<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Earth Sci.</journal-id>
<journal-title>Frontiers in Earth Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Earth Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-6463</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1357759</article-id>
<article-id pub-id-type="doi">10.3389/feart.2024.1357759</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Earth Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Modeling the morphological response of gravel&#x2013;bed rivers subject to repeated sediment bypass tunnel operations</article-title>
<alt-title alt-title-type="left-running-head">Facchini et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/feart.2024.1357759">10.3389/feart.2024.1357759</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Facchini</surname>
<given-names>Matteo</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2608154/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Vetsch</surname>
<given-names>David F.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2608041/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Boes</surname>
<given-names>Robert M.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1021899/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Siviglia</surname>
<given-names>Annunziato</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1417297/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Eastern Alps River Basin District Authority</institution>, <addr-line>Trento</addr-line>, <country>Italy</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Laboratory of Hydraulics</institution>, <institution>Hydrology and Glaciology (VAW)</institution>, <institution>ETH Zurich</institution>, <addr-line>Zurich</addr-line>, <country>Switzerland</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Department of Civil</institution>, <institution>Environmental and Mechanical Engineering</institution>, <institution>University of Trento</institution>, <addr-line>Trento</addr-line>, <country>Italy</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/92023/overview">Tianming Huang</ext-link>, Chinese Academy of Sciences (CAS), China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2624478/overview">Tetsuya Sumi</ext-link>, Kyoto University, Japan</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1175327/overview">Alessandro Stocchino</ext-link>, Hong Kong Polytechnic University, Hong Kong SAR, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2636913/overview">Claudia Ylla Arbos</ext-link>, Delft University of Technology, Netherlands</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: David F. Vetsch, <email>dvetsch@ethz.ch</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>05</day>
<month>04</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1357759</elocation-id>
<history>
<date date-type="received">
<day>19</day>
<month>12</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>15</day>
<month>03</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Facchini, Vetsch, Boes and Siviglia.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Facchini, Vetsch, Boes and Siviglia</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>
<bold>Introduction:</bold> Sediment bypass tunnels (SBTs) are operated during floods to re-establish sediment continuity past dams with the ultimate goal of decreasing accumulation in reservoirs. Depending on the location of the intake structure and operational conditions, SBTs may release either bed load-laden or bed load-free flows in the form of sudden floods (pulses). When evaluating the impacts of SBTs, the traditional approach has been to monitor channel changes in the years following project completion. However, by the time these impacts become evident, mitigation is challenging.</p>
<p>
<bold>Methods:</bold> This research adopts a forward-thinking methodology, emphasizing the prediction of potential impacts during the project&#x2019;s early planning stages. We use a one-dimensional morphodynamic model to forecast the potential morphological response of a gravel-bed river with an idealized geometry to a series of SBT operations characteristic of projects in the Swiss Alps.</p>
<p>
<bold>Results:</bold> The morphological response of the downstream reach over the medium to long term is influenced by the alternating pattern of bed load-laden and bed load-free flow releases. This pattern is instrumental in forming a transition region near the SBT outlet hydrograph boundary layer (HBL).</p>
<p>
<bold>Discussion:</bold> Within the (HBL), fluctuations in bed elevation and slope are localized, while variations in grain-size distribution persist throughout the downstream region.</p>
</abstract>
<kwd-group>
<kwd>sediment bypass tunnel</kwd>
<kwd>gravel-bed</kwd>
<kwd>sediment transport</kwd>
<kwd>1D numerical modeling</kwd>
<kwd>basement</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Hydrosphere</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Interruption of sediment continuity in gravel bed rivers due to dam construction causes reservoir sedimentation in the upstream reach and sediment starvation in the downstream reach. In turn, this affects downstream river morphology, resulting in channel narrowing and incision (<xref ref-type="bibr" rid="B72">Williams and Wolman, 1984</xref>; <xref ref-type="bibr" rid="B12">Brandt, 2000</xref>) and coarsening of the bed material (<xref ref-type="bibr" rid="B20">Dietrich et al., 1989</xref>), diminishing the ecological integrity of downstream rivers (<xref ref-type="bibr" rid="B34">Kondolf, 1997</xref>; <xref ref-type="bibr" rid="B52">Pitlick and Wilcock, 2001</xref>).</p>
<p>In the last decades, approaches for sediment management aiming at restoring natural sediment regimes have been widely discussed by the scientific community. Consequently, many sediment management techniques have been studied and implemented, whereby sediments are 1) mechanically removed and added downstream of dams (<xref ref-type="bibr" rid="B14">Bunte, 2004</xref>), 2) mobilized by induced side erosion caused by mechanical channel reconfiguration (<xref ref-type="bibr" rid="B19">Die Moran et al., 2013</xref>), 3) mobilized with artificial floods (<xref ref-type="bibr" rid="B37">Kondolf and Wilcock, 1996</xref>), or 4) routed around or through reservoirs (<xref ref-type="bibr" rid="B62">Sumi et al., 2012</xref>; <xref ref-type="bibr" rid="B26">Hager et al., 2020</xref>).</p>
<p>Implementation of such strategies significantly influences long-term changes in the downstream reach, leading to a variety of morphological configurations. Among them, sediment bypass tunnels (SBTs) have been designed with the main aim of routing sediment around reservoirs, thereby avoiding accumulation of new or even removing previous sediment (<xref ref-type="bibr" rid="B68">Vischer et al., 1997</xref>; <xref ref-type="bibr" rid="B35">Kondolf et al., 2014</xref>; <xref ref-type="bibr" rid="B8">Boes, 2015</xref>; <xref ref-type="bibr" rid="B61">Sumi, 2017</xref>; <xref ref-type="bibr" rid="B15">Chang and Lai, 2019</xref>; <xref ref-type="bibr" rid="B26">Hager et al., 2020</xref>). Investigations on existing SBTs showed that in some cases, up to 98% of the incoming sediments are diverted to the downstream river reach, considerably increasing the estimated lifespan of interested reservoirs (<xref ref-type="bibr" rid="B5">Auel et al., 2016</xref>; <xref ref-type="bibr" rid="B1">Albayrak et al., 2019</xref>; <xref ref-type="bibr" rid="B26">Hager et al., 2020</xref>). Most SBTs are located in mountainous regions, with small (&#x3c;1 million m<sup>3</sup>) to medium-sized (&#x3c;10 million m<sup>3</sup>) reservoirs, where a considerable amount of coarse material is entrained and transported as bed load (<xref ref-type="bibr" rid="B4">Auel and Boes, 2011</xref>; <xref ref-type="bibr" rid="B9">Boes et al., 2014</xref>). Typically, SBTs are solely operated during flood events with significant bed load transport (<xref ref-type="fig" rid="F1">Figure 1A</xref>), when the incoming bed load-laden flows are diverted to the downstream reach (see the planar view and cross-section of the Solis SBT in Switzerland in <xref ref-type="fig" rid="F1">Figure 1B</xref>). The frequency and duration of an SBT in operation is thus generally a function of the incoming flow. They are generally operated for a few hours, a few times a year (<xref ref-type="bibr" rid="B45">Mueller-Hagmann, 2017</xref>), imposing cycled hydrographs with or without synchronized sedimentographs as boundary conditions to the downstream reach. We distinguish two different types of SBTs (<xref ref-type="bibr" rid="B4">Auel and Boes, 2011</xref>; <xref ref-type="bibr" rid="B26">Hager et al., 2020</xref>), namely, type A with intakes at the reservoir head and type B with intakes positioned in the reservoir and therefore submerged, unless the reservoir level is substantially drawn down (see the 3D view of typical SBT systems in <xref ref-type="fig" rid="F1">Figure 1C</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Example and functioning of sediment bypass tunnels. <bold>(A)</bold> Outlet structure of the Solis SBT (canton of Grisons, Switzerland) during an SBT operation; <bold>(B)</bold> layout and cross-section of the Solis SBT (view against flow direction); <bold>(C)</bold> 3D view of type A (intake structure at the upstream end of the reservoir) and type B (intake structure within the reservoir) SBT systems.</p>
</caption>
<graphic xlink:href="feart-12-1357759-g001.tif"/>
</fig>
<p>While numerous studies have focused on the design, construction, and operation of SBTs (<xref ref-type="bibr" rid="B8">Boes, 2015</xref>; <xref ref-type="bibr" rid="B61">Sumi, 2017</xref>; <xref ref-type="bibr" rid="B15">Chang and Lai, 2019</xref>) and the ecological state of downstream reaches post-construction (<xref ref-type="bibr" rid="B6">Auel et al., 2017</xref>; <xref ref-type="bibr" rid="B43">Mart&#xed;n et al., 2017</xref>; <xref ref-type="bibr" rid="B55">Serrana et al., 2018</xref>), few have examined the morphological impacts of SBTs. In particular, research on the morphological responses of rivers to repeated SBT operations remains scarce. For instance, <xref ref-type="bibr" rid="B25">Fukuda et al. (2012)</xref> noted that Japan&#x2019;s Asahi Dam&#x2019;s downstream riffle&#x2013;pool morphology and grain-size distribution (GSD) were recovered after more than 15 years of SBT operations, yet highlighting the uncertainties of long-term effects of SBTs. Similarly, <xref ref-type="bibr" rid="B33">Kobayashi et al. (2018)</xref> showcased the potential of SBTs in mirroring upstream eco-morphology downstream based on in-field results. The current understanding of SBT impacts largely stems from monitoring post-SBT operations, but these are temporally and spatially constrained. Few studies, like <xref ref-type="bibr" rid="B22">Facchini (2018)</xref>, have adopted a sediment-focused, process-based approach to investigate these morphological effects. Given the rising global implementation of SBTs, their morphological impacts might become pronounced before detection. Modeling of these operations can shed light on the varied time-scaled responses of downstream river segments, offering a glimpse into both their immediate and long-term consequences.</p>
<p>Regardless of the type, SBTs exert a strong control on rates of water and sediment supply being released to the downstream reach. Essentially, SBTs act as boundary conditions for the downstream reach, where they impose particular cyclic flow and sediment regimes, mostly depending on 1) the SBT type (particularly in terms of their intake location, see <xref ref-type="fig" rid="F1">Figure 1C</xref>), 2) the magnitude of the event during which they are operated (i.e., the release scenario), and 3) other operational conditions related, for instance, to activities occurring in the upstream reach, such as gravel mining (<xref ref-type="bibr" rid="B7">Baumer and Radogna, 2015</xref>) or release of bed load-free flows (<xref ref-type="bibr" rid="B45">Mueller-Hagmann, 2017</xref>). Morphological effects in response to cycled hydrographs are well-documented and have been studied using laboratory experiments and by means of simplified mathematical models (<xref ref-type="bibr" rid="B75">Wong and Parker, 2006</xref>; <xref ref-type="bibr" rid="B50">Parker et al., 2008</xref>; <xref ref-type="bibr" rid="B3">An et al., 2017b</xref>). Adjustment of the bed surface size distribution of gravel-bed rivers in response to cycled hydrographs, similar to that of SBT operations, shows that the riverbed cyclically degrades and coarsens at a high flow and aggrades and becomes finer at a low flow (<xref ref-type="bibr" rid="B50">Parker et al., 2008</xref>). These variations of riverbed level and composition caused by repeated hydrographs under a constant sediment feed rate are found to be restricted to a short reach downstream of the feed point, i.e., a hydrograph boundary layer (HBL), where riverbed level and composition changes are confined. Downstream of the HBL, the variations are found in bed load composition and rate, while riverbed level and composition are invariant to the hydrograph (<xref ref-type="bibr" rid="B50">Parker et al., 2008</xref>). The dynamic related to the formation and evolution of the HBL described above, originally proposed by <xref ref-type="bibr" rid="B75">Wong and Parker (2006)</xref> and <xref ref-type="bibr" rid="B50">Parker et al. (2008)</xref> and confirmed by <xref ref-type="bibr" rid="B2">An et al. (2017a)</xref> breaks down if a sufficiently wide grain-size distribution is considered (<xref ref-type="bibr" rid="B3">An et al., 2017b</xref>). In this case, <xref ref-type="bibr" rid="B3">An et al. (2017b)</xref> demonstrates the persistence of forced grain-scale oscillations of the riverbed composition over the entire reach. Similarly, cyclic SBT operations might trigger the formation of an HBL. However, given the huge variability in the feed rate and composition caused by the particular operational conditions of the existing SBT operations (see Tables 2.1 and 2.2 in <xref ref-type="bibr" rid="B22">Facchini (2018)</xref>), it is unclear whether the very specific conditions needed for hydrograph boundary layer formation, i.e., the ones described by <xref ref-type="bibr" rid="B50">Parker et al. (2008)</xref>, can be matched by SBT operations or if these have the power to break down the mechanism of HBL formation. In this case, the effects of SBT operations, concerning in particular the riverbed composition, might be substantial even far from the feeding point (<xref ref-type="bibr" rid="B3">An et al., 2017b</xref>).</p>
<p>Herein, we use a numerical model to describe the potential impacts of a range of SBT release scenarios, which are typical of existing Swiss and Japanese projects, under different operational conditions such as riverbed slope and bed material grain size. We performed one-dimensional (1D) morphodynamic simulations of a river reach where bed load-laden or bed load-free flows are fed to the downstream reach in the form of repeated hydrographs and potentially combined sedimentographs. By changing the boundary conditions in our model to better resemble SBT release scenarios, we aim to produce a matrix of the predicted changes in the riverbed level and composition as a function of the imposed boundary conditions. While a classical geomorphologist or river engineer can qualitatively describe the investigated processes, only a modeling approach can provide quantitative predictions, thereby offering a more accurate assessment of the potential impacts of forthcoming SBT operations and thus inform operators to optimize the SBT operations. The framework produced can be further used to extend the scope of our modeling, e.g., to model in detail the effects of specific projects or to predict the effects of SBT releases with more advanced modeling techniques (e.g., two-dimensional).</p>
</sec>
<sec id="s2">
<title>2 Sediment bypass tunnel release scenarios</title>
<p>Morphological changes on the downstream reach are determined by SBT operations, i.e., the sequence of cycled hydrographs and potentially combined sedimentographs released by the SBT. In this section, a conceptual framework is developed to identify the possible scenarios under which SBTs can be operated. Each scenario is characterized by 1) the peak water discharge of the hydrograph and 2) the peak of the combined sedimentograph.</p>
<sec id="s2-1">
<title>2.1 Identification of possible release scenarios</title>
<p>The volumes of water and bed load conveyed through the SBT are linked to the discharge and bed load transport rate of the upstream (us) river reach. The relationship between the water discharge <italic>Q</italic>
<sub>
<italic>w</italic>
</sub> and the bed load transport rate <italic>Q</italic>
<sub>
<italic>b</italic>
</sub> (i.e., the bed load rating curve, BRC) in the upstream reach can be established under the hypothesis of uniform flow. Therefore, given the slope and the grain-size distribution of the upstream river reach, it is possible to determine the upstream bed load rating curve (line BRC<sub>us</sub> in <xref ref-type="fig" rid="F2">Figure 2</xref>). In general, SBTs are designed with a larger transport capacity than those of the upstream reach (<xref ref-type="bibr" rid="B4">Auel and Boes, 2011</xref>). Since SBTs are mostly operated under open channel flow conditions, the SBT bed load rating curve can be derived considering the slope of the SBT and the grain-size distribution of the sediment mixture entering the tunnel (line BRC<sub>SBT</sub> in <xref ref-type="fig" rid="F2">Figure 2</xref>).</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Sediment bypass tunnel release scenarios: bed load rating curves for the upstream reach (BRC<sub>us</sub>, solid blue line) and for the SBT (BRC<sub>SBT</sub>, solid red line). Black symbols represent the values of hydrograph and sedimentograph peaks used as input for the numerical runs under OC1, OC2a, and OC2b, and to each run, a specific symbol is given.</p>
</caption>
<graphic xlink:href="feart-12-1357759-g002.tif"/>
</fig>
<p>The minimum value of water discharge for which the SBT is first put into operation is <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>m</italic>
</sub>, which defines, on the BRC<sub>us</sub>, the corresponding minimum bed load discharge transported by the tunnel (<italic>Q</italic>
<sub>
<italic>b</italic>,<italic>m</italic>
</sub>). SBTs are designed according to a given water discharge capacity, i.e., a design water discharge <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>d</italic>
</sub>. On the BRC<sub>SBT</sub>, <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>d</italic>
</sub> identifies the maximum bed load discharge that can be carried by the SBT (<italic>Q</italic>
<sub>
<italic>b</italic>,<italic>M</italic>
</sub>). The <italic>Q</italic>
<sub>
<italic>w</italic>
</sub> needed to carry the maximum bed load discharge <italic>Q</italic>
<sub>
<italic>b</italic>,<italic>M</italic>
</sub> in the upstream reach is <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>M</italic>
</sub>. Then, we identified four possible scenarios based on the intensity of the flood events in the upstream reach causing (or not) the opening of the SBT. We define scenarios from I to IV in terms of water and bed load discharges released through the SBT to the downstream river reach (see <xref ref-type="fig" rid="F2">Figure 2</xref>) as follows:<list list-type="simple">
<list-item>
<p>1 For small events, i.e., water discharges in the upstream reach below the minimum discharge (<italic>Q</italic>
<sub>
<italic>w</italic>,<italic>us</italic>
</sub> &#x3c; <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>m</italic>
</sub>, point 1), we define <italic>scenario I</italic> for which the SBT is not operated and sediments carried by the upstream reach are not conveyed to the downstream reach.</p>
</list-item>
<list-item>
<p>2 For values of water discharge in the design range of the SBT (<italic>Q</italic>
<sub>
<italic>w</italic>,<italic>m</italic>
</sub> &#x2264; <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>us</italic>
</sub> &#x3c; <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>d</italic>
</sub>, points 2 and 3), we define <italic>scenario II</italic>, where the entire amount of sediments coming from upstream is conveyed through the SBT, and the possible SBT operations are identified by the points lying on the BRC<sub>us</sub> line (points 2 and 3).</p>
</list-item>
<list-item>
<p>3 In case of large floods exceeding the water discharge capacity of the SBT (<italic>Q</italic>
<sub>
<italic>w</italic>,<italic>d</italic>
</sub> &#x2264; <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>us</italic>
</sub> &#x2264; <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>M</italic>
</sub>, points 4&#x2013;10), we define <italic>scenario III</italic>, where the water and bed load discharges being delivered to the downstream reach range between the SBT design discharges and the maximum discharges possible. For scenario III, two extreme conditions can be considered. The first one occurs when the water discharge fed to the downstream reach flows only through the SBT, and thus it is kept constant and equal to its design discharge <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>d</italic>
</sub>. In this case, the water surplus (when <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>us</italic>
</sub> exceeds <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>d</italic>
</sub>) is stored in the reservoir, and the bed load discharge carried through the SBT ranges between the design and maximum one (vertical segment between points 4 and 6) since, for a given water discharge, the SBT can transport larger volumes of bed load as compared to the upstream reach, due to its larger transport capacity (<italic>BRC</italic>
<sub>
<italic>SBT</italic>
</sub> &#x3e; <italic>BRCus</italic>). The second condition occurs when the surplus water (when <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>us</italic>
</sub> exceeds <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>d</italic>
</sub>) flows inside the reservoir, but water is also released from the dam outlets. In this case, the delivered water and bed load discharges range both between the design and the maximum values, and possible release discharges can be identified on the BRC<sub>us</sub> (points 4, 7, and 8). Between these two situations, a number of other release scenarios are also possible, which can be identified generally inside the green area delimited by points 4, 6, and 8.</p>
</list-item>
<list-item>
<p>4 In case of very large floods exceeding the maximum water discharge needed to transport the maximum bed load discharge through the SBT (<italic>Q</italic>
<sub>
<italic>w</italic>,<italic>us</italic>
</sub> &#x3e; <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>M</italic>
</sub>), the bed load discharge fed to the downstream reach is constant and equal to the maximum transport capacity of the SBT (points 11 and 12). The water discharge exceeding the maximum discharge is released from the dam outlets.</p>
</list-item>
</list>
</p>
</sec>
<sec id="s2-2">
<title>2.2 Definition of possible operational conditions</title>
<p>The release scenarios presented in the previous section describe how SBTs can be operated based only on the intensity of flood events. However, SBTs can be operated in different ways (<xref ref-type="fig" rid="F2">Figure 2</xref>), i.e., the operational conditions (OCs) of the SBT depend on its design and on the transport conditions in the upstream reach. In this work, hydrographs and sedimentographs fed to the domain resemble real-event hydrographs, such as the ones relative to the Solis SBT (represented in <xref ref-type="fig" rid="F3">Figures 3A,B</xref>) with steep rising and falling limbs and a rather long release duration. The magnitude of the outflow hydrographs resembles that of natural floods up to a 10-year return period. Furthermore, we consider three different operational conditions, OC1, OC2, and OC3, among the most common adopted worldwide (<xref ref-type="bibr" rid="B61">Sumi, 2017</xref>). For OC1 and OC2a and b, the peak values of hydrographs and sedimentographs (if any) correspond to the numbered points in <xref ref-type="fig" rid="F2">Figure 2</xref>, while for OC3, the peak of the sedimentograph is reduced as specified hereafter.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Hydrographs and sedimentographs released through sediment bypass tunnels. <bold>(A)</bold> and <bold>(B)</bold> examples of hydrographs from real SBT operations at Solis, Switzerland, 13 August 2014 <bold>(A)</bold> and 11 June 2016 <bold>(B)</bold>. Solid and dashed lines refer to the water discharge conveyed through the SBT and the total water discharge delivered to the downstream reach, respectively. <bold>(C)</bold> SBT sediment release patterns (cycles) under OC1, OC2, and OC3. Under OC1 and OC3 at each SBT operation, bed load (red) is released (at different magnitudes) together with water (blue); under OC2a, bed load-laden and bed load-free flow releases are alternating; and under OC2b, three bed load-free flow releases follow a bed load-laden release.</p>
</caption>
<graphic xlink:href="feart-12-1357759-g003.tif"/>
</fig>
<sec id="s2-2-1">
<title>2.2.1 OC1: reference condition</title>
<p>The reference condition is relative to the conditions under which the SBT releases water and bed load during each operation. Thus, during each SBT operation, a hydrograph and a synchronized sedimentograph are released through the SBT. The grain-size distribution of the released bed load is identical to the distribution of the bed load transported in the upstream reach, and no bed load-free flow releases occur between two consecutive SBT operations (see <xref ref-type="fig" rid="F3">Figure 3C</xref>). This typically holds for type A SBTs (see <xref ref-type="fig" rid="F1">Figure 1C</xref>) operated during floods.</p>
</sec>
<sec id="s2-2-2">
<title>2.2.2 OC2: alternate bed load-laden and bed load-free flow releases</title>
<p>In practical situations, i.e., most commonly, SBTs may be used outside their reference (or design) conditions. Between two consecutive releases, certain SBTs may be opened just to route bed load-free flow (<xref ref-type="bibr" rid="B22">Facchini, 2018</xref>). To describe this case of alternating bed load-laden and bed load-free flow releases, we consider two release patterns for the bed load supply rate, one where bed load-laden and bed load-free flow releases are alternated (OC2a) and one where one bed load-laden flow release is followed by three bed load-free flow releases (OC2b) (see <xref ref-type="fig" rid="F3">Figure 3C</xref>).</p>
</sec>
<sec id="s2-2-3">
<title>2.2.3 OC3: upstream mining of coarse material</title>
<p>Gravel may be extracted upstream of the SBT intake structure, typically at the upstream end of the reservoir where backwater begins. This is the case for most Japanese SBTs, where, to avoid SBT invert abrasion (see <xref ref-type="bibr" rid="B34">Kondolf, 1997</xref>), coarse bed load is trapped and mined upstream of the SBT inlet structure and sold as a construction material (<xref ref-type="bibr" rid="B7">Baumer and Radogna, 2015</xref>). In this work, we consider that only the coarsest gravel fraction is removed, causing in turn a decrease in the volume of the released bed load (see <xref ref-type="fig" rid="F3">Figure 3C</xref>).</p>
</sec>
</sec>
</sec>
<sec id="s3">
<title>3 Numerical model</title>
<p>To predict the downstream changes in the riverbed slope and grain-size distribution, we run 1D numerical simulations using BASEMENT (<xref ref-type="bibr" rid="B18">Vetsch et al., 2017a</xref>). The hydrodynamic component of the model is solved using the Saint-Venant equations, while morphodynamic processes are described using the Exner equation, assuming a cohesionless bed composed of mixed-size sediment. The grain-size distribution on the riverbed surface and size stratification are described by the active layer modeling developed by <xref ref-type="bibr" rid="B29">Hirano (1971, 1972)</xref>. Further insights into the numerical solution and management of size stratification, including its vertical storage, in BASEMENT, are available in the software reference manual (<xref ref-type="bibr" rid="B66">Vetsch et al., 2017b</xref>). Our numerical investigations are restricted to the case of bed load transport of gravel with some admixtures of sand, typical for gravel-bed rivers (see <xref ref-type="sec" rid="s13">Supplementary Appendix A1.3</xref>). To account for the effect of sand (i.e., <italic>d</italic>
<sub>
<italic>k</italic>
</sub> &#x3c; 2 mm) on coarser grains, we adopt the <xref ref-type="bibr" rid="B70">Wilcock and Crowe (2003)</xref> transport formula (see <xref ref-type="sec" rid="s13">Supplementary Appendix A1.3</xref>).</p>
<p>The mathematical model resolved by BASEMENT is a standard model that has been widely used to describe stream-wise sorting and its relationship to bed elevation in gravel-bed rivers (<xref ref-type="bibr" rid="B49">Parker, 2008</xref>). It has been effective in describing the morphodynamic evolution in gravel-bed rivers, as demonstrated in laboratory experiments (<xref ref-type="bibr" rid="B16">Cui et al., 2003</xref>; <xref ref-type="bibr" rid="B47">Nelson et al., 2009</xref>; <xref ref-type="bibr" rid="B67">Viparelli et al., 2011</xref>; <xref ref-type="bibr" rid="B59">Stecca et al., 2016</xref>; <xref ref-type="bibr" rid="B2">An et al., 2017a</xref>). Consequently, it is expected to be well-suited for predicting potential morphodynamic changes in the downstream reach due to variations in sediment supply. Moreover, the lack of data relative to SBT-affected river reaches hinders the possibility of calibrating the model. However, the absence of a calibration does not affect the soundness of the results of this work or the possibility of extending them to specific SBT operations or to different modeling approaches. More details on the numerical modeling approaches are given in <xref ref-type="sec" rid="s13">Supplementary Appendix A1</xref>.</p>
</sec>
<sec id="s4">
<title>4 Definition of and input to the numerical model runs</title>
<p>We constructed our model using an idealized river geometry that captures the key characteristics of the Albula River, located just downstream of Tiefencastel in the canton of Grisons, Switzerland. Although the SBT operations of our model are hypothetical, they are based on the characteristics observed in the SBT of the Solis reservoir, which was constructed in 2012 (see <xref ref-type="fig" rid="F1">Figures 1A,B</xref>) (<xref ref-type="bibr" rid="B22">Facchini, 2018</xref>). A schematic representation of the modeled reach with input hydrographs and sedimentographs is given in <xref ref-type="fig" rid="F4">Figure 4</xref>.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Schematic representation of the modeled reach with input hydrographs (<italic>Q</italic>
<sub>
<italic>w</italic>
</sub>) and sedimentographs (<italic>Q</italic>
<sub>
<italic>b</italic>
</sub>).</p>
</caption>
<graphic xlink:href="feart-12-1357759-g004.tif"/>
</fig>
<sec id="s4-1">
<title>4.1 Upstream reach and SBT bed load rating curves</title>
<p>The bed load rating curve for the upstream reach (BRC<sub>us</sub> line in <xref ref-type="fig" rid="F2">Figure 2</xref>) is calculated considering uniform flow conditions in a rectangular channel 1.5% steep and 15 m wide where the friction term is calculated using a Gauckler&#x2013;Strickler parameter <italic>k</italic>
<sub>
<italic>s</italic>
</sub>&#x3d;32 m<sup>1/3</sup>s<sup>&#x2212;1</sup> obtained with <xref ref-type="sec" rid="s13">Supplementary Appendix Eq. 4</xref> considering d<sub>50</sub> &#x3d; 80 mm, relative to a sample from the Albula River upstream of the Solis reservoir. The riverbed grain-size distribution is similar to the that of the material sampled in the vicinity of the SBT inlet structure and is represented by a bimodal mixture with geometric mean size <italic>d</italic>
<sub>
<italic>g</italic>
</sub> &#x3d; 16.22 mm and geometric standard deviation <italic>&#x3c3;</italic>
<sub>
<italic>g</italic>
</sub> &#x3d; 7.37 mm. It is composed of 25% sand mode (d&#x3c;2 mm) and a 75% gravel mode (d &#x2265;2 mm) (<xref ref-type="bibr" rid="B23">Facchini et al., 2015</xref>). The full description of GSD is given in <xref ref-type="sec" rid="s13">Supplementary Tables A1, A2</xref>.</p>
<p>Usually, sediment transport calculations for SBTs differ from those for alluvial rivers because of low relative roughness (<italic>k</italic>
<sub>
<italic>e</italic>
</sub>/<italic>h</italic> &#x3c;&#x3c; 0.1, where <italic>k</italic>
<sub>
<italic>e</italic>
</sub> is the equivalent sand roughness height), a non-movable bed, and an initial planar bed. Therefore, the bed load rating curve for the SBT (BRC<sub>SBT</sub> line in <xref ref-type="fig" rid="F2">Figure 2</xref>) is calculated using the <xref ref-type="bibr" rid="B57">Smart and Jaeggi (1983)</xref> formula for the fixed planar bed (<xref ref-type="bibr" rid="B10">Boes et al., 2017</xref>) considering uniform open channel flow conditions and the characteristics of the Solis SBT, which has a slope of 1.9% and a channel width of 4.4 m. For the concrete invert, the roughness is characterized by a coefficient <italic>k</italic>
<sub>
<italic>s</italic>
</sub>&#x3d;65 m<sup>1/3</sup>s<sup>&#x2212;1</sup> (<xref ref-type="bibr" rid="B65">VAW, 2010</xref>). These conditions are used to evaluate the maximum bed load discharge carried by the SBT (<italic>Q</italic>
<sub>
<italic>b</italic>,<italic>M</italic>
</sub> in <xref ref-type="fig" rid="F2">Figure 2</xref>). The minimum (<italic>Q</italic>
<sub>
<italic>w</italic>,<italic>m</italic>
</sub>), the design (<italic>Q</italic>
<sub>
<italic>w</italic>,<italic>d</italic>
</sub>), and the maximum (<italic>Q</italic>
<sub>
<italic>w</italic>,<italic>M</italic>
</sub>) water discharges of the Solis SBT, considered in this work for running the numerical simulations, are 50, 170, and 275 m<sup>3</sup>s<sup>&#x2212;1</sup>, respectively, and the corresponding volumetric bed load rates (<italic>Q</italic>
<sub>
<italic>b</italic>,<italic>m</italic>
</sub>, <italic>Q</italic>
<sub>
<italic>b</italic>,<italic>d</italic>
</sub>, <italic>Q</italic>
<sub>
<italic>b</italic>,<italic>M</italic>
</sub>) are 0.18, 0.67, and 1.06 m<sup>3</sup>s<sup>&#x2212;1</sup>, respectively.</p>
</sec>
<sec id="s4-2">
<title>4.2 Downstream reach</title>
<p>Numerical simulations in the downstream reach are conducted under unsteady flow conditions in a straight rectangular channel with a width of 15 m, an initial constant slope of <italic>S</italic> &#x3d; 1.4%, and a Gauckler&#x2013;Strickler parameter of <italic>k</italic>
<sub>
<italic>s</italic>
</sub> &#x3d; 32 m<sup>1/3</sup>s<sup>&#x2212;1</sup>, as in the upstream reach. The total length of the channel is 15 km, even though for the results, only the first 10 km are considered to neglect the possible effects of the downstream boundary condition. The channel width is assumed to be sufficiently wide compared to the flow depth that sidewall effects can be neglected. Floodplains are not considered and deemed negligible due to the morphology of river stretches downstream of SBTs. The initial bed composition is the same as the upstream reach (see <xref ref-type="sec" rid="s13">Supplementary Tables A1</xref>). The numerical domain consists of cross-sections that are spaced out evenly (100 m), and the active layer thickness <italic>L</italic>
<sub>
<italic>a</italic>
</sub> has been considered constant in time and twice as thick as the <italic>d</italic>
<sub>90</sub> of its initial composition, i.e., considering the dimensionless constant of order 1 <italic>n</italic>
<sub>
<italic>a</italic>
</sub> &#x3d; 2 in <xref ref-type="sec" rid="s13">Supplementary Appendix Eq. 6</xref> (<xref ref-type="bibr" rid="B48">Parker, 2004</xref>). A schematic representation of the 1D channel used for the numerical simulations is given in <xref ref-type="fig" rid="F4">Figure 4</xref>. For all runs, initial conditions (<italic>t</italic> &#x3d; <italic>t</italic>
<sub>0</sub>) are obtained by setting constant discharge, i.e., uniform flow conditions both at the inlet and the outlet of the numerical domain and by running simulations under fixed-bed conditions until the steady state is reached.</p>
</sec>
<sec id="s4-3">
<title>4.3 Boundary conditions</title>
<p>Water and bed load are fed into the downstream reach at its inlet in the form of cycled hydrographs with or without synchronized sedimentographs. At the channel outlet (x &#x3d; 15 km), bed elevation is held constant throughout the calculation. The shape and duration of the hydrograph and sedimentograph and of the rising and falling limbs used as boundary conditions for the runs resemble the ones that occurred in 2014 and 2016 at the Solis SBT (<xref ref-type="fig" rid="F3">Figures 3A,B</xref>). Details of the SBT operations at Solis can be found in <xref ref-type="bibr" rid="B45">Mueller-Hagmann (2017)</xref>. The geometrical characteristics of the channel and the hydrograph features used in the numerical study are common to all runs and are represented in <xref ref-type="fig" rid="F4">Figure 4</xref> and summarized in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Summary of the geometrical characteristics of the upstream and downstream reach and of the hydrograph and sedimentograph characteristics used for the numerical study. The downstream reach values refer to the topography of the Albula River downstream of the Solis SBT (canton of Grisons, Switerland).</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Parameter</th>
<th align="left">Downstream</th>
<th align="left">SBT</th>
<th align="left">Upstream</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left"/>
<td align="left">Reach</td>
<td align="left"/>
<td align="left">Reach</td>
</tr>
<tr>
<td align="left">Total channel length [m]</td>
<td align="left">15,000</td>
<td align="left">&#x2a;</td>
<td align="left">&#x2a;</td>
</tr>
<tr>
<td align="left">Considered channel length <italic>L</italic> [m]</td>
<td align="left">10,000</td>
<td align="left">&#x2a;</td>
<td align="left">&#x2a;</td>
</tr>
<tr>
<td align="left">Channel width <italic>b</italic> [m]</td>
<td align="left">15</td>
<td align="left">4.4</td>
<td align="left">15</td>
</tr>
<tr>
<td align="left">Initial bed slope [%]</td>
<td align="left">1.4</td>
<td align="left">1.9</td>
<td align="left">1.5</td>
</tr>
<tr>
<td align="left">Strickler coefficient <italic>k</italic>
<sub>
<italic>s</italic>
</sub> [m<sup>1/3</sup>s<sup>&#x2212;1</sup>]</td>
<td align="left">32</td>
<td align="left">65</td>
<td align="left">32</td>
</tr>
<tr>
<td align="left">Hydrograph and sedimentograph duration [h]</td>
<td align="left">12</td>
<td align="left">&#x2a;</td>
<td align="left">&#x2a;</td>
</tr>
<tr>
<td align="left">Rising and falling limb duration [h]</td>
<td align="left">1</td>
<td align="left">&#x2a;</td>
<td align="left">&#x2a;</td>
</tr>
<tr>
<td align="left">Hydrograph lower discharge [m<sup>3</sup>s<sup>&#x2212;1</sup>]</td>
<td align="left">5</td>
<td align="left">&#x2a;</td>
<td align="left">&#x2a;</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>&#x2a;: Non relevant because the calculations have been developed with uniform open channel flow formulas.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Model runs are grouped into four sets based on the operational condition considered (i.e., OC1, OC2a, OC2b, and OC3). For each set, we run 12 numerical simulations numbered 1 to 12 (see the symbols in <xref ref-type="fig" rid="F2">Figure 2</xref>). For each run, the relative single release pattern (see <xref ref-type="fig" rid="F3">Figure 3C</xref>) is cyclically repeated. If any of the configurations is sustained for a sufficient amount of time (ca. 10<sup>3</sup> events for this study), a mobile-bed equilibrium is finally reached, which is not characterized by constant values of riverbed level and composition, but is instead dynamic. Thus, the equilibrium can strictly be defined only in relation to a cycle of events that is repeated over time (<xref ref-type="bibr" rid="B49">Parker et al., 2008</xref>). In this study, we refer to a cycle of four events, regardless of the OC, considered to encompass all operational conditions and to highlight the different dynamics developed when bed load-free flows are (OC2a and b) or are not (OC1 and OC3) released between two consecutive bed load-laden flow releases (see <xref ref-type="fig" rid="F3">Figure 3C</xref>).</p>
<p>For OC1, OC2a, and OC2b, the set of two peak values (<italic>Q</italic>
<sub>
<italic>w</italic>
</sub> and <italic>Q</italic>
<sub>
<italic>b</italic>
</sub>) are calculated as follows: 1) for run 1, <italic>Q</italic>
<sub>
<italic>b</italic>
</sub> is 0; 2) for runs 2, 3, 4, and 7, the peak values correspond to the ones lying on the upstream bed load rating curve; 3) simulations 5, 9, and 7 are runs using increasing values of the water discharge and using the same sediment discharge corresponding to the value found on the BRC<sub>us</sub> for a water discharge equal to the average of <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>d</italic>
</sub> and <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>M</italic>
</sub>; iv) runs 6, 10, 8, 11, and 12 are characterized by increasing values of water discharge and the same bed load rate equal to the maximum discharge that can be carried by the SBT (i.e.,.<italic>Q</italic>
<sub>
<italic>b</italic>,<italic>M</italic>
</sub>), which is calculated using the bed load rating curve of the SBT. The grain-size distribution of the material fed to the downstream reach is the mixture transported to the upstream reach and is calculated by running numerical simulations with constant discharge using the geometry of the upstream reach, i.e., 15 m width, 1.5% slope, and 32 m<sup>1/3</sup>s<sup>&#x2212;1</sup> Strickler <italic>k</italic>
<sub>
<italic>s</italic>
</sub> (see the <xref ref-type="sec" rid="s13">Supplementary Appendix A1.1</xref>). Peak values for OC3 runs are calculated considering that the coarsest fraction (<italic>d</italic> &#x3d; 256 mm) is removed from the sediment transport of the upstream reach before entering the SBT. As a consequence, the feeding GSD changes and becomes finer. We define two bimodal mixtures of sand and gravel used for our runs: one composed by 31% of the sand mode and 69% of the gravel mode (OC1 and OC2a and b) and one composed by 38% of the sand mode and 62% of the gravel mode (OC3). The peak values and the mean geometrical diameter used for the 4 &#xd7; 12 &#x3d; 48 numerical runs are given in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Parameters defining numerical model runs. <italic>Q</italic>
<sub>
<italic>w</italic>
</sub> &#x3d; hydrograph peak value; <italic>Q</italic>
<sub>
<italic>b</italic>
</sub> &#x3d; sedimentograph peak value; <italic>Q</italic>
<sub>
<italic>b</italic>
</sub>/<italic>Q</italic>
<sub>
<italic>b</italic>,<italic>us</italic>
</sub> &#x3d; ratio between the sedimentograph peak value and the bed load discharge in the upstream river reach under uniform flow conditions and water discharge <italic>Q</italic>
<sub>
<italic>w</italic>
</sub>; <italic>V</italic>
<sub>
<italic>b</italic>4</sub> &#x3d; bed load volume released over a cycle of four releases; <italic>d</italic>
<sub>
<italic>g</italic>,<italic>f</italic>
</sub> &#x3d; mean geometric diameter of the grain-size distribution of the feeding; <italic>d</italic>
<sub>
<italic>g</italic>,<italic>f</italic>
</sub>/<italic>d</italic>
<sub>
<italic>g</italic>
</sub> &#x3d; ratio between the mean geometric diameter of the grain-size distribution of the feeding and the upstream river reach.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th colspan="13" align="center">OC1, OC2a, and OC2b</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">run &#x23;</td>
<td align="left">1</td>
<td align="left">2</td>
<td align="left">3</td>
<td align="left">4</td>
<td align="left">5</td>
<td align="left">6</td>
<td align="left">7</td>
<td align="left">8</td>
<td align="left">9</td>
<td align="left">10</td>
<td align="left">11</td>
<td align="left">12</td>
</tr>
<tr>
<td align="left">
<italic>Q</italic>
<sub>
<italic>w</italic>
</sub> [<italic>m</italic>
<sup>3</sup>/<italic>s</italic>]</td>
<td align="left">30</td>
<td align="left">50</td>
<td align="left">100</td>
<td align="left">170</td>
<td align="left">170</td>
<td align="left">170</td>
<td align="left">223</td>
<td align="left">275</td>
<td align="left">197</td>
<td align="left">222</td>
<td align="left">428</td>
<td align="left">623</td>
</tr>
<tr>
<td align="left">
<italic>Q</italic>
<sub>
<italic>b</italic>
</sub> [<italic>m</italic>
<sup>3</sup>/<italic>s</italic>]</td>
<td align="left">0</td>
<td align="left">0.18</td>
<td align="left">0.39</td>
<td align="left">0.67</td>
<td align="left">0.87</td>
<td align="left">1.06</td>
<td align="left">0.87</td>
<td align="left">1.06</td>
<td align="left">0.87</td>
<td align="left">1.06</td>
<td align="left">1.06</td>
<td align="left">1.06</td>
</tr>
<tr>
<td align="left">
<italic>Q</italic>
<sub>
<italic>b</italic>
</sub>/<italic>Q</italic>
<sub>
<italic>b</italic>,<italic>us</italic>
</sub> [-]</td>
<td align="left">-</td>
<td align="left">1.0</td>
<td align="left">1.0</td>
<td align="left">1.0</td>
<td align="left">1.3</td>
<td align="left">1.58</td>
<td align="left">1.0</td>
<td align="left">1.0</td>
<td align="left">1.13</td>
<td align="left">1.22</td>
<td align="left">0.68</td>
<td align="left">0.50</td>
</tr>
<tr>
<td align="left">&#x2a;<italic>V</italic>
<sub>
<italic>b</italic>4</sub> &#x22c5; 10<sup>3</sup> [<italic>m</italic>
<sup>3</sup>]</td>
<td align="left">-</td>
<td align="left">28</td>
<td align="left">61.1</td>
<td align="left">106</td>
<td align="left">138</td>
<td align="left">167</td>
<td align="left">138</td>
<td align="left">167</td>
<td align="left">138</td>
<td align="left">167</td>
<td align="left">167</td>
<td align="left">167</td>
</tr>
<tr>
<td align="left">
<italic>d</italic>
<sub>
<italic>g</italic>,<italic>f</italic>
</sub> [<italic>mm</italic>]</td>
<td align="left">-</td>
<td align="left">7.6</td>
<td align="left">8.8</td>
<td align="left">9.6</td>
<td align="left">10.0</td>
<td align="left">10.3</td>
<td align="left">10.0</td>
<td align="left">10.3</td>
<td align="left">10.0</td>
<td align="left">10.3</td>
<td align="left">10.8</td>
<td align="left">11.2</td>
</tr>
<tr>
<td align="left">
<italic>d</italic>
<sub>
<italic>g</italic>,<italic>f</italic>
</sub>/<italic>d</italic>
<sub>
<italic>g</italic>
</sub> [-]</td>
<td align="left">-</td>
<td align="left">0.47</td>
<td align="left">0.54</td>
<td align="left">0.59</td>
<td align="left">0.62</td>
<td align="left">0.64</td>
<td align="left">0.62</td>
<td align="left">0.64</td>
<td align="left">0.62</td>
<td align="left">0.64</td>
<td align="left">0.67</td>
<td align="left">0.69</td>
</tr>
<tr>
<td colspan="13" align="center">OC3</td>
</tr>
<tr>
<td align="left">run &#x23;</td>
<td align="left">1</td>
<td align="left">2</td>
<td align="left">3</td>
<td align="left">4</td>
<td align="left">5</td>
<td align="left">6</td>
<td align="left">7</td>
<td align="left">8</td>
<td align="left">9</td>
<td align="left">10</td>
<td align="left">11</td>
<td align="left">12</td>
</tr>
<tr>
<td align="left">
<italic>Q</italic>
<sub>
<italic>w</italic>
</sub> [<italic>m</italic>
<sup>3</sup>/<italic>s</italic>]</td>
<td align="left">30</td>
<td align="left">50</td>
<td align="left">100</td>
<td align="left">170</td>
<td align="left">170</td>
<td align="left">170</td>
<td align="left">223</td>
<td align="left">275</td>
<td align="left">197</td>
<td align="left">222</td>
<td align="left">428</td>
<td align="left">623</td>
</tr>
<tr>
<td align="left">
<italic>Q</italic>
<sub>
<italic>b</italic>
</sub> [<italic>m</italic>
<sup>3</sup>/<italic>s</italic>]</td>
<td align="left">0</td>
<td align="left">0.02</td>
<td align="left">0.05</td>
<td align="left">0.08</td>
<td align="left">0.11</td>
<td align="left">0.13</td>
<td align="left">0.11</td>
<td align="left">0.13</td>
<td align="left">0.11</td>
<td align="left">0.13</td>
<td align="left">0.13</td>
<td align="left">0.13</td>
</tr>
<tr>
<td align="left">
<italic>Q</italic>
<sub>
<italic>b</italic>
</sub>/<italic>Q</italic>
<sub>
<italic>b</italic>,<italic>us</italic>
</sub> [-]</td>
<td align="left">-</td>
<td align="left">0.11</td>
<td align="left">0.13</td>
<td align="left">0.12</td>
<td align="left">0.16</td>
<td align="left">0.19</td>
<td align="left">0.13</td>
<td align="left">0.12</td>
<td align="left">0.14</td>
<td align="left">0.15</td>
<td align="left">0.08</td>
<td align="left">0.06</td>
</tr>
<tr>
<td align="left">
<italic>V</italic>
<sub>
<italic>b</italic>4</sub> &#x22c5; 10<sup>3</sup> [<italic>m</italic>
<sup>3</sup>]</td>
<td align="left">-</td>
<td align="left">3.5</td>
<td align="left">7.64</td>
<td align="left">13.2</td>
<td align="left">17.2</td>
<td align="left">20.9</td>
<td align="left">17.2</td>
<td align="left">20.9</td>
<td align="left">17.2</td>
<td align="left">20.9</td>
<td align="left">20.9</td>
<td align="left">20.9</td>
</tr>
<tr>
<td align="left">
<italic>d</italic>
<sub>
<italic>g</italic>,<italic>f</italic>
</sub> [<italic>mm</italic>]</td>
<td align="left">-</td>
<td align="left">4.8</td>
<td align="left">5.4</td>
<td align="left">5.9</td>
<td align="left">6.1</td>
<td align="left">6.2</td>
<td align="left">6.1</td>
<td align="left">6.2</td>
<td align="left">6.1</td>
<td align="left">6.2</td>
<td align="left">6.5</td>
<td align="left">6.7</td>
</tr>
<tr>
<td align="left">
<italic>d</italic>
<sub>
<italic>g</italic>,<italic>f</italic>
</sub>/<italic>d</italic>
<sub>
<italic>g</italic>
</sub> [-]</td>
<td align="left">-</td>
<td align="left">0.30</td>
<td align="left">0.33</td>
<td align="left">0.36</td>
<td align="left">0.38</td>
<td align="left">0.38</td>
<td align="left">0.38</td>
<td align="left">0.38</td>
<td align="left">0.38</td>
<td align="left">0.38</td>
<td align="left">0.40</td>
<td align="left">0.41</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>&#x2a;: Volumes are relative to OC1, used here as a reference. OC2a and b volumes are half and a quarter of the reference, respectively.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
</sec>
<sec sec-type="results" id="s5">
<title>5 Results</title>
<p>Our results are presented and later discussed in terms of longitudinal riverbed changes (channel slope and riverbed level) and evolution of the active layer composition. As a reference configuration, we assume the slope and the grain-size distribution of the reach upstream of the reservoir. In particular, to assess the ability of SBT operations in reproducing the upstream morphological conditions in the downstream reach, we define two dimensionless variables. Namely, we define a dimensionless slope as <italic>S</italic>
<sup>&#x2a;</sup> &#x3d; <italic>S</italic>/<italic>S</italic>
<sub>
<italic>us</italic>
</sub> and a dimensionless mean geometric size as <inline-formula id="inf1">
<mml:math id="m1">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>A</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>, where <italic>S</italic> is the riverbed slope and <italic>d</italic>
<sub>
<italic>g</italic>
</sub> is the geometric mean of the grain-size distribution in the downstream reach. The subscripts <sub>
<italic>us</italic>
</sub> and <sub>
<italic>AL</italic>
</sub> stand for &#x201c;upstream&#x201d; (i.e., reference) and &#x201c;active layer&#x201d; (i.e., actual dg in the active layer of the simulated downstream reach), respectively. Where not specifically related to the upstream reach with the relative subscript, the variables refer to the downstream reach. <italic>S</italic>
<sup>&#x2a;</sup> and <inline-formula id="inf2">
<mml:math id="m2">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula> &#x3d; 1 indicate values of slope and riverbed composition in the downstream reach equal to the upstream ones, respectively, while <italic>S</italic>
<sup>&#x2a;</sup> &#x3d; 0.93 indicates the initial condition. Thus, if <italic>S</italic>
<sup>&#x2a;</sup> is smaller than 0.93, incision occurs, and if it is larger, deposition occurs. Similarly, if <italic>S</italic>
<sup>&#x2a;</sup> and <inline-formula id="inf3">
<mml:math id="m3">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula> are smaller (larger) than 1, the resulting slopes and riverbed composition are milder (steeper) and finer (coarser) than the upstream ones. To analyze the changes in the riverbed level, referring to the longitudinal profile of the channel, we define a deviatoric riverbed level (or deviatoric elevation) as &#x394;<italic>&#x3b7;</italic>(<italic>x</italic>, <italic>t</italic>) &#x3d; <italic>&#x3b7;</italic>(<italic>x</italic>, <italic>t</italic>) &#x2212; <italic>&#x3b7;</italic>(<italic>x</italic>, 0). Negative (positive) values of &#x394;<italic>&#x3b7;</italic> indicate incision (aggradation), while &#x394;<italic>&#x3b7;</italic> &#x3d; 0 represents the initial condition.</p>
<sec id="s5-1">
<title>5.1 Riverbed slope and composition at mobile-bed equilibrium</title>
<p>The results in <xref ref-type="fig" rid="F5">Figure 5</xref> show that the values of the riverbed slope and composition at the mobile-bed equilibrium are highly dependent on the volume of water that is released through the SBT, i.e., on the unbalance between the transport capacity in the downstream reach and the bed load feeding from the SBT. The results in <xref ref-type="fig" rid="F5">Figure 5</xref> are displayed in terms of dimensionless slope <italic>S</italic>
<sup>&#x2a;</sup> and mean geometric sizes <inline-formula id="inf4">
<mml:math id="m4">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula>, as a function of the dimensionless water discharge <inline-formula id="inf5">
<mml:math id="m5">
<mml:msubsup>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>, where <italic>Q</italic>
<sub>
<italic>w</italic>,<italic>d</italic>
</sub> is the design water discharge of the SBT. They are relative to the configuration reached after the dynamic mobile-bed equilibrium has been attained (after thousands to tens of thousands of releases) and refer precisely to the end of a cycle of four events (see <xref ref-type="fig" rid="F3">Figure 3C</xref>). Since both the downstream riverbed level and grain-size distribution are not uniformly distributed along the reach, we present <italic>S</italic>
<sup>&#x2a;</sup> and <inline-formula id="inf6">
<mml:math id="m6">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula> results relative to the cross-section located 10 km downstream of the feeding point, chosen because it represents the location least affected by the disturbance caused by both boundary conditions. Furthermore, at mobile-bed equilibrium, the dynamics is, to a high degree of approximation, everywhere the same outside a short inlet region. In that region, the dynamics is more affected by the cyclic fluctuations imposed by repeated hydrographs combined or not with synchronized sedimentographs than outside. In the following sections, we discuss the results that will support this choice.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Long-term results concerning dimensionless riverbed slope <italic>S</italic>
<sup>&#x2a;</sup> (panels <bold>(A, C, E,</bold> and <bold>G)</bold>) and composition <inline-formula id="inf7">
<mml:math id="m7">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula> (panels <bold>(B, D, F,</bold> and <bold>H)</bold>) at mobile-bed equilibrium considering four different SBT operational conditions, i.e., OC1, OC2a and b, and OC3, plotted against a non-dimensional water discharge relative to the SBT design discharge, i.e., <inline-formula id="inf8">
<mml:math id="m8">
<mml:msubsup>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>. Riverbed slopes and compositions are calculated 10 km downstream of the feeding point. The dashed line in panels <bold>(A)</bold>, <bold>(C)</bold>, <bold>(E)</bold>, and <bold>(G)</bold> indicates the initial condition of the runs (<italic>S</italic>
<sup>&#x2a;</sup>&#x3d;0.93). Dotted lines separate different scenarios, and in panels <bold>(A)</bold>, <bold>(C)</bold>, <bold>(E)</bold>, and <bold>(G)</bold>, dashed lines connect results obtained with runs having the same feeding rate.</p>
</caption>
<graphic xlink:href="feart-12-1357759-g005.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F5">Figure 5</xref> shows that, for events where water is released only from the dam outlets and the SBT is not in operation (scenario I: no bed load at the inlet), the riverbed slope and composition at the mobile-bed equilibrium are smaller (<italic>S</italic>
<sup>&#x2a;</sup> &#x3d; 1.0, <xref ref-type="fig" rid="F5">Figures 5A,C,E,G</xref>) and larger (<inline-formula id="inf9">
<mml:math id="m9">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula>&#x3d; 6.6, <xref ref-type="fig" rid="F5">Figures 5B,D,F,H</xref>) than the reference, respectively. Furthermore, since <italic>S</italic>
<sup>&#x2a;</sup> is smaller than the initial value (i.e., <italic>S</italic>
<sup>&#x2a;</sup> &#x3d; 0.93, dashed horizontal line in <xref ref-type="fig" rid="F5">Figures 5A,C,E,G</xref>), riverbed incision occurs. This behavior is similar for all operational conditions. If the SBT is operational (scenarios II, III, and IV, see <xref ref-type="fig" rid="F2">Figure 2</xref>), there is high variability among equilibrium slopes (<xref ref-type="fig" rid="F5">Figures 5A,C,E,G</xref>) and compositions (<xref ref-type="fig" rid="F5">Figures 5B,D,F,H</xref>). For the runs with the same <italic>Q</italic>
<sub>
<italic>b</italic>
</sub> (scenarios III and IV), <italic>S</italic>
<sup>&#x2a;</sup> decreases with increasing <italic>Q</italic>
<sub>
<italic>w</italic>
</sub> (dashed blue lines in the <xref ref-type="fig" rid="F5">Figures 5A,C,E,G</xref>). Under OC1 (<xref ref-type="fig" rid="F5">Figure 5A</xref>), <italic>S</italic>
<sup>&#x2a;</sup>~&#x3e;1.0 if the bed load supply is larger than the transport capacity (runs 5, 6, 9, and 10). If it is equal to the transport capacity (runs 2, 3, 4, 7, and 8), <italic>S</italic>
<sup>&#x2a;</sup> &#x3d; 1.0, and if it is smaller (runs 11 and 12), <italic>S</italic>
<sup>&#x2a;</sup>; &#x3c;1.0. Reductions in the volume of the total bed load delivered by the SBT due to operational conditions (see values of <italic>V</italic>
<sub>
<italic>b</italic>4</sub> in <xref ref-type="table" rid="T2">Table 2</xref>) cause a further generalized decrease in the equilibrium slopes (<xref ref-type="fig" rid="F5">Figures 5C,E,G</xref>); that is, for each scenario under OC2a, OC2b, and OC3, there are relevant reductions in the equilibrium slope, with <italic>S</italic>
<sup>&#x2a;</sup> values even smaller than 0.2 in the most extreme cases (e.g., scenario IV in <xref ref-type="fig" rid="F5">Figure 5G</xref>). Concerning the riverbed surface composition (excluding scenario I), the results show that at mobile-bed equilibrium under OC1 (<xref ref-type="fig" rid="F5">Figure 5B</xref>), the downstream riverbed has the same composition as the upstream one. On the contrary, bed load-free releases (OC2a and b) and reduced sediment input (OC3) shift the downstream riverbed composition away from the upstream one (<xref ref-type="fig" rid="F5">Figures 5D,F,H</xref>). In the first two cases (<xref ref-type="fig" rid="F5">Figures 5D,F</xref>), the more bed load-free flow is released from the SBT, the more material from the downstream riverbed surface is entrained in the bed load transport. As a consequence, mostly coarser particles are left on the downstream riverbed surface, causing its composition to become coarser as compared to the upstream one. In the other case (<xref ref-type="fig" rid="F5">Figure 5H</xref>), releasing finer material and reducing the released volume significantly induces a remarkable fining of the riverbed material, mostly due to the replacement of the initial composition of the downstream riverbed surface with the feeding one, lacking the coarsest grain size.</p>
<p>Our equilibrium results are in good agreement with the widely applied, qualitative relationship for channel adjustments proposed by <xref ref-type="bibr" rid="B39">Lane (1955)</xref>, which can be extended to infer that gravel-bed rivers adjust prevailing channel characteristics (e.g., slope and riverbed composition) to provide, with the available water discharge, just the velocity required for the transportation of the load supplied from upstream (<xref ref-type="bibr" rid="B41">Mackin, 1948</xref>). This suggests that our model supports the expected physics of a river reach subject to changes in water and bed load supply regimes.</p>
</sec>
<sec id="s5-2">
<title>5.2 Riverbed level and composition fluctuations at the inlet of the downstream reach</title>
<p>The results in <xref ref-type="fig" rid="F6">Figure 6</xref> show that both riverbed level and composition might fluctuate, depending on the operational conditions. The first fluctuates in a short region close to the SBT outlet, i.e., at the upstream end of the considered river reach, while the second fluctuates along the whole reach. For representing the general trend of both riverbed level and composition, we average the results at the mobile-bed equilibrium over a four-event cycle (see <xref ref-type="fig" rid="F3">Figure 3C</xref>), regardless of the operational conditions. For representing the fluctuations, we calculate their minimum and maximum as well as the 25th and 75th percentiles of the distribution over the whole cycle. We refer to the equilibrium characterized by values averaged over a cycle as <italic>macroscopic</italic> mobile-bed equilibrium (<xref ref-type="bibr" rid="B3">An et al., 2017b</xref>). Results in <xref ref-type="fig" rid="F6">Figure 6</xref> show that under OC2a and b, both riverbed level and composition present larger average values close to the upstream end (black solid lines in <xref ref-type="fig" rid="F6">Figures 6C</xref>&#x2013;<xref ref-type="fig" rid="F6">F</xref>). The main cause of the pronounced deposition (and consequent coarsening) occurring during alternate bed load-laden&#x2013;bed load-free flow releases is the imbalance between bed load supply and transport capacity evaluated over a four-event cycle of SBT operations. During the releases of bed load-laden flow, the riverbed aggrades, and during the following clear-water releases, the aggraded material is transported downstream, returning the morphology of the channel back to the equilibrium state at the end of each cycle. The more bed load-free flow is released after a bed load-laden flow release, the larger the deposition close to the inlet is. Furthermore, over a cycle of events, both riverbed level and composition may or may not fluctuate, depending on the conditions under which the SBT is operated. As for the average values, if the feeding rate and transport capacity are significantly unbalanced and bed load-free flows are released between two consecutive bed load-laden flow releases (OC2a and OC2b), fluctuations in riverbed level (<xref ref-type="fig" rid="F6">Figures 6C,E</xref>) and composition (<xref ref-type="fig" rid="F6">Figures 6D,F</xref>) are clearly observable. On the contrary, if bed load-free flows are never released (OC1 and OC3), results show no significant fluctuations in the deviatoric riverbed level (<xref ref-type="fig" rid="F6">Figures 6A,G</xref>) and small fluctuations of the riverbed composition (<xref ref-type="fig" rid="F6">Figures 6B,H</xref>). Under all circumstances, the disturbance of the average values and the fluctuations of the riverbed level are always confined within the first 1 to 2 km downstream of the feeding point (see <xref ref-type="fig" rid="F6">Figures 6C,E</xref>), i.e., they are limited to the HBL (<xref ref-type="bibr" rid="B50">Parker et al., 2008</xref>). In addition to the findings of <xref ref-type="bibr" rid="B50">Parker et al. (2008)</xref>, our results indicate that oscillations of the riverbed composition extend beyond the HBL (see <xref ref-type="fig" rid="F6">Figures 6D,F</xref>) due to the different upstream boundary conditions imposed in this work, that is, the original idea of the HBL breaks down in this case, mostly due to the fact that alternating bed load-laden and bed load-free flow releases generate an imbalance between sediment feeding and transport capacity, which propagates the oscillations of the riverbed composition along the reach. The oscillations are anyway more intense inside the HBL due to the fact that the riverbed level, and thus the slope, is also strongly oscillating.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Long-term fluctuations of deviatoric bed elevation &#x394;<italic>&#x3b7;</italic> [m] (left) and geometric mean size of riverbed composition <italic>d</italic>
<sub>
<italic>g</italic>
</sub> during a cycle of four events for run 3 (Scenario II). Panels <bold>(A)</bold> and <bold>(B)</bold> refer to OC1, <bold>(C)</bold> and <bold>(D)</bold> refer to OC2a, <bold>(E)</bold> and <bold>(F)</bold> refer to OC2b, and <bold>(G)</bold> and <bold>(H)</bold> refer to OC3. In all panels, the black line represents the average over the cycle (macroscopic equilibrium), the light gray area represents the minimum&#x2013;maximum oscillation, and the dark gray area represents the oscillation between the 25% and 75% percentiles. In panels <bold>(B)</bold>, <bold>(D)</bold>, <bold>(F)</bold>, and <bold>(H)</bold>, the black dashed line represents the geometric mean size of the feeding.</p>
</caption>
<graphic xlink:href="feart-12-1357759-g006.tif"/>
</fig>
<p>The results in <xref ref-type="fig" rid="F7">Figure 7</xref> expand the analysis on the effect of alternating bed load-laden and bed load-free flow releases on the dynamics of the armoring ratio. The latter fluctuates more intensely close to the feeding point (<xref ref-type="fig" rid="F7">Figures 7A,C,E,G</xref>) and if bed load-free flows are released between two consecutive bed load-laden flow releases (<xref ref-type="fig" rid="F7">Figures 7C,E</xref>). In these cases, fluctuations of the armoring ratio follow a pattern, which follows, in turn, the release cycle of bed load-laden and bed load-free flows from the SBT (refer to <xref ref-type="fig" rid="F3">Figure 3C</xref> for the duration of the cycles under different operational conditions). In <xref ref-type="fig" rid="F7">Figure 7</xref>, two cross sections are considered, namely, 1 km downstream of the feeding point (i.e., &#x201c;upstream end,&#x201d; <xref ref-type="fig" rid="F7">Figures 7A,C,E,G</xref>) and 10 km downstream of the feeding point (i.e., &#x201c;downstream end,&#x201d; <xref ref-type="fig" rid="F7">Figures 7B,D,F,H</xref>). The dynamics of fining and coarsening is different depending on the operational conditions, i.e., if no, one or more bed load-free flow releases are carried out between two consecutive bed load-laden flow releases, that is, both at the upstream and downstream end of the domain, there are no significant fluctuations under OC1 and OC3 (<xref ref-type="fig" rid="F7">Figures 7A,B,G,H</xref>, respectively). However, a small increase in the armoring ratio occurs at the upstream end during the rising limb of the hydrograph (hours 0 to 1 in the inset in <xref ref-type="fig" rid="F7">Figure 7A</xref>), followed by a fining during the first hour of the peak (hours 1 to 2 in the inset in <xref ref-type="fig" rid="F7">Figure 7A</xref>). Then, the armoring ratio remains constant (hours 2 to 11 in the inset in <xref ref-type="fig" rid="F7">Figure 7A</xref>), followed by coarsening during the falling limb of the hydrograph (hours 11 to 12 in the inset in <xref ref-type="fig" rid="F7">Figure 7A</xref>). At the downstream end, the coarsening&#x2013;fining process is delayed by about 5 h (inset in <xref ref-type="fig" rid="F7">Figure 7B</xref>), which is the time needed by the perturbations in the average diameter of the mixture to propagate until the downstream end. As previously described in the literature (<xref ref-type="bibr" rid="B54">Ribberink, 1987</xref>; <xref ref-type="bibr" rid="B58">Stecca et al., 2014</xref>; <xref ref-type="bibr" rid="B56">Siviglia et al., 2017</xref>), such perturbations, which are associated with changes in the grain-size distribution in the active layer, travel slower than hydraulic perturbations. In this case, the flow velocity is 4.5 ms<sup>&#x2212;1</sup>, i.e., it takes ca. 35 min for the hydraulic perturbation to travel to the downstream end of the reach.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Long-term evolution of the armoring ratio during a cycle of four events at equilibrium. Panels <bold>(A)</bold> and <bold>(B)</bold> refer to OC1, and the two insets are a zoom in to the first event. Panels <bold>(C)</bold> and <bold>(D)</bold> refer to OC2a, <bold>(E)</bold> and <bold>(F)</bold> to OC2b, and <bold>(G)</bold> and <bold>(H)</bold> to OC3. Panels <bold>(A)</bold>, <bold>(C)</bold>, <bold>(E)</bold>, and <bold>(G)</bold> refer to a cross-section 1 km downstream of the feeding point (upstream end), while panels <bold>(B)</bold>, <bold>(D)</bold>, <bold>(F)</bold>, and <bold>(H)</bold> refer to a cross-section 10 km downstream of the feeding point (downstream end). Both the armoring ratios during bed load-laden (solid lines) and bed load-free (dashed lines) water releases are plotted. Black dashed vertical lines define the beginning and the end of the peak of the hydrograph, while the black dash&#x2013;dotted vertical line divides two consecutive events. Arrows indicate the positive marching of time during each release, i.e., releases 3 and 4 also span over 24 h like releases 1 and 2. They are represented like this to highlight their cyclic pattern.</p>
</caption>
<graphic xlink:href="feart-12-1357759-g007.tif"/>
</fig>
<p>Close to the inlet, under OC2a and b, the armoring is removed within the first hour (rising limb of the hydrograph, hours 0 to 1 in <xref ref-type="fig" rid="F7">Figures 7C,E</xref>) of the bed load-laden flow release (solid lines in <xref ref-type="fig" rid="F7">Figures 7C,E</xref>), and then coarsening occurs during bedload-free flow releases (dashed lines in <xref ref-type="fig" rid="F7">Figures 7C,E</xref>). Under OC2a, two identical cycles composed of two events follow one another, while under OC2b, coarsening occurs during the three bed load-free flow releases following the first bed load-laden release. However, more than 80% of the coarsening process under OC2b occurs during the first bed load-free flow release, where the armoring ratio increases from 1.06 to 2.72. Similarly, the armoring ratio increases from 1.04 to 2.54 during the bed load-free flow release under OC2a.</p>
</sec>
<sec id="s5-3">
<title>5.3 Riverbed level and composition fluctuations for different scenarios at the inlet of the downstream reach</title>
<p>Results in <xref ref-type="fig" rid="F8">Figure 8</xref> refer to all scenarios and show that the riverbed composition evolution relative to scenario II presented in <xref ref-type="fig" rid="F6">Figures 6</xref>, <xref ref-type="fig" rid="F7">7</xref> is similar in each run, i.e., results at mobile-bed equilibrium appear to be similar regardless of the scenario. In <xref ref-type="fig" rid="F8">Figure 8</xref>, the results are plotted for the armoring ratio <inline-formula id="inf10">
<mml:math id="m10">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula> <italic>versus</italic> the dimensionless water discharge <inline-formula id="inf11">
<mml:math id="m11">
<mml:msubsup>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula> and are presented in terms of the average values of the riverbed slope and composition calculated over a cycle of four releases. Two cross-sections are considered here, identical to <xref ref-type="fig" rid="F7">Figure 7</xref>. Moreover, to quantify the oscillation of the riverbed composition during a cycle, we represent the range between the 25th and the 75th percentile of its distribution with whiskers (similar to the dark gray area in <xref ref-type="fig" rid="F6">Figure 6</xref>). The armoring ratio <inline-formula id="inf12">
<mml:math id="m12">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula> lies, under almost all circumstances, inside the region of unarmored conditions, as defined in section A1.5, both upstream (<xref ref-type="fig" rid="F8">Figures 8A,C,E,G</xref>) and downstream (<xref ref-type="fig" rid="F8">Figures 8B,D,F,H</xref>), with the only exception being the upstream end composition of the riverbed under OC2b (<xref ref-type="fig" rid="F8">Figure 8E</xref>). Fluctuations of the armoring ratio (represented by the whiskers) occur only for OC2a and OC2b since they are caused by the imbalance between the bed load supply and the transport capacity. Similarly, similar fluctuations are observed for each scenario also for the riverbed level [results not shown, see Section 5.4 in <xref ref-type="bibr" rid="B22">Facchini (2018)</xref>]. Differently from the results presented in <xref ref-type="fig" rid="F6">Figure 6</xref>, the fluctuations of the armoring ratio are restricted to the HBL since <inline-formula id="inf13">
<mml:math id="m13">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula> incorporates both oscillations of the active layer and substrate composition, while in <xref ref-type="fig" rid="F6">Figure 6,</xref> only the active layer composition is represented. Moreover, the more bed load-free flow is released (OC2b), the higher the imbalance and thus the more armored the riverbed surface results, that is, while the mean values of <inline-formula id="inf14">
<mml:math id="m14">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula> under OC2a are below 2.0 for each scenario, they are all larger than 2.0 under OC2b. Results in <xref ref-type="fig" rid="F8">Figure 8</xref> also show that under OC1 (<xref ref-type="fig" rid="F8">Figure 8A</xref>) and OC3 (<xref ref-type="fig" rid="F8">Figure 8G</xref>), there are no fluctuations in the riverbed composition due to the balance between bed load supply and transport capacity. Under these operational conditions, the armoring ratio is close to 1.0 (unarmored conditions), with the only exception of OC3 at the downstream end of the domain.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Long-term results concerning the armoring ratio <inline-formula id="inf15">
<mml:math id="m15">
<mml:msubsup>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula> at mobile-bed equilibrium under OC1 (panels <bold>(A,B)</bold>), OC2a (panels <bold>(C,D)</bold>), OC2b (panels <bold>(E, F)</bold>), and OC3 (panels <bold>(G,H)</bold>), plotted against a non-dimensional water discharge relative to the SBT design discharge, i.e., <inline-formula id="inf16">
<mml:math id="m16">
<mml:msubsup>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>. Panels <bold>(A)</bold>, <bold>(C)</bold>, <bold>(E)</bold>, and <bold>(G)</bold> refer to a cross-section 1 km downstream of the feeding point (upstream end), while panels <bold>(B)</bold>, <bold>(D)</bold>, <bold>(F)</bold>, and <bold>(H)</bold> refer to a cross-section 10 km downstream of the feeding point (downstream end). Dotted lines separate different scenarios.</p>
</caption>
<graphic xlink:href="feart-12-1357759-g008.tif"/>
</fig>
<p>With respect to the magnitude and extension of the riverbed fluctuations, we represent their amplitude and length in <xref ref-type="fig" rid="F9">Figure 9</xref>, in terms of dimensionless variables, i.e., <italic>A</italic>
<sup>&#x2a;</sup> and <italic>L</italic>
<sup>&#x2a;</sup> as defined in section A1.7.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Long-term results concerning the dimensionless hydrograph boundary layer (HBL) amplitude panels <bold>(A,B)</bold> and length panels <bold>(C,D)</bold> under OC2a and OC2b at mobile-bed equilibrium plotted against the dimensionless water discharge <inline-formula id="inf17">
<mml:math id="m17">
<mml:msubsup>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula>.</p>
</caption>
<graphic xlink:href="feart-12-1357759-g009.tif"/>
</fig>
<p>The aggradation (i.e., positive <italic>A</italic>
<sub>HBL</sub>) induced by SBT operations at the mobile-bed equilibrium is between 5 and 35 times the <italic>d</italic>
<sub>90</sub> (<xref ref-type="fig" rid="F9">Figures 9A,B</xref>). Despite the large amplitude of the upstream aggradation caused by the release of bed load-laden flow (as for results in <xref ref-type="fig" rid="F6">Figure 6</xref>), this remains confined ca. within the first 1.5 km downstream of the feeding point for all scenarios (see <xref ref-type="fig" rid="F9">Figures 9C,D</xref>).</p>
</sec>
<sec id="s5-4">
<title>5.4 Short-term effects of SBT releases</title>
<p>Results in <xref ref-type="fig" rid="F10">Figure 10</xref> are relative to a cycle of four SBT releases after 12 operations, i.e., after three complete cycles, under OC2a and OC2b. As compared to the results in <xref ref-type="fig" rid="F6">Figures 6C,E</xref>, the results in <xref ref-type="fig" rid="F10">Figures 10A,C</xref> show that the riverbed level is oscillating during a cycle close to the upstream end, while it is adapting to the boundary conditions imposed by the SBT (incision occurs). The general trend of erosion is confirmed by results in <xref ref-type="fig" rid="F5">Figure 5</xref>, where initial values of the slope are equal to 1.4% for each run (S<sup>&#x2a;</sup> &#x3d; 0.93, dashed line in <xref ref-type="fig" rid="F5">Figures 5C,E</xref>), while the equilibrium slopes for run 3 (upward-pointing triangle) are 1% and 0.7% for OC2a and b (S<sup>&#x2a;</sup> &#x3d; 0.67 and 0.47), respectively. Results in <xref ref-type="fig" rid="F10">Figures 10B,D</xref> are quite similar to the results in <xref ref-type="fig" rid="F6">Figures 6D,F</xref>, showing that the evolution of the riverbed composition is faster than the one of the riverbed level, that is, the riverbed composition is close to the equilibrium configuration and it oscillates around the average (solid black line in <xref ref-type="fig" rid="F10">Figures 10B,D</xref>) along the whole reach already after a few SBT operations.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Short-term fluctuations of deviatoric bed elevation &#x394;<italic>&#x3b7;</italic> [m] (left) and geometric mean size of riverbed composition <italic>d</italic>
<sub>
<italic>g</italic>
</sub> [mm] during a cycle of four events for run 3 (Scenario II), after the first three cycles (i.e., 12 events). Panels <bold>(A)</bold> and <bold>(B)</bold> refer to OC2a and <bold>(C)</bold> and <bold>(D)</bold> refer to OC2b. In all panels, the black line represents the average over the cycle (macroscopic equilibrium), the light gray area represents the minimum&#x2013;maximum oscillation, and the dark gray area the oscillation between the 25% and 75% percentiles. In panels <bold>(B)</bold> and <bold>(D)</bold>, the black dashed line represents the geometric mean size of the feeding.</p>
</caption>
<graphic xlink:href="feart-12-1357759-g010.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="discussion" id="s6">
<title>6 Discussion</title>
<sec id="s6-1">
<title>6.1 Long-term effects of SBT operations on riverbed slope and composition</title>
<p>SBTs release water, and potentially bedload, in the form of cyclic hydrographs combined or not with synchronized sedimentographs. The characteristics of the hydrograph, i.e., flow magnitude, duration, shape, and sequence, are relevant to the bedload flux only in terms of their contribution to the total transported volume (<xref ref-type="bibr" rid="B51">Phillips et al., 2018</xref>), but they greatly affect the river morphology (riverbed level and composition). Compared to less intense hydrographs with more gradual rising and falling limbs (<xref ref-type="bibr" rid="B31">Humphries et al., 2012</xref>), intense hydrographs with steep rising and falling limbs may affect the morphology far from the inlet and cause coarsening of the riverbed (e.g., <xref ref-type="bibr" rid="B42">Mao, 2012</xref>). In addition to that, the limitation of sediment supply has been proven to be a first-order control on the formation of an armored layer (<xref ref-type="bibr" rid="B27">Hassan et al., 2006</xref>), like long periods of low flow between floods, which favor the consolidation and stabilization of the riverbed (<xref ref-type="bibr" rid="B53">Reid et al., 1985</xref>; <xref ref-type="bibr" rid="B69">Waters and Curran, 2015</xref>). Similarly, short periods between two consecutive floods reduce the time over which fine sediment can be entrained and transported as bed load, thus favoring downstream fining (<xref ref-type="bibr" rid="B24">Ferrer-Boix and Hassan, 2015</xref>).</p>
<p>In this work, we fed the domain with hydrographs having steep rising and falling limbs as compared to the duration of the peak and no low-flow periods between two consecutive hydrographs. Furthermore, we limit the sediment feeding only during bed load-free flow releases. Results in <xref ref-type="fig" rid="F5">Figure 5</xref> show that if the feeding grain-size distribution is the same as the upstream one (this is always the case, except for OC3) and the feed rate is the one of the upstream bedload rating curves (runs 2, 3, 4, 7, and 8 under OC1 lying on the BRC<sub>us</sub> in <xref ref-type="fig" rid="F2">Figure 2</xref>), the SBT can be considered a means to make the reservoir <italic>transparent</italic> since the volume and composition of the sediment delivered downstream correspond exactly to what is transported by the upstream reach. In this case, the upstream and downstream reaches are perfectly connected by the SBT, and the resulting riverbed composition and slopes will be the same as the reference SBT. If the feed rate or grain-size distribution are different from those described above, i.e., if the scenario or the operational conditions are changed, the volume and composition of the sediment delivered downstream do not correspond to what is transported upstream; thus, the reservoir acts as a <italic>disturbance</italic> which cannot be countered by the SBT. As a consequence, the downstream reach will adapt to the modified conditions, and its mobile-bed equilibrium will diverge from the reference.</p>
<p>Bed load-free flow releases play a major role in creating unbalanced conditions. Under OC2a (alternating bed load-laden and bed load-free releases) and OC2b (one bed load-laden release followed by three bed load-free releases), the equilibrium slopes and riverbed composition are milder and coarser than the reference ones (see <xref ref-type="fig" rid="F5">Figures 5C&#x2013;F</xref>), respectively. This is a consequence of the bed load-free flow releases during which the sediment is picked up from the riverbed and entrained into the bed load (see <xref ref-type="fig" rid="F7">Figure 7</xref>). Under OC3, the volume of the delivered sediment is further reduced, and the feeding grain-size distribution is deprived of its coarsest fractions (refer to <italic>V</italic>
<sub>
<italic>b</italic>4</sub> and <italic>d</italic>
<sub>
<italic>g</italic>,<italic>f</italic>
</sub> in <xref ref-type="table" rid="T2">Table 2</xref>); thus, reductions in the resulting slope and riverbed composition are to be expected (see <xref ref-type="fig" rid="F5">Figures 5G,H</xref>).</p>
<p>Mobile-bed equilibrium conditions develop after sustaining an operational SBT regime for enough time, regardless of the initial conditions, i.e., of the downstream reach morphological state at the beginning of the process. The mobile-bed equilibrium developed after sustaining an SBT operational regime for a sufficient amount of time is a fluctuating equilibrium, which consists of a state in which the bed degrades and coarsens at high-flow stages and during bed load-free releases, and aggrades and becomes finer at the low-flow stages or during bed load-laden releases (<xref ref-type="fig" rid="F6">Figures 6</xref>, <xref ref-type="fig" rid="F7">7</xref>). A cycle can be identified depending on the boundary conditions sustained to reach the equilibrium, at the end of which the riverbed level and composition are the same as at the beginning. In this study, the cycle is characterized by a varying number of events to represent a pattern, depending on the operational conditions. For OC1 and OC3, a cycle is represented by one event; for OC2a and OC2b, it is represented by two and four events, respectively (see <xref ref-type="fig" rid="F3">Figure 3B</xref>). Given the different equilibrium slopes, the amplitude of the riverbed level changes, i.e., the amplitude of the hydrograph boundary layer (<xref ref-type="fig" rid="F9">Figure 9</xref>), is different depending on the operational conditions. During a single bed load-laden release, the release of the same volume of sediment causes larger depositions for milder slopes (OC2a and OC2b, <xref ref-type="fig" rid="F5">Figures 5C,E</xref>) than for steeper slopes (OC1, <xref ref-type="fig" rid="F5">Figure 5A</xref>). Furthermore, during the first bed load-laden release under OC2a and OC2b, the bed aggrades and becomes finer, reaching its maximum level (upper limit of the light gray area in <xref ref-type="fig" rid="F6">Figures 6C,E</xref>) and its finest riverbed composition (lower limit of the light gray area in <xref ref-type="fig" rid="F6">Figures 6D,F</xref>). During the following one (OC2a) or three (OC2b) bed load-free releases, the riverbed gradually degrades and becomes coarser, reaching the initial riverbed level (&#x394;<italic>&#x3b7;</italic> &#x3d; 0 in <xref ref-type="fig" rid="F6">Figures 6C,E</xref>) and composition (upper limit of the light gray area in <xref ref-type="fig" rid="F6">Figures 6D,F</xref>). Under OC1 and OC3, the whole aggradation&#x2013;fining and degradation&#x2013;coarsening dynamics is completed within one event, i.e., during a single hydrograph, as shown in <xref ref-type="bibr" rid="B50">Parker et al. (2008)</xref>, <xref ref-type="bibr" rid="B75">Wong and Parker (2006)</xref>, and <xref ref-type="bibr" rid="B3">An et al. (2017b)</xref> since every release is equal to the previous and following releases.</p>
</sec>
<sec id="s6-2">
<title>6.2 Hydrograph boundary layer formation, rework of the riverbed composition, and short-term effects</title>
<p>On a short temporal scale (say a few cycles), results reported in <xref ref-type="fig" rid="F10">Figure 10</xref> indicate that even far from the mobile-bed equilibrium conditions, riverbed slope and grain-size distribution cyclically fluctuate during every single-cycle SBT operation. This means that the mobile-bed equilibrium is slowly attained through a repeated sequence of hydrographs and sedimentographs, which at every cycle of SBT operations results in cyclic variations of both riverbed slope and grain-size distribution.</p>
<p>As shown in <xref ref-type="fig" rid="F6">Figures 6</xref>, <xref ref-type="fig" rid="F7">7</xref>, oscillations of riverbed level and composition are mostly induced by the presence of bed load-free flow releases. The oscillations over a single bed load-laden release, i.e., over the single hydrograph (<xref ref-type="fig" rid="F7">Figures 7A,B</xref>), are negligible and occur mostly because the sedimentograph is synchronized with the hydrograph (see <xref ref-type="fig" rid="F3">Figure 3</xref>) and not with the time distribution of the transport capacity, that is, the sediment discharge relative to the hydrographs lags behind the feeding rate (<xref ref-type="bibr" rid="B22">Facchini, 2018</xref>). Such a time lag between the end of the rising limb of the hydrograph and of the sedimentograph has been observed under unsteady flows (<xref ref-type="bibr" rid="B11">Bombar et al., 2011</xref>; <xref ref-type="bibr" rid="B42">Mao, 2012</xref>; <xref ref-type="bibr" rid="B32">Khosravi et al., 2019</xref>). Fluctuations of riverbed levels shown in <xref ref-type="fig" rid="F6">Figure 6</xref> are realized only in a relatively short region, i.e., the hydrograph boundary layer (<xref ref-type="bibr" rid="B50">Parker et al., 2008</xref>; <xref ref-type="bibr" rid="B3">An et al., 2017b</xref>), downstream of the feed point. Unlike the findings of <xref ref-type="bibr" rid="B50">Parker et al. (2008)</xref>, here oscillations of the riverbed composition extend along the whole reach and are not limited to the HBL (see <xref ref-type="fig" rid="F6">Figures 6D,F</xref>). Thus, the original idea of the HBL breaks down in this case, mostly due to the fact that alternating bed load-laden and bed load-free flow releases propagate the oscillations of the riverbed composition along the whole reach. The grain-size distribution oscillations are anyway more intense inside the hydrograph boundary layer due to the fact that the riverbed level, and thus the slope, is strongly oscillating inside it. Bed load-free releases produce fluctuations of riverbed level and composition during each event, even after the first few SBT operations (see <xref ref-type="fig" rid="F10">Figure 10</xref>), and induce the system to evolve toward the equilibrium configuration, repeating the same cycle over and over. Results in <xref ref-type="fig" rid="F10">Figures 10A,C</xref> show that the amplitude (ca. 3 m) as well as the length of the hydrograph boundary layer (less than 2 km) are comparable to the ones shown in <xref ref-type="fig" rid="F6">Figures 6C,E</xref>. However, the average deviatoric riverbed level (solid black line in <xref ref-type="fig" rid="F10">Figures 10A,C</xref>) shows an erosion trend that confirms the general trajectory predicted by the results shown in <xref ref-type="fig" rid="F5">Figure 5</xref>, where to all runs under OC2a and b (<xref ref-type="fig" rid="F5">Figures 5C,E</xref>), an equilibrium slope is associated, which is milder than the initial runs (compare the results with the dashed line in <xref ref-type="fig" rid="F5">Figure 5</xref>). The fluctuations of riverbed composition caused by bed load-free flow releases are the same after few events (<xref ref-type="fig" rid="F10">Figures 10B,D</xref>) or at equilibrium (<xref ref-type="fig" rid="F6">Figures 6D,F</xref>), that is, riverbed composition adapts to the imposed boundary conditions faster than the riverbed level.</p>
<p>The cycle over which the fluctuations around the equilibrium state occur is represented in <xref ref-type="fig" rid="F7">Figure 7</xref>, where it is shown how the system returns to the equilibrium conditions after each bed load-laden flow release under OC1 (<xref ref-type="fig" rid="F7">Figures 7A,B</xref>) and OC3 (<xref ref-type="fig" rid="F7">Figures 7G,H</xref>), after two releases (alternate bed load-laden and bed load-free) under OC2a (<xref ref-type="fig" rid="F7">Figures 7C,D</xref>), and after four releases (one bed load-laden followed by three bed load-free) under OC2b (<xref ref-type="fig" rid="F7">Figures 7E,F</xref>), that is, the reworking of the riverbed surface is negligible under OC1 and OC3, where hydrographs and sedimentographs are always synchronized, while it is not under OC2a and b, where bed load-free releases cause an imbalance between transport capacity and sediment feeding. The dynamics imposed by bedload-free releases is repeated from the first few events on (not shown here), which confirms that the riverbed composition is reworked along the whole reach, starting with the initial SBT operations. The trend of the riverbed composition close to the upstream end (<xref ref-type="fig" rid="F7">Figures 7C&#x2013;F</xref>) is in line with the one relative to the bed load rate during flood events described in the literature (<xref ref-type="bibr" rid="B42">Mao, 2012</xref>). That is, during the initial stages of the rising limb, mostly fine particles are entrained, while coarse particles of the feeding and those on the riverbed surface are not, hence causing coarsening. During the following stages of the rising limb, the transport capacity increases together with the capacity of entrained coarse particles into bed load, leaving finer particles on the riverbed surface, hence causing fining. On the contrary, during the falling limb and the following bed load-free water release(s), coarse particles are deposited first, causing the coarsening of the riverbed.</p>
</sec>
<sec id="s6-3">
<title>6.3 Implications for river management</title>
<p>The construction of SBTs is among the few techniques applied to restore water and sediment continuity without removing man-made obstructions. Projects restoring natural sediment and water regimes in rivers aim to promote morphological complexity, hydraulic diversity, and ecosystem health (<xref ref-type="bibr" rid="B40">Lisle, 1989</xref>; <xref ref-type="bibr" rid="B44">Merz and Ochikubo Chan, 2005</xref>; <xref ref-type="bibr" rid="B21">Elosegi et al., 2010</xref>; <xref ref-type="bibr" rid="B17">Czapiga et al., 2022</xref>). Studies on the ecological effects of SBT operations suggest that micro- and macroinvertebrate density and richness (<xref ref-type="bibr" rid="B6">Auel et al., 2017</xref>; <xref ref-type="bibr" rid="B43">Mart&#xed;n et al., 2017</xref>; <xref ref-type="bibr" rid="B55">Serrana et al., 2018</xref>) as well as fish habitat suitability (<xref ref-type="bibr" rid="B62">Sumi et al., 2012</xref>) are positively affected by SBT operations. This notwithstanding, the outcome of such operations is strongly dependent on the frequency of the events (<xref ref-type="bibr" rid="B43">Mart&#xed;n et al., 2017</xref>; <xref ref-type="bibr" rid="B55">Serrana et al., 2018</xref>). Extreme magnitudes or frequencies of SBT operations could have severe consequences on the river ecosystem (<xref ref-type="bibr" rid="B43">Mart&#xed;n et al., 2017</xref>). However, since SBTs are operated only few times a year (<xref ref-type="bibr" rid="B45">Mueller-Hagmann, 2017</xref>; <xref ref-type="bibr" rid="B61">Sumi, 2017</xref>), they impose an army schedule on the downstream reach: &#x201c;long periods of relative quiet punctuated by short bursts of frenetic activity&#x201d; (<xref ref-type="bibr" rid="B73">Wohl, 2013</xref>). That is, when SBTs are not releasing bed load-laden or bed load-free flows to the downstream reach, the latter is fed only with the environmental flow and no sediment (<xref ref-type="bibr" rid="B43">Mart&#xed;n et al., 2017</xref>; <xref ref-type="bibr" rid="B55">Serrana et al., 2018</xref>) and can be considered morphologically inactive during non-operational times (not counting potential feeding from tributaries). A release pattern is repeated every year or so when SBTs are operated, but when they are not, the downstream reach keeps the morphological configuration reached at the end of the SBT operation. Our study confirms the importance of the frequency of SBT operations and proves additionally that the effects of SBT operations are strongly dependent on how (not only how often) an SBT is operated, i.e., on the operational conditions and scenarios. Results in <xref ref-type="fig" rid="F5">Figures 5</xref>&#x2013;<xref ref-type="fig" rid="F10">10</xref> show that releasing bed load-free flows can induce the formation of an HBL and the rework of the riverbed composition in the downstream reach even at the early stages. The riverbed level oscillates in a region close to the upstream end, while the riverbed composition is reworked along the entire reach (see <xref ref-type="fig" rid="F6">Figure 6</xref>), even after the first few SBT operations (see <xref ref-type="fig" rid="F10">Figure 10</xref>). Except for OC2b, the average armoring ratio at equilibrium is always smaller than 2, i.e., SBT operations cause the riverbed to be in unarmored conditions at equilibrium (<xref ref-type="bibr" rid="B27">Hassan et al., 2006</xref>). Furthermore, riverbed composition changes occur, almost in all cases, at the event scale (<xref ref-type="fig" rid="F7">Figure 7</xref>) and are not greatly affected by the initial state of the riverbed surface (as discussed in the <xref ref-type="sec" rid="s13">Supplementary Appendix A1.2</xref>). Improved sediment mobility and the presence of gravel and cobbles on the riverbed have tremendous ecological importance as a habitat for benthic macroinvertebrates and fish (<xref ref-type="bibr" rid="B36">Kondolf and Matthews, 1991</xref>; <xref ref-type="bibr" rid="B38">Kondolf et al., 2008</xref>; <xref ref-type="bibr" rid="B28">Hauer, 2015</xref>).</p>
<p>Results in <xref ref-type="fig" rid="F7">Figure 7</xref> show that the downstream riverbed composition oscillates during a cycle of bed load-free and bed load-laden releases, causing the riverbed to be either in unarmored or in mobile-armor conditions. That is, depending on whether the last SBT release is bed load-laden or bed load-free, the resulting riverbed composition is different (see <xref ref-type="fig" rid="F7">Figures 7C,E</xref>). Moreover, SBTs are operating only a few times a year and are therefore unable to restore natural hydrological processes (<xref ref-type="bibr" rid="B74">Wohl et al., 2015</xref>) since during non-operational times, no sediment and environmental flows are fed to the downstream reach (<xref ref-type="bibr" rid="B43">Mart&#xed;n et al., 2017</xref>; <xref ref-type="bibr" rid="B55">Serrana et al., 2018</xref>). Thus, knowing the pattern of the SBT releases allows predicting the riverbed composition in the downstream reach, which may last for several months during non-operational times.</p>
<p>Our analysis demonstrates the critical role of managing SBT operations in shaping the downstream river morphology. Specifically, to reduce incision, or rather favor deposition to revert river damming effects, prioritizing bed load-laden water releases over bed load-free ones is essential. Bed load-free releases facilitate the transport of sediment deposited near the SBT outlet, accelerating the reworking of the riverbed texture, with a minimal impact on riverbed elevation. Moreover, bedload-free water releases favor the formation of the HBL, which causes strong oscillations of riverbed level and composition in a confined region. This may have negative effects on 1) infrastructures due to the large variation in the riverbed level (<xref ref-type="bibr" rid="B63">Surian and Rinaldi, 2003</xref>), 2) fish due to stranding (<xref ref-type="bibr" rid="B46">Nagrodski et al., 2012</xref>), and 3) macroinvertebrates due to drift (<xref ref-type="bibr" rid="B13">Bruno et al., 2013</xref>). To avoid such strong oscillations and counteract the formation of the HBL, bedload-laden releases should be favored against bedload-free releases. However, the release of only bedload-laden waters (OC1 and OC3) can also have negative effects such as an increase in the downstream slope due to large floods (scenario III) where the SBT is operated above the transport capacity of the downstream river reach. Overall, coupling bedload-laden releases with controlled water releases may be the most effective strategy to reduce sediment deposition and promote downstream transport.</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s7">
<title>7 Conclusion</title>
<p>The rapid increase in SBT development worldwide has outpaced our understanding of their potential morphological effects. Although SBTs are commonly assumed to be beneficial in restoring sediment continuity past dams, this assumption has not been rigorously evaluated. The primary method of gaining an insight into the effects of SBTs has been through long-term field monitoring. However, due to the duration of time it takes for these effects to be observed in the field, we may be experiencing the full effects of SBT dams before they are fully understood. To proactively understand the geomorphic impacts of SBTs, we modeled the effects on grain-size distribution and bed elevation associated with typical (still idealized) SBT release scenarios in the Alps. Using a 1-D morphodynamic model, we simulated different SBT release scenarios based on a simplified geometry resembling the reach of the Albula River downstream of the Solis SBT. Our model results can be summarized as follows:<list list-type="simple">
<list-item>
<p>&#x2022; The more bed load-free flow is released, the milder the riverbed slopes become, and the coarser their composition turns out to be. If the released material&#x2019;s grain-size distribution becomes finer in relation to a reduction in the released bedload volume, the downstream riverbed composition becomes finer, contrasting with the coarser upstream composition.</p>
</list-item>
<list-item>
<p>&#x2022; Repetitive SBT operations at mobile-bed equilibrium yield cyclic riverbed changes. The riverbed undergoes aggradation and fining during bedload-laden releases and degradation and coarsening with bed load-free releases. Riverbed level oscillations mainly occur upstream near the inlet, forming a hydrograph boundary layer (HBL). The HBL&#x2019;s longitudinal extension and amplitude strongly depend on the frequency of bed load-free releases. However, riverbed grain-size distribution changes span the entire downstream length with the strongest intensity within the HBL.</p>
</list-item>
<list-item>
<p>&#x2022; In the short term, the riverbed displays oscillations after only a few SBT operations.</p>
</list-item>
<list-item>
<p>&#x2022; Our insights particularly pertain to type-B SBTs, where bed load-free releases are more prevalent than in type-A SBTs.</p>
</list-item>
</list>
</p>
<p>Although an experienced geomorphologist could qualitatively describe most of the processes happening in the downstream river reach of an SBT, our quantitative, objective framework offers a way to evaluate the potential influences of upcoming SBT developments. Present SBT project monitoring, due to its limited scope, faces challenges in providing significant insights into the detailed description of the morphological response of the river downstream from SBT. While our method should not replace regular monitoring of the morphological impact, it aids in formulating testable hypotheses for more targeted monitoring.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s8">
<title>Data availability statement</title>
<p>The data can be made available by the first author upon reasonable request.</p>
</sec>
<sec id="s9">
<title>Author contributions</title>
<p>MF: conceptualization, data curation, formal analysis, investigation, methodology, writing review&#x2013;editing. DV: funding acquisition, project administration, resources, software, validation, and writing&#x2013;review and editing. RB: funding acquisition, project administration, resources, validation, and writing&#x2013;review and editing. AS: conceptualization, methodology, project administration, resources, software, supervision, validation, writing&#x2013;original draft, and writing&#x2013;review and editing.</p>
</sec>
<sec sec-type="funding-information" id="s10">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was co-funded by the Swiss Federal Office for the Environment (FOEN) and open access funding by ETH Zurich.</p>
</sec>
<ack>
<p>The authors would like to thank the Swiss Federal Office for the Environment (FOEN) for their financial support in the scope of the project &#x201c;Sediment transport and habitat dynamics&#x201d; within the research program &#x201c;Hydraulic Engineering and Ecology,&#x201d; and the Swiss Competence Center for Energy Research&#x2014;Supply of Electricity (SCCER-SoE), to which the project was further affiliated. They extend their gratitude to C. Oertli and T. Ziegler of the electric power company of Zurich (Elektrizit&#xe4;tswerk der Stadt Z&#xfc;rich, ewz) and D. Rickenmann of the Swiss Federal Institute for Forest, Snow, and Landscape Research (WSL) for providing data relative to the Solis SBT and the upstream reach of the Albula River, respectively. They would also like to thank M. Saletti of the Department of Fisheries and Oceans (Government of Canada) for accurate proofreading and the valuable comments that greatly helped them improve the paper. Eventually, they kindly acknowledge the three reviewers for their careful reading of the manuscript and their many insightful comments and suggestions that helped them in improving the paper, too.</p>
</ack>
<sec sec-type="COI-statement" id="s11">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
<p>The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.</p>
</sec>
<sec sec-type="disclaimer" id="s12">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s13">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/feart.2024.1357759/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/feart.2024.1357759/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet1.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Albayrak</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Mueller-Hagmann</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2019</year>). &#x201c;<article-title>Efficiency evaluation of swiss sediment bypass tunnels</article-title>,&#x201d; in <conf-name>Proceedings of 3rd International Workshop on Sediment Bypass Tunnels</conf-name>. <pub-id pub-id-type="doi">10.3929/ethz-b-000337945</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>An</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Cui</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Parker</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2017a</year>). <article-title>gravel-bed river evolution in earthquake-prone regions subject to cycled hydrographs and repeated sediment pulses</article-title>. <source>Earth Surf. Process. Landforms</source> <volume>42</volume>, <fpage>2426</fpage>&#x2013;<lpage>2438</lpage>. <pub-id pub-id-type="doi">10.1002/esp.4195</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>An</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Parker</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2017b</year>). <article-title>Effect of grain sorting on gravel-bed river evolution subject to cycled hydrographs: bedload sheets and breakdown of the hydrograph boundary layer</article-title>. <source>J. Geophys. Res. Earth Surf.</source> <volume>122</volume>, <fpage>1513</fpage>&#x2013;<lpage>1533</lpage>. <pub-id pub-id-type="doi">10.1002/2016JF003994</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Auel</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2011</year>). &#x201c;<article-title>Sediment bypass tunnel design &#x2014; review and outlook</article-title>,&#x201d; in <source>ICOLD Symposium - dams under changing challenges</source>. Editors <person-group person-group-type="editor">
<name>
<surname>Schleiss</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
</person-group> (<publisher-loc>London, UK</publisher-loc>: <publisher-name>Taylor and Francis</publisher-name>), <fpage>403</fpage>&#x2013;<lpage>412</lpage>.</citation>
</ref>
<ref id="B5">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Auel</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Kantoush</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Sumi</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2016</year>). &#x201c;<article-title>Positive effects of reservoir sedimentation management on reservoir life - Examples from Japan</article-title>,&#x201d; in <conf-name>Proc. 84th ICOLD Annual meeting</conf-name>, <conf-loc>Johannesburg, South Africa</conf-loc>, <conf-date>May, 4.11&#x2013;14</conf-date>, <fpage>20</fpage>.</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Auel</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Kobayashi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Takemon</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Sumi</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Effects of sediment bypass tunnels on grain size distribution and benthic habitats in regulated rivers</article-title>. <source>Int. J. River Basin Manag.</source> <volume>15</volume>, <fpage>433</fpage>&#x2013;<lpage>444</lpage>. <pub-id pub-id-type="doi">10.1080/15715124.2017.1360320</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Baumer</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Radogna</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2015</year>). &#x201c;<article-title>Rehabilitation of the Palagnedra sediment bypass tunnel (2011-2013)</article-title>,&#x201d; in <conf-name>First International Workshop on Sediment Bypass Tunnels</conf-name>, <conf-loc>ETH Zurich, Switzerland</conf-loc>. Editor <person-group person-group-type="editor">
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
</person-group>, <fpage>235</fpage>&#x2013;<lpage>245</lpage>. <comment>Laboratory of Hydraulics, Hydrology and Glaciology (VAW)</comment>.</citation>
</ref>
<ref id="B8">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2015</year>). &#x201c;<article-title>First international workshop on sediment bypass tunnels</article-title>,&#x201d; in <source>VAW-Mitteilung 232</source>. Editor <person-group person-group-type="editor">
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
</person-group> <comment>Laboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich, Switzerland</comment>.</citation>
</ref>
<ref id="B9">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Auel</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Hagmann</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Albayrak</surname>
<given-names>I.</given-names>
</name>
</person-group> (<year>2014</year>). &#x201c;<article-title>Sediment bypass tunnels to mitigate reservoir sedimentation and restore sediment continuity</article-title>,&#x201d; in <source>Reservoir sedimentation</source>. Editors <person-group person-group-type="editor">
<name>
<surname>Schleiss</surname>
<given-names>A. J.</given-names>
</name>
<name>
<surname>de Cesare</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Franca</surname>
<given-names>M. J.</given-names>
</name>
<name>
<surname>Pfister</surname>
<given-names>M.</given-names>
</name>
</person-group> (<publisher-loc>London, UK</publisher-loc>: <publisher-name>Taylor and Francis Group</publisher-name>), <fpage>221</fpage>&#x2013;<lpage>228</lpage>.</citation>
</ref>
<ref id="B10">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Beck</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Lutz</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Lais</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Albayrak</surname>
<given-names>I.</given-names>
</name>
</person-group> (<year>2017</year>). &#x201c;<article-title>Hydraulics of water, air-water and sediment flow in downstream-controlled sediment bypass tunnels</article-title>,&#x201d; in <conf-name>Proc. 2nd Intl. Workshop on Sediment Bypass Tunnels</conf-name>, <conf-loc>Kyoto, Japan</conf-loc>, <fpage>1</fpage>&#x2013;<lpage>19</lpage>.</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bombar</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>El&#xe7;i</surname>
<given-names>&#x15e;.</given-names>
</name>
<name>
<surname>Tayfur</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>G&#xfc;ney</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Bor</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Experimental and numerical investigation of bed-load transport under unsteady flows</article-title>. <source>J. Hydraulic Eng.</source> <volume>137</volume>, <fpage>1276</fpage>&#x2013;<lpage>1282</lpage>. <pub-id pub-id-type="doi">10.1061/(asce)hy.1943-7900.0000412</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brandt</surname>
<given-names>S. A.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Classification of geomorphological effects downstream of dams</article-title>. <source>Catena</source> <volume>40</volume>, <fpage>375</fpage>&#x2013;<lpage>401</lpage>. <pub-id pub-id-type="doi">10.1016/S0341-8162(00)00093-X</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bruno</surname>
<given-names>M. C.</given-names>
</name>
<name>
<surname>Siviglia</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Carolli</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Maiolini</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Multiple drift responses of benthic invertebrates to interacting hydropeaking and thermopeaking waves</article-title>. <source>Ecohydrology</source> <volume>6</volume>, <fpage>511</fpage>&#x2013;<lpage>522</lpage>. <pub-id pub-id-type="doi">10.1002/eco.1275</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bunte</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>State of the science review, gravel mitigation and augmentation below hydroelectric dams: a geomorphological perspective</article-title>. <comment>Tech. rep</comment>.</citation>
</ref>
<ref id="B15">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Chang</surname>
<given-names>T.-J.</given-names>
</name>
<name>
<surname>Lai</surname>
<given-names>J.-S.</given-names>
</name>
</person-group> (<year>2019</year>). <conf-name>Proceedings of the third International Workshop on Sediment Bypass Tunnels</conf-name>, <conf-loc>Taipei, Taiwan</conf-loc>. <publisher-name>National Taiwan University</publisher-name>.</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cui</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Parker</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Pizzuto</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lisle</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Sediment pulses in mountain rivers: 2. Comparison between experiments and numerical predictions</article-title>. <source>Water Resour. Res.</source> <volume>39</volume>. <pub-id pub-id-type="doi">10.1029/2002wr001805</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Czapiga</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Blom</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Viparelli</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Sediment nourishments to mitigate channel bed incision in engineered rivers</article-title>. <source>J. Hydraulic Eng.</source> <volume>148</volume>. <pub-id pub-id-type="doi">10.1061/(asce)hy.1943-7900.0001977</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>[Dataset] Vetsch</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Siviglia</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ehrbar</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Facchini</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kammerer</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Koch</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2017a</year>). <source>Basement &#x2014; basic simulation environment for computation of environmental flow and natural hazard simulation</source>.</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Die Moran</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>El Kadi Qbderrezzak</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Mosselman</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Habersack</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Lebert</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Aelbrecht</surname>
<given-names>D.</given-names>
</name>
<etal/>
</person-group> (<year>2013</year>). <article-title>Physical model experiments for sediment supply to the old Rhine through induced bank erosion</article-title>. <source>Int. J. Sediment Res.</source> <volume>28</volume>, <fpage>431</fpage>&#x2013;<lpage>447</lpage>. <pub-id pub-id-type="doi">10.1016/S1001-6279(14)60003-2</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dietrich</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Kirchner</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ikeda</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Iseya</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>Sediment supply and the development of the coarse surface layer in gravel-bedded rivers</article-title>. <source>Nature</source> <volume>342</volume>, <fpage>215</fpage>&#x2013;<lpage>217</lpage>. <pub-id pub-id-type="doi">10.1038/340215a0</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elosegi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>D&#xed;ez</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Mutz</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Effects of hydromorphological integrity on biodiversity and functioning of river ecosystems</article-title>. <source>Hydrobiologia</source> <volume>657</volume>, <fpage>199</fpage>&#x2013;<lpage>215</lpage>. <pub-id pub-id-type="doi">10.1007/s10750-009-0083-4</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Facchini</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2018</year>). &#x201c;<article-title>Downstream moprhological effects of sediment bypass tunnels (ph.D. Thesis)</article-title>,&#x201d; in <source>VAW-Mitteilung 243</source>. Editor <person-group person-group-type="editor">
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
</person-group> <comment>Laboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich</comment>.</citation>
</ref>
<ref id="B23">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Facchini</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Siviglia</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2015</year>). &#x201c;<article-title>Downstream morphological impact of a sediment bypass tunnel: preliminary results and forthcoming actions</article-title>,&#x201d; in <conf-name>First International Workshop on Sediment Bypass Tunnels</conf-name>, <conf-loc>ETH Zurich, Switzerland</conf-loc>. Editor <person-group person-group-type="editor">
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
</person-group>, <fpage>137</fpage>&#x2013;<lpage>146</lpage>. <comment>Laboratory of Hydraulics, Hydrology and Glaciology (VAW)</comment>.</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ferrer-Boix</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Hassan</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Channel adjustments to a succession of water pulses in gravel bed rivers</article-title>. <source>Water Resour. Res.</source> <volume>51</volume>, <fpage>8773</fpage>&#x2013;<lpage>8790</lpage>. <pub-id pub-id-type="doi">10.1002/2015wr017664</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Fukuda</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Yamashita</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Osada</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Fukuoka</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2012</year>). &#x201c;<article-title>Study on flushing mechanism of dam reservoir sedimentation and recovery of riffle-pool in downstream reach by a flushing bypass tunnel</article-title>,&#x201d; in <conf-name>International Symposium on Dams for a Changing World</conf-name>, <fpage>1</fpage>&#x2013;<lpage>6</lpage>.</citation>
</ref>
<ref id="B26">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Hager</surname>
<given-names>W. H.</given-names>
</name>
<name>
<surname>Schleiss</surname>
<given-names>A. J.</given-names>
</name>
<name>
<surname>Boes</surname>
<given-names>R. M.</given-names>
</name>
<name>
<surname>Pfister</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2020</year>). <source>Hydraulic engineering of dams</source>. <publisher-loc>London</publisher-loc>: <publisher-name>Taylor and Francis</publisher-name>. <pub-id pub-id-type="doi">10.1201/9780203771433</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hassan</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Egozi</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Parker</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Experiments on the effect of hydrograph characteristics on vertical grain sorting in gravel bed rivers</article-title>. <source>Water Resour. Res.</source> <volume>42</volume>, <fpage>1</fpage>&#x2013;<lpage>15</lpage>. <pub-id pub-id-type="doi">10.1029/2005WR004707</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hauer</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Review of hydro-morphological management criteria on a river basin scale for preservation and restoration of freshwater pearl mussel habitats</article-title>. <source>Limnologica</source> <volume>50</volume>, <fpage>40</fpage>&#x2013;<lpage>53</lpage>. <pub-id pub-id-type="doi">10.1016/j.limno.2014.11.002</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hirano</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>1971</year>). <article-title>River bed degradation with armoring</article-title>. <comment>Transactions of the Japan Society of Civil Engineers</comment> <volume>3</volume>, <fpage>194</fpage>&#x2013;<lpage>195</lpage>.</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hirano</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>1972</year>). <article-title>Studies on variation and equilibrium state of a river bed composed of nonuniform material</article-title>. <comment>Transactions of the Japan Society of Civil Engineers</comment> <volume>4</volume>, <fpage>128</fpage>&#x2013;<lpage>129</lpage>.</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Humphries</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Venditti</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Sklar</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wooster</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Experimental evidence for the effect of hydrographs on sediment pulse dynamics in gravel-bedded rivers</article-title>. <source>Water Resour. Res.</source> <volume>48</volume>. <pub-id pub-id-type="doi">10.1029/2011wr010419</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Khosravi</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Chegini</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Binns</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Daggupati</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Mao</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Difference in the bed load transport of graded and uniform sediments during floods: an experimental investigation</article-title>. <source>Hydrology Res.</source> <volume>50</volume>, <fpage>1645</fpage>&#x2013;<lpage>1664</lpage>. <pub-id pub-id-type="doi">10.2166/nh.2019.078</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kobayashi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Fukuroi</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Sumi</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Takemon</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2018</year>). &#x201c;<article-title>Sediment derivation by bypass tunnel restores downstream environment</article-title>,&#x201d; in <source>Twenty-sixth international congress on large dams</source> (<publisher-loc>Vienna, Austria</publisher-loc>, <publisher-name>CRC Press</publisher-name>), <fpage>3464</fpage>&#x2013;<lpage>3465</lpage>. <pub-id pub-id-type="doi">10.1201/9780429465086</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kondolf</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>PROFILE: hungry water: effects of dams and gravel mining on river channels</article-title>. <source>Environ. Manag.</source> <volume>21</volume>, <fpage>533</fpage>&#x2013;<lpage>551</lpage>. <pub-id pub-id-type="doi">10.1007/s002679900048</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kondolf</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Annandale</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Morris</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2014</year>). <article-title>Sustainable sediment management in reservoirs and regulated rivers: experiences from five continents</article-title>. <source>Earth&#x2019; s Future</source> <volume>2</volume>, <fpage>256</fpage>&#x2013;<lpage>280</lpage>. <pub-id pub-id-type="doi">10.1002/2013ef000184</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kondolf</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Matthews</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>1991</year>). &#x201c;<article-title>Management of coarse sediment in regulated rivers of California</article-title>,&#x201d; in <source>UC berkley technical completion reports</source> (<publisher-loc>Davis, California</publisher-loc>: <publisher-name>University of California Water Resources Center</publisher-name>), <volume>80</volume>, <fpage>102</fpage>.</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kondolf</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Wilcock</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>1996</year>). <article-title>The flushing flow problem: defining and evaluating objectives</article-title>. <source>Water Resour. Res.</source> <volume>32</volume>, <fpage>2589</fpage>&#x2013;<lpage>2599</lpage>. <pub-id pub-id-type="doi">10.1029/96wr00898</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Kondolf</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Williams</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Horner</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Milan</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2008</year>). &#x201c;<article-title>Assessing physical quality of spawning habitat</article-title>,&#x201d; in <conf-name>American Fisheries Society Symposium</conf-name>. Editors <person-group person-group-type="editor">
<name>
<surname>Sear</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>DeVries</surname>
<given-names>P.</given-names>
</name>
</person-group>, <fpage>1</fpage>&#x2013;<lpage>26</lpage>.</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lane</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>1955</year>). <article-title>The importance of fluvial morphology in hydraulic engineering. Hydraulic Laboratory Report No. 372, Denver, Colorado, United States: Engineering Laboratories of the US Department of the Interior, Bureau of Reclamation</article-title>. <volume>81</volume>, <fpage>1</fpage>&#x2013;<lpage>17</lpage>.</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lisle</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>Sediment transport and resulting deposition in spawning gravels, north coastal California</article-title>. <source>Water Resour. Res.</source> <volume>25</volume>, <fpage>1303</fpage>&#x2013;<lpage>1319</lpage>. <pub-id pub-id-type="doi">10.1029/wr025i006p01303</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mackin</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>1948</year>). <article-title>Concept of the graded river</article-title>. <source>Geol. Soc. Am. Bull.</source> <volume>59</volume>, <fpage>463</fpage>&#x2013;<lpage>512</lpage>. <pub-id pub-id-type="doi">10.1130/0016-7606(1948)59[463:cotgr]2.0.co;2</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mao</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>The effect of hydrographs on bed load transport and bed sediment spatial arrangement</article-title>. <source>J. Geophys. Res. Earth Surf.</source> <volume>117</volume>, <fpage>1</fpage>&#x2013;<lpage>16</lpage>. <pub-id pub-id-type="doi">10.1029/2012JF002428</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mart&#xed;n</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Doering</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Robinson</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Ecological assessment of a sediment by-pass tunnel on a receiving stream in Switzerland</article-title>. <source>River Res. Appl.</source> <volume>33</volume>, <fpage>925</fpage>&#x2013;<lpage>936</lpage>. <pub-id pub-id-type="doi">10.1002/rra.3145</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Merz</surname>
<given-names>J. E.</given-names>
</name>
<name>
<surname>Ochikubo Chan</surname>
<given-names>L. K.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Effects of gravel augmentation on macroinvertebrate assemblages in a regulated California River</article-title>. <source>River Res. Appl.</source> <volume>21</volume>, <fpage>61</fpage>&#x2013;<lpage>74</lpage>. <pub-id pub-id-type="doi">10.1002/rra.819</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Mueller-Hagmann</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2017</year>). &#x201c;<article-title>Hydroabrasion in high speed flows at sediment bypass tunnels (Ph.D. thesis)</article-title>,&#x201d; in <source>VAW-Mitteilung 239</source>. Editor <person-group person-group-type="editor">
<name>
<surname>Boes</surname>
<given-names>R.</given-names>
</name>
</person-group> <comment>Laboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich</comment>.</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nagrodski</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Raby</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Hasler</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Taylor</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Cooke</surname>
<given-names>S. J.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Fish stranding in freshwater systems: sources, consequences, and mitigation</article-title>. <source>J. Environ. Manag.</source> <volume>103</volume>, <fpage>133</fpage>&#x2013;<lpage>141</lpage>. <pub-id pub-id-type="doi">10.1016/j.jenvman.2012.03.007</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nelson</surname>
<given-names>P. A.</given-names>
</name>
<name>
<surname>Venditti</surname>
<given-names>J. G.</given-names>
</name>
<name>
<surname>Dietrich</surname>
<given-names>W. E.</given-names>
</name>
<name>
<surname>Kirchner</surname>
<given-names>J. W.</given-names>
</name>
<name>
<surname>Ikeda</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Iseya</surname>
<given-names>F.</given-names>
</name>
<etal/>
</person-group> (<year>2009</year>). <article-title>Response of bed surface patchiness to reductions in sediment supply</article-title>. <source>J. Geophys. Res. Earth Surf.</source> <volume>114</volume>. <pub-id pub-id-type="doi">10.1029/2008jf001144</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Parker</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2004</year>). <source>1D Sediment transport morphodynamics with applications to rivers and turbidity currents</source> <publisher-loc>Champaign, Illinois, USA</publisher-loc>: <publisher-name>University of Illinois - Ven Te Chow Hydrosystems Laboratory</publisher-name>.</citation>
</ref>
<ref id="B49">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Parker</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Transport of gravel and sediment mixtures</article-title>. In <source>Sedimentation engineering: processes, measurements, modeling, and practice</source>, ed. <person-group person-group-type="editor">
<name>
<surname>Garcia</surname>
<given-names>M. H.</given-names>
</name>
</person-group> <publisher-name>ASCE American Society of Civil Engineers</publisher-name>, <fpage>165</fpage>&#x2013;<lpage>251</lpage>.</citation>
</ref>
<ref id="B50">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Parker</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Hassan</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wilcock</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2008</year>). &#x201c;<article-title>Adjustment of the bed surface size distribution of gravel-bed rivers in response to cycled hydrographs</article-title>,&#x201d; in <source>Gravel-bed rivers VI: from process understanding to river restoration</source>, <fpage>1</fpage>&#x2013;<lpage>47</lpage>.</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Phillips</surname>
<given-names>C. B.</given-names>
</name>
<name>
<surname>Hill</surname>
<given-names>K. M.</given-names>
</name>
<name>
<surname>Paola</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Singer</surname>
<given-names>M. B.</given-names>
</name>
<name>
<surname>Jerolmack</surname>
<given-names>D. J.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Effect of flood hydrograph duration, magnitude, and shape on bed load transport dynamics</article-title>. <source>Geophys. Res. Lett.</source> <volume>45</volume>, <fpage>8264</fpage>&#x2013;<lpage>8271</lpage>. <pub-id pub-id-type="doi">10.1029/2018GL078976</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pitlick</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wilcock</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Relations between streamflow, sediment transport, and aquatic habitat in regulated rivers</article-title>. <source>Water Sci. Appl.</source> <volume>4</volume>, <fpage>185</fpage>&#x2013;<lpage>198</lpage>. <pub-id pub-id-type="doi">10.1029/WS004p0185</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Reid</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Frostick</surname>
<given-names>L. E.</given-names>
</name>
<name>
<surname>Layman</surname>
<given-names>J. T.</given-names>
</name>
</person-group> (<year>1985</year>). <article-title>The incidence and nature of bedload transport during flood flows in coarse-grained alluvial channels</article-title>. <source>Earth Surf. Process. Landforms</source> <volume>10</volume>, <fpage>33</fpage>&#x2013;<lpage>44</lpage>. <pub-id pub-id-type="doi">10.1002/esp.3290100107</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Ribberink</surname>
<given-names>J. S.</given-names>
</name>
</person-group> (<year>1987</year>). <source>Mathematical modelling of one-dimensional morphological changes in rivers with non-uniform sediment</source>. <comment>Ph.D. thesis</comment>. <comment>TU Delft, The Netherlands</comment>.</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Serrana</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yaegashi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Kondoh</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Robinson</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Watanabe</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Ecological influence of sediment bypass tunnels on macroinvertebrates in dam-fragmented rivers by DNA metabarcoding</article-title>. <source>Sci. Rep.</source> <volume>8</volume>, <fpage>10185</fpage>&#x2013;<lpage>10210</lpage>. <pub-id pub-id-type="doi">10.1038/s41598-018-28624-2</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Siviglia</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Stecca</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Blom</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2017</year>). &#x201c;<article-title>Modelling of mixed&#x2013;sediment morphodynamics in gravel&#x2013;bed rivers using the active&#x2013;layer approach</article-title>,&#x201d; in <source>Gravel-bed rivers: processes and disasters</source>. Editors <person-group person-group-type="editor">
<name>
<surname>Tsutsumi</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Laronne</surname>
<given-names>J.</given-names>
</name>
</person-group> (<publisher-loc>Chichester, UK</publisher-loc>: <publisher-name>John Wiley and Sons, Ltd). John wiley edn.</publisher-name>), <fpage>703</fpage>&#x2013;<lpage>728</lpage>. <pub-id pub-id-type="doi">10.1002/9781118971437.ch26</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Smart</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Jaeggi</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>1983</year>). &#x201c;<article-title>Sediment transport on steep slopes</article-title>,&#x201d;. <source>VAW-Mitteilung</source>. Editor <person-group person-group-type="editor">
<name>
<surname>Vischer</surname>
<given-names>D.</given-names>
</name>
</person-group> (<publisher-loc>Zurich, Switzerland</publisher-loc>: <publisher-name>Laboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich</publisher-name>), <volume>64</volume>, <fpage>1</fpage>&#x2013;<lpage>191</lpage>.</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stecca</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Siviglia</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Blom</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Mathematical analysis of the Saint-Venant-Hirano model for mixed-sediment morphodynamics</article-title>. <source>Water Resour. Res.</source> <volume>50</volume>, <fpage>7563</fpage>&#x2013;<lpage>7589</lpage>. <pub-id pub-id-type="doi">10.1002/2014WR015251</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stecca</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Siviglia</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Blom</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>An accurate numerical solution to the Saint-Venant-Hirano model for mixed-sediment morphodynamics in rivers</article-title>. <source>Adv. Water Resour.</source> <volume>93 Part A</volume>, <fpage>39</fpage>&#x2013;<lpage>61</lpage>. <pub-id pub-id-type="doi">10.1016/j.advwatres.2015.05.022</pub-id>
</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Strickler</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>1923</year>). <article-title>Beitr&#xe4;ge zur Frage der Geschwindigkeitsformel und der Rauhigkeitszahlen f&#xfc;r Str&#xf6;me, Kan&#xe4;le und geschlossene Leitungen</article-title>. <comment>No. 16 in Amt f&#xfc;r Wasserwirtschaft: Mitteilungen des Amtes f&#xfc;r Wasserwirtschaft (Sekretariat des Eidgen&#xf6;ssischen Amtes f&#xfc;r Wasserwirtschaft)</comment>
</citation>
</ref>
<ref id="B61">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Sumi</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2017</year>). <source>Proceedings of the second international workshop on sediment Bypass tunnels</source>. <publisher-loc>Japan</publisher-loc>: <publisher-name>DPRI, University of Kyoto</publisher-name>.</citation>
</ref>
<ref id="B62">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Sumi</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Kantoush</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Suzuki</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2012</year>). &#x201c;<article-title>Performance of Miwa Dam sediment bypass tunnel: evaluation of upstream and downstream state and bypassing efficiency</article-title>,&#x201d; in <conf-name>International Symposium on Dams for a Changing World (Kyoto: International COmmission on Large Dams (ICOLD))</conf-name>, <fpage>576</fpage>&#x2013;<lpage>596</lpage>.</citation>
</ref>
<ref id="B63">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Surian</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Rinaldi</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Morphological response to river engineering and management in alluvial channels in Italy</article-title>. <source>Geomorphology</source> <volume>50</volume>, <fpage>307</fpage>&#x2013;<lpage>326</lpage>. <pub-id pub-id-type="doi">10.1016/S0169-555X(02)00219-2</pub-id>
</citation>
</ref>
<ref id="B64">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Toro-Escobar</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Paola</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Parker</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>1996</year>). <article-title>Transfer function for the deposition of poorly sorted gravel in response to streambed aggradation</article-title>. <source>J. Hydraulic Res.</source> <volume>34</volume>, <fpage>35</fpage>&#x2013;<lpage>53</lpage>. <pub-id pub-id-type="doi">10.1080/00221689609498763</pub-id>
</citation>
</ref>
<ref id="B65">
<citation citation-type="book">
<collab>VAW</collab> (<year>2010</year>). <source>Solis sediment Bypass tunnel &#x2014; hydraulic study (geschiebeumleitstollen Solis &#x2014; hydraulische modellversuche)</source>. <comment>Tech. Rep. 4269, Laboratory of Hydraulics, Hydrology and Glaciology (VAW), [in German], ETH Zurich, Zurich, Switzerland</comment>.</citation>
</ref>
<ref id="B66">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Vetsch</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Siviglia</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ehrbar</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Facchini</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kammerer</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Koch</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2017b</year>). <article-title>System manuals of BASEMENT</article-title>. <source>Laboratory of hydraulics, glaciology and hydrology (VAW)</source>. <publisher-loc>Zurich, Switzerland</publisher-loc>: <publisher-name>ETH Zurich</publisher-name>. <comment>Version 2.7</comment>.</citation>
</ref>
<ref id="B67">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Viparelli</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Gaeuman</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Wilcock</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Parker</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>A model to predict the evolution of a gravel bed river under an imposed cyclic hydrograph and its application to the Trinity River</article-title>. <source>Water Resour. Res.</source> <volume>47</volume>. <pub-id pub-id-type="doi">10.1029/2010wr009164</pub-id>
</citation>
</ref>
<ref id="B68">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Vischer</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Hager</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Casanova</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Joos</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Lier</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Martini</surname>
<given-names>O.</given-names>
</name>
</person-group> (<year>1997</year>). &#x201c;<article-title>Bypass tunnels to prevent reservoir sedimentation</article-title>,&#x201d; in <source>
<italic>International Congress on large dams</italic> (florence: international COmmission on large dams (ICOLD))</source>, <fpage>605</fpage>&#x2013;<lpage>624</lpage>.</citation>
</ref>
<ref id="B69">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Waters</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Curran</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Linking bed morphology changes of two sediment mixtures to sediment transport predictions in unsteady flows</article-title>. <source>Water Resour. Res.</source> <volume>51</volume>, <fpage>2724</fpage>&#x2013;<lpage>2741</lpage>. <pub-id pub-id-type="doi">10.1002/2014wr016083</pub-id>
</citation>
</ref>
<ref id="B70">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wilcock</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Crowe</surname>
<given-names>J. C.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Surface-based transport model for mixed-size sediment</article-title>. <source>J. Hydraulic Eng.</source> <volume>129</volume>, <fpage>120</fpage>&#x2013;<lpage>128</lpage>. <pub-id pub-id-type="doi">10.1061/(asce)0733-9429(2003)129:2(120)</pub-id>
</citation>
</ref>
<ref id="B71">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wilcock</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>MacArdell</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>1993</year>). <article-title>Surface&#x2010;based fractional transport rates: mobilization thresholds and partial transport of a sand&#x2010;gravel sediment</article-title>. <source>Water Resour. Res.</source> <volume>29</volume>, <fpage>1297</fpage>&#x2013;<lpage>1312</lpage>. <pub-id pub-id-type="doi">10.1029/92wr02748</pub-id>
</citation>
</ref>
<ref id="B72">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Williams</surname>
<given-names>G. P.</given-names>
</name>
<name>
<surname>Wolman</surname>
<given-names>M. G.</given-names>
</name>
</person-group> (<year>1984</year>). <source>Downstream effects of dams on alluvial rivers</source>, <volume>1286</volume>. <publisher-loc>Washington, DC, USA</publisher-loc>: <publisher-name>U.S. Geological Surveys Professional Papers (U.S. Geological Surveys</publisher-name>), <fpage>83</fpage>. <pub-id pub-id-type="doi">10.1126/science.277.5322.9j</pub-id>
</citation>
</ref>
<ref id="B73">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Wohl</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2013</year>). <source>Wide rivers crossed: the south platte and the Illinois of the American prairie</source>. <publisher-loc>Boulder, CO, USA</publisher-loc>: <publisher-name>University Press of Colorado</publisher-name>.</citation>
</ref>
<ref id="B74">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wohl</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Lane</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wilcox</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>The science and practice of river restoration</article-title>. <source>Water Resour. Res.</source> <volume>51</volume>, <fpage>5974</fpage>&#x2013;<lpage>5997</lpage>. <pub-id pub-id-type="doi">10.1002/2014WR016874</pub-id>
</citation>
</ref>
<ref id="B75">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wong</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Parker</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>One-dimensional modeling of bed evolution in a gravel bed river subject to a cycled flood hydrograph</article-title>. <source>J. Geophys. Res. Earth Surf.</source> <volume>111</volume>. <pub-id pub-id-type="doi">10.1029/2006JF000478</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>