170 likes | 360 Views
4D-var assimilation of CRISTA-NF H20 and MLS retrievals with the high resolution SACADA system. Ketevan Kasradze kk@riu.uni-koeln.de Supervisor: PD Dr. Hendrik Elbern Rhenish Institute for Environmental Research at the University of Cologne Germany. Atmospheric layers. 3 /18. Why?.
E N D
4D-var assimilation of CRISTA-NF H20 and MLS retrievals with the high resolution SACADA system KetevanKasradze kk@riu.uni-koeln.de Supervisor: PD Dr. HendrikElbern Rhenish Institute for Environmental Research at the University of Cologne Germany
Atmospheric layers 3/18
Why? 4/18
Data assimilation „using all the available information, to determine as accurately as possible the state of the atmospheric (or oceanic) flow” Talagrand (1997) available information: • observations • background 5/18
Data assimilation Background Analysis Xb Xa Xb 1st day 2nd day t 6/18
Quasi-Newton method L-BFGS SACADA assimilation-system Model operator Background Projection operator Cost function Observation error covariance matrix Gradient Adjungiertes Modell Vector of observations Background error covariance matrix 54x23 042 ~ 1x10^6 BECM ~ 4x10^11 for 51 species ~ 80 Terrabyte 7/18
SACADA assimilation-system Background Meteorological ECMWF analyses Trace gas observations PREP SACADA DWD GME CTM CTMad Diffusion L-BFGS Analysis 8/18
Horizontal GME Grid • ~147km between the grid points, • 23 042 grid points pro Model layer 9/18
Additional refinement troposphere/lower stratosphere SACADA Vertical Grid 54 layer MLS CRISTA-NF 10/18
Chemistry • 51 species • 200 chemical reaction equations • Reaction rates: JPL 2006 recommendations • KPP (Kinetic Preprocessor) • ROS2 • Chemistry Module • Adjoint Chemistry Module 11/18
SCOUT-O3campaignStratospheric-Climate Links with Emphasis on the UTLS - O3November-December 2005 AMMA-campaignAfricanMonsoonMultidisciplinaryAnalyses 29.07.2006 -17.08.2006 12/18
Assimilation Results • Background • Derived from SOCRATES (=2D) model results (1-st day) • At the following day: Analyses of the day before • Observations • MLS retrievals • CRISTA-NF retrievals • Analyses • Ozone analyses from SACADA for 7.11.2005 • Pressure level 30 (~ 137hPa) • H2O profiles from SACADA for 29.07.2006 13/18
Ozone O3 7.11.2005 ~137hPa ~14km 12 hUTC Analysis Observation Background 14/18
20km 20km 10km 10km CRISTA-NF - Background MLS Asm: H2O 29.07.2006 15/18
Computationalcosts Fine grid (ni=48, lev=54) with new chemistry mechanism, timestep: 800sec. 1 day assimilation: • ~50 GB temporary storage • ~10 clock hour (16 AMD OPTERON processor) 16/18
Acknowledgements to: Garmany. University of Cologne. Rhenish Institute for Environmental Research at the University of Cologne. Prof A. Ebel, Prof A. Wahner, Dr. J. Schwinger, all my co-workers and especially my supervisor Dr. H. Elbern. Tbilisi State University. Applied mathematics and computer sciences department stuff. Prof. R. Botchorishvili, Dr. M. Menteshashvili.
Thank you! Danke schön! დიდი მადლობა!