Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data inference”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Inference of the ring current ion composition by means of charge exchange decay

The analysis of the measured ion fluxes during the several day storm recovery period and the assumption that beside hydrogen other ions were present and that the decays were exponential in nature, it was possible to establish three separate lifetimes for the ions. These fitted decay lifetimes are in excellent agreement with the expected charge exchange decay lifetimes for H(+), O(+), and He(+) in the energy and L-value range of the data. This inference technique, thus, establishes the presence of measurable and appreciable quantities of oxygen and helium ions as well as protons in the storm-time ring current. Indications that He(+) may also be present under these same conditions were found.

Smith, P. H.↗

Describing hadronization via histories and observables for Monte-Carlo event reweighting

We introduce a novel method for extracting a fragmentation model directly from experimental data without requiring an explicit parametric form, called Histories and Observables for Monte-Carlo Event Reweighting (HOMER), consisting of three steps: the training of a classifier between simulation and data, the inference of single fragmentation weights, and the calculation of the weight for the full hadronization chain. We illustrate the use of HOMER on a simplified hadronization problem, a q\bar{q} q q ‾ string fragmenting into pions, and extract a modified Lund string fragmentation function f(z) f ( z ) . We then demonstrate the use of HOMER on three types of experimental data: (i) binned distributions of high-level observables, (ii) unbinned event-by-event distributions of these observables, and (iii) full particle cloud information. After demonstrating that f(z) f ( z ) can be extracted from data (the inverse of hadronization), we also show that, at least in this limited setup, the fidelity of the extracted f(z) f ( z ) suffers only limited loss when moving from (i) to (ii) to (iii). Public code is available at https://gitlab.com/uchep/mlhad.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Uncertainty quantification for high explosive reactant and product equations of state

Equations of state (EOSs) are typically represented as physics-informed models with tunable parameters that are adjusted to replicate calibration data as closely as possible. Uncertainty quantification (UQ) allows for the development of an ensemble of EOS parameters that are consistent with the calibration data instead of a single EOS. In this work, we perform UQ for the reactant and product EOSs for a variety of high explosives (HEs). In doing so, we demonstrate a strategy for dealing with heterogeneous (both experimental and calculated) data. We also use a statistical distance metric to quantify the differences between the various HEs using the UQ results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Preliminary Results on Bayesian Inverse UQ for OECD/NEA WPNCS Subgroup 14 Benchmark Exercise for Error Recovery and Experimental Coverage

The Organization for Economic Cooperation and Development (OECD) Working Party on Nucelar Criticality Safety (WPNCS) has proposed a benchmark exercise representative of neutronic behavior in criticality experiments. Here, the goal is to develop confidence in data assimilation techniques used to adjust nuclear data. Participants are given synthetic experimental models with associated measured data and asked to estimate the model parameters given the model and measurements as well as provide predictions for separate application models. In this work, we performed data assimilation using Bayesian inverse Uncertainty Quantification (UQ) with machine learning surrogate models to produce posterior parameter distributions for the requested parameters and posterior predictive distributions for the requested responses. Several experimental models are shown to insufficiently inform the posterior parameter distributions for the applications involved. However, given sufficient experimental data, posterior parameter estimates yielded reduced uncertainty in the response predictions of interest while covering the experimental data.

Bayesian Inference↗

Neutron Leakage Spectra Sensitivities for ICSBEP Benchmarks

Neutron leakage spectra have been measured, simulated, and investigated by many groups. These spectra have many uses, including determining shielding requirements and calculating dose for radiation protection purposes, validating nuclear data, and determining material composition. It has even been proposed to use measurements of neutron leakage spectra to determine the soil composition of Mars. As part of the Los Alamos National Laboratory (LANL) Experiments Underpinned by Computational Learning for Improvements in nuclear Data (EUCLID) Laboratory Directed Research Development (LDRD) project, the authors are working on developing methods to calculate sensitivities to neutron leakage spectra. This may allow nuclear data evaluators to better use neutron leakage spectra data to constrain or adjust nuclear data. This may also be particularly useful as many neutron leakage spectra measurements do not require fissile material and can be performed with well characterized neutron sources. Neutron leakage spectra measurements are usually performed by placing a detector at the outside of a nuclear system. This may be a reactor, a neutron generator, or a neutron source. Certain detection systems can use pulse height data to infer neutron energies, other detector systems rely on other ways of determining neutron energy (for example, a Bonner sphere with multiple moderator thicknesses can be used to measure neutron spectrum). These measurements typically rely on some unfolding of the measured results. For pulsed systems like the "Livermore Pulsed Spheres," time-of-flight information can also help to constrain the neutron energy spectra data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gravity wave structure between 60 and 90 km inferred from Space Shuttle reentry data

Density fluctuations obtained along seven Space Shuttle reentry tracks are used to examine the horizontal structure and the vertical distribution of density variance in the mesosphere and lower thermosphere. The tracks lie primarily over open ocean at middle and low latitudes and represent the only measurements of horizontal atmospheric structure at these heights available to date. The density fluctuations are interpreted in terms of gravity-wave motions and reveal significant density (and velocity) variance at horizontal scales ranging from about 10 to 1000 km. Fluctuation amplitudes are used to infer corresponding velocity perturbations and characteristic vertical scales and frequencies of the wave spectrum. Results suggest that the mean velocity variance is smaller over the Pacific ocean than over major land masses, and that the variance increases with height in a manner consistent with that expected in the present of wave saturation processes.

Fritts, David C.↗

Midlatitude Lightning NOx Production Efficiency Inferred From OMI and WWLLN Data

Oxides of nitrogen are critical trace gases in the troposphere and are precursors for nitrate aerosol and ozone, which is an important pollutant and greenhouse gas. Lightning is the major source of NOx(NO + NO2) in the mid-to upper troposphere. We estimate the production efficiency (PE) of lightning NOx (LNOx) using satellite data from the Ozone Monitoring Instrument (OMI) and the ground-based WorldWide Lightning Location Network (WWLLN) in three northern midlatitude, primarily continental regions that include much of North America, Europe and East Asia. Data were obtained over 5 boreal summers, 2007 –2011 and comprise the largest number of midlatitude convective events to date for estimating the LNOx PE with satellite NO2and ground-based lightning measurements. In contrast to some previous studies, the algorithm assumes no minimum flash-rate threshold and estimates freshly produced LNOxby subtracting a background of aged NOx estimated from the OMI dataset itself. We infer an average value of 180 ± 100 moles LNOx produced per lightning flash. We also show evidence of a dependence of PE on lightning flash rate and find an approximate empirical power function relating moles LNOxto flashes. PE decreases by an order of magnitude for a 2-order of magnitude increase in flash rate. This phenomenon has not been reported in previous satellite LNOxstudies but is consistent with ground-based observations suggesting an inverse relationship between flash rate and size.

LNOx↗

Spike-and-Slab Shrinkage Priors for Structurally Sparse Bayesian Neural Networks

Network complexity and computational efficiency have become increasingly significant aspects of deep learning. Sparse deep learning addresses these challenges by recovering a sparse representation of the underlying target function by reducing heavily overparameterized deep neural networks. Specifically, deep neural architectures compressed via structured sparsity (e.g., node sparsity) provide low-latency inference, higher data throughput, and reduced energy consumption. In this article, we explore two well-established shrinkage techniques, Lasso and Horseshoe, for model compression in Bayesian neural networks (BNNs). To this end, we propose structurally sparse BNNs, which systematically prune excessive nodes with the following: 1) spike-and-slab group Lasso (SS-GL) and 2) SS group Horseshoe (SS-GHS) priors, and develop computationally tractable variational inference, including continuous relaxation of Bernoulli variables. We establish the contraction rates of the variational posterior of our proposed models as a function of the network topology, layerwise node cardinalities, and bounds on the network weights. Furthermore, we empirically demonstrate the competitive performance of our models compared with the baseline models in prediction accuracy, model compression, and inference latency.

97 MATHEMATICS AND COMPUTING↗

