Relationships between deformation and mixing processes in lava flows: a case study from Salina (Aeolian Islands,Tyrrhenian Sea)

 The massive unit of a lava flow from Porri volcano (Salina, Aeolian Islands) displays many unusual structures related to the physical interaction between two different magmas. The magma A represents approximately 80% of the exposed lava surface; it has a crystal content of 51 vol.% and a dacitic glass composition (SiO₂=63–64 wt.%). The magma B has a basaltic-andesite glass composition (SiO₂=54–55 wt.%) and a crystal content of approximately 18 vol.%. It occurs as pillow-like enclaves, banding, boudin-like and rolling structures which are hosted in magma A. Structural analysis suggests that banding and boudin-like structures are the result of the deformation of enclaves at different shear strain. The linear correlation between strain and stratigraphic height of the measured elements indicates a single mode of deformation. We deduce that the component B deformed according to a simple shear model. Glass analyses of the A–B boundary indicate that A and B liquids mix together at high shear strain, whereas only mingling occurs at low shear strain. This suggests that the amount of deformation (i.e. forced convection) plays an important role in the formation of hybrid magmas. High shear strain may induce stretching, shearing and rolling of fluids which promote both forced convection and dynamical diffusion processes. These processes allow mixing of magmas with large differences in their physical properties. Received: 15 July 1995 / Accepted: 30 May 1996     

Cooling and crystallization of lava in open channels, and the transition of Pāhoehoe Lava to 'A'ā

 Samples collected from a lava channel active at Kīlauea Volcano during May 1997 are used to constrain rates of lava cooling and crystallization during early stages of flow. Lava erupted at near-liquidus temperatures (∼1150  °C) cooled and crystallized rapidly in upper parts of the channel. Glass geothermometry indicates cooling by 12–14  °C over the first 2 km of transport. At flow velocities of 1–2 m/s, this translates to cooling rates of 22–50  °C/h. Cooling rates this high can be explained by radiative cooling of a well-stirred flow, consistent with observations of non-steady flow in proximal regions of the channel. Crystallization of plagioclase and pyroxene microlites occurred in response to cooling, with crystallization rates of 20–50% per hour. Crystallization proceeded primarily by nucleation of new crystals, and nucleation rates of ∼10⁴/cm³s are similar to those measured in the 1984 open channel flow from Mauna Loa Volcano. There is no evidence for the large nucleation delays commonly assumed for plagioclase crystallization in basaltic melts, possibly a reflection of enhanced nucleation due to stirring of the flow. The transition of the flow surface morphology from pāhoehoe to 'a'ā occurred at a distance of 1.9 km from the vent. At this point, the flow was thermally stratified, with an interior temperature of ∼1137  °C and crystallinity of ∼15%, and a flow surface temperature of ∼1100  °C and crystallinity of ∼45%. 'A'ā formation initiated along channel margins, where crust was continuously disrupted, and involved tearing and clotting of the flow surface. Both observations suggest that the transition involved crossing of a rheological threshold. We suggest this threshold to be the development of a lava yield strength sufficient to prevent viscous flow of lava at the channel margin. We use this concept to propose that 'a'ā formation in open channels requires both sufficiently high strain rates for continued disruption of surface crusts and sufficient groundmass crystallinity to generate a yield strength equivalent to the imposed stress. In Hawai'i, where lava is typically microlite poor on eruption, these combined requirements help to explain two common observations on 'a'ā formation: (a) 'a'ā flow fields are generated when effusion rates are high (thus promoting crustal disruption); and (b) under most eruption conditions, lava issues from the vent as pāhoehoe and changes to 'a'ā only after flowing some distance, thus permitting sufficient crystallization. Received: 3 September 1998 / Accepted: 12 April 1999     

Degassing during magma ascent in the Mule Creek vent (USA)

 The structures and textures of the rhyolite in the Mule Creek vent (New Mexico, USA) indicate mechanisms by which volatiles escape from silicic magma during eruption. The vent outcrop is a 300-m-high canyon wall comprising a section through the top of a feeder conduit, vent and the base of an extrusive lava dome. Field relations show that eruption began with an explosive phase and ended with lava extrusion. Analyses of glass inclusions in quartz phenocrysts from the lava indicate that the magma had a pre-eruptive dissolved water content of 2.5–3.0 wt% and, during eruption, the magma would have been water-saturated over the vertical extent of the present outcrop. However, the vesicularity of the rhyolite is substantially lower than that predicted from closed-system models of vesiculation under equilibrium conditions. At a given elevation in the vent, the volume fraction of primary vesicles in the rhyolite increases from zero close to the vent margin to values of 20–40 vol.% in the central part. In the centre the vesicularity increases upward from approximately 20 vol.% at 300 m below the canyon rim to approximately 40 vol.% at 200 m, above which it shows little increase. To account for the discrepancy between observed vesicularity and measured water content, we conclude that gas escaped during ascent, probably beginning at depths greater than exposed, by flow through the vesicular magma. Gas escape was most efficient near the vent margin, and we postulate that this is due both to the slow ascent of magma there, giving the most time for gas to escape, and to shear, favouring bubble coalescence. Such shear-related permeability in erupting magma is supported by the preserved distribution of textures and vesicularity in the rhyolite: Vesicles are flattened and overlapping near the dense margins and become progressively more isolated and less deformed toward the porous centre. Local zones have textures which suggest the coalescence of bubbles to form permeable, collapsing foams, implying the former existence of channels for gas migration. Local channelling of gas into the country rocks is suggested by the presence of sub-horizontal syn-eruptive rhyolitic tuffisite veins which depart from the vent margin and invade the adjacent country rock. In the central part of the vent, similar local channelling of gas is indicated by steep syn-eruption tuffisite veins which cut the rhyolite itself. We conclude that the suppression of explosive eruption resulted from gas separation from the ascending magma and vent structure by shear-related porous flow and channelling of gas through tuffisite veins. These mechanisms of gas loss may be responsible for the commonly observed transition from explosive to effusive behaviour during the eruption of silicic magma. Received: 24 May 1995 / Accepted: 13 March 1996     

The changing morphology of an open lava channel on Mt. Etna

An open channel lava flow on Mt. Etna (Sicily) was observed during May 30–31, 2001. Data collected using a forward looking infrared (FLIR) thermal camera and a Minolta-Land Cyclops 300 thermal infrared thermometer showed that the bulk volume flux of lava flowing in the channel varied greatly over time. Cyclic changes in the channel's volumetric flow rate occurred over several hours, with cycle durations of 113–190 min, and discharges peaking at 0.7 m³ s⁻¹ and waning to 0.1 m³ s⁻¹. Each cycle was characterized by a relatively short, high-volume flux phase during which a pulse of lava, with a well-defined flow front, would propagate down-channel, followed by a period of waning flow during which volume flux lowered. Pulses involved lava moving at relatively high velocities (up to 0.29 m s⁻¹) and were related to some change in the flow conditions occurring up-channel, possibly at the vent. They implied either a change in the dense rock effusion rate at the source vent and/or cyclic-variation in the vesicle content of the lava changing its bulk volume flux. Pulses would generally overspill the channel to emplace pāhoehoe overflows. During periods of waning flow, velocities fell to 0.05 m s^–1. Blockages forming during such phases caused lava to back up. Occasionally backup resulted in overflows of slow moving ‘a‘ā that would advance a few tens of meters down the levee flank. Compound levees were thus a symptom of unsteady flow, where overflow levees were emplaced as relatively fast moving pāhoehoe sheets during pulses, and as slow-moving ‘a‘ā units during backup. Small, localized fluctuations in channel volume flux also occurred on timescales of minutes. Volumes of lava backed up behind blockages that formed at constrictions in the channel. Blockage collapse and/or enhanced flow under/around the blockage would then feed short-lived, wave-like, down-channel surges. Real fluctuations in channel volume flux, due to pulses and surges, can lead to significant errors in effusion rate calculations. Editorial responsibility: A. Woods     

Cone-building block-and-ash flows: the Senyama volcanic products of O’e Takayama volcano,SW Japan

The Senyama volcanic products of the late Pliocene to early Pleistocene O’e Takayama volcano overlie a 100-m-thick, late Pliocene coastal quartz-sandstone and are intruded by an early Pleistocene dacite dome. The Senyama volcanic products are the remains of a cone that retains a basal part 1.5 km across and 150–250 m high from the substrate. The cone comprises dacite block-and-ash flow deposits and minor base-surge deposits occur at the base. Single beds of the block-and-ash flow deposits are 1–16 m thick and dip inward 20–40° at the base of the cone and inward or outward 10–20° at the summit. Juvenile fragments in the block-and-ash flow deposits are non- to poorly vesicular and commonly have curviplanar surfaces and prismatic joints extending inward from the surfaces, which imply quenching and brittle fracturing of dacite lava. They are variably hydrothermally altered. Nevertheless, juvenile blocks appear to retain a uniform direction of the magnetization vector residual during thermal demagnetization between 280°C and 625°C. At the time of the eruption, the well-sorted sand of the substrate was at the coast and a good aquifer that facilitated explosive interaction of water and the ascending dacite lava. The mechanism of the explosion perhaps involved thermal contraction cracking of the dacite lava, water-inflow into the interior of the lava, and explosive expansion of the water. Initial phreatomagmatic explosions opened the vent. Succeeding phreatomagmatic or phreatomagmatic–vulcanian explosions produced block-and-ash flow deposits around the vent. Hydrothermal silver-ore deposits and manganese-oxide deposits occur in the Senyama volcanic products and the underlying sandstone, respectively. They could represent post-eruptive activity of the hydrothermal system developed in and around the cone.     

Eruptive style of the young high-Mg basaltic-andesite Pelagatos scoria cone,southeast of México City

The eruption of the Pelagatos scoria cone in the Sierra Chichinautzin monogenetic field near the southern suburbs of Mexico City occurred less than 14,000 years ago. The eruption initiated at a fissure with an effusive phase that formed a 7-km-long lava flow, and continued with a phase of alternating and/or simultaneous explosive and effusive activity that built a 50-m-high scoria cone on the western end of the fissure and formed a compound lava flow-field near the vent. The eruption ended with the emplacement of a short lava flow that breached the cone and was accompanied by weak explosions at the crater. Products consist of a microlite-rich high-Mg basaltic andesite. Samples were analyzed to determine the magma’s initial properties as well as the effects of degassing-induced crystallization on eruptive style. Although distal ash fallout deposits from this eruption are not preserved, a recent quarry exposes a large section of the scoria cone. Detailed study of exposed layers allows us to elucidate the mode of cone-building activity. Petrological and textural data, combined with models calibrated by experimental work and melt-inclusion analyses of similar magmas elsewhere, indicate that the magma was initially hot (>1,200°C), gas-rich (up to 5 wt.% H₂O), crystal-poor (~10 vol.% Fo₉₀ olivine phenocrysts) and thus poorly viscous (40–80 Pa s). During the early phase, low magma ascent velocity at the fissure vent allowed low-viscosity magma to degas and crystallize during ascent, producing lava flows with elevated crystal contents at T < 1,100°C, and blocky surfaces. Later, the closure of the fissure by cooling dikes focused the magma flow at a narrow section of the fissure. This led to an increased magma ascent velocity. Rapid and shallow degassing (<3 km deep) triggered ~40 vol.% microlite crystallization. Limited times for gas-escape and higher magma viscosity (6 × 10⁵–4 × 10⁶ Pa s) drove strong explosions of highly (60–80 vol.%) and finely vesicular magma. Coarse clasts broke on landing, which implies brittle behavior due to complete solidification. This requires sufficient time to cool and in turn implies ejection heights of over 1 km, which is much higher than “normal” Strombolian activity. Hence, magma viscosity significantly impacts eruption style at monogenetic volcanoes because it affects the kinetics of shallow degassing. The long-lasting eruptions of Jorullo and Paricutin, which produced similar magmas in western México, were more explosive. This can be related to higher magma fluxes and total erupted volumes. Implications of this study are important because basaltic andesites are commonly erupted to form monogenetic scoria cones of the Trans-Mexican Volcanic Belt.     

