Challenges in Modeling Ice Floods Printer-friendly Version Interactive Discussion Hydrology and Earth System Sciences Challenges in Modeling Ice Floods on the Ningxia-inner Mongolia Reach of the Yellow River, China Hessd Challenges in Modeling Ice Floods Printer-friendly Version Interactive Discussi

Discussions This discussion paper is/has been under review for the journal Hydrology and Earth System Sciences (HESS). Please refer to the corresponding final paper in HESS if available. Abstract During winter the Yellow River in China is frequently subjected to ice flood disasters. Possible dike-breaking due to ice floods poses a serious threat to the part of the region located along the river, in particular the Ning-Meng reach (including Ningxia Hui and the Inner Mongolia Autonomous Region). Due to its special geographical location and 5 river flow direction, the ice dams and jams lead to dike-breaking and overtopping on the embankment, which has resulted in huge casualties and property losses throughout history. Therefore, there is a growing need to develop capability in forecasting and analysing river ice floods. Research into ice floods along the river is taking place at the Yellow River Conservancy Commission (YRCC). A numerical model is one of the 10 essential parts of the current research going on at the YRCC, which can be used to supplement the inadequacies in the field and lab studies which are being carried out to help understand the physical processes of river ice on the Yellow River. Based on the available data about the Ning-Meng reach of the Yellow River, the YRCC River Ice Dynamic Model (YRIDM) has been tested for capabilities to conduct ice flood forecasting. 15 The YRIDM can be applied to simulate water level, discharge, water temperature, and ice cover thickness under unsteady-state conditions. Different scenarios were designed to explore the model uncertainty for two bounds (5 % and 95 %) and probability distribution. The YRIDM is an unsteady-state flow model that can show the basic regular pattern of ice floods; hence it can be used as an important tool to support decision-20 making. The recommendation is that data and research should be continued in order to support the model and to measure improvements.


Introduction
River ice is a natural phenomenon that can be commonly seen in the cold regions of the world.River ice plays an important role in cold regions, being the main transportation means in wintertime: the ice roads and ice bridges caused by river ice in the northern Figures regions of Canada, Russia, and the USA (Alaska), where the population is sparse, are a positive effect of the ice formation in cold regions (Petrow et al., 2007;Rojas et al., 2011).However, on the negative side, river ice can cause ice flooding, hamper hydropower generation, threaten hydraulic structures, hinder water supply and river navigation, and has many other important detrimental effects.
When it comes to flood caused by river ice, the Ning-Meng reach of the Yellow River, China, cannot be ignored.The ice flood of the Yellow River basin is one of the most dangerous situations of all the Chinese rivers.According to historical data, different ice flood events have been recorded in the past during the winter on the main stream and the tributaries of the Yellow River.These floods have caused huge casualties and property losses, especially on the Ning-Meng reach of the Yellow River.Since 1949 there have been 29 ice flood hazards, 9 of them causing the levee to breach in the Ning-Meng reach.
Ice regime forecasting is an efficient method for ice flood control, regulation and decision making, and mathematical modeling is playing an essential part in the development of river ice research.Mathematical models can supplement the inadequacies of field and lab studies to help understand the physical processes of river ice, and at the same time, they can also be a tool to help design and plan engineering projects.They have become valuable tools for exploring the research area of river ice, such as understanding the physical processes and simulating river ice phenomena, and even forecasting ice floods (Dahkle et al., 2012).
The development of the economy and society with the resulting climate change and human activities have changed the characteristics of ice regimes, especially ice disasters during the freezing or breaking-up periods, attracting more and more attention from water authorities and local governments.Hence, it is very important to know ice regime characteristics and use mathematical models to enable ice forecasting, ice flood prevention and ice flood alleviation.
At the YRCC a numerical model named YRIDM has been developed to simulate ice regimes.It also contains data on observed data including hydrological, meteorological, Introduction

Conclusions References
Tables Figures

