Engineering PapersSearch

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 55 records · Page 3

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

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

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

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

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)

Separate the Role of Southern and Northern Extra‐Tropical Pacific in Tropical Pacific Climate Variability

Abstract Observational and modeling studies have elucidated the influential role played by the southern and northern extratropical Pacific (SEP and NEP) forcing in shaping dynamics of tropical Pacific climate variability. However, the relative importance of the NEP and SEP and the timescale on which they impact the tropics remain unclear. Using a linear inverse model (LIM) that selectively incorporates or excludes tropical‐extratropical coupling, we find a reduction in tropical interannual variability (∼40%) and low‐frequency (sub‐decadal to decadal) variability in the southeastern tropical Pacific region (∼70%) in the absence of SEP. Conversely, the absence of NEP yields no significant impact on tropical interannual variability but markedly diminishes low‐frequency variability in the central tropical Pacific region (∼70%). LIM and statistic diagnostics on CMIP6 models show the low‐frequency to total variability ratio in the tropical Pacific depending on their NEP and SEP representation. Models with more (less) low‐frequency power tend to show stronger NEP (SEP) dynamics.

Geology

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

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

Subdecadal Holocene Warm-Season Temperature Variability in Central Europe Recorded by Biochemical Varves

Paleoclimate data provide important information about the character of natural climate variability. However, records with sufficient length and resolution to resolve high-frequency (decadal-scale) variability across the Holocene are scarce. We present a 10,800-year reconstruction of spring and summer temperature at three-year resolution based on biochemical varves from Lake $\dot{Z}abi$$\acute{n}skie$, Poland. The reconstruction is based on Ca/Ti ratio, which are significantly correlated with instrumental spring and summer temperature spanning 240 years. Major climate events of the Holocene period are represented in the reconstruction, including the Holocene Thermal Maximum, 8.2 ka Event, Medieval Climate Anomaly, and Little Ice Age. A low-frequency 8,000-year decreasing trend in warm-season temperatures is driven by declining summer insolation. Temperature variability is highest during the early Holocene, likely related to warmer and drier conditions. The rate of warming during the past 90 years is extremely unusual, if not unprecedented for the Holocene, based on our reconstruction.

58 GEOSCIENCES

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

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

No constraint on long-term tropical land carbon-climate feedback uncertainties from interannual variability

Unraveling drivers of the interannual variability of tropical land carbon cycle is critical for understanding land carbon-climate feedbacks. Here we utilize two generations of factorial model experiments to show that interannual variability of tropical land carbon uptake under both present and future climate is consistently dominated by terrestrial water availability variations in Earth system models. The magnitude of this interannual sensitivity of tropical land carbon uptake to water availability variations under future climate shows a large spread across the latest 16 models (2.3 ± 1.5 PgC/yr/Tt H 2 O), which is constrained to 1.3 ± 0.8 PgC/yr/Tt H 2 O using observations and the emergent constraint methodology. However, the long-term tropical land carbon-climate feedback uncertainties in the latest models can no longer be directly constrained by interannual variability compared with previous models, given that additional important processes are not well reflected in interannual variability but could determine long-term land carbon storage. Our results highlight the limited implication of interannual variability for long-term tropical land carbon-climate feedbacks and help isolate remaining uncertainties with respect to water limitations on tropical land carbon sink in Earth system models.

54 ENVIRONMENTAL SCIENCES

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

Multiscale Wildfire Simulation Framework and Remote Sensing

Wildfire as one type of climate extreme events causes huge socioeconomic losses and damages. Large wildfires (i.e., generated pyrocumulonimbus (PyroCb)) can inject tremendous amounts of smoke into the stratosphere, where black carbon and organic carbon aerosols can persist months to years and influence climate by imposing a significant reduction in the radiative forcing like that associated with large volcanic eruptions or proposed via climate interventions such as geoengineering. Both observations and numerical modeling results clearly indicate an increasing trend in wildfire frequency and intensity in many regions during the recent decades with climate change. However, current understanding of wildfire remains largely uncertain owing to limitations of modeling capabilities in representing the multiscale wildfire physics and dynamics and a scarcity of observations constraining important wildfire and environmental variables. This study primarily aims to improve the wildfire simulation capabilities in the state-of-the-art climate model by filling in two major gaps: (1) model resolution is typically too coarse to resolve fine scale processes associated with fires, and (2) chemistry and aerosol processes in fire smoke are poorly represented.

54 ENVIRONMENTAL SCIENCES

Hybrid-Energy-Powered Electrochemical Ocean Alkalinity Enhancement Model: Plant Operation, Cost, and Profitability

