The Yellow River (YR) is famous for its ‘hanging river’ due to the high sedimentation rate in the lower reaches. Reservoir sediment retention (RSR) and water–sediment regulation (WSR) of the Xiaolangdi Reservoir (XLDR) contribute to mitigating the siltation and adjusting the bankfull channel geometry in the lower Yellow River (LYR). However, few studies have been conducted to quantify the contribution of RSR and WSR to medium-sized channel shaping. In this paper, we analyzed the temporal and spatial variations of bankfull discharge (BFD) in the LYR, established empirical expressions between the BFD and the previous five-year average discharge above 2,000 m3/s and the incoming sediment coefficient during flood seasons, and then constructed three different water and sediment scenarios to quantify the contribution of the WSR of the XLDR to medium-sized channel shaping in the LYR. Calculated results indicated that WSR accounted alone for 68%–94% of the medium-sized channel shaping, with RSR contributing 6%–32%. The contribution of WSR to the recovery of medium-sized channels increased with time and became stabilized eventually. Along the LYR, such contribution also increased from 68% to 97% and eventually became stabilized in narrow reaches, reaching over 90%. The results will lay a solid technical foundation for the WSR of the XLDR, river regulation downstream of the reservoir, and sediment disposal of the YR.

  • The contributions of RSR and WSR to the medium-sized channel shaping of the LYR were 88% and 12% respectively in the period 2002–2015.

  • Contribution rate of WSR to medium-sized channel shaping increased first and then became stabilized with time.

  • The influence of WSR on medium-sized channel shaping increased and eventually became stabilized in narrow reaches along the LYR.

The Yellow River (YR) is well known in the world due to its high sediment concentration and rapid sedimentation rate in its lower reach (Wang & Li 2011). The riverbed in the lower Yellow River (LYR) has been dramatically aggrading for many years, and is now up to 10 m higher than the overbank plains in some places, forming a so-called ‘suspended river’ or ‘perched river’ due to its high sedimentation rate (Wu et al. 2008a). The main-channel width in the LYR varied from 0.9 to 1.4 km after the 1999 flood season, and there are extensive floodplains utilized by local inhabitants on both sides of the LYR, the area of which accounts for over 80% of the whole river (Xia et al. 2014a; Cheng et al. 2021). Therefore, studying the adjustment characteristics and influencing factors of bankfull discharge (BFD) (herein referring to the flow at which the in-channel water level is about to overtop its banks onto the active floodplain) would provide significant insights into additional information to guide further river management and flood prevention in the LYR.

Through long-term practical exploration of YR management, reservoir sediment retention (RSR) and water and sediment regulation (WSR) have been proposed for disposal and utilization of the sediment. Before the operation of the Xiaolangdi Reservoir (XLDR), the LYR was characterized by obvious shrinkage of the main channel accompanied by a sharp decrease in the flood discharge capacity (Bi et al. 2019), and later the situation was improved due to the reservoir operation. Upon its completion in 1999, the XLDR has been used to retain most of the inflowing sediment to the LYR. To further desilt the lower reaches and recover the medium-sized channels, a water–sediment regulation (WSR) project was initiated with the XLDR as the core in 2002. Remarkable effects of the WSR on the LYR have been reported (Wang et al. 2017), including changes in the river channels (Xu & Si 2009; Chen et al. 2012; Ma et al. 2012), hydrological characteristics (Lu et al. 2022), and sedimentation features (Kong et al. 2022), as well as changes in the evolution of the Yellow River Delta (Kong et al. 2015; Wang et al. 2017).

WSR was accompanied with a large-flow discharge process that was conducive to the sediment transport of the lower reaches and the recovery and maintenance of the medium-sized channels (Zhang et al. 2009; Xia et al. 2014a). The BFD at Huayuankou Station increased from 3,700 m3/s in 2000 to 6,900 m3/s in 2012 (Qi et al. 2013) and, overall, the LYR channel has switched from a sediment-deposition state to an erosive state (Xu et al. 2005). Some studies have been conducted on the operation of the XLDR and the variation of BFD of the medium-sized channels in the LYR. Wu et al. (2008a) pointed out that the variation of BFD could be predicted by using the changes in the incoming discharge and sediment load in the LYR. Li et al. (2015), Xia et al. (2010a, 2010b) and Chen et al. (2018) studied the BFD in the LYR, and analyzed the variation of BFD and its relation with water and sediment. With the operation of the XLDR, the flow and sediment regime entering the LYR has been altered dramatically, and the LYR channel is experiencing continuous degradation, leading to a significant variation in the bankfull channel geometry (Xia et al. 2010a, 2010b, 2014a, 2014b). Hu et al. (2012) analyzed the effect of various water and sediment combinations released from the XLDR on the shaping of the main channel and found that under recent incoming flow and sediment conditions, it is possible to shape and maintain a medium-sized channel with a BFD of approximately 4,000 m3/s. Miao et al. (2016) analyzed the impact of WSR of the XLDR on the scouring and silting of the LYR. Bi et al. (2019) analyzed the scouring and silting characteristics of the LYR before and after the operation of the XLDR. Chen et al. (2021) put forward the operation mode of the XLDR based on the maintenance of the medium-sized channels in the LYR. Cheng et al. (2021) estimated the adjustment characteristics of various bankfull parameters over the period 1986–2015 in a braided reach of the LYR. All of these studies have gone into the operation of the XLDR, the scouring and silting of the middle and lower reaches, and the evolution of the medium-sized channels. However, no research has been done into the quantitative contributions of RSR and WSR to medium-sized channel shaping. Whether WSR or RSR plays a greater role in shaping the medium-sized channels in the LYR remains to be explored. In the existing literature, the scale of the restored medium-sized channels is often attributed to the role of WSR, ignoring the contribution of RSR, resulting in misunderstandings, which are necessary to clarify.

In this paper, based on the analysis of the variation of BFD of the medium-sized channels in the LYR, we established the response relationship between the BFD and the inflow sediment coefficient in the flood season and the large water discharge within the year. On this basis, we further constructed different water and sediment scenarios to analyze and calculate the quantitative contribution of WSR of the reservoir to the medium-sized channel shaping in the LYR. In addition, spatial and temporal analyses of the contribution of current WSR on the LYR were conducted. The results will lay a solid technical foundation for WSR of the XLDR, the river regulation downstream of the reservoir, and the sediment disposal of the YR as well as other large rivers with severe sedimentation throughout the world.

Study area

The YR, originating from the Qinghai–Tibetan Plateau and extending toward the Bohai Sea, is the second longest river in China (5,464 km) and the largest river in the world in terms of sediment concentration (Zhang et al. 2021a, 2021b). The mean annual suspended load measured at Sanmenxia Station is 1.6 billion tonnes. The average sediment concentration is 35 kg/m3, and hyperconcentrated floods with sediment concentrations over 100 kg/m3 occur quite often (Wu et al. 2008a, 2008b; Xia et al. 2010a, 2010b). The river basin is mostly arid and semiarid with an annual average temperature of 8–14 °C and an annual average precipitation of about 442 mm (Chen et al. 2005; Cheng et al. 2021). For the purpose of this paper, the research covers the reaches downstream of Tongguan, mainly including XLDR, Sanmenxia Reservoir and the LYR (Figure 1).
Figure 1

Map of the middle and lower reaches of the YR and locations of the main hydrologic stations.

Figure 1

Map of the middle and lower reaches of the YR and locations of the main hydrologic stations.

Close modal

The Xiaolangdi Hydroproject is located at the exit of the last canyon in the middle reaches of the Yellow River, 130 km away from Sanmenxia Hydroproject on the upper reaches and 131.9 km away from Huayuankou Hydrologic Station on the lower reaches. As a backbone project to control the WSR system of the YR (Zhang et al. 2021a, 2021b), the Xiaolangdi Hydroproject has been developed primarily for flood control (including ice prevention) and desilting, and secondarily for water supply, irrigation and power generation, with the aim to create favorable conditions to fend off natural disasters and boost comprehensive utilization of resources. The reservoir has a designed normal storage level of 275 m and a total storage of 12.65 billion m³, including 7.55 billion m³ for sediment retention and 1 billion m³ for WSR. The gate was closed for impoundment in October 1999, and the reservoir is currently in the first phase of the later sediment-retention stage (Chen et al. 2016).

The Sanmenxia Reservoir is the first large-scale hydroproject built on the main stream of the middle Yellow River, developed primarily for flood control, and secondarily for ice prevention, irrigation, power generation and water supply. Since the gating for impoundment in September 1960, the reservoir has been reconstructed twice, and experienced three operation stages of ‘storing water and retaining sediment’, ‘detaining flood and discharging sediment ‘ and ‘storing clean water and discharging muddy flow’. To maintain active storage and control sediment deposition from affecting the elevation of river reaches at Tongguan (Yang et al. 2020), the reservoir is limited to operate at a low water level of 305 m in the flood season, with the corresponding storage capacity of only about 50 million m3, rendering its water and sediment regulating capacity weak. The reservoir area is basically under an equilibrium between scouring and deposition.

The LYR, with a length of around 740 km, is defined as the reach between Mengjin in Henan province and Lijin in Shandong province. The bed elevation in the LYR drops significantly by about 94 m, with the gradient ranging from 0.265‰ to 0.1‰. Here, the level of the riverbed rising above cities and farmland along most of the LYR is as a result of the heavy sedimentation in the river channel, a phenomenon known as a ‘suspended river’ or ‘perched river’ (Figure 2). Notably wide floodplains and water channels can be found in the lower reaches. Different from other river channels, the floodplain in the LYR is densely populated. The LYR is usually divided into three distinct reaches based on their geomorphological characteristics (Wu et al. 2008a). (i) The reach from Mengjin to Gaocun has a typical braided channel pattern with a sinuosity as low as 1.15. It extends a length of 284 km and has a large width of 5–20 km between the left and right levees, a main channel width of 0.5–3.0 km and a longitudinal channel slope of 0.1‰–0.5‰. The channel depth is only about 2 m for flows under bankfull conditions. (ii) The reach roughly from Aishan to Lijin, with a length of 301 km and an average sinuosity of 1.21, is a well-restricted and stable meandering reach. The average channel slope is 0.09‰, and the distance between levees is 0.5 and 5.0 km. The main channel has a narrow and deep section as compared with the braided reach, with a width of 0.3 km to 0.7 km and a depth of 2–3 m. (iii) The reach between Gaocun and Taochengpu has a transitional channel pattern from braided to meandering with a length of 155 km. The reach is characterized by a width between two levees of 1.4–8.5 km, a main-channel width of 0.4–1.2 km, a main-channel depth of 2–3 m, and a longitudinal channel slope of 0.11‰ (Wu et al. 2008a; Xia et al. 2014b).
Figure 2

A generalized diagram of a typical cross-section of the LYR.

Figure 2

A generalized diagram of a typical cross-section of the LYR.

Close modal

Data source

