Earth Observation Data Group

Abstract:

On 15 January 2022, the Hunga Tonga–Hunga Ha’apai (HTHH) eruption injected 146 MtH2O and 0.42 MtSO2 into the stratosphere. This large water vapour perturbation means that HTHH will probably increase the net radiative forcing, unusual for a large volcanic eruption, increasing the chance of the global surface temperature anomaly temporarily exceeding 1.5 °C over the coming decade. Here we estimate the radiative response to the HTHH eruption and derive the increased risk that the global mean surface temperature anomaly shortly exceeds 1.5 °C following the eruption. We show that HTHH has a tangible impact of the chance of imminent 1.5 °C exceedance (increasing the chance of at least one of the next 5 years exceeding 1.5 °C by 7%), but the level of climate policy ambition, particularly the mitigation of short-lived climate pollutants, dominates the 1.5 °C exceedance outlook over decadal timescales.

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:

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.