Automatic generation of large ensembles for air quality forecasting using the Polyphemus system

Abstract. This paper describes a method to automatically generate a large ensemble of air quality simulations. Such an ensemble may be useful for quantifying uncertainty, improving forecasts, evaluating risks, identifying process weaknesses, etc. The objective is to take into account all sources of uncertainty: input data, physical formulation and numerical formulation. The leading idea is to build different chemistry-transport models in the same framework, so that the ensemble generation can be fully controlled. Large ensembles can be generated with a Monte Carlo simulations that address at the same time the uncertainties in the input data and in the model formulation. This is achieved using the Polyphemus system, which is flexible enough to build various different models. The system offers a wide range of options in the construction of a model: many physical parameterizations, several numerical schemes and different input data can be combined. In addition, input data can be perturbed. In this paper, some 30 alternatives are available for the generation of a model. For each alternative, the options are given a probability, based on how reliable they are supposed to be. Each model of the ensemble is defined by randomly selecting one option per alternative. In order to decrease the computational load, as many computations as possible are shared by the models of the ensemble. As an example, an ensemble of 101 photochemical models is generated and run for the year 2001 over Europe. The models' performance is quickly reviewed, and the ensemble structure is analyzed. We found a strong diversity in the results of the models and a wide spread of the ensemble. It is noteworthy that many models turn out to be the best model in some regions and some dates.

[1]  M. Wesely.  Parameterization of surface resistances to gaseous dry deposition in regional-scale numerical models , 1989 .

[2]  Rohit Mathur,et al.  Assessment of an ensemble of seven real-time ozone forecasts over eastern North America during the summer of 2004 , 2005 .

[3]  Bruno Sportisse,et al.  MICS Asia Phase II - Sensitivity to the aerosol module , 2007, 0704.3190.

[4]  Robert Vet,et al.  A revised parameterization for gaseous dry deposition in air-quality models , 2003 .

[5]  S. Madronich.  Photodissociation in the atmosphere: 1. Actinic flux and the effects of ground reflections and clouds , 1987 .

[6]  G. Evensen.  Sequential data assimilation with a nonlinear quasi‐geostrophic model using Monte Carlo methods to forecast error statistics , 1994 .

[7]  Vivien Mallet,et al.  Investigation of some numerical issues in a chemistry-transport model: Gas-phase simulations , 2007 .

[8]  Willem Hundsdorfer,et al.  Numerical time integration for air pollution models , 1998 .

[9]  Edouard Debry,et al.  Atmospheric Chemistry and Physics Technical Note : A new SIze REsolved Aerosol Model ( SIREAM ) , 2007 .

[10]  J. Louis.  A parametric model of vertical eddy fluxes in the atmosphere , 1979 .

[11]  M. Sofiev,et al.  Ensemble dispersion forecasting—Part I: concept, approach and indicators , 2004 .

[12]  Steven E. Peckham,et al.  Application of dynamic linear regression to improve the skill of ensemble-based deterministic ozone forecasts , 2006 .

[13]  F. Kirchner,et al.  A new mechanism for regional atmospheric chemistry modeling , 1997 .

[14]  H. Christopher Frey,et al.  Uncertainties in predicted ozone concentrations due to input uncertainties for the UAM-V photochemical grid model applied to the July 1995 OTAG domain , 2001 .

[15]  Vivien Mallet,et al.  Development and validation of a fully modular platform for numerical modelling of air pollution: POLAIR , 2004 .

[16]  Matthias Beekmann,et al.  Monte Carlo uncertainty analysis of a regional‐scale transport chemistry model constrained by measurements from the Atmospheric Pollution Over the Paris Area (ESQUIF) campaign , 2003 .

[17]  W. Stockwell,et al.  Comparison of the EMEP, RADM2 and RACM Mechanisms , 2003 .

[18]  M. Smith,et al.  The sea spray generation function , 1998 .

[19]  D. E. Spiel,et al.  A Model of Marine Aerosol Generation Via Whitecaps and Wave Disruption , 1986 .

[20]  Philippe Thunis,et al.  Evaluation of long-term ozone simulations from seven regional air quality models and their ensemble , 2007 .

[21]  Willem Hundsdorfer,et al.  Numerical time integration for air pollution problems , 2002 .

[22]  I. Troen,et al.  A simple model of the atmospheric boundary layer; sensitivity to surface evaporation , 1986 .

[23]  C. N. Hewitt,et al.  Inventorying emissions from nature in Europe , 1999 .

[24]  Anne Grete Straume,et al.  A More Extensive Investigation of the Use of Ensemble Forecasts for Dispersion Model Evaluation , 2001 .

[25]  S. Hanna,et al.  Monte carlo estimates of uncertainties in predictions by a photochemical grid model (UAM-IV) due to uncertainties in input variables , 1998 .

[26]  V. Mallet,et al.  Uncertainty in a chemistry-transport model due to physical parameterizations and numerical approximations: An ensemble approach applied to ozone modeling , 2006 .

[27]  Vivien Mallet,et al.  Ensemble‐based air quality forecasts: A multimodel approach applied to ozone , 2006 .

[28]  W. Stockwell,et al.  The second generation regional acid deposition model chemical mechanism for regional air quality modeling , 1990 .

[29]  Arjo Segers,et al.  Data Assimilation in Atmospheric Chemistry Models Using Kalman Filtering , 2002 .

[30]  D. Byun.  Science algorithms of the EPA Models-3 community multi-scale air quality (CMAQ) modeling system , 1999 .

[31]  Vivien Mallet,et al.  Ozone ensemble forecast with machine learning algorithms , 2009 .

[32]  Lin Wu,et al.  Technical Note: The air quality modeling system Polyphemus , 2007 .

[33]  Paulette Middleton,et al.  A three‐dimensional Eulerian acid deposition model: Physical concepts and formulation , 1987 .