R. S. J. Sparks

Quantitative models of the fallout and dispersal of tephra from volcanic eruption columns

A theoretical model of clast fallout from convective eruption columns has been developed which quantifies how the maximum clast size dispersal is determined by column height and wind strength. An eruption column consists of a buoyant convecting region which rises to a heightHB where the column density equals that of the atmosphere. AboveHB the column rises further to a heightHT due to excess momentum. BetweenHT andHB the column is forced laterally into the atmosphere to form an upper umbrella region. Within the eruption column, the vertical and horizontal velocity fields can be calculated from exprimental and theoretical studies and consideration of mass continuity. The centreline vertical velocity falls as a nearly linear function over most of the columns height and the velocity decreases as a gaussian function radially away from the centreline. Both column height and vertical velocity are strong functions of magma discharge rate. From calculations of the velocity field and the terminal fall velocity of clasts, a series of particle support envelopes has been constructed which represents positions where the column vertical velocity and terminal velocity are equal for a clast of specific size and density. The maximum range of a clast is determined in the absence of wind by the maximum width of the clast support envelope.The trajectories of clasts leaving their relevant support envelope at its maximum width have been modelled in columns from 6 to 43 km high with no wind and in a wind field. From these calculations the shapes and areas of maximum grain size contours of the air-fall deposit have been predicted. For the no wind case the theoretical isopleths show good agreement with the Fogo A plinian deposit in the Azores. A diagram has been constructed which plots, for a particular clast size, the maximum range normal to the dispersal axis against the downward range. From the diagram the column height (and hence magma discharge rate) and wind velocity can be determined. Historic plinian eruptions of Santa Maria (1902) and Mount St. Helens (1980) give maximum heights of 34 and 19 km respectively and maximum wind speeds at the tropopause of m/s and 30 m/s respectively. Both estimates are in good agreement with observations. The model has been applied to a number of other plinian deposits, including the ultraplinian phase of theA.D. 180 Taupo eruption in New Zealand which had an estimated column height of 51 km and wind velocity of 27 m/s.

The dimensions and dynamics of volcanic eruption columns

Eruption columns can be divided into three regimes of physical behaviour. The basal gas thrust region is characterized by large velocities and decelerations and is dominated by momentum. This region is typically a few hundred metres in height and passes upwards into a much higher convective region where buoyancy is dominant. The top of the convective region is defined by the level of neutral density (heightHB) where the column has a bulk density equal to the surrounding atmosphere. Above this level the column continues to ascend to a heightHT due to its momentum. The column spreads horizontally and radially outwards between heightHT andHB to form an umbrella cloud. Numerical calculations are presented on the shape of eruption columns and on the relationships between the heightHB and the mass discharge rate of magma, magma temperature and atmospheric temperature gradients. Spreading rate of the column margins increases with height principally due to the decrease in the atmospheric pressure. The relationship between column height and mass discharge rate shows good agreement with observations. The temperature inversion above the tropopause is found to only have a small influence on column height and, eruptions with large discharge rates can inject material to substantially greater heights than the inversion level. Approximate calculations on the variation of convective velocities with height are consistent with field data and indicate that columns typically ascend at velocities from a few tens to over 200 m/s. In very large columns (greater than 30 km) the calculated convective velocities approach the speed of sound in air, suggesting that compressibility effects may become important in giant columns. Radial velocities in the umbrella region where the column is forced laterally into the atmosphere can be substantial and exceed 55 m/s in the case of the May 18th Mount St. Helens eruption. Calculations on motions in this region imply that it plays a major role in the transport of coarse pyroclastic fragments.

The significance of vitric-enriched air-fall ashes associated with crystal-enriched ignimbrites

Abstract A distinctive type of fine-grained air-fall ash is found intimately associated with many ignimbrites. They have crystal/glass ratios systematically lower than artificially crushed pumice from the same ignimbrites. The crystal enrichment found in crystal-bearing ignimbrites indicates substantial losses of the vitric component, amounting to an average of at least 35% by weight of the original juvenile material, and this lost material is believed to occur in the ash-fall deposits. These ashes thus complement ignimbrite, and are here called “co-ignimbrite ashes”. The loss is believed to take place during ignimbrite eruptions as a result of: (1) the escape of fine ash and gas above a collapsing eruptive column; (2) the preferential entry of fine vitric ash into an upper turbulent cloud when (immediately following column collapse) the segregation of a dense pyroclastic flow from an initially highly turbulent, low-concentation density flow takes place; (3) the elutriation of fine vitric ash (generated in part within the pyroclastic flow) from the fluidised flow. Ash from all three mechanisms would be expected to rise to a great height in convective plumes and be dispersed by winds to produce extensive, vitric-enriched ash-fall deposits. The data indicate that the co-ignimbrite ashes must have volumes comparable with those of ignimbrites, and examples are given of particularly large ash-fall deposits (including some found in deep-sea cores) associated with large ignimbrites which may be of this type rather than fall-out from a preceding plinian phase as hitherto thought.

Tephra fallout in the eruption of Soufrière Hills Volcano, Montserrat

