Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “METRICS”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Methods development towards automated, physics-informed, quantitative quality control of TRISO-SiC

Tristructural-isotropic (TRISO) fuel particles have been developed as a high-performance fuel for use in high-temperature gas-cooled reactor (HTGR) systems due to their high efficiency and stability under both normal and off-normal conditions. Broader deployment of this technology in advanced nuclear applications may benefit from quantitative quality assurance and quality control (QA/QC) methods that directly link TRISO properties to downstream performance. A key TRISO property is the SiC layer microstructure, which influences fission product retention during irradiation. However, existing QA/QC for the TRISO-SiC microstructure comprises only a qualitative visual inspection; therefore, there is a clear opportunity for the development of quantitative methods for TRISO QA/QC. Here, to this end, previous work has demonstrated an image processing approach to grain boundary (GB) identification and subsequent extraction of microstructural metrics; however, extensive twinning within the SiC layer complicates such analyses because twin GBs significantly influence microstructural metrics but are not expected to contribute to fission product transport. This study presents the initial development, training, and testing of an ML-based image segmentation algorithm designed to identify and remove twin GBs from standard backscattered electron micrographs, providing an industrially applicable, quantitative, and physically meaningful QA/QC approach for the TRISO-SiC microstructure. Although pixel-wise performance metrics for the twin predictions are low, the change in grain area and the number of GB pixels after twin removal predicted by the ML workflow are within 1% of the true values calculated using crystallographic data. This suggests that the model is well capable of predicting overall twin boundary structures and grain morphology, and continued advancement of this approach could enable automated, scalable, and physics-informed QA/QC for TRISO-SiC microstructures, supporting the reliable qualification of coated particle fuels for next-generation reactor systems.

Computer vision↗
DOE OSTI2020

A lightweight method for evaluating in situ workflow efficiency

Performance evaluation is crucial to understanding the behavior of scientific workflows. In this study, we target an emerging type of workflow, called in situ workflows. These workflows tightly couple components such as simulation and analysis to improve overall workflow performance. To understand the tradeoffs of various configurable parameters for coupling these heterogeneous tasks, namely simulation stride, and component placement, separately monitoring each component is insufficient to gain insights into the entire workflow behavior. Through an analysis of the state-of-the-art research, we propose a lightweight metric, derived from a defined in situ step, for assessing resource usage efficiency of an in situ workflow execution. By applying this metric to a synthetic workflow, which is parameterized to emulate behaviors of a molecular dynamics simulation, we explore two possible scenarios (Idle Simulation and Idle Analyzer) for the characterization of in situ workflow execution. In addition to preliminary results from a recently published study [11], we further exploit the proposed metric to evaluate a practical in situ workflow with a real molecular dynamics application, i.e., GROMACS. Here, experimental results show that the in transit placement (analytics on dedicated nodes) sustains a higher frequency for performing in situ analysis compared to the helper-core configuration (analytics co-allocated with simulation).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effect of alloying on intrinsic ductility in WTaCrV high entropy alloys

Tungsten (W) exhibits desirable properties for extreme applications, such as the divertor in magnetic fusion reactors, but its practicality remains limited due to poor formability and insufficient irradiation resistance. In this work, we study the intrinsic ductility of body-centered cubic WTaCrV based high entropy alloys (HEAs), which are known to exhibit excellent irradiation resistance. The ductility evaluations are carried out using a criterion based on the competition between the critical stress intensity factors for emission (K Ie ) and cleavage (K Ic ) in the {110} slip planes and {110} crack planes, which are evaluated within the linear elastic fracture mechanics framework and computed using density functional theory. The results suggest that increasing the alloying concentrations of V and reducing the concentrations of W can significantly improve the ductility in these HEAs. The elastic anisotropy for these HEAs is analyzed using the Zener anisotropy ratio and its correlation with the concentration of W in the alloys is studied. Results indicate that these alloys tend to be fairly isotropic independently from the concentration of W in them. The computed data for the elastic constants of these HEAs is also compared against available experimental data. The results are in good agreement, validating the robustness and accuracy of the computational methods. Multiple phenomenological ductility metrics were also computed and analyzed against the analytical model. Some metrics, mainly the surrogate D parameter, show a good correlation with the Rice model. The potential of these empirical metrics to serve as surrogate screening models for optimizing the compositional space is also discussed.

