Search

This study uses a chemistry–climate model CHASER (MIROC) to explore the roles of heterogeneous reactions (HRs) in global tropospheric chemistry. Three distinct HRs of N2O5, HO2, and RO2 are considered for surfaces of aerosols and cloud particles. The model simulation is verified with EANET and EMEP stationary observations; R/V Mirai ship-based data; ATom1 aircraft measurements; satellite observations by OMI, ISCCP, and CALIPSO-GOCCP; and reanalysis data JRA55. The heterogeneous chemistry facilitates improvement of model performance with respect to observations for NO2, OH, CO, and O3, especially in the lower troposphere. The calculated effects of heterogeneous reactions cause marked changes in global abundances of O3 (−2.96 %), NOx (−2.19 %), CO (+3.28 %), and global mean CH4 lifetime (+5.91 %). These global effects were contributed mostly by N2O5 uptake onto aerosols in the middle troposphere. At the surface, HO2 uptake gives the largest contributions, with a particularly significant effect in the North Pacific region (−24 % O3, +68 % NOx, +8 % CO, and −70 % OH), mainly attributable to its uptake onto clouds. The RO2 reaction has a small contribution, but its global mean negative effects on O3 and CO are not negligible. In general, the uptakes onto ice crystals and cloud droplets that occur mainly by HO2 and RO2 radicals cause smaller global effects than the aerosol-uptake effects by N2O5 radicals (+1.34 % CH4 lifetime, +1.71 % NOx, −0.56 % O3, +0.63 % CO abundances). Nonlinear responses of tropospheric O3, NOx, and OH to the N2O5 and HO2 uptakes are found in the same modeling framework of this study (R>0.93). Although all HRs showed negative tendencies for OH and O3 levels, the effects of HR(HO2) on the tropospheric abundance of O3 showed a small increment with an increasing loss rate. However, this positive tendency turns to reduction at higher rates (>5 times). Our results demonstrate that the HRs affect not only polluted areas but also remote areas such as the mid-latitude sea boundary layer and upper troposphere. Furthermore, HR(HO2) can bring challenges to pollution reduction efforts because it causes opposite effects between NOx (increase) and surface O3 (decrease).

Ensembles of chemistry-climate models (CCMs) are fundamental for the exploration of the chemistry-climate system. A particular focus of chemistry-climate modelling is stratospheric ozone, whose concentrations have been decreased by anthropogenic releases of ozone depleting substances. In conjunction with observational data, CCM ensembles have been relied upon to simulate historic effects of ozone depletion and to project future ozone recovery. However, many widely used ensemble analysis methods are simplistic and are based upon incorrect assumptions about the design of the ensemble. Multi-model means used to construct future ozone projections do not account for variable model performance or similarity and therefore give biased and inaccurate projections. Similarly, simplistic linear regression methods used to infill historic ozone records underestimate interannual variability and are inaccurate in regions of sparse data coverage. Moreover, given advances in machine learning and data science and their increased use in environmental science, it is timely to apply more advanced tools to CCM ensembles. To address this methodological deficit, this thesis presents a set of novel tools to improve the predictions and projections from CCM ensembles of stratospheric ozone. A process-based weighted mean is developed which accounts for model performance and similarity in CCM ensembles. This improvement over pre-existing methods was used to generate accurate ozone hole recovery projections. This thesis also developed a Bayesian neural network (BNN) which fuses together CCMs with observational data to produce accurate and uncertainty-aware predictions. The BNN framework was used to produce historic continuous datasets of total ozone column and vertically resolved ozone, and represents a significant improvement in methods used to ensemble models. Though designed for CCM ensembles these flexible tools have the potential to be applied to other environmental modelling disciplines to improve the accuracy of projections, better understand uncertainty and to make better use of historic observations.

This study developed a next‐generation atmospheric aerosol/chemistry‐climate model, the BCC_AGCM_CUACE2.0. Then, the performance of the model for nitrate was evaluated, and the nitrate direct radiative forcing (DRF) and effective radiative forcing (ERF) due to aerosol‐radiation interactions were simulated for the present day (2010), near‐term future (2030), and middle‐term future (2050) under the Representative Concentration Pathway 4.5, 6.0, and 8.5 scenarios relative to the preindustrial era (1850). The model reproduced the distributions and seasonal changes in nitrate loading well, and simulated surface concentrations matched observations in Europe, North America, and China. Current global mean annual loading of nitrates was predicted to increase by 1.50 mg m−2 relative to 1850, with the largest increases occurring in East Asia (9.44 mg m−2), Europe (4.36 mg m−2), and South Asia (3.09 mg m−2). The current global mean annual ERF of nitrates was −0.28 W m−2 relative to 1850. Due to global reductions in pollutant emissions, the nitrate ERF values were predicted to decrease to −0.17, −0.20, and −0.24 W m−2 in 2030 and −0.07, −0.18, and −0.19 W m−2 in 2050 for Representative Concentration Pathway 4.5, 6.0, and 8.5 relative to 1850, respectively. Although global mean nitrate values showed a declining trend, future nitrate loading remained high in East Asia and South Asia.

This study describes the simulation characteristics of a newly developed Global/Regional Integrated Model system‐Chemistry Climate Model (GRIMs‐CCM), which is listed in Chemistry Climate Model Initiative as a participating model. The GRIMs‐CCM was run using a standard set of forcings, and historical sea surface temperatures and sea ice concentration from the Hadley Centre were prescribed. The simulation results of GRIMs‐CCM were compared to the National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis data set 1. The GRIMs‐CCM satisfactorily simulated the climatological (1960–2010) atmospheric features and atmospheric teleconnections resulting from tropical sea surface temperature forcing. However, the GRIMs‐CCM also had some regional biases; for instance, particularly, the temperature bias over the Antarctic was noticeable. We further analyzed physical processes that caused such biases and the influence of coupled chemistry‐climate processes in the GRIMs‐CCM, which may provide further guidance to improve an earlier version of the GRIMs‐CCM and other climate‐chemistry models participating in the Chemistry Climate Model Initiative.

With the gradual yet unequivocal phasing out of ozone depleting substances (ODSs), the environmental crisis caused by the discovery of an ozone hole over the Antarctic has lessened in severity and a promising recovery of the ozone layer is predicted in this century. However, strong volcanic activity can also cause ozone depletion that might be severe enough to threaten the existence of life on Earth. In this study, a transport model and a coupled chemistry-climate model were used to simulate the impacts of super volcanoes on ozone depletion. The volcanic eruptions in the experiments were the 1991 Mount Pinatubo eruption and a 100 × Pinatubo size eruption. The results show that the percentage of global mean total column ozone depletion in the 2050 RCP8.5 100 × Pinatubo scenario is approximately 6% compared to two years before the eruption and 6.4% in tropics. An identical simulation, 100 × Pinatubo eruption only with natural source ODSs, produces an ozone depletion of 2.5% compared to two years before the eruption, and with 4.4% loss in the tropics. Based on the model results, the reduced ODSs and stratospheric cooling lighten the ozone depletion after super volcanic eruption.

