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

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

An Effective Method for InSAR Mapping of Tropical Forest Degradation in Hilly Areas

Current satellite remote sensing methods struggle to detect and map forest degradation, which is a critical issue as it is likely a major and growing source of carbon emissions and biodiveristy loss. TanDEM-X InSAR phase height is a promising variable for measuring forest disturbances, as it is closely related to the mean canopy height, and thus should decrease if canopy trees are removed. However, previous research has focused on relatively flat terrains, despite the fact that much of the world's remaining tropical forests are found in hilly areas, and this inevitably introduces artifacts in sideways imaging systems. In this paper, we find a relationship between InSAR phase height and aboveground biomass change in four selectively logged plots in a hilly region of central Gabon. We show that minimising multilooking prior to the calculation of InSAR phase height on a pixel-by-pixel basis. This shows that TanDEM-X InSAR can measure the magnitude of degradation, and that topographic effects can be mitigated if data from multiple SAR viewing geometries are available

Reliably Mapping Low-intensity Forest Disturbance Using Satellite Radar Data

In the last decades tropical forests have experienced increased fragmentation due to a global growing demand for agricultural and forest commodities. Satellite remote sensing offers a valuable tool for monitoring forest loss, thanks to the global coverage and the temporal consistency of the acquisitions. In tropical regions, C-band Synthetic Aperture Radar (SAR) data from the Sentinel-1 mission provides cloud-free and open imagery on a 6- or 12-day repeat cycle, offering the unique opportunity to monitor forest disturbances in a timely and continuous manner. Despite recent advances, mapping subtle forest losses, such as those due to small-scale and irregular selective logging, remains problematic. A Cumulative Sum (CuSum) approach has been recently proposed for forest monitoring applications, with preliminary studies showing promising results. Unfortunately, the lack of accurate in-situ measurements of tropical forest loss has prevented a full validation of this approach, especially in the case of low-intensity logging. In this study, we used high-quality field measurements from the tropical Forest Degradation Experiment (FODEX), combining unoccupied aerial vehicle (UAV) LiDAR, Terrestrial Laser Scanning (TLS), and field-inventoried data of forest structural change collected in two logging concessions in Gabon and Peru. The CuSum algorithm was applied to VV-polarized Sentinel-1 ground range detected (GRD) time series to monitor a range of canopy loss events, from individual tree extraction to forest clear cuts. We developed a single change metric using the maximum of the CuSum distribution, retrieving location, time, and magnitude of the disturbance events. A comparison of the CuSum algorithm with the LiDAR reference map resulted in a 78% success rate for the test site in Gabon and 65% success rate for the test site in Peru, for disturbances as small as 0.01 ha in size and for canopy height losses as fine as 10 m. A correlation between the change metric and above ground biomass (AGB) change was found with R2 = 0.95, and R2 = 0.83 for canopy height loss. From the regression model we directly estimated local AGB loss maps for the year 2020, at 1 ha scale and in percentages of AGB loss. Comparison with the Global Forest Watch (GFW) Tree Cover Loss (TCL) product showed a 61% overlap between the two maps when considering only deforested pixels, with 504 ha of deforestation detected by CuSum vs. 348 ha detected by GFW. Low intensity disturbances captured by the CuSum method were largely undetected by GFW and by the SAR-based Radar for Detecting Deforestation (RADD) Alert System. The results of this study confirm this approach as a simple and reproducible change detection method for monitoring and quantifying fine-scale to high intensity forest disturbances, even in the case of multi-storied and high biomass forests

Table_1_Reliably mapping low-intensity forest disturbance using satellite radar data.XLSX

