scholarly journals Visualizing Volcanic Ash Forecasts: Scientist and Stakeholder Decisions Using Different Graphical Representations and Conflicting Forecasts

2017 ◽  
Vol 9 (3) ◽  
pp. 333-348 ◽  
Author(s):  
Kelsey J. Mulder ◽  
Matthew Lickiss ◽  
Natalie Harvey ◽  
Alison Black ◽  
Andrew Charlton-Perez ◽  
...  

Abstract During volcanic eruptions, Volcanic Ash Advisory Centres issue ash advisories for aviation showing the forecasted outermost extent of the ash cloud. During the 2010 Icelandic volcano Eyjafjallajökull eruption, the Met Office produced supplementary forecasts of quantitative ash concentration, due to demand from airlines. Additionally, satellite retrievals of estimated volcanic ash concentration are now available. To test how these additional graphical representations of volcanic ash affect flight decisions, whether users infer uncertainty in graphical forecasts of volcanic ash, and how decisions are made when given conflicting forecasts, a survey was conducted of 25 delegates representing U.K. research and airline operations dealing with volcanic ash. Respondents were more risk seeking with safer flight paths and risk averse with riskier flight paths when given location and concentration forecasts compared to when given only the outermost extent of the ash. Respondents representing operations were more risk seeking than respondents representing research. Additionally, most respondents’ hand-drawn no-fly zones were larger than the areas of unsafe ash concentrations in the forecasts. This conservatism implies that respondents inferred uncertainty from the volcanic ash concentration forecasts. When given conflicting forecasts, respondents became more conservative than when given a single forecast. The respondents were also more risk seeking with high-risk flight paths and more risk averse with low-risk flight paths when given conflicting forecasts than when given a single forecast. The results show that concentration forecasts seem to reduce flight cancellations while maintaining safety. Open discussions with the respondents suggested that definitions of uncertainty may differ between research and operations.

2018 ◽  
Vol 18 (6) ◽  
pp. 4019-4038 ◽  
Author(s):  
Alejandro Marti ◽  
Arnau Folch

Abstract. Volcanic ash modeling systems are used to simulate the atmospheric dispersion of volcanic ash and to generate forecasts that quantify the impacts from volcanic eruptions on infrastructures, air quality, aviation, and climate. The efficiency of response and mitigation actions is directly associated with the accuracy of the volcanic ash cloud detection and modeling systems. Operational forecasts build on offline coupled modeling systems in which meteorological variables are updated at the specified coupling intervals. Despite the concerns from other communities regarding the accuracy of this strategy, the quantification of the systematic errors and shortcomings associated with the offline modeling systems has received no attention. This paper employs the NMMB-MONARCH-ASH model to quantify these errors by employing different quantitative and categorical evaluation scores. The skills of the offline coupling strategy are compared against those from an online forecast considered to be the best estimate of the true outcome. Case studies are considered for a synthetic eruption with constant eruption source parameters and for two historical events, which suitably illustrate the severe aviation disruptive effects of European (2010 Eyjafjallajökull) and South American (2011 Cordón Caulle) volcanic eruptions. Evaluation scores indicate that systematic errors due to the offline modeling are of the same order of magnitude as those associated with the source term uncertainties. In particular, traditional offline forecasts employed in operational model setups can result in significant uncertainties, failing to reproduce, in the worst cases, up to 45–70 % of the ash cloud of an online forecast. These inconsistencies are anticipated to be even more relevant in scenarios in which the meteorological conditions change rapidly in time. The outcome of this paper encourages operational groups responsible for real-time advisories for aviation to consider employing computationally efficient online dispersal models.


2021 ◽  
Vol 14 (1) ◽  
pp. 409-436
Author(s):  
Andrew T. Prata ◽  
Leonardo Mingari ◽  
Arnau Folch ◽  
Giovanni Macedonio ◽  
Antonio Costa

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.


2019 ◽  
Author(s):  
Isabelle A. Taylor ◽  
Elisa Carboni ◽  
Lucy J. Ventress ◽  
Tamsin A. Mather ◽  
Roy G. Grainger

