Next Article in Journal
Groundwater Arsenic-Attributable Cardiovascular Disease (CVD) Mortality Risks in India
Next Article in Special Issue
Field Study on Wood Accumulation at a Bridge Pier
Previous Article in Journal
Comparative Assessment of Methods for Coupling Regional and Local Groundwater Flow Models: A Case Study in the Beijing Plain, China
Previous Article in Special Issue
Flume Experiments Evaluating the Efficacy of a Large Wood Trap Featuring Horizontal Rods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Wood Retention at Inclined Bar Screens: Effect of Wood Characteristics on Backwater Rise and Bedload Transport

by
Isabella Schalko
1,*,
Virginia Ruiz-Villanueva
2,
Fiona Maager
1 and
Volker Weitbrecht
1
1
ETH Zurich, Laboratory of Hydraulics, Hydrology and Glaciology (VAW), 8093 Zurich, Switzerland
2
Institute of Earth Surface Dynamics (IDYST), Faculty of Geosciences and the Environment, University of Lausanne, 1015 Lausanne, Switzerland
*
Author to whom correspondence should be addressed.
Water 2021, 13(16), 2231; https://doi.org/10.3390/w13162231
Submission received: 29 June 2021 / Revised: 11 August 2021 / Accepted: 12 August 2021 / Published: 16 August 2021
(This article belongs to the Special Issue Impact of Large Wood on River Ecosystems)

Abstract

:
In forested mountain catchment areas, both bedload and large wood (LW) can be transported during ordinary flows. Retention structures such as sediment traps or racks are built to mitigate potential hazards downstream. Up to now, the design of these retention structures focuses on either LW or bedload. In addition, the majority of LW retention racks tend to retain both LW and bedload, while bedload transport continuity during ordinary flows is an important aspect to be considered in the design. Therefore, a series of flume experiments was conducted to study the effect of LW accumulations at an inclined bar screen with a bottom clearance on backwater rise and bedload transport. The main focus was put on testing different LW characteristics such as LW size, density, fine material, and shape (branches and rootwads), as well as a sequenced flood. The results demonstrated that a few logs (wood volume of ≈ 7 m 3 prototype scale with a model scale factor of 30) are sufficient to reduce the bedload transport capacity to below 75% compared to the condition without LW. Fine material and smaller wood sizes further reduced bedload transport and increased backwater rise. In contrast, LW density and LW shape had a negligible effect. The test focusing on a sequenced flood highlighted the need for maintenance measures to avoid self-flushing of the bed material. The results of this study further indicate that an inclined bar screen may need to be adapted by considering LW characteristics in the design of the bottom clearance to enable bedload continuity during ordinary flows.

1. Introduction

Potential hazards associated with intense bedload transport in steep mountain streams have traditionally been managed by building check dams and sediment trap structures [1]. The design of such retention structures has evolved through time, and nowadays a variety of models exists, ranging from solid dams to slit and other types of open check dams [2,3]. The design depends on both the purpose and function of the structure, which can be grouped in (i) stabilization and consolidation to prevent or stop bed incision or lateral erosion, (ii) retention and storage of water and/or sediment, (iii) dosing and filtering to regulate peak flows and to sort selected sediment fractions, and (iv) energy dissipation to reduce stream power and break debris flows [1,4]. An important aspect for designing sediment traps is related to the impact on the stream ecosystem. By retaining flow and sediment, the natural regimes downstream of the structure are altered, leading to morphological and physical changes, such as incision, lateral erosion, or armoring, with important consequences in terms of habitat loss and changes in ecological characteristics and biodiversity [5,6,7]. Therefore, it is important to design retention structures that are effective in retaining sediment during potentially hazardous events (high flows), but are permeable during ordinary flows. With this aim, a novel design of semipermeable check dams has recently been proposed by combining a hydraulic with a mechanical barrier and a guiding channel [8,9,10,11]. This semipermeable structure proved to be successful in terms of maintaining bedload transport continuity during ordinary flows and bedload retention during higher flows. However, mountain streams in forested catchments can also transport large wood (LW; herein defined as logs with a length ≥1 m and a diameter ≥0.1 m [12,13,14]), such as downed trees, trunks, rootwads, and branches, during both ordinary and high flows. The choice of retention structure type should, therefore, also consider the transport and retention of LW [15,16]. Design guidelines for retention structures in mountain streams that properly manage both bedload and LW are still lacking. A novel rack design for low-gradient (small Froude number) river sections is the bypass retention system. This rack is installed at an outer river bend parallel to the flow and enables the separation of bedload and LW due to the secondary currents induced by the river curvature, thereby allowing the structure to trap LW without limiting the bedload transport continuity [17]. As this rack design requires both smaller slope gradients and a river bend, it may not be efficient in steeper, mountain streams. The majority of sediment traps has been designed without considering LW or the conducted experiments were limited to a specific case study [11,18,19]. In addition, common wood retention structures are usually designed to only retain LW, whereas they tend to also retain bedload [20]. Flume experiments with inclined racks showed that both backwater rise and local scour decreases with increasing rack inclination [21]. First experiments at an inclined bar screen with a lower gap demonstrated that LW volume blocking 20% of the flow cross section is sufficient to reduce bedload transport by 50% [22] and the effect of LW on bedload transport should, thus, not be overlooked.
Based on a preliminary proof-of-concept by Schalko and Weitbrecht [22], this paper aims at testing the effects of LW accumulation at inclined bar screens on the bedload transport continuity during ordinary flows. The functionality of a complete sediment trap during large floods was not the focus of this study. The flume experiments were performed to reproduce morphologically effective ordinary flows, so bedload transport occurs (following the work in [23]). This effective discharge was assumed to be equivalent to a bankfull discharge, which can be linked to return periods of 1 to 3 years for most rivers [24], but could increase to 30 to 50 years for steep mountain streams [8,25]. Following the procedure tested by Schalko and Weitbrecht [22], we broadened the tested parameter range and focused on wood characteristics. Specifically, we studied the impact of different LW sizes, shapes, densities, and organic fine material (i.e., fir branches and leaves [26,27,28]) on backwater rise upstream of the rack and on bedload transport continuity. In this paper, we hypothesized that (H1) LW retention at the inclined bar screen affects the behaviour of the structure, triggering backwater rise and leading to bedload retention; (H2) LW characteristics (i.e., dimension, shape, and density) influence backwater rise and bedload transport continuity; (H3) organic fine material (i.e., fir branches and leaves) within an accumulation results in increased backwater rise and bedload deposition; (H4) the maintenance of the structure requires both LW and bedload management to avoid bedload self-flushing, especially during sequenced floods; and (H5) the design of inclined bar screens, including the lower gap, is not only a function of characteristic grain size diameters, but also of LW dimensions.

2. Materials and Methods

2.1. Experimental Setup

Physical model tests were performed at the Laboratory of Hydraulics, Hydrology and Glaciology (VAW) at ETH Zurich. The experiments were conducted according to Froude similitude with a model scale factor of λ = 30 in an 8 m long, 0.4 m wide, and 0.7 m deep glass-walled channel with a slope of S = 3.3%. The test setup and notation are illustrated in Figure 1. The intake is 0.66 m long and equipped with two flow straighteners to suppress secondary currents and surface waves. Approximately 4 m downstream from the intake an inclined bar screen with aluminum circular bars (5 mm diameter) and a lower gap was placed in the flume to mimic the mechanical control of a semipermeable check-dam within a guiding channel. The inclined bar screen was designed based on the recommendations by Schwindt et al. [9] and Roth et al. [10] with an inclination of 2:1. The concrete orifice dam of the Schwindt et al. [9] design was neglected, as the focus of this study was on the interaction between bedload transport and LW at the inclined bar screen during ordinary flows. The open bar spacing b = 1.2 d 84 and the bottom clearance height f m   d 84 (with f m defined as a pre-factor) are both a function of the characteristic grain size diameter d 84 and were kept constant during the experiments with b = 7 mm and f m   d 84 = 19 mm (similar relative bottom clearance height f m   d 84 / h o = 0.56 to Roth et al. [10]). The model bar diameter d b was set to 5 mm, in the same order of magnitude as d 84 . Steady-state discharge with Q = 11.6 l/s was used to reproduce morphologically effective ordinary flows. The tested Q and S correspond to supercritical flow conditions with an approach flow (subscript o) Froude number F o = v/ g h ≥ 1, with v flow velocity, g gravitational acceleration, and h water depth. The discharge was measured with an electromagnetic flow meter (MID) with ±0.1% accuracy and the water depth was measured with two Ultrasonic Distance Sensors (UDS) positioned 1.2 m and 0.3 m upstream and one UDS 1.2 m downstream of the inclined bar screen with ±0.3 mm accuracy. For all tests, backwater rise ( Δ h ) was calculated as the difference between upstream and downstream water depth to Δ h = h h o with relative backwater rise defined as Δ h / h o .

2.2. Model Large Wood and Sediment

Natural cylindrical wooden logs were used to model LW. In addition, logs with branches, rootwads, and model fine material (using plastic fir tree made of polyvinyl chloride) were also tested (Figure 2). The log diameter varied between d L = 0.2 and 2 cm and log length L L = 4.7 to 20 cm. LW density ρ L was varied by watering the logs for several days until the respective density of 500 to 1100 kg/m 3 was reached. The fine material and rootwads were added to the solid LW volume V s as a volume percentage with 10, 30, and 50% of fine material and 25 and 50% of rootwads. A total of 200 g of LW was added in four wood packages to the flow 0.5 m downstream of the inlet. Each wood package consisted of a solid wood mass M s = 50 g (i.e., 6–7 logs for V1a–c). Based on the solid wood mass M s and the wood density ρ L , the respective solid wood volume V s can be calculated for each wood package as well as the relative wood volume V s , r e l , defined as the wood volume required to block an idealized rectangular cuboid V s , r e l = V s / ( B h o 2 ) , with B channel width. Considering the different tested LW characteristics, V s , r e l = 0.14 to 0.30 for each wood package. The wood trapping ratio T R = V s / V s , a l l is the ratio of the wood volume blocked at the inclined bar screen to the entire wood volume added to the flow.
Uniform bed material was used as model sediment with the characteristic grain size diameters of d m = 5.5 mm (subscript m as mean), d 90 = 6.0 mm (0.18 m in prototype with λ = 30 to model coarse gravel bed rivers), d 84 = 5.9 mm, and d 16 = 5.2 mm (see in [22]). The sorting coefficient of the grain size distribution was σ g = 1.1 and verified with a sieving analysis. Prior to the addition of LW, the initial condition of bedload transport was defined as apparent transport across channel width with the initial (subscript IC) non-dimensional bed shear stress θ I C ≥ 0.11 and larger than the threshold value for full transport in gravel bed rivers θ c r = 0.06 [29]. The corresponding bedload transport input rate G b , i n in (g/s) was calculated using the bedload transport equation by Rickenmann [30], resulting in G b , i n = 88 g/s for the tested flow and bedload conditions. The calculated value was then verified during the experiments and added continuously to the flow over the entire test duration using a sediment feeder. The experiments were conducted for 100% transport capacity, so bedload transport input equals output rate. The bedload transport output rate G b , o u t was manually measured for a defined duration (≈5 min) using a filter basket, and then drying in the oven and weighing the sediment with ±0.5 g/s accuracy. The bedload transport reduction factor ξ was computed as the ratio between the bedload output rate and bedload input rate to ξ = G b , o u t / G b , i n . A value of ξ c r = 0.75 (i.e., the bedload output is reduced to 75% of the bedload supply and 25% of the input bedload is retained due to the inclined bar screen) was defined as the critical threshold for maintaining bedload transport continuity [10,22]. During the experiments, ξ was determined after each LW package addition.

2.3. Test Program and Procedure

