Next Article in Journal
Carbon Emission Accounting and Reduction Evaluation in Sponge City Residential Areas
Previous Article in Journal
Impact of Chemical Oxygen Demand/Total Nitrogen Ratio on Shifting Autotrophic Partial Nitrification to Heterotrophic Nitrification and Aerobic Denitrification in High-Strength Ammonium Wastewater Treatment
Previous Article in Special Issue
Economic Evaluation of Water Management Alternatives in the Upper Green River Basin of Wyoming
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on the Threshold of the Transverse Gradient of the Floodplain in the Lower Yellow River Based on a Flood Risk Assessment Model

1
Henan Yellow River Engineering and Consulting Co., Ltd., Zhengzhou 450003, China
2
Henan Engineering Research Center for Protection and Governance of Yellow River, Zhengzhou 450003, China
3
School of Water Conservancy and Transportation, Zhengzhou University, Zhengzhou 450001, China
4
Yellow River Henan Bureau, YRCC, Zhengzhou 450003, China
*
Authors to whom correspondence should be addressed.
Water 2024, 16(17), 2533; https://doi.org/10.3390/w16172533
Submission received: 26 July 2024 / Revised: 28 August 2024 / Accepted: 3 September 2024 / Published: 6 September 2024
(This article belongs to the Special Issue Socio-Economics of Water Resources Management)

Abstract

:
Due to the influence of water and sediment conditions, engineering projects, channel erosion and siltation, river-related factors, and human activities (such as adjustments in floodplain production structures and village construction), there have been significant variations in the transverse gradient of the floodplain in the lower Yellow River. An irrational transverse gradient can lead to the rapid conversion of gravitational potential energy into kinetic energy during the flood evolution process, resulting in increased flow velocity and inundated areas. Exploring reasonable transverse gradients can provide technical support for floodplain management. Using “flood risk assessment” as a keyword, research papers from the Web of Science core database and CNKI published in the past five years were collected. Through a VOS viewer analysis of indicators, a flood risk assessment model based on the “Source–Path–Receptor–Consequence–Resilience” framework was established. A two-dimensional water and sediment model was used to simulate flood inundation scenarios with different transverse gradients in the same flood event, evaluate flood risks in the floodplain, and determine the optimal transverse gradient based on flood risk levels. The results indicate that, compared to low transverse gradients, moderate and high transverse gradients have a more significant driving effect on flood inundation, increasing flood risk opportunities for floodplains. Lower transverse gradients (i.e., TG = 10LG = 1.25‰) are the most favorable for flood protection in the floodplain after flood inundation.

1. Introduction

The lower Yellow River is a “hanging river” crossing the heartland of China, and flood safety is crucial for the stability of people’s lives in the floodplain [1]. The floodplain within the river channel serves as a place for major floods to flow, be detained, and deposit sediment while also being the homeland for millions of people, playing a dual role in flood safety along the Yellow River and social and economic development in the region [2]. The floodplains in the lower Yellow River are widespread in meandering and transitional river sections. Transitional sections have highly curved channels, poor connectivity, and significant changes in sediment transport with imbalanced water–sediment relationships due to a large sediment load, leading to massive sediment deposition and severe channel infilling. It has been nearly two decades since large-scale flooding occurred in the downstream floodplain, allowing residents to live and work in peace. However, human activities have significantly altered the conditions of the floodplain [3]. The downstream Yellow River represents a typical compound river channel, where the production embankment areas on both sides enclose the river channel, with vast floodplains between the main embankments [4]. Since the beginning of this century, changes in water and sediment flow and the regulation of the Xiaolangdi Reservoir have caused new alterations in the water and sediment characteristics entering the downstream channel. There has been a substantial decrease in sediment loads and a noticeable increase in low-flow processes [5]. Erosion in the downstream channel has intensified, with maximum average erosion depths reaching around 3 m. Inadequate alignment between channel management projects and water–sediment processes has led to substantial adjustments in local river stretches, resulting in deformed river profiles and instances of bank collapses and village destruction [6]. Despite not experiencing significant floods in over twenty years, changes are occurring on the floodplain due to human activities [7].
The transverse gradient of the floodplain surface is an important indicator reflecting the development of a “secondary hanging river” [8]. Similar to the water surface transverse gradient, the transverse gradient of the floodplain refers to the ratio of the elevation difference between the floodplain lip and the embankment toe to the horizontal distance between them [9]. The variation in transverse gradient largely depends on water and sediment conditions, channel and floodplain management projects, cross-sectional morphology, river-related factors, erosion and silting changes due to channel flooding, changes in floodplain land use, and flood routing along embankment toes. Influenced by water and sediment conditions, engineering projects, channel erosion and siltation, river-related parameters, and human activities, significant changes also occur in the transverse gradient of the floodplain [10]. This has a major impact on flood inundation, where a higher transverse gradient can lead to the rapid conversion of gravitational potential energy into kinetic energy during the flood evolution process, resulting in increased flow velocity and inundated areas. An appropriate transverse gradient not only facilitates material and energy exchange between the floodplain and the river channel and promotes vegetation growth in the floodplain but also helps avoid flood routing along the embankment during flood events. Chinese scholars have utilized measured cross-sectional data to analyze trends in transverse gradient variations, conducting spatial and temporal studies on the degree of change in transverse gradients. Zhang Jinliang [11] conducted a temporal and spatial analysis of changes in the transverse gradient on both sides of the High Cun-Taochengpu river section. They found that since 2000, there has been an increase in the transverse gradient, with the ratios concentrated between 0.5‰ and 2‰ along the course, where the transverse gradient of the left bank was slightly larger than that of the right bank. Liu Yan weighted the difference in height between the floodplain and channel along the Dongbatou-Gaocun section and found that the “secondary hanging river” appeared in 1972 and has persisted since then. It is widely recognized among domestic scholars that the development of secondary hanging rivers in the Gaocun-Taochengpu river section is severe [12]. Many scholars have investigated the impact of various factors on transverse gradients. Xu Linjuan studied flood processes and the development of secondary hanging rivers in meandering river sections, finding that a sand influx coefficient less than 0.1 and a flooding coefficient greater than 0.29 can effectively restrict the formation of secondary hanging rivers [13]. Yang Jishan, analyzing longitudinal cross-sectional data and embankment conditions, highlighted how embankments limit sediment accumulation areas and play a significant role in increasing the transverse gradient. In general, an increase in the transverse gradient is the result of adverse water and sediment conditions combined with boundary conditions [14]. Changes in water and sediment conditions are the primary reasons for a significant increase in the transverse gradient, while human activities within the river channel, particularly the construction of embankments, alter flow conditions, leading to a drastic rise in the transverse gradient, thus exacerbating the development of secondary hanging rivers [15]. The specific manifestations of transverse gradients include deep channels, low floodplains, and depressions at the base of embankments. During floods, this can easily trigger the formation of transverse flows, oblique flows, and rolling currents, resulting in flood overtopping of embankments and causing flood routing along embankments, posing a serious threat to flood defenses [14]. Moreover, the swift flow and deep water near embankment toes can lead to dangerous situations, such as seepage and piping at the back of large embankments. Sun Dongpo established a two-dimensional mathematical model for water and sediment transport, studying the effects of different suspended sediment loads on flood evolution in the floodplain, flow structures in the flood channel, and flood routing along embankments [16]. The results indicated that increasing suspended sediment loads could heighten the risk of embankment erosion. Once they exceeded critical values, they could significantly affect flow velocities along embankments. The results of an experimental project for controlling secondary hanging rivers in 2003 revealed that reducing the transverse gradient to some extent decreased the likelihood of flood inundation, flood routing along embankments, transverse and oblique flows, and rolling currents occurring in the study river section [17].
The evolution of floods refers to the process of flood waves propagating and evolving along river channels (or diversion areas or detention areas) [18]. After flooding on the floodplain, there are roughly two modes of evolution on the floodplain surface: the triangular floodplain exchange mode or the elongated floodplain exchange mode. The former indicates that water first enters from the triangular floodplain area, undergoes sediment deposition and accumulation on the floodplain, and then returns to the river channel. The latter mode involves water flowing into the floodplain at the breach connecting the river engineering structures and embankments, gradually advancing toward the near-bank floodplain, and resulting in flood routing along the embankment or longitudinal flooding along the floodplain. It then returns to the main river at the next weak point in the engineering structure or embankment. Conducting research on floodplain evolution requires calculations and studies on water and sediment transport processes in compound river channels [19]. To determine the characteristics of floodplain water flow evolution, scholars both domestically and internationally have conducted in-depth research. Wang Shudong, starting from the turbulent Reynolds transport equation, combining turbulence theory with floodplain water flow characteristics, derived analytical solutions for the two-dimensional distribution of point flow velocities [20]. Shinon and Knight proposed the SKM method, which integrates the turbulent Reynolds transport equation along the water depth. They derived the motion equations of floodplain water flow based on the Boussinesq assumption and the Darcy–Weisbach formula [21,22]. Zong Hucheng, in a study on solving floodplain water flow velocities, incorporated the factor of the transverse gradient into the model, validating the driving role of the transverse gradient in the floodplain evolution process [23]. Li Xizhi simulated the floodplain evolution process at a specific location in Langyuan using the Delft3D model. They observed that there was a buffer zone near the river channel before widespread floodplain inundation. Floodwaters first gathered in the vicinity of the river channel before advancing extensively onto the floodplain, entering a converging state prior to flood occurrence. At the middle stage of flood evolution, floodwaters spread along a certain trajectory within the floodplain. In the later stages of floodplain inundation, the floodwaters were insufficient to support the expansion of the flooded area [24]. Zhang Xiaolei constructed a two-dimensional mathematical model to explore the evolution of floodwaters from breached embankments on the floodplain with different roughness values. The results indicated that the time taken for the flood wave front to reach various measurement points in the floodplain increased gradually with higher roughness values and decreased gradually with increasing flood levels. Flow velocities at various measurement points in the floodplain showed an initial increase followed by a decrease, with the rate of increase in velocity being significantly greater than the rate of decrease [25]. Guo Peng established a two-dimensional mathematical model to simulate the progress of floodwaters resulting from embankment breaches on the floodplain [26]. The study demonstrated that, driven by the transverse gradient, water discharged from breaches spread laterally across the floodplain, eventually reaching the main levee of the Yellow River, posing a threat to its safety. When foreign scholars simulate flood evolution in floodplains, they generally do not consider sediment transport and bed deformation processes, with most of them utilizing diffusive wave models [27,28]. Caleffi used a two-dimensional shallow-water equation to simulate flood evolution in the Toce River [29]. Bates, based on high-precision digital elevation data, calculated the flood evolution process in the Severn River floodplain using a diffusive wave model [30]. Currently, research on floodplain evolution driven by the transverse gradient mainly focuses on theoretical derivation and experimental simulation. In terms of theory, studies have deduced the movement of floodwaters on the floodplain under the influence of the transverse gradient and analyzed flow velocity distributions.
Previous research has made some progress in understanding the lateral gradient of the floodplain and the flood evolution within the floodplain driven by this gradient. However, the terrain and geomorphology of the Yellow River’s downstream floodplain are heavily influenced by human activities, leading to complex variations. Most existing studies focus on the spatiotemporal distribution characteristics of the lateral gradient, its driving factors, and its impacts on flood management. Few scholars have explored how changes in the lateral gradient affect the progression of overbank flooding. Further research is needed to determine how to select an appropriate lateral gradient to accommodate changes in the floodplain during flood events. Flood risk assessment involves establishing an evaluation system with a multi-level, multi-indicator structure. This system standardizes each indicator value, assigns weights to the indicators, and performs superposition operations on them. Evaluating the differences in flood risk levels in different areas of the study region after these processes becomes a crucial method for assessing the extent of losses caused by floods in disaster-affected areas [31,32]. This study established a comprehensive and reliable flood risk assessment model to explore the impact of changes in the lateral gradient of the floodplain on the progression of overbank flooding. By utilizing the flood inundation results within the floodplain, this study provides theoretical support and a decision-making basis for optimizing the threshold values of the lateral gradient. The objective of this study is as follows: (1) establish a two-dimensional water and sediment model for the Gaocun–Sunkou river section to simulate flood evolution with different transverse gradients; (2) develop a flood risk assessment model suitable for floodplains and apply it to the Gaocun–Sunkou river section to calculate the spatial distribution of flood risk under different transverse gradients; (3) based on the results of the flood risk assessment, select a more reasonable transverse gradient.

2. Materials and Methods

2.1. An Overview of the Study Area

The Gaosun River section (located between 115°39′ to 116°08′ E and 35°43′ to 36°02′ N) spans the region between the Gaocun and Sunkou cross-sections, covering parts of Henan and Shandong provinces (Figure 1). This river segment is 118.2 km long with a channel slope of 1.11‰. It represents a transitional river section from a meandering to a sinuous morphology and is characterized by wide floodplains. The river consists of a compound channel formed by both the main channel and the floodplain. The evolution of the river shows dual characteristics [33], preserving the basic features of a meandering river, including scattered flow, significant swing amplitudes, rapid changes, and relatively straight channels. Simultaneously, it exhibits characteristics of sinuous rivers with single-channel bends and minor variations in the main channel’s curvature in the plan view. The river’s planform alternates between straight and curved sections, featuring mostly straight and elongated river segments such as Liuzhuang-Susizhuang and Yingfang-Penglou, each exceeding a length of 8 km. Due to its typical compound channel features, moderate to small flows pass through the main channel, while during major flood events, floodwaters need to inundate the floodplain, with overland flow occupying around 20% of the entire cross-section due to the presence of broad floodplains.
In the Gaosun River section, the channel is close to the main embankment, and there is a dense population of channel engineering structures [34]. The average distance between embankments on both sides is 4.5 km. The floodplain area between the production embankment and the main embankment is approximately 540 km2, with the left bank and right bank floodplain areas covering 321 km2 and 139 km2, respectively. The floodplain at the base of the embankments is low-lying, contains numerous sandbars, and has uneven surfaces with many depressions (Figure 2). Some areas of the floodplain have difficulty draining, leading to the potential formation of stagnation zones during inundation events, showcasing significant water storage capacity [2]. The average transverse gradients of the left and right bank floodplains are 9.80‰ and 10.39‰, respectively [35].

2.2. Data Processing

2.2.1. Natural Geographic Data

