Projected Changes in Hydrological Variables in the Agricultural Region of Alberta, Canada

The responses of regional hydrological variables to climate change are of prime concern for agricultural water resources planning and management. Therefore, the seasonal (April–September) and annual (January–December) evolution of precipitation, temperature, evapotranspiration (ET), soil moisture (SM), deep aquifer recharge (DA), and water yield (WYLD) was investigated using established statistical techniques for the historical, near and far future (1983–2007: His, 2010–2034: NF, 2040–2064: FF) in the agricultural region of Alberta, Canada. Previously calibrated and validated agro-hydrological models (Soil and Water Assessment Tool) were used to generate these variables. Future changes were investigated under two representative concentration pathways, i.e., RCP 2.6 and RCP 8.5, projected by nine global climate models (GCM). Results revealed that Alberta had become warmer and drier during the His period. The future projection showed an increase in precipitation, SM, DA, and WYLD, in turn, indicated more water resources. Precipitation and temperature were projected to increase between 1 to 7% and 1.21 to 2.32 ◦C, respectively. Seasonal precipitation showed a higher trend magnitude than that of annual precipitation. The temperature generally had an increasing trend in the future with a maximum in the southern Alberta. Monthly average ET was likely to increase and decrease in the rising and falling limbs of the bell-shaped curve with the peak in July. A comparison of water demand from two land use types (dominant land use and barley) during the His period showed that water deficit existed in July and August. The results of this study could help in understanding anticipated changes in hydrological variables and decision-making regarding the regional agricultural water resources management.


Introduction
Global climate warming is undoubtedly happening due to the anthropogenic actions which include burning fossil fuels, deforestation and rapid urbanization [1]. It modifies the intensity and frequency of precipitation, in turn, changes the hydrological cycle. This change is likely to increase in the future [2] and will continue to pose a threat to the water availability for different sectors, e.g., agriculture [3]. Soil water availability in an agricultural region is one of the primary inputs to raise agricultural production to feed a growing world population. Precipitation is a key driver controlling soil water in most of the rainfed croplands, whereas temperature is particularly important in most snow-dominated and semi-arid areas [3] as in Canada. The sensitivity of semi-arid regions to changes in precipitation and temperature increases regional vulnerability to potential effects of climate change on water resources and agriculture. Therefore, it is essential to investigate the changing properties of agro-hydrological variables including precipitation, temperature, evapotranspiration, soil water, deep aquifer recharge, and water yield due to climate change. The province has 17 river basins and 19 census divisions (CD) of agriculture with 43,234 agricultural farms. The agricultural region of Alberta is spread out in 11 out of 17 river basins ( Figure 1) that are chosen to assess the spatial variability of hydrological parameters in this study. A total of seven counties from seven agricultural CD were selected based on the number of farms for trend analysis (Table 1). First, the CD with the total number of farms in the division more than 2000 was selected. Next, the county in each CD with a total number of farms >900 was selected for trend analysis. No county had been selected having less than 900 farms, even if the CD had more than 2000 farms ( Table 1). The name of these counties are-(a) Grande Prairie from Grande Prairie CD, (b) LacStAnne from Barrhead/Athabasca CD, (c) Leduc from Edmonton CD, (d) Lethbridge from Lethbridge CD, (e) Mountain View from Calgary CD, (f) Red Deer from Red Deer CD, and (g) Vermilion from Camrose/Vermilion CD (Figure 1). The selected counties represent various hydro-climate, geospatial, and management conditions in the province's agricultural lands.
Based on the best available land use-land cover map (see Section 2.3), the dominant land use types were different across counties. For example, Agricultural Land-Row Crops (AGRR) and Forest-Deciduous (FRSD) land use types were found to be dominant in the Grande Prairie county, whereas Pasture (PAST) and Spring Wheat (SWHT) were considered as the main land use types in LacStAnne, Leduc and Lethbridge County, respectively ( Table 1). The dominant land use type in the rest of the counties were AGRR. The overall dominant land use type in the agricultural region of Alberta was AGRR followed by FRSD. The province has 17 river basins and 19 census divisions (CD) of agriculture with 43,234 agricultural farms. The agricultural region of Alberta is spread out in 11 out of 17 river basins ( Figure 1) that are chosen to assess the spatial variability of hydrological parameters in this study. A total of seven counties from seven agricultural CD were selected based on the number of farms for trend analysis (Table 1). First, the CD with the total number of farms in the division more than 2000 was selected. Next, the county in each CD with a total number of farms >900 was selected for trend analysis. No county had been selected having less than 900 farms, even if the CD had more than 2000 farms ( Table 1). The name of these counties are-(a) Grande Prairie from Grande Prairie CD, (b) LacStAnne from Barrhead/Athabasca CD, (c) Leduc from Edmonton CD, (d) Lethbridge from Lethbridge CD, (e) Mountain View from Calgary CD, (f) Red Deer from Red Deer CD, and (g) Vermilion from Camrose/Vermilion CD (Figure 1). The selected counties represent various hydro-climate, geospatial, and management conditions in the province's agricultural lands.
Based on the best available land use-land cover map (see Section 2.3), the dominant land use types were different across counties. For example, Agricultural Land-Row Crops (AGRR) and Forest-Deciduous (FRSD) land use types were found to be dominant in the Grande Prairie county, whereas Pasture (PAST) and Spring Wheat (SWHT) were considered as the main land use types in LacStAnne, Leduc and Lethbridge County, respectively ( Table 1). The dominant land use type in the rest of the counties were AGRR. The overall dominant land use type in the agricultural region of Alberta was AGRR followed by FRSD.

