Next Article in Journal
Phosphorus Release from Sediments in a Raw Water Reservoir with Reduced Allochthonous Input
Previous Article in Journal
Retrieving Surface Soil Moisture over Wheat-Covered Areas Using Data from Sentinel-1 and Sentinel-2
Previous Article in Special Issue
Assessing Hydrological Vulnerability to Future Droughts in a Mediterranean Watershed: Combined Indices-Based and Distributed Modeling Approaches
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of SWAT Using Snow Data and Detecting Climate Change Impacts in the Mountainous Eastern Regions of Turkey

1
Department of Civil Engineering, Istanbul University-Cerrahpaşa, Istanbul 34320, Turkey
2
Department of Civil Engineering, Eskisehir Technical University, Eskisehir 26555, Turkey
*
Author to whom correspondence should be addressed.
Water 2021, 13(14), 1982; https://doi.org/10.3390/w13141982
Submission received: 21 June 2021 / Revised: 11 July 2021 / Accepted: 17 July 2021 / Published: 19 July 2021

Abstract

:
In recent years, the potential impacts of climate change on water resources and the hydrologic cycle have gained importance especially for snow-dominated mountainous basins. Within this scope, the Euphrates-Tigris Basin, a snow-fed transboundary river with several large dams, was selected to investigate the effects of changing climate on seasonal snow and runoff. In this study, two headwater basins of the Euphrates River, ranging in elevation between 1500–3500 m, were assigned and SWAT was employed as a hydrological modeling tool. Model calibration and validation were conducted in a stepwise manner for snow and runoff consecutively. For the snow routine, model parameters were adjusted using MODIS daily snow-covered area, achieving hit rates of more than 95% between MODIS and SWAT. Other model parameters were calibrated successively and later validated according to daily runoff, reaching a Nash-Sutcliffe efficiency of 0.64–0.82 in both basins. After the modeling stage, the focus was drawn to the impacts of climate change under two different climate scenarios (RCP4.5 and RCP8.5) in two 30-year projection periods (2041–2070 and 2071–2099). From the results, it is estimated that on average snow water equivalent decreases in the order of 30–39% and snow-covered days shorten by 37–43 days for the two basins until 2099. In terms of runoff, a slight reduction of at most 5% on average volume is projected but more notably, runoff center-time is expected to shift 1–2 weeks earlier by the end of the century.

1. Introduction

Mountainous regions serve as a lifeline in socio-hydrological systems, as the origin of most rivers stems from such elevations [1,2,3,4]. Snow at high altitudes, sometimes called ‘water towers’ [3,5], contributes significantly to discharge downstream, and plays a major role in the hydrological cycle. Therefore, knowledge of the amount and spatial distribution of snow cover are essential [6]. However, snow measurements are particularly difficult in such areas, and to overcome the limitations of data scarcity, the utilization of satellite imagery has been an effective method for many hydrological studies since the late 1970s [7,8,9,10,11,12,13,14,15]. Because snow-dominated areas are very susceptible to increasing temperatures, snow is also a strong indicator of climate change. According to the Intergovernmental Panel on Climate Change (IPCC) report [16], there is no doubt that the climate system has been changing due to human influence; therefore, in recent years, many researchers have focused their attention on climate change impacts in snow-dominated mountainous basins [2,17,18,19]. They commonly emphasize the need for analyzing future water availability in such regions, usually lacking sufficient ground observations.
The Soil and Water Assessment Tool (SWAT) [20] has had widely successful applications in hydrological and environmental issues over many regions of the world [21,22,23,24]. Generally, hydrologic modelling in mountainous areas is challenging; moreover, [25] reported that SWAT had inadequate performance on high-elevated catchments until a new algorithm was adopted by [26], incorporating elevation bands that discretize the snowmelt process based on basin topographic controls.
Along with this development, different studies have been applied at various parts of the world on the purpose of observing orographic effects using snow-melt parameters [27,28,29,30]. In some snow-dominated basins, runoff and snow processes of different models are compared with SWAT [31,32]. Additionally, remote sensing data like MODIS images are used along with SWAT in various mountainous regions [14,33]. In these studies, ground or satellite snow data are exploited in order to validate the model representation of the snow extent both spatially and temporally. Furthermore, SWAT is also made use of as an effective hydrologic simulator when assessing the impact of climate change on different regions of the globe [34,35,36,37]. Overall, SWAT produces acceptable results and has been shown to be an efficient tool for the planning and management of water resources, including snow-dominated catchments, when compared to other models [24,38,39,40,41].
In the current study, SWAT is utilized for the simulation of two contiguous headwater basins located in the mountainous eastern Turkey, with ArcSWAT, a GIS-based interface. There have been previous applications of SWAT in Turkey covering different fields: water quality [42,43], optimal water-management strategies [43,44], the effect of climate change on groundwater [45], and water resource potential [46,47], but the novelty of this study is its use of remotely sensed and in situ measured snow data for SWAT model calibration/validation. Hence, the novel contribution and main objectives of the current research are: (i) the application of SWAT in snow-dominated mountainous basins as a pioneer study in Turkey, (ii) utilizing remote sensing and ground snow data for model calibration and validation besides runoff, and (iii) evaluating the impacts of climate change scenarios on the mountainous headwater basins.

2. Materials and Data

2.1. Study Area

For the area of application, Murat (E21A022, 39.54° N–42.78° E) and Karasu (E21A054, 39.94° N–40.76° E) basins, both located at the headwaters of Euphrates River in the eastern part of Turkey, were selected (Figure 1). Murat and Karasu basins have areas of 5900 km2 and 2890 km2 and mean slopes of 15.7% and 16.5%, respectively. Elevations range from 1559 m to 3516 m for Murat and from 1675 m to 3156 m for Karasu with a similar hypsometric mean of around 2100 m. According to the European Environment Agency—Coordination of Information on the Environment (EEA-CORINE) [48] Land Cover dataset, the major land cover types for each basin are rangelands (55.4% Murat and 52.4% Karasu) and agricultural areas (36.7% Murat and 38.8% Karasu). As for the soil map produced by the Food and Agriculture Organization of the United Nations Educational, Scientific and Cultural Organization (FAO-UNESCO) [49], different types of kastanozems and leptosols are the major soil groups in the two basins (100% Murat and 93% Karasu).
Annual average discharges (Q) recorded at stream gauging stations of the basins are 55.03 m3/s (294.2 mm) for Murat and 22.38 m3/s (244.2 mm) for Karasu. Agri (1648 m) and Erzurum (1758 m) are two meteorological stations within the basin boundaries (Figure 1). Annual total precipitation (P) and mean temperature (T) are around 490 mm and 6.7 °C for Agri station and 470 mm and 5.4 °C for Erzurum station. Yearly hydro-meteorological conditions for both basins are presented in Figure 2.

2.2. SWAT Model

The Soil and Water Assessment Tool (SWAT) is a comprehensive, semi-distributed, continuous, process-based hydrological model developed in the early 1990s [20]. Initially, the objective in model development was to predict the impact of management on water, sediment, and agricultural chemical yields in large ungauged basins. The current version of SWAT [50] has been extensively developed since its origin based on the Simulator for Water Resources in Rural Basins (SWRRB) model [51]. A detailed history of SWAT model development can be found in [52,53,54].
In recent years, SWAT has become a well-known and convenient model for many applications, from hydrological to environmental issues, and its area of utilization is widespread throughout the world [22,23,40,41]. On the other hand, the model’s setup and calibration procedures are demanding because of several physical model inputs and parameters.
In SWAT, a basin is divided into subbasins based on topography; moreover, hydrological response units (HRUs) for each subbasin are defined as unique combinations of slope, land use, and soil. In this study, ArcSWAT GIS-based interface is used to identify subbasins and HRUs, as well as handling model input files. Major hydrological processes that include surface runoff, lateral flow, evapotranspiration, recharge, revap, return flow, infiltration, and percolation are computed in each HRU within a water-balance calculation based on four reservoirs (snow, soil, and shallow and deep aquifer) [55,56,57].
It has been proven in several studies that using elevation bands with associated snow parameters and lapse rates increases the model’s success for snow-dominated mountainous basins [6,14,26,58,59,60,61]. Among these studies, there is no common idea or recommendation on the number of elevation bands selected, which can be defined as a function of elevation or area. For this study, ten elevation bands are preferred in modelling to elaborate the model snow outputs in detail.

2.3. Model Inputs

Table 1 shows model inputs under four major classes: HRU definition, climate, calibration/validation, and climate projection data. SWAT needs some physical data in order to properly define HRUs. Shuttle Radar Topography Mission (SRTM) data [62] is utilized as a digital elevation model (DEM) in delineating the basin and river network. CORINE and FAO data are employed for generating land-use and soil type, respectively.
For the climate data, daily precipitation and temperature were obtained from the Turkish Meteorological Service (MGM) gauge measurements. The optional relative humidity, solar radiation, and wind speed data were gathered from global weather data derived from Climate Forecast System Reanalysis (CFSR) [63], which is provided as an option under SWAT.
For model calibration/validation, discharge and snow water equivalent (SWE) data were collected from the Turkish Hydraulic Works (DSI). Daily discharge at basin outlets and monthly snow tube SWE measurements from snow stations located in the vicinity of the basins are utilized. As for remote sensing data, MODerate Resolution Imaging Spectroradiometer (MODIS) daily snow-cover images (M*D10A1) were downloaded from NASA [64]. MODIS (Terra/Aqua) snow-cover extent was processed through a series of cloud-filtering algorithms (including a combination of Terra/Aqua, spatio-temporal, elevation, and seasonal filters) to obtain the daily cloud-free snow extent [65].
Climate projection data were gathered from the Turkish Meteorological Service (MGM) using the output of the MPI-ESM-MR Global Circulation Model whose climate parameters are dynamically downscaled to 20 km domain.

3. Modelling Studies and Results

3.1. Base Model Setup

For setting up the base model without any calibration, ArcSWAT 2012.10_4.19 desktop version was used as an interface. Firstly, in the watershed delineation step, the minimum subbasin threshold area (STA) values were defined credibly [66,67]. When a certain STA was selected for each basin (Table 2), SWAT automatically discretized 45 subbasins and 663 HRUs for Murat, along with 41 subbasins and 462 HRUs for Karasu basin. Ten elevation bands were preferred for each basin. Finally, the selected timeline for running the model in terms of water year was chosen as: three years (1999–2001) for warm-up, six years (2002–2007) for calibration, and four years (2008–2011) for validation.

3.2. Snow Parameters Fitting Procedure

SWAT has several model parameters for calibration and for such processes; an external automatic calibration software called SWAT-CUP (SWAT-Calibration and Uncertainty Program) [68] is preferred by SWAT-users for ease and efficiency [57]. When all the model parameters were set for calibration with runoff only in SWAT-CUP, unrealistic values were achieved, especially for some of the snow parameters, as also experienced by [61]. These abnormalities were especially seen in SMFMN (minimum melt factor in the Northern Hemisphere) parameter being greater than SMFMX (maximum melt factor in the Northern Hemisphere), as well as in high/low PLAPS/TLAPS rates. These results could be attributed to either calibrating too many parameters at the same time or fitting the model parameters to flow only. As stated by [69], according to the parametrization procedure, snow parameters and lapse rates should be fitted first before flow calibration. For this reason, a temporal methodology was undertaken in this study to compare the SWE-SCA (SWAT derived Snow Water Equivalent in mm and MODIS derived Snow-Covered Area in percent) relationship.
SWE calculations in SWAT are HRU-based; therefore, values were first transformed into subbasins and later converted as basin-wide SWE outputs. On the other hand, SCA values were derived from MODIS daily cloud-free snow extent for each basin. Although these two snow components (SWE vs SCA) have different units, a time series comparison could be carried out [12,14,70,71,72,73,74,75,76].
In our fitting procedure, the consistency between SWE and SCA was achieved by manually adjusting the SWAT snow parameters within physical limits using previous investigations on the basins [77,78,79]. Figure 3 illustrates an example of the SWE-SCA comparison for one year (2004) in each basin, including default and calibrated snow parameter sets. In this comparison, the important elements are the timing match between the beginning and end dates of snow cover as well as the accumulation and depletion phases. Figure 4 depicts SWAT-derived SWE using default/calibrated snow parameters and MODIS-derived SCA in both Murat and Karasu basins for the whole modeling period (2002–2011). A simple hit rate performance measure was utilized to assess the comparison success as given in Equation (1). In order not to degrade the match performance from small fluctuations (due to improper satellite snow-cloud discrimination or ground partial snow cover at the beginning and end periods of the snow season), a limit of 1 mm for SWE and 5% for SCA was set.
Hit rate = [(a + d)/(a + b + c + d)] × 100
where ‘a’ is the number of snow days for both SWE and SCA, ‘d’ is the number of no snow days for both SWE and SCA, ‘b’ is the number of snow days for SWE and no snow days for SCA, and ‘c’ is the number of no snow days for SWE and snow days for SCA.
Table 3 displays the SWE-SCA comparison hit-rate scores for the calibrated model in each basin and water year. Overall, hit rate scores are above 95% both in the calibration and validation periods, indicating a good fit in terms of temporal comparison. Table 4 tabulates the suitable snow parameters and lapse rates fitted after applying the SWE-SCA temporal methodology. These values were manually calibrated, and tended to suit both basins (since the two basins are quite similar in terms of climate, topography, and land use), hence they may not be global optimums but rather a physically meaningful and trustable parameter set.

3.3. Flow-Sensitivity Analysis, Calibration, Validation and Water Balance

After manually fitting the snow parameters, the calibration procedure for the remaining model parameters was conducted utilizing SWAT-CUP according to observed runoff. Before starting the calibration, to determine the critical parameters, a sensitivity analysis is performed by SWAT-CUP using the Sequential Uncertainty Fitting Algorithm—version2 (SUFI-2) algorithm, which relies on the Latin hypercube sampling [80]. A global (all-at-a-time) sensitivity analysis procedure as well as a large number of runs (around 1000 were proposed in order to identify the impact of each parameter on the objective function) was preferred to achieve more reliable results described by [23]. Table 5 tabulates the parameters considered in sensitivity analysis and their statistical results for each basin.
There are two statistics employed for the assessment of sensitive parameters in SUFI-2 algorithm, named t-stat and p-value. The larger the absolute value of t-stat, the more sensitive a parameter becomes; on the other hand, the significance of parameter sensitivity increases as p-value approaches zero. At the end of the sensitivity analysis, ten common model parameters: CN2, ALPHA_BF, SOL_Z, SOL_AWC, SOL_K, SOL_BD, GW_DELAY, RCHRG_DP, GWQMN, and REVAPMN were selected as being more sensitive at each basin and set forth for the calibration process, as shown in Table 5.
For the automatic calibration phase on streamflow, the SUFI-2 algorithm of SWAT-CUP was again utilized with the selected sensitive model parameters. Suitable change methods, relative (r_) or replace (v_) were assigned for each parameter along with the upper and lower limits. The term “r” is used for the relative adjustment of a parameter within a given range and the term “v” means directly replacing the parameter value with the assigned value. Two iterations with 1500 runs each were performed with Nash-Sutcliffe Efficiency (NSE) set as the objective function. The results of the calibration process are provided in Table 6 showing the parameter value/range for each basin.
During the SWAT-CUP calibration process, the degree of uncertainty was also evaluated using two statistical measures referred to as p-factor and r-factor. The p-factor is defined as the percentage of measure data captured by 95% prediction uncertainty (95 PPU), calculated at the 2.5% and 97.5% levels of cumulative distribution of an output variable obtained through Latin hypercube sampling, and the r-factor shows the average thickness of the 95 PPU band. The p-factor and r-factor are 0.94, 0.38 for Murat and are 0.79, 0.21 for Karasu basin, respectively.
In order to evaluate the goodness of fit, Nash-Sutcliffe Efficiency (NSE), Coefficient of Determination (R2) and Percent Bias (PBIAS) values were determined. Table 7 shows performance statistics for SWAT modeling considering daily calibration/validation results in each basin. Although there are no explicit standards, model outcomes could be assessed according to the ratings adapted by [81]. Accordingly, almost all categories are within a “good” to “very good” performance rating, confirming a trustable model setup. For Murat basin, calibration results are slightly better than the validation period, whereas this is reversed for Karasu basin. A pictorial representation of the hydrographs for each basin is shown in Figure 5. All periods considering the calibration (2002–2007) and validation (2008–2011) are seen in common charts for Murat and Karasu basins and the (dis)agreement can easily be detected between observed and simulated flow data for each year.
SWAT provided water balance components as model outputs on monthly and annual basis. These outputs were first compared annually in terms of total water yield (WYLD) between the simulated and observed values considering 2002–2011 water years for each basin (Table 8) where the error margin seems to be within 1–2%. Afterwards, each modelled water balance component was analyzed annually as shown in Table 9. Most of the input components appear to be similar between the two basins, although output constituents show important differences. Finally, Figure 6 depicts some of the principal water balance components (precipitation, water yield, and evapotranspiration) in monthly terms for Murat and Karasu basins. Since they are mountainous basins, most of the streamflow is generated in the spring months due to snowmelt, and evapotranspiration is high during summer.

