An official website of the United States government
Abstract Mineral chemistry records the pressure and temperature conditions of lithospheric processes. Active tectonic margins, however, are subjected to non‐hydrostatic stresses wherein stress magnitudes vary directionally, and the impact of non‐hydrostatic stress on mineral chemistry is uncertain. The work of materials scientists F. Larché and J. Cahn provides a framework for quantifying how stress affects mineral chemistry. Crystallographically and mechanically anisotropic, multicomponent minerals will have different compositions as a function of their orientation under a fixed stress meaning that grain‐to‐grain compositional variation can be used to estimate stress. We develop two “orientation piezometry” methods that use the chemistry and orientations of multicomponent, anisotropic minerals to estimate stress. The first method uses chemistry and orientation (“coupled orientation piezometry”) whereas the second method uses composition alone (“decoupled orientation piezometry”). We apply the methods to clinopyroxene and feldspar solid solutions using synthetic data sets. The first method determines the full stress tensor whereas the second method can only determine the differential stress magnitude unless additional a priori information is specified. Plausible scenarios for orientation piezometry include minerals undergoing diffusion creep, recrystallized grains formed during dislocation creep, and minerals grown statically under stress. Preliminary application of the decoupled piezometer to the famous eclogite facies shear zones on Holsnøy, Norway, suggests differential stresses in the range of 300–900 MPa, broadly consistent with previous estimates from the area. Thus, orientation piezometry techniques may provide valuable constraints on geodynamic processes and insights into long‐standing geological problems such as the relationship between pressure and depth.
more »
« less
Quantifying the Effects of Non‐Hydrostatic Stress on Multi‐Component Minerals
Abstract Mineral compositions are used to infer pressures, temperatures, and timescales of geological processes. The thermodynamic techniques underlying these inferences assume a uniform, constant pressure. Nonetheless, convergent margins generate significant non‐hydrostatic (unequal) stresses, violating the uniform pressure assumption and creating uncertainty. Materials scientists F. Larché and J. Cahn derived an equation suitable for non‐hydrostatically stressed geologic environments that links stress and equilibrium composition in elastic, multi‐component crystals. However, previous works have shown that for binary solid solutions with ideal mixing behavior, hundreds of MPa to GPa‐level stresses are required to shift mineral compositions by a few hundredths of a mole fraction, limiting the equation's applicability. Here, we apply Larché and Cahn's equation to garnet, clinopyroxene, and plagioclase solid solutions, incorporating for the first time non‐ideal mixing behavior and more than two endmembers. We show that non‐ideal mixing increases predicted stress‐induced composition changes by up to an order of magnitude. Further, incorporating additional solid solution endmembers changes the predicted stress‐induced composition shifts of the other endmembers being considered. Finally, we demonstrate that Larché and Cahn's approach yields positive entropy production, a requirement for any real process to occur. Our findings reveal that stresses between tens and a few hundred MPa can shift mineral compositions by several hundredths of a mole fraction. Consequently, mineral compositions could plausibly be used to infer stress states. We suggest that stress‐composition effects could develop via intracrystalline diffusion in any high‐grade metamorphic setting, but are most likely in hot, dry, and strong rocks such as lower crustal granulites.
more »
« less
- Award ID(s):
- 2208229
- PAR ID:
- 10373474
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 127
- Issue:
- 9
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Most theoretical descriptions of stresses induced by freezing are rooted in the (generalized) Clapeyron equation, which predicts the pressure that a solid can exert as it cools below its melting temperature. This equation is central for topics ranging beyond glaciology to geomorphology, civil engineering, food storage and cryopreservation. However, it has inherent limitations, requiring isotropic solid stresses and conditions near bulk equilibrium. Here, we examine when the Clapeyron equation is applicable by providing a rigorous derivation that details all assumptions. We demonstrate the natural extension for anisotropic stress states, and we show how the temperature and pressure ranges for validity depend on well-defined material properties. Finally, we demonstrate how the range of applicability of the (linear) Clapeyron equation can be extended by adding higher-order terms, yielding results that are in good agreement with experimental data for the pressure melting of ice.more » « less
-
Rhabdophane is a hydrous phosphate that commonly replaces monazite as a weathering product in critical mineral deposits during the alteration of rare earth elements (REE) bearing carbonatites and alkaline igneous complexes. It is an important host to the light (L)REE (i.e., La to Gd) but the stability and structure of binary solid solutions between the Ce and the other LREE endmembers have not yet been determined experimentally. Here we present room temperature calorimetric experiments that were used to measure the enthalpy of precipitation of rhabdophane (Ce1−xREExPO4·nH2O; REE = La, Pr, Nd, Sm, Eu, and Gd). The solids were characterized using X-ray diffraction, scanning electron microscopy, Raman spectroscopy, and the role of water in the rhabdophane structure was further determined using thermogravimetric analysis coupled with differential scanning calorimetry. The calorimetric experiments indicate a non-ideal behavior for all of the binary solid solutions investigated with an excess enthalpy of mixing (ΔHex) described by a 2- to 3-term Guggenheim parameters equation. The solid solutions were categorized into three groups: (1) binary Ce-La and Ce-Pr which display positive ΔHex values with a slight asymmetry; (2) binary Ce-Nd and Ce-Sm which display negative ΔHex values with a nearly symmetric shape; (3) Ce-Eu and Ce-Gd which display both negative and positive ΔHex values with nearly symmetric shape. The excess Gibbs energy (ΔGex) of the solid solutions was further investigated using a thermodynamic analysis approach of aqueous-solid solution equilibria and the optimization programs GEMS and GEMSFITS. The resulting ΔGex values combined with the calorimetric ΔHex values indicate that there is likely an excess entropy contribution implying important short-range structural modifications in the solid solutions dependent on the deviation of the REE ionic radii from the size of Ce3+. These observations corroborate with the trends in the Raman v1 stretching bands of the PO4-site. The excess molar volumes determined from X-ray diffraction analysis further indicate an overall asymmetric behavior in all of the studied binary solid solutions, which becomes increasingly important from La to Gd. The pronounced short-range order–disorder occurring in groups 2 and 3 solid solutions mimics some of the behavior observed from previous studies in anhydrous monazite solid solutions. This study highlights the potential to use the chemistry and the structural modifications of rhabdophane as potential indicators of formation conditions in geologic systems and permits improving our modeling capabilities of REE partitioning in critical minerals systems.more » « less
-
Abstract. We present new models for the activity of iron (γFe) in solid face-centered cubic (fcc) and liquid FePt alloy at high temperature and pressure to facilitate their use as sliding buffer redox sensors under extreme conditions. Numerous experimental studies of γFe in FePt alloy at 100 kPa have produced a wide spread of values. By favoring high-temperature studies that are more likely to have produced equilibrium measurement and excluding experiments for compositions and temperatures that probably encountered ordered or unmixed low-temperature phases, we regress an asymmetric Margules activity–composition model with parameters WFePtfcc=-121.5±2.1 kJ mol−1 and WPtFefcc=-93.3±4.3 kJ mol−1. These values are close to the widely used model of Kessel et al. (2001), but for Pt-rich compositions they predict larger Fe activities and correspondingly more reduced oxygen fugacities. Activity–composition relations in liquid FePt are calibrated from direct measurements of activities and, most sensitively, from the trace of the Fe–Pt liquidus. Together, these yield asymmetric Margules parameters of WFePtliq=-124.5 kJ mol−1 and WPtFeliq=-94.0 kJ mol−1. The effects of pressure on both fcc and liquid FePt alloy are considered from excess-volume relations. Both solid and liquid alloy display significant positive excess volumes of mixing. Extraction of the excess volume of mixing for fcc FePt alloy requires filtering data for ordered low-temperature phases and corrections for the effects of magnetostriction on Fe-rich compositions which exhibit “Invar” behavior. Applied at high temperatures and pressures, both solid and liquid FePt alloys have strongly negative deviations from ideality at low pressure, which become closer to ideal at high pressure. These models provide a provisional basis for the calculation of aFe in high-temperature, high-pressure experiments that, when combined with estimates of aFeO, allow characterization of fO2 under conditions relevant to magma oceans, core formation, and differentiation processes in the lower mantle of Earth or on other terrestrial planets. Improvements in these models require new constraints on the equation of state of FePt fcc alloy and documentation of the high-pressure melting relations in the system Fe–Pt.more » « less
-
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 estimatedP‐Tconditions. 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
