Search

Previous research has shown that stratospheric aerosol causes only a small temperature change per unit forcing because they produce stronger cooling in the tropical Indian Ocean and the western Pacific Ocean than in the global mean. The enhanced temperature change in this so-called “warm-pool” region activates strongly negative local and remote feedbacks, which dampen the global mean temperature response. This paper addresses the question of why stratospheric aerosol forcing affects warm-pool temperatures more strongly than CO2 forcing, using idealized MPI-ESM simulations. We show that the aerosol's enhanced effective forcing at the top of the atmosphere (TOA) over the warm pool contributes to the warm-pool-intensified temperature change but is not sufficient to explain the effect. Instead, the pattern of surface effective forcing, which is substantially different from the effective forcing at the TOA, is more closely linked to the temperature pattern. Independent of surface temperature changes, the aerosol heats the tropical stratosphere, accelerating the Brewer–Dobson circulation. The intensified Brewer–Dobson circulation exports additional energy from the tropics to the extratropics, which leads to a particularly strong negative forcing at the tropical surface. These results show how forced circulation changes can affect the climate response by altering the surface forcing pattern. Furthermore, they indicate that the established approach of diagnosing effective forcing at the TOA is useful for global means, but a surface perspective on the forcing must be adopted to understand the evolution of temperature patterns.

Volcanic eruptions impact the climate and environment. The volcanic forcing is determined by eruption source parameters, including the mass and composition of volcanic volatiles, eruption season, eruption latitude, and injection altitude. Moreover, initial atmospheric conditions of the climate system play an important role in shaping the volcanic forcing and response. However, our understanding of the combination of these factors, the distinctions between tropical and extratropical volcanic eruptions, and the co-injection of sulfur and halogens remains limited. Here, we perform ensemble simulations of volcanic eruptions at 15 and 64° N in January, injecting 17 Mt of SO2 together with HCl and HBr at 24 km altitude. Our findings reveal that initial atmospheric conditions control the transport of volcanic volatiles from the first month and modulate the subsequent latitudinal distribution of sulfate aerosols and halogens. This results in different volcanic forcing, surface temperature and ozone responses over the globe and Northern Hemisphere extratropics (NHET) among the model ensemble members with different initial atmospheric conditions. NH extratropical eruptions exhibit a larger NHET mean volcanic forcing, surface cooling and ozone depletion compared with tropical eruptions. However, tropical eruptions lead to more prolonged impacts compared with NH extratropical eruptions, both globally and in the NHET. The sensitivity of volcanic forcing to varying eruption source parameters and model dependency is discussed, emphasizing the need for future multi-model studies to consider the influence of initial conditions and eruption source parameters on volcanic forcing and subsequent impacts.

Perturbations in stratospheric aerosol due to explosive volcanic eruptions are a primary contributor to natural climate variability. Observations of stratospheric aerosol are available for the past decades, and information from ice cores has been used to derive estimates of stratospheric sulfur injections and aerosol optical depth over the Holocene (approximately 10 000 BP to present) and into the last glacial period, extending back to 60 000 BP. Tephra records of past volcanism, compared to ice cores, are less complete but extend much further into the past. To support model studies of the potential impacts of explosive volcanism on climate variability across timescales, we present here an ensemble reconstruction of volcanic stratospheric sulfur injection (VSSI) over the last 140 000 years that is based primarily on terrestrial and marine tephra records. VSSI values are computed as a simple function of eruption magnitude based on VSSI estimates from ice cores and satellite observations for identified eruptions. To correct for the incompleteness of the tephra record, we include stochastically generated synthetic eruptions assuming a constant background eruption frequency from the ice core Holocene record. While the reconstruction often differs from ice core estimates for specific eruptions due to uncertainties in the data used and reconstruction method, it shows good agreement with an ice-core-based VSSI reconstruction in terms of millennial-scale cumulative VSSI variations over the Holocene. The PalVol reconstruction provides a new basis to test the contributions of forced vs. unforced natural variability to the spectrum of climate and the mechanisms leading to abrupt transitions in the palaeoclimate record with low- to high-complexity climate models. The PalVol volcanic forcing reconstruction is available at https://doi.org/10.26050/WDCC/PalVolv1 (Toohey and Schindlbeck-Belo, 2023).

Here we investigate how uncertainties in the solar and volcanic forcing records of the past millennium affect the large-scale temperature response using a two-box impulse response model. We use different published solar forcing records and present a new volcanic forcing ensemble that accounts for random uncertainties in eruption dating and sulfur injection amount. The simulations are compared to proxy reconstructions from PAGES 2k and Northern Hemispheric tree ring data. We find that low solar forcing is most consistent with all the proxy reconstructions, even when accounting for volcanic uncertainty. We also find that the residuals are in line with CMIP6 control variability at centennial timescales. Volcanic forcing uncertainty induces a significant spread in the temperature response, especially at periods of peak forcing. For individual eruptions and superposed epoch analyses, volcanic uncertainty can strongly affect the agreement with proxy reconstructions and partly explain known proxy–model discrepancies.

The injection of sulfur into the stratosphere by volcanic eruptions is the dominant driver of natural climate variability on interannual to multidecadal timescales. Based on a set of continuous sulfate and sulfur records from a suite of ice cores from Greenland and Antarctica, the HolVol v.1.0 database includes estimates of the magnitudes and approximate source latitudes of major volcanic stratospheric sulfur injection (VSSI) events for the Holocene (from 9500 BCE or 11 500 years BP to 1900 CE), constituting an extension of the previous record by 7000 years. The database incorporates new-generation ice-core aerosol records with a sub-annual temporal resolution and a demonstrated sub-decadal dating accuracy and precision. By tightly aligning and stacking the ice-core records on the WD2014 chronology from Antarctica, we resolve long-standing inconsistencies in the dating of ancient volcanic eruptions that arise from biased (i.e., dated too old) ice-core chronologies over the Holocene for Greenland. We reconstruct a total of 850 volcanic eruptions with injections in excess of 1 teragram of sulfur (Tg S); of these eruptions, 329 (39 %) are located in the low latitudes with bipolar sulfate deposition, 426 (50 %) are located in the Northern Hemisphere extratropics (NHET) and 88 (10 %) are located in the Southern Hemisphere extratropics (SHET). The spatial distribution of the reconstructed eruption locations is in agreement with prior reconstructions for the past 2500 years. In total, these eruptions injected 7410 Tg S into the stratosphere: 70 % from tropical eruptions and 25 % from NH extratropical eruptions. A long-term latitudinally and monthly resolved stratospheric aerosol optical depth (SAOD) time series is reconstructed from the HolVol VSSI estimates, representing the first Holocene-scale reconstruction constrained by Greenland and Antarctica ice cores. These new long-term reconstructions of past VSSI and SAOD variability confirm evidence from regional volcanic eruption chronologies (e.g., from Iceland) in showing that the Early Holocene (9500–7000 BCE) experienced a higher number of volcanic eruptions (+16 %) and cumulative VSSI (+86 %) compared with the past 2500 years. This increase coincides with the rapid retreat of ice sheets during deglaciation, providing context for potential future increases in volcanic activity in regions under projected glacier melting in the 21st century. The reconstructed VSSI and SAOD data are available at https://doi.org/10.1594/PANGAEA.928646 (Sigl et al., 2021).

The 852/3 CE eruption of Mount Churchill, Alaska, was one of the largest first-millennium volcanic events, with a magnitude of 6.7 (VEI 6) and a tephra volume of 39.4–61.9 km3 (95 % confidence). The spatial extent of the ash fallout from this event is considerable and the cryptotephra (White River Ash east; WRAe) extends as far as Finland and Poland. Proximal ecosystem and societal disturbances have been linked with this eruption; however, wider eruption impacts on climate and society are unknown. Greenland ice core records show that the eruption occurred in winter 852/3 ± 1 CE and that the eruption is associated with a relatively moderate sulfate aerosol loading but large abundances of volcanic ash and chlorine. Here we assess the potential broader impact of this eruption using palaeoenvironmental reconstructions, historical records and climate model simulations. We also use the fortuitous timing of the 852/3 CE Churchill eruption and its extensively widespread tephra deposition of the White River Ash (east) (WRAe) to examine the climatic expression of the warm Medieval Climate Anomaly period (MCA; ca. 950–1250 CE) from precisely linked peatlands in the North Atlantic region. The reconstructed climate forcing potential of the 852/3 CE Churchill eruption is moderate compared with the eruption magnitude, but tree-ring-inferred temperatures report a significant atmospheric cooling of 0.8 ∘C in summer 853 CE. Modelled climate scenarios also show a cooling in 853 CE, although the average magnitude of cooling is smaller (0.3 ∘C). The simulated spatial patterns of cooling are generally similar to those generated using the tree-ring-inferred temperature reconstructions. Tree-ring-inferred cooling begins prior to the date of the eruption suggesting that natural internal climate variability may have increased the climate system's susceptibility to further cooling. The magnitude of the reconstructed cooling could also suggest that the climate forcing potential of this eruption may be underestimated, thereby highlighting the need for greater insight into, and consideration of, the role of halogens and volcanic ash when estimating eruption climate forcing potential. Precise comparisons of palaeoenvironmental records from peatlands across North America and Europe, facilitated by the presence of the WRAe isochron, reveal no consistent MCA signal. These findings contribute to the growing body of evidence that characterises the MCA hydroclimate as time-transgressive and heterogeneous rather than a well-defined climatic period. The presence of the WRAe isochron also demonstrates that no long-term (multidecadal) climatic or societal impacts from the 852/3 CE Churchill eruption were identified beyond areas proximal to the eruption. Historical evidence in Europe for subsistence crises demonstrate a degree of temporal correspondence on interannual timescales, but similar events were reported outside of the eruption period and were common in the 9th century. The 852/3 CE Churchill eruption exemplifies the difficulties of identifying and confirming volcanic impacts for a single eruption, even when the eruption has a small age uncertainty.

