Search NASASearch

SEARCH · Search NASA

Results for “Variable Importance”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

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

Leveraging Observations of Untrained Panelists to Screen for Quality of Fresh-Cut Romaine Lettuce

Fresh-cut romaine lettuce’s high perishability challenges ready-to-eat (RTE) salad production. Selecting cultivars less prone to browning and decay is crucial for extending shelf life. Traditional quality evaluation methods using instrumentation and trained panelists are time-consuming and logistically complex. This study investigated the effectiveness of untrained volunteers in assessing fresh-cut romaine lettuce quality. Given that the average consumer in the USA is familiar with the flavor characteristics of romaine lettuce, this study proposed to investigate the value of having untrained volunteers discern the quality of fresh-cut romaine lettuce. Therefore, six romaine lettuce accessions (Green Forest, King Henry, Parris Island Cos, PI 491224, SM13-R2, and Sun Valley) were assessed for sensory quality attributes (browning, green color, decay, and overall quality) and compared with instrumentation analyses (gas composition including O2 and CO2, electrolyte leakage, and color). The results showed significant quality differences (p < 0.05) among the accessions, with some seasonal variability. Very importantly, the consumers’ (n = 159) assessments revealed similar results to those produced by either instrumentation or a trained panel. The consumers provided sensory scores that allowed for the grouping of accessions based on their postharvest quality, which efficiently matched their pedigree relationship. In conclusion, ad hoc consumer panels can be an effective way to characterize the quality of romaine lettuce for RTE salads.

Agriculture

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

SetEnvironment

SAND2025-11727O SetEnvironment software is a utility library module that uses ConfigParserEnhanced, also a library module, to read specifically formatted .ini files for setting the user's environment. It creates a consistent environment with an expected set of variables and module loads, which is important for testing. Environment variable operations include set, append, prepend, unset, remove, and others and uses the Python os.environ methods. The module operations are performed with a wrapper. Systems using the Lmod package can generate a module function, which can also perform module operations specified in the text file. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

McLendon III, William

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

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

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

54 ENVIRONMENTAL SCIENCES

Estimating irrigation water use from remotely sensed evapotranspiration data: Accuracy and uncertainties at field, water right, and regional scales

Irrigated agriculture is the dominant user of water globally, but most water withdrawals are not monitored or reported. As a result, it is largely unknown when, where, and how much water is used for irrigation. Here, we evaluated the ability of remotely sensed evapotranspiration (ET) data, integrated with other datasets, to calculate irrigation water withdrawals and applications in an intensively irrigated portion of the United States. We compared irrigation calculations based on an ensemble of satellite-driven ET models from OpenET with reported groundwater withdrawals from hundreds of farmer irrigation application records and a statewide flowmeter database at three spatial scales (field, water right group, and management area). At the field scale, we found that ET-based calculations of irrigation agreed best with reported irrigation when the OpenET ensemble mean was aggregated to the growing season timescale (bias = 1.6–4.9%, R 2 = 0.53–0.74), and agreement between calculated and reported irrigation was better for multi-year averages than for individual years. At the water right group scale, linking pumping wells to specific irrigated fields was the primary source of uncertainty. At the management area scale, calculated irrigation exhibited similar temporal patterns as flowmeter data but tended to be positively biased with more interannual variability. Disagreement between calculated and reported irrigation was strongly correlated with annual precipitation, and calculated and reported irrigation agreed more closely after statistically adjusting for annual precipitation. The selection of an ET model was also an important consideration, as variability across ET models was larger than the potential impacts of conservation measures employed in the region. From these results, we suggest key practices for working with ET-based irrigation data that include accurately accounting for changes in soil moisture, deep percolation, and runoff; careful verification of irrigated area and well-field linkages; and conducting application-specific evaluations of uncertainty.

59 BASIC BIOLOGICAL SCIENCES

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

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

58 GEOSCIENCES

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

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

54 ENVIRONMENTAL SCIENCES

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

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

09 BIOMASS FUELS

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

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

Data-driven discovery

Evaluation of average leaf inclination angle quantified by indirect optical instruments in crop fields

Average leaf inclination angle ($\overline{θ}$ L ) is an important canopy structure variable that influences light regime, photosynthesis, and evapotranspiration of plants. $\overline{θ}$ L can be measured through direct methods (e.g., protractor), which are labor-intensive and time-consuming, or through indirect optical instruments, which are more efficient than the direct methods. However, uncertainties of different indirect optical instruments for quantifying $\overline{θ}$ L remain largely unquantified. In this study, we evaluated and compared the performances of three major indirect optical instruments: (1) LAI-2200, (2) 30°-tilted camera, and (3) digital hemispherical photography (DHP), in different crop fields over a growing season, benchmarked with direct measurements. LAI-2200 and 30°-tilted camera showed higher agreement with direct $\overline{θ}$ measurements (R 2 = 0.54, RMSE = 7.37°; R 2 = 0.58, RMSE = 8.08°) than DHP (R 2 = 0.14, RMSE = 13.96°). Different performances of indirect optical instruments could be attributed to the accuracy of gap fraction measurement and the performance of the $\overline{θ}$ L quantification algorithms. When using the LAI-2200 algorithm, larger gap fraction gradients over view zenith angles led to larger $\overline{θ}$ L values, and smaller gap fraction gradients led to smaller $\overline{θ}$ L values. Such error propagation was larger in sparse canopy than in dense canopy. The Wilson G function of the LAI-2200 algorithm performed better in estimating $\overline{θ}$ L than the G function based on the ellipsoidal LAD function used by the CAN_EYE algorithm. We also proposed a modification of the LAI-2200 algorithm, which further improved the performance of LAI-2200 and 30°-tilted cameras in estimating $\overline{θ}$ L . We envision that the low-cost 30°-tilted cameras provide a promising sensor solution to continuously monitor canopy structure for various ecosystems.

30°-tilted camera

Machine learning of factors for improving oyster hatchery production

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

Vishwakarma, Srishti [Oak Ridge National Laborator

Modeling Oxygen Partial Pressure in Solid Oxide Electrolysis Cells: The Microstructure Effect

Oxygen partial pressure is an important thermodynamic state variable that affects both the performance and degradation of solid oxide electrolysis cells. In this work, a 3D model developed from Virkar’s 1D model has been applied to reconstructed and synthetic microstructures of Ni-YSZ-GDC-LSCF cell. The effect of the microstructures, including the thickness of the YSZ and GDC layer, and the compositions of the hydrogen and oxygen electrode, on the distribution of oxygen partial pressure has been investigated. The results show that the maximum oxygen partial pressure may occur on the interface between oxygen electrode and GDC layer or between GDC layer and YSZ layer depending on the rate of oxygen ion exchange between GDC and YSZ. A thicker GDC layer lowers the maximum oxygen partial pressure in the cell, while a thicker YSZ layer lowers the maximum oxygen partial pressure in the hydrogen electrode. In addition, the Ni:YSZ ratio and porosity also affects the maximum partial pressure. These findings provide insights on mitigating degradation in solid oxide electrolysis cell by tuning the microstructures.

Lei, Yinkai

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

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

54 ENVIRONMENTAL SCIENCES

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