Back Close
Full and ice regime data, covering about fifty years of measurements.The available model and collected field data provide an opportunity to examine the possibility of ice flood forecasting.
The purpose of this paper is to apply the YRIDM to the Ning-Meng reach of the Yellow River to examine the accuracy of the model.A secondary objective is to examine the possibility of applying the YRIDM to forecasting ice flooding in support of decision making.
This paper first describes the current state of the art in ice flood modeling, followed by a description of the case study.The modelling application is detailed next.Modeling of ice floods on the Ning-Meng reach is discussed with respect to model set-up, calibration, verification and uncertainty analysis.This paper ends with a conclusion section.

Ice flood research
Ice floods in rivers are generated due to ice jams, ice dams, and ice-snow melt (Liu et al., 2000).The ice flood process cannot be separated from the whole river ice process, although it mainly involves the ice-covered period and breaking-up period.
As summarized by Shen (2006), during the past fifty years engineering and environmental issues have largely driven the development of river ice research, and significant achievement has been made during this period.However, Beltaos (2008) thought that although there has been remarkable progress in understanding and quantifying the complex river ice processes, many problems concerning river ice still remain today or are only partially solved.Shen (2006) classifies the river ice research area into two parts: firstly, the energy budget methods and water temperature distribution calculations before and during the freezing-up period and secondly, the evolution of frazil ice, frazil floc, anchor ice and Figures

Back Close
Full ice dams.Currently the water temperature distribution is well understood (e.g.Shen and Chiang, 1984), and the mechanisms of super cooling and frazil ice formation are also relatively well understood (Osterkemp, 1978;Daly, 1984).However, the evolution of frazil ice, frazil floc, anchor ice and ice dams needs to be studied further (Ye et al., 2004), the transitional conditions among different ice run regimes are not clearly understood (Hammar et al., 2002), and also knowledge on the mechanism of ice pans and ice floe formation is limited.In addition, a complete analytical formulation of the mechanical breakup should be developed.Beltaos (2008) showed the challenges and opportunities for research into river ice processes.The main challenge is to avoid or decrease the negative influence of river ice processes and to make sure that the positive influence is not affected by human activities (Beltaos, 2008).In order to meet this challenge, it is necessary that there is a good qualitative understanding of river ice processes.However, there still remain serious gaps, such as research into anchor ice, breakup, ice jamming, and the influence of climate change on the river ice process.The main opportunities lie in how to use new technologies to understand the river ice processes.These new technologies include instrumentation, numerical modeling, mitigation and prediction of climate impact on river ice processes.
Research into ice cover has developed from static ice cover research, such as Beltaos and Wong (1986), into dynamic ice cover research due to the fact that the original research did not take ice dynamics into consideration.Dynamic ice cover research has been further divided into one-dimensional dynamic ice cover research (Shen et al., 1990) and two-dimensional dynamic ice cover research, such as the DynaRICE model (Shen et al., 2000), in order to take into account the frictional resistance of riverbanks and channel bottoms, irregular cross-sections of river channels, and unsteady flow state in reality.
Research into ice transportation under the ice covers and ice jams started from a critical velocity criterion or Froude Number criterion (Kivislid, 1959;Tesaker, 1975).However, this method could not provide a way to compute an accurate value of the Figures

