Next Article in Journal
A GPU-Accelerated and LTS-Based Finite Volume Shallow Water Model
Previous Article in Journal
Scaling an Artificial Neural Network-Based Water Quality Index Model from Small to Large Catchments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing Lake Response to Extreme Climate Change Using the Coupled MIKE SHE/MIKE 11 Model: Case Study of Lake Zazari in Greece

1
Greek Biotope/Wetland Centre, The Goulandris Natural History Museum, 57001 Thessaloniki, Greece
2
Department of Hydraulics, Soil Science and Agricultural Engineering, Faculty of Agriculture, Forestry and Natural Environment, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
*
Author to whom correspondence should be addressed.
Water 2022, 14(6), 921; https://doi.org/10.3390/w14060921
Submission received: 17 January 2022 / Revised: 11 March 2022 / Accepted: 12 March 2022 / Published: 15 March 2022
(This article belongs to the Section Hydrology)

Abstract

:
Lakes, either artificial or natural, are greatly important as a component in their catchments’ hydrology, but also as ecosystem service providers. However, due to climate change, they have begun to face numerous problems with their water quality and quantity. Furthermore, general circulation models (GCMs) show future climate change with a reduction in rainfall and increase in temperature. The aim of the current study is to present an application where GCMs and state-of-the-art hydrological modelling system MIKE SHE/MIKE 11 are combined for assessing the response of a Greek lake in terms of its water balance and water level under climate change. Four general circulation models (GCMs; GFDL-CM3, MIROC-ESM-CHEM, MIROC-ESM, IPSL-CM5A-LR) for the extreme climate change scenario of RCP8.5 were used in the basin of Lake Zazari in Greece as a case study. Results showed that, by keeping the irrigated demands (the main water user) unchanged in the future, the lake exhibited a lower water level for all GCMs, fluctuating from −0.70 to −1.8 m for the mean (min) water level and from −0.30 to −1.20 m for the mean (max) water level. Instead of the above and by preserving the amount of withdraw water n from the lake at a certain percentage of inflows, the irrigated area should be reduced from 54.1% to 64.05% depending on the circulation model.

1. Introduction

Natural and artificial lakes are greatly important due to their great contribution to the supply of high-value ecosystem services [1]. The impact of climate change and human activities on their water quality and quantity is thoroughly investigated [2,3,4,5,6,7,8,9]. Among the main reasons is to support management actions related to their ecosystem protection, water supply, flood, and water shortage phenomena [10]. Such actions are expected to face great challenges for setting appropriate measures to preserve or restore lake conditions [11,12,13,14], especially in identified climate change hot spots such as the Mediterranean basin [15] where Greece is located, with more than 50 lakes or reservoirs in its territory.
A solution for developing management strategies is the application of hydrological models with future climate data based on predictions of general circulation models (GCMs). In this way, the water balance components of lakes and their response under extreme climate conditions can be investigated. Depending on the purposes of modelling and the specific attributes of a lake and its catchment, there are different types of hydrological model structures and levels of model complexities that can be selected. For example, there are more complex and data-demanding models, such as MIKE SHE/MIKE 11 [16], WASP [17], CE-QUAL-W2 [18], GSFLOW [19], WATLAC [20], PCLake [21], DYRESM-CAEDYM [22], or simpler ones [6,10,23,24]. The MIKE SHE modelling system is among the most integrated and it was used for analyzing climate change impact on water resources [25], and more specifically on lakes [26,27]. The MIKE SHE model consists of various compartments for different hydrologic analysis purposes. The most integrated combination is its coupling with MIKE 11, which can be used to perform detailed hydrologic and hydraulic analysis at catchment scale applications based on the hydraulic characteristics and dimensions of rivers and lakes [28].
As far as the Greek territory is concerned, there are few studies that have been conducted for wetlands and specifically for natural lakes. There are studies based on either statistical or mathematical models, but these refer to streams that flow into lakes [29] or analyze the impact of the whole hydrological response of the basin into lakes [30,31]. For example, the authors in [32] studied Lake Karla, one of the most interesting lakes in Greece that was drained in 1962 and then restored by the authorities in 1999 due to its environmental importance but also its vital role in irrigated agriculture. The purpose of the study was the water resources management of the basin. In that direction, the authors in [33] studied the role of adaptation measures on agriculture to water resources under climate change, while the authors in [34] studied the socioeconomic impact that a catchment can have with a lake. There is no clear approach that focuses on lakes apart from an effort conducted in [35] to simulate Lakes Volvi and Koronia using hydrological models such as MIKE SHE and MIKE Hydro. However, apart from the Greek territory, scientists around the word are studying the whole response of a catchment that contains a lake or a reservoir in climate change and not the lake itself and the lake’s behavior [36,37,38].
The current study investigates Zazari Lake, situated in northern Greece, due to its environmental importance. There are few studies that focused on the origin of Zazari. The authors in [39] investigated the effect of using annual updated land use maps derived from satellite images instead of a static coarse land use map in a SWAT model for the estimation of sediments and nitrates that load in Zazari Lake. In [40], the authors investigated the internal trends of surface temperature, precipitation, and wind speed over 18 lakes of Greece using ERA5-Land monthly averaged data from 1981–2020. Moreover, the authors in [29] investigated the impacts of climate change on the hydrology of Lakes Cheimaditida and Kerkini by using monthly data and a conceptual model (MIKE BASIN-NAM). The basin of Zazari Lake is part of the wider basin of Cheimaditida Lake. Lastly, the authors in [41] studied the mean monthly water balance of the wider basin of Cheimaditida Lake using the Thorthwaite and Turc methods for actual evapotranspiration estimation.
The aim of this study is to present an application for assessing the response of a lake to climate change using the coupled MIKE SHE/MIKE 11 modelling system in the catchment of Lake Zazari in Greece, thereby creating a fully appropriate system for wetland hydrological applications. Although the coupled MIKE SHE/MIKE 11 model is not new, the novelty of this work is the fact that this coupled MIKE SHE/MIKE 11 model approach has not been used so far for investigating changes in the lake levels due to climate change in a daily time step in Greek lakes. The study was performed for Lake Zazari in Greece considering the four worst climate change predictions of four respective general circulation models (GCMs) for the climate change scenario of RCP8.5. Inflows and outflows from the lake, and its water level fluctuation were used to assess climate change pressure on the lake’s hydrology.

2. Materials and Methods

2.1. Study Area

The study area is the hydrologic basin of Zazari Lake (40°37′30″ N, 21°32′50″ E; northwestern Greece; Figure 1), which covers an area of 98.4 km2 and belongs to the River Basin District of Western Macedonia–Greece. According to the updated world map of the Köppen–Geiger climate classification of [42], the region belongs to the Csa climatic zone (hot summer Mediterranean climate). Annual precipitation during the five-year period of the study ranged from 600 to 800 mm, and mean temperature from 4 °C during winter, and 22.6 °C during summer.
Τhe mountainous area of the basin is mainly composed of low-permeability rocks of a crystalline background, while alluvial deposits cover its lower part. In these quaternary sediments of about 120–130 m thickness [43], two main hydraulically independent aquifer systems were formed, a surface shallow one with aquifer thickness of up to 45 m that grows in coarse-grained surface formations, and a deeper confined aquifer. The shallow aquifer is hydraulically connected with the lake [44,45].
The lake is supplied with water from both its own subcatchment and the Sclithros torrent. The latter, before outflowing into the lake, receives water from the Asprogia and Aetos streams (Figure 1a). Aetos stream, just after it leaves the mountainous area, becomes a drainage ditch of nonpermanent flow.
A sharp crested weir of 17 m width that is located at the outlet of the lake regulates its surface water level. Overspill starts at 599.4 m above sea level (a.s.l.). At this elevation, the lake covers an area of ~1.95 km2 with a maximal depth of 6 m.
The main water user in the area is agriculture. The lake supplies with water the irrigation network of the Limnohori area (220 ha, 1.19 × 106 m3 mean yearly abstraction from May to September), which is located out of its basin, and 60 more hectares of individual plots close to its shore. From the agricultural area that is located close to the settlement of Sclithros, 364 ha are irrigated from Sclithros River, Asprogia stream, and pumping wells. Moreover, 480 ha of agricultural land that is located south of the Aetos and Agrapidies settlements are irrigated by groundwater.

