Explore the vast range of titles available.

Explosive volcanic eruptions affect climate, but how climate change affects the stratospheric volcanic sulfate aerosol lifecycle and radiative forcing remains unexplored. We combine an eruptive column model with an aerosol-climate model to show that the stratospheric aerosol optical depth perturbation from frequent moderate-magnitude tropical eruptions (e.g. Nabro 2011) will be reduced by 75% in a high-end warming scenario compared to today, a consequence of future tropopause height rise and unchanged eruptive column height. In contrast, global-mean radiative forcing, stratospheric warming and surface cooling from infrequent large-magnitude tropical eruptions (e.g. Mt. Pinatubo 1991) will be exacerbated by 30%, 52 and 15% in the future, respectively. These changes are driven by an aerosol size decrease, mainly caused by the acceleration of the Brewer-Dobson circulation, and an increase in eruptive column height. Quantifying changes in both eruptive column dynamics and aerosol lifecycle is therefore key to assessing the climate response to future eruptions. How climate change influences the lifecycle of stratospheric volcanic aerosols and the associated radiative forcing is unknown. Here, the authors present model experiments suggesting that climate change amplifies the forcing of large-magnitude tropical eruptions but reduces the forcing of moderate-magnitude tropical eruptions.

This paper quantifies the pre-industrial (1850) to present-day (2014) effective radiative forcing (ERF) of anthropogenic emissions of NOX, volatile organic compounds (VOCs; including CO), SO2, NH3, black carbon, organic carbon, and concentrations of methane, N2O and ozone-depleting halocarbons, using CMIP6 models. Concentration and emission changes of reactive species can cause multiple changes in the composition of radiatively active species: tropospheric ozone, stratospheric ozone, stratospheric water vapour, secondary inorganic and organic aerosol, and methane. Where possible we break down the ERFs from each emitted species into the contributions from the composition changes. The ERFs are calculated for each of the models that participated in the AerChemMIP experiments as part of the CMIP6 project, where the relevant model output was available. The 1850 to 2014 multi-model mean ERFs (± standard deviations) are −1.03 ± 0.37 W/sq.m for SO2 emissions, −0.25 ± 0.09 W/sq.m for organic carbon (OC), 0.15 ± 0.17 W/sq.m for black carbon (BC) and −0.07 ± 0.01 W/sq.m for NH3. For the combined aerosols (in the piClim-aer experiment) it is −1.01 ± 0.25 W/sq.m. The multi-model means for the reactive well-mixed greenhouse gases (including any effects on ozone and aerosol chemistry) are 0.67 ± 0.17 W/sq.m for methane (CH4), 0.26 ± 0.07 W/sq.m for nitrous oxide (N2O) and 0.12 ± 0.2 W/sq.m for ozone-depleting halocarbons (HC). Emissions of the ozone precursors nitrogen oxides (NOx), volatile organic compounds and both together (O3) lead to ERFs of 0.14 ± 0.13, 0.09 ± 0.14 and 0.20 ± 0.07 W/sq.m respectively. The differences in ERFs calculated for the different models reflect differences in the complexity of their aerosol and chemistry schemes, especially in the case of methane where tropospheric chemistry captures increased forcing from ozone production.

Biogenic volatile organic compounds (BVOCs) affect climate via changes to aerosols, aerosol-cloud interactions (ACI), ozone and methane. BVOCs exhibit dependence on climate (causing a feedback) and land use but there remains uncertainty in their net climatic impact. One factor is the description of BVOC chemistry. Here, using the earth-system model UKESM1, we quantify chemistry’s influence by comparing the response to doubling BVOC emissions in the pre-industrial with standard and state-of-science chemistry. The net forcing (feedback) is positive: ozone and methane increases and ACI changes outweigh enhanced aerosol scattering. Contrary to prior studies, the ACI response is driven by cloud droplet number concentration (CDNC) reductions from suppression of gas-phase SO 2 oxidation. With state-of-science chemistry the feedback is 43% smaller as lower oxidant depletion yields smaller methane increases and CDNC decreases. This illustrates chemistry’s significant influence on BVOC’s climatic impact and the more complex pathways by which BVOCs influence climate than currently recognised. The modelling of BVOC chemistry strongly affects how doubling of BVOC emissions affects climate. Lower oxidant depletion with state-of-science chemistry leads to 43% smaller positive forcing from smaller methane increases and cloud albedo decreases.

The evolution of volcanic sulfur and the resulting radiative forcing following explosive volcanic eruptions is well understood. Petrological evidence suggests that significant amounts of halogens may be co-emitted alongside sulfur in some explosive volcanic eruptions, and satellite evidence indicates that detectable amounts of these halogens may reach the stratosphere. In this study, we utilise an aerosol–chemistry–climate model to simulate stratospheric volcanic eruption emission scenarios of two sizes, both with and without co-emission of volcanic halogens, in order to understand how co-emitted halogens may alter the life cycle of volcanic sulfur, stratospheric chemistry, and the resulting radiative forcing. We simulate a large (10 Tg of SO2) and very large (56 Tg of SO2) sulfur-only eruption scenario and a corresponding large (10 Tg SO2, 1.5 Tg HCl, 0.0086 Tg HBr) and very large (56 Tg SO2, 15 Tg HCl, 0.086 Tg HBr) co-emission eruption scenario. The eruption scenarios simulated in this work are hypothetical, but they are comparable to Volcanic Explosivity Index (VEI) 6 (e.g. 1991 Mt Pinatubo) and VEI 7 (e.g. 1257 Mt Samalas) eruptions, representing 1-in-50–100-year and 1-in-500–1000-year events, respectively, with plausible amounts of co-emitted halogens based on satellite observations and volcanic plume modelling. We show that co-emission of volcanic halogens and sulfur into the stratosphere increases the volcanic effective radiative forcing (ERF) by 24 % and 30 % in large and very large co-emission scenarios compared to sulfur-only emission. This is caused by an increase in both the forcing from volcanic aerosol–radiation interactions (ERFari) and composition of the stratosphere (ERFclear,clean). Volcanic halogens catalyse the destruction of stratospheric ozone, which results in significant stratospheric cooling, offsetting the aerosol heating simulated in sulfur-only scenarios and resulting in net stratospheric cooling. The ozone-induced stratospheric cooling prevents aerosol self-lofting and keeps the volcanic aerosol lower in the stratosphere with a shorter lifetime. This results in reduced growth by condensation and coagulation and a smaller peak global-mean effective radius compared to sulfur-only simulations. The smaller effective radius found in both co-emission scenarios is closer to the peak scattering efficiency radius of sulfate aerosol, and thus co-emission of halogens results in larger peak global-mean ERFari (6 % and 8 %). Co-emission of volcanic halogens results in significant stratospheric ozone, methane, and water vapour reductions, resulting in significant increases in peak global-mean ERFclear,clean (> 100 %), predominantly due to ozone loss. The dramatic global-mean ozone depletion simulated in large (22 %) and very large (57 %) co-emission scenarios would result in very high levels of UV exposure on the Earth's surface, with important implications for society and the biosphere. This work shows for the first time that co-emission of plausible amounts of volcanic halogens can amplify the volcanic ERF in simulations of explosive eruptions. It highlights the need to include volcanic halogen emissions when simulating the climate impacts of past or future eruptions, as well as the necessity to maintain space-borne observations of stratospheric compounds to better constrain the stratospheric injection estimates of volcanic eruptions.

The South Asian summer monsoon supplies over 80 % of India's precipitation. Industrialization over the past few decades has resulted in severe aerosol pollution in India. Understanding monsoonal sensitivity to aerosol emissions in general circulation models (GCMs) could improve predictability of observed future precipitation changes. The aims here are (1) to assess the role of aerosols in India's monsoon precipitation and (2) to determine the roles of local and regional emissions. For (1), we study the Precipitation Driver Response Model Intercomparison Project experiments. We find that the precipitation response to changes in black carbon is highly uncertain with a large intermodel spread due in part to model differences in simulating changes in cloud vertical profiles. Effects from sulfate are clearer; increased sulfate reduces Indian precipitation, a consistency through all of the models studied here. For (2), we study bespoke simulations, with reduced Chinese and/or Indian emissions in three GCMs. A significant increase in precipitation (up to ∼20 %) is found only when both countries' sulfur emissions are regulated, which has been driven in large part by dynamic shifts in the location of convective regions in India. These changes have the potential to restore a portion of the precipitation losses induced by sulfate forcing over the last few decades.