The Digital Elevation Model (DEM) data are sourced from the Shuttle Radar Topography Mission (SRTM) data by the National Aeronautics and Space Administration (NASA). It can be downloaded from the Chinese Academy of Sciences Geographic Spatial Data Cloud website (http://www.csdb.cn/index.html#, accessed on 1 September 2023), with a grid size of 30 m × 30 m. The data type is raster data.

2.2.2. Land Use Data

The land use data used in this study are based on the National Resource and Environment Database of the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences. It utilizes Landsat satellite remote sensing image data from the United States as the primary information source and establishes a multi-period land use/land cover remote sensing monitoring database (CNLUCC) through manual visual interpretation.

2.2.3. Remote Sensing Image Data

The Landsat series data are from the joint program by NASA and the U.S. Geological Survey, providing the longest continuous space-based observation records globally. The spatial resolution is 30 m, and it can be downloaded from the U.S. Geological Survey website (https://earthexplorer.usgs.gov/, accessed on 1 September 2023).

2.2.4. Hydrometeorological Data

To monitor changes in riverbeds, the Yellow River Conservancy Commission established fixed cross-sections along the lower reaches of the Yellow River in 1950. Professional personnel conduct cross-sectional measurements before and after the flood season each year. Additionally, they regularly monitor the water and sediment conditions of the river channel during both the flood and non-flood seasons according to different standards, documenting and organizing the data. The Yellow River Conservancy Commission publishes technical standards and checks the data to ensure measurement quality, focusing mainly on water level, flow rate, and sediment concentration data at the Gaocun and Sunkou cross-sections.

2.2.5. Cross-Sectional Data

A total of 38 hydrological cross-sections are established within the study area. Actual measured cross-sectional data before and after the flood season for each section are collected to reflect changes in the cross-sectional morphology at different times. The method of cross-section analysis was employed to calculate erosion and deposition within the river section.

2.3. Calculation of Transverse Gradient

In order to facilitate calculations, the river channel cross-sections were generalized for the Gaosun River, a compound channel characterized by a narrow and deep main channel with broad floodplains and a significant transverse gradient. Referencing the longitudinal gradient formula for rivers, the transverse gradient of the floodplain was computed as shown in Figure 3. In Figure 3, b1 represents the width of the river, and b2,L and b2,R are the widths of the left and right bank floodplains, respectively, while HL and b2,R indicate the elevation differences between the lip of the floodplain and the base of the embankment on the left and right banks. To minimize errors in transverse gradient calculations near the lip of the floodplain and the base of the embankment due to local variations, an average elevation of 100m was used for these locations. Therefore, the transverse gradient of the floodplain can be calculated based on Equations (1) and (2):
T G L = h L b 2 , L
T G R = h R b 2 , L

2.4. Flood Risk Assessment Model

2.4.1. Selection of Model Indicators

In order to construct a flood risk assessment model specifically for the Gaosun River section’s floodplain, the Web of Science (WOS) core collection database was used as the platform for collecting English literature research data, while the CNKI database served as the platform for collecting Chinese literature research data. The search term “Flood Risk Assessment” was employed, and the time span from 2018 to 2023 was selected to ensure the accuracy of the chosen indicators. Academic papers published since 2018 were filtered, with “Article” as the document type. Through manual screening and the elimination of less relevant subject entries, a total of 183 articles were gathered as the data source for knowledge graph analysis. Utilizing the professional literature analysis tool VOSviewer v1.6.21, this study conducted a visual analysis of the selected indicators from the retrieved literature to identify evaluation indicators used in flood risk assessment studies. By integrating the results of clustering and co-occurrence analyses of hotspots in the literature, appropriate evaluation indicators were selected to establish an assessment system tailored to the floodplain of the Gaosun River section. The data content related to the evaluation indicators was derived from the measured data of the Yellow River Institute of Hydraulic Research. It was preliminarily screened and deemed relatively detailed and reliable.
Based on the indicators drawn in the heat map shown in Figure 4, it is evident that current research on flood assessment is mainly focused on urban flood disasters, watershed floods, and coastal storm surges. These indicators can generally be categorized into four classes: source, path, acceptor, and consequence. The indicators related to flood propagation paths often pertain to natural attributes, such as topography and landforms, while those linked to flood acceptors mainly revolve around societal attributes like population structure. However, the Tan District of Gaosun River involves three major socioeconomic industries, and its post-disaster resilience, to some extent, reflects the level of flood risk in the district. Considering that the study area is the Tan District of Gaosun River, relevant indicators of resilience were incorporated into the aspects of source, flood propagation paths, flood acceptors, and consequences. Building upon the indicators presented in the heat map, additional indicators with characteristics unique to the district were included to establish an indicator system. A total of 14 indicators were selected, and the constructed indicator system is shown in Table 1.

2.4.2. Determination of Indicator Weights

Currently, there are three main types of comprehensive quantitative evaluation methods commonly used in domestic and international contexts [36]. These include methods based on deterministic and uncertain information for quantitative or qualitative evaluations [37,38], such as the Analytic Hierarchy Process (AHP) and Gray Relational Analysis. Statistical-law-based methods [39,40] include Principal Component Analysis and Cluster Analysis, among others. Methods based on goal programming [41,42] include the Technique for Order Preference by Similarity to an Ideal Solution (TOPSIS). Ultimately, the three types can be classified, according to objective and subjective factors, into subjective weighting methods, objective weighting methods, and combined weighting methods. The AHP is selected for subjective weighting as the subjective approach, and the CRITIC method (Criteria Importance Though Inter criteria Correlation) is chosen for objective weighting as the objective approach, followed by using game theory methods for combined weighting [43]. Utilizing this comprehensive weighting method not only leverages the information inherent in the data but also integrates the expertise of industry specialists, resulting in more credible weights [44].
This study employed a normalization method to standardize the raw data, aiming to reduce the differences between various indicator variables [45].
If x j is a positive indicator (“the larger, the better” type of indicator), then
x i j = n i j m i n ( n j ) m a x ( n j ) m i n ( n j )
If x j is a negative indicator (“the smaller, the better” type of indicator), then
x i j = max ( n j ) n i j max ( n j ) min ( n j )
The formula for the AHP is as follows:
A = a i j n × n = a 11 a 1 n a n 1 a n n
W i ¯ = j = 1 n a i j n i , j = 1,2 , n
W i = W i ¯ i = 1 n W i ¯ i = 1,2 , n
where a i j is the ratio scale, indicating the relative importance of element i compared to element j at the same level; W i ¯ represents the n-th root of the product of the ratios; and Wi is the weight value.
The AHP consistency test is performed as follows:
C I = λ m a x n 1
C R = C I R I
where λ m a x denotes the maximum eigenvalue of the judgment matrix; CI is the consistency index; CR is the consistency ratio; and RI is the average random consistency index. When CR < 0.1, the judgment matrix A is considered to have satisfactory agreement.
The formula for CRITIC is as follows:
σ j = i = 1 e ( x i j x j ¯ ) e 1
f j = i = 1 e ( 1 r i j )
C j = σ j f j
w j = C j j = 1 f C j
In the formulas, σ j represents the distinctiveness of the j-th criterion, f j denotes the discordance of criterion j with the other criteria, r i j is the Pearson correlation coefficient between criteria i and j, C j signifies the information load of the j-th criterion, and w j denotes the weight.
The process of the integrated weighting calculation based on game theory is as follows:
W ¯ = a 1 w 1 + a 2 w 2 ,     a 1 + a 2 = 1 ,   a 1 > 0 ,   a 2 > 0
m i n i = 1 2 a i w i T w j T | | 2
w 1 w 1 T w 1 w 2 T w 2 w 1 T w 2 w 2 T a 1 a 2 = w 1 w 1 T w 2 w 2 T
where a 1 and a 2 represent the combination coefficients, and W denotes the comprehensive weight of the indicators.

2.5. Two-Dimensional Water–Sediment Model

In a flood risk assessment system, indicator data related to the source, the path, the acceptor, and resilience can be obtained from remote sensing images and statistical yearbooks. These data can be processed using spatial analysis techniques to convert them into raster data for flood risk assessment. Indicator data related to the consequences of flooding require conducting flood inundation simulations through the construction of a two-dimensional water–sediment model [46]. Below are the basic control equations and solution methods for the two-dimensional water–sediment model, along with an explanation of how key issues in the model are addressed, aiming to provide technical support for flood risk assessment in the Tan District of Gaosun River.

2.5.1. Hydrological Data

The Yellow River has not witnessed a large-scale flooding event in the floodplain for nearly 20 years. In recent years, a significant flooding incident occurred during the autumn flood of 2003, mainly affecting the Landon floodplain area. Detailed records are available for this event. In 1996, a major flood inundated the downstream floodplain areas of the Yellow River, impacting as many as 804 villages, with 830,000 affected individuals and 1.46 million mu (approximately 243,333 acres) of flooded farmland, resulting in significant economic losses to the floodplain region. The hydrological data collected during the 1996 flood event are quite comprehensive, with complete flow, water level, and sediment concentration data available at the Gao Village and Sunkou sections. These datasets meet the requirements for conducting two-dimensional water–sediment modeling. Additionally, Landsat 5 satellite imagery offers high-precision remote sensing images that vividly depict the inundation situation in the Gaosun River section of the floodplain region during this flood event. This imagery can be used to validate the model accuracy. Please refer to Figure 5 for the satellite image data.

2.5.2. Roughness Settings

According to previous research findings and the “Hydrological Yearbook of the People’s Republic of China” (hydrological data for the Yellow River Basin), Table 2 lists the roughness coefficients for various land use types.

2.5.3. Grid Division

In numerical simulations, the computational grid typically represents the discretization of irregular physical calculation regions. The number of grids, geometric shapes, grid quality, density, etc., all directly affect the accuracy and efficiency of numerical analysis. Based on the shape generated by generating 2D planar grids, they can be divided into quadrilateral grids and triangular grids. Triangular grids are suitable for complex regions. In the model, structured grids or unstructured grids, or quadrilateral grids or triangular grids, can be used based on the complexity of the physical region’s boundaries. Considering the morphology of the Gao Sun River section, triangular grids were selected. The determination of the grid size needs to be comprehensively considered from the perspective of both model accuracy and computational performance. Smaller grids result in higher model accuracy but longer simulation times, requiring higher computational performance. Due to different simulation focuses, different grid sizes were chosen for local refinement for the river channel and floodplain in the study area. The grid is shown in Figure 6.

2.6. Optimal Cross-Slope Selection Scheme

In order to express the severity of the transverse gradient of the beach surface, experts and scholars often use multiple relationships between the transverse gradient and longitudinal gradient. To differentiate between the impacts of changes in the transverse gradient on hydrological connectivity in the Gao Sun River section, this study categorizes the transverse gradient into three levels, namely, a low transverse gradient, a medium transverse gradient, and a high transverse gradient, based on the ratio relationship between the longitudinal gradient and transverse gradient of the river section. The classification criteria are shown in Table 3.
Calculations showed that the longitudinal gradient of the Gao Sun River section is 0.125‰. Based on the calculated results of the transverse gradient in the Gao Sun River section, the ratio relationship between the transverse gradient and the river channel’s longitudinal gradient was computed. Different simulation scenarios were then set up based on this ratio relationship. As shown in Table 4, during the flood inundation simulation process in the Gao Sun River section, the transverse gradient was sequentially set to low, medium, and high levels according to the ratio relationship between the transverse gradient and the longitudinal gradient for flood simulations.

3. Results

3.1. The Trend of TG Variation in the Floodplain

Based on the measured cross-sectional data of the Gaosun River before and after the flood, the TG of the river section from 1986 to 2021 was calculated. The calculation results are shown in Figure 7. The average transverse gradients before and after the flood over the past 35 years were 1.271‰ and 1.272‰, respectively. The long-term average transverse gradient before the flood is slightly lower than after the flood. Before the flood, the long-term average transverse gradients in the early, middle, and late stages of the Xiaolangdi Reservoir construction were 1.255‰, 1.057‰, and 1.369‰, respectively. After the dam was built, the average transverse gradient increased by 9.08% compared to before the dam construction. After the flood, the average transverse gradients in the early, middle, and late stages of the Xiaolangdi Reservoir construction were 1.221‰, 1.083‰, and 1.384‰ respectively, with an increase of 13.35% in the transverse gradient after dam construction compared to before. The operation of the Xiaolangdi Reservoir changed the water and sediment characteristics of the river section, affecting the erosion and sedimentation evolution of the river section and reducing the probability of flooding in the floodplain. This results in more sediment deposition at the flood lip, leading to an increase in the transverse gradient of the floodplain.
The operation of the Xiaolangdi Dam has increased the transverse gradient of the Gaosun River section. The change cycles and trends of the transverse gradient before and after the flood are quite consistent. Before 1997, it was in a state of fluctuating decrease, but it rapidly increased and stabilized thereafter. The main difference between the transverse gradients before and after the flood is that the interannual variability of the transverse gradient after the flood is more unstable. This is because the transverse gradient after the flood is more easily influenced by flood erosion and sediment deposition.

3.2. Optimal Selection of TG of the Floodplain

3.2.1. Flood Inundation Results

Based on the low, medium, and high transverse gradient scenarios described in Table 4, a two-dimensional water–sediment model was constructed to simulate flood inundation in different scenarios. Figure 8, Figure 9 and Figure 10 show the distribution of the total water depth during the evolution of flood inundation under different transverse gradient conditions. For comparison purposes, only the flooded areas are displayed.
Based on the simulation results, compared to the low transverse gradient scenario, the inundation area with a higher water depth slightly increases in the medium transverse gradient scenario, while the inundation area with a lower water depth decreases. Additionally, the inundation water depth at the embankment toe is significantly larger in the medium transverse gradient scenario compared to the floodplain results in the low transverse gradient scenario. On the other hand, the increase in the inundation area is more pronounced in the high transverse gradient scenario, especially in the area where the water depth exceeds 1.8 m. The areas with higher inundation water depths driven by medium and high transverse gradients are mainly concentrated at the embankment toe. It is easier for flooding to occur along the embankment in the medium and high transverse gradient scenarios compared to the low transverse gradient scenario. This is because, after floodplain inundation, the gravitational driving force accelerates the flow due to topographical effects, leading to higher flow velocities.

3.2.2. Determination of Indicator Weights

The combined weights of each evaluation indicator were calculated based on the calculation steps described earlier, as shown in Table 5.

3.2.3. Flood Risk Zoning

For different transverse gradients, the flood risk is mainly analyzed by changing the transverse gradient to analyze flood inundation. This involves altering the inundation depth, duration, and extent in the results, achieving the flood risk assessment of the floodplain with different transverse gradients for the same event (i.e., the autumn flood of 2003). Visualizing the flood risk level in the floodplain provides a basis for optimizing the threshold of the transverse gradient. Based on the comprehensive weights mentioned above and actual data, the flood risk assessment results were calculated using a raster calculator, as shown in Figure 11, Figure 12 and Figure 13.
To evaluate the flood risk in low, medium, and high scenarios, it is essential to assess the area proportions of different risk levels and analyze the changes in the areas of medium- to high-risk zones to determine the optimal transverse gradient for the three scenarios. By considering these two perspectives, one can make informed decisions on the most suitable transverse gradient for each scenario based on the distribution of risk levels and changes in medium- to high-risk areas.
(1) Assessment of flood risk with low TG
With a low transverse gradient, medium- and high-risk levels are mainly distributed in the river channel, while the flood risk level in the floodplain is mostly low or below. According to the statistical results, in the Gao Sun River segment, there are 43 blocks classified as high risk and 3 blocks as very high risk, accounting for only 7.64% of the total river segment area. These high-risk areas are primarily located within the river channel, with lower flood risks near the embankments, resulting in relatively low exposure to flood risks. The majority of the study area falls within the medium risk level or below, with the largest area classified as low risk, covering 55.98% of the river segment area. During the autumn flood of 2003, the flood risk faced by the floodplain with a low transverse gradient was relatively low. Combined with the land use conditions mentioned in Section 3.2.4, areas where residents live and work, such as developed land and farmland, face minimal threats from floods, which is beneficial for the economic development of the floodplain region.
(2) Assessment of flood risk with medium TG
Compared to the floodplain conditions with a low transverse gradient, the areas classified as extremely low risk and low risk decrease with a moderate gradient. Meanwhile, the areas classified as medium risk and high risk increase by 18.73%. With a moderate transverse gradient, the low-risk and medium-risk areas constitute the largest proportions in the total river segment area, accounting for 44.11% and 42.16%, respectively, while the high-risk area proportion is 11.27%, indicating an increase of 4.23% compared to the low transverse gradient scenario. Apart from the river channel, significant areas of high-risk and very high-risk zones are observed in the Gao Cun section and near the Sun Kou section. The area near the Sun Kou section has densely developed land, making it susceptible to flood threats. Under the influence of a moderate transverse gradient, after flooding overflows into the floodplain, some areas experience increased flood risks or even transition to higher-risk categories. This situation impacts flood control measures in the vicinity of the Gao Cun section’s and the Sun Kou section’s floodplain areas.
(3) Assessment of flood risk with high TG
With the high transverse gradient, the area of high-risk zones in the river segment dramatically increases, accounting for 26.35% of the total area. The high-risk areas are 3.74 times greater with high transverse gradients compared to those with low transverse gradients. High-risk zones and very high-risk zones also appear in the mid-section floodplain of the river segment. These high-risk areas are characterized by flat terrain, close proximity to the river channel, high population density, developed socioeconomic conditions, and predominantly agricultural land use. As these areas serve as flood buffer zones, they are highly susceptible to adverse effects from flooding. Consequently, they face significant flood control pressures.

3.2.4. The Optimization of the TG for the Floodplain

The flood risk assessment results in different scenarios indicate that, compared to a low transverse gradient, the gravitational driving force for flood inundation is more pronounced with moderate and high transverse gradients, increasing the flood risk faced by the floodplain. For floods of similar magnitudes, the areas classified as low risk and below with moderate and high transverse gradients decreased by 18.73% and 27.66%, respectively, compared to those with a low transverse gradient. The areas classified as medium risk increased by 14.50% and 8.02% with moderate transverse gradients compared to those with low transverse gradients. Additionally, the high-risk and above areas increased by 4.23% and 19.64% with high transverse gradients compared to those with low transverse gradients. Therefore, considering the different scenarios analyzed during the research process, it is most advantageous for floodplain flood safety to have a low transverse gradient (i.e., TG = 10LG = 1.25‰) in place when flash flooding occurs.

4. Conclusions

Based on historical cross-sectional data from the downstream Gao Sun River segment of the Yellow River, calculations and analyses of the transverse gradient trends on the floodplain were conducted. A flood risk assessment model for the Gao Sun River segment was constructed with a two-dimensional hydro-sediment model as the core. Different scenarios of transverse gradients were established to simulate floods and evaluate flood risk levels in the floodplain areas. A comparative analysis of flood risk areas with low, moderate, and high transverse gradients was carried out to determine the most favorable transverse gradient for flood protection in the floodplain. The results indicate the following:
(a) Over the past 35 years, the average pre-flood season and post-flood season transverse gradients on the floodplain were 1.271‰ and 1.272‰, respectively. Due to the operation of the Xiaolangdi Reservoir, the transverse gradient increased by 9.08% before the flood season and increased by 13.35% after the flood season.
(b) The mechanism by which water–sediment conditions and sedimentation influence the transverse gradient differs from the impact of river channel sedimentation control projects. Water overflows onto the floodplain, causing a rapid decrease in sediment-carrying capacity, leading to significant deposition near the flood lip. Subsequently, the difference between the sediment concentration and the sediment-carrying capacity gradually decreases, resulting in reduced sediment deposition. Water–sediment conditions and channel sedimentation directly impede sediment exchange on the floodplain, while sediment control projects alter the path of water–sediment evolution on the floodplain, consequently changing local sedimentation patterns and influencing transverse gradient changes.
(c) With moderate to high transverse gradients, there is a higher likelihood of floodwaters flowing along embankments after overflowing onto the floodplain, posing a threat to levee safety. For floods of similar magnitude, the areas classified as low risk and below decreased by 18.73% and 27.66% with moderate and high transverse gradients compared to those with low transverse gradients. The medium-risk areas increased by 14.50% and 8.02%, while high-risk and above areas increased by 4.23% and 19.64% with moderate and high transverse gradients. A low transverse gradient (TG = 10LG = 1.25‰) is deemed more favorable for flood protection.
This paper focuses on the Gaosun River section of the lower Yellow River, analyzing the trends in lateral gradient variation in the river’s floodplains and their driving factors. Subsequently, considering flood prevention factors in the floodplain area, a threshold optimization scheme for the lateral gradient based on a flood risk assessment model is proposed. The scheme is evaluated by comparing the methodology itself and the optimization results, which appear to be relatively reasonable. However, due to limitations in time and resources, there are certain constraints in the current research. Future studies could improve on these results by incorporating flood risk assessments for different flood events and lateral gradient scenarios, thereby enhancing the continuity of scenario settings in the research process. In practical applications, it is necessary to conduct threshold studies of the lateral gradient along different sections of the lower Yellow River. In future research, the lateral gradient of major cross-sections can replace the average lateral gradient of the river sections. Assessing floodplain spillover risk using different flood frequencies allows for timely forecasting and warning. Flood risk assessments across various ranges of the lateral gradient can more precisely identify the most effective gradients conducive to flood prevention in floodplain areas.

Author Contributions

Conceptualization, L.Q. and Z.Z.; methodology, Z.Z.; software, C.N. and M.L.; validation, Z.Z.; formal analysis, Z.Z.; investigation, L.Q.; resources, C.N. and C.H.; data curation, L.Q. and M.L.; writing—original draft preparation, Z.Z.; writing—review and editing, C.N., M.L. and G.A.; visualization, G.A. and M.L.; supervision, C.H.; project administration, C.H.; funding acquisition, Z.Z. and C.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “Outstanding Youth scientific research project of Yellow River Conservancy Commission of the Ministry of Water Resources, grant number HQK-202302”, “State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research, grant number IWHR-SKL-202104”, and “the National Natural Science Foundation of China, grant number No. U2243219”.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

Author Zheng Zhao was employed by the company Henan Yellow River Engineering and Consulting Co., Ltd., Zhengzhou, Henan Engineering Research Center for Protection and Governance of Yellow River, Zhengzhou. Author Ai Guangzhang was employed by the company Yellow River Henan Bureau, YRCC, Zhengzhou The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Xia, J.; Wang, Y.; Zhou, M.; Deng, S.; Li, Z.; Wang, Z. Variations in channel centerline migration rate and intensity of a braided reach in the Lower Yellow River. Remote Sens. 2021, 13, 1680. [Google Scholar] [CrossRef]
  2. Hu, C. Evolution of channel roughness and influencing factors in the wide floodplain of the Lower Yellow River in recent 20 years. J. Basic Sci. Eng. 2024, 32, 984–999. [Google Scholar]
  3. Zheng, L.; Wang, Y.; Li, J. Quantifying the spatial impact of landscape fragmentation on habitat quality: A multi-temporal dimensional comparison between the Yangtze River Economic Belt and Yellow River Basin of China. Land Use Policy 2023, 125, 106463. [Google Scholar] [CrossRef]
  4. Wang, Y.; Xia, J.; Deng, S.; Zhou, M.; Wang, Z.; Xu, X. Numerical simulation of bank erosion and accretion in a braided reach of the Lower Yellow river. Catena 2022, 217, 106456. [Google Scholar] [CrossRef]
  5. Miao, C.; Kong, D.; Wu, J.; Duan, Q. Functional degradation of the water–sediment regulation scheme in the lower Yellow River: Spatial and temporal analyses. Sci. Total Environ. 2016, 551, 16–22. [Google Scholar] [CrossRef]
  6. Peng, J.; Chen, S. Response of delta sedimentary system to variation of water and sediment in the Yellow River over past six decades. J. Geogr. Sci. 2010, 20, 613–627. [Google Scholar] [CrossRef]
  7. Dong, J.; Xia, X.; Wang, M.; Lai, Y.; Zhao, P.; Dong, H.; Zhao, Y.; Wen, J. Effect of water–sediment regulation of the Xiaolangdi Reservoir on the concentrations, bioavailability, and fluxes of PAHs in the middle and lower reaches of the Yellow River. J. Hydrol. 2015, 527, 101–112. [Google Scholar] [CrossRef]
  8. Wu, J.; Shang, K. Regional flood resilience evaluation of the wandering river basin of the lower Yellow River in Henan Province, China—Based on interval intuitionistic fuzzy numbers and the discrete Choquet integrals. J. Taiwan Inst. Chem. Eng. 2023, 105161. [Google Scholar] [CrossRef]
  9. Annis, A.; Karpack, M.; Morrison, R.R.; Nardi, F. On the influence of river Basin morphology and climate on hydrogeomorphic floodplain delineations. Adv. Water Resour. 2022, 159, 104078. [Google Scholar] [CrossRef]
  10. Wang, N.; Zhang, C.; Xiao, Y.; Jin, G.; Li, L. Transverse hyporheic flow in the cross-section of a compound river system. Adv. Water Resour. 2018, 122, 263–277. [Google Scholar] [CrossRef]
  11. Jinliang, Z.; Junzheng, L.; Yuchuan, B.; Haiyu, X.; Yan, L. Spatiotemporal water-sediment variations and geomorphological evolution in wide-floodplain transitional reach of lower Yellow River. J. Hydroelectr. Eng. 2021, 40, 1–12. [Google Scholar] [CrossRef]
  12. Yan, L.; Junhua, L.; Xiang, Z. Experimental Study on Low Land Area Control in the Secondary Suspended River Reach of the Lower Yellow River. Yellow River 2022, 44, 41–44. [Google Scholar]
  13. Xu, L.; Li, J.; Xu, H.; Zhang, X.; Lai, R.; Zhang, X.; Gao, X. Evolution and drivers of secondary suspended rivers in typical wandering sections of the lower Yellow River from 1960–2021. Front. Ecol. Evol. 2023, 11, 1330749. [Google Scholar] [CrossRef]
  14. Jishan, Y.; Jiongxin, X.; Jianhua, L. The Process of Secondary Suspended Channel in the Lower Yellow River under Differ ent Conditions of Runoff and Sediment Load. Acta Geogr. Sin. 2006, 61, 66–76. [Google Scholar]
  15. Zheng, Z.; Wu, T.; Cui, T. Application of 2-D hydrodynamic mathematical model in the regulation project of secondary suspended river. In Proceedings of the 3rd International Conference on Mechatronics, Robotics and Automation, Shenzhen, China, 14–15 May 2015; pp. 209–212. [Google Scholar]
  16. Dongpo, S.; Xiaoping, L.; Hai, X.; Pegtao, W.; Xiolong, L. Numerical simulation of fluvial process in the Lower Yellow River with “secondary perched river”. J. Hydroelectr. Eng. 2008, 27, 136–141. [Google Scholar]
  17. Yongwei, Z. “Second Hangs the River” Prevented Flood to the Low Reaches of the Yellow River Influence and the Countermeasure; Wu Han University: Wuhan, China, 2004. [Google Scholar]
  18. Qiu, W.; Li, Y.; Zhang, Y.; Wen, L.; Wang, T.; Wang, J.; Sun, X. Numerical investigation on the evolution process of cascade dam-break flood in the downstream earth-rock dam reservoir area based on coupled CFD-DEM. J. Hydrol. 2024, 635, 131162. [Google Scholar] [CrossRef]
  19. Kejun, Y.; Xingnian, L.; Shuyou, C.; Zhixiang, Z. Turbulence characteristics of overbank flow in compound river channel with vegetated floodplain. J. Hydraul. Eng. 2005, 36, 1263–1268. [Google Scholar] [CrossRef]
  20. Hanxiang, X. Simplified Calculation of Diffuse Beach Flow. Res. Waterborne Transp. 1982, 28, 84–92. [Google Scholar] [CrossRef]
  21. Shudong, W. Two-Dimensional Velocity Distribution and Hydraulic Calculation of Diffuse Beach Flow. J. Hydraul. Eng. 1986, 36, 51–59. [Google Scholar] [CrossRef]
  22. Shiono, K.; Knight, D.W. Turbulent open-channel flows with variable depth across the channel. J. Fluid Mech. 1991, 222, 617–646. [Google Scholar] [CrossRef]
  23. Rathor, S.K.; Mohanta, A.; Patra, K. Validation of Computational Fluid Dynamics Approach of Lateral Velocity Profile Due to Curvature Effect on Floodplain Levee of Two-stage Meandering Channel. Water Resour. Manag. 2022, 36, 5495–5520. [Google Scholar] [CrossRef]
  24. Hucheng, Z.; Enhui, J.; Lianjun, Z.; Xiaoxue, Z. The esearch of Two Dimensional Analytical Solution for Overbank Flow Velocity on the Floodplain Transverse Slope. Yellow River 2015, 37, 45–49. [Google Scholar] [CrossRef]
  25. Xizhi, L.; Huan, Z.; Zhenshan, L.; Peng, H. Analysis Method of Inflow and Outflow Evolution of Flood Inundation in the Lower Yellow River Floodplain. Yellow River 2023, 45, 77–84. [Google Scholar] [CrossRef]
  26. Xiaolei, Z.; Junqiang, X.; Peng, G.; Qian, C. Experimental and Numerical Investigations of Farm Dike-break Induced Floods. Adv. Eng. Sci. 2018, 50, 71–81. [Google Scholar] [CrossRef]
  27. Peng, G.; Junqiang, X.; Meirong, Z. Numerical modelling of farm dike lateral breach induced by flooding processes. J. Sediment Res. 2022, 47, 15–22. [Google Scholar] [CrossRef]
  28. Horritt, M.S.; Bates, P.D.; Mattinson, M.J. Effects of mesh resolution and topographic representation in 2D finite volume models of shallow water fluvial flow. J. Hydrol. 2006, 329, 306–314. [Google Scholar] [CrossRef]
  29. Yu, D.; Lane, S.N. Urban fluvial flood modelling using a two-dimensional diffusion-wave treatment, part 1: Mesh resolution effects. Hydrol. Process. Int. J. 2006, 20, 1541–1565. [Google Scholar] [CrossRef]
  30. Caleffi, V.; Valiani, A.; Zanni, A. Finite volume method for simulating extreme flood events in natural channels. J. Hydraul. Res. 2003, 41, 167–177. [Google Scholar] [CrossRef]
  31. Bates, P.D.; Wilson, M.D.; Horritt, M.S.; Mason, D.C.; Holden, N.; Currie, A. Reach scale floodplain inundation dynamics observed using airborne synthetic aperture radar imagery: Data analysis and modelling. J. Hydrol. 2006, 328, 306–318. [Google Scholar] [CrossRef]
  32. Cossart, É.; Fressard, M. Assessment of structural sediment connectivity within catchments: Insights from graph theory. Earth Surf. Dyn. 2017, 5, 253–268. [Google Scholar] [CrossRef]
  33. Turnbull, L.; Wainwright, J.; Brazier, R.E. A conceptual framework for understanding semi-arid land degradation: Ecohydrological interactions across multiple-space and time scales. Ecohydrol. Ecosyst. Land Water Process Interact. Ecohydrogeomorphol. 2008, 1, 23–34. [Google Scholar] [CrossRef]
  34. Yuchuan, B.; Yan, L.; Jinliang, Z.; Yang, B.; Ziqing, J. Energy dissipation of boundary resistance and stability analysis of riverbed of Gaocun to Tao-chengpu reach in lower Yellow River. J. Hydraul. Eng. 2020, 51, 1165–1174. [Google Scholar] [CrossRef]
  35. Hongwu, Z. Sanmenxia Reservoir on Yellow River Should Shine Again in the New Era. Yellow River 2022, 44, 1–4. [Google Scholar] [CrossRef]
  36. Wenyi, Y.; Wenxue, L.; Zhijun, H.; Wenhua, C. Channel shrinkage and its disaster-causing mechanism in the Lower Yellow River. J. Hydraul. Eng. 2005, 36, 257–264. [Google Scholar]
  37. Zhanlin, P.; Qiang, Z.; Shanlin, Y. Overview of Comperhensive Evaluation Theory and Methodology. Chin. J. Manag. Sci. 2015, 23, 245–256. [Google Scholar]
  38. Aidinidou, M.T.; Kaparis, K.; Georgiou, A.C. Analysis, prioritization and strategic planning of flood mitigation projects based on sustainability dimensions and a spatial/value AHP-GIS system. Expert. Syst. Appl. 2023, 211, 118566. [Google Scholar] [CrossRef]
  39. Han, S.; Li, D.; Li, K.; Wu, H.; Gao, Y.; Zhang, Y.; Yuan, R. Analysis and Study of Transmission Line Icing Based on Grey Correlation Pearson Combinatorial Optimization Support Vector Machine. Measurement 2024, 236, 115086. [Google Scholar] [CrossRef]
  40. Liu, J.; Kang, H.; Tao, W.; Li, H.; He, D.; Ma, L.; Tang, H.; Wu, S.; Yang, K.; Li, X. A spatial distribution—Principal component analysis (SD-PCA) model to assess pollution of heavy metals in soil. Sci. Total Environ. 2023, 859, 160112. [Google Scholar] [CrossRef]
  41. Zha, S.; Jin, Y.; Wheeler, R.; Bosarge, E. A mixed-method cluster analysis of physical computing and robotics integration in middle-grade math lesson plans. Comput. Educ. 2022, 190, 104623. [Google Scholar] [CrossRef]
  42. Sahoo, D.; Parida, P.K.; Pati, B. Efficient fuzzy multi-criteria decision-making for optimal college location selection: A comparative study of min–max fuzzy TOPSIS approach. Results Control Optim. 2024, 15, 100422. [Google Scholar] [CrossRef]
  43. Nafei, A.; Azizi, S.P.; Edalatpanah, S.A.; Huang, C.-Y. Smart TOPSIS: A Neural Network-Driven TOPSIS with Neutrosophic Triplets for Green Supplier Selection in Sustainable Manufacturing. Expert Syst. Appl. 2024, 255, 124744. [Google Scholar] [CrossRef]
  44. Cheng, G.; Li, G.; Pu, X.; Chen, C.; He, Y. Advancing coupling coordination simulation in the social-human-ecological system of the Three Gorges Reservoir Area: A multi-scenario system dynamics approach. Ecol. Indic. 2024, 158, 111504. [Google Scholar] [CrossRef]
  45. Ruiying, W.; Huaiwei, S.; Dong, Y.; Hui, T.; Weihong, L.; Gaorui, C.; Dongwei, G. Evaluation of flood disaster risk in China-Pakistan Economic Corridor by combination weighting based on improved game theory and grid data. Trans. Chin. Soc. Agric. Eng. 2021, 37, 145–154. [Google Scholar] [CrossRef]
  46. Li, M.; Niu, C.; Li, X.; Quan, L.; Li, W.; Liu, C.; Shi, C.; Soomro, S.-e.-h.; Song, Q.; Hu, C. Modeling and Evaluating the Socio-Economic–Flood Safety–Ecological System of Landong Floodplain Using System Dynamics and the Weighted Coupling Coordination Degree Model. Water 2024, 16, 2366. [Google Scholar] [CrossRef]
Figure 1. Topography and distribution of key cross-sections in the Gaosun River section.
Figure 1. Topography and distribution of key cross-sections in the Gaosun River section.
Water 16 02533 g001
Figure 2. A flow structure diagram of the floodplain channel in the lower Yellow River.
Figure 2. A flow structure diagram of the floodplain channel in the lower Yellow River.
Water 16 02533 g002
Figure 3. Generalized cross-section diagram of the Gaosun River section.
Figure 3. Generalized cross-section diagram of the Gaosun River section.
Water 16 02533 g003
Figure 4. Indicator hotspot map based on VOSviewer.
Figure 4. Indicator hotspot map based on VOSviewer.
Water 16 02533 g004
Figure 5. Image of the submerged area in the Gao Sun River section in 1996.
Figure 5. Image of the submerged area in the Gao Sun River section in 1996.
Water 16 02533 g005
Figure 6. The 2D water–sediment model grid.
Figure 6. The 2D water–sediment model grid.
Water 16 02533 g006
Figure 7. The results of the calculation for the transverse gradient of the floodplain from 1986 to 2021.
Figure 7. The results of the calculation for the transverse gradient of the floodplain from 1986 to 2021.
Water 16 02533 g007
Figure 8. The distribution map of the total water depth during flooding in the low TG.
Figure 8. The distribution map of the total water depth during flooding in the low TG.
Water 16 02533 g008
Figure 9. The distribution map of the total water depth during flooding in the medium TG.
Figure 9. The distribution map of the total water depth during flooding in the medium TG.
Water 16 02533 g009
Figure 10. The distribution map of the total water depth during flooding in the high TG.
Figure 10. The distribution map of the total water depth during flooding in the high TG.
Water 16 02533 g010
Figure 11. Flood risk assessment of low TG.
Figure 11. Flood risk assessment of low TG.
Water 16 02533 g011
Figure 12. Flood risk assessment of medium TG.
Figure 12. Flood risk assessment of medium TG.
Water 16 02533 g012
Figure 13. Flood risk assessment of high TG.
Figure 13. Flood risk assessment of high TG.
Water 16 02533 g013
Table 1. Construction table of evaluation indicators.
Table 1. Construction table of evaluation indicators.
CriteriaIndicatorExplanationAttributes
SourceWaterbodyIdentifies whether the water body within the area is classified as a water body to determine whether it would serve as a source of flooding during inundation. The value can be either 1 or 0.Positive
Distance from channelBy calculating the distance of each point within the area to the river channel, the likelihood of inundation at each point can be determined based on the proximity to the river. Generally, the farther a point is from the river channel, the less likely it is to be inundated.Negative
PathElevationBy assessing the elevation of each point within the area, it can be determined whether it is susceptible to flooding. Generally, lower-elevation points are more prone to inundation.Negative
RoughnessA dimensionless parameter reflecting the influence on water flow resistance. The rougher the boundary surface, the higher the roughness coefficient, resulting in slower water flow; conversely, the smoother the boundary surface, the lower the roughness coefficient, leading to faster water flow.Positive
AcceptorNDVIA standardized index used to generate images displaying the vegetation amount (relative biomass). It can reflect losses in inundated areas during floods.Positive
NDWIThe NDWI is typically used to extract water body information from images, reflecting water bodies within inundated areas and displaying flood losses.Positive
Population densityPopulation density is the number of people per unit area of land and is an important indicator for measuring the distribution of population in inundated areas.Positive
ImperviousnessA crucial indicator for identifying impermeable surfaces and can reflect the infiltration situation in inundated areas after floodplain inundation.Positive
GDP per unitGross Domestic Product (GDP) per unit area in the study area represents the economic status of the inundated area.Positive
Nighttime light indexThe nighttime light index is based on the image data of human nighttime activities and production extracted using satellite remote sensing and data analysis techniques. Economically developed and densely populated areas often shine brightly at night, resulting in a high nighttime light index. Conversely, economically underdeveloped and sparsely populated areas exhibit dim or no nighttime lights, leading to a low nighttime light index.Positive
ConsequenceFloodwater depthThe data derived from the results of a two-dimensional hydro-sediment model can output the water depth across the entire computation area. In reality, the water depth is the difference between the calculated water surface elevation and the elevation in the Digital Elevation Model (DEM) below it. The greater the submerged water depth, the more significant the resulting damages.Positive
Submergence durationData exported from the simulation results of a two-dimensional water–sediment model can output the submergence duration of the flooded area. The longer the submergence duration, the greater the resulting damage.Positive
Flood flow velocityData exported from the simulation results of a two-dimensional water–sediment model can output the flood flow velocity of the inundated area. The higher the flow velocity, the more dangerous the inundated area becomes, resulting in greater losses.Positive
ResilienceNDBIA remote sensing index used to identify the distribution of buildings in urban areas, capable of identifying the distribution of buildings.Negative
Table 2. The roughness coefficients for different land use types.
Table 2. The roughness coefficients for different land use types.
TypeRange
Farm0.02–0.06
Forest0.03–0.2
Grassland0.02–0.05
Water0.02–0.035
Floodplain0.02–0.038
Building0.025–0.07
Unuse0.02–0.06
Table 3. Criteria for dividing transverse gradient levels.
Table 3. Criteria for dividing transverse gradient levels.
Low TGMedium TGHigh TG
RelationshipTG   12 LG12LG < TG 17LG17LG < TG
Note: TG represents the transverse gradient, while LG indicates the longitudinal gradient.
Table 4. Flood simulation scenario settings.
Table 4. Flood simulation scenario settings.
NumberType of FloodTGExplanation
1Flood of 1996TG of 2000Simulating a real flood to validate the model’s feasibility for calculating the inundation of floodplains.
2Low TGSetting the TG to be 10 times the LG of the river channel, simulating the flooding and free evolution process toward a two-dimensional plane. This helps determine the driving effect of low TG on the flood evolution process.
3Medium TGSetting the TG of the floodplain to be 14 times the LG of the river channel, simulating the flooding and free evolution process toward a two-dimensional plane. This helps determine the driving effect of medium TG on the flood evolution process.
4High TGSetting the TG of the floodplain to be 18 times the LG of the river channel, simulating the flooding and free evolution process toward a two-dimensional plane. This helps determine the driving effect of high TG on the flood evolution process.
Table 5. Combined weights of each evaluation indicator.
Table 5. Combined weights of each evaluation indicator.
Criteria
Source
Indicator
Waterbody
Subjective Weights
(%)
Objective Weights
(%)
Weight CoefficientsCombined Weights
(%)
α1α2
PathDistance from channel11.5429.6370.4280.57210.453
Elevation5.1576.9336.173
AcceptorRoughness2.5585.2134.077
NDVI2.2186.7274.797
ConsequenceNDWI1.3525.3413.634
Population density9.943.0345.990
Imperviousness12.4661.5526.223
GDP per unit3.32810.1847.250
Nighttime light index3.872.7823.247
Floodwater depth3.3052.2272.688
ResilienceSubmergence duration13.25410.73411.813
Flood flow velocity13.25414.83614.159
NDBI13.2546.5209.402
CriteriaIndicator4.50314.28010.095
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zheng, Z.; Li, M.; Quan, L.; Ai, G.; Niu, C.; Hu, C. Research on the Threshold of the Transverse Gradient of the Floodplain in the Lower Yellow River Based on a Flood Risk Assessment Model. Water 2024, 16, 2533. https://doi.org/10.3390/w16172533

AMA Style

Zheng Z, Li M, Quan L, Ai G, Niu C, Hu C. Research on the Threshold of the Transverse Gradient of the Floodplain in the Lower Yellow River Based on a Flood Risk Assessment Model. Water. 2024; 16(17):2533. https://doi.org/10.3390/w16172533

Chicago/Turabian Style

Zheng, Zhao, Ming Li, Liyu Quan, Guangzhang Ai, Chaojie Niu, and Caihong Hu. 2024. "Research on the Threshold of the Transverse Gradient of the Floodplain in the Lower Yellow River Based on a Flood Risk Assessment Model" Water 16, no. 17: 2533. https://doi.org/10.3390/w16172533

APA Style

Zheng, Z., Li, M., Quan, L., Ai, G., Niu, C., & Hu, C. (2024). Research on the Threshold of the Transverse Gradient of the Floodplain in the Lower Yellow River Based on a Flood Risk Assessment Model. Water, 16(17), 2533. https://doi.org/10.3390/w16172533

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop