Model validation and data insertion with FALL3D-8.0: exploiting geostationary satellite retrievals of volcanic ash and SO2

The new version of FALL3D has recently been released with several new features and improvements in model physics, solving algorithms, code accuracy and performance [1]. Among the new features are a data insertion scheme and the ability to simulate volcanic SO2 clouds. The data insertion scheme enables users to initialise model runs from satellite retrievals. This modelling approach is useful for removing uncertainties associated with source term parameters such as the mass flow rate, plume height, source duration and start time. Here we demonstrate and validate the new data insertion scheme in FALL3D-8.0 using geostationary satellite retrievals of volcanic ash and SO2

Atmospheric processes affecting the separation of volcanic ash and SO2 in volcanic eruptions: inferences from the May 2011 Grímsvötn eruption

The separation of volcanic ash and sulfur dioxide (SO2) gas is sometimes observed during volcanic eruptions. The exact conditions under which separation occurs are not fully understood but the phenomenon is of importance because of the effects volcanic emissions have on aviation, on the environment, and on the earth’s radiation balance. The eruption of Grímsvötn, a subglacial volcano under the Vatnajökull glacier in Iceland during 21–28 May 2011 produced one of the most spectacular examples of ash and SO2 separation, which led to errors in the forecasting of ash in the atmosphere over northern Europe. Satellite data from several sources coupled with meteorological wind data and photographic evidence suggest that the eruption column was unable to sustain itself, resulting in a large deposition of ash, which left a low-level ash-rich atmospheric plume moving southwards and then eastwards towards the southern Scandinavian coast and a high-level predominantly SO2 plume travelling northwards and then spreading eastwards and westwards. Here we provide observational and modelling perspectives on the separation of ash and SO2 and present quantitative estimates of the masses of ash and SO2 that erupted, the directions of transport, and the likely impacts. We hypothesise that a partial column collapse or “sloughing” fed with ash from pyroclastic density currents (PDCs) occurred during the early stage of the eruption, leading to an ash-laden gravity intrusion that was swept southwards, separated from the main column. Our model suggests that water-mediated aggregation caused enhanced ash removal because of the plentiful supply of source water from melted glacial ice and from entrained atmospheric water. The analysis also suggests that ash and SO2 should be treated with separate source terms, leading to improvements in forecasting the movement of both types of emissions

A novel and inexpensive method for measuring volcanic plume water fluxes at high temporal resolution

© 2017 by the authors.Water vapour (H2O) is the dominant species in volcanic gas plumes. Therefore,measurements of H2O fluxes could provide valuable constraints on subsurface degassing and magmatic processes. However, due to the large and variable concentration of this species in the background atmosphere, little attention has been devoted to monitoring the emission rates of this species from volcanoes. Instead, the focus has been placed on remote measurements of SO2, which is present in far lower abundances in plumes, and therefore provides poorer single flux proxies for overall degassing conditions. Here, we present a new technique for the measurement of H2O emissions at degassing volcanoes at high temporal resolution (≈1 Hz), via remote sensing with low cost digital cameras. This approach is analogous to the use of dual band ultraviolet (UV) cameras for measurements of volcanic SO2 release, but is focused on near infrared absorption by H2O. We report on the field deployment of these devices on La Fossa crater, Vulcano Island, and the North East Crater of Mt. Etna, during which in-plume calibration was performed using a humidity sensor, resulting in estimated mean H2O fluxes of ≈15 kg·s-1 and ≈34 kg·s-1, respectively, in accordance with previously reported literature values. By combining the Etna data with parallel UV camera and Multi-GAS observations, we also derived, for the first time, a combined record of 1 Hz gas fluxes for the three most abundant volcanic gas species: H2O, CO2, and SO2. Spectral analysis of the Etna data revealed oscillations in the passive emissions of all three species, with periods spanning ≈40-175 s, and a strong degree of correlation between the periodicity manifested in the SO2 and H2O data, potentially related to the similar exsolution depths of these two gases. In contrast, there was a poorer linkage between oscillations in these species and those of CO2, possibly due to the deeper exsolution of carbon dioxide, giving rise to distinct periodic degassing behaviour

A satellite chronology of plumes from the April 2021 eruption of La Soufrière, St Vincent

Satellite instruments play a valuable role in detecting, monitoring and characterising emissions of ash and gas into the atmosphere during volcanic eruptions. This study uses two satellite instruments, the Infrared Atmospheric Sounding Interferometer (IASI) and the Advanced Baseline Imager (ABI), to examine the plumes of ash and sulfur dioxide (SO2) from the April 2021 eruption of La Soufrière, St Vincent. The frequent ABI data have been used to construct a 14 d chronology of a series of explosive events at La Soufrière, which is then complemented by measurements of SO2 from IASI, which is able to track the plume as it is transported around the globe. A minimum of 35 eruptive events were identified using true, false and brightness temperature difference maps produced with the ABI data. The high temporal resolution images were used to identify the approximate start and end times, as well as the duration and characteristics of each event. From this analysis, four distinct phases within the 14 d eruption have been defined, each consisting of multiple explosive events with similar characteristics: (1) an initial explosive event, (2) a sustained event lasting over 9 h, (3) a pulsatory phase with 25 explosive events in a 65.3 h period and (4) a waning sequence of explosive events. It is likely that the multiple explosive events during the April 2021 eruption contributed to the highly complex plume structure that can be seen in the IASI measurements of the SO2 column amounts and heights. The bulk of the SO2 from the first three phases of the eruption was transported eastwards, which based on the wind direction at the volcano implies that the SO2 was largely in the upper troposphere. Some of the SO2 was carried to the south and west of the volcano, suggesting a smaller emission of the gas into the stratosphere, there being a shift in wind direction around the height of the tropopause. The retrieved SO2 heights show that the plume had multiple layers but was largely concentrated between 13 and 19 km, with the majority of the SO2 being located in the upper troposphere and around the height of the tropopause, with some emission into the stratosphere. An average e-folding time of 6.07±4.74 d was computed based on the IASI SO2 results: similar to other tropical eruptions of this magnitude and height. The SO2 was trackable for several weeks after the eruption and is shown to have circulated the globe, with parts of it reaching as far as 45∘ S and 45∘ N. Using the IASI SO2 measurements, a time series of the total SO2 mass loading was produced, with this peaking on 13 April (descending orbits) at 0.31±0.09 Tg. Converting these mass values to a temporally varying SO2 flux demonstrated that the greatest emission occurred on 10 April with that measurement incorporating SO2 from the second phase of the eruption (sustained emission) and the beginning of the pulsatory phase. The SO2 flux is then shown to fall during the later stages of the eruption: suggesting a reduction in eruptive energy, something also reflected in ash height estimates obtained with the ABI instrument. A total SO2 emission of 0.63±0.5 Tg of SO2 has been derived, although due to limitations associated with the retrieval, particularly in the first few days after the eruption began, this, the retrieved column amounts and the total SO2 mass on each day should be considered minimum estimates. There are a number of similarities between the 1979 and 2021 eruptions at La Soufrière, with both eruptions consisting of a series of explosive events with varied heights and including some emission into the stratosphere. These similarities highlight the importance of in-depth investigations into eruptions and the valuable contribution of satellite data for this purpose; as these studies aid in learning about a volcano's behaviour, which may allow for better preparation for future eruptive activity

Geostationary aerosol retrievals of extreme biomass burning plumes during the 2019–2020 Australian bushfires

