Land cover change and its response to water level around Tonle Sap Lake in 1988-2020

地理学报(英文版) ›› 2024, Vol. 34 ›› Issue (2) : 329-354.

PDF(5826 KB)
PDF(5826 KB)
地理学报(英文版) ›› 2024, Vol. 34 ›› Issue (2) : 329-354. DOI: 10.1007/s11442-024-2207-1
研究论文

作者信息 +

Land cover change and its response to water level around Tonle Sap Lake in 1988-2020

Author information +
文章历史 +

Abstract

The transboundary influence of environmental change is a critical issue in the Lancang-Mekong region. As the largest river-connected lake in the lower Mekong, the ecological change and influence of Tonle Sap Lake have received widespread attention and discussion, especially after 2008, when the hydrological regime of the Lancang-Mekong River mainstream underwent distinct changes. However, the linkage and coupling mechanism between the lake riparian environment and mainstream water level change are still unclear. In this study, the interannual spatiotemporal changes in land cover in the Tonle Sap Lake riparian zone (TSLRZ) and their relationship with mainstream water levels were analysed. The results showed that the expansion of farmland was the most notable change in 1988-2020. After 2008, the land cover changes intensified, manifested as accelerated farmland expansion, intensified woodland fragmentation and significant water body shrinkage. Furthermore, the responses of the water body, degraded land, wasteland and grassland areas to the mainstream water levels weakened after 2008. Evidently, the land cover changes in the TSLRZ in the last 30 years were less related to the mainstream water level change than to local reclamation and logging. These results can offer a new scientific basis for the transboundary influence analysis of hydrological change.

Key words

riparian zone land cover change / mainstream water level change and transboundary influence / Tonle Sap Lake / Lancang-Mekong River

引用本文

导出引用
ZHANG Jing, MA Kai, FAN Hui, HE Daming. [J]. Journal of Geographical Sciences, 2024, 34(2): 329-354 https://doi.org/10.1007/s11442-024-2207-1
ZHANG Jing, MA Kai, FAN Hui, HE Daming. Land cover change and its response to water level around Tonle Sap Lake in 1988-2020[J]. Journal of Geographical Sciences, 2024, 34(2): 329-354 https://doi.org/10.1007/s11442-024-2207-1

1 Introduction

The transboundary influence of environmental change is a critical issue in the context of global change, especially in the Lancang-Mekong River Basin, which flows through six countries, including China, Myanmar, Laos, Thailand, Cambodia and Vietnam (He et al., 2016; Wang et al., 2021; Zhong et al., 2021). Tonle Sap Lake (abbreviated as TSL) in Cambodia, connected to the Lancang-Mekong River mainstream by the Tonle Sap River, is the largest freshwater lake in Southeast Asia and the largest fishery base in Asia, and it is an important habitat for endangered and rare species, with very important ecological and economic value (Lamberts, 2008; Ziv et al., 2012).
In recent decades, with increasingly obvious global climate change and intensified human activities, the ecological environment of the TSL area has deteriorated, showing effects such as lake shrinkage, rare species habitat destruction and biodiversity loss and receiving widespread attention and discussion (Arias, 2013; Lin and Qi, 2017; Pokhrel et al., 2018; Mahood et al., 2020; Morovati et al., 2023). In this area, the ecological changes are directly affected by local climate change and human activities and are also related to the hydrological regime (such as water level and runoff) of the mainstream since the majority of the supply of lake water comes from the mainstream on an annual scale (the supply proportion may fluctuate over time) (Kummu et al., 2014; Ji et al., 2018; Chua et al., 2022; Morovati et al., 2023). Therefore, the attribution of the ecological changes in the TSL area is extremely complex and sensitive.
Land cover changes affect resource availability, environmental conditions, habitat structure and spatial configuration, closely linking human socio-economic activities with natural ecological processes (Lambin et al., 2001; Turner et al., 2007; Sterling et al., 2013; He et al., 2022; Wang et al., 2022). The Tonle Sap Lake riparian zone (abbreviated as TSLRZ) is flat and broad and has various land cover types (JICA, 1999; Eng and Ouch, 2006; Mahood et al., 2020). Compared with the instantaneous change in the lake water body (indicated by the permanent water body and its associated water body areas, abbreviated as LWB), the land cover changes in the TSLRZ indicate a more stable and comprehensive manifestation of various impacts, including local human activities and mainstream water levels, which can more accurately represent the status and quality of the eco-environment. Thus, the analysis of the relationship between the land cover change in the TSLRZ and the mainstream water level change can further reveal their hydroecological linkage and coupling mechanism.
To date, there have been numerous studies on the hydrological regime changes of the Lancang-Mekong River (abbreviated as LMR), especially after the reservoir impoundment of the Xiaowan mega dam in 2008, which was the first mega dam in the mainstream of Lancang River (upper Mekong) (Li et al., 2017; Binh et al., 2020; Yun et al., 2020; Guan and Zheng, 2021; Lu and Chua, 2021; Chua et al., 2022; Dang et al., 2022; Morovati et al., 2023). Most of these studies showed that the hydrological regime of the mainstream experienced distinct changes after 2008, indicating that the water level or flow increased in the dry season and decreased in the wet season, although whether the dominant influencing factor is climate change or mega dam operation remains controversial. Moreover, studies on the changes in the various land cover types in the TSL area have made some progress (Zhan et al., 2002; Trung et al., 2013; Kim et al., 2019; Zhao et al., 2019; Huang et al., 2020). In particular, the studies by Meng et al. (2018), Mahood et al. (2020) and Chen et al. (2022) focused on these changes in the area around the TSLRZ (abbreviated as A-TSLRZ) based on Landsat images or European Space Agency Climate Change Initiative (ESA CCI) land cover products, which have temporal scales of 26 years, 5 years and 9 years, respectively. These previous studies showed that the area of farmland increased and other land cover types generally shrank in the A-TSLRZ from 1990 to 2016 (Meng et al., 2018), from 1993 to 2018 (Mahood et al., 2020) and from 1992 to 2019 (Chen et al., 2022), and they further discussed the influence of local precipitation, reclamation and logging on land cover changes. In addition, since 2008, scholars have paid more attention to the relationship between the land cover change in the TSL area and the hydrological regime of the mainstream and the lake (Arias, 2013; Trung et al., 2013; Ji et al., 2018; Wang et al., 2020; Gu et al., 2021). For instance, Trung et al. (2013) analysed the seasonal change in these land cover types (forest, shrub, grass and agriculture) in the northwestern part of the TSL during the period from 2007 to 2010 based on Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) data and Advanced Land Observation Satellite (ALOS) data and its response to the lake water level. Arias (2013) proposed a spatial modelling framework to link the habitat, flood duration and animal species in the A-TSLRZ. Ji et al. (2018) explored the relationship between the LWB and mainstream runoff during the period from 2000 to 2014 at an annual scale based on MODIS images. Gu et al. (2021) discussed the relationship between the LWB and the mainstream water level during the period from 1987 to 2017 at monthly and annual scales based on Landsat images.
In these studies, the overall characteristics of the land cover changes in the TSLRZ and the relationship between the change in the LWB and the hydrological regime of the mainstream were clarified, providing a great foundation for our study. However, the time scales for the change analyses of the other land cover types excluding water bodies in these studies were relatively broad due to different research objectives, while the water levels of the TSL and the LMR mainstream varied and fluctuated greatly every year. On different time scales, the relationship between the two is difficult to clarify and still unclear. In addition, examining broad time scales makes it difficult to understand the detailed spatiotemporal changes in land cover in the TSLRZ, especially after 2008, when the hydrological regime of the mainstream underwent distinct changes.
In this study, we combined land cover classification maps extracted from Landsat images and hydrologic measurements to study the changes in land cover and its major driving factors in the TSLRZ, with two main objectives: (1) to create an annual scale land cover classification dataset and document interannual spatiotemporal changes in land cover in the TSLRZ during the period from 1988 to 2020, focusing on the changes after 2008; and (2) to clarify the relationship between these land cover changes and the water level of the LMR mainstream.

2 Materials and methods

2.1 Study area

Tonle Sap Lake is located on the alluvial plain of the lower Lancang-Mekong River (Figure 1). The lake displays a distinct dichotomy between dry and wet seasons due to the influence of the tropical monsoon climate and the changes in the recharge of the LMR (Tangdamrongsub et al., 2016; Siev et al., 2018; Morovati et al., 2023). During the wet season, the water level of the LMR is higher than that of the TSL, reversing the flow in the Tonle Sap River and discharging into the lake, which results in a maximum lake area greater than 15,000 km2. During the dry season, the water level of the LMR becomes lower, and the lake water pours into the river when the water level of the river falls below that of the lake, which results in a minimum lake area of less than 2500 km2. The annual mean temperature of the lake area is 28-33℃, and the annual precipitation is between 1300 and 1500 mm (Siev et al., 2018; Daly et al., 2020).
Figure 1 Location and spatial pattern of the study area. The false-colour maps of the Tonle Sap Lake riparian zone (TSLRZ) are composed of the 5, 4 and 3 bands of Landsat images, which were collected in 1988 and 2020. The extent of Tonle Sap Lake (permanent water body) was extracted from Landsat images during 1988-2020.

Full size|PPT slide

