scholarly journals EnKF and 4D-Var Data Assimilation with a Chemistry Transport Model

2016 ◽  
Author(s):  
Sergey Skachko ◽  
Richard Menard ◽  
Quentin Errera ◽  
Yves Christophe ◽  
Simon Chabrillat
Keyword(s):  
Data Assimilation ◽  
Transport Model ◽  
Model Error ◽  
The Other ◽  
Full Description ◽  
Observation Error ◽  
Tracer Transport ◽  
Chemistry Transport Model ◽  
Stratospheric Chemistry ◽  
Chemical Tracer

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.

scholarly journals EnKF and 4D-Var data assimilation with chemical transport model BASCOE (version 05.06)

2016 ◽  
Vol 9 (8) ◽  
pp. 2893-2908 ◽  
Author(s):  
Sergey Skachko ◽  
Richard Ménard ◽  
Quentin Errera ◽  
Yves Christophe ◽  
Simon Chabrillat
Keyword(s):  
Data Assimilation ◽  
Transport Model ◽  
Chemical Transport ◽  
Error Variance ◽  
Model Error ◽  
Chemical Processes ◽  
Observation Error ◽  
Chemical Transport Model ◽  
Error Covariance ◽  
Chemical Tracer

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) data assimilation, using a comprehensive stratospheric chemistry transport model (CTM). This 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 the full chemistry model with EnKF are taken into account. One issue is a large number of error variance parameters that need to be optimized. We estimate an observation error variance parameter as a function of pressure level for each observed species using the Desroziers method. For comparison purposes, we apply the same estimate procedure in the 4D-Var data assimilation, where both scale factors of the background and observation error covariance matrices are estimated using the Desroziers method. However, in EnKF the background error covariance is modelled using the full chemistry model and a model error term which is tuned using an adjustable parameter. We found that it is adequate to have the same value of this parameter based on the chemical tracer formulation that is applied for all observed 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 noise in the cross-covariance between species that occurs when species are weakly chemically related 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 8-month long assimilation of limb sounding observations from EOS Aura MLS. This 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 the presence of model error or observation biases. If the erroneous chemical modelling is associated with moderately fast chemical processes, but whose lifetimes are longer than the model time step, then EnKF performs better, while 4D-Var develops spurious increments in the chemically related species. If, however, the observation biases are significant, then 4D-Var is more robust and is able to reject erroneous observations while EnKF does not.

scholarly journals First steps towards the assimilation of IASI ozone data into the MOCAGE-PALM system

2009 ◽  
Vol 9 (2) ◽  
pp. 6691-6737 ◽  
Author(s):  
S. Massart ◽  
C. Clerbaux ◽  
D. Cariolle ◽  
A. Piacentini ◽  
S. Turquety ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Total Ozone ◽  
Transport Model ◽  
Polar Regions ◽  
Observation Error ◽  
Random Errors ◽  
Chemistry Transport Model ◽  
Ozone Data ◽  
The One ◽  
New Generation

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.

scholarly journals Diagnosing the transition layer at extratropical latitudes using MLS O<sub>3</sub> and MOPITT CO analyses

2013 ◽  
Vol 13 (14) ◽  
pp. 7225-7240 ◽  
Author(s):  
J. Barré ◽  
L. El Amraoui ◽  
P. Ricaud ◽  
W. A. Lahoz ◽  
J.-L. Attié ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Transition Region ◽  
Transition Layer ◽  
Transport Model ◽  
Global Scale ◽  
British Isles ◽  
Chemistry Transport Model ◽  
Free Model ◽  
The Relationship

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°.

scholarly journals Assimilation of surface NO<sub>2</sub> and O<sub>3</sub> observations into the SILAM chemistry transport model

2015 ◽  
Vol 8 (2) ◽  
pp. 191-203 ◽  
Author(s):  
J. Vira ◽  
M. Sofiev
Keyword(s):  
Transport Model ◽  
Time Of Day ◽  
Covariance Matrices ◽  
Atmospheric Composition ◽  
Integrated Modelling ◽  
Daily Maximum ◽  
Observation Error ◽  
Data Set ◽  
Chemistry Transport Model ◽  
Error Covariance

