Analyzing the Effects of Climate Change for the Water Tower of Mesopotamia, Türkiye

. This study aims to analyze future snowmelt runoff and snow extent for the Upper Euphrates Basin (41 109 km 2 ) which is the most upstream basin of the water tower Mesopotamia and located in the Eastern mountainous part of Türkiye. The future characterization of snowmelt runoff and snow water equivalent is investigated taking the regional climate model projections into account through hydrological model applications. The basin is divided into two major sub-basins Karasu Basin and Murat Basin. The conceptual hydrological models, HBV and HEC-HMS have been utilized to establish a rainfall-runoff relationship. The models are calibrated and validated with observed daily total precipitation and average temperature data for the 1980–2011 period. Future projections are selected for the 2025–2100 periods using EU-CORDEX data with various global circulation models. MPI-ESM-LR and CNRM-CM5 data sets and RCP4.5 and RCP8.5 emission scenarios are used for future projections. The results indicate a decrease in snow cover extent and an increase in snowmelt runoff. The overall assessment will help us to understand how these changes will affect the operations of water resources systems (downstream reservoirs) in terms of ﬂood control, energy production, irrigation


Introduction
Global climate change greatly affects climatological and hydrological elements like temperature, precipitation regime, groundwater levels and hydrological response (Gleick, 1989).Although floods and droughts are/will be the main expected consequences in the future, also it is expected that global change will have different effects on each region (Tabari, 2020).Therefore, the assessment of potential impacts of climate change on regional water resources is receiving growing attention from various researchers.The Mediterranean Basin has been identified as one of the world's most sensitive regions to climate change consequences (URL-1, 2022).Fujihara et al. (2008) explored the potential impacts of climate change on the Seyhan River Basin in Türkiye and concluded an increase in temperature and a reduction in precipitation for future projections.Bozkurt and Sen (2013) investigated future climate change in the Euphrates-Tigris Basin using dynamically downscaled outputs of different Global Circulation Models (GCMs) and emissions scenario simulations for two 30 year periods, 2041-2070 and 2071-2099.According to their analyses, significant decreases in the annual surface runoff is expected.
GCMs are generally used in large-scale climate change impact studies and play a key role in recognizing future global climatic changes (Fujihara et al., 2008).The model outputs necessitate a proper hydrological modeling background to produce future projections.The eastern part of Türkiye is covered with full or partially snow during winter and large rivers such as Euphrates and Tigris are mainly fed by the snowmelt.This study aims to project future runoff and snow components in the snow-dominated basins of the Euphrates River basin by integrating climate model data sets and two well-known conceptual hydrological models.

Data and Methodology
This section briefly describes the data used for the set-up and the structures of hydrological modelling approaches together G. Uysal et al.: Analyzing the Effects of Climate Change for the Water Tower of Mesopotamia, Türkiye with the climate data and their assessment for the future projections.Hydrological modelling is accomplished with two conceptual models: HBV-Light and HEC-HMS.After calibrating and validating models, GCM data set based runs are employed to simulate the historical and future changes.The changes over the reference period are calculated for runoff and snow components.

Hydro-Meteorological and Climate Data
Türkiye has a complex topography and variable climate.Two upstream snow-dominated catchments (Fig. 1) are selected; Karasu Basin (15 330 km 2 ) and Murat Basin (25 779 km 2 ) having the elevation ranges between 860-3565 and 840-4045 m, respectively.Annual total precipitation is 429.3 and 593.9 mm yr −1 , annual average temperature is 9.1 and 8.2 °C and annual average discharge is 153.5 and 239.0 m 3 s −1 between 1980-2011 for Karasu and Murat, respectively.These basins play an important role in downstream water structures.The region is also the concern of various national and international projects ( Şensoy and Uysal, 2012;Uysal et al., 2016;Şorman et al., 2019;Uysal and Şorman, 2021).
Historical simulations and projections of future changes in precipitation and temperature are provided by the EURO-CORDEX platform (URL-2, 2021).MPI-ESM-LR and CNRM-CM5 GCMs data sets having RCP4.5 and RCP8.5 concentration scenarios with 112.5 km spatial resolution are selected in the study.Data over the basin are derived from Regional Climate Model (RCM) ModeRegCM4.3.4 to provide a series of high-resolution inputs (20 km resolution).Linear bias correction is applied for precipitation and temperature data sets.Stationary climate simulations are used for reference  and future (2025-2100) periods.The future projection period is evaluated under 2025-2050, 2051-2075, and 2076-2100 periods.