Anisotropy↗
DOE OSTI2023

Automated reactor physics analysis framework of High Flux Isotope Reactor low-enriched uranium silicide dispersion fuel designs

The High Flux Isotope Reactor (HFIR) is a versatile research reactor that provides one of the highest steady-state neutron fluxes of any reactor in the world. The HFIR reactor physics team investigated the conversion of the current 93 wt% highly enriched uranium U 3 O 8 -Al dispersion fuel to a 19.75% low-enriched uranium (LEU) U 3 Si 2 -Al dispersion fuel. The team continuously develops a Python module to streamline the analysis steps required for an LEU core design to ensure reproducible and agile design iteration. The Python module automates the data processing between analysis steps and automates the input perturbation for branch calculations and design changes. The automated framework has proven to significantly increase the efficiency and reproducibility of the reactor physics team to design High Flux Isotope Reactor (HFIR) LEU cores and thoroughly analyze performance metrics, safety metrics, and thermal safety margins. Consequently, the team can now respond rapidly to fuel fabrication engineer and thermal-hydraulic-structural analyst requests. Numerous combinations of LEU fuel designs are explored, of which two LEU fuel designs are presented here in this paper: a low density silicide design, and a high-density silicide design. Results show that both designs meet or exceed safety and performance metrics with exception for minor differences caused by the hardened spectrum from LEU.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Measurement and analysis of GPU-accelerated applications with HPCToolkit

To address the challenge of performance analysis on the US DOE’s forthcoming exascale supercomputers, Rice University has been extending its HPCToolkit performance tools to support measurement and analysis of GPU-accelerated applications. To help developers understand the performance of accelerated applications as a whole, HPCToolkit’s measurement and analysis tools attribute metrics to calling contexts that span both CPUs and GPUs. To measure GPU-accelerated applications efficiently, HPCToolkit employs a novel wait-free data structure to coordinate monitoring and attribution of GPU performance. To help developers understand the performance of complex GPU code generated from high-level programming models, HPCToolkit constructs sophisticated approximations of call path profiles for GPU computations. To support fine-grained analysis and tuning, HPCToolkit uses PC sampling and instrumentation to measure and attribute GPU performance metrics to source lines, loops, and inlined code. To supplement fine-grained measurements, HPCToolkit can measure GPU kernel executions using hardware performance counters. To provide a view of how an execution evolves over time, HPCToolkit can collect, analyze, and visualize call path traces within and across nodes. Finally, on NVIDIA GPUs, HPCToolkit can derive and attribute a collection of useful performance metrics based on measurements using GPU PC samples. Here, we illustrate HPCToolkit’s new capabilities for analyzing GPU-accelerated applications with several codes developed as part of the Exascale Computing Project.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Abiotic influences on continuous conifer forest structure across a subalpine watershed

Understanding the abiotic drivers of high-elevation forest physiognomy is essential for forecasting how mountain ecosystems will respond to emerging environmental pressures. Most prior studies of these relationships have relied on small samples of the full landscape, resulting in limited power to detect dominant covariates and their interactions. Here we report the first evaluation of abiotic influences on a complement of accurate, wall-to-wall estimates of conifer forest structure and composition at the watershed scale. In a subalpine conifer domain in the Colorado Rocky Mountains (USA), we developed a novel method for deriving stand structure metrics from waveform LiDAR data, which showed high fidelity with field inventory. We quantified the relationships between structural and compositional metrics and climate, topographic, edaphic, and geologic factors. Our results showed that peak snow water equivalent (SWE), snow disappearance rate, and elevation explained most of the variation in forest structure. The highest stand density, basal area, maximum canopy height, and quadratic mean diameter occurred in sites with SWE around one standard deviation below mean, but with long snow residence times. Stand density decreased linearly with elevation, while other metrics peaked between 3000 m and 3200 m. Substrate properties had weaker influence. Continuous mapping of through-canopy forest structure enabled our novel findings of the dominant role of snowpack in explaining structural and compositional variation, and of elevation thresholds. Our reproducible approach facilitates assessment of forest-topoclimate relationships in other conifer-dominated landscapes and improves understanding of the baseline patterns controlling forest structure, which is needed for predicting long-term ecological change.

forest ecology↗
DOE OSTI2022

Exposure pathways and bioaccumulation of per- and polyfluoroalkyl substances in freshwater aquatic ecosystems: Key considerations

