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

Polynomial growth of volume of balls for zero-entropy geodesic systems

The aim of this paper is to state and prove polynomial analogues of the classical Manning inequality relating the topological entropy of a geodesic flow with the growth rate of the volume of balls in the universal covering. To this aim we use two numerical conjugacy invariants, the {\em strong polynomial entropy $h_{pol}$} and the {\em weak polynomial entropy $h_{pol}^*$}. Both are infinite when the topological entropy is positive and they satisfy $h_{pol}^*\leq h_{pol}$. We first prove that the growth rate of the volume of balls is bounded above by means of the strong polynomial entropy and we show that for the flat torus this inequality becomes an equality. We then study the explicit example of the torus of revolution for which we can give an exact asymptotic equivalent of the growth rate of volume of balls, which we relate to the weak polynomial entropy.Comment: 22 page

Adaptive density estimation for stationary processes

We propose an algorithm to estimate the common density $s$ of a stationary process $X_1,...,X_n$. We suppose that the process is either $\beta$ or $\tau$-mixing. We provide a model selection procedure based on a generalization of Mallows' $C_p$ and we prove oracle inequalities for the selected estimator under a few prior assumptions on the collection of models and on the mixing coefficients. We prove that our estimator is adaptive over a class of Besov spaces, namely, we prove that it achieves the same rates of convergence as in the i.i.d framework

Combined assimilation of IASI and MLS observations to constrain tropospheric and stratospheric ozone in a global chemical transport model

Accurate and temporally resolved fields of free-troposphere ozone are of major importance to quantify the intercontinental transport of pollution and the ozone radiative forcing. We consider a global chemical transport model (MOdèle de Chimie Atmosphérique à Grande Échelle, MOCAGE) in combination with a linear ozone chemistry scheme to examine the impact of assimilating observations from the Microwave Limb Sounder (MLS) and the Infrared Atmospheric Sounding Interferometer (IASI). The assimilation of the two instruments is performed by means of a variational algorithm (4D-VAR) and allows to constrain stratospheric and tropospheric ozone simultaneously. The analysis is first computed for the months of August and November 2008 and validated against ozonesonde measurements to verify the presence of observations and model biases. Furthermore, a longer analysis of 6 months (July–December 2008) showed that the combined assimilation of MLS and IASI is able to globally reduce the uncertainty (root mean square error, RMSE) of the modeled ozone columns from 30 to 15% in the upper troposphere/lower stratosphere (UTLS, 70–225 hPa). The assimilation of IASI tropospheric ozone observations (1000–225 hPa columns, TOC – tropospheric O₃ column) decreases the RMSE of the model from 40 to 20% in the tropics (30° S–30° N), whereas it is not effective at higher latitudes. Results are confirmed by a comparison with additional ozone data sets like the Measurements of OZone and wAter vapour by aIrbus in-service airCraft (MOZAIC) data, the Ozone Monitoring Instrument (OMI) total ozone columns and several high-altitude surface measurements. Finally, the analysis is found to be insensitive to the assimilation parameters. We conclude that the combination of a simplified ozone chemistry scheme with frequent satellite observations is a valuable tool for the long-term analysis of stratospheric and free-tropospheric ozone

Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and Rural Africa

peer reviewe

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

Forecasting global atmospheric CO_2

A new global atmospheric carbon dioxide (CO_2) real-time forecast is now available as part of the pre-operational Monitoring of Atmospheric Composition and Climate – Interim Implementation (MACC-II) service using the infrastructure of the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS). One of the strengths of the CO_2 forecasting system is that the land surface, including vegetation CO_2 fluxes, is modelled online within the IFS. Other CO_2 fluxes are prescribed from inventories and from off-line statistical and physical models. The CO_2 forecast also benefits from the transport modelling from a state-of-the-art numerical weather prediction (NWP) system initialized daily with a wealth of meteorological observations. This paper describes the capability of the forecast in modelling the variability of CO_2 on different temporal and spatial scales compared to observations. The modulation of the amplitude of the CO_2 diurnal cycle by near-surface winds and boundary layer height is generally well represented in the forecast. The CO_2 forecast also has high skill in simulating day-to-day synoptic variability. In the atmospheric boundary layer, this skill is significantly enhanced by modelling the day-to-day variability of the CO_2 fluxes from vegetation compared to using equivalent monthly mean fluxes with a diurnal cycle. However, biases in the modelled CO_2 fluxes also lead to accumulating errors in the CO_2 forecast. These biases vary with season with an underestimation of the amplitude of the seasonal cycle both for the CO_2 fluxes compared to total optimized fluxes and the atmospheric CO_2 compared to observations. The largest biases in the atmospheric CO_2 forecast are found in spring, corresponding to the onset of the growing season in the Northern Hemisphere. In the future, the forecast will be re-initialized regularly with atmospheric CO_2 analyses based on the assimilation of CO_2 products retrieved from satellite measurements and CO_2 in situ observations, as they become available in near-real time. In this way, the accumulation of errors in the atmospheric CO_2 forecast will be reduced. Improvements in the CO_2 forecast are also expected with the continuous developments in the operational IFS

Discriminant Analysis of Geographical Origin of Cork Planks and Stoppers by Near Infrared Spectroscopy

The objective of this study was to assess the potential of visible and near infrared spectroscopy (VIS+NIRS) combined with multivariate analysis for identifying the geographical origin of cork. The study was carried out on cork planks and natural cork stoppers from the most representative cork-producing areas in the world. Two training sets of international and national cork planks were studied. The first set comprised a total of 479 samples from Morocco, Portugal, and Spain, while the second set comprised a total of 179 samples from the Spanish regions of Andalusia, Catalonia, and Extremadura. A training set of 90 cork stoppers from Andalusia and Catalonia was also studied. Original spectroscopic data were obtained for the transverse sections of the cork planks and for the body and top of the cork stoppers by means of a 6500 Foss-NIRSystems SY II spectrophotometer using a fiber optic probe. Remote reflectance was employed in the wavelength range of 400 to 2500 nm. After analyzing the spectroscopic data, discriminant models were obtained by means of partial least square (PLS) with 70% of the samples. The best models were then validated using 30% of the remaining samples. At least 98% of the international cork plank samples and 95% of the national samples were correctly classified in the calibration and validation stage. The best model for the cork stoppers was obtained for the top of the stoppers, with at least 90% of the samples being correctly classified. The results demonstrate the potential of VIS + NIRS technology as a rapid and accurate method for predicting the geographical origin of cork plank and stopper

Sequential schemes for frequentist estimation of properties in statistical model checking

National Research Foundation (NRF) Singapor

Comparison of multivariate calibration techniques applied to experimental NIR data sets

The present study compares the performance of different multivariate calibration techniques applied to four near-infrared data sets when test samples are well within the calibration domain. Three types of problems are discussed: the nonlinear calibration, the calibration using heterogeneous data sets, and the calibration in the presence of irrelevant information in the set of predictors. Recommendations are derived from the comparison, which should help to guide a nonchemometrician through the selection of an appropriate calibration method for a particular type of calibration data. A flexible methodology is proposed to allow selection of an appropriate calibration technique for a given calibration problem.54460862