SEARCH · Search NASA

Results for “surface model”

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.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Daily evapotranspiration changes during heatwaves at 32 NEON sites, 2019-2021

This dataset provides partitioned evapotranspiration (ET, the combined loss of water from soil and plant surfaces) anomalies during heatwave events—soil evaporation (E) and transpiration (T)—for 268 heatwave events across 32 National Ecological Observatory Network (NEON) flux sites in the contiguous United States from 2019–2021. Using an ensemble of four high-frequency turbulence methods (Flux-variance Similarity, Conditional Eddy Covariance [CEC], CEC with Water-Use Efficiency, and Conditional Eddy Accumulation; see Zahn and Bou-Zeid 2024), half-hourly transpiration-to-evapotranspiration (T/ET) ratios were derived from 20 hertz (Hz, cycles per second) eddy covariance measurements of carbon dioxide (CO₂) and water vapor (H₂O) concentrations. The dataset spans six vegetation types including evergreen and deciduous forests, grasslands, cultivated crops, shrublands, and emergent herbaceous wetlands. Data Package Contents: The dataset includes a single CSV (comma-separated values) file containing daily anomalies (deviations from baseline conditions) for transpiration (Delta_T), evaporation (Delta_E), total evapotranspiration (Delta_ET), and T/ET ratio (Delta_T_ET) during each day of identified heatwave events. The file also includes site codes, dates, heatwave event identifiers, and day-of-heatwave indicators. The CSV file can be opened with spreadsheet software (Microsoft Excel, Google Sheets) or programming environments (Python, R, MATLAB). This resource enables researchers to investigate ecosystem-specific responses to thermal extremes, validate land surface model partitioning of ET fluxes, and examine feedbacks between water cycling and surface energy balance during heatwaves. The dataset is particularly valuable for studies linking vegetation hydraulic strategies to climate resilience, as it captures the divergent responses of shallow-rooted versus deep-rooted ecosystems. Potential applications include improving drought early warning systems, informing irrigation management strategies, and advancing our mechanistic understanding of land-atmosphere interactions under extreme heat conditions.

Day of Heatwave↗

DOE OSTI2025

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗

DOE OSTI2025

Characterizing How Meteorological Forcing Selection and Parameter Uncertainty Influence Community Land Model Version 5 Hydrological Applications in the United States

Despite the increasing use of large-scale Land Surface Models (LSMs) in predicting hydrological responses in extreme conditions, there's a critical gap in understanding the uncertainties in these predictions. This study addresses this gap through a detailed diagnostic evaluation of the uncertainties arising from meteorological forcing selection and model parametrization in hydrological simulations of the Community Land Model version 5 (CLM5). CLM5 is configured at a spatial scale of about 12-km to simulate runoff processes for 464 headwater watersheds, selected from the Catchment Attributes for Large-Sample Studies (CAMELS) dataset to be representative of physiographic and climatic gradients across the conterminous United States. For each watershed, CLM5 is driven by five commonly used gridded forcing datasets in combination with a large ensemble (> 1200) of key CLM5 hydrologic parameters. Our results suggest that uncertainty in CLM5 runoff simulations resulting from both forcing and parametric sources is markedly higher in arid regions, e.g., Great Plains and Midwest regions. Uncertainty in low flow is dominated by parametric uncertainty, while the selection of meteorological forcing contributes more dominantly to high flow and seasonal flows during fall and spring. Our analysis also demonstrates that the selection of forcing datasets and the metrics used to calibrate CLM5 significantly impact the model’s predictive accuracy in extreme event severity for both floods and droughts. Overall, the results from this study highlight the need to understand and account for forcing and parametric uncertainties in CLM5 simulations, particularly for hazard and risk assessments addressing hydrologic extremes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Interfacial Ice Density Fluctuations Inform Surface Ice-Philicity

