scholarly journals Aerosol microphysics simulations of the Mt. Pinatubo eruption with the UKCA composition-climate model

2014 ◽  
Vol 14 (2) ◽  
pp. 2799-2855 ◽  
Author(s):  
S. S. Dhomse ◽  
K. M. Emmerson ◽  
G. W. Mann ◽  
N. Bellouin ◽  
K. S. Carslaw ◽  
...  

Abstract. We have enhanced the capability of a microphysical aerosol-chemistry module to simulate the atmospheric aerosol and precursor gases for both tropospheric and stratospheric conditions. Using the Mount Pinatubo eruption (June 1991) as a test case, we evaluate simulated aerosol properties in a composition-climate model against a range of satellite and in-situ observations. Simulations are performed assuming an injection of 20 Tg SO2 at 19–27 km in tropical latitudes, without any radiative feedback from the simulated aerosol. In both quiescent and volcanically perturbed conditions, simulated aerosol properties in the lower stratosphere show reasonable agreement with the observations. The model captures the observed timing of the maximum aerosol optical depth (AOD) and its decay timescale in both tropics and Northern Hemisphere (NH) mid-latitudes. There is also good qualitative agreement with the observations in terms of spatial and temporal variation of the aerosol effective radius (Reff), which peaks 6–8 months after the eruption. However, the model shows significant biases against some observational data sets. Simulated AOD and Surface Area Density (SAD) in the tropics are substantially higher than the gap-filled satellite data products during the first 6 months after the eruption. The model shows consistently weaker enhancement in Reff compared to satellite and in-situ measurements. Simulated aerosol particle size distribution is also compared to NH mid-latitude in-situ balloon sounding measurements of size-resolved number concentrations. Before the eruption, the model captures the observed profiles of lower stratospheric particle number concentrations with radii larger than 5, 150 and 250 nm (N5, N150 and N250) very well. However, in the first 6 months after the eruption, the model shows high bias in N5 concentrations in the lower stratosphere, suggesting too strong nucleation. Following particle growth via condensation and coagulation, this bias in the finest particles propagates into a factor 2 high bias in N150. Our comparison suggests that new particle formation in the initial phase of large eruptions, and subsequent particle growth to optically-active sizes, might be playing an important role in determining the magnitude of the climate impacts from volcanoes like Pinatubo.

2014 ◽  
Vol 14 (20) ◽  
pp. 11221-11246 ◽  
Author(s):  
S. S. Dhomse ◽  
K. M. Emmerson ◽  
G. W. Mann ◽  
N. Bellouin ◽  
K. S. Carslaw ◽  
...  

Abstract. We use a stratosphere–troposphere composition–climate model with interactive sulfur chemistry and aerosol microphysics, to investigate the effect of the 1991 Mount Pinatubo eruption on stratospheric aerosol properties. Satellite measurements indicate that shortly after the eruption, between 14 and 23 Tg of SO2 (7 to 11.5 Tg of sulfur) was present in the tropical stratosphere. Best estimates of the peak global stratospheric aerosol burden are in the range 19 to 26 Tg, or 3.7 to 6.7 Tg of sulfur assuming a composition of between 59 and 77 % H2SO4. In light of this large uncertainty range, we performed two main simulations with 10 and 20 Tg of SO2 injected into the tropical lower stratosphere. Simulated stratospheric aerosol properties through the 1991 to 1995 period are compared against a range of available satellite and in situ measurements. Stratospheric aerosol optical depth (sAOD) and effective radius from both simulations show good qualitative agreement with the observations, with the timing of peak sAOD and decay timescale matching well with the observations in the tropics and mid-latitudes. However, injecting 20 Tg gives a factor of 2 too high stratospheric aerosol mass burden compared to the satellite data, with consequent strong high biases in simulated sAOD and surface area density, with the 10 Tg injection in much better agreement. Our model cannot explain the large fraction of the injected sulfur that the satellite-derived SO2 and aerosol burdens indicate was removed within the first few months after the eruption. We suggest that either there is an additional alternative loss pathway for the SO2 not included in our model (e.g. via accommodation into ash or ice in the volcanic cloud) or that a larger proportion of the injected sulfur was removed via cross-tropopause transport than in our simulations. We also critically evaluate the simulated evolution of the particle size distribution, comparing in detail to balloon-borne optical particle counter (OPC) measurements from Laramie, Wyoming, USA (41° N). Overall, the model captures remarkably well the complex variations in particle concentration profiles across the different OPC size channels. However, for the 19 to 27 km injection height-range used here, both runs have a modest high bias in the lowermost stratosphere for the finest particles (radii less than 250 nm), and the decay timescale is longer in the model for these particles, with a much later return to background conditions. Also, whereas the 10 Tg run compared best to the satellite measurements, a significant low bias is apparent in the coarser size channels in the volcanically perturbed lower stratosphere. Overall, our results suggest that, with appropriate calibration, aerosol microphysics models are capable of capturing the observed variation in particle size distribution in the stratosphere across both volcanically perturbed and quiescent conditions. Furthermore, additional sensitivity simulations suggest that predictions with the models are robust to uncertainties in sub-grid particle formation and nucleation rates in the stratosphere.