The test program is listed in Table 1. The focus of the experiments in this paper was put on the wood characteristics. Therefore, the rack geometry was kept constant with a bottom clearance height of f m   d 84 = 1.9 cm as well as the initial conditions with discharge Q = 11.6 l/s, approach flow depth h o = 3.4 cm, flow Froude number F o = 1.55, channel slope S = 3.3%, non-dimensional bed shear stress θ = 0.11, and bedload transport rate G b , i n = 88 g/s. The reproducibility of the experiments was evaluated in tests V1a–c. In tests V1–V5, the effect of log density was studied. Different log dimensions were investigated in tests V6–V8. Model fine material was added to the LW accumulation in tests V9–V11. In tests V12–V14, the effect of branches and rootwads was studied. For test V15, the test procedure was adapted to model a sequenced flood and consisted of two runs. Similar to the previous tests, four LW packages (i.e., wood supply 1) were added to the flow during the first run. After the addition of the last LW package, the accumulated LW was removed and the second run started with the addition of a new set of four LW packages (i.e., wood supply 2) to the flow without changing flow and bedload conditions. The time between wood removal and the beginning of the second run (wood supply 2) amounted to 3 min, similar to the time interval between the addition of LW packages.
The test procedure can be summarized by the following steps. First, the approach flow conditions were set by gradually increasing Q with S = 0% at test start. The channel slope was then increased to S = 3.3% and the bedload addition started. The initial bedload transport condition was evaluated by measuring h o with UDS and calculating θ I C . Equilibrium of bedload transport ( G b , i n = G b , o u t ; 100% transport capacity) was validated by manually measuring G b , o u t . Second, the inclined bar screen was inserted to the flume with the defined bottom clearance height f m   d 84 , and both water depth h and bedload transport output rate G b , o u t were measured to check the influence of the inclined bar screen on flow and bedload transport continuity. Prior to rack addition, ξ = 99.2% (±1.5) compared to ξ = 99.0% (±1.3) with the rack, demonstrating that the rack itself had a negligible effect on ξ . Third, LW was added to the flow in random orientation and distributed across the flume width. After each LW package, h and G b , o u t were measured after an equilibrium was reached with respect to h (≈2 min). After the addition of the respective LW volume (i.e., four LW packages), S was reduced to 0%, bedload supply stopped, and Q reduced. LW that accumulated at the inclined bar screen was extracted and dried.

3. Results and Discussion

3.1. Test Reproducibility and Wood Trapping Ratio

Test reproducibility was confirmed by tests V1a–c (Figure 3). Relative backwater rise Δ h / h o increased and relative bedload transport reduction factor ξ decreased with increasing relative solid wood volume V s , r e l . The mean relative standard error amounted to 6.2% for Δ h / h o and 7.2% for ξ (Table 2). The differences between the three tests were larger for V s , r e l < 0.5. The LW accumulation process is highly stochastic, leading to different log arrangements for the respective tests. With increasing V s , r e l and number of logs, the differences in log arrangement (i.e., blocking ratio) decreased and likewise the differences in the resulting backwater rise and bedload transport reduction factor. For the later analysis, the average of tests V1a–c (V1*) was used to compare the effects of wood characteristics on backwater rise and bedload transport continuity. For all tests, the trapping ratio T R of LW at the inclined bar screen was on average 89.6%, with a minimum T R of 59.2% for test V7 due to the small log sizes.

3.2. Large Wood Density

A photo series of resulting backwater rise and bedload deposition at the inclined bar screen is illustrated for test V1b with ρ L = 500 kg/m 3 in Figure 4a,b and test V4 with ρ L = 900 kg/m 3 in Figure 4c,d. We observed that with increasing ρ L and in particular for ρ L > 1000 kg/m 3 , logs tend to be transported closer to the channel bottom. In contrast, logs with ρ L < 1000 kg/m 3 were transported at the water surface (i.e., floating). This affected the blocking position of individual logs at the inclined bar screen as well as the possibility of dense logs being transported through the bottom clearance further downstream. The trapping ratio T R decreased with increasing ρ L and was 96% for ρ L = 500 kg/m 3 (V1*) compared to 71% for ρ L = 1100 kg/m 3 (V5).
According to Figure 5a, relative backwater rise Δ h / h o increased with increasing relative solid LW volume V s , r e l . With increasing V s , r e l , the open flow cross section decreases, leading to a decrease in flow velocity and, thus, increase in resulting flow depth h. However, there was a negligible effect of ρ L on Δ h / h o , as the mean relative standard error amounted to 7.2% and is comparable to the mean relative standard error of the reproducibility tests (6.2%) and no systematic trend was detected. The LW accumulation at the rack and the resulting increase in Δ h / h o will also have implications on the bedload transport capacity. With increasing Δ h / h o , the respective flow velocity and non-dimensional bed shear stress are reduced. Given θ < θ c r , bedload will deposit upstream of the inclined bar screen (Figure 4). As soon as the first LW package was retained at the inclined bar screen, the formation of a bedload deposition was observed. Initially, the bedload deposition was formed upstream of the inclined bar screen and evolved in the upstream direction in a succession of quasi-equilibrium states [31,32]. The change in bottom slope may also affect the flow conditions, and further enhance bedload deposition and reduce bedload transport downstream. Depending on the extent, bedload transport continuity will be significantly reduced and the bedload transport reduction factor will fall below the threshold value ξ < ξ c r (herein defined as ξ c r = 0.75). Similar to the results for Δ h / h o , no systematic trend of the effect of ρ L on ξ can be deduced based on Figure 5b. Given logs with ρ L = 700 and 800 kg/m 3 (V2 and V3), ξ was reduced to 0.55 and 0.34, respectively, after the addition of V s , r e l = 0.18 and V s , r e l = 0.21, respectively. The variation in ξ in Figure 5b resulted due to the differences in log arrangement at the inclined bar screen and the stochastic LW accumulation process. For all tests, ξ < ξ c r after the addition of V s , r e l > 0.50, corresponding to ≈150 g of LW (model scale). The results for the tests with varying ρ L confirmed our first hypothesis (H1) that LW retention at the inclined bar screen affects the behaviour of the structure by triggering backwater rise and leading to bedload retention. The second hypothesis (H2) was not confirmed as the effect of ρ L on both backwater rise and bedload retention was negligible.

3.3. Large Wood Size

Figure 6 shows a photo series of resulting backwater rise and bedload deposition at the inclined bar screen for test V7 with d L = 0.5 cm (Figure 6a,b) compared to test V8 with d L = 0.2 cm (Figure 6c,d). Similar to previous studies on backwater rise due to LW accumulations at rack structures, Δ h / h o increased with decreasing LW size. Due to smaller d L , the flow within an accumulation is more heterogeneous, leading to increased friction losses and increased Δ h / h o [26,28]. Given V s , r e l = 0.5, Δ h / h o increased by a factor of ≈2.7 for smaller wood dimensions ( d L = 0.2 cm and L L = 5 cm; V8) compared to the reference test ( d L = 1 cm and L L = 15 cm; V1*) (Figure 7a). The increase in Δ h / h o led to a reduction in bedload transport, and ξ decreased with increasing V s , r e l and decreasing LW size (Figure 7b). After the addition of the first LW package ( V s , r e l ≈ 0.15), bedload transport was significantly reduced to ξ = 0.25 for the smallest tested wood size (V8) compared to ξ = 0.85 for the reference test (V1*) (Figure 7b).
For all tests, ξ decreased below ξ c r = 0.75 for V s , r e l > 0.38. Due to the decrease in ξ , bedload deposited upstream of the rack and the bedload deposition increased with increasing V s , r e l and decreasing LW size (i.e., with increasing Δ h / h o ) (Figure 6). As observed for the tests with varying LW density, the bedload deposition upstream of the rack was extending further upstream with increasing V s , r e l . According to Figure 7b, the effect of LW accumulation on bedload transport was not linear and different for the tested LW sizes. During the experiments, we noted that the log position can shift, thereby affecting the open flow cross section and also the resulting bedload transport. The results of tests V6 and V8 further show that ξ increased for V s , r e l > 0.50 due to partial self-flushing of the bedload deposition upstream of the inclined bar screen. The self-flushing may be initiated due to the changes in log position at the bar screen or due to the increasing slope of the bedload deposition upstream of the bar screen. The experimental results for tests with different LW size confirmed both the first and second hypotheses (H1 and H2), as the LW size had a governing effect on backwater rise and bedload transport continuity.

3.4. Model Fine Material

The presence of fine material ( F M ) such as fir branches and leaves in the LW accumulation significantly affected both backwater rise and bedload transport continuity (Figure 8 and Figure 9). Previous studies [26,27,28] demonstrated that F M in an accumulation lead to a more heterogeneous flow within the accumulation, increased friction losses, and, thus, increased backwater rise (Figure 9a). To show the effect of F M , the volume V s , r e l was determined with the logs only.
The addition of 30% and 50% F M increased Δ h / h o by a factor of 1.9 and 2.3, respectively, compared to the reference test without F M (V1*) for V s , r e l ≈ 0.2. By conducting flume experiments with bedload supply, we noted that F M significantly decreased bedload transport (Figure 9b). If we compare ξ for V s , r e l ≈ 0.5, bedload transport capacity is reduced to ξ = 0.18 for 50% F M (V11), ξ = 0.30 for 30% F M (V10), and ξ = 0.47 for 10% F M (V8) compared to ξ = 0.65 without F M (V1*) (Figure 9b). The resulting decrease in bedload transport capacity leads to bedload deposition upstream of the rack, as illustrated in Figure 8. Based on the results according to Figure 9b, F M significantly affects bedload transport capacity and, thus, the functionality of the inclined bar screen to enable bedload transport in the presence of LW. The results, therefore, confirmed the first and third hypotheses (H1 and H3), as the consideration of F M in a LW accumulation had a governing effect on backwater rise and bedload retention.

3.5. Large Wood Shape

For selected tests, the accumulation and resulting bedload deposition was investigated for different LW shapes including logs with branches and rootwads (Figure 10). In the prototype, especially freshly recruited wood may still include rootwads and branches. In contrast to LW size such as log length and diameter, the shape of LW did not have a governing effect on backwater rise and bedload transport capacity (Figure 11). In addition, the differences in both Δ h / h o and ξ are minimal for test V13 with 25% R W compared to test V14 with 50% R W . According to Figure 11b, ξ decreased below ξ c r = 0.75 for V s , r e l > 0.50. The results of tests V13 and V14 further show that ξ increased for V s , r e l > 0.75 due to partial self-flushing of the bedload deposition upstream of the inclined bar screen, as also observed for tests V6 and V8. To improve the process understanding of self-flushing, continuous monitoring of the bedload transport rate is required. The results for the tests with varying LW shape also confirmed the first hypothesis (H1), while the second hypothesis (H2) was not confirmed as the effect of LW shape on both backwater rise and bedload retention was negligible.

3.6. Sequenced Floods

For the last test (V15), the test procedure was adapted by modeling a sequence of floods with two consecutive wood supplies. For the first run, LW was added to the flow in four packages similar to the previous tests (Figure 12b–e). After the addition of the last LW package (Figure 12e), the retained LW at the bar screen was removed and the second run started (Figure 12f) by repeating the procedure and adding four LW packages to the flow (Figure 12g–j). Similar to the previous tests, Δ h / h o increased with increasing V s , r e l . Due to the bedload deposition during the first run, Δ h / h o = 0.36 after the removal of the first wood supply and prior to the second wood supply (Figure 13a). Given V s , r e l = 0.5, Δ h / h o increased by a factor of 1.7 during the second run compared to the first run (Figure 13a). The LW removal at the inclined bar screen led to bedload self-flushing indicated by an increase in bedload transport with ξ ≈ 1.5 prior to the wood supply in the second run (Figure 13b). For a LW accumulation of V s , r e l ≈ 0.6, ξ was larger by a factor of 1.5 in the first run compared to the second run (Figure 13b).
The results demonstrated that the consecutive LW supply can affect both backwater rise and bedload deposition. They confirmed the fourth hypothesis (H4), highlighting that the maintenance of inclined bar screens requires a combined management of both LW and bedload in order to avoid bedload deposition during ordinary flows, bedload self-flushing, as well as backwater rise.

3.7. Interaction of Large Wood Accumulation and Bedload Transport Continuity