This paper provides initial results from a multi-model ensemble analysis based on the volc-pinatubo-full experiment performed within the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP) as part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). The volc-pinatubo-full experiment is based on an ensemble of volcanic forcing-only climate simulations with the same volcanic aerosol dataset across the participating models (the 1991–1993 Pinatubo period from the CMIP6-GloSSAC dataset). The simulations are conducted within an idealized experimental design where initial states are sampled consistently across models from the CMIP6-piControl simulation providing unperturbed preindustrial background conditions. The multi-model ensemble includes output from an initial set of six participating Earth system models (CanESM5, GISS-E2.1-G, IPSL-CM6A-LR, MIROC-E2SL, MPI-ESM1.2-LR and UKESM1). The results show overall good agreement between the different models on the global and hemispheric scales concerning the surface climate responses, thus demonstrating the overall effectiveness of VolMIP's experimental design. However, small yet significant inter-model discrepancies are found in radiative fluxes, especially in the tropics, that preliminary analyses link with minor differences in forcing implementation; model physics, notably aerosol–radiation interactions; the simulation and sampling of El Niño–Southern Oscillation (ENSO); and, possibly, the simulation of climate feedbacks operating in the tropics. We discuss the volc-pinatubo-full protocol and highlight the advantages of volcanic forcing experiments defined within a carefully designed protocol with respect to emerging modelling approaches based on large ensemble transient simulations. We identify how the VolMIP strategy could be improved in future phases of the initiative to ensure a cleaner sampling protocol with greater focus on the evolving state of ENSO in the pre-eruption period.

Volcanic fallout in polar ice sheets provides important opportunities to date and correlate ice-core records as well as to investigate the environmental impacts of eruptions. Only the geochemical characterization of volcanic ash (tephra) embedded in the ice strata can confirm the source of the eruption, however, and is a requisite if historical eruption ages are to be used as valid chronological checks on annual ice layer counting. Here we report the investigation of ash particles in a Greenland ice core that are associated with a volcanic sulfuric acid layer previously attributed to the 79 CE eruption of Vesuvius. Major and trace element composition of the particles indicates that the tephra does not derive from Vesuvius but most likely originates from an unidentified eruption in the Aleutian arc. Using ash dispersal modeling, we find that only an eruption large enough to include stratospheric injection is likely to account for the sizable (24–85 µm) ash particles observed in the Greenland ice at this time. Despite its likely explosivity, this event does not appear to have triggered significant climate perturbations, unlike some other large extratropical eruptions. In light of a recent re-evaluation of the Greenland ice-core chronologies, our findings further challenge the previous assignation of this volcanic event to 79 CE. We highlight the need for the revised Common Era ice-core chronology to be formally accepted by the wider ice-core and climate modeling communities in order to ensure robust age linkages to precisely dated historical and paleoclimate proxy records.

We present a transient simulation of global vegetation and climate patterns of the mid- and late Holocene using the MPI-ESM (Max Planck Institute for Meteorology Earth System Model) at T63 resolution. The simulated vegetation trend is discussed in the context of the simulated Holocene climate change. Our model captures the main trends found in reconstructions. Most prominent are the southward retreat of the northern treeline that is combined with the strong decrease of forest in the high northern latitudes during the Holocene and the vast increase of the Saharan desert, embedded in a general decrease in precipitation and vegetation in the Northern Hemisphere monsoon margin regions. The Southern Hemisphere experiences weaker changes in total vegetation cover during the last 8000 years. However, the monsoon-related increase in precipitation and the insolation-induced cooling of the winter climate lead to shifts in the vegetation composition, mainly between the woody plant functional types (PFTs). The large-scale global patterns of vegetation almost linearly follow the subtle, approximately linear, orbital forcing. In some regions, however, non-linear, more rapid changes in vegetation are found in the simulation. The most striking region is the Sahel–Sahara domain with rapid vegetation transitions to a rather desertic state, despite a gradual insolation forcing. Rapid shifts in the simulated vegetation also occur in the high northern latitudes, in South Asia and in the monsoon margins of the Southern Hemisphere. These rapid changes are mainly triggered by changes in the winter temperatures, which go into, or move out of, the bioclimatic tolerance range of individual PFTs. The dynamics of the transitions are determined by dynamics of the net primary production (NPP) and the competition between PFTs. These changes mainly occur on timescales of centuries. More rapid changes in PFTs that occur within a few decades are mainly associated with the timescales of mortality and the bioclimatic thresholds implicit in the dynamic vegetation model, which have to be interpreted with caution. Most of the simulated Holocene vegetation changes outside the high northern latitudes are associated with modifications in the intensity of the global summer monsoon dynamics that also affect the circulation in the extra tropics via teleconnections. Based on our simulations, we thus identify the global monsoons as the key player in Holocene climate and vegetation change.

