An official website of the United States government
Abstract A comprehensive analysis of experimental data and theoretical simulations on the partial molar volume of water in silicate melt indicates that finite strain theory successfully describes the compression of the H2O component dissolved in silicate melt at high pressures and temperatures. However, because of the high compressibility of the water component, a fourth order equation of state fit is required to accurately simulate experimental results on water's volume in silicate melts at a deep upper mantle, transition zone, and lower mantle pressures. Data from previous shock compression experiments on hydrous minerals in which melting occurs along the Hugoniot are used to provide an experimental constraint on the partial molar volume of water in silicate melt at deep mantle temperatures and pressures. The equation of state of the water component indicates that, depending on elastic averaging technique, the amount of water that could be present in neutrally or negatively buoyant mafic/ultramafic melts above the 410 km seismic discontinuity is upper‐bounded at 5.6 wt%: smaller than previously inferred, and consistent with melt being confined to a narrow depth range above the 410 km discontinuity. If melt is predominantly distributed along grain boundaries in low aspect ratio films, extents of melting as low as 2% could produce observed seismic velocity reductions. The ability of the lowermost mantle to contain negatively buoyant hydrous liquids hinges on the trade‐off between iron content and hydration: at these depths, substantially higher degrees of hydration could be present within partial melts.
more »
« less
This content will become publicly available on December 1, 2026
Low melt viscosity enables melt doublets above the 410-km discontinuity
Abstract Seismic and magnetotelluric studies suggest hydrous silicate melts atop the 410 km discontinuity form 30–100 km thick layers. Importantly, in some regions, two layers are observed. These stagnant layers are related to their comparable density to the surrounding mantle, but their formation mechanisms and detailed structures remain unclear. Here we report a large decrease of silicate melt viscosity at ~14 GPa, from 96(5) to 11.7(6) mPa⋅s, as water content increases from 15.5 to 31.8 mol% H₂O. Such low viscosities facilitate rapid segregation of melt, which would typically prevent thick layer accumulation. Our 1D finite element simulations show that continuous dehydration melting of upwelling mantle material produces a primary melt layer above 410 km and a secondary layer at the depth of equal mantle-melt densities. These layers can merge into a single thick layer under low density contrasts or high upwelling rates, explaining both melt doublets and thick single layers.
more »
« less
- Award ID(s):
- 2022492
- PAR ID:
- 10597570
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 16
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Density of silicate melt dictates melt migration and establishes the gross structure of Earth's interior. However, due to technical challenges, the melt density of relevant compositions is poorly known at deep mantle conditions. Particularly, water may be dissolved in such melts in large amounts and can potentially affect their density at extreme pressure and temperature conditions. Here we perform first‐principles molecular dynamics simulations to evaluate the density of Fe‐rich, eutectic‐like silicate melt (Emelt) with varying water content up to about 12 wt %. Our results show that water mixes nearly ideally with the nonvolatile component in silicate melt and can decrease the melt density significantly. They also suggest that hydrous melts can be gravitationally stable in the lowermost mantle given its likely high iron content, providing a mechanism to explain seismically slow and dense layers near the core‐mantle boundary.more » « less
-
Abstract Silicate melts have served as transport agents in the chemical and thermal evolution of Earth. Molecular dynamics simulations based on a deep neural network potential trained byab initiodata show that the viscosity of MgSiO3melt decreases with increasing pressure at low pressures (up to ∼6 GPa) before it starts to increase with further compression. The melt electrical conductivity also behaves anomalously; first increasing and then decreasing with pressure. The melt accumulation implied by the viscosity turnover at ∼23 GPa along mantle liquidus offers an explanation for the low‐velocity zone at the 660‐km discontinuity. The increase in electrical conductivity up to ∼50 GPa may contribute to the steep rise of Earth's electrical conductivity profiles derived from magnetotelluric observations. Our results also suggest that small fraction of melts could give rise to detectable bulk conductivity in deeper parts of the mantle.more » « less
-
Abstract By taking advantage of the recent availability of a broadband seismic data set from Networks NR and BX covering the entire country of Botswana, we conduct a systematic receiver function investigation of the topography of the 410 and 660 km discontinuities beneath the incipient Okavango rift zone (ORZ) in northern Botswana and its adjacent Archean‐Proterozoic tectonic provinces in southern Africa. Similar to a previous mantle transition zone (MTZ) discontinuity study using data from a 1‐D profile traversing the ORZ, a normal MTZ thickness is observed in most parts of the study area. This is inconsistent with the existence of widespread positive thermal anomalies in the MTZ and further implies that active thermal upwelling from the lower mantle plays an insignificant role in the initiation of continental rifting. The results also suggest that cold temperature presumably associated with thick cratonic keels has indiscernible influence on the thermal structure of the MTZ. The expanded data set reveals several isolated areas of slight (~10 km or smaller) MTZ thinning. The largest of such areas has a NE‐SW elongated shape and is mostly caused by relative deepening of the 410 km discontinuity rather than shallowing of the 660 km discontinuity. These characteristics are different from those expected for a typical mantle plume. We speculate that the thinner‐than‐normal MTZ may be induced by minor thermal upwelling associated with late Mesozoic‐early Cenozoic lithospheric delamination, a recently proposed mechanism that might be responsible for the high elevation of southern Africa.more » « less
-
Abstract Surface meltwater runoff dominates present-day mass loss from the Greenland Ice Sheet. In Greenland’s interior, porous firn can limit runoff by retaining meltwater unless perched low-permeability horizons, such as ice slabs, develop and restrict percolation. Recent observations suggest that such horizons might develop rapidly during extreme melt seasons. Here we present radar sounding evidence that an extensive near surface melt layer formed following the extreme melt season in 2012. This layer was still present in 2017 in regions up to 700 m higher in elevation and 160 km further inland than known ice slabs. We find that melt layer formation is driven by local, short-timescale thermal and hydrologic processes in addition to mean climate state. These melt layers reduce vertical percolation pathways, and, under appropriate firn temperature and surface melt conditions, encourage further ice aggregation at their horizon. Therefore, the frequency of extreme melt seasons relative to the rate at which pore space and cold content regenerates above the most recent melt layer may be a key determinant of the firn’s multi-year response to surface melt.more » « less
