Integrating stream gage data and Landsat imagery to complete time-series of surface water extents in Central Valley, California

Jessica J. Walker; Christopher E. Soulard; Roy E. Petrakis

doi:10.1016/j.jag.2019.101973

Home/Archive/USGS

USGS · 70205929

Integrating stream gage data and Landsat imagery to complete time-series of surface water extents in Central Valley, California

Jessica J. Walker·Christopher E. Soulard·Roy E. Petrakis

View original source ↗Export metadata ↓

Abstract

Accurate monitoring of surface water location and extent is critical for the management of diverse water resource phenomena. The multi-decadal archive of Landsat satellite imagery is punctuated by missing data due to cloud cover during acquisition times, hindering the assembly of a continuous time series of inundation dynamics. This study investigated whether streamflow volume measurements could be integrated with satellite data to fill gaps in monthly surface water chronologies for the Central Valley region of California, USA, from 1984 to 2015. We aggregated measurements of maximum monthly water extent within each of the study area’s 50 8-digit hydrologic unit code [HUC] watersheds from two Landsat-derived datasets: the European Commission’s Joint Research Centre (JRC) Monthly Water History and the U.S. Geological Survey Dynamic Surface Water Extent (DSWE). We calculated Spearman rank correlation coefficients between water extent values in each HUC and streamflow discharge data. Linear regression fits of the water extent/streamflow data pairs with the highest correlations served as the basis for interpolation of missing imagery surface water values on a HUC-wise basis. Results show strong (ρ > 0.7) maximum correlations in 11 (22.4%) and 25 (51.0%) HUCs for the DSWE and JRC time series, respectively, when comparisons were restricted to imagery and gages co-located in each HUC. Strong maximum correlations occurred in 39 (79.6%; DSWE) and 42 (85.7%; JRC) HUCs when imagery was paired with discharge data from any study area gage, providing a solid basis for reconstruction of water extent values. We generated continuous time series of 30+ years in 35 HUCs, demonstrating that this technique can provide quantitative estimates of historical surface water extents and elucidate flooding or drought events over the period of data collection. Results of a non-parametric trend analysis of the long-term time series on an annual, seasonal, and monthly basis varied among HUCs, though most trends indicate an increase in surface water over the past 30 years.

Explore related subjects

International Journal of Applied Earth Observation and Geoinformation Article

Source-reported bounding extent: 34.70549341022544° to 41.5579215778042° latitude; -123.53027343749999° to -116.65283203124999° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Jessica J. Walker, Christopher E. Soulard, Roy E. Petrakis. 2019. Integrating stream gage data and Landsat imagery to complete time-series of surface water extents in Central Valley, California. https://doi.org/10.1016/j.jag.2019.101973

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

KEEP EXPLORING

Related USGS reports

USGS2026

Satellite embeddings for crop type classification: A comparative examination

Embedding datasets encode complex relationships among multiple sources of Earth observation data into a compact format. Here, we evaluated the utility of a 10-m global Satellite Embedding product (SE) for classifying crop types in central California for the year 2020. We compared the classification accuracy of a random forest model based exclusively on the SE layer to an existing random forest model with multiple imagery inputs. Our results showed the SE-based classification had higher agreement with the reference dataset (California Department of Water Resources crop map) than the classification based on Landsat and National Agricultural Imagery Program inputs (94.7% versus. 91.9% overall accuracy, respectively). The performance of individual crop types was consistent across models, ranging from high agreement for rice (98.4% versus 98% accuracy) to lower agreement for pasture, grain, and fallow/young perennial classes (< 65% accuracy in both models). The SE-based workflow used three times less cloud-based computational resources and represented substantial savings of predictor development time. Geospatial embedding products can aid classification efforts by reducing predictor development time and processing demands while maintaining classification accuracy.

California↗

USGS2026

The impacts of cover crop biomass on satellite-based detectability of cover crops in Maryland

Cover crop adoption in the U.S. has increased over the past decades, increasing the need to quantify their performance and environmental benefits. While remote sensing (RS)-based approaches for detecting cover crop presence have been developed, there has been limited research on how cover crops with varied biomass and management practices influence detectability. Using unique field-level cover crop presence and biomass datasets in Maryland, U.S., we investigated how RS-based detectability changes for cover crops with varied aboveground biomass, planting, and termination dates. Specifically, we proposed a time-integrated satellite-based greenness feature from Harmonized Landsat-8 and Sentinel-2 (HLS) time series from 2017 to 2021 to estimate biomass of cover crops and evaluated their detectability using a phenology-based cover crop detection framework. The impacts of cover crop planting and termination dates on cover crop biomass were also analyzed. Our results demonstrated that Normalized Difference Vegetation Index (NDVI) estimated biomass with higher accuracy compared to other vegetation indices, and the time-integrated model estimated biomass with higher accuracy than the single-date “snapshot” linear model (R2 from 0.53 to 0.66 and RMSE from 922 kg/ha to 747 kg/ha). While the snapshot models were species sensitive, the time-integrated models showed strong robustness across different species. Detectability increased with cover crop biomass, as detected cover crops averaged 963.3 ±719.5 kg/ha compared to 297.2 ± 209.0 kg/ha for non-detected cover crops. Detection accuracy reached 96.1% for fields exceeding 500 kg/ha, compared with an overall accuracy of 62.7%. Earlier planting and later termination increased biomass and detectability, with biomass rising by 4.14 kg/ha/day (p < 0.01). This study demonstrates how management practices affect cover crop biomass and detectability via satellite time series and provides insights that can inform management of cover crops and monitoring of their effects on agroecosystems.

Maryland↗

USGS2026

Aligning legacy NLCD land cover maps based on Landsat Collection 1 to Collection 2

The transition from Landsat Collection 1 to Collection 2 introduced significant improvements in radiometric and geometric accuracy. However, the improvements cause location misalignment between the existing Landsat-derived land cover products and the new collection. The legacy National Land Cover Database (NLCD) has been used as a cornerstone land cover source for a variety of research. Therefore, a method aligning the legacy NLCD product to Collection 2 is required to ensure its continuity and consistency of service. We developed a strategy to not only align legacy NLCD to match new Collection 2 geometric locations but also improve land cover labeling in the region that was affected by the geometric shifts. The method identifies boundary pixels of homogeneous land cover patches as potential problem areas that are likely impacted by geometric shifts and generates candidate labels from 3 × 3 window with the target pixel at the center and segmentation-derived majority label. Standard phenology patterns of each candidate land cover type are established based on the random samples except boundary pixels within a 1000-pixels × 1000-pixels processing window region. The phenological distance to each standard land cover type pattern is calculated through a penalty dynamic time warping (DTW) method for each target pixel in the boundary region. Finally, the method determines the most suitable label based on the phenological distance from the candidate labels. Both visual and accuracy assessment results demonstrate that the alignment preserves the overall land cover patterns in the original legacy NLCD product while reducing the spatial discrepancies between the Landsat Collection 2 and land cover. In addition, it enhances the accuracy of land cover labeling of boundary pixels. The overall accuracy (OA) was increased by 7% in the land cover boundary regions after alignment. The quality and confusion matrix comparison between the alignment results and the original legacy NLCD confirm the reliability of the method. Our alignment method has the potential to serve as a framework for aligning other Landsat-derived land cover products to future collections.

conterminous United States↗