A laboratory model of surface crust formation and disruption on lava flows through non-uniform channels

Crust formation on basaltic lava flows dictates conditions of both flow cooling and emplacement. For this reason, flow histories are dramatically different depending on whether lava is transported through enclosed lava tubes or through open channels. Recent analog experiments in straight uniform channels (Griffiths et al. J Fluid Mech 496:33–62, 2003) have demonstrated that tube flow, dictated by a stationary surface crust, can be distinguished from a mobile crust regime, where a central solid crust is separated from channel walls by crust-free shear zones, by a simple dimensionless parameter ϑ, such that ϑ<25 produces tube flow and ϑ>25 describes the mobile crust regime. ϑ combines a previously determined parameter ψ, which describes the balance between the formation rate of surface solid and the shear strain that disrupts the solid crust, with the effects of thermal convection (described by the Rayleigh number Ra).Here we explore ways in which ϑ can be used to describe the behavior of basaltic lava channels. To do this we have extended the experimental approach to examine the effects of channel irregularities (expansions, contractions, sinuosity, and bottom roughness) on crust formation and disruption. We find that such changes affect local flow behavior and can thus change channel values of ϑ. For example, gradual widening of a channel results in a decrease in flow velocity that causes a decrease in ϑ and may allow a down-flow transition from the mobile crust to the tube regime. In contrast, narrowing of the channel causes an increase in flow velocity (increasing ϑ), thus inhibiting tube formation.We also quantify the fraction of surface covered by crust in the mobile crust regime. In shallow channels, variations in crust width (d _c) with channel width (W) are predicted to follow d _c∼W ^5/3. Analysis of channelized lava flows in Hawaii shows crustal coverage consistent with this theoretical result along gradually widening or narrowing channel reaches. An additional control on crustal coverage in both laboratory and basaltic flows is disruption of surface crust because of flow acceleration through constrictions, around bends, and over breaks in slope. Crustal breakage increases local rates of cooling and may cause local blockage of the channel, if crusts rotate and jam in narrow channel reaches. Together these observations illustrate the importance of both flow conditions and channel geometry on surface crust development and thus, by extension, on rates and mechanisms of flow cooling. Moreover, we note that this type of analysis could be easily extended through combined use of FLIR and LiDAR imaging to measure crustal coverage and channel geometry directly.Editorial responsibility: A. Harris     

The morphology and evolution of the Stromboli 2002–2003 lava flow field: an example of a basaltic flow field emplaced on a steep slope

The use of a hand-held thermal camera during the 2002–2003 Stromboli effusive eruption proved essential in tracking the development of flow field structures and in measuring related eruption parameters, such as the number of active vents and flow lengths. The steep underlying slope on which the flow field was emplaced resulted in a characteristic flow field morphology. This comprised a proximal shield, where flow stacking and inflation caused piling up of lava on the relatively flat ground of the vent zone, that fed a medial–distal lava flow field. This zone was characterized by the formation of lava tubes and tumuli forming a complex network of tumuli and flows linked by tubes. Most of the flow field was emplaced on extremely steep slopes and this had two effects. It caused flows to slide, as well as flow, and flow fronts to fail frequently, persistent flow front crumbling resulted in the production of an extensive debris field. Channel-fed flows were also characterized by development of excavated debris levees in this zone (Calvari et al. 2005). Collapse of lava flow fronts and inflation of the upper proximal lava shield made volume calculation very difficult. Comparison of the final field volume with that expecta by integrating the lava effusion rates through time suggests a loss of ~70% erupted lava by flow front crumbling and accumulation as debris flows below sea level. Derived relationships between effusion rate, flow length, and number of active vents showed systematic and correlated variations with time where spreading of volume between numerous flows caused an otherwise good correlation between effusion rate, flow length to break down. Observations collected during this eruption are useful in helping to understand lava flow processes on steep slopes, as well as in interpreting old lava–debris sequences found in other steep-sided volcanoes subject to effusive activity.     

A quantitative look at the demise of a basaltic vent: the death of Kupaianaha,Kilauea Volcano,Hawai'i

 The Kupaianaha vent, the source of the 48th episode of the 1983-to-present Pu'u 'O'o–Kupaianaha eruption, erupted nearly continuously from July 1986 until February 1992. This investigation documents the geophysical and geologic monitoring of the final 10 months of activity at the Kupaianaha vent. Detailed very low frequency (VLF) electromagnetic profiles across the single lava tube transporting lava from the vent were used to determine the cross-sectional area of the molten lava within the tube. Combined with measurements of lava velocity, these data provide an estimate of the lava output of Kupaianaha. In addition, lava temperatures (calculated from analysis of quenched glass) and bulk-rock chemistry were obtained for samples taken from the tube at the same site. The combined data set shows the lava flux from Kupaianaha vent declining linearly from 250 000 m³/day in April 1991 to 54 000 m³/day by November 1991. During that time surface breakouts of lava from weak points along the tube occurred progressively closer to the vent, consistent with declining efficiency in lava transport. There were no significant changes in lava temperature or in bulk MgO content during this period. Another eruptive episode (the 49th) began uprift of Kupaianaha on 8 November 1991 and erupted lava concurrently with Kupaianaha for 18 days. Lava flux from Kupaianaha decreased in response to this new episode, but the response was delayed by approximately 1 day. After 14 November 1991, lava velocities were no longer measurable in the tube because the lava stream beneath the skylight had crusted over; however, the VLF-derived electrical conductances documented the decreasing flux of molten lava through the tube. Kupaianaha remained active, but output continued to decrease until early February 1992 when the last active surface flows were seen. In November 1991 we used the linearly decreasing effusion rate to accurately predict the date for the death of the Kupaianaha vent. The linear nature of the decline in lava tube conductance and the delayed and slow response of the Waha'ula tube conductances to the 49th eruptive episode led us to speculate that (a) the Kupaianaha vent shut down because of a decrease in driving pressure and not because of a freeze-up of the vent, and (b) that Pu'u 'O'o, episode 49, and Kupaianaha were fed nearly vertically from a source deep within the rift zone. Received: 29 September 1995 / Accepted: 21 November 1995     

Rootless shield and perched lava pond collapses at Kīlauea Volcano,Hawai'i

Effusion rate is a primary measurement used to judge the expected advance rate, length, and hazard potential of lava flows. At basaltic volcanoes, the rapid draining of lava stored in rootless shields and perched ponds can produce lava flows with much higher local effusion rates and advance velocities than would be expected based on the effusion rate at the vent. For several months in 2007–2008, lava stored in a series of perched ponds and rootless shields on Kīlauea Volcano, Hawai'i, was released episodically to produce fast-moving 'a'ā lava flows. Several of these lava flows approached Royal Gardens subdivision and threatened the safety of remaining residents. Using time-lapse image measurements, we show that the initial time-averaged discharge rate for one collapse-triggered lava flow was approximately eight times greater than the effusion rate at the vent. Though short-lived, the collapse-triggered 'a'ā lava flows had average advance rates approximately 45 times greater than that of the pāhoehoe flow field from which they were sourced. The high advance rates of the collapse-triggered lava flows demonstrates that recognition of lava accumulating in ponds and shields, which may be stored in a cryptic manner, is vital for accurately assessing short-term hazards at basaltic volcanoes.     

Viscosity of microlite-bearing rhyolitic obsidians: an experimental study

 To investigate the influence of microlites on lava flow rheology, the viscosity of natural microlite-bearing rhyolitic obsidians of calc-alkaline and peralkaline compositions containing 0.1–0.4 wt.% water was measured at volcanologically relevant temperatures (650–950  °C), stresses (10³–10⁵ Pa) and strain rates (10^–5 to 10^–7 s^–1). The glass transition temperatures (T _g ) were determined from scanning calorimetric measurements on the melts for a range of cooling/heating rates. Based on the equivalence of enthalpic (calorimetric) and shear (viscosity) relaxation, we calculated the viscosity of the melt in crystal-bearing samples from the T _g data. The difference between the calculated viscosity of the melt phase and the measured viscosity for the crystal-bearing samples is interpreted to be the physical effect of microlites on the measured viscosity. The effect of <5 vol.% rod-like microlites on the melt rheology is negligible. Microlite-rich and microlite-poor samples from the same lava flow and with identical bulk chemistry show a difference of 0.6 log₁₀ units viscosity (Pa s), interpreted to be due to differences in melt chemistry caused by the presence of microlites. The only major differences between measured and calculated viscosities were for two samples: a calc-alkaline rhyolite with 1 vol.% branching crystals, and a peralkaline rhyolite containing crystal-rich bands with >45 vol.% crystals. For both of these samples a connectivity factor is apparent, with, for the latter, a close packing framework of crystals which is interpreted to influence the apparent viscosity. Received: 14 March 1996 / Accepted: 30 May 1996     

