Next Article in Journal
Development of a New Testing Approach for Decentralised Technical Sustainable Drainage Systems
Previous Article in Journal
Urban Groundwater Contamination by Non-Steroidal Anti-Inflammatory Drugs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

How Can Be Lotic Ecosystem Size More Precisely Estimated? Comparing Different Approximations in Pre-Pyrenean and Pyrenean Mountains

1
Section of Quality, Environmental Evaluation and Environment—TRAGSATEC-SEPI, C/Julian Camarillo, 6B, Sector D, 28037 Madrid, Spain
2
Department of Animal Science—Wildlife Section, University of Lleida, Av. Alcalde Rovira Roure 191, 25198 Lleida, Spain
3
Forest Science and Technology Centre of Catalonia-CTFC, 25280 Solsona, Spain
4
Iberian Society of Ichthyology, Department of Environmental Biology, University of Navarra, C/Irunlarrea 1, 31008 Pamplona/Iruña, Spain
*
Author to whom correspondence should be addressed.
Water 2021, 13(5), 721; https://doi.org/10.3390/w13050721
Submission received: 31 December 2020 / Revised: 17 February 2021 / Accepted: 1 March 2021 / Published: 6 March 2021
(This article belongs to the Section Water Quality and Contamination)

Abstract

:
Rivers are among the most biodiverse and endangered ecosystems on earth. In Europe, concern over their conservation promoted the development of legal instruments for habitat and species conservation, the Habitats Directive, and water resource management, the Water Framework Directive. This legal protection demanded the estimate of river ecosystem surface for different purposes. Different approaches allow river surface to be measured at a low cost. Some accurate techniques like satellite images or LiDAR (Light Detection and Ranging) do not always work at a large scale or for streams and small rivers. We discuss here the use of the traditional hydraulics relationship between drainage area and bankfull width as a good approach to river surface estimation. We confirm that the use of this cheap and simple method could be a good approach to estimate river surface. However, we also proved that the development of regional curves, i.e., to establish the empirical relationship based on study area data, constitutes an essential improvement to estimation.

1. Introduction

Freshwater comprises about only 2.5–2.8% of the Earth’s total volume of water, and running waters are only a small fraction (~0.0001–0.006%) of the world’s freshwater [1,2,3,4], occupying approximately 0.8% of the Earth’s surface [5]. They host a huge richness of habitats, and as a result, around 100,000–150,000 [5,6,7] out of approximately 1.81 million species [8], 5.55% of the total number [5,9], inhabit freshwater ecosystems, even though freshwater covers only around 0.8% of the Earth’s surface [7]. Accordingly, these ecosystems are considered biodiversity hotspots [7].
Rivers support some of the most biodiverse ecosystems in the world and provide essential ecosystem services to society [10], but, at the same time, freshwater ecosystems may be the most endangered ecosystems in the world [6], and rivers are between the most threatened [11,12] due to overexploitation, pollution, regulation, climate change, drainage-basin degradation, non-native species, and synergistic impacts [6,13]. As a consequence, this extensive habitat deterioration is causing a decline of biodiversity in freshwater ecosystems far greater than in the most affected terrestrial ecosystems [14].
Flow and its regime are considered as the main structuring factors of stream communities [15,16], and several studies have demonstrated a relationship between the hydrological regime and biological communities [17,18,19,20,21]. However, it is estimated that around 50% of the river volume is currently altered by either flow regulation or fragmentation, and it is expected to increase to 93% with the construction of new major hydropower dams [22]. Expanding hydropower is considered as one of the 12 pressing and emerging threats to freshwater biodiversity [11]. In Europe, most rivers are largely impacted and significantly modified by anthropogenic structures like dams, weirs, and barrages [23], negatively impacting the ecological processes and aquatic species [24]. The rivers also suffer from fragmentation by barriers to free flow [7,25], causing loss of habitat, altering the natural flow regimen, and reducing the connection between the main channel and floodplain due to flow regulation [23]. Hydrological variation caused by weirs and barrages is a substantial contributor to aquatic habitat alteration and community disturbance by altering the longitudinal continuity of lotic ecosystems, the lateral interaction with floodplains, the interchange with adjacent groundwaters, and the natural variation in their relative magnitude over time [26,27,28]. As a result, European rivers suffer from intense morphological habitat degradation, which is considered one of the most important river ecosystem alterations in Europe [29]. Pyrenean rivers are not an exception, as they are highly affected by reservoirs that impound about one-third of the mean annual runoff [30].
In this context, our societies have called for increased conservation of our ecosystems. Responding to these requirements, the Habitats Directive (HD), the main legislative work of the European Union’s nature conservation policy, and the Water Framework Directive (WFD), the main legislative work regarding water resource management, have been developed to conserve our ecosystems. Both directives, from different perspectives and approximations, have common objectives and links [31]. These legal responsibilities demand the estimation of river size and its evolution at different scales. On one hand, the Habitats Directive requires surface estimation for the different running water habitats every six years [32,33]. On the other hand, the Water Framework Directive demands the size estimation of each WFD management unit [34], and by evaluating the hydromorphology of ecosystems the connectivity and continuity of these ecosystems are indirectly measured. Therefore, a method for river size estimation for all streams and rivers at a large scale is mandatory in order to achieve the legal requirements. However, the lack of a common method and some interpretation and implementation problems have led to calls for an established protocolized method for monitoring running water habitats.
River size is a key parameter related to several river ecosystem processes such as flux exchange with the atmosphere [35,36], species richness [37], and nutrient removal [38]. Many different measures and approximations have been used to estimate river size [39], e.g., length [2,37,40], drainage area [36,41], and river surface. At a large scale, easy measures like drainage area or length have been traditionally widely used [39]; however, new techniques like satellite images can be used to measure river surface [35,42]. As of yet, these estimation methods are limited by pixel resolution, allowing measurements for rivers wider than 30 m [35,42]. Other techniques, like LiDAR and photogrammetry, are limited at a large scale by their computing processing time [43,44]. Thus, several methods have been proposed to estimate river width and river surface [35,36], but none of them seems to be widely applicable with satisfactory results and with an optimal cost–benefit ratio.
When the river ecosystem surface is estimated, it is necessary to delimitate its lateral extent, which can refer to the water channel, the bankfull, the riparian zone, or even the floodplain. Estimation of the bankfull and water channel at a large scale has been widely developed [35,36,42,45]; however, riverine and floodplain surface estimation has only been developed at a local scale [46,47]. In this study, we considered the bankfull width as the lateral extent for the river ecosystem assessment due to two main reasons: (1) Ecosystem classifications usually differ riparian and floodplain ecosystem from running water ecosystems, e.g., habitat types of the HD [48]; (2) although the river channel and floodplain could be considered as components of a single dynamic system [27], the main material and energy of most European rivers are derived through local production and the riparian zone [49]. Therefore, bankfull width is indeed a good proxy for river lateral extent.
Bankfull width has been measured by different methods in stations or small sections including field measurements [50,51], high-resolution digital elevation models, and photogrammetry [43,52]. However, these measures only have been realized in stations or small sections. At a large scale, the most widely accepted estimation for bankfull channel geometry relates river width to discharge or drainage area in a power law model [42,53]. These power laws have been largely estimated in the USA at a regional scale and for the entire country [54,55,56]. However, few studies have estimated them in Europe and only at a regional or local scale [51]. A global-scale bankfull estimation method was recently developed but only for rivers wider than 30 m [42].
Therefore, the main aim of this study was to analyze the size of freshwater lotic ecosystems from different approximations and methodologies to determine which of the existing methodologies are the most accurate and adequate by identifying the advantages and disadvantages of these approximations and, as a consequence, quantifying the size of the freshwater lotic ecosystem to achieve the mandatory requirements of the HD and WFD by providing a proved tool that can be adapted to local contexts, enabling easy implementation.

2. Materials and Methods

2.1. Study Area and Data

This study was developed for Pre-Pyrenean and Pyrenean rivers, an ecosystem type defined in the methodologies for the monitoring of the conservation of habitat types in Spain [57]. This ecosystem type is defined based on physical, chemical, and biological variables and covers rivers along approximately 390 km of the southern slopes of the Pyrenees [57] (Figure 1). They are small rivers and streams located in altitudes between 700 and 900 m with a mean discharge between 1.37 and 4.93 m3 · s−1 and a drainage area between 4.5 and 3800 km2. A total of 232 management water bodies defined according to the WFD are classified as this type of ecosystem.
Bankfull widths for 76 river sections, or reaches, of 1–2 km were collected from the national cartography of legal boundaries delimitation [58]. This digital maps establish the bankfull channel limits based on historical and actual orthophotos, digital elevation models, and bankfull discharge [59]. However, they have only been developed for a few rivers. For each of the 76 water bodies where these guidelines have been established, we randomly selected a point and created a buffer of 1 km to delimitate and standardize a river section. For every section, we obtained the bankfull limits via the national guidelines of legal river boundaries (Figure 1).

2.2. Methods: Size Estimation

The surface of each polygon delimitated by the bankfull channel banks was calculated (Figure 1). This surface was considered as the real bankfull channel surface of each section (Sb). We extracted the centerline of each section using a skeletonization procedure of the bankfull channel bank lines. The length (L) of each centerline was calculated as an approximation of the river ecosystem size, a classical and probably the most widely used approximation of river size. The bankfull width of each section (Wb) was calculated as the median of the measurements of Wb every 50 m within each reach. The drainage area of each section (A) was calculated as the mean value of all sections. We used the drainage area raster developed by the Spanish government for its estimation [60]. All these operations were conducted using ArcMap 10.5 [61].
At a large scale, the most widely accepted estimation method for bankfull channel geometry relates river width (W) to discharge (Q) or drainage area (A) in a power law model (Equations (1) and (2)) with two empirical factors, α and β [51,54,55,56,62,63,64,65,66,67,68,69,70,71,72,73]. All the power laws found in the bibliography were collected and published (see Data Availability Statement) to analyze the difference between them and the regional power law of our study area.
W = α · Q β
W = α · A β
First, we estimated the power law that relates bankfull width and drainage area for our study area. We used the median value of the bankfull width (Wb) and the mean drainage area (A) of each section to estimate the α’ and β’ coefficients. To develop the power law for our study, we log-transformed both variables and estimated α’ and β’ through simple linear regression, i.e., Y = β0 + β1 X, where Y = log10 (Wb), X = log10 (A), and β0 and β1 are regression parameters (intercept and slope, respectively). Thus, the coefficient and the exponent of the power law relation (Equation (2)) are α’ = 10β0 and β’ = β1, respectively.
To explore the relationship between bankfull channel width and drainage area, we used multiple linear regression models in combination with mean annual precipitation, channel slope (primary controls on discharge and channel morphology), and elevation as explanatory variables of channel width [56]. Mean annual precipitation in mm (MAP) data were collected from WordlCilm2 [74]. Mean altitude (MA) in m and mean channel slope (MS) were obtained using a digital elevation model of 5 m pixels. All variables, except elevation, were log-transformed, yielding a multiplicative effects model of the form [58]. We ranked the models using the corrected Akaike information criterion (AICc) and selected the model with the lowest AICc (∆AICc < 2).
To estimate the bankfull channel surface for each water body, the average bankfull width (Wbe) was calculated for each section using our regional power law, considering the mean drainage area of the section (A). Then, we estimated the bankfull surface (Sbe) for a water body as the product of the centerline length (L) and Wbe, considering the section as a parallelogram. This operation was repeated for all the power laws found in the bibliography to estimate bankfull surface (Wbe) for each one.

2.3. Methods: Statistical Analysis

We used multiple linear regression to explore the differences between the real bankfull surface (Sb) of each section and the estimated surface (Sbe) according to our regional power law. We included the length of the centerline (L) and median width (Wb) of the section, as well as the sinuosity index (SI) [2] (primary controls on size, width and length, and channel complexity, sinuosity index). We ranked the models using the corrected Akaike information criterion (AICc) and selected the model with the lowest AICc (∆AICc < 2). All data and results are shown in Supplementary Materials 1.
Finally, we compared the differences between the real bankfull surface (Sb) and the estimated surfaces (Sbe) via linear regressions. The Sb and Sbe values were previously log-transformed to achieve residual normality, linearity, and homoscedasticity. The resultant linear models had the expression log10 (Sb) = m + n log10 (Sbe), where Sb and Sbe are the real and the estimated bankfull surface, and m and n are the model parameters (intercept and slope, respectively). If we remove the logarithms, the resulting expression is Sb = 10m Sben. Thus, when m is nearer to 0, and n to 1, the estimated values of Sb are more similar to the real Sbe.
All statistical analyses were carried out by using R version 4.0.2 [75].

3. Results

The relationship between the logarithms of Wb and A showed a significant increase in log10 (Wb) with the increase of log10 (A) (Figure 2) (F1,74 = 62.89, p < 0.001, r2adjusted = 0.45). With the linear model coefficients, we calculated α and β values of Equation (2) for our regional data. The regional power law relationship is as follows:
Wbe = 5.531 · A 0.362
The multiple linear regression including in the relationship log10 (MAP), log10 (MS), and MA showed an improvement in the relation between log10 (Wb) and log10 (A) when log10 (MS) was included in the model (Table 1), although the model with only log10 (A) as an explanatory variable was roughly equivalent to the best model, and it was the only significant explanatory variable for all models (Table 1).
The multiple linear regressions developed to assess the differences between the estimated bankfull channel surface (Sbe) and real bankfull surface (Sb) showed that size and river complexity influence the estimation error. As shown in Table 2, the most parsimonious model only included Wb and SI (F2,73 = 94.35, p < 0.001, r2adjusted = 0.72). The difference between Sb and Sbe was positively correlated with Wb and SI. The surface was overestimated (Sbe > Sb or Sb − Sbe < 0) in sections with small Wb and underestimated (Sbe < Sb or Sb − Sbe > 0) where Wb was high (Figure 3a). SI showed a slightly positive correlation with the estimation error. The surface of complex rivers (SI > 1.5) was more frequently overestimated than the surface of less complex rivers (SI ≈ 1) (Figure 3b).
Lastly, we compared the results of the bankfull channel surface estimation obtained by our regional curve with the results obtained by the other power laws. Figure 4 shows the regional curve as one of the most similar to the real values. Approximately 50% of the linear models between log10 (Sb) and log10 (Sbe) had a higher r2adjusted than the model in which Sbe was estimated by our regional curve (F1,74 = 51.24, p < 0.001, r2adjusted = 0.40) (the highest r2adjusted is 0.41). However, the mean differences between Sb and Sbe estimated by the regional curve were the third-smallest.

4. Discussion

As expected, and confirming previous similar studies, the drainage area is intensively related to width, and our results have shown that the estimated bankfull channel surface is directly related to the real bankfull surface, and this estimation significantly improves when the mean channel slope is included in the model. Moreover, the development of a regional power law based on real data improved the estimation. However, the surface was overestimated and underestimated due to river size and complexity. On one hand, upstream sections were normally overestimated, whereas downstream sections were underestimated (Figure 3a). On the other hand, sections highly sinuous were more frequently underestimated (Figure 3a).
We demonstrate the utility of power laws relating drainage area and bankfull width as a fast assessment of river ecosystem size for streams and small rivers at a large scale, confirming this method as a tool of easy implementation. Other studies indicated the importance of developing power law relationships between bankfull channel width and drainage areas for preliminary channel design for restoration projects [73] or modeling habitat availability and suitability for fish [76]. According to our results, these estimations are robust and can be consolidated as an optimal approximation of river size with a low cost–benefit ratio, at least until LiDAR and remote sensing can be optimized for these purposes. We found that the best obtained results corresponded to the power law established with our data within the regional context of the Pre-Pyrenean and Pyrenean range. Thanks to these results, we can state that using regional power law models is the most optimal approximation, and in case comparing with other regions or models is required, models can be easily compared by recalculating the values. As outlined, using a regional power law provides the best results. However, in case it is not possible to use a regional power law due to the lack of digital terrain models or real data, it is also possible to use a power law model from similar ecoregions.
Our reported value for α could be considered as high compared with the values found in the bibliography, the 13th largest α detected/historically detected, whereas the β value is near the median value of all β values registered, 0.36 (see Data Availability Statement). High values of α, such as those estimated in our study area, are related to high pluviometry [56]. The mean pluviometry for all ecosystem sizes is 1273.84 ± 213.84 [57], a value that could be considered as relatively high and could be related to the fact that the Pyrenees is a mountain range that links the Atlantic Ocean with the Mediterranean Sea with a clear west–east gradient of humidity and pluviometry. The β value is lower the less the width grows with the watershed area, but this relationship can be affected by the degree of alteration of the river basin. It has been proved that human impacts can reduce β values [56,69,70]. On the other hand, high β values have also been related to climatic conditions, e.g., arid zones [56]. In our case, the high estimated β value is coherent with other studies [56,70] because most of our studied rivers can be considered as altered by dams, hydropower stations, etc. [30]. Moreover, the hydrological alteration and flow regulation caused by these perturbances can influence the estimated regional power law, increasing the β value [70].
To improve our model, we included in the power law other factors that are related to discharge (mean annual precipitation) and channel morphology (mean channel slope and mean altitude). We found only slight improvement when we included mean channel slope in the model (∆AICc < 2). This result is coherent since the chosen ecosystem type was defined partly based on physical variables such as discharge, altitude, or slope [57]. Therefore, it is logical that the relationship was not significantly improved based on these variables because they are intrinsically included.
As previously mentioned, the regional power law provides the most optimal results. However, this is an imperfect approximation, and the comparison between the real and the estimated bankfull channel area (Figure 3) shows an underestimation in wider sections that are located in downstream reaches and an overestimation in small lotic ecosystems (Figure 3). Degradation of the riverbank often causes a widening processes [77], which could explain the observed underestimation, i.e., the downstream reaches are wider than expected. Most of the areas surrounding rivers in the low parts of our study area are occupied by cultures, whereas the high areas are mostly forest in a mountainous area where rivers flow through narrow valleys where widening is not possible and overestimation of Wbe is explained. According to these results, altitudinal gradient and orography could explain these over- and underestimation results. By including mean channel slope (MS) in our model, we obtained a good model (ER = 1) for explaining river width (Table 2). This is coherent with the mentioned under- and overestimation problems because upstream rivers have a higher MS than lower reaches. The surface was slightly underestimated in more complex sections with higher SI. Complex rivers had a higher heterogeneity of width [78], which was not captured in the estimations producing underestimations.
As is shown in Figure 4, the variation in the empirical values of the power laws triggered a huge variability in surface estimation. This variation has been reported in several studies developed in the USA [54,55,56], and it has been demonstrated that different factors influence this relationship, e.g., human alteration [69,70], riparian vegetation [69], riverbed material [56], and climate and topographic factors [56]. Therefore, it is important to develop a power law at a regional scale for river size estimation when possible. However, all estimated bankfull channel surfaces are significantly related to the real surface and could be useful for comparison of the river size.
These results confirm that regional power law models are superior to other models that have been developed for specific hydrogeological contexts. Of course, even in a regional context, we can expect variability and errors in the estimation as we have observed. For this reason, the use of regional power law models is conditioned by the characteristics of the region; therefore, these regions must be internally homogenous with similar climate, geology, pluviometry, and, of course, orography. However, this fast approach to river size might be replaced in the future by more accurate new techniques such as LiDAR or satellite images that currently present limitations [44]. Indeed, recent studies have demonstrated the precision of satellite images to estimate bankfull and water channel at a large scale for rivers wider than 30 m [35,42] or accurate measurements derived from LiDAR [79]. However, these methods still present limitations that restrain their use at a large scale for narrower rivers or sections. As a consequence, using regional power laws could be considered as an optimal temporary solution until these more accurate techniques can be applied to most of the lotic ecosystems. With the application of this method, we can address the lack of surface data from lotic ecosystems since the Habitats Directive was approved 30 years ago. Classical methods like river length also should be calculated because of historical reasons and due to the rivers mainly comprising a lineal ecosystem where the most important processes and patterns take place in longitudinal axes [80,81,82].

