SEARCH · Search NASA

Results for “value–chain estimation”

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.

SourceResearch connectionFrom yearThrough year

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

424 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2026

The 10 September 2025 M w 4.1 Earthquake in Northeastern Utah, United States: An Archetypal Continental Mantle Event

The 10 September 2025 M w 4.1 earthquake in northeastern Utah, United States, had a focal depth 68 km beneath sea level, which is ∼20–25 km greater than estimates of local crustal thickness, making it a rare example of a continental mantle earthquake (CME). The focal depth is well resolved from arrival-time inversion (nearest station ∼13 km away) and moment tensor inversion of regional waveforms. Similar to other CMEs in the Intermountain West, there were no obvious aftershocks or foreshocks, and the waveforms were enriched in high-frequency energy. Spectral modeling gives a stress drop of ∼80 MPa and a radiation efficiency of ∼0.08, albeit with large uncertainties. The high stress drop and low radiation efficiency are consistent with a dissipative source process such as thermal runaway. Also similar to previous Intermountain West CMEs, the event occurred along the boundary of the Archean Wyoming craton, where pressure–temperature conditions favor ductile deformation. We hypothesize that edge-driven or regional-scale mantle convection produces increased strain rates near the craton boundary that make either conventional brittle failure or thermal runaway feasible at relatively high pressure–temperature conditions. High conductivity inferred around the edge of the craton may suggest that fluids also contribute to CME occurrence.

Koper, Keith D. [Univ. of Utah, Salt Lake City, UT↗

NASA NTRS2026

Hot Water, Cold Reality: Feasibility Assessment of Iodine Removal in Heated Spacecraft Potable Water Systems

The current eXploration Potable Water Dispenser (xPWD) design removes iodine upstream of the heated leg due to concerns with the Activated Carbon and Ion Exchange (ACTEX) functionality in hot water, leaving the downstream volume without residual biocide. The NESC determined that this non-iodinated volume is a concern for microbial growth during exploration missions and proposed 33 biocide architecture options for future missions that could address this concern. The top-ranked architecture out of the report was Option 1: moving the iodine removal media as close to the dispensing needle as possible to minimize the wetted components without biocide in the xPWD. Three main challenges were identified with this proposed configuration. First, the hot water at 175 ± 25 °F is a concern for the potential physical degradation of the ion exchange resin and lowered adsorption capacity in activated carbon. Second, moving the ACTEX or alternative sorption media closer to the dispense needle increases the unheated volume downstream of the heater, challenging the ability for dispensed water to meet temperature requirements. Finally, bubbles evolved from dissolved gas coming out of solution in the heater could clog or reduce the efficiency of the sorption media. To address the first challenge, more thermally robust ion exchange resins were identified and adsorption capacity tests were planned and will be discussed in a companion ICES paper (ICES-2026-5). To address the dispense temperature concerns, allowable bed size and architectural configuration changes are proposed. The value of adding phase separators to remove bubbles and potential implementation schemes are discussed. These findings support the development of potable water systems resilient to microbial risks during long-duration space missions.

Biocide↗

NASA NTRS2026

Hot Water, Cold Reality: Feasibility Assessment of Iodine Removal in Heated Spacecraft Potable Water Systems

PWD↗

NASA NTRS2026

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

mineral dust↗

DOE OSTI2026

Robustness of inertial fusion energy relevant implosions to low-mode asymmetries

Inertial confinement fusion implosions for power production need to be robust. Degradation sources need to be understood and considered for an economic power plant. In this work, the impact of low-mode asymmetries on yield is studied using 2D radiation hydrodynamics simulations for inertial fusion energy relevant implosions. First, a high-gain indirectly driven capsule design is developed. This design is subsequently degraded with low-mode drive asymmetries. The results show that yield is constant with increasing amounts of asymmetry until a cliff is reached where the yield quickly drops below ignition. Robustness arises because the degraded hotspot still has enough time to get to ignition before capsule decompression. This robustness occurs even though the asymmetry couples less energy to the hotspot and causes decompression to occur sooner. The robustness is linked to the generalized Lawson parameter, χ⁠. The value of χ steadily decreases with increasing asymmetry. Once χ is close to one, the hotspot is no longer able to ignite before the capsule decompresses. A simple piston model is used to predict where the robustness cliff is for any design. The trade-off between yield and robustness is also demonstrated.

Physics - Plasma physics↗

DOE OSTI2024

Determining Stellar Elemental Abundances from DESI Spectra with the Data-driven Payne

Abstract Stellar abundances for a large number of stars provide key information for the study of Galactic formation history. Large spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI) and LAMOST take median-to-low-resolution (R≲ 5000) spectra in the full optical wavelength range for millions of stars. However, the line-blending effect in these spectra causes great challenges for elemental abundance determination. Here we employDD-Payne, a data-driven method regularized by differential spectra from stellar physical models, to the DESI early data release spectra for stellar abundance determination. Our implementation delivers 15 labels, including effective temperatureT eff , surface gravity log g , microturbulence velocityv mic , and the abundances for 12 individual elements, namely C, N, O, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, and Ni. Given a spectral signal-to-noise ratio of 100 per pixel, the internal precisions of the label estimates are about 20 K forT eff , 0.05 dex for log g , and 0.05 dex for most elemental abundances. These results agree with the theoretical limits from the Crámer–Rao bound calculation within a factor of 2. The majority of the accreted halo stars contributed by the Gaia–Enceladus–Sausage are discernible from the disk and in situ halo populations in the resultant [Mg/Fe]–[Fe/H] and [Al/Fe]–[Fe/H] abundance spaces. We also provide distance and orbital parameters for the sample stars, which spread over a distance out to ∼100 kpc. The DESI sample has a significantly higher fraction of distant (or metal-poor) stars than the other existing spectroscopic surveys, making it a powerful data set for studying the Galactic outskirts. The catalog is publicly available.