The “1809 eruption” is one of the most recent unidentified volcanic eruptions with a global climate impact. Even though the eruption ranks as the third largest since 1500 with a sulfur emission strength estimated to be 2 times that of the 1991 eruption of Pinatubo, not much is known of it from historic sources. Based on a compilation of instrumental and reconstructed temperature time series, we show here that tropical temperatures show a significant drop in response to the ∼ 1809 eruption that is similar to that produced by the Mt. Tambora eruption in 1815, while the response of Northern Hemisphere (NH) boreal summer temperature is spatially heterogeneous. We test the sensitivity of the climate response simulated by the MPI Earth system model to a range of volcanic forcing estimates constructed using estimated volcanic stratospheric sulfur injections (VSSIs) and uncertainties from ice-core records. Three of the forcing reconstructions represent a tropical eruption with an approximately symmetric hemispheric aerosol spread but different forcing magnitudes, while a fourth reflects a hemispherically asymmetric scenario without volcanic forcing in the NH extratropics. Observed and reconstructed post-volcanic surface NH summer temperature anomalies lie within the range of all the scenario simulations. Therefore, assuming the model climate sensitivity is correct, the VSSI estimate is accurate within the uncertainty bounds. Comparison of observed and simulated tropical temperature anomalies suggests that the most likely VSSI for the 1809 eruption would be somewhere between 12 and 19 Tg of sulfur. Model results show that NH large-scale climate modes are sensitive to both volcanic forcing strength and its spatial structure. While spatial correlations between the N-TREND NH temperature reconstruction and the model simulations are weak in terms of the ensemble-mean model results, individual model simulations show good correlation over North America and Europe, suggesting the spatial heterogeneity of the 1810 cooling could be due to internal climate variability.

The Stratosphere-troposphere Processes and their Role in Climate (SPARC) Data Initiative (SPARC, 2017) performed the first comprehensive assessment of currently available stratospheric composition measurements obtained from an international suite of space-based limb sounders. The initiative's main objectives were (1) to assess the state of data availability, (2) to compile time series of vertically resolved, zonal monthly mean trace gas and aerosol fields, and (3) to perform a detailed intercomparison of these time series, summarizing useful information and highlighting differences among datasets. The datasets extend over the region from the upper troposphere to the lower mesosphere (300–0.1 hPa) and are provided on a common latitude–pressure grid. They cover 26 different atmospheric constituents including the stratospheric trace gases of primary interest, ozone (O3) and water vapor (H2O), major long-lived trace gases (SF6, N2O, HF, CCl3F, CCl2F2, NOy), trace gases with intermediate lifetimes (HCl, CH4, CO, HNO3), and shorter-lived trace gases important to stratospheric chemistry including nitrogen-containing species (NO, NO2, NOx, N2O5, HNO4), halogens (BrO, ClO, ClONO2, HOCl), and other minor species (OH, HO2, CH2O, CH3CN), and aerosol. This overview of the SPARC Data Initiative introduces the updated versions of the SPARC Data Initiative time series for the extended time period 1979–2018 and provides information on the satellite instruments included in the assessment: LIMS, SAGE I/II/III, HALOE, UARS-MLS, POAM II/III, OSIRIS, SMR, MIPAS, GOMOS, SCIAMACHY, ACE-FTS, ACE-MAESTRO, Aura-MLS, HIRDLS, SMILES, and OMPS-LP. It describes the Data Initiative's top-down climatological validation approach to compare stratospheric composition measurements based on zonal monthly mean fields, which provides upper bounds to relative inter-instrument biases and an assessment of how well the instruments are able to capture geophysical features of the stratosphere. An update to previously published evaluations of O3 and H2O monthly mean time series is provided. In addition, example trace gas evaluations of methane (CH4), carbon monoxide (CO), a set of nitrogen species (NO, NO2, and HNO3), the reactive nitrogen family (NOy), and hydroperoxyl (HO2) are presented. The results highlight the quality, strengths and weaknesses, and representativeness of the different datasets. As a summary, the current state of our knowledge of stratospheric composition and variability is provided based on the overall consistency between the datasets. As such, the SPARC Data Initiative datasets and evaluations can serve as an atlas or reference of stratospheric composition and variability during the “golden age” of atmospheric limb sounding. The updated SPARC Data Initiative zonal monthly mean time series for each instrument are publicly available and accessible via the Zenodo data archive (Hegglin et al., 2020).

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.

The atmospheric CO2 concentration increased by about 20 ppm from 6000 BCE to the pre-industrial period (1850 CE). Several hypotheses have been proposed to explain mechanisms of this CO2 growth based on either ocean or land carbon sources. Here, we apply the Earth system model MPI-ESM-LR for two transient simulations of climate and carbon cycle dynamics during this period. In the first simulation, atmospheric CO2 is prescribed following ice-core CO2 data. In response to the growing atmospheric CO2 concentration, land carbon storage increases until 2000 BCE, stagnates afterwards, and decreases from 1 CE, while the ocean continuously takes CO2 out of the atmosphere after 4000 BCE. This leads to a missing source of 166 Pg of carbon in the ocean–land–atmosphere system by the end of the simulation. In the second experiment, we applied a CO2 nudging technique using surface alkalinity forcing to follow the reconstructed CO2 concentration while keeping the carbon cycle interactive. In that case the ocean is a source of CO2 from 6000 to 2000 BCE due to a decrease in the surface ocean alkalinity. In the prescribed CO2 simulation, surface alkalinity declines as well. However, it is not sufficient to turn the ocean into a CO2 source. The carbonate ion concentration in the deep Atlantic decreases in both the prescribed and the interactive CO2 simulations, while the magnitude of the decrease in the prescribed CO2 experiment is underestimated in comparison with available proxies. As the land serves as a carbon sink until 2000 BCE due to natural carbon cycle processes in both experiments, the missing source of carbon for land and atmosphere can only be attributed to the ocean. Within our model framework, an additional mechanism, such as surface alkalinity decrease, for example due to unaccounted for carbonate accumulation processes on shelves, is required for consistency with ice-core CO2 data. Consequently, our simulations support the hypothesis that the ocean was a source of CO2 until the late Holocene when anthropogenic CO2 sources started to affect atmospheric CO2.

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.

The eruption of Mt. Tambora in 1815 was the largest volcanic eruption of the past 500 years. The eruption had significant climatic impacts, leading to the 1816 year without a summer, and remains a valuable event from which to understand the climatic effects of large stratospheric volcanic sulfur dioxide injections. The eruption also resulted in one of the strongest and most easily identifiable volcanic sulfate signals in polar ice cores, which are widely used to reconstruct the timing and atmospheric sulfate loading of past eruptions. As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), five state-of-the-art global aerosol models simulated this eruption. We analyse both simulated background (no Tambora) and volcanic (with Tambora) sulfate deposition to polar regions and compare to ice core records. The models simulate overall similar patterns of background sulfate deposition, although there are differences in regional details and magnitude. However, the volcanic sulfate deposition varies considerably between the models with differences in timing, spatial pattern and magnitude. Mean simulated deposited sulfate on Antarctica ranges from 19 to 264 kg km−2 and on Greenland from 31 to 194 kg km−2, as compared to the mean ice-core-derived estimates of roughly 50 kg km−2 for both Greenland and Antarctica. The ratio of the hemispheric atmospheric sulfate aerosol burden after the eruption to the average ice sheet deposited sulfate varies between models by up to a factor of 15. Sources of this inter-model variability include differences in both the formation and the transport of sulfate aerosol. Our results suggest that deriving relationships between sulfate deposited on ice sheets and atmospheric sulfate burdens from model simulations may be associated with greater uncertainties than previously thought.

The injection of sulfur into the stratosphere by explosive volcanic eruptions is the cause of significant climate variability. Based on sulfate records from a suite of ice cores from Greenland and Antarctica, the eVolv2k database includes estimates of the magnitudes and approximate source latitudes of major volcanic stratospheric sulfur injection (VSSI) events from 500 BCE to 1900 CE, constituting an update of prior reconstructions and an extension of the record by 1000 years. The database incorporates improvements to the ice core records (in terms of synchronisation and dating) and refinements to the methods used to estimate VSSI from ice core records, and it includes first estimates of the random uncertainties in VSSI values. VSSI estimates for many of the largest eruptions, including Samalas (1257), Tambora (1815), and Laki (1783), are within 10 % of prior estimates. A number of strong events are included in eVolv2k which are largely underestimated or not included in earlier VSSI reconstructions, including events in 540, 574, 682, and 1108 CE. The long-term annual mean VSSI from major volcanic eruptions is estimated to be ∼ 0.5 Tg [S] yr−1, ∼ 50 % greater than a prior reconstruction due to the identification of more events and an increase in the magnitude of many intermediate events. A long-term latitudinally and monthly resolved stratospheric aerosol optical depth (SAOD) time series is reconstructed from the eVolv2k VSSI estimates, and the resulting global mean SAOD is found to be similar (within 33 %) to a prior reconstruction for most of the largest eruptions. The long-term (500 BCE–1900 CE) average global mean SAOD estimated from the eVolv2k VSSI estimates including a constant background injection of stratospheric sulfur is ∼ 0.014, 30 % greater than a prior reconstruction. These new long-term reconstructions of past VSSI and SAOD variability give context to recent volcanic forcing, suggesting that the 20th century was a period of somewhat weaker than average volcanic forcing, with current best estimates of 20th century mean VSSI and SAOD values being 25 and 14 % less, respectively, than the mean of the 500 BCE to 1900 CE period. The reconstructed VSSI and SAOD data are available at https://doi.org/10.1594/WDCC/eVolv2k_v2.

The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data).

