Accessing the Difference in the Climate Elasticity of Runoff across the Poyang Lake Basin , China

Understanding the effects of climate and catchment properties’ changes on water yield is a challenging component in assessments of future water resources. Here, we spatially applied the water-energy balance equation, based on the widely-used Budyko framework, to quantify the temporal and spatial differences of the climate elasticity of runoff in the Poyang Lake Basin (PYLB), highlighting the influence of the catchment properties’ parameter n variation on the climate elasticity and runoff prediction. By using Sen’s slope and the Mann–Kendall method, we found that, for the whole study period (1960–2010), annual temperature in PYLB significantly increased at a rate of 1.44% per decade. Basin-wide wind speed and net radiation had been declining at 0.17 m/s and 46.30 MJ/m2 per decade. No significant trend was detected in precipitation and relative humidity. The moving average method was applied to evaluate the temporal pattern of n. The results showed that the calibrated catchment properties’ parameter and the derived elasticities were not constant during the past 50 years. We found that in most sub-basins, the n values increased during 1970–1980, followed by a decreasing trend in the period from 1980 to 1990, whereas the n value in Fuhe sub-basin kept increasing for the almost the whole study period. In addition, the climate elasticity is highly correlated with the n value, indicating that the catchment properties’ parameter was the dominant factor influencing climate elasticity in PYLB in the past 50 years. We also attempted to predict the runoff trend with the consideration of trends in n. However, in some sub-basins, there were still considerable differences between the predicted runoff trend and the observed one. The method used here to evaluate the temporal pattern of n should be an extension of the existing literature and will provide a better understanding of elasticity in the regional hydrological cycle.