The water and sediment into the Sanmenxia Reservoir and the XLDR mainly come from the main stream, as there are no large tributaries in these reservoir areas. The Tongguan Hydrologic Station is located at the inlet of the Sanmenxia Reservoir. The Sanmenxia Hydrologic Station is located at the exit of the Sanmenxia Reservoir and the inlet of the XLDR. The Xiaolangdi Hydrologic Station is located at the outlet of the XLDR, and serves as a water and sediment control station for the main stream joining the LYR. In addition, there are two major tributaries merging into the LYR, i.e. Qinhe River and Yiluo River, with the Wuzhi and Heishiguan Hydrologic Stations for inflowing water and sediment control. There are many hydrologic stations on the main stream of the LYR (see Figure 1), including Huayuankou, Gaocun, Sunkou, Aishan, and Lijin, which are 132 km, 310 km, 431 km and 766 km away from the Xiaolangdi Dam, respectively. These stations are mainly tasked to observe such data as discharge, sediment discharge, sediment concentration, water level, cross-section and topography, which are then compiled and issued by the hydrological authorities. The research in this paper is based on these measured data. The BFD in the cross-section of each hydrologic station on the LYR is obtained from the measured stage–discharge relationship with reference to the position of the cross-section at the hydrologic station (Xia et al. 2010b), as shown in Figure 3.
Figure 3

Measured bankfull discharge in 2005 using the mean daily stage–discharge relation at Gaocun: (a) determination of bankfull stage; (b) obtaining of bankfull discharge using the stage–discharge curve.

Figure 3

Measured bankfull discharge in 2005 using the mean daily stage–discharge relation at Gaocun: (a) determination of bankfull stage; (b) obtaining of bankfull discharge using the stage–discharge curve.

Close modal

Methods

We first established expressions to describe the response relationship between the BFD and the water and sediment using the measured data. On this basis, we then calculated the variation of BFD according to the water and sediment process in three different scenarios. Finally, we obtained the quantitative contribution of WSR to the medium-sized channel shaping by comparing the variation of BFD between these three different scenarios.

Calculation and verification of BFD

A medium-sized channel is the main channel for flood and sediment transport, and its discharge capacity can be expressed by BFD (Hu et al. 2006). The variation of BFD is not only associated with the conditions of inflowing water and sediment to the LYR in the corresponding year, but also with the conditions in the early stage. It is generally the synthetic result of the water and sediment conditions over the years. To analyze the response relationship between BFD and water and sediment, the following have been considered: (1) for an alluvial channel, the water and sediment conditions in the flood season generally have an obvious influence on the shaping of medium-sized channels. Therefore, factors including water discharge and inflow sediment coefficient in the flood season can be taken as the influencing factors of BFD. To reflect the influence of water and sediment conditions in the non-flood season on BFD, such factors as annual average water discharge and annual average inflow sediment coefficient can also be selected as the influencing factors of BFD, with due consideration given to the delayed response of BFD to water and sediment (Wu et al. 2008a). (2) The shaping of medium-sized channels is not only related to the water and sediment inflow, but also to the discharge process. The larger the discharge is, the greater the role of sediment transport in medium-sized channel shaping will be. Under the same water inflow, the greater the number of days is with a larger discharge, and the more favorable it will be to shape and maintain the medium-sized channels. Based on previous research results, the water discharge above a given value is selected as the influencing factor of BFD. Given all this, in this paper, we analyzed and calculated the correlation coefficient R2 between the BFD at each station and its influencing factors using the data during the period from 1974 to 2015 (42 years in total). The calculation results indicate that the BFD at each station is non-linearly related to the two factors, i.e. water discharge above 2,000 m3/s in the five years before discharge to the LYR and the inflow sediment coefficient in the flood season, with the strongest correlation (Figure 4 shows the variation of the correlation coefficient between the minimum BFD and these influencing factors). Therefore, these two factors are selected as the influencing factors of BFD.
Figure 4

Variation of correlation coefficient between minimum BFD and influencing factors.

Figure 4

Variation of correlation coefficient between minimum BFD and influencing factors.

Close modal
The non-linear relationship is expressed as follows:
(1)
where,
  • – BFD (m3/s);

  • , and – coefficients to be determined;

  • – water discharge above 2,000 m3/s in the LYR in the first five years (100 million m3);

  • – inflow sediment coefficient of the LYR in the flood season.

Through regression analysis of the measured data about the LYR, the response relationships between the BFD at Huayuankou, Gaocun, Sunkou, and Lijin and the minimum BFD and the water and sediment are expressed through calibration as follows:
(2)
(3)
(4)
(5)
(6)
where,
  • – BFD at Huayuankou (m3/s);

  • – BFD at Gaocun (m3/s);

  • – BFD at Sunkou (m3/s);

  • – BFD at Lijin (m3/s);

  • – minimum BFD of the LYR (m3/s).

In Equations (2)–(6), the comprehensive correlation coefficient is above 0.8. The BFD calculated as per these expressions is close to the measured value. Figure 5 is the correlation scatter diagram of the calculated and measured BFDs at Gaocun Hydrologic Station. It can be seen from the figure that the scatters are centralized near the 45° line, indicating that the calibrated expressions of response relationships between the BFD and the water and sediment are credible.
Figure 5

Comparison between calculated and measured BFDs at Gaocun Hydrologic Station.

Figure 5

Comparison between calculated and measured BFDs at Gaocun Hydrologic Station.

Close modal

Scenario setting

Upon its completion, the XLDR has successfully reduced the inflowing sediment to the LYR, and played a significant role in scouring and silting of river channels and shaping of medium-sized channels. In order to separately analyze the effect of WSR on shaping the medium-sized channels in the LYR, it is necessary to first analyze the role of WSR in controlling the inflowing water and sediment to the LYR. To this end, three scenarios in this regard are developed according to the locations of reservoirs and hydrologic stations and the WSR of reservoirs:

  • Scenario 1: The water and sediment processes at Tongguan, Heishiguan and Wuzhi Hydrologic Stations represent the water and sediment process without RSR or WSR (i.e. without reservoir regulation).

  • Scenario 2: The water and sediment processes at Xiaolangdi, Heishiguan and Wuzhi Hydrologic Stations represent the water and sediment process with RSR and WSR.

  • Scenario 3: The water process at Tongguan Hydrologic Station, the sediment process at Xiaolangdi Hydrologic Station, and the water and sediment processes at Heishiguan and Wuzhi Hydrologic Stations represent the water and sediment process with only RSR.