Abstract. This paper describes the assimilation of trace gas observations into the chemistry transport model SILAM (System for Integrated modeLling of Atmospheric coMposition) using the 3D-Var method. Assimilation results for the year 2012 are presented for the prominent photochemical pollutants ozone (O3) and nitrogen dioxide (NO2). Both species are covered by the AirBase observation database, which provides the observational data set used in this study. Attention was paid to the background and observation error covariance matrices, which were obtained primarily by the iterative application of a posteriori diagnostics. The diagnostics were computed separately for 2 months representing summer and winter conditions, and further disaggregated by time of day. This enabled the derivation of background and observation error covariance definitions, which included both seasonal and diurnal variation. The consistency of the obtained covariance matrices was verified using χ2 diagnostics. The analysis scores were computed for a control set of observation stations withheld from assimilation. Compared to a free-running model simulation, the correlation coefficient for daily maximum values was improved from 0.8 to 0.9 for O3 and from 0.53 to 0.63 for NO2.

scholarly journals Assimilation of surface NO<sub>2</sub> and O<sub>3</sub> observations into the SILAM chemistry transport model

2014 ◽  
Vol 7 (4) ◽  
pp. 5589-5621 ◽  
Author(s):  
J. Vira ◽  
M. Sofiev
Keyword(s):  
Transport Model ◽  
Model Simulation ◽  
Time Of Day ◽  
Covariance Matrices ◽  
Daily Maximum ◽  
Observation Error ◽  
Chemistry Transport Model ◽  
Error Covariance ◽  
Seasonal And Diurnal Variation ◽  
Free Running

Abstract. This paper describes assimilation of trace gas observations into the chemistry transport model SILAM using the 3D-Var method. Assimilation results for year 2012 are presented for the prominent photochemical pollutants ozone (O3) and nitrogen dioxide (NO2). Both species are covered by the Airbase observation database, which provides the observational dataset used in this study. Attention is paid to the background and observation error covariance matrices, which are obtained primarily by iterative application of a posteriori diagnostics. The diagnostics are computed separately for two months representing summer and winter conditions, and further disaggregated by time of day. This allows deriving background and observation error covariance definitions which include both seasonal and diurnal variation. The consistency of the obtained covariance matrices is verified using χ2 diagnostics. The analysis scores are computed for a control set of observation stations withheld from assimilation. Compared to a free-running model simulation, the correlation coefficient for daily maximum values is improved from 0.8 to 0.9 for O3 and from 0.53 to 0.63 for NO2.

scholarly journals Representing model error in ensemble data assimilation

2014 ◽  
Vol 21 (5) ◽  
pp. 971-985 ◽  
Author(s):  
C. Cardinali ◽  
N. Žagar ◽  
G. Radnoti ◽  
R. Buizza
Keyword(s):  
Data Assimilation ◽  
Covariance Matrix ◽  
Model Error ◽  
The Other ◽  
Error Statistics ◽  
Model Uncertainties ◽  
Background Error ◽  
Error Covariance ◽  
Ensemble Data ◽  
Assimilation System

