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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Signatures of aerosol-cloud interactions in GiOcean: a coupled global reanalysis with two-moment cloud microphysics

Aerosols influence the Earth's radiative balance through direct interactions with radiation and by affecting cloud properties. Anthropogenic aerosols have led to cooling during the industrial era through aerosol–cloud interactions (ACI), including aerosol effects on cloud microphysical properties and the subsequent adjustments. However, large uncertainties remain in Earth system models (ESMs) regarding the magnitude of this cooling. In part, ESMs substantially disagree on cloud properties, thermodynamics, the hydrological cycle, and general circulation. Reanalysis provides a useful avenue for exploring the impact of ACI on clouds and radiation because its atmosphere is forced to match realistic conditions through the assimilation of observations. Here, we explore the impact of ACI on clouds in the GiOcean reanalysis – the first to incorporate aerosol-cloud interactions. We contrast variables important for ACI between GiOcean and satellite observations and develop 2-dimensional lookup tables of ACI for both using a source-sink budget perspective to attribute the changes in cloud droplet number (Nd) and liquid water path (LWP) to aerosol and meteorology. A compositing analysis using lookup tables shows that GiOcean captures key aspects of aerosol–cloud–precipitation interactions, including (1) activation of aerosol into cloud droplets, (2) effective precipitation scavenging of Nd, (3) suppression of precipitation by high Nd in regions with heavy aerosol emissions. In contrast, satellite observations do not exhibit clear patterns for processes (2) and (3). Random Forest analysis shows that interannual variability in Nd and LWP over the Northern Hemisphere ocean in GiOcean is primarily driven by precipitation, consistent with satellite observations.

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

Optimization based process modeling of an anaerobic membrane bioreactor system: Application to swine wastewater

To maintain current levels of consumption in the economy with the dwindling supply of non-renewable material and energy, alternative resource streams more traditionally viewed as waste streams must be considered. Fermentation of high-strength wastewaters is one such pathway that allows for the recovery of energy, nitrogen, phosphorus, and carbon compounds. Anaerobic membrane bioreactors (AnMBRs) are an emerging technology that allow for the digestion of wastewater in a much smaller footprint than traditional anaerobic digesters. Adoption of this technology into industry has been limited by membrane capital and cleaning costs, but these costs may be offset through the recovery of valuable products. To evaluate the viability of AnMBR technology in the context of swine wastewater treatment, an optimization-based process model built upon Anaerobic Digestion Model No. 1 (ADM1) has been developed. Modeling results show that a swine wastewater stream provides potential for net positive energy generation from the AnMBR system in most cases. Sensitivity analyses around important variables were conducted to determine focus areas for future research into AnMBR technology and evaluate the robustness of the model to microbial variables that may change with different microbial communities.

09 BIOMASS FUELS↗
DOE OSTI2025

Exploring the Potential of Rainwater Harvesting for Cooling Towers: A Systematic Review and Regional Feasibility Assessment

Water scarcity has driven interest in rainwater harvesting, especially for U.S. industries like power and manufacturing, which dedicate a large percentage of their water needs to cooling towers. This paper provides a comprehensive systematic review of studies on the use of rainwater in cooling tower applications, along with an assessment of the regional feasibility of integrating RWH in U.S. manufacturing cooling towers. The systematic review examines the technical, economic, environmental, and policy feasibility of utilizing harvested rainwater, and the regional feasibility analysis evaluates the practical implementation of RWH in manufacturing cooling towers, considering factors such as regional RWH potential, water costs, state policies, and the industrial water use of manufacturing facilities per state. This study supports the hypothesis that RWH for cooling towers is technically feasible, economically viable, and environmentally beneficial. Harvested rainwater is naturally less conductive and soft, and rainwater reuse minimizes the ecological footprint. Supportive state policies, regional RWH potential, and rising water costs across the U.S. are important variables that may impact RWH adoption. The review highlights rooftop RWH as the most studied method and notes that implementing RWH requires infrastructure changes and filtration techniques. While initial investment costs may be high, operational and maintenance costs are low, making RWHS economically feasible over time. Regions with higher water costs and supportive policies are more likely to benefit from RWH adoption. The study provides a foundation for understanding the potential for using RWH in industrial cooling towers.

