Tropical regions such as Java, Indonesia, still lack publication of soil water retention (SWR) information, particularly at upper Citarum watershed. The SWR is one of the critical elements in water storage and movement in the soil and very important to solve ecological and environmental problems. However, getting the access requires a lot of laboratory meas-urement that is time-consuming and expensive. Therefore, utilizing pedotransfer functions (PTFs) to estimate the water in the soil is needed. This study aims to define soil properties related to the SWR and to evaluate the performance of existing PTFs in predicting SWR. The study was carried out at agroforestry land system soil at upper Citarum watershed, Indonesia. Ten point and two continuous existing PTFs developed for tropical regions were applied in this study. Pearson's correlation (r), mean error (ME), root mean square error (RMSE), and modelling efficiency (EF) were used for evaluation. Cation exchange capacity (CEC), organic carbon (OC), bulk density (BD), and clay were considered as potential soil properties for soil water retention prediction. The performance of PTFs by MINASNY, HARTEMINK [2011] at matric potential of –10 kPaand BOTULA [2013] at matric potential of –33 kPa and –1500 kPa were recommended for point PTFs, while PTFs by HODNETT, TOMASELLA [2002] was for continuous PTFs in predicting SWR. The accuracy of the point PTFs is almost better than the continuous PTFs in predicting SWR in agroforestry land system soil at upper Citarum watershed, Indonesia.

Exploring the drivers of changes in ecosystem services is crucial to maintain ecosystem functionality, especially in the diverse Central Citarum watershed. This study utilises the integrated valuation of ecosystem service and trade-offs (InVEST) model and multiscale geographically weighted regression (MGWR) model to examine ecosystem services patterns from 2006 to 2018. The InVEST is a hydrological model to calculate water availability and evaluate benefits provided by nature through simulating alterations in the amount of water yields driven by land use/cover changes. Economic, topographic, climate, and vegetation factors are considered, with an emphasis on their essential components. The presence of a geographical link between dependent and explanatory variables was investigated using a multiscale geographic weighted regression model. The MGWR model is employed to analyse spatial impacts. The integration of both models simplified the process and enhanced its understanding. The findings reveal the following patterns: 1) decreasing land cover and increasing ecosystem services demand in the watershed, along with a decline in water yield, e.g. certain sub-districts encounter water scarcity, while others have abundant water resources; 2) the impact of natural factors on water yield shifts along vegetation > climate > topography (2006) changes to climate > vegetation > topography (2018).