We report that due to the bioaccumulative behavior, toxicity, and recalcitrance to degradation, per- and polyfluoroalkyl substances (PFAS) are a focus for many researchers investigating freshwater aquatic ecosystems. PFAS are a diverse set of chemicals that accumulate and transport quite differently in the environment depending on the length of their fluoroalkyl chains and their functional groups. This diversity in PFAS chemical characteristics combined with varying environmental factors also impact the bioaccumulation of these compounds in different organisms. In this review, we evaluate environmental factors (such as organic carbon, proteins, lipids, and dissolved cations) as well as PFAS characteristics (head group, chain-length, and concentration) that contribute to the significant variation seen in the literature of bioaccumulation metrics reported for organisms in aquatic ecosystems. Of the factors evaluated, it was found that PFAS concentration, dissolved organic matter, sediment organic matter, and biotransformation of precursor PFAS tended to significantly impact reported bioaccumulation metrics the most. Based on this review, it is highly suggested that future studies provide sufficient details of important environmental factors, specific organism traits/ behavior, and PFAS concentrations/compounds when reporting on bioaccumulation metrics to further fill data gaps and improve our understanding of PFAS in aquatic ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A holey cavity for single-transducer 3D ultrasound imaging with physical optimization

Within the compressive sensing (CS) framework, one effective way to increase the likelihood of successful signal reconstruction is to employ random processes in the construction of the sensing matrix. This work presents a 3D holey cavity, with diverse frequency modes, to spectrally code, that is, randomize, the ultrasound wave fields. The simulated results show that the use of such a cavity enables imaging simple or complex targets, such as spheres or the letter E, by only a single transceiver—something that is not possible without the use of a coding structure like the cavity. The effect of noise on imaging results and the size of the targets on the first-order Born approximation (BA) are also investigated. Moreover, this study attempts to optimize the cavity, based on a single numerical metric, such as the sum of singular values (SSV) or mutual coherence (MC). Yet, it will be shown that neither of these metrics can consistently compare the norm-one imaging performance between two cavities of different materials or hole sizes. This leaves finding a quantitative metric for these kinds of optimizations an open problem.

42 ENGINEERING↗
DOE OSTI2025

Environmental and socio-economic Pareto-front trade-off analysis of U.S. PET packaging material in a circular economy

Various recycling technologies are emerging to implement circular economy in plasticssupply chain systems. However, the environmental and socio-economic trade-offs of in circular economy are not well understood at a systems level. Particularly, quantifying these trade-offs as a function of end-of-life (EOL) management decisions, including transition of recycling technologies, systems level metrics such as circularity, recycled content, and the need for fossil-derived plastics are not well understood. Here, the present study addressed these research gaps by applying a systems analysis modeling approach that utilizes material flow analysis, life cycle assessment, socioeconomic data, and system optimization techniques for polyethylene terephthalate (PET) packaging supply chains in the United States. Pareto-front trade-offs between conflicting environmental and socio-economic impacts as well as those between socioeconomic impacts and circularity were explored using the epsilon constraint method. The Pareto-front trade-off analysis revealed the transition of EOL management strategies for PET packaging systems, including changes in selection of recycling technologies, to aid decision making process by quantifying studied system metrics. Transitioning from environmentally optimal to socio-economically optimal systems led to increased employment (by 17%), wages (by 26%), and revenues (by 6%) but also led to increased global warming potential (GWP; by 65%), energy consumption (by 59%), and reliance on fossil PET in the system (by 78%). Finally, the results show that there is not a unique set of recycling technologies to achieve a sustainable circular economy of PET packaging system, instead it depends on the decision maker’s objectives and targeted metrics of the system.

54 - ENVIRONMENTAL SCIENCES/GLOBAL CLIMATE CHANGE ↗
DOE OSTI2025

A physiologically based pharmacokinetic (PBPK) model to align dosimetry of the isobutyl metabolic series in rats and humans