Cooling towers↗
DOE OSTI2025

Controlling matric potential in microfluidics to examine microbial dynamics in unsaturated porous media

The use of microfluidics for the study of soil microbial ecology is an emerging field. Most microfluidic studies of biological systems, however, have been performed under fully saturated conditions that are not representative of natural soil. Therefore, while microfluidics offer many unique capabilities that other methodologies cannot, they are not currently suited to address the effects of matric potential, an important variable defining the microbial moisture niche. Here, a methodology is presented that allows the user to control the aqueous conditions within microfluidic networks by manipulating matric potential using a hanging water column. The method relies on hydrophilic surface treatment of the microfluidic device using polyvinyl alcohol (PVA) and incorporating a bed of small pores at the network boundaries, which serve as a porous ceramic plate analogue (PPA). The method was validated on a simple capillary bundle and then on a more complex pore network. A water retention curve, exhibiting hysteresis, was generated for the pore network over a narrow matric potential range of 0 to – 5 kPa. Both the drainage and wetting curves were reproducible, as were the spatial configuration and the number of fragmented moisture niches in the pore network, particularly on the drainage curve. In contrast, the wetting curve exhibited greater variability in spatial configuration due to the “ink bottle effect,” where capillarity was interrupted by wider pore bodies. Ultimately, the methodology provides realistic pore-scale moisture conditions that can be easily manipulated and maintained, enabling new opportunities to explore soil biophysics and microbial biogeography in unsaturated porous media. As a brief example, images showing the localization of fluorescently tagged Pantoea sp. YR343 at −4.3 kPa are presented, highlighting bacterial distributions in water films and air-water interfaces.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Sculpting 2D Crystals via Membrane Contractions before and during Solidification

When phospholipids crystallize within the otherwise fluid membranes of giant unilamellar vesicles, the resulting molecularly thin “2D” solids exhibit great variety in their morphology evolution. For example, within membranes containing moderate amounts of the crystallizing component, crystals grow with a fixed morphology depending on vesicle size. Conversely for membranes containing large amounts of the crystallizing species, we find small compact crystals on vesicles of all sizes. However, on large vesicles, growing crystals sprout flower petals that lengthen progressively. These behaviors result from two combined mechanisms: first, like other 2D solids, the shear rigidity of phospholipid crystals renders them intolerant to morphologies with nonzero Gaussian curvature. As a result and especially at elevated membrane tension, the cost of bending elasticity is reduced at the expense of line energy by the formation of flowers as opposed to compact crystals. Second, the composition-dependent tension rise during cooling relaxes via water permeation of the membrane with a time constant scaling as R2. The amount of crystal formed for a small decrease in temperature determines this composition-dependent increase in stress from thermal contractions versus solidification. Surface Evolver computations were motivated using the predicted tension evolution to develop a processing space that maps to experimental observations for initial and growing crystal morphology. Important variable groups are identified, including a scaled ratio of bending to line energy, a vesicle-size-independent group for membrane contractions, and a time constant for stress relaxation. Though processing stresses ultimately relax, the crystal morphology persists well beyond the processing window.

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

Temperature and Water Levels Collectively Regulate Methane Emissions From Subtropical Freshwater Wetlands

