Changes in water temperature in lakes and reservoirs is one of the key processes affecting the spatial and temporal distribution of chemical compounds and biological elements such as plankton and fish [Elci 2008; Genova et al. 2010; Yu et al. 2010]. The water temperature depends on geographic conditions (e.g., bathymetry, elevation, latitude and interactions with groundwater) and seasonal changes in meteorological conditions (e.g., solar radiation, cloud cover, wind speed and direction) [Lee et al. 2013; Tuan et al. 2009; Tarasiuk et al. 2015]. The temperaturebiochemistry interdependences are often bidirectional, and the concentration of suspended solids or plankton inhibiting the penetration of solar radiation, for example, may impact the water temperature [Kumagai et al. 2000; Mazumder et al. 1990]. Temperature affects the chemical and biological status directly, as particular processes occur more intensively at a given range of temperature and indirectly influence the water density, and thus, water mixing or stratification [Elci 2008]. In addition to the effects of temperature, water mixing processes also depend on the inflow and outflow currents, wind and rainfall. The strength of all these factors depends regionally on the climate and geographical conditions (e.g., rainfalls and currents can dominate in monsoon seasons and inflow currents can be the main factor during heavy rainfalls in mountains or foothills areas) [Heungsoo et al. 2013].
During summers, the warm water of epilimnion (upper well-mixed zone) is unable to mix with more dense water in hypolimnion (lower, colder zone) because of density differences that cannot be disturbed by the energy of winds, inflow or outflow currents and rainfall [Boehrer, Schultze 2008; Elci 2008; Hodges 2000; Lee et al.2013; Yu et al. 2010]. When the density differences arise in water column, the upper layer is usually influenced by winds and therefore well mixed. In contrast, the lower zone is cut-off from atmospheric oxygen and light. Photosynthesis may be inhibited, and the oxygen consumed by bacteria and other organisms. Near the bottom, anoxic conditions may occur. In such conditions, a solubility of nutrients, heavy metals and hydrogen sulphides increases and poses a risk of release from bottom sediments [Boehrer, Schultze 2008; Hudson, Kirschner 1997; Lee et al. 2013]. In addition to release from sediments, the organic matter falling from the epilimnion increases nutrients accumulation [Yu et al. 2010]. During the fall turnover, water mixing status changes as the solar radiation energy is outweighed by winds, currents and other factors causing vertical water flows. During the more intensive water mixing, substances accumulated in deep (hypolimnietic) water become available near the surface [Boehrer, Schultze 2008]. In addition to the seasonal changes in reservoir thermodynamics, shallow ones can stratify and destratify on a daily basis [Branco, Torgersen 2009].
Analyses of hydrodynamics, thermodynamics and impacts on water quality and ecosystems can be supported by mathematical modelling. There is a variety of numerical models describing the physical processes affecting water mixing and heat transport or coupled hydrodynamic and water quality models. Jensen et al. (2015) reported over 1500 such tools based on a literature review. There are numerous one-dimensional models (e.g., CE-QUAL-R1, DYRESM, DUFLOW, GLM, GOTM, LIMNMOD, MINLAKE, Mylake, PROTECH, SIMSTRAT) widely applied and in most cases, allowing for analyses of vertical water mixing and/or quality [Jensen et al. 2015]. These models were often a basis for two dimensional models such as CE-QUAL-W2, which can be used to simulate changes along rivers and narrow reservoirs, however, with the assumption that horizontal variations in water temperature and quality are of minor importance (not valid for analyses of local algal blooms) [Romero et al. 2004; Lee et al. 2013]. The most advanced tools are three dimensional models (e.g., AEM3D, CAEDYM, ELCOM, GEMSS, GETM) or modelling systems that allow to choose the number of dimensions according to the needs of application (e.g., Delft3D, EFDC, WASP) [Jansen et al. 2015]. Analyses of advantages and disadvantages of various models and modelling approaches were a subject of numerous studies (e.g., Saloranta et al. 2004; Gao, Li 2014; Jensen et al. 2015).
In the presented study, the ELCOM model was applied and calibrated (in terms of the water temperature) in order to simulate the summer and fall changes in water temperature of the Goczałkowice Reservoir and to prepare a basis for the comparisons of the reservoir’s thermodynamics and its chemical and ecological status. Mathematical model used and described in this study was applied as a part of task 5: ‘Development, application and verification of models for lakes physics, chemistry and ecosystems’ of the ZiZOZAP project.
Goczałkowice Reservoir is the biggest dam reservoir in the south of Poland. It covers over 32 km2 and has two main inflows: the Small Vistula River (the main inflow, approx. 80%) and the Bajerka River with total basin area of 530 km2 (Figure 1).
The reservoir was created in 1955 and serves as a major part of the system supplying the Upper Silesian agglomeration (approx. 3.4 million inhabitants) with potable water. It is also a storage reservoir protecting downstream areas from floods and droughts. Additionally, the reservoir being a part of the Natura 2000 system helps to protect a wide range of habitats and species [Dabioch et al. 2013; Młynarczyk et al. 2013; Polak et al. 2011], The reservoir is included in the national monitoring system and since 2010, extensive research monitoring has been carried out in the framework of the ZiZOZap project [ZiZOZap 2010], Both the operational monitoring system and real-time research measurements indicate strong changes in the water quality posing a risk to the water treatment plants and biodiversity dynamics (parameters of concern include: nutrients, sediments, dissolved oxygen and heavy metals). Some of the water quality parameters (e.g., phosphorus, humic acids and heavy metals) were recently analysed in relation to their abilities to be released from sediments [Dabioch et al. 2013; Młynarczyk et al. 2013; Polak et al. 2011]. However, the reservoir hydrodynamics and thermodynamics were not analysed in order to identify their impact on concentrations of chlorophyll
In this study, the Estuary, Lake and Coastal Ocean Model (ELCOM) was applied. ELCOM is a three-dimensional hydrodynamics model for lakes and reservoirs and is used to predict the variation of water temperature and salinity in space and time. Heat exchange through the water’s surface is governed by standard bulk transfer models. The energy transfer across the free surface is separated into non-penetrative components of long-wave radiation, sensible heat transfer, and evaporative heat loss, complemented by penetrative shortwave radiation. Non-penetrative effects are introduced as sources of temperature in the surfacemixed layer, whereas penetrative effects are introduced as source terms in one or more grid layers on the basis of an exponential decay and an extinction coefficient [Hodges, Dallimore 2013].
The model of Goczałkowice Reservoir presented here was applied (tested and calibrated) for a period of 6 months: June – November 2010. The time step used in simulations was set to 5 minutes. The starting time was determined by the setup of the real-time monitoring system, which was fully operational in June. The six-month period of analyses was chosen in order to be: 1) long enough to include processes crucial from the ecological and management points of view (e.g., algal bloom, high concentrations of suspended solids) and 2) as short as possible to enable multiple iterations required for the model calibration in an acceptable time.
The model of Goczałkowice Reservoir consists of 20 layers of thickness varying from 0.5 to 1.25 m. The thickness of 0.5 was applied to the surface water layers (up from 255.25 m a.s.l.), the thickness of 1.25 m was applied to the bottom layer (244–245.25 m a.s.l.) and the thickness of 1 m was used for ten intermediate layers. Twelve of these 20 model layers are usually wet, with the average water table elevation at 255.5 m a.s.l. The next three layers of 0.5 m thickness are reserved for higher water levels, with the maximum impoundment level at 257 m a.s.l. The top five layers should be always dry and were left only for calibration purposes. Horizontal resolutions of 25, 50, 100 and 250 m were tested and finally a resolution of 100 m was chosen as a compromise on the level of details and the calculation time (one ELCOM simulation of 180 days lasted 35 minutes for 250 m resolution, 5 hours for 100 m resolution and 9 days for the resolution of 25 m). The chosen resolution resulted in 119 columns, 60 rows and 41,683 total calculation cells. The model includes 7 inflows: the Vistula River (the main inflow), the Bajerka River and five pumping stations transferring the water excess from depressed agricultural and grassland areas around the western part of the reservoir (Figure 1). The average inflows in the analysed period are as follows: Vistula 7.98 m3/s, Bajerka 0.40 m3/s and pumping stations 0.41 m3/s. Outflows include: intake (average 2.02 m3s), spillway (average 7.07 m3/s) and bottom outflow required for the protection of ecosystems downstream (0.6 m3s). The temporal resolution of all inflow and outflow data is one day and all inflows include the information on water temperature.
The model included inputs from one meteorological station located above the water surface near the Bajerka river inflow. Meteorological data are of hourly resolution and are presented in Figure 2. The real-time water monitoring system launched in 2010 includes hourly measurements of water temperature in one water column with one-meter intervals (Figure 2). These measurements were the main input for the model calibration and verification. Additionally, the water temperature was measured periodically (approx. once a month) in 7 research monitoring points and in two monitoring points of the national monitoring system (Figure 3). These data were used only for model verification purposes.
The surface heat exchange includes shortwave radiation, longwave radiation, sensible heat flux and evaporation. The latter three are constituents of the non-penetrative energy density deposited in the surface layer. Solar radiation is one of the most important factors affecting the lake hydrodynamics and ecosystem. The ELCOM model divides the shortwave radiation into four bands: Photosynthetically Active Radiation (PAR), Near Infrared (NIR) and two Ultra Violet bands: A (UVA) and B (UVB). Percentages for each band in the whole shortwave radiation were set in the model as follows: PAR = 45%, NIR = 51%, UVA = 3.5% and UVB = 0.5%. The depth of the shortwave penetration depends in the model on the extinction coefficient, which is set separately for each band (Table 1).
Parameters used in ELCOM model. Photosynthetically Active Radiation (PAR) Near Infrared (NIR) Ultra Violet A (UVA) Ultra Violet B (UVB)Parameter Value Mean albedo of the water for shortwave radiation 0.08 Mean albedo of the water for long wave radiation 0.03 Drag coefficient on bottom cells 0.005 Initial wind speed for domain 0 m/s Initial water temperature for domain 11.4°C Light extinction coefficients for shortwave radiation bands used to calculate the shortwave penetration: 0.2-0.4 m–1 1 m–1 1-1.8 m–1 2.5-2.8 m–1 Sediments reflectivity 0.9 Surface heat transfer coefficient 0.0013-0.0015
The penetration was calculated according to the Beer-Lambert law [Hodges, Dallimore 2013]. In case of an excess of shortwave energy at the bottom of water columns, 90% of the energy is allowed to propagate back. For the longwave energy, two options were tested for the Goczalkowice Reservoir model: first based on the albedo and longwave radiation energy density deposited into the surface layer of water, and second, based on the albedo, cloud cover fraction and the temperature of the air and surface water layer [Hodges, Dallimore 2013]. The sensible heat loss from the surface of the reservoir was calculated based on the wind speed (at 10 m reference height above the water surface), the difference between the surface water temperature and air temperature, and the sensible heat transfer coefficient (Table 1). Finally, the last component of non-penetrative energy, the evaporation, was calculated based on atmospheric pressure and density, vapour pressure, saturated vapour pressure, water surface temperature, wind speed and heat transfer coefficient. The variability in the sensible and latent (evaporation) heat transfer described above can be affected by air column stability and water roughness [Imberger, Patterson 1990]. These factors have the effect of altering heat exchange coefficients, especially if meteorological parameters are measured over the water surface, as it is in the case of Goczalkowice Reservoir. ELCOM includes an atmospheric stability correction procedure applied after Hicks [1975] [Hodges, Dallimore 2013].
The surface thermodynamics described above are applied in ELCOM to the upper grid cells changing their temperature. Heat flux is added to the uppermost grid cell while solar radiation is added to the water column using exponential decay over depth. The heat transfer changes the density stratification, which is also affected by the concentration of suspended solids. Once the new density field is calculated, the mixing process is modelled on a layer-by-layer basis through each water column [Hodges, Dallimore 2013].
The model of Goczalkowice Reservoir excludes: the Coriolis term, reservoir generated tides and rainfall input. Whilst the excluding of the Coriolis term and tides may be easily explained by their minor importance for the analysed problem, the same could not be said of the rainfall input. Reasons for ignoring the rainfall effect on the reservoir thermodynamics are: 1) the geographical characteristics of the reservoir catchment (in this foothill area, the inflow currents and winds are the most dominant factors affecting the water mixing process), and 2) the observed low intensity of rainfalls in the analysed period. More about the factors determining the reservoir hydrodynamics and thermodynamics can be found in the ‘Results and discussion’ section.
ELCOM (and its successor – AEM3D) is a widely used model with teens of scientific papers confirming its applicability worldwide. Numerous applications are reported in Europe [Papadimitrakis, Karalis 2009; Lang et al. 2010; Imberger et al. 2017], and several in Poland, where it was used mainly to analyse the shallow dam reservoirs [Kliś et al. 2014; Karpiński, Łozowski 2017; Woźnica et al. 2017; 2017b]. Calibration of the Goczałkowice Reservoir model included: 1) modification of heat transfer coefficients, light extinction coefficients and initial water temperature in the reservoir or inflows, 2) using different meteorological data sources and 3) using different methods of longwave energy input to the surface layer (as described in the section above). The calibration process included over 400 iterations and the simulated water temperatures were compared to the output from the real-time monitoring system (Figure 4). It is a typical approach to calibrate the model’s thermodynamics prior to the application of the model’s results to the analyses of water quality or ecosystems [Zhang et al. 2008]. Such calibration and subsequent validation are often (but not as a rule, e.g., Leon et al. 2006) concluded using a statistical assessment of the model performance [Abbasi et al. 2016]. In case of the Goczalkowice Reservoir, model validation was based on both real-time and periodical monitoring and was evaluated using the determination coefficient (
From the beginning of summer (middle of June) to late summer (end of August), according to the monitoring data, the average water temperature in the surface layer is approximately 3.3°C higher than at the bottom of the reservoir, whereas the maximum difference (
The simulation of water temperature was followed by an analysis of relations between water temperature and observed concentration of chlorophyll
The presented study confirms the applicability of the ELCOM model to simulate thermodynamics of a dammed reservoir using the Goczalkowice reservoir as an example. Simulated and measured water temperatures are matched with high accuracy (
Basing on the high correlation between water temperature (water mixing processes) and concentration of chlorophyll