The propensity of a surface to nucleate ice or bind to ice is governed by its ice-philicity─its relative preference for ice over liquid water. However, the relationship between the features of a surface and its ice-philicity is not well understood, and for surfaces with chemical or topographical heterogeneity, such as proteins, their ice-philicity is not even well-defined. In the analogous problem of surface hydrophobicity, it has been shown that hydrophobic surfaces display enhanced low water-density (vapor-like) fluctuations in their vicinity. To interrogate whether enhanced ice-like fluctuations are similarly observed near ice-philic surfaces, here we use molecular simulations and enhanced sampling techniques. Using a family of model surfaces for which the wetting coefficient, k , has previously been characterized, we show that the free energy of observing rare interfacial ice-density fluctuations decreases monotonically with increasing k . By utilizing this connection, we investigate a set of fcc systems and find that the (110) surface is more ice-philic than the (111) or (100) surfaces. By additionally analyzing the structure of interfacial ice, we find that all surfaces prefer to bind to the basal plane of ice, and the topographical complementarity of the (110) surface to the basal plane explains its higher ice-philicity. Using enhanced interfacial ice-like fluctuations as a measure of surface ice-philicity, we then characterize the ice-philicity of chemically heterogeneous and topologically complex systems. In particular, we study the spruce budworm antifreeze protein (sbwAFP), which binds to ice using a known ice-binding site (IBS) and resists engulfment using nonbinding sites of the protein (NBSs). We find that the IBS displays enhanced interfacial ice-density fluctuations and is therefore more ice-philic than the two NBSs studied. We also find the two NBSs are similarly ice-phobic. By establishing a connection between interfacial ice-like fluctuations and surface ice-philicity, our findings thus provide a way to characterize the ice-philicity of heterogeneous surfaces.

crystals↗

DOE OSTI2025

TxDOT Road Elevation Model Dataset

This dataset provides three formats of Road Elevation Model (REM) data: 3D road line/polygon GeoPackage (GPKG), road lidar LAZ and COPC LAZ, and road digital surface model (DSM) GeoTIFF. Data are produced from the ~50TB TxGIO (formerly TNRIS) state lidar collections. This dataset is currently organized by maintenance section in each TxDOT district. Computation is done on GPU computing resources at Oak Ridge National Laboratory (ORNL), through a Strategic Partnership Project with UT Austin and an NSF ACCESS computing allocation award that enables fast massive data movement between TACC Corral and ORNL CADES/OLCF using Globus. In addition to this release from ORNL, a copy of this dataset can also be downloaded at https://web.corral.tacc.utexas.edu/nfiedata/road3d/.

13 HYDRO ENERGY↗

DOE OSTI2025

A unified ensemble soil moisture dataset across the continental United States

Abstract A unified ensemble soil moisture (SM) package has been developed over the Continental United States (CONUS). The data package includes 19 products from land surface models, remote sensing, reanalysis, and machine learning models. All datasets are unified to a 0.25-degree and monthly spatiotemporal resolution, providing a comprehensive view of surface SM dynamics. The statistical analysis of the datasets leverages the Koppen-Geiger Climate Classification to explore surface SM’s spatiotemporal variabilities. The extracted SM characteristics highlight distinct patterns, with the western CONUS showing larger coefficient of variation values and the eastern CONUS exhibiting higher SM values. Remote sensing datasets tend to be drier, while reanalysis products present wetter conditions. In-situ SM observations serve as the basis for wavelet power spectrum analyses to explain discrepancies in temporal scales across datasets facilitating daily SM records. This study provides a comprehensive soil moisture data package and an analysis framework that can be used for Earth system model evaluations and uncertainty quantification, quantifying drought impacts and land–atmosphere interactions and making recommendations for drought response planning.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

A multi‐variable framework for selecting WRF physics configurations at convection‐permitting scales: An Amazon wet‐season case study

