Next Article in Journal
Importance of Flood Samples for Estimating Sediment and Nutrient Loads in Mediterranean Rivers
Next Article in Special Issue
Effect of Rainfall Regime on Rainwater Harvesting Tank Sizing for Greenhouse Irrigation Use
Previous Article in Journal
What Is the Contribution of Urban Trees to Mitigate Pluvial Flooding?
Previous Article in Special Issue
Incorporating aSPI and eRDI in Drought Indices Calculator (DrinC) Software for Agricultural Drought Characterisation and Monitoring
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Procedure for Estimating Drought Duration and Magnitude at the Uniform Cutoff Level of Streamflow: A Case of the Weekly Flows of Canadian Rivers

Department of Civil Engineering, Lakehead University, Thunder Bay, ON P7B 5E1, Canada
*
Author to whom correspondence should be addressed.
Submission received: 28 May 2022 / Revised: 10 June 2022 / Accepted: 14 June 2022 / Published: 16 June 2022
(This article belongs to the Special Issue Drought and Water Scarcity: Monitoring, Modelling and Mitigation)

Abstract

:
At times, hydrological drought is defined using Q90 or Q95 (90% or 95% flows equaling or exceeding) or even at higher levels, such as Q75 as the cutoff level regardless of their seasonal variation (i.e., truncation at the uniform flow level). In the past, the estimation of drought length and magnitude at the aforesaid uniform cutoff levels of flow has been a challenging issue. A procedure is presented to first estimate the drought magnitude (M), which then forms the basis for estimating the drought duration or length (L). The drought magnitude (M) and the length of the critical period (Lcr) are estimated using the concept of behavior analysis prevalent in the hydrologic design of reservoirs. This information is used for estimating the drought length (LT-e′, the estimated value of drought length for the return period of T weeks) involving a Markov chain model on the standardized weekly flow sequences. A weighted average of Lcr and LT-e′ (=0.60 Lcr + 0.40 LT-e′) results in the estimate of drought length, which is compatible to the observed counterpart. The performance of the procedure to estimate drought length was found to be satisfactory up to the truncation level of Q75, whereas the estimation of drought magnitude was rated as good.

1. Introduction

The annual, monthly, and daily identification and modeling of hydrological droughts has been the subject of investigations for which significant research material is available in the literature. At times the hydrological droughts have been defined when the flow drops below a uniform predefined level such as Q95 (95% flows are equal to or exceeding this level of flow or the probability of drought, q = 5% or 0.05), Q90, Q80, or Q70 [1,2,3,4], or even the mean or median (=Q50) and 75% of the mean [5], generally in case of daily flow sequences. Lately, an attempt has also been made to quantify drought duration and deficit volume on a weekly basis by Sharma and Panu [6,7]. A weekly time scale is a good compromising resolution between a long period of a month and a short period of a day. Put another way, a weekly analysis of drought is simple and tractable in a statistical sense, whereas a daily analysis might catapult into a complex relationship, which at times tends to be intractable. Furthermore, a weekly analysis is vindicated by the outcome of low flow studies in which a minimum annual 7-day average flow is found to be a rational descriptor of drought conditions by the task committee for low flow conditions ((TCLFE) [8] in USA and Gustard et al. [9] in UK). One major advantage of a weekly analysis is that it can easily be scaled down to daily values in terms of drought duration (or length). For example, if a 50-year drought extends for 10 weeks based on the cutoff level of Q90, then a drought period spanning 70 (=10 × 7) days can be inferred.
In a majority of hydrological drought studies dwelling on a seasonal basis, the truncation level is variable in the sense that statistics of the respective month, season, or week are used to truncate the series. For instance, if we are dealing with monthly flows, then droughts are identified with respect to monthly means or medians [6,10,11] or some statistic of respective monthly flow such as Q70 of the respective month [12]. In a statistical sense, the variable flow level simplifies the analysis by which the flow series can be easily made stationary and the resultant series can be analyzed using the known concepts of stochastic theory. Such a stationarization can be achieved through week-by-week standardization of the weekly flows, with the resultant series termed as the standardized hydrological index, SHI series [6]. The SHI series can be regarded as analogous to a standardized precipitation index (SPI) series [13,14] that is used in the ambit of meteorological drought. In this kind of analysis, although the truncation level is uniform in terms of probability, the uniformity of the truncation in terms of the flow is not preserved. The term uniform probability implies that when the analysis is done in the standardized domain of flows (SHI sequences), then the truncation is done at the uniform value of SHIx, with the uniform counterpart drought probability, qx. In doing so, constancy of qx is preserved, but the constancy of the cutoff flow is distorted.
A recent study [15] has used meteorological and hydrological (streamflow) anomaly indices to assess water needs during the hydrological drought on a monthly basis. The above study is based on the truncation level approach and the theory of runs. There has been an introduction of copula and entropy based approaches to characterize the drought duration and the deficit volumes [16,17], and these approaches are gaining prominence, though their efficacy is still to be evaluated in comparison to the traditional approaches advanced by Dracup et al.[10] and Yevjevich et al. [11].
In practice, the information on drought parameters (i.e., duration and magnitude) at the uniform flow level of truncation is of greater importance, which is being addressed in this paper. The parameters of the hydrologic drought, viz., T-week drought duration, LT and magnitude, MT (standardized deficit volume), are of practical utility at uniform truncation levels of flow for the design and the operation of drought alleviating facilities and structures. The standardized deficit volume is named here as drought magnitude MT (=DT/σ) in which DT is the deficit volume in volumetric units and σ is a characteristic standard deviation of the flow sequence, thus MT is a dimensionless entity and suffix T stands for the return period (week). In the traditional approach, the drought length is estimated first (say by the Markov Chain methodology or the extreme number theorem), which is transformed into drought magnitude through the linkage relationship (drought magnitude = drought intensity × drought length [10]). In the present analysis, firstly the drought magnitude is determined, which then forms the basis for estimating the drought length. The concepts of behavior analysis [18,19] used in the hydrologic design of reservoirs, simple and conditional probabilities, and the Markov chain using the SHI sequences are invoked to establish T-week drought parameters (duration and magnitude). The weekly flows have been subjected to drought analysis with a uniform truncation level of Q95 and Q90 up to Q75.

2. Background of the Model

