1. Introduction
Lakes are an important part of the earth’s hydrosphere and are a link that connects all spheres of the earth [
1]. They play an important role in the surface water cycle [
2]. In recent decades, China has been attaching great importance to the study of lake resources. Since the 1960s, China has conducted two national lake surveys and many local lake surveys [
3,
4,
5,
6,
7]. This had led to the achievement of milestones of research progress [
8,
9]. These efforts have made important contributions to the construction and development of Chinese lake databases and provide valuable information to comprehensively understand the dynamic changes of Chinese lakes and their response to climatic change and human activities.
Lakes are sensitive indicators of climatic change [
10,
11]. As a sensitive area of climate change and a vulnerable ecosystem area, the dynamics of lakes on the Qinghai-Tibet Plateau has become an important indicator for studying global climate change. Existing studies show that lakes on the Qinghai-Tibet Plateau are undergoing rapid expansion and that an increase in precipitation is a major reason for the expansion of lakes, followed by the melting of glaciers and the degradation of frozen soil [
12]. The Qaidam Basin is located northeast of the Qinghai-Tibet Plateau and is a typical arid inland basin in China. In contrast to other regions of the Qinghai-Tibet Plateau, this region has a dry climate and limited precipitation [
13]. Since lakes are among the scarce surface water resources, they play an important role in regulating the climate of arid regions and maintaining local economic development [
14]. Studies have shown that the air temperature in the Qaidam Basin changed by 0.53 °C/10a from 1961 to 2010 [
15], which is much higher than the trend of changes in other parts of the Qinghai-Tibet Plateau. Thus, the Qaidam Basin has become a key area for the study of climatic change.
There are various types of lakes in the Qaidam Basin, and most of them are terminal, salt water, and salt lakes. The change in terminal lakes is the result of the water volume balance in the basin and is very sensitive to the climatic change in the basin and even the world [
16]. Salt lakes are important economic resources in the basin. The unique salt lake industry in the basin is the result of the most densely populated area of salt lakes in the world [
17]. Thus, the lakes in the basin are an important foundation for local economic development. However, the development of salt lakes is often accompanied by intensive human activities [
18]. Unreasonable development could not only affect the sustainable use of resources but can also cause irreversible damage to the regional microclimate. In view of the dual characteristics of lakes as both ecological and economic resources in the basin, it is necessary to comprehensively study the changes of lakes in the basin and their responses to climatic change and human activities.
Existing studies have shown that lakes in the Qaidam Basin respond significantly to changes in climate and land cover and human activities. The climatic factors that affect the change in lakes are regional as well, and this has been shown by lake precipitation, air temperature, and evaporation [
19]. Analyses of changes in the land cover have proven that changes in vegetation and the cropland area are also important factors that affect the changes in lakes in the basin [
17]. Single factor correlations have been studied by singly exploring the relationship between each factor and lake change [
17,
19]. However, lake change is a comprehensive result of the influence of multiple factors. These studies cannot compare the contribution of various factors to lake change and are not conducive to forming an in-depth understanding of the driving forces in lake change. A multiple linear regression (MLR) method can evaluate the impact of different independent variables on dependent variables by establishing the relationship between dependent variables and multiple independent variables. Thus, these types of analyses are applicable to analyze the driving forces involved in lake changes. Currently, this method is widely applied in the fields of economics and sociology, but it is rarely considered in analyzing the driving forces of lake changes [
20,
21]. Human activities that impact changes in lakes in the basin include the exploitation of salt lakes and the construction of dams [
18]. However, owing to the difficulty of quantifying human activities, it is problematic for current studies to directly analyze their correlation with lake changes. The lake area is currently the main dependent variable in the analysis of drivers of lake change, but due to differences in topography around the lakes, the lake area does not fully reflect the true variability of lake water [
22]. In addition, current studies in the Qaidam Basin are primarily based on analyses of individual lake changes, and there is a lack of understanding of the lake changes on the whole basin scale [
17]. Some studies have conducted analyses on the scale of the whole basin [
18]. However, they still adopted the lake area as the indicator, and this approach cannot adequately analyze the differences in changes across subregions.
The objective of this study was to analyze the lake water volume changes and their driving forces in the Qaidam Basin by combining multi-source remote sensing, model simulations and data from historical sources. First, the spatiotemporal pattern of climatic change in the basin under the background of global warming was analyzed with lake water volume as the main indicator. Changes in the lake water volume in different regions of the basin in response to climate change and human activities are then comprehensively discussed. Finally, the MLR model was used to delineate the main factors that affect the change in lake water volume in different regions. This study reveals the characteristics of changes in the water volume and driving forces of the lakes in inland arid basins in China and is intended to provide a theoretical basis for managing climatic change and the rational utilization of lake resources in inland arid regions in China.
2. Materials and Methods
2.1. Overview of the Study Area
The Qaidam Basin (34°37′57”~39°18′42”, 90°27′56”~99°25′51”) is an enormous tectonic subsidence basin northeast of the Qinghai-Tibet Plateau that is approximately 800 km long from east to west and 300 km wide from north to south with a total area of 276,000 km
2. The basin connects the Tibet Autonomous Region in the south, reaches Gansu Province in the north, borders the Xinjiang Uygur Autonomous Region in the west, and borders the Haibei and Hainan Tibetan Autonomous Prefectures of Qinghai Province in the east. It is the central hub among Qinghai, Gansu, Xinjiang, and Tibet. The basin is surrounded by high mountains, including the Kunlun Mountains in the south, the Qilian Mountains in the northeast, and the Altun Mountains and Qimantag Mountains in the northwest. The highest altitude in the basin is 6826 m, while the lowest is 2676~3200 m. As one of the most famous inland arid basins in the world, the Qaidam Basin has a dry climate and rare precipitation. Many lakes, salt lakes, and marshes are also characteristic of the Qaidam Basin. There are 31 lakes larger than 10 km
2 in the basin, and most of them are terminal lakes in the center of the basin. These lakes only have water inflow and no water outflow. Moreover, owing to powerful evaporation, mineralization is very high, and most lakes are saltwater or saline. Alternatively, inflow-outflow lakes contain a large amount of freshwater resources and are primarily distributed in the mountainous areas with an altitude of more than 4000 m north of the Kunlun Mountains. Among them, the Donggi Conag Lake is the largest alternately inflow-outflow lake in the basin. According to the water resources zoning scheme of the first and second national comprehensive water resources assessment of China, the Qaidam Basin was divided into 16 level 4 water resources zones (WRZs), and we analyzed the spatial distribution of lakes by analyzing the area of lakes within different WRZs [
23]. The names and ID of WRZ and lakes in the Qaidam are provided in
Table 1.
2.2. Data
Three types of publicly available data were examined in this study, including lake storage variables, meteorological elements, and land cover (
Table 2).
There are two types of data on lake storage variables that were used in this study. Among them, QB_LWSC (Qaidam Basin Lake Water Storage Change) is a lake storage variable that was estimated using multi-source remote sensing images and satellite altimetry data combined with a pyramid model (
https://github.com/zhangch-ss/Monitoring_of_lake_dynamics, accessed on 13 May 2022) [
24]. This dataset has a time span of 1990–2020, and a time resolution of 1a that is used to analyze changes in the time series and driving forces. TPDC_LWSC (Tibet Plateau Science Data Center Lake Water Storage Change) was released by the Qinghai-Tibet Plateau Science Data Center (
https://data.tpdc.ac.cn/zh-hans/, accessed on 13 May 2022) and formulated using remote sensing images, Shuttle Radar Topography Mission (SRTM), and the pyramid model [
25]. This dataset spans from 1990 to 2019 and has a time resolution of 5a. However, the small amount of data available is unsuitable for a correlation analysis. Therefore, it was primarily used to compare and verify the QB_LWSC in this study. The lakes used to compare and verify the TPDC_LWSC and QB_LWSC were lakes 1, 2, 3, 5, 7, 9, 16, 21, 25 and 31 (
Figure 1).
Meteorological elements include precipitation, air temperature, and evaporation, which originated from MSWEP and GLDAS-Noah, respectively. MSWEP is a set of multi-source integrated precipitation products formed by integrating the advantages of stations, remote sensing, and reanalysis data (
http://www.gloh2o.org/mswep/, accessed on 20 May 2022) with a spatial resolution of 0.1°. This data is well suited to evaluate the lakes in inland China and is highly applicable for areas that lack data [
26]. GLDAS-Noah is the global land assimilation data generated using the Noah model based on satellite and ground observation data (
https://disc.gsfc.nasa.gov/datasets, accessed on 20 May 2022) and has been proven to perform well in the arid regions of northwest China [
27,
28]. The evaporation and air temperature data provided in GLDAS-Noah were the only values used in this study. However, the V2 version only provides data since 2000, while the V1 version provides data from 1984 to 2014. Therefore, the two sets of data of the overlapping period (2000–2014) were subjected to deviation correction for this study, so that they met the analytical requirements for a long time series.
The land cover data was used to analyze the impact of vegetation and cropland on lake water volume. CLCD data is based on the Google Earth Engine Platform, which adopted the random forest method to classify the Landsat data to generate yearly land cover data of 30 m per year from 1985 to 2020 (
https://zenodo.org/record/5210928#.YrREA3pBztV, accessed on 25 May 2022) [
29]. This data has a long time series and has been widely verified in China. Its mapping accuracy is higher than that of MCD12Q1, GlobeLand30 and ESACCI_LC data of the same type.
2.3. Analysis of Change Trend
The Mann-Kendall (MK) test was adopted to analyze the trend of changes and the abrupt change characteristics of each element factor. This method is a nonparametric statistical test method, which does not rely on the characteristics of normal distribution and is insensitive to a few outliers. It has been widely used to analyze the trend of changes of meteorological and hydrological elements [
30,
31]. Assuming that the
variables of a time series element
are
, the statistics
can be defined as follows:
variance
is as follows:
where
is the number of groups of repeating elements in
, and
is the number of repeating elements in Group
.
The test statistics
are as follows:
During the test, the original hypothesis was set to no linear trend, and the alternative hypothesis was considered a linear trend. The two-sided test was adopted. Under the significance level of , was rejected, and was used in case of . Among them, if the conditions were met, indicates an increasing trend, while indicates a decreasing trend; .
The abrupt change point was further determined for the time series. An order column
was built to count the cumulative number of
in the time series:
Assuming that the time series were random and independent, the statistics
were defined as follows:
where
;
.
The time series was arranged in a reverse order, and the process above was repeated, which resulted in the following equation:
At the 0.05 significance level, the two curves and and two straight lines were drawn. > 0 and < 0 means increasing trend, on the contrary, a decreasing trend. The intersection (if any) between two straight lines and was the abrupt change point.
2.4. Correlation Analysis
The multiple linear regression (MLR) method was adopted to explore the correlation between lake water volume and multiple elements in this study as shown in
Figure 2.
First, a linear regression model between lake storage variables
and other elements (
) was constructed:
where
is the regression constant, and
was the regression coefficient.
The ordinary least square (OLS) method was adopted to estimate the parameters to minimize the square sum of the error between the estimated value and actual value of the lake storage variable:
To determine the minimum value of
, the partial derivative of each element was set to 0 as shown in equation 1.6. The equation can be solved to complete the estimation of multiple linear regression constants and regression coefficients.
The F-test was then adopted to evaluate the overall significance of the multiple linear regression model to ensure the overall significance of the model, and a t-test was used to evaluate the significance of each element on the lake storage variables.
Finally, the backward selection method was adopted to eliminate the elements with significance <0.05. Moreover, previous research indicated that the negative correlation between evaporation and lake storage variable (if taken as the dependent variable) is explainable, namely, more lake water volume can be maintained owing to the reduction in evaporation. If there was a positive correlation, the increase in volume of lake water could lead to an increase in evaporation. Although this indicates the impact of lake expansion, the evaporation in a positive correlation cannot be used to characterize the contribution to lake water volume. Therefore, the evaporation elements, which do not meet the conditions, are eliminated in the backward selection method. In addition, the original elements are standardized in this study, and the coefficient of the final multiple linear regression model is the standardized coefficient Beta. It is possible to determine the contribution or influence of different elements on the lake storage variables by comparing the size of Beta.
3. Results and Discussion
3.1. Comparison and Verification of Lake Water Volume Change
Comparisons showed that the average mean absolute error (MAE) of the change of lake storage volume in TPDC_LWSC and QB_LWSC was between 0.003 km
3 and 0.155 km
3, and Tuosu Lake (2, 0.136) and Dabuxun Lake (9, 0.155) had higher MAE values (
Figure 3). The maximum relative error (AE) is 0.316 km
3, which appears in the comparison of storage variables in Tuosu Lake from 2015 to 2019. The average mean relative error (MRE) was higher and ranged from 47.71% to 1899.91%, primarily because of the difference between the two types of data at low values. For example, TPDC_LWSC and QB_LWSC show that the storage variables of Gasikule Lake (ID:7) were 0.00073 km
3 and 0.081 km
3 from 1990 to 1995, respectively, and the AE was −0.080 km
3. However, the relative error (RE) was as high as 11,101.82%, which resulted in the increase in overall MRE. These values reflect the great uncertainty of both types of data in characterizing small lake storage variables. The root mean square error (RMSE) was between 0.003 km
3 and 0.175 km
3, and the characteristics of its distribution were similar to those of the MAE.
In summary, owing to the difference in data sources and the substantial time span (5a), the lake storage variables obtained from TPDC_LWSC and QB_LWSC have some differences, but almost all show a trend of “similar increase and decrease.” The most obvious relative error of TPDC_LWSC and QB_LWSC occurs at low values, and the AE at the low value is small. Thus, TPDC_LWSC and QB_LWSC are highly consistent. QB_LWSC has a higher time resolution and complete time series compared with TPDC_LWSC and is more suitable to analyze the characteristics of time series and perform a correlation analysis.
3.2. Spatiotemporal Change Characteristics of the Lake Water Volume
Table 3 and
Figure 4 show the characteristics of changes in water volume in major lakes in the Qaidam Basin from 1990 to 2020 over a time series. The change and trend of change clearly show that the total water volume of major lakes in the basin increased by 3.81 km
3 with a rate of increase of 0.79 km
3/a, which presented a trend of a significant increase (QB, Z > 1.96,
p < 0.1). A total of 73.33% of the lakes tended to significantly increase and contributed 4.22 km
3 of the water volume (In
Supplementary Materials Figure S1), which was 110.76% of the net increase. In contrast, the water volume of North Huobuxun Lake (ID 19, Z < −1.96,
p > 0.1) and West Taijinar Lake (ID 24, Z < −1.96,
p > 0.1) decreased by 0.17 km
3 and 0.01 km
3, respectively, which both tended to significantly decrease. However, the changes in Dabuxun Lake (ID 9, Z < 1.96,
p >0.1) and East Taijinar Lake (ID 26, Z < 1.96,
p >0.1) were substantial at 0.36 km
3 and −0.21 km
3, respectively, but the trends of their changes were not significant.
An analysis of the characteristics of changes in the time series showed that the total volume of lake water in the basin changed very slightly from 1990 to 2000, with a decrease of only 0.72 km3 and a rate of change of 0.07 km3/a. In contrast, since the 21st century, the volume of total lake water has increased by 4.53 km3, with a rate of change (0.23 km3/a) that was 3.29-fold of that between 1990 and 2000, and the increment and rates of change were very significant. The average changes in the volume of lake water in the 1990s, 2000s and 2010s were −0.26 km3, −0.71 km3, and 1.36 km3, respectively. It is readily apparent that the volume of total lake water in the basin has been steadily increasing since 2010.
It is apparent from the abrupt changes in the test results, that the total volume of lake water decreased and increased in 1996 and 2014, respectively, but it was not significant at the 0.05 significance level. In terms of specific lakes, some of the lakes abruptly changed points around 1995 and 2015. This was particularly true for Tuosu Lake, Keke Salt Lake and Dabuxun Lake.
In conclusion, the volume of water in the lakes primarily increased from 1990 to 2020 in the Qaidam Basin, and only a few lakes had a trend of changing insignificantly or significantly decreasing. The volume of water in most lakes decreased abruptly around 1995 but increased abruptly around 2015.
Large and medium-sized lakes are primarily distributed in the arid region in the middle of the Qaidam Basin, and small lakes are primarily distributed in the mountains around the basin (
Figure 5). It is apparent from the perspective of spatial pattern in the change in volume of lake water that WRZ I (two lakes) northeast of the basin had the most dramatic change in lake water volume. In 1990–2000, 2000–2010 and 2010–2020, the total lake storage variables in WRZ I accounted for 54.36%, 28.69% and 53.54% of the total lake storage variables in the basin, respectively. The sum of lake storage variables in WRZ VIII, X and XIV (eight lakes) in the middle of the basin accounted for 27.4%, 50.84% and 20.79%, respectively. Thus, there were substantial differences in the changes of lake water volume in different regions of the basin.
3.3. Characteristics of Spatiotemporal Changes of Meteorological Elements
Table 4 and
Figure 6 show the characteristics of change in precipitation in the basin from 1990 to 2020. The trend of change in precipitation clearly showed that the precipitation in the whole basin significantly increased from 1990 to 2020 (QB, z > 1.96,
p < 0.1). Except for WRZ V and XII that tended to increase and WRZ VIII that lacked a significant trend, the other areas tended to significantly increase. The change in precipitation showed that the decades average precipitation of the basin was 162.08 mm in the 1990s. The change increased slightly to 167.47 mm in the 2000s. The decades average precipitation in the 2010s increased significantly to 235.83 mm, which was 1.46- and 1.41-fold of those during the previous two periods, respectively.
The abrupt test results of the changes (the area names correspond to the numbers of WRZ I-XVI from left to right and from top to bottom) showed that there was an abrupt change in the decrease in precipitation in the basin from 1995 to 2000, while an abrupt increasing change appeared from 2015 to 2020. These changes were consistent with the change in total lake water volume in the basin. Among all the WRZ (In
Supplementary Materials Figure S2), except for WRZ III (Chaikai Lake-Keke Salt Lake Water System) and WRZ VI (Dulan Lake Water System), the precipitation in other areas abruptly changed to varying degrees, and all passed the significance test at the level of 0.05. Moreover, the abrupt increase in precipitation in different regions appeared after 2000, particularly in 2015. These results show that the predominant change in the precipitation in the Qaidam Basin was an increase. Particularly since the 21st century, the trend of increase in precipitation has been very significant in the basin.
The spatial pattern of precipitation indicated that the precipitation in the basin was characterized by more precipitation in the southeast, less in the northwest, more in the periphery, and less in the middle, which corresponds to the altitude characteristics of high elevation around the basin and low elevation in the middle (
Figure 7). This is owing to the barrier action of the Qinghai-Tibet Plateau and the obstruction of the Hengduan Mountains to the warm and humid airflow of the southwest monsoons, which results in the formation of a unique spatial distribution pattern of precipitation in the Qaidam Basin.
Table 5 and
Figure 8 show the characteristics of changes over a time series of the actual evaporation in the basin from 1990 to 2020. The trend of changes in evaporation clearly showed that the actual evaporation of the basin increased significantly (QB, Z > 1.96,
p < 0.1), and except for WRZ XIV (Nalinggele Wutumeiren Area), the other WRZ tended to significantly increase (In
Supplementary Materials Figure S3). In terms of change, the annual averages of the actual evaporation of the basin in 1990s, 2000s and 2010s were 93.93 mm, 109.98 mm, and 119.35 mm, respectively, and the growth of the actual evaporation remained relatively stable for 30 years.
There was no abrupt change in the actual evaporation of the basin from 1990 to 2020. With the exceptions of WRZs VI (Dulan Lake Water System), XIV (Nalinggele Wutumeiren Area) and XV (Quanji River/Xiezuo Lake Water System), there was no abrupt change among the WRZs.An analysis of the spatial pattern (
Figure 9) indicated that the spatial distribution of actual evaporation in the basin was consistent with precipitation, but it did not show the obvious characteristics of less in the middle and more in the periphery owing to the fact that the evaporation in the middle was more significant. The evaporation in the middle of basin was significantly greater than the precipitation (
Figure 7 and
Figure 9), indicating the remarkable characteristics of arid/semiarid regions.
The air temperature in the basin did not change dramatically from 1990 to 2020 compared with the precipitation and actual evaporation. However, it still significantly trended upward (QB, Z > 1.96,
p < 0.1) (
Table 6). Except that the changes in air temperature in WRZ V, IX and XI were not significant, and there was a slight uptrend in WRZ X, the other areas tended to significantly increase (In
Supplementary Materials Figure S4). The annual average air temperature in the 1990s, 2000s and 2010s changed by 0.25, 0.55 and 0.80 °C, respectively, with an increasing rate of 0.28 °C/10a. The whole basin and some areas obviously increased in 1995 and 2014.
The results of an abrupt change test indicated that there was no abrupt change in the air temperature of the basin, only a stable uptrend (
Figure 10). However, the air temperature in WRZ II (Chahansilatu Dry Saline Area), WRZ V (Dalangtan-Fengnan Dry Salina Area), WRZ XI (Laomangya Lake Water System) and WRZ XIV (Nalinggele Wutumeiren Area) among each WRZ all abruptly changed and were significant at the 0.05 level (In
Supplementary Materials Figure S4). The abrupt change was concentrated around 1995 and 2015 in which the air temperature of WRZ II and XIV decreased around 1995 and increased around 2015, while the air temperature of WRZs V and Xi reacted in the opposite manner.
The spatial pattern of air temperature (
Figure 11) indicated that the average annual air temperature in Qaidam Basin from 1990 to 2020 was between −9.42℃ and 7.59℃ and were spatially distributed with high air temperature in the middle and low air temperature in the periphery. The rate of change in the air temperature tended to decrease, while that of other areas increased to varying degrees, except for WRZ VII (Gasikule Lake Water System) and XI (Laomangya Lake Water System). There was a relatively high increase in the air temperature of WRZs I, III, VI and XVI northeast of the basin, which were 0.80 °C/10a, 0.61 °C/10a, 0.53 °C/10a and 0.64 °C/10a, respectively. These values were significantly higher than the average change in the rate of air temperature of the basin. There was a relatively low rate of change in the air temperature in the middle and west of the basin.
3.4. Spatiotemporal Changes of Land Cover
Table 7 and
Figure 12 show the characteristics of changes of vegetation area in the Qaidam Basin from 1990 to 2020. The area of vegetation of the basin fluctuated between 66,754.62 km
2 (1996) and 79,883.36 km
2 (2011) from 1990 to 2020. The vegetation area has increased by 1832.38 km
2 over the past 30 years, with a rate of increase of only 2.48%. However, there was still a significant trend of increase overall (QB, Z > 1.96,
p < 0.1).
There was a diverse change in the distribution of vegetation (In
Supplementary Materials Figure S5), and the area of vegetation in 56.25% of the WRZs (nine in total) tended to significantly increase, particularly since approximately 1995. The change in area of vegetation in 31.25% of the areas (five in total) was not significant, and only two areas of vegetation significantly decreased. The annual average vegetation areas in the 1990s, 2000s and 2010s were 70,268.65, 74,805.98 and 78,762.90 km
2, respectively, and the increase was relatively stable.
The results of abrupt change test indicated that there was no abrupt change in the vegetation area of the whole basin and most of the WRZ. The only areas that had significantly abrupt changes in vegetation around 1995 and 2015 were VIII (Golmud Area), XI (Laomangya Lake Water System), XIII (Mahai Lake Water System) and XV (Quanji River/Xiezuo Lake Water System).
Figure 13 and
Table 8 show that the characteristics of cropland area changed from 1990 to 2020 in the Qaidam Basin. The cropland area from 1990 to 2020 fluctuated between 878.73 km
2 (2020) and 1,143.77 km
2 (1990) with a significant tendency to decrease (Z < -1.96,
p > 0.01). Only 62.50% (10) of the cropland area among the 16 WRZs decreased to varying degrees (In
Supplementary Materials Figure S6). There were only significant trends of increase in WRZ II (Chahansilatu Dry Saline Area) and VI (Dulan Lake Water System). In summary, the area of cropland primarily decreased from 1990 to 2020. The annual average of the area of cropland in the 1990s was 1066.53 km
2, while it was 968.19 km
2 in the 2000s, which was a large decrease (9.22%). It increased slightly in the 2010s.
The results of abrupt change test were significant decrease. There was no obvious abrupt change in the overall cropland area of the basin from 1990 to 2020. There were six WRZs that abruptly changed, and the patterns of abrupt change in WRZ II, IV, VI and XIV are relatively similar, which occurred around 1995 and 2015. Abrupt changes in the other areas mostly arose in 2000.
3.5. Responses of Lakes in the Qaidam Basin to Climate Change and Human Activities
3.5.1. Single Factor Correlation of Lake Storage, Climate and Change in Land Cover
Precipitation and evaporation are the primary supply and loss terms of lake water, respectively. Temperatures, on the other hand, indirectly alter the amount of water in lakes by affecting the melting of permafrost and glaciers. Changes in cropland and vegetation area also have potential impacts on changes in lake water storage. Examples include WRZ I, X and XIV of the basin. The water resources in freshwater lakes are used directly for irrigation. The up-stream runoff of salt water or salt lakes is intercepted to irrigate the farmland, and the change in farmland area could have a potential impact on the lake storage. An increase in the area of vegetation can not only alleviate the warming effect of the surface and reduce water evaporation but also reduce the ineffective infiltration and loss of precipitation. This enables the vegetation to contribute to the generation and continuous replenishment of runoff, affecting the runoff into lake.
As shown in
Figure 14. The elements of the whole basin that are associated with lake storage are precipitation (R = 0.80), vegetation (R = 0.66), and evaporation (R = 0.65). In addition, there is a strong positive correlation (
p < 0.001) between evaporation and precipitation (R = 0.66) and evaporation and vegetation (R = 0.74). For each WRZ (In
Supplementary Materials Figure S7), except for the Gasikule Lake Water System and Golmud Area, the precipitation and vegetation in the other areas significantly correlated with the lake water volume, indicating that it was significantly affected by the direct supply of precipitation and vegetation. The evaporation of the Chaikai Lake–Keke Salt Lake Water System (R = 0.64), the Dachaidan Lake Water System (R = 0.74), the Huobuxun Lake Water System (R = 0.60) and the Xiaochaidan Lake Water System (R = 0.68) highly correlated with the lake water volume, which was owing to the increase in lake water volume that led to the expansion of lake water surface, resulting in more intense evaporation. There was a significant positive correlation between the air temperature and lake water volume in the Bayin River/Delingha Area (R = 0.46), Chaikai Lake—Keke Salt Lake Water System (R = 0.63), Nalinggele Wutumeiren Area (R = 0.42) and Xiaochaidan Lake Water System (R = 0.54), and the increase in lake water volume was related to the increase in air temperature. The lakes in these areas are supplied by melted water from glaciers. Combined with the analytical results in
Section 3, the cause was the increase in water from melting glaciers and the frozen soil caused by the increase in air temperature. There are no glaciers distributed in the Gasikule Lake Water System, and its air temperature negatively correlated with the lake water volume (R = −0.36) because the decrease in air temperature inhibited the evaporation of water. There was a strong negative correlation between the cropland of Xiaochaidan Lake Water System (R = −0.74) and the lake water volume. Combined with the analytical results in
Section 3, the reduction in area of cropland led to a reduction in the consumption of irrigation water, which could be the reason for the increase in runoff into the lake.
3.5.2. Degree of the Influence of Climate and Changes in Land Cover on the Lake Water Volume
The correlation analysis of a single factor can only show the influence of a single element on the lake water volume, while the change in lake water volume under particular natural conditions is the result of a comprehensive influence of multiple factors. To further evaluate the influence of different elements on lake water volume, a model of the relationship between lake water volume and multiple factors was established using the method described in
Section 2.3.
Table 9 shows the elements that significantly correlated with the lake storage variables in the basin and their standardized coefficients. Combined with the analytical results in
Section 3, the increase in total lake water volume in the basin significantly correlated with the increase in precipitation and area of vegetation. Precipitation was the main factor (Beta = 0.64), while vegetation was the secondary factor (Beta = 0.40). No significant elements that affect the water volume of the lakes in WRZ VII and VIII have been identified. Precipitation is still an important factor that affects the lake water volume in the other WRZs, and the precipitation in WRZ I (Beta = 0.71) and X (Beta = 0.58) contributes the most to the change in lake water volume. However, evaporation is not a significant element that affects the lake water volume, and it did not pass the significance test in all the WRZs owing to the marked increase in evaporation (
Section 3.3). Although the amount of evaporation can indicate the loss of lake storage to some extent, the increase in precipitation is more significant, and the actual water volume of most of the lakes also increased. Therefore, the “side effects” owing to the increase in evaporation cannot be effectively reflected. The lake water volume in WRZ III and IX increased owing to air temperature. The increase in lake water volume in WRZ III positively correlated with air temperature, but there is no glacier in this area. When the impact of temperature on the water volume of lakes is considered, the increase in water volume of lakes in this area could be correlated with the melting of frozen soil caused by the increase in air temperature. In contrast, the increase in lake water volume in Area IX was related to the decrease in air temperature, which could inhibit the evaporation of lake water. However, the reduction in temperature plays a limited role (Beta = 0.27).
The change in land cover somewhat contributed to the change in lake water volume, but the change in cropland area in each WRZ did not significantly contribute to the lake water volume. The expansion of vegetation area significantly contributed to the lake water volume in 66.67% of the area, and the change in vegetation in the other areas was the main factor except for WRZ X. This shows that when the amount of precipitation increases, an increase in the area of vegetation in the basin helps the precipitation to form runoff, which then promotes the effective supply of lake water volume. In summary, precipitation has made a major contribution to the increase in lake water volume in the whole basin, while vegetation plays an important role in the increase in lake water volume in each WRZ.
3.5.3. Influence of Human Activities on Lake Changes
The lakes in the central Qaidam Basin are the most significantly affected by human activities. Industrial activities, such as brine mining, salt field construction and damming, have led to dramatic changes in the spatial pattern of lakes in the central basin for decades. An analysis of the time sequence showed that there was relatively little industrial activity in the central basin from 1990 to 2000. Except for the slight increase in the water volume of East Taijinar Lake, the water volume of other lakes decreased to varying degrees. The shrinkage of water surface of West Taijinar Lake was the most obvious, which was caused by underground brine mining and the diversion of lake water to East Taijinar Lake.
The total lake water volume in the central basin increased significantly from 2000 to 2010. CITIC Guoan Science and Technology Development Co., Ltd. (Gelmud, China) (2003) and Qinghai Lithium Co., Ltd. (2008) dammed the east and west sides of West and East Taijinar, respectively, to maintain brine mining, resulting in the rapid increase in the water volume of Yahu Lake (25) by 0.33 km
3 [
18]. Part of the lake water flows into the northern area of West Taijinar Lake and forms another new lake (Yiliping Lake, 23). This new lake was the largest in the basin in 2010, with an increase of 0.29 km
3 in water volume. In addition, dams have been built in the west of Senie Lake (28) and Dabsan Lake (9), which formed the two lakes Noname2 (27) and Noname3 (8), respectively.
China Minmetals Corporation (Beijing, China) (2013) dammed south of Yiliping Lake from 2010 to 2020 to block water and runoff into the lake to mine underground brine and build salt fields. Since then, Yiliping Lake has rapidly dried up, and the lake water volume decreased by 0.29 km3 in one year. Except for the western salt field, Yiliping Lake has almost completely dried up. Simultaneously, the dam to the west of Yahu Lake blocks a large amount of incoming water, and its water volume increased rapidly, with a rate of increase of 0.08 km3/a and an increase of 0.81 km3. In contrast, the damming and salt field construction around Dabsan Lake had little impact on its water volume, indicating that the salt lake mining in this area is relatively reasonable for the area.
In summary, brine mining and damming activities in the central basin are the primary factors that affect the change in lake water volume. Lakes not only function to regulate the regional microclimate; they are also important resources for local economic development. The development of salt lakes will inevitably bring challenges to the sustainable development of lake resources and the stability of regional climate. The question of how to coordinate ecological protection and economic development and formulate reasonable policies for salt lake development from the perspective of sustainable development are major challenges faced by the local authorities.