The annual average water discharges under scenarios 1, 2 and 3 are 26.175 billion m³, 27.176 billion m³ and 26.175 billion m³, respectively (as shown in Table 1), in which the water discharge in the flood season accounts for 48.5%, 36.9% and 48.5%, respectively, and the annual average water discharge above 2,000 m³/s accounts for 16.7%, 23.0% and 16.7%, respectively. The annual sediment loads are 237 million tonnes, 71 million tonnes and 71 million tonnes, respectively, in which the sediment loads in the flood season account for 81.1%, 94.8% and 94.8%, respectively. It is thus clear that the WSR of the reservoir has increased the large-flow discharge to the LYR. The annual average number of days with discharge above 2,000 m³/s has increased from 16.93 days to 24.64 days, and the corresponding water discharge from 4.379 billion m³ to 6.262 billion m³. Through reservoir sediment retention, the annual average sediment inflow to the LYR has dropped from 237 million tonnes to 71 million tonnes.

Table 1

Inflowing water and sediment to the LYR under different scenarios

Scenario no.Annual average water discharge (100 million m³)Annual average sediment load (100 million tonnes)Water discharge in flood season (100 million m³)Sediment load in flood season (100 million tonnes)Number of days with discharge over 2,000 m³/s (days)Annual average water discharge over 2,000 m³/s (100 million m³)
261.75 2.37 126.94 1.92 16.93 43.79 
271.76 0.71 100.30 0.67 24.64 62.62 
261.75 0.71 126.94 0.67 16.93 43.79 
Scenario no.Annual average water discharge (100 million m³)Annual average sediment load (100 million tonnes)Water discharge in flood season (100 million m³)Sediment load in flood season (100 million tonnes)Number of days with discharge over 2,000 m³/s (days)Annual average water discharge over 2,000 m³/s (100 million m³)
261.75 2.37 126.94 1.92 16.93 43.79 
271.76 0.71 100.30 0.67 24.64 62.62 
261.75 0.71 126.94 0.67 16.93 43.79 

BFD in the LYR under different scenarios

Based on the process of inflowing water and sediment to the LYR in different scenarios, the variation of BFD at each station in these scenarios is calculated as per the established expressions of relationships between BFD and its influencing factors (i.e. water and sediment), and the results are shown in Table 2 and Figure 6.
Table 2

Calculated BFD at each station under different scenarios

ItemHuayuankou
Gaocun
Sunkou
Lijin
Minimum BFD
Scenario 1Scenario 2Scenario 3Scenario 1Scenario 2Scenario 3Scenario 1Scenario 2Scenario 3Scenario 1Scenario 2Scenario 3Scenario 1Scenario 2Scenario 3
2002 3,075 3,890 3,717 1,717 2,244 1,984 1,679 1,985 1,747 2,540 2,902 2,690 1,630 1,953 1,738 
2003 3,967 4,785 4,687 2,763 3,317 3,136 2,583 2,835 2,675 3,345 3,641 3,518 2,454 2,738 2,596 
2004 3,875 5,544 5,078 2,561 3,988 3,145 2,378 3,221 2,516 3,188 4,006 3,459 2,278 3,127 2,496 
2005 3,778 5,682 4,811 2,332 4,395 2,802 2,140 3,596 2,251 3,002 4,264 3,230 2,074 3,449 2,250 
2006 3,856 5,965 4,751 2,208 4,790 2,588 1,957 3,880 2,045 2,876 4,479 3,063 1,927 3,709 2,068 
2007 3,921 6,167 4,442 1,886 5,038 2,074 1,564 4,041 1,605 2,564 4,605 2,663 1,594 3,861 1,663 
2008 3,889 5,949 4,385 1,825 4,564 1,999 1,506 3,644 1,544 2,510 4,332 2,603 1,542 3,514 1,606 
2009 4,560 6,683 5,485 2,508 4,924 2,885 2,027 3,673 2,106 3,025 4,448 3,198 2,039 3,603 2,170 
2010 4,551 6,698 5,345 2,411 5,017 2,725 1,926 3,759 1,992 2,944 4,505 3,090 1,951 3,677 2,060 
2011 4,832 6,805 5,524 2,609 5,076 2,888 2,039 3,770 2,097 3,069 4,526 3,195 2,069 3,695 2,164 
2012 4,960 6,901 5,471 2,658 5,365 2,864 2,047 4,014 2,089 3,090 4,690 3,183 2,084 3,908 2,154 
2013 4,982 7,047 5,426 2,646 5,727 2,823 2,027 4,307 2,064 3,077 4,886 3,158 2,068 4,165 2,129 
2014 3,980 7,247 4,040 1,351 6,027 1,367 988 4,407 991 2,021 5,036 2,030 1,079 4,265 1,085 
2015 4,152 7,247 4,382 1,470 6,027 1,531 1,069 4,407 1,081 2,123 5,036 2,158 1,162 4,265 1,183 
ItemHuayuankou
Gaocun
Sunkou
Lijin
Minimum BFD
Scenario 1Scenario 2Scenario 3Scenario 1Scenario 2Scenario 3Scenario 1Scenario 2Scenario 3Scenario 1Scenario 2Scenario 3Scenario 1Scenario 2Scenario 3
2002 3,075 3,890 3,717 1,717 2,244 1,984 1,679 1,985 1,747 2,540 2,902 2,690 1,630 1,953 1,738 
2003 3,967 4,785 4,687 2,763 3,317 3,136 2,583 2,835 2,675 3,345 3,641 3,518 2,454 2,738 2,596 
2004 3,875 5,544 5,078 2,561 3,988 3,145 2,378 3,221 2,516 3,188 4,006 3,459 2,278 3,127 2,496 
2005 3,778 5,682 4,811 2,332 4,395 2,802 2,140 3,596 2,251 3,002 4,264 3,230 2,074 3,449 2,250 
2006 3,856 5,965 4,751 2,208 4,790 2,588 1,957 3,880 2,045 2,876 4,479 3,063 1,927 3,709 2,068 
2007 3,921 6,167 4,442 1,886 5,038 2,074 1,564 4,041 1,605 2,564 4,605 2,663 1,594 3,861 1,663 
2008 3,889 5,949 4,385 1,825 4,564 1,999 1,506 3,644 1,544 2,510 4,332 2,603 1,542 3,514 1,606 
2009 4,560 6,683 5,485 2,508 4,924 2,885 2,027 3,673 2,106 3,025 4,448 3,198 2,039 3,603 2,170 
2010 4,551 6,698 5,345 2,411 5,017 2,725 1,926 3,759 1,992 2,944 4,505 3,090 1,951 3,677 2,060 
2011 4,832 6,805 5,524 2,609 5,076 2,888 2,039 3,770 2,097 3,069 4,526 3,195 2,069 3,695 2,164 
2012 4,960 6,901 5,471 2,658 5,365 2,864 2,047 4,014 2,089 3,090 4,690 3,183 2,084 3,908 2,154 
2013 4,982 7,047 5,426 2,646 5,727 2,823 2,027 4,307 2,064 3,077 4,886 3,158 2,068 4,165 2,129 
2014 3,980 7,247 4,040 1,351 6,027 1,367 988 4,407 991 2,021 5,036 2,030 1,079 4,265 1,085 
2015 4,152 7,247 4,382 1,470 6,027 1,531 1,069 4,407 1,081 2,123 5,036 2,158 1,162 4,265 1,183 
Figure 6

