Stratigraphy and eruptive dynamics of a pulsating Plinian eruption of Somma-Vesuvius: the Pomici di Mercato (8900 years B.P.)

New volcanological studies allow reconstruction of the eruption dynamics of the Pomici di Mercato eruption (ca 8,900 cal. yr B.P.) of Somma-Vesuvius. Three main Eruptive Phases are distinguished based on two distinct erosion surfaces that interrupt stratigraphic continuity of the deposits, indicating that time breaks occurred during the eruption. Absence of reworked volcaniclastic deposits on top of the erosion surfaces suggests that quiescent periods between eruptive phases were short perhaps lasting only days to weeks. Each of the Eruptive Phases was characterised by deposition of alternating fall and pyroclastic density current (PDC) deposits. The fallout deposits blanketed a wide area toward the east, while the more restricted PDC deposits inundated the volcano slopes. Eruptive dynamics were driven by brittle magmatic fragmentation of a phonolitic magma, which, because of its mechanical fragility, produced a significant amount of fine ash. External water did not significantly contribute either to fragmentation dynamics or to mechanical energy release during the eruption. Column heights were between 18 and 22 km, corresponding to mass discharge rates between 1.4 and 6 × 10⁷ kg s⁻¹. The estimated on land volume of fall deposits ranges from a minimum of 2.3 km³ to a maximum of 7.4 km³. Calculation of physical parameters of the dilute pyroclastic density currents indicates speeds of a few tens of m s⁻¹ and densities of a few kg m⁻³ (average of the lowermost 10 m of the currents), resulting in dynamic pressures lower than 3 kPa. These data suggest that the potential impact of pyroclastic density currents of the Pomici di Mercato eruption was smaller than those of other Plinian and sub-Plinian eruptions of Somma-Vesuvius, especially those of 1631 AD and 472 AD (4–14 kPa), which represent reference values for the Vesuvian emergency plan. The pulsating and long-lasting behaviour of the Pomici di Mercato eruption is unique in the history of large explosive eruptions of Somma-Vesuvius. We suggest an eruptive scheme in which discrete magma batches rose from the magma chamber through a network of fractures. The injection and rise of the different magma batches was controlled by the interplay between magma chamber overpressure and local stress. The intermittent discharge of magma during a large explosive eruption is unusual for Somma-Vesuvius, as well as for other volcanoes worldwide, and yields new insights for improving our knowledge of the dynamics of explosive eruptions.     

The Pomici di Avellino eruption of Somma-Vesuvius (3.9 ka bp). Part I: stratigraphy, compositional variability and eruptive dynamics

The stratigraphic succession of the Pomici di Avellino Plinian eruption from Somma-Vesuvius has been studied through field and laboratory data in order to reconstruct the eruption dynamics. This eruption is particularly important in the Somma-Vesuvius eruptive history because (1) its vent was offset with respect to the present day Vesuvius cone; (2) it was characterised by a distinct opening phase; (3) breccia-like very proximal fall deposits are preserved close to the vent and (4) the pyroclastic density currents generated during the final phreatomagmatic phase are among the most widespread and voluminous in the entire history of the volcano. The stratigraphic succession is, here, divided into deposits of three main eruptive phases (opening, magmatic Plinian and phreatomagmatic), which contain five eruption units. Short-lived sustained columns occurred twice during the opening phase (H_t of 13 and 21.5 km, respectively) and dispersed thin fall deposits and small pyroclastic density currents onto the volcano slopes. The magmatic Plinian phase produced the main volume of erupted deposits, emplacing white and grey fall deposits which were dispersed to the northeast. Peak column heights reached 23 and 31 km during the withdrawal of the white and the grey magmas, respectively. Only one small pyroclastic density current was emplaced during the main Plinian phase. In contrast, the final phreatomagmatic phase was characterised by extensive generation of pyroclastic density currents, with fallout deposits very subordinate and limited to the volcano slopes. Assessed bulk erupted volumes are 21 × 10⁶ m³ for the opening phase, 1.3–1.5 km³ for the main Plinian phase and about 1 km³ for the final phreatomagmatic phase, yielding a total volume of about 2.5 km³. Pumice fragments are porphyritic with sanidine and clinopyroxene as the main mineral phases but also contain peculiar mineral phases like scapolite, nepheline and garnet. Bulk composition varies from phonolite (white magma) to tephri-phonolite (grey magma).     

The Pomici di Avellino eruption of Somma–Vesuvius (3.9 ka BP). Part II: sedimentology and physical volcanology of pyroclastic density current deposits

Pyroclastic density currents (PDCs) generated during the Plinian eruption of the Pomici di Avellino (PdA) of Somma–Vesuvius were investigated through field and laboratory studies, which allowed the detailed reconstruction of their eruptive and transportation dynamics and the calculation of key physical parameters of the currents. PDCs were generated during all the three phases that characterised the eruption, with eruptive dynamics driven by both magmatic and phreatomagmatic fragmentation. Flows generated during phases 1 and 2 (EU1 and EU3pf, magmatic fragmentation) have small dispersal areas and affected only part of the volcano slopes. Lithofacies analysis demonstrates that the flow-boundary zones were dominated by granular-flow regimes, which sometimes show transitions to traction regimes. PDCs generated during eruptive phase 3 (EU5, phreatomagmatic fragmentation) were the most voluminous and widespread in the whole of Somma–Vesuvius’ eruptive history, and affected a wide area around the volcano with deposit thicknesses of a few centimetres up to more than 25 km from source. Lithofacies analysis shows that the flow-boundary zones of EU5 PDCs were dominated by granular flows and traction regimes. Deposits of EU5 PDC show strong lithofacies variation northwards, from proximally thick, massive to stratified beds towards dominantly alternating beds of coarse and fine ash in distal reaches. The EU5 lithofacies also show strong lateral variability in proximal areas, passing from the western and northern to the eastern and southern volcano slopes, where the deposits are stacked beds of massive, accretionary lapilli-bearing fine ash. The sedimentological model developed for the PDCs of the PdA eruption explains these strong lithofacies variations in the light of the volcano’s morphology at the time of the eruption. In particular, the EU5 PDCs survived to pass over the break in slope between the volcano sides and the surrounding volcaniclastic apron–alluvial plain, with development of new flows from the previously suspended load. Pulses were developed within individual currents, leading to stepwise deposition on both the volcano slopes and the surrounding volcaniclastic apron and alluvial plain. Physical parameters including velocity, density and concentration profile with height were calculated for a flow of the phreatomagmatic phase of the eruption by applying a sedimentological method, and the values of the dynamic pressure were derived. Some hazard considerations are summarised on the assumption that, although not very probable, similar PDCs could develop during future eruptions of Somma–Vesuvius.     

The 1919–1920 eruption of Mauna Iki,Kilauea: chronology,geologic mapping,and magma transport mechanisms

Utilizing historical accounts, field mapping, and photogeology, this paper presents a chronology of, and an analysis of magma transport during, the December 1919 to August 1920 satellitic shield eruption of Mauna Iki on the SW rift zone of Kilauea Volcano, Hawaii. The eruption can be divided into four stages based on the nature of the eruptive activity. Stage 1 consisted of the shallow injection of a dike from the summit region to the eventual eruption site 10 km downrift. During stage 2, a low ridge of pahoehoe formed in the vent area; later a large a'a flow broke out of this ridge and flowed 8.5 km SW at an average flow front velocity of 0.5 km/day. The eruption continued until mid-August producing almost exclusively pahoehoe, first as gas-rich overflows from a lava pond (stage 3), and later as denser tube-fed lava (stage 4) that reached almost 8 km from the vent at an average flow-front velocity of 0.1 km/day. Magma transport during the Mauna Iki eruption is examined using three criteria: (1) eruption characteristics and volumetric flow rates; (2) changes in the surface height of the Halemaumau lava lake; and (3) tilt measurements made at the summit of Kilauea. We find good correlation between Halemaumau lake activity and the eruptive stages. Additionally, the E-W component of summit tilt tended to mimic the lake activity. The N-S component, however, did not. Multiple storage zones in the shallow summit region probably accounted for the decoupling of E-W and N-S tilt components. Analysis of these criteria shows that at different times during the eruption, magma was either emplaced into the volcano without eruption, hydraulically drained from Halemaumau to Mauna Iki, or fed at steady-state conditions from summit storage to Mauna Iki. Volume calculations indicate that the supply rate to Kilauea during the eruption was around 3 m³/s, similar to that calculated during the Mauna Ulu and Kupaianaha shield-building eruptions, and consistent with previously determined values of long-term supply to Kilauea.     

Volcanological evolution of the Rivi–Capo Volcanic Complex at Salina,Aeolian Islands: magma storage processes and ascent dynamics

Critical to understanding explosive eruptions is establishing how accurately representative pyroclasts are of processes during magma vesiculation and fragmentation. Here, we present data on densities, and vesicle size and number characteristics, for representative pyroclasts from six silicic eruptions of contrasting size and style from Raoul volcano (Kermadec arc). We use these data to evaluate histories of bubble nucleation, coalescence, and growth in explosive eruptions and to provide comparisons with pumiceous dome carapace material. Density/vesicularity distributions show a scarcity of pyroclasts with ～65–75 % vesicularity; however, pyroclasts closest to this vesicularity range have the highest bubble number density (BND) values regardless of eruptive intensity or style. Clasts with vesicularities greater than this 65–75 % “pivotal” vesicularity range have decreasing BNDs with increasing vesicularities, interpreted to reflect continuing bubble growth and coalescence. Clasts with vesicularities less than the pivotal range have BNDs that decrease with decreasing vesicularity and preserve textures indicative of processes such as stalling and open system degassing prior to vesiculation in a microlite-rich magma, or vesiculation during slow ascent of degassing magma. Bubble size distributions (BSDs) and BNDs show variations consistent with 65–75 % representing the vesicularity at which vesiculating magma is most likely to undergo fragmentation, consistent with the closest packing of spheres. We consider that the observed vesicularity range may reflect the development of permeability in the magma through shearing as it flows through the conduit. These processes can act in concert with multiple nucleation events, generating a situation of heterogeneous bubble populations that permit some regions of the magma to expand and bubbles to coalesce with other regions in which permeable networks are formed. Fragmentation preserves the range in vesicularity seen as well as any post-fragmentation/pre-quenching expansion which may have occurred. We demonstrate that differing density pyroclasts from a single eruption interval can have widely varying BND values corresponding to the degree of bubble maturation that has occurred. The modal density clasts (the usual targets for vesicularity studies) have likely undergone some degree of bubble maturation and are therefore may not be representative of the magma at the onset of fragmentation.     

The Holocene Pucón eruption of Volcán Villarrica, Chile: deposit architecture and eruption chronology

The Pucón eruption was the largest Holocene explosive outburst of Volcán Villarrica, Chile. It discharged >1.0 km³ of basaltic-andesite magma and >0.8 km³ of pre-existing rock, forming a thin scoria-fall deposit overlain by voluminous ignimbrite intercalated with pyroclastic surge beds. The deposits are up to 70 m thick and are preserved up to 21 km from the present-day summit, post-eruptive lahar deposits extending farther. Two ignimbrite units are distinguished: a lower one (P1) in which all accidental lithic clasts are of volcanic origin and an upper unit (P2) in which basement granitoids also occur, both as free clasts and as xenoliths in scoria. P2 accounts for ∼80% of the erupted products. Following the initial scoria fallout phase, P1 pyroclastic flows swept down the northern and western flanks of the volcano, magma fragmentation during this phase being confined to within the volcanic edifice. Following a pause of at least a couple of days sufficient for wood devolatilization, eruption recommenced, the fragmentation level dropped to within the granitoid basement, and the pyroclastic flows of P2 were erupted. The first P2 flow had a highly turbulent front, laid down ignimbrite with large-scale cross-stratification and regressive bedforms, and sheared the ground; flow then waned and became confined to the southeastern flank. Following emplacement of pyroclastic surge deposits all across the volcano, the eruption terminated with pyroclastic flows down the northern flank. Multiple lahars were generated prior to the onset of a new eruptive cycle. Charcoal samples yield a probable eruption age of 3,510 ± 60 ¹⁴C years BP.     

The November 2002 eruption of Piton de la Fournaise,Réunion: tracking the pre-eruptive thermal evolution of magma using melt inclusions

The November 2002 eruption of Piton de la Fournaise in the Indian Ocean was typical of the activity of the volcano from 1999 to 2006 in terms of duration and volume of magma ejected. The first magma erupted was a basaltic liquid with a small proportion of olivine phenocrysts (Fo₈₁) that contain small numbers of melt inclusions. In subsequent flows, olivine crystals were more abundant and richer in Mg (Fo_83–84). These crystals contain numerous melt and fluid inclusions, healed fractures, and dislocation features such as kink bands. The major element composition of melt inclusions in this later olivine (Fo_83–84) is out of equilibrium with that of its host as a result of extensive post-entrapment crystallization and Fe²⁺ loss by diffusion during cooling. Melt inclusions in Fo₈₁ olivine are also chemically out of equilibrium with their hosts but to a lesser degree. Using olivine–melt geothermometry, we determined that melt inclusions in Fo₈₁ olivine were trapped at lower temperature (1,182 ± 1°C) than inclusions in Fo_83–84 olivine (1,199–1,227°C). This methodology was also used to estimate eruption temperatures. The November 2002 melt inclusion compositions suggest that they were at temperatures between 1,070°C and 1,133°C immediately before eruption and quenching. This relatively wide temperature range may reflect the fact that most of the melt inclusions were from olivine in lava samples and therefore likely underwent minor but variable amounts of post-eruptive crystallization and Fe²⁺ loss by diffusion due to their relatively slow cooling on the surface. In contrast, melt inclusions in tephra samples from past major eruptions yielded a narrower range of higher eruption temperatures (1,163–1,181°C). The melt inclusion data presented here and in earlier publications are consistent with a model of magma recharge from depth during major eruptions, followed by storage, cooling, and crystallization at shallow levels prior to expulsion during events similar in magnitude to the relatively small November 2002 eruption.     

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.     

Insights into eruption dynamics from textural analysis: the case of the May, 2008, Chaitén eruption

The May, 2008, Chaitén (southern Chile) eruption was characterized by several explosive events, each associated with plumes which reached up to about 19?km above sea level on May 6. A study of the textural and physical features of the juvenile clasts erupted during the climactic phase of the 2008 eruption of Chaitén is presented. Pumice clasts show unimodal density distribution (main mode at 600?kg/m³), average vesicularity of about 69?%, a glassy groundmass with no microcrystals, and vesicles with dimension between ～1?μm and ～2?mm. They also show a unimodal vesicle size distribution with most frequent vesicle size in the range 0.05–0.08?mm and an estimated vesicle number density of 1.3?±?0.5?×?10⁵?mm^?3 related to a rapid nucleation event produced during the late phases of magma rise. This is confirmed by the absence of microcrystals that could otherwise have delayed vesicle formation and allowed the magma to maintain a low viscosity and a supersaturation in volatiles. Vesiculation and fragmentation were triggered by a sudden decompression of the melt associated with the opening of the volcanic conduit (～10?MPa?s^?1).     

Melt inclusion analysis of the Unzen 1991–1995 dacite: implications for crystallization processes of dacite magma

Dacitic magma, a mixture of high-temperature (T) aphyric magma and low-T crystal-rich magma, was erupted during the 1991–1995 Mount Unzen eruptive cycle. Here, the crystallization processes of the low-T magma were examined on the basis of melt inclusion analysis and phase relationships. Variation in water content of the melt inclusions (5.1–7.2 wt% H₂O) reflected the degassing history of the low-T magma ascending from deeper levels (250 MPa) to a shallow magma chamber (140 MPa). The ascent rate of the low-T magma decreased markedly towards the emplacement level as crystal content increased. Cooling of magma as well as degassing-induced undercooling drove crystallization. With the decreasing ascent rate, degassing-induced undercooling decreased in importance, and cooling became more instrumental in crystallization, causing local and rapid crystallization along the margin of the magma body. Some crystals contain scores of melt inclusions, whereas there are some crystals without any inclusions. This heterogeneous distribution suggests the variation in the crystallization rate within the magma body; it also suggests that cooling was dominant cause for melt entrapment. Numerical calculations of the cooling magma body suggest that cooling caused rapid crystal growth and enhanced melt entrapment once the magma became a crystal-rich mush with evolved interstitial melt. The rhyolitic composition of melt inclusions is consistent with this model.Editorial responsibility: H Shinohara     

The settling dynamics of flocculating mud and sand mixtures: part 2—numerical modelling

Estuarine and coastal sediment transport is characterised by the transport of both sand-sized particles (of diameter greater than 63?μm) and muddy fine-grained sediments (silt, diameter less than 63?μm; clay, diameter less than 2?μm). These fractions are traditionally considered as non-cohesive and cohesive, respectively, because of the negligible physico-chemical attraction that occurs between sand grains. However, the flocculation of sediment particles is not only caused by physico-chemical attraction. Cohesivity of sediment is also caused by biology, in particular the sticky extra-cellular polymeric substances secreted by diatoms, and the effect of biology in binding sediment particles can be much larger than that of physico-chemical attraction. As demonstrated by Manning (2008) and further expanded in part 1 of this paper (Manning et al., submitted), the greater binding effect of biology allows sand particles to flocculate with mud. In many estuaries, both the sand and fine sediment fractions are transported in significant quantities. Many of the more common sediment transport modelling suites now have the capability to combine mud and sand transport. However, in all of these modelling approaches, the modelling of mixed sediment transport has still essentially separated the modelling of sand and mud fractions assuming that these different fractions do not interact except at the bed. However, the use of in situ video techniques has greatly enhanced the accuracy and reliability of settling velocity measurements and has led to a re-appraisal of this widely held assumption. Measurements of settling velocity in mixed sands presented by Manning et al. (2009) have shown strong evidence for the flocculation of mixed sediments, whilst the greater understanding of the role of biology in flocculation has identified mechanisms by which this mud-sand flocculation can occur. In the first part of this paper (Manning et al., submitted), the development of an empirical flocculation model is described which represents the interaction between sand and mud particles in the flocculation process. Measurements of the settling velocity of varying mud-sand mixtures are described, and empirical algorithms governing the variation of settling velocity with turbulence, suspended sediment concentration and mud-sand content are derived. The second part of this paper continues the theme of examination of the effects of mud-sand interaction on flocculation. A 1DV mixed transport model is developed and used to reproduce the vertical transport of mixed sediment fractions. The 1DV model is used to reproduce the measured settling velocities in the laboratory experiments described in the part 1 paper and also to reproduce measurements of concentration of mixed sediments in the Outer Thames. In both modelling exercises, the model is run using the algorithms developed in part 1 and repeated using an assumption of no interaction between mud and sand in the flocculation process. The results of the modelling show a significant improvement in the ability of the 1DV to reproduce the observed sediment behaviour when the empirical equations are used. This represents further strong evidence of the interaction between sand and mud in the flocculation process.     

Eruption dynamics of Hawaiian-style fountains: the case study of episode 1 of the Kīlauea Iki 1959 eruption

Hawaiian eruptions are characterized by fountains of gas and ejecta, sustained for hours to days that reach tens to hundreds of meters in height. Quantitative analysis of the pyroclastic products from the 1959 eruption of Kīlauea Iki, Kīlauea volcano, Hawai‘i, provides insights into the processes occurring during typical Hawaiian fountaining activity. This short-lived but powerful eruption contained 17 fountaining episodes and produced a cone and tephra blanket as well as a lava lake that interacted with the vent and fountain during all but the first episode of the eruption, the focus of this paper. Microtextural analysis of Hawaiian fountaining products from this opening episode is used to infer vesiculation processes within the fountain and shallow conduit. Vesicle number densities for all clasts are high (10⁶–10⁷ cm⁻³). Post-fragmentation expansion of bubbles within the thermally-insulated fountain overprints the pre-fragmentation bubble populations, leading to a reduction in vesicle number density and increase in mean vesicle size. However, early quenched rims of some clasts, with vesicle number densities approaching 10⁷ cm⁻³, are probably a valid approximation to magma conditions near fragmentation. The extent of clast evolution from low vesicle-to-melt ratio and corresponding high vesicle number density to higher vesicle-to-melt ratio and lower vesicle-number density corresponds to the length of residence time within the fountain.     

The May 2005 eruption of Fernandina volcano,Galápagos: The first circumferential dike intrusion observed by GPS and InSAR

The May 2005 eruption of Fernandina volcano, Galápagos, occurred along circumferential fissures parallel to the caldera rim and fed lava flows down the steep southwestern slope of the volcano for several weeks. This was the first circumferential dike intrusion ever observed by both InSAR and GPS measurements and thus provides an opportunity to determine the subsurface geometry of these enigmatic structures that are common on Galápagos volcanoes but are rare elsewhere. Pre- and post- eruption ground deformation between 2002 and 2006 can be modeled by the inflation of two separate magma reservoirs beneath the caldera: a shallow sill at ~1 km depth and a deeper point-source at ~5 km depth, and we infer that this system also existed at the time of the 2005 eruption. The co-eruption deformation is dominated by uplift near the 2005 eruptive fissures, superimposed on a broad subsidence centered on the caldera. Modeling of the co-eruption deformation was performed by including various combinations of planar dislocations to simulate the 2005 circumferential dike intrusion. We found that a single planar dike could not match both the InSAR and GPS data. Our best-fit model includes three planar dikes connected along hinge lines to simulate a curved concave shell that is steeply dipping (~45–60°) toward the caldera at the surface and more gently dipping (~12–14°) at depth where it connects to the horizontal sub-caldera sill. The shallow sill is underlain by the deep point source. The geometry of this modeled magmatic system is consistent with the petrology of Fernandina lavas, which suggest that circumferential eruptions tap the shallowest parts of the system, whereas radial eruptions are fed from deeper levels. The recent history of eruptions at Fernandina is also consistent with the idea that circumferential and radial intrusions are sometimes in a stress-feedback relationship and alternate in time with one another.     

Complex dynamics of small-moderate volcanic events: the example of the 2011 rhyolitic Cordón Caulle eruption,Chile

After decades of repose, Puyehue-Cordón Caulle Volcano (Chile) erupted in June 2011 following a month of continuously increasing seismic activity. The eruption dispersed a large volume of rhyolitic tephra over a wide area and was characterized by complex dynamics. During the initial climactic phase of the eruption (24–30 h on 4–5 June), 11–14-km-high plumes dispersed most of the erupted tephra eastward towards Argentina, reaching as far as the Atlantic Ocean. This first eruptive phase was followed by activity of lower intensity, leading to the development of a complex stratigraphic sequence, mainly due to rapid shifts in wind direction and eruptive style. The resulting tephra deposits consist of 13 main layers grouped into four units. Each layer was characterized based on its dispersal direction, sedimentological features, and on the main characteristics of the juvenile fraction (texture, density, petrography, chemistry). The lowest part of the eruptive sequence (Unit I), corresponding to the tephra emitted between 4 and 5 June, is composed of alternating lapilli layers with a total estimated volume of ca. 0.75 km³; these layers record the highest intensity phase, during which a bent-over plume dispersed tephra towards the southeast-east, with negligible up-wind sedimentation. Products emitted during 5–6 June (Unit II) signaled an abrupt shift in wind direction towards the north, leading to the deposition of a coarse ash deposit in the northern sector (ca. 0.21 km³ in volume), followed by a resumption of easterly directed winds. A third phase (Unit III) began on 7 June and resulted in tephra deposits in the eastern sector and ballistic bombs around the vent area. A final phase (Unit IV) started after 15 June and was characterized by the emission of fine-grained white tephra from ash-charged plumes during low-level activity and the extrusion of a viscous lava flow. Timing and duration of the first eruptive phases were constrained based on comparison of the dispersal of the main tephra layers with satellite images, showing that most of the tephra was emitted during the first 72 h of the event. The analyzed juvenile material tightly clusters within the rhyolitic field, with negligible chemical variations through the eruptive sequence. Textural observations reveal that changes in eruption intensity (and consequently in magma ascent velocity within the conduit) and complex interactions between gas-rich and gas-depleted magma portions during ascent resulted in vesicular clasts with variable degrees of shear localization, and possibly in the large heterogeneity of the juvenile material.     

Estimation of the heat stored by the magma chamber of El’brus Volcano in the host rocks and how it can be extracted

The results of geological and geophysical studies are analyzed to throw light on the presence of an unsolidified magma chamber beneath El’brus Volcano in the Caucasus, as well as its depth and approximate dimensions. We provide an upper and a lower bound on the estimates of the heat stored by the chamber in the host rocks, which were heated by the volcano’s magma chamber from its origination until now, incorporating the variation in the dimensions of the magma chamber during its evolution and the heat storage in it. We examine the geological and geophysical conditions that favor the use of thermal energy in the hot rocks that surround the magma chamber of El’brus Volcano.