Download the full PDF ➡
Read the full paper ⬇
ARTDOI: 10.1590/2317-4889201920180118
Brazilian Journal of Geology
BJ GEO
Experimental bedforms by saline density currents
Débora Koller1*
Abstract
Bedforms are sedimentary features that can be generated on the seafloor by the interaction between density currents and mobile beds. Developing knowledge about the hydraulic and sedimentary processes involved during these events is in the interest of research groups and oil companies. Because of the magnitude of the density currents in its natural environment and the challenge in collecting data, studies in laboratory are of great value. We present results of 29 experiments focusing in the bedform development generated by saline currents, testing two different sediment types and three grain sizes: melamine (245 μm-plastic) and sands (206 and 485 μm). We analyzed the cur- rent velocity and fractional density profiles as well as pictures taken during and after each run. Results have showed classical velocity and concentration profiles, for 8 subcritical and 21 supercritical currents, with densimetric Froude numbers (Frd) between 0.5 and 2.2. Some
correlations were identified, such as the decrease of the velocity peak height and increase in the mean velocity (with consequent reduction
of the current thickness), due to an increase of the concentration and/or flume slope. The occurrence of bedforms was more likely for high discharge and concentration rates of current density, which directly influence the Frd. Bedforms were classified according to the shear stresses values applied by the current to the bed, resulting in the generation of lower plane bed, ripples and dunes. Dunes and ripples were observed in supercritical flow conditions, which is a hydraulic scenario of bedforms generation not predicted by fluvial models. Thus, this
study demonstrated the existence of differences in generation and, consequently, the classification approach for density current bedforms, compared to those generated by river flows. To this fact is attributed the hydrodynamic (velocity and concentration profiles) and sediment transport differences between fluvial flows and density currents. Further studies may be carried out in order to constructing new concepts of bedforms generation by density currents.
KEYWORDS: Physical modelling; density currents; mobile bed; bedforms.
INTRODUCTION
Gravity flows are a type of flow whose movement takes place because of relative differences in the density between flow and ambi- ent water. These flows are called turbidity currents when density excess is promoted by the presence of sediments in turbulent sus- pension (Simpson 1982, Middleton 1993), and are responsible for transporting large amounts of sediment to the deep sea. Turbidity currents may be correlated with the formation of submarine can- yons and fans, and the sediment deposits created bythese flowscan form important hydrocarbon reservoirs (Meiburg & Kneller 2010, Talling et al. 2012). The interaction between these currents and mobilebeds,througherosion,transportanddepositionofsediments, may result in the generation and development of bedforms. These sedimentary structures give us information about current parame- ters and help us in the interpretation of ancient geological records. Fluvial plane bed, ripples, dunes and antidunes are the most studied bedforms and are associated with currents energy and
1Universidade do Rio Grande do Sul – Porto Alegre (RS), Brazil. E-mails: debora.koller@ufrgs.br, rafaelmanica@googlemail.com, alborges@ipg.ufrgs.br
2ExxonMobil – Houston (TX), United States. E-mail: juan.j.fedele@exxonmobil.com
*Corresponding author.
© 2019 The autors. This is an open access article distributed under the terms of the Creative Commons license.
mobile bed compositions. They can be developed under subcrit- ical, critical and supercritical hydraulic conditions, established by the Froude number (Simons & Richardson 1961, Kennedy 1969, Engelund & Fredsφe 1982, Van Rijn 1984), which con- siders the role of inertial and gravitational forces in a flow.
Monitoring the mechanisms of generation and development of these deep-sea processes is quite difficult due to the density currents magnitude (Heezen & Ewing 1952, Gorsline et al. 2000, Azpiroz- Zabala et al. 2017) and its capacity to destroy equipment used in data measurement. Therefore, it is extremely important to improve and to develop experimental studies in bedform generation by density currents. Many studies have attempted to gather and compare information on these processes (Parker et al. 1987, Raudkivi 1997, Puhl 2012, Cartigny & Postma 2016, Fedele et al. 2016), but there are still some information gaps with respect to bedform prediction and bringing flow fluvial concepts into the density currents realm can be problematic.
With the purpose of contributing to knowledge about bed- form generation by density currents in marine environments, the present study developed an experimental methodology for the generation of saline density currents and investigated the associated bedforms developed from the interaction between these flows and the mobile bed. Experimental saline currents are widely used to model sediment gravity flows, because of the complexity of turbidity currents (non-linear interaction of mixing, sediment entrainment and suspension, besides water-column stratification) (Parsons et al. 2007). Saline flows
Bed | S (º) | Experiment | Q (L min-1) | ρmix (kg m-3) |
0.5Mq1* | 380 | 1,016 | ||
0.5Mq2* | 384 | 1,026 | ||
0.5Mq3 | 383 | 1,042 | ||
0.5 | 0.5MQ1* | 506 | 1,016 | |
0.5MQ2 | 517 | 1,031 | ||
0.5MQ3 | 479 | 1,039 | ||
0.5MQ4 474 1,045 |
do not interact with their boundary and therefore conserve their buoyancy flux as they propagate downslope, and hence their parameters are more easily measured. Moreover, Sequeiros et al. (2010) assert that dissolved salt is capable of reproducing hydraulic and sedimentological processes of fine-grained tur- bidity currents. Furthermore, correlations on hydrodynamic, depositional and granulometric characteristics between the experimental currents and bedforms were performed.
generation by density currents. The experiments were named with information about slope, bed composition (melamine — M, beach sand — B or fluvial sand — F), mixture discharge (q ~ 380 L min-1 and Q ~ 505 L min-1) and density (1 ~ 1,015, 2 ~1,025, and 3 ~1,040 kg m-3).
Table 1. Experimental data: bed composition, slope (S, º), discharge (Q, L min-1) and mixture density (ρmix, kg m-3).
EXPERIMENTAL PROCEDURE
Melamine
Experiments were carried out at the Universidade Federal do Rio Grande do Sul, Brazil (at the hydraulic laboratory NECOD) and reported in detail by Koller (2016). The exper- iments were conducted in a 20 cm wide acrylic-sided flume with an 18 m long tilting section, submerged inside a long tank (Fig. 1) and filled with fresh water.
For each run, dense mixtures of water, salt and red dye were fed into the upstream part of the flume (Fig. 1D), where a 0.6 m long and 0.1 m high compartment was built. A 3 cm spillway was placed at the end of this compartment, distributing and guid- ing the mixture to the bottom of the flume. The water column above the entrance of the mixture varied from, approximately,
0.85 to 1.15 m, depending on the slope. This distance increases throughout the flume, reaching 1.39 cm of water column at its end (set 0.1 m from the bottom of the tank and approximately
1.5
0.5
1.5Mq1 381 1,016
1.5Mq2 383 1,024
1.5Mq3 383 1,041
1.5MQ1 519 1,016
1.5MQ2 520 1,025
1.5MQ3 519 1,043
0.5BQ1 504 1,010
0.5BQ2 503 1,024
4 m away from the long tank wall, to avoid flow reflection).
Mixtures were prepared in 5 m3 mixing tank (Fig. 1A) and flow discharges were measured at a rate of 4 Hz for all experiments using a magnetic flow meter (Siemens Sitrans
0.5BQ3 506 1,041
Beach sand
1.5Bq1 364 1,015
1.5Bq2 364 1,026
1.5Bq3 383 1,045
Mag 5100 W) installed in the supply line (Fig. 1C), after the
pump (KSB Megabloc, 32-160, 7.5 CV) (Fig. 1B). The pump flow was controlled by a register, so as to obtain the required constant flow rate throughout the experiment, which was obtained with a datalogger MyPCLab® connected to a com- puter. In addition, the water temperature of both mixture and tank water was then recorded before the experiment started.
1.5
0.5
1.5Bq4 325 1,041
1.5BQ1 510 1,015
1.5BQ2 516 1,025
1.5BQ3 461 1,042
0.5FQ1 509 1,016
0.5FQ2 504 1,025
In order to investigate several scenarios of bedform develop- ment, a wide range of different parameters were tested, includ- ing saline mixture density (1,010 < ρ (kg m-3) < 1,045), current discharge (325 < Q (L min-1) < 520), mobile bed composi- tion (sand and melamine) and flume slope (0.5 < S (º) < 1.5). The combination of these parameters totalized 29 experiments (Tab. 1) and enabled a more complete investigation of bedform
0.5FQ3 508 1,041
Fluvial sand
1.5Fq1 381 1,016
1.5 1.5Fq2 381 1,026
1.5Fq3 379 1,041
1.5Fq4 377 1,041
*Experiments realized with “melamine*”
Figure 1. Not to scale. Experimental schematic view: (A) mixture tank, (B) pump (C), flowmeter, (D) density current flume inlet, (E) mobile bed, (F) lateral pictures, (G) ultrasound velocity profiler (UVP), (H) siphons, (I) flume outlet valve.
Flow velocities were measured at 14 m from the current inlet with an Ultrasonic Velocity Perfilator (UVP — Duo MetFlow AS) at an acquisition rate of 12.5 Hz. The UVP profile con- sisted of 10 vertically distributed probes disposed at 0.8, 2.15, 4.95, 7.85, 10.75, 13.65, 18, 22.35, 26.7, 31.05 cm from the
mobile bed (Fig. 1G). The instrumental error associated to the UVP technique was below 1%, once the measurement condi- tions were with satisfactory signal-to-noise ratio (SNR – good echo from the flow).
Samples of the saline current were taken with a six-tube siphon system and used to extract vertical density distribu- tions (Fig. 1H). The tubes were settled up at 2, 5, 10, 13, 18 and 21 cm from the mobile bed and 15.5 m from the current
sands are sub-rounded with moderated to high spheric- ity. Melamine, on the other hand, was angular and with low sphericity (Krumbein 1963).
Regarding melamine, it was found after some experiments that much of its fine fraction was transported to the exit region of the flume. Because of that, it was decided to perform a new melamine bed sampling and grain size analysis, resulting in a larger average diameter (d50 = 310 μm). Melamine with- out this finer fraction was named melamine* and is also pre- sented in Table 2.
All three types of sediment used were classified according to Folk & Ward (1957) as moderately sorted (eqs. 1 and 2) with σ(Φ) from 0.53 to 0.73.
inlet. The excess of density was then measured from these sam-
ples using a refractometer ATAGO S28E 2~28%, which were converted to concentration values (C) through a calibration curve, with a mean error of 4.3%.
σ (Φ):
%84 −%16
4
× %95 −%5
6,6
(1)
The flow and bedform behavior could be observed during the entire experiment through a glass window, located 12 to 16 m from the flow inlet (Fig. 2).
Bedform (Fig. 1E) generation and development were recorded, as well as the current thickness, every second from the sidewall of the flume, using a Nikon D5000 camera (Fig. 1F). In order to observe bedform crests in plan view, the flume was slowly drained (Fig. 1I) after the end of each experiment and the deposit photographed from above, as described by Koller (2016). To avoid possible image distortion caused by the camera lens, approximately 15% of each side of the pic- ture were rejected in the image assembly. In addition, the scale provided the preserve of its real dimension, even with distor- tion of the image.
Melamine (SG: 1.5) and two types of sand, named here beach and fluvial sands (SG: 2.5), were used as erodible bed. Particle size distributions of the sediments were analyzed with a laser particle size analyzer, Cilas 1180, which reads grain sizes between 0.04 and 2,500 mm (Fig. 3 and Tab. 2).
Sediment morphoscopic properties (sphericity and round- ness) were assessed to reinforce grain-size properties of the material used as mobile bed. Both beach and fluvial
In which:
Φ:−log2 (d ) (2)
In which:
σ = the degree of selection values in relation to parameter Φ; Φ = the scale parameter;
%84 = the percentile 84;
%16 = the percentile 16;
%95 = the percentile 95;
%5 = the percentile 5;
d = the diameter of the particle (μm).
Figure 3. Particle size distribution curve from sands and from melamine, as well as median grain size diameter (d50, μm).
Table 2. Grain size data from beach and fluvial sands, and melamine, tested in this study.
Bed | ρS (kg m-3) | d10 | d50 d90 (μm) | d* | |
Melamine | 1,500 | 165 | 245 | 410 | 3.9 |
Melamine* | 1,500 | 169 | 310 | 487 | 4.9 |
Beach sand | 2,600 | 131 | 206 | 324 | 4.8 |
U: current mean velocity; H: thickness; ρCD: density.
Fluvial sand 2,600 208 480 790 11.3
Figure 2. Experimental saline current (experiment 0.5Mq1*, U = 0.14 m s-1, H = 0.26 m and ρDC = 1,006 kg m-3).
d*: dimensionless diameter; Melamine*: melanine sampled after some experiments and which grain size was re-analyzed.
DATA ANALYSIS
Velocity, concentration
in a fluid, where U is the mean current velocity (U, m s-1), H thickness (m), and υ is the kinematic viscosity (m s- 2).
UH
and depth calculation
Velocity (u) and saline concentration (c) profiles were mea-
Re = υ
(9)
sured in the flow direction in different vertical elevations from bed (z) during each experimental run and used to compute mean vertically averaged current parameters, as velocity (U, m s-1), concentration (C, g L-1) and thickness (H, m). These parameters were determined by Ellison & Turner (1959) and are presented in equations 3, 4 and 5.
Dimensionless grain size
This parameter is commonly used to compare sediments with different densities and is described by Van Rijn (1984) in Eq. 10, where d* is a dimensionless particle parameter, ρS is
the density of the sediment (kg m-3), ρDC is the mean density of
the current (kg m- 3), d50 is the mean grain size (m), and υ is the
∞ m−1 u +u
kinematic viscosity coefficient (m s-2), as shown in Equation 10.
UH = ∫ udz =∑ i i+1 (zi+1 − zi )
(3)
0 i=1
∞
2
m−1 u2+u 2
d* =
(10)
(ρ
S ρ
DC
−1 gd 3
)
13
50
υ2
U 2H = ∫ u2 dz =∑ i i+1 (zi+1 − zi )
(4)
0 i=1 2
∞ m−1 u c +u c
Mobile bed analysis
UCH = ∫ uc dz =∑
i i i+1 i+1 (zi+1 − zi )
(5)
Besides the density current and bed parameters previously
0 i=1 2
Shear velocity values (u*) were estimated from experimen- tal velocity data collected in the lower region of the profile, between bed and maximum velocities. The logarithmic profile method (Eq. 6) was used (Altinakar et al. 1996, Manica 2009).
presented, pictures taken laterally (during the experiments) and from above (after the experiment) were also used to assist in the recognition of the different bedform types generated by the experimental currents (Fig. 4), from which the wave- lengths and crest types of the bedforms could be measured.
u = 1 ln z u* k z0
(6)
Bedforms were classified with respect to three factors:
- sediment transport near the bed and the presence of sedi- ment in suspension during bedform development (verified
- In which:
- z = the vertical elevation (m); z0 = the bottom reference (m); u* = the shear velocity (m s-1);
- k = the Von Kàrman constant (0.41).
- Shear velocity (u*) and near-bed saline density of the cur- rent (ρDCb) then obtained were used to estimate bed shear stress values (τb), as shown in Equation 7.
- with the pictures taken during the experiments);
- bedform wavelengths;
- near-bed shear stress values.
RESULTS AND DISCUSSION
This study generated experimental currents in the sub- and supercritical regimes, as well as three different bedforms types, lower plane bed, ripples and dunes. Lower plane beds and rip- ples presented no visible amounts of suspended sediments
τb
ρDCb
u* =
(7)
during its evolution, as expected. Dunes (1.35 < Frd < 1.66), and in some cases ripples (1.07 < Frd < 1.70), were identi- fied in supercritical flows conditions, unlike the expected for
Densimetric Froude number and Reynolds number The parameters U, C and H allow the computation of the densimetric Froude number (Frd) of the flow (Eq. 8), which
is the ratio of inertial to gravitational forces.
fluvial bedforms but already observed in other experimental studies with density currents (Hand 1974, Fedele et al. 2016, Koller 2016). It seems that, if compared to fluvial bedforms,
ρDC − ρamb gH
U
ρamb
Frd =
(8)
In which:
g = the gravity acceleration (m s-2);
ρDC = the mean density of the current (kg m-3); ρamb = the density of the ambient water (kg m-3).
All density currents were fully turbulent, presenting Reynolds numbers (Eq. 9) between 21,903 and 61,882. This number represents the ratio of inertial to viscous forces
U: current mean velocity; H: thickness; ρCD: density.
Figure 4. Plan view of bedforms generated by experimental saline current (experiment 1.5Mq3, U: 0.28 m s-1, H: 0.19 m and ρDC: 1019 kg m-3).
the bedforms generated by density currents require larger flow velocities and, consequently, higher energy to be developed. Sediment composition also played an important role in bedform generation, due to its density and mean grain size. Meanwhile, these characteristics have a direct effect on the mobility of the particle, that is, erosion, suspension and depo- sition processes during the density current flow. Among the three types of tested bed material in this study (melamine and sands), melamine showed high mobility and seems more likely to enter and remain in suspension, either by its low density or
its lower fall velocity in relation to the sands.
Density current parameters
As previously presented in the Data Analysis section, gath- ering information on the vertical distribution of velocities and concentrations assists in the understanding on currents hydrodynamics and on the stresses applied by the flow on the bed, as shown below.
Although there have been challenges in maintaining the reproduction of the pre-defined parameters (flow rates — high and low — and densities — low, medium and high), there seems to have been no significant influences on the mean velocities and densities of the turbidity currents at the measured point (14 m from the current inlet). For high flow rates (Q ~ 505 L min- 1) the standard deviation (σ) of the mean flow velocities was 0.062 cm s-1 and, at low flow rates (~ 380 L min-1), 0.058 cm s-1. Moreover, low (~ 1,015 kg m-3), medium (~ 1,025 kg m-3) and high (~ 1,040 kg m-3) variations in input flow densities presented, respectively, 1.7, 3.0 and
5.3 kg m-3 of σ in the near bed density, in the measured point (15.5 m from the current inlet).
A summary of the more significant hydraulic parameters is given in Table 3.
Velocity profiles
Experimental saline currents from this study developed mean velocities ranging from 0.10 to 0.36 m s-1.
As previously documented (Fabian 2002, Sequeiros 2012, Puhl 2012), velocity profiles (Fig. 5) present lower values near the mobile bed due to current-bed interaction and a positive velocity gradient until a maximal point, defining the inner region, which is similar to a boundary layer flow (Fedele et al. 2016). Above the maximum velocity, the values con- tinue to decrease until they reach the mixing region (Stacey & Bowen 1988), where there is a more pronounced ambient water incorporation.
Moreover, Figure 5 (A and B) shows the pronounced influ- ence of the slope (0.5 and 1.5°) on the development of the flows and on the increase of the velocities along the vertical. These differences have an influence on the calculated mean velocities and hydraulic parameters.
Dimensionless profiles were also created using H and umax scaly parameters (Fig. 6), in order to better evaluate and com- pare vertical velocity distributions of each experiment. Figure 6
also shows a slight trend of most bedforms located in the right and bottom region of the graph, indicating that their occurrence is more likely for currents with lower distances
between the highest velocity and bed (0.08 < z (m) < 0.11) (thinner boundary layer), smaller thickness (z/H), and higher flume slopes (1.5º). This trend is approximately indicated in Figure 6 by a translucent light gray region.
Density profiles
Currents density profiles from each experiment are pre- sented in Figure 7 and indicate similar behavior to typical den- sity profiles for density currents of previous studies (Fabian 2002, Puhl 2012, Sequeiros 2012).
Experimental saline currents from this study developed densities reached values between 1,002.9 and 1,027.6 kg m-3. The highest current density values (Fig. 7) were found near the mobile bed. Along the vertical profile, these values tended to decrease until the current-water interface, where turbulence and current dilution were more intense
Experiments performed on a 0.5º flume slope developed currents with mean concentration values approximately constant
Table 3. Currents mean velocity (U, m s-1), density (ρCD, kg m-3), thickness (H, m) and Reynolds number (Re, -).
Exp | U (m s-1) | ρCD (kg m-3) | H (m) | Re (m) |
0.5Mq1* | 0.14 | 1,006.0 | 0.26 | 32,615 |
0.5Mq2* | 0.15 | 1,009.5 | 0.23 | 30,271 |
0.5Mq3 | 0.18 | 1,024.3 | 0.23 | 36,962 |
0.5MQ1* | 0.20 | 1,006.3 | 0.20 | 36,170 |
0.5MQ2 | 0.20 | 1,018.0 | 0.25 | 44,849 |
0.5MQ3 | 0.19 | 1,010.6 | 0.25 | 40,606 |
0.5MQ4 | 0.13 | 1,024.9 | 0.25 | 28,725 |
1.5Mq1 | 0.10 | 1,005.2 | 0.24 | 21,903 |
1.5Mq2 | 0.22 | 1,009.8 | 0.20 | 37,828 |
1.5Mq3 | 0.28 | 1,019.2 | 0.19 | 47,159 |
1.5MQ1 | 0.22 | 1,005.4 | 0.25 | 48,285 |
1.5MQ2 | 0.26 | 1,011.6 | 0.23 | 53,901 |
1.5MQ3 | 0.36 | 1,027.6 | 0.19 | 61,882 |
0.5BQ1 | 0.13 | 1,002.9 | 0.29 | 32,978 |
0.5BQ2 | 0.20 | 1,014.0 | 0.29 | 50,567 |
0.5BQ3 | 0.24 | 1,032.4 | 0.26 | 55,431 |
1.5Bq1 | 0.16 | 1,002.9 | 0.20 | 28,481 |
1.5Bq2 | 0.22 | 1,009.1 | 0.16 | 31,596 |
1.5Bq3 | 0.25 | 1,018.1 | 0.14 | 31,251 |
1.5Bq4 | 0.25 | 1,017.3 | 0.16 | 34,682 |
1.5BQ1 | 0.22 | 1,003.7 | 0.20 | 40,600 |
1.5BQ2 | 0.25 | 1,007.8 | 0.23 | 49,925 |
1.5BQ3 | 0.30 | 1,020.8 | 0.17 | 45,267 |
0.5FQ1 | 0.16 | 1,004.8 | 0.27 | 37,754 |
0.5FQ2 | 0.15 | 1,008.6 | 0.26 | 34,495 |
0.5FQ3 | 0.25 | 1,031.6 | 0.25 | 53,801 |
1.5Fq1 | 0.15 | 1,004.0 | 0.24 | 32,130 |
1.5Fq2 | 0.24 | 1,008.7 | 0.17 | 35,692 |
1.5Fq3 | 0.29 | 1,017.9 | 0.16 | 41,932 |
1.5Fq4 0.28 1,016.1 0.16 39,316
along the vertical profile, which differs from observations of Kneller & Buckee (2000) and Manica (2009) e.g., who tested turbidity currents. We attribute this fact to the hydraulic dif- ferences between conservative and non-conservative currents (saline flows do not interact with their boundary and therefore conserve their buoyancy flux as they propagate downslope). On the other hand, 1.5º slopes promoted concentration profiles
with a more apparent concentrated layer near the bottom and values rapidly decreasing toward the ambient water, probably due to the higher current velocities (between 0.36–0.22 m s-1 for Q group ~ 505 L min-1 and 0.29–0.15 m s-1 for q group
~ 380 L min-1) developed in this slope when compared to 0.5° (0.25–0.13 m s-1 to Q and 0.14–0.18 m s-1 to q) and due to the more intense current-ambient mixing.
Figure 7 shows a tendency of the higher the currents concen- tration, the more likely the presence of bedforms (filled points in the profiles). We observed that the higher the concentra- tion closest to the bed are the maximal velocities and, conse- quently, the higher are the shear stresses applied to the bed (Tab. 4). Indeed, 78% of the high concentration experiments generated bedforms and solely two of them did not modify the mobile bed.
Furthermore, dunes appeared only for slopes of S= 1.5° (a condition that influences the flow velocity and the stresses near the bed) and, more commonly in conditions of high density (from 1,025 to 1,045 kg m- 3). Differently, low con- centration density currents generated bedforms in only three experiments (0.5MQ1*, 1.5Mq1 and 1.5MQ1), while six pre- sented flat beds.
Densimetric Froude number (Fr ) and shear stress ( )
d b
Figure 5. Velocity profiles from all experiments, separated by slope. Note the influence of the slope in the vertical velocity distribution. Slopes of 0.5º resulted in profiles with maximal velocities between
0.15 and 0.31 m s-1, whereas 1.5º caused an increment in the maximal velocities (0.21 m s-1 < umax < 0.48 m s-1) and a decrease in the distance between these maximal velocities and the bed (~ 0.14 m to 0.5º against ~ 0.078 m to 1.5º).
Figure 6. Dimensionless velocity profiles presented according to the bed composition used in each experiment (melamine — circle, beach sand — triangle, and fluvial sand — square) and generated bedforms (black — dunes, medium gray — ripples, and light gray
— plane bed). Note that dunes (black points) are more likely to occur for high maximal velocity values (umax) and thinner boundary layers (distance between maximal velocities and bed).
The densimetric Froude number results (Tab. 4) ranged
from 0.5 to 2.1 (8 occurred in subcritical and 21 in supercrit- ical flow conditions), and the bottom shear stresses between
0.26 and 3.28 N m-2.
In order to compare the experiments composed with three types of sediments, we plotted the densimetric Froude num- ber (Frd) and the dimensionless particle parameter (d*) in a
diagram (Fig. 8).
In general terms, bedforms occured to high Frd values, regardless of the type of sediment used. It demonstrates the need of higher flow energies to obtain bedforms, if compared
with the subaerial bedform formation, where lower plane beds, ripples and dunes, for example, are formed under subcritical subaerial flows (Fr < 1) (Simon & Richardson 1961, Kennedy 1969, Engelund & Fredsφe 1982).
Figure 7. Saline concentration profiles. Filled points refer to currents that were capable to form ripples or dunes, whereas unfilled points refer to lower plane beds. Gray and black colors represent slopes of, respectively, 0.5 and 1.5º.
Also, the diagram shows that the Frd and d* results obtained for lower plane bed, ripples and dunes tend to cluster these three bedforms in similar regions of the diagram, which may aid in the classification of this sedimentary features. However, experiments 1.5Fq2 and 1.5BQ1 presented high Frd values (1.84 and 2.10, respectively) (upper part of the diagram) and were classified as lower plane beds because of the visual absence of bedload and suspended sediment transport during the exper- iments (as seem in the next section – Bedforms).
Dunes and, in some cases, ripples were observed for supercritical flows, although they are not found in the subaer- ial realm. Hand (1974) identified the coexistence of experi- mental upstream-migrating ripples and downstream-migrat- ing antidunes, under supercritical flow conditions. More recently, Spinewine et al. (2009) reached similar conclusions, indicating greater diversity of bedforms for density currents when compared to subaerial flows. In addition, Fedele et al. (2009, 2016) have discovered a new type of bedform whose
Table 4. Bedforms according to the Frd and shear stress values (τb), and bedforms dimensions (wavelength — λ and amplitude — η).
genesis and development resembles small wavelength down- stream-migration antidunes (opposite to what usually occurs in fluvial flows), due to the interface interaction between flow and bedforms.
Bedforms
The classification of bedforms generated in this study and their relationship with the mobile bed composition are described below.
Lower Plane Bed
Lower plane bed (Fig. 9) were generated in 13 of 29 exper- iments. In five of those (0.5Mq1*, 0.5Mq2*, 0.5FQ2, 0.5FQ3 and 0.5BQ2), flows were subcritical, but very near to critical (0.92 < Frd < 0.98). The remaining experiments (0.5FQ1,
1.5Fq1, 1.5Fq2, 0.5BQ1, 0.5BQ3, 1.5Bq1, 1.5Bq2 and 1.5BQ1)
resulted in supercritical currents (1.04 < Frd < 2.10).
During these experiments, and based on the picture’s
analysis, we did not notice near bed sediment suspension or movement, as expected for this bed state.
Upper plane bed is another type of plane bed that nor- mally occurs along with increasing flow energy and the sup- pression of ripples and dunes due to the high shear stresses applied over the mobile bed.
Experiments 1.5BQ1 and 1.5Fq2 presented high values of bottom shear stress (1.65 and 3.28 N m-2, respectively) and
Figure 8. Densimetric Froude (Frd) number versus dimensionless grain size (d*), and the classification of likely bedforms occurrence. Filled shapes indicate the occurrence of bedforms (ripples and dunes), whereas unfilled shapes refer to lower plane beds. Bed material were melamine (circle), beach sand (triangle), and fluvial sand (square).
Bedform Exp | Frd (-) | τb (N m-2) | λ (cm) | η (cm) |
0.5Mq1* | 0.98 | 0.17 | NO | NO |
0.5Mq2* | 0.92 | 1.06 | NO | NO |
0.5BQ1 | 1.12 | 0.34 | NO | NO |
0.5BQ2 | 0.94 | 0.64 | NO | NO |
0.5BQ3 | 1.04 | 1.24 | NO | NO |
1.5Bq1 | 1.71 | 0.72 | NO | NO |
LPB 1.5Bq2 | 1.65 | 1.07 | NO | NO |
1.5BQ1 | 2.10 | 1.65 | NO | NO |
0.5FQ1 | 1.22 | 1.56 | NO | NO |
0.5FQ2 | 0.92 | 0.72 | NO | NO |
0.5FQ3 | 0.97 | 0.74 | NO | NO |
1.5Fq1 | 1.32 | 0.88 | NO | NO |
1.5Fq2 | 1.84 | 3.28 | NO | NO |
0.5Mq3 | 0.79 | 0.35 | 9.4 | 0.73 |
0.5MQ1* | 1.55 | 1.70 | 9.0 | < 0.5 |
0.5MQ2 | 0.91 | 0.29 | 8.9 | < 0.5 |
0.5MQ3 | 1.07 | 0.46 | 9.0 | 0.5 |
R 0.5MQ4 | 0.51 | 0.65 | 8.1 | < 0.5 |
1.5Mq1 | 0.81 | 0.33 | 9.5 | 1.02 |
1.5Bq3 | 1.54 | 0.67 | 5.9 | < 0.5 |
1.5Bq4 | 1.44 | 0.74 | 5.6 | 0.70 |
1.5BQ2 | 1.70 | 1.76 | 6.5 | < 0.5 |
1.5BQ3 | 1.56 | 1.48 | 5.4 | 0.54 |
1.5Mq2 | 1.45 | 0.84 | 11.0 | 1.02 |
1.5Mq3 | 1.35 | 2.93 | 17.5 | 1.42 |
1.5MQ1 | 1.62 | 0.94 | 9.3 | 0.62 |
D 1.5MQ2 | 1.5 | 1.21 | 13.0 | 1.5 |
1.5MQ3 | 1.53 | 1.64 | 8.2 | 1.5 |
1.5Fq3 | 1.66 | 2.25 | 9.6 | < 0.5 |
Q: discharge, ρ : density; Fr : densimetric Froude; τ : shear stress
CD d b
value.
*Experiments realized with “melamine*”
Figure 9. Lower plane bed from experiment 1.5Bq1 (Q = 364 L min-1, ρDC = 1,002.9 kg m-3, Frd = 1.71, τb= 0.72 N m-2). These bedforms were observed as result of low shear stresses and no sediment motion or transport.
densimetric Froude number (2.18 and 1.84, respectively) when compared to the other obtained lower plane beds. These val- ues could indicate a possible transport of sediments and the misclassification of these forms in upper plane bed.
However, it should be noted that no sediment transport (neither bedload nor suspension) was observed during these tests and, because of that, 1.5BQ1 and 1.5Fq2 were also clas- sified as lower plane beds.
Ripples
During the formation of this bedform type, the interaction between current and mobile bed did not produced bed sediment suspension. The grains were transported as bedload, developing forms with gentle upstream faces and more steep downstream faces, like already described by Fedele et al. (2009). All bedforms presented downstream migration, wavelength (λ) between 5.4 and 13 cm, and height (η) values between 0.5 and 1.5 cm(Tab. 4).
Ripples (Fig. 10) were observed in 10 experiments, of which four (0.5Mq3, 0.5MQ2, 0.5MQ4 and 1.5Mq1) occurred under subcritical flow conditions (0.51 < Frd < 0.91). The
six remaining experiments (0.5MQ1*, 0.5MQ3, 1.5Bq3,
1.5Bq4, 1.5BQ2 and 1.5BQ3) occurred in supercritical flow condition (1.07 < Frd < 1.70).
Dunes
Dunes were observed in six experiments, 1.5Mq2, 1.5Mq3, 1.5MQ1, 1.5MQ2, 1.5MQ3 and 1.5Fq3. These bedforms were developed under high shear stress values, which are applied by the flow over the sedimentary bottom (0.84 < (N m-2)
< 2.93, Tab. 4). Because of that, during the generation of these bedforms, erosion, suspension and transport of the sediment present in the bed were visible during the experiment (Fig. 11) and along the observation window.
Figure 11 shows dune generation and evolution from exper- iment 1.5Mq3, where it is possible to note suspension of the melamine present in the bed and its transport along the flume.
At the lee side of the bedforms, the movement of the grains is more disordered. The flow separates at the dune crest, gen- erating a shear layer that plays a crucial role in the turbulent transport of momentum and energy. It was verified by the observation of bed thickness decrease during the experiment and by the large amount of sediment present at the exit of the flume after the end of these experiments.
Then, the flow reattachment seems to occur in the beginning of the next bedform (Fig. 11), where there is bed removal with posterior sediment release in the next stoss side of the dune. Grain movement mechanisms involved in the bedform development, as reported above, are similar to those cited in the literature for subaerial dunes (Simons & Richardson 1961,
Kennedy 1963, 1969).
Figures 12 and 13 show pictures taken after experiments
1.5BQ3 and 1.5Mq3. These experiments produced ripples and dunes, respectively, and it is possible to observe differ- ences between its wavelengths (λ) (Figs. 12 and 13) and heights (η) (Fig. 11).
Dunes formed in this study are larger than ripples (Fig. 12), with λbetween 8.2 and 17.5 cm and ηbetween 1.02 and 1.42 cm.
CONCLUSIONS
This paper reported 29 experimental saline currents in order to simulate bedforms generation and investigate its relation- ship with hydraulic, grain size and morphological information about the current and mobile bed. The results demonstrate the influence of mean velocity, concentration and slope angle in
Q: discharge, ρCD: density; Frd: densimetric Froude; τb: shear stress value. Figure 11. Dunes from experiment 1.5Mq3 (Q = 383 L min-1, ρDC = 1,019.2 kg m-3, Frd = 1.35, τb = 2.93 N m-2). High shear stresses are applied at the stoss-side of the bedform, where there is sediment erosion. The flow detaches at the dune crest and at the leeside of these bedforms, the movement of the grains is more disordered.
Figure 12. Plan view from experiments 1.5BQ3 (ripples) and
1.5Mq2 (dunes). While ripples crests are commonly perpendicular (transverse) to the flow, dunes present sinuous crest.
Frd: densimetric Froude; ρCD: density; τb: shear stress value.
Figure 10. Ripples from experiments 1.5BQ2 (Frd = 1.7, ρDC = 1,007.8 kg m-3 and τb= 1.76 N m- 2) and 1.5BQ3 (Frd = 1.56, ρDC = 1,020.8 kg m-3, and τb = 1.48 N m-2)
Figure 13. Lateral view from experiments 1.5BQ2 (ripples) and
1.5Mq3 (dunes). Ripples were generated under lower shear stresses when compared to dunes and presented wavelength of ~ 6 cm, whereas dunes ~ 19 cm.
bedform generation, which occur more frequently according to the following conditions:
- High values of density current concentration;
- High values of density current maximal velocity (umax), relative to the mean velocity;
- Steep slopes have shown strong influence over the mean velocity of the current;
- Ripples and dunes were found to be developed in super- critical densimetric Froude flows;
- Sediment material (ρs) and grain size (d50) have distinct influences on the bedforms generation;
The experimental generation of ripples and dunes under sub and supercritical densimetric Froude condi- tions show the versatility of these bedforms occurrence and may aggregate knowledge in the interpretation of bedforms observed in natural environments, particularly in deepwater realms.
ARTICLE INFORMATION
Manuscript ID: 20180118. Received on: 10/23/2018. Approved on: 03/06/2019.
D. K. wrote the first draft of the manuscript, prepared all the figures and tables and carried out the corrections requested by the reviewers;
R. M. provided advisorship regarding hydraulic processes and sediment transport, and also improved the manuscript through corrections and suggestions. A. B. brought to light several questions about the mechanism of bedform’s development and its relationship with hydraulic parameters. J. F. cooperated to compile all the initial ideas brought by the D. K. and helped to improve the paper after the corrections requested by the reviewers.
Competing interests: The authors declare no competing interests.
REFERENCES
Altinakar M.S., Graf W.H., Hopfinger E.J. 1996. Flow structure in turbidity currents. Journal Hydraulic Research, 34(5):713-718. https://doi. org/10.1080/00221689609498467
Azpiroz-Zabala M., Cartigny M.J.B., Talling P.J., Parsons D.R., Sumner E.J., Clare M.A., Simmons S.M., Cooper C., Pope E.L. 2017. Newly recognized turbidity current structure can explain prolonged flushing of submarine canyons. Science Advances, 3(10):e1700200. https://doi. org/10.1126/sciadv.1700200
Cartigny M.J.B., Postma G. 2016. Turbidity current bedforms. In: Guillén J., Acosta-Yepes J., Chiocci F.L., Palanques A. (Eds.). Atlas of bedforms in the Western Mediterranean. New York, Springer.
Ellison T.H., Turner J.S. 1959. Turbulent entrainment in stratified flows. Journal of Fluid Mechanics, 6(3):423-448. https://doi.org/10.1017/ S0022112059000738
Engelund F., Fredsφe J. 1982. Sediment ripples and Dunes. Annual Review of Fluid Mechanics, 14:13-37. https://doi.org/10.1146/annurev. fl.14.010182.000305
Fabian S. 2002. Modelagem Física de Correntes de Densidade Conservativas em Canal de Declividade Variável. 2002. 107f. Dissertation, Instituto de Pesquisas Hidráulicas, Universidade Federal do Rio Grande do Sul, Porto Alegre.
Fedele J.J., Guentzel K., Hoyal D.C. 2009. Experiments on Bedforms Created by Density Currents. In: River, Coastal and Estuarine Morphodynamics Symposium, 9., Boca Raton. Annals… Florida: CRC Press: 2009, p. 833-840.
Fedele J.J., Hoyal D.C., Barnaal Z., Tulenko J., Awalt S. 2016. Bedforms created by gravity flows. In: Budd D.A., Hajek E.A., Purkis S.J. (Eds.). Autogenic Dynamics and Self-Organization in Sedimentary Systems. SEPM Special Publication, n. 106.
Folk R.L., Ward W.C. 1957. Brazos River bar: a study in the significance of gran size parameters. Journal of Sedimentary Petrology, 27(1):3-26. http:// dx.doi.org/10.1306/74D70646-2B21-11D7-8648000102C1865D
Gorsline D.S., de Diego T., Nava-Sanchez E.H. 2000. Seismically triggered turbidites in small margin basins: Alfonso Basin, Western Gulf of California and Santa Monica Basin, California Borderland. Sedimentary Geology, 135(1):21-35. http://dx.doi.org/10.1016/S0037-0738(00)00060-9
Hand B.M. 1974. Supercritical flow in density currents. Journal of Sedimentary Research, 44(3):637-648. https://doi. org/10.1306/74D72AB3-2B21-11D7-8648000102C1865D
Heezen B.C. & Ewing M. 1952. Turbidity currents and submarine slumps and the 1929 Grand Banks earthquake. American Journal of Science, 250(12):849-873. https://doi.org/10.2475/ajs.250.12.849
Kennedy J.F. 1963. The mechanics of dunes and antidunes in erodible- bed channels. Journal of Fluid Mechanics, 16(4):521-544. https://doi. org/10.1017/S0022112063000975
Kennedy J.F. 1969. The Formation of Sediment Ripples, Dunes, and Antidunes. Annual Review of Fluid Mechanics, 1:147-168. https://doi. org/10.1146/annurev.fl.01.010169.001051
Kneller B.C., Buckee C. 2000. The Structure and Fluid Mechanics of Turbidity Currents: a Review of Some Recent Studies and Their Geological Implications. Sedimentology, 47(s1):62-94. https://doi. org/10.1046/j.1365-3091.2000.047s1062.x
Koller D.K. 2016. Estudo experimental de formas de fundo por correntes de densidade salina em canal de fundo móvel. Dissertation, Programa de Pós- Graduação em Recursos Hídricos e Saneamento Ambiental, Instituto de Pesquisas Hidráulicas, Universidade Federal do Rio Grande do Sul. Porto Alegre.
Krumbein W.C. 1963. Stratigraphy and Sedimentation. 2. ed. San Francisco, Freeman, 660 p.
Manica R. 2009. Geração das correntes de turbidez de alta densidade: condicionantes hidráulicos e deposicionais. Thesis, Instituto de Pesquisas Hidráulicas, Universidade Federal do Rio Grande do Sul, Porto Alegre.
Meiburg E. & Kneller B. 2010. Turbidity Currents and Their Deposits. Annual Review of Fluid Mechanics, 42:135-156. https://doi.org/10.1146/ annurev-fluid-121108-145618
Middleton, G. V., 1993. Sediment Deposition from Turbidity Currents. Annual Review Earth Planet. Sci., 21:89-114. https://doi. org/10.1146/annurev.ea.21.050193.000513
Parker G., Garcia M., Fukushima Y., Yu W. 1987. Experiments on turbidity currents over an erodible bed. Journal of Hydraulic Research, 25(1):123-147. https://doi.org/10.1080/00221688709499292
Parsons J.D., Friedrichs C.T., Traykovski P.A., Mohrig D., Imran J., Syvitski J.P.M., Parker G., Puig P., Buttles J.L., García M.H. 2007. The Mechanics of Marine Sediment Gravity Flows. In: Nittrouer C., Austin J., Field M., Kravitz J., Syvitski J., Wiberg P. (Eds.). Continental Margin Sedimentation: From Sediment Transport to Sequence Stratigraphy, v. 37, p. 275-334. IAS Special Publication. https://doi.org/10.1002/9781444304398.ch6
Puhl E. 2012. Morfodinâmica e condição de equilíbrio do leito sob a ação de correntes de turbidez. Thesis, Instituto de Pesquisas Hidráulicas, Universidade Federal do Rio Grande do Sul, Porto Alegre, 155 f.
Raudkivi A.J. 1997. Ripples on stream bed. Journal of Hydraulic Engineering, 123(1):58-64. https://doi.org/10.1061/(ASCE)0733-9429(1997)123:1(58)
Sequeiros O.E. 2012. Estimating turbidity current conditions from channel morphology: a Froude number approach. Journal of Geophysical Research, 117(C4). https://doi.org/10.1029/2011JC007201
Sequeiros O.E., Spinewine B., Beaubouef R., Sun T., Garcia M., Parker
G. 2010. Bedload transport and bed resistance associated with density and turbidity currents. Sedimentology, 57(6):1463-1490. https://doi. org/10.1111/j.1365-3091.2010.01152.x
Simons D.B., Richardson E.V. 1961. Forms of Bed Roughness in Alluvial Channels. Journal of the Hydraulics Division, 87(3):87-105.
Simpson J.E. 1982. Gravity currents in the laboratory, atmosphere, and ocean. Annual Review Fluid Mechanics, 14(2):213-234. https://doi. org/10.1146/annurev.fl.14.010182.001241
Spinewine B., Sequeiros O., Garcia M., Beaubouef R., Sun T., Savoye B., Parker G. 2009. Experiments on wedge-shaped deep sea sedimentary
deposits in mini basins and/or on channel levees emplaced by turbidity currents. Part II. Morphodynamic evolution of the wedge and of the associated bedforms. Journal of Sedimentary Research, 79(8):608-628. https://doi.org/10.2110/jsr.2009.065
Stacey M.W., Bowen A.J. 1988. Vertical structure of density and turbidity currents: Theory and observations. Journal Geophysical Research, 93(C4):3528-3542. https://doi.org/10.1029/JC093iC04p03528
Talling P.J., Masson D.G., Sumner E.J., Malgesini G. 2012. Subaqueous sediment density flows: Depositional processes and deposit types. Sedimentology, 59(7):1937-2003. https://doi. org/10.1111/j.1365-3091.2012.01353.x
Van Rijn L.C. 1984. Sediment transport, part III: bed forms and alluvial roughness. Journal of Hydraulic Engineering, 110(12):1733-1754. https:// doi.org/10.1061/(ASCE)0733-9429(1984)110:12(1733)