Halogenated very short-lived substances (VSLSs) are naturally produced in the ocean and emitted to the atmosphere. When transported to the stratosphere, these compounds can have a significant influence on the ozone layer and climate. During a research cruise on RV Sonne in the subtropical and tropical west Indian Ocean in July and August 2014, we measured the VSLSs, methyl iodide (CH3I) and for the first time bromoform (CHBr3) and dibromomethane (CH2Br2), in surface seawater and the marine atmosphere to derive their emission strengths. Using the Lagrangian particle dispersion model FLEXPART with ERA-Interim meteorological fields, we calculated the direct contribution of observed VSLS emissions to the stratospheric halogen burden during the Asian summer monsoon. Furthermore, we compare the in situ calculations with the interannual variability of transport from a larger area of the west Indian Ocean surface to the stratosphere for July 2000–2015. We found that the west Indian Ocean is a strong source for CHBr3 (910 pmol m−2 h−1), very strong source for CH2Br2 (930 pmol m−2 h−1), and an average source for CH3I (460 pmol m−2 h−1). The atmospheric transport from the tropical west Indian Ocean surface to the stratosphere experiences two main pathways. On very short timescales, especially relevant for the shortest-lived compound CH3I (3.5 days lifetime), convection above the Indian Ocean lifts oceanic air masses and VSLSs towards the tropopause. On a longer timescale, the Asian summer monsoon circulation transports oceanic VSLSs towards India and the Bay of Bengal, where they are lifted with the monsoon convection and reach stratospheric levels in the southeastern part of the Asian monsoon anticyclone. This transport pathway is more important for the longer-lived brominated compounds (17 and 150 days lifetime for CHBr3 and CH2Br2). The entrainment of CHBr3 and CH3I from the west Indian Ocean to the stratosphere during the Asian summer monsoon is lower than from previous cruises in the tropical west Pacific Ocean during boreal autumn and early winter but higher than from the tropical Atlantic during boreal summer. In contrast, the projected CH2Br2 entrainment was very high because of the high emissions during the west Indian Ocean cruise. The 16-year July time series shows highest interannual variability for the shortest-lived CH3I and lowest for the longest-lived CH2Br2. During this time period, a small increase in VSLS entrainment from the west Indian Ocean through the Asian monsoon to the stratosphere is found. Overall, this study confirms that the subtropical and tropical west Indian Ocean is an important source region of halogenated VSLSs, especially CH2Br2, to the troposphere and stratosphere during the Asian summer monsoon.

Stratospheric sulfate aerosols from volcanic eruptions have a significant impact on the Earth's climate. To include the effects of volcanic eruptions in climate model simulations, the Easy Volcanic Aerosol (EVA) forcing generator provides stratospheric aerosol optical properties as a function of time, latitude, height, and wavelength for a given input list of volcanic eruption attributes. EVA is based on a parameterized three-box model of stratospheric transport and simple scaling relationships used to derive mid-visible (550 nm) aerosol optical depth and aerosol effective radius from stratospheric sulfate mass. Precalculated look-up tables computed from Mie theory are used to produce wavelength-dependent aerosol extinction, single scattering albedo, and scattering asymmetry factor values. The structural form of EVA and the tuning of its parameters are chosen to produce best agreement with the satellite-based reconstruction of stratospheric aerosol properties following the 1991 Pinatubo eruption, and with prior millennial-timescale forcing reconstructions, including the 1815 eruption of Tambora. EVA can be used to produce volcanic forcing for climate models which is based on recent observations and physical understanding but internally self-consistent over any timescale of choice. In addition, EVA is constructed so as to allow for easy modification of different aspects of aerosol properties, in order to be used in model experiments to help advance understanding of what aspects of the volcanic aerosol are important for the climate system.

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.

Observations and simple theoretical arguments suggest that the Northern Hemisphere (NH) stratospheric polar vortex is stronger in winters following major volcanic eruptions. However, recent studies show that climate models forced by prescribed volcanic aerosol fields fail to reproduce this effect. We investigate the impact of volcanic aerosol forcing on stratospheric dynamics, including the strength of the NH polar vortex, in ensemble simulations with the Max Planck Institute Earth System Model. The model is forced by four different prescribed forcing sets representing the radiative properties of stratospheric aerosol following the 1991 eruption of Mt. Pinatubo: two forcing sets are based on observations, and are commonly used in climate model simulations, and two forcing sets are constructed based on coupled aerosol–climate model simulations. For all forcings, we find that simulated temperature and zonal wind anomalies in the NH high latitudes are not directly impacted by anomalous volcanic aerosol heating. Instead, high-latitude effects result from enhancements in stratospheric residual circulation, which in turn result, at least in part, from enhanced stratospheric wave activity. High-latitude effects are therefore much less robust than would be expected if they were the direct result of aerosol heating. Both observation-based forcing sets result in insignificant changes in vortex strength. For the model-based forcing sets, the vortex response is found to be sensitive to the structure of the forcing, with one forcing set leading to significant strengthening of the polar vortex in rough agreement with observation-based expectations. Differences in the dynamical response to the forcing sets imply that reproducing the polar vortex responses to past eruptions, or predicting the response to future eruptions, depends on accurate representation of the space–time structure of the volcanic aerosol forcing.

