Decomposition and Attribution Analysis of Runo ﬀ Alteration of the Dongting Lake in China

: The runo ﬀ process in the Dongting Lake has been inﬂuenced by climate change and human activities in recent decades. To manage the Dongting Lake e ﬃ ciently and exploit water resources properly under the background of water shortage, it is desired to detect the factors of runo ﬀ change in the Dongting Lake. Hydro-meteorological data from 1961 to 2019 are analyzed to reveal the climate change and runo ﬀ alteration of the Dongting Lake comprehensively. Mutation test is used to detect the change points of runo ﬀ depth series, ﬁnding that 1984 and 2005 are change points and therefore 1961–1983, 1984–2004, and 2005–2019 are regarded as baseline period (BP), period 1 (P1), and period 2 (P2), respectively. Eight methods are used to quantitatively assess the relative contribution of human activities and climate change on runo ﬀ variation. It reveals that climate change especially precipitation change plays the dominant role (climate change makes runo ﬀ depth increase 70.14–121.51 mm, human activities make runo ﬀ depth decrease 51.98–103.35 mm) in runo ﬀ alteration in P1 while human activities play a prime role (account for 88.47–93.17%) in P2. Human activities such as reservoir construction, water consumption, and land-use (land-cover) change may be the main factors that inﬂuence the runo ﬀ in the Dongting Lake in P2. According to the sensitivity analysis, runo ﬀ in the Dongting Lake is more sensitive to climate change in P2 compared with that in P1, and no matter in P1 or P2, runo ﬀ is more sensitive to change in precipitation than the change in potential evapotranspiration. Combined with climate forecast, the results of sensitivity analysis can be used to estimate runo ﬀ change caused by climate change in the future.


Introduction
Runoff, an important source of social water use, has changed in many regions of the world in recent decades [1]. The factors for runoff change can be divided into two groups: climate change and human activities [2][3][4]. Therefore, a lot of studies have attempted to estimate the effects of climate change and human activities on runoff [5][6][7]. There are four kinds of methods commonly used to separate the effects of human activities and climate change, namely empirical statistics, elasticity-based method, hydrological modeling, and conceptual model [8,9]. Empirical statistics usually establishes a linear relationship between runoff and precipitation. If there is no influence of human activities on runoff, the relationship between runoff series and precipitation series will not change. The most frequently used empirical statistics methods are the simple linear regression (SLR) method and the double-mass curve (DMC) method [10,11]. It requires long-time hydrological and meteorological data. The elasticity-based method was proposed by Schaake [12]. It relays on the climate elasticity coefficients to represent the sensitivity of runoff to variations in meteorological factors, then runoff

Study Area
The Dongting Lake (28°30′-30°20′N, 111°40′-113°10′E), with an area about 2625 km 2 , is the second-largest freshwater lake of China located in north-eastern Hunan province, and on the south bank of the middle reaches of Yangtze River (Figure 1). It is composed of three sub-lakes, named the Eastern, Western, and Southern Dongting Lakes, respectively. Among them, the Eastern Dongting Lake is the largest one. The river-lake system of the Yangtze River and the Dongting Lake is complex. The water of the Yangtze River pours into the Dongting Lake by the Three Outlets (Songzi, Taiping, and Ouchi) which are natural waterways. Besides, there are four rivers (Xiangjiang, Zishui, Yuanjiang, and Lishui) flow directly into the lake. The water of the Dongting Lake naturally flows into the Yangtze River at Chenglingji. Characterized by the subtropical monsoon climate, the Dongting Lake is hot and rainy in summer, but warm and dry in winter. The annual mean precipitation and temperature in the Dongting Lake are about 1350 mm and 17.4 °C, respectively. Precipitation in this area exhibits remarkable seasonal variation, the seasonal distribution of discharge in the Dongting Lake is timed simultaneously with precipitation, thus causing frequent flood and drought. The drought has occurred more than 10 times and flood occurred about 8 times in recent decades [28].

