Assessing Changes of Water Yield in Qinghai Lake Watershed of China

: Water yield is an important ecosystem service, which is directly related to human welfare and a ﬀ ects the sustainable development. Using the integrated valuation of environmental services and tradeo ﬀ s model (InVEST model), we simulated the dynamic change of water yield in Qinghai lake watershed, Qinghai, China, and veriﬁed the simulation results. This paper emphatically explored how precipitation change and land use / land cover change (LUCC) a ﬀ ected the change of water yield on the spatial and temporal scales. Before 2004, the areas of cultivated land and unused land showed a dramatic increasing tendency, while forestland and water area presented a decreasing trend. After 2004 cultivated land changed slowly, unused land decreased. Grassland revealed a general trend of decline during 1977–2018, while built-up land basically presented a linear increase. The results show that water yield ﬂuctuated and increased during 1977–2018. From 1977 to 2000, the mean water yield in each sub-watershed showed an increasing trend and afterward a decreasing one. After 2000, the sub-watersheds basically showed an increasing tendency. There was a strong correlation, with a correlation coe ﬃ cient of 0.954 ** (** correlation is signiﬁcant at the 0.01 level), between precipitation and water yield. Land use / land cover change can change the hydrological state of inﬁltration, evapotranspiration, and water retention. Meanwhile, the correlation between built-up land and water yield was the highest, with a correlation coe ﬃ cient of 0.932, followed by forestland, with a correlation coe ﬃ cient of 0.897. Through the analysis of di ﬀ erent scenarios, we found that compared with land use / land cover change, precipitation played a more dominant role in a ﬀ ecting water yield.


Introduction
Ecosystem services, which represent the benefits that human beings derive from ecosystem, are related to socioeconomic development and human wellbeing [1,2]. Ecosystem services include provisioning, regulatory, supporting, and cultural services [3,4]. Many ecosystem services are crucial for human survival, especially the water-related ecosystem services [5,6], such as water yield, which plays a basic role in agriculture, industry, and quality of life for humans [7,8]. Abundant water provision means the sustainable development of the regional economy and society. Quantitative evaluation and visualization of water yield is helpful to understand the trends of water supply function of ecosystem, and is beneficial to reveal the relationship between human beings and water resources, which is of great significance for scientific management and utilization of water resources [9].
Models play a significant role in assessment of water yield, because they can estimate and simulate the spatial distribution of water yield [10]. However, more and more attention has been paid to the prevent the west from desert spreading from the east, and is an important node to maintain ecological security in the northeast of Qinghai-Tibet Plateau [34,35]. In general, the temperature increases from northwest to southeast in the watershed. The watershed has a typical plateau continental climate with average annual temperatures of about −0.7 °C and a mean annual precipitation of 381 mm [36,37]. The topography of the watershed decreases from the northwest to the southeast. There are three main vegetation types in the study area, they are Potentilla fruticosa shrub, Kobresia meadow, and Achnatherum splendens steppe. The ecological environment of Qinghai lake watershed largely reflects the overall ecological change trend of Qinghai-Tibet plateau. Because of its unique geographical location and environment, the watershed has become an ecologically fragile area.

Water Yield Module
The InVEST model was developed in 2007 by Stanford University, the World Wildlife Fund for Nature and the Nature Conservancy to assess the ecosystem services and support environmental decision-making. The water yield module in the InVEST model is based on the Budyko curve and the annual average precipitation [38,39]. The model algorithm is as follows: where x Y is the water yield of pixel x , x AET is the annual actual evapotranspiration for pixel x , and x P is the annual precipitation on pixel x .
x x AET P approximates the Budyko curve [40]. 1 1 1 where x PET is the potential evapotranspiration for pixel x ,  is the empirical parameter and its expression is:

Water Yield Module
The InVEST model was developed in 2007 by Stanford University, the World Wildlife Fund for Nature and the Nature Conservancy to assess the ecosystem services and support environmental decision-making. The water yield module in the InVEST model is based on the Budyko curve and the annual average precipitation [38,39]. The model algorithm is as follows: where Y x is the water yield of pixel x, AET x is the annual actual evapotranspiration for pixel x, and P x is the annual precipitation on pixel x. AET x P x approximates the Budyko curve [40].
where PET x is the potential evapotranspiration for pixel x, ω is the empirical parameter and its expression is: where AWC x is the plant available water content and Z is empirical parameter. Potential evapotranspiration is defined as: where K c ( x ) is the plant evapotranspiration coefficient associated with LUCC on pixel, and ET 0 (x) is the reference evapotranspiration.

