Hydrogen, helium, oxygen, iron, silicon, and magnesium are the main constituent elements for planets. Depending on the ratio of these elements, a wide compositional range of planets can be built, from rocky planets to ice and gas giants. Although the internal structures of the rocky planets in the inner solar system have been relatively well studied through modeling and laboratory studies8-10, those of the outer-solar system planets (i.e. Uranus and Neptune) rich in volatile elements remain under-explored. Yet, understanding the geophysics and geochemistry of water-rich planets are increasingly more important as sub-Neptunes, many of which are likely rich in H2O, appear to be the most common type of planets in our galaxy1.
In conventional models, water-rich planets (water-worlds) have been assumed to have separate layers of atmosphere, ice/fluid, rocky mantle, and metallic core, positioning the layer interface between H2O and the rocky layers at high pressure and temperature11. Although hydrous minerals have been widely studied at subducting slab conditions under H2O-controlled (or limited) environments12, chemical interaction and thermodynamic processes of major rock-forming minerals at the H2O-rock interface conditions of water-rich planetary interiors have been scarcely explored particularly at high pressures. Based on the balance between the equilibrium temperature and effective temperature6, the existence of heavy elements in the ice/fluid layer of Uranus and Neptune has been proposed13-14, challenging the conventional view of the well differentiated planetary body and a sharp H2O-rock interface. To provide further insights into the interior models of these water-rich planets and exoplanets, it is essential to understand the phase relations in the H2O-MgO-FeO-SiO2 system at high pressure and temperature (P–T) conditions.
Although previous studies on MgO, SiO2, H2O, and combined MgO-SiO2-H2O system12 have led to the discovery of various hydrous silicates and hydroxide phases at high pressure and temperatures, such as phase D and phase H15-16, experimental conditions have been focused on temperatures expected in the subducting slabs of the Earth17-18. In order to understand the interior of water-rich planets, we have investigated typical rock-forming minerals, olivine ((Mg0.9,Fe0.1)2SiO4) and ferropericlase ((Mg0.9,Fe0.1)O), under H2O-rich conditions. We found that Mg, the major rock-forming cation, becomes highly soluble in dense H2O environments at high P–T, requiring a revision of interior models such as the geochemical cycle and thermal evolution of water-rich planets.
Under a H2O pressure medium, the diffraction peaks from San Carlos olivine disappear within 5–10 seconds upon heating to 1770 K at 33 GPa, and the diffraction pattern becomes dominated by peaks from stishovite (SiO2) with weaker peaks from brucite (Mg(OH)2) and ε-FeOOH (Fig. 1A). This result is in contrast to observations under anhydrous conditions at the same P–T conditions: bridgmanite and ferropericlase form from olivine. Although we observed transient peaks from other Mg-bearing phases (such as phase D and akimotoite, see SI text #1) during heating, they mostly disappeared near the end of heating and temperature quenching at high pressures. Compared to the expected amounts of Mg-rich phases based on the mass balance assuming the starting material and stishovite as a reaction product, the peak intensity of brucite was much weaker than expected. After temperature quenching at high pressures, the diffraction lines of brucite gained intensity, but not to the extent expected considering those of stishovite. The intensities of brucite peaks, however, increased dramatically after quenching the pressure to 1 bar (Fig. 1A). Rietveld refinements of the diffraction pattern yield a 1:1.7 molar ratio between stishovite and brucite (Extended Data Fig. 1), which is in reasonable agreement with the expected molar ratio of 1:1.8 from the breakdown reaction of the San Carlos olivine ((Mg0.9,Fe0.1)2SiO4 + 2H2O → 1.8Mg(OH)2 + SiO2 + 0.2FeOOH), assuming all iron is consumed to produce ε-FeOOH19.
At high temperatures in the 27–50 GPa range, Mg-rich phases were either systematically absent or presented in minor amounts. Although brucite is observed with weak intensities in the diffraction patterns during some heating runs, it is unlikely stable at these conditions because our temperatures were 200–900 K above the known dehydration temperatures at the same pressures20-21. We interpret the observation of brucite as originating from the low-temperature areas between liquid-H2O and diamond anvils (second image in Fig. 2E). Periclase (MgO), as a dehydration product of brucite, was not observed (Extended Data Table 1), supporting to the conclusion that the missing MgO content is no longer crystalline. Our preferred interpretation is that the MgO component exists in a dissolved form in liquid H2O during laser heating at high pressures. The crystallization of a large amount of brucite after quenching to 1 bar and 300 K (also confirmed by an infrared (IR) spectroscopy study, see SI text #2) can thus be explained by the reduced solubility of MgO component in H2O with decreases in pressure and temperature and subsequent precipitation of the dissolved MgO-component into Mg(OH)2.
To gain further insights into the processes occurred during the laser heating at high pressures, we analyzed the laser-heated spots of the recovered samples using scanning electron microscopy (SEM) combined with the focused (Ga+) ion beam (FIB) and energy dispersive X-ray spectroscopy (EDS). Above the laser-heated spots, new dome-like structures were found (Fig. 2A). The vertical cross sections of the structure reveal that the dome is highly porous with an empty space inside (Figs. 2B and 2C). The domes are rich in Mg whereas below the domes, in place of the original starting material, is concentrated in Si (Fig. 2D). The minerals of the dome structures show cleavage planes, which are characteristic of the morphology of a layered crystal structure, such as brucite (white box in Fig. 2C). Based on the chemical compositions combined with the sample morphology and the synchrotron XRD observations, we propose that the domes are composed of brucite precipitates originating from MgO dissolved in H2O from the starting olivine sample. The Si released during decomposition resulted consequently in the formation stishovite during high-pressure laser heating.
Because H2O medium does not couple directly with near IR laser beams, the heat from the sample, which does couple, does provide sufficient heat to melt the surrounding H2O-ice (Fig. 2E). At high temperatures, pockets of molten H2O would form and contain dissolved MgO originating from olivine. Upon quenching from high temperature at high pressures, only a marginal increase in the peak intensities of brucite would suggest that high-pressure ice could still store most of the dissolved MgO (Fig. 1A). Finally, after decompression to almost 1 bar, most of the dissolved MgO precipitates as brucite because of its low solubility in water at pressures. Such processes from high P–T to 1 bar, and 300 K must have resulted in the formation of porous dome structures as observed in the electron microscopy (Fig. 2) and an abrupt increase in the brucite intensity as observed in the diffraction patterns (Fig. 1A) and IR spectra (Extended Data Fig. 3, see SI text #3). In conclusion, the following chemical reaction can be constructed for the decomposition of olivine when subjected to high P–T conditions under H2O saturated conditions:
Mg2SiO4 + 2H2O → 2Mg2+ (aq) + 4OH– (aq) + SiO2 (s).
When ferropericlase ((Mg0.9,Fe0.1)O) starting material was heated to 1660 K at 36 GPa under a H2O medium, peaks from brucite and ε-FeOOH, though weak, appeared while ferropericlase peaks disappeared (Fig. 1B, see SI text #4). Similar to the results from the olivine experiment, the intensity of brucite peaks increased in steps after quenching to 300 K and decompression to 1 bar. The SEM images of the recovered samples also showed highly porous features in the heated spot (Extended Data Fig. 4). Unlike the olivine + water experiments, however, the dome-like structures (Figs. 2A–2C) were not observed from the ferropericlase + water experiment. Ferropericlase does not contain SiO2, hence undissolving materials (e.g. stishovite) into H2O were not left at the center of heating spot. In other words, only the leaching process would have proceeded at the center while brucite recrystallized at low-temperature areas. Therefore, we observed the porous texture at the center instead of dome structure.
At pressures above 60 GPa, bridgmanite and ferropericlase became the dominant phases in the olivine + H2O experiment (Extended Data Fig. 5). Our diffraction patterns measured at high P–T conditions and upon temperature quench do not show any clear evidence for brucite, suggesting the decrease in the solubility of MgO in H2O above 60 GPa (Extended Data Table 1). We also found a similar trend for the ferropericlase + water experiment at pressures above 60 GPa: brucite peaks were not observed during laser heating and in the quenched sample at high pressures (Extended Data Fig. 6, see SI text #5). The SEM images of the recovered samples from these experiments still showed porous textures in the heated spots (Extended Data Fig. 4), indicating certain level of solubility of MgO in H2O persists at high pressures above 60 GPa albeit smaller.
An increase in the content of MgO in the fluid has been reported from the recovered sample in a similar Mg-silicate system from a multi-anvil press experiment at pressures below 10.5 GPa22, although performed at lower temperatures (up to 1500 K) than our laser-heated DAC experiments. MgO in the fluid (wt%) increased as a function of temperature and pressure, suggesting pressure-promoted solubility of MgO increase in H2O. While there is a gap in pressure between the multi-anvil press study and our DAC experiment, the general trends are in agreement supporting the enhanced solubility of MgO at high pressures and temperatures.
The amount of dissolved MgO in H2O at high pressures is difficult to estimate from the XRD patterns alone. In the SEM images (Figs. 2B-2D), the volume of the heated spot appears to expand by approximately 3-4 times from the initial sample thickness. Based on the solid residue and the pore volume, we estimate the volume ratio between MgO and H2O involved in the reaction to be between 1:5 and 1:10. Considering molar volumes at the reaction conditions at 30 GPa, the solubility of MgO in H2O could be then between 200 and 400 g/L (see SI text #6). This amount is comparable to the solubility of NaCl in pure H2O at 1 bar, i.e., 360 g/L.
In an attempt to estimate variation in MgO solubility with pressure, we measured the intensities of selected ferropericlase peaks as a function of pressure. The degree of the decrease in the Bragg peak intensity would be proportional to the amount of MgO loss from the crystalline phase, and hence the amount of MgO dissolved in H2O. At 24–38 GPa and 1400–1850 K, ferropericlase peak intensity decreases by up to 87–100% upon heating (Fig. 3A), meaning most of MgO is in dissolved state. The intensity only decreased to 10–40% of the original peak as pressure increases to 80 GPa. However, at pressures above 55 GPa, our temperatures are below the H2O melting when considering data reported by Schwager and Boehler (2008)23 (Fig. 3B). Therefore, it is feasible that the decreased solubility of MgO above 55 GPa is due to a solid-solid reaction. On the other hand, the melting temperatures of H2O at such high pressures are highly controversial, and if the melting curve of Lin et al. (2005)24 is used, our temperatures estimates lay above the melting.
While water-rich planets are common in our galaxy (including Uranus and Neptune in our solar system), their internal structures and geochemical cycles are not well understood13. Our new data here provide important insights into constraining the internal structures of the water-rich planets. While a wide range of sizes exists in the sub-Neptune class1, we focus here on two sizes where the water-rock interfaces are within our experimental pressure range: (1) 0.7-1.4MEarth and 1.0-1.1REarth planets (such as TRAPPIST-1c and -1f)7 and (2) 6.6MEarth and 2.7REarth planets (such as GJ 1214b)25. The water-rock interface of these planets would exist at approximately 25 GPa26 and 40-140 GPa27, respectively, assuming that the planets are composed of H2O and rocky layers with a metallic iron core. We then consider two different thermal gradients within the planet: (1) sufficiently high temperatures above the melting curve of H2O, and (2) sufficiently low temperatures below the melting curve of H2O. The former can be considered as young water-rich planets and the latter represents an older version (Fig. 4). For the sake of simplicity, we only assume pure H2O for the outer layer and Earth’s mantle like composition for the rock layer, although the thermally insulating atmosphere could also play a role in containing the internal heat of the water-rich planets2, 28.
For Earth-size planets (such as TRAPPIST-1c and -1f), the pressures of the water-rock interface would be close to where we identify the peak solubility of MgO in liquid H2O at high temperatures. If such a planet is sufficiently warm to stay above the H2O melting temperature, a large amount of MgO could be dissolved into the liquid H2O layer. Furthermore, if the temperature of the planet is high enough for the hydrated rock layer to be molten29, causing a vigorous convection in the magma ocean30, the topmost rock layer would be continuously rejuvenated, supplying a large amount of MgO into the H2O layer. In exchange, H2O would also be dissolved into the magma ocean because H2O solubility in magma increases dramatically with pressure31-33. Finally, if H2O and silicate magma are fully miscible at sufficient high pressures34, the boundary between the rock and H2O layers would be fuzzy. On the other hand, if the rock layer below the liquid H2O layer is solidified, MgO would still remain dissolved in the H2O layer and therefore the solid near the H2O-rock interface would be enriched in SiO2, potentially resulting in a silica-rich layer above the remaining rocky mantle that will have an effect on the structure of the planet (Fig. 4B, see SI text #7 for details).
Because MgO solubility in H2O-ice would be somewhat lower than that in H2O fluid as observed in our experiment, some amounts of Mg-silicates and/or Mg-bearing hydrous minerals would precipitate upon the solidification of H2O and deposit at the top of the silica-rich (hydrated) layer above the rocky mantle (yellow dots in Fig. 4B). The density of hydrous minerals, such as brucite and phase D, are much lower by 20-40% compared with silica and bridgmanite. Therefore, it would be difficult for the Mg-rich hydrous mineral layer to be mixed with the rest of the rock layer, especially in cooler planets with reduced convection in the rock layer.
Now let us consider the larger water-rich planets (such as GJ 1214b: 2.68REarth) that are warm enough to dissolve MgO into H2O layer. These planets would also have an MgO-dissolved H2O layer based on our higher pressure experiments. If MgO and H2O liquids are miscible at high pressure34, much MgO would exist in the H2O layer, while the magma ocean below would be hydrated due to the earlier described exchange through the interface. Thus, the boundary would be undefined over a large depth range. Once the rock layer is solidified by cooling, MgO content at greater depths in the H2O layer would be lowered as we observed reduced reaction between olivine and H2O above 60 GPa. Consequently, the thickness of the silica-rich layer at the top of the rock layer could be smaller in comparison to the smaller water rich-planets. In addition, the topmost silica layer would not be as MgO-deficient as in the smaller Earth-sized water-rich planet as the reaction is less extensive at higher pressure conditions (Fig. 4B).
Because the density of silica at 100 GPa is higher than bridgmanite and pyrolite, similar to the case of smaller Earth-sized water-rich planets, the lower temperature expected for the silica layer would further increase the density contrast at the interface. However, silica hydration is likely to be promoted by pressure35. For example, 13 wt% H2O can reduce the density of silica by 10%. Therefore, the silica layer is less likely to be mixed with the rock layer below, potentially acting as a compositional barrier between the water and rock layers in the large water-rich planets. Similar to the smaller Earth-size water-rich planets, hydrous minerals would precipitate at the top of the rock layer during the cooling of a large water-rich planet while high-pressure ice can also store certain amount of MgO. It is therefore feasible that an Mg-rich hydrated layer would form at the top of the rocky mantle (Fig. 4B).
Our data can also help to understand Uranus, a water-rich planet in our solar system that is cold (its atmospheric temperature at 49 K)5 and shows very low heat flux6 of 0.042x10-4 W/cm2. Its low luminosity may also suggest rapid cooling of the outer envelope while much heat might be trapped at greater depths. It has been hypothesized that some form of a compositional barrier exists in Uranus's upper layer preventing the heat transfer from the deeper interior3-4. Our experimental results support this idea because we found that the solubility of MgO in H2O may peak at 20-40 GPa especially considering that a large amount of rocky materials may have been mixed with water to shallow depths during large impacts in early Uranus36. Thus, during cooling, the rock components would settle to the greater depths because of their high density. However, at the depth range where the solubility is high, MgO may remain in water layer, resulting in a chemically stabilized compositional barrier as has been proposed.
Our experiments also predict compositional layering at the interface between the H2O and rock layers in Uranus, i.e., Mg-rich hydrated layer and Si-rich hydrated layer (Fig. 4B). If hydrated as supported by recent experiments35, the layer would be difficult to mix with the rock layer below and therefore likely form a dynamically stable compositional gradient at the interface between H2O and rock layers. However, if sufficient heat is trapped in Uranus’ deep interior due to the thermal barrier, supported by our results, the rocky portions could be molten2, 36, facilitating vigorous mixing instead of layering at the interface. Therefore, the deeper compositional gradients of larger water-rich planets remain open question.