A linear CO chemistry parameterization in a chemistry-transport model: evaluation and application to data assimilation

This paper presents an evaluation of a new linear parameterization valid for the troposphere and the stratosphere, based on a first order approximation of the carbon monoxide (CO) continuity equation. This linear scheme (hereinafter noted LINCO) has been implemented in the 3-D Chemical Transport Model (CTM) MOCAGE (MOdèle de Chimie Atmospherique Grande Echelle). First, a one and a half years of LINCO simulation has been compared to output obtained from a detailed chemical scheme output. The mean differences between both schemes are about ±25 ppbv (part per billion by volume) or 15% in the troposphere and ±10 ppbv or 100% in the stratosphere. Second, LINCO has been compared to diverse observations from satellite instruments covering the troposphere (Measurements Of Pollution In The Troposphere: MOPITT) and the stratosphere (Microwave Limb Sounder: MLS) and also from aircraft (Measurements of ozone and water vapour by Airbus in-service aircraft: MOZAIC programme) mostly flying in the upper troposphere and lower stratosphere (UTLS). In the troposphere, the LINCO seasonal variations as well as the vertical and horizontal distributions are quite close to MOPITT CO observations. However, a bias of ~−40 ppbv is observed at 700 Pa between LINCO and MOPITT. In the stratosphere, MLS and LINCO present similar large-scale patterns, except over the poles where the CO concentration is underestimated by the model. In the UTLS, LINCO presents small biases less than 2% compared to independent MOZAIC profiles. Third, we assimilated MOPITT CO using a variational 3D-FGAT (First Guess at Appropriate Time) method in conjunction with MOCAGE for a long run of one and a half years. The data assimilation greatly improves the vertical CO distribution in the troposphere from 700 to 350 hPa compared to independent MOZAIC profiles. At 146 hPa, the assimilated CO distribution is also improved compared to MLS observations by reducing the bias up to a factor of 2 in the tropics. This study confirms that the linear scheme is able to simulate reasonably well the CO distribution in the troposphere and in the lower stratosphere. Therefore, the low computing cost of the linear scheme opens new perspectives to make free runs and CO data assimilation runs at high resolution and over periods of several years

How realistic are air quality hindcasts driven by forcings from climate model simulations?

Predicting how European air quality could evolve over the next decades in the context of changing climate requires the use of climate models to produce results that can be averaged in a climatologically and statistically sound manner. This is a very different approach from the one that is generally used for air quality hindcasts for the present period; analysed meteorological fields are used to represent specifically each date and hour. Differences arise both from the fact that a climate model run results in a pure model output, with no influence from observations (which are useful to correct for a range of errors), and that in a "climate" set-up, simulations on a given day, month or even season cannot be related to any specific period of time (but can just be interpreted in a climatological sense). Hence, although an air quality model can be thoroughly validated in a "realistic" set-up using analysed meteorological fields, the question remains of how far its outputs can be interpreted in a "climate" set-up. For this purpose, we focus on Europe and on the current decade using three 5-yr simulations performed with the multiscale chemistry-transport model MOCAGE and use meteorological forcings either from operational meteorological analyses or from climate simulations. We investigate how statistical skill indicators compare in the different simulations, discriminating also the effects of meteorology on atmospheric fields (winds, temperature, humidity, pressure, etc.) and on the dependent emissions and deposition processes (volatile organic compound emissions, deposition velocities, etc.). Our results show in particular how differing boundary layer heights and deposition velocities affect horizontal and vertical distributions of species. When the model is driven by operational analyses, the simulation accurately reproduces the observed values of O<sub>3</sub>, NO<sub>x</sub>, SO<sub>2</sub> and, with some bias that can be explained by the set-up, PM<sub>10</sub>. We study how the simulations driven by climate forcings differ, both due to the realism of the forcings (lack of data assimilated and lower resolution) and due to the lack of representation of the actual chronology of events. We conclude that the indicators such as mean bias, mean normalized bias, RMSE and deviation standards can be used to interpret the results with some confidence as well as the health-related indicators such as the number of days of exceedance of regulatory thresholds. These metrics are thus considered to be suitable for the interpretation of simulations of the future evolution of European air quality

Equilibrium of sinks and sources of sulphate over Europe: comparison between a six-year simulation and EMEP observations

Sulphate distributions were simulated with a global chemistry transport model. A chemical scheme describing the sulphur cycle and the parameterisations of the main sinks for sulphate aerosols were included in the model. A six-year simulation was conducted from the years 2000 to 2005, driven by the ECMWF operational analyses. Emissions come from an inventory representative of the year 2000. This paper focuses on the analysis of the sulphate sinks and sources over Europe for the entire period of simulation. The Sulphate burden shows a marked annual cycle, which is the result of the annual variations of the aqueous and gaseous chemistry. Regionally, the monthly mean aerosol burden can vary by a factor of 2 from one year to another, because of different weather conditions, driving chemistry, transport and wet deposition of sulphate aerosols. Sulphate ground concentrations, scavenging fluxes and precipitation modelled were compared with observations. The model represents quite well sulphate fields over Europe, but has a general tendency to overestimate sulphate ground concentrations, in particular over Northern Europe. We assume that it is linked to the representation of the scavenging fluxes, which are underestimated. We suggest that uncertainties in modelled precipitation explain only partially the underestimation of the scavenging fluxes in the model

