skip to main content


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. Abstract

    Airborne light detection and ranging is an emerging measurement tool for snowpack estimation, and data are now emerging to better assess multiscale snow depth patterns. We used airborne light detection and ranging measurements from four sites in the southern Sierra Nevada to determine how snow depth varies with canopy structure and the interactions between canopies and terrain. We processed the point clouds into snow depth rasters at 0.5×0.5‐m2resolution and performed statistical analysis on the processed snow depth data, terrain attributes, and vegetation attributes, including the individual tree bole locations, canopy crown area, and canopy height. We studied the snow depth at such a fine scale due in part to the spatial heterogeneity introduced by canopy interception and enhanced melting caused by tree trunks in forested areas. We found that the dominant direction of a tree well, the area around the tree bole that has shallower snowpack, is correlated with the local aspect of the terrain, and the gradient of the snow surface in a tree well is correlated with the tree's crown area. The regression‐tree based XGBoost model was fitted with the topographic variables and canopy variables, and about 71% of snow depth variability can be explained by the model.

     
    more » « less
  3. 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
  4. 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
  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