Next Article in Journal
Development of Boosted Machine Learning Models for Estimating Daily Reference Evapotranspiration and Comparison with Empirical Approaches
Previous Article in Journal
Introducing a Calculator for the Environmental and Financial Potential of Drain Water Heat Recovery in Commercial Kitchens
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impacts of Agricultural Areas on Spatio-Temporal Variability of Daily Minimum Extreme Flows during the Transitional Seasons (Spring and Fall) in Southern Quebec

1
Department of Environmental Sciences and Research Centre for Watershed-Aquatic Ecosystem Interactions (RIVE, UQTR), University of Quebec at Trois-Rivières, 3351 Boulevard des Forges, Trois-Rivières, QC G9A 5H7, Canada
2
Higher National School of Hydraulics, Blida 09000, Algeria
3
Centre d’Applications et de Recherches en Télédétection (CARTEL), Université de Sherbrooke, Sherbrooke, QC J1K 2R1, Canada
*
Author to whom correspondence should be addressed.
Water 2021, 13(24), 3487; https://doi.org/10.3390/w13243487
Submission received: 29 October 2021 / Revised: 24 November 2021 / Accepted: 2 December 2021 / Published: 7 December 2021
(This article belongs to the Section Hydrology)

Abstract

:
Several statistical methods were used to analyze the spatio-temporal variability of daily minimum extreme flows (DMEF) in 17 watersheds—divided into three homogenous hydroclimatic regions of southern Quebec—during the transitional seasons (spring and fall), during the 1930–2019 period. Regarding spatial variability, there was a clear difference between the south and north shores of the St. Lawrence River, south of 47° N. DMEF were lower in the more agricultural watersheds on the south shore during transitional seasons compared to those on the north shore. A correlation analysis showed that this difference in flows was mainly due to more agricultural areas ((larger area (>20%) on the south than on the north shore (<5%)). An analysis of the long-term trend of these flows showed that the DMEF of south-shore rivers have increased significantly since the 1960s, during the fall (October to December), due to an increase in rainfall and a reduction in cultivated land, which increased the infiltration in the region. Although there was little difference between the two shores in the spring (April to June), we observed a decrease in minimum extreme flows in half (50%) of the south-shore rivers located north of 47° N.

1. Introduction

In cold temperate regions, the impacts of global warming on minimum river flows are relatively complex. This complexity results from the different impacts of warming on the two types of river flow feeding regimes: snowmelt and rainfall. Many studies have already shown a significant decrease in snowfall in many regions, in the current context of global warming [1,2,3,4,5]. As for rainfall, the impacts of warming generally result in increased rainfall in winter (increase in temperature) and in summer (increase in evapotranspiration) [6,7,8]. Despite numerous studies on the temporal variability of minimum flows in the context of global warming [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26], few have focused on analyzing the impacts of these changes in seasonal precipitation regimes on minimum flows. Furthermore, the impacts of global warming on minimum flows can be amplified or mitigated by land use, particularly agriculture, as observed in a number of watersheds in certain agricultural regions of the USA [13,27].
Like other cold temperate regions, the climate in Quebec regions consists primarily of two seasons: winter (cold and snowy) and summer (hot and/or rainy), and these longer seasons are bisected by two transitional seasons, spring and fall. Spring marks the transition from cold (winter) to hot (summer), while fall marks the opposite transition from hot (summer) to cold (winter). From a hydrological perspective, the sequence of the seasons strongly influences river flow regimes during the two transitional seasons. Spring flows are primarily influenced by snowmelt from the winter, while fall flows are primarily influenced by rainfall. As such, the factors affecting spatio-temporal flow variability during these two transitional seasons are not necessarily the same, and are not affected by global warming in the same way.
Studies in Quebec on the impacts of global warming on temperature and precipitation have shown that, in winter, these impacts mean rising temperatures and decreasing snowfall [28,29,30]. In summer, these impacts also mean a rise in the temperature, combined with increasing rainfall [31,32]. The hydrological impact of these changes in temperature and precipitation regimes were studied in winter and summer only [33,34,35,36]. Research conducted during the spring and fall has so far focused on the temporal variability of floods [37,38,39]. However, unlike the two previous seasons, these two transitional seasons are characterized by periods of high and low flows, in Quebec. In spring, low flows occur before or after snowmelt, and, in fall, after a period without rain. It is important to stress that these periods of low flow do not necessarily correspond to periods of hydrological drought. However, despite the existence of these periods of low flows, there are still no studies dedicated specifically to daily minimum extreme flows during the transitional seasons (spring and fall), despite their vital role in the function and evolution of aquatic ecosystems and human activity [40,41]. In Quebec, for example, certain indices (standards) for calculating the instream flows, to protect aquatic habitats, take these minimum flows into account in spring and in fall [42]. From a seasonal perspective, almost all previous research conducted in cold temperate regions analyzed minimum daily flows in winter and/or summer exclusively. There are no studies on minimum flows in the transitional seasons (spring and fall), and none that analyze the factors affecting these flows from a spatial perspective. Given these considerations, our study focused on the following three objectives:
  • Identify factors that influence spatial variability in daily minimum extreme flows in southern Quebec in spring (April to June) and fall (October to December). This objective is based on the hypothesis that, due to the differences in the mechanisms and processes affecting daily minimum extreme flows in the spring and fall, the factors affecting spatial variability also differ during these seasons in Quebec. It is important to remember that minimum extreme flows in spring are mainly associated with recession (surface runoff), whereas in fall they are generally associated with the direct recharge of bedrock aquifers during low-water events.
  • Compare the temporal variability, or long-term trend (stationarity of hydrological series), of daily minimum extreme flows in spring and fall. This objective is based on the hypothesis that daily minimum extreme flows decrease in spring, due to a decrease in snowfall in winter and early snowmelt, and in fall, due to increased evapotranspiration caused by higher temperatures in summer and early fall.
  • Finally, the 2010s is considered the wettest decade on record in North America. Given that this moisture caused several catastrophic floods in Quebec over the course of the decade, its impact on the evolution of daily minimum extreme flows during these seasons should be analyzed. We will test the hypothesis of whether this moisture caused a significant increase in daily minimum extreme flows in the spring and fall, in addition to flooding.

2. Materials and Methods

2.1. Description of the Watersheds and Data Sources

To test these three hypotheses, we analyzed the daily minimum extreme flows in 17 watersheds in spring and fall from 1930 to 2019. Human activity has very little impact on flow measurements. The 17 watersheds were divided into three hydroclimatic regions, defined in southern Quebec [33,37]: two on the south shore of the St. Lawrence River, and only one on the north shore (Figure 1 and Table 1). The eastern hydroclimatic region, located on the south shore, north of 47° N, is characterized by a cold temperate maritime climate, whereas the climate of the southeastern hydroclimatic region, just south-east of this parallel, is mixed temperate (part continental and part maritime). The north shore hydroclimatic region is characterized by a cold temperate continental climate. The rivers on the south shore run through the Appalachians and the St. Lawrence Lowlands, becoming more scarce heading northward. These two geological formations consist mainly of sedimentary rocks. The rivers on the north shore flow almost entirely on the Canadian Shield, which consists mainly of metamorphic rocks.
In hydrogeological terms, on the north shore of the St. Lawrence River, on the Canadian Shield, aquifers are located in cracks and faults in the bedrock, as well as in superficial deposits originating mainly from glaciers and rivers. On the south shore, they are found in crystalline, metamorphic, and sedimentary rock, as well as in lacustrine, glacial, marine, and river deposits covering these deep rocks. It is important to note that no studies have been done on the hydrogeological characteristics of these aquifers for southern Quebec as a whole. Existing studies only examined a few watersheds on both shores.
Physiographic variables were summarized per basin; they include mean watershed slope, drainage density, and land cover and use. Land cover and use (forest, wetland, agricultural and urban relative areas) variables were obtained from the 1.1 km resolution Canada land-cover database, derived from an advanced very high resolution radiometer (AVHRR) sensor operating on board the United States National Oceanic and Atmospheric Administration (NOAA) satellites. The classification was based on approximately 45 cloud free summer images, spanning the period 1988–1991.
As for watershed morphology, two key variables were measured: the mean slope of the watershed, and drainage density. These two variables were calculated using the 400 m resolution Canadian Digital Elevation Model (CDEM), and the drainage network produced by the Ministry of Energy and Natural Resources of Quebec (MERN). These two variables significantly influence the runoff and infiltration processes that impact minimum flows. Mean slope values of the watersheds analyzed varied between 0.46 and 3.85 m/km on the south shore, and 2.2 and 3.75 m/km on the north shore. Drainage density values varied between 0.33 and 0.53 km/km2 on the south shore, and 0.37 and 0.49 km/km2 on the north shore. Other morphological variables, such as drainage length, were strongly correlated with the surface area and/or mean slope of the watersheds.
With respect to land use, the proportion of agricultural area was greater than 10% on the south shore south of 47° N, and below 10% on the north shore (Table 1). However, the proportion of wetland surface area was less than 5% on the south shore, but greater than 8% on the north shore (except for the watershed of the L’Assomption River, which was the most agricultural). Overall, the south shore had a smaller forest area than the north shore. There was a significant difference between the two shores, south of 47° N. Finally, in all of the watersheds analyzed, the amount of urbanized surface area was less than 1%, with the exception of the watersheds of the Châteauguay River (43%), on the south shore, and L’Assomption River (5%), on the north shore. There was little water extraction from aquifers for irrigation or industry, even in the highly agricultural watersheds on the south shore.
Flow data were taken from the website of the Ministry of Environment and the Fight against Climate Change’s Water expertise Centre of Quebec (https://www.cehq.gouv.qc.ca/index_en.asp, accessed on 20 February 2020). Data on physiographic variables were also kindly provided to us by the Centre. They include mean watershed slope and drainage density, as well as percentages of forest, wetland, agricultural and urbanized areas. Climate data was taken from the Environment Canada website (https://climat.meteo.gc.ca/climate_normals/index_f.html, accessed on 18 June 2021). These are the climatic norms, measured at a station in each watershed during the period of 1971–2000. These climatic norms were compared to those of 1941–1970 and 1981–2010, when values were available. In the case of a difference between the three, we calculated their means. The climatic variables analyzed were seasonal precipitation and temperature.
For each of the 17 rivers under study, we created two seasonal hydrological series: spring (April to June) and fall (October to December) daily minimum extreme flows. This involved measuring the lowest flow values of one (1) day per season, each year between 1930 and 2019. We calculated the arithmetic mean (Mean) for each hydrological series and determined the maximum (Max) and minimum (Min) values.

2.2. Statistical Analysis of the Hydrological Series

To analyze the spatial variability of daily minimum extreme flows, we compared the arithmetic means of the hydrological series using parametric and non-parametric tests (ANOVA and Kruskal–Wallis respectively). The arithmetic means, as well as maximum and minimum flow values, were then correlated with the watersheds’ physiographic and climatic factors. To ensure that watershed size would have no influence on the results, we converted these three flow variables into specific flows (L/s/km2).
We used two types of tests to analyze the temporal variability of daily minimum extreme flows. The first was to analyze the long-term trend in the hydrological series using six different statistical tests based on the Mann–Kendall (MK) test. The purpose of some of these tests was to eliminate the effects of short-term persistence (STP) by removing the influence of autocorrelation (MMK-PW and TFPW), or correcting the variance introduced by autocorrelation (MMKY and MMKH). Other tests were designed to eliminate the effects of long-term persistence (LTP). Only one test (MMK-LTP) was used for this purpose. To avoid weighing down the text with formulas, all tests have been summarized in Table 2 and supplemented with references containing detailed descriptions. For more information on their mathematical descriptions and application in hydroclimatology, we direct readers to this section.
In the final stage of the statistical analysis, we conducted three more tests to detect shifts in the mean of the hydrological series: the Pettitt (PT) test, the modified Pettitt test (TPM), and the Lombard test (LT). While the first two tests exclusively detect abrupt shifts, the final test can detect both abrupt and gradual (smooth) shifts. Unlike the original Pettitt (PT) test [49], the modified Pettitt test (TPM) and the Lombard test eliminate the influence of autocorrelation on the detection of shift in mean [50,51,52].

3. Results

3.1. Comparison of Spatial Variability of Daily Minimum Extreme Flows in Spring and Fall

Spring and fall daily minimum extreme flow values are presented in Table 3 and Figure 2. In spring, average flow values were higher on the north shore than on the south shore. The difference between the two shores was most significant south of 47° N, between the rivers of the southwestern hydroclimatic region (>7 L/s/km2) and those of the southeastern hydroclimatic region (<6 L/s/km2). The same applies to minimum values. This difference between the flow values of the two hydroclimatic regions of the two shores was also observed in the fall, even for maximum values.
Based on the analysis of correlation coefficients, calculated between these values and the physiographic and climatic variables, the mean and minimum values of the daily minimum extreme flows were significantly better correlated with main slope (positive correlation) and agricultural areas (negative correlation) during the two seasons (Table 4). There were also positively correlated with wetland areas, in spring only. In fall, this factor was correlated only with the minimum value of flows. Forest areas were positively correlated with the mean value of flows in spring. Maximum values were not significantly correlated with any physiographic factor during the two seasons. Finally, no climatic factor was significantly correlated with daily minimum extreme flows during the two seasons.

3.2. Comparison of Temporal Variability (Long-Term Trend) of Daily Extreme Flows in Spring and Fall

The results of six trend analysis tests are presented in Table 5 (spring) and Table 6 (fall). Except for the MMKY test, the other five tests yielded consistent results in spring. In total, a statistically significant long-term trend in daily extreme minimum flows was only observed in six rivers (35% of the rivers under analysis): five on the south shore (three to the south and two to the north of 47° N) and only one (L’Assomption) on the north shore, in spring. This trend was also observed in fall. Once again, except for the MMKY, a significant long-term trend in daily minimum extreme flows was observed in eight series: seven on the south shore and only one on the north shore. Unlike in spring, the trend was observed in almost all rivers (four out of six) south of 47° N. While long-term trends tended to be negative (most Z values were negative) on both shores in spring, except for south-shore rivers located south of 47° N, they tended to become positive in fall, except for south-shore rivers located north of 47° N (Eastern hydroclimatic region). As such, daily minimum extreme flows tended to increase in the southeastern hydroclimatic region in both spring and fall.
These long-term trends were partially confirmed by the three tests to detect shifts in means, the results of which were nearly consistent. In spring, these three tests detected breaks in the mean of flow of four rivers, three of which were located north of 47° N. However, the dates of these shifts were not synchronous, occurring in different decades (Table 7). In fall, breaks in mean were detected in seven rivers, five of which were located south of 47° N (Table 8). Surprisingly, these breaks in mean were nearly synchronous, and occurred during the 1960s (Figure 3) in these five rivers, which clearly indicates that their consistency is regional. In the two rivers north of this parallel, these shifts occurred later.
In the 1980s and 1990s, none of the three tests detected shifts in the mean of daily minimum flows on the north shore in either season, which suggests that there were also no shifts in the mean of daily minimum extreme flows in the rivers on the north shore in either season.
Rivers in the southeastern hydroclimatic region were further analyzed in the fall to determine if a shift in mean occurred subsequent to the break in the 1960s. The MK tests, and the tests to detect changes in shifts in mean, were repeated, and revealed no other statistically significant changes in either the long-term trend or the means of the hydrological series after the 1960s. Consequently, the moisture of the 2010s did not cause a significant increase in daily minimum extreme flows during the two transitional seasons in southern Quebec. However, Figure 3 shows that, even in the most agricultural south shore watersheds south of 47° N, minimum flows decreased over the course of the decade, despite a substantial increase in flows in the 1960s.

4. Discussion

Analysis of the spatial variability of minimum extreme flows in spring and fall revealed a clear difference between the rivers on the south shore of the St. Lawrence River (south of 47° N) and those on the north shore. The means of these flows tended to be greater on the north shore (>5 L/s/km2) than on the south shore (<5 L/s/km2). The correlation analysis between these flows and physiographic factors revealed that two main factors influence the spatial variability of flows: agricultural areas and mean slope, which, respectively, were negatively and positively correlated. In the spring, these flows are also influenced by wetland areas (positive correlation). The south shore watersheds, located south of 47° N, tend to have relatively large agricultural areas (>20%), but very small wetlands (<4%). Conversely, the watersheds on the north shore have larger wetlands (>8%) but very small agricultural areas (<5%). The impacts of land use on river flows are well documented in the scientific literature, even in Canada [53,54]. With respect to agriculture, previous research in Quebec has shown that agricultural land use reduces minimum flows due to soil sealing, which limits infiltration into aquifers, the main source of streamflow during low-water periods in winter, summer, and fall [55,56,57,58].
This soil sealing also accelerates spring recession in agricultural watersheds [51], which explains the low minimum extreme flows observed in spring and fall in the most agricultural watersheds on the south shore. As for the impacts of wetlands on these minimum extreme flows, there is still no research on the topic in Quebec, and very little internationally. In their summary on the hydrological impacts of wetlands, Bullock and Acreman [52] highlighted the complexity of the relationship between wetlands, aquifers, and minimum flows. This complexity led to divergent results: some studies showed that minimum flows decreased in wetlands, while others showed an increase. In Quebec, it is quite clear that wetlands cause an increase in extreme minimum flows in spring. In this season, minimum flows are mainly associated with recession caused by snowmelt. Thus, wetlands slow this recession, as minimum flows are maintained at higher values than those observed in agricultural watersheds on the south shore (the “sponge effect”), where discharge from snowmelt and/or rainwater is much faster. A recent study showed that spring freshets caused by snowmelt don’t last as long in agricultural watersheds on the south shore than in those on the north shore, as runoff is more significant and faster in the first watersheds (Assani, 2021; unpublished work). In fall, minimum flows are mainly influenced by groundwater. In this season, precipitation mainly falls as rain. But, unlike freshets and rain-induced flood flows, the “sponge effect” of wetlands on the daily minimum extreme flows (July to November) acts tacitly, and results in a slower decrease in flows during low-flows events, probably due to the poor connectivity between these wetlands and aquifers [59]. This may partly explain the significant positive correlation observed between wetlands and the lowest value of the daily minimum extreme flows in fall.
In addition to agricultural areas and wetlands, there was a significant positive correlation between the mean slope and daily minimum extreme flows in both spring and fall. This correlation is explained by the fact that this physiographic factor, in the very agricultural south shore watersheds, is generally less steep than on the north shore, because a large share of the south shore watersheds’ surface area lies in the St. Lawrence Lowlands, a geological formation with a relatively flat topography. By contrast, nearly all the north shore watersheds are located on the Canadian Shield, which has a relatively irregular topography (consisting of hills and valleys). Thus, this correlation between these two variables does not express a direct causal link, but a covariation due to agricultural areas of the watersheds.
Analysis of the temporal variability of daily minimum extreme flows also revealed a marked difference between south-shore and north-shore rivers south of 47° N, particularly in fall. In fact, at that time of year, south-shore rivers located south of this parallel tended to see a significant increase in daily minimum extreme flows, while no significant changes were observed in the north-shore rivers. This increase in flows happened in near-perfect synchronicity for most rivers during the 1960s. This increase can be explained by two factors:
  • The modernization of agricultural practices in southern Quebec since 1950, which led to a significant reduction of cultivated land. Cultivated land was left uncropped or reforested, the amount of which nearly tripled in area between 1950 and 1981 [60]. From a hydrological perspective, this type of land use change would promote the infiltration of water into the soil, increasing fall minimum extreme flows supplied primarily by groundwater. In addition to reduction in the cultivated area, the type of the crop has also changed. As a result, forage crop acreage was reduced by over half to make room for oats, corn, and soybeans, particularly between 1950 and 2011 [60]. This change in crop type could change the conditions for runoff and infiltration of snowmelt and rainwater to varying degrees. The increase in minimum flows resulting from changes in agricultural practices has been observed in many watersheds in agricultural regions of the United States [27].
  • Increased rainfall in summer and fall in the northeastern part of North America. Several studies have noted a significant increase in rainfall, particularly in fall, in the northeastern part of North America north of 44° N [13,17,22,25], and on the south shore of Quebec south of 47° N [61,62]. This increase occurred during the 1960s and early 1970s.
The reduction of cultivated land in favor of fallow land, in addition to the increase in precipitation, may therefore explain the increase in daily minimum extreme flows observed in fall in the southeastern hydroclimatic region (south of 47° N)—the most agricultural region in Quebec. This increase was also observed in rivers in the northeastern part of the United States [29,32].
The predominant fact to emerge from our analysis of temporal variability in spring was the tendency of daily minimum extreme flows to decrease on the south shore north of 47° N (Eastern hydroclimatic region), in contrast with the region just south of this parallel. This decrease may be linked to snowfall in winter [28,36], although this decrease in snow cover was more extensive on the north shore than on the south shore [36]. However, no significant changes in the means of flows were observed on the north shore. It follows that the decrease in amount of snow does not seem to affect the spring daily minimum extreme flows. The impacts of this decrease in the amount of snow would probably be mitigated by wetlands in the region. This hypothesis deserves to be explored.

5. Conclusions

Analysis of the spatial and temporal variability of daily minimum extreme flows during the transitional seasons (spring and fall) between 1930 and 2019 revealed a marked difference between rivers on the north and south shores of the St. Lawrence River south of 47° N. In terms of spatial variability, rivers on the south shore have lower daily minimum extreme flows than those on the north shore. This difference was mainly due to agricultural areas. South shore watersheds have a larger agricultural area than those on the north shore. Agriculture reduces water infiltration in fall and accelerates runoff in spring.
With respect to temporal variability, the same difference was observed between the north and south shores in fall. Daily minimum extreme flows have increased significantly in south-shore rivers since the 1960s. This increase is due to a reduction in cultivated land following the modernization of agricultural practices in Quebec since the 1960, as well as an increase in rainfall in the region since the 1970s. Decreased snowfall in Quebec does not yet appear to significantly affect the long-term trend of daily minimum extreme flows in spring, despite the observation of a downward trend north of 47° N on both the south and north shores. In fall, the overall increase in temperature in summer and early fall did not lead to a decrease in fall daily minimum extreme flows, as the increase in rainfall offsets the increase in evapotranspiration resulting from this summer–fall warming.
According to climate models, warmer temperatures in Quebec will result in a significant decrease in the amount of snow in spring, which will in turn reduce runoff and infiltration. This will result in a significant decrease in the magnitude of spring freshets and low-water levels. In fall, these climate models predict a significant increase in temperature that will be accompanied by more intense rainfall. This scenario will result in a continuous increase in low-water flows. There will therefore be a continuous decrease in low-water flows in spring but a continuous increase in low-water flows in fall. It is this trend that will be followed in the future.

Author Contributions

Conceptualization, A.A.A.; methodology, A.A.A. and A.Z.; formal analysis, A.A.A. and A.Z.; investigation, A.A.A., A.Z., A.R. and C.K.; writing—original draft preparation, A.A.A., A.Z., A.R. and C.K.; writing—review and editing, A.A.A., A.Z., A.R. and C.K.; funding acquisition, A.A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Sciences and Engineering Research Council of Canada (NSERC), grant number “261274/2019”.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Danco, J.F.; DEAngilils, A.; Raney, B.K.; Broccoli, A.J. Effects of a warming climate on daily snowfall events in the Northern Hemisphere. J. Clim. 2016, 29, 6295–6318. [Google Scholar] [CrossRef]
  2. Kapnick, S.B.; Delworth, T.L. Controls of global snow under a changed climate. J. Clim. 2013, 26, 5537–5562. [Google Scholar] [CrossRef]
  3. Krasting, J.P.; Broccoli, A.J.; Dixon, K.W.; Lanzante, J.R. Future changes in Northern Hemisphere snowfall. J. Clim. 2013, 26, 7813–7828. [Google Scholar] [CrossRef] [Green Version]
  4. O’Gorman, P. Contrasting responses of mean and extreme snowfall to climate change. Nature 2014, 512, 416–418. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Janoski, T.P.; Broccoli, A.J.; Kapnick, S.B.; Johnson, N. Effects of climate change on wind-driven heavy-snowfall events over Eastern North America. J. Clim. 2018, 31, 9037–9054. [Google Scholar] [CrossRef]
  6. Ingram, W. Extreme precipitation: Increases all round. Nat. Clim. Chang. 2016, 6, 443–444. [Google Scholar] [CrossRef]
  7. Tabari, H. Climate change impact on flood and extreme precipitation increase with water availibility. Nat. Sci. Rep. 2020, 10, 13768. [Google Scholar] [CrossRef] [PubMed]
  8. Marx, A.; Kumar, R.; Thober, S.; Zink, M.; Wanders, N.; Wood, E.F.; Pan, M.; Sheffield, J.; Samaniogo, L. Climate change alters low flows in Europe under a 1.5, 2, and 3 degree global warming. Hydrol. Earth. Syst. Sci. Discuss. 2017, 1–24. [Google Scholar] [CrossRef] [Green Version]
  9. Assefa, K.; Moges, M.A. Low Flow Trends and Frequency Analysis in the Blue Nile Basin, Ethiopia. J. Water Res. Protec. 2018, 10, 182–203. [Google Scholar] [CrossRef] [Green Version]
  10. Blum, A.G.; Archfield, S.A.; Hirsch, R.M.; Vogel, R.M.; Kiang, J.E.; Dudley, W. Updating estimates of low-streamflow statistics to account for possible trends. Hydrol. Sci. J. 2019, 64, 1404–1414. [Google Scholar] [CrossRef]
  11. Ehsanzadeh, E.; Adamowski, K. Trends in timing of low stream flows in Canada: Impact of autocorrelation and long-term persistence. Hydrol. Process. 2010, 24, 970–980. [Google Scholar] [CrossRef]
  12. Déry, S.J.; Hernández-Henríquez, M.A.; Owens, P.N.; Parkes, M.W.; Petticrew, E.L. A century of hydrological variability and trends in the Fraser Rive basin. Environ. Res. Lett. 2012, 7, 024019. [Google Scholar] [CrossRef]
  13. Dudley, R.W.; Hirsh, R.M.; Archfiled, S.A.; Blum, A.G.; Renard, B. Low trends at human-impacted and reference basins in the United States. J. Hydrol. 2020, 580, 124254. [Google Scholar] [CrossRef]
  14. Gain, A.K.; Immerzeel, W.W.; Sperna-Weiland, F.C.; Bierkens, M.F.P. Impact of climate change on the stream flow of lower Brahmaputra: Trends in high and low flows based on discharge-weighted ensemble modelling. Hydrol. Earth Syst. Sci. 2011, 8, 365–390. [Google Scholar] [CrossRef] [Green Version]
  15. Hammond, J.C.; Fleming, B.J. Evaluating low flow patterns, drivers and trends in the Delaware river basin. J. Hydrol. 2021, 598, 126246. [Google Scholar] [CrossRef]
  16. Hannaford, J. Climate-driven changes in UK river flows: A review of the evidence. Progr. Phys. Geogr. 2015, 39, 29–48. [Google Scholar] [CrossRef] [Green Version]
  17. Hodgkins, G.A.; Dudley, R.W. Historical summer base flow and stromflow trends for New England rivers. Water Resourc. Res. 2011, 47, W07528. [Google Scholar] [CrossRef]
  18. Khaliq, M.N.; Ouarda, T.B.M.J.; Gachon, P. Identification of temporal trends in annual and seasonal low flows occuring in Canadian rivers: The effect of short- and long-term persistence. J. Hydrol. 2009, 369, 183–197. [Google Scholar] [CrossRef]
  19. Khaliq, M.N.; Ouarda, T.B.M.J.; Gachon, P.; Sushama, L. Temporal evolution of low-flow regimes in Canadian rivers. Water Resourc. Res. 2008, 44, W08436. [Google Scholar] [CrossRef] [Green Version]
  20. Kay, A.L.; Bell, V.A.; Guillod, B.P.; Jones, R.G. Rudd, A.C. National-scale analysis of low flow frequency: Historical trends and potential future changes. Clim. Chang. 2018, 147, 585–599. [Google Scholar] [CrossRef] [Green Version]
  21. Kormos, P.R.; Luce, C.H.; Wenger, S.J.; Berghuijs, W.R. Trends and sensitivities of low streamflow extremes to discharge timing and magnitude in Pacific Northwest mountain streams. Water Resour. Res. 2016, 52, 4990–5007. [Google Scholar] [CrossRef] [Green Version]
  22. Sadri, S.; Kam, J.; Sheffield, J. Nonstationarity of low flows and their timing in the eastern United States. Hydrol. Earth Syst. Sci. 2016, 20, 633–649. [Google Scholar] [CrossRef] [Green Version]
  23. Sahoo, B.B.; Jha, R. Assessment of low flow trends and change point detection in Mahanadi River basin, India. Sustain. Water Resour. Manag. 2020, 6, 81. [Google Scholar] [CrossRef]
  24. Sawaske, S.P.; Freyberg, D.L. An analysis of trends in baseflow recession and low-flows in rain-dominated coastal streams of the pacific coast. J. Hydrol. 2014, 519, 599–610. [Google Scholar] [CrossRef]
  25. Small, D.; Islam, S.; Vogel, R.M. Trends in precipitation and streamflow in the eastern U.S.: Paradox or perception? Geophys. Res. Lett. 2006, 13, L03403. [Google Scholar] [CrossRef]
  26. Vlach, V.; Ledvinka, O.; Matouskova, M. Changing low flow and streamflow drought seasonality in Central European headwaters. Water 2020, 12, 3575. [Google Scholar] [CrossRef]
  27. Zhang, Y.-K.; Schilling, K.E. Increasing streamflow and base flow in Mississippi River since the 1940s: Effect of land use change. J. Hydrol. 2006, 324, 412–422. [Google Scholar] [CrossRef]
  28. Brown, R.D. Analysis of snow cover variability and change in Québec, 1948–2005. Hydrol. Process. 2010, 24, 1929–1954. [Google Scholar] [CrossRef]
  29. Fortin, G.; Hétu, B. Estimating winter trends in climatic variables in the Chic-Chocs Mountains, Canada (1970–2009). Int. J. Climatol. 2013, 34, 3078–3088. [Google Scholar] [CrossRef]
  30. Yagouti, A.; Boulet, G.; Vincent, L.; Vescovi, L.; Mekis, E. Observed changes in daily temperature and precipitation indices for Southern Québec, 1960–2005. Atmos. Ocean. 2008, 46, 243–256. [Google Scholar] [CrossRef]
  31. Assani, A.A. Comparison of the temporal variability of maximum daily temperatures for summer months in relation to El Nino events in Southern Québec. In Exteme Weather; Sallis, P.J., Ed.; IntechOpen: London, UK, 2008. [Google Scholar] [CrossRef]
  32. Assani, A.A.; Maloney-Dumont, V.; Pothier-Champagne, A.; Kinnard, C.; Quéessy, J.-F. Comparison of the temporal variability of summer temperature and rainfall as it relates to climate indices in southern Quebec (Canada). Theor. Appl. Climatol. 2019, 137, 2425–2435. [Google Scholar] [CrossRef]
  33. Assani, A.A.; Landry, R.; Quessy, J.-F.; Clément, F. Temporal Variability of Rain-Induced Floods in Southern Quebec, Climate Change. In Geophysical Foundations and Ecological Effects; Blanco, J., Kheradmand, H., Eds.; IntechOpen: London, UK, 2011. [Google Scholar] [CrossRef] [Green Version]
  34. Assani, A.A.; Landry, R.; Laurencelle, M. Comparison of interannual variability modes and trends of seasonal precipitation and streamflow in Southern Quebec (Canada). River Res. Appl. 2012, 28, 1740–1752. [Google Scholar] [CrossRef]
  35. Beauchamp, M.; Assani, A.A.; Landry, R.; Massicotte, P. Temporal variability of the magnitude and timing of winter maximum daily flows in southern Quebec (Canada). J. Hydrol. 2015, 520, 410–417. [Google Scholar] [CrossRef]
  36. Guerfi, N.; Assani, A.A.; Mesfioui, M.; Kinnard, C. Comparison of the temporal variability of winter daily extreme temperatures and precipitations in southern Quebec (Canada) using the Lombard and copula methods. Int. J. Climatol. 2015, 35, 4237–4246. [Google Scholar] [CrossRef]
  37. Assani, A.A.; Charron, S.; Matteau, M.; Mesfioui, M.; Quessy, J.F. Temporal variability modes of floods for catchments in the St. Lawrence Watershed (Quebec, Canada). J. Hydrol. 2010, 385, 292–299. [Google Scholar] [CrossRef]
  38. Mazouz, R.; Assani, A.A.; Quessy, J.-F.; Légaré, G. Comparison of the interannual variability of spring heavy floods characteristics of tributaries of St. Lawrence River in Quebec (Canada). Adv. Water Res. 2012, 35, 110–120. [Google Scholar] [CrossRef]
  39. Mazouz, R.; Assani, A.A.; Rodriguez, M. Application of redundancy analysis to hydroclimatology: A case study of spring heavy floods in southern Québec (Canada). J. Hydrol. 2013, 496, 187–194. [Google Scholar] [CrossRef]
  40. Smakhtin, V.U. Low flow hydrology: A review. J. Hydrol. 2001, 240, 147–186. [Google Scholar] [CrossRef]
  41. Rolls, R.J.; Leigh, C.; Sheldon, F. Mechanistic effects of low-flow hydrology on riverine ecosystems: Ecological principles and consequences of alteration. Freshw. Sci. 2012, 31, 1163–1186. [Google Scholar] [CrossRef] [Green Version]
  42. Belzile, L.; Bérubé, P.; Hoang, V.D.; Leclerc, M. Méthode Écohydrologique de Détermination des débits Réservés pour la Protection des Habitats du Poisson dans les Rivières du Québec. Rapport Présenté par l’INRS-Eau et le Groupe-Conseil Génivar inc. au Ministère de l’Environnement et de la Faune et à Pêches et Océans Canada. 1997; 83p. Available online: http://www.obvcapitale.org/plans-directeurs-de-leau-2/2e-generation/diagnostic/section-3-problematiques-associees-a-la-quantite-et-a-la-securite/3-4-debits-reserves (accessed on 20 December 2020).
  43. Sneyers, R. On the Statistical Analysis of Series of Observations; Technical Notes N°143; World Meteorological Organization: Geneva, Switzerland, 1990; 192p. [Google Scholar]
  44. Von Storch, H. Misuses of statistical analysis in climate research. In Analysis of Climate Variability; von Storch, H., Navarra, A., Eds.; Springer: Dordrecht, The Netherlands, 1995; pp. 11–26. [Google Scholar]
  45. Yue, S.; Pilon, P.; Cavadias, G. Power of the Mann-Kendall and Spearman’s rho tests for detecting monotonic trends in hydrological series. J. Hydrol. 2002, 259, 254–271. [Google Scholar] [CrossRef]
  46. Yue, S.; Wang, C.Y. The Mann-Kendall test modified by effective sample size to detect trend in serially correlated hydrological series. Water Resourc. Manag. 2004, 18, 201–218. [Google Scholar] [CrossRef]
  47. Hamed, K.H.; Rao, A.R. A modified Mann-Kendall trend test for autococorrelated data. J. Hydrol. 1998, 204, 182–196. [Google Scholar] [CrossRef]
  48. Hamed, K.H. Enhancing the effectiveness of prewhitening in trend analysis of hydrologic data. J. Hydrol. 2009, 368, 143–155. [Google Scholar] [CrossRef]
  49. Pettitt, A.N. A non-parametric approach to the change-point problem. J. R. Statist. Soc. Ser. C 1979, 28, 126–135. [Google Scholar] [CrossRef]
  50. Serinaldi, F.; Kilsby, C. The importance of prewhiting in change point analysis under persitence. Stoch. Environ. Res. Risk Assess. 2016, 30, 763–777. [Google Scholar] [CrossRef] [Green Version]
  51. Lombard, F. Rank tests for changepoint problems. Biometrika 1987, 74, 615–624. [Google Scholar] [CrossRef]
  52. Quessy, J.F.; Favre, A.C.; Saïd, M.; Champagne, M. Statistical inference in Lombard’S smooths-change model. Environmetrics 2011, 22, 882–893. [Google Scholar] [CrossRef]
  53. Mao, D.; Cherkauer, K.A. Impacts of land-use change on hydrologic responses in the Great Lakes region. J. Hydrol. 2009, 374, 71–82. [Google Scholar] [CrossRef]
  54. Tan, X.; Gan, T.Y. Contribution of human and climate change impacts to changes in streamflow of Canada. Sci. Rep. 2015, 5, 17767. [Google Scholar] [CrossRef]
  55. Assani, A.A.; Landry, R.; Kinnard, C.; Azouaoui, O.; Demers, C.; Lacasse, K. Comparison of the spatiotemporal variability of temperature, precipitation, and maximum daily spring flows in two watersheds in Quebec characterized by different land use. Adv. Meteorol. 2016, 2016, 3746460. [Google Scholar] [CrossRef] [Green Version]
  56. Muma, M.; Assani, A.A.; Landry, R.; Quessy, J.-F.; Mesfioui, M. Effects of the change from forest to agriculture land use on the spatial variability of summer extreme daily flow charactersitics in southern Quebec (Canada). J. Hydrol. 2011, 407, 153–163. [Google Scholar] [CrossRef]
  57. Quilbé, R.; Rousseau, A.; Moquet, J.-S.; Savary, S.; Ricard, S.; Garbouj, M.S. Hydrological response of a watershed to historical land use evolution and future land use scenarios under climate change conditions. Hydrol. Earth Syst. Sci. 2008, 12, 101–110. [Google Scholar] [CrossRef] [Green Version]
  58. Landry, R.; Assani, A.A.; Biron, S.; Quessy, J.-F. The management modes of seasonal floods and their impact on the relatioship between climate and streamflow downstream from dams in Quebec (Canada). River Res. Appl. 2014, 30, 287–298. [Google Scholar] [CrossRef]
  59. Bullock, A.; Acreman, M.C. The role of wetlands in the hydrological cycle. Hydrol. Earth Syst. Sci. 2003, 7, 358–389. [Google Scholar] [CrossRef] [Green Version]
  60. Ruiz, J. Modernization agriculture and agricultural land cover in Quebec (1951–2011). Cahiers Géographie Québec 2019, 63. in press (In French) [Google Scholar]
  61. Perrault, L.; Haché, M.; Slivitzky, M.; Bobée, B. Detection of changes in precipitation and runoff over eastern Canada and U.S. using a Bayesian approach. Stoch. Environ. Res. Risk Assess. 1999, 13, 201–216. [Google Scholar] [CrossRef]
  62. Assani, A.A.; Lajoie, F.; Vadnais, M.-E.; Beauchamp, G. Influence of the Arctic oscillation on the interannual variability of precipitation in the Saint-François river watershed (Quebec, Canada) as determined by canonical correlation analysis. J. Water Sci. (Rev. Sci. Eau) 2008, 21, 21–33. (In French) [Google Scholar]
Figure 1. Location of river stations. SW = Southwestern hydroclimatic Region; SE = Southeastern Hydroclimatic Region; E = Eastern Hydroclimatic Region; the circles represent major towns.
Figure 1. Location of river stations. SW = Southwestern hydroclimatic Region; SE = Southeastern Hydroclimatic Region; E = Eastern Hydroclimatic Region; the circles represent major towns.
Water 13 03487 g001
Figure 2. Comparison of average daily minimum extreme flow values in spring (blue bars) and fall (red bars) between 1930 and 2019.
Figure 2. Comparison of average daily minimum extreme flow values in spring (blue bars) and fall (red bars) between 1930 and 2019.
Water 13 03487 g002
Figure 3. Example of the interannual variability in the daily minimum extreme flows in fall of three rivers in the southeast hydroclimatic region on the south shore. The dashed vertical green line indicates the year of shift in mean of the hydrological series. Beaurivage River (blue curve); Etchemin River (red curve); Nicolet SW River (black curve). The shift in mean of flows occurred almost the same year for all three rivers.
Figure 3. Example of the interannual variability in the daily minimum extreme flows in fall of three rivers in the southeast hydroclimatic region on the south shore. The dashed vertical green line indicates the year of shift in mean of the hydrological series. Beaurivage River (blue curve); Etchemin River (red curve); Nicolet SW River (black curve). The shift in mean of flows occurred almost the same year for all three rivers.
Water 13 03487 g003
Table 1. Rivers analyzed.
Table 1. Rivers analyzed.
RiversCodeDrainage Area (Km2)Latitude
(N)
Longitude
(W)
Agricultural Area (%)
Southeastern Hydroclimatic Region (South Shore)
ChateaugaySE1249245°19′49″73°45′44″41.9
EatonSE264645°28′05″71°39′18″10.7
Nicolet SWSE356245°47′30″71°58′05″26.3
EtcheminSE4115246°39′25″71°39′18″25.5
BeaurivageSE570846°39′25″71°17′20″34.7
Du SudSE682146°49′22″70°45′22″10.4
Eastern Hydroclimatic Region (South Shore)
OuelleE179647°22′52″67°57′14″2.8
Du LoupE2104247°36′43″69°38′41″10.2
Trois-PistolesE393048°05′21″69°11′43″16.5
RimouskiE4161548°24′46″68°33′18″8.8
MataneE5166548°46′25″67°32′25″8.9
BlancheE622348°47′20″67°41′51″30.2
Southernwest Hydroclimatic Region (North Shore)
Petite NationSW1133145°47′27″75°05′22″0.7
Du NordSW2116345°31′08″74°20′11″0.4
L’AssomptionSW3128646°02′45″73°26′19″8.6
MatawinSW4138746°40′50″73°55′00″0
VermillonSW5266247°39′20″72°57′44″0
Table 2. Summary of Mann–Kendall tests applied to the hydrological series under analysis.
Table 2. Summary of Mann–Kendall tests applied to the hydrological series under analysis.
Name of the TestAcronymPurpose of the TestReferences
Elimination of STP (Short Term Persistence) Effects
Original Mann–KendallMKDetect the long-term trend (does not take STP or LTP into account).[43]
Prewhitenning methodMMK-PWEliminate the influence of autocorrelation by prewhitenning data.[44]
Trend Free Prewhitening MethodTFPWEliminate autocorrelation by prewhitening data.[45]
Modified Mann–Kendall Test1MMKYEliminate autocorrelation by variance correction.[46]
Modified Mann–Kendall Test2MMKHEliminate autocorrelation by variance correction.[47]
Elimination of LTP (Long Term persistence) Effects
Long Term Persistence MKLTPEliminate LTP effects.[48]
Table 3. Comparison of seasonal minimum extreme daily flows (L/s/km2) in spring and fall between 1930 and 2019.
Table 3. Comparison of seasonal minimum extreme daily flows (L/s/km2) in spring and fall between 1930 and 2019.
RiversSpringFall
MaxMeanMinMaxMeanMin
Southeastern Hydroclimatic Region
Chateaugay9.23.3 (1.61)0.569.32.6 (1.74)0.28
Eaton7.03.1 (1.30)0.888.13.2 (1.62)0.78
Nicolet SW10.82.8 (1.95)0.258.33.4 (1.98)0.53
Etchemin14.13.6 (2.67)0.8818.14.5 (2.30)1.00
Beaurivage47.33.6 (4.98)0.569.53.3 (1.94)0.16
Du Sud13.55.5 (2.52)0.3013.04.6 (2.08)0.12
Eastern Hydroclimatic Region
Ouelle16.13.9 (2.38)0.547.32.9 (1.74)0.35
Du Loup13.15.8 (2.42)0.7111.74.3 (2.43)0.27
Trois-Pistoles11.44.1 (2.14)1.43.21.3 (0.59)0.37
Rimouski25.46.3 (3.95)1.715.94.7 (2.76)0.70
Matane21.48.0 (4.22)1.7315.94.7 (2.77)0.70
Blanche21.15.2 (3.43)2.018.34.1 (3.75)0.38
Southernwest Hydroclimatic Region
Petite Nation24.39.8 (4.01)3.620.94.8 (3.26)1.2
Du Nord16.58.1 (2.79)3.915.27.0 (2.56)2.7
L’Assomption14.37.9 (2.83)2.512.95.4 (2.80)0.72
Matawin27.28.9 (4.33)2.714.26.6 (2.45)2.0
Vermillon15.67.6 (3.09)2.214.46.5 (2.64)1.7
(1.61) = standard deviation.
Table 4. Correlation coefficients calculated between physio-climatic variables and magnitude (L/s/km2) of daily minimum extreme flows in spring and fall between 1930 and 2019.
Table 4. Correlation coefficients calculated between physio-climatic variables and magnitude (L/s/km2) of daily minimum extreme flows in spring and fall between 1930 and 2019.
VariablesSpringFall
MaxMeanMinMaxMeanMin
Physiographic Variables
Drainage density (km/km2)0.472−0.119−0.147−0.014−0.127−0.216
Mean slope (m/km)−0.0300.671 **0.550 **0.4690.669 **0.549 **
Forests surface area (%)0.0020.484 **0.3880.1500.4630.381
Agricultural surface area (%)0.081−0.687 **−0.648 **−0.231−0.580 **−0.572 **
Wetlands surface area (%)0.2660.666 **0.789 **0.4470.4260.528 **
Climatic Variables
Seasonal (spring/fall) total rainfall (mm)−0.0820.161−0.0820.094−0.146−0.088
Seasonal (spring/fall) total snowfall (cm)0.253−0.271−0.407−0.052−0.399−0.370
Seasonal (spring/fall) total precipitation (mm)0.0150.0710.1380.007−0.399−0.351
Winter-spring/ summer-fall total rainfall (mm)−0.0980.0990.2540.109−0.087−0.160
Winter-spring/summer-fall snowfall (cm)0.058−0.139−0.282−0.052−0.399−0.369
Winter-spring/summer-fall total precipitation (mm)0.001−0.018−0.0390.052−0.284−0.321
Seasonal (Winter/Spring) daily mean maximum temperature (°C)−0.4000.1220.274−0.238−0.118−0.003
Winter-Spring/Summer-Fall daily mean maximum temeperature (°C)−0.4230.0560.252−0.0840.0880.121
** = significant values over the 5% threshold.
Table 5. Results of the various Mann–Kendall tests applied to daily minimum flow series in spring between 1930 and 2019.
Table 5. Results of the various Mann–Kendall tests applied to daily minimum flow series in spring between 1930 and 2019.
RiversMKMMK-PWTFPWMMKYMMKHLTP
Zp-ValueZp-ValueZp-ValueZp-ValueZp-ValueZp-Value
Southeastern Hydroclimatic Region
Chateaugay1.1150.2651.3930.1641.2790.2012.621 **0.0091.1020.2711.1900.234
Eaton−0.3170.751−0.5940.553−0.5650.572−0.9660.334−0.3510.726−0.4410.667
Nicolet SW2.2660.0231.4880.1161.3860.6662.884 **0.0041.708 **0.0881.700 *0.089
Etchemin1.4600.1441.4880.1371.6690.0952.791 **0.0051.3640.1721.1100.267
Beaurivage5.3090.0005.1420.0005.5890.0009.559 **0.0004.723 **0.0004.737 **0.000
Du Sud0.2200.8260.6200.5340.6410.5210.5020.6150.1790.8580.1830.855
Eastern Hydroclimatic Region
Ouelle−4.392 **0.000−3.81 **0.000−4.4 **0.000−11.9 **0.000−4.15 **0.000−4.32 **0.000
Du Loup0.5400.5891.0240.3061.1520.2490.9920.3210.5230.6010.2840.777
Trois-Pistoles−2.18 **0.029−1.92 *0.055−2.0 **0.048−3.74 **0.000−2.31 **0.021−2.10 **0.037
Rimouski−1.0560.291−1.1230.261−1.1940.232−1.3190.187−1.0500.294−1.1300.262
Matane−2.27 **0.023−2.69 **0.007−2.7 **0.006−3.84 **0.000−2.46 **0.014−2.92 **0.004
Blanche−0.3840.701−0.0490.961−0.1180.906−0.5960.551−0.4500.653−0.2190.830
Southwestern Hydroclimatic Region
Petite Nation1.0630.2881.4140.1571.2650.2062.4010.0161.0630.2881.2570.209
Du Nord−0.3930.694−0.4580.647−0.4070.684−0.8760.381−0.7610.447−0.6120.548
L’Assomption−1.5060.132−1.4490.147−1.2860.198−5.57 **0.000−1.726 *0.084−2.97 **0.003
Matawin−1.5060.132−0.1570.876−0.1570.876−0.2310.817−0.1210.904−0.0910.931
Vermillon0.8210.4171.2580.2081.0450.2963.5130.0000.9200.3581.0620.288
** = significant values over the 5% threshold are shown in red/blue bold; * = significant values over the 10% threshold are shown in red/blue bold. Blue = negative trends; red = positive trends.
Table 6. Results of the various Mann–Kendall tests applied to daily minimum flow series in fall between 1930 and 2019.
Table 6. Results of the various Mann–Kendall tests applied to daily minimum flow series in fall between 1930 and 2019.
RiversMKMMK-PWTFPWMMKYMMKHLTP
Zp-ValueZp-ValueZp-ValueZp-ValueZp-ValueZp-Value
Southeastern Hydroclimatic Region
Chateaugay1.3620.1730.1733.3993.694 **0.0008.468 **0.0003.831 **0.0002.947 **0.003
Eaton1.1000.2170.9320.3331.0670.2861.968 **0.0491.1000.2710.7220.470
Nicolet SW3.275 **0.0012.671 **0.0083.254 **0.0015.888 **0.0003.744**0.0002.085 **0.037
Etchemin2.754 **0.0062.112 **0.0352.676 **0.0074.545 **0.0001.992**0.0461.2250.221
Beaurivage3.820 **0.0002.662 **0.0083.600 **0.0006.502 **0.0003.574**0.0001.948 *0.051
Du Sud1.99 **0.0441.0630.2881.3950.1632.232 **0.0260.9980.6680.6680.504
Eastern Hydroclimatic Region
Ouelle−0.4720.637−0.3080.758−0.3590.719−0.8070.419−0.4540.650−0.2440.810
Du Loup−0.9910.322−0.8360.403−0.7990.424−2.036 ** 0.042 −0.9910.318−0.7260.471
Trois-Pistoles−2.141 ** 0.032 −1.921 * 0.055 −1.957 * 0.050 −4.26 ** 0.000 −2.14 ** 0.032 −1.6180.107
Rimouski−0.4930.622−0.2560.798−0.3140.754−0.8990.369−0.5050.614−0.2930.773
Matane−2.265 ** 0.024 −1.770 * 0.077 −2.05 ** 0.040 −6.46 ** 0.000 −2.27 ** 0.024 −1.687 * 0.093
Blanche0.7520.4520.5370.5910.5830.5601.5040.1330.7520.4520.4980.620
Southwestern Hydroclimatic Region
Petite Nation1.3620.1731.5030.1331.5320.1252.995 **0.0031.0460.1631.5310.126
Du Nord0.8720.3830.9680.3330.7920.4282.881 **0.0041.1560.2481.860 *0.063
L’Assomption0.9110.3620.9180.3580.7820.4342.455 **0.0140.8690.3851.0850.278
Matawin0.3580.7200.2470.8050.1940.8471.0220.3070.4440.6570.4960.620
Vermillon0.0920.9270.2560.7980.2490.8040.2000.8410.0790.9300.0940z.925
** = significant values over the 5% threshold are shown in red/blue bold; * = significant values over the 10% threshold are shown in red/blue bold. Blue = negative trends; red = positive trends.
Table 7. Results of the various Pettitt and Lombard tests applied to daily minimum flow series in spring between 1930 and 2019.
Table 7. Results of the various Pettitt and Lombard tests applied to daily minimum flow series in spring between 1930 and 2019.
Pettitt TestModified Pettitt TestLombard Test
Kp-ValueTKp-ValueTSnT1–T2
Southeastern Hydroclimatic Region
Chateaugay5450.178-5800.129-0.017-
Eaton2821.001-3170.834-0.003-
Nicolet SW4920.247-5460.177-0.021-
Etchemin6470.066-5910.117-0.029-
Beaurivage1298 **0.00019681219 **0.00019680.258 **1967–68
Du Sud3760.633-3760.633-0.002-
Eastern Hydroclimatic Region
Ouelle984 ** 0.001 1986 889 ** 0.003 1986 0.176 ** 1930–99
Du Loup5300.203-5360.193-0.006-
Trois-Pistoles5930.114-6060.101-0.0489 ** 1951–52
Rimouski3630.684-3410.776-0.010-
Matane765 ** 0.017 1963 856 ** 0.005 1963 0.061 ** 1961–62
Blanche3890.488-2970.879-0.003-
Southwestern Hydroclimatic Region
Petite Nation5250.212-6080.099-0.014-
Du Nord2691.088-2381.242-0.005-
L’Assomption3780.625-4150.492-0.018-
Matawin4820.247-5020.207-0.002-
Vermillon3870.591-4970.268-0.008-
** = significant values at the 5% threshold are shown in red bold; T = year of the shift mean; T1–T2 = year before (T1) and after (T2) of shift mean; K = calculated values of Pettitt test; Sn = calculated values of Lombard test. Blue = drop in the means after the shift; Red bold= Rise in the means after the shift.
Table 8. Results of the two Pettitt and Lombard tests applied to the daily minimum flow series in fall between 1930 and 2019.
Table 8. Results of the two Pettitt and Lombard tests applied to the daily minimum flow series in fall between 1930 and 2019.
Pettitt TestModified Pettitt TestLombard Test
Kp-ValueTKp-ValueTSnT1–T2
Southeastern Hydroclimatic Region
Chateaugay1023 **0.0001964930 **0.00119620.156 **1961–62
Eaton5720.105-5100.192-0.010-
Nicolet SW1137 **0.0001965948 **0.00119640.107 **1963–64
Etchemin966 **0.0011965784 **0.00919650.075 **1964–65
Beaurivage1114 **0.0001964892 **0.00219640.121 **1963–64
Du Sud814 **0.0441970752 **0.01119680.060 **1966–68
Eastern Hydroclimatic Region
Ouelle4310.397-3730.596-0.003-
Du Loup4560.327-4300.400-0.013-
Trois-Pistoles656 * 0.080 1988 660 * 0.051 1988 0.049 ** 1992–93
Rimouski5020.240-4880.270-0.004-
Matane619 * 0.080 1976 556 * 0.148 - 0.055 ** 1982–83
Blanche3080.800-3020.829-0.004-
Southwestern Hydroclimatic Region
Petite Nation4250.415-4690.295-0.010-
Du Nord2980.923-3850.550-0.004-
L’Assomption4420.386-5260.195-0.009-
Matawin3140.798-4090.420-0.002-
Vermillon3760.632-4320.416-0.002-
** (red/blue bold) = significant values over the 5% threshold; * (red bold/blue bold) = significant value over the 10% threshold. T = year of the shift in mean; T1–T2 = year before (T1) and after (T2) the shift in mean; K = calculated values of Pettitt test; Sn = calculated values of Lombard test. Blue = drop in the means after the shift; red = rise in the means after the shift.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Assani, A.A.; Zeroual, A.; Roy, A.; Kinnard, C. Impacts of Agricultural Areas on Spatio-Temporal Variability of Daily Minimum Extreme Flows during the Transitional Seasons (Spring and Fall) in Southern Quebec. Water 2021, 13, 3487. https://doi.org/10.3390/w13243487

AMA Style

Assani AA, Zeroual A, Roy A, Kinnard C. Impacts of Agricultural Areas on Spatio-Temporal Variability of Daily Minimum Extreme Flows during the Transitional Seasons (Spring and Fall) in Southern Quebec. Water. 2021; 13(24):3487. https://doi.org/10.3390/w13243487

Chicago/Turabian Style

Assani, Ali Arkamose, Ayoub Zeroual, Alexandre Roy, and Christophe Kinnard. 2021. "Impacts of Agricultural Areas on Spatio-Temporal Variability of Daily Minimum Extreme Flows during the Transitional Seasons (Spring and Fall) in Southern Quebec" Water 13, no. 24: 3487. https://doi.org/10.3390/w13243487

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