Abstract. The paper investigates a method to represent model error in the ensemble data assimilation (EDA) system. The ECMWF operational EDA simulates the effect of both observations and model uncertainties. Observation errors are represented by perturbations with statistics characterized by the observation error covariance matrix whilst the model uncertainties are represented by stochastic perturbations added to the physical tendencies to simulate the effect of random errors in the physical parameterizations (ST-method). In this work an alternative method (XB-method) is proposed to simulate model uncertainties by adding perturbations to the model background field. In this way the error represented is not just restricted to model error in the usual sense but potentially extends to any form of background error. The perturbations have the same correlation as the background error covariance matrix and their magnitude is computed from comparing the high-resolution operational innovation variances with the ensemble variances when the ensemble is obtained by perturbing only the observations (OBS-method). The XB-method has been designed to represent the short range model error relevant for the data assimilation window. Spread diagnostic shows that the XB-method generates a larger spread than the ST-method that is operationally used at ECMWF, in particular in the extratropics. Three-dimensional normal-mode diagnostics indicate that XB-EDA spread projects more than the spread from the other EDAs onto the easterly inertia-gravity modes associated with equatorial Kelvin waves, tropical dynamics and, in general, model error sources. The background error statistics from the above described EDAs have been employed in the assimilation system. The assimilation system performance showed that the XB-method background error statistics increase the observation influence in the analysis process. The other EDA background error statistics, when inflated by a global factor, generate analyses with 30–50% smaller degree of freedom of signal. XB-EDA background error variances have not been inflated. The presented EDAs have been used to generate the initial perturbations of the ECMWF ensemble prediction system (EPS) of which the XB-EDA induces the largest EPS spread, also in the medium range, leading to a more reliable ensemble. Compared to ST-EDA, XB-EDA leads to a small improvement of the EPS ignorance skill score at day 3 and 7.

scholarly journals The global chemistry transport model TM5: description and evaluation of the tropospheric chemistry version 3.0

2010 ◽  
Vol 3 (3) ◽  
pp. 1009-1087 ◽  
Author(s):  
V. Huijnen ◽  
J. E. Williams ◽  
M. van Weele ◽  
T. P. C. van Noije ◽  
M. C. Krol ◽  
...  
Keyword(s):  
Seasonal Variation ◽  
Nitrogen Dioxide ◽  
Transport Model ◽  
Surface Ozone ◽  
Oxidative Capacity ◽  
Tropospheric Chemistry ◽  
Full Description ◽  
Ozone Monitoring Instrument ◽  
Chemistry Transport Model ◽  
Spatial And Seasonal Variation

Abstract. We present a comprehensive description and benchmark evaluation of the tropospheric chemistry version of the global chemistry transport model TM5 (Tracer Model 5, version TM5-chem-v3.0). A full description is given concerning the photochemical mechanism, the interaction with aerosol, the treatment of the stratosphere, the wet and dry deposition parameterizations, and the applied emissions. We evaluate the model against a suite of ground-based, satellite, and aircraft measurements of components critical for understanding global photochemistry for the year 2006. The model exhibits a realistic oxidative capacity at a global scale. The methane lifetime is ~8.9 years with an associated lifetime of methyl chloroform of 5.86 years, which is similar to that derived using an optimized hydroxyl radical field. The seasonal cycle in observed carbon monoxide (CO) is well simulated at different regions across the globe. In the Northern Hemisphere CO concentrations are underestimated by about 20 ppbv in spring and 10 ppbv in summer, which is related to missing chemistry and underestimated emissions from higher hydrocarbons, as well as to uncertainties in the seasonal variation of CO emissions. The model also captures the spatial and seasonal variation in formaldehyde tropospheric columns as observed by SCIAMACHY. Positive model biases over the Amazon and eastern United States point to uncertainties in the isoprene emissions as well as its chemical breakdown. Simulated tropospheric nitrogen dioxide columns correspond well to observations from the Ozone Monitoring Instrument in terms of its seasonal and spatial variability (with a global spatial correlation coefficient of 0.89), but TM5 fields are lower by 25–40%. This is consistent with earlier studies pointing to a high bias of 0–30% in the OMI retrievals, but uncertainties in the emission inventories have probably also contributed to the discrepancy. TM5 tropospheric nitrogen dioxide profiles are in good agreement (within ~0.1 ppbv) with in situ aircraft observations from the INTEX-B campaign over (the Gulf of) Mexico. The model reproduces the spatial and seasonal variation in background surface ozone concentrations and tropospheric ozone profiles from the World Ozone and Ultraviolet Radiation Data Centre to within 10 ppbv, but at several tropical stations the model tends to underestimate ozone in the free troposphere. The presented model results benchmark the TM5 tropospheric chemistry version, which is currently in use in several international cooperation activities, and upon which future model improvements will take place.

Sign in / Sign up

Export Citation Format

Share Document