Next Article in Journal
Investigation of the Sloshing Behavior Due to Seismic Excitations Considering Two-Way Coupling of the Fluid and the Structure
Previous Article in Journal
Urban Spatial Configuration and Functional Runoff Connectivity: Influence of Drainage Grid Density and Landscape Metrics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantifying the Variability in Flow Competence and Streambed Mobility with Water Discharge in a Gravel-Bed Channel: River Esva, NW Spain

1
INDUROT, University of Oviedo, Campus de Mieres, s/n 33600 Mieres, Spain
2
Department of Geology, University of Oviedo, c\Jesús Arias de Velasco, s/n 33005 Oviedo, Spain
*
Author to whom correspondence should be addressed.
Water 2019, 11(12), 2662; https://doi.org/10.3390/w11122662
Submission received: 16 November 2019 / Revised: 12 December 2019 / Accepted: 13 December 2019 / Published: 17 December 2019
(This article belongs to the Section Hydraulics and Hydrodynamics)

Abstract

:
Streambed mobility in gravel-bed rivers is largely controlled by the rate at which particles with different grain sizes are recruited from the riverbed into the bed load. In this paper, we present a study case in which we explored this question, based on combining field observations using painted plots and the grain size analysis of a large flood sediment deposit in the River Esva, northwest Spain, and the generalized threshold model (GTM) competence model developed by Recking. The main aim was to accomplish a complete characterization of streambed mobility in this river. The obtained results suggest the large potential of the GTM model compared to previous competence models when searching for the quantitative description of particle entrainment and streambed mobility in the River Esva. We observed how the grain size of the bed load in the River Esva tended to be closer to that of the sub-armour bed material during large floods, while moderate magnitude flows tended to carry a relatively fine bed load. Additionally, we compared our results with previously published field observations on flow competence. This comparison outlined the large degree of site specificity in the links between grain size of the bed load and that of the bed material.

1. Introduction

Unravelling the complex interplays between flow strength, grain-size distribution (GSD) of the bed load, and that of the bed material is a central issue in river morphodynamics [1,2]. This topic has been classically approached in fluvial geomorphology and sediment transport studies through the analysis of the thresholds for bed particle entrainment as long as the ability of the stream to displace riverbed particles controls the GSD of the bed load in alluvial channels [2,3,4,5]. In this regard, the inception of motion for streambed particles also exerts significant control over the area of bed surface that is actively involved in bed load motion and the vertical extent of scour and fill processes. Thus, particle entrainment is closely related to the overall streambed stability and the spatial extent of bed disturbance caused by sediment erosion and deposition during floods [6], which ultimately represent important drivers of aquatic habitat quality in river ecosystems [7,8,9,10].
Additionally, particle entrainment is closely related to the “flow competence” concept, introduced first by Gilbert and Murphy (1914) [11] and defined based on the size of the largest clast moved by a given flow. The idea of river competence relies on assuming that the coarsest clast transported by a streamflow provides a measure of the stream’s ability to entrain the bed sediment. Consequently, the measure of the largest clast size of fluvial deposits has been a common and important tool for most geologists. This procedure was frequently used as an instrument to infer velocities, discharges, mean stresses, and depths of flows that transported the particles found in fluvial conglomerates [12,13,14,15,16,17,18,19], or in the evaluation of extreme flood hydraulics [20]. However, the use of the largest moved clasts in this way, i.e., as a technique to infer differences in flow strength and stream’s ability to convey sediment between different sites, involves the implicit assumption that finer grains are always entrained at lower flows than coarser grains and over a relatively broad range of flow discharges [21,22]. This implies ignoring important physical aspects of sediment entrainment that nuance the (in principle) expectable “weight selective” assumption for particle entrainment [2,23,24].
Consequently, there is a clear interest in better understanding the links between the GSD of the bed material, the streambed surface, the bed load, and the spatial–temporal patterns of streambed mobility in coarse bed rivers. Previous works have already explored the relationship between the grain size of the bed load, sampled using sediment traps and/or different samplers, and the grain size of the bed material [1,13,18,25]. However, in many cases, these observations were made before some of the more modern bed load transport and flow competence models were developed (i.e., Parker (1990) [26]; Wilcock and Crowe (2003) [27]; Recking, (2016) [2]), so the former could still benefit from new field testing. Additionally, available field measures of the GSD of the bed load were carried out in many cases during moderate and regular bed load transport episodes (when sampling bed load with Helley-Smith and other kinds of samplers is easier), whereas information about relatively strong events is still scarce. In summary, more field studies are needed in order to better grasp several questions related to streambed mobility, such as understanding whether the GSD of the bed material proxies that of the bed load or not, or how the bed state modulates the relationship between particle entrainment and flow strength.
In this regard, the current paper presents a case study accomplished in a gravel bed stream from northwest Spain, where an extreme flood episode occurred in January 2019. It left important sediment deposits, and its GSD was sampled in the field. Additionally, painted plots were employed in order to monitor streambed mobility during regular, low-magnitude floods occurring before this major flood episode. Both sources of data were used in order to calibrate a streambed mobility model, based on Recking (2016) [2], which, in turn, was used for modelling the GSD of the average annual bed load and the temporal patterns of streambed mobility. The obtained results highlight how the grain size of the bed load in the River Esva tends to be closer to that of the sub-armour bed sediment during large floods, while moderate magnitude flows carry a relatively finer bed load. These results show the potential of coupling a competence model (such as the GTM model [2]) with field observations in order to predict streambed mobility. We also compared our measures and estimates with the bed load information available in Milhous (1973) [28], Andrews (1994) [29], Lisle (1995) [25] and Whittaker and Potts (2007) [1]. This comparison outlined how the relations between grain size of the bed load and that of the bed material can be largely variable with flow strength and among different sites.

2. Rationale

2.1. Physical Basis for Flow Competence Models

Flow competence and streambed mobility are closely related to particle entrainment and the hydraulic thresholds for incipient sediment motion. In this regard, classical research tried to relate shear stress and/or flow velocity to the size of the largest particles set in motion [30,31,32,33]. Initially, threshold stresses for sediment entrainment were related to particle size [30]:
τ c = τ c * × ( ρ s ρ ) × g × D
where ρ is the density of water, ρs that of sediment, g is the gravity acceleration, D the grain size of the bed sediment, τc is the critical bed shear stress for incipient motion (dimensional), and τc* is the critical dimensionless Shields stress, for which a large range of values have been proposed, normally between 0.03 and 0.06 [23,34]. Once known the value of this parameter, the grain size of the particles that are mobilized with increasing flow discharges (Dmax) are easy to estimate:
D m a x = τ τ c * × ( ρ s ρ ) × g
The previous model Equation (2) holds well in the ideal case of uniform sized particles. In the case of a graded mixture, some works revealed how the effective stress acting on a given grain is strongly influenced by the relative size of the grain within the mixture [35,36]. Sand content seems to influence gravel mobility by reducing the flow resistance and intergranular friction [11,27,37]. Clast arrangements and imbrications may also have a large impact on streambed mobility [38]. As a result, a large variability of situations has been documented regarding the relative mobility of each size class of a sediment mixture in gravel-bed rivers with several studies reporting size selective mobility over a significant range of flow [13,15,28,39,40,41], while others suggested that all sizes in a sediment mixture may begin moving over a range of narrow flow conditions [42,43,44,45,46]. So-called hiding functions were proposed in order to mathematically describe this range of situations:
τ c i * = τ c R e f * × ( D i D r e f ) α
where τci is the critical Shields parameter for incipient motion of a given diameter Di and τcref the critical Shields parameter for incipient motion of a reference particle diameter Dref; α is a coefficient whose value varies from 0 in the case of a total size selective behaviour (threshold conditions only depend on particle size) to −1 in case of perfect equal mobility conditions (all sizes are entrained at the same threshold stress, independent of their size). Rearranging Equation (3), we can estimate again the size of the largest clast moved during floods:
D m a x = D r e f × ( τ c i * τ c R e f * ) 1 1 + α
In spite of the great improvement hiding functions (such as Equations (3) and (4)) provide for modelling flow competence and streambed mobility, there are still some open issues. For instance, some researchers have found that a simple power law (Equations (3) and (4)) is not enough to represent the complete relation for any particle size critical stresses in gravel mixtures [34,47]. In addition, gravel-bed rivers are systematically armored [48,49,50,51,52], i.e., the streambed surface is commonly coarser than the underlying sub-armor deposit. The potential changes in the shape of the hiding parameter α with armor break-up are not considered when using a single hiding function. Additionally, the entire bed load GSD cannot be computed directly from Equation (4) (which only supplies the coarsest mobile size) but must be deduced from the bed load rates computed for each size class [26,27,34,49]. This commonly requires the use of an indirect scaling approach [3], implying that the same bed load function can capture (on average) the transport of each individual size class for a given shear stress, an idea that still needs field validation [2].

2.2. Recking’s (2016) GTM Model

