Internal water storage capacity of terrestrial planets and the effect of hydration on the M-R relation

Understanding the chemical interactions between water and Mg-silicates or iron is essential to constrain the interiors of water-rich planets. Hydration effects have, however, been mostly neglected by the astrophysics community so far. As such effects are unlikely to have major impacts on theoretical mass-radius relations this is justified as long as the measurement uncertainties are large. However, upcoming missions, such as the PLATO mission (scheduled launch 2026), are envisaged to reach a precision of up to $\approx 3 \%$ and $\approx 10 \%$ for radii and masses, respectively. As a result, we may soon enter an area in exoplanetary research where various physical and chemical effects such as hydration can no longer be ignored. Our goal is to construct interior models for planets that include reliable prescriptions for hydration of the cores and the mantles. These models can be used to refine previous results for which hydration has been neglected and to guide future characterization of observed exoplanets. We have developed numerical tools to solve for the structure of multi-layered planets with variable boundary conditions and compositions. Here we consider three types of planets: dry interiors, hydrated interiors and dry interiors + surface ocean where the ocean mass fraction corresponds to the mass fraction of $\rm H_2 O$ equivalent in the hydrated case. We find H/OH storage capacities in the hydrated planets equivalent to $0-6 \rm \ wt\% \ \rm H_{2}O$ corresponding to up to $\approx 800 \rm \ km$ deep ocean layers. In the mass range $0.1 \leq M/M_\oplus \leq 3$ the effect of hydration on the total radius is found to be $\leq 2.5\%$ whereas the effect of differentiation into an isolated surface ocean is $\leq 5 \ \%$. Furthermore, we find that our results are very sensitive to the bulk composition.


Introduction
Water is abundant in the Universe; it is thought to have played an essential role in the onset of biochemistry on Earth and is indispensable for the sustainability of life as we know it (Kotwicki 1991;Kasting et al. 1993;Franck et al. 2001;Mottl et al. 2007;Jones et al. 2010;Güdel et al. 2014). These notions have led to a widely used definition of the classical habitable zone around a star. This zone is defined as the orbital region within which liquid water could exist on a planet's surface in terms of the stellar flux it receives and assuming simple greenhouse effects mediated mainly by CO 2 and H 2 O and maintained by the carbon cycle (Huang 1959;Rasool et al. 1970;Kasting et al. 1993;Kaltenegger et al. 2011;Ramirez 2018). This crude definition, although very useful as a first approach, cannot capture the full complexity of possible habitats as many other effects may be equally important for habitability. These include more diverse atmospheric compositions and chemical interactions, complex geochemical recycle processes, protection against harmful high energy particles from the host star via for example magnetic fields, and influences of the galactic neighbourhood (Ramirez 2018;Pierrehumbert et al. 2011;Gonzalez et al. 2001;Gonzalez 2005;López-Morales et al. 2011;Gallet et al. 2017;Javaux et al. 2010). Nevertheless, the presence of liquid water on a planet's surface is still commonly used as a proxy to constrain the orbital region around a star where life could exist. Hence, water is of the utmost interest for both astrophysics and astrobiology and has attracted increasing attention from the planetary science community in recent decades (for example Kuchner 2003;Léger et al. 2004;Seager et al. 2007;Sotin et al. 2007;Noack et al. 2016;Alibert 2014Alibert , 2016Fu et al. 2010;Kitzmann et al. 2015;Noack et al. 2017;Levi et al. 2017). It has been hypothesized that planets that formed beyond the ice line in the protoplanetary disk could migrate inwards into the habitable zone, where the accreted water ices would melt and form gigantic liquid surface oceans (Kuchner 2003;Léger et al. 2004). Indeed, in recent decades a large number of planets in the mass-range M 10M ⊕ have been detected with mean densities that are too low to correspond to a purely rocky composition (Rogers 2015;Lozovsky et al. 2018). These objects are not massive enough to have accreted substantial H-He envelopes in their past (Selsis et al. 2007;Sotin et al. 2007;Seager et al. 2007;Léger et al. 2004). This may suggest that their low densities are the results of extensive amounts of water on their surfaces. Such oceans could be maintained over timescales of several gigayears, long enough to allow for the emergence of life (Kuchner 2003;Noack et al. 2016). However, to sustain a biosphere, stable long-term climate Article number, page 1 of 34 arXiv:2012.06455v2 [astro-ph.EP] 16 Dec 2020 A&A proofs: manuscript no. aanda conditions are required. For instance, it has been proposed by Laskar et al. 1993 that the stabilization of the Earth's obliquity is crucial to avoid dramatic changes in climate that could have fatal consequences for life. These results were challenged later on by Lissauer et al. 2011, who found that the presence of a large moon has only a minor influence on the evolution of the obliquity for timescales of hundreds of millions of years. Furthermore, Armstrong et al. 2014 found large variations of obliquity to even extend the outer edge of the habitable zone. Consequently, on Earth, long term climate stability is probably more strongly supported by a negative feedback mediated by the carbon cycle between the atmosphere and the mantle (Walker et al. 1981;Post et al. 1990; Kasting et al. 1993). The presence of such a geochemical cycle has often been assumed to be a requirement for the long term sustainability of life on other worlds as well. However, the interior dynamics of ocean planets could be quite different from those of water depleted planets like the Earth. It is still debated whether geochemical cycles, such as the carbon cycle, could be significant on water worlds. The mantles of these planets are expected to be isolated from the gas envelopes by liquid water layers and further shielded by high pressure ice layers at the bottom of the ocean. As a result, the presence of massive surface oceans could inhibit chemical transport mechanisms between the atmosphere and the interior (Alibert 2014;Kitzmann et al. 2015;Noack et al. 2016;Abbot et al. 2012;Wordsworth et al. 2013;Nakayama et al. 2019).
Apart from its astrobiological significance, water is also an important ingredient in many solar system bodies including the ice giants Uranus and Neptune and the icy satellites of Jupiter and Saturn. The exact role of water in internal processes and its distribution in the interiors of these objects are subjects of active research efforts and there are many remaining open questions. The interior structures of Uranus and Neptune, for instance, are still unknown (see review by Helled et al. 2020). A number of theoretical studies show that various structure models with different assumed compositions can match the gravitational moments measured by the Voyager 2 flybys (for example Helled et al. 2020 and references therein). These studies further suggest that it is possible that the 'ice giants' contain substantially lower ice to silicate ratios than is commonly assumed, and that their interiors could be rock dominated (Helled et al. 2011). Uranus could even have solar ice to rock ratios (Nettelmann et al. 2016). It is debated whether the rocks and ices are separated into differentiated layers or form a mixed interior with gradual composition gradients (Helled et al. 2011;Nettelmann et al. 2013). In the mixed case chemical interactions between the Mg-silicates in the mantle or the iron in the core and the water ices could be important factors in further constraining the composition and thermal profiles of these planets. However, such hydration effects in the interiors are generally not included in such studies. A detailed investigation of these aspects could significantly improve our understanding of the interiors of water rich solar system planets and would allow for a more reliable exoplanet characterization in the future.
In the context of exoplanets, interior characterization is aided by rather limited data. The planetary radii and masses along with some stellar properties are the primary parameters that are remotely accessible for planets outside the solar system. From these measurements it is possible to gain insight into their compositions and hence formation and evolution histories. Early efforts to link observed masses and radii with theoretical models for the interior compositions have been presented for example by Léger et al. 2004and Sotin et al. 2007 over a decade ago. It is also worth noting the pioneering work by Zapolsky et al. 1969 long before the discovery of the first exoplanet. These authors have developed simple one-dimensional structure models to compute the total radii of planets as a function of mass, composition, and surface conditions. It has become common practice over the past years to use such models to interpret measured planetary masses and radii and predict possible interior structures and compositions. While such early approaches are quite powerful for inferring the general landscape of planetary properties, these models are highly degenerate. That is, many different internal compositions and profiles could match an observed M-R pair. This degeneracy has been noticed early on (for example Adams et al. 2007) and many researchers have devoted their work to minimize the possible parameter space for planetary interior properties for a given set of observed quantities (for example Rogers et al. 2010;Dorn et al. 2015Dorn et al. , 2017Grasset et al. 2009;Lozovsky et al. 2018). In these studies terrestrial planets are assumed to be fully differentiated into a number of compositionally distinct layers with water being modelled as an isolated surface layer. However, numerous studies over the past decades have shown that water can be incorporated in many minerals that are likely to be major constituents in the Earth's mantle. The upper mantle mineral (Mg, Fe) 2 SiO 4 (Olivine polymorphs) has been found to incorporate up to several wt% of water (see references in Tables E.1-E.3) in the correspondingly relevant temperature and pressure regime. At pressures ≈ 25-30 GPa, Ringwoodite (Rw), a high pressure polymorph of Olivine, dissociates into (Mg, Fe)O (Magnesiowüstite) and (Mg, Fe)SiO 3 (Perovskite) (Wang et al. 1997). MgO can react with water to form Mg(OH) 2 (Brucite) containing as much as ≈ 31 wt % H 2 O at moderate temperatures and pressures of 1500 K and 35 GPa (Hermann et al. 2016). Frost 1998 has determined the stability of the dense hydrous Mg-silicate Phase D and found it to be stable up to 50 GPa. This phase can incorporate up to 18 wt% of water. Nishi et al. 2014 have shown that Phase D transforms to an assemblage with another hydrous Mg-silicate, Phase H, at pressures above ≈ 48 GPa. They concluded that this phase could transport water in the form of hydrates into the lowermost mantle. Furthermore, Perovskite, which is thought to play a minor role for water storage in the lower mantle of Earth (Inoue et al. 2010;Bolfan-Casanova 2005), transforms into post-Perovskite (pPv) between ≈ 100-130 GPa and remains stable up to ≈ 0.8 − 0.9 TPa (Umemoto et al. 2011), approximately the relevant pressure regime in the interiors of Uranus and Neptune. Although the hydration behaviour of pPv remains unknown, density-functional theory (DFT) simulations suggest that stable configurations containing at least 2-3 wt% water exist up to pressures of 150 GPa and possibly higher (Townsend et al. 2015). Finally, Silica (SiO 2 ), which is stable over the entire pressure range relevant for the interiors of small to giant planets (Umemoto et al. 2011), has been reported to exhibit significant saturation contents of up to ≈ 3.2 wt% below 10 GPa and even up to ≈ 8.4 wt% for up to ≈ 100 GPa (Nisr et al. 2020 and references therein). In addition to the hydration of Mg-Silicates in the mantle, hydrogen solubility in the iron core might be equally important for constraining the total internal storage capacity of H 2 O equivalent in terrestrial planets. Indeed, it has been found that the solubility of hydrogen into iron is strongly enhanced at high pressures, relevant for the interiors of these objects (Fukai et al. 1983;Sugimoto et al. 1992). Although hydration effects under high pressure conditions have been actively studied by many researchers over the past decades, the Article number, page 2 of 34 O. Shah et al.: Internal water storage capacity of terrestrial planets and the effect of hydration on the M-R relation bulk water content of the mantle remains one of the most poorly constrained compositional parameters of the Earth (Townsend et al. 2015;Cowan et al. 2014).
The effects of hydration on the equations of state of the core and mantle materials could affect the thermal profiles and total radii of the planets. More importantly, the capacity of storing large amounts of H and OH in planetary interiors could lead to the distribution of the total amount of water between an internal reservoir and a surface ocean. For an assumed total mass fraction of H 2 O equivalent this would mean that the surface ocean, and hence the pressure at the mantle-ocean interface, could be considerably smaller than one would expect on a fully differentiated planet. Indeed, Cowan et al. 2014 have derived steady-state solutions to the water partitioning between hydrosphere and mantle on terrestrial planets and find, based on simplified assumptions for the mineral hydration, that super-Earths tend to store large amounts of water in their interiors. This can have relevant implications for the results of previous studies. For example, Sotin et al. 2007 have presented a model that takes the atomic ratios of Mg, Si and Fe in the interior as input parameters to generate mass-radius relations for given surface conditions. Grasset et al. 2009 have used this model to quantify the effect of a surface water ocean on the total radius of a planet with given bulk composition. They inferred that bulk composition plays only a secondary role for the total radius in comparison to the effect of the ocean mass fraction. The authors estimated that if the planetary mass and radius are precisely known, it is possible to determine the total amount of water with a standard deviation of 4.5 %. Alibert 2014 used a very similar model to compute the maximum radius of habitable planets assuming that habitability is hindered on water worlds by the formation of high pressure ice layers. It was concluded that for planets with 1-12 M ⊕ this maximum radius is limited to 1.7-2.2 R ⊕ and that larger planets are likely to be not habitable. These studies assume dry mantles and isolated surface oceans. However, in the light of the aforementioned hydration behaviour of Mg-silicates and iron, it is possible that these results are affected to a non-negligible extent when hydration is included. Understanding and quantifying these effects is essential given the increasing precision in mass and radius measurements for exoplanets. Although the mass and radius will never be precisely known, the uncertainties have decreased drastically over the past decades and will decrease further. To date, one of the most precise measurements of the radius of an exoplanet has been achieved for Kepler-93b (Ballard et al. 2014). Its radius has been determined from transit observations to with an uncertainty of ≈ ±1.3 %. The PLATO mission, scheduled for launch in 2026, is expected to determine the radii and masses of terrestrial planets with precisions up to ≈ 3 % and ≈ 10 %, respectively (Rauer et al. 2016). Hence, second-order effects that were justified to be excluded in the past are expected to become necessary ingredients for future planetary characterization.
In this study we build upon the models presented by Sotin et al. 2007, Grasset et al. 2009and Alibert 2014 and include the hydration of the mantle and the core as two separate reservoirs for H and OH and assess upper bounds for the corresponding effects on theoretical M-R relations. Furthermore, we apply our model to planets up to 3 M ⊕ with simplified bulk compositions and show how the effects of hydration on M-R relations can be estimated. The mass limit of 3 M ⊕ is dictated by the pressure range over which the core hydration model is valid. The generalization to more complex bulk compositions is briefly discussed. This paper is structured as follows: In Section 2 we describe how the internal profiles of the planets are constructed for variable boundary conditions and compositions. In Section 3 we explain how the H and OH content in the Mg-silicates and the core is calculated as a function of pressure and temperature. The resulting estimates on the effect of hydration and differentiation on total planetary radii for different bulk compositions as a function of total planetary mass are discussed in Section 4. Caveats arising from some simplifying assumptions about the composition and hydration in our model are outlined in Section 5 and a summary of the work is provided in Section 6.

