Neutron spectroscopy is used to determine bulk water abundances in the near surface of planetary bodies. The Dynamic Albedo of Neutrons (DAN) instrument on the Mars Science Laboratory (MSL) rover, Curiosity, is able to determine the depth distribution of water and neutron absorbers in the top ~50 cm of the subsurface. In this dissertation, I focus on answering significant geologic questions by interpreting DAN results in the geologic context provided by other MSL and orbital datasets. This approach enabled me to investigate significant outstanding questions in Gale crater geology, with implications for the evolution and habitability of Mars.I mapped an extensive silicic volcaniclastic layer in the subsurface, the first identified and mapped on Mars. This layer served as a silica source for other silica-rich features. But unlike those features, this layer contains abundant rhyolitic glass, indicating an evolved volcanic origin. Similar material on Earth is produced by plate tectonics, so this layer has important implications for the evolution of Mars, which has no evidence of plate tectonics. One of the primary motivations for exploring Gale crater is a distinct clay mineral signature from orbital data of the Compact Reconnaissance Imaging Spectrometer at Mars (CRISM), which has also identified a corresponding hydration signature. I compared DAN and CRISM hydration results and found that CRISM hydration results are biased by the presence of regolith, indicating that this regolith is either more hydrated or has a different grain size texture than bedrock. Clay minerals are primary binding sites for organics on Earth, and most organic-mineral binding mechanisms involve either water or hydroxyl. This makes hydrated clays the most efficient hosts for organic preservation, but clays are normally dehydrated when measured by MSL. However, my DAN-derived water abundances are greater in the most clay-rich unit of Gale crater, suggesting that clay minerals may be hydrated in the subsurface. I developed a new amorphous component analysis method that simultaneously constrains clay mineral hydration and abundances of various hydrated amorphous phases. I found a strong correlation between “excess” water and smectites (expandable clay minerals), indicating that these clay minerals are hydrated in the subsurface.

Information about the elemental composition of a planetary surface can be determined using nuclear instrumentation such as gamma-ray and neutron spectrometers (GRNS). High-energy Galactic Cosmic Rays (GCRs) resulting from cosmic super novae isotropically bombard the surfaces of planetary bodies in space. When GCRs interact with a body’s surface, they can liberate neutrons in a process called spallation, resulting in neutrons and gamma rays being emitted from the planet’s surface; how GCRs and source particles (i.e. active neutron generators) interact with nearby nuclei defines the nuclear environment. In this work I describe the development of nuclear detection systems and techniques for future orbital and landed missions, as well as the implications of nuclear environments on a non-silicate (icy) planetary body. This work aids in the development of future NASA and international missions by presenting many of the capabilities and limitations of nuclear detection systems for a variety of planetary bodies (Earth, the Moon, metallic asteroids, icy moons). From bench top experiments to theoretical simulations, from geochemical hypotheses to instrument calibrations—nuclear planetary science is a challenging and rapidly expanding multidisciplinary field. In this work (1) I describe ground-truth verification of the neutron die-away method using a new type of elpasolite (Cs2YLiCl6:Ce) scintillator, (2) I explore the potential use of temporal neutron measurements on the surface of Titan through Monte-Carlo simulation models, and (3) I report on the experimental spatial efficiency and calibration details of the miniature neutron spectrometer (Mini-NS) on board the NASA LunaH-Map mission. This work presents a subset of planetary nuclear science and its many challenges in humanity's ongoing effort to explore strange new worlds.

Characterizing the surface mineralogy of asteroids is critical to constraining their formation history and provides insight into the processes of planetary formation. One method of determining the surface mineralogy of asteroids is comparison of their visible to near-infrared reflectance (VNIR) spectra with laboratory spectra from meteorites and minerals. Subsequent in-situ investigation of these asteroids by spacecraft can supplement or supersede interpretations derived from Earth-based observations.I investigated a suite of aubrites, sulfide minerals, and metal-rich chondrites in a variety of forms (hand samples, powders, and slabs) to identify similarities with ‘spectrally featureless’ asteroids. I collected VNIR spectra and powder X-ray diffraction patterns of these samples and compared their overall reflectance and spectral slope with X-complex and T-, L-, and D-type asteroid spectra. The Psyche Mission will orbit asteroid (16) Psyche beginning in 2026. I provide a pre-flight assessment of the surface composition of Psyche by comparing spectra of Psyche to a large spectral library of possible surface analog materials (e.g., iron meteorites, mesosiderites, pallasites, sulfides, enstatite, ordinary, and metal-rich chondrites, endmember silicates, and mixtures of silicates, metal, and sulfides). Spectra of Psyche are generally consistent with iron meteorite powder, mixtures of iron meteorite powder and low-Fe, low-Ca pyroxene, sulfide minerals, and the CH/CBb chondrite Isheyevo. Next, I demonstrate some anticipated capabilities of the Psyche Multispectral Imager by comparing spectral parameters derived from Imager-convolved data to those from high resolution laboratory spectra. I offer preliminary strategies for classifying surface composition based on Imager filter ratios and overall reflectance. Last, I present an assessment of a benchtop, commercial-off-the-shelf (COTS) version of the Psyche Imager. The COTS Imager uses the same model CCD and a similar f-number commercial camera lens. I measured the gain, full well, linearity, read noise, quantum efficiency, and modulation transfer function to compare with eventual calibration data from the flight Imager. I validate the results of a radiometric model developed for the flight Imager with signal measurements from the COTS Imager. This work demonstrates that the COTS Imager is an effective testbed for validating Imager requirements and developing software and procedures for eventual calibration of the flight instrument.

This work is concerned with how best to reconstruct images from limited angle tomographic measurements. An introduction to tomography and to limited angle tomography will be provided and a brief overview of the many fields to which this work may contribute is given.

