Estimating European volatile organic compound emissions using satellite observations of formaldehyde from the Ozone Monitoring Instrument

Emission of non-methane Volatile Organic Compounds (VOCs) to the atmosphere stems from biogenic and human activities, and their estimation is difficult because of the many and not fully understood processes involved. In order to narrow down the uncertainty related to VOC emissions, which negatively reflects on our ability to simulate the atmospheric composition, we exploit satellite observations of formaldehyde (HCHO), an ubiquitous oxidation product of most VOCs, focusing on Europe. HCHO column observations from the Ozone Monitoring Instrument (OMI) reveal a marked seasonal cycle with a summer maximum and winter minimum. In summer, the oxidation of methane and other long-lived VOCs supply a slowly varying background HCHO column, while HCHO variability is dominated by most reactive VOC, primarily biogenic isoprene followed in importance by biogenic terpenes and anthropogenic VOCs. The chemistry-transport model CHIMERE qualitatively reproduces the temporal and spatial features of the observed HCHO column, but display regional biases which are attributed mainly to incorrect biogenic VOC emissions, calculated with the Model of Emissions of Gases and Aerosol from Nature (MEGAN) algorithm. These "bottom-up" or a-priori emissions are corrected through a Bayesian inversion of the OMI HCHO observations. Resulting "top-down" or a-posteriori isoprene emissions are lower than "bottom-up" by 40% over the Balkans and by 20% over Southern Germany, and higher by 20% over Iberian Peninsula, Greece and Italy. We conclude that OMI satellite observations of HCHO can provide a quantitative "top-down" constraint on the European "bottom-up" VOC inventories

STRATAQ: A three-dimensional Chemical Transport Model of the stratosphere

International audienceA three-dimensional (3-D) Chemical Transport Model (CTM) of the stratosphere has been developed and used for a test study of the evolution of chemical species in the arctic lower stratosphere during winter 1996/97. This particular winter has been chosen for testing the model's capabilities for its remarkable dynamical situation (very cold and strong polar vortex) along with the availability of sparse chlorine, HNO3 and O3 data, showing also very low O3 values in late March/April. Due to those unusual features, the winter 1996/97 can be considered an excellent example of the impact of both dynamics and heterogeneous reactions on the chemistry of the stratosphere. Model integration has been performed from January to March 1997 and the resulting long-lived and short-lived tracer fields compared with available measurements. The model includes a detailed gas phase chemical scheme and a parameterization of the heterogeneous reactions occurring on liquid aerosol and polar stratospheric cloud (PSC) surfaces. The transport is calculated using a semi-lagrangian flux scheme, forced by meteorological analyses. In such form, the STRATAQ CTM model is suitable for short-term integrations to study transport and chemical evolution related to "real" meteorological situations. Model simulation during the chosen winter shows intense PSC formation, with noticeable local HNO3 capture by PSCs, and the activation of vortex air leading to chlorine production and subsequent O3 destruction. The resulting model fields show generally good agreement with satellite data (MLS and TOMS), although the available observations, due to their limited number and time/space sparse nature, are not enough to effectively constraint the model. In particular, the model seems to perform well in reproducing the rapid processing of air inside the polar vortex on PSC converting reservoir species in active chlorine. In addition, it satisfactorily reproduces the morphology of the continuous O3 decline as shown by the satellite during the investigated period, with a tendency, however, to underestimate the total column values inside the polar vortex during late winter. As possible causes of this model/observation difference we suggest an incorrect estimation of the vertical transport and of the tropospheric contribution

Assimilation of stratospheric ozone in the chemical transport model STRATAQ

International audienceWe describe a sequential assimilation approach useful for assimilating tracer measurements into a three-dimensional chemical transport model (CTM) of the stratosphere. The numerical code, developed largely according to Kha00, uses parameterizations and simplifications allowing assimilation of sparse observations and the simultaneous evaluation of analysis errors, with reasonable computational requirements. Assimilation parameters are set by using ?2 and OmF (Observation minus Forecast) statistics. The CTM used here is a high resolution three-dimensional model. It includes a detailed chemical package and is driven by UKMO (United Kingdom Meteorological Office) analyses. We illustrate the method using assimilation of Upper Atmosphere Research Satellite/Microwave Limb Sounder (UARS/MLS) ozone observations for three weeks during the 1996 antarctic spring. The comparison of results from the simulations with TOMS (Total Ozone Mapping Spectrometer) measurements shows improved total ozone fields due to assimilation of MLS observations. Moreover, the assimilation gives indications on a possible model weakness in reproducing polar ozone values during springtime

Quinpi: Integrating Conservation Laws with CWENO Implicit Methods

Many interesting applications of hyperbolic systems of equations are stiff, and require the time step to satisfy restrictive stability conditions. One way to avoid small time steps is to use implicit time integration. Implicit integration is quite straightforward for first-order schemes. High order schemes instead also need to control spurious oscillations, which requires limiting in space and time also in the linear case. We propose a framework to simplify considerably the application of high order non-oscillatory schemes through the introduction of a low order implicit predictor, which is used both to set up the nonlinear weights of a standard high order space reconstruction, and to achieve limiting in time. In this preliminary work, we concentrate on the case of a third-order scheme, based on diagonally implicit Runge Kutta (DIRK) integration in time and central weighted essentially non-oscillatory (CWENO) reconstruction in space. The numerical tests involve linear and nonlinear scalar conservation laws

Fundamental diagrams in traffic flow: the case of heterogeneous kinetic models

