Assessing Proxy-Based Grassland Gross Primary Productivity Using Machine Learning Approaches and Multi-Source Remote Sensing

Gross Primary Productivity (GPP) in grassland ecosystems is a fundamental eco-biophysical indicator for assessing carbon cycling, grazing capacity, and ecosystem responses to climatic stress. However, robust estimation of GPP in arid and semi-arid rangelands remains challenging because of pronounced spatial heterogeneity, strong climate variability, and inherent uncertainties associated with remotely sensed observations. Together, these factors constrain both modeling performance and out-of-sample generalization beyond the training domain. In this dryland grassland context, this study compares the performance of machine learning (ML) models for grassland GPP proxy-based characterization, downscaling, and predictive agreement using a multivariate dataset that integrates Sentinel-2-derived spectral and phenological features, a Moderate-Resolution Imaging Spectroradiometer (MODIS)-derived GPP proxy, and complementary climatic and geographic information. Pixel-level observations spanning multiple years are analyzed, with ordinary linear regression used as a baseline benchmark and ensemble decision-tree models, including Random Forest, Gradient Boosting, and Histogram-based Gradient Boosting (HGB), compared. Instead of relying solely on random cross-validation, model performance is systematically assessed using a combination of spatially structured validation and a leave-one-year-out scheme to explicitly examine spatial and temporal generalization. The results indicate that ensemble tree-based models outperform linear approaches, with the HGB model showing the strongest agreement with the MODIS-derived GPP proxy (R 2 = 0.95, RMSE = 0.035 on the test set) and maintaining stable performance across spatial and temporal validations (R 2 = 0.86–0.96 across years). Taken together, the findings demonstrate that integrating multi-source remote sensing data with climatic information within a rigorous validation framework enables a more reliable assessment of model generalization and gap-filling consistency with respect to a remote-sensing-based proxy target, rather than an absolute validation against ground-based measurements, thereby supporting sustainability-relevant monitoring of arid grassland ecosystems.