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.
Recurring discoveries of abiotic methane in gas seeps and springs in ophiolites and peridotite massifs worldwide raised the question of where, in which rocks, methane was generated. Answers will ...impact the theories on life origin related to serpentinization of ultramafic rocks, and the origin of methane on rocky planets. Here we document, through molecular and isotopic analyses of gas liberated by rock crushing, that among the several mafic and ultramafic rocks composing classic ophiolites in Greece, i.e., serpentinite, peridotite, chromitite, gabbro, rodingite and basalt, only chromitites, characterized by high concentrations of chromium and ruthenium, host considerable amounts of
C-enriched methane, hydrogen and heavier hydrocarbons with inverse isotopic trend, which is typical of abiotic gas origin. Raman analyses are consistent with methane being occluded in widespread microfractures and porous serpentine- or chlorite-filled veins. Chromium and ruthenium may be key metal catalysts for methane production via Sabatier reaction. Chromitites may represent source rocks of abiotic methane on Earth and, potentially, on Mars.
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
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.
The petrological study of volcanic products emitted during the paroxysmal events of December 2015 from the summit craters of Mount Etna allow us to constrain T-P-XH
2
O phase stability, ...crystallization conditions, and mixing processes along the main open-conduit feeding system. In this study, we discuss new geochemical, thermo-barometric data and related Rhyolite-MELTS modelling of the eruptive activity that involved the concomitant activation of all summit craters. The results, in comparison with the previous paroxysmal events of the 2011–2012, reinforce the model of a vertically extended feeding system and highlight that the activity at the New South-East Crater was fed by magma residing at a significantly shallower depth with respect to the Central Craters (CC) and North-East Crater (NEC), even if all conduits were fed by a common deep (
P
= 530–440 MPa) basic magmatic input. Plagioclase dissolution, resorption textures, and the Rhyolite-MELTS stability model corroborate its dependence on H
2
O content; thus, suggesting that further studies on the effect that flushing from fluids with different H
2
O/CO
2
ratio are needed to understand the eruption-triggering mechanisms for high energetic strombolian paroxysmal episodes.
The aim of this study was to investigate the fate of the conjugated forms of the three most common natural estrogens in the municipal aqueous environment. Levels of conjugated and free estrogens in ...(1) female urine; (2) a septic tank collecting domestic wastewater; (3) influents and effluents of six activated sludge sewage treatment plants (STPs) were measured. The analytical method was based on solid-phase extraction by using a Carbograph 4 cartridge and Liquid Chromatography-tandem Mass Spectrometry. On average, a group of 73 women selected to represent a typical cross section of the female inhabitants of a Roman condominium, excreted 106, 14 and 32 μg/day of conjugated estriol (E
3), estradiol (E
2) and estrone (E
1), respectively. Apart from some E
3 in pregnancy urine, free estrogens were never detected in urine samples. Estrogen sulfates represented 21% of the total conjugated estrogens. This situation changed markedly in the condominium collecting tank. Here, significant amounts of free estrogens were observed and the estrogen sulfate to estrogen glucuronated ratio rose to 55/45. A laboratory biodegradation test confirmed that glucuronated estrogens are readily deconjugated in unmodified domestic wastewater, presumably due to the large amounts of the β-glucuronidase enzyme produced by fecal bacteria (
Escherichia coli). Deconjugation continued in sewer transit. At the STP entrance, free estrogens and sulfated estrogens were the dominant species. The sewage treatment completely removed residues of estrogen glucuronates and with good efficiency (84–97%) the other analytes, but not E
1 (61%) and estrone-3-sulfate (E
1-3S) (64%). Considering that (1) E
1 has half the estrogenic potency of E
2, (2) the amount of the former species discharged from STPs into the receiving water was more than ten times larger than the latter one and (3) a certain fraction of E
1-3S could be converted to E
1 in the aquatic environment, E
1 appears to be the most important natural endocrine disrupter.
In this study we parameterize the textural attributes of plagioclase phenocrysts and microlites from nineteen pyroclasts ejected during mild to violent explosions at Stromboli over a timespan of ...∼18 years, from 2003 to 2021. By allying kinetic and crystal size distribution principles, we document that the morphological stability of large-sized, euhedral phenocrysts is superimposed on an internal textural heterogeneity due to growth-dissolution phenomena associated with the input rate of hot, H2O-rich recharge magmas rising from depth. As a result, the volumetric plagioclase proportion, dominant size, and number of phenocrysts per unit volume decrease from mild to violent explosions responding to a more efficient magma mixing process via sustained injections of mafic magmas into the shallow reservoir. On the other hand, the crystallization of anhedral plagioclase microlites is controlled by fast growth kinetics taking place in the uppermost part of the conduit during magma acceleration towards the surface. Under such highly dynamic crystallization conditions, the microlite number density closely depends on the increase of melt liquidus temperature via magma decompression and H2O exsolution. This mutualism allows to model the degassing rate and ascent velocity of magma under open-conduit flow regimes for the different eruptive styles, thereby supporting the idea that violent explosions at Stromboli are driven by sustained influxes of recharge magmas favoring strong acceleration (∼12–27 m/s), decompression (∼0.25–0.49 MPa/s), and H2O exsolution (∼0.005–0.01 wt%/s) before magma discharge at the vent.
•Plagioclase crystallization dominates scoria clasts erupted at Stromboli.•Plagioclase textural attributes markedly change as a function of the eruptive style.•Phenocrysts record the renewal rate of the shallow magmatic reservoir.•Microlites reflect conduit dynamics before magma discharge at the vent.
Petrological studies of active volcanoes typically focus on eruptive phenomena occurring over long timescales of the order of days to years, aiming at identifying major changes in the ...physico-chemical state of magma during ascent towards the surface. Exceptionally, we present results from an integrated petrological and statistical approach based on the compilation of ∼5300 major and trace element data for glass and crystals, in combination with volcanological data on eruptive events occurred over timescales of minutes at Stromboli volcano (Sicily). On May 11, 2019, we had the rare opportunity to collect individual fresh fallout ash products from eighteen mostly consecutive explosions, erupted in a 2-h time span and, at the same time, to acquire continuous high frequency (50 Hz) infrared thermal data of the same explosions. Through video analysis, we observe that explosions were more frequent and ash-dominated at the southwestern crater area (SCA, 8–10 events/h) than at the northeastern crater area (NCA, 3–5 events/h), where coarser material was ejected. The statistical analysis of glass and plagioclase compositions reveals differences in the products erupted from the two crater areas. SCA explosions tapped less differentiated magmas (Mg#∼42–46, ∼257–365 LaN, ∼0.7–0.9 Eu/Eu*) in equilibrium with more anorthitic plagioclase cores (An∼72–88), whereas NCA area explosions are more differentiated (Mg#∼40–44, ∼286–387 LaN, ∼0.6–0.8 Eu/Eu*) and in equilibrium with less anorthitic plagioclase cores (An∼68–82). Thermometric calculations based on major and trace element clinopyroxene-plagioclase-melt equilibrium modeling highlight that the SCA explosions were statistically fed by hotter magmas in comparison to NCA explosions. Plagioclase-based diffusion modeling also indicates longer timescales for the dynamic ascent of NCA magmas, leading to preferential groundmass crystallization at the conduit walls and transition from sideromelane to tachylite groundmass textures. The final emerging picture is that in May 2019, concurrent normal eruptions from different crater areas at Stromboli were heralds of compositionally and thermally diverse magmas rising at different rates within the uppermost branched part of the conduit region. High frequency petrological investigations aided by statistical treatment of data have the potential to constrain dynamic conduit processes related to transient, explosive eruptions in persistently active volcanoes, thereby offering new insights on the interplay between magma dynamics, ascent timescales, and eruptive behavior.
•High-frequency eruption sampling at Stromboli reveals subtle chemical changes.•Sample compositions and textures differ between the two crater areas.•Thermal and compositional gradients are reconstructed in the shallow mush system.•Crystallization is promoted by slower magma ascent rates at the conduit walls.
Phonolite-trachyte associations are a common feature of alkaline volcanoes in intraplate settings, and their coexistence challenges closed-system magmatic differentiation scenarios. Here we have ...investigated the mineralogical and petrochemical features of dikes, lavas, pyroclastic deposits, and comagmatic crystal-rich enclaves outcropping at Dunedin Volcano (Otago region, southern New Zealand). These alkaline magmatic products show both highly and mildly alkaline affinities, trending towards phonolitic and trachytic end-members, respectively. Intermediate rocks are phonotephrites + tephriphonolites (highly alkaline series) and mugearites + benmoreites (mildly alkaline series) with a phenocryst assemblage of clinopyroxene + plagioclase ± amphibole formed at low to mid-crustal levels (i.e., ~29–16 km). Phonolites are porphyritic rocks characterized by alkali feldspar ± amphibole ± clinopyroxene. Their whole-rock compositions are highly enriched in incompatible elements, with variable Ba + Sr contents. A weak negative to slightly positive Eu anomaly is also associated with 87Sr/86Sr ratios of 0.7028–0.7031, which are comparable to those of parental magmas. Geochemical models indicate that phonolites originate as interstitial melts that are generated via abundant alkali feldspar crystallization from a shallow crystalline mush (i.e., ~14–5 km). Strong melt differentiation and extraction is testified by crystal-rich enclaves, as remnants of the mush region. On the other hand, trachytes are phenocryst-poor products strongly depleted in Ba + Sr and with a marked negative Eu anomaly. Trachytes are characterized by 87Sr/86Sr ratios of 0.7040–0.7060, which are different from intermediate rocks and phonolites, and trend towards crustal isotopic compositions. Integrated mass balance, trace element, and energy-constrained modeling confirm that trachytes originate from mildly alkaline magmas interacting with the country rock during feldspar fractionation. We interpret the transition from trachyte to phonolite formation and eruption resulting from the maturation of the plumbing system through accumulation, cooling, and degassing of both highly and mildly alkaline magmas.
•Thermobarometric constraints on the architecture of the magmatic plumbing system•Phonolites and trachytes differentiate from intermediate magmas in the shallow crust.•Trachytes originate by assimilation and fractional crystallization.•Phonolites represent interstitial liquids of upper crustal mush regions.•Trachyte to phonolite transition relates to maturity of the magmatic plumbing system.
Mafic alkaline magmas, such as those feeding the persistent eruptive activity of Stromboli and Mt. Etna volcanoes in Italy, are dominated by the crystallization of plagioclase via cooling and ...degassing phenomena related to the dynamics of shallow crustal reservoirs and eruptive conduits. Because plagioclase textures and compositions are extremely sensitive to the changes of intensive variables in subvolcanic plumbing systems, the phenomenological variability of erupted crystals preserves detailed evidence of complex growth histories. From this point of view, we reappraise the textural maturation and compositional complexity of plagioclase by allying thermodynamic and kinetic principles to natural and experimental observations, with the purpose of drawing up guidelines for reconstructing magma dynamics in mafic alkaline volcanic settings. A multifaceted statistical method is adopted to parameterize the decay of crystal growth rate with increasing crystallization time, as relaxation kinetics prevails over melt supersaturation effects. This model parameterization is combined with the textural analysis of natural plagioclase crystals to quantify the residence time of phenocrysts in equilibrium with magmas at Stromboli and Mt. Etna and/or the timescale of rapid microlite growth during disequilibrium ascent of magmas within the conduit. The role played by temperature and melt-water content on plagioclase components and major cation substitution mechanisms is also evaluated under both isobaric-isothermal and decompression conditions. The emerging paradigm is that the influence of dissolved water on anorthite-albite exchange between plagioclase and melt is overwhelmingly mitigated by changes in temperature at conditions of P = 30–300 MPa, T = 1050–1150 °C, fO2 = NNO + 1.9-NNO + 2.3, and melt-H2O = 0.6–4.4 wt%. As a corollary, anorthite and albite melt activities are almost fully encapsulated in the variation of anhydrous melt components as the crystallization of plagioclase proceeds during magma cooling. Following this line of reasoning, we propose an integrated modeling approach to decipher complex zoning patterns in natural plagioclase phenocrysts from mafic alkaline eruptions. Key findings from our re-assessment of equilibrium, thermometric, and hygrometric models indicate that temperature and dissolved water can be iteratively estimated for different plagioclase textural patterns if crystals are sufficiently strongly zoned and probability-based criteria are applied to determine the maximum probability distribution from kernel density analysis.