A line-integral-based method to partition climate and catchment effects on runoff

It is a common task to partition the synergistic impacts of drivers in the environmental sciences. However, there is no mathematically precise solution to this partition task. Here I present a line-integral-based method, which addresses the sensitivity to the drivers throughout the drivers’ evolutionary paths so as to ensure a precise partition. The method reveals that the partition depends on both the change magnitude and pathway (timing of the change) but not on the magnitude alone unless used for a linear system. To illustrate this method, I applied the Budyko framework to partition the effects of climatic and catchment conditions on the temporal change in the runoff for 19 catchments from Australia and China. The proposed method reduces to the decomposition method when assuming a path in which climate change occurs first, followed by an abrupt change in catchment properties. The proposed method re-defines the widely used sensitivity at a point as the path-averaged sensitivity. The totaldifferential and the complementary methods simply concern the sensitivity at the initial and/or the terminal state, so they cannot give precise results. Although the path-averaged sensitivities varied greatly among the catchments, they can be readily predicted within the Budyko framework. As a mathematically accurate solution, the proposed method provides a generic tool for conducting quantitative attribution analyses.


Introduction
The impacts of certain drivers on observed changes of interest often require quantification in environmental sciences.In the hydrology community, both climate and human activities have posed global-scale impact on hydrologic cycle and water resources (Barnett et al., 2008;Xu et al., 2014;Wang and Hejazi, 2011).Diagnosing their relative contributions to runoff is of considerable relevance to the researchers and managers.Unfortunately, performing a quantitative attribution analysis of runoff changes remains a challenge (Wang and Hejazi, 2011;Berghuijs and Woods, 2016;Zhang et al., 2016); this is to a considerable degree due to a lack of a mathematically precise method of decoupling synergistic and often confounding impacts of climate change and human activities.
Numerous studies have detected the long-term variability in runoff and attempted to partition the effects of climate change and human activities through various methods (Dey and Mishra, 2017); these include the paired-catchments method and the hydrological modeling method.The pairedcatchment method can filter the effect of climatic variability and thus isolate the runoff change induced by vegetation changes (Brown et al., 2005).However, this method is capital intensive; moreover, it generally involves small catchments and experiences difficulties when extrapolating to large catchments (Zhang et al., 2011).Physically based hydrological models often have limitations such as a high data requirement, labor-intensive calibration and validation processes, and inherent uncertainty and interdependence in parameter estimations (Binley et al., 1991;Wang et al., 2013;Liang et al., 2015).Conceptual models such as Budyko-type equations (see Sect. 2.1) have consequently gained interest in recent years.
Within the Budyko framework, studies (Roderick and Farquhar, 2011;Zhang et al., 2016) have used the total differential of runoff as a proxy for the runoff change and the partial derivatives as the sensitivities (hereafter called the total-differential method).The total differential, however, is simply a first-order approximation of the observed change (Fig. 1a).This approximation has caused an error in the calculation of climate impact on runoff, with the deviation ranging from 0 to 20 × 10 −3 m (or −118 to 174 %) in China (Yang et al., 2014).The elasticity method proposed by Schaake (1990) is also based on the total-differential expression (Sankarasubramanian et al., 2001;Zheng et al., 2009).The method uses the "elasticity" concept to assess the climate sensitivity of runoff.The elasticity coefficients, however, have been estimated in an empirical way and are not physically sound (Roderick and Farquhar, 2011;Liang et al., 2015).
The so-called decomposition method developed by Wang and Hejazi (2011) has also been widely used.The method assumes that climate changes cause a shift along a Budyko curve and then human interferences cause a vertical shift from one Budyko curve to another (Fig. 2).Under this assumption, the method extrapolates the Budyko models that are calibrated using observations of the reference period, in which human impacts remain minimal, to determine the human-induced runoff changes that occur during the evaluation period.
Recently, Zhou et al. (2016) established a Budyko complementary relationship for runoff and further applied it to partitioning the climate and catchment effects.Superior to the total-differential method, the complementary method culminates by yielding a no-residual partition.Nevertheless, this method depends on a given weighted factor that is determined in an empirical but not a precise way.Furthermore, Zhou et al. (2016) argued that the partition is not unique in the Budyko framework because the path of the climate and catchment changes cannot be uniquely identified.
Obtaining a precise partition remains difficult, even when giving a precise mathematical model.This difficulty can be illustrated by using a precise hydrology model R = f (x, y), where R represents runoff, and x and y represent the climate factors and catchment characteristics, respectively.We assumed that R changes by R when x changes by x and y changes by y; i.e., R = f (x + x, y + y) − f (x, y).
To determine the effect of x on R -i.e., R x -a common practice is to assume that y remains constant when x changes by x.We thus obtain R x = f (x + x, y)−f (x, y).Similarly, we can obtain R y = f (x, y+ y)−f (x, y).Although this derivation seems quite reasonable, it is problematic as R x + R y = R.A further examination shows that a variable's effect on R seems to differ depending on the chang-ing path (timing of the change).For example, R x = f (x + x, y) − f (x, y) and R y = f (x + x, y + y) − f (x + x, y) if x changes first and y subsequently changes (note that the partition is precise with R x + R y = R at this moment).If y changes first and x subsequently changes, the partition then becomes R x = f (x + x, y + y)−f (x, y + y) and R y = f (x, y + y)−f (x, y).In the case of x and y changing simultaneously, unfortunately, the current literature seems not to provide a mathematically precise solution.
The aim of this study is to propose a mathematically precise method to conduct a quantitative attribution to drivers.The method is based on the line integer (called the LI method hereafter) and takes into account the sensitivity throughout the evolutionary path of the drivers rather than at a point as the total-differential method does.To present and evaluate the proposed method, I decomposed the relative influences of climate and catchment conditions on runoff within the Budyko framework using data from 19 catchments from Australia and China.