Extreme biomass burning (BB) events, such as those seen during the 2019-2020 Australian bushfire season, are becoming more frequent and intense with climate change. Ground-based observations of these events can provide useful information on the macro-and micro-physical properties of the plumes, but these observations are sparse, especially in regions which are at risk of intense bushfire events. Satellite observations of extreme BB events provide a unique perspective, with the newest generation of geostationary imagers, such as the Advanced Himawari Imager (AHI), observing entire continents at moderate spatial and high temporal resolution. However, current passive satellite retrieval methods struggle to capture the high values of aerosol optical thickness (AOT) seen during these BB events. Accurate retrievals are necessary for global and regional studies of shortwave radiation, air quality modelling and numerical weather prediction. To address these issues, the Optimal Retrieval of Aerosol and Cloud (ORAC) algorithm has used AHI data to measure extreme BB plumes from the 2019-2020 Australian bushfire season. The sensitivity of the retrieval to the assumed optical properties of BB plumes is explored by comparing retrieved AOT with AErosol RObotic NETwork (AERONET) level-1.5 data over the AERONET site at Tumbarumba, New South Wales, between 1 December 2019 at 00:00UTC and 3 January 2020 at 00:00UTC. The study shows that for AOT values >2, the sensitivity to the assumed optical properties is substantial. The ORAC retrievals and AERONET data are compared against the Japan Aerospace Exploration Agency (JAXA) Aerosol Retrieval Product (ARP), Moderate Resolution Imaging Spectroradiometer (MODIS) Deep Blue over land, MODIS MAIAC, Sentinel-3 SYN and VIIRS Deep Blue products. The comparison shows the ORAC retrieval significantly improves coverage of optically thick plumes relative to the JAXA ARP, with approximately twice as many pixels retrieved and peak retrieved AOT values 1.4 times higher than the JAXA ARP. The ORAC retrievals have accuracy scores of 0.742-0.744 compared to the values of 0.718-0.833 for the polar-orbiting satellite products, despite successfully retrieving approximately 28 times as many pixels over the study period as the most successful polar-orbiting satellite product. The AHI and MODIS satellite products are compared for three case studies covering a range of BB plumes over Australia. The results show good agreement between all products for plumes with AOT values ≤2. For extreme BB plumes, the ORAC retrieval finds values of AOT >15, significantly higher than those seen in events classified as extreme by previous studies, although with high uncertainty. A combination of hard limits in the retrieval algorithms and misclassification of BB plumes as cloud prevents the JAXA and MODIS products from returning AOT values significantly greater than 5

Lidar ratios of stratospheric volcanic ash and sulfate aerosols retrieved from CALIOP measurements

We apply a two-way transmittance constraint to nighttime CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) observations of volcanic aerosol layers to retrieve estimates of the particulate lidar ratio (Sp) at 532 nm. This technique is applied to three volcanic eruption case studies that were found to have injected aerosols directly into the stratosphere. Numerous lidar observations permitted characterization of the optical and geometric properties of the volcanic aerosol layers over a time period of 1–2 weeks. For the volcanic ash-rich layers produced by the Puyehue-Cordón Caulle eruption (June 2011), we obtain mean and median particulate lidar ratios of 69 ± 13 sr and 67 sr, respectively. For the sulfate-rich aerosol layers produced by Kasatochi (August 2008) and Sarychev Peak (June 2009), the means of the retrieved lidar ratios were 66 ± 19 sr (median 60 sr) and 63 ± 14 sr (median 59 sr), respectively. The 532 nm layer-integrated particulate depolarization ratios (δp) observed for the Puyehue layers (δp = 0.33 ± 0.03) were much larger than those found for the volcanic aerosol layers produced by the Kasatochi (δp = 0.09 ± 0.03) and Sarychev (δp = 0.05 ± 0.04) eruptions. However, for the Sarychev layers we observe an exponential decay (e-folding time of 3.6 days) in δp with time from 0.27 to 0.03. Similar decreases in the layer-integrated attenuated colour ratios with time were observed for the Sarychev case. In general, the Puyehue layers exhibited larger colour ratios (χ′ = 0.53 ± 0.07) than what was observed for the Kasatochi (χ′ = 0.35 ± 0.07) and Sarychev (χ′ = 0.32 ± 0.07) layers, indicating that the Puyehue layers were generally composed of larger particles. These observations are particularly relevant to the new stratospheric aerosol subtyping classification scheme, which has been incorporated into version 4 of the level 2 CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) data products

The 2019 Raikoke eruption as a testbed used by the Volcano Response group for rapid assessment of volcanic atmospheric impacts

The 21 June 2019 Raikoke eruption (48° N, 153° E) generated one of the largest amounts of sulfur emission to the stratosphere since the 1991 Mt. Pinatubo eruption. Satellite measurements indicate a consensus best estimate of 1.5 Tg for the sulfur dioxide (SO2) injected at an altitude of around 14–15 km. The peak Northern Hemisphere (NH) mean 525 nm stratospheric aerosol optical depth (SAOD) increased to 0.025, a factor of 3 higher than background levels. The Volcano Response (VolRes) initiative provided a platform for the community to share information about this eruption which significantly enhanced coordination efforts in the days after the eruption. A multi-platform satellite observation subgroup formed to prepare an initial report to present eruption parameters including SO2 emissions and their vertical distribution for the modeling community. It allowed us to make the first estimate of what would be the peak in SAOD 1 week after the eruption using a simple volcanic aerosol model. In this retrospective analysis, we show that revised volcanic SO2 injection profiles yield a higher peak injection of the SO2 mass. This highlights difficulties in accurately representing the vertical distribution for moderate SO2 explosive eruptions in the lowermost stratosphere due to limited vertical sensitivity of the current satellite sensors (±2 km accuracy) and low horizontal resolution of lidar observations. We also show that the SO2 lifetime initially assumed in the simple aerosol model was overestimated by 66 %, pointing to challenges for simple models to capture how the life cycle of volcanic gases and aerosols depends on the SO2 injection magnitude, latitude, and height. Using a revised injection profile, modeling results indicate a peak NH monthly mean SAOD at 525 nm of 0.024, in excellent agreement with observations, associated with a global monthly mean radiative forcing of −0.17 W m−2 resulting in an annual global mean surface temperature anomaly of −0.028 K. Given the relatively small magnitude of the forcing, it is unlikely that the surface response can be dissociated from surface temperature variability

The 2019 Raikoke eruption as a testbed used by the Volcano Response group for rapid assessment of volcanic atmospheric impacts

The 2019 Raikoke eruption (Kamchatka, Russia) generated one of the largest emissions of particles and gases into the stratosphere since the 1991 Mt. Pinatubo eruption. The Volcano Response (VolRes) initiative, an international effort, provided a platform for the community to share information about this eruption and assess its climate impact. The eruption led to a minor global surface cooling of 0.02 C in 2020 which is negligible relative to warming induced by human greenhouse gas emissions

Treetop: A Shiny-based application and R package for extracting forest information from LiDAR data for ecologists and conservationists

Individual tree detection (ITD) and crown delineation are two of the most relevant methods for extracting detailed and reliable forest information from LiDAR (Light Detection and Ranging) datasets. However, advanced computational skills and specialized knowledge have been normally required to extract forest information from LiDAR.The development of accessible tools for 3D forest characterization can facilitate rapid assessment by stakeholders lacking a remote sensing background, thus fostering the practical use of LiDAR datasets in forest ecology and conservation. This paper introduces the treetop application, an open-source web-based and R package LiDAR analysis tool for extracting forest structural information at the tree level, including cutting-edge analyses of properties related to forest ecology and management.We provide case studies of how treetop can be used for different ecological applications, within various forest ecosystems. Specifically, treetop was employed to assess post-hurricane disturbance in natural temperate forests, forest homogeneity in industrial forest plantations and the spatial distribution of individual trees in a tropical forest.treetop simplifies the extraction of relevant forest information for forest ecologists and conservationists who may use the tool to easily visualize tree positions and sizes, conduct complex analyses and download results including individual tree lists and figures summarizing forest structural properties. Through this open-source approach, treetop can foster the practical use of LiDAR data among forest conservation and management stakeholders and help ecological researchers to further understand the relationships between forest structure and function.The authors thank Nicholas L. Crookston for co‐developing the web‐LiDAR treetop tool, and the two anonymous reviewers for their helpful suggestions on the first version of the manuscript. This study is based on the work supported by the Department of Defence Strategic Environmental Research and Development Program (SERDP) under grants No. RC‐2243, RC19‐1064 and RC20‐1346 and USDA Forest Service (grand No. PRO00031122