Back Close
Full critical velocity and Froude Number.Hence, Shen and Wang (1995) developed a concept of ice transport capacity to simulate ice deposition under the ice cover, which was demonstrated by field observation data.Nevertheless Beltaos (2008) still thought that ice transportation under the ice cover was only partially understood, and further research should be continued in the future.
The process of the thermal growth and erosion of ice cover is relatively clearly understood.Firstly, research focused on the ice cover without layers of snow ice, snow slush or black ice (Shen and Chiang, 1984).Secondly, the layers of snow ice, snow slush and black ice were taken into consideration (Calkins, 1979;Shen and Lal, 1986).Finally, based on the previous research, a lot of models were produced to simulate the processes.The most classical one was the degree-day method proposed by Stefan (1889) that has been used to simulate ice cover growth for a long time.Shen and Yapa (1985) refined the classical degree-day method which they then developed into a modified degree-day method, which was demonstrated on the St. Lawrence River.
The breakup of ice cover can be classified as mechanical breakup or thermal breakup.A mechanical breakup can exert a negative influence on the hydraulic facilities and the safety of people along the river.However, the ability to simulate a mechanical breakup is still limited (Shen, 2006).
Only a few researchers have tried to simulate the propagation of ice jam release waves.Most of them have used a one-dimensional model without considering the effect of ice (e.g.Blackburn and Hicks, 2003).Field investigations were conducted by Jasek (2003), who found out that release wave celerity seemed to change with different ice conditions.Based on this research output, Liu and Shen (2004) started to take into consideration the effect of ice on wave propagation.She and Hicks (2006) built a model, named River-1-D, to look into the effect of ice on ice jam release waves, which was tested by the release event on the Saint John River in 1993, on the Saint John River in 2002, and on the Athabasca River in 2002.At the same time, statistical methods and ANN are also used to forecast the breakup of river ice.The potential for using Fuzzy Expert Systems to forecast the potential risk of ice jams was discussed Figures  (Mahabir et al., 2002), and the systems identified five years when high water levels could occur, including four years when the predicted high water level did occur.Based on this potential, more research has been carried out (Mahabir et al., 2006;Wang et al., 2008) which took more related parameters into consideration.

Ice flood modeling
Among all the river ice research, mathematical modeling is an essential part of the progress.Ice flooding models have developed from 1-D steady state to 1-D unsteady state and then into 2-D models, and data-driven models are also being applied to forecast ice floods.The 1-D steady state models are based on the static ice jam theory, namely that the formulas to determine the final thickness of water surface ice can be deduced according to the static balance of internal and external forces on a floating ice block.The models based on the static ice jam theory are HEC-2 (US Army, 1990), ICETHK (Tuthill et al., 1998), and HEC-RAS (US Army, 1998).The basic assumption is that the flow is steady, gradually varied, and one-dimensional, and that river channels have small slopes (less than 1 : 10).The basic equation is the force and energy balance equations; a standard step method is used to solve the equation.However, the limitation of 1-D steady state models is that they ignore ice dynamic conditions.Due to the limitation of 1-D steady state models, the development of 1-D unsteady state models, such as RICE (Lal and Shen, 1991), RICEN (Shen et al., 1995), and the Comprehensive River Ice Simulation System (CRISSP1D) model (Chen et al., 2006), was promoted.These models are based on the assumption of unsteady state flow; the governing equations are 1-D Saint Venant equations, i.e. mass and momentum conservation equations with floating ice, and they can be used to simulate the entire ice process in rivers during the winter season.The basic assumptions of the RICE, RICEN, and CRISSP1D models are that they ignore the following: (1) the effect on the water body mass balance due to the change in ice phase; (2) the river has a floating ice cover; (3) two-layer ice transportation theory; (4) suspension ice is full of the suspension ice 12299 Introduction

Conclusions References
Tables Figures

