To accurately understand and investigate carbon fluxes in cropland ecosystems, this study adopted a machine learning ensemble model for estimation. Focusing on the Jinzhou station of the ChinaFLUX, we integrated eddy covariance carbon flux observations with multi-source satellite remote sensing data to construct a machine learning-based cropland carbon flux estimation model. For environmental driver selection, a strategy combining correlation analysis with ecological mechanism understanding was employed to screen LST, NDVI, and NDMI as model input variables, effectively avoiding multicollinearity issues. Using footprint-weighted integrated data from 2005 to 2014 for model training and validation, a Stacking ensemble model was constructed with the RF model serving as the meta-learner to stack the predictions of RF, CART, and GBM. The ensemble model further reduced the prediction error (RMSE = 39.82), maintaining an R2 > 0.9 in most years and effectively improving predictive performance during anomalous years where single models underperformed. Based on these findings, the model was applied to analyze the spatiotemporal evolution of NEE in Jinzhou croplands from 2005 to 2014. The analysis revealed that while the region functioned overall as a carbon sink, it exhibited significant spatiotemporal heterogeneity. Spatially, the distribution followed a pattern of “strong intensity in the northeast and center, and weak intensity in the northwest and southwest.” Temporally, the sink intensity underwent significant interannual oscillations characterized by a “strengthening–weakening–re-strengthening–declining” trajectory. The high-precision prediction method proposed in this study is of great significance for revealing spatiotemporal variations in carbon sources/sinks, guiding green agricultural development, and supporting relevant policy formulation.
Jiang et al. (2026) studied this question.