Search

The recovery and re-calibration of a dataset of vertical aerosol extinction profiles of the 1963/64 stratospheric aerosol layer measured by a searchlight at 32° N in New Mexico, US, is reported. The recovered dataset consists of 105 aerosol extinction profiles at 550 nm that cover the period from December 1963 to December 1964. It is a unique record of the portion of the aerosol cloud from the March 1963 Agung volcanic eruption that was transported into the Northern Hemisphere subtropics. The data-recovery methodology involved re-digitizing the 105 original aerosol extinction profiles from individual Figures within a research report, followed by the re-calibration. It involves inverting the original equation used to compute the aerosol extinction profile to retrieve the corresponding normalized detector response profile. The re-calibration of the original aerosol extinction profiles used Rayleigh extinction profiles calculated from local soundings. Rayleigh and aerosol slant transmission corrections are applied using the MODTRAN code in transmission mode. Also, a best-estimate aerosol phase function was calculated from observations and applied to the entire column. The tropospheric aerosol phase function from an AERONET station in the vicinity of the searchlight location was applied between 2.76 to 11.7 km. The stratospheric phase function, applied for a 12.2 to 35.2 km altitude range, is calculated from particle-size distributions measured by a high-altitude aircraft in the vicinity of the searchlight in early 1964. The original error estimate was updated considering unaccounted errors. Both the re-calibrated aerosol extinction profiles and the re-calibrated stratospheric aerosol optical depth magnitudes showed higher magnitudes than the original aerosol extinction profiles and the original stratospheric aerosol optical depth, respectively. However, the magnitudes of the re-calibrated variables show a reasonable agreement with other contemporary observations. The re-calibrated stratospheric aerosol optical depth demonstrated its consistency with the tropics-to-pole decreasing trend, associated with the major volcanic eruption stratospheric aerosol pattern when compared to the time-coincident stratospheric aerosol optical depth lidar observations at Lexington at 42° N. [ABSTRACT FROM AUTHOR]

A growing number of general circulation models are adapting interactive sulfur and aerosol schemes to improve the representation of relevant physical and chemical processes and associated feedbacks. They are motivated by investigations of climate response to major volcanic eruptions and potential solar geoengineering scenarios. However, uncertainties in these schemes are not well constrained. Stratospheric sulfate is modulated by emissions of sulfur-containing species of anthropogenic and natural origin, including volcanic activity. While the effects of volcanic eruptions have been studied in the framework of global model intercomparisons, the background conditions of the sulfur cycle have not been addressed in such a way. Here, we fill this gap by analyzing the distribution of the main sulfur species in nine global atmospheric aerosol models for a volcanically quiescent period. We use observational data to evaluate model results. Overall, models agree that the three dominant sulfur species in terms of burdens (sulfate aerosol, OCS, and SO 2) make up about 98 % stratospheric sulfur and 95 % tropospheric sulfur. However, models vary considerably in the partitioning between these species. Models agree that anthropogenic emission of SO 2 strongly affects the sulfate aerosol burden in the northern hemispheric troposphere, while its importance is very uncertain in other regions, where emissions are much lower. Sulfate aerosol is the main deposited species in all models, but the values deviate by a factor of 2. Additionally, the partitioning between wet and dry deposition fluxes is highly model dependent. Inter-model variability in the sulfur species is low in the tropics and increases towards the poles. Differences are largest in the dynamically active northern hemispheric extratropical region and could be attributed to the representation of the stratospheric circulation. The differences in the atmospheric sulfur budget among the models arise from the representation of both chemical and dynamical processes, whose interplay complicates the bias attribution. Several problematic points identified for individual models are related to the specifics of the chemistry schemes, model resolution, and representation of cross-tropopause transport in the extratropics. Further model intercomparison research is needed with a focus on the clarification of the reasons for biases, given the importance of this topic for the stratospheric aerosol injection studies. [ABSTRACT FROM AUTHOR]

