Significance: We assessed the effects of air pollution and greenhouse gases on public health, climate, and the hydrologic cycle. We combined a global atmospheric chemistry–climate model with air pollution exposure functions, based on an unmatched large number of cohort studies in many countries. We find that fossil-fuel-related emissions account for about 65% of the excess mortality rate attributable to air pollution, and 70% of the climate cooling by anthropogenic aerosols. We conclude that to save millions of lives and restore aerosol-perturbed rainfall patterns, while limiting global warming to 2 °C, a rapid phaseout of fossil-fuel-related emissions and major reductions of other anthropogenic sources are needed.

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

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

The self-cleaning or oxidation capacity of the atmosphere is principally controlled by hydroxyl (OH) radicals in the troposphere. Hydroxyl has primary (P) and secondary (S) sources, the former mainly through the photodissociation of ozone, the latter through OH recycling in radical reaction chains. We used the recent Mainz Organics Mechanism (MOM) to advance volatile organic carbon (VOC) chemistry in the general circulation model EMAC (ECHAM/MESSy Atmospheric Chemistry) and show that S is larger than previously assumed. By including emissions of a large number of primary VOC, and accounting for their complete breakdown and intermediate products, MOM is mass-conserving and calculates substantially higher OH reactivity from VOC oxidation compared to predecessor models. Whereas previously P and S were found to be of similar magnitude, the present work indicates that S may be twice as large, mostly due to OH recycling in the free troposphere. Further, we find that nighttime OH formation may be significant in the polluted subtropical boundary layer in summer. With a mean OH recycling probability of about 67 %, global OH is buffered and not sensitive to perturbations by natural or anthropogenic emission changes. Complementary primary and secondary OH formation mechanisms in pristine and polluted environments in the continental and marine troposphere, connected through long-range transport of O₃, can maintain stable global OH levels.

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 (1950–2011), hindcast simulations with specified dynamics (1979–2013), i.e. nudged towards ERA-Interim reanalysis data, and combined hindcast and projection simulations (1950–2100). 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 set-ups, 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.