A Simplified Concept to Test Erosion or Sedimentation Potential along a Tidal River in the Ganges-Brahmaputra-Meghna Delta ()
1. Introduction
Estuarine environment of deltaic river systems is always attractive for human settlements and economic activities because of its valuable ecosystems, natural resources and cost-effective communications, although, there exists a number of inherent threats and challenges experienced by multiple natural hazards such as flooding, erosion, storm surges and salinity coupled with anthropogenic stresses (Ericson et al., 2006; Syvitski et al., 2009). The proximity of the river banks is the most suitable locations for such developments and ensuring the safety against the hazards, specially, the land erosion is one of the important challenges (Akter et al., 2019; Dastagir, 2014; Rahman et al., 2020) for the sustained development activities. Moreover, with the growth of population and economic expansion land reclamation in these environments are widely chosen using a wide range of methods including cross-dams to ensure the emerging required spaces (Rahman et al., 2022). The elevation of the lower region of the Ganges-Brahmaputra-Meghna (GBM) delta, has an average height of 1 - 2 m above mean sea level and large areas of land mass are exposed to natural inundation due to higher tidal amplitude in the tide dominated delta (tidal amplitude ranges from 2 - 7 m from the West to East) that often propagates up to 100 km towards landward direction (Hussain & Tajima, 2017; Nicholls et al., 2018). The estuarine networks in the southwest region of the GBM delta are recognized as an eastern, central, and western estuarine system (henceforth EES, CES, and WES, respectively), connected through several cross channels (Rahman et al., 2022). There exist strong seasonal and spatial variations of water-sediment regime that causes regional and temporal changes of erosion-sedimentation within the estuarine system (Dasgupta et al., 2014; Haque et al., 2016).
Generally, erosion is occurring at locations where transport capacity exceeds the availability of the sediment flux, while sedimentation occurring under the reverse condition when sediment load exceeds its transport capacity. The former case we define here as the condition of “sediment deficit system,
SDS
” and the latter case, we define as “sediment surplus system,
SSS
”. In semidiurnal tidal environment, the flow-sediment regime changes from flood tide to ebb tide or ebb tide to flood tide that eventually lead towards the changes from
SDS
to
SSS
or
SSS
to
SDS
, seasonally, fortnightly or even daily, at a fixed location. On the other hand, the development sites generally in need for stable land mass and sustainable navigational routes to facilitate smooth waterway communications. To meet the above requirement, it is important to identify whether the site is dominated by or sediment surplus to facilitating the selection of appropriate sites. Through such processes, the prevailing natural processes can be utilized and the development efforts, repair and maintenance will be at the possible minimum level.
As part of wider development activities (Bangladesh Planning Commission, 2018), the Government of Bangladesh (GoB) has planned to establish a large build up area on a sedimentary deposit along the Payra river, developed in the last 40 years at Lebukhali which is about 60 km upstream from the mouth of Bishkhali/Burishwar estuary (Figure 1). The natural formative processes of this land mass (sand bars) through development of bends and cut-offs of the meandering reach are discussed earlier (Islam et al., 2018a) which is further demonstrated
Figure 1. Overall study area, about 60 km upstream from the mouth of Bishkhali/Burishwar estuary. (a) Central estuarine system (CES) of the southern estuarine region of Bangladesh, (b) study site with surrounding districts and river system, (c) evolution of the study area from 1975-2016, and (d) Study site with closer look, 22.274490˚N, 90.173525˚E to 22.301957˚N, 90.214919˚E.
briefly in Figure 1(c). The site is located along the right bank of the Payra river that is tidal in nature. It is located within the CES, is dominated by fresh water, and flows between the districts of Barisal and Patuakhali which are considered as food basket of the country in terms of fisheries and agriculture. The river is linked with the Lower-Meghna through the Tentulia in the north-east, flows to the south and finally falls to the Bay of Bengal through the Bishkhali/Burishwar estuarine mouth. Around 2 km reach of 6.5 km long sedimentary deposit is taken as the study site at the downstream of the existing ferry ghat at which Lebukhai bridge is constructed recently (Figure 1(d)).
The river reaches in this location transports dominantly fresh water and sediment under the complex interactions of fresh water input from the lower Meghna system, tidal exchange, and other fluvial and marine processes. Sediment transport and associated erosion-deposition processes in the study reach is much complicated as it is characterized by the variation of tidal duration, peak velocities and discharges during daily ebb and flood cycles and fortnightly neap-spring cycles as well as slack time during flow reversal. Although the existing four major tidal creeks are kept active during the development work, the study site is still experiencing intermittent bank erosion (Islam et al., 2018b) and the authority needed to manage river erosion in order to continue the planned development activities. Within the above backdrop, the objective of this paper is to identify the reasons of erosion through a confirmative study whether the site is existing within
SDS
or
SSS
so that the authority can take appropriate steps to ensure the stability of the developing site.
2. Methodology
The study site does not experience any significant wave action except the seasonal wind induced wave and therefore, only current related suspended sediment transport is considered here. As the sedimentary deposit along the right bank of the study site experiences intermittent erosion as demonstrated by Figure 2, with greater intensity during dry season and ebb tide phase, the aim of this paper is to identify the diagnostic reason for the observed erosion through relating the changes of transport capacity at different tidal phases.
2.1. Conceptualization of Sediment Surplus or Sediment Deficit
As the flow induced bed shear velocity (
u∗ ) exceeds a critical value for sediment transport (
u∗,cr ), bed particles start moving by sliding, rolling or jumping, usually known as bed load transport. As the vertical component bed shear velocity (
u∗,cr,sus, Equation (1) & (2)) exceeds the particle fall velocity (
ws, Equation (3)), the particle moves in suspension towards the vertical direction (Bagnold, 1966). This movement continues in the direction of flow as if horizontal velocity component is strong enough to carry them away (Guo et al., 2018). The above facts of sediment transport processes are schematically conceptualized in Figure 3 within the tidal phases, namely during flood and ebb tides.
Critical bed shear velocity (
u∗,cr,sus ) for the initiation of suspension can be calculated using the relations of critical Shields parameter for suspension (
θcr,sus ) as follows (Soulsby, 1997; van Rijn, 2019):
θcr,sus=u2∗,cr,sus(s−1)gD (1)
Figure 2. Intermittent erosion-accretion along the bank line of the study site observed during dry season ebb time. Images of 2006, 2012 and 2016. Image: 22.264607˚N, 90.164659˚E to 22.293240˚N, 90.203933˚E.
Figure 3. Schematic diagram for sediment transport in a semidiurnal tidal river. Bed load transport takes place in flow direction by sliding or rolling when
u∗>u∗c,sr and by jumping when
u∗>u∗c,js. Bank and bed material goes for suspension when
u∗>u∗c,sus or
ws (Hinze, 1975). These suspended particles plus wash load of river and marine flows move horizontally when horizontal component of flow velocity is strong enough to transport them away.
and
θcr,sus=0.31+D∗+0.1[1−exp(−0.05D∗)] (2)
where, the dimensionless grain size parameter
D∗=[(s−1)gυ²]1/3D.
Particle fall velocity for marine environment can be estimated as follows:
ws=νD[(10.362+1.049D3∗)12−10.36] (3)
where, D is the mean grain size assumed as
D50.
However, for the equilibrium transportation of suspended sediment, the volume of deposition caused by gravity must be equal to the volume of upward sediment diffusion, typically conceptualize as
qv,downward=qv,upward. If the fall velocity ws of particles of uniform size having the time-averaged suspended sediment concentration (SSC)
C(z) at a vertical position z above the bed, the downward volume of flux can be expressed as
qv,downward=wsC(z), while the upward
volume flux of particles can be expressed as
qv,upward=εsdC(z)dz (Figure 4)
where
εs is the vertical diffusion coefficient (Keetels, Goeree, & van Rhee, 2018).
Equating the two fluxes gives an expression for the vertical distribution of SSC in equilibrium over the depth (Nie et al., 2017; van Rijn, 1993, 2012) as
wsC(z)=−εsdC(z)dz
Or,
wsC(z)+εsdC(z)dz=0 (4)
In Equation (4), the diffusion coefficient
εs is assumed proportional to the
kinematic eddy viscosity,
vt as,
εs=βvt with,
vt=κu∗zh(1−zh) (Tsujimoto, 2010) and the Equation can be solved as (Rouse, 1937):
Figure 4. Suspended sediment transport. Volume of suspension
qv,upward and deposition
qv,downward at bottom boundary layer (
z0=a ) for suspension will be equal i.e.
wsC(z)=εsdC(z)dz.
C(z)C(a)=[(h−zz)(ah−a)]N with
N=wsβκu∗ (5)
where, C(a) = reference concentration at reference level a above river bed, h= flow depth, z= height above the bed level. Here, N is known as suspension number or Rouse number by which the relative distribution of the suspended sediment transport is determined. κ = Karman’s constant = 0.4, and β = dimensionless proportionality coefficient related to diffusion of sediment particles.
Equation (5) is widely known as Rouse equation which is used to estimate equilibrium profiles of vertical sediment concentration in tidal and non-tidal rivers (Van Rijn, 1984; Wu et al., 2017). Quantitatively, if the actual sediment concentration exceeds the value of the estimated equilibrium concentration, the hydraulic condition is recognized here as sediment surplus, while the condition is considered as sediment deficit if the actual sediment concentration is less than the equilibrium value. However, in Equation (6), the vertical concentration profile depends significantly on the selection of appropriate reference height a and reference sediment concentration C(a). It is found in several previous researches (Garde & Ranga Raju, 1977; Parker, 1990; Singh et al., 2004) that value of a can be estimated as a function of either flow depth or particle diameter. For plain bed condition, the selection of a based on particle diameter might be applicable but for rough bed, the estimation of a can be done related to flow depth or bed roughness height. The real life field condition is usually representing rough bed which is also valid in the study site. As the measurement of roughness height is relatively difficult in field condition, the choice of appropriate value of a would be evaluated as a function of flow depth (0.01h to 0.05h) as prescribed in previous researches (Engelund & Fredsoe, 1976; Garcia & Parker, 1991; Shibayama & Winyu, 1993; Van Rijn, 1984; Zyserman & Fredsoe, 1994). To estimate another important parameter C(a) at reference level, either the numerical expression (Equation (6)) prescribed by van Rijn (1984, 2007) can be used or the value can be measured in the field.
C(a)=0.015d50T1.5ad0.3∗ (6)
where, T = transport stage parameter
=u2*−u2∗,cru2∗,cr.
2.2. Measurement of Flow and Sediment Parameters
To evaluate the parameters in Equations (2)-(6), required arrangements to collect flow and sediment data were established in the field and data collection continued over the period of 5 years from May 2016 to January 2020. Different sensors coupled with total station/RTK GPS were deployed to measure the spatial and temporal changes of flow, depth, sediment concentration and related parameters (Islam et al., 2021). ADCP is used to measure the velocity and discharge data across the river at different sections of the study site. Vertical velocity profiles are then retrieved from this ADCP records which are used for subsequent analyses. Peak ebb velocities are found varying from 0.59 to 0.79 m/s and peak flood velocities varying from 0.52 to 0.70 m/s in the upstream (VIP Ghat) to downstream location (OBM Ghat). Eco-sounder is used for river bathymetry, and Total Station/RTK GPS is used for horizontal location, digital elevation, and water level data. OBS-3A sensor was mounted on a steel frame and placed 2 m above the bed level at both locations to measure the turbidity from 7 am to 6 pm, in a specified day. The suspended solid concentration at the same locations were also measured. The measured suspended solid concentration and OBS turbidity reading were then calibrated. The calibrated equation was then used to calculate the SSC. The bathymetry and water level data were measured covering the entire study reach, while detailed of sediment concentration and flow velocities were measured at two locations, the OMB Ghat and VIP Ghat (Figure 1) at different vertical points within different tidal phases. Grain size analysis was done to estimate the sediment characteristics as well as transport parameters of the river reach.
3. Results and Discussions
3.1. Estimated Parameters
Peak flow velocities for flood and ebb tides, grain size, fluid properties were used to calculate the various hydrodynamic parameters. Summary of the estimated parameters in Equations (1)-(6) are shown in Table 1, whereas, details of measured data and estimation of required parameters are discussed in the supplementary.
Here, the reference concentration C(a) is taken at reference height, a = 0.15 to 0.165 m (which is the average of 0.01h to 0.05h), while h varies between 5 to 6 m.
Parameters | VIP Ghat | OBM Ghat | Equations |
Current related shear velocity, (m/s) | FT = 0.054 FE = 0.063 | FT = 0.074 FE = 0.085 | |
Current related shear stress (Nm−2) | FT = 2.94 FE = 4.12 | FT = 5.57 FE = 7.41 | |
Dimensionless grain size, | 2.38 ≈ 2.4 | 2.39 ≈ 2.4 | |
Critical shear velocity (m/s) | 0.0128 | 0.0128 | |
Settling velocity, ws (m/s) | 0.0073 | 0.0073 | |
Rouse number, N | FT = 0.34 FE = 0.29 | FT = 0.25 FE = 0.22 | |
Reference concentration, Ca (mg/l) at a = 0.155 m | FT 690.0 FE 406.7 | FT 560.0 FE 156.2 | Field data analyzed in laboratory |
Transport stage parameter, T | FT = 16.8 FE = 23.0 | FT = 32.4 FE = 43.0 | or |
翻译:
Table 1. Summary of different hydrodynamic parameters.
The Rouse number N, during the peak ebb and peak flood flow varies 0.25 to 0.22 at OBM Ghat and the above values are relatively higher in the
VIP
Ghat (0.29 to 0.34). While the corresponding transport stage parameter, T, during the peak flood and peak ebb flow varies 32 to 43 at OBM Ghat and the above values are relatively lower in the
VIP
Ghat (16 to 23). It is important to note that the higher value of N indicates the lower level of equilibrium sediment concentration. As N values in
VIP
Ghat are higher than OBM Ghat, the transport capacity in the
VIP
Ghat is lower as compared with the OBM Ghat. Moreover, the transport stage during the ebb tide is higher as compared with flood tide in both locations which indicates the sediment transport is higher during the ebb tide and residual transport occurred along the ebb direction.
3.2. Sensitivity of Rouse Number on Vertical Profiles of Sediment Concentration
Using C(a) is equal to 690 mg/l and a is equal to 0.15 m, and the sensitivity of N on the vertical concentration of sediment profiles are show in Figure 5.
It is seen that the larger the value of N, the more rapidly the suspended-sediment concentration decreases with height above the reference level, a. The value of N is a function of flow induced shear velocity and particle fall velocity. The particle fall velocity is primarily governed by the particle characteristic and fluid property which generally remain constant for a particular river segment, whereas, the value of flow velocity can be modified often due to anthropogenic interventions. Thereby, the value of N can be increased or decreased to meet the specific requirements of the site. If a site needs sedimentation for land reclamation, the
Figure 5. Vertical distribution of SSC for different N values. The smaller the N value, the larger the sediment concentration.
required value for N can be adjusted through suitable interventions. Reverse situation can be also ensured where lowering of riverbed is essential, for example in navigational routes or landing ports.
3.3. Testing for Sediment Deficit System or Sediment Surplus System
Using the required parameters calculated from measured data, sediment concentration profiles are plotted and compared with the equilibrium profiles which are estimated using Equation (5) for both the upstream (VIP Ghat) and downstream (OBM Ghat) locations. The equilibrium profiles are then compared with the measured concentration of those locations during full tide and full ebb. The selection of reference level a and corresponding reference concentration C(a) are important to have an accurate distribution profile. Criteria for selecting reference height are already discussed. Since measured sediment concentrations of the study site are available, the estimated reference sediment concentration obtained using Equation (6) is not used here. SSC were measured for four different time periods of a tidal cycle; half tide, full tide, half ebb and full ebb (shown in supplementary). Considering the sensitivity of reference concentration on the distribution profile, these profiles are plotted considering the specific reference concentration of different tidal periods as shown in Figure 6.
The profiles show that the differences between equilibrium concentration and measured concentration are different in full tide and full ebb situations at both upstream and downstream locations. The difference between measured concentration and equilibrium concentration is less in the upstream location but more in the downstream location. This implies downstream location exhibits higher SD and hence more liable to erosion than upstream location. Again, in the upstream or downstream locations, concentration profiles are different during tide and ebb phases. In the upstream, measured concentration nearly merges with
Figure 6. Depth averaged SSC profiles: measured vs. equilibrium. (a) Upstream VIP Ghat location during full tide (FT), (b) upstream VIP Ghat locations full ebb (FE), (c) downstream OBM Ghat location during full tide and (b) downstream OBM Ghat locations during full ebb. The difference between measured concentration and equilibrium concentration is less in the upstream location (a & b) but more in the downstream location (c & d). Again, within upstream or downstream locations, difference between measured concentration and equilibrium concentration is more during FE than FT.
the equilibrium concentration during full tide but lags behind the equilibrium concentration during full ebb. This implies the site is nearly in the equilibrium condition during full tide but remain in SD condition during full ebb and hence more prone to erosion during ebb. In the downstream, measured concentration is found significantly less than the equilibrium concentration during both full tide and full ebb but the difference is more during ebb phase than tide phase. This means the downstream location remains in SD condition and is liable for erosion throughout the tidal cycle. Moreover, the profiles show that SD is more during ebb than during flood which implies that the ebb flow has more sediment transport capacity than the tidal flow. Therefore, intensity of erosion is more during ebb phase than during tide phase. To attain the equilibrium condition, additional sediment supply is required which would come from bank or bed erosion. Therefore, the analysis reveals that the overall status of the study site is sediment deficit condition and thus the site is subject to riverbank or bed erosion. This phenomenon matches quite well with the field condition as it is noticed during field observation and data collection. The field scenario reveals that bank line of 7 - 10 m width particularly in the downstream location had been washed away in last 4 - 5 years period with greater intensity during ebb phase. The above field observations support the hypothesis that system is SD and erosion may continue in the future. Again the time series plot of bank lines in Figure 7 (reproduced from google earth), confirm that during the period of 2006-2012, the total land loss was around 0.56 km2, while during the period of 2012-2016, the lost land was around 0.33 km2 which translate a consistent erosion rate of around 0.09 km2/year. While the rate of accretion translates within the above period with a much lower value which is around 0.01 km2/year (around 1/10th). The above quantitative analysis of images, validates that the site is erosion dominating that need continuous attention to protect land.
3.4. Implications to Sediment Management
Sediment management is one of the big challenges and less explored issue in the GBM delta which is very important to materialize the inherent opportunities with efficient use of incoming sediment as resource in these environments. The practiced methods are focusing either to accelerate the siltation or de-siltation based on the requirements. For accelerating siltation, cross-dams are typically applied in the coastal zone to trap the sediment for land reclamation which has some inherent consequences to the adjacent areas to initiate or accelerate fresh erosion (Rahman et al., 2022). Tidal river management (TRM) are usually implemented to elevate low lying land inside polders with additional benefit of increased conveyance capacity of the adjacent rivers (Talchabhadel et al., 2020). Bandal-like structures are adopted to increase the conveyance capacity of the navigational routes and ports and enhance near bank sedimentation (Rahman et al.,
Figure 7. Image analysis for the validation of SDS. The site was in SDS throughout the tidal cycle. That is why, the quantity of erosion in found significant.
2020; Kibriya, 2020). Dredging, on the other hand, is adopted for the maintenance of navigational depths in rivers and recently been applied to solve waterlogging problems within polders (WARPO & BUET, 2020). As the working principles of the methods are different and are applied in isolation depending on the objectives, the interactions within each of the methods are also relatively unknown (Rahman et al., 2022). Therefore, the identification of SDS and SSS will enable to adopt site specific appropriate measures and/or its combinations for sediment management. Through the application of remote sensing and
GIS
, the real time distribution of sediment concentration is possible to evaluate that can be further utilized to identify whether site is sediment surplus or sediment deficit. Such an insight on the prevailing morpho-dynamic behaviour will be very useful in devising an appropriate solution for river channel stabilization and maintenance of navigational channel and ports.
There exists plenty of examples where river management projects are unsuccessful in Bangladesh, specially, in the estuarine environment because of the inappropriate sediment management options. The improved understanding of whether it is sediment deficit or sediment surplus, may help to avoid the further negative consequences. Further, early identification would provide warning for subsequent erosion or accretion that can enable to take timely countermeasures and thereby could minimize the losses and save the important installations.
4. Conclusion
The selection of appropriate sites in terms of physical stability in the tidal environment is very important to reduce the challenges related to required repair and maintenance against the changes of river channels. The undesirable erosion-scouring-siltation often leads to the unsustainability of the project site leading to the complete failure. Therefore, methodology adopted in the present study to evaluate sediment deficit or sediment surplus can identify the river reach whether it will be an eroding or siltation dominated site which in turn will create the ground for technical justification in the selection processes of a development project site. The study site is selected for the establishment of major infrastructures with the understanding that the sedimentary deposit will remain stable. However, the findings from this study reveal that the site, in general, is representing the SDS with residual sediment transport towards the ebb direction that will experience continuous riverbank erosion. More specifically, the downstream reach with higher depth exhibits more transport capacity than upstream and the ebb phase experiences more sediment deficit than tide phase. Consequently, the concerned river reach needs appropriate bank protection measures to ensure the physical stability of the project.
Further field validation is required to test this method in sediment deficit region in the future through which the above method can be replicable to other sites to understand whether a certain reach of a river system exists within the sediment deficit or sediment surplus and can be an important technique for the selection of stable project site. Moreover, it can be regarded as the precursor of likely erosion or sedimentation and appropriate measures can be adapted accordingly. Real time flow-sediment data would complement in developing a sound sediment management for implementable and environmentally adaptable development plan which needs to be explored in the future.
Supplementary: Measured Data and Estimation of Required Parameters
Velocity profile
Velocity data are measured using ADCP across the river at different sections. Vertical velocity profiles of two locations about 20 m away from the bank line were retrieved from ADCP records for both ebb and flood peaks. The vertical velocity profiles are shown in FigureS.1.
Suspended sediment data
Sediment sample were collected from two locations (VIP Ghat and OBM Ghat, 0.6 km away from each other) using OBS and manually using improvised sample collector. The results are shown in TableS.1 and TableS.2 below respectively for two locations 0.6 km away from each other. Samples were collected with respect to local measurement unit feet and then converted to metric unit.
Bed soil data
River bed sample were collected from both VIP Ghat and OBM Ghat locations of the study site and grain size distributions are done. The results are shown in TableS.3 and TableS.4.
Figure S.1. Vertical velocity profile. Upstream location (left), downstream location (right).
Depth (ft) | Half Tide | Full Tide | Half Ebb | Full Ebb |
19 | | 152.2 | | |
18 | | 154.5 | 32.2 | |
17.5 | 35.1 | | | |
16 | | | | 45.8 |
15.5 | 40.2 | 161 | 35 | |
12.5 | 48 | 190.2 | 50 | 46.4 |
9.5 | 67.5 | 218.9 | 70.5 | 45 |
6.5 | 86 | 256.1 | 77.9 | 51.5 |
3.5 | 109 | 288.6 | 106.1 | 73 |
0.5 | 270 | 560 | 430 | 156.2 |
bottom | | | | |
翻译:
Table S.1. Sediment concentration at VIP Ghat.
Depth (ft) | Half Tide | Full Tide | Half Ebb | Full Ebb |
21 | - | 51.3 | - | - |
20 | 32.7 | | - | - |
19 | | | 36.6 | - |
18.5 | 34 | 50 | | |
18 | | | | 35.5 |
15.5 | 33 | 57 | 39.3 | 36 |
12.5 | 38 | 64.5 | 47 | 46 |
9.5 | 52 | 81.9 | 61.5 | 38.3 |
6.5 | 118.8 | 130.5 | 72.7 | 50.1 |
3.5 | 128.8 | 157.5 | 110.4 | 95.5 |
0.5 | 670 | 690 | 560 | 406.7 |
(bottom) | | | | |
翻译:
Table S.2. Sediment concentration at OBM Ghat.
| Gravel (>4.75 mm) | Sand (4.75 - 0.075 mm) 69.0% | Silt (0.075 - 0.005 mm) | Clay (0.005 - 0.001 mm) |
Coarse sand (4.75 - 2.00 mm) | Medium sand (2.00 - 0.425 mm) | Fine sand (0.425 - 0.075 mm) |
OBM Ghat | 0.0% | 0.0% | 2.0% | 67.0% | 22.0% | 8.0% |
VIP Ghat | 1.7% | 0.0% | 1.3% | 70% | 25% | 1.0% |
翻译:
Table S.3. Grain size limit classifications (ASTM D 7928-17), OBM Ghat and VIP Ghat.
| D95 | D84 | D60 | D50 | D30 | D16 | D10 | D | |
OBM Ghat | 0.3000 | 0.1800 | 0.1234 | 0.1042 | 0.0767 | 0.0300 | 0.0062 | 0.1047 | 0.082 |
VIP Ghat | 0.3000 | 0.2000 | 0.1208 | 0.1049 | 0.0790 | 0.0400 | 0.0214 | 0.1049 | 0.073 |
翻译:
Table S.4. Grain size parameters (mm).
Sediments become cohesive when the clay and silt contents are over 3% - 5% (Van Ledden, Van Kesteren, & Winterwerp, 2004). Here the clay particles are found 8%. Again, the sample may be termed as weakly cohesive if it contains 60% - 80% of sand (>62 µm) (van Rijn et al., 2007). Moreover, sediment will behave as cohesive when mud fraction (<0.05 mm) are 30% (van Rijn, 2012). Result shows that particle size < 0.05 mm for both locations are below 30%. Therefore, the sediment of the study site is weakly cohesive. Sorting of both samples are less than 0.350 indicating very well sorted and have undergone much transport (Folk, 1974).
Flow induced shear stress (τb) and shear velocity (u*)
The current related bed shear velocity (
u∗ ) and bed shear stress (τb) can be obtained from following equations:
u(z)=u∗kln(zz0) (A.1)
and
τb=ρwu2∗ (A.2)
where z = vertical coordinate representing water depth, k = Von Karman’s constant, usually taken as 0.4, and z0 = vertical level with zero velocity, also often referred to as bed roughness. For hydraulically rough flow z0 = 0.033ks where, ks is bed roughness (Soulsby, 1997). During field measurement velocity was observed zero nearly at 0.1 m depth. Therefore, z0 is reasonably taken as 0.1 m.
For OBM Ghat: for ebb condition u at z= 4 m is 0.785 m/s, z0 = 0.1 m. Thus
using Equations (A.1) & (A.2),
u∗=0.085 m/s. [
0.785=u∗0.4ln(40.1) ] and corresponding
τb=1023×(0.078)2=7.41 N⋅m−2.
Similarly for flood condition u at z = 4.5 m is 0.702 m/s. Thus
u∗=0.074 m/s and
τb=5.57 N⋅m−2.
VIP Ghat: for ebb condition u at z = 4 m is 0.585 m/s, z0 = 0.1 m. Thus
u∗=0.063 m/s.
τb=1023×(0.063)2=4.12 N⋅m−2.
For flood condition u at z = 5 m is 0.524 m/s. Thus
u∗=0.054 m/s. and
τb=1023×(0.054)2=2.94 N⋅m−2.
Critical shear velocity (u*,cr,sus) and shear stress (τ*,cr,suspension) for suspension
For 20˚C and for
s=2.61,
ν=1.0×10−6 m2⋅s−1.
ν is kinematic viscosity
that varies with temperature
=1.79×10−61+0.03369T+0.00021T2 with temperature T
in ˚C. Data was collected during winter season during day time in Bangladesh. The average water temperature was considered as 15˚C. For this temperature,
ν=1.15×10−6 m2⋅s−1.
Here,
s=ρs/ρw=2.6 which is determined from laboratory test of bed material (
ρs=2650 kg/m3,
ρw=1023 kg/m3 ).
Thus for OBM Ghat,
D∗=[(s−1)gυ²]1/3D50=2.38, when
D50=0.1042 mm and for VIP Ghat
D∗=[(s−1)gυ²]1/3D50=2.39, when
D50=0.1049 mm.
Critical Shield parameter for suspension obtained using Equation (1) is
θcr,suspension=0.099.
Then the critical velocity component Equation (2),
u∗,cr,suspension=0.0128 m/s. and
τ∗,cr,suspension=ρw×u2∗,cr,suspension=1023×1.56×10−4=0.167 N⋅m−2
Particle fall velocity ws
Particle fall velocity for marine environment is estimated using Equation (3) as follows:
ws=1.07×10−61.042×10−4[(10.362+1.049×2.383)12−10.36]=0.0073 m/s
From grain size analysis,
D50=0.1042 mm and
D=D84+D50+D163=0.1047 mm.
Rouse Number, N
OBM Ghat: peak flood
N=wsβκu∗=0.00731×0.4×0.074=0.247 and peak ebb
N=0.00731×0.4×0.085=0.215.
VIP Ghat: peak flood
N=wsβκu∗=0.00731×0.4×0.054=0.34 and peak ebb
N=0.00731×0.4×0.063=0.29.
Estimated reference concentration, Ca and Transport Parameter T
OBM Ghat
T(Peakebb)=τb−τcrτcr=u2∗−u2∗,cru2∗,cr=0.0852−0.012820.01282=43 and
T(Peakflood)=0.0742−0.012820.01282=32.4.
VIP Ghat
T(Peakebb)=u2∗−u2∗,cru2∗,cr=0.0632−0.012820.01282=23 and
T(Peakflood)=0.0542−0.012820.01282=16.8.
At OBM Ghat when T = 43,
Ca=0.015D50T1.5aD0.3∗=0.0150.0001042×431.50.16×2.390.3=2.12×10−3×2650=5.62 kg/m3=5620 mg/L.
At VIP Ghat when T = 23,
Ca=2198 mg/L and when T = 16,
Ca=1275 mg/L.
Evidence of rough bed condition
Analysis of measured data of the study site reveals that the current related bed shear velocity (0.074 - 0.085 m/s) is larger by 10% - 20% than critical velocity (0.0096 m/s). The median particle size of the bed material is 0.105 mm which is smaller than 0.5 mm. Water depth at the study site varies from 3 - 10 m which is >0.02 m. Transport stage parameter
T=23>3 and
D∗=2.39 which falls between 1 to 10. All these imply the presence of rough bed mostly with ripple formation. Flow depth where the velocity and sediment measurement were taken varies from 5 - 6 m. Therefore, for this study reference concentration C(a) is taken at average of 0.05h and 0.01h between 0.15 - 0.165 m.