The structure equations
The internal temperature, pressure, and density profiles need to be constructed to compute theoretical mass-radius relations. To this end the well known one-dimensional structure equations for spherical objects in hydrostatic equilibrium are integrated: Here, P(r), T (r), ρ(r), and m(r) are the pressure, temperature, density, and enclosed mass at radial distance r, respectively, G is the gravitational constant and (dT/dP) ad is the adiabatic gradient. To compute the density as a function of pressure and temperature adequate equations of state need to be employed (see Appendix B). The hydration of the mantle and the core, hereafter referred to in terms of mass fraction X H 2 O or the molar hydrogen content ξ H in the Mg-silicates and the iron is incorporated into the EoS (see Section 3).
The system of differential equations given by eq. (1) -(3) can readily be solved using a standard 4th order Runge-Kutta scheme with adaptive integration step size dr: where r is the refinement parameter. For an optimal trade off between stability and performance we recommend to use 0.1 ≤ r ≤ 0.5. All results presented here have been obtained with r = 0.25. The gradients have to be updated in each integration step to define the subsequent integration step size. With this procedure the planets are typically divided into a total of ∼ 100 − 1000 individual shells. For simplicity, we set the water content X H 2 O (P, T ) constant in each integration step. This induces a small error as P and T are of course different at the bottom of a shell and the top of the same shell. However, this effect can be controlled by the refinement parameter and can be neglected for our choice of r . Type-1 Type-2 Type-3 Fig. 1. Schematic overview of the different planet types modelled in this study. All objects consist of a pure iron or iron hydride core and a dry or hydrated Mg-silicate mantle. We have estimated the effect of hydration by comparing the dry planets (Type-1) with their hydrated counterparts (Type-2) at fixed composition and boundary conditions. We have then extracted the water content in the Type-2 planets and moved it into an isolated surface layer (Type-3) to assess the effect of ocean separation on the total radius at fixed composition and boundary conditions. The oceans of the Type-3 planets can consist of an upper liquid part and a lower solid part if the pressure at the bottom of the oceans is high enough to allow for the formation of high pressure ices.

Thermal model
The adiabatic gradient in eq. (3) is related to the material properties according to: where γ(T, P) and K S (T, P) are the Grüneisen parameter and the adiabatic bulk modulus, respectively. Eq. (6) is commonly written in terms of the logarithmic temperature derivative ∇ ad as: with The adiabatic bulk modulus can be related to the isothermal bulk modulus K T via the Grüneisen parameter γ and the thermal expansion coefficient α th using: K T and α th can be computed from the EoS (see Appendix B) directly: The P and T dependance of the Grüneisen parameter can be expressed in terms of the density: In eq. (6) -(12), T and P are functions of the radial distance T (r) and P(r), which we have omitted to write explicitly for clarity. Eq. (6) is only used for the iron in the core and the Mg-silicates in the mantle as for the pure water the adiabatic gradient (dT/dP) ad can self-consistently be extracted directly from the EoS used in this study (see Appendix B for details).
The corresponding values for the Grüneisen parameter for the different layers are summarized in Table 2.
In the case of mixtures of N materials denoted by the index i the density is obtained using linear mixing: where ξ i is the mole fraction of material i. The mole fractions are normalized to one, such that i ξ i = 1. The isothermal bulk modulus and thermal expansion are then given by: where we have omitted to write the T − P dependence of ξ i for clarity. Here we do not account for composition gradients and hence the derivatives ∂ξ i /∂P and ∂ξ i /∂T vanish.

Composition
The model planets are composed of the elements Fe, Si, Mg, O and H. Minor elements such as C, Ni, Ca, S or Al are excluded since they would significantly increase the complexity of the composition beyond the scope of this study (see also Sotin et al. 2007). We point out that the presence of these elements could affect the hydration behaviour of the Mg-silicates and therefore change their water content (see also Section 5). Based on the assumed elemental composition the modelled objects are divided into two to four layers: a pure iron core, an upper and lower Mgsilicate mantle with variable Fe content and possibly a surface water ocean. The transition from the upper to the lower mantle is defined via the dissociation of Ringwoodite (Rw) into Magnesiowüstite (Mw) and Perovskite (Pv) occurring at ≈ 25 − 30 GPa (Wang et al. 1997): If the pressure at the bottom of the mantle does not exceed this transition pressure, only an upper mantle is present. We do not account for a possible differentiation of the core into a liquid and solid part and neglect the possible presence of lighter elements other than hydrogen in the core for simplicity.
For the purpose of illustrating how our model can be applied to estimate the water contents in planetary interiors, we imposed three simplifying assumptions: (1) The mantles of the modelled objects consist of pure (Mg (1−ξ Fe ) , Fe ξ Fe ) 2 SiO 4 or (Mg 1−ξ Fe , Fe ξ Fe )SiO 3 + (Mg 1−ξ Fe , Fe ξ Fe )O in the upper or lower mantle, respectively. The choice of these minerals was motivated by the fact that they are major constituents in the Earth's mantle (Nunez-Valdez et al. 2013;Yao et al. 2012;Jacobsen 2006;Sinmyo et al. 2014;Li et al. 2014). Furthermore, Olivine (and its high pressure polymorphs) can incorporate considerable amounts of water and MgO can react with H 2 O to form Mg(OH) 2 (see references in Table E.1-E.3 and also Section 3). (2) The minerals in the interiors of the hydrated planets are fully saturated. (3) Since higher Fe content in Olivine polymorphs generally leads to an increase in the water saturation level, we chose ξ Fe = 0.25 as it roughly marks the upper limit for which our hydration model is valid (see Section 3). We point out that, for the simplified bulk composition, the density in the lower mantle can be lower than in the upper mantle in the fully water saturated case if only little Fe is present. This is due to the high water storage capacity of ≈ 31 wt% in Brucite (Mg(OH) 2 (Br)). This further motivates the choice of high Fe contents in the mantle. For our choice of ξ Fe this spurious behaviour could be avoided for all modelled cases and would naturally vanish in the Fe-free case if more Si-rich phases were present in the lower mantle.
In order to estimate the effects of hydration and ocean separation on M-R relations we consider here three types of planets: Dry planets (Type-1), hydrated planets (Type-2) and ocean planets (Type-3). The surface oceans are assumed to consist of pure water. An overview of the different types is provided in Fig. 1. The Type-1 planets have fully OH and H depleted interiors and no surface oceans. For the Type-2 planets the same boundary conditions, that is bulk composition and surface conditions, as for the Type-1 planets have been employed, but the mantles and cores are hydrated. A comparison between Type-1 and Type-2 planets allows for the assessment of the effect of hydration on total planetary radii. To estimate the effect of ocean differentiation we computed the amount of H 2 O equivalent in the Type-2 planets and added it as an isolated surface ocean on top of a dry mantle for the Type-3 planets. Since we are interested in estimating maximum effects in this study, intermediate cases where the water is partially distributed between an internal reservoir and a surface ocean, were not considered.