Methodology
2.1 Budyko framework and the MCY equation Budyko (1974) argued that mean annual evapotranspiration (E) is largely determined by the water and energy balance of a catchment.Using precipitation (P ) and potential evapotranspiration (E 0 ) as proxies for water and energy availabilities, respectively, the Budyko framework relates evapotranspiration losses to the aridity index defined as the ratio of E 0 over P .The Budyko framework has gained wide acceptance in the hydrology community (Berghuijs and Woods, 2016;Sposito, 2017).In recent decades, several equations have been developed to describe the Budyko framework.Among them, the Mezentsev-Choudhury-Yang equation (Mezentsev, 1955;Choudhury, 1999;Yang et al., 2008) (called the MCY equation hereafter) has been widely accepted and was used in this study: where n ∈ (0, ∞) is an integration constant that is dimensionless and represents catchment properties.Equation (1) requires a relatively long timescale whereby the water storage of a catchment is negligible and the water balance equation reduces to be R = P − E.Here I adopted a "tuned" n value that can obtain an exact accordance between the E calculated by Eq. (1) and that actually encountered (P − R).
The partial differentials of R with respect to P , E 0 , and n are given as  Point A denotes the initial state (the reference period), and point C denotes the terminal state (the evaluation period).The decomposition method assumes that the catchment state first evolves from Point A to B along the Budyko curve for the reference period, then jumps to Point C. R 2 represents the mean annual runoff of the evaluation period, and R 2 the mean annual runoff given the climate conditions of the evaluation period and the catchment conditions of the reference period.See Sect.2.4 for details.

∂R ∂E
Figure 3.A schematic plot illustrating the LI method.