We provide an overview of the REF-C1SD specified-dynamics experiment that was conducted as part of phase 1 of the Chemistry-Climate Model Initiative (CCMI). The REF-C1SD experiment, which consisted of mainly nudged general circulation models (GCMs) constrained with (re)analysis fields, was designed to examine the influence of the large-scale circulation on past trends in atmospheric composition. The REF-C1SD simulations were produced across various model frameworks and are evaluated in terms of how well they represent different measures of the dynamical and transport circulations. In the troposphere there are large (~40 %) differences in the climatological mean distributions, seasonal cycle amplitude, and trends of the meridional and vertical winds. In the stratosphere there are similarly large (~50 %) differences in the magnitude, trends and seasonal cycle amplitude of the transformed Eulerian mean circulation and among various chemical and idealized tracers. At the same time, interannual variations in nearly all quantities are very well represented, compared to the underlying reanalyses. We show that the differences in magnitude, trends and seasonal cycle are not related to the use of different reanalysis products; rather, we show they are associated with how the simulations were implemented, by which we refer both to how the large-scale flow was prescribed and to biases in the underlying free-running models. In most cases these differences are shown to be as large or even larger than the differences exhibited by free-running simulations produced using the exact same models, which are also shown to be more dynamically consistent. Overall, our results suggest that care must be taken when using specified-dynamics simulations to examine the influence of large-scale dynamics on composition.

Here we present results from an evaluation of model simulations from the International Hemispheric Transport of Air Pollution Phase II (HTAPII) and Chemistry Climate Model Initiative (CCMI) model inter-comparison projects against a comprehensive series of ground-based, aircraft and satellite observations of ozone mixing ratios made at various locations across India. The study focuses on the recent past (observations from 2008 to 2013, models from 2009–2010) as this is most pertinent to understanding the health impacts of ozone. To our understanding this is the most comprehensive evaluation of these models' simulations of ozone across the Indian subcontinent to date. This study highlights some significant successes and challenges that the models face in representing the oxidative chemistry of the region. The multi-model range in area-weighted surface ozone over the Indian subcontinent is 37.26–56.11 ppb, whilst the population-weighted range is 41.38–57.5 ppb. When compared against surface observations from the Modelling Atmospheric Pollution and Networking (MAPAN) network of eight semi-urban monitoring sites spread across India, we find that the models tend to simulate higher ozone than that which is observed. However, observations of NOx and CO tend to be much higher than modelled mixing ratios, suggesting that the underlying emissions used in the models do not characterise these regions accurately and/or that the resolution of the models is not adequate to simulate the photo-chemical environment of these surface observations. Empirical orthogonal function (EOF) analysis is used in order to identify the extent to which the models agree with regards to the spatio-temporal distribution of the tropospheric ozone column, derived using OMI-MLS observations. We show that whilst the models agree with the spatial pattern of the first EOF of observed tropospheric ozone column, most of the models simulate a peak in the first EOF seasonal cycle represented by principle component 1, which is later than the observed peak. This suggests a widespread systematic bias in the timing of emissions or some other unknown seasonal process. In addition to evaluating modelled ozone mixing ratios, we explore modelled emissions of NOx, CO, volatile organic compounds (VOCs) and the ozone response to the emissions. We find a high degree of variation in emissions from non-anthropogenic sources (e.g. lightning NOx and biomass burning CO) between models. Total emissions of NOx and CO over India vary more between different models in the same model intercomparison project (MIP) than the same model used in different MIPs, making it impossible to diagnose whether differences in modelled ozone are due to emissions or model processes. We therefore recommend targeted experiments to pinpoint the exact causes of discrepancies between modelled and observed ozone and ozone precursors for this region. To this end, a higher density of long-term monitoring sites measuring not only ozone but also ozone precursors including speciated VOCs, located in more rural regions of the Indian subcontinent, would enable improvements in assessing the biases in models run at the resolution found in HTAPII and CCMI.

