scholarly journals The 2019 Raikoke volcanic eruption: Part 1 Dispersion model simulations and satellite retrievals of volcanic sulfur dioxide

2020 ◽  
Author(s):  
Johannes de Leeuw ◽  
Anja Schmidt ◽  
Claire S. Witham ◽  
Nicolas Theys ◽  
Isabelle A. Taylor ◽  
...  

Abstract. Volcanic eruptions can cause significant disruption to society and numerical models are crucial for forecasting the dispersion of erupted material. Here we assess the skill and limitations of the Met Office’s Numerical Atmospheric-dispersion Modelling Environment (NAME) in simulating the dispersion of the sulfur dioxide (SO2) cloud from the 21–22 June 2019 eruption of the Raikoke volcano (48.3° N, 153.2° E). The eruption emitted around 1.5 ± 0.2 Tg of SO2, which represents the largest volcanic emission of SO2 into the stratosphere since the 2011 Nabro eruption. We simulate the temporal evolution of the volcanic SO2 cloud across the Northern Hemisphere (NH) and compare our model simulations to high-resolution SO2 measurements from the Tropospheric Monitoring Instrument (TROPOMI) and the Infrared Atmospheric Sounding Interferometer (IASI) satellite SO2 products. We show that NAME accurately simulates the observed location and horizontal extent of the SO2 cloud during the first 2–3 weeks after the eruption, but is unable, in its standard configuration, to capture the extent and precise location of very high-concentration regions within the volcanic cloud. Using the Fractional Skill Score as metric for model skill, NAME shows skill in simulating the horizontal extent of the cloud for 12–17 days after the eruption where vertical column densities (VCD) of SO2 (in Dobson Units, DU) are above 1 DU. For SO2 VCDs above 20 DU, which are predominantly observed as small-scale features within the SO2 cloud, the model shows skill on the order of 2–4 days only. The lower skill for these high-concentration regions is partly explained by the model-simulated SO2 cloud in NAME being too diffuse compared to TROPOMI retrievals. Reducing the standard diffusion parameters used in NAME by a factor of four results in a slightly increased model skill during the first five days of the simulation, but on longer timescales the simulated SO2 cloud remains too diffuse when compared to TROPOMI measurements. We find that the temporal evolution of the NH-mean SO2 mass burden simulated by NAME strongly depends on the fraction of SO2 mass emitted into the lower stratosphere, which is uncertain for the 2019 Raikoke eruption. When emitting 0.9–1.1 Tg of SO2 into the lower stratosphere (11–18 km) and 0.4–0.7 Tg into the upper troposphere (8–11 km), both NAME and TROPOMI show a similar peak in SO2 mass burden (1.4–1.6 Tg of SO2) with an average SO2 e-folding time of 14–15 days in the NH. Our work demonstrates the large potential of using high-resolution satellite retrievals to identify and rectify limitations in dispersion models like NAME, which will ultimately help to improve dispersion modelling efforts of volcanic SO2 clouds.

2021 ◽  
Vol 21 (14) ◽  
pp. 10851-10879
Author(s):  
Johannes de Leeuw ◽  
Anja Schmidt ◽  
Claire S. Witham ◽  
Nicolas Theys ◽  
Isabelle A. Taylor ◽  
...  

Abstract. Volcanic eruptions can cause significant disruption to society, and numerical models are crucial for forecasting the dispersion of erupted material. Here we assess the skill and limitations of the Met Office's Numerical Atmospheric-dispersion Modelling Environment (NAME) in simulating the dispersion of the sulfur dioxide (SO2) cloud from the 21–22 June 2019 eruption of the Raikoke volcano (48.3∘ N, 153.2∘ E). The eruption emitted around 1.5±0.2 Tg of SO2, which represents the largest volcanic emission of SO2 into the stratosphere since the 2011 Nabro eruption. We simulate the temporal evolution of the volcanic SO2 cloud across the Northern Hemisphere (NH) and compare our model simulations to high-resolution SO2 measurements from the TROPOspheric Monitoring Instrument (TROPOMI) and the Infrared Atmospheric Sounding Interferometer (IASI) satellite SO2 products. We show that NAME accurately simulates the observed location and horizontal extent of the SO2 cloud during the first 2–3 weeks after the eruption but is unable, in its standard configuration, to capture the extent and precise location of the highest magnitude vertical column density (VCD) regions within the observed volcanic cloud. Using the structure–amplitude–location (SAL) score and the fractional skill score (FSS) as metrics for model skill, NAME shows skill in simulating the horizontal extent of the cloud for 12–17 d after the eruption where VCDs of SO2 (in Dobson units, DU) are above 1 DU. For SO2 VCDs above 20 DU, which are predominantly observed as small-scale features within the SO2 cloud, the model shows skill on the order of 2–4 d only. The lower skill for these high-SO2-VCD regions is partly explained by the model-simulated SO2 cloud in NAME being too diffuse compared to TROPOMI retrievals. Reducing the standard horizontal diffusion parameters used in NAME by a factor of 4 results in a slightly increased model skill during the first 5 d of the simulation, but on longer timescales the simulated SO2 cloud remains too diffuse when compared to TROPOMI measurements. The skill of NAME to simulate high SO2 VCDs and the temporal evolution of the NH-mean SO2 mass burden is dominated by the fraction of SO2 mass emitted into the lower stratosphere, which is uncertain for the 2019 Raikoke eruption. When emitting 0.9–1.1 Tg of SO2 into the lower stratosphere (11–18 km) and 0.4–0.7 Tg into the upper troposphere (8–11 km), the NAME simulations show a similar peak in SO2 mass burden to that derived from TROPOMI (1.4–1.6 Tg of SO2) with an average SO2 e-folding time of 14–15 d in the NH. Our work illustrates how the synergy between high-resolution satellite retrievals and dispersion models can identify potential limitations of dispersion models like NAME, which will ultimately help to improve dispersion modelling efforts of volcanic SO2 clouds.