Tropical convection over rainforests modulates atmospheric circulation and the energy and hydrological cycles across multiple scales. However, the scarcity of observations still limits our understanding of these processes. Although numerical models are utilized to investigate atmospheric physical processes, their performance depends on the choice of parameterizations and grid resolution. Here, in this work, we introduce and apply a multi‐variable, multi‐physics framework to evaluate and rank the Advanced Research Weather Research and Forecasting (WRF–ARW) configurations at 1‐km resolution over the central Amazon during the wet season. A 48‐member ensemble combines three land‐surface models (LSMs), four planetary boundary layer (PBL), and four microphysics (MP) schemes. Model performance for seven convection‐related near‐surface and boundary‐layer variables is assessed using Taylor diagrams and the Taylor skill score (TSS), analysis of variance (ANOVA)‐based sensitivity metrics, and non‐parametric rank tests. We then construct a combined, weighted TSS to identify configurations that are comparatively robust across variables. Results showed that most configurations reproduce near‐surface temperature, sensible and latent heat fluxes, and boundary‐layer height reasonably well, whereas humidity and rainfall remain challenging. LSM choice has the strongest impact on the surface fluxes and a secondary influence on near‐surface temperature and humidity, PBL schemes dominate boundary‐layer height, and MP schemes exert the largest control on rainfall. No single configuration is optimal for all variables, but the combination of Noah (LSM), Yonsei University (PBL), and Morrison (MP) emerges as the most robust configuration for this case, with the WRF single‐moment six‐class scheme (WSM6) providing a competitive, computationally cheaper MP alternative. The framework is general and can be applied to other regions, seasons, and convective regimes.

Amazon rainforest↗

DOE OSTI2024

Carbon-phosphorus cycle models overestimate CO 2 enrichment response in a mature Eucalyptus forest

The importance of phosphorus (P) in regulating ecosystem responses to climate change has fostered P-cycle implementation in land surface models, but their CO 2 effects predictions have not been evaluated against measurements. Here, we perform a data-driven model evaluation where simulations of eight widely used P-enabled models were confronted with observations from a long-term free-air CO 2 enrichment experiment in a mature, P-limited Eucalyptus forest. We show that most models predicted the correct sign and magnitude of the CO 2 effect on ecosystem carbon (C) sequestration, but they generally overestimated the effects on plant C uptake and growth. We identify leaf-to-canopy scaling of photosynthesis, plant tissue stoichiometry, plant belowground C allocation, and the subsequent consequences for plant-microbial interaction as key areas in which models of ecosystem C-P interaction can be improved. Together, this data-model intercomparison reveals data-driven insights into the performance and functionality of P-enabled models and adds to the existing evidence that the global CO 2 -driven carbon sink is overestimated by models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

More concentrated precipitation decreases terrestrial water storage

Terrestrial water availability is a key determinant of human and ecosystem well-being. Apart from mean precipitation and evaporation changes, it is unknown how daily-scale precipitation concentration into fewer, heavier events affects hydrologic partitioning and the land water balance. Here we show observationally that more concentrated precipitation decreases land water availability across all climates globally, a drying effect as strong in magnitude as the wetting effect of increased total precipitation. Simple and complex land-surface models recover the observed effect, whereas idealized simulations show that it arises from enhanced evaporation caused by hydrologic partitioning changes at the land surface. Projected terrestrial water storage impacts of warming-driven precipitation concentration at about 2 °C of warming shift the land surface to abnormally dry conditions (≥0.5 standard deviation) for 27% of the global population, independent of any total precipitation or irrigation changes. Our results show new key determinants of the land water balance, highlighting its sensitivity to the temporal distribution of precipitation, with broad implications for future water availability.

climate change↗

DOE OSTI2026

Distinct Kinetic Signatures of Photodesorption from Metal Nanoparticles

