ABSTRACT The objective of this study is to develop a methodological approach to assess ecosystem functioning in terms of carbon uptake and sequestration. To achieve this, a model was developed to forecast Net Primary Productivity (NPP) and tested in a rural area, the Thrace region (NE Greece). Monthly net photosynthesis data from July 2015 to February 2024 was estimated using the MODerate Resolution Imaging Spectroradiometer (MODIS) Aqua MYD17A2H v061 dataset for the study area. This dataset was used as the dependent variable in an Autoregressive Integrated Moving Average (ARIMA) model. The predictive ability of monthly Land Cover (LC) and Soil Moisture (SM) data at various time lags was examined, and variables found to be significant predictors were introduced into the model as external predictors. The study found that crop areas and SM conditions exert the greatest influence on NPP, with areas containing trees exerting a lesser impact. Additionally, NPP was estimated for seven different LC and SM scenarios from March 2023 to February 2024.