skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on September 1, 2026

Title: Spatial heterogeneity of agricultural evapotranspiration as quantified by satellite and flux tower sources
Satellite-based evapotranspiration (ET) products inform decision-making regarding water availability and plant water use from sub-field to watershed scales. These products are validated using eddy covariance flux towers, where observations are subject to the influence of landscape heterogeneity due to a constantly shifting contributing area, or flux footprint, governed by surface and atmospheric drivers. In agricultural regions with multiple crop types, local heterogeneity may amplify the importance of appropriate grid cell selection for satellite ET comparisons. We evaluate the extent to which different satellite ET products capture both ET magnitudes and spatial variability due to landscape heterogeneity. We compare ECOSTRESS 70m instantaneous and daily ET products to flux tower observations in central Illinois with measurements at two heights, representing different but overlapping flux footprint areas. We estimate satellite ET based on a flux footprint weighting, gridded averages around the tower, and for crop-specific corn or soybean grid cells. In this region, the disALEXI daily ET product has better overall performance relative to the PT-JPL daily product, which is more sensitive to overpass times and tends to overestimate instantaneous ET. However, the PT-JPL model reflects more variability at high spatial resolution. Specifically, a footprint-derived ET improves the data-model comparison relative to disALEXI, and PT-JPL more closely replicates crop-specific and field-specific differences as inferred from the 2-height experimental setup. This study highlights differences in how models integrate spatial inputs, which lead to different representations of spatial variability for the same nominal resolution. This can also have important implications for understanding and predicting field-level differences in land-atmosphere fluxes.  more » « less
Award ID(s):
2012850
PAR ID:
10630772
Author(s) / Creator(s):
; ; ;
Editor(s):
NA
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Agricultural and Forest Meteorology
Volume:
372
Issue:
C
ISSN:
0168-1923
Page Range / eLocation ID:
110608
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Spatial heterogeneity in methane (CH 4 ) flux requires a reliable upscaling approach to reach accurate regional CH 4 budgets in the Arctic tundra. In this study, we combined the CLM-Microbe model with three footprint algorithms to scale up CH 4 flux from a plot level to eddy covariance (EC) tower domains (200 m × 200 m) in the Alaska North Slope, for three sites in Utqiaġvik (US-Beo, US-Bes, and US-Brw), one in Atqasuk (US-Atq) and one in Ivotuk (US-Ivo), for a period of 2013–2015. Three footprint algorithms were the homogenous footprint (HF) that assumes even contribution of all grid cells, the gradient footprint (GF) that assumes gradually declining contribution from center grid cells to edges, and the dynamic footprint (DF) that considers the impacts of wind and heterogeneity of land surface. Simulated annual CH 4 flux was highly consistent with the EC measurements at US-Beo and US-Bes. In contrast, flux was overestimated at US-Brw, US-Atq, and US-Ivo due to the higher simulated CH 4 flux in early growing seasons. The simulated monthly CH 4 flux was consistent with EC measurements but with different accuracies among footprint algorithms. At US-Bes in September 2013, RMSE and NNSE were 0.002 μmol m −2  s −1 and 0.782 using the DF algorithm, but 0.007 μmol m −2  s −1 and 0.758 using HF and 0.007 μmol m −2  s −1 and 0.765 using GF, respectively. DF algorithm performed better than the HF and GF algorithms in capturing the temporal variation in daily CH 4 flux each month, while the model accuracy was similar among the three algorithms due to flat landscapes. Temporal variations in CH 4 flux during 2013–2015 were predominately explained by air temperature (67–74%), followed by precipitation (22–36%). Spatial heterogeneities in vegetation fraction and elevation dominated the spatial variations in CH 4 flux for all five tower domains despite relatively weak differences in simulated CH 4 flux among three footprint algorithms. The CLM-Microbe model can simulate CH 4 flux at both plot and landscape scales at a high temporal resolution, which should be applied to other landscapes. Integrating land surface models with an appropriate algorithm provides a powerful tool for upscaling CH 4 flux in terrestrial ecosystems. 
    more » « less
  2. Abstract Understanding tree transpiration variability is vital for assessing ecosystem water‐use efficiency and forest health amid climate change, yet most landscape‐level measurements do not differentiate individual trees. Using canopy temperature data from thermal cameras, we estimated the transpiration rates of individual trees at Harvard Forest and Niwot Ridge. PT‐JPL model was used to derive latent heat flux from thermal images at the canopy‐level, showing strong agreement with tower measurements (R2 = 0.70–0.96 at Niwot, 0.59–0.78 at Harvard at half‐hourly to monthly scales) and daily RMSE of 33.5 W/m2(Niwot) and 52.8 W/m2(Harvard). Tree‐level analysis revealed species‐specific responses to drought, with lodgepole pine exhibiting greater tolerance than Engelmann spruce at Niwot and red oak showing heightened resistance than red maple at Harvard. These findings show how ecophysiological differences between species result in varying responses to drought and demonstrate that these responses can be characterized by deriving transpiration from crown temperature measurements. 
    more » « less
  3. Eddy covariance measurements quantify the magnitude and temporal variability of land-atmosphere exchanges of water, heat, and carbon dioxide (CO 2 ) among others. However, they also carry information regarding the influence of spatial heterogeneity within the flux footprint, the temporally dynamic source/sink area that contributes to the measured fluxes. A 25 m tall eddy covariance flux tower in Central Illinois, USA, a region where drastic seasonal land cover changes from intensive agriculture of maize and soybean occur, provides a unique setting to explore how the organized heterogeneity of row crop agriculture contributes to observations of land-atmosphere exchange. We characterize the effects of this heterogeneity on latent heat ( LE ), sensible heat ( H ), and CO 2 fluxes ( F c ) using a combined flux footprint and eco-hydrological modeling approach. We estimate the relative contribution of each crop type resulting from the structured spatial organization of the land cover to the observed fluxes from April 2016 to April 2019. We present the concept of a fetch rose, which represents the frequency of the location and length of the prevalent upwind distance contributing to the observations. The combined action of hydroclimatological drivers and land cover heterogeneity within the dynamic flux footprint explain interannual flux variations. We find that smaller flux footprints associated with unstable conditions are more likely to be dominated by a single crop type, but both crops typically influence any given flux measurement. Meanwhile, our ecohydrological modeling suggests that land cover heterogeneity leads to a greater than 10% difference in flux magnitudes for most time windows relative to an assumption of equally distributed crop types. This study shows how the observed flux magnitudes and variability depend on the organized land cover heterogeneity and is extensible to other intensively managed or otherwise heterogeneous landscapes. 
    more » « less
  4. Abstract This study developed and evaluated 30‐m daily evapotranspiration (ET) estimates using the Priestley‐Taylor Jet Propulsion Laboratory (PT‐JPL) model with ECOSTRESS, Moderate MODIS, harmonized Landsat Sentinel‐2 (HLS) imagery, ERA5‐Land reanalysis, and eddy covariance measurements. The new daily 30‐m ET showed significantly improved performance (overall,r = 0.8, RMSE = 1.736, KGE = 0.466) at 145 EC sites over contiguous United States compared to the current 70‐m ECOSTRESS ET (overall,r = 0.485, RMSE = 4.696, KGE = −0.841). A deep neural network postprocessing model trained with ET measurements from EC sites further improved the performance on test sites that were not used for model training (overall,r = 0.842, RMSE = 0.88, KGE = 0.792). The 30‐m ET estimation biases were significantly related to the biases in the upwelling longwave (RUL) and downwelling shortwave radiation (RDS) inputs, with ET estimates driven by MODIS radiation showing higher biases compared to those driven by ERA5‐Land radiation. The error diagnosis using random forest indicates that ET biases tend to be larger under higher ET estimates, andRULandRDSwere the primary contributors to the high bias at the higher ET ranges, with partial dependence plots revealing that the estimation biases tend to be higher under more humid environment, denser vegetation covers, and high net radiation conditions. In conclusion, higher spatial resolution satellite imagery of vegetation characteristics and higher temporal resolution radiation data, combined with continent‐wide EC measurements and deep learning, provided substantial added value for improving ET estimations at the field scale (30‐m). 
    more » « less
  5. Surface albedo is of crucial interest in land–climate interaction studies, since it is a key parameter that affects the Earth’s radiation budget. The temporal and spatial variation of surface albedo can be retrieved from conventional satellite observations after a series of processes, including atmospheric correction to surface spectral bi-directional reflectance factor (BRF), bi-directional reflectance distribution function (BRDF) modelling using these BRFs, and, where required, narrow-to-broadband albedo conversions. This processing chain introduces errors that can be accumulated and then affect the accuracy of the retrieved albedo products. In this study, the albedo products derived from the multi-angle imaging spectroradiometer (MISR), moderate resolution imaging spectroradiometer (MODIS) and the Copernicus Global Land Service (CGLS), based on the VEGETATION and now the PROBA-V sensors, are compared with albedometer and upscaled in situ measurements from 19 tower sites from the FLUXNET network, surface radiation budget network (SURFRAD) and Baseline Surface Radiation Network (BSRN) networks. The MISR sensor onboard the Terra satellite has 9 cameras at different view angles, which allows a near-simultaneous retrieval of surface albedo. Using a 16-day retrieval algorithm, the MODIS generates the daily albedo products (MCD43A) at a 500-m resolution. The CGLS albedo products are derived from the VEGETATION and PROBA-V, and updated every 10 days using a weighted 30-day window. We describe a newly developed method to derive the two types of albedo, which are directional hemispherical reflectance (DHR) and bi-hemispherical reflectance (BHR), directly from three tower-measured variables of shortwave radiation: downwelling, upwelling and diffuse shortwave radiation. In the validation process, the MISR, MODIS and CGLS-derived albedos (DHR and BHR) are first compared with tower measured albedos, using pixel-to-point analysis, between 2012 to 2016. The tower measured point albedos are then upscaled to coarse-resolution albedos, based on atmospherically corrected BRFs from high-resolution Earth observation (HR-EO) data, alongside MODIS BRDF climatology from a larger area. Then a pixel-to-pixel comparison is performed between DHR and BHR retrieved from coarse-resolution satellite observations and DHR and BHR upscaled from accurate tower measurements. The experimental results are presented on exploring the parameter space associated with land cover type, heterogeneous vs. homogeneous and instantaneous vs. time composite retrievals of surface albedo. 
    more » « less