Abstract. Ash clouds are a geographically far reaching hazard associated with volcanic eruptions. To minimise the risk that these pose to aircraft and to limit disruption to the aviation industry, it is important to closely monitor the emission and atmospheric dispersion of these plumes. The altitude of the plume is an important consideration and is an essential input into many models of ash cloud propagation. CO2 slicing is an established technique for obtaining the top height of meteorological clouds and previous studies have demonstrated that there is potential for this method to be used for volcanic ash. In this study, the CO2 slicing technique has been adapted for volcanic ash and applied to spectra obtained from the Infrared Atmospheric Sounding Interferometer (IASI). Simulated ash spectra are first used to select the most appropriate channels and then demonstrate that the technique has merit for determining the altitude of the ash. These results indicate a strong match between the true heights and CO2 slicing output with a root mean square error (RMSE) of less than 800 m. Following this, the technique was applied to spectra obtained with IASI during the Eyjafjallajökull and Grimsvötn eruptions in 2010 and 2011 respectively, both of which emitted ash clouds into the troposphere, and which have been extensively studied with satellite imagery. The CO2 slicing results were compared against those from an optimal estimation scheme, also developed for IASI, and a satellite borne LiDAR is used for validation. Overall, the CO2 slicing tool performs better than the optimal estimation scheme. The CO2 slicing heights returned a RMSE value of 2.2 km when compared against the LiDAR. This is lower than the RMSE for the optimal estimation scheme (2.8 km). The CO2 slicing technique is a relatively fast tool and the results suggest that this method could be used to get a first approximation of the ash cloud height, potentially for use for hazard mitigation, or as an input for other retrieval techniques or models of ash cloud propagation.


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.


2013 ◽  
Vol 6 (4) ◽  
pp. 6215-6248 ◽  
Author(s):  
M. Montopoli ◽  
G. Vulpiani ◽  
D. Cimini ◽  
E. Picciotti ◽  
F. S. Marzano

Abstract. The important role played by ground-based microwave weather radars for the monitoring of volcanic ash clouds has been recently demonstrated. The potential of microwaves from satellite passive and ground-based active sensors to estimate near-source volcanic ash cloud parameters has been also proposed, though with little investigation of their synergy and the role of the radar polarimetry. The goal of this work is to show the potentiality and drawbacks of the X-band Dual Polarization radar measurements (DPX) through the data acquired during the latest Grímsvötn volcanic eruptions that took place on May 2011 in Iceland. The analysis is enriched by the comparison between DPX data and the observations from the satellite Special Sensor Microwave Imager/Sounder (SSMIS) and a C-band Single Polarization (SPC) radar. SPC, DPX, and SSMIS instruments cover a large range of the microwaves spectrum, operating respectively at 5.4, 3.2, and 0.16–1.6 cm wavelengths. The multi-source comparison is made in terms of Total Columnar Concentration (TCC). The latter is estimated from radar observables using the "Volcanic Ash Radar Retrieval for dual-Polarization X band systems" (VARR-PX) algorithm and from SSMIS brightness temperature (BT) using a linear BT–TCC relationship. The BT–TCC relationship has been compared with the analogous relation derived from SSMIS and SPC radar data for the same case study. Differences between these two linear regression curves are mainly attributed to an incomplete observation of the vertical extension of the ash cloud, a coarser spatial resolution and a more pronounced non-uniform beam filling effect of SPC measurements (260 km far from the volcanic vent) with respect to the DPX (70 km from the volcanic vent). Results show that high-spatial-resolution DPX radar data identify an evident volcanic plume signature, even though the interpretation of the polarimetric variables and the related retrievals is not always straightforward, likely due to the possible formation of ash and ice particle aggregates and the radar signal depolarization induced by turbulence effects. The correlation of the estimated TCCs derived from DPX and SSMIS BTs reaches −0.73.


2015 ◽  
Vol 57 ◽  
Author(s):  
Frank S. Marzano ◽  
Luigi Mereu ◽  
Mario Montopoli ◽  
Domenico Cimini ◽  
Giovanni Martucci

<div class="page" title="Page 1"><div class="layoutArea"><div class="column"><p><span>Volcanic ash plumes are formed during explosive volcanic eruptions. After advection over several thousands of kilometers, volcanic ash particles are highly fragmented, dispersed and aged with micron- sized sorting. This Annex describes the ash microphysical modeling and the simulated radar and lidar signatures. [...]</span></p></div></div></div>


2017 ◽  
Vol 17 (4) ◽  
pp. 2865-2879 ◽  
Author(s):  
Tianfeng Chai ◽  
Alice Crawford ◽  
Barbara Stunder ◽  
Michael J. Pavolonis ◽  
Roland Draxler ◽  
...  

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.


2016 ◽  
Vol 16 (9) ◽  
pp. 5705-5720 ◽  
Author(s):  
Dimitris Balis ◽  
Maria-Elissavet Koukouli ◽  
Nikolaos Siomos ◽  
Spyridon Dimopoulos ◽  
Lucia Mona ◽  
...  