Theory of the line-integral-based method
We start by considering an example of a two-variable function z = f (x, y) and assuming that x and y are independent.The function has continuous partial derivatives ∂z/∂x = f x (x, y) and ∂z/∂y = f y (x, y).Suppose that x and y vary along a smooth curve L (e.g., in Fig. 3) from the initial state (x 0 , y 0 ) to the terminal state (x N , y N ), and z co-varies from z 0 to z N .Let z = z N − z 0 , x = x N − x 0 , and y = y N − y 0 .
Our goal is to determine a mathematical solution that quantifies the effects of x and y on z, i.e., z x and z y .z x and z y should be subject to the constraint z x + z y = z.
As shown in Fig. 3, points M 1 (x 1 , y 1 ), . . ., www.hydrol-earth-syst-sci.net/24/2365/2020/ Hydrol.Earth Syst.Sci., 24, 2365-2378, 2020 For each segment, z i can be approximated as dz i : We thus obtain the following respective approximation of z x and z y : f y (x i , y i ) y i .Next, define τ as the maximum length of the N segments.The smaller the value of τ , the closer the value of dz i is to z i and thus the more accurate the approximations are.The approximations become exact in the limit τ → 0. Taking the limit τ → 0 then converts the sum into integrals and gives a precise expression (this is an informal derivation; please see Appendix A for a formal one): f y (x i , y i ) y i denote the line integral of f x and f y along L (termed integral path) with respect to x and y, respectively.L f x (x, y)dx and L f y (x, y)dy exist provided that f x and f y are continuous along L. We thus obtain a precise evaluation of z x and z y : Unlike the total-differential method, the sum of z x and z y persistently equals z (Appendix B).If f (x, y) is linear, then f x and f y are constant.Let f x (x, y) and f y (x, y) remain constant at C x and C y , respectively; then z x = C x x and z y = C y y. z x and z y are thus independent of L. If f (x, y) is nonlinear, however, both z x and z y vary with L, as is exemplified in Appendix C. Hence, the initial and the terminal states, together with the path connecting them, determine the resultant partition unless f (x, y) is linear.
The mathematical derivation above applies to a threevariable function as well.By doing the line integrals for the MCY equation, we obtain the desired results: where R p , R E 0 , and R n denote the effects on runoff change of P , E 0 , and n, respectively.The sum of R p and R E 0 represents the effect of climate change, and R n is often related to human activities although it probably includes the effects of other factors, such as climate seasonality (Roderick and Farquhar, 2011;Berghuijs and Woods, 2016).L denotes a three-dimensional curve along which climate and catchment changes have occurred.I approximated L by a series of line segments.R p , R E 0 , and R n were finally determined by summing up the integrals along each of the line segments (see Sect. 2.3).

2.3
Using the LI method to determine R p , R E 0 , and R n within the Budyko framework 1. Determining R p , R E 0 , and R n assuming a linear integral path: A curve can always be approximated as a series of line segments.Hence, we can first handle the case of a linear integral path.Given two consecutive periods and assuming that the catchment state has evolved from (P 1 , E 01 , n 1 ) to (P 2 , E 02 , n 2 ) along a straight line L, let P = P 2 − P 1 , E 0 = E 02 − E 01 , and n = n 2 − n 1 ; then the line L is given by parametric equations: . Given these equations, Eq. ( 2) becomes a univariate function of t, i.e., ∂R/∂P = R p (t), ∂R/∂E 0 = R E 0 (t), and ∂R/∂n = R n (t).Then, R P , R E 0 , and R n can be evaluated as 2. Dividing the evaluation period into a number of subperiods: I first determined a change point and divided the whole observation period into the reference and evaluation periods.To determine the integral path, the evaluation period was further divided into a number of subperiods.The Budyko framework assumes a steady-state condition of a catchment and therefore requires no change in soil water storage.Over a time period of 5-10 years, it is reasonable to assume that changes in soil water storage will be sufficiently small (Zhang et al., 2001).Here, I divided the evaluation period into a number of 7-year subperiods with the exception for the final subperiod, which varied from 7 to 13 years in length depending on the length of the evaluation period.
3. Determining R p , R E 0 , and R n by approximating the integral path as a series of line segments: For a short period, the integral path L can be considered as linear, which implies a temporally invariant change rate.For a long period in which the change rate varies over time, L can be fitted using a number of line segments.Given a reference period and an evaluation period comprising N subperiods, the catchment state was assumed to evolve from (P 0 , E 00 , n 0 ), . . ., (P i , E 0i , n i ), . . ., to (P N , E 0N , n N ), where the subscript "0" denotes the reference period, and "i" and "N" denote the ith and the final subperiods of the evaluation period, respectively.I used a series of line segments L 1 , L 2 , . . ., L N to approximate the integral path L, where L i connects points Then R p , R E 0 , and R n were evaluated as the sum of the integrals along each of the line segments, which were calculated using Eq. ( 6).