Astronomy & Astrophysics↗

DOE OSTI2026

Coalition for Community-Supported Affordable Geothermal Energy Systems (C2SAGES)

The C2SAGES project evaluated the feasibility of a community geothermal system for the planned Windy Ridge affordable housing development in Hinesburg, Vermont. Led by GTI Energy with Vermont Gas Systems, LN Consulting, NREL, and Frontier Energy, the work assessed technical design, energy performance, costs, business models, community engagement, maintenance, workforce development, and permitting. The proposed system was designed to serve 100% of the development’s heating, cooling, and domestic hot water loads. Compared with a baseline using air-source heat pumps and natural gas water heating, the geothermal system was estimated to reduce HVAC and domestic hot water energy use by about 45% to 48%, lower operating and maintenance costs, and reduce 30-year life-cycle costs by 37% for Phase 1 and 10% for Phase 2. Technical testing and modeling indicated that the Windy Ridge site is suitable for a community-scale geothermal system. The project also developed borehole field layouts, piping concepts, pump house designs, controls, maintenance plans, and supporting engineering drawings. The business model analysis found that first cost, ownership structure, and customer affordability remain major deployment challenges. Utility-led maintenance and operation were viewed favorably, but traditional utility cost-recovery models may require subsidy or revised financing structures to be practical for affordable housing. Community engagement highlighted the need for clear public education, transparent financing, reliable long-term maintenance, trained technicians, and the potential to pair geothermal systems with weatherization. Overall, the report concludes that community geothermal is technically feasible and offers meaningful energy, emissions, and life-cycle cost benefits, but broader deployment will depend on workable financing models and workforce readiness.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Sustainable Production of Biomass‐Derived Graphite and Graphene Conductive Inks from Biochar

Abstract Graphite is a commonly used raw material across many industries and the demand for high‐quality graphite has been increasing in recent years, especially as a primary component for lithium‐ion batteries. However, graphite production is currently limited by production shortages, uneven geographical distribution, and significant environmental impacts incurred from conventional processing. Here, an efficient method of synthesizing biomass‐derived graphite from biochar is presented as a sustainable alternative to natural and synthetic graphite. The resulting bio‐graphite equals or exceeds quantitative quality metrics of spheroidized natural graphite, achieving a RamanI D /I G ratio of 0.051 and crystallite size parallel to the graphene layers (L a ) of 2.08 µm. This bio‐graphite is directly applied as a raw input to liquid‐phase exfoliation of graphene for the scalable production of conductive inks. The spin‐coated films from the bio‐graphene ink exhibit the highest conductivity among all biomass‐derived graphene or carbon materials, reaching 3.58 ± 0.16 × 10 4 S m −1 . Life cycle assessment demonstrates that this bio‐graphite requires less fossil fuel and produces reduced greenhouse gas emissions compared to incumbent methods for natural, synthesized, and other bio‐derived graphitic materials. This work thus offers a sustainable, locally adaptable solution for producing state‐of‐the‐art graphite that is suitable for bio‐graphene and other high‐value products.

Chemistry↗

NASA NTRS1999