Experimental studies on vehicular traffic provide data on quantities like density, flux, and mean speed of the vehicles. However, the diagrams relating these variables (the fundamental and speed diagrams) show some peculiarities not yet fully reproduced nor explained by mathematical models. In this paper, resting on the methods of kinetic theory, we introduce a new traffic model which takes into account the heterogeneous nature of the flow of vehicles along a road. In more detail, the model considers traffic as a mixture of two or more populations of vehicles (e.g., cars and trucks) with different microscopic characteristics, in particular different lengths and/or maximum speeds. With this approach we gain some insights into the scattering of the data in the regime of congested traffic clearly shown by actual measurements

Lidar observations of mesospheric sodium over Italy

We have developed a lidar system for the measurements of the mesospheric atomic sodium (Na) number density profile. The performances and the main characteristics of the system are evaluated showing the first observations: Na number density profiles can be obtained with a resolution of 300 m in altitude and of 200 s in time, the relative uncertainty is between 5% and 10%. These features also allow to discuss the capability of the system to detect the time and space variability of Na density profiles, which is typical of the gravity wave propagation. The preliminary analysis of the sodium profile dynamics along a 7 hours measurement session, started from 19:30UT on 10 April 1997, shows the presence of wave-like perturbations with characteristics typical of propagating gravity waves

Preliminary results of a Mesoscale Model for MARS

The Earth atmospheric circulation has been studied for long time using both GCM (General Circulation Models) and Mesoscale Models or LAM (Limited Area Models). The latter have been widely applied to study local circulation at high resolution and for weather forecasting. In the last years, the Martian atmosphere arouse the interest of the scientific community, both for supporting the landing of Beagle 2 lander and for studying and assessing similarities/differences with the Earth atmosphere. To this aim, GCM have been successfully used. Recently, also Earth LAMs have been changed to simulate the Mars atmosphere, showing good results. The scarce availability of observations did not allowed for validating these models. In this work an attempt is made to validate the newly developed MARS-MM5 against GCM. The model simulation produced using a data base on the basis of output from multi-annual integration of two CGM (see Lewis S. R. et al., J. Geophys. Res., 104 (E10) (1999) 177) is used for statistically evaluates MARS-MM5. The preliminary results suggest that MARS-MM5 is able to correctly reproduce the Mars atmosphere, indeed either the horizontal and the vertical structure of temperature produced by MARS-MM5 is in good agreement with the ones produced by GCM. A few discrepancies are found in the PBL, probably produced by a different parameterization

Validation and analysis of regional present-day climate and climate change simulations over Europe

In the European Commission (EC) project "Regionalization of Anthropogenic Climate Change Simulations, RACCS, recently terminated, 11 European institutions have carried out tests of dynamical and statistical regionalization techniques. The outcome of the "dynamical part" of the project, utilizing a series of high resolution LAMs and a variable resolution global model (all of which we shall refer to as RCMs, Regional Climate Models), is presented here. The per- formance of the dqterent LAMs had first, in a preceding EC project, been tested with "perfect" boundary forcing fields (ECMWF analyses) and also multi-year present-day climate simula- tions with AMIP "perfect ocean " or mixed layer ocean GCM boundary conditions had been validated against available climatological data. The present report involves results of vali- dation and analysis of RCM present-day climate simulations and anthropogenic climate change experiments. Multi-year (5 - 30 years) present-day climate simulations have been per- formed with resolutions between 19 and 70 km (grid lengths) and with boundary conditions from the newest CGCM simulations. The climate change experiments involve various 2xCO2 - ]xCO2 transient greenhouse gas experiments and in one case also changing sulphur aerosols. A common validation and inter-comparison was made at the coordinating institution, MPIfor Meteorology. The validation of the present-day climate simulations shows the importance of systematic errors in the low level general circulation. Such errors seem to induce large errors in precipitation and surface air temperature in the RCMs as well as in the CGCMs providing boundary conditions. Over Europe the field of systematic errors in the mean sea level pressure (MSLP) usually involve an area of too low pressure, often in the form of an east-west trough across Europe with too high pressure to the north and south. New storm-track analyses confirm that the areas of too low pressure are caused by enhanced cyclonic activity and similarly that the areas of too high pressure are caused by reduced such activity. The precise location and strength of the extremes in the MSLP error field seems to be dependent on the physical param- eterization package used. In model pairs sharing the same package the area of too low pressure is deepened further in the RCM compared to the corresponding CGCM, indicating an increase of the excessive cyclonic activity with increasing resolution. From the experiments performed it seems not possible to decide to what extent the systematic errors in the general circulation are the result of local errors in the physical parameterization schemes or remote errors trans- mitted to the European region via the boundary conditions. Additional errors in precipitation and temperature seems to be due to direct local effects of errors in certain parameterization schemes and errors in the SSTs taken from the CGCMs. For all seasons many biases are fOund to be statistically significant compared to estimates of the internal model variability of the time- slice mean values. In the climate change experiments statistically significant European mean temperature changes which are large compared to the corresponding biases are found. How- ever, the changes in the deviations from the European mean temperature as well as the changes in precipitation are only partly sign wcan ce and are of the same order of magnitude or smaller than the corresponding biases found in the present-day climate simulations. Cases of an inter- action between the systematic model errors and the radiative forcing show that generally the errors are not canceling out when the changes are computed. Therefore, reliable regional cli- mate changes can only be achieved after model improvements which reduce their systematic errors sufficiently. Also in future RCM experiments sujiciently long time-slices must be used in order to obtain statistically sign ijicant climate changes on the sub-continental scale aimed at with the present regionalization technique