An extremely large magnitude eruption close to the Plio-Pleistocene boundary: reconstruction of eruptive style and history of the Ebisutoge-Fukuda tephra, central Japan

An extremely large magnitude eruption of the Ebisutoge-Fukuda tephra, close to the Plio-Pleistocene boundary, central Japan, spread volcanic materials widely more than 290,000 km² reaching more than 300 km from the probable source. Characteristics of the distal air-fall ash (>150 km away from the vent) and proximal pyroclastic deposits are clarified to constrain the eruptive style, history, and magnitude of the Ebisutoge-Fukuda eruption.Eruptive history had five phases. Phase 1 is phreatoplinian eruption producing >105 km³ of volcanic materials. Phases 2 and 3 are plinian eruption and transition to pyroclastic flow. Plinian activity also occurred in phase 4, which ejected conspicuous obsidian fragments to the distal locations. In phase 5, collapse of eruption column triggered by phase 4, generated large pyroclastic flow in all directions and resulted in more than 250–350 km³ of deposits. Thus, the total volume of this tephra amounts over 380–490 km³. This indicates that the Volcanic Explosivity Index (VEI) of the Ebisutoge-Fukuda tephra is greater than 7. The huge thickness of reworked volcaniclastic deposits overlying the fall units also attests to the tremendous volume of eruptive materials of this tephra.Numerous ancient tephra layers with large volume have been reported worldwide, but sources and eruptive history are often unknown and difficult to determine. Comparison of distal air-fall ashes with proximal pyroclastic deposits revealed eruption style, history and magnitude of the Ebisutoge-Fukuda tephra. Hence, recognition of the Ebisutoge-Fukuda tephra, is useful for understanding the volcanic activity during the Pliocene to Pleistocene, is important as a boundary marker bed, and can be used to interpret the global environmental and climatic impact of large magnitude eruptions in the past.     

Volcanic plumes and wind: Jetstream interaction examples and implications for air traffic

Volcanic plumes interact with the wind at all scales. On smaller scales, wind affects local eddy structure; on larger scales, wind shapes the entire plume trajectory. The polar jets or jetstreams are regions of high [generally eastbound] winds that span the globe from 30 to 60° in latitude, centered at an altitude of about 10 km. They can be hundreds of kilometers wide, but as little as 1 km in thickness. Core windspeeds are up to 130 m/s. Modern transcontinental and transoceanic air routes are configured to take advantage of the jetstream. Eastbound commercial jets can save both time and fuel by flying within it; westbound aircraft generally seek to avoid it.Using both an integral model of plume motion that is formulated within a plume-centered coordinate system (BENT) as well as the Active Tracer High-resolution Atmospheric Model (ATHAM), we have calculated plume trajectories and rise heights under different wind conditions. Model plume trajectories compare well with the observed plume trajectory of the Sept 30/Oct 1, 1994, eruption of Kliuchevskoi Volcano, Kamchatka, Russia, for which measured maximum windspeed was 30–40 m/s at about 12 km. Tephra fall patterns for some prehistoric eruptions of Avachinsky Volcano, Kamchatka, and Inyo Craters, CA, USA, are anomalously elongated and inconsistent with simple models of tephra dispersal in a constant windfield. The Avachinsky deposit is modeled well by BENT using a windspeed that varies with height.Two potentially useful conclusions can be made about air routes and volcanic eruption plumes under jetstream conditions. The first is that by taking advantage of the jetstream, aircraft are flying within an airspace that is also preferentially occupied by volcanic eruption clouds and particles. The second is that, because eruptions with highly variable mass eruption rate pump volcanic particles into the jetstream under these conditions, it is difficult to constrain the tephra grain size distribution and mass loading present within a downwind volcanic plume or cloud that has interacted with the jetstream. Furthermore, anomalously large particles and high mass loadings could be present within the cloud, if it was in fact formed by an eruption with a high mass eruption rate. In terms of interpretation of tephra dispersal patterns, the results suggest that extremely elongated isopach or isopleth patterns may often be the result of eruption into the jetstream, and that estimation of the mass eruption rate from these elongated patterns should be considered cautiously.     

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.     

A model for the numerical simulation of tephra fall deposits