Three different methods of specifying ozone in an atmosphere‐only version of the HadGEM3‐A global circulation model are compared to the coupled chemistry configuration of this model. These methods include a specified zonal‐mean ozone climatology, a specified 3‐D ozone climatology, and a calculated‐asymmetry scheme in which a specified zonal‐mean ozone field is adapted online to be consistent with dynamically produced zonal asymmetries. These simulations all use identical boundary conditions and, by construction, have the same climatological zonal‐mean ozone, that of the coupled chemistry configuration of the model. Prescribing ozone, regardless of scheme, results in a simulation which is 3–4 times faster than the coupled chemistry‐climate model (CCM). Prescribing climatological zonal asymmetries leads to a vortex which is the correct intensity but which is systematically displaced over regions with lower prescribed ozone. When zonal asymmetries in ozone are free to evolve interactively with model dynamics, the modeled wintertime stratospheric vortex shape and mean sea level pressure patterns closely resemble that produced by the full CCM in both hemispheres, in terms of statistically significant differences. Further, we separate out the two distinct pathways by which zonal ozone asymmetries influence modeled dynamics. We present this interactive‐ozone zonal‐asymmetry scheme as an inexpensive tool for accurately modeling the impacts of dynamically consistent ozone fields as seen in a CCM which ultimately influence mean sea level pressure and tropospheric circulation (particularly during wintertime in the Northern Hemisphere, when ozone asymmetries are generally largest), without the computational burden of simulating interactive chemistry., Key Points Dynamical feedbacks resulting from ozone zonal asymmetries are important for accurately representing Northern Hemisphere circulation patterns during wintertimeThree simulations using different methods of specifying ozone mixing ratios are compared with a fully coupled CCM simulation in order to examine the pathways by which zonal ozone asymmetries influence modeled dynamicsDynamical feedbacks occurring in a coupled CCM can be reproduced, at a greatly reduced computational cost, by appropriately placing a prescribed ozone climatology onto a dynamically evolving coordinate

In the extratropical atmosphere, Rossby waves (RWs) and internal gravity waves (GWs) propagating from the troposphere mediate a coupling with the middle atmosphere by influencing the dynamics herein. In the current generation chemistry-climate models (CCMs), RW effects are well resolved while GW effects have to be parameterized. Here, we analyze orographic GW (OGW) interaction with resolved dynamics in a comprehensive CCM on the time scale of days. For this, we apply a recently developed method of strong OGW drag event composites for the three strongest northern hemisphere OGW hotspots. We show that locally-strong OGW events considerably alter the properties of resolved wave propagation into the middle atmosphere, which subsequently influences zonal winds and RW transience. Our results demonstrate that the influence of OGWs is critically dependent on the hotspot region, which underlines the OGW-resolved dynamics interaction being a two-way process.

Ozone feedback effects on the quasi-biennial oscillation (QBO) were investigated with a chemistry–climate model (CCM) by modifying ozone abundance in the radiative process. Under a standard run for 50 years, the CCM could realistically reproduce the QBO of about a 28-month period for wind and ozone. Five experiment runs were made for 20 years through varying ozone abundance only in the equatorial stratosphere from 100 to 10 hPa by −40, −20, −10, +10, and +20%, respectively, after the chemistry module and transferring the resultant ozone to the radiation calculation. It was found that the modification of ozone abundance in the radiation substantially changed the period of the QBO but slightly influenced the amplitude of the QBO. The 10% and 20% increase runs led to longer QBO periods (31 and 34 months) than that of the standard run, i.e., lengthening by 3 and 6 months, while the 10%, 20%, and 40% decrease runs resulted in shorter periods (24, 22, and 17 months), i.e., shortening by 4, 6, and 11 months. These substantial changes in the QBO period in the experiment runs indicate that the ozone feedback significantly affects the QBO dynamics through the modulation in solar heating.

Arctic ozone amount in winter to spring shows large year-to-year variation. This study investigates Arctic spring ozone in relation to the phase of quasi-biennial oscillation (QBO)/the 11-year solar cycle, using satellite observations, reanalysis data, and outputs of a chemistry climate model (CCM) during the period of 1979–2011. For this duration, we found that the composite mean of the Northern Hemisphere high-latitude total ozone in the QBO-westerly (QBO-W)/solar minimum (Smin) phase is slightly smaller than those averaged for the QBO-W/Smax and QBO-E/Smax years in March. An analysis of a passive ozone tracer in the CCM simulation indicates that this negative anomaly is primarily caused by transport. The negative anomaly is consistent with a weakening of the residual mean downward motion in the polar lower stratosphere. The contribution of chemical processes estimated using the column amount difference between ozone and the passive ozone tracer is between 10–20% of the total anomaly in March. The lower ozone levels in the Arctic spring during the QBO-W/Smin years are associated with a stronger Arctic polar vortex from late winter to early spring, which is linked to the reduced occurrence of sudden stratospheric warming in the winter during the QBO-W/Smin years.

Reconstructions of volcanic aerosol forcing and its climatic impacts are undermined by uncertainties in both the models used to build these reconstructions as well as the proxy and observational records used to constrain those models. Reducing these uncertainties has been a priority and in particular, several modelling groups have developed interactive stratospheric aerosol models. Provided with an initial volcanic injection of sulfur dioxide, these models can interactively simulate the life cycle and optical properties of sulfate aerosols, and their effects on climate. In contrast, most climate models that took part in the Coupled Model Intercomparison Project Phase 5 and 6 (CMIP6) directly prescribe perturbations in atmospheric optical properties associated with an eruption. However, before the satellite era, the volcanic forcing dataset used for CMIP6 mostly relies on a relatively simple aerosol model and a volcanic sulfur inventory derived from ice-cores, both of which have substantial associated uncertainties.In this study, we produced a new set of historical simulations using the UK Earth System Model UKESM1, with interactive stratospheric aerosol capability (referred to as interactive runs hereafter) instead of directly prescribing the CMIP6 volcanic forcing dataset as was done for CMIP6 (standard runs, hereafter). We used one of the most recent volcanic sulfur inventories as input for the interactive runs, in which aerosol properties are consistent with the model chemistry, microphysics and atmospheric components. We analyzed how the stratospheric aerosol optical depth, the radiative forcing and the climate response to volcanic eruptions differed between interactive and standard runs, and how these compare to observations and proxy records. In particular, we investigate in detail the differences in the response to the large-magnitude Krakatoa 1883 eruption between the two sets of runs. We also discuss differences for the 1979-2015 period where the forcing data in standard runs is directly constrained from satellite observations. Our results shed new light on uncertainties affecting the reconstruction of past volcanic forcing and highlight some of the benefits and disadvantages of using interactive stratospheric aerosol capabilities instead of a unique prescribed volcanic forcing dataset in CMIP’s historical runs.