Abstract Wetlands are the largest and most climate‐sensitive natural sources of methane. Accurately estimating wetland methane emissions involves reconciling inversion (“top‐down”) and process‐based (“bottom‐up”) models within the global methane budget. However, estimates from these two model types are inherently interdependent and often reveal substantial discrepancies. To enhance the reliability of both approaches, we need a comprehensive understanding of wetland methane emissions and an independent high‐resolution long‐term flux data set. Here, we employed a data‐driven random forest approach to identify key variables influencing methane emissions from subtropical freshwater wetlands in the Southeastern United States. The model‐estimated monthly mean methane fluxes fit well with measured methane fluxes ( R 2 = 0.67) at four representative FLUXNET‐CH4 wetland sites across the region. Variable importance analysis highlighted the sensitivity of subtropical freshwater wetland methane emissions to variations in both temperature and water levels. High temperatures facilitate methanogenesis by enhancing microbial activities, while elevated water levels maintain anaerobic conditions necessary for methane production. Notably, the response of methane emissions to water level fluctuations is contingent on temperature conditions, and vice versa. Moreover, we constructed the first high‐spatial‐resolution (∼1 km × 1 km) and long‐term (1982–2010) gridded regional wetland methane flux product for the Southeastern United States, estimating annual methane emissions from subtropical freshwater wetlands in the region at 4.93 ± 0.11 Tg CH 4 yr −1 for 1982–2010. This new benchmark product holds promise for validating and parameterizing uncertain wetland methane emission processes in bottom‐up models and provides improved prior information for top‐down models.

He, Keqi [Earth and Climate Sciences Nicholas Scho↗
DOE OSTI2024

Use of Satellite, Surface Observations and Numerical Weather Prediction Model Data to Improve Cloud Base Height and Cloud Base Vertical Velocity Estimation

Cloud base height (CBH) and cloud base vertical velocity (CBVV) are important variables that impact the overall climate in a region as they influence the formulation, longevity, and evolution of clouds. Retrieval of both parameters have long used ground instrumentation (e.g., Doppler lidar (DL), ground base radar); however, retrieving CBH from satellites is particularly challenging given that space-based instruments only observe cloud tops. In this manuscript, CBH is retrieved using a multi-linear regression equation, while CBVV used a random forests model. Both retrievals combine satellite and numerical weather prediction data. The satellite data used are the Visible Infrared Imaging Radiometer Suite imagery, while measurements of CBH and CBVV include DL and radiosonde data at the Southern Great Plains (SGP) Atmospheric Radiation Measurement observatory. Data from 83 summer days (May-August) in 2018–2021 featuring cumulus clouds forced by solar heating were examined and used to train the models, with years 2022–2023 used for validation. Various spatial domains were defined with one large (2.4° longitude by 2.0° latitude) SGP domain being split into smaller sections (smallest being 0.99° and 0.61° longitude and latitude respectably). CBH and CBVV values obtained from the DL as compared to the models show root mean square errors between 150 and 200 m, with CBVV values between 0.45 and 1 ms -1 . Finally, it was found that the CBH formulation performs well over all domains, while the CBVV retrievals become less accurate due to more turbulence being introduced into the observations as the number of DL stations decreases in the smaller domains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Machine learning and process-based modeling of spatiotemporal changes in active layer thickness across Alaska

Permafrost degradation poses a growing threat to infrastructure stability and ecosystem resilience in the rapidly warming Arctic. We investigated the spatiotemporal dynamics of active layer thickness (ALT) across Alaska by integrating field observations, environmental datasets, a physically based Stefan model, and machine learning (ML) techniques. Using weather projections from the Coupled Model Intercomparison Project Phase 6 under two Shared Socioeconomic Pathways (SSP 2-4.5 and SSP 5-8.5), we assessed ALT sensitivity to projected future weather conditions. The random forest (RF) model outperformed the Stefan approach in predicting ALT on the training dataset (R² = 0.84 vs. 0.53) but demonstrated lower generalizability on the test dataset (R² = 0.24 vs. 0.54). The root mean square error (RMSE) for the RF model for training and testing ranged from 14 to 22 cm, compared to 17 and 18 cm for the Stefan model. Variable importance analysis revealed that mean annual temperature and slope angle were the strongest predictors of ALT, accounting for 19% and 18% of the variance, respectively, followed by sediment transport index (14%) and stream power index (11%). Comparative analysis of baseline ALT predictions showed the Stefan model tended to project a thicker active layer (mean ± SD: 65 ± 16 cm), compared to the RF model (mean ± SD: 59 ± 8.8) cm). Both models indicated a latitudinal gradient in ALT, with shallower depths at higher latitudes. Projected ALT increases by 2100 were estimated at 3.3 ± 2.2 cm under SSP 2-4.5 and 5.9 ± 4.0 cm under SSP 5-8.5 for the ML model, whereas the Stefan model projected substantially larger increases of 13 ± 2.6 cm (SSP 2-4.5) and 28 ± 4.4 cm (SSP5-8.5). Spatial analysis showed the greatest ALT increases in northern Alaska, with relatively smaller changes in southern regions. These findings highlight the complex, multifactorial nature of ALT dynamics and the value of hybrid modeling approaches. As rising temperatures accelerate permafrost thaw, changes in ALT can disrupt ecosystems, damage infrastructures, and enhance the release of stored soil carbon, highlighting the urgent need for improved predictive capabilities to inform adaptation strategies in the Arctic.

Climate sciences↗
DOE OSTI2025

Exploring the role of clouds in offshore wind potential off the US West Coast in a changing climate

To meet US goals of deploying additional wind energy as part of the decarbonization strategy, wind plants are being planned for the deep water offshore the western US. The wind flow in that region is complex due to the proximity to the coast, cold water upwelling, and persistent stratiform clouds that interact with radiation in ways that have the potential to destabilize the atmosphere. That flow has the potential to change with a changing climate. To address these issues, we assess the flow and the clouds in that region using downscaled climate model data, under both historic climate (1975–2005) and projected future (2025–2055) conditions. We note that the climate simulations agree fairly well with the cloud patterns observed by satellite data in the nearshore and offshore regions. We then assess the projected changes in clouds, wind speed, and other important variables, noting that our simulations project that the predominant north/northwesterly low-level jet is expected to strengthen and clouds are likely to be commensurately enhanced, although projected changes are within about 10% of current conditions. Our examination of the dynamics associated with the changes in the climate simulations provides confidence in the dynamical consistency of these projected changes.

17 WIND ENERGY↗
DOE OSTI2025

Enhancing Multi-Step Reservoir Inflow Forecasting: A Time-Variant Encoder–Decoder Approach

Accurate reservoir inflow forecasting is vital for effective water resource management. Reliable forecasts enable operators to optimize storage and release strategies to meet competing sectoral demands—such as water supply, irrigation, and hydropower scheduling—while also mitigating flood and drought risks. To address this need, in this study, we propose a novel time-variant encoder–decoder (ED) model designed specifically to improve multi-step reservoir inflow forecasting, enabling accurate predictions of reservoir inflows up to seven days ahead. Unlike conventional ED-LSTM and recursive ED-LSTM models, which use fixed encoder parameters or recursively propagate predictions, our model incorporates an adaptive encoder structure that dynamically adjusts to evolving conditions at each forecast horizon. Additionally, we introduce the Expected Baseline Integrated Gradients (EB-IGs) method for variable importance analysis, enhancing interpretability of inflow by incorporating multiple baselines to capture a broader range of hydrometeorological conditions. The proposed methods are demonstrated at several diverse reservoirs across the United States. Our results show that they outperform traditional methods, particularly at longer lead times, while also offering insights into the key drivers of inflow forecasting. These advancements contribute to enhanced reservoir management through improved forecasting accuracy and practical decision-making insights under complex hydroclimatic conditions.

58 GEOSCIENCES↗
DOE OSTI2025

Characterization of Arsenic and Selenium in Coal Fly Ash to Improve Evaluations for Disposal and Reuse Potential (Final Technical Report)

