Inefficient volatile loss from the Moon-forming disk: reconciling the giant impact hypothesis and a wet Moon
Miki Nakajima David J. Stevenson
































mnakajima@carnegiescience.edu

^1Department of Terrestrial Magnetism, Carnegie Institution for Science, 5241 Broad Branch Rd NW, Washington, DC 20015, USA. 

^2Division of Geological and Planetary Sciences, California Institute of Technology, 1200 E California Blvd, MC 150-21, Pasadena, CA 91125, USA.




The Earth's Moon is thought to have formed from a circumterrestrial disk generated by a giant impact between the proto-Earth and an impactor approximately 4.5 billion years ago.
Since this impact was energetic, the disk would have been hot (4000-6000 K) and partially vaporized (20-100 % by mass).
This formation process is thought to be responsible for the geochemical observation that the Moon is depleted in volatiles (e.g., K and Na). 
To explain this volatile depletion, some studies suggest the Moon-forming disk was rich in hydrogen, which was dissociated from water, and it escaped from the disk as a hydrodynamic wind accompanying heavier volatiles (hydrodynamic escape). 
This model predicts that the Moon should be significantly depleted in water, but this appears to contradict some of the recently measured lunar water abundances and D/H ratios that suggest that the Moon is more water-rich than previously thought.
Alternatively, the Moon could have retained its water if the upper parts (low pressure regions) of the disk were dominated by heavier species because hydrogen would have had to diffuse out from the heavy-element rich disk, and therefore the escape rate would have been limited by this slow diffusion process (diffusion-limited escape). 
To identify which escape the disk would have experienced and to quantify volatiles loss from the disk, we compute the thermal structure of the Moon-forming disk considering various bulk water abundances (100-1000 ppm) and mid-plane disk temperatures (2500-4000 K). Assuming that the disk consists of silicate (SiO_2 or Mg_2SiO_4) and water and that the disk is in the chemical equilibrium, our calculations show that the upper parts of the Moon-forming disk are dominated by heavy atoms or molecules (SiO and O at T_ mid > 2500 - 2800 K and H_2O at T_ mid < 2500 - 2800 K) and hydrogen is a minor species. This indicates that hydrogen escape would have been diffusion-limited, and therefore the amount of lost water and hydrogen would have been small compared to the initial abundance assumed. 
This result indicates that the giant impact hypothesis can be consistent with the water-rich Moon.
Furthermore, since the hydrogen wind would have been weak, the other volatiles would not have escaped either.
Thus, the observed volatile depletion of the Moon requires another mechanism. 


  





Moon; giant impact; volatiles; lunar water; volatile loss; hydrodynamic escape

https://doi.org/10.1016/j.epsl.2018.01.026














§ INTRODUCTION


It is widely accepted that the Earth's Moon formed by a collision between the proto-Earth and an impactor approximately 4.5 billion years ago <cit.>. This impact created a partially vaporized disk around the planet, from which the Moon accreted. 
In the standard version of this hypothesis, the impactor was approximately Mars-sized and the impact velocity was close to the escape velocity <cit.>.
This model has been favored because it can explain several observed aspects of the Earth-Moon system, such as its angular momentum, the lunar mass, and the small iron core of the Moon. 
However, the model cannot easily explain the observation that the Earth and Moon have identical or strikingly similar isotopic ratios (e.g., oxygen and tungsten, ) given that the impact simulations indicate that most of the disk materials originate from the impactor, which presumably had different isotopic ratios from the Earth.
It may be, however, possible that the impactor happened to have similar isotopic ratios to those of Earth because the inner solar system may have been well-mixed <cit.>. This idea might bolster the standard model, but it still requires an explanation for the identical tungsten isotopic ratios. 

Alternatively, <cit.> suggest that a small impactor hit a rapidly rotating Earth while <cit.> suggests that two half Earth-sized objects collided. 
In these models, the composition of the disk is similar to that of the Earth, and therefore, the isotopic similarities could be naturally explained. 
These new models are promising alternatives, but they may predict that the Earth's mantle becomes mixed by the energetic impact <cit.>. This appears to contradict observed anomalies of short-lived isotopes indicating that the Earth has never been completely mixed (e.g., ^182Hf-^182W,  and noble gases, ). A recent suggestion that the Moon could have formed as the outcome of merging of smaller Moons and multiple impacts <cit.> revives an old idea requiring some specific, perhaps unlikely dynamical conditions to be acceptable.

In addition to the isotopic ratios of the Earth-Moon system, the chemical compositions provide further essential information.
The giant impact has been thought to be at least partly responsible for the observation that the Moon is depleted in volatiles, such as K, Rb, Na, and other volatile elements <cit.>. 
Smaller K/Th and K/U ratios of the Moon than those of the Earth also indicate that the Moon is depleted in volatiles (K is more volatile than Th and U) <cit.>.
The Moon-forming disk was hot and partially vaporized (up to 4000-5000 K and 20-30% for the standard case and 6000-7000 K and 80-90% for the recent models, ). 
Hydrogen, which would have been dissociated from water at this high temperature, may have escaped from the hot Moon-forming disk as a wind (hydrodynamic escape) together with heavier atoms and molecules <cit.>. 