Supplementary Materials

Supplementary material is available at https://doi.org/10.6084/m9.figshare.14169530.

Author Contributions

Conceptualization, F.C.S., and J.R.S.-G.; methodology, F.C.S. and J.R.S.-G.; formal analysis, F.C.S. and J.R.S.-G.; investigation, F.C.S. and J.R.S.-G.; resources, F.C.S.; data curation, F.C.S.; writing—original draft preparation, F.C.S. and J.R.S.-G.; writing—review and editing, F.C.S., F.C. and J.R.S.-G.; su-pervision, F.C. and J.R.S.-G.; project administration, F.C.; funding acquisition, F.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Ecological Transition and Demographic Challenge, project number 3081147. This research was supported by the Government of Catalonia (CERCA Program).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data used in this study is available at https://doi.org/10.6084/m9.figshare.13490208.

Acknowledgments

This manuscript would not be possible without the support and the motivation of Rafael Hidalgo, Olga Lamas Murúa, and Elena Bermejo Bermejo. This manuscript was born within the context of the Project “Establishment of a national monitoring system of the Conservation Status of the Habitat of Community Interest” 3081147 and supported by the Government of Spain and the Ministry for the Ecological Transition and Demographic Challenge.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Downes, B.J.; Barmuta, L.A.; Fairweather, P.G.; Faith, D.P.; Keough, M.J.; Lake, P.S.; Mapstone, B.D.; Quinn, G.P. Monitoring Ecological Impacts. Concepts and Practice in Flowing Waters; Cambridge University Press: Cambridge, UK, 2004; ISBN 978-0-521-06529-0. [Google Scholar]
  2. Allan, J.D. Stream Ecology. Structure and Function of Running Waters, 1st ed.; Chapman and Hall: London, UK, 1995; ISBN 978-0-412-35530-1. [Google Scholar]
  3. Shiklomanov, I.A. World Freshwater Resources. In Water in Crisis: A Guide to World’s Freshwater Resources; Gleick, P.H., Ed.; Oxford University Press Inc.: New York, NY, USA, 1993; pp. 13–24. ISBN 0-19-507628-1. [Google Scholar]
  4. Malmqvist, B.; Rundle, S. Threats to the running water ecosystems of the world. Environ. Conserv. 2002, 29, 134–153. [Google Scholar] [CrossRef]
  5. Gleick, P.H. Water Resources. In Encyclopedia of Climate and Weather; Shneider, S.H., Ed.; Oxford University Press: New York, NY, USA, 1996; pp. 817–823. [Google Scholar]
  6. Dudgeon, D.; Arthington, A.H.; Gessner, M.O.; Kawabata, Z.-I.; Knowler, D.J.; Lévêque, C.; Naiman, R.J.; Prieur-Richard, A.-H.; Soto, D.; Stiassny, M.L.J.; et al. Freshwater biodiversity: Importance, threats, status and conservation challenges. Biol. Rev. Camb. Philos. Soc. 2006, 81, 163–182. [Google Scholar] [CrossRef]
  7. Strayer, D.L.; Dudgeon, D. Freshwater biodiversity conservation: Recent progress and future challenges. J. North Am. Benthol. Soc. 2010, 29, 344–358. [Google Scholar] [CrossRef] [Green Version]
  8. Hawksworth, D.J.; Kalin-Arroyo, M.T. Magnitude and Distribution of Biodiversity. In Global Biodiversity Assessment; Heywood, V.H., Ed.; Cambridge University Press: Cambridge, UK, 1995; pp. 107–191. [Google Scholar]
  9. Species 2000 & ITIS Catalogue of Life, 2019 Annual Checklist. Available online: www.catalogueoflife.org/annual-checklist/2019 (accessed on 19 December 2020).
  10. Grizzetti, B.; Liquete, C.; Pistocchi, A.; Vigiak, O.; Zulian, G.; Bouraoui, F.; De Roo, A.; Cardoso, A.C. Relationship between ecological condition and ecosystem services in European rivers, lakes and coastal waters. Sci. Total Environ. 2019, 671, 452–465. [Google Scholar] [CrossRef] [PubMed]
  11. Reid, A.J.; Carlson, A.K.; Creed, I.F.; Eliason, E.J.; Gell, P.A.; Johnson, P.T.J.; Kidd, K.A.; Maccormack, T.J.; Olden, J.D.; Ormerod, S.J.; et al. Emerging threats and persistent conservation challenges for freshwater biodiversity. Biol. Rev. 2019, 94, 849–873. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Grzybowski, M.; Glińska-Lewczuk, K. Principal threats to the conservation of freshwater habitats in the continental biogeographical region of Central Europe. Biodivers. Conserv. 2019, 28, 4065–4097. [Google Scholar] [CrossRef] [Green Version]
  13. Dudgeon, D. Anthropocene Extinction: Global Threats to Riverine Biodiversity and the Tragedy of the Freshwater Commons. In River Conservation. Challenges and Opportunities; Sabater, S., Elosegi, A., Eds.; Bilbao: Bloomington, IN, USA, 2013; pp. 129–167. [Google Scholar]
  14. Sala, O.E.; Chapin, F.S.; Armesto, J.J.; Berlow, E.; Bloomfield, J.; Dirzo, R.; Huber-Sanwald, E.; Huenneke, L.F.; Jackson, R.B.; Kinzig, A.; et al. Global biodiversity scenarios for the year 2100. Science 2000, 287, 1770–1774. [Google Scholar] [CrossRef]
  15. Poff, N.L.; Allan, J.D.; Bain, M.B.; Karr, J.R.; Prestegaard, K.L.; Richter, B.D.; Sparks, R.E.; Stromberg, J.C. The natural flow regime. Bioscience 1997, 47, 769–784. [Google Scholar] [CrossRef]
  16. Poff, N.L. Beyond the natural flow regime? Broadening the hydro-ecological foundation to meet environmental flows challenges in a non-stationary world. Freshw. Biol. 2018, 63, 1011–1021. [Google Scholar] [CrossRef]
  17. Bunn, S.E.; Arthington, A.H. Basic Principles and Ecological Consequences of Altered Flow Regimes for Aquatic Biodiversity. Environ. Manag. 2002, 30, 492–507. [Google Scholar] [CrossRef] [Green Version]
  18. Monk, W.A.; Wood, P.J.; Hannah, D.M.; Wilson, D.A.; Extence, C.A.; Chadd, R.P. Flow variability and macroinvertebrate community response within riverine systems. River Res. Appl. 2006, 22, 595–615. [Google Scholar] [CrossRef] [Green Version]
  19. Belmar, O.; Velasco, J.; Gutiérrez-Cánovas, C.; Mellado-Díaz, A.; Millán, A.; Wood, P.J. The influence of natural flow regimes on macroinvertebrate assemblages in a semiarid Mediterranean basin. Ecohydrology 2013, 6, 363–379. [Google Scholar] [CrossRef] [Green Version]
  20. Solans, M.A.; Mellado-Díaz, A. A Landscape-Based Regionalization of Natural Flow Regimes in the Ebro River Basin and Its Biological Validation. River Res. Appl. 2015, 31, 457–469. [Google Scholar] [CrossRef]
  21. Mellado-Díaz, A.; Sánchez-González, J.R.; Guareschi, S.; Magdaleno, F.; Toro Velasco, M. Exploring longitudinal trends and recovery gradients in macroinvertebrate communities and biomonitoring tools along regulated rivers. Sci. Total Environ. 2019, 695, 133774. [Google Scholar] [CrossRef] [PubMed]
  22. Barbarossa, V.; Schmitt, R.J.P.; Huijbregts, M.A.J.; Zarfl, C.; King, H.; Schipper, A.M. Impacts of current and future large dams on the geographic range connectivity of freshwater fish worldwide. Proc. Natl. Acad. Sci. USA 2020, 117, 3648–3655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Tockner, K.; Uehlinger, U.; Robinson, T.R. Rivers of Europe, 1st ed.; Elsevier: London, UK, 2009; ISBN 978-0-12-369449-2. [Google Scholar]
  24. Gough, P.; Fernandez Garrido, P.; Van Herk, J. Dam Removal. A Viable Solution for the Future of our European Rivers. Available online: https://damremoval.eu/wp-content/uploads/2018/07/Dam-Removal-Europe-Report-2018-DEF-1.pdf (accessed on 19 December 2020).
  25. Grill, G.; Lehner, B.; Thieme, M.; Geenen, B.; Tickner, D.; Antonelli, F.; Babu, S.; Borrelli, P.; Cheng, L.; Crochetiere, H.; et al. Mapping the world’s free-flowing rivers. Nature 2019, 569, 215–221. [Google Scholar] [CrossRef] [PubMed]
  26. Poff, N.L. A river that flows free connects up in 4D. Nature 2019, 569, 201–202. [Google Scholar] [CrossRef] [Green Version]
  27. Junk, W.J.; Bayley, P.B.; Sparks, R.E. The flood pulse concept in river & floodplain systems. Can. Spec. Pub. Aquat. Sci. 1989, 106, 110–127. [Google Scholar]
  28. Tockner, K.; Malard, F.; Ward, J.V. An extension of the Flood Pulse Concept. Hydrol. Process. 2000, 14, 2861–2883. [Google Scholar] [CrossRef]
  29. Belletti, B.; Garcia de Leaniz, C.; Jones, J.; Bizzi, S.; Börger, L.; Segura, G.; Castelletti, A.; van de Bund, W.; Aarestrup, K.; Barry, J.; et al. More than one million barriers fragment Europe’s rivers. Nature 2020, 588, 436–441. [Google Scholar] [CrossRef]
  30. Batalla, R.J.; Gómez, C.M.; Kondolf, G.M. Reservoir-induced hydrological changes in the Ebro River basin (NE Spain). J. Hydrol. 2004, 290, 117–136. [Google Scholar] [CrossRef]
  31. DG Environment. Links between the Water Framework Directive (WFD 2000/60/EC) and Nature Directives (Birds Directive 2009/147/EC and Habitats Directive 92/43/EEC) Frequently Asked Questions; EC Europa: Brussels, Belgium, 2011. [Google Scholar]
  32. European Community. Council Directive 92/43/EEC of 21 May 1992 on the Conservation of Natural Habitats and of Wild Fauna and Flora; European Community: Brussels, Belgium, 1992; pp. 568–583. [Google Scholar]
  33. European Commission. Reporting under Article 17 of the Habitats Directive Explanatory Notes and Guidelines for the Period 2013–2018; European Commision: Brussels, Belgium, 2017. [Google Scholar]
  34. European Commission. Directive 2000/60/EC oftThe European Parlamient And of the Council of 23 October 2000 Establishing a Framework for Community Action in the Field of Water Policy; European Commision: Brussels, Belgium, 2000; pp. 1–71. [Google Scholar]
  35. Allen, G.H.; Pavelsky, T. Global extent of rivers and streams. Science 2018, 361, 585–588. [Google Scholar] [CrossRef] [Green Version]
  36. Downing, J.A.; Cole, J.J.; Duarte, C.M.; Middelburg, J.J.; Melack, J.M.; Prairie, Y.T.; Kortelainen, P.; Striegl, R.G.; McDowell, W.H.; Tranvik, L.J. Global abundance and size distribution of streams and rivers. Inl. Waters 2012, 2, 229–236. [Google Scholar] [CrossRef]
  37. Oberdorff, T.; Guilbert, E.; Lucchetta, J.C. Patterns of fish species richness in the Seine River basin, France. Hydrobiologia 1993, 259, 157–167. [Google Scholar] [CrossRef]
  38. Wollheim, W.M.; Vörösmarty, C.J.; Peterson, B.J.; Seitzinger, S.P.; Hopkinson, C.S. Relationship between river size and nutrient removal. Geophys. Res. Lett. 2006, 33, 2–5. [Google Scholar] [CrossRef] [Green Version]
  39. Dynesius, M.; Nilsson, C. Fragmentation and flow regulation of river systems in the northern third of the world. Science 1994, 266, 753–762. [Google Scholar] [CrossRef] [PubMed]
  40. Hauer, F.R.; Lamberti, G.A. Methods in Stream Ecology, 2nd ed.; Hauer, F.R., Lamberti, G.A., Eds.; Academic Press: San Diego, CA, USA, 2007; ISBN 9780123329080. [Google Scholar]
  41. Leopold, L.B.; Maddock, T.J. The Hydraulic Geometry of Stream Channels and Some Physiographic Implications; US Government Printing Office: Washington, DC, USA, 1953; Volume 252.
  42. Lin, P.; Pan, M.; Allen, G.H.; de Frasson, R.P.; Zeng, Z.; Yamazaki, D.; Wood, E.F. Global Estimates of Reach—Level Bankfull River Width Leveraging Big Data Geospatial Analysis. Geophys. Res. Lett. 2020, 47, 1–12. [Google Scholar] [CrossRef]
  43. Demarchi, L.; Bizzi, S.; Piégay, H. Hierarchical object-based mapping of riverscape units and in-stream mesohabitats using lidar and VHR imagery. Remote Sens. 2016, 8, 97. [Google Scholar] [CrossRef] [Green Version]
  44. Tomsett, C.; Leyland, J. Remote sensing of river corridors: A review of current trends and future directions. River Res. Appl. 2019, 35, 779–803. [Google Scholar] [CrossRef]
  45. Andreadis, K.M.; Schumann, G.J.P.; Pavelsky, T. A simple global river bankfull width and depth database. Water Resour. Res. 2013, 49, 7164–7168. [Google Scholar] [CrossRef]
  46. Thoms, M.C. Floodplain–river ecosystems: Lateral connections and the implications of human interference. Geomorphology 2003, 56, 335–349. [Google Scholar] [CrossRef]
  47. Scown, M.W.; Thoms, M.C.; De Jager, N.R. An index of floodplain surface complexity. Hydrol. Earth Syst. Sci. 2016, 20, 431–441. [Google Scholar] [CrossRef] [Green Version]
  48. Evans, D. The habitats of the European union habitats directive. Biol. Environ. 2006, 106, 167–173. [Google Scholar] [CrossRef] [Green Version]
  49. Thorp, J.H.; Delong, M.D. Dominance of autochthonous autotrophic carbon in food webs of heterotrophic rivers. Oikos 2002, 96, 543–550. [Google Scholar] [CrossRef]
  50. US Environmental Protection Agency. Wadeable Streams Assessment: A Collaborative Survey of the Nation’s Streams; US Environmental Protection Agency: Washington, DC, USA, 2006; p. 113.
  51. Vianello, A.; D’Agostino, V. Bankfull width and morphological units in an alpine stream of the dolomites (Northern Italy). Geomorphology 2007, 83, 266–281. [Google Scholar] [CrossRef]
  52. Mount, N.J.; Louis, J.; Teeuw, R.M.; Zukowskyj, P.M.; Stott, T. Estimation of error in bankfull width comparisons from temporally sequenced raw and corrected aerial photographs. Geomorphology 2003, 56, 65–77. [Google Scholar] [CrossRef]
  53. Gleason, C.J. Hydraulic geometry of natural rivers: A review and future directions. Prog. Phys. Geogr. 2015, 39, 337–360. [Google Scholar] [CrossRef]
  54. Blackburn-Lynch, W.; Agouridis, C.T.; Barton, C.D. Development of Regional Curves for Hydrologic Landscape Regions (HLR) in the Contiguous United States. J. Am. Water Resour. Assoc. 2017, 53, 903–928. [Google Scholar] [CrossRef]
  55. Bieger, K.; Rathjens, H.; Allen, P.M.; Arnold, J.G. Development and Evaluation of Bankfull Hydraulic Geometry Relationships for the Physiographic Regions of the United States. J. Am. Water Resour. Assoc. 2015, 51, 842–858. [Google Scholar] [CrossRef] [Green Version]
  56. Faustini, J.M.; Kaufmann, P.R.; Herlihy, A.T. Downstream variation in bankfull width of wadeable streams across the conterminous United States. Geomorphology 2009, 108, 292–311. [Google Scholar] [CrossRef] [Green Version]
  57. Sánchez-González, J.-R.; Mellado-Díaz, A. Descripción de un Procedimiento para Generar una Tipología de Hábitats Lóticos Existentes en España; Ministerio para la Transición Ecológica: Madrid, Spain, 2019.
  58. MITERD. Cauces con DPH Cartográfico. Available online: https://www.miteco.gob.es/es/cartografia-y-sig/ide/descargas/agua/dph-cartografico.aspx (accessed on 1 August 2020).
  59. MARM. Guía Metodológica para el Desarrollo del Sistema Nacional de Cartografía de Zonas Inundables; Ministerio para la Transición Ecológica y el Reto Demográfico: Madrid, Spain, 2011.
  60. MITERD. Superficie de Cuenca Vertiente a Cada Celda de 25 × 25 m. Available online: https://www.miteco.gob.es/es/cartografia-y-sig/ide/descargas/agua/superficie-cuenca-25x25.aspx (accessed on 1 August 2020).
  61. Environmental System Research Institute (ESRI). ArcGIS Desktop: Release 10.5; Environmental Systems Research Institute: Redlands, CA, USA, 2016. [Google Scholar]
  62. Moody, J.A.; Meade, R.H.; Jones, D.R. Lewis and Clarke’s Observations and Measurements of Geomorphology and Hidrology, and Changes with Time; U.S. Geological Survey Circular 1246: Reston, VA, USA, 2003. [Google Scholar]
  63. Castro, J.M.; Jackson, P.L. Bankfull discharge recurrence interval and regional hidraulic geometric relationship: Patterns in the pacific northwest, USA. J. Am. Water Resour. Assoc. 2001, 37, 1249–1262. [Google Scholar] [CrossRef]
  64. Metcalf, C.K.; Wilkerson, S.D.; Harman, W.A. Bankfull regional curves for north and Northwest Florida streams. J. Am. Water Resour. Assoc. 2009, 45, 1260–1272. [Google Scholar] [CrossRef]
  65. Krstolic, J.L.; Chaplin, J.J. Bankfull Regional Curves for Streams in the Non-Urban, Non-Tidal Coastal Plain Physiographic Province, Virginia and Maryland Scientific Investigations Report 2007–5162; US Geological Survey: Washington, DC, USA, 2007.
  66. Brockman, R.R.; Agouridis, C.T.; Workman, S.R.; Ormsbee, L.E.; Fogle, A.W. Bankfull regional curves for the inner and the outer bluegrass regions of kentucky. J. Am. Water Resour. Assoc. 2012, 48, 391–406. [Google Scholar] [CrossRef]
  67. Mohamoud, Y.M.; Parmar, R.S. Estimating streamflow and associated hydraulic geometry, the mid-Atlantic region, USA. J. Am. Water Resour. Assoc. 2006, 5017, 755–768. [Google Scholar] [CrossRef]
  68. Lee, W.H.; Choi, H.S. Characteristics of Bankfull Discharge and its Estimation using Hydraulic Geometry in the Han River Basin. J. Civ. Eng. 2018, 22, 2290–2299. [Google Scholar] [CrossRef]
  69. Hession, W.C.; Pizzuto, J.E.; Johnson, T.E.; Horwitz, R.J. Influence of bank vegetation on channel morphology in rural and urban watersheds. Geology 2003, 31, 147–150. [Google Scholar] [CrossRef]
  70. Doll, B.A.; Wise-Frederick, D.E.; Buckner, C.M.; Wilkerson, S.D.; Harman, W.A.; Smith, R.E.; Spooner, J. Hydraulic geometry relationships for urban streams throughout the piedmont of North Carolina. J. Am. Water Resour. Assoc. 2002, 38, 641–651. [Google Scholar] [CrossRef]
  71. Sweet, W.V.; Geratz, J.W. Bankfull hydraulic geometry relationship and recurrence intervals for North Carolina’s coasta plain 1. JAWRA 2003, 39, 861–871. [Google Scholar]
  72. Golden, L.A.; Springer, G.S. Channel geometry, median grain size, and stream power in small mountain streams. Geomorphology 2006, 78, 64–76. [Google Scholar] [CrossRef]
  73. Johnson, P.A.; Fecko, B.J. Regional channel geometry equations: A statiscial comparison for physiographic provinces in the Eastern US. River Res. Appl. 2008, 24, 823–834. [Google Scholar] [CrossRef]
  74. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  75. R Development Core Team. R: A Language and Environment for Statistical Computing. R version 4.0.2 “Taking Off Again”; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  76. Rosenfeld, J.S.; Post, J.; Robins, G.; Hatfield, T. Hydraulic geometry as a physical template for the River Continuum: Application to optimal flows and longitudinal trends in salmonid habitat. Can. J. Fish. Aquat. Sci. 2007, 64, 755–767. [Google Scholar] [CrossRef] [Green Version]
  77. ASCE Task Committee on Hydraulics, Bank Mechanics, and Model of River Width Adjust. River Width Adjustment. Processes and Mechanisms. J. Hydraul. Eng. 1998, 124, 881–902. [Google Scholar] [CrossRef]
  78. Eke, E.C.; Czapiga, M.J.; Viparelli, E.; Shimizu, Y.; Imran, J. Coevolution of width and sinuosity in meandering rivers. J. Fluid Mech. 2014, 127–174. [Google Scholar] [CrossRef]
  79. Demarchi, L.; Bizzi, S.; Piégay, H. Regional hydromorphological characterization with continuous and automated remote sensing analysis based on VHR imagery and low-resolution LiDAR data. Earth Surf. Process. Landf. 2017, 42, 531–551. [Google Scholar] [CrossRef]
  80. Allan, J.D.; Castillo, M.M. Stream Ecology. Structure and Function of Running Waters; Springer: Dordrecht, The Netherlands, 2007; p. 436. ISBN 9781402055829. [Google Scholar]
  81. Vannote, R.L.; Minshall, G.W.; Cummins, K.W.; Sedell, J.R.; Cushing, C.E. The River Continuum Concept. Can. J. Fish. Aquat. Sci. 1980, 37, 130–137. [Google Scholar] [CrossRef]
  82. Newbold, J.D.; Elwood, J.W.; O’Neill, R.V.; Van Winkle, W. Measuring Nutrient Spiralling in Streams. Can. J. Fish. Aquat. Sci. 1981, 38, 860–863. [Google Scholar] [CrossRef]
