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.

The Bekasi River Basin is highly vulnerable to severe and recurrent flooding, as evidenced by significant infrastructure and environmental damage during major events. This study investigates the catastrophic floods of 2016, 2020, 2022, and 2025 by implementing the Rainfall-Runoff-Inundation (RRI) model to simulate key hydrological processes. After validation using historical water level data, the model performed effectively, achieving the highest coefficient of determination (R² = 0.75) and lowest root mean square error (RMSE = 0.66) at Cileungsi Station. In contrast, the lowest R² = 0.02, and the highest RMSE = 3.74 at Pondok Gede Permai (PGP) Station. The results reveal a concerning trend of worsening 5-year flood events, with the 2025 flood reaching a peak inundation depth exceeding 3 m and affecting an area of 2.97 km², caused by a rainfall threshold of more than 180 mm/day. Furthermore, the model shows a rapid hydrological response, with a time lag of approximately 7 h or less between peak rainfall and flood onset across three monitoring stations. Analysis indicates these severe floods were primarily triggered by heavy rainfall combined with significant land cover changes. The findings provide valuable insights for flood prediction and mitigation strategies in this vulnerable region.

Water conservation (WC) is a vital ecosystem service (ES) that plays an essential role in the sustainable management of water resources and ensures ecological security. This research examines the WC capacity of the Citarum watershed in West Java Province, Indonesia, from 2010 to 2020. The specific objectives of this research are as follows: (1) to assess the 10-year temporal and spatial variations of WC using the Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) and topographic index model; (2) to analyze the temporal interchange between WC and its influencing factors through the Spatial Elastic Coefficient Trajectory Model (SECTM); and (3) to determine the driving factors (DFs) affecting WC by applying Multiscale Geographically Weighted Regression (MGWR). Key findings of this study reveal significant trends in WC from 2010 to 2020: the average WC in the Citarum watershed decreased from 513.96 mm/ha to 453.64 mm/ha (11.74%), indicating a concerning decline in ES capacity. This study also noted that regions implementing comprehensive regulations dominated the study area, covering approximately 72.70% of the total area (502,295 hectares). This illustrates that the implementation of rules plays a crucial role in the WC. Lastly, the MGWR analysis found that WC in the Citarum watershed positively correlated with topography, climate, and vegetation while negatively correlating with socioeconomic factors. This indicates that WC levels are generally lower in areas with higher human activity and economic growth, highlighting the impact of anthropogenic pressures on natural resources. This framework helps stakeholders plan to ensure sustainable development in the area, as it provides valuable insights into the interactions between the ecological and socioeconomic factors affecting WC.

Population growth is increasing food demand, while watershed food systems face intense pressure from land-use/land-cover (LULC) change. Therefore, this study aimed to identify the critical threshold for food sufficiency in the Cimanuk Watershed. The methods used include dynamic analysis and prediction of LULC changes in line with four scenarios using Support Vector Machine (SVM)–Markov model, ecosystem service-based food production modeling through Crop Production module in Integrated Valuation of Ecosystem Services and Trade-offs (InVEST) model, food demand estimation based on population projections, and calculation of Food Sufficiency Index (FSI) to identify the critical threshold at which the food system shifts from surplus to deficit. The results showed that paddy rice production is projected to continue meeting food demand across all scenarios through 2042. However, food demand is expected to increase as the population grows, specifically in the Accelerated Population Growth (APGS) scenario. FSI analysis indicated that the Cimanuk Watershed remained in food surplus (FSI > 1) up to 2042. The watershed may transition to a food-deficit condition after 2062 in the absence of intervention measures. Therefore, this study recommends protecting productive paddy fields, conserving forests, and managing population growth to maintain long-term food sufficiency.