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.

[1] Ozone changes and associated climate impacts in the Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations are analyzed over the historical and future (2006-2100) period under four Representative Concentration Pathways (RCP). In contrast to CMIP3, where half of the models prescribed constant stratospheric ozone, CMIP5 models all consider past ozone depletion and future ozone recovery. Multimodel mean climatologies and long-term changes in total and tropospheric column ozone calculated from CMIP5 models with either interactive or prescribed ozone are in reasonable agreement with observations. However, some large deviations from observations exist for individual models with interactive chemistry, and these models are excluded in the projections. Stratospheric ozone projections forced with a single halogen, but four greenhouse gas (GHG) scenarios show largest differences in the northern midlatitudes and in the Arctic in spring (~20 and 40 Dobson units (DU) by 2100, respectively). By 2050, these differences are much smaller and negligible over Antarctica in austral spring. Differences in future tropospheric column ozone are mainly caused by differences in methane concentrations and stratospheric input, leading to~10 DU increases compared to 2000 in RCP 8.5. Large variations in stratospheric ozone particularly in CMIP5 models with interactive chemistry drive correspondingly large variations in lower stratospheric temperature trends. The results also illustrate that future Southern Hemisphere summertime circulation changes are controlled by both the ozone recovery rate and the rate of GHG increases, emphasizing the importance of simulating and taking into account ozone forcings when examining future climate projections.

Abstract. Lineshaped contrails were detected with the research aircraft Falcon during the CONCERT -CONtrail and Cirrus ExpeRimenT -campaign in October/November 2008. The Falcon was equipped with a set of instruments to measure the particle size distribution, shape, extinction and chemical composition as well as trace gas mixing ratios of sulfur dioxide (SO 2 ), reactive nitrogen and halogen species (NO, NO y , HNO 3 , HONO, HCl), ozone (O 3 ) and carbon monoxide (CO). During 12 mission flights over Europe, numerous contrails, cirrus clouds and a volcanic aerosol layer were probed at altitudes between 8.5 and 11.6 km and at temperatures above 213 K. 22 contrails from 11 different aircraft were observed near and below ice saturation. The observed NO mixing ratios, ice crystal and soot number densities are compared to a process based contrail model. On 19 November 2008 the contrail from a CRJ-2 aircraft was penetrated in 10.1 km altitude at a temperature of 221 K. The contrail had mean ice crystal number densities of 125 cm −3 with effective radii r eff of 2.6 µm. The presence of particles with r>50 µm in the less than 2 min old contrail suggests that natural cirrus crystals were entrained in the contrail. Mean HONO/NO (HONO/NO y ) ratios of 0.037 (0.024) and the fuel sulfur conCorrespondence to: C. Voigt (christiane.voigt@dlr.de) version efficiency to H 2 SO 4 ( S↓ ) of 2.9 % observed in the CRJ-2 contrail are in the range of previous measurements in the gaseous aircraft exhaust. On 31 October 2010 aviation NO emissions could have contributed by more than 40% to the regional scale NO levels in the mid-latitude lowest stratosphere. The CONCERT observations help to better quantify the climate impact from contrails and will be used to investigate the chemical processing of trace gases on contrails.

Abstract.A community diagnostics and performance metrics tool for the evaluation of Earth system models (ESMs) has been developed that allows for routine comparison of single or multiple models, either against predecessor versions or against observations. The priority of the effort so far has been to target specific scientific themes focusing on selected essential climate variables (ECVs), a range of known systematic biases common to ESMs, such as coupled tropical climate variability, monsoons, Southern Ocean processes, continental dry biases, and soil hydrology-climate interactions, as well as atmospheric CO 2 budgets, tropospheric and stratospheric ozone, and tropospheric aerosols. The tool is being developed in such a way that additional analyses can easily be added. A set of standard namelists for each scientific topic reproduces specific sets of diagnostics or performance metrics that have demonstrated their importance in ESM evaluation in the peer-reviewed literature. The Earth System Model Evaluation Tool (ESMValTool) is a community effort open to both users and developers encouraging open exchange of diagnostic source code and evaluation results from the Coupled Model Intercomparison Project (CMIP) ensemble. This will facilitate and improve ESM evaluation beyond the stateof-the-art and aims at supporting such activities within CMIP and at individual modelling centres. Ultimately, we envisage running the ESMValTool alongside the Earth System Grid Federation (ESGF) as part of a more routine evaluation of CMIP model simulations while utilizing observations available in standard formats (obs4MIPs) or provided by the user.

Abstract.A quasi chemistry-transport model mode (QCTM) is presented for the numerical chemistry-climate simulation system ECHAM/MESSy Atmospheric Chemistry (EMAC). It allows for a quantification of chemical signals through suppression of any feedback between chemistry and dynamics. Noise would otherwise interfere too strongly. The signal is calculated from the difference of two QCTM simulations, a reference simulation and a sensitivity simulation. In order to avoid the feedbacks, the simulations adopt the following offline chemical fields: (a) offline mixing ratios of radiatively active substances enter the radiation scheme, (b) offline mixing ratios of nitric acid enter the scheme for re-partitioning and sedimentation from polar stratospheric clouds, (c) and offline methane oxidation is the exclusive source of chemical water-vapor tendencies. Any set of offline fields suffices to suppress the feedbacks, though may be inconsistent with the simulation setup. An adequate set of offline climatologies can be produced from a non-QCTM simulation using the setup of the reference simulation. Test simulations reveal the particular importance of adequate offline fields associated with (a). Inconsistencies from (b) are negligible when using adequate fields of nitric acid. Acceptably small inconsistencies come from (c), but should vanish for an adequate prescription of chemical water vapor tendencies. Toggling between QCTM and non-QCTM is done via namelist switches and does not require a source code re-compilation.