2.2. Hydrological Modeling System

In the context of this study, the coupled MIKE SHE/MIKE 11 modelling system was used. MIKE SHE is a physical, deterministic, fully distributed hydrological model [16,46] that can simulate the major processes of water cycles on both surface and subsurface [47]. Evapotranspiration, infiltration, overland flow, unsaturated and saturated flow, irrigation and drainage, and overland–groundwater exchange were simulated in this study using the corresponding MIKE SHE modules.
MIKE 11 is a 1-dimensional hydraulic modelling system for rivers, lakes, and estuaries that is based on the complete dynamic wave formulation of the Saint Venant equations. Furthermore, it can represent a wide range of hydraulic structures including sluice gates, weirs, gates, bridges, and culverts [28,48].
The two systems, MIKE SHE and MIKE 11, were dynamically coupled [49,50], exchanging information after each computational time step [47,48]. The main streams in the basin and Lake Zazari were included in MIKE 11 for the purposes of this study.

2.3. Data

In the context of this study, MIKE SHE/MIKE 11 was calibrated on the basis of the conditions of the period of 25 July 2012–31 December 2016, and validated on the conditions of the period of 1 January 2017–31 December 2017. The observed data were the variation in lake surface elevation (LSE), with maximal and minimal water levels of 599.98 and 598.56 m, respectively, and the daily temperature and precipitation data from the Limnochori meteorological station (located nearby the lake, 40°37′20.70″ N, 21°33′19.27″ E) for the aforementioned period. Reference evapotranspiration (ETo) was calculated on a daily time step by the equation in [51], taking into consideration the local value from the global map of revised coefficients provided by [52] for achieving equivalent estimations of FAO-56/ASCE reference crop evapotranspiration ETo for short grass (Figure 2) [53,54]. On the basis of ETo, MIKE SHE estimates the actual evapotranspiration from land surface by taking into account the LAI and RD of the land use, the available ponded water, and the soil moisture content in each cell. In the case of flooded areas such as lakes, MIKE SHE assumes that actual evaporation equals to ETo, unless a different approach is adopted [47].
Τhe digital elevation model for the study area was obtained by the EU-DEM v1.1 database with a spatial resolution of 25 m, in which lake bathymetry was incorporated [55].
The land uses of the basin according to the Corine Land Cover 2012 database [56] consist of 32.6% croplands, 46.8% forests, 8.7% pastures/grasslands, 8.3% shrubs/sclerophyllous vegetation, 2% lake, 1.3% artificial surfaces, and 0.3% bare soils (Figure 1b). The general soil properties obtained from the ESDB-ESDAC soil database [57,58] showed that the study area is dominated by loamy sand and sandy loam soils at 36.9% and 44.2% coverage, respectively.
For the analysis of climate change, the predictions of four GCMs {GFDL-CM3 (gf), MIROC-ESM-CHEM (mi), MIROC-ESM (mr), IPSL-CM5A-LR (ip)} for the worst scenario of the highest greenhouse-gas emissions RCP8.5 (mean conditions of 2061–2080) were used. These 4 GCMs were selected among 19 GCMs provided by the WorldClim database [59] because they represented the worst conditions in terms of rainfall reduction and temperature increase compared to the current conditions at the position of the lake. The data are in raster form with spatial resolution of 30 arc-s (1 km2).

2.4. Model Setup and Parameterization

The model area included the whole basin of the lake. Two grid sizes, 400 × 400 m and 200 × 200 m were tested [60] for its discretization. The 400 × 400 m size did not allow for an accurate representation of the river network and lake shape, so it was abandoned before any calibration attempt. Instead, the grid size of 200 × 200 m gave a finer and more accurate representation of the river and lake shape, and of the spatial distribution of physical parameters in the basin. This discretization resulted in 2709 cells (70 rows × 85 columns) and was adopted.
The overall maximal time step was set at 6 h, with maximal time step for the sheet overland flow, 0.5 h; water flow in rivers, 30 s, in the unsaturated zone, 2 h, and 24 h in the saturated zone.
Land use spatial extent was further corrected, taking into account Google Earth images. Leaf area index (LAI) and root depth distribution (RDD) over time were acquired from [61] for each class of CLC 2012.
For the simulation of the irrigation in the basin, eight irrigation command areas were adopted (areas that receive water from the same source). Their spatial extent and water sources are given in Figure 3. In all command areas, irrigation starts at wilting point and ends at field capacity during the period for which LAI and root depth are greater than zero. The irrigation method was sprinklers in all cases.
The hydraulic properties of soil types (Figure 4) were estimated using Van Genuchten equations according to the dataset of [62] (Table 1). The full Richard’s equation was used with the computations to be carried out in all grid sells (soil columns) of the model. For the functional relationship of moisture–retention curve and effective conductivity, the formulas of Van Genuchten were selected using default values for their parameters (Table 1). The unsaturated zone was considered to be uniform to the vertical direction, and its discretization was chosen to vary from 10 cm for the first 0.5 m (5 nodes), 50 cm from 0.5 to 9.0 m (17 nodes) depth, and 200 cm from 9.0 to 45 m depth (18 nodes).
According to [44,45], the saturated zone was simulated by adopting a single surface aquifer seated on an impermeable layer that extended across the catchment. In addition, the aquifer is extended beneath the lake, allowing for a dynamic exchange of water between them. Both the depth of the aquifer and hydraulic conductivities were subjected to calibration. Initial values were set to 5 m depth at the mountainous area, and 40 m at the lower part of the catchment. Following [63], almost 90% of the aquifers under exploitation in Greece consist mainly of sand and gravel with hydraulic conductivities between 10−6 and greater than 10−4 m/s. Since no further information regarding this parameter was found for the study area, the initial value of 5 × 10−5 m/s was adopted for both horizontal (Kx) and vertical (Ky) hydraulic conductivity.
The outer boundary condition at the periphery of the aquifer was set as zero flux, except for a section east of the catchment. For the latter, a permanent hydraulic gradient of −0.0025 was adopted that followed the natural slope of the ground surface. This boundary condition drains ground water downstream out of the catchment.
The dense network of streams in the mountainous areas and the artificial drainage network in the agricultural land were simulated by using the drainage module of MIKE SHE. The drain level was set at −1.5 m below ground surface for the whole catchment, while the time constant (leakage factor) was set to be 1 × 10−7 for the whole catchment except for the valleys, where it was set to be 1 × 10−6. For drainage routing, grid codes were adopted in order for the drain flow to be routed to the nearest river link. The subcatchment of the lake was excluded from the drainage (drain codes set to 0) due to the absence of significant streams.
For the simulation of water flow in Sklithros River, and the Asprogia and Aetos streams (Figure 1a), MIKE 11 was used. Only the sections of the above water bodies located in the lowlands of the catchment were included in the model (Figure 5). In this way, steep slopes of the mountainous branches and accompanied instabilities for the model were avoided. Lake Zazari was simulated as part of Sklithros River by inserting several cross-sections following the shape of its shores (Figure 5). The whole length of the 11 MIKE branches was coupled with MIKE SHE.
Cross-sections for each branch were specified from field investigations. The Manning coefficient for all branches was set to be n = 0.04 [64], except for the lake section, which was set to be n = 0.03. Leakage coefficients were initially set to be 1 × 10−7, but they were subjected to calibration for all branches.
As far as the initial conditions of the model are concerned, two conditions were set at the beginning of the simulation, the first was water depth along the branches that was set to be to 0.01 m, and the second was the initial water depth of the lake, which was set to be equal to the observed water level at the beginning of the simulation.
A small discharge of 0.05 m3/s for the first few days that became zero for the rest of the simulation period was set to be the boundary condition at the upper free ends of the branches. At the outlet of the lake, a Q/h boundary condition was set representing the functioning of the crested weir (Figure 6). Further, the water abstraction from the lake for the supply of the Limnohori collective irrigation network was set as a discharge boundary condition.

