The ability to quantify evapotranspiration (ET) is crucial for smart agriculture and sustainable groundwater management. Efficient ET estimation strategies often rely on the vertical‐flow assumption to assimilate data from soil‐moisture sensors. While adequate in some large‐scale applications, this assumption fails when the horizontal component of the local flow velocity is not negligible due to, for example, soil heterogeneity or drip irrigation. We present novel implementations of the ensemble Kalman filter (EnKF) and the maximum likelihood estimation (MLE), which enable us to infer spatially varying ET rates and root water uptake profiles from soil‐moisture measurements. While the standard versions of EnKF and MLE update the predicted soil moisture prior to computing ET, ours treat the ET sink term in Richards' equation as an updatable observable. We test the prediction accuracy and computational efficiency of our methods in a setting representative of drip irrigation. Our strategies accurately estimate the total ET rates and root‐uptake profiles and do so up to two‐orders of magnitude faster than the standard EnKF.
Evapotranspiration is arguably the least quantified component of the hydrologic cycle. We propose two complementary strategies for estimation of evapotranspiration rates and root water uptake profiles from soil‐moisture sensor‐array data. One is our implementation of ensemble Kalman filter (EnKF); it treats the evapotranspiration sink term in the Richards equation, rather than soil moisture, as the observable to update. The other is a maximum likelihood estimator (MLE) applied to the same observable; it is supplemented with the Fisher information matrix to quantify uncertainty in its predictions. We use numerical experiments to demonstrate the accuracy and computational efficiency of these techniques. We found our EnKF implementation to be two orders of magnitude faster than either the standard EnKF or MLE, and our MLE procedure to require an order of magnitude fewer iterations to converge than its counterpart applied to soil moisture. These findings render our methodologies a viable and practical tool for estimation of the root water uptake profiles and evaporation rates, with the MLE technique to be used when the prior knowledge about evapotranspiration at the site is elusive.
more » « less- PAR ID:
- 10446300
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 57
- Issue:
- 11
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
In the semiarid, water-limited deserts of the Southwest United States, soil moisture is a crucial factor influencing atmospheric, hydrologic, and ecological processes. These dynamics are driven by infrequent yet significant precipitation events that redistribute moisture and establish hydrologic connectivity across the landscape. The Chihuahuan Desert, particularly within its endorheic basins, exemplifies these large-scale interactions where a complex balance of hydrological fluxes is maintained within a closed system. These basins receive most of their precipitation in upland regions, from which surface runoff can lead to downstream connectivity. This connectivity is influenced by the local water balance, including interactions among precipitation, leakage, and evapotranspiration, which are essential for understanding soil moisture variability. Additionally, soil moisture is affected by soil profile characteristics, vegetation, and atmospheric conditions. Field-scale methods like Cosmic-Ray Neutron Sensing (CRNS) are more appropriate than point-scale in situ sensors for quantifying hydrologic connectivity between upland and downstream regions, as CRNS reliably captures soil moisture temporal dynamics over several hectares. This study examines these dynamics within the endorheic Jornada Basin of the Chihuahuan Desert, focusing on two contrasting sites: an Upland Watershed (UW) on a piedmont slope and a Downstream Playa (DP) in a valley bottom. Using CRNS and complementary water balance instrumentation, I compared soil moisture dynamics at these two sites from July 2022 to February 2024. My analysis centered on a significant precipitation event early in the study period that generated surface runoff and playa inundation, followed by an extended dry period. Although temporal variations in leakage and evapotranspiration are similar at both sites, their rates differ significantly. The UW experienced a higher drying rate, necessitating greater plant water uptake from the subsurface. This led to an increased upward leakage to sustain vegetation, resulting in a leakage value of -205 mm, indicating vertical plant water uptake. Conversely, at the DP, the inundation event was formed by 228 mm of surface runoff, supplementing water inputs from precipitation. This additional water reduced the need for upward soil water movement to sustain plant water uptake, resulting in a leakage value of -97 mm. These findings enhance our understanding of hydrologic fluxes within endorheic basins and improve the applicability of hydrological models and the downscaling of remotely sensed soil moisture products.more » « less
-
Crop growth depends on the root-zone soil moisture (RZSM) (~top 1m). Accurate estimation of RZSM is vital to optimize irrigation management for saving water and energy while sustaining crop yield. The High-Resolution Land Assimilation System (HRLDAS) from NCAR can generate RZSM at field scales for irrigation management. The soil moisture data from various agriculture sites in the AmeriFlux network, U.S. Climate Reference Network (USCRN), and Soil Climate Analysis Network (SCAN) are used to verify the soil moisture products generated by HRLDAS. Although the HRLDAS products is not location specific and could be applied nationwide, this study will focus on Nebraska for evaluation, validation, and further calibration. We also compared NASA’s SMAP surface soil moisture products to HRLDAS surface layer soil moisture. Since the accuracy of the SMAP product is known, this comparison directly validates the HRLDAS surface soil moisture product and indirectly validate its RZSM products. Results from these two validation methods show a good accuracy of HRLDAS soil moisture products. The conspicuous differences between HRLDAS and SMAP products indicate that HRLDAS omits the irrigation activities as its simulation is based on weather variables and energy balance. It’s hard for HRLDAS to consider and include the irrigation actions in its results, while as SMAP products remotely sense the soil moisture as it is, the changes caused by irrigation are clearly reflected. Therefore, a simple calibration is applied to the HRLDAS products by including irrigation amount as its variables.more » « less
-
Abstract The subsurface processes that mediate the connection between evapotranspiration and groundwater within forested hillslopes are poorly defined. Here, we investigate the origin of diel signals in unsaturated soil water, groundwater, and stream stage on three forested hillslopes in the H.J. Andrews Experimental Forest in western Oregon, USA, during the summer of 2017, and assess how the diurnal signal in evapotranspiration (ET) is transferred through the hillslope and into these stores. There was no evidence of diel fluctuations in upslope groundwater wells, suggesting that tree water uptake in upslope areas does not directly contribute to the diel signal observed in near‐stream groundwater and streamflow. The water table in upslope areas resided within largely consolidated bedrock, which was overlain by highly fractured unsaturated bedrock. These subsurface characteristics inhibited formation of diel signals in groundwater and impeded the transfer of diel signals in soil moisture to groundwater because (1) the bedrock where the water table resides limited root penetration and (2) the low unsaturated hydraulic conductivity of the highly fractured rock weakened the hydraulic connection between groundwater and soil/rock moisture. Transpiration‐driven diel fluctuations in groundwater were limited to near‐stream areas but were not ubiquitous in space and time. The depth to the groundwater table and the geologic structure at that depth likely dictated rooting depth and thus controlled where and when the transpiration‐driven diel fluctuations were apparent in riparian groundwater. This study outlines the role of hillslope hydrogeology and its influence on the translation of evapotranspiration and soil moisture fluctuations to groundwater and stream fluctuations.
-
Abstract Soil hydrology provides important background for understanding the fate of organic carbon (OC) buried by geomorphic processes as well as the influence of runoff, infiltration, and plant root uptake on long‐term erosion and landscape evolution. We modeled the hydrology of a 4.5‐m loess‐paleosol sequence on an eroding tableland in the U.S. central Great Plains using Hydrus 1D, a numerical unsaturated flow model, parameterized with high resolution measurements of the soil water retention and hydraulic conductivity curves, which were distinct for the loess and paleosols. We hypothesized that (a) the connection of paleosols to modern climate depends on their burial depth, (b) paleosols in the root zone would have broader pore‐size distributions than unweathered loess, and (c) this broader pore‐size distribution increased root water uptake and made vegetation more resilient to drought, increasing the stability of loess tablelands despite high erodibility and high local relief. Four years with varying total annual precipitation were simulated for the observed profile and two hypothetical profiles, one without paleosols and another with a shallow, strongly developed paleosol. In these simulations, soil moisture in shallow paleosols responds quickly to precipitation while a deeply buried paleosol is largely disconnected from the modern climate, contributing to buried OC preservation. Contrary to our expectation, the presence of paleosols did not increase root uptake relative to unweathered loess in either wet or dry years. The unweathered coarse loess we studied may have an optimal pore‐size distribution for root uptake, providing an alternative hypothesis for why highly erodible loess tablelands persist.