USGS ScienceSearch

USGS · 70230528

Mapping actual evapotranspiration using Landsat for the conterminous United States: Google Earth Engine implementation and assessment of the SSEBop model

Abstract

The estimation and mapping of actual evapotranspiration (ETa) is an active area of applied research in the fields of agriculture and water resources. Thermal remote sensing-based methods, using coarse resolution satellites, have been successful at estimating ETa over the conterminous United States (CONUS) and other regions of the world. In this study, we present CONUS-wide ETa from Landsat thermal imagery-using the Operational Simplified Surface Energy Balance (SSEBop) model in the Google Earth Engine (GEE) cloud computing platform. Over 150,000 Landsat satellite images were used to produce 10 years of annual ETa (2010–2019) at unprecedented scale. The accuracy assessment of the SSEBop results included point-based evaluation using monthly Eddy Covariance (EC) data from 25 AmeriFlux stations as well as basin-scale comparison with annual Water Balance ETa (WBET) for more than 1000 sub-basins. Evaluations using EC data showed generally mixed performance with weaker (R 2 < 0.6) correlation on sparsely vegetated surfaces such as grasslands or woody savanna and stronger correlation (R 2 > 0.7) over well-vegetated surfaces such as croplands and forests, but location-specific conditions rather than cover type were attributed to the variability in accuracy. Croplands performed best with R 2 of 0.82, root mean square error of 29 mm/month, and average bias of 12%. The WBET evaluation indicated that the SSEBop model is strong in explaining the spatial variability (up to R 2 > 0.90) of ETa across large basins, but it also identified broad hydro-climatic regions where the SSEBop ETa showed directional biases, requiring region-specific model parameter improvement and/or bias correction with an overall 7% bias nationwide. Annual ETa anomalies over the 10-year period captured widely reported drought-affected regions, for the most part, in different parts of the CONUS, indicating their potential applications for mapping regional- and field-scale drought and fire effects. Due to the coverage of the Landsat Path/Row system, the availability of cloud-free image pixels ranged from less than 12 (mountainous cloud-prone regions and U.S. Northeast) to more than 60 (U.S. Southwest) per year. However, this study reinforces a promising application of Landsat satellite data with cloud-computing for quick and efficient mapping of ETa for agricultural and water resources assessments at the field scale.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 25.08° to 49.38905° latitude; -124.68721° to -66.96466° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Gabriel B. Senay, MacKenzie Friedrichs, Charles Morton, Gabriel Edwin Lee Parrish, Matthew Schauer, Kul Bikram Khand, Stefanie Kagone, Olena Boiko, Justin Huntington. 2022. Mapping actual evapotranspiration using Landsat for the conterminous United States: Google Earth Engine implementation and assessment of the SSEBop model. https://doi.org/10.1016/j.rse.2022.113011

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

On-demand global Landsat evapotranspiration product: Development, evaluation, and dissemination

Global actual evapotranspiration (ET) is one of the essential climate variables needed to understand and manage the relationships among food, energy, and water resources. The U.S. Geological Survey Earth Resources Observation and Science (EROS) Center launched a provisional ET product in 2020, offering on-demand, field-scale global coverage derived from Landsat data through the EROS Science Processing Architecture (ESPA) platform. The ESPA interface provides ET data for cloud-free Landsat overpasses starting in 1982 with Landsat 4 through the current Landsat 9. The ET data are delivered as a Provisional Level-3 Science product created using the Operational Simplified Surface Energy Balance (SSEBop) model. Landsat surface temperature and reference ET are the main model drivers along with vegetation index and net radiation for model parameterization. A large volume of Landsat-based ET orders (e.g., over 1,200,000 images from June 2020 through December 2025) around the world indicate increasing awareness and application of the ET data. The ESPA platform enables land and water resource managers and researchers to access a first-order ET product without requiring advanced knowledge of remote sensing technology or evapotranspiration modeling. We present the methodology and workflow of the on-demand Landsat ET product and its performance evaluations over diverse hydro-climatic settings. The product can help estimate field-scale consumptive water use and thus quickly and consistently assess historical water use, allocation, and budget to inform water management under changing environments. Future ET data aggregated to monthly and seasonal time scales are expected to enhance integration with decision-making tools and procedures.

