1. Introduction
Land use and land cover (LULC) change has become a popular topic in recent decades because change in the landscape is linked to so many vital environmental processes [
1,
2]. Ecosystem function, biological cycles, and biodiversity can be profoundly changed by alteration of land cover [
3]. The most common and widespread LULC changes across the globe are deforestation, agricultural expansion, and urbanization [
4,
5,
6]. The trend in LULC has been toward cash crop plantations due to various government policies, growing physical infrastructure, social, and economic development in various parts of mountainous areas in Southeast Asia. Forests have become the main focus for conversion to agriculture plantation [
7].
The mountainous forested areas of Cambodia have been transformed into agricultural land through economic land concessions (ELC). ELC is the process of long-term (usually 70 to 90 years) granting the state-owned land to private sector operations for economic development through agricultural and industrial–agricultural operations, including large-scale plantations, stock rearing, and factory construction [
8]. According to the Ministry of Agriculture, Forestry, and Fishery (MAFF) of Cambodia, a total land area of 1,204,750 hectares of ELC was granted to 118 companies [
9] over the last decade. Land use in mountain areas has seen a rapid transformation from native forests to cash crops, including rubber, acacia, cassava, sugarcane, jatropha, and other crops [
10]. These changes have undoubtedly altered the distribution of biotic populations and pathways of ecosystem function [
11] in steepland terrain.
It is essential to have a good understanding of the historic trends and drivers of LULC change because these can help to establish practical and effective land conservation and management, as well as facilitate modeling of future LULC change [
2]. Remote sensing (RS) technology provides an effective way to extract a time series of LULC information [
12]. A variety of land use models have been developed to simulate and predict the future of LULC change [
13], mainly by coupling RS and geographic information system (GIS) data. Land use change models are tools used to analyze and provide a better understanding of the causes and consequences of land use changes. They can be used to investigate complex driving forces that affect the spatiotemporal patterns of LULC changes [
14]. They have been widely applied as a basis for land use planning, policy-making, and environmental impact assessment [
15]. Land use models have been used by researchers to analyze and forecast LULC change, employing techniques such as Markov chain [
16,
17], artificial neural network [
18], cellular automata [
19], and logistic regression [
20]. However, each model has its own limitation, and these have been discussed extensively in the literature [
21]. The ANN and Markov chain limitations were addressed in this study. For instance, MC is a qualitative model that does not incorporate various variables of land use change [
22], and ANN mode is the “black-box” that limits understanding of land use change [
21].
One strategy for dealing with these shortcomings has been to develop an integrated model that leverages the strengths of the various techniques [
23]. The land-use change modeler (LCM) tool is an integrated LULC change model that employs both multi-Layer perceptron neural network (MLPNN) and Markov chain (MC) methods. The MLPNN consists of three layers, named ‘input’, ‘hidden’, and ‘output’; these are used to define complex non-linear relationships between LULC changes and driving variables of change. This technique has been shown to produce accurate estimates of transition potential [
24,
25]. According to Roodposhti et al. [
26], “transition potential represents the behavioral propensities of the actors determining land-use change and is defined based on the inferred logic from a set of transition rules which offer an easily understandable and transparent way to find the most reliable land-use class allocation”. The MC method is applied to calculate the transition probability matrix of predicted LULC change, based on the transition results of MLPNN. This LCM has been successfully applied in many studies for analysis and prediction of urban growth around the world such as in Guanzhong, China [
11]; Patna district, India [
23]; Cairo Region, Egypt [
27]; Banda Aceh, Indonesia [
28]; Kailash Sacred Landscape, Nepal [
29]; Dhaka city, Bangladesh [
30]; and Asmara, Eritrea [
31]. However, only a few studies have carried out LULC change analysis and forecasting at the watershed scale. Two examples are the Neka Watershed in Iran [
32] and the Des Plaines River watershed in the USA [
33].
In this study, the LCM tool was applied to simulate and project LULC change in the upper Siem Reap basin, Cambodia, a catchment that has experienced considerable change over the last few decades. There has been no previous modelling work of LULC change in the catchment, and the LCM has not previously been applied in this monsoon-dominated mountain environment. This environment is characterized by two seasons: A dry season from November to April, and a rainy season from May to October [
34]. It is important to note that the fundamental assumption in our modelling framework is that state change is represented strictly as a function of probability—there is no consideration of contingencies that may be introduced by social, economic, or administrative factors.
The objective of this study is to evaluate LULC change trend and annual rate of deforestation from 1988 to 2018 using supervised classification, and to predict future LULC maps for years 2058 and 2098 using the LCM tool. The results of this paper will be used in future studies to further assess the impacts of LULC and future climate changes on river flows that maintain the Angkor temple complex, a World Heritage Site. The Angkor temple is sited on sand substrate [
35] and requires a constant water supply to avoid land subsidence and collapse. As outlined in the following section, the water supply to the temple complex and the municipality is vital for a number of economic, social, and cultural reasons, and thus it is of utmost importance to understand the factors that might impact this resource in the future.
2. Study Area
Siem Reap province, northwestern Cambodia, is the site of numerous ancient temples, including Angkor Wat temple. More than one million tourists visit Angkor Wat annually, and the sustainability of this vital cultural and spiritual landmark depends critically on the supply of water from the upper source catchment. Thus, our study focuses on the upper Siem Reap River watershed, located about 20 km from the Siem Reap downtown area and about 374 km from Phnom Penh capital city as in
Figure 1. It is located mainly within the Banteay Srei and Svay Leu districts, in the Siem Reap Province.
The hydrological regime of the catchment is driven by two seasons, a rainy season and a dry season. According to the Siem Reap Provincial Department of Water Resources and Meteorology (PDWRAM), the annual average temperature is 28 °C, with an average minimum of 25 °C in the rainy season, an average maximum of 34 °C in the dry season, and the annual average rainfall is 1580 mm in the rainy season in Siem Reap province.
The upper Siem Reap watershed includes large swaths of lowland evergreen natural forests and agricultural land. Much of the catchment has a fairly gentle slope, and is thus attractive for agricultural land use, especially in areas above the Angkor temple complex up to Kulen mountainous range and the eastern part of the upper Kulen Mountain. The landscape has changed significantly over the last few decades, most prominently by deforestation. There are several factors driving this change, including a loosely enforced land management framework, as well as socio-economic, demographic, and environmental factors (see
Figure 2). The Cambodian government framework on land management therefore plays a critical role in regulating LULC change in the catchment.
According to established boundaries, reflected in Ministry of Environment maps, about 54% of the study watershed should be a protected area/national park with forest cover, but recent fieldwork (early 2019), reveals significant portions of the protected areas have been disturbed, owing to poor enforcement of regulations. In periods running up to political elections, as politicians look for support from the local community, they may overlook illegal encroachments upon forest land. Consequently, these areas become legally titled to the occupiers following the election.
Although there are a few tourism sites in this upstream watershed, the majority of people who live there are farmers and some younger workers who seek employment in Siem Reap town. Typical family size is increasing, and thus, they require more land for farming and crop cultivation: Deforestation is a common option. According to the local authority of Anlong Thom village [
36] in the watershed area (early 2019), in some cases, there are migrants from other provinces who come into the watershed and clear parts of the forest for plantation and then resell the land.
Given the extent of these recent changes, the trend of land use and land cover is likely to affect the runoff characteristics of the basin. According to Bruijnzeel [
37] and Foley et al. [
1] LULC change can have significant influence on runoff and groundwater in the watershed. The conversion from natural forest to agriculture crop in this catchment might change the streamflow patterns and increase peak flows, owing to reduced interception by forest canopy. It can also influence weather factors, through changing surface temperature, increasing evaporation, and decreasing transpiration rates. Moreover, studies have found that there is a clear positive correlation between forest cover and water yield [
38]. Wang et al. [
39], for instance, found that mean annual runoff declined by 2.3% following an increase of 25% in the proportion of forestland in a mountainous watershed. Khoi and Suetsugi [
40] discovered an increase in streamflow, from 0.2 to 0.4%, was associated with a decline of 16.3% in forestland catchment. Since LULC has changed considerably over the last few decades, it is reasonable to assume that this will exert some influence on water yield and streamflow peaks in the upstream catchment of Siem Reap River, and thus may have significant effects on the Angkor temple complexes and the people in the watershed area.
4. Results
The results of the analysis are broken into two main parts, for clarity: (1) The LULC trend detection and annual rate of deforestation analysis, including the reasons behind the change trends and accuracy assessment and (2) land change modeler results, which consist of land change analysis from 1998 to 2018, simulation of LULC map year 2018 using MLPNN and MC methods, validation by comparing the simulated and the actual classification 2018, and prediction of LULC maps 2058 and 2098 using LULC maps 2008–2018 and calibration driver parameters that were adapted from the best results in the 2018 simulation.
4.1. LULC Trend Detection and Annual Rate of Deforestation Analysis
Based on the information from the classified images for the last three decades in the upstream catchment of the Siem Reap River, there is a trend of decreasing forest cover with a corresponding, inverse trend in non-forest cover. It can be seen in
Figure 5 that forest cover increased very slightly in the first decade (1988–1998). Following this, there was a decreasing trend in the last two decades (1998–2008 and 2008–2018); the reasons behind those changes are discussed at the end of this section.
According to this analysis, the annual deforestation rate for the last 30 years was 1.22%, which is similar to rates of forest cover changes in the Greater Mekong Sub-region, particularly Cambodia and Myanmar at (−1.20%) from 1990 to 2015 [
48]. While the net gain in forest cover was 0.57% annually for the first decade, over the last two decades (1998–2008 and 2008–2018) there was a net decrease of 2.32% and 1.91%, respectively. It emerges that deforestation proceeded at a remarkable rate −36.61% between 1988 and 2018. Thus, non-forest areas became the majority cover type in the watershed. Forest covered areas, inversely, decreased by almost half, from 81.08% (331.87 km
2) to 44.47% (182.02 km
2) over the same epoch (see
Table 2). The extent of forested land retreated to more upland or mountainous areas, or sites with increasingly marginal potential for cultivation and other development. Pixels with water cover type (natural lakes, wide river reaches) made up less than 0.08% of the catchment for 17 years (1988, 1990, 1992, 1993, 1999, 2000, 2004, 2005, 2007, and from 2009 to 2018), and there was no water detected in the remaining 12 years.
There were two main factors behind the low rate of forest cover change in the first decade: (1) Limited access to the road network, and (2) landmine hazard, since this area was a stronghold of the Khmer Rouge. In addition to this, cultivators followed traditional practices and remained within their farming boundaries, with little additional land clearance. There was less illegal logging and agricultural expansion, only small-scale activities, and there was firewood and charcoal production for household usage only. Therefore, the disturbance areas were primarily in the lower elevation areas close to the road and village environs.
In the second decade, the proportion of forest cover declined more quickly: More than one-third of the catchment area was stripped. The most likely cause of this relatively rapid deforestation can be traced to improved access to road infrastructure and landmine clearance. These improvements may have facilitated illegal logging activities, firewood and charcoal production, and cultivation. Gaughan et al. [
49] underlined that land use changes from 1989–2005 in the Siem Reap watershed could be related to these factors, which were triggered by explosive growth in tourism.
In the third decade, the extent of forested area gradually dropped almost 10% with some ups and downs due to fluctuations in the costs of agricultural production. Farms often change their agriculture crops, and some years they do not grow any crops; hence, natural regeneration could be improved. A notable peak in 2008 might relate to the global economic crisis, which heavily affected agriculture, tourism, and other sectors. The main reason behind the decrease in forest coverage in this last decade could be further improvement of road infrastructure and large-scale agricultural expansion, including shifting cultivation that has been a significant forest clearance in the Mount Kulen range.
The key drivers were determined based on informal interviews/meetings with regional and local officials (see
Section 3.2.3, above), and these are summarized in
Table 3.
4.2. Accuracy Assessment
While the average overall accuracy of the classification of all images (1988–2018) was 91.46%, the kappa statistic was 82.48%. The minimum overall accuracy and kappa were 86% and 70.76%, respectively, and the maximum values were 97.37% and 94.92%. There was no individual class with an overall accuracy below 70%, which was deemed satisfactory. Annual overall accuracy and kappa statistic for each frame over the last three decades are shown in
Table 4.
4.3. The Land Change Modeler (LCM)
4.3.1. LULC Change Analysis
The LULC maps from 1998, 2008, and 2018 were selected to perform the change analysis and cross-tabulate analysis, thus the surface areas (gains and losses) and spatial changes of each class were assessed for Epoch 1 (1998–2008), Epoch 2 (2008–2018), and Epoch 3 (1998–2018). In this analysis, there were four possible transitions: Non-Forest remained the same (NF-NF), non-forest changed to forest or regeneration (NF-F), forest stayed the same (F-F), and forest changed to non-forest (F-NF). During Epoch 1, while forest covered areas had decreased 104.89 km
2 (25.62%) and become non-forest areas (F-NF), areas that regenerated from non-forest to forest (NF-F) constituted 10 km
2 (2.44%), thus the net deforestation was 94.89 km
2 (23.1%). During Epoch 2, there appeared to be a decrease of 93.20 km
2 (22.77%) of forest cover (F-NF) and the forest recovery (NF-F) was 14.84 km
2 (3.63%), hence the net deforestation was 78.36 km
2 (19.14%). During Epoch 3, the conversion from forest cover to non-forest (F-NF) was 179.31 km
2 (43.81%) and the regrowth forest (NF-F) was 6.06 km
2 (1.48%). The regeneration of forest from one epoch to another was not cumulative because it was disturbed or wiped out in the later epoch. Therefore, the net deforestation was 173.25 km
2 (42.08%). The pattern of gains and losses for the LULC classes in all epochs is shown in
Figure 6a–c.
It can be seen from the land change maps that there were substantial losses of forest cover and reciprocal gains in non-forest cover in all periods. The LULC change and gains and losses maps of each class (forest and non-forest) from 1998–2018 are shown in
Figure 7a–c.
4.3.2. The Spatial Trend of Change
A generalized trend surface, representing the spatial weighting of the state transitions within a given epoch, is a useful technique for visualizing the geographic ‘hot-spots’ of change. This tool was developed by Eastman [
25]. The level of detail can be generalized by using a lower polynomial order, while higher-order (up to 9 terms) surfaces reveal more spatial details of the trend. In this study, a ninth-order surface was selected for a detailed visualisation of the trend of land use change across the whole catchment. The spatial trend maps of both forest and non-forest changes in Epoch 1 (1998–2008), 2 (2008–2018), and 3 (1998–2018) are shown in
Figure 8a–f.
In
Figure 8a,c, it can be seen that the spatial trend of change in forest cover in Epochs 1 and 3 occurred in similar locations. A notable change happened around the west of the map (dense red colour), along the main road (red and yellow colours), in the east of the catchment area (yellow colour), and in the lower right (south-east) of the map area (the Mount Kulen range). These changes were closely linked to the landmine clearance and development of access road infrastructure. Greater road access contributed to increased firewood and charcoal production, illegal logging activities, shifting cultivation, and agricultural expansion.
During Epoch 2, significant forest change took place in a few locations, such as along the main road (most notably in the north), in the north-west and the south-east of map (
Figure 8b). Trends in forest change continued the pattern established in Epoch 1, with further incursion into the former forest. There are three driving factors that may be responsible for this: Further development of access roads such as the main road (under construction on Mount Kulen), secondary roads, village roads, and trekking roads (around the watershed areas), large-scale agricultural expansion, and shifting cultivation (mainly cashew plantation).
Deforestation began at lower elevations and proceeded to higher elevations. While in the first epoch (1998–2008) deforestation occurred in areas with elevations from 40 to 100 m above mean sea level (a.m.s.l.), in the second epoch (2008–2018), further encroachment of forest took place at elevations of 100 to around 120 a.m.s.l. On the other hand, there was a different pattern of deforestation along the top of Mount Kulen. The disturbances are scattered throughout the mountain range, since farmers have been shifting cultivation for more than six decades. However, the density of these scattered farming plots has increased considerably over the last two decades (both Epochs 1 and 2: 1998–2008 and 2008–2018) with the elevation from 320 to 420 a.m.s.l.
The trend surfaces within non-forested areas in the epoch of 1998–2008 (
Figure 8d) and 1998–2018 (
Figure 8f) exhibit a similar form, with hot-spots of change in the southwest part of the map and southeast part of the map or in the Mount Kulen range. However, a majority of change happened in the southwest part of the map (close to the lower main road), the north part (close to upper main road), the northwest and the southeast of the map (
Figure 8e) during 2008–2018. Forest regeneration generally occurred in the mountainous areas such as the north part (close to upper main road), the northwest and the southeast of the map, and in areas close to the protected area of Apsara Authority or in temple areas in the southwest part of the map (close to the lower main road). There are often encroachments upon this government-owned land due to lax enforcement of land laws.
4.3.3. Multi-Layer Perceptron Markov chain (MLP-MC) analysis
MLPNN-MC methods were applied to simulate LULC change to year 2018 using LULC maps years 1998 and 2008 (Epoch 1) and drivers of change. A further aim was to predict LULC maps out to years 2058 and 2098 using LULC maps years 2008 to 2018 (Epoch 2) and successful calibration of drivers setting of simulated LULC map 2018.
The MLP-MC method, as other spatial models, normally carried data, model, and future predictions uncertainty [
50]. Noise in land use data input can cause uncertainty, but it was reduced by a validation process by comparing classified images with ground truth [
51] and field survey. Further reduction of uncertainty of input data was carried out by MLPNN-MC, which is capable of handling data uncertainties [
51,
52,
53]. During training and testing, the model worked with drivers and hidden layers under 10,000 interactions to generate the expected transition of Epoch 1 by minimizing the root mean square (RMS) (training RMS 0.165 and testing RMS 0.164). Moreover, sensitivity analysis can lower the model uncertainty [
53,
54]. Sensitivity analysis of drivers for LULC changes between 1998 to 2008 can be tested through MLPANN. Each driver was fed in the model separately to identify their influential factors contributing to land use change in the study catchment. The relative influence of each driver was evaluated by the model as shown in
Table 5.
The most influential factors were elevation, distance from main road, and distance from villages, and the least influence drivers were distance from minor road, slope, and aspect. Elevation was one of the most effective discriminant variables for LULC change because deforestation happens mainly at lower elevations across the catchment. Similarly, distance from main road was a good predictor, since roads provided easy access for logging. Another significant driver factor is distance from villages. It can clearly be seen that villages in the catchment are loci for forest clearance, particularly those villages located along the main road. As the population grows, rates of deforestation surrounding villages will continue to grow.
A combination of all seven drivers of change and LULC maps in Epoch 1 were fed in the MLPNN, and it produced the potential transition maps (F-NF and NF-F) with an accuracy rate of 93.29%, which was satisfactory and easily surpassed the acceptance threshold of 80% [
45].
The MC technique was then used with the potential transition results from MLPNN to generate the transition probability matrix and simulate LULC map 2018. As shown in
Table 6, the probability of non-forest areas remaining as non-forest areas (NF-NF) was 81.51%, and the probability of change in the non-forest land to forest land (NF-F) was only 18.49% in year 2018. The probability of forested areas staying the same (F-F) was 70.48% and the probability of change from forest areas to non-forest areas (F-NF) was only 29.52% in the same year.
After simulation of LULC map 2018 using MLPANN-MC, a validation process was carried out by comparing pixels of the simulated and actual LULC maps 2018 (
Figure 9a,b) using kappa index statistics. Based on this comparison, the statistical results show that the general kappa statistic or Cohen’s kappa statistic, K
standard was 0.8112, K
no value was 0.8349, K
location was 0.8560, and K
locationStrata was 0.8560. According to this statistic, all kappa index values were greater than 0.70, showing good agreement between simulated and actual LULC change. It indicated that the simulated and observed land use were very similar.
Furthermore, it is difficult to overcome the uncertainty of future LULC maps since it is inherently spatial and cannot be represented in model [
51,
55]. Therefore, it is assumed that LULC change in Epoch 2 (2008–2018) would follow the same trend in the next four and eight decades.
4.3.4. Predicting Future LULC Maps 2058 and 2098
After successful validation of modelled and actual LULC maps 2018, simulation of LULC maps 2058 and 2098 was performed using the same setting of driver parameters and LULC maps in Epoch 2. The transition probability matrixes of simulated LULC maps 2058 and 2098 are shown in
Table 7 and
Table 8.
Based on
Table 7, the probability of non-forest areas remaining in the same class (NF-NF) was 81.11% and the probability of change in the non-forest land to forest land (NF-F) was 19.89% in the year 2058. The probability of forested areas staying the same (F-F) was 28.55% and the probability of change from forest areas to non-forest areas (F-NF) was only 71.45% in the same year. According to
Table 8, the probability of non-forest areas remaining the same class (NF-NF) was 78.39% and the probability of change in the non-forest land to forest land (NF-F) was 21.61% in the year 2098. The probability of forested areas staying the same (F-F) was 22.36% and the probability of change from forest areas to non-forest areas (F-NF) was only 77.64% in the same year.
The LULC maps 2058 and 2098 images were generated as shown in
Figure 10. The prediction revealed that there would be an increase in non-forest areas, from 227.30 km
2 (55.53%) to 313.12 km
2 (76.49%) between 2018 and 2058. Conversely, the total forest area would decrease from 182.02 km
2 (44.47%) to 96.20 km
2 (23.50%) in the same period. Furthermore, the non-forest covered areas would also continue to increase from 313.12 km
2 (76.49%) to 320.44 km
2 (78.28%) between 2058 and 2098. The forest covered areas would decrease from 96.20 km
2 (23.50%) to 88.88 km
2 (21.71%) in the same period. There was a small proportion of water in both predicted years, accounting for less than 0.01%. While the non-forest covered class could be expanded up to 22.75% of the catchment from 2018 to 2098, the forest cover naturally could be dropped the same proportion in upstream of Siem Reap River, Cambodia. LULC maps 2058 and 2098 are presented in
Figure 10a,b.
The results from this model seem to indicate that there would be a continuous upward trend of non-forest and downward trend of forest cover in the future at the upstream of Siem Reap River. However, this model did not include a number of factors that might affect the system trajectory, including government policies and framework, socio-economic, and demographic information.
5. Conclusions and Discussion
LULC in the Siem Reap watershed has changed considerably over the last few decades. Given these highly dynamic conditions, and the great importance of water supply to the Angkor Wat complex, determining the current and future LULC is crucial. Based on the analysis, there has been a strongly decreasing trend in forest cover over the last three decades (1988–2018). The deforestation over that interval has been remarkable, at 36.61%. The average annual deforestation rate is 1.22%. Non-Forest areas have thus become the majority cover type in the watershed. The main reasons behind these trends are landmine clearance and road infrastructure development, which have contributed significantly to the expansion of agriculture, shifting cultivation, illegal logging, and firewood and charcoal production.
According to the LCM analysis, the majority of forest loss occurred within a low band of terrain elevation, between 40 to 120 m a.m.s.l., as well as in areas along the main road from south to north, in the west part of the catchment, and in some parts of Kulen National park at elevations of 320 to 420 a.m.s.l. The remaining forest cover is mainly located on the steep slopes of mountainous areas. This is consistent with the relative weighting of driver influence (
Table 5), showing that elevation, distance from major roads, and distance from villages are more important than distance from minor road and river, slope, and aspect. The model also reveals that there would be a continuous upward trend of non-forest and downward trend of forest cover in the future (year 2058 and 2098) at the upstream of Siem Reap River, Cambodia. This forecasting results might have high possibility according to the field visit on the early of 2019 because it is observed that a main road and secondary roads are under construction along Mount Kulen, which provide an easy access and promote further deforestation in Mount Kulen range, but an actual spatial and quantity change might be different from the simulation LULC change (can be less or more than the prediction results). Overall, this model works satisfactorily and acceptable to simulate and forecast the LULC change in this study.
This model can be improved by including other missing parameters of change such as government policies and framework on land management, socio-economic, and demographic information, which have significant impact on the land use changes and future change. The LCM may be an imperfect model, but it can help to simulate the possible future of LULC maps with sufficient and visible driver factors of change. In addition, although GIS and LCM enable researchers to analyse and predict future LULC maps scenarios, high-resolution images should be used in the model in order to achieve further accuracy and prediction results.
Furthermore, since the indigenous communities have lived within the Mount Kulen National Park boundary, the scattered plots of regrowth and deforestation have appeared over the period of time. This evidence shows the shifting cultivation which is commonly practiced by people live in the upland regions. However, those scattered areas of deforestation have been converted into permanent agriculture since 2009, and those scatter plots have been increased significantly in the last decade. Therefore, the Cambodian government as well as related agencies should have visible plans either to identify the protected areas and shifting cultivation boundaries or even relocate the indigenous communities outside the protected boundary. Moreover, reforestation and forest-regrow should be taken into account by Cambodian government to restore forest cover and sustain the river flow since this watershed is important in providing water to maintain the Angkor complex, particularly Angkor Temple, which is a World Heritage Site.
The results of this study will contribute to the sustainable land management and development managers and related agencies involved in the decision-making process. These outcomes will also be an important step for further study on the effect of LULC and climate changes on hydrology because the continuous decrease of forest cover in the watershed may influence the hydrological components such as infiltration rate, runoff, and base flow [
56], as well as water yield, which could threaten the unique characteristics of the river, including sustainable of Angkor temple, drinking supplies, and irrigation.