A simple semianalytical model to simulate ash dispersion and deposition produced by sustained Plinian and sub-Plinian eruption columns based on the 2D advection–dispersion equation was applied. The eruption column acts as a vertical line source with a given mass distribution and neglects the complex dynamics within the eruption column. Thus, the use of the model is limited to areas far from the vent where the dynamics of the eruption column play a minor role. Vertical wind and diffusion components are considered negligible with respect to the horizontal ones. The dispersion and deposition of particles in the model is only governed by gravitational settling, horizontal eddy diffusion, and wind advection. The model accounts for different types and size classes of a user-defined number of particle classes and changing settling velocity with altitude. In as much as wind profiles are considered constant on the entire domain, the model validity is limited to medium-range distances (about 30–200 km away from the source).The model was used to reconstruct the tephra fall deposit from the documented Plinian eruption of Mt. Vesuvius, Italy, in 79 A.D. In this case, the model was able to broadly reproduce the characteristic medium-range tephra deposit. The results support the validity of the model, which has the advantage of being simple and fast to compute. It has the potential to serve as a simple tool for predicting the distribution of ash fall of hypothetical or real eruptions of a given magnitude and a given wind profile. Using a statistical set of frequent wind profiles, it also was used to construct air fall hazard maps of the most likely affected areas around active volcanoes where a large eruption is expected to occur.     

Doppler radar sounding of volcanic eruption dynamics at Mount Etna

Besides their common use in atmospheric studies, Doppler radars are promising tools for the active remote sensing of volcanic eruptions but were little applied to this field. We present the observations made with a mid-power UHF Doppler radar (Voldorad) during a 7-h Strombolian eruption at the SE crater of Mount Etna on 11–12 October 1998. Main characteristics of radar echoes are retrieved from analysis of Doppler spectra recorded in the two range gates on either side of the jet axis. From the geometry of the sounding, the contribution of uprising and falling ejecta to each Doppler spectrum can be discriminated. The temporal evolution of total power backscattered by uprising targets is quite similar to the temporal evolution of the volcanic tremor and closely reproduces the overall evolution of the eruption before, during and after its paroxysm. Moreover, during the sharp decrease of eruptive activity following the paroxysm, detailed analysis of video (from camera recording), radar and seismic measurements reveals that radar and video signals start to decrease simultaneously, approximately 2.5 min after the tremor decline. This delay is interpreted as the ascent time through a magma conduit of large gas slugs from a shallow source roughly estimated at about 500 m beneath the SE crater. Detailed analysis of eruptive processes has been also made with Voldorad operating in a high sampling rate mode. Signature of individual outburst is clearly identified on the half part of Doppler spectra corresponding to rising ejecta: temporal variations of the backscattered power exhibit quasi periodic undulations, whereas the maximum velocity measured on each spectrum displays a sharp peak at the onset of each outburst followed by a slow decay with time. Periodicity of power variations (between 3.8 and 5.5 s) is in agreement with the occurrence of explosions visually observed at the SE vent. Maximum vertical velocities of over 160 m s^–1 were measured during the paraoxysmal stage and the renewed activity. Finally, by using a simplified model simulating the radar echoes characteristics, we show that when Voldorad is operating in high sampling rate mode, the power and maximum velocity variations are directly related to the difference in size and velocity of particles crossing the antenna beam.Editorial responsibility: A. Woods     

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.     

Plinian activity during the early eruptive history of the Sabatini Volcanic District, Central Italy

The early activity of the Sabatini Volcanic District (SVD; central Italy) was characterised by highly explosive eruptions that produced widespread subplinian and plinian fall deposits. In this study, four major eruptive units—informally named as units A, B, C and D—were recognised in the 514–449 ka age interval. In particular, unit D was emplaced during the early phase of the 449 ka Tufo Rosso a Scorie Nere pyroclastic flow-forming eruption, the most important event in the whole SVD activity history. Estimates of relevant eruptive parameters indicate tephra fall volumes up to 4 km³ for individual units, peak eruption column heights in the range of 14–29 km and corresponding mass eruption rates of 7.8×10⁶–1.3×10⁸ kg/s. Isopach and isopleth maps of fallout deposits—as well as the distribution of the proximal lag-breccia of the Tufo Rosso a Scorie Nere—consistently indicate a common vent area, which does not correspond to any volcanic centre identified up to now in the SVD. This was located along NE–SW-trending tectonic lineaments that also controlled the location of the other major volcanic centres of the SVD. The characterisation by means of field aspects, grain size, componentry and density and chemical composition of juvenile clasts, renders the studied fall deposits as valuable stratigraphic markers for the SVD and well beyond it. In fact, their wide areal dispersals toward the E and SE may allow correlations on a regional scale for the Quaternary successions of intermountain basins of the Central Apennine and the Adriatic Sea basin successions. Finally, the correct identification of distal tephra from plinian and co-ignimbrite plumes and their attribution to specific explosive eruptions of the SVD and the other volcanic districts of the Roman Province—rather than to local intra-Apennine centres—provides crucial implications for geodynamic reconstructions.     

Volcanic risk ranking for Auckland, New Zealand. II: Hazard consequences and risk calculation

In a companion paper, a methodology for ranking volcanic hazards and events in terms of risk was presented, and the likelihood and extent of potential hazards in the Auckland Region, New Zealand investigated. In this paper, the effects of each hazard are considered and the risk ranking completed. Values for effect are proportions of total loss and, as with likelihood and extent, are based on order of magnitude.Two outcomes were considered – building damage and loss of human life. In terms of building damage, tephra produces the highest risk by an order of magnitude, followed by lava flows and base surge. For loss of human life, risk from base surge is highest. The risks from pyroclastic flows and tsunami are an order of magnitude smaller. When combined, tephra fall followed by base surge produces the highest risk. The risks from lava flows and pyroclastic flows are an order of magnitude smaller. For building damage, the risk from Mt. Taranaki volcano, 280 km from the Auckland CBD, is largest, followed by Okataina volcanic centre, an Auckland volcanic field eruption centred on land, then Tongariro volcanic centre. In terms of human loss, the greatest risk is from an Auckland eruption centred on land. The risks from an Auckland eruption centred in the ocean, Okataina volcanic centre, and Taupo volcano are more than an order of magnitude smaller. When combined, the risk from Mt. Taranaki remains highest, followed by an Auckland eruption centred on land. The next largest risks are from the Okataina and Tongariro volcanic centres, followed by Taupo volcano.Three alternative situations were investigated. As multiple eruptions may occur from the Auckland volcanic field, it was assumed that a local event would involve two eruptions. This increased risk of a local eruption occurring on land so that it was equal to that of an eruption from Mt. Taranaki. It is possible that a future eruption may be of a similar, or larger size, to the previous Rangitoto eruption. Risk was re-calculated for local eruptions based on the extent of hazards from Rangitoto. This increased the risk of lava flow to greater than that of base surge, and the risk from an Auckland land eruption became greatest. The relative probabilities used for Mt. Taranaki volcano and the Auckland volcanic field may only be minimum values. When the probability of these occurring was increased by 50%, there was no change in either ranking.Editorial responsibility: J. S. Gilbert     

The interaction between volcanic gases and tephra: Fluorine adhering to tephra of the 1970 hekla eruption

The mass distribution and sorting of tephra produced in the plinian phase of the 1970 Hekla eruption was controlled by the particle size distribution, the height of the eruption column, and velocity of transport. Near the volcano the mass distribution of soluble fluorine was controlled by particle size of the deposits, but approaches the mass distribution of the tephra at longer distances. Adsorbed soluble fluorine reaches a maximum at a distance from the volcano determined by the velocity of the transporting medium.SEM studies show the soluble fluorine to be chemically adsorbed on the surface of tephra particles. The adsorption is shown by experiment to occur at temperatures below 600°C in the cooling eruption column. Evaluation of reactions in the eruption column leads to the conclusion that formation of water soluble compounds adhering to tephra is principally controlled by environmental factors and to a lesser degree by the composition of the volcanic gas phase.     

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.     

Nature and significance of small volume fall deposits at composite volcanoes: Insights from the October 14, 1974 Fuego eruption,Guatemala