Data Sources and Preparation
Water yield module requires annual precipitation, average annual reference evapotranspiration, root restricting layer depth, plant available water content (PAWC), LUCC, watersheds, sub-watersheds, biophysical table, and seasonality factor (Z).
The annual precipitation and temperature data were acquired from the China Meteorological Administration. The original data were the daily meteorological observation data of meteorological station and were interpolated into 30 m × 30 m grid data by the Australian National University Splines (ANUSPLIN) package. It should be pointed out that DEM data were used as auxiliary data in the process of meteorological data interpolation. Meteorological stations within and around the Qinghai lake watershed were selected, including 16 stations, such as Tianjun station, Gangcha station, and Haiyan station.
The Penman-Monteith equation is the best equation for estimating ET 0 , because it can be used globally without any local calibration [41]. This method is limited in application by the lack of input parameters. With this in mind, Allen suggested to use the Hargreaves equation [42] for estimating ET 0 . So in this study, ET 0 is calculated based on the Modified-Hargreaves equation, its expression is: where T max and T min are the monthly highest temperature and lowest temperature respectively. R a is the monthly radiation from the sun's upper atmosphere and is obtained by Allen equation. We used the values of C, E, and T recommended by Hu on the Qing-Tibet plateau [43]. Root restricting layer depth is replaced by reference soil depth as a proxy. PAWC is defined as the difference between the fraction of volumetric field capacity and permanent wilting point. Zhou developed a method for PAWC estimation using physical and chemical properties of soil [44], its expression is: where PAWC is the plant available water fraction; sand%, silt%, clay%, OM% are the proportion of sand, silt, clay, and organic matter in soil. The soil data comes from the Harmonized World Soil Database version 1.2 (HWSD) [45]. The data set was provided by Cold and Arid Regions Sciences Data Center at Lanzhou. The LUCC data included the LUCC data in 1977,1987,2000,2005,2010,2018. The first five phases of land use data were interpreted from Landsat Thematic Mapper (TM) and/or Enhanced Thematic Mapper (ETM) images. The last phase of the natural features was interpreted by Landsat ETM images, while the artificial features were interpreted by Gaofen-2 remote sensing images. The LUCC data were classified into 22 classes at 30 m spatial resolution. Considering the difference between the data source of the first five phases and the last phase, we collected the LUCC of Qinghai lake watershed in 2010, which was derived from Data Center for Eco-Environment Protection in the Qinghai Lake Basin, and also obtained the LUCC with spatial resolution of 100 m in 2018, which was taken from the Data Center of the Chinese Academy of Sciences. Those datasets were used for cross-validation of our interpretation results.
The sub-watersheds data were generated from the 30 m digital elevation model (DEM) by Hydrology Analyst Tools in ArcGIS 10.5. Biophysical table is a table that contains the biophysical coefficients used in water yield module. Each column contains a different attribute. Lucode is the unique code corresponding to each LUCC class. LUCC_desc is to describe the name of each LUCC class. LUCC_veg contains the information on which AET Equation (1) to use. Root_depth is often given as the depth as which 95% of a vegetation type's root biomass occurs. Kc is the plant evapotranspiration coefficient for each LUCC class.
Seasonality factor (Z) is an empirical parameter which describes the rainfall intensity and topography characteristics, with typical values ranging from 1 to 30. There are three methods to calculate Z [46]: (i) Using the observed streamflow data; (ii) as a function of the average annual number of rain events; (iii) using global estimations of ω. In our study, we used the third approach to estimate Z, and its expression, that is: where P and AWC are the average annual values of precipitation and available water capacity, respectively [8]. According to this method, the calculated Z of the Qinghai lake watershed was equal to 8.5. We obtained the volume of water yield of Qinghai lake watershed over the years from the Qinghai water resources bulletin, which was published by Qinghai water resources information website.

Dynamic Change of LUCC
LUCC data were divided into 6 level (I) classes and 22 level (II) classes. The interpretation results were verified with the field survey data in 2018, indicating that the overall accuracies of level (I) class and level (II) class were 94.5% and 86.4%, respectively. We compared the results of our interpretation with the collected LUCC data. The results showed that our interpretation was 92% consistent with the published data in 2010. We used field survey data to evaluate the accuracy of published data, and we found that our interpretation accuracy was 5% higher than the data collected in 2018. Meanwhile, it is further explained that different data sources do not affect the final interpretation result.
The mainly land cover types in the Qinghai lake watershed were cultivated land, forestland, grassland, water area, built-up land, and unused land. We divided 22 secondary subcategories under the six primary land types. It can be seen from the spatial distribution of LUCC that the area of grassland coverage was the highest, accounting for approximately 63% of the total area of the watershed. That is to say, the watershed is well covered by vegetation. The area of water area, unused land, forestland, and cultivated land accounted for about 18%, 15%, 1.6%, and 1.5% of the total area, respectively. Forestland is mainly distributed in the Gangcha county. Sandy land is mostly located at the eastern and western parts of Qinghai lake. Wetland is mainly distributed in Tianjun county and Gangcha county. Built-up land is mainly concentrated in the surrounding area of Qinghai lake ( Figure 2).
The area of cultivated land increased sharply before 2000 and slowly changed after 2000. The area of cultivated land had increased by 163 km 2 in the past 40 years. From 1977 to 2000, the area of forest land decreased promptly. After 2000, it began to rise. In 2018, the area of forestland was basically the same as that in 1970s. Grassland basically showed a downward trend, which is closely related to regional human activities. Grassland area decreased by 207 km 2     The proportion of non-cultivated land conversion from cultivated land was less than 3%. Cultivated land was mostly converted into built-up land, and the proportion had been growing. Forestland was mainly converted into grassland. The percentage of forestland converted to non-forestland increased from 5.47% to 6.72% between 1987 and 2010. Grassland was mostly converted into cultivated land and unused land, and the grassland area remained above 97.76% of the original area. Water area was largely converted into unused land, and proportions of water area were 96.64%, 96.30%, 94.88%, and 95.17% in 1987, 2000, 2004, and 2010, respectively. The percentage of built-up land that had not been transformed into other types of land increased year by year and remained above 99%. Unused land was mainly converted into grassland. and finally decreased sharply ( Figure 3).
Taking 1977 as the baseline year, we calculated the land cover change transfer probability matrix of 1987, 2000, 2004, 2010, and 2018. The results indicated that the transformation rules of land types were basically the same in 1987, 2000, 2004, and 2010. The proportion of non-cultivated land conversion from cultivated land was less than 3%. Cultivated land was mostly converted into builtup land, and the proportion had been growing. Forestland was mainly converted into grassland. The percentage of forestland converted to non-forestland increased from 5.47% to 6.72% between 1987 and 2010. Grassland was mostly converted into cultivated land and unused land, and the grassland area remained above 97.76% of the original area. Water area was largely converted into unused land, and proportions of water area were 96.64%, 96.30%, 94.88%, and 95.17% in 1987, 2000, 2004, and 2010, respectively. The percentage of built-up land that had not been transformed into other types of land increased year by year and remained above 99%. Unused land was mainly converted into grassland. From the transition probability matrix of LUCC during 1977 to 2018, the percentage of cultivated land to forestland, grassland, and built-up land was 3.69%, 2.78%, and 2.81%, respectively. The proportion of forestland to grassland was 8.37%, in addition, there were 2.83% and 0.11% forestland converted to unused land and cultivated land. The proportion of non-grassland conversion from grassland was 3.62%. And 1.80% grassland was converted to unused land, it implied that grassland in these areas had degraded. Water area was basically converted into grassland and unused land, with the proportion of 2.49% and 3.00% respectively. 13.25% of the built-up land was converted to cultivated land, while 1.40% of the built-up land was transformed into grassland. It is notable that 3.54% of built-up land was converted to unused land. The percentage of unused land to grassland and water was 6.08% and 1.11% respectively (Table 1).  From the transition probability matrix of LUCC during 1977 to 2018, the percentage of cultivated land to forestland, grassland, and built-up land was 3.69%, 2.78%, and 2.81%, respectively. The proportion of forestland to grassland was 8.37%, in addition, there were 2.83% and 0.11% forestland converted to unused land and cultivated land. The proportion of non-grassland conversion from grassland was 3.62%. And 1.80% grassland was converted to unused land, it implied that grassland in these areas had degraded. Water area was basically converted into grassland and unused land, with the proportion of 2.49% and 3.00% respectively. 13.25% of the built-up land was converted to cultivated land, while 1.40% of the built-up land was transformed into grassland. It is notable that 3.54% of built-up land was converted to unused land. The percentage of unused land to grassland and water was 6.08% and 1.11% respectively (Table 1). From 1977 to 2004, the proportion of cultivated land converted to built-up land was 1.8%. An obvious phenomenon was the degradation of forestland into grassland and the transformation of water into unused land ( Table 2). From 2004 to 2018, cultivated land was transformed into forestland and grassland. Forestland was transformed into grassland and unused land. The percentage of grassland to water area was the largest. The water area was mainly transformed into grassland and unused land. The proportion of built-up land converted to cultivated land was 10.25%. Unused land was mainly transformed into grassland and water area (Table 3).

Simulation Verification
It should be pointed out that no data were published before 2000, and the water yield data in 2000 refers to the water yield in the surrounding area of Qinghai lake (QHLS). Simulation error was calculated by published data and simulation results. It turned out that simulation errors were 6.93%, 3.45% and 0.26% in 2004, 2010, and 2018, respectively. The year with the largest simulation error was 2004, and the year with the smallest was 2018 (Table 4). Therefore, it can be proved that the InVEST model can be well used to simulate the spatial-temporal change of water yield in Qinghai lake watershed. Meanwhile, it is proved that the input parameters of the model are in accordance with the geographical characteristics of the study area.

Dynamic Change of Water Yield
We focused on the temporal and spatial dynamic changes of water yield in 1977-2018. The results show that water yields tended to increase in Qinghai lake watershed. Volumes of total annual water yield, which were calculated by InVEST model during 1977 to 2018, are showed in Table 2. The lowest water yield value was in 1977, and the highest water yield value was in 2018 (Figure 4). Comparing the water yield and average annual precipitation, we found that two trends were basically the same.
Water 2020, 12, x FOR PEER REVIEW 9 of 18 lowest water yield value was in 1977, and the highest water yield value was in 2018 (Figure 4). Comparing the water yield and average annual precipitation, we found that two trends were basically the same. The watershed was delineated into 12 sub-watersheds ( Figure 5). We calculated the mean water yield in each sub-watershed. QHL-1 was the largest sub-watershed, which contributed the most to the water yield of Qinghai lake watershed. It should be noted that the mean values of sub-watershed QHL-8 and sub-watershed QHL-9 were close to 0, so no analysis was performed. The lowest mean water yield was sub-watershed QHL-2, and the highest mean water yield was sub-watershed QHL-  As can be seen from the spatial distribution, the distribution of water yield had significant spatial heterogeneity. The spatial pattern of high and low values of water yield distribution was basically consistent during 1977 to 2018. Water yield increased gradually from southeast to northwest. We can find that areas of low water yield were mainly concentrated in the surrounding area of Qinghai lake. Water yield in mountain area was obviously more than that in the surrounding area of Qinghai lake. Areas with high water yield were located in the northern, western mountainous areas and edge of the watershed, especially in the central part of Tianjun county. Regional water yields showed an The watershed was delineated into 12 sub-watersheds ( Figure 5). We calculated the mean water yield in each sub-watershed. QHL-1 was the largest sub-watershed, which contributed the most to the water yield of Qinghai lake watershed. It should be noted that the mean values of sub-watershed QHL-8 and sub-watershed QHL-9 were close to 0, so no analysis was performed. The lowest mean water yield was sub-watershed QHL-2, and the highest mean water yield was sub-watershed QHL-5. Water 2020, 12, x FOR PEER REVIEW 9 of 18 lowest water yield value was in 1977, and the highest water yield value was in 2018 (Figure 4). Comparing the water yield and average annual precipitation, we found that two trends were basically the same. The watershed was delineated into 12 sub-watersheds ( Figure 5). We calculated the mean water yield in each sub-watershed. QHL-1 was the largest sub-watershed, which contributed the most to the water yield of Qinghai lake watershed. It should be noted that the mean values of sub-watershed QHL-8 and sub-watershed QHL-9 were close to 0, so no analysis was performed. The lowest mean water yield was sub-watershed QHL-2, and the highest mean water yield was sub-watershed QHL-  As can be seen from the spatial distribution, the distribution of water yield had significant spatial heterogeneity. The spatial pattern of high and low values of water yield distribution was basically consistent during 1977 to 2018. Water yield increased gradually from southeast to northwest. We can find that areas of low water yield were mainly concentrated in the surrounding area of Qinghai lake. Water yield in mountain area was obviously more than that in the surrounding area of Qinghai lake. Areas with high water yield were located in the northern, western mountainous areas and edge of the watershed, especially in the central part of Tianjun county. Regional water yields showed an As can be seen from the spatial distribution, the distribution of water yield had significant spatial heterogeneity. The spatial pattern of high and low values of water yield distribution was basically consistent during 1977 to 2018. Water yield increased gradually from southeast to northwest. We can find that areas of low water yield were mainly concentrated in the surrounding area of Qinghai lake. Water yield in mountain area was obviously more than that in the surrounding area of Qinghai lake. Areas with high water yield were located in the northern, western mountainous areas and edge of the watershed, especially in the central part of Tianjun county. Regional water yields showed an obvious increasing trend in Tianjun County and Gangcha County between 1977 and 2018. Water yield was very high in the area of alpine tundra. Water yields fluctuated between 0 mm and 780.68 mm in study area. The maximum water yield in 2018 was abnormally higher than before ( Figure 6). In fact, the population was mainly concentrated in the surrounding area of Qinghai lake. In contrast, areas with high water yield had less population distribution.
Water 2020, 12, x FOR PEER REVIEW 10 of 18 obvious increasing trend in Tianjun County and Gangcha County between 1977 and 2018. Water yield was very high in the area of alpine tundra. Water yields fluctuated between 0 mm and 780.68 mm in study area. The maximum water yield in 2018 was abnormally higher than before ( Figure 6). In fact, the population was mainly concentrated in the surrounding area of Qinghai lake. In contrast, areas with high water yield had less population distribution.

