skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The Viscosity and Atomic Structure of Volatile-Bearing Melilititic Melts at High Pressure and Temperature and the Transport of Deep Carbon
Understanding the viscosity of mantle-derived magmas is needed to model their migration mechanisms and ascent rate from the source rock to the surface. High pressure–temperature experimental data are now available on the viscosity of synthetic melts, pure carbonatitic to carbonate–silicate compositions, anhydrous basalts, dacites and rhyolites. However, the viscosity of volatile-bearing melilititic melts, among the most plausible carriers of deep carbon, has not been investigated. In this study, we experimentally determined the viscosity of synthetic liquids with ~31 and ~39 wt% SiO2, 1.60 and 1.42 wt% CO2 and 5.7 and 1 wt% H2O, respectively, at pressures from 1 to 4.7 GPa and temperatures between 1265 and 1755 °C, using the falling-sphere technique combined with in situ X-ray radiography. Our results show viscosities between 0.1044 and 2.1221 Pa·s, with a clear dependence on temperature and SiO2 content. The atomic structure of both melt compositions was also determined at high pressure and temperature, using in situ multi-angle energy-dispersive X-ray diffraction supported by ex situ microFTIR and microRaman spectroscopic measurements. Our results yield evidence that the T–T and T–O (T = Si,Al) interatomic distances of ultrabasic melts are higher than those for basaltic melts known from similar recent studies. Based on our experimental data, melilititic melts are expected to migrate at a rate ~from 2 to 57 km·yr−1 in the present-day or the Archaean mantle, respectively.  more » « less
Award ID(s):
1761388 1732256
PAR ID:
10140589
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Minerals
Volume:
10
Issue:
3
ISSN:
2075-163X
Page Range / eLocation ID:
267
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Many lines of evidence from high P–T experiments, thermodynamic models, and natural observations suggest that slab-derived aqueous fluids, which flux mantle wedges contain variable amounts of dissolved carbon. However, constraints on the effects of H2O–CO2 fluids on mantle melting, particularly at mantle wedge P–T conditions, are limited. Here, we present new piston cylinder experiments on fertile and depleted peridotite compositions with 3.5 wt.% H2O and XCO2 [= molar CO2 / (CO2 + H2O)] of 0.04–0.17. Experiments were performed at 2–3 GPa and 1350°C to assess how temperature, peridotite fertility, and XCO2 of slab-derived fluid affects partial melting in mantle wedges. All experiments produce olivine + orthopyroxene +7 to 41 wt.% partial melt. Our new data, along with previous lower temperature data, show that as mantle wedge temperature increases, primary melts become richer in SiO2, FeO*, and MgO and poorer CaO, Al2O3, and alkalis when influenced by H2O–CO2 fluids. At constant P–T and bulk H2O content, the extent of melting in the mantle wedge is largely controlled by peridotite fertility and XCO2 of slab-fluid. High XCO2 depleted compositions generate ~7 wt.% melt, whereas, at identical P–T, low XCO2 fertile compositions generate ~30 to 40 wt.% melt. Additionally, peridotite fertility and XCO2 have significant effects on peridotite partial melt compositions. At a constant P–T–XCO2, fertile peridotites generate melts richer in CaO and Al2O3 and poorer in SiO2, MgO + FeO, and alkalis. Similar to previous experimental studies, at a constant P–T fertility condition, as XCO2 increases, SiO2 and CaO of melts systematically decrease and increase, respectively. Such distinctive effects of oxidized form of dissolved carbon on peridotite partial melt compositions are not observed if the carbon-bearing fluid is reduced, such as CH4-bearing. Considering the large effect of XCO2 on melt SiO2 and CaO concentrations and the relatively oxidized nature of arc magmas, we compare the SiO2/CaO of our experimental melts and melts from previous peridotite + H2O ± CO2 studies to the SiO2/CaO systematics of primitive arc basalts and ultra-calcic, silica-undersaturated arc melt inclusions. From this comparison, we demonstrate that across most P–T–fertility conditions predicted for mantle wedges, partial melts from bulk compositions with XCO2 ≥ 0.11 have lower SiO2/CaO than all primitive arc melts found globally, even when correcting for olivine fractionation, whereas partial melts from bulk compositions with XCO2 = 0.04 overlap the lower end of the SiO2/CaO field defined by natural data. These results suggest that the upper XCO2 limit of slab-fluids influencing primary arc magma formation is 0.04 < XCO2 < 0.11, and this upper limit is likely to apply globally. Lastly, we show that the anomalous SiO2/CaO and CaO/Al2O3 signatures observed in ultra-calcic arc melt inclusions can be reproduced by partial melting of either CO2-bearing hydrous fertile and depleted peridotites with 0 < XCO2 < 0.11 at 2–3 GPa, or from nominally CO2-free hydrous fertile peridotites at P > 3 GPa. 
    more » « less
  2. Both seismic observations of dense low shear velocity regions and models of magma ocean crystallization and mantle dynamics support enrichment of iron in Earth’s lowermost mantle. Physical properties of iron-rich lower mantle heterogeneities in the modern Earth depend on distribution of iron between coexisting lower mantle phases (Mg,Fe)O magnesiowüstite, (Mg,Fe)SiO3 bridgmanite, and (Mg,Fe)SiO3 post-perovskite. The partitioning of iron between these phases was investigated in synthetic ferrous-iron-rich olivine compositions (Mg0.55Fe0.45)2SiO4 and (Mg0.28Fe0.72)2SiO4 at lower mantle conditions ranging from 33–128 GPa and 1900–3000 K in the laser-heated diamond anvil cell. The resulting phase assemblages were characterized by a combination of in situ X-ray diffraction and ex situ transmission electron microscopy. The exchange coefficient between bridgmanite and magnesiowüstite decreases with pressure and bulk Fe# and increases with temperature. Thermodynamic modeling determines that incorporation and partitioning of iron in bridgmanite are explained well by excess volume associated with Mg-Fe exchange. Partitioning results are used to model compositions and densities of mantle phase assemblages as a function of pressure, FeO-content and SiO2-content. Unlike average mantle compositions, iron-rich compositions in the mantle exhibit negative dependence of density on SiO2-content at all mantle depths, an important finding for interpretation of deep lower mantle structures. 
    more » « less
  3. Xu, Hongwu (Ed.)
    Abstract We have measured the sound velocities and elasticity of synthetic polycrystalline β-Mg2SiO4 containing 1.2 wt% H2O to 10 GPa and 600 K using ultrasonic interferometry with synchrotron X-radiation. We determined sample length at high pressure and temperature using the sample’s X-radiographic image and applied travel times bond corrections appropriate to the experimental cell assembly configuration. Fitting the entire moduli data to third-order finite strain equations yields the adiabatic bulk [KS0 = 153.3(4) GPa] and shear [G0 = 101.8(2) GPa] moduli, their pressure derivatives (∂KS/∂P)T = 5.15(6) and (∂G/∂P)T = 1.68(3) and temperature derivatives (∂KS/∂T)P = −0.0179(9) GPa/K and (∂G/∂T)P = −0.0151(7) GPa/K. Comparing the bulk sound velocity contrast between the new hydrous wadsleyite data and olivine (0.38 wt% H2O) with seismic bulk sound velocity contrasts of 3.5% and 4.0% yields 53% and 60% olivine content, respectively, assuming an iso-chemical mantle model of the Earth. The results suggest that a hydrous mantle transition zone with a pyrolite model composition could explain the 410 km seismic velocity jump. 
    more » « less
  4. Abstract A collection of quaternary, high-MgO (≤13.4 wt%) basanite and minette cinder and lava cones, with an enhanced arc geochemical signature, are located along the northern margin of the N–S Colima rift in western Mexico. The Colima rift overlies the lithospheric suture between the Jalisco block and Guerrero terrane, as well as the tear between the Rivera and Cocos subducting oceanic plates. From the literature, volatile analyses of olivine-hosted melt inclusions in the Colima cone samples document notably high concentrations of dissolved H2O in the melt (≤ 7.0 wt%) as well as degassing-induced phenocryst growth over a range of depths ≤25 km. In this study, it is shown that the high-MgO character of the Colima suite reflects liquid compositions, consistent with evidence for their rapid transit to the surface, without stalling in a crustal magma chamber. The most Mg-rich olivine analyzed in each sample matches the equilibrium composition at the liquidus based on olivine-melt Mn–Mg and Fe2+–Mg exchange coefficients. Application of both a Mg- and Ni-based olivine-melt thermometer, calibrated on the same experimental data set, to the Colima cone suite provides the temperature and dissolved H2O content at the liquidus. Because the Ni thermometer is insensitive to dissolved H2O in the melt, it gives the actual temperature at the onset of olivine phenocryst growth. For the nine Colima samples that range from 13.4–9.2 wt% MgO, resulting temperatures range from 1221°C to 1056°C (± 6–11°C). In contrast, the Mg thermometer is sensitive to dissolved H2O in the melt, and its application (without a correction of H2O) gives the temperature of olivine crystallization under anhydrous conditions. When the Mg- and Ni-based temperatures are paired, the depression of the liquidus (∆T = TMg–TNi) due to dissolved H2O can be obtained. For the high-MgO (>9 wt%) Colima samples, ∆T values range from 188°C to 109°C. Corrections for the effect of pressure (i.e. evidence that phenocryst growth began at ~700 MPa), increase ∆T by ~21°C. An updated model calibration (on experiments from the literature) that relates ∆T with dissolved H2O in the melt shows that the best fit (R2 = 0.95) is linear, wt% H2O = 0.047*∆T, with a standard error of ±0.5 wt%. Although the experimental data set spans a wide range of melt composition (e.g. 47–58 wt% SiO2, 4.4–10.2 wt% MgO, 1.3–4.9 wt% Na2O, 0.1–5.0 wt% K2O, 0.3–5.3 wt% H2O), no dependence on anhydrous melt composition is resolved. Application of this updated model to the Colima suite gives H2O contents of 5.1–8.8 wt% H2O, consistent with those analyzed in olivine-hosted MIs from the literature. When the thermometry and hygrometry results for the Colima cone suite are compared to those for the adjacent calc–alkaline basalts from the Tancítaro Volcanic Field (TVF) in Michoacán, two distinct linear trends in a plot of wt% H2O vs. temperature are found, indicative of different mantle sources. It is proposed that the high-MgO (>11 wt%) Colima cone melts were derived from a phlogopite-bearing harzburgitic mantle at the base of the Jalisco block lithosphere, whereas both TVF and Colima melts with ≤10 wt% MgO were derived from the asthenosphere (i.e. arc mantle wedge). In both mantle sources, slab-derived fluids were an important source of H2O. 
    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