Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial heterogeneity”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulation of dilatancy-controlled gas migration processes in saturated bentonite using a coupled multiphase flow and elastoplastic H 2 M model

Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories. This is particularly the case for bentonite, an often-preferred barrier material. Gas flow in preferential pathways is characterized by localization and spontaneous behavior, which is challenging to simulate in numerical models due to strong hydro-mechanical coupling. To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project, this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations. The proposed model utilized hydro-mechanical spatial distributions, namely Young's modulus and gas entry pressure, and elastoplasticity combined with a linear swelling model. A conceptual strain-dependent permeability approach simulated dilatancy-controlled gas flow based on hydro-mechanical coupling. To test the effectiveness of the presented approach, a gas injection test in a compacted, saturated bentonite sample was simulated using the open-source code OpenGeoSys 5.8 and compared with experimental observations. The presented methodology is capable of simulating localized gas flow in preferential pathways. The spatial distributions of Young's modulus and gas entry pressure affect the swelling pressure, relative permeability and, in combination with the strain-dependent permeability model, also the intrinsic permeability.

42 ENGINEERING↗
DOE OSTI2020

An Unmanned Aerial System (UAS) for concurrent measurements of solar-induced chlorophyll fluorescence and hyperspectral reflectance toward improving crop monitoring

Unmanned aerial system (UAS)-based remote sensing can provide high spatial- and temporal-resolution crop monitoring for precision management. Existing crop monitoring UAS primarily use conventional techniques such as multi-spectral broadband vegetation indices (VIs) to track canopy structure/biomass. Recently developed lightweight hyperspectral sensors can track crop physiology and performance via complementary signals such as solar-induced chlorophyll fluorescence (SIF) and photochemical reflectance index (PRI) that can be derived from hyperspectral reflectance. However, few existing UAS can acquire high-quality SIF and hyperspectral reflectance. Here, we designed a novel UAS that simultaneously captures far-red SIF and hyperspectral reflectance, using a single bifurcated fiber and motorized arm to measure downwelling and upwelling irradiance. The UAS was tested over (1) a heterogeneous corn field with spatially varying crop yield potential and (2) a set of soybean and corn plots under differing nutrient treatments. Our UAS can maintain stability under low-to-moderate winds (2–3 m/s) with <0.4 m geolocation accuracy, <0.1 m altitude drift, and <1.5° error of the fiber optic from nadir position during scanning. Weekly seasonal campaigns reveal that SIF outperforms conventional VIs such as the normalized difference vegetation index (NDVI) for distinguishing plots with different crop yield potential, especially after canopy closure when NDVI tends to saturate. The UAS can capture the diurnal dynamics of SIF and PRI for both corn and soybean. At the seasonal scale, SIF acquired from the UAS correlates well with that from a fixed SIF tower system ( R 2 = 0.81), which measured the same target albeit with different footprints, demonstrating the capability of UAS in characterizing the seasonal progression of crop activity. In conclusion, our newly developed UAS provides high-quality SIF and hyperspectral reflectance, facilitating mechanistic understanding of the physiological control on photosynthesis dynamics. Our findings also imply that this UAS could enable extension of fixed tower-based systems to monitor multiple targets from diurnal to seasonal scales and improve monitoring of heterogeneous fields.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Advances in mass spectrometry-enabled multiomics at single-cell resolution

We report biological organisms are multifaceted, intricate systems where slight perturbations can result in extensive changes in gene expression, protein abundance and/or activity, and metabolic flux. These changes occur at different timescales, spatially across cells of heterogeneous origins, and within single-cells. Hence, multimodal measurements at the smallest biological scales are necessary to capture dynamic changes in heterogeneous biological systems. Of the analytical techniques used to measure biomolecules, mass spectrometry (MS) has proven to be a powerful option due to its sensitivity, robustness, and flexibility with regard to the breadth of biomolecules that can be analyzed. Recently, many studies have coupled MS to other analytical techniques with the goal of measuring multiple modalities from the same single-cell. It is with these concepts in mind that we focus this review on MS-enabled multiomic measurements at single-cell or near-single- cell resolution.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Data used in Wainwright, H.M. et al. 2021, “Watershed zonation through hillslope clustering for tractably quantifying above- and belowground watershed heterogeneity and functions”

