Professor Roy Grainger

Abstract:

Field-portable Open Path Fourier Transform Infrared (OP-FTIR) spectrometers can be used to remotely measure the composition of volcanic plumes using absorption spectroscopy, providing invaluable data on total gas emissions. Quantifying the temporal evolution of gas compositions during an eruption helps develop models of volcanic processes and aids in eruption forecasting. Absorption measurements require a viewing geometry which aligns infrared source, plume, and instrument, which can be challenging. Here, we present a fast retrieval algorithm to estimate quantities of gas, ash and sulphate aerosols from thermal emission OP-FTIR measurements, and the results from two pilot campaigns on Stromboli volcano in Italy in 2019 and 2021. We validate the method by comparing time series of SO2 slant column densities retrieved using our method with those obtained from a conventional UV spectrometer, demonstrating that the two methods generally agree to within a factor of 2. The algorithm correctly identifies ash-rich plumes and gas bursts associated with explosions and quantifies the mass column densities and particle sizes of ash and sulphate aerosols (SA) in the plume. We compare the ash sizes retrieved using our method with the particle size distribution (PSD) of an ash sample collected during the period of measurements in 2019 by flying a Remotely Piloted Aircraft System into the path of a drifting ash plume and find that both modes of the bimodal PSD (a fine fraction with diameter around 5–10 μm and a coarse fraction around 65 μm) are identified within our datasets at different times. We measure a decrease in the retrieved ash particle size with distance downwind, consistent with settling of larger particles, which we also observed visually. We measure a decrease in the SO2/SA ratio as the plume travels downwind, coupled with an increase in measured SA particle size (range 2–6 μm), suggesting rapid hygroscopic particle growth and/or SO2 oxidation. We propose that infrared emission spectroscopy can be used to examine physical and chemical changes during plume transport and opens the possibility of remote night-time monitoring of volcanic plume emissions. These ground-based analyses may also aid the refinement of satellite-based aerosol retrievals.

Abstract:

A series of explosive eruptions began at La Soufrière, on St Vincent, on 9th April 2021, with the last event occurring on 22nd April. The eruption was observed by the Advanced Baseline Imager (ABI) on the Geostationary Operational Environmental Satellite East (GOES-East) platform. The instrument measures in 16 channels in the visible, near-infrared and infrared, across a large area which includes the Carribean. During the La Soufrière eruption, the ABI instrument made measurements for the full disc every 10-minutes, and for part of the eruption, a mesoscale region was moved over the volcano, providing data every minute. In this study, data from this satellite instrument has been used to construct four images: (a) a true colour image, (b) a false colour image (12.3 - 10.3 µm, 11.2 - 8.4 µm and 10.3 µm channels assigned to red, green and blue respectively), (c) 10.3 - 11.2 µm brightness temperature difference, and (d) 11.2 - 12.3 µm brightness temperature difference. These have been used to study the plumes produced during the April 2021 eruption, allowing the identification of at least 35 explosive events. In this repository, an animation of the created images has been provided for each of the explosive events identified. The animation names are constructed to indicate the volcano (LaSoufriere), the instrument (ABI), the type of data (i.e. full disc or meso), the event number, the start date/time (yyyymmddhhmm) and the end date/time (yyyymmddhhmm). For example: ‘LaSoufriere_ABI_full_event01_202104091250_202104091740.mp4’ Note that the start and end times indicate the start time of measurement period and some adjustment is needed to ascertain the time over the volcano for the full disc (~ +243 seconds). This has been accounted for in the title shown in each frame of the full disc animations. Also note that two animations are provided for event 5 as measurements from the mesoscale region began during this period. We are grateful to the NOAA Big Data Program through which the GOES-16 data are available. These animations accompany a paper: Taylor, I.A., Grainger, R.G., Prata, A.T., Proud, S,R, Mather, T.A., Pyle, D.M.: A satellite chronology of plumes from the April 2021 eruption of La Soufrière, St Vincent, Accepted for publication in Atmospheric, Chemistry and Physics.

Abstract:

Unsieved ash density measurements from 23 unsieved raw ash samples originating from 15 volcanoes, including mass of sample measured with a digital scale and volume measured with a nitrogen gas pycnometer. Data measured in July-August 2022.

Abstract:

Uncertainty-bounded satellite retrievals of volcanic ash cloud properties such as ash cloud-top height, effective radius, optical depth and mass loading are needed for the robust quantitative assessment required to warn aviation of potential hazards. Moreover, there is an imperative to improve quantitative ash cloud estimation due to the planned move towards quantitative ash concentration forecasts by the Volcanic Ash Advisory Centers. Here we apply the Optimal Retrieval of Aerosol and Cloud (ORAC) algorithm to Advanced Himawari Imager (AHI) measurements of the ash clouds produced by the June 2019 Raikoke (Russia) eruption. The ORAC algorithm uses an optimal estimation technique to consolidate a priori information, satellite measurements and associated uncertainties into uncertainty-bounded estimates of the desired state variables. Using ORAC, we demonstrate several improvements in thermal infrared volcanic ash retrievals applied to broadband imagers. These include an improved treatment of measurement noise, accounting for multi-layer cloud scenarios, distinguishing between heights in the troposphere and stratosphere, and the retrieval of a wider range of effective radii sizes than existing techniques by exploiting information from the 10.4 µm channel. Our results indicate that 0.73 ± 0.40 Tg of very fine ash (radius ≤ 15 µm) was injected into the atmosphere during the main eruptive period from 21 June 18:00 UTC to 22 June 10:00 UTC. The total mass of very fine ash decreased from 0.73 to 0.10 Tg over ∼ 48 h, with an e-folding time of 20 h. We estimate a distal fine ash mass fraction of 0.73 % ± 0.62 % based on the total mass of very fine ash retrieved and the ORAC-derived height–time series. Several distinct ash layers were revealed by the ORAC height retrievals. Generally, ash in the troposphere was composed of larger particles than ash present in the stratosphere. We also find that median ash cloud concentrations fall below peak ash concentration safety limits (< 4 mg m−3) 11–16 h after the eruption begins, if typical ash cloud geometric thicknesses are assumed. The ORAC height retrievals for the near-source plume showed good agreement with GOES-17 side-view height data (R=0.84; bias = −0.75 km); however, a larger negative bias was found when comparing ORAC height retrievals for distal ash clouds against Cloud-Aerosol Lidar with Orthogonal Polarisation (CALIOP) measurements (R=0.67; bias = −2.67 km). The dataset generated here provides uncertainties at the pixel level for all retrieved variables and could potentially be used for dispersion model validation or be implemented in data assimilation schemes. Future work should focus on improving ash detection, improving height estimation in the stratosphere and exploring the added benefit of visible channels for retrieving effective radius and optical depth in opaque regions of nascent ash plumes.

Abstract:

Estimates of the anthropogenic effective radiative forcing (ERF) trend have increased by 50% since 2000 (+0.4W/m2/decade in 2000-2009 to +0.6W/m2/decade in 2010-2019), the majority of which is driven by changes in the aerosol ERF trend, due to aerosol emissions reductions. Here we study the extent to which observations of the climate system agree with these ERF assumptions. We use a large ERF ensemble from IPCC’s Sixth Assessment Report (AR6) to attribute the anthropogenic contributions to global mean surface temperature (GMST), top-of-atmosphere radiative flux, and aerosol optical depth observations. The GMST trend has increased from +0.18°C/decade in 2000-2009 to +0.35°C/decade in 2010-2019, coinciding with the anthropogenic warming trend rising from +0.19°C/decade in 2000-2009 to +0.24°C/decade in 2010-2019. This, and observed trends in top-of-atmosphere radiative fluxes and aerosol optical depths support the claim of an aerosol-induced temporary acceleration in the rate of warming. However, all three observation datasets additionally suggest smaller aerosol ERF trend changes are compatible with observations since 2000, since radiative flux and GMST trends are significantly influenced by internal variability over this period. A zero-trend-change aerosol ERF scenario results in a much smaller anthropogenic warming acceleration since 2000, but is poorly represented in AR6’s ERF ensemble. Short-term ERF trends are difficult to verify using observations, so caution is required in predictions or policy judgments that depend on them, such as estimates of current anthropogenic warming trend, and the time remaining to, or the outstanding carbon budget consistent with, 1.5°C warming. Further systematic research focused on quantifying trends and early identification of acceleration or deceleration is required.