Ash fallout scenarios at Vesuvius: Numerical simulations and implications for hazard assessment

Volcanic ash fallout subsequent to a possible renewal of the Vesuvius activity represents a serious threat to the highly urbanized area around the volcano. In order to assess the relative hazard we consider three different possible scenarios such as those following Plinian, Sub-Plinian, and violent Strombolian eruptions. Reference eruptions for each scenario are similar to the 79 AD (Pompeii), the 1631 AD (or 472 AD) and the 1944 AD Vesuvius events, respectively. Fallout deposits for the first two scenarios are modeled using HAZMAP, a model based on a semi-analytical solution of the 2D advection–diffusion–sedimentation equation. In contrast, fallout following a violent Strombolian event is modeled by means of FALL3D, a numerical model based on the solution of the full 3D advection–diffusion–sedimentation equation which is valid also within the atmospheric boundary layer. Inputs for models are total erupted mass, eruption column height, bulk grain-size, bulk component distribution, and a statistical set of wind profiles obtained by the NCEP/NCAR re-analysis. We computed ground load probability maps for different ash loadings. In the case of a Sub-Plinian scenario, the most representative tephra loading maps in 16 cardinal directions were also calculated. The probability maps obtained for the different scenarios are aimed to give support to the risk mitigation strategies.     

Tephra fallout hazard assessment at the Campi Flegrei caldera (Italy)

Tephra fallout associated with renewal of volcanism at the Campi Flegrei caldera is a serious threat to the Neapolitan area. In order to assess the hazards related with tephra loading, we have considered three different eruption scenarios representative of past activity: a high-magnitude event similar to the 4.1 ka Agnano-Monte Spina eruption, a medium-magnitude event, similar to the ∼3.8 ka Astroni 6 eruption, and a low-magnitude event similar to the Averno 2 eruption. The fallout deposits were reconstructed using the HAZMAP computational model, which is based on a semi-analytical solution of the two-dimensional advection–diffusion–sedimentation equation for volcanic tephra. The input parameters into the model, such as total erupted mass, eruption column height, and bulk grain-size and components distribution, were obtained by best-fitting field data. We carried out tens of thousands simulations using a statistical set of wind profiles, obtained from NOAA re-analysis. Probability maps, relative to the considered scenarios, were constructed for several tephra loads, such as 200, 300 and 400 kg/m². These provide a hazard assessment for roof collapses due to tephra loading that can be used for risk mitigation plans in the area.     

An automatic procedure to forecast tephra fallout

Tephra fallout constitutes a serious threat to communities around active volcanoes. Reliable short-term forecasts represent a valuable aid for scientists and civil authorities to mitigate the effects of fallout on the surrounding areas during an episode of crisis. We present a platform-independent automatic procedure with the aim to daily forecast transport and deposition of volcanic particles. The procedure builds on a series of programs and interfaces that automate the data flow and the execution and subsequent postprocess of fallout models. Firstly, the procedure downloads regional meteorological forecasts for the area and time interval of interest, filters and converts data from its native format, and runs the CALMET diagnostic model to obtain the wind field and other micro-meteorological variables on a finer local-scale 3-D grid defined by the user. Secondly, it assesses the distribution of mass along the eruptive column, commonly by means of the radial averaged buoyant plume equations depending on the prognostic wind field and on the conditions at the vent (granulometry, mass eruption rate, etc). All these data serve as input for the fallout models. The initial version of the procedure includes only two Eulerian models, HAZMAP and FALL3D, the latter available as serial and parallel implementations. However, the procedure is designed to incorporate easily other models in a near future with minor modifications on the model source code. The last step is to postprocess the outcomes of models to obtain maps written in standard file formats. These maps contain plots of relevant quantities such as predicted ground load, expected deposit thickness and, for the case of or 3-D models, concentration on air or flight safety concentration thresholds.     

A quantitative uncertainty assessment of eruptive parameters derived from tephra deposits: the example of two large eruptions of Cotopaxi volcano,Ecuador

