Search

Accurately quantifying volcanic impacts on climate is a key requirement for robust attribution of anthropogenic climate change. Here we use the Unified Model – United Kingdom Chemistry and Aerosol (UM-UKCA) composition–climate model to simulate the global dispersion of the volcanic aerosol clouds from the three largest eruptions of the 20th century: 1963 Mt Agung, 1982 El Chichón, and 1991 Mt Pinatubo. The model has interactive stratospheric chemistry and aerosol microphysics, with coupled aerosol–radiation interactions for realistic composition–dynamics feedbacks. Our simulations align with the design of the Interactive Stratospheric Aerosol Model Intercomparison (ISA-MIP) "Historical Eruption SO2 Emissions Assessment". For each eruption, we perform three-member ensemble model experiments for upper, mid-point, and lower estimates of SO2 emission, each re-initialised from a control run to approximately match the observed transition in the phase of the quasi-biennial oscillation (QBO) in the 6 months after the eruptions. With this experimental design, we assess how each eruption's emitted SO2 translates into a tropical reservoir of volcanic aerosol and analyse the subsequent dispersion to mid-latitudes. We compare the simulations to the volcanic forcing datasets (e.g. Space-based Stratospheric Aerosol Climatology (GloSSAC); , and) that are used in historical integrations for the two most recent Coupled Model Intercomparison Project (CMIP) assessments. For Pinatubo and El Chichón, we assess the vertical extent of the simulated volcanic clouds by comparing modelled extinction to the Stratospheric Aerosol and Gas Experiment (SAGE-II) v7.0 satellite measurements and to 1964–1965 Northern Hemisphere ground-based lidar measurements for Agung. As an independent test for the simulated volcanic forcing after Pinatubo, we also compare simulated shortwave (SW) and longwave (LW) top-of-the-atmosphere radiative forcings to the flux anomalies measured by the Earth Radiation Budget Experiment (ERBE) satellite instrument. For the Pinatubo simulations, an injection of 10 to 14 Tg SO2 gives the best match to the High Resolution Infrared Sounder (HIRS) satellite-derived global stratospheric sulfur burden, with good agreement also with SAGE-II mid-visible and near-infra-red extinction measurements. This 10–14 Tg range of emission also generates a heating of the tropical stratosphere that is consistent with the temperature anomaly present in the ERA-Interim reanalysis. For El Chichón, the simulations with 5 and 7 Tg SO2 emission give best agreement with the observations. However, these simulations predict a much deeper volcanic cloud than represented in the GloSSAC dataset, which is largely based on an interpolation between Stratospheric Aerosol Measurements (SAM-II) satellite and aircraft measurements. In contrast, these simulations show much better agreement during the SAGE-II period after October 1984. For 1963 Agung, the 9 Tg simulation compares best to the forcing datasets with the model capturing the lidar-observed signature of the altitude of peak extinction descending from 20 km in 1964 to 16 km in 1965. Overall, our results indicate that the downward adjustment to SO2 emission found to be required by several interactive modelling studies when simulating Pinatubo is also needed when simulating the Agung and El Chichón aerosol clouds. This strengthens the hypothesis that interactive stratospheric aerosol models may be missing an important removal or re-distribution process (e.g. effects of co-emitted ash) which changes how the tropical reservoir of volcanic aerosol evolves in the initial months after an eruption. Our model comparisons also identify potentially important inhomogeneities in the CMIP6 dataset for all three eruption periods that are hard to reconcile with variations predicted in the interactive stratospheric aerosol simulations. We also highlight large differences between the CMIP5 and CMIP6 volcanic aerosol datasets for the Agung and El Chichón periods. Future research should aim to reduce this uncertainty by reconciling the datasets with additional stratospheric aerosol observations. [ABSTRACT FROM AUTHOR]