The TSLRZ (Figure 1) includes the area from the lakeside to the lake road, with a total area of approximately 15,149.72 km2 (Figure 1). There are three Biosphere Reserve core areas belonging to the United Nations Educational, Scientific and Cultural Organization Programme on Man and the Biosphere (UNESCO, 2006) and one Ramsar site defined under the Ramsar Convention (RAMSAR, 2010) in the TSLRZ. The land cover types in the study area mainly include woodland, wasteland, water body, and grassland (JICA, 1999; Eng and Ouch, 2006; Mahood et al., 2020; Chen et al., 2022) and man-made land use, such as rice-growing farmland and fragmented woodland, due to human interference (ADB, 2019; Mahood et al., 2020; Chen et al., 2022). In recent years, the TSL area has been vigorously reclaimed to increase rice production, with the use of high-yielding rice varieties and an expansion of irrigation infrastructure (ADB, 2019; Mahood et al., 2020). In addition, large net forest loss is occurring in this area due to logging and agricultural farming (Chen et al., 2022).

2.2 Data

(1) Landsat dataset
The surface reflectance (SR) products of Landsat images, with a spatial resolution of 30 m, were used for land cover classification in the TSLRZ during the period from 1988 to 2020. These products were generated based on Landsat-5 TM, Landsat-7 ETM+ and Landsat-8 OLI images. All products were subjected to atmospheric correction and geometric fine correction, and the cloud cover of a single-scene image was less than 5%. The 78 SR products were downloaded from the USGS website (https://www.usgs.gov/land-resources/nli/landsat). The clear images mainly appeared in January, February, March and December in the TSLRZ due to cloudy and rainy weather throughout the year. To ensure that the various land cover types were identifiable and that the acquired data had similar hydrological and phenological conditions in each year, land cover classification in our study was performed using images from January to March in each year. Notably, we did not gather clear images covering the entire study area for all the seven years (1990, 1997, 2004, 2008, 2010, 2012 and 2013) because the cloud cover was too high. Three Landsat images covered the entire study area, with paths/rows of 127/051, 126/051, and 126/052. Specific information on the satellite data is listed in Table 1.
Table 1 Specific information on the satellite data collected for land cover classification
Year Image name, Cloud cover (%) Image name, Cloud cover (%) Image name, Cloud cover (%)
1988 LT51270511988044BKT01, 0.02 LT05L1TP12605119880222, 0.99 LT05L1TP12605219880309, 0.35
1989 LT41270511989054XXX02, 0.15 LT05L1TP12605119890224, 1.33 LT05L1TP12605219890224, 0.04
1991 LT51270511991036BKT00, 0.14 LT05L1TP12605119910129, 0.01 LT05L1TP12605219910129, 0
1992 LT51270511992087BKT01, 1.39 LT05L1TP12605119920304, 0.01 LT05L1TP12605219920304, 0.01
1994 LT51270511994028BKT01, 0.01 LT05L1TP12605119940105, 1.31 LT05L1GS12605219940121, 0.46
1995 LT51270511995031BKT00, 0 LT05L1TP12605119950209, 0.03 LT05L1TP12605219950209, 0
1996 LT51270511996002CLT00, 0.03 LT05L1TP12605119960127, 0.39 LT05L1TP12605219960127, 0
1998 LT51270511998007BKT00, 0.04 LT05L1TP12605119980321, 0.59 LT05L1TP12605219980305, 2
1999 LT51270511999026BKT00, 0.4 LT05L1TP12605119990119, 0.88 LT05L1TP12605219990220, 0
2000 LT51270512000013BKT00, 1.08 LT05L1TP12605120000326, 1.74 LT05L1TP12605220000326, 0
2001 LT51270512001047BKT00, 0.07 LT05L1TP12605120010124, 1.88 LT05L1TP12605220010124, 0
2002 LE71270512002010SGS00, 0.02 LE71260512002051BKT00, 0.02 LE07L1TP12605220020220, 0
2003 LE71270512003029SGS00, 0.01 LE71260512003038SGS00, 0.31 LE07L1TP12605220030207, 0
2005 LT51270512005042BKT00, 0.7 LT05L1TP12605120050119, 0.14 LT05L1TP12605220050119, 0.02
2006 LT51270512006045BKT00, 2.02 LT05L1TP12605120060207, 0.14 LT05L1TP12605220060207, 0
2007 LT51270512007032BKT00, 0.01 LT05L1TP12605120070125, 0.13 LT05L1TP12605220070210, 0.07
2009 LT05L1TP12705120090121, 0.98 LT05L1TP12605120090114, 0.12 LT05L1TP12605220090114, 0
2011 LT05L1TP12705120110127, 1.57 LT05L1TP12605120110120, 6.18 LT05L1TP12605220110120, 3.12
2014 LC81270512014035LGN00, 2.16 LC81260512014028LGN00, 0.09 LC08L1TP12605220140128, 0.17
2015 LC81270512015022LGN00, 0.01 LC81260512015015LGN00, 0 LC08L1TP12605220150115, 0.01
2016 LC81270512016089LGN00, 3.67 LC81260512016034LGN00, 0 LC08L1TP12605220160203, 0
2017 LC08L1TP12705120170212, 0.04 LC08L1TP12605120170205, 0.18 LC08L1TP12605220170205, 0.78
2018 LC08L1TP12705120180215, 0.22 LC08L1TP12605120180312, 2.63 LC08L1TP12605220180312, 1.12
2019 LC08L1TP12705120190117, 0.01 LC08L1TP12605120190211, 0.07 LC08L1TP12605220190211, 0.01
2020 LC81270512020020LGN00, 0.04 LC81260512020013LGN00, 0.01 LC81260522020013LGN00, 0.03
(2) Dataset for generating training and validation samples
One hundred and twenty-four Gaofen-1 (GF-1) and GF-2 images (30 in 2015, 28 in 2016, 38 in 2017, and 28 in 2018) were acquired from Suzhou Zhongke Tianqi Remote Sensing Technology Co., Ltd. The GF-1 images included panchromatic images with a 2 m spatial resolution and multispectral images with an 8 m spatial resolution. The GF-2 images included panchromatic images with a 1 m spatial resolution and multispectral images with a 4 m spatial resolution. Sixty-two fused GF images covering the entire study area, with a resolution of at least 2 m, were produced from these images through principal component analysis (PCA) (Pal et al., 2007) and were used to generate training and validation samples for land cover classification.
Two thematic maps of land cover classification (JICA, 1999; Eng and Ouch, 2006; Arias, 2013) created by the Japan International Cooperation Agency in 1996 and the Spatial Analysis and Information Laboratory, Ministry of Agriculture, Forestry, and Fisheries, Cambodia, in 2005 were used to generate the training and validation samples, respectively. The former was created based on SPOT images with a spatial resolution of 10 m and Landsat images with a spatial resolution of 30 m. The latter was created based on aerial images at a 1:5,000 scale. Notably, the two maps were unpublished data and sourced from the work of Arias (2013).
Google Earth historical images, with a spatial resolution of 0.61-10 m, were another useful data source for generating training and validation samples for these years without thematic land cover classification maps and fused GF images. These historical images were produced by integrating aerospace and aerial images (such as QuickBird, IKONOS, SPOT and others). The time sequence of the Google Earth historical images from 1984 to the present and the historical images can be achieved through the Google Earth (GE) Platform.
(3) Digital elevation model (DEM)
The DEM used in this study was the Shuttle Radar Topography Mission (SRTM) DEM product (version 3.0). These data have a spatial resolution of 30 m and were obtained from https://gdex.cr.usgs.gov/gdex/. In the follow-up study, the elevation, slope and aspect of the TSLRZ were extracted from the DEM for land cover classification.
(4) Measured water level data
The water level changes of the Kratie gauging station (Figure 1) were used to represent the water level changes of the LMR mainstream. The selection of the Kratie gauging station mainly hinges on its historically statistical correlation within the TSL’s inundation range prior to the development of upstream mega dams (Wang et al., 2020; Gu et al., 2021), a timeframe that establishes the baseline for this study. In addition, it has continuous and complete water level data and is a nearer station to the TSL among all gauging stations in the LMR mainstream. The daily water level data of the Kratie station during the period from 1987 to 2020 were obtained from http://ffw.mrcmekong.org/. In subsequent analyses (trend analysis and correlation analysis), the daily water level was processed into the average water level during the dry season (abbreviated as WLd), the average water level during the wet season (abbreviated as WLw) and the annual average water level (abbreviated as WLa). Notably, before August 31, 2006, water level measurements per day were taken once, while after this date, two measurements per day were recorded, which may have led to the occasional occurrence of distinct fluctuations within a single day. These fluctuations on the daily scale have been ignored in this study, considering that the time scales of hydrological indicators analysed (WLd, WLw and WLa) are longer than the daily scale.

2.3 Methods

2.3.1 Remote sensing-based mapping of land cover

(1) Land cover classification system in the TSLRZ
The land cover classification system of the study area (Table 2) was mainly determined based on the real-world conditions in the TSLRZ, referring to the two abovementioned thematic maps and the study by Mahood et al. (2020). The land cover types include farmland, woodland, wasteland, degraded land, water body and grassland. Notably, farmland includes rice fields and village crops that are under cultivation or not under cultivation, and wasteland is defined as a land cover type that is suitable for agricultural farming, which is sometimes cultivated and sometimes uncultivated. Woodland refers to natural vegetation areas with vegetation coverage greater than or equal to 60%. Degraded land is fragmented woodland and refers to natural vegetation areas with a vegetation coverage of less than 60%. The sample areas are marked with red wireframes (Table 2). The base maps were composed of the 5, 4 and 3 bands of Landsat images collected in 1988.
Table 2 Land cover classification system in the Tonle Sap Lake riparian zone
Land cover types Sample areas Land cover types Sample areas
Water body Woodland
Wasteland Degraded land
Grassland Farmland
(2) Land cover classification based on random forest (RF)
The RF model is an ensemble learning algorithm that has been widely used for land cover classification (Breiman, 2001; Liaw and Wiener, 2002; Yigez et al., 2023). A random forest is a combination of many trees. The number of trees (ntree) in the forest and the number of predictor variables in the random subset at a node of a tree (mtry) are the only two user-defined parameters needed for an RF model (Breiman, 2001; Liaw and Wiener, 2002). In the present study, b1, b2, b3, b4, b5, b6, longitude, latitude, elevation, slope, aspect and vegetation coverage were input as independent variables, which were extracted from the six bands of the SR products of Landsat images, graticules, DEM and other auxiliary data, respectively. The land cover type was used as the dependent variable. The calculation of vegetation coverage was performed as described in previous studies (Li, 2003; Wang et al., 2023). The RF model was used for land cover classification in this study, as explained below.
First, a point dataset with 1,000,000 points covering the study area was generated, and then the sample dataset with 3000 points was constructed through two consecutive random samplings from the point dataset, which was conducted separately in each year. Second, the land cover type of each point in the sample dataset was assigned using visual interpretation based on the reference images, and its corresponding independent variables were linked to the point by the map coordinates. Third, all independent variables from each year were merged by utilizing 30×30 m grids covering the study area and then converted into a matrix as the predicted dataset. Fourth, the RF model was trained and validated over the sample dataset, and the verified model was applied to the predicted dataset. The value of ntree was set as 1000, and mtry was equal to the square root of the number of independent variables. Fifth, the classification results were postprocessed via a majority filter and visual interpretation. Sixth, the accuracy of the classification results was evaluated using producer’s accuracy and user’s accuracy based on a confusion matrix (Foody, 2002).

2.3.2 Measured metrics of land cover change

(1) Trend analysis
The trend line based on the univariate linear regression equation and the least squares method (Wen et al., 2006) was used to analyse the change trends in the land cover type areas. The formula of the linear slope of the trend line (k) is as follows:
k=i=1nxi×i1n(i=1nxi)(i=1ni)i=1ni21n(i=1ni)2
(1)
where a positive k indicates an increasing trend and a negative k indicates a decreasing trend. xi denotes the area of a land cover type in year i. n denotes the total number of years. The significance and the significance level are denoted by p and α, respectively, and α was set to 0.05.
Considering that the operation schedule of the mega dams in the mainstream of the LMR and the hydrological conditions (water level, runoff) of the mainstream showed distinct changes in 2008 (Ji et al., 2018; Yun et al., 2020), we focused on the changes in land cover in two periods: the prechange period (1988-2007) and the postchange period (2009-2020). Importantly, 2008 was not included in the study periods due to the lack of clear Landsat images available to generate a land cover classification map.
(2) Transition matrix
A transition matrix (Takada et al., 2010) is used to reflect the mutual conversion process between various land cover types from the initial stage to the final stage in an area. The matrix is as follows:
Sij=[s11s12s1ns21s22s2nsn1sn2snn]
(2)
where s is the area and n is the number of land cover types. i and j (i, j=1, 2,..., n) are the land cover types at the initial and final stages, respectively. sij is the area of the initial i land type converted to the final j land type. Notably, the transition matrix was calculated based on land cover classification maps in vector format and then presented as a Sankey diagram in this study.
(3) Hot spot analysis
Hot spot analysis (Mitchel, 2010) is a spatial clustering method that can measure the change dynamics of the land cover in a given area and period. High values and low values with statistical significance are called hot spots and cold spots, respectively. The hot spots and cold spots correspond to areas with dramatic land cover changes and areas with minor land cover changes, respectively. In addition, the change frequency map of the land cover was used as the base map for hot spot analysis in this study.

2.3.3 Analysis of mainstream water level change and its relationship with the area of various land cover types

The trend line (see Formula 1 for details) was also used to analyse the trend of the WLd, WLw and WLa. As 80%-90% of the discharge of the LMR occurs from May to October (Delgado et al., 2010), we considered the wet period to be from May to October, and the other months were considered the dry season.
The Spearman correlation coefficient (Spearman, 2010) was used to analyse the relationship between the areas of different land cover types and the mainstream water level. The Spearman correlation coefficient does not require the original variables to have a certain distribution, and its value is between -1 and 1. The greater the absolute value is, the stronger the correlation. Suppose the datasets X(xi) (i=1, 2, 3, …, n) and Y(yi) (i=1, 2, 3, …, n) represent the land cover type area and mainstream water level, respectively. The datasets Q(qi) and S(si) are the level datasets of datasets X and Y, respectively. The formula to calculate the Spearman coefficient (ρ) is as follows:
ρ=i=1n(qiq¯)(sis¯)(i=1n(qiq¯)2)(i=1n(sis¯)2)
(3)
where ρ indicates the degree of correlation. qi and si denote the levels of samples xi and yi in the datasets, respectively. q¯ and s¯ denote the means of the level datasets Q and S, respectively. n denotes the total number of samples. The significance and the significance level are denoted by p and α, respectively, and α was set to 0.05.
Notably, in the correlation analysis, the areas of various land cover types in the current year corresponded to the average water levels of the mainstream (WLd, WLw and WLa) of the previous year, indicating a time lag response of the lake water level and land cover types to the mainstream water level, as well as the acquisition time of the Landsat images at the beginning of each year.

3 Results

3.1 The accuracy evaluation of the annual scale land cover classification dataset

According to the evaluation results (Figure 2), the producer’s accuracy and user’s accuracy of a single land cover type were greater than 80% for all land cover classification maps during the period from 1988 to 2020. Among these land cover types, the two precision indicators in farmland and water bodies were generally greater than those in other land use types, followed by woodland and grassland, and the indicators in wasteland and degraded land were relatively low.
Figure 2 The accuracy evaluation results of the annual scale land cover classification maps during 1988-2020

Full size|PPT slide

3.2 Spatiotemporal changes in land cover

3.2.1 Spatial dynamics of land cover

As shown in Figure 3, the distribution of land cover in the TSLRZ was characterized by an obvious ringed belt structure. Woodland, wasteland and farmland were distributed in the first, second and third belts, respectively, from the inside to the outside around the TSL, and the other land cover types were scattered among the three belts. The spatial pattern of land cover in the TSLRZ during the period from 1988 to 2020 showed the following changes.
Figure 3 Spatial patterns of land cover in the Tonle Sap Lake riparian zone during 1988-2020. The last picture is a schematic diagram of the belt distribution.

Full size|PPT slide

First, with the expansion of farmland towards the lakeside, the continuous woodland and grassland areas, especially in the southeastern part, and the wasteland in the second belt were obviously fragmented and gradually disappeared. Second, interactive conversion often occurred between wasteland and farmland. After 2006, the conversion gradually decreased. Third, the conversion of woodland and degraded land was often interactive. Clearly, more woodland was converted to degraded land in 2016 and 2020 than in other years. In particular, most of the degraded land located on both sides of the lake along the northwest‒southeast axis, which was converted from woodland in 2017, was not restored to woodland until 2020. Fourth, the distribution of the water body areas fluctuated greatly, and its area in 1996, 2001, 2009 and 2014 was larger than that in other years. Notably, some water body areas were reclaimed and converted to farmland after 2015, especially around the Tonle Sap River.
Figure 4 shows that in the period from 1988 to 2020, the area of land cover change hot spots in the TSLRZ was 4337.47 km2, accounting for 28.63% of the study area, and most of these hot spots were located in the second belt, the lakeside and the southeastern part. In the prechange period, the area of the hot spots was 3773.93 km2, accounting for 24.91% of the study area, and these hot spots were mostly located in the second belt. In the postchange period, the area of the hot spots was 4367.48 km2, accounting for 28.83% of the study area, and these hot spots were mainly located along the lakeside and in the southeastern part.
Figure 4 Hot spots of land cover changes in the Tonle Sap Lake riparian zone in the whole study period (1988-2020), the prechange period (1988-2007) and postchange period (2009-2020)

Full size|PPT slide

3.2.2 Temporal dynamics of land cover

Figure 5 shows that the proportion of farmland area in the study area has always been higher than that of other land use types, with values between 39.30% and 58.17%, followed by that of woodland. The area proportions of these two land cover types were greater than 70% in all years except 1996. In 1988, the proportions of farmland, woodland, wasteland, degraded land, water body and grassland areas in the study area were 41.72%, 32.61%, 13.14%, 5.42%, 4.33% and 2.78%, respectively. By 2020, the proportions of farmland and degraded land areas had increased by 16.46% and 4.40%, respectively, while the proportions of wasteland, woodland, grassland and water body areas decreased by 10.23%, 7.98%, 2% and 0.66%, respectively.
Figure 5 Structure of land cover in the Tonle Sap Lake riparian zone during 1988-2020

Full size|PPT slide

Figure 6 shows that during the period from 1988 to 2020, the area of farmland significantly increased, the areas of wasteland, grassland and woodland significantly decreased, and the areas of water bodies and degraded land did not significantly change. Compared with their areas in the prechange period, the increase in farmland and the decrease in grassland in the postchange period became significant, the change in the area of water body transitioned from a significant increase to a significant decrease, and the change in the area of degraded land transitioned from a nonsignificant decrease to a nonsignificant increase.
Figure 6 Change trends of various land cover types in the Tonle Sap Lake riparian zone. k denotes the linear slope. Subscripts 1, 2 and 3 correspond to the prechange period (1988-2007), postchange period (2009-2020) and the whole study period (1988-2020), respectively. p denotes significance, the significance level is set to 0.05, and * indicates p ≤ 0.05.

Full size|PPT slide

3.2.3 Conversion of land cover types

As shown in Figure 7a, from 1988 to 2020, the land cover type with the largest added area was farmland, the land cover type with the largest reduction was woodland, and the conversion type with the largest area was the conversion from wasteland to farmland. Figures 7b and 7c show that, first, the land cover type with the largest added area was farmland for both periods. However, the added farmland in the postchange period was approximately 2.02 times that in the prechange period. Second, the land cover type with the largest reduction was woodland, and 47.64% and 41.58% of the reduced part was converted to water bodies and degraded land in the prechange period. The land cover type with the largest reduction was water body, and 55.75% and 26.89% of the reduced part was converted to farmland and woodland in the postchange period. Third, the conversion type with the largest area was from wasteland to farmland in the prechange period, while it was from water body to farmland in the postchange period. Fourth, the area of woodland converted to degraded land was much larger than the area of woodland converted to farmland in each period. In addition, the area of woodland converted to degraded land in the postchange period was 1.78 times that in the previous period.
Figure 7 Sankey diagrams for conversions of the areas of various land cover types in the Tonle Sap Lake riparian zone from 1988 to 2020 (a), from 1988 to 2007 (b), and from 2009 to 2020 (c)

Full size|PPT slide

3.3 Changes in mainstream water level and its correlation with the area of land cover type

3.3.1 Changes in mainstream water level

Figure 8a shows that the highest daily water level was 21.94 m in 1996, and the lowest daily water level was 4.24 m in 1993. In addition, the daily water levels in 1998, 2015, 2016, 2019 and 2020 were lower than those in other years. Figures 8b-8d show that WLd and WLa significantly increased in the period from 1987 to 2007, while they did not significantly change in the period from 2008 to 2020. WLw showed no significant change trend in these two periods. Additionally, compared with that in the period from 1987 to 2007, the multiyear average of WLd and WLa in the period from 2008 to 2020 increased by 6.31% and 0.46%, respectively, and the multiyear average of WLw decreased by 2.72%. The changes in the multiyear averages of WLd and WLw indicated that the water level of the mainstream increased in the dry season and decreased in the wet season, and this phenomenon was related to the regulating effect of the upstream mega dams. These results can be supported by previous studies (Lu and Chua, 2021; Chua et al., 2022), although there have been some differences in the hydrological indicators analysed and time series.
Figure 8 Changes in the water level of the Lancang-Mekong River mainstream. k denotes the linear slope. Subscripts 1 and 2 correspond to the two periods of 1987-2007 and 2008-2020, respectively, which indicate the period before and the period after impoundment of the Xiaowan hydropower station, respectively. p denotes significance, the significance level is set to 0.05, and * indicates p ≤ 0.05.

Full size|PPT slide

3.3.2 Correlation between land cover type area and mainstream water level

The results of the Spearman correlation analysis (Figure 9 and Table 3) showed that the water body area was significantly related to WLd, WLw and WLa; the area of degraded land was significantly related to WLd and WLa; and the areas of two other land cover types (wasteland and grassland) were significantly related only to WLd. Compared with the prechange period, in the postchange period, the correlations between the areas of the four land cover types (water body, degraded land, wasteland and grassland) and WLd, WLw and WLa were reduced to varying degrees and did not reach significant levels. The results further indicated a weak linkage between the interannual change in local land cover and the mainstream water level and the operation of upstream mega dams.
Figure 9 Scatter plots of the areas of various land cover types in the Tonle Sap Lake riparian zone and the mainstream water levels. WLd, WLw and WLa represent the same meanings as those in Figure 8.

Full size|PPT slide

Table 3 Spearman correlation coefficient of the areas of various land cover types in the Tonle Sap Lake riparian zone and the mainstream water levels
Water
levels
Statistics Study
periods
Land cover types
Water body Wasteland Grassland Woodland Degraded land Farmland
WLd r 1988-2007 0.75* -0.64* -0.54* -0.15 -0.52* -0.02
2009-2020 0.29 -0.11 -0.18 0.21 -0.38 -0.29
p 1988-2007 0.00* 0.01* 0.02* 0.57 0.03* 0.93
2009-2020 0.44 0.78 0.65 0.59 0.31 0.44
WLw r 1988-2007 0.68* -0.37 -0.29 -0.16 -0.42 -0.09
2009-2020 0.28 -0.18 -0.13 0.30 -0.62 -0.17
p 1988-2007 0.00* 0.14 0.26 0.53 0.09 0.72
2009-2020 0.46 0.64 0.73 0.43 0.08 0.67
WLa r 1988-2007 0.75* -0.44 -0.36 -0.13 -0.52* -0.13
2009-2020 0.28 -0.18 -0.22 0.22 -0.53 -0.17
p 1988-2007 0.00* 0.08 0.16 0.63 0.03* 0.63
2009-2020 0.46 0.64 0.58 0.58 0.14 0.67
Note: WLd, WLw and WLa denote the average water level during the dry season and the wet season and the annual average water level; r denotes the Spearman correlation coefficient; p denotes significance, the significance level is set to 0.05, and * indicates p ≤ 0.05.

4 Discussion

4.1 Notable land cover changes after 2008

Three notable land cover changes were observed in the TSLRZ after 2008 based on the annual scale land cover classification dataset, namely, accelerated expansion of farmland, intensified fragmentation of woodland and significant shrinkage of water bodies (Figures 6 and 7). First, the increasing trend of the farmland area became significant in the postchange period (Figure 6), and the added farmland in the postchange period was approximately 2.02 times that in the prechange period (Figure 7). This result is supported by the study of Mahood et al. (2020), who observed that the added area of farmland was 441 km2 in the A-TSLRZ from 1993 to 2008, while it reached 2356 km2 from 2008 to 2018. Second, the area of woodland converted to degraded land in the postchange period was 1.78 times that in the previous period (Figure 7). This noticeable change was in sound agreement with the results of the study by Chen et al. (2022), who reported that forests have been rapidly fragmented in the A-TSLRZ since 1992 and that hardly any intact forests remained after 2010 due to logging and agricultural farming. Third, the change in the area of the water body significantly decreased in the postchange period (Figure 6). This is consistent with the results of several previous studies (Ji et al., 2018; Wang et al., 2020; Ng and Park, 2021). For example, Ji et al. (2018) found that the multiyear average area of the LWB in the TSL decreased by 1.5% in the period from 2008 to 2014 compared with the period from 2000 to 2007. Wang et al. (2020) also reported that the monthly average area of the LWB in the TSL changed only slightly in the period from 1988 to 2000 but showed a significant decreasing trend in the period from 2000 to 2018. Notably, compared with the results of these studies, the water body change detected in this study was more significant because our study did not include the smallest lake body, which is regarded as fixed.
Furthermore, some special phenomena have been documented based on the annual scale dataset. First, the interactive conversion of farmland and wasteland occurred often before 2006, and the areas of occurrence were unstable (Figures 3 and 7). This result may have been related to the local industrial structure. When local tourism and fishing are booming, labour tends to shift from agriculture to these industries, and vice versa. After 2006, this phenomenon was rare, which further reflected the acceleration of farmland expansion. Second, the conversion of woodland and degraded land was interactive, and the conversion was more likely to occur in the drier period. For example, 2015 and 2019 were typical drought years in the TSL area (Frappart et al., 2018; Sok et al., 2021). Correspondingly, more woodland was converted to degraded land in 2016 and 2020 than in other years (Figure 3). In addition, in contrast to the above interactive conversion caused by drought, most of the degraded land located on both sides of the lake, which was converted from woodland in 2017, was not restored to woodland until 2020. This difference can be ascribed to fires. According to the report by Lovgren (2020), massive fires burned as much as one-third of the 750,000 acres of flooded forest in the TSL area in 2016.

4.2 Response of land cover change to mainstream water level

The correlations between the water body area in the TSLRZ and the mainstream water levels (WLd, WLw and WLa) were statistically significant during the prechange period (1988-2007) but insignificant during the postchange period (2009-2020) (Figure 9 and Table 3). This result indicated that the water body shrinkage in the last 10 years (Figure 6) was not dominated by the mainstream water level and was still affected by other factors. This result is supported by the study of Ji et al. (2018), who reported that the decrease in the water area of the TSL during the dry season after 2008 could not be attributed solely to the change in mainstream runoff. The changes in the correlations between the water body area and the mainstream water levels during the postchange period can be attributed to the following three aspects. First, according to previous studies (Morovati et al., 2021; 2023; Chua et al., 2022), the flood pulse of the LMR has declined in recent years due to the development and operation of upstream mega dams in the upstream, climate variations, local irrigation and channel incision (caused by declining sediments combined with sand-mining operations). This may be an important reason for the shrinkage of water body areas. Second, the components of the water balance model for this lake-floodplain system mainly include the LMR mainstream, TSL’s own tributaries, rainfall, evaporation, groundwater infiltration and overland flow via floodplain (Kummu et al., 2014; Morovati et al., 2023), which suggested that the change in other components (excluding the mainstream) can also lead to the change in the water body areas in the TSLRZ. Daly et al. (2020) noted that the average temperature in the TSL increased by 0.03°C/year in the period from 1988 to 2018, which meant that evaporation increased in the area. The study of Ng and Park (2021) showed that tributary discharge and rainfall, which provide approximately 50% of the supply of lake water (the proportion of the supply may vary in different periods), did not display any significant trend over the investigated period (1980 to 2018) and were not the primary factors contributing to the shrinking of the lake. Third, based on the results of our study, the water body areas have been subjected to greater human interference during the postchange period. Under natural conditions, water body areas will be transformed into natural land cover types, such as woodland and grassland. However, approximately 55.75% of the reduced water body areas in 2009 were converted into farmland in 2020 (Figure 7), especially around the Tonle Sap River (Figure 3). This reflects that local agricultural farming interfered with the original response mode between the water body area of the TSLRA and the mainstream water level. Moreover, in the postchange period, the area of added farmland was approximately 2.02 times that in the previous period (Figure 7), which meant that the irrigation water consumption increased accordingly and would also interfere with the original response mode. Based on the above results, it can be concluded that the changes in water body areas in the TSLRA were dominated by mainstream water levels during the prechange period but were influenced by multiple factors together during the postchange period.
The other land cover types (farmland, woodland, wasteland, degraded land, and grassland) in the TSLRZ indirectly responded to the mainstream water level through the conduction of the lake water level. In the prechange period, the area of degraded land was significantly related to WLd and WLa, and the areas of two other land cover types (wasteland and grassland) were significantly related only to WLd (Figure 9 and Table 3), in which the correlation between the area of wasteland and WLd was strong. However, as shown in Figure 7b, most changed wasteland was converted into farmland; therefore, the causality of this statistically significant correlation was weak. In addition, the correlations between the areas of two land cover types (degraded land and grassland) and mainstream water levels were weak during the prechange period but reduced to varying degrees and did not reach significant levels during the postchange period. These fluctuations in both the areas of these land cover types and the mainstream water levels may relate to the more intense human activities, as well as more frequent and higher extreme changes in the climate in the postchange period.
According to the above analysis, the land cover changes in the last 30 years were mainly the conversion of nonagricultural land to agricultural land by local reclamation and logging, which stressed the importance of agriculture and the imperative of food security in Cambodia. The main driving force of land cover changes was the ever-growing population in Cambodia. The population of Cambodia changed from 7,975,597 in 1987 to 16,767,842 in 2022 (data sourced from the World Bank public database), and the demand for food continues to increase accordingly. Another driving force may be international market demand. The policy document on promoting rice production and rice export issued by Cambodia (Royal Government of Cambodia, 2010) highlighted how to turn Cambodia into a “rice basket” and a major global rice exporter and proposed an export target of 1 million tons in 2015. Moreover, Cambodia’s rice exports have accelerated since the launch of the China-ASEAN Free Trade Area in 2010 (Prom, 2016). In addition, although the reduction in the flood pulse of the LMR was not a driving force, it indirectly provided convenient conditions for agriculture activities to a certain extent. Usually, human activities such as reclamation and logging are constrained by higher flood levels and wider inundation ranges. After the flood pulse weakened, the range of water level changes became relatively narrow. This change was conducive to the development and utilization of land in the TSLRZ, thereby increasing the area of farmland.

4.3 Limitations and way forward

The relationship between the water level change of the mainstream and the land cover change in the TSLRZ, as identified in this study, represents a linkage relationship rather than a quantitative direct relationship. This is due to the complexity of the indirect effect of the mainstream water level on the land cover and the lack of relevant hydrological factor data in the TSLRZ. A quantitative relationship between them needs further investigation through sequential analysis of the two quantitative relationships. The first is between the land cover change in the TSLRZ and the hydrological factor change in the TSLRZ, while the second is between the hydrological factor change and the water level change of the mainstream. The hydrological factors, which include but are not limited to inundation frequency, inundation period, and humidity, are located in the same two-dimensional space as land cover and directly affect its dynamics.
The interaction mechanism between LMR and TSL is indeed complex due to the influence of the annual reverse flow phenomenon. The reverse flow occurs when the water level of the LMR is higher than that of the TSL, which will enlarge the TSL inundation range and impact the land cover changes in the TSLRZ. Currently, scholars have deeply studied the interaction mechanism (Kummu et al., 2014; Li et al., 2019; Morovati et al., 2021; 2023; Chua et al., 2022). In contrast to previous studies, the focus of this study is to explore the changes in the ecological environment of the TSL when the hydrological regime of the LMR mainstream underwent distinct changes and its relationship with the water level of the LMR. In this analysis, the land cover changes in the TSL were used to represent its ecological changes, while the water level changes at the Kratie gauging station were used to represent the water level changes in the LMR mainstream. The selection of the Kratie gauging station mainly hinges on its historically statistical correlation within the TSL’s inundation range prior to the development of upstream mega dams, a timeframe that establishes the baseline for this study. Notably, the water level fluctuations at the gauging station do not necessarily indicate water body area changes in the TSLRZ, and the flood pulse at this station does not necessarily trigger the onset of reverse flow. Therefore, the linkage relationship obtained in this study can only reflect the localized aspects of the interaction mechanism between LMR and TSL. Future research to comprehensively decipher this mechanism necessitates a multifaceted approach and integrating additional hydrological indicators.
While we do not advocate ignoring the future implications of the mainstream hydrological regime, the land cover changes that we have documented and their correlations with the mainstream water level showed that the influences of local reclamation and logging on the land cover changes were far greater than those of the mainstream water level, greatly threatening the ecological environment in the TSLRZ. The continuous expansion of farmland may be unsustainable in the long run, endangering regional food security. Therefore, it is of great importance for government agencies, relevant investors and development agencies to take urgent, coordinated action to address the issue of the conversion of nonagricultural land to agricultural land. By doing so, it is conducive to mitigating possible risk in the future of the TSL and the livelihoods of those residing in the region.
We propose three measures to address the issue identified here. First, areas that are already zoned for protection, such as the three Biosphere Reserve core areas and one Ramsar site, must be prohibited for reclamation and logging and strictly protected so that they can provide ecosystem services to the surrounding communities. Second, secondary protected areas should be built in some areas that have already experienced extensive reclamation and logging, such as the areas located on both sides of the lake along the northwest‒southeast axis and areas around the Tonle Sap River, and implement ecological restoration actions such as reverting reclaimed land to its natural state, water source protection and waterway dredging. Third, local government agencies should establish a data sharing platform. Specifically, this platform should include different types of data, such as hydrological data, meteorological data, precious species data and grain production and export data. Through data sharing, more scholars, research institutions and environmental protectors can participate in the ecological protection of the TSL.

5 Conclusions

In this study, we created an annual scale land cover classification dataset with a 30 m spatial resolution in the TSLRZ during the period from 1988 to 2020, and the producer’s accuracy and user’s accuracy of a single land cover type were greater than 80%. On this foundation, this study clarified the interannual spatiotemporal changes in land cover in the TSLRZ and assessed their relationship with mainstream water level changes in the LMR. The main conclusions of the study were as follows:
(1) Farmland was the largest land cover type, and its expansion towards the lakeside was the most notable change in the TSLRZ in the last 30 years. During this period, the farmland area showed a significant increase, the wasteland, grassland and woodland areas all showed significant decreases, and the water body and degraded land areas showed stage-specific characteristics.
(2) There were three notable land cover changes in the TSLRZ after 2008, namely, accelerated farmland expansion, intensified woodland fragmentation and significant water body shrinkage. During the postchange period (2009-2020), the water body area significantly decreased, and the areas of added farmland and fragmented woodland were 2.02 and 1.78 times those in the prechange period (1988-2007), respectively.
(3) The responses of the water body, degraded land, wasteland and grassland areas to the mainstream water levels weakened in the TSLRZ after 2008. Their correlations were reduced and did not reach significant levels, indicating a weak linkage between the interannual change in local land cover and the mainstream water level and the operation of upstream mega dams.
(4) The land cover changes in the TSLRZ mainly resulted from local reclamation and logging and were driven by the ever-growing population in Cambodia and the international market demand.
The results above can provide a new scientific basis for the transboundary influence analysis of hydrological change, and they can be used to support resource utilization and ecological protection in the riparian zone of river-connected lakes under climate change and cascade dam building.

参考文献

[1]
Arias M E, 2013. Impacts of hydrological alterations in the Mekong Basin to the Tonle Sap ecosystem[D]. Christchurch, New Zealand: University of Canterbury.
[2]
Asian Development Bank (ADB), 2019. Kingdom of Cambodia: Preparing the irrigated agriculture improvement project. https://www.adb.org/sites/default/files/project-documents/51159/51159-001-tacr-en.pdf.
[3]
Binh D V, Kantoush S A, Saber M et al., 2020. Long-term alterations of flow regimes of the Mekong River and adaptation strategies for the Vietnamese Mekong Delta. Journal of Hydrology: Regional Studies, 32: 100742.
[4]
Breiman L, 2001. Random Forest. Machine Learning, 45: 5-32.
[5]
Chen A F, Chen A P, Varis O et al., 2022. Large net forest loss in Cambodia’s Tonle Sap Lake protected areas during 1992-2019. Ambio, 51(8): 1889-1903.
Historical land-use practices have caused forest loss in Cambodia’s Tonle Sap Lake area (TSLA), the largest freshwater lake in Southeast Asia. However, it remains unclear if this deforestation trend had continued since 2001 when the land was designated as protected areas. Using satellite imagery, we investigated forest conversion flows and fragmentation patterns in the TSLA for 1992–2001, 2001–2010, and 2010–2019, respectively. Results show substantial forest losses and fragmentations occurring at the lower floodplain where the protected areas are located until 2010, with some forest regain during 2010–2019. The land conversions indicated that forest clearing and agricultural farming were the primary causes for observed extensive forest loss during 1992–2010. Hence, despite the creating of protected areas in 2001, our findings reveal the persistence of alarming forest loss in the TSLA until 2010. On the other hand, while net forest loss has stopped after 2010, forest regain during 2010–2019 is way too small to restore the region’s total forest area to even the level when the protected areas were established. Thus, more effective planning and implementations of forest management and restoration policies are needed for the TSLA.
[6]
Chua S D X, Lu X X, Oeurng C et al., 2022. Drastic decline of flood pulse in the Cambodian floodplains (Mekong River and Tonle Sap system). Hydrology and Earth System Sciences, 26(3): 609-625.
. The Cambodian floodplains experience a yearly flood pulse that is\nessential to sustain fisheries and the agricultural calendar. Sixty years of data, from 1960–2019, are used to track the changes to the flood pulse there.\nWe find that minimum water levels over 2010–2019 increased by up to 1.55 m at Kratie and maximum water levels decreased by up to 0.79 m at Prek\nKdam when compared to 1960–1991 levels, causing a reduction of the annual\nflood extent. Concurrently, the duration of the flooding season has\ndecreased by about 26 d (Kampong Cham) and 40 d (Chaktomuk), with the\nseason starting later and ending much earlier. Along the Tonle Sap River,\nthe average annual reverse flow from the Mekong to the Tonle Sap Lake has\ndecreased by 56.5 %, from 48.7 km3 in 1962–1972 to 31.7 km3 in\n2010–2018. As a result, wet-season water levels at Tonle Sap Lake\ndropped by 1.05 m in 2010–2019 compared to 1996–2009, corresponding to a 20.6 %\nshrinkage of the lake area. We found that upstream contributors such as\ncurrent hydropower dams cannot fully account for the observed decline in\nflood pulse. Instead, local anthropogenic causes such as irrigation and\nchannel incision are important drivers. We estimate that water withdrawal in\nthe Cambodian floodplains is occurring at a rate of (2.1 ± 0.3) km3 yr−1. Sediment decline and ongoing sand-mining operations\nhave also caused channel erosion. As the flood pulse is essential for the\necological habitats, fisheries and livelihoods of the region, its reduction\nwill have major implications throughout the basin, from the Tonle Sap system\nto the Vietnamese Mekong Delta downstream.\n
[7]
Daly K, Ahmad S K, Bonnema M et al., 2020. Recent warming of Tonle Sap Lake, Cambodia: Implications for one of the world’s most productive inland fisheries. Lakes and Reservoirs: Research and Management, 25(2): 133-142.
[8]
Dang H, Pokhrel Y, Shin S et al., 2022. Hydrologic balance and inundation dynamics of Southeast Asia’s largest inland lake altered by hydropower dams in the Mekong River basin. Science of the Total Environment, 831: 154833.
[9]
Delgado J M, Apel H, Merz B, 2010. Flood trends and variability in the Mekong River. Hydrology and Earth System Sciences, 14(3): 407-418.
. Annual maximum discharge is analyzed in the Mekong river in Southeast Asia with regard to trends in average flood and trends in variability during the 20th century. Data from four gauging stations downstream of Vientiane, Laos, were used, covering two distinct hydrological regions within the Mekong basin. These time series span through over 70 years and are the longest daily discharge time series available in the region. The methods used, Mann Kendal test (MK), ordinary least squares with resampling (OLS) and non-stationary generalized extreme value function (NSGEV), are first tested in a Monte Carlo experiment, in order to evaluate their detection power in presence of changing variance in the time series. The time series are generated using the generalized extreme value function with varying scale and location parameter. NSGEV outperforms MK and OLS, both because it resulted in less type II errors, but also because it allows for a more complete description of the trends, allowing to separate trends in average and in variability. Results from MK, OLS and NSGEV agreed on trends in average flood behaviour. However, the introduction of a time-varying scale parameter in the NSGEV allowed to isolate flood variability from the trend in average flood and to have a more complete view of the changes. Overall, results showed an increasing likelihood of extreme floods during the last half of the century, although the probability of an average flood decreased during the same period. A period of enhanced variance in the last quarter of the 20th century, estimated with the wavelet power spectrum as a function of time, was identified, which confirmed the results of the NSGEV. We conclude that the absence of detected positive trends in the hydrological time series was a methodological misconception due to over-simplistic models.\n
[10]
Eng C, Ouch V, 2006. TSBR Land cover map derived from the orthophoto map (unpublished data from the Spatial Analysis and Information Laboratory, Ministry of Agriculture, Forestry, and Fisheries,
[11]
Foody G M, 2002. Status of land cover classification accuracy assessment. Remote Sensing of Environment, 80(1): 185-201.
[12]
Frappart F, Biancamaria S, Normandin C et al., 2018. Influence of recent climatic events on the surface water storage of the Tonle Sap Lake. Science of The Total Environment, 636: 1520-1533.
[13]
Gu Z K, Zhang Y, Fan H, 2021. Mapping inter- and intra-annual dynamics in water surface area of the Tonle Sap Lake with Landsat time-series and water level data. Journal of Hydrology, 601: 126644.
[14]
Guan Y, Zheng F, 2021. Alterations in the water-level regime of Tonle Sap Lake. Journal of Hydrologic Engineering, 26(1): 05020045.
[15]
He C, Zhang J, Liu Z, et al., 2022. Characteristics and progress of land use/cover change research during 1990-2018. Journal of Geographical Sciences, 32(3): 537-559.

Land use/cover change (LUCC) is the foundation and frontier for integrating multiple land surface processes. This paper aims to systematically review LUCC research from 1990 to 2018. Based on qualitative and quantitative analyses, we delineated the history of LUCC research and summarized their characteristics and major progress at different stages. We also identified the main challenges and proposed future directions for LUCC research. We found that the number of publications on LUCC research and their total citations grew exponentially. The research foci shifted from the process of LUCC during 1990-2004 to the impact of LUCC during 2005-2013 and then to the sustainability of LUCC from 2014 onwards. Currently, LUCC research is facing theoretical, methodological and practical challenges ranging from integrating the framework of sustainability science, adopting emerging technologies to supporting territorial spatial planning. To move forward, LUCC research should be closely integrated with landscape sustainability science and geodesign and take the leading role in territorial spatial planning to achieve the related Sustainable Development Goals.

[16]
He D M, Liu H, Feng Y et al., 2016. Perspective on theories and methods study of transboundary water resources under the global change. Advances in Water Science, 27(6): 928-934. (in Chinese)
[17]
Huang Z Y, Wang J D, Zhang Y B, 2020. A brief analysis of the relationship between landscape change and human activities in Tonle Sap Lake Basin, Cambodia. China Rural Water and Hydropower, (4): 30-38. (in Chinese)
[18]
Ji X, Li Y, Luo X et al., 2018. Changes in the lake area of Tonle Sap: Possible linkage to runoff alterations in the Lancang River? Remote Sensing, 10(6): 866.
[19]
JICA, 1999. Cambodia reconnaissance survey digital data (unpublished data from the Japan International Cooperation Agency).
[20]
Kim S, Sohn H G, Kim M K et al., 2019. Analysis of the relationship among flood severity, precipitation, and deforestation in the Tonle Sap Lake Area, Cambodia using multi-sensor approach. KSCE Journal of Civil Engineering, 23: 1330-1340.
[21]
Kummu M, Tes S, Yin S et al., 2014. Water balance analysis for the Tonle Sap Lake-floodplain system. Hydrological Processes, 28(4): 1722-1733.
[22]
Lamberts D, 2008. The unintended role of the local private sector in biodiversity conservation in the Tonle Sap Biosphere Reserve, Cambodia. Local Environment, 13(1): 43-54.
[23]
Lambin E F, Turner B L, Geist H J et al., 2001. The causes of land-use and land-cover change: Moving beyond the myths. Global Environmental Change, 11(4): 261-269.
[24]
Li C W, You Z Q, Li A Q et al., 2019. Variation characteristics, influencing factors and hydrological conditions of the reverse flow from Mekong River to Tonle Sap Lake. Applied Ecology and Environmental Research, 17(6): 13875-13895.
[25]
Li D, Long D, Zhao J et al., 2017. Observed changes in flow regimes in the Mekong River basin. Journal of Hydrology, 551: 217-232.
[26]
Li M M, 2003. The method of vegetation fraction estimation by remote sensing[D]. Beijing, China: Institute of Remote Sensing Applications, Chinese Academy of Sciences. (in Chinese)
[27]
Liaw A, Wiener M, 2002. Classification and regression by random forest. R News, 2(3): 18-22.
[28]
Lin Z, Qi J, 2017. Hydro-dam-A nature-based solution or an ecological problem: The fate of the Tonle Sap Lake. Environmental Research, 158: 24-32.
[29]
Lovgren S, 2020. Cambodia’s biggest lake is running dry, taking forests and fish with it. https://www.nationalgeographic.com/science/article/cambodia-tonle-sap-lake-running-dry-taking-flooded-forest-fish.
[30]
Lu X X, Chua S, 2021. River discharge and water level changes in the Mekong River: Droughts in an era of mega-dams. Hydrological Processes, 35(7): e14265.
[31]
Mahood S P, Poole C M, Watson J E M et al., 2020. Agricultural intensification is causing rapid habitat change in the Tonle Sap floodplain, Cambodia. Wetlands Ecology and Management, 28(5): 713-726.
[32]
Meng Q J, Zang S Y, Song K S et al., 2018. A comparison of conditions of Tonle Sap Lake wetlands in dry and wet seasons in 1990 with those in 2016. Wetland Science, 16(6): 801-807. (in Chinese)
[33]
Mitchel A, 2010. The ESRI Guide to GIS Analysis (Volume 2): Spatial Measurements and Statistics. RedLands: ESRI Press.
[34]
Morovati K, Nakhaei P, Tian F et al., 2021. A machine learning framework to predict reverse flow and water level: A case study of Tonle Sap Lake. Journal of Hydrology, 603: 127168.
[35]
Morovati K, Tian F, Kummu M et al., 2023. Contributions from climate variation and human activities to flow regime change of Tonle Sap Lake from 2001 to 2020. Journal of Hydrology, 616: 128800.
[36]
Ng W X, Park E, 2021. Shrinking Tonle Sap and the recent intensification of sand mining in the Cambodian Mekong River. Science of the Total Environment, 777: 146180.
[37]
Pal S K, Majumdar T J, Bhattacharya A K, 2007. ERS-2 SAR and IRS-1C LISS III data fusion: A PCA approach to improve remote sensing based geological interpretation. ISPRS Journal of Photogrammetry and Remote Sensing, 61(5): 281-297.
[38]
Pokhrel Y, Burbano M, Roush J et al., 2018. A review of the integrated effects of changing climate, land use, and dams on Mekong River hydrology. Water, 10(3): 266.
The ongoing and proposed construction of large-scale hydropower dams in the Mekong river basin is a subject of intense debate and growing international concern due to the unprecedented and potentially irreversible impacts these dams are likely to have on the hydrological, agricultural, and ecological systems across the basin. Studies have shown that some of the dams built in the tributaries and the main stem of the upper Mekong have already caused basin-wide impacts by altering the magnitude and seasonality of flows, blocking sediment transport, affecting fisheries and livelihoods of downstream inhabitants, and changing the flood pulse to the Tonle Sap Lake. There are hundreds of additional dams planned for the near future that would result in further changes, potentially causing permanent damage to the highly productive agricultural systems and fisheries, as well as the riverine and floodplain ecosystems. Several studies have examined the potential impacts of existing and planned dams but the integrated effects of the dams when combined with the adverse hydrologic consequences of climate change remain largely unknown. Here, we provide a detailed review of the existing literature on the changes in climate, land use, and dam construction and the resulting impacts on hydrological, agricultural, and ecological systems across the Mekong. The review provides a basis to better understand the effects of climate change and accelerating human water management activities on the coupled hydrological-agricultural-ecological systems, and identifies existing challenges to study the region’s Water, Energy, and Food (WEF) nexus with emphasis on the influence of future dams and projected climate change. In the last section, we synthesize the results and highlight the urgent need to develop integrated models to holistically study the coupled natural-human systems across the basin that account for the impacts of climate change and water infrastructure development. This review provides a framework for future research in the Mekong, including studies that integrate hydrological, agricultural, and ecological modeling systems.
[39]
Prom S, 2016. Research on rice export trade of Cambodia[D]. Guilin, China: Guangxi Normal University. (in Chinese)
[40]
RAMSAR, 2010. The list of wetlands of international importance. http://www.ramsar.org/.
[41]
Royal Government of Cambodia, 2010. Policy paper on promotion of paddy rice production and export of milled rice. Supreme National Economic Council, Phnom Penh.
[42]
Siev S, Yang H, Sok T et al., 2018. Sediment dynamics in a large shallow lake characterized by seasonal flood pulse in Southeast Asia. The Science of the Total Environment, 631: 597-607.
[43]
Sok S, Chhinh N, Hor S et al., 2021. Climate change impacts on rice cultivation: A comparative study of the Tonle Sap and Mekong River. Sustainability, 13(16): 8979.
Climate change is unequivocal. Farmers are increasingly vulnerable to floods and drought. In this article, the negative impact of climate hazards on rice cultivation in the Tonle Sap and Mekong River influenced by climatic variability between 1994 and 2018 are analyzed. A cohort of 536 households from four Cambodian districts participated in household surveys designed to consider how various vulnerability factors interacted across this time series. It was found that: (i) The major climate hazards affecting rice production between 1994 and 2018 were frequent and extreme flood and drought events caused by rainfall variability; (ii) In 2018, extreme flood and drought occurred in the same rice cultivation cycle. The impact caused by each hazard across each region were similar; (iii) An empirical model was used to demonstrate that drought events tend to limit access to irrigation, impact rice production, and result in an increased prevalence of water-borne diseases. Flood events cause reduced rice production, damage to housing, and impede children from accessing education. The impact of drought events on rice production was found to be more severe than flood events; however, each climatic hazard caused physical, economic, social, and environmental vulnerabilities. It is recommended that sufficient human and financial resources are distributed to local authorities to implement adaptation measures that prepare rice farmers for flood and drought events and promote equitable access to water resources.
[44]
Spearman C, 2010. The proof and measurement of association between two things. International Journal of Epidemiology, 39(5): 1137-1150.
[45]
Sterling S M, Ducharne A, Polcher J, 2013. The impact of global land-cover changes on the terrestrial water cycle. Nature Climate Change, 3(4): 385-390.
[46]
Takada T, Miyamoto A, Hasegawa S F, 2010. Derivation of a yearly transition probability matrix for land-use dynamics and its applications. Landscape Ecology, 25: 561-572.
[47]
Tangdamrongsub N, Ditmar P G, Steele-Dunne S C et al., 2016. Assessing total water storage and identifying flood events over Tonle Sap basin in Cambodia using GRACE and MODIS satellite observations combined with hydrological models. Remote Sensing of Environment, 181: 162-173.
[48]
Trung N V, Choi J H, Won J S, 2013. A land cover variation model of water level for the floodplain of Tonle Sap, Cambodia, derived from ALOS PALSAR and MODIS data. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 6(5): 2238-2253.
[49]
Turner B L, Lambin E F, Reenberg A, 2007. The emergence of land change science for global environmental change and sustainability. Proceedings of the National Academy of Sciences of the United States of America, 104(52): 20666-20671.
Land change science has emerged as a fundamental component of global environmental change and sustainability research. This interdisciplinary field seeks to understand the dynamics of land cover and land use as a coupled human-environment system to address theory, concepts, models, and applications relevant to environmental and societal problems, including the intersection of the two. The major components and advances in land change are addressed: observation and monitoring; understanding the coupled system-causes, impacts, and consequences; modeling; and synthesis issues. The six articles of the special feature are introduced and situated within these components of study.
[50]
UNESCO, 2006. Biosphere Reserves: World Network. http://www.unesco.org/mab/.
[51]
Wang D, Liu Y, Zheng L, et al., 2023. Growing impacts of low-flow events on vegetation dynamics in hydrologically connected wetlands downstream Yangtze River Basin after the operation of the Three Gorges Dam. Journal of Geographical Sciences, 33(4): 885-904.

Wetland vegetation is intimately related to floodplain inundations, which can be seriously affected by dam operation. Poyang Lake is the largest floodplain wetland in China and naturally connected with the Yangtze River and the Three Gorges Dam (TGD) upstream. To understand the potential impacts of TGD on Poyang Lake wetlands, we collected remote sensing imagery acquired during dry season from 1987 to 2020 and extracted vegetation coverage data in the Ganjiang Northern-branch Delta (GND) and the Ganjiang Southern-branch Delta (GSD), using the Object-oriented Artificial Neural Network Regression. Principal components analysis, correlation analysis, and the random forest model were used to explore the interactions between vegetation extent in the two deltas and 33 hydrological variables regarding magnitude, duration, timing, and variation. The implementation of the TGD advanced and extended the low-flow periods in Poyang Lake. Vegetation coverage in the GND and GSD increased at the rates of 0.39 and 0.22 km2/year, respectively. The reservoir storage at the end of September accelerated the runoff recession in the GND and the GSD, making low-flow events more influential for vegetation dynamics and shortening the response time of vegetation to the water regime. This study provides an important reference for evaluating the impacts of dam engineering on downstream wetlands.

[52]
Wang H W, Qi Y, Lian X H et al., 2022. Effects of climate change and land use/cover change on the volume of the Qinghai Lake in China. Journal of Arid Land, 14(3): 245-261.

Qinghai Lake is the largest saline lake in China. The change in the lake volume is an indicator of the variation in water resources and their response to climate change on the Qinghai-Tibetan Plateau (QTP) in China. The present study quantitatively evaluated the effects of climate change and land use/cover change (LUCC) on the lake volume of the Qinghai Lake in China from 1958 to 2018, which is crucial for water resources management in the Qinghai Lake Basin. To explore the effects of climate change and LUCC on the Qinghai Lake volume, we analyzed the lake level observation data and multi-period land use/land cover (LULC) data by using an improved lake volume estimation method and Integrated Valuation of Ecosystem Services and Trade-offs (InVEST) model. Our results showed that the lake level decreased at the rate of 0.08 m/a from 1958 to 2004 and increased at the rate of 0.16 m/a from 2004 to 2018. The lake volume decreased by 105.40×108 m3 from 1958 to 2004, with the rate of 2.24×108 m3/a, whereas it increased by 74.02×108 m3 from 2004 to 2018, with the rate of 4.66×108 m3/a. Further, the climate of the Qinghai Lake Basin changed from warm-dry to warm-humid. From 1958 to 2018, the increase in precipitation and the decrease in evaporation controlled the change of the lake volume, which were the main climatic factors affecting the lake volume change. From 1977 to 2018, the measured water yield showed an "increase-decrease-increase" fluctuation in the Qinghai Lake Basin. The effects of climate change and LUCC on the measured water yield were obviously different. From 1977 to 2018, the contribution rate of LUCC was -0.76% and that of climate change was 100.76%; the corresponding rates were 8.57% and 91.43% from 1977 to 2004, respectively, and -4.25% and 104.25% from 2004 to 2018, respectively. Quantitative analysis of the effects and contribution rates of climate change and LUCC on the Qinghai Lake volume revealed the scientific significance of climate change and LUCC, as well as their individual and combined effects in the Qinghai Lake Basin and on the QTP. This study can contribute to the water resources management and regional sustainable development of the Qinghai Lake Basin.

[53]
Wang T, Liu C L, Du D B, 2021. Spatio-temporal dynamics of international freshwater conflict events and relations from 1948 to 2018. Acta Geographica Sinica, 76(7): 1792-1809. (in Chinese)

With global climate change and the rapid development of human society and economy, the contradiction between water supply and demand has become increasingly prominent in recent years, and the freshwater conflicts in international river basins have intensified, which has aroused widespread concern in academia. Here we analyzed the spatio-temporal dynamics of global freshwater conflicts (GFCs) over the last 70 years from the "event-relations" perspective, and establish a spatio-temporal database of GFCs from 1948 to 2018 based on data mining method and spatial analysis. The results show that: (1) The evolution of GFCs is a non-monotonic dynamic process with multi-dimensional characteristics of trend, mutation and volatility. The GFCs showed a general trend of fluctuating growth, with an obvious sudden change around 1987. (2) The GFCs are mainly composed of low-intensity conflicts, and the hydrological intervention and contention for resource ownership are the focus of conflicts. The number of conflicts caused by the construction of dams and other water conservancy projects increases significantly. South Asia, West Asia and East Africa are the leading forces driving the evolution of GFCs. (3) The pattern of GFCs has changed from single-center to multi-center, and there is a clear trend of spatial spread. However, the overall distribution pattern with more conflicts in the northern and eastern hemispheres and the pattern with less conflicts in the southern and western hemispheres is relatively stable. Along 30-degree north latitude, a dense zone of freshwater conflicts covering high water stress basins in South Asia, Central Asia, West Asia, and East Africa has formed. (4) International freshwater conflict has gradually become more ubiquitous, complicated and networked, and the basin communities of freshwater conflict network have increased significantly. But the "Matthew effect" of freshwater conflicts among countries are obvious, and its polarized distribution pattern is relatively stable. A "path-locking" effect has been formed among the major conflictive countries. There is a certain spatial mismatch between the quantity relationship and intensity relationship of GFCs.

[54]
Wang Y, Feng L, Liu J G et al., 2020. Changes of inundation area and water turbidity of Tonle Sap Lake: Responses to climate changes or upstream dam construction? Environmental Research Letters, 15(9): 0940a1.
Using long-term Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat satellite observations, the inundation changes of Tonle Sap Lake between 1988 and 2018 were investigated. The results show that the inundation area was stable before 2000, followed by a significant shrinking trend between 2000 and 2018. Quantitative remote sensing retrievals for concentrations of the total suspended sediments (TSS) also demonstrate an evident increasing trend (7.92 mg l−1yr−1) since 2000. A strong correlation (R2 = 0.67) was found between the annual mean inundation area and concurrent precipitation in a region located in the lower basin of the Mekong River (mostly outside the drainage basin of Tonle Sap Lake). A multiple general linear model (GLM) regression further pointed to the precipitation variation as a major contributor (76.1%) to the interannual fluctuation of the inundation area, while the dams constructed in China only contributed to 6.9%. The limited impacts of Chinese dams on the inundation area of the lake could be revealed through the limited fraction of water discharge from the Mekong River within China (∼17%). The analysis also found significant impacts of inundation changes on the recent lake turbidity increase in the dry seasons. We clearly revealed that the contribution of dam construction in China to the recent lake shrinkage was insignificant when compared with the impacts of the precipitation decrease. The results of this study provide important scientific evidence for settling water volume-related transboundary disputes regarding the control of the inundation area and water turbidity of Tonle Sap Lake.
[55]
Wen G, Wang Z, Xia S et al., 2006. Least-squares fitting of multiple M-dimensional point sets. Visual Computer, 22(6): 387-398.
[56]
Yigez B, Xiong D, Zhang B et al., 2023. Dynamics of soil loss and sediment export as affected by land use/cover change in Koshi River Basin, Nepal. Journal of Geographical Sciences, 33(6): 1287-1312.

How the dynamics in soil loss (SL) and sedimentation are affected by land use/cover change (LULCC) has long been one of the most important issues in watershed management worldwide, especially in fragile mountainous river basins. This study aimed to investigate the impact of LULCC on SL and sediment export (SE) in eastern regions of the Koshi River basin (KRB), Nepal, from 1990 to 2021. The Random Forest classifier in the Google Earth Engine platform was employed for land use/land cover (LULC) classification, and the Integrated Valuation Ecosystem Services and Trade-offs (InVEST) Sediment Delivery Ratio model was used for SL and SE modeling. The results showed that there was a pronounced increase in forest land (4.12%), grassland (2.35%), and shrubland (3.68%) at the expense of agricultural land (10.32%) in KRB over the last three decades. Thus, the mean SL and SE rates decreased by 48% and 60%, respectively, from 1990 to 2021. The conversion of farmland to vegetated lands has greatly contributed to the decrease in SL and SE rates. Furthermore, the rates of SL and SE showed considerable spatiotemporal variations under different LULC types, topographic factors (slope aspect and gradient), and sub-watersheds. The higher rates of SL and SE in the study area were observed mostly in slope gradient classes between 8° and 35° (accounting for 83%-91%) and sunny and semi-sunny slope aspects (SE, S, E, and SW) (accounting for 57%-65%). Although the general mean rate of SL presented a decreasing trend in the study area, the current mean SL rate (23.33 t ha-1 yr-1) in 2021 is still far beyond the tolerable SL rate of both the global (10 Mg ha-1 yr-1) and the Himalayan region (15 t ha-1 yr-1). Therefore, landscape restoration measures should be integrated with other watershed management strategies and upscaled to hotspot areas to regulate basin sediment flux and secure ecosystem service sustainability.

[57]
Yun X, Tang Q, Wang J et al., 2020. Impacts of climate change and reservoir operation on streamflow and flood characteristics in the Lancang-Mekong River Basin. Journal of Hydrology, 590: 125472.
[58]
Zhan X, Sohlberg R A, Townshend J R G et al., 2002. Detection of land cover changes using MODIS 250 m data. Remote Sensing of Environment, 83(1/2): 336-350.
[59]
Zhao J C, Yang K, Zhu Y H et al., 2019. The spatial-temporal changes of wetland in the Tonle Sap Lake Basin from 1988 to 2009. Journal of Southwest Forestry University, 39(6): 130-136. (in Chinese)
[60]
Zhong R D, Zhao T T G, Chen X H, 2021. Evaluating the tradeoff between hydropower benefit and ecological interest under climate change: How will the water-energy-ecosystem nexus evolve in the upper Mekong basin? Energy, 237: 121518.
[61]
Ziv G, Baran E, Nam S et al., 2012. Trading-off fish biodiversity, food security, and hydropower in the Mekong River Basin. Proceedings of the National Academy of Sciences of the United States of America, 109(15): 5609-5614.
The Mekong River Basin, site of the biggest inland fishery in the world, is undergoing massive hydropower development. Planned dams will block critical fish migration routes between the river's downstream floodplains and upstream tributaries. Here we estimate fish biomass and biodiversity losses in numerous damming scenarios using a simple ecological model of fish migration. Our framework allows detailing trade-offs between dam locations, power production, and impacts on fish resources. We find that the completion of 78 dams on tributaries, which have not previously been subject to strategic analysis, would have catastrophic impacts on fish productivity and biodiversity. Our results argue for reassessment of several dams planned, and call for a new regional agreement on tributary development of the Mekong River Basin.
PDF(5826 KB)

632

Accesses

0

Citation

Detail

段落导航
相关文章

/