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 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.

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 Characterizing the Atmospheric Conditions Leading to Large Error Growth in Volcanic Ash Cloud Forecasts

2018 ◽  
Vol 57 (4) ◽  
pp. 1011-1019 ◽  
Author(s):  
H. F. Dacre ◽  
N. J. Harvey
Keyword(s):  
Flow Separation ◽  
Volcanic Ash ◽  
Dispersion Model ◽  
Forecast Accuracy ◽  
Large Error ◽  
Ensemble Member ◽  
Forecast Errors ◽  
Atmospheric Conditions ◽  
Wind Fields ◽  
Ensemble Forecasts

ABSTRACTVolcanic ash poses an ongoing risk to safety in the airspace worldwide. The accuracy with which volcanic ash dispersion can be forecast depends on the conditions of the atmosphere into which it is emitted. In this study, meteorological ensemble forecasts are used to drive a volcanic ash transport and dispersion model for the 2010 Eyjafjallajökull eruption in Iceland. From analysis of these simulations, the authors determine why the skill of deterministic-meteorological forecasts decreases with increasing ash residence time and identify the atmospheric conditions in which this drop in skill occurs most rapidly. Large forecast errors are more likely when ash particles encounter regions of large horizontal flow separation in the atmosphere. Nearby ash particle trajectories can rapidly diverge, leading to a reduction in the forecast accuracy of deterministic forecasts that do not represent variability in wind fields at the synoptic scale. The flow‐separation diagnostic identifies where and why large ensemble spread may occur. This diagnostic can be used to alert forecasters to situations in which the ensemble mean is not representative of the individual ensemble‐member volcanic ash distributions. Knowledge of potential ensemble outliers can be used to assess confidence in the forecast and to avoid potentially dangerous situations in which forecasts fail to predict harmful levels of volcanic ash.

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 Assessing the value meteorological ensembles add to dispersion modelling using hypothetical releases

2022 ◽  
Vol 22 (1) ◽  
pp. 577-596
Author(s):  
Susan J. Leadbetter ◽  
Andrew R. Jones ◽  
Matthew C. Hort
Keyword(s):  
Volcanic Ash ◽  
Dispersion Model ◽  
Weather Prediction ◽  
Skill Score ◽  
Decision Makers ◽  
Radioactive Material ◽  
Ensemble Prediction ◽  
Ensemble Simulations ◽  
Model Predictions ◽  
Flight Level

Abstract. Atmospheric dispersion model output is frequently used to provide advice to decision makers, for example, about the likely location of volcanic ash erupted from a volcano or the location of deposits of radioactive material released during a nuclear accident. Increasingly, scientists and decision makers are requesting information on the uncertainty of these dispersion model predictions. One source of uncertainty is in the meteorology used to drive the dispersion model, and in this study ensemble meteorology from the Met Office ensemble prediction system is used to provide meteorological uncertainty to dispersion model predictions. Two hypothetical scenarios, one volcanological and one radiological, are repeated every 12 h over a period of 4 months. The scenarios are simulated using ensemble meteorology and deterministic forecast meteorology and compared to output from simulations using analysis meteorology using the Brier skill score. Adopting the practice commonly used in evaluating numerical weather prediction (NWP) models where observations are sparse or non-existent, we consider output from simulations using analysis NWP data to be truth. The results show that on average the ensemble simulations perform better than the deterministic simulations, although not all individual ensemble simulations outperform their deterministic counterpart. The results also show that greater skill scores are achieved by the ensemble simulation for later time steps rather than earlier time steps. In addition there is a greater increase in skill score over time for deposition than for air concentration. For the volcanic ash scenarios it is shown that the performance of the ensemble at one flight level can be different to that at a different flight level; e.g. a negative skill score might be obtained for FL350-550 and a positive skill score for FL200-350. This study does not take into account any source term uncertainty, but it does take the first steps towards demonstrating the value of ensemble dispersion model predictions.

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 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