2.5. Criteria for Evaluation of Lake Zazari Basin Hydrological Model

Model performance was evaluated by comparing the simulated lake surface elevation (LSE) against its observed measurements, since it was the only parameter that is monitored in the catchment almost at a monthly time step since 2012 [65]. The statistical criteria that were used, the range of their expected values, and values that indicated the best fitting of the model results are described in Table 2 [47,66].

2.6. Climate Change Scenarios

For the analysis of climate change, the predictions of four GCMs {GFDL-CM3 (gf), MIROC-ESM-CHEM (mi), MIROC-ESM (mr), IPSL-CM5A-LR (ip)} for the worst scenario of highest greenhouse-gas emissions RCP8.5 (2061–2080) were used. These 4 GCMs were selected among 19 GCMs provided by the WorldClim database [59] because they represented the worst conditions in terms of precipitation reduction and temperature increase compared to the current conditions at the position of the lake. As the coupled MIKE SHE/MIKE 11 needs daily data for the simulation, the following was assumed in order for the model to be able to simulate in future time: historical data series taken from the weather station of Limnochori were assumed to have the same future pattern. So, the average values of meteorological (temperature) or aggregated (precipitation) parameters in the historical series for each calculated month and the respective values of each month from the Fick and Hijmans dataset were used and compared, and the percentage of increase or decrease of each parameter was calculated. Then, future time series were generated by multiplying each value of the historical dataset with the percentage of increase or decrease that had been calculated in the previous step. After that, the coupled models were again applied using the four new climate datasets that represented the four worst future climate conditions in order to assess the response of the lake.

3. Results and Discussion

3.1. Model Calibration and Validation

The calibration and validation of the model were conducted against observed measurements of the lake surface elevation (LSE) for the period of 25 July 2012–31 December 2016 and 1 January–31 December 2017, respectively.
Simulated vs. observed LSE data for both of the above-mentioned phases are given in Figure 7, which shows that the model satisfactorily followed the seasonal pattern of the lake’s water level fluctuation.
The values of the statistical criteria described in Table 2 were calculated and are given in Table 3, which shows that, in both the calibration and validation phases, ME and MAE were in the order of less than 10 cm when the LSE fluctuated between 598.56 and 599.98 m a.s.l. (in the range of 1.42 m). RMSE and STDres had values close to zero, which is the best value for these criteria. The correlation coefficient (R) and Nash–Sutcliffe model efficiency (R2), also suggested a very good predictive capability of the model, which is also reflected in the validation phase.
Table 4 shows the final values of the parameters that were subjected to calibration. Hence, for the mountainous areas, the depth of the aquifer was 12 m, while for the lowlands, a depth from 15 to 25 m was adopted. Beneath the lake, the depth of the aquifer was 6 m. The final values for hydraulic conductivity after calibration were 1.5 × 10−5 m/s for the mountainous areas, 2 × 10−4 m/s for the lowlands, and 1.5 × 10−5–2 × 10−4 beneath the lake. The leakage coefficient for MIKE 11 branches was in the range of 1 × 10−7 to 1 × 10−6 s−1. Lower values do not allow for the basin to drain, leading to a continuous increase in the water that is stored in the saturated zone, while greater values lead to the continuous drainage of the basin.

3.2. Response of Lake Zazari to Climate Change

For the assessment of Lake Zazari’s response to climate change, the hydrological model of the basin was rerun for each future GCM, and pressures on the lake’s water balance and state were quantified. Data of daily reference evapotranspiration (ETo) and precipitation (P) that had been generated from each GCM were used in the hydrological model. The same initial conditions were adopted as those at the current period. In addition, the irrigated area in the basin and the abstraction of water to meet the needs of the Limnohori irrigation network were kept the same. In this way, four sets of results, one for each GCM, were produced for the future period.
Parameters that were used to quantify the pressure on the lake was precipitation (mm) and evaporation (ETo, mm) from its surface, inflow from Sklithros (mean Qin, m3) and outflow as overspill (mean Qout, m3) as well as the mean minimal (mean (min)) and maximal (mean (max)) LSE (m), and are given in Table 5.
Inflow from Sklithros was estimated using the components of Table 6 following the equation:
Qin = ORR + DR + BR − (IR − PG) − RB
where ORR is the amount of overland runoff to the river (mm), DR is the amount of water from the saturated zone that drains to rivers (mm), BR is the amount of water from the baseflow that flows to the river (mm), IR is the amount needed for irrigation (mm), PG is the pumping amount of water abstracted from the groundwater (mm), and RB is the amount of water that flows from the river to the baseflow (mm).
Table 6 shows the components of the water balance of the area (sub-basins of Sklithros torrent, Asprogia and Aetos streams, 9116 ha) that contributes with runoff to torrent Sklithros and hence to Lake Zazari (Table 5). In Table 6, the irrigation component includes both pumping water from groundwater and extractions from surface water (Sklithros torrent and Asprogia stream).
The values of Table 5 and Table 6 represent mean yearly values of the five-year period of 1 January 2013–31 December 2017 and future conditions of the period of 2061–2080 of each climatic projection. Pressure was evaluated by comparing the values of the above parameters under the current period against future ones.
In Figure 8, the observed values of LSE for the period of 1 January 2013–31 December 2017 are compared against the estimation from the hydrological model values of LSE for the future period and for each GCM {GFDL-CM3 (gf), MIROC-ESM-CHEM (mi), MIROC-ESM (mr), IPSL-CM5A-LR (ip)}. Under all GCMs, the LSE shifted to lower levels, with that of the gf-model being the lowest. The mean (max) water level is expected to drop from −0.30 to −1.20 m, while the mean (min) from −0.70 to −1.80 m, according to each GCM.
Changes in the lake’s level fluctuation are due to changes caused in the components of its water balance, that is, evaporation from and precipitation on the lake’s surface, water inflows and outflows, and direct abstractions.
Evaporation from the lake is increasing, as is indicated by the values of ETo (Table 5). Hence, greater amounts of water (from 18.1% to 9.4%) directly leave the lake.
In all GCMs, mean yearly precipitation substantially dropped, from 544 mm (−25.6%) under the mi-model up to 405 mm (−44.7%) under the gf-model, instead of 731 mm under current climatic conditions. This leads to less direct rainfall on the lake’s surface and reduced amount of water inflows into the lake as runoff from its basin.
Table 5 shows that the mean yearly discharge under the four GCMs from Sklithros River into Zazari drops by more than its half (−54.11% to −64.15%) compared to current climatic conditions. This can be seen from Table 6, which shows the water balance components of the Sklithros river. The lake’s water balance and thereby its level fluctuation are greatly affected, since the irrigated area was considered to be the same.
The reduction in lake level results in smaller quantities of water that overspill downstream. Figure 8 shows that, just after the second year of the future period, the LSE drops below the crest of the weir (599.4 m a.s.l.), and water stops overspilling downstream.
In the above future climatic projections, the area of the irrigation network of Limnohori (220 ha) was considered to be unchanged. Under current climatic conditions, 1.19 × 106 m3 (5400 m3/ha) are directly withdrawn every year from the lake for the irrigation of this area. Further to the above, 60 ha more of individual plots located near the shores of the lake are irrigated directly from the lake (0.324 × 106 m3). Hence, the total amount of water that is yearly abstracted is estimated to be Qir = 1.512 × 106 m3 (Table 5). Rir = 0.26 (or 26%) of the mean yearly inflows mean (Qin) to the lake under current climatic conditions (Table 5 and Table 7). However, under future climatic projections, the above percentage (Rir) increases due to the reduced inflows to the lake, and ranges from 0.73 to 0.57. Keeping the percentage of abstractions (and hence of irrigation pressure) unchanged, and equal to that under current conditions (26%), the available amount of water for irrigation varies from 693,989.5 to 543,581.2 m3 (Table 7). These quantities of water (and under the assumption that irrigation needs per hectare remain the same) allow for the irrigation of 100.66, 121.19, 128.52, and 118.25 ha for each of the future projections instead of 280 ha under current climatic conditions. That means that, in order for irrigation pressure to remain at the same level as it is under current climatic conditions, a reduction in irrigated land from 64.05% to 54.10% (Table 7) following the climatic projections should be considered. In fact, this reduction in irrigated land is expected to be greater, since the rise in temperature leads to higher values of evapotranspiration and hence to greater amounts for irrigation per hectare.
According to the results above, evapotranspiration increases within the basin because temperature is expected to increase, while precipitation and runoff face serious yearly reduction. The above trends can be observed in all of the four climate change models (scenario RCP 8.5) that were used. Similar findings with regard to the trends of the above parameters under different climate change models and scenarios were reported from many other researchers [67,68] for the Eastern Mediterranean region.
Furthermore, it is obvious that Lake Zazari faces a serious future problem as the lake shrinks, and available water volume is decreased. This is something that can be seen intensively in the Mediterranean area [69,70]. Lake Zazari is small compared to its basin (1:50). In addition, more than 25% of its yearly inflow is withdrawn during summer for irrigation purposes. This study showed that a reduction in rainfall of more than 25% in the basin, which is expected under severe climate change scenarios, leads to a substantial decrease in LSE. Lakes located at similar latitudes or climatic zones, especially those with similar hydrology and exploitation patterns, are expected to be subjected to analogous pressure due to climate change. As it regards Greece, there are at least 8 lakes located in the northern part of the country that have similar climatic conditions and exploitation patterns. Irrigation is one of the main aspects in the basin that would be affected, in accordance to other studies as well [71,72,73]. The results of this study can offer a quantifiable framework to the competent authorities for the design of adaptation measures related to irrigated agriculture for mitigating the impact of climate change on lake surface elevation (LSE).

