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.


Title: High-Resolution Snow-Covered Area Mapping in Forested Mountain Ecosystems Using PlanetScope Imagery
Improving high-resolution (meter-scale) mapping of snow-covered areas in complex and forested terrains is critical to understanding the responses of species and water systems to climate change. Commercial high-resolution imagery from Planet Labs, Inc. (Planet, San Francisco, CA, USA) can be used in environmental science, as it has both high spatial (0.7–3.0 m) and temporal (1–2 day) resolution. Deriving snow-covered areas from Planet imagery using traditional radiometric techniques have limitations due to the lack of a shortwave infrared band that is needed to fully exploit the difference in reflectance to discriminate between snow and clouds. However, recent work demonstrated that snow cover area (SCA) can be successfully mapped using only the PlanetScope 4-band (Red, Green, Blue and NIR) reflectance products and a machine learning (ML) approach based on convolutional neural networks (CNN). To evaluate how additional features improve the existing model performance, we: (1) build on previous work to augment a CNN model with additional input data including vegetation metrics (Normalized Difference Vegetation Index) and DEM-derived metrics (elevation, slope and aspect) to improve SCA mapping in forested and open terrain, (2) evaluate the model performance at two geographically diverse sites (Gunnison, Colorado, USA and Engadin, Switzerland), and (3) evaluate the model performance over different land-cover types. The best augmented model used the Normalized Difference Vegetation Index (NDVI) along with visible (red, green, and blue) and NIR bands, with an F-score of 0.89 (Gunnison) and 0.93 (Engadin) and was found to be 4% and 2% better than when using canopy height- and terrain-derived measures at Gunnison, respectively. The NDVI-based model improves not only upon the original band-only model’s ability to detect snow in forests, but also across other various land-cover types (gaps and canopy edges). We examined the model’s performance in forested areas using three forest canopy quantification metrics and found that augmented models can better identify snow in canopy edges and open areas but still underpredict snow cover under forest canopies. While the new features improve model performance over band-only options, the models still have challenges identifying the snow under trees in dense forests, with performance varying as a function of the geographic area. The improved high-resolution snow maps in forested environments can support studies involving climate change effects on mountain ecosystems and evaluations of hydrological impacts in snow-dominated river basins.  more » « less
Award ID(s):
1947875 2117834
PAR ID:
10406223
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Remote Sensing
Volume:
14
Issue:
14
ISSN:
2072-4292
Page Range / eLocation ID:
3409
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Mountain snowpack provides critical water resources for forest and meadow ecosystems that are experiencing rapid change due to global warming. An accurate characterization of snowpack heterogeneity in these ecosystems requires snow cover observations at high spatial resolutions, yet most existing snow cover datasets have a coarse resolution. To advance our observation capabilities of snow cover in meadows and forests, we developed a machine learning model to generate snow-covered area (SCA) maps from PlanetScope imagery at about 3-m spatial resolution. The model achieves a median F1 score of 0.75 for 103 cloud-free images across four different sites in the Western United States and Switzerland. It is more accurate (F1 score = 0.82) when forest areas are excluded from the evaluation. We further tested the model performance across 7,741 mountain meadows at the two study sites in the Sierra Nevada, California. It achieved a median F1 score of 0.83, with higher accuracy for larger and simpler geometry meadows than for smaller and more complexly shaped meadows. While mapping SCA in regions close to or under forest canopy is still challenging, the model can accurately identify SCA for relatively large forest gaps (i.e., 15m < DCE < 27m), with a median F1 score of 0.87 across the four study sites, and shows promising accuracy for areas very close (>10m) to forest edges. Our study highlights the potential of high-resolution satellite imagery for mapping mountain snow cover in forested areas and meadows, with implications for advancing ecohydrological research in a world expecting significant changes in snow. 
    more » « less
  2. Accurate measurements of terrain elevation are crucial for many ecological applications. In this study, we sought to assess new global three-dimensional Earth observation data acquired by the spaceborne Light Detection and Ranging (LiDAR) missions Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) and Global Ecosystem Dynamics Investigation (GEDI). For this, we examined the “ATLAS/ICESat-2 L3A Land and Vegetation Height”, version 5 (20 × 14 m and 100 × 14 m segments) and the “GEDI Level 2A Footprint Elevation and Height Metrics”, version 2 (25 m circle). We conducted our analysis across four land cover classes (bare soil, herbaceous, forest, savanna), and six forest types (temperate broad-leaved, temperate needle-leaved, temperate mixed, tropical upland, tropical floodplain, and tropical secondary forest). For assessment of terrain elevation estimates from spaceborne LiDAR data we used high resolution airborne data. Our results indicate that both LiDAR missions provide accurate terrain elevation estimates across different land cover classes and forest types with mean error less than 1 m, except in tropical forests. However, using a GEDI algorithm with a lower signal end threshold (e.g., algorithm 5) can improve the accuracy of terrain elevation estimates for tropical upland forests. Specific environmental parameters (terrain slope, canopy height and canopy cover) and sensor parameters (GEDI degrade flags, terrain estimation algorithm; ICESat-2 number of terrain photons, terrain uncertainty) can be applied to improve the accuracy of ICESat-2 and GEDI-based terrain estimates. Although the goodness-of-fit statistics from the two spaceborne LiDARs are not directly comparable since they possess different footprint sizes (100 × 14 m segment or 20 × 14 m segment vs. 25 m circle), we observed similar trends on the impact of terrain slope, canopy cover and canopy height for both sensors. Terrain slope strongly impacts the accuracy of both ICESat-2 and GEDI terrain elevation estimates for both forested and non-forested areas. In the case of GEDI the impact of slope is, however, partly caused by horizontal geolocation error. Moreover, dense canopies (i.e., canopy cover higher than 90%) affect the accuracy of spaceborne LiDAR terrain estimates, while canopy height does not, when considering samples over flat terrains. Our analysis of the accuracy and precision of current versions of spaceborne LiDAR products for different vegetation types and environmental conditions provides insights on parameter selection and estimated uncertainty to inform users of these key global datasets. 
    more » « less
  3. Abstract Located at northern latitudes and subject to large seasonal temperature fluctuations, boreal forests are sensitive to the changing climate, with evidence for both increasing and decreasing productivity, depending upon conditions. Optical remote sensing of vegetation indices based on spectral reflectance offers a means of monitoring vegetation photosynthetic activity and provides a powerful tool for observing how boreal forests respond to changing environmental conditions. Reflectance‐based remotely sensed optical signals at northern latitude or high‐altitude regions are readily confounded by snow coverage, hampering applications of satellite‐based vegetation indices in tracking vegetation productivity at large scales. Unraveling the effects of snow can be challenging from satellite data, particularly when validation data are lacking. In this study, we established an experimental system in Alberta, Canada including six boreal tree species, both evergreen and deciduous, to evaluate the confounding effects of snow on three vegetation indices: the normalized difference vegetation index (NDVI), the photochemical reflectance index (PRI), and the chlorophyll/carotenoid index (CCI), all used in tracking vegetation productivity for boreal forests. Our results revealed substantial impacts of snow on canopy reflectance and vegetation indices, expressed as increased albedo, decreased NDVI values and increased PRI and CCI values. These effects varied among species and functional groups (evergreen and deciduous) and different vegetation indices were affected differently, indicating contradictory, confounding effects of snow on these indices. In addition to snow effects, we evaluated the contribution of deciduous trees to vegetation indices in mixed stands of evergreen and deciduous species, which contribute to the observed relationship between greenness‐based indices and ecosystem productivity of many evergreen‐dominated forests that contain a deciduous component. Our results demonstrate confounding and interacting effects of snow and vegetation type on vegetation indices and illustrate the importance of explicitly considering snow effects in any global‐scale photosynthesis monitoring efforts using remotely sensed vegetation indices. 
    more » « less
  4. Abstract Changes in vegetation productivity based on normalized difference vegetation index (NDVI) have been reported from Arctic regions. Most studies use very coarse spatial resolution remote sensing data that cannot isolate landscape level factors. For example, on Yamal Peninsula in West Siberia enhanced willow growth has been linked to widespread landslide activity, but the effect of landslides on regional NDVI dynamics is unknown. Here we apply a novel satellite-based NDVI analysis to investigate the vegetation regeneration patterns of active-layer detachments following a major landslide event in 1989. We analyzed time series data of Landsat and very high-resolution (VHR) imagery from QuickBird-2 and WorldView-2 and 3 characterizing a study area of ca. 35 km2. Landsat revealed that natural regeneration of low Arctic tundra progressed rapidly during the first two decades after the landslide event. However, during the past decade, the difference between landslide shear surfaces and surrounding areas remained relatively unchanged despite the advance of vegetation succession. Time series also revealed that NDVI generally declined since 2013 within the study area. The VHR imagery allowed detection of NDVI change ‘hot-spots’ that included temporary degradation of vegetation cover, as well as new and expanding thaw slumps, which were too small to be detected from Landsat satellite data. Our study demonstrates that landslides can have pronounced and long-lasting impacts on tundra vegetation. Thermokarst landslides and associated impacts on vegetation will likely become increasingly common in NW Siberia and other Arctic regions with continued warming. 
    more » « less
  5. Abstract Patterns of alpine plant productivity are extremely variable in space and time. Complex topography drives variations in temperature, wind, and solar radiation. Surface and subsurface flow paths route water between landscape patches. Redistribution of snow creates scour zones and deep drifts, which drives variation in water availability and growing season length. Hence, the distribution of snow likely plays a central role in patterns of alpine plant productivity. Given that these processes operate at sub‐1 m to sub‐10 m spatial scales and are dynamic across daily to weekly time scales, historical studies using manual survey techniques have not afforded a comprehensive assessment of the influence of snow distribution on plant productivity. To address this knowledge gap, we used weekly estimates of normalised difference vegetation index (NDVI), snow extent, and peak snow depth, acquired from drone surveys at 25 cm resolution. We derived six snowpack‐related and topographic variables that may influence vegetation productivity and analysed these with respect to the timing and magnitude of peak productivity. Peak NDVI and peak NDVI timing were most highly correlated with maximum snow depth, and snow‐off‐date. We observed up to a ~30% reduction in peak NDVI for areas with deeper and later snow cover, and a ~11‐day delay in the timing of peak NDVI in association with later snow‐off‐date. Our findings leverage a novel approach to quantify the importance of snow distribution in driving alpine vegetation productivity and provide a space for time proxy of potential changes in a warmer, lower snow future. 
    more » « less