Data Preparation
Yueyang, Yuanjiang, and Changde meteorological stations are selected as the representative station of the Dongting Lake. Daily meteorological data from 1961-2019 consisted of the temperature, precipitation, actual evapotranspiration, sunshine duration, relative humidity, wind speed, etc. are obtained from the China Meteorological Data Sharing Service System (http://data.cma.cn). The potential evapotranspiration is calculated by the Penman-Monteith equation [29]. The data of the Dongting Lake is calculated by the Tessellation polygon method [30].

Methods
The total change of the runoff depth includes the changes caused by human activities and climate change (Equation (1)). Moreover, the total change of annual mean runoff depth is the

Data Preparation
Yueyang, Yuanjiang, and Changde meteorological stations are selected as the representative station of the Dongting Lake. Daily meteorological data from 1961-2019 consisted of the temperature, precipitation, actual evapotranspiration, sunshine duration, relative humidity, wind speed, etc. are obtained from the China Meteorological Data Sharing Service System (http://data.cma.cn). The potential evapotranspiration is calculated by the Penman-Monteith equation [29]. The data of the Dongting Lake is calculated by the Tessellation polygon method [30].

Methods
The total change of the runoff depth includes the changes caused by human activities and climate change (Equation (1)). Moreover, the total change of annual mean runoff depth is the difference of the annual mean runoff depth of baseline period and variation period (baseline period is the reference, Equation (2). The relative contribution of climate change and human activities in runoff depth variation can be calculated by Equations (3) and (4).
where ∆Q total is the total change of the annual mean runoff depth (mm); ∆Q variation is the annual runoff depth during variation period (mm); ∆Q baseline is the annual runoff depth during the baseline period (mm); ∆Q climate is the runoff depth change caused by climate change (mm); ∆Q human is the runoff depth change caused by human activities (mm); µ climate and µ human is the contribution rate of the climate change and human activities to runoff depth change.

Simple Linear Regression Method
Simple linear regression (SLR) is a statistical method that establishes the linear relation between runoff and precipitation, which does not consider the physical processes [31]. The required meteorological data of this method is limited, many researchers have applied SLR to estimate the relative contributions of climate change and human activities to changes in runoff or sediment [10,31,32].
The regression equation between runoff depth and precipitation in the baseline period can be expressed as Equation (5). If there is no human activity, the constants (a and b) will not change with time, therefore the reconstructed runoff depth only affected by climate during the variation can be calculated by Equation (6).
where a and b are constant, determined by runoff depth and precipitation in baseline period; Q baseline is annual mean runoff depth in the baseline period (mm); P baseline is annual mean precipitation in the baseline period (mm); Q reconstruct is annual mean runoff depth only effected by the climate in the variation period (mm), and P variation is annual mean precipitation in the variation period (mm). The runoff depth change caused by human activities and climate change can be estimated by Equations (7) and (8).
where ∆Q human refers to the runoff depth changes caused by human activities during the variation period (mm); ∆Q climate refers to the runoff depth changes caused by climate change during the variation period (mm); Q variation is the observed annual mean runoff depth during the variation period (mm); Q reconstruct is annual mean runoff depth during variation period only affected by climate change (mm); the contribute rates of climate change and human activities are calculated by Equations (3) and (4).

Double-Mass Curve Method
The double-mass curve (DMC) is a simple, visual, and practical method, and it is highly applied in the consistency test of long-term hydro-meteorological data [11]. Recently, it has been applied to assess the impacts of climate change and human activities to runoff and sediment changes [10,12,32].
The DMC is composed of cumulative values of two parameters plotted against one another over a certain time span [33]. The theory of the DMC is based on the fact that a plot of the two cumulative quantities during the same period exhibits a straight line so long as the proportionality between the two remains unchanged [11].
The linear relationship between cumulative runoff depth and cumulative precipitation during the baseline period can be expressed as Equation (9), and the reconstructed cumulative runoff depth in the variation period, which is only affected by climate change, is calculated by Equation (10).
where c and d are constant, which can be calculated by  (7) and (8), and contribute rate is calculated by Equations (3) and (4).

Budyko-Based Methods
Given independent variables of precipitation and potential evapotranspiration which are considered as complete climate-controlled variables, so changes of river runoff due to climate change could be described as [34][35][36] where Q is the runoff depth (mm); P is the precipitation (mm); E 0 is the potential evapotranspiration (mm); ∆Q, ∆P, and ∆E 0 are the change of runoff depth, precipitation, and potential evapotranspiration, respectively (mm). The climate elasticity of runoff (ε) can be defined as the ratio of the runoff variation rate to the variation rate of a certain climate factor [13]: where X represents precipitation or potential evapotranspiration in this paper. Therefore, Equation (11) can be written as where ε P and ε E0 are the runoff elasticity to precipitation and potential evapotranspiration, respectively. The elasticity coefficients are calculated by the Budyko-based method. According to the Budyko framework [16], the ratio of actual evapotranspiration to precipitation (E/P, where E is the actual evapotranspiration (mm)) is a function of the aridity index (φ = E 0 /P), which is expressed as Water 2020, 12, 2729 6 of 16 Combined with the water balance equation (Equation (15)), the elasticity coefficients of runoff to precipitation and potential evapotranspiration can be calculated as Equations (16) and (17) [15].
where ∆S is the change in water storage (mm), which is assumed to be 0 over a long time (e.g., five to ten years) [27,37].
Several pieces of research have proposed different functions to describe F. Six popular functions are used in this study ( Table 1). The parameters (m in Fu and ω in Zhang) represent catchment properties which can influence runoff, and these parameters should be calibrated using the data in the baseline period.

Meteorological Data and Hydrological Data Analysis
According to the water balance equation (Equation (15)), the runoff depth can be simplified as the difference between precipitation and actual evapotranspiration [27,28,36]. The Mann-Kendall (MK) test [45,46] is used to analyze the temporal variation of the runoff depth data and detect the change points of runoff depth series [12]. As presented in Figure 2, the mutations occurred around 1984 and 2005. Hence the period of 1961-2019 is divided into three parts. The period of 1961-1983 is the natural time which is regarded as the baseline period (BP). The period of 1984-2004 is variation period 1 (P1), the period of 2005-2019 is variation period 2 (P2). Annual temperature, actual evapotranspiration, precipitation, potential evapotranspiration, and runoff depth series are shown in Figure 3a. Box plots are used to investigate the characteristics of those data in BP, P1, and P2, which include Q1 (the median of the lower half of the sample), Q2 (median of the sample), Q3 (the median of the upper half of the sample), minimum (Q1 − 1.5 × IQR, where IQR = Q3 − Q1), maximum (Q3 + 1.5 × IQR) and outliers (lower than minimum or higher than maximum) (Figure 3b). In addition, the Cv (Coefficient of variation) of these variables are calculated (Table 2). Note: T = temperature; E = actual evapotranspiration; P = precipitation; E 0 = potential evapotranspiration; Q = runoff depth.
K) test [45,46] is used to analyze the temporal variation of the runoff depth data and detect th nge points of runoff depth series [12]. As presented in Figure 2, the mutations occurred aroun 4 and 2005. Hence the period of 1961-2019 is divided into three parts. The period of 1961-1983 i natural time which is regarded as the baseline period (BP). The period of 1984-2004 is variatio iod 1 (P1), the period of 2005-2019 is variation period 2 (P2). Annual temperature, actua potranspiration, precipitation, potential evapotranspiration, and runoff depth series are shown i ure 3a. Box plots are used to investigate the characteristics of those data in BP, P1, and P2, whic lude Q1 (the median of the lower half of the sample), Q2 (median of the sample), Q3 (the media the upper half of the sample), minimum (Q1 − 1.5  IQR, where IQR = Q3 − Q1), maximum (Q1  IQR) and outliers (lower than minimum or higher than maximum) (Figure 3b). In addition, th (Coefficient of variation) of these variables are calculated (Table 2). Water 2020, 9 FOR PEER REVIEW 7

Meteorological Data and Hydrological Data Analysis
According to the water balance equation (Equation (15)), the runoff depth can be simplified as the difference between precipitation and actual evapotranspiration [27,28,36]. The Mann-Kendall (MK) test [45,46] is used to analyze the temporal variation of the runoff depth data and detect the change points of runoff depth series [12]. As presented in Figure 2, the mutations occurred around 1984 and 2005. Hence the period of 1961-2019 is divided into three parts. The period of 1961-1983 is the natural time which is regarded as the baseline period (BP). The period of 1984-2004 is variation period 1 (P1), the period of 2005-2019 is variation period 2 (P2). Annual temperature, actual evapotranspiration, precipitation, potential evapotranspiration, and runoff depth series are shown in Figure 3a. Box plots are used to investigate the characteristics of those data in BP, P1, and P2, which include Q1 (the median of the lower half of the sample), Q2 (median of the sample), Q3 (the median of the upper half of the sample), minimum (Q1 − 1.5  IQR, where IQR = Q3 − Q1), maximum (Q1 + 1.5  IQR) and outliers (lower than minimum or higher than maximum) (Figure 3b). In addition, the Cv (Coefficient of variation) of these variables are calculated ( Table 2).   As shown in Figure 3, the average runoff depth in BP is 565.14 mm. The average runoff depth in P1 is 583.30 mm, which is 18.17 mm higher than that in BP. The average runoff depth in P2 is 440.84 mm, which is 124.30 mm less than that in BP. The runoff depth increased by 3.21% in P1 and decreased by 22.00% in P2 compared with BP. The Cv of runoff depth in P2 is larger than other periods, which means the inter-annual change of runoff depth in P2 is greater (Table 2).
Different from runoff depth, the annual temperature is on the rise as a whole (Figure 3a). In the background of global warming, which began in the 1970s [47], the temperature of the Dongting Lake began to increase slowly in the 1970s and rose rapidly from 1997. The average temperature in BP is 17.04 • C. The average temperature in P1 is 17.41 • C, which is 0.37 • C. higher than that in BP. The average temperature in P2 is 18.06 • C, which is 1.02 • C. higher than that in BP. It increased by 2.17% and 6.03% in P1 and P2 compared with BP, respectively. The temperature in BP is distributed mainly in the low-temperature zone, and that in P2 is distributed mainly in the high-temperature zone (Figure 3b). The inter-annual change of temperature in P1 is greater than the other two periods ( Table 2).
Similar to temperature, there are distinct stepwise increases in actual evapotranspiration over the three periods. The average actual evapotranspiration in BP is 729.36 mm, and that in P1 and P2 are 817.06 mm and 859.00 mm, which are 87.70 mm and 129.64 mm more than that in BP, respectively. The actual evapotranspiration increased by 12.02% and 17.78% in P1 and P2 compared with BP, respectively (Figure 3a). However different from temperature, the actual evapotranspiration distributed mainly in low actual evapotranspiration in all three periods (Figure 3b), and the data in P2 distributed most concentrated, which means the inter-annual change in P2 is smaller than BP and P1 ( Table 2).
The average precipitation in BP is 1294.50 mm. The average precipitation in P1 is 1400.37 mm, which is 105.87 mm more than that in BP, and the average precipitation in P2 is 1299.84 mm, which is 5.34 mm higher than that in BP. It increased by 8.18% and 0.41% in P1 and P2, respectively (Figure 3a). The inter-annual change of precipitation in P2 is greater than the other two periods (Table 2), and the data in P2 distributed mainly in a low rainfall zone (Figure 3b).
The average potential evapotranspiration in BP is 1021.77 mm. The average potential evapotranspiration in P1 is 1011.20 mm, which is 10.57 mm less than that in BP. The average potential evapotranspiration in P2 is 1052.71 mm, which is 30.94 mm more than that in BP. It decreased by 1.03% at first and then increased by 3.03% (Figure 3a). Therefore, the evaporation capacity of the Dongting Lake reduced in P1 and improved in P2 compared with BP. The inter-annual change of potential evapotranspiration in BP is the greatest in three periods ( Table 2). Located in the humid area, the aridity index (φ) of the Dongting Lake is less than 1 in all periods, the energy supply limits actual evapotranspiration [17].

Assessing Contribution of Climate Change and Human Activity to Runoff Alteration
Runoff depth alterations caused by climate change and human activities, calculated with eight methods, as well as the relative contribution rates of human activities and climate change to runoff depth change in P1 and P2 are shown in Table 3. The first group of methods we used is the statistic method including the SLR and the DMC. Both the results of two methods suggest that climate change and human activities have opposite effects on runoff depth change in P1. Climate change causes the runoff depth to increase 70.14-121.51 mm, while human activities make it decrease 51.98-103.35 mm. Thus, the effect of climate change on runoff depth change is greater than that of human activities. As for P2, both the results of two methods indicate that climate change and human activities have opposite effects on runoff depth change. Climate change causes the runoff depth to increase 6.12-22.91 mm, while human activities make it decrease 130.42-147.21 mm. Different from P1, the effect of human activities on runoff depth change is greater than that of climate change. The change of runoff depth is dominated by human activities, while climate change plays a complementary role.
The second group of methods are Budyko-based methods. Different studies propose different methods based on the climate elastic model; six calculation formulations are chosen in this study. Calibrated with the data in BP, m = 2.5 and ω = 1.1. As shown in Table 3, these methods indicate that climate change and human activities influence the runoff depth oppositely in P1. The increase of the runoff depth caused by climate change ranges from 83.81 mm to 105.81 mm, the decrease of the runoff depth caused by human activities ranges from 65.65 mm to 87.64 mm. The effect of climate change on runoff depth change is greater than that of human activities. As for P2, these six methods have the same results that both climate change and human activities make runoff depth decrease and human activities play the dominant role. The decrease caused by climate change ranges from 8.49 mm to 14.34 mm, accounting for 6.83% to 11.53% of total runoff depth decrease. The decrease caused by human activities ranges from 109.96 mm to 115.81 mm, accounting for 88.47-93.17% of total runoff depth decrease. So human activities play a dominant role in runoff depth alteration, while climate change plays a complementary role in P2.
Generally, for P1, all of these methods show that climate change has the opposite effect on human activities on runoff depth change. Climate change causes runoff depth increase while human activities make it decreased. The decline of runoff depth caused by human activities ranges from 51.98 mm to 103.35 mm (average 75.94 mm), the increase of runoff depth caused by climate change ranges from 70.14 mm to 121.51 mm (average 94.11 mm). Climate change plays a dominant role and leads to an increase in runoff depth. As for P2, all of these methods reveal the human activities play a chief part in reducing runoff depth. However, two groups of methods come to the opposite conclusion in the influence of climate change in runoff depth, the SLR and DMC indicate that climate change leads to the runoff depth increasing while the Budyko-based methods show that climate change leads to the runoff depth decrease. As to the SLR and DMC, they consider precipitation as the only climatic factor influences runoff. While the Budyko-based method considered both precipitation and potential evapotranspiration as climatic factors. Hence, in this study, the Budyko-based methods are considered more reliable in partitioning the impact of climate change and human activities on runoff. Based on the second group of methods, the decline of runoff depth caused by human activities accounts for 88.47-93.17% (average 91%), which is the primary reason for runoff depth variation. The SLR and MDC's performances in quantifying the impact of human activities and climate change on runoff alteration are different from region to region [9,31], so both study areas and data should be considered when choosing a method for assessing the effect of climate change and human activities on runoff change.

Influence Factors for Runoff
For P1, climate change makes runoff depth increase while human activities make it decrease, and the runoff depth increased on the whole. The precipitation has increased in P1, which is a completely climatic factor. Actual evapotranspiration also increased in P1, which may be affected by climate change and human activities. However, the effect of increased rainfall is greater. Therefore, precipitation is the main factor that influences the runoff in P1. For P2, both climate change and human activities decrease runoff. The precipitation in P2 increases very slightly, and temperature raise causes actual evapotranspiration increase. The influence of temperature might greater than precipitation, so climate change causes runoff decrease as a whole. At the same time, human activities bring about a runoff decrease. Human activities that change the runoff can be divided into two groups: (1) direct impact such as the construction of reservoirs and water consumption; and (2) indirect impacts such as the activities which may change the land-use or land-cover and thereby change the conditions for runoff generation [48]. Firstly, as presented in Figure 4 [49], the midsize reservoirs whose capacity is between 10 and 100 million m 3 in Hunan province increased greatly in 2005 (up to 276) and 2017 (up to 362), and the water storage of these reservoirs increased to 3.808 billion m 3 in 2018. The large reservoirs whose capacity is more than 100 million m 3 increased gradually from 2001 to 2016 and increased greatly in 2018 (up to 45), and the water storage of large reservoirs increased to 20.388 billion m 3 in 2018. Water consumption has also increased in P2 ( Figure 5) [49,50]. Secondly, land-use or land-cover around the Dongting Lake have changed during 1993-2010: percent coverage of mudflat and meadow decreased, wood coverage rate increased, while the coverage of reed changed little [51]. Although there is little data before 1993, the data mentioned above partly indicated the change of land-use types. Hence land-use or land-cover change may account for the runoff alteration caused by human activities. Besides, the GD and TGD are in the main stream of the Yangtze River upper reaches of the Dongting Lake. However, the mechanisms for these human activities to change the runoff of the lake and how much they each change the runoff still need further research. The underground water also should be considered in further research.

Influence Factors for Runoff
For P1, climate change makes runoff depth increase while human activities make it decrease, and the runoff depth increased on the whole. The precipitation has increased in P1, which is a completely climatic factor. Actual evapotranspiration also increased in P1, which may be affected by climate change and human activities. However, the effect of increased rainfall is greater. Therefore, precipitation is the main factor that influences the runoff in P1. For P2, both climate change and human activities decrease runoff. The precipitation in P2 increases very slightly, and temperature raise causes actual evapotranspiration increase. The influence of temperature might greater than precipitation, so climate change causes runoff decrease as a whole. At the same time, human activities bring about a runoff decrease. Human activities that change the runoff can be divided into two groups: (1) direct impact such as the construction of reservoirs and water consumption; and (2) indirect impacts such as the activities which may change the land-use or land-cover and thereby change the conditions for runoff generation [48]. Firstly, as presented in Figure 4 [49], the midsize reservoirs whose capacity is between 10 and 100 million m 3 in Hunan province increased greatly in 2005 (up to 276) and 2017 (up to 362), and the water storage of these reservoirs increased to 3.808 billion m 3 in 2018. The large reservoirs whose capacity is more than 100 million m 3 increased gradually from 2001 to 2016 and increased greatly in 2018 (up to 45), and the water storage of large reservoirs increased to 20.388 billion m 3 in 2018. Water consumption has also increased in P2 ( Figure  5) [49,50]. Secondly, land-use or land-cover around the Dongting Lake have changed during 1993-2010: percent coverage of mudflat and meadow decreased, wood coverage rate increased, while the coverage of reed changed little [51]. Although there is little data before 1993, the data mentioned above partly indicated the change of land-use types. Hence land-use or land-cover change may account for the runoff alteration caused by human activities. Besides, the GD and TGD are in the main stream of the Yangtze River upper reaches of the Dongting Lake. However, the mechanisms for these human activities to change the runoff of the lake and how much they each change the runoff still need further research. The underground water also should be considered in further research.   [7]. Yanhe River, a first-order tributary of the Yellow River, belongs to a semi-arid area in China, while the Dongting Lake belongs to the humid area. Although the study period of the  [7]. Yanhe River, a first-order tributary of the Yellow River, belongs to a semi-arid area in China, while the Dongting Lake belongs to the humid area. Although the study period of the Yanhe River is not the same with this paper, the conclusions of these two studies might concur with the regularity found in the study of Li et al. [36] which indicated that influences of human activities were found considerable in the grassland covered semi-humid area, while climate change functioned remarkably in the semi-arid area with sparse vegetation.

Sensitivity Analysis
As shown in Equation (13), the change in runoff caused by climate change is a function of change in precipitation and change in potential evapotranspiration. In Budyko-based methods, ε P and ε E0 in that function can be regarded as sensitivity coefficients. Therefore, the sensitivity coefficients of the two periods are calculated with six formulations in Table 1, the results are shown in Table 4. In P1, ε P ranges from 1.76 to 2.19 (average 1.95), which means 1% change in precipitation results in 1.76% to 2.19% (average 1.95%) change in runoff; ε E0 ranges from −1.19 to −0.76 (average −0.95), which means 1% increase in precipitation results in 0.76% to 1.19% (average −0.95%) decrease in runoff. While in P2, ε P ranges from 1.80 to 2.24 (average 1.99), which means 1% change in precipitation results in 1.80% to 2.24% (average 1.99%) change in runoff; ε E0 ranges from −1.24 to −0.80 (average −0.99), which means 1% increase in precipitation results in 0.80% to 1.24% (average −0.99%) decrease in runoff. Although the values of the climate elasticity calculated by various methods are different, all climate elasticity in P2 is larger than that in P1, which reveals runoff are more sensitive to climate change in P2 compared with P1. Besides, no matter in P1 or P2, runoff is more sensitive to change in precipitation than the change in potential evapotranspiration. In other words, the effect of precipitation to runoff is greater than the evaporation capacity. The average sensitivity coefficients in two periods are used to draw contours with which runoff change can be calculated when there is any change in precipitation or potential evapotranspiration ( Figure 6). There is a slight difference between the two periods. With Figure 6b, the change of runoff in the future can be rough estimated when the change of precipitation and potential evapotranspiration is known. Moreover, ε P and ε E0 in the Yanhe River calculated by the same method are larger than those in the Dongting Lake [7]. Combined with the analysis of 4.1, it can be speculated that the change of the underlying surface might account for that runoff being more sensitive to climate change in P2 compared with P1. Therefore, climate change and human activities are not fully independent factors, which need further research. the two periods. With Figure 6b, the change of runoff in the future can be rough estimated when the change of precipitation and potential evapotranspiration is known. Moreover, P and E0 in the Yanhe River calculated by the same method are larger than those in the Dongting Lake [7]. Combined with the analysis of 4.1, it can be speculated that the change of the underlying surface might account for that runoff being more sensitive to climate change in P2 compared with P1. Therefore, climate change and human activities are not fully independent factors, which need further research.

Conclusions
In this study, the meteorological and hydrological data from 1961 to 2019 of the Dongting Lake are analyzed; eight methods are used to decompose the effects of climate change and human activities on runoff alteration, and the factors for runoff changing are analyzed; sensitivity of the runoff to changes in climatic variables is analyzed. The key findings are summarized as follows: (1) According to the MK test, the mutation of runoff occurred in 1984 and 2005, so the study period is divided into three periods: 1961-1983 (BP), 1984-2004 (P1), and 2005-2019 (P2). Precipitation increases in P1 and decreases in P2 compared with BP, while potential evapotranspiration decreases in P1 and increases in P2 compared with BP. Temperature and actual evapotranspiration have the same changing tendency, they increase with time in three periods. Runoff increases in P1 and decreases in P2.
(2) The results of eight methods indicate that climate change causes runoff increases 386% to 669% (average 518%) while human activities make it decrease 286% to 569% (average 418%) in P1. Climate change plays a prime role in altering runoff in P1. Based on the Budyko-based methods, the decline of runoff caused by human activities account for 88% to 93% (average 91%), while the decline of runoff caused by climate change accounts for 7% to 12% in P2. Human activities are the main factors for runoff change in P2.
(3) For P1, a climatic factor especially precipitation is the main factor that influences the runoff. For P2, human activities such as reservoirs construction, water consumption, and land-use (landcover) change are the main factors that influence the runoff in the Dongting Lake. The mechanisms for human activities to alter runoff and how much they each change the runoff need further research.
(4) According to the sensitivity analysis, runoff in the Dongting Lake is more sensitive to climate change in P2 compared with P1. Besides, no matter in P1 or P2, runoff is more sensitive to change in precipitation than a change in potential evapotranspiration. Combined with climate forecast, the results of sensitivity analysis can be used to estimate runoff change caused by climate change in the future.

Conclusions
In this study, the meteorological and hydrological data from 1961 to 2019 of the Dongting Lake are analyzed; eight methods are used to decompose the effects of climate change and human activities on runoff alteration, and the factors for runoff changing are analyzed; sensitivity of the runoff to changes in climatic variables is analyzed. The key findings are summarized as follows: (1) According to the MK test, the mutation of runoff occurred in 1984 and 2005, so the study period is divided into three periods: 1961-1983 (BP), 1984-2004 (P1), and 2005-2019 (P2). Precipitation increases in P1 and decreases in P2 compared with BP, while potential evapotranspiration decreases in P1 and increases in P2 compared with BP. Temperature and actual evapotranspiration have the same changing tendency, they increase with time in three periods. Runoff increases in P1 and decreases in P2.
(2) The results of eight methods indicate that climate change causes runoff increases 386% to 669% (average 518%) while human activities make it decrease 286% to 569% (average 418%) in P1. Climate change plays a prime role in altering runoff in P1. Based on the Budyko-based methods, the decline of runoff caused by human activities account for 88% to 93% (average 91%), while the decline of runoff caused by climate change accounts for 7% to 12% in P2. Human activities are the main factors for runoff change in P2.
(3) For P1, a climatic factor especially precipitation is the main factor that influences the runoff. For P2, human activities such as reservoirs construction, water consumption, and land-use (land-cover) change are the main factors that influence the runoff in the Dongting Lake. The mechanisms for human activities to alter runoff and how much they each change the runoff need further research.
(4) According to the sensitivity analysis, runoff in the Dongting Lake is more sensitive to climate change in P2 compared with P1. Besides, no matter in P1 or P2, runoff is more sensitive to change in precipitation than a change in potential evapotranspiration. Combined with climate forecast, the results of sensitivity analysis can be used to estimate runoff change caused by climate change in the future.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.