This model predicts that the Moon also lost a significant amount of water, but this appears to be inconsistent with some of the measured lunar water abundances.
Determining the bulk lunar water abundance is an active area of research; 
based on these measurements and modeling of the lunar interior evolution <cit.>, the bulk water content of the Moon has been estimated to range from <10 ppm <cit.> to a few hundred ppm <cit.>. 
<cit.> suggest that the crustal thickness estimated by GRAIL (34-43 km, ) can be consistent with an initially deep (∼ 700 km) lunar magma ocean with presence of water (270 - 1650 ppm).
On the other hand,  work using Cl and F in addition to H in apatite suggests that the water content of the Moon can be much lower <cit.>.
These results cover a wide range: the Moon could be “wet”, which indicates here that the Moon is as water-rich as Earth (Earth's bulk water abundance is estimated as a couple of hundred ppm, ) or could be drier (< 100 ppm).
Needless to say, the possibility that water is heterogeneously distributed within the Moon <cit.> makes it even more difficult to estimate the bulk lunar water abundance based on a small set of samples. 



If a significant amount of water escaped from the disk, the lunar D/H ratios should be more enhanced than that of the Earth because H is lighter and would have escaped more efficiently than D. However, analyses of pristine lunar water suggest that lunar D/H ratios may be comparable to the terrestrial values, which may indicate that water loss was insignificant  <cit.>.
It should be noted that measuring the bulk content of the indigenous water and D/H ratio is a very challenging task because the available lunar samples are limited and because a number of processes, including fractional crystallization, degassing, solar wind irradiation, and cosmic-ray spallation, would likely alter the original values.

Thus, the Moon is depleted in some volatiles, but it may or may not be depleted in water. 
To understand the history of the lunar volatiles, we propose to reevaluate the water loss mechanism. <cit.> suggest that hydrodynamic escape could have occurred and blew off the disk atmosphere when the disk temperature is 2000 K and the mean molecular weight of the disk is m̅=6 g mol^-1 (i.e., water in the disk was dissociated to 2H and O). Conventionally, hydrodynamic escape from a planetary atmosphere occurs when the Jeans parameter λ≡ GM_⊕m̅/RTr' is smaller than ∼ 2 <cit.> (the exact number of this criterion can vary depending on the atmospheric composition, , and the geometry, ).
Here, G is the gravitational constant, M_⊕ is the Earth mass, R is the gas constant, T is the temperature, and r' is the distance from the planet.
The work done by <cit.> is certainly insightful, but an important aspect here is that this criterion of λ has been developed for a gas that behaves as a material with a single molecular weight (for example the solar wind, which is primarily hydrogen), and it is necessary to understand if this model is applicable to the specific system of interest.
If the disk were dominated by heavier elements that were gravitationally bound (i.e., the escape fluxes of the heavy elements were negligible), the hydrodynamic escape model is no longer valid.
For hydrogen to escape from a disk dominated by heavier elements, it must diffuse out from the heavy elements that are gravitationally bound to the planet-disk system.
Thus, the hydrogen escape rate is limited by this diffusion process and this is called diffusion-limited escape, which is much slower than the “blow off” hydrodynamic escape. This type of hydrogen escape likely occurred from early planetary atmospheres <cit.>.



As an example, consider a disk that is dominated by a light element i and heavy element j.
We assume here that i and j are hydrogen and oxygen, respectively, and that their mole fractions are f_i=n_i/(n_i+n_j)=f_ H=0.1, where n_i and n_j are the number densities of the element i and j, respectively. 
The upper limit of the diffusion-limited escape rate is described as <cit.>,

    ϕ_l =  b_ij f_i (1/H_j - 1/H_i)∼ b_ij f_i/H_j,
 
where b_ij is the binary collision parameter between elements i and j, f_i is the mole fraction of the element i, H_i(=RT/m_i g, where m_i is the molecular weight of the element i) is the scale height of the element i. The subscripts i and j represent the elements i and j, respectively. R is the gas constant, and g is the gravity.  The last approximation is valid when H_i ≫ H_j. Here, we are assuming that the heavy element j is not escaping from the system.

Under the hard-sphere approximation, b_ij is described as <cit.>, 

    b_ij=3/64Q(2π RTm_i+m_j/m_i m_j)^1/2,
 
where

    Q=π/16(σ_i+σ_j)^2.
 
σ is the collision diameter.
Assuming i is atomic hydrogen and j is atomic oxygen, the hydrogen escape flux becomes 1.83 ×10^15 atoms m^-2s^-1 (σ_i=2 ×53 pm, σ_j=2 ×60 pm, m_i=1 g mol^-1, m_j=16 g mol^-1, T=2000 K, f_i=0.1, r=3 R_⊕, z=3 R_⊕, r'=√(r^2 + z^2)=4.2 R_⊕ where r is the horizontal distance from the planetary spin axis and z is the vertical distance from the disk mid-plane, and g=GM_⊕ z/r'^3. The choices of these parameters are discussed in Section <ref>). 
Assuming the surface area of the disk is 2π((5 R_⊕)^2 - R_⊕^2) and the disk life time is 1000 years, the total amount of lost hydrogen is 5.86× 10^14 kg and the equivalent amount of water is 5.86  ×10^14×(18/2) = 5.27 ×10^15 kg. It should be noted that this surface area is likely an upper limit because part the disk outside of the Roche radius would fragment.  
If the total mass of the disk is 1.5 lunar masses and the disk contains 100 ppm of water, then the mass fraction of the lost water with respect to the total water (i.e., the water loss mass fraction) is 4.78 × 10^-4. 
This is too small to have a significant effect on the interpretation of measurements of the water abundance or D/H ratio of the Moon.


Thus, determining the escape mechanism is highly important for estimating the volatile loss from the Moon-forming disk.
In this paper, we determine the structure of the disk and find that the upper parts of the disk are dominated by heavy atoms and molecules (Sections <ref> and <ref>).
This indicates that the hydrogen escape is likely diffusion-limited and that the amount of hydrogen and other volatiles escaping from the disk is too small to be observed.
This may indicate that the observed volatile loss (e.g., K, Na, and Rb) would require another explanation as discussed in Section <ref>.







§ MODEL


We assume that the disk has a liquid layer in the mid-plane that is sandwiched by vapor layers (this picture is similar to Figure 3 in , but we assume that the disk is isolated from the Earth's atmosphere).
We assume that the disk consists of water and silica (SiO_2 except Section <ref>, where Mg_2SiO_4 is considered). 
To estimate the hydrogen abundance in the upper parts of the disk, we simply investigate the vertical disk structure at a certain radial location (at r=3R_⊕ except Figure <ref>) given a mid-plane temperature T_ mid instead of modeling the whole disk structure. 
We make this simplification mainly because our outcome is insensitive to the detailed radial disk structure and partly because the time-dependent disk structure is not well known.
The effect of r on the water loss is considered in Section <ref>.




The surface density of the disk is assumed to be 5× 10^7 kg m^-2 based on previous work <cit.>.
We also test 10^8 kg m^-2 (not shown), but the outcome is similar to the case presented here. 
The model parameters are the bulk water content of the disk (100, 500, and 1000 ppm) and mid-plane temperature T_ mid (2500 - 4000 K).
This range can be thought of as representing different stages in the cooling of the disk or different giant impact models (or both).

To estimate the amount of water loss, we first determine the mixing ratio of silicate vapor and water vapor in the upper parts of the disk without considering dissociation of molecules (Section <ref>).
Based on the pressure and temperature ranges obtained, we estimate the mole fractions of the molecules and atoms considering dissociation given that the system in the chemical equilibrium (Section <ref>). If the disk is dominated by heavy elements, hydrogen escape is diffusion-limited. The amount of lost hydrogen and water is estimated using Equation (<ref>).





 §.§ Boundary condition at the liquid-vapor interface

At the liquid-vapor interface (at the mid-plane), the partial pressure of water is given as <cit.> 

    p_ H_2O=(y_ H_2O (wt%)/2.08× 10^-4)^1/0.54 (Pa),
 
where y_ H_2O is the mass fraction of water in the liquid.
The saturation vapor pressure of pure SiO_2 liquid (2000-6000 K) is written as (the units are modified from )

    p^*_ SiO_2=p_0 exp(-L/RT),
 
where p_0=1.596× 10^13 Pa and L=4.96× 10^5 J mol^-1.
The total pressure at the interface becomes p=(1-x^l_ H_2O)p^*_ SiO_2+p_ H_2O, where x^l_ H_2O is the mole fraction of water in the liquid (x^l_ H_2O=y_ H_2O/18 ( y_ H_2O/18+100-y_ H_2O/60)^-1).






 §.§ Vertical structure of the disk


We assume that the disk is in the radiative-convective equilibrium, which implies that lower part of the disk (i.e., small z) is convective and the upper part (i.e., large z) is radiative.  
In previous studies, the Moon-forming disk is assumed to be convective <cit.> or isothermal in the vertical direction <cit.>, but we modify this assumption in order to determine the temperature of the upper part of the disk. 
Provided the radial flow is highly subsonic, the disk is in hydrostatic equilibrium, i.e., dp/dz=-ρ g where p is the pressure, and ρ is the density. The disk self-gravity is ignored.



The convective region of the disk follows the moist pseudoadiabat curve, which assumes that the partial pressure of silicate vapor is equal to its saturation vapor pressure.
The amount of water extracted by solution in the silicate rain is small because of the modest water vapor pressure. This implies that only silicate rains out and this rain-out is efficient (if it were not then our conclusion that the water escape is minor would be even stronger since the upper atmosphere would contain droplets or particles that impede outflow). 
A moist pseudoadiabatic lapse rate is described as <cit.>,

    (∂ T/∂ p)_s= RT/pC_p, H_2O+x_s/x_ H_2OL/pC_p, H_2O/x_ H_2O+x_sC_p,s/C_p, H_2O+x_s/x_ H_2OL^2/RT^2 C_p, H_2O,
 
where x_s is the mole fraction of saturated silicate (=p_ SiO_2^*(T)/p) and x_ H_2O(=1-x_s) is the mole fraction of water. The subscripts, s and H_2O represent the parameters of SiO_2 and water, respectively.
It should be noted that dissociation of SiO_2 and H_2O are not considered in this section and it would be considered only in Section <ref>.
The effect of this simplification on the temperature profile is further discussed in Section <ref>).

C_p is the specific heat and C_p, s=62 JK^-1mol^-1, and C_p,  H_2O=55.7 JK^-1mol^-1 (at 1 bar at 3000K, ).


The optical depth throughout the disk is described as <cit.>

    dτ=(κ_s x_s m_s+κ_ H_2O x_ H_2O m_ H_2O)dp/m_ aveg,
 
where τ is the optical depth, κ is the absorption coefficient, m is the molecular weight, and m_ ave is the average molecular weight (m_s=m_ SiO_2=60 g mol^-1 and m_ H_2O=18 g mol^-1).
Here, κ_s=0.1 m^2kg^-1 <cit.> and κ_ H_2O=0.01 m^2kg^-1 are used <cit.>.
We assume that the disk is treated as a gray atmosphere and the upward radiation flux at given τ is written as 


    F_↑(τ)    =3/2∫_τ^τ_ midπ B(τ')exp[-3/2 (τ'-τ) ]dτ' 
        +π B(τ_ mid)exp[3/2 (τ-τ_ mid)],
 
where τ_ mid is the optical depth at the mid-plane (the liquid-vapor boundary).
Likewise, the downward radiation flux is written as


    F_↓(τ)=3/2∫_0^τπ B(τ')exp[-3/2 (τ'-τ) ]dτ' .
 
The net upward flux is written as F(τ)=F_↑(τ)-F_↓(τ). Here, π B=σ T^4 where σ is the Stefan-Boltzmann constant.
In the radiative part of the disk, the temperature profile follows the relationship

    π B=σ T^4(τ)=1/2F_↑ top(3/2τ+1),
 
where F_↑ top is the radiation flux from the top of the atmosphere (where τ=0).
The transition (tropopause) between the convective lower disk and radiative upper disk is iteratively determined.
First, assuming that the whole atmosphere is convective, we compute the vertical structure of the disk until the condition described by Equation (<ref>) is met.
This provides τ_ tp, T_ tp, and F_ tp where F_ tp=F_↑ top.
We recalculate Equations (<ref>) and (<ref>) and adjust the value of F_↑ top assuming that the lower parts of the disk are convective and the upper parts of the disk are radiative (generally speaking, this correction is very minor).





 §.§ Dissociation of molecules


At high temperatures and low pressures, molecules in the disk can dissociate. 
For example, assuming that the system is in the chemical equilibrium, the dissociation of SiO_2= SiO + O is described as

    K_ th=p_ SiO p_ O/p_ SiO_2=exp(-Δ G^0 /RT),
 
