Monitoring tropical forest succession at landscape scales despite uncertainty in Landsat time series

Long-term annual estimation of forest above ground biomass, canopy cover, and height from airborne and spaceborne sensors synergies in the Iberian Peninsula

The Mediterranean Basin has experienced substantial land use changes as traditional agriculture decreased and population migrated from rural to urban areas, which have resulted in a large forest cover increase. The combination of Landsat time series, providing spectral information, with lidar, offering three-dimensional insights, has emerged as a viable option for the large-scale cartography of forest structural attributes across large time spans. Here we develop and test a comprehensive framework to map forest above ground biomass, canopy cover and forest height in two regions spanning the most representative biomes in the peninsular Spain, Mediterranean (Madrid region) and temperate (Basque Country). As reference, we used lidar-based direct estimates of stand height and forest canopy cover. The reference biomass and volume were predicted from lidar metrics. Landsat time series predictors included annual temporal profiles of band reflectance and vegetation indices for the 1985-2023 period. Additional predictor variables including synthetic aperture radar, disturbance history, topography and forest type were also evaluated to optimize forest structural attributes retrieval. The estimates were independently validated at two temporal scales, i) the year of model calibration and ii) the year of the second lidar survey. The final models used as predictor variables only Landsat based metrics and topographic information, as the available SAR time-series were relatively short (1991-2011) and disturbance information did not decrease the estimation error. Model accuracies were higher in the Mediterranean forests when compared to the temperate forests (R2 = 0.6-0.8 vs. 0.4-0.5). Between the first (1985-1989) and the last (2020-2023) decades of the monitoring period the average forest cover increased from 21 ± 2% to 32 ± 1%, mean height increased from 6.6 ± 0.43 m to 7.9 ± 0.18 m and the mean biomass from 31.9 ± 3.6 t ha-1 to 50.4 ± 1 t ha-1 for the Mediterranean forests. In temperate forests, the average canopy cover increased from 55 ± 4% to 59 ± 3%, mean height increased from 15.8 ± 0.77 m to 17.3 ± 0.21m, while the growing stock volume increased from 137.8 ± 8.2 to 151.5 ± 3.8 m3 ha-1. Our results suggest that multispectral data can be successfully linked with lidar to provide continuous information on forest height, cover, and biomass trends.

A Novel Method for Estimating Spatial Distribution of Forest Above-Ground Biomass Based on Multispectral Fusion Data and Ensemble Learning Algorithm

Optical remote sensing technology has been widely used in forest resources inventory. Due to the influence of satellite orbits, sensor parameters, sensor errors, and atmospheric effects, there are great differences in vegetation spectral information captured by different satellite sensor images. Spectral fusion technology can couple the advantages of different multispectral sensor images to produce new multispectral data with high spatial and spectral resolution, it has great potential for improving the spectral sensitivity of forest vegetation and alleviating the spectral saturation. However, how to quickly and effectively select the multi-spectral fusion data suitable for forest above-ground biomass (AGB) estimation is a very critical issue. This study proposes a scheme (RF-S) to comprehensively evaluate multispectral fused images and develop the appropriate model for forest AGB estimation, on the basis of random forest (RF) and the stacking ensemble algorithm. First, four classic fusion methods are used to fuse the preprocessed GaoFen-2 (GF-2) multispectral image with Sentinel-2 image to generate 12 fused Sentinel-like images. Secondly, we apply a comprehensive evaluation method to quickly select the optimal fused image for the follow-up research. Subsequently, two feature combination optimization methods are used to select feature variables from the three feature sets. Finally, the stacking ensemble algorithm based on model dynamic integration and hyperparameter automatic optimization, as well as some classic machine learners, are used to construct the forest AGB estimation model. The results show that the fused image NND_B3 (based on nearest neighbor diffusion pan sharpening method and Band3_Red) selected by the evaluation method proposed in this study has the best performance in AGB estimation. Using the stacking ensemble method and NND_B3 image, we get the highest estimation accuracy, with the adjusted R2 and relative root mean square error (RMSEr) of 0.6306 and 15.53%, respectively. The AGB estimation RMSEr of NND_B3 is 19.95% and 24.90% lower than those of GF-2 and Sentinel-2, respectively. We also found that the multi-window texture factor has better performance in the area with low AGB, and it can suppress the overestimation significantly. The AGB spatial distribution estimated using the NND_B3 image matches the field observations well, indicating that the multispectral fusion image combined with the Stacking algorithm can increase the accuracy and saturation of the AGB estimates.