Parallel 4D-Variational Data Assimilation for an Eulerian Chemistry Transport Model

Author(s):  
H. Elbern ◽  
H. Schmidt ◽  
A. Ebel
2013 ◽  
Vol 13 (14) ◽  
pp. 7225-7240 ◽  
Author(s):  
J. Barré ◽  
L. El Amraoui ◽  
P. Ricaud ◽  
W. A. Lahoz ◽  
J.-L. Attié ◽  
...  

Abstract. The behavior of the extratropical transition layer (ExTL) is investigated using a chemistry transport model (CTM) and analyses derived from assimilation of MLS (Microwave Limb Sounder) O3 and MOPITT (Measurements Of Pollution In The Troposphere) CO data. We firstly focus on a stratosphere–troposphere exchange (STE) case study that occurred on 15 August 2007 over the British Isles (50° N, 10° W). We evaluate the effect of data assimilation on the O3–CO correlations. It is shown that data assimilation disrupts the relationship in the transition region. When MLS O3 is assimilated, CO and O3 values are not consistent between each other, leading to unphysical correlations at the STE location. When MLS O3 and MOPITT CO assimilated fields are taken into account in the diagnostics the relationship happens to be more physical. We then use O3–CO correlations to quantify the effect of data assimilation on the height and depth of the ExTL. When the free-model run O3 and CO fields are used in the diagnostics, the ExTL distribution is found 1.1 km above the thermal tropopause and is 2.6 km wide (2σ). MOPITT CO analyses only slightly sharpen (by −0.02 km) and lower (by −0.2 km) the ExTL distribution. MLS O3 analyses provide an expansion (by +0.9 km) of the ExTL distribution, suggesting a more intense O3 mixing. However, the MLS O3 analyses ExTL distribution shows a maximum close to the thermal tropopause and a mean location closer to the thermal tropopause (+0.45 km). When MLS O3 and MOPITT CO analyses are used together, the ExTL shows a mean location that is the closest to the thermal tropopause (+0.16 km). We also extend the study at the global scale on 15 August 2007 and for the month of August 2007. MOPITT CO analyses still show a narrower chemical transition between stratosphere and troposphere than the free-model run. MLS O3 analyses move the ExTL toward the troposphere and broaden it. When MLS O3 analyses and MOPITT CO analyses are used together, the ExTL matches the thermal tropopause poleward of 50°.


2009 ◽  
Vol 9 (2) ◽  
pp. 6691-6737 ◽  
Author(s):  
S. Massart ◽  
C. Clerbaux ◽  
D. Cariolle ◽  
A. Piacentini ◽  
S. Turquety ◽  
...  

Abstract. The Infrared Atmospheric Sounding Interferometer (IASI) is one of the five European new generation instruments carried by the polar-orbiting MetOp-A satellite. Data assimilation is a powerful tool to combine these data with a numerical model. This paper presents the first steps made towards the assimilation of the total ozone columns from the IASI measurements into a chemistry transport model. The IASI ozone data used are provided by an inversion of radiances performed at the LATMOS (Laboratoire Atmosphères, Milieux, Observations Spatiales). As a contribution to the validation of this dataset, the LATMOS-IASI data are compared to a four dimensional ozone field, with low systematic and random errors compared to ozonesondes and OMI-DOAS data. This field results from the combined assimilation of ozone profiles from the MLS instrument and of total ozone columns from the SCIAMACHY instrument. It is found that on average, the LATMOS-IASI data tends to overestimate the total ozone columns by 2% to 8%. The random observation error of the LATMOS-IASI data is estimated to about 6%, except over polar regions and deserts where it is higher. Using this information, the LATMOS-IASI data are then assimilated, combined with the MLS data. This first LATMOS-IASI data assimilation experiment shows that the resulting analysis is quite similar to the one obtained from the combined MLS and SCIAMACHY data assimilation.


2016 ◽  
Author(s):  
Sergey Skachko ◽  
Richard Menard ◽  
Quentin Errera ◽  
Yves Christophe ◽  
Simon Chabrillat