Atmosphere ◽  
2020 ◽  
Vol 11 (9) ◽  
pp. 900 ◽  
Author(s):  
Zhe Wang ◽  
Itsushi Uno ◽  
Kazuo Osada ◽  
Syuichi Itahashi ◽  
Keiya Yumimoto ◽  
...  

Atmospheric ammonia (NH3) plays an important role in the formation of secondary inorganic aerosols, the neutralization of acid rain, and the deposition to ecosystems, but has not been well understood yet, especially over East Asia. Based on the GEOS-Chem model results, the IASI satellite retrievals, the in-site surface observations of a nationwide filter pack (FP) network over Japan and the long-term high resolution online NH3 measurements at Fukuoka of western Japan, the spatio-temporal distributions of atmospheric NH3 over East Asia was analyzed comprehensively. A significant seasonal variation with a summer peak was found in all datasets. Comparison between the satellite retrievals and model simulations indicated that the IASI NH3 vertical column density (VCD) showed good consistency with GEOS-Chem results over North and central China, but had large differences over South China due to the effect of clouds. Over the Japan area, GEOS-Chem simulated NH3 concentrations successfully reproduced the spatio-temporal variations compared with in-situ observations, while IASI NH3 VCD retrievals were below or near the detection limit and difficult to obtain a reasonable correlation for with model results. The comprehensive analysis indicated that there were still some differences among different datasets, and more in-situ observations, improved satellite retrievals, and high-resolution model simulations with more accurate emissions are necessary for better understanding the atmospheric NH3 over East Asia.


2017 ◽  
Vol 2017 ◽  
pp. 1-13 ◽  
Author(s):  
Jiali Luo ◽  
Jiayao Song ◽  
Hongying Tian ◽  
Lei Liu ◽  
Xinlei Liang

We use ERA-Interim reanalysis, MLS observations, and a trajectory model to examine the chemical transport and tracers distribution in the Upper Troposphere and Lower Stratosphere (UTLS) associated with an east-west oscillation case of the anticyclone in 2016. The results show that the spatial distribution of water vapor (H2O) was more consistent with the location of the anticyclone than carbon monoxide (CO) at 100 hPa, and an independent relative high concentration center was only found in H2O field. At 215 hPa, although the anticyclone center also migrated from the Tibetan Mode (TM) to the Iranian Mode (IM), the relative high concentration centers of both tracers were always colocated with regions where upward motion was strong in the UTLS. When the anticyclone migrated from the TM, air within the anticyclone over Tibetan Plateau may transport both westward and eastward but was always within the UTLS. The relative high concentration of tropospheric tracers within the anticyclone in the IM was from the east and transported by the westward propagation of the anticyclone rather than being lifted from surface directly. Air within the relative high geopotential height centers over Western Pacific was partly from the main anticyclone and partly from lower levels.


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.


1980 ◽  
Vol 2 ◽  
Author(s):  
Fernando A. Ponce

ABSTRACTThe structure of the silicon-sapphire interface of CVD silicon on a (1102) sapphire substrate has been studied in crøss section by high resolution transmission electron microscopy. Multibeam images of the interface region have been obtained where both the silicon and sapphire lattices are directly resolved. The interface is observed to be planar and abrupt to the instrument resolution limit of 3 Å. No interfacial phase is evident. Defects are inhomogeneously distributed at the interface: relatively defect-free regions are observed in the silicon layer in addition to regions with high concentration of defects.


2012 ◽  
Vol 9 (6) ◽  
pp. 3521-3566 ◽  
Author(s):  
R. M. A. Caldeira ◽  
X. Couvelard ◽  
E. Casella ◽  
A. Vetrano