BFD at each station under different scenarios: (a) Huayuankou, (b) Gaocun, (c) Sunkou, (d) Lijin.

Figure 6

BFD at each station under different scenarios: (a) Huayuankou, (b) Gaocun, (c) Sunkou, (d) Lijin.

Close modal

It can be seen from the table that the BFD in the three scenarios shows different variations with time. In scenarios 1 and 3, the calculated BFD fluctuates, showing a trend of increasing first and then decreasing, except for the Huayuankou Hydrologic Station, where the BFD at the beginning of the period is greater than that at the end of the period. In scenario 2, the calculated BFD at all stations shows an overall increasing trend, with the BFD at the beginning less than that at the end. The variation of calculated BFD at the beginning and end of the period in the three scenarios is compared as follows.

In scenario 1, the calculated BFD during the period from 2002 to 2015 increased from 3,075 m³/s to 4,152 m³/s at Huayuankou Station, but decreased from 1,717 m³/s to 1,470 m³/s at Gaocun Station, from 1,679 m³/s to 1,069 m³/s at Sunkou Station, and from 2,540 m³/s to 2,123 m3/s at Lijin Station. The minimum BFD during this period also decreased from 1,630 m³/s to 1,162 m³/s. It is thus clear that in this scenario, the BFD decreased at all other stations except at Huayuankou Station, so this scenario represents the variation trend of BFD without WSR of the reservoir.

In scenario 2, the calculated BFD during the period from 2002 to 2015 increased from 3,890 m³/s to 7,247 m³/s at Huayuankou Station, from 2,244 m³/s to 6,027 m³/s at Gaocun Station, from 1,985 m³/s to 4,407 m³/s at Sunkou Station, and from 2,902 m³/s to 5,036 m³/s at Lijin Station. The minimum BFD during this period also increased from 1,953 m³/s to 4,265 m³/s. It is thus clear that in this scenario, the BFD at all stations increased, and the variation of calculated BFD was consistent with the measured one in the lower Yellow River, so this scenario represents the variation of BFD under the combined effect of RSR and WSR under natural water and sediment conditions.

In scenario 3, the calculated BFD during the period from 2002 to 2015 increased from 3,717 m³/s to 4,382 m³/s at Huayuankou Station, but decreased from 1,984 m³/s to 1,531 m³/s at Gaocun Station, from 1,747 m³/s to 1,081 m³/s at Sunkou Station, and from 2,690 m³/s to 2,518 m³/s at Lijin Station. The minimum BFD during this period also decreased from 1,738 m³/s to 1,183 m³/s. It is thus clear that in this scenario, the BFD decreased at all other stations except at Huayuankou Station close to the dam, indicating that it is impossible to recover the medium-sized channels in a large scale from a long distance only by RSR.

Longitudinal variations in contribution rate of WSR to medium-sized shaping