Substantial and prolonged enhancements in stratospheric water vapor (SWV) have occurred after large‐magnitude explosive tropical volcanic eruptions, with modified tropopause entry caused by aerosol‐absorptive heating. Here, we analyze the timing and longevity of heating‐driven post‐eruption SWV changes within CMIP6‐VolMIP short‐term climate‐response experiments with the UK Earth System Model (UKESM1). We find aerosol‐absorptive heating causes peak SWV increases of 17% (∼1 ppmv) and 10% (0.5 ppmv) at 100 and 50 hPa, at ∼18 and ∼23 months after a Pinatubo‐like eruption, respectively. We track the temperature response in the tropical lower stratosphere and identify the main SWV increase occurs only after the descending aerosol heating reaches the tropopause, suggesting a key role for aerosol microphysical processes (sedimentation rate). We explore how El Niño–Southern Oscillation variability modulates this effect. Post‐eruption SWV increases are ∼80% stronger for the La Nina phase compared to the ensemble mean. Tropical upwelling strongly mediates this effect. Plain Language Summary: Strong volcanic eruptions, such as the 1991 eruption of Mt Pinatubo, inject a large amount of SO2 directly into the stratosphere, thereby enhancing the stratospheric aerosol layer and causing a short‐term climatic perturbation. Another substantial part of the climatic influence is the change in stratospheric water vapor (SWV), which affects the chemical processes and the radiative budget of the atmosphere. Along with near‐instantaneous injection of water vapor into the stratosphere, volcanic eruptions can indirectly enhance the entry of water vapor into the stratosphere through aerosol‐induced tropopause heating. This work analyses Earth system model experiments designed to explore how volcanic impacts combine with internal climate variability. We find that peak SWV entry mixing ratios occur only within the second post‐eruption year, consistent with the substantially lagged timing of SWV increase seen in post‐Pinatubo satellite measurements. This analysis provides a new perspective on the temporal evolution of the observed post‐Pinatubo SWV increase and an improved quantification of its impacts. Key Points: Aerosol‐induced absorptive‐heating increases stratospheric water vapor (SWV) by up to 17% at 18 months post‐eruption in a Pinatubo‐like experimentAnalyzing simulations by El Niño–Southern Oscillation (ENSO) variability show an 80% larger peak SWV increase occurs if an eruption is followed by a La Niña phaseThe timing of peak SWV increase occurs when volcanic‐aerosol‐induced heating reaches the tropopause, with ENSO modulation of tropical upwelling [ABSTRACT FROM AUTHOR]

A previous model intercomparison of the Tambora aerosol cloud has highlighted substantial differences among simulated volcanic aerosol properties in the pre-industrial stratosphere and has led to questions about the applicability of global aerosol models for large-magnitude explosive eruptions prior to the observational period. Here, we compare the evolution of the stratospheric aerosol cloud following the well-observed June 1991 Mt. Pinatubo eruption simulated with six interactive stratospheric aerosol microphysics models to a range of observational data sets. Our primary focus is on the uncertainties regarding initial SO 2 emission following the Pinatubo eruption, as prescribed in the Historical Eruptions SO 2 Emission Assessment experiments (HErSEA), in the framework of the Interactive Stratospheric Aerosol Model Intercomparison Project (ISA-MIP). Six global models with interactive aerosol microphysics took part in this study: ECHAM6-SALSA, EMAC, ECHAM5-HAM, SOCOL-AERv2, ULAQ-CCM, and UM-UKCA. Model simulations are performed by varying the SO 2 injection amount (ranging between 5 and 10 Tg S) and the altitude of injection (between 18–25 km). The comparisons show that all models consistently demonstrate faster reduction from the peak in sulfate mass burden in the tropical stratosphere. Most models also show a stronger transport towards the extratropics in the Northern Hemisphere, at the expense of the observed tropical confinement, suggesting a much weaker subtropical barrier in all the models, which results in a shorter e-folding time compared to the observations. Furthermore, simulations in which more than 5 Tg S in the form of SO 2 is injected show an initial overestimation of the sulfate burden in the tropics and, in some models, in the Northern Hemisphere and a large surface area density a few months after the eruption compared to the values measured in the tropics and the in situ measurements over Laramie. This draws attention to the importance of including processes such as the ash injection for the removal of the initial SO 2 and aerosol lofting through local heating. [ABSTRACT FROM AUTHOR]

