2.1. Study Site and Data Sources
The Manouane River is sourced in its namesake lake and flows into the Peribonka River, which is the main tributary of Lake Saint-Jean. The Manouane River watershed is entirely comprised within the Canadian Shield, covering 9483 km
2 (
Figure 1). Two water diversion structures were built on the Manouane River: one in 1961 by Alcan, to feed into Peribonka Lake, and the second in 2003 by Hydro-Québec, to increase the capacity of the Pipmuacan reservoir, which supplies the Bersimis-1 and Bersimis-2 hydroelectric plants on the Betsiamites River.
The diversion dam built by Hydro-Québec on the Manouane River is 9 m high and 90 m long. The maximum surface covered by the related reservoir is 21 km2. Water diverted from the Manouane River to the Betsiamites River is first routed through a 7 km long channel then through the des Hirondelles Brook. Mean diverted flow is 30 m3/s and maximum flow through the channel is 80 m3/s. At the level of the dam, the Manouane River watershed covers 1717 km2 in surface area.
The station at which Manouane River flows are measured is located downstream from the diversion dam. Flow measurements at that station started towards the end of 1999 and continue uninterrupted to this day. The data were taken from the website of [
22]. Between the station and the dam, the main tributary of the Manouane River, the Petite Manouane River (1445 km
2) flows into it as well as several small tributaries, such that the watershed surface area at the level of the measurement station is 3600 km
2. It should be kept in mind that the effects of the diversion dam on flows measured at the station are attenuated by inputs from these natural tributaries. This study only looks at the impacts produced by the latter diversion because no flow data are available for the period prior to the first diversion in 1961. It should also be pointed out that there are no temperature and precipitation measuring stations in the Manouane River watershed. The meteorological station of reference for this watershed is at the Bagotville military airport (48°20′ N; 71°00′ W).
The Manouane River was selected because of the availability of flow measurement data collected before and after construction of the diversion dam. Flow measurements began in late 1979 and are ongoing to this day. They therefore cover the periods both before (1980–2002) and after (2004–2014) construction of the dam.
2.2. Statistical Analysis
Three methods were used to analyze hydrologic impacts induced by the diversion dam: the monitoring method, which is based on a comparison of mean values of flows before (1980–2002) and after (2004–2014) dam construction using the non-parametric Kruskal-Wallis method; the analysis of the long-term trend of flows using the classic non-parametric Mann-Kendall method; and the non-parametric Lombard method, which allows the detection of shifts in mean values. The last two methods, which are briefly described below, were applied over the whole period of interest (1980–2014).
The Mann-Kendall test is widely used in hydrology. Given a sample
of values that are independent from a random variable
X for which the stationarity or long-term trend must be assessed, the Mann-Kendall statistic is defined as follows:
where
Xi and
Xj are sequential values of
X and n is the sample size. The test statistic is obtained by counting, for each
pair, the number of cases where the second value is greater than the first, and the number of cases where the second value is less than the first, then subtracting these two numbers. The presence of a statistically significant trend is assessed using the
Z score value as follows:
A positive (negative) Z score reflects an increasing (a decreasing) long-term trend, and its significance is compared with the critical value or significance threshold for the test. The critical Z score values when using a 95% confidence level are −1.96 and +1.96 standard deviations. The p-value associated with a 95% confidence level is 0.05. If Z score is between −1.96 and +1.96, p-value will be larger than 0.05, null hypothesis cannot be rejected.
The Lombard method can be used to detect the nature and timing of shifts in the mean and/or variance of a statistical series. Unlike other methods commonly used in hydrology (e.g., the Pettitt method), the Lombard method can distinguish between sharp and gradual shifts in mean and variance. It is therefore a more general method than the other methods, which only detect sharp changes. The mathematical basis for the Lombard method has been described in detail by [
23] and [
24]. In the present study, only the most important statistical aspects of the methods are described (e.g., [
25,
26]). Given a series of independent observations
where
Xi is the observation taken at time
A question that is often of interest is to assess whether the mean of this series has changed at some unknown time. To this end, one considers as a possible pattern for the mean of these observations the smooth-change model introduced by [
23], where the mean of
Xi is defined by
In other words, the mean changes gradually from to between the times and . As a special case, one has the usual abrupt-change model when
In order to test formally that the mean of a univariate series is stable, or on the contrary that is follows model (3), one can use the statistical procedure by [
23]. To this end, let
denote the rank of
among
X1,…,
Xn and define the rank score of
Xi by
where
is Wilcoxon’s score function, while
Lombard’s test statistic is
where
A value of
Sn greater than 0.0403 derived for a series of observations indicates that there is a shift in the mean value of this series at the 5% probability level. This value of 0.0403 corresponds to the asymptotic theoretical (critical) value as obtained by [
23]. Note that the test proposed by [
23] in order to detect multiple abrupt changes in the mean was also performed; here again, the test concluded to a significant change in the mean. It is important to note that the independence assumption among the observations is necessary for the validity Lombard’s test [
23,
24].
The rate of variation
R of the mean values before and after shifts detected using the Lombard method was derived using the following formula:
where
M1 and
M2 are, respectively, the mean values of flow before and after the shift,
R being expressed as a percentage.
The three methods were applied to the following two types of series for the period from 1980 to 2014:
The series of annual maximum and minimum daily flows, which consist of the highest (maximum) and lowest (minimum) daily flow values measured each year (October to September) from 1980 to 2014;
The series of seasonal maximum and minimum daily flows, which consist of the highest (maximum) and lowest (minimum) daily flow values measured each season (winter: January to March; spring: April to June; summer: July to September; winter: October to December) from 1980 to 2014.