Physical parameters of explosive eruptions are typically derived from tephra deposits. However, the characterization of a given eruption relies strongly on the quality of the dataset used, the strategy chosen to obtain and process field data and the particular model considered to derive eruptive parameters. As a result, eruptive parameters are typically affected by a certain level of uncertainty and should not be considered as absolute values. Unfortunately, such uncertainty is difficult to assess because it depends on several factors and propagates from field sampling to the application and interpretation of dispersal models. Characterization of explosive eruptions is made even more difficult when tephra deposits are poorly exposed and only medial data are available. In this paper, we present a quantitative assessment of the uncertainty associated with the characterization of tephra deposits generated by the two largest eruptions of the last 2,000 years of Cotopaxi volcano, Ecuador. In particular, we have investigated the effects of the determination of the maximum clast on the compilation of isopleth maps, and, therefore, on the characterization of plume height. We have also compared the results obtained from the application of different models for the determination of both plume height and erupted volume and for the eruption classification. Finally, we have investigated the uncertainty propagation into the calculation of mass eruption rate and eruption duration. We have found that for our case study, the determination of plume height from isopleth maps is more sensitive to the averaging techniques used to define the maximum clast than to the choice of dispersal models used (i.e. models of Carey and Sparks 1986; Pyle 1989) and that even the application of the same dispersal model can result in plume height discrepancies if different isopleth lines are used (i.e. model of Carey and Sparks 1986). However, the uncertainties associated with the determination of erupted mass, and, as a result, of the eruption duration, are larger than the uncertainties associated with the determination of plume height. Mass eruption rate is also associated with larger uncertainties than the determination of plume height because it is related to the fourth power of plume height. Eruption classification is also affected by data processing. In particular, uncertainties associated with the compilation of isopleth maps affect the eruption classification proposed by Pyle (1989), whereas the VEI classification is affected by the uncertainties resulting from the determination of erupted mass. Finally, we have found that analytical and empirical models should be used together for a more reliable characterization of explosive eruptions. In fact, explosive eruptions would be characterized better by a range of parameters instead of absolute values for erupted mass, plume height, mass eruption rate and eruption duration. A standardization of field sampling would also reduce the uncertainties associated with eruption characterization.     

Numerical inversion and analysis of tephra fallout deposits from the 472 AD sub-Plinian eruption at Vesuvius (Italy) through a new best-fit procedure

A simple semi-analytical model for ash-fall deposit was applied to reconstruct the tephra deposits of the sub-Plinian 472 AD eruption of Vesuvius, Italy, which is of the scale of the reference eruptive scenario for the emergency planning, at Vesuvius. Applying a novel least-squares method, the bulk grain-size distribution, the total mass, and the eruption column height were obtained by fitting the computed ground load and granulometries with the observed ones. The analysis of the effect of three different weighting factors in the minimization procedure was also performed. Results showed that the statistical weighting factor produced the minimum bias. The best correlation between calculated and measured deposit was found, even though the quantity of the input data was not very high, as it commonly occurs for several ancient eruptions. Model results were also in agreement with estimations provided by other independent methods.     

Reconstructing eruptive source parameters from tephra deposit: a numerical study of medium-sized explosive eruptions at Etna volcano

Since the 1970s, multiple reconstruction techniques have been proposed and are currently used, to extrapolate and quantify eruptive parameters from sampled tephra fall deposit datasets. Atmospheric transport and deposition processes strongly control the spatial distribution of tephra deposit; therefore, a large uncertainty affects mass derived estimations especially for fall layer that are not well exposed. This paper has two main aims: the first is to analyse the sensitivity to the deposit sampling strategy of reconstruction techniques. The second is to assess whether there are differences between the modelled values for emitted mass and grainsize, versus values estimated from the deposits. We find significant differences and propose a new correction strategy. A numerical approach is demonstrated by simulating with a dispersal code a mild explosive event occurring at Mt. Etna on 24 November 2006. Eruptive parameters are reconstructed by an inversion information collected after the eruption. A full synthetic deposit is created by integrating the deposited mass computed by the model over the computational domain (i.e., an area of 7.5 × 10⁴ km ²). A statistical analysis based on 2000 sampling tests of 50 sampling points shows a large variability, up to 50 % for all the reconstruction techniques. Moreover, for some test examples Power Law errors are larger than estimated uncertainty. A similar analysis, on simulated grain-size classes, shows how spatial sampling limitations strongly reduce the utility of available information on the total grain size distribution. For example, information on particles coarser than ?(?4) is completely lost when sampling at 1.5 km from the vent for all columns with heights less than 2000 m above the vent. To correct for this effect an optimal sampling strategy and a new reconstruction method are presented. A sensitivity study shows that our method can be extended to a wide range of eruptive scenarios including those in which aggregation processes are important. The new correction method allows an estimate of the deficiency for each simulated class in calculated mass deposited, providing reliable estimation of uncertainties in the reconstructed total (whole deposit) grainsize distribution.     