Abstract Four mechanisms caused tephra fallout at Soufrière Hills Volcano, Montserrat, during the 1995-1999 period: explosive activity (mainly of Vulcanian type), dome collapses, ash-venting and phreatic explosions. The first two mechanisms contributed most of the tephra-fallout deposits (minimum total dense-rock equivalent volume of 23 x 106 m3), which vary from massive to layered and represent the amalgamation of the deposits from a large numbers of events. The volume of co-pyroclastic-flow fallout tephra is in the range 4-16° of the associated pyroclastic flow deposits. Dome-collapse fallout tephra is characterized by ash particles generated by fragmentation in the pyroclastic flows and by elutriation of fines. Vulcanian fallout tephra is coarser grained, as it is formed by magma fragmentation in the conduit and by elutriation from the fountain-collapse flows and initial surges. Vulcanian fallout tephra is typically polymodal, whereas dome-collapse fallout tephra is predominantly unimodal. Polymodality is attributed to: overlapping of fallout tephra of different types, premature fallout of fine particles, multiple tephra-fallout sources, and differences in density and grain-size distribution of different components. During both dome collapses and explosions, ash fell as aggregates of various sizes and types. Accretionary lapilli grain size is independent of their diameter and is characterized by multiple subpopulations with a main mode at 5ø. Satellite data indicate that very fine ash can stay in a volcanic cloud for several hours and show that exponential thinning rates observed in proximal areas cannot apply in distal areas.

The Los Chocoyos Ash, guatemala: A major stratigraphic marker in middle America and in three ocean basins

Abstract The Los Chocoyos Ash, having erupted from vents near the Lake Atitlan caldera, Guatemala, is perhaps the largest Quaternary silicic pyroclastic unit in Central America. It consists of an underlying H-tephra member and an overlying ash-flow member. One-hundred-and-five samples of ash from the Guatemalan Highlands and deep-sea cores in the equatorial Pacific and Gulf of Mexico were analyzed by neutron activation and/or electron microprobe. Glass shard chemistry, determined by microprobe, is useful for distinguishing several very widespread, distinct, deep-sea ash layers, but needs support from trace-element data when applied on land to distinguish between many individual eruptions from the same province. Data from this study support the correlation of the Worzel ‘D’ layer and the Los Chocoyos Ash proposed by Hahn et al. (1979) and Bowles et al. (1973). Chemical data from this study are used to correlate the Y-8 ash layer of the Gulf of Mexico with the Los Chocoyos Ash. The recognition of the Los Chocoyos Ash in the Gulf of Mexico and equatorial Pacific increases the known areal extent of the unit to more than 6 × 106 km2 and allows an age of 84,000 yr B.P. to be assigned to the formation on the basis of oxygen-isotope stratigraphy, biostratigraphy, and Pa-Th-isotope data. Trace-element data obtained from seven other ash layers in the Gulf of Mexico and the equatorial Pacific, when combined with new land-based data, should allow further correlation and dating of ash units in Central America.

The volcanological significance of deep-sea ash layers associated with ignimbrites

Many volcanic ash layers preserved in deep-sea sediments are the products of large magnitude ignimbrite eruptions. The characteristics of such co-ignimbrite ash-fall deposits are illustrated by two layers from the Eastern Mediterranean: the Minoan ash, Santorini, and the Campanian ash, Italy. These layers are divisible into a coarse lower unit and a fine upper unit in proximal cores. Both layers also show striking bimodal grain size distributions in more distal cores. The coarser mode decreases in median diameter with distance from source whereas the finer mode shows no lateral variation. These features are interpreted in terms of a model for ignimbrite formation by eruption column collapse. Comparable volumes of ignimbrite and associated air-fall ejecta are produced.

Dynamics of giant volcanic ash clouds from supervolcanic eruptions

[1]xa0The largest explosive volcanic eruptions that have occurred on Earth generated giant ash clouds from rising plumes that spread in the stratosphere around a height of neutral buoyancy, with estimated supply rates that are in the range 1011 to 1013 m3/s. These giant ash clouds are controlled by a balance between gravity and Coriolis forces, forming spinning bodies of nearly fixed proportions after a few hours and are initially insensitive to stratospheric winds. In contrast, volcanic plumes from eruptions with small to intermediate magnitude spread as inertial intrusions under the influence of stratospheric winds, with the Earths rotation being unimportant. In the largest eruptions the giant spinning ash clouds typically develop diameters greater than 600 km and up to 6000 km in the most powerful super-eruptions, thus explaining why areas of continental size can be covered with volcanic ash. The radial expansion and spinning velocities are calculated at tens of metres per second and increase with eruption intensity. Higher spreading velocities carry larger ash grain sizes to a given radius, so that grain size at a given distance from the source increases with eruption intensity, consistent with geological observations.

Sedimentology of the Minoan deep-sea tephra layer in the Aegean and Eastern Mediterranean

