Forest stock volume (FSV) is a basic data source for estimating forest carbon sink. It is also a crucial parameter that reflects the quality of forest resources and forest management level. The use of remote sensing data combined with a support vector regression (SVR) algorithm has been widely used in FSV estimation. However, due to the complexity and spatial heterogeneity of the forest biological community, in the FSV high-value area with dense vegetation, the optical re-mote sensing variables tend to be saturated, and the sensitivity of synthetic aperture radar (SAR) backscattering features to the FSV is significantly reduced. These factors seriously affect the ac-curacy of the FSV estimation. In this study, Landsat 8 (L8) Operational Land Imager multispectral images and C-band Sentinel-1 (S1) hyper-temporal SAR images were used to extract three re-mote sensing feature datasets: spectral variables (L8), backscattering coefficients (S1), and inter-ferometric SAR factors (S1-InSAR). We proposed a feature selection method based on SVR (FS-SVR) and compared the FSV estimation performance of FS-SVR and stepwise regression analysis (SRA) on the aforementioned three remote sensing feature datasets. Finally, an estima-tion model of coniferous FSV was constructed using the SVR algorithm in Wangyedian Forest Farm, Inner Mongolia, China, and the spatial distribution map of coniferous FSV was predicted. The experimental results show the following: (1) The coherence amplitude and DSM data ob-tained based on S1 images contain information relat-ed to forest canopy height, and the hy-per-temporal S1 image data significantly enrich the diversity of S1-InSAR feature factors. There-fore, the S1-InSAR dataset has a better FSV response than remote sensing factors such as the S1 backscattering coefficient and L8 vegetation index, and the corresponding root mean square er-ror (RMSE) and relative RMSE (rRMSE) values reached 47.6 m³/ha and 20.9%, respectively. (2) The integrated dataset can provide full play to the synergy of the L8, S1, and S1-InSAR remote sensing data. Its RMSE and rRMSE values are 44.3 m³/ha and 19.4% respectively. (3) The proposed FS-SVR method can better select remote sensing variables suitable for FSV estimation than SRA. The average value of the rRMSE (23.17%) based on the three datasets was 13.8% lower than that of the SRA method (26.87%). This study provides new insights into forest FSV retrieval based on active and passive multisource remote sensing joint data.

Unaffected by cloud cover and solar illumination, synthetic aperture radar (SAR) images coupled with quad-polarimetric techniques have significant potential for mapping forest aboveground biomass (AGB) in the mountains of southern China. To improve the accuracy of mapping forest AGB, it is necessary to accurately interpret and evaluate the sensitivity of polarimetric features related to polarimetric response in complex forests. In this study, several rotated polarimetric features were extracted from L-band quad-polarimetric ALOS PALSAR-2 images based on uniform polarimetric matrix rotation theory. In addition, the sensitivity of rotated polarimetric features with forest parameters was evaluated by the Pearson correlation coefficient, sensitivity index (SI), and saturation levels. Ultimately, the forest AGB was mapped with various combinatorial feature sets by a proposed feature selection method based on the sensitivity index. The results illustrated that rotated polarimetric features extracted from the rotational domain have higher sensitivity with various forest parameters and higher saturation levels for mapping forests than other traditional features. After using the proposed feature selection method and combinatorial feature sets, the rRMSE of mapped forest AGB ranged from 22.5% to 33.9% for two acquired images, and the best result was obtained from the combination of three types of polarimetric features (BC + C4 + Ro). It is also confirmed that different types of features extracted from quad-polarimetric SAR images have better compensation effects and the accuracy of mapped forest AGB is significantly improved.

As one of the important types of forest resources, mapping forest stock volume (FSV) in larch (Larix decidua) forests holds significant importance for forest resource management, carbon cycle research, and climate change monitoring. However, the accuracy of FSV mapping using common spectral and texture features is often limited due to their failure in fully capturing seasonal changes and growth cycle characteristics of vegetation. Phenological features can effectively provide essential information regarding the growth status of forests. In this study, multi-temporal Sentinel-2 satellite imagery were initially acquired in the Wangyedian Forest Farm in Chifeng City, Inner Mongolia. Subsequently, various phenological features were extracted from time series variables constructed by Gaussian Process Regression (GPR) using Savitzky–Golay filters, stepwise differentiation, and Fourier transform techniques. The alternative features were further refined through Pearson’s correlation coefficient analysis and the forward selection algorithm, resulting in six groups of optimal subsets. Finally, four models including the Random Forest (RF), K-Nearest Neighbors (KNN), Support Vector Machine (SVM), and Multiple Linear Regression (MLR) algorithms were developed to estimate FSV. The results demonstrated that incorporating phenological features significantly enhanced model performance, with the SVM model exhibiting the best performance—achieving an R² value of 0.77 along with an RMSE value of 46.36 m³/hm² and rRMSE value of 22.78%. Compared to models without phenological features, inclusion of these features led to a 0.25 increase in R² value while reducing RMSE by 10.40 m³/hm² and rRMSE by 5%. Overall, integration of phenological feature variables not only improves the accuracy of larch forest FSV mapping but also has potential implications for delaying saturation phenomena.