This data package contains spatial data layers and processing scripts used in Wainwright, H.M. et al. 2021, “Watershed zonation approach for tractably quantifying above-and- belowground watershed heterogeneity and functions”. The purpose of the data and paper is to develop a watershed zonation approach for characterizing watershed organization and function in a tractable manner by applying clustering methods to multiple spatial data layers. The data package contains the geotiff files of spatial data layers, and the processed data values corresponding to the figures in the paper. The Data_description file describe each file in details. The spatial data sets (geotiff) are included in the zip files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Machine Learning Inference of Random Medium Properties

Earth materials are heterogeneous across a range of spatial scales, but the resolvability of small structures is limited by sparse data coverage, noise, bandlimitedness, and other difficulties. In practice, heterogeneities below a certain size cannot be recovered from seismic data except through statistical medium descriptions, which even then can be difficult to uniquely determine. To improve the characterization of such heterogeneities, we develop a novel supervised machine learning (ML) model that provides insight about the recoverability of statistical medium properties from elastic waveform data and succeeds despite cycle-skipping and other challenges well known from elastic waveform inversion. We demonstrate the approach using random media generated by superimposing self-affine random variations on homogeneous and layered background structures. After training on sparsely-recorded, high-frequency waveforms from hundreds of different random medium realizations, we show the ability of our ML model to recover correlation lengths and other statistical properties of interest to near-surface and crustal seismology, among other fields. For frequency passbands and spatial offsets encountered in seismology, Gaussian correlation lengths and the amplitude of the random variations relative to the background model are recovered even in challenging scenarios involving unknown medium parameters, complex crustal structures, and low signal-to-noise ratio. In comparison, von Kármán correlation lengths, which are related to larger-wavelength variations of the medium than Gaussian correlation lengths, are not as well recovered. These results provide one of the first and most systematic investigations of the recoverability of statistical properties of heterogeneities below the resolution limit of deterministic seismic tomography, and suggest practical ML strategies for high-frequency waveform seismology.

58 GEOSCIENCES↗
DOE OSTI2025

An Atomistic Study of Reactivity in Solid-State Electrolyte Interphase Formation for Li/Li7P3S11

Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.

Diffusion↗
DOE OSTI2024

Correlating activities and defects in (photo)electrocatalysts using in-situ multi-modal microscopic imaging

Photo(electro)catalysts use sunlight to drive chemical reactions such as water splitting. A major factor limiting photocatalyst development is physicochemical heterogeneity which leads to spatially dependent reactivity. To link structure and function in such systems, simultaneous probing of the electrochemical environment at microscopic length scales and a broad range of timescales (ns to s) is required. Here, we address this challenge by developing and applying in- situ (optical) microscopies to map and correlate local electrochemical activity, with hole lifetimes, oxygen vacancy concentrations and photoelectrode crystal structure. Using this multi-modal approach, we study prototypical hematite (α-Fe 2 O 3 ) photoelectrodes. We demonstrate that regions of α-Fe 2 O 3 , adjacent to microstructural cracks have a better photoelectrochemical response and reduced back electron recombination due to an optimal oxygen vacancy concentration, with the film thickness and extended light exposure also influencing local activity. Our work highlights the importance of microscopic mapping to understand activity, in even seemingly homogeneous photoelectrodes.

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

Soil gas probes for monitoring trace gas messengers of microbial activity