3.4. Snow Validation with Ground and Satellite Data

After the model had been calibrated and validated using flow data, additional validation was carried out with different types of snow data. The first of these utilized manually measured discrete snow observations. There are six snow stations located in and around the two basins under study (Figure 1). Manual snow observations with a snow tube were conducted at these snow stations once/twice a month. The station locations, names, elevations and the corresponding subbasin numbers are shown in Figure 7. Being located at different elevations, each snow station corresponded to a certain elevation zone of the nearest subbasin. In order to test the representativeness of the snow station versus elevation band of the subbasin, Figure 8 illustrates discrete SWE measurements on the continuous model SWE for selected years.
Looking at Figure 8, it is seen that in Murat basin, Eleskirt and Haciomer stations show a good correspondence with their elevation band, whereas Dogangun station gives a higher observed SWE compared to its model. In Karasu basin, Hacimahmut station gives a fine agreement to its elevation; on the other hand, Guzelyayla and Yesildere stations seem to slightly overestimate their corresponding elevation range in the snow season. These outcomes are commonly observed in many situations when point measurements are assumed to represent a certain spatial/elevation extent. Overall, it can be deduced that the match between discretely measured and continuously modelled SWE may show differences from year to year and station to station, although the general trend agrees quite well.
The second validation assessment using snow data compared SWAT-modelled SWE with MODIS satellite SCA in spatial terms. In the calibration process, a temporal analysis was conducted, resulting in a very high correspondence (over 95% match). This time a spatial analysis was undertaken to relate modelled SWE and remotely sensed SCA during a snow season. In the literature, this relation (SWE-SCA threshold, meaning the amount of modelled SWE in mm over which there is complete snow cover by satellite) has been tested by different authors, sometimes keeping the threshold constant [14,71] and variable at others [73,76] during the modeling period. In this study, HRU-based SWE calculations of SWAT were converted into subbasins and compared with MODIS cloud-free SCA for each subbasin unit. This methodology was evaluated for an example snow season of 2006 where most of the snow conditions (SWAT SWE and MODIS SCA accumulation, peak and depletion value/timing) were on average terms (Figure 4).
Since it is not practical to visualize the results in daily steps, example periods representing snow accumulation (4 December 2005) and snow depletion (8 April 2006) were selected for each basin as depicted in Figure 9. Two different assumptions are tested here, one keeping the SWE-SCA threshold constant throughout the year, and the other varying the threshold for accumulation and depletion periods until reaching the best possible agreement between SWE and SCA grids. Figure 9 parts (a,b) and (e,f) represent the two basins in the accumulation and depletion periods, respectively, with a 15 mm SWE-SCA constant threshold as employed by [14]. On the other hand, Figure 9 parts (c,d) and (g,h) show the same intervals utilizing variable thresholds in each subbasin to have the best possible visual match between SWAT-modelled SWE and MODIS-derived SCA as proposed by [76]. The results for the daily 2006 snow season are summarized in Table 10. It is seen that even though the SWE-SCA threshold ranges are a little different for the two basins, most probably due to topographic and climatic factors, there is a clear variation between the accumulation and depletion periods. The thresholds were higher during the depletion phase as compared to the accumulation. Not only time intervals but also aspect differences were spotted for each basin. The subbasins comprising of the north mountain ridges (facing south) have lower SWE-SCA thresholds compared to the south ridges (facing north) in both basins. These results indicate that besides snow season accumulation/depletion differences, aspect variations may also need to be considered directly or indirectly while producing the HRUs within the snow routine of the model. More details on this topic can be assessed in the future, when there is a good set of ground snow data representing different elevation bands along with the availability of finer resolution satellite snow cover for a smaller-sized basin.

3.5. Climate Change Impacts

The climate is a pivotal factor that directly affects all elements of the hydrological cycle. As mentioned in the Fifth Assessment Report of the IPCC [16], changing precipitation and temperatures are altering hydrological systems, especially in mountainous regions. In the same report, four different scenarios referred to as Representative Concentration Pathways (RCPs—RCP2.6, RCP4.5, RCP6, and RCP8.5) are emphasized depending on the greenhouse gas concentration trajectories.
In the recent years, the Turkish Meteorological Office (MGM) has conducted climate projections to reveal the impacts of climate change in different regions of the country under the project “Climate Projections with New Scenarios for Turkey” [82]. Three Global Circulation Models (GCMs), HadGEM2-ES, MPI-ESM-MR, and GFDL-ESM2M, were exploited, and climate parameters were further dynamically downscaled using RegCM4.3.4 (Regional Climate Model) based on the RCP4.5 and RCP8.5 scenarios. While 1971–2000 was used as the reference period, the projection period covers 2016–2099.
In this study, the precipitation and temperature values of MPI-ESM-MR projections using the two scenarios (RCP4.5 and RCP8.5) were utilized in SWAT, as they were seen to have the best match with the reference observation datasets [83]. The projections are treated in two 30-year sub-periods: 2041–2070 (1st period) and 2071–2099 (2nd period). According to the MPI-ESM-MR projections based on two emission scenarios during the selected two periods, the model input variables showed changes as given in Table 11 with respect to the reference period (1971–2000). For the two basins under study, the variations in temperature (ΔTmp) are all in the positive direction (+1.8–+4.5 °C), increasing with the more pessimistic scenario as well as period. When the change in precipitation (ΔPcp) is considered, there are more fluctuations within periods and basins. Some of the positive changes may be attributed to the Black Sea effect (several climate change projections indicate an increase or no change in precipitation especially in the eastern Black Sea region) close to where Karasu basin is situated, as seen in Figure 1. For Murat basin, precipitation shows a small decrease almost unrelated to the type of scenario. Using these projections, climate change impacts were assessed with the calibrated SWAT model over the two basins in terms of snow and discharge.
For snow comparison, SWAT-averaged SWE outputs were plotted for reference, first and second periods of the two emission scenarios (RCP4.5 and RCP8.5) in Murat and Karasu basins as shown in Figure 10. It is clearly seen that for long-term averages, the change in mean SWE values (ΔSWE) is negative, which means less snow will accumulate on the ground compared to what currently occurs. From the figure, it is also detectable that by the end of the century, the snow season will narrow as snow starts to accumulate in mid-November instead of late October, and it will deplete quicker and disappear at the end of May rather than June. In terms of numbers, the change in the snow-covered days (ΔSnow days) is anticipated to decrease by 23 to 43 days in Murat basin and 22 to 37 days for Karasu basin in the two projection periods, respectively (Table 12). For average SWE values, a decrease of 21% and 30% is expected to occur in Murat basin during the first and second periods, respectively, whereas these values reach 26% and 39% for Karasu basin. Figure 11 depicts SWE change rates for different elevation zones in the second period (2071–2099) of the RCP8.5 scenario, when the impact is most evident. It is clearly observed that the highest variations are in the lower elevations below 1800 m.
Considering the impacts on discharge (ΔQ), the volume of runoff does not seem to be affected as much as snow, showing a decrease of at most 5% (Table 13). This is probably closely related to the change in precipitation between the periods and scenarios (Table 11) that show varying trends in each basin. For discharges, a center-time (CT) concept [17] is used which marks the day of the water year when 50% of annual total runoff is reached. There is a close relation between CT and snow melt in snow-dominated basins, because a temporal shift in the CT may indicate the impact of climate change. Table 13 shows the days representing CT for the reference time, first and second projection intervals of RCP4.5 and RCP8.5 scenarios in Murat and Karasu basins. Although the reference CTs for the two basins are almost a week apart, from the numbers it is evident that CT is shifting to earlier dates in the future (Figure 12). Since the region of interest is mountainous headwaters and the climate change scenarios project a gradual increase in temperatures, as expected more dramatical shifts are foreseen in the second period of the RCP8.5 scenario. With these results, CT for Murat is predicted to shift towards mid-April, whereas in Karasu this time shift shows the end of April. Overall, a 9- to 11-day CT translation is expected for the area of study by 2099. These results agree with the outcomes of other studies in the literature [18,19,84] stating that a 1–2 week CT translation is expected for regions above 2000 m and 2–4 week shifts are anticipated for areas between 1000–2000 m, thus predicting that no permanent snow cover will exist below 1000 m by the end of the century.