4. Conclusions

This study presented an application for assessing the response of a lake to climate change using the coupled MIKE SHE/MIKE 11 modelling system in the catchment of Lake Zazari in Greece, thereby creating a fully appropriate system for wetland hydrological applications. Overall analysis and the provided results of the Lake Zazari hydrological model showed that MIKE SHE/MIKE 11 accurately simulated hydrology in the basin where both surface and ground water not only coexist, but also interact. The modelling system was also able to incorporate hydraulic structures and water management practices, and simulate irrigation at the basin scale through the combined use of surface and groundwater resources.
The calibrated model provided a detailed evaluation of the lake’s water level and its water balance components, taking into consideration the climatic conditions of the current studied period (1 January 2012–31 December 2017), but also of future climate predictions.
The four climate change GCMs showed that precipitation dropped by almost half, while temperature and thereby evaporation increased. Under these climatic conditions, runoff and hence inflow discharge into the lake decrease.
In case that irrigation demands remain the same, the lake faces diminished water level. In addition, fewer quantities of water overspill from the lake, in turn affecting the hydrology in the downstream basin.
On the other hand, any attempt to maintain the abstracted amounts of water for irrigation at the same level as it is under the current climatic conditions would lead to a reduction in irrigated area of more than half.

Author Contributions

Conceptualization, D.P. (Dimitris Papadimos) and K.D.; methodology, D.P. (Dimitris Papadimos) and K.D.; software, D.P. (Dimitris Papadimos) and K.D.; validation, D.P. (Dimitris Papadimos) and K.D.; formal analysis, D.P. (Dimitris Papadimos) and K.D.; investigation, D.P. (Dimitris Papadimos) and K.D.; resources, D.P. (Dimitris Papadimos) and K.D.; data curation, D.P. (Dimitris Papadimos) and K.D.; writing—original draft preparation, D.P. (Dimitris Papadimos), K.D. and D.P. (Dimitris Papamichail); writing—review and editing, D.P. (Dimitris Papadimos), K.D. and D.P. (Dimitris Papamichail); visualization, D.P. (Dimitris Papadimos) and K.D.; supervision, D.P. (Dimitris Papamichail). All authors have read and agreed to the published version of the manuscript.

Funding

This research was financed by the European Union Cohesion Fund (Partnership Agreement 2014–2020) Act MIS 5001204. Data used in this research came from Act MIS 371140 funded by the European Regional Development Fund (National Strategic Reference Framework 2007–2013).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The present study was based on freely available data provided by seven international databases that are fully described in Section 2: https://land.copernicus.eu/imagery-in-situ/eu-dem/eu-dem-v1.1 (accessed on 1 October 2021), https://repository.biodiversity-info.gr/handle/11340/1935 (accessed on 1 October 2021), https://land.copernicus.eu/pan-european/corine-land-cover/clc-2012 (accessed on 1 October 2021), https://esdac.jrc.ec.europa.eu/content/european-soil-database-derived-data (accessed on 1 October 2021), https://doi.pangaea.de/10.1594/PANGAEA.868808 (accessed on 1 October 2021), https://www.worldclim.org/ (accessed on 1 October 2021), https://earthexplorer.usgs.gov/ (accessed on 1 October 2021).

Acknowledgments

