Mapping the Forest Canopy Height in Northern China by Synergizing ICESat-2 with Sentinel-2 Using a Stacking Algorithm

Mapping global mangrove canopy height by integrating Ice, Cloud, and Land Elevation Satellite-2 photon-counting LiDAR data with multi-source images

Large-scale and precise measurement of mangrove canopy height is crucial for understanding and evaluating wetland ecosystems' condition, health, and productivity. This study generates a global mangrove canopy height map with a 30 m resolution by integrating Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) photon-counting light detection and ranging (LiDAR) data with multi-source imagery. Initially, high-quality mangrove canopy height samples were extracted using meticulous processing and filtering of ICESat-2 data. Subsequently, mangrove canopy height models were established using the random forest (RF) algorithm, incorporating ICESat-2 canopy height samples, Sentinel-2 data, TanDEM-X DEM data and WorldClim data. Furthermore, a global 30 m mangrove canopy height map was generated utilizing the Google Earth Engine platform. Finally, the global map's accuracy was evaluated by comparing it with reference canopy heights derived from both space-borne and airborne LiDAR data. Results indicate that the global 30 m resolution mangrove height map was found to be consistent with canopy heights obtained from space-borne (r = 0.88, Bisa = -0.07 m, RMSE = 3.66 m, RMSE% = 29.86 %) and airborne LiDAR (r = 0.52, Bisa = -1.08 m, RMSE = 3.39 m, RMSE% = 39.05 %). Additionally, our findings reveal that mangroves worldwide exhibit an average height of 12.65 m, with the tallest mangrove reaching a height of 44.94 m. These results demonstrate the feasibility and effectiveness of using ICESat-2 data integrated with multi-source imagery to generate a global mangrove canopy height map. This dataset offers reliable information that can significantly support government and organizational efforts to protect and conserve mangrove ecosystems.

Continuous mapping of forest canopy height using ICESat-2 data and a weighted kernel integration of multi-temporal multi-source remote sensing data aided by Google Earth Engine

Forest Canopy Height (FCH) is a crucial parameter that offers valuable insights into forest structure. Spaceborne LiDAR missions provide accurate FCH measurements, but a significant challenge is their point-based measurements lacking spatial continuity. This study integrated ICESat-2's ATL08-derived FCH values with multi-temporal and multi-source remote sensing (RS) datasets to generate continuous FCH maps for northern forests in Iran. Sentinel-1/2, ALOS-2 PALSAR-2, and FABDEM datasets were prepared in Google Earth Engine (GEE) for FCH mapping, each possessing unique spatial and geometrical characteristics that differ from those of the ATL08 product. Given the importance of accurately representing the geometrical characteristics of the ATL08 segments in modeling FCH, a novel Weighted Kernel (WK) approach was proposed in this paper. The WK approach could better represent the RS datasets within the ATL08 ground segments compared to other commonly used resampling approaches. The correlation between all RS data features improved by approximately 6% compared to previously employed approaches, indicating that the RS data features derived after convolving the WK approach are more predictive of FCH values. Furthermore, the WK approach demonstrated superior performance among machine learning models, with random forests outperforming other models, achieving a coefficient of determination (R2) of 0.71, root mean square error (RMSE) of 4.92 m, and mean absolute percentage error (MAPE) of 29.95%. Furthermore, in contrast to previous studies using only summer datasets, this study included spring and autumn data from Sentinel-1/2, resulting in a 6% increase in R2 and a 0.5-m decrease in RMSE. The proposed methodology filled the research gaps and improved the accuracy of FCH estimations.

Integrating spaceborne LiDAR and Sentinel-2 images to estimate forest aboveground biomass in Northern China

BACKGROUND