Coal fly ash is a high volume waste material that is discarded in landfills and surface water impoundments across the U.S. and is also widely recycled for a variety of applications. The leaching of potential of contaminants of concern, such as arsenic (As) and selenium (Se), is often the driver of risk assessments for coal ash disposal and reuse. The extent of leachable As and Se depends on several factors related to environmental conditions and fly ash characteristics. Previous studies employed various methods to delineate the concentration, chemical form, and distribution of As and Se in fly ash materials. However, few studies have attempted to directly correlate these properties to mobilization parameters relevant to disposal and reuse. Instead, the coal residuals industries often rely upon standardized leaching protocols that can be laborious or involve hazardous chemicals. The goals of the project were to: 1) Develop and evaluate a characterization protocol that can be used to screen fly ash samples for leachability of As and Se; 2) Characterize As, Se, and associated constituents of fly ash particles at multiple length scales (nanometer to micrometer) to determine if elemental associations differ as a function of the resolution of characterization; and 3) Establish a predictive model for the chemical composition of coal ash produced annually at major U.S. coal fired power facilities on 50-year national coal supply records. For the first objective, we performed leaching experiments with 52 fly ash samples collected from 15 different U.S. power plants and representing coal feedstocks from the three major domestic coal regions. For this work, we assessed the mobilization potential of As and Se in fly ash based on standardized leaching protocols and performed multivariate and lasso regression analyses to explore correlations of leachable As and Se contents with characteristics such as major element contents, loss on ignition (LOI) and pH. The results of regression models indicated that major elements (Fe, Ca, Al) for a wide range of fly ashes can serve as predictor variables for the leaching potential of As, but not for Se. LOI and pH were not important predictive variables in the models. Both regression approaches resulted in relatively strong fits for leachable As (correlation coefficient R 2 = 0.78 for both models) compared to models for leachable Se (R 2 = 0.49). Overall, these results suggest that correlation models combined with on-site elemental analysis with portable analyzers may enable a screening method for leachable As in coal ash. For the second objective, we utilized nanoscale 2-D imaging (30-50 nm spot size) with the Hard X-ray Nanoprobe (HXN) in combination with microprobe X-ray capabilities (~5 µm resolution) to determine As and Se elemental associations in fly ash particles. Speciation of As and Se was also measured at the nano- to microscale with X-ray absorption spectroscopy. The enhanced resolution of HXN showed As and Se that were diffusely located around or comingled with Ca- and Fe-rich particles. The results also showed nanoparticles of Se attached to the surface of fly ash grains. Overall, a comparison of As and Se species across scales highlights the heterogeneity and complexity of chemical associations for these trace elements of concern in coal fly ash. For the final objective, we developed a predictive model for major element composition of coal ash in reserve at disposal sites of major U.S. coal fired power plants. This model was constructed from coal purchase records of 705 power stations from 1973-2022 and was trained on coal ash composition data showing that coal ash elemental composition is strongly associated with the source of feedstock coal. The model showed regional shifts in the major element contents of ash produced by power plants in the last 50 years, particularly for calcium and iron (expressed as %CaO and %Fe 2 O 3 ), as coal-fired power stations changed their source of coal over this time frame. Our approach enables an estimation of coal ash chemical composition that is stored in waste impoundments at individual power stations. Such information can help delineate the regional market potential for material applications that would utilize coal ash harvested from disposal sites across the U.S.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Artificial Intelligence in Nuclear Safeguards; Evaluating Safeguards and Security Risks and Benefits for Advanced and Small Modular Reactor Deployments