Thermal Modeling and Testing of High-Temperature Refractory Ceramic Insulation Felts

Heat transfer in high-temperature, high-porosity, flexible refractory ceramic fibrous insulation felts is investigated. Heat transfer in these insulation materials consists of combined gas conduction, solid conduction, and radiation modes, with the precise theoretical modeling of the latter two modes being formidable. A semi-empirical model that requires inverse methods and steady-state thermal test data to infer some of the required model parameters is further developed in this study, and applied to five insulation materials for temperatures between 300 K and 1900 K. The steady-state thermal test setup at NASA Langley Research Center with recent modifications to increase its testing capability to 1900 K is discussed. Design considerations to ensure one-dimensional heat transfer in the test setup are described. Test data and corresponding thermal models for alumina and zirconia-based fibrous insulation felts are presented. Furthermore, test data and thermal models on two fibrous insulation samples containing additives to further suppress either radiation or gas conduction modes of heat transfer are presented. Previously published alumina-based insulation data are also re-processed using the updated modeling methodology. The significance of various heat transfer modes in typical insulation samples is discussed and used to provide general guidance on optimum insulation layups.

Kamran Daryabeigi↗

Quantum critical point followed by Kondo-like behavior due to Cu substitution in the itinerant antiferromagnet La 2 ⁢(Cu 𝑥 ⁢Ni 1−𝑥 ) 7

La 2 ⁢Ni 7 is an itinerant magnetic system with a small ordered moment of ∼ 0.1⁢µ 𝐵 /Ni and a series of antiferromagnetic (AFM) transitions at 𝑇 1 = 61.0K, 𝑇 2 = 56.5K, and 𝑇 3 = 42.2K. 𝑀 ⁡(𝐻), and 𝜌 ⁡(𝐻) isotherms as well as constant field 𝑀⁡ (𝑇) and 𝜌 ⁡(𝑇) measurements on single-crystalline samples manifest a complex, anisotropic 𝐻−𝑇 phase diagram with multiple phase lines. Here, in this study, we present the growth and characterization of single crystals of the La 2 ⁢(Cu 𝑥 ⁢Ni 1−𝑥 ) 7 series for 0 ≤ 𝑥 ≤ 0.181. We measured powder x-ray diffraction and composition, as well as anisotropic temperature- and field-dependent resistivity, temperature- and field-dependent magnetization, and temperature-dependent heat capacity on these single crystals. Using the measured data we infer a transition temperature-composition (𝑇−𝑥) phase diagram for this system to study the evolution of the AFM ordering upon Cu substitution. For 0 ≤ 𝑥 ≤ 0.097, the system remains magnetically ordered at base temperature with 𝑥 ≤ 0.012, showing signs of multiple AFM ordering temperatures. For the higher substitution levels 0.125 ≤ 𝑥 ≤ 0.181, there are no signatures of magnetic ordering, but anomalous features in resistance and heat capacity data are observed which are consistent with the Kondo effect in this system. The intermediate 𝑥 = 0.105 sample lies between the magnetic ordered and the Kondo regime and is in the vicinity of the AFM quantum critical point (QCP). Thus, La 2⁢ (Cu 𝑥 ⁢Ni 1−𝑥 ) 7 is an example of a small moment system that can be tuned through a QCP. Given these data combined with the fact that the La 2⁢ Ni 7 structure has kagomelike, Ni sublattices running perpendicular to the crystallographic 𝑐 axis, and a predicted 3⁢𝑑-electron flat band that contributes to the density of states near the Fermi energy, La 2 ⁢(Cu 𝑥 ⁢Ni 1−𝑥 ) 7 becomes a promising system to host and study exotic physics.

