Investigation of Long-Term Climate and Streamflow Patterns in Ontario

To develop mitigation and adaptation strategies for undesired consequences of climate change, it is important to understand the changing hydrological and climatological trends in the past few decades. Although the changing climate is a cause of concern for the entire planet, its effects can vary significantly on a regional scale. Canada has experienced a rapid rise in the annual mean surface air temperature in the past decades. The current study aims to investigate trends in monthly mean precipitation, rainfall, snowfall, maximum and minimum temperature, as well as baseflow, surface runoff, and total streamflow values for the province of Ontario, Canada. To the best of the author’s knowledge, a similar study involving rural and urban watersheds, that quantifies the impact of changing climate on temperature and other hydrological processes over a period ranging from 1968 to 2017, has not yet been conducted for Ontario. Man-Kendall trend test was used to analyze trends in the above mentioned climatic and hydrometric parameters for rural and urban watersheds situated in the northern and southern parts of Ontario. The results of this study indicate that the mean monthly minimum temperatures for rural watersheds situated in southern Ontario have increased significantly for the winter and summer months, which may have caused an increase in snowmelt and consequently the streamflow for the winter months in the region. Unlike the watersheds in southern Ontario, the northern watersheds witnessed relatively fewer instances of significant changes in mean monthly temperatures, and in some cases, declining rates have been noted. Similarly, only a few watersheds in the north saw a substantial drop in baseflow over the summer months. For nearly all the months, the average monthly minimum and maximum temperatures were found to increase for urban watersheds. The streamflow, baseflow, and surface runoff increased, likely due to rapid urbanization, resulting in a lower infiltration rate. These results will contriHow to cite this paper: Azarkhish, A., Rudra, R., Daggupati, P., Dhiman, J., Dickinson, T., & Goel, P. (2021). Investigation of Long-Term Climate and Streamflow Patterns in Ontario. American Journal of Climate Change, 10, 467-489. https://doi.org/10.4236/ajcc.2021.104024 Received: June 28, 2021 Accepted: December 14, 2021 Published: December 17, 2021 Copyright © 2021 by author(s) and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access


Introduction
Climate change has a wide range of potential social, economic, and environmental impacts. Global concerns related to rising sea levels (Church et al., 2013), increase in the frequency and intensity of heat waves (Peterson et al., 2013), and frequent flooding and drought events (Dai, 2013), makes it important to study the change in trends of climatic and hydrological data. However, the impacts of climate change along with land use changes on a regional scale may vary considerably based on geographic locations (Hamilton & Keim, 2009). For instance, in a study conducted by Shahid et al. (2018) it was found that climate change along with land use change played an important role in decrease of runoff from the Soan basin, Pakistan. Thus, understanding the short and long-term effects of climate change on parameters such as temperature, precipitation, baseflow, surface runoff, and total streamflow at a regional scale becomes essential, especially for the local authorities and decision-makers (Shahid & Rahman, 2021;Wang et al., 2016). The response of a local catchment to the climatic changes (e.g., variation in precipitation and temperature) is represented in its hydrologic regime. Also, the geomorphologic evolutions of catchments are much slower than the possible climatic variations being experienced. Thus, in the pristine watersheds with the unregulated flows, a stream's hydrologic regime can be a good indicator of the effects of climate change (Zhang et al., 2001).
The rate of increase in the annual mean temperature is almost twice the global warming rate in the northern regions of the planet (ACIA, 2005;AMAP-SWIPA, 2011). For the past few decades, Canada has experienced a rapid rise in the annual mean surface air temperature; there has been an increase of 1.5˚C between years 1950 and 2010 (Vincent et al., 2012). A warming shift has also been observed for the Hudson Bay area from the year 1998 (Fazel-Rastgar, 2020). The province of Ontario, which is the most populous and the largest economy among all Canadian provinces, has been impacted considerably by regional climate change. For instance, in July 2013, a flood in Toronto caused the largest natural disaster ever recorded in the history of the province (Wang et al., 2014).
Based on the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change, heavy precipitation followed an increasing trend between the mid-20th century and early 21st century in North America (Barros et 2014). Studies conducted towards the end of the 20th century showed that in the last five decades of the century, the changes in annual precipitation within Canada have been between −10% to +35%. The northern regions have experienced the strongest changes, whereas in south-eastern Canada, a significant decrease in winter precipitation was observed (Zhang et al., 2001). In higher latitudes, basins with snowmelt-flooding events are expected to witness changes in their flood frequencies (Burn & Elnur, 2002).
For the past few decades, many researchers have studied the impacts of climate change on hydrology around the globe and within Canada. However, the results of these studies vary a lot depending on the location of the study; these results include both negative and positive trends in the hydrologic and climatic patterns (Nalley et al., 2012). Several studies have indicated increasing mean annual stream flows in Canada (Burn & Elnur, 2002;McBean & Motiee, 2006;Ehsanzadeh et al., 2007;Vincent et al., 2015), corresponding to increasing trends in temperature and precipitation (Barros et al., 2014;Wang et al., 2016;Whitfield & Cannon, 2000). However, only a few studies have examined the long-term changes in surface runoff and baseflow of the streams in Ontario; these two hydrologic parameters correspond to different characteristics within each watershed, and their analyses can present a better view of the long-term impacts of changing climate on the hydrology of the watershed. There is also a need to investigate long term trends of climatic and hydrologic parameters in the urban watersheds.
Global climate change presents a real threat to human health, ecosystems, and water management systems; there is an urgent need to develop mitigation and adaptation strategies at local and global levels (Upadhyay, 2020). To be able to take any further steps in preparation and implementation of new strategies to mitigate the undesired consequences of the ongoing regional climate change in Ontario, it is essential to perform an updated and comprehensive study that analyzes the historical trends in hydrology and climatology of the province for the past few decades; the changing climate and its effects on hydrological processes can be assessed through a trend analysis of the location's historical climatic and hydrologic data (Shahid & Rahman, 2021). Thus, the objective of this study is to understand the effects of the changing climate on hydrological process for rural and urban watersheds situated in northern and southern parts of Ontario.

Data Selection
In this study, we have looked at 13 rural and 2 urban watersheds within Ontario.
Due to rapid urbanization in major Canadian cities, investigating the impacts of land use changes in urban watersheds on hydrologic responses in Ontario is one of the objectives of this study. For this purpose, two watersheds in the Greater Toronto Area (GTA) were also selected. The gross drainage areas for all the se- tions and are considered pristine with less than 5% of their area been altered (Cunderlik & Burn, 2004;Cunderlik & Ouarda, 2009 (Burn & Elnur, 2002). In this research, the minimum duration of the time period for which the data is available was 36 years; however, for most of the stations, 60 years of data was considered.

Baseflow Separation
To have a better understanding of variations in the streamflow, it is essential to separate the baseflow from the surface runoff. Various strategies can be followed to separate baseflow from surface runoff; these include graphical separation methods, filtering methods, frequency analysis methods, and recession analysis methods. One of the most commonly used methods for separating baseflow from the surface runoff is the one-parameter filter method proposed by Lyne and Hollick (1979). The results of a study performed by Rudra et al. (2015), showed that the digital filter method can generate more realistic baseflow hydrographs in the studied watersheds in southern Ontario. The equation for the one-parameter digital filter method used in the present study is given by where, y represents the total streamflow, b represents baseflow, and k represents the time-step number. Parameter a is the filter parameter with a suggested value of 0.925; using this value in this equation gives realistic results (Arnold et al., 1995;Nathan & McMahon, 1990).

Trend Analysis
The linear regression method is one of the simplest techniques for estimating the slope of a trend line in the data series. Linear least-square method is the most used technique for linear regression. However, very low and very high (extreme) values can considerably affect the results; this method, therefore, is not suitable for trend analysis of hydrologic and climatic parameters where the data consists of these low and high values. However, due to the simplicity of this method, it is one of the simplest ways to visualize the presence of trends in the data-series charts.
Man-Kendall test (Kendall, 1975;Mann, 1945) is also one of the most widely used trend tests to estimate nonparametric linear trends in hydrologic and climatic time-series (Cunderlik & Ouarda, 2009 (Cunderlik & Ouarda, 2009;Zhang et al., 2001). One of the most frequently used approaches to overcome this problem is the trend free pre-whitening (TFPW) method. The method was used to remove serial correlation by using Equation (2) (Cunderlik & Ouarda, 2009;Von Storch & Navarra, 1995;Zhang et al., 2001): where, t x′ is the pre-whitened value of the time-series, x t is the original value of the time-series, and AC 1 is the lag-1 autocorrelation coefficient. The pre-whitening method is applied when there is a serial correlation in the time-series at a significance level of 0.05.
Sen's non-parametric method was applied to estimate the slope of time-series.
This method assumes that the data time-series has a linear trend. The slope (Q) of the trend can be calculated based on the slopes of all pairs of data. A positive value for Q shows an increasing trend and a negative value for Q represents a decreasing slope in the data-series. The Excel program MAKESENS (Salmi et al., 2002) was used to find the Sen's non-parametric slope of the data-series.

Results
Data from all flow, precipitation and temperature parameters collected from the hydrometric and climate stations were analyzed to identify trends. The results of the analysis for each of the rural and urban watersheds are presented in Table 2 and Table 3 respectively; statistically significant (p < 0.05) trends are indicated within the table. The time periods for the trend analyses were selected in a way that covered the maximum common duration for the data collected from different stations. Some of the climate stations provided data which did not lie between the selected duration; such data was not included in the current analysis.
However, due to availability of the data, it is important to note that the time period selected for the comparison of the flow from different watersheds are slightly different from the common time period selected for comparing the climatic trends within each watershed.

Rural Watersheds
A total of 13 watersheds in the rural area across Ontario were investigated. Table   2 presents the results of the trend analysis for these selected watersheds. The watershed corresponding to the Sydenham river near Alvinston has two climate stations. However, as per the common time duration of the data provided by the hydrometric station of this watershed, only the Strathroy climate station was selected for analysis. Only the month of October was found to have an increasing trend in precipitation, and the rest of the months presented no significant trends. Also, no significant trends in the hydrometric parameters were observed  The upward arrow (↑) signifies a significant increasing trend, and the downward arrow (↓) signifies a significant decreasing trend (α = 0.05). The asterisk symbol (*) signifies presence of considerable number of "zero" values in the data that did not permit us to perform the trend analysis for the specified month.
A. Azarkhish et al. The upward arrow (↑) signifies a significant increasing trend, and the downward arrow (↓) signifies a significant decreasing trend (α = 0.05). The asterisk symbol (*) signifies presence of considerable number of "zero" values in the data that did not permit us to perform the trend analysis for the specified month.
for all the months. It was found that, in general, the minimum and maximum temperature was found to increase in the summer months for this watershed. It should be noted that although this station is not a part of the Reference Hydro- for the watershed of Goulais river near Searchmont. Unlike other watersheds, data from the hydrometric station 02BF002 suggested no increasing trend for the baseflow for this watershed during the winter months. On the contrary, the mean monthly baseflow was found to decrease during the summer (Jun-Aug). Unfortunately, climate data was not available for this watershed, for the selected duration. The decadal average of monthly streamflow for the period of 2008 to 2017 is more than the average of monthly streamflow for the period of 1998 to 2007 for all of the months except for the month of April.
The mean monthly surface flow for the watershed of Neebing River near Thunder Bay was found to increase during the month of July and decrease significantly for the month of December. For the same watershed, the minimum temperature for the month of October was found to decrease over the period of past 36 years, and consequently the amount of mean monthly snowfall ws found to increase for the same month. This was not in accordance with the results for other watersheds considered for this study. A duration of 36 years  was chosen to perform the trend analysis. The decadal average monthly streamflow for the period of 2008 to 2017 is more than the average of monthly streamflow observed between 1998 to 2007 for all of the months except for the months of April and October. The average monthly streamflow and baseflow were found to decrease during the fall months (Aug-Oct) for the watershed of Missinaibi river at Mattice across a period of 1958 to 2017. The average monthly surface flow was found to significantly increase in March and decrease during the months of May of September. The observed surface runoff values are mostly "Zero" during the winters (Dec-Mar) due to icy conditions. No climate data was available for the watershed for the given time period. The decadal average monthly streamflow for the period of 2008 to 2017 is more than the average of monthly streamflow for the period of 1998 to 2007 for all of the months except for the months of April and October.
Trend analysis of various parameters considered across 38 years  for the watershed of Nagagami river at Highway No.11 revealed that the average monthly streamflow and baseflow decreased during the month of June. Also, a decreasing trend in the monthly average maximum temperature for October was observed. Unlike the other watersheds, the decadal average monthly streamflow for the period of 2008 to 2017 was not found to be higher than the average monthly streamflow of the decade before that. Also, as compared to other watersheds considered within the study, the decadal average monthly streamflow was found to be higher in the month of May; the month of April had the highest decadal average monthly streamflow for other watersheds.
The mean monthly streamflow and baseflow were found to increase during   (Figure 3). Figure 4 shows the decadal average monthly streamflow for the urban watersheds; even with the huge amount of land use change in the past decades, the Etobicoke

Urban Watersheds
Creek station still receives unregulated flow (Environment and Climate Change Canada, 2012) and thus, it can be used as a measure to study the effects of urbanization on hydrometric and climatic parameters. No significant increasing or decreasing trends in the average monthly precipitation for the past fifty years (1968 to 2017) was observed within the watershed (Figure 4(a)). However, the average monthly minimum temperature increased significantly for almost every month. Also, the streamflow, baseflow and the surface runoff seemed to increase during the months of May to October. One of the possible reasons for this significant increase in the surface runoff during these months is the effect of urbanization, leading to reduced infiltration rate. In order to investigate the significant increase in baseflow, without an increase in precipitation, snowfall, and rainfall for the watershed, baseflow separation from the surface runoff was studied using the one-parameter filter method. This technique separates the slow flow from the quick flow using the filter method. The major components of the slow flow in urban watersheds are considered to be: 1) the classical baseflow from groundwater inflow to the stream, and 2) the drainage flow from the detention ponds that are created for storm-water management purposes (Liu et al., 2013).
The Greater Toronto Area (GTA) has developed rapidly during the last 50 years, and a complex network of storm-water management ponds and canals have been designed to control flooding in the newly developed urban area. Before urbanization, the infiltrated water was eventually discharged into the streams, lakes, or neighboring watersheds. But, after development as the permeable layer decreased at the surface, larger amounts of water that was supposed to infiltrate, stayed at the surface. To control flooding, many storm-water management ponds A. Azarkhish et al. have been designed, thus changing the shape of the hydrograph in a way that the slow flow increases.
The Don River at Todmorden watershed is 319 km 2 in size and its station 02HC024 receives regulated flow from the river. As the streamflow and the baseflow data have a significant lag 1 autocorrelation, the flow data from this sta-  Figure 4(b) shows the decadal average of the monthly streamflow of the Don River at Todmorden; with passing decades, the average monthly streamflow was found to increase significantly.
the trends of climatic and hydrometric parameters can be found for watersheds within a given region of the province. For investigating these similarities, all watersheds were sorted according to their latitude and land use. Significant trends in the streamflow, baseflow, and surface runoff, were then calculated for a common period of time  for all of the watersheds. Table 4 shows the significant streamflow trends for rural watersheds located in the southern region of the province (below 45 degrees latitude, north). According to the trend analysis, all of the RHBN rural watersheds exhibited a significant increase in the streamflow during the month of January. Correspondingly, the baseflow for these watersheds was also found to increase during the month of January, and for some cases in February. The baseflow of all of the rural watersheds is significantly increasing in the month of January and for some cases in February. Since the surface runoff has can be impacted by changes in climatic parameters more than baseflow, the pattern of significant trends of surface runoff in this region is slightly different than the pattern of significant baseflow for the same period. Surface runoff was mostly found to increase during the month of January.   Table 2.

