SEARCH · Engineering Papers

Results for “root”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Effect of the integral length scales of turbulent inflows on wind turbine loads

As wind turbines become larger, the fluctuations in the inflow become increasingly influential in the turbine structural loading. These fluctuations are characterized by the integral length scale, which represents the average size of the largest energy-containing turbulent eddies. Current design standards neglect the varying integral length scales that characterize inflows of wind turbines in operation. Using large-eddy simulations, we generate turbulent inflows of varying integral length scales and quantify the loads of the IEA 15-MW reference wind turbine. Results illustrate that the impact of turbulence on rotor and tower loads is up to 10 times greater than the impact of the mean shear profile. Increasing integral length scales from 0.3x to 0.5x the rotor diameter reduces blade root flapwise moments and rotor and tower loads. Increasing integral length scales from 0.5x to 1.7x the rotor diameter increases the rotor aerodynamic thrust force and the blade root flapwise shear loads and decreases the tower base torsional moment and the tilting and yawing rotor aerodynamic moments. Additionally, turbulence intensity has a greater impact on wind turbine loads than integral length scales. In conclusion, findings indicate that design standards should consider varying integral length scales for accurate wind turbine loading characterization in turbulent inflow conditions.

17 WIND ENERGY↗

DOE OSTI2021

A decreasing carbon allocation to belowground autotrophic respiration in global forest ecosystems

Belowground autotrophic respiration (RAsoil) depends on carbohydrates from photosynthesis flowing to roots and rhizospheres, and is one of the most important but least understood components in forest carbon cycling. Carbon allocation plays an important role in forest carbon cycling and reflects forest adaptation to changing environmental conditions. However, carbon allocation to RAsoil has not been fully examined at the global scale. To fill this knowledge gap, first, the spatio-temporal patterns of RAsoil from 1981 to 2017 were predicted by a Random Forest (RF) algorithm using the most updated Global Soil Respiration Database (v5) with global environmental variables; second, carbon allocation from photosynthesis to RAsoil (CAB), was calculated as the ratio of RAsoil to gross primary production; and its temporal and spatial patterns were assessed in global forest ecosystems. . Globally, mean RAsoil from forests was 8.9 ± 0.08 Pg C yr-1 (mean ± standard deviation) from 1981 to 2017 with strong spatial variabilities. Temporally, RAsoil increased at a rate of 0.0059 Pg C yr-2, paralleling broader soil respiration changes and indicating increasing carbon respired by roots. Mean CAB was 0.243 ± 0.016 and decreased over time. The temporal trend of CAB varied greatly in space, reflecting uneven responses of CAB to environmental changes. This study is the first attempt to predict global CAB and analyze its temporal and spatial patterns. With the linkage of carbon use efficiency, the developed CAB offers an completely independent approach to quantify global aboveground autotropic respiration spatially and temporally, which could provide crucial insights into carbon flux partition and global carbon cycling under climate change.

Tang, Xiaolu↗

DOE OSTI2021

Belowground allocation and dynamics of recently fixed plant carbon in a California annual grassland

Plant roots and the organisms that surround them are a primary source for stabilized soil organic carbon (SOC). While grassland soils have a large capacity to store organic carbon (C), few field-based studies have quantified the amount of plant-fixed C that moves into soil and persists belowground over multiple years. Yet this characteristic of the soil C cycle is critical to C storage, soil water holding capacity, nutrient provisions, and the management of soil health. Here, we tracked the fate of plant-fixed C following a five-day 13 CO 2 labeling of a Northern California annual grassland, measuring C pools starting at the end of the labeling period, at three days, four weeks, six months, one year, and two years. Soil organic carbon was fractionated using a density-based approach to separate the free-light fraction (FLF), occluded-light fraction (OLF), and heavy fraction (HF). Using isotope ratio mass spectrometry, we measured 13 C enrichment and total C content for plant shoots, roots, soil, soil dissolved organic carbon (DOC), and the FLF, OLF, and HF. The chemical nature of C in the HF was further analyzed by solid state 13 C nuclear magnetic resonance (NMR) spectroscopy.

