Effects of Transverse Groynes on Meso-Habitat Suitability for Native Fish Species on a Regulated By-Passed Large River: A Case Study along the Rhine River

: River regulations ultimately degrade ﬂuvial forms and morphodynamics and simplify riparian and aquatic habitats. For several decades, river restoration actions have been performed to recover geomorphic processes and diversify these habitats to enhance both river biodiversity and ecosystem services. The objective of this study is to provide quantitative feedback on the experimental restoration of a large regulated and by-passed river (the Upper Rhine downstream of the Kembs Dam, France/Germany). This restoration consisted of the construction of two transverse groynes and the removal of bank protection. A monitoring framework composed of topo-bathymetric surveys as well as ﬂow velocity and grain size measurements was established to assess the channel morphodynamic responses and evaluate their effects on habitat suitability for ﬁve native ﬁsh species using habitat models. A riverscape approach was used to evaluate the landscape changes in terms of both the conﬁguration and the composition, which cannot be considered with classic approaches (e.g., Weighted Usable Area). Our results show that the two transverse groynes and, to a lesser extent, bank erosion, which was locally enhanced by the two groynes, increased habitat diversity due to the creation of new macroforms (e.g., pools and mid-bars) and ﬁning of the bed grain size. Using a riverscape approach, our ﬁndings highlight that the restoration improved eel and juvenile nase species due to slowing down of the current and the deposition of ﬁne sediments downstream of both groynes. As a consequence, the restoration improved the habitat suitability of the studied reach for more ﬁsh species compared with the pre-restoration conditions. This study also demonstrates that the salmon habitats downstream of the restored reach were improved due to ﬁning of the bed grain size. This ﬁnding highlights that, for restorations aimed at ﬁsh habitats, the grain size conditions must be taken into consideration along with the ﬂow conditions. Furthermore, the implementation of groynes, while not a panacea in terms of functional restoration, can be a strategy for improving ﬁsh habitats on highly regulated rivers, but only when more functional and natural options are impossible due to major constraints.


