Accurately estimating forest aboveground biomass (AGB) is imperative for comprehending carbon cycling, calculating carbon budgets, and formulating sustainable forest management plans. Currently, random forest (RF) and other machine learning models are widely used to estimate forest AGB, as they can effectively handle nonlinear relationships. However, by constructing a global model using all the samples collected from a study area, these models fail to account for the spatial heterogeneity in the AGB and cannot correct the prediction biases, thereby constraining the estimation accuracy. To overcome these limitations, we proposed a novel approach termed geographical random forest and empirical Bayesian kriging (GRFEBK). This hybrid model combines the localized modeling capability of geographical random forest (GRF) with the bias correction strength of empirical Bayesian kriging (EBK). GRF adapts RF to account for the spatial heterogeneity of the AGB, while EBK utilizes the spatial autocorrelation of residuals to correct the prediction deviations. This study was conducted in Hainan Island, utilizing spectral bands, vegetation indices, tasseled cap components derived from Landsat-8 imagery, backscattering coefficients from ALOS-2 synthetic aperture radar, topographic features, and the forest canopy height as the explanatory variables. A total of 195 forest aboveground biomass (AGB) samples were collected for modeling and assessing the predictive accuracy. The results demonstrate that, among the tested models, including GRFEBK, RF, support vector machine (SVM), k-nearest neighbor (KNN), geographically weighted regression (GWR), GRF, and EBK, GRFEBK attains the highest ${\mathrm{R}}^2$ (0.78) and the lowest RMSE (36.04 Mg/ha) and RRMSE (22.87%), significantly outperforming the conventional models and using GRF or EBK alone. These results demonstrate that by accounting for local non-stationarity in AGB and correcting prediction biases, GRFEBK achieves significantly higher accuracy than conventional RF and other models. While the results are promising, the computational cost of GRFEBK and its performance under varying geographical conditions warrant further investigation at larger scales to assess its broader applicability. Nevertheless, GRFEBK provides an innovative and more reliable approach for accurate forest AGB estimation with great potential to support global forest resource monitoring.

The accurate monitoring of crop areas is essential for food security and agriculture, but accurately extracting multiple-crop distribution over large areas remains challenging. To solve the above issue, in this study, the Pixel-based One-dimensional convolutional neural network (PB-Conv1D) and Pixel-based Bi-directional Long Short-Term Memory (PB-BiLSTM) were proposed to identify multiple-crop cultivated areas using time-series NaE (a combination of NDVI and EVI) as input for generating a baseline classification. Two approaches, Snapshot and Stochastic weighted averaging (SWA), were used in the base-model to minimize the loss function and improve model accuracy. Using an ensemble algorithm consisting of five PB-Conv1D and seven PB-BiLSTM models, the temporal vegetation index information in the base-model was comprehensively exploited for multiple-crop classification and produced the Pixel-Based Conv1D and BiLSTM Ensemble model (PB-CB), and this was compared with the PB-Transformer model to validate the effectiveness of the proposed method. The multiple-crop cultivated area was extracted from 2005, 2010, 2015, and 2020 in North China by using the PB-Conv1D combine Snapshot (PB-CDST) and PB-CB models, which are a performance-optimized single model and an integrated model, respectively. The results showed that the mapping results of the multiple-crop cultivated area derived by PB-CDST (OA: 81.36%) and PB-BiLSTM combined with Snapshot (PB-BMST) (OA: 79.40%) showed exceptional accuracy compared to PB-Transformer combined with Snapshot and SWA (PB-TRSTSA) (OA: 77.91%). Meanwhile, the PB-CB (OA: 83.43%) had the most accuracy compared to the pixel-based single algorithm. The MODIS-derived PB-CB method accurately identified multiple-crop areas for wheat, corn, and rice, showing a strong correlation with statistical data, exceeding 0.7 at the municipal level and 0.6 at the county level.

