water
Article Understanding Morphodynamic Changes of a Tidal River Confluence through Field Measurements and Numerical Modeling
Qiancheng Xie 1,* , James Yang 2,3 , Staffan Lundström 1 and Wenhong Dai 4
1 Division of Fluid and Experimental Mechanics, Luleå University of Technology (LTU), 97187 Luleå, Sweden; [email protected] 2 Division of Resources, Energy and Infrastructure, Royal Institute of Technology (KTH), 10044 Stockholm, Sweden; [email protected] 3 Vattenfall AB, Research and Development (R & D), Älvkarleby Laboratory, 81426 Älvkarleby, Sweden 4 College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing 210098, China; [email protected] * Correspondence: [email protected]; Tel.: +46-72-2870-381
Received: 20 August 2018; Accepted: 8 October 2018; Published: 11 October 2018
Abstract: A confluence is a natural component in river and channel networks. This study deals, through field and numerical studies, with alluvial behaviors of a confluence affected by both river run-off and strong tides. Field measurements were conducted along the rivers including the confluence. Field data show that the changes in flow velocity and sediment concentration are not always in phase with each other. The concentration shows a general trend of decrease from the river mouth to the confluence. For a given location, the tides affect both the sediment concentration and transport. A two-dimensional hydrodynamic model of suspended load was set up to illustrate the combined effects of run-off and tidal flows. Modeled cases included the flood and ebb tides in a wet season. Typical features examined included tidal flow fields, bed shear stress, and scour evolution in the confluence. The confluence migration pattern of scour is dependent on the interaction between the river currents and tidal flows. The flood tides are attributable to the suspended load deposition in the confluence, while the ebb tides in combination with run-offs lead to erosion. The flood tides play a dominant role in the morphodynamic changes of the confluence.
Keywords: tidal river confluence; flow features; morphological changes; field measurements; numerical simulations
1. Background A river confluence is a key feature of a drainage basin in terms of hydrology and geomorphology, for geological records, as well as from a habitat point of view [1]. In the confluence, two merging run-off streams often result in enhanced turbulent mixing. This has a bearing on transported sediment and its amount delivered downstream. In a long-term perspective, the flow patterns govern the morphology changes of the confluence, e.g., the scouring and sediment deposition [2]. Many studies focused on hydrodynamic patterns and morphology changes in run-off confluences. Mosley [3] was a pioneer in the research, which was further developed by Best who defined six distinct hydraulic regions in a confluence [2,4]. Those included areas of flow stagnation, flow deflection, flow separation, maximum velocity, gradual flow recovery, and shear layers. With advanced instrumentation and novel experimental design, research of river confluences evolved and focused on the separation zone and the shear layer [5–8]. Yuan et al. [9] made a review of the state of the art in hydraulic research of run-off confluences. Best [2] examined principal morphological features
Water 2018, 10, 1424; doi:10.3390/w10101424 www.mdpi.com/journal/water Water 2018, 10, 1424 2 of 21 such as deep scour holes, bank-attached lateral bars, tributary-mouth bars, and a region of sediment accumulation in a river confluence. Other similar studies looked at the sediment and morphological aspects at non-tidal alluvial confluences [10–12]. They contributed to the understanding of the confluence scouring. A literature review shows that limited attention is drawn to tidal confluences with bi-directional flows [4]. In tidal environments, the confluence is also affected by tidal currents. As a result, the alluvial process in terms of erosion and deposition is different, which is an issue of concern for many practical applications, especially if the confluence is in an urban development area. Pittaluga et al. [13] investigated the morphodynamic equilibrium of alluvial estuaries, where river flows and tides meet. The complexity in a tidal confluence does not draw much attention [4,14,15]. In bi-directional flows, the shift of the dominant processes between river run-off and tides, featuring periodical changes in both magnitude and direction, induces more degrees of complexity in terms of flow patterns, sediment transport, and bed morphology change in the confluence. Ferrarin et al. [16] identified 29 scour holes at tidal channel confluences by examining their geomorphological characteristics and comparing them with scours in rivers. As a consequence of changes in the flow regime, their findings revealed, on a century scale, the morphological dynamics of scouring. Field studies [10–12] and laboratory experiments [17–20] are major tools for studying sediment transport and the hydromorphic process in a confluence. In some studies, field measurements were made with tides; their purpose was to examine transport of the bed load [21–23]. The study of suspended load in tidal confluences is limited. One reason is attributable to the fact that it is not easy to, in a controlled manner in the laboratory, produce flow conditions of suspended load in combination with tides. Due to the complexity of the process, erosion and deposition are still not well understood [18]. Hypotheses are usually made that the origin of mid-stream scour is related to high flow velocity, strong turbulence, and the effect of shear layer or curvature-induced helical circulation [2,18–20]. Numerical modeling allows duplications of complex boundary conditions and prediction of different scenarios in a short time [24–29], which is also true for the study of tidal river confluences. Previous attempts were made with three-dimensional (3D) models for simulations of secondary circulations and flow variations in the vertical direction [30,31]. For a shallow confluence with a large width-to-depth ratio, a depth-averaged model is an acceptable compromise if necessary corrections of cross-circulatory motion are made [32–35]. In this research, field measurements of flow and sediment were made in the study area including the confluence. A two-dimensional (2D) morpho-dynamic model was set up, in which the cohesive suspended sediment transport was simulated. The objective was, by means of field measurements and simulations, to provide insight into the physical phenomenon that governs flow features of the tidal confluence, to describe circulatory patterns of suspended load transport, and to predict the scour-hole evolution. This study reveals the relationship between the velocity and suspended sediment movement influenced by both the run-off and the tides. The results provide reference to behaviors of tidal currents, sediment transport patterns, and fluvial process in similar situations. The paper includes a description of the study area, field measurements of flow and sediment and data analyses, numerical formulation, model set-up with calibration and validation, major flow features, and morphological changes of the confluence.
2. Study Area The confluence in question is called Sanjiangkou, formed by the Fenghua and the Yao River in southeast China. The stream after the confluence is called the Yong River. It is approximately 26 km upstream from the mouth into the Pacific Ocean (Figure1). The Yao River runs into the junction roughly at a 90◦ angle with the other two rivers. The confluence is significantly affected by combined actions of river run-off and tidal currents. The bed load is negligibly small; the sediment transport in the rivers is mainly in form of suspended load, a common feature of many fluvial rivers [36]. WaterWater 20182018,, 1010,, 1424x FOR PEER REVIEW 33 of of 2120
Measured at the normal water surface, the width of the Fenghua River ranges from 90 to 180 m, and itsMeasured average at bed the normalslope is water 0.81%. surface, The reach the width near ofthe the confluence Fenghua River is almost ranges straight; from 90 its to 180cross-section m, and its averageis nearly bed U-shaped. slope is The 0.81%. yearly The averaged reach near run-off the confluence is 53.6 m3 is/s; almost the annual straight; sediment its cross-section discharge isis nearly4.35 × U-shaped.104 tons [36]. The yearly averaged run-off is 53.6 m3/s; the annual sediment discharge is 4.35 × 104 tons [36]. The normal width of the Yao River ranges from 180180 to 230 m. About 500 m m before before the the confluence, confluence, a local constriction exists, with it itss water-surface width expanding from 140 to 180 m at at the the confluence. confluence. The dailydaily run-off,run-off, as wellwell asas sedimentsediment transport,transport, is controlled by the sluicesluice gates located 3.3 km upstream ofof thethe confluence;confluence; thethe annualannual sedimentsediment dischargedischarge isis comparativelycomparatively small.small.
Figure 1. The water system, locations of the river confluence, confluence, and hydrological stations of water levels (WL)(WL) andand ofof flowflow andand sedimentsediment (FS).(FS).
The Yong RiverRiver runsruns roughlyroughly inin thethe west–east direction,direction, withwith its normal width ranging from 150 3 toto 250250 m.m. TheThe averageaverage river-bedriver-bed slopeslope isis 0.117%.0.117%. ItsIts annualannual meanmean run-offrun-off is is approximately approximately 92 92 m m/s3/s × 9 3 (annual(annual averageaverage run-offrun-off 2.9122.912 × 109 m3).). The The peak peak discharge discharge of of run-off run-off and tides occurs normallynormally 3 during the secondsecond halfhalf ofof June,June, amountingamounting toto aboutabout 18001800 mm3/s./s. The The field field recording stations for river waterwater levels (WL) (WL) and and of of flow flow and and sediment sediment (FS) (FS) are are also also marked marked in Figure in Figure 1. At1. the At theriver river mouth, mouth, e.g., e.g.,at station at station WL6, WL6, the theannual annual mean mean tidal tidal range range is 1. is91 1.91 m; m; the the flood flood and and ebb ebb durations durations are are almost the same, approximatelyapproximately 6 6 h. h. The The tidal tidal asymmetry asymmetry aggregates aggregates from from the the river river mouth mouth to theto the upstream. upstream. In the In confluence,the confluence, the ebbthe durationebb duration is 40 minis 40 longer min longer than the than flood the duration. flood duration. At station At WL7, station the WL7, duration the differenceduration difference becomes 60becomes min [36 60]. min On the[36]. Yao On River,the Yao the River, sluice the gates sluice stop gates the stop tidal the propagation tidal propagation further upstream.further upstream. As for theAs for location the location of the upstreamof the upst limitream of limit current of current reversals reversals along thealong Fenghua the Fenghua River, theRiver, tides the affect tides approximatelyaffect approximately 15 km 15 upstream km upstream of WL7/FS7; of WL7/FS7; no records no records are availableare available to show to show the run-offthe run-off influence. influence. The The sediment sediment in the in riverthe river is mainly is ma frominly from the coastal the coastal area andarea is and carried is carried by the by tides. the tides.Local scouring is a typical morphological feature of the river confluence. Our concern of the bed morphologyLocal scouring is is its a scour-holetypical morphological evolution (Figure feature1 ).of the Previous river confluence. field measurements Our concern show of the that bed a scourmorphology hole, like is its a narrowscour-hole and evolution deep crater (Figure along 1) the. Previous Fenghua–Yong field measurements River, exists show in the that confluence a scour andhole, its like deepest a narrow point and is moredeep thancrater 10 along m below the theFenghua–Yong confluence bed River, elevation. exists in During the confluence the earlier and years, its scouringdeepest point dominated is more the than river 10 sedimentation m below the process. confluence Since bed the 1980s,elevation. the studyDuring area the was earlier affected years, by ascouring number dominated of factors, the such river as the sedimentation construction process. of the sluice Since gatesthe 1980s, on the the Yao study River area and was other affected human by activities.a number Theseof factors, changes such modified as the construction the hydraulic of the conditions sluice gates and affectedon the Yao the River erosion and potential other human in the wateractivities. system, These including changes modified the confluence. the hydraulic Bathymetric conditions surveys, and althoughaffected the irregular erosion and potential fragmental, in the show that the morphology tends to shift from erosion to deposition.
Water 2018, 10, x FOR PEER REVIEW 4 of 20
Water 2018, 10, 1424 4 of 21 water system, including the confluence. Bathymetric surveys, although irregular and fragmental, show that the morphology tends to shift from erosion to deposition. 3. Field Measurements 3. Field Measurements 3.1. Data Collection 3.1. Data Collection To map the river and confluence topography and to record the tidal hydrological data, field surveysTo map were the carried river out and during confluence two major topography periods, i.e.,and Juneto record 2015 andthe Januarytidal hydrological 2016. The formerdata, field was surveysused for were this study. carried The out river during bathymetry two major used periods, in the simulationsi.e., June 2015 was and mapped January from 2016. June The 2015–January former was used2016, for which this wasstudy. achieved The river using bathymetry an HY1600 used bathymetric in the simulations profiler (SunNavwas mapped Technology from June Co., 2015 Ltd.,‒ Tianjin,January China).2016, which The hydrologicalwas achieved data using included an HY1600 water bathymetric level, flow velocity,profiler (SunNav flow discharge, Technology sediment Co., Ltd.,concentration, Tianjin, China). grain-size The distribution, hydrological water data quality, included and water salinity. level, flow velocity, flow discharge, sedimentThe waterconcentration, levels were grain-size monitored distribution, at seven cross-sectionswater quality, (WL1–WL7),and salinity. five of which were along the YongThe River.water Tolevels measure were flowmonitored velocity at andseven suspended cross-sections sediment, (WL1 seven‒WL7), corresponding five of which cross-sections were along (FS1–FS7)the Yong River. were arranged, To measure each flow with velocity three plumb and suspended lines (Figure sediment,2). Their seven distances corresponding to the confluence cross- (measuredsections (FS1 along‒FS7) the were river arranged, centerline) each are with given three in Table plumb1. FS5 lines and (Figure FS6 are 2). a fewTheir meters distances apart to from the confluenceeach other and(measured are treated along as the sameriver centerline) section. Along are given each line, in Table the sampling 1. FS5 and was FS6 made are a at few six depthsmeters apartfrom thefrom water each surface,other and i.e., areh =treated 0, 0.2H as0, the 0.4 Hsame0, 0.6 section.H0, 0.8H Along0, and each 1.0H line,0, where the samplingH0 (m) is was the watermade atdepth six depths at each from line. Allthe thewater data surface, were recordedi.e., h = 0, in 0.2 one-hourH0, 0.4H intervals.0, 0.6H0, 0.8H0, and 1.0H0, where H0 (m) is the water depth at each line. All the data were recorded in one-hour intervals. Table 1. Distance of field measurement stations to the confluence. Table 1. Distance of field measurement stations to the confluence. Water Level Station WL1 WL2 WL3 WL4 WL5 WL6 WL7 DistanceWater to confluenceLevel Station (km) 2.20 WL1 0.25WL2 6.00WL3 14.80WL4 20.30WL5 25.20WL6 WL7 2.80 FlowDistance and Sedimentto confluence Station (km) FS1 2.20 FS2 0.25 FS3 6.00 FS4 14.80 20.30 FS5 (FS6) 25.20 2.80FS7 Flow and Sediment Station FS1 FS2 FS3 FS4 FS5 (FS6) FS7 Distance to confluence (km) 2.20 0.25 9.70 18.10 24.70 3.20 Distance to confluence (km) 2.20 0.25 9.70 18.10 24.70 3.20
With a four-beamfour-beam 600/1200-kHz600/1200-kHz RDI RDI Workhorse Workhorse acoustic acoustic Doppler Doppler current current profilers profilers (ADCPs) (Nortek group, Rud, Norway),Norway), water-flowwater-flow velocitiesvelocities andand dischargesdischarges werewere measured.measured. They were attached to the measurement vessels that werewere anchoredanchored onon land.land. The The uncertainty uncertainty of of the flow flow measurements was below ±±5%.5%. The The YJD-1-type YJD-1-type pres pressuresure sensors sensors (Tekscan, (Tekscan, Inc., Inc., South South Boston, Boston, MA, MA, USA) were used for thethe water-levelwater-level measurementsmeasurements andand theirtheir accuracyaccuracy isis ±±11 cm. Major efforts were made to sample the suspended sediment, using point point-integrative-integrative water water samplers samplers (Hoskin (Hoskin Scientific, Scientific, Ltd., Saint-Laurent, QC, Canada). Samples Samples of of the the bed load, although limited in amount, were also taken with Shipek grab samplers samplers (Envco, (Envco, Auckland Auckland,, New New Zealand) Zealand) and and their their amount amount was was calculated. calculated.
Figure 2. Sketch of measurement arra arrangementngement at a cross-section.
The grain-size distribution of the suspended lo loadad was analyzed using an automatic sieving device (SFY-D)(SFY-D) (Zhonghu (Zhonghu Scientific Scientific Ltd., Ltd., Nanjing, Nanjing, China) China) and an and automated an automated laser particle-size laser particle-size analyzer (Mastersizer2000)analyzer (Mastersizer2000) (Malvern (Malvern Panalytical Panalytical Ltd., Worcestershire, Ltd., Worcestershire, UK). Particle UK). sizes Particle falling sizes in thefalling range in thebetween range 0.0002between and 0.0002 2 mm and were 2 mm identified. were identified. The obtained The obtained data are data well aresuited well suited for calibration for calibration and validationand validation of numerical of numerical models. models. The fieldThe datafield acquireddata acquired during during the wet the season, wet season, i.e., the i.e., second the second half of Junehalf of 2015, June were 2015, analyzed were analyzed to determine to determine the sediment the sediment features in features the study in areathe study including area the including confluence. the Theyconfluence. are also They used are for also calibration used for and calibration validation. and validation.
Water 2018, 10, 1424 5 of 21
3.2. Features of Suspended Sediment For each WL and FS station, the time series of the raw measurement data were analyzed. For each time period, the average value was first obtained for each plumb line with the six points. Based on the results of the three lines, the cross-sectionally averaged value was achieved using the weighted average method. According to the grain-size distribution, the sediment is classified as sand (0.05–2 mm), silt (0.005–0.05 mm), and clay (<0.005 mm) [37]. The field measurements show that approximately 95% of the river sediment is suspended load of cohesive silt and clay, most of which is carried into the river by the tides from the coastal area, as shown later. Table2 shows their median grain sizes ( D50) from the field data. In the spring tide, the D50 values vary from 0.005–0.009 mm for the suspended load, and from 0.008–0.017 mm for the bed load. In the neap tide, the corresponding ranges are 0.006–0.009 mm and 0.009–0.020 mm, implying that the D50 values are slightly larger.
Table 2. Median grain sizes (D50) of suspended and bed load during the wet season (second half of June 2015).
Spring Tide, D50 (mm) Neap Tide, D50 (mm) Station Suspended Load Bed Load Suspended Load Bed Load FS1 0.005 0.008 0.008 0.009 FS2 0.007 0.013 0.007 0.016 FS3 0.007 0.011 0.008 0.012 FS4 0.008 0.019 0.006 0.020 FS5 (FS6) 0.008 0.011 0.007 0.011 FS7 0.009 0.017 0.009 0.013
The study area experiences a semi-diurnal tide—two nearly equal high and low tides each day, belonging to the category of incomplete standing waves. The second half of June 2015 includes 30 semi-diurnal tides, with 30 flood and ebb tides. S (kg/m3) denotes the mean sediment concentration at a cross-section. To look at its spatial changes during the period, Table3 summarizes the time-averaged S values at FS1–FS7. The following features were observed: (1) Going upstream from FS4, including the confluence, the suspended sediment experiences a trend of decrease in concentration during the semi-diurnal tides. The maximum value of S along the Yong River always occurs at FS4, with a peak value of 2.341 kg/m3 during the spring tides. (2) At the same location during the semi-diurnal tides, the sediment concentration during the spring tide is 3–8 times higher than that during the neap tide. (3) For either the spring or neap tide at the same location, the sediment concentration of the flood tides differs from that of the ebb tides; the spring tides feature higher values than the neap tides. This implies that the spring tides govern the transport of the suspended sediment from the coastal area. (4) At FS2, the flood tides carry more sediment than the neap tides, which means that the former dominates the sediment to the confluence.
Table 3. Time-averaged mean sediment concentration (S) values during the wet season (second half of June 2015).
Spring Tide (kg/m3) Neap Tide (kg/m3) Station Flood Tide Ebb Tide Semi-Diurnal Tide Flood Tide Ebb Tide Semi-Diurnal Tide FS1 - - 0.229 - - 0.069 FS2 0.780 0.541 0.677 0.121 0.089 0.106 FS3 1.247 1.615 1.441 0.200 0.215 0.208 FS4 2.341 2.117 2.207 0.232 0.318 0.282 FS5 (FS6) 1.673 1.476 1.578 0.314 0.214 0.265 FS7 0.895 0.920 0.909 0.142 0.141 0.141 Water 2018, 10, x FOR PEER REVIEW 6 of 20
Water 2018, 10, x FOR PEER REVIEW 6 of 20 Water 2018, 10, 1424FS7 0.895 0.920 0.909 0.142 0.141 0.141 6 of 21 FS7 0.895 0.920 0.909 0.142 0.141 0.141 To further unveil the streamwise influence of the tides on the sediment, cross-sections FS5 (FS6) andTo FS2To further further were selected unveil unveil the theon streamwisethestreamwise Yong River. influence influence The former of of the the istides tides close on on to the the sediment, sediment,river mouth cross-sections cross-sections and the latter FS5 FS5 to(FS6) the (FS6) andandconfluence. FS2 FS2 were were Figures selected selected 3 and onon the the4 compare, YongYong River. for the The spring former former and is is closeneap close totides, to the the theirriver river relationshipmouth mouth and and the between the latter latter toS andthe tothe confluence.confluence.V, where V Figures (m/s)Figures refers3 and3 and to4 compare,the4 compare, mean flow for for the velocitythe spring spring at and aand river neap neap cross-section, tides, tides, their their relationship positiverelationship toward betweenbetween the sea SS [38].andand V, whereVA, negativewhereV (m/s) V (m/s)V value refers refers implies to to the the tidal mean mean reversal flow flow velocity velocity of the current. at a a ri riverver cross-section, cross-section, positive positive toward toward the thesea sea[38]. [ 38]. AA negative negativeV Vvalue value implies implies tidaltidal reversal of the the current. current.
Figure 3. Station FS5 (FS6), mean flow velocity vs. mean sediment concentration (V–S) relationship for spring and neap tides: (a) Spring tide; (b) Neap tide. FigureFigure 3. 3.Station Station FS5 FS5 (FS6), (FS6), mean mean flow flow velocity velocity vs. vs. mean mean sediment sedimentconcentration concentration ((VV––SS)) relationshiprelationship for springfor spring and neapand neap tides: tides: (a) Spring(a) Spring tide; tide; (b) ( Neapb) Neap tide. tide.
FigureFigure 4. 4.Station Station FS2, FS2,V V––SS relationshiprelationship for for spring spring and and neap neap tides: tides: (a ()a Spring) Spring tide; tide; (b) ( bNeap) Neap tide. tide. Figure 4. Station FS2, V–S relationship for spring and neap tides: (a) Spring tide; (b) Neap tide. TheTheV Vand andS Svalues values inin thethe springspring tides are are larger larger than than in in the the neap neap tides. tides. S variesS varies significantly significantly duringduringThe a a tidal tidalV and period, period, S values with with in twotwo thepeaks, peaks,spring implyingtides are thlarger thatat the the than sediment sediment in the concentration neap concentration tides. S isvaries dominated is dominated significantly by bythe the tidesduringtides along along a tidal the the period, river. river. The Thewith floodflood two peaks, tidaltidal V implying and SS areare th alwaysat always the sediment out out of of phase phaseconcentration with with one one isanother; dominated another; the the peakby peakthe S S appearstidesappears along during during the theriver. the floodflood The tides.floodtides. tidalAlong V andthe the river,S river, are alwaysthe the peak peak out of of ofSS phasealmostalmost with synchronizes synchronizes one another; with with the that peak that of VS of V exceptappearsexcept for for during during during the the the flood neap neap tides. tidetide atatAlong thethe riverriverthe river, mouth. the This Thispeak implies impliesof S almost that that the synchronizes the sediment sediment concentration with concentration that of Vis is subjectedexceptsubjected for to to during modifications modifications the neap byby tide thethe at tidestides the at river the the mouth.rive riverr mouth mouth This andimplies and other other that oceanic oceanic the sediment processes processes concentration like like coastal coastal up- is up- andsubjectedand downwelling, downwelling, to modifications storm storm surges, surges, by the etc. etc.tides at the river mouth and other oceanic processes like coastal up- andDuring downwelling,During the the spring spring storm tide, tide, surges, it isit the is etc. strongthe strong tidal tidal currents currents that dominatethat dominate the river the flowriver andflow that and transport that thetransport sedimentDuring the towardthe sediment spring the toward confluence.tide, it theis theconfluence. During strong the tidal Duri neap currentsng tide, the neap the that tide tide, dominate is the comparatively tide the is comparativelyriver weak;flow and its weak;S thatvalues aretransportits much S values lower the are sediment (Figuresmuch lower toward3b and (Figures4 theb). confluence. Judging 3b and from4b). Duri Judging this,ng onethe from neap can this, saytide, one that the can thetide say run-off is thatcomparatively the plays run-off a major weak; plays role. Thisitsa major S leads values role. to are the This much non-similarity leads lower to the (Figures non-similarity of S and3b andV between 4b). of S Judging and the V between springfrom this,and the one spring neap can tides. sayand that neap The the tides. tidal run-off currentsThe plays tidal are currents are essential for stirring sediment, modifying its peak duration, as well as its transport. If essentiala major forrole. stirring This leads sediment, to the non-similarity modifying its of peak S and duration, V between as the well spring as its and transport. neap tides. If there The tidal was no there was no tide, S would be directly proportional to the river discharge and sediment would only tide,currentsS would are essential be directly for proportionalstirring sediment, to the modifying river discharge its peak and duration, sediment as well would as its only transport. be diverted If be diverted downstream [39]. downstreamthere was no [39 tide,]. S would be directly proportional to the river discharge and sediment would only be diverted downstream [39]. 4. Numerical Modeling 4. Numerical Modeling 4. NumericalBased on Modeling the Delft3D package [40], a 2D depth-averaged model was used to help understand Based on the Delft3D package [40], a 2D depth-averaged model was used to help understand the complex flow features and morphology changes of the confluence. Delft3D-Flow is a separate the complexBased on flow the features Delft3D andpackage morphology [40], a 2D changes depth-averaged of the confluence. model was Delft3D-Flowused to help understand is a separate module in the package, in which sediment motion and morphology change are coupled with the flow modulethe complex in the package,flow features in which and morphology sediment motion changes and of morphology the confluence. change Delft3D-Flow are coupled is witha separate the flow moduleto simulate in the flow package, patterns in andwhich morphological sediment motion changes. and morphology change are coupled with the flow to simulate flow patterns and morphological changes. to simulate flow patterns and morphological changes. 4.1. Mathematical Formulation The governing equations are based on the Navier–Stokes equations, with Leibniz integration in the vertical direction, to obtain depth-averaged flow parameters. The vertical flow acceleration Water 2018, 10, 1424 7 of 21 is neglected, leading to the hydrostatic pressure. The turbulence shear stress is solved by the k-ε turbulence model. The river-bed stability coefficient and resistance are parameters that govern the sediment scouring and deposition. The governing equations include, therefore, mass continuity, flow motion, sediment transport, and bed deformation. The depth-averaged continuity equation is ∂ξ 1 ∂ HUξ Gηη 1 ∂ HUη Gξξ + + = q, (1) ∂t Gξξ Gηη ∂ξ Gξξ Gηη ∂η where H is water depth, Uξ(ξ, η) and Uη(ξ, η) (m/s) are the depth-averaged velocities in the ξ and η coordinate system, Gξξ and Gηη (m) are coefficients for transformation of the curvilinear to orthogonal coordinates, q (m/s) is the global source or sink term per unit area, and t (s) is time. The momentum equations in the ξ and η directions are
2 ∂Uξ Uξ ∂Uξ Uη ∂Uξ Uη ∂Gηη UξUη ∂Gξξ 1 + + − + − f Uη = − Pξ + Fξ + Mξ, (2) ∂t Gξξ ∂ξ Gηη ∂η Gξξ Gηη ∂ξ Gξξ Gηη ∂η ρ0Gξξ
2 ∂Uη Uξ ∂Uη Uη ∂Uη Uξ ∂Gξξ UξUη ∂Gηη 1 + + − + + f Uξ = − Pη + Fη + Mη, (3) ∂t Gξξ ∂ξ Gηη ∂η Gξξ Gηη ∂η Gξξ Gηη ∂ξ ρ0Gηη 3 where ρ0 (kg/m ) is water density (except in the baroclinic pressure terms, variations in ρ0 are 2 2 neglected), Pξ and Pη (kg/(m s )) are pressure gradients, f (1/s) is the Coriolis parameter (i.e., inertial 2 frequency), Mξ and Mη (m/s ) refer to contributions from external sources or sinks of momentum 2 (within the computational domain Mξ = Mη = 0), Fξ and Fη (m/s ) are horizontal Reynolds stresses ∂τ ∂τ based on the eddy viscosity concept and change in both space and time, and F = 1 ξξ + 1 ξη , ξ Gξξ ∂ξ Gηη ∂η 1 ∂τηξ 1 ∂τηη 2 Fη = + , and τ ,τηη, τ , and τ (N/m ) are turbulence shear stresses. Gξξ ∂ξ Gηη ∂η ξξ ξη ηξ As noted from the in situ data, the bed load amount is small (below 5%). For the sake of simplification, only the cohesive suspended sediment is considered in the sediment model. The 2D transport equations for suspended load are given by