Back Close
Full layer; and (5) the thickness of the ice layer on the water surface is equal to the thickness of the ice block floating on the water surface.However, the limitation of these models is that they lack detailed consideration of complex flow patterns and river geometry (Shen, 2010).The RICE and RICEN models have been used to simulate the ice regime for a few rivers; for example, the RICE model was used for the Upper St. Lawrence River near New York and the Peace River in northern British Columbia (Li et al., 2002) and the RICEN model was used for ice jam simulation in 1995 on the Niagara River.
Based on the limitations of 1-D unsteady state flow, 2-D models have also been developed, such as DynaRICE (Shen et al., 2000) and CRISSP2D (Liu and Shen, 2006) that can be used to simulate the ice regime.These models solve the 2-D depthintegrated hydrodynamic equations for shallow water flow.The two basic assumptions are that the movement of the surface ice layer is continuous, and ice is a kind of continuous medium.A finite element method is used to solve the equations.However, the limitation is that they lack detailed consideration of the third dimension.
DynaRICE has been applied to understand the ice jam evolution, ice boom design and navigation structure design, on several rivers.DynaRICE was applied in the Niagara Power Project to study both the ice control and ice-period operation, on the Missouri-Mississippi River to study ice jam formation, and on the Shokotsu River in Hokkaido, Japan, to study breakup jams.CRISSP2D has also been applied to several rivers and lakes to simulate the ice regime.It has been applied to the Nelson River to simulate freeze-up ice conditions (Malenchak et al., 2008); on the Red River near Netley Cut (Haresign and Clark, 2011) to model ice formation; and to test the potential for anchor ice growth after the Conawapa Power Generation Station was constructed and operational (Morris et al., 2008).
The potential for using Fuzzy Expert Systems to forecast the potential risk of ice jam was discussed by Mahabir et al. (2002).Based on this research, more research has since been carried out (Mahabir et al., 2006;Wang et al., 2008).Data-driven modeling has been applied to the ice regime simulation several times and one example is the 12300 Introduction

Conclusions References
Tables Figures

Back Close
Full Yellow River (Chen et al., 2012).However, the limitation of using such methods is that ice processes cannot be understood by data-driven models.
3 Case study description

The Ning-Meng reach
The Yellow River is the second longest river in China and is also one of the most famous rivers in the world.Due to the fact that it is the cradle of the Chinese civilization, the river is also called the "Mother River of China".The Yellow River originates from the Bayanhar Mountain in Qinghai Province and flows through the nine provinces of Qinghai, Sichuan, Gansu, Ningxia, Inner Mongolia, Shaanxi, Shanxi, Henan, and Shandong, and finally flows into the Bohai Sea in Shandong Province.It has a total length of 5464 km and a basin area of 752 443km 2 .
The Ning-Meng reach of the Yellow River (Fig. 1) is where the ice flood mainly happens.The Ning-Meng reach length is 1237 km and it consists of two consecutive parts, the Ningxia and Inner Mongolia reaches.The starting point for it is Nanchangtan, Zhongwei County in the Ningxia Hui Autonomous Region, and the end point is Yushuwan, Mazha Town, Zhungeerqi County in the Inner Mongolia Autonomous Region.
The Qingtongxia and Sanshenggong reservoirs are located on this reach.The total length of the Ningxia reach is 397 km and it flows from southwest to northeast.On this reach one part, from Nanchangtan to Zaoyuan, does not usually freeze because the bed slope is steep and the water velocity is large, and it only freezes in a very cold winter.The rest of the reach, from Zaoyuan to Mahuanggou, flows from south to north and it often freezes because the slope is gentle and the water velocity is low.The Inner Mongolia reach is located on the north of the Yellow River basin between 106 twists and turns.Although it is located in the middle and lower reaches, the slope of the reach that is under study is close to that at the Yellow River Estuary.

