Документ взят из кэша поисковой машины. Адрес оригинального документа : http://selena.sai.msu.ru/Bus/Publications/ms38/ms063.pdf
Дата изменения: Mon Oct 6 18:09:59 2008
Дата индексирования: Thu Feb 27 20:57:37 2014
Кодировка:
Microsymposium 38, MS063, 2003 (A4 format)

FORMATION OF HYDRATED SILICATES IN EDGEWORTH-KUIPER BELT OBJECTS. A. B. Makalkin, Institute of Earth Physics, RAS, Moscow, RF (e-mail: makalkin@uipe-ras.scgis.ru); Dorofeeva, V. A. Vernadsky Institute of Geochemisry, (RAS), Moscow, RF (e-mail: dorofeeva@geokhi.ru); V. V. Busarev, Sternberg State Astronomical Institute, Moscow University, RF; (e-mail: busarev@sai.msu.ru). Introduction: Visible-range absorption bands at 600-750 nm are recently detected on two EdgeworthKuiper Belt (EKB) objects [1]. Most probably the spectral features may be attributed to hydrated silicates. In this paper we indicate possible conditions of the phyllosilicate formation in the EKB bodies. The model and results: The EKB objects orbit the Sun outwards of the Neptune orbit, 30 AU to 50 AU, and are possibly the most primitive large solid bodies in the solar system. According to [2,3], the EKB objects formed in situ though some part of their material could come with projectiles from the formation zones of giant planets, mainly of Neptune and Uranus. Contemporary models of the solar nebula [4,5] yield very low temperatures and pressures of T = 15-30 K and P = 10-9-10-10 bar at the radial distance (r) of 30-50 AU and the nebula age of about 106-107 yr, when the EKB bodies of sub-planetary size were formed. At these T-P conditions all volatiles excluding H2, He and Ne were in the solid state (mostly ices and some organics), and the abundance of rocky (silicate) dust component was lower than that of ices in accordance with the solar ratios of elements [6]. Phyllosilicates may form in EKB bodies during aqueous alteration [7]. A liquid state of water requires a considerable elevation of temperature in the bodies' interior. Plausible mechanisms of heating were decay of radionuclides (short-lived 26Al and long-lived 40K, 235 U, 238U and 232Th) in the bodies and mutual collisions between the bodies. In this paper we consider the former mechanism only. As shown in model calculations [8], the long-lived radioisotopes were insufficient for total melting of ice fraction in the icy satellites of giant planets of radii up to 800 km, although partial melting was possible. Was the concentration of captured 26Al (half-life 7.2 x 105 yr) sufficient for melting of water ice in the EKB objects? If the time of EKB bodies' formation was substantially larger than the half-life of 26Al, then independently of the isotope concentration it couldn't heat the EKB bodies with rather high efficiency. The formation time of hundreds-km-sized EKB bodies was from about one million year [9] to several tens of million years [10]. Taking into account the model of cometary body formation by [9], accretion of bodies up to 100 km in radius at the EKB distances 3550 AU within (1­1.5)в106 yr seems to be possible, though this time is near the lower limit of accretion timescales. In this consideration we suppose that formation of planetesimals at r of the EKB could begin several 105 yr after the collapse of the protosolar cloud. If the accretion of bodies of radius R = 100 km had completed no later than a few 26Al half-life times, the decay of this isotope provides enough heat to melt the water ice in the interiors of these bodies. To check this conclusion we adopt the mass fraction of silicates (silicates /(ices + dust) ) in EKB parent bodies of 0.30 [11]. The rock component with chondritic (solar) abundances of refractory elements contains 1.3 wt. % of aluminum. We also adopt the 26Al/27Al ratio of 1в10-5 which is obtained from the "canonical" initial 26 Al/27Al ratio of 5в10-5 and accretion time of a EKB body 1.6 Myr (after CAIs). This time possibly but not necessarily coincides with the age of the solar nebula (from the collapse stage). The above figures, giving the 26Al abundance, should be added with the decay energy of 26Al = 3 MeV per atom and its decay constant = 9.63в10-7 year-1 to yield the heat production rate Q = 0.4 J kg-1 yr-1. The time m required to heat a large EKB body to the water-ice melting point and to melt the ice in its interiors can be estimated from the equation





