assessment of hydrological drought in northern ontario using standardized streamflow index
abstract
momentum is gaining for the use of a standardized streamflow index (ssi) for the assessment
of hydrological drought, which is an extension of the well-established and popularly used
meteorological drought index, standardized precipitation index (mckee et al, (1993). drought
characterization using ssi consists of various data treatments to transform streamflow data to
standard normal z-score values, where the mean is zero and the standard deviation is one. the
underlying assumption when using z-score in early warning systems and drought monitoring
programs is that the scores associated probabilities may be represented by the normal (gaussian)
distribution. applying data treatments that fail to achieve the conditions of normal distribution
results inaccurate ssi drought assessments. there is an opportunity in northern ontario to
conduct ssi drought assessments in early warning systems and drought monitoring programs
using the abundance of long and continuous streamflow records; however, widespread
applications would only be achieved if the resulting assessments are accurate, equitably and
reliable. the data treatments investigated in the thesis include: untreated and treated (i.e., log
normal transformation, and the fitted gamma distribution).
monthly assessments for a total of 40 rivers from across northern ontario with an
average record length of ≈ 45 years were utilized in the analysis. historical droughts sample
periods 1976/77 and 1998 were used in evaluations as they were determined to pertain to the
well-defined regional drought events with confirmation of historically significant impacts in
addition to containing seasonal influences. the log normal transformation followed by the
fitted gamma distribution exhibited the most consistency in not rejecting normality assumption
of the data set when using the shapiro wilks and the anderson darling tests. these tests also
show the best performance of the fitted gamma distribution for the winter season (november to
may) and of the log normal transformation for the summer season (june to october). using
these two data treatments for the respective seasons permitted the assumption of normality to be
rejected 15.7% of the time. the use of best performing data treatments and removal of data for
the month of march further reduced the assumption of normality being rejected only 12.0% of
the time. in evaluating the frequency of occurrence for the severe and extreme z-score values (-
1.64 and -2.33), the fitted gamma distribution followed by the log normal transformation
demonstrated the most acceptable scoring distribution. untreated data performed less efficiently.
all data treatments underperformed for the expected frequency at the extreme z-score of -2.33
(i.e., the extreme drought).
the threshold level q80 was determined to be the most appropriate trigger of monthly
assessments for northern ontario because it was found to equitably delineate severe regional
drought perceived by the respectively impacted sectors. the threshold level q80 represents these
sectors demands by identifying the point where streamflows are exceeded 80% of time. since it
is representative of the impacted sectors; it is used to assess the performance of ssi. significant
differences were recorded when comparing q80 to the typical moderate drought classification
with a z-score equal to - 1. untreated data tacitly refers to the assumption of the normal distribution of
the monthly flows and correspondingly of ssi values. thus, the analysis based on the assumption of the
normality of the monthly ssi data turned out to be less reliable. all the data treatments underperformed
for the expected frequency at extreme z-score of -2.33. however, when applying q80 to the
theoretical equivalent z-score of - 0.84; identical results were achieved 76% and 66% of the
time, respectively for the log normal transformation and the fitted gamma distribution. for
the case of q80, the average difference for the drought initiation, termination, and duration were found to range approximately 1.3 months with the log normal distribution to 3.0 months for the
lesser performing data treatments, such as the untreated data.