Past ice floods
According to historical data for the Yellow River, ice disasters appeared frequently in the past and every year between 1855 and 1949, when dams were destroyed by ice floods 27 times.In the period between 1951 and 1955 the dams on the Lijing breach were also destroyed by ice flood.In addition, there were 28 ice flood seasons with ice disasters between 1951 and 2005 (Rao et al., 2012).
Especially on the Ning-Meng reach, ice jams, frazil jams, and other ice flood disasters have happened frequently.According to the statistics for ice flood disasters on the Ning-Meng reach there were 13 ice flood events from 1901 to 1949, almost one event every 4 yr.Even after the Liujiaxia Reservoir and the Qingtongxia Reservoir started operation in 1968 and 1960 respectively, ice flood disasters still occurred in 1993, 1996, 2003, and 2008(Gao et al., 2012)).
The Yellow River basin is under the influence of cold air from the vicinity of Siberia The governing equations that are used by YRIDM are: (3) (4) The reach between Bayangaole station and Toudaoguai station was chosen to be modelled due to the fact that ice flood has happened frequently here and data was available on the reach.The simulation reach has a total length of 475.6 km.A hydrological station named Sanhuhekou is positioned in the middle of the reach, i.e. 205.6 km away from Bayangaole station.
The available data includes river hydraulics, meteorological, ice regime, and cross sectional at the four hydrometric stations, namely Shizuishan, Bayangaole, Sanhuhekou, Toudaoguai stations, and it covers the period of ten winters, from 2001 to 2011.River hydraulic data includes water level and discharge with daily measuring frequency.Meteorological data includes air and water temperature with a daily measurement of the water temperature and twice a day measurement of air temperature (the daily highest and lowest temperature per day).Ice regime data includes ice run data, freeze-up date, break-up date, and ice cover thickness.The measured frequency is per winter except that the measured frequency of ice cover thickness is every 5 days.and Toudaoguai stations, and cross sections and bed elevation were used to build the model.

Model calibration and verification
During the winter of 2008/2009, daily mean discharge at the Bayangaole station was input into the model as the upstream boundary condition, and daily mean water level at the Toudaoguai station was set as the downstream boundary condition.The air temperature at the Sanhuhekou and Toudaoguai stations were input into the model, and the air temperature at the Sanhuhekou station was used to cover the reach from the Bayangaole station to the Sanhuhekou station, and at the same time, the air temperature at the Toudaoguai station was used to cover the reach from the Sanhuhekou station to the Toudaoguai station.The daily lowest and highest temperatures were interpolated by sine function.Daily mean water temperature at the Bayangaole station was input into the model as the upstream boundary condition of water temperature.The initial conditions were water temperature, water level, and discharge.Based on the literature review and the reference book of the numerical model from YRCC, the sensitivity parameters are listed in Table 1.The One-At-A-Time sensitivity measure method was used to conduct the sensitivity analysis, which means the value of one parameter is changed from the minimum value to the maximum value while at the same time other parameters are kept constant at the mean value, and then the variation of the model output is checked.
After the sensitivity analysis, scenarios of the model calibration can be designed based on the sensitivity analysis results (Table 2).If the model is sensitive to a certain parameter, then the parameter needs to be calibrated carefully, otherwise the default value is used.In this case, bed roughness, end ice roughness, and decay constant are sensitive to the water level at the Sanhuhekou station; heat exchange coefficient between ice and air, and heat exchange coefficient between water and air are sensitive to ice cover thickness at the Sanhuhekou station.Therefore these 5 parameters should be calibrated carefully; heat exchange coefficient between water and ice is not sensitive 12305 Introduction

Conclusions References
Tables Figures

Back Close
Full

Model uncertainty analysis
Based on the sensitivity analysis results, the water level at the Sanhuhekou station is sensitive to the Manning coefficient of river bed, decay constant, and end ice roughness; and the ice cover thickness at the Sanhuhekou station is sensitive to the heat exchange coefficient between ice and air, and heat exchange coefficient between water and air.Hence, the uncertainty analysis is divided into an uncertainty analysis of the water level at the Sanhuhekou station and the uncertainty analysis of ice cover thickness at the Sanhuhekou station, and the Monte Carlo simulation is used to conduct the parametric uncertainty analysis.
Due to the limitation of time, it was impossible to run the model several times.When conducting the uncertainty analysis of the water level at the Sanhuhekou station, the related four parameters were the Manning coefficients of the river bed at the upstream 12306 Figures