DS-2 Mars Microprobe Battery

In January of 1999 the NM DS-2 Mars microprobe will be launched to impact on Mars in December. The technical objectives of the missions are to demonstrate: key technologies, a passive atmospheric entry, highly integrated microelectronics which can withstand both low temperatures and high decelerations, and the capability to conduct in-situ, surface and subsurface science data acquisition. The scientific objectives are to determine if ice is present below the Martian surface, measure the local atmospheric pressure, characterize the thermal properties of the martian subsurface soil, and to estimate the vertical temperature gradient of the Martian soil. The battery requirements are 2-4 cell batteries, with voltage of 6-14 volts, capacity of 550 mAh at 80C, and 2Ah at 25C, shelf life of 2.5 years, an operating temperature of 60C and below, and the ability to withstand shock impact of 80,000 g's. The technical challenges and the approach is reviewed. The Li-SOCL2 system is reviewed, and graphs showing the current and voltage is displayed, along with the voltage over discharge time. The problems encountered during the testing were: (1) impact sensitivity, (2) cracking of the seals, and (3) delay in voltage. A new design resulted in no problems in the impact testing phase. The corrective actions for the seal problems involved: (1) pre weld fill tube, (2) an improved heat sink during case to cover weld and (3) change the seal dimensions to reduce stress. To correct the voltage delay problem the solutions involved: (1) drying the electrodes to reduce contamination by water, (2) assemblage of the cells within a week of electrode manufacture, (3) ensure electrolyte purity, and (4) provide second depassivation pulse after landing. The conclusions on further testing were that the battery can: (1) withstand anticipated shock of up to 80,000 g, (2) meet the discharge profile post shock at Mars temperatures, (3) meet the required self discharge rate and (4) meet environmental requirements.

Frank, H.↗

DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

MARVEL Reactor Digital Engineering Developments

The MARVEL reactor project has served to introduce a new generation of engineers to the processes required to transform a reactor design from simply an idea on paper into what will be an approved, constructed, and operational nuclear power system. Much as there have been advances in materials, analysis, and evaluation methodologies over the 50 years since the last reactor was built at INL, so too has the technology for managing the engineering process itself advanced. Digital Engineering tools and methods provide improved coordination between previously siloed engineering disciplines, reduced burdens of non-value-added data transcription processes and bring forward insights and improvements that might otherwise fall later in the design stage, where changes are much more costly. While the tools and techniques to support the full digital engineering vision are not yet complete, the MARVEL design processes provide valuable demonstrations and validations of key aspects and illuminate further areas for implementation by subsequent projects.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2026

Cost and Carbon Intensity Implications of Coprocessing Sustainable Aviation Fuel at Petroleum Refineries

Sustainable aviation fuel (SAF) will play a critical role in decarbonizing the aviation industry. Among SAF production pathways, alcohol-to-jet (ATJ) stands out for its scalability, supported by abundant feedstock availability and a well-established bioethanol industry. However, significant reductions in SAF carbon intensity (CI) require the use of future feedstocks (e.g., cellulosic) whose adoption is hindered by high capital costs for feedstock processing and ethanol upgrading. Here, we evaluate the financial viability and environmental implications of integrating an ATJ SAF biorefinery within a petroleum refinery, utilizing miscanthus and switchgrass as example feedstocks. Three scenarios are evaluated: standalone (benchmark), colocated, and repurposing (coprocessing SAF within the petroleum refinery). Results show repurposing reduces baseline capital costs by 36% and SAF minimum selling price (MSP) by 12% to 8.14 USD·gal −1 ; the superior performance of repurposing is consistent across both feedstocks. Integration has a limited effect on SAF CI, which remains stable across scenarios, whereas using cellulosic feedstocks reduces CI by over 70% relative to corn, with baseline values of 17.01 g CO 2 e· MJ −1 for miscanthus and 12.23 g CO 2 e·MJ −1 for switchgrass. Global sensitivity analysis reveals MSP declines with greater coprocessing levels.

09 BIOMASS FUELS↗

DOE OSTI2026

In Situ UV–vis Quantification of Chlorobasicity in Binary Molten Chlorides