Das, Atreyee [Ames Laboratory, and Iowa State Univ↗

Dynamic Harris current sheet thickness from Cluster current density and plasma measurements

We use the first accurate measurements of current densities in the plasma sheet to calculate the half-thickness and position of the current sheet as a function of time. Our technique assumes a Harris current sheet model, which is parameterized by lobe magnetic field B(o), current sheet half-thickness h, and current sheet position z(sub o). Cluster measurements of magnetic field, current density, and plasma pressure are used to infer the three parameters as a function of time. We find that most long timescale (6-12 hours) current sheet crossings observed by Cluster cannot be described by a static Harris current sheet with a single set of parameters B(sub o), h, and z(sub o). Noting the presence of high-frequency fluctuations that appear to be superimposed on lower frequency variations, we average over running 6-min intervals and use the smoothed data to infer the parameters h(t) and z(sub o)(t), constrained by the pressure balance lobe magnetic field B(sub o)(t). Whereas this approach has been used in previous studies, the spatial gnuhen& now provided by the Cluster magnetometers were unavailable or not well constrained in earlier studies. We place the calculated hdf&cknessa in a magnetospheric context by examining the change in thickness with substorm phase for three case study events and 21 events in a superposed epoch analysis. We find that the inferred half-thickness in many cases reflects the nominal changes experienced by the plasma sheet during substorms (i.e., thinning during growth phase, thickening following substorm onset). We conclude with an analysis of the relative contribution of (Delta)B(sub z)/(Delta)X to the cross-tail current density during substorms. We find that (Delta)B(sub z)/(Delta)X can contribute a significant portion of the cross-tail c m n t around substorm onset.

Thompson, S. M.↗

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: measurement of the BAO and growth rate of structure of the emission line galaxy sample from the anisotropic power spectrum between redshift 0.6 and 1.1

We analyse the large-scale clustering in Fourier space of emission line galaxies (ELG) from the Data Release 16 of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey. The ELG sample contains 173,736 galaxies covering 1,170 square degrees in the redshift range 0.6 < z < 1.1. We perform a BAO measurement from the post-reconstruction power spectrum monopole, and study redshift space distortions (RSD) in the first three even multipoles. Furthermore, photometric variations yield fluctuations of both the angular and radial survey selection functions. Those are directly inferred from data, imposing integral constraints which we model consistently. The full data set has only a weak preference for a BAO feature (1.4σ). At the effective redshift zeff = 0.845 we measure $D_{\rm V}(z_{\rm eff})/r_{\rm drag} = 18.33_{-0.62}^{+0.57}$, with DV the volume-averaged distance and rdrag the comoving sound horizon at the drag epoch. In combination with the RSD measurement, at zeff = 0.85 we find $f\sigma _8(z_{\rm eff}) = 0.289_{-0.096}^{+0.085}$, with f the growth rate of structure and σ8 the normalisation of the linear power spectrum, $D_{\rm H}(z_{\rm eff})/r_{\rm drag} = 20.0_{-2.2}^{+2.4}$ and DM(zeff)/rdrag = 19.17 ± 0.99 with DH and DM the Hubble and comoving angular distances, respectively. These results are in agreement with those obtained in configuration space, thus allowing a consensus measurement of fσ8(zeff) = 0.315 ± 0.095, $D_{\rm H}(z_{\rm eff})/r_{\rm drag} = 19.6_{-2.1}^{+2.2}$ and DM(zeff)/rdrag = 19.5 ± 1.0. This measurement is consistent with a flat ΛCDM model with Planck parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty

Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. We develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)-based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one-dimensional problem and demonstrate it numerically for a three-dimensional aquifer model. Here we design a generalized preconditioned Crank-Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low-rank-based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two-dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.

54 ENVIRONMENTAL SCIENCES↗

On the Distribution and Implications of Mantled and Exhumed Terrains on Mars

The intent of this paper is to show the relationships for Mars among albedo, thermal inertia, roughness inferred from MOLA pulse width spread data, and geology inferred from photogeological analyses. Mapping of surface units using these parameters and approaches, in combination with analysis of hyperspectral image data from ISM, TES, OMEGA, and CRISM observations, will maximize our understanding of the distribution and nature of surface units on the red planet. Results will directly impact the selection of landing sites that exhibit geological records needed to understand planetary habitability.

Arvidson, R. E.↗

Mineralogy, fabric and deformation domains in D" across the southwestern border of the African LLSVP

