An official website of the United States government
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
Quantifying the Effects of Non‐Hydrostatic Stress on Single‐Component Polymorphs
Abstract Gibbs free energy, the fundamental thermodynamic potential used to calculate equilibrium mineral assemblages in geological systems, does not apply to non‐hydrostatically stressed solids. Consequently, there is debate over the significance of non‐hydrostatic stress in petrological and geophysical processes. To help resolve this debate, we consider the effects of non‐hydrostatic stress on the polymorph pairs kyanite/sillimanite, graphite/diamond, calcite/aragonite, and quartz/coesite. While these polymorphs are most relevant to metamorphic processes, the concepts developed are applicable to any single‐component solid reaction. We quantitatively show how stress variations normal to an interface alter equilibrium temperatures of polymorph pairs by approximately two orders of magnitude more than stress variations parallel to an interface. Thus, normal stress controls polymorph stability to first order. High‐pressure polymorphs are expected to preferentially nucleate normal to and grow parallel to the maximum stress and low‐pressure polymorphs, the minimum stress. Nonetheless, stress variations parallel to an interface allow for the surprising possibility that a high‐pressure polymorph can become more stable relative to a low‐pressure polymorph as stress decreases. The effects of non‐hydrostatic stress on mineral equilibrium are unlikely to be observed in systems with interconnected, fluid‐filled porosity, as fluid‐mediated reactions yield mineral assemblages at approximately constant pressures. In dry systems, however, reactions can occur directly between elastic solids, facilitating the direct application of non‐hydrostatic thermodynamics. Non‐hydrostatic stress is likely to be important to the evolution of metamorphic systems, as preferential orientations of polymorphic reactions can generate seismicity and may influence fundamental rock properties such as porosity and seismic anisotropy.
more »
« less
- Award ID(s):
- 2208229
- PAR ID:
- 10633571
- Publisher / Repository:
- Journal of Geophysical Research: Solid Earth
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 126
- Issue:
- 5
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Porosity generated during fluid–rock reaction can facilitate fluid transport and metasomatism in low permeability high-pressure metamorphic rocks. Evidence for reaction-induced porosity is found in an eclogite-facies clinopyroxene + apatite vein in an undeformed eclogitized Fe–Ti metagabbro from the Monviso Ophiolite (W. Alps) with a distinct garnet-rich selvage. Vein-forming fluids were sourced from adjacent metagabbros and reaction with the host rock removed Ca and P from the selvage and added Fe, REE, Pb and Cr. Textures at the selvage–host rock interface and in the host rock record local heterogeneity in reactivity and porosity during metasomatism linked to variable initial lawsonite abundance. These features reflect a hierarchy of pervasive-to-channelized porosity structures that facilitated widespread metasomatism of the host rock. Development of this metasomatic system in response to locally derived fluids suggests large-scale externally derived fluid transport is not required to drive extensive fluid–rock exchange. The production of porosity during metasomatic reactions could be important in facilitating further fluid–rock reaction and fluid transport in subducting slabs where permeability is low.more » « less
-
Abstract. Modelling the pressure in the Earth's interior is a common problem in Earth sciences. In this study we propose a method based on the conservation of the momentum of a fluid by using a hydrostatic scenario or a uniformly moving fluid to approximate the pressure. This results in a partial differential equation (PDE) that can be solved using classical numerical methods. In hydrostatic cases, the computed pressure is the lithostatic pressure. In non-hydrostatic cases, we show that this PDE-based approach better approximates the total pressure than the classical 1D depth-integrated approach. To illustrate the performance of this PDE-based formulation we present several hydrostatic and non-hydrostatic 2D models in which we compute the lithostatic pressure or an approximation of the total pressure, respectively. Moreover, we also present a 3D rift model that uses that approximated pressure as a time-dependent boundary condition to simulate far-field normal stresses. This model shows a high degree of non-cylindrical deformation, resulting from the stress boundary condition, that is accommodated by strike-slip shear zones. We compare the result of this numerical model with a traditional rift model employing free-slip boundary conditions to demonstrate the first-order implications of considering “open” boundary conditions in 3D thermo-mechanical rift models.more » « less
-
Bone mineralization is critical to maintaining tissue mechanical function. The application of mechanical stress via exercise promotes bone mineralization via cellular mechanotransduction and increased fluid transport through the collagen matrix. However, due to its complex composition and ability to exchange ions with the surrounding body fluids, bone mineral composition and crystallization is also expected to respond to stress. Here, a combination of data from materials simulations, namely density functional theory and molecular dynamics, and experimental studies were input into an equilibrium thermodynamic model of bone apatite under stress in an aqueous solution based on the theory of thermochemical equilibrium of stressed solids. The model indicated that increasing uniaxial stress induced mineral crystallization. This was accompanied by a decrease in calcium and carbonate integration into the apatite solid. These results suggest that weight-bearing exercises can increase tissue mineralization via interactions between bone mineral and body fluid independent of cell and matrix behaviours, thus providing another mechanism by which exercise can improve bone health. This article is part of a discussion meeting issue ‘Supercomputing simulations of advanced materials’.more » « less
-
This study discusses turbulent suspension flows of non-Brownian, non-colloidal, neutrally buoyant and rigid spherical particles in a Newtonian fluid over porous media with particles too large to penetrate and move through the porous layer. We consider suspension flows with the solid volume fraction $${{\varPhi _b}}$$ ranging from 0 to 0.2, and different wall permeabilities, while porosity is constant at 0.6. Direct numerical simulations with an immersed boundary method are employed to resolve the particles and flow phase, with the volume-averaged Navier–Stokes equations modelling the flow within the porous layer. The results show that in the presence of particles in the free-flow region, the mean velocity and the concentration profiles are altered with increasing porous layer permeability because of the variations in the slip velocity and wall-normal fluctuations at the suspension-porous interface. Furthermore, we show that variations in the stress condition at the interface significantly affect the particle near-wall dynamics and migration toward the channel core, thereby inducing large modulations of the overall flow drag. At the highest volume fraction investigated here, $${{\varPhi _b}}= 0.2$$ , the velocity fluctuations and the Reynolds shear stress are found to decrease, and the overall drag increases due to the increase in the particle-induced stresses.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1029/2020JB021594
