Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variable Importance”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Data-driven emulation of modal aerosol microphysics via neural operator-based modeling

The complexity and the small characteristic scales of aerosol microphysical processes pose a big challenge for accurate and efficient Earth system simulations at regional and global scales. In this work, we construct and evaluate a surrogate model: the aerosol deep operator network (ADON), a physics-inspired dual-net architecture for emulating the aerosol microphysics parameterization suite in the version 2 of the Energy Earth System Model (E3SMv2). The current version of the surrogate model is trained on a dataset comprising 9.8 million samples obtained from a global E3SMv2 simulation with the horizontal resolution of about one degree under cloud-free conditions. Incorporating domain spatial and temporal coordinates, as well as principle components extracted from training data, the dual-net surrogate model effectively captures the intricate representations of aerosol and the relationship with atmospheric state variables, achieving an R-squared score over $$95.7\%$$ for all the lognormal aerosol modes in the extrapolated regime. The validated model provides feature importance of input variables and their impact on the predictive capacity of the surrogate model in relation to the E3SM. The computational cost of online inference time deployed on CPUs and GPUs with lower precisions highlights ADON’s efficiency and potential in robust predictive modeling for large-scale Earth system computations.

Bai, Zhe↗
DOE OSTI2025

Temporal Variability in Reservoir Surface Area Is an Important Source of Uncertainty in GHG Emission Estimates

Ebullitive methane (CH 4 ) emissions in lentic ecosystems tend to concentrate at river-lake interfaces and within shallow littoral zones. However, inconsistent definitions of the littoral zone and static representations of the lake or reservoir surface area contribute to major uncertainties in greenhouse gas (GHG) emissions estimates, particularly in reservoirs with large water-level fluctuations. This study examines temporal variation in littoral and total surface areas of US reservoirs and demonstrates how different methods and data sources lead to discrepencies in reservoir GHG emissions at large scales and over time. We also explore variability in remotely sensed water occurrence according to maximum surface area, reservoir purposes, and hydrologic regions. Notably, the largest relative variability in surface area is exhibited by small reservoirs with a maximum surface area <1 km 2 and non-hydroelectric reservoirs. Additionally, we use a case study of measured CH 4 emissions from the southeastern United States (Douglas Reservoir) to illustrate the effects of varying surface area on reservoir-wide GHG estimates. Upscaled CH 4 emissions in Douglas Reservoir differed by nearly two-fold depending on the source of total surface area data and whether estimates accounted for seasonal fluctuations in surface area. During seasonal drawdown in Douglas Reservoir, relative littoral area varies non-linearly; periods of lower pool elevation (and thus larger relative littoral area) likely contribute disproportionately high CH 4 emission rates compared to the commonly sampled summer season when water levels are at full-pool elevation. Improved GHG monitoring and upscaling techniques require accounting for temporal variability in reservoir surface extent and littoral area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

How Rains and Floods Become Groundwater: Understanding Recharge Pathways With Stable and Cosmogenic Isotopes

Anthropogenic climate change leads to increased precipitation intensity and exacerbated droughts in California, challenging the reliability and drought resiliency of water supply. Storing floodwater underground via managed aquifer recharge can mitigate these effects through direct infiltration or streambed infiltration. Seasonally dry streams (arroyos) already play an important part in managing groundwater recharge to the Livermore basin (CA). Understanding how, when and where stormwater and arroyo water infiltrate is critical to effectively utilise this strategy. To track water from recent storms (water year 2022–2023, WY23) into the Livermore Valley Groundwater Basin, we analysed stable water isotopes (δ 18 O and δ 2 H) in combination with naturally occurring radioactive isotopic tracers, sulphur-35 ( 35 S, t ½ = 87 days) and tritium ( 3 H, t ½ = 12.3 years). By comparing measurements of δ 18 O, 35 S and 3 H in arroyos to precipitation and groundwater, we classified the relative age and identified source of recharge to 16 wells near two arroyos. Two wells contained water with recent recharge (from WY23) from local precipitation. One well had recent recharge from variable (precipitation and imported water) sources. One well contained imported water recharge. Three wells contained water from mixed recent and older (pre-WY23) waters, from local precipitation sources. Two wells contained recent recharge from local mine settling ponds. Seven wells had older recharge from local precipitation sources. This combination of isotopes allows us to delineate where local and imported water recharges in this highly managed basin and identify locations where managed aquifer recharge is contributing to rapid groundwater infiltration. Our combined interpretation of isotopic water ages and sources in the context of land use shows that local infiltration of precipitation in open spaces is an important recharge mechanism, in addition to the managed arroyo recharge. Finally, a broader familiarity with 35 S will enable more extensive research on the infiltration of urban floodwaters.