Through comparison of the calculation results in these three scenarios, we can obtain the quantitative contributions made by RSR and WSR to the recovery of medium-sized channels, To be specific, the difference in the calculation results between scenarios 1 and 2 shows that the quantitative contributions made by RSR and WSR to the BFD at Huayuankou, Gaocun, Sunkou and Lijin, and to the minimum BFD are 2,016 m³/s, 2,540 m³/s, 1,829 m³/s, 1,571 m³/s and 1,712 m³/s, respectively (as shown in Table 3). The difference in the calculation results between scenarios 2 and 3 shows that the quantitative contributions of WSR to the BFD at Huayuankou, Gaocun, Sunkou and Lijin, and to the minimum BFD are 1,361 m³/s, 2,263 m³/s, 1,766 m³/s, 1,437 m³/s and 1,612 m³/s, respectively. The difference in the calculation results between scenarios 1 and 3 shows that the quantitative contributions of RSR to the BFD at Huayuankou, Gaocun, Sunkou and Lijin, and to the minimum BFD are 655 m³/s, 277 m³/s, 63 m³/s, 134 m3/s and 100 m³/s, respectively. It is thus clear that RSR and WSR have increased the average BFD at each station. To be specific, the average BFD at Huayuankou is increased by 2,016 m³/s, including 1,361 m³/s (68%) contributed by WSR and 655 m³/s (32%) by RSR. That at Gaocun is increased by 2,540 m³/s, including 2,263 m³/s (89%) contributed by WSR and 277 m³/s (11%) by RSR. That at Sunkou is increased by 1,829 m³/s, including 1,766 m³/s (97%) contributed by WSR and 63 m³/s (3%) by RSR. That at Lijin is increased by 1,571 m³/s, including 1,437 m³/s (92%) contributed by WSR and 134 m³/s (8%) by RSR. The average minimum BFD is increased by 1,712 m³/s, including 1,612 m³/s (94%) contributed by WSR and 100 m³/s (6%) by RSR.

Table 3

Contribution of reservoir operation to medium-sized channel shaping at each station in different scenarios

ItemAverage (m3/s)
Comparative difference (m3/s)
Contribution of reservoir operation to medium-sized channel shaping
Scenario 1 ①Scenario 2 ②Scenario 3 ③②-①②-③③-①WSRRSR
Huayuankou 4,170 6,186 4,825 2,016 1,361 655 68% 32% 
Gaocun 2,210 4,750 2,487 2,540 2,263 277 89% 11% 
Sunkou 1,852 3,681 1,915 1,829 1,766 63 97% 3% 
Lijin 2,812 4,383 2,946 1,571 1,437 134 92% 8% 
Minimum BFD 1,854 3,566 1,954 1,712 1,612 100 94% 6% 
ItemAverage (m3/s)
Comparative difference (m3/s)
Contribution of reservoir operation to medium-sized channel shaping
Scenario 1 ①Scenario 2 ②Scenario 3 ③②-①②-③③-①WSRRSR
Huayuankou 4,170 6,186 4,825 2,016 1,361 655 68% 32% 
Gaocun 2,210 4,750 2,487 2,540 2,263 277 89% 11% 
Sunkou 1,852 3,681 1,915 1,829 1,766 63 97% 3% 
Lijin 2,812 4,383 2,946 1,571 1,437 134 92% 8% 
Minimum BFD 1,854 3,566 1,954 1,712 1,612 100 94% 6% 

The results revealed that the WSR is more efficient for shaping and maintaining a medium-sized channel. After the implementation of the WSR, the downstream main channel obtained comprehensive erosion, which reduced the deposition of the LYR, increased the BFD of the main channel, and improved the flow capacity (Xu & Si 2009; Hu et al. 2012; Miao et al. 2016). However, it is noteworthy that the scale of the restored medium-sized channels is not only attributed to the role of WSR; RSR also contributed to approximately 6%–32% of the medium-sized channel shaping in the LYR.

Figure 7 shows the variation of the contribution of reservoir operation to the recovery of the medium-sized channels in the LYR. The contribution of WSR increases with the increase of distance from the dam, and tends to become stable in the narrow reaches downstream of Sunkou, reaching over 90%, while the contribution of RSR decreases with the increase of distance from the dam, reaching below 10%. This law of contribution conforms to the actual conditions of water and sediment movement and scouring and silting in the river channel. The reservoir retains sediment and discharges clear water, so the sediment load in the LYR can be controlled. The closer the river channel is to the reservoir, the greater the effect of scouring and the more significant the effect of recovering the medium-sized channels will be. From 1999 to 2015, 1.918 billion tonnes of sediment were scoured in the lower reaches, including 1.404 billion tonnes in the 310-km-long reaches upstream of Gaocun, accounting for 73.2%, with the BFD in these reaches recovered to more than 6,000 m3/s. WSR is aimed at increasing the discharge flow, coordinating the equilibrium between water and sediment, improving the sediment transport capacity in water flow, and realizing the scouring in the lower reaches, especially the reaches downstream of Gaocun. During the 19 water–sediment regulations, the total amount of sediment scoured in the LYR is 430 million tonnes, including 162 million tonnes in the reaches from Gaocun to Aishan and 111 million tonnes from Aishan to Lijin, accounting for 41.0% and 30.0% of the total amount of scouring in the corresponding reaches since the operation of the reservoir, with the scouring efficiency during WSR being 7.1 times and 1.9 times that in other periods. Miao et al. (2016) noted that the cumulative channel erosion increased gradually and was then stable from the braided reach, through the transitional reach to the meandering reach. In the braided reach, the average annual sediment erosion was about 0.20 × 106 t per kilometre, whereas the average annual sediment erosion dropped to 0.02 × 106 t per kilometre in the meandering reach. Bi et al. (2019) also proved that the Gaocun–Aishan transitional channel reach witnessed the maximum depositional and erosional flux and intensity, whereas the deposition–erosion flux and intensity in the meandering channel reach were lower. Therefore, the contribution of WSR to medium-sized shaping increased and eventually became stabilized along the LYR.
Figure 7

Variation of contribution of reservoir operation to the recovery of medium-sized channel in the LYR.

Figure 7

Variation of contribution of reservoir operation to the recovery of medium-sized channel in the LYR.

Close modal
Figure 8

Variation of contribution of WSR in recovering medium-sized channel in the LYR.

Figure 8

Variation of contribution of WSR in recovering medium-sized channel in the LYR.

Close modal