The inclined bar screen was designed based on recommendations by Schwindt et al. [9] and Roth et al. [10], where the bottom clearance height is a function of the characteristic grain size diameter d 84 . In order to also apply this novel rack structure to sites with LW transport during ordinary flows, it is important to study how LW deposits at the structure alter bedload transport capacity. To determine the effect of LW accumulation on bedload transport continuity, we defined that bedload transport is impaired, if the bedload transport output rate G b , o u t is reduced by a minimum of 25% compared to the bedload transport input rate G b , i n , i.e., ξ c r = 0.75 [10]. Similar to the results of Roth et al. [10], bedload transport continuity was not affected by the rack structure prior to the addition of LW with a bedload transport reduction factor ξ ≈ 1 (Figure 14). By plotting the bedload transport reduction factor ξ versus the relative solid LW volume V s , r e l for all tests with varied LW characteristics (V1*–V14), we can differentiate between two areas (gray- vs. white-shaded area in Figure 14). Due to the stochastic LW accumulation process, there are significant deviations in ξ for the different LW characteristics. However, the white-shaded area illustrates that for V s , r e l > 0.5, all tests fell below the threshold value of ξ c r and bedload transport capacity is impaired. The results, therefore, demonstrate that a comparably small LW volume is sufficient to reduce bedload transport. Similar to the results of Schalko and Weitbrecht [22], a LW accumulation of ≈10–20 logs (100–150 g model scale) with V s , r e l > 0.5 at an inclined bar screen may be large enough to already impair bedload transport. The results plotted in Figure 14 support the last hypothesis (H5) that the design of inclined bar screens and specifically the design of the bottom clearance, should also include LW dimensions to enable bedload transport continuity during ordinary flows. However, additional experiments are required to define the bottom clearance considering LW size and to confirm this hypothesis (H5).

4. Summary and Outlook

Flume experiments were conducted to study the effect of LW accumulations at an inclined bar screen on both backwater rise and bedload transport continuity with a special focus on varying LW characteristics. The results demonstrated that the inclined bar screen does not affect bedload transport prior to the addition of LW with a bedload reduction coefficient of ξ ≈ 1. However, a relative solid LW volume of V s , r e l > 0.5 or ≈10–20 logs (corresponding to a wood volume of ≈6.8 m 3 prototype scale with a model scale factor of 30) are sufficient to reduce the bedload transport capacity to below 75% compared to the condition without LW. The results further demonstrated that smaller LW sizes (i.e., log diameter and length), the addition of fine material F M such as fir branches and leaves, and the simulation of a sequenced flood significantly increased backwater rise Δ h / h o and reduced the bedload transport reduction factor ξ compared to the reference test with regular logs. In contrast, LW density and LW shape such as logs with branches and the addition of rootwads resulted in similar backwater rise and bedload transport capacity as the reference test. However, in prototype, decayed wood (i.e., small density) may be transported during smaller flows compared to dense wood.
The results further confirm our hypotheses that a LW accumulation at an inclined bar screen affects the behavior of the structure (H1) and that LW characteristics such as fine material and LW size significantly influence backwater rise and bedload retention (H2, H3). The test focusing on the sequenced flood highlighted the need for improved maintenance measures in case an inclined bar screen is used for a combined LW and bedload retention to avoid self-flushing of the bed material (H4). Furthermore, the design of an inclined bar screen may need to be adapted for a combined LW and bedload retention to enable bedload transport continuity during ordinary flows (H5). As the presented experiments did not vary the design of an inclined bar screen, i.e., bottom clearance, the last hypothesis (H5) cannot yet be confirmed.
Additional flume experiments are required focusing on (1) design alterations of the inclined bar screen as a function of LW size, (2) unsteady flow conditions, and (3) different bed material. Variation of the bottom clearance needs to be tested in combination with large floods to guarantee the functionality of inclined bar screens as well as the entire sediment trap system during typical flood hydrographs.

Author Contributions

Conceptualization, I.S.; methodology, I.S. and V.R.-V.; investigation, V.R.-V. and F.M.; formal analysis, I.S. and V.R.-V.; writing—original draft preparation, I.S. and V.R.-V.; writing—review and editing, I.S., V.R.-V., F.M., and V.W.; visualization, I.S.; supervision, V.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study is available in Schalko et al. [34].

Acknowledgments