0

Q exp(-t )dt =



Tm

T0

c p dT + L f m w ,

(1)

where T0 = 30 K is the adopted value for the initial temperature of the body, Tm = 273 K is the melting temperature of water ice (a good approximation to at P< 25 MPa, characteristic for interiors of a EKB body of radius R 300 km), Lf = 3.34в105 J kg-1 is the latent heat of fusion for H2O, mw = 0.38 is adopted for the H2O mass fraction, cp is the thermal capacity at constant pressure per unit mass for the body's material. With some overestimation of cp at temperatures from 30 to 150 K we can take the temperature dependence of specific heat for all main components similar to that for water ice: cpw=7.67 T J kg-1 K-1 [12]. In this approximation the values of thermal capacities (all in J kg-1K-1) is: cp r=3.1 T for rocks (mainly silicates), cp CHON = 5.7 T for refractory organics, and cp vol =10 T for volatile organics and gases (the approximation for gases is most crude, but this has little effect due to their low content). We use also mass fraction of CHON mCHON = 0.22 and combined mass fraction of volatile organics and gases mvol+g = 0.10. With these values we obtain the thermal capacity for the mixture cp cp0 T , where cp0 = 6.1 J kg-1K-2. After substitution of this value in Equation (1) and integration we have the majoring estimation for the time m : m = --1 ln{1 - [c p 0 (Ti2 - T02 ) / 2 + L f mw ] / Q} 1.9 в 106 yr . (2) Thus the water ice in the bodies can be melted in less than 2 million years after the body formation and, hence, at the age of the solar nebula of 3.5 Myr. During this time only a surface layer of thickness R 10 km could be frozen, as follows from the simple estimation:
R ~ , (3) where is the thermal diffusivity related to the thermal conductivity k as = k /( c p ) . The temperature


dependence of for water ice I is = 0T -2 , where 0 9.1 в 10 -2 m 2 K 2 s -1 [13]. However, the porosity of ices p 0.5 decreases the thermal conductivity 5 to 50 times [14]. The porosity is the maximal at the surface and reduces to the low values at the bottom of the layer. Thus the reasonable estimate for the thermal diffusivity of the layer is ~ 10 -6 m 2 s -1 . The consequences of water ice heating are much more important. First, huge amount of water ice evaporated at low pressures in the porous medium should recondense in the upper layers of the bodies, substantially reducing their porosity. As a result of insulation of the interiors from the outer space, the pressure below the upper layer of thickness R would become higher than 1 bar and melting of water ice should occur when heated to T > Tm 270 K. Probable admixture of volatile organics might slightly decrease this temperature. Thus, as follows from Equations (1) and (2), internal water ocean in the young EKB bodies forms at their age 1.9 Myr, that is after 3.5 Myr after CAIs and solar nebula formation. Consider the evolution of the internal water ocean in a young EKB body of r = 100­300 km. The thermal convection in the ocean should be vigorous, if the Rayleigh number Ra is much higher than its critical value Racr 103. We can estimate the value -4 -1 Ra = g d 3 T /( ) , where 10 K is the volumetric thermal expansion coefficient of the mixture, dominated by liquid water, g is the gravitational acceleration ( g ~ 4 G R ), d 0.8­0.9 R is the convective layer thickness, T is the temperature difference across the layer, ~ 10 -7 m 2 s -1 and ~ 10 -5 m 2 s -1 are the thermal diffusivity and kinematic viscosity of the water­solids mixture. At d = 70 km and the very low value of T = 1 K we nevertheless obtain a very high value Ra 1021. The Nusselt number (Nu), which is the ratio of the total heat flow (including convective one) to the conductive flow is related to Ra by [15] Nu 0.2Ra 1 / 3 . With these data we can estimate the time of heat transport through the convective water ocean c by relation (3) where R and is substituted for d 0.8 R and the molecular thermal diffusivity is substituted for the effective thermal diffusivity e which accounts for convection, with e = Nu . At the above parameters we obtain c 103 years. The time is very short relative to the thermal evolution time scale of the bodies m 106 years. The latter is also the time scale for heat transport through the outer body's shell (thermal boundary layer) of thickness R. Owing to the rapid radial heat transport through the water ocean its T is stabilized near the temperature of maximum water density 277 K (the adiabatic compression for hundreds-km-sized bodies is negligible) and probably never exceeds 280 K. After a lapse of time a continuing decrease of radiogenic heat production yields the freezing of the internal ocean beginning (as in a usual terrestrial ocean) from the upper layer.