58 GEOSCIENCES↗
DOE OSTI2024

Identifying recharge sources and their impacts on a North Central New Mexico shallow aquifer using unsupervised machine learning

In this article, shallow aquifers are important but highly variable resources in arid to semi-arid regions. Limited shallow aquifer volume results in high sensitivity to recharge fluctuations, which can impact the local fauna and flora, and transport of contaminants in the aquifer or vadose zone. Aquifer response to external forcing (e.g., precipitation) is usually solved by estimating aquifer parameters and running physics-based models to match known fluctuations of hydraulic head. However, this technique is time and computationally expensive. Furthermore, high aquifer complexity decreases precision in physics-based models. Alternatively supervised machine learning is used to predict aquifer dynamics. However, these techniques rely on input data and struggle to interpret aquifer response for missing sources (i.e., snowpack data). To counter these problems, we propose an unsupervised machine learning technique (NMFk) to estimate the impact of different sources on aquifer recharge. NMFk is used to understand the influence of external forcing on shallow aquifer recharge in the Pajarito Plateau (Los Alamos, NM, USA). The results show how NMFk can be used to reduce the data dimension in a complex field dataset to three recharge signals that cause fluctuations within the field data. Here, the source signals are interpreted as rainfall, snowmelt, and a delayed aquifer response to the previous two signals. These results evidence how heterogeneous aquifers delimited by canyons incised into the Pajarito Plateau respond in similar ways to the source signals identified by NMFk. Furthermore, results show the importance of the local geology where faults act as sinks, and anthropogenic disturbances can facilitate infiltration amplifying the interpreted signal.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Variation in Biomass Yield and Cell Wall Composition in Switchgrass Natural Variants Under Two Nitrogen Regimes

Switchgrass (Panicum virgatum) is a promising lignocellulosic biofuel crop for which biomass and processing quality are important. Inherent plant variability across genotypes and environments challenges uniformity and product quality. In this study, the impact of nitrogen (N) application on switchgrass yield and quality was examined under field conditions using a highly diverse switchgrass panel over a 4-year period at Knoxville, TN. Overall, biomass production was correlated between low (0 kg of added N/ha) and moderate (135 kg of added N/ha) nitrogen treatments, suggesting that the N impact is largely uniform across the genotypes. Nonetheless, high biomass genotypes were identified with high nitrogen-use efficiency ; biomass was congruent or even higher (up to 9-fold) in the low N treatment. Genotypes were also identified with up to 94% nitrogen-remobilization efficiency. Furthermore, nitrogen application appeared to impact lignin content in whole tillers but was neutral to lignin monomer syringyl-to-guaiacyl (S/G) ratio. The same panel grown under natural conditions in Watkinsville, GA, produced significantly more biomass than the Tennessee panel in the first 3 years, but biomass was similar across both sites and treatments in year 4. Top performing genotypes overlapped between sites by 20-37%. There were low correlations in lignin content in whole tillers across the two field sites, but moderate correlations were observed for S/G ratios. The high yielding genotypes from the low N plot identified in this study can be used in breeding programs and management strategies in switchgrass to evade adverse environmental and economic effects.

09 BIOMASS FUELS↗
DOE OSTI2025

Constrained or unconstrained? Neural-network-based equation discovery from data

Throughout many fields, practitioners often rely on differential equations to model systems. Yet, for many applications, the theoretical derivation of such equations and/or the accurate resolution of their solutions may be intractable. Instead, recently developed methods, including those based on parameter estimation, operator subset selection, and neural networks, allow for the data-driven discovery of both ordinary and partial differential equations (PDEs), on a spectrum of interpretability. The success of these strategies is often contingent upon the correct identification of representative equations from noisy observations of state variables and, as importantly and intertwined with that, the mathematical strategies utilized to enforce those equations. Specifically, the latter has been commonly addressed via unconstrained optimization strategies. Representing the PDE as a neural network, we propose to discover the PDE (or the associated operator) by solving a constrained optimization problem and using an intermediate state representation similar to a physics-informed neural network (PINN). The objective function of this constrained optimization problem promotes matching the data, while the constraints require that the discovered PDE is satisfied at a number of spatial collocation points. We present a penalty method and a widely used trust-region barrier method to solve this constrained optimization problem, and we compare these methods on numerical examples. Our results on several example problems demonstrate that the latter constrained method outperforms the penalty method, particularly for higher noise levels or fewer collocation points. This work motivates further exploration into using sophisticated constrained optimization methods in scientific machine learning, as opposed to their commonly used, penalty-method or unconstrained counterparts. For both of these methods, we solve these discovered neural network PDEs with classical methods, such as finite difference methods, as opposed to PINNs-type methods relying on automatic differentiation. Here, we briefly highlight how simultaneously fitting the data while discovering the PDE improves the robustness to noise and other small, yet crucial, implementation details.

Data-driven discovery↗
DOE OSTI2026

Machine learning of factors for improving oyster hatchery production

Oyster aquaculture and restoration in the Chesapeake Bay are vital, yet hatcheries frequently struggle with inconsistent larval growth and sudden mass mortality events. Unpredictable disruptions in larval production cause large economic losses, represent a perceived risk to growers, and impede industry expansion. To better understand associations between production yield and its potential predictors, we applied machine learning (random forest, and neural network) and statistical (generalized additive model) models to a comprehensive dataset of environmental, water quality, and operational parameters from a Maryland oyster hatchery, aiming to identify key yield predictors and develop a robust forecasting tool. We used recursive Boruta algorithm for variable selection, pinpointing critical predictors, and employed cross-validation to fine-tune model settings. Shapley value analysis offered crucial insights into model interpretations, highlighting week number, Normalized Difference Vegetation Index, salinity, turbidity, and fecundity as primary drivers of yield variability. For low-yield cases, salinity-related variables were particularly important. Our findings provide an early warning system for potential production downturns, empowering hatchery operators to make data-driven decisions for optimizing water conditions, feeding schedules, and broodstock management. By boosting predictability and efficiency, this research directly supports economic stability of the oyster industry and ecological health of the Chesapeake Bay.

