A Hydrological Modeling Approach for Assessing the Impacts of Climate Change on Runoff Regimes in Slovakia

.


Introduction
Understanding the impact of climate change on the water balance is one of the main scientific interests in hydrology.Possible changes in runoff conditions caused by changing climatic conditions are one of the significant sources of uncertainty in water resources management and flood protection.Climate change can cause a decrease in water resources, changes in the regimes of runoff and snow accumulation, and increases in the extremes of floods and droughts (e.g., [1][2][3][4][5]).Changes in hydrological time series in the observation period in Slovakia were detected by Ďurigová, et al. [6,7], and changes in the water balance and the soil moisture regime due a changing climate were analyzed by Rončák at al. [8,9].Therefore, it is important to estimate the possible impacts of future climate change on water resources.
Rainfall-runoff (r-r) models, especially conceptual r-r models, are widely used for climate impact assessments (e.g., [10][11][12]).There are numerous regional [3,5,[13][14][15][16] and national (e.g., [17][18][19][20][21][22][23][24][25]) studies that have evaluated the impact of climate change on water resources using r-r models, which are forced with different global or regional climate circulation models, and which consider various emission scenarios.Most of them have shown that climate change will play an important role in future runoff conditions.For example, Vyshnevskyi et al. [25] assessed the possible climate change impact on the runoff regime at 11 meteorological stations in the Ukrainian Carpathians.Their results showed a tendency for runoff to increase in a winter period.Wang et al. [18] showed that increases in precipitation lead to greater changes in runoff.Their research was performed in the Liao River basin in China.Wu et al. [22] showed a reduction (11.6%) in annual runoff in the Chinese Yellow River basin according to their hydrological model's projection.This is in line with Huntington [17], who showed that climate warming could reduce runoff in 38 forested watersheds in New England (the United States).According to their results, the largest reductions in runoff were observed during April and May.Shahvari et al. [24] evaluated the impact of climate change on water resources through the use of a SWAT hydrological model under three scenarios (RCP 2.6, RCP 6, RCP 8.5).The study was performed for three basins (the Lavasanat, Damavand, and Varamin) in Iran.Their results showed that runoff will increase in the fall/spring and decrease in the winter/summer under future time horizons.Aich et al. [20] assessed the possible impact of climate change on the runoff regime in three African river basins (the Niger, the Upper Blue Nile, the Ouban).Their results indicated the impact of climate change on runoff, especially on high/low flows and mean discharges.The impact of future climate changes on runoff was also presented by Vano et al. [21], who used the variable infiltration capacity (VIC) hydrological model in the Pacific Northwest (USGS Hydrologic Region 17 plus the Canadian portion of the Columbia).Their results indicated seasonal shifts in runoff in response to warming during the cold season.
Most of these results are comparable with the results of previous studies performed in various Slovak regions (e.g., [5,8,[13][14][15][16]).For example, Kohnová et al. [5] evaluated the impact of climate change on extreme runoff regimes in two Slovak basins (the Váh and Laborec basins).They predicted a general increase in mean monthly runoff, especially during the winter period.These results are supported by the results of Rončák et al. [3] and Rončák et al. [16], who indicated that overall monthly runoff will increase (mainly during the winter period) with changes in climate as a result of increases in the precipitation/air temperature.
Despite the fact that quantifying the impact of climate change on runoff processes has been extensively investigated over recent decades, this topic is still challenging (because of complex runoff mechanisms, etc.) for hydrologists.The scientific objective of this study was to develop a calibration strategy based on the selection of the most suitable decade in the observation period for the parameterization of a rainfall-runoff model for modelling runoff in future climate conditions.The parameterized model was then applied to estimate the possible effect of climate change on the runoff regime in the eight selected basins in Slovakia, which were divided into two groups according to their mean basin elevation.In extending previous studies performed in these Slovak regions, this paper concentrates on an approach for studying the effects of climate change on the whole territory of Slovakia.Due to a strong climate signal in the recent decades (e.g., [26]), this region is well suited for such analyses.We applied a lumped hydrological model based on the HBV (Hydrologiska Byråns Vattenbalansavdelning) concept, which is widely used for climate impact assessments.The climate projections were represented by the model runs forced by the KNMI (Koninklijk Nederlands Meteorologisch Instituut) and MPI (Max Planck Institute) regional circulation models with the A1B emission scenario.
The paper is structured as follows.The study basins, data, and climatic scenarios are presented in Section 2. Section 3 describes the rainfall-runoff model.Section 4 contains the final results of the study, and Section 5 contains the discussion, followed by the conclusions in Section 6.

Data
We used the daily basin data (mean daily runoff, daily precipitation totals, mean daily air temperature) from the period 1981-2019.This data set was provided by the Slovak Hydrometeorological Institute (SHMI).Visual inspection showed that the data from all the river basins selected do not have any observation gaps.Point precipitation data measurements were used to interpolate the basins' mean areal precipitation using the inverse distance weighting (IDW) method [27], and the point air temperatures were processed to the basins' mean areal values by the altitude gradient method.The daily potential evapotranspiration was calculated by a modified Blaney-Criddle method [28].
Figure 2 and Table 2 show the changes in the climatic characteristics over four decades for all the basins studied, which are divided in two groups according to their mean basin elevations.We can see that the medians of the long-term mean annual air temperatures (T) show an increasing trend over the four decades.In the basins with a lower mean elevation, the median of T increased from 7.8 °C to 9.4 °C between the periods 1981-1990 and 2011-2019, i.e., by +1.6 °C; in the basins with a higher mean elevation, the median of T increased from 5.6 °C to 7.2 °C between the periods 1981-1990 and 2011-2019, i.e., by +1.6 °C.
The long-term mean annual precipitation (P) median values in the basins with a lower mean elevation increased from 657 to 750 mm/year between the periods 1981-1990 and 2001-2010, and from 657 to 661 mm/year between the periods 1981-1990 and 2011-2019, i.e., by +4 mm/year.In the basins with a higher mean basin elevation, the median values of the long-term mean annual precipitation (P) increased from 794 to 899 mm/year between the periods 1981-1990 and 2001-2010, and from 794 to 847 mm/year between the periods 1981-1990 and 2011-2019, i.e., by +53 mm/year.
The medians of the long-term mean annual runoff (Q) in the basins with a lower mean elevation show slightly increasing or decreasing trends over the decades and an overall decrease from 219 to 184 mm/year between the periods 1981-1990 and 2011-2019, i.e., by −35 mm/year.In the basins with a higher mean elevation over the four decades, a slight increase in the runoff median from 368 to 382 mm/year between the periods 1981-1990 and 2011-2019, i.e., + 14 mm/year, was observed.
Table 2.The medians (50%) and percentiles (25% and 75%) of the long-term mean annual precipitation (P), air temperature (T), and runoff (Q) averaged over four decades (a) for the river basins with a mean elevation < 435 m a.s.l., (b) for the river basins with a mean elevation > 435 m a.s.l.Changes in the hydroclimatic characteristics (long-term mean annual precipitation, air temperature, and runoff) over four decades.A line represents the medians of these characteristics for each of the four periods (i.e., 1981-1990, 1991-2000, 2001-2010, and 2011-2019) and for the two groups of the basins.The scatter (i.e., the 75 and 25% percentiles) indicates the variability between the basins.

Climate Change Scenarios
The data (i.e., the mean daily air temperatures and daily precipitation totals) generated from the regional circulation models (RCMs), i.e., KNMI (Koninklijk Nederlands Meteorologisch Instituut) and MPI (Max Planck Institute), under the A1B (moderate) scenario were applied in the hydrological model to project the future runoff conditions.Both RCMs (i.e., the KNMI and MPI) are based on the German ECHAM5 general circulation model (GCM) boundary conditions [26], and represent a detailed integration of the atmospheric and oceanic dynamic equations with a grid point resolution of about 25 × 25 km.The outputs from both climate models for the region studied were downscaled to the climatic stations and interpolated using the IDW (Inverse Distance Weighted) and altitude gradient methods on the basin's data.The values of the mean daily air temperatures and daily precipitation totals were used as the daily basin data for the period between 2011 and 2100.
Figure 3 and Table 3 show the long-term mean monthly values of the precipitation (P) and air temperature (T) for the reference period (1981-2010) and three future time horizons (2011-2040, 2041-2070, and 2071-2100) according to the KNMI and MPI climate scenarios, separately for the river basins with lower and higher mean elevations.The lines in the graphs indicate the averages for the basins.We can see that the changes in the long-term mean monthly climatic characteristics (T and P) according to both climate models used are expected to be more or less the same for both groups of the river basins.During the spring, autumn, and winter periods, an increase in precipitation can be expected.Conversely, a decrease in precipitation can be observed in the summer periods.Based on the KNMI model, the largest increase in precipitation is expected to be in September (period 2071-2100).This may also be caused due to the changes in seasonality when the maximum of the long-term mean monthly precipitation shifted from July to September in all future time periods.The MPI model predicts the largest increase to be in March (period 2071-2100).The long-term mean monthly air temperatures will rise with the future time horizons.The KNMI climate model predicts an increase in the air temperature during the winter months by 3 °C, and in the summer seasons by 4 °C in the future time horizons.In the spring and autumn seasons, the air temperature is expected to increase by about 3 °C by 2100.According to the KNMI scenario, the largest increase in air temperature will be observed in July.The highest increase in air temperatures according to the MPI scenario is assumed to occur in July/August.

Hydrological Model
To study the possible effects of climate change on the hydrological regime in the selected Slovak basins, a conceptual rainfall-runoff (r-r) model (the TUW model, [29]) has been chosen.The TUW model is a lumped model which follows the structure of the widely used Swedish Hydrologiska Byråns Vattenbalansavdelning (HBV) model [30,31].The model has proven to be suitable for the analysis of the impact of climate change on water resources (e.g., [10][11][12]).The model uses the basin values of mean daily air temperatures, daily precipitation totals, and potential evapotranspiration as the inputs.
The model consists of three routines, namely: (1) the snow routine, which uses a simple degree-day method for calculating snow accumulation and snow melt processes in the river basin; (2) the soil moisture routine, which simulates the processes in the soil profile (e.g., the evapotranspiration of water from a soil profile, the infiltration of water from melted snow into the soil, the distribution and accumulation of water in soil, and the formation of surface, subsurface/groundwater runoff); and (3) the runoff generation routine, which is used to transform the total runoff (surface runoff, subsurface runoff and baseflow) through a river basin.The stream routing uses a triangular weighting function for the runoff propagation to the basin outlet.
The model has 15 parameters that need to be calibrated.According to the three routines and their calibration ranges, the model parameters are listed in Table 4.The calibration ranges were specified based on the existing studies (e.g., [32,33]).
More details about the model's structure and the routines/parameters used to calculate the individual runoff components are given, e.g., by Parajka et al. [29], Sleziak et al. [11,34], and Széles et al. [35].In this study, the TUW model is calibrated by using the differential evolution algorithm (Deoptim) automatic calibration procedure [36,37].The selected algorithm uses biologically inspired processes (i.e., crossover, mutation, and selection on a population) in searching for the best results.The optimization process used in the Deoptim algorithm is based on the repeated evaluation of the objective function in order to move the initial population toward a global minimum.For further information about the Deoptim algorithm, see, e.g., Ardia et al. [36] and Mullen et al. [37].As the objective function (OF), we apply a combination of the Nash-Sutcliffe coefficient (NSE) and the logarithmic Nash-Sutcliffe coefficient (logNSE).Whereas the NSE accentuates the high flows, the logNSE emphasize the low flows more.The OF is mathematically defined as: where the NSE and logNSE criteria are estimated as: , are the simulated mean daily runoffs on day i; i obs Q , are the observed mean daily runoffs; and obs Q is the average of the observed runoffs over the calibration/validation period of n days.The NSE and logNSE coefficients vary between −∞ (a poor match between the observed and simulated values) and 1 (a perfect match).
The performance of the model is judged by comparing the observed and simulated runoff in terms of the above-defined model performance metrics (OF, NSE, logNSE) for the calibration/validation periods.The second measure of the model's performance used in this study is the volume error (VE), which is the measure of the bias between the simulated and observed runoff [22].VE = 0 indicates no bias; VE < 0 implies an underestimation of the total runoff volume; and VE > 0 represents an overestimation of the total runoff volume.The VE is estimated as follows: We used the differential split-sample test (DSST) proposed by Klemeš [38] for the testing of the model's performance.The model was calibrated/validated consecutively for four 10-year periods between 1981-2019.For further analysis, i.e., for simulating the future runoff conditions, we choose parameters from the recent calibration period (i.e., 2011-2019; Table 5) because we assume that this period should be similar (mainly in terms of the mean daily air temperatures) to a recent/warmer climate.

Assessment of the Model's Performance
The results of the model's calibration and validation in terms of the four different model performance metrics (the objective function (OF), the logarithmic Nash-Sutcliffe efficiency (logNSE), the Nash-Sutcliffe efficiency (NSE), and the volume error (VE)) for all the basins studied are shown in Figure 4. Generally, in terms of adequate model performance, the OF, NSE, and logNSE should be close to 1, and the VE should be close to 0.
The results indicate that the model performed better in the calibrations compared to the validations.These expected findings were caused mainly by the differences among climatic and hydrological characteristics in the four decades in the observation period 1981-2019.During the calibration period (2011-2019), the OF values for all the basins ranged from 0.72 (the Váh River basin) to 0.84 (the Hron River basin).The volume errors from -0.5 % (the Turiec River basin) to 4.9 % (the Váh River basin) show a trend to slightly overestimate the runoff volume in the calibrations.The OF values in the validations are as follows: in the validation period of 1981-1990, the OF values were between 0.55 (the Myjava River basin) to 0.77 (the Hron River basin).In the validation period of 1991-2000, the OF values varied from 0.49 (the Myjava River basin) to 0.84 (the Hron River basin).In the validation period of 2001-2010, the OF values ranged from 0.41 (the Myjava River basin) to 0.78 (the Hron River basin).Concerning all of the model performance metrics (the OF/NSE/VE values), the best model fit was found for the Hron River basin.Conversely, poorer simulations were obtained for the Myjava River basin.For this small basin, the VE indicated a more than 20% overestimation of the runoff volumes in the validation period of 1991-2000, and a more than 30% overestimation in the validation period of 2001-2010.
Due to the differences in OF in the calibration and validation periods, for simulating the future runoff, we chose parameters from a recent calibration period (i.e., 2011-2019) because we assume that this period should be similar (mainly in terms of the mean daily air temperatures) to a future/warmer climate.
The simulated discharges of the hydrological model driven with the two climate models were compared with the observed discharges to visually check the model's performance for the reference period (Figure 5).The visual comparison indicates that the model provided a reasonable simulation of the runoff.We can see that the seasonal pattern of the runoff dynamics observed was well captured.However, there are some differences in representing high flows, which can especially be seen in the months of March, April, and May.For example, the largest increase in runoff for the Hron River basin is observed in April (Qobs = 70 mm/month).Simulations (Qsim) obtained for this month are 64 mm/month (Qsim), 55 mm/month (Qsim considering the KNMI scenario), and 68 mm/month (Qsim considering the MPI scenario).A similar pattern (i.e., differences in representing high flows) has been observed for most of the river basins studied.
A comparison of the different scenarios points that they are similar in terms of simulating an increase of high and low runoff.Overall, it is obvious that the simulated runoff using downscaled inputs from both RCM models for the reference period represents the simulated hydrological regime from the measured data well, and, at the same time, the model parameters from the measured inputs are usable for an acceptable reproduction of the hydrological regime of the reference period.Therefore, we were able to use these data as a basis for comparing the changes in runoff according to the scenarios.

Impact of Climate Change on Runoff
The predicted future changes in runoff due to climate change were evaluated by comparing the simulated long-term mean monthly runoff for the reference period (1981-2010) and the modelled scenarios in three time periods (2011-2040, 2041-2070, and 2071-2100) (Figures 6 and 7).A comparison of the simulated long-term mean monthly runoff between the reference period and projected scenarios is presented in percentages (Table 6).The values that decreased by 10% or more than 10% are marked in gray.
The results indicate that changes in the runoff regime in all the basins analyzed could be expected in the future, and are assumed to be more or less the same according to both climate scenarios used.Both climate scenarios foresee an increase in the winter runoff (December-February) with a maximum increase in the period 2071-2100.The highest increase in runoff is expected for the Topľa River basin in January (219% for the KNMI scenario and 141% for the MPI scenario).This could be caused by increases in the air temperature (snowfalls will change into rain).On the other hand, the runoff will react in the reverse in the summer season (June-August).According to the KNMI scenario, the runoff will decrease by 28% up to 52%.
The highest percentage decrease in the summer runoff (56% in July) according to the KNMI scenario is assumed for the Myjava River basin.
For the Nitra River basin, we can expect a maximum increase in runoff in the winter months of December and January (according to the KNMI scenario, up to 94%; and according to the MPI scenario, up to 78%).Decreasing runoff (up to 43% for the KNMI scenario) in this basin mainly will occur from June to August.For the Laborec River basin, the most significant decrease (up to 43%) under the scenarios is shown in the spring and summer months.According to both scenarios, the runoff is expected to increase from December to February.
A similar pattern (i.e., a decrease in runoff in the summer months and an increase in the winter months) is observed for the other basins.For example, according to the MPI scenario, the runoff in the Hron River basin may amount to a more than 100% increase in January/February in the last time horizon (2017-2100).According to the KNMI scenario, the runoff will increase by 96-104% in the months of January to February.Conversely, the runoff will decrease by 18% up to 38% in the months of May to August.
A decrease in the summer runoff (50% in August for the KNMI scenario) is projected for the Poprad River basin.The MPI scenario also presents a decrease, but a little bit lower (from 19-24%).
For the Váh River basin, an increase in runoff is projected from December to March (this is true for the climatic data from both climatic models).This is a high-altitude basin, and an increasing air temperature here could result in less snow accumulation and more intensive runoff from snow melt.
Similar to the other basins, an increase in winter runoff (from December to March) can be seen in the Turiec River basin.For this river basin, the highest increases in the runoff are observed in January (according to the KNMI scenario, 101% for the 2071-2100 horizon; and according to the MPI scenario, 98% for the period 2071-2100).Decreasing runoff is mainly observed during the summer months.
According to the results presented, we can conclude that for the entire territory of Slovakia, both climate change scenarios (the KNMI and MPI) predict similar seasonality changes.They foresee an increase in runoff in the winter/spring and a decrease in the summer due to changes in the precipitation regime and an increase in the air temperature, with the intensity of the changes increasing towards the future horizons.The results indicate that the medians of the maximum annual daily runoff tend to be higher in 2071-2100 for the climatic data from the KNMI model and the river basins with lower mean basin elevations.According to the MPI scenario, the medians of the maximum annual daily runoff show a slightly increasing trend until the year 2040, with an increase of about 0.68 mm/day between the periods 1981-2010 and 2011-2040.The maximum annual daily runoff does not change much during the four periods.
From the analysis, it is also evident that in the basins with lower mean basin elevations, the values of maximum annual daily runoff are higher in size and frequency compared to the basins with higher mean basin elevations.On the other hand, a few very extreme values of maximum annual daily runoff were simulated in the basins with higher mean basin elevations.These extremes may be caused by an increase in the air temperature in the future and extreme runoff from the intensive melting of snow here.Table 6.Simulated long-term mean monthly runoff (mm/month) and its changes (%) for the three future time horizons: (a) in the river basins with a mean elevation < 435 m a.s.l.; (b) in the river basins with a mean elevation > 435 m a.s.l.

Discussion
This study evaluates the potential impacts of climate change on the runoff regime in eight selected basins located in the whole territory of Slovakia, divided into two groups according to their mean basin elevations.The comparison of changes in the climatic and hydrological characteristics over four decades in the observation period showed differences in these river basins.We analyzed the same increase in the medians of the long-term mean annual air temperatures in the latest decade, 2011-2019, compared to the first decade, 1981-1990, for both groups of the basins.In the river basins with a lower mean elevation, we can see only a very slight increase in the long-term mean annual precipitation and a slight decrease in the long-term mean annual runoff in the latest period, 2011-2019, compared to the first period, 1981-1990.On the contrary, in the river basins with a higher mean basin elevation, an increase in both the long-term mean precipitation and in the long-term mean annual runoff in the latest period compared to the first period is evident.
The climate projections were represented by the TUW rainfall-runoff model runs forced by the KNMI and MPI climate models with the A1B emission scenario.Because we cannot calibrate the rainfall-runoff model parameters for the future time horizons, the model was calibrated and validated for four decades in the observation period 1981-2019.The differences in objective functions in the calibration and validation periods confirmed that the model parameters highly depend on input climate data which changed over the four decades.Therefore, for simulating the future runoff, we chose parameters from a recent calibration period 2011-2019, which should be the most similar, mainly in terms of the mean daily air temperatures to a future warmer climate.
According to the KNMI and MPI models evaluated in the reference period (1981-2010) and three future time horizons (2011-2040, 2041-2070, 2071-2100) for all the river basins, the simulated long-term mean monthly runoff shows the known river basin response to climate change from previous works (e.g., [3,5,[13][14][15][16]25]). The runoff will rise in the winter/spring, but drop in the summer, throughout the future time horizons.The variability of the runoff can be attributed to the variability of the rainfall and increasing air temperatures towards the future horizons (e.g., [3,16]).Higher air temperatures could increase winter snow-melt runoff due to less snow accumulation.The monthly precipitation totals may not change significantly, but the extremes may increase (alternating low-precipitation periods with periods featuring more extreme precipitation).
Several studies on climate change impacts on the runoff regime have been conducted in recent years (e.g., [23][24][25]39]).Some of these studies (e.g., [23,24,39]) observed an increase in future runoff, whereas others (e.g., [17,[40][41][42]) observed a decrease in runoff.For instance, Worqlul et al. [23], in their study, showed that the runoff will rise by up to 64% in dry periods and decrease by 19% in wet periods.Their research was performed in two subcatchments in the upper Blue Nile catchment (Ethiopia).This is in line with Liersch et al. [39], who, for the same study area (upper Blue Nile catchment), presented the increasing pattern in average annual discharges.The authors also observed a shift in seasonal discharge patterns (increase from August to November, and decrease in June and July).In another study, Shahvari et al. [24] showed that runoff will increase in the fall/spring and decrease in the winter/summer according to their hydrological model's projection.Results from these as well as many other studies suggest that climate change will play an important role in future runoff conditions.
Our results indicate that for most of the basins, studied increases in runoff could occur in the December to February period.On the other hand, decreases in runoff could occur from March/April to September.Such results have also been documented, e.g., Štefunková et al. [14] and Rončák et al. [3].Those studies were conducted for three and five Slovak basins, respectively.Another study on the impact of climate change on discharge seasonality for selected Slovak basins can be found, e.g., in Hlavčová et al. [13], Kohnová et al. [5], and Rončák et al. [16].We used a larger sample of basins (eight basins) in our investigation; therefore, the results can be generalized for the whole territory of Slovakia.In contrast to previous studies undertaken in Slovak basins, we also evaluated the maximum annual daily runoff for the reference period and future horizons.Our evaluation showed that in the reference period and future time horizons, in the basins with lower mean basin elevations, the values of maximum annual daily runoff are higher in size and frequency compared to the basins with higher mean basin elevations.On the other hand, a few very extreme values of maximum annual daily runoff were simulated in the basins with higher mean basin elevations.Generally, the maximum annual daily runoff could be more extreme for later time horizons (according to the KNMI scenario for 2071-2100).
The results of our study show how the hydrological regime might change in the basins studied in the future.This analysis could help develop optimal strategies for better management of available water resources in the territory of Slovakia.

Conclusions
For estimating the impact of climate change on the runoff regime in Slovakia, eight pilot river basins spread throughout the whole territory of Slovakia were selected and divided into two groups according to their mean basin elevations.Analyzing the climatic characteristics over four decades in the observation period, we can see that the medians of the long-term mean annual air temperatures show an increasing trend.The median values of the long-term mean annual precipitation in the river basins with a lower mean elevation increased between the first and the last decade by 4 mm/year, and in the basins with a higher mean basin elevation by 53 mm/year.The long-term mean annual runoff medians in the river basins with a lower mean elevation show an overall decrease over the four decades by 35 mm/year and a slight increase by 14 mm/year in the basins with a higher mean elevation.
A lumped hydrological model based on the HBV concept was applied for modelling runoff in new climate conditions.For the parameterization of the model for future climate, we used a calibration strategy based on selecting the most suitable decade in the observation period.The differences in objective functions in the calibration and validation periods in four decades in the observation periods confirmed differences in the model parameters due to different climate characteristics in these decades.Therefore, for simulating the future runoff, we chose parameters from a recent calibration period (i.e., 2011-2019) because we assume that this period should be similar (mainly in terms of the mean daily air temperatures) to a future/warmer climate.
The climate projections were represented by the model runs forced by the KNMI and MPI regional circulation models with the A1B emission scenario.The model's performance metrics and a visual inspection showed that the simulated runoff using downscaled inputs from both RCM models for the reference period represents the simulated hydrological regimes well.
The outcomes of the runoffs simulated under changing climate conditions show that changes in the seasonal distribution of the long-term monthly runoff can be expected.The monthly runoff is predicted to increase during winter periods and decrease during summer seasons.These changes are expected to be more extreme for the later time horizons.The maximum annual daily runoff also could be more extreme for the last time horizon.
Our results have shown that climate change will impact the runoff regimes of the eight basins studied, with consequences for water resource planning/management.The results from this study can be used to support decision-making processes, for example, in developing optimal strategies for adapting to the changing climate predicted.

Figure 1 .
Figure 1.Topography of Slovakia and locations of the eight river basins selected.

Figure 2 .
Figure 2. Changes in the hydroclimatic characteristics (long-term mean annual precipitation, air temperature, and runoff) over four decades.A line represents the medians of these characteristics for each of the four periods (i.e.,1981-1990, 1991-2000, 2001-2010, and 2011-2019)  and for the two groups of the basins.The scatter (i.e., the 75 and 25% percentiles) indicates the variability between the basins.

Figure 3 .
Figure 3. Seasonal variability of the long-term mean monthly air temperature and precipitation totals for the KNMI and MPI climate scenarios: (a) for the river basins with a mean elevation < 435 m a.s.l.; (b) for the river basins with a mean elevation > 435 m a.s.l.The lines show the basin averages.

Figure 4 .
Figure 4. Values of the objective function (OF), the logarithmic Nash-Sutcliffe efficiency (logNSE), the Nash-Sutcliffe efficiency (NSE), and the volume error (VE) in the calibration/validation periods of the eight basins studied.

Figure 5 .
Figure 5.Comparison of the measured and simulated long-term mean monthly runoffs for the reference period (1981-2010) and for the two groups of river basins.

Figure 6 .
Figure 6.Long-term mean monthly runoff projected for the reference period (1981-2010) and three future time horizons (2011-2040, 2041-2070, 2071-2100) according to the KNMI and MPI scenarios for the river basins with a mean elevation < 435 m a.s.l.

Figure 7 .
Figure 7.Long-term mean monthly runoff projected for the reference period (1981-2010) and three future time horizons (2011-2040, 2041-2070, 2071-2100) according to the KNMI and MPI scenarios for the river basins with a mean elevation > 435 m a.s.l.

Figure 8
Figure8compares the maximum annual daily runoff for the reference period(1981- 2010)  and three future time horizons(2011-2040, 2041-2070, and 2071-2100).The results indicate that the medians of the maximum annual daily runoff tend to be higher in 2071-2100 for the climatic data from the KNMI model and the river basins with lower mean basin elevations.According to the MPI scenario, the medians of the maximum annual daily runoff show a slightly increasing trend until the year 2040, with an increase of about 0.68 mm/day between the periods 1981-2010 and 2011-2040.The maximum annual daily runoff does not change much during the four periods.From the analysis, it is also evident that in the basins with lower mean basin elevations, the values of maximum annual daily runoff are higher in size and frequency compared to the basins with higher mean basin elevations.On the other hand, a few very extreme values of maximum annual daily runoff were simulated in the basins with higher mean basin elevations.These extremes may be caused by an increase in the air temperature in the future and extreme runoff from the intensive melting of snow here.

Figure 8 .
Figure 8.Comparison of the simulated maximum annual daily runoff for the KNMI and MPI climate change scenarios for the reference period (1981-2010) and future horizons (2011-2040, 2041-2070, and 2071-2100.The lines show the basin averages.The blue lines represent the median of the maximum annual mean runoff.The red lines indicate the median for each of the four periods.The scatter (i.e., 75 and 25 % percentiles) indicates the variability between the basins.

Table 1 .
Elevation (m a.s.l.) and size (km 2 ) of the eight Slovak basins.

Table 3 .
The long-term mean monthly values of the precipitation (P) and air temperature (T) for the reference period (1981-2010) and three future time horizons (2011-2040, 2041-2070, 2071-2100) according to the KNMI and MPI climate scenarios: (a) for the river basins with a mean elevation < 435 m a.s.l.; (b) for the river basins with a mean elevation > 435 m a.s.l.

Table 4 .
15TUW model parameters sorted due to the routines and their calibration ranges.

Table 5 .
The optimum HBV model parameters for eight studied basins.