Introduction
Most large rivers have been regulated for the purposes of flood protection, navigation, and agricultural development, as well as electricity production. However, these regulations have induced the degradation of both geomorphological and hydrological conditions (e.g., the fragmentation of river networks and a deficit of sediment) [1,2]. As a result, river habitats have been heavily disturbed, altering aquatic and riverine communities and biodiversity [3,4].
Over the past several decades, an increasing number of river restorations has been undertaken, notably in Europe, where this rise in restoration programs has been fostered by the European Water Framework Directive (WFD) [5]. The main objective of such restoration actions is to improve the functional processes of rivers and thus recover ecosystem services [6][7][8][9][10]. River restoration actions can be implemented for one or more of the seven following categories: continuity restoration, pollution control, recovery of sediment dynamics, land-use management, adaptation to climate change, implementation of ecological discharges, and habitat restoration [11].
Numerous restorations targeted the recovering habitats of some native fish species due to their emblematic value and habitat requirements. For example, since the 1960s, restorations of rivers in California, USA were performed to recover spawning areas for salmon by using gravel augmentations, which were sometimes combined with artificial riffles [12]. Either fish passes were launched or low dams were removed to restore river continuity [13]. More recently, several engineering strategies were tested along the Danube River, such as (i) riprap removal and groyne shortening, (ii) riprap and groyne field removal, and (iii) side channel reconnection, to increase the abundance of rheophilic fishes [14]. Historically, groynes were structures composed by wood and/or rocks and used to promote bed erosion in the main channel and maintain a minimum water level for navigation purposes [15]. Recently, longitudinal training dams (LTDs) were tested along the Rhine river to simultaneously improve navigation during low flows, enhance safety for human activities and hydraulic infrastructures against floods, and create refugia for native fish at their juvenile stages [16]. On the channelized Godwin Creek River, stones were added to extend 18 small groynes to promote channel diversification and enhance fish habitats [17]; this research highlighted that the extension of groynes promoted habitat heterogeneity by increasing pool habitats and promoting grain size fining, following which the composition and relative abundance of fish were similar to those of the non-incised reference site. Lacey and Millar [18] in their research used hydrodynamic models on the Chiliwack River to assess the effect of large wood within the stream and the implementation of rock groyne habitat structures (both simple and double rock groynes); their modeling results showed that these actions promoted an increase in habitat suitability for two salmonids (Oncorhynchus kisutch, Oncorhynchus mykiss) during high flows. Fish species can be viewed as umbrella species, i.e., species whose conservation is expected to confer protection to a large number of naturally co-occurring species [19]. Hence, a restoration implemented for a native fish community is usually good for the whole ecosystem. Fish react promptly to human modifications of an ecosystem and thus can be used as indicators of river recovery and, more broadly, the ecological status of waterbodies [20,21]. Thorough knowledge of the biology of fish species and the requirements for those species to complete their life cycle makes it possible to predict which habitats will be able to host a given species or a given functional category of fish. Consequently, fish habitat models are frequently employed to determine minimal in-stream discharges to improve the habitat suitability for native fish species on by-passed reaches [22][23][24], to assess habitat changes induced by hydropeaking or weir and reservoir removal [25], and to evaluate the effects of restoration actions [26,27]. The most commonly used modeling method is the in-stream flow incremental methodology (IFIM) developed by the U.S. Fish and Wildlife Service during the 1970s [28]. This methodology evaluates habitat suitability through the habitat suitability index (HSI), derived from fish species preferences and related to numerous environmental variables (e.g., water depth, flow velocity, and grain size of the bed substrate) [20,[29][30][31].
Most fish habitat models are based on quantitative estimations of habitat suitability using microhabitat approaches. However, the composition and configuration of riverscape units and the associated meso-habitats at a larger scale are rarely taken into account, although ecological processes at the population and community levels are affected by these factors [32]. The use of a meso-habitat rather than a micro-habitat approach allows these changes to be quantified [33] and is potentially useful for detecting riverscape changes following natural or human disturbances [34]. Based on this approach, a new ecological hydrograph relating flow to habitat quality for Spinibarbus hollandi on the Lijiang River was proposed by Li et al. [33]. The advantage of using the riverscape approach rather than the micro-scale approach to assess the effects of low dams on fish habitat quality was shown on the Neosho River by Hitchman et al. [34].
A riverscape approach was used to quantitatively assess the effects of channel restoration based on the implementation of two transverse groynes on habitat diversification and a rise in potential fish abundance using riverscape metrics. Two hypotheses were tested: (i) the restoration action promoted the diversification of aquatic habitats, and (ii) the diversification of aquatic habitats improved the habitat suitability for the native fish studied. Geomorphological monitoring based on a before-aftercontrol-impact (BACI) design [35] was achieved to evaluate the channel responses following the restoration, allowing us to calibrate two-dimensional hydraulic models to build fish habitat models for five native species based on empirical preferences determined along French rivers [36,37].

Study Area
The Rhine River is the third largest river in Europe with a length and drainage basin equal to 1250 km and 185,000 km 2 , respectively. The Upper Rhine, a 300 km long section located between Basel and Bingen-am-Rhein, shows a nivo-glacial hydrological regime (Figure 1a). The mean annual discharge has been estimated at 1059 m 3 /s at the Basel gauging station (during the period 1891-2008) [38]. Since the beginning of the 19th century, the Upper Rhine has been profoundly impacted by several engineering works. These works induced the severe degradation of both geomorphological conditions (e.g., simplification of channel patterns, narrowing and stabilization of the main channel, bed incision, and armoring) and ecological function (losses of habitat diversity and dynamics). These effects were notably observed along a 50 km by-passed reach parallel to the Grand Canal d'Alsace (GCA) downstream of Basel [39] (Figure 1b). Since 2010, several restoration actions, including an increase in instream flow [40], gravel augmentations [41,42], and controlled bank erosion with the implementation of two transverse groynes [43], have been performed along this by-passed reach to enhance the bedload transport and morphodynamics therein and thus diversify aquatic and riverine habitats. In December 2010, the instream flow entering the Old Rhine increased from a quasi-constant discharge of 20-30 m 3 /s to a seasonally varying discharge of 52-150 m 3 /s in accordance with the life cycle of native fish species [44]. An initial experiment of controlled bank erosion was implemented in 2013 over a span of 300 m on the upstream part of the Old Rhine ; April-May 2013). This ecological action consisted of removing the bank protection and three old groynes built during the beginning of the 20th century ( Figure 1d) and implementing two new transverse groynes within the flow channel to maximize flow velocities and shear stresses along the bank. The location and shape of these groynes had been previously defined through numerical and physical modeling to promote bank erosion, deliver sediment into the channel, and avoid any damage to human construction due to uncontrolled lateral erosion [45].
The monitoring was composed of four states (S i ), namely, an initial state (S 0 ) and three periods (P i ). Since this restoration, the flood peaks reached 2600 m 3 /s, 1400 m 3 /s and 1590 m 3 /s during P 1 , P 2 and P 3 , respectively ( Figure 2). The flood discharges were measured at the Rheinweiler gauging station (located 5 km upstream of the restored site). No tributaries are located between this gauging station and the studied reach. Following these flood events, only few bank erosion events occurred. However, significant and unexpected geomorphological changes, of the bathymetry and bed grain size, were observed that show that transverse groynes had significant effects on aquatic habitats. The latter motivated the implementation of our research strategy.   [46]. Stars correspond to topo-bathymetric surveys, squares correspond to the collection of physical parameters, and the grey arrow denotes the end of the restoration.

Field Monitoring
Field monitoring was performed following a before-after-control-impact (BACI) framework [35] on three river sections: the upstream, impacted, and downstream sections ( Figure 1). The upstream and downstream sections were used as control sections, and the impacted section was judged according to the impact of the restoration action.
In March 2013 (S 0 ), a bathymetric survey was performed using a Tritech PA500 single beam echo-sounder (Tritech International Limited, Aberdeenshire, Scotland) to survey the channel bathymetry. Since March 2014 (S 1 ), yearly airborne topo-bathymetric LiDAR surveys have been performed to systematically characterize the channel geometry with a high resolution in both above-water and underwater conditions along the first upstream 20 km of the by-passed channel. For S 1 , a Riegl VQ-820-G (RIEGL Laser Measurement Systems, Horn, Austria) sensor was used, whereas an Optech Titan sensor (Teledyne Optech, Vaughan Canada) was used during S 2 (March 2016) and S 3 (March 2017) (Figure 3). High-resolution digital elevation models (DEMs; cell size of 0.25 m 2 ) were produced using Arcmap software v.10.3 (ESRI, Redlands, California, USA). The channel elevation changes for each period (P i ) were calculated by the DEMs of difference (DoDs) between two surveys. Field measurements of water depths and flow velocities and visual estimates of the bed grain size were conducted along cross-sections (n = 18) at a spacing of 100 m with 12 to 19 sampling points for each cross-section, and the measurements were collected by a private company [47][48][49][50]. All sampling points were georeferenced by a global positioning system (Promark 700). The same sampling points were recorded during each survey to evaluate the changes in physical parameters as accurately as possible ( Figure 1c). A rod and a Nautilus C2000 current meter were used to measure the water depth and the flow velocity, respectively, at 0.4 times the water depth from the channel bottom using a boat. A bathyscope was used to qualitatively estimate the bed grain size according to the ordinal classification of [51]. The discharges during data collection for S 1 , S 2 , S 3 , and S 4 were approximately 67 m 3 /s, 52 m 3 /s, 70 m 3 /s, and 70 m 3 /s, respectively ( Figure 2).

Grain Size Quantification
The qualitative estimates of the bed grain size were transformed into quantitative values on the basis of the EVHA codification [52]: Grain size value = 0.2 * coarsest substrate + 0.4 * dominant substrate 1 + 0.4 * dominante substrate 2 (1) In the case of only one dominant substrate: The bed grain size along the entire study reach was interpolated following the triangulation interpolation method.

Hydraulic Models
A two-dimensional hydraulic model was performed for each state (S i ) by using HEC-RAS software (v.5.0.3). The topo-bathymetric data for the models were provided by high-resolution DEMs. Each model run was calibrated with the field measurements of the water depths and flow velocities performed at each S i and correlated with the daily discharges provided by the Rheinweiler gauging station (Figures 1b and 2). Sensitivity analyses of the Manning coefficient were performed to determine the best calibration values for modeling purposes. The values of the Manning coefficient varied over the four S i from 0.025 to 0.030, similar to the values proposed by Chow [53] and Benson et al. [54] for rivers similar to this by-passed reach. For each state, the water depths and flow velocities were modeled for the studied reach at 52 m 3 /s because this discharge has been the most frequent daily discharge recorded since December 2010.

Estimation of Aquatic Habitat Heterogeneity
Habitat heterogeneity was estimated using the eco-hydro-morphological index of diversity (EHMID) developed by [55] following the hydro-morphological index of diversity (HMID) proposed by [56]. The EHMID was slightly modified because the suitability curves following the grain size were based on the EVHA codification rather than the median grain size, as proposed by [55]. The metric was normalized by the section width, the resulting metric, nEHMID m , was calculated as follows: ) 2 correspond to the partial diversity of the flow velocities (m.s −1 ), water depths (m), and bed grain sizes (codification), respectively. The cross-section diversity (CSD) developed by Gostner et al. [56] was calculated for each crosssection to evaluate the relation between the physical parameters changes and the channel geometry changes. This metric was calculated with the following formula: where |∆Y i | is the absolute height difference in the channel bottom between two consecutive points on a cross-section (m) and X i is the distance between those two points (m).