A new version of the CNRM Chemistry-Climate Model, CNRM-CCM: description and improvements from the CCMVal-2 simulations

This paper presents a new version of the Météo-France CNRM Chemistry-Climate Model, so-called CNRM-CCM. It includes some fundamental changes from the previous version (CNRM-ACM) which was extensively evaluated in the context of the CCMVal-2 validation activity. The most notable changes concern the radiative code of the GCM, and the inclusion of the detailed stratospheric chemistry of our Chemistry-Transport model MOCAGE on-line within the GCM. A 47-yr transient simulation (1960–2006) is the basis of our analysis. CNRM-CCM generates satisfactory dynamical and chemical fields in the stratosphere. Several shortcomings of CNRM-ACM simulations for CCMVal-2 that resulted from an erroneous representation of the impact of volcanic aerosols as well as from transport deficiencies have been eliminated. <br><br> Remaining problems concern the upper stratosphere (5 to 1 hPa) where temperatures are too high, and where there are biases in the NO<sub>2</sub>, N<sub>2</sub>O<sub>5</sub> and O<sub>3</sub> mixing ratios. In contrast, temperatures at the tropical tropopause are too cold. These issues are addressed through the implementation of a more accurate radiation scheme at short wavelengths. Despite these problems we show that this new CNRM CCM is a useful tool to study chemistry-climate applications

Midlatitude stratosphere - troposphere exchange as diagnosed by MLS O3 and MOPITT CO assimilated fields

International audienceThis paper presents a comprehensive characterization of a very deep stratospheric intrusion which occurred over the British Isles on 15 August 2007. The signature of this event is diagnosed using ozonesonde measurements over Lerwick, UK (60.14 N, 1.19 W) and is also well characterized using meteorological analyses from the global operational weather prediction model of Météo-France, ARPEGE. Modelled as well as assimilated fields of both ozone (O3) and carbon monoxide (CO) have been used in order to better document this event. O3 and CO from Aura/MLS and Terra/MOPITT instruments, respectively, are assimilated into the three-dimensional chemical transport model MOCAGE of Météo-France using a variational 3-DFGAT (First Guess at Appropriate Time) method. The validation of O3 and CO assimilated fields is done using selfconsistency diagnostics and by comparison with independent observations such as MOZAIC (O3 and CO), AIRS (CO) and OMI (O3). It particularly shows in the upper troposphere and lower stratosphere region that the assimilated fields are closer to MOZAIC than the free model run. The O3 bias between MOZAIC and the analyses is −11.5 ppbv with a RMS of 22.4 ppbv and a correlation coefficient of 0.93, whereas between MOZAIC and the free model run, the corresponding values are 33 ppbv, 38.5 ppbv and 0.83, respectively. In the same way, for CO, the bias, RMS and correlation coefficient between MOZAIC and the analyses are −3.16 ppbv, 13 ppbv and 0.79, respectively, whereas between MOZAIC and the free model run, the corresponding values are 33 ppbv, 38.5 ppbv and 0.83, respectively. In the same way, for CO, the bias, RMS and correlation coefficient between MOZAIC and the analyses are −3.16 ppbv, 13 ppbv and 0.79, respectively, whereas between MOZAIC and the free model they are 6.3 ppbv, 16.6 ppbv and 0.71, respectively. The paper also presents a demonstration of the capability of O3 and CO assimilated fields to better describe a stratosphere-troposphere exchange (STE) event in comparison with the free run modelled O3 and CO fields. Although the assimilation of MLS data improves the distribution of O3 above the tropopause compared to the free model run, it is not sufficient to reproduce the STE event well. Assimilated MOPITT CO allows a better qualitative description of the stratospheric intrusion event. The MOPITT CO analyses appear more promising than the MLS O3 analyses in terms of their ability to capture a deep STE event. Therefore, the results of this study open the perspectives for using MOPITT CO in the STE studies

A new chemistry-climate tropospheric and stratospheric model MOCAGE-Climat: evaluation of the present-day climatology and sensitivity to surface processes