Accurately mapping the forest canopy height is vital for conserving forest ecosystems. Employing the forest height measured by satellite light detection and ranging (LiDAR) systems as ground samples to establish forest canopy height extrapolation (FCHE) models presents promising opportunities for mapping large-scale wall-to-wall forest canopy height. However, despite the potential to provide more samples and alleviate the stripe effect by synergistically using the data from two existing LiDAR datasets, Global Ecosystem Dynamics Investigation (GEDI) and Ice, Cloud, and land Elevation Satellite-2 (ICESat-2), the fundamental differences in their operating principles create measurement biases, and thus, there are few studies combining them for research. Furthermore, previous studies have typically employed existing regression algorithms as FCHE models to predict forest canopy height, without customizing a model that achieves optimal performance based on the current samples. These shortcomings constrain the accuracy of predicting forest canopy height using satellite LiDAR data. To surmount these difficulties, we proposed a genetic programming (GP) guided method for mapping forest canopy height by combining the GEDI and ICESat-2 LiDAR data with Sentinel-1/2, terrain, and climate data. In this method, GP autonomously constructs the fusion model of the GEDI and ICESat-2 datasets (hereafter GIF model) and the optimal FCHE model based on the explanatory variables for the specific study area. The outcomes demonstrate that the fusion of the GEDI and ICESat-2 data shows high consistency (${\mathrm{R}}^2$ = 0.85, RMSE = $2.2$$\mathrm{m}$, pRMSE = 11.24%). The synergistic use of the GEDI and ICESat-2 data, coupled with the optimization of the FCHE model, substantially improves the precision of forest canopy height predictions, and finally achieves ${\mathrm{R}}^2$, RMSE, and pRMSE of 0.64, $3.38$$\mathrm{m}$, and 16.08%, respectively. In summary, our research presents a reliable approach to accurately estimate forest canopy height using remote sensing data by addressing measurement biases between the GEDI and ICESat-2 data and overcoming the limitations of traditional FCHE models.

Spatiotemporal vegetation changes serve as a key indicator of regional ecological environmental quality and provide crucial guidance for developing strategies for regional ecological protection and sustainable development. Currently, vegetation change studies in the Yangtze River Basin primarily rely on the Normalized Difference Vegetation Index (NDVI). However, the NDVI is susceptible to atmospheric and soil conditions and exhibits saturation phenomena in areas with high vegetation coverage. In contrast, the kernel NDVI (kNDVI) demonstrates significant advantages in suppressing background noise and improving saturation thresholds through nonlinear kernel transformation, thereby enhancing sensitivity to vegetation changes. To elucidate the spatiotemporal characteristics and driving mechanisms of vegetation changes in the Yangtze River Basin, this study constructed a temporal kNDVI using MOD09GA data from 2000 to 2022. Considering sectional heterogeneity, rather than analyzing the entire region as a whole as in previous studies, this research examined spatiotemporal evolution characteristics by sections using four statistical metrics. Subsequently, Partial Least Squares Path Modeling (PLSPM) was innovatively introduced to quantitatively analyze the influence mechanisms of topographic, climatic, pedological, and socioeconomic factors. Compared to traditional correlation analysis and the geographical detector method, PLSPM, as a theoretically driven statistical method, can simultaneously process path relationships among multiple latent variables, effectively revealing the intensity and pathways of driving factors’ influences, while providing more credible and interpretable explanations for kNDVI variation mechanisms. Results indicate that the overall kNDVI in the Yangtze River Basin exhibited an upward trend, with the midstream demonstrating the most significant improvement with minimal interannual fluctuations, the upstream displaying an east-increasing and west-stable spatial pattern, and the downstream demonstrating coexisting improvement and degradation characteristics, with these trends expected to persist. Driving mechanism analysis reveals that the upstream was predominantly influenced by the climatic factor, the midstream was dominated by terrain, and the downstream displayed terrain–soil coupling effects. Based on these findings, it is recommended that the upstream focus on enhancing vegetation adaptation management to climate change, the midstream need to coordinate the relationship between topography and human activities, and the downstream should concentrate on controlling the negative impacts of urban expansion on vegetation.