The present study was conducted in the Frame of the Greek National Water Monitoring Network according to JMD 140384/2011, implemented by The Goulandris Natural History Museum, Greek Biotope/Wetland Centre (EKBY). The network is supervised by the General Directorate for Waters of the Ministry of Environment and Energy. We thank the reviewers and the editor for their comments and recommendations that improved our article.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Wong, C.P.; Jiang, B.; Bohn, T.J.; Lee, K.N.; Lettenmaier, D.P.; Ma, D.; Ouyang, Z. Lake and wetland ecosystem services measuring water storage and local climate regulation. Water Resour. Res. 2017, 53, 3197–3223. [Google Scholar] [CrossRef] [Green Version]
  2. Mooij, W.M.; Hülsmann, S.; De Senerpont Domis, L.N.; Nolet, B.A.; Bodelier, P.L.E.; Boers, P.C.M.; Dionisio Pires, L.M.; Gons, H.J.; Ibelings, B.W.; Noordhuis, R.; et al. The impact of climate change on lakes in the Netherlands: A review. Aquat. Ecol. 2005, 39, 381–400. [Google Scholar] [CrossRef]
  3. Altunkaynak, A. Forecasting Surface Water Level Fluctuations of Lake Van by Artificial Neural Networks. Water Resour. Manag. 2006, 21, 399–408. [Google Scholar] [CrossRef]
  4. Moss, B. Cogs in the endless machine: Lakes, climate change and nutrient cycles: A review. Sci. Total Environ. 2012, 434, 130–142. [Google Scholar] [CrossRef] [PubMed]
  5. Demertzi, K.A.; Papamichail, D.M.; Georgiou, P.E.; Karamouzis, D.N.; Aschonitis, V.G. Assessment of rural and highly seasonal tourist activity plus drought effects on reservoir operation in a semi-arid region of Greece using the WEAP model. Water Int. 2014, 39, 23–34. [Google Scholar] [CrossRef]
  6. Demertzi, K.; Papadimos, D.; Aschonitis, V.; Papamichail, D. A Simplistic Approach for Assessing Hydroclimatic Vulnerability of Lakes and Reservoirs with Regulated Superficial Outflow. Hydrology 2019, 6, 61. [Google Scholar] [CrossRef] [Green Version]
  7. Abdel-Fattah, S.; Krantzberg, G. A review: Building the resilience of Great Lakes beneficial uses to climate change. Sustain. Water Qual. Ecol. 2014, 3, 3–13. [Google Scholar] [CrossRef]
  8. Zhang, C.; Lai, S.; Gao, X.; Liu, H. A review of the potential impacts of climate change on water environment in lakes and reservoirs. J. Lake Sci. 2016, 28, 691–700. [Google Scholar] [CrossRef] [Green Version]
  9. Gaglio, M.; Aschonitis, V.; Pieretti, L.; Santos, L.; Gissi, E.; Castaldelli, G.; Fano, E.A. Modelling past, present and future Ecosystem Services supply in a protected floodplain under land use and climate changes. Ecol. Model. 2019, 403, 23–34. [Google Scholar] [CrossRef]
  10. Dessie, M.; Verhoest, N.E.C.; Pauwels, V.R.N.; Adgo, E.; Deckers, J.; Poesen, J.; Nyssen, J. Water balance of a lake with floodplain buffering: Lake Tana, Blue Nile Basin, Ethiopia. J. Hydrol. 2015, 522, 174–186. [Google Scholar] [CrossRef]
  11. Cooke, G.D.; Welch, E.B.; Peterson, S.A.; Nicholson, S.A. Restoration and Management of Lakes and Reservoirs; CRC Press: Boca Raton, FL, USA, 2005; p. 616. [Google Scholar]
  12. Doulgeris, C.; Georgiou, P.; Papadimos, D.; Papamichail, D. Ecosystem approach to water resources management using the MIKE 11 modeling system in the Strymonas River and Lake Kerkini. J. Environ. Manag. 2012, 94, 132–143. [Google Scholar] [CrossRef]
  13. Cookey, P.E.; Darnsawasdi, R.; Ratanachai, C. A conceptual framework for assessment of governance performance of lake basins: Towards transformation to adaptive and integrative governance. Hydrology 2016, 3, 12. [Google Scholar] [CrossRef] [Green Version]
  14. Jeppesen, E.; Søndergaard, M.; Liu, Z. Lake restoration and management in a climate change perspective: An introduction. Water-SUI 2017, 9, 122. [Google Scholar] [CrossRef] [Green Version]
  15. Loizidou, M.; Giannakopoulos, C.; Bindi, M.; Moustakas, K. Climate change impacts and adaptation options in the Mediterranean basin. Reg. Environ. Chang. 2016, 16, 1859–1861. [Google Scholar] [CrossRef] [Green Version]
  16. Abbott, M.B.; Bathurst, J.C.; Cunge, J.A.; O’Connell, P.E.; Rasmussen, J. An introduction to the European Hydrological System—Systeme Hydrologique Europeen, “SHE”, 1: History and philosophy of a physically-based, distributed modelling system. J. Hydrol. 1986, 87, 45–59. [Google Scholar] [CrossRef]
  17. Ambrose, R.B.; Wool, T.A.; Matrin, J.L. Part A: Model Documentation. In The Water Quality Analysis Simulation Program, WASP5, Version 5.10; U.S. Environmental Protection Agency: Athens, GA, USA, 1993; p. 251. [Google Scholar]
  18. Cole, T.M.; Buchak, E.M. Draft User Manual, Instruction Report EL-95-1. In CE-QUAL-W2: A Two-Dimensional, Laterally Averaged Hydrodynamic and Water Quality Model, Version 2.0; Corps of Engineers Waterways Experiment Station Vicksburg: Vicksburg, MI, USA, 1995; p. 357. [Google Scholar]
  19. Markstrom, S.L.; Niswonger, R.G.; Regan, R.S.; Prudic, D.E.; Barlow, P.M. GSFLOW-Coupled Ground-Water and Surface-Water FLOW Model Based on the Integration of the Precipitation-Runoff Modeling System (PRMS) and the Modular Ground-Water Flow Model (MODFLOW-2005); U.S. Geological Survey Techniques and Methods 6-D1: Reston, VA, USA, 2008; p. 240. [Google Scholar]
  20. Zhang, Q. Development and application of an integrated hydrological model for lake watersheds. Proc. Environ. Sci. 2011, 10, 1630–1636. [Google Scholar] [CrossRef] [Green Version]
  21. Hu, F.; Bolding, K.; Bruggeman, J.; Jeppesen, E.; Flindt, M.R.; van Gerven, L.; Janse, J.H.; Janssen, A.B.G.; Kuiper, J.J.; Mooij, W.M.; et al. FABM-PCLake—Linking aquatic ecology with hydrodynamics. Geosci. Model. Dev. 2016, 9, 2271–2278. [Google Scholar] [CrossRef] [Green Version]
  22. Luo, L.; Hamilton, D.; Lan, J.; McBride, C.; Trolle, D. Autocalibration of a one-dimensional hydrodynamic-ecological model (DYRESM 4.0-CAEDYM 3.1) using a Monte Carlo approach: Simulations of hypoxic events in a polymictic lake. Geosci. Model. Dev. 2018, 11, 903–913. [Google Scholar] [CrossRef] [Green Version]
  23. Bracht-Flyr, B.; Istanbulluoglu, E.; Fritz, S. A hydro-climatological lake classification model and its evaluation using global data. J. Hydrol. 2013, 486, 376–383. [Google Scholar] [CrossRef] [Green Version]
  24. Yang, K.; Lu, H.; Yue, S.; Zhang, G.; Lei, Y.; La, Z.; Wang, W. Quantifying recent precipitation change and predicting lake expansion in the Inner Tibetan Plateau. Clim. Chang. 2018, 147, 149–163. [Google Scholar] [CrossRef]
  25. Keilholz, P.; Disse, M.; Halik, Ü. Effects of Land Use and Climate Change on Groundwater and Ecosystems at the Middle Reaches of the Tarim River Using the MIKE SHE Integrated Hydrological Model. Water-SUI 2015, 7, 3040–3056. [Google Scholar] [CrossRef] [Green Version]
  26. Dimitriou, E.; Moussoulis, E. Hydrological and nitrogen distributed catchment modeling to assess the impact of future climate change at Trichonis Lake, western Greece. Hydrogeol. J. 2010, 18, 441–454. [Google Scholar] [CrossRef]
  27. Singh, C.R.; Thompson, J.R.; French, J.R.; Kingston, D.G.; MacKay, A.W. Modelling the impact of prescribed global warming on runoff from headwater catchments of the Irrawaddy River and their implications for the water level regime of Loktak Lake, northeast India. Hydrol. Earth Syst. Sci. 2010, 14, 1745–1765. [Google Scholar] [CrossRef]
  28. Thompson, J.R.; Sørensen, H.R.; Gavin, H.; Refsgaard, A. Application of the coupled MIKE SHE/MIKE 11 modelling system to a lowland wet grassland in southeast England. J. Hydrol. 2004, 293, 151–179. [Google Scholar] [CrossRef]
  29. Doulgeris, C.; Papadimos, D.; Kapsomenakis, J. Impacts of climate change on the hydrology of two Natura 2000 sites in Northern Greece. Reg. Environ. Chang. 2016, 16, 1941–1950. [Google Scholar] [CrossRef]
  30. Antonopoulos, V.Z.; Gianniou, S.K.; Antonopoulos, A.V. Artificial neural networks and empirical equations to estimate daily evaporation: Application to Lake Vegoritis, Greece. Hydrol. Sci. J. 2016, 61, 2590–2599. [Google Scholar] [CrossRef] [Green Version]
  31. Sidiropoulos, P.; Tziatzios, G.; Vasiliades, L.; Mylopoulos, N.; Loukas, A. Groundwater Nitrate Contamination Integrated Modeling for Climate and Water Resources Scenarios: The Case of Lake Karla Over-Exploited Aquifer. Water 2019, 11, 1201. [Google Scholar] [CrossRef] [Green Version]
  32. Panagopoulos, Y.; Dimitriou, E. A Large-Scale Nature-Based Solution in Agriculture for Sustainable Water Management: The Lake Karla Case. Sustainability 2020, 12, 6761. [Google Scholar] [CrossRef]
  33. Kourgialas, N.N. A critical review of water resources in Greece: The key role of agricultural adaptation to climate-water effects. Sci. Total Environ. 2021, 775, 145857. [Google Scholar] [CrossRef]
  34. Zafeiriou, E.; Andrea, V.; Tampakis, S.; Karanikola, P. Wetlands Management in Northern Greece: An Empirical Survey. Water 2020, 12, 3181. [Google Scholar] [CrossRef]
  35. Malamataris, D.; Kolokytha, E.; Loukas, A. Integrated hydrological modelling of surface water and groundwater under climate change:the case of the Mygdonia basin in Greece. J. Water Clim. Chang. 2020, 11, 1429–1454. [Google Scholar] [CrossRef]
  36. Althoff, D.; Rodrigues, L.N.; Da Silva, D.D. Impacts of climate change on the evaporation and availability of water in small reservoirs in the Brazilian savannah. Clim. Chang. 2020, 159, 215–232. [Google Scholar] [CrossRef]
  37. Xu, D.; Lyon, S.W.; Mao, J.; Dai, H.; Jarsjö, J. Impacts of multi-purpose reservoir construction, land-use change and climate change on runoff characteristics in the Poyang Lake basin, China. J. Hydrol. Reg. Stud. 2020, 29, 100694. [Google Scholar] [CrossRef]
  38. Gardner, E.; Burningham, H.; Thompson, J.R. Impacts of climate change and hydrological management on a coastal lake and wetland system. Irish Geogr. 2019, 1, 52. [Google Scholar] [CrossRef]
  39. Samarinas, N.; Tziolas, N.; Zalidis, G. Improved Estimations of Nitrate and Sediment Concentrations Based on SWAT Simulations and Annual Updated Land Cover Products from a Deep Learning Classification Algorithm. ISPRS Int. J. Geo-Inf. 2020, 9, 576. [Google Scholar] [CrossRef]
  40. Stefanidis, K.; Varlas, G.; Papadopoulos, A.; Dimitriou, E. Four Decades of Surface Temperature, Precipitation, and Wind Speed Trends over Lakes of Greece. Sustainability 2021, 13, 9908. [Google Scholar] [CrossRef]
  41. Gudulas, K.; Voudouris, K.; Soulios, G.; Dimopoulos, G. Comparison of different methods to stimate actual evapotranspiration and hydrologic balance. Desalination Water Treat. 2013, 51, 13–15. [Google Scholar] [CrossRef]
  42. Peel, M.C.; Finlayson, B.L.; McMahon, T.A. Updated world map of the Köppen-Geiger climate classification. Hydrol. Earth Syst. Sci. 2007, 11, 1633–1644. [Google Scholar] [CrossRef] [Green Version]
  43. Soulios, G.; Dimopoulos, G.; Mountrakis, D.; Psilovikos, A.; Pennas, P.; Chatzidimitriadis, E.; Vafiadis, P. Research on the hydrological balance of basins in Greece (Example from the basin of Sklithros, Florina). B Geol. Soc. Greece 1989, XXVI, 129–153. [Google Scholar]
  44. Gudulas, K. Investigation of the Operating Mechanisms of the Lakes Cheimaditida and Zazari in the Amyntaio basin, of Florina County, from a Hydrogeological and Environmental Point of View. Ph.D. Thesis, Aristotle University, Thessaloniki, Greece, 2012. (In Greek). [Google Scholar]
  45. Pashos, P.; Nikolaou, N.; Stamos, A. Report on technical and geological identification of landslides in the local settlements of Valtonera, Pedinou, Rodona and Fanou in municipality of Amynteo. In Hellenic Survey of Geology and Mineral Exploration; Directorate of Technical Geology, Regional Unit of Epirus, Regional Unit of Dytiko: Macedonia, Greece, 2018. (In Greek) [Google Scholar]
  46. Graham, D.N.; Butts, M.B. Flexible, integrated watershed modelling with MIKE SHE. In Watershed Models; Singh, V.P., Frevert, D.K., Eds.; CRC Press: Boca Raton, FL, USA, 2005; pp. 245–272. ISBN 0849336090. [Google Scholar]
  47. DHI. MIKE SHE Water Movement: User Guide, Volume 1; Danish Hydraulic Institute: Hørsholm, Denmark, 2017; p. 420. [Google Scholar]
  48. DHI. MIKE 11 A Modelling System for Rivers and Channels User Guide; Danish Hydraulic Institute: Hørsholm, Denmark, 2017; p. 510. [Google Scholar]
  49. Sørensen, H.R.; Klucovska, J.; Topolska, J.; Clausen, T.; Refsgaard, J.C. An engineering case study—Modelling the influences of Gabcikovo hydropower plant on the hydrology and ecology in the Slovakian part of the river branch system of ZitnyOstrov. In Distributed Hydrological Modelling; Refsgaard, J.C., Ed.; Kluwer: Dordrecht, The Netherlands, 1996; pp. 233–253. [Google Scholar]
  50. Refsgaard, J.C.; Sørensen, H.R. Water management of the Gabcikovo Scheme for balancing the interest of hydropower and environment. In Operational Water Management; Refsgaard, J.C., Karalis, E.A., Eds.; Balkema: Rotterdam, The Netherlands, 1997; pp. 365–372. [Google Scholar]
  51. Hargreaves, G.H.; Samani, Z.A. Estimating potential evapotranspiration. J. Irrig. Drain. 1982, 108, 225–230. [Google Scholar] [CrossRef]
  52. Aschonitis, V.G.; Papamichail, D.; Demertzi, K.; Colombani, N.; Mastrocicco, M.; Ghirardini, A.; Castaldelli, G.; Fano, E. High-resolution global grids of revised Priestley-Taylor and Hargreaves-Samani coefficients for assessing ASCE-standardized reference crop evapotranspiration and solar radiation. Earth Syst. Sci. Data 2017, 9, 615–638. [Google Scholar] [CrossRef] [Green Version]
  53. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements; Irrigation and Drainage Paper 56; Food and Agriculture Organization of the United Nations: Rome, Italy, 1998. [Google Scholar]
  54. Allen, R.G.; Walter, I.A.; Elliott, R.; Howell, T.; Itenfisu, D.; Jensen, M. The ASCE Standardized Reference Evapotranspiration Equation; Allen, R.G., Walter, I.A., Elliott, R., Howell, T., Itenfisu, D., Jensen, M., Eds.; Final Report (ASCE-EWRI); Environmental and Water Resources Institute, Task Committee on Standardization of Reference Evapotranspiration of the Environmental and Water Resources Institute: Reston, VA, USA, 2005. [Google Scholar]
  55. Apostolakis, A.; Papadimos, D.; Zervas, D. Bathymetry Map of Zazari Lake. Greek Biotope/Wetland Centre (ΕΚΒΥ), 2015, Scale 1:12.000. Available online: http://hdl.handle.net/11340/1935 (accessed on 1 October 2021).
  56. Corine Land Cover 2012. Available online: https://land.copernicus.eu/pan-european/corine-land-cover/clc-2012 (accessed on 1 October 2021).
  57. Hiederer, R. Mapping Soil Properties for Europe—Spatial Representation of Soil Database Attributes. In EUR26082EN Scientific and Technical Research Series; Publications Office of the European Union: Luxembourg, 2013; pp. 2013–2047. [Google Scholar] [CrossRef]
  58. Hiederer, R. Mapping Soil Typologies—Spatial Decision Support Applied to European Soil Database. In EUR25932EN Scientific and Technical Research Series; Publications Office of the European Union: Luxembourg, 2013; pp. 2013–2147. ISSN 1831-9424. [Google Scholar] [CrossRef]
  59. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  60. Vazquez, R.F.; Feyen, L.; Feyen, J.; Refsgaard, J.C. Effect of grid size on effective parameters and model performance of the MIKE-SHE code. Hydrol. Proces. 2002, 16, 355–372. [Google Scholar] [CrossRef]
  61. Ministry of Development. Development of Water Resources Management Systems and Tools for the Water Districts of Makedonia and Thraki; Special Issue I; ENM Consulting Engineers: Athens, Greece, 2006; p. 55. (In Greek) [Google Scholar]
  62. Carsel, R.F.; Parrish, R.S. Developing joint probability distributions of soil water retention characteristics. Water Resour. Manag. 1988, 24, 755–769. [Google Scholar] [CrossRef] [Green Version]
  63. Kallergis, G. Hydrology of geological formations, groundwater exploration, development and management, hydrothermal-thermal waters, mathematical models. In Applied Environmental Hydrogeology, 2nd ed.; Technical Chamber of Greece: Athens, Greece, 2001; p. 432. (In Greek) [Google Scholar]
  64. Chow, V.T. Open-Channel Hydraulics; International Student Edition; McGraw-Hill: New York, NY, USA, 1959; p. 680. [Google Scholar]
  65. Greek Biotope/Wetland Center, Water Level Data of Zazari Lake. 2018. Available online: http://www.biodiversity-info.gr/index.php/el/lakes-data#!Oze20170904_IMG_4511 (accessed on 1 October 2021).
  66. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  67. Coppens, J.; Trolle, D.; Jeppesen, E.; Beklioglu, M. The impact of climate change on a Mediterranean shallow lake: Insights based on catchment and lake modelling. Reg. Environ. Chang. 2020, 20, 62. [Google Scholar] [CrossRef]
  68. Bucak, T.; Trolle, D.; Andersen, H.E.; Thodsen, H.; Erdogan, S.; Levi, E.E.; Filiz, N.; Jeppesen, E.; Beklioglu, M. Future water availability in the largest freshwater Mediterranean lake is at great risk as evidenced from simulations with the SWAT model. Sci. Total Environ. 2017, 581, 413–425. [Google Scholar] [CrossRef] [PubMed]
  69. Gorguner, M.; Kavvas, M.L. Modeling impacts of future climate change on reservoir storages and irrigation water demands in a Mediterranean basin. Sci. Total Environ. 2020, 748, 141246. [Google Scholar] [CrossRef] [PubMed]
  70. Rocha, J.; Carvalho-Santos, C.; Diogo, P.; Beca, P.; Keizer, J.J.; Nunes, J.P. Impacts of climate change on reservoir water availability, quality and irrigation needs in a water scarce Mediterranean region (southern Portugal). Sci. Total Environ. 2020, 736, 139477. [Google Scholar] [CrossRef]
  71. Ahmadaali, J.; Barani, G.-A.; Qaderi, K.; Hessari, B. Analysis of the Effects of Water Management Strategies and Climate Change on the Environmental and Agricultural Sustainability of Urmia Lake Basin, Iran. Water 2018, 10, 160. [Google Scholar] [CrossRef] [Green Version]
  72. Woolway, R.I.; Kraemer, B.M.; Lenters, J.D.; Merchant, C.; O’Reilly, C.M.; Sharma, S. Global lake responses to climate change. Nat. Rev. Earth Environ. 2020, 1, 388–403. [Google Scholar] [CrossRef]
  73. Van Der Schriek, T.; Giannakopoulos, C.; Varotsos, K.V. The impact of future climate change on bean cultivation in the Prespa Lake catchment, northern Greece. Euro-Mediterr. J. Environ. Integr. 2020, 5, 14. [Google Scholar] [CrossRef]