International audienceWe present the chemistry-climate configuration of the Météo-France Chemistry and Transport Model, MOCAGE-Climat. MOCAGE-Climat is a state-of-the-art model that simulates the global distribution of ozone and its precursors (82 chemical species) both in the troposphere and the stratosphere, up to the mid-mesosphere (~70 km). Surface processes (emissions, dry deposition), convection, and scavenging are explicitly described in the model that has been driven by the ECMWF operational analyses of the period 2000–2005, on T21 and T42 horizontal grids and 60 hybrid vertical levels, with and without a procedure that reduces calculations in the boundary layer, and with on-line or climatological deposition velocities. Model outputs have been compared to available observations, both from satellites (TOMS, HALOE, SMR, SCIAMACHY, MOPITT) and in-situ instrument measurements (ozone sondes, MOZAIC and aircraft campaigns) at climatological timescales. The distribution of long-lived species is in fair agreement with observations in the stratosphere putting apart shortcomings linked to the large-scale circulation. The variability of the ozone column, both spatially and temporarily, is satisfactory. However, the too fast Brewer-Dobson circulation accumulates too much ozone in the lower to mid-stratosphere at the end of winter. Ozone in the UTLS region does not show any systematic bias. In the troposphere better agreement with ozone sonde measurements is obtained at mid and high latitudes than in the tropics and differences with observations are the lowest in summer. Simulations using a simplified boundary layer lead to ozone differences between the model and the observations up to the mid-troposphere. NOx in the lowest troposphere is in general overestimated, especially in the winter months over the northern hemisphere, which might result from a positive bias in OH. Dry deposition fluxes of O3 and nitrogen species are within the range of values reported by recent inter-comparison model exercises. The use of climatological deposition velocities versus deposition velocities calculated on-line had greatest impact on HNO3 and NO2 in the troposphere

A geostationary thermal infrared sensor to monitor the lowermost troposphere: O₃ and CO retrieval studies

This paper describes the capabilities of a nadir thermal infrared (TIR) sensor proposed for deployment onboard a geostationary platform to monitor ozone (O3) and carbon monoxide (CO) for air quality (AQ) purposes. To assess the capabilities of this sensor we perform idealized retrieval studies considering typical atmospheric profiles of O3 and CO over Europe with different instrument configuration (signal to noise ratio, SNR, and spectral sampling interval, SSI) using the KOPRA forward model and the KOPRA-fit retrieval scheme. We then select a configuration, referred to as GEO-TIR, optimized for providing information in the lowermost troposphere (LmT; 0–3 km in height). For the GEO-TIR configuration we obtain ~1.5 degrees of freedom for O3 and ~2 for CO at altitudes between 0 and 15 km. The error budget of GEO-TIR, calculated using the principal contributions to the error (namely, temperature, measurement error, smoothing error) shows that information in the LmT can be achieved by GEO-TIR. We also retrieve analogous profiles from another geostationary infrared instrument with SNR and SSI similar to the Meteosat Third Generation Infrared Sounder (MTG-IRS) which is dedicated to numerical weather prediction, referred to as GEO-TIR2. We quantify the added value of GEO-TIR over GEO-TIR2 for a realistic atmosphere, simulated using the chemistry transport model MOCAGE (MOd`ele de Chimie Atmospherique `a Grande Echelle). Results show that GEO-TIR is able to capture well the spatial and temporal variability in the LmT for both O3 and CO. These results also provide evidence of the significant added value in the LmT of GEO-TIR compared to GEO-TIR2 by showing GEO-TIR is closer to MOCAGE than GEO-TIR2 for various statistical parameters (correlation, bias, standard deviation)

Evaluation of the MOCAGE Chemistry Transport Model during the ICARTT/ITOP Experiment

We evaluate the Meteo-France global chemistry transport 3D model MOCAGE (MOdele de Chimie Atmospherique a Grande Echelle) using the important set of aircraft measurements collected during the ICARRT/ITOP experiment. This experiment took place between US and Europe during summer 2004 (July 15-August 15). Four aircraft were involved in this experiment providing a wealth of chemical data in a large area including the North East of US and western Europe. The model outputs are compared to the following species of which concentration is measured by the aircraft: OH, H2O2, CO, NO, NO2, PAN, HNO3, isoprene, ethane, HCHO and O3. Moreover, to complete this evaluation at larger scale, we used also satellite data such as SCIAMACHY NO2 and MOPITT CO. Interestingly, the comprehensive dataset allowed us to evaluate separately the model representation of emissions, transport and chemical processes. Using a daily emission source of biomass burning, we obtain a very good agreement for CO while the evaluation of NO2 points out incertainties resulting from inaccurate ratio of emission factors of NOx/CO. Moreover, the chemical behavior of O3 is satisfactory as discussed in the paper

Ensemble forecasts of air quality in eastern China – Part 1: Model description and implementation of the MarcoPolo–Panda prediction system, version 1

An operational multi-model forecasting system for air quality including nine different chemical transport models has been developed and provides daily forecasts of ozone, nitrogen oxides, and particulate matter for the 37 largest urban areas of China (population higher than 3 million in 2010). These individual forecasts as well as the mean and median concentrations for the next 3 days are displayed on a publicly accessible website (http://www.marcopolo-panda.eu, last access: 7 December 2018). The paper describes the forecasting system and shows some selected illustrative examples of air quality predictions. It presents an intercomparison of the different forecasts performed during a given period of time (1–15 March 2017) and highlights recurrent differences between the model output as well as systematic biases that appear in the median concentration values. Pathways to improve the forecasts by the multi-model system are suggested.</p