The traditional tomographic image reconstruction approach involves Fourier domain representations. The classic Filtered Back Projection algorithm will be discussed and used for comparison throughout the work. Bayesian statistics and information entropy considerations will be described. The Maximum Entropy reconstruction method will be derived and its performance in limited angular measurement scenarios will be examined.

Many new approaches become available once the reconstruction problem is placed within an algebraic form of Ax=b in which the measurement geometry and instrument response are defined as the matrix A, the measured object as the column vector x, and the resulting measurements by b. It is straightforward to invert A. However, for the limited angle measurement scenarios of interest in this work, the inversion is highly underconstrained and has an infinite number of possible solutions x consistent with the measurements b in a high dimensional space.

The algebraic formulation leads to the need for high performing regularization approaches which add constraints based on prior information of what is being measured. These are constraints beyond the measurement matrix A added with the goal of selecting the best image from this vast uncertainty space. It is well established within this work that developing satisfactory regularization techniques is all but impossible except for the simplest pathological cases. There is a need to capture the "character" of the objects being measured.

The novel result of this effort will be in developing a reconstruction approach that will match whatever reconstruction approach has proven best for the types of objects being measured given full angular coverage. However, when confronted with limited angle tomographic situations or early in a series of measurements, the approach will rely on a prior understanding of the "character" of the objects measured. This understanding will be learned by a parallel Deep Neural Network from examples.

The lunar poles have hydrated materials in their permanently shadowed regions (PSRs), also known as lunar cold traps. These cold traps exist because of the Moon’s slight tilt of 1.5, which consequently creates these PSRs. In these shadows, the temperature remains cold enough to prevent the sublimation of volatile materials for timescales spanning that of geologic times [Hayne et. al 2015]. PSRs are significant because they create an environment where water ice can exist within the first meter of regolith at the lunar poles, where many cold traps are present. These volatile materials can be observed through a process called neutron spectroscopy. Neutron spectroscopy is a method of observing the neutron interactions caused by galactic and extragalactic cosmic ray proton collisions. Neutron interactions are more sensitive to hydrogen than other elements found in the regolith, and thus are a good indicator of hydrated materials. Using neutron spectroscopy, it is possible to detect the hydrogen in these cold traps up to a meter deep in the regolith, thus detecting the presence of hydrated materials, water, or ice.
For this study, we used the Monte Carlo Neutral Particle Transport Code (MCNP6) to create a homogenous sphere that represented the PSRs on Moon, and then modeled five differing water contents for the lunar regolith ranging from 0-20 percent weight. These percent weights were modeled after the estimates for Shackleton crater, data from Lunar Reconnaissance Orbiter (LRO) mission, and data from Lunar Orbiter Laser Altimeter (LOLA).
This study was created with the LunaH-Map mission as motivation, seeking to exhibit what neutron data might be observed. The LunaH-Map mission is an array of mini-Neutron Spectrometers that will orbit the Moon 8-20 km away from the lunar surface and map the spatial
distribution of hydrogen at the lunar poles. The plots generated show the relationship between neutron flux and energy from the surface of the Moon as well as from 10km away. This data provides insight into the benefits of collecting orbital data versus surface data, as well as illustrating what LunaH-Map might observe within a PSR.

Part I – I analyze a database of Smoothed Particle Hydrodynamics (SPH) simulations of collisions between planetary bodies and use the data to define semi-empirical models that reproduce remant masses. These models may be leveraged when detailed, time-dependent aspects of the collision are not paramount, but analytical intuition or a rapid solution is required, e.g. in ‘N-body simulations’. I find that the stratification of the planet is a non-negligible control on accretion efficiency. I also show that the absolute scale (total mass) of the collision may affect the accretion efficiency, with larger bodies more efficiently disrupting, as a function of gravitational binding energy. This is potentially due to impact velocities above the sound speed. The interplay of these dependencies implies that planet formation, depending on the dynamical environment, may be separated into stages marked by differentiation and the growth of planets more massive than the Moon.

Part II – I examine time-resolved neutron data from the Dynamic Albedo of Neutrons (DAN) instrument on the Mars Science Laboratory (MSL) Curiosity rover. I personally and independently developed a data analysis routine (described in the supplementary material in Chapter 2) that utilizes spectra from Monte Carlo N-Particle Transport models of the experiment and the Markov-chain Monte Carlo method to estimate bulk soil/rock properties. The method also identifies cross-correlation and degeneracies. I use data from two measurement campaigns that I targeted during remote operations at ASU. I find that alteration zones of a sandstone unit in Gale crater are markedly elevated in H content from the parent rock, consistent with the presence of amorphous silica. I posit that these deposits were formed by the most recent aqueous alteration events in the crater, since subsequent events would have produced matured forms of silica that were not observed. I also find that active dunes in Gale crater contain minimal water and I developed a Monte Carlo phase analysis routine to understand the amorphous materials in the dunes.

With the natural resources of earth depleting very fast, the natural resources of other celestial bodies are considered a potential replacement. Thus, there has been rise of space missions constantly and with it the need of more sophisticated spectrometer devices has increased. The most important requirement in such an application is low area and power consumption.

To save area, some scintillators have been developed that can resolve both neutrons and gamma events rather than traditional scintillators which can do only one of these and thus, the spacecraft needs two such devices. But with this development, the requirements out of the readout electronics has also increased which now need to discriminate between neutron and gamma events.

This work presents a novel architecture for discriminating such events and compares the results with another approach developed by a partner company. The results show excellent potential in this approach for the neutron-gamma discrimination and the team at ASU is going to expand on this design and build up a working prototype for the complete spectrometer device.