NUCLEAR SPECTROSCOPY OF IRREGULAR BODIES

46th Lunar and Planetary Science Conference (2015)
1501.pdf
NUCLEAR SPECTROSCOPY OF IRREGULAR BODIES: COMPARISON OF VESTA AND PHOBOS.
T. H. Prettyman1 and J. S. Hendricks2, 1Planetary Science Institute (1700 East Fort Lowell, Suite 106, Tucson, AZ
85719, [email protected]), 2TechSource, Inc. (1475 Central Ave., Suite 250, PO Box 988, Los Alamos, NM 87544,
[email protected]).
Introduction: Gamma rays and neutrons are steadily produced by the interaction of galactic cosmic rays
(GCR) within the surfaces of airless planetary bodies.
Gamma rays are also made by the decay of long-lived,
radioelements K, Th, and U found in regolith materials.
The elemental composition of planetary regoliths, to
depths of a few decimeters, can be determined by analyzing emitted radiation measured from orbit or in situ.
Nuclear spectroscopy is effective within within about 1
body radius of planetary surfaces.
Nuclear spectrometers have been deployed on missions to the Moon, terrestrial planets, and asteroids Eros
and Vesta [e.g. 1]. Since the Moon and planets are are
round, with minimal topographic relief, corrections for
measurement geometry are often based on the first-order assumption that the flux of radiation arriving at the
detector varies in proportion to the solid angle subtended by a sphere with the same radius as the planet.
The assumption of solid-angle proportionality (ASAP)
is strictly valid, regardless of the shape of the planetary
body, when the flux of emitted radiation is isotropic
(e.g. for uncollided gamma rays made by the decay of
radioelements). For “smooth,” spherical bodies, observed departures from solid-angle proportionality may
result from anisotropic emission distributions for GCR
secondary radiation and, in the case of neutrons, loss
due to beta decay [2].
Vesta is the smallest body for which elemental composition has been mapped using nuclear methods. Vesta
has a triaxial shape with high relief in comparison to the
Moon, Mercury and Mars. For orbital mapping studies
using data acquired by NASA Dawn Mission’s Gamma
Ray and Neutron Detector (GRaND) [3], the flux of
gamma rays and neutrons was assumed to vary in proportion to the solid angle of Vesta, calculated at each
measurement location from a shape model [3,4]. The
validity of ASAP was verified empirically by examining the variation of counts with solid angle at high latitudes, which were well sampled in Dawn’s low altitude
mapping orbit [4-8]. Variations in solid angle were removed from the measurements to reveal Vesta’s surface
composition. Since the correction for solid angle was
by far the largest made to GRaND counting data, validation of ASAP is important in assessing sources of systematic error in mapped data.
For GCR-induced radiation, the relationship between measured counts and surface composition may
not be fully explained by the solid angle. For example,
it is easy to imagine that secondary radiation production
A
B
Fig. 1 Shape models of A) Vesta with 261-km mean volumetric radius (MVR) and B) Phobos with 11-km MVR.
by GCRs would be suppressed in shielded regions (e.g.,
within basins, canyons and craters) or that radiation
emission could be enhanced by scattering with topographic features such as basin walls. If so, then removal
of solid angle might introduce artifacts of shape and topography that could obscure variations in surface chemical composition.
Summary: In this study, we quantify the errors introduced by ASAP using detailed simulations of GCR
interactions with irregular bodies. For these simulations,
we modified the Monte Carlo radiation transport code
MCNPX [9] to treat complex shapes. Preliminary results indicate that spatial artifacts introduced by ASAP
are minimal for low-altitude measurements by Dawn at
Vesta: Corrections based upon MCNPX calculations do
not significantly change the results of previous mapping
studies. In contrast, ASAP breaks down over large altitude ranges, and artifacts of topography are observed
close proximity to Phobos. Improved correction methods are needed for missions proposed to small solar system bodies and for the analysis of low-altitude any data
set for which the instrument field-of-view is comparable
to topographic relief.
Modeling: MCNPX 2.7.5 was modified to include
a new surface type: a shape model specified by a triangle-vertex list. A voxel-based traversal algorithm was
used to speed-up intersections with surfaces containing
a large number of triangles, enabling fast simulations of
complex shapes. Shapes of Phobos and Vesta, determined by stereophotoclinometry, are represented by interconnected quadrilaterals (ICQs) [10, 11]. For our
study, the ICQ mesh was decimated and converted to
triangle-vertex format for ingestion into MCNPX. The
decimated models contained about 5×104 triangles (Fig.
1). MCNPX was also modified to treat neutron decay,
which can affect the flux of neutrons arriving at high
altitudes.
46th Lunar and Planetary Science Conference (2015)
GCR protons and alpha particles with energy distributions given by [12] were modeled. The gamma ray
and neutron leakage currents (energy-angle distribution) were tallied on concentric spheres spanning a
range of orbital altitudes. The spheres were subdivided
spatially into circular longitude and latitude regions using cones. For comparison, the average solid angle was
calculated for each longitude-latitude region. The composition of both bodies was assumed to be spatially uniform, with Vesta modeled as H-free howardite and Phobos as H-rich CM chondrite.
Results and Discussion: The total neutron leakage
current calculated by MCNPX is plotted against solid
angle and equivalent altitude (heq) in Fig. 2A for Vesta
(black) and Phobos (red). The points cluster into five
groups for corresponding to five concentric spheres on
which current was tallied (radii indicated on chart).
For Vesta, the radii span Dawn’s low altitude mapping orbit. For Phobos, tally spheres were selected to
represent a range of altitudes that could be sampled by
a future mission [e.g. 13]. The leakage current for Phobos is suppressed due to the presence of H in CM, which
moderates neutrons to low energies where they are preferentially absorbed. In both cases, current is strongly
correlated with solid angle; however, residuals (top
panel, Fig. 2A) show departures in linearity over the
wide range of altitudes sampled for Phobos.
If ASAP was valid, then maps of the ratio of neutron
current to solid angle should be featureless when composition is spatially uniform. Maps of this ratio are
shown for the innermost tally spheres in Fig. 2B&C.
Artifacts of relief (e.g. associated with the giant impact
basin Rheasilvia) are absent in Vesta’s ratio map. The
4% min-max variation in the ratio is small compared to
that of GRaND measurements with solid angle removed
(e.g. ~20% dynamic range for epithermal neutron
counts [4]). In comparison, the ratio map shown for
Phobos has higher contrast (~20%), which would bias
compositional maps determed using ASAP. The lowest
ratios are found in and around Stickney crater, which
also corresponds to the lowest sampled altitudes.
Improved methods: A model-based approach that
treats deviations from ASAP (e.g. caused by anisotropy,
shielding, and secondary scattering) is needed for closeup measurements of small, irregular bodies. MCNPX
provides a framework for investigating the these effects
and can be used directly to correct for shape and topography. A combination of MCNPX and simplified models [3] is used to characterize each of the aforementioned effects. For example, topographic shielding of
GCRs is not expected to be a large effect: GCRs with
low incidence angle, which would be blocked by basin
walls, contribute relatively little to radiation production
compared to those with normal incidence.
1501.pdf
A
Representative
uncertainty
Vesta
R2=0.998
420 km
445 km
470 km
495 km
520 km
Radii of tally spheres
14.5 km
17.1 km
Phobos
R2=0.993
19.8 km
22.4 km
25 km
B. Vesta
Ratio
1.02
Rheasilvia
1.00
0.98
C. Phobos
Ratio
1.06
0.97
Stickney
0.88
-180E
0
180E
Fig. 2. See Results and Discussion (left).
Acknowledgements: This work was funded by the
Dawn Mission under contract with NASA/JPL and by
NASA/NIAC grant NNX14AT50G.
References: [1] Prettyman T. H. (2014), in Encyclopedia of the Solar System ISBN:9780124158450.
[2] Maurice S. et al. (2004) JGR, 109, E7. [3] Prettyman
T.H. et al. (2011), Space Sci. Rev., 163, 371-459.
[4] Prettyman T.H. et al. (2012) Science, 338(6104),
242-246. [5] Lawrence D.J. et al. (2013), MAPS, 22712288. [6] Prettyman et al. (2013), MAPS, 48, 22112236. [7] Peplowski P.N. et al. (2013), MAPS, 48, 22522270. [8] Yamashita et al. (2013), MAPS, 48, 22372251. [9] Pelowitz, D.B., Ed. (2011), LA-CP-11-00438.
[10] Gaskell, R.W. priv. comm. [11] Gaskell, R.W.,
Gaskell Phobos Shape Model V1.0. VO1-SAVISA/VISB-5-PHOBOSSHAPE-V1.0. NASA/PDS,
2011. [12] O’Neill P. M. (2010), IEEE TNS, 57(6),
3148-3152. [13] Raymond, C. A. et al. (2015), AstroRecon, Abstract #6029.