Back Close
Full and downstream of the Sanhuhekou station, decay constant, and end ice roughness.
The scenarios of the uncertainty analysis of the water level at the Sanhuhekou station were designed based on the calibrated parameters, the range was calculated by increasing and decreasing the calibrated value by 20 %, the sample generation was uniformly random, and the number of simulation was 500.When conducting the uncertainty analysis of the ice cover thickness at the Sanhuhekou station, the two related parameters were the heat exchange coefficient between ice and air, and the heat exchange coefficient between water and air.The scenarios of the uncertainty analysis about ice cover thickness at the Sanhuhekou station were designed based on the calibrated parameters, the range was calculated through increasing and decreasing the calibrated value by 20 %, the sample generation was uniformly random, and the number of simulations was 400.Based on the above case designs, the parameters were input into the model, the model was run and the result were stored, the distribution and quartile of the output were analysed, namely the PDF (Probability Density Function) at one time step and two bounds (5 % and 95 %).

Results and discussion
The results of the sensitivity analysis show that bed roughness, end ice roughness, and the decay constant are sensitive to the water level at the Sanhuhekou station, and that the heat exchange coefficient between ice and air, and the heat exchange coefficient between water and air, are sensitive to ice cover thickness at the Sanhuhekou station.
Based on the results of the sensitivity analysis, the model was calibrated, and the values of calibrated parameters can be seen in Table 3.
Figure 3  with the observed results, which is because in the beginning, the simulated water level is affected by the initial conditions.From day 125 to the end, the simulated result is not accurate and the trend is even adverse, which is because the model cannot simulate an ice jam breakup during the breakup period.The data observed during the period clearly shows that the water level decreases suddenly, which means that the ice jam at the downstream of the Sanhuhekou station collapses, but the model cannot simulate this phenomenon, which is why the simulated results are not so accurate.During the rest of the period, the results are acceptable.
Figure 4 shows a comparison between the simulated discharge data and observed discharge data at the Sanhuhekou station during the winter of 2008/2009.The RMSE is 160.5 compared with the peak discharge of 1400 m 3 s −1 , so the error rate is 11.2 %, and hence the simulated result is acceptable and reasonable.From day 1 to day 25, the simulated results are not accurate, but the trend in the simulated results is consistent with the observed trend, which is because in the beginning, the simulated water level is affected by the initial conditions.From day 25 to day 50, the simulated results are also not accurate.The discharge decreases suddenly; this is because during the freeze-up period, the channel storage capacity increases when the water has been changed into ice, then the amount of water decreases, which can result in a sudden decrease in discharge, but the discharge can increase again when the channel storage capacity of the river continues to be stable.However, there is an assumption that the model ignores the effect on the water body mass balance due to changes in the ice phase; this is the reason why the results are not good enough.During the other period, the result is acceptable.
Figure 5 shows the comparison between the simulated data and observed data for the water temperature at the Sanhuhekou station during the winter of 2008/2009; the RMSE is 0.854.The simulated results are acceptable and reasonable.From day 1 to day 15, the simulated results are not good, which is because in the beginning, the simulated water temperature is affected by the initial water temperature.From day 55 to day 95, the water temperature is below 0 • C; this is the super cooling phenomenon, 12308 Introduction

Conclusions References
Tables Figures

Back Close
Full which means before the formation of ice cover the water body loses heat very quickly, which can result in the negative value for the water temperature.However, once the ice cover forms, it can prevent heat exchange between the water body and the air; this is the reason why the water temperature keeps steady at 0 • C between day 100 and day 130.During the rest of the period, the results are acceptable.
Figure 6 shows the comparison between the simulated data and observed data for ice cover thickness at the Sanhuhekou station during the winter of 2008/2009.Because the ice cover thickness data was measured per five days, hence the observed ice cover thickness data is not continuous and the RMSE between the simulated data and observed data is 0.109 m.The simulated results can describe the variation trend for ice cover thickness, and the simulated maximum value of the ice cover thickness is 0.5786 m on 28 January 2009, which is 0.0036 m higher than the observed measurement (0.575 m) and 4 days earlier than the observed measurement (1 February 2009), hence the result is acceptable.
The prepared data during the winter of 2007/2008 were input into the model to verify the calibrated model.However, the YRIDM cannot work.According to the reference book for the model, the model cannot work when the water level exceeds the height of the embankment.In the winter of 2007/2008, the water level exceeded the height of embankment at Duguitalakuisu county, and resulted in a dike breach.Hence, the YRIDM should be improved so that it has the ability to deal with the problem.
After running the model, the uncertainty analysis results can be summarized in Fig. 7. Figure 7 shows the observed data, 5 % percentile bound, 95 % percentile bound, and the results of 500 cases.Day 90 (Fig. 8) and Day 130 (Fig. 9) are chosen to show the probability distribution.When it comes to the probability distribution of the uncertainty analysis results, if the distribution looks like a normal distribution, the uncertainty analysis results are good.According to Figs. 8 and 9, the probability distribution of the water level on day 90 looks good, which is because it looks like the normal distribution.Although the probability distribution of the water level on day 130 is a skewed normal 12309 Introduction

Conclusions References
Tables Figures

Back Close
Full distribution, at least it shows the basic shape of normal distribution, and hence the uncertainty analysis results are reasonable.After running the model, the uncertainty analysis results can be summarized in Fig. 10. Figure 10 shows the observed data, 5 % percentile bound, 95 % percentile bound, and the results of 400 cases.Day 60 (Fig. 11) and Day 120 (Fig. 12) are chosen to show the probability distribution.According to Figs. 11 and 12, the uncertainty analysis results are not good, due to the fact that the probability distribution of the ice cover thickness on these two days is not a normal distribution.This is because the number of cases designed for the uncertainty analysis is 400, which is not sufficient to show the characteristics of the distribution.

Conclusions
In this paper, a one-dimensional model named YRIDM was applied to the Ning-Meng reach of the Yellow River to simulate the river ice processes.After the model had been set up, a sensitivity analysis was conducted by the One-At-A-Time sensitivity measure method, and the results showed that bed roughness, end ice roughness, and decay constant are sensitive to water level at the Sanhuhekou station; and that both the heat exchange coefficient between ice and air and the heat exchange coefficient between water and air are sensitive to ice cover thickness at the Sanhuhekou station.Based on the sensitivity analysis results, the model was calibrated, and the results showed the acceptable capacity of the YRIDM to simulate the water level, discharge, water temperature, and ice cover thickness.However, the model verification cannot be conducted due to the model weakness, namely the model cannot work when the water level exceeds the height of embankments.Nevertheless reasonable results about uncertainty analysis can be achieved.Hence the model can be applied to the Ning-Meng reach to simulate its ice regime.Once the model has been calibrated, it can be used to forecast the ice regime to support decision making, such as on artificial ice-breaking and reservoir regulation.The limitation of the YRIDM model is that the effect on the water body mass balance due to ice phase change is ignored because of the basic premise by the model which ignores the effect on the water body mass balance due to ice phase's change, hence the model cannot simulate a mechanical breakup during the breakup period.
The recommendation is that the measured items and frequency should be increased, and then more research can be conducted.Based on the data and research, the model can be further improved, for example, the one-dimensional models extended to two-dimensional models, and an ice flood early warning system built to support the decision-making.Introduction

Conclusions References
Tables Figures

