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:

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.

Abstract:

Due to the remote location of many volcanoes, there is substantial uncertainty about the timing, amount and vertical distribution of volcanic ash released when they erupt. One approach to determine these properties is to combine prior estimates with satellite retrievals and simulations from atmospheric dispersion models to create posterior emission estimates, constrained by both the observations and the prior estimates, using a technique known as source inversion. However, the results are dependent not only on the accuracy of the prior assumptions, the atmospheric dispersion model and the observations used, but also on the accuracy of the meteorological data used in the dispersion simulations. In this study, we advance the source inversion approach by using an ensemble of meteorological data from the Met Office Global and Regional Ensemble Prediction System to represent the uncertainty in the meteorological data and apply it to the 2019 eruption of Raikoke. Retrievals from the Himawari-8 satellite are combined with NAME dispersion model simulations to create posterior emission estimates. The use of ensemble meteorology provides confidence in the posterior emission estimates and associated dispersion simulations that are used to produce ash forecasts. Prior mean estimates of fine volcanic ash emissions for the Raikoke eruption based on plume height observations are more than 15 times higher than any of the mean posterior ensemble estimates. In addition, the posterior estimates have a different vertical distribution, with 27 %–44 % of ash being emitted into the stratosphere compared to 8 % in the mean prior estimate. This has consequences for the long-range transport of ash, as deposition to the surface from this region of the atmosphere happens over long timescales. The posterior ensemble spread represents uncertainty in the inversion estimate of the ash emissions. For the first 48 h following the eruption, the prior ash column loadings lie outside an estimate of the error associated with a set of independent satellite retrievals, whereas the posterior ensemble column loadings do not. Applying a risk-based methodology to an ensemble of dispersion simulations using the posterior emissions shows that the area deemed to be of the highest risk to aviation, based on the fraction of ensemble members exceeding predefined ash concentration thresholds, is reduced by 49 %. This is compared to estimates using an ensemble of dispersion simulations using the prior emissions with ensemble meteorology. If source inversion had been used following the eruption of Raikoke, it would have had the potential to significantly reduce disruptions to aviation operations. The posterior inversion emission estimates are also sensitive to uncertainty in other eruption source parameters and internal dispersion model parameters. Extending the ensemble inversion methodology to account for uncertainty in these parameters would give a more complete picture of the emission uncertainty, further increasing confidence in these estimates.

Abstract:

Due to the remote location of many volcanoes, there is substantial uncertainty about the timing, amount and vertical distribution of volcanic ash released when they erupt. One approach to determine these properties is to combine prior estimates with satellite retrievals and simulations from atmospheric dispersion models to create posterior emission estimates, constrained by both the observations and the prior estimates, using a technique known as source inversion. However, the results are dependent not only on the accuracy of the prior assumptions, the atmospheric dispersion model and the observations used, but also on the accuracy of the meteorological data used in the dispersion simulations. In this study, we advance the source inversion approach by using an ensemble of meteorological data from the Met Office Global and Regional Ensemble Prediction System to represent the uncertainty in the meteorological data and apply it to the 2019 eruption of Raikoke. Retrievals from the Himawari-8 satellite are combined with NAME dispersion model simulations to create posterior emission estimates. The use of ensemble meteorology provides confidence in the posterior emission estimates and associated dispersion simulations that are used to produce ash forecasts. Prior mean estimates of fine volcanic ash emissions for the Raikoke eruption based on plume height observations are more than 15 times higher than any of the mean posterior ensemble estimates. In addition, the posterior estimates have a different vertical distribution, with 27 %-44 % of ash being emitted into the stratosphere compared to 8 % in the mean prior estimate. This has consequences for the long-range transport of ash, as deposition to the surface from this region of the atmosphere happens over long timescales. The posterior ensemble spread represents uncertainty in the inversion estimate of the ash emissions. For the first 48 h following the eruption, the prior ash column loadings lie outside an estimate of the error associated with a set of independent satellite retrievals, whereas the posterior ensemble column loadings do not. Applying a risk-based methodology to an ensemble of dispersion simulations using the posterior emissions shows that the area deemed to be of the highest risk to aviation, based on the fraction of ensemble members exceeding predefined ash concentration thresholds, is reduced by 49 %. This is compared to estimates using an ensemble of dispersion simulations using the prior emissions with ensemble meteorology. If source inversion had been used following the eruption of Raikoke, it would have had the potential to significantly reduce disruptions to aviation operations. The posterior inversion emission estimates are also sensitive to uncertainty in other eruption source parameters and internal dispersion model parameters. Extending the ensemble inversion methodology to account for uncertainty in these parameters would give a more complete picture of the emission uncertainty, further increasing confidence in these estimates. Copyright: