Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://selena.sai.msu.ru/Kozl/Publications/ASR%20Origin%20and%20Stability.pdf Äàòà èçìåíåíèÿ: Mon Apr 22 21:21:11 2013 Äàòà èíäåêñèðîâàíèÿ: Thu Feb 27 20:41:44 2014 Êîäèðîâêà: Windows-1251 Ïîèñêîâûå ñëîâà: astro-ph

Available online at www.sciencedirect.com

Advances in Space Research 50 (2012) 1638-1646 www.elsevier.com/locate/asr

Origin and stability of lunar polar volatiles
A.A. Berezhnoy , E.A. Kozlova, M.P. Sinitsyn, A.A. Shangaraev, V.V. Shevchenko
Sternberg Astronomical Institute, Moscow State University, Universitetskij pr., 13, Moscow, Russia Available online 24 March 2012

Abstract Temperature regime at the LCROSS impact site is studied. All detected species in the Cabeus crater as well as CH4 and CO clathrate hydrates except H2, CO, and CH4 are stable against evaporation at the LCROSS impact site. CO and CH4 can be chemisorbed at the surface of the regolith particles and exist in the form of clathrate hydrates in the lunar cold traps. Flux rates of delivery of volatile species by asteroids, micrometeoroids, O-rich, C-rich, and low-speed comets into the permanently shadowed regions are estimated. Significant amounts of H2O, CO, H2, H2S, SO2, and CO2 can be impact-produced during collisions between asteroids and O-rich comets with the Moon while CH3OH, NH3 and complex organic species survive during low-speed comet impacts as products of disequilibrium processes. C-rich comets are main sources of CH4, and C2H4. Ó 2012 COSPAR. Published by Elsevier Ltd. All rights reserved.
Keywords: Impact processes; Cold traps; Comets; Polar ices; Atmosphere; The Moon

1. Introduction The possible existence of water ice in permanently shadowed regions near the lunar poles was considered by Watson et al. (1961) and Arnold (1979). Ice-like radar echoes were detected there by the Clementine-Deep Space Network bistatic experiment (Nozette et al., 2001), but subsequent radar study of south pole of the Moon are consistent with the ice being present only as disseminated grains in the lunar regolith (Campbell et al., 2006). Based on Lunar Prospector neutron flux data Feldman et al. (2000) estimated the mass of water ice to be 2 Â 1011 kg in the south pole region, and its mass fraction as about 1.5% in the south polar caps. The possibility of polar volatiles other than water ice has been considered theoretically by Sprague et al. (1995) for
Corresponding author. Tel.: +7 095 939 1029; fax: +7 095 932 8841.

E-mail addresses: ber@sai.msu.ru (A.A. Berezhnoy), Hrulis@yandex. ru (E.A. Kozlova), mik2756@yandex.ru (M.P. Sinitsyn), art-shangaraev@ yandex.ru (A.A. Shangaraev), vladislav_shevch@mail.ru (V.V. Shevchenko).

the case of S-bearing species, Reed (1999) for the case of Hg, and Duxbury et al. (2001) for the case of clathrates. However, only active impact experiments and sample return missions from the lunar cold traps are able to give more detailed information about the chemical composition of lunar polar volatiles. During active Lunar Crater Observation and Sensing Satellite (LCROSS) impact experiment several volatile compounds including H-, C-, N-, S-bearing species were detected in the Cabeus crater near the south pole of the Moon (Colaprete et al., 2010; Gladstone et al., 2010). Properties of the LCROSS impact-produced cloud were studied by Hurley (2011). According to Lunar Reconnaissance Orbiter (LRO) Lunar Exploration Neutron Detector (LEND) data the LCROSS impact site inside the Cabeus crater demonstrates the highest hydrogen concentration in the lunar south polar region, corresponding to about 0.5-4 wt.% of water ice (Mitrofanov et al., 2010). In this paper we re-evaluate comet hypothesis of the origin of lunar polar volatiles (Shevchenko, 1999) based on new data from LCROSS impact experiment and consider the thermal stability of volatile compounds detected in the Cabeus crater.

0273-1177/$36.00 Ó 2012 COSPAR. Published by Elsevier Ltd. All rights reserved. http://dx.doi.org/10.1016/j.asr.2012.03.019


A.A. Berezhnoy et al. / Advances in Space Research 50 (2012) 1638-1646

1639

2. Temperature regime and thermal stability of volatiles at the south pole of the Moon 2.1. Description of used model To calculate the lighting conditions and temperature of the lunar surface we used data from Lunar Orbiter Laser Altimeter (LOLA) (PDS, 2012), working on board the LRO spacecraft. The data were taken with step of 0.05œ for latitude and with step of 0.5œ for longitude. The investigated part of the lunar surface was divided into areas and for each area we have determined on the basis of altimeter data the height, the slope angle and the orientation of the area with respect to other areas. To investigate the illumination regime for each area we determined the azimuths and the angular heights of all the surrounding areas in order to get real picture of the horizon. The temperature of any element of the surface was calculated in accordance with approach previously described by (Carruba and Corradini, 1999; Vasavada et al., 1999): T j ? ?ðF 3 þð1 À aÞÂ?F 2 þ F 4 þ F 5 Þ=e  r1=4 ; ð1Þ
T, K

40 35 30 25 20 15 10 5 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Part of th e day

Fig. 1. Temperature regime at the LCROSS impact site in the Cabeus crater.