Abstract. The vulnerability of the European airspace to volcanic eruptions was brought to the attention of the public and the scientific community by the 2010 eruptions of the Icelandic volcano Eyjafjallajökull. As a consequence of this event, ash concentration thresholds replaced the “zero tolerance to ash” rule, drastically changing the requirements on satellite ash retrievals. In response to that, the ESA funded several projects aiming at creating an optimal end-to-end system for volcanic ash plume monitoring and prediction. Two of them, namely the SACS-2 and SMASH projects, developed and improved dedicated satellite-derived ash plume and sulfur dioxide level assessments. The validation of volcanic ash levels and height extracted from the GOME-2 and IASI instruments on board the MetOp-A satellite is presented in this work. EARLINET lidar measurements are compared to different satellite retrievals for two eruptive episodes in April and May 2010. Comparisons were also made between satellite retrievals and aircraft lidar data obtained with the UK's BAe-146-301 Atmospheric Research Aircraft (managed by the Facility for Airborne Atmospheric Measurements, FAAM) over the United Kingdom and the surrounding regions. The validation results are promising for most satellite products and are within the estimated uncertainties of each of the comparative data sets, but more collocation scenes would be desirable to perform a comprehensive statistical analysis. The satellite estimates and the validation data sets are better correlated for high ash optical depth values, with correlation coefficients greater than 0.8. The IASI retrievals show a better agreement concerning the ash optical depth and ash layer height when compared with the ground-based and airborne lidar data.


2017 ◽  
Author(s):  
Alejandro Marti ◽  
Arnau Folch

Abstract. Volcanic ash modeling systems are used to simulate the atmospheric dispersion of volcanic ash and to generate forecasts that quantify the impacts from volcanic eruptions on infrastructures, air quality, aviation, and climate. The efficiency of response and mitigation actions is directly associated to the accuracy of the volcanic ash cloud detection and modeling systems. Operational forecasts build on off-line coupled modeling systems where meteorological variables are updated at the specified coupling intervals. Despite the concerns from other communities regarding the accuracy of this strategy, the quantification of the systematic errors and shortcomings associated to the off-line modeling systems has received no attention. This paper employs the NMMB-MONARCH-ASH model to quantify these errors by employing different quantitative and categorical evaluation scores. The skills of the off-line coupling strategy are compared against those from an on-line forecast considered to be the best estimate of the true outcome. Case studies are considered for a synthetic eruption with constant eruption source parameters and for two historical events, which suitably illustrate the severe aviation disruptive effects of European (2010 Eyjafjallajökull) and South-American (2011 Cordón Caulle) volcanic eruptions. Evaluation scores indicate that systematic errors credited to off-line modeling are of the same order of magnitude that those associated to the source term uncertainties. In particular, traditional off-line forecasts employed in operational model setups can result in significant uncertainties, failing to reproduce, in the worst cases, up to 45–70 % of the ash cloud of an on-line forecast. These inconsistencies are anticipated to be even more relevant in scenarios where the meteorological conditions change rapidly in time. The outcome of this paper encourages operational groups responsible for real‐time advisories for aviation to consider employing computationally efficient on-line dispersal models.


2019 ◽  
Vol 12 (7) ◽  
pp. 3853-3883 ◽  
Author(s):  
Isabelle A. Taylor ◽  
Elisa Carboni ◽  
Lucy J. Ventress ◽  
Tamsin A. Mather ◽  
Roy G. Grainger

Abstract. Ash clouds are a geographically far-reaching hazard associated with volcanic eruptions. To minimise the risk that these pose to aircraft and to limit disruption to the aviation industry, it is important to closely monitor the emission and atmospheric dispersion of these plumes. The altitude of the plume is an important consideration and is an essential input into many models of ash cloud propagation. CO2 slicing is an established technique for obtaining the top height of aqueous clouds, and previous studies have demonstrated that there is potential for this method to be used for volcanic ash. In this study, the CO2 slicing technique has been adapted for volcanic ash and applied to spectra obtained from the Infrared Atmospheric Sounding Interferometer (IASI). Simulated ash spectra are first used to select the most appropriate channels and then demonstrate that the technique has merit for determining the altitude of the ash. These results indicate a strong match between the true heights and CO2 slicing output with a root mean square error (RMSE) of less than 800 m. Following this, the technique was applied to spectra obtained with IASI during the Eyjafjallajökull and Grímsvötn eruptions in 2010 and 2011 respectively, both of which emitted ash clouds into the troposphere, and which have been extensively studied with satellite imagery. The CO2 slicing results were compared against those from an optimal estimation scheme, also developed for IASI, and a satellite-borne lidar is used for validation. The CO2 slicing heights returned an RMSE value of 2.2 km when compared against the lidar. This is lower than the RMSE for the optimal estimation scheme (2.8 km). The CO2 slicing technique is a relatively fast tool and the results suggest that this method could be used to get a first approximation of the ash cloud height, potentially for use for hazard mitigation, or as an input for other retrieval techniques or models of ash cloud propagation.


Sign in / Sign up

Export Citation Format

Share Document