Representivity of incompletely sampled fall deposits in estimating eruption source parameters: a test using the 12–13 January 2011 lava fountain deposit from Mt. Etna volcano,Italy

The Southeast Crater (SEC) of Mt. Etna, Italy, is renowned for its high activity, mainly long-lived eruptions consisting of sequences of individual paroxysmal episodes which have produced more than 150 eruptive events since 1998. Each episode typically forms eruption columns followed by tephra fallout over distances of up to about 100 km from the vent. One of the last sequences consisted of 25 lava fountaining events, which took place between January 2011 and April 2012 from a pit-vent on the eastern flank of the SEC and built a new scoria cone renamed New Southeast Crater. The first episode on 12–13 January 2011 produced tephra fallout which was unusually dispersed toward to the South extending out over the Mediterranean Sea. The southerly deposition of tephra permitted an extensive survey at distances between ~1 and ~100 km, providing an excellent characterization of the tephra deposit. Here, we document the stratigraphy of the 12–13 January fallout deposit, draw its dispersal, and reconstruct its isopleth map. These data are then used to estimate the main eruption source parameters. The total erupted mass (TEM) was calculated by using four different methodologies which give a mean value of 1.5?±?0.4?×?10⁸ kg. The mass eruption rate (MER) is 2.5?±?0.7?×?10⁴ kg/s using eruption duration of 100 min. The total grain-size (TGS) distribution, peaked at ?3 phi, ranges between ?5 and 5 phi and has a median value of ?1.4 phi. Further, for the eruption column height, we obtained respective values of 6.8–13.8 km by using the method of Carey and Sparks (1986) and 3.4?±?0.3 km by using the methods of Wilson and Walker (1987), Mastin et al. (2009), and Pistolesi et al. (2011) and considering the mean value of MER from the deposit. We also evaluated the uncertainty and reliability of TEM and TGS for scenarios where the proximal and distal samples are not obtainable. This is achieved by only using a sector spanning the downwind distances between 6 and 23 km. This scenario is typical for Etna when the tephra plume is dispersed eastward, i.e., in the prevailing wind direction. Our results show that, if the analyzed deposit has poorer sample coverage than presented in this study, the TEM (3.4?×?10⁷ kg) is 22 % than the TEM obtained from the whole deposit. The lack of the proximal (<6 km) deposit may cause more significant differences in the TGS estimations.     

Deformation of Mount Etna, 1971–1974

Electro-optical distance measurements made on the summit of Mt. Etna from 1971 to 1974 show evidence of large surface deformation of the volcano. This deformation cannot be satisfactorily analysed in terms of the models of subsurface magma reservoirs of various geometries that have been previously used, as they have, for instance, on Kilauea in Hawaii. A model that gives a better fit between the observed and computed data involves horizontal, radial strain about an open, cylindrical magma column. In this model, strain is inversely proportional to the square of the distance from the centre of the deformation. This strain pattern is probably confined to the immediate vicinity of the summit vents and is of a different nature lower down the volcano. Tiltmeter, precise levelling and distance measurement data collected over the period of a small flank eruption in January–March 1974 indicate that the eruption was fed by magma through a conduit from the summit reservoir system of the Chasm and Bocca Nuova. Inflation of the summit around the Northeast Crater, which had been measured since 1971, continued despite the flank eruption, and eruptive activity was resumed at the Northeast Crater in September 1974.     

Shifting styles of basaltic explosive activity during the 2002–03 eruption of Mt. Etna, Italy

