scholarly journals Using Satellite Data to Determine Empirical Relationships between Volcanic Ash Source Parameters

Atmosphere ◽  
2020 ◽  
Vol 11 (4) ◽  
pp. 342 ◽  
Author(s):  
Meelis J. Zidikheri ◽  
Chris Lucas
Keyword(s):  
Case Studies ◽  
Volcanic Ash ◽  
Dispersion Model ◽  
Source Parameters ◽  
Difficult Problem ◽  
Emission Rates ◽  
Empirical Relationships ◽  
Mass Load ◽  
Satellite Retrievals ◽  
Mass Loads

Poor knowledge of dispersion model source parameters related to quantities such as the total fine ash mass emission rate, its effective spatial distribution, and particle size distribution makes the provision of quantitative forecasts of volcanic ash a difficult problem. To ameliorate this problem, we make use of satellite-retrieved mass load data from 14 eruption case studies to estimate fine ash mass emission rates and other source parameters by an inverse modelling procedure, which requires multidimensional sampling of several thousand trial simulations with different values of source parameters. We then estimate the dependence of these optimal source parameters on eruption height. We show that using these empirical relationships in a data assimilation procedure leads to substantial improvements to the forecasts of ash mass loads, with the use of empirical relationships between parameters and eruption height having the added advantage of computational efficiency because of dimensional reduction. In addition, the use of empirical relationships, which encode information in satellite retrievals from past case studies, implies that quantitative forecasts can still be issued even when satellite retrievals of mass load are not available in real time due to cloud cover or other reasons, making it especially useful for operations in the tropics where ice and water clouds are ubiquitous.

Ensemble-based volcanic ash forecasts using satellite retrievals for quantitative verification

2021 ◽  
Author(s):  
Antonio Capponi ◽  
Natalie J. Harvey ◽  
Helen F. Dacre ◽  
Keith Beven ◽  
Mike R. James
Keyword(s):  
Volcanic Ash ◽  
Dispersion Model ◽  
Source Parameters ◽  
Forecast Accuracy ◽  
Satellite Observations ◽  
Ensemble Simulations ◽  
Planning Decisions ◽  
Satellite Retrievals ◽  
Initial Ensemble ◽  
Quantitative Verification

<p>Volcanic ash poses a significant hazard for aviation. If an ash cloud forms as result of an eruption, it forces a series of flight planning decisions that consider important safety and economic factors. These decisions are made using a combination of satellite retrievals and volcanic ash forecasts issued by Volcanic Ash Advisory Centres.  However, forecasts of ash hazard remain deterministic, and lack quantification of the uncertainty that arises from the estimation of eruption source parameters, meteorology and uncertainties within the dispersion model used to perform the simulations. Quantification of these uncertainties is fundamental and could be achieved by using ensemble simulations. Here, we explore how ensemble-based forecasts — performed using the Met Office dispersion model NAME — together with sequential satellite retrievals of ash column loading, may improve forecast accuracy and uncertainty characterization.</p><p>We have developed a new methodology to evaluate each member of the ensemble based on its agreement with the satellite retrievals available at the time. An initial ensemble is passed through a filter of verification metrics and compared with the first available set of satellite observations. Members far from the observations are rejected. The members within a limit of acceptability are used to resample the parameters used in the initial ensemble, and design a new ensemble to compare with the next available set of satellite observations. The filtering process and parameter resampling are applied whenever new satellite observations are available, to create new ensembles propagating forward in time, until all available observations are covered.</p><p>Although the method requires the run of many ensemble batches, and it is not yet suited for operational use, it shows how combining ensemble simulations and sequential satellite retrievals can be used to quantify confidence in ash forecasts. We demonstrate the method by applying it to the recent Raikoke (Kurii Islands, Russia) eruption, which occurred on the 22<sup>nd</sup> July 2019. Each ensemble consists of 1000 members and it is evaluated against 6-hourly HIMAWARI satellite ash retrievals.</p>

scholarly journals Refining an ensemble of volcanic ash forecasts using satellite retrievals: Raikoke 2019

2021 ◽  
Author(s):  
Antonio Capponi ◽  
Natalie J. Harvey ◽  
Helen F. Dacre ◽  
Keith Beven ◽  
Cameron Saint ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Volcanic Ash ◽  
Dispersion Model ◽  
Source Parameters ◽  
Aviation Industry ◽  
Model Parameters ◽  
Flight Operations ◽  
Operational Forecasts ◽  
The Pacific ◽  
Satellite Retrievals

