Search NASA⌕ Search

SEARCH · Search NASA

Results for “evapotranspiration”

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 127 records · Page 7

X-BASE: the first terrestrial carbon and water flux products from an extended data-driven scaling framework, FLUXCOM-X

Mapping in situ eddy covariance measurements of terrestrial land–atmosphere fluxes to the globe is a key method for diagnosing the Earth system from a data-driven perspective. We describe the first global products (called X-BASE) from a newly implemented upscaling framework, FLUXCOM-X, representing an advancement from the previous generation of FLUXCOM products in terms of flexibility and technical capabilities. The X-BASE products are comprised of estimates of CO 2 net ecosystem exchange (NEE), gross primary productivity (GPP), evapotranspiration (ET), and for the first time a novel, fully data-driven global transpiration product (ETT), at high spatial (0.05°) and temporal (hourly) resolution. X-BASE estimates the global NEE at −5.75 ± 0.33 Pg C yr −1 for the period 2001–2020, showing a much higher consistency with independent atmospheric carbon cycle constraints compared to the previous versions of FLUXCOM. The improvement of global NEE was likely only possible thanks to the international effort to increase the precision and consistency of eddy covariance collection and processing pipelines, as well as to the extension of the measurements to more site years resulting in a wider coverage of bioclimatic conditions. However, X-BASE global net ecosystem exchange shows a very low interannual variability, which is common to state-of-the-art data-driven flux products and remains a scientific challenge. With 125 ± 2.1 Pg C yr −1 for the same period, X-BASE GPP is slightly higher than previous FLUXCOM estimates, mostly in temperate and boreal areas. X-BASE evapotranspiration amounts to 74.7×10 3 ± 0.9×10 3 km 3 globally for the years 2001–2020 but exceeds precipitation in many dry areas, likely indicating overestimation in these regions. On average 57 % of evapotranspiration is estimated to be transpiration, in good agreement with isotope-based approaches, but higher than estimates from many land surface models. Despite considerable improvements to the previous upscaling products, many further opportunities for development exist. Pathways of exploration include methodological choices in the selection and processing of eddy covariance and satellite observations, their ingestion into the framework, and the configuration of machine learning methods. For this, the new FLUXCOM-X framework was specifically designed to have the necessary flexibility to experiment, diagnose, and converge to more accurate global flux estimates.

Nelson, Jacob A.↗

Data and Scripts associated with a manuscript on ecosystem responses to wildfires in the Columbia River Basin

This data package is associated with the publication “Ecosystem leaf area, gross primary production, and evapotranspiration responses to wildfire in the Columbia River Basin” submitted to Biogeosciences (Shi et al., 2024; doi: 10.22541/au.171053013.30286044/v1). In this research, data products, leaf area index (LAI), gross primary production (GPP), and evapotranspiration (ET), from the Moderate Resolution Imaging Spectroradiometer (MODIS) are used to quantify the resistance and resilience of different ecosystem types in the Columbia River Basin (CRB). A machine learning algorithm, random forest (RF), was used to examine the impacts of precipitation, vapor pressure deficit (VPD), and burn severity from Monitoring Trends in Burn Severity (MTBS) on ecosystem resilience. The data package includes the processed MODIS data products, precipitation, VPD, and burn severity in 138 fire regions in CRB and the input files for RF model training. This data package includes six folders. The MODIS products are included in three MODIS_* folders with shell scripts for data clipping and *ncl files for data processing: (1) “/MODIS_LAI_CRB”; (2) “/MODIS_GPP_CRB”; and (3) “/MODIS_ET_CRB”. All the processed data for each fire event are NetCDF formatted. The MTBS burn severity data and the shell and *ncl scripts used for data processing are in the folder named (4) “MTBS_fire”. The ERA meteorological fields and the data processing scritps are in (5) “ERA_Var_CR”. All the scripts for figure development are in the format of *ncl and in the folder (6) “paper_scripts”. See the file ending in “flmd.csv” for a list of all files contained in this data package and descriptions for each. Tabular column headers and units are described in the data dictionary file ending in “dd.csv”.

54 ENVIRONMENTAL SCIENCES↗

Meteorological Variables and Energy Fluxes at the Pumphouse Site, Crested Butte, CO 2017-2019

This data contains output from the pumphouse eddy covariance tower that includes shortwave radiation, longwave radiation, net radiation, air temperature, relative humidity, as well as sensible, latent, and ground heat fluxes. Also included is calculated evapotranspiration from the latent heat flux and the latent heat of vaporization. All data are on a daily timestep and displayed in Mountain Time. The data has been processed, and Quality Assurance / Quality Control (QA/QC) was done, but any daily gaps in the data have not been filled in. This research was funded by the Department of Energy and performed as part of the Watershed Function Scientific Focus Area. This research aimed to constrain evapotranspiration in a high-elevation catchment.The dataset includes one comma-separated values (CSV) data file (EddyCovariance_MeteorlogicalVariables_CrestedButtePumphouse.csv). Additionally, three metadata CSV files are included: (1) location metadata file (locations.csv), which contains location metadata and coordinates; (2) a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and (3) a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

Contribution of Tibetan Plateau ecosystems to local and remote precipitation through moisture recycling

The ecosystems of the Tibetan Plateau (TP) provide multiple important ecosystem services that benefit both local populations and those beyond, such as through climate regulation services on precipitation for East Asia and China. However, the precipitation regulation service of the TP ecosystems for supplying moisture and maintaining precipitation is yet to be evaluated. In this study, we used the moisture recycling framework and a moisture tracking model to quantify the precipitation regulation services of TP ecosystems for their contribution to precipitation. We found TP ecosystems contributed substantially to local and downwind precipitation, with a declined contribution from 221 mm/year for the TP and neighboring areas, to <10 mm/year for eastern China and other surrounding countries. Among ecosystem types, grassland contributed most to precipitation, followed by barren and snow lands, forests, and shrublands. In terms of seasonality, precipitation contribution from TP ecosystems was greater in summer months than in non-summer months for western China, while the opposite was true for eastern China – though the magnitude was much smaller. Over the past two decades, the significant evapotranspiration increases in TP translated to a widespread increase in precipitation contribution for TP and downwind beneficiary regions from 2000 to 2020. Our study provides a quantitative way to understand the precipitation regulation services of TP ecosystems through moisture recycling, substantiating their key role to safeguard precipitation and the water cycle for downwind regions– effectively acting as an ecological security shelter that could be perceived by the public.

Li, Yan↗

Data used in manuscript "Quantifying Subsurface Flow and Solute Transport in a Snowmelt-recharged Hillslope with Multiyear Water Balance", submitted to Water Resources Research on 5-26-2022

Data in this package were used to calculate subsurface flow and solute (specifically dissolved nitrogen) transport in the Pumphouse Lower Montane (PLM) hillslope. These data are included in separate CSV files associated with each figure presented in the manuscript. The files describe site locations and elevations, weather data (air temperatures and precipitation), modeled evapotranspiration, time dependent groundwater elevations, calculated annual subsurface flow, hydraulic conductivity values, time trends in subsurface flow within different zones, pore water concentrations of total dissolve nitrogen (TDN), rates of TDN transport, soil water matric potentials following snowmelt, soil water contents during drainage, and evapotranspiration calculated by soil water mass balance versus by CLM4 model.

54 ENVIRONMENTAL SCIENCES↗

Energy‐Water Asynchrony Principally Determines Water Available for Runoff From Snowmelt in Continental Montane Forests

ABSTRACT Changes in the volume, rate, and timing of the snowmelt water pulse have profound implications for seasonal soil moisture, evapotranspiration (ET), groundwater recharge, and downstream water availability, especially in the context of climate change. Here, we present an empirical analysis of water available for runoff using five eddy covariance towers located in continental montane forests across a regional gradient of snow depth, precipitation seasonality, and aridity. We specifically investigated how energy‐water asynchrony (i.e., snowmelt timing relative to atmospheric demand), surface water input intensity (rain and snowmelt), and observed winter ET (winter AET) impact multiple water balance metrics that determine water available for runoff (WAfR). Overall, we found that WAfR had the strongest relationship with energy‐water asynchrony (adjusted r 2 = 0.52) and that winter AET was correlated to total water year evapotranspiration but not to other water balance metrics. Stepwise regression analysis demonstrated that none of the tested mechanisms were strongly related to the Budyko‐type runoff anomaly (highest adjusted r 2 = 0.21). We, therefore, conclude that WAfR from continental montane forests is most sensitive to the degree of energy‐water asynchrony that occurs. The results of this empirical study identify the physical mechanisms driving variability of WAfR in continental montane forests and are thus broadly relevant to the hydrologic management and modelling communities.

Webb, Ryan William↗

Systematic Evaluation of Atmospheric Forcing, Surface Datasets, and Mesh Effects on Kilometer-Scale Land Surface and River Modeling

Earth system models are advancing toward kilometer-scale resolution to capture local climate impacts and extremes. High-resolution land and river modeling depends on multiple factors, including mesh, surface datasets, and atmospheric forcing, but their relative effects at kilometer scales remain unquantified. We evaluated five Energy Exascale Earth System Model land and river configurations over the Mid-Atlantic region using two mesh (1/8° structured versus variable-resolution unstructured mesh), two surface datasets (default versus newly developed), and three atmospheric forcings (NLDAS2, MSWX, GSWP). Evaluation against satellite, reanalysis, and in situ benchmarks across water, energy, and carbon cycles quantifies how these factors affect model performance. Forcing selection produces the largest bias reductions (12-99% across variables), followed by surface datasets (7-75%) and mesh (up to 21%). Forcing effects vary by variable, with MSWX reducing biases for snow water equivalent, evapotranspiration, albedo, temperature, and gross primary productivity, GSWP for snow cover and runoff, and NLDAS for soil moisture and streamflow. The use of newly developed surface datasets improves gross primary productivity (58% bias reduction) and evapotranspiration but increase soil moisture and albedo biases due to current modeling limitations. Variable-resolution unstructured mesh improves the simulation of small-basin streamflow through better capturing drainage networks, though mesh minimally affects other land variables. These findings provide important guidance for high-resolution modeling development and actionable science.

Land and River modeling↗

Exploring Micro-Environmental Conditions of Urban Agrivoltaics: Advancing Sustainable Green Spaces and Agriculture in Temperate Cities

Agrivoltaics, the integrated land use combining renewable energy production and agriculture, can potentially address key challenges faced by urban communities, including limited access to vacant land, fresh produce, and electricity. However, urban agrivoltaics has rarely been considered a viable solution, primarily due to a lack of experimental systems demonstrating the functionality of urban food production or green spaces alongside solar energy production. To evaluate the potential of solar photovoltaics to help mitigate heat stress on urban agriculture and green spaces in a temperate climate, we examined the microclimatic changes introduced by an urban agrivoltaics system (UrAV) when combined with vegetable crops and turfgrass, respectively. Accordingly, we installed an experimental setup instrumented with environmental sensors to compare full-sun conditions with those of an agrivoltaics system designed for urban environments. We found minimal differences for both land cover types between the control and agrivoltaics sites in air and soil temperature. However, during periods of intensified summer heat, temperatures beneath the panels cooled due to reductions in wind speed and relative humidity, which impeded the movement of hot, moist air and reduced reference evapotranspiration. Further, soil moisture in UrAV was highly spatially heterogeneous, influenced by the layout of the photovoltaic panels and their support structure, which redistributed rainfall and controlled where solar radiation could penetrate and drive evapotranspiration. Overall, our results suggest that PV-induced environmental changes in temperate climates are compatible with cultivating turfgrass or appropriate crops. These findings can help planners and designers integrate agrivoltaics into community gardens, farms, and green spaces in temperate cities.

14 SOLAR ENERGY↗

A Retrospective and Prospective Examination of the 1960s U.S. Northeast Drought