Introduction
Recently, the climate has changed considerably [1], and many of the observed changes are unprecedented over decades to millennia [2].Hence, understanding the change in regional water availability and hydrological cycle under a warming climate is of great importance for sustainable development of human society, ecosystem services, as well as the regional water resources' management [3][4][5][6].Planning for potential future changes in water yield requires not only an understanding of possible climate change, but also understandings of how and to what extent a catchment's hydrology might respond to such change.However, estimations of future climatic conditions are currently poorly suited to hydrological prediction [7], while climate elasticity is considered as an important indicator for quantifying how the runoff would change in response to climate change [5,8,9].The climate elasticity was defined by Fu [8] as the proportion of the change in runoff to the change in climate variables.Following that, various methods have been developed to compute the climate elasticity.Sankarasubramanian et al. [6] well classified these methods into five categories: (1) varying the input of a calibrated conceptual watershed model to observe the resulting runoff change; (2) analytically derivation; (3) fitting multivariate regional hydrological models; (4) empirical estimation based on historical data; and (5) applying multivariate statistical methods.The model-based approach might be physically sound, but requires major efforts on model calibration and can lead to remarkably different results because of the uncertainty in model structure and parameter estimation [10][11][12][13].Moreover, given that climatic conditions are not stationary [14], such hindcasting approaches are problematic for prediction.In contrast, the Budyko hypothesis ingeniously couples hydrological cycle and terrestrial energy balance [3,[15][16][17], which contains not only the regional climate conditions, but also the parameter reflecting underlying characteristics [18][19][20][21][22]. Based on the Budyko hypothesis, more generic methods of elasticity have been developed [6,16], amongst which, the analytical derivation approach, depending less on historical data, is clear in theory and, hence, has captured much attention [3,4,9,16,18].
In light of the derivation approach, a unified framework for assessing the elasticity of runoff to changes in average precipitation P and potential evaporation E 0 has been recently established [3,18,23,24].This framework renders climate elasticity analyses highly accessible, such that it can be applied wherever long-term average precipitation and other climate data exist.Additionally, the framework can be used to assess different hydrological responses likely across different regions of a catchment or basin, which is important for management decisions relating to possible future hydrological conditions.Consequently, it has been widely applied in many studies to evaluate the impact of climate change on runoff [3,4,18,20,23].Although the simple framework is powerful for assessing the impacts of average precipitation P and potential evaporation E 0 on runoff, the variations in potential evaporation are hardly measured due to its inconvenience in calculation, which requires many data [25].Is it possible to separate the influences of potential evaporation into the influences of more accessible climate variables (e.g., temperature, relative humidity, wind speed, etc.)?Since these climate variables' influences cannot be derived simply based on the water-energy balance equation introduced by previous studies (the M-C-Y equation) [16,26,27], more efforts are needed to figure it out.Yang and Yang [3] recently introduced an analytical solution to this problem by combining the original M-C-Y equation with the Penman equation [28], which made it possible to relate the elasticity of runoff to net radiation, temperature, wind speed and relative humidity, providing a guide for separating the effects of different climate variables on runoff.
Due to its high accessibility, the above-mentioned framework has been widely used.Many previous studies, however, attempted to derive the climate elasticity based on a single (maybe slightly different in forms) Budyko-type curve [5,6,9,23,24].In accordance with the derived climate elasticity by Yang and Yang [3], we could expect that not only the changes in climate conditions, but also the variation in catchment properties will impact the derived elasticity.To our knowledge, Zheng et al. [23] appears to be the first to demonstrate the temporal variation of climate elasticity using a moving average method, but he did not pay attention to the temporal variation of catchment properties.The inclusion of changes in catchment properties over time seems to have been first considered by Yang et al. [21], but they assumed that the catchment properties would remain constant in that study, and no further calculation of the associated elasticity was made.Following that, Roderick and Farquhar [18] first introduced a qualitative way to explain why the catchment properties would vary over time, but their work also did not consider how the catchment properties' parameter might change over time.Recently, more and more attention has been paid to the variation of catchment properties and its impact on the climate elasticity.As no definitive process-based understanding of what determines the catchment properties' parameter is currently available, it is hard to give an a priori estimate of the parameter [4,29].Consequently, some studies tried to relate the catchment properties to other climate factors, such as climate seasonality and mean storm depth [29], while some attempted to link it to vegetation coverage or effective rotting depth and plant root characteristics [20,29].The key papers related to the Budyko-based framework are summarized in Table 1.As the key factors influencing the catchment properties' parameter remain vague, the evaluation of the climate elasticity variation with the consideration of the change in catchment properties should be an extension of the existing literature and will provide a better understanding of elasticity in the regional hydrological cycle.
Table 1.Summary of the relevant literature on the climate elasticity of runoff studies based on the Budyko hypothesis in relation to the aims of the present study.In the second column, the consideration of the catchment properties' parameter is identified by the code: (a) consider the spatial variation of the parameter over different catchments; (b) consider the temporal variation of the parameter.We add the current study for completeness.The catchment properties' parameter has a significant regional pattern and represents the land surface conditions.The Poyang Lake Basin (PYLB) provides a good example where hydrological characteristics vary across space and where management strategies need to be carefully targeted.The objectives of this paper are to (1) evaluate the trend in climatic and hydrologic variables across PYLB during the past 50 years; (2) figure out how the catchment properties' parameter would change over time in PYLB during the past 50 years; (3) determine the spatial and temporal pattern of climate elasticity of runoff across PYLB; and (4) evaluate the impact of catchment properties' parameter variation on the derived elasticities and runoff prediction.

Overview of the Poyang Lake Basin
Located on the middle and lower reaches of the Yangtze River on the south shore, Poyang Lake (Figure 1) is the largest freshwater lake and an important ecological function protected area in China [32].The lake has a drainage basin area of 162,225 km 2 , covering about 97% of Jiangxi province [33].The topography of the Poyang Lake Basin (PYLB) varies from highly mountainous regions (maximum elevation of 2200 m above sea level) to alluvial plains in lower reaches of the primary water resources [34,35].Poyang Lake exchanges water with the Yangtze River through a narrow outlet located in Hukou [35], while receiving water from the five sub-drainage basins of Ganjiang River, Fuhe River, Xinjiang River, Raohe River and Xiushui River [36,37].Headwater of these rivers is located in boundaries of the east, south and west of Jiangxi Province that are surrounded by high mountains.The stream gradient decreases as these rivers flow onto the relatively flat region surrounding Poyang Lake.The flat central area of the Poyang Lake Basin is an important region for farming (major crops are rice and oilseed rape) [35], not only in Jiangxi province, but also in China [34].The PYLB belongs to a subtropical wet climate zone with an annual mean precipitation of 1680 mm and annual mean temperature of 17.5 • C. The hydrological processes in PYLB is complex [36,38], leading to dramatic seasonal water-level fluctuations of Poyang Lake.In normal years, Poyang Lake can expand to a large water surface of 4000 km 2 with a volume of 3.2 × 10 10 m 3 and will shrink to little more than a river during the dry season, shaping the unique landscape that has been described as "flooding like the sea, drying like thread" [34,39].

Data Preparation
Observed discharges (daily runoff) of seven gauging station (Figure 1) on the lower reaches of the five rivers were obtained from Management Bureau of the Yangtze River (MBYR) catchment for the period from 1960 to 2010.For sub-basins of the Ganjiang, Fuhe and Xinjiang, only one gauging station at each main stream was selected to represent basin runoff, contributing more than 90% of the total inflow of the five rivers; while for sub-basins of the Xiushui and Raohe, gauging stations at tributaries (Wanjiabu station for the Xiushui River, Hushan and Dufengken stations for the Raohe River) were included to account for the runoff contributions to the main stream, contributing less than 10% of the total inflow of the five rivers.The basic features of these gauging stations are listed in Table 2.The total drainage area of these gauging stations is 137,143 km 2 , leaving an area of 25,082 km 2 (15.5% of the whole basin area, including the lake surface) that is not gauged.The meteorological data from 14 National Meteorological Observatory (NMO) stations inside the PYLB were obtained from the National Climate Centre of China Meteorological Administration (CMA), including daily observations of air temperature, wind speed, relative humidity, sunshine hours and absolute vapor pressure, among others.The period of records used here is 1960-2010.Based on the meteorological data, potential evapotranspiration (PET or E 0 ) of each weather station was then estimated following the Penman-Monteith method recommended by the Food and Agriculture Organization of the United Nations (FAO), and the computation of all of the data required for the potential evaporation followed the method and procedure given in Chapter 3 of FAO Paper 56 [40].
All of the daily runoff records and climate data provided by MBYR and CMA had gone through a standard quality control process before delivery [34,41,42].No missing data were found in the variables used in this study, except for daylight hours.The missing daylight hours data were further interpolated using data from the nearest three meteorological stations.
As the Penman-Monteith equation [40] is highly non-linear, consequently, instead of applying the Penman-Monteith equation [40] with basin average variables, we calculated the daily E 0 for each weather station using daily hydro-meteorological records.Following that, all of the hydro-meteorological records, as well as the calculated E 0 were aggregated from daily to monthly and to yearly.In addition, the large degree of variation in topography and the uneven distribution of weather stations across the PYLB are noticeable and should be take into consideration when estimating the basin average data.Consequently, Ye et al. [34] suggested that an area-based weighing method (i.e., Thiessen polygon method [43]) should be used in calculating the average data.In the present study, we followed the same procedure used by Ye et al. [34] when calculating the average value of climate variables and potential evapotranspiration for the whole PYLB, as well as the five sub-basins.

Trend Analysis for Hydro-Meteorological Variables
In the present study, two non-parametric methods (Mann-Kendall test and Sen's slope estimator) were used to detect the trend of climatic variables in PYLB, using temporally-and spatially-aggregated data, namely annual data averaged over five sub-basins.

Mann-Kendall Test
The Mann-Kendall test statistics S [44] is calculated as: where n is the number of data points.x i and x j are the data values in times series i and j (j > i), respectively, and sgn(x j − x i ) is the sign function defined as: The variance is calculated as: where n is the number of data points, m is the number of tied groups and t i denotes the number of ties of extent i.A tied group is a set of sample data having the same value.In cases where the sample size n > 10, the standard normal test statistic Z S is computed as follows: Positive values of Z S indicate an increasing trend, while negative Z S values mean a decreasing trend.Testing the trend is done at the specific α significance level.When |Z S | > Z 1−α/2 , the null hypothesis is rejected, and a significant trend exists in the time series.Z 1−α/2 can be obtained from the standard normal distribution table.In this study, significance levels α = 0.05 and α = 0.01 were used.At the 5% significance level, the null hypothesis of no trend is rejected if |Z S | > 1.96 and rejected if |Z S | > 2.56 at the 1% significance level.

Sen's Slope Estimator
Sen [45] developed the non-parametric procedure for evaluating the slope of trend in the sample of N pairs of data: where x j and x k are the data values at times j and k (j > k), respectively.dt is the desired time interval.
If there is only one datum in each time period, then N = n(n−1)

2
, where n is the number of time periods.Otherwise, N < n(n−1)

2
, where n is the total number of observations.The N values of Q i are sorted from smallest to the largest, and the median of the slope or the Sen's slope estimator is then computed as: where the sign of Q med reflect the data trend pattern, while its value indicates the steepness of the trend.Once the slope Q med has been calculated, the regression line can then been determined from the sample points by setting the intercept b to be the median of the values of y i − Q med × x i [46].

Calibration of the Catchment Properties Parameter
Budyko [15] quantified mean annual water balance in terms of the ratio of mean annual evaporation to mean annual precipitation, E/P, demonstrating that E/Pis determined, to first order, by the ratio of mean annual potential evaporation to mean annual precipitation, E 0 /P (the so-called climatic dryness index), a measure of climatic aridity.Many attempts have been made to formulate the mean annual water-energy balance.Bagrov [47] first derived the analytical equation for the mean annual water-energy balance by introducing dE/dP = 1 − (E/E 0 ) n .A modification of that formula was derived by Mezentsev [27] and was expressed as dE/dP = [1 − (E/E 0 ) n ] (n+1/n) .More efforts have been made for that equation, and it was finally expressed in a generalized form (the M-C-Y equation) [16,26,27]: where E 0 is the potential evaporation, E is the actual evaporation P is the precipitation and n (dimensionless) is a catchment-related parameter, which modifies the partitioning of P between E and Q.
In terms of Equation (7), n is the main model parameter that needs to be calibrated.The values of n are typically in the range of 0.6-3.6, and most are in a smaller range from 1.5 to 2.6 [3,16,18,26].A default value of n is recommended by Choudhury [26] as 1.8.In this study, the inhomogeneous characteristics of each sub-basin have been taken into consideration when estimating n.Apart from the spatial pattern of n, we also aim to evaluate the temporal variation of n and its subsequent effects on climate elasticity estimation.In fact, the n value can be estimated according to Equation ( 7) whenever the annual mean values of climatic and hydrologic variables are available.Note that the Budyko framework used here assumes steady state conditions and therefore requires a time scale whereby changes in catchment storage are small relative to the magnitude of fluxes (i.e., P, E, R) [4].In practice, this requires averages over at least one year [18], and previous studies often applied the Budyko framework with averages over 10 years [23,34], 30 years [18] and even longer periods [3].In the present study, we assume that it would be appropriate to apply the Budyko framework based on the 10-year average.In addition, with that assumption in mind, we further assume that the Budyko framework can be further applied with 10-year moving averages.Following the assumption, the annual values of the related climatic and hydrologic variables were first calculated using the moving average method.The moving window length used in the present study was 10-year.After that, we estimated the parameter n based on R (mean annual runoff) and P (mean annual precipitation) for each sub-basin within each moving window in order to obtain exact agreement with observations.The calibration was done using a Newton-Raphson method with a desired accuracy of 0.001 [48].The correlation coefficient and mean absolute error (MAE) between the simulated and observed runoff within each moving window were further calculated to evaluate the calibration results of n [34].

Derivation of Climate Elasticity
Climate elasticity was derived following the same framework developed by Yang and Yang [3], except for the potential evaporation equation used.As mentioned before, the potential evaporation was estimated following the Penman-Monteith method recommended by FAO [40].By combining Equation (7), which can be denoted as f (E0, P, n), and the Penman-Monteith equation, which is denoted as v(Rn, T, U, RH, ...), the changes in E due to changes in climate and the catchment properties can be derived.To first order, the change in E is: When evaluating the trend in annual evaporation or runoff for a long period (10 years in this study), the water balance equation can be simplified as R = P − E, which leads to the differential form dR = dP − dE.By direct substitution into Equation (8), we obtain: and the relative change in R is: The terms in square brackets can be called the sensitivity coefficients and are analytical expressions formally equivalent to the "elasticity" concept [3,18,49].The calculation was done using the R statistical software 3.3.2[50].

Runoff Prediction with Climate Elasticity
By combining the trend results (Section 3.1) and the derived elasticities, the predicted trend in runoff can be calculated.As mentioned before, most previous studies often predict the runoff trend using climate elasticities with no consideration of changes in catchment properties.In other words, dn is treated as zero, and the trend in runoff can be denoted as: However, it is problematic if the temporal variation in n cannot be neglected.Hence, in addition to Equation (11), we also predicted the runoff trend with consideration of n variation using Equation (10).Following that, the absolute differences between the two predicted runoff trends, namely trends from Equations ( 10) and (11), and the observed trend (Section 3.1) were calculated to evaluate the influence of n on runoff prediction.

Trend in Climatic and Hydrologic Variables
Following Section 2.2, the average values of the climatic and hydrologic variables across the PYLB were calculated.The results are listed in Table 3.As shown in Table 3 for Xiushui, Xinjiang, Fuhe, Ganjiang and Raohe, respectively.The maximum value of R is in Xinjiang sub-basin, with a value of 1151.8 mm•year −1 , followed by those of Raohe (1010.0mm•year −1 ), Xiushui (938.4 mm•year −1 ), Ganjiang (844.0 mm•year −1 ) and Fuhe (778.9 mm•year −1 ), respectively.No sufficient differences in other climate variables can be found among the sub-basins in PYLB. Figure 2 presents the trends of the annual climatic and hydrologic variable series.The change rate is defined as the ratio of the trend slope to the mean annual value.Generally, the PYLB had been warming during the study period, with annual temperature significantly increasing at a rate of 1.44% (i.e., 0.19 • C) per decade.In addition, basin-wide wind speed (U) and net radiation (R n ) had been declining at 0.17 m/s and 46.30MJ/m 2 per decade.No significant trend was detected in precipitation and relative humidity (except for Raohe sub-basin).Besides the trend of climatic variables, the runoff trend in the PYLB, as well as the five sub-basins was also calculated, and its results are presented in Figure 2. As presented in Figure 2, no significant trend can be detected in runoff across the PYLB over the whole study period due to its huge inter-annual fluctuations.The cell without asterisks denotes that the trend is not statistically significant at both the 1% and 5% significance levels.

Temporal and Spatial Pattern of the Catchment Properties' Parameter
Figure 3 shows the temporal variation of the catchment properties' parameters for each sub-basin, as well as the whole PYLB.As shown in Figure 3, on a 10-year moving window scale, the ranges of the n value are 0.85-1.65,0.86-1.16,1.56-4.77,1.13-1.49and 1.01-1.56for Xiushui, Xinjiang, Fuhe, Ganjiang and Raohe sub-basin, respectively.The basin-wide n value is between 1.22 and 1.55.In sub-basins except for Fuhe, the n values increased during 1970-1980, followed by a sharp decrease in period from 1980 to 1990 and ended with a slight upward trend in the 2000s.However in Fuhe sub-basin, the n value kept increasing for the period from 1960 to the 2000s, which almost persists for the whole study period.In order to evaluate the results of the n calibration, the correlation coefficient and the mean absolute error (MAE) between the simulated and observed runoff within each 10-year moving window were further calculated and summarized in Figure 4.As shown in Figure 4, despite of the great fluctuations, most of the correlation coefficient values are greater than 0.7 with small mean absolute error (MAE), indicating that the simulated results are acceptable.1960 1970 1980 1990 2000 20101960 1970 1980 1990 2000 20101960 1970 1980 1990

Climatic Elasticity in PYLB
After estimating the parameter n, we calculated the climate elasticity based on the annual data for the period from 1960 to 2010.Following the framework developed by Yang and Yang [3], the climate elasticity can be estimated according to the mean annual climatic variables.The mean elasticity results for each sub-basin, as well as the whole PYLB for the entire study period are listed in Table 4.As shown in Table 4, the value of precipitation elasticity (ε P ) is much larger than other climatic variable elasticities, indicating that runoff is more sensitive to precipitation than to other climatic variables (i.e., radiation, temperature, wind speed and humidity) across the PYLB.Specifically, a 1% P increase will result in A 1.67% runoff increase; a 1% increase in R n and T will lead to A 0.50% and A 0.22% decrease in runoff, respectively; a 1% T increase will induce A 0.22% decrease in runoff; and RH appears to have little influence on runoff.Despite the climate variables, the catchment properties' parameter also influences runoff; a 1% increase in n will lead to A 0.27% decrease in runoff.Additionally, the climate elasticities in Fuhe sub-basin ARE relative higher than those in other sub-basins, and no sufficient differences are found among the other four sub-basins.The temporal variation of runoff elasticities in PYLB is shown in Figure 5.As expected, the climate elasticities across the PYLB keep changing over time.However, the variation shows no monotonic trend except for n in Fuhe sub-basin.The dominant pattern to emerge from the analysis of runoff elasticities in Figure 5 is that, in the PYLB, the temporal variation of runoff elasticities to climate variables turns out to be a W-shaped curve, which is to say the elasticities decreased in the beginning of the period and then comes in a small period of recovery followed by another decrease and an eventual long-term upward trend.In addition, the temporal pattern of runoff elasticities to the catchment properties' parameter in the PYLB is different; specifically, the catchment properties' parameter elasticity ( n ) values in sub-basins except for Fuhe fluctuate for the period from 1960 to 2000 with a slight downward trend followed by a sharp upward trend in the end of the study period.In contrast, a sharp dip in n is found for the period from 1960 to 2000 in Fuhe sub-basin and ends with a slight upward trend.

Influence of the Catchment Properties' Parameter on Climate Elasticities and Runoff Prediction
In Sections 4.2 and 4.3, we evaluated the temporal variation of the catchment properties' parameter and the derived climate elasticities, and the temporal pattern of the two were similar, which indicates a relationship between n and climate elasticities.The relationships between n and derived elasticities are shown in Figure 6.As expected, the catchment properties' parameter is significantly correlated with the derived elasticities.Specifically, precipitation elasticity, wind speed elasticity and relative humidity elasticity showed significant positive correlations with n, whereas temperature elasticity and net radiation elasticity were significantly negatively correlated with catchment properties parameter.Figure 7 shows the differences of two (i.e., dn = 0 and dn calculated from the moving window method) predicted runoff trends compared to the observed trend according to Section 4.1.As shown in Figure 7, the predicted runoff trend would get closer to the observed one when taking n variation into consideration in Xinjiang, Fuhe and the whole PYLB.However, in sub-basins like Xiushui, Ganjiang and Roahe, the prediction error could become larger when n was considered.

Potential Factors Influencing the Catchment Properties' Parameter
In this study, with the assumption that the Budyko framework could be applied with moving averages, we estimated the temporal variation of the catchment properties' parameter n in PYLB, as well as the five sub-basins using a 10-year moving average method.Our results showed that the parameter n in Equation ( 7) ranges from 1.0 to 2.0 for most sub-basins in PYLB, consistent with previous findings for the same region [16,51,52]; whereas the n value in Fuhe sub-basin was notably higher.A comparable result was given in Ye et al. [34], who found that the value of the catchment properties' parameter in Fuhe sub-basin was three times as much as n values in other sub-basins in PYLB.
Among the four factors above, vegetation characteristics have always been regarded as the the important factors for n variation in many previous studies [51,52,55], and it is also the case in the PYLB.Due to the rapid development of urbanization, industrialization and agricultural cultivation [56], PYLB, like many other regions in China, has undergone intensive human activities since 1950, including afforestation and deforestation, land reclamation, river regulation, agriculture intensification and extensive infrastructure construction, which may have exerted considerable impacts on catchment hydrology [34].Guo et al. [37] applied the SWAT model in PYLB and demonstrated that when agriculture land changed into forest, accounting for up to 23.3% of the catchment area, it would result in a decrease of annual runoff by up to 3.2%.Figure 8 shows the change of forest coverage and soil erosion area in PYLB for the period from 1960 to 2010 [57].Remarkably, the forest coverage in PYLB increased dramatically since the 1970s.In addition, other studies also demonstrated that the Normalized Difference Vegetation Index (NDVI) in PYLB presented an increasing trend since the 1980s [35,58,59].Both historical and remote sensing evidence highlight the increase of vegetation coverage in PYLB, which may play an important role in reducing runoff.The causes of vegetation change include both climate and human factors, such as: (1) the "Grain-for-Green" revegetation project in not only PYLB, but also the entirety of China [60,61]; (2) irrigation and fertilization [59]; and (3) the impact on vegetation due to increasing CO 2 concentration from human activities [62][63][64].
A previous study [34] also pointed out the importance of soil erosion in the catchment.As shown in Figure 8, an upward trend can be found in the soil erosion area in PYLB during 1960-1990, followed by a dramatic decreasing trend from 1990-2010.Erosion will reduce the water-holding capacity because of rapid water runoff [65], which may increase the amount of runoff in the catchment and, thus, reduce the n value.
Despite the vegetation coverage and soil erosion area, extensive water utilization is also considered as the main factor that reduces the runoff in PYLB.As demonstrated by Ye et al. [34], the volume of water storage in PYLB increased steadily from 1970 to 2010, especially after the 1990s.The reasons for this are the rapid changes in society and social life, population, economic forces and technological development, which increased the demand and the ability of water supply for industrial and domestic consumptions.
In summary, the increase of n values in PYLB during the period from 1970 to 1980 could possible be the result of the increased water storage volume due to the rapid development in the catchment; whereas the dramatic increase in soil erosion area in 1980-1990 should enhance the runoff generation in the catchment and could possibly account for the decline in n for that period.Finally, the upward trend in n in the 2000s may be the sum effects of decreased erosion area, increased vegetation coverage and water storage volume.In addition, the different changing patterns in n values in Fuhe sub-basin since the 1990s (Figure 3) as compared with other sub-basins may be the result of high water-use efficiency.In Fuhe sub-basin, a previous study [34] suggested that it is the only catchment where the reduction in runoff should be mainly attributed to the intensive human activities.The biggest irrigation farmland, as well as the irrigation systems of Jiangxi province are also located in the middle and lower reach of Fuhe, which will notably increase the water utilization and directly decrease runoff.Hence, in the Fuhe sub-basin, a relatively lower runoff Q, which means a higher E (calculated as R = P − E) when P varies little among all of the sub-basins, results in the higher estimated value of n.In addition, the intensive human activities may possibly contribute to the increase of n in Fuhe sub-basin.

Uncertainty and the Limitations of This Study
In the present study, we assumed that the long-term water balance of a catchment will follow the Budyko curve if climate change occurs and further estimated the temporal variation of climate elasticities across the PYLB.Remarkable, as shown in Equation ( 7), the catchment water balance will follow one Budyko curve only when n keeps constant.Hence, to evaluate the temporal variation of climate elasticities across the PYLB, we treated parameter n as a changeable constant, that is the parameter will keep constant in each relatively shorter sub-period (five or 10 years) and will change as sub-periods change.Following the simple assumption, the temporal variation of climate elasticities can be estimated.
However, when applying the n trend results in runoff prediction, the results were not as expected.There were still considerable differences between calculated and observed trend in runoff.The possible causes are the following: (1) some climatic variables or runoff itself have relatively large variabilities, but no significant trends.(2) under the framework developed by Yang and Yang [3], when calculating the runoff trend, the elasticity itself should keep constant, which is not the case in reality.This highlights a shortcoming in the analyses presented here and of the framework as a whole, which is that values of n are not easily estimated a priori for a given area.In practical, the values of n have to be either estimated or, if runoff data exist, fitted to observational data, thus limiting the framework to gauged catchments and to historical conditions.Additionally, the framework used in this study only assesses the changes in steady state conditions.In other words, when the effects of a 10% change in P, for example, have been discussed, it is assumed that such a change has no bearing on soil storage dynamics.
In addition, other possible influences of data error on the results are yet to be investigated.Meteorological data from 14 weather stations in this study area (PYLB) might not be sufficient coverage for such a lager-scale catchment.The influence of human activities on the estimation of mean annual E is not considered.Although the calibrated catchment properties' parameter n in this study well reflects the average catchment condition in most periods, these could still affect the final results to some extent.
Consequently, we should further make n as a function of climate, vegetation and other potential variables, i.e., n = f (C, V, ...), which should be able to quantify the temporal variation of n in a finer time scale (e.g., yearly) and make the n a priori estimate.By combing this function with the original Budyko hypothesis, the climate impact on runoff should include not only the impacts from P and E 0 change, but also the impacts from the changes in other climate characteristics (such as climate seasonality and mean storm depth) and vegetation change caused by climate change.

Conclusions
As an important part of the hydrological cycle, the change of runoff can significantly affect water resources, society's safety and ecosystem health [66].Understanding the effects of climate and land cover changes on water yield is a challenging component in the assessments of future water resources [53].In this study, we used the M-C-Y equation [16,26,27], based on the widely-used Budyko framework, to quantify the temporal and spatial differences of climate elasticity of runoff in PYLB.By using the moving average method, we highlighted the influence of the catchment properties' parameter (n) variation on the climate elasticity, as well as runoff prediction.From the results obtained in this study, the following conclusions can be made: 1.
Changes in climatic variables and runoff were found using the Mann-Kendall test and Sen's slope estimator.Annual temperature in PYLB significantly increased at a rate of 1.44% (i.e., 0.19 • C) per decade.Basin-wide wind speed (U) and net radiation (R n ) had been declining at 0.17 m/s and 46.30MJ/m 2 per decade.No significant trend was detected in precipitation and relative humidity.

2.
The catchment properties' parameter is not constant during the whole study period.As we evaluated, in sub-basins, except for the Fuhe, a slight upward trend can be found during 1970-1980, followed by a decrease trend in the period from 1980 to 1990.However, the n value in Fuhe sub-basin kept increasing for the period from 1970 to 2010, which is almost persistent for the whole study period.In addition, the derived climate elasticities were significantly correlated with the catchment properties' parameter, indicating that the catchment properties' parameter was the dominant factor influencing climate elasticity in PYLB in the past 50 years.

3.
The moving window method presented in this study is relatively simple, but it is a feasible method to detect the temporal variation of climate elasticity and catchment properties.Taking the variation of the catchment properties' parameter into consideration when predicting future runoff may enhance the accuracy.

Figure 1 .
Figure 1.Topography and river networks of PYLB.Locations of the 14 meteorological stations and seven hydrological stations are presented using black triangles and red dots, respectively.

Figure 2 .
Figure 2.Changes in climatic variables and runoff in the PYLB.R, RH, U, T, R n and P represent runoff, relative humidity, wind speed, temperature, net radiation and precipitation, respectively.The shade of the cell color represents the change rate per decade.The cell labeled with two asterisks (**) denotes the statistically-significant trend at the 1% significance level.The cell without asterisks denotes that the trend is not statistically significant at both the 1% and 5% significance levels.

Figure 3 .
Figure 3. Temporal variation of the catchment properties' parameter n for PYLB and each sub-basin during the period from 1960 to 2010.The annual value of n was calibrated by adapting the Budyko framework with a 10-year moving average, including the five antecedent and precedent years.

Figure 4 .
Figure 4. Correlation between observed and estimated annual runoff within each moving window for the period 1960-2010.The y-axis value of each dot in the figure denotes the coefficient of correlation between the simulated and observed runoff in a 10-year moving window, and the shade of the dot represents the mean absolute error (MAE) between the observed and the simulated annual runoff in each 10-year moving window.

Figure 5 .
Figure 5. Temporal variation of climate elasticities across PYLB during 1960-2010.The annual value of elasticity was calculated by adapting the Budyko framework with a 10-year moving average, including the five antecedent and precedent years.P , Rn , T , U , RH and n represent precipitation elasticity, net radiation elasticity, temperature elasticity, wind speed elasticity, relative humidity elasticity and catchment properties' parameter elasticity, respectively.

Figure 6 .
Figure 6.Correlation between the calibrated n and climate elasticities in PYLB, as well as the five sub-basins.Blank denotes no significant correlation between n and climate elasticities (at the 95% significance level).

Figure 7 .
Figure 7. Absolute differences between predicted and observed runoff trend.P0 and P1 stand for prediction without consideration of n and prediction with consideration of n, respectively.

Figure 8 .
Figure 8. Changes of forest coverage and soil erosion area in PYLB from 1960 to 2010.

Table 2 .
List of hydrological gauging stations used in this study.

Table 3 .
Basic characteristics of the PYLB.P, Ē0 , R, T, Ū, RH and Rn represent the mean annual values of precipitation, evaporation, runoff depth, temperature, wind speed, relative humidity and net radiation.