Abstract. A community diagnostics and performance metrics tool for the evaluation of Earth System Models (ESMs) has been developed that allows for routine comparison of single or multiple models, either against predecessor versions or against observations. The priority of the effort so far has been to target specific scientific themes focusing on selected Essential Climate Variables (ECVs), a range of known systematic biases common to ESMs, such as coupled tropical climate variability, monsoons, Southern Ocean processes, continental dry biases and soil hydrology-climate interactions, as well as atmospheric CO2 budgets, tropospheric and stratospheric ozone, and tropospheric aerosols. The tool is being developed in such a way that additional analyses can easily be added. A set of standard namelists for each scientific topic reproduces specific sets of diagnostics or performance metrics that have demonstrated their importance in ESM evaluation in the peer-reviewed literature. The Earth System Model Evaluation Tool (ESMValTool) is a community effort open to both users and developers encouraging open exchange of diagnostic source code and evaluation results from the CMIP ensemble. This will facilitate and improve ESM evaluation beyond the state-of-the-art and aims at supporting such activities within the Coupled Model Intercomparison Project (CMIP) and at individual modelling centres. Ultimately, we envisage running the ESMValTool alongside the Earth System Grid Federation (ESGF) as part of a more routine evaluation of CMIP model simulations while utilizing observations available in standard formats (obs4MIPs) or provided by the user.

Abstract. Four simulations with the ECHAM/MESSy Atmospheric Chemistry (EMAC) model have been evaluated with the Earth System Model Validation Tool (ESMValTool) to identify differences in simulated ozone and selected climate parameters that resulted from (i) different setups of the EMAC model (nudged vs. free-running) and (ii) different boundary conditions (emissions, sea surface temperatures (SSTs) and sea ice concentrations (SICs)). To assess the relative performance of the simulations, quantitative performance metrics are calculated consistently for the climate parameters and ozone. This is important for the interpretation of the evaluation results since biases in climate can impact on biases in chemistry and vice versa. The observational data sets used for the evaluation include ozonesonde and aircraft data, meteorological reanalyses and satellite measurements. The results from a previous EMAC evaluation of a model simulation with nudging towards realistic meteorology in the troposphere have been compared to new simulations with different model setups and updated emission data sets in free-running time slice and nudged quasi chemistrytransport model (QCTM) mode. The latter two configurations are particularly important for chemistry-climate projections and for the quantification of individual sources (e.g., the transport sector) that lead to small chemical perturbations of the climate system, respectively. With the exception of some specific features which are detailed in this study, no large differences that could be related to the different setups (nudged vs. free-running) of the EMAC simulations were found, which offers the possibility to evaluate and improve the overall model with the help of shorter nudged simulations. The main differences between the two setups is a better representation of the tropospheric and stratospheric temperature in the nudged simulations, which also better reproduce stratospheric water vapor concentrations, due to the improved simulation of the temperature in the tropical tropopause layer. Ozone and ozone precursor concentrations, on the other hand, are very similar in the different model setups, if similar boundary conditions are used. Different boundary conditions however lead to relevant differences in the four simulations. Biases which are common to all simulations are the underestimation of the ozone hole and the overestimation of tropospheric column ozone, the latter being significantly reduced when lower lightning emissions of nitrogen oxides are used. To further investigate possible other reasons for such bias, two sensitivity simulations with an updated scavenging routine and the addition of a newly proposed HNO 3 -forming channel of the HO 2 + NO reaction were performed. The update in the scavenging routine resulted in a slightly better representation of ozone compared to the reference simulation. The introduction of the new HNO 3 -forming channel significantly reduces the overestimation of tropospheric ozone. Therefore, including the new reaction rate could potentially...

Tear resistance at the edge of a slab is an important parameter controlling the evolution of subduction zones. However, compared with other subduction parameters such as plate strength, plate viscosity, plate thickness and trench width, the dynamics of tearing are poorly understood. Here we obtain a first-order understanding of the dynamics and morphology of subduction zones to resistance during tear propagation, by developing and using a novel computational modelling technique for subducting slabs, with side boundaries described by visco-plastic weak zones, developing into tear faults. Our 3D model is based upon a visco-plastic slab that sinks into the less dense mantle, generating poloidal and toroidal flows. The asthenospheric mantle field is static and only develops flow due to the subducting slab. We use the finite element code eScript/Finley and the level set method to describe the lithosphere to solve this fluid dynamics problem. Our results show the importance of tear resistance for the speed of trench migration and for shaping the final geometry of subduction systems. We show that slab tearing along a weak layer can result in a relatively straight slab hinge shape, while increasing the strength in the weak layer results in the curvature of the hinge increasing substantially. High tear resistance at the slab edges may hinder rollback to the extent that the slab becomes stretched and recumbently folded at the base of the domain. Tear resistance also controls whether the subducting lithosphere can experience accelerating rollback velocities or a constant rollback velocity.