Methane (CH4), the second most important greenhouse gas in terms of radiative forcing, is on the rise; but there are extensive opportunities for mitigation with existing technologies. Anthropogenic emissions account for around 60% of the global methane source, and the recent atmospheric methane growth rate puts us on a trajectory comparable to the most extreme future methane scenarios in the sixth Coupled Model Intercomparison Project (CMIP6). We use a new methane emissions-driven configuration of the UK Earth System Model (UKESM1) to explore the role of anthropogenic methane in the earth system. The full methane cycle is represented, including surface deposition, chemistry and interactive wetland emissions. As a baseline scenario we used Shared Socioeconomic Pathway 3-7.0 (SSP3-7.0) – the highest methane emissions scenario in CMIP6. In an idealised experiment, all anthropogenic methane emissions were instantaneously stopped from 2015 onwards in a coupled atmosphere-ocean simulation running from 2015-2050, to make a net-zero anthropogenic methane emissions scenario. Within a decade, significant changes can be seen in atmospheric composition and climate, compared to SSP3-7.0. The atmospheric methane burden declines to below pre-industrial levels within 12 years, and by the late 2030s reaches a constant level around 44% below that of the present day (2015). The tropospheric ozone burden and surface mean ozone concentrations decreased by 12% and 15% respectively by 2050 – key in terms of limiting global warming as well as improving air quality and human health. By 2050 the net-zero anthropogenic methane scenario results in a global mean surface temperature (GMST) 1˚C lower than the baseline, a significant value in the context of climate goals such as the Paris Agreement. Through decomposition of the radiation budget, the change in climate can be directly attributed to the reduction in methane and indirectly to the resulting changes in ozone, clouds and ozone precursors such as CO. In addition, the changes in climate result in impacts on the interactive wetland emissions via changes in temperature and wetland extent, highlighting the coupled nature of methane in the earth system. Cessation of anthropogenic methane emissions has profound impacts on near-term warming and on tropospheric ozone, but ultimately cannot single-handedly achieve the necessary reductions for meeting Paris goals.

We present the first results of modelling of the short-living cosmogenic isotope 7Be production, deposition, and transport using the chemistry-climate model SOCOLv3.0 aimed to study solar-terrestrial interactions and climate changes. We implemented an interactive deposition scheme, based on gas tracers with and without nudging to the known meteorological fields. Production of 7Be was modelled using the 3D time-dependent Cosmic Ray induced Atmospheric Cascade (CRAC) model. The simulations were compared with the real concentrations (activity) and depositions measurements of 7Be in the air and water at Finnish stations. We have successfully reproduced and estimated the variability of the cosmogenic isotope 7Be produced by the galactic cosmic rays (GCR) on time scales longer than about a month, for the period of 2002–2008. The agreement between the modelled and measured data is very good (within 12%) providing a solid validation for the ability of the SOCOL CCM to reliably model production, transport, and deposition of cosmogenic isotopes, which is needed for precise studies of cosmic-ray variability in the past.

10 pags., 4 figs., Many Chemistry-Climate Models (CCMs) include a simplified treatment of brominated very short-lived (VSL) species by assuming CHBr as a surrogate for VSL. However, neglecting a comprehensive treatment of VSL in CCMs may yield an unrealistic representation of the associated impacts. Here, we use the Community Atmospheric Model with Chemistry (CAM-Chem) CCM to quantify the tropospheric and stratospheric changes between various VSL chemical approaches with increasing degrees of complexity (i.e., surrogate, explicit, and full). Our CAM-Chem results highlight the improved accuracy achieved by considering a detailed treatment of VSL photochemistry, including sea-salt aerosol dehalogenation and heterogeneous recycling on ice-crystals. Differences between the full and surrogate schemes maximize in the lowermost stratosphere and midlatitude free troposphere, resulting in a latitudinally dependent reduction of ∼1–7 DU in total ozone column and a ∼5%–15% decrease of the OH/HO ratio. We encourage all CCMs to include a complete chemical treatment of VSL in the troposphere and stratosphere., This study has been funded by the European Union's Horizon 2020 Re-search and Innovation program (Project ‘ERC-2016-COG 726349 CLIMAHAL’), and supported by the Consejo Superior de Investigaciones Científicas of Spain. Computing resources and support are provided by the Computational and Information System Laboratory (CISL,2017). R. P. Fernandez would like to thank financial support from CONICET, ANPCyT (PICT 2015-0714), UNCuyo (SeCTyP M032/3853) and UTN (PID 4920-194/2018). NCAR is sponsored by NSF under grant number 1852977. R. J. Salawitch appreciates support from the NASA (grant ACMP 80NSSC19K0983). The authors thank two anonymous reviewers for their constructive comments

The atmospheric chemistry of halogenated species (Cl, Br, I) participates in the global chemical sink of tropospheric ozone and perturbs the oxidising capacity of the troposphere, notably by influencing the atmospheric lifetime of methane. Global chemistry–climate models are commonly used to assess the global budget of ozone and its sensitivity to emissions of its precursors, as well as to project its long-term evolution. Here, we report on the implementation of tropospheric sources and chemistry of halogens in the chemistry–climate model LMDZ-INCA (Laboratoire de Météorologie Dynamique general circulation model, LMDZ, and Interactions with Chemistry and Aerosols, INCA, version with Non-Methane HydroCarbon chemistry, vNMHC) and evaluate halogen effects on the tropospheric ozone budget. Overall, the results show that the model simulates satisfactorily the impact of halogens on the photo-oxidising system in the troposphere, in particular in the marine boundary layer. To quantify the effects of halogen chemistry in LMDZ-INCA, standard metrics representative of the behaviour of the tropospheric chemical system (Ox, HOx, NOx, CH4 and non-methane volatile organic compounds – NMVOCs) are computed with and without halogens. The addition of tropospheric halogens in the LMDZ-INCA model leads to a decrease of 22 % in the ozone burden, 8 % in OH and 33 % in NOx. Sensitivity simulations show for the first time that the inclusion of halogen chemistry makes ozone more sensitive to perturbations in CH4, NOx and NMVOCs. Consistent with other global model studies, the sensitivity of the tropospheric ozone burden to changes from pre-industrial to present-day emissions is found to be ∼20 % lower when tropospheric halogens are taken into account.