where a is the albedo of the surface assumed to be 0.11 (Braden et al., 2011), e = 0.95 is the coefficient of radiation (Greenhagen and Paige, 2006). F2 is the flux of solar radiation which drops on the illuminated element of the surface, F3 is the lunar heat flux estimated as 0.016 W/m2 (Langseth et al., 1976), F4 is the reflected flux from the other illuminated elements of the surface calculated in accordance with (Ueno et al., 1991), and F5 is the flux of infrared radiation that drops on the element j calculated in accordance with (Ueno et al., 1991). The approach of (Ueno et al., 1991) does not take into account the thermal flux from the neighboring elements of the surface but we can neglect this flux due to low thermal conductivity of the lunar regolith, about 10À2 W mÀ1 KÀ1 (Langseth et al., 1976). We also neglect influence of the solar wind while solar wind particles are able to penetrate to permanently shadowed craters and change the chemical composition of lunar polar volatiles (Zimmerman et al., 2011). 2.2. Areas of thermal stability and masses of deposits of volatile species at the south pole of the Moon The spatial resolution of obtained temperature maps is about 1.5 km. Maximal temperatures in the Cabeus crater were estimated during "summer" lunar day, when illuminations are maximal. Permanently shadowed region is located at the north part of the Cabeus crater. Maximal day-time surface temperature at the LCROSS impact site is about 37 K, this value is in agreement with observations (Hayne et al., 2010), while the night-time surface temperature is about 20 K (see Fig. 1). Several volatile compounds such as CO, H2, Ca, Hg, Mg (Gladstone et al., 2010), H2O, H2S, NH3, SO2, C2H4, CO2, CH3OH, CH4, OH (Colaprete et al., 2010), and Na (Killen

et al., 2010) were detected during the LCROSS impact experiment in the Cabeus crater. It is adopted that species are stable on the surface against thermal evaporation if the evaporation rate is less than 10À10 cm/year (Zhang and Paige, 2009). The evaporation rates of surface volatile species as a function of temperature were calculated according to approach of (Schorghofer and Taylor, 2007). Volatility temperatures at 2 cm depth were calculated assuming that evaporation rates of volatile species at 2 cm depth is in 530 times less than at the surface in accordance with (Schorghofer and Taylor, 2007). Obtained results of estimations of surface and subsurface volatility temperatures and calculations of areas of thermal stability are presented at Table 1. Areas of thermal stability of volatile species are estimated also from Fig. 2 of (Paige et al., 2010). Subsurface temperature regime favors survival of subsurface polar volatiles; evaporation rate of subsurface volatiles is significantly lower than that of volatiles at the surface. Both factors lead to increasing of areas of thermal stability of subsurface ice deposits in comparison with that of surface volatiles (Siegler et al., 2011). Calculated volatility temperatures for pure surface volatiles are low limits because physical and chemical interactions with other volatiles and regolith are not considered. OH-containing minerals are stable against thermal vaporization even in the equatorial regions of the Moon at daytime temperatures of about 350 K (Pieters et al., 2009). Comparing temperature regime of the crater Cabeus with volatility temperatures of detected species shows that H2, CO, and CH4 are unstable against thermal vaporization. It is therefore reasonable to consider possibility of the existence of such volatile species as H2, CO, and CH4 in the Cabeus crater as not pure solids in the cold traps, but as chemisorbed at the surface of the regolith particles or in the form of clathrate hydrates. Masses of deposits of volatile species at the south pole of the Moon were estimated assuming that depth of such deposits is 1 m, areas of existence of these deposits are equal to areas of thermal stability of these species in the cold traps at 2 cm depth (Paige et al., 2010), and mass fraction of species in areas of its thermal stability is the same as


1640

A.A. Berezhnoy et al. / Advances in Space Research 50 (2012) 1638-1646 References to studies of temperature regime: 1 - regions above 80œ (Paige et al., 2010), 2 - this work for the case of vicinities of the crater Cabeus, 3 - regions above 80œ (Zhang and Paige, 2010). References to vapor pressure data: 3 - (Zhang and Paige, 2010), 4 - (Lunine and Stevenson, 1985), 5 - (Fray and Schmitt, 2009), 6-(Huber et al., 2009), 7 - (Ferreira and Lobo, 2011), 8 - (Schins et al., 1971), 9 - from (Moser et al., 1998), 10 - (Brutti et al., 2001).

Table 1 Volatility temperatures and areas of thermal stability of volatile species at the south pole of the Moon.

Area of regions of thermal stability at the surface (km2) (Diviner data)1

in the LCROSS impact site. It is assumed that evaporation rate of about 10À3 cm/s is required for sublimation of volatiles from LCROSS impact site; it corresponds to vapor pressure of about 10À5 bar. This value of vapor pressure achieves at 105 K for CO2, 215 K for H2O, 315 K for Hg, 380 K for S, 700 K for Mg, and 870 K for Ca while temperature of grains of the lunar regolith excavated by LCROSS impact to the exosphere and heated by direct solar light is about 300 K. Water ice evaporates from the LCROSS impact-produced cloud with mass of about 3150 kg (Colaprete et al., 2010). Mass fractions of CH4, C2H4, H2S, CO2, NH3, SO2, CH3OH, and H2O in the regolith were estimated from (Colaprete et al., 2010) and assuming that CO, CH4, C2H4, H2S, CO2, NH3, SO2, CH3OH, H2O, Hg, and S should be released from 50,000, 40,000, 35,000, 30,000, 25,000, 15,000, 12,000, 4000, 3150, 3000, and 2000 kg of regolith, respectively, based on relative volatility of these species. Mass fractions of Hg, CO and H2 were estimated from masses of H2 and CO in the LCROSS plume (Gladstone et al., 2011) while H2 should be released from 250,000 kg of regolith (Hurley et al., 2012). Let us note that we use value of mass of regolith for CO in 3.6 times less than (Hurley et al., 2012) because CO is not so extremely volatile in comparison with H2. Estimated masses of volatiles in the lunar cold traps are presented at Table 2. 3. Origin of lunar polar volatiles 3.1. Origin of metal atoms delivered to the lunar exosphere by the LCROSS impact Several elements such as Na, K, Ca, Mg already detected at the Cabeus crater have lunar origin because these elements are abundant on the Moon and at typical velocities of collisions between celestial bodies and the Moon the target-to-impactor mass ratio is high, about 10-50 (Cintala, 1992). Thus, only volatile elements such as H, C, N, and Hg which are almost absent on the Moon can be delivered to the cold traps by impactors. Na and K atoms can be released from the surface by UV solar photons and can migrate toward the lunar poles by its ballistic hops above the surface. Detection of Mg, Ca, and Hg in the LCROSS impact experiment can be explained that atoms of these elements like Na and K atoms can migrate toward the lunar poles. Atoms of Mg, Ca, and Hg may be released from the lunar surface in the equatorial regions by high-energetic processes such as micrometeoroid bombardment. The maximal temperature in the LCROSS impact site is about 1000 K (Hayne et al., 2010), it is high enough for thermal evaporation of Na, Mg, and Ca deposits (see Section 2.2). Vapor pressure above MgO and CaO reaches 10À5 bar at about 2000 and 2100 K, respectively, while evaporation occurs in the form of Ca, Mg, O2, and O (Jacobson, 1989). Vapor pressure above Na2O reaches 10À5 bar at about 1050 K while evaporation occurs in the form of Na, NaO2, and O2 (Steinberg and Schofield,

Area of regions of thermal stability at the surface (km2) (this work)2 Area of regions of thermal stability at 2 cm depth (km2) (model)1 Area of regions of thermal stability at the surface (km2) (model)1 Volatility temperature at 2 cm depth (K) Surface volatility temperature (K) Compound

CO CH4 COÁ5.75H2O C2H4 CH4Á7H2O H2S CO2 NH3 SO2 CH3OH H2O Hg S Na Mg Ca

16.85 223 33.43 405 41.83 47.85 53.45 635 70.55 905 100.85 1326 181 + 77 2268 3039 378 + 1210

18.95 253 36.53 44.65 46.53 53.75 59.55 70.45 78.35 1005 112.75 1486 202 + 87 2568 3429 426 + 1410

0 0 0 800 1300 3900 7800 18300 29000 59000 70000 90000 126000 190000 349000 349000

0 0 0 0 0 140 1400 6700 13000 37000 53000 90000 130000 190000 349000 349000

0 0 350 4500 6100 13500 21400 39000 53000 100000 133000 300000 349000 349000 349000 349000

0 0 500 900 1430 4600 8600 20100 32000 63000 72000 90000 130000 188000 349000 349000


A.A. Berezhnoy et al. / Advances in Space Research 50 (2012) 1638-1646

1641

1991).Temperature in the LCROSS impact-produced cloud is about 1000 K (Hayne et al., 2010), it is not sufficient for evaporation of Mg, Ca, and, possibly, Na from the silicates. Thus, Mg and Ca were delivered to the LCROSS impact-produced cloud and presented in the polar regolith mainly in the form of atoms while Na could be presented in Na-bearing silicates in the polar regolith. Based on constants of reactions between CaO, MgO and H2O, CO2 (NIST, 2012) and estimating temperature and pressure in the LCROSS impact-produced cloud from (Hurley, 2011) it was found that chemical reactions with participation of Mg-, and Ca- bearing species are too slow for significant changes of the chemical composition of Mg-, Ca- bearing species in the cloud. If Mg and Ca were released to the LCROSS impact-produced plume mainly in the form of molecules then these molecules can be quickly destroyed by solar photons because photolysis lifetimes of Mg- and Ca-bearing species are quite short, about 10 s; MgO and CaO photolysis leads to formation of photolysis-generated metal atoms with excess energies of about 0.5 eV (Berezhnoy, 2010). However, the quality of Mg and Ca observations (Hurley et al., 2012) is not high enough for careful estimation of velocities of Mg and Ca atoms and fraction of hot photolysis-generated Mg and Ca atoms in the LCROSS impact-produced cloud. 3.2. Sources of lunar polar volatiles Big amounts of volatiles are delivered to the Moon by low-speed short-period comets (Berezhnoy et al., 2003). However, volatile-rich asteroids are also important sources of lunar polar volatiles (Ong et al., 2010). Let us consider other types of impactors. Average velocity of collisions of interstellar comets with the Moon is about 70 km/s. For such high-speed impacts the amount of cometary material captured by the Moon is very low, about 10À3 (Ong et al., 2010). Thus, interstellar comets are just minor sources of lunar polar volatiles because the frequency of collisions of such comets is lower than that of short-period comets and the amount of cometary matter captured after impacts of high-speed interstellar comets is lower than that after impacts of low-speed short-period comets. To explain mass of hydrogen at the poles of the Moon (Feldman et al., 2000) about 100 impacts of 20 km/s mini-comets per year are required. This frequency of impacts is unrealistically high. Thus, mini-comets are just minor sources of lunar polar volatiles. Solar wind is considered as a main source of hydrogencontaining species such as H2, OH, and H2O at the poles of the Moon (Crider and Vondrak, 2002). However, the content of other volatile elements such as N, C, and S is very small as in the solar wind as in the equatorial lunar regolith. For example, in the solar wind S/H, N/H, and C/H ratios are equal to 2.6 Â 10À5, 4.2 Â 10À5, and 3.5 Â 10À4, respectively (von Steiger et al., 2000). Thus, solar wind cannot be considered as an important source of S-, N-, and C-containing species in the cold traps.

Outgassing of lunar interiors was studied by Fegley (1991), but outgassing rate is still poorly known. 3.3. Chemistry and physics of impact processes Let us consider the behavior of H-, C-, N-, and S- bearing species detected in the Cabeus crater during collisions between comets and volatile-rich asteroids with the Moon. Impact flux F of volatile species to the cold traps can be estimated as F ði; jÞ ? F
imp

ði; jÞF ret ði; jÞF

mass

ði; jÞF

cap

ði; jÞ;

ð 2Þ