The Minoan eruption of Santorini resulted in deposition of pyroclastic material over a large area of the Aegean Sea and Eastern Mediterranean. The eruptive activity commenced with a plinian phase, which was followed by a phreatomagmatic phase, and ended with numerous pyroclastic flows. We present a sedimentological study of the Minoan ash deposits found in cores taken throughout the region. The cores are divided into three groups: (1) those from the sub-marine slopes immediately surrounding Santorini; (2) those from the Aegean Sea where the topography consists of many steep-sided basins; and (3) those from the Eastern Mediterranean where the topography is relatively flat and smooth. This last group is designated the Distal Ash Zone. The deposits from the Aegean Sea have been remobilised as turbidity currents. Grain-size distribution analysis reveals that these deposits are bimodal, the Mdo of the coarse population increasing away from source, whilst that of the fine population remains uniform irrespective of distance from source. The bimodality is caused by a combination of aggregation of fine particles in ash clouds produced during the last two phases of the eruptive activity and mixing of the coarse plinian and fine phreatomagmatic and co-ignimbrite deposits within the turbidity current. n nAn unusual feature found in two of the cores is a coarse, density sorted layer beneath the Bouma “a” division of several turbidite units. This is interpreted as an example of a flow head deposit as described by Sparks and Wilson (1983). The deposits from the Distal Ash Fall Zone are fine grained and unimodal, containing phreatomagmatic and co-ignimbrite material. They show only very minor remobilisation of the original air-fall ash, due to the flat topography of the area. The Mdo of the deposit remains uniform irrespective of distance from source and is very similar to that of the fine population found in the bimodal Aegean deposits. This is interpreted as the result of aggregation of fines in the phreatomagmatic and co-ignimbrite ash clouds. The deposits in both areas have been extensively bioturbated. A study of dispersed ash in the cores to reconstruct the original thickness of the ash layers suggests that ash layers less than 0.62 cm thick will be completely dispersed in the surrounding sediment, whilst those thicker than 2.5 cm will be preserved as discrete layers, the base being unaffected by bioturbation. The deposits on the submarine slopes of Santorini were not originally deposited as air-fall ash, but are derived from the erosion of coastal exposures of non-welded ignimbrite. They form storm-sand layers and are typical of the type of deposit that may be formed by rapid erosion of fresh pyroclastic material after the eruption of an island-arc volcano.

Late quaternary explosive silicic volcanism on St Lucia, West Indies

Many explosive eruptions of dacitic magmas have occurred on St Lucia during the late Quaternary. These have produced widespread aprons and fans of pumice flow and ash flow deposits radiating around the central highlands, with co-eruptive air-fall and surge layers interbedded with palaeosols and epiclastic deposits. Vents in the highlands have not been located because of the dense tropical jungle but we suspect they are now plugged by lava domes surrounded by aprons of block and ash flow deposits. Young magmatically related dacitic lava domes have been extruded in the Qualibou depression. The pumice succession can be divided into older quartz-poor deposits forming the Choiseul Pumice and younger crystal-rich deposits with abundant large quartz which are called the Belfond Pumice . The Choiseul Pumice groups together scattered remnants of the products of many eruptions of different low-silica dacitic magma types. The Belfond Pumice is the product of several eruptions of a high-silica magma type and 14 C ages have dated these between 20900 to 34200 years B.P. The pumice flow deposits occur as small-volume valley fills. A granulometric study of Belfond pumice flow deposits shows them to be strongly depleted in finer ash and vitric components. It is suggested that the narrow, winding and vegetated valleys on the island locally induced turbulence and the flows moved with large, highly fluidized and inflated heads, resulting in substantial loss of fine vitric ash. One ash flow deposit which is extremely rich in crystals and carbonized vegetation is highly depleted in fines and shows enhanced vitric losses. This flow may have been a much more violent ash hurricane or blast which surmounted topography ingesting large amounts of lush vegetation. Ignition of this released the large quantities of gas needed to elutriate most of the fines. A model is suggested for the recent volcanic activity on St Lucia in which separate batches of silicic magma, each having a distinctive petrological and chemical character, rose into high level chambers over a large area. Eruptions of volatile-rich magma led to highly explosive pumice-forming activity from vents in the central highlands. Degassed and more crystal-rich magma was extruded later from the same vents or in the attenuated flank of the Qualibou depression to from lava dome complexes.

Morphology and structure of the 1999 lava flows at Mount Cameroon Volcano (West Africa) and their bearing on the emplacement dynamics of volume-limited flows

The morphology and structure of the 1999 lava flows at Mount Cameroon volcano are documented and discussed in relation to local and source dynamics. Structures are analysed qualitatively and more detailed arguments are developed on the processes of levee formation and systematic links between flow dynamics and levee–channel interface geometry. The flows have clear channels bordered by four main types of levees: initial, accretionary, rubble and overflow levees. Thermally immature pahoehoe lava units with overflow drapes define the proximal zone, whereas rubble and accretionary levees are common in the distal region bordering thermally mature aa clinker or blocky aa flow channels. Pressure ridges, squeeze-ups and pahoehoe ropes are the prevalent compressive structures. Standlines displayed on clinkery breccias are interpreted to represent levee–channel interactions in response to changing flow levels. These data complement previous knowledge on lava flow morphology, thus far dominated by Etnean and Hawaiian examples.