We describe the baseline model configuration and simulation characteristics of the Geophysical Fluid Dynamics Laboratory (GFDL)'s Atmosphere Model version 4.1 (AM4.1), which builds on developments at GFDL over 2013–2018 for coupled carbon‐chemistry‐climate simulation as part of the sixth phase of the Coupled Model Intercomparison Project. In contrast with GFDL's AM4.0 development effort, which focused on physical and aerosol interactions and which is used as the atmospheric component of CM4.0, AM4.1 focuses on comprehensiveness of Earth system interactions. Key features of this model include doubled horizontal resolution of the atmosphere (~200 to ~100 km) with revised dynamics and physics from GFDL's previous‐generation AM3 atmospheric chemistry‐climate model. AM4.1 features improved representation of atmospheric chemical composition, including aerosol and aerosol precursor emissions, key land‐atmosphere interactions, comprehensive land‐atmosphere‐ocean cycling of dust and iron, and interactive ocean‐atmosphere cycling of reactive nitrogen. AM4.1 provides vast improvements in fidelity over AM3, captures most of AM4.0's baseline simulations characteristics, and notably improves on AM4.0 in the representation of aerosols over the Southern Ocean, India, and China—even with its interactive chemistry representation—and in its manifestation of sudden stratospheric warmings in the coldest months. Distributions of reactive nitrogen and sulfur species, carbon monoxide, and ozone are all substantially improved over AM3. Fidelity concerns include degradation of upper atmosphere equatorial winds and of aerosols in some regions.

We report airborne measurements of acetaldehyde (CH(3)CHO) during the first and second deployments of the National Aeronautics and Space Administration (NASA) Atmospheric Tomography Mission (ATom). The budget of CH(3)CHO is examined using the Community Atmospheric Model with chemistry (CAM-chem), with a newly-developed online air-sea exchange module. The upper limit of the global ocean net emission of CH(3)CHO is estimated to be 34 Tg a(−1) (42 Tg a(−1) if considering bubble-mediated transfer), and the ocean impacts on tropospheric CH(3)CHO are mostly confined to the marine boundary layer. Our analysis suggests that there is an unaccounted CH(3)CHO source in the remote troposphere and that organic aerosols can only provide a fraction of this missing source. We propose that peroxyacetic acid (PAA) is an ideal indicator of the rapid CH(3)CHO production in the remote troposphere. The higher-than-expected CH(3)CHO measurements represent a missing sink of hydroxyl radicals (and halogen radical) in current chemistry-climate models.

One of the key questions in the air quality and climate sciences is how tropospheric ozone concentrations will change in the future. This will depend on two factors: changes in stratosphere-to-troposphere transport (STT) and changes in tropospheric chemistry. Here we aim to identify robust changes in STT using simulations from the Chemistry Climate Model Initiative (CCMI) under a common climate change scenario (RCP6.0). We use two idealized stratospheric tracers to isolate changes in transport: stratospheric ozone (O3S), which is exactly like ozone but has no chemical sources in the troposphere, and st80, a passive tracer with fixed volume mixing ratio in the stratosphere. We find a robust increase in the tropospheric columns of these two tracers across the models. In particular, stratospheric ozone in the troposphere is projected to increase 10 %–16 % by the end of the 21st century in the RCP6.0 scenario. Future STT is enhanced in the subtropics due to the strengthening of the shallow branch of the Brewer–Dobson circulation (BDC) in the lower stratosphere and of the upper part of the Hadley cell in the upper troposphere. The acceleration of the deep branch of the BDC in the Northern Hemisphere (NH) and changes in eddy transport contribute to increased STT at high latitudes. These STT trends are caused by greenhouse gas (GHG) increases, while phasing out of ozone-depleting substances (ODS) does not lead to robust transport changes. Nevertheless, the decline of ODS increases the reservoir of ozone in the lower stratosphere, which results in enhanced STT of O3S at middle and high latitudes. A higher emission scenario (RCP8.5) produces stronger STT trends, with increases in tropospheric column O3S more than 3 times larger than those in the RCP6.0 scenario by the end of the 21st century.

With low concentrations of tropospheric aerosol, the Southern Ocean offers a "natural laboratory" for studies of aerosol–cloud interactions. Aerosols over the Southern Ocean are produced from biogenic activity in the ocean, which generates sulfate aerosol via dimethylsulfide (DMS) oxidation, and from strong winds and waves that lead to bubble bursting and sea spray emission. Here, we evaluate the representation of Southern Ocean aerosols in the Hadley Centre Global Environmental Model version 3, Global Atmosphere 7.1 (HadGEM3-GA7.1) chemistry–climate model. Compared with aerosol optical depth (AOD) observations from two satellite instruments (the Moderate Resolution Imaging Spectroradiometer, MODIS-Aqua c6.1, and the Multi-angle Imaging Spectroradiometer, MISR), the model simulates too-high AOD during winter and too-low AOD during summer. By switching off DMS emission in the model, we show that sea spray aerosol is the dominant contributor to AOD during winter. In turn, the simulated sea spray aerosol flux depends on near-surface wind speed. By examining MODIS AOD as a function of wind speed from the ERA-Interim reanalysis and comparing it with the model, we show that the sea spray aerosol source function in HadGEM3-GA7.1 overestimates the wind speed dependency. We test a recently developed sea spray aerosol source function derived from measurements made on a Southern Ocean research voyage in 2018. In this source function, the wind speed dependency of the sea spray aerosol flux is less than in the formulation currently implemented in HadGEM3-GA7.1. The new source function leads to good agreement between simulated and observed wintertime AODs over the Southern Ocean; however, it reveals partially compensating errors in DMS-derived AOD. While previous work has tested assumptions regarding the seawater climatology or sea–air flux of DMS, we test the sensitivity of simulated AOD, cloud condensation nuclei and cloud droplet number concentration to three atmospheric sulfate chemistry schemes. The first scheme adds DMS oxidation by halogens and the other two test a recently developed sulfate chemistry scheme for the marine troposphere; one tests gas-phase chemistry only, while the second adds extra aqueous-phase sulfate reactions. We show how simulated sulfur dioxide and sulfuric acid profiles over the Southern Ocean change as a result and how the number concentration and particle size of the soluble Aitken, accumulation and coarse aerosol modes are affected. The new DMS chemistry scheme leads to a 20% increase in the number concentration of cloud condensation nuclei and cloud droplets, which improves agreement with observations. Our results highlight the importance of atmospheric chemistry for simulating aerosols and clouds accurately over the Southern Ocean.