Rapidly growing interest in advanced and small modular reactor (A/SMR) technologies presents challenges as well as opportunities for implementing international safeguards and security. A/SMR deployments are expected to be more numerous, more geographically dispersed, and more varied in their designs, placing new demands on the data systems and analytical tools used to support oversight (Alberti et al., 2023; Canadian Nuclear Safety Commission et al., 2024). Because of this variability, the importance and reliance on data systems for A/SMR deployments is expected to be higher than for previous reactor generations. Artificial Intelligence and Machine Learning (AI/ML) offer potential capabilities to address the high variability inherent in A/SMR technology. The beneficiaries of AI-assisted tools include facility operators, government regulators, IAEA inspectors, and A/SMR vendors. This report analyzes how AI/ML-assisted technologies can strengthen the implementation of IAEA safeguards and security measures. It also identifies AI-assisted tools to strengthen operator, facility, and regulator knowledge management practices and examines the potential risks AI/ML-based tools may introduce to IAEA safeguards and security efforts. It concludes with a set of hypothetical, standards-style requirements for AI/ML systems used in safeguards contexts, grounded in an inspector-centric view of system verification. Despite the potential benefits of AI/ML systems, understanding potential intentional and unintentional failure modes is critical for ensuring adequate protection of nuclear materials and facilities. Unique features of A/SMRs including sealed cores, remote and novel paradigms of operation, off-site reactor fabrication, novel fuel forms, and varied refueling requirements, introduce challenges for traditional safeguards technological approaches (Pensado et al., 2024; Federation of American Scientists, 2025). AI/ML systems deployed to address these challenges may introduce new risks requiring systematic evaluation rooted in both AI-specific risk frameworks, such as the NIST AI Risk Management Framework (NIST AI RMF), and established cyber risk management standards such as NIST SP 800-30 (National Institute of Standards and Technology [NIST], 2023; NIST, 2012).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Editorial: Applications of spectroscopy and chemometrics in nuclear materials analysis

Optical analysis techniques, including spectroscopy and image analysis, have many advantages when applied to the study of nuclear materials. They require small sample sizes, can be performed remotely, and can be proceduralized through consistent practice. Most importantly, they provide a wealth of information by generating multivariate data. For example, ultraviolet–visible–near-infrared absorbance spectroscopy of actinides in aqueous and organic solutions is dependent on the oxidation state, anionic complexation, and temperature. These variables are important for solution-based separation processes, and sensitivity to these factors, combined with online monitoring, can drive the efficiency and control of these processes. The morphology and chemical composition of actinide particles can also provide a vital clue to the mechanisms by which the particles were formed, providing forensic information on the origins of the particles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

An efficient hybrid downscaling framework to estimate high-resolution river hydrodynamics

Flow depth and velocity are the most important hydrodynamic variables that govern various river functions, including water resources, navigation, sediment transport, and biogeochemical cycling. Existing high-resolution flow depth simulations rely on either computationally expensive river hydrodynamic models (RHMs) or data-driven models with formidable training costs, whereas data-driven modeling of flow velocity has rarely been explored. Here, using the hybrid Low-fidelity, Spatial analysis, and Gaussian process learning (LSG) model, we developed a downscaling approach to construct high-resolution flow depth and velocity from a two-dimensional (2-D) RHM simulation at coarse resolution. The LSG models were trained and tested in an urban watershed in Houston using two different hurricane-driven flood events. The high-resolution (as fine as 30 m resolution) and low-resolution (mostly 1000 m resolution) meshes include 664 724 and 14 536 grid cells, respectively. The results showed that through downscaling, the simulation errors were reduced to less than one-fourth and one-third of the errors of the low-resolution 2-D RHM for flow depth and velocity, respectively. Our analysis further revealed that the dominant uncertainty sources of the downscaled hydrodynamics are different, with flow velocity dominated by the dimensionality reduction error, which we reduced by using a regionalized training procedure. The downscaling approach achieves an 84-fold acceleration in computational time compared to the high-resolution 2-D RHM, making high-fidelity ensemble flood modeling feasible. More importantly, the developed method provides an opportunity to couple large-scale hydrodynamical processes with local physical, chemical, and biological processes in river models.