As the most severe drought over the Northeastern United States (NEUS) in the past century, the 1960s drought had pronounced socioeconomic impacts. Although it was followed by a persisting wet period, the conditions leading to the 1960s extreme drought could return in the future, along with its challenges to water management. To project the potential consequences of such a future drought, pseudo-global warming simulations using the Weather Research and Forecasting Model are performed to simulate the dynamical conditions of the historical 1960s drought, but with modified thermodynamic conditions under the shared socioeconomic pathway SSP585 scenario in the early (2021–2027), middle (2041–2047), and late (2091–2097) 21st century. Our analysis focuses on essential hydroclimatic variables including temperature, precipitation, evapotranspiration, soil moisture, snowpack, and surface runoff. In contrast to the historical 1960s drought, similar dynamical conditions will generally produce more precipitation, increased soil moisture and evapotranspiration, and reduced snowpack. However, we also find that although wet months get much wetter, dry months may become drier, meaning that wetting trends are most significant in wet months but are essentially negligible for extremely dry months with negative monthly mean net precipitation. For these months, the trend toward wetting conditions provides little relief from the effects of extreme dry months. These conditions may even aggravate water shortages due to an increasingly rapid transition from wet to dry conditions. Other challenges emerge for residents and stakeholders in this region, including more extreme hot days, record-low snow pack, frozen ground degradation and subsequent decreases in surface runoff.

54 ENVIRONMENTAL SCIENCES↗

Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in situ and remote sensing data sets

The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).

Cheng, Yanyan↗

High Vapor Pressure Deficit Decreases the Productivity and Water Use Efficiency of Rain‐Induced Pulses in Semiarid Ecosystems

Abstract Intermittent rain events drive dynamic pulses of carbon and water exchange in many arid and semiarid ecosystems. Although soil moisture is known to control these pulses, the effect of atmospheric dryness on pulses is not well documented. Here we hypothesized that vapor pressure deficit ( VPD ) modulates net ecosystem production ( NEP ) and ecosystem‐scale water use efficiency ( WUE ) during pulse events due to its effects on canopy stomatal conductance and evapotranspiration. We quantified relationships between VPD and carbon and water exchange during growing season rain events and tested their generality across four semiarid flux sites with varied vegetation in the southwest United States. Across grassland, shrubland, and savanna sites, we found that high VPD during pulses suppressed ecosystem photosynthesis and surface conductance to a greater degree than respiration or evapotranspiration, particularly when soil moisture was high. Thus, periods of high VPD were associated with a 13–64% reduction in NEP and an 11–25% decrease in WUE , relative to moderate VPD conditions. Sites dominated by shrubs with the C3 photosynthetic pathway were more sensitive to VPD than sites dominated by C4 grasses. We found that a 1 kPa increase in VPD reduced the average NEP of pulse events by 13–56%, which illustrates the potential for projected increases in atmospheric demand to reduce the net productivity of semiarid ecosystems.

Roby, Matthew C.↗

Evaluation of Methods for Causal Discovery in Hydrometeorological Systems

Understanding causal relations is of utmost importance in hydrology and climate research for systems identification, prediction, and understanding systems behavior in a changing climate. Traditionally, researchers in hydrometeorology attempted to study causal questions by conducting controlled experiments using numerical models. This approach, however, in most cases of interest provides uncertain results because the models are approximate representation of the natural system. An alternative approach that has recently drawn significant attention in several fields is to infer causal relations from purely observational data. It possesses several traits to its utility particularly in hydrometeorology due to the rapid accumulation of in situ and remotely sensed data records. The first objective of this study is to present a brief description of four causal discovery methods (Granger causality, Transfer Entropy, graph-based algorithms, and Convergent Cross Mapping) with special emphasis on the assumptions on which they are built. Second, using synthetic data generated from a hydrological model, we assess their performance in retrieving causal information taking into account sensitivity to sample size and presence of noise. Last, we use causal analysis to examine and formulate hypotheses on causal drivers of evapotranspiration in a shrubland region during summer and winter seasons. An interpretation of the hypotheses based on canopy seasonal dynamics and evapotranspiration processes is presented. It is hoped that the results presented here can be useful in guiding researchers studying hydrometeorological systems as to which causal method is most appropriate to the characteristics of the system under study.

54 ENVIRONMENTAL SCIENCES↗

Permafrost Promotes Shallow Groundwater Flow and Warmer Headwater Streams