HBV Hydrological Model
The Swedish Meteorological and Hydrological Institute (SMHI) developed HBV as a conceptual and continuous hydrological model (Bergström, 1995).The HBV-Light version (Seibert and Vis, 2012) was utilized in this study as it offers a user-friendly graphical interface and an automatic calibration opportunity.The HBV model distributes precipitation, temperature, and potential evapotranspiration inputs based on elevation and vegetation zones.The model calculates snow and soil moisture layers based on zones, as well as runoff and routing processes for the entire basin, to mimic daily flow.

HEC-HMS Hydrological Model
HEC-HMS (Hydrologic Modeling System) model is developed by the U.S. Army Corps of Engineers Hydrologic Engineering Center (USACE) and is designed to simulate the rainfall-runoff process of dendritic watershed systems as continuous or event-based and to be applicable in a wide range of geographic areas (Scharffenberg and Fleming, 2000).Inputs are precipitation, temperature, and potential evapotranspiration.In this study, soil moisture accounting (SMA), linear reservoir and Clark unit hydrograph method are chosen respectively as a loss, baseflow and transform methods.SMA module considers continuous situation-based soil moisture variation which is carried out by calculating the amount of water stored in and moving through vegetation, soil surface, profile and groundwater layers.Linear reservoir method is more compatible with SMA.

Assessment Metrics
Nash-Sutcliffe Efficiency (NSE), Killing-Gupta Efficiency (KGE) and Pearson correlation coefficient (r) are used to test hydrological model performances (Eqs.1-4): where x is observation, y is modelled value, x is average observation, y is average modelled value, S x is the standard deviation of the observed data, S y is the standard deviation of the modelled data, n is the number of data sets.Contingency tables evaluate the consistency by using two main indicators i.e., Probability of Detection (POD) (Eq.5) and False Alarm Rate (FAR) (Eq.6).

POD
where "a" (detected by both ground and GMC), "b" (detected by only GCM), "c" (detected by only ground) and "d" (no detection in both).

Hydrological Modelling Performances
The conceptual daily hydrological models are calibrated and validated from 1980-2011.Tables 1 and 2 show that the HBV model has a higher model efficiency for both basins and periods.The performance of HBV is slightly better than the HEC-HMS model in terms of NSE, mainly associated with peak flow accuracy.KGE accounts for various components and both models show similar performances in Karasu.: 1980: -1999: , Val: 1991: -2011. .

Raw and Bias-Corrected Climate Data Performances
Bias correction is applied for temperature and precipitation GCM data sets.Performances of the raw and bias-corrected GCMs data sets are compared with gauge measurements for the reference period.Pearson-R values of monthly average bias-corrected temperature against gauge measurement are around 0.96 for MPI-ESM-LR and CNRM-CM5 in both basins.Figure 2 depicts both basins' POD and FAR values with monthly (mo) threshold precipitations.POD values increased in both products after the bias correction application.CNRM-CM5 performs slightly better in comparison with MPI-ESM-LR.