Because of some of the issues stated above, Recking (2016) [2] proposed a new flow competence model (named “generalized threshold model”, GTM), as an alternative approach for characterizing streambed mobility and computing the GSD of the bed load. The GTM model computes the GSD of the bed load after modelling which fraction of sediment in each size range is entrained from and which fraction remains static on the riverbed. It does so without assuming a particular shape for the bed load function, in contrast to common approaches based on the hiding function (i.e., Parker (1990) [26]; Wilcock and Crowe (2003) [27]). The GTM workflow operates in two steps. First, it computes the maximum competent size from the following expression (see more details in Recking [2]):
M = 84 × ( τ 84 * τ c 84 * ) β
where M is the Mth percentile of the surface GSD corresponding to the largest mobile class, τ84* is the dimensionless Shields stress estimated based on the 84th percentile of the surface GSD (D84), τc84* is the critical Shields stress for entraining the surface D84. The use of the 84-th percentile as reference parameter in this model is meaningful, insofar as the inception of motion of the coarser particles defines the moment when the armour layer starts to become disorganized [53,54]. Then, a transport stage (τ84*/τc84*) close to 1 defines the threshold conditions when the armour layer destabilizes. In this regard, the β exponent establishes the rate at which progressively coarser grains are set in motion with increasing flow strengths. It may take a different value depending on whether or not the armour layer destabilizes. A value of ~0.25 is suggested by Recking [2] when the armour layer is broken (τ84*c84* > 1) which implies assuming full mobility conditions at τ84*c84 > 2 as Wilcock and McArdell (1997) [46] reported in the flume. In cases where the armour layer is not broken (τ84*c84 < 1), there is more uncertainty in choosing the right value for β as it probably depends on the degree of bed armoring/pavement. High β simulates an abrupt increase in competent size when armour breaks (i.e., a well-paved stream), whereas a low value models the case of a bed with a more progressive transition from partial to full mobility conditions.
The second step of the GTM model is calculating the fraction (φ) of bed particles mobilized in individual size ranges (i) of the bed surface:
φ i = 0.01 + ( 1 ( i M ) γ )
where γ can be estimated from:
γ = 0.1 + 0.5 × ( τ 84 * τ c 84 * ) γ 2
where γ2 is a parameter that needs to be calibrated. The combination of Equations (6) and (7) is very flexible and allows to model a broad range of streambed behaviors, such as progressive mobilizations of surface particles or more abrupt transitions from stability to mobility. Thus, a right γ2 parameter, based on previous empirical information about riverbed behavior, allows the proper modelling of a given study reach. This is one of the strongest points of the GTM model: the fact of identifying the key parameters (β and γ2) driving streambed behavior and constraining the competence model. The definition of these parameters from field observations permits to estimate the GSD of the bed load with the following expression:
F i = φ i j = i M φ j
In summary, there are three main advantages of the GTM model compared to previous approaches to flow competence: (1) its flexibility, which allows to model a large range of particle entrainment and streambed stability conditions; (2) the model is well-rooted into our current knowledge on sediment transport mechanics and considers some situations that were not adequately grasped by previous flow competence models; i.e., the change in entrainment conditions with armour break up; and (3) while most of the previous flow competence models focused mainly on computing the mobilized grain sizes, the GTM model also allows the estimation of the mobile fraction in each size range which is at least as important as predicting the maximum transported grain in terms of evaluating streambed disturbance (interesting, for example in many ecological applications). In spite of all these advantages, there have been no previous attempts to take advantage of the GTM model for the characterization of streambed mobility in specific field cases (Vázquez-Tarrío et al. (2019) [55] being the only exception to the best of our knowledge).

2.3. Proposed Workflow

The GTM model provides a well-suited practical framework for streambed mobility studies. Hence, we propose a field workflow adapted to the characterization of streambed mobility in gravel-bed rivers, which combines the flexibility of the GTM model with field observations (Figure 1). It combines two different sources of field data. Painted tracers, which can be employed in order to characterize streambed mobility at low transport conditions, and, as consequence, to obtain the best value for the β exponent in Equation (5) when τ84*/τc84* < 1. The second source of data comes from collecting in the field the GSD information of sediment deposited during large floods. This kind of data helps to calibrate Equation (6) by forward modelling: by changing the γ parameter in Equation (6) until obtaining an adequate value to fit the field measures. Finally, the γ parameter obtained from flood deposits could be combined with the field identification of incipient motion conditions in order to calibrate the γ2 parameter in Equation (7). This paper illustrates the implementation, potential and capabilities of this workflow, using the River Esva (NW Spain) as study case.

3. Study Site

The River Esva is located in the northern watershed of the Cantabrian Mountains, in NW Spain (Figure 2). The climate is temperate-oceanic, with average monthly rainfalls ranging between 52 and 146 mm. The River Esva drains a surface of 464 km2. The flood regime is perennial and rainfall dominated, with larger discharges in winter. Average annual discharge is 10.5 m3/s, with average maximum and minimum annual discharges of 19.2 and 4.5 m3/s, respectively.
We selected a ~1 km study reach located close to the village of Trevias (Figure 3). It is a single and straight river reach with an average ~0.003 channel slope and 25 m width. This reach was embanked in 1993 along its entire length in order to protect Trevias from overbank flows.
A lateral gravel bar develops on the left bank in the downstream part of the study reach. Streambed sediment is mainly composed of sub-angular and rod shaped sandstone gravels and cobbles. One square plot of blue-painted stones was deployed on this gravel bar in November 2016 (hereinafter called PP1 painted plot, for more details see Section 4). A second forced bar deposit is located on the left bank in the upstream part of the study reach, where a second plot of yellow-painted stones was deployed in May 2017 (hereinafter named PP2 painted plot, more details in Section 2) (Figure 3).
Additionally, a gauging station located at the downstream part of the study reach (Figure 3) provided water level and flow-discharge records during the whole study period (November 2016 to May 2019) (Figure 4).

4. Materials and Methods

4.1. Painted Plots

As outlined above, two plots of painted stones were deployed in the study site (Figure 3). The painted plots were deployed on the main body of two lateral gravel-bars. As far as our main aim with the plots was monitoring incipient motion conditions for coarse sediment particles, doing so with painted stones is easier in those areas of the channel that are not permanently under water, such as bars.
The first plot of tagged stones (PP1) covers a surface of 1 m × 1 m square on the riverbed. The stones were painted blue in situ, without moving or disturbing in any moment their position. Each grain was labelled and measured with a ruler. Additionally, we took several photographs to record the initial state of the painted stones. In May 2017, 100 stones were collected from the riverbed and taken to the lab, painted yellow, measured, and labelled. We seeded them on the riverbed. Several deployment options were considered for this second generation of tagged stones such as a square plot, a transverse line on the channel bed, or the random distribution of the stones on the streambed. Finally, we decided to seed the stones defining a square of tracers (1 m × 1 m) on the streambed (PP2) (Figure 3). Three were the main reasons: (1) to assure that all the stones were submitted to similar hydraulic conditions; (2) to keep coherence with the strategy previously followed with PP1; and (3) the identification in the field of incipient displacements and stone disturbance is easier with a painted plot
We made regular visits to the field (once every month, on average, and after each major peak of flow), in order to monitor the changes experienced by the painted plots. Each time the painted stones were displaced, we measured the grain size of the largest moved clast. To help this task, we compared photographs of the painted plots taken before and after the transport episode. We also measured the downstream transport distance travelled by the recovered clasts and registered the amount of painted stones recovered and lost after each flood.

4.2. Grain Size Analysis of Flood Deposits

In January 2019, a high-magnitude flow event occurred in the study river. Aerial photographs were available before and after the flood episode, thanks to a monitoring program with an unmanned aerial vehicle (UAV-drone platform). This helped the identification of a conspicuous gravel body deposited on the top of a gravel bar located in the study reach, undoubtedly related to the large flood. This sediment accumulation consists in an elongated gravel-sheet body overlapping a lateral bank deposit along the left margin of the river (Figure 5). Field visits were carried out to identify this deposit and to sample its GSD. We also sampled the GSD of the streambed surface and subsurface.
In the case of the streambed surface and the flood deposit, we used the Wolman (1954) [56] pebble count method. Each pebble count consisted of 200 grains collected along two ~50 m sampling lines spaced ~1–5 m apart. The sampling of grains was done systematically, extracting them at every 1 m intersection along a tape (around twice the largest grain size visually estimated in the field). To minimize the operator’s bias, all the grains were selected and measured by the same person. Metallic templates were used to measure the b-axis of grains > 8 mm and to classify them into half-Ψ size classes. Smaller grains were classified into two groups: grains between 4 and 8 mm and grains < 4 mm.
To sample the GSD of the subarmoured bed material, we used a variant of the pebble count method [56] adapted for the subsurface sediment [57,58]. A surface of 1.5 × 1.5 m was selected on the streambed, and the most superficial and coarse sediment layer was retired. We then mixed the underlying sediment using a shovel. Afterwards, a purpose-designed 10 × 10 cm wide sampling grid was deployed over the 1.5 × 1.5 m square surface. We measured the grain size of each pebble falling in each node of the grid. We repeated this operation at three different emplacements, measuring the size of the b-axis on ~100 subsurface particles per site. Finally, we integrated the three measures into a single GSD for the subarmour bed material. The sampled surface was 100 times larger than the largest clast, so we satisfied the Diplas and Fripp’s (1992) [59] criterion for a representative sample.