Kristina Milojevic is acknowledged for her laboratory work during her project thesis at ETH Zurich.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Piton, G.; Recking, A. Design of sediment traps with open check dams. II: Woody debris. J. Hydraul. Eng. 2016, 142, 04015046. [Google Scholar] [CrossRef]
  2. Armanini, A.; Dellagiacoma, F.; Ferrari, L. From the check dam to the development of functional check dams. In Fluvial Hydraulics of Mountain Regions; Springer: Berlin/Heidelberg, Germany, 1991; pp. 331–344. [Google Scholar]
  3. Wehrmann, H.; Hübl, J.; Holzinger, G. Classification of dams in torrential watersheds. In Proceedings of the INTERPRAEVENT Conference Citeseer, Niigata, Japan, 25–29 September 2006; pp. 829–838. [Google Scholar]
  4. Rudolf-Miklau, F.; Suda, J. Design criteria for torrential barriers. In Dating Torrential Processes on Fans and Cones; Springer: Berlin/Heidelberg, Germany, 2013; pp. 375–389. [Google Scholar]
  5. Boix-Fayos, C.; Barberá, G.G.; López-Bermúdez, F.; Castillo, V.M. Effects of check dams, reforestation and land-use changes on river channel morphology: Case study of the Rogativa catchment (Murcia, Spain). Geomorphology 2007, 91, 103–123. [Google Scholar] [CrossRef]
  6. Bombino, G.; Boix-Fayos, C.; Gurnell, A.M.; Tamburino, V.; Zema, D.A.; Zimbone, S.M. Check dam influence on vegetation species diversity in mountain torrents of the Mediterranean environment. Ecohydrology 2014, 7, 678–691. [Google Scholar] [CrossRef]
  7. Galia, T.; Škarpich, V.; Ruman, S.; Macurová, T. Check dams decrease the channel complexity of intermediate reaches in the Western Carpathians (Czech Republic). Sci. Total. Environ. 2019, 662, 881–894. [Google Scholar] [CrossRef] [PubMed]
  8. Schwindt, S. Hydro-Morphological Processes through Permeable Sediment Traps at Mountain Rivers. Ph.D. Thesis, EPFL, Lausanne, Switzerland, 2017. [Google Scholar]
  9. Schwindt, S.; Franca, M.J.; Reffo, A.; Schleiss, A.J. Sediment traps with guiding channel and hybrid check dams improve controlled sediment retention. Nat. Hazards Earth Syst. Sci. 2018, 18, 647–668. [Google Scholar] [CrossRef] [Green Version]
  10. Roth, A.; Jafarnejad, M.; Schwindt, S.; Schleiss, A. Design optimization of permeable sediment traps for fluvial bed load transport. E3S Web of Conferences. Edp Sci. 2018, 40, 03009. [Google Scholar]
  11. Moldenhauer-Roth, A.; Piton, G.; Schwindt, S.; Jafarnejad, M.; Schleiss, A.J. Design of sediment detention basins: Scaled model experiments and application. Int. J. Sediment Res. 2021, 36, 136–150. [Google Scholar] [CrossRef]
  12. Keller, E.A.; Swanson, F.J. Effects of large organic material on channel form and fluvial processes. Earth Surf. Process. 1979, 4, 361–380. [Google Scholar] [CrossRef]
  13. Nakamura, F.; Swanson, F.J. Effects of coarse woody debris on morphology and sediment storage of a mountain stream system in western Oregon. Earth Surf. Process. Landforms 1993, 18, 43–61. [Google Scholar] [CrossRef]
  14. Piégay, H.; Gurnell, A. Large woody debris and river geomorphological pattern: Examples from SE France and S. England. Geomorphology 1997, 19, 99–116. [Google Scholar] [CrossRef]
  15. Rimböck, A. Schwemmholzrückhalt in Wildbächen: Grundlagen zu Planung und Berechnung von Seilnetzsperren. Ph.D. Thesis, Technische Universität München, Munich, Germany, 2003. [Google Scholar]
  16. Piton, G.; Recking, A. Design of sediment traps with open check dams. I: Hydraulic and deposition processes. J. Hydraul. Eng. 2016, 142, 04015045. [Google Scholar] [CrossRef]
  17. Schmocker, L.; Weitbrecht, V. Driftwood: Risk Analysis and Engineering Measures. J. Hydraul. Eng. 2013, 139, 683–695. [Google Scholar] [CrossRef]
  18. Piton, G.; Mano, V.; Richard, D.; Evin, G.; Laigle, D.; Tacnet, J.M.; Rielland, P.A. Design of a debris retention basin enabling sediment continuity for small events: The Combe de Lancey case study (France). In Proceedings of the International Conference on Debris-Flow Hazards Mitigation: Mechanics, Prediction, and Assessment, Golden, CO, USA, 10–13 June 2019; pp. 1019–1026. [Google Scholar]
  19. Piton, G.; Horiguchi, T.; Marchal, L.; Lambert, S. Open check dams and large wood: Head losses and release conditions. Nat. Hazards Earth Syst. Sci. 2020, 20, 3293–3314. [Google Scholar] [CrossRef]
  20. Lange, D.; Bezzola, G.R. Schwemmholz: Probleme und Lösungsansätze. VAW-Mitteilungen 2006, 188, 1–135. [Google Scholar]
  21. Schalko, I. Wood retention at inclined racks: Effects on flow and local bedload processes. Earth Surf. Process. Landforms 2020, 45, 2036–2047. [Google Scholar] [CrossRef]
  22. Schalko, I.; Weitbrecht, V. Wood blockage and sediment transport at inclined bar screens. J. Hydraul. Res. 2021, 1–9. [Google Scholar] [CrossRef]
  23. Wolman, M.G.; Miller, J.P. Magnitude and frequency of forces in geomorphic processes. J. Geol. 1960, 68, 54–74. [Google Scholar] [CrossRef] [Green Version]
  24. Wohl, E.E. Mountain Rivers; American Geophysical Union: Washington, DC, USA, 2000. [Google Scholar]
  25. Hassan, M.A.; Brayshaw, D.; Alila, Y.; Andrews, E. Effective discharge in small formerly glaciated mountain streams of British Columbia: Limitations and implications. Water Resour. Res. 2014, 50, 4440–4458. [Google Scholar] [CrossRef]
  26. Schalko, I.; Schmocker, L.; Weitbrecht, V.; Boes, R.M. Backwater Rise due to Large Wood Accumulations. J. Hydraul. Eng. 2018, 144, 04018056. [Google Scholar] [CrossRef]
  27. Schalko, I.; Lageder, C.; Schmocker, L.; Weitbrecht, V.; Boes, R.M. Laboratory flume experiments on the formation of spanwise large wood accumulations Part I: Effect on backwater rise. Water Resour. Res. 2019, 55, 4854–4870. [Google Scholar] [CrossRef]
  28. Follett, E.; Schalko, I.; Nepf, H. Momentum and Energy Predict the Backwater Rise Generated by a Large Wood Jam. Geophys. Res. Lett. 2020, 47. [Google Scholar] [CrossRef]
  29. 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–2029. [Google Scholar] [CrossRef] [Green Version]
  30. Rickenmann, D. Bedload Transport Capacity of Slurry Flows at Steep Slopes. Ph.D. Thesis, ETH Zurich, Zurich, Switzerland, 1990. [Google Scholar]
  31. Armanini, A.; Larcher, M. Rational Criterion for Designing Opening of Slit-Check Dam. J. Hydraul. Eng. 2001, 127, 94–104. [Google Scholar] [CrossRef]
  32. Campisano, A.; Cutore, P.; Modica, C. Improving the Evaluation of Slit-Check Dam Trapping Efficiency by Using a 1D Unsteady Flow Numerical Model. J. Hydraul. Eng. 2014, 140, 04014024. [Google Scholar] [CrossRef]
  33. Schalko, I.; Ruiz-Villanueva, V.; Weitbrecht, V. Effect of wood accumulation on sediment continuity at permeable sediment traps. In Proceedings of the 10th Conference on Fluvial Hydraulics, Delft, The Netherlands, 7–10 July 2020; CRC Press: Boca Raton, FL, USA, 2020; pp. 1–6. [Google Scholar]
  34. Schalko, I.; Ruiz-Villanueva, V.; Maager, F.; Weitbrecht, V. Dataset: Wood Retention at Inclined Bar Screens: Effect of Wood Characteristics on Backwater Rise and Bedload Transport. Zenodo. Available online: https://doi.org/10.5281/zenodo.5040461 (accessed on 30 June 2021). [CrossRef]
Figure 1. Test setup and notation with Froude number F o [-], discharge Q [l/s], approach flow depth h o [cm], resulting water depth with wood blockage h [cm], backwater rise Δ h [cm], and bottom clearance height f m   d 84 [mm]. Ultrasonic distance sensors (UDS) to measure h. LW accumulation with log length L L [cm], log diameter d L [cm], and solid wood volume V s [ m 3 ]. Bedload transport input rate G b , i n [g/s] and output rate G b , o u t [g/s] with mean sediment diameter d m [mm] and sorting coefficient of the grain size distribution σ g [-]. Potential bedload deposition is illustrated upstream of the inclined bar screen.
Figure 1. Test setup and notation with Froude number F o [-], discharge Q [l/s], approach flow depth h o [cm], resulting water depth with wood blockage h [cm], backwater rise Δ h [cm], and bottom clearance height f m   d 84 [mm]. Ultrasonic distance sensors (UDS) to measure h. LW accumulation with log length L L [cm], log diameter d L [cm], and solid wood volume V s [ m 3 ]. Bedload transport input rate G b , i n [g/s] and output rate G b , o u t [g/s] with mean sediment diameter d m [mm] and sorting coefficient of the grain size distribution σ g [-]. Potential bedload deposition is illustrated upstream of the inclined bar screen.
Water 13 02231 g001
Figure 2. Model wood, fine material, and sediment.
Figure 2. Model wood, fine material, and sediment.
Water 13 02231 g002
Figure 3. (a) Relative backwater rise Δ h / h o versus relative solid wood volume V s , r e l for tests V1a–c, (b) bedload transport reduction factor ξ versus V s , r e l for V1a–c, and (c) wood trapping ratio T R = V s / V s , a l l for all tests.
Figure 3. (a) Relative backwater rise Δ h / h o versus relative solid wood volume V s , r e l for tests V1a–c, (b) bedload transport reduction factor ξ versus V s , r e l for V1a–c, and (c) wood trapping ratio T R = V s / V s , a l l for all tests.
Water 13 02231 g003
Figure 4. LW accumulation and bedload deposition at the inclined bar screen for (a) test V1b ( ρ L = 500 kg/m 3 ) with V s , r e l = 0.58 and (b) V s , r e l = 0.74, compared to (c) test V4 ( ρ L = 900 kg/m 3 ) with V s , r e l = 0.48 and (d) V s , r e l = 0.64.
Figure 4. LW accumulation and bedload deposition at the inclined bar screen for (a) test V1b ( ρ L = 500 kg/m 3 ) with V s , r e l = 0.58 and (b) V s , r e l = 0.74, compared to (c) test V4 ( ρ L = 900 kg/m 3 ) with V s , r e l = 0.48 and (d) V s , r e l = 0.64.
Water 13 02231 g004
Figure 5. (a) Relative backwater rise Δ h / h o and (b) bedload transport reduction factor ξ versus V s , r e l for tests V1* ( ρ L = 500 kg/m 3 ), V2 ( ρ L = 700 kg/m 3 ), V3 ( ρ L = 800 kg/m 3 ), V4 ( ρ L = 900 kg/m 3 ), and V5 ( ρ L = 1100 kg/m 3 ).
Figure 5. (a) Relative backwater rise Δ h / h o and (b) bedload transport reduction factor ξ versus V s , r e l for tests V1* ( ρ L = 500 kg/m 3 ), V2 ( ρ L = 700 kg/m 3 ), V3 ( ρ L = 800 kg/m 3 ), V4 ( ρ L = 900 kg/m 3 ), and V5 ( ρ L = 1100 kg/m 3 ).
Water 13 02231 g005
Figure 6. LW accumulation and bedload deposition at the inclined bar screen for (a) test V7 ( d L = 0.5 cm and L L = 5–10 cm) with V s , r e l = 0.17 and (b) V s , r e l = 0.34, compared to (c) test V8 ( d L = 0.2 cm and L L = 5 cm) with V s , r e l = 0.15 and (d) V s , r e l = 0.30.
Figure 6. LW accumulation and bedload deposition at the inclined bar screen for (a) test V7 ( d L = 0.5 cm and L L = 5–10 cm) with V s , r e l = 0.17 and (b) V s , r e l = 0.34, compared to (c) test V8 ( d L = 0.2 cm and L L = 5 cm) with V s , r e l = 0.15 and (d) V s , r e l = 0.30.
Water 13 02231 g006
Figure 7. (a) Relative backwater rise Δ h / h o and (b) bedload transport reduction factor ξ versus V s , r e l for tests V1* ( d L = 1 cm and L L = 15–20 cm), V6 ( d L = 2 cm and L L = 15–20 cm), V7 ( d L = 0.5 cm and L L = 5–10 cm), and V8 ( d L = 0.2 cm and L L = 5 cm).
Figure 7. (a) Relative backwater rise Δ h / h o and (b) bedload transport reduction factor ξ versus V s , r e l for tests V1* ( d L = 1 cm and L L = 15–20 cm), V6 ( d L = 2 cm and L L = 15–20 cm), V7 ( d L = 0.5 cm and L L = 5–10 cm), and V8 ( d L = 0.2 cm and L L = 5 cm).
Water 13 02231 g007
Figure 8. LW accumulation and bedload deposition at the inclined bar screen for test V11 with 50% F M and increasing V s , r e l : (a) V s , r e l = 0.42 and (b) V s , r e l = 0.83 (adapted from the work in [33]).
Figure 8. LW accumulation and bedload deposition at the inclined bar screen for test V11 with 50% F M and increasing V s , r e l : (a) V s , r e l = 0.42 and (b) V s , r e l = 0.83 (adapted from the work in [33]).
Water 13 02231 g008
Figure 9. (a) Relative backwater rise Δ h / h o and (b) bedload transport reduction factor ξ versus V s , r e l for tests V1* (no F M ), V9 (10% F M ), V10 (30% F M ), and V11 (50% F M ).
Figure 9. (a) Relative backwater rise Δ h / h o and (b) bedload transport reduction factor ξ versus V s , r e l for tests V1* (no F M ), V9 (10% F M ), V10 (30% F M ), and V11 (50% F M ).
Water 13 02231 g009
Figure 10. LW accumulation and bedload deposition at the inclined bar screen for test V14 with 50% R W and increasing V s , r e l : (a) side view and (b) zoom with V s , r e l = 0.52 and (c) side view and (d) zoom with V s , r e l = 1.04.
Figure 10. LW accumulation and bedload deposition at the inclined bar screen for test V14 with 50% R W and increasing V s , r e l : (a) side view and (b) zoom with V s , r e l = 0.52 and (c) side view and (d) zoom with V s , r e l = 1.04.
Water 13 02231 g010
Figure 11. (a) Relative backwater rise Δ h / h o and (b) bedload transport reduction factor ξ versus V s , r e l for tests V1* (regular logs), V12 (logs with branches), V13 (25% R W ), and V14 (50% R W ).
Figure 11. (a) Relative backwater rise Δ h / h o and (b) bedload transport reduction factor ξ versus V s , r e l for tests V1* (regular logs), V12 (logs with branches), V13 (25% R W ), and V14 (50% R W ).
Water 13 02231 g011
Figure 12. LW accumulation and bedload deposition at the inclined bar screen for test V15 with increasing V s , r e l . Panels (ae) illustrate first run and panels (fj) second run. (a) Prior to first wood supply, (b) V s , r e l = 0.16, (c) V s , r e l = 0.31, (d) V s , r e l = 0.47, (e) V s , r e l = 0.62, (f) after wood removal of first run and prior to second wood supply, (g) V s , r e l = 0.15, (h) V s , r e l = 0.30, (i) V s , r e l = 0.45, (j) V s , r e l = 0.60.
Figure 12. LW accumulation and bedload deposition at the inclined bar screen for test V15 with increasing V s , r e l . Panels (ae) illustrate first run and panels (fj) second run. (a) Prior to first wood supply, (b) V s , r e l = 0.16, (c) V s , r e l = 0.31, (d) V s , r e l = 0.47, (e) V s , r e l = 0.62, (f) after wood removal of first run and prior to second wood supply, (g) V s , r e l = 0.15, (h) V s , r e l = 0.30, (i) V s , r e l = 0.45, (j) V s , r e l = 0.60.
Water 13 02231 g012
Figure 13. (a) Relative backwater rise Δ h / h o versus V s , r e l for test V1* (addition of 4 LW packages; white circles), V15 with first run (wood supply 1: addition of 4 LW packages; gray squares), and V15 with second run (wood supply 2: LW removal from first run and then re-addition of 4 LW packages; black squares); (b) bedload transport reduction factor ξ versus V s , r e l for test V15 (first and second run).
Figure 13. (a) Relative backwater rise Δ h / h o versus V s , r e l for test V1* (addition of 4 LW packages; white circles), V15 with first run (wood supply 1: addition of 4 LW packages; gray squares), and V15 with second run (wood supply 2: LW removal from first run and then re-addition of 4 LW packages; black squares); (b) bedload transport reduction factor ξ versus V s , r e l for test V15 (first and second run).
Water 13 02231 g013
Figure 14. Bedload transport reduction factor ξ versus V s , r e l for all tests. Circles correspond to different LW density, diamonds to LW size, triangles to fine material, and squares to LW shapes including rootwads.
Figure 14. Bedload transport reduction factor ξ versus V s , r e l for all tests. Circles correspond to different LW density, diamonds to LW size, triangles to fine material, and squares to LW shapes including rootwads.
Water 13 02231 g014
Table 1. Test program with bottom clearance height f m   d 84 = 19 mm, discharge Q = 11.6 l/s, approach flow depth h o = 3.4 cm, flow Froude number F o = 1.55, channel slope S = 3.3%, non-dimensional bed shear stress θ = 0.11, and bedload transport input rate G b , i n = 88 g/s. F M = Fine material and R W = rootwads.
Table 1. Test program with bottom clearance height f m   d 84 = 19 mm, discharge Q = 11.6 l/s, approach flow depth h o = 3.4 cm, flow Froude number F o = 1.55, channel slope S = 3.3%, non-dimensional bed shear stress θ = 0.11, and bedload transport input rate G b , i n = 88 g/s. F M = Fine material and R W = rootwads.
TestLog Diameter d L
[cm]
Log Length L L [cm]Wood Shape [-]Wood Density
ρ L [kg/m 3 ]
V1a115–20Logs500
V1b115–20Logs500
V1c115–20Logs500
V2115–20Logs700
V3115–20Logs800
V4115–20Logs900
V5115–20Logs1100
V6215–20Logs620
V70.55–10Logs620
V80.25Logs500
V9115–20Logs + 10% F M 500
V10115–20Logs + 30% F M 500
V11115–20Logs + 50% F M 500
V12115–20Logs with branches660
V13115–20Logs + 25% R W 500
V14115–20Logs + 50% R W 500
V15115–20Logs500
Table 2. Relative standard error and mean relative standard error of reproducibility tests (V1a–c).
Table 2. Relative standard error and mean relative standard error of reproducibility tests (V1a–c).
TestParameterRelative Standard Error [%]Mean Relative Standard Error [%]
V1a–c Δ h / h o 2.2–11.76.2
V1a–c ξ 5.2–8.57.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Schalko, I.; Ruiz-Villanueva, V.; Maager, F.; Weitbrecht, V. Wood Retention at Inclined Bar Screens: Effect of Wood Characteristics on Backwater Rise and Bedload Transport. Water 2021, 13, 2231. https://doi.org/10.3390/w13162231

AMA Style

Schalko I, Ruiz-Villanueva V, Maager F, Weitbrecht V. Wood Retention at Inclined Bar Screens: Effect of Wood Characteristics on Backwater Rise and Bedload Transport. Water. 2021; 13(16):2231. https://doi.org/10.3390/w13162231

Chicago/Turabian Style

Schalko, Isabella, Virginia Ruiz-Villanueva, Fiona Maager, and Volker Weitbrecht. 2021. "Wood Retention at Inclined Bar Screens: Effect of Wood Characteristics on Backwater Rise and Bedload Transport" Water 13, no. 16: 2231. https://doi.org/10.3390/w13162231

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