We report the recovery and re-calibration of an extensive dataset of vertical profile measurements of the 1963/64 stratospheric aerosol layer measured from a two-site searchlight measurement facility at White Sands missile base and Sacramento Peak observatory, in New Mexico, US. The recovered dataset comprises 105 profiles of 550nm aerosol extinction (βp(z)) and is part of a longer program of measurements with the US Air Force Cambridge Research Laboratories (AFCRL) searchlight facility that began in February 1963. The recovered series of βp(z) profiles span the 13-month period December 1963 to December 1964 and provide a unique record of the altitude and vertical extent of the Northern Hemisphere dispersed portion of the aerosol cloud from the March 1963 Agung volcanic eruption. The data recovery methodology involved first re-digitizing the 105 original βp(z) profiles (βp(z)) from individual Figures within an AFCRL research report (Elterman, 1966a). The re-calibration involves inverting the original equation used to compute βp(z) in Elterman (1966a; 1966b) to retrieve a normalized detector response (Er(z) Er(35)) profile for each of the 105 re-digitized βp 0(z) profiles. An iterative procedure was then used to compute the re-calibration βp(z) profiles (βp Ra(z)), with the molecular extinction profile calculated with the corresponding daily molecular extinction profile calculated from local soundings, rather than the US Standard Atmosphere 1962 in the original dataset. Two-way molecular and aerosol transmittance corrections are applied using the MODTRAN (MODerate resolution atmospheric TRANsmission) code in transmission mode, applying a best-estimate aerosol phase function calculated from measurements, applied for the entire 2.76 to 35.2 km column. For the tropospheric aerosol transmittance, the AERONET aerosol phase function from White Sands High Energy Laser Systems Test Facility (HELSTF) was applied (2.76 to 10.7 km), a separate stratospheric phase function applied between 11.2 and 35.2 km, calculated from a set of particle size distributions measured by the U-2 high-altitude aircraft over a US region in the vicinity of White Sands in early 1964. Errors were estimated taking as a reference the errors determined in the computation of βp 0(z). Using available tabulated data from the original procedure the errors in the re-digitalization of βp 0(z) and in the retrieval of the Er(z) Er(35) procedures were calculated and later added to the original estimates. Both the βp R(z) and the stratospheric aerosol optical depth magnitudes showed higher magnitudes than βp 0(z) and the original stratospheric aerosol optical depth, however their magnitudes show a reasonable agreement with other contemporary observations. Both the original and recalibration datasets are being submitted to PANGAEA data repository for its storage and public access. [ABSTRACT FROM AUTHOR]

As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), several climate modeling centers performed a coordinated pre-study experiment with interactive stratospheric aerosol models simulating the volcanic aerosol cloud from an eruption resembling the 1815 Mt. Tambora eruption (VolMIP-Tambora ISA ensemble). The pre-study provided the ancillary ability to assess intermodel diversity in the radiative forcing for a large stratospheric-injecting equatorial eruption when the volcanic aerosol cloud is simulated interactively. An initial analysis of the VolMIP-Tambora ISA ensemble showed large disparities between models in the stratospheric global mean aerosol optical depth (AOD). In this study, we now show that stratospheric global mean AOD differences among the participating models are primarily due to differences in aerosol size, which we track here by effective radius. We identify specific physical and chemical processes that are missing in some models and/or parameterized differently between models, which are together causing the differences in effective radius. In particular, our analysis indicates that interactively tracking hydroxyl radical (OH) chemistry following a large volcanic injection of sulfur dioxide (SO2) is an important factor in allowing for the timescale for sulfate formation to be properly simulated. In addition, depending on the timescale of sulfate formation, there can be a large difference in effective radius and subsequently AOD that results from whether the SO2 is injected in a single model grid cell near the location of the volcanic eruption, or whether it is injected as a longitudinally averaged band around the Earth. [ABSTRACT FROM AUTHOR]

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]

The radiative forcing caused by a volcanic eruption is dependent on several eruption source parameters such as the mass of sulfur dioxide (SO2) emitted, the eruption column height, and the eruption latitude. General circulation models with prognostic aerosol and chemistry schemes can be used to investigate how each parameter influences the volcanic forcing. However, the range of multidimensional parameter space that can be explored is restricted because such simulations are computationally expensive. Here we use statistical emulation to explore the radiative impact of eruptions over a wide covarying range of SO2 emission magnitudes, injection heights, and eruption latitudes based on only 30 simulations. We use the emulators to build response surfaces to visualize and predict the sulfate aerosol e‐folding decay time, the stratospheric aerosol optical depth, and net radiative forcing of thousands of different eruptions. We find that the volcanic stratospheric aerosol optical depth and net radiative forcing are primarily determined by the mass of SO2 emitted, but eruption latitude is the most important parameter in determining the sulfate aerosol e‐folding decay time. The response surfaces reveal joint effects of the eruption source parameters in influencing the net radiative forcing, such as a stronger influence of injection height for tropical eruptions than high‐latitude eruptions. We also demonstrate how the emulated response surfaces can be used to find all combinations of eruption source parameters that produce a particular volcanic response, often revealing multiple solutions. Key Points: We demonstrate the feasibility and value of using statistical emulation to quantify the radiative impact of volcanic eruptionsEmulated response surfaces illustrate the dependencies of model output such as net radiative forcing on eruption source parametersEmulated response surfaces can also be used to constrain the eruption source parameters for a particular volcanic response [ABSTRACT FROM AUTHOR]