4.3. Modelling Streambed Mobility in the River Esva

Data from the painted plots were used with the purpose of defining the maximum mobile grain size at different discharges and to calibrate Equation (5) at low magnitude flows (size-selective conditions). Similarly, the GSD of the January 2019 flood deposit was used for the calibration of Equation (6) at high-magnitude flows (equal mobility conditions). The parameters β (in Equation (5)) and γ (in Equation (6)) were iteratively modified until obtaining those values that minimize the root square mean of differences between the observed (in the field) and the modelled competent sizes (Equation (5)) and GSD of the January 2019 flood deposit (in Equation (6)), respectively. Finally, Equation (7) was calibrated based on the value of γ computed for the January 2019 flood and the γ = 0 assigned to the May 2017 transport episode (close to the incipient motion conditions, see Section 5.3).
The resulting GTM model was applied to characterize streambed mobility in the River Esva. More specifically, we used the calibrated model to answer three questions: i. how the maximum and median competent sizes evolve with flow strength in the Esva river; ii. how the mobile fraction of fine particles (<8 mm) changes with flow stage; and iii. how the number of mobile particles grows with flow discharges. The maximum and median competent sizes were estimated using Equations (5)–(8). The fraction of bed load finer than 8 mm was estimated with Equation (6). The increase in the number of mobile particles was computed based on the mobile fraction of particles in each i-size range:
f s = φ i × p i p i
where fs is the portion of the streambed particles that are entrained and pi the relative fraction of each i-size class in the surface GSD.
Applying the GTM model requires the computation of shear stresses from flow discharge. To achieve this, we used the Rickenmann and Recking (2011) [60] fit to the Ferguson (2007) [61] flow resistance equation and estimated flow depths (d) from flow discharge using:
d = 0.015 × D 84 × Q *   2 p p 2.5
where:
Q * = Q / w × g S D 84 3
where Q is the water discharge, Q* the dimensionless water discharge, w the channel width, S the channel slope, and p = 0.24 if Q* < 100 and 0.31 otherwise. Then, assuming a rectangular cross-section, we estimated the hydraulic radius (R) from:
R = d × w ( 2 d + w )
and, finally, shear stresses were directly computed from the hydraulic radius–slope product:
τ = ρ × g × S × R
The GTM model also needs a value for the critical threshold of motion of the D84 (τc84*). We based it on the slope dependent relation proposed by Recking (2009) [62]:
τ c 84 * = ( 1.32 × S + 0.037 ) × ( D 84 D 50 ) 0.93
Finally, we computed the average GSD of the bed load mobilized during the entire study period (November 2016–July 2019). We used the 15 min discharge record available from the gauging station for the study period. Then, we applied the calibrated GTM model and Equation (8) to estimate the evolution of the GSD of the bed load at each 15 min time step of the flood hydrograph.
Then, the time integrated bed load volumes (Vi) in each size class, during the entire study period, can be estimated from:
V i = T = T 0 T = T f q s T · p p i T × T
where T0 and Tf are the departure and final times of the study period (respectively), qsT is the bed load rate at each time step T, and ppiT is the proportion of the i-size class in the bed load. Bed load rates were estimated using Recking’s (2013) [63] sediment transport equation, which appears to work well in coarse-bed rivers with the same range of slope and bed conditions as the River Esva [63,64]. Once we computed Vi for the different size classes, we calculated the percentage of each size class in the total bed load volumes (VT) mobilized during the study period:
% i = V i V T

4.4. Comparison to Previous Studies

We compared the GSD of the flood deposit sampled in the field to both the surface and subsurface GSD, in order to see how the grain size of the carried load relates to the grain size of the bed material. Later, in order to nourish the discussion on how the grain size of the carried load may or may not correlate to flow strength and the GSD of the streambed material, we re-analyzed previously published information on the grain size of the bed load available in Milhous (1973) [28], Andrews (1994) [29], Lisle (1995) [25] and Whittaker and Potts (2007) [1] (Table 1) and we compared them to our field observations.

5. Results

5.1. Painted Plots

Four floods were monitored with the painted stones from November 2016 to March 2019 (Table 2).
The first of the studied flood episodes occurred on 17 January 2017. It had a 41.93 m3/s peak discharge which corresponds to a frequent, regular peak of flow (<1 year return period discharge). This episode was able to move 23% of the seeded painted stones, but displacements were not very important, i.e., the maximum measured transport distances were relatively low (~0.2 m). The b-axis of the largest moved clast was 55 mm. These field observations can be considered as reliable, as approximately 85% of the tracers initially painted were retrieved. They reflect a partial mobility regime and close to the threshold of motion for the median-sized particles in the streambed.
The second transport episode, which occurred on 2 March 2017, was also monitored with the PP1 painted plot. It corresponded to a 74.63 m3/s discharge (<1 year return period flow). During this peak of discharge, a larger amount of painted tracers was disturbed (55%), but the displacements were again not very important; the maximum measured downstream transport distance was around 4 m. The size of the largest moved clast was 90 mm, coarser than in the previous episode. All these observations suggest partial mobility conditions and bed shear stresses close to the threshold of motion for the coarser particles in the riverbed. We recovered a large number of the initially painted stones (~82%), so all our observations were reliable.
The third transport episode took place on 29 May 2017. It was a 31.34 m3/s peak of discharge that was monitored with both the PP1 and PP2 painted plots. The PP1 painted stones did not experience any apparent movement. Conversely, the stones from the PP2 painted plot experienced some kind of reorganization, consisting in rolling and pivoting of some clasts and a very short downstream dispersion (<0.1 m). The PP2 was painted only six days before this event happened, while the PP1 painted stones had already experienced two previous peaks of flow. Consequently, differences in clast structuration may explain why the positions of PP2 stones were reorganized while those of PP1 clasts remained undisturbed.
The fourth and relative high-magnitude transport episode occurred on 12 November 2017 and was surveyed with both PP1 and PP2 painted stones. It had a peak discharge of 183.95 m3/s, which corresponds to a ~2 year return period flow. This episode displaced 86% of the painted stones from PP1, with some of the retrieved tracers experiencing downstream displacements as far as 35 m. The b-axis of the largest moved clast was 210 mm. We were able to recover 82% of the initially PP1 painted stones, so again the observations made with PP1 tracers could be considered as trustworthy. In contrast, we could only recover three of the initially painted stones from PP2; the rest were lost. The largest clast of the three PP2 retrieved stones measured 158 mm.

5.2. GSD of a Large Flood Deposit

An exceptionally large flood episode occurred in the River Esva between the 19th and 29th of January 2019. Peak discharge during this event reached 545 m3/s, which according to the available gauging records corresponds to a ~50 year return period flow discharge. A large sediment deposit was accumulated during this flood over a gravel bar in our study site (Figure 5). This deposit was first identified with the help of a UAV drone. Later, we recognized it in the field and we measured its GSD using the Wolman (1954) [56] method (Figure 6).
We also compared the GSD of the flood deposit to those of the surface and subsurface streambed material (Figure 6 and Table 3). The GSD of the flood deposit was very similar to that of the subsurface GSD, and both were considerably finer than the surface GSD. The similarity between the GSD of the January 2019 flood and the subsurface GSD suggests that the GSD of the bed material was mainly sculpted by large floods along the River Esva. Nevertheless, the flood deposit showed a larger percentage of sediment finer than 16 mm than the subsurface GSD. This difference might be related to winnowing of fine sediment during low floods, coarsening the topmost layer of the streambed and washing relatively fine particles out of the riverbed. However, potential bias when sampling GSD of the subsurface bed material cannot be discarded. Actually, when the three GSD are truncated at 8 mm, the GSD of the January 2019 flood deposit and the subsurface bed material seem to collapse into a very similar curve. This reinforces the idea that the subsurface sediment is depleted in fine sediment compared to the bed load carried during large floods.

5.3. Calibrating the GTM Model

The results from our field observations were used to calibrate the GTM competence model [2] for the River Esva. Painted plots allowed the calibration of Equation (5) which describes how the maximum mobile percentiles evolve with flow stage. We obtained an optimum ~3 − β parameter for transport stages below 1. In the case of transport stages larger than 1, we noticed how the painted plot data were in accordance with the 0.25 − β parameter initially suggested by Recking (2016) [2] that involves full mobility conditions when transport stages are larger than 2 (Figure 7A).
To calibrate Equation (6) (quantifying how mobile fraction in each size evolves with flow strength) at high flows (equal mobility conditions) we based on the GSD of the January 2019 flood deposit. We searched for the value of the γ parameter (in Equation (6)) that best minimizes the root mean square of differences between observed and modelled GSD (Figure 7B). We obtained a 2.58 − γ value. Nevertheless, according to Reference [2], the value of the γ parameter may vary with flow strength (Equation (7)). In order to better constrain our calibrated model, based on the data of 29 May peak of flow (when painted tracers practically did not move), we considered that the fraction of mobile particles in all size classes decreased to almost 0 at transport stages close to 0.6; this involves a γ parameter ~0. According to this and the ~2.6 γ value for transport stages close to 3.3, the best fit for Equation (7) was obtained with a 1.60 γ2 parameter (Figure 7C).

5.4. Characterizing Streambed Mobility in the River Esva

Using the GTM model calibrated for our study site, we accomplished an analysis on how streambed mobility evolves in the River Esva with flow stage. We estimated four parameters: (1) the median size of the entrained particles; (2) the size of the coarsest mobile clast; (3) the fraction of fine sediment (<8 mm) in the bed load at every flow stage; and (4) the number of bed surface particles that are entrained and actively incorporated into bed load motion with increasing discharges.
The modelled values clearly show how all these four parameters changed with flow strength (Figure 8). There was a remarkable change around ~75 m3/s where both the grain size and the number of entrained particles experienced a sudden jump. This discharge corresponded to a transport stage close to 1; that is to say, close to the moment when the coarser clasts in the riverbed started to incorporate into motion and the armour layer destabilized and broke up.

5.5. Modelling the Long-Term Averaged GSD of the Bed Load

Using the flow discharge records available for the study site, we estimated the GSD of the bed load integrated for the whole study period. This covers almost three entire hydrological years (2016–2017, 2017–2018, and 2018–2019), so we could consider it somewhat a proxy of the annual average GSD of the bed load. The obtained GSD is considerably finer than the GSD of the January 2019 flood deposit and the subsurface bed material (Figure 6). Indeed, this difference persists even when the GSDs are truncated at 8 mm (Figure 6B), so this is not just a consequence of a larger amount of fine sediment: the time integrated GSD of the bed load is still finer than the subsurface bed material, even when the fine sediment bed material load is not considered.

5.6. Comparison to Previous Studies

In Figure 9 and Figure 10 we plotted the River Esva data together with data from previous studies [1,25,28,29]. This comparison shows a statistically significant (p-value < 0.05) correlation between the largest clast size of the bed load and the transport stage ratio. The grain size of the largest transported clast approximates the diameter of the D84 when transport stages were close to 1 (Figure 9A and Figure 10A). This confirms how once the thresholds of motion of the coarser grains are reached, these coarse particles become incorporated into the bed load mass.
Conversely, in the case of the median sizes of the bed load, we did not observe any connection with the flow strength (Figure 9B and Figure 10B). Data from the different study sites appear well segregated which outlines the high degree of site specificity in the relations between GSD of the bed load and flow strength. This highlights the importance of site-specific controls on streambed mobility. Additionally, data plots corresponding to transport stages higher than 1 show an important amount of scatter. This brings into question the idea of a simple relation between the GSD of the subarmour material and that of the bed load during large floods.
Moreover, correlation between the size of the largest clast and flow strength is stronger when the largest moved clast sizes are normalized by the D84 of the surface sediment than when they are normalized using the D84 of the subsurface bed material. This may suggest how the remobilization of the coarser grains from the riverbed surface may largely control the size of the largest clast in motion.

6. Discussion

6.1. General Discussion

Our work illustrates the great potential and interest of combining a simple and open model such as the GTM [2] with relatively easy-to-take field measures in order to characterize the largest mobile size, fine sediment mobility or overall streambed mobility/stability behavior (Figure 8). Indeed, our field observations were in good accordance with the general shape of the GTM model [2] which is partially due to the large flexibility of the GTM model compared to previous competence models (Figure 11).
One main advantage of the GTM model relies is in its ability to model armour break-up and its adaptability to a broad range of armour destabilization conditions. The Esva calibrated GTM model shows two break points at ~75 and ~275 m3/s (Figure 11), corresponding to the flow discharges when the coarser grains started to be entrained and when full mobility conditions were reached, respectively. According to Figure 11, the GTM was able to model a situation where we have size-selective entrainment at low flows, a progressive armour destabilization at moderate flows and a sudden final breakup of the armour layer at very large flows. In fact, this involves some complex combination of size-selective and equal mobility conditions, very suitable for the River Esva (Figure 11). Previous competence models based on power-law hiding functions (Equation (3)) seem unable to grasp such behavior (Figure 11), as they are only able to model a progressive mobilization of the different grain sizes or a sudden mobilization of all the grain sizes at a narrow range of flows. The Komar (1987) [13] fit to Carling (1983) [15] data (involving a α-parameter ~0.85) was the one closest to the calibrated GTM model.
In this work, we did not consider the possibility of temporal changes in surface GSD when calibrating of the GTM model. However, changes in surface GSD could be expected in gravel-bed rivers, driven by variable flow discharges and temporal fluctuations in sediment supply [51,65]. Nevertheless, we accomplished several pebble counts (~100 counts per sample) in the study reach all along the time period of study (Figure 12). Between 2016 and the January 2019 flood, we did not observe significant grain size changes: we observed a ~15% variability around the mean value for the D50 which is comparable to the precision expected for pebble counts based on 100–200 clast-size measures [58,66]. Thus, we could consider that surface GSD remained quite stable between 2016 and 2019. However, after the January 2019 flood, we reported significant streambed fining. This was probably related to the liberation of fines after widespread armour disorganization during this high-magnitude flow.
Consequently, major floods in the River Esva seem to be associated with a temporal replenishment in fines of the streambed surface. Hence, we could expect differences in the GSD of the bed load after major floods, which will persist until the streambed exhausts again the fine sediment. Thus, we repeated the GTM model based on the surface GSD measured after the January 2019 flood and we compared to our previous estimations (Figure 13). In general, we observed how both the median grain size and the largest mobile clast during moderate flow discharges are finer after the January 2019 flood than before; conversely, differences tend to disappear for higher water discharges. On this point, seasonal and inter-annual changes in bed load rating curves have been documented by several authors in gravel-bed rivers [67,68]. We believe that the GTM model has the potential to grasp this kind of seasonal and annual bed load hysteresis patterns. Changes in in-channel sediment supply conditions should have an imprint on bed surface texture, grain size and streambed structuration [38,66]. Hence, if these changes in surface GSD are adequately monitored in the field, then they could be accounted for in the workflow here proposed (Figure 1) and the hysteresis behavior in flow competence could be approached based on the GTM model.

6.2. GSD of the Bed Load and the Streambed Sediment: Implications for Paleo-Hydrological Analysis

According to our results, the GSD of the time-integrated bed load is considerably finer than the GSD of the subsurface bed material. In addition, the comparison between the GSD of the January flood deposit and the GSD of the subsurface bed material suggests that the GSD of the subarmour material proxies that of the bed load carried during large floods in the River Esva. Church and Hassan (2005) [69] also compared the grain size of the bed load to that of the subsurface bed material in Harris Creek (Canada), based on bed load volumes collected with a pit trap. They also observed how the grain size of the bed load approaches that of the subsurface bed material in periods of large floods, and a time averaged bed load finer than the subsurface sediment.
Low-magnitude peaks of flows and recession limbs of major flows seem able to mobilize a large amount of fine sediment, which passes over the Esva’s riverbed and is not stocked into the bed material. Additionally, during large floods, much of the sand fraction of the bed material can escape as suspension load. Both issues help explain why the GSD of the time integrated bed load is finer than the GSD of the sub-armour sediment, although some potential methodological bias cannot be neglected. We defined the GSD of the time integrated bed load based on bed load rates computed using Recking’s (2013) [63] (updated by Recking et al. (2016) [70]) bed load equation, which is said to perform well in gravel-bed rivers [64]. However, at low flows, fine sediment exhaustion from the riverbed surface can give rise to a supply-limited situation, compromising the applicability of bed load formulae. Thus, the amount of fine sediment in the time integrated bed load may be overestimated.
In summary, the GSD of the bed load during large floods tends to resemble the GSD of the subarmour bed material. However, re-analysis of previously published data suggests that this is not a general trend in gravel-bed rivers, and it could be largely site specific. A wide scatter is present in the available data (Figure 9 and Figure 10) and the different data sets appear well-segregated, so site-specific conditions control the existing relations between the GSD of the subsurface bed material and that of bed load during large floods. Consequently, patterns of armour development and destabilization are largely site-specific, influencing the links between GSD of the bed load and that of the bed material.
In this regard, the size of the largest clast has been a classical tool in the evaluation of flow competence by many geologists and sedimentologists interested in the paleohydrological and paleohydraulical interpretation of old fluvial conglomerate deposits, who used it (for instance) for paleo-velocity or paleo-slope analysis [12,13,14,15,16,17,18,19,20,71,72,73,74,75]. However, translating flow competence measures, based on the largest clast size, into paleohydraulic conditions is not direct; as we have stated above, there is no single and general function linking the grain size of the bed load to flow strength. In addition, particle weight is not the only control on particle mobility, but packing, imbrication, clast interlocking, fine sediment content and a broad range of exposure and hiding effects distort the links between particle size and entrainment conditions. On this point, some authors documented equal mobility conditions for all the grain sizes present in the bed, whilst others have reported a much more size-selective behavior.
Moreover, gravel-bed rivers are systematically armored [48,49,50,51,52], i.e., the streambed surface is commonly coarser than the underlying subarmour deposit. In this regard, the topmost layer of the streambed represents the local source of sediment particles during bed load transport and surface coarsening results from the adjustment of the streambed’s surface to the bed load transfers during the dominant channel-forming flows. Thus, the differences between the surface and subsurface GSDs may potentially provide some interesting clues for evaluating streambed mobility and/or paleo-hydrological analysis. Unfortunately, the topmost layer of the streambed is rarely preserved in old fluvial deposits [75]. Consequently, the grain size of the subsurface sediment is the only available information for the stratigraphic application of flow competence. Additionally, grain size of fluvial deposits is probably driven mainly by hydraulic conditions during deposition [76,77] rather than by particle entrainment.
With the above in mind, some doubts arise about the stratigraphic interpretation of flow competences based on the largest grain measured in subarmour deposits, whether they may be a good proxy of large magnitude and less frequent transport episodes, or instead they provide an average picture of the more frequent channel-forming flows. Subsurface grain-size distribution has been commonly considered a good proxy for the grain size of the average annual bed load [53,78,79], but our results suggest that there is not a straightforward relation between the GSD of the subsurface material, the time integrated GSD of the bed load and/or that of the bed load during moderate to large floods. These relations are largely controlled by the specific streambed conditions (armoring/pavement, packing, imbrications, etc.) at a given river reach, so the paleo-hydrological interpretation of large clast measures in fluvial conglomerates could not be straightforward.

6.3. Frequency and Intensity of Streambed Mobility: Implications for Streambed Stability Analysis

The timing and intensity of riverbed mobility are of great interest in river ecology [6], since they determine the evolution of bed surface structure, which is the physical support of habitat for many aquatic organisms. Therefore, it is not surprising that stream ecologists have also used the flow competence concept in order to evaluate stream stability during flows [8,22,80,81] in relation to several issues, such as streambed perturbation [6] or the response of benthic macroinvertebrates to floods [82,83,84,85].
The magnitude and frequency of coarse particle motion contributes to shape habitat for macroinvertebrates that live on the riverbed surface and defines the bed conditions that render river gravel more stable through time between large floods. A reductionist view of gravel riverbed mobility based on the flow competence concept would suggest that the timing and intensity of streambed disturbance is driven by the ability of the channel to set in motion the coarser clasts. However, this overlooks the importance of other questions. For instance, the fact that a clast in a given class size is entrained does not mean that all clasts in that size are actually moving (partial mobility conditions [86]). We believe that the workflow outlined here is well suited for the characterization of streambed mobility in gravel-bed rivers, always bearing in mind the site-specific character of each riverbed.
Our results show an abrupt change in the streambed behavior at discharges able to destabilize the armored surface. The size of the largest moved clast and the median size of the bed load show a sudden increase. Furthermore, the percentage of fine sediment in the bed load decreases. This relates to a larger amount of coarser clasts set in motion that decreases the relative proportion of fine sediment. The percentage of mobile particles in the riverbed also increases, reaching a limiting value of ~80% during large floods as full mobility conditions are entirely never reached. This implies that an important number of clasts remain static during floods. This contrasts with observations made by Lisle et al. (2000) [87] and Haschenburger and Wilcock (2003) [88] who found that only a small portion of the bed is inactive during bankfull floods. Nevertheless, other authors have also observed that a large amount of the riverbed could be immobile during large floods [6,89,90], in agreement with our observations.

7. Conclusions

In this work, we combined field observations of particle entrainment using painted plots and grain size analysis of a large flood deposit, with the GTM model, recently developed [2], in order to quantify streambed mobility in the River Esva. The GTM model proved to be a very flexible tool for the characterization of patterns of particle entrainment and armour destabilization in our study site, and we illustrated the potential of this workflow to address these complex issues.
Additionally, we compared our results with field observations from previous studies that we compiled from the literature. We observed how the grain size of the bed load tended to be closer to that of the sub-armour bed material during large floods, while moderate magnitude flows tended to carry a relatively fine bed load. However, there was a huge scatter in the data. Hence, inferring the grain size of the bed load or paleo-flow conditions from the grain size analysis of flood deposits should be done with care.

Author Contributions

Conceptualization, D.V.-T.; methodology, D.V.-T., E.F.-I. and M.F.G.; formal analysis, D.V.-T.; investigation, D.V.-T. and E.F.-I.; resources, E.F.-I., M.F.G. and D.V.-T.; data curation, E.F.-I., M.F.G. and D.V.-T.; writing—original draft preparation, D.V.-T.; writing—review and editing, D.V.-T.; visualization, D.V.-T.; supervision, D.V.-T. and J.M.; project administration, J.M. and D.V.-T.; funding acquisition, J.M. and D.V.-T.

Funding

The present work was possible thanks to the financial support provided by the grant ACB17-44, co-funded by the post-doctoral “Clarín” program-FICYT (Government of the Principality of Asturias) and the Marie Curie Co-Fund. The authors acknowledge the support from the project RIVERCHANGES-CGL2015-68824-R (MINECO/FEDER, UE).

Acknowledgments