Figure 1. (a) Study site of Zazari basin and (b) land uses of the basin.
Figure 1. (a) Study site of Zazari basin and (b) land uses of the basin.
Water 14 00921 g001
Figure 2. Monthly precipitation, reference evapotranspiration, and temperature of Zazari basin for the period from July 2012 to December 2017.
Figure 2. Monthly precipitation, reference evapotranspiration, and temperature of Zazari basin for the period from July 2012 to December 2017.
Water 14 00921 g002
Figure 3. Spatial extent of irrigation command areas and their water source in the hydrological model of Lake Zazari basin.
Figure 3. Spatial extent of irrigation command areas and their water source in the hydrological model of Lake Zazari basin.
Water 14 00921 g003
Figure 4. Spatial extent of the soil types in the hydrological model of Lake Zazari basin.
Figure 4. Spatial extent of the soil types in the hydrological model of Lake Zazari basin.
Water 14 00921 g004
Figure 5. MIKE 11 river network and cross-sections that were included in the hydrological model of Lake Zazari basin.
Figure 5. MIKE 11 river network and cross-sections that were included in the hydrological model of Lake Zazari basin.
Water 14 00921 g005
Figure 6. Q/LSE relation of crested weir at the outlet of Zazari Lake.
Figure 6. Q/LSE relation of crested weir at the outlet of Zazari Lake.
Water 14 00921 g006
Figure 7. Observed and simulated surface elevation (LSE, m a.s.l.) of Lake Zazari for the calibration period (25 July 2012–31 December 2016) and the validation period (1 January–31 December 2017).
Figure 7. Observed and simulated surface elevation (LSE, m a.s.l.) of Lake Zazari for the calibration period (25 July 2012–31 December 2016) and the validation period (1 January–31 December 2017).
Water 14 00921 g007
Figure 8. Water level of Lake Zazari 1 January 2013–31 December 2017 and future climate change conditions.
Figure 8. Water level of Lake Zazari 1 January 2013–31 December 2017 and future climate change conditions.
Water 14 00921 g008
Table 1. Soil type, saturated moisture content (θs, cm3/cm3), residual moisture content (θr, cm3/cm3), suction at field capacity (pFfc), suction at wilting point (pFw), saturated hydraulic conductivity (Ks, m/s), empirical constants in Van Genuchten formulas a (cm−1) and n (-), used in the hydrological model of Lake Zazari basin.
Table 1. Soil type, saturated moisture content (θs, cm3/cm3), residual moisture content (θr, cm3/cm3), suction at field capacity (pFfc), suction at wilting point (pFw), saturated hydraulic conductivity (Ks, m/s), empirical constants in Van Genuchten formulas a (cm−1) and n (-), used in the hydrological model of Lake Zazari basin.
ParameterClayLoamSiltLoamy SandSandy Clay LoamSandy Loam
θs0.380.430.460.410.390.41
θr0.0680.0780.0340.0570.10.065
pFfc 12.82.22.41.92.52
pFw 14.24.24.24.24.24.2
Ks5.55 × 10−72.89 × 10−66.944 × 10−74.05 × 10−53.64 × 10−61.23 × 10−5
a0.0080.0360.0160.1240.0590.075
n1.091.561.372.281.481.89
1 pF: log cm of H2O column.
Table 2. Statistical criteria for evaluating the performance of the Lake Zazari basin hydrological model.
Table 2. Statistical criteria for evaluating the performance of the Lake Zazari basin hydrological model.
CriterionEquationRangeBest Value
Mean error M E j = t O b s i t S i m i t n −∞ ~ +∞0
Mean absolute error M A E j = t O b s i , t S i m i , t n 0 ~ +∞0
Root mean square error R M S E j = t O b s i , t S i m i , t 2 n 0 ~ +∞0
Standard deviation of residuals S T D r e s j = t O b s i , t S i m i , t t O b s i , t S i m i , t n 2 n 0 ~ +∞0
Correlation coefficient R j = t S i m i , t S i m i , t ¯ · O b s i , t O b s i , t ¯ t S i m i , t S i m i , t ¯ 2 · t O b s i , t O b s i , t ¯ 2 0 ~ 11
Nash–Sutcliffe correlation coefficient R j 2 = 1 t O b s i , t S i m i , t 2 t O b s i , t O b s i , t ¯ 2 −∞ ~ 1.01
Obs: observed values at a location i and at a time t, Sim: simulated values at a location i and at a time t, n: number of observations, O b s i , t ¯ : means of observations at location i, S i m i , t ¯ : means of calculations at location i.
Table 3. Values of statistical criteria used for the evaluation of the Lake Zazari basin hydrological model performance.
Table 3. Values of statistical criteria used for the evaluation of the Lake Zazari basin hydrological model performance.
MEMAERMSESTDresRR2
Calibration
25 July 2012–31 December 2016
0.01930.07190.09420.09220.97200.9331
Validation
1 January–31 December 2017
0.07780.09820.12480.09750.98990.8733
Table 4. Final values of calibrated parameters in the coupled MIKE SHE/MIKE 11 hydrological model of Lake Zazari basin.
Table 4. Final values of calibrated parameters in the coupled MIKE SHE/MIKE 11 hydrological model of Lake Zazari basin.
ParametersMountainous AreaLowlandsLake
MIKE SHE
Aquifer depth (m)1215–256
Hydraulic conductivty, Kx, Ky (m/s)1.5 × 10−52 × 10−41.5 × 10−5–2 × 10−4
MIKE 11
Leakage coefficient (s−1) 1 × 10−7–1 × 10−6
Table 5. Comparison of Lake Zazari water balance components assessed by MIKE SHE/MIKE 11 under current climatic conditions and future climatic projections.
Table 5. Comparison of Lake Zazari water balance components assessed by MIKE SHE/MIKE 11 under current climatic conditions and future climatic projections.
Current ConditionsFuture Predictions
Observedgfipmimr
mmmmΔ (%)mmΔ (%)mmΔ (%)mmΔ (%)
Precipitation731405−44.7475−35.1544−25.6520−28.9
ETo961113618.110529.4109513.9111716.2
DischargeQ (m3)Q (m3)Δ (%)Q (m3)Δ (%)Q (m3)Δ (%)Q (m3)Δ (%)
Mean (Qin)5.73 × 1062.06 × 106−64.152.48 × 106−56.702.63 × 106−54.112.42 × 106−57.75
Mean (Qout)4.02 × 1066.81 × 104−98.302.60 × 105−93.542.99 × 105−92.552.11 × 105−94.76
Water levelLSE (m)LSE (m)Δ (%)LSE (m)Δ (%)LSE (m)Δ (%)LSE (m)Δ (%)
Mean (max)599.70598.50−0.20599.20−0.08599.40−0.05599.10−0.10
Mean (min)598.70596.90−0.30597.80−0.15598.00−0.12597.70−0.17
Fluctuation1.001.60−60.001.40−40.001.40−40.001.40−40.00
Table 6. Water balance components of sub-basins (9116 ha) upstream the lake Zazari under current (obs) (1 January 2013–3 December 2017) and future (2061–2080) climatic conditions.
Table 6. Water balance components of sub-basins (9116 ha) upstream the lake Zazari under current (obs) (1 January 2013–3 December 2017) and future (2061–2080) climatic conditions.
Water Balance ComponentsFuture Climatic Projections
obs (mm)gf (mm)ip (mm)mi (mm)mr (mm)
Precipitation731.3404.7474.6544.5520.1
Evapotranspiration593.2432.0460.3513.3510.6
Overland runoff to river (ORR)2.70.20.30.70.3
Irrigation (IR)31.340.740.341.340.3
Pumping from groundwater (PG)29.538.938.539.138.4
Saturated zone drain to rivers (DR)44.29.812.012.811.7
Baseflow to river (BR)17.814.516.617.316.4
Baseflow from river (RB)0.10.10.10.10.1
Sklithros River outflow, mm6323272826
Sklithros River outflow, m35.73 × 1062.06 × 1062.48 × 1062.59 × 1062.41 × 106
Table 7. Pressure analysis of irrigation under current climatic conditions and future climate projections.
Table 7. Pressure analysis of irrigation under current climatic conditions and future climate projections.
Current ConditionsFuture Climatic Projections
Observedgfipmimr
Qir (m3)1,512,000.01,512,000.01,512,000.01,512,000.01,512,000.0
Mean (Qin) (m3)5,730,000.02,060,000.02,480,000.02,630,000.02,420,000.0
Rir = Qir/Qin0.260.730.610.570.62
Qir under stable
Rir = 26%
1,512,000.0543,581.2654,408.4693,989.5638,575.9
Irrigated area (ha)280.00100.66121.19128.52118.25
Reduced irrigated area (%) 64.0556.7254.1057.77
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Papadimos, D.; Demertzi, K.; Papamichail, D. Assessing Lake Response to Extreme Climate Change Using the Coupled MIKE SHE/MIKE 11 Model: Case Study of Lake Zazari in Greece. Water 2022, 14, 921. https://doi.org/10.3390/w14060921

AMA Style

Papadimos D, Demertzi K, Papamichail D. Assessing Lake Response to Extreme Climate Change Using the Coupled MIKE SHE/MIKE 11 Model: Case Study of Lake Zazari in Greece. Water. 2022; 14(6):921. https://doi.org/10.3390/w14060921

Chicago/Turabian Style

Papadimos, Dimitris, Kleoniki Demertzi, and Dimitris Papamichail. 2022. "Assessing Lake Response to Extreme Climate Change Using the Coupled MIKE SHE/MIKE 11 Model: Case Study of Lake Zazari in Greece" Water 14, no. 6: 921. https://doi.org/10.3390/w14060921

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