Optical basicity measurements provide a basis for determining the relative chlorobasicity of a given molten salt composition. This study uses Bi3+ as a spectroscopic probe to investigate optical basicity across three binary molten chloride systems (LiCl-KCl, MgCl2–KCl, and ZnCl2–KCl) through in situ UV–vis spectroscopy measurements of the 1S0→3P1 transition. LiCl-KCl systems exhibit the highest optical basicity with minimal compositional dependence, while both divalent systems show stronger, yet nonlinear, dependence that is reminiscent of acid–base titration curves. The influence of cation electronegativity and polarizing power (Li+ < Mg2+ < Zn2+) directly reflects the basicity curves, with ZnCl2–KCl displaying the steepest gradients and lowest basicity values. Nonlinear dependencies in divalent systems indicate formation of chloride-bridged network structures with coordination states dictated by the amount of KCl in a salt composition, which was further investigated and validated by molecular dynamics simulations of polymeric structures formed in ZnCl2 and MgCl2-based salts using machine learning force fields. The temperature dependence of chlorobasicity in terms of optical basicity was investigated, and a loose, positive correlation between the two was revealed. These findings establish Bi3+ as a sensitive probe for understanding electronic environments in molten chloride systems and provide insights for predicting the chemical behavior of molten salts for use in high-temperature environments.