Accurate quantification of the effects of volcanic eruptions on climate is a key requirement for better attribution of anthropogenic climate change. Here we use the UM-UKCA composition-climate model to simulate the atmospheric evolution of the volcanic aerosol clouds from the three largest eruptions of the 20th century: 1963 Agung, 1982 El Chichón and 1991 Pinatubo. The model has interactive stratospheric chemistry and aerosol microphysics, with coupled aerosol-radiation interactions for realistic composition-dynamics feedbacks. Our simulations align with the design of the Interactive Stratospheric Aerosol Model Intercomparison (ISA-MIP) "Historical Eruption SO2 Emissions Assessment". For each eruption, we perform 3-member ensemble model experiments with upper, mid-point and lower estimates for SO2 emission, each initialised to a meteorological state to match the observed phase of the quasi-biennial oscillation (QBO) at the times of the eruptions. We assess how each eruption's emitted SO2 evolves into a tropical reservoir of volcanic aerosol and analyse the subsequent dispersion to mid-latitudes. We compare the simulations to the three volcanic forcing datasets used in historical integrations for the two most recent Coupled Model Intercomparison Project (CMIP) assessments: the Global Space-based Stratospheric Aerosol Climatology (GloSSAC) for CMIP6, and the Sato et al. (1993) and Ammann et al. (2003) datasets used in CMIP5. We also assess the vertical extent of the volcanic aerosol clouds by comparing simulated extinction to Stratospheric Aerosol and Gas Experiment II (SAGE-II) v7.0 satellite aerosol data (1985-1995) for Pinatubo and El Chichón, and to 1964-65 northern hemisphere ground-based lidar measurements for Agung. As an independent test for the simulated volcanic forcing after Pinatubo, we also compare to the shortwave (SW) and longwave (LW) Top-of-the-Atmosphere flux anomalies measured by the Earth Radiation Budget Experiment (ERBE) satellite instrument. For the Pinatubo simulations, an injection of 10 to 14 Tg SO2 gives the best match to the High Resolution Infrared Sounder (HIRS) satellite-derived global stratospheric sulphur burden, with good agreement also to SAGE-II mid-visible and near-infrared extinction measurements. This 10-14 Tg range of emission also generates a heating of the tropical stratosphere that is comparable with the temperature anomaly seen in the ERA-Interim reanalyses. For El Chichón the simulations with 5 Tg and 7 Tg SO2 emission give best agreement with the observations. However, these runs predict a much deeper volcanic cloud than present in the CMIP6 data, with much higher aerosol extinction than the GloSSAC data up to October 1984, but better agreement during the later SAGE-II period. For 1963 Agung, the 9 Tg simulation compares best to the forcing datasets with the model capturing the lidar-observed signature of peak extinction descending from 20 km in 1964 to 16 km in 1965. Overall, our results indicate that the downward adjustment to previous SO2 emission estimates for Pinatubo as suggested by several interactive modelling studies is also needed for the Agung and El Chichón aerosol clouds. This strengthens the hypothesis that interactive stratospheric aerosol models may be missing an important removal or redistribution process (e.g. effects of co-emitted ash) which changes how the tropical reservoir of volcanic aerosol evolves in the initial months after an eruption. Our analysis identifies potentially important inhomogeneities in the CMIP6 dataset for all three periods that are hard to reconcile with variations predicted by the interactive stratospheric aerosol model. We also highlight large differences between the CMIP5 and CMIP6 volcanic aerosol datasets for the Agung and El Chichón periods. Future research should aim to reduce this uncertainty by reconciling the datasets with additional stratospheric aerosol observations. [ABSTRACT FROM AUTHOR]

Here we present a description of the UKCA StratTrop chemical mechanism, which is used in the UKESM1 Earth system model for CMIP6. The StratTrop chemical mechanism is a merger of previously well-evaluated tropospheric and stratospheric mechanisms, and we provide results from a series of bespoke integrations to assess the overall performance of the model. We find that the StratTrop scheme performs well when compared to a wide array of observations. The analysis we present here focuses on key components of atmospheric composition, namely the performance of the model to simulate ozone in the stratosphere and troposphere and constituents that are important for ozone in these regions. We find that the results obtained for tropospheric ozone and its budget terms from the use of the StratTrop mechanism are sensitive to the host model; simulations with the same chemical mechanism run in an earlier version of the MetUM host model show a range of sensitivity to emissions that the current model does not fall within. Whilst the general model performance is suitable for use in the UKESM1 CMIP6 integrations, we note some shortcomings in the scheme that future targeted studies will address. [ABSTRACT FROM AUTHOR]