where K_ th is the equilibrium constant, and p_ SiO, p_ O, and p_ SiO_2 are the partial pressures of SiO, O and SiO_2 in bar. G^0 is the Gibbs free energy under the standard conditions (Δ G^0=Δ H^0 - TΔ S^0). 
Here, Δ H^0 is the change in the Helmholtz energy and Δ S^0 is the change in the entropy. 
We assume that Δ H^0 and Δ S^0 are not sensitive to temperature and pressure.
The rest of the reactions and the thermal constants are listed in Table <ref>.













 §.§ Homopause location


The escape flux of hydrogen is determined by thermal properties at the homopause, where eddy diffusion K equals molecular diffusion D (the molecular diffusion coefficient is described as D=b_ij/N where N is the number density). 
At lower parts of the disk (small z), eddy diffusion is more efficient than molecular diffusion (K>D), and therefore the disk is homogenized (i.e., the mole fractions of atoms and molecules are constant).
On the other hand, in the upper parts of the disk (large z), molecular diffusion becomes more dominant (K<D), and therefore each molecule or atom has its own scale height (i.e. light elements are more abundant at large z). 
If a hydrogen atom (or any light element) is present above the homopause, the atom can easily escape from the disk. 
However, this escape rate cannot exceed the hydrogen supply below the homopause. 
In other words, the supply rate is determined by how fast a hydrogen atom can diffuse from a heavy-element rich disk.
This is the definition of diffusion-limited escape as we briefly describe in Section <ref>.
We estimate that the homopause pressure range is ∼ 10^-4 - 10^1 Pa, as we discuss in Section <ref>. Even though the pressure range of homopause is not well defined, fortunately, our result is insensitive to this parameter.











§ RESULTS







 §.§ Vertical structure of the disk



The vertical structure of the disk is shown in Figure <ref>. The left panels show the temperature-pressure structure and right panels show the mixing ratio of water (x_ H_2O). 
The top, middle, and bottom panels correspond to the cases when the bulk water abundance is  (a) 100 ppm, (b) 500 ppm, and (c) 1000 ppm.
The location of the homopause is indicated by the shade (a large uncertainty, as discussed in Section <ref> and <ref>).
The lower parts of the disk  (i.e., large p and small z) are in the convection regime while the upper parts (i.e., small p and small z) are in the radiative regime.
The location of this transition (i.e., tropopause) is indicated by the kink in the temperature profile.
This kink appears because the disk temperature continues to decrease in the convective regime as z increases, whereas it is nearly constant in the radiative regime.



