The Modular Earth Submodel System (MESSy) is an open, multi-institutional project providing a strategy for developing comprehensive Earth System Models (ESMs) with highly flexible complexity. The first version of the MESSy infrastructure and process submodels, mainly focusing on atmospheric chemistry, has been successfully coupled to an atmospheric General Circulation Model (GCM) expanding it into an Atmospheric Chemistry GCM (AC-GCM) for nudged simulations and into a Chemistry Climate Model (CCM) for climate simulations. <br><br> Here, we present the second development cycle of MESSy, which comprises (1) an improved and extended infrastructure for the basemodel independent coupling of process-submodels, (2) new, highly valuable diagnostic capabilities for the evaluation with observational data and (3) an improved atmospheric chemistry setup. With the infrastructural changes, we place the headstone for further model extensions from a CCM towards a comprehensive ESM. The new diagnostic submodels will be used for regular re-evaluations of the continuously further developing model system. The updates of the chemistry setup are briefly evaluated

Abstract.The new Modular Earth Submodel System (MESSy) describes atmospheric chemistry and meteorological processes in a modular framework, following strict coding standards. It has been coupled to the ECHAM5 general circulation model, which has been slightly modified for this purpose. A 90-layer model setup up to 0.01 hPa was used at spectral T42 resolution to simulate the lower and middle atmosphere. With the high vertical resolution the model simulates the Quasi-Biennial Oscillation. The model meteorology has been tested to check the influence of the changes to ECHAM5 and the radiation interactions with the new representation of atmospheric composition. In the simulations presented here a Newtonian relaxation technique was applied in the tropospheric part of the domain to weakly nudge the model towards the analysed meteorology during the period 1998-2005. This allows an efficient and direct evaluation with satellite and in-situ data. It is shown that the tropospheric wave forcing of the stratosphere in the model suffices to reproduce major stratospheric warming events leading e.g. to the vortex split over Antarctica in 2002. Characteristic features such as dehydration and denitrification caused by the sedimentation of polar stratospheric cloud particles and ozone depletion during winter and spring are simulated well, although ozone loss in the lower polar stratosphere is slightly underestimated. The model realistically simulates stratosphere-troposphere exchange processes as indicated by comparisons with satellite and in situ measurements. The evaluation of tropospheric chemistry presented here focuses on the distributions of ozone, hydroxyl radicals, carbon monoxide and reactive nitrogen compounds. In spite of minor shortcomings, mostly related to the relatively Correspondence to: P. Jöckel (joeckel@mpch-mainz.mpg.de) coarse T42 resolution and the neglect of inter-annual changes in biomass burning emissions, the main characteristics of the trace gas distributions are generally reproduced well. The MESSy submodels and the ECHAM5/MESSy1 model output are available through the internet on request.

Abstract. Three types of reference simulations, as recommended by the Chemistry-Climate Model Initiative (CCMI), have been performed with version 2.51 of the European Centre for Medium-Range Weather Forecasts -Hamburg (ECHAM)/Modular Earth Submodel System (MESSy) Atmospheric Chemistry (EMAC) model: hindcast simulations , hindcast simulations with specified dynamics , i.e. nudged towards ERA-Interim reanalysis data, and combined hindcast and projection simulations . The manuscript summarizes the updates of the model system and details the different model set-ups used, including the on-line calculated diagnostics. Simulations have been performed with two different nudging setups, with and without interactive tropospheric aerosol, and with and without a coupled ocean model. Two different vertical resolutions have been applied. The on-line calculated sources and sinks of reactive species are quantified and a first evaluation of the simulation results from a global perspective is provided as a quality check of the data. The focus is on the intercomparison of the different model set-ups. The simulation data will become publicly available via CCMI and the Climate and Environmental Retrieval and Archive (CERA) database of the German Climate Computing Centre (DKRZ). This manuscript is intended to serve as an extensive reference for further analyses of the Earth System Chemistry integrated Modelling (ESCiMo) simulations.

