Inversion of Grassland Aboveground Biomass in the Three Parallel Rivers Area Based on Genetic Programming Optimization Features and Machine Learning

Aboveground biomass (AGB) in grasslands is a vital metric for assessing ecosystem functioning and health. Accurate and efficient AGB estimation is essential for the scientific management and sustainable use of grassland resources. However, achieving low-cost, high-efficiency AGB estimation via remote sensing remains a key challenge. This study integrates Sentinel-1 and Sentinel-2 imagery to derive 38 multi-source feature variables, including backscatter coefficients, texture, spectral reflectance, vegetation indices, and topographic factors. These features are combined with AGB data from 112 field plots in the Three Parallel Rivers area. Feature selection was performed using Pearson correlation, Random Forest (RF), and SHAP values to identify optimal variable sets. Genetic Programming (GP) was then applied for nonlinear optimization of the selected features. Three machine learning models—RF, GBRT, and KNN—were used to estimate AGB and generate spatial distribution maps. The results revealed notable differences in model accuracy, with RF performing best overall, outperforming GBRT and KNN. After GP optimization, all models showed improved performance, with the RF model based on RF-selected features achieving the highest accuracy (R2 = 0.90, RMSE = 0.31 t/ha, MAE = 0.23 t/ha), improving R2 by 0.03 and reducing RMSE and MAE by 0.05 and 0.03 t/ha, respectively. Spatial mapping showed the AGB ranged from 0.41 to 3.59 t/ha, with a mean of 1.39 t/ha, closely aligned with the actual distribution characteristics. This study demonstrates that the RF model, combined with multi-source features and GP optimization, provides an effective approach to grassland AGB estimation and supports ecological monitoring in complex areas.