Here, we developed a physiologically based pharmacokinetic (PBPK) model in rats and humans for the isobutyl metabolic series including isobutyl acetate, isobutanol, isobutyraldehyde, and isobutyric acid. Chemical manufactures routinely use these compounds as solvents, for chemical synthesis, as potential biofuels, de-icing fluids, and additives for food and/or fragrance in consumer products. Human exposure to isobutyl compounds can occur through inhalation or oral routes. We previously developed a PBPK model for the propyl metabolic series and utilized it as a framework to create the isobutyl PBPK model due to the chemical similarities between the two series. To support model development, we measured in vitro metabolism of isobutyl acetate in rat and human blood and liver S9 fractions. Compared to rats, humans demonstrated faster isobutyl acetate hydrolysis in liver S9 fractions, while the hydrolysis rates in blood were similar between the two species. We used concentrations of isobutyl compounds measured in air and blood from rats exposed to isobutyl acetate and isobutanol as well as other published data to further parameterize the model. Following exposure to either isobutyl acetate or isobutanol, we observed isobutanol concentrations highest among the isobutyl compounds in the blood of rats. Overall, the model accurately predicts measured time course concentrations of isobutyl acetate, isobutanol, and isobutyric acid in available data in rats and humans. Sensitivity analyses identified alveolar ventilation rates, isobutyl metabolism rates, and cardiac output as the most sensitive parameters affecting concentrations of isobutyl compounds in blood. The isobutyl PBPK model enables comparisons of internal dose metrics across various isobutyl compound exposures and species and allows for calculation of equivalent external exposures that result in the same dose metric. Regulators can employ this PBPK model to predict and align internal dose metrics of isobutyl compounds for risk assessment purposes.

2-methyl-1-propanol↗
DOE OSTI2025

MS25: Materials Science-Focused Benchmark Data Set for Machine Learning Interatomic Potentials

Here, we present MS25, a benchmark data set for evaluating machine learning interatomic potentials (MLIPs) across diverse materials-relevant systems including MgO surfaces, liquid water, zeolites, a catalytic Pt surface reaction, high-entropy alloys (HEAs), and disordered Zr-oxides. Five MLIP architectures (MACE, NequIP, Allegro, MTP, and Torch-ANI) are trained and tested, focusing not only on traditional metrics (energies, forces, and stresses) but also explicitly validating derived physical observables such as lattice constants, volumes, and reaction barriers. We find that most models reach comparable accuracy on standard error metrics across the simple systems, although equivariant MLIPs offer 1.5–2× improvements over nonequivariant MLIPs in energy and force error for structurally complex or compositionally disordered environments such as HEAs and Zr–O systems. Our analysis highlights that low errors in energy and force predictions do not guarantee reliable observables, emphasizing the necessity of explicit validation. We demonstrate limitations in cross-framework transferability, as models trained on one zeolite framework (CHA) fail to reliably generalize to predictions of structurally distinct frameworks (e.g., MFI). Size-extensive tests show some dependence on system size for MgO, resulting from forced periodicity. The HEA and Zr–O data sets are identified as challenging tests for future benchmarks and MLIP model architecture developments as they show significant differentiation in error between MLIP architectures and are still relatively difficult at 1000 training images. Moving forward, we recommend that benchmarking efforts shift their focus from marginal accuracy improvements in energy and force errors toward identifying and understanding model failure modes, rigorously assessing transferability, and evaluating how their errors affect observable predictions. For researchers looking to choose an MLIP architecture, we suggest selecting equivariant MLIP architectures if the complexity of the system is a challenge. For simple materials problems, auxiliary features such as integration with molecular dynamics engines, trade-offs between computational data set generation cost vs MLIP inference speed, and framework integration may play a more important decision factor than small differences in error metrics that are unlikely to matter for production-level research.

chemical structure↗
DOE OSTI2021

Dynamics-Based Peptide–MHC Binding Optimization by a Convolutional Variational Autoencoder: A Use-Case Model for CASTELO

An unsolved challenge in the development of antigen-specific immunotherapies is determining the optimal antigens to target. Comprehension of antigen–major histocompatibility complex (MHC) binding is paramount toward achieving this goal. Here, we apply CASTELO, a combined machine learning-molecular dynamics (ML-MD) approach, to identify per-residue antigen binding contributions and then design novel antigens of increased MHC-II binding affinity for a type 1 diabetes-implicated system. We build upon a small-molecule lead optimization algorithm by training a convolutional variational autoencoder (CVAE) on MD trajectories of 48 different systems across four antigens and four HLA serotypes. We develop several new machine learning metrics including a structure-based anchor residue classification model as well as cluster comparison scores. ML-MD predictions agree well with experimental binding results and free energy perturbation-predicted binding affinities. Moreover, ML-MD metrics are independent of traditional MD stability metrics such as contact area and root-mean-square fluctuations (RMSF), which do not reflect binding affinity data. Finally, our work supports the role of structure-based deep learning techniques in antigen-specific immunotherapy design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Bottom-Up Simulation, Reconstruction, and Quantification of Macromolecule Sequences from Experimental Polymerizations