Abstract. A high-resolution ocean circulation modelling system forced with a high-resolution numerical wind product was used to study the mesoscale and sub-mesoscale eddy population of the North-Western Mediterranean Sea, contrasting eddy-activity between the Tyrrhenian and Ligurian sub-basins. Numerical solutions reproduced some of the known regional dynamics, namely the occurrence and oceanic implications of Mistral events, the convective cell leeward of the Gulf of Lion, as well as the Balearic frontal system. Calculated transport across the Corsica Channel followed a similar trend, when compared to the transport computed from a moored current meter. The analysis of the results showed that surface eddy activity is mostly confined to the boundary-currents, whereas in the deeper layers most eddies are concentrated on the central-deeper part of the basins. The Liguro-Provençal basin shows a much higher concentration of intermediate and deep-water eddies, when compared to the Tyrrhenian basin. Sub-mesoscale surface eddies tend to merge and migrate vertically onto intermediate waters. Intense eddy activity in the boundary-current surrounding the Liguro-Provençal Gyre, concentrate high-productivity, manifested by higher concentrations of mean sea surface chlorophyll, in the central part of the gyre, defined herein as the Ligurian Productive Pool (LPP). On average, the Tyrrhenian was mostly oligotrophic except for a small productive vortice in the south-eastern (leeward) side of Corsica. The transport in the Tyrrhenian Gyre, and across the basin is one order of magnitude higher than the transport calculated for the Liguro-Provençal basin. A high concentration of eddies in the passage between the Balearic Archipelago and Sardinia suggests retention and longer residence times of nutrient rich water in the "Ligurian pool", compared to a "fast draining" Tyrrhenian basin. Previous studies support the cyclonic gyre circulation generated in the Liguro-Provençal basin but more studies are needed to address the surface and deep mesoscale activity of the Tyrrhenian basin.


2021 ◽  
Author(s):  
Meike Rotermund ◽  
Vera Bense ◽  
Martyn Chipperfield ◽  
Andreas Engel ◽  
Jens-Uwe Grooß ◽  
...  

<p>We report on measurements of total bromine (Br<sup>tot</sup>) in the upper troposphere and lower stratosphere (UTLS) taken from the German High Altitude and LOng range research aircraft (HALO) over the North Atlantic, Norwegian Sea and north-western Europe in September/ October 2017 during the WISE (Wave-driven ISentropic Exchange) research campaign. Br<sup>tot</sup> is calculated from measured total organic bromine (Br<sup>org</sup>) (i.e., the sum of bromine contained in CH<sub>3</sub>Br, the halons and the major very short-lived brominated substances) added to inorganic bromine (Br<sub>y</sub><sup>inorg</sup>), evaluated from measured BrO and photochemical modelling. Combining these data, the weighted mean [Br<sup>tot</sup>] is 19.2 ± 1.2 ppt in the extratropical lower stratosphere (Ex-LS) of the northern hemisphere. The inferred average Br<sup>tot</sup> for the Ex-LS is slightly smaller than expected for the middle stratosphere in 2016 (~19.6 ppt (ranging from 19-20 ppt) as reported by the WMO/UNEP Assessment (2018)). However, it reflects the expected variability in Br<sup>tot</sup> in the Ex-LS due to influxes of shorter lived brominated source and product gases from different regions of entry. A closer look into Br<sup>org</sup> and Br<sub>y</sub><sup>inorg</sup> as well as simultaneously measured transport tracers (CO, N<sub>2</sub>O, ...) and an air mass lag-time tracer (SF<sub>6</sub>), suggests that a filament of air with elevated Br<sup>tot</sup> protruded into the extratropical lowermost stratosphere (Ex-LMS) from 350-385 K and between equivalent latitudes of 55-80˚N (high bromine filament – HBrF). Lagrangian transport modelling shows the multi-pathway contributions to Ex-LMS bromine. According to CLaMS air mass origin simulations, contributions to the HBrF consist of predominantly isentropic transport from the tropical troposphere (also with elevated [Br<sup>tot</sup>] = 21.6 ± 0.7 ppt) as well as a smaller contribution from an exchange across the extratropical tropopause which are mixed into the stratospheric background air. In contrast, the surrounding LS above and below the HBrF has less tropical tropospheric air, but instead additional stratospheric background air. Of the tropical tropospheric air in the HBrF, the majority is from the outflow of the Asian monsoon anticyclone and the adjacent tropical regions, which greatly influences concentrations of trace gases transported into the Ex-LMS in boreal summer and fall. The resulting increase of Br<sup>tot</sup> in the Ex-LMS and its consequences for ozone is investigated through the TOMCAT/SLIMCAT model simulations. However, more extensive monitoring of total stratospheric bromine in more aged air (i.e., in the middle stratosphere) as well as globally and seasonally is required in addition to model simulations to fully understand its impact on Ex-LMS ozone and the radiative forcing of climate.</p>


Sign in / Sign up

Export Citation Format

Share Document