Remote Sensing of Environment

A framework for integrating spatiotemporal deep learning methods with landsat for annual land cover and impervious surface mapping

Land cover information is essential for understanding Earth’s surface dynamics and how vegetation, water, soil, climate, and terrain interact. The National Land Cover Database (NLCD) has been the authoritative source for consistent U.S. land cover mapping. To extend NLCD’s temporal resolution and reduce production latency, we developed the Land Cover Artificial Mapping System (LCAMS)—a prototype spatiotemporal deep learning framework piloted as the foundation for the new Annual NLCD. LCAMS builds on concepts from legacy NLCD and the U.S. Geological Survey Land Change Monitoring, Assessment, and Projection (LCMAP) initiatives. It employs a loosely coupled two-stage architecture consisting of independent but functionally interdependent spatial and temporal models. Spatial models extract per-year information from Landsat data, while the temporal models refine the spatial outputs to enforce inter-annual consistency—critical for reliable land change monitoring. LCAMS produces annual 30 m resolution land cover and impervious surface outputs, with region-specific fine-tuning to generalize across diverse landscapes and temporal dynamics. Validation was conducted using an independent dataset of 1925 randomly sampled plots from five U.S. Landsat Analysis Ready Data (ARD) tiles spanning 1985-2021, selected for spatial and temporal variability. This dataset was used consistently to evaluate LCAMS, Legacy NLCD, and LCMAP. Using the NLCD legend, LCAMS achieved 72.1 ± 1.60% overall agreement, compared to 71.1 ± 1.7% agreement for Legacy NLCD. Using the LCMAP legend, LCAMS achieved 83.4 ± 1.22% agreement, compared to 84.6 ± 1.11% agreement for LCMAP. Overall, LCAMS delivers comparable accuracy while offering higher thematic resolution, longer temporal coverage, and automated production of annual 30 m CONUS land cover.

Remote Sensing of Environment

Towards global mapping of dynamic surface water extents using Sentinel-1 SAR data

We introduce a fully automated and scalable method for mapping surface water extents from single-acquisition Sentinel-1 synthetic aperture radar (SAR) imagery. This approach integrates adaptive thresholding of radiometric terrain-corrected SAR backscatter data, fuzzy-logic classification, region growing, dark land estimation, and a bimodality test to minimize false positives in low-backscattering areas and false negatives in high-backscattering areas. By combining these steps, the algorithm achieves classification accuracies exceeding 85% in detecting surface water extents across diverse environmental conditions. Accuracy was first assessed at meter scale using 52 PlanetScope scenes acquired worldwide in September–October 2019; the algorithm achieved 93% overall accuracy, 86% user's accuracy, and 94% producer's accuracy. Global robustness was then evaluated by processing every Sentinel-1 acquisition from 1 to 12 November 2023 and cross-comparing the resulting maps with 6561 temporally matched observational products for end-users from remote sensing analysis (OPERA) dynamic surface water extent from Harmonized Landsat and Sentinel-2 (DSWx-HLS) products. This large-scale test yielded 90% user's and 94% producer's accuracies, confirming reliable performance at continental extent. Additional case studies demonstrate the algorithm's ability to handle surface water extent in sand-dominated deserts, to track seasonal amplitude in Folsom Lake (California), drought-induced loss in Cerro Prieto Reservoir (Mexico), and rapid filling of the Grand Ethiopian Renaissance Dam. These results show that the method scales across local to global domains and maintains high accuracy, providing a practical tool for near-real-time monitoring of floods, droughts, and water-resource management. Because the approach is sensor-agnostic, it can be ported to forthcoming L- and S-band missions such as NASA-ISRO synthetic aperture radar (NISAR), broadening its applicability to future hydrologic observations.

Remote Sensing of Environment