Boundary conditions
After defining the compositions of the layers the interior structure was determined by integrating from the centre to the outermost layer using initial guesses for the central pressure P C and temperature T C . For differentiated planets, the layer transitions must be assessed numerically using one of the integrated variables Q as a constraint, where Q is either the pressure, temperature, enclosed mass or radial extent of a layer or the total magnesium number of the planet. We used the core mass M Core to define the core mantle boundary (CMB) and the surface temperature T S to constrain the central temperature T C . The boundary conditions used in this study are: The total magnesium number Mg# ≡ [Mg]/([Mg] + [Fe]) including core and mantle, and the surface conditions P S , T S and M. The surface pressure is automatically matched as it defines the surface of the planet. T S , M and Mg# were iteratively probed by adjusting T C , P C and M Core until the desired precision for the boundary conditions was achieved. Since these quantities are not independent, we employed the following prescription for predicting a new set of initial values S in each iteration: with Q ≡ (1, Mg#, log(P S ), T S ), S ≡ (log(M Core /M), log(P C ), T C ).
We found eq. 17 to be sufficient in most cases to obtain the desired precision for all parameters and have hence omitted mixed terms in the Q i . In each iteration the coefficients a i j are obtained by solving the following matrix equation for k, i = 0, 1, 2, 3 and j = 0, 1, 2: M ki · a i j = x k j (20) The subscript k refers to the kth data point that is needed to constrain all a i j . In our case, four independent data points are required before the iterative prediction for the subsequent x k j can be carried out. These data points were obtained by roughly estimating the values S j for a given set of boundary conditions Q from a pre-constructed grid and then adopting slight variations to the S j . This procedure results in a very robust and efficient evaluation of the values S j with very high accuracy for a given Q. In the case of hydrated planets, the hydrogen content in the core ξ H is a priori unknown since it depends on the water content in the mantle and the temperature at the core mantle boundary (see Section 3.2 for details). However, as ξ H is a required input for the FeH x EoS in the core, it must be estimated prior to the integration. Here, we treated ξ H as a passive parameter in addition to the actively probed parameters S j . In each iteration the two previous data points were used to predict the new value for ξ H as a function of the new value for T C solely. With the procedure described above, ξ H could be matched with high accuracy if T S was probed with sufficient accuracy. We note that in general a less sophisticated iteration scheme, where the values S j are probed alternately via simple bisection, can be sufficient to achieve the same accuracy for the dry planets and the ocean planets. However, it failed for a considerable fraction of models for the hydrated planets. Furthermore, the procedure outlined above was found to be robust and more efficient than bisection. To determine when the quantity Q i has reached the desired precision we evoked the following criterion: where Q i is the desired precision for the parameter Q i . The same condition is used to probe the layer transitions using either the enclosed mass, the transition pressure or transition temperature as constraints. As the layer transitions in all our models depend on only one parameter, we probed the transition using simple bisection during the structure integration. The core-mantle transition has been probed using a fixed value for M Core , the upperlower mantle transition via a fixed value for the dissociation pressure of Olivine, P UM→LM and the mantle-ocean transition using the enclosed mass M − M Ocean at the bottom of the ocean. The corresponding values for the different parameters that have been used to define the boundary conditions and layer transitions are summarized in Table 1. We set higher accuracy for M of the Type-3 planets and T S of the Type-2 planets. This was necessary to achieve the desired precision for M Ocean /M and ξ H , respectively. The transition criteria for the different layers are listed in Table 2. The specified accuracy was reached in about 80 − 90 % of the modelled cases. For the rest of the planets the achieved accuracy for M, Mg#, ξ H , and M Ocean /M were at least 10 −3 , 5 · 10 −3 , 3 · 10 −2 , and 2 · 10 −2 , respectively. These numerical uncertainties have, however, minor effect on the main results presented in this study. To obtain the mass-radius relations it is sufficient to vary the central pressure P C to cover a desired mass range. However, since we aim to compare different mass-radius curves directly we explicitly fixed M via P C for each of the modelled planets individually. This allows a direct comparison of the planetary properties for a given mass, composition, and surface conditions. In the case of dry mantles, where the composition in all layers is constant with depth, the core mass fraction can be computed analytically as a function of the total bulk and mantle compositions Mg# and ξ Fe , the total mass M and the ocean mass M Ocean : where m Fe the molar mass of iron.m is the mass of a portion of the material in the dry mantle containing one mole of Fe. It is defined as:m The oxygen number O# is given by: The mantle mass fraction can then simply be computed for a given total mass and ocean mass fraction as In these cases the iteration for the core mass to match Mg# can be omitted. The core mass fractions and corresponding mantle mass fractions as a function of the mantle composition are plotted in Fig. 2 for different bulk compositions and ocean mass fractions.
For differentiated planets each individual layer can have a different convection behaviour based on the composition and the P − T -profiles. Hence, the heat transport via convection from the interior can be partially interrupted at a layer transition. This leads to a temperature drop upon transiting from a convective layer into a less convective one. Stixrude 2014 have presented models for the temperature profiles of super-Earths considering thermal regulation controlled by silicate melting. It was found that the temperature drop at the CMB denoted as ∆T CMB , can be represented by a simple scaling law as a function of the total planetary mass M: Here, we use eq. 27 to compute ∆T CMB over the modelled mass range and for all bulk compositions. We did not imposed additional temperature jumps in the mantle transition zone (MTZ). We treat the temperature drop ∆T TBL in the thermal boundary layer (TBL) at the top of the Mg-silicate mantle (mantle-surface or mantle-ocean interface) as a free parameter and have considered values of 200 K, 700 K, 1200 K and 1700 K. The value ∆T TBL = 1200 K roughly corresponds to modern day Earth with T S ≈ 300 K (for example Sotin et al. 2007). Since the water saturation in the mantle Mg-silicates drops with increasing temperature (see Section 3), the maximum effect at given composition and surface conditions is reached when the temperature discontinuity in the MTZ vanishes. We note, however, that a non-vanishing temperature jump in the MTZ could in principle increase the hydration of the iron core (see Section 3.2). Furthermore, we only consider one value T S = 300 K because different values for T S could, in principle, be interpreted as different values for ∆T TBL . For instance, T S = 300 K and ∆T TBL = 1700 K could also be interpreted as T S = 800 and ∆T TBL = 1200 K. Furthermore, at elevated surface temperatures the definition of the Type-3 analogues would be less straightforward as part of the ocean mass fraction would enter the vapor phase. Since we do not account for atmospheres in this study, we explicitly fix T S = 300 K for the Type-3 planets.
In this study, we computed the mass-radius curves for fixed surface pressures and temperatures of P S = 1 bar and T S = 300 K. The Mg# has been varied between 0.2 − 0.7. This is roughly the range covered by estimates for elemental ratios of planet hosting stars presented by Grasset et al. 2009 (see their Fig. 2). We point out that for their presented value range of Mg/Si and Fe/Si, the magnesium number could be as high as ≈ 0.8. However, in this study Mg# must be strictly smaller than 1 − ξ Fe = 0.75 for the assumed Fe content in the mantles of ξ Fe = 0.25 if a core is present.