An important source of polar stratospheric clouds (PSCs), which play a crucial role in controlling polar stratospheric ozone depletion, is from the temperature fluctuations induced by mountain waves, enabling stratospheric temperatures to fall below the threshold value for PSC formation in the cold phases of these waves even if the synoptic-scale temperatures are too high. However, this formation mechanism is usually missing in chemistry–climate models because these temperature fluctuations are neither resolved nor parameterised. Here, we investigate the representation of parameterised stratospheric mountain-wave-induced temperature fluctuations over the Antarctic Peninsula from a 30-year run of the global chemistry-climate configuration of the UM-UKCA model against climatologies of Atmospheric Infrared Sounder (AIRS) radiance measurements and high-resolution radiosonde temperature soundings from Rothera. The results demonstrate that the local mountain wave-induced cooling phases computed by the scheme are in relatively good agreement with both sets of observations. For example, the scheme is able to capture the observed probability distribution of the temperature fluctuations, particularly the cold tails of the distribution that are critical for exceeding the temperature threshold for PSC formation. Further analysis shows that the increased stratospheric cooling induced by the scheme results in a large increase in total PSC ‘pseudo-volume’ of the area over the Antarctic Peninsula where the model temperature exceeds the temperature threshold of formation of PSCs.

The organic aerosols (OA) contribute significantly to fine particulate mass in the atmosphere, however, most global climate models do not include elaborate treatment associated with the production of secondary organic aerosols (SOA) involving complex chemical processes to save computational time. As a result, the concentrations of SOA simulated by these climate models are often highly uncertain. Moreover, very limited research has been done on SOA and its precursors, particularly on the contribution of individual sources towards the SOA concentrations across India. In this study, we investigate the sensitivity of the production of SOA from different VOC sources and different atmospheric oxidants by the Community Atmospheric Model version 4 coupled with an extensive interactive atmospheric chemistry module (CAM4-Chem). The main objective of our present research is to understand the contribution of individual sources of VOCs towards the production and distribution of SOA across the Indian region. We carried out a series of systematically designed simulations using the CAM4-Chem model to understand the sensitivity of simulated SOA over the Indian region to changes in only emissions of VOCs from anthropogenic, biogenic, and biomass burning emissions from preindustrial (PI) to present-day (PD) period. In order to avoid the influence of changes in meteorology from PI to PD on the production of SOA, all simulations are performed for the same period from 2004 to 2014 with identical meteorology prescribed to the model based on MERRA2 data, while the VOC emissions from anthropogenic, biogenic, and biomass burning sources are allowed to change from PI to PD in different simulations. Our results show that the simulated distribution of SOA over the Indian region in PD is linked to the significant changes in the emissions of VOCs from anthropogenic, biogenic, and biomass burning emissions sources from PI to PD. We find that the changes in emissions of VOCs from biogenic sources from PI to PD associated with land use and land cover changes contribute significantly along with the changes in emissions from anthropogenic sources towards the total changes in SOA distribution over the Indian region over the same period. The global annual mean burden of SOA from our sensitivity simulations vary in the range of 0.65Tg to 0.80Tg due to variations in emission of different VOCs that are precursors to the production of SOA in the atmosphere. These sensitivity simulations improve our understanding of atmospheric chemistry and specifically about the formation of SOA from different precursor gases originating from diverse anthropogenic, biogenic, and biomass burning emissions sources. More results with greater detail will be presented.

Urban areas are recognised as a significant source of greenhouse gas emissions (GHG), such as carbon dioxide (CO2) and methane (CH4). The total amount of urban GHG emissions, especially for CH4, however, is not well quantified. Here we report on airborne in situ measurements using a Picarro G1301-m analyser aboard the DLR Cessna Grand Caravan to study GHG emissions downwind of the German capital city Berlin. In total, five aircraft-based mass balance experiments were conducted in July 2018 within the Urban Climate Under Change [UC]2 project. The detection and isolation of the Berlin plume was often challenging because of comparatively small GHG signals above variable atmospheric background concentrations. However, on July 20th enhancements of up to 4 ppm CO2 and 21 ppb CH4 were observed over a horizontal extent of roughly 45 to 65 km downwind of Berlin. These enhanced mixing ratios are clearly distinguishable from the background and can partly be assigned to city emissions. The estimated CO2 emission flux of 1.39 ± 0.75 t s-1 is in agreement with current inventories, while the CH4 emission flux of 5.20 ± 1.61 kg s-1 is almost two times larger than the highest reported value in the inventories. We localized the source area with HYSPLIT trajectory calculations and the high resolution numerical model MECO(n) (down to ~1 km), and investigated the contribution from sewage-treatment plants and waste deposition to CH4, which are treated differently by the emission inventories. Our work highlights the importance of a) strong CH4 sources in the surroundings of Berlin and b) a detailed knowledge of GHG inflow mixing ratios to suitably estimate emission rates.

Ozone loss pathways and their rates in the ozone quasi-biennial oscillation (QBO), which is simulated by a chemistry-climate model developed by the Meteorological Research Institute of Japan, are evaluated using an ob- jective pathway analysis program (PAP). The analyzed chemical system contains catalytic cycles caused by NOx , HOx , ClOx , Ox , and BrOx . PAP quantified the rates of all significant catalytic ozone loss cycles, and evaluated the partitioning among these cycles. The QBO amplitude of the sum of all cycles amounts to about 4 and 14 % of the annual mean of the total ozone loss rate at 10 and 20 hPa, respectively. The contribution of catalytic cycles to the QBO of the ozone loss rate is found to be as follows: NOx cycles contribute the largest fraction (50 – 85 %) of the QBO amplitude of the total ozone loss rate; HOx cycles are the second-largest (20 – 30 %) below 30 hPa and the third-largest (about 10 %) above 20 hPa; Ox cycles rank third (5 – 20 %) below 30 hPa and second (about 20 %) above 20 hPa; ClOx cycles rank fourth (5 – 10 %); and BrOx cycles are almost negligible. The relative contribution of the NOx and Ox cycles to the QBO amplitude of ozone loss differs by up to 10 % and 20 %, respectively, from their contribution to the annual mean ozone loss rate. The ozone QBO at 20 hPa is mainly driven by ozone transport, which then alters the ozone loss rate. In contrast, the ozone QBO at 10 hPa is driven chemically by NOx and the temperature dependence of [O]/[O3], which results from the temperature dependence of the reaction O + O2 + M → O3 + M. In addition, the ozone QBO at 10 hPa is influenced by the overhead ozone column, which affects [O]/[O3] (through ozone photolysis) and the ozone production rate (through oxygen photolysis).