4. Conclusions

This research investigates a pioneer application of SWAT in the mountainous eastern regions of Turkey. Since the study comprises high elevated areas, snow is the most significant element of the hydrological cycle. However, it is well-known that hydrological modeling in snow-dominated basins is challenging due to the scarcity of observed data, considering complex topography and harsh climatic conditions.
In this context, SWAT was chosen to model two snow-dominated headwater basins of the transboundary Euphrates River in Turkey. A stepwise calibration procedure was applied whereby firstly, model snow parameters were manually adjusted within physical limits exploiting a satellite-assisted temporal methodology. Over 95% hit rate was achieved comparing MODIS daily cloud-free snow cover images with SWAT snow outputs for the modeling period. Such a calibration method showed a valuable example of adjusting model snow parameters using snow observations rather than runoff which could be applicable for poorly gaged areas. The second part of the calibration procedure continued by applying the SUFI-2 algorithm within SWAT-CUP software. Initially, sensitive model parameters were determined and later, they were automatically calibrated according to NSE objective function on daily runoff terms. NSE calibration/validation results ranged between 0.64–0.82 for daily timescales, corresponding to 0.74–0.89 for monthly terms, respectively. These values indicate that SWAT could produce acceptable results also in mountainous catchments.
Apart from validating the model with runoff only, validation with respect to snow was also undertaken. A satisfactory trend was observed between continuous model and discrete point snow measurements during the snow season, although some point observations deviated from their corresponding elevation bands at specific periods. This could be improved by collecting more manual points or even better by setting up automatic stations measuring continuous SWE or snow depth for a more valuable direct comparison. In terms of spatially distributed snow validation, it was displayed that instead of a single SWE-SCA threshold, varying limits were demonstrated especially for accumulation and depletion periods. Furthermore, a noteworthy aspect effect was highlighted as a detailed future study that could be assessed while producing HRUs within the model setup.
Once a trustable hydrological model was attained, the impacts of climate change on the basins were investigated. Regionally downscaled MPI-ESM-MR model clearly indicated an increasing temperature trend with different emission scenarios (RCP4.5 and RCP8.5) and projection periods (2041–2070 and 2071–2099), but a variable precipitation change. Accordingly, a more drastic variation is expected to occur in terms of snow amount (30–39%) and shortening of snow days (37–43 days) rather than runoff volume. The greatest SWE change is foreseen to occur in the lower elevations (<1800 m) of the region. Although runoff volume seems to be affected less, the shift in streamflow timing to earlier dates (1–2 weeks) is a notable indication of changing climate for the headwater basins.
The results of this study indicate that climate change will have a significant impact on snow amounts and runoff regimes for the two headwater basins which may continue to increase downstream to the large storage reservoirs by the end of the century. Multi-hydrologic models including multi-climate projections expanding to the water-scarce lower areas should be conducted in order to draw more concrete conclusions on the outcomes. Overall, the consequences of climate change promote decision-makers to consider new plans and more sustainable management techniques in terms of water resources, agriculture, livestock, and winter tourism in the region.

Author Contributions

I.B.P. performed the model simulations and drafted the first manuscript. A.A.S. performed the data and satellite image compilation, and made contributions to simulation results and editing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was partly funded by the Eskisehir Technical University Scientific Research Fund (Project No: 19ADP077).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data source details are given in “Model Inputs Section” and summarized in Table 1.

Acknowledgments