The water mixing ratio at the homopause is most sensitive to the mid-plane temperature.
As Figure <ref> a2-c2 shows, at high mid-plane temperatures (T_ mid > 2500 - 2800 K) the homopause is dominated by silicate vapor (x_ H_2O<0.5) because the partial pressure of the silicate vapor is larger than that of water vapor. In contrast, it is dominated by water vapor (x_ H_2O>0.5) at lower mid-plane temperatures (T_ mid < 2500 - 2800 K) because most of the silicate condenses into liquid and the partial pressure of silicate vapor becomes negligible.
As the bulk water abundance increases, the water mixing ratio increases, but this is a relatively weak effect compared to the mid-plane temperature.


When the disk is dominated by silicate vapor, the temperature of the disk closely follows the saturation vapor pressure of SiO_2 (i.e., it is thermodynamically determined rather than determined by the heat flow). 
This is a universal curve (independent of disk mass, heat flow, and water provided x_ H_2O≪ 1). 
In contrast, when the disk is dominated by water (x_ H_2O∼ 1), the disk temperature profile is determined by the heat flow (i.e., the radiative boundary condition) and the temperature can be less than the value it would have if it were determined by vapor pressure equilibrium with silicate because there is no longer significant silicate present. The homopause height is z ∼ 3 R_⊕ for a silicate-rich disk, whereas z ∼ 5 R_⊕ for a water-rich disk.




 §.§ Atoms and molecules present in the disk


The atoms and molecules present at the homopause in the disk at a constant temperature are shown in Figures <ref> and <ref>.
In Figure <ref> (a) and (b), the water mixing ratio x_ H_2O is 0.1 and 0.5, respectively while x_ H_2O is 1 in Figure <ref>. 
Figure <ref> corresponds to cases where the mid-plane temperature is high (i.e. T_ mid> 2500 - 2800 K) and the disk is dominated by silicate vapor, while Figure <ref>  corresponds to the cases where the mid-plane temperature is relatively small (i.e. T_ mid< 2500 - 2800 K) and the homopause is dominated by water vapor.

In Figure <ref> (a), the main reactions are described as SiO_2 = SiO + O and H_2O = 2H + O at p < 10^-3 Pa.
Hydrogen mole fraction f_ H at this pressure range is ∼0.1 (this is approximately estimated by  0.9 SiO_2 + 0.1 H_2O = 0.9 SiO + 1O + 0.2 H. f_ H∼ 0.2/(0.9+1+0.2)=0.095). 
A similar argument can be made for Figure <ref> (b), where 0.5 SiO_2 + 0.5 H_2O = 0.5 SiO + 1O + 1H, therefore f_ H = 1/(0.5+2)=0.4. 

This f_ H is likely to be an overestimate because the homopause temperature is below 2000 K at x_H_2O=0.5 (for example, the homopause temperature is 1660 K at T_ mid=2500 K and p=10^-2 Pa with 100 ppm of water as shown in a1, Figure <ref>). x_ H_2O>0.5 can occur when the water abundance is 1000 ppm (c2, Figure <ref> at T_ mid=2800 K), but 1000 ppm is likely to be larger than the actual bulk lunar water abundance (Section <ref>). 


When the homopause is dominated by water (x_ H_2O∼ 1), the homopause temperature is approximately 1600 K or below at p < 1 (Pa) (Figure <ref>). 
Under these circumstances, H_2O exists as its molecular form, and thus, hydrogen H is not abundant due to this relatively low temperature as shown in Figure <ref> (a). If the homopause temperature is close to 2000 K (Figure <ref>b), H_2O is dissociated to 2H and O, and hydrodynamic escape would be expected from such a disk as suggested by <cit.>. However, we suggest that this is unlikely to be the case for the lunar disk because the homopause temperature is smaller (∼ 1600 K). 






 §.§ Water loss from the disk


Hydrogen loss from the disk is estimated from Equation (<ref>).
First, we consider the case when the disk temperature is high and the disk is dominated by silicate vapor (f_ H < 0.4 as discussed in Section <ref>). 
This is likely the case shortly after the Moon-forming impact when the disk temperature is close to the initial disk temperature (4000-6000 K, ). 
Assuming the disk is dominated by O with a minor amount of H (f_ H=0.1, f_ O=0.9), the lost water by diffusion-limited escape is estimated as 5.27 ×10^15 kg and the mass fraction of lost lunar water is 4.78 × 10^-4 as discussed in Section <ref>. 

Now we consider the case when the homopause is dominated by water. Given the homopause temperature is 1600 K, f_ H=0.3, m_2=18 g mol^-1, σ_2 = 265 pm, z=5R_⊕ with bulk 100 ppm of water, the hydrogen escape flux is 1.64 × 10^15 atoms m^-2s^-1, the lost water mass is 4.73×10^15 kg and the mass fraction of lost water is 4.28× 10^-4. 