The first of four successive pulses of the 1974 explosive eruption of Fuego volcano, Guatemala, produced a small volume (∼0.02 km³ DRE) basaltic sub-plinian tephra fall and flow deposit. Samples collected within 48 h after deposition over much of the dispersal area (7–80 km from the volcano) have been size analyzed down to 8 φ (4 μm). Tephra along the dispersal axis were all well-sorted (σ _φ = 0.25–1.00), and sorting increased whereas thickness and median grain size decreased systematically downwind. Skewness varied from slightly positive near the vent to slightly negative in distal regions and is consistent with decoupling between coarse ejecta falling off the rising eruption column and fine ash falling off the windblown volcanic cloud advecting at the final level of rise. Less dense, vesicular coarse particles form a log normal sub-population when separated from the smaller (Md_φ < 3φ or < 0.125 mm), denser shard and crystal sub-population. A unimodal, relatively coarse (Md_φ = 0.58φ or 0.7 mm σ _φ = 1.2) initial grain size population is estimated for the whole (fall and flow) deposit. Only a small part of the fine-grained, thin 1974 Fuego tephra deposit has survived erosion to the present day. The initial October 14 pulse, with an estimated column height of 15 km above sea level, was a primary cause of a detectable perturbation in the northern hemisphere stratospheric aerosol layer in late 1974 to early 1975. Such small, sulfur-rich, explosive eruptions may substantially contribute to the overall stratospheric sulfur budget, yet leave only transient deposits, which have little chance of survival even in the recent geologic record. The fraction of finest particles (Md_φ = 4–8φ or 4–63 μm) in the Fuego tephra makes up a separate but minor size mode in the size distribution of samples around the margin of the deposit. A previously undocumented bimodal–unimodal–bimodal change in grain size distribution across the dispersal axis at 20 km downwind from the vent is best accounted for as the result of fallout dispersal of ash from a higher subplinian column and a lower “co-pf” cloud resulting from pyroclastic flows. In addition, there is a degree of asymmetry in the documented grain-size fallout pattern which is attributed to vertically veering wind direction and changing windspeeds, especially across the tropopause. The distribution of fine particles (<8 μm diameter) in the tephra deposit is asymmetrical, mainly along the N edge, with a small enrichment along the S edge. This pattern has hazard significance.     

Mortality in England during the 1783–4 Laki Craters eruption

1783/4 has been recognised as a mortality crisis year in the population history of England. This demographic incident coincides with the Laki Craters eruption, Iceland, which began in June 1783 and fumigated many parts of Europe with volcanic gases and particles. Many reports and proxy climate records implicate the volcanic cloud in meteorological anomalies, including notably hot 1783 summer conditions in England and a severe subsequent winter. We present here a detailed analysis of the geographical and temporal trends in English mortality data, and interpret them in the light of the climatological records and observations of the pollutant cloud. We show that there were two distinct crisis periods: in August-September 1783, and January-February 1784, which together accounted for ~20,000 extra deaths. In both cases, the East of England was the worst affected region. Possible causes for the two crisis periods are considered and we conclude that the timing and magnitude of the winter mortality peak can be explained by the severe cold of January 1784. The late summer mortality followed 1–2 months after the very hot July of 1783 and may also have been related to the weather, with the time lag reflecting the relatively slow spread of enteric disease or the contraction of malaria. However, it is hard to explain the entire late summer anomaly by these high temperature causes. We therefore consider that fine acid aerosol and/or gases in the volcanic haze may also have contributed to the unusual August-September mortality. Given that complex radiative and dynamical effects of the volcanic cloud are implicated in the climatic anomalies in 1783–4, it is likely that the Laki Craters eruption did play a role in the English mortality crises of the same period.Editorial responsibility: R. Cioni     

Explosive silicic volcanism in Iceland and the Jan Mayen area during the last 6 Ma: sources and timing of major eruptions

Fifty-three major explosive eruptions on Iceland and Jan Mayen island were identified in 0–6-Ma-old sediments of the North Atlantic and Arctic oceans by the age and the chemical composition of silicic tephra. The depositional age of the tephra was estimated using the continuous record in sediment of paleomagnetic reversals for the last 6 Ma and paleoclimatic proxies (δ¹⁸O, ice-rafted debris) for the last 1 Ma. Major element and normative compositions of glasses were used to assign the sources of the tephra to the rift and off-rift volcanic zones in Iceland, and to the Jan Mayen volcanic system. The tholeiitic central volcanoes along the Iceland rift zones were steadily active with the longest interruption in activity recorded between 4 and 4.9 Ma. They were the source of at least 26 eruptions of dominant rhyolitic magma composition, including the late Pleistocene explosive eruption of Krafla volcano of the Eastern Rift Zone at about 201 ka. The central volcanoes along the off-rift volcanic zones in Iceland were the source of at least 19 eruptions of dominant alkali rhyolitic composition, with three distinct episodes recorded at 4.6–5.3, 3.5–3.6, and 0–1.8 Ma. The longest and last episode recorded 11 Pleistocene major events including the two explosive eruptions of Tindfjallajökull volcano (Thórsmörk, ca. 54.5 ka) and Katla volcano (Sólheimar, ca. 11.9 ka) of the Southeastern Transgressive Zone. Eight major explosive eruptions from the Jan Mayen volcanic system are recorded in terms of the distinctive grain-size, mineralogy and chemistry of the tephra. The tephra contain K-rich glasses (K₂O/SiO₂>0.06) ranging from trachytic to alkali rhyolitic composition. Their normative trends (Ab–Q–Or) and their depleted concentrations of Ba, Eu and heavy-REE reflect fractional crystallisation of K-feldspar, biotite and hornblende. In contrast, their enrichment in highly incompatible and water-mobile trace elements such as Rb, Th, Nb and Ta most likely reflect crustal contamination. One late Pleistocene tephra from Jan Mayen was recorded in the marine sequence. Its age, estimated between 617 and 620 ka, and its composition support a common source with the Borga pumice formation at Sør Jan in the south of the island.     

The 1991 eruption of Hekla,Iceland

The eruption that started in the Hekla volcano in South Iceland on 17 January 1991, and came to an end on 11 March, produced mainly andesitic lava. This lava covers 23 km² and has an estimated volume of 0.15 km³. This is the third eruption in only 20 years, whereas the average repose period since 1104 is 55 years. Earthquakes, as well as a strain pulse recorded by borehole strainmeters, occurred less than half an hour before the start of the eruption. The initial plinian phase was very short-lived, producing a total of only 0.02 km³ of tephra. The eruption cloud attained 11.5 km in height in only 10 min, but it became detached from the volcano a few hours later. Several fissures were active during the first day of the eruption, including a part of the summit fissure. By the second day, however, the activity was already essentially limited to that segment of the principal fissure where the main crater subsequently formed. The average effusion rate during the first two days of the eruption was about 800 m³ s^–1. After this peak, the effusion rate declined rapidly to 10–20 m³ s^–1, then more slowly to 1 m³ s^–1, and remained at 1–12 m³ s^–1 until the end of the eruption. Site observations near the main crater suggest that the intensity of the volcanic tremor varied directly with the force of the eruption. A notable rise in the fluorine concentration of riverwater in the vicinity of the eruptive fissures occurred on the 5th day of the eruption, but it levelled off on the 6th day and then remained essentially constant. The volume and initial silica content of the lava and tephra, the explosivity and effusion rate during the earliest stage of the eruption, as well as the magnitude attained by the associated earthquakes, support earlier suggestions that these parameters are positively related to the length of the preceeding repose period. The chemical difference between the eruptive material of Hekla itself and the lavas erupted in its vicinity can be explained in terms of a density-stratified magma reservoir located at the bottom of the crust. We propose that the shape of this reservoir, its location at the west margin of a propagating rift, and its association with a crustal weakness, all contribute to the high eruption frequency of Hekla.     

The caldera-forming eruption of Ksudach volcano about cal. A.D. 240: the greatest explosive event of our era in Kamchatka, Russia

The largest Plinian eruption of our era and the latest caldera-forming eruption in the Kuril-Kamchatka region occurred about cal. A.D. 240 from the Ksudach volcano. This catastrophic explosive eruption was similar in type and characteristics to the 1883 Krakatau event. The volume of material ejected was 18–19 km³ (8 km³ DRE), including 15 km³ of tephra fall and 3–4 km³ of pyroclastic flows. The estimated height of eruptive column is 22–30 km. A collapse caldera resulting from this eruption was 4 × 6.5 km in size with a cavity volume of 6.5–7 km³. Tephra fall was deposited to the north of the volcano and reached more than 1000 km. Pyroclastic flows accompanied by ash-cloud pyroclastic surges extended out to 20 km. The eruption was initially phreatomagmatic and then became rhythmic, with each pulse evolving from pumice falls to pyroclastic flows. Erupted products were dominantly rhyodacite throughout the eruption. During the post-caldera stage, when the Shtyubel cone started to form within the caldera, basaltic-andesite and andesite magma began to effuse. The trigger for the eruption may have been an intrusion of mafic magma into the rhyodacite reservoir. The eruption had substantial environmental impact and may have produced a large acidity peak in the Greenland ice sheet.     

Changes in diatom, pollen, and chironomid assemblages in response to a recent volcanic event in Lake Galletué (Chilean Andes)