The modeling study presented here aims to estimate how uncertainties in global hydroxyl radical (OH) distributions, variability, and trends may contribute to resolving discrepancies between simulated and observed methane (CH4) changes since 2000. A multi-model ensemble of 14 OH fields was analyzed and aggregated into 64 scenarios to force the offline atmospheric chemistry transport model LMDz (Laboratoire de Meteorologie Dynamique) with a standard CH4 emission scenario over the period 2000–2016. The multi-model simulated global volume-weighted tropospheric mean OH concentration ([OH]) averaged over 2000–2010 ranges between 8.7×105 and 12.8×105 molec cm−3. The inter-model differences in tropospheric OH burden and vertical distributions are mainly determined by the differences in the nitrogen oxide (NO) distributions, while the spatial discrepancies between OH fields are mostly due to differences in natural emissions and volatile organic compound (VOC) chemistry. From 2000 to 2010, most simulated OH fields show an increase of 0.1–0.3×105 molec cm−3 in the tropospheric mean [OH], with year-to-year variations much smaller than during the historical period 1960–2000. Once ingested into the LMDz model, these OH changes translated into a 5 to 15 ppbv reduction in the CH4 mixing ratio in 2010, which represents 7 %–20 % of the model-simulated CH4 increase due to surface emissions. Between 2010 and 2016, the ensemble of simulations showed that OH changes could lead to a CH4 mixing ratio uncertainty of >±30 ppbv. Over the full 2000–2016 time period, using a common state-of-the-art but nonoptimized emission scenario, the impact of [OH] changes tested here can explain up to 54 % of the gap between model simulations and observations. This result emphasizes the importance of better representing OH abundance and variations in CH4 forward simulations and emission optimizations performed by atmospheric inversions.

Various forms of geoengineering have been proposed to counter anthropogenic climate change. Methods which aim to modify the Earth's energy balance by reducing insolation are often subsumed under the term solar radiation management (SRM). Here, we present results of a standard SRM modelling experiment in which the incoming solar irradiance is reduced to offset the global mean warming induced by a quadrupling of atmospheric carbon dioxide. For the first time in an atmosphere–ocean coupled climate model, we include atmospheric composition feedbacks for this experiment. While the SRM scheme considered here could offset greenhouse gas induced global mean surface warming, it leads to important changes in atmospheric composition. We find large stratospheric ozone increases that induce significant reductions in surface UV-B irradiance, which would have implications for vitamin D production. In addition, the higher stratospheric ozone levels lead to decreased ozone photolysis in the troposphere. In combination with lower atmospheric specific humidity under SRM, this results in overall surface ozone concentration increases in the idealized G1 experiment. Both UV-B and surface ozone changes are important for human health. We therefore highlight that both stratospheric and tropospheric ozone changes must be considered in the assessment of any SRM scheme, due to their important roles in regulating UV exposure and air quality.

The amount of dioxygen (O2) in the atmosphere may have varied from as little as 5 % to as much as 35 % during the Phanerozoic eon (54 Ma–present). These changes in the amount of O2 are large enough to have led to changes in atmospheric mass, which may alter the radiative budget of the atmosphere, leading to this mechanism being invoked to explain discrepancies between climate model simulations and proxy reconstructions of past climates. Here, we present the first fully 3-D numerical model simulations to investigate the climate impacts of changes in O2 under different climate states using the coupled atmosphere–ocean Hadley Centre Global Environmental Model version 3 (HadGEM3-AO) and Hadley Centre Coupled Model version 3 (HadCM3-BL) models. We show that simulations with an increase in O2 content result in increased global-mean surface air temperature under conditions of a pre-industrial Holocene climate state, in agreement with idealised 1-D and 2-D modelling studies. We demonstrate the mechanism behind the warming is complex and involves a trade-off between a number of factors. Increasing atmospheric O2 leads to a reduction in incident shortwave radiation at the Earth's surface due to Rayleigh scattering, a cooling effect. However, there is a competing warming effect due to an increase in the pressure broadening of greenhouse gas absorption lines and dynamical feedbacks, which alter the meridional heat transport of the ocean, warming polar regions and cooling tropical regions. Case studies from past climates are investigated using HadCM3-BL and show that, in the warmest climate states in the Maastrichtian (72.1–66.0 Ma), increasing oxygen may lead to a temperature decrease, as the equilibrium climate sensitivity is lower. For the Asselian (298.9–295.0 Ma), increasing oxygen content leads to a warmer global-mean surface temperature and reduced carbon storage on land, suggesting that high oxygen content may have been a contributing factor in preventing a “Snowball Earth” during this period of the early Permian. These climate model simulations reconcile the surface temperature response to oxygen content changes across the hierarchy of model complexity and highlight the broad range of Earth system feedbacks that need to be accounted for when considering the climate response to changes in atmospheric oxygen content.