The 2002–03 flank eruption of Etna was characterized by two months of explosive activity that produced copious ash fallout, constituting a major source of hazard and damage over all eastern Sicily. Most of the tephra were erupted from vents at 2750 and 2800 m elevation on the S flank of the volcano, where different eruptive styles alternated. The dominant style of explosive activity consisted of discrete to pulsing magma jets mounted by wide ash plumes, which we refer to as ash-rich jets and plumes. Similarly, ash-rich explosive activity was also briefly observed during the 2001 flank eruption of Etna, but is otherwise fairly uncommon in the recent history of Etna. Here, we describe the features of the 2002–03 explosive activity and compare it with the 2001 eruption in order to characterize ash-rich jets and plumes and their transition with other eruptive styles, including Strombolian and ash explosions, mainly through chemical, componentry and morphology investigations of erupted ash. Past models explain the transition between different styles of basaltic explosive activity only in terms of flow conditions of gas and liquid. Our findings suggest that the abundant presence of a solid phase (microlites) may also control vent degassing and consequent magma fragmentation and eruptive style. In fact, in contrast with the Strombolian or Hawaiian microlite-poor, fluidal, sideromelane clasts, ash-rich jets and plumes produce crystal-rich tachylite clasts with evidence of brittle fragmentation, suggesting that high groundmass crystallinity of the very top part of the magma column may reduce bubble movement while increasing fragmentation efficiency.     

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

Ash clouds are one of the major hazards that result from volcanic eruptions. Once an eruption is reported, volcanic ash transport and dispersion (VATD) models are used to forecast the location of the ash cloud. These models require source parameters to describe the ash column for initialization. These parameters include: eruption cloud height and vertical distribution, particle size distribution, and start and end time of the eruption. Further, if downwind concentrations are needed, the eruption mass rate and/or volume of ash need to be known. Upon notification of an eruption, few constraints are typically available on many of these source parameters. Recently, scientists have defined classes of eruption types, each with a set of pre-defined eruption source parameters (ESP). We analyze the August 18, 1992 eruption of the Crater Peak vent at Mount Spurr, Alaska, which is the example case for the Medium Silicic eruption type. We have evaluated the sensitivity of two of the ESP – the grain size distribution (GSD) and the vertical distribution of ash – on the modeled ash cloud. HYSPLIT and Puff VATD models are used to simulate the ash clouds from the different sets of source parameters. We use satellite data, processed through the reverse absorption method, as reference for computing statistics that describe the modeled-to-observed comparison. With the grain size distribution, the three options chosen, (1) an estimated distribution based on past eruption studies, (2) a distribution with finer particles and (3) the National Oceanic and Atmospheric Administration HYSPLIT GSD, have little effect on the modeled ash cloud. For the initial vertical distribution, both linear (uniform concentration throughout the vertical column) and umbrella shapes were chosen. For HYSPLIT, the defined umbrella distribution (no ash below the umbrella), apparently underestimates the lower altitude portions of the ash cloud and as a result has a worse agreement with the satellite detected ash cloud compared to that with the linear vertical distribution for this particular eruption. The Puff model, with a Poisson function to represent the umbrella cloud, gave similar results as for a linear distribution, both having reasonable agreement with the satellite detected cloud. Further sensitivity studies of this eruption, as well as studies using the other source parameters, are needed.     

Simulating the dispersal of tephra from the 1991 Pinatubo eruption: Implications for the formation of widespread ash layers