Abstract The presence of permafrost influences the flow paths of water through Arctic landscapes and thereby has the potential to impact stream discharge and thermal regimes. Observations from 11 headwater streams in Alaska showed that July water temperatures were higher in catchments with more near‐surface permafrost. We apply a fully coupled cryohydrology model to investigate if the impact of permafrost on flow path depth could cause the same pattern in temperatures of groundwater discharging from hillslopes to streams. The model simulates surface energy and water balances, snow, and subsurface water and energy balances for two‐dimensional hillslope model cases with varying permafrost extent. We find that hillslopes with continuous permafrost have more shallow flow paths and twice as high rates of evapotranspiration, compared to hillslopes with no permafrost. For our simulated cases, 6.7% of the horizontal water flux moves through the top organic soil layers when there is continuous permafrost, while only 0.5% moves through organic layers without permafrost. The deeper flow paths in permafrost‐free simulations buffer seasonal temperature extremes, so that summer groundwater discharge temperatures are highest with continuous permafrost. Our results suggest that permafrost thawing alters groundwater flow paths and can lead to decreases in summer stream temperatures and reductions in evapotranspiration in headwater catchments. These changes are of potential importance for stream biotic components of ecosystems, however, the full impact remains unknown.

54 ENVIRONMENTAL SCIENCES↗

Predicting Characteristics of the Water Cycle From Scaling Relationships

Over multi-decadal time scales, assuming that changes in subsurface water storage are negligible, the continental precipitative water flux, P, can be divided into two principal components, Q (run-off, including soil infiltration and groundwater recharge) and ET (evapotranspiration). Taking into account a broadly applied Budyko's phenomenology to describe the relationship of ET/P as a function of PET/P, where PET is the potential evapotranspiration, we propose a theoretical framework for predicting characteristics of the water cycle from scaling relationships. In this framework, the ecosystem net primary productivity is expressed in terms of soil formation and vegetation growth, which is mathematically optimized with respect to the water partitioning, generating directly the value ET/P. Here, the mathematical optimization is based on the general ecological principle that dominant ecosystems tend to be those that, for any given conditions, maximize conversion of atmospheric carbon to biomass. It is shown that application of the results of mathematical optimization to water-limited ecosystems is possible by applying the optimization only to a vegetation covered portion of the surface. For energy-limited ecosystems, the optimization can be applied only to a portion of the precipitation equal to PET, assuming that the remaining P simply runs off. We use theoretical and actual values of plant root fractal dimensionalities, d f , to predict ranges of ET/P as a function of PET/P for 0 ≤ PET/P ≤ 1 and compare with annual and multi-decadal means of ET/P. By comparing the developed approach with a large amount of data collected from the literature, we demonstrate its successful applications to both water- and energy-limited systems.

54 ENVIRONMENTAL SCIENCES↗

Differentiable, Learnable, Regionalized Process-Based Models With Multiphysical Outputs can Approach State-Of-The-Art Hydrologic Prediction Accuracy

Predictions of hydrologic variables across the entire water cycle have significant value for water resources management as well as downstream applications such as ecosystem and water quality modeling. Recently, purely data-driven deep learning models like long short-term memory (LSTM) showed seemingly insurmountable performance in modeling rainfall runoff and other geoscientific variables, yet they cannot predict untrained physical variables and remain challenging to interpret. Here, we show that differentiable, learnable, process-based models (called δ models here) can approach the performance level of LSTM for the intensively observed variable (streamflow) with regionalized parameterization. We use a simple hydrologic model HBV as the backbone and use embedded neural networks, which can only be trained in a differentiable programming framework, to parameterize, enhance, or replace the process-based model's modules. Without using an ensemble or post-processor, δ models can obtain a median Nash-Sutcliffe efficiency of 0.732 for 671 basins across the USA for the Daymet forcing data set, compared to 0.748 from a state-of-the-art LSTM model with the same setup. For another forcing data set, the difference is even smaller: 0.715 versus 0.722. Meanwhile, the resulting learnable process-based models can output a full set of untrained variables, for example, soil and groundwater storage, snowpack, evapotranspiration, and baseflow, and can later be constrained by their observations. Both simulated evapotranspiration and fraction of discharge from baseflow agreed decently with alternative estimates. The general framework can work with models with various process complexity and opens up the path for learning physics from big data.

54 ENVIRONMENTAL SCIENCES↗

Impacts of Topography‐Driven Water Redistribution on Terrestrial Water Storage Change in California Through Ecosystem Responses