Total-differential, decomposition, and complementary methods
To evaluate the LI method, I compared it with the existing methods, including the decomposition method, the totaldifferential method, and the complementary method.The total-differential method approximated R as dR: where λ P = ∂R/∂P , λ E 0 = ∂R/∂E 0 , and λ n = ∂R/∂n, representing the sensitivity coefficient of R with respect to P , E 0 , and n, respectively.Within the totaldifferential method, R P = λ P P , R E 0 = λ E 0 E 0 , and R n = λ n n.I used the forward approximation, i.e., substituting the observed mean annual values of the reference period into Eq.( 2), to estimate λ P , λ E 0 , and λ n , as is standard in most studies (Roderick and Farquhar, 2011;Yang and Yang, 2011;Sun et al., 2014).
The decomposition method (Wang and Hejazi, 2011) calculated R n as follows: where R 2 , P 2 , and E 2 represent the mean annual runoff, precipitation, and evapotranspiration of the evaluation period, respectively; R 2 and E 2 represent the mean annual runoff and evapotranspiration, respectively, given the climate conditions of the evaluation period and the catchment conditions of the reference period (Fig. 2).Both E 2 and E 2 were calculated by Eq. ( 1) but using n values of the evaluation period and the reference period, respectively.The complementary method (Zhou et al., 2016) uses a linear combination of the complementary relationship for runoff to determine R p , R E 0 , and R n : where the subscripts "1" and "2" denotes the reference and the evaluation periods, respectively.a is a weighting factor and varies from 0 to 1.As suggested by Zhou et al. (2016), I set a = 0.5.Equation ( 9) thus gave an estimation of R p , R E 0 , and R n as follows: ) R n = 0.5 ∂R ∂P (P 1 + P 2 )  91 1952-1996 1997-2011 2004-2011 a R, P , and E 0 represent the mean annual runoff, precipitation, and potential evaporation, all in 10 −3 m yr −1 .n (dimensionless) is the parameter representing catchment properties in the MCY equation.AI is the dimensionless aridity index (AI = E 0 /P ).Data of catchments 1-12 were derived from Zhang et al. (2010).Data of catchments 13-16 were from Sun et al. (2014).Data of catchments 17, 18, and 19 were from Zheng et al. (2009), Jiang et al. (2015), and Gao et al. (2016), respectively.I used the change points given in the literatures to divide the observation period into the reference and elevation periods.The LI method further divides the evaluation period into a number of subperiods.
The column "The final subperiod" denotes the final subperiod, which was used as the evaluation period for the total-differential method, the decomposition method, and the complementary method.

Data
I collected runoff and climate data from 19 selected catchments evaluated in previous studies (Table 1).The changepoint years given in these studies were directly used to determine the reference and evaluation periods for the LI method.
As mentioned above, the LI method further divides the evaluation period into a number of subperiods.For the sake of comparison, the final subperiod of the evaluation period was used as the evaluation period for the decomposition, the total-differential, and the complementary methods (it can be equally considered that all of the four methods used the final subperiod as the evaluation period, but the LI method cares about the intermediate period between the reference and the evaluation periods, and the other methods do not).Eight of the 19 catchments had a reference period comprising only one subperiod (Table 1), and the others had two to seven subperiods.
The 19 selected catchments have diverse climates and landscapes, with 12 from Australia and seven from China (Table 1).The catchments span from tropical to subtropical and temperate areas and from humid to semi-humid and semiarid regions, with the mean annual rainfall varying from 506 to 1014 × 10 −3 m and potential evaporation from 768 to 1169 × 10 −3 m.The dryness index ranges between 0.86 and 1.91.The catchment areas vary by 5 orders of magnitude from 1.95 to 121 972 with a median 606 × 10 6 m 2 .The key data include annual runoff, precipitation, and potential evaporation.The record length varied between 19 and 76 with a median of 39 years.All the catchments experienced changes in climate and catchment properties over the observation periods.The precipitation changes from the reference to the evaluation period ranged between −153 and 79 × 10 −3 m yr −1 , and between −35 and 41 × 10 −3 m yr −1 for potential evaporation (Table 2).The coeval change in the parameter n of the MCY equation ranged between −0.2 and 1.4.The mean annual streamflow reduced for all catchments, ranging from 0.4 to 169 with a median 38×10 −3 m yr −1 .The change in catchment properties mainly refers to the vegetation cover or land use change.More details on the data and the catchments can be found in Zhang et al. (2010Zhang et al. ( , 2011)), Sun et al. (2014), Zheng et al. (2009), Jiang et al. (2015), and Gao et al. (2016).The subscript "1" denotes the reference period, and "2" denotes the evaluation period.X=X2 − X1 (X as a substitute for R, P , E 0 , and n).

Catchment LI method Decomposition
Total-differential Complementary no.method method method  n 1 f n (P 2 , E 02 , n)dn, where the subscript "1" denotes the reference period and "2" denotes the final subperiod of the evaluation period.

Results
Table 3 lists the resultant values of R p , R E 0 , and R n from the LI method and the three other methods.Please see the Supplement for detailed calculation steps.Figure 4a compares the resultant R n of the LI method and the decomposition method.Although they are quite similar, the discrepancies can be >20×10 −3 m yr −1 .The decomposition method assumes that climate change occurs first and then human interferences cause a sudden change in catchment properties (Fig. 2).Such a fictitious path is identical to the path ABC in Fig. 3, provided that x represents climate factors and y catchment properties.When ABC was adopted as the integral path, the LI method yielded the same results as the decomposition method did (Fig. 4b).Hence, the decom-position method can be considered as a case of the LI method that uses a special integral path.
The total-differential method is predicated on an approximate equation, i.e., Eq. ( 7).The LI method reveals that the precise form of the equation is R = λ P P + λ E 0 E 0 + λ n n (Appendix D), where λ P , λ E 0 , and λ n denote the pathaveraged sensitivity of R to P , E 0 , and n, respectively.All points along the path have the same weight in determining λ P , λ E 0 , and λ n .To determine them, the total-differential and the complementary methods utilize only the initial and/or the terminal states.Neglecting the intermediate states results in an imprecise partition, as was illustrated in Figs. 5 and 6 as well as in Fig. 1 using a univariate function, and even a reverse trend estimation (see R E 0 for catchment no. 1 in Table 3).As with the LI method, the complementary method produced R P , R E 0 , and R n that exactly summed up to R.Although its resultant R p , R E 0 , and R n values were all in accordance with the LI method (Fig. 6), the LI method often yielded values beyond the bounds given by the complementary method (Fig. 7); this is because the maximum or minimum sensitivities do not necessarily occur at the initial or terminal states.
λ P , λ E 0 , and λ n imply the average runoff change induced by a unit change in P , E 0 , and n, respectively (Appendix D). λ P , λ E 0 , and λ n all varied by up to several times or even 10-fold between the studied catchments (Table S4). Figure 8 shows that Eq. (2) reproduced λ P , λ E 0 , and λ n very well when taking the long-term means of P , E 0 , and n as inputs, a reflection of the fact that the dependent variable approached its average when the independent variables were set to be their averages.This finding is of relevance to the spatial prediction of λ P , λ E 0 , and λ n .