Updating the shell contents
We assume the hydrogen content in the core to be homogeneous. For a given composition the water content in the mantle, however, is a function of temperature and pressure (and hence of radius). The water content was updated in each shell for the  hydrous phases before integration according to the water storage capacity of the Mg-silicates (see Section 3). A change in the molar water content, ξ H 2 O , leads to a change of the ratios Si/Mg. Therefore, in order to keep Si/Mg at a fixed value in each shell, the molar fractions of the different minerals need to be updated according to the water content. Here, this is irrelevant for the upper mantle since it consists of only one stoichiometry, the Olivine polymorphs, and hence its molar abundance, denoted by ξ Ol , is always equal to unity. In the lower mantle two minerals are present: Pv or pPv and Mw with molar Fe content ξ Fe . If water is added, the mole fractions of Pv or pPv and Mw ξ Pv and ξ Mw must be updated. The molar abundance of Mw is given by: The molar abundance of Pv/pPv is then simply given by: Here ξ H 2 O,Mw is either 0 or (1 − ξ Fe )/2, depending on the location in the phase diagram of Mg(OH) 2 (see Section 3). For Pv/pPv we assume constant water content of 0.1 wt% (Pv) and 3.0 wt% (pPv) (see Section 3). We have assumed the hydrogen content in the core to be homogeneous with depth and hence ξ H remains constant throughout the integration of the core.
2.6. Extracting the deep water reservoirs 2.6.1. In the core The total amount of H 2 O equivalent in the core is given by the hydrogen content in FeH x , which is assumed to be homogeneous in the entire core. We discuss how the parameter x can be computed in Section 3.2. Once the stoichiometric hydrogen content x in the core is known, it can easily be converted to a molar abundance ξ H as: Article number, page 7 of 34 A&A proofs: manuscript no. aanda The total atomic hydrogen content in the core is then given by: The H 2 O equivalent mass in the core can then easily be obtained: 2.6.2. In the mantle The local water weight fraction X H 2 O (P, T ) in the mantle depends on the radial distance via P(r) and T (r). As outlined in Section 2.5, we assumed X H 2 O to be constant in each shell. The total mantle reservoir in a planet is then simply given by summation over all water containing shells. For a given shell, the water content M H 2 O,Shell was extracted from the shell properties using the following statement: where N H 2 O,Shell and m H 2 O are the molar amount of water in the shell and the molar mass of water, respectively. N H 2 O,Shell is related to the shell mass M Shell and the molar water contents ξ H 2 O,i in the Mg-silicates (hereafter denoted as Mg-Si): with the relative molar masses given by: Here, n Mg−Si denotes the number of coexisting Mg-Si phases in the shell. ξ H 2 O,i can directly be computed from X H 2 O,i by invoking mass balance: and hence:

The mantle
Here we present a simple model for water saturation in Fe bearing Olivine polymorphs (upper mantle) and Brucite and Perovskite or post-Perovskite (lower mantle). We choose these minerals because they are known to be major constituents of the Earth's mantle, where they are thought to play significant roles as hydrated Mg-silicate (see Section 2.3). The upper limit of hydration in the Mg-silicates is in general given by the saturation water content, which must be either measured in experiments or computed using molecular dynamics simulations.
In an astrophysical framework, however, the extent of internal hydration of a planet cannot be assessed directly. That is, we cannot constrain it from stellar composition and total mass and radius alone. Thus, the water content in the mantle is written as ] is a free model parameter. It can be constrained from geophysical studies of terrestrial planets. Indeed, numerous authors have investigated the influence of water on the interior evolution of small objects and super-earths and how the internal water content in mantle Mg-silicates evolves under the influence of a variety of geochemical and physical processes. For Mars, for instance, it has been shown that water loss from the interior over its evolution has been minor and that most of the initial water content in the mantle should still be present today (Hauck et al. 2002;Ruedas et al. 2013). Nakagawa et al. 2015 modelled water circulation mechanisms for the Earth and found that the water content in an initially anhydrous mantle increases over time as a result of hydration of subducting slabs. They found that a balance between regassing and degassing at shallow depths is reached, which dictates the total water content in the mantle. Nakagawa 2017 showed that in the case of strongly water-dependent viscosity water transport into the lower mantle is more effective, which can lead to higher internal water contents. Tikoo et al. 2017 modelled magma ocean solidification with variable initial water contents for the young Earth and showed that the water can be retained at saturation levels in the mantle upon cooling. Furthermore, Cowan et al. 2014 investigated water cycling between mantle and ocean and found that super-Earths tend to absorb more water in their mantles as the seafloor pressure is proportional to the surface gravity. These studies suggest that whenever a sufficient amount of water is present during planetary evolution, it is likely that the Mg-silicate mantle of the planet either remains or becomes strongly hydrated, which means H 2 O ≈ 1. We are mainly interested in ocean planets that contain between a few and tens of wt% of water for which the condition of 'sufficient' amount of water is likely met. For the cases of lower water contents we still adopted H 2 O = 1 to estimate the maximum possible effects of hydration. Therefore, whenever we refer to hydrated planets we shall assume fully water saturated Mg-silicates in the mantle, that is X H 2 O (P, T ) = X H 2 O,sat (P, T ), unless otherwise stated.
Our prescription for water saturation in Olivine polymorphs in the upper mantle is based on a collection of experimental data that has been acquired by various groups over the last decades. We construct a simple tri-linear model from a least squares fit to compute the saturation water content X H 2 O,sat (T, P, ξ Fe ) as a function of temperature T , pressure P and iron content ξ Fe : A full list of the collected data that has been used for the fit, along with the corresponding references, is provided in Hydrogen in the Core ∆T TBL = 1700 K Fig. 3. Partial hydrogen pressure P H 2 (left column) and resulting hydrogen content (right column) in the cores of the modelled Type-2 planets. We have limited our models to M ≤ 3 M ⊕ to ensure P H 2 1 GPa as the hydrogen solubility from eq. 41 has been found to be considerably underestimated at higher pressures. The hydrogen content is computed from chemical equilibrium between the core and the mantle. It depends on the temperature at the CMB and the partial hydrogen pressure in the lowermost part of the mantle (see text for details).
the α-phase. However, since the water content in the α-phase spans a large range depending primarily on pressure and iron content, the effect of ξ Fe on X H 2 O,sat is misrepresented at elevated pressures if c 011 and c 101 are excluded from the fit. Furthermore, at very high temperatures or very low pressures our fit can yield negative values for X H 2 O,sat . We therefore invoke an artificial lower threshold value of X H 2 O,sat ≥ 0 to avoid a non-physical behaviour.
In the lower mantle we chose Brucite (Mg(OH) 2 ) and Perovskite (MgSiO 3 ) as the water bearing phases. We fix the water content in hydrous Perovskite at 0.1 wt% and neglect effects on the EoS for such low water contents (Inoue et al. 2010;Bolfan-Casanova 2005). For post-Perovskite it has been found that configurations containing more than 2 wt% of water are stable (Townsend et al. 2015). However, the pressure and temperature dependence of the water content has not been constrained. For this reason we adopt a constant value of 3 wt% (see also Appendix B.4). The stability field of Brucite is limited to rather low temperatures and pressures. It dissociates into MgO + H 2 O above ≈ 30 − 35 GPa (Hermann et al. 2016). We extracted the stability field of Mg(OH) 2 in the P − T plane from Fig. 2 in Hermann et al. 2016. As a result, the lower mantle can itself be divided into two parts: an upper part where Mg(OH) 2 is present and a lower part where only MgO is present.

The core
We assume that the main source of hydrogen that is dissolved in iron and eventually deposited into the core comes from the Mgsilicates and enters the iron via the reaction (Wu et al. 2018): In an early magma ocean hydrogen dissolves in liquid iron droplets and is incorporated into the core during the segregation process. If we assume that iron hydride, FeH x , droplets are in chemical equilibrium with the Mg-silicates at all times, then the hydrogen content in the iron is dictated by the conditions at the CMB. This notion allows us to compute the hydrogen content in the core by adopting chemical equilibrium between the core and the lowermost mantle. If the water content in the mantle changes with time, either through external sources of water or the thermal evolution of the mantle, the hydrogen content in the core would adapt in order to maintain chemical equilibrium. Here we employ a simple model for the H-solubility as a function of the temperature at the CMB and the partial pressure of H 2 in the Article number, page 9 of 34 A&A proofs: manuscript no. aanda Type-1   lowermost part of the silicate mantle following Wu et al. 2018. These authors have used Sieverts' law for gas solubility in metals to estimate the hydrogen concentration ξ H in the iron cores of planetary embryos. We refer the interested reader to their Section 2.3 for further details. For our purpose it is sufficient to state the prescription for the hydrogen solubility in pure iron, which is given by: where P 0 = 1 bar is the reference pressure, P H 2 is the partial pressure of H 2 in the silicates at the CMB, R is the universal gas constant and ∆H 0 = +31.8 kJ mol −1 and ∆S 0 = −38.1 kJ mol −1 K −1 are the enthalpy and entropy for dissolution of hydrogen into iron, respectively. P H 2 is given by the amount of hydrogen present in the silicate mantle at the CMB. Here we used the van-der-Waals (VDW) EoS to relate the density of H 2 to the partial pressure at given temperature. The Sieverts law is well suited to compute gas solubility in metals over an extended pressure range. However, Sugimoto et al. 1992 suggested that the Sieverts law underestimates the hydrogen solubility in iron at partial pressure P H 2 > 1 GPa. Therefore, we restrict our study to planets with M ≤ 3 M ⊕ , which keeps P H 2 strictly below ≈ 1 GPa for all modelled cases (see. Fig. 3).

Results and discussion
In the following we present the results for the total storage capacities of H 2 O equivalent (including core-and mantle reservoirs) in the mass range 0.1 ≤ M/M ⊕ ≤ 3 for different bulk compositions, 1 bar surface pressure, 300 K surface temperature and ∆T TBL between 200 K and 1700 K. We also compare the different types of planets as described in Fig. 1 and estimate the net effect on the total radii as well as the internal density, pressure and temperature profiles. Finally, we discuss possible implications of our results to exoplanet characterization.