Abstract Soil microbes vigorously produce and consume gases that reflect active soil biogeochemical processes. Soil gas measurements are therefore a powerful tool to monitor microbial activity. Yet, the majority of soil gases lack non-disruptive subsurface measurement methods at spatiotemporal scales relevant to microbial processes and soil structure. To address this need, we developed a soil gas sampling system that uses novel diffusive soil probes and sample transfer approaches for high-resolution sampling from discrete subsurface regions. Probe sampling requires transferring soil gas samples to above-ground gas analyzers where concentrations and isotopologues are measured. Obtaining representative soil gas samples has historically required balancing disruption to soil gas composition with measurement frequency and analyzer volume demand. These considerations have limited attempts to quantify trace gas spatial concentration gradients and heterogeneity at scales relevant to the soil microbiome. Here, we describe our new flexible diffusive probe sampling system integrated with a modified, reduced volume trace gas analyzer and demonstrate its application for subsurface monitoring of biogeochemical cycling of nitrous oxide (N 2 O) and its site-specific isotopologues, methane, carbon dioxide, and nitric oxide in controlled soil columns. The sampling system observed reproducible responses of soil gas concentrations to manipulations of soil nutrients and redox state, providing a new window into the microbial response to these key environmental forcings. Using site-specific N 2 O isotopologues as indicators of microbial processes, we constrain the dynamics of in situ microbial activity. Unlocking trace gas messengers of microbial activity will complement -omics approaches, challenge subsurface models, and improve understanding of soil heterogeneity to disentangle interactive processes in the subsurface biome.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2025

Structural basis of differential gene expression at eQTLs loci from high-resolution ensemble models of 3D single-cell chromatin conformations

Abstract Motivation Techniques such as high-throughput chromosome conformation capture (Hi-C) have provided a wealth of information on nucleus organization and genome important for understanding gene expression regulation. Genome-Wide Association Studies have identified numerous loci associated with complex traits. Expression quantitative trait loci (eQTL) studies have further linked the genetic variants to alteration in expression levels of associated target genes across individuals. However, the functional roles of many eQTLs in noncoding regions remain unclear. Current joint analyses of Hi-C and eQTLs data lack advanced computational tools, limiting what can be learned from these data. Results We developed a computational method for simultaneous analysis of Hi-C and eQTL data, capable of identifying a small set of nonrandom interactions from all Hi-C interactions. Using these nonrandom interactions, we reconstructed large ensembles (×105) of high-resolution single-cell 3D chromatin conformations with thorough sampling, accurately replicating Hi-C measurements. Our results revealed many-body interactions in chromatin conformation at the single-cell level within eQTL loci, providing a detailed view of how 3D chromatin structures form the physical foundation for gene regulation, including how genetic variants of eQTLs affect the expression of associated eGenes. Furthermore, our method can deconvolve chromatin heterogeneity and investigate the spatial associations of eQTLs and eGenes at subpopulation level, revealing their regulatory impacts on gene expression. Together, ensemble modeling of thoroughly sampled single-cell chromatin conformations combined with eQTL data, helps decipher how 3D chromatin structures provide the physical basis for gene regulation, expression control, and aid in understanding the overall structure-function relationships of genome organization. Availability and implementation It is available at https://github.com/uic-liang-lab/3DChromFolding-eQTL-Loci.

Du, Lin (ORCID:0009000289869812)↗
DOE OSTI2021

Statistical upscaling of ecosystem CO 2 fluxes across the terrestrial tundra and boreal domain: Regional patterns and uncertainties

Abstract The regional variability in tundra and boreal carbon dioxide (CO 2 ) fluxes can be high, complicating efforts to quantify sink‐source patterns across the entire region. Statistical models are increasingly used to predict (i.e., upscale) CO 2 fluxes across large spatial domains, but the reliability of different modeling techniques, each with different specifications and assumptions, has not been assessed in detail. Here, we compile eddy covariance and chamber measurements of annual and growing season CO 2 fluxes of gross primary productivity (GPP), ecosystem respiration (ER), and net ecosystem exchange (NEE) during 1990–2015 from 148 terrestrial high‐latitude (i.e., tundra and boreal) sites to analyze the spatial patterns and drivers of CO 2 fluxes and test the accuracy and uncertainty of different statistical models. CO 2 fluxes were upscaled at relatively high spatial resolution (1 km 2 ) across the high‐latitude region using five commonly used statistical models and their ensemble, that is, the median of all five models, using climatic, vegetation, and soil predictors. We found the performance of machine learning and ensemble predictions to outperform traditional regression methods. We also found the predictive performance of NEE‐focused models to be low, relative to models predicting GPP and ER. Our data compilation and ensemble predictions showed that CO 2 sink strength was larger in the boreal biome (observed and predicted average annual NEE −46 and −29 g C m −2 yr −1 , respectively) compared to tundra (average annual NEE +10 and −2 g C m −2 yr −1 ). This pattern was associated with large spatial variability, reflecting local heterogeneity in soil organic carbon stocks, climate, and vegetation productivity. The terrestrial ecosystem CO 2 budget, estimated using the annual NEE ensemble prediction, suggests the high‐latitude region was on average an annual CO 2 sink during 1990–2015, although uncertainty remains high.

