skip to main content


Title: WISTFUL: Whole‐Rock Interpretative Seismic Toolbox for Ultramafic Lithologies
Abstract

To quantitatively convert upper mantle seismic wave speeds measured into temperature, density, composition, and corresponding and uncertainty, we introduce theWhole‐rockInterpretativeSeismicToolboxForUltramaficLithologies (WISTFUL). WISTFUL is underpinned by a database of 4,485 ultramafic whole‐rock compositions, their calculated mineral modes, elastic moduli, and seismic wave speeds over a range of pressure (P) and temperature (T) (P = 0.5–6 GPa,T = 200–1,600°C) using the Gibbs free energy minimization routine Perple_X. These data are interpreted with a toolbox of MATLAB® functions, scripts, and three general user interfaces:WISTFUL_relations, which plots relationships between calculated parameters and/or composition;WISTFUL_geotherms, which calculates seismic wave speeds along geotherms; andWISTFUL_inversion, which inverts seismic wave speeds for best‐fit temperature, composition, and density. To evaluate our methodology and quantify the forward calculation error, we estimate two dominant sources of uncertainty: (a) the predicted mineral modes and compositions, and (b) the elastic properties and mixing equations. To constrain the first source of uncertainty, we compiled 122 well‐studied ultramafic xenoliths with known whole‐rock compositions, mineral modes, and estimatedPTconditions. We compared the observed mineral modes with modes predicted using five different thermodynamic solid solution models. The Holland et al. (2018,https://doi.org/10.1093/petrology/egy048) solution models best reproduce phase assemblages (∼12 vol. % phase root‐mean‐square error [RMSE]) and estimated wave speeds. To assess the second source of uncertainty, we compared wave speed measurements of 40 ultramafic rocks with calculated wave speeds, finding excellent agreement (VpRMSE = 0.11 km/s). WISTFUL easily analyzes seismic datasets, integrates into modeling, and acts as an educational tool.

 
more » « less
Award ID(s):
1844340
NSF-PAR ID:
10372952
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geochemistry, Geophysics, Geosystems
Volume:
23
Issue:
8
ISSN:
1525-2027
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    It is well known that elasticity is a key physical property in the determination of the structure and composition of the Earth and provides critical information for the interpretation of seismic data. This study investigates the stress‐induced variation in elastic wave velocities, known as the acoustoelastic effect, in San Carlos olivine. A recently developed experimental ultrasonic acoustic system, the Directly Integrated Acoustic System Combined with Pressure Experiments (DIASCoPE), was used with the D‐DIA multi‐anvil apparatus to transmit ultrasonic sound waves and collect the reflections. We use the DIASCoPE to obtain longitudinal (P) and shear (S) elastic wave velocities from San Carlos olivine at pressures ranging from 3.2–10.5 GPa and temperatures from 450–950°C which we compare to the stress state in the D‐DIA derived from synchrotron X‐ray diffraction. We use elastic‐plastic self‐consistent (EPSC) numerical modeling to forward model X‐ray diffraction data collected in D‐DIA experiments to obtain the macroscopic stress on our sample. We can observe the relationship between the relative elastic wave velocity change (ΔV/V) and macroscopic stress to determine the acoustoelastic constants, and interpret our observations using the linearized first‐order equation based on the model proposed by Hughes and Kelly (1953),https://doi.org/10.1103/physrev.92.1145. This work supports the presence of the acoustoelastic effect in San Carlos olivine, which can be measured as a function of pressure and temperature. This study will aid in our understanding of the acoustoelastic effect and provide a new experimental technique to measure the stress state in elastically deformed geologic materials at high pressure conditions.

     
    more » « less
  2. Abstract

    We address in situ serpentinization and mineral carbonation processes in oceanic lithosphere using integrated field magnetic measurements, rock magnetic analyses, superconducting quantum interference device (SQUID) microscopy, microtextural observations, and energy dispersive spectroscopy phase mapping. A representative suite of ultramafic rock samples were collected, within the Atlin ophiolite, along a 100‐m long transect across a continuous outcrop of mantle harzburgite with several alteration fronts: serpentinite, soapstone (magnesite + talc), and listvenite (magnesite + quartz). Strong correlations between changes in magnetic signal strengths and amount of alteration are shown with distinctive contrasts between serpentinite, transitional soapstone, and listvenite that are linked to the formation and breakdown of magnetite. While previous observations of the Linnajavri ultramafic complex indicated that the breakdown of magnetite occurred during listvenite formation from the precursor soapstone (Tominaga et al., 2017,https://doi.org/10.1038/s41467-017-01610-4), results from our study suggest that magnetite destabilization already occurred during the replacement of serpentinite by soapstone (i.e., at lower fluid CO2concentrations). This difference is attributed to fracture‐controlled flow of sulfur‐bearing alteration fluid at Atlin, causing reductive magnetite dissolution in thin soapstone zones separating serpentinite from sulfide‐mineralized listvenite. We argue that magnetite growth or breakdown in soapstone provides insight into the mode of fluid flow and the composition, which control the scale and extent of carbonation. This conclusion enables us to use magnetometry as a viable tool for monitoring the reaction progress from serpentinite to carbonate‐bearing assemblages in space and time with a caution that the three‐dimensionality of magnetic sources impacts the scalability of measurements.

     
    more » « less
  3. Abstract

    Mineral inclusions are ubiquitous in metamorphic rocks and elastic models for host‐inclusion pairs have become frequently used tools for investigating pressure–temperature (P–T) conditions of mineral entrapment. Inclusions can retain remnant pressures () that are relatable to their entrapmentP–Tconditions using an isotropic elastic model andP–T–Vequations of state for host and inclusion minerals. Elastic models are used to constrainP–Tcurves, known as isomekes, which represent the possible inclusion entrapment conditions. However, isomekes require a temperature estimate for use as a thermobarometer. Previous studies obtained temperature estimates from thermometric methods external of the host‐inclusion system. In this study, we present the firstP–Testimates of quartz inclusion entrapment by integrating the quartz‐in‐garnet elastic model with titanium concentration measurements of inclusions and a Ti‐in‐quartz solubility model (QuiG‐TiQ). QuiG‐TiQ was used to determine entrapmentP–Tconditions of quartz inclusions in garnet from a quartzofeldspathic gneiss from Goodenough Island, part of the (ultra)high‐pressure terrane of Papua New Guinea. Raman spectroscopic measurements of the 128, 206, and 464 cm−1bands of quartz were used to calculate inclusion pressures using hydrostatic pressure calibrations (), a volume strain calculation (), and elastic tensor calculation (), that account for deviatoric stress.values calculated from the 128, 206, and 464 cm−1bands’ hydrostatic calibrations are significantly different from one another with values of 1.8 ± 0.1, 2.0 ± 0.1, and 2.5 ± 0.1 kbar, respectively. We quantified elastic anisotropy using the 128, 206 and 464 cm−1Raman band frequencies of quartz inclusions and stRAinMAN software (Angel, Murri, Mihailova, & Alvaro, 2019, 234:129–140). The amount of elastic anisotropy in quartz inclusions varied by ~230%. A subset of inclusions with nearly isotropic strains gives an averageandof 2.5 ± 0.2 and 2.6 ± 0.2 kbar, respectively. Depending on the sign and magnitude, inclusions with large anisotropic strains respectively overestimate or underestimate inclusion pressures and are significantly different (<3.8 kbar) from the inclusions that have nearly isotropic strains. Titanium concentrations were measured in quartz inclusions exposed at the surface of the garnet. The average Ti‐in‐quartz isopleth (19 ± 1 ppm [2σ]) intersects the average QuiG isomeke at 10.2 ± 0.3 kbar and 601 ± 6°C, which are interpreted as theP–Tconditions of quartzofeldspathic gneiss garnet growth and entrapment of quartz inclusions. TheP–Tintersection point of QuiG and Ti‐in‐quartz univariant curves represents mechanical and chemical equilibrium during crystallization of garnet, quartz, and rutile. These three minerals are common in many bulk rock compositions that crystallize over a wide range ofP–Tconditions thus permitting application of QuiG‐TiQ to many metamorphic rocks.

     
    more » « less
  4. Abstract

    Elevated seismic noise for moderate‐size earthquakes recorded at teleseismic distances has limited our ability to see their complexity. We develop a machine‐learning‐based algorithm to separate noise and earthquake signals that overlap in frequency. The multi‐task encoder‐decoder model is built around a kernel pre‐trained on local (e.g., short distances) earthquake data (Yin et al., 2022,https://doi.org/10.1093/gji/ggac290) and is modified by continued learning with high‐quality teleseismic data. We denoise teleseismic P waves of deep Mw5.0+ earthquakes and use the clean P waves to estimate source characteristics with reduced uncertainties of these understudied earthquakes. We find a scaling of moment and duration to beM0 ≃ τ4, and a resulting strong scaling of stress drop and radiated energy with magnitude ( and ). The median radiation efficiency is 5%, a low value compared to crustal earthquakes. Overall, we show that deep earthquakes have weak rupture directivity and few subevents, suggesting a simple model of a circular crack with radial rupture propagation is appropriate. When accounting for their respective scaling with earthquake size, we find no systematic depth variations of duration, stress drop, or radiated energy within the 100–700 km depth range. Our study supports the findings of Poli and Prieto (2016,https://doi.org/10.1002/2016jb013521) with a doubled amount of earthquakes investigated and with earthquakes of lower magnitudes.

     
    more » « less
  5. The granitic water-saturated solidus (G-WSS) is the lower temperature limit of magmatic mineral crystallization. The accepted water-saturated solidus for granitic compositions was largely determined >60 years ago1. More recent advances in experimental petrology, improved analytical techniques, and recent observations that granitic systems can remain active or spend a significant proportion of their lives at conditions below the traditional G-WSS2–5 necessitate a careful experimental investigation of the near-solidus regions of granitic systems. Natural and synthetic starting materials were melted at 10 kbar and 900°C with 48 wt% H2O to produce hydrous glasses for subsequent experiments at lower PT conditions used to locate the G-WSS. We performed crystallization experiments and melting experiments at temperatures ranging from 575 to 800°C and 1, 6, 8, and 10 kbar on 12 granitoid compositions. First, we ran a series of isothermal crystallization experiments along each isobar at progressively lower temperatures until runs completely crystallized to identify apparent solidus temperatures. Geochemical analyses of quenched glass compositions demonstrate that progressive crystallization drives all starting compositions towards silica-rich, water-saturated rhyolitic/granitic melts (e.g., ~7578 wt% SiO2). After identifying the apparent solidus temperatures at which the various compositions crystallized, we then ran series of reversal-type melting experiments. With the goal of producing rocks with hydrous equilibrium microstructures, we crystallized compositions at temperatures ~10°C below the apparent solidus identified in crystallization experiments, and then heated isobarically to conditions that produced ~20% melt during the crystallization experiments. Importantly, crystallization experiments and heating experiments at the same PT conditions produced similar proportions of melt, crystals, and vapor. A time-series of experiments 230 days at PT conditions previously identified to produce ~10% to 20% melt did not reveal any kinetic effects on melt crystallization. Experiments at 6 to 10 kbar crystallized/melted at temperatures close to the published G-WSS. However, at lower pressures where the published G-WSS is strongly curved in PT space, all compositions investigated contained melt to temperatures ~75 to 100°C below the accepted G-WSS. The similarity of crystallization temperatures for the higher-pressure experiments to previously published results, similar phase proportions in melting and crystallization experiments, and the lack of kinetic effects on crystallization collectively suggest that our lower pressure constraints on the G-WSS are accurate. The new experimental results demonstrating that the lower-pressure G-WSS is significantly lower than unanimously accepted estimates will help us to better understand the storage conditions, evolution, and potential for eruption in mid- to upper-crustal silicic magmatic systems. (1) Tuttle, O.; Bowen, N. Origin of Granite in the Light of Experimental Studies in the System NaAlSi3O8–KAlSi3O8–SiO2–H2O; Geological Society of America Memoirs; Geological Society of America, 1958; Vol. 74. https://doi.org/10.1130/MEM74. (2) Rubin, A. E.; Cooper, K. M.; Till, C. B.; Kent, A. J. R.; Costa, F.; Bose, M.; Gravley, D.; Deering, C.; Cole, J. Rapid Cooling and Cold Storage in a Silicic Magma Reservoir Recorded in Individual Crystals. Science 2017, 356 (6343), 1154–1156. https://doi.org/10.1126/science.aam8720. (3) Andersen, N. L.; Jicha, B. R.; Singer, B. S.; Hildreth, W. Incremental Heating of Bishop Tuff Sanidine Reveals Preeruptive Radiogenic Ar and Rapid Remobilization from Cold Storage. Proceedings of the National Academy of Sciences 2017, 114 (47), 12407–12412. https://doi.org/10.1073/pnas.1709581114. (4) Ackerson, M. R.; Mysen, B. O.; Tailby, N. D.; Watson, E. B. Low-Temperature Crystallization of Granites and the Implications for Crustal Magmatism. Nature 2018, 559 (7712), 94–97. https://doi.org/10.1038/s41586-018-0264-2. (5) Glazner, A. F.; Bartley, J. M.; Coleman, D. S.; Lindgren, K. Aplite Diking and Infiltration: A Differentiation Mechanism Restricted to Plutonic Rocks. Contributions to Mineralogy and Petrology 2020, 175 (4). https://doi.org/10.1007/s00410-020-01677-1. 
    more » « less