Abstract. We compare two optimized chemical data assimilation systems, one based on the ensemble Kalman filter (EnKF) and the other based on four-dimensional variational (4D-Var), using a comprehensive stratospheric chemistry transport model (CTM). The work is an extension of the Belgian Assimilation System for Chemical ObsErvations (BASCOE), initially designed to work with a 4D-Var data assimilation. A strict comparison of both methods in the case of chemical tracer transport was done in a previous study and indicated that both methods provide essentially similar results. In the present work, we assimilate observations of ozone, HCl, HNO3, H2O and N2O from EOS Aura-MLS data into the BASCOE CTM with a full description of stratospheric chemistry. Two new issues related to the use of full chemistry model with EnKF are taken into account. One issue concerns to a large number of error variance parameters that need to be optimized. We estimate an observation error parameter as function of pressure level for each observed species using the Desroziers' method. For comparison reasons, we apply the same estimate procedure in the 4D-Var data assimilation, where we keep both estimates: the background and observation error variances. However in EnKF, the background error covariance is modelled using the full chemistry model and a model error term. We found that it is adequate to have a single model error based on the chemical tracer formulation that is applied for all species. This is an indication that the main source of model error in chemical transport model is due to the transport. The second issue in EnKF with comprehensive atmospheric chemistry models is the sampling errors between species. When species are weakly chemically related, cross-species sampling noise errors occur at the same location. These errors need to be filtered out, in addition to a localization based on distance. The performance of two data assimilation methods was assessed through an eight-month long assimilation of limb sounding observations from EOS Aura-MLS. The paper discusses the differences in results and their relation to stratospheric chemical processes. Generally speaking, EnKF and 4D-Var provide results of comparable quality but differ substantially in presence of model error or observation biases. If the erroneous chemical modelling is associated with not too small chemical life-times, then EnKF performs better, while 4D-Var develops spurious increments in the chemically related species. If, on the other hand, the observation biases are significant, then 4D-Var is more robust and is able to reject erroneous observations, while EnKF does not.


2013 ◽  
Vol 6 (1) ◽  
pp. 1-16 ◽  
Author(s):  
J. D. Silver ◽  
J. Brandt ◽  
M. Hvidberg ◽  
J. Frydendall ◽  
J. H. Christensen

Abstract. Data assimilation is the process of combining real-world observations with a modelled geophysical field. The increasing abundance of satellite retrievals of atmospheric trace gases makes chemical data assimilation an increasingly viable method for deriving more accurate analysed fields and initial conditions for air quality forecasts. We implemented a three-dimensional optimal interpolation (OI) scheme to assimilate retrievals of NO2 tropospheric columns from the Ozone Monitoring Instrument into the Danish Eulerian Hemispheric Model (DEHM, version V2009.0), a three-dimensional, regional-scale, offline chemistry-transport model. The background error covariance matrix, B, was estimated based on differences in the NO2 concentration field between paired simulations using different meteorological inputs. Background error correlations were modelled as non-separable, horizontally homogeneous and isotropic. Parameters were estimated for each month and for each hour to allow for seasonal and diurnal patterns in NO2 concentrations. Three experiments were run to compare the effects of observation thinning and the choice of observation errors. Model performance was assessed by comparing the analysed fields to an independent set of observations: ground-based measurements from European air-quality monitoring stations. The analysed NO2 and O3 concentrations were more accurate than those from a reference simulation without assimilation, with increased temporal correlation for both species. Thinning of satellite data and the use of constant observation errors yielded a better balance between the observed increments and the prescribed error covariances, with no appreciable degradation in the surface concentrations due to the observation thinning. Forecasts were also considered and these showed rather limited influence from the initial conditions once the effects of the diurnal cycle are accounted for. The simple OI scheme was effective and computationally feasible in this context, where only a single species was assimilated, adjusting the three-dimensional field for this compound. Limitations of the assimilation scheme are discussed.


Sign in / Sign up

Export Citation Format

Share Document