Inter-annual variability in contribution rate of WSR to medium-sized channel shaping

In addition, in 2002, the first year of WSR, the contributions of WSR to the BFD at Huayuankou, Gaocun, Sunkou and Lijin and to the minimum BFD were 21%, 49%, 78%, 59% and 67%, respectively (as shown in Figure 8). In subsequent years, the contribution showed an upward trend. The contribution to the BFD of the reaches downstream of Huayuankou was gradually stabilized after 2006, keeping at over 90% in most years. The Huayuankou Station was closer to the dam, so it was greatly affected by the sediment retention and clear water discharge of the reservoir. The riverbed was easily scoured in the early stage, and the discharge of clear water during the sediment-retention period would definitely be conducive to recovering the medium-sized channel. Nevertheless, with the gradual coarsening of riverbed sediment, the effect of scouring by clear water, especially at a small discharge, would inevitably decrease (Shen et al. 2020). Miao et al. (2016) have suggested that the current WSR is no longer as effective as it was initially for the LYR, exhibiting a tendency towards functional degradation. Hence, the contribution of WSR to the BFD at Huayuankou generally showed a slow increase, accounting for more than 90% in 2014.

Through analysis of the temporal and spatial variations of BFD in the LYR and its influencing factors (i.e. water and sediment), it is clear that there is a good response relationship between the BFD at Huayuankou, Gaocun, Sunkou, Lijin and the minimum flow section and the inflow sediment coefficient in the flood season and the water discharge above 2,000 m³/s during a year. Considering the delayed response of BFD to water and sediment, expressions of response relationships are established, to fully describe the impact of the equilibrium between water and sediment and the large-flow discharge process on the variation of BFD in the lower Yellow River.

Three water and sediment scenarios are constructed based on the WSR of the XLDR. On this basis, the quantitative contributions of WSR of the reservoir to the medium-sized channel shaping in the LYR are calculated as per the expressions of response relationships. From 2002 to 2015, the contributions of the reservoir operation to the discharge capacity recovery of the medium-sized channels at Huayuankou, Gaocun, Sunkou, Lijin and the minimum flow section in the LYR were 2,016 m³/s, 2,540 m³/s, 1,829 m³/s, 1,571 m³/s and 1,712 m³/s, respectively. The contributions of WSR to the recovery of the medium-sized channels gradually increased with time and became stabilized eventually, i.e., for the Sunkou section, which increased from 63.5% in 2003 to 98.3% in 2007, and remained at around 97.8% during 2008–2015. With the increase of the distance from the dam, such contribution also increased and eventually became stabilized in narrow reaches, reaching over 90%. Along the LYR, such contribution also increased from 68% to 97% and eventually became stabilized in narrow reaches, reaching over 90%. The research results indicate that the large discharge process during WSR of the XLDR plays a significant role in shaping the medium-sized channels in the LYR.

With the changes in the natural environment of the YR Basin and the increasingly stringent requirements on basin protection and management, and given the fact that the YR has a small amount of water resources and there is a tension between supply and demand, it is recommended to focus research on water discharge for WSR in the future, in order to optimize the discharge process through reservoir operation, thus improving the efficiency of sediment transport and medium-sized channel shaping during WSR.

The current investigation was financially supported by Young Elite Scientists Sponsorship Program by Henan Association for Science and Technology (grant number 2022HYTP022), the National Key R&D Program of China (grant number 2016YFC0402503 and 2018YFC1508404) and Science-Technology Program for Excellent Talents in Yellow River Conservancy Commission of the Ministry of Water Resources (HQK-202319).

Data cannot be made publicly available; readers should contact the corresponding author for details.

We declare that we are not affiliated with or involved with any organization or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this paper.