Several lakes in Chile are near important volcanic areas where eruption impacts can limit the quality of lacustrine sediments for reconstructing past environmental changes. In this study, we report changes in diatoms, pollen, and chironomids assemblages after a tephra deposition in Lake Galletué (Chilean Andes). A sediment core obtained from Lake Galletué (40 m water depth) was sliced in 1 cm intervals and subsamples were taken to analyze each proxy. ²¹⁰Pb and ¹³⁷Cs activities were measured to obtain the geochronology and mineralogical analyses were performed to determine the mineral composition of the tephra. Diatom species composition and productivity were modified when the lake received the tephra; Aulacoseira granulata decreased and was later replaced by Cyclotella af. glomerata. After the tephra input, Aulacoseira granulata abundance increased to pre-disturbance levels and Cyclotella af. glomerata decreased. These changes seem to suggest a momentary increase in lake nutrient levels after the tephra deposition. Chironomid assemblages also decreased in head capsules just after the tephra deposition, but the most important change was the replacement of Ablabesmyia by Parakiefferiella, probably due to the sedimentological changes produced by the input of coarse tephra grains. Furthermore, unlike other studies, chironomid assemblages in Lake Galletué did not show a decrease drastically in diversity within the tephra layer. The pollen analysis indicated that, prior to the volcanic event, the vegetal community was dominated by Nothofagus sp., Araucaria araucana, and Blechnum sp.-type. After the tephra deposition, the same taxa are dominant, indicating that the volcanic event seems not produce changes in the vegetation. Nevertheless, within the tephra layer it is possible to see an increase in Poaceae, which represent – due to the percolation process – the effect of eruption on the vegetation. According to our results, diatoms were the most sensitive proxy for describing the changes produced by tephra deposition into the aquatic ecosystem and, despite the noticeable changes in its sedimentological properties; the lake seems to have a high resilience capacity, allowing it to return to pre-tephra input conditions.     

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.     

The May 1915 eruptions of Lassen Peak, II: May 22 volcanic blast effects, sedimentology and stratigraphy of deposits, and characteristics of the blast cloud

The May 22, 1915 eruptions of Lassen Peak involved a volcanic blast and the emplacement of three geographically and temporally distinct lahar deposits. The volcanic blast occurred when a Vulcanian explosion at the summit unroofed a shallow magma source, generating an eruption cloud that rose to an estimated height of 9 km above sea level. The blast cloud was probably caused by the collapse of a small portion of the eruption column; absence of a flank vent associated with these eruptions argues against it originating as an explosion that has been directed by vent geometry or location. The volcanic blast devasted 7 km² of the northeast flank of the volcano, and emplaced a deposit of juvenile tephra and accidental lithic and mineral fragments. Decrease in blast deposit thickness and median grain size with increasing distance from the vent suggests that the blast cloud lost transport competence as it crossed the devastated area. Scanning electron microscope examination of pyroclasts from the blast deposit indicates that the blast cloud was a dry, turbulent suspension that emplaced a thin deposit which cooled rapidly after deposition. Lahar deposits were emplaced primarily in Lost Creek, with minor lahars flowing down gullies on the west, northwest and north flanks of the volcano. The initial lahar was apparently triggered early in the eruption when the blast cloud melted the residual snowpack as it moved down the northeast flank of the peak. The event that triggered the later lahars is enigmatic; the presence of approximately five times more juvenile dacite bombs on the surface of the later lahars suggests that they may have been triggered by a change in eruption style or dynamics.     

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.     

Possible effects of volcanic eruptions on stratospheric minor constituent chemistry

Although stratosphere penetrating volcanic eruptions have been infrequent during the last half century, periods have existed in the last several hundred years when such eruptions were significantly more frequent. Several mechanisms exist for these injections to affect stratospheric minor constitutent chemistry, both on the long-term average and for short-term perturbations. These mechanisms are reviewed and, because of the sensitivity of current models of stratospheric ozone to chlorine perturbations, quantitative estimates are made of chlorine injection rates. It is found that, if chlorine makes up as much as 0.5 to 1% of the gases released and if the total gases released are about the same magnitude as the fine ash, then a major stratosphere penetrating eruption could deplete the ozone column by several percent. The estimate for the Agung eruption of 1963 is just under 1% an amount not excluded by the ozone record but complicated by the peak in atmospheric nuclear explosions at about the same time. The long-term contribution to stratospheric CIX by volcanic eruptions is estimated as 0.1 ppbv for the period 1900–60 and 1 ppbv for the much more volcanically active period 1780–1840. All of the estimates are subject to large uncertainties, perhaps a factor of 2 or 3 on the high side and a factor of 10 or more on the low side.Paper presented at the IAGA/IAMAP Joint Assembly, Seattle, WA, U.S.A., August 1977.