PUFF and HAZMAP, two tephra dispersal models developed for volcanic hazard mitigation, are used to simulate the climatic 1991 eruption of Mt. Pinatubo. PUFF simulations indicate that the majority of ash was advected away from the source at the level of the tropopause (~ 17 km). Several eruptive pulses injected ash and SO₂ gas to higher altitudes (~ 25 km), but these pulses represent only a small fraction (~ 1%) of the total erupted material released during the simulation. Comparison with TOMS images of the SO₂ cloud after 71 and 93 h indicate that the SO₂ gas originated at an altitude of ~ 25 km near the source and descended to an altitude of ~ 22 km as the cloud moved across the Indian Ocean. HAZMAP simulations indicate that the Pinatubo tephra fall deposit in the South China Sea was formed by an eruption cloud with the majority of the ash concentrated at a height of 16–18 km. Results of this study demonstrate that the largest concentration of distal ash was transported at a level significantly below the maximum eruption column height (~ 40 km) and at a level below the calculated height of neutral buoyancy (~ 25 km). Simulations showed that distal ash transport was dominated by atmospheric circulation patterns near the regional tropopause. In contrast, the movement of the SO₂ cloud occurred at higher levels, along slightly different trajectories, and may have resulted from gas/particle segregations that took place during intrusion of the Pinatubo umbrella cloud as it moved away from source.     

Reconstruction and analysis of sub-plinian tephra dispersal during the 1530 A.D. Soufrière (Guadeloupe) eruption: Implications for scenario definition and hazards assessment

The last magmatic eruption of Soufrière of Guadeloupe dated at 1530 A.D. (Soufrière eruption) is characterized by an onset with a partial flank-collapse and emplacement of a debris-avalanche that was followed by a sub-plinian VEI 2–3 explosive short-lived eruption (Phase-1) with a column that reached a height between 9 and 12 km producing about 3.9 × 10⁶ m³ DRE (16.3 × 10⁶ m³ bulk) of juvenile products. The column recurrently collapsed generating scoriaceous pyroclastic flows in radiating valleys up to a distance of 5–6 km with a maximum interpolated bulk deposit volume of 11.7 × 10⁶ m³ (5 × 10⁶ m³ DRE). We have used HAZMAP, a numerical simple first-order model of tephra dispersal [Macedonio, G., Costa, A., Longo, A., 2005. A computer model for volcanic ash fallout and assessment of subsequent hazard. Comput. Geosci. 31, 837–845] to reconstruct to a first approximation the potential dispersal of tephra and associated tephra mass loadings generated by the sub-plinian Phase 1 of the 1530 A.D. eruption. We have tested our model on a deterministic average dry season wind profile that best-fits the available data as well as on a set of randomly selected wind profiles over a 5 year interval that allows the elaboration of probabilistic maps for the exceedance of specific tephra mass load thresholds. Results show that in the hypothesis of a future 1530 A.D. scenario, populated areas to a distance of 3–4 km west–southwest of the vent could be subjected to a static load pressure between 2 and 10 kPa in case of wet tephra, susceptible to cause variable degrees of roof damage. Our results provide volcanological input parameters for scenario and event-tree definition, for assessing volcanic risks and evaluating their impact in case of a future sub-plinian eruption which could affect up to 70 000 people in southern Basse-Terre island and the region. They also provide a framework to aid decision-making concerning land management and development. A sub-plinian eruption is the most likely magmatic scenario in case of a future eruption of this volcano which has shown, since 1992, increasing signs of low-energy seismic, thermal, and acid degassing unrest without significant deformation.     

Sedimentation of tephra by volcanic plumes. Part 2: controls on thickness and grain-size variations of tephra fall deposits

A model for sedimentation from turbulent suspensions predicts that tephra concentration decreases exponentially with time in an ascending volcanic column and in the overlying umbrella cloud. For grain-size distributions typical of plinian eruptions application of the model predicts for thickness variations in good agreement with the exponential thinning observed in tephra fall deposits. The model also predicts a proximal region where fallout from the plume margins results in a more rapid decrease in thickness so that the deposit shows two segments on a thickness versus distance plot. Several examples of deposits with two segments are known. The distance at which the two segments intersect is a measure of eruption column height. The thickness half-distance ( equivalent to the dispersal index of Walker) is strongly correlated with column height, but is also weakly dependent on grain-size distribution of the ejecta. For a dispersal index of 500 km² (the plinian/subplinian boundary of Walker) column heights between 14 and 18 km are calculated. For ultraplinian deposits with D>50000 km² column heights of at least 45 km are implied. Model grain-size distributions of the deposits have sorting values comparable to those observed in tephra fall deposits formed from eruption columns in a weak or negligible cross-wind. Median diameter decreases exponentially with distance as is observed. Sorting () improves with distance as is observed in plinian deposits in a weak wind. However, tephra fall deposits formed in strong winds do not show improved sorting with distance and proximal deposits are typically somewhat better sorted than the model calculations. Differences are attributed to the influence of wind which disperses particles further than predicted in our model and which has an increasing influence as particle size decreases.     