Fish Habitat Models and Metrics
In this study, the suitability curves were based on statistical relationships defined for six large rivers in France according to field observations (Table 1). Suitability curves were established for three categories of physical parameters: water depth, flow velocity, and bed grain size ( Figure 4). We estimated the habitat quality for five native fish observed in the study reach (mostly for the adult and juvenile life stages): barbel (Barbus barbus), nase (Chondrostoma nasus), grayling (Thymallus thymallus), salmon (Salmo salar), and eel (Anguilla anguilla). Table 1. References for the suitability curves used in this study for each native fish. (a) indicates studies devoted to by-passed reaches.

Common and Scientific Names Life Stages Source Studied Reaches in the Fish Surveys
Eel (Anguilla anguilla) [36] Rhône river (a); Ain river; Ardèche river; Drôme river; Loire river; Garonne river Barbel (Barbus barbus) Adult and juvenile [36] Rhône river (a); Ain river; Ardèche river; Drôme river; Loire river; Garonne river Nase (Chondrostoma nasus) Adult and juvenile [36] Rhône river (a); Ain river; Ardèche river; Drôme river; Loire river; Garonne river; Grayling (Thymallus thymallus) Adult and juvenile [37] Ain river Salmon (Salmo salar) Adult and juvenile Expert curves revised in 1997 Habitat suitability was estimated according to the suitability index (SI) calculated at the cell scale (0.25 m 2 for each physical variable varying from 0 (bad habitat quality) to 1 (excellent habitat quality) ( Figure 3). According to [29,57,58], the combined suitability index (CSI), which corresponds to the product among the SI values of the n physical variables was calculated as follows: where i, j, and n correspond to the cell in which the computation is occurring, the physical parameter being multiplied, and the number of physical parameters, respectively.
From the CSI, the habitat suitability index (HSI) was calculated to estimate the habitat suitability at the section scale and calculated it as follows [29,59,60]: where a i is the surface area of cell i and A is the cross-sectional area of the section being studied (upstream, impacted, and downstream sections). All indices were calculated by using ArcGIS v10.3 software based on Python scripts and calculated from the results of the two-dimensional hydraulic model runs. In addition, the following three metrics were calculated using FRAGSTATS software v.4.0 to take into account both the configuration and the composition of the riverscape (Figure 4): i The highly suitable habitat proportion (HSHP): ii The largest patch index (LPI): iii The landscape division index (LDI): where a ij is equal to the surface area of patch ij and A is the cross-sectional area of the section being studied. These three metrics were calculated for the highly suitable habitat, which corresponds to a CSI equal to or greater than 0.5. All of the metrics were estimated for the upstream, impacted, and downstream sections and for each S i (Figures 1c and 3). Example of the landscape evolution for a given HSI value and the dynamics of the landscape metrics. HSI, LPI, LDI, and HSHP correspond to the habitat suitability index, the largest patch index of a highly suitable habitat, the landscape division index, and the highly suitable habitat proportion, respectively. LPI, LDI, and HSHP were calculated for patches that were judged to be highly suitable (HSI ≥ 0.5) and are represented in black.
Furthermore, a synthetic index varying from 0 (bad) to 1 (excellent) was developed to quantify habitat suitability; this index is called the synthetic habitat suitability index (SHSI), which combines the four previous metrics with identical weighting, and is calculated as follows: Then, we calculated the temporal change ratio (TCR) of the SHSI for each period as follows:

Channel Responses
During P 1 , contrasting channel adjustments occurred following the intense flood of May-June 2013 (Q ix < 15 years). The bed was scoured by 0.1 m to 0.5 m in the upstream section, whereas both bed scouring (max. of 1 m) and bed deposition (max. of 0.7 m) occurred on the two riffles located in the downstream section (Figure 5a). In the impacted section, new macroforms appeared with the formation of two mid-bars downstream of the two artificial groynes, where the deposition thickness ranged from 0.1 m to 2.5 m (Figure 5a-c), and scouring (max. of 4 m) occurred in two pools at the eastern extremity of the artificial groynes. During P 2 , the upstream and downstream sections morphologically evolved less than during P 1 , with less than 0.5 m of scouring and deposition occurring. Moreover, in the impacted section, the bed was globally eroded. During P 3 , generalized scour of the bed occurred along the whole study reach, reaching a maximum of 2 m locally in the upstream pool. Nevertheless, macroforms were clearly still present during the whole monitoring period (Figure 5d). These channel responses induced a significant increase in the CSD values within the impacted section during the whole monitoring period, whereas the CSD values in the upstream and downstream sections remained globally stable (Figure 6e).

Precision of the Hydraulic Models
The hydraulic models reliably simulated the two hydraulic variables, namely, the water depth and flow velocity. The mean absolute error (MAE) values for the simulated water depths are equal to 0.18 m, 0.05 m, 0.08 m, and 0.11 m for S 0 , S 1 , S 2 , and S 3 , respectively (Figure 6a-d). Likewise, the MAE values for the simulated flow velocities are equal to 0.17 m/s, 0.14 m/s, 0.19 m/s, and 0.15 m/s for S 0 , S 1 , S 2 , and S 3 , respectively (Figure 6e-h). The lower accuracy of the simulations related to the first survey during S 0 , especially for the water depths, is probably due to the lower quality of the DEM [29], which was produced using a single-beam echosounder rather than airborne LiDAR topo-bathymetric surveying, which is characterized by a high density of points per square metre.

Temporal Evolution of Habitat Heterogeneity
During the whole monitoring period, significant changes were observed throughout the upstream section with an increase in the water depth at S 1 and a reduction in the flow velocity at S 3 . Within the downstream section, significant fining of the bed grain size was observed over the whole monitoring period, and slightly significant changes in the water depth and flow velocity were observed at S 1 (Figure 7a-c). In the impacted section, significant changes in all three physical parameters were observed over the whole monitoring period with decreases in the flow velocity and bed grain size and an increase in the water depth (Figure 7a-c). The results of the paired Wilcoxon test performed between each state and the initial state (S 0 ) for each physical variable are indicated as follows: * = p < 0.05, ** = p < 0.01, *** = p < 0.001. The number of sampling points in the upstream, impacted and downstream sections is 11, 74 and 122, respectively. Boxplots (d) of the nEHMID m and (e) of the CSD index calculated for each cross-section, state and section. The number of cross-sections for the upstream, impacted, and downstream sections is 1, 6, and 11, respectively. Blue, red, and green boxplots correspond to the upstream, impacted, and downstream sections, respectively. S 0 , S 1 , S 2 , and S 3 correspond to each monitoring state.
Changes in the physical parameters along the impacted section induced progressive increases in the CSD and nEHMID m values. The mean CSD was equal to 0.034, 0.037, 0.043, and 0.060 for S 0 , S 1 , S 2 , and S 3 , respectively, and the mean nEHMID m was equal to 0.039, 0.046, 0.051, and 0.078 for S 0 , S 1 , S 2 , and S 3 , respectively. Relatively small changes occurred in both the CSD and the nEHMID m values along the downstream section with mean CSD values equal to 0.038, 0.032, 0.031, and 0.039 for S 0 , S 1 , S 2 , and S 3 , respectively, and mean nEHMID m values equal to 0.045, 0.039, 0,038, and 0.047 for S 0 , S 1 , S 2 , and S 3 , respectively. The upstream section showed a singular pattern comprising a decrease followed by an increase (Figure 7e). Statistically significant linear relationships were found (at a 95% confidence interval) between the partial diversity of the water depth and flow velocity and the CSD index, whereas lower significance was found for the partial diversity of the grain size (Figure 8a-c). The relationship between the nEHMID m and CSD was almost stronger than these relationships (Figure 8d).

Temporal Evolution of Fish Habitat Suitability
For the whole monitoring period and throughout the whole study reach, the highest values of the SHSI (>0.50) were observed for barbel at both the adult and the juvenile life stages and for nase at the adult life stage (Figure 9a). In contrast, the lowest values of the SHSI (<0. 45) were observed for all other species on the three sections, except for salmon at the juvenile stage in the downstream section, for which the SHSI ranged from 0.45 to 0.66.
Significant changes in the TCR values were calculated for both the impacted and the downstream sections (Figure 9b). In the impacted section, the TCR values were positive for both eel and juvenile nase species with a maximum at S 3 equal to +149% and juvenile nase with a maximum at S 1 equal to +57%. At S 2 , positive TCR values in the downstream section were related to juvenile and adult salmon species with maxima equal to +78% and +48%, respectively, and adult grayling with a maximum of +57% at S 3 ( Figure 9b). Globally, during the monitoring period, the TCR values calculated within the upstream section were stable, except for juvenile grayling during S 2 with a TCR value equal to −52% (Figure 9b).

Channel Responses and Habitat Heterogeneity
The results show that the implementation of the two groynes induced significant channel responses. Two mid-bars were created downstream of the two groynes due to the reduction in the flow velocities and the formation of countercurrents downstream of the groynes. Sediments were generated by the scour of the two pools and by bank erosion, both following the reduction in the area of the wetted section induced by the two groynes, while other sediments were produced from residual bedload transport within the main channel upstream of the restoration site [39]. These channel responses are in accordance with previous modeling and field reports [45,61,62]. Field monitoring showed that the flow velocities were reduced along the impacted section during the whole monitoring period (four years) while maintaining the same range of flow velocities as measured during the pre-restoration period for both the upstream and the downstream sections (Figure 7b). Concerning the bed grain size, we observed significant fining of the bed grain size in both the impacted and the downstream sections. Finer grain sizes may be related to the deposition of fine sediments provided by bank erosion [43] and/or by the scouring of the two pools (Figure 7c). These finer particles may also originate from the upstream reaches due both to probable local bank erosion and to break-up of the pre-existing bed armor layer following the intense flood of May-June 2013 (Q 15 ) [42]. Indeed, the break-up of this armor layer induces the entrainment of the subsurface layer, which contains particles that are finer than those within the surface armor layer [63]. Along the regulated lower Ebro River, the authors in [64] showed that the bed armor layer was broken up during intensive flood events and was subsequently re-established during lower flood events. The authors in [65] and [39] further demonstrated that residual geomorphological processes occurred along this by-passed section during intensive flood events, as shown by bed elevation changes.
This study indicates that the transverse groynes significantly increased the physical habitat heterogeneity, whereas few changes occurred in the upstream and downstream sections (Figures 5 and 7). The empirical relationship found between the nEHMID m and CSD indicates that habitat heterogeneity increases as the channel cross-section becomes more diversified (Figure 8d). This finding validates the first hypothesis that the implementation of transverse groynes would enhance habitat heterogeneity in the impacted section. The increase in the nEHMID m values along the impacted section suggests that the geomorphological effects of this restoration are still effective after several floods and may persist as long as the two groynes remain, probably over several decades (Figure 5d).
In such large and heavily impacted by-passed systems within which habitat improvement is critical, such experiences question the opportunity of introducing new groynes in specific reaches to enhance habitat suitability. Such active restoration measures have shown positive effects on fish responses [17,18]. These results raise questions regarding environmental ethics related to the increased use of artificial features to enhance ecological conditions. To restore geomorphic processes in a more sustainable way and achieve improvements to aquatic and riverine habitats efficiently, the authors in [39,42,66] recommended combining gravel augmentations and channel widening. More specifically, channel widening should be performed downstream of gravel augmentations, promoting the deposition of some injected sediments and potentially diversifying aquatic and riverine habitats [41,67]. Several methods may be tested to perform channel widening, such as the removal of bank protection, the partial removal of historical groyne fields and/or the excavation of floodplain areas connected to the channel. In the last option, excavated sediments should be stored temporarily on an island and injected upstream periodically depending on the hydrological regime. However, such operations, which must avoid any risks to human activities (e.g., navigation and flood control), depend strongly on the financial resources available and local technical characteristics. Furthermore, some small interventions may also be tested to enhance local fish habitats [14,16] to reconfigure existing historical groynes and/or introduce wood structures [68].

Temporal Evolution of Fish Habitat Suitability
From four landscape metrics, a synthetic metric, namely, the SHSI, was developed to estimate the quality of fish habitat suitability by integrating changes in the riverscape in terms of both the composition and the configuration, which are rarely studied. It is essential to estimate these factors, notably in the context of river restoration, to integrate all possible riverscape changes because ecological processes at the population and community levels are highly affected by the composition and configuration of a riverscape [32]. Our results highlight that the SHSI metric can be adopt to evaluate habitat changes because it integrates physical variables that are key indicators of habitat diversity (Figures 7e and 9a).
This study demonstrates that the habitat suitability of the studied reach was higher for a larger fish species than before the restoration work. Indeed, an improvement for species that like low flow velocities and low bed grain size conditions occured during the monitoring period on the impacted section (Figures 3 and 9). These results are explained by the reduction in flow velocities induced by the two artificial groynes (Figure 7b) and the deposition of fine sediments downstream of those groynes [69]. Taking into account the whole study reach, the restoration clearly increased the habitat suitability for a larger number of native fish species than in the pre-restoration conditions ( Figure 9).
Indeed, the results of electric fishing conducted along the impacted section showed that the fish communities captured were composed mostly of rheophilic species [43], and eel individuals were present along the impacted section for the first time after the first post-restoration flood of May-June 2013. Interestingly, the results of field fish surveys are in accordance with the results obtained from fish habitat models. Furthermore, these results should be consolidated by further field surveys of fish communities.
An unexpected result was the improvement in the fish habitat suitability for salmon in the downstream section due to sustainable fining of the bed grain size over the four years (Figures 7c and 9a,b). This finding highlights that the bed grain size is as important as the hydraulic conditions for improving the habitat suitability of targeted species. Therefore, in some cases, heavy restorations could be avoided by focusing on the bed grain size rather than only on the hydraulic conditions. Chardon et al. (submitted) [70] showed the ability of gravel augmentations to restore potential spawning areas for rheophilic species due to the preferential deposition of the injected sediments on riffles. This finding also shows that it is important to consider grain size issues in habitat suitability approaches because the bed grain size is a critical ecological parameter that notably concerns lithophilous species and can be improved by specific restoration measures, such as providing a coarse sediment supply from artificial gravel augmentations or bank erosion [26,71].

Eco-Hydraulic Modelling, a Tool for River Restoration Design
Currently, a large amount of financial resources is dedicated to river restoration. The authors in [10] estimated that the restoration of 1 ha of a European river costs, on average, 310,000 e. Furthermore, despite the high costs of restoration actions, several authors highlighted that restoration projects have no clear goals and suffer from a lack of planning [72]. With the introduction of or recent developments in remote sensing methods such as airborne topo-bathymetric LiDAR and hydraulic open-source software, it is now possible to establish two-dimensional hydraulic models and characterize fluvial habitats at finer scales along larger reaches than was previously possible. Hydro-ecological models are a powerful tool for restoration planning to (i) evaluate the degradation states of river channels, (ii) test several restoration designs, and (iii) define the best among those designs according to restoration objectives [73].
Nevertheless, these models have some shortcomings. Most studies of fish habitat modeling focus on hydraulic parameters because the flow conditions are generally considered the most important factors that influence the distributions of fish populations and control, at least partly, some other environmental parameters (e.g., water temperature, water oxygenation, and pollutant concentrations) [74,75]. However, these parameters may limit other factors because the influences of physical parameters decrease in the downstream direction [76], and their influences may increase in the future in the context of global change and increasing land use pressure. Although fish species are able to adapt to climate change by modifying their phenology, physiology, and distribution [77,78], a delay can be observed [79]. Moreover, several species may not be able to adapt to environmental changes [78], which could lead to declines in population and profound community modifications. In the same way, water physico-chemistry and some biotic factors (e.g., food availability, competition, predation, parasitism, and invasion) may also be significant drivers of the fish population distribution and abundance [80][81][82]. Supplementary research must be conducted to better understand the impacts of abiotic, biotic, and chemical factors and their inter-relations with the repartition and abundance of fish populations to sustain them in a changing world.

Conclusions
This paper highlights the advantages of using a riverscape approach rather than a classic approach in fish-habitat modeling. This riverscape approach allowed us to evaluate landscape dynamics in terms of the changes in their configuration and composition, which are crucial in the context of river restoration. The restoration work induced a rise in habitat suitability for a large number of native fish species compared with the pre-restoration conditions. In addition, this study reveals that the improvement of grain size conditions can be as efficient as improvements to flow conditions at enhancing the habitat suitability for some targeted fish species, notably for salmon species. Finally, the implementation of groynes could be a strategy for improving fish habitats, but only where more functional and natural restoration options are impossible.