The authors would like to thank the Turkish Meteorological Service (MGM) and Turkish Hydraulic Works (DSI) for providing the data. We also thank Anonymous Reviewers for their constructive and helpful comments.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Viviroli, D.; Weingartner, R. The Significance of Mountains as Sources of the World’s Fresh Water. GAIA Ecol. Perspect. Sci. Soc. 2002, 11, 182–186. [Google Scholar] [CrossRef]
  2. Adam, J.C.; Hamlet, A.F.; Lettenmaier, D.P. Implications of global climate change for snowmelt hydrology in the twenty-first century. Hydrol. Process. 2009, 23, 962–972. [Google Scholar] [CrossRef]
  3. Viviroli, D.; Dürr, H.H.; Messerli, B.; Meybeck, M.; Weingartner, R. Mountains of the world, water towers for humanity: Typology, mapping, and global significance. Water Resour. Res. 2007, 43. [Google Scholar] [CrossRef] [Green Version]
  4. Milano, M.; Reynard, E.; Köplin, N.; Weingartner, R. Climatic and anthropogenic changes in Western Switzerland: Impacts on water stress. Sci. Total Environ. 2015, 536, 12–24. [Google Scholar] [CrossRef] [PubMed]
  5. Immerzeel, W.W.; Lutz, A.F.; Andrade, M.; Bahl, A.; Biemans, H.; Bolch, T.; Hyde, S.; Brumby, S.; Davies, B.; Elmore, A.C.; et al. Importance and vulnerability of the world’s water towers. Nat. Cell Biol. 2020, 577, 364–369. [Google Scholar] [CrossRef]
  6. Pradhanang, S.M.; Anandhi, A.; Mukundan, R.; Zion, M.S.; Pierson, D.C.; Schneiderman, E.M.; Matonse, A.; Frei, A. Application of SWAT model to assess snowpack development and streamflow in the Cannonsville watershed, New York, USA. Hydrol. Process. 2011, 25, 3268–3277. [Google Scholar] [CrossRef]
  7. Butt, M.J.; Bilal, M. Application of snowmelt runoff model for water resource management. Hydrol. Process. 2011, 25, 3735–3747. [Google Scholar] [CrossRef]
  8. Seidel, K.; Brusch, W.; Steinmeier, C. Experiences from real time runoff forecasts by snow cover remote sensing. IEEE Int. Geosci. Remote Sens. Symp. 2002, 4, 2090–2093. [Google Scholar] [CrossRef]
  9. Rango, A. Spaceborne remote sensing for snow hydrology applications. Hydrol. Sci. J. 1996, 41, 477–494. [Google Scholar] [CrossRef]
  10. Grayson, R.B.; Bloschl, G.; Western, A.W.; McMahon, T.A. Advances in the use of observed spatial patterns of catchment hydrological response. Adv. Water Resour. 2002, 25, 1313–1334. [Google Scholar] [CrossRef]
  11. Parajka, J.; Bloschl, G. Spatio-temporal combination of MODIS images—Potential for snow cover mapping. Water Resour. Res. 2008, 44, 03406. [Google Scholar] [CrossRef]
  12. Şorman, A.A.; Sensoy, A.; Tekeli, A.E.; Sorman, A.U.; Akyürek, Z. Modelling and forecasting snowmelt runoff process using the HBV model in the eastern part of Turkey. Hydrol. Process. 2009, 23, 1031–1040. [Google Scholar] [CrossRef]
  13. Hancock, S.; Baxter, R.; Evans, J.; Huntley, B. Evaluating global snow water equivalent products for testing land surface models. Remote Sens. Environ. 2013, 128, 107–117. [Google Scholar] [CrossRef] [Green Version]
  14. Grusson, Y.; Sun, X.; Gascoin, S.; Sauvage, S.; Raghavan, S.; Anctil, F.; Sáchez-Pérez, J.-M. Assessing the capability of the SWAT model to simulate snow, snow melt and streamflow dynamics over an alpine watershed. J. Hydrol. 2015, 531, 574–588. [Google Scholar] [CrossRef]
  15. Yılmaz, Y.A.; Aalstad, K.; Sen, O.L. Multiple Remotely Sensed Lines of Evidence for a Depleting Seasonal Snowpack in the Near East. Remote Sens. 2019, 11, 483. [Google Scholar] [CrossRef] [Green Version]
  16. IPCC. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Pachauri, R.K., Meyer, L.A., Eds.; IPCC: Geneva, Switzerland, 2014; p. 151. [Google Scholar]
  17. Stewart, I.T.; Cayan, D.R.; Dettinger, M.D. Changes toward Earlier Streamflow Timing across Western North America. J. Clim. 2005, 18, 1136–1155. [Google Scholar] [CrossRef]
  18. Bozkurt, D.; Sen, O.L. Climate change impacts in the Euphrates–Tigris Basin based on different model and scenario simulations. J. Hydrol. 2013, 480, 149–161. [Google Scholar] [CrossRef]
  19. Yucel, I.; Güventürk, A.; Sen, O.L. Climate change impacts on snowmelt runoff for mountainous transboundary basins in eastern Turkey. Int. J. Clim. 2014, 35, 215–228. [Google Scholar] [CrossRef]
  20. Arnold, J.G.; Srinivasan, R.; Muttiah, R.S.; Williams, J.R. Large area hydrologic modeling and assessment part I: Model development. J. Am. Water Resour. Assoc. 1998, 34, 73–89. [Google Scholar] [CrossRef]
  21. Abbaspour, K.C.; Rouholahnejad, E.; Vaghefi, S.; Srinivasan, R.; Yang, H.; Kløve, B. A continental-scale hydrology and water quality model for Europe: Calibration and uncertainty of a high-resolution large-scale SWAT model. J. Hydrol. 2015, 524, 733–752. [Google Scholar] [CrossRef] [Green Version]
  22. Gassman, P.W.; Reyes, M.R.; Green, C.H.; Arnold, J.G. The Soil and Water Assessment Tool: Historical Development, Applications, and Future Research Directions. Trans. ASABE 2007, 50, 1211–1250. [Google Scholar] [CrossRef] [Green Version]
  23. Gassman, P.W.; Sadeghi, A.M.; Srinivasan, R. Applications of the SWAT Model Special Section: Overview and Insights. J. Environ. Qual. 2014, 43, 1–8. [Google Scholar] [CrossRef] [PubMed]
  24. Tan, M.L.; Gassman, P.W.; Yang, X.; Haywood, J. A review of SWAT applications, performance and future needs for simulation of hydro-climatic extremes. Adv. Water Resour. 2020, 143, 103662. [Google Scholar] [CrossRef]
  25. Arnold, J.G.; Srinivasan, R.; Muttiah, R.S.; Allen, P.M. Continental scale simulation of the hydrologic balance. J. Am. Water Resour. Assoc. 1999, 35, 1037–1051. [Google Scholar] [CrossRef]
  26. Fontaine, T.; Cruickshank, T.; Arnold, J.; Hotchkiss, R. Development of a snowfall—Snowmelt routine for mountainous terrain for the soil water assessment tool (SWAT). J. Hydrol. 2002, 262, 209–223. [Google Scholar] [CrossRef]
  27. Wang, X.; Melesse, A. Evaluation of the swat model’s snowmelt hydrology in a northwestern minnesota watershed. Trans. ASAE 2005, 48, 1359–1376. [Google Scholar] [CrossRef]
  28. Ahl, R.S.; Woods, S.W.; Zuuring, H.R. Hydrologic Calibration and Validation of SWAT in a Snow-Dominated Rocky Mountain Watershed, Montana, USA. J. Am. Water Resour. Assoc. 2008, 44, 1411–1430. [Google Scholar] [CrossRef]
  29. Kang, K.; Lee, J.H. Hydrologic modelling of the effect of snowmelt and temperature on a mountainous watershed. J. Earth Syst. Sci. 2014, 123, 705–713. [Google Scholar] [CrossRef] [Green Version]
  30. Tuo, Y.; Marcolini, G.; Disse, M.; Chiogna, G. Calibration of snow parameters in SWAT: Comparison of three approaches in the Upper Adige River basin (Italy). Hydrol. Sci. J. 2018, 63, 657–678. [Google Scholar] [CrossRef]
  31. Troin, M.; Caya, D. Evaluating the SWAT’s snow hydrology over a Northern Quebec watershed. Hydrol. Process. 2013, 28, 1858–1873. [Google Scholar] [CrossRef]
  32. Dhami, B.; Himanshu, S.K.; Pandey, A.; Gautam, A.K. Evaluation of the SWAT model for water balance study of a mountainous snowfed river basin of Nepal. Environ. Earth Sci. 2018, 77, 21. [Google Scholar] [CrossRef]
  33. Stehr, A.; Debels, P.; Arumí, J.L.; Romero, F.; Alcayaga, H. Combining the Soil and Water Assessment Tool (SWAT) and MODIS imagery to estimate monthly flows in a data-scarce Chilean Andean basin. Hydrol. Sci. J. 2009, 54, 1053–1067. [Google Scholar] [CrossRef]
  34. Abbaspour, K.C.; Faramarzi, M.; Ghasemi, S.S.; Yang, H. Assessing the impact of climate change on water resources in Iran. Water Resour. Res. 2009, 45, 10434. [Google Scholar] [CrossRef] [Green Version]
  35. Grusson, Y.; Anctil, F.; Sauvage, S.; Pérez, J.M.S. Coevolution of Hydrological Cycle Components under Climate Change: The Case of the Garonne River in France. Water 2018, 10, 1870. [Google Scholar] [CrossRef] [Green Version]
  36. Čerkasova, N.; Ertürk, A.; Zemlys, P.; Denisov, V.; Umgiesser, G. Curonian Lagoon drainage basin modelling and assessment of climate change impact. Oceanologia 2016, 58, 90–102. [Google Scholar] [CrossRef] [Green Version]
  37. Saddique, N.; Usman, M.; Bernhofer, C. Simulating the Impact of Climate Change on the Hydrological Regimes of a Sparsely Gauged Mountainous Basin, Northern Pakistan. Water 2019, 11, 2141. [Google Scholar] [CrossRef] [Green Version]
  38. Borah, D.K.; Bera, M. Watershed-scale hydrologic and nonpoint-source pollution models: Review of mathematical bases. Trans. ASAE 2003, 46, 1553–1566. [Google Scholar] [CrossRef] [Green Version]
  39. Golmohammadi, G.; Prasher, S.; Madani, A.; Rudra, R. Evaluating Three Hydrological Distributed Watershed Models: MIKE-SHE, APEX, SWAT. Hydrology 2014, 1, 20–39. [Google Scholar] [CrossRef] [Green Version]
  40. Krysanova, V.; Arnold, J.G. Advances in ecohydrological modelling with SWAT—A review. Hydrol. Sci. J. 2008, 53, 939–947. [Google Scholar] [CrossRef]
  41. Krysanova, V.; White, M. Advances in water resources assessment with SWAT—An overview. Hydrol. Sci. J. 2015, 60, 1–13. [Google Scholar] [CrossRef] [Green Version]
  42. Ekdal, A.; Gurel, M.; Guzel, C.; Erturk, A.; Tanik, A.; Gonenc, I.E. Application of WASP and SWAT models for a mediterranean coastal lagoon with limited seawater exchange. J. Coast. Res. 2011, S1, 1023–1027. [Google Scholar]
  43. Özcan, Z.; Kentel, E.; Alp, E. Evaluation of the best management practices in a semi-arid region with high agricultural activity. Agric. Water Manag. 2017, 194, 160–171. [Google Scholar] [CrossRef]
  44. Gungor, O.; Göncü, S. Application of the soil and water assessment tool model on the Lower Porsuk Stream Watershed. Hydrol. Process. 2013, 27, 453–466. [Google Scholar] [CrossRef]
  45. Ertürk, A.; Ekdal, A.; Gürel, M.; Karakaya, N.; Guzel, C.; Gönenç, E. Evaluating the impact of climate change on groundwater resources in a small Mediterranean watershed. Sci. Total Environ. 2014, 499, 437–447. [Google Scholar] [CrossRef]
  46. Cuceloglu, G.; Abbaspour, K.C.; Ozturk, I. Assessing the Water-Resources Potential of Istanbul by Using a Soil and Water Assessment Tool (SWAT) Hydrological Model. Water 2017, 9, 814. [Google Scholar] [CrossRef] [Green Version]
  47. Cuceloglu, G.; Ozturk, I. Assessing the Impact of CFSR and Local Climate Datasets on Hydrological Modeling Performance in the Mountainous Black Sea Catchment. Water 2019, 11, 2277. [Google Scholar] [CrossRef] [Green Version]
  48. CORINE Land Cover. Copernicus Land Monitoring Service; CORINE Land Cover: Copenhagen, Denmark, 2000; Available online: https://land.copernicus.eu/pan-european/corine-land-cover (accessed on 3 May 2021).
  49. Food and Agriculture Organization of the United Nations. FAO Digital Soil Map of the World (DSMW). Available online: http://www.fao.org/geonetwork/srv/ (accessed on 24 October 2020).
  50. Bieger, K.; Arnold, J.G.; Rathjens, H.; White, M.J.; Bosch, D.D.; Allen, P.M.; Volk, M.; Srinivasan, R. Introduction to SWAT+, A Completely Restructured Version of the Soil and Water Assessment Tool. J. Am. Water Resour. Assoc. 2016, 53, 115–130. [Google Scholar] [CrossRef]
  51. Arnold, J.G.; Williams, J.R. Validation of SWRRB—Simulator for Water Resources in Rural Basins. J. Water Resour. Plan. Manag. 1987, 113, 243–256. [Google Scholar] [CrossRef]
  52. Williams, J.R.; Arnold, J.G.; Kiniry, J.R.; Gassman, P.W.; Green, C.H. History of model development at Temple, Texas. Hydrol. Sci. J. 2008, 53, 948–960. [Google Scholar] [CrossRef] [Green Version]
  53. Douglas-Mankin, K.R.; Srinivasan, R.; Arnold, J.G. Soil and Water Assessment Tool (SWAT) Model: Current Developments and Applications. Trans. ASABE 2010, 53, 1423–1431. [Google Scholar] [CrossRef]
  54. Tuppad, P.; Douglas-Mankin, K.R.; Lee, T.; Srinivasan, R.; Arnold, J.G. Soil and Water Assessment Tool (SWAT) Hydrologic/Water Quality Model: Extended Capability and Wider Adoption. Trans. ASABE 2011, 54, 1677–1684. [Google Scholar] [CrossRef]
  55. Arnold, J.G.; Fohrer, N. SWAT2000: Current capabilities and research opportunities in applied watershed modelling. Hydrol. Process. 2005, 19, 563–572. [Google Scholar] [CrossRef]
  56. Neitsch, S.L.; Arnold, J.G.; Kiniry, J.R.; Williams, J.R. Soil and Water Assessment Tool Theoretical Documentation Version 2009; Texas Water Resources Institute: Forney, TX, USA, 2011; pp. 1–618. [Google Scholar]
  57. Arnold, J.G.; Moriasi, D.N.; Gassman, P.W.; Abbaspour, K.C.; White, M.J.; Srinivasan, R.; Santhi, C.; Harmel, R.D.; van Griensven, A.; Van Liew, M.W.; et al. SWAT: Model Use, Calibration, and Validation. Trans. ASABE 2012, 55, 1491–1508. [Google Scholar] [CrossRef]
  58. Zhang, X.; Srinivasan, R.; Debele, B.; Hao, F. Runoff Simulation of the Headwaters of the Yellow River Using The SWAT Model With Three Snowmelt Algorithms1. J. Am. Water Resour. Assoc. 2008, 44, 48–61. [Google Scholar] [CrossRef]
  59. Luo, Y.; Arnold, J.; Allen, P.; Chen, X. Baseflow simulation using SWAT model in an inland river basin in Tianshan Mountains, Northwest China. Hydrol. Earth Syst. Sci. 2012, 16, 1259–1267. [Google Scholar] [CrossRef] [Green Version]
  60. Rahman, K.; Maringanti, C.; Beniston, M.; Widmer, F.; Abbaspour, K.; Lehmann, A. Streamflow Modeling in a Highly Managed Mountainous Glacier Watershed Using SWAT: The Upper Rhone River Watershed Case in Switzerland. Water Resour. Manag. 2012, 27, 323–339. [Google Scholar] [CrossRef] [Green Version]
  61. Tuo, Y.; Marcolini, G.; Disse, M.; Chiogna, G. A multi-objective approach to improve SWAT model calibration in alpine catchments. J. Hydrol. 2018, 559, 347–360. [Google Scholar] [CrossRef]
  62. Shuttle Radar Topography Mission (SRTM). SRTM 90m DEM Digital Elevation Database. Available online: https://srtm.csi.cgiar.org/ (accessed on 21 June 2021).
  63. Fuka, D.R.; Walter, M.T.; MacAlister, C.; DeGaetano, A.T.; Steenhuis, T.S.; Easton, Z.M. Using the Climate Forecast System Reanalysis as weather input data for watershed models. Hydrol. Process. 2014, 28, 5613–5623. [Google Scholar] [CrossRef]
  64. MODerate Resolution Imaging Spectroradiometer (MODIS). Available online: https://modis.gsfc.nasa.gov/ (accessed on 21 June 2021).
  65. Sorman, A.A.; Yamankurt, E. Modified satellite products on snow covered area in upper Euphrates basin, Turkey. Geophys. Res. Abstr. 2011, 13, EGU2011-7887. [Google Scholar]
  66. Jha, M.; Gassman, P.W.; Secchi, S.; Gu, R.; Arnold, J. Effect of watershed subdivision on SWAT flow, sediment, and nutrient predictions. J. Am. Water Resour. Assoc. 2004, 40, 811–825. [Google Scholar] [CrossRef] [Green Version]
  67. Winchell, M.; Srinivasan, R.; Di Luzio, M.; Arnold, J. ArcSWAT Interface for SWAT2012: User’s Guide; Blackland Research and Extension Center, Texas AgriLife Research, and USDA-ARS Grassland, Soil, and Water Research Laboratory: Temple, TX, USA, 2013. [Google Scholar]
  68. Abbaspour, K.C.; Vejdani, M.; Haghighat, S. SWAT-CUP Calibration and Uncertainty Programs for SWAT. In Proceedings of the Fourth International SWAT Conference 2007, Delft, The Netherlands, 4–6 July 2007; pp. 1596–1602. [Google Scholar]
  69. Abbaspour, K.C.; Vaghefi, S.A.; Srinivasan, R. A Guideline for Successful Calibration and Uncertainty Analysis for Soil and Water Assessment: A Review of Papers from the 2016 International SWAT Conference. Water 2017, 10, 6. [Google Scholar] [CrossRef] [Green Version]
  70. Dressler, K.A.; Leavesley, G.H.; Bales, R.C.; Fassnacht, S.R. Evaluation of gridded snow water equivalent and satellite snow cover products for mountain basins in a hydrologic model. Hydrol. Process. 2006, 20, 673–688. [Google Scholar] [CrossRef]
  71. Shamir, E.; Georgakakos, K.P. Estimating snow depletion curves for American River basins using distributed snow modeling. J. Hydrol. 2007, 334, 162–173. [Google Scholar] [CrossRef]
  72. Molotch, N.P.; Margulis, S.A. Estimating the distribution of snow water equivalent using remotely sensed snow cover data and a spatially distributed snowmelt model: A multi-resolution, multi-sensor comparison. Adv. Water Resour. 2008, 31, 1503–1514. [Google Scholar] [CrossRef]
  73. Franz, K.; Karsten, L.R. Calibration of a distributed snow model using MODIS snow covered area data. J. Hydrol. 2013, 494, 160–175. [Google Scholar] [CrossRef]
  74. Skaugen, T.; Randen, F. Modeling the spatial distribution of snow water equivalent, taking into account changes in snow-covered area. Ann. Glaciol. 2013, 54, 305–313. [Google Scholar] [CrossRef] [Green Version]
  75. Rittger, K.; Bair, E.H.; Kahl, A.; Dozier, J. Spatial estimates of snow water equivalent from reconstruction. Adv. Water Resour. 2016, 94, 345–363. [Google Scholar] [CrossRef] [Green Version]
  76. Riboust, P.; Thirel, G.; Le Moine, N.; Ribstein, P. Revisiting a Simple Degree-Day Model for Integrating Satellite Data: Implementation of Swe-Sca Hystereses. J. Hydrol. Hydromech. 2019, 67, 70–81. [Google Scholar] [CrossRef] [Green Version]
  77. Uysal, G.; Sensoy, A.; Şorman, A.A. Improving daily streamflow forecasts in mountainous Upper Euphrates basin by multi-layer perceptron model with satellite snow products. J. Hydrol. 2016, 543, 630–650. [Google Scholar] [CrossRef]
  78. Sorman, A.A.; Sensoy, A.; Coskun, C. Use of different satellite snow products in short term hydrological forecasting. In Proceedings of the EU-METSAT Meteorological Satellite Conference 2017, Rome, Italy, 2–6 October 2017. [Google Scholar]
  79. Şorman, A.A.; Uysal, G.; Şensoy, A. Probabilistic Snow Cover and Ensemble Streamflow Estimations in the Upper Euphrates Basin. J. Hydrol. Hydromech. 2019, 67, 82–92. [Google Scholar] [CrossRef] [Green Version]
  80. McKay, M.D.; Beckman, R.J.; Conover, W.J. Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output from a Computer Code. Technometrics 1979, 21, 239–245. [Google Scholar] [CrossRef]
  81. Moriasi, D.N.; Arnold, J.G.; Van Liew, M.W.; Bingner, R.L.; Harmel, R.D.; Veith, T.L. Model Evaluation Guidelines for Systematic Quantification of Accuracy in Watershed Simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar] [CrossRef]
  82. Demircan, M.; Demir, Ö.; Atay, H.; Yazici, B.; Eskioğlu, O.; Tuvan, A.; Akçakaya, A. Climate change projections for Turkey with new scenarios. In Proceedings of the Climate Change and Climate Dynamics Conference—CCCD 2014, Istanbul, Turkey, 8–10 September 2014. [Google Scholar]
  83. Demircan, M.; Gürkan, H.; Eskioğlu, O.; Arabacı, H.; Coşkun, M. Climate Change Projections for Turkey: Three Models and Two Scenarios. Turk. J. Water Sci. Manag. 2017, 1, 22–43. [Google Scholar] [CrossRef]
  84. Önol, B.; Bozkurt, D.; Turuncoglu, U.U.; Şen, Ö.L.; Dalfes, H.N.; Dalfes, N. Evaluation of the twenty-first century RCM simulations driven by multiple GCMs over the Eastern Mediterranean–Black Sea region. Clim. Dyn. 2013, 42, 1949–1965. [Google Scholar] [CrossRef]