The authors would like to thank Daniel Grace and Hodei Uzkeda Apestegia for their review of the English of this manuscript. We are also thankful to the three anonymous reviewers that helped improve the final version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Whitaker, A.C.; Potts, D.F. Analysis of flow competence in an alluvial gravel bed stream, Dupuyer Creek, Montana. Water Resour. Res. 2007, 43, W07433. [Google Scholar] [CrossRef] [Green Version]
  2. Recking, A. A generalized threshold model for computing bed load grain-size distribution. Water Resour. Res. 2016, 52, 9274–9289. [Google Scholar] [CrossRef] [Green Version]
  3. Parker, G.; Klingeman, P.C. On why gravel bed streams are paved. Water Resour. Res. 1982, 18, 1409–1423. [Google Scholar] [CrossRef] [Green Version]
  4. Brayshaw, A.C. Bed microtopography and entrainment thresholds in gravel-bed rivers. GSA Bull. 1985, 96, 218–223. [Google Scholar] [CrossRef]
  5. Wilcock, P.R. Critical shear stress of natural sediment. J. Hydraul. Eng. 1993, 119, 491–505. [Google Scholar] [CrossRef]
  6. Haschenburger, J.K. Streambed Disturbance over a Long Flood Series. River Res. Appl. 2017, 33, 753–765. [Google Scholar] [CrossRef]
  7. Lancaster, J.; Belyea, L.R. Nested hierarchies and scale dependence of mechanisms of flow refugium use. J. N. Am. Benthol. Soc. 1997, 16, 221–238. [Google Scholar] [CrossRef]
  8. Townsend, C.R.; Scarsbrook, M.R.; Doledec, S. The intermediate disturbance hypothesis, refugia, and biodiversity in streams. Limnol. Oceanogr. 1997, 42, 939–949. [Google Scholar] [CrossRef] [Green Version]
  9. Matthaei, C.D.; Townsend, C.R. Long-Term effects of local disturbance history on mobile stream invertebrates. Oecologia 2000, 125, 119–126. [Google Scholar] [CrossRef]
  10. Pfeiffer, A.M.; Finnegan, N.J. Regional variation in gravel riverbed mobility, controlled by hydrologic regime and sediment supply. Geophys. Res. Lett. 2018, 45, 3097–3106. [Google Scholar] [CrossRef]
  11. Gilbert, G.K.; Murphy, E.C. The Transportation of Debris by Running Water; US Government Printing Office: Washington, DC, USA, 1914.
  12. Komar, P.D. The competence of turbidity current flow. Geol. Soc. Am. Bull. 1970, 81, 1555–1562. [Google Scholar] [CrossRef]
  13. Komar, P.D. Selective gravel entrainment and the empirical evaluation of flow competence. Sedimentology 1987, 34, 1165–1176. [Google Scholar] [CrossRef]
  14. Baker, V.R.; Ritter, D.F. Competence of rivers to transport coarse bedload material. Geol. Soc. Am. Bull. 1975, 86, 975–978. [Google Scholar] [CrossRef]
  15. Carling, P.A. Bedload transport in two gravel bedded streams. Earth Surf. Process. Landf. 1983, 14, 27–39. [Google Scholar] [CrossRef]
  16. Costa, J.E. Paleohydraulic reconstruction of flash-flood peaks from boulder deposits in the Colorado Front Range. Geol. Soc. Am. Bull. 1983, 94, 986–1004. [Google Scholar] [CrossRef]
  17. Williams, G.P. Paleohydrological methods and some Examples from Swedish fluvial environments. Geogr. Ann. 1983, 65, 227–243. [Google Scholar] [CrossRef]
  18. Komar, P.D.; Carling, P.A. Grain sorting in gravel-bed streams and the choice of particle sizes for flow-competence evaluations. Sedimentology 1991, 38, 489–502. [Google Scholar] [CrossRef]
  19. O’Connor, J.E. Hydrology, Hydraulics and Geomorphology of the Bonneville Flood; Geological Society of America: Boulder, CO, USA, 1991. [Google Scholar]
  20. Bradley, W.C.; Mears, A.I. Calculations of flows needed to transport coarse fraction of Boulder Creek alluvium at Boulder, Colorado. Geol. Soc. Am. Bull. 1993, 86, 1057–1090. [Google Scholar]
  21. Wilcock, P.R. Flow competence: A criticism of a classic concept. Earth Surf. Process. Landf. 1992, 17, 289–298. [Google Scholar] [CrossRef]
  22. Lorang, M.; Hauer, F. Flow competence and streambed stability: An evaluation of technique and application. J. N. Am. Benthol. Soc. 2003, 22, 475–491. [Google Scholar] [CrossRef]
  23. Parker, G. Transport of gravel and sediment mixtures. In Sedimentation Engineering: Processes, Measurements, Modeling, and Practice. The ASCE Manuals and Reports in Engineering Practice, No. 110; García, M., Ed.; American Society of Civil Engineers: Reston, VA, USA, 2008; pp. 165–264. [Google Scholar]
  24. Juez, C.; Hassan, M.A.; Franca, M.J. The origin of fine sediment determines the observations of suspended sediment fluxes under unsteady flow conditions. Water Resour. Res. 2018, 54, 1–16. [Google Scholar] [CrossRef]
  25. Lisle, T.E. Particle size variations between bed load and bed material in natural gravel bed channels. Water Resour. Res. 1995, 31, 1107–1118. [Google Scholar] [CrossRef]
  26. Parker, G. Surface-Based bedload transport relation for gravel rivers. J. Hydraul. Res. 1990, 28, 417–428. [Google Scholar] [CrossRef]
  27. Wilcock, P.R.; Crowe, J.C. Surface-Based transport model for mixed-size sediment. J. Hydraul. Eng. 2003, 129, 120–128. [Google Scholar] [CrossRef]
  28. Milhous, R.T. Sediment Transport in a Gravel Bottom Stream. Ph.D. Thesis, Oregon State University, Corvallis, OR, USA, 1973. [Google Scholar]
  29. Andrews, E.D. Marginal bed load transport in a gravel bed stream, Sagehen Creek, California. Water Resour. Res. 1994, 30, 2241–2250. [Google Scholar] [CrossRef]
  30. Shields, A. Application of Similarity Principles and Turbulence Research to Bedload Movement. Ph.D. Theisis, California Institute of Technology, Pasadena, CA, USA, 1936. Originally published in German in 1936, English translation by W.P. Ott and J.C. van Uchelen, available in Hydrodynamics Laboratory Publication No. 167. Hydrodynamics Laboratory. [Google Scholar]
  31. Hjülstrom, F. Studies of morphological activity of rivers as illustrated by the River Fyris, Upsala Mineral. Geologis. Inst. Bull. 1935, 25, 221–527. [Google Scholar]
  32. Sundborg, A. The river Klarälven a study of fluvial processes. Geogr. Ann. 1956, 38, 125–316. [Google Scholar]
  33. Buffington, J.M.; Montgomery, D.R. A systematic analysis of eight decades of incipient motion studies, with special reference to gravel-bedded rivers. Water Resour. Res. 1997, 33, 1993–2027. [Google Scholar] [CrossRef] [Green Version]
  34. Parker, G. 1D Sediment Transport Morphodynamics with Applications to Rivers and Turbidity Currents. 2004. Available online: http://hydrolab.illinois.edu/people/parkerg//morphodynamics_e-book.htm (accessed on 26 February 2019).
  35. Einstein, H.A.; Chien, N. Transport of Sediment Mixtures with Large Range of Grain Size; University of California, Institute of Engineering Research: Oakland, CA, USA, 1953. [Google Scholar]
  36. Egiazaroff, I.V. Calculation of nonuniform sediment concentrations. J. Hydraul. Div. 1965, 91, 225–247. [Google Scholar]
  37. Wilcock, P.R. Toward a practical method for estimating sediment transport rates in gravel-bed rivers. Earth Surf. Process. Landf. 2001, 26, 1395–1408. [Google Scholar] [CrossRef]
  38. Church, M.; Hassan, M.A.; Wolcott, J.F. Stabilizing self-organized structures in gravel-bed stream channels: Field and experimental observations. Water Resour. Res. 1998, 34, 3169–3179. [Google Scholar] [CrossRef]
  39. Ashworth, P.J.; Ferguson, R.I. Size-selective entrainment of bed load in gravel bed streams. Water Resour. Res. 1989, 25, 627–634. [Google Scholar] [CrossRef]
  40. Lenzi, M.A.; D’Agostino, V.; Billi, P. Bedload transport in the instrumented catchment of the Rio Cordon, part I: Analysis of bedload records, conditions and threshold of bedload entrainment. Catena 1989, 36, 171–190. [Google Scholar] [CrossRef]
  41. Powell, D.M.; Reid, I.; Laronne, J.B. Evolution of bed load grain-size distribution with increasing flow strength and the effect of flow duration on the caliber of bed load sediment yield in ephemeral gravel bed rivers. Water Resour. Res. 2001, 37, 1463–1474. [Google Scholar] [CrossRef]
  42. Parker, G.; Klingeman, P.C.; McLean, D.G. Bed load and size distribution in gravel-bed streams. J. Hydraul. Div. Am. Soc. Civ. Eng. 1982, 108, 544–571. [Google Scholar]
  43. Andrews, E.D. Entrainment of gravel from naturally sorted riverbed material. Geol. Soc. Am. Bull. 1983, 94, 1225–1231. [Google Scholar] [CrossRef]
  44. Wilcock, P.R.; Southard, J.B. Experimental study of incipient motion in mixed-size sediment. Water Resour. Res. 1988, 24, 1137–1151. [Google Scholar] [CrossRef] [Green Version]
  45. Andrews, E.D.; Smith, J.D. A theoretical model for calculating marginal bedload transport rates of gravel. In Dynamic of Gravel Bed Rivers; Billi, P., et al., Eds.; John Wiley: Chichester, UK, 1992; pp. 41–52. [Google Scholar]
  46. Wilcock, P.R.; McArdell, B.W. Partial transport of a sand/gravel sediment. Water Resour. Res. 1997, 33, 235. [Google Scholar] [CrossRef] [Green Version]
  47. Proffitt, G.T.; Sutherland, A.J. Transport of non uniform sediments. J. Hydraul. Res. 1983, 21, 33–43. [Google Scholar] [CrossRef]
  48. Gomez, B. Temporal variations in bedload transport rates: The effect of progressive bed armoring. Earth Surf. Process. Landf. 1983, 8, 41–54. [Google Scholar] [CrossRef]
  49. Parker, G.; Sutherland, A.J. Fluvial armor. J. Hydraul. Res. 1990, 28, 529–544. [Google Scholar] [CrossRef]
  50. Richards, K.; Clifford, N.J. Fluvial geomorphology: Structured beds in gravelly rivers. Prog. Phys. Geogr. 1991, 15, 407–422. [Google Scholar] [CrossRef]
  51. Hassan, M.A.; Egozi, R.; Parker, G. Experiments on the effect of hydrograph characteristics on vertical grain sorting in gravel bed rivers. Water Resour. Res. 2006, 42, W09408. [Google Scholar] [CrossRef]
  52. Venditti, J.G.; Nelson, P.A.; Bradley, R.W.; Haught, D.; Gitto, A.B. Bedforms, structures, patches, and sediment supply in gravel-bed rivers. In Gravel-Bed Rivers: Processes and Disasters; Tsutsumi, D., Laronne, J.B., Eds.; John Wiley: Chichester, UK, 2017; pp. 439–466. [Google Scholar]
  53. Recking, A. An analysis of nonlinearity effects on bed load transport prediction. J. Geophys. Res. Earth Surf. 2013, 118, 1264–1281. [Google Scholar] [CrossRef] [Green Version]
  54. MacKenzie, L.G.; Eaton, B.C.; Church, M. Breaking from the average: Why large grains matter in gravel-bed streams. Earth Surf. Process. Landf. 2018, 43, 3190–3196. [Google Scholar] [CrossRef]
  55. Vázquez-Tarrío, D.; Tal, M.; Camenen, B.; Piegay, H. Effects of continuous embankments and sucessive run-of-the-river dams on bedload transport capacities along the Rhône River, France. Sci. Total. Environ. (STOTEN) 2019, 658, 1375–1389. [Google Scholar] [CrossRef] [PubMed]
  56. Wolman, M. A method for sampling coarse river-bed material. Trans. Am. Geophys. Union 1954, 35, 951–956. [Google Scholar] [CrossRef]
  57. Buffington, J.M. An alternative method for determining subsurface grain-size distributions of gravel-bedded rivers, Supplement to EOS. AGU Trans. 1996, 77, 250. [Google Scholar]
  58. Bunte, K.; Abt, S.R. Sampling Surface and Subsurface Particle-Size Distributions in Wadable Gravel-and Cobble-Bed Streams for Analyses in Sediment Transport, Hydraulics, and Streambed Monitoring; Gen. Tech. Rep. RMRS-GTR-74; U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2001; p. 428.
  59. Diplas, P.; Fripp, J.B. Properties of various sediment sampling procedures. J. Hydraul. Eng. 1992, 118, 955–970. [Google Scholar] [CrossRef]
  60. Rickenmann, D.; Recking, A. Evaluation of flow resistance in gravel-bed rivers through a large field data set. Water Resour. Res. 2011, 47, W07538. [Google Scholar] [CrossRef] [Green Version]
  61. Ferguson, R. Flow resistance equations for gravel and boulder bed streams. Water Resour. Res. 2007, 43, W05427. [Google Scholar] [CrossRef] [Green Version]
  62. Recking, A. Theoretical development on the effects of changing flow hydraulics on incipient bed load motion. Water Resour. Res. 2009, 45, W04401. [Google Scholar] [CrossRef]
  63. Recking, A. Simple method for calculating reach-averaged bed-load transport. J. Hydraul. Eng. 2013, 139, 70–75. [Google Scholar] [CrossRef]
  64. Hinton, D.; Hotchkiss, R.H.; Cope, M. Comparison of Calibrated Empirical and Semi-Empirical Methods for Bedload Transport Rate Prediction in Gravel Bed Streams. J. Hydraul. Eng. 2018, 144, 04018038. [Google Scholar] [CrossRef]
  65. Plumb, B.D.; Juez, C.; Annable, W.K.; McKie, C.W.; Franca, M.J. The impact of hydrograph variability and frequency on sediment transport dynamics in a gravel-bed flume. Earth Surf. Process. Landf. 2019. [Google Scholar] [CrossRef]
  66. Rice, S.; Church, M. Sampling surficial fluvial gravels: the precision of size distribution percentile estimates. J. Sediment. Res. 1996, 66, 654–665. [Google Scholar] [CrossRef]
  67. Moog, D.B.; Whiting, P.J. Annual hysteresis in bed load rating curves. Water Resour. Res. 1998, 34, 2393–2399. [Google Scholar] [CrossRef]
  68. Hassan, M.A.; Church, M. Sensitivity of bed load transport in Harris Creek: Seasonal and spatial variation over a cobble-gravel bar. Water Resour. Res. 2001, 37, 813–825. [Google Scholar] [CrossRef]
  69. Church, M.; Hassan, M. Upland gravel-bed rivers with low sediment transport. In Developments in Earth Surface Processes 7; Garcia, C., Batalla, R.J., Eds.; Elsevier: Amsterdam, The Netherlands, 2005; Volume 7, pp. 141–168. [Google Scholar]
  70. Recking, A.; Piton, G.; Vázquez-Tarrío, D.; Parker, G. Quantifying the morphological print of bedload transport. Earth Surf. Process. Landf. 2016, 41, 809–822. [Google Scholar] [CrossRef]
  71. Reid, I.; Frostick, L.E. Particle interaction and its effect on the thresholds of initial and final bedload motion in coarse alluvial channels. In Sedimentology of Gravels and Conglomerates; Koster, E.H., Steel, R.J., Eds.; Canadian Society of Petroleum Geologists: Calgary, AB, Canada, 1984; pp. 61–68. [Google Scholar]
  72. Church, M.; Wolcott, J.; Maizels, J.K. Paleovelocity: A parsimonious proposal. Earth Surf. Process. Landf. 1990, 15, 475–480. [Google Scholar] [CrossRef]
  73. Rotnicki, A. Retroduction and palaeodischarges of meandering and sinuous alluvial rivers and its palaeohydroclimatic implications. In Temperate Palaeohydrology; Starkel, L., Gregory, K.J., Thorne, J.B., Eds.; John Wiley: Chichester, UK, 1991; pp. 431–471. [Google Scholar]
  74. Paola, C.; Morhrig, D. Paleohydraulics revisited: Paleoslope estimation in coarse grained braided rivers. Basin Res. 1996, 8, 243–254. [Google Scholar] [CrossRef]
  75. Jones, S.J.; Frostick, L.E. Inferring bedload transport from stratigraphic successions: Examples from cenozoic and pleistocene rivers, south central pyrenees, spain. In Landscapes Evolution: Denudation, Climate and Tectonics over Different Time and Space Scales; Gallagher, K., Jones, S.J., Wainwright, J., Eds.; Special Publications of the Geological Society: London, UK, 2008; pp. 129–145. [Google Scholar]
  76. Reid, I.; Frostick, L.E. Role of settling, entrainment and dispersive equivalence and of interstice trapping in placer formation. J. Geol. Soc. 1985, 142, 739–746. [Google Scholar] [CrossRef]
  77. Pähtz, T.; Durán, O. The cessation threshold of nonsuspended sediment transport across aeolian and fluvial environments. J. Geophys. Res. Earth Surf. 2018, 123, 1638–1666. [Google Scholar]
  78. Kuhnle, R.A.; Willis, J.C. Mean size distribution of bed load on Goodwin Creek. J. Hydraul. Eng. 1992, 118, 1443–1446. [Google Scholar] [CrossRef]
  79. Segura, C.; Pitlick, J. Coupling fluvial-hydraulic models to predict gravel transport in spatially variable flows. J. Geophys. Res. Earth Surf. 2015, 120, 834–855. [Google Scholar] [CrossRef]
  80. Downes, B.J.; Glaister, A.; Lake, P.S. Spatial variation in the force required to initiate rock movement in 4 upland streams: Implications for estimating disturbance frequencies. J. N. Am. Benthol. Soc. 1997, 16, 203–220. [Google Scholar] [CrossRef]
  81. Duncan, M.J.; Suren, A.M.; Brown, S.L.R. Assessment of streambed stability in steep, bouldery streams: Development of a new analytical technique. J. N. Am. Benthol. Soc. 1999, 18, 445–456. [Google Scholar] [CrossRef]
  82. Scrimgeour, G.J.; Davidson, R.J.; Davidson, J.M. Recovery of benthic macroinvertebrate and epilithic communities following a large flood, in an unstable, Braided, New Zealand River. N. Z. J. Mar. Freshw. Res. 1998, 22, 337–344. [Google Scholar] [CrossRef]
  83. Robinson, C.T.; Aebischer, S.; Uehlinger, U. Immediate and habitat-specific responses of macroinvertebrates to sequential, experimental floods. J. N. Am. Benthol. Soc. 2004, 23, 853–867. [Google Scholar] [CrossRef]
  84. Gibbins, C.; Vericat, D.; Batalla, R.J. When is stream invertebrate drift catastrophic? The role of hydraulics and sediment transport in initiating drift during flood events. Freshw. Biol. 2007, 52, 2369–2384. [Google Scholar] [CrossRef]
  85. Power, M.E.; Parker, M.S.; Dietrich, W.E. Seasonal reassembly of a river food web: Floods, droughts, and impacts of fish. Ecol. Monogr. 2008, 78, 263–282. [Google Scholar] [CrossRef] [Green Version]
  86. Wilcock, P.R.; Mc Ardell, B.W. Surface-Based fractional transport rates: Mobilization thresholds and partial transport of a sand-gravel sediment. Water Resour. Res. 1993, 29, 1297–1312. [Google Scholar] [CrossRef]
  87. Lisle, T.E.; Nelson, J.M.; Pitlick, J.; Madej, M.A.; Barkett, B.L. Variability of bed mobility in natural gravel-bed channels and adjustments to sediment load at local and reach scales. Water Resour. Res. 2000, 36, 3743–3755. [Google Scholar] [CrossRef]
  88. Haschenburger, J.K.; Wilcock, P.R. Partial transport in a natural gravel bed channel. Water Resour. Res. 2003, 39, 1020. [Google Scholar] [CrossRef]
  89. May, C.L.; Pryor, B.; Lisle, T.E.; Lang, M. Coupling hydrodynamic modeling and empirical measures of bed mobility to predict the risk of scour and fill of salmon redds in a large regulated river. Water Resour. Res. 2009, 45. [Google Scholar] [CrossRef] [Green Version]
  90. McKean, J.; Tonina, D. Bed stability in unconfined gravel bed mountain streams: With implications for salmon spawning viability in future climates. J. Geophys. Res. Earth Surf. 2013, 118, 1227–1240. [Google Scholar] [CrossRef]
Figure 1. Workflow followed in the present research with the aim of quantitatively characterizing streambed mobility in the River Esva.
Figure 1. Workflow followed in the present research with the aim of quantitatively characterizing streambed mobility in the River Esva.
Water 11 02662 g001
Figure 2. Location of the study river in NW Spain.
Figure 2. Location of the study river in NW Spain.
Water 11 02662 g002
Figure 3. Aerial image of the study reach and location of the two monitored painted plots.
Figure 3. Aerial image of the study reach and location of the two monitored painted plots.
Water 11 02662 g003
Figure 4. Flow hydrograph of the River Esva during the study period. Vertical dashed lines represent dates of field surveys.
Figure 4. Flow hydrograph of the River Esva during the study period. Vertical dashed lines represent dates of field surveys.
Water 11 02662 g004
Figure 5. (A,B) UAV-aerial image of the River Esva before (A) and after (B) January 2019 flood. (C,D) gravel sheet body studied and sampled in the field, days after the flood. Blue arrows indicates flow direction
Figure 5. (A,B) UAV-aerial image of the River Esva before (A) and after (B) January 2019 flood. (C,D) gravel sheet body studied and sampled in the field, days after the flood. Blue arrows indicates flow direction
Water 11 02662 g005
Figure 6. The grain-size distribution (GSD) curves for the River Esva: (A) non-truncated and (B) truncated at 8 mm.
Figure 6. The grain-size distribution (GSD) curves for the River Esva: (A) non-truncated and (B) truncated at 8 mm.
Water 11 02662 g006
Figure 7. Results of the GTM model calibration for the River Esva. (A). Calibrating Equation (5) using data from painted plots. (B): Calibrating Equation (6) using the GSD of the January 2019 flood deposit. (C): Calibrating Equation (7) using data derived from painted plots and the GSD of the January 2019 flood deposit.
Figure 7. Results of the GTM model calibration for the River Esva. (A). Calibrating Equation (5) using data from painted plots. (B): Calibrating Equation (6) using the GSD of the January 2019 flood deposit. (C): Calibrating Equation (7) using data derived from painted plots and the GSD of the January 2019 flood deposit.
Water 11 02662 g007
Figure 8. Streambed mobility model for the River Esva, derived from the calibrated GTM model.
Figure 8. Streambed mobility model for the River Esva, derived from the calibrated GTM model.
Water 11 02662 g008
Figure 9. Comparison of the results obtained for the River Esva with those of previous studies focused on flow competence. In this figure, we compare the grain size of the bed load to that of the streambed surface, accounting for flow strength. (A) Maximum mobile size compared to surface D84. (B) Median mobile size compared to surface D50. References: Oak creek [28]; Sagehen creek [29]; Dupuyer creek [1]; Other [25].
Figure 9. Comparison of the results obtained for the River Esva with those of previous studies focused on flow competence. In this figure, we compare the grain size of the bed load to that of the streambed surface, accounting for flow strength. (A) Maximum mobile size compared to surface D84. (B) Median mobile size compared to surface D50. References: Oak creek [28]; Sagehen creek [29]; Dupuyer creek [1]; Other [25].
Water 11 02662 g009
Figure 10. Comparison of the results obtained for the River Esva with those of previous studies focused on flow competence. In this figure, we compare the grain size of the bed load to that of the streambed subsurface, accounting for flow strength. (A) Maximum mobile size compared to subsurface D84. (B) Median mobile size compared to subsurface D50. References: Oak creek [28]; Sagehen creek [29]; Dupuyer creek [1]; Other [25].
Figure 10. Comparison of the results obtained for the River Esva with those of previous studies focused on flow competence. In this figure, we compare the grain size of the bed load to that of the streambed subsurface, accounting for flow strength. (A) Maximum mobile size compared to subsurface D84. (B) Median mobile size compared to subsurface D50. References: Oak creek [28]; Sagehen creek [29]; Dupuyer creek [1]; Other [25].
Water 11 02662 g010
Figure 11. Comparison between the GTM model (calibrated for the River Esva) and several other flow competence models applied to the River Esva. In general, previous competence models simulated either a progressive increase in the size of the mobile clast with flow discharge (size-selective behavior, e.g., Komar [13], Andrews [42], Shields [30]) or a sudden incorporation of all the particle sizes at a narrow range of discharges (equal mobility, e.g., Parker et al. [41] and Recking [62]). However, the GTM model was able to simulate a combination of size-selective behaviors at low flows and equal mobility once the armour layer destabilized entirely which seems to agree better with our field observations in the River Esva.
Figure 11. Comparison between the GTM model (calibrated for the River Esva) and several other flow competence models applied to the River Esva. In general, previous competence models simulated either a progressive increase in the size of the mobile clast with flow discharge (size-selective behavior, e.g., Komar [13], Andrews [42], Shields [30]) or a sudden incorporation of all the particle sizes at a narrow range of discharges (equal mobility, e.g., Parker et al. [41] and Recking [62]). However, the GTM model was able to simulate a combination of size-selective behaviors at low flows and equal mobility once the armour layer destabilized entirely which seems to agree better with our field observations in the River Esva.
Water 11 02662 g011
Figure 12. Comparison among the different pebble counts accomplished in the study site during the time period spanned by the present research (2016–2019). Calibration of the GTM model was based in the pebble count accomplished in September 2018.
Figure 12. Comparison among the different pebble counts accomplished in the study site during the time period spanned by the present research (2016–2019). Calibration of the GTM model was based in the pebble count accomplished in September 2018.
Water 11 02662 g012
Figure 13. Comparison of the GTM model estimations for the River Esva before and after the January 2019 flood.
Figure 13. Comparison of the GTM model estimations for the River Esva before and after the January 2019 flood.
Water 11 02662 g013
Table 1. Data of grain size of the bed load, compiled from previous studies and compared to our own data for the River Esva. S: Slope. W: Channel width (in m). D50s and D50ss: Median size of the surface and subsurface grain-size distribution (in mm), respectively. D50bl: Median size of the bed load. Dmax: Largest mobile grain (in mm). Qbkf: Water discharge at bankfull (m3/s). Q: Water discharge (m3/s) during the sampled sediment transport episodes. Sample: Bed load sediment transport method.
Table 1. Data of grain size of the bed load, compiled from previous studies and compared to our own data for the River Esva. S: Slope. W: Channel width (in m). D50s and D50ss: Median size of the surface and subsurface grain-size distribution (in mm), respectively. D50bl: Median size of the bed load. Dmax: Largest mobile grain (in mm). Qbkf: Water discharge at bankfull (m3/s). Q: Water discharge (m3/s) during the sampled sediment transport episodes. Sample: Bed load sediment transport method.
RiverSourceSWD50sD50ssD50blDmaxQbkfQSample
Oak creek[28]0.00943.754203–278–983.40.2–3.4Vortex sampler
Sagehen creek[29]0.01022.6583012–6531–1260.30.2–0.6HS
Bambi creek[25]0.00823.6-93221.72.5HS
Goodwin creek[25]0.003312.912882490.081HS
Jacoby creek[25]0.006317.22784259.034.2HS
NF Caspar creek[25]0.01304.41594803.16.2Pit
Redwood creek 1[25]0.014044.654130430.0374.1HS
Tanana river[25]0.0008315.01311722-2040HS
Tom Mac Donald Creek[25]0.00606.119.8118223.612.6HS
Turkey brook[25]0.00863.022161110013.016.9Pit
Dupuyer creek[1]0.01008.054-47–6088–1556.54.6–10.3Basket sampler
River EsvaThis study0.003025.25922-16–208-31.3–183.9Painted stones
Table 2. Summary of results of the field surveys in the River Esva accomplished here.
Table 2. Summary of results of the field surveys in the River Esva accomplished here.
DatePainted PlotQ (m3/s)τ*c*Dmax% RecoveryL (m)
17/01/2017PP141.930.7354.586~0.2
02/03/2017PP174.630.9890824.0
29/05/2017PP131.340.630820.0
12/11/2017PP1183.951.552088235.2
29/05/2017PP231.340.631696<0.1
12/11/2017PP2183.951.55158319
Q: Peak flow discharge. τ*c*: Transport stage ratio (at peak discharge). Dmax: Maximum mobile size (mm). L: Maximum distance travelled by the retrieved painted stones.
Table 3. Comparison between the grain-size distribution (GSD) of the surface sediment, the subarmour bed material, and the January 2019 flood deposit in the River Esva. Values between brackets correspond to the grain size metrics estimated for the truncated (<8 mm) GSD.
Table 3. Comparison between the grain-size distribution (GSD) of the surface sediment, the subarmour bed material, and the January 2019 flood deposit in the River Esva. Values between brackets correspond to the grain size metrics estimated for the truncated (<8 mm) GSD.
GSD% Particles < 8 mmD16 (mm)D50 (mm)D84 (mm)
Surface034 (34)59 (59)107 (107)
Subsurface811 (20)22 (30)47 (58)
January 2019
(flood deposit)
223 (20)22 (30)42 (49)

Share and Cite

MDPI and ACS Style

Vázquez-Tarrío, D.; Fernández-Iglesias, E.; Fernández García, M.; Marquínez, J. Quantifying the Variability in Flow Competence and Streambed Mobility with Water Discharge in a Gravel-Bed Channel: River Esva, NW Spain. Water 2019, 11, 2662. https://doi.org/10.3390/w11122662

AMA Style

Vázquez-Tarrío D, Fernández-Iglesias E, Fernández García M, Marquínez J. Quantifying the Variability in Flow Competence and Streambed Mobility with Water Discharge in a Gravel-Bed Channel: River Esva, NW Spain. Water. 2019; 11(12):2662. https://doi.org/10.3390/w11122662

Chicago/Turabian Style

Vázquez-Tarrío, Daniel, Elena Fernández-Iglesias, María Fernández García, and Jorge Marquínez. 2019. "Quantifying the Variability in Flow Competence and Streambed Mobility with Water Discharge in a Gravel-Bed Channel: River Esva, NW Spain" Water 11, no. 12: 2662. https://doi.org/10.3390/w11122662

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

Article Metrics

Back to TopTop