We document the development of the first version of the U.K. Earth System Model UKESM1. The model represents a major advance on its predecessor HadGEM2‐ES, with enhancements to all component models and new feedback mechanisms. These include a new core physical model with a well‐resolved stratosphere; terrestrial biogeochemistry with coupled carbon and nitrogen cycles and enhanced land management; tropospheric‐stratospheric chemistry allowing the holistic simulation of radiative forcing from ozone, methane, and nitrous oxide; two‐moment, five‐species, modal aerosol; and ocean biogeochemistry with two‐way coupling to the carbon cycle and atmospheric aerosols. The complexity of coupling between the ocean, land, and atmosphere physical climate and biogeochemical cycles in UKESM1 is unprecedented for an Earth system model. We describe in detail the process by which the coupled model was developed and tuned to achieve acceptable performance in key physical and Earth system quantities and discuss the challenges involved in mitigating biases in a model with complex connections between its components. Overall, the model performs well, with a stable pre‐industrial state and good agreement with observations in the latter period of its historical simulations. However, global mean surface temperature exhibits stronger‐than‐observed cooling from 1950 to 1970, followed by rapid warming from 1980 to 2014. Metrics from idealized simulations show a high climate sensitivity relative to previous generations of models: Equilibrium climate sensitivity is 5.4 K, transient climate response ranges from 2.68 to 2.85 K, and transient climate response to cumulative emissions is 2.49 to 2.66 K TtC−1. Plain Language Summary: We describe the development and behavior of UKESM1, a novel climate model that includes improved representations of processes in the atmosphere, ocean, and on land. These processes are inter‐related: For example, dust is produced on the land and blown up into the atmosphere where it affects the amount of sunlight falling on Earth. Dust can also be dissolved in the ocean, where it affects marine life. This in turn changes both the amount of carbon dioxide absorbed by the ocean and the material emitted from the surface into the atmosphere, which has an affect on the formation of clouds. UKESM1 includes many processes and interactions such as these, giving it a high level of complexity. Ensuring realistic process behavior is a major challenge in the development of our model, and we have carefully tested this. UKESM1 performs well, correctly exhibiting stable results from a continuous pre‐industrial simulation (used to provide a reference for future experiments) and showing good agreement with observations toward the end of its historical simulations. Results for some properties—including the degree to which average surface temperature changes with increased amounts of carbon dioxide in the atmosphere—are examined in detail. Key Points: UKESM1 represents a major advance over its predecessor HadGEM2‐ES, both in the complexity of its components and its internal couplingThe complex coupling presents challenges to the model development; we document the tuning process employed to obtain acceptable performanceUKESM1 performs well, having a stable pre‐industrial state and showing good agreement with observations in a wide variety of contexts [ABSTRACT FROM AUTHOR]

Here we present a description of the UKCA StratTrop chemical mechanism which is used in the UKESM1 Earth System Model for CMIP6. The StratTrop chemical mechanism is a merger of previously well evaluated tropospheric and stratospheric mechanisms and we provide results from a series of bespoke integrations to assess the overall performance of the model. We find that the StratTrop scheme performs well when compared to a wide array of observations. The analysis we present here focuses on key components of atmospheric composition, namely the performance of the model to simulate ozone in the stratosphere and troposphere and constituents that are important for ozone in these regions. We find that the results obtained from the use of the StratTrop mechanism are sensitive to the host model; simulations with the same chemical mechanism run in an earlier version of the MetUM host model show a range of sensitivity to emissions that the current model does not fall within. Whilst the general model performance is suitable for use in the UKESM1 CMIP6 integrations, we note some shortcomings in the scheme that future targeted studies will address. [ABSTRACT FROM AUTHOR]

The vertical profile of aerosol is important for its radiative effects, but weakly constrained by observations on the global scale, and highly variable among different models. To investigate the controlling factors in one particular model, we investigate the effects of individual processes in HadGEM3-UKCA and compare the resulting diversity of aerosol vertical profiles with the inter-model diversity from the AeroCom Phase II control experiment. In this way we show that (in this model at least) the vertical profile is controlled by a relatively small number of processes, although these vary among aerosol components and particle sizes. We also show that sufficiently coarse variations in these processes can produce a similar diversity to that among different models in terms of the global-mean profile and, to a lesser extent, the zonal-mean vertical position. However, there are features of certain models' profiles that cannot be reproduced, suggesting the influence of further structural differences between models. In HadGEM3-UKCA, convective transport is found to be very important in controlling the vertical profile of all aerosol components by mass. In-cloud scavenging is very important for all except mineral dust. Growth by condensation is important for sulfate and carbonaceous aerosol (along with aqueous oxidation for the former and ageing by soluble material for the latter). The vertical extent of biomass-burning emissions into the free troposphere is also important for the profile of carbonaceous aerosol. Boundary-layer mixing plays a dominant role for sea salt and mineral dust, which are emitted only from the surface. Dry deposition and below-cloud scavenging are important for the profile of mineral dust only. In this model, the microphysical processes of nucleation, condensation and coagulation dominate the vertical profile of the smallest particles by number (e.g. total CN >3 nm), while the profiles of larger particles (e.g. CN>100 nm) are controlled by the same processes as the component mass profiles, plus the size distribution of primary emissions. We also show that the processes that affect the AODnormalised radiative forcing in the model are predominantly those that affect the vertical mass distribution, in particular convective transport, in-cloud scavenging, aqueous oxidation, ageing and the vertical extent of biomass-burning emissions. [ABSTRACT FROM AUTHOR]