Electrochemical ocean alkalinity enhancement is a form of marine carbon dioxide removal, a rapidly growing industry that is powered by efficient onshore or offshore energy sources. As more and larger deployments are being planned, it is important to consider how variable energy sources like tidal energy can impact plant performance and costs. An open-source Python-based generalizable model for electrodialysis-based ocean alkalinity enhancement has been developed that can capture key system-level insights of the electrochemistry, ocean chemistry, acid disposal, and co-product creation of these plants under various conditions. The model additionally accounts for hybrid energy system performance profiles and costs via the National Laboratory of the Rockies’ H2Integrate tool. The model was used to analyze an example theoretical plant deployment in North Admiralty Inlet, including how the plant is impacted by the available energy sources in the region and the scale at which plant costs are covered by the co-products it generates, such as recycled concrete aggregates, without requiring carbon credits. The results show that the example plant could be profitable without carbon credits at commercial scales of 100,000 to 1 million tons of carbon dioxide removal per year, so long as it uses low-cost electricity sources and either sells acid or recovers recycled concrete aggregates with about 1 molar acid concentrations, though more research is needed to confirm these results.

hybrid energy

Mesoscale variations of the chemical and electronic landscape on the surface of the Weyl semimetal Co 3 ⁢Sn 2⁢ S 2 visualized by ARPES and XPS

The multiple crystalline terminations in magnetic Weyl semimetal Co 3 Sn 2 S 2 display distinct topological and trivial surface states, which have successfully been distinguished experimentally. However, a model of pure terminations is known to be inadequate because these surfaces exhibit a high degree of spatial heterogeneity and point disorder. Here, in this study, we perform a spectromicroscopy study of the surface chemistry and surface electronic structure using photoemission measurements in combination with first-principles calculations of core levels. We identify an intermediate region with properties distinct from both the sulfur and tin terminations, and demonstrate that the spectral features in this region can be associated with a disordered termination with a varying density of surface tin vacancies. This work establishes heuristics for identifying variable surface disorder using photoemission, an important prerequisite to experimentally establishing the behavior of momentum-space topological surface features subject to variable surface disorder on a single cleave.

Sreedhar, Sudheer Anand [Univ. of California, Davi

Influence of simple terrain on the spatial variability of a low-level jet and wind farm performance in the AWAKEN field campaign

In wind energy research, scientific challenges are often associated with complex terrain sites, where orography, vegetation, and buildings disrupt flow uniformity. However, even sites characterized as simple terrain can exhibit significant spatial variability in wind speed, particularly during stable boundary layers (SBLs) and low-level jets (LLJs). This study investigates these terrain interactions using both simulations and observations from the American WAKe ExperimeNt (AWAKEN). We employ a multiscale Weather Research and Forecasting (WRF) model simulation, integrating mesoscale forcing in the coarse domains and representing three rows of turbines from the King Plains wind farm as generalized actuator disks (GAD) in the large-eddy simulation (LES) domains. During a nocturnal LLJ event on 3 April 2023, the downstream, wake-affected turbine rows outperformed the upstream, unwaked row by 25 %–51 %. This counterintuitive result arises from terrain-induced streamwise variations in hub-height wind speed of approximately 4 m s −1 over 5 km – equivalent to ∼50 % of the upstream reference speed. This enhancement outweighs the wake-induced reduction in mean wind speed (∼12 %) and global blockage effects reported in the literature (∼1 %–3.4 %). The multiscale simulations capture the intra-farm spatial variability in power performance observed in SCADA data. Terrain-induced vertical displacement of the LLJ, coupled with large wind shear below the jet maximum, drives the substantial streamwise acceleration within the wind farm. These findings underscore the importance of accounting for spatial variability related to terrain, even in simple landscapes, particularly during LLJ conditions. Incorporating such effects into reduced-order modeling frameworks for wind farm design and control could significantly enhance their effectiveness.

17 WIND ENERGY

Unique Temperature Trend Pattern Associated With Internally Driven Global Cooling and Arctic Warming During 1980–2022

Abstract Diagnosing the role of internal variability over recent decades is critically important for both model validation and projections of future warming. Recent research suggests that for 1980–2022 internal variability manifested as Global Cooling and Arctic Warming (i‐GCAW), leading to enhanced Arctic Amplification (AA), and suppressed global warming over this period. Here we show that such an i‐GCAW is rare in CMIP6 large ensembles, but simulations that do produce similar i‐GCAW exhibit a unique and robust internally driven global surface air temperature (SAT) trend pattern. This unique SAT trend pattern features enhanced warming in the Barents and Kara Sea and cooling in the Tropical Eastern Pacific and Southern Ocean. Given that these features are imprinted in the observed record over recent decades, this work suggests that internal variability makes a crucial contribution to the discrepancy between observations and model‐simulated forced SAT trend patterns.

54 ENVIRONMENTAL SCIENCES