Lightning and convection parameterisations ? uncertainties in global modelling

International audienceThe simulation of convection, lightning and consequent NOx emissions with global atmospheric chemistry models is associated with large uncertainties since these processes are heavily parameterised. Each parameterisation by itself has deficiencies and the combination of these substantially increases the uncertainties compared to the individual parameterisations. In this study several combinations of state-of-the-art convection and lightning parameterisations are used in simulations with the global atmospheric chemistry general circulation model ECHAM5/MESSy, and are evaluated against lightning observations. A wide range in the spatial and temporal variability of the simulated flash densities is found, attributed to both types of parameterisations. Some combinations perform well, whereas others are hardly applicable. In addition to resolution dependent rescaling parameters, each combination of lightning and convection schemes requires individual scaling to reproduce the observed flash frequencies. The resulting NOx profiles are inter-compared, however definite conclusions about the most realistic profiles can currently not be drawn

Quantifying atmospheric transport, chemistry, and mixing using a new trajectory-box model and a global atmospheric-chemistry GCM

We present a novel method for the quantification of transport, chemistry, and mixing along atmospheric trajectories based on a consistent model hierarchy. The hierarchy consists of the new atmospheric-chemistry trajectory-box model CAABA/MJT and the three-dimensional (3-D) global ECHAM/MESSy atmospheric-chemistry (EMAC) general circulation model. CAABA/MJT employs the atmospheric box model CAABA in a configuration using the atmospheric-chemistry submodel MECCA (M), the photochemistry submodel JVAL (J), and the new trajectory submodel TRAJECT (T), to simulate chemistry along atmospheric trajectories, which are provided offline. With the same chemistry submodels coupled to the 3-D EMAC model and consistent initial conditions and physical parameters, a unique consistency between the two models is achieved. Since only mixing processes within the 3-D model are excluded from the model consistency, comparisons of results from the two models allow to separate and quantify contributions of transport, chemistry, and mixing along the trajectory pathways. Consistency of transport between the trajectory-box model CAABA/MJT and the 3-D EMAC model is achieved via calculation of kinematic trajectories based on 3-D wind fields from EMAC using the trajectory model LAGRANTO. The combination of the trajectory-box model CAABA/MJT and the trajectory model LAGRANTO can be considered as a Lagrangian chemistry-transport model (CTM) moving isolated air parcels. The procedure for obtaining the necessary statistical basis for the quantification method is described as well as the comprehensive diagnostics with respect to chemistry. <br><br> The quantification method presented here allows to investigate the characteristics of transport, chemistry, and mixing in a grid-based 3-D model. The analysis of chemical processes within the trajectory-box model CAABA/MJT is easily extendable to include, for example, the impact of different transport pathways or of mixing processes onto chemistry. Under certain prerequisites described here, the results can be used to complement observations with detailed information about the history of observed air masses

Technical note: Implementation of prescribed (OFFLEM), calculated (ONLEM), and pseudo-emissions (TNUDGE) of chemical species in the Modular Earth Submodel System (MESSy)

International audienceWe present the submodels OFFLEM, ONLEM, and TNUDGE for the Modular Earth Submodel System (MESSy). Prescribed emissions from input files are handled by OFFLEM. ONLEM deals with online-calculated emissions, i.e., emissions that are calculated during the simulation. The submodel TNUDGE uses the "tracer nudging" technique for pseudo-sources and -sinks. For species with highly uncertain emission fluxes and/or with sufficiently long lifetimes, e.g., CH4, it is common to create such pseudo-fluxes by prescribing the observed mixing ratio of the species at a given boundary (e.g., the mixing ratio of methane at the surface, or the ozone mixing ratio at the tropopause). All three submodels substantially simplify the inclusion of emissions into a model. Specific emissions can easily be switched on or off. New prescribed emissions can be included without rewriting any code. New online emissions only require one additional subroutine containing the new parameterization. A major advantage is that input fields at arbitrary resolution can be used. The problem of incompatible grids between emission data and model is overcome by utilizing the MESSy data import interface. To further simplify the creation of new offline emission data, the preprocessing program EDGAR2NC is provided. EDGAR2NC transforms files from the EDGAR format into the netCDF format which is required by OFFLEM. The presented routines are a part of the community modeling project MESSy and can be made available for use to the atmospheric modeling community