13C-NMR↗

DOE OSTI2023

Enzymatically hydrolyzed fluorescence-based chemical probe enables in situ mapping of chitinase activity in the rhizosphere

Rhizosphere microbiomes and root exudates play a pivotal biochemical role in helping to catalyze chitin catabolism. Chitin is a recalcitrant and ubiquitous soil biopolymer, estimated to be the second most abundant organic soil biopolymer on Earth. Despite its abundance, role as a source of C and N in soil, and importance to ecosystem function, the biochemical mechanisms controlling chitin fate in the rhizosphere are elusive and poorly understood. To enable spatial mapping of chitinase activity in the rhizosphere, we designed and synthesized an enzymatically activated fluorogenic substrate, chitotriose-TokyoGreen (chitotriose-TG), by incorporating a fluorescein derivative (TG) onto the trimeric unit of chitin. This non-fluorescent substrate is selectively hydrolyzed by chitinase to release TG and yield a high fluorescence signal, which can be used to spatially image and measure chitinase activity in the rhizosphere. To demonstrate the application of this technique, we grew switchgrass (Panicum virgatum) in rhizoboxes amended with a horizontal layer of chitin. We extracted mobile proteins from the rhizobox using a nitrocellulose membrane blotting technique which offer non-destructive enzyme extraction while preserving the 2D spatial position of the enzymes. We then subjected these membranes to the synthesized chitotriose-TG stain to spatially visualize the distribution of chitinase activity within the rhizosphere. Furthermore, we observed increased chitinase activity near plant roots and higher activity within the soil zone enriched in chitin, showing an adaptive response of chitinase activity with spatial focusing in areas of higher chitin abundance. Thus, the enzyme extraction and visualization strategy we describe here can help enlighten efforts to better understand spatial controls on chitin breakdown in rhizosphere, further elucidating the role of chitin as a C and N source in these systems.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Aligning theoretical and empirical representations of soil carbon-to-nitrogen stoichiometry with process-based terrestrial biogeochemistry models

Soil carbon-nitrogen (C:N) stoichiometry acts as a control over decomposition and soil organic matter formation and loss, making it a key soil property for understanding ecosystem dynamics and projected ecosystems responses to global environmental change. However, the controls of soil C:N and how they respond to increasing pressures from global change agents are not fully understood. The “foundational” controls on soil C:N, namely plant and microbial C:N, have been used to predict soil C:N, but fail to accurately simulate all ecosystems and may be insufficient for predictions under global environmental change. Here we present an “emerging” representation of controls of soil C:N that includes plant-microbe-mineral feedbacks that have been shown to regulate soil C:N. We argue that including representation of these emerging drivers in process-based terrestrial biogeochemistry models, which include biological N fixation, mycorrhizae, priming, root exudation of organic acids, and mineralogy (including soil texture, mineral composition, and aggregation), will improve mechanistic representation of soil C:N and associated processes. Such improvements will produce models that will better simulate a variety of ecological states and predict soil C:N when global changes modify plant-microbe-mineral interactions. Here, we align our empirical understanding of controls of soil C:N with those controls represented in models, identifying contexts where emerging drivers might be particularly important to represent (e.g., priming and root exudation in nutrient-limited conditions) and areas of future work. Additionally, we show that implementing emerging drivers of soil C:N results in different simulated outcomes at steady state and in response to elevated atmospheric CO2. Our review and preliminary simulations support the need to incorporate emerging drivers of soil C:N into process-based terrestrial biogeochemistry models, allowing for both theoretical exploration of mechanisms and potentially more accurate predictions of land biogeochemical responses to global change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Shrub and sedge rhizosphere communities display distinct affinities toward exudates and soil organic matter degradation: a quantitative stable isotope probing analysis