Motivated by the canonical sequence–structure–function paradigm, tools to characterize chemical patterning in natural biomacromolecules, from proteins to nucleic acids, have grown exponentially in recent years. However, analogous strategies for synthetic macromolecules remain in nascent stages, complicated by sequence polydispersity and analytical limitations. To address this, we have developed a comprehensive and open-source Python package, PRISM (polymer rate insights and sequence modeling), an end-to-end workflow that provides a path from experimental kinetics measurements to quantitative and qualitative metrics for describing chemical patterning in stochastic polymers. First, a numerical integration strategy was constructed to simulate and fit experimental data from reversible addition–fragmentation chain transfer (RAFT) polymerization kinetics, enabling the facile estimation of relevant reactivity ratios. These ratios were then used in a mechanism-specific stochastic kinetic simulation strategy to simulate sequence ensembles corresponding to model systems spanning experimental copolymers, classes of statistical polymers (e.g., alternating, block, and gradient), and multiblock copolymers. Lastly, inspired by sequence homology metrics from bioinformatics, we introduce visualization strategies and quantitative metrics to facilitate comparisons of different sequence ensembles. As the sequence–structure–function paradigm becomes increasingly central in de novo design of synthetic macromolecules, this toolkit provides a first step toward accurate and representative sequence description and featurization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Data Assimilative Optimization of WSA Source Surface and Interface Radii using Particle Filtering

TheWang-Sheeley-Arge (WSA) model estimates solar wind speed and interplanetary magnetic field polarity in the inner heliosphere using global photospheric magnetic field maps. WSA employs the Potential Field Source Surface (PFSS) and Schatten Current Sheet (SCS) models to determine the Sun's global coronal magnetic field. The PFSS and SCS models are connected through two radial parameters, the source surface and interface radii, which specify the overlap region between the inner SCS and outer PFSS models. Though both radii values are adjustable, they have typically been fixed to 2.5 solar radii. Our work highlights how solar wind predictions improve when the radii are allowed to vary over time. Data assimilation using particle filtering (sequential Monte Carlo) is used to infer optimal values over a fixed time window. Solar wind model predictions and satellite observations are compared with a newly developed quality-of-agreement prediction metric. The agreement metric between the model and observations is assumed to correspond to the probability of the two key WSA model parameters, the source surface and interface radii, where the highest metric value implies the optimal radii.We find that the optimal particle filter values of solar radii can perform twice as well as standard values for an exploratory period during Carrington Rotation 1901, with these values also reducing nonphysical kinking effects seen in solar magnetic field lines. Data assimilation choices of input realization and time frame have implications for variation in the solar wind over time. We present this work's theoretical context and practical applications for prediction accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Investigating Tropical Versus Extratropical Influences on the Southern Hemisphere Tropical Edge in the Unified Model

Abstract Since the late 1970s, observations have shown a widening of the tropical Hadley cell (HC) circulation. State‐of‐the‐art climate models reproduce the general trend along with a projected continuous expansion. Discrepancies in expansion rates of observation‐ and model‐based studies have been attributed to differences in applied methods, natural variability and model shortcomings. Furthermore, the driving influence of tropical or extratropical processes on these changes is not well understood. All of this highlights the dynamical mechanisms and the region of origin controlling the tropical width are still insufficiently understood. Here we examine the influence of systematic model biases of the atmosphere‐only Unified Model (UM) onto the simulation of the Southern Hemisphere (SH) tropical edge. We utilize nudged experiments with prescribed sea surface temperatures, where potential temperature and horizontal winds are relaxed back to ERA‐Interim reanalysis for a 20‐year period in selected regions. Correcting model biases in the tropics and extratropics separately allows us to dissect the dominant remote impacts of present model errors onto the SH tropical edge simulation. The experiments are applied to established tropical width metrics ranging from near‐surface to upper‐level metrics capturing the poleward flank of the HC. We find both regions work remotely to reduce errors in the UM fields and location of the tropical edge. Surprisingly, correcting the extratropical biases, south of 45°S, more consistently improves the tropical width across the metrics and seasons than nudging the tropics (10°N–10°S). These findings demonstrate the substantial role of extratropical influences in locating the SH tropical edge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Framework for Multisector Scenarios of Outcomes for Well-Being and Resilience