Recent advances in seismic anisotropy studies that jointly use reflections and shear wave splitting have proven to place tight constraints on the plausible anisotropic and deformation scenarios in the D" region. We apply this novel methodology to a large area of the D" region beneath the South Atlantic, in proximity to and within the African large low seismic velocity province (LLSVP). This area of the mantle is characterized by a transition from fast to slow seismic velocity anomalies and it is thought to be the location of deep-seated plumes responsible for hotspot volcanism. Attempting to probe mantle composition and deformation along the LLSVP borders may provide key information on mantle dynamics. By analysing seismic phases sampling this region, we detect a D" discontinuity over a large area beneath the South Atlantic, with inferred depth ranges ~170 to ~240 km above the core–mantle boundary. We find evidence for a D" reflector within the area of the LLSVP. Shear wave splitting observations suggest that anisotropy is present in this region of the mantle, in agreement with previous studies that partially sampled this region. We model the observations considering lattice- and shape-preferred orientation of materials expected in the D" region. A regional variation of mineralogy, phase transition boundaries, and deformation direction is required to explain the data. We infer two distinct domains of mineralogy and deformation: aligned post-perovskite outside the LLSVP and aligned bridgmanite within the LLSVP. Finally, the scenario depicted by this study agrees well with the current hypotheses for the composition of the LLSVP and with the prevalence of vertical deformation directions expected to occur along the LLSVPs borders.

58 GEOSCIENCES↗

The interpretation of crustal dynamics data in terms of plate motions and regional deformation near plate boundaries

During our participation in the NASA Crustal Dynamics Project under NASA contract NAS-27339 and grant NAG5-814 for the period 1982-1991, we published or submitted for publication 30 research papers and 52 abstracts of presentations at scientific meetings. In addition, five M.I.T. Ph.D. students (Eric Bergman, Steven Bratt, Dan Davis, Jeanne Sauber, Anne Sheehan) were supported wholly or in part by this project during their thesis research. Highlights of our research progress during this period include the following: application of geodetic data to determine rates of strain in the Mojave block and in central California and to clarify the relation of such strain to the San Andreas fault and Pacific-North American plate motions; application of geodetic data to infer post seismic deformation associated with large earthquakes in the Imperial Valley, Hebgen Lake, Argentina, and Chile; determination of the state of stress in oceanic lithosphere from a systematic study of the centroid depths and source mechanisms of oceanic intraplate earthquakes; development of models for the state of stress in young oceanic regions arising from the differential cooling of the lithosphere; determination of the depth extent and rupture characteristics of oceanic transform earthquakes; improved determination of earthquake slip vectors in the Gulf of California, an important data set for the estimation of Pacific-North American plate motions; development of models for the state of stress and mechanics of fold-and-thrust belts and accretionary wedges; development of procedures to invert geoid height, residual bathymetry, and differential body wave travel time residuals for lateral variations in the characteristic temperature and bulk composition of the oceanic upper mantle; and initial GPS measurements of crustal deformation associated with the Imperial-Cerro Prieto fault system in southern California and northern Mexico. Full descriptions of the research conducted on these topics may be found in the Semi-Annual status Reports submitted regularly to NASA over the course of this project and in the publications listed.

Solomon, Sean C.↗

The warped neutral sheet and plasma sheet in the near-earth geomagnetic tail

An analysis of ISEE 2 plasma and magnetic field data indicates that the plasma sheet and neutral sheet in the near-earth magnetotail are warped in such a manner that in summer (winter) the neutral sheet rises above (dips below) the solar magnetospheric equatorial plane near the center of the tail, but dips below (rises above) the equatorial plane along the tail flanks. In the near tail, the neutral sheet crosses the equatorial plane at about 12 earth radii from the aberrated X axis, considerably closer to the center of the tail than has been inferred from data obtained farther downstream. This increase in the warp with decreasing distance from the earth is consistent with theoretical predictions. In the near tail, the warp is sufficiently strong when the dipole tilt angle is large that even in quiet times the upper or lower edge of the plasma sheet can be found close to the solar magnetospheric equatorial plane along the tail flanks. The seasonal dependence of the warp can produce certain dawn-dusk asymmetries in satellite data which are more apparent than real.

Gosling, J. T.↗

Mineralogical, Elemental, and Tomographic Reconnaissance Investigation for CLPS (METRIC)

METRIC is a robotic science laboratory that can determine the mineralogy, elemental chemistry, micromorphology, and thermophysical properties of planetary regolith. The METRIC suite comprises METRIC XRD/F, an X-ray diffraction/X-ray fluorescence instrument that can determine the mineralogy and elemental chemistry of regolith samples; METRIC XCT, a micro X-ray computed tomography instrument that can be used to evaluate grain/crystallite sizes and textures; METRIC IRS, an imaging spectrometer mounted on a rover that can determine mineralogy and thermophysical properties at the landing site; and a pneumatic sample collection, processing, distribution system developed by Honeybee Robotics. The payload elements could be deployed on a static lander or a rover. Data returned from the METRIC payload would inform origin, formation, and evolution of rocky planetary bodies. METRIC XRD/F draws on heritage from the CheMin instrument on the Mars Science Laboratory (MSL) Curiosity rover [1], with a few important improvements. Like CheMin, METRIC XRD/F operates in transmission geometry and uses piezoelectric actuators on sample cells in a tuning fork geometry to induce convective grain motion of the regolith to create a randomly oriented powder. MSL CheMin uses an energy-sensitive CCD to collect XRD patterns and XRF spectra simultaneously from the same sample cell, resulting in qualitative XRF data. METRIC XRD/F uses two different sample cells, one optimized for XRD and one optimized for XRF, and a silicon drift detector to detect fluoresced X-rays. This improvement to the XRF capabilities provides quantitative geochemical data of major elements down to Z = 11 and allows for the detection of minor and trace elements that are critical for evaluating geologic evolution of the Moon (e.g., P and Th). Modest improvements to the XRD geometry and hardware allow for better angular resolution and the ability to distinguish between members of the pyroxene group. METRIC XCT uses the same geometry and much of the same hardware as METRIC XRD/F, where a CCD would capture images of a regolith sample in a 3 mm diameter sample tube that is rotated 360° in steps <1°. Image brightness can be used to infer compositional data, where brighter materials indicate a higher Z, much like scanning electron microscopy. Data from METRIC XCT complement those from METRIC XRD/F. Particle size, shape, and texture can provide petrologic and provenance information, whereas vesicle size and morphology in volcanic or impact melt lithologies can inform cooling rates. METRIC IRS is a hyperspectral thermal imager that can be mounted to a lander or rover to provide mineralogical data from the broader landing site and help determine whether the samples analyzed by METRIC XRD/F and XCT are representative. The METRIC IRS spectral range (8–14 μm) and resolution (10.8 cm-1) allow for quantitative mineralogy from modelling Reststrahlen bands of major rock-forming minerals (e.g., silicates, phosphates). Radiance cubes can be processed and modelled with an onboard high-performance computer to determine mineral abundances of plagioclase, high-Ca pyroxene, pigeonite, orthopyroxene, olivine, and glass. Regolith samples can be acquired, processed, and delivered to the X-ray instruments via multiple sample handling systems, but the pneumatic sampling systems developed by Honeybee Robotics [e.g., 2] are best suited for relatively low-cost missions that are being competed for the Moon (e.g., NASA’s Payloads and Research Investigations for the Surface of the Moon program). There are pneumatic sampling systems that collect surface material and other systems that pneumatically drill up to ~1 m below the surface, providing material that has not been space weathered and has not been affected by the lander’s exhaust. [1] Blake, D. F., Vaniman, D., Achilles, C., Anderson, R., Bish, D., et al. (2012). Space Sci. Rev. 170, 341-478. https://doi.org/10.1007/s11214-012-9905-1. [2] Zacny, K., Betts, B., Hedlund, M., Long, P., Gramlich, M., Tura, K., Chu, P., Jacob, A., Garcia, A. (2014). IEEE Aerospace Conference, 3-7 March 2014, Big Sky, MT, U.S.A.

X-ray diffraction↗