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.

This paper evaluates the current status of global modeling of the organic aerosol (OA) in the troposphere and analyzes the differences between models as well as between models and observations. Thirty-one global chemistry transport models (CTMs) and general circulation models (GCMs) have participated in this intercomparison, in the framework of AeroCom phase II. The simulation of OA varies greatly between models in terms of the magnitude of primary emissions, secondary OA (SOA) formation, the number of OA species used (2 to 62), the complexity of OA parameterizations (gas-particle partitioning, chemical aging, multiphase chemistry, aerosol microphysics), and the OA physical, chemical and optical properties. The diversity of the global OA simulation results has increased since earlier AeroCom experiments, mainly due to the increasing complexity of the SOA parameterization in models, and the implementation of new, highly uncertain, OA sources. Diversity of over one order of magnitude exists in the modeled vertical distribution of OA concentrations that deserves a dedicated future study. Furthermore, although the OA/OC ratio depends on OA sources and atmospheric processing, and is important for model evaluation against OA and OC observations, it is resolved only by a few global models. The median global primary OA (POA) source strength is 56 Tg a(-1) (range 34-144 Tg a(-1)) and the median SOA source strength (natural and anthropogenic) is 19 Tg a(-1) (range 13-121 Tg a(-1)). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a(-1) (range 16-121 Tg a(-1)), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a(-1); range 13-20 Tg a(-1), with one model at 37 Tg a(-1)). The median atmospheric burden of OA is 1.4 Tg (24 models in the range of 0.6-2.0 Tg and 4 between 2.0 and 3.8 Tg), with a median OA lifetime of 5.4 days (range 3.8-9.6 days). In models that reported both OA and sulfate burdens, the median value of the OA/sulfate burden ratio is calculated to be 0.77; 13 models calculate a ratio lower than 1, and 9 models higher than 1. For 26 models that reported OA deposition fluxes, the median wet removal is 70 Tg a(-1) (range 28-209 Tg a(-1)), which is on average 85% of the total OA deposition. Fine aerosol organic carbon (OC) and OA observations from continuous monitoring networks and individual field campaigns have been used for model evaluation. At urban locations, the model-observation comparison indicates missing knowledge on anthropogenic OA sources, both strength and seasonality. The combined model-measurements analysis suggests the existence of increased OA levels during summer due to biogenic SOA formation over large areas of the USA that can be of the same order of magnitude as the POA, even at urban locations, and contribute to the measured urban seasonal pattern. Global models are able to simulate the high secondary character of OA observed in the atmosphere as a result of SOA formation and POA agi...

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.