Isothermal and undercooling experiments were conducted on one of the most primitive trachybasalts from Mt. Etna volcano in order to examine the crystallization mechanisms controlling the textural and ...compositional variability of clinopyroxene. Experiments were performed at 400–800 MPa, 1050–1200 °C, 0–4 wt.% H2O and at oxygen fugacity 2 log units above the Ni-NiO + 2 buffer. In isothermal experiments, the final resting temperature is approached from room temperature and clinopyroxene growth is dominated by an interface-controlled mechanism, leading to the formation of small (∼10 µm) and euhedral crystals with homogeneous compositions. Conversely, in undercooling experiments, the final resting temperature is approached after annealing at temperature above the liquidus, imposing an effective degree undercooling (ΔT) to the system. In presence of undercooling, the crystallization of clinopyroxene is dominated by a diffusion-controlled mechanism that determines the formation of large (>100 µm) crystals, constituted by two compositionally distinct domains, enriched in Al2O3 + TiO2 and SiO2 + MgO, respectively. The maximum growth rate (Gmax) decreases progressively from ∼10−7 to ∼10−8 cm/s as the degree of undercooling increases from ∼20 to ∼230 °C, due to the increase in nucleation rate. At low to moderate degrees of undercooling (ΔT = 23–41 °C) clinopyroxene is prevalently euhedral to subhedral, whereas at high degrees of undercooling, the crystal shape changes from prevalently subhedral (ΔT = 73–123 °C) to skeletal and dendritic (ΔT = 132–233 °C). Hourglass sector zoning similar to that documented for natural phenocrysts from eruptions at Mt. Etna volcano is observed only at low degrees of undercooling (ΔT = 23–32 °C). This type of zoning develops in the form of the cation exchange Si + Mg{-111} ↔ Al + Ti{100} and demonstrates that hourglass sector zoning is an effective indicator of sluggish kinetic effects caused by relatively low degrees of undercooling. In contrast, at increasing degrees of undercooling (ΔT > 32 °C), strong melt supersaturation determines the early formation of Al2O3 + TiO2-rich dendritic crystals and further SiO2 + MgO-rich overgrowths, as the bulk system attempts to return to a near-equilibrium state between the advancing crystal surface and the feeding melt.
The experimentally-determined relationship between ΔT and clinopyroxene chemistry is used to reconstruct the crystallization conditions of natural clinopyroxenes from 1974 and 2002–2003 eccentric eruptions at Mt. Etna volcano. Clinopyroxene rims record much higher degrees of undercooling (up to ∼110 °C) than crystal mantles associated with magma recharge at depth (mostly 0–40 °C). Hence, the rims track decompression-induced degassing and cooling during the ascent of magma towards the surface.
We present the variation in trace element partition coefficients measured at the interface between rapidly cooled clinopyroxene crystals and co-existing melts. Results indicate that, as the cooling ...rate is increased, clinopyroxene crystals are progressively depleted in Si, Ca and Mg counterbalanced by enrichments in Al (mainly tetrahedral Al
iv
), Na and Ti. Partition coefficients (Ds) for rare earth elements (REE), high field strength elements (HFSE) and transition elements (TE) increase with increasing cooling rate, in response to clinopyroxene compositional variations. The entry of REE into the M2 site is facilitated by a coupled substitution where either Na substitutes for Ca on the M2 site or Al
iv
substitutes for Si in the tetrahedral site. The latter substitution reflects an increased ease of locally balancing the excess charge at M2 as the number of surrounding Al
iv
atoms increases. Due to the lower concentration of Ca in rapidly cooled clinopyroxenes, divalent large ion lithophile elements (LILE) on M2 decrease at the expense of monovalent cations. Conversely, higher concentrations of HFSE and TE on the M1 site are facilitated as the average charge on this site increases with the replacement of divalent-charged cations by Al
vi
. Although crystallization kinetics modify clinopyroxene composition, deviations from equilibrium partitioning are insufficient to change the tendency of a trace element to be compatible or incompatible. Consequently, there are regular relationships between ionic radius, valence of the trace element and D. At both equilibrium and cooling rate conditions, Ds for isovalent cations define parabola-like curves when plotted against ionic radius, consistent with the lattice strain model, demonstrating that the partitioning of trace elements is driven by charge balance mechanisms; cation substitution reactions can be treated in terms of the energetics of the various charge-imbalanced configurations.
Strombolian activity is common in low-viscosity volcanism. It is characterised by quasi-periodic, short-lived explosions, which, whilst typically weak, may vary greatly in magnitude. The current ...paradigm for a strombolian volcanic eruption postulates a large gas bubble (slug) bursting explosively after ascending a conduit filled with low-viscosity magma. However, recent studies of pyroclast textures suggest the formation of a region of cooler, degassed, more-viscous magma at the top of the conduit is a common feature of strombolian eruptions. Following the hypothesis that such a rheological impedance could act as a ‘viscous plug’, which modifies and complicates gas escape processes, we conduct the first experimental investigation of this scenario. We find that: 1) the presence of a viscous plug enhances slug burst vigour; 2) experiments that include a viscous plug reproduce, and offer an explanation for, key phenomena observed in natural strombolian eruptions; 3) the presence and extent of the plug must be considered for the interpretation of infrasonic measurements of strombolian eruptions. Our scaled analogue experiments show that, as the gas slug expands on ascent, it forces the underlying low-viscosity liquid into the plug, creating a low-viscosity channel within a high-viscosity annulus. The slug's diameter and ascent rate change as it enters the channel, generating instabilities and increasing slug overpressure. When the slug reaches the surface, a more energetic burst process is observed than would be the case for a slug rising through the low-viscosity liquid alone. Fluid-dynamic instabilities cause low and high viscosity magma analogues to intermingle, and cause the burst to become pulsatory. The observed phenomena are reproduced by numerical fluid dynamic simulations at the volcanic scale, and provide a plausible explanation for pulsations, and the ejection of mingled pyroclasts, observed at Stromboli and elsewhere.
•We present laboratory experiments in which gas slugs ascend a pipe plugged with a viscous liquid.•The presence of a viscous plug enhances burst vigour.•Slug–plug interaction causes pulsatory bursting resembling natural strombolian eruptions.•Experiments indicate slug passage promotes mingling between vertically stratified magmas.
The explosive activity of the 2021 Tajogaite eruption eludes pigeonholing into well‐defined eruption styles, with a variety of pyroclast ejection modes occurring both alternately and simultaneously ...at multiple vents. Visually, we defined four endmembers of explosive activity, referred to as fountaining, spattering, ash‐poor jets and ash‐rich jets. To capture the physical parameters of these activities, we deployed a camera array including one high‐speed camera and three high‐definition cameras in two field campaigns. Transitions between and fluctuations within activity occurred at the time scale of minutes to hours, likely driven by the same shallow conduit and vent processes controlling Strombolian activity at other volcanoes, but at higher gas and magma fluxes. From a physical standpoint, mean pyroclast rise velocity ranged 5–50 m/s, maximum ejection velocity 10–220 m/s, and sub‐second mass flux of lapilli to bomb‐sized pyroclasts at the vent 0.2–200 × 103 kg/s. The largest mass flux occurred during fountaining, which contributed by far more than other activities to cone building. All explosive activity exhibited well‐defined pyroclast ejection pulses, and we found a positive correlation between the occurrence rate of ejection pulses and maximum pyroclast ejection velocity. Despite orders of magnitude variations, physical parameters shift gradually with no boundary from one activity endmember to another. As such, attributing this explosive activity specifically to any currently defined style variations is arbitrary and potentially misleading. The highly variable explosive activity of the Tajogaite eruption recalls previous definitions of violent Strombolian eruptions, an eruption style whose pyroclast ejection dynamics, however, were so far largely undefined.
Plain Language Summary
The 2021 Tajogaite volcanic eruption offered a rare opportunity to study in detail the physical properties and the controlling factors of explosive activity driven by basaltic magmas. The activity lasted almost uninterrupted for almost 3 months and had visually different manifestations occurring simultaneously and alternating at different volcanic vents. To study the explosive activity, we used one high‐speed camera, taking short, slow motion videos, and three commercial grade high‐definition camcorders recording for many hours. We found that the activity changed in features and intensity at the time scale of minutes to hours, largely controlled by changes in the size and debris cover of the vent, magma viscosity, and magma flux and gas content. The ejection velocity of large volcanic particles ranged 5–220 m/s, with mean values around 10–50 m/s. The mass flux of particles erupted reached peaks of 200 metric tons per second. Particle ejection was never steady but always proceeded in pulses, which were more frequent if the ejection velocity was higher. Our measurements show that the current classification schemes for explosive eruptions of basaltic magmas do not adequately describe the activity of the Tajogaite eruption, which represents a type of eruption that was not yet measured in detail.
Key Points
High‐definition and high‐speed imaging record the velocity, size, and mass flux of pyroclasts
Activity shifted in location, nature and vigor at the time scale of hours and progressed in ejection pulses at the time scale of seconds
Physical parameters of explosive activity vary gradually between apparently different activity styles, without any clear boundary
Open‐conduit conditions characterize several of the most hazardous and active volcanic systems of basaltic composition worldwide, persistently refilled by magmatic inputs. Eruptive products with ...similar bulk compositions, chemically buffered by continual mafic inputs, nevertheless exhibit heterogeneous glass compositions in response to variable magma mixing, crystallization, and differentiation processes within different parts of the plumbing system. Here, we document how multivariate statistics and magma differentiation modeling based on a large data set of glass compositions can be combined to constrain magma differentiation and plumbing system dynamics. Major and trace elements of matrix glasses erupted at Stromboli volcano (Italy) over the last 20 years provide a benchmark against which to test our integrated petrological approach. Principal component analysis, K‐means cluster analysis, and kernel density estimation reveal that trace elements define a multivariate space whose eigenvectors are more readily interpretable in terms of petrological processes than major elements, leading to improved clustering solutions. Comparison between open‐ and closed‐system differentiation models outlines that steady state magma compositions at constantly replenished and erupting magmatic systems approximate simple fractional crystallization trends, due to short magma residence times. Open‐system magma evolution is associated with magma storage crystallinities that are lower than those associated with closed‐system scenarios. Accordingly, open‐system dynamics determine the efficient crystal‐melt separation toward the top of the reservoir, where eruptible melts continuously supply the ordinary activity. Conversely, a mush‐like environment constitutes the bottom of the reservoir, where poorly evolved magmas result from mixing events between mush residual melts and primitive magmas injected from deeper crustal levels.
Plain Language Summary
Volcanoes characterized by continuous eruptive activity are typified by constant replenishment of new magma, rising from deeper regions of the crust. The volcanic glass (supercooled silicate melt), represents the residual liquid of magma crystallization, and is found as the intracrystalline matrix of eruptive products. The study of its chemical composition may provide insight into the processes occurring at depths beneath the volcanic vent, where magma compositional changes result from crystallization and mixing with new magma rising from depth. We combine statistical analyses and analytical equations based on the chemical composition of the matrix glasses from Stromboli volcano, in order to constrain the processes which produce their chemical variations, identifying different environments where magmas are stored at depth. Our results also show that when magma is stored for a short period of time, the chemical changes to which the magma is subjected in a constantly replenished system are similar to those occurring in a system which is closed to new inputs of magma.
Key Points
The combination of multivariate statistics with geochemical modeling provides new constraints on magma differentiation processes
Multivariate statistics based on trace elements allow better retrieval of petrological information than those based on major elements
Magma differentiation in open systems approximates that occurring in closed systems when magma residence timescales are short
For a better understanding of the final solidification history of eruptions at Mt. Etna volcano (Sicily, Italy), we have investigated cation redistributions at the interface between ...sub-millimetre-sized clinopyroxene and titanomagnetite crystal rims and coexisting melts. The studied products were scoria clasts from lava fountains and rock samples from pahoehoe and aa lava flows. Our data indicate that scoria clasts from lava fountaining were rapidly quenched at the contact with the atmosphere, preserving the original crystal textures and compositions inherited during magma dynamics within the plumbing system. Kinetics and energetics of crystallization were instantaneously frozen-in and post-eruptive effects on mineral chemistry were negligible. The near-equilibrium compositions of clinopyroxene and titanomagnetite indicate that lava fountain episodes were supplied by high-temperature, H2O-rich magmas ascending with velocities of 0.01–0.31m/s. In contrast, magmas feeding lava flow eruptions underwent a more complex solidification history where the final stage of the crystal growth was mostly influenced by volatile loss and heat dissipation at syn- and post-eruptive conditions. Due to kinetic effects associated with magma undercooling, clinopyroxenes and titanomagnetites formed by crystal attachment and agglomeration mechanisms leading to intricate intergrowth textures. The final compositions of these minerals testify to closure temperatures and melt–water concentrations remarkably lower than those estimated for lava fountains. Kinetically-controlled cation redistributions at the crystal–melt interface suggest that the solidification of magma was driven by degassing and cooling processes proceeding from the uppermost part of the volcanic conduit to the surface.
•Cation redistributions have been measured in clinopyroxenes and titanomagnetites at the final stage of crystal growth.•These compositional variations reflect the cooling and thermal path of magmas ascending from the conduit to the surface.•Remarkable implications are found for the final solidification history of lava fountains and lava flows at Mt. Etna.
Centimeter to meter‐sized volcanic ballistic projectiles from explosive eruptions jeopardize people and properties kilometers from the volcano, but they also provide information about the past ...eruptions. Traditionally, projectile trajectory is modeled using simplified ballistic theory, accounting for gravity and drag forces only and assuming simply shaped projectiles free moving through air. Recently, collisions between projectiles and interactions with plumes are starting to be considered. Besides theory, experimental studies and field mapping have so far dominated volcanic projectile research, with only limited observations. High‐speed, high‐definition imaging now offers a new spatial and temporal scale of observation that we use to illuminate projectile dynamics. In‐flight collisions commonly affect the size, shape, trajectory, and rotation of projectiles according to both projectile nature (ductile bomb versus brittle block) and the location and timing of collisions. These, in turn, are controlled by ejection pulses occurring at the vent. In‐flight tearing and fragmentation characterize large bombs, which often break on landing, both factors concurring to decrease the average grain size of the resulting deposits. Complex rotation and spinning are ubiquitous features of projectiles, and the related Magnus effect may deviate projectile trajectory by tens of degrees. A new relationship is derived, linking projectile velocity and size with the size of the resulting impact crater. Finally, apparent drag coefficient values, obtained for selected projectiles, mostly range from 1 to 7, higher than expected, reflecting complex projectile dynamics. These new perspectives will impact projectile hazard mitigation and the interpretation of projectile deposits from past eruptions, both on Earth and on other planets.
Key Points
Volcanic Ballistic Projectiles (VBPs) in volcanic deposits, theory, and direct observations are reviewed
High‐speed imaging and measurements of VBPs spinning, deforming, fragmenting, colliding, and impacting with the ground are provided
In‐flight fragmentation, collisions, and spinning are important for VBPs dynamics, and apparent drag coefficient can be higher than expected
Plain Language Summary
Explosive volcanic eruptions launch incandescent fragments, sometimes partially molten, to distances of up to several kilometers from the volcano. The largest fragments, from the size of an apple to that of a van, travel in air following the same laws that control the flight of artillery shells and, on landing, may cause the same harmful consequences. To protect people and properties from these volcanic projectiles, their occurrence in volcanic rocks is documented, and their motion is simulated by computer models. However, both field studies and computer models require validation, but in‐flight observation of the projectiles have been sparse, so far. We used state‐of‐the‐art high‐speed cameras, filming volcanic projectiles in slow motion to understand and measure the processes that control their flight dynamics. We found that the in‐flight deformation, rotation, and collision of the projectiles have a deep impact on their trajectory. We also measured the size of craters left by the projectiles on landing, and we derived specific parameters that are essential to model projectiles flight. We found that currently used models often do not account for all the in‐flight dynamics. Our findings will improve interpreting the motion of the projectiles and mitigating their hazard.
The intense explosive and effusive volcanic activity of the last 1000 years at La Fossa volcano (Vulcano Island, Italy) was characterized by the eruption of magmas ranging in composition from latites ...to trachytes and rhyolites, as well as K-rich trachytes. Evidence of syn-eruptive mixing among these magmas is frequently observed in the form of magmatic enclaves and bands in lava flows and pyroclastic products. The petrological and volcanological diversity of the erupted materials suggests that complex differentiation processes occurred in the shallow part of the plumbing system. With the aim to reconstruct the magmatic feeding system and to identify the differentiation processes behind such a petrologic complexity, we analysed lavas and pyroclastic products representative of the recent eruptive sequences at La Fossa and combined the petro-chemical features with thermo-barometric calculations, geochemical modelling and temperature gradient experiments. Thermo-barometric calculations indicate that the K-rich trachytic magma crystallized at lower pressure (160 ± 54 MPa) compared to the latitic (307 ± 47 MPa) and trachytic (208 ± 30 MPa) magmas. Differentiation modelling suggests that both trachytic and rhyolitic compositions can be obtained through differentiation of a common latitic magma, essentially by varying the plagioclase/sanidine ratio. Temperature gradient experiments, performed at the conditions inferred for the shallow plumbing system of La Fossa volcano (150 MPa and 1050–900 °C), indicate different paths of melt differentiation that overall produce an increase of the SiO
2
/K
2
O ratio with the increasing H
2
O in the system (from 0 to 4 wt.%). This is consistent with the origin of K-rich trachytes at lower pressure and lower H
2
O content. In turn, the formation of crystal-poor rhyolites is explained by the segregation of the interstitial melt formed in a latitic–trachytic crystal mush, favoured by the second boiling of the melt and consequent exsolution of a fluid phase.
Abstract
Systematic variations in the crystal cargo and whole-rock isotopic compositions of mantle-derived basalts in the intraplate Dunedin Volcano (New Zealand) indicate the influence of a complex ...mantle-to-crust polybaric plumbing system. Basaltic rocks define a compositional spectrum from low-alkali basalts through mid-alkali basalts to high-alkali basalts. High-alkali basalts display clinopyroxene crystals with sector (hourglass) and oscillatory zoning (Mg#61–82) as well as Fe-rich green cores (Mg#43–69), whereas low-alkali basalts are characterized by clinopyroxenes with unzoned overgrowths (Mg#69–83) on resorbed mafic cores (Mg#78–88), coexisting with reversely zoned plagioclase crystals (An43–68 to An60–84 from core to rim). Complex magma dynamics are indicated by distinctive compositional variations in clinopyroxene phenocrysts, with Cr-rich zones (Mg#74–87) indicating continuous recharge by more mafic magmas. Crystallization of olivine, clinopyroxene and titanomagnetite occurred within a polybaric plumbing system extending from upper mantle to mid-crustal depths (485–1059 MPa and 1147–1286°C), whereas crystallization of plagioclase with subordinate clinopyroxene and titanomagnetite proceeded towards shallower crustal levels. The compositions of high-alkali basalts and mid-alkali basalts resemble those of ocean island basalts and are characterized by FOZO-HIMU isotopic signatures (87Sr/86Sri = 0.70277–0.70315, 143Nd/144Ndi = 0.51286–0.51294 and 206Pb/204Pb = 19.348–20.265), whereas low-alkali basalts have lower incompatible element abundances and isotopic compositions trending towards EMII (87Sr/86Sri = 0.70327–70397, 143Nd/144Ndi = 0.51282–0.51286 and 206Pb/204Pb = 19.278–19.793). High- and mid-alkali basalt magmas mostly crystallized in the lower crust, whereas low-alkali basalt magma recorded deeper upper mantle clinopyroxene crystallization before eruption. The variable alkaline character and isotope composition may result from interaction of low-alkaline melts derived from the asthenosphere with melts derived from lithospheric mantle, possibly initiated by asthenospheric melt percolation. The transition to more alkaline compositions was induced by variable degrees of melting of metasomatic lithologies in the lithospheric mantle, leading to eruption of predominantly small-volume, high-alkali magmas at the periphery of the volcano. Moreover, the lithosphere imposed a filtering effect on the alkalinity of these intraplate magmas. As a consequence, the eruption of low-alkali basalts with greater asthenospheric input was concentrated at the centre of the volcano, where the plumbing system was more developed.