Figure 1. Study area. (a) Location of the analyzed rivers, Pre-Pyrenean and Pyrenean, in Europe. (b,c) Examples of sections, in red, of legal bankfull river boundaries according to national cartography [58]. (d) Locations of all Pre-Pyrenean and Pyrenean rivers are marked with dark blue, and the locations of all sections in red. Coordinates in meters, EPSG: 25831.
Figure 1. Study area. (a) Location of the analyzed rivers, Pre-Pyrenean and Pyrenean, in Europe. (b,c) Examples of sections, in red, of legal bankfull river boundaries according to national cartography [58]. (d) Locations of all Pre-Pyrenean and Pyrenean rivers are marked with dark blue, and the locations of all sections in red. Coordinates in meters, EPSG: 25831.
Water 13 00721 g001
Figure 2. Relationship between the logarithm of the drainage area (A) in km2 and the logarithm of the real bankfull width (Wb) in m at the 76 studied sections. The blue line shows linear regression (F1,74 = 62.73, p < 0.005, r2adjusted = 0.46) with 95% CI (grey).
Figure 2. Relationship between the logarithm of the drainage area (A) in km2 and the logarithm of the real bankfull width (Wb) in m at the 76 studied sections. The blue line shows linear regression (F1,74 = 62.73, p < 0.005, r2adjusted = 0.46) with 95% CI (grey).
Water 13 00721 g002
Figure 3. (a) Relationship between the differences between Sb and Sbe in m2 and the real bankfull width (Wb) in m at the 76 studied sections. The blue line shows linear regression (F1,74 = 158.4, p < 0.005, r2adjusted = 0.68) with 95% CI (grey). (b) Relationship between the differences between Sb and Sbe in m2 and the sinuosity index (SI) in m at the 76 studied sections. The blue line shows linear regression (F1,74 = 1.52, p > 0.05, r2adjusted = 0.01) with 95% CI (grey).
Figure 3. (a) Relationship between the differences between Sb and Sbe in m2 and the real bankfull width (Wb) in m at the 76 studied sections. The blue line shows linear regression (F1,74 = 158.4, p < 0.005, r2adjusted = 0.68) with 95% CI (grey). (b) Relationship between the differences between Sb and Sbe in m2 and the sinuosity index (SI) in m at the 76 studied sections. The blue line shows linear regression (F1,74 = 1.52, p > 0.05, r2adjusted = 0.01) with 95% CI (grey).
Water 13 00721 g003
Figure 4. Relationships between log10 (Sb) and log10 (Sbe) in m2 at the 76 studied sections for the Sbe estimated by the different regional curves. The black line shows the linear regression where Sbe was estimated by our regional curve (F1,74 = 51.24, p < 0.001, r2adjusted = 0.40). The red line shows the linear regression considering the Sb values for Sbe. Grey lines show all the different relationships between Sb and the Sbe estimated for the different curves.
Figure 4. Relationships between log10 (Sb) and log10 (Sbe) in m2 at the 76 studied sections for the Sbe estimated by the different regional curves. The black line shows the linear regression where Sbe was estimated by our regional curve (F1,74 = 51.24, p < 0.001, r2adjusted = 0.40). The red line shows the linear regression considering the Sb values for Sbe. Grey lines show all the different relationships between Sb and the Sbe estimated for the different curves.
Water 13 00721 g004
Table 1. Comparison of Akaike information criterion (AICc) deltas and weights for linear models testing the effects of drainage area (A), mean channel slope (MS), mean annual precipitation (MP), and mean altitude (MA) and response (bankfull Width, Wb). Shown are the number of parameters (K), the AICc, the ΔAICc, AICc weight (wAICc), and the evidence ratios (ERs). The best model and models roughly equivalent to the best (ΔAICc ≈ 2 or lower) are highlighted in bold. The only significant explanatory variable in all models is log10 (A) (p < 0.005). The sample size is 76 for all models.
Table 1. Comparison of Akaike information criterion (AICc) deltas and weights for linear models testing the effects of drainage area (A), mean channel slope (MS), mean annual precipitation (MP), and mean altitude (MA) and response (bankfull Width, Wb). Shown are the number of parameters (K), the AICc, the ΔAICc, AICc weight (wAICc), and the evidence ratios (ERs). The best model and models roughly equivalent to the best (ΔAICc ≈ 2 or lower) are highlighted in bold. The only significant explanatory variable in all models is log10 (A) (p < 0.005). The sample size is 76 for all models.
Candidate ModelKr2adjAICc∆AICcwAICcER
log10(Wb) ~ log10(A) + log10(MS) + log10(MP) + MA60.45−232.263.830.106.80
log10(Wb) ~ log10(A) + log10(MS) + log10(MP)50.45−233.992.100.232.86
log10(Wb) ~ log10(A) + log10(MS)40.46−236.6600.671
log10(Wb) ~ log10(A)30.45−236.030.070.651.03
Table 2. Comparison of AICc deltas and weights for linear models testing the effects of bankfull width (Wb), sinuosity index (SI), river length (L), and bankfull channel Surface (Sb) and response (differences between real and estimated bankfull channel surface, Sb - Sbe). Shown are the number of parameters (K), the AICc, the ΔAICc, AICc weight (wAICc), and the evidence ratios (ERs). The best model and models roughly equivalent to the best (ΔAICc ≈ 2 or lower) are highlighted in bold.
Table 2. Comparison of AICc deltas and weights for linear models testing the effects of bankfull width (Wb), sinuosity index (SI), river length (L), and bankfull channel Surface (Sb) and response (differences between real and estimated bankfull channel surface, Sb - Sbe). Shown are the number of parameters (K), the AICc, the ΔAICc, AICc weight (wAICc), and the evidence ratios (ERs). The best model and models roughly equivalent to the best (ΔAICc ≈ 2 or lower) are highlighted in bold.
Candidate ModelKr2adjAICc∆AICcwAICcER
(Sb−Sbe) ~ Wb + SI + L + Sb60.711482.244.260.088.43
(Sb−Sbe) ~ Wb *** + SI * + L50.711480.102.120.232.89
(Sb−Sbe) ~ Wb *** + SI **40.711477.9700.671
(Sb−Sbe) ~ Wb ***30.681485.797.810.0149.70
(Sb−Sbe) ~ SI30.011571.2393.250.00>100
Significance code for variables: ‘***’ 0.001, ‘**’ 0.01, ‘*’ 0.05, ‘.’ 0.1, ‘’ 1.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Coello Sanz, F.; Casals, F.; Sánchez-González, J.R. How Can Be Lotic Ecosystem Size More Precisely Estimated? Comparing Different Approximations in Pre-Pyrenean and Pyrenean Mountains. Water 2021, 13, 721. https://doi.org/10.3390/w13050721

AMA Style

Coello Sanz F, Casals F, Sánchez-González JR. How Can Be Lotic Ecosystem Size More Precisely Estimated? Comparing Different Approximations in Pre-Pyrenean and Pyrenean Mountains. Water. 2021; 13(5):721. https://doi.org/10.3390/w13050721

Chicago/Turabian Style

Coello Sanz, Fernando, Frederic Casals, and Jorge Rubén Sánchez-González. 2021. "How Can Be Lotic Ecosystem Size More Precisely Estimated? Comparing Different Approximations in Pre-Pyrenean and Pyrenean Mountains" Water 13, no. 5: 721. https://doi.org/10.3390/w13050721

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

Article Metrics

Back to TopTop