Future Projections
The performances of the hydrological models, GMCs data are assessed and the change in future projections (runoff and SWE) are demonstrated against the reference period.Future projections are assessed for RCP8.5 as the worst-case scenario (Fig. 3).Both models demonstrate a possible increase in the runoff for two mountainous basins (more emphasized for Karasu Basin) for the first two projection periods and a significant decrease for the last projection period according https://doi.org/10.5194/piahs-385-465-2024Proc.IAHS, 385, 465-470, 2024 These mismatches might be assigned to variations of precipitation in each GCM in each period.Figure 4 depicts the change in snow water equivalent (SWE) in terms of percentages in the future projections.Regardless of the hydrological model and GCM, the reduction is expected in future SWE values.Fewer reduction is expected for the 2024-2054 period.The reductions are exacerbated along with a longer lead time period.This might be attributed to the increase in temperature in the future.

Discussion of the Results
A bias correction is required for GCM data sets, however the improvement for precipitation data is limited compared to the temperature data.Both hydrological models (HBV and HEC-HMS) provide satisfactory calibration and validation performances.The change in runoff is variable depending on GCM data set.For example, MPI-ESM-LR based runoff present no change or slight increase for change in percent runoff between 2024-2054 and 2055-2084, however an increase is expected for 2085-2099.On the other hand, the change for CNRM-CM5 is either vice versa or less.The future projection of the SWE is decreasing related to the increase in temperature data.Climate change will alter snowfall to rainfall, thus will have various impacts on surface runoff and snow extent.

Conclusions and Recommendations
The study investigates the future implications of snowmelt runoff against the changing climate for the upstream part of the Euphrates River.EU-CORDEX GCM data sets MPI-ESM-LR and CNRM-CM5 are accounted to estimate future projections of discharge and SWE.The study accounts different hydrological models (HBV & HEC) in the assessment.Both hydrological models provide satisfactory validation NSE values (0.84-0.90).The bias correction is required for GCM data sets before using them in the models.A possible increase in the runoff for the first two projection period, and a significant decrease for the last projection period especially according to MPI-ESM-LR are expected.These results indicate a drastic and continuous decrease in SWE due to changes in the temperature.These changes will negatively affect the operation of reservoirs, the agriculture policy, and tourism management.Nonetheless of model selection, change in precipitation holds more uncertainty which might cause the variation in streamflow runoff.Changes in snow line elevation and lag time are being assessed in future studies.Considering the variation in GCMs, increasing the number of climate models would provide more information to interpret the results.Author contributions.GU analyzed the data, wrote the draft; YOD and HSC performed the model simulations; AS designed the research; GU, YOD, HSC, AS, AAS reviewed and edited the manuscript.
Competing interests.The contact author has declared that none of the authors has any competing interests.

Figure 1 .
Figure 1.Snow cover across Türkiye from IMS satellite (30 March 2009) and the selected catchments on the Euphrates River.

Figure 3 .
Figure 3.The change in runoff for three future projections (a) Karasu Basin, (b) Murat basin.

Figure 4 .
Figure 4.The change in snow water equivalent (SWE) for three future projections (a) Karasu Basin, (b) Murat basin.
Disclaimer.Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper.While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.Special issue statement.This article is part of the special issue "IAHS2022 -Hydrological sciences in the Anthropocene: Variability and change across space, time, extremes, and interfaces".It is a result of the XIth Scientific Assembly of the International Association of Hydrological Sciences (IAHS 2022), Montpellier, France, 29 May-3 June 2022.Acknowledgements.Special thanks are extended to the State Hydraulics Works of Türkiye (DSI) as well as the Turkish State Meteorological Service (MGM).Financial support.This research has been supported by Eskisehir Technical University (Scientific Research Projects 22DRP190 and 22ADP051).https://doi.org/10.5194/piahs-385-465-2024Proc.IAHS, 385, 465-470, 2024

Table 1 .
Model performances for Karasu basin.

Table 2 .
Model performances for Murat basin.
* Cal Data availability.Publicly available free climate datasets (https:// esg-dn1.nsc.liu.se/search/cordex/,ESGF, 2024) are analyzed in this study.The hydro-meteorological data set is not public or open to third parties.