where i is the number of considered mechanism of delivery, j is the number of studied compound, Fimp is the impact flux, Fret is the retaining fraction of the impact-produced cloud captured by the Moon, Fmass is the mass fraction of considered species in the impact-produced cloud after quenching of chemical reactions, Fcap is the capture probability of considered species. Impact flux and retaining mass fraction of volatile-rich asteroids to the Moon are assumed to be 1.5 Â 108 kg/year and 17%, respectively (Ong et al., 2010). Impact flux and retaining mass fraction of O-rich comets are assumed to be 1.2 Â 107 kg/year and 6.5%, respectively (Ong et al., 2010). Impact flux and retaining mass fraction of degassed C-rich comets are assumed to be the same as for O-rich short-period comets. The behavior of volatiles in the impact-produced cloud is determined by the elemental composition of impactors because the Moon is almost dry. The elemental composition of O-rich short-period comets is assumed to be that of the comet Halley (Delsemme, 1988). The elemental composition of volatile-rich asteroids is assumed to be that of CM chondrites (Lodders and Fegley, 1998). Inactive degassed C-rich comets can be modeled by substantial removal of water molecules from the comet Halley. However, during degassing of comets other volatile species such as CO, CO2, and CH4 evaporate. For this reason existence of comets with extremely low O/C ratio, less than 0.6 is still in doubt. The elemental composition of C-rich comets is assumed to be that of dust of comet Halley (Krueger, Kissel, 1987). The elemental composition of sporadic meteoroids is assumed to be that of CI chondrites (Lodders and Fegley, 1998). The quenching theory was already applied for study of chemical processes during collisions between comets and the Moon (Berezhnoi and Klumov, 1998; Berezhnoy et al., 2003). Namely, the initial temperatures and pressures in the impact-produced vapor cloud are so high that typical time scales of chemical reactions are shorter than the hydrodynamic time scale and the chemical composition of such a cloud is in equilibrium. During expansion of the cloud into vacuum temperatures and pressure are rapidly decrease while chemical time scales increase. Quenching of the chemical composition of impact-produced cloud occurs when chemical time scales are comparable with hydrodynamic time scale.


1642 Table 2 Origin of volatile species at the poles of the Moon. Content of species detected at the LCROSS impact-produced plume, in comets, and impact-produced cloud at the quenching of chemical reactions is given in units of number of molecules normalized to water content assumed to be equal to 100%. Compound Content in the LCROSS plume 6005 Mass fraction in polar regolith (wt.%) 0.05 Mass in the cold traps at the south pole (kg) 4 Â 108 Content in comets Product of photolysis 100 Product of photolysis 2-207 5-109 0.2-68 0.2-1.57 0.38 0.12-0.67 0.110 Product of photolysis 0.1-1.69 0.0212
12

Photolysis lifetime (s1) 6.7 Â 106

Main sources Asteroids, O-rich comets Asteroids Photolysis Comets, asteroids Asteroids Low-speed comets C-rich comets C-rich comets Asteroids Asteroids Meteoroids Low-speed comets Asteroids, meteoroids Comets

O-rich comets (kg/year2) 3 Â 104

Asteroids (kg/year2) 6 Â 104

C-rich comets (kg/year2) 7 Â 103

Low-speed comets (kg/ year) -

Micrometeoroids (kg/year3) A.A. Berezhnoy et al. / Advances in Space Research 50 (2012) 1638-1646 104

H4 2

H2O OH4 CO CO2 CH3OH CH
4

100 0.036 155 2.176 1.556 0.656 3.126 16.756 3.196 <0.65 6.036 Undetected 0.65

5.6 - 0.08 0.04 0.1 0.003 0.02 0.2 0.2 < 0.3 0.07 - 0.4

10 -

13

8.3 Â 104 5 Â 104 1.3 Â 106 5 Â 105 8.3 Â 104 1.3 Â 105 2.1 Â 104 3.1 Â 103 4.8 Â 103 106 5.6 Â 103 106 -

105 2 Â 10À3 2 Â 105 105 2 Â 10À7 10À 10À
2

2 Â 106 5 Â 10À2 2 Â 105 1.5 Â 106 2 Â 10À8 2 Â 10À4 10À
13

4 Â 10À3 10À10 2 Â 105 0.02 4 Â 10À8 2 Â 104 104 0.02 10À18 3 Â 10À7 0.01 2 Â 103 20

2 Â 10 - 2 Â 10 4 Â 10 4 Â 10 70 60 30 30 - 60 - 3

4

105 104 7 Â 103 7 Â 104 5 Â 10À7 2 Â 10À7 10À16 0.2 4 Â 103 4 Â 105 10À3 104 0.2

5 Â 108 2 Â 1010 2 Â 1011 2 Â 108 2 Â 109 4 Â 1010 1011 <2 Â 1012 5 Â 1010 - 10

3

3 2

C2H4 H2S SO2 S11 NH3 N2
4

9

5 Â 103 3 0.05 0.2 3 Â 104 20

105 2 Â 104 4 0.1 4 Â 104 3

Hg13

10À

6

References and remarks: 1 - for quiet Sun at 1 AU according to Huebner et al. (1992); 2 - quenching of the chemical composition of impact-produced cloud occurs at Tq = 1200 K and Pq = 0.03 bar; 3 - quenching of the chemical composition of impact-produced cloud occurs at Tq = 3000 K and Pq = 10 bar; 4 - impact flux to the whole Moon is given, capture of delivered species by the cold traps is not considered; 5 - quenching of the chemical composition of impact-produced cloud occurs at Tq = 3000 K and Pq = 10 bar; 3 - ratios given by Gladstone et al. (2011) are normalized to water content assumed to be equal 5.6 wt.% (Colaprete et al., 2010); 6 - (Colaprete et al., 2010); 7 - (Crovisier, 2007); 8 - (Altwegg et al., 1999), 9 - (Kim et al., 1997); 10 - (Ootsubo et al., 2010); 11 - photoionization lifetime is given, impact flux from meteoproids was calculated at assumption that all atoms of sulfur produced during photolysis of S-bearing species reach the cold traps; 12 - (Wyskoff et al., 1991), 13 - Hg content in comets is assumed to be 3 ppm.


