Sentinel-1 Shadows Used to Quantify Canopy Loss from Selective Logging in Gabon

Abstract

Selective logging is a major cause of forest degradation in the tropics, but its precise scale, location and timing are not known as wide-area, automated remote sensing methods are not yet available at this scale. This limits the abilities of governments to police illegal logging, or monitor (and thus receive payments for) reductions in degradation. Sentinel-1, a C-band Synthetic Aperture Radar satellite mission with a 12-day repeat time across the tropics, is a promising tool for this due to the known appearance of shadows in images where canopy trees are removed. However, previous work has relied on optical satellite data for calibration and validation, which has inherent uncertainties, leaving unanswered questions about the minimum magnitude and area of canopy loss this method can detect. Here, we use a novel bi-temporal LiDAR dataset in a forest degradation experiment in Gabon to show that canopy gaps as small as 0.02 ha (two 10 m × 10 m pixels) can be detected by Sentinel-1. The accuracy of our algorithm was highest when using a timeseries of 50 images over 20 months and no multilooking. With these parameters, canopy gaps in our study site were detected with a false alarm rate of 6.2%, a missed detection rate of 12.2%, and were assigned disturbance dates that were a good qualitative match to logging records. The presence of geolocation errors and false alarms makes this method unsuitable for confirming individual disturbances. However, we found a linear relationship (r2=0.74) between the area of detected Sentinel-1 shadow and LiDAR-based canopy loss at a scale of 1 hectare. By applying our method to three years’ worth of imagery over Gabon, we produce the first national scale map of small-magnitude canopy cover loss. We estimate a total gross canopy cover loss of 0.31 Mha, or 1.3% of Gabon’s forested area, which is a far larger area of change than shown in currently available forest loss alert systems using Landsat (0.022 Mha) and Sentinel-1 (0.019 Mha). Our results, which are made accessible through Google Earth Engine, suggest that this approach could be used to quantify the magnitude and timing of degradation more widely across tropical forests.