The lifetime of the water ocean (till the beginning of its freezing) in the early EKB body of radius R can be estimated by comparing the heat flow F1 underneath the solid shell (thermal boundary layer) of thickness R and the heat flow F2 in the shell. The flow F1 is generated in the interiors being heated by the 26Al decay and quickly transferred to the lithosphere. Thus we can write F1 1 3 ( R - R )Q exp(-t ) (4) 26 on the assumption of Al homogeneous distribution, where = 1.4 в 10 3 kg m -3 is the mean density of the body (calculated at the above fractions of components). Flow F2 can be written as F2 = k T / R , (5) where k is the thermal conductivity of shell, T = 273 ­30 240 K. We assume k = 2 W m -1K -1 , taking into account the empirical relation for crystalline ice k (T ) = 567 / T W m -1K -1 and compensating effect of increasing porosity from the base to the surface of the shell [16]. The freezing of the water ocean begins when the incoming flow from interiors F1 becomes lower than the flow F2 coming from the shell. By equating two flows from (4) and (5) we obtain the estimate of the lifetime of the ocean of liquid water as l 0.6­1.8 Myr for the bodies of radius R = 100­300 km respectively. This lifetime is quite sufficient for silicates to form phyllosilicates by reaction with water. This research was supported, in part, by the Russian Foundation for Basic Research Grant 01-0564318. References: [1] Boehnhardt, H., Delsanti, A., Hainaut, O. et al. (2002) Proceedings of Asteroids, Comets, Meteors (ACM 2002), 47-50, ESA-SP-500, Berlin, Germany. [2] Safronov, V. S. (1996) Solar Sys. Res., 30, 251-257. [3] Farinella, P., Davis, D.R., and Stern, S.A. (2000) In: Protostars and Planets IV (V. Mannings, A.P. Boss, and S.S. Russell, eds.), 12551282, Univ. of Arizona Press, Tucson, USA. [4] Makalkin, A. B., Dorofeeva, V. A. (1996) Solar Sys. Res., 30, 440-455. [5] Mousis, O., Gautier, D., BockeleeMorvan, D. et al. (2000) Icarus, 148, 513-525. [6] Lodders, K., Fegley, B. (1998) The planetary scientist's companion. Oxford Univ. Press. New York, Oxford, USA. [7] Prinn, R. G. and Fegley, B., Jr. (1989) In Origin and Evolution of Planetary and Satellite Atmospheres (Atreya , S. K. et al., eds.), p. 78-136, Univ. of Arizona, Tucson, USA. [8] Prialnik, D., Bar-Nun, A. (1990) Astrophys. J., 355, 281-286. [9] Weidenschilling, S. J. (1997) Icarus, 127, 290-306. [10] Kenyon, S. J. and Luu, J. X. (1998) Astron. J., 115, 2136-2160. [11] Busarev V. V., Dorofeeva, V. A., Makalkin, A. B. (2003) Earth, Moon, and Planets, 93 (in press). [12] Hobbs, P.V. (1974) Ice Physics. Clarendon, Oxford. [13] Kirk, R.L., Stevenson, D.J. (1987) Icarus, 69, 91-134. [14] Shoshany, Y., Prialnik, D., Podolack, M. (2002) Icarus, 157, 219-227. [15] Schubert, G., Cassen, P., and Young, R.E.: 1979, Icarus, 38, 192-211. [16] Spohn, T., Schubert, G.: 2003, Icarus, 161, 456-467.