Evaluating suitability of a tephra dispersal model as part of a risk assessment framework

In volcanic risk assessment it is necessary to determine the appropriate level of sophistication for a given predictive model within the contexts of multiple sources of uncertainty and coupling between models. A component of volcanic risk assessment for the proposed radioactive waste repository at Yucca Mountain (Nevada, USA) involves prediction of dispersal of contaminated tephra during violent Strombolian eruptions and the subsequent transport of that tephra toward a hypothetical individual via surface processes. We test the suitability of a simplified model for volcanic plume transport and fallout tephra deposition (ASHPLUME) coupled to a surface sediment-transport model (FAR) that calculates the redistribution of tephra, and in light of inherent uncertainties in the system. The study focuses on two simplifying assumptions in the ASHPLUME model: 1) constant eruptive column height and 2) constant wind speed and direction during an eruption. Variations in tephra dispersal resulting from unsteady column height and wind conditions produced variations up to a factor of two in the concentration of tephra in sediment transported to the control population. However, the effects of watershed geometry and terrain, which control local remobilization of tephra, overprint sensitivities to eruption parameters. Because the combination of models used here shows limited sensitivity to the actual details of ash fall, a simple fall model suffices to estimate tephra mass delivered to the hypothetical individual.     

Petrology of lavas from the 2004–2005 flank eruption of Mt. Etna,Italy: inferences on the dynamics of magma in the shallow plumbing system

Following the 2001 and 2002–2003 flank eruptions, activity resumed at Mt. Etna on 7 September 2004 and lasted for about 6 months. This paper presents new petrographic, major and trace element, and Sr–Nd isotope data from sequential samples collected during the entire 2004–2005 eruption. The progressive change of lava composition allowed defining three phases that correspond to different processes controlling magma dynamics inside the central volcano conduits. The compositional variability of products erupted up to 24 September is well reproduced by a fractional crystallization model that involves magma already stored at shallow depth since the 2002–2003 eruption. The progressive mixing of this magma with a distinct new one rising within the central conduits is clearly revealed by the composition of the products erupted from 24 September to 15 October. After 15 October, the contribution from the new magma gradually becomes predominant, and the efficiency of the mixing process ensures the emission of homogeneous products up to the end of the eruption. Our results give insights into the complex conditions of magma storage and evolution in the shallow plumbing system of Mt. Etna during a flank eruption. Furthermore, they confirm that the 2004–2005 activity at Etna was triggered by regional movements of the eastern flank of the volcano. They caused the opening of a complex fracture zone extending ESE which drained a magma stored at shallow depth since the 2002–2003 eruption. This process favored the ascent of a different magma in the central conduits, which began to be erupted on 24 September without any significant change in eruptive style, deformation, and seismicity until the end of eruption.     

Hazard assessment of far-range volcanic ash dispersal from a violent Strombolian eruption at Somma-Vesuvius volcano, Naples, Italy: implications on civil aviation

Long-range dispersal of volcanic ash can disrupt civil aviation over large areas, as occurred during the 2010 eruption of Eyjafjallaj?kull volcano in Iceland. Here we assess the hazard for civil aviation posed by volcanic ash from a potential violent Strombolian eruption of Somma-Vesuvius, the most likely scenario if eruptive activity resumed at this volcano. A Somma-Vesuvius eruption is of concern for two main reasons: (1) there is a high probability (38?%) that the eruption will be violent Strombolian, as this activity has been common in the most recent period of activity (between AD 1631 and 1944); and (2) violent Strombolian eruptions typically last longer than higher-magnitude events (from 3 to 7?days for the climactic phases) and, consequently, are likely to cause prolonged air traffic disruption (even at large distances if a substantial amount of fine ash is produced such as is typical during Vesuvius eruptions). We compute probabilistic hazard maps for airborne ash concentration at relevant flight levels using the FALL3D ash dispersal model and a statistically representative set of meteorological conditions. Probabilistic hazard maps are computed for two different ash concentration thresholds, 2 and 0.2?mg/m³, which correspond, respectively, to the no-fly and enhanced procedure conditions defined in Europe during the Eyjafjallaj?kull eruption. The seasonal influence of ash dispersal is also analysed by computing seasonal maps. We define the persistence of ash in the atmosphere as the time that a concentration threshold is exceeded divided by the total duration of the eruption (here the eruption phase producing a sustained eruption column). The maps of averaged persistence give additional information on the expected duration of the conditions leading to flight disruption at a given location. We assess the impact that a violent Strombolian eruption would have on the main airports and aerial corridors of the Central Mediterranean area, and this assessment can help those who devise procedures to minimise the impact of these long-lasting low-intensity volcanic events on civil aviation.     