Technical Note: Simulation of detailed aerosol chemistry on the global scale using MECCA-AERO

International audienceWe present the MESSy submodel MECCA-AERO, which simulates both aerosol and gas phase chemistry within one comprehensive mechanism. Including the aerosol phase into the chemistry mechanism increases the stiffness of the resulting set of differential equations. The numerical aspects of the approach followed in MECCA-AERO are presented. MECCA-AERO requires input of an aerosol dynamical/microphysical model to provide the aerosol size and particle number information of the modes/bins for which the chemistry is explicitly calculated. Additional precautions are required to avoid the double counting of processes, especially for sulphate in the aerosol dynamical and the chemistry model. This coupling is explained in detail. To illustrate the capabilities of the new aerosol submodel, examples for species usually treated in aerosol dynamical models are shown. The aerosol chemistry as provided by MECCA-AERO is very sumptuous and not readily applicable for long-term simulations, though it provides a reference to evaluate simplified approaches

Technical note: A new comprehensive SCAVenging submodel for global atmospheric chemistry modelling

International audienceWe present the new scavenging scheme SCAV, simulating the removal of trace gases and aerosol particles by clouds and precipitation in global atmospheric chemistry models. The scheme is quite flexible and can be used for various purposes, e.g. long term chemistry simulations as well as detailed cloud and precipitation chemistry calculations. The presence of clouds can substantially change the chemical composition of the atmosphere. We present a new method of mechanistically coupling gas phase, aerosol, cloud and precipitation chemistry, which enables studies of feedbacks between multiphase chemistry and transport processes

Technical Note: An implementation of the dry removal processes DRY DEPosition and SEDImentation in the Modular Earth Submodel System (MESSy)

International audienceWe present the submodels DRYDEP and SEDI for the Modular Earth Submodel System (MESSy). Dry deposition of gases and aerosols is calculated within DRYDEP, whereas SEDI deals with aerosol particle sedimentation. Dry deposition velocities depend on the near-surface turbulence and the physical and chemical properties of the surface cover (e.g. the roughness length, soil pH or leaf stomatal exchange). The dry deposition algorithm used in DRYDEP is based on the big leaf approach and is described in detail within this Technical Note. The sedimentation submodel SEDI contains two sedimentation schemes: a simple upwind zeroth order scheme and a first order approach

Global cloud and precipitation chemistry and wet deposition: tropospheric model simulations with ECHAM5/MESSy1

International audienceThe representation of cloud and precipitation chemistry and subsequent wet deposition of trace constituents in global atmospheric chemistry models is associated with large uncertainties. To improve the simulated trace gas distributions we apply the new submodel SCAV, which includes detailed cloud and precipitation chemistry and present results of the atmospheric chemistry general circulation model ECHAM5/MESSy1. A good agreement with observed wet deposition fluxes for species causing acid rain is obtained. The new scheme enables prognostic calculations of the pH of clouds and precipitation, and these results are also in accordance with observations. We address the influence of detailed cloud and precipitation chemistry on trace constituents based on sensitivity simulations. The results confirm previous results from regional scale and box models, and we extend the analysis to the role of aqueous phase chemistry on the global scale. Some species are directly affected through multiphase removal processes, and many also indirectly through changes in oxidant concentrations, which in turn have an impact on the species lifetime. While the overall effect on tropospheric ozone is relatively small (3 can reach ?20%, and several important compounds (e.g., H2O2, HCHO) are substantially depleted by clouds and precipitation

Technical Note: Coupling of chemical processes with the Modular Earth Submodel System (MESSy) submodel TRACER

International audienceThe implementation of processes related to chemistry into Earth System Models and their coupling within such systems requires the consistent description of the chemical species involved. We provide a tool (written in Fortran95) to structure and manage information about constituents, hereinafter referred to as tracers, namely the Modular Earth Submodel System (MESSy) generic (i.e., infrastructure) submodel TRACER. With TRACER it is possible to define a multitude of tracer sets, depending on the spatio-temporal representation (i.e., the grid structure) of the model. The required information about a specific chemical species is split into the static meta-information about the characteristics of the species, and its (generally in time and space variable) abundance in the corresponding representation. TRACER moreover includes two submodels. One is TRACER_FAMILY, an implementation of the tracer family concept. It distinguishes between two types: type-1 families are usually applied to handle strongly related tracers (e.g., fast equilibrating species) for a specific process (e.g., advection). In contrast to this, type-2 families are applied for tagging techniques. Tagging means the artificial decomposition of one or more species into parts, which are additionally labelled (e.g., by the region of their primary emission) and then processed as the species itself. The type-2 family concept is designed to conserve the linear relationship between the family and its members. The second submodel is TRACER_PDEF, which corrects and budgets numerical negative overshoots that arise in many process implementations due to the numerical limitations (e.g., rounding errors). The submodel therefore guarantees the positive definiteness of the tracers and stabilises the integration scheme. As a by-product, it further provides a global tracer mass diagnostic. Last but not least, we present the submodel PTRAC, which allows the definition of tracers via a Fortran95 namelist, as a complement to the standard tracer definition by application of the TRACER interface routines in the code. TRACER with its submodels and PTRAC can readily be applied to a variety of models without further requirements. The code and a documentation are included in the electronic supplement

A refined method for calculating equivalent effective stratospheric chlorine

Chlorine and bromine atoms lead to catalytic depletion of ozone in the stratosphere. Therefore the use and production of ozone-depleting substances (ODSs) containing chlorine and bromine is regulated by the Montreal Protocol to protect the ozone layer. Equivalent effective stratospheric chlorine (EESC) has been adopted as an appropriate metric to describe the combined effects of chlorine and bromine released from halocarbons on stratospheric ozone. Here we revisit the concept of calculating EESC. We derive a refined formulation of EESC based on an advanced concept of ODS propagation into the stratosphere and reactive halogen release. A new transit time distribution is introduced in which the age spectrum for an inert tracer is weighted with the release function for inorganic halogen from the source gases. This distribution is termed the "release time distribution". We show that a much better agreement with inorganic halogen loading from the chemistry transport model TOMCAT is achieved compared with using the current formulation. The refined formulation shows EESC levels in the year 1980 for the mid-latitude lower stratosphere, which are significantly lower than previously calculated. The year 1980 is commonly used as a benchmark to which EESC must return in order to reach significant progress towards halogen and ozone recovery. Assuming that – under otherwise unchanged conditions – the EESC value must return to the same level in order for ozone to fully recover, we show that it will take more than 10 years longer than estimated in this region of the stratosphere with the current method for calculation of EESC. We also present a range of sensitivity studies to investigate the effect of changes and uncertainties in the fractional release factors and in the assumptions on the shape of the release time distributions. We further discuss the value of EESC as a proxy for future evolution of inorganic halogen loading under changing atmospheric dynamics using simulations from the EMAC model. We show that while the expected changes in stratospheric transport lead to significant differences between EESC and modelled inorganic halogen loading at constant mean age, EESC is a reasonable proxy for modelled inorganic halogen on a constant pressure level

Simulating organic species with the global atmospheric chemistry general circulation model ECHAM5/MESSy1: a comparison of model results with observations

The atmospheric-chemistry general circulation model ECHAM5/MESSy1 is evaluated with observations of different organic ozone precursors. This study continues a prior analysis which focused primarily on the representation of atmospheric dynamics and ozone. We use the results of the same reference simulation and apply a statistical analysis using data from numerous field campaigns. The results serve as a basis for future improvements of the model system. ECHAM5/MESSy1 generally reproduces the spatial distribution and the seasonal cycle of carbon monoxide (CO) very well. However, for the background in the Northern Hemisphere we obtain a negative bias (mainly due to an underestimation of emissions from fossil fuel combustion), and in the high latitude Southern Hemisphere a yet unexplained positive bias. The model results agree well with observations of alkanes, whereas severe problems in the simulation of alkenes and isoprene are present. For oxygenated compounds the results are ambiguous: The model results are in good agreement with observations of formaldehyde, but systematic biases are present for methanol and acetone. The discrepancies between the model results and the observations are explained (partly) by means of sensitivity studies