1. Introduction
Occasionally, the water stored in a dam reservoir can be released suddenly due to a dam burst, creating a large-scale burst flood [
1,
2,
3,
4,
5]. These large-scale outburst floods have significant entrainment effects on riverbeds [
6]. For example, they can erode and scour both bedrock and soils and lead to alterations in the topography of the original river course. These effects can even cause the disappearance of certain landforms present in the original river channel. In addition, these outburst floods can transport sediments from upstream, which are mainly bed loads, which, once they reach the terraces of riverbanks, will accumulate and change the topography and landforms of the riverbanks [
7,
8,
9].
A boulder bar is an accumulated river landform that appears in the river channel and along the riverbanks. The sizes or shapes of boulder bars can be changed by erosion or deposition around them due to outburst floods [
10]. Therefore, a boulder bar is a reference that reflects the impact of outburst floods on the topography and geomorphology of river courses. The sizes of boulder bars are directly affected by the velocities and water levels of floods; that is, boulder bars can be altered slightly by outburst floods with low velocities and low water levels. Meanwhile, under conditions of rapid currents and high-water levels, boulder bars will be strongly eroded and even destroyed due to the action of bursting floods. Furthermore, new boulder bars may be formed downstream due to the same action [
11]. In addition, the sizes of the sediment particles that are carried by burst floods also affect the sizes of boulder bars. In general, the transported sediments with coarser particles have a dominant effect on the shapes and sizes of boulder bars, while smaller particles are more likely to settle and accumulate around the boulder bars, increasing their sizes [
12,
13,
14].
Floods that are caused by natural dam failures are a common type of flooding in river basins. Many cases have shown that outburst floods triggered by natural dam bursts can change the sizes and shapes of boulder bars downstream [
10,
11]. For example, Wu et al. [
11] selected a research area located 17 km downstream of Yigong Barrier Lake and found that outburst floods were responsible for a reduction in the area of boulder bars in the study area. However, these studies analyzed only the areas, lengths and widths of boulder bars, and the quantitative relationships among these different parameters were not examined. In particular, there is a lack of comparative analyses of the sizes of boulder bars before and after natural dam failures. In addition, the variation modes of boulder bars after a natural failure are also an interesting topic, but they have not been reported in the literature.
The topography of fluvial erosion has often been cited as an example of self-similarity in nature [
15,
16]. This indicates that riverine geometric structures conserve self-similar shapes for long-term fluviation. Transient river dynamics, such as a sudden and dominant impact of outburst floods on river systems, followed by an accumulation of sediments or adjustment of bed topography, typically exhibit both translation and dispersion [
17]. Whether self-similarity still exists in the stream bed topography after this transient process is still unclear. Hence, the objectives of our research are to show how a natural dam failure influences downstream riverine morphological features and whether the bedform features take self-similar shapes, as found in previous research. For this, a 173 km section of the downstream river channel of the Sedongpu natural dam was taken as the research area.
The Sedongpu natural dam occurred in the Sedongpu Basin, with a landslide volume of 3000 × 10
4 m
3, on the 29 October 2018. The fundamental cause of this event was the collapse of a glacier. The collapse occurred approximately 6 km upstream in the Sedongpu Basin at an elevation of 4070 m [
18]. During the collapse process, part of the glacier detached and carried rocks and debris from the sides of the valley and the base of the glacier. Carrying ice and debris from the glacier, the mixture rushed down the valley. The snow and debris moved rapidly towards the bottom of the valley, traveling a long distance before entering the Yarlung Tsangpo River and colliding with the right bank, ultimately blocking the Yarlung Tsangpo River and forming a dam. The dam was approximately 2.5 km long, with a width ranging from 415 to 620 m and an area of approximately 1.29 km
2, with a volume ranging from 33 to 45 × 10
6 m
3 [
18]. The Shedongpu natural dam was mainly composed of gravel and fine soil, with a small amount of snow, and a few boulders. The backwater in the upper reaches of the dam seriously threatened the road water conservancy facilities, electric power communication facilities and cultivated land in Gara Village, Blunt Village and Chibai Village in Tibet. The ratio of soil to rock was approximately 8:2 [
19]. The natural dam breached in approximately 2 days [
18]. The dam failure seriously damaged downstream Medog County and the Medog Yarang Hydropower Station. And 16,600 people were affected, 0.34 km
2 of arable land was affected, 7 km of road was damaged, and the preliminary estimate of economic loss was more than CNY 3 × 10
8.
In this study, the sizes of boulder bars, such as their lengths, widths, areas and perimeters, were compared before and after the failure of the Sedongpu natural dam combined with remote sensing images of the study area. And the number variation of boulder bars after the dam failure was studied. In addition, the evolutionary pattern of boulder bars in the study area was investigated and summarized. Furthermore, the most important finding is self-similarity of riverine morphological features before and after the natural dam failure. This finding supplement previous observations, as the dam failure had a sudden and significant impact on the river system.
3. Results
3.1. Morphological Characteristics of Boulder Bars
Before the dam failure, the lengths of the boulder bars ranged from 62.7 m to 2195 m, the widths of the boulder bars ranged from 15.69 m to 238 m, the perimeters of the boulder bars ranged from 163 m to 4450 m, the areas of the boulder bars ranged from 1211 m2 to 228,447 m2, and the flatness of the boulder bars ranged from 3.5% to 97.7%. After the dam failure, the lengths of the boulder bars ranged from 77.23 m to 2357 m, the widths of the boulder bars ranged from 24.5 m to 203 m, the perimeters of the boulder bars ranged from 183.88 m to 4904 m, the areas of the boulder bars ranged from 1994.2 m2 to 335,213 m2, and the flatness of the boulder bars ranged from 36.6% to 95.8%.
The shapes of the boulder bars in the study area can be divided into three types: (1) The sickle type: the boulder bar bodies bend towards the convex bank of the river, the boulder bar heads (upstream of the boulder bar) are slightly blunt, and the boulder bar tails (downstream of the boulder bar) are sharpened like a sickle. (2) The bamboo leaf type: the edges of the boulder bar are smooth, and the boulder bar heads are blunt, short and thick, while the boulder bar tails are sharp and elongated, similar to bamboo leaves. (3) The oval type: the boulder bars have smooth borders that are wide in the middle and narrow at both ends, similar to an ellipse (
Figure 2).
The remote sensing images indicate that the sickle-shaped boulder bars are mostly located in the curved parts of the river courses. The bamboo-leaf boulder bars are relatively large in length and width. These are a stable type of boulder bar. They are mostly located on the edges of the river. Elliptical boulder bars generally have small areas and are mostly located near dams and in the middle of river channels. They can be considered as the early form of boulder bars.
3.2. The Relationship between the Lengths and Widths of Boulder Bars
Figure 2 shows that before the failure of the Sedongpu natural dam, the boulder bars in the study area were within an interval of
λ ∈ [1.04, 43.74]. Moreover, 95% of the boulder bars were less than 1000 m in length and less than 100 m in width. From a data distribution perspective, the lengths and widths of the boulder bars are linearly positively correlated.
Three
SR intervals, namely (
SR < −1), (−1 ≤
SR ≤ 1) and (1 <
SR), were used to discriminate the dispersion degree of the data from the fitting curve in
Figure 2. The
λ values at the boundaries (
λb) of the three internals were 2.46 and 13.09. The boulder bars with a small
λ ∈ [1.04, 2.46], falling in interval Ⅰ, accounted for 10.5%, some of which were located within 1/3
LT from the Sedongpu natural dam, and some within 2/3
LT–
LT. The boulder bars at
λ ∈ [2.46, 13.09] falling in interval Ⅱ, which accounted for 77%, were distributed throughout the entire study area. Boulder bars with
λ(
λ ∈ [13.09, 43.74]), accounting for 12.5% and falling in interval Ⅲ, were mainly distributed within 1/8
LT–
LT downstream of the Sedongpu natural dam, and most of them were present in the curved parts of river channels.
After the Sedongpu natural dam failure, the λ values were found to be in the range of [1.58, 23.86], and 94% of the boulder bars were less than 1000 m in length and less than 100 m in width. Similar to those before the failure of the Sedongpu natural dam, the lengths and widths of boulder bars have a linear positive correlation. As noted above, the boulder bars with small λ values ([1.58, 2.66]), falling in interval Ⅰ, accounted for 7.5%, all of which were concentrated within the range of 1/4 LT–LT downstream of the Sedongpu natural dam, and most of them were boulder bars that were newly formed after the outburst flood. The boulder bars at λ ∈ [2.66, 14.92], accounting for 70% and falling in interval Ⅱ, were distributed throughout the entire study area. The boulder bars had large λ values within the range of ([14.92, 23.86]), accounting for 22.5% and falling in interval Ⅲ, and were distributed within 1/4 LT–LT of the downstream reach length of the natural dam in the curved part of the river channel.
When comparing the size characteristics of the boulder bars before and after the failure of the Sedongpu natural dam, it can be found that the λ values of the boulder bars in the whole study area changed little after the dam failure, except for the large λ values in the purple area in the figure. In addition, the fitting lines are approximately the same (slope and intercept of the line) for all the data before and after dam failure. This means that the outburst flood did not change the relationship between the lengths and widths of boulder bars.
3.3. The Relationship between the Perimeters and Lengths of the Boulder Bars
Based on the
λb values for the three intervals in
Figure 2, the relationship of
α-λ, which is the relationship between the perimeters and lengths of the boulder bars, was also divided into three intervals before and after the Sedongpu natural dam failure (
Figure 3). The figure shows that except for interval III,
α changed little before and after the dam failure for the corresponding intervals. The perimeters in interval III decreased after the dam failure. However, the percentage of boulder bars in this interval was much smaller than the sum of the other two intervals.
The linear relationship between
λ and
α was the same before and after the failure (i.e., slope of the line and the intercept at the origin), although the upper and lower limits of
α and
λ of different intervals shifted towards a bit higher value after the failure. This means that the relationship between the perimeters and lengths was not affected by the outburst flood and follow a linear relationship in the form of Equation (6):
The Equation (6) suggests that the boulder bars were “self-similar” both before and after the natural dam failure. Specifically, the result indicates that the boulder bars conserved self-similar morphology features after the sudden outburst flood. This adds to the previous observations of the lack of transient river dynamics records. To further investigate self-similarity, the relationship between the number of boulder bars whose perimeters were larger than
P and
P was investigated (
Figure 4) [
20]. This showed that there is a power–law relationship between the number of boulder bars with perimeters greater than
P and
P. This power distribution is equivalent to a fractal distribution [
20], which proves the self-similarity of the boulder bars.
3.4. Relationship between Boulder Bar Areas and Lengths
Based on the above analysis, it found that the lengths of the boulder bars determine their perimeters. To explore whether the boulder bar lengths are also a key parameter affecting the boulder bar area (A), we used the correlation coefficient to analyze the degree of correlation between the areas, lengths and widths of the boulder bars.
According to Equation (4), before and after the dam failure, the correlation coefficient PXY for areas, lengths and widths were high. This means that the boulder bar areas before and after the dam failure were more related to the lengths.
The relationship between
β and
λ is shown in
Figure 5, which reflects the relationship between the boulder bar areas and lengths. Three partitions are defined based on the values of
λb before and after the failure of the Sedongpu natural dam, as above in
Figure 2. The figure shows that except for interval III,
β changed little before and after the dam failure for the corresponding intervals. The boulder bar areas in interval III decreased after the dam failure. However, the percentage of boulder bars in this interval was much smaller than the sum of the other two intervals. There is a nonlinear relationship between
β and
λ before and after the dam failure, which can be described by a power exponential function, namely Equation (7):
This equation is applicable to
λ values with a range of 0–45. The equation provides a quantitative method for calculating the boulder bar areas. In addition, the power exponent equation between
β and
λ before and after the dam failure remained unchanged, which suggests a fractal relationship. This indicates the self-similar shapes of the boulder bars before and after the natural dam failure, which coincides with the results in
Section 3.3.
To analyze the characteristics of the distributions of the number of boulder bars in the study area, we divided the river course into ten sections with lengths of 17.3 km each (
Figure 6). Before the dam failure, the number of boulder bars was the largest at 0–17.3 km, i.e., 24 within the section. In addition, the number was the smallest, i.e., 2, at 69.2–86.5 km within the section. After the dam failure, the number of boulder bars was the largest, i.e., 17, at 86.5–103.9 km, and the number was the smallest, i.e., 2, at 121.2–138.5 km.
In addition, the number of boulder bars after the dam break was significantly smaller than that before the dam break, with a total decrease of 37. This indicates that the number of boulder bars was affected by the outburst flood. By comparison, the number of boulder bars decreased the most, by 18, within 0–17.3 km. Therefore, there is a pattern: as the distance from the dam increased, the decrease in the number of boulder bars became smaller. This suggests that the impact of outburst floods on the number of boulder bars is strongest near the dam.
3.5. Evolution Modes of Boulder Bars
The evolution of boulder bars that are affected by outburst floods can generally be categorized into four modes, as shown in
Figure 7. Mode 1 shows that the original boulder bars disappeared after the outburst floods. This mode mostly occurs in the straight segments of the channel. When the outburst flood flows through a straight channel portion, some sediments in the boulder bars are carried away, causing the disappearance of boulder bars.
Mode 2 suggests that the boulder bars were displaced by the outburst flood. In this mode, the shapes and sizes of the boulder bars were similar before and after the dam failure, and the variations in area were less than 10% of the original boulder bars. Mode 3 means that the boulder bars were synthesized from multiple boulder bars into one by the outburst flood. This mode mainly occurred on the convex banks of river bends. Mode 4 indicates that the sizes of the boulder bars decreased due to outburst floods. However, in this mode, the reduction of area of the boulder bars affected by floods was small, i.e., generally less than 10% of the original boulder bars.