Abstract. Volcanic ash advisories are produced by specialised forecasters who combine several sources of observational data and volcanic ash dispersion model outputs based on their subjective expertise. These advisories are used by the aviation industry to make decisions about where it is safe to fly. However, both observations and dispersion model simulations are subject to various sources of uncertainties that are not represented in operational forecasts. Quantification and communication of these uncertainties are fundamental for making more informed decisions. Here, we develop a data assimilation technique which combines satellite retrievals and volcanic ash transport and dispersion model (VATDM) output, considering uncertainties in both data sources. The methodology is applied to a case study of the 2019 Raikoke eruption. To represent uncertainty in the VATDM output, 1000 simulations are performed by simultaneously perturbing the eruption source parameters, meteorology and internal model parameters (known as the prior ensemble). The ensemble members are filtered, based on their level of agreement with Himawari satellite retrievals of ash column loading, to produce a posterior ensemble that is constrained by the satellite data and its uncertainty. For the Raikoke eruption, filtering the ensemble skews the values of mass eruption rate towards the lower values within the wider parameters ranges initially used in the prior ensemble (mean reduces from 1 Tg h−1 to 0.1 Tg h−1). Furthermore, including satellite observations from subsequent times increasingly constrains the posterior ensemble. These results suggest that the prior ensemble leads to an overestimate of both the magnitude and uncertainty in ash column loadings. Based on the prior ensemble, flight operations would have been severely disrupted over the Pacific Ocean. Using the constrained posterior ensemble, the regions where the risk is overestimated are reduced potentially resulting in fewer flight disruptions. The data assimilation methodology developed in this paper is easily generalisable to other short duration eruptions and to other VATDMs and retrievals of ash from other satellites.

scholarly journals Atmospheric Dispersion Modelling at the London VAAC: A Review of Developments since the 2010 Eyjafjallajökull Volcano Ash Cloud

Atmosphere ◽  
2020 ◽  
Vol 11 (4) ◽  
pp. 352 ◽  
Author(s):  
Frances M. Beckett ◽  
Claire S. Witham ◽  
Susan J. Leadbetter ◽  
Ric Crocker ◽  
Helen N. Webster ◽  
...  
Keyword(s):  
Volcanic Ash ◽  
Dispersion Model ◽  
Source Parameters ◽  
Atmospheric Dispersion ◽  
Lessons Learned ◽  
Dispersion Modelling ◽  
Atmospheric Dispersion Modelling ◽  
Volcanic Ash Cloud ◽  
Atmospheric Dispersion Model ◽  
Ash Cloud

It has been 10 years since the ash cloud from the eruption of Eyjafjallajökull caused unprecedented disruption to air traffic across Europe. During this event, the London Volcanic Ash Advisory Centre (VAAC) provided advice and guidance on the expected location of volcanic ash in the atmosphere using observations and the atmospheric dispersion model NAME (Numerical Atmospheric-Dispersion Modelling Environment). Rapid changes in regulatory response and procedures during the eruption introduced the requirement to also provide forecasts of ash concentrations, representing a step-change in the level of interrogation of the dispersion model output. Although disruptive, the longevity of the event afforded the scientific community the opportunity to observe and extensively study the transport and dispersion of a volcanic ash cloud. We present the development of the NAME atmospheric dispersion model and modifications to its application in the London VAAC forecasting system since 2010, based on the lessons learned. Our ability to represent both the vertical and horizontal transport of ash in the atmosphere and its removal have been improved through the introduction of new schemes to represent the sedimentation and wet deposition of volcanic ash, and updated schemes to represent deep moist atmospheric convection and parametrizations for plume spread due to unresolved mesoscale motions. A good simulation of the transport and dispersion of a volcanic ash cloud requires an accurate representation of the source and we have introduced more sophisticated approaches to representing the eruption source parameters, and their uncertainties, used to initialize NAME. Finally, upper air wind field data used by the dispersion model is now more accurate than it was in 2010. These developments have resulted in a more robust modelling system at the London VAAC, ready to provide forecasts and guidance during the next volcanic ash event.

scholarly journals FALL3D-8.0: a computational model for atmospheric transport and deposition of particles, aerosols and radionuclides – Part 2: Model validation

2021 ◽  
Vol 14 (1) ◽  
pp. 409-436
Author(s):  
Andrew T. Prata ◽  
Leonardo Mingari ◽  
Arnau Folch ◽  
Giovanni Macedonio ◽  
Antonio Costa
Keyword(s):  
Case Studies ◽  
Model Validation ◽  
Volcanic Ash ◽  
Atmospheric Transport ◽  
Test Suite ◽  
Structure Amplitude ◽  
Satellite Retrievals ◽  
Mt Etna ◽  
Ash Cloud ◽  
Tephra Fallout