Back Close
Full  Full  Full       12329 Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | • 10 E-112 • 50 E, 37 • 35 N-41 • 50 N.The total length of the main stream is 830 km, and the total height difference is 162.5 m.The reach is wide with a gentle slope and meandering 12301 Discussion Paper | Discussion Paper | Discussion Paper | and Mongolia; the mean air temperature during the whole winter is below 0 • , and this lasts for 4 to 5 months.Due to the special geographical location and difference in latitude of the Ning-Meng reach on the Yellow River, together with the special river flow direction, the downstream part both freezes and breaks up earlier than the upstream part.Hence it is very easy for ice jams or backwater to be formed, which can result in ice flood disasters such as dams being destroyed and dike breaks.Based on the literature review about ice flood disasters, one can conclude that on the Yellow River the main problems regarding ice floods are caused by ice dam floods and ice jam floods which can result in dike breaks and overtopping on the embankment.This is the key problem to be solved on the Ning-Meng reach of the Yellow River.Present research uses the YRIDM, which is a 1-D unsteady state model.The YRIDM is based on the RICE and RICEN models, and designed taking river ice processes into account.It was designed to simulate the ice flood with unsteady-state flow, and its logical flowchart is presented in Fig.2.It consists of three main components, the river hydrodynamics, thermodynamics, and ice dynamic modules.The advantage of this model is that it can be subdivided even further into the following modules: river hydraulics, heat exchange, water temperature and ice concentration distributions, ice cover formation, ice transport and cover progression, undercover deposition and erosion, and thermal growth and decay of ice covers.
e )e −α n T (5) Where, Q = discharge (m 3 s −1 ); A = net flow cross-sectional area (m 2 ); x = distance (m); t = time (s); g = gravitational acceleration (m s −2 ); S f = friction slope; n c = equivalent roughness; n i = roughness of ice cover; n b = roughness of river bed; u = average velocity (m s −1 ); R = river hydraulic radius (m); n i, i = initial roughness Discussion Paper | Discussion Paper | Discussion Paper | of ice cover; n i, e = final roughness of ice cover at the end of an ice-covered period (0.008 ∼ 0.012); T = freeze-up time (day); α n = decay coefficient.Based on the available data analysis for the Ning-Meng reach of the Yellow River, and literature review about ice models, the YRIDM has been selected to conduct the modeling for ice flood middle-term prediction (10-15 days ahead).The ice regime on the Ning-Meng reach is mainly determined by thermal factors, dynamic factors, channel course conditions and human activities.During the whole research period between 2001 and 2011, the most serious ice flood occurred during 2 winters, namely the winter from 2007 to 2008 and the winter from 2008 and 2009.The winter from 2008 to 2009 was chosen as the simulation winter to calibrate the model, and the winter from 2007 to 2008 was used to verify the calibrated model.
Data for four cross sections for 2009 is available.During the winter of 2008/2009 the discharge and water temperature at Bayangaole station, the water level at Toudaoguai station, air temperature at Sanhuhekou 12304 Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | shows the water level comparison at the Sanhuhekou station during the winter of 2008/2009 between the simulated data and observed data.The resulting RMSE is 0.464.The simulated result is acceptable and reasonable.From day 1 to day 40 the simulated result is not good, but the trend in the simulated results is consistent Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Figure 1 .
Figures and tables refered in text
to the output, hence the default value can be set.The model can be calibrated under different ice regime conditions, which means that some parameters can be calibrated under open water conditions, such as the Manning coefficient of river beds and the heat exchange coefficient between water and air.Some parameters can be calibrated under ice conditions, such as the Manning coefficient under ice cover, end ice roughness, and decay constant.Hence the model calibration procedure is divided into a model calibration under open water conditions and a model calibration under ice conditions.In the calibration procedure, RMSE (Root Mean Square Error) is used as the criterion to check the model performance.After calibration of YRIDM during the winter of 2008/2009, the discharge and water temperature at the Bayangaole station, the water level at the Toudaoguai station, air temperature at the Sanhuhekou and Toudaoguai stations, and cross sections and bed elevation during the winter of 2007/2008 were used to build the model.The simulated water level, discharge, water temperature, and ice cover thickness at the Sanhuhekou station during the winter of 2007/2008 were used to compare with the observed values to verify the calibrated model.

Table 1 .
Sensitivity parameters of the YRCC River Ice Dynamic model.