2017 ◽  
Author(s):  
Laura Revell ◽  
Andrea Stenke ◽  
Beiping Luo ◽  
Stefanie Kremser ◽  
Eugene Rozanov ◽  
...  

Abstract. To simulate the impacts of volcanic eruptions on the stratosphere, chemistry-climate models that do not include an online aerosol module require temporally and spatially resolved aerosol size parameters for heterogeneous chemistry and aerosol radiative properties as a function of wavelength. For phase 1 of the Chemistry-Climate Model Initiative (CCMI-1) and, later, for phase 6 of the Coupled Model Intercomparison Project (CMIP6) two such stratospheric aerosol data sets were compiled, whose functional capability and representativeness are compared here. For CCMI-1, the SAGE-4λ data set was compiled, which hinges on the measurements at four wavelengths of the SAGE (Stratospheric Aerosol and Gas Experiment) II satellite instrument and uses ground-based Lidar measurements for gap-filling immediately after the Mt. Pinatubo eruption, when the stratosphere was optically opaque for SAGE II. For CMIP6, the new SAGE-3λ data set was compiled, which excludes the least reliable SAGE II wavelength and uses CLAES (Cryogenic Limb Array Etalon Spectrometer) measurements on UARS, the Upper Atmosphere Research Satellite, for gap-filling following the Mt. Pinatubo eruption instead of ground-based Lidars. Here, we performed SOCOLv3 (Solar Climate Ozone Links version 3) chemistry-climate model simulations of the recent past (1986–2005) to investigate the impact of the Mt. Pinatubo eruption in 1991 on stratospheric temperature and ozone and how this response differs depending on which aerosol data set is applied. The use of SAGE-4λ results in heating and ozone loss being overestimated in the lower stratosphere compared to observations in the post-eruption period by approximately 3 K and 0.2 ppmv, respectively. However, less heating occurs in the model simulations based on SAGE-3λ, because the improved gap-filling procedures after the eruption lead to less aerosol loading in the tropical lower stratosphere. As a result, simulated temperature anomalies in the model simulations based on SAGE-3λ for CMIP6 are in excellent agreement with MERRA and ERA-Interim reanalyses in the post-eruption period. Less heating in the simulations with SAGE-3λ means that the rate of tropical upwelling does not strengthen as much as it does in the simulations with SAGE-4λ, which limits dynamical uplift of ozone and therefore provides more time for ozone to accumulate in tropical mid-stratospheric air. Ozone loss following the Mt. Pinatubo eruption is overestimated by 0.1 ppmv in the model simulations based on SAGE-3λ, which is a better agreement with observations than in the simulations based on SAGE-4λ. Overall, the CMIP6 stratospheric aerosol data set, SAGE-3λ, allows SOCOLv3 to more accurately simulate the post-Pinatubo eruption period.


2010 ◽  
Vol 10 (9) ◽  
pp. 4295-4317 ◽  
Author(s):  
D. Wurl ◽  
R. G. Grainger ◽  
A. J. McDonald ◽  
T. Deshler

Abstract. Stratospheric aerosol particles under non-volcanic conditions are typically smaller than 0.1 μm. Due to fundamental limitations of the scattering theory in the Rayleigh limit, these tiny particles are hard to measure by satellite instruments. As a consequence, current estimates of global aerosol properties retrieved from spectral aerosol extinction measurements tend to be strongly biased. Aerosol surface area densities, for instance, are observed to be about 40% smaller than those derived from correlative in situ measurements (Deshler et al., 2003). An accurate knowledge of the global distribution of aerosol properties is, however, essential to better understand and quantify the role they play in atmospheric chemistry, dynamics, radiation and climate. To address this need a new retrieval algorithm was developed, which employs a nonlinear Optimal Estimation (OE) method to iteratively solve for the monomodal size distribution parameters which are statistically most consistent with both the satellite-measured multi-wavelength aerosol extinction data and a priori information. By thus combining spectral extinction measurements (at visible to near infrared wavelengths) with prior knowledge of aerosol properties at background level, even the smallest particles are taken into account which are practically invisible to optical remote sensing instruments. The performance of the OE retrieval algorithm was assessed based on synthetic spectral extinction data generated from both monomodal and small-mode-dominant bimodal sulphuric acid aerosol size distributions. For monomodal background aerosol, the new algorithm was shown to fairly accurately retrieve the particle sizes and associated integrated properties (surface area and volume densities), even in the presence of large extinction uncertainty. The associated retrieved uncertainties are a good estimate of the true errors. In the case of bimodal background aerosol, where the retrieved (monomodal) size distributions naturally differ from the correct bimodal values, the associated surface area (A) and volume densities (V) are, nevertheless, fairly accurately retrieved, except at values larger than 1.0 μm2 cm−3 (A) and 0.05 μm3 cm−3 (V), where they tend to underestimate the true bimodal values. Due to the limited information content in the SAGE II spectral extinction measurements this kind of forward model error cannot be avoided here. Nevertheless, the retrieved uncertainties are a good estimate of the true errors in the retrieved integrated properties, except where the surface area density exceeds the 1.0 μm2 cm−3 threshold. When applied to near-global SAGE II satellite extinction measured in 1999 the retrieved OE surface area and volume densities are observed to be larger by, respectively, 20–50% and 10–40% compared to those estimates obtained by the SAGE~II operational retrieval algorithm. An examination of the OE algorithm biases with in situ data indicates that the new OE aerosol property estimates tend to be more realistic than previous estimates obtained from remotely sensed data through other retrieval techniques. Based on the results of this study we therefore suggest that the new Optimal Estimation retrieval algorithm is able to contribute to an advancement in aerosol research by considerably improving current estimates of aerosol properties in the lower stratosphere under low aerosol loading conditions.


2020 ◽  
Author(s):  
Christoph Brühl ◽  
Hans Schlager ◽  
Ralf Weigel ◽  
Oliver Appel ◽  
Stephan Borrmann ◽  
...  

<p>Results from a transient 28 year simulation with the chemistry climate model EMAC with interactive modal aerosol scheme nudged to observed tropospheric meteorology (ERA-Interim) which includes about 500 volcanic SO<sub>2</sub> injections are compared with in situ aircraft observations in the UT/LS in the Asian Monsoon anticyclone. Enhanced SO<sub>2</sub> observed by STRATOMAS and enhanced sulfate aerosol observed by ERICA in the LS point to impact of several explosive eruptions of the Indonesian volcano Sinabung during summer 2017 seen by the OSIRIS satellite instrument. This is supported by freshly nucleated particles observed by COPAS in the UTLS. We present several sensitivity studies with EMAC with different assumptions on the injection patterns in comparison to the observations in July/August 2017.  <br>The monsoon dynamics distributes the volcanic material together with Asian pollution into the global lower stratosphere.</p>


2017 ◽  
Vol 17 (21) ◽  
pp. 13139-13150 ◽  
Author(s):  
Laura E. Revell ◽  
Andrea Stenke ◽  
Beiping Luo ◽  
Stefanie Kremser ◽  
Eugene Rozanov ◽  
...  

Abstract. To simulate the impacts of volcanic eruptions on the stratosphere, chemistry–climate models that do not include an online aerosol module require temporally and spatially resolved aerosol size parameters for heterogeneous chemistry and aerosol radiative properties as a function of wavelength. For phase 1 of the Chemistry-Climate Model Initiative (CCMI-1) and, later, for phase 6 of the Coupled Model Intercomparison Project (CMIP6) two such stratospheric aerosol data sets were compiled, whose functional capability and representativeness are compared here. For CCMI-1, the SAGE-4λ data set was compiled, which hinges on the measurements at four wavelengths of the SAGE (Stratospheric Aerosol and Gas Experiment) II satellite instrument and uses ground-based lidar measurements for gap-filling immediately after the 1991 Mt Pinatubo eruption, when the stratosphere was too optically opaque for SAGE II. For CMIP6, the new SAGE-3λ data set was compiled, which excludes the least reliable SAGE II wavelength and uses measurements from CLAES (Cryogenic Limb Array Etalon Spectrometer) on UARS, the Upper Atmosphere Research Satellite, for gap-filling following the Mt Pinatubo eruption instead of ground-based lidars. Here, we performed SOCOLv3 (Solar Climate Ozone Links version 3) chemistry–climate model simulations of the recent past (1986–2005) to investigate the impact of the Mt Pinatubo eruption in 1991 on stratospheric temperature and ozone and how this response differs depending on which aerosol data set is applied. The use of SAGE-4λ results in heating and ozone loss being overestimated in the tropical lower stratosphere compared to observations in the post-eruption period by approximately 3 K and 0.2 ppmv, respectively. However, less heating occurs in the model simulations based on SAGE-3λ, because the improved gap-filling procedures after the eruption lead to less aerosol loading in the tropical lower stratosphere. As a result, simulated tropical temperature anomalies in the model simulations based on SAGE-3λ for CMIP6 are in excellent agreement with MERRA and ERA-Interim reanalyses in the post-eruption period. Less heating in the simulations with SAGE-3λ means that the rate of tropical upwelling does not strengthen as much as it does in the simulations with SAGE-4λ, which limits dynamical uplift of ozone and therefore provides more time for ozone to accumulate in tropical mid-stratospheric air. Ozone loss following the Mt Pinatubo eruption is overestimated by up to 0.1 ppmv in the model simulations based on SAGE-3λ, which is a better agreement with observations than in the simulations based on SAGE-4λ. Overall, the CMIP6 stratospheric aerosol data set, SAGE-3λ, allows SOCOLv3 to more accurately simulate the post-Pinatubo eruption period.


2021 ◽  
Author(s):  
Mariam Tidiga ◽  
Gwenaël Berthet ◽  
Fabrice Jegou ◽  
Adriana Bossolasco ◽  
Corinna Kloss ◽  
...  

<p>The cumulative impacts of frequent moderate-magnitude eruptions on stratospheric aerosols were identified among the factors in recent decadal climate trends. Moderate volcanic eruptions are a recurrent source of sulfur dioxide (SO2) in the Upper Troposphere and Lower Stratosphere (UTLS) region and the resulting formation of sulfuric acid aerosol particles from the SO2 emitted provides sites for chemical reactions leading to enhancement of stratospheric optical depth (SAOD) and ozone depletion. Modelling properly the volcanic aerosol content and its evolution in this region is important for radiative impact issues. In this work, we explore the variability of the tropical UTLS aerosol content between 2013 and 2019, a period which was particularly impacted by moderate tropical and mid-latitude volcanic eruptions. For that purpose, space-borne observations from OMPS (version 2, datasets from GES DISC), and IASI, together with simulations by the Whole Atmosphere Community Climate Model (WACCM) coupled with the Community Aerosol and Radiation Model for Atmospheres (CARMA), are used. Different model sensitive experiments, particularly for the injection altitude and timing, have been conducted to evaluate how the model captures the aerosol plume in terms of content, optical and microphysical properties, transport and residence time. We find that the decay of the Calbuco and Kelud plumes observed by OMPS version 2 is well reproduced by the model. Comparisons with unique datasets in the tropical southern hemisphere from the NDACC Maïdo observatory (Reunion Island, France, 20.5°S, 55.5°E) show good agreement between the lidar SAOD observations and WACCM-CARMA SAOD simulations although we observe a difference in the altitude of the maximum aerosol concentration between the model and the in situ profile after Calbuco eruption in April 2015. A particular focus is also made on recent eruptions like Raikoke, Ambae and Ulawun. The plume of the Ambae volcano (15°S, 167°E) which erupted in July 2018 is shown to propagate to the northern hemisphere with some influence until summer 2019 in the Asian monsoon region. For the year 2019, we investigate how the Ulawun (5°S, 151°E; ~0.14 Tg of SO2) tropical eruption and the Raikoke mid-latitude eruption (48°N, 153°E; ~1.5Tg of SO2), have influenced the aerosol burden in the tropics.</p>