Unaffected by cloud cover and solar illumination, synthetic aperture radar (SAR) images have great capability to map forest growing stem volume (GSV) in complex biophysical environments. Up to now, c-band dual-polarization Gaofen-3 (GF-3) SAR images, acquired by the first Chinese civilian satellite equipped with multi-polarized modes, are rarely applied in mapping forest GSV. To evaluate the capability of dual-polarization GF-3 SAR images in mapping forest GSV, several proposed derived features were initially extracted by mathematical operations and applied to obtain optimal feature sets by different feature sorting methods and feature selection methods. Then, the maps of GSV in an evergreen coniferous forest were inverted by various machine learning algorithms and stacking ensemble learning methods with different strategies. The results implied that backscattering coefficients and partially proposed derived features showed high sensitivity to the forest GSV, and the saturation phenomenon also obviously occurred once the forest GSV was larger than 300 m³/ha. Furthermore, the results showed that the accuracy of the mapped GSV was significantly improved using the stacking ensemble learning methods. Using various optimal feature sets and base models (MLR, KNN, SVM, and RF), the rRMSE values mainly ranged from 30% to 40%. After using the stacking ensemble learning methods, the values of rRMSE ranged from 16.71% to 20.51%. This confirmed that dual-polarization GF-3 images have great potential to map forest GSV in evergreen coniferous forests.

Polarimetric Synthetic Aperture Radar (PolSAR) images with dual polarization modes have great potential to map forest stock volume (FSV) by excellent penetration capabilities and distinct microwave scattering processes. However, the response of these SAR data to FSV is still uncertain in the non-growing season. To further interpret the response of FSV to different dual polarization SAR images, three types of dual polarization SAR images (GF-3, Sentinel-1, and ALOS-2) were initially acquired in coniferous planted forest in the non-growing season. Then, sensitivity between FSV and all alternative features extracted from each type of SAR image was analyzed to express the response of FSV to dual polarization SAR images with bands and polarization modes in the non-growing season in deciduous (Larch) and evergreen (Chinese pine) forests. Finally, mapped FSV using single and combined dual polarization images were derived by optimal feature sets and four machine learning models, respectively. The combined effects were also analyzed to clarify the difference of bands and polarization modes in deciduous and evergreen forests in the non-growing season. The results demonstrated that the backscattering energy from different sensors is significantly different in Chinese pine, and the difference is gradually reduced in Larch forests. It is also implied that the polarization mode is more important than penetration capability in mapping forest FSV in deciduous forest in the non-growing season. By comparing the accuracy of mapped FSV using single and combined images, combined images have more capability to improve the accuracy and reliability of mapped FSV. Meanwhile, it is confirmed that compensation effects with bands and polarization modes not only have great potential to delay the saturation phenomenon, but also have the capability to reduce errors caused by overestimation.

SAR images with two polarizations show strong potential for mapping forest stock volume (FSV) combined with limited samples. However, accurately mapping FSV still presents challenges in selecting the optimal acquisition date to obtain the SAR images during specific phenological stages within the annual forest cycle (growth and dormant stages). To clarify the impacts of phenological stages within the annual cycle on FSV mapping, SAR images with various polarization models and bands (Sentinel-1(S), GaoFen-3(GF-3 (G)) and ALOS-2(A)) were acquired within the growth and dormant stages of an annual cycle in a boreal evergreen coniferous forest (Chinese pine) and a deciduous coniferous forest (Larch). Subsequently, single-band (G, S, and A) and multi-band combined alternative variable sets (A + G, A + S, S + G, and A + S + G) were extracted within the same stage, respectively. Finally, the forward selection approach was utilized in conjunction with four different models (MLR, KNN, RF, and SVR) to obtain the most suitable variable sets and generate FSV mapping. The results demonstrated a strong correlation between the intensity of backscattering coefficients and the phenological stages of the forest. Within the dormant stage, there was a significant decrease in the gaps of backscattering coefficients obtained from the same polarization compared to those within the growth stage. Furthermore, the results also revealed that more signals from inside the canopy could be detected during the dormant stage in both evergreen coniferous forests and deciduous coniferous forests. Subsequently, the accuracy in mapping FSV obtained from single-band SAR images within the dormant stage are slightly higher than that within the growth stage, and the accuracy was still significantly affected by both overestimation and underestimation. Moreover, the combined effects of different bands significantly improve the reliability of mapped FSV. The rRMSE values in four multi-band combinations ranged from 22.37% to 29.40% for Chinese pine forests and from 21.27% to 34.38% for Larch forests, and the optimal result was observed from combinations of A + S + G acquired within the dormant stage. It is confirmed that SAR signal and their sensitivity to FSV depends on the stages of forest annual growth cycle. In comparison to the growth period, dual-polarization SAR data acquired during the dormant stage is more suitable for estimating FSV in boreal forests.