Virkkala, Anna‐Maria↗
DOE OSTI2023

Changing temperature profiles and the risk of dengue outbreaks

As temperatures change worldwide, the pattern and competency of disease vectors will change, altering the global distribution of both the burden of infectious disease and the risk of the emergence of those diseases into new regions. To evaluate the risk of potential summer dengue outbreaks triggered by infected travelers under various climate scenarios, we develop an SEIR-type model, run numerical simulations, and conduct sensitivity analyses under a range of temperature profiles. Our model extends existing theoretical frameworks for studying dengue dynamics by introducing temperature dependence of two key parameters: the mosquito extrinsic incubation period and the lifespan of mosquitoes, which empirical data suggests are both highly temperature dependent. We find that changing temperature significantly alters dengue risk in an inverted U-shape, with temperatures in the range 27-31°C producing the highest risk. As temperatures increase beyond 31°C, the determinants of dengue risk begin to shift from mosquito biting rate and carrying capacity to the duration of the human infectious period, suggesting that changing temperatures not only alter dengue risk but also the potential efficacy of control measures. To illustrate the role of spatial and temporal temperature heterogeneity, we select five US cities where the primary dengue vector, the mosquito Aedes aegypti , has been observed, and which have had dengue cases in the past: Los Angeles, Houston, Miami, Brownsville, and Phoenix. Our analysis suggests that an increase of 3°C leads to an approximate doubling of the risk of dengue in Los Angeles and Houston, but a reduction of risk in Miami, Brownsville, and Phoenix due to extreme heat.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Engineering of Nanoscale Heterogeneous Transition Metal Dichalcogenide–Au Interfaces

Engineering the transition metal dichalcogenide (TMD)–metal interface is critical for the development of two-dimensional semiconductor devices. By directly probing the electronic structures of WS 2 –Au and WSe 2 –Au interfaces with high spatial resolution, we delineate nanoscale heterogeneities in the composite systems that give rise to local Schottky barrier height modulations. Photoelectron spectroscopy reveals large variations (>100 meV) in TMD work function and binding energies for the occupied electronic states. Characterization of the composite systems with electron backscatter diffraction and scanning tunneling microscopy leads us to attribute these heterogeneities to differing crystallite orientations in the Au contact, suggesting an inherent role of the metal microstructure in contact formation. Here we then leverage our understanding to develop straightforward Au processing techniques to form TMD–Au interfaces with reduced heterogeneity. Our findings illustrate the sensitivity of TMDs’ electronic properties to metal contact microstructure and the viability of tuning the interface through contact engineering.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Comparison of automated chemical-guided segmentation and human annotation of soil organic matter in X-ray microcomputed tomography imaging in contrasted soil types

Soil organic matter (OM) formation and persistence is strongly influenced by the spatial distribution of organic substrates and microscale soil heterogeneity by dictating OM accessibility to microorganisms. However, traditional size and/or density fractionation techniques disrupt aggregate architecture, eliminating spatial information needed to fully understand intra-aggregate OM distribution. To quantify three-dimensional OM spatial distribution and automate segmentation in X-ray microcomputed tomography (µCT) imaging without human annotation bias, we developed an iodine gas vapor (I2) based staining workflow that eliminates labor-intensive manual annotation while maintaining segmentation accuracy, using aggregates from four taxonomically diverse soils (Xerofluvent, Haploxeroll Sphagnofibrist, Palehumult) with an 8-fold range of soil organic carbon. Human annotation of 10 µCT slices by the experienced and inexperienced annotators resulted in variations up to 3% in the Dice similarity coefficient (DSC), reflecting a degree of inherent subjectivity of manual labeling. Such inconsistencies are expected to compound as the number of manually annotated slices increases. Dual-energy µCT imaging at 33.1 keV (below the iodine (I) K-edge) and 33.2 keV (above the I K-edge) was used to resolve aggregate microstructure following I2 staining. The automated image subtraction pipeline identified OM regions by the I Kedge induced brightness increases, achieving DSC values of 0.58–0.83 relative to an experienced annotator. Sensitivity analyses revealed that the reconstruction alpha value—optimized via the open-source tool TomocuPy—and the 3D registration slice count were the primary determinants of accuracy, providing a novel benchmark for dual-energy soil imaging. The pipeline without GPU acceleration achieved 9.6 to 43.2 times faster than manual annotation. Using GPU-accelerated image post-processing and affine transformation matrices, the pipeline successfully segmented OM elements for large-scale datasets (3232×3232 pixel, 2048 slices) within ~5200 s from raw file acquisition to segmented output. The high-throughput approach enables the quantification of OM spatial distribution across diverse and heterogeneous soil.

Soil microbial biomass↗
DOE OSTI2020

Mapping of Heterogeneous Catalyst Degradation in Polymer Electrolyte Fuel Cells

Pt catalysts in polymer electrolyte fuel cells degrade heterogeneously as the catalyst particles are exposed to local variations throughout the catalyst layer during operation. State-of-the-art analytical techniques for studying degradation of Pt catalysts do not possess fine spatial resolution to elucidate such non-uniform degradation behavior at a large electrode level. A new methodology is developed to spatially resolve and quantify the heterogeneous Pt catalyst degradation over a large area (several cm ) of aged MEAs based on synchrotron X-ray microdiffraction. PEFC single cells are aged using voltage cycling as an accelerated stress test and the degradation heterogeneity at a micrometer length scale is visualized by mapping Pt catalyst particle size after voltage cycling. It is demonstrated in detail that the Pt catalyst particle size growth is non-uniform and follows the flow field geometry. The Pt particle size growth is greater in the area under the flow field land, while it is minimal in the area under the flow field channel. Additional non-uniformity is observed with the Pt particle size increasing more rapidly at the air outlet area than the Pt particle size at the inlet area.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Deep Learning Modeling Framework to Capture Mixing Patterns in Reactive-Transport Systems

Prediction and control of chemical mixing are vital for many scientific areas such as subsurface reactive transport, climate modeling, combustion, epidemiology, and pharmacology. Due to the complex nature of mixing in heterogeneous and anisotropic media, the mathematical models related to this phenomenon are not analytically tractable. Numerical simulations often provide a viable route to predict chemical mixing accurately. However, contemporary modeling approaches for mixing cannot utilize available spatial-temporal data to improve the accuracy of the future prediction and can be compute-intensive, especially when the spatial domain is large and for long-term temporal predictions. To address this knowledge gap, in this work we will present in this paper a deep learning (DL) modeling framework applied to predict the progress of chemical mixing under fast bimolecular reactions. This framework uses convolutional neural networks (CNN) for capturing spatial patterns and long short-term memory (LSTM) networks for forecasting temporal variations in mixing. By careful design of the framework—placement of non-negative constraint on the weights of the CNN and the selection of activation function, the framework ensures non-negativity of the chemical species at all spatial points and for all times. Our DL-based framework is fast, accurate, and requires minimal data for training. The time needed to obtain a forecast using the model is a fraction (≈ O(-6)) of the time needed to obtain the result using a high-fidelity simulation. To achieve an error of 10% (measured using the infinity norm) for capturing local-scale mixing features such as interfacial mixing, only 24% to 32% of the sequence data for model training is required. To achieve the same level of accuracy for capturing global-scale mixing features, the sequence data required for model training is 64% to 70% of the total spatial-temporal data. Hence, the proposed approach—a fast and accurate way to forecast long-time spatial-temporal mixing patterns in heterogeneous and anisotropic media—will be a valuable tool for modeling reactive-transport in a wide range of applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A Spatially Resolved Evaluation of Accelerated Environmental Aging on Emerging Polypropylene-Based Photovoltaic Backsheets Using Raman Spectroscopy

