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
Monitoring tropical forest succession at landscape scales despite uncertainty in Landsat time series
Forecasting rates of forest succession at landscape scales will aid global efforts to restore tree cover to millions of hectares of degraded land. While optical satellite remote sensing can detect regional land cover change, quantifying forest structural change is challenging. We developed a state‐space modeling framework that applies Landsat satellite data to estimate variability in rates of natural regeneration between sites in a tropical landscape. Our models work by disentangling measurement error in Landsat‐derived spectral reflectance from process error related to successional variability. We applied our modeling framework to rank rates of forest succession between 10 naturally regenerating sites in Southwestern Panama from about 2001 to 2015 and tested how different models for measurement error impacted forecast accuracy, ecological inference, and rankings of successional rates between sites. We achieved the greatest increase in forecasting accuracy by adding intra‐annual phenological variation to a model based on Landsat‐derived normalized difference vegetation index (NDVI). The best‐performing model accounted for inter‐ and intra‐annual noise in spectral reflectance and translated NDVI to canopy height via Landsat–lidar fusion. Modeling forest succession as a function of canopy height rather than NDVI also resulted in more realistic estimates of forest state during early succession, including greater confidence in rank order of successional rates between sites. These results establish the viability of state‐space models to quantify ecological dynamics from time series of space‐borne imagery. State‐space models also provide a statistical approach well‐suited to fusing high‐resolution data, such as airborne lidar, with lower‐resolution data that provides better temporal and spatial coverage, such as the Landsat satellite record. Monitoring forest succession using satellite imagery could play a key role in achieving global restoration targets, including identifying sites that will regain tree cover with minimal intervention.
more » « less- NSF-PAR ID:
- 10378163
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Ecological Applications
- Volume:
- 31
- Issue:
- 1
- ISSN:
- 1051-0761
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Many studies have used time series of satellite-derived vegetation indices to identify so-called greening and browning trends across the northern high-latitudes and to suggest that the productivity of Arctic-Boreal ecosystems is changing in response to climate forcing at local and continental scales. However, disturbances that alter land cover are prevalent in Arctic-Boreal ecosystems, and changes in Arctic-Boreal land cover, which complicate interpretation of trends in vegetation indices, have mostly been ignored in previous studies. Here we use a new land cover change dataset derived from Landsat imagery to explore the extent to which land cover and land cover change influence trends in the normalized difference vegetation index (NDVI) over a large (3.76 M km2) area of NASA’s Arctic Boreal Vulnerability Experiment, which spans much of northwestern Canada and Alaska. Between 1984 and 2012, 21.2% of the study domain experienced land cover change and 42.7% had significant NDVI trends. Land cover change occurred in 27.6% of locations with significant NDVI trends during this period and resulted in greening and browning rates 48%–128% higher than in areas of stable land cover. While the majority of land cover change areas experienced significant NDVI trends, more than half of areas with stable land cover did not. Further, the extent and magnitude of browning and greening trends varied substantially as a function of land cover class and land cover change type. Forest disturbance from fire and timber harvest drove over one third of statistically significant NDVI trends and created complex mosaics of recent forest loss (as browning) and post-disturbance recovery (as greening) at both landscape and continental scale. Our results demonstrate the importance of land cover changes in highly disturbed high-latitude ecosystems for interpreting trends of NDVI and productivity across multiple spatial scales.
-
more » « less
Abstract Here, we evaluate the ecosystem functioning and the ecosystems services supply of different vegetation types (grasslands, shrublands and woodlands) under contrasting management regimes by comparing a protected area with the surrounding landscape, which has been subjected to human disturbance in the Eastern Hills of Uruguay. We propose, based on functional attributes and vegetation physiognomy, a State and Transition Model for the dynamics of the grassland–woodland mosaic. We used remote sensing techniques to: (i) develop a land‐cover map of the study area based on supervised Landsat imagery classification, and (ii) compare attributes of the ecosystem functioning (productivity and seasonality) and service supply derived from the Normalized Difference Vegetation Index (NDVI) images provided by the moderate resolution imaging spectroradiometer (MODIS) sensor. The land‐cover map showed that grasslands and shrublands were the most extensive land covers in the study area. These vegetation types presented higher productivity, seasonality and ecosystem service supply, outside the protected area than inside it. On the other hand, woodlands showed higher productivity, ecosystem service supply and lower seasonality inside the protected area than outside of it. Two axes represented the grassland–woodland mosaic dynamic: (i) the mean annual and (ii) the intra‐annual coefficient of variation of the NDVI. Our results highlight that conservation of grasslands, shrublands and woodlands require different management strategies based on particular disturbance regimes like moderate grazing and controlled burns. Moderate disturbances may help to preserve ecosystem services provisioning in grasslands and shrublands. On the contrary, woodland conservation requires a more rigorous regime of protection against disturbances.
-
ArcticDEM provides the public with an unprecedented opportunity to access very high-spatial resolution digital elevation models (DEMs) covering the pan-Arctic surfaces. As it is generated from stereo-pairs of optical satellite imagery, ArcticDEM represents a mixture of a digital surface model (DSM) over a non-ground areas and digital terrain model (DTM) at bare grounds. Reconstructing DTM from ArcticDEM is thus needed in studies requiring bare ground elevation, such as modeling hydrological processes, tracking surface change dynamics, and estimating vegetation canopy height and associated forest attributes. Here we proposed an automated approach for estimating DTM from ArcticDEM in two steps: (1) identifying ground pixels from WorldView-2 imagery using a Gaussian mixture model (GMM) with local refinement by morphological operation, and (2) generating a continuous DTM surface using ArcticDEMs at ground locations and spatial interpolation methods (ordinary kriging (OK) and natural neighbor (NN)). We evaluated our method at three forested study sites characterized by different canopy cover and topographic conditions in Livengood, Alaska, where airborne lidar data is available for validation. Our results demonstrate that (1) the proposed ground identification method can effectively identify ground pixels with much lower root mean square errors (RMSEs) (<0.35 m) to the reference data than the comparative state-of-the-art approaches; (2) NN performs more robustly in DTM interpolation than OK; (3) the DTMs generated from NN interpolation with GMM-based ground masks decrease the RMSEs of ArcticDEM to 0.648 m, 1.677 m, and 0.521 m for Site-1, Site-2, and Site-3, respectively. This study provides a viable means of deriving high-resolution DTM from ArcticDEM that will be of great value to studies focusing on the Arctic ecosystems, forest change dynamics, and earth surface processes.more » « less
-
more » « less
Abstract In post‐fire Siberian larch forests, where tree density can vary within a burn perimeter, shrubs constitute a substantial portion of the vegetation canopy. Leaf area index (LAI), defined as the one‐sided total green leaf area per unit ground surface area, is useful for characterizing variation in plant canopies. We estimated LAI with allometry for trees and tall shrubs (>0.5 and <1.5 m) across 26 sites with varying tree stem density (0.05–3.3 stems/m2) and canopy cover (4.6%–76.9%) in a uniformly‐aged mature Siberian larch forest that regenerated following a fire ∼75 years ago. We investigated relationships between tree density, tree LAI, and tall shrub LAI, and between LAI and satellite observations of Normalized Difference and Enhanced Vegetation Indices (NDVI and EVI). Across the density gradient, tree LAI increases with increasing tree density, while tall shrub LAI decreases, exhibiting no patterns in combined tree‐shrub LAI. We also found significant positive relationships between tall shrub LAI and NDVI/EVI from PlanetScope and Landsat imagery. These findings suggest that tall shrubs compensate for lower tree LAI in tree canopy gaps, forming a canopy with contiguous combined tree‐shrub LAI across the density gradient. Our findings suggest that NDVI and EVI are more sensitive to variation in tall shrub canopies than variation in tree canopies or combined tree‐shrub canopies in these ecosystems. The results improve our understanding of the relationships between forest density and tree and shrub leaf area and have implications for interpreting spatial variability in LAI, NDVI, and EVI in Siberian boreal forests.