Effect of Precipitation on Water Yield
Meteorological stations we analyzed were mainly three meteorological stations within the Qinghai lake watershed and the nearby Haiyan station. The names of the meteorological stations are Tianjun station, Gangcha station, Haiyan station, and Jiang Xigou station, respectively. We calculated

Effect of Precipitation on Water Yield
Meteorological stations we analyzed were mainly three meteorological stations within the Qinghai lake watershed and the nearby Haiyan station. The names of the meteorological stations are Tianjun station, Gangcha station, Haiyan station, and Jiang Xigou station, respectively. We calculated the cumulative annual precipitation for each station (Figure 7). We used the moving average method to fit the data of each meteorological station, in order to reflect the detailed change characteristics. Then, the linear fitting method was utilized to find a line to show the overall trend of data changes. The results showed that, the annual precipitation of Tianjun station had fluctuated from 250 mm to 450 mm, with an overall upward trend since 1955. The highest annual precipitation was 563.6 mm, and the lowest annual precipitation was 211.6 mm. Annual precipitation of Gangcha station had fluctuated from 270 mm to 450 mm, with an overall upward trend. The highest annual precipitation was 572.3 mm, and the lowest annual precipitation was 260.0 mm. It should be pointed out that Haiyan station only had precipitation data since 1976. Minimum precipitation of Haiyan station was 248.2 mm, and the maximum value was 522.3 mm. It showed a general upward tendency. Although the data of Jiang Xigou station was rare, the precipitation data of Jiang Xigou station from 1973 to 1998 only showed a slow rising trend. The highest annual precipitation was 635.7 mm, and the lowest annual precipitation was 325.2 mm.
Water 2020, 12, x FOR PEER REVIEW 11 of 18 the cumulative annual precipitation for each station (Figure 7). We used the moving average method to fit the data of each meteorological station, in order to reflect the detailed change characteristics. Then, the linear fitting method was utilized to find a line to show the overall trend of data changes. The results showed that, the annual precipitation of Tianjun station had fluctuated from 250 mm to 450 mm, with an overall upward trend since 1955. The highest annual precipitation was 563.6 mm, and the lowest annual precipitation was 211.6 mm. Annual precipitation of Gangcha station had fluctuated from 270 mm to 450 mm, with an overall upward trend. The highest annual precipitation was 572.3 mm, and the lowest annual precipitation was 260.0 mm. It should be pointed out that Haiyan station only had precipitation data since 1976. Minimum precipitation of Haiyan station was 248.2 mm, and the maximum value was 522.3 mm. It showed a general upward tendency. Although the data of Jiang Xigou station was rare, the precipitation data of Jiang Xigou station from 1973 to 1998 only showed a slow rising trend. The highest annual precipitation was 635.7 mm, and the lowest annual precipitation was 325.2 mm. We calculated the average precipitation of the three meteorological stations in 1977, 1987, 2000, 2004, 2010, and 2018. Pearson correlation analysis was subsequently performed using these values and the annual water yields of the watershed. The results of this analysis indicated that the Pearson correlation coefficient between annual water yield and annual average precipitation was 0.954 ** (** correlation is significant at the 0.01 level), with a significance of 0.001. Taking regional variability into account, we carried out the Pearson correlation analysis on the annual precipitation of each station and the water yield of the sub-watershed where the station was located. It turned out that the Pearson correlation coefficient between annual water yield and annual precipitation was 0.899 * (* correlation is significant at the 0.05 level), with a significance of 0.006 in Tianjun station. And the Pearson correlation coefficient was 0.929 **, with a significance of 0.002. When precipitation changed with time period, water yield of the watershed also changed significantly, and there was a high correlation between the two variables. Therefore, precipitation became a very important factor which affected the changes of water yield.
Geographical distribution of the ratio of water yield to precipitation, which can reflect the consistency between the two, is showed in Figure 8. Comparing the spatial distribution of water yield depth per pixel to annual precipitation, the water yield rises with high precipitation. It can be seen from the spatial and temporal dynamic changes of the ratio of water yield to precipitation that precipitation had a high contribution rate to water yield. For instance, in 2018, the highest water yield was located in the northern, western mountainous areas ( Figure 6). The ratio of water yield to precipitation also tended to be the highest in the same places, correspondingly (Figure 8). We calculated the average precipitation of the three meteorological stations in 1977, 1987, 2000, 2004, 2010, and 2018. Pearson correlation analysis was subsequently performed using these values and the annual water yields of the watershed. The results of this analysis indicated that the Pearson correlation coefficient between annual water yield and annual average precipitation was 0.954 ** (** correlation is significant at the 0.01 level), with a significance of 0.001. Taking regional variability into account, we carried out the Pearson correlation analysis on the annual precipitation of each station and the water yield of the sub-watershed where the station was located. It turned out that the Pearson correlation coefficient between annual water yield and annual precipitation was 0.899 * (* correlation is significant at the 0.05 level), with a significance of 0.006 in Tianjun station. And the Pearson correlation coefficient was 0.929 **, with a significance of 0.002. When precipitation changed with time period, water yield of the watershed also changed significantly, and there was a high correlation between the two variables. Therefore, precipitation became a very important factor which affected the changes of water yield.
Geographical distribution of the ratio of water yield to precipitation, which can reflect the consistency between the two, is showed in Figure 8. Comparing the spatial distribution of water yield depth per pixel to annual precipitation, the water yield rises with high precipitation. It can be seen from the spatial and temporal dynamic changes of the ratio of water yield to precipitation that precipitation had a high contribution rate to water yield. For instance, in 2018, the highest water yield was located in the northern, western mountainous areas ( Figure 6). The ratio of water yield to precipitation also tended to be the highest in the same places, correspondingly (Figure 8). Water 2020, 12, x FOR PEER REVIEW 12 of 18

Effect of LUCC on Water Yield
The Pearson correlation was carried out to study the influence of LUCC on the changes of water yield ( Table 5). The 22 LUCC classes were aggregated into six land cover types for each year. The Pearson correlation coefficients were calculated between changes of water yield and changes of cultivated land, forestland, grassland, water area, built-up land, and unused land. The highest positive correlation was observed between the change of water yield and the change of built-up land area, and the correlation coefficient was 0.932 *. The Pearson correlation coefficient between changes of water yield and changes of built-up land was 0.897 *. It should be noted that the changes of water yield were negatively correlated with the changes of unused land and cultivated land. In conclusion,

Effect of LUCC on Water Yield
The Pearson correlation was carried out to study the influence of LUCC on the changes of water yield ( Table 5). The 22 LUCC classes were aggregated into six land cover types for each year. The Pearson correlation coefficients were calculated between changes of water yield and changes of cultivated land, forestland, grassland, water area, built-up land, and unused land. The highest positive correlation was observed between the change of water yield and the change of built-up land area, and the correlation coefficient was 0.932 *. The Pearson correlation coefficient between changes of water yield and changes of built-up land was 0.897 *. It should be noted that the changes of water yield were negatively correlated with the changes of unused land and cultivated land. In conclusion, the change in the forestland, built-up land and unused land made major contribution to the change in the water yield.  The average water yield of unused land was basically above 115 mm, and the upper quartile was above 300 mm from 1977 to 2010, and the median was above 160 mm in 2018. The area of wetlands accounted for approximately 8.2% of the whole watershed. Wetlands which absorbed water during wet periods and released it during dry seasons could play an important role in hydrological process. In fact, water yield of wetlands was not high from 1977 to 2010, but increased in 2018. While, water yield of unused land increased from 1977 to 2018. In terms of spatial distribution, we found that the high value of water yield of unused land was mainly in the area of alpine tundra. The region of alpine tundra was less affected by human activities.
Water 2020, 12, x FOR PEER REVIEW 13 of 18 the change in the forestland, built-up land and unused land made major contribution to the change in the water yield.  The average water yield of unused land was basically above 115 mm, and the upper quartile was above 300 mm from 1977 to 2010, and the median was above 160 mm in 2018. The area of wetlands accounted for approximately 8.2% of the whole watershed. Wetlands which absorbed water during wet periods and released it during dry seasons could play an important role in hydrological process. In fact, water yield of wetlands was not high from 1977 to 2010, but increased in 2018. While, water yield of unused land increased from 1977 to 2018. In terms of spatial distribution, we found that the high value of water yield of unused land was mainly in the area of alpine tundra. The region of alpine tundra was less affected by human activities.