Abstract. This paper presents model validation results for the latest version release of the FALL3D atmospheric transport model. The code has been redesigned from scratch to incorporate different categories of species and to overcome legacy issues that precluded its preparation towards extreme-scale computing. The model validation is based on the new FALL3D-8.0 test suite, which comprises a set of four real case studies that encapsulate the major features of the model; namely, the simulation of long-range fine volcanic ash dispersal, volcanic SO2 dispersal, tephra fallout deposits and the dispersal and deposition of radionuclides. The first two test suite cases (i.e. the June 2011 Puyehue-Cordón Caulle ash cloud and the June 2019 Raikoke SO2 cloud) are validated against geostationary satellite retrievals and demonstrate the new FALL3D data insertion scheme. The metrics used to validate the volcanic ash and SO2 simulations are the structure, amplitude and location (SAL) metric and the figure of merit in space (FMS). The other two test suite cases (i.e. the February 2013 Mt. Etna ash cloud and associated tephra fallout deposit, and the dispersal of radionuclides resulting from the 1986 Chernobyl nuclear accident) are validated with scattered ground-based observations of deposit load and local particle grain size distributions and with measurements from the Radioactivity Environmental Monitoring database. For validation of tephra deposit loads and radionuclides, we use two variants of the normalised root-mean-square error metric. We find that FALL3D-8.0 simulations initialised with data insertion consistently improve agreement with satellite retrievals at all lead times up to 48 h for both volcanic ash and SO2 simulations. In general, SAL scores lower than 1.5 and FMS scores greater than 0.40 indicate acceptable agreement with satellite retrievals of volcanic ash and SO2. In addition, we show very good agreement, across several orders of magnitude, between the model and observations for the 2013 Mt. Etna and 1986 Chernobyl case studies. Our results, along with the validation datasets provided in the publicly available test suite, form the basis for future improvements to FALL3D (version 8 or later) and also allow for model intercomparison studies.

scholarly journals Improving volcanic ash predictions with the HYSPLIT dispersion model by assimilating MODIS satellite retrievals

2016 ◽  
Author(s):  
Tianfeng Chai ◽  
Alice Crawford ◽  
Barbara Stunder ◽  
Michael Pavolonis ◽  
Roland Draxler ◽  
...  
Keyword(s):  
Volcanic Ash ◽  
Dispersion Model ◽  
Temporal Variations ◽  
Current Case ◽  
Cost Functional ◽  
Free Region ◽  
Model Predictions ◽  
Satellite Retrievals ◽  
Moderate Resolution Imaging Spectroradiometer ◽  
Ash Cloud

Abstract. Currently NOAA's National Weather Service (NWS) runs the HYSPLIT dispersion model with a unit mass release rate to predict the transport and dispersion of volcanic ash. The model predictions provide information for the Volcanic Ash Advisory Centers (VAAC) to issue advisories to meteorological watch offices, area control centers, flight information centers, and others. This research aims provide quantitative forecasts of ash distributions generated by objectively and optimally estimating the volcanic ash source strengths, vertical distribution and temporal variations using an observation-modeling inversion technique. In this top-down approach, a cost functional is defined to mainly quantify the differences between model predictions and the satellite measurements of column integrated ash concentrations, weighted by the model and observation uncertainties. Minimizing this cost functional by adjusting the sources provides the volcanic ash emission estimates. As an example, MODIS (MOderate Resolution Imaging Spectroradiometer) satellite retrievals of the 2008 Kasatochi volcanic ash clouds are used to test the HYSPLIT volcanic ash inverse system. Because the satellite retrievals include the ash cloud top height but not the bottom height, three options for matching the model concentrations to the observed mass loadings are tested. Although the emission estimates vary significantly with different options the subsequent model predictions with the different release estimates all show decent skill when evaluated against the unassimilated satellite observations at later times. Among the three options, integrating over three model layers yields slightly better results than integrating from the surface up to the volcanic ash cloud top or using a single model layer. Inverse tests also show that including the ash-free region to constrain the model is not beneficial for the current case. In addition, extra constraints to the source terms can be given by explicitly enforcing ``no-ash'' for the atmosphere columns above or below the observed ash cloud top height. However, in this case such extra constraints are not helpful for the inverse modeling. It is also found that simultaneously assimilating observations at different times produces better hindcasts than only assimilating the most recent observations.

Estimation of optimal dispersion model source parameters using satellite detections of volcanic ash

2017 ◽  
Vol 122 (15) ◽  
pp. 8207-8232 ◽  
Author(s):  
Meelis J. Zidikheri ◽  
Christopher Lucas ◽  
Rodney J. Potts
Keyword(s):  
Volcanic Ash ◽  
Dispersion Model ◽  
Source Parameters ◽  
Model Source