For this work, accelerated aging was used to assess environmental degradation in emerging co-extruded polypropylene (PP)-based backsheets under three different environmental conditions (65°C/20% relative humidity (RH), 75°C/20% RH, and 75°C/50% RH). Although differential scanning calorimetry did not measure crystallinity changes with exposure, spatially resolved Raman spectroscopy identified crystallinity increases in the core layer of aged samples, indicating a heterogeneous postcrystallization process. The Raman results were in agreement with synchrotron-based microfocused wide-angle X-ray scattering measurements. Cross-sectional nanoindentation was used to correlate localized crystallinity shifts with changes in Young's modulus. A similar trend was found where increased modulus was measured in the core layer, supporting the relationship between modulus and crystallinity. Finally, dielectric characterization was used to assess the impact of these material property changes on performance. While changes in the backsheet material properties and dielectric performance were observed with accelerated aging, these shifts generally equilibrated with time, indicating overall stability in response to environmental stressors. Additionally, the identified heterogeneous material property changes indicate that spatially resolved crystallinity measurements may be a valuable early failure indicator to be used in the assessment of PV backsheet long-term durability.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nonempirical Range-Separated Hybrid Functional with Spatially Dependent Screened Exchange

Electronic structure calculations based on density functional theory (DFT) have successfully predicted numerous ground-state properties of a variety of molecules and materials. However, exchange and correlation functionals currently used in the literature, including semilocal and hybrid functionals, are often inaccurate to describe the electronic properties of heterogeneous solids, especially systems composed of building blocks with large dielectric mismatch. Here, in this study, we present a dielectric-dependent range-separated hybrid functional, screened-exchange range-separated hybrid (SE-RSH), for the investigation of heterogeneous materials. We define a spatially dependent fraction of exact exchange inspired by the static Coulomb-hole and screened-exchange (COHSEX) approximation used in many-body perturbation theory, and we show that the proposed functional accurately predicts the electronic structure of several nonmetallic interfaces, three- and two-dimensional, pristine, and defective solids and nanoparticles.

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

Integrated parameter and process learning for hydrologic and biogeochemical modules in Earth System Models

Focus area: Primary focal area #2; secondary focal area #3: Learning about parameters and processes of land surface hydrologic and biogeochemical models in Earth System models by integrating machine learning, physics, and big data. Science challenges: How do we maximally leverage big-data observations to improve hydrobiogeochemical process description and parameterization so that such modules more realistically capture hydrologic and vegetation responses and feedbacks under the future climate? For example, how can we leverage physics, limited observations of vegetation and streamflow to better estimate evapotranspiration, and, relatedly, net primary productivity, especially for drought areas? Vegetation plays a critical role in regional and global water cycles; however, existing vegetation models have failed to predict vegetation response to droughts (McDowell & Xu, 2017) , arctic greening (Keenan & Riley, 2018) , and critical transitions between forest and savanna (Hirota et al., 2011) . These studies suggest that when we build process-based models (PBM) parameterized from regional and global plant traits, we tend to poorly describe plant adaptation and local-scale competition processes. The models and their associated parameters assigned for different regions in the world are not capturing essential heterogeneity in vegetation responses at finer spatial scales. Many parameters of the land surface models control hydrology and vegetation dynamics at the same time. The heterogeneity in vegetation response is a function of (i) plant type, (ii) plant size, (iii) competition and succession, (iv) environmental controls, and (v) local variations due to the unique ecological community that are very difficult to describe (e.g., the size of gaps resulting from fire that facilitated the coexistence of pioneering species). In the demographic models, only factors (i) and (iv) were captured, and plant types were generally described only by leaf phenology and climate zones. With current demographic models, we generally consider more traits to define plant types (i) and calibrate these traits to consider factors (ii), (iii) and (iv); however, it is substantially challenging to scale to regional and global simulations due to trait variations across space (Ali et al., 2016). Moreover, it has been noted that hillslope processes, including ridge-to-valley flow and sunny vs. shady slopes are primary organizers of water, energy, and vegetation (Clark et al., 2015; Fan et al., 2019) . Although gradual improvements in the hydrologic model component in earth system models may reduce this error (at a remarkably slow pace), the long-term, gradual impact of hydrology on plant traits are not well captured. Recent work showed that the hydrologic controls exerted by groundwater and lateral flow are primary regulators of rooting depth (Fan et al., 2017) . Such hydrologic controls have seldom been reflected in vegetation model parameterizations.

54 ENVIRONMENTAL SCIENCES↗