Data Collection
The study period was divided into three different time horizons (a) historical (His: 1983 to 2007), (b) near future (NF: 2010 to 2034), and (c) far future (FF: 2040 to 2064). The climate data were collected from 300 meteorological stations at a daily time step from Environment and Climate Change Canada (http://climate.weather.gc.ca/) for His period. For NF and FF, climate change data for nine Global Climate Models (GCMs) namely-CanESM2, CCSM4, CNRM-CM5, CSIRO-MK3, GFDL-ESM2G, HADGEM2, MIROC5, MPI-ESMLR, and MRI-CGCM3, were collected from the Pacific Climate Impacts Consortium (PCIC; https://www.pacificclimate.org/data) under two RCP scenarios, i.e., RCP 2.6 and 8.5. The selected GCMs are the best model combination out of the majority of IPCC collection that represents the largest plausible future in Canada [22]. The calibrated hydrological model [20,23] was used to simulate hydrologic water balance variables for each climate model and scenario. Owing to uncertainties in the detailed representation of many complex climate processes, individual climate models vary in their representation of these processes and will have biases of various kinds [24]. Hence, it is suggested to use multi-model ensemble average of projections, which produces smaller errors than any individual model [24,25]. Therefore, results of this study were presented as the ensemble average of nine GCMs in this study. In this study, months from April to September were considered to explain the results in 'seasonal' scale, whereas January to December was considered for the 'annual' scale.

The SWAT Model
The Soil and Water Assessment Tool (SWAT) is a physical process-based, time-continuous hydrological model [26]. This model simulates hydrological variables within small or large watersheds under different management, land use change, and climate change scenarios for impact assessment. Based on topographic data, soil characteristics and land use type, the model delineates watersheds into unique spatial units and predicts future hydrological variables using a water balance equation in each unit ( Figure 2) [27]: where, SW t = final soil water content (mm); SW 0 = initial soil water content on day i (mm); t = time (days); R day = amount of precipitation on day i (mm); Q sur f = amount of surface runoff on day i (mm); E a = amount of evapotranspiration on day i (mm); w seep = amount of water entering the vadose zone from the soil profile on day i (mm); Q gw = amount of return flow on day i (mm). The SWAT model also simulates potential evapotranspiration (PET). Three methods have incorporated into SWAT including the Penman-Monteith, Priestley-Taylor and Hargreaves method to calculate PET. We used the Penman-Monteith method as follow [28]: where is the latent heat flux density (MJ m −2 d −1 ), is the depth rate of evaporation (mm day −1 ), Δ is the slope of the saturation vapor pressure-temperature curve, / (kPa °C −1 ), is the net radiation (MJ m −2 day −1 ), is the heat flux density to the ground (MJ m −2 day −1 ), is the air density (kg m −3 ), is the specific heat at constant pressure (MJ kg −1 °C −1 ), is the saturation vapor pressure of air at height z (kPa), is the water vapor pressure of air at height z (kPa), is the psychrometric constant (kPa °C −1 ), is the plant canopy resistance (s m −1 ), and is the diffusion resistance of the air layer (aerodynamic resistance) (s m −1 ). The model simulates actual evapotranspiration (AET) based on PET, soil water availability, and the maximum amount of transpiration depending on the plant type and related daily above-ground and below-ground biomass production. This AET is used in this study and referred to evapotranspiration (ET) afterward.
The model works on the concept that the daily precipitation can generate surface runoff, and a fraction can infiltrate into the soil depending on land use and soil characteristics. Further, a fraction of daily soil water becomes available for plant evapotranspiration. It will also partition into a subsurface flow to stream and recharge to the shallow aquifer that will partially feed streamflow. The recharge can be routed to the deep aquifer. The amount of water will be diverted from the shallow aquifer due to percolation to the deep aquifer is the total recharge entering the aquifer multiplied by the aquifer percolation coefficient. The water yield in the SWAT model is the amount of water leaving a hydrologic response unit and entering the main channel during the simulated time-step. For more details see Neitsch et al. [28].
The model was built using the digital elevation model at 10-m resolution [29], land use-land cover maps and related physical parameters (http://www.geobase.ca/geobase/en/data/landcover/ csc2000v/description.html) and soil map with the soil physical properties embodied in the map (http://sis.agr.gc.ca/cansis/nsdb/slc/index.html) as well as agro-hydrological data [20,23]. A total of 2255 sub-basins were delineated for Alberta [23], however, in this study we extracted the sub-basins that are located within the agricultural region. An extensive calibration and validation of the The SWAT model also simulates potential evapotranspiration (PET). Three methods have incorporated into SWAT including the Penman-Monteith, Priestley-Taylor and Hargreaves method to calculate PET. We used the Penman-Monteith method as follow [28]: where λE is the latent heat flux density (MJ m −2 d −1 ), E is the depth rate of evaporation (mm day −1 ), ∆ is the slope of the saturation vapor pressure-temperature curve, de/dt (kPa • C −1 ), H net is the net radiation (MJ m −2 day −1 ), G is the heat flux density to the ground (MJ m −2 day −1 ), ρ air is the air density (kg m −3 ), c p is the specific heat at constant pressure (MJ kg −1 • C −1 ), e o z is the saturation vapor pressure of air at height z (kPa), e z is the water vapor pressure of air at height z (kPa), γ is the psychrometric constant (kPa • C −1 ), r c is the plant canopy resistance (s m −1 ), and r a is the diffusion resistance of the air layer (aerodynamic resistance) (s m −1 ). The model simulates actual evapotranspiration (AET) based on PET, soil water availability, and the maximum amount of transpiration depending on the plant type and related daily above-ground and below-ground biomass production. This AET is used in this study and referred to evapotranspiration (ET) afterward.
The model works on the concept that the daily precipitation can generate surface runoff, and a fraction can infiltrate into the soil depending on land use and soil characteristics. Further, a fraction of daily soil water becomes available for plant evapotranspiration. It will also partition into a subsurface flow to stream and recharge to the shallow aquifer that will partially feed streamflow. The recharge can be routed to the deep aquifer. The amount of water will be diverted from the shallow aquifer due to percolation to the deep aquifer is the total recharge entering the aquifer multiplied by the aquifer percolation coefficient. The water yield in the SWAT model is the amount of water leaving a hydrologic response unit and entering the main channel during the simulated time-step. For more details see Neitsch et al. [28].
The model was built using the digital elevation model at 10-m resolution [29], land use-land cover maps and related physical parameters (http://www.geobase.ca/geobase/en/data/landcover/ csc2000v/description.html) and soil map with the soil physical properties embodied in the map (http://sis.agr.gc.ca/cansis/nsdb/slc/index.html) as well as agro-hydrological data [20,23]. A total of 2255 sub-basins were delineated for Alberta [23], however, in this study we extracted the sub-basins that are located within the agricultural region. An extensive calibration and validation of the hydrologic and crop models were conducted [20,21,23]. Two statistical measures were used to quantify the goodness-of-fit and model output uncertainty using the Sequential Uncertainty Fitting (SUFI-2) program. The SUFI-2 program is linked to SWAT and provides the basis for parallel processing of calibration and large-scale parameterization schemes [20,21,23]. Two statistical measures are p-factor, which is the percentage of observed data bracketed by the model outputs uncertainty, quantified as 95% prediction uncertainty (95PPU), and r-factor, which is the average thickness of the 95PPU band, were used to check the model performance during calibration and validation. The 95PPU is calculated at 2.5% and 97.5% levels of the cumulative distribution of an output variable [20,21,23]. The r-factor is the average width of the 95PPU band divided by the standard deviation of the observed variable. In an ideal condition, a p-factor of hundred percent and r-factor of zero is the simulation that exactly brackets the observed data. Calibrated-validated models were used to simulate historical and future water balance components for further analysis in this study.

Assessment of Hydrological Variables
Simulated daily records of agro-hydrological variables (temperature, precipitation, evapotranspiration, soil moisture, deep aquifer recharge, and water yield) were aggregated into seasonal (April-September) and annual (January-December) means. Long-term monthly mean also used to examine monthly changes in hydrological variables. Long-term mean has been used by hydrologists and climatologists in different regions of the world [30,31].

Trends of Hydrological Variables
A variety of statistical methods can detect the trends of hydrological variables at different temporal scale. Both non-parametric (Mann-Kendall test) [32,33] and parametric (linear regression analysis) procedures are used to identify the annual and seasonal trends. The Mann-Kendall (MK) statistical test is used worldwide to detect the trend of precipitation, temperature, and relative humidity [14,[34][35][36][37]. The non-parametric tests are preferred as they are more suitable for non-normally distributed data, which are frequently encountered in hydro-meteorological time series [38]. Although the original MK test has been widely used and recommended by the World Meteorological Organization, it fails to deal with the issue of autocorrelation in the hydro-meteorological time series. Hamed and Rao [38] improved the MK test by taking the lag-i autocorrelation into consideration; and the modified MK (MMK) has been shown to be robust in capturing the trends of the hydro-meteorological time series [39]. Chen et al. [15] also found that the MMK test's performance is better than the original MK test if data is autocorrelated. In this study, the trend obtained by the MMK test was assessed at the significance level of α = 5%.
In addition to the trend detection, it is necessary to estimate the trend magnitude. Sen's nonparametric method [38] is used to estimate the true slope of an existing trend (as change per year). Sen's slope method gives a robust estimation of the trend and is used in various trend analyses of hydrological variables [38]. Hence, the non-parametric statistical trend analysis is recommended to assess the impact of climate change on the hydrological variables. A 100 (1 − α)% two-sided confidence interval about the slope estimate is obtained by the nonparametric technique based on the normal distribution. A positive value of trend estimation indicates an upward trend and a negative value indicates a downward trend in the time series.

Potential Land Use Type and Water Source
The ET was simulated by SWAT agro-hydrologic models for dominant land use (ET_DLU) and barley crop (ET_barley) for His period [20,21]. These two ET were used for comparison of the water demand differences in a county for two land use types. We intended to get a sense of how large the water demand of current dominant land use type is in a county as compared to a potential crop type. We also calculated water deficit (WD) for His period by taking the difference between precipitation and ET for each land use type (ET_DLU and ET_barley). The WD can provide useful information on how the water needs of an existing land use type are compensated by precipitation. For example, in a given county, large negative delta (i.e., precipitation minus ET) for a given land use type indicates that the plant water needs were from either irrigation or other sources than the natural soil water from precipitation.

The SWAT Model Performance
We used outputs from the SWAT hydrological and crop simulation models in this study. A multi-gauge and multi-objective calibration using river discharges and crop yield ensured proper apportioning of precipitation and soil water into surface runoff, actual evapotranspiration, and groundwater recharge. Altogether, measured discharge data of 130 hydrometric stations and annual yield data of 67 counties were used to calibrate and validate the agro-hydrologic model. The model performance was satisfactory for most of the hydrometric stations and counties. Overall, 63% and 71% of the observed streamflow data, and 88% and 85% of observed crop yield data were captured by the model uncertainty prediction for calibration and validation periods, respectively. The average r-factor was about 1.04 and 1.43 for discharge and 4.48 and 5.35 for crop yields during the calibration and validation, respectively at the Alberta scale (Table 2). For yield simulation in the best performing counties, the model captured 79% of observed data with an r-factor of 2.51, while the model captured 95% of observed data with an r-factor of 6.85 in the worst case. Larger r-factor in the crop model might be due to climate variability and unreliable precipitation distribution in Alberta, and this has increased standard deviation as well as the thickness of the uncertainty band. Detailed performance results of the hydrological and crop simulation models can be found in Faramarzi et al. [20] and Masud et al. [21], respectively.

Seasonal and Annual Variation of Hydrological Variables
We examined seasonal and annual changes of hydrological variables for the agricultural region of Alberta. The study was performed for the historical (1983-2007) and future (2010-2064) periods projected by nine GCMs under two RCP scenarios. In this section, we presented and discussed results based on the ensemble average of all GCMs. Figure 3 shows historical and projected future seasonal and annual precipitation distribution. Historical precipitation distribution showed that southeastern part was drier than any other area of the agricultural region. The precipitation was high in the west and followed by the central part of the region. Annual and seasonal precipitation ranged from 280 to 745 mm and 208 to 493 mm, respectively. Projected changes in annual and seasonal precipitation in both NF and FF indicated a likely increase compare to the His, which was higher for FF followed by NF. Most of the sub-basins were likely to gain higher precipitation during NF as compared to His, while some of them were likely to receive less precipitation. In the FF, precipitation was projected to increase uniformly over the study area having a maximum in the eastern part. projected to increase with high precipitation. These findings are in agreement with some small-scale studies in southern Alberta [41,43,45]. Spatial distribution of deep aquifer recharge and water yield is shown in the supplementary materials ( Figures S3 and S4). According to Falkenmark and Rockstrom [46], the summation of water yield (WYLD) and deep aquifer recharge (DA) is termed as blue water. The spatial pattern of blue water components was similar to that of precipitation. Out of the 440-mm mean annual precipitation in His period, about 64 mm was renewable blue water resources in the agricultural region. Based on the projected changes, the blue water resources were expected to marginally increase in most of the sub-basins.  Annual precipitation was projected to increase by 3% and 7% in most of the sub-basins during the NF and FF. Similarly, the seasonal precipitation could increase by 1% and 5% for two future horizons. It is worth to mention that projected increase in southern sub-basins will not result in a significant amount of precipitation since the historical precipitation is minimum in this area.
The variation of seasonal and annual mean temperature is presented in Figure 4, and the maximum and minimum temperature maps are shown in the Supplementary Materials ( Figures S1 and S2). The southern region was hotter with the lowest amount of precipitation for all periods while the northwestern part was colder with the highest amount of precipitation. The average annual mean, maximum, and minimum temperatures were found to be 5, 9, and −0.2 • C during the His period. Likewise, the average seasonal mean, maximum, and minimum temperatures were found to be 13, 19, and 7 • C. Compare to the His period, the temperature was rising in NF and FF with a maximum increase in FF (Figure 4). The spatial distribution of the temperature rise was different for the annual and seasonal case. The largest temperature increase was projected for the eastern and southern part of the region for the annual and seasonal case, respectively. In general, the average rise of mean annual and the seasonal temperature was 1.33 and 1.21 • C in NF, while it was 2.32 and 2.14 • C in FF, respectively. The spatial patterns of maximum and minimum temperature for His and projected changes for the NF and FF were similar to the spatial pattern of mean temperature. This GCM projected precipitation and temperature changes agreed with previous studies [40][41][42][43][44]. These results are also supported by the recent results reported by Environment and Climate Change Canada [24] in a country-wide study.
The spatial distribution of ET and soil moisture (SM) is shown in Figures 5 and 6. The projection suggested larger increase (decrease) in ET (SM) in the sub-basins, where the temperature was also projected to increase with high precipitation. These findings are in agreement with some small-scale studies in southern Alberta [41,43,45]. Spatial distribution of deep aquifer recharge and water yield is shown in the Supplementary Materials ( Figures S3 and S4). According to Falkenmark and Rockstrom [46], the summation of water yield (WYLD) and deep aquifer recharge (DA) is termed as blue water. The spatial pattern of blue water components was similar to that of precipitation. Out of the 440-mm mean annual precipitation in His period, about 64 mm was renewable blue water resources in the agricultural region. Based on the projected changes, the blue water resources were expected to marginally increase in most of the sub-basins.

Seasonal and Annual Trends of Hydrological Variables
Trend of seasonal and annual precipitation and temperature were analyzed for the current and future periods. The magnitude of the trend and its significance in the time series were determined using the Sen's slope estimator and the MMK test ( Table 3). The trend analysis was performed separately for RCP 2.6 and 8.5 in seven selected counties (Table 1) in the Alberta agricultural region for His, NF and FF.
During His period, annual and seasonal precipitation showed a decreasing trend in all counties,

Seasonal and Annual Trends of Hydrological Variables
Trend of seasonal and annual precipitation and temperature were analyzed for the current and future periods. The magnitude of the trend and its significance in the time series were determined using the Sen's slope estimator and the MMK test (Table 3). The trend analysis was performed separately for RCP 2.6 and 8.5 in seven selected counties (Table 1) in the Alberta agricultural region for His, NF and FF.
During His period, annual and seasonal precipitation showed a decreasing trend in all counties, except in Lethbridge and Red Deer, where seasonal precipitation showed an increasing trend (Table 3). However, Grande Prairie county only observed the statistically significant negative trend. Likewise, non-significant positive and negative annual trends were found for NF and FF, which was similar to the His period in most of the counties. Only Grande Prairie observed a statistically significant negative trend in NF for RCP 2.6. In general, seasonal precipitation showed an upward or downward trend with a higher magnitude than annual precipitation for both RCPs.
The trend magnitude of both annual and seasonal mean temperature varied slightly among different counties, and it ranged from 0.03-0.05 • C (annual) and 0.01-0.02 • C (seasonal). Table 3 also indicated that the trend of mean temperature generally increased from southern to northern Alberta. In other words, the region with lower mean temperature ( Figure 4) in His period had a higher warming rate in the future. The increasing trend of mean temperature was observed in the NF and FF. However, most of the counties showed non-significant trend except Mountain View, where the trend was found to be significant. The trend magnitude of annual mean temperature was usually higher than the seasonal mean temperature. Similar trends were found for both maximum and minimum temperature for the historical as well as future periods (Table 3). Jiang et al. [42] investigated historical and future precipitation and temperature of Alberta using the Special Report on Emissions Scenarios (A1B, A2, and B1) from an ensemble of GCMs. They found future June-July-August precipitation had downward and temperature had upward trends. The magnitude of precipitation and temperature are close to our findings. Zhang et al. [4] also found similar results for Alberta at their Canada-wide study for the historical period of the 20th century. These changes of temperature and precipitation are attributed, at least in part, to the increasing atmospheric greenhouse gases concentration which is mainly due to anthropogenic activities [11,47].
Trend analysis of other hydrological variables in a changing climate is vital to assess climate changes and suggests appropriate water resource management strategies, especially in the agricultural region for the future. Trend analysis of ET, SM, DA, and WYLD are given in Table 4 for the seven counties, and results of these variables are somewhat similar for both RCPs. ET showed a non-significant decreasing trend in all counties during His period, and this trend continued in the future. However, there was a positive trend for annual ET in Grande Prairie county. Similar to the temperature, the trend magnitude of annual ET was higher than that of seasonal magnitudes during the His, NF, and FF. The annual SM during His period showed a decreasing trend in most of the counties except Leduc and Red Deer. For seasonal SM in the His period, most of the counties showed positive trend except Grande Prairie and Mountain View, where the trend was negative. In the future, the annual SM showed a negative trend for both RCPs. In contrast, the seasonal trend in the future was found to be positive for most of the counties. There were some counties which observed a statistically significant trend in annual and seasonal SM for RCP 8.5 only, including Lethbridge (annual SM), Red Deer and Vermilion (seasonal SM). The trend magnitude of DA was very mild. The slope was zero for some of the counties and time horizons. These results indicate a relatively constant level of DA with respect to other hydrological variables during the study periods. In all time periods (His, NF and FF), Grande Prairie showed a non-significant negative trend of DA. Seasonal and annual WYLD trends showed a positive magnitude in most of the counties during the His period (Table 4).  Figure 7 shows the long-term average monthly variation of hydrological variables. The results showed that monthly precipitation, temperature, ET, and SM had a distinct and similar pattern for current and future periods. Precipitation, temperature, and ET followed the bell-shaped curve with the maximum during the agricultural growing season (May-August). Precipitation was projected to increase during the winter and spring seasons (December to May). However, it was likely to show a mixed pattern for the summer and fall seasons (June to November). Projected changes showed a likely uniform increase in temperature throughout the year (Figure 7). Overall, ET was likely to increase and decrease in the ascending and recession limbs of the bell-shaped curve, respectively except for Lethbridge county (Figure 7). SM was projected to decrease in all selected counties except in Lac Ste. Anne and Lethbridge where it was projected to increase concerning the historical. This might be the results of complex in-situ interaction of other hydrological variables such as precipitation, temperature, and ET. Overall, monthly SM is likely to decrease in the future at the scale of entire agricultural region, except in the RCP 8.5 scenario where the results showed an increasing SM. Monthly DA and WYLD were expected to increase in the future, except in Lac Ste. Anne and Lethbridge where monthly variations showed a mixed pattern ( Figure S5).  Figure 7 shows the long-term average monthly variation of hydrological variables. The results showed that monthly precipitation, temperature, ET, and SM had a distinct and similar pattern for current and future periods. Precipitation, temperature, and ET followed the bell-shaped curve with the maximum during the agricultural growing season (May-August). Precipitation was projected to increase during the winter and spring seasons (December to May). However, it was likely to show a mixed pattern for the summer and fall seasons (June to November). Projected changes showed a likely uniform increase in temperature throughout the year (Figure 7). Overall, ET was likely to increase and decrease in the ascending and recession limbs of the bell-shaped curve, respectively except for Lethbridge county (Figure 7). SM was projected to decrease in all selected counties except in Lac Ste. Anne and Lethbridge where it was projected to increase concerning the historical. This might be the results of complex in-situ interaction of other hydrological variables such as precipitation, temperature, and ET. Overall, monthly SM is likely to decrease in the future at the scale of entire agricultural region, except in the RCP 8.5 scenario where the results showed an increasing SM. Monthly DA and WYLD were expected to increase in the future, except in Lac Ste. Anne and Lethbridge where monthly variations showed a mixed pattern ( Figure S5).

Implications for Agriculture
We compared ET for seven selected counties under two land use conditions (ET_DLU and ET_barley) to analyze their water demand and water deficits. We considered ET as the consumptive water use of plants, and WD as water deficit in the soil that is meant to be compensated by irrigation from other sources (i.e., precipitation minus ET) allowing the plants to reach their actual yields (e.g., barley) or biomass production (e.g., DLU). Figure 8 demonstrates both ET and WD for the selected counties. In a semi-arid region, where agriculture is not the primary land use type, this ET and WD comparison creates useful information to convert the land use for a suitable agricultural crop instead of DLU; in turn, may help to feed growing population. However, we may need to consider the environmental consequences of converting land use types. Because, the change in land use type can alter the atmospheric concentration of CO 2 , the principal heat-trapping gas and local climate by changing the energy balance [48]. The WD results showed an apparent pattern for all counties with a negligible WD during May and June, and a considerable amount in July and August for some counties. The maximum negative WD was found in Lethbridge County, indicating a large soil water deficit for the growing season of plants (e.g., barley and DLU) for this county. This result was expected, since historically the Lethbridge County receives less precipitation with high temperature, and most water needs of crop are supplied from irrigation. This information about the ET and WD for DLU and barley (as an example) may help agricultural managers and planners to adopt different strategies to cope with climate change.
Water 2018, 10, x FOR PEER REVIEW 16 of 20

Implications for Agriculture
We compared ET for seven selected counties under two land use conditions (ET_DLU and ET_barley) to analyze their water demand and water deficits. We considered ET as the consumptive water use of plants, and WD as water deficit in the soil that is meant to be compensated by irrigation from other sources (i.e., precipitation minus ET) allowing the plants to reach their actual yields (e.g., barley) or biomass production (e.g., DLU). Figure 8 demonstrates both ET and WD for the selected counties. In a semi-arid region, where agriculture is not the primary land use type, this ET and WD comparison creates useful information to convert the land use for a suitable agricultural crop instead of DLU; in turn, may help to feed growing population. However, we may need to consider the environmental consequences of converting land use types. Because, the change in land use type can alter the atmospheric concentration of CO2, the principal heat-trapping gas and local climate by changing the energy balance [48]. The WD results showed an apparent pattern for all counties with a negligible WD during May and June, and a considerable amount in July and August for some counties. The maximum negative WD was found in Lethbridge County, indicating a large soil water deficit for the growing season of plants (e.g., barley and DLU) for this county. This result was expected, since historically the Lethbridge County receives less precipitation with high temperature, and most water needs of crop are supplied from irrigation. This information about the ET and WD for DLU and barley (as an example) may help agricultural managers and planners to adopt different strategies to cope with climate change. Increased food demand due to a growing population has to be met by agriculture intensification [49]. Climate change is posing a major challenge to the agricultural sector through increasing concentration of greenhouse gases [3]. In the agriculture sector, climate change may lead to a major shift and extension of croplands favoring proper environment for water demand of crop growths. Results of this study showed clear evidence of increasing temperature in the Alberta agricultural region. This rising temperature trend could increase the risk of plant diseases, insects Increased food demand due to a growing population has to be met by agriculture intensification [49]. Climate change is posing a major challenge to the agricultural sector through increasing concentration of greenhouse gases [3]. In the agriculture sector, climate change may lead to a major shift and extension of croplands favoring proper environment for water demand of crop growths.
Results of this study showed clear evidence of increasing temperature in the Alberta agricultural region. This rising temperature trend could increase the risk of plant diseases, insects infestations, and invasive weeds [50]. Warming temperature also results in longer frost-free seasons which is undoubtedly favorable for the Alberta agriculture. However, the extreme temperature could increase the likelihood of heat-wave related reduction in crop production and consequently in other agricultural industry such as beef-cattle and milk production in the dairy industries in Alberta.

•
Results revealed that the climate in the agricultural region of Alberta had become warmer and drier during the His period. The climate condition is expected to be similar in future periods. Seasonal and annual precipitation is expected to increase by 1% and 3% in the NF while they are projected to increase by 5% and 7% in the FF. The mean seasonal and annual temperature is likely to increase by 1.21 and 1.33 • C in the NF while they are expected an increase by 2.14 and 2.32 • C in the FF, respectively. ET and SM distribution in the future has a resemblance with temperature and precipitation distribution. For instance, a region with high temperature is projected to have high ET and low SM. The blue water resources (DA and WYLD) is likely to increase in the future.

•
Trend analysis showed that magnitude of increase and decrease in seasonal precipitation is higher than that of annual precipitation. Mean temperature generally has a higher trend magnitude in the southern part than the north, and a region with a low mean temperature has a higher warming rate. ET shows decreasing trends in the historical as well as in the future periods. SM does not indicate an apparent trend to conclude in the selected counties. However, it is likely to have an increasing trend for the study area. DA and WYLD show very mild trend both in the historical and future periods. • Long-term average monthly variation of precipitation is expected to increase in winter and spring seasons. The temperature is likely to increase all the year round. ET is expected to increase and decrease in the ascending and recession limbs of the bell-shaped curve having the peak in July. The SM is projected to decrease considering the entire agricultural region, while blue water resources are projected to increase in the future. • Comparison of water demand (ET) and water deficit (WD) for DLU and barley (as an example crop) indicated that there was no water deficit in May and June, while water deficit existed in July and August in some counties during the His period, that was compensated by irrigation.
It is noteworthy that the causes behind the changes of hydrological variables such as increasing atmospheric greenhouse gases due to anthropogenic actions, natural climate variability, and the effect of management and dam operations on water availability were not explicitly addressed in this study. The results of trend analysis could be influenced by the time period used for the analysis. Hence, the results need to be interpreted carefully. However, the information on projected changes will undoubtedly be useful in adaptation-related decision-making.

Conflicts of Interest:
The authors declare no conflicts of interest.