Three-dimensional ozone data analysis with an air quality model over the Paris area

[1] We present a description, an evaluation, and a comparison of four methods designed to produce objective and physically consistent maps of ozone concentration fields. These methods are based on the use of a chemistry-transport model (CTM) and available observations. In most existing analysis systems, the error covariance is modeled using assumptions of homogeneity and isotropy. However, these assumptions may fail in the case of strongly heterogeneous terrain or emission patterns. Therefore we propose a simple method for specifying an anisotropic and heterogeneous background error covariance model in a statistical interpolation. We illustrate the positive impact of the implementation of this model. Since the covariance model is independent of the state of the atmosphere and invariant in time, we simultaneously test kriging techniques that are generally used for spatial interpolations. Kriging is applied to the observed values or to the differences between simulated and observed concentrations, namely the innovations. We perform an objective statistical validation of all methods for data recorded over an entire summer season, in contrast to other evaluations of data assimilation methods made in air quality modeling. We demonstrate that the RMS error of the ozone analyses at the surface is 30–50% smaller than the one from simulations, regardless of the method used. Most of the time, the kriging method applied to the innovations gives results equivalent to that of the anisotropic statistical interpolation, in spite of very different formulations. In addition, we show that this method outperforms the classical kriging method applied to only observations. The information provided by the CTM is therefore essential to a good-quality representation of ozone patterns such as city plumes or urban/suburban gradients. We also use the Atmospheric Pollution Over the Paris Area (ESQUIF) airborne measurements to demonstrate that the anisotropic method efficiently corrects CTM fields in altitude.

[1]  Arjo Segers,et al.  Data assimilation of ozone in the atmospheric transport chemistry model LOTOS , 2000, Environ. Model. Softw..

[2]  Gérald Desroziers,et al.  A Coordinate Change for Data Assimilation in Spherical Geometry of Frontal Structures , 1997 .

[3]  Michael E. Chang,et al.  On using inverse methods for resolving emissions with large spatial inhomogeneities , 1997 .

[4]  A. Russell,et al.  Estimation of emission adjustments from the application of four-dimensional data assimilation to photochemical air quality modeling , 2001 .

[5]  Henk Eskes,et al.  Assimilation of total ozone satellite measurements in a three‐dimensional tracer transport model , 1999 .

[6]  P. Woodward,et al.  Application of the Piecewise Parabolic Method (PPM) to meteorological modeling , 1990 .

[7]  Robert Vautard,et al.  A comparison of simulated and observed ozone mixing ratios for the summer of 1998 in Western Europe , 2001 .

[8]  R. Vautard,et al.  Measurements and modelling of atmospheric pollution over the Paris area: an overview of the ESQUIF Project , 2000 .

[9]  Valerii V. Fedorov,et al.  Kriging and other Estimators of Spatial Field Characteristics (With Special Reference to Environmental Studies) , 1989 .

[10]  Montserrat Fuentes,et al.  A high frequency kriging approach for non‐stationary environmental processes , 2001 .

[11]  Robert Vautard,et al.  Validation of a hybrid forecasting system for the ozone concentrations over the Paris area , 2001 .

[12]  O. Talagrand,et al.  4D-variational data assimilation with an adjoint air quality model for emission analysis , 2000, Environ. Model. Softw..

[13]  P. L. Finkelstein The Spatial Analysis of Acid Precipitation Data. , 1984 .

[14]  J. M. Lewis,et al.  The use of adjoint equations to solve a variational adjustment problem with advective constraints , 1985 .

[15]  R. E. Kalman,et al.  New Results in Linear Filtering and Prediction Theory , 1961 .

[16]  Eric Mayer,et al.  Analytical determination and classification of pollutant concentration fields using air pollution monitoring network data: Methodology and application in the Paris area, during episodes with peak nitrogen dioxide levels , 2000, Environ. Model. Softw..

[17]  Lars Peter Riishojgaard,et al.  A direct way of specifying flow‐dependent background error correlations for meteorological analysis systems , 1998 .

[18]  M. Ghil,et al.  Data assimilation in meteorology and oceanography , 1991 .

[19]  R. E. Kalman,et al.  A New Approach to Linear Filtering and Prediction Problems , 2002 .

[20]  S. Cohn,et al.  Ooce Note Series on Global Modeling and Data Assimilation Construction of Correlation Functions in Two and Three Dimensions and Convolution Covariance Functions , 2022 .

[21]  Jan G. Verwer,et al.  Gauss-Seidel Iteration for Stiff ODES from Chemical Kinetics , 1994, SIAM J. Sci. Comput..

[22]  J. Hoelzemann,et al.  PSAS and 4D-var data assimilation for chemical state analysis by urban and rural observation sites , 2001 .

[23]  B. Hicks,et al.  Some factors that affect the deposition rates of sulfur dioxide and similar gases on vegetation , 1977 .

[24]  R. Kalman,et al.  New results in linear prediction and filtering theory Trans. AMSE , 1961 .

[25]  Matthias Beekmann,et al.  Effect of biogenic volatile organic compound emissions on tropospheric chemistry during the Atmospheric Pollution Over the Paris Area (ESQUIF) campaign in the Ile-de-France region , 2003 .

[26]  J. Erisman,et al.  Parametrization of surface resistance for the quantification of atmospheric deposition of acidifying pollutants and ozone , 1994 .

[27]  Pierre H. Flamant,et al.  A synthesis of the Air Pollution Over the Paris Region (ESQUIF) field campaign , 2003 .

[28]  J. O'Brien,et al.  A Note on the Vertical Structure of the Eddy Exchange Coefficient in the Planetary Boundary Layer , 1970 .

[29]  Noel A Cressie,et al.  Statistics for Spatial Data, Revised Edition. , 1994 .

[30]  R. Daley Atmospheric Data Analysis , 1991 .

[31]  A. Raftery,et al.  Model Validation and Spatial Interpolation by Combining Observations with Outputs from Numerical Models via Bayesian Melding , 2001 .

[32]  P. Rasch,et al.  MOZART, a global chemical transport model for ozone , 1998 .

[33]  M. Ghil,et al.  A stochastic-dynamic model for the spatial structure of forecast error statistics , 1983 .

[34]  Michael E. Chang,et al.  Inverse modeling of biogenic isoprene emissions , 1996 .

[35]  Anthony Hollingsworth,et al.  The statistical structure of short-range forecast errors as determined from radiosonde data , 1986 .

[36]  Hendrik Elbern,et al.  Ozone episode analysis by four-dimensional variational chemistry data assimilation , 2001 .

[37]  Rainer Friedrich,et al.  Paris emission inventory diagnostics from ESQUIF airborne measurements and a chemistry transport model , 2003 .

[38]  Mireille Lattuati,et al.  Impact des emissions europeennes sur le bilan de l'ozone tropospherique a l'interface de l'europe et de l'atlantique nord : apport de la modelisation lagrangienne et des mesures en altitude , 1997 .

[39]  A. Gilliland,et al.  A sensitivity study of the discrete Kalman filter (DKF) to initial condition discrepancies , 2001 .