The ability to accurately monitor and anticipate changes in consumptive water use associated with changing land use and land management is critical to developing sustainable water management strategies in water-limited climatic regions. In this paper, we present an application of a remote sensing data fusion technique for developing high spatiotemporal resolution maps of evapotranspiration (ET) at scales that can be associated with changes in land use. The fusion approach combines ET map timeseries developed using an multi-scale energy balance algorithm applied to thermal data from Earth observation platforms with high spatial but low temporal resolution (e.g., Landsat) and with moderate resolution but frequent temporal coverage (e.g., MODIS (Moderate Resolution Imaging Spectroradiometer)). The approach is applied over the Sacramento-San Joaquin Delta region in California—an area critical to both agricultural production and drinking water supply within the state that has recently experienced stresses on water resources due to a multi-year (2012–2017) extreme drought. ET “datacubes” with 30-m resolution and daily timesteps were constructed for the 2015–2016 water years and related to detailed maps of land use developed at the same spatial scale. The ET retrievals are evaluated at flux sites over multiple land covers to establish a metric of accuracy in the annual water use estimates, yielding root-mean-square errors of 1.0, 0.8, and 0.3 mm day⁻¹ at daily, monthly, and yearly timesteps, respectively, for all sites combined. Annual ET averaged over the Delta changed only 3 mm year⁻¹ between water years, from 822 to 819 mm year⁻¹, translating to an area-integrated total change in consumptive water use of seven thousand acre-feet (TAF). Changes were largest in areas with recorded land-use change between water years—most significantly, fallowing of crop land presumably in response to reductions in water availability and allocations due to the drought. Moreover, the time evolution in water use associated with wetland restoration—an effort aimed at reducing subsidence and carbon emissions within the inner Delta—is assessed using a sample wetland chronosequence. Region-specific matrices of consumptive water use associated with land use changes may be an effective tool for policymakers and farmers to understand how land use conversion could impact consumptive use and demand.

As a primary flux in the global water cycle, evapotranspiration (ET) connects hydrologic and biological processes and is directly affected by water and land management, land use change and climate variability. Satellite remote sensing provides an effective means for diagnosing ET patterns over heterogeneous landscapes; however, limitations on the spatial and temporal resolution of satellite data, combined with the effects of cloud contamination, constrain the amount of detail that a single satellite can provide. In this study, we describe an application of a multi-sensor ET data fusion system over a mixed forested/agricultural landscape in North Carolina, USA, during the growing season of 2013. The fusion system ingests ET estimates from the Two-Source Energy Balance Model (TSEB) applied to thermal infrared remote sensing retrievals of land surface temperature from multiple satellite platforms: hourly geostationary satellite data at 4 km resolution, daily 1 km imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS) and biweekly Landsat thermal data sharpened to 30 m. These multiple ET data streams are combined using the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) to estimate daily ET at 30 m resolution to investigate seasonal water use behavior at the level of individual forest stands and land cover patches. A new method, also exploiting the STARFM algorithm, is used to fill gaps in the Landsat ET retrievals due to cloud cover and/or the scan-line corrector (SLC) failure on Landsat 7. The retrieved daily ET time series agree well with observations at two AmeriFlux eddy covariance flux tower sites in a managed pine plantation within the modeling domain: US-NC2 located in a mid-rotation (20-year-old) loblolly pine stand and US-NC3 located in a recently clear-cut and replanted field site. Root mean square errors (RMSEs) for NC2 and NC3 were 0.99 and 1.02 mm day⁻¹, respectively, with mean absolute errors of approximately 29 % at the daily time step, 12 % at the monthly time step and 0.7 % over the full study period at the two flux tower sites. Analyses of water use patterns over the plantation indicate increasing seasonal ET with stand age for young to mid-rotation stands up to 20 years, but little dependence on age for older stands. An accounting of consumptive water use by major land cover classes representative of the modeling domain is presented, as well as relative partitioning of ET between evaporation (E) and transpiration (T) components obtained with the TSEB. The study provides new insights about the effects of management and land use change on water yield over forested landscapes.

The phenological response of vegetation to ongoing climate change may have great implications for hydrological regimes in the eastern United States. However, there have been few studies that analyze its resultant effect on catchment discharge dynamics, separating from dominant climatic controls. In this study, we examined the net effect of phenological variations on the long‐term and interannual gross primary production (GPP) and evapotranspiration (ET) fluxes in a temperate deciduous forest, as well as on the catchment discharge behavior in a mixed deciduous‐conifer forest catchment. First, we calibrated the spring and autumn leaf phenology models for the Harvard Forest in the northeastern United States, where the onsets of greenup and senescence have been significantly advanced and delayed, 10.3 and 6.0 days respectively, over the past two decades (1992–2011). We then integrated the phenology models into a mechanistic watershed ecohydrological model (RHESSys), which improved the interannual and long‐term simulations of both the plot‐scale daily GPP and ET fluxes and the catchment discharge dynamics. We found that the phenological changes amplified the long‐term increases in GPP and ET driven by the climatic controls. In particular, the earlier greenup onsets resulted in increases in annual ET significantly, while the delayed senescence onsets had less influence. Consequently, the earlier greenup onsets reduced stream discharge not only during the growing season but also during the following dormant season due to soil water depletion. This study highlights the importance of understanding vegetation response to ongoing climate change in order to predict the future hydrological nonstationarity in this region.

The energy delivered to the land surface via insolation is a primary driver of evapotranspiration (ET)—the exchange of water vapor between the land and atmosphere. Spatially distributed ET products are in great demand in the water resource management community for real-time operations and sustainable water use planning. The accuracy and deliverability of these products are determined in part by the characteristics and quality of the insolation data sources used as input to the ET models. This paper investigates the practical utility of three different insolation datasets within the context of a satellite-based remote sensing framework for mapping ET at high spatiotemporal resolution, in an application over the Sacramento–San Joaquin Delta region in California. The datasets tested included one reanalysis product: The Climate System Forecast Reanalysis (CFSR) at 0.25° spatial resolution, and two remote sensing insolation products generated with geostationary satellite imagery: a product for the continental United States at 0.2°, developed by the University of Wisconsin Space Sciences and Engineering Center (SSEC) and a coarser resolution (1°) global Clouds and the Earth’s Radiant Energy System (CERES) product. The three insolation data sources were compared to pyranometer data collected at flux towers within the Delta region to establish relative accuracy. The satellite products significantly outperformed CFSR, with root-mean square errors (RMSE) of 2.7, 1.5, and 1.4 MJ·m⁻²·d⁻¹ for CFSR, CERES, and SSEC, respectively, at daily timesteps. The satellite-based products provided more accurate estimates of cloud occurrence and radiation transmission, while the reanalysis tended to underestimate solar radiation under cloudy-sky conditions. However, this difference in insolation performance did not translate into comparable improvement in the ET retrieval accuracy, where the RMSE in daily ET was 0.98 and 0.94 mm d⁻¹ using the CFSR and SSEC insolation data sources, respectively, for all the flux sites combined. The lack of a notable impact on the aggregate ET performance may be due in part to the predominantly clear-sky conditions prevalent in central California, under which the reanalysis and satellite-based insolation data sources have comparable accuracy. While satellite-based insolation data could improve ET retrieval in more humid regions with greater cloud-cover frequency, over the California Delta and climatologically similar regions in the western U.S., the CFSR data may suffice for real-time ET modeling efforts.