Simulation of the 1980 eruption of Mount St. Helens using the ash-tracking model PUFF

The dispersal of volcanic ash from the May 18, 1980 eruption of Mount St. Helens (MSH) has been simulated using the Lagrangian ash-tracking model PUFF. Previous applications of the model were limited to smaller, short-lived eruptions with ash dispersal occurring mainly within the troposphere. Two high-resolution atmospheric reanalysis datasets (ERA-40 and NCEP/NCAR-40) allowed MSH ash cloud dispersal to be simulated up to 30 km elevation. The 1980 eruption was divided into two distinct eruptive phases, (1) an initial, relatively short-lived blast/surge phase that injected ash up to 30 km and (2) a subsequent nine-hour plinian phase that maintained an average eruption column height of 16 km. Using PUFF, the two phases of the MSH eruption were modeled separately based on a range of individual input parameters and then combined to produce an integrated simulation of the entire eruption. The trajectory and areal extent of the modeled atmospheric ash cloud best match the actual distribution of MSH ash when input parameters are set to values inferred from satellite and radar data collected on May 18, 1980. The prevailing wind field exerts the strongest control on the advection and ultimate position of the modeled ash cloud, making the maximum column height and the vertical distribution of ash the most sensitive of the PUFF input parameters for this event. The results indicate that the PUFF model works well at simulating the dispersal of ash injected well into the lower stratosphere from a moderate, relatively long-lived eruption, such as MSH. However, attempts to use PUFF to recreate some granulometric aspects of the MSH fallout deposit, such as the maximum particle size as a function of distance from source, were not successful. PUFF consistently predicts much greater fallout distances for small ash particles (< 500 µm) than actually observed in the MSH deposit. The effective settling velocities used by the PUFF model appear to be too slow to accurately predict fallout distances of small ash particles. As a consequence the PUFF model may overestimate the duration of ash loading in the atmosphere associated with the distal fine ash component of explosive eruptions.     

A Plinian treatment of fallout from Hawaiian lava fountains

A new model is presented which simulates the dispersal and deposition of material from a Hawaiian eruption column. The model treats the Hawaiian column as a coarse-grained Plinian column and uses a modified version of the Wilson and Walker [Wilson, L., Walker, G.P.L., 1987. Explosive volcanic eruptions: VI. Ejecta dispersal in Plinian eruptions: the control of eruption conditions and atmospheric properties. Geophys. J. R. Astron. Soc. 89, 657–679.] Plinian pyroclast dispersal model to simulate the fall out of material during a Hawaiian eruption. The model results are found to be in good agreement with independent estimates of various parameters made for the 1959 Kilauea Iki eruption of Kilauea volcano. The close agreement between the model results and these independent estimates shows that, dynamically, Hawaiian eruptions are indistinguishable from Plinian eruptions. The major differences in the styles and deposits of these two types of eruptions are accounted for by differences in the mass fluxes and gas contents of the erupting magmas and, most fundamentally, by differences in the grainsize distribution of the erupted clasts. Plume heights predicted by the model are greater than those found for previous models of Hawaiian eruptions. This is because previous models did not allow for the progressive fall out of particles from the plume and, more importantly, made no correction for the velocity disequilibrium between gas and clasts when the grainsize distribution is coarse.