Thermal structure and heat loss at the summit crater of an active lava dome

Forward-Looking Infrared (FLIR) nighttime thermal images were used to extract the thermal and morphological properties for the surface of a blocky-to-rubbley lava mass active within the summit crater of the Caliente vent at Santiaguito lava dome (Guatemala). Thermally the crater was characterized by three concentric regions: a hot outer annulus of loose fine material at 150–400°C, an inner cold annulus of blocky lava at 40–80°C, and a warm central core at 100–200°C comprising younger, hotter lava. Intermittent explosions resulted in thermal renewal of some surfaces, mostly across the outer annulus where loose, fine, fill material was ejected to expose hotter, underlying, material. Surface heat flux densities (radiative + free convection) were dominated by losses from the outer annulus (0.3–1.5 × 10⁴  s⁻¹m⁻²), followed by the hot central core (0.1–0.4 × 10⁴ J s⁻¹m⁻²) and cold annulus (0.04–0.1 × 10⁴ J s⁻¹m⁻²). Overall surface power output was also dominated by the outer annulus region (31–176 MJ s⁻¹), but the cold annulus contributed equal power (2.41–7.07 MJ s⁻¹) as the hot central core (2.68–6.92 MJ s⁻¹) due to its greater area. Cooled surfaces (i.e. the upper thermal boundary layer separating surface temperatures from underlying material at magmatic temperatures) across the central core and cold annulus had estimated thicknesses, based on simple conductive model, of 0.3–2.2 and 1.5–4.3 m. The stability of the thermal structure through time and between explosions indicates that it is linked to a deeper structural control likely comprising a central massive plug, feeding lava flow from the SW rim of the crater, surrounded by an arcuate, marginal fracture zone through which heat and mass can preferentially flow.     

Explosive lava–water interactions I: architecture and emplacement chronology of volcanic rootless cone groups in the 1783–1784 Laki lava flow, Iceland

To determine the relationships between rootless cone emplacement mechanisms, morphology, and spatial distribution, we mapped the Hnúta and Hrossatungur groups of the 1783–1784 Laki lava flow in Iceland. We based our facies maps on Differential Global Positioning System (DGPS) measurements, photogeological interpretations, and supporting field observations. The study area covers 2.77 km² and includes 2216 explosion sites. To establish the timing of rootless cone formation we incorporated tephrochronological constraints from eighty-eight stratigraphic sections and determined that the Hnúta and Hrossatungur groups are composite structures formed by the emplacement of six geographically and chronologically discrete domains. Rootless eruptions initiated in domain 1 on the first day of the Laki eruption (June 8, 1783) and lasted 1–2 days. The second episode of rootless activity began in domain 2 on June 11 and lasted 1–3 days. The four domains of the Hrossatungur group dominantly formed after June 14 and exhibit a complex emplacement sequence that reflects interactions between the Laki lava, contemporaneously emplaced rootless cones, and an existing topographic ridge. In the study area, we identify three distinct rootless cone archetypes (i.e., recurring morphological forms) that are related to tube-, channel-, and broad sheet lobe-fed eruptions. We assert that emplacement of lava above compressible substrates (e.g., unconsolidated sediments) may trigger rootless eruptions by causing subsidence-induced flexure and failure of the basal crust, thereby allowing molten lava (fuel) to come into direct contact with groundwater (coolant) and initiating analogs to explosive molten fuel–coolant interactions (MFCIs).     