2021 ◽  
Vol 164 (3-4) ◽  
Author(s):  
Seshagiri Rao Kolusu ◽  
Christian Siderius ◽  
Martin C. Todd ◽  
Ajay Bhave ◽  
Declan Conway ◽  
...  

AbstractUncertainty in long-term projections of future climate can be substantial and presents a major challenge to climate change adaptation planning. This is especially so for projections of future precipitation in most tropical regions, at the spatial scale of many adaptation decisions in water-related sectors. Attempts have been made to constrain the uncertainty in climate projections, based on the recognised premise that not all of the climate models openly available perform equally well. However, there is no agreed ‘good practice’ on how to weight climate models. Nor is it clear to what extent model weighting can constrain uncertainty in decision-relevant climate quantities. We address this challenge, for climate projection information relevant to ‘high stakes’ investment decisions across the ‘water-energy-food’ sectors, using two case-study river basins in Tanzania and Malawi. We compare future climate risk profiles of simple decision-relevant indicators for water-related sectors, derived using hydrological and water resources models, which are driven by an ensemble of future climate model projections. In generating these ensembles, we implement a range of climate model weighting approaches, based on context-relevant climate model performance metrics and assessment. Our case-specific results show the various model weighting approaches have limited systematic effect on the spread of risk profiles. Sensitivity to climate model weighting is lower than overall uncertainty and is considerably less than the uncertainty resulting from bias correction methodologies. However, some of the more subtle effects on sectoral risk profiles from the more ‘aggressive’ model weighting approaches could be important to investment decisions depending on the decision context. For application, model weighting is justified in principle, but a credible approach should be very carefully designed and rooted in robust understanding of relevant physical processes to formulate appropriate metrics.


2021 ◽  
Vol 13 (9) ◽  
pp. 1846
Author(s):  
Vivek Kumar ◽  
Isabel M. Morris ◽  
Santiago A. Lopez ◽  
Branko Glisic

Estimating variations in material properties over space and time is essential for the purposes of structural health monitoring (SHM), mandated inspection, and insurance of civil infrastructure. Properties such as compressive strength evolve over time and are reflective of the overall condition of the aging infrastructure. Concrete structures pose an additional challenge due to the inherent spatial variability of material properties over large length scales. In recent years, nondestructive approaches such as rebound hammer and ultrasonic velocity have been used to determine the in situ material properties of concrete with a focus on the compressive strength. However, these methods require personnel expertise, careful data collection, and high investment. This paper presents a novel approach using ground penetrating radar (GPR) to estimate the variability of in situ material properties over time and space for assessment of concrete bridges. The results show that attributes (or features) of the GPR data such as raw average amplitudes can be used to identify differences in compressive strength across the deck of a concrete bridge. Attributes such as instantaneous amplitudes and intensity of reflected waves are useful in predicting the material properties such as compressive strength, porosity, and density. For compressive strength, one alternative approach of the Maturity Index (MI) was used to estimate the present values and compare with GPR estimated values. The results show that GPR attributes could be successfully used for identifying spatial and temporal variation of concrete properties. Finally, discussions are presented regarding their suitability and limitations for field applications.


Sign in / Sign up

Export Citation Format

Share Document