Water storage capacity and mass-radius relations
Here we address three main questions to illustrate the application of our model: 1. What is the maximum amount of H 2 O equivalent that a terrestrial object of a given size, composition, and surface conditions can store in its interior?     To answer these questions we modelled three types of planets in the mass range 0.1 ≤ M/M ⊕ ≤ 3 by solving the structure equations for adequate boundary conditions (see Section 2 for details). An overview of the different types of planets is given in Fig. 1. We modelled a total of 1728 planets for four different values of ∆T TBL , six different Mg#, and twenty four values for M uniformly distributed in logarithmic space in the range 0.1 ≤ M/M ⊕ ≤ 3. Fig. 4-6 show the central temperatures T C (first column), central pressures P C (second column) and the mass-radius relations (third column) obtained over the entire range of compositions and masses. The mass-radius relations for pure H 2 O, MgO and Fe at P S = 1 bar, T S = 300 K are also shown for reference. We find that while the central temperature changes significantly between the different types, the central pressure is relatively unaffected. This is because the adiabatic temperature gradient is steeper in hydrated Mg-silicates than in the anhydrous case, which leads to a net increase in the internal temperature profile for the Type-2 planets with respect to the Type-1 and Type-3 planets. Furthermore, Fig. 5 shows a slight change in the slope of the central temperature between ≈ 0.1 − 0.3 M ⊕ corresponding to the positions of local maximum in the internal water content (see also discussion below and Fig.  7). Larger core mass fractions generally lead to larger central temperatures because the temperature gradient is steeper in the core than in the Mg-silicates. A similar behaviour is obtained for the central pressure. The central pressure decreases with Mg# for all planets because the pressure gradient is steeper in the core than in the Mg-silicate mantle (see also Section 4.2). That is, for a given surface pressure and total mass, increasing the core mass fraction (or equivalently decreasing Mg#) must be compensated by an increase in the central pressure. The local maximum of the total water content are also reflected in the curves for the central pressure by slight changes of the slope between ≈ 0.1 − 0.3 M ⊕ . The mass-radius relations are very similar for all planet types and show consistently larger radii for larger Mg# due to the lower density of the Mg-silicates with respect to the iron in the core. Furthermore, in Fig. 6 6. Central temperature (left column), central pressure (middle column), and M-R relation (right column) for the Type-3 planets. The different rows correspond to different assumed values for ∆T TBL (see annotations at the right). While the core temperature increases significantly for higher ∆T TBL the central pressure and radius remain vastly unaffected. The kinks in the mass-radius curves correspond to the Pv-pPv transition, which leads to a drastic increase in the core reservoirs in the Type-2 planets. Accordingly, the ocean mass fractions in the Type-3 planets increase in the same mass range, which decreases the mean density of the palents. The mass-radius relations for pure H 2 O, MgO and Fe are shown for reference.
planets starts to increase drastically (see also next paragraph). This leads to a significant increase in the ocean mass fraction and hence reduction of the mean density in the Type-3 planets in the same mass range.

Comparison between Type-1 and Type-2 planets
In Fig. 7 we show the total water content (left column), the coreand mantle reservoirs (middle column) of the Type-2 planets and the resulting relative change in radius between Type-1 and Type-2 planets (right column). As expected, the fractional water content in the mantle increases strongly with the total Mg#. This is because larger mantle mass fractions have the effect of increasing the pressure range in the mantle, which leads to higher water saturation contents in Olivine polymorphs (see Fig. B.2). In addition, higher mantle mass fractions increase the amount of Mg-silicates with respect to the total planetary mass. Furthermore, the mantle reservoirs scale with the planetary mass only as long as an upper mantle is present. If the mass is sufficiently large that the transition pressure Ol → Pv + Mw is reached, more material from the lower mantle with lower water content is added. This is reflected in the distinct peaks in the total water content and the mantle reservoirs between 0.1 M/M ⊕ 0.3. The peak is highest for Mg# = 0.7 and ∆T TBL = 200 K at ≈ 3.8 wt% and decreases strongly for lower Mg# or higher ∆T TBL . In contrast, at higher masses the total water content is dominated by the core reservoir, which leads to an opposite trend. This is because the core reservoir scales with the core mass fraction and hence scales inversely with the Mg#. The core reservoirs themselves exhibit local maximum between 0.1 M/M ⊕ 0.3 where the transition from the water rich upper mantle to the water poor lower mantle occurs. These maxima coincide with the local maxima in the mantle reservoirs (we note that this is only roughly true in the Fig. 7 because of the finite resolution of the curves). At higher masses 0.6 M/M ⊕ 2 the pressure at the bottom of the lower mantle exceeds the transition pressure between Pv and pPv leading to higher water contents in the Mg-Silicates at the CMB and hence higher equilibrium hydrogen content in the core. This is reflected in an increase in both core and mantle reservoirs for reservoirs at 3 M ⊕ of up to ≈ 3 wt% for ∆T TBL = 200 K and ≈ 6 wt% for ∆T TBL = 1700 K. The relative change in radius upon hydration is shown in the third column of Fig. 7 and exhibits three major trends: First, independent of temperature, the effect of hydration on the mean density is always strongest roughly in the range 0.1 M/M ⊕ 0.3 where the mantle reservoirs peak despite the fact that for higher temperatures, the total fractional water content is much larger at higher masses. This means that hydration manifests more strongly for the upper mantle minerals than the iron core.
Second, the heights of the peaks clearly depend on ∆T TBL . Between 200 K and 700 K the maximum δR 1 /R 1 decreases from ≈ 0.012 to ≈ 0.01, remains roughly constant between 700 K and 1200 K and strongly increases between 1200 K and 1700 K up to ≈ 0.025. To understand this behaviour several effects must be taken into account: (a) The water content in the upper mantle typically decreases with temperature. (b) The hydrogen content in the core increases with the water content in the upper mantle, provided there is no lower mantle. (c) The hydrogen solubility in the core increases with temperature. (d) Hydration tends to steepen the adiabatic temperature gradient leading to generally hotter interiors for the same boundary conditions, which increases the radius. Effects (a) and (b) lead to a decrease in the total fractional water content with temperature, which generally leads to a smaller change in the radius. Effect (c), on the other hand, leads to an increase in the total fractional water content with temperature. Hence, both effects (c) and (d) contribute positively to the relative change in radius while effects (a) and (b) have a negative contribution. The observed trend in the heights of the peaks therefore means that below 700 K effects (a) and (b) dominate while above 1200 K effects (c) and (d) are more profound. In the intermediate range between ∼700-1200 K, the negative and positive contributions balance each other.
Finally, the third trend is that the Mg# for which δR 1 /R 1 is largest at given mass and temperature corresponds to the value of Mg# for which the total fractional water content is largest. Exceptions of this trend can be seen for ∆T TBL = 200 K below ≈ 0.8 M ⊕ . However, these small fluctuations are likely to have a numerical rather than a physical origin.

Comparison between Type-2 and Type-3 planets
The comparison of Type-2 and Type-3 planets is summarized in Fig. 8. The ocean mass fractions (first column) for the Type-3 planets are calculated from the total water content in the Type-2 planets (compare first column in Fig. 7). The corresponding ocean depths are shown in the second column of Fig. 8. As expected, the ocean depths exhibit overall similar behaviour as the ocean mass fractions. But for lower masses a fixed ocean mass fraction results in generally lower ocean depths in comparison to the same ocean mass fraction at higher total mass. This is because ρ R 3 ∝ M for spherical objects (see Appendix C for further discussion). We find that the estimated maximum H 2 O equivalent mass for all modelled planets corresponds to ocean layers with depths up to ≈ 800 km for M/M ⊕ = 3 and Mg# = 0.2. The effect of ∆T TBL on the ocean depth is twofold: higher temperatures reduce the water content in the mantle but can increase the core reservoir of the Type-2 planets (see Fig. B.2 and 3). Because the water content in the low mass planets is dominated by the upper mantles, this results in a downwards trend of the ocean depth with increasing ∆T TBL . Since X H 2 O in pPv is constant, at higher masses only the positive effect of ∆T TBL on the core reservoir manifests. This explains why the ocean depths scale upwards monotonously with temperature for M/M ⊕ 0.8 but not for lower masses. However, we remind the reader that the water content in pPv is in reality likely to change with P and T . A more realistic model for the hydration of pPv could therefore change the aforementioned behaviour.
The relative difference in radius δR 2 /R 2 ≡ (R 3 − R 2 )/R 2 between the hydrated Type-2 planets and the fully differentiated Type-3 planets is in general larger than δR 1 /R 1 for a given mass. The maximum δR 2 /R 2 can reach up to ≈ 5 % for M/M ⊕ = 3 and Mg# = 0.2. For temperatures ∆T TBL > 700 K the δR 2 /R 2 can in our model framework even become negative. This means, that differentiation into a surface ocean results in a net decrease in the planetary radius. Although counterintuitive at first glance, this feature does in fact arise intrinsically from the properties of the EoS used in this study. The temperature gradient in pure water tends to be shallower than in the hydrated silicates. As a result there is a considerable reduction of the interior temperature profiles between the Type-2 and Type-3 planets of up to a few thousand Kelvin (compare Figs. 5 and 6 or see Section 4.2).  ∆T TBL = 1700 K Fig. 9. Internal profiles as functions of radial distance from the centre for M/M ⊕ = 0.1 for all modelled ∆T TBL and compositions for the temperature (left column), pressure (middle column), and density (right). The solid curves correspond to Type-1 planets, the dashed curves to Type-2 planets and the dotted curves to Type-3 planets (see text for details). Since hydration tends to steepen the temperature gradient in the mantle, the Type-2 planets exhibit generally higher temperatures in the interiors than their Type-1 and Type-3 analogues.
The densities of the cores and mantles of the Type-3 planets are hence enhanced due to the much lower temperatures and the absence of H or OH, respectively. For the 200 K and 700 K planets, this effect is, however, very small, explaining why in this case all Type-3 planets are indeed inflated with respect to their Type-2 counterparts due to the presence of the low density surface ocean. At higher values of ∆T TBL the effect becomes large enough to counteract the decreasing effect on the mean density of the surface ocean. This can lead to a decrease in the total radius unless enough water is present. For large ocean mass fractions the decreasing contribution of the low density ocean dominates and the Type-3 planets are always larger than the Type-2 planets in the super-Earth regime. It is interesting to note that the change of sign of δR 2 /R 2 roughly lies at 1 M ⊕ for solar composition (Mg# 0.5). This is because the Pv-pPv transition occurs at roughly the temperature and pressure conditions near the bottom of the Earth's lower mantle.
Our findings indicate that the change in the total radius due to hydration is very likely below the detection limit of current or  ∆T TBL = 1700 K Fig. 10. Internal profiles as functions of radial distance from the centre for M/M ⊕ = 0.6 for all modelled ∆T TBL and compositions for the temperature (left column), pressure (middle column), and density (right column). The solid curves correspond to Type-1 planets, the dashed curves to Type-2 planets and the dotted curves to Type-3 planets (see text for details). Since hydration tends to steepen the temperature gradient in the mantle, the Type-2 planets exhibit generally higher temperatures in the interiors than their Type-1 and Type-3 analogues.
near future instrumentation for characterizing exoplanets. Furthermore, the maximum effect of hydration on the radius occurs in a mass range that is well below 1 M ⊕ . The partitioning of water between internal and surface reservoirs, on the other hand, can indeed affect a planet's radius within plausible limits of observations in the foreseeable future. As this effect is enhanced for larger planetary masses, it is of particular interest for future characterization of planets in the super-Earth regime. However, here we used a relatively simple model for the hydration of the cores and mantles and have assumed a rather simple mantle composition. We discuss the limitations of our model in Section 5.

Internal profiles
In Figs. 9-10 we show the internal temperature (left columns), pressure (middle columns), and density (right columns) profiles as functions of radial distance from the centre for three selected planetary masses. The solid, dashed, and dotted lines correspond to the Type-1, Type-2, and Type-3 analogues, respectively. The ∆T TBL = 1700 K Fig. 11. Internal profiles as functions of radial distance from the centre for M/M ⊕ = 3 for all modelled ∆T TBL and compositions for the temperature (left column), pressure (middle column), and density (right column). The solid curves correspond to Type-1 planets, the dashed curves to Type-2 planets and the dotted curves to Type-3 planets (see text for details). Since hydration tends to steepen the temperature gradient in the mantle, the Type-2 planets exhibit generally higher temperatures in the interiors than their Type-1 and Type-3 analogues. For larger masses, the temperature profiles of the Type-3 planets can, however, exceed the temperature of the Type-2 planets for low values of ∆T TBL (see top row).
pressure profiles are very similar for all planet types. They are shifted to slightly lower pressures between Type-1 and Type-2 planets but are essentially indistinguishable by eye between Type-1 and Type-3 planets. The reason is the steeper pressure gradient in the core. As the mantle is hydrated, the ratio between total mantle mass (Mg-silicates plus water) and the core mass becomes larger. Therefore, for a fixed total mass and surface pressure, the central pressure needs to be reduced. Since the pressure difference in the surface oceans is very small in compar-ison to the central pressures no similar effect occurs between the Type-1 and Type-3 planets. The temperature profiles are much more sensitive to the model type and this sensitivity increases with increasing mass. The profiles shift to higher temperatures between Type-1/Type-2 planets. This is because hydration tends to steepen the temperature gradient in the mantle. Furthermore, the temperature profiles of the Type-1 and Type-3 planets are indistinguishable for 0.  ). For large masses, the ocean mass is sufficient to reach pressures at the bottom of the oceans that exceed the transition pressure from liquid water to high pressure ices. This can also be seen in the density profiles in Fig. 11. and become deeper for higher masses. As deeper oceans generally lead to larger temperature differences between the surface and the oceans, the internal temperature profiles at given surface temperature are more affected for deeper oceans (compare also Fig. 12). An interesting feature in Fig. 11 is that the temperature profiles in the Type-3 planets are higher for ∆T TBL = 200 K than even in the Type-2 planets. For higher values of ∆T TBL , however, the difference between the temperature profiles of the Type-3 planets and the Type-2 or Type-1 planets decreases. At very high temperatures the profiles of the Type-3 planets eventually fall even below the Type-1 planets. Considering the effect of the surface ocean on the temperature profiles alone cannot explain this behaviour. It would be intuitive to assume that larger oceans (larger temperature differences between surface and mantle) would monotonically increase the temperature in the mantle and the core. However, large ocean mass fractions also increase the pressure at the ocean-mantle interface. While the temperature gradient in the mantle becomes generally steeper for higher temperatures, it becomes shallower for higher pressures. For lower values for ∆T TBL the temperature effect dominates, rendering the temperature gradient steeper in the interior and hence resulting in somewhat higher temperature. For larger values of ∆T TBL , however, the pressure effect takes over, leading to a net decrease in the temperature gradient. The density profiles are very similar for all planet types. The density of the mantles and cores of the Type-2 planets are reduced compared to the Type-1 planets because of the hydration effects on the EoS. The Type-1 and Type-3 planets have, apart from the surface oceans, nearly identical density profiles. Fig. 12 shows the pressure at the bottom of the oceans (left column) and the temperature difference between the surface and the bottom of the oceans (right column) for Type-3 planets. The ocean pressures reach up to ∼ 0.01 − 10 GPa where P Ocean 1 GPa is only reached for M 1 M ⊕ . This means that for the majority of the modelled objects the oceans are fully liquid (compare also Fig. B.1). For higher masses both the pressure and the temperature at the bottom of the ocean significantly increase. The temperature increase between the surface and the bottom of the ocean reaches up to ≈ 500 K, corresponding to temperatures at the bottom of the ocean of ≈ 800 K. The transition pressures from liquid water to the high pressure ices in the temperature range ≈ 300 − 800 K roughly lie in the range ≈ 1 − 10 GPa. This suggests that for larger masses the oceans of the Type-3 planets can posses an upper liquid layer and a lower solid layer (this is also evident in the density profiles shown in Fig. 11). These results indicate that the amount of H 2 O that can be stored in the interiors of terrestrial planets could reach up to ocean mass fractions sufficient for the formation of high pressure ice layers below the liquid parts of the oceans. This means, we find that large total water mass Article number, page 18 of 34 O. Shah et al.: Internal water storage capacity of terrestrial planets and the effect of hydration on the M-R relation fractions of up to a few wt % could in fact easily be partitioned into interior and surface reservoirs so that the remaining surface oceans, that would in isolation be sufficiently large to form high pressure ice layers, would be liquid throughout. For low enough surface temperatures it is in principle also possible to have an ocean that is frozen throughout, regardless of the total ocean mass fraction. This scenario has, however, not been considered in this study.

Implications for previous studies
We have found that planets with considerable water content of at least a few wt% could still have only very small or no surface oceans. This could have important implications in the context of previously published studies. The maximum radii for habitable planets derived by Alibert 2014, for instance, are likely to be underestimated. That is, the maximum total amount of H 2 O equivalent for which no high pressure ice layer forms at the bottom of the ocean could be considerably larger. This would lead to larger maximum radii for habitable planets.
Furthermore, our results could alter the findings of Grasset et al. 2009. These authors have estimated that the total ocean mass fraction could be derived from the bulk composition, the mass and the radius of a planet alone with a standard deviation of no more than 4.5 % if the mass and radius were to be measured precisely. However, our study shows that the partitioning of the total water content into internal and surface reservoirs can have considerable effects on the radius of a planet. In particular, for super-Earths the partitioning between surface and internal reservoirs would decrease the total radius with respect to a fully differentiated ocean planet. That means that the mean density of such a planet would become more similar to the mean density of a fully water depleted planet and it becomes more difficult to distinguish between these objects. In this light, it has to be assumed that the determination of the water mass fraction present on an observed exoplanet from mass-radius measurements is less precise than suggested by Grasset et al. 2009.

Caveats
First, a number of simplifying assumptions for the bulk composition of the planets have been made (see Section 2.3). Including a larger variety of mantle minerals could induce a different overall hydration behaviour. For instance, Pyroxene plays a dominant role along with Olivine in the upper mantle of the Earth (Jacobsen 2006). The coexistence of different minerals can change the water content in the end member phases (Grant et al. 2006). Furthermore, SiO 2 has been found to have large water storage capacity over a very wide pressure range and could be a major agent for water storage under lower mantle conditions. Nisr et al. 2020 have studied hydrous SiO 2 samples up to ≈ 100 GPa and find water contents of 0.4 − 8.4 wt%. We have not included the possible presence of SiO 2 in our models. In addition to the exclusion of some minerals that are thought to be important constituents in the earths mantle, the hydration behaviour of some of the minerals are very weakly constrained. Most importantly, the saturation content in pPv is only poorly studied and to our knowledge, no data on the hydration behaviour as a function of pressure and temperature is currently available. The simplified assumption of constant water content of 3 wt% over the entire pressure and temperature range is likely to be unrealistic. The water content in pPv dictates the equilibrium hydrogen content in the core and is therefore a key parameter in our models and presumably one of the major sources of uncertainty at this stage.
These simplifying assumptions on the composition limit the interpretability of our results with respect to specific objects. For example, when assuming roughly solar [Fe]/[Mg] ≈ 0.5 for the bulk Earth, our results indicate that a maximum of ≈ 1.4 wt% of water could be stored in the Earth's mantle (see Fig. 7). This is much higher than standard estimates for the actual water content. From isotopic mass balance the total amount of water in the Earth's interior at present is found to be ≈ 0.06 − 0.3 wt% (see review by Marty 2012). Geochemical and geophysical considerations lead to generally somewhat lower estimates for the deep water reservoirs of ≈ 0.006 − 0.07 wt% (see review by Ohtani 2019). Furthermore, Cowan et al. 2014 have developed simple two-box models for the water cycling between ocean and mantle and derived the internal water content of terrestrial planets based on the equilibrium between degassing and regassing. They reported a steady state solution for Earth with an internal water content of ≈ 0.08 wt%. One reason for the large discrepancy between our model and these results is, of course, the fact that we assumed water saturated Mg-silicates in the mantle. However, the Earth's mantle is likely not fully water saturated (Ohtani 2019) and hence we do not expect our model to reproduce the actual water content but rather to give an upper limit. Another important factor is the fact that the ratio [Si]/[Mg] takes a fixed value of [Si]/[Mg]=2/3 for our simplified composition. This is much lower than the solar value of [Si]/[Mg] ≈ 0.9 (Sotin et al. 2007). Adjusting the composition to roughly solar [Si]/[Mg] requires the inclusion of more Si-enriched phases in the mantle. We performed preliminary tests by adding dry (Mg, Fe) 2 Si 2 O 6 to the upper mantle and increase the amount of (Mg, Fe)SiO 3 in the lower mantle accordingly in order to increase the ratio [Si]/[Mg] to roughly solar. Furthermore, we have used a lower Fe content in the Mg-silicates of ξ Fe ≈ 0.1, which is thought to approximately represent the average Fe content in the Earth's mantle (Sotin et al. 2007). Finally, to mimic the Earth as closely as possible we adopted temperature discontinuities in the mantle transition zone (MTZ) at the top of the upper mantle of ∆T MTZ = 300 K (Sotin et al. 2007). With this the maximum water content for an Earth analogue decreases by about a factor of 50 to ≈ 0.028 wt%, in good agreement with the value range obtained from Ohtani 2019 but below the suggested range from Marty 2012 and the equilibrium water content reported in Cowan et al. 2014. However, since we ignored the hydration of Pyroxene for this preliminary test, the actual predicted value for the maximum water content in the Earth's mantle from our model would be somewhat higher. Furthermore, we have not included an SiO 2 phase in the mantle composition, which could further increase the overall water storage capacity. In addition, the water content should be expected to be rather sensitive to ∆T MTZ and ∆T TBL as these parameters strongly affect the thermal profile in the mantle. These aspects indicate that the water content in the planetary interior is extremely sensitive to the bulk composition and the thermal state of the interior. Therefore, our model is suited for estimating the upper limits for the hydration effects for a given planetary mass and for Si/Mg = 2/3. Currently our model cannot be used to accurately predict the maximum water reservoirs in the interior of an object with a different composition. To this end, more diverse compositions need to be considered.
Second, in addition to excluding a number of mineral phases in the mantle that could play an import role for water storage, we Article number, page 19 of 34 A&A proofs: manuscript no. aanda also have not included possible effects on hydration of other elements s.a. C, Ni, Ca, S or Al. Various studies show, however, that water solubility can be rather sensitive to the incorporation of some of these elements. For example, Al has been found to affect the stability field of hydrous phases (Hermann et al. 2016) and could lead to the stability of hydrous phases up to higher pressures. It furthermore has been found to increase the water solubility in Olivine polymorphs and Pyroxene (Grant et al. 2006;Férot et al. 2012).
Accounting for more complex compositions as outlined above would also rise the question of how adequate our assumption of uniform composition in the mantle is (see secs. 2.3 and 2.5). For example the Al content in Olivine is found to decrease significantly with increasing pressure and decreasing temperature (Férot et al. 2012). As the Al content can strongly influence the hydration behaviour, its pressure and temperature dependence might be an important factor. Similar issues could arise for the Fe content but has not been further investigated by the others at this point. However, we argue that the generalization to more complex compositions could inevitably lead to the necessity of allowing for composition gradients at least in the mantles of the planets. This would require a more sophisticated geochemical model for the phase assemblages as functions of pressure and temperature and the corresponding partitioning of water between the coexisting phases. We plan to address some of these topics in a subsequent study. Third, our model is restricted to M ≤ 3 M ⊕ due to the simple model for the chemical equilibrium between the core and mantle. As a result, for higher masses the partial hydrogen pressure in the lower parts of the mantle can exceed 1 GPa where the Sieverts law underestimates the hydrogen solubility in iron (Sugimoto et al. 1992). In the context of exoplanets, however, the super-Earth regime 1 M/M ⊕ 10 is of great interest. Therefore, it is desirable to employ more sophisticated chemical models in the future so the model can be extended to higher masses.
Fourth, we treated ∆T TBL and ξ Fe as free parameters. Furthermore, the estimate of ∆T CMB via eq. 27 might not be universally applicable to planets of different sizes, compositions, and evolution histories. This renders these three aforementioned parameter as important sources of uncertainty. In Appendix A we present a brief sensitivity study of our model with respect to ξ Fe and ∆T CMB . Both parameters could have significant effects on the H and OH reservoirs and the resulting effects on the mean densities of the planets. In the super-Earth regime, the role of ξ Fe is difficult to estimate at this stage. This is because the effect of the iron content on the hydration behaviour of pPv could not be included in our model due to a lack of available data. Furthermore, the roles of pressure and temperature on the water content in pPv are, to our knowledge, not further constrained neither by experiments nor numerical simulations. However, the exact behaviour of X H 2 O (P, T, ξ Fe ) in pPv could significantly affect the predicted mantle and core reservoirs of super-Earths. We therefore emphasize the importance of further investigation of these minerals to allow more reliable predictions of internal H and OH reservoirs of terrestrial planets.

Summary and conclusions
We have presented planetary structure models that account for the hydration of the Mg-silicate mantles and the iron cores. This allows for assessing the corresponding effects of water-interior interactions and ocean differentiation on the mass-radius relation and the partitioning of water between surface and internal reservoirs.
We considered three types of planets with simple bulk compositions (see Fig. 1): Dry planets (Type-1), hydrous planets (Type-2) and fully differentiated ocean planets (Type-3). For a first set of models we only considered Olivine polymorphs, Perovskite, Magnesiowüstit and Brucite in the mantle with a uniform Fe content of 25 mol%. The core-to-mantle mass ratio was derived from the total magnesium number Mg#, which is a free parameter. We then computed mass-radius curves for different bulk compositions at 300 K surface temperature and 1 bar surface pressure in the mass range M = 0.1 − 3 M ⊕ and for different values for the temperature drop in the thermal boundary layer. Our main conclusions can be summarized as follows: 1. The storage capacity of H 2 O equivalent in Type-2 planets is ≈ 0 − 6 wt% for 0.2 ≤ Mg# ≤ 0.7 and M ≤ 3 M ⊕ . This corresponds to ocean depths in Type-3 planets of ≈ 0 − 800 km.

The effect of hydration on the planetary radius peaks at
3. The radii of super-Earths can be considerably smaller if water is stored in the interior instead of an isolated surface ocean.
5. The temperature drop at the TBL has a great influence on the H and OH storage capacity and the resulting total radii.
6. The results are very sensitive to the bulk composition and planet mass.
These results reveal the necessity of accounting for internal hydration for a more reliable characterization of the compositions and structures of exoplanets, in particular in the super-Earth regime. This will become even more essential with the increased precision of the mass-radius measurements that are expected in the near future. Here we have laid the foundations for further guiding the interpretation of observational data by incorporating these aspects into theoretical predictions of mass-radius relations.
Appendix A: The sensitivity of the results to ∆T CMB and ξ Fe In this section we investigate the sensitivity of our results to the assumed temperature drop at the CMB (∆T CMB ) and the iron content in the Mg-silicates in the mantle (ξ Fe ). We ran three comparison cases: ∆T CMB = 0.8∆T CMB,nom , ∆T CMB = 1.2∆T CMB,nom and ξ Fe = 0, where the subscript 'nom' refers to the nominal case as presented in the main text. For simplicity, we only ran the bracketing cases of Mg# = 0.2, 0.7 and ∆T TBL = 200, 1700 K. ξ Fe is kept at the nominal value when ∆T CMB was changed and vice versa. Furthermore, the convergence criteria for the planets were somewhat relaxed as we only aim to give a rough estimation of the effects of these parameters. Similarly to Figs. 7 and 8 of the main text, Figs. A.1 -A.6 show the comparisons between the Type-1 and Type-2 planets and the Type-2 and Type-3 planets. The thick and thin curves correspond to the nominal and comparison cases, respectively. We find that changing ∆T CMB by ±20 % has a rather minor effect on the overall results.
The differences between the comparison and nominal cases are most prominent at higher masses where the core reservoir dominates. This is expected since for fixed surface conditions the temperature drop at the CMB predominantly affects the temperature in the core and not in the mantle. Therefore, the mantle reservoirs are rather unaffected by ∆T CMB . Furthermore, an increase or decrease in ∆T CMB leads to an increase or decrease in the hydrogen content in the core, which is consistent with eq. 41.
The effects of ∆T CMB on δR 1 /R 1 are very weak while δR 2 /R 2 is shifted to slightly smaller/larger values for smaller/larger values of ∆T CMB for high masses and low Mg#. This indicates that ∆T CMB plays an important role for iron-rich super-Earths. The iron content in the mantle ξ Fe has a considerable effect on the main results, especially at lower masses, where the mantle reservoir dominates. This reflects the fact that the iron content in the mantle strongly influences the saturation water content in the upper mantle. If no iron is present in the mantle, the mantle reservoir is significantly decreased.
Furthermore, less iron in the mantle leads to larger cores at fixed bulk composition, which is reflected in slightly larger core reservoirs. In addition, the Pv-pPv transition is shifted to somewhat higher masses as the core to mantle ratio is enhanced for lower iron contents in the mantle. This is reflected in an offset of the core reservoirs between the nominal and the comparison case ξ Fe = 0. The effects on δR 1 /R 1 and δR 2 /R 2 are most prominent at lower masses and diminish at higher masses. It is still, however, unclear how reliably this behaviour reflects reality. This is because the water content in the lower mantle is intrinsically independent of ξ Fe in our model because no sophisticated model for the hydration of pPv could be employed. In reality, the water storage capacity in pPv could, however, very well be influenced by the iron content, which would render our results more sensitive to ξ Fe also at higher masses. These comparisons indicate that the iron content in the mantle is a crucial factor for determining internal H and OH reservoirs for low mass planets but its role is somewhat unclear for super-Earths at this stage.
The opposite is true for ∆T CMB , which mainly affects the core reservoirs of super-Earths. In the context of exoplanets most of the observed terrestrial planets are in the super-Earth regime. For these planets both ∆T CMB and ξ Fe could play an important role for determining the internal H and OH reservoirs and corresponding effects of hydration and ocean differentiation on massradius relations.

