Icedammed lake outburst floods in the Altai mountains, siberia – a review with links for further readings
Reconstruction of the size of the flood formed during the break of the icepond lake. Study of the data on the flow rates and velocities of the water corresponding to the time of passing the peak of the flood and the various stages of its downward phase.
Рубрика  Иностранные языки и языкознание 
Вид  статья 
Язык  английский 
Дата добавления  24.11.2018 
Размер файла  587,3 K 
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
150 Jьrgen Herget
Размещено на http://www.allbest.ru/
УДК 574
University of Bonn (Bonn, Germany)
ICEDAMMED LAKE OUTBURST FLOODS IN THE ALTAI MOUNTAINS, SIBERIA  A REVIEW WITH LINKS FOR FURTHER READINGS
Jьrgen Herget
Area of investigation and evidence of floods
The headwaters of the river Ob, which form the centre of this study, have their origin close to the border to Mongolia in the AltaiMountains, located in southwestern Siberia (fig. 1). During the Pleistocene, several mountain glaciers reached the valleybottom of the Chuja River, the main tributary of the Katun River, which is the major source of discharge to the Ob River. Near the village of Aktash, the extended glaciers blocked the course of the river Chuja and formed an icedammed lake, which attained a volume of 607 kmі. The water surface of the lake reached a maximum altitude of 2 100 m which is indicated by strandlines and icerafted debris in the intramountainous basins of the Kuray and Chuja.
Fig. 1
At the time of icedam failure, giant floods passed through the valleys of the Chuja and Katun Rivers. The collapse was probably caused by overtopping and rapid incision into the glaciers. The flow transported gravel in suspension [1, 2] deposited this gravel in the form of giant bars at local valley expansions and tributary mouths along the flood route [3] (fig. 2).
At the mouths of some tributaries along the flood route the giant bars blocked valleys and generated secondary lakes, which lasted over different periods of time. These lakes are indicated by lacustrine deposits, such as at Injushka in a tributary valley to the Katun Valley at Inja, where lacustrine sediments are interbedded with suspension gravel indicating repeated outburst floods from the icedammed lake upstream. Other floodrelated features are gravel dunes and deposits of erratic boulders. In combination with exposures of characteristic suspension gravels at the current bedrock valleybottom, the surfaces of the giant bars reveal depths of flow of up to 400 m in Chuja Valley, close to the downstream extent of the icedam. In the wider valley of Katun River, the depth of flow reached about 250 m. At several locations, drawdown levels along the slopes of the giant bars indicate phases of temporary stable water surface levels during the waning flow (fig. 3).
Fig. 2
Fig. 3. Longitudinal profile along Katun and Chuja River downstream of the former icedam with locations and elevations of giant bars and related levels
Runup sediments are deposited in front of local valley obstructions, as a thin layer of suspension gravel on more or less steep bedrock surfaces and can be used to estimate local flow velocity. They indicated a rising water surface level when flow velocity abruptly decreased at the valley obstruction caused by energy conservation [4].
This review focuses on selected attempts to reconstruct the icedammed lake outburst flood, while several previous publications present general reviews based on the state of knowledge available at the time of publication [2, 48]. Other studies focus on special aspects, for example, Parnachov [1] and Carling et al. [3] on the geology of giant bars, Carling [910] on the gravel dunes Reuter et al. [11], on the age of the outburst floods, or Borodavko [12] and Carling et al. [13] on the icedammed lake itself. Recent investigations will provide more detailed information regarding the age of the floods. In the review papers mentioned above, the majority of dated floodrelated features indicate ages between 40 and ~17 ka. Some of the dates are based on luminescence methods, for which methodological problems make older applications doubtful (Carling, personal comm.).
Different methods of reconstruction
Previous estimations of the discharge of the outburst floods were carried out by different researchers. Butvilovski [6 and written comm.] applied a number of different empirical equations to estimate the mean velocity of the flow, based on geometrical parameters of the valleys and basins. Discharges were calculated for crosssectional areas consisting of the entire basins resulting in estimated peak discharges of up to 50 000 000 mіs1. Baker et al. [5] applied the waterlevel calculation software HEC2 to estimate the discharge based on the elevation of floodrelated features. In this approach the assumed discharge values are iteratively optimised until the model calculation of the waterlevel matches the altitude of the waterlevel indicators in the field. They chose a location within the extension of the former icedam and estimated a peak discharge of 18 000 000 mіs1 based on a hydraulic jump located in the modelled reach of the valley. Carling [9] investigated the hydraulic implications of the gravel dunes in Kuray Basin, which was inundated by the icedammed lake. He found a discharge of the magnitude of 750 000 mіs1 for the volume of flow passing over the dunefield which value does not represent the outflow from the Kuray Basin, following the icedam failure. This estimated discharge could not be related to peak discharge of the outburst flood downstream of the icedam as the dune field is located upstream of the valley blockage. Contrary to Butvilovski [6 and written comm.], Carling's attempted to estimate the local discharge of an unspecified stage associated with the decreasing flood in the Kuray Basin appears more plausible. Based on my own investigations of floodrelated features along the Chuja and Katun Valleys, seven independent attempts to gain information about the discharge of the outburst floods were made. The main findings of the study are presented in this paper. Additional details are discussed elsewhere [4].
Conveyance  slope method
Based on the continuity equation and the empirical formula for the estimation of the mean flow velocity in channels by Manning, the discharge Q of the outburst flood can be expressed as Q = A n1 R2/3 S1/2 where A is the crosssection area, n is the roughness coefficient, R is the hydraulic radius and S is the slope of the energy line [14]. Data for the geometric parameters were obtained by surveys of the giant bars along the flood's pathway. The location and the elevations were derived by using GPS as well as calibrated altimeter measurements, or were obtained from 1:50 000 topographic maps (the crosssection area and the hydraulic radius). It is assumed that the slope S of the energy line is, more or less, similar to the slope of the water surface along the direction of the flow (indicated by the slope between the giant bar surfaces along the valleys). The roughness coefficient n was also estimated on the basis of earlier studies. Previous studies on floods of comparable magnitude, such as the drainage of Pleistocene Lake Bonneville [15, 16] and the outburst flood of Lake Missoula [17, 18], the most representative values for the roughness coefficient were determined in the range of 0,04< n <0,07. Hence, for the calculation of the outburst floods in the Altai Mountains, the upper and the limits were considered as extremes.
Six locations of the highest bar surfaces along the Chuja and Katun Valley were selected to estimate peak discharge of the outburst flood. The elevation of the waterlevel indicated by the giant bars, the slope between the bar surfaces, and the recent slope of the rivers are shown in fig. 4. Considering the unknown slope of the flood flow water surface, the surfaces of the giant bars and the current mean river waterlevel are taken as the slope parameter S. Based on the investigations of hydraulic conditions of the flow in steep gradient river flow [19], it is assumed that the subcritical flow conditions are dominant along the outburst flood flow. For the Pleistocene Lake Missoula and Lake Bonneville Floods also, subcritical flow conditions were found to be predominant [16; 17. P. 14]. Therefore, locations and hydraulic parameters associated with flows with Fr >1 (table 1) were excluded from further analysis.
Considering these assumptions, a peak discharge was estimated to be in the range of 10 000 000 mіs1 to 14 000 000 mіs1.
Runup sediments indicating velocity head of flow
Runup sediments form a relatively thin layer of suspension gravel deposited in front of local valley obstructions and are the uppermost deposits related to the outburst flood (fig. 3, fig. 5). They could be considered an indicator of locally elevated water surfaces, due to the physical law of energy conservation. Based on the energy equation after Bernoulli, the total energy H is constant along the channel [14, 20]. As pressure energy is not of relevant for free surface flow, the energy equation can be expressed as
H = (y + z) + (vІ / 2 g),
where, (y + z) is the potential energy, y is the depth of flow, z is the elevation above datum, (vІ/2 g) is the kinetic energy, v is the velocity of flow, and g is the acceleration due to gravity.
While total energy H, elevation above datum z and acceleration of gravity g remain constant, an abrupt decreased in the flow velocity leads to an increased depth of flow, hence, the waterlevel rises on the upstream side of an obstruction.
Fig. 4. Giant bar surfaces and terrace levels along Chuja and Katun valley
Fig. 5
This phenomenon can be observed, for instance, at bridge piers where the elevated water surface on the upstream side of the obstruction occurs and it is used to measure the flow velocity by velocityhead rods [21, 22]. The deposited suspension load indicates the risen water surface level in front of the obstruction, hence provides evidence of the energy head of the peak discharge of the flood. Due to nonuniform distribution of velocity within the open channel flow, an energy coefficient б is introduced to calculate the change in height of the waterlevel from the total loss of kinetic energy [14]. Except in the case of uniform flow, the energy coefficient б is always greater than unity and increases with the steepness of the channel. The equation for the mean velocity head vh in a channel can be written as vh = б (vІ / 2 g).
Depending on the amount of irregularities in channels, the values of б vary between 1 and 2, rarely exceeding 4,70 [23]. Available literature indicates that the average value is about 1,3. However, experience shows that the effect of an increased velocity head vh by б is small compared with other uncertainties. Hence, the energy coefficient is frequently considered to be equal to unity, as the data are sparse and are found to be inconsistent [14, 24, 25].
Assuming that the runup sediments and uppermost giant bar surfaces located nearby, were generated during the same stage of the flood, flow velocity v can be determined by taking velocity vh as the difference in elevations of the runup sediments and the giant bar surfaces and assuming б = 1,3.
While the elevated waterlevel upstream of obstructions is indicated by suspension gravel being deposited in the uppermost position, the depth of the flow above the giant bar surfaces is unknown. The maximum value for velocity head vh is given by the assumption that the undisturbed waterlevel of the flood was identical to that of the giant bar's surfaces. Therefore, calculated flow velocity also reaches maximum values. To check the sensitivity of this parameter, an undisturbed waterlevel half way between the giant bar and the related runup sediments is considered.
Fig. 6. Locations of runup sediments and related bar surfaces indicated velocity heads along Chuja and Katun valley (The maximum possible velocity head vh max is limited by the difference in elevation of the runup sediments and the bar surfaces. Its value is given in the figure.
The velocity head vh assumed for the assumed water level above the bar's surfaces is half of this value) For several locations along the Chuja and Katun Valley (fig. 3), velocities of the flow can be estimated on the basis of the relation between runup sediments and related giant bar surfaces nearby. It is important to note, that the runup sediments are deposited only during the peak discharge due the elevated waterlevel on the upstream side of local obstructions, but also during lower stages of the decreasing flood. Figure 6 shows values of local velocity heads vh for different locations at two elevations considered for the undisturbed waterlevels. A maximum value of vh = 137 m is indicated by field data.
Figure 7 shows the calculated mean flow velocities based on the velocity heads given in Figure 6 and an assumed value of the energy coefficient б = 1,3. Velocity heads would decrease for higher waterlevels, hence the undisturbed flow velocity would also decrease. Flow velocities estimated from both waterlevels are within the expected range considering the magnitude of the flood as estimated previously by uniform flow calculations (table 1). For vh max flow velocities vary between 24 and 45 ms1, for the higher waterlevel with vh assumed 17 ms1 < v < 32 ms1 are calculated.
Fig. 7. Estimated mean undisturbed flow velocities derived from velocity heads indicated by runup sediments for two local water levels along Chuja and Katun valley
These current conditions of higher kinetic energy (= higher flow velocity) can be expressed by Froude number (Fr). For lower waterlevels, Fr was found to be in the range of 0.48 and 1.17, whereas for the higher waterlevel the number was determined as 0,32 < Fr < 0,69. The broad bedload terrace filling the valleybottom is assumed to be deposited up to the current thickness later during the recession limb of the flood and also during the Holocene. Hence, the depth of flow used to calculate Froude number reaches from the current river channels (bedrock) up to the waterlevel of the flood stage, forming the runup sediments.
To obtain subcritical conditions of flow, the energy coefficient б must be modified up to a value of 1,8. This gave Froude numbers between 0,40 and 0,99 for the lower waterlevel, and between 0,28 and 0,58 for the higher one. For this value of the energy coefficient, flow velocities vary between 21 and 39 ms1 for the lower and between 15 and 27 ms1 for the higher water surface. None of these values appear critical as all the data are within the expected magnitude. This sensitivity analysis reveals that, considering the scale of the flood, the calculated flow velocities and Froude numbers are less sensitive to modifications of the energy coefficient б. Hence, the value for б is not a critical in the estimation of largemagnitude floods. Nevertheless, supercritical conditions of flow might have occurred at locations of valley obstructions.
The calculation of discharge Q based on the continuity equation Q = v. A, for estimating the flow velocities v is occasionally a difficult task, because of the problems related with the calculation of the crosssection area A. As some of the runup sediments are located at prominent valley obstructions and extended into ineffective areas of flow, the incident crosssection area cannot be determined. Similar problems occur at other locations, where the runup sediments are deposited within the mouth of tributaries, and where flow separation phenomena must have occurred. Here too, related crosssection areas cannot be accurately estimated.
The calculated discharges (table 2) are within the expected range for peak discharge of the outburst flood. The relatively small difference between discharges of higher and of lower waterlevel elevations is remarkable. The lower waterlevel with reduced crosssection area is compensated by higher flow velocities to reach the height of the runup sediments from a lower water surface elevation. Only occasionally, the flow velocity increases by an amount that results in a higher discharge for the lower water surface elevation. This effect might be enhanced by steep valleyslopes, which result into lower decrease of crosssection area at a decreased waterlevel than for valley sections with gentle slopes. On the other hand, the differences are in many cases within the range of accuracy. Therefore, the need for assuming the height of the water surface seems to be less crucial for the estimation of discharge.
The method to estimate flow velocities for discharges from investigated data of the velocity head appears to be a valid approach to determine the magnitude of the flood. This is due to the fact that the results are within the expected range of the peak flood associated with the giant bars and runup sediments, to give a peak discharge of about 10 000 000 mіs1. The assumed parameters of energy coefficient б and the surface of the undisturbed waterlevel appear to be less critical in the estimation of the discharge. Occasionally, depending on the local topography with extended areas of ineffective flow in tributary valleys, the determination of the related crosssection area is complicated and precludes the accurate calculation of discharge from the estimated velocity of flow. Due to the nature of the obstructions  such as the ridges at tributary confluences causing the waterlevel to rise  the valley crosssections show a significant modification in width which varies the crosssectional area of flow for a given water surface elevation by several hundred percent.
Table 2 Discharges calculated by velocity heads for selected runup sediment locations
Location valleykm 
vh (assumed) : altitude of waterlevel flow velocity v crosssection area A discharge Q 
vh (max) : altitude of waterlevel flow velocity v crosssection area A discharge Q 

KezekDzhala km 9,35 
850 m v = 26 m/s A = 222 000 mІQ = 5,8Ч106 mі/s 
806 m v = 36 m/s A = 133 000 mІQ = 4,7Ч106 mі/s 

Anijakh km 9,65 
843 m v = 22 m/s A = 205 000 mІQ = 4,5Ч106 mі/s 
810 m v = 31 m/s A = 139 000 mІQ = 4,3Ч106 mі/s 

Little Jaloman Village km 14,6 
908 m v = 29 m/s A = 306 000 mІQ = 8,9Ч106 mі/s 
850 m v = 42 m/s A = 215 000 mІQ = 9,0Ч106 mі/s 

Giant bar between Inja and Little Jaloman Village km 20,65 
949 m v = 32 m/s A = 384 000 mІQ = 12,3Ч106 mі/s 
880 m v = 45 m/s A = 225 000 mІQ = 10,1Ч106 mі/s 

Lower Chuja Valley km 35,7 
1040 m v = 28 m/s A = 362 000 mІQ = 10,1Ч106 mі/s 
989 m v = 39 m/s A = 256 000 mІQ = 10,0Ч106 mі/s 

Lower YalbakTash km 39 
1057 m v = 32 m/s A = 300 000 mІQ = 9,6Ч106 mі/s 
990 m v = 45 m/s A = 186 000 mІQ = 8,4Ч106 mі/s 

Upper YalbakTash km 41,4 
1077 m v = 31 m/s A = 342 000 mІQ = 10,6Ч106 mі/s 
1015 m v = 43 m/s A = 240 000 mІQ = 10,3Ч106 mі/s 

The roughness of the surface of the valleyslopes might lead to an underestimation of the kinetic energy of the flow. Parts of the energy necessary to raise the waterlevel at local obstructions might be dissipated to overcome the roughness of the obstruction's surface. This effect, in turn, might lead to deposition of the suspension load at lower elevation, implying that the extent of the rise in the water surface is determined by the surface roughness, which cannot be easily quantified. Validation of the method is desirable to estimate it accurately. Wilm and Storey [21] found 10% errors in the flow velocity estimated on the basis of velocityhead rod method. However, this value was derived for flows of much lower magnitude and partly from controlled flume experiments.
Regression drained lake volumepeak discharge
A wellestablished method to estimate the peak discharge of modern icedammed lake outburst floods is the empirical regression of the drained volume of water and the peak discharge. Even though there is no physical relationship between these parameters and significant differences are expected occur along the floods' pathway due to valley topography and slope and retention effects, the relationships work remarkably well. Several attempts have been made to derive regression equations, in spite of the criticism related to the selection of the events, differences in outflow mechanism and the triggering mechanism of the outbursts. The plots given in figure 7 show that the regression lines obtained by various workers are comparable.
The application of such relationships for largescale events is limited. Only Clague and Mathews [26] considered estimations of peak discharge of the outburst flood from the Pleistocene icedammed Lake Missoula. Therefore, a new regression was derived for largescale events. An arbitrarily threshold drained volume of 1 kmі (= 109 mі) was chosen to identify largescale events and separate them from the smaller events. Eighteen instances of icedammed lake outburst floods from Pleistocene and modern times were taken from the literature to derive the relationship. The instances are listed in Table 3. Note that wellknown examples that were either not exclusively related to icedams and those that were related to some kind of uncertainties, were not considered. For a drained volume V given in kmі a regression equation of Q = 6 645 V0,98 was determined, Q is given in mіs1 and rІ = 0,93. ice pond lake flood
Assuming that the entire volume of the icedammed lake in the Kuray and Chuja Basin, with a volume of 607 kmі (for a lake level at 2100 m), was drained during the outburst flood, the estimated peak discharge is 3 500 000 mіs1.
A comparison of this estimate with the earlier estimates indicates that regression equation underestimates the peak discharge. A closer examination of the examples given in table 3, reveals that Lake Missoula Flood is the only outburst flood with an estimated drained volumes of more than 100 kmі. Even after taking into account the fact that it is incorrect to include two different values for the estimated peak discharge from Lake Missoula (based on different estimation attempts) and that the complex topography of the flood's pathway have a profound ponding effect, the amount of data to develop a regression equation is insufficient. Hence, it could be concluded that such an estimate of the peak discharge is not reliable.
Table 3 Large icedammed lake outburst floods in historic and prehistoric times
Location 
Date 
Lake volume, kmі 
Peak discharge, mі/s 
Reference 

Missoula,Montana USA 
Pleistocene 
2,184 
17 000 000 
O'Connor and Baker (1992) 

Missoula,Montana USA 
Pleistocene 
2,184 
10 000 000 
Baker and Costa (1987) 

Nedre Glеmsjш,Norway 
Pleistocene 
99 
170 000 
Berthling and Sollid (1999) 

Lake Alsek,Canada 
Holocene 
30 
470 000 
after Clarke from Clague and Evans (1994) 

Kjцlur, Iceland 
9500 BP 
25 
200 000 
Tуmasson (2002) 

Hubbart Glacier,Alaska USA 
1986 
5,3 
104 500 
Mayo (1987), Seitz et al. (1986),Mayo (1989) 

Lake Alsek,Canada 
ca. 1850 
4,7 
30 000 
after Clarke from Clague and Evans (1994) 

Lake Elk, Canada 
Pleistocene 
4 
20 000 
Clague (1973) 

Lake George, Alaska USA 
1958 
2,2 
10 160 
after US Geological Survey from Fahnestock and Bradley (1973) 

Lake George, Alaska USA 
1961 
1,73 
10 050 
after US Geological Survey from Fahnestock and Bradley (1973) 

Graenalon, Iceland 
1935 
1,5 
5 800 
Thorarinsson (1939) 

Graenalon, Iceland 
1939 
1,5 
5 000 
Thorarinsson (1939) 

Lake Batal,Himalaya 
Pleistocene 
1,496 
24 000 
Coxon et al. (1996) 

Lake George, Alaska USA 
1960 
1,48 
9 280 
after US Geological Survey from Fahnestock and Bradley (1973) 

Van Cleve Lake, Alaska USA 
1992 
1,4 
4 500 
after Brabets from Walder and Costa (1996) 

Chong Kumdan(Shyo), India 
1929 
1,35 
22 650 
after Gunn or Mason et al. From Hewitt (1982) 

Lake George, Alaska USA 
1959 
1,11 
6 310 
after US Geological Survey from Fahnestock and Bradley (1973) 

Lake George, Alaska USA 
1965 
1,11 
6 680 
after US Geological Survey from Fahnestock and Bradley (1973) 

Further attempts
Additional attempts to estimate discharge related to the outburst floods in Altai Mountains are applied. Due to the limitation of the space, only a short review of the attempts and related problems is given.
The waterlevel calculation software HECRAS (Hydrologic Engineering Center of the US Army Corps of Engineers  River Analysis System) [27, 28, http://www.hec.usace.army.mil/software/hecras/  27th Sept. 2011] is applied to estimate discharge based on the palaeostage, represented by giant bars' surfaces [29]. The software package allows calculations of discharge and water surface level for onedimensional steady and unsteady flow under sub and supercritical conditions. Energy losses are evaluated by friction (Manning equation) and contraction or expansion by coefficients multiplied by the change in velocity head. The momentum equation is applied in situations where the water surface profile is rapidly varied, including mixed flow regime calculations with hydraulic jumps, hydraulics of bridges and profiles at river confluences.
Fig. 8. Selected published relations of drained icedammed lake volume V and peak discharge Q
To accomplish the presumption of gradually varied flow between the 244 crosssections derived from detailed topographic maps, the modelled reaches of the valleys must be relatively homogenous. Hence, three units were separated and modelled independently, because obvious changes in valley characteristics require this: the broad and less steep Katun Valley up to the confluence with the Chuja River (km 030; 2,2‰), the straight and relatively steeper lower part of the Chuja Valley (km 3044; 6,0‰) and the winding upper part of the Chuja Valley (km 4485; 6,0‰). The area of the confluence itself is left out for modelling as the conditions of flow change abruptly in this section. Considering, the available literature the following loss coefficients of roughness, contraction and expansion were chosen: n = 0,04, c = 0,2 and e = 0,4. Regarding runup sediments as indicator of maximum waterlevel and giant bar surfaces as minimum elevations, peak discharge for subcritical flow conditions can be estimated as between 8Ч106 mіs1 and 12Ч106 mіs1 for all subreaches. Depth of flow varies between 280 and 400 m, with mean flow velocity of 937 ms1. Froude numbers of 0,200,85 confirm subcritical flow conditions along the entire valley.
In the meantime, improvements of the HECRAS software allow unsteady flow simulations which were not possible due to numerical instabilities with previous versions. Rudoy and Zemtsov [30] modelled the unsteady outburst flood principally, but the modelled flood levels do not reach the palaeostage indicators like the giant bar surfaces and peak discharge is less than found by other approaches. Carling et al. [31] used alternative software to model the outburst flood as unsteady, one and twodimensional flow and found peak discharges confirming previous magnitudes.
Between the villages of Inja and Little Jaloman in the Katun Valley (fig. 2), boulders, which were obviously transported by the flood are deposited. The mean diameter of the five largest boulders is 11,3 m, while the largest one has a diameter of 13,5 m. To interpret the competence of the flow indicated by the dimension of the boulders, is an apparent attempt but limited by the dimension of the transported boulders. As described previously, for example, by Baker [32], transport of boulders with diameters of more than about 23 m is related to macroturbulence effects, which create an uplift force due to modification of the current system around the boulder. Due to the magnitude of the related flow conditions, hydraulic interpretation of this effect of nonlinearity is not well understood [33] and theorybased attempts are not available.
In spite of these problems, a rough estimation of the flow conditions required to transport these large boulders, based on relationship derived for smaller dimension boulder deposits available in the literature, result in flow velocities of about 20 ms1 with related discharges of 3 000 000 mіs1. As these values could be overestimations, due to the macroturbulence effects mentioned above, the transport of the boulders may not be necessarily related to the peak discharge of the outburst flood. Even the magnitude of the required discharge is significantly smaller than peak discharge estimated by other attempts. Hence, the transport of the boulder occurred at an unspecified stage associated with falling flood stage.
Fluvial gravel dunes are formed as bedform features at the previous lake bottom and along the flood's pathway down to the opening of the Katun Valley towards the western Siberian Plain, where traces of the event fade out. As mentioned above, the sedimentology of several gravel dune fields is subject of detailed investigations by Carling [9, 10], who also carried out hydraulic interpretations, which result in estimated flow velocities of 1,58 ms1 and related discharges of 20 000750 000 mіs1 over the span of the dune fields. The largest gravel dunes were found at the eastern margin of the Chuja Basin with heights up to 23 m and wavelengths up to 320 m.
Based on a) empirical relationships for the geometry of fluvial dunes; b) the depth of flow in combination with estimations of the flow competence indicated by the gravel sizes of the investigated dunes, the estimation of flow conditions forming the residual dune structures can be made. Depending up on the different sizes of the dunes, mean flow velocities of 4,810 ms1 and discharges in the range of 10 000 1 340 000 mіs1 were determined. As graveldunes grow during the rising stage of a flood and decrease at the end, the residual dimensions cannot be related to the peak discharge but again to an unspecified stage at the termination of the outflow.
Similar problems are encountered in case of the obstacle marks, found in different sizes at the boulder field between Inja and Little Jaloman, and around a bedrock hill near the village of ChaganUzun at the eastern margin of the Chuja Basin. As scour holes refill during the decreasing flood [34], they are dynamic features like graveldunes.
The boulder between Inja and Little Jaloman is about 10,8 m wide with a length of 13,8 m and a maximum height of 3,1 m above the surrounding surface. The scour hole is filled by loess deposits blown in after the outburst floods, but can easily be separated from the gravelpaved surface around by changes in the vegetation cover. The width of the depression is 24,3 m with a maximum observed depth of 2,3 m. The scour hole around the bedrock hill near ChaganUzun has a measured depth of 8,1 m. Length along the direction of flow is 91,5 m and the width is nearly 400 m, while the hill itself is 50 m high. The maximum elevation of the deposits on leeward side of the hill is about 4 m above the area, while its counterpart at the boulder reaches the height is 0,5 m.
Calculations of scour hole dimensions and dynamics are a classical engineering task as bridges and oil platforms, for example, are in potential danger to fail by scour around the piers. Unfortunately experiences gained by related investigations cannot be transferred, as the current pattern is different, because piers engineering structures are not usually submerged by the flow. Similarly, the results from flume experiments are not applicable directly, because earlier investigations derived mainly qualitative results, while current studies are mainly carried out on a fluid mechanical base in flumes under welldefined environmental conditions.
Finally, the estimations of flow conditions based on the characteristics of the obstacle marks are based on minimum values of the depth of flow derived from the scour hole dimensions and flow velocity assessments related to the flow competence indicated by grain sizes of sediments at the outer scour hole walls. The scour hole around the boulder near Inja seems to be finally formed by a current with a velocity of 1,4 ms1 at a discharge of about 800 mіs1. The larger depression around the bedrock hill contains evidence of a flow of about 11 500 mіs1 with flow velocity of approximately 1,7 ms1. Obviously, these discharges represent the very final stage of the flood, and had previously reduced the scour holes dimension by refilling with transported sediments.
Results
The different attempts to estimate flow conditions of the Pleistocene outburst floods in Altai Mountains result in data for flow velocity and discharge of the flood's peak and unspecified stages of the waning flood. Based on the conservative estimate of the peak discharge to be about 10 000 000 mіs1 and assuming that the entire lake basin, with a volume of 607 kmі, drained during the outburst a hydrograph can be generated. The hydrograph indicates a flood that did not last for more than 23 days (fig. 9). The duration of the flood was determined by the integral curve based on the limitations of the drained volume (= integral) and the peak discharge. Advanced modelling of unsteady twodimensional flow confirms this magnitude. From a negligible baseflow, an abrupt rise of the hydrograph represents the outburst flood wave, which reached peak discharge values almost immediately. It should be noted that the hypothetical hydrograph does not consider ponding effects along the flood's pathway. Consequently, a tendency towards underestimation of the duration of the flood is inherent.
Fig. 9. Theoretical hydrograph of the icedammed lake outburst flood based on peak discharge and drained volume
Compared with the earlier estimates of the peak discharge of the outburst flood [5] the new calculations give a slightly higher discharge magnitude. Considering the fact that the present study is based on larger number of palaeostage indicators and also includes the calculation of the discharge downstream of the former icedam, the new estimations appear reliable and contain several elements of conservative assessment.
Acknowledgements
The author appreciates fruitful cooperation with Paul Carling and Pavel Borodavko for more than 10 years, scientific cooperators became friends during this time. Peter Martini, Sergei Parnachov and Heike Agatz provided valuable comments and friendly cooperation during fieldwork. A fruitful cooperation with Greg Balco, Dr. Barbara Mauz and Anne Reuter produced successful dating of several samples of flood related features. Kirsten Hennrich kindly improved the language of the contribution. The studies are financially supported by Deutsche Forschungsgemeinschaft DFG, Tomsk State University, British Research Council and the Department of Geography, RuhrUniversitдt Bochum.
References
1. Parnachov S.V. Geology of giant Altay terraces. Tomsk : Tomsk Polytechnic University, 1999, 137 p. (in Russian)
2. Carling P.A., Kirkbride A.D., Parnachov S.V., Borodavko P.S., Berger G.W. Late Quaternary catastrophic flooding in the Altai Mountains of southcentral Siberia: a Synoptic overview and an introduction to the flood deposit sedimentology / tds. P.I. Martini, V.R. Baker, G. Garson // Flood and megaflood processes and deposits: resend and ancient examples : Int. Ass. of Sedimentologists. Oxford, England, 2002. Spec. Publ. 32. P. 1735.
3. Carling P.A., Martini P., Herget J., Borodavko P., Parnachov S. Megaflood sedimentary valley fill: Altai Mountains, Siberia / tds. D. Burr, V.R. Baker, P.A. Carling // Megaflooding on Earth and Mars. UK : Cambridge University Press, 2009. P. 243266.
4. Herget J. Reconstruction of icedammed lake outburst floods in the AltaiMountains, Siberia // Geological Society of America. Special Publication. 2005. 386 p.
5. Baker V.R., Benito G. and Rudoy A.N. Palaeohydrology of Late Pleistocene superflooding, Altay Mountains, Siberia // Science. 1993. Vol. 259. P. 348350.
6. Butvilovsky V.V. Last glaciation and Holocene palaeogeography of the Altai: a cataclysmic model. Tomsk : Tomsk State University, 1993, 252 p. (in Russian)
7. Rudoy A.N. Mountain icedammed lakes of Southern Siberia and their influence on the development and regime of the intracontinental runoff systems of North Asia in the late Pleistocene / eds. G. Benito, V.R. Baker, K.J. Gregory // Palaeohydrology and environmental change. Chichester: John Wiley & Sons Ltd., 1998. P. 215234.
8. Rudoy A.N., Baker V.R. Sedimentary effects of cataclysmic late Pleistocene glacial outburst flooding, Altay Mountains, Siberia // Sedimentary Geology. 1993. Vol. 85. P. 5362.
9. Carling P.A. A preliminary palaeohydraulic model applied to late Quaternary gravel dunes: Altai Mountains, Siberia / eds. J. Branson, A.G. Brown, K.J. Gregory // Global Continental Changes: the context of palaeohydrology // Geological Society. Special Publication. 1996. Vol. 115. P. 165179.
10. Carling P.A. Morphology, sedimentology and palaeohydraulic significance of large gravel dunes, Altai Mountains, Siberia // Sedimentology. 1996. Vol. 43. P. 647664.
11. Reuter A., Herget J., IvyOchs S., Borodavko P., Kubik P.W., Heine K. Constraining the timing of the most recent cataclysmic flood event from icedammed lakes in the Russian AltaiMountains, Siberia, using cosmogenic in situ 10Be // Geology. 2006. Vol. 34. P. 913916.
12. Borodavko P.S. LatePleistocene glacial lakes in southeastern Altai // Voprosy Geografii Sibili. 1999. Vol. 23. P. 8794 (in Russian).
13. Carling P.A., Knaapen M., Borodavko P., Herget J. et al. Palaeoshorelines of glacial lake KurayChuja, southcentral Siberia: form, sediments and process / tds. I.P. Martini, H.M. French, A. Pйrez Alberti // IceMarginal and Periglacial Processes and Sediments. London: Geological Society, Special Publications. 2011. Vol. 354. P. 111128.
14. Chow V.T. Openchannel hydraulics. Tokyo : McGrawHill, 1959. 680 p.
15. Jarrett R.D., Malde H.E. Palaeodischarge of the late Pleistocene Bonneville Flood, Snake River, Idaho, computed from new evidence // Geological Society of America Bulletin. 1987. Vol. 99. P. 127134.
16. O'Connor J.E. Hydrology, hydraulics, and geomorphology of the Bonneville Flood // Geological Society of America. Special Paper. 1993. Vol. 274. 83 p.
17. O'Connor J.E., Baker V.R. Magnitudes and implications of peak discharges from glacial Lake Missoula // Geological Society of America Bulletin. 1992. Vol. 104. P. 267279.
18. Benito G. Energy expenditure and geomorphic work of the cataclysmic Missoula flooding in the Columbia River gorge, USA // Earth Surface Processes and Landforms. 1997. Vol. 22. P. 457472.
19. Jarrett R.D. Hydraulics of highgradient streams // Journal of Hydraulic Engineering. 1984. Vol. 110. P. 15191539.
20. Chanson H. The hydraulics of open channel flow. London : Arnold, 1999. 495 p.
21. Wilm H.G., Storey H.C. Velocityhead rod calibrated for measuring stream flow // Civil Engineering. 1944. Vol. 14. P. 475476.
22. Peterson M., Cromwell C.F. Measuring irrigation water in a ditch, stream or reservoir. Columbia, University of Missouri, Department of Agricultural Engineering, Agricultural Publication, 1993. Vol. G01681. 7 p.
23. Hulsing H., Smith W., Cobb E.D. Velocityhead coefficients and open channels // U.S.
Geological Survey WaterSupply Paper. 1966. Vol. 1869C. 21 p.
24. Henderson F.M. Open channel flow. N.Y. : Macmillan, 1966. 522 p.
25. Barnes H.H., Davidian J. Indirect methods / ed. R.W. Herschy // Hydrometry  principles and practises. Chichester : Wiley, 1978. P. 149204.
26. Clague J.J., Mathews W.H. The magnitude of jцkulhlaups // Journal of Glaciology. 1973. Vol. 12. P. 501504.
27. Brunner G.W. HECRAS River Analysis System  User's manual, Version 3.0. Davis, California, US Army Corps of Engineers Hydrologic Engineering Center (HEC), 2001. 320 p.
28. Brunner G. WHECRAS River Analysis System  Hydraulic reference manual. Davis, California, US Army Corps of Engineers Hydrologic Engineering Center (HEC), 2001. 262 p.
29. Herget J., Agatz H. Modelling icedammed lake outburst floods in AltaiMountains (Siberia) with HECRAS / eds. V.R. Thorndycraft, G. Benito, M. Barriendos, and M.C. Llasat // Palaeofloods, historical data & climatic variability: Applications in flood risk assessment. Madrid, Proceedings of the International PHEFRA Workshop held in Barcelona, 1619th October 2002. Centro de Ciencias Medioambientales, 2003. P. 177181.
30. Rudoy A.N., Zemtsov V.A. Modelling of hydraulic parameters of diluvial floods of the ChuyaKuray icedammed lake in the Quaternary // Ice and Snow. 2010. Vol. 109. P. 111 118 (in Russian).
31. Carling P.A., Villanueva I., Herget J. et al. Unsteady 1D and 2D hydraulic models with ice dam break for Quaternary megaflood, Altai Mountains, southern Siberia // Global and Planetary Change. 2010. Vol. 70. P. 2434.
32. Baker V.R. Palaeohydrology and sedimentology of Lake Missoula flooding in eastern Washington // Geological Society of America. Special Paper. 1973. Vol. 144. 79 p.
33. Jackson R.G. Sedimentological and fluiddynamic implications of the turbulent bursting phenomenon in geophysical flows // Journal of Fluid Mechanics. 1976. Vol. 77. P. 531560. 34. Hoffmanns G.J.C.M., Verheij H.J. Scour manual. Rotterdam : Balkema, 1997. 205 p. Received November 13, 2011
Abstract
This review focuses on selected attempts to reconstruct the icedammed lake outburst flood in the valleybottom of the Chuja River, the main tributary of the Katun River, which is the major source of discharge to the Ob River, while several previous publications present general reviews based on the state of knowledge available at the time of publication. Other studies focus on special aspects, for example, Parnachov and Carling et al. on the geology of giant bars, Carling on the gravel dunes, Reuter et al., on the age of the outburst floods, or Borodavko and Carling et al. on the icedammed lake itself. Recent investigations will provide more detailed information regarding the age of the floods. In the review papers mentioned above, the majority of dated floodrelated features indicate ages between 40 and ~17 ka. Some of the dates are based on luminescence methods, for which methodological problems make older applications doubtful (Carling, personal comm.).
The different attempts to estimate flow conditions of the Pleistocene outburst floods in Altai Mountains result in data for flow velocity and discharge of the flood's peak and unspecified stages of the waning flood. Based on the conservative estimate of the peak discharge to be about 10 000 000 mіs1 and assuming that the entire lake basin, with a volume of 607 kmі, drained during the outburst a hydrograph can be generated. The hydrograph indicates a flood that did not last for more than 23 days. The duration of the flood was determined by the integral curve based on the limitations of the drained volume (= integral) and the peak discharge. Advanced modelling of unsteady twodimensional flow confirms this magnitude. From a negligible baseflow, an abrupt rise of the hydrograph represents the outburst flood wave, which reached peak discharge values almost immediately. It should be noted that the hypothetical hydrograph does not consider ponding effects along the flood's pathway. Consequently, a tendency towards underestimation of the duration of the flood is inherent.
Compared with the earlier estimates of the peak discharge of the outburst flood the new calculations give a slightly higher discharge magnitude. Considering the fact that the present study is based on larger number of palaeostage indicators and also includes the calculation of the discharge downstream of the former icedam, the new estimations appear reliable and contain several elements of conservative assessment.
Key words: icedammed lake; outburst flood
Размещено на Allbest.ru
...Подобные документы
A Phase Diagram Approach to Microemulsions. Percolating Phenomenon Microemulsions: Effect of External Entity. Influethylene of Polyethylene Glycols and Polyethylene Glycol Dimethyl Ethers upon the Internal Dynamics of Water in Oil Microemulsions.
презентация [8,7 M], добавлен 07.12.2010Charles Dickens and his autobiography. Charles Dickens as the most popular and productive novelist combined social Critism with comedy and sentiment, creating a tone that the world identifies as Victorian lake chancer and Shakespeare before him.
курсовая работа [50,5 K], добавлен 21.07.2009When we have time for leisure, we usually need something that can amuse and entertain us. Some people find that collecting stamps, badges, model cars, planes or ships, bottles, or antiques are relaxing hobbies. Free time is organized in many schools.
сочинение [5,2 K], добавлен 04.02.2009The results of theoretical analysis and computer simulation of the amplitude and phase errors of the narrowband signal. Vector representation of input and output signals. Standard deviation of the phase. Probability distribution laws of the phase error.
реферат [469,7 K], добавлен 06.04.2011Introduction to geographical location, population size, state of the industry, energy resources, transportation infrastructure in Alaska. Study location, swimming pools, demographics, and the main attractions of California  one of the states of America.
презентация [387,4 K], добавлен 05.11.2010The main theories in the field of human origin, their basic content and direction of research. Basic stages of human development from the primitive to the modern form of the form. Character change erectus skeleton human time frame of the process.
презентация [614,1 K], добавлен 26.09.2014Comparative analysis of acronyms in English business registers: spoken, fiction, magazine, newspaper, nonacademic, misc. Productivity acronyms as the most difficult problem in translation. The frequency of acronym formation in British National Corpus.
курсовая работа [145,5 K], добавлен 01.03.2015Air, water and soil as necessity for existence of all living things. Importance of solving the environmental problems that endanger people's lives. Water and air pollution. Pesticides, rubbish and poisonbeware. Reduction of pollution. Drainage systems.
доклад [27,1 K], добавлен 08.01.2011Definition and general characteristics of the wordgroup. Study of classification and semantic properties of the data units of speech. Characteristics of motivated and unmotivated wordgroups; as well as the characteristics of idiomatic phrases.
реферат [49,3 K], добавлен 30.11.2015Contextual and functional features of the passive forms of grammar in English. Description of the rules of the time in the passive voice. Principles of their translation into Russian. The study of grammatical semantics combinations to be + Participle II.
курсовая работа [51,9 K], добавлен 26.03.2011