Chen
J.
,
Wang
F.
,
Meybeck
M.
,
He
D.
,
Xia
X.
&
Zhang
L.
2005
Spatial and temporal analysis of water chemistry records (1958–2000) in the Huanghe (Yellow River) basin
.
Global Biogeochemical Cycles
19
,
GB3016
.
http://dx.doi.org/10.1029/2004GB002325
.
Chen
J.
,
Zhou
W.
&
Sun
G.
2016
Studies on operation scheme and water–sediment regulation of Xiaolangdi reservoir
.
Journal of Sediment Research
2016 (4),
1
8
(in Chinese). https://doi.org/10.16239/j.cnki.0468-155x.2016.04.001
.
Chen
L.
,
Hu
C.
&
Chen
X.
2018
Relationship between bank-full discharges and processes of flow-sediment in lower Yellow River
.
Journal of Sediment Research
43
(
4
),
1
7
(in Chinese). https://doi.org/10.16239/j.cnki.0468-155x.2018.04.001
.
Chen
C.
,
Lu
J.
,
Wu
M.
,
Wei
S.
&
Cui
Z.
2021
Study on the operation of Xiaolangdi Reservoir based on maintaining the river channel of median water of the lower Yellow River
.
Yellow River
43
(
11
),
35
39
(in Chinese)
.
Cheng
Y.
,
Xia
J.
,
Zhou
M.
,
Deng
S.
&
Li
Z.
2021
Adjustment of flood discharge capacity with varying boundary conditions in a braided reach of the Lower Yellow River
.
Journal of Geographical Sciences
31
(
11
),
1615
1632
.
Hu
C.
,
Guo
Q.
,
Chen
J.
&
Chen
X.
2006
Measures of reforming and maintaining the river channel for medium flood in lower Yellow River
.
Journal of Hydraulic Engineering
37
(
4
),
381
388
(in Chinese). https://doi.org/10.13243/j.cnki.slxb.2006.04.001
.
Hu
C.
,
Chen
J.
&
Guo
Q.
2012
Shaping and maintaining a medium-sized main channel in the Lower Yellow River
.
International Journal of Sediment Research
27
(
3
),
259
270
.
Kong
D.
,
Miao
C.
,
Borthwick
A. G. L.
,
Duan
Q.
,
Liu
H.
,
Sun
Q.
, Ye, A., Di, Z. &
Gong
W.
2015
Evolution of the Yellow River Delta and its relationship with runoff and sediment load from 1983 to 2011
.
Journal of Hydrology
520
,
157
167
.
Kong
D.
,
Miao
C.
,
Li
J.
&
Zheng
H.
2022
Full-stream erosion in the lower Yellow River: feasibility, sustainability and opportunity
.
Science of the Total Environment
807
,
150810
.
Li
J.
,
Xia
J.
,
Zhang
X.
&
Deng
S.
2015
Relationships between bankfull discharge-area and flow-sediment condition in Lower Yellow River under quasi-equilibrium
.
Journal of Sediment Research
2015 (
10
),
37
43
(in Chinese). https://doi.org/10.16239/j.cnki.0468-155x.2015.05.007
.
Ma
Y.
,
Huang
H.
,
Nanson
G. C.
,
Li
Y.
&
Yao
W.
2012
Channel adjustments in response to the operation of large dams: the upper reach of the lower Yellow River
.
Geomorphology
147–148
,
35
48
.
Miao
C.
,
Kong
D.
,
Wu
J.
&
Duan
Q.
2016
Functional degradation of the water–sediment regulation scheme in the lower Yellow River: spatial and temporal analyses
.
Science of the Total Environment
551–552
,
16
22
.
Qi
P.
,
Qi
H.
&
Tian
S.
2013
Great changes on the Lower Yellow River channel since 2000 and future prospective
.
Advances in Science and Technology of Water Resources
33
(
6
),
23
28
(in Chinese)
.
Shen, G., Zhang, Y. & Zhang, M. 2020 Characteristics on floods and sediment transport in the Lower Yellow River since the Xiaolangdi Reservoir impoundment. Journal of Sediment Research 45 (6), 59–66 (in Chinese).
Wang
H.
,
Wu
X.
,
Bi
N.
, Li, S., Yuan, P., Wang, A., Syvitski, J. P. M., Saito, Y., Yang, Z., Liu, S. &
Nittrouer
J.
2017
Impacts of the dam-orientated water–sediment regulation scheme on the lower reaches and delta of the Yellow River (Huanghe): a review
.
Global and Planetary Change
157
,
93
113
.
Wu
B.
,
Wang
G.
,
Xia
J.
,
Fu
X.
&
Zhang
Y.
2008a
Response of bankfull discharge to discharge and sediment load in the Lower Yellow River
.
Geomorphology
100
(
3
–4),
366
376
.
Wu
B.
,
Xia
J.
,
Fu
X.
,
Zhang
Y.
&
Wang
G.
2008b
Effect of altered flow regime on bankfull area of the Lower Yellow River, China
.
Earth Surface Processes and Landforms
33
(
10
),
1585
1601
.
Xia
J.
,
Wu
B.
,
Wang
Y.
& Li, W.
2010a
Estimating the bankfull discharge in the Lower Yellow River and analysis of its variation processes
.
Journal of Sediment Research
2010 (
2
),
6
14
(in Chinese). https://doi.org/10.16239/j.cnki.0468-155x.2010.02.003
.
Xia
J.
,
Wu
B.
,
Wang
G.
&
Wang
Y.
2010b
Estimation of bankfull discharge in the Lower Yellow River using different approaches
.
Geomorphology
117
(
1
–2),
66
77
.
Xia
J.
,
Li
X.
,
Zhang
X.
&
Li
T.
2014b
Recent variation in reach-scale bankfull discharge in the Lower Yellow River
.
Earth Surface Processes & Landforms
39
(
6
),
723
734
.
Xu
G.
&
Si
C.
2009
Effect of water and sediment regulation on lower Yellow River
.
Transactions of Tianjin University
15
(
2
),
113
120
.
Xu
G.
,
Zhang
J.
&
Lian
J.
2005
Effect of water–sediment regulation of the Yellow River on the lower reach
.
Advances in Water Science
16
(
4
),
518
523
(in Chinese)
.
Yang
G.
,
Wu
B.
,
Zhang
R.
&
Shang
H.
2020
Influence of ‘318’ controlled operation of the Sanmenxia Reservoir on variations of the Tongguan elevation
.
Journal of Sediment Research
45
(
3
),
38
45
(in Chinese). https://doi.org/10.16239/j.cnki.0468-155x.2020.03.006
.
Zhang
M.
,
Zhang
J.
,
Liu
M.
,
Zhao
H.
&
Jiang
G.
2009
Impact on discharge capacity by sediment–runoff regulation in downstream of Yellow River
.
Journal of Water Resources and Water Engineering
20
(
3
),
140
142
.
Zhang
J.
,
Shang
Y.
,
Liu
J.
,
Lu
J.
,
Wei
S.
,
Wan
Z.
, Luo, Q., Chen, C., Tong, L., Wang, Q. & Fu, J.
2021a
Optimisation of reservoir operation mode to improve sediment transport capacity of silt-laden rivers
.
Journal of Hydrology
594
,
125951
.
Zhang
J. L.
,
Lu
J.
, Wei, S., Luo, Q. & Wan, Z.
2021b
Causes and countermeasures of insufficient follow-up power for water and sediment regulation in Xiaolangdi Reservoir
.
Yellow River
43
(
1
),
5
9
(in Chinese)
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY-NC-ND 4.0), which permits copying and redistribution for non-commercial purposes with no derivatives, provided the original work is properly cited (http://creativecommons.org/licenses/by-nc-nd/4.0/).