Vishwakarma, Srishti [Oak Ridge National Laborator↗
DOE OSTI2025

Key Environmental and Ecological Variables of Wetland CH 4 and CO 2 Fluxes Change With Warming

Wetlands are important ecosystems for the global carbon cycle, impacting regional and global methane (CH 4 ) and carbon dioxide (CO 2 ) budgets. This study examines how environmental and ecological variables impact wetland CH 4 flux and net ecosystem exchange of CO 2 (NEE) across 17 sites globally. We also quantified the importance of variables for each wetland type and site at monthly scale under normal and warm temperatures using dominance analysis. We identified soil and air temperature (TS, TA, respectively) as key variables influencing wetland CH4, and latent heat (LE) and shortwave radiation (SW) for NEE under normal and warm conditions. However, the importance of some variables shifted with warming. For predicting the variability of wetland CH4 flux under warming, gross primary productivity (GPP) and LE, replacing wind direction (WD), were dominant variables for tropical swamps, while NEE was important for high-latitude fens and bogs under warm temperatures. For wetland NEE, the role of TA and TS decreased across all wetland types with warming, while vapor pressure deficit (VPD) became more important for mid and high-latitude wetlands. Our results reveal the complex responses of wetland carbon flux to environmental and ecological variables with warming and provide new insights into improving wetland models by incorporating additional variables and accounting for the changing roles of variables in carbon flux under warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Explainable Machine Learning for Functional Data

Black-box machine learning models are recognized as useful tools for prediction applications, but the algorithmic complexity of some models causes interpretation challenges. Explainability methods have been proposed to provide insight into these models, but there is little research focused on supervised modeling with functional data inputs. We argue that, especially in applications of high consequence, it is important to explicitly model the functional dependence in a black-box analysis to not obscure or misrepresent patterns in explanations. As such, we propose the V ariable importance E xplainable E lastic S hape A nalysis (VEESA) pipeline for training supervised machine learning models with functional inputs. The pipeline is an analysis process that includes the data preprocessing, modeling, and post-hoc explanations. The preprocessing is done using elastic functional principal components analysis, which accounts for vertical and horizontal variability in functional data and, ultimately, allows for explanations in the original data space that identify the important functional variability without bias due to correlated variables. Here, we demonstrate the pipeline on two high-consequence applications: explosives classification for national security and inkjet printer identification in forensic science. The applications exhibit the VEESA pipeline’s ability to provide an understanding of the characteristics of the functional data useful for prediction. Code for implementing the pipeline is available in the veesa R package (and supplemental python code).

Elastic Shape Analysis↗
DOE OSTI2026

The relative importance of wind and hydroclimate drivers in modulating the interannual variability of dust emissions in Earth system models

Windblown dust emissions are controlled by near-surface wind speed and sediment erodibility, the latter modulated by hydroclimate and land-use conditions. Accurate representations of these drivers are critical for reproducing historical dust variability and projecting future dust changes in Earth system models (ESMs). This study examines the discrepancies among 21 ESMs in the relative importance of wind speed versus five hydroclimate drivers in explaining the historical (1980–2014) variability of dust emissions from global drylands. In hyperarid areas, models show poor agreement in the simulated dust variability, with only 9 % out of 210 inter-model comparisons exhibiting significant positive correlations. In contrast, arid and semiarid areas exhibit a dual pattern driven by a “double-edged sword” effect of land surface memory: models with coherent hydroclimate variability show better agreement, whereas those with divergent hydroclimate representations show larger disagreement. While the ESMs capture the dominant role of wind speed in hyperarid areas, they diverge markedly in the relative contributions of wind and hydroclimate drivers in arid and semiarid areas. Replacing the Zender et al. (2003) dust scheme with the Kok et al. (2014) scheme in CESM and E3SM generally strengthens hydroclimate influences while reducing wind speed contributions to simulated dust variability. MERRA-2 reanalysis produces stronger wind influences than most ESMs across all dryland regions. These results underscore the need for improved near-surface wind simulations in hyperarid areas and more realistic land surface and hydroclimate representations in arid and semiarid areas to reduce uncertainties in global dust emission simulations.

Li, Xinzhu [Michigan Technological University, Hou↗
DOE OSTI2025

Evapotranspiration Partitioning Using Flux Tower Data in a Semi-Arid Ecosystem

Information about evapotranspiration (ET) and its components, that is, evaporation and transpiration, is crucial for a wide range of water and ecosystem management applications. However, partitioning ET into its two components is often challenging because of their spatiotemporal variabilities and lack of process understanding. This study developed a machine learning (ML) framework to shed light on ET processes and assess the relative importance of different drivers by incorporating hydrometeorology and biomass productivity variables. The Shapley Additive Explanations (SHAP) approach was applied to enhance explainability and rank the importance of ET drivers and their components. A total of 62 variables covering hydrometeorological and biomass productivity dimensions were considered from the Reynolds Creek Critical Zone Observatory (CZO) station in Idaho. The variable importance assessment identified the leading drivers individually for evaporation, transpiration and ET (soil water content for evaporation, vapour pressure deficit for transpiration and soil water content for ET). The results further highlighted the value of combining hydrometeorological and biomass productivity variables to achieve better predictability of ET processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Convergence Criteria for Multiphysics Simulations

The behavior of engineered systems is often influenced by multiple physical phenomena, such as mechanical deformation, heat transfer, and chemical species transport and reactions. There are often strong interactions between these phenomena, and there is increasing interest in applying coupled-physics models to improve understanding of physical behavior under complex environmental conditions. Multiple simulation frameworks that facilitate coupled-physics simulations are in widespread use, and these employ a variety of techniques to account for interactions between those physics. Many frameworks solve the physics models independently and transfer results between them. Alternatively, a single monolithic system of equations for every physics model can be formed and solved. Each of these approaches has its benefits and drawbacks, and the optimal approach varies depending on the nature of the problem. The open-source MOOSE framework was developed targeting solution of large-scale multiphysics problems. Although it provides options for all these coupling approaches, its standard approach for multiphysics solutions is to form and solve a single monolithic system of equations containing the unknowns for all physics models. MOOSE provides a streamlined approach for users to define the solution variables, the terms in the partial differential equations pertaining to each variable, and interactions between solution variables. One aspect of the monolithic solution approach that can be problematic, however, is defining appropriate convergence criteria for the nonlinear system. A standard approach is to determine convergence is to simply take a norm of the residual vector corresponding to the full vector of unknowns. However, if the residual vector contains variables for multiple physics models, the magnitudes of those variables can differ significantly, and the variables can converge at significantly different rates from each other. It is important to ensure that the variables for each of the physics are converged, and also ensure that the convergence criteria are not excessively stringent in cases when there is little change in the solution. This talk presents representative multiphysics problems to highlight these issues, and shows strategies for convergence criteria in MOOSE that are robust for multiphysics models under a variety of conditions.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Using feature importance as an exploratory data analysis tool on Earth system models

Abstract. Machine learning (ML) models are commonly used to generate predictions, but these models can also support the discovery of new science. Generating accurate predictions necessitates that a model captures the structure of the underlying data. If the structure is properly extracted, ML could be a useful exploratory and evidential tool. In this paper, we present a case study that demonstrates the use of ML for exploratory data analysis (EDA) in the climate space. We apply the ML explainability method of spatiotemporal zeroed feature importance (stZFI) to understand how climate-variable associations evolve over space and time. Our analyses focus on data from ensembles of Earth system models (ESMs) which provide data on different climate states and conditions. We elect to work with ESM ensembles since they allow us to compare feature importance across alternative scenarios not available with observed data. The ensembles also account for natural variability so that we can distinguish between signal and noise due to natural climate variability when computing feature importance. The use of perturbed initial condition ensembles introduces variability mimicking the natural variability in the atmosphere; thus the signals emerging using feature importance (FI) can be evaluated against the natural variability in the climate system. For our analyses, we consider the 1991 volcanic eruption of Mount Pinatubo, which was a large stratospheric aerosol injection. We explore the climate pathway associated with the eruption from aerosols to radiation to temperature at both the near-surface and stratospheric levels. In addition to applying the method to data generated from two different ESMs, we apply stZFI to reanalysis data to compare the associations identified by stZFI. We show how stZFI tracks the importance of aerosol optical depth over time on forecasting temperatures. This case study illustrates usefulness of an ML tool (stZFI) for EDA on a well-studied climate exemplar.

Ries, Daniel (ORCID:0000000250294647)↗
DOE OSTI2025

Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles

In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.

Adsorption-desorption process↗
DOE OSTI2025

Impact of a dynamic grid mix and climate on operational carbon emissions modeling for different building typologies and climate zones

Calculating operational carbon emissions through a building’s lifecycle is complex due to the dynamic nature of influencing factors such as climate and energy grid mix. This paper introduces a novel methodology for modeling 30-year operational carbon impacts of buildings and applies this method to mid-rise office and residential typologies across various US climate zones. The method accounts for these temporal variabilities using new and scarcely cited data sources. Key findings indicate that future changes in the climate, while impactful, play a relatively modest role in operational carbon emissions compared to significant reductions with modeling scenarios using the projected decarbonization of the electricity grid. Here, the study also finds that using annual, month-hourly, or hourly grid emission factors have a minimal impact on carbon accounting, except in certain climates and program types where emission patterns do not align with a building’s energy consumption. Warmer climates like Miami, Florida and Tucson, Arizona, which rely heavily on cooling, demonstrate larger variations in carbon emissions when using higher temporal resolution emission factors. Ultimately, this study underscores the critical role of grid decarbonization in reducing long-term emissions and the importance of incorporating this variable in life cycle assessment (LCA) modeling.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

The Freundlich isotherm equation best represents phosphate sorption across soil orders and land use types in tropical soils of Puerto Rico

Biomass production in the lowland wet tropical forest is greater than in any other biome, and it is typically limited by soil phosphorus (P) availability. However, the mechanisms involved in the P cycle remain poorly represented in Earth System Models (ESMs). Soil P sorption processes are key in the P cycle and for understanding the extent of P limitation for plant productivity. Currently, a few ESMs include isotherm equations to model these processes. Although the Langmuir equation is widely cited, other isotherm equations may better describe sorption in tropical soils. Here, we use a diverse range of soil samples from Puerto Rico to test the validity of the Langmuir, Freundlich, and Temkin equation. We found that across four soil orders (Inceptisols, Mollisols, Oxisols, Ultisols), and forested and cultivated land use types, the Freundlich equation best represented soil P sorption. Furthermore, the Langmuir and the Temkin equations poorly represent soil P adsorption, especially at low P concentrations. Specifically, the Langmuir equation underestimated soil P adsorption by 40% and the Temkin equation overestimated adsorption by 76%. We also found, as expected, that soil clay content and pH were the most important parameters explaining the variability of the Freundlich (K f ) constant. Greater clay content and lower pH, common in highly weathered Ultisols and Oxisols which are abundant in the tropics, led to greater K f values. Overall, our results suggest that a diversity of soils can prompt underestimation of P sorption when using the Langmuir isotherm, which leads to an overestimation of available P that can have repercussions on ESM predictions of the P cycle and tropical forest productivity.

Earth system models↗
DOE OSTI2026

Systematic Evaluation of Atmospheric Forcing, Surface Datasets, and Mesh Effects on Kilometer-Scale Land Surface and River Modeling

Earth system models are advancing toward kilometer-scale resolution to capture local climate impacts and extremes. High-resolution land and river modeling depends on multiple factors, including mesh, surface datasets, and atmospheric forcing, but their relative effects at kilometer scales remain unquantified. We evaluated five Energy Exascale Earth System Model land and river configurations over the Mid-Atlantic region using two mesh (1/8° structured versus variable-resolution unstructured mesh), two surface datasets (default versus newly developed), and three atmospheric forcings (NLDAS2, MSWX, GSWP). Evaluation against satellite, reanalysis, and in situ benchmarks across water, energy, and carbon cycles quantifies how these factors affect model performance. Forcing selection produces the largest bias reductions (12-99% across variables), followed by surface datasets (7-75%) and mesh (up to 21%). Forcing effects vary by variable, with MSWX reducing biases for snow water equivalent, evapotranspiration, albedo, temperature, and gross primary productivity, GSWP for snow cover and runoff, and NLDAS for soil moisture and streamflow. The use of newly developed surface datasets improves gross primary productivity (58% bias reduction) and evapotranspiration but increase soil moisture and albedo biases due to current modeling limitations. Variable-resolution unstructured mesh improves the simulation of small-basin streamflow through better capturing drainage networks, though mesh minimally affects other land variables. These findings provide important guidance for high-resolution modeling development and actionable science.

Land and River modeling↗
DOE OSTI2025

Development of a military-specific mesh-type computational phantom library and its application to internal dosimetry and in-field radiological triage screening

Estimates of organ-absorbed and committed doses to individuals exposed to radioactive materials via acute inhalation often rely on internal dose coefficients and detector responses from reference human computational models. To achieve more accurate dose assessments to United States Armed Forces service members exposed in-field, computational models with varying morphometric parameters representative of this population are necessary. The International Commission on Radiological Protection (ICRP) Publication 145 provides detailed mesh reference computational phantoms (MRCPs) for adult males and females, with morphometric parameters matched to the 50th percentile. Previously, these phantoms were 2D and 3D scaled to match desired height, mass, and secondary anthropomorphic parameters in the creation of the University of Florida / Memorial Sloan Kettering (UF/MSK) computational phantom library. To achieve body fat percentage targets required for accession into the US Armed Forces, muscle and fat volumes were adjusted accordingly, thus, creating the UF/Department of Defence computational phantom library presented in this study. A comprehensive library of mesh-type computational human phantoms was created, including 57 adult males and 49 adult females with morphometric parameters aligned with United States Armed Forces service members. Phantoms were restricted to a body mass index between 19 and 27.5, with body fat percentages below 26% for males and 36% for females. Specific absorbed fractions were computed for selected source and target combinations, demonstrating how variations in height and body mass influence energy absorption in target regions relative to the ICRP MRCPs. Radiation detector responses were also computed, revealing that higher body masses resulted in decreased registered counts in the detection volume. These findings highlight the importance of incorporating morphometric variability in computational phantoms to achieve more accurate dose assessments and radiation detection responses for United States Armed Forces service members who inhale radioactive materials in-field.

computational phantoms↗