Warming temperatures are accelerating permafrost thaw and changing tundra vegetation, where woody shrubs are displacing sedges. Shrubs, such as Betula nana, and sedges, such as Eriophorum vaginatum, exhibit distinct life strategies including unique root-associated, or rhizosphere microbial communities. As permafrost thaws it unlocks previously unavailable carbon and nutrient sources resulting in deeper roots and a translocation of rhizosphere communities. Because permafrost microbial communities contain lower diversity and biomass than rhizosphere communities, the coalescence of rhizosphere and permafrost microbial communities could alter soil organic matter (SOM) degradation rates and increase greenhouse gas emissions. To identify metabolic strategies across distinct rhizosphere and permafrost microbial communities we conducted an isotope tracing incubation experiment. We inoculated thawed permafrost with shrub and sedge rhizosphere communities while adding exudates or water daily and compared this to an uninoculated control. After 46 days, we spiked samples with 18O enriched water or 13C enriched exudates and measured isotope incorporation into microbial DNA with quantitative stable isotope probing (qSIP). Our results indicate that exudate additions had little effect on uninoculated permafrost communities but the addition of exudates and rhizosphere inoculants had a compounding effect on respiration rates. We found that soils inoculated with shrub rhizosphere communities contained a mixture of exudate and SOM degraders while soils inoculated with sedge rhizosphere communities contained mainly SOM degraders. Finally, we found that individual microbial taxa exhibited maximum growth rates under specific combinations of microbial inoculant communities and exudate addition treatments. Our results reveal that microbial niches are strongly influenced by substrate preferences and community context, and suggest that a reduction in sedges and an expansion of shrubs may provide a mechanism by which permafrost carbon losses are mitigated through corresponding shifts in microbial communities and their substrate preferences.

