Lunar Lava Tube Radiation Safety Analysis

G. De Angelis1,2,3, J.W.Wilson3, M.S.Clowdsley4, J.E.Nealy1, D.Humes3 and J.M.Clem5
1Old Dominion University, Norfolk, VA; 2Istituto Superiore di Sanita', Rome, Italy; 3NASA Langley Research Center, Hampton, VA; 4National Research Council, Washington, DC; 5Bartol Research Institute, Newark, DE

An analysis of radiation safety issues on lunar lava tubes as potential habitats has been performed. Lava tubes are basically formed when an active low viscosity lava flow develops a continuous and hard crust, which thickens and forms a roof above the still flowing lava stream, and at the end of the extrusion period an empty flow channel now free from molten magma is left [1]. Lava tubes are commonly observed on the Earth [2], on basaltic volcanic terrains, with typical sizes of the order of 1-2 km of length, and few meters for cross-sectional parameters (i.e. height and width). Under lunar conditions (lower gravity field, absence of atmosphere), lava channels and tubes are at least an order of magnitude larger in each size dimension [3], i.e. hundreds of meters wide by hundred of meters or more deep and tenths of kilometers long. For many years it has been [1-3] suggested that these natural cavities on the Moon could provide an ideal location for a manned lunar base, by providing shelter from various natural hazards, such as cosmic rays radiation, meteorites, micrometeorites impacts, and impact crater ejecta for example, and also providing a natural environmental control, with a nearly constant temperature of -20°C unlike that of the lunar surface showing extreme variation in its diurnal cycle. The purpose of this work is an assessment of the lunar lava tube radiation characteristics and an evaluation of the their actual safety features.

The formation of lava tubes is generally associated to the formation of "sinuous rilles" [4], frequently observed on the lunar basalt surface, especially in the maria floors, which formed from high extrusion and very low viscosity magma which filled the existing basins. In contrast to the so numerous flow channels in the form of sinuous rilles, real lava tubes cannot be easily observed on the Moon, for the reason of being subsurface objects, therefore unobservable in surface imagery, and only those with at least a partially collapsed roof are observable. Moreover, lunar surface imagery is at most at medium resolution, so rilles or tubes smaller than few meters wide are not observable with present lunar imagery. A catalog of lava tube candidates has been obtained by analyzing Lunar Orbiter and Apollo imagery along lunar rilles on the lunar nearside [5], and more than 90 candidates were identified in maria, namely Oceanus Procellarum, Mare Imbrium, Mare Serenitatis and Mare Tranquillitatis, as discontinuous rilles alternating between open lava channel segments and roofed-over segments. An estimation of the cross-sectional size of the observed lava tubes was performed by projecting the walls of the adjacent rille segments all along the roofed-over segments, whereas the length were measured directly from the imagery and the roof thickness was estimated through the craters superimposed to the uncollapsed roof. This catalog provided a large lunar lava tube data set, from which parameters typical for minimum, average and maximum values for lunar lava tube size have been extracted. The "minimum" values are such with respect to the currently available imagery, with tubes with a roof thickness of e.g. 3 m being currently unobservable.

The analysis has been performed considering radiation from galactic cosmic rays (GCR) interacting with the lunar surface. The surface has been modeled as a 5 m regolith layer, followed by rock. The regolith density profile has been obtained by combining data from groundbased radiophysical measurements and from in-situ analysis data from the Luna, Surveyor and Apollo missions, as modeled in [6], whereas for the rock layer a constant value of 3.3 g/cm3 has been used as typical of mare basalt rock. The composition has been adopted as the same for both surface and rock layers, and has been chosen as an average of the Apollo 12 surface samples [7, 8, 9], taken at the Oceanus Procellarum landing site, the region with the largest number of lava tube candidates in the catalog. Two different scenarios have been considered, namely a Lunar Night (Tsurface = 100 K) and a Lunar Day (Tsurface = 400 K) scenario, with temperature profiles for regolith and rock as from data from the Apollo 15 and Apollo 17 landing sites measurements [10, 11, 12]. The range of describing parameters provided by the existing database of lunar lava tubes has been incorporated into the transport calculation.

For the initial conditions a primary spectrum of GCR (p, a, HZE) for Solar Minimum conditions modulated at 510 MV has been adopted as background radiation, and a spectrum with particle fluxes equivalent to four times the intensity of the 29 September 1989 event has been adopted for Solar Particle Events (p). All primary particles heavier than protons have been approximated as individual nucleons, e.g. a He4 nuclei have been transported as 4 individual protons. Radiation profiles given by natural and induced radioactivity (a, b, g) have been taken into account. Particles (p, n, K+, K-, p+, p-, m+, m-, e-, e+, g) have been transported with the three-dimensional Monte Carlo transport code FLUKA [13]. The evaluation of the radiation safety-related quantities, namely the Effective Dose (E) and the Ambient Dose Equivalent (H*10), from particle fluence has been performed with the conversion coefficients from [14], whereas the physical quantity Absorbed Dose (D) has been obtained by rescaling with the use of the ICRP60 [15] radiation weighting factors wr. Both downward and upward (backscattering) particles have been considered in the dose evaluation, so the results obtained are to be taken as radiation exposure upper limits.

No significant differences in the results have been observed between the Lunar Night and the Lunar Day scenarios. After 6 m of depth no effects of radiation due to or induced by GCRs in both quiet and disturbed scenarios are observable in the simulation, and after far less than 1 m no effects of radiation due to or induced by SPE particles are observable. Natural and induced radioactivity seems not to play a significant role in the lava tube exposures. Mesons can interact in a quite dense material like lunar rocks before decay, and this is the reason why this usual component, mainly m+m-, is not present at larger depths. As a by-product of the transport results, the particle fluence from arriving GCR particles and from upward backscattering just at Moon surface and the relative dose equivalents have been obtained. Also in the very shallow and presently unobservable lunar lava tubes, with roof thickness of the order of 1-2 m, the doses are well below the montly, annual and career limits given by NCRP 132 [16]. The radiation safety of lunar lava tubes environments has been demonstrated.

References:

[01] Horz, F., 'Lava Tubes: Potential Shelters for Habitats', in Lunar Bases and Space Activities of the 21st Century (W.W.Mendell, Ed.), Lunar and Planetary Institute, Houston TX, p. 405, 1985.

[02] Cruikshank, D.P., and C.A.Woods, 'Lunar Rilles and Hawaiian Volcanic Features: Possible Analogues', Moon, Vol. 3, p. 412, 1972.

[03] Wilson, I., and J.W. Head., 'Ascent and Eruption of Basaltic Magma on the Earth and Moon', Journal of Geophysical Research, Vol. 86, p. 2971, 1981.

[04] Oberbeck, V.R., W.L. Quaide and R.G. Greeley, 'On the Origin of Lunar Sinuous Rilles', Modern Geology, Vol. 1, p. 75, 1969.

[05] Coombs, C.R., and B.R. Hawke, 'A Search for Intact Lava Tubes on the Moon: Possible Lunar Base Habitats', in The Second Conference on Lunar Bases and Space Activities of the 21stCentury (W.W.Mendell, Ed.), NASA CP-3166, Vol. I, p. 219, 1992.

[06] Cherkasov, I.I., and V.V. Shvarev, 'Lunar Soil Science' (in Russian), Nauka Publishers, Moscow, 1975.

[07] 'Apollo 12 Preliminary Science Report', NASA SP-235, 1970.

[08] Warner, J., 'Apollo 12 Lunar Sample Information', NASA TRR-353, 1970.

[09] Lodders, K., and B. Fegley, Jr., 'The Planetary Scientist's Companion', Oxford University Press, NY, 1998.

[10] 'Apollo 15 Preliminary Science Report', NASA MSC SP-289, 1972.

[11] 'Apollo 17 Preliminary Science Report', NASA JSC SP-330, 1972.

[12] Zharkov, V.N., 'Internal Structure of Earth and Planets' (in Russian), Nauka Publishers, Moscow, 1983.

[13] Fasso', A., A. Ferrari, A. Ranft, P. R. Sala, G.R. Stevenson and J.M.Zazula, `A comparison of FLUKA simulations with measurements of fluence and dose in calorimeter structures', Nuclear Instruments and Methods of the Physical Research, Sect. A, Vol. 332, p. 459, 1993.

[14] Pelliccioni, M., 'Protection Quantities and Conversion Coefficients for Use in Radiation Shielding', Journal of Nuclear Science and Technology, Suppl. 1, p. 854, 2000.

[15] ICRP, 'Recommendations of the International Commission on Radiological Protection', ICRP Publication N. 60, Annals of the ICRP, Vol. 21, Pergamon Press, Elmsford NY, 1991.

[16] NCRP, 'Radiation Protection Guidance for Activities in Low-Earth Orbit', National Council on Radiation Protection and Measurements, Publication N. 132, Bethesda MD, 2001.