Visible photon fluxes can influence the rate and selectivity of heterogeneously catalyzed reactions on metal nanoparticle surfaces. Models describing the influence of photon fluxes have typically introduced photon flux dependent apparent thermal kinetic parameters (reaction orders, activation energies, binding energies, etc.). This has relied on empirical fitting of reaction rate data, making mechanistic interpretations of how photon fluxes influence elementary step rates challenging and inconsistent with fundamental descriptions of photochemistry on metal surfaces developed from surface science studies. Using the CO adsorption–desorption quasi-equilibrium reaction on Pt/Al 2 O 3 catalysts as a model system, we measured steady state adsorbed CO (CO*) coverages under isothermal and isobaric (1 mbar CO) conditions as a function of temperature (473–573 K) and of 440 nm photon flux ((0.1–5.2) × 10 3 # hv Pt site –1 s –1 ) using in situ IR spectroscopy. Steady state CO* coverage on Pt was photon flux dependent with increasing photon flux causing decreasing coverage, consistent with photons driving CO* desorption rates faster than thermal CO* desorption rates. However, photon flux dependent CO* coverages were essentially temperature independent, inconsistent with models that describe photon effects using perturbations to apparent thermal kinetic parameters. Instead, 120 steady state CO* coverages as a function of temperature and photon flux are quantitatively described by a kinetic model in which the overall desorption rate is a summation of independent thermal and photon induced CO* desorption rates. Site-resolved analysis reveals distinct kinetic parameters for photon driven desorption of CO* from well-coordinated, under-coordinated, and highly under-coordinated Pt sites, with temperature-dependent apparent quantum efficiencies (AQE) consistent with temperature dependence of vibrational quanta distribution of adsorbed CO. The rigorous kinetic rate laws for independent photon and thermal driven pathways allow for predictive modeling of the influence of photon fluxes on the rates of CO* desorption under catalytic conditions. Further, the analysis provides evidence that steady state continuous wave photon fluxes can drive desorption/adsorption reactions on metal surfaces out of thermal equilibrium, reconciling surface science observations of molecular photodesorption with applied catalysis. The work establishes a general kinetic framework to be considered for photon driven processes on metals, and defines catalyst, reaction, and photon flux characteristic design principles for breaking Sabatier limitations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Development of a plant carbon–nitrogen interface coupling framework in a coupled biophysical-ecosystem–biogeochemical model (SSiB5/TRIFFID/DayCent-SOM v1.0)

Plant and microbial nitrogen (N) dynamics and N availability regulate the photosynthetic capacity and capture, allocation, and turnover of carbon (C) in terrestrial ecosystems. Studies have shown that a wide divergence in representations of N dynamics in land surface models leads to large uncertainties in the biogeochemical cycle of terrestrial ecosystems and then in climate simulations as well as the projections of future trajectories. In this study, a plant C–N interface coupling framework is developed and implemented in a coupled biophysical-ecosystem–biogeochemical model (SSiB5/TRIFFID/DayCent-SOM v1.0). The main concept and structure of this plant C–N framework and its coupling strategy are presented in this study. This framework takes more plant N-related processes into account. The dynamic ratio (CNR) for each plant functional type (PFT) is introduced to consider plant resistance and adaptation to N availability to better evaluate the plant response to N limitation. Furthermore, when available N is less than plant N demand, plant growth is restricted by a lower maximum carboxylation capacity of RuBisCO (V c,max ), reducing gross primary productivity (GPP). In addition, a module for plant respiration rates is introduced by adjusting the respiration with different rates for different plant components at the same N concentration. Since insufficient N can potentially give rise to lags in plant phenology, the phenological scheme is also adjusted in response to N availability. All these considerations ensure a more comprehensive incorporation of N regulations to plant growth and C cycling. This new approach has been tested systematically to assess the effects of this coupling framework and N limitation on the terrestrial carbon cycle. Long-term measurements from flux tower sites with different PFTs and global satellite-derived products are employed as references to assess these effects. The results show a general improvement with the new plant C–N coupling framework, with more consistent emergent properties, such as GPP and leaf area index (LAI), compared to the observations. The main improvements occur in tropical Africa and boreal regions, accompanied by a decrease in the bias in global GPP and LAI by 16.3 % and 27.1 %, respectively.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Understanding the Active Site Structures and Achieving Catalytic Activity Tuning of Atomically Dispersed FeN 4 Sites for Oxygen Reduction Reaction

Atomically dispersed Fe–—N—C catalysts with high oxygen reduction reaction (ORR) activity have attracted great attention since the last decade. Due to comparable ORR activity and low material cost, they are promising platinum group metal (PGM)-free catalysts that can replace the commercialized Pt/C materials; furthermore, it can facilitate the efficiency of the fuel cell technologies and mitigate dependence on fossil fuels. Great advancements have been made to experimentally optimize the synthesis approach of the Fe–—N—C catalysts, enhance the ORR activity, and improve the catalyst stability. Similarly, recent theoretical studies also provide enriched understanding of the active site structures, properties, and reaction mechanisms. In this review, discussions are made upon utilizing combined experimental and computational spectroscopy to reveal the active site structures, employing mechanistic studies to investigate reaction thermodynamics and kinetics, as well as developing scaling relationships to assist the design and development of future PGM-free catalyst materials. Furthermore, recent advances in studying Fe–—N—C catalysts utilizing electrified surface models and explicit solvation models are also discussed. Not only can these aspects improve the accuracy of theoretical simulation and predictions but also deepen the understanding of the catalyst properties and reaction mechanisms under the effect of surface charges and solvent molecules.

Fe single-atom catalysts↗

DOE OSTI2025

Shrub Expansion Can Counteract Carbon Losses From Warming Tundra

Arctic warming is causing substantial compositional, structural, and functional changes in tundra vegetation including shrub and tree-line expansion and densification. However, predicting the carbon trajectories of the changing Arctic is challenging due to interacting feedbacks between vegetation composition and structure, and surface characteristics. We conduct a sensitivity analysis of the current-date to 2100 projected surface energy fluxes, soil carbon pools, and CO 2 fluxes to different shrub expansion rates under future emission scenarios (intermediate—RCP4.5, and high—RCP8.5) using the Arctic-focused configuration of E3SM Land Model (ELM). We focus on Trail Valley Creek (TVC), an upland tundra site in the western Canadian Arctic, which is experiencing shrub densification and expansion. We find that shrub expansion did not significantly alter the modeled surface energy and water budgets. However, the carbon balance was sensitive to shrub expansion, which drove higher rates of carbon sequestration as a consequence of higher shrubification rates. Thus, at low shrub expansion rates, the site would become a carbon source, especially under RCP8.5, due to higher temperatures, which deepen the active layer and enhance soil respiration. At higher shrub expansion rates, TVC would become a net CO 2 sink under both Representative Concentration Pathway scenarios due to higher shrub productivity outweighing temperature-driven respiration increase. Our simulations highlight the effect of shrub expansion on Arctic ecosystem carbon fluxes and stocks. We predict that at TVC, shrubification rate would interact with climate change intensity to determine whether the site would become a carbon sink or source under projected future climate.

Yazbeck, Theresia [The Ohio State Univ., Columbus,↗

DOE OSTI2025

Aerial and Processed Model Data Representing As-built Conditions in Coastal Port Arthur, Texas in May 2025

This dataset was collected by the Co-Design Team of the Southeast Texas Urban Integrated Field Lab, a research initiative led by the University of Texas at Austin and funded by the U.S. Department of Energy. The broader project focuses on developing climate-resilient design solutions for the Beaumont–Port Arthur region, with more information available at www.setx-uifl.org. Our team conducted aerial surveys of the Port Arthur coastal neighborhood in May 2025, before the start of construction scheduled for Summer 2026. These pre-construction datasets are designed to facilitate comparative analyses, including pre- and post-construction assessments and simulated inundation scenario evaluations. Aerial images were captured using DroneDeploy autonomous flight systems, with imagery processed through the DroneDeploy engine. All original aerial photographs are provided in JPG format and organized in zipped folders by area. The processed data package includes: 3D surface models Orthomosaics Geospatial and topographic mappings Point clouds For guidance on file contents, structure, and recommended usage, please refer to the included README file.

2D mapping↗

DOE OSTI2024

Proton-Coupled Electron Transfer Mechanisms for CO 2 Reduction to Methanol Catalyzed by Surface-Immobilized Cobalt Phthalocyanine

Immobilized cobalt phthalocyanine (CoPc) is a highly promising architecture for the six-proton, six-electron reduction of CO 2 to methanol. This electroreduction process relies on proton-coupled electron transfer (PCET) reactions that can occur by sequential or concerted mechanisms. Immobilization on a conductive support such as carbon nanotubes or graphitic flakes can fundamentally alter the PCET mechanisms. We use density functional theory (DFT) calculations of CoPc adsorbed on an explicit graphitic surface model to investigate intermediates in the electroreduction of CO 2 to methanol. Our calculations show that the alignment of the CoPc and graphitic electronic states influences the reductive chemistry. These calculations also distinguish between charging the graphitic surface and reducing the CoPc and adsorbed intermediates as electrons are added to the system. This analysis allows us to identify the chemical transformations that are likely to be concerted PCET, defined for these systems as the mechanism in which protonation of a CO 2 reduction intermediate is accompanied by electron abstraction from the graphitic surface to the adsorbate without thermodynamically stable intermediates. Furthermore, this work establishes a mechanistic pathway for methanol production that is consistent with experimental observations and provides fundamental insight into how immobilization of the CoPc impacts its CO 2 reduction chemistry.

Adsorption↗

DOE OSTI2025

Data Files for Runoff Evaluation in an Earth System Land Model for Permafrost Regions

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce.This data set includes all files that were produced and applied in the paper Runoff Evaluation in an Earth System Land Model for Permafrost Regions [Xiang et al. in review]. The paper is in review as of July 1 2025 in Geoscientific Model Development (GMD). In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-rich integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM’s parameterized representation of total runoff. This dataset contains 2 figure image files (*.png, *jpg) that describe the study site and methods, as well as folders (Figure*.zip) that contain the associated data files (*.csv, *.dat) and python code notebooks (*.ipynb) for figures 3-7 in the paper. Jupyter notebook (*.ipynb) files that produce the figure files using the associated data files will run within a python environment configured with Jupyter Lab or Notebook packages.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Best practices in software development for robust and reproducible geoscientific models based on insights from the Global Carbon Budget's dynamic vegetation models

Computational models play an increasingly vital role in scientific research by enabling the numerical simulation of complex processes. Such models are also fundamental in geosciences. For instance, they offer critical insights into the impacts of global change on the Earth system today and in the future. Beyond their value as research tools, models are also software products and should therefore adhere to certain established software engineering standards. However, scientists are rarely trained as software developers, which can lead to potential deficiencies in software quality like unreadable, inefficient, or erroneous code. The complexity of models, coupled with their integration into broader workflows, also often makes it challenging to reproduce results, evaluate processes, and build upon them. In this paper, we review the state and current practices of the development processes of the state-of-the-art land surface models used by the Global Carbon Budget. We combine the experience of modelers from the respective research groups with the expertise of software engineers from tech companies to outline key principles and tools for improving software quality in research. We explore four main areas: (1) model testing and validation, (2) scientific, technical, and user documentation, (3) version control, continuous integration, and code review, and (4) the portability and reproducibility of workflows. Our review reveals that while modeling communities are incorporating many best practices, significant room for improvement remains in areas such as automated testing, automated documentation, and reproducibility. Therefore, we here identify and promote essential software engineering practices, including numerous examples of practices from within the community that can serve as guidelines for other models and could help streamline processes across the entire community. We conclude with an open-source example implementation of these principles, demonstrating portable and reproducible data flows, a continuous integration setup, and web-based visualizations. This example may serve as a practical resource for model developers, users, and all scientists engaged in scientific programming.

Gregor, Konstantin [Technical Univ. of Munich (Ger↗

DOE OSTI2025

A surprising proliferation of detwinning in β -tin at extreme loading rates

Integrating data from dynamic compression experiments of condensed matter across three national laboratories has led to insight and quantitative calibration of materials strength over decades of loading rate. For many materials, a single strength model (such as PTW) is sufficient to capture the flow-stress strain rate relationship which is monotonic. Here, we show here that β -tin, a tetragonal metal, exhibits dramatic deviations from this behavior. Naive fitting to a single PTW model is insufficient to capture the behavior; indeed, such resulting inferred flow stress versus strain exhibits a non-monotonic behavior. We suggest a resolution to this by proposing that in β -tin there are important Bauschinger effects arising from favorable conditions for twinning and detwinning. A simple yield surface model when paired with PTW hardening captures the experimental data.

36 MATERIALS SCIENCE↗