Watersheds in Northern Ontario
Based on the discussion of the results, an urgent need for local authorities and decision makers to propose and develop new water management strategies emerges. These strategies should be designed by taking into account the regional variations of the impact of climate change on hydrology of watersheds within Ontario.

Conclusion
Significant upward trends were observed for total streamflow, baseflow and surface runoff during the winter and summer months in watersheds situated in southern Ontario. This was mainly attributed to the increasing trends for the monthly average minimum temperatures for this time period. The observed in-A. Azarkhish et al. American Journal of Climate Change creasing trend in the minimum and maximum temperatures for the winter months can speed up the process of snowmelt and alter the hydrological processes for these watersheds. Unlike the southern Ontario watersheds, very few instances of significant changes in the mean monthly temperatures (minimum and maximum) were observed for the watersheds located in northern part of Ontario; in some cases, decreasing trends were also observed. Similarly, very few instances of significant decrease in the baseflow during summer months were also observed for some of the watersheds in the north (except North French River watershed). The results of this study also indicate that in the selected urban watersheds, the average monthly minimum and maximum temperatures are increasing for almost all of the months. Also, the streamflow, baseflow, and surface runoff were found to increase for the Etobicoke watershed during the summer and fall season, possibly due to rapid urbanization, causing reduced infiltration rate.
According to the findings of this study, local and regional decision-makers must soon propose and adopt new water management strategies, particularly in light of changing climate and land use patterns. These methods should be developed with geographic differences in the effects of climate change on the hydrology of Ontario's watersheds in mind.