Discussion
The LI method highlights the role of the evolutionary path in determining the resultant partition.Yet it seems that no studies have accounted for the path issue while evaluating the relative influences of drivers.The limit of the LI method is high data requirements for obtaining the evolutionary path.When the path data are unavailable, the complementary method can be considered as an alternative.The complementary method is free of residuals; moreover, it employs data from both the reference and the evaluation periods, thereby generally yielding sensitivities closer to the path-averaged results than the total-differential method.
While using the Budyko models, a reasonable timescale is relevant to meet the assumption that changes in catchment water storage are small relative to the magnitude of fluxes of P , R, and E (Donohue et al., 2007;Roderick and Farquhar, 2011).A 7-year timescale was used in the present study, as most studies have suggested that a time period of 5-10 years (Zhang et al., 2001(Zhang et al., , 2016;;Wu et al., 2017a, b;Li et al., 2017) or even 1 year (Roderick and Farquhar, 2011;Sivapalan et al., 2011;Carmona et al., 2014;Ning et al., Figure 8. Performances of Eq. ( 2) to be used to predict λ P , λ E 0 , and λ n with the long-term mean values of P , E 0 , and n as inputs.
is the mean absolute error, where O and P are values that were actually encountered (given in Table S4) and predicted using Eq. ( 2), respectively, and N is the number of selected catchments.2017) is reasonable.Nevertheless, some studies have argued that the time period should be longer than 10 years (Li et al., 2016;Dey and Mishra, 2017).Using the Gravity Recovery and Climate Experiment (GRACE) satellite gravimetry, Zhao et al. (2011) detected the water storage variations for the three largest river basins of China, namely, the Yellow, Yangtze, and Zhujiang.The Yellow River mostly drains an arid and semiarid region (P , 450 × 10 −3 m; R, 70 × 10 −3 m; E, 380 × 10 −3 m), and the Yangtze (P , 110 × 10 −3 m; R, 550 × 10 −3 m; E, 550 × 10 −3 m) and the Zhujiang River basins (P , 1400×10 −3 m; R, 780×10 −3 m; E, 620×10 −3 m) are humid.The amplitude of the water storage variations between years was 7, 37.2, and 65×10 −3 m for the three rivers, respectively, at 1 order of magnitude smaller than the fluxes of P , R and E.Although the observations cannot be directly extrapolated to other regions, the possibility seems remote that the use of a 7-year aggregated time period strongly violates the assumption of the steady-state condition.
The mutual independence of the drivers is crucial for a valid partition.In the present study, although annual P and E 0 exhibited significant correlation for most catchments (p<0.05), the aggregated P , E 0 , and n over a 7-year period showed minimal correlation (mostly p>0.1).The interdependence between the drivers can considerably confound the resultant partitions of the LI method and other existing methods.
The LI method revises the concept of sensitivity at a point as the path-averaged sensitivity.Mathematically, the LI method is unrelated to a functional form and hence applies to communities other than just hydrology.For example, identifying the carbon emission budgets (an allowable amount of anthropogenic CO 2 emission consistent with a limiting warming target) is crucial for global efforts to mitigate climate change.The LI method suggested that the emission budgets depends on both the emission magnitude and pathway (timing of emissions), which is in line with a recent study by Gasser et al. (2018).An optimal pathway would facilitate an elevated carbon budget unless the carbon-climate system behaves in a linear fashion.
This study presented the LI method using time-series data, but it applies equally to the case of spatial series of data.Given a model that relates fluvial or eolian sediment load to the influencing factors (e.g., rainfall and topography), for example, the LI method can be used to separate their contributions to the sediment-load change along a river or in the along-wind direction.

Conclusions
Based on the line integral, I created a mathematically precise method to partition the synergistic effects of several factors that cumulatively drive a system to change from one state to the other.The method is relevant for quantitative assessments of the relative roles of the factors behind the change in the system state.I applied the LI method to partition the effects of climatic and catchment conditions on runoff within the Budyko framework.The method reveals that, in addition to the change magnitude, the change pathways of climatic and catchment conditions also play a role.Instead of using the runoff sensitivity at a point, the LI method uses the pathaveraged sensitivity, thereby ensuring a mathematically precise partition.As a mathematically accurate scheme, the LI method has the potential to be a generic attribution approach in the environmental sciences.

Figure 1 .
Figure1.For a nonlinear function z = f (x), the total-differential method (a) and the complementary method (b) do not accurately estimate the effect ( z x ) of x on z when x changes by x, but the LI method (c) does.For a univariate function, the z change is exclusively driven by x, so that z x should be equal to z. z x = z in (c) but not in (a) and (b).λ x in (c) represents the average sensitivity along the curve AC and λ x = z/ x; see Appendix E for details.

Figure 2 .
Figure2.A schematic plot illustrating the decomposition method.Point A denotes the initial state (the reference period), and point C denotes the terminal state (the evaluation period).The decomposition method assumes that the catchment state first evolves from Point A to B along the Budyko curve for the reference period, then jumps to Point C. R 2 represents the mean annual runoff of the evaluation period, and R 2 the mean annual runoff given the climate conditions of the evaluation period and the catchment conditions of the reference period.See Sect.2.4 for details.

Figure 4 .
Figure 4. Comparisons between the LI method and the decomposition method.(a) Comparison of the estimated contributions to the runoff changes from the catchment changes ( R n ).(b) The decomposition method is equivalent to the LI method that assumes a sudden change in catchment properties following climate change.In this case, the integral path of the LI method can be considered as the path ABC in Fig. 3 (x represents climate factors, and y catchment properties, i.e., n) and R n = AB+BC ∂R ∂n dn = AB ∂R ∂n dn + BC ∂R ∂n dn = 0 + BC ∂R ∂n dn =

Figure 5 .
Figure 5. Comparisons of the estimated contribution to runoff from the changes in (a) precipitation ( R P ), (b) potential evapotranspiration ( R E 0 ), and (c) catchment properties ( R n ) between the LI method and the total-differential method.

Figure 6 .
Figure 6.Comparisons of (a) R P , (b) R E 0 , and (c) R n between the LI method and the complementary method (a = 0.5).

Figure 7 .
Figure 7. Comparisons of (a) R P , (b) R E 0 , and (c) R n by the LI method with the upper and lower bounds given by the complementary method.According to Zhou et al. (2016), R P , R E 0 , and R n reach their bounds when a is 0 or 1.

Table 1 .
Summary of the long-term hydrometeorological characteristics of the selected catchments a .