A statistical approach to evaluate the tephra deposit and ash concentration from PUFF model forecasts

Abstract In this paper we present a new statistical approach which provides tephra deposit load and ash concentration using PUFF, a Lagrangian model widely used to forecast volcanic ash dispersal during volcanic crisis. We perform a parametric study in order to analyze the influence of each input parameter on model outputs. For this test, we simulate two eruptive scenarios similar to the 2001 ( Scenario 1 ) and 1998 ( Scenario 2 ) Etna eruptions using high resolution weather data and a domain of 170 × 170 km. Results show that for both scenarios, we are able to calculate the tephra deposit load and ash concentration but the use of millions of particles is required. Specifically, up to 33 and 220 millions of particles were necessary to accurately predict the tephra deposit and ash concentration in air, respectively. This is approximately two orders of magnitude larger than the values typically considered running PUFF. The parametric study shows that the horizontal diffusion coefficient, the time step of the simulations, the topography and the standard deviation of the particle distribution greatly affect the model outputs. We also validate the model by best-fit procedures. Results show a good comparison between field data of the 2001 Etna eruption and PUFF simulations, being inside 5 and 1/5 times the observed data, comparable with results of Eulerian models. This work will allow to reliably outline the areas of contaminated airspace using PUFF or any other Lagrangian model in order to define the No Fly Zone and ensure the safety aviation operations as required after the Eyjafjallajokull eruption.

[1]  G. Macedonio,et al.  A model for the numerical simulation of tephra fall deposits , 2005 .

[2]  Barbara J. B. Stunder,et al.  Preliminary sensitivity study of eruption source parameters for operational volcanic ash cloud transport and dispersion models — A case study of the August 1992 eruption of the Crater Peak vent, Mount Spurr, Alaska , 2009 .

[3]  C. Bonadonna,et al.  Numerical modelling of tephra fallout associated with dome collapses and Vulcanian explosions: application to hazard assessment on Montserrat , 2002, Geological Society, London, Memoirs.

[4]  I. Charpentier.  Variational coupling of Plinian column models and data: Application to El Chichón volcano , 2008 .

[5]  Jan-Peter Muller,et al.  Assessment of MISR and MODIS cloud top heights through inter‐comparison with a back‐scattering lidar at SIRTA , 2004 .

[6]  David G. Stork,et al.  Pattern Classification , 1973 .

[7]  Costanza Bonadonna,et al.  Modeling tephra sedimentation from a Ruapehu weak plume eruption , 2005 .

[8]  William I. Rose,et al.  Hydrometeor-enhanced tephra sedimentation: Constraints from the 18 May 1980 eruption of Mount St. Helens , 2009 .

[9]  Simona Scollo,et al.  Tephra fallout of 2001 Etna flank eruption: Analysis of the deposit and plume dispersion , 2007 .

[10]  Larry G. Mastin,et al.  A multidisciplinary effort to assign realistic source parameters to models of volcanic ash-cloud transport and dispersion during eruptions , 2009 .

[11]  Augusto Neri,et al.  High performance computing simulations of pyroclastic flows , 2005, Comput. Phys. Commun..

[12]  T. Casadevall,et al.  Volcanic ash and aviation safety; proceedings of the First international symposium on Volcanic ash and aviation safety , 1994 .

[13]  William I. Rose,et al.  Quantitative shape measurements of distal volcanic ash , 2003 .

[14]  Saltelli Andrea,et al.  Global Sensitivity Analysis: The Primer , 2008 .

[15]  J. Merrill,et al.  Simulating the dispersal of tephra from the 1991 Pinatubo eruption: Implications for the formation of widespread ash layers , 2009 .

[16]  R. S. J. Sparks,et al.  Thickness variations and volume estimates of tephra fall deposits: the importance of particle Reynolds number , 1998 .

[17]  Nicola Spinelli,et al.  Applications of the PUFF model to forecasts of volcanic clouds dispersal from Etna and Vesuvio , 2009, Comput. Geosci..

[18]  K. Dean,et al.  PUFF: A high-resolution volcanic ash tracking model , 1998 .

[19]  S. Scollo,et al.  3D volcanic aerosol dispersal: a comparison between misr data and numerical simulations , 2010 .

[20]  S. Scollo,et al.  Three‐dimensional volcanic aerosol dispersal: A comparison between Multiangle Imaging Spectroradiometer (MISR) data and numerical simulations , 2010 .

[21]  S. Carey,et al.  Simulations of tephra dispersal from the 1991 explosive eruptions of Hudson volcano, Chile , 2010 .

[22]  J. Pelletier,et al.  Evaluating suitability of a tephra dispersal model as part of a risk assessment framework , 2008 .

[23]  Stefano Tarantola,et al.  Sensitivity analysis and uncertainty estimation for tephra dispersal models , 2008 .

[24]  Simona Scollo,et al.  Monitoring ash emission episodes at Mt. Etna : The 16 November 2006 case study , 2009 .

[25]  Giancarlo Neri,et al.  Satellite analysis and PUFF simulation of the eruptive cloud generated by the Mount Etna paroxysm of 22 July 1998 , 2002 .

[26]  Arnau Folch,et al.  A three-dimensional Eulerian model for transport and deposition of volcanic ashes , 2006 .

[27]  Y. Ishimine.  A simple integral model of buoyancy-generating plumes and its application to volcanic eruption columns , 2007 .

[28]  Arlene Laing,et al.  Probabilistic modeling of tephra dispersal: Hazard assessment of a multiphase rhyolitic eruption at Tarawera, New Zealand , 2005 .

[29]  R. Sparks,et al.  The eruption of Soufrière Hills Volcano, Montserrat, from 1995 to 1999 , 2002 .

[30]  Simona Scollo,et al.  Characterization of shape and terminal velocity of tephra particles erupted during the 2002 eruption of Etna volcano, Italy , 2008 .

[31]  J. Turner,et al.  Buoyancy Effects in Fluids , 1973 .

[32]  H. Tanaka,et al.  Numerical simulation of volcanic plume dispersal from Usu volcano in Japan on 31 March 2000 using PUFF model , 2002 .

[33]  Carlo Cavazzoni,et al.  An automatic procedure to forecast tephra fallout , 2008 .

[34]  Arnau Folch,et al.  A parametric and comparative study of different tephra fallout models , 2008 .

[35]  Jonathan Dehn,et al.  Predicting Regions Susceptible to High Concentrations of Airborne Volcanic Ash in the North Pacific Region (Abstract Only) , 2004 .

[36]  Maria Teresa Pareschi,et al.  Role of local wind circulation in plume monitoring at Mt. Etna volcano (Sicily): Insights from a mesoscale numerical model , 2004 .

[37]  M. Coltelli,et al.  Monitoring and forecasting Etna volcanic plumes , 2009 .

[38]  Antonella Longo,et al.  A computer model for volcanic ash fallout and assessment of subsequent hazard , 2005, Comput. Geosci..

[39]  S. Barnard.  Results of a reconnaissance trip to Mt. Etna, Italy: , 2004 .

[40]  Augusto Neri,et al.  The VOL‐CALPUFF model for atmospheric ash dispersal: 2. Application to the weak Mount Etna plume of July 2001 , 2008 .

[41]  B. Voight,et al.  The eruption of Soufrière Hills Volcano, Montserrat from 1995 to 1999 , 2002 .

[42]  J. Merrill,et al.  Simulation of the 1980 eruption of Mount St. Helens using the ash-tracking model PUFF , 2008 .

[43]  Jonathan Dehn,et al.  Integrated satellite observations of the 2001 eruption of Mt. Cleveland, Alaska , 2004 .