Heavy ice jams during spring floods are common on rivers in the northern European Russia. Climate changes of the last decades have caused significant changes in the ice regime of rivers. Continuous freezing periods, followed by sludge ice jams and late freezing at high water levels, occur with increasing frequency. Winter break-ups have become more frequent. The most recent catastrophic flood in the area under study took place in the spring of 2013. The probability of floods and ice jamming is evaluated; the duration and depth of inundation at various water levels are calculated with an assessment of the accompanying economic losses in the Northern Russia.
The use of water resources is inevitably accompanied by an increase in the risk of exposure of the public safety, as well as the social and industrial facilities on riverbanks due to the effect of hydrological processes. The details of the analysis of hazardous hydrological phenomena and consequently the assessment of the natural risk depend on the spatial scale of the study (Frolova et al., 2004). This scale determines the structure of the source information, the methods of its processing, and the choice of appropriate indicators and the features of the phenomena to be displayed on the map. A method for solving this problem is described here as applied to some developed reaches of the Northern Dvina River from the confluence of the Sukhona and Yug rivers, at Velikiy Ustyug one of the oldest cities in the northern Russia to Kotlas City at the inflow of the Vychegda River (Fig. 1). It stretches from the south to the north almost along meridian, so powerful ice jams form during spring floods on the Sukhona, Yug, and Northern Dvina practically every year. The control of accompanying inundations has been an ongoing issue there for more than a century. The major ice jams are often 10–20 km in length, and their lifetimes range from several hours to 3–5 days. Nine catastrophic floods, which caused enormous economic damage, were recorded within this reach from 1877 to 2013. The assessment of the possible flood risk required thorough analysis of historical, geographical, hydrometeorological, social, and economic information; field studies; the development of methods for calculation and prediction of hazardous ice phenomena in view of possible climate change; simulation of manifestations of dangerous hydrological processes based on different scenarios; the evaluation of the probable damage caused by the adverse impact of water, including the analysis of the efficiency of measures against ice and sludge ice jams; and eventually, the development of recommendations regarding protective measures.
The Northern Dvina is a typical river of the northern European Russia (Vuglinsky, 2002).The reach of the Northern Dvina under study, about 90 km in length, is located from the confluence of the Sukhona and Yug rivers (Velikiy Ustyug) to the inflow of the Vychegda R. (Kotlas). Five water stage gauges record water levels within this reach, and two of them (upstream) record water discharges. The rivers are mostly snow-fed; they show high spring floods and low winter low-water periods. No statistically significant changes were identified in recent years in the rate of stream flow and the character of its seasonal distribution.
The study area.
The probability of a dangerous water rise during high spring floods due to snowmelt inflow and the backwater effects of ice jams in the river reach as of 1946–2013 is about 30–35 % (Fig. 2).
An attempt was made to identify changes in river ice regime on the basis of hydrometeorological data covering two distinct periods: 1961–1990 (baseline) and 1991–2013 (present). The results of this study allowed trends to be identified and conclusions on possible qualitative and quantitative changes in river ice hydrological characteristics to be made. However, the quantitative change of various river ice regime characteristics is not so obvious and depends on the ice regime phase. The results of statistical analysis of the dates of floating ice occurrence and freezing up show, for example, that the shift in mean freeze-up dates in 1991–2013 increased by seven days as compared to the period of 1961–1990 (Fig. 3). In many cases, the change in freeze-up duration occurs due to later freezing (Fig. 3). Upward trend in dates and duration of freeze up period is statistically significant. Trend's significance has been evaluated using Spearman's rank correlation coefficient.
Such changes are associated not only with increase of air temperature in autumn prior to freezing, but also with the specifics of rivers regime. When the temperature regime for the study region becomes milder, the low water flows at the beginning of autumn are observed more frequently, which leads to early freeze-up dates. Air temperature fluctuations during autumn may cause a return of positive temperatures in November and result in floods due to snow melting and rain. High water levels during this period stimulate later freeze-up, increase of autumn ice and sludge drifts duration and formation of ice jams. Later freeze-up shortens the period of exploitation of ice crossing – an important feature of the north of Russia.
Changes in the flow component
Changes in freeze-up date
The milder winter conditions due to climate change result in the decreasing sum of negative air temperatures, increasing volume of solid precipitation and increasing number of positive air temperatures events that cause thawing spells (for many weather stations, the analysis of meteorological time series points out to the existence of statistically meaningful trends in these characteristics). The change in the maximum ice thickness during 1991–2013 compared to 1961–1990 is found to exceed 10 cm, which is statistically significant. The air temperature in winter often becomes positive and changes in sum of negative temperatures and ice thickness are more evident. On one hand, thawing spells lead to winter break-ups that are accompanied by catastrophic ice jams due to increased discharge and associated volume of sludge along river channels. Ice jams developed during such periods may cause subsequent spring ice jams. On the other hand, longer thaws result in thin ice cover and decreasing snow supply. High water levels during ice jam periods and large volumes of sludge form favorable conditions for hazardous ice jams during break-ups. Another example of adverse river ice regime changes is the reduction of ice crossings' exploitation periods due to winter break-ups and reduced thickness and strength of the ice.
The changes of the river ice regime during spring are also ambiguous. Spring
processes begin with the date of steady temperature transition above
0
The years of 2006 and 2007 showed examples of possible future ice regime changes. Despite ice appearance and freeze-up at the dates close to long-term mean ones, lengthy thawing spells in November and December 2006 contributed to ice cover destruction and break-up. The developed ice jams caused catastrophic increase of water level and flooding of developed territories in the Northern Dvina basin and created a risk of destruction of a motor road bridge at Kotlas city. Despite an early break-up in spring and low ice thickness, ice jams developed in the upstream area around Velikiy Ustyug town. Negative temperatures led to stabilization of the ice jam which lasted for 12 days despite of active measures to destroy it.
Measures to reduce water level rises caused by ice jams on the Northern Dvina were initiated since the 1950s. Works for ice weakening (ice blackening and cutting, as well as preventive explosive works on the most ice-jam-prone zones) are performed before break-ups. Ice jams were directly destroyed by explosions (including bombing) or weakened by preventive destruction of ice fields by icebreakers. However, the analysis of long-term fluctuations of the maximum levels induced by ice jams within the study reach over the recent decade does not show a significant positive effect of expensive ice-jam-control measures on the river rifts located downstream. Ice jams form on them almost every year and only in case of favorable climatic conditions (a relatively thin and weak ice cover, a low flood wave and slow snowmelt, and the passage of ice drifts on the Sukhona and Yug with a time lag), the spring flood on the Northern Dvina passes without significant damages.
Assessment of the flooded area is fulfilled by means of different
hydrodynamic models (Pietroniro et al., 2001; Beltaos, 2003). The available
observational data of 1896–2013 were used to calculate the maximum levels of
the flow and ice jam origin and maximum water discharges of 1 and 5 %
recurrence at the boundaries of the reach under study. Two-dimensional
Saint-Venant equations were taken as a basis for a mathematical model of
flooding. It has been developed by V. V. Belikov with colleagues since the
1990s. It is based on approaches of irregular hybrid computational meshes and
original methods of interpolation (Belikov and Semenov, 2000). The system of
Saint-Venant equations is most suitable for this class of problems, as it
takes into account a real bed bathymetry and floodplain topography, thus
allowing one to obtain the distribution of discharges over the branches of a
braided riverbed and to evaluate flow velocities and depths at any point of
the area. The simulation was performed using “STREAM-2D” software complex
(Belikov et al., 2002). The two-dimensional equations were discretized on a
triangular–quadrilateral mesh adapted to the channel network. The input data
included the results of field studies with depth measurements; the slopes of
water surface; measured water discharges and flow velocities, morphometric
cross-sections through river valleys, the locations and elevations of road
embankments, the presence of discharge openings; topographic maps with
The values of depths
In the calculations, the economic damage from the flooding of territories was divided into direct and indirect. The direct damage is associated with the direct physical contact of water with economic facilities, the amount of such damage being determined by the cost of restoration of the economy or the current market value of the destroyed (or damaged) economic facilities. The indirect economic losses are due to shutdown of facilities that process the production of plants, damaged by floods or inundations. The indirect damage also includes payments for shutdown periods, losses on transport, lost profits, adverse effects of flushed pollutants and damage to sewage treatment plants, manure storage structures, etc. The initial data for the assessment of damage to productive assets and finished products of companies, housing stock, and private property are data on the population in the flooded zone and elevations of the populated localities, as well as the type of residential and industrial buildings, which were classified according to the type of settlement. The lands subjected to flooding and waterlogging are defined in the riverside with the width corresponding to the rise of maximum levels of 1 and 5 % recurrence. The method of assessing the economic damage caused by floods, which is used in the Russian practice, was adjusted according to the available statistical information and economic and geographical characteristics of the territory (Shpagina et al., 2005). The calculations allowed estimating the total potential damage in the area of possible flooding. For example, in the case of a flood of 5 % recurrence, it turned out that 70 % of all costs are the costs of housing reconstruction and compensation to the population; 20 % are the costs of the restoration of production facilities; 10 % are the costs of rebuilding the social infrastructure. However, a comparison of the costs of eliminating the consequences of floods and the cost of construction of dams shows that the construction of flood control dams for the cities located on the banks of the Northern Dvina are effective only at a very high flood (1 % recurrence). The efficiency ratio (the ratio of the annual loss to the cost of protective structures) in this case is about 4–5.
The above methods of prevention of adverse impacts of water are basic and based on the major hydrotechnical intervention in natural processes of water and streamflow regime. They also involve significant environmental risks. The measures aimed to reduce the ice jam component during the passage of the wave of a spring flood, are not as effective but more economical and “soft” for river ecosystems. In particular, the prevention of flooding of settlements is optimal when hydraulic engineering is combined with annual dredging.
The analysis of current trends in the streamflow and hydrological regime of the Northern Dvina and the use of GIS technology and mathematical modeling of the flooding of floodplain–channel systems enabled formulating recommendations for preventing the inundation of populated localities and carrying out their critical evaluation. The results of model calculations provided a basis for recommendations and proposals for hydrotechnical engineering protection of populated localities by flood control dams and the assessment of their possible effectiveness and the ratio of the potential prevented damage to investments in the construction.
The study was supported by the Russian Science Foundation (Project no. 14-17-00700).