Research Paper
Villemonte's approach: a general method for modeling uniform and nonuniform performance in stepped fishways
L'approche de Villemonte : une méthode générale de modélisation de la performance uniforme et non uniforme dans les échelles à poissons
^{1}
Centre for Biorobotics, Tallinn University of Technology,
Akadeemia tee 15A111,
12618
Tallinn, Estonia
^{2}
Department of Hydraulics and Hydrology, ETSIIAA, University of Valladolid,
Avenida de Madrid 44, Campus La Yutera,
34004
Palencia, Spain
^{⁎} Corresponding author: juan.fuentes@ttu.ee
Stepped fishways are the most popular solutions to enable the free movement of fish fauna through weirs and dams. Given the flow variation of rivers throughout the year, successful fish migration through stepped fishways relies on the accurate discharge calculation and their modeling under variable scenarios. This study aims to propose a general method for flow and water level calculation of stepped fishways, unifying different findings in specialized literature. To achieve this purpose, the relation defined by Villemonte is used and tested under laboratory and field case studies. This study shows that the hydraulic behavior of a wide range of stepped fishway typologies can be explained based on a single equation, as well as the need of calibration of the coefficients involved in this equation for different subtypes. Furthermore, the proposed method enables the water level modeling under variable boundary conditions, which in turns allows the analysis of stepped fishways hydraulic performance under different river scenarios. The comparison of the hydraulic parameters in the fishways with the physical capacities and preferences of fish will contribute to the fulfillment of their main objective: allow free movement of fish fauna.
Résumé
Les échelles à poissons sont les solutions les plus populaires pour permettre la libre circulation de la faune piscicole à travers les seuils et les barrages. Compte tenu de la variation de débit des cours d'eau tout au long de l'année, la migration réussie des poissons par ces passes à poissons repose sur le calcul exact de la décharge et sa modélisation dans des scénarios variables. Cette étude a pour but de proposer une méthode générale de calcul du débit et du niveau d'eau des échelles à poissons, synthétisant différents résultats de la littérature spécialisée. Pour atteindre cet objectif, la relation définie par Villemonte est utilisée et testée dans des études de cas en laboratoire et sur le terrain. Cette étude montre que le comportement hydraulique d'une large gamme de typologies de poissons en gradins peut être expliqué sur la base d'une seule équation, avec besoin d'étalonnage des coefficients impliqués dans cette équation pour différents soustypes. De plus, la méthode proposée permet de modéliser le niveau de l'eau dans des conditions aux limites variables, ce qui permet, à son tour, d'analyser la performance hydraulique des passes de poissons en escaliers selon différents scénarios de rivières. La comparaison des paramètres hydrauliques dans les passes à poissons avec les capacités physiques et les préférences des poissons contribuera à la réalisation de leur objectif principal : permettre la libre circulation de la faune aquatique.
Key words: discharge coefficient / flow measurement / vertical slot fishway / pool and weir fishway / naturelike fishway
Mots clés : coefficient de décharge / mesure de débit / passe à fentes verticales / passe en bassins / passe de type rivière artificielle
© J.F. FuentesPérez et al., Published by EDP Sciences 2017
This is an Open Access article distributed under the terms of the Creative Commons Attribution License CCBYND (http://creativecommons.org/licenses/bynd/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. If you remix, transform, or build upon the material, you may not distribute the modified material.
1 Introduction
Humans have found in rivers a source to satisfy many of their basic necessities, which have resulted in their geomorphological and ecological alteration (Nilsson et al., 2005). One of the most notable alteration in rivers is the installation of crosssectional structures (e.g. weirs and dams) to satisfy water and energy requirements or for flood control. These structures fragment the stream and can block the movement of some animals such as fish, which require different environments to complete their life cycles (Lucas et al., 2001; Branco et al., 2013). In recent years, ensuring undisturbed fish migration has become a key component of watershed restoration (Santos et al., 2014) and the installation of fishways is one of the most widely adopted solution in order to achieve this objective.
There are many types of fishways. The most common ones consist of a succession of crosswalls in a sloped channel, namely, stepped fishways or fish ladders (Fig. 1) (e.g. vertical slot fishway (VSF) (Rajaratnam et al., 1986; Larinier, 2002b; FuentesPérez et al., 2014), pool and weir fishway (PWF) (Rajaratnam et al., 1988; Larinier, 2002b; FuentesPérez et al., 2016), and steppool naturelike fishway (SPNF) (FAO/DVWK, 2002; Wang and Hartlieb, 2011)). These structures divide the total height of the obstacle (H) in smaller drops (ΔH) in each crosswall to ensure that the hydraulic conditions inside are in the range of the physical capacities of fish fauna and, thus, enable their passage.
Depending on the type of fishway, they can have different kinds (slots, notches, or orifices) and number of connections in the crosswalls between pools, from a single slot, like for example in some VSFs, to multiple combinations in SPNFs.
The discharge and performance of fishways can be modeled in many different ways, which enables the classification of these hydraulic calculations into two big groups: dimensionless relationships (Rajaratnam et al., 1986, 1989; Ead et al., 2004; Yagci, 2010) and classical weir equations (Larinier, 1992; Clay, 1995; Martínez de Azagra, 1999; Boiten and Dommerholt, 2006; Krüger et al., 2010). The first group of equations is only useful when ΔH is equal to the topographic difference between pools (ΔZ) (i.e. same water depth in all pools) (Rajaratnam et al., 1986). This is known as uniform regimen, which is difficult to achieve in field conditions due to the temporal variability of river flow (FuentesPérez et al., 2016; Marriner et al., 2016). Regarding classical weir equations, the most commonly used equations are the orifice equations derived from Torricelli's law (Torricelli, 1644) as well as Poleni's weir equation (Poleni, 1717). An accurate selection of the discharge coefficients for these equations is vital in order to achieve precise results, both under uniform (ΔH = ΔZ) as well as under nonuniform (ΔH ≠ ΔZ) conditions. In order to use these equations under both conditions, it will be necessary to consider the water level upstream () and downstream () of the crosswall (FuentesPérez et al., 2014).
Nonuniform performance in fishways is usually classified in two main different water level profiles (Rajaratnam et al., 1986): M1 or backwater profile, which produces higher mean depth (h_{0}) and smaller drops (ΔH < ΔZ) in the downstream pools of the fishway, and M2 or drawdown profile, which contrary, generates lower h_{0} and higher drops (ΔH > ΔZ) in downstream pools. It is likely to occur in field conditions due to the variable hydrological regimen of river or small deviances in the construction of the fishway (FuentesPérez et al., 2016; Marriner et al., 2016). This will modify the hydrodynamics of the flow inside the fishway, which may lead to an incompatibility with fish fauna preferences or capabilities, affecting the efficiency of fish passage (FuentesPérez et al., 2016; SanzRonda et al., 2016). For instance, M1 profiles may improve the passability due to the lower velocities in the crosswalls or the reduction on the volumetric power dissipation (VPD) within the pools, but it may decrease the attractivity due to reduction of velocity in the most downstream crosswalls. Contrary, M2 profiles may led to more attractive scenarios but may generate too demanding drops to be surpassed or conditions with too high VPD in downstream pools.
Thus, nonuniformity must be considered in fishway research. However, the simplification of the study of fishways is a widespread problem and, in many cases, their performance is modeled according to uniform conditions (Rajaratnam et al., 1986; Wu et al., 1999; Puertas et al., 2004; Bermúdez et al., 2010; Wang et al., 2010; Tarrade et al., 2011; among others).
Using some of the most extended fishway guidelines (Larinier, 1992; FAO/DVWK, 2002; Krüger et al., 2010) it is possible to explain, at least partly, the nonuniform performance. Likewise, some recent works have shown that it is possible to use classical weir equations together with the submerged weir discharge coefficient (C_{s}) proposed by Villemonte (1947) to model fishway performance under uniform and nonuniform scenarios in specific subtypes of VSF (FuentesPérez et al., 2014) and PWF (FuentesPérez et al., 2016).
In the present paper, we study the use of Villemonte's equation as a general discharge coefficient definition for flow and water level calculation of all stepped fishway types under field conditions. This is achieved by calibrating the discharge coefficient for the different fishway types (PWF, VSF and SPNF) and subtypes (different morphologies within the types) studied in the field and in specialized literature, which will enable us to validate this general methodology. The main contributions of this paper are: to (i) prove the usefulness of Villemonte's equation as a general method for the estimation of stepped fishway discharge and water levels under uniform and nonuniform scenarios; (ii) adjust the coefficients of this equation for flow measurement in stepped fishway types and subtypes proposed by the specialized literature as well as field cases; and (iii) show that this equation is able to unify findings in different scientific references.
Fig. 1 Examples of sections of stepped fishways. Each type (a) vertical slot fishway, (b) pool and weir fishway and (c) steppool naturelike fishway can have different subtypes according to their morphology and connections. 
2 Material and methods
2.1 Theoretical background: recommended discharge equations
The performance and discharge of stepped fishways will be explained by means of classical weir equations. The equation for weirs proposed by Poleni (1717) (Eq. (1)) will be used to describe the flow (Q) through notches (Larinier, 1992; Martínez de Azagra, 1999; FAO/DVWK, 2002) and slots (FAO/DVWK, 2002; Krüger et al., 2010). On the other hand, the orifice equation derived from Torricelli's law (Torricelli, 1644) (Eq. (2)) will be used to describe the Q through orifices (Larinier, 1992; Martínez de Azagra, 1999; FAO/DVWK, 2002; Boiten and Dommerholt, 2006). (1) (2) where g stands for the gravity, h_{1} is the water level upstream the crosswall deducting the sill height (p), b is the width of notches, slots or orifices, a is the height of the orifice and C (C_{s} and C_{o}) stands for the discharge coefficients.
Multiple fishways studies have demonstrated that C_{o} can be considered independent of the different water level variables involved (ΔH, h_{2} and h_{1}, where h_{2} is the water level downstream the crosswall deducting p) (Brater et al., 1996; Boiten and Dommerholt, 2006; FuentesPérez et al., 2016). This property is derived from equation (2) and, for this reason, it can be considered that C_{o} will only vary according to the geometrical variables of the orifice (dimensions, shape and crosswall thickness).
For slots and notches, Villemonte's submerged coefficient will be used as a generic discharge coefficient (Villemonte, 1947). Villemonte derived his coefficient assuming that Q is the difference between the freeflow discharge due to h_{1} (Q_{1}) and the counterflow discharge due to h_{2} (Q_{2}) (Fig. 2 and Eq. (3)). (3)
The experiments carried out by Villemonte confirmed this hypothesis, Q/Q_{1} (C_{s}) is related functionally to 1 − Q_{2}/Q_{1}. Considering this and replacing Q_{1} and Q_{2} by equation (1), the final form of C_{s} is obtained (Eq. (4)). (4) where β_{0} and β_{1} are coefficients which depend on the geometry of the slot or notch, pool dimensions, and the discharge equation used. Villemonte concluded that β_{1} was equal to 0.385 in all his experiments for submerged weirs, however he recommended its calibration for different setups. β_{0} can be considered as a constant approximation of a free discharge coefficient.
Fig. 2 Sketch of Villemonte's assumption for submerged flows. Q_{1} is the free flow discharge in positive direction controlled by h_{1} and Q_{2} is the counterflow discharge due to h_{2}. p = 0 in case of slots. 
2.2 Hypothesis testing through case studies
Data from peerreviewed scientific studies, grey literature and from our own field experiments were used to test the hypothesis: Villemonte's equation can be adapted to model different types of stepped fishway under uniform and nonuniform conditions.
A systematic search was carried out following the PRISMA (Preferred Reporting Items for Systematic Reviews and MetaAnalyses) guidelines (Moher et al., 2009) (S1 Appendix). This systematic search was carried out in Google Scholar as well as in Web of Science and Scopus from conception to September 2016 for articles and grey literature on VSFs, PWFs and SPNFs. To standardize and minimize errors only studies containing raw data and appropriate geometrical descriptions were considered for calibrating purpose. This is rather difficult to find and, aside data from our own field tests, only references for VSFs were found useful (Rajaratnam et al., 1986, 1992; FAO/DVWK, 2002; Puertas et al., 2004) (Tab. 1). Likewise, as only English has been used for the search, this might have introduced a language bias.
Regarding field studies (8 in total), all the fishways are located in the basin of the Duero River, in NorthCentral Spain (Fig. 3). The Department of Hydraulics and Hydrology of the University of Valladolid was involved in the design of all of them. The information about location, characteristics and geometry of each fishway can be found in S2 Appendix.
Summary of collected data through references and field experiments.
Fig. 3 Field cases location. Three types of fishways were studied in field: vertical slot fishway (VSF), pool and weir fishway (PWF) and steppool naturelike fishway (SPNF). Localization coordinates are defined in S2 Appendix. 
2.3 Experimental arrangement
The following experimental procedure was used for all field experiments. The geometrical parameters of each structure were measured by topographic surveying with a total station Leica TC307 to a precision of 0.001 m. Smaller details, such as orifice, notch and slot dimensions were measured by means of metal rulers to the same precision level.
The discharge of each scenario (different studied water level configurations, Tab. 2) through the structures was obtained by chemical gaging (precision of 0.004 m^{3}/s) using Rhodamine WT as tracer (Martínez, 2001). In cases where was possible, boundary conditions were artificially modified in the fishways in order to increase the number of studied scenarios. For example, when there was a gate situated upstream, different discharges were measured and, if possible, the water level downstream was also increased and decreased by modifying the section of the last control structures, creating new scenarios (Tab. 2) (see S2 Appendix for the geometrical information of the field cases and S3 Data Appendix for discharge, water levels and boundary conditions of each studied scenario).
Water levels were measured to a millimetric precision in each of the pools, by means of metal rulers installed where water surface was found more stable (downstream of the crosswall and 0.40 m from the sidewall opposite to the notch (0.15B) or slot (0.25B)). Water level oscillations were recorded for 8 s using a camera (Canon EOS 600D) with a sampling rate of 25 Hz; in all cases after ≤2 s (50 samples) a stable mean value for the water level was obtained. For SPNF, due to the naturelike morphology and, thus, the absence of clear reference points for the metal rulers, water levels were measured in each of the corners of the pools using the total station. The mean value was used for calculations.
Summary of the field scenarios for each fishway type (see S3 Data Appendix for detailed information about each studied scenario).
2.4 Data fitting and validation
Data from the literature review and field experiments were combined for each of the different stepped fishway types considered (VSF, PWF or SPNF) and within them for the different possible subtypes (e.g. VSF designs from No. 1 to No. 18 of Rajaratnam et al., 1986). Fishway subtypes are defined as fishways that, within the same type, have substantial differences in the configurations which result in a modification of the flow pattern. A global data table was created for each subtype with the necessary variables to study the proposed relations (i.e. type, subtype, h_{2}, h_{1}, ΔZ, Q and b and a, for each of the connections (Fig. 1) if applicable). Curve fitting to calculate Villemonte's coefficients (β_{0} and β_{1} in Eq. (4)) for each study case was undertaken in R, release 3.3.1 (R Core Team, 2016).
Following this, all fitted coefficients for each type and subtype of fishway were evaluated using root mean square errors (RMSEs) and determination coefficient (R^{2}) as well as graphically. The proposed method was compared with other available methods of calculation for each fishway type. When different number of fitting parameters were involved between two different fishway calculation methods Bayesian Information Criterion (BIC) was used to compare the performance (Priestley, 1981; Kass and Wasserman, 1995).
In VSFs analysis, as an example, to demonstrate the applicability of the method here proposed, fitted equations were used to predict water level distributions (predicted data). In these cases, mean relative error (MRE) with respect to the measured real data (observed data), was used as descriptor of the achieved accuracy.
3 Results and discussion
3.1 Vertical slot fishways
Considering the number of references by fishway, VSFs seem to be one of the most studied and extended fishway types in the world. In 1986 and 1992 Rajaratnam et al. conducted the largest serial study to date of this type of structures. The VSF subtypes defined in these papers (18 in total) have been used as a reference for successive works (Puertas et al., 2004; FuentesPérez et al., 2014) (for a complete description of the geometry of these subtypes please see the original references (Rajaratnam et al., 1986, 1992)). These authors defined dimensionless relationships to explain VSF performance (Eqs. (S4.1) and (S4.2) in S4 Appendix), including slope (S) as a variable and only mean water level in the pool (h_{0}) as a water level descriptor. This has been proven to be useful for uniform scenarios; however, it is invalid to estimate flow and water level distributions in nonuniform scenarios (Rajaratnam et al., 1986; FuentesPérez et al., 2014).
Table 3 summarizes the number of observations, the origin of data as well as the fitted coefficients for equation (4) for each subtype of VSF defined by Rajaratnam et al. (1986, 1992) after literature and experimental data combination. In all cases, the use of equation (4) as discharge coefficient enables to explain precisely the observed variability as shown by the R^{2} value. It is worth mentioning that the lower R^{2} for subtype 1 may be explained by the higher variances of the four field studies included, compared to the laboratory experiments (Rajaratnam et al., 1986, 1992; Puertas et al., 2004) (Fig. 4a).
Figures 4a and b show the fits for some common VSF subtypes using Villemonte's equation. This equation is able to unify findings and observations of different studies. In addition, as shown in Figure 4a, the same fit could be used in subtypes 2 and 1 despite the presence of a small sill (p = 0.15 m) in the former.
Although in equation (4), S is not considered, the relation between h_{2} and h_{1} alone is able to account for the observed variance caused by different S configurations in the same subtype. This can be explained by the fact that S will determine the relation between both variables (h_{2} and h_{1}) under uniform conditions (S∙(L + e) = ΔZ = ΔH = h_{1}−h_{2}, where L is the length of the pool and e is the thickness of the crosswall) (FuentesPérez et al., 2014). For this reason, it is possible to use the equation independently of the selected S, considering only the subtype and the observed water level distributions. It is worth mentioning that in cases where S = 0 (fishways without slope) (Bice et al., 2017), the performance will always be nonuniform (Musall et al., 2015) and ΔH will only depend on the boundary conditions used. Although equation (4) allows its calculation, dimensionless relationships (Eqs. (S4.1) and (S4.2) in S4 Appendix) (Rajaratnam et al., 1986, 1992) will lead to an indetermination in these cases, as S is a divider in equation (S4.1).
Figure 5a shows the potential of the proposed equation by correctly estimating the performances observed by Rajaratnam et al. (1986). In contrast to dimensionless relationships proposed by Rajaratnam et al. (1986), as h_{2} and h_{1} are involved in the proposed method, considering the geometrical parameters of the fishway, the bottomup resolution of equation (1) and equation (4) (i.e. between two consecutive crosswalls, i and i+1 (cf. Figs. S2.2 and S2.4 in S2 Appendix): and (FuentesPérez et al., 2014)) can explain the water level distribution in each of the pools under uniform and nonuniform conditions accurately (Tab. 4).
The correct estimation of water level distribution is vital to assess the performance of fishways, as from this, other variables which are directly related with fish preferences and swimming abilities (Larinier, 2002a, 2002b), such as volumetric power dissipation (VPD = ρ · g · ΔH · Q/(B · L · h_{0}), where ρ is the water density and B is the width of the pool) (FAO/DVWK, 2002; Towler et al., 2015) or maximum velocity in the crosswall () (Rajaratnam et al., 1986; Larinier, 2002b) can be estimated. Therefore, this method will enable the analysis of the possible consequences in fish passage of different hydraulic scenarios and, thus, it will allow the preliminary evaluation of the passage and detection of possible problems during the whole hydrological period.
In recent years, different equations or equation combinations have been proposed to explain the performance of the discharge coefficient in VSFs under uniform and nonuniform scenarios (Krüger et al., 2010; FuentesPérez et al., 2014; Marriner et al., 2016). However, all of them present limitations that can be overcome with the use of equation (4).
For instance, Krüger et al. (2010) proposed the use of a similar equation to equation (4) but allowing the constant exponent 1.5 to be adjusted for each VSF subtype (Eq. (S4.3) in S4 Appendix). This only provides slightly higher R^{2} in the 46% of cases (in the 54% of cases Eq. (4) provides higher R^{2}); however, in the 96% of cases 1.5 is in the 95% of confident interval of the exponent, which consequently seems to indicate that the fit of this variable is unnecessary. Likewise, if the number of free parameters (three against two) is taken into account and the collected data is applied to fit both equations, equation (4) will produce better fits in 67% of studied cases (BIC_{Eq.(4)} < BIC_{Eq.(S4.3)}) and in 40% of all cases the difference will not be significant (BIC_{Eq.(4)} − BIC_{Eq.(S4.3)} < 2), which seems to indicate an overfitting problem in the equation proposed by Krüger et al. (2010). Moreover, Villemonte's approach enables the comparison of different fishways, as β_{0} and β_{1} in equation (4) provide, directly, information about hydrodynamic performance. β_{0} provides information about the horizontal asymptote to be reached by the fit and β_{1} (indicative of the curvature of the fit) provides information about the evolution of the coefficient with the submergence.
In FuentesPérez et al. (2014), as an alternative to the combination of equations (1) and (4), equation (2) together with equation (4) was suggested for flow estimation in VSFs (e.g. Figs. 4c and 5c). Equation (2) has been widely recommended by specialized references as it has a physical justification (Clay, 1995; Martínez de Azagra, 1999; Larinier, 2002b; Marriner et al., 2016) and as it collects part of the variability caused by the nonuniform performance itself (without Eq. (4)). This allows a better flow calculation than with equation (1) in the first states of nonuniformity when a constant discharge coefficient is used (Fig. 5b). However, as can be seen from Figures 4c, 5b, 5c and Table 4, the use of equation (4) is also necessary to collect all the variability produced by nonuniformity and to accurately calculate the flow and water level distribution. When equation (4) is used, good results are expected with both equations (Tab. 4). Nevertheless, the use of equation (2) instead of equation (1), can also reduce the interpretability of equation (4), as the variability due to nonuniformity is distributed in two different equations (Eqs. (2) and (4)), instead of only one (Eq. (4)).
Fig. 4 Examples of discharge coefficient fits for most relevant (in terms of data points) VSF subtypes. (a) Subtypes 1–3 (fit coefficients and R^{2} in Tab. 3). (b) Subtype 6 (fit coefficients and R^{2} in Tab. 3). (c) Fit of equation (4) for subtype 1 to be used in equation (2) (β_{0} = 0.509, β_{1} = −0.248, R^{2} = 0.587, RMSE = 0.0458). 
Fig. 5 Nonuniform performance in a VSF of subtype 3, scenario observed (grey) by Rajaratnam et al. (1986) and simulated (black) by different equation combinations. M1: h_{2} = 2.712 m and Q = 0.66 m^{3}/s; M2: h_{2} = 0.931 m and Q = 0.66 m^{3}/s; uniform: h_{2} = 1.416 m and Q = 0.66 m^{3}/s. (a) Equations (1) and (4). (b) Equation (2) and constant coefficient (C = 0.56). (c) Equations (2) and (4). 
3.2 Pool and weir fishways
Villemonte's equation was originally proposed for submerged weirs and, therefore, it has been widely used to calculate the flow over PWF. Although Villemonte adjusted the equation for sharpcrested weirs with specific experimental conditions and suggested new fits for different conditions, in practice, the originally proposed coefficients have been indiscriminately used in PWFs (Larinier, 1992; Martínez de Azagra, 1999; FAO/DVWK, 2002). The utility of equation (4) and the convenience of specific fits for a PWF subtype has been previously shown in FuentesPérez et al. (2016).
The notch of PWFs can perform under two different flow regimen conditions, plunging (or free flow) and streaming (or submerged flow), depending on the influence of h_{2} over h_{1} (Rajaratnam et al., 1988). In the first case, h_{2} = 0 (, see Fig. 2) allowing the free flow, while in the later performance h_{2} > 0 (, see Fig. 2) and, thus, h_{2} slowdowns the flow due to h_{1}. These performances can be found independently or mixed within a fishway. Therefore, the proposed equations must be able to describe both possible conditions.
Usually, and for practical reasons, a constant discharge coefficient is used for notch in plunging performance (Larinier, 1992; Martínez de Azagra, 1999; FAO/DVWK, 2002). This is because the coefficient only changes significantly for low water levels over the sill which is not a typical performance of fishways (Rehbock, 1929; Ead et al., 2004; FuentesPérez et al., 2016). For this reason, β_{0} in equation (4) can be considered as a constant approximation of the free discharge coefficient (when free flow h_{2}/h_{1} = 0 consequently C_{s} = β_{0}) and, thus, equation (4) can be used as a general descriptor of plunging and streaming performances. For studies focused on low discharges over fishways, and usually outside their operation range, the use of a variable β_{0} is recommended (cf. FuentesPérez et al. (2016)).
PWFs besides a notch also tend to have an orifice, as an alternative passage for fish fauna (Silva et al., 2009; Branco et al., 2013) or to facilitate fishway cleaning (Larinier, 2002b). Thus, the discharge of the orifice must be considered. As previously shown, the C_{o} for orifices will be constant. Therefore, 0.876 has been used for C_{o} in all field scenarios to calculate the parameters of equation (4). This value is within the coefficient range recommended by specialized references (Martínez de Azagra, 1999; FAO/DVWK, 2002; Larinier, 2002b; Boiten and Dommerholt, 2006) and was calculated for the subtype under study in a systematic field study (FuentesPérez et al., 2016). Nevertheless, appropriate specific constant C_{o} values for submerged orifices can be found in Brater et al. (1996).
Considering all of the above, i.e. constant C_{o} and β_{0}, Figure 6 shows the C_{s} distribution for the studied PWFs cases. As shown in the figure, considering the β_{1} proposed by Villemonte (β_{1} = 0.385), as suggested by specialized references, and fitting only β_{0}, provides a good result; however as recommended by Villemonte, a specific fit offers significant better results (R^{2} of 0.755 versus 0.683, RMSE 0.0474 versus 0.0538 and BIC_{Eq.(4)} < BIC_{Eq.(4),β1 = 0.385}).
The studied field cases involve a b value for the notch in the range of 0.20–0.40 m and a side length for the square orifice in the range of 0.15–0.25 m (Tab. S2.2 in S2 Appendix), which demonstrates the utility of the fit for different notch dimensions. However, all the studied cases have similar configuration (i.e. the flow pattern is similar and therefore they can be considered the same subtype) and, thus, as for the VSF, new fits are recommended for different possible subtypes.
A common alternative calculation to the one described above are the dimensionless relations proposed by Rajaratnam et al. (1988) (Eqs. (S4.4) and (S4.5) in S4 Appendix). These equations have been studied and validated in a number of studies (Ead et al., 2004; Yagci, 2010). However, each of the equations proposed by Rajaratnam et al. (1988) is defined for a possible regimen (plunging or streaming), not allowing mixed scenarios. Furthermore, as the dimensionless relations proposed by the same authors for VSFs, they depend on S and they only consider h_{0} as water level descriptor, which does not allow their use under nonuniform conditions and generate indeterminations in some cases (e.g. S = 0).
Fig. 6 C_{s} fit for studied PWF and alternative methods proposed by specialized references. In total 180 data points were collected. 
3.3 Steppool naturelike fishway
Most of hydraulic studies on SPNFs have focused on the evaluation of velocity, dissipated energy or friction factors (Pagliara and Chiavaccini, 2006; Wang and Hartlieb, 2011; Oertel and Schlenkhoff, 2012) and, there are no specific studies on discharge and discharge coefficients.
Despite the lack of published studies, design guidelines recommend the use of discharge and discharge coefficient equations. For instance, FAO/DVWK (2002) and Larinier et al. (2006) recommend the use of the equation for submerged flow in broadcrested weirs. This consists of the use of an equation similar to equation (1) together with a drownedflow reduction factor (σ) (Schröder, 1994) (Eqs. (S4.6) and (S4.7) in S4 Appendix) and a reduction constant coefficient (μ) that depends on the geometry of the stones (range from 0.5 to 0.8 if they are sharpedged or rounded stones respectively (FAO/DVWK, 2002)). In the same way, Hegberg et al. (2001) proposed the use of the equation for sharpcrested weirs of Villemonte (Eqs. (1) and (4) with β_{1} = 0.385 and a calibrated β_{0}). Likewise, due to the small sills that are presented in notches in these kind of configurations (mean p = 0.166 m for the studied case), it is a common practice (FAO/DVWK, 2002; Larinier et al., 2006) to consider for all type of connections same discharge coefficient (i.e. same exponents and coefficients in the discharge coefficient equation).
Figure 7 shows both methods fitted for the field case data in a SPNF together with the fit of Villemonte's equation, where both parameters involved (β_{0} and β_{1}) have been recalculated given the different settings from the ones originally used by the author. As shown by the fit for the broadcrested weir method (FAO/DVWK, 2002; Larinier et al., 2006), it assumes that for situations with h_{2}/h_{1} < 0.614 there is no influence of h_{2} in the flow, which does not fit the observed data distribution (RMSE = 0.3433).
When comparing the β_{1} fitted with that proposed in Villemonte's experiments, as recommended by Hegberg et al. (2001), fitting both parameters provide better significant results (RMSE of 0.0573 versus 0.0583 and BIC_{Eq.(4)} < BIC_{Eq.(4),β1 = 0.385}) (Fig. 7). However, the difference is smaller than the one observed for other fishway types (Figs. 4 and 6). This could be explained by the smaller slopes, longer pool length and less turbulence of SPNFs compared with other fishway types, which are closer to the experimental conditions used by Villemonte (without slope (S = 0) and h_{1} and h_{2} measured outside the influence of turbulence). It is noteworthy that the presented result needed to be further studied as only data from one field study has been considered.
In addition, further considerations should be taken into account when overtopping flows occur. In these cases, the crosswalls are completely submerged and the fishway behaves like a channel with continuous roughness in the bottom. Under these conditions, the flow is usually explained by means of relative energy dissipation and friction factors in SPNFs (Oertel and Schlenkhoff, 2012). These variables will depend on the configuration of the boulder arrangements (boulder size, shape and concentration), the slope (Oertel and Schlenkhoff, 2012; Cassan and Laurens, 2016) and also on the submergence ratio (Oertel and Schlenkhoff, 2012; Baki et al., 2014). Likewise, more studies under these experimental conditions are necessary to evaluate the performance of the proposed equations. Nevertheless, these conditions can always be controlled with an adequate design or a control gate upstream of the facility.
Fig. 7 Fits for the studied SPNF and alternative methods proposed by specialized references. In total 45 data points were collected. 
4 Summary and conclusions
This article shows how a single equation, Villemonte's equation, can be used to calculate the flow and water level distribution of the most common type of stepped fishways under different hydrodynamic conditions, which provides a standardized method for the design and calibration of the performance of these structures. In addition, compared to more extended methods, this equation has shown a better performance not only due to its simplicity, better fit or capacity to unify the different references, but also due to its ability to explain uniform and nonuniform scenarios with a single equation.
In general, in order to obtain an optimal solution, it seems necessary to estimate both free parameters involved in equation (4) (β_{0} and β_{1}), which is understandable given the different experimental conditions used by Villemonte in his original studies.
The possibility of modeling water level distribution in variable hydrodynamic scenarios provides us with the opportunity to study the variation of hydraulic parameters (e.g. VPD and V_{max}) within the fishway and ensures the fulfillment of fish fauna preferences in different flow conditions and during the whole hydrological period. Furthermore, if incompatibilities with fish fauna are detected, possible solutions can be tested using the proposed equations. This possibility of optimization in both design and operation will contribute to ensure the longitudinal connectivity in rivers and, thus, facilitate fish migration and species conservation.
In future studies, it would be interesting to validate Villemonte's methodology in other subtypes of fishway, for instance those with nonrectangular shape connections (e.g. triangular, parabolic, etc.).
Nomenclature
b: notch, slot or orifice width (m)
C_{o}: discharge coefficient for orifices (dimensionless)
C_{s}: discharge coefficient for slots and notches (dimensionless)
g: acceleration due to gravity (m/s^{2})
H: total height of the transversal obstacle (m)
h_{0}: mean water level of the flow in the pool (m)
h_{1}: mean water level of the flow in the pool upstream of the crosswall measured from the sill (m)
h_{2}: mean water level of the flow in the pool downstream of the crosswall measured from the sill (m)
: mean water level of the flow in the pool upstream of the crosswall () (m)
: mean water level of the flow in the pool downstream of the crosswall () (m)
Q_{1}: freeflow discharge due to h_{1} (m^{3}/s)
Q_{2}: counterflow discharge due to h_{2} (m^{3}/s)
Q_{design}: theoretical discharge or flow, as calculated in the fishway design project (m^{3}/s)
Q: discharge or flow rate (m^{3}/s)
R^{2}: determination coefficient (dimensionless)
V_{max}: maximum velocity (m/s)
VPD: volumetric power dissipation (W/m^{3})
β_{0}, β_{1}: dimensionless coefficients for equation (4)
ΔH: water level difference between pools or head drop (ΔH = h_{1} – h_{2}) (m)
ΔZ: topographic difference between crosswalls (m)
ρ: density of water (kg/m^{3})
σ: reduction factor for broadcrested weirs (dimensionless)
μ: reduction constant coefficient for broadcrested weirs (dimensionless)
Acknowledgements
The authors would like to thank the owners and staff of the hydropower plants and fish farms associated to field cases (SAVASA, Iberdrola, Molinos de Castilla, Comunidad de Regantes de El Carracillo, Junta de Castilla y León, Piscifactoría de Campoo) for their collaboration during the experimental procedures, as well as Dr. Sara FuentesPérez for her editorial advice and active participation in the revision of this paper.
References
 Baki ABM, Zhu DZ, Rajaratnam N. 2014. Mean flow characteristics in a rockramptype fish pass. J Hydraul Eng 140: 156–168. [CrossRef] (In the text)
 Bermúdez M, Puertas J, Cea L, Pena L, Balairón L. 2010. Influence of pool geometry on the biological efficiency of vertical slot fishways. Ecol Eng 36: 1355–1364. (In the text)
 Bice CM, Zampatti BP, MallenCooper M. 2017. Paired hydraulically distinct verticalslot fishways provide complementary fish passage at an estuarine barrier. Ecol Eng 98: 246–256. (In the text)
 Boiten W, Dommerholt A. 2006. Standard design of the Dutch pool and orifice fishway. Int J River Basin Manag 4: 219–227. [CrossRef] (In the text)
 Branco P, Santos JM, Katopodis C, Pinheiro A, Ferreira MT. 2013. Pooltype fishways: two different morphoecological cyprinid species facing plunging and streaming flows. PLOS ONE 8: e65089. (In the text)
 Brater EF, King HW, Lindell JE, Wei CY. 1996. Handbook of hydraulics. Boston, MA, USA: McGraw Hill. (In the text)
 Cassan L, Laurens P. 2016. Design of emergent and submerged rockramp fish passes. Knowl Manag Aquat Ecosyst 45: 417. (In the text)
 Clay CH. 1995. Design of fishways and other fish facilities. Ottawa, Canada: CRC Press. (In the text)
 Ead SA, Katopodis C, Sikora GJ, Rajaratnam N. 2004. Flow regimes and structure in pool and weir fishways. J Environ Eng Sci 3: 379–390. [CrossRef] (In the text)
 FAO/DVWK. 2002. Fish passes: design, dimensions, and monitoring. Rome, Italy: FAO, 119 p. (In the text)
 FuentesPérez JF, SanzRonda FJ, Martínez de Azagra Paredes A, GarcíaVega A. 2014. Modeling waterdepth distribution in verticalslot fishways under uniform and nonuniform scenarios. J Hydraul Eng 140: 6014016. [CrossRef] (In the text)
 FuentesPérez JF, SanzRonda FJ, Martínez de AzagraParedes A, GarcíaVega A, Martínez de Azagra A, GarcíaVega A. 2016. Nonuniform hydraulic behavior of poolweir fishways: a tool to optimize its design and performance. Ecol Eng 86: 5–12. [CrossRef] (In the text)
 Hegberg CH, Jacobs S, Schlindwein PA, Cohen SJ. 2001. Fish passage in urban streams. Part 2: Natural channel design methods. Wetl Eng River Restor 1–12. (In the text)
 Kass RE, Wasserman L. 1995. A reference Bayesian test for nested hypotheses and its relationship to the Schwarz criterion. J Am Stat Assoc 90: 928–934. [CrossRef] [MathSciNet] (In the text)
 Krüger F, Heimerl S, Seidel F, Lehmann B. 2010. Ein Diskussionsbeitrag zur hydraulischen Berechnung von Schlitzpässen. Wasser Wirtschaft 3: 31–36. (In the text)
 Larinier M. 1992. Passes à bassins successifs, prébarrages et rivières artificielles. Bull Fr Pêche Piscic 45–72: 326–327. (In the text)
 Larinier M. 2002a. Biological factors to be taken into account in the design of fishways, the concept of obstructions to upstream migration. Bull Fr Pêche Piscic 364: 28–38. [CrossRef] (In the text)
 Larinier M. 2002b. Pool fishways, prebarrages and natural bypass channels. Bull Fr Pêche Piscic 364: 54–82. [CrossRef] (In the text)
 Larinier M, Courret D, Gomes P. 2006. Guide technique pour la conception des passes à poissons “naturelles”. Rapport GHAAPPE RA.06.05V1, Agence de l'Eau Adour Garonne, Compagnie Nationale du Rhône, Conseil Supérieur de la Pêche, Groupe d'Hydraulique Appliquée aux Aménagements Piscicoles et à la Protection de l'Environnement, Institut de mécanique des fluides de Toulouse, 66 p. (In the text)
 Lucas MC, Baras E, Thom TJ, Duncan A, Slavík O. 2001. Migration of freshwater fishes. Oxford, UK: Wiley Online Library, 440 p. [CrossRef] (In the text)
 Marriner BA, Baki ABM, Zhu DZ, Cooke SJ, Katopodis C. 2016. The hydraulics of a vertical slot fishway: a case study on the multispecies VianneyLegendre fishway in Quebec, Canada. Ecol Eng 90: 190–202. [CrossRef] (In the text)
 Martínez E. 2001. Hidráulica fluvial: principios y práctica. Madrid, Spain: Bellisco. (In the text)
 Martínez de Azagra A. 1999. Escalas para peces. Palencia, Spain: ETSIIAA, Universidad de Valladolid. (In the text)
 Moher D, Liberati A, Tetzlaff J, Altman DG, The PRISMA Group. 2009. Preferred reporting items for systematic reviews and metaanalyses: The PRISMA statement. PLoS Med 6: e1000097. [CrossRef] (In the text)
 Musall M, Oberle P, CarbonellBaeza R, FuentesPérez JF, Tuhtan JA, Nestmann F. 2015. Beitrag zu detaillierten Analysen der Hydraulik von Schlitzpässen. Wasser Wirtschaft 7/8: 67–71. (In the text)
 Nilsson C, Reidy CA, Dynesius M, Revenga C. 2005. Fragmentation and flow regulation of the world's large river systems. Science 308: 405–408. [CrossRef] [PubMed] (In the text)
 Oertel M, Schlenkhoff A. 2012. Crossbar block ramps: flow regimes, energy dissipation, friction factors, and drag forces. J Hydraul Eng 138: 440–448. [CrossRef] (In the text)
 Pagliara S, Chiavaccini P. 2006. Flow resistance of rock chutes with protruding boulders. J Hydraul Eng 132: 545–552. [CrossRef] (In the text)
 Poleni G. 1717. De motu aquae mixto libri duo, Iosephi Comini, Patavii. (In the text)
 Priestley MB. 1981. Spectral analysis and time series. London, UK: Academic Press, 890 p. (In the text)
 Puertas J, Pena L, Teijeiro T. 2004. Experimental approach to the hydraulics of vertical slot fishways. J Hydraul Eng 130: 10–23. [CrossRef] (In the text)
 R Core Team. 2016. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. (In the text)
 Rajaratnam N, Van der Vinne G, Katopodis C. 1986. Hydraulics of vertical slot fishways. J Hydraul Eng 112: 909–927. [CrossRef] (In the text)
 Rajaratnam N, Katopodis C, Mainali A. 1988. Plunging and streaming flows in pool and weir fishways. J Hydraul Eng 114: 939–944. [CrossRef] (In the text)
 Rajaratnam N, Katopodis C, Mainali A. 1989. Poolorifice and poolorificeweir fishways. Can J Civil Eng 16: 774–777. [CrossRef] (In the text)
 Rajaratnam N, Katopodis C, Solanki S. 1992. New designs for vertical slot fishways. Can J Civil Eng 19: 402–414. [CrossRef] (In the text)
 Rehbock T. 1929. Wassermessung mit scharfkantigen Überfallwehren. Z Ver dtsch Ing 73: 817–823. (In the text)
 Santos JM, Branco P, Katopodis C, Ferreira T, Pinheiro A. 2014. Retrofitting poolandweir fishways to improve passage performance of benthic fishes: effect of boulder density and fishway discharge. Ecol Eng 73: 335–344. [CrossRef] (In the text)
 SanzRonda FJ, BravoCordoba FJ, FuentesPérez JF, CastroSantos T. 2016. Ascent ability of brown trout, Salmo trutta, and two Iberian cyprinids – Iberian barbel, Luciobarbus bocagei, and northern straightmouth nase, Pseudochondrostoma duriense – in a vertical slot fishway. Knowl Manag Aquat Ecosyst 417: 10. [CrossRef] (In the text)
 Schröder RCM. 1994. Technische Hydraulik: Kompendium für den Wasserbau. Berlin, Heidelberg: SpringerVerlag, 308 p. (In the text)
 Silva AT, Santos JM, Franco AC, Ferreira MT, Pinheiro AN. 2009. Selection of Iberian barbel Barbus bocagei (Steindachner, 1864) for orifices and notches upon different hydraulic configurations in an experimental pooltype fishway. J Appl Ichthyol 25: 173–177. [CrossRef] (In the text)
 Tarrade L, Pineau G, Calluaud D, Texier A, David L, Larinier M. 2011. Detailed experimental study of hydrodynamic turbulent flows generated in vertical slot fishways. Environ Fluid Mech 11: 1–21. [CrossRef] (In the text)
 Torricelli E. 1644. Opera geometrica. Florence, Italy: Amadoro Massa & Lorenzo de Landis. (In the text)
 Towler B, Mulligan K, Haro A. 2015. Derivation and application of the energy dissipation factor in the design of fishways. Ecol Eng 83: 208–217. [CrossRef] (In the text)
 Villemonte JR. 1947. Submergedweir discharge studies. Eng News Rec 139: 866–869. (In the text)
 Wang RW, Hartlieb A. 2011. Experimental and field approach to the hydraulics of naturelike pooltype fish migration facilities. Knowl Manag Aquat Ecosyst 400: 1–18. (In the text)
 Wang RW, David L, Larinier M. 2010. Contribution of experimental fluid mechanics to the design of vertical slot fish passes. Knowl Manag Aquat Ecosyst 396: 1–21. (In the text)
 Wu S, Rajaratnam N, Katopodis C. 1999. Structure of flow in vertical slot fishway. J Hydraul Eng 125: 351–360. [CrossRef] (In the text)
 Yagci O. 2010. Hydraulic aspects of poolweir fishways as ecologically friendly water structure. Ecol Eng 36: 36–46. [CrossRef] (In the text)
Cite this article as: FuentesPérez JF, GarcíaVega A, SanzRonda FJ, Martínez de Azagra Paredes A. 2017. Villemonte's approach: a general method for modeling uniform and nonuniform performance in stepped fishways. Knowl. Manag. Aquat. Ecosyst., 418, 23.
Supplementary Material
Appendix S1. Search strategy.
Appendix S2. Field case studies.
Appendix S3. Data Appendix.
Appendix S4. Auxiliary equations.
(Access here)All Tables
Summary of the field scenarios for each fishway type (see S3 Data Appendix for detailed information about each studied scenario).
All Figures
Fig. 1 Examples of sections of stepped fishways. Each type (a) vertical slot fishway, (b) pool and weir fishway and (c) steppool naturelike fishway can have different subtypes according to their morphology and connections. 

In the text 
Fig. 2 Sketch of Villemonte's assumption for submerged flows. Q_{1} is the free flow discharge in positive direction controlled by h_{1} and Q_{2} is the counterflow discharge due to h_{2}. p = 0 in case of slots. 

In the text 
Fig. 3 Field cases location. Three types of fishways were studied in field: vertical slot fishway (VSF), pool and weir fishway (PWF) and steppool naturelike fishway (SPNF). Localization coordinates are defined in S2 Appendix. 

In the text 
Fig. 4 Examples of discharge coefficient fits for most relevant (in terms of data points) VSF subtypes. (a) Subtypes 1–3 (fit coefficients and R^{2} in Tab. 3). (b) Subtype 6 (fit coefficients and R^{2} in Tab. 3). (c) Fit of equation (4) for subtype 1 to be used in equation (2) (β_{0} = 0.509, β_{1} = −0.248, R^{2} = 0.587, RMSE = 0.0458). 

In the text 
Fig. 5 Nonuniform performance in a VSF of subtype 3, scenario observed (grey) by Rajaratnam et al. (1986) and simulated (black) by different equation combinations. M1: h_{2} = 2.712 m and Q = 0.66 m^{3}/s; M2: h_{2} = 0.931 m and Q = 0.66 m^{3}/s; uniform: h_{2} = 1.416 m and Q = 0.66 m^{3}/s. (a) Equations (1) and (4). (b) Equation (2) and constant coefficient (C = 0.56). (c) Equations (2) and (4). 

In the text 
Fig. 6 C_{s} fit for studied PWF and alternative methods proposed by specialized references. In total 180 data points were collected. 

In the text 
Fig. 7 Fits for the studied SPNF and alternative methods proposed by specialized references. In total 45 data points were collected. 

In the text 