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

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.

Atmosphere ◽  
2020 ◽  
Vol 11 (4) ◽  
pp. 342 ◽  
Author(s):  
Meelis J. Zidikheri ◽  
Chris Lucas

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.


2021 ◽  
Author(s):  
Antonio Capponi ◽  
Natalie J. Harvey ◽  
Helen F. Dacre ◽  
Keith Beven ◽  
Mike R. James

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


2011 ◽  
Vol 11 (9) ◽  
pp. 4333-4351 ◽  
Author(s):  
A. Stohl ◽  
A. J. Prata ◽  
S. Eckhardt ◽  
L. Clarisse ◽  
A. Durant ◽  
...  

Abstract. The April–May, 2010 volcanic eruptions of Eyjafjallajökull, Iceland caused significant economic and social disruption in Europe whilst state of the art measurements and ash dispersion forecasts were heavily criticized by the aviation industry. Here we demonstrate for the first time that large improvements can be made in quantitative predictions of the fate of volcanic ash emissions, by using an inversion scheme that couples a priori source information and the output of a Lagrangian dispersion model with satellite data to estimate the volcanic ash source strength as a function of altitude and time. From the inversion, we obtain a total fine ash emission of the eruption of 8.3 ± 4.2 Tg for particles in the size range of 2.8–28 μm diameter. We evaluate the results of our model results with a posteriori ash emissions using independent ground-based, airborne and space-borne measurements both in case studies and statistically. Subsequently, we estimate the area over Europe affected by volcanic ash above certain concentration thresholds relevant for the aviation industry. We find that during three episodes in April and May, volcanic ash concentrations at some altitude in the atmosphere exceeded the limits for the "Normal" flying zone in up to 14 % (6–16 %), 2 % (1–3 %) and 7 % (4–11 %), respectively, of the European area. For a limit of 2 mg m−3 only two episodes with fractions of 1.5 % (0.2–2.8 %) and 0.9 % (0.1–1.6 %) occurred, while the current "No-Fly" zone criterion of 4 mg m−3 was rarely exceeded. Our results have important ramifications for determining air space closures and for real-time quantitative estimations of ash concentrations. Furthermore, the general nature of our method yields better constraints on the distribution and fate of volcanic ash in the Earth system.


Atmosphere ◽  
2020 ◽  
Vol 11 (4) ◽  
pp. 352 ◽  
Author(s):  
Frances M. Beckett ◽  
Claire S. Witham ◽  
Susan J. Leadbetter ◽  
Ric Crocker ◽  
Helen N. Webster ◽  
...  

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.


2019 ◽  
Author(s):  
Soledad Osores ◽  
Juan Ruiz ◽  
Arnau Folch ◽  
Estela Collini

Abstract. Quantitative volcanic ash cloud forecasts are prone to uncertainties coming from the source term quantification (e.g. eruption strength or vertical distribution of the emitted particles), with consequent implications on operational ash impact assessment. We present an ensemble-based data assimilation and forecast system for volcanic ash dispersal and deposition aimed at reducing uncertainties related to eruption source parameters. The FALL3D atmospheric dispersal model is coupled with the Ensemble Transform Kalman Filter (ETKF) data assimilation technique by combining ash mass loading observations with ash dispersal simulations in order to obtain a better joint estimation of 3D ash concentration and source parameters. The ETKF-FALL3D data assimilation system is evaluated performing Observation System Simulation Experiments (OSSE) in which synthetic observations of fine ash mass loadings are assimilated. The evaluation of the ETKF-FALL3D system considering reference states of steady and time-varying eruption source parameters shows that the assimilation process gives both better estimations of ash concentration and time-dependent optimized values of eruption source parameters. The joint estimation of concentrations and source parameters leads to a better analysis and forecast of the 3D ash concentrations. Results show the potential of the methodology to improve volcanic ash cloud forecasts in operational contexts.


2017 ◽  
Vol 3 (3) ◽  
pp. 343
Author(s):  
Cindy Greenman ◽  
Javad Gorjidooz

<p><em>The World Bank Group (WBG) has provided grant and credits to support projects related to aviation and air transport including capacity building, policy and regulation, safety and security, infrastructure rehabilitation, and institutional strengthening for over 60 years. Today, the WBG remains actively engaged on aviation related project in every region of the world including the Pacific Island Countries. The purpose of this paper is to investigate the economic impact of WBG investments in aviation infrastructures in Pacific Island Countries.</em></p><p><em>The preliminary results of this study show a serious deficiency in the areas of aviation infrastructure, safety and security, management efficiency, and airport environmental sustainability in the Pacific Island countries. Without a dynamic aviation industry, most developing countries would see a drastic reduction in tourism, resulting in an economic slowdown. In Pacific Island countries, if these aviation infrastructure deficiency continued to go unaddressed, they would likely lead to the end of international flight operations to the airports in this region. Therefore, the WBG investments in aviation infrastructure is necessary in this region to support tourism, economic development, job creation, and higher standard of living in for Pacific Islanders.</em></p>


2021 ◽  
Author(s):  
Leonardo Mingari ◽  
Arnau Folch ◽  
Andrew T. Prata ◽  
Federica Pardini ◽  
Giovanni Macedonio ◽  
...  

Abstract. Modelling atmospheric dispersal of volcanic ash and aerosols is becoming increasingly valuable for assessing the potential impacts of explosive volcanic eruptions on infrastructures, air quality, and aviation. Management of volcanic risk and reduction of aviation impacts can strongly benefit from quantitative forecasting of volcanic ash. However, an accurate prediction of volcanic aerosol concentrations using numerical modelling relies on proper estimations of multiple model parameters which are prone to errors. Uncertainties in key parameters such as eruption column height, physical properties of particles or meteorological fields, represent a major source of error affecting the forecast quality. The availability of near-real-time geostationary satellite observations with high spatial and temporal resolutions provides the opportunity to improve forecasts in an operational context by incorporating observations into numerical models. Specifically, ensemble-based filters aim at converting a prior ensemble of system states into an analysis ensemble by assimilating a set of noisy observations. Previous studies dealing with volcanic ash transport have demonstrated that a significant improvement of forecast skill can be achieved by this approach. In this work, we present a new implementation of an ensemble-based Data Assimilation (DA) method coupling the FALL3D dispersal model and the Parallel Data Assimilation Framework (PDAF). The FALL3D+PDAF system runs in parallel, supports online-coupled DA and can be efficiently integrated into operational workflows by exploiting high-performance computing (HPC) resources. Two numerical experiments are considered: (i) a twin experiment using an incomplete dataset of synthetic observations of volcanic ash and, (ii) an experiment based on the 2019 Raikoke eruption using real observations of SO2 mass loading. An ensemble-based Kalman filtering technique based on the Local Ensemble Transform Kalman Filter (LETKF) is used to assimilate satellite-retrieved data of column mass loading. We show that this procedure may lead to nonphysical solutions and, consequently, conclude that LETKF is not the best approach for the assimilation of volcanic aerosols. However, we find that a truncated state constructed from the LETKF solution approaches the real solution after a few assimilation cycles, yielding a dramatic improvement of forecast quality when compared to simulations without assimilation.


2016 ◽  
Author(s):  
Tianfeng Chai ◽  
Alice Crawford ◽  
Barbara Stunder ◽  
Michael Pavolonis ◽  
Roland Draxler ◽  
...  

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.


2017 ◽  
Vol 122 (15) ◽  
pp. 8207-8232 ◽  
Author(s):  
Meelis J. Zidikheri ◽  
Christopher Lucas ◽  
Rodney J. Potts

Sign in / Sign up

Export Citation Format

Share Document