To begin with, the weekly flow data for the Upper Humber River (#23, Table 1) were analyzed for the reservoir volume (VR) as well as DT (drought deficiency volume) at the truncation (demand level) of Q75. It should be noted that the Upper Humber River was chosen because of its size (2210 km2) and its long uninterrupted length of flow record of 68 years during which the persistence in the weekly SHI sequences has been modestly high (ρ = 0.48), with the coefficient of variation (=1.07). Most importantly, the watershed is least provoked by abstractions, and flows can be deemed entirely in a natural state. For this river, the Q75 flow based on the weekly flow duration curve of the historical data spanning from 1953 to 2020 (68 years = 3536 weeks) was found to be equal to 24.16 m3/s. In the behavior analysis [18], the weekly flow sequence was truncated at the uniform flow level =24.16 m3/s (demand level) and several episodes of full reservoir conditions (assumed to be above the cutoff level) were tracked. Behavior analysis was chosen as it has been found to be a simple procedure to estimate reservoir capacity, and it can be applied to the data of any time interval. The procedure takes into account autocorrelation, seasonality, and other flow parameters imbued in the historical data [18]. Further calculations were done using the water balance equation (i.e., behavior analysis), Vt = Vt−1 + Qt − Q75 (all losses including evaporation were neglected for the sake of simplicity in the drought analysis). In the above water balance equation, Vt is the reservoir volume at the end of week t, Vt−1 is the volume at the end of the week t−1, Qt is the river inflow during week t, and Q75 is the outflow needed to meet the demand during the week t. In the process of counting, Vt < Q75 indicates no fill condition (deficit conditions) and Vt ≥ Q75 indicates the reservoir full condition (surpluses). Under the condition, Vt > 0, the reservoir will spill over the water from the reservoir say through a spillway. The calculations are then begun with Vt−1 = 0, and all deficits are counted below zero as magnitudes with a negative sign. All the values of Vt ≥ 0 are set to zero. The deficit spell begins with zero, drops down to the minimum level and then begins to recover and ends again with zero after full recovery. All the values during this dry spell are negative and the minimum negative value represents the maximum deficit volume (VR) during this dry period, i.e., critical period, Lcr. By applying the above water balance equation and a counting scheme, the maximum value of VR was found = 221.99 m3/s-week (absolute value) with the Lcr equal to 21 weeks.
In the drought length and magnitude based analysis [7], the flow series of 3536 weeks is chopped at a uniform flow level of 24.16 m3/s and several episodes of deficit (marked by 0) and surplus (marked by 1) conditions erupted. The episode with the largest length LT and DT (T = 3536) was identified, and these entities were estimated by a counting method. The LT (largest drought length) and the DT (largest deficit volume) were found = 17 weeks and 221.99 m3/s-week, respectively. It can be seen that both VR and DT are equal but LT (=17) is not equal to Lcr (=21). This behavior of equality of DT and VR was found in almost all the rivers at the chosen truncation levels. The above calculations were done using the historical observed data. So, the terms LT and DT are designated as LT-o and DT-o (subscript “o” stands for observed) from here onward. When LT and DT are estimated by a modeling scheme, they are designated by LT-e and DT-e (subscript “e” stands for estimated).
Following the aforesaid calculations, the behavior of VR, DT-o, Lcr, and LT-o was investigated for all rivers [Table 1]. It was found that on a weekly basis, the value of DT-o using flow sequences was almost equal to VR at the truncation levels from Q95 to Q75 for the Canadian rivers having ρ less than 0.90 [Table 2, Table 3 and Table 4]. That means the estimates of VR can be obtained using behavior analysis and the estimate of VR thus obtained can be standardized and designated as VR′. The value of VR′ can be equated to the standardized DT-o (termed as drought magnitude MT-o). Though the estimation of drought magnitude DT can be easily accomplished by estimating VR by behavior analysis, the estimation of LT requires further investigation. The estimate of Lcr was generally found > LT-o, so a need arose to determine other estimates of LT-e that could be combined with Lcr, with the resulting value equivalent to LT-o. A method to estimate the LT-e to be equivalent to LT-o is the subject matter of the present paper.

Estimation of Drought Length, LT-e

It is mentioned in the foregoing section that Lcr is one estimate of LT-e, which needs further work because of the discrepancy between LT-o and Lcr. To this end, recourse was taken for the estimation of LT-e by a Markov chain (MC) based relationship. To make use of the MC relationship, the weekly flow sequences must be stationary (i.e., stationarity up to the second order). The historical weekly flow sequences can be stationarized using the week-by-week standardization of flow sequences thus forming the SHI sequences. These SHI sequences form the platform for applying the MC based modeling methodology to derive an estimate of LT-e [6].
Since MT-e is a standardized value of the deficit volume, DT-e, therefore DT-e = σav × MT-e, where σav is the average of 52 values of the weekly standard deviations. It should be noted that σav was found to be the best estimator of the standard deviation for the above transformation (Sharma and Panu [20,21]). Further, the value of MT-o from the non-standardized weekly sequences is computed as MT-o = [VRav] = VR′ = [DT-oav], subject to the condition VR = DT-o, which is true for the range of uniform cutoff, Qx such as Q90, considered in the present analysis.
Since VR′ = MT-o, to find LT-e, VR′ is set = MT-e. In the process of equating VR′ to MT-e, the following steps are pursued:
Step-1: Choose a suitable cutoff level corresponding to flow level, Qx (say Q90). The standardized version of the cutoff level (denoted as SHIx) can be obtained in the following fashion.
Step-2: Let (i) = SHIo = (Qx − µo)/σo; (ii) = SHImax = (Qx − µo)/σmax; (iii) = SHIav = (Qx − µo)/σav; or (iv) = SHIgm = (Qx − µo)/σgm, where σmax is the maximum value, σav is the arithmetic mean, and σgm is the geometric mean of the 52 values of the weekly standard deviations. It has been mentioned earlier in the text that µo and σo are the overall mean and the standard deviations of the weekly historical (non-standardized) flow sequences. At times, cutoff levels [average of (i) and (ii) denoted as (ia)] or [average of (iii) and (i), denoted as (iiia)] can be combined to improve the results of the analysis.
To elucidate the point, there shall be six cutoff levels: (iv), (iii), (iiia), (i), (ia), and (ii). These levels, in general, in the ascending order are ((ii) being the highest), the most frequent levels being (iiia) and (i), whereas levels (ii) and (iv) are rare. Based on the chosen cutoff level SHIx, the probabilities q1, qq, and qp are computed by a counting method. At the Q95 level (drought probability q = 0.05), the q1 is computed by the chosen cutoff level and it should be close to 0.05. Likewise, at the Q75 level, q1 should be close to 0.25. So, the identification of an appropriate cutoff level is a trial and error procedure. It should be noted that q1 is the drought probability computed by a counting method and q is the drought probability based on the flow duration curve, and they can be designated as qx corresponding to cutoff flow Qx.
Step-3: Once an appropriate cutoff level is chosen, obtain an estimate of LT-e′ by using the MC1 [7] representation of drought lengths, expressed as follows.
L T - e = 1 log   F   T 1 q 1 q p log q q
where, F is the factor to account for the plotting position in the empirical estimation of the exceedance probability. In this analysis, the plotting position formula [22] developed for Canadian rivers has been used. The formula evaluates the exceedance probability (=0.75/(T + 0.25)), so F (=1.33 (1 + 0.25/T)) ≈ 1.33 as T is generally large. The term qq stands for the conditional probability of the present period being drought given the previous period was also a drought and likewise qp stands for the present period being drought given the previous period was wet. The probabilities qq and qp can be computed from the counting method. It should also be borne in mind that at low cutoff levels, such as Q95, LT-e′ may be adequately represented by MC0 in which case qq = qp = q1.
Step-4: Obtain the initial estimate of MT-e by following the linear linkage relationship between drought magnitude (M) = drought intensity (Id) × drought length (L) [10]. Thus, the MT-e can be expressed as
MT-e = µd × Lc
in which Lc is a characteristic drought length obeying the following relationship (Sharma and Panu [20,21]).
L c =   Φ   L m + 1 Φ     L T - e
The value of the parameter Φ varies between 0 and 1, and its proper value can be estimated by a trial and error procedure, which will be explained in the following text.
The mean drought length ( L m )   can be expressed as Sen [23]
L m = 1 1 q q
It is tacitly assumed that µd, mean of the drought intensity, is approximated as independent of T, and it can be expressed by the truncated normal pdf of drought intensities as follows [23,24].
µ d = exp 0.5   z 0 2     q n 2 π z 0
In Equation (5), qn is the drought probability corresponding to z0 (standard normal variate). The z0 can be regarded as equivalent to the cutoff level (or truncation level), SHIx. The qn can be determined from the standard normal tables or the polynomial function [25]. In strict normal pdf conditions of SHI sequences, the value q1 by the counting procedure is equal to qn. Otherwise, minor deviation may occur. However, in the calculations, qn based on the polynomial function should be used to ensure consistency.
Equation (2), involving µd, can be reduced to the following form.
M T e =   abs exp 0.5   z 0 2     q n 2 π z 0 L c
If the parameter Φ = 0, then Equation (6) takes the following form.
M T - e = abs exp 0.5   z 0 2     q n 2 π z 0 L T - e
After having chosen the value of the cutoff level (z0 = SHIx) and estimating the parameters (q1, qq, and qp) by the counting method, the LT-e′ and the MT-e′ are computed. MT-e′ (means Φ = 0) value so computed may be compared to VR′. If MT-e′ ≥ VR′, then the model structure (MC1) is acceptable, and MT-e can be made equal to VR′ by adjusting the value of Φ within its range 0 and 1. If MT-e′ < VR′, then matching should be done between MT-e′ (based on LT-e′ with MC0 structure of the model) and VR′. Again, MT-e can be made = VR′ by adjusting the value of Φ within the range 0 and 1. This situation was found to occur with flow levels of Q95 and at Q90 (on a few occasions), as can be seen in Table 3 in the context of the Beaverbank River. The value of LT-e′ so obtained should be combined with Lcr to arrive at an estimate of the LT-e as was discussed earlier in the text. The detailed calculations based on the modeling methodology are presented in the next section along with a flow diagram.

3. Data and Methods of Analysis

3.1. Data Acquisition

Data for the analysis constituted natural (i.e., unregulated) and uninterrupted flow records of 24 rivers across Canada (Figure 1), which are listed in Table 1. Twenty rivers were used for fitting the model whereas four rivers (written in italics) were used for validation of the model. Daily flow data for these 24 rivers were extracted from the Canadian Hydrologic Data Base, HYDAT, Environment Canada [26]. Selected rivers are representative of a wide range of drainage basins (37 to 32,400 km2) and a long period of the historical database (1911 to 2020). Daily flows were transformed into weekly flows such that each of the first 51 weeks would be composed of 7 days, while week 52 would contain the remaining days. That is, the last week of the year would comprise 8 or 9 days (in case of a leap year).

3.2. Computation of Flow Statistics and Probabilities

The statistics such as μo (the overall mean), cvo, cvmx, cvav, cvgm, and ρ (lag-1 autocorrelation from the SHI sequences), were computed from historical weekly flow data, as shown in Table 1. Four variants of standard deviations were computed: σo, the overall standard deviation without standardizing the weekly flows; σmx, the maximum value among the 52 standard deviations; σav the arithmetic mean of 52 standard deviations; and σgm the geometric mean of these 52 standard deviations. The corresponding coefficients of variation were computed as cvo (=σoo), cvmx (=σmxo), and so on.
The observed values of drought duration (LT-o) and drought deficiency volume (DT-o) were computed using the historical weekly flow data. The computations are illustrated using three rivers, viz., the Athabasca (River #2), Goulis (River #10), and the Upper Humber (River #23) in Table 1. Likewise, behavior analysis as described briefly in Section 2 was used to estimate the VR at the cutoff level of Q95, Q90, --- Q75. The length of the critical period, Lcr was also computed from behavior analysis. In the majority of instances, DT-o was found equal to VR. The entities DT-o and VR were standardized using σav as a common denominator such that MT-o (=DT-oav) and VR′ (=VRav). The standard deviation σav was found to be the best estimator of MT in the modeling mode as stated before [20,21], which is why it was used as a common denominator in the standardization process. For almost all the conditions of cutoff levels, VR′ was found equal to MT-o (Table 2, column 2). In other words, VR′ can be taken as an estimate of drought magnitude, which should compare well with the estimated value (denoted by MT-e) by the modeling methodology.
The parameters q1, qp, and qq were estimated by a counting method [6,7] by truncating the SHI series at the desired truncation (cutoff) level of SHIx. Values above the truncation level were designated as surplus ‘1’ and below the truncation level as deficit ‘0’. Consider a number of scenarios involving occurrences of ‘0’ as follows. The occurrences of 0 are, isolated single zeros (say n1), and consecutive two 0’s, i.e., ‘00’ (say n2). Knowing that the total sample size is n, hence q1 (=n1/n), qq (=n2/n1). Since there are only two numerals, namely 1 and 0, the number of 1’s = (n − n1). The number of pairs ‘11’ occurring in succession (i.e., 1 preceded by 1) is counted (say n2). Therefore, pp (=n2/(n − n1)) and hence qp (= (1 − pp)). The probabilities q1, qp, and qq were computed at each of the chosen cutoff levels (i.e., Q95, Q90, Q85, Q80, and Q75).

3.3. Identification of Cutoff Level

Identification of the proper cutoff level (SHIx) plays a crucial role. A value of SHIx could take one of six forms, as discussed in Section 2. In the process of choosing the right value of SHIx, one can begin from the level (iii) SHIx = SHIav (= (Qx − µo)/σav)) and compute the values of q1, qq, qp, LT-e′, and MT-e′. At the chosen cutoff level, VR′ should be compared to MT-e′ and the value of MT-e′ should be ≥ VR′. Otherwise, the next cutoff level should be considered until the criterion MT-e′ ≥ VR′ has been satisfied. It should be borne in mind that MT-e′ is a value of MT-e that is based on Φ = 0. The value of Φ (=0) results in the highest value of MT-e, whereas Φ (=1) yields the lowest value. Thus, MT-e′ can be computed with MC0 or MC1 based LT-e′, and matching can be done between VR′ and MT-e′ to confirm the adequacy of the order of MC applicable at a particular Qx. This is particularly significant at Q95 and Q90 flow levels. All the above steps of calculations are depicted in the flow chart shown in Figure 2.
For elucidation, an example of the Upper Humber River (Table 2) is presented, at the Q75 level with cutoff level (iii) and MC1 based drought lengths, VR′ (=4.18) (behavior analysis), and MT-e′ (=2.84). Thus, the criterion (MT-e′ ≥ VR′) was not satisfied, so prima-facie, cutoff level (iii) seemed to be inadequate. To meet the criterion, MT-e′ should have a larger value, therefore, the next higher cutoff level (iiia) was tried. Using this cutoff level, the corresponding probability values (column 4), MT-e′ and LT-e′ (column 5) were computed as shown in Table 2. For these computations, the criterion MT-e′ ≥ VR′ was satisfied, so cutoff level (iiia) seemed adequate. However, looking at the values of probabilities in column 4, it can be seen that q1 (=0.117) is much lower than the expected value of q = 0.25 at Q75. So, there is further scope for adjustment in the cutoff level, therefore, the next level (i) = [(Qx − µo)/σo] was attempted. Cutoff level (i) turned out to be appropriate in the sense that criterion MT-e′ ≥ VR′ was satisfied as well as q1 (=0.242) resulted in being much closer to q = 0.25. Using cutoff level (i), the estimates of q1, qq, and qp, (column 4) and LT-e′ (column 5) were obtained. Likewise, by trial and error, Φ was found= 0.47 (column 7), which rendered MT-e (=4.20) closer to VR′ (=4.18). Thus, the calculations resulted in LT-e′ = 11 weeks and Lcr = 21 weeks. These two values were combined to obtain (a) LT-e = (0.5Lcr + 0.50 LT-e′) = 16; (b) = (0.75 Lcr + 0.25 LT-e′) = 19; and (c) = (0.60 Lcr + 0.40 LT-e′) = 17, which are shown in column 8. The value based on the weighted averaging option (c) compared with the observed value (LT-o = 17 in column 2) perfectly. Likewise, at Q80, Q85, and Q90 flow levels, similar calculations were repeated and cutoff level (i) was found adequate with MC1 for drought lengths. The relevant calculations are shown in Table 2.
Because, at the Q95 demand flow level, the values of VR′ and Lcr are, respectively, 0.70 and 9, there arose a need to revisit cutoff level (i) with MC1 for the drought lengths. This combination resulted in MT-e′ (=4.99) >> VR′ (=0.70), meaning that MC1 for length is an over-fit. Whereas for the MC0 condition MT-e′ = 2.20, which compares relatively closer to VR′ (=0.70). Although, with MC0 for the drought length, the criterion of drought magnitude was satisfied, yet the q1 value was equal to 0.138 (column 4), which is much higher than the expected value of 0.05 (Q95). This called for a need to lower the cutoff level to (iiia). By doing so, q1 reduced to 0.043 (column 4) with MT-e′ (=1.27), which is > VR′ (=0.70). The calculations thus revealed that at the Q95 flow level, cutoff level (iiia) with MC0 for drought length is more appropriate. Using the revised value of q1 and LT-e′, the appropriate value of Φ was sought, which turned out to be 0.97 (column 7) and MT-e (=0.70) (column 7). Accordingly, the final three values of LT-e with three options of averaging are shown in column 8. At the Q95 cutoff level, option (a), i.e., the arithmetic averaging of LT-e߰ and Lcr matched LT-o perfectly (column 8).
Following the calculations for the Upper Humber River, all rivers (Table 1) were analyzed. The results based on the analysis of three typical rivers are summarized in Table 3. The three typical rivers from Table 1 are: the Athabasca River (#2 from Canadian Prairies), the Goulis River (#11 from northern Ontario), and the Beaverbank River (#21 from Atlantic Canada. These rivers represent a typical scenario in terms of the MC structure of drought lengths and cutoff levels. In a majority of cases, it was found that the optimal level of cutoff level yielded q1 values close to q (=qx) at Q75 or Q80 flow levels. Therefore, it is advisable to commence the analysis from Q75 downward to the Q95 level. At Q95 and Q90 flow levels, the values of q1 are very sensitive to SHIx and so q1 may not correspond closely to q or qx. Under such an ambiguous scenario, a SHIx level that is one step lower than at Q75 and Q80 can be adopted (Table 3).

4. Results and Discussion

4.1. Estimation of MT-e and LT-e: Fitting the Model Structure

The model structure as outlined in the aforesaid section was fitted on 20 rivers identified by bold numerals (Table 1). The remaining four rivers with italicized numerals and identity (#3 Bow River in Canadian prairies, #4 Pipestone River, western Ontario, #14 Chippewa River in northern Ontario, and #24 Lepreau River in Atlantic Canada) were used for model validation. The adequacy and the quality of fitting was adjudged by the Nash–Sutcliffe (NSE) efficiency and associated mean error (MER) criteria [20].
The most important element discovered from the analysis was that in almost all cases VR′ = MT-o (column 2, Table 3 and Table 4). This means one can safely use the equality VR′ = MT-e as a pivotal point to calculate the values of Lcr and LT-e′ for a given condition of the uniform cutoff level, Qx. Thus, at a given Qx, VR′ and Lcr remain unchanged with respect to any value of corresponding SHIx [i.e., (i) for (Qx − µo)/σo, (iii) for (Qx − µo)/σav, etc.]. However, MT-e′ does change with respect to the aforesaid SHIx levels because of a change in the values of probabilities and thereby in LT-e′ and µd (Equation (7)). Therefore, compute VR′ (fixed value) from the behavior analysis and compare it with MT-e′ computed from Equation (7). So, while moving the cutoff level SHIx up and down, the corresponding probabilities, q1, qq, and qp, are computed by the counting method and MT-e (with Φ, a value of which is found by the trial and error) is calculated. There are two estimates of MT-e, namely, (1) MC0 based on LT-e′ and µd and (2) MC1 based LT-e′ and µd. The combination that yields the closest correspondence between VR′ and MT-e is chosen.
For illustrative purposes, in the three selected rivers, the appropriate combination for various Qx levels is shown in Table 3. It can be noted that SHIx levels vary from river to river, though the frequent levels used are (iiia) and (i). Generally, it was observed that at the Q95 level, the appropriate MC order was found to be MC0 and from the Q90 to the Q75 levels, MC1 (with few exceptions at the Q90 level). Beyond the Q75 level, mostly MC1 did not seem to work, probably because in the higher range of Qx the LT is better simulated by MC2, as demonstrated by Sharma and Panu [6].
The values of Lcr and LT-e′ so computed are shown in column 5 (Table 3). It can be seen that Lcr and LT-e′ display a large divergence between themselves and none of them matched to LT-o, shown in bold numeral in column 2. So, the first thought that came to mind was to compute the average of Lcr and LT-e′ [option (a)] as an estimator of LT-e and compare it with LT-o. The LT-e values so obtained are shown in column 6, which shows some correspondence with LT-o values presented in column 2.
A plot (Figure 3A) between LT-o and LT-e based on option (a) was drawn on a 1:1 line, and statistics NSE and MER were evaluated. The plot resulted in NSE ≈ 84% but with a significant under-prediction (MER = −8.86%). To ameliorate the under-prediction, option (b), i.e., LT-e was computed as LT-e (=0.75 Lcr + 0.25 LT-e′). This estimate of LT-e using option (b) proved less meaningful as NSE dropped to ≈ 80% accompanied by a significant over-prediction (MER ≈ 9%), as is shown in Figure 3B. Therefore, a more compromising estimate of LT-e, option (c) was considered as LT-e (=0.60 Lcr + 0.40 LT-e′). This new estimate of LT-e turned out better with an MER ≈ −1%, with an acceptable NSE ≈ 83% as is evinced by the plot in Figure 4. The values of LT-e so obtained are shown in column 9 in bold letters. In short, the weighted average of LCR and LT-e′ compared quite well with the LT-o in the process of model fitting.

4.2. Validation of the Model Structure

To investigate the adequacy of the model, the LT-e values for the four rivers (Table 1, italicized identification) were estimated and compared with the observed (LT-o) counterparts, and they are shown in Table 4. The validation affirmed that the modeling worked with the weighted averaging option (c) to be the best one. In this type of averaging, the relative difference between LT-e and LT-o was found to vary from −15 to 21%, with a mean deviation of 3% and a standard error of 10% (column 9, Table 4). The simple averaging tended to be next in line with mean deviation = −4% with the standard error of 11% (column 5, Table 4). The weighted averaging of option (b) was less satisfactory, as the mean deviation turned out to be excessively large = 12%, with the standard error of 12% (column 7, Table 4). Based on the foregoing discussion, one can easily conjure that the suggested procedure and model in this paper are robust for providing dependable performance in estimating drought magnitudes and lengths.
For weekly flow sequences at the uniform cutoff level, it may be worthwhile to compare the present methodology with that used by the authors [7]. In their paper, a traditional approach was used to first estimate the LT, using weekly SHI sequences by truncating at a pulsating SHIx, with two cutoff levels, viz. Q90 and Q95. The MT was predicted by using the relationship (MT = µd × LT). The NSE and the MER for this kind of formulation turned out to be low (NSE = 72% for the LT and 76% for MT). In specific terms, there was acute under-prediction for the high end of MT values. Using a different approach, the present paper proposes a predictive model with five cutoff levels (Q95 to Q75), and it has yielded predicted values of MT (i.e., MT-e) almost the same as the observed values (i.e., MT-o). Likewise, the NSE for LT turned out to be 83%, which falls within the acceptable range.
In summary, the model fitting and validation with 24 rivers revealed that MT-e can be simulated adequately by taking it equal to VR′, which in turn allowed the determination of the length of the critical period (Lcr) via behavior analysis. In other words, at low uniform cutoff levels from Q95 up to Q75, a combination of Lcr with MC based LT-e′ on SHI sequences can be used to estimate LT-e.

4.3. An Illustrative Example for the Estimation of MT-e and LT-e

Consider the Neebing River (near Thunder Bay), northern Ontario, Canada for which one needs to estimate the drought duration (weeks) and the magnitude (m3) for a return period of 65 years (T = 3380 weeks = 52 × 65), with the following statistics: μo = 1.62 m3/s; σo = 3.04 m3/s; σav = 1.79 m3/s; and ρ = 0.63. A severe drought is defined when flow drops below the Q90 (= 0.089 m3/s) level. One may investigate the additional requirement of water when the river has already receded below the Q90 level. The drought probabilities from the historical weekly data were computed as q = 0.10 (Q90), qp = 0.050, and qq = 0.654. The drought lengths were adequately simulated by MC0. These parameters can also be estimated using the equations involving ρ and a normal pdf of the SHI sequences [6].
A value of a 65-year drought duration LT-e′ at the Q90 level can be computed to be approximately = 4 weeks (MC0), using Equation (1) and plugging in values of relevant parameters (T = 3380, q = 0.10, qp = 0.10, qq = 0.10 because of MC0). Based on the historical data, behavior analysis yields Lcr = 36 weeks. Therefore, LT-e (=0.60 × 36 + 0.40 × 4 =) is 23 weeks. MT-e was found to be = 0.58. When MT-e was converted to the total deficit volume (volumetric units) DT-e = MT-e × σav = 0.58 × 1.79 × (7 × 24 × 3600) m3 = 0.63 × 106 m3 per week (note that the duration is in weeks and hence flow is to be converted on a weekly basis that brings in the multiplier 7 × 24 × 3600). Since drought is lasting for 23 weeks computed above, the total deficit volume at the Q90 level = 0.63 × 106 × 23 m3 ≈ 14.50 × 106 m3. That means only 14.50×106 m3 is available in the river, because the river is in the grip of a severe drought.
If the cutoff level for the drought-free condition is defined as Q75, then similar calculations (as done at Q90) will result in LT-e = 35 weeks (Lcr = 48 and LT-e′ = 16) and MT-e = 2.16. At Q75, the drought lengths fell in the regime of MC1 (with q = 0.25, qq = 0.70, and qp = 0.12). Therefore, the deficit volume would amount to 2.16 × 1.79 × 7 × 24 × 3600 × 35 = 81.84 × 106 m3. This volume of water met all water demands from the river, when the drought was not there.
To cater for needs during severe drought (Q90 level), we need the additional water to shore up to the drought-free level of Q75. Thus, the additional volume of water needed is 81.84 − 14.50 = 67.34 × 106 m3, which must be arranged to restore the health of the river and to assuage other demands that the river has been coping with in the past. This could be met from external sources by diverting the water from reservoirs or by pumping from the lakes or similar water bodies. It can be appreciated that under the severe drought conditions (Q90), to restore the river flows to the modest drought free level (Q75), nearly 4.6 times more water is needed, i.e., (67.34/14.50). Likewise, to deal with extreme drought, at the Q95 level, the requirement would further shore up and would work out = (81.84 − 6.28) = 75.56×106 m3, (i.e., 12 times = 75.56/6.28). At the Q95 level, calculations can be done with the MC0 regime of drought lengths with Lcr = 15, LT-e = 10, and MT-e = 0.58, and the available volume will work out to be 6.28 × 106 m3.

4.4. A Comment on the Present State of Drought Magnitude Assessment and Cutoff Levels at the Uniform Flow Levels

The truncation level approach in assessing hydrological droughts originated in the 1960’s [11], which used the mean or the median of the flows of the respective seasons (months) as the truncation levels. In doing so, the non-stationary series of drought variables is rendered stationary in the statistical sense. The analysis of stationary series is amenable to statistical analysis, using the theory of runs and Markovian processes. The stationary series, when standardized is tantamount to the standardized precipitation index series, SPI [13,14] used in the context of meteorological drought. The parameters of drought, viz. duration, magnitude and intensity, have been synthesized using the aforesaid variable cutoff levels (on a monthly or seasonal basis) and the results have been deemed satisfactory for drought mitigation purposes. One outcome of such analyses has catapulted in identifying interactions between meteorological and hydrological drought characteristics [27]. The other outcome has resulted in the estimation of drought deficit volume as an aid in planning and in evaluating the performance of the water storage systems [15,16,17,27,28].
A major challenge of drought assessment in the statistical sense is spurred by the desire for uniformity of the cutoff level throughout the flow regime, such as Q90, Q70, etc. In other words, if the demand level of water is relatively constant irrespective of seasons, then the analysis using the stochastic theory applicable to stationary time series turns out to be tedious. The current approach for handling such a problem on daily flow sequences is confined to the frequency distribution fitting of durations and magnitudes empirically abstracted from the historical data by truncating at the desired level [1,2,3,4]. In such an analysis, the linkage between the parameters of frequency functions of duration and magnitude, and that of the drought variable (such as streamflow series) is still to be defined. Ideally, the endeavor should be directed to the assessment of drought parameters based on the readily available statistics (μ, cv, and ρ) of the drought variable at the uniform cutoff level. This paper has demonstrated some success in meeting the above goal in which drought parameters are being predicted using the estimates of probabilities derived from the data. Most importantly, the estimates of reservoir volume and the length of the critical period (Lcr) have been used as key elements in the analysis.

5. Conclusions

Drought parameters, viz. MT-e and LT-e on a weekly basis at the uniform truncation level, such as from Q95 to Q75 flow levels, can be estimated from the analysis of weekly flow sequences using behavior analysis used in the estimation of reservoir volumes. The behavior analysis provided an estimate of reservoir volume, VR′ (≈MT-o) and that of Lcr. The Markov chain models (of order 1: MC1) and (of zero order: MC0) performed satisfactorily in simulating LT-e′ by using the estimates of simple and conditional drought probabilities from weekly SHI sequences. The most important element in estimating the probabilities was the appropriate identification of the cutoff level (SHIx). The frequently occurring cutoff levels were found to be [(Qx − µo)/σav], [(Qx − µo)/σo] and the average of these two entities. Using the hypothesis that drought magnitude = drought intensity × drought length, the scaling parameter Φ for the characteristic drought length was evaluated to fine-tune the estimates of MT-e to be equal to VR. The LT-e′ from the MC analysis when combined with Lcr from the behavior analysis (i.e., LT-e = 0.60Lcr + 0.40LT-e′) yielded the drought length LT-e, which was found compatible with the observed drought length, LT-o.

Author Contributions

T.C.S. and U.S.P. collaboratively accomplished various tasks related to this manuscript, such as draft preparation and the finalization of the manuscript for journal publication, including the conceptualization and the development of the methodology; data collection; and data analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The partial financial support of the Natural Sciences and Engineering Research Council of Canada for this paper is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zelenhasic, E.; Salvai, A. A method of streamflow drought analysis. Water Resour. Res. 1987, 23, 156–168. [Google Scholar] [CrossRef]
  2. Tallaksen, L.M.; Madsen, H.; Clausen, B. On definition and modeling of streamflow drought duration and deficit volume. Hydrol. Sci. J. 1997, 42, 15–33. [Google Scholar] [CrossRef]
  3. Tallaksen, L.M.; Van Lanen, H.A. Hydrological Drought: Processes and Estimation Methods for Streamflow and Groundwater, Developments in Water Science; Elsevier: Amsterdam, The Netherlands, 2004; Volume 48. [Google Scholar]
  4. Laaha, G. Hydrological drought severity explained by climate and catchment characteristics. J. Hydrol. 2015, 526, 3–14. [Google Scholar] [CrossRef] [Green Version]
  5. Clausen, B.; Pearson, C.P. Regional frequency analysis of annual maximum streamflow drought. J. Hydrol. 1995, 173, 111–130. [Google Scholar] [CrossRef]
  6. Sharma, T.C.; Panu, U.S. Analytical procedures for weekly hydrological droughts: A case of Canadian rivers. Hydrol. Sci. J. 2010, 55, 79–92. [Google Scholar] [CrossRef]
  7. Sharma, T.C.; Panu, U.S. Predicting drought durations and magnitudes at weekly time scale: Constant flow as a truncation level. Open J. Atmos. Clim. Chang. 2014, 1, 1–16. [Google Scholar] [CrossRef]
  8. TCLFE (Task Committee on Low Flow Evaluation). Characteristics of low flows. ASCE J. Hydraulic Div. 1980, 106, 717–731. [Google Scholar] [CrossRef]
  9. Gustard, A.; Bullock, A.; Dixon, J.M. Low Flow Estimation in the United Kingdom; Report no.108; Institute of Hydrology: Wallingford, UK, 1992. [Google Scholar]
  10. Dracup, J.A.; Lee, K.S.; Paulson, E.G., Jr. On identification of droughts. Water Resour. Res. 1980, 16, 289–296. [Google Scholar] [CrossRef]
  11. Yevjevich, V. Methods for determining statistical properties of droughts. In Coping with Droughts; Yevjevich, V., da Cunha, L., Vlachos, E., Eds.; Water Resources Publications: Littleton, CO, USA, 1983; pp. 22–43. [Google Scholar]
  12. Kjeldsen, T.R.; Lundorf, A.; Rosbjerg, D. Use of a two-component exponential distribution in partial duration modeling of hydrological droughts in Zimbabwean rivers. Hydrol. Sci. J. 2000, 45, 285–298. [Google Scholar] [CrossRef]
  13. McKee, T.B.; Doesen, N.J.; Kleist, J. The relationship of drought frequency and duration to time scales. In Proceedings of the 8th Conference on Applied Climatology, Anaheim, CA, USA, 17–22 January 1993; pp. 179–184. [Google Scholar]
  14. Hayes, M. Drought Indices. Available online: www.drought.unl.edu/whatis/indices.pdf.lastaccesed (accessed on 15 May 2022).
  15. Wu, J.; Chen, X.; Love, C.A.; Yao, H.; Chen, X.; AghaKouchak, A. Determination of water required to recover from hydrological drought: Perspective for drought propagation and non-standardized indices. J. Hydrol. 2020, 590, 125227. [Google Scholar] [CrossRef]
  16. Borgomeo, E.; Pflug, G.; Hall, J.W.; Hochrainer-Stigler, S. Assessing water resources system vulnerability to understand hydrological drought using copulas to characterize drought duration and deficit. Water Resour. Res. 2015, 51, 8927–8948. [Google Scholar] [CrossRef] [Green Version]
  17. Rad, A.M.; Ghahraman, B.; Khalili, D.; Ghabremani, Z.; Ardakani, S.A. Integrated meteorological and hydrological drought model: A management tool for proactive water resources planning in semi-arid regions. Adv. Water Resour. 2017, 107, 336–353. [Google Scholar] [CrossRef]
  18. McMahon, T.A.; Adeloye, A.J. Water Resources Yield; Water Resources Publications: Littleton, CO, USA, 2005; pp. 80–88. [Google Scholar]
  19. McMahon, T.A.; Pegram, G.G.S.; Vogel, R.M.; Peel, M.C. Revisiting Reservoir Storage–Yield Relationships using a Global Streamflow Database. Adv. Water Resour. 2007, 30, 1858–1872. [Google Scholar] [CrossRef]
  20. Sharma, T.C.; Panu, U.S. Reservoir sizing at the draft level of 75% of mean annual flow using drought magnitude method on Canadian rivers. Hydrology 2021, 8, 79. [Google Scholar] [CrossRef]
  21. Sharma, T.C.; Panu, U.S. Compatibility of drought magnitude based method with SPA for assessing reservoir volume: Analysis using Canadian river flows. J. Water Resour. Protection 2022, 14, 1–20. [Google Scholar] [CrossRef]
  22. Adamowski, K. Plotting formula for flood frequency. Water Resour. Bull. 1981, 17, 197–202. [Google Scholar] [CrossRef]
  23. Sen, Z. Statistical analysis of hydrological critical droughts. ASCE J. Hydraul. Eng. 1980, 106, 99–115. [Google Scholar]
  24. Sen, Z. Applied Drought Modeling, Prediction and Mitigation; Elsevier Inc.: Amsterdam, The Netherlands, 2015. [Google Scholar]
  25. Chow, V.T.; Maidment, D.R.; Mays, L.W. Applied Hydrology; McGraw-Hill: New York, NY, USA, 1988; p. 367. [Google Scholar]
  26. Environment Canada. HYDAT CD-ROM Version 96-1.04 and HYDAT CD-ROM User’s Manual. In Surface Water and Sediment Data; Water Survey of Canada: Gatineau, QC, Canada, 2020. [Google Scholar]
  27. Lorenzo-Lacruz, J.; Vicente-Serrano, S.M.; Gonzalez-Hidalgo, J.C.; Lopez-Moreno, J.I.; Cortesi, N. Hydrological drought response to meteorological drought in the Iberian Peninsula. Clim. Res. 2013, 58, 117–131. [Google Scholar] [CrossRef] [Green Version]
  28. Lopez-Moreno, J.I.; Vicente-Serrano, S.M.; Beguria, S.; Garcia-Ruiz, J.M.; Portela, M.M.; Almeida, A.B. Dam effect on drought magnitude and duration in a transboundary basin: The lower River Tagus, Spain and Portugal. Water Resour. Res. 2009, 45, W02405. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Map of Canada showing hydrometric gauging stations (source: Environment Canada).
Figure 1. Map of Canada showing hydrometric gauging stations (source: Environment Canada).
Hydrology 09 00109 g001
Figure 2. Flow chart of calculations to compute MT-e and LT-e.
Figure 2. Flow chart of calculations to compute MT-e and LT-e.
Hydrology 09 00109 g002
Figure 3. Comparison of estimated (LT-e) and observed (LT-o) drought lengths.
Figure 3. Comparison of estimated (LT-e) and observed (LT-o) drought lengths.
Hydrology 09 00109 g003
Figure 4. Comparison of estimated (LT-e) and observed (LT-o) drought lengths.
Figure 4. Comparison of estimated (LT-e) and observed (LT-o) drought lengths.
Hydrology 09 00109 g004
Table 1. Summary of statistical properties of weekly flows of the rivers under consideration.
Table 1. Summary of statistical properties of weekly flows of the rivers under consideration.
Numeric Identifier of the River in Figure 1 with Name and Gauging Station IdentityData Size (Years)Area (km2)μocvocvmxcvavcvgmρ
[1] Fraser at Shelley, BC08KB00170 (1951–2020)32,400817.340.900.820.360.280.75
[2] Athabasca River at Athabasca, AB07BE00169 (1952–2020)74,600429.190.981.320.430.260.81
[3] Bow River at Banff, AB05BB001110 (1911–2020)221039.241.111.220.310.140.72
[4] pipestone River at Karl lake, ON04DA00154 (1967–2020)596059.051.041.560.630.410.89
[5] Neebing at Thunder Bay, ON02AB00866 (1954–2019)1871.621.873.841.100.680.63
[6] Pic River near Marathon, ON02BB00350 (1971–2020)427050.101.242.160.710.480.74
[7] Pagwachaun at highway#11, ON04JD00553 (1968–2020)202053.081.452.770.790.480.74
[8] Nagamgami at highway#11, ON04JC00270 (1951–2020)241024.561.111.660.550.400.87
[9] Batchawana at Batchawana, ONBF00150 ((1971–2020)119022.381.382.750.740.520.62
[10] Goulis near Searchmont, ON02FB00253 (1968–2020)116018.371.322.690.750.550.69
[11] Whitson at Chemsford, ON02CF00760 (1961–2020)2433.061.503.620.780.570.68
[12] North French near Mouth, ON04MF00154 (1967–2020)119095.721.292.440.710.460.72
[13] Labase River at North Bay, ON02DD01354 (1975–2018)70.40.911.493.240.960.790.44
[14] Chippewa Creek at North Bay, ON02DD01454 (1975–2018)37.30.621.112.100.810.670.43
[15] Commanda at Commanda, ON02DD01546 (1975–2020)1061.761.222.310.770.660.58
[16] N. Magnetwan at Pickerel Lake, ON02EA01052 (1969–2020)1492.861.262.540.820.700.51
[17] Becancour A Lyster, QC02PL00146 (1923–1968)141030.621.322.460.820.690.62
[18] Beaurivage A. Sainte Entiene, QC02PJ00775 (1926–2000)70914.211.472.670.900.770.49
[19] Lepreau River at Lepreau, NB01AQ001101 (1919–2019)2397.431.082.010.870.800.49
[20] Carruther at Saint Anthony, PE01CA00359 (1962–2020)46.80.971.332.890.820.620.48
[21] Bevearbank River at Kinsac, NS01DG00388 (1922–2019)973.041.091.530.900.840.43
[22] N. Margaree at Margaree valley, NS01FB00190 (1929–2020)36817.010.961.600.680.590.46
[23] Upper Humber at Reidville, NF02YL00168 ((1953–2020)221080.291.071.570.660.600.48
[24] Torrent River at Bristol Pool, NF02YC00161 (1960–2020)62424.901.072.210.620.530.57
Table 2. Summary of calculations in the identification of the cutoff level for the estimation of drought length (LT-e), Upper Humber River, Newfoundland, Canada.
Table 2. Summary of calculations in the identification of the cutoff level for the estimation of drought length (LT-e), Upper Humber River, Newfoundland, Canada.
River IdentityQx, (q), VR′, LcrCutoff
SHIx = Z0
q1 *, qq, qpLT-e′, µd, MT-eDecisionModel
Φ, MT-e
LT-e
12345678
Upper Humber
(#23, Table 1)
T = 3536
σav = 53.10
ρ = 0.48
Q75 = 24.16, (q = 0.25), 4.18, 21
LT-o = 17
(iii) = −1. 06
(iiia) = −0.86
(i) = −0.66
0.046, 0.389, 0.030
0.117, 0.528, 0.063
0.242, 0.575, 0.136
6, 0.51, 2.84
9, 0.55, 4.87
11, 0.60, 6.68
MT-e < VR′, next
MT-e > VR′, q1
is low, next
cutoff (i) is fine
MC1
0.47, 4.20
A = 16
b = 19
c = 17
Q80 = 20.47, (q = 0.20), 3.06, 20
LT-o = 16
(iii) = −1.13
(iiia) = −0.91
(i) = −0.70
0.030, 0.358, 0.020
0.091, 0.497, 0.051
0.210, 0.552, 0.119
5, 0.50, 2.41
8, 0.54, 4.25
10, 0.59, 6.06
MT-e < VR′, next
MT-e > VR′, q1
is low, next
cutoff (i) fine
MC1
0.63, 3.07
a = 15
b = 18
c = 16
Q85 = 16.78 (q = 0.15), 2.02, 17
LT-o = 15
(i) = −0.740.185, 0.530, 0.107
At cutoff levels (iiia) and (iii) q = 0.071 and 0.02, very low
10, 0.58,5.54Cutoff(i) fine, Level (iii) and (iiia) were rejected due to low q1 valuesMC1
0.82, 2.01
a = 13
b = 15
c = 14
Q90 = 13.57 (q = 0.10), 1.27, 16
LT-o = 14
(i) = −0.780.158, 0.538, 0.087
At cutoff levels (iiia) and (iii) q = 0.056 and 0.014, very low.
9, 0.57, 5.42Cutoff Level (i) is fine. Cutoff (iiia) and (iii) resulted in very low values of q1, so were rejectedMC1
0.99, 1.27
a = 13
b = 14
c = 13
Q95 * = 10.48 (q = 0.05), 0.70, 15
LT-o = 9
(i) = −0.81
(iiia) = −1.06
0.138, 0.518, 0.077
0.043, 0.379, 0.028
4, 0.51, 2.20
(9, 0.56, 4.99)
3, 0.51, 1.27
(5, 0.51, 2.75)
Cutoff level (i) is too high—so next lower
cutoff level (iiia) is fine
MC0
0.97, 0.70
a = 9
b = 11
c = 10
Asterisk (*) denotes that at Q95 level drought length resembles MC0 because the values of VR′ better correspond to the MC0 based values of MT-e′. The values in italics show the relevant values with MC1 based lengths. However, from Q90 to Q75, MC1 yielded MT-e′ values in better correspondence to VR′. Note that all entities have no units except Qx, which is in m3/s, and LT-e and LT-e′ are in weeks.
Table 3. Summary of the results based on the calculations for estimated values, LT-e and MT-e at varying cutoff levels.
Table 3. Summary of the results based on the calculations for estimated values, LT-e and MT-e at varying cutoff levels.
River IdentityVRMT-oLT-oQx *q1 *qqqpCutoffMT-eΦMT-eLcrLT-eMC typeLT-e **
abc
123456
Athabasca (#2, Table 1)
T = 3588
σav = 183.97
ρ = 0.81
0.930.931569.730.0220.5000.011iiia0.930.000.93202MC0111613
1.991.901979.770.0250.5380.012iiia3.000.501.992417MC1152017
3.373.372191.920.2020.7270.069i3.580.093.362817MC1182320
4.564.5622102.320.2110.7280.073i9.580.674.572817MC1222524
5.705.7023112.350.2240.7380.076i10.130.565.692918MC1232624
Goulis
(#11, Table 1)
T = 2756
σav = 13.84
ρ = 0.69
1.501.50312.2040.0790.6470.031iiia5.921.001.533411MC1222825
4.104.10343.2880.0950.6540.037iiia6.390.484.093712MC1243127
5.905.90354.0100.2590.7700.081i12.920.685.863920MC1303432
7.367.36354.5870.2750.7820.083i13.750.597.333922MC1313532
8.898.89355.1900.2920.7850.089i13.970.478.863922MC1313532
Bevearbank
(#21, Table 1)
T = 5096
σav = 2.74
ρ = 0.43
0.240.24120.0650.0650.4130.041i1.551.000.58163MC091311
0.800.80190.1740.1450.5150.082ia2.310.930.79214MC0121714
1.791.79190.3240.1680.5180.097ia5.380.871.77249MC1172018
3.203.20210.5110.2030.5410.117ia5.970.593.222710MC1182320
4.634.63220.6960.2550.5660.148ia6.710.404.593411MC1222825
Asterisk (*) indicates Qx values at Q95, Q90, Q85, Q80, and Q75 levels for each river in column 1. The values of q (=qx) are 0.05, 0.10, 0.15, 0.20, and 0.25, respectively, based on the Qx flow levels (flow duration curves). Asterisk (**) Indicates LT-e values (weeks) are computed for 3 options (a), (b), and (c) discussed in the text.
Table 4. Summary of the results in the validation process in terms of estimated values, LT-e against LT-o at varying cutoff levels.
Table 4. Summary of the results in the validation process in terms of estimated values, LT-e against LT-o at varying cutoff levels.
River IdentityQx *VRMT-0LcrLT-0CutoffΦMT-eModel orderLT-eLT-eLT-e
(a) Ϯ(%) ҂(b) ϮϮ(%)҂(c) ϮϮϮ(%) ҂
123456789
Lepreau River (#19, Table 1)
T = 3588
σav = 183.97, ρ = 0.49
0.580.780.571812iii0.780.78MC010−16.671416.67120.00
1.031.781.782415iiia0.551.78MC1150.001926.671713.33
1.402.702.702719i0.562.72MC118−5.262215.79190.00
1.783.753.683219i0.423.78MC12110.532636.842321.05
2.195.175.173423i0.275.20MC122−4.352821.74258.70
Pipestone River
(#4,Table 1), T = 2808
σav = 37.31
ρ = 0.89
9.200.910.911613iiia0.600.91MC09−30.77130.0011−15.38
10.841.521.521716iiia0.021.52MC11718.751718.751718.75
12.672.282.282018iiia1.002.54MC1195.56195.56195.56
14.463.083.082220iiia0.933.04MC1200.00215.00200.00
16.664.164.162422i0.964.22MC1260.00254.55250.00
Chippewa River
(#14, Table 1)
T = 2288
σav = 0.48,
ρ = 0.43
0.120.450.4397iiia10.59MC06−14.2970.0070.00
0.140.860.721110iiia10.98MC19−10.00100.009−10.00
0.161.231.231611iiia0.951.24MC1129.091427.271318.18
0.191.621.621813iiia0.871.63MC1130.001623.08147.69
0.212.272.272014iiia0.742.27MC1157.141721.431614.29
Bow River
(#3, Table 1)
T = 5720
σav = 12.05
ρ = 0.72
7.111.771.772018i0.991.77MC117−5.56195.56180.00
7.672.672.672220i0.862.69MC118−10.00200.0019−5.00
8.183.553.552421i0.743.56MC119−9.52224.7620−4.47
8.634.344.342622i0.634.37MC120−9.09234.5521−4.55
9.165.275.272624i0.505.28MC120−10.6723−4.1721−12.56
Mean
Stan. error
−4.0%
11.4%
12.0%
11.9%
3%
10.4%
Asterisk (*) indicates Qx values at Q95, Q90, Q85, Q80, and Q75 levels in 5 rows for each river in column 1. (a) Ϯ indicates option-a, (b) ϮϮ indicates option-b, (c) ϮϮϮ indicates option-c, and (%) ҂ indicates percent deviation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sharma, T.C.; Panu, U.S. A Procedure for Estimating Drought Duration and Magnitude at the Uniform Cutoff Level of Streamflow: A Case of the Weekly Flows of Canadian Rivers. Hydrology 2022, 9, 109. https://0-doi-org.brum.beds.ac.uk/10.3390/hydrology9060109

AMA Style

Sharma TC, Panu US. A Procedure for Estimating Drought Duration and Magnitude at the Uniform Cutoff Level of Streamflow: A Case of the Weekly Flows of Canadian Rivers. Hydrology. 2022; 9(6):109. https://0-doi-org.brum.beds.ac.uk/10.3390/hydrology9060109

Chicago/Turabian Style

Sharma, Tribeni C., and Umed S. Panu. 2022. "A Procedure for Estimating Drought Duration and Magnitude at the Uniform Cutoff Level of Streamflow: A Case of the Weekly Flows of Canadian Rivers" Hydrology 9, no. 6: 109. https://0-doi-org.brum.beds.ac.uk/10.3390/hydrology9060109

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop