Bubble nucleation is the critical first step during magma degassing. The resultant number density of bubbles provides a record of nucleation kinetics and underlying eruptive conditions. The rate of bubble nucleation is strongly dependent on the surface free energy associated with nucleus formation, making the use of bubble number density for the interpretation of eruptive conditions contingent upon a sound understanding of surface tension. Based on a suite of nucleation experiments with up to >1016bubbles per unit volume of melt, and using numerical simulations of bubble nucleation and growth during each experiment, we provide a new formulation for surface tension during homogeneous nucleation of H2O bubbles in rhyolitic melt. It is based on the Tolman correction with a Tolman length of
Magnetite (Fe3O4) is an iron ore mineral that is globally mined especially for steel production. It is denser (5.15 g/cm3) than Earth’s crust (~2.7 g/cm3) and is expected to accumulate at the bottom of melt-rich magma reservoirs. However, recent studies revealed heterogeneous fluid bubble nucleation on oxide minerals such as magnetite during fluid degassing in volcanic systems. To test if the attachment on fluid bubbles is strong enough to efficiently float magnetite in silicate magma, decompression experiments were conducted at geologically relevant magmatic conditions with subsequent annealing to simulate re-equilibration after decompression. The results demonstrate that magnetite-bubble pairs do ascend in silicate melt, accumulating in an upper layer that grows during re-equilibration. This outcome contradicts the paradigm that magnetite must settle gravitationally in silicate melt.
more » « less- NSF-PAR ID:
- 10153613
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 9
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract δ = 0.32 nm, which implies an increase in surface tension of bubbles with decreasing nucleus size. Our model results indicate that experiments encompass two distinct nucleation regimes, distinguishable by the ratio of the characteristic diffusion time of water,τ diff, to the decompression time,t d. Experiments with >1013 m−3bubbles are characterized byτ diff/t d≪ 1, wherein the nucleation rate predominantly depends on the interplay between decompression and diffusion rates. Nucleation occurs over a short time interval with nucleation rate peaks at high values. For experiments with comparatively low bubble number density the average distance between adjacent bubbles and the diffusion timescale is large. Consequently,τ diff/t d≫ 1 and nucleation is nearly unaffected by diffusion and independent of decompression rate, with bubbles nucleating at an approximately constant rate until the sample is quenched. -
more » « less
Abstract Despite its increasing application to estimate magma decompression rates for explosive eruptions, the embayment speedometer has long awaited critical experimental evaluation. We present the first experimental results on the fidelity of natural quartz‐hosted embayments in rhyolitic systems as recorders of magma decompression. We conducted two high pressure‐temperature isobaric equilibrium experiments and 13 constant‐rate, continuous isothermal decompression experiments in a cold‐seal pressure vessel where we imposed rates from 0.005 to 0.05 MPa s−1in both H2O‐saturated and mixed‐volatile (H2O + CO2)‐saturated systems. In both equilibrium experiments, we successfully re‐equilibrated embayment melt to new fluid compositions at 780°C and 150 MPa, confirming the ability of embayments to respond to and record changing environmental conditions. Of the 32 glassy embayments recovered, seven met the criteria previously established for successful geospeedometry and were thus analyzed for their volatile (H2O ± CO2) concentrations, with each producing a good model fit and recovering close to the imposed decompression rate. In one H2O‐saturated experiment, modeling H2O concentration gradients in embayments from three separate crystals resulted in best‐fit decompression rates ranging from 0.012 to 0.021 MPa s−1, in close agreement with the imposed rate (0.015 MPa s−1) and attesting to the reproducibility of the technique. For mixed‐volatile experiments, we found that a slightly variable starting fluid composition (2.4–3.5 wt.% H2O at 150 MPa) resulted in good fits to both H2O + CO2profiles. Overall our experiments provide confidence that the embayment is a robust recorder of constant‐rate, continuous decompression, with the model successfully extracting experimental conditions from profiles representing nearly an order of magnitude variation (0.008–0.05 MPa s−1) in decompression rate.
-
more » « less
Abstract Magma from Plinian volcanic eruptions contains an extraordinarily large numbers of bubbles. Nucleation of those bubbles occurs because pressure decreases as magma rises to the surface. As a consequence, dissolved magmatic volatiles, such as water, become supersaturated and cause bubbles to nucleate. At the same time, diffusion of volatiles into existing bubbles reduces supersaturation, resulting in a dynamical feedback between rates of nucleation due to magma decompression and volatile diffusion. Because nucleation rate increases with supersaturation, bubble number density (BND) provides a proxy record of decompression rate, and hence the intensity of eruption dynamics. Using numerical modeling of bubble nucleation, we reconcile a long-standing discrepancy in decompression rate estimated from BND and independent geospeedometers. We demonstrate that BND provides a record of the time-averaged decompression rate that is consistent with independent geospeedometers, if bubble nucleation is heterogeneous and facilitated by magnetite crystals.
-
more » « less
Abstract Interpreting unrest at silicic volcanoes requires knowledge of the magma storage conditions and dynamics that precede eruptions. The Laguna del Maule volcanic field, Chile, has erupted ~40 km3of rhyolite over the last 20 ka. Astonishing rates of sustained surface inflation at >25 cm/year for >12 years reveal a large, restless system. Integration of geochronologic, petrologic, geomorphic, and geophysical observations provides an unusually rich context to interpret ongoing and prehistoric processes. We present new volatile (H2O, CO2, S, F, and Cl), trace element, and major element concentrations from 109 melt inclusions hosted in quartz, plagioclase, and olivine from seven eruptions. Silicic melts contain up to 8.0 wt. % H2O and 570 ppm CO2. In rhyolites melt inclusions track decompression‐driven fractional crystallization as magma ascended from ~14 to 4 km. This mirrors teleseismic tomography and magnetotelluric findings that reveal a domain containing partial melt spanning from 14 to 4 km. Ce and Cl contents of rhyolites support the generation of compositionally distinct domains of eruptible rhyolite within the larger reservoir. Heat, volatiles, and melt derived from episodic mafic recharge likely incubate and grow the shallow reservoir. Olivine‐hosted melt inclusions in mafic tephra contain up to 2.5 wt. % H2O and 1,140 ppm CO2and proxy for the volatile load delivered via recharge into the base of the silicic mush at ~14 to 8 km. We propose that mafic recharge flushes deeper reaches of the magma reservoir with CO2that propels H2O exsolution, upward accumulation of fluid, pressurization, and triggering of rhyolitic eruptions.
-
null (Ed.)Magnetite is the most important iron ore in iron oxide-apatite (IOA) deposits which represent the Cu-poor endmember of the iron oxide-copper–gold (IOCG) clan. Magnetite chemistry has been used as a petrogenetic indicator to identify the geological environment of ore formation and as a fingerprint of the source reservoir of iron. In this study, we present new textural and microanalytical EPMA and LA-ICP-MS data of magnetite from Carmen, Fresia, Mariela and El Romeral IOA deposits located in the Cretaceous Coastal Cordillera of northern Chile. We also provide a comprehensive summary and discussion of magnetite geochemistry from Andean IOAs including Los Colorados, Cerro Negro Norte, El Romeral (Chilean Iron Belt) and the Pliocene El Laco IOA deposit located in the Central Volcanic Zone of the Chilean Andes. Microtextures coupled with geochemical data were used to define and characterize the occurrence of different magnetite types. Magnetite exhibits a variety of textural features including oscillatory zoning, colloform banding, re-equilibration textures, exsolution lamellae and symplectites. The magmatic vs. hydrothermal origin of the different magnetite types and the evolution of IOA deposits can be assessed using diagrams based on compatible trace elements. However, magnetite is very susceptible to hydrothermal alteration and to both textural and compositional re-equilibration during magmatic and superimposed hydrothermal events. Based on the data presented here, we conclude that V and Ga are possibly the most reliable compatible elements in magnetite to trace ore-forming processes in the Andean IOA deposits. Magnetite chemistry reveals different conditions/events of formation for each IOA deposit ranging from high-temperature, low-oxygen fugacity (ƒO2), purely magmatic (> 600 °C) conditions; to lower temperature and higher ƒO2 magmatic-hydrothermal (300–600 °C) to low-temperature hydrothermal (< 200–300 °C) conditions. Specifically, a continuous transition from high-temperature, low- ƒO2 conditions in the deepest portions of the deposits to low-temperature, relatively higher ƒO2 conditions towards surface are described for magnetite from El Laco. The new and compiled magnetite data from IOA deposits from the Chilean Iron Belt and El Laco are consistent with a transition from magmatic to hydrothermal conditions. The flotation model plausibly explains such features, which result from the crystallization of magnetite microlites from a silicate melt, nucleation and coalescence of aqueous fluid bubbles on magnetite surfaces, followed by ascent of a fluid-magnetite suspension along reactivated transtensional faults or through fissures formed during the collapse of the volcanic structure (El Laco). The decompression of the coalesced fluid-magnetite aggregates during ascent promotes the continued growth of magnetite microlites from the Fe-rich magmatic-hydrothermal fluid. As with any general genetic model, the flotation model allows variation and the definition of different styles or subtypes of IOA mineralization. The deeper, intrusive-like Los Colorados deposit shows contrasting features when compared with the Cerro Negro Norte hydrothermal type, the pegmatitic apatite-rich deposits of Carmen, Fresia and Mariela, and the shallow, subaerial deposits of El Laco. These apparent differences depend fundamentally on the depth of formation, the presence of structures and faults that trigger decompression, the composition of the host rocks, and the source and flux rate of hydrothermal fluids.more » « less