The mass fraction of water loss with various parameters (T=1600 K and 2000 K, f_ H=0.1, 0.3, and 0.5, r=3-7 R_⊕) is shown in Figure <ref>, which confirms that water loss from the disk is minor (< a few 10^-3 under these conditions). The range of r is taken from previous work (r≤ 5-7 R_⊕, ).













 
























§ DISCUSSION




 §.§ Mg_2SiO_4 disk


In Section <ref>, we only consider a disk that contains SiO_2 and H_2O, but the disk composition would be better modeled by Mg_2SiO_4. 
The disk temperature-pressure profile would not be significantly different given that the saturation vapor pressure of the bulk silicate Earth composition is similar to that of SiO_2 (Figure 3,  ), whereas species in the disk would be different. 
Figure <ref> shows the mole fractions of the species present in a disk of Mg_2SiO_4 and H_2O.
At x_ H_2O = 0.5, the main reactions are described as 0.5 Mg_2SiO_4=0.5 (2Mg + SiO + 3O) and 0.5 H_2O= 0.5 (2 H + O), which leads to f_ H=0.22. 
The presence of Mg decreases f_ H and makes it even more difficult to lose hydrogen from the disk. 
f_ H may decrease even further when more species (such as Al_2O_3 and CaO) that are expected to be in the disk are considered <cit.>. We do not consider effects of Mg on the solubility of water in the magma even though it could dissolve as Mg(OH)_2 (for detailed discussion, see ).
Furthermore, some OH can be present in the disk ( and personal communications with J. Melosh), which is not considered in this paper. It is possible that presence of OH can decrease f_ H and the amount of water loss.




 §.§ Escape of other volatiles


In the previous sections, we only consider the escape of hydrogen, but other volatile elements, such as K and Na, can be considered in a similar framework.
<cit.> consider the escape of a minor species with the presence of two major constituents.
The escape flux of a minor constituent k, ϕ_k, is approximately described as,


    ϕ_k=F_kϕ_ H(1-m_k-m_ H/m_2-m_ Hb_ Hk/b_ H2+m_2-m_k/m_2-m_ Hb_ Hk/b_ H2F_2/1+b_ Hk/b_2kF_2),
 

where the subscripts H, 2, and k represent hydrogen, heavy element that does not escape (corresponding to species such as SiO and O in the Moon-forming disk), and the element k, respectively.  
Here, F_k=n_k/n_H and F_2=n_2/n_H (the definition of F_k is slightly different from f_k=n_k/n_H+n_i+n_k as described in Section <ref>).


Assumptions of this model are (1) F_k  ≪ F_2, 1, and (2) F_2 is small. 
If the element k is potassium and the element 2 is atomic oxygen, given that σ_2=2× 60 pm, σ_k=2× 220 pm, m_2=16 g mol^-1, m_k=39.1 g mol^-1, and F_2=1, and the equation above becomes ϕ_k=0.0683 F_k ϕ_ H. Assuming that F_ k =  1, and ϕ_H∼ 2 × 10^15 atoms m^-2s^-1, ϕ_k=1.37 × 10^14 atoms m^-2s^-1, which corresponds to 1.73 ×10^15 kg. This is much smaller than the total potassium mass of the Moon is 2.64× 10^19 kg if the potassium abundance of the Moon-forming disk (1.5 lunar masses) is the same as terrestrial values  (240 ppm, ).
Equation (<ref>) may not necessarily be useful for estimating other volatiles such as Na, because Na may be a dominant constituent instead of a minor constituent <cit.>. However, it is still unlikely that this weak hydrogen flux would drag other heavier volatile elements.


Given that diffusion-limited escape is not an efficient mechanism for removing water or volatiles from the Moon-forming disk, the next question becomes how the Moon lost its volatiles or failed to efficiently accrete the volatiles.
A potential explanation is to accrete these volatiles onto Earth. 
<cit.> suggest that some volatiles that were initially present in the disk could have been preferentially accreted onto the Earth at the end of the Moon accretion process. 
The amount is not a significant part of the total volatile budget of Earth and may only have a small effect on the isotopic ratios for the Moon because high disk temperature may not cause significant isotopic fractionation especially if it is described as equilibrium isotope fractionation.
Slightly higher lunar potassium isotopic ratios than those of the Earth may also indicate that the Moon forms from liquid part of the disk <cit.>.
An alternative is to lose volatiles directly from the Moon as or after it accretes (Section <ref>). 






 §.§ Model validation


First, in Section <ref>, we assume that heavy elements, such as oxygen O, would not escape from the disk, but this is an approximation. Some oxygen would escape, but the amount is likely to be limited. The escape regime of oxygen would be diffusion-limited because O would be too heavy to hydrodynamically escape. Consider a disk that consists of SiO and O, and f_ O=0.5, T=2000 K, σ_1=2× 60 pm, σ_2=320 pm, z=3R_⊕ then the escape flux of O becomes 1.88× 10^15 molecules m^-2 s^-1. This is too small to affect the oxygen abundance in the Moon. 