Tan, Zeli [Pacific Northwest National Laboratory (↗
DOE OSTI2025

Supersaturation, Nucleation, and Phase Separation of Mesoscopic Systems

Supersaturation, nucleation, and phase separation are ubiquitous phenomena of great interest in both science and industry. However, a unified, quantitative understanding of these phenomena has yet to be achieved for mesoscopic systems. Here, we present a set of general equations that determine the monomer saturation degree, the size distribution, and the free energy of mesoscopic systems, as well as their phase-transition conditions. These equations reveal that, under supersaturation, the largest cluster size (LCS) is an important state variable; the supersaturation degree decreases with the LCS, approaching unity in the macroscopic limit. We identify the critical supersaturation, at which the nuclei undergo the phase transition to form large crystals. Below this critical supersaturation, the nucleus size distribution is either a unimodal function or a monotonically decreasing function of size, depending on the system and temperature. We also predict the most probable nucleus size and the direction of spontaneous changes of the LCS. Our theory provides a unified, quantitative explanation of the nucleus-size-distribution across six different systems, including nanoparticles and biological condensates. This work serves as a general theoretical framework useful for understanding and designing nucleation and phase transitions of mesoscopic systems.

Kang, Jingyu↗
DOE OSTI2025

Nonlocal, Pattern-Aware Response and Feedback Framework for Regional Climate Response

We devise a pattern-aware feedback framework for representing the forced climate response using a suite of Green’s function experiments with solar radiation perturbations. By considering the column energy balance, a comprehensive linear response function (CLRF) for important climate variables and feedback quantities such as moist static energy, sea surface temperature, albedo, cloud optical depth, and lapse rate is learned from Green’s function data. The learned CLRF delineates the effects of the energy diffusion in both the ocean and atmosphere and the pattern-aware feedbacks from the aforementioned radiatively active processes. The CLRF can then be decomposed into forcing–response mode pairs, which are in turn used to construct a reduced-order model describing the dominant dynamics of climate responses. These mode pairs capture nonlocal effects and teleconnections in the climate and thus make the reduced-order model apt for capturing regional features of climate response. A key observation is that the CLRF captures the polar-amplified response as the most excitable mode of the climate system, and this mode is explainable in the data-learned pattern-aware feedback framework. The reduced-order model can be used for predicting the response for a given forcing and for reconstructing the forcing from a given response; we demonstrate these capabilities for multiple independent forcing scenarios.

Feedback↗
DOE OSTI2026

VIC-Global Parameter Dataset Sensitivity with the Variable Infiltration Capacity Model: Evaluating the importance of dynamic land surface parameters when using the VIC-Global parameter dataset

Accurate prediction of runoff is essential to water resources management, flood risk assessment, and ecosystem protection. However, many hydrological models still have relatively substantial limitations when representing the influence of land use and land cover (LULC) on runoff generation and routing. Changes in LULC, such as deforestation, urban expansion, agricultural intensification, and wetland loss, have been shown to alter the water balance at the land surface through fundamental hydrologic processes (e.g., interception, infiltration, evapotranspiration, and soil storage). However, it remains an open question what the exact magnitude and timing of these impacts are for the spatial and temporal scales commonly used in engineering applications. In this analysis we focus on one aspect of recent LULC change for assessing human impacts, which is urbanization. Specifically we seek to determine the impacts of urbanization on the magnitude and timing of surface runoff and baseflow in HUC-12 basins in Clark County, Nevada which has experienced rapid urbanization. We use the Variable Infiltration Capacity (VIC) hydrology model with a widely used off-the-shelf dataset of land surface parameters, VIC-Global, both of which have been commonly used in the past for water and energy balance modeling for large scale hydrologic studies. We examine two scenarios where the first scenario removes all urbanized land cover and parameterizes those areas of the basins as barren or open shrubland. The second scenario tests the opposite case where all areas of the basins are classified as urban regardless of their present classification. The results from the VIC model show there is a low sensitivity for daily surface runoff between scenarios. The daily baseflow values indicate similar low sensitivity to the classification change during specific periods, but then have substantial differences during other period when large precipitation events are occurring. This is likely due to the assumed parameter values for the urban land cover classification made by the VIC-Global dataset. Using a static land cover parameterization is reasonable for large domain hydrology models that are being used for near-term planning horizons (<30 years). However, longer planning horizons where feedbacks between the atmosphere and land surface are important, especially in transient climate situations, considerations for how to update land surface parameters should be incorporated.

42 ENGINEERING↗
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↗