Mechanism for LUCC and Climate Change Influencing Water Yield
LUCC and climate change are key drivers to the change of water yield in the Qinghai lake watershed. To further explore which factor is more significant for water yield changes in the study area, two different scenarios were designed: under the scenario without climate change and under the scenario without LUCC. Under the scenario without climate change, there was no change in the input meteorological data, which was always the meteorological data from 1977. However, the LUCC data in 1977,1987,2000,2004,2010, and 2018 were treated as the input data for the corresponding years. Under the scenario without LUCC, the LUCC conditions were assumed to be the same as 1977. Finally, the two scenarios and the actual conditions were compared to reveal the impacts of LUCC and climate change on water yield. It should be acknowledged that the impact of LUCC was far smaller than that of climate change (Table 6). Under the scenario without LUCC, the water yield in the Qinghai lake watershed significantly increased, which was basically consistent with the real situation. Precipitation tends to quickly form runoff, which reduces soil moisture while increasing the water yield. Thus, precipitation directly impacts the water yield. Under the scenario without climate change, the water yield showed an increasing trend and afterward a decreasing one.

Mechanism for LUCC and Climate Change Influencing Water Yield
LUCC and climate change are key drivers to the change of water yield in the Qinghai lake watershed. To further explore which factor is more significant for water yield changes in the study area, two different scenarios were designed: under the scenario without climate change and under the scenario without LUCC. Under the scenario without climate change, there was no change in the input meteorological data, which was always the meteorological data from 1977. However, the LUCC data in 1977,1987,2000,2004,2010, and 2018 were treated as the input data for the corresponding years. Under the scenario without LUCC, the LUCC conditions were assumed to be the same as 1977. Finally, the two scenarios and the actual conditions were compared to reveal the impacts of LUCC and climate change on water yield. It should be acknowledged that the impact of LUCC was far smaller than that of climate change (Table 6). Under the scenario without LUCC, the water yield in the Qinghai lake watershed significantly increased, which was basically consistent with the real situation. Precipitation tends to quickly form runoff, which reduces soil moisture while increasing the water yield. Thus, precipitation directly impacts the water yield. Under the scenario without climate change, the water yield showed an increasing trend and afterward a decreasing one.
To explore which one had the most impact on changes of water yield, the Pearson correlation was carried out between LUCC and the change of water yield under the scenario without climate change ( Table 7). The highest positive correlation was observed between the change of water yield and the change of built-up land area, and the correlation coefficient was 0.979 **. The change of water yield was negatively correlated with the change of grassland and forestland, and the correlation coefficients were −0.952 * and −0.885 *, respectively. Correlation results indicate again that built-up land, grassland, and forestland made major contribution to the change in the water yield. Note: * correlation is significant at the 0.05 level, ** correlation is significant at the 0.01 level.

Discussion
LUCC can change the hydrological state of infiltration, evapotranspiration, and water retention [47]. The mechanism for LUCC influencing water yield was far more complex. Cultivated land had a capacity to hold the water in plants and soil. Meanwhile, crops need water to grow. During the stage from 1977 to 2000, the area of grassland decreased, which was consistent with the previous study [48]. Vegetation land produced lower water yield, because vegetated land had the ability to convert water into soil. This conclusion has been proved in existing studies [49,50]. Water yield of forestland increased, however, forestland decreased before 2010, and increased from 2010 to 2018. This was at odds with most previous studies [51], which suggested that forestland was responsible for the most water yield loss. A possible reason is that area of forestland in Qinghai lake watershed was only about 1.6% of the total area, and up to 85% of forest were sparse forest and shrubs, which have a lower water retention capacity. Built-up land is usually covered with asphalt, cement and concrete, forming an impermeable layer, which reduced the infiltration and concentration time [52]. Increases in built-up land led to increases in water yield. The conclusion was in line with the existing conclusions [53].
The finding that grassland and forestland produce lower water yield while built-up land generates more water yield was consistent with the existing studies [33,54]. But it is worth noting that the processes of LUCC are more complex. Meanwhile, conversions of LUCC lead to both negative and positive influence. Therefore, the overall consequence is not clear. In contrast, climate change is straightforward, but requires a long time to have a noticeable impact.
However, our study still has some limitations. InVEST represents bio-physical processes in a simplified manner, the model assumes all water yield from a pixel reaches a point. So the model does not distinguish between surface and subsurface water. At the same time, the limited ability of the model accounts for intra-annual variation in water supply. The retention and transport of water are both neglect in the model. Therefore, it is necessary to make some further improvement in future research. Although with these uncertainties, the results of this study can still reflect the general tendency of water yield, and reveal the relationship between water yield and LUCC and climate change. It can provide some reference information for scientific management and utilization of water resources and conserving the ecological environment.

Conclusions
Water yield is an important component of ecosystem services, which is related to regional ecological security and sustainable development of resources. In this study, we evaluated the dynamic change of water yield in Qinghai lake watershed from 1977 to 2018, and verified the simulation results with the data published by Qinghai water resources department. We also explored the influence of precipitation and LUCC on the water yield.
Our results demonstrated that water yield of the whole watershed tended to increase from 18.17 × 10 8 m 3 in 1977 to 45.45 × 10 8 m 3 in 2018. The distribution of water yield had obviously spatial heterogeneity. The water yield of mountain area was larger than that of the surrounding lake area.
The results show that precipitation change and LUCC both had important influence on water yield during 1977 to 2018. Compared with land use/land cover change, precipitation played a more dominant role in affecting water yield. In particular, we focused on the specific impact of each individual LUCC type on the change of water yield. The water yield of each LUCC type showed an increasing trend. The change in the built-up land, forestland, and unused land made major contribution to the change in the water yield.
In conclusion, InVEST model has good applicability and performance in Qinghai lake watershed. There are still some uncertainties in the simulation, but the results provide decision foundation for the scientific management of water resources and ecological environment protection.