Abstract Lateral subsurface flow plays an essential role in sustaining the terrestrial ecosystem, but it is not explicitly represented in most Earth System Models. In this study, we implemented an explicit lateral saturated flow model into the E3SM land model (ELM). The model explicitly describes lateral flow in the saturated zone by representing, for each model grid, an idealized hillslope consisting of five hydrologically connected soil columns. We conducted three model experiments driven by 0.125° atmospheric forcing data during 1980–2015 over California using models of the default ELM, a modified version of ELM to enhance infiltration, and the model with the lateral saturated flow model. The simulated runoff, evapotranspiration, and terrestrial water storage anomaly (TWSA) from the three simulations were evaluated against available observations, and the model explicitly representing lateral flow performs best. The new model produces greater gridcell‐averaged evapotranspiration especially over the mountainous regions with moderate relief and seasonally dry climates. Most importantly, it improves the modeled seasonal variations, interannual variabilities, and the recent decadal decline of TWSA. Many of these improvements can be attributed to the enhanced ecosystem resilience to droughts as demonstrated by transpiration increases caused by lateral flow. Model sensitivity experiments suggest that subsurface runoff is most sensitive to the ratio between horizontal and vertical saturated hydraulic conductivity, followed by hillslope planforms (convergent, divergent, and uniform), number of columns, and lower boundary conditions. Future work should effectively characterize hillslopes in global models and explore the long‐term influences of lateral water movement on modeled biogeochemical cycle.

54 ENVIRONMENTAL SCIENCES↗

The Role of Vapor Pressure Deficit in the CLM Simulated Interaction Between Urban Heat Islands and Heat Waves Over CONUS

Abstract Although both urban and rural temperatures are expected to increase under heat waves (HWs), whether the urban heat island (UHI) intensity becomes stronger under HWs remains unknown especially at the daily mean and large spatial scales. Using an urbanized land surface model, we quantify the interactions between UHIs and HWs over the Contiguous United States (CONUS). Synergistic interactions (i.e., increased UHI intensities under HWs) are observed over the eastern and western U.S. However, negative interactions are found in the Central U.S. due to the stronger inhibition of rural evapotranspiration by vapor pressure deficit (VPD) stresses. The interactions between UHIs and HWs in the Central U.S. will be further reduced along with the elevated VPD stresses in a hotter future. The results highlight the importance of properly parameterizing the sensitivity of urban and rural evapotranspiration to various environmental stresses in climate and earth system models.

Wang, Linying [National Key Laboratory of Earth Sy↗

Quantifying Groundwater Response and Uncertainty in Beaver‐Influenced Mountainous Floodplains Using Machine Learning‐Based Model Calibration

Abstract Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model‐data integration workflow to address uncertainty in groundwater response to beaver‐induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning‐based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down‐valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down‐valley underflow, from 2 during dry pond periods to 20 during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down‐valley underflow, minimizes the effect of beaver‐induced inundation on water quality. We emphasize the need for field‐scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response. Plain Language Summary Beavers change the flow of water in river corridors by creating ponds, expanding wetlands, and flooding floodplains. This increases surface water area, promotes plant growth, and enhances biodiversity. However, the impact of this flooding on groundwater flow is not well understood, especially in mountainous areas with gravel layers where water moves easily beneath soil. In this study, we used numerical modeling to investigate how beaver ponds influence groundwater in a mountainous floodplain of the Upper Colorado River Basin. We adapted a machine learning method to validate our numerical models using multiple field data sets. Our findings show that beaver ponds significantly increase vertical water flow from the soil to the gravel during wet periods, compared to when the ponds are fully drained. The study also highlights the importance of floodplain structure in controlling both water flow in gravel layers along the river direction and vertical flow from the soil to the gravel with the presence of beavers. To reduce uncertainty in groundwater response, we emphasize the need for more field‐scale measurements of floodplain structure, hydraulic properties, and evapotranspiration changes. Key Points Floodplain structures and hydraulic conductivities are important for groundwater response with beaver ponds in mountainous floodplains Large down‐valley underflow in permeability‐stratified floodplains reduces beaver‐induced impacts on groundwater storage and water quality Machine learning‐based model calibration methods are effective for estimating posterior distributions of groundwater model parameters

Wang, Lijing↗