Climatologies of atmospheric observations are often produced by binning measurements according to latitude and calculating zonal means. The uncertainty in these climatological means is characterised by the standard error of the mean (SEM). However, the usual estimator of the SEM, i.e., the sample standard deviation divided by the square root of the sample size, holds only for uncorrelated randomly sampled measurements. Measurements of the atmospheric state along a satellite orbit cannot always be considered as independent because (a) the time-space interval between two nearest observations is often smaller than the typical scale of variations in the atmospheric state, and (b) the regular time-space sampling pattern of a satellite instrument strongly deviates from random sampling. We have developed a numerical experiment where global chemical fields from a chemistry climate model are sampled according to real sampling patterns of satellite-borne instruments. As case studies, the model fields are sampled using sampling patterns of the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) and Atmospheric Chemistry Experiment Fourier-Transform Spectrometer (ACE-FTS) satellite instruments. Through an iterative subsampling technique, and by incorporating information on the random errors of the MIPAS and ACE-FTS measurements, we produce empirical estimates of the standard error of monthly mean zonal mean model O3 in 5° latitude bins. We find that generally the classic SEM estimator is a conservative estimate of the SEM, i.e., the empirical SEM is often less than or approximately equal to the classic estimate. Exceptions occur only when natural variability is larger than the random measurement error, and specifically in instances where the zonal sampling distribution shows non-uniformity with a similar zonal structure as variations in the sampled field, leading to maximum sensitivity to arbitrary phase shifts between the sample distribution and sampled field. The occurrence of such instances is thus very sensitive to slight changes in the sampling distribution, and to the variations in the measured field. This study highlights the need for caution in the interpretation of the oft-used classically computed SEM, and outlines a relatively simple methodology that can be used to assess one component of the uncertainty in monthly mean zonal mean climatologies produced from measurements from satellite-borne instruments.

Simulations of tropical volcanic eruptions using a general circulation model with coupled aerosol microphysics are used to assess the influence of season of eruption on the aerosol evolution and radiative impacts at the Earth's surface. This analysis is presented for eruptions with SO2 injection magnitudes of 17 and 700 Tg, the former consistent with estimates of the 1991 Mt. Pinatubo eruption, the later a near-"super eruption". For each eruption magnitude, simulations are performed with eruptions at 15° N, at four equally spaced times of year. Sensitivity to eruption season of aerosol optical depth (AOD), clear-sky and all-sky shortwave (SW) radiative flux is quantified by first integrating each field for four years after the eruption, then calculating for each cumulative field the absolute or percent difference between the maximum and minimum response from the four eruption seasons. Eruption season has a significant influence on AOD and clear-sky SW radiative flux anomalies for both eruption magnitudes. The sensitivity to eruption season for both fields is generally weak in the tropics, but increases in the mid- and high latitudes, reaching maximum values of ~75 %. Global mean AOD and clear-sky SW anomalies show sensitivity to eruption season on the order of 15–20 %, which results from differences in aerosol effective radius for the different eruption seasons. Smallest aerosol size and largest cumulative impact result from a January eruption for Pinatubo-magnitude eruption, and from a July eruption for the near-super eruption. In contrast to AOD and clear-sky SW anomalies, all-sky SW anomalies are found to be insensitive to season of eruption for the Pinatubo-magnitude eruption experiment, due to the reflection of solar radiation by clouds in the mid- to high latitudes. However, differences in all-sky SW anomalies between eruptions in different seasons are significant for the larger eruption magnitude, and the ~15 % sensitivity to eruption season of the global mean all-sky SW anomalies is comparable to the sensitivity of global mean AOD and clear-sky SW anomalies. Our estimates of sensitivity to eruption season are larger than previously reported estimates: implications regarding volcanic AOD timeseries reconstructions and their use in climate models are discussed.

Vertical profiles of NO2 and NO have been obtained from solar occultation measurements by the Atmospheric Chemistry Experiment (ACE), using an infrared Fourier Transform Spectrometer (ACE-FTS) and (for NO2) an ultraviolet-visible-near-infrared spectrometer, MAESTRO (Measurement of Aerosol Extinction in the Stratosphere and Troposphere Retrieved by Occultation). In this paper, the quality of the ACE-FTS version 2.2 NO2 and NO and the MAESTRO version 1.2 NO2 data are assessed using other solar occultation measurements (HALOE, SAGE II, SAGE III, POAM III, SCIAMACHY), stellar occultation measurements (GOMOS), limb measurements (MIPAS, OSIRIS), nadir measurements (SCIAMACHY), balloon-borne measurements (SPIRALE, SAOZ) and ground-based measurements (UV-VIS, FTIR). Time differences between the comparison measurements were reduced using either a tight coincidence criterion, or where possible, chemical box models. ACE-FTS NO2 and NO and the MAESTRO NO2 are generally consistent with the correlative data. The ACE-FTS and MAESTRO NO2 volume mixing ratio (VMR) profiles agree with the profiles from other satellite data sets to within about 20% between 25 and 40 km, with the exception of MIPAS ESA (for ACE-FTS) and SAGE II (for ACE-FTS (sunrise) and MAESTRO) and suggest a negative bias between 23 and 40 km of about 10%. MAESTRO reports larger VMR values than the ACE-FTS. In comparisons with HALOE, ACE-FTS NO VMRs typically (on average) agree to ±8% from 22 to 64 km and to +10% from 93 to 105 km, with maxima of 21% and 36%, respectively. Partial column comparisons for NO2 show that there is quite good agreement between the ACE instruments and the FTIRs, with a mean difference of +7.3% for ACE-FTS and +12.8% for MAESTRO.

