We provide a new algorithm for mass‐balance calculations in petrology and geochemistry based on the log‐ratio approach championed initially by John Aitchison (e.g., Aitchison, 1982,
To quantitatively convert upper mantle seismic wave speeds measured into temperature, density, composition, and corresponding and uncertainty, we introduce the
- Award ID(s):
- 1844340
- PAR ID:
- 10372952
- 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
-
more » « less
Abstract https://doi.org/10.1111/j.2517-6161.1982.tb01195.x ; Aitchison, 1984,https://doi.org/10.1007/bf01029316 ) along with the underlying principles, mathematical frameworks, and data requirements. Log‐ratio Inversion of Mixed End‐members (LIME) is written in MATLAB and calculates phase proportions in an experiment or rock given a bulk composition, the composition of each phase, and the associated compositional uncertainties. An important advantage of LIME is that performing the mass‐balance calculation in inverse log‐ratio space constrains phase proportions to be between 0 and 100 wt.%. Further, the resulting LIME phase proportions provide realistic estimates of uncertainty regardless of data distribution. These two characteristics of LIME improve upon standard multiple linear regression techniques, which may yield negative values for phase proportions if non‐constrained or report oversimplified symmetric errors. Primary applications of LIME include estimating phase abundances, calculating melting and metamorphic reaction stoichiometries, and checking for open system behavior in phase equilibria experiments. The technique presented here covers whole‐rock analysis, mineralogy, and phase abundance, but could be extended to isotopic tracers, trace element modeling, and regolith component un‐mixing. We highlight the importance of uncertainty estimations for phase abundances to the fields of petrology and geochemistry by comparing our results from LIME to previously published mass‐balance calculations. Furthermore, we present case studies that demonstrate the role of mass‐balance calculations in determining magma crystallinity and defining melting reactions. -
more » « less
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 beM 0 ≃τ 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
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
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.