Buttice, Kailee [Department of Nuclear Engineering↗

NASA NTRS1993

Enhanced Decoding for the Galileo S-Band Mission

A coding system under consideration for the Galileo S-band low-gain antenna mission is a concatenated system using a variable redundancy Reed-Solomon outer code and a (14,1/4) convolutional inner code. The 8-bit Reed-Solomon symbols are interleaved to depth 8, and the eight 255-symbol codewords in each interleaved block have redundancies 64, 20, 20, 20, 64, 20, 20, and 20, respectively (or equivalently, the codewords have 191, 235, 235, 235, 191, 235, 235, and 235 8-bit information symbols, respectively). This concatenated code is to be decoded by an enhanced decoder that utilizes a maximum likelihood (Viterbi) convolutional decoder; a Reed Solomon decoder capable of processing erasures; an algorithm for declaring erasures in undecoded codewords based on known erroneous symbols in neighboring decodable words; a second Viterbi decoding operation (redecoding) constrained to follow only paths consistent with the known symbols from previously decodable Reed-Solomon codewords; and a second Reed-Solomon decoding operation using the output from the Viterbi redecoder and additional erasure declarations to the extent possible. It is estimated that this code and decoder can achieve a decoded bit error rate of 1 x 10 -7 at a concatenated code signal-to-noise ratio of 0.76 dB. By comparison, a threshold of 1.17 dB is required for a baseline coding system consisting of the same (14,1/4) convolutional code, a (255,223) Reed-Solomon code with constant redundancy 32 also interleaved to depth 8, a one-pass Viterbi decoder, and a Reed Solomon decoder incapable of declaring or utilizing erasures. The relative gain of the enhanced system is thus 0.41 dB. It is predicted from analysis based on an assumption of infinite interleaving that the coding gain could be further improved by approximately 0.2 dB if four stages of Viterbi decoding and four levels of Reed-Solomon redundancy are permitted. Confirmation of this effect and specification of the optimum four-level redundancy profile for depth-8 interleaving is currently being done.

S Dolinar↗

NASA NTRS2026

CH3OH and HCN in Interstellar Comet 3I/ATLAS Mapped with the ALMA Atacama Compact Array: Distinct Outgassing Behaviors and a Remarkably High CH3OH/HCN Production Rate Ratio

We report the detection of methanol (CH3OH) toward interstellar comet 3I/ATLAS using the At-acama Compact Array of the Atacama Large Millimeter/Submillimeter Array (ALMA) on UT 2025 August 28, September 18 and 22, and October 1, and of hydrogen cyanide (HCN) on September 12 and 15. These observations spanned pre-perihelion heliocentric distances (r H ) of 2.6 – 1.7 au. The molecules showed outgassing patterns distinct from one another, with HCN production being depleted in the sunward hemisphere of the coma, whereas CH3OH was enhanced in that direction. Statistical analysis of molecular scale lengths in 3I/ATLAS indicated that CH3OH included production from coma sources at Lp > 258 km at 99% confidence. However low signal-to-noise on long baselines, which sample emission on small spatial scales closest to the nucleus, prevented definitively ruling out CH3OH as purely a parent species. In contrast, HCN production was indistinguishable from direct nucleus sublimation. The CH3OH production rate increased sharply from August through October, including an uptick near the inner edge of the H2O sublimation zone at rH = 2 au. Compared to comets measured to date at radio wavelengths, the derived CH3OH/HCN ratios in 3I/ATLAS of 124+30−34 and 79+11−14on September 12 and 15, respectively, are among the most enriched values measured in any comet, surpassed only by anomalous solar system comet C/2016 R2 (PanSTARRS)

Comae↗

NASA NTRS2026

Space Mission Options for Reconnaissance and Mitigation of Asteroid 2024 YR4

Near-Earth asteroid 2024 YR 4 was discovered on 2024-12-27 and its probability of Earth impact in December 2032 peaked at about 3% on 2025-02-18. Additional observations ruled out Earth impact by 2025-02-23. However, the probability of lunar impact in December 2032 then rose, reaching about 4% by the end of the apparition in May 2025. James Webb Space Telescope (JWST) observations on 2025-03-26 estimated the asteroid's diameter at 60 +/- 7 m. Studies of 2024 YR 4 's potential lunar impact effects suggest lunar ejecta could increase micrometeoroid debris flux in low Earth orbit up to 1000 times above background levels over just a few days, possibly threatening astronauts and spacecraft. In this work, we present options for space missions to 2024 YR 4 that could be utilized if lunar impact is confirmed. We cover flyby & rendezvous reconnaissance, deflection, and robust disruption of the asteroid. We examine both rapid-response and delayed launch options through 2032. We evaluate chemical and solar electric propulsion, various launch vehicles, optimized deep space maneuvers, and gravity assists. Re-tasking extant spacecraft and using built spacecraft not yet launched are also considered. The best reconnaissance mission options launch in late 2028, leaving only approximately three years for development at the time of this writing in August 2025. Deflection missions were assessed and appear impractical. However, kinetic robust disruption missions are available with launches between April 2030 and April 2032. Nuclear robust disruption missions are also available with launches between late 2029 and late 2031. Finally, even if lunar impact is ruled out there is significant potential utility in deploying a reconnaissance mission to characterize the asteroid.

Asteroid Disruption↗

DOE OSTI2026

Sparse non-Markovian Noise Modeling of Transmon-Based Multi-Qubit Operations

The influence of noise on quantum dynamics is one of the main factors preventing current quantum processors from performing accurate quantum computations. Sufficient noise characterization and modeling can provide key insights into the effect of noise on quantum algorithms and inform the design of targeted error protection protocols. However, constructing effective noise models that are sparse in model parameters, yet predictive can be challenging. In this work, we present an approach for effective noise modeling of multi-qubit operations on transmon-based devices. Through a comprehensive characterization of seven devices offered by the IBM Quantum Platform, we show that the model can capture and predict a wide range of single- and two-qubit behaviors, including non-Markovian effects resulting from spatiotemporally correlated noise sources. The model’s predictive power is further highlighted through multi-qubit dynamical decoupling demonstrations and an implementation of the variational quantum eigensolver. As a training proxy for the hardware, we show that the model can predict expectation values within a relative error of 0.5%; this is a sevenfold improvement over default hardware noise models. Through these demonstrations, we highlight key error sources in superconducting qubits and illustrate the utility of reduced noise models for predicting hardware dynamics.

open quantum systems & decoherence↗

NASA NTRS1977

Magnetic and Dielectric Properties of Lunar Samples

The dielectric properties of lunar soil and rock samples show a systematic character when careful precautions are taken to ensure there is no moisture present during measurement. The dielectric constant (K) above 10 5 Hz is directly dependent on density according to the formula K = (1.93 ± 0.17) p where p is the density in g/cc. The dielectric loss tangent is only slightly dependent on density and has values less than 0.005 for typical soils and 0.005-0.03 for typical rocks. In addition to a density dependence, the loss tangent appears to be directly related to the metallic ilmenite content. These results are in good agreement with the results of the Surface Electrical Properties Experiment carried on Apollo 17. It showed a surface layer of dielectric constant 3.8 and a loss tangent of 0.008. This is presumed to be a layer of compact soil about 7 m thick, overlying a medium with a dielectric constant of roughly 7.5 and a loss tangent interpreted to be 0.035, The medium is presumed to be bedrock. These results are compatible with seismic results. The magnetic properties of lunar samples can be used to study the distribution of metallic and ferrous iron which shows systematic variations from soil-type to soil-type and a general tendency to increasing Fe°/Fe** distribution in the more highly rewelded breccias. The other magnetic characteristics can also be used to determine the distribution of grain sizes. There are a number of ways of interpreting the origin of the stable remanent magnetization in lunar samples, but several lines of evidence suggest that it is of thermal origin and was acquired at a time when the igneous rocks and breccias cooled from above 800°C in the presence of an ancient field.

D W Strangway↗