In the last decades tropical forests have experienced increased fragmentation due to a global growing demand for agricultural and forest commodities. Satellite remote sensing offers a valuable tool for monitoring forest loss, thanks to the global coverage and the temporal consistency of the acquisitions. In tropical regions, C-band Synthetic Aperture Radar (SAR) data from the Sentinel-1 mission provides cloud-free and open imagery on a 6- or 12-day repeat cycle, offering the unique opportunity to monitor forest disturbances in a timely and continuous manner. Despite recent advances, mapping subtle forest losses, such as those due to small-scale and irregular selective logging, remains problematic. A Cumulative Sum (CuSum) approach has been recently proposed for forest monitoring applications, with preliminary studies showing promising results. Unfortunately, the lack of accurate in-situ measurements of tropical forest loss has prevented a full validation of this approach, especially in the case of low-intensity logging. In this study, we used high-quality field measurements from the tropical Forest Degradation Experiment (FODEX), combining unoccupied aerial vehicle (UAV) LiDAR, Terrestrial Laser Scanning (TLS), and field-inventoried data of forest structural change collected in two logging concessions in Gabon and Peru. The CuSum algorithm was applied to VV-polarized Sentinel-1 ground range detected (GRD) time series to monitor a range of canopy loss events, from individual tree extraction to forest clear cuts. We developed a single change metric using the maximum of the CuSum distribution, retrieving location, time, and magnitude of the disturbance events. A comparison of the CuSum algorithm with the LiDAR reference map resulted in a 78% success rate for the test site in Gabon and 65% success rate for the test site in Peru, for disturbances as small as 0.01 ha in size and for canopy height losses as fine as 10 m. A correlation between the change metric and above ground biomass (AGB) change was found with R2 = 0.95, and R2 = 0.83 for canopy height loss. From the regression model we directly estimated local AGB loss maps for the year 2020, at 1 ha scale and in percentages of AGB loss. Comparison with the Global Forest Watch (GFW) Tree Cover Loss (TCL) product showed a 61% overlap between the two maps when considering only deforested pixels, with 504 ha of deforestation detected by CuSum vs. 348 ha detected by GFW. Low intensity disturbances captured by the CuSum method were largely undetected by GFW and by the SAR-based Radar for Detecting Deforestation (RADD) Alert System. The results of this study confirm this approach as a simple and reproducible change detection method for monitoring and quantifying fine-scale to high intensity forest disturbances, even in the case of multi-storied and high biomass forests.</p

Image_3_Reliably mapping low-intensity forest disturbance using satellite radar data.JPEG

In the last decades tropical forests have experienced increased fragmentation due to a global growing demand for agricultural and forest commodities. Satellite remote sensing offers a valuable tool for monitoring forest loss, thanks to the global coverage and the temporal consistency of the acquisitions. In tropical regions, C-band Synthetic Aperture Radar (SAR) data from the Sentinel-1 mission provides cloud-free and open imagery on a 6- or 12-day repeat cycle, offering the unique opportunity to monitor forest disturbances in a timely and continuous manner. Despite recent advances, mapping subtle forest losses, such as those due to small-scale and irregular selective logging, remains problematic. A Cumulative Sum (CuSum) approach has been recently proposed for forest monitoring applications, with preliminary studies showing promising results. Unfortunately, the lack of accurate in-situ measurements of tropical forest loss has prevented a full validation of this approach, especially in the case of low-intensity logging. In this study, we used high-quality field measurements from the tropical Forest Degradation Experiment (FODEX), combining unoccupied aerial vehicle (UAV) LiDAR, Terrestrial Laser Scanning (TLS), and field-inventoried data of forest structural change collected in two logging concessions in Gabon and Peru. The CuSum algorithm was applied to VV-polarized Sentinel-1 ground range detected (GRD) time series to monitor a range of canopy loss events, from individual tree extraction to forest clear cuts. We developed a single change metric using the maximum of the CuSum distribution, retrieving location, time, and magnitude of the disturbance events. A comparison of the CuSum algorithm with the LiDAR reference map resulted in a 78% success rate for the test site in Gabon and 65% success rate for the test site in Peru, for disturbances as small as 0.01 ha in size and for canopy height losses as fine as 10 m. A correlation between the change metric and above ground biomass (AGB) change was found with R2 = 0.95, and R2 = 0.83 for canopy height loss. From the regression model we directly estimated local AGB loss maps for the year 2020, at 1 ha scale and in percentages of AGB loss. Comparison with the Global Forest Watch (GFW) Tree Cover Loss (TCL) product showed a 61% overlap between the two maps when considering only deforested pixels, with 504 ha of deforestation detected by CuSum vs. 348 ha detected by GFW. Low intensity disturbances captured by the CuSum method were largely undetected by GFW and by the SAR-based Radar for Detecting Deforestation (RADD) Alert System. The results of this study confirm this approach as a simple and reproducible change detection method for monitoring and quantifying fine-scale to high intensity forest disturbances, even in the case of multi-storied and high biomass forests.</p