An official website of the United States government
Hydrogen is a rapidly diffusing monovalent cation in nominally anhydrous minerals (NAMs, such as olivine, orthopyroxene, and clinopyroxene), which is potentially re-equilibrated during silicate melt-rock and aqueous fluid-rock interactions in massif and abyssal peridotites. We apply a 3D numerical diffusion modeling technique to provide first-order timescales of complete hydrogen re-equilibration in olivine, clinopyroxene, and orthopyroxene over the temperature range 600-1200°C. Model crystals are 1-3 mm along the c-axis and utilize H+ diffusion coefficients appropriate for Fe-bearing systems. Two sets of models were run with different boundary compositions: 1) “low-H models” are constrained by mineral-melt equilibrium partitioning with a basaltic melt that has 0.75 wt% H2O and 2) “high-H models,” which utilize the upper end of the estimated range of mantle water solubility for each phase. Both sets of models yield re-equilibration timescales that are identical and are fast for all phases at a given temperature. These timescales have strong log-linear trends as a function of temperature (R2 from 0.97 to 0.99) that can be used to calculate expected re-equilibration time at a given temperature and grain size. At the high end of the model temperatures (1000-1200°C), H+ completely re-equilibrates in olivine, orthopyroxene, and clinopyroxene within minutes to hours, consistent with previous studies. These short timescales indicate that xenolith NAM mantle water contents are likely to be overprinted prior to eruption. The models also resolve the decoupled water-trace element relationship in Southwest Indian Ridge peridotites, in which peridotite REE abundances are reproduced by partial melting models whereas the relatively high NAM H2O contents require later re-equilibration with melt. At temperatures of 600-800°C, which correspond to conditions of hydrothermal alteration of pyroxene to amphibole and talc, H+ re-equilibration typically occurs over a range of timescales spanning days to years. These durations are well within existing estimates for the duration of fluid flow in oceanic hydrothermal systems, suggesting that peridotite NAM water contents are susceptible to diffusive overprinting during higher temperature hydrothermal alteration. Thus, diffusion during aqueous fluid-rock interactions may also explain NAM H2O contents that are too high to reflect residues of melting. These relatively short timescales at low temperatures suggest that the origin of water contents measured in peridotite NAMs requires additional constraints on sample petrogenesis, including petrographic and trace element analyses. Our 3D model results also hint that H+ may diffuse appreciably during peridotite serpentinization, but diffusion coefficients at low temperature are unconstrained and additional experimental investigations are needed.
more »
« less
This content will become publicly available on November 1, 2025
Distribution of REE between amphibole and pyroxenes in the lithospheric mantle: An assessment from the lattice strain model
Abstract Amphibole and pyroxenes are the main reservoirs of rare earth elements (REEs) in the lithospheric mantle that has been affected by hydrous metasomatism. In this study, we developed semi-empirical models for REE partitioning between orthopyroxene and amphibole and between clinopyroxene and amphibole. These models were formulated on the basis of parameterized lattice strain models of mineral-melt REE partitioning for orthopyroxene, clinopyroxene, and amphibole, and they were calibrated using major element and REE data of amphibole and pyroxenes in natural mantle samples from intraplate settings. The mineral-melt REE partitioning models suggest that amphibole is not in equilibrium with coexisting pyroxenes in the mantle samples and that the amphibole crystallized at a lower temperature than that of the pyroxenes. We estimated the apparent amphibole crystallization temperature using major element compositions of the amphibole and established temperature- and composition-dependent models that can be used to predict apparent pyroxene-amphibole REE partition coefficients for amphibole-bearing peridotite and pyroxenite from intraplate lithospheric mantle. Apparent pyroxene-amphibole REE partition coefficients predicted by the models can be used to infer REE contents of amphibole from REE contents of coexisting pyroxenes. This is especially useful when the grain size of amphibole is too small for trace element analysis.
more »
« less
- Award ID(s):
- 2147598
- PAR ID:
- 10574493
- Publisher / Repository:
- Mineralogical Society of America
- Date Published:
- Journal Name:
- American Mineralogist
- Volume:
- 109
- Issue:
- 11
- ISSN:
- 0003-004X
- Page Range / eLocation ID:
- 1921 to 1933
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Amphibole is a common hydrous mineral in mantle rocks. To better understand processes leading to the formation of amphibole‐bearing peridotites and pyroxenites in the lithospheric mantle, we conducted experiments by juxtaposing a lherzolite against hydrous basaltic melts in Au‐Pd capsules. Two melts were examined, a basaltic andesite and a basalt, each containing 4 wt% of water. The experiments were run at 1200°C and 1 GPa for 3 or 12 h, and then cooled to 880°C and 0.8 GPa over 49 h. The reaction at 1200°C produced a melt‐bearing orthopyroxenite‐dunite sequence. Crystallization of the partially reacted melts during cooling lead to the formation of an amphibole‐bearing gabbronorite‐orthopyroxenite‐peridotite sequence. Orthopyroxene in the peridotite and orthopyroxenite has a poikilitic texture enclosing olivines and spinels. Amphibole in the peridotite occurs interstitial to olivine, orthopyroxene, clinopyroxene, and spinel. Comparisons of texture and mineral compositions in the experimental products with those from field observations allow a better understanding of hydrous melt‐rock reaction in the lithospheric mantle. Amphibole‐bearing pyroxenite veins (or dikes) can be formed in the lithospheric mantle or at the crust‐mantle boundary by interaction between hydrous melt and peridotite and subsequent crystallization. Hornblendite or amphibole gabbronorite can be formed in the veins when the flux of hydrous melt is high. Differences in reacting melt and peridotite compositions are responsible for the variation in amphibole composition in mantle xenoliths from different tectonic settings. The extent of melt‐rock reaction is a factor that control amphibole composition across the amphibole‐bearing vein and the host peridotite.more » « less
-
Abstract Pyroxenite veins and dikes are commonly observed in the mantle section of ophiolites. Because of their mantle occurrence, these pyroxenites are free from crustal contamination and offer a unique opportunity for studying mantle compositions and melt–rock interaction processes. We conducted an integrated petrological and geochemical study of a suite of composite orthopyroxenite, websterite, and pyroxene-bearing dunite veins from the Xiugugabu ophiolite located on the western segment of Yarlung–Zangbo Suture Zone. The dunite is separated from the host peridotite by a layer of pyroxenite, forming a composite vein system. Systematic variations in major, minor, and trace element compositions in minerals across the composite veins are observed. Two generations of orthopyroxenes in the pyroxenites are characterized by high Mg#, low TiO2 concentrations, and depleted patterns of incompatible trace elements. Clinopyroxenes in the pyroxenites are characterized by high Mg#, low contents of TiO2 and Na2O, spooned shaped REE patterns, and a negative Zr anomaly. Through major and trace element modeling, we showed that both orthopyroxene and clinopyroxene were in equilibrium with melts with different compositions. This hypothesis is further confirmed by distinct initial Nd and Hf isotope ratios in the two pyroxenes. A model for the formation of composite pyroxenite veins is developed, whereby hydrous and silica-rich melts percolate along the margins of a dunite channel. The orthopyroxenite was formed by the reaction between a hydrous, silica-rich melt and the surrounding peridotite. The websterite is formed by reactive crystallization of a hybrid melt produced by mixing silica-rich melt and the melt formed by remelting of previously depleted peridotite in the deeper part of the mantle column. The extremely enriched Nd–Hf isotope compositions of the pyroxenite veins (εNd = −20.3 to +11.5 and εHf = −13.2 to +25.3, 125 million years ago) can be explained by the addition of ancient, recycled sediments to the mantle source in a supra-subduction setting. Based on the low-Cr# spinel in the Xiugugabu dunites (Cr# = 19–50) and the depleted nature of the parental melt of the Xiugugabu pyroxenites, we deduced that the formation of pyroxenites postdate the formation of the Xiugugabu ophiolite at ~125–130 Ma. Collectively, results from this study have provided support to the hypothesis that the Xiugugabu ophiolite experience a two-stage evolution, i.e., firstly formed in a mid-ocean ridge setting and subsequently modified in a supra subduction zone.more » « less
-
null (Ed.)Mantle heterogeneity has a first-order control on the petrological and geochemical differences of erupted mafic lavas across the globe. It is debated whether this heterogeneity reflects only chemical variability or also lithological differences in source regions. Because of their various partitioning behaviors between mantle minerals, First Row Transition Elements (FRTEs) have been identified as potential lithological tracers. Here, we investigate the various parameters that control FRTE partitioning between common mantle phases through a comparison of partition coefficients calculated from natural pyroxenites obtained from the Earthchem database with previous partitioning experiments and new electron microprobe analyses. Using naturally occurring pyroxenites from alpine massifs and xenoliths provides the opportunity to explore the behavior of FRTEs on a much larger range of compositions and temperatures than covered by experimental studies. Our preliminary results show that natural partition coefficients for Fe and Mn depend on temperature and vary distinctly between lithologies. The effect of composition, however, is difficult to resolve and will require further inspection. Natural exchange coefficients, or Kd’s (mineral/mineral) for Mn/Fe, largely match previous experimental data across peridotite and pyroxenite compositions for garnet/clinopyroxene(cpx), orthopyroxene/cpx, and olivine/cpx. However, natural samples often present evidence of chemical disequilibrium and/or secondary alteration which can significantly increase the scatter in analyses. Importantly, despite the larger uncertainty on the natural Kd’s than on experimental ones, natural exchange coefficients show distinct values between the various pairs of minerals. These distinctions, and the fact that Kd’s do not seem to be influenced by temperature, make the bulk Mn/Fe ratio in lavas a good lithological tracer, supporting previous claims. Hence, we show that natural compositions can be used to expand trends in FRTE distribution behavior across a wider range of temperatures (500-1500°C) and compositions than determined previously by experiments alone.more » « less
-
Abstract A major challenge in mantle geochemistry is determining the source composition and melt fraction involved in melting. We provide a new Rare‐Earth Element (REE) inverse model that provides source concentration, source and melt mineral modes, and melt fraction based on the difference between separate determinations of bulk distribution coefficients and constrained by boundary conditions. An analytical inverse of the batch melting equation provides expressions for source, , and bulk distribution coefficient of the mantle, , with two unknowns, the initial concentration of La in the mantle, , andPi, the bulk distribution coefficient of the melt. We traverse through a range of steps and examine thousands of melt modes,Pi, at each step. Thousands of trial melt modes fail by generating that are inconsistent with partition coefficients. Many surviving trials cannot be inverted to estimate a mantle mode. Other boundary conditions eliminate even more trials. Surviving trials are ordered by the difference between calculated from the REE data of a lava suite and calculated from partition coefficients and mantle mode. We select the solution with the closest fit that passes all the boundary conditions. We tested our new model with lava suites from Hawaii where different lines of evidence suggest that they melted from different mantle sources, Mauna Kea representing shield‐stage lava and submarine Kiekie representing rejuvenated stage lava. Our inverse determination of mantle composition and melting parameters was consistent with earlier models based on assumptions of HREE composition.more » « less