Schaefer, Sean R [Natural Resources and Earth Syst↗

DOE OSTI2023

Toward accurate prediction of partial-penetration laser weld performance informed by three-dimensional characterization – Part II: μCT based finite element simulations

The mechanical behavior of partial-penetration laser welds exhibits significant variability in engineering quantities such as strength and apparent ductility. Understanding the root cause of this variability is important when using such welds in engineering designs. In Part II of this work, we develop finite element simulations with geometry derived from micro-computed tomography (μCT) scans of partial-penetration 304L stainless steel laser welds that were analyzed in Part I. We use these models to study the effects of the welds’ small-scale geometry, including porosity and weld depth variability, on the structural performance metrics of weld ductility and strength under quasi-static tensile loading. We show that this small-scale geometry is the primary cause of the observed variability for these mechanical response quantities. Additionally, we explore the sensitivity of model results to the conversion of the μCT data to discretized model geometry using different segmentation algorithms, and to the effect of small-scale geometry simplifications for pore shape and weld root texture. The modeling approach outlined and results of this work may be applicable to other material systems with small-scale geometric features and defects, such as additively manufactured materials.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Computing the viscous effect in early-time drop impact dynamics

The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. ( Phys. Rev. Lett. , vol. 112, no. 13, 2014 b , p. 134501) show that the liquid–air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier–Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square-root scaling relationship. We also predict the effect of interfacial tension at the liquid–gas interface on the drop impact, showing that it mediates the lift-off behaviour.

42 ENGINEERING↗

DOE OSTI2024

Soil Carbon Stocks Not Linked to Aboveground Litter Input and Chemistry of Old-Growth Forest and Adjacent Prairie

The long-standing assumption that aboveground plant litter inputs have a substantial influence on soil organic carbon storage (SOC) and dynamics has been challenged by a new paradigm for SOC formation and persistence. We tested the importance of plant litter chemistry on SOC storage, distribution, composition, and age by comparing two highly contrasting ecosystems: an old-growth coast redwood (Sequoia sempervirens) forest, with highly aromatic litter, and an adjacent coastal prairie, with more easily decomposed litter. We hypothesized that if plant litter chemistry was the primary driver, redwood would store more and older SOC that was less microbially processed than prairie. Total soil carbon stocks to 110 cm depth were higher in prairie (35 kg C m —2 ) than redwood (28 kg C m —2 ). Radiocarbon values indicated shorter SOC residence times in redwood than prairie throughout the profile. Higher amounts of pyrogenic carbon and a higher degree of microbial processing of SOC appear to be instrumental for soil carbon storage and persistence in prairie, while differences in fine-root carbon inputs likely contribute to younger SOC in redwood. We conclude that at these sites fire residues, root inputs, and soil properties influence soil carbon dynamics to a greater degree than the properties of aboveground litter.

13C-NMR spectroscopy↗

DOE OSTI2021

Impacts of Anisotropic Porosity on Heat Transfer and Off-Gassing during Biomass Pyrolysis

The pore structure of biogenic materials imbues the ability to deliver water and nutrients through a plant from root to leaf. This anisotropic pore granularity can also play a significant role in processes such as biomass pyrolysis that are used to convert these materials into useful products like heat, fuel, and chemicals. Evolutions in modeling of biomass pyrolysis as well as imaging of pore structures allow for further insights into the concerted physics of phase change-induced off-gassing, heat transfer, and chemical reactions. In this work, we report a biomass single particle model which incorporates these physics to explore the impact of implementing anisotropic permeability and diffusivity on the conversion time and yields predicted for pyrolysis of oak and pine particles. Simulation results showed that anisotropic permeability impacts predicted conversion time more than 2 times when the Biot number is above 0.1 and pyrolysis numbers (Py 1 , Py 2 ) are less than 20. Pore structure significantly impacts predicted pyrolytic conversion time (>8 times) when the Biot number is above 1 and the pyrolysis number is below 1, i.e., the “conduction controlled” regime. Therefore, these nondimensional numbers reflect that when internal heat conduction limits pyrolysis performance, internal pyrolysis off-gassing further retards effective heat transfer rates as a closely coupled phenomenon. Overall, this study highlights physically meaningful opportunities to improve particle-scale pyrolysis modeling and experimental validation relevant to a variety of feedstock identities and preparations, guiding the future design of pyrolyzers for efficient biomass conversion.

09 BIOMASS FUELS↗

DOE OSTI2022

Understanding the Origin of the Nonpassivating Behavior of Si-Based Anodes during the Initial Cycles

In this contribution, we combined electrochemical cycling and X-ray photoelectron spectroscopy (XPS) to understand the nonpassivating behavior of the solid electrolyte interphase (SEI) on Si anodes during the first cycles. Based on galvanostatic measurements, we show that the irreversible capacity loss is reduced after the first cycle, and it stays almost constant from the second cycle onwards. XPS was used to determine the root causes of the Coulombic inefficiency, showing that the rate of decomposition of the organic solvents strongly decreased after the first cycle, whereas the rate of salt decomposition is almost unchanged between cycles. We determine that the inhibition of the decomposition reaction of the organic solvent is responsible for the lower Coulombic loss during the second electrochemical cycle in comparison to the first, whereas the nonpassivating behavior toward the salt decomposition is one of the main causes of capacity loss upon cycling. Here, we further revisit the role of cracking in contributing to capacity loss. Whereas high volumetric expansion remains an issue plaguing the performance of Si anodes, our chronoamperometry studies reveal that the SEI formed on Si anodes does not passivate even when the electrode is fully expanded, and no additional surface is exposed. Overall, our work establishes the need to address the chemical and electrochemical instability of the SEI on the Si anode in addition to the more notorious issue of cracking.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Flash NanoPrecipitation as an Agrochemical Nanocarrier Formulation Platform: Phloem Uptake and Translocation after Foliar Administration

The increasing severity of pathogenic and environmental stressors that negatively affect plant health has led to interest in developing next-generation agrochemical delivery systems capable of precisely transporting active agents to specific sites within plants. In this work, we adapt Flash NanoPrecipitation (FNP), a scalable nanocarrier (NC) formulation technology used in the pharmaceutical industry, to prepare organic core–shell NCs and study their efficacy as foliar or root delivery vehicles. NCs ranging in diameter from 55 to 200 nm, with surface zeta potentials from -40 to +40 mV, and with seven different shell material properties were prepared and studied. Shell materials included synthetic polymers poly(acrylic acid), poly(ethylene glycol), and poly(2-(dimethylamino)ethyl methacrylate), naturally occurring compounds fish gelatin and soybean lecithin, and semisynthetic hydroxypropyl methylcellulose acetate succinate (HPMCAS). NC cores contained a gadolinium tracer for tracking by mass spectrometry, a fluorescent dye for tracking by confocal microscopy, and model hydrophobic compounds (alpha tocopherol acetate and polystyrene) that could be replaced by agrochemical payloads in subsequent applications. After foliar application onto tomato plants with Silwet L-77 surfactant, internalization efficiencies of up to 85% and NC translocation efficiencies of up to 32% were observed. Significant NC trafficking to the stem and roots suggests a high degree of phloem loading for some of these formulations. Results were corroborated by confocal microscopy and synchrotron X-ray fluorescence mapping. NCs stabilized by cellulosic HPMCAS exhibited the highest degree of translocation, followed by formulations with a significant surface charge. The results from this work indicate that biocompatible materials like HPMCAS are promising agrochemical delivery vehicles in an industrially viable pharmaceutical nanoformulation process (FNP) and shed light on the optimal properties of organic NCs for efficient foliar uptake, translocation, and delivery.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Persistence and Plasticity in Conifer Water–Use Strategies

The selective use of seasonal precipitation by vegetation is critical to understanding the residence time and flow path of water in watersheds, yet there are limited datasets to test how climate alters these dynamics. Here, we use measurements of the seasonal cycle of tree ring δ 18 O for two widespread conifer species in the Rocky Mountains of North America to provide a multi-decadal depiction of the seasonal origins of forest water use. The results show that while the conifer tree stands had a dominant preference for use of snowmelt, there were multi-annual periods over the last four decades when use of summer precipitation was preferential. Utilization of summer rain emerged during years with increased snowfall and tree growth, suggesting that summer rain enhanced the transpiration stream only during the periods of highest water use. We hypothesize this could be explained through shallowing of the root profile during wetter periods and/or through the influence of changing water table depths on the residence time of summer precipitation in the root zone. We suggest the tree ring proxy approach used here could be applied in other watersheds to provide critical insight into the temporal dynamics of plant water use that could not be inferred from short measurement campaigns. Furthermore, these data on the seasonal origins of forest water are critical for understanding forest vulnerability to drought, the processes that affect precipitation pathways and residence time in watersheds and the interpretation of tree ring proxy data.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2019

Spatial Predictions and Associated Uncertainty of Annual Soil Respiration at the Global Scale

Abstract Soil respiration (Rs), the soil‐to‐atmosphere CO 2 flux produced by microbes and plant roots, is a critical but uncertain component of the global carbon cycle. Our current understanding of the variability and dynamics is limited by the coarse spatial resolution of existing estimates. We predicted annual Rs and associated uncertainty across the world at 1‐km resolution using a quantile regression forest algorithm trained with observations from the global Soil Respiration Database spanning from 1961 to 2011. This model yielded a global annual Rs estimate of 87.9 Pg C/year with an associated global uncertainty of 18.6 (mean absolute error) and 40.4 (root mean square error) Pg C/year. The estimated annual heterotrophic respiration (Rh), derived from empirical relationships with Rs, was 49.7 Pg C/year over the same period. Predicted Rs rates and associated uncertainty varied widely across vegetation types, with the greatest predicted rates of Rs in evergreen broadleaf forests (accounting for 20.9% of global Rs). The greatest prediction uncertainties were in northern latitudes and arid to semiarid ecosystems, suggesting that these areas should be targeted in future measurement campaigns. This study provides predictions of Rs (and associated prediction uncertainty) at unprecedentedly high spatial resolution across the globe that could help constrain local‐to‐global process‐based models. Furthermore, it provides insights into the large variability of Rs and Rh across vegetation classes and identifies regions and vegetation types with poor model performance that should be prioritized for future data collection.

Warner, D. L.↗

DOE OSTI2020

Parameterizing Perennial Bioenergy Crops in Version 5 of the Community Land Model Based on Site‐Level Observations in the Central Midwestern United States

Abstract With projected expansion of biofuel production at a global scale, there is a pressing need to develop adequate representation of bioenergy crops in land surface models to help effectively quantify the biogeophysical and biogeochemical effects of its associated land use changes. This study implements two new perennial bioenergy crops, Miscanthus and switchgrass, into the Community Land Model Version 5 based on site‐level observations from the midwestern United States by modifying parameters associated with photosynthesis, phenology, allocation, decomposition, and carbon cost of nitrogen uptake and integrating concomitantly land management practices. Sensitivity analyses indicate that carbon and energy fluxes of the perennial crops are most sensitive to photosynthesis and phenology parameters. Validation of simulated fluxes against site‐level measurements demonstrates that the model is capable of capturing the overall patterns of energy and carbon fluxes, as well as physiological transitions from leaf emergence to senescence. Compared to annual crops, perennial crops feature longer growing season, greater leaf areas, and higher productivity, leading to increased transpiration, lower annual runoff, and larger carbon uptake. The model simulations suggest that with higher CO 2 assimilation rates and lower demands for nutrients and water, high‐yielding perennial crops are promising alternatives of bioenergy feedstocks compared to traditional annual crops not only for mitigating climate change but also for environmental conservation purposes by reducing fertilizer application and therefore alleviating surface‐ and ground‐water contaminations. Although the local‐scale simulations shed light on potential benefits of using perennial grasses as bioenergy feedstocks, quantifying consequences of their plantations at larger scales warrants additional investigation.

Cheng, Yanyan↗

DOE OSTI2021

Topographical Controls on Hillslope-Scale Hydrology Drive Shrub Distributions on the Seward Peninsula, Alaska

Observations indicate shrubs are expanding across the Arctic tundra, mainly on hillslopes and primarily in response to climate warming. However, the impact topography exerts on hydrology, nutrient dynamics, and plant growth can make untangling the mechanisms behind shrub expansion difficult. We examined the role topography plays in determining shrub expansion by applying a coupled transect version of a mechanistic ecosystem model (ecosys) in a tundra hillslope site in the Seward Peninsula, Alaska. Modeled biomass of the dominant plant functional types agreed well with field measurements (R 2 = 0.89) and accurately represented shrub expansion over the past 30 years inferred from satellite observations. In the well-drained crest position, canopy water potential and plant nitrogen (N) uptake was modeled to be low from plant and microbial water stress. Intermediate soil water content in the mid-slope position enhanced mineralization and plant N uptake, increasing shrub biomass. The deciduous shrub growth in the mid-slope position was further enhanced by symbiotic N 2 fixation primed by increased root carbon allocation. The gentle slope in the poorly drained lower-slope position resulted in saturated soil conditions that reduced soil O 2 concentrations, leading to lower root O 2 uptake and lower nutrient uptake and plant biomass. Additionally, a simulation that removed topographical interconnectivity between grid cells resulted in (1) a 28% underestimate of mean shrub biomass and (2) over or underestimated shrub productivity at the various hillslope positions. Our results indicate that land models need to account for hillslope-scale coupled surface and subsurface hydrology to accurately predict current plant distributions and future trajectories in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Crustal Structure of the Mesopotamian Plain, East of Iraq

The crustal structure of Iraq was investigated through analyzing teleseismic data from 12 new seismic stations. Three seismic stations are located within the Zagros Fold-Thrust Belt, and nine are within the Mesopotamian Plain. Joint inversion of P wave receiver function and Rayleigh wave dispersion data were employed to resolve S wave velocity structure models beneath each station. Combining these models with available Moho depths within and around the study area reveals that the Moho depth is smoothly increasing from the Arabian platform toward the Zagros Mountains. Further, an exceptional deeper root is observed near the eastern edge of the Mesopotamian Plain where sedimentary pile is the thickest. This root was interpreted as a structure inherited from the successive Mesozoic rifting of the NE Arabian platform enhanced by sedimentary loading and progressive subsidence. A low-velocity uppermost mantle was observed beneath the Arabian Foredeep attesting the existence of a low-strength lithospheric mantle beneath the region southwest from the Zagros deformation front. The weak rheology of the uppermost mantle may have allowed local sinking of the crust and deepening of the Moho boundary due to vertical loading and the Late Miocene lateral contraction.

58 GEOSCIENCES↗

DOE OSTI2025

Baseflow Identification via Explainable AI With Kolmogorov‐Arnold Networks

Abstract Hydrological models often involve constitutive laws that may not be optimal in every application. We propose to replace such laws with the Kolmogorov‐Arnold networks (KANs), a class of neural networks designed to identify symbolic expressions. We demonstrate KAN's potential on the problem of baseflow identification, a notoriously challenging task plagued by significant uncertainty. KAN‐derived functional dependencies of the baseflow components on the aridity index outperform their original counterparts; they demonstrate that water availability, rather than potential evapotranspiration, drives baseflow by constraining actual evapotranspiration under arid conditions. On a test set, they increase the Nash‐Sutcliffe efficiency (NSE) by 65%, decrease the root mean squared error by 29%, and increase the Kling‐Gupta efficiency by 34%. This superior performance is achieved while reducing the number of fitting parameters from three to two. Next, we use data from 378 catchments across the continental United States to refine the water‐balance equation at the mean‐annual scale. The KAN‐derived equations based on the refined water balance outperform both the current aridity index model, with up to a 105% increase in NSE, and the KAN‐derived equations based on the original water balance. While the performance of our model and tree‐based machine learning methods is similar, KANs offer the advantage of simplicity and transparency and require no specific software or computational tools. This case study focuses on the aridity index formulation, but the approach is flexible and transferable to other hydrological processes. Plain Language Summary Equations used in hydrologic model are often suboptimal, resulting in reduced prediction accuracy and efficiency. We implemented Kolmogorov‐Arnold networks (KAN), a machine learning algorithm for deriving symbolic formulations, to estimate groundwater recharge and showed that it outperforms an existing state‐of‐the‐art semi‐empirical formulation. In hydrology, Nash‐Sutcliffe efficiency (NSE), root mean squared error (RMSE), and Kling‐Gupta efficiency (KGE) are commonly used to evaluate model performance. Higher NSE and KGE values indicate better performance, while lower RMSE values are preferable. Our results show that NSE increased by 71%, RMSE decreased by 32%, and KGE improved by 25%. In addition, KAN identifies an optimal functional form and can be used to derive new analytical formulas using the prior knowledge. The KAN‐inspired equation outperformed the original formulation and reduced the fitting parameters. Furthermore, we refined the water‐balance equation at the mean‐annual scale and showed that, based on the new water‐balance equation, KAN can derive new formulations that are superior to the original aridity index formulations (up to 105% increase in NSE) and KAN‐derived equations based on the original water balance. These findings highlight the significant potential of KAN to advance the scientific understanding of a wide range of hydrologic processes. Key Points Kolmogorov‐Arnold networks (KANs) enhance interpretability of machine‐learned hydrological models KAN‐derived symbolic formulations outperform state‐of‐the‐art semi‐empirical aridity indices KAN‐identified functional form yields an analytical index with fewer fitting parameters and improved performance

baseflow↗