Abstract
Baseflow plays a crucial role in sustaining the alpine ecosystem during rainless or cold periods. Despite its importance, information on how and why baseflow has changed in the source region of the Yangtze River (SRYR) is sparse. In our study, statistical analysis and the elastic coefficient method were used to identify the dynamic characteristics of baseflow and the underlying causes. The results show that monthly baseflow contributed 62–97% of runoff with a mean value of 75%, and they followed remarkable increasing trends from 1957 to 2020. The contributions of precipitation, temperature, evapotranspiration, and ecological conservation programs (ECPs) on baseflow variations were 86, 53, −15, and −24%, respectively. However, their contributions differed across months. During the warm months of May to September, precipitation played a dominant role, followed by evapotranspiration. In contrast, during other colder months, temperature was dominant; meanwhile, the effect of precipitation was almost absent. Moreover, climatic change had a hysteretic effect on baseflow variation, with a maximum lag time of 10 months. Our results highlighted critical roles of both precipitation and temperature, and indicated that climate change, rather than ECPs, dominated the variation in baseflow in the SRYR.
HIGHLIGHTS
Monthly baseflow contributed 62–97% of runoff with a mean value of 75%.
Contributions of precipitation, temperature, evapotranspiration, and ecological conservation to the increasing baseflow were 86, 53, −15, and −24%, respectively.
During warm months, precipitation played a dominant role, while during colder months, the impact of temperature was dominant.
Climatic change had a hysteretic effect on baseflow variation.
INTRODUCTION
The health of river ecosystems is undoubtedly associated with water cycle processes at the basin level (Kayitesi et al. 2022). Groundwater flow feeding streams, referred to as baseflow, are an important portion of the watershed hydrologic cycle. Baseflow is the reason many streams do not dry up during dry periods, and it can provide a persistent source of water to streams, sustaining life living in and around them (Guo et al. 2022). However, climate variability is expected to intensify or accelerate the hydrologic cycle (Yin et al. 2023a, 2023b, 2023c), thereby causing an alteration of the baseflow. At the same time, anthropogenic activity is likely to modify the baseflow regime of watersheds (Chen & Teegavarapu 2021). Prior studies have noted that the baseflow trends of many streams across the globe have experienced profound changes (Ficklin et al. 2016). Thus, ensuring sufficient groundwater availability is greatly challenged by climate change and human activities (Scanlon et al. 2023). Studying baseflow dynamics and the underlying causes will be important in understanding the evolution of the hydrologic cycle and groundwater storage in changing environments (Brutsaert 2008).
Accurate baseflow estimation is critical for the quantitative management of watershed ecological protection/conservation and the utilization of water resources (Smakhtin 2001). However, few techniques exist to help water managers directly monitor baseflow changes in streams. Numerous baseflow separation approaches have been developed to estimate the baseflow from the streamflow time series. They can be broadly classified into four major categories: hydrograph separation (e.g., HYSEP, flow recession analysis; Ficklin et al. 2016; Cheng et al. 2017), hydrological model (e.g., SPARROW; Miller et al. 2016), numerical simulation (e.g., two-parameter digital filtering method; Eckhardt 2005), and water and chemical mass balance (e.g., geochemical tracers; Tetzlaff & Soulsby 2008). The HYSEP method, as a traditional baseflow segmentation method, is simple but labor-intensive when dealing with long-term series data, and its result is highly subjective (Sloto & Crouse 1996). The physical interpretations of the hydrological model are relatively clear, but its model construction and parameter calibration are complicated and tedious, and the results have a certain degree of uncertainty (Kayitesi et al. 2022). The chemical mass balance method separates baseflow by analyzing the composition of stable isotopes in rainfall, surface water, and groundwater, and its results are more reliable. However, it requires large amounts of manpower and material resources and is difficult to widely use (Lott & Stewart 2016). Numerical simulations have been applied worldwide to the baseflow separation of long-term hydrologic series because they are easily operated, efficient, practical, and objective, avoiding the arbitrariness of manual methods (Ahiablame et al. 2017). For example, Murray et al. (2023) used the Lyne–Hollick one-parameter digital filter method to estimate baseflow and characterized the impacts of climatic change on monthly baseflow trends in Canadian watersheds. However, the baseflow separation results obtained by the above methods may vary, and sometimes the differences are significant. To reduce the uncertainties of the results, it is preferable that multiple baseflow separation methods are adopted simultaneously and qualitatively compared with each other based on their hydrologic plausibility. This serves to select the most suitable method.
With such leaps in understanding the importance of baseflow to stream water dynamics, more research has focused on exploring the factors that control baseflow and how baseflow responds to said factors. In general, the discharge rate of groundwater mainly depends on basin characteristics, climate regime, and human intervention (Segura et al. 2019). Watershed characteristics experience minute changes over the course of decades, and as a result, they are largely unnoticeable. Temporal variations of baseflow yield are mainly impacted by climate variability and human intervention (land cover change, increased water utilization, etc.). For example, the reduction of baseflow, which was considered to be the direct reason for the summer drying of the Ganges river, was found to be triggered by groundwater depletion (Mukherjee et al. 2018). Groundwater withdrawals were likely accelerating the reduction of baseflow in the heavily irrigated Sao Francisco River basin in Brazil (Lucas et al. 2020). Zomlot et al. (2015) used the WetSpass water balance model to estimate the recharge rate of 67 hydrologic stations in Flanders, discovering the prominent role of land use changes on recharge variation. Ficklin et al. (2016) researched baseflow for 674 sites across the United States, finding that baseflow displayed strong associations with precipitation trends, followed by evapotranspiration. At one point, there was a hypothesis that the baseflow fraction generated by precipitation is primarily attributed to the aridity index, but great differences occurred in humid and arid catchments (Gnann et al. 2019). This influence mechanism can be attributed to distinct antecedent wetness. A great many trends in baseflow across Canada were not related to any climatic predictors (Murray et al. 2023), likely due to intense anthropogenic activities, such as catchment urbanization, reservoir operation, and river abstractions for industrial, and agricultural, or municipal purposes (Smakhtin 2001). Through summarizing previous studies, we found that many of them were limited to qualitative evaluation about the cause of baseflow change based on correlation analysis. The quantitative attributions of baseflow relative to climate variables and human activities are still rarely discussed, particularly for basins in alpine regions, where temperature plays an important role in driving baseflow variation.
The source region of the Yangtze River (SRYR), located in the hinterland of the Qinghai–Tibet Plateau, is an important part of the Asian Water Tower and China's ecological security barrier. Over the last several decades, the SRYR has suffered unprecedented climate warming and intensified human activities (Cui et al. 2023). Consequently, climate change will likely lead to changes in the water cycle processes of the basin, which will then have adverse effects on river ecosystems in the SRYR (Rounce et al. 2023). The consistently rising trend in runoff and its related causes have been widely reported (Liu et al. 2019; Li et al. 2023a, 2023b). However, in-depth studies on the dynamic variations and attributions of baseflow have not been adequately conducted, and the understanding of influence mechanisms in alpine areas remains in the exploratory stage. The primary objectives of this study are to (1) compare the results of multiple baseflow separation methods and select the most appropriate one, (2) clarify the temporal variations of baseflow at various time scales and investigate the associations of baseflow with trends in specific climates, and (3) quantify the relative influences of climate variables (e.g., precipitation, evapotranspiration, temperature) and human activities on baseflow. This study is designed to advance the understanding of baseflow variation and its causes in alpine areas, and to guide future activities regarding water resources management and ecosystem restoration/protection within the SRYR.
MATERIALS AND METHODS
Study area
The vegetation that covers the headwaters primarily consists of three types: alpine meadow, grassland, and marsh. Under the joint influences of climate change and anthropogenic disturbance, problems like accelerated permafrost melting, grassland degradation, soil erosion, and land desertification directly threaten the stability of the rivers' ecosystems and has recently been recognized as a serious social concern (Li et al. 2023a, 2023b). Moreover, supra-permafrost water has become a major source of surface water at the Zhimenda and Tuotuohe stations (Li et al. 2020), implying that the Yangtze River Water Tower has become increasingly unstable. To prevent further degradation, in 2005, China's State Council invested 7.5 billion yuan RMB to implement more than 22 ecological conservation programs (ECPs) in the Sanjiangyuan Nature Reserve, covering the headwaters of the Yangtze River, the Yellow River, and the Lancang River (Luo 2022).
Data and source
The SRYR occupies a large area of inaccessible land, so the distribution of hydrological and climate gauges is sparse. Until now, only six hydrological stations have been distributed in the SRYR: the Zhimenda, Tuotuohe, Xinzhai, Yanshiping, Qumalaio, and Longbaotan stations. Unfortunately, all stations, except Zhimenda, have produced patchy observation records. As Zhimenda's hydrometric station is located at the outlet of the SRYR, its daily discharge records were selected for analysis of baseflow characteristics. This station had available data dating back 64 years (1957–2020), and these data were gathered from the Qinghai Hydrological and Water Resources Measurement and Reporting Centre. In addition, we collected meteorological records from five other stations: Tuotuohe, Wudaoliang, Qumalai, Yushu, and Chenduo. These data were provided by the National Meteorological Information Center. These climate data included daily precipitation, relative humidity, and maximum, mean, and minimum air temperatures. The daily potential evapotranspiration data were estimated using the Penman–Monteith approach. The five climate stations mentioned above are representative of the SRYR, and their weighted averages were calculated using the Thiessen polygon. These climate observations provided context for discussing the potential influences of climate change on baseflow. The SRYR topographic map, along with the hydrological and meteorological stations, is presented in Figure 1.
Furthermore, the thawing index was chosen to predict the seasonal ground freeze–thaw and permafrost thaw in the SRYR (Woo et al. 2004). The thawing index (°C days) is calculated for each year by summing up the mean daily air temperatures for all the days exceeding 0 °C (Duan et al. 2017). Observed groundwater level records from 2010 to 2020 at the Zhimenda station were selected to reflect groundwater storage variability, and glacier area change data referenced from Li et al. (2023a, 2023b) were used to reflect the changes in warming-induced glacier meltwater in the SRYR.
Baseflow separation approaches
A great many baseflow separation methods have been developed to derive baseflow from the observed streamflow time series. However, applying a single method is risky. Thus, the use of multiple methods to calculate the baseflow is recommended (Tan et al. 2020). To obtain accurate estimates of baseflow, seven methods were adopted in our study, including one digital filter method, two baseflow index (BFI) methods, three HYSEP methods, and one PART method. These methods are detailed in the following sections.
Baseflow separation approaches
BFI methods
The BFI procedure was developed by the Institute of Hydrology, Wallingford, United Kingdom, and then programmed by Wahl & Wahl (1995). The BFI method, also known as the minimum sliding method, uses the BFI as a weighing factor to calculate the baseflow. The BFI method mainly includes the standard BFI(f) method and the modified BFI(k) method. Typically, rivers that are recharged by groundwater have a BFI close to 1, while intermittent rivers have a BFI value tending toward 0.
HYSEP methods
The HYSEP software (Sloto & Crouse 1996) uses three methods originally developed by Pettyjohn & Henning (1979) to separate baseflow and runoff components from streamflow hydrographs: the fixed interval method (HYSEP1), the sliding interval method (HYSEP2), and the local minimum method (HYSEP3). These methods use different algorithms to draw connecting lines (baseflow hydrograph) between low points of the streamflow hydrograph. The detailed description of these methods can be seen in Pettyjohn & Henning (1979).
PART methods
The PART method is a procedure proposed by the U.S. Geological Survey for baseflow separation. Baseflow is considered equal to streamflow on days that are designated as being unaffected by surface runoff. The PART method requires daily hydrological series with no less than one hydrologic year. The separation is based on the receding of the antecedent streamflow, which first arranges the daily flow data into a single-dimensional array and then selects the values that meet the recession requirements. If the daily recession has less than a 0.1 log period, the daily streamflow value is used as the baseflow. During the remaining time, the baseflow values are obtained by linear interpolation.
Baseflow evaluation criterion
Since the applicability of the baseflow separation methods differed in various regions, it was necessary to evaluate their practicalities in the SRYR. Considering that baseflow is relatively stable, the reliability of annual separated baseflow can be assessed using the standard deviation and extreme value ratio of the BFI series. For the purposes of this study, the ratio of annual maximum minus annual minimum to annual minimum was calculated as an extreme value ratio. In addition, strict baseflow points were selected as baseflow references (considered to be true values). We evaluated the accuracies of seven baseflow separation methods using the Nash–Sutcliffe efficiency (NSE) and Kling–Gupta efficiency (KGE). Strict baseflow points can be obtained using the following four rules: (1) Erase all the data points of daily streamflow with , where ; (2) Eliminate the previous two points before points with , as well as the next three points; (3) Eliminate five points after major events that were identified by flood peaks greater than the 90th quantile of all streamflow observations; (4) Exclude data points followed by a data point with smaller , namely, (Xie et al. 2020). Detailed explanations of the four rules can be found in Xie et al. (2020).
Trend and breakpoints detecting
The trend in baseflow and its corresponding predictors at annual/seasonal scales from 1957 to 2020 can be estimated using the nonparametric Mann–Kendall test, Sen's Slope test, and linear regression analysis. The Mann–Kendall test was employed to detect whether there are monotonic trends with a significance level of p = 0.05 in such long-term time series. The slope of the trend was calculated using Sen's Slope test.
Pettitt's test is widely used for detecting mutations. In our study, Pettitt's test, as well as the Mann–Kendall test, was applied to assess if there were any mutation point in the time series. The results from the two test methods were validated by each other.
Correlation analyses
The correlation coefficient r ranges from −1 to +1. A negative r indicates a negative linear relationship between x and y, and a positive r indicates a positive linear relationship between x and y. The greater the absolute value of r, the stronger the correlation.
Attribution analysis of baseflow
RESULTS
Selection of baseflow separation methods
In summary, the Eckhardt method was the most suitable and reliable method for baseflow separation in the SRYR. The upcoming sections will discuss the baseflow and BFI time series derived from the Eckhardt method employed to further analyze the temporal change in baseflow and its potential causes.
Temporal variability of annual and seasonal baseflow and BFI
Notes: NS = non-significant trend. March to May, June to August, September to November, and December to February were selected to represent spring, summer, fall, and winter, respectively. April to September and October to March were chosen to represent warm and cold months, respectively.
Though the annual BFI, in both warm and cold months, had gradual increasing trends from 1957 to 2020, their trends were not statistically significant. Spring and autumn BFI showed almost no pronounced monotonic trends, while the BFI in summer and winter exhibited statistically insignificant increasing trend. In summary, BFI had a much weaker upward trend compared to baseflow.
The mutation detecting process showed that the abrupt year in annual baseflow was 1998, and those in spring, summer, autumn, and winter were 2004, 1998, 2002, and 2004, respectively. Detecting the temporal pattern in BFI was unsuccessful; only abrupt points in annual and winter baseflow could be detected by Pettitt's test, and they were found in the years 2000 and 2006. The change of BFI obviously lagged behind that of baseflow, with a delay time of 2 years.
Intra-annual distribution of baseflow and BFI
By contrast, BFI in summer and spring were nearly identical and lower than that found in autumn and winter. Compared to the cold months, warm months had a lower BFI. In addition, the BFI had two peaks, one of 0.74 in mid-May and the other of 0.99 in late October. Clearly, the most streamflow in cold months was discharged by groundwater. The mean BFI on April 19 and June 16 reached the two relatively lowest values, 0.56 and 0.55, respectively.
Temporal variability in meteorological factors
Pettitt's test demonstrated that precipitation and temperature exhibited significant changes in 1998 and 1997, close to the year of the mutation in baseflow and BFI. Nevertheless, a much earlier breakpoint of evapotranspiration was detected in 1969. It can be inferred from the facts above that precipitation and temperature, not evapotranspiration, played an important role in the long-term evolution of the baseflow in the SRYR.
Correlation of baseflow with meteorological factors
Correlations between baseflow and climate parameters varied across seasons and months. All year round, precipitation, temperature, and evapotranspiration played different roles in controlling baseflow. For spring and summer, precipitation made the greatest contribution with correlation coefficients of 0.54 and 0.77, followed by evapotranspiration (−0.38 and −0.58) and then temperature (0.17 and 0.21). For autumn, precipitation and temperature made the top two contributions with correlation coefficients of 0.60 and 0.40, followed by evapotranspiration (−0.13). During winter, temperature made the greatest contribution, with a correlation coefficient of 0.50, followed by evapotranspiration (0.33). However, precipitation had no relationship to the fluctuation of baseflow. Specifically, from October to April of the following year, temperature played a dominant role in driving the increase in baseflow. With April as the exception, the influence of precipitation was almost nonexistent during this period. However, during May to September, precipitation was dominant, followed by evapotranspiration, while temperature's impact was negligible. In other words, during the warm months, precipitation (0.84) made more contribution than evapotranspiration (−0.50) and temperature (0.40). In the cold months, temperature (0.58) made the greatest contribution, compared to evapotranspiration (0.29) and precipitation (0.19).
As seen in Figure 7, there was a response delay in the baseflow caused by meteorological factors, and this response process varied for different periods. There was a significant correlation between baseflow and precipitation in the months of April to September, but their strongest correlation occurred in the second month. This hysteresis of baseflow to precipitation can last several months. For example, the precipitation in September maintained a large influence on the baseflow in the following 10 months. Evapotranspiration exhibited a moderate relationship with baseflow in those same months, and its strong relation was found in the second month, but soon ceased in the third month. Evidently, the delay of baseflow to evapotranspiration lasted only for 1–2 months. Moreover, during January to April and September to October, the strong relationship between temperature and baseflow did not cease until May, and in November, the strong relationship persisted until the second month.
Contribution of meteorological factors to increase in baseflow
Based on the results of attribution analysis, climate variability contributed the most, with a percentage of 124%; precipitation contributed 86%, temperature 53%, and evapotranspiration −15%. In contrast, the impact of human activity had a small proportion of −24%.
DISCUSSION
Critical role of groundwater to streamflow
Groundwater is a vital source of freshwater, and as such, maintaining and enhancing groundwater storage and discharge is central to climate adaptation (Stigter et al. 2023; Yin et al. 2023a, 2023b, 2023c). Ahiablame et al. (2017) pointed out that 60% of streamflow in the Missouri River basin is derived from baseflow. The simulation results from Miller et al. (2016) show that 56% of the surface water in the Upper Colorado River Basin originated as baseflow. Aquifers beneath our own study area contain a large amount of solid and liquid groundwater, but its hidden role remains unclear. Compared to the rivers in temperate regions, the SRYR has a much greater proportion of groundwater discharge. Although the streamflow is composed of multiple potential sources, groundwater is the major contributor to runoff generation. Groundwater contributed 62–97% of the monthly mean streamflow, and 75% of the annual mean, indicating a significant role in sustaining the discharge of rivers in the SRYR. Gonzales et al. (2009) found a much higher proportion of groundwater in the streamflow in cold regions. For a small catchment in the central part of the Netherlands, approximately 90% of the total runoff comes from groundwater. Our findings corroborated the ideas of Guo et al. (2022), who suggested that groundwater contribution is much more important to the glaciated alpine watersheds on the Tibetan Plateau than previously thought. As an underlying system mediator, groundwater can provide relatively stable water sources against extremes, especially in the face of global climate change (Condon et al. 2020).
Lagged effects of baseflow in response to meteorological factors
Baseflow itself is a delayed source, sustaining streamflow in an extended period after precipitation or during periods of drought (Ayers et al. 2021). Baseflow is directly associated with precipitation, but there are apparent time lags between the two. Changes in precipitation may drive baseflow discharge by altering the gradient between surface water and groundwater levels. However, the points at which precipitation inputs arrive at surface water and groundwater bodies are vastly different due to their discrepancies of recharging routing. This reflects, in part, the different time scales over which surface water and groundwater operate. This is why there was a delay between baseflow and precipitation in our study. Meanwhile, the hysteresis effect of baseflow on precipitation is diverse across regions. For example, the hysteretic influence of precipitation on baseflow was evident in the US Pacific Northwest and Midwest but was not the case in the remaining areas (Ayers et al. 2022). In addition, the precipitation hysteretic effect of baseflow may vary by month. Our results, as shown in Figure 7, also reflected this diversity. The influence of precipitation on baseflow in the SRYR was large during April to September and could be extended for no less than 3 months, but that extension cannot be applied for other months.
Our results also revealed that the delay effect for temperature was similar to that for precipitation, but their influencing mechanisms were completely different. When the daily maximum temperature is greater than 0 °C, the melting of frozen soils or glaciers begin to intensify. However, it will take months for the soils or glaciers to melt, due to slow heat conduction through soils and ice, as well as production and confluence processes. This causes baseflow to lag behind temperature. However, it is noted that the correlation and hysteresis effect of baseflow and temperature were absent from May to August because of the disturbance of precipitation inputs.
The influence mechanism of evapotranspiration on baseflow was similar to that of precipitation, but the impact of evapotranspiration was not significant and lasted for only 1–2 months. This was likely perturbed by the combination of precipitation and temperature.
Climatic contribution of baseflow variations
A great amount of knowledge has been accumulated on baseflow in the past two decades. Despite this, our understanding of baseflow-generating mechanisms remains rather limited. This limitation is particularly pronounced with alpine areas, where the influence of climate change on baseflow variations is much more intricate (Sun et al. 2020).
Precipitation generally plays a dominant role in driving the baseflow variability in unregulated rivers in a variety of climate regions (Miller et al. 2016). Globally, precipitation was the most prominent cause for temporal changes in mean runoff for 83% of land grid cells. Here, we have identified that increased precipitation was the most important factor in driving baseflow variations, with a contribution of 86%. This is somewhat smaller than that in other more temperate regions, likely because only some months out of the year experienced precipitation domination during the formation of baseflow. The role of precipitation inputs in the increasing trend of baseflow was lost in the cold months, because during these cold months, the dominant form of precipitation was snow, and the incidence of precipitation recharging baseflow was reduced. It is clear that precipitation made minimal contributions to baseflow during the winter months, but the long-term mean winter baseflow trend was somewhat influenced by precipitation inputs (Paznekas & Hayashi 2016). Certainly, our results were not always consistent with previous studies. For example, the investigation of Paznekas & Hayashi (2016) showed that precipitation inputs had a strong correlation with the long-term mean winter flow in the Canadian Rockies, reflecting the dominant role of precipitation during wintertime.
A number of researchers have reported that precipitation and evaporation cannot fully explain changes in the baseflow in alpine areas (Ahiablame et al. 2017). Temperature is an important predictor in high elevations (Ayers et al. 2021). For rivers with melting water from glaciers and permafrost as the main source of recharge, temperature often serves as the primary factor influencing runoff variations (Lin et al. 2020). It is worth noting that temperature may also influence baseflow variations by altering the spatiotemporal patterns of precipitation (Yin et al. 2023a, 2023b, 2023c), but this is beyond the scope of our study. Many studies have suggested that solid water is the most important source of river water during ablation episodes in the SRYR. Li et al. (2020) pointed that the proportion of precipitation for rivers in the SRYR permafrost area was only 20.79%, whereas those of supra-permafrost water and glacier/snow meltwater were 69.54 and 9.67%, respectively. Our study area experienced more rapid changes in temperature than environments at lower elevations (Pepin et al. 2015). Our results also show that upward trends in baseflow coincided with an increased temperature. Thus, another reason why baseflow was always increasing is that temperature experienced a marked rising trend. In our study, temperature was a secondary contributor. By comparison, Yi et al. (2021) did not agree with our conclusion. They thought that rising temperature was the primary cause. Regardless, temperature played an important role in the generation of baseflow in the SRYR.
By comparison, evapotranspiration was the least important factor to negatively impact baseflow variations, with a small contribution of −11%. That being said, evapotranspiration may sometimes act as a dominant contributor (Smakhtin 2001). For example, as a predictor of baseflow, terrestrial water storage throughout the Arctic experienced a declining trend, and Suzuki et al. (2018) largely attributed it to an increase in evapotranspiration driven by rising summer air temperatures.
Possible role of anthropogenic activities
Perhaps the greatest obstacle in assessing baseflow response to climate change is the confounding factor of concurrent anthropogenic activity. Anthropogenic activities in the SRYR included the implementation of ECPs, water withdrawal for domestic purposes, and grazing-induced land cover change. However, this source area was recently a depopulated zone because of poverty and harsh environments. More and more residents moved away from this area after 2000 when massive ecological migrations were initiated. Considering that the population in headwaters was sparse, the influences of water withdrawal and grazing were approaching zero and can consequently be neglected in the attribution of baseflow. Therefore, the ECPs may represent much of the human activity in this area.
The ECPs in this area, including Wetland Protection Program and Desertification Land Control Program, greatly changed land use and land cover (LULC), producing significant ecological effects. The evidence that vegetation coverage is continuously increasing in the SRYR in recent years has been widely reported. According to the Statistical Bulletin of National Economic and Social Development of Yushu Prefecture, more than 130 million sporadic plants and 200,000 hectares of grasslands have been developed. Rodents have been wiped out in more than 5,000,000 hectares since the implementation of ECPs. As of now, water conservation has increased by more than 6% annually, and grassland coverage and grass production have risen by more than 11 and 30%, respectively, compared to 10 years ago. Certainly, vegetation cover could be an integrated index (Tan et al. 2020), and its development was prompted not only by the implementation of ECPs but also by the rising temperatures, precipitation, and atmospheric CO2 (Zhang et al. 2022). The LULC change reflected by these statistics may lead to the alteration of the water cycle and water retention capacity, thus directly affecting the baseflow-generating processes. Our findings show that the ECPs essentially representing the impact of human activities can explain −24% of the increasing baseflow trends, indicating a negative relationship between baseflow and ECPs. This is somewhat consistent with results from other studies (Brown et al. 2005). A large number of studies on the impact of land cover change on baseflow revealed a statistically significant negative association between changes in baseflow and basin vegetation (Tan et al. 2020). Essentially, evapotranspiration was limited by water availability (Babst et al. 2019). The additional grasslands may improve water retention by trapping precipitation in the vegetation canopy, litter, and soil layers, resulting in a persistent increase in the actual evaporation rate (Trancoso et al. 2017). Soil drying would be exacerbated, ultimately leading to a reduction in baseflow. Additional vegetation cover may improve the transpiration rate, thereby greatly consuming shallow groundwater (Condon et al. 2020). Thus, this negative relationship between catchment vegetation cover and baseflow volume may be due to greater interception and evapotranspiration rates associated with vegetation cover.
CONCLUSIONS AND FUTURE PROSPECTS
Conclusions
The main purpose of this study was to investigate dynamic characteristics of baseflow variations and identify possible climatic and anthropogenic drivers across the SRYR. Sixty-four years of hydrological records were gathered from the Zhimenda gauging station at the catchment outlet to investigate the changes in annual and monthly baseflow and BFI across the whole watershed. Statistical monotonic trends and corresponding trend slopes were determined using the Mann–Kendall test and Sen's Slope test, and breakpoints were detected using Pettitt's test. Correlation analysis and the elasticity approach were carried out to study the contributions of climate change and human activity to changes in baseflow.
The Eckhardt method was assessed to be the most suitable approach in the SRYR. Annual and seasonal trends in baseflow and temperature were consistently increasing, which was not the case for precipitation and evapotranspiration; annual precipitation and evapotranspiration had a rising trend, but only springtime precipitation and wintertime evapotranspiration increased. In addition, time lags between baseflow and climate parameters were found, and the hysteretic effects varied for different climate parameters and months.
We found that the roles of precipitation, temperature, and evapotranspiration were distinct in spring, summer, autumn, and winter. Of the three climate predictors, precipitation was the most critical driver in contributing to the change in baseflow, followed by temperature, and then evapotranspiration. The contributions of precipitation, temperature, and evapotranspiration on baseflow were 86, 53, and −15%, respectively. However, across the four seasons, the three predictors played different roles in influencing baseflow. Precipitation made the greatest contribution in spring and summer, followed by evapotranspiration and temperature. In autumn, precipitation and temperature made the greatest contributions, followed by evapotranspiration. In winter, temperature made the greatest contribution, followed by evapotranspiration. The contribution of precipitation was so minute that it can be overlooked. On a monthly scale, from October to April of the following year, temperature played a dominant role in the evolution of baseflow, while the effect of precipitation was almost absent (except in April). During the period from May to September, precipitation was dominant, followed by evapotranspiration, while the impact of temperature was negligible. In addition, ECP-induced vegetation cover greening made a negative contribution of −24% to the increase in baseflow.
All in all, the glacier/permafrost melting water and snowmelt associated with increasing temperature, in addition to the increasing precipitation inputs, drove the increase in groundwater discharge, albeit under the negative influence of evapotranspiration and ECP implementations. The variation of baseflow was the result of the combined effect of the aforementioned factors. Clearly, the increased temperature superimposed changes upon the precipitation effect. The hydrologic effect of precipitation could be exacerbated by temperature change in the future.
Future prospects
Although our results highlighted temporal patterns in baseflow and its attributions, the spatial variability was less seriously considered. The baseflow-generating mechanisms vary with elevation, and patterns of baseflow variability may be of small difference across the source region. In the future, a distributed hydrological model coupling the melted ground and glacial snow should be constructed to investigate the spatial pattern of baseflow.
Furthermore, with the steady progress of the pilot program for the Sanjiangyuan National Park and the comprehensive implementation of the Qinghai–Tibet Plateau Ecological Protection Law, the vegetation in the SRYR will experience long-term positive development. This will lead to a further enhancement of water source conservation capacity, profoundly affecting the future trends in baseflow changes. Moreover, according to the Climate Change 2022: Impacts, Adaptation and Vulnerability released by IPCC, under the backdrop of climate change, future runoff variations are uncertain and complexities are heightened, emphasizing the more pronounced role of baseflow in stabilizing runoff. Hence, accurate baseflow prediction is crucial, necessitating comprehensive consideration of vegetation cover and the impacts of climate change.
AUTHOR CONTRIBUTIONS
GW contributed to the methodology, formal analysis, writing – original draft, writing – review & editing, and funding acquisition. JZ contributed to the methodology, conceptualization, and supervision. YL contributed to the data curation and validation. YL and MY contributed to the investigation and validation. HR contributed to the visualization.
ACKNOWLEDGEMENTS
This research was jointly funded by the National Key Research and Development Programs of China (No. 2022YFC3201700), the National Natural Science Foundation of China (No. 52009006), and the Fundamental Research Funds for Central Public Welfare Research Institutes (No. CKSF2023315/TB). We also acknowledge the Qinghai Hydrology and Water Resources Survey Center for data support.
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.
CONFLICT OF INTEREST
The authors declare there is no conflict.