For the two cases discussed in Section <ref> (silicate-rich and water-rich disks), the Jeans parameter λ=GM_⊕m̅/RTr' is estimated as 26.9 (m̅=30 g mol^-1, T=2000 K, and r'=4.2R_⊕) and 14.6 (m̅=18 g mol^-1, T=1600 K, and r'=5.8R_⊕), both of which are much smaller than 2, suggesting that hydrogen escape is diffusion-limited. It should be noted, however, that the criterion (i.e., hydrodynamic escape at λ <2) is developed for an isotropic atmosphere with a single component, and therefore it may not be directly applicable to a Moon-forming disk which is neither isotropic nor a single component. 
Nevertheless, λ being much larger than 2 further supports our model that the escape is diffusion-limited.


Furthermore, we use the formula of diffusion-limited escape developed for spherical geometry, but this geometry is not exactly applicable to the Moon-forming disk.
Nevertheless, we argue that the geometry becomes similar to a sphere given that the disk is extended in the vertical direction (z∼ 3 - 5 R_⊕).
Furthermore, the physics of diffusion-limited escape is not sensitive to the geometry.  
We have not explicitly modeled the non-thermal aspect of the disk evolution, such as the decrease in surface density as disk material spreads outward as well as disk materials accreting to the Moon and Earth simultaneously.
However, much of this effect is simply equivalent to changing surface density and the results are not sensitive to this. Of course, aggregating material at large radius or onto Earth can only decrease the escape. 
Moreover, here we assume that the disk is isolated, but if the interaction of the disk and atmosphere is considered, the escape may become even smaller because volatiles in the disk may be replenished from Earth.

In this paper, we focus on diffusion-limited escape, but volatile escape can be limited by other sources. 
Some of the previous studies investigate energy-limited escape, where the escape is caused by energy input from the stellar EUV, and therefore the rate is limited by the stellar flux.
In that case, the escape rate is written as ϕ_ EL=ε F^* r'/GMm^', where ε is the efficiency factor, m^' is the averaged molecular weight, F^* is the globally averaged EUV flux and Lyman α flux, 
which could have been  10^-2-10^-1 W m^-2 during this time period <cit.>. At ε = 0.3 <cit.>, m^'=12.9 g mol^-1 (mixture of 0.3 H and 0.7 H_2O), F^* = 0.1 Wm^-2, r'=4.2 R_⊕, this becomes 9.41 × 10^16 atoms m^-2s^-1 assuming the same disk surface area and disk life time discussed in Section <ref>. This value is much larger than the diffusion-limited case.
However, this is equivalent to losing 4.61 ×10^17 kg of water, which is still relatively small compared to the lunar bulk water (1.10 × 10^19 kg at 100 ppm).  
This energy-limited escape can be more efficient if the escape occurs at higher z, where f_ H is larger and m^' is smaller, than the homopause, but then the rate would be ultimately limited by the hydrogen supply from below, and therefore the escape rate may not exceed the diffusion-limited escape rate. Thus, it is possible that energy-limited escape would contribute to water loss, but the extent is likely to be limited.
We further discuss model assumptions in Section <ref>.



















 §.§ Comparison among different impact models

Our calculations show that hydrogen escape is minor at T_ mid=2500-4000 K, where 4000 K is the initial disk temperature estimated for the canonical Moon-forming impact.
T_ mid can be higher than 4000 K in the fast-spinning Earth and half-Earths models (T_ mid=6000-7000 K, ).
Nevertheless, even under this temperature rage, escape is still inefficient because f_ H at the homopause is small.
Thus, water escape would be inefficient even outside of the temperature range considered in this paper.






 §.§ Lunar volatiles after the Moon formation


In this paper, we focus on volatile loss during the disk phase, and it is possible that some volatiles were lost or added before and after this phase.
For example, the giant impact may have induced a vapor jet <cit.>, which may have removed some volatiles from the Moon-forming disk. 
This extent is difficult to quantify with the conventional impact method called smoothed particle hydrodynamics (SPH), where a fluid is expressed as a collection of spherical particles.
This is because the current SPH method cannot treat physics of two-phase (liquid-vapor) flow. 
Furthermore, some volatiles may have been added to or removed from the Moon before the lunar crust formation <cit.>, which depends on the impact flux and crust formation time scale (∼ 10 - 10^2 Myrs, ). 

We will further investigate evolution of lunar volatiles in the future.






§ CONCLUSIONS


We estimate the upper limit for hydrogen and volatile loss by thermal escape from the Moon-forming disk under various disk mid-plane temperatures (2500 - 4000 K) and the bulk water abundances (100, 500 and 1000 ppm). 
When the mid-plane disk temperature is large (> 2500-2800 K), the disk is dominated by silicate vapor. The major species in the upper part of the disk are SiO and O and the hydrogen mole fraction is small. In contrast, under low disk mid-plane temperature (< 2500-2800 K), the upper part of the disk is dominated by water and its temperature is  ∼ 1600 K or below. 
In this temperature range, water stays in its molecular form (H_2O) and the hydrogen mole fraction is small as well. Since hydrogen is not the major element and other heavy elements (such as O, SiO, SiO_2, and H_2O) are the dominant species, hydrogen would need to diffuse out from this heavy-element rich disk. This escape regime is called diffusion-limited escape and it is an inefficient escape process. We estimate the total mass of lost water and volatiles, such as potassium, and find that the escape is inefficient and that it would not remove water or other volatiles from the disk to the extent that the loss is measurable in lunar rock samples. 
To remove volatiles from the Moon-forming disk or from the Moon, another mechanism, such as losing volatiles from the disk to the Earth or degassing from the lunar surface, would be required.

 















§ ACKNOWLEDGEMENT


This work is supported by NASA Headquarters under the NASA Earth and Space Science Fellowship Program Grant NNX14AP26H and the Carnegie DTM Postdoctoral Fellowship. We would like to thank Masahiro Ikoma, Francis Nimmo, Cheng Li, Jay Melosh, Steve Desch, Erik Hauri, Alycia Weinberger, Peng Ni, and anonymous reviewers for helpful discussions.











§ SUPPLEMENTARY INFORMATION






 §.§ Homopause location


The location of the homopause is defined where K=D. The value of K is highly uncertain because it is determined by dynamical processes that are imperfectly understood, such as upward propagating waves that break or other fluid dynamical instabilities (for example, magnetorotational instabilities, MRI, in the Moon-forming disk have been previously discussed, ). 
An upper bound in the convective region follows from consideration of heat flow (essentially the bound provided by convective vigor). For a characteristic fluid velocity v_ conv of (F_ conv/ρ)^1/3∼ 10^2 m s^-1, where F_ conv is the convective flux and ρ is the density, and the scale height of the disk H̅ of hundreds of km, we expect K∼ v_ convH̅∼ 10^7 m^2 s^-1, but values in upper stably stratified regions can easily be many orders of magnitude lower. Simultaneously, the values can also be higher since the density is much lower (meaning that less energy per unit volume is required for mixing). 

The values of K for planetary atmospheres are typically 10^2 - 10^6 m^2 s^-1 <cit.> but these are for much less energetic systems and for smaller characteristic length scales, and therefore a higher value is possible in our case. Fortunately, the value of K does not matter much, as we now explain.


The upper limit of the escape flux for a light element i (hydrogen), ϕ_i is approximately described as ϕ_i ≤ϕ_l ∼  b_ij f_i/H_j (Equation <ref>) as discussed in Section <ref>.
Importantly, the parameters in this equation change rather little even as the pressure and number density at the homopause change by many orders of magnitude. It should be noted that Equation (<ref>) has a simple physical interpretation. Except for factors of order unity, it indicates that the escape flux is bounded above by ∼ nc, where n is the number density of the atom or molecule in question and c is the sound speed for that species, with the number density being evaluated at the place where the mean free path for the dominant species is of order H (this is the so-called exobase).


We can approximately estimate the homopause location based on the following two scenarios. If the region of the disk is dominated by oxygen atoms and has a small fraction of hydrogen atoms at 2000K, then D∼ K at p∼ 10^-3-10^1 Pa
(m_i=1 g mol^-1, m_j=16 g mol^-1, σ_i = 2 × 53 pm, and σ_j = 2×60 pm). 
If the disk is dominated by SiO at 2000K, this condition is met at p∼ 10^-4-10^0 Pa (m_i=1 g mol^-1, m_j=30 g mol^-1, σ_i = 2×53 pm, and σ_j = 320 pm). Thus, the homopause pressure range is estimated as ∼ 10^-4 - 10^1 Pa as discussed in Section <ref>.




 §.§ Further model validation


There are further assumptions in the model in addition to the points raised in Section <ref>.
One of the underlying assumptions of hydrodynamic escape is that once the flow reaches its sound speed as it expands outwards and eventually escapes from the disk and Earth.
According to previous studies <cit.>, this approximation is likely to be valid if the exobase is above the critical point. 
The critical point is where the velocity reaches the sound velocity.
This indicates that the criterion is l_c<H_1c, where l_c is the mean free path and H_1c is the scale height of a light element and the subscript c describes the critical point. 
This is rewritten as <cit.>


    ϕ_l≥GM_⊕/r'^2m_1/RT4b_11/√(π)
 

This becomes 5.69 × 10^15 atoms m^-2s^-1  (r'=4.24 R_⊕, T=2000 K,  m_1=1×10^-3 g mol^-1), which is comparable to the hydrogen escape flux we estimate.
It is possible that the escape flux becomes smaller than this value and thus the condition above may not be met.
In this case, the hydrogen flux is very weak and heavy elements would not be dragged to space <cit.>.
In other words, the diffusion-limited hydrogen escape rate is considered the upper limit and this makes our argument even stronger that hydrogen and volatile loss from the disk are minor.

Another assumption is that we ignore effects of dissociation on the lapse rate for simplicity (Section <ref>). The reaction of SiO_2= SiO + O is endothermic and the absolute value of the enthalpy is comparable to that of SiO_2 condensation. This indicates that this reaction would affect the lapse rate, but this would still not change the final outcome because our result that the water escape is not efficient holds under the wide temperature range (2500-4000 K).










































































model2-names



























