ketevan kasradze [email protected] supervisor: pd dr. hendrik elbern

Tbilisi, 6.08.2012

www.riu.uni-koeln.de Ketevan Kasradze

Ketevan [email protected]

Supervisor: PD Dr. Hendrik Elbern

Rhenish Institute for Environmental Researchat the University of Cologne

Germany

4D-var assimilation of CRISTA-NF H20 and MLS retrievals with the high resolution SACADA system

mailto:[email protected]

Tbilisi, 6.08.2012


Atmospheric layers

3/18

Tbilisi, 6.08.2012


Why?

4/18

Tbilisi, 6.08.2012


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

Tbilisi, 6.08.2012


Data assimilation

1st day 2nd day t

Xa

Xb

BackgroundAnalysis

Xb

6/18

Tbilisi, 6.08.2012


N

iii

Tibx yHxRHMxxBxJ

0

1*0

10 )(

0

N

iii

Tiib

Tb yxMHRyxMHxxBxxxJ

00

100

100 2

121)(

Background error covariance matrix54x23 042 ~ 1x10^6BECM ~ 4x10^11 for 51 species ~ 80 Terrabyte

Cost function

Vector of observations

Observation error covariance matrix

Projection operator

BackgroundModel operator

GradientAdjungiertes Modell

)(),(, 000 0xJxJx x 0xQuasi-Newton method

L-BFGS

SACADA assimilation-system

7/18

Tbilisi, 6.08.2012


SACADA assimilation-system

Background Meteorological ECMWF analyses

Trace gas observations

Analysis

SACADA

PREPDWD GMECTM

CTMadDiffusion

L-BFGS

8/18

Tbilisi, 6.08.2012


Horizontal GME Grid

• ~147km between the grid points,• 23 042 grid points pro Model layer

9/18

Tbilisi, 6.08.2012

www.riu.uni-koeln.de Ketevan KasradzeAd

ditio

nal r

efine

men

t tro

posp

here

/low

er

stra

tosp

here

SACADA Vertical Grid54 layer

CRISTA-NF

MLS

10/18

Tbilisi, 6.08.2012


Chemistry

H2O2 H2O NO3

N2O N2O5 HNO3 HNO4 Cl2O2

OClOCl2

HClHOCl ClONO2

BrClHBr

HOBrBrONO2

O(1D) O(3P)

O3 H

OHHO2

NONO2 Cl ClO

BrBrOH2

CH4

CH3ClCCl4

CFCl3

CF2Cl2

CHF2ClCF2ClCFCl2

CH3CCl3

CH3Br CF2ClBr

CF3BrCH2O CH3O

CH3O2 CH3OOH

CH3 HCO CO O2

N2

51 species 200 chemical reaction equations

Reaction rates: JPL 2006 recommendations

KPP (Kinetic Preprocessor)• ROS2• Chemistry Module• Adjoint Chemistry Module

11/18

Tbilisi, 6.08.2012


SCOUT-O3 campaign Stratospheric-Climate Links with Emphasis on the UTLS - O3

November-December 2005

AMMA-campaign African Monsoon Multidisciplinary Analyses

29.07.2006 -17.08.2006

12/18

Tbilisi, 6.08.2012


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

Tbilisi, 6.08.2012


A

naly

sis

Obs

erva

tion

Ba

ckgr

ound

Ozone O3

7.11.2005

~137hPa ~14km

12 h UTC

14/18

Tbilisi, 6.08.2012


CRISTA-NF - Background MLS

29.07.2006

Asm: H2O20km

10km

20km

10km

15/18

Tbilisi, 6.08.2012


Computational costs

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

Tbilisi, 6.08.2012


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.

Tbilisi, 6.08.2012


Thank you!Danke schön!

!დიდი მადლობა

ketevan kasradze [email protected] supervisor: pd dr. hendrik elbern

Documents