Shared community scenarios of societal and environmental system changes have underpinned a broad range of research and assessment studies over the past several decades. These scenarios have largely aimed to address specific questions within broad issue areas like climate change or biodiversity and generally provided information on the drivers of change. The consequences of those drivers, such as impacts on society and policy responses, have tended to be left to the research community to investigate, using scenarios of drivers as inputs to their studies, producing projections of a disparate set of relevant output metrics. While this approach has had many benefits, it has fallen short of producing a robust, comparable literature describing outcomes across studies in common metrics. We argue that new scenarios are needed that extend current approaches to be organized around common outcome metrics for the well-being and resilience of society and ecosystems. We propose an approach that would focus on agreed upon outcomes for well-being and resilience as well as critical drivers of change, cut across issues and scales in multiple sectors, and draw on new systematic methods of scenario generation and discovery to highlight scenarios that are most critical in understanding societal risks and responding to them. Research derived from this outcome-based scenario development approach would facilitate improved assessment of risks of and responses to a range of stressors and the multi-sector interactions they generate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Fungal Spore Seasons Advanced Across the US Over Two Decades of Climate Change

Abstract Phenological shifts due to climate change have been extensively studied in plants and animals. Yet, the responses of fungal spores—organisms important to ecosystems and major airborne allergens—remain understudied. This knowledge gap limits our understanding of their ecological and public health implications. To address this, we analyzed a long‐term (2003–2022), large‐scale (the continental US) data set of airborne fungal spores collected by the US National Allergy Bureau. We first pre‐processed the spore data by gap‐filling and smoothing. Afterward, we extracted 10 metrics describing the phenology (e.g., start and end of season) and intensity (e.g., peak concentration and integral) of fungal spore seasons. These metrics were derived using two complementary but not mutually exclusive approaches—ecological and public health approaches, defined as percentiles of total spore concentration and allergenic thresholds of spore concentration, respectively. Using linear mixed‐effects models, we quantified annual shifts in these metrics across the continental US. We revealed a significant advancement in the onset of the spore seasons defined in both ecological (11 days, 95% confidence interval: 0.4–23 days) and public health (22 days, 6–38 days) approaches over two decades. Meanwhile, total spore concentrations in an annual cycle and in a spore allergy season tended to decrease over time. The earlier start of the spore season was significantly correlated with climatic variables, such as warmer temperatures and altered precipitations. Overall, our findings suggest possible climate‐driven advanced fungal spore seasons, highlighting the importance of climate change mitigation and adaptation in public health decision‐making.

Environmental Sciences & Ecology↗
DOE OSTI2024

Surface Variability Mapping and Roughness Analysis of the Moon Using a Coarse–Graining Decomposition

The lunar surface contains a wide variety of topographic shapes and features, each with different distributions and scales, and any analysis technique to objectively measure roughness must respect these qualities. Coarse-graining is a naturally scale-dependent filtering technique that preserves scale-dependent symmetries and produces coarse elevation maps that gradually erase the smaller features from the original topography. In this study of the lunar surface, we present two surface variability metrics obtained from coarse-graining lunar topography: fine elevation and coarse curvature. Both metrics are isotropic, deterministic, slope-independent, and coordinate-agnostic. Fine (detrended) elevation is acquired by subtracting the coarse elevation from the original topography and contains features that are smaller than the coarse-graining length-scale. Coarse curvature is the Laplacian of coarsened topography, and naturally quantifies the curvature at any scale and indicates whether a location is elevated or depressed relative to its neighborhood at that scale. We find that highlands and maria have distinct roughness characteristics at all length-scales. Our topographic spectra reveal four scale-breaks that mark characteristic shifts in surface roughness: 100, 300, 1,000, and 4,000 km. Comparing fine elevation distributions between maria and highlands, we show that maria fine elevation is biased toward smaller-magnitude elevations and that the maria–highland discrepancies are more pronounced at larger length-scales. Here, we also provide local examples of selected regions to demonstrate that these metrics can successfully distinguish geological features of different length-scales.

58 GEOSCIENCES↗