A Study on the Spatio-Temporal Land-Use Changes and Ecological Response of the Dongting Lake Catchment

: Catchments support the survival and development of humans in a region and investigat-ing the mechanism of land-use changes and ecological responses in catchments is of great significance for improving watershed ecological service functions. Taking the Dongting Lake catchment as the study area, this study used spatial analysis, an ecosystem service value (ESV)-equivalent factor method, grid method, and other spatial analysis methods to explore land-use changes and the corresponding ecological service value response from 1990 to 2015, to provide an important theo-retical reference for ecological service management, regional planning and ecological service function improvement in the Dongting Lake Basin. Our findings are as follows: (1) apart from a trend of notable expansion in construction land, the land-use types in the Dongting Lake catchment did not change significantly. (2) Grassland had the fastest transfer-out rate; forest land were cultivated land comparably transferred to each other with a larger area, where both were simultaneously and continuously transferred out as construction land; water areas, wetlands, and construction land were all transferred in, where construction land had the fastest transfer rate. (3) The total ESV of the watershed first increased and then decreased, but the overall change was small; spatially, the wetlands and water areas had a higher ESV, whereas construction land and cultivated land had lower ESVs. (4) Soil protection, gas exchange, climate regulation, biodiversity, and water conservation are always the main ecosystem service functions of a catchment, where the service function of the catchment ecosystem is greater than the productive function. The increase in construction land was the main factor for the increased differences between the spatial distributions of the soil, raw materials, biodiversity protection, and gas exchange. on Land-Use Changes and Ecological Response


Introduction
Spatial quantification and collaborative research on land use/cover change (LUCC) are important methods for examining global changes, which have become a research hotspot in geography and ecology [1][2][3]. LUCC affects the structure, quality, process and other aspects of ecosystems and promotes changes in ecosystem services [4][5][6][7][8]. Ecosystem

Study Area
The Dongting Lake Basin is located in the southern part of central China (ranging from 24°38′-30°26′ N and 107°13′-114°18′ E) and includes the Hunan, Hubei, Guangxi, Guizhou, and Chongqing provinces ( Figure 1). It covers a total area of approximately 26.28 km 2 , of which Hunan Province accounts for 78%, Guizhou Province accounts for 11.6%, and Guangxi, Chongqing, and Hubei provinces together account for 10.4%. The basin is surrounded by mountains to the east, south, and west; the northern part of the basin is the Dongting Lake Plains; and the central area consists of hills and basin, forming a unique "horseshoe-shaped" pattern. The basin has a typical subtropical monsoon climate zone with alternating winter and summer monsoons. The precipitation in the basin is unevenly distributed with large inter-annual variability and a high frequency of floods and droughts. The water system in the basin is complex, including four large rivers: the Xiangjiang, Zijiang, Yuanjiang, and Lishui rivers. Dongting Lake is the second largest lake in China and plays an important role in ecosystem protection of the entire Yangtze River Basin.

Data Resources and Processing
The data used in this study mainly include three types of data, including basic geographic data, meteorological data, and socio-economic data.
(1) The basic geographic data included administrative division data, land-use data, and digital elevation models (DEM), which were mainly obtained from the Resource and Environment Science and Data Center of the Chinese Academy of Sciences (https://www.resdc.cn/, accessed on 27 July 2021). The administrative division data were the 2015 national, provincial, and municipal vector boundaries. The 1:100,000 land-use data were obtained from the interpretation of Landsat series images and included six years: 1990, 1995, 2000, 2005, 2010, and 2015. The land-use data were divided into two levels. The first level was divided into six categories: cultivated land, forest land, grassland, water area, construction land, and unused land; the second level was divided into 26 categories [29]. This set of data has been extensively applied in relevant studies with a higher classification accuracy [30][31][32][33]. To study land-use in the Dongting Lake Basin, we adjusted for the first-level classification structure that merges the second level tidal flats, beaches, and marshes into wetlands, together with the other six first-level categories. Referring to the relevant literatures [34][35][36], compared with other DEM data, the application effect of SRTM DEM was relatively good; thus, the 30 m resolution DEM of SRTM 3 was mosaicked and clipped to provide the DEM data for the Dongting Lake Basin area. (2) Meteorological data, including rainfall and temperature, were obtained from the National Meteorological Information Centre of China (http://www.nmic.cn/, accessed on 27 July 2021). (3) The socio-economic data included 10 indicators: GDP; population; output values of the primary, secondary, and tertiary industries; fixed asset input; per capita income; agricultural product output; retail sales of social goods; and fiscal expenditure. These were all obtained from the Hunan Province Statistical Yearbook.

Land-Use Change
The study aimed at quantifying the ESV of Dongting Lake Basin and to explore the mechanism of its response to LUCC. Thus, we first analyzed the land use changes, and then adopted the 1 km grid method to spatially quantify the ESV of the study area, and associate its spatiotemporal change process with the LUCC changes. The methods used in the study mainly included land cover dynamic index, conversion matrix, ESV-equivalent factor method, etc., which almost were implemented in ArcGIS10.4 software.

Land Cover Dynamic Index
The land cover dynamic index ( ) and bilateral LCDI ( ) are good indicators of the intensity of land-use change, in which LCDI can characterize the direction and speed of land-use change effectively and BLCDI can quantize the comprehensive speed of landuse change well [37]. These are calculated as follows: where is the area of land cover type at the beginning of the study period in the study area, is the area of type at the end of the period in the study area, ∑ is the total area of type converted into other types, ∑ is the total area of type converted from other types, is the study period, and and are the dynamic index and bilateral dynamic index of land cover type , respectively [37,38].

Conversion Matrix
A conversion matrix is the main quantitative approach for examining the number and direction of different types of converted land-use. It can specifically reflect the structural characteristics of land-use change and the corresponding directions of conversion. Table 1 expresses the mathematical formulas for the conversion matrix.
In Table 1, and represent the land-use types at the beginning and end of the study period, respectively, represents the unchanged area of land type during the study period, represents the area of land-use type at the beginning of the study period, represents land-use type at the end of the study period, − represents the reduction in the area of land type during the study period, and − represents the increase in the area of land type during the study period.   [23,40]. We used the ESV equivalent per unit area of the Chinese ecosystem as estimated by Xie et al. [23] to evaluate the ecosystem services of the study area. To reflect the regional differences, we calculated the ESV equivalent as follows [41,42]: where is the value of the ESV-equivalent factor (RMB·km −2 ·a −1 ), P is the national average food price (RMB·kg −1 ), Q is the unit grain output in the study area (kg·km −2 ), and n is the number of years. To facilitate comparison, the national grain unit price in 2015 was selected as 2.36 RMB·kg −1 . The unit grain output in Dongting Lake Basin was calculated as 527,156 kg·km −2 by multiplying the average annual grain yield and the area proportion of Hunan, Guizhou, Hubei, Guangxi, and Chongqing. After regional correction, the value of the equivalent factor was 303,346 RMB·km −2 ·a −1 . The equivalent coefficients are shown in Table 2. The value of each ecological service function of the construction land was set to 0. The ESV was calculated as follows: = ∑ ( × ) (4) = ∑( × ) (5) where the ESV is the ecosystem service value, is the area of land-use type k (km 2 ), is the ecological value coefficient (RMB·km −2 ·a −1 ), is the value of the f-th ecosystem service function, and is the f-th ecological value coefficient of land-use type k (RMB·km −2 ·a −1 ). Owing to the advantages of using the 1 km grid [43], the above ESV calculations were combined with a 1 km grid for spatialization. Finally, the 1 km grid data for the ESVs of Dongting Lake Basin were obtained, which was used to analyze the spatiotemporal changes of ESVs and to compare with an analysis results of land-use changes through spatial analysis function of ArcGIS 10.4 to explore the driving mechanism of ESV changes.

Principal Component Analysis
Principal component analysis (PCA) is a class of mathematical statistical means to remove data redundancy by transforming the relevant multivariate spatial data into a few unrelated comprehensive indicators through the rotation of the original spatial axis to simplify the data structure [44]. Based on the principle of variance maximization, PCA can use as few principal components as possible to replace the original index while retaining as much information of it as possible, simplifying the complex problem. Through PCA, the information is concentrated into the first few principal components, usually with a cumulative contribution rate above 95%, which can be used as representatives to reflect the original information. In this study, the PCA is used to explore and identify the main drive factors to the spatial differences of ESVs with IBM SPSS 23.0.

Spatio-Temporal Analysis
Damages to the ecosystem functions of land caused by LUCC is the main reason for regional ecological and environmental effects [45,46]. The results were divided into timeseries and spatial-series analyses. The time-series analysis (Table 3) showed a decrease in the area of arable land, forest land, and grassland in the basin. Among these, the area of arable land decreased the most (0.88%). The areas of construction land, water area, and wetlands increased, where construction land had the largest increase (1.23%). As the proportion of unused land was only 0.1% of the area, it was calculated without a specific analysis. Based on the spatial-series analysis ( Figure 2), forest land was found to be the main land type in Dongting Lake Basin, with an area proportion of more than 60%. Forest land was widely distributed and concentrated in the hilly and mountainous areas to the west and south of the study area. The distribution of arable land was relatively concentrated, mainly in the hinterland plains in the middle of the study area and the plains near Dongting Lake in the northeast. Grassland was mainly distributed in the mountainous areas of the west, northwest, and southwest. The vast and mixed water system was distributed throughout the basin; the largest water body was Dongting Lake in the northeast. Compared with the other land types, construction land was relatively scattered; the trend of outward expansion was relatively notable. Large patches were mainly concentrated in the eastern and central areas of the basin.

Land-Use Transfer Matrix Analysis
The land-use transfer matrix can quantitatively analyze important information such as the direction and speed of changes in each land-use type, which is conducive to further understanding the mechanism of land-use change. Combined with the land cover dynamic index (Formulas (4) and (5)), the land-use transfer matrix was calculated based on land-use data from 1990 to 2015 (Table 4, Figure 3), which characterizes the changed laws of land-use more effectively. The results showed that arable land was mainly transferred out, with the largest area being transferred to construction land, followed by forest land, and compared to spatial distribution, the reduced area of arable is closed to urban areas, especially Changsha, Zhuzhou and Xiangtan ( Figure 3a). Forest land was also mainly transferred out with an LCDI degree of −0.01%; the area converted to construction land was also the largest, followed by that converted to arable land; in spatial, as shown in Figure 3b, the eastern cities, such as Yiyang, Yueyang, Changsha, Zhuzhou, Xiangtan, Loudi, Shaoyang, and Hengyang, have many reduced areas of forest land, and many western cities, like Zhangjiajie, Xiangxi, Huaihua and so on, have many increased areas, so the differences between east and west changes in forest land are more obvious. In a word, during the study period, the forest and arable land were transformed into each other; meanwhile, the transfer to construction land was always prominent, mainly in the suburbs and adjoining areas. The grassland was mainly transferred out with an LCDI degree of −0.28%, indicating the fastest decrease, and the areas with a significant decrease were mainly in the western regions (Figure 3c), such as Huaihua and Zhangjiajie. The water area and wetland were mainly transformed into each other; part of the area was converted to arable land, especially in the Dongting Lake area, which is closely related to agricultural economic development (Figure 3d,e). The trend of external expansion of construction land was obvious based on the central urban area of each city (Figure 3f), and it was mainly transferred in with an LCDI degree of 4.28%, having the fastest increased speed with an BLCDI degree of 4.79%. The increase in construction land was significantly greater in the eastern part of the study area than in the western part. The reasons for this include the influence of the topography-high west and low east-and the differences in economic development between the east and west. Overall, the forest land and arable land transform to each other comparably from 1990 to 2015 and both also transferred out as construction land with a larger area. Grassland was transferred out at the highest rate, and water areas and wetlands were transformed into each other. The expansion of construction land was obvious and had the fastest conversion speed.
Note: LCDI, land cover dynamic index; BLCDI, bilateral land cover dynamic index.

Analysis of ESV of the Dongting Lake Basin
Although the overall land-use changes in Dongting Lake Basin were not significant during the time-series, the spatial differences were substantial. It is worth noting that land-use changes can only indirectly reflect certain mechanisms of change in the ecological environment of a basin; the ESV best reflects the final response result. The ESV is one of the most important means to quantitatively analyze the service functions of regional ecosystems and is a relatively mature technique. Thus, this study evaluated the changes in the ecosystem functions of the basin at two levels: the overall ESV and individual ESV.

Time Change Characteristics of the Total ESV
We combined Formula (4) with land-use data to calculate the ESV of Dongting Lake Basin ( Table 5). The results showed that the total ESV of Dongting Lake Basin presented an initial increase followed by an overall decrease, in which the ESV reached its maximum in 2010 while the minimum was in 2015 and the changing amplitude of ESV between 2010 and 2015 was the largest. Compared with the ESV of various land-use types, the total ESV of arable land and forest land accounted for approximately 90% of that of the entire basin. Thus, these were the main bodies supporting the ecosystem services of Dongting Lake Basin. During the study period, the ESV of arable land, forest land, and grassland decreased to varying degrees, which was consistent with the corresponding changes in the land-use area; the ESV of grassland had the largest reduction. The water areas and wetlands accounted for only 3.5% of the study area, but their ESVs accounted for over 10% of the total, and compensated for the loss in the ESV caused by reductions in the arable land, forest land, and grassland. Since the implementation of the Program to Return Farmland to Forests in 2000, the ESV of arable land has decreased and the ESV of woodland has partially increased to a certain degree. Moreover, except for grassland, the ESVs of the water areas and wetlands have also significantly increase, especially wetlands, which had the largest increase (17.42%). Overall, the ESV of the basin decreased; in the other words, in terms of the time-series, land-use changes in the basin were not conducive to the maintenance of the ecological service function and ESV.   (4), a 1 km grid was used to spatialize the ESV of the Dongting Lake Basin, and the ESV of each period from 1990 to 2015 was averaged to obtain the average distribution across the study area. Combined with the natural break point method, the ESV of the Dongting Lake Basin was divided into five levels: the low value area (0-3.72 × 10 6 RMB·km −2 ), lower value area (3.72 × 10 6 -5.58 × 10 6 RMB·km −2 ), medium value area (5.58 × 10 6 -7.73 × 10 6 RMB·km −2 ), higher value area (7.73 × 10 6 -11.28 × 10 6 RMB·km −2 ), and high value area (11.28 × 10 6 -14.38 × 10 6 RMB·km −2 ) (Table 2, Figure 4). The contribution of different land-use types to the regional ESVs varied significantly. The ESVequivalent coefficients of the water areas and wetlands were relatively large and the ESVs were notably higher, thus indicating that these are high-value areas. The ESV-equivalent coefficient of construction land was small and the ESV was low, indicating a low-value area. The ESV-equivalent coefficient of arable land was slightly higher than that of construction land, indicating a lower value area. The ESV-equivalent coefficient of forest land and grassland was intermediate between that of water and arable land, and the coverage was wide, indicating a generally higher value area. Overall, the areas to the northeast of the basin near Dongting Lake, southwestern areas, and central urban areas of various cities had low ESVs. The mountainous areas to the west of the middle of the basin, the northwestern and southeastern margins, and the lake areas such as Dongting Lake had relatively high ESVs. The spatial difference of total ESV may be closely related to the topography and regional economic development level. In general, farmland and construction land are located in plain areas, which often have lower ESVs due to smaller ESV-equivalent coefficients. At the same time, forest land and grassland are often in mountainous areas, such that these areas generate higher ESVs with larger ESV-equivalent coefficients. Moreover, the expansion of construction land has direct links to the regional economic development level, i.e., the more developed the regional economy, the lower the regional ESV. In the study area, the development level of cities surrounding Dongting Lake is better than other cities, especially the western cities, and the ESVs of these areas are relatively low. Thus, we need to consider more the ESV promotion in these developed cities or region.

Temporal Characteristics of Individual Ecological Services
According to Formula (5), the value of individual ecological services can be calculated, as listed in Table 6. From the perspective of the value of an individual ecological service, soil formation and protection, water conservation, gas exchange, climate regulation, and biodiversity protection were the main ESVs of the basin, where the sum of the proportions exceeded 70% of the total value. The contributions of food production, entertainment, and leisure were small, whose sum was less than 10%. From 1990 to 2015, the value of the ecological services increased in functions such as climate regulation, water conservation, waste treatment, and entertainment, among which the greatest increase was in the water conservation function (0.16%), followed by waste treatment (0.12%). Since the large flood in 1998, the ecological function of Dongting Lake has received increasing attention. The protection of water resources in the lake area has continuously increased, and various lake and wetland resources have been repaired and protected by policies, especially the implementation of the Program to Return Farmland to Forests in 2000 and Rural Environmental Remediation, which effectively promoted the function of water conservation and waste treatment. The ESVs decreased mainly in functions such as gas exchange, soil formation and protection, biodiversity protection, food production, and raw materials. Among these, the greatest decrease was in soil formation and protection (−0.15%), followed by gas exchange (−0.07%). This was mainly because the area of arable land and grassland-both with higher contributions to the two functions-had been significantly reduced, whereas construction land, with a lower contribution to the functions, had increased rapidly and could not effectively supplement the loss of its functional value. Moreover, the majority of the Dongting Lake Basin belong to the red soil region, and in recent years, frequent extreme climate events and an increased rainfall intensity have also reduced the soil protection function of soil formation to some extent. In a word, the service function of the basin ecosystem was higher than the productive function. At the same time, the value of the productive function showed a downward trend. The value of each ecological service was significantly different from the spatial distribution ( Figure 5). For gas exchange (Figure 5a), the high-value region of its ESV was mainly located in mountainous areas and closely associated with woodlands. For climatic regulation and headwaters conservation (Figure 5b,c), the high-value region was in Dongting Lake and its hydrographic net, respectively, which are affected by the water area. For soil formation and protection (Figure 5d), most of the high-value regions were located in woodlands. In waste disposal (Figure 5e), the high-value regions were also influenced by water bodies. For biodiversity protection (Figure 5f), the high-value regions were mainly distributed in woodlands and at the Dongting Lake, while the low-value regions were located on construction land and farmland. Farmland had the most significant effect on food production within the low-value regions located (Figure 5g). The ESVs of raw material were also dominated by woodlands (Figure 5h). Entertainment (Figure 5i), which mainly refers to providing opportunities for recreational, cultural, and artistic values, was affected by water areas and wetlands with high-value regions in the Dongting Lake region.
In general, water areas contributed the most to the ecological service function of water conservation, followed by wetlands. Simultaneously, water areas and wetlands also significantly contributed to waste treatment, climate regulation, and entertainment, ultimately leading to a correspondingly higher ESV with significant differences in the spatial distribution. Arable land contributed most to the food production ecological service function, and the unit value was considerably higher than that of the other land-use types. The ESVs of raw materials, biodiversity protection, and gas exchange were generally higher in forest land and were significantly lower in arable land and construction land. Therefore, differences in the spatial distribution of the value of individual ecological services were markedly related to the land-use type. The increase in construction land was the main factor associated with the increase in the spatial distribution of the ESVs, such as soil formation, raw materials, biodiversity protection, and gas exchange.

The Exploration of Drive Factors
Although the land-use patterns and changes have a decisive driving process for the ESVs in the study area, in a fundamental sense, the spatial difference of ESVs was the result of a comprehensive effect of natural and human activities. Therefore, based on data accessibility and spatialization, this study attempted to use the principal component analysis (PCA) method and to select relevant drivers from the scope of nature and social economy to further explore the drivers and mechanisms that caused the spatial differences in the ESVs.
We divided the factors into two categories: natural factors and socio-economic factors. Based on the typicality and comprehensiveness of the factors, the natural factors were the DEM, rainfall, and temperature, whereas the socio-economic factors were the population, GDP, and the output values of primary and secondary industries ( Table 7). The average value of each factor for a total of six periods from 1990 to 2015 was calculated using the county as a unit. The spatial differences in the ESVs driving the natural and socioeconomic factors were explored using the principal component factor method and regression analysis. The multiple regression model was as follows: = −2.057 × 10 − 0.416 − 0.168 + 0.101 ( 2 = 0.211, P < 0.05) (6) where represents the ESV after standardization, and , , and represent the three principal component factors extracted by PCA after the standardization of each factor. The cumulative contribution rate was 89.43%. Combined, the regression model and PCA (Table 7) showed that principal component F1 reflected the GDP, population, secondary and tertiary industries, per capita income, retail of social goods, and fiscal expenditure, among others; in other words, socio-economic information was more prominent. The coefficient of was negative, indicating that the social factors had a negative correlation with the ESV. Principal component reflected information regarding the output of the primary industry and agricultural products, characterized by a negative correlation with the ESV. Principal component F3 reflected information regarding natural factors, such as the temperature, DEM, and rainfall. The temperature and rainfall had a negative correlation with the ESV, whereas the DEM had a positive correlation. In conclusion, the socio-economic factors exerted more effects on the ESVs of Dongting Lake Basin than the natural factors.

Discussion
The Dongting Lake basin covers a wide area, and its ecological service function is of great significance to the whole middle and upper reaches of the Yangtze River. This study used a 1 km grid to conduct spatial quantitative research on the ESV of the Dongting Lake basin, which was conducive to explore the fundamental driving forces and to identify the weak area of ecological service function in the Dongting Lake Basin, having important reference significance for realizing watershed management and sustainable development.

The Determination of Evaluation Method and the Comparison of the Results
In fact, it is difficult to unify the results with so many evaluation methods of ESVs. Currently, the spatialized estimation methods of ESVs can be roughly divided into two categories: a "material quantitative method" and "equivalent factor method". The material quantitative method is that the total value is derived from the material amount of ecosystem services and the unit price based on it. Through this method, a simulation of ecosystem service functions in study areas was performed by establishing production equations between a single service function and local eco-environmental variables. The method requires more parameters with a complex computational process, and the evaluation method and parameter criteria for each type of ESV are also difficult to unify. For example, Kong et al. (2019) evaluated five ecological service functions in the Dongting Lake catchment, namely water conservation, soil conservation, carbon fixation, flood regulation and biodiversity protection, and put forward suggestions for protection and restoration [25]. Ma et al. (2019) estimated the material amount of water yield, soil conservation, carbon storage and schistosomiasis prevention of Dongting Lake wetland, and valued them through the unit price [47]. Liu et al. (2019) evaluated the water conservation function in the upper reaches of the Minjiang River, and estimated its ecological service value by a market value method [48]. The equivalent factor method is a method to construct the value equivalent of various ecosystem service functions based on distinguishing different types of ecosystem service functions, and then evaluates them combined with the distribution area of each ecosystem. Compared with the material quantitative method, the equivalent factor method is applied widely due to the intuitive and small amount of the data required, and it is particularly suitable for the ESV evaluation at the regional and global scale [39,40]. Therefore, the equivalent factor method was chosen to evaluate the ESVs in this study. Obviously, even in the same region, the results are diverse and not comparable due to different methods. While comparing them with the results of the regions contained within the watershed by the same method, such as the evaluation results of ESVs in Hunan Province conducted by Xiong et al. (2018) [49] and the results of the Dongting Lake wetland studied by Deng et al. (2019) [50], the results of this study were small and the main reason for this was because the study revised per unit area yield with the area occupied by multiple provinces, which were relatively small, so the ESV of the whole basin is relatively low. Still, the changed laws of ESVs in the corresponding region are consistent, which was conducive to grasp the changed laws of ecological service function in local regions or to identify the weakened areas.

The ESV at the County Scale
In order to further explore the changed laws of the ESV in local regions, this study calculated the ESV in each county by spatial statistical analysis of ArcGIS 10.4, as shown in Figure 6. At the county scale, the counties within the Dongting Lake had the highest ESV, owing to influences from lakes and wetlands. Among the four major sub-basins, the ESVs of the counties adjacent to the Xiangjiang River Basin were relatively low, whereas those of counties adjacent to the Yuanjiang River were significantly higher. There were two reasons for this dichotomy. First, the topography of the Xiangjiang River Basin is lower than that of the Yuanjiang River Basin. Second, the Xiangjiang River flows through the central urban areas of many cities, such as Changsha and Zhuzhou, among others. These areas are economically developed, densely populated, and subject to a considerable interference from human activity. Therefore, we must focus on the protection and restoration of the ecological system in and around the Xiangjiang River Basin, especially in cities located in the upstream area. We must also strengthen the protection and restoration of high-ESV areas, such as wetlands, water areas, and forest lands, maintain and strengthen the continuity and integrity of ecosystem patterns, and guide the development of land-use in the direction of ESV preservation or appreciation.

Limitations of the Study
In this study, the factors selection is not enough and perfect due to the fact it only considers the availability of data. Thereinto, natural factors, such as evaporation and wind speed, need to be added in and social-economic factors can also include policy factors such as returning the farmland to forest land. Meanwhile, although the principal component analysis was used to preliminarily explore the driving factors of the ESV in the study area, the analysis of the internal driving mechanism process was far from enough, and at least the in-depth analysis of the social-economic factors was required, which also was closely related to the availability of other types of data. Moreover, the R 2 was not sufficiently large enough, at only 0.211, possibly because the Dongting Lake Basin is too large. Thus, further subregional analyses on the Dongting Lake Basin should be performed in future research.

Conclusions
The Dongting Lake Basin is one of the most important sub-basins in the middle reaches of the Yangtze River. The ecological services of water conservation, biodiversity protection, food production, and other functions are highly important for regional development. In this study, we used land-use data from 1990 to 2015 and combined it with an ESV-equivalent exchange algorithm to analyze the temporal and spatial changes in the land use and ESVs in the Dongting Lake Basin. The results of this study can be a great reference point and could guide significance for watershed management, regional planning, and river basin ecological service function improvement of Dongting Lake Basin. The main conclusions are as follows: (1) Apart from construction land, the land-use changes in the Dongting Lake Basin did not show significant changes in other land-use types. The main land-use type in the basin was forest land, with an area proportion of ≥ 60%. Construction land had the largest increase at 1.23%, with a notable trend of outward expansion. The area of arable land decreased the most, at only 0.88%. (2) During the study period, the conversions of various land types were markedly different. Forest and arable land were transformed into each other with comparable area proportions. Construction land occupied most of the share of forest and arable lands that were transferred out, resulting in a change in these two land-use types. A large area of grassland was transformed to forest and arable land, which was characterized by a rapid conversion speed. Water areas and wetlands were mainly transformed in. Construction land was mainly transformed in with the fastest conversion speed. (3) The overall ESV of the basin first increased and then decreased. However, the magnitude of the change was small. Forest and arable lands have continuously supported the ecosystem services of the Dongting Lake Basin, but there was a small reduction (to varying degrees) during the study period. The ESVs of the water areas and wetlands increased significantly, partially compensating for the overall ESV loss. In terms of the spatial distribution, the ESV was higher in areas adjacent to wetlands and water areas and were lower in areas adjacent to construction and arable lands. The ESV differentiation was notable. (4) The main ecological service functions included soil formation and protection, water conservation, gas exchange, climate regulation, and biodiversity protection. The value of the ecological service functions was generally greater than that of the productive functions. The differences in the distribution of various ESVs were closely related to the land-use types. The increase in construction land was the main factor driving the increase in the spatial distribution of the ecological service functions, such as soil formation, raw materials, biodiversity protection, and gas exchange.
In this study, land-use data based on a Landsat data source has great advantages in spatial resolution, but the time resolution is poor. In the future, it is necessary to monitor the ecological service value of the basin from a higher temporal resolution, to dynamically grasp the ecological environmental changes of the Dongting Lake basin.
Author Contributions: Nan Yang devised the idea and wrote the paper; Wenbo Mo performed the analysis; Maohuang Li, Xian Zhang and Min Chen conducted the data processing; Feng Li and Wanchao Gao helped in data collection. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available from the first author (Y.N.) on reasonable request.