Figure 1. Geographic location of Murat and Karasu basins.
Figure 1. Geographic location of Murat and Karasu basins.
Water 13 01982 g001
Figure 2. Hydro-meteorological conditions of (a) Murat and (b) Karasu basin.
Figure 2. Hydro-meteorological conditions of (a) Murat and (b) Karasu basin.
Water 13 01982 g002
Figure 3. Critical parts of SWE-SCA comparison chart in an example year for (a) Murat and (b) Karasu basin.
Figure 3. Critical parts of SWE-SCA comparison chart in an example year for (a) Murat and (b) Karasu basin.
Water 13 01982 g003
Figure 4. SWE (default/calibrated)-SCA comparison chart for (a) Murat and (b) Karasu basin in the modeling period.
Figure 4. SWE (default/calibrated)-SCA comparison chart for (a) Murat and (b) Karasu basin in the modeling period.
Water 13 01982 g004
Figure 5. Hydrographs for (a) Murat and (b) Karasu basin in daily time-step.
Figure 5. Hydrographs for (a) Murat and (b) Karasu basin in daily time-step.
Water 13 01982 g005
Figure 6. Average monthly values of precipitation, water yield and evapotranspiration for (a) Murat and (b) Karasu basin.
Figure 6. Average monthly values of precipitation, water yield and evapotranspiration for (a) Murat and (b) Karasu basin.
Water 13 01982 g006
Figure 7. Location of snow stations with corresponding subbasins.
Figure 7. Location of snow stations with corresponding subbasins.
Water 13 01982 g007
Figure 8. SWE comparison charts between SWAT elevation bands and snow stations for Dogangun in (a) 2002 and (b) 2003, Eleskirt in (c) 2007 and (d) 2010, Haciomer in (e) 2003 and (f) 2007, Guzelyayla in (g) 2002 and (h) 2003, Hacimahmut in (i) 2005 and (j) 2006, Yesildere in (k) 2002 and (l) 2003 hydrological year.
Figure 8. SWE comparison charts between SWAT elevation bands and snow stations for Dogangun in (a) 2002 and (b) 2003, Eleskirt in (c) 2007 and (d) 2010, Haciomer in (e) 2003 and (f) 2007, Guzelyayla in (g) 2002 and (h) 2003, Hacimahmut in (i) 2005 and (j) 2006, Yesildere in (k) 2002 and (l) 2003 hydrological year.
Water 13 01982 g008aWater 13 01982 g008b
Figure 9. SWAT-MODIS overlapped snow extent visualizations for Murat and Karasu basins using a constant 15 mm SWE-SCA threshold on (a,b) 4 December 2005, (e,f) 8 April 2006, and using variable SWE-SCA values on (c,d) 4 December 2005, (g,h) 8 April 2006.
Figure 9. SWAT-MODIS overlapped snow extent visualizations for Murat and Karasu basins using a constant 15 mm SWE-SCA threshold on (a,b) 4 December 2005, (e,f) 8 April 2006, and using variable SWE-SCA values on (c,d) 4 December 2005, (g,h) 8 April 2006.
Water 13 01982 g009
Figure 10. Average SWE values for reference, 1st (2041–2070) and 2nd (2071–2099) periods in RCP4.5 and RCP8.5 emission scenarios for (a) Murat and (b) Karasu basins.
Figure 10. Average SWE values for reference, 1st (2041–2070) and 2nd (2071–2099) periods in RCP4.5 and RCP8.5 emission scenarios for (a) Murat and (b) Karasu basins.
Water 13 01982 g010
Figure 11. Percent SWE change rates in different elevation zones for 2071–2099 (2nd period) of RCP8.5 emission scenario.
Figure 11. Percent SWE change rates in different elevation zones for 2071–2099 (2nd period) of RCP8.5 emission scenario.
Water 13 01982 g011
Figure 12. Center-time changes for (a) Murat and (b) Karasu basin.
Figure 12. Center-time changes for (a) Murat and (b) Karasu basin.
Water 13 01982 g012
Table 1. Data description, source, and scale/resolution used in SWAT.
Table 1. Data description, source, and scale/resolution used in SWAT.
Data TypeData SourceScale/Resolution
HRU Definition Data
Digital Elevation Model (DEM)Shuttle Radar Topography Mission (SRTM) Grid cell 90 × 90 m
Land useEuropean Environment Agency CORINE Land Cover (year 2000) Grid cell 100 × 100 m
SoilFAO-UNESCO Global Soil MapScale 1:5,000,000
Climate Data
Precipitation, max./min. temperatureTurkish Meteorological Service (MGM)Ground station
Rel. hum., solar rad., wind speedClimate Forecast System Reanalysis (CFSR)Grid cell ~38 km
Calibration/Validation Data
DischargeTurkish Hydraulic Works (DSI) Ground station
Snow Water Equivalent (SWE)Turkish Hydraulic Works (DSI) Ground station
Snow Cover Area (SCA)MODerate Resolution Imaging Spectroradiometer (MODIS) Grid cell 500 × 500 m
Climate Projection Data
Temperature/precipitationTurkish Meteorological Service (MGM)Grid cell 20 × 20 km
Table 2. Base model setup descriptions for each basin.
Table 2. Base model setup descriptions for each basin.
DescriptionMuratKarasu
Minimum subbasin threshold area10,000 Ha5000 Ha
HRU threshold (slope/land use/soil)0/0/0
Subbasin number4541
Slope classes0–12/12–25/>25 (%)
HRU number663462
Elevation band number10
Warm-up period1999–2001 (3 years)
Calibration period2002–2007 (6 years)
Validation period2008–2011 (4 years)
Table 3. Hit rates for SWE-SCA comparison charts for calibrated snow parameters.
Table 3. Hit rates for SWE-SCA comparison charts for calibrated snow parameters.
Hydrological
Year
Murat BasinKarasu Basin
abcdHit Rate
(%)
abcdHit Rate
(%)
20022137514096.7120016214795.07
20031798916995.3418423314690.41
20042095215098.362090814998.08
200517971016995.341836916795.89
200617102616892.882030315999.18
20071932516598.081930916397.53
Total Cal.
Period
1144295796196.071172246493195.98
200815721918894.52180112315290.96
20091775617796.991782817797.26
20101986315897.531944815996.71
201115911119496.71179241115190.41
Total Val.
Period
691143971796.37731415063993.77
Total18354396167896.25189375114157095.15
Table 4. Fitted snow parameters and lapse rates.
Table 4. Fitted snow parameters and lapse rates.
ParameterFitting Value
SFTMP1 °C
SMTMP0.5 °C
SMFMX2.5 mm H2O/°C-day
SMFMN0.5 mm H2O/°C-day
TIMP1
SNOCOVMX55 mm H2O
SNO50COV0.55
PLAPS175 mm H2O /km
TLAPS−5.5 °C/km
Table 5. Statistical results of global sensitivity analysis.
Table 5. Statistical results of global sensitivity analysis.
ParametersDescriptionMurat BasinKarasu Basin
t-Statp-Valuet-Statp-Value
CN2SCS runoff curve number2.480.01−11.730.00
ALPHA_BFBase flow alpha factor35.050.006.180.00
GW_DELAYGroundwater delay−6.880.004.570.00
GWQMNThreshold water depth in the shallow aquifer for return flow to occur−1.290.192.180.02
REVAPMNThreshold water depth in the shallow aquifer for ‘revap’ to occur2.120.031.830.06
RCHRG_DPDeep aquifer percolation fraction1.410.159.810.00
GW_REVAPGroundwater ‘revap’ coefficient−0.430.661.060.28
ESCOSoil evaporation compensation factor−0.110.90−0.680.49
EPCOPlant evaporation compensation factor−0.570.56−0.430.66
CH_K2Effective hydraulic conductivity in main channel alluvium0.650.510.450.65
CH_N2Manning’s “n” value for the main channel−0.190.841.010.31
CANMXMaximum canopy storage−0.710.470.170.86
LAT_TTIMELateral flow travel time−1.010.300.910.36
SLSOILSlope length for lateral subsurface flow0.350.720.320.74
SURLAGSurface runoff lag coefficient0.670.500.520.04
FFCBInitial soil water storage0.060.941.360.17
SOL_ZDepth from soil to bottom of layer−1.660.0927.540.00
SOL_AWCAvailable water capacity of the soil layer−1.280.206.270.00
SOL_BDSoil moist bulk density4.430.002.090.12
SOL_KSaturated hydraulic conductivity−11.950.004.620.00
Table 6. Calibrated parameters with change methods and initial ranges.
Table 6. Calibrated parameters with change methods and initial ranges.
ParametersChange MethodDefault RangeInitial RangeCalibrated Value
Lower LimitUpper LimitLower LimitUpper LimitMuratKarasu
CN2r_3598−0.30.3−0.28−0.16
ALPHA_BFv_01010.730.75
SOL_Zr_03000−0.30.3−0.210.30
SOL_AWCr_01−0.30.30.290.19
SOL_Kr_02000−0.30.3−0.190.30
SOL_BDr_0.92.5−0.30.30.100.12
GW_DELAYv_0500510010.697.53
RCHRG_DPv_010.20.50.360.40
GWQMNv_0500051000172.85155.04
REVAPMNv_010005001000770.45750.10
Table 7. Performance statistics of SWAT modeling (daily time-step).
Table 7. Performance statistics of SWAT modeling (daily time-step).
Objective FunctionMurat BasinKarasu Basin
CalibrationValidationCalibrationValidation
NSE0.730.670.640.82
R20.740.760.630.82
PBIAS−4.3−14.00 8.5−3.00
Table 8. Annual total water yield of SWAT.
Table 8. Annual total water yield of SWAT.
Murat BasinKarasu Basin
Simulated WYLD (mm)304.2246.3
Observed WYLD (mm)294.2244.2
Table 9. Annual water balance components of SWAT.
Table 9. Annual water balance components of SWAT.
Water Balance ComponentVolume (mm)
Murat BasinKarasu Basin
Precipitation508.1529.2
Snowfall183.2173.7
Snowmelt170.1168.3
Sublimation17.010.0
Surface runoff11.92.26
Lateral flow5.15.6
Ground water flow (Shallow aquifer)216.2122.6
Ground water flow (Deep aquifer)70.9115.5
Revap18.214.5
Recharge from deep aquifer70.9115.5
Recharge from total aquifer294.1238.2
Total water yield304.2246.3
Percolation294.1238.1
Evapotranspiration274.8290.9
Potential evapotranspiration925.7737.4
Table 10. SWE-SCA threshold ranges of subbasin groups for each basin.
Table 10. SWE-SCA threshold ranges of subbasin groups for each basin.
PeriodSWE Threshold Range (mm)
Murat Subbasin GroupsKarasu Subbasin Groups
North SideSouth SideNorth SideSouth Side
Accumulation Period (Oct-Dec)7–8 25–3020–2535–40
100% Snow Cover (Jan-Feb)20–25
Recession Period (Mar-May)10–1550–6025–3045–50
Table 11. Climate change data (temperature and precipitation) for each basin.
Table 11. Climate change data (temperature and precipitation) for each basin.
MuratKarasu
ΔTmp
(°C)
ΔPcp
(%)
ΔTmp
(°C)
ΔPcp
(%)
1st PeriodRCP4.5+1.86−1.10+1.79+2.84
RCP8.5+2.74−4.40+2.67−1.82
2nd PeriodRCP4.5+2.33−3.70+2.29−4.40
RCP8.5+4.43−1.46+4.34+1.31
Table 12. Detecting climate change impacts for snow outputs in different scenarios and periods.
Table 12. Detecting climate change impacts for snow outputs in different scenarios and periods.
MuratKarasu
ΔSWE
(%)
ΔSnow Days
(day)
ΔSWE
(%)
ΔSnow Days
(day)
1st PeriodRCP4.5−18.70−19−19.64−18
RCP8.5−20.95−23−25.78−22
2nd PeriodRCP4.5−25.56−22−31.82−22
RCP8.5−29.84−43−38.83−37
Table 13. Detecting climate change impacts for runoff outputs in different scenarios and periods.
Table 13. Detecting climate change impacts for runoff outputs in different scenarios and periods.
MuratKarasu
Reference CT: 27-AprilReference CT: 5-May
ΔQ (%)CT/ShiftΔQ (%)CT/Shift
1st PeriodRCP4.5−0.5722-April/5 days+3.8304-May/1 day
RCP8.5−4.2319-April/8 days−1.528-April/7 days
2nd PeriodRCP4.5−2.7323-April/4 days−2.6502-May/3 days
RCP8.5−1.1816-April/11 days+1.5626-April/9 days
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peker, I.B.; Sorman, A.A. Application of SWAT Using Snow Data and Detecting Climate Change Impacts in the Mountainous Eastern Regions of Turkey. Water 2021, 13, 1982. https://doi.org/10.3390/w13141982

AMA Style

Peker IB, Sorman AA. Application of SWAT Using Snow Data and Detecting Climate Change Impacts in the Mountainous Eastern Regions of Turkey. Water. 2021; 13(14):1982. https://doi.org/10.3390/w13141982

Chicago/Turabian Style

Peker, Ismail Bilal, and Ali Arda Sorman. 2021. "Application of SWAT Using Snow Data and Detecting Climate Change Impacts in the Mountainous Eastern Regions of Turkey" Water 13, no. 14: 1982. https://doi.org/10.3390/w13141982

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