Spatio-Temporal Dynamics of Land Use and Land Cover Change and Ecosystem Service Value Assessment in Citarum Watershed, Indonesia: A Multi-Scenario and Multi-Scale Approach

Rapid land use and land cover (LULC) changes in densely populated watersheds pose serious challenges to the sustainability of ecosystem services (ES), yet their spatially explicit economic consequences remain insufficiently understood. This study analyzes the spatio-temporal dynamics of LULC and ecosystem service values (ESVs) in the Citarum Watershed, Indonesia, one of the country’s most critical and intensively transformed watersheds. Multi-temporal Landsat imagery from 2003, 2013, and 2023 was classified using a Random Forest algorithm, while future LULC conditions for 2043 were projected using a Multi-layer Perceptron–Markov Chain (MLP–MC) model under three scenarios: Business-as-Usual (BAU), Protecting Paddy Field (PPF), and Protecting Forest Area (PFA). ESVs were quantified at multiple spatial scales (county, 250 m grids, and 100 m grids) using both the Traditional Benefit Transfer (TBT) method and a Spatial Benefit Transfer (SBT) approach that integrates biophysical indicators with socio-economic variables. The contribution of LULC transitions to ESV dynamics was further assessed using the Ecosystem Service Change Intensity (ESCI) index. The results reveal substantial historical forest and shrubland losses, alongside rapid expansion of settlements and dryland agriculture, indicating intensifying anthropogenic pressure on watershed functions. Scenario analysis shows continued degradation under BAU, limited mitigation under PPF, and improved forest retention under PFA; although settlement expansion persists across all scenarios. Total ESV declined from USD 2641.33 million in 2003 to USD 1585.01 million in 2023, representing a cumulative loss of 46.13%. Projections indicate severe ESV losses under BAU and PPF by 2043, while PFA substantially reduces, but does not eliminate economic degradation. ESCI results identify forest and shrubland conversion to settlements and dryland agriculture as the dominant drivers of ESV decline. These findings demonstrate that integrating multi-scenario LULC modeling with spatially explicit ESV assessment provides a more robust basis for ecosystem-based spatial planning and supports sustainable watershed management under increasing development pressure.