Thick lava flows of Karisimbi Volcano, Rwanda: insights from SIR-C interferometric topography

 We use a digital elevation model (DEM) derived from interferometrically processed SIR-C radar data to estimate the thickness of massive trachyte lava flows on the east flank of Karisimbi Volcano, Rwanda. The flows are as long as 12 km and average 40–60 m (up to >140 m) in thickness. By calculating and subtracting a reference surface from the DEM, we derived a map of flow thickness, which we used to calculate the volume (up to 1 km³ for an individual flow, and 1.8 km³ for all the identified flows) and yield strength of several flows (23–124 kPa). Using the DEM we estimated apparent viscosity based on the spacing of large folds (1.2×10¹² to 5.5×10¹² Pa s for surface viscosity, and 7.5×10¹⁰ to 5.2×10¹¹ Pa s for interior viscosity, for a strain interval of 24 h). We use shaded-relief images of the DEM to map basic flow structures such as channels, shear zones, and surface folds, as well as flow boundaries. The flow thickness map also proves invaluable in mapping flows where flow boundaries are indistinct and poorly expressed in the radar backscatter and shaded-relief images. Received: 6 September 1997 / Accepted: 15 May 1998     

Hydrothermal hexahydrite spherules erupted during the 2008–2010 summit eruption of Kīlauea Volcano,Hawai`i

Small (1–3 mm), hollow spherules of hexahydrite have been collected falling out of the magmatic gas plume downwind of Kīlauea’s summit vent. The spherules were observed on eight separate occasions during 2009–2010 when a lake of actively spattering lava was present ~150–200 m below the rim of the vent. The shells of the spherules have a fine bubbly foam structure less than 0.1 mm thick, composed almost entirely of hexahydrite [MgSO₄·6H₂O] Small microspherules of lava (<5 μm across) along with mineral and rock fragments from the magmatic plume adhered to the outside of the hexahydrite spherules. Phase relationships and the particulate matter in the magmatic plume indicate that the spherules originated as a bubbly solution injected into and mixed with the magmatic plume. The most likely mechanism for production of hexahydrite spherules is boiling of MgSO₄-saturated meteoric water in the walls of the conduit above the surface of the lava lake. Solfataric sulfates may thus be recycled and reinjected into the plume, creating particulates of sulfate minerals that can be distributed far from their original source.     

Voluminous lava flows at Oldoinyo Lengai in 2006: chronology of events and insights into the shallow magmatic system

The largest natrocarbonatite lava flow eruption ever documented at Oldoinyo Lengai, NW Tanzania, occurred from March 25 to April 5, 2006, in two main phases. It was associated with hornito collapse, rapid extrusion of lava covering a third of the crater and emplacement of a 3-km long compound rubbly pahoehoe to blocky aa-like flow on the W flank. The eruption was followed by rapid enlargement of a pit crater. The erupted natrocarbonatite lava has high silica content (3% SiO₂). The eruption chronology is reconstructed from eyewitness and news media reports and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data, which provide the most reliable evidence to constrain the eruption’s onset and variations in activity. The eruption products were mapped in the field and the total erupted lava volume estimated at 9.2 ± 3.0 × 10⁵ m³. The event chronology and field evidence are consistent with vent construct instability causing magma mixing and rapid extrusion from shallow reservoirs. It provides new insights into and highlights the evolution of the shallow magmatic system at this unique natrocarbonatite volcano.     

Effusive eruption of viscous silicic magma triggered and driven by recharge: a case study of the Cerro Chascon-Runtu Jarita Dome Complex in Southwest Bolivia

 The Cerro Chascon-Runtu Jarita Complex is a group of ten Late Pleistocene (∼85 ka) lava domes located in the Andean Central Volcanic Zone of Bolivia. These domes display considerable macroscopic and microscopic evidence of magma mixing. Two groups of domes are defined chemically and geographically. A northern group, the Chascon, consists of four lava bodies of dominantly rhyodacite composition. These bodies contain 43–48% phenocrysts of plagioclase, quartz, sanidine, biotite, and amphibole in a microlite-poor, rhyolitic glass. Rare mafic enclaves and selvages are present. Mineral equilibria yield temperatures from 640 to 750  °C and log ƒO₂ of –16. Geochemical data indicate that the pre-eruption magma chamber was zoned from a dominant volume of 68% to minor amounts of 76% SiO₂. This zonation is best explained by fractional crystallization and some mixing between rhyodacite and more evolved compositions. The mafic enclaves represent magma that intruded but did not chemically interact much with the evolved magmas. A southern group, the Runtu Jarita, is a linear chain of six small domes (<1 km³ total volume) that probably is the surface expression of a dike. The five most northerly domes are composites of dacitic and rhyolitic compositions. The southernmost dome is dominantly rhyolite with rare mafic enclaves. The composite domes have lower flanks of porphyritic dacite with ∼35 vol.% phenocrysts of plagioclase, orthopyroxene, and hornblende in a microlite-rich, rhyodacitic glass. Sieve-textured plagioclase, mixed populations of disequilibrium plagioclase compositions, xenocrystic quartz, and sanidine with ternary composition reaction rims indicate that the dacite is a hybrid. The central cores of the composite domes are rhyolitic and contain up to 48 vol.% phenocrysts of plagioclase, quartz, sanidine, biotite, and amphibole. This is separated from the dacitic flanks by a banded zone of mingled lava. Macroscopic, microscopic, and petrologic evidence suggest scavenging of phenocrysts from the silicic lava. Mineral equilibria yield temperatures of 625–727  °C and log ƒO₂ of –16 for the rhyolite and 926–1000  °C and log ƒO₂ of –9.5 for the dacite. The rhyolite is zoned from 73 to 76% SiO₂, and fractionation within the rhyolite composition produced this variation. Most of the 63–73% SiO₂ compositional range of the lava in this group is the result of mixing between the hybrid dacite and the rhyolite. Eruption of both groups of lavas apparently was triggered by mafic recharge. A paucity of explosive activity suggests that volatile and thermal exchanges between reservoir and recharge magmas were less important than volume increase and the lubricating effects of recharge by mafic magmas. For the Runtu Jarita group, the eruption is best explained by intrusion of a dike of dacite into a chamber of crystal-rich rhyolite close to its solidus. The rhyolite was encapsulated and transported to the surface by the less-viscous dacite magma, which also acted as a lubricant. Simultaneous effusion of the lavas produced the composite domes, and their zonation reflects the subsurface zonation. The role of recharge by hotter, more fluid mafic magma appears to be critical to the eruption of some highly viscous silicic magmas. Received: 23 August 1998 / Accepted: 10 March 1999     

Lava channel roofing,overflows, breaches and switching: insights from the 2008–2009 eruption of Mt. Etna

During long-lived basaltic eruptions, overflows from lava channels and breaching of channel levées are important processes in the development of extensive 'a'ā lava flow-fields. Short-lived breaches result in inundation of areas adjacent to the main channel. However, if a breach remains open, lava supply to the original flow front is significantly reduced, and flow-field widening is favoured over lengthening. The development of channel breaches and overflows can therefore exert strong control over the overall flow-field development, but the processes that determine their location and frequency are currently poorly understood. During the final month of the 2008–2009 eruption of Mt. Etna, Sicily, a remote time-lapse camera was deployed to monitor events in a proximal region of a small ephemeral lava flow. For over a period of ~10 h, the flow underwent changes in surface elevation and velocity, repeated overflows of varying vigour and the construction of a channel roof (a required prelude to lava tube formation). Quantitative interpretation of the image sequence was facilitated by a 3D model of the scene constructed using structure-from-motion computer vision techniques. As surface activity waned during the roofing process, overflow sites retreated up the flow towards the vent, and eventually, a new flow was initiated. Our observations and measurements indicate that flow surface stagnation and flow inflation propagated up-flow at an effective rate of ~6 m h⁻¹, and that these processes, rather than effusion rate variations, were ultimately responsible for the most vigorous overflow events. We discuss evidence for similar controls during levée breaching and channel switching events on much larger flows on Etna, such as during the 2001 eruption.