51ÁÔÆæÈë¿Ú

Skip to main content
Department Of Physics text logo
  • Research
    • Our research
    • Our research groups
    • Our research in action
    • Research funding support
    • Summer internships for undergraduates
  • 51ÁÔÆæÈë¿Ú
    • Undergraduates
    • Postgraduates
  • Engage
    • For young people
    • For teachers
    • For the public
    • For alumni
    • For business
  • Support
Menu
Juno Jupiter image

Professor Roy Grainger

Reader in Atmospheric Physics

Research theme

  • Climate physics

Sub department

  • Atmospheric, Oceanic and Planetary Physics

Research groups

  • Earth Observation Data Group
Don.Grainger@physics.ox.ac.uk
Telephone: 01865 (2)72888
Robert Hooke Building, room S47
  • About
  • Publications

Mie scattering from optically levitated mixed sulfuric acid-silica core- shell aerosols: observation of core-shell morphology for atmospheric science

Physical Chemistry Chemical Physics Royal Society of Chemistry 24 (2022) 5813-5822

Abstract:

Sulfuric acid is shown to form a core–shell particle on a micron-sized, optically-trapped spherical silica bead. The refractive indices of the silica and sulfuric acid, along with the shell thickness and bead radius were determined by reproducing Mie scattered optical white light as a function of wavelength in Mie spectroscopy. Micron-sized silica aerosols (silica beads were used as a proxy for atmospheric silica minerals) were levitated in a mist of sulfuric acid particles; continuous collection of Mie spectra throughout the collision of sulfuric acid aerosols with the optically trapped silica aerosol demonstrated that the resulting aerosol particle had a core–shell morphology. Contrastingly, the collision of aqueous sulfuric acid aerosols with optically trapped polystyrene aerosol resulted in a partially coated system. The light scattering from the optically levitated aerosols was successfully modelled to determine the diameter of the core aerosol (±0.003 μm), the shell thickness (±0.0003 μm) and the refractive index (±0.007). The experiment demonstrated that the presence of a thin film rapidly changed the light scattering of the original aerosol. When a 1.964 μm diameter silica aerosol was covered with a film of sulfuric acid 0.287 μm thick, the wavelength dependent Mie peak positions resembled sulfuric acid. Thus mineral aerosol advected into the stratosphere would likely be coated with sulfuric acid, with a core–shell morphology, and its light scattering properties would be effectively indistinguishable from a homogenous sulfuric acid aerosol if the film thickness was greater than a few 100 s of nm for UV-visible wavelengths.

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

Atmospheric Chemistry and Physics Copernicus Publications 21:14 (2021) 10851-10879

Authors:

Johannes de Leeuw, Anja Schmidt, Claire S Witham, Nicolas Theys, Isabelle A Taylor, Roy Grainger, Richard J Pope, Jim Haywood, Martin Osborne, Nina I Kristiansen

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.

The Independent Volcanic Eruption Source Parameter Archive (IVESPA, version 1.0): a new observational database to support explosive eruptive column model validation and development

Journal of Volcanology and Geothermal Research Elsevier 417 (2021) 107295

Authors:

Tj Aubry, S Engwell, C Bonadonna, G Carazzo, S Scollo, Ar Van Eaton, Ia Taylor, D Jessop, J Eychenne, M Gouhier, Lg Mastin, Kl Wallace, S Biass, M Bursik, Roy Grainger, Am Jellinek, A Schmidt

Abstract:

Eruptive column models are powerful tools for investigating the transport of volcanic gas and ash, reconstructing past explosive eruptions, and simulating future hazards. However, the evaluation of these models is challenging as it requires independent estimates of the main model inputs (e.g. mass eruption rate) and outputs (e.g. column height). There exists no database of independently estimated eruption source parameters (ESPs) that is extensive, standardized, maintained, and consensus-based. This paper introduces the Independent Volcanic Eruption Source Parameter Archive (IVESPA, ivespa.co.uk), a community effort endorsed by the International Association of Volcanology and Chemistry of the Earth's Interior (IAVCEI) Commission on Tephra Hazard Modelling. We compiled data for 134 explosive eruptive events, spanning the 1902-2016 period, with independent estimates of: i) total erupted mass of fall deposits; ii) duration; iii) eruption column height; and iv) atmospheric conditions. Crucially, we distinguish plume top versus umbrella spreading height, and the height of ash versus sulphur dioxide injection. All parameter values provided have been vetted independently by at least two experts. Uncertainties are quantified systematically, including flags to describe the degree of interpretation of the literature required for each estimate. IVESPA also includes a range of additional parameters such as total grain size distribution, eruption style, morphology of the plume (weak versus strong), and mass contribution from pyroclastic density currents, where available. We discuss the future developments and potential applications of IVESPA and make recommendations for reporting ESPs to maximize their usability across different applications. IVESPA covers an unprecedented range of ESPs and can therefore be used to evaluate and develop eruptive column models across a wide range of conditions using a standardized dataset.

Observations of plumes from the 2019 Raikoke eruption with the Infrared Atmospheric Sounding Interferometer (IASI)

(2021)

Authors:

Isabelle A Taylor, Roy G Grainger, Tamsin A Mather

Abstract:

<p>Raikoke, a remote volcano in the Kuril Islands, erupted on the 21<sup>st</sup> June 2019. The eruption injected significant quantities of SO<sub>2</sub> into the atmosphere along with volcanic ash. These plumes have been studied with tools developed for the Infrared Atmospheric Sounding Interferometer (IASI) by the Earth Observation Data Group (EODG) at the University of Oxford. IASI is a hyperspectral sensor onboard of three meteorological satellites (Metop A, B and C). Each instrument obtains near global coverage twice a day and has a spectral range which includes sensitivity to both SO<sub>2</sub> and ash: making them useful for studying the Raikoke plumes. A fast linear SO<sub>2</sub> retrieval was first applied to flag pixels with elevated amounts of SO<sub>2</sub>. With this tool it was possible to follow the Raikoke plume as it circulated the northern hemisphere above 30 degrees, with parts of the plume still visible around 2 months after the eruption took place. Next an iterative SO<sub>2</sub> retrieval was used to quantify the amount and height of the SO<sub>2</sub> in each pixel. In the first few days after the eruption took place, very high column amounts are recorded, in some cases exceeding 600 DU. Using this retrieval, a preliminary estimate of 1.6 Tg was obtained for the total amount of SO<sub>2 </sub>emitted (measured on the 23<sup>rd</sup> of June). Height information from this technique shows that there were probably multiple injection heights during the eruption and that SO<sub>2</sub> was emitted into both the troposphere and stratosphere. The tropospheric plume remains visible for just a few days after the eruption, while the stratospheric portion of the plume persists for several weeks.</p>

Extremely fast retrieval of volcanic SO2 layer heights from UV satellite data using inverse learning machines

Copernicus Publications (2021)

Authors:

Pascal Hedelt, MariLiza Koukouli, Konstantinos Michaelidis, Taylor Isabelle, Dimitris Balis, Don Grainger, Dmitry Efremenko, Diego Loyola

Abstract:

Precise knowledge of the location and height of the volcanic sulfur dioxide (SO2) plume is essential for accurate determination of SO2 emitted by volcanic eruptions, however so far not available in operational near-real time UV satellite retrievals. The FP_ILM algorithm (Full-Physics Inverse Learning Machine) enables for the first time to extract the SO2 layer height information in a matter of seconds for current UV satellites and is thus applicable in NRT environments.

The FP_ILM combines a principal component analysis (PCA) and a neural network approach (NN) to extract the information about the volcanic SO2 layer height from high-resolution UV satellite backscatter measurements. So far, UV based SO2 layer height retrieval algorithms were very time-consuming and therefore not suitable for near-real-time applications like aviation control, although the SO2 LH is essential for accurate determination of SO2 emitted by volcanic eruptions.

In this presentation, we will present the latest FP_ILM algorithm improvements and show results of recent volcanic eruptions.

The SO2 layer height product for Sentinel-5p/TROPOMI is developed in the framework of the SO2 Layer Height (S5P+I: SO2 LH) project, which is part of ESA Sentinel-5p+ Innovation project (S5P+I). The S5P+I project aims to develop novel scientific and operational products to exploit the potential of the S5P/TROPOMI capabilities. The S5P+I: SO2 LH project is dedicated to the generation of an SO2 LH product and its extensive verification with collocated ground- and space-born measurements.

Pagination

  • First page First
  • Previous page Prev
  • …
  • Page 3
  • Page 4
  • Page 5
  • Page 6
  • Current page 7
  • Page 8
  • Page 9
  • Page 10
  • Page 11
  • …
  • Next page Next
  • Last page Last

Footer Menu

  • Contact us
  • Giving to the Dept of Physics
  • Work with us
  • Media

User account menu

  • Log in

Follow us

FIND US

Clarendon Laboratory,

Parks Road,

Oxford,

OX1 3PU

CONTACT US

Tel: +44(0)1865272200

Department Of Physics text logo

© University of Oxford - Department of Physics

Cookies | Privacy policy | Accessibility statement

  • Home
  • Research
  • 51ÁÔÆæÈë¿Ú
  • Engage
  • Our people
  • News & Comment
  • Events
  • Our facilities & services
  • 51ÁÔÆæÈë¿Ú
  • Giving to Physics