The Stratospheric Sulfur and its Role in Climate (SSiRC) Interactive Stratospheric Aerosol Model Intercomparison Project (ISA-MIP) explores uncertainties in the processes that connect volcanic emission of sulfur gas species and the radiative forcing associated with the resulting enhancement of the stratospheric aerosol layer. The central aim of ISA-MIP is to constrain and improve interactive stratospheric aerosol models and reduce uncertainties in the stratospheric aerosol forcing by comparing results of standardized model experiments with a range of observations. In this paper we present four co-ordinated inter-model experiments designed to investigate key processes which influence the formation and temporal development of stratospheric aerosol in different time periods of the observational record. The Background (BG) experiment will focus on microphysics and transport processes under volcanically quiescent conditions, when the stratospheric aerosol is controlled by the transport of aerosols and their precursors from the troposphere to the stratosphere. The Transient Aerosol Record (TAR) experiment will explore the role of small- to moderate-magnitude volcanic eruptions, anthropogenic sulfur emissions, and transport processes over the period 1998-2012 and their role in the warming hiatus. Two further experiments will investigate the stratospheric sulfate aerosol evolution after major volcanic eruptions. The Historical Eruptions SO2 Emission Assessment (HErSEA) experiment will focus on the uncertainty in the initial emission of recent large-magnitude volcanic eruptions, while the Pinatubo Emulation in Multiple models (PoEMS) experiment will provide a comprehensive uncertainty analysis of the radiative forcing from the 1991 Mt Pinatubo eruption. [ABSTRACT FROM AUTHOR]

The enhancement of the stratospheric aerosol layer by volcanic eruptions induces a complex set of responses causing global and regional climate effects on a broad range of timescales. Uncertainties exist regarding the climatic response to strong volcanic forcing identified in coupled climate simulations that contributed to the fifth phase of the Coupled Model Intercomparison Project (CMIP5). In order to better understand the sources of these model diversities, the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP) has defined a coordinated set of idealized volcanic perturbation experiments to be carried out in alignment with the CMIP6 protocol. VolMIP provides a common stratospheric aerosol data set for each experiment to minimize differences in the applied volcanic forcing. It defines a set of initial conditions to assess how internal climate variability contributes to determining the response. VolMIP will assess to what extent volcanically forced responses of the coupled ocean-atmosphere system are robustly simulated by state-of-the-art coupled climate models and identify the causes that limit robust simulated behavior, especially differences in the treatment of physical processes. This paper illustrates the design of the idealized volcanic perturbation experiments in the VolMIP protocol and describes the common aerosol forcing input data sets to be used. [ABSTRACT FROM AUTHOR]

The Stratospheric Sulfur and its Role in Climate (SSiRC) interactive stratospheric aerosol model intercomparison project (ISA-MIP) explores uncertainties in the processes that connect volcanic emission of sulphur gas species and the radiative forcing associated with the resulting enhancement of the stratospheric aerosol layer. The central aim of ISA-MIP is to constrain and improve interactive stratospheric aerosol models and reduce uncertainties in the stratospheric aerosol forcing by comparing results of standardized model experiments with a range of observations. In this paper we present 4 co-ordinated inter-model experiments designed to investigate key processes which influence the formation and temporal development of stratospheric aerosol in different time periods of the observational record. The "Background" (BG) experiment will focus on microphysics and transport processes under volcanically quiescent conditions, when the stratospheric aerosol is controlled by the transport of aerosols and their precursors from the troposphere to the stratosphere. The "Transient Aerosol Record" (TAR) experiment will explore the role of small- to moderate-magnitude volcanic eruptions, anthropogenic sulphur emissions and transport processes over the period 1998-2012 and their role in the warming hiatus. Two further experiments will investigate the stratospheric sulphate aerosol evolution after major volcanic eruptions. The "Historical Eruptions SO2 Emission Assessment" (HErSEA) experiment will focus on the uncertainty in the initial emission of recent large-magnitude volcanic eruptions, while the "Pinatubo Emulation in Multiple models" (PoEMS) experiment will provide a comprehensive uncertainty analysis of the radiative forcing from the 1991 Mt. Pinatubo eruption. [ABSTRACT FROM AUTHOR]