Climate impact of aviation is expected to increase further. Aircraft routings are an important measure for climate impact reductions. To find an effective aircraft routing strategy for reducing the impact, the first version of the submodel AirTraf has been developed; this submodel can simulate global air traffic in the ECHAM/MESSy Atmospheric Chemistry (EMAC) model. This paper describes the updated submodel AirTraf 2.0. Seven new aircraft routing options are introduced, including contrail avoidance, minimum economic costs, and minimum climate impact. Example simulations of the new routing options are presented by using around 100 north-Atlantic flights of an Airbus A330 aircraft for a typical winter day. The results clearly show that the family of optimum flight trajectories (three-dimensional) varies according to the routing options. The comparison of the results for various routing options reveals characteristics of the routing with respect to air traffic performances. The minimum cost option obtains a trade-off solution between the minimum time and the minimum fuel solutions. The aircraft routings for contrail avoidance and minimum climate impact reduce the potential climate impact, which is estimated by using algorithmic Climate Change Functions, whereas these two routings increase flight operating costs. A trade-off between the aircraft operating costs and the climate impact is confirmed. The simulation results are compared with literature data and the consistency of the submodel AirTraf 2.0 is verified.

Specified dynamics (SD) schemes relax the circulation in climate models toward a reference meteorology to simulate historical variability. These simulations are widely used to isolate the dynamical contributions to variability and trends in trace gas species. However, it is not clear if trends in the stratospheric overturning circulation are properly reproduced by SD schemes. This study assesses numerous SD schemes and modeling choices in the Community Earth System Model (CESM) Whole Atmosphere Chemistry Climate Model (WACCM) to determine a set of best practices for reproducing interannual variability and trends in tropical stratospheric upwelling estimated by reanalyses. Nudging toward the reanalysis meteorology as is typically done in SD simulations expectedly changes the model’s mean upwelling compared to its free-running state, but does not accurately reproduce upwelling trends present in the underlying reanalysis. In contrast, nudging to anomalies from the climatological winds or from the zonal mean winds and temperatures preserves WACCM’s climatology and better reproduces trends in stratospheric upwelling. An SD scheme’s performance in simulating the acceleration of the shallow branch of the mean meridional circulation from 1980–2017 hinges on its ability to simulate the downward shift of subtropical lower stratospheric wave momentum forcing. Key to this is not nudging the zonal-mean temperature field. Gravity wave momentum forcing, which drives a substantial fraction of the upwelling in WACCM, cannot be constrained by nudging and presents an upper-limit on the performance of these schemes.

31 pags., 12 figs., 6 tabs. -- Open Access funded by Creative Commons Atribution Licence 4.0. -- jame20925-sup-0001_Supporting_Information.pdf, Current chemistry climate models do not include polar emissions and chemistry of halogens. This work presents the first implementation of an interactive polar module into the very short-lived (VSL) halogen version of the Community Atmosphere Model with Chemistry (CAM-Chem) model. The polar module includes photochemical release of molecular bromine, chlorine, and interhalogens from the sea-ice surface, and brine diffusion of iodine biologically produced underneath and within porous sea-ice. It also includes heterogeneous recycling of inorganic halogen reservoirs deposited over fresh sea-ice surfaces and snow-covered regions. The polar emission of chlorine, bromine, and iodine reach approximately 32, 250, and 39 Gg/year for Antarctica and 33, 271, and 4 Gg/year for the Arctic, respectively, with a marked seasonal cycle mainly driven by sunlight and sea-ice coverage. Model results are validated against polar boundary layer measurements of ClO, BrO, and IO, and satellite BrO and IO columns. This validation includes satellite observations of IO over inner Antarctica for which an iodine “leapfrog” mechanism is proposed to transport active iodine from coastal source regions to the interior of the continent. The modeled chlorine and bromine polar sources represent up to 45% and 80% of the global biogenic VSL and VSL emissions, respectively, while the Antarctic sea-ice iodine flux is ~10 times larger than that from the Southern Ocean. We present the first estimate of the contribution of polar halogen emissions to the global tropospheric halogen budget. CAM-Chem includes now a complete representation of halogen sources and chemistry from pole-to-pole and from the Earth's surface up to the stratopause., This study has been funded by the European Research Council Executive Agency under the European Union′s Horizon 2020 Research and Innovation program (Project “ERC‐2016‐COG 726349 CLIMAHAL”) and supported by the Consejo Superior de Investigaciones Científicas (CSIC) of Spain. Computing resources, support, and data storage are provided and maintained by the Computational and Information System Laboratory from the National Center of Atmospheric Research (CISL,2017). R. P. F. would like to thank CONICET, ANPCyT (PICT 2015‐0714), UNCuyo (SeCTyP M032/3853), and UTN (PID 4920‐194/2018) for the financial support. Partial funding for this work was provided by the Korea Polar Research Institute (KOPRI) project (PE18200). The contributions of the University of Bremen have been supported by the State of Bremen, the German Research Foundation (DFG), the German Aerospace (DLR), and the European Space Agency (ESA). We gratefully acknowledge the funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) —Projektnummer 268020496—TRR 172, within the Transregional Collaborative Research Center “ArctiC Amplification: Climate Relevant Atmospheric and SurfaCe Processes,and Feedback Mechanisms (AC)3 ” in subproject C03 as well as the support by the University of Bremen Institutional Strategy Measure M8 in the framework of the DFG Excellence Initiative.

We use the Goddard Earth Observing System Chemistry-Climate Model (GEOSCCM), a contributor to both the 2010 and 2014 WMO Ozone Assessment Reports, to show that inclusion of 5 parts per trillion (ppt) of stratospheric bromine (Bry) from very short-lived substances (VSLS) is responsible for about a decade delay in ozone hole recovery. These results partially explain the significantly later recovery of Antarctic ozone noted in the 2014 report, as bromine from VSLS was not included in the 2010 Assessment. We show multiple lines of evidence that simulations that account for VSLS Bry are in better agreement with both total column BrO and the seasonal evolution of Antarctic ozone reported by the Ozone Monitoring Instrument (OMI) on NASA's Aura satellite. In addition, the near zero ozone levels observed in the deep Antarctic lower stratospheric polar vortex are only reproduced in a simulation that includes this Bry source from VSLS.

To achieve better agreement of simulated Arctic winter stratospheric dynamic with observations assimilation procedure nudging was incorporated in CCM SOCOL. Trajectories based on SOCOL output winds demonstrate the reasonable agreement with trajectories based on reanalysis data inside the polar vortex and can be used for analysis and forecast of ozone related processes in winter-spring seasons. Obtained results of several recent major Arctic SSW events analysis show that CCM SOCOL could be used for SSW forecast over the period up to 8 days.

A new module, ORACLE 2-D, simulating organic aerosol formation and evolution in the atmosphere has been developed and evaluated. The module calculates the concentrations of surrogate organic species in two-dimensional space defined by volatility and oxygen-to-carbon ratio. It is implemented into the EMAC global chemistry–climate model, and a comprehensive evaluation of its performance is conducted using an aerosol mass spectrometer (AMS) factor analysis dataset derived from almost all major field campaigns that took place globally during the period 2001–2010. ORACLE 2-D uses a simple photochemical aging scheme that efficiently simulates the net effects of fragmentation and functionalization of the organic compounds. The module predicts not only the mass concentration of organic aerosol (OA) components, but also their oxidation state (in terms of O : C), which allows for their classification into primary OA (POA, chemically unprocessed), fresh secondary OA (SOA, low oxygen content), and aged SOA (highly oxygenated). The explicit simulation of chemical OA conversion from freshly emitted compounds to a highly oxygenated state during photochemical aging enables the tracking of hygroscopicity changes in OA that result from these reactions. ORACLE 2-D can thus compute the ability of OA particles to act as cloud condensation nuclei and serves as a tool to quantify the climatic impact of OA.

Primary biological aerosol particles (PBAPs) may impact human health and aerosol-climate interactions. The role of PBAPs in the earth system is associated with large uncertainties, related to source estimates and atmospheric transport. We used a chemistry-climate model to simulate PBAPs in the atmosphere including bacteria, fungal spores and pollen. Three fungal spore emission parameterizations have been evaluated against an updated set of spore counts synthesized from observations reported in the literature. The comparison indicates an optimal fit for the emission parameterization proposed by Heald and Spracklen (2009), although the model significantly over-predicts PBAP concentrations in some locations. Additional evaluation was performed by comparing our combined bacteria and fungal spore simulations to a global dataset of fluorescent biological aerosol particle (FBAP) concentrations. The model predicts the sum total of measured PBAP concentrations relatively well, with an over- or under-prediction of less than a factor of 2 compared to FBAP. The ratio of bacteria to fungal spores reflects a greater difference, however, and the simulated bacteria concentrations outnumber the simulated fungal spore concentrations in almost all locations. Further, the modeled fungal spore results under-predict the FBAP concentrations, which are used here as a rough proxy for spores. Uncertainties related to technical aspects of the FBAP and direct-counting spore measurements challenge the ability to further refine quantitative comparison on this scale. We estimate that the global PBAPs mass concentration (apart from desert dust and sea salt aerosols), i.e. of fungal spores and pollen, amounts to 19 % and 52 % of the total aerosol mass, respectively.

Bromine from very short lived substances (VSLS), primarily from natural oceanic sources, contributes substantially to the stratospheric bromine loading. This source of stratospheric bromine has so far been ignored in most chemistry climate model calculations of stratospheric ozone trends. Here we present a transient simulation with the chemistry climate model EMAC for the period 1960–2005 including emissions of the five brominated VSLS CHBr3, CH2Br2, CH2BrCl, CHBrCl2, and CHBr2Cl. The emissions lead to a realistic stratospheric bromine loading of about 20 pptv for present-day conditions. Comparison with a standard model simulation without VSLS shows large differences in modeled ozone in the extratropical lowermost stratosphere and in the troposphere. Differences in ozone maximize in the Antarctic Ozone Hole, resulting in more than 20% less ozone when VSLS are included. Even though the emissions of VSLS are assumed to be constant in time, the model simulation with VSLS included shows a much larger ozone decrease in the lowermost stratosphere during the 1979–1995 period and a faster ozone increase during 1996–2005, in better agreement with observed ozone trends than the standard simulation without VSLS emissions.

Within the framework of the global chemistry climate model ECHAM-HAMMOZ a novel explicit coupling between the sectional aerosol model HAM-SALSA and the chemistry model MOZ was established to form isoprene derived secondary organic aerosol (iSOA). Isoprene oxidation in the chemistry model MOZ is described by a semi-explicit scheme consisting of 147 reactions, embedded in a detailed atmospheric chemical mechanism with a total of 779 reactions. Low volatile compounds (LVOC) produced during isoprene photooxidation are identified and explicitly partitioned by HAM-SALSA. A group contribution method was used to estimate their evaporation enthalpies and corresponding saturation vapor pressures, which are used by HAM-SALSA to calculate the saturation concentration of each LVOC. With this method, every single precursor is tracked in terms of condensation and evaporation in each aerosol size bin. This approach lead to the identification of ISOP(OOH)2 as a main contributor to iSOA formation. Further, reactive uptake of isoprene epoxidiols (IEPOX) and isoprene derived glyoxal were included as iSOA sources. The parameterization of IEPOX reactive uptake includes a dependency on aerosol pH value. This model framework connecting semi-explicit isoprene oxidation with explicit treatment of aerosol tracers leads to a global, annual isoprene SOA yield of 16 % relative to the primary oxidation of isoprene by OH, NO3, and ozone. With 445 Tg (392 TgC) isoprene emitted, an iSOA source of 148 Tg (61 TgC) is simulated. The major part of iSOA in ECHAM-HAMMOZ is produced by IEPOX (24.4 TgC) and ISOP(OOH)2 (28.3 TgC). The main sink process is particle wet deposition which removes 143 Tg (59 TgC). The iSOA burden reaches 1.6 Tg (0.7 TgC) in the year 2012.

Recent research has demonstrated that the concentrations of anthropogenic halocarbons have decreased in response to the worldwide phaseout of ozone depleting substances. Yet in 2015 the Antarctic ozone hole reached a historical record daily average size in October. Model simulations with specified dynamics and temperatures based on a reanalysis suggested that the record size was likely due to the eruption of Calbuco but did not allow for fully coupled dynamical or thermal feedbacks. We present simulations of the impact of the 2015 Calbuco eruption on the stratosphere using the Whole Atmosphere Community Climate Model with interactive dynamics and temperatures. Comparisons of the interactive and specified dynamics simulations indicate that chemical ozone depletion due to volcanic aerosols played a key role in establishing the record-sized ozone hole of October 2015. The analysis of an ensemble of interactive simulations with and without volcanic aerosols suggests that the forced response to the eruption of Calbuco was an increase in the size of the ozone hole by 4.5 × 10⁶ km².

We have improved an ozone DIfferential Absorption Lidar (DIAL) system, originally developed in March 2010. The improved DIAL system consists of a Nd:YAG laser and a 2 m Raman cell filled with 8.1 × 105 Pa of CO2 gas which generate four Stokes lines (276, 287, 299, and 312 nm) of stimulated Raman scattering, and two receiving telescopes with diameters of 49 and 10 cm. Using this system, 44 ozone profiles were observed in the 1–6 km altitude range over Saga (33.24° N, 130.29° E) in 2012. High-ozone layers were observed at around 2 km altitude during April and May. Ozone column amounts within the 1–6 km altitude range were almost constant (19.1 DU on average) from January to March, and increased to 26.7 DU from late April to July. From mid-July through August, ozone column amounts decreased greatly to 14.3 DU because of exchanges of continental and maritime air masses. Then in mid-September they increased again to 22.1 DU within 1−6 km, and subsequently decreased slowly to 17.3 DU, becoming almost constant by December. The Meteorological Research Institute's chemistry–climate model version 2 (MRI-CCM2) successfully predicted most of these ozone variations with the following exceptions. MRI-CCM2 could not predict the high-ozone volume mixing ratios measured at around 2 km altitude on 5 May and 11 May, possibly in part because emissions were assumed in the model to be constant (climatological data were used). Ozone volume mixing ratios predicted by MRI-CCM2 were low in the 2–6 km range on 7 July and high in the 1–4 km range on 19 July compared with those measured by DIAL.

Recent observations reveal a seasonally occurring layer of aerosol located from 0° to 100°E, 20° to 45°N and extending vertically from about 13 km to 18 km; this has been termed the Asian tropopause aerosol layer (ATAL), and its existence is closely associated with the Asian summer monsoon circulation. Observational studies argue that the ATAL is a recent phenomenon, as the layer is not observed in the satellite record prior to 1998. This suggests that the ATAL may be of anthropogenic origin associated with a shift in the dominant regional emission of sulfur dioxide (SO2) to China and India in the late 1990s. Here we test the hypothesis that SO2emitted from Asia led to the formation of the ATAL using an aerosol microphysical model coupled to a global chemistry climate model. This is the first modeling study to specifically examine the ATAL and its possible origin. From our results, we conclude that the ATAL is most likely due to anthropogenic emissions, but its source cannot solely be attributed to emissions from Asia. Specifically, the results indicate that Chinese and Indian emissions contribute ∼30% of the sulfate aerosol extinction in the ATAL during volcanically quiescent periods. We also show that even small volcanic eruptions preclude our ability to make any conclusions about the existence of the ATAL before 1998 with observations alone.

Tropospheric ozone is a potent greenhouse gas, biological irritant and significant source of highly reactive hydroxyl radicals. Simulations with a chemistry climate model suggest that shifts in atmospheric circulation can account for the seasonally dependent trends in tropospheric ozone levels observed at Mauna Loa, Hawaii, over the past three decades.

We present an overview of state-of-the-art chemistry–climate and chemistry transport models that are used within phase 1 of the Chemistry–Climate Model Initiative (CCMI-1). The CCMI aims to conduct a detailed evaluation of participating models using process-oriented diagnostics derived from observations in order to gain confidence in the models' projections of the stratospheric ozone layer, tropospheric composition, air quality, where applicable global climate change, and the interactions between them. Interpretation of these diagnostics requires detailed knowledge of the radiative, chemical, dynamical, and physical processes incorporated in the models. Also an understanding of the degree to which CCMI-1 recommendations for simulations have been followed is necessary to understand model responses to anthropogenic and natural forcing and also to explain inter-model differences. This becomes even more important given the ongoing development and the ever-growing complexity of these models. This paper also provides an overview of the available CCMI-1 simulations with the aim of informing CCMI data users.

The fractional release factor (FRF) gives information on the amount of a halocarbon that is released at some point into the stratosphere from its source form to the inorganic form, which can harm the ozone layer through catalytic reactions. The quantity is of major importance because it directly affects the calculation of the ozone depletion potential (ODP). In this context time-independent values are needed which, in particular, should be independent of the trends in the tropospheric mixing ratios (tropospheric trends) of the respective halogenated trace gases. For a given atmospheric situation, such FRF values would represent a molecular property.We analysed the temporal evolution of FRF from ECHAM/MESSy Atmospheric Chemistry (EMAC) model simulations for several halocarbons and nitrous oxide between 1965 and 2011 on different mean age levels and found that the widely used formulation of FRF yields highly time-dependent values. We show that this is caused by the way that the tropospheric trend is handled in the widely used calculation method of FRF.Taking into account chemical loss in the calculation of stratospheric mixing ratios reduces the time dependence in FRFs. Therefore we implemented a loss term in the formulation of the FRF and applied the parameterization of a mean arrival time to our data set.We find that the time dependence in the FRF can almost be compensated for by applying a new trend correction in the calculation of the FRF. We suggest that this new method should be used to calculate time-independent FRFs, which can then be used e.g. for the calculation of ODP.

We present the third generation of the coupled chemistry–climate model (CCM) SOCOL (modeling tools for studies of SOlar Climate Ozone Links). The most notable modifications compared to the previous model version are (1) the dynamical core has been updated with the fifth generation of the middle-atmosphere general circulation model MA-ECHAM (European Centre/HAMburg climate model), and (2) the advection of the chemical species is now calculated by a mass-conserving and shape-preserving flux-form transport scheme instead of the previously used hybrid advection scheme. The whole chemistry code has been rewritten according to the ECHAM5 infrastructure and transferred to Fortran95. In contrast to its predecessors, SOCOLvs3 is now fully parallelized. The performance of the new SOCOL version is evaluated on the basis of transient model simulations (1975–2004) with different horizontal (T31 and T42) resolutions, following the approach of the CCMVal-1 model validation activity. The advanced advection scheme significantly reduces the artificial loss and accumulation of tracer mass in regions with strong gradients that was observed in previous model versions. Compared to its predecessors, SOCOLvs3 generally shows more realistic distributions of chemical trace species, especially of total inorganic chlorine, in terms of the mean state, but also of the annual and interannual variability. Advancements with respect to model dynamics are for example a better representation of the stratospheric mean state in spring, especially in the Southern Hemisphere, and a slowdown of the upward propagation in the tropical lower stratosphere. Despite a large number of improvements model deficiencies still remain. Examples include a too-fast vertical ascent and/or horizontal mixing in the tropical stratosphere, the cold temperature bias in the lowermost polar stratosphere, and the overestimation of polar total ozone loss during Antarctic springtime., Geoscientific Model Development, 6 (5), ISSN:1991-9603, ISSN:1991-959X