skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Evaluating the Role of Titanomagnetite in Bubble Nucleation: Novel Applications of Low Temperature Magnetic Analysis and Textural Characterization of Rhyolite Pumice and Obsidian From Glass Mountain, California
Nucleation of H2O vapor bubbles in magma requires surpassing a chemical supersaturation threshold via decompression. The threshold is minimized in the presence of a nucleation substrate (heterogeneous nucleation, <50 MPa), and maximized when no nucleation substrate is present (homogeneous nucleation, >100 MPa). The existence of explosively erupted aphyric rhyolite magma staged from shallow (<100 MPa) depths represents an apparent paradox that hints at the presence of a cryptic nucleation substrate. In a pair of studies focusing on Glass Mountain eruptive units from Medicine Lake, California, we characterize titanomagnetite nanolites and ultrananolites in pumice, obsidian, and vesicular obsidian (Brachfeld et al., 2024,https://doi.org/10.1029/2023GC011336), calculate titanomagnetite crystal number densities, and compare titanomagnetite abundance with the physical properties of pumice to evaluate hypotheses on the timing of titanomagnetite crystallization. Titanomagnetite crystals with grain sizes of approximately 3–33 nm are identified in pumice samples from the thermal unblocking of low‐temperature thermoremanent magnetization. The titanomagnetite number densities for pumice are 10^18 to 10^20 m^−3, comparable to number densities in pumice and obsidian obtained from room temperature methods (Brachfeld et al., 2024,https://doi.org/10.1029/2023GC011336'>https://doi.org/10.1029/2023GC011336). This range exceeds reported bubble number densities (BND) within the pumice from the same eruptive units (average BND ∼4 × 10^14 m^−3). The similar abundances of nm‐scale titanomagnetite crystals in the effusive and explosive products of the same eruption, together with the lack of correlation between pumice permeability and titanomagnetite content, are consistent with titanomagnetite formation having preceded the bubble formation. Results suggest sub‐micron titanomagnetite crystals are responsible for heterogeneous bubble nucleation in this nominally aphyric rhyolite magma.  more » « less
Award ID(s):
1839230 1839313 0521069 0948262 2153786
PAR ID:
10502181
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
American Geophysical Union
Date Published:
Journal Name:
Geochemistry, Geophysics, Geosystems
Volume:
25
Issue:
4
ISSN:
1525-2027
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Nanometer‐scale titanomagnetite crystals have been detected in nominally aphyric rhyolite pumice, but whether they are numerous enough to impact bubble nucleation in explosive silicic volcanism was unresolved. This study examines sub‐micron crystals using rock magnetic techniques, Rhyolite‐MELTS modeling, and physical characterization. We analyzed pumice from four eruptions spanning wide ranges in intensity, storage depth, and bubble number density (1016to 1013 m−3liquid): 1060 CE Glass Mountain, 1912 CE Novarupta, 232 CE Taupo, and 0.45 Ma Pudahuel. Calculations assuming monospecific assemblages of 10 and 1,000 nm cubic particles yield titanomagnetite number densities of 1021to 1013m−3dense rock equivalent, respectively. In all cases, titanomagnetite is thermodynamically stable at pre‐eruptive storage conditions and magnetic susceptibility (χLF) is independent of vesicularity and permeability, indicating that crystals likely formed prior to vesiculation. The existence of nm‐scale Fe‐Ti oxides in four diverse cases suggests that heterogeneous bubble nucleation is a general feature of explosive rhyolite volcanism. 
    more » « less
  2. This dataset archived with the Magnetics Information Consortium contains rock-magnetic data for rhyolitic pumice and obsidian from Glass Mountain, Medicine Lake, California, USA. Data were generated at Montclair State University and include magnetic susceptibility measured at 976Hz and 3904Hz, magnetic susceptibility vs. temperature, anhysteretic remanent magnetization (ARM), and magnetic hysteresis measurements. This dataset accompanies the publication Brachfeld, S., McCartney, K., Hammer, J.E., Shea, T., Giachetti, T., Evaluating the role of titanomagnetite in bubble nucleation: Rock magnetic detection and characterization of nanolites and ultra-nanolites in rhyolite pumice and obsidian from Glass Mountain, California, Geochemistry Geophysics Geosystems, https://doi.org/10.1029/2023GC011336. 
    more » « less
  3. We document the presence, composition, and number density (TND) of titanomagnetite nanolites and ultra‐nanolites in aphyric rhyolitic pumice, obsidian, and vesicular obsidian from the 1060 CE Glass Mountain volcanic eruption of Medicine Lake Volcano, California, using magnetic methods. Curie temperatures indicate compositions of Fe2.40Ti0.60O4 to Fe3O4. Rock‐magnetic parameters sensitive to domain state, which is dependent on grain volume, indicate a range of particle sizes spanning superparamagnetic (<50–80 nm) to multidomain (>10 μm) particles. Cylindrical cores drilled from the centers of individual pumice clasts display anisotropy of magnetic susceptibility with prolate fabrics, with the highest degree of anisotropy coinciding with the highest vesicularity. Fabrics within a pumice clast require particle alignment within a fluid, and are interpreted to result from the upward transport of magma driven by vesiculation, ensuing bubble growth, and shearing in the conduit. Titanomagnetite number density (TND) is calculated from titanomagnetite volume fraction, which is determined from ferromagnetic susceptibility. TND estimates for monospecific assemblages of 1,000 nm–10 nm cubes predict 10^12 to 10^20 m^−3 of solid material, respectively. TND estimates derived using a power law distribution of grain sizes predict 10^18 to 10^19  m^−3. These ranges agree well with TND determinations of 10^18 to 10^20  m^−3 made by McCartney et al. (2024), and are several orders of magnitude larger than the number density of bubbles in these materials. These observations are consistent with the hypothesis that titanomagnetite crystals already existed in extremely high number‐abundance at the time of magma ascent and bubble nucleation. 
    more » « less
  4. Abstract In their article entitled “Trapdoor Fault Activation:A Step Toward Caldera Collapse at Sierra Negra,Galapagos,Ecuador” Shreve and Delgado (2023,https://doi.org/10.1029/2023jb026437) examine co‐eruptive deformation during the 2018 eruption of Sierra Negra Volcano. One of their major conclusions is that the 2018 eruption, and specifically co‐eruptive faulting, represents the initial stages of caldera collapse. They reach this conclusion because they focus their analysis solely on co‐eruptive deformation, and do not investigate the total (net) deformation for the 2005 to 2018 eruption cycle. Bell, La Famina, et al. (2021,https://doi.org/10.1038/s41467‐021‐21596‐4) investigated both the pre‐ and co‐eruptive phases of the 2018 eruption and showed that net deformation was one of caldera resurgence, not subsidence. In this comment, we demonstrate that the conclusion of collapse, or even initiation of collapse, is attributable to not accounting for pre‐eruptive deformation on the intra‐caldera Trapdoor Fault system and incorrectly assuming that the volcano‐tectonic dynamics of Sierra Negra mimic those of other basaltic calderas. 
    more » « less
  5. Abstract We examine the behavior of natural basaltic and trachytic samples during paleointensity experiments on both the original and laboratory‐acquired thermal remanences and characterize the samples using proxies for domain state including curvature (k) and the bulk domain stability parameters of Paterson (2011,https://doi.org/10.1029/2011JB008369) and Paterson et al. (2017,https://doi.org/10.1073/pnas.1714047114), respectively. A curvature value of 0.164 (suggested by Paterson, 2011,https://doi.org/10.1029/2011JB008369) as a critical threshold that separates single‐domain‐like remanences from multidomain‐like remanances on the original paleointensity data was used to separate samples into “straight” (single‐domain‐like) and “curved” (multidomain‐like) groups. Specimens from the two sample sets were given a “fresh” thermal remanent magnetization in a 70 μT field and subjected to an infield‐zerofield, zerofield‐infield (IZZI)‐type (Yu et al., 2004,https://doi.org/10.1029/2003GC000630) paleointensity experiment. The straight sample set recovered the laboratory field with high precision while the curved set had much more scattered results (70.5 ± 1.5 and 71.9 ± 5.2 μT, respectively). The average intensity of both sets for straight and curved was quite close to the laboratory field of 70 μT, however, suggesting that if experiments contain a sufficient number of specimens, there does not seem to be a large bias in the field estimate. We found that the dependence of the laboratory thermal remanent magnetization on cooling rate was significant in most samples and did not depend on domain states inferred from proxies based on hysteresis measurements and should be estimated for all samples whose cooling rates differ from that used in the laboratory. 
    more » « less