Fast and accurate forest aboveground biomass (AGB) estimation and mapping is the basic work of forest management and ecosystem dynamic investigation, which is of great significance to evaluate forest quality, resource assessment, and carbon cycle and management. The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2), as one of the latest launched spaceborne light detection and ranging (LiDAR) sensors, can penetrate the forest canopy and has the potential to obtain accurate forest vertical structure parameters on a large scale. However, the along-track segments of canopy height provided by ICESat-2 cannot be used to obtain comprehensive AGB spatial distribution. To make up for the deficiency of spaceborne LiDAR, the Sentinel-2 images provided by google earth engine (GEE) were used as the medium to integrate with ICESat-2 for continuous AGB mapping in our study. Ensemble learning can summarize the advantages of estimation models and achieve better estimation results. A stacking algorithm consisting of four non-parametric base models which are the backpropagation (BP) neural network, k-nearest neighbor (kNN), support vector machine (SVM), and random forest (RF) was proposed for AGB modeling and estimating in Saihanba forest farm, northern China.

RESULTS

The results show that stacking achieved the best AGB estimation accuracy among the models, with an R² of 0.71 and a root mean square error (RMSE) of 45.67 Mg/ha. The stacking resulted in the lowest estimation error with the decreases of RMSE by 22.6%, 27.7%, 23.4%, and 19.0% compared with those from the BP, kNN, SVM, and RF, respectively.

CONCLUSION

Compared with using Sentinel-2 alone, the estimation errors of all models have been significantly reduced after adding the LiDAR variables of ICESat-2 in AGB estimation. The research demonstrated that ICESat-2 has the potential to improve the accuracy of AGB estimation and provides a reference for dynamic forest resources management and monitoring.

Performance of GEDI Space-Borne LiDAR for Quantifying Structural Variation in the Temperate Forests of South-Eastern Australia

Monitoring forest structural properties is critical for a range of applications because structure is key to understanding and quantifying forest biophysical functioning, including stand dynamics, evapotranspiration, habitat, and recovery from disturbances. Monitoring of forest structural properties at desirable frequencies and cost globally is enabled by space-borne LiDAR missions such as the global ecosystem dynamics investigation (GEDI) mission. This study assessed the accuracy of GEDI estimates for canopy height, total plant area index (PAI), and vertical profile of plant area volume density (PAVD) and elevation over a gradient of canopy height and terrain slope, compared to estimates derived from airborne laser scanning (ALS) across two forest age-classes in the Central Highlands region of south-eastern Australia. ALS was used as a reference dataset for validation of GEDI (Version 2) dataset. Canopy height and total PAI analyses were carried out at the landscape level to understand the influence of beam-type, height of the canopy, and terrain slope. An assessment of GEDI’s terrain elevation accuracy was also carried out at the landscape level. The PAVD profile evaluation was carried out using footprints grouped into two forest age-classes, based on the areas of mountain ash (Eucalyptus regnans) forest burnt in the Central Highlands during the 1939 and 2009 wildfires. The results indicate that although GEDI is found to significantly under-estimate the total PAI and slightly over-estimate the canopy height, the GEDI estimates of canopy height and the vertical PAVD profile (above 25 m) show a good level of accuracy. Both beam-types had comparable accuracies, with increasing slope having a slightly detrimental effect on accuracy. The elevation accuracy of GEDI found the RMSE to be 10.58 m and bias to be 1.28 m, with an R2 of 1.00. The results showed GEDI is suitable for canopy densities and height in complex forests of south-eastern Australia.

Forest Canopy Height Mapping by Synergizing ICESat-2, Sentinel-1, Sentinel-2 and Topographic Information Based on Machine Learning Methods

Spaceborne LiDAR has been widely used to obtain forest canopy heights over large areas, but it is still a challenge to obtain spatio-continuous forest canopy heights with this technology. In order to make up for this deficiency and take advantage of the complementary for multi-source remote sensing data in forest canopy height mapping, a new method to estimate forest canopy height was proposed by synergizing the spaceborne LiDAR (ICESat-2) data, Synthetic Aperture Radar (SAR) data, multi-spectral images, and topographic data considering forest types. In this study, National Geographical Condition Monitoring (NGCM) data was used to extract the distributions of coniferous forest (CF), broadleaf forest (BF), and mixed forest (MF) in Hua’ nan forest area in Heilongjiang Province, China. Accordingly, the forest canopy height estimation models for whole forest (all forests together without distinguishing types, WF), CF, BF, and MF were established, respectively, by Radom Forest (RF) and Gradient Boosting Decision Tree (GBDT). The accuracy for established models and the forest canopy height obtained based on estimation models were validated consequently. The results showed that the forest canopy height estimation models considering forest types had better performance than the model grouping all types of forest together. Compared with GBDT, RF with optimal variables had better performance in forest canopy height estimation with Pearson’s correlation coefficient (R) and the root-mean-squared error (RMSE) values for CF, BF, and MF of 0.72, 0.59, 0.62, and 3.15, 3.37, 3.26 m, respectively. It has been validated that a synergy of ICESat-2 with other remote sensing data can make a crucial contribution to spatio-continuous forest canopy height mapping, especially for areas covered by different types of forest.

Improving Leaf Area Index Retrieval Using Multi-Sensor Images and Stacking Learning in Subtropical Forests of China

The leaf area index (LAI) is a key indicator of the status of forest ecosystems that is important for understanding global carbon and water cycles as well as terrestrial surface energy balances and the impacts of climate change. Machine learning (ML) methods offer promising ways of generating spatially explicit LAI data covering large regions based on optical images. However, there have been few efforts to analyze the LAI in heterogeneous subtropical forests with complex terrain by fusing high-resolution multi-sensor data from the Sentinel-1 Synthetic Aperture Radar (SAR), Sentinel-2 Multi Spectral Instrument (MSI), and Advanced Land Observing Satellite-1 digital elevation model (DEM). Here, forest LAI mapping was performed by integrating the MSI, SAR, and DEM data using a stacking learning (SL) approach that incorporates distinct predictions from a set of optimized individual ML algorithms. The method’s performance was evaluated by comparison to field forest LAI measurements acquired in Xingguo and Gandong of subtropical China. The results showed that the addition of the SAR and DEM images using the SL model compared to the inputs of only optical images reduced the mean absolute error (MAE) and root mean square error (RMSE) by 26% and 18%, respectively, in Xingguo, and by 12% and 8%, respectively, in Gandong. Furthermore, the combination of all images had the best prediction performance. SL was found to be more robust and accurate than conventional individual ML models, while the MAE and RMSE were decreased by 71% and 64%, respectively, in Xingguo, and by 68% and 59%, respectively, in Gandong. Therefore, the SL model using the three-source data combination produced satisfied prediction accuracy with the coefficients of determination (R2), MAE, and RMSE of 0.96, 0.17, and 0.28, respectively, in Xingguo and 0.94, 0.30, and 0.47, respectively, in Gandong. This study revealed the potential of the SL algorithm for retrieving the forest LAI using multi-sensor data in areas with complex terrain.

The Effect of Synergistic Approaches of Features and Ensemble Learning Algorith on Aboveground Biomass Estimation of Natural Secondary Forests Based on ALS and Landsat 8

Although the combination of Airborne Laser Scanning (ALS) data and optical imagery and machine learning algorithms were proved to improve the estimation of aboveground biomass (AGB), the synergistic approaches of different data and ensemble learning algorithms have not been fully investigated, especially for natural secondary forests (NSFs) with complex structures. This study aimed to explore the effects of the two factors on AGB estimation of NSFs based on ALS data and Landsat 8 imagery. The synergistic method of extracting novel features (i.e., COLI1 and COLI2) using optimal Landsat 8 features and the best-performing ALS feature (i.e., elevation mean) yielded higher accuracy of AGB estimation than either optical-only or ALS-only features. However, both of them failed to improve the accuracy compared to the simple combination of the untransformed features that generated them. The convolutional neural networks (CNN) model was much superior to other classic machine learning algorithms no matter of features. The stacked generalization (SG) algorithms, a kind of ensemble learning algorithms, greatly improved the accuracies compared to the corresponding base model, and the SG with the CNN meta-model performed best. This study provides technical support for a wall-to-wall AGB mapping of NSFs of northeastern China using efficient features and algorithms.