The Atmospheric Chemistry Experiment (ACE), also known as SCISAT, was launched on 12 August 2003, carrying two instruments that measure vertical profiles of atmospheric constituents using the solar occultation technique. One of these instruments, the ACE Fourier Transform Spectrometer (ACE-FTS), is measuring volume mixing ratio (VMR) profiles of nitrous oxide (N2O) from the upper troposphere to the lower mesosphere at a vertical resolution of about 3–4 km. In this study, the quality of the ACE-FTS version 2.2 N2O data is assessed through comparisons with coincident measurements made by other satellite, balloon-borne, aircraft, and ground-based instruments. These consist of vertical profile comparisons with the SMR, MLS, and MIPAS satellite instruments, multiple aircraft flights of ASUR, and single balloon flights of SPIRALE and FIRS-2, and partial column comparisons with a network of ground-based Fourier Transform InfraRed spectrometers (FTIRs). Between 6 and 30 km, the mean absolute differences for the satellite comparisons lie between −42 ppbv and +17 ppbv, with most within ±20 ppbv. This corresponds to relative deviations from the mean that are within ±15%, except for comparisons with MIPAS near 30 km, for which they are as large as 22.5%. Between 18 and 30 km, the mean absolute differences for the satellite comparisons are generally within ±10 ppbv. From 30 to 60 km, the mean absolute differences are within ±4 ppbv, and are mostly between −2 and +1 ppbv. Given the small N2O VMR in this region, the relative deviations from the mean are therefore large at these altitudes, with most suggesting a negative bias in the ACE-FTS data between 30 and 50 km. In the comparisons with the FTIRs, the mean relative differences between the ACE-FTS and FTIR partial columns (which cover a mean altitude range of 14 to 27 km) are within ±5.6% for eleven of the twelve contributing stations. This mean relative difference is negative at ten stations, suggesting a small negative bias in the ACE-FTS partial columns over the altitude regions compared. Excellent correlation (R=0.964) is observed between the ACE-FTS and FTIR partial columns, with a slope of 1.01 and an intercept of −0.20 on the line fitted to the data.

The Atmospheric Chemistry Experiment (ACE) satellite was launched on 12 August 2003. Its two instruments measure vertical profiles of over 30 atmospheric trace gases by analyzing solar occultation spectra in the ultraviolet/visible and infrared wavelength regions. The reservoir gases HNO3, ClONO2, and N2O5 are three of the key species provided by the primary instrument, the ACE Fourier Transform Spectrometer (ACE-FTS). This paper describes the ACE-FTS version 2.2 data products, including the N2O5 update, for the three species and presents validation comparisons with available observations. We have compared volume mixing ratio (VMR) profiles of HNO3, ClONO2, and N2O5 with measurements by other satellite instruments (SMR, MLS, MIPAS), aircraft measurements (ASUR), and single balloon-flights (SPIRALE, FIRS-2). Partial columns of HNO3 and ClONO2 were also compared with measurements by ground-based Fourier Transform Infrared (FTIR) spectrometers. Overall the quality of the ACE-FTS v2.2 HNO3 VMR profiles is good from 18 to 35 km. For the statistical satellite comparisons, the mean absolute differences are generally within ±1 ppbv ±20%) from 18 to 35 km. For MIPAS and MLS comparisons only, mean relative differences lie within±10% between 10 and 36 km. ACE-FTS HNO3 partial columns (~15–30 km) show a slight negative bias of −1.3% relative to the ground-based FTIRs at latitudes ranging from 77.8° S–76.5° N. Good agreement between ACE-FTS ClONO2 and MIPAS, using the Institut für Meteorologie und Klimaforschung and Instituto de Astrofísica de Andalucía (IMK-IAA) data processor is seen. Mean absolute differences are typically within ±0.01 ppbv between 16 and 27 km and less than +0.09 ppbv between 27 and 34 km. The ClONO2 partial column comparisons show varying degrees of agreement, depending on the location and the quality of the FTIR measurements. Good agreement was found for the comparisons with the midlatitude Jungfraujoch partial columns for which the mean relative difference is 4.7%. ACE-FTS N2O5 has a low bias relative to MIPAS IMK-IAA, reaching −0.25 ppbv at the altitude of the N2O5 maximum (around 30 km). Mean absolute differences at lower altitudes (16–27 km) are typically −0.05 ppbv for MIPAS nighttime and ±0.02 ppbv for MIPAS daytime measurements.

Low-resolution atmospheric thermal emission spectra collected by balloon-borne radiometers over the time span of 1990–2002 are used to retrieve vertical profiles of HNO3, CFC-11 and CFC-12 volume mixing ratios between approximately 10 and 35 km altitude. All of the data analyzed have been collected from launches from a Northern Hemisphere mid-latitude site, during late summer, when stratospheric dynamic variability is at a minimum. The retrieval technique incorporates detailed forward modeling of the instrument and the radiative properties of the atmosphere, and obtains a best fit between modeled and measured spectra through a combination of onion-peeling and optimization steps. The retrieved HNO3 profiles are consistent over the 12-year period, and are consistent with recent measurements by the Atmospheric Chemistry Experiment-Fourier transform spectrometer satellite instrument. We therefore find no evidence of long-term changes in the HNO3 summer mid-latitude profile, although the uncertainty of our measurements precludes a conclusive trend analysis.