A.A. Berezhnoy et al. / Advances in Space Research 50 (2012) 1638-1646

1643

Let us assume that quenching temperature Tq for the case of collisions between comets and volatile-rich asteroids with the Moon is the same as for collision between comet Shoemaker-Levy 9 with Jupiter, about 1200 K (Berezhnoi et al., 1996). For initial temperature and pressure in the impact-produced cloud, about 10,000 K and 10,000 bar (Berezhnoy, 2010), and c = 1.2 it corresponds to quenching pressure Pq equal to about 0.03 bar. Quenching parameters of chemical composition of clouds produced during collisions of meteoroids with the Moon are assumed to be 3000 K and 10 bar (Herzog et al., 2009). At the time of quenching of the composition of impact-produced cloud H2O, CO, H2, H2S, and CO2 molecules are abundant as for impacts of O-rich comets as for impacts of asteroids while SO2 is abundant only for impacts of volatile-rich asteroids (see Fig. 2). However, for impacts of O-rich comets and asteroids the content of other species such as CH4, C2H4, CH3OH, and NH3 is very low in comparison with the observed values at the LCROSS impactproduced plume (see Fig. 2). Other species such as CO, CH4, and C2H4 can be delivered to the atmosphere of the Moon by C-rich comets (see Fig. 3). For the case of low-speed impacts temperature in the impact-produced cloud is too low for full destruction of complex organic compounds contained in comets and at

1.E+06 1.E+04 1.E+02

[X]/[H2O], %
C2H
4

CO CH
4

CO NH
3

H
2

2

H2S 0.4 0.6 0.8 1

1.E+00 0 1.E-02 1.E-04 1.E-06 CH3OH OH 0.2

O(free)/C
1.2 SO
2

1.4

Fig. 3. Equilibrium content of species detected at the LCROSS impact site normalized to H2O content at 1200 K and 0.03 bar versus O(free)/C ratio in impacting comets. The O(free) value is the amount of oxygen atoms which are not connected with Mg, Si, Al, Fe, Ca, Na atoms. The elemental composition of comets is assumed to be that of the comet Halley (Delsemme, 1988). The O(free)/C ratio in the comet Halley is 1.3. Decreasing of O(free)/C ratio is performed by removing of H2O molecules from the elemental composition of the comet Halley.

a
1.E+02

[X]/[H2O], %

H

2 2

CO H2S

CO

1.E+00

SO

2

1.E-02 NH 1.E-04 800
3

CH

OH
4

1200

T, K

1600

2000

such impacts these species can survive. Based on study of survival of amino acids during low-speed 10 km comet impacts on Earth (Pierazzo and Chyba, 1999) it is assumed that fraction of surviving species of cometary origin is equal to 1, 0.1, 10À2, and 10À3 for the case of 7, 11, 15, and 20 km/s impacts, respectively. For this case about 5% and 15% of species originally presented in impactors are survived during impacts of near-Earth asteroids and short-period comets, respectively. Taking into account content of considered species in comets (see Table 2), velocity distribution of low-speed comet impacts (Jeffers et al., 2001), and its retaining factors (Ong et al., 2010) impact fluxes of species captured by the Moon and survived during low-speed comet impacts were estimated. Content of volatile species in asteroids is unknown. If content of considered species in asteroids is in 30 times lower than that in comets than comets are main sources of species for the case of low-speed impacts. 3.4. Capture of temporal impact-produced atmosphere in the cold traps

b
1.E+02

[X]/[H2O], %
CO SO
2 2

1.E+00

H

2

H2S CO

1.E-02

OH

1.E-04 800 1200 1600 2000

T, K
Fig. 2. Equilibrium content of species detected in the crater Cabeus during adiabatic cooling of impact-produced vapor cloud. Water vapor content is assumed to be equal to 100%. Initial temperature is 10,000 K, initial pressure is 10,000 bar, c = 1.2. The elemental composition of the vapor cloud is assumed to be that of the comet Halley for the case (a) (Delsemme, 1988), and mixture of equal masses of CI chondrites and ferroan anorthosites (Lodders and Fegley, 1998) for the case (b).

After an impact temporal atmosphere is formed over the whole Moon during about 2 Â 104 s (Berezhnoy et al., 2003). At typical height scale of temporary impact-produced atmosphere of about 100 km and cross sections of photolysis of considered species of about 10À17 cmÀ2 (Huebner et al., 1992) photolysis lifetimes of molecules decreases with increasing column density at column densities greater than 1010 cmÀ3 corresponding to masses of atmosphere greater than 109 kg or diameters of impactors greater than about 200 m. Without taking into account adsorption of solar photons by photolyzed species molecules at impact-produced lunar atmosphere are stable against photolysis at column densities greater than 1013 cmÀ3 (Berezhnoy et al., 2003). For the case of Martian atmosphere molecules are stable against photolysis at


1644

A.A. Berezhnoy et al. / Advances in Space Research 50 (2012) 1638-1646

column densities greater than 109 cmÀ3 (Fox and Dalgarno, 1979). Majority of volatiles are delivered by biggest impacting asteroids and comets with typical diameters of about 20-100 km (Ong et al., 2010), such impactors can produce temporal atmosphere with surface column density of about 1016 cmÀ3. For this reason it was assumed that dense impact-produced atmosphere protects species against photolysis. Ratio of typical time of capture of temporal atmosphere in cold traps td and lifetime of atmosphere ta (assuming that it is determined by photolysis and products of photolysis escape the Moon) can be estimated in accordance with (Berezhnoi and Klumov, 1998) as td =ta $ ðS m =S tr Þd c =tph v;
7 2

2 Â 10À3 for H2S, 5 Â 10À3 for C2H4, 7 Â 10À3 for NH3, 0.015 for CH4 and SO2, 0.03 for CO, 0.3 for H2O, 0.4 for CH3OH, 0.45 for CO2, and about unity for S and Hg. In fact capture probabilities increase with increasing mass of atmosphere impactors at the transition case between two ultimate regimes described by Eqs. (3) and (4). Fluxes of species of different origin to the cold traps estimated in accordance with adopted mass fluxes, retaining factors, chemical composition of the impact-produced cloud at the time of quenching of chemical reactions, and capture probabilities are given at Table 2. 4. Discussion and conclusion The LCROSS impact experiment gives us information about existence of volatile species in the Cabeus crater considered as one of the coldest lunar polar craters. Stability of these species against evaporation is studied; areas of thermal stability of studied species at the south pole of the Moon are estimated. Comparison of estimated masses of deposits of lunar polar volatiles and its flux rates shows that considered sources of majority of lunar polar volatiles can deliver required amounts of volatiles to the poles of the Moon for 104-106 year. It means that only small fraction of delivered volatiles survives during existence of the cold traps (about 109 year). Maximal required delivery time, about 108-109 year, is estimated for Hg, CH3OH, and NH3. Hg can be delivered to the poles from the upper layers of lunar regolith during meteoroid bombardment (Reed, 1999). Molecular nitrogen is the main nitrogen-bearing compound for almost all studied sources of volatiles. Formation of NH3 molecules from implanted protons and nitrogen atoms in the cold caps upon action of ionization radiation is possible; this mechanism is already discussed by Colaprete et al. (2010). Our estimates of flux rates of molecular nitrogen show that amounts of nitrogen of asteroid's and cometary origin are comparable while isotopic studies show that chondritic micrometeoroids are main sources of non-solar nitrogen on the Moon (Furi et al., Å 2012). Methanol can be formed also during reactions at the surface of cold grains of the polar regolith. Asteroids are main sources of H2, H2O, CO, CO2, and H2S while O-rich comets are responsible for delivery of H2 and CO. C-rich comets are main sources of CO, CH4, and C2H4. Low-speed comets deliver NH3, CH3OH, and complex organic species originally presented in comets while meteoroids deliver elemental sulfur to the poles of the Moon. Estimation of the isotopic composition of lunar polar ices will be useful for determination of the main sources of polar volatiles because each considered source has unique isotopic composition. Thermodynamic calculations of equilibrium chemical composition of impact-produced cloud at the time of quenching of chemical reactions, about 3000 K and 10 bar, formed after impact of micrometeoroids in cold traps shows that main impact-produced species are H2O,

ð3Þ

where Sm = 3.8 Â 10 km is the total area of the lunar surface, Str is the area of cold traps assumed to be equal to area of thermal stability of species at 2 cm depth (Paige et al., 2010), dc $ 1 km is the characteristic depth of permanently shadowed craters, tph is the photolysis lifetime taken from (Huebner et al., 1992), v $ 300 m/s is the thermal velocity of species at 200 K. Based on estimated td/ta ratio capture probabilities were estimated as unity for CO2, CH3OH, H2O, and S; about 0.7 for CO, CH4, and NH3; about 0.4 for SO2 and C2H4; and 0.2 for H2S. These estimations are upper limits of values of capture probabilities because photolysis of species below the exobase level, which may be important for 100 m - 1 km impactors, was not taken into account. During photolysis of impact-produced molecules in the temporal atmosphere preferable escape of atomic and molecular hydrogen into space occurs due to its low masses and high velocities, it leads to depletion of hydrogen content in species captured by lunar cold traps in comparison with that in the impact-produced cloud at the time of quenching. For the case of smaller impactors collisional atmosphere is not produced. Transition from collisionless to collisional atmosphere occurs at column densities of about 108 cmÀ3 and mass of atmosphere of about 5 Â 106 kg (Vondrak, 1974). During micrometeoroid bombardment temporary impact-produced collisional atmosphere is not produced and photolysis of species in the lunar exosphere must be taken into account. Based on measurements of terrestrial cosmic dust flux impact flux of sporadic meteoroids is assumed to be 2 Â 106 kg/year (Love and Brownlee, 1993). The retaining fraction Fret is estimated with usage of Maxwell velocity distribution function at 3000 K. For the case of micrometeoroid bombardment capture probabilities of species to the cold traps were estimated in accordance with (Berezhnoy et al., 2003) as F
cap

$ ctph S tr =t

mig S m

;

ð4Þ

where c $ 1 is the capture probability of molecules during their collisions with the surface of cold traps, tmig is the mean duration of hops. Using Eq. (5) capture probabilities to cold traps at the south pole of the Moon are estimated as


A.A. Berezhnoy et al. / Advances in Space Research 50 (2012) 1638-1646

1645

H2, CO, CO2, HSOH, SO, SO2, H2S, and N2 while other species such as NH3, C2H4, CH3OH, and CH4 can be destroyed during such impacts. For this reason preferable survival of NH3, C2H4, CH3OH, and CH4 is possible in subsurface layers better protected from meteoroid bombardment in comparison with surface. Different sources are responsible for delivery of simple molecules such as H2O, CO2, and SO2 and complex organic molecules to the poles of the Moon. Progress in our knowledge about sources of lunar polar volatiles can be achieved after delivery of samples of lunar polar deposits to Earth. Further study of thermal regime and remote sensing mapping of polar deposits are required for careful estimation of amount and spatial distribution of lunar polar volatiles. Acknowledgments This research is supported by RFBR Grant 11-0290440-Ukr_f_a. The authors wish to thank quest editor and three anonymous reviewers for useful comments and suggestions. References
Altwegg, K., Balsiger, H., Geiss, J. Composition of the volatile material in Halley's coma from in situ measurements. Space Sci. Rev. 90, 3-18, 1999. Arnold, J.R. Ice in the lunar polar regions. J. Geophys. Res. 84, 5659- 5668, 1979. Berezhnoi, A.A., Klumov, B.A., Fortov, V.E., et al. Collision of a comet with Jupiter: determination of fragment penetration depths from the molecular spectra. JETP Lett. 63, 405-410, 1996. Berezhnoi, A.A., Klumov, B.A. Lunar ice: can its origin be determined? JETP Lett. 68, 163-167, 1998. Berezhnoy, A.A., Hasebe, N., Hiramoto, T., et al. Possibility of the presence of S, SO2, and CO2 at the poles of the Moon. Publ. Astron. Soc. Jpn. 55, 859-870, 2003. Berezhnoy, A.A. Meteoroid bombardment as a source of the lunar exosphere. Adv. Space Res. 45, 70-76, 2010. Braden, S.E., Robinson, M.S., Denevi, B.W., et al. Reflectance of Mercury and the Moon. EPSC-DPS Joint Meeting 1737, 2011. Brutti, S., Ciccioli, A., Balducci, G., et al. Thermodynamic stabilities of intermediate phases in the Ca-Si system. J. Alloy. Compd., 317-318, 525-531, 2001. Campbell, D.B., Campbell, B.A., Carter, L.M., et al. No evidence for thick deposits of ice at the lunar south pole. Nature 443, 835-837, 2006. Carruba, V., Corradini, A. Lunar cold traps: effects of double shielding. Icarus 142, 402-413, 1999. Cintala, M. Impact-induced thermal effects in the lunar and Mercurian regoliths. J. Geophys. Res. 97, 947-973, 1992. Colaprete, A., Schultz, P., Heldmann, J., et al. Detection of water in the LCROSS ejecta plume. Science 330, 463-468, 2010. Crider, D.H., Vondrak, R.R. Hydrogen migration to the lunar poles by solar wind bombardment of the Moon. Adv. Space Res. 30, 1869- 1874, 2002. Crovisier, J. Cometary diversity and cometary families. Available from: arxiv:astro-ph/0703785, 2007. Delsemme, A.H. The chemistry of comets. Royal Soc. Philos. Trans., Ser. A. 325, 509-523, 1988. Duxbury, N.S., Nealson, K.H., Romanovsky, V.E. On the possibility of clathrate hydrates on the Moon. J. Geophys. Res. 106, 27811-27814, 2001.

Fegley, Br. Thermodynamic models of the chemistry of lunar volcanic gases. Geophys. Res. Lett. 18, 2073-2076, 1991. Feldman, W.C., Lawrence, D.J., Elphic, R.C., et al. Polar hydrogen deposits on the Moon. J. Geophys. Res. 105, 4175-4196, 2000. Ferreira, A.G.M., Lobo, L.Q. The low-pressure phase diagram of sulfur. J. Chem. Thermodyn. 43, 95-104, 2011. Fox, J.L., Dalgarno, A. Ionization, luminosity and heating of upper atmosphere of Mars. J. Geophys. Res. 84, 7315-7333, 1979. Fray, N., Schmitt, B. Sublimation of ices of astrophysical interest: a bibliographic review. Planet. Space Sci. 57, 2053-2080, 2009. Furi, E., Marty, B., Assonov, S.S. Constraints on the flux of meteoritic Å and cometary water on the Moon from volatile element (N-Ar) analyses of single lunar soil grains, Luna 24 core. Icarus 218, 220-229, 2012. Gladstone, G.R., Hurley, D.M., Retherford, K.D., et al. LRO-LAMP observations of the LCROSS impact plume. Science 330, 472-476, 2010. Gladstone, G.R., Hurley, D.M., Retherford, K.D., et al. Corrections and clarifications: LRO-LAMP observations of the LCROSS impact plume. Science 333, 1703, 2011. Greenhagen, B.T., Paige, D.A. Mapping lunar surface petrology using the mid-infrared emissivity maximum with the LRO diviner radiometer. Lunar Planet. Sci. Conf. 37, 2046, 2006. Hayne, P.O., Greenhagen, B.T., Foote, M.C., et al. Diviner lunar radiometer observations of the LCROSS impact. Science 330, 477- 479, 2010. ` Herzog, G.F., Moynier, F., Albarede, F., et al. Isotopic and elemental ? abundances of copper and zinc in lunar samples, Zagami, Pele's hairs, and a terrestrial basalt. Geochimica Cosmochimica Acta 73, 5884- 5904, 2009. Huber, M.L. Laesecke, A., Friend, D.G. The vapor pressure of mercury, NIST, NISTIR 6643, p. 59, 2009. Huebner, W.F., Keady, J.J., Lyon, S.P. Solar photo rates for planetary atmospheres and atmospheric pollutants. Astrophys. Space Sci. 195 (1-289), 291-294, 1992. Hurley, D.M. Modeling of the vapor release from the LCROSS impact: parametric dependencies. J. Geophys. Res. 116, E10007, http:// dx.doi.org/10.1029/2010JE003793, 2011. Hurley, D.M., Hibbitts, C., Ernst, C., et al. Modeling of the vapor release from the LCROSS impact: II. Observations from LAMP. J. Geophys. Res. 117, E00H07, http://dx.doi.org/10.1029/2011JE003841, 2012. Jacobson, N.S. Thermodynamic properties of some metal oxide-zirconia systems. NASA Technical Memorandum. NASA TM-102351 62, 1989. Jeffers, S.V., Manley, S.P., Bailey, M.E., et al. Near-Earth object velocity distributions and consequences for the Chicxulub impactor. Mon. Not. R. Astron. Soc. 327, 126-132, 2001. Killen, R.M., Potter, A.E., Hurley, D.M., et al. Observations of the lunar impact plume from the LCROSS event. Geophys. Res. Lett. 37, CiteID L23201, 2010. Kim, S.J., Bockelee-Morvan, D., Crovisier, J., et al. Abundances of SO and SO2 in comet Hale-Bopp (C/1995 O1). Earth Moon Planet. 78, 65- 66, 1997. Krueger, F.R., Kissel, J. The chemical composition of the dust of comet P/ Halley as measured by "PUMA" on board Vega-1. Naturwissenschaften 74, 312-316, 1987. Langseth, M.G., Keihm, S.J., Peters, K. Revised lunar heat-flow values, in: Proc. 7th Lunar Science Conference 3, 3143-3171, New York, Pergamon Press, 1976. Lodders, K., Fegley, B. The Planetary Scientist Companion. Oxford University Press, New York, 1998. Love, S.G., Brownlee, D.E. A direct measurement of the terrestrial mass accretion rate of cosmic dust. Science 262, 550-553, 1993. Lunine, J.I., Stevenson, D.J. Thermodynamics of clathrate hydrate at low and high pressures with application to the outer solar system. Astrophys. J. Suppl. Ser. 58, 493-531, 1985. Mitrofanov, I.G., Sanin, A.B., Boynton, W.V., et al. Hydrogen mapping of the lunar south pole using the LRO neutron detector experiment LEND. Science 330, 483-486, 2010.


1646

A.A. Berezhnoy et al. / Advances in Space Research 50 (2012) 1638-1646 Sprague, A.L., Hunten, D.M., Lodders, K. Sulfur at Mercury, elemental at the poles and sulfides in the regolith. Icarus 118, 211-215, 1995. Steinberg, M., Schofield, R. A reevaluation of the vaporization behavior of sodium oxide and the bond strengths of NaO and Na2O: implications for the mass spectrometric analyses of alkali/oxygen systems. J. Chem. Phys. 94, 3901-3907, 1991. Ueno, S., Mukai, T., Azuma H. Temperature distribution in the permanently shadowed area in lunar polar regions, in: Proc. 24th ISAS Lunar and Planetary Symposium. ISAS, Tokyo, pp. 190-196, 1991. Vasavada, A.R., Paige, D.A., Wood, S.E. Near-surface temperatures on Mercury and the Moon and the stability of polar ice deposits. Icarus 141, 179-193, 1999. Vondrak, R.R. Creation of an artificial lunar atmosphere. Nature 248, 657-659, 1974. von Steiger, R., Geiss, J., Fisk, L.A., et al. Composition of quasistationary solar wind flows from SWIGS/Ulysses. J. Geophys. Res. 105, 27217-27238, 2000. Watson, K., Murray, B.C., Brown, H. The behavior of volatiles on the lunar surface. J. Geophys. Res. 66, 3033-3045, 1961. Wyskoff, S., Tegler, S.C., Engel, L. Nitrogen abundance in comet Halley. Astrophys. J. 367, 641-648, 1991. Zhang, J.A., Paige, D.A. Cold-trapped organic compounds at the poles of the Moon and Mercury: implications for origins. Geophys. Res. Lett. 36, CiteID L16203, 2009. Zhang, J.A., Paige, D.A. Correction to "Cold-trapped organic compounds at the poles of the Moon and Mercury: implications for origins". Geophys. Res. Lett. 37, CiteID L03203, 2010. Zimmerman, M.I., Farrell, W.M., Stubbs, T.J., et al. Solar wind access to lunar polar craters: feedback between surface charging and plasma expansion. Geophys. Res. Lett. 38, CiteID L19202, 2011.

Moser, Z., Zakulski, W., Rzyman, K., et al. New thermodynamic data for liquid aluminum-magnesium alloys from emf, vapor pressures, and calorimetric studies. J. Phase Equilibria 19, 38-47, 1998. NIST. Available from: , 2012. Nozette, S., Spudis, P.D., Robinson, M.S., et al. Integration of lunar polar remote-sensing data sets: evidence for ice at the lunar south pole. J. Geophys. Res. 106, 23253-23266, 2001. Ong, L., Asphaug, E.I., Korycansky, D., et al. Volatile retention from cometary impacts on the Moon. Icarus 207, 578-589, 2010. Ootsubo, T., Usui, F., Kawakita, H., et al. Detection of parent H2O and CO2 molecules in the 2.5-5 lm spectrum of comet C/2007 N3 (Lulin) observed with AKARI.. Astrophys. J. Lett. 717, L66-L70, 2010. Paige, D.A., Siegler, M.A., Zhang, J.A., et al. Diviner lunar radiometer observations of cold traps in the Moon's south polar region. Science 330, 479-482, 2010. Pierazzo, E., Chyba, C.F. Amino acid survival in large cometary impacts. Meteorit. Planet. Sci. 34, 909-918, 1999. Pieters, C.M., Goswami, J.N., Clark, R.N., et al. Character and spatial distribution of OH/H2O on the surface of the Moon seen by M3 on Chandrayaan-1. Science 326, 568-572, 2009. PDS. Available from: , 2012. Reed Jr., G.W. Don't drink the water. Meteorit. Planet. Sci. 34, 809-812, 1999. Schins, H.E.J., Van Wijk, R.W.M., Dorpema, B. Heat-pipe boiling-point method and vapor pressure of 12 metallic elements in range 10-104 torr. Zeitschrift fur Metallkunde 62, 330-336, 1971. Schorghofer, N., Taylor, G.J. Subsurface migration of H2O at lunar cold traps. J. Geophys. Res. 112, E02010, 2007. Siegler, M.A., Bills, Br.G., Paige, D.A. Effects of orbital evolution on lunar ice stability. J. Geophys. Res., 116, CiteID E03010, 2011. Shevchenko, V.V. On the cometary origin of the lunar ice. Sol. Syst. Res. 33, 400-408, 1999.