Appendix B: Thermal equation of state
Here we describe in detail the equations of state used for the different materials in the core, mantle, and surface ocean. All results acquired in this study have been obtained using pregenerated EoS tables to optimize computational performance. We use piece-wise logarithmic spacing for both the temperature and the pressure axis, that is, we divided each decade into 10 n grid points where we fixed n = 1 for all materials except for pure water where n = 2 was used. During the structure integration the corresponding parameters (s.a. ρ, K T or α th ) at given P and T were extracted from the tables by bi-linear interpolation using the four adjacent grid points in the P − T -plane. (Mg, Fe) 2 SiO 4 exhibits three major polymorphic structures denoted by α-Ol, β-Ol and γ-Ol and hereafter generally referred to as Olivine polymorphs. To model the phase transitions of these phases, we performed linear fits from a set of experimental data for the transition pressure as a function of temperature (see Table  E.4). The transition pressure is given as: The corresponding values for the fit coefficients a and b are given in Table B.5. These transitions depend on the Fe content in Olivine polymorphs and can be shifted by up to ≈ 0.5-2 GPa for Fe contents up to 25 mol %. The stability field of the γ-phase reaches up to at least ≈ 25 GPa and likely higher (Miyahara et al. 2011(Miyahara et al. , 2016 and depends on the temperature. For simplicity we omit the effect of the Fe content on the phase transitions and assume a temperature independent stability field for γ-Ol up to 30 GPa to ensure that our prescription yields a strict upper bound for the water content. In light of the other simplifying assumptions for the composition (Section 2.3) these approximations are justified. Furthermore, we point out that for super-Earths the mantles and hence the mantle reservoirs are dominated by the lower mantles. Therefore, the exact value for the dissociation pressure of Olivine becomes irrelevant for larger planetary masses.
Appendix B.1.2: Thermoelastic properties In order to compute thermoelastic properties of water saturated Olivine polymorphs for different Fe contents ξ Fe relevant for planetary interior modelling an adequate EoS is required. A large number of studies have investigated in depth the EoS parameters for the anhydrous case and some efforts have been made to prescribe hydration effects on these parameters (see references in Table E.1-E.3). The most commonly used analytical form of the EoS used in such studies is a third order Birch-Murnaghan equation. It relates the pressure P to the temperature T and the density ρ according to:  The reference bulk modulus and the reference density are functions of temperature and given by: where we have adopted linear correction terms to account for hydration. Here we neglect the effect of Fe content on a T , b T and c T . Furthermore we do not account for any effects of Fe and water content on K T,0 . The reference parameters are computed following Mao et al. 2016 for the density and bulk modulus: Here, ξ Fe and X H 2 O are to be inserted in mole and weight fraction, respectively. The parameter values for eq.(B.3), (B.4) and (B.5) used here are summarized in Table B.1 for the different phases along with the corresponding literature references. The coefficients for eq. (B.6) and (B.7) are listed in Table B.2.
All required parameters in eq. (B.2) can then be computed as functions of water and Fe content. The density, bulk modulus, thermal expansion coefficient, and water saturation content X H 2 O,sat are plotted for the relevant P − T range in  3 we discuss how the saturation water content X H 2 O,sat (P, T ) is computed.

Appendix B.2: Hydration model for Olivine polymorphs
The pressure in the compiled data set for water saturated Olivine polymorphs covers the entire stability range for each phase individually within the given temperature ranges. The latter is lim-ited to the relevant regime for the Earth's upper mantle (up to 2000 K). The colour maps in Fig. B.2 provide an overview of the density, bulk modulus, water content and thermal expansion for the Fe free case and for ξ Fe = 0.25. The overall effect of increasing temperature is a relatively strong decrease in the saturation water content X H 2 O,sat for all three phases. The water content reaches zero for T < 2500 K for all phases. This means that at higher temperatures the water content diminishes to zero   and extrapolation to higher temperatures should therefore be unproblematic with regard to the hydration model. At low pressures P 5 GPa and high temperatures T 1500 K the fit predicts an increase in water content with temperature. This behaviour might not be physical. We point out that the experimental data do not reach temperatures below ≈ 1000 K. This is because they have been performed in the context of the Earth's upper mantle for which temperatures between ≈ 1000 − 2000 K are relevant. Therefore, to model low values of ∆T TBL , it was inevitable to

Appendix B.3: Brucite
For Brucite (Br) we use experimental fits for the 3rd order Birch-Murnaghan EoS (eq. (B.2)-(B.5)). The correction terms ∂a T /∂X H 2 O , ∂b T /∂X H 2 O , and ∂c T /∂X H 2 O vanish in these cases. The EoS parameters used in this study along with the corresponding references are summarized in Table B.3.

Appendix B.4: Perovskite, Periclase, and Wüstite
For Perovskite (Pv), Periclase (Per), and Wüstite (Ws) we use experimental fits for the Mie-Grüneisen-Debye EoS. In this case the pressure as a function of density and temperature is given as: where: The Grüneisen parameter γ(T, P) at elevated temperatures and pressures is computed via eq. 12. The Debye temperature is given as: The EoS parameters used in this study along with the corresponding references are summarized in Table B.3. We model (Mg, Fe)O (Magnesiowüstite) in the lower mantle as a mixture of MgO and FeO using linear mixing (eq. 13). For (Mg, Fe)SiO 3 (Perovskite and post-Perovskite) we give the parameter values for the Fe free case in Table B.3. The effect of the Fe content on the thermoelastic properties is modelled by adjusting ρ 0 , K 0 , θ D,0 and γ 0 . For ρ 0 we taken the results for 3.8 mol% Fe and 6.1 mol% Fe from Sinmyo et al. 2014 Table  B.5). Townsend et al. 2015 have found stable configurations of hydrous pPv containing at up to 2.31 wt% of water but have not investigated the saturation content as a function of pressure and temperature. Here we assume a constant maximum water content in pPv of 3 wt%.
Appendix B.5: Pure Fe and FeH x For the iron in the cores of the modelled objects we proceed in analogy to Alibert 2014 and use the EoS fit to molecular dynamics calculations from Belonoshko 2010. It has a similar form than the 3rd order Birch-Murnaghan EoS but instead of the temperature dependance of K 0 , K 0 and ρ 0 an additional term is intro-duced to account for the thermal pressure: P(ρ, T ) =P th (ρ, T ) with P th (ρ, T ) = 3Rγ(T − T 0 )m/ρ. (B.14) Here, γ is the Grüneisen parameter, given by eq. (12), R the gas constant and m the molar mass of the modelled substance, here Fe. The EoS parameters are T 0 = 300 K, K 0 = 174 GPa, K 0 = 5.3, ρ 0 = 8334 kg/m 3 , m = m Fe = 0.0558 kg, γ 0 = 2.434 and q = 0.489. The effect of the hydrogen content in FeH x on the EoS parameters of pure Fe are taken from Thompson et al. 2018. For convenience we rewrite FeH x = ξ Fe Fe + (1 − ξ Fe )H where ξ Fe = 1 − ξ H and ξ H = x/(1 + x). We then express the change of density and bulk modulus as:   (French et al. 2015) include ices VII and X. Other high pressure ice phases are not considered at this stage. The solid-liquid phase boundary is computed following Wagner et al. 2002 andFrench et al. 2016 for the low and high pressure regimes, respectively.
To obtain the parameters d and k the change in density and bulk modulus between the pure Fe EoS and the values for FeH x from Table 2 in Thompson et al. 2018 have been fitted as linear functions of the pressure and the hydrogen content. From this we extracted the values d ≈ −1620 kg m −3 and k ≈ −40.2 GPa. With this simple model fit the densities in Table 2 of Thompson et al. 2018 for low values of x, which are relevant for this study, can be reproduced with relative deviations of 1.2 %. The hydrogen content in the core itself, x, is obtained from chemical equilibrium between the core and the Mg-silicate mantle and is discussed in more detail in Section 3.2.
Appendix B.5.1: Pure water To model the most relevant phases of pure water in the surface oceans we employ different EoS for different P − T regions. At pressures below ≈ 1 GPa, the IAPWS formulation for liquid water (Wagner et al. 2002) and water ice Ih (Feistel et al. 2006) allow for as very reliable prescriptions of all relevant thermodynamic parameters. At higher pressures, where less or no experimental data are available, we switch to results obtained from molecular dynamics simulations. For the solid phase we use the thermodynamic potentials from French et al. 2015. Their fit includes water ices VII and X, which are the most important high pressure water ice polymorphs. The transition from low pressure ice Ih to the high pressure ices occurs at ≈ 0.2 GPa regardless of the temperature. Hence we switch from Wagner et al. 2002to French et al. 2015 at 0.2 GPa. At high temperatures and pressures we switch to the results of Mazevet et al. 2018. Their fit was performed such that it is consistent with the IAPWS at low pressures P 1 GPa. We switch from the IAPWS to Mazevet at 0.2 GPa. The transition between solid and liquid phases at all pressures are given by the water solidus. The solidus line is taken from Wagner et al. 2002 up to 1 GPa. For P > 10 GPa we take 2). Shown are the Fe free case, that is ξ Fe = 0 (left) and ξ Fe = 0.25 (right). The dependence of the phase transitions on the Fe and water content are neglected here. The pressure and temperature range in the plot roughly marks the validity range of the model fit. We note that we have assumed a stability field of the γ-phase of up to 30 GPa, which overestimates the dissociation pressure of Ringwoodite. We imposed this assumption to ensure that our model yields strict upper bounds on the water content.
the results from French et al. 2016 and interpolate between 1-10 GPa. A summary of the phase diagram of pure water along with the EoS reference for each region used in this study is provided in Fig. B.1. For all phase regions the EoS is given in the form of a thermodynamic potential. This allows computing the adiabatic gradient directly from the EoS using: dT dP ad = Pα th c P ρ , (B.17) where ρ, c P , P and α th are the density, specific isobaric heat, pressure, and thermal expansion. The isobaric heat is computed