1. Introduction
The importance of preserving hydrological variability in environmental flow regimes has been described by several authors, e.g., [
1,
2], and methods to assess and preserve flow regime variability in environmental flows are suggested [
3]. The described problems with preserving hydrological variability in environmental flow regimes are valid for all kinds of encroachments on river regimes. Norway produces more than 96% of its electrical energy from hydropower and has several large-scale hydropower developments that have wide ranging effects in the affected river systems. Processes related to the implementation of the European Union Water Framework Directive and the upcoming relicensing of existing hydropower systems drive the need for better information on the impact of regulation and the development of better methods for mitigating changes in flow regimes. Numerous methods of different types and on different scales are developed for assessing impacts and designing environmental flows, e.g., [
4,
5]. One method used is the computation of hydrological indices that describes variability in the hydrological regime [
2]. The various indices can be linked to environmental performance [
6], and effects of changes in flow on the environment can be assessed. Applying these indices before (natural conditions) and then after flow changes provides a measure of change and a foundation for finding mitigation, and a basis for further and more detailed analysis of the impacts. Probably the best known method is the Index of Hydrological Alteration (IHA) described by Richter et al. [
7], a total of 33 indexes describing components of the flow regime relevant for ecological assessment.
Table 1 shows the category of indexes suggested by Richter. This method has gained some use and has been developed into a software tool for assessment and evaluation of the before and after situation in regulated rivers. An adapted version of this method is used as a component of the environmental flow methodology recently developed in Norway [
8], and also in the methodology for hydro-morphological assessment with regard to the implementation of the water framework directive [
9].
For assessment in cold regions, it is important to note that the common description of hydrological indexes of alteration [
7,
10] does not consider the significant effect formation of river ice on the flow conditions during winter. The formation and break-up of river ice can have a significant effect on the water level and flow regime in rivers due to the storage and release of water when ice forms and breaks up. Particularly, during freeze-up, ice can lead to large variability in water levels along the river under stable discharge conditions. A consequence of this is that the in-stream habitat can differ between summer and winter for the same discharge. For cold climate regions, knowledge of ice should therefore be combined with the open water hydrological variability analysis to ensure that relevant ice effects on the flow regime are taken into account in environmental flow analysis. This is particularly the case for rivers regulated for hydropower, since production releases have a significant effect on the ice regime after regulation [
11]. Research and data on discharge and river ice is most common in larger rivers from Canada and USA, but recent research has provided insights on the formation of ice and the effects ice has on the flow and related hydraulic parameters in rivers of the size most commonly found in Norway [
12,
13,
14,
15,
16]. Stickler, Alfredsen, Linnansaari, and Fjeldstad [
14] show how formation of anchor ice dams in the river Sokna in Norway transform the hydraulic conditions, and thereby the available habitat, under relatively constant discharge conditions. Shallow, fast riffle areas were transformed into a succession of deeper pools that eventually formed a continuous ice cover over the reach. During the formation process, the water level in this reach was considerably higher than what it would have been with the same discharge under open water conditions, and as a consequence anchor ice damming at this site would lead to ice-induced low flows further downstream. Similar results were found by Turcotte and Morse [
15], who also compared three sections of a different size and documented the similarities and differences in ice formation between scales from river to creek. Turcotte, Morse, Dube, and Anctil [
16] quantified the ice formation in small rivers combining an energy balance model with river hydraulics providing a method to extend the assessment to sections where direct observations don’t exist. Lind, Alfredsen, Kuglerova, and Nilsson [
12] used a statistical modelling approach to link ice formation on 25 river-reaches in Sweden to parameters easily obtainable from meteorological and geographical information, thereby providing tools for assessing ice formation over larger spatial scales. As discussed above, ice can induce high water conditions and even floods during freeze-up. In addition, ice can also induce flood levels during break-up, particularly in association with ice jamming [
17]. For a further discussion on low flows and the corresponding channel storage effects, see [
18].
There are both direct and indirect effects on the environment from ice formation which should be incorporated in an environmental assessment. It is therefore necessary to link the relevant hydrological indices to the response of the environment. For salmonid fish, ice formation and break-up can influence habitat availability, behaviour, and apparent survival, while a stable ice cover is often considered beneficial. A reduced or removed surface ice cover due to hydropower regulation in rivers where a stable ice cover is the normal can have a negative impact on Atlantic salmon [
19]. Summaries of the general impact of ice on fish can be found in reviews by Brown et al. [
20,
21], and Heggenes et al. The authors of [
22] recently summarized the knowledge on the relations between fish, hydropower development, and ice conditions. Formation of ice may also have effects on the terrestrial environment and the river floodplains. Lind and Nilsson [
23] found a higher diversity of riparian plants in rivers with flooding related to anchor ice formation, and also a difference in aquatic plants where anchor ice formed. Lind et al. [
24] review the relationship between ice dynamics on riparian and aquatic vegetation and find evidence of both a direct physical impact and the physiological effects of ice on plant communities. Lesack and Marsh [
25] show the relationship between flood peaks and water renewal in adjacent lakes on the Mackenzie delta, and how ice jams during the spring flood influence this process. Ice-induced effects on the flow regime can have a negative impact on the river environment, but also beneficial effects, which should be addressed in an environmental assessment. There is, therefore, a need to develop tools to do this, since omitting ice effects would eliminate effects like ice-induced flooding during freeze up, flooding during break up, the duration of the ice-covered period, and ice-induced low flows during freeze up, which are known to appear during the winter period in cold climates. In a study of the Peace–Athabasca delta, Canada, Monk et al. [
26] added cold climate relevant parameters to a subset of the original IHA parameters to assess environmental change also related to the ice processes. This work was later developed into the Cold-regions Hydrological Indicators of Change (CHIC) [
27]. This is a suite of indexes that characterizes the impact of ice on the flow regime, and, combined with a selection of IHA indexes, it covers the important features of a cold climate flow regime. This development work used the Mackenzie river, Canada as a case study.
Table 2 shows the CHIC revised ice indexes.
To compute the ice-related indexes, Peters, Monk, and Baird [
27] used data from Environment Canada that carried a flag in the meta data, which showed if the measurement was influenced by ice or not. For Norway, such data is not available from the national discharge database, and data on ice duration is just available for a few rivers and has uncertainties that are difficult to assess [
28]. On the other hand, ice-corrected data based on manual correction procedures are available from the national database on a daily scale and provide an alternative to data marked as ice-influenced.
The objective of this paper is to describe the effect of river ice on hydrological flow indices and suggest additional indices to measure ice effects, particularly those that focus on smaller rivers with short duration ice events. This paper also describes how indices can be used to both describe ice effects on flow in unregulated gauges and how to compare ice-induced hydrological variability before and after flow regulation. Further, this paper discusses temporal and spatial issues related to ice indices and how indices can be used to improve our understanding of winter flow regimes. Lack of data is a common issue related to winter and ice, and this paper discusses strategies for evaluating ice using air temperature, together with the flow record, to assess periods where ice effects on the flow can occur. This paper is an extended version of initial work presented at the CRIPE 2017 conference [
29].
3. Results
Figure 2 shows a hydrograph with ice-reduced data (blue), with a time resolution of one day and raw data (red) with an hourly resolution taken from the regulated Syrstad gauge. The large variability of flow over the winter period is caused by the operation of the hydropower system. Key indices are marked on the figure, blue circles indicate high peaks, green circles indicate low peaks, and open circles indicate peaks in the hourly series that are considered to be instantaneous and subsequently lost in the daily data series because of averaging. These are not included in the following analysis. The duration is computed as the width of each peak within the detection threshold, and counting the high and low peaks provides the frequency of events.
The formation of the identified ice-induced peaks on
Figure 2 coincides with periods of cold weather and shows that these peaks are driven by frazil and anchor ice accumulation at the gauge site. It is also worth noting that the peaks seem to appear independent of the volume flow in the river, which, for this winter, clearly shows the effect of hydropower operation. The indices for the high events registered in
Figure 2 (marked by blue circles) are summarized in
Table 3. All event data is presented in m
3·s
−1; this includes the rise and fall of the hydrograph (in m
3·s
−1·h
−1), in order to be consistent with the indices hydrological alteration commonly used, and since discharge is usually more readily available for these sites than water level. Data can easily be converted to water level in meters above sea level by utilizing the stage–discharge curve for the gauge. For evaluation purposes, e.g., purposes related to rises and falls, a conversion to water level can be a useful for computing, e.g., loss of wetted area or loss of habitat at the site.
Over the winters from 2000/2001 to 2015/2016, the maximum ice-induced peak was 289.7 m3·s−1, while the maximum difference between ice corrected discharge and the ice-induced peak was 238.9 m3·s−1, with an average of 35.3 m3·s−1. The median duration of all events longer than 5 h was 11 h. The average duration was 84 h, a number strongly affected by two cold years in which nearly the entire record was corrected. In such cases, the raw flow data would significantly overestimate the inflow to the river, while on the other hand the corrected flow record would underestimate the actual environmental conditions. It is particularly for events of long duration that an evaluation of the ice effects is important.
Figure 3 shows a comparison between flow records from the Syrstad gauge in panel A and the Hugdal bru gauge in panel B. In the Syrstad case, several peaks are seen forming in cold periods, raising the discharge above the correct discharge in the river, which is controlled by releases from the power plant. For the Hugdal Bru case (panel B), two freeze-up events are detected in November and December, separated by what is, most likely, drainage of the anchor ice dam formed at the gauge site (indicated by increase in temperature and fast drainage). The two other peaks occur during milder periods and could be effects of temporary ice jams at the gauge site, showing another type of event in which ice processes alter the discharge in the river outside of the natural stage–discharge relationship. The ice-controlled data for Hugdal bru shows a natural recession into a low winter flow.
For many practical comparisons, hourly data is not available, particularly in cases in which older discharge records are needed for the evaluation of pre- and post-regulation effects on the flow. It is, therefore, necessary to use ice adjusted daily data for the assessment.
Figure 4 shows a comparison for the Syrstad gauge between daily and hourly data.
Based on the data in
Figure 4, we see that the ice signal is also present in the daily averaged data, albeit somewhat reduced, but is still marked enough for analysis on a daily time scale to give an indication of the effects of ice on the discharge conditions. In cases in which an early spring flood affects the maximum peak during the winter season, the period of analysis could be shortened to eliminate this from the analysis. Other ice-induced peaks that occur before the start of the spring flood will still be detected. This can be seen by summarizing the number of events larger than 5 h for the winters from 2000/2001 to 2015/2016 (
Table 4), noting that winters 2004, 2006, and 2011 all show several peaks.
To compare changes in effects before and after hydropower regulation, indexes were computed for the pre-regulation and post-regulation period on daily data and the results were compared. All measurements are taken from the Syrstad gauge, but it is worth noting that the measurement site pre-regulation was originally a kilometer downstream from the current gauge and was then moved, since the old measurement site was affected by the construction of the Bjørset dam and intake. The series is later spliced by NVE and adjusted to the Syrstad site.
Figure 5 shows the 1-day peak discharge before and after the regulation. There is a significant difference (
p < 0.05) between peak sizes before and after the regulation. From the figure, we can see that ice had an effect on the flow also before the regulation, and that this has increased after the regulation due to loss of ice cover and a more dynamic ice regime. The effect of ice on the natural flow regime regarding timing of occurrence here is similar to what is seen in panel B of
Figure 3.
In
Figure 6, a similar comparison is made between the size of daily rises and falls in ice-induced peaks identified in the hydrograph. The difference between drops and rises are also significantly different (
p < 0.05), with larger drops and rises in discharge after the regulation. Regarding the duration of the peaks, there is no significant difference before and after the regulation. We get more and larger peaks, but the duration of each event is not very different in the regulated and unregulated case. The plot in
Figure 7 shows the winter mean index commonly used in the application of IHA in Norway for 10 years before regulation and 16 years after regulation. The plot shows the difference between uncorrected data and ice-corrected data on an important parameter in cold climate impact assessment, and it also illustrates the effect hydropower regulation typically has on the flow regime through increased winter discharge.
4. Discussion
For many evaluations of flow regimes in rivers during winter, open water data is used even if it is well known that formation and break-up of ice can have a large impact on the flow regime. This paper presents a set of indices that can be used to describe the effect ice has on flow in winter in rivers, information that may be important in environmental analysis during the cold season. The indices are mainly developed for assessing effects of flow regulations in winter, which is known to influence ice formation and ice production [
35]. The work builds on the work by Richter, Baumgartner, Powell, and Braun [
7] and follows the principles of the Index of Hydrological Alteration. The method presented will add functions to the first cold regions indices, the Cold Region Hydrological Indicators of Change, developed by Peters, Monk, and Baird [
27]. The indices presented here adds measures handling duration and frequency of ice-induced peaks and the rise and fall rates of such peaks. These are measures that are particularly relevant for smaller streams with dynamic ice regimes. The evaluation of the differences between recorded raw data and ice-corrected data as shown also highlights the importance of using the right winter data for the analysis and that winter data, without a clear status regarding if they are ice-influenced or not should, should be applied with care.
The derivation of the indexes and the analysis carried out here utilizes the official quality-controlled and ice-corrected data available for many Norwegian gauges operated by Norwegian Water Resources and Energy Directorate (NVE). The ice correction is carried out by NVE to provide discharge series that shows the true inflow from the catchment during winter. The analysis approach is an alternative to the method used by Peters, Monk, and Baird [
27], who utilized meta data available in the Canadian discharge database to identify ice-influenced flow periods. Since such information is not available for Norway; the approach used here is to combine the corrected data with the raw data showing the effect of ice formation. The analysis is therefore dependent on the accuracy in the ice corrections, and uncertainties in the corrections will lead to uncertainties in the computed indices. For the study site at Syrstad, the ice corrections should be precise, since both releases from upstream hydropower stations and the production discharge from a downstream station are available. In case such data is missing, a crude approximation can be made from temperature to identify periods in which ice may have an impact on the flow regime. As seen from the results presented here, this will mostly be an approximation, and little quantitative data can be reliably derived from this approach. An alternative approach in the case corrections or meta data on ice is missing—that is, applying a hydrological model to derive the natural inflow to the catchment as a basis for evaluating the magnitude of ice impacts. In this case, the simulated inflow at the gauge could be used as a replacement for the ice-corrected flow data used in this study. Independent of the method used, data for winter length and, particularly, ice cover stability may be difficult to assess without specific ice data for the sites.
The results from the analysis capture the ice dynamics of the flow regime as described in recent literature [
13,
14,
15], and the evaluation of differences between regulated and natural streams also confirms previous studies on extended ice formation period in rivers regulated from hydropower due to releases of warmer reservoir water into the downstream rivers. The results from this analysis show that, for the rivers studied, data with an hourly time resolution captures more of the dynamics of the ice-flow regime regarding the number of peaks, peak size, and the rate of change than what is seen in the data with a daily resolution. When data are subsequently aggregated into daily time resolution, the signal from short duration ice events are reduced or even removed from the data. On the other hand, since ice-corrected data with an hourly time resolution is not available, one should be careful not to attribute natural instantaneous peaks in the hourly data to ice driven events. For some sites, finer time resolution than one day could be important information related to environmental effects, e.g., Stickler et al. [
36] and Alfredsen et al. [
37] showed sub-daily responses in Atlantic salmon during several ice formation and release events in the river Orkla. The loss of detail in the daily ice data is, in many ways, similar to the loss of accuracy experienced when analyzing hydro-peaking flow data, e.g., Bevelhimer et al. [
38] found that a sub-daily time step is necessary to capture the dynamics of a peaking flow regime. As discussed, it is a drawback that the assessment of hourly peaking relies on daily ice-corrected data. Further, it can also be a problem that hourly data mainly are available for the recent decades, and such data are typically not available for the time periods that are often needed to define the pre-regulation flow conditions. However, as shown in
Figure 4 and
Figure 7, there is still an ice signal in the daily data, and it is possible to draw information on both timing and the effects ice has on the flow regime.
In a natural river, the winter season often can be divided into a freeze-up period, a stable mid-winter period, and a breakup period, of which each has distinct processes. The typical freeze-up effects are shown in
Figure 3B, and the same figure also shows effects of releases of ice on the flow measurements. In a regulated river with releases of production water from reservoirs, the mid-winter, stable ice period is often reduced or completely removed and replaced with a more dynamic ice formation [
21], as is seen in
Figure 3A. The method described here will catch the temporal distribution of ice effects, and as seen in
Table 3, events are timed and identified through the winter period.
A source of uncertainty associated with station-based analysis of hydrologic variability, as presented here, is if the method properly catches the spatial variability of the hydrologic regime in the river. This is even more important for assessment of ice effects than for the open water hydrological analysis, since the spatial variability of ice is larger than the variability of the flow alone, and is also driven by local hydraulic and morphological conditions and the energy balance of the river section. We can assume that observed, ice-induced low flows in a river probably mean that there are extensive ice formations and, likely also, ice-induced high flow at another location in the same river reach. This can be remedied by having a network of discharge gauges along the river, but, at least for Norway, there are few rivers that have the necessary spatial detail in their measurement networks. An option to overcome this is to model ice formation along the river reach and use this to develop a spatially distributed hydrological index. Timalsina, Charmasson, and Alfredsen [
11] showed the possibility of simulating ice development along the Orkla river using the MIKE-ICE model, and the findings could have been used to improve the hydrological assessment. The drawbacks of this approach is the significant data needed and the potential problems with setting up the model in smaller and steeper rivers with an even more dynamic ice regime. Turcotte, Morse, Dube, and Anctil [
16] presented an energy balance model that can be used to describe ice processes and quantify ice formation in steep rivers, which could be used as a tool to establish the ice regime and thereby better understand the spatial variability of the hydrological indices. This approach is an interesting way to better describe the ice regime of a river, but it also requires detailed data on climate and morphology with a spatial resolution that may not be available for many sites. A third option would be to use river classification or river mapping system on meso-scale, which is applied in many environmental studies, e.g., [
39,
40]. This is a rapid method of classifying river sections according to simple morphological and hydraulic parameters, and it could be combined with knowledge on ice processes to map probable sites where ice will have an effect on the flow regime, and then this information could be used to assess the spatial distribution of the ice effects measured at the gauge site. Typically, the classification systems’ map riffles, runs, and pools, which will have different ice formation processes, and this could then be related to the type of river found at the gauge site.
Figure 8 shows a section of the river Orkla mapped into river types. The Syrstad gauge is indicated on the map in a river type (shallow glide), which can have large frazil ice accumulations in winter, which explains the ice effects on the flow, as seen, e.g., in
Figure 2. By comparing river types, the extent of possible sites with similar ice regimes can be found, and, thereby, the spatial distribution of similar ice effects, as observed at the gauge site, can be evaluated. This method can be improved by more field data to confirm the ice formation processes. The linkage between river types and ice formation can be further developed, e.g., by utilizing a statistical model as described by Lind, Alfredsen, Kuglerova, and Nilsson [
12] to link simple parameters on river type to ice types and ice processes, and if more detail is needed this could further be developed into a spatial model using a strategy similar to the method described by Lindenschmidt and Chun [
41]. This could lead to a method that provides some modelling capability at the same time as it preserves the desired simplicity of the hydrological index method. A similar strategy could be to explore the possibility of linking the hydrological ice indexes to the conceptual ice model developed by Turcotte and Morse [
42]. This method links ice cover formation with ice formation processes using simple parameters like channel type, channel size, and winter intensity that can easily be found for any river. This model could also be combined with data from the river classification shown in
Figure 8 to give an overview of both the spatial and temporal aspect of ice formation. An important factor to consider when these models are applied to rivers like Orkla at Syrstad is the effect of the hydropower regulation, and a system like the one in Orkla with several reservoir and power plants will have a large effect on the flow regime and the spatial distribution of ice. This must be considered when existing tools developed for natural conditions are used in regulated rivers.
The indexes presented here will be directly useful as a complement to the existing tools when measuring the before and after situation in rivers with a regulated discharge. In addition, in the recent measures for hydro-morphological effects made for the water framework directive in Norway [
9], ice effects are categorized on a scale from natural to severely modified based on percentage changes in ice cover and/or anchor ice formation. This is data that can be difficult to find directly from measurements, so the derivation of data needed for this assessment must in many cases be based on proxy information. The indices of ice-induced flow effects could be used to evaluate the categories in this classification. The frequency, duration, and temporal distribution of ice-induced peaks could be related to reduced stability of the ice regime, and they are direct effects of the formation of frazil and anchor ice in the river section.
The use of hydrological indexes presented here is mainly focused on river regulations and a situation in which we compare the flow situation before and after the regulation to determine changes that can be critical to the environment and as a basis for finding and evaluating mitigation measures. From this application, it follows that the method also can be used to identify ice effects on a hydrograph and as a measure of the impact of ice on discharge and water levels also in natural rivers, e.g., as a tool for assessing effects of vegetation as shown by Lind and Nilsson [
23]. The method is also applicable to other types of studies in which measures of changes in the hydrology are useful. Bin Ashraf et al. [
43] use the index of hydrologic alteration to evaluate climate change effects on flow regimes in rivers in Finland. If ice is a case in climate change studies, the combination of changes in temperature and the current indexes of ice-induced flow changes could be used to assess future ice effects. Linnansaari et al. [
44] showed that ice had an effect on how Atlantic salmon utilized an artificial habitat in a regulated river, and unpublished data from the same river shows significant ice-driven modification to the water level and thereby the function of the habitat in the same river. The method presented here could be used to evaluate ice effects in such reaches by combining local measurements with a base flow assessment from an upstream unregulated gauge.