Abstract. Bacteria are constantly being transported through the atmosphere, which may have implications for human health, agriculture, cloud formation, and the dispersal of bacterial species. We simulate the global transport of bacteria, represented as 1 μm and 3 μm diameter spherical solid particle tracers in a general circulation model. We investigate factors influencing residence time and distribution of the particles, including emission region, cloud condensation nucleus activity and removal by ice-phase precipitation. The global distribution depends strongly on the assumptions made about uptake into cloud droplets and ice. The transport is also affected, to a lesser extent, by the emission region, particulate diameter, and season. We find that the seasonal variation in atmospheric residence time is insufficient to explain by itself the observed seasonal variation in concentrations of particulate airborne culturable bacteria, indicating that this variability is mainly driven by seasonal variations in culturability and/or emission strength. We examine the potential for exchange of bacteria between ecosystems and obtain rough estimates of the flux from each ecosystem by using a maximum likelihood estimation technique, together with a new compilation of available observations described in a companion paper. Globally, we estimate the total emissions of bacteria-containing particles to the atmosphere to be 7.6×1023–3.5×1024 a−1, originating mainly from grasslands, shrubs and crops. We estimate the mass of emitted bacteria- to be 40–1800 Gg a−1, depending on the mass fraction of bacterial cells in the particles. In order to improve understanding of this topic, more measurements of the bacterial content of the air and of the rate of surface-atmosphere exchange of bacteria will be necessary. Future observations in wetlands, hot deserts, tundra, remote glacial and coastal regions and over oceans will be of particular interest.

Abstract. The new Modular Earth Submodel System (MESSy) describes atmospheric chemistry and meteorological processes in a modular framework, following strict coding standards. It has been coupled to the ECHAM5 general circulation model, which has been slightly modified for this purpose. A 90-layer model version up to 0.01 hPa was used at T42 resolution (~2.8 latitude and longitude) to simulate the lower and middle atmosphere. The model meteorology has been tested to check the influence of the changes to ECHAM5 and the radiation interactions with the new representation of atmospheric composition. A Newtonian relaxation technique was applied in the tropospheric part of the domain to weakly nudge the model towards the analysed meteorology during the period 1998–2005. It is shown that the tropospheric wave forcing of the stratosphere in the model suffices to reproduce the Quasi-Biennial Oscillation and major stratospheric warming events leading e.g. to the vortex split over Antarctica in 2002. Characteristic features such as dehydration and denitrification caused by the sedimentation of polar stratospheric cloud particles and ozone depletion during winter and spring are simulated accurately, although ozone loss in the lower polar stratosphere is slightly underestimated. The model realistically simulates stratosphere-troposphere exchange processes as indicated by comparisons with satellite and in situ measurements. The evaluation of tropospheric chemistry presented here focuses on the distributions of ozone, hydroxyl radicals, carbon monoxide and reactive nitrogen compounds. In spite of minor shortcomings, mostly related to the relatively coarse T42 resolution and the neglect of interannual changes in biomass burning emissions, the main characteristics of the trace gas distributions are generally reproduced well. The MESSy submodels and the ECHAM5/MESSy1 model output are available through the internet on request.

Abstract. We 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.

Abstract. The development of a comprehensive Earth System Model (ESM) to study the interactions between chemical, physical, and biological processes, requires coupling of the different domains (land, ocean, atmosphere, ...). One strategy is to link existing domain-specific models with a universal coupler, i.e. an independent standalone program organizing the communication between other programs. In many cases, however, a much simpler approach is more feasible. We have developed the Modular Earth Submodel System (MESSy). It comprises (1) a modular interface structure to connect to a , (2) an extendable set of such for miscellaneous processes, and (3) a coding standard. MESSy is therefore not a coupler in the classical sense, but exchanges data between a and several within one comprehensive executable. The internal complexity of the is controllable in a transparent and user friendly way. This provides remarkable new possibilities to study feedback mechanisms (by two-way coupling). Note that the MESSy and the coupler approach can be combined. For instance, an atmospheric model implemented according to the MESSy standard could easily be coupled to an ocean model by means of an external coupler. The vision is to ultimately form a comprehensive ESM which includes a large set of submodels, and a base model which contains only a central clock and runtime control. This can be reached stepwise, since each process can be included independently. Starting from an existing model, process submodels can be reimplemented according to the MESSy standard. This procedure guarantees the availability of a state-of-the-art model for scientific applications at any time of the development. In principle, MESSy can be implemented into any kind of model, either global or regional. So far, the MESSy concept has been applied to the general circulation model ECHAM5 and a number of process boxmodels.

We present version 3.0 of the atmospheric chemistry box model CAABA/MECCA. In addition to a complete update of the rate coefficients to the most recent recommendations, a number of new features have been added: chemistry in multiple aerosol size bins; automatic multiple simulations reaching steady-state conditions; Monte-Carlo simulations with randomly varied rate coefficients within their experimental uncertainties; calculations along Lagrangian trajectories; mercury chemistry; more detailed isoprene chemistry; tagging of isotopically labeled species. Further changes have been implemented to make the code more user-friendly and to facilitate the analysis of the model results. Like earlier versions, CAABA/MECCA-3.0 is a community model published under the GNU General Public License