scholarly journals Improving Ensemble Volcanic Ash Forecasts by Direct Insertion of Satellite Data and Ensemble Filtering

Atmosphere ◽  
2021 ◽  
Vol 12 (9) ◽  
pp. 1215
Author(s):  
Meelis J. Zidikheri ◽  
Chris Lucas
Keyword(s):  
Volcanic Ash ◽  
Weather Prediction ◽  
Source Parameters ◽  
Volcanic Eruptions ◽  
Forecast Skill ◽  
Aviation Industry ◽  
Transport Modelling ◽  
Volcanic Source ◽  
Mass Load ◽  
Ensemble Filtering

Improved quantitative forecasts of volcanic ash are in great demand by the aviation industry to enable better risk management during disruptive volcanic eruption events. However, poor knowledge of volcanic source parameters and other dispersion and transport modelling uncertainties, such as those due to errors in numerical weather prediction fields, make this problem very challenging. Nonetheless, satellite-based algorithms that retrieve ash properties, such as mass load, effective radius, and cloud top height, combined with inverse modelling techniques, such as ensemble filtering, can significantly ameliorate these problems. The satellite-retrieved data can be used to better constrain the volcanic source parameters, but they can also be used to avoid the description of the volcanic source altogether by direct insertion into the forecasting model. In this study we investigate the utility of the direct insertion approach when employed within an ensemble filtering framework. Ensemble members are formed by initializing dispersion models with data from different timesteps, different values of cloud top height, thickness, and NWP ensemble members. This large ensemble is then filtered with respect to observations to produce a refined forecast. We apply this approach to 14 different eruption case studies in the tropical atmosphere. We demonstrate that the direct insertion of data improves model forecast skill, particularly when it is used in a hybrid ensemble in which some ensemble members are initialized from the volcanic source. Moreover, good forecast skill can be obtained even when detailed satellite retrievals are not available, which is frequently the case for volcanic eruptions in the tropics.

scholarly journals Improving volcanic ash predictions with the HYSPLIT dispersion model by assimilating MODIS satellite retrievals

2017 ◽  
Vol 17 (4) ◽  
pp. 2865-2879 ◽  
Author(s):  
Tianfeng Chai ◽  
Alice Crawford ◽  
Barbara Stunder ◽  
Michael J. Pavolonis ◽  
Roland Draxler ◽  
...  
Keyword(s):  
Volcanic Ash ◽  
Dispersion Model ◽  
Temporal Variations ◽  
Current Case ◽  
Cost Functional ◽  
Free Region ◽  
Model Predictions ◽  
Satellite Retrievals ◽  
Moderate Resolution Imaging Spectroradiometer ◽  
Ash Cloud

Abstract. Currently, the National Oceanic and Atmospheric Administration (NOAA) National Weather Service (NWS) runs the HYSPLIT dispersion model with a unit mass release rate to predict the transport and dispersion of volcanic ash. The model predictions provide information for the Volcanic Ash Advisory Centers (VAAC) to issue advisories to meteorological watch offices, area control centers, flight information centers, and others. This research aims to provide quantitative forecasts of ash distributions generated by objectively and optimally estimating the volcanic ash source strengths, vertical distribution, and temporal variations using an observation-modeling inversion technique. In this top-down approach, a cost functional is defined to quantify the differences between the model predictions and the satellite measurements of column-integrated ash concentrations weighted by the model and observation uncertainties. Minimizing this cost functional by adjusting the sources provides the volcanic ash emission estimates. As an example, MODIS (Moderate Resolution Imaging Spectroradiometer) satellite retrievals of the 2008 Kasatochi volcanic ash clouds are used to test the HYSPLIT volcanic ash inverse system. Because the satellite retrievals include the ash cloud top height but not the bottom height, there are different model diagnostic choices for comparing the model results with the observed mass loadings. Three options are presented and tested. Although the emission estimates vary significantly with different options, the subsequent model predictions with the different release estimates all show decent skill when evaluated against the unassimilated satellite observations at later times. Among the three options, integrating over three model layers yields slightly better results than integrating from the surface up to the observed volcanic ash cloud top or using a single model layer. Inverse tests also show that including the ash-free region to constrain the model is not beneficial for the current case. In addition, extra constraints on the source terms can be given by explicitly enforcing no-ash for the atmosphere columns above or below the observed ash cloud top height. However, in this case such extra constraints are not helpful for the inverse modeling. It is also found that simultaneously assimilating observations at different times produces better hindcasts than only assimilating the most recent observations.

Sign in / Sign up

Export Citation Format

Share Document