Next Article in Journal
Fixed-Bed and Mobile-Bed Resistance of Channels with Steep Gradients in Mountainous Areas
Next Article in Special Issue
Diets and Trophic Structure of Fish Assemblages in a Large and Unexplored Subtropical River: The Uruguay River
Previous Article in Journal
Effect of Pilot-Points Location on Model Calibration: Application to the Northern Karst Aquifer of Qatar
Previous Article in Special Issue
Periods of Extreme Shallow Depth Hinder but Do Not Stop Long-Term Improvements of Water Quality in Lake Apopka, Florida (USA)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Beta Diversity Partitioning and Drivers of Variations in Fish Assemblages in a Headwater Stream: Lijiang River, China

1
College of Environmental Science and Engineering, Guilin University of Technology, Guilin 541004, China
2
Guangxi Key Laboratory of Environmental Pollution Control Theory and Technology, Guilin 541004, China
3
Collaborative Innovation Center for Water Pollution Control and Water Safety in Karst Area, Guilin University of Technology, Guilin 541004, China
4
Department of Bioscience, Aarhus University, Vejlsøvej 25, 8600 Silkeborg, Denmark
5
Sino-Danish Centre for Education and Research, Beijing 100049, China
6
Department of Hydrology and Water Resources Management, Kiel University, Kiel 24118, Germany
7
Department of Health and Environmental Sciences, Xi’an Jiaotong-Liverpool University, Suzhou 215123, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(4), 680; https://doi.org/10.3390/w11040680
Submission received: 21 January 2019 / Revised: 9 February 2019 / Accepted: 1 April 2019 / Published: 2 April 2019
(This article belongs to the Special Issue Trophic Interactions in Warm Freshwater Ecosystems)

Abstract

:
Beta diversity partitioning has currently received much attention in research of fish assemblages. However, the main drivers, especially the contribution of spatial and hydrological variables for species composition and beta diversity of fish assemblages are less well studied. To link species composition to multiple abiotic variables (i.e., local environmental variables, hydrological variables, and spatial variables), the relative roles of abiotic variables in shaping fish species composition and beta diversity (i.e., overall turnover, replacement, and nestedness) were investigated in the upstream Lijiang River. Species composition showed significant correlations with environmental, hydrological, and spatial variables, and variation partitioning revealed that the local environmental and spatial variables outperformed hydrological variables, and especially abiotic variables explained a substantial part of the variation in the fish composition (43.2%). The overall species turnover was driven mostly by replacement (87.9% and 93.7% for Sørensen and Jaccard indices, respectively) rather than nestedness. Mantel tests indicated that the overall species turnover (ßSOR and ßJAC) and replacement (ßSIM and ßJTU) were significantly related to hydrological, environmental, and spatial heterogeneity, whereas nestedness (ßSNE or ßJNE) was insignificantly correlated with abiotic variables (P > 0.05). Moreover, the pure effect of spatial variables on overall species turnover (ßSOR and ßJAC) and replacement (ßSIM and ßJTU), and the pure effect of hydrological variables on replacement (ßSIM and ßJTU), were not important (P > 0.05). Our findings demonstrated the relative importance of interactions among environmental, hydrological, and spatial variables in structuring fish assemblages in headwater streams; these fish assemblages tend to be compositionally distinct, rather than nested derivatives of one another. Our results, therefore, indicate that maintaining natural flow dynamics and habitat continuity are of vital importance for conservation of fish assemblages and diversity in headwater streams.

1. Introduction

Fish assemblages in streams and rivers are, to a large degree determined by the prevailing spatial and environmental conditions, including physical, chemical, and biological characteristics [1,2]. Several studies have focused on the spatial variation of fish assemblages in rivers and their relationship to abiotic variables [3,4], and most of studies support the river continuum concept, description of the changes in fish composition from upstream to downstream [5,6]. However, temporal variations of fish assemblages occur in some sites [3,7], but not in others [5,8]. Although the relationship between fish assemblage and abiotic variables has been intensively investigated, the relative role of different variables contributing to fish assemblage variations remains poorly elucidated [9,10]. Some studies found that hydromorphologic variables (i.e., channel geometry, substrate, flow, vegetation, and land use in the riparian zone) were dominant regulating factors [11,12]. However, local environmental variables (e.g., nutrients, dissolved oxygen, pH, conductivity) and climate are more often considered to be the key determining factors [5,10]. It is, therefore, uncertain which of the abiotic variables (i.e., environmental, hydrological, and spatial variables) best explain the fish assemblage patterns on the regional scale [13,14]. Moreover, knowledge of the influence of regional constraints and historical processes on species contribution, abundance, and diversity remain superficial, although it has been considered critical for an advance in stream fish ecology and conservation [15,16].
Beta diversity (ß-diversity) describes the variation in species composition among communities, measured by ecological dissimilarities between sampling units based on Jaccard or Sørensen indices [17,18]. ß-diversity can be further partitioned into two additive components accounting for (i) the dissimilarity associated with the replacement of some species by others between assemblages (i.e., species replacement) and (ii) the dissimilarity associated with species losses in which an assemblage is a strict subset of another (i.e., nestedness-resultant dissimilarity of assemblages [19,20]. Some studies have revealed that replacement contributes more than nestedness to the overall species turnover (ß-diversity) of fish assemblages in rivers and streams [21,22], while others have shown an equal contribution of replacement and nestedness to ß-diversity [13,23]. Thus, the contribution of replacement or nestedness to ß-diversity of fish assemblages may be area-specific.
Most Chinese studies of fish assemblages in rivers have focused on fish composition and its distribution, spatially and temporally, and on the role of environmental variables [24,25,26,27]. Recent studies indicate that fish assemblages correlate significantly with abiotic variables, such as river width, discharge, velocity, water temperature, dissolved oxygen, and altitude [27,28]. However, the question of which variables (i.e., physical variables, chemical variables, hydrological variables, and spatial variables) contribute most to fish assemblage variations is not wellelucidated. In addition, in China, α-diversity within fish assemblages (i.e., species richness, Shannon index, Simpson’s index, etc.) is widely studied [5,29,30], while ß-diversity is poorly documented, especially ß-diversity partitioning [30]. Exploring the two additive components of ß-diversity, replacement and nestedness, can help with the acquisition of information about the underlying ecological mechanisms at different environmental conditions [20].
In the present study, we attempt to quantify the role of different abiotic variables in determining changes in fish assemblages through partitioning of ß-diversity in the upstream Lijiang River, an area where information on freshwater fish diversity and distribution is currently poor. We address three questions: (i) What is most important for fish species composition: spatial variables, environmental variables, or hydrological variables? (ii) What are the major drivers of species and ß-diversity of fish assemblages? (iii) Which factors contribute most to ß-diversity: replacement or nestedness? The hypotheses were that (i) environmental, hydrological, and spatial variables interacted to determine species composition and ß-diversity, and that (ii) spatial variables would be a key driver of ß-diversity, but it was uncertain whether replacement or nestedness would prove to be of superior importance.

2. Materials and Methods

2.1. Study Area

Lijiang River, which lies in Guangxi Zhuang Autonomous Region, southwest China (Figure 1), originates from the main peak of Mao’er Mountain (maximum altitude 2141.5 m) in the northeast of Guilin, and has a total length of 214 km and drainage area of 2173.29 km2. Lijiang River has typical rain source characteristics and the bankfull discharge is 7810 m3/s, usually occurring in May and June [31]. The river has an extremely uneven flow during the year, mainly due to the huge heavy rainfall during the southeast monsoon in the upstream of the river flowing in a steep watershed. Ludong Stream, with a mean slope of 7.38‰, is one of the upstream tributaries to Lijiang River. The total length of Ludong Stream is 46 km with a total of drainage area of 314 km2. All stream sections in the study area are mountainous with a turbulent waterflow and a large hydraulic gradient (31.7‰). The average annual discharge is 5.99 × 108 m3 and the average annual flow is 19 m3/s. However, the Fuzikou Reservoir was established in the Ludong Stream by means of a dam (with a total length of 239 m and a maximum height of 76.5 m) in January 2018 for the purpose of flood control, ecological hydration of Lijiang River, and power generation. The reservoir has a water storage capacity of 1.88 × 108 m3 with an impoundment level of 267 m, which can compensate water discharge to Lijiang River with a total of 1.114 × 108 m3 annually. The vegetation types on both sides of the stream are mixed, mostly consisting of coniferous forest, broad-leaved forest, and bamboo forest [32]. The fish community in most upstream areas of Lijiang River is dominated by Zacco platypus, Pseudogastromyzon fangi, Erromyzon sinensis, and Opsariichthys bidens, with little temporal variation in assemblages [33]. No floating and submerged vegetation was observed.

2.2. Field Sample Collection

Fish samples were collected from 26 sampling sites (Figure 1) from September to November 2016. All fish in each sampling reach, which was about 150 m long, were collected by a backpack electro-fisher (CWB-2000P, China; 12V import, 250 V export) downstream to upstream, with each ending of the reach blocked by block net. Fish specimens longer than 20 mm were identified to species, counted, and then returned to the river. Fish with total length smaller than 20 mm were not included and were released on site [34,35].
About 15 abiotic variables were measured at each sampling site, including physical and chemical parameters, such as water temperature, pH, dissolved oxygen, conductivity, turbidity, substrate size, and habitat structure character. Water temperature, pH, dissolved oxygen and conductivity were measured using a portable water quality analysis device (HACH HQ 40 d), and turbidity was measured applying a portable turbidimeter (HACH 2100 Q). Velocity was recorded at middle depth using a printing flow velocity meter (LJD-10 A). Latitude, longitude, and altitude were obtained from GPS (Garmin eTrex 10), and the distances of each sampling site to the dam of Fuzikou Reservoir were obtained from Google Earth. Measurements of stream wetted width and depth were taken at the widest and deepest section, and habitat percentage of habitat was estimated by walking along the stream reach after sampling and recording the types and relative proportions of the habitats (i.e., cascade, riffle, pool, step pool, run, glide). The same procedure was applied for substrate size, which was calculated by measuring the median particle diameter of thirty randomly chosen particles in each sampling reach [36].

2.3. Biotic and Abiotic Datasets

Two biotic datasets were applied in this study: (i) A species abundance data and species presence/absence data—a total of all 37 fish species recorded and their abundance were included and (ii) beta diversities (ß-diversity). To calculate pairwise dissimilarities, we used the Sørensen dissimilarity index (ßSOR) and the Jaccard dissimilarity index (ßJAC) and two additive components (i.e., ßSIM, ßNES, and ßSOR for the Sørensen index and ßJAC, ßJTU, and ßJNE for the Jaccard index) [19,37]. Abbreviations used in the paper are explained in Table 1.

2.4. Abiotic Datasets

Three abiotic datasets were compiled as well including hydrological variables (Hv), environmental variables (Ev), and spatial variables (Sv). Hydrological variables (Hv) were in situ measured stream wet width, maximum depth, current velocity at the sampling sites, stream order, altitude, and the distance from the sampling sites to the dam of Fuzikou Reservoir. Environmental variables (Ev) were pH, water temperature (T), dissolved oxygen (DO), conductivity, turbidity and substrate size, hydrostatic habitats (%), fast-flow habitats (%), and slow-flow habitats (%). Except for the coordinates (X: latitude, Y: longitude), Moran’s eigenvector maps were used to generate the spatial variables (Sv) representing geographical positions and dispersal across the streams. A Euclidean distance matrix between each pair of sampling sites was calculated using the earth.dist function in the R package fossil (version 0.3.7). Then a principal coordinates of neighborhood matrix (PCNM) analysis was applied to compute the spatial variables representing geographical positions through the pcnm function in the R package vegan (version 2.5-3). The generated eigenvectors were considered as spatial variables (i.e., PCNMs), which could reflect unmeasured broadscale variation in the modern environment or historic factors, for example, natural dispersal-generated patterns demonstrating internal local-scale dispersal dynamics or regional-scale migration history [38]. Finally, a total of 17 PCNMs were generated, which were used in the following analyses (Table 2).

2.5. Data Analysis

To assess the species richness in the study area, species accumulation curves based on samples and individuals were applied. We used the specpool function in the R package vegan to get extrapolated species richness estimates based on presence/absence data, and species rarefaction curves were plotted based on abundance data with the specaccum function in the R package vegan.
To explore the potential impacts of spatial variables, hydrological variables and environmental variables on fish assemblages (question i), preliminary data analyses were conducted. Firstly, species composition with abundance was Hellinger-transformed using the function decosdtand in the R package vegan. Secondly, a detrended correspondence analysis (DCA) using the function decorana on the Hellinger-transformed species data produced the longest gradient length of 2.13 along the first axis, suggesting that redundancy analysis (RDA) was appropriate for further analyses [39]. We performed RDA with the rda function and tested the significance using the anova function. If significant, a forward selection using the function ordistep could be utilized to obtain reduced variables for the spatial variables, hydrological variables, and environmental variables, respectively. To eliminate factors with high collinearity, all abiotic variables were log(x+1) transformed, and variance inflation factors (VIF) using the vif function in the R package car were computed, and any factor with a VIF > 10 and a Spearman’s rank correlation coefficient (|r| > 0.75) was removed from the analysis [40]. All the selected variables (i.e., Ev, Hv, and Sv) were used as explanatory variables and Hellinger-transformed fish composition data were used for the variation partitioning analysis using the function varpart in the R package vegan [41]. The significance of the fractions of variation explained was tested with 9999 permutations [42].
Next, we ran Mantel tests in order to examine the changes in species composition along spatial, hydrological, and environmental gradients (question ii). The Mantel test has been utilized as a distance-based approach to study ß-diversity in relation to distance matrices [40]. The significance of this distancedecay relationship, which measures how dissimilarity decays with increasing distance between pairwise sites, was determined using Mantel test with 9999 permutations [42]. We first calculated the total species turnover (ßSOR or ßJAC) and its two additive components (replacement (ßSIM or ßJTU) and nestedness (ßNES or ßJNE)) based on presence/absence data with the function beta.multi in the R package betapart (version 1.5.1), and we used the function beta.pair to compute three dissimilarity matrices for ß-diversity (i.e., ßSIM, ßNES, and ßSOR for index.family = “sørensen”, and ßJAC, ßJTU, and ßJNE for index.family = “jaccard”). Dissimilarity matrices based on Euclidean distances for the spatial, hydrological, and environmental variables (i.e., Svdis, Hvdis, and Evdis) were constructed. In addition to simple Mantel tests using six matrices, we applied partial Mantel tests to ease apart the pure effects of spatial, hydrological, and environmental variables on ß-diversity matrices, and the significance was assessed using 9999 permutations. Mantel and partial Mantel tests were run using the functions mantel and mantel.partial, respectively, in the R package vegan [40,43]. All analyses were performed with the R software (version 3.5.1, R Development Core Team 2018).

3. Results

3.1. Variability of Abiotic Variables

The hydrological and environmental variables varied widely in the river (Table 2). For instance, stream wet width ranged from 2.4 m to 63.5 m (mean: 14.2 m), current velocity from 14.2 cm/s to 59.4 cm/s (mean: 31.2 cm/s), water temperature from 18.8 °C to 25.1 °C (mean: 22.1 °C), and conductivity from 16.1 μS/cm to 157.4 μS/cm (mean: 62.4 μS/cm). However, the spatial variables showed a small variation with latitude ranging from 25.70 to 25.88°N and longitude from 110.50 to 110.60°E, which reflects the small catchment (314 km2) studied.

3.2. Fish Assemblage Composition

A total of 2968 fish individuals were collected, comprising 37 fish species belonging to 11 families, most of which were native to the drainage area (Table 3). Cyprinidae, Balitoridae and Gobiidae contributed significantly to the typical fish assemblage composition, the proportions of individuals being 40.84%, 39.08%, and 14.82%, respectively. The dominant species were Pseudogastromyzon fangi (24.8%), Zacco platypus (21.0%), Acrossocheilus parallens (12.8%), and Erromyzon sinensis (10.5%). Of these species, 10 species were singleton species and 5 were doubleton species. Species richness ranged from 2 to 18 per site, and mean species richness was 11.19. In addition, the individuals ranged from 27 to 223 with an average of 114.15 individuals per site.
The species rarefaction curve showed no apparent asymptote, but a trend towards stabilization (Figure 2). Species richness in the areas was estimated to 47 ± 8 using the function specpool in the R package vegan, indicating that 78.7% of the expected richness were collected in this study.

3.3. Driver of Species Composition

In the RDA analysis for species composition, environmental, hydrological, and spatial variables all showed a significant relationship with composition. Two environmental variables, two hydrological variables, and four spatial variables were selected by forward selection (Table 4). According to variation partitioning analysis, the three abiotic variable sets explained 43.2% of species composition (Figure 3). The pure effect of environmental variables (7.1%) and spatial variables (4.0%) contributed more than the pure effect of hydrological variables (0.7%) on species variance, the joint effect of the three sets being the largest with 18.8%. The variation partitioning showed that fish composition was determined by the interaction of hydrological, environmental and spatial variables, thus supporting our first hypothesis.

3.4. Main Driver of ß-Diversity

The partitioning total dissimilarity for the Sørensen and Jaccard indices showed that replacement (ßSIM or ßJTU) contributes more to the overall species turnover than nestedness (ßSNE or ßJNE). For example, the Sørensen indices (ßSIM = 0.767; ßSNE = 0.096; ßSOR = 0.863) suggested that replacement alone was responsible for 87.9% of the overall turnover of the species composition, and only 11.1% was accounted for by nestedness. Furthermore, the Jaccard indices (ßJTU = 0.868; ßJNE = 0.058; ßJAC = 0.926) showed that replacement alone was responsible for 93.7% of the overall turnover of the species composition, with only 6.3% being accounted for by nestedness.
Mantel tests showed that the overall turnover of species composition, based on both the Sørensen and Jaccard indices (ßSOR and ßJAC), increased significantly with hydrological and environmental distances, and the relationships between the overall species turnover (ßSOR and ßJAC) and hydrological distances were similar to the relationships between the overall turnover (ßSOR and ßJAC) and environmental distances. The weakest relationships were found for spatial distances (Table 5). The replacement based on both indices (ßSIM and ßJTU) increased significantly with hydrological, spatial, and environmental distances (excepting ßJTU) (Figure 4 and Figure 5), and the relationships between replacement (ßSIM and ßJTU) and hydrological distances were similar to the relationships between replacement (ßSIM and ßJTU) and spatial distances. The weakest relationships were those for environmental distances (Figure 4 and Figure 5). In contrast, there were no significant relationships between nestedness (ßSNE or ßJNE) and environmental, hydrological, and spatial distances (P > 0.05) (Table 5).
Based on partial Mantel tests, the pure effects of hydrological and environmental distances on the overall species turnover (ßSOR and ßJAC) were significant, whereas the pure effect of spatial distances was not significant (P > 0.05). The pure effect of environmental distances on replacement (ßSIM and ßJTU) was significant (P < 0.05), while the pure effects of hydrological and spatial distances were insignificant. Similarly, there were no significant relationships between nestedness (ßSNE or ßJNE) and environmental, hydrological, and spatial distances (P > 0.05) (Table 5).

4. Discussion

An important task in ecology is to explore the variables controlling the abundance and distribution patterns of aquatic organisms and the underlying causes [40]. Our results suggest that environmental filtering is a significant force structuring fish assemblage. The environmental, hydrological and spatial variables together explained approximately 43.2% of the total variation in fish composition. The most significant factors determining fish assemblages in upstream Lijiang River were the interactions among environmental variables, hydrological variables, and spatial variables accounting in total for 18.8% of the variation. Moreover, the pure environmental variables contributed much more than the pure spatial and hydrological variables (Figure 3). This concurs with other studies showing environmental variables to be more important than spatial variables in explaining fish composition. Earlier studies have revealed that the amount of variation explained by environmental variables, including elevation and stream width, among others, ranged from 17.3% to 34.9% [14,44,45,46]. However, we found the amount of variation explained by environmental variables to be lower (i.e., 7.1%) than recorded in these previous studies, which might reflect that we separated various hydrological variables (i.e., stream width, depth, flow velocity, altitude, etc.) from environmental variables, potentially increasing the overall interactions among the three different abiotic sets.
Longitudinal patterns in fish species richness and assemblage structure have been well studied in American and European streams [47,48], as well as in China [5,24,49,50]. In wadable headwater streams, local fish species richness usually increases downstream, which is due to the general increase in stream size and habitat diversity and complexity downstream [5,47,51]. We found that fish species richness, abundance, and Shannon index increased downstream, but exhibiting maximum values in the second stream order followed by a decrease towards third stream order (Figure 6), as found also in other Chinese headwater streams [50]. A subsequently lower fish species richness, abundance and Shannon index values in the third stream order might reflect anthropogenic disturbances, most importantly construction of the Fuzikou Reservoir dam(i.e., three sampling sites in the third stream order being located less than 4 km from the dam). Besides, many rare species were found in the present study area, resulting in no apparent asymptote in species rarefaction curves (Figure 2) [52,53]. A total of 14 fish species were collected from just one or two locations, representing as much as 37.8% of all species found.
We used two approaches to calculate the overall species turnover and its two additive components (replacement and nestedness) based on the Sørensen and the Jaccard indices. Replacement was more important than nestedness in explaining the overall species turnovers, implying that species replacement was responsible for the pattern in species turnover in the upstream Lijiang River. This indicates that the variation in species between sampling sites and the compositional difference in individual assemblages were the main drivers of diversity at a watershed scale [14,43]. The pattern in species turnover concurs with other recent findings from studies on ß-diversity of stream fish assemblages [14,22,30].
Additionally, identifying the mechanisms underlying the spatial patterns of biodiversity is a major task in community ecology, especially with regards to biodiversity conservation and stream restoration [30,40,54]. The relative roles of different distances showed considerable variability for ß-diversity and the additive components. For instance, Mantel tests suggested that the following rank of importance of the variables accounting for the among-site differences in species based on both indices (ßSOR or ßJAC): Environmental variables ≈ hydrological variables > spatial variables, while the importance for replacement in species (ßSIM or ßJTU) was: Environmental variables ≈ spatial variables > hydrological variables. In contrast, nestedness (ßSNE or ßJNE) was not significantly correlated with the environmental, hydrological, and spatial variables (Table 5). Our results revealed a clear distancedecay of fish assemblage dissimilarity with increasing environmental and hydrological distances (Figure 4 and Figure 5), which indicated that maintaining various habitats continuously in streams are of vital importance.
The Fuzikou Reservoir started to be impounded in January 2018 and when filled it will have a storage capacity of 1.88 × 108 m3. By then, all sampling sites in the third stream order and most sampling sites in the second stream order will potentially be affected, which may reduce the fish species richness due to the increase of permanent lentic habitats and the decrease in hydrological connectivity [28,55,56,57]. Moreover, the dam may weaken the influence of stream spatial position on fish species richness [28], which could lead to the replacement of species turnover by nestedness [22,55].
In conclusion, the present study reveals that the importance of different abiotic variables in constructing fish assemblage composition and ß-diversity patterns. Maintenance of naturally existing ß-diversity patterns is of vital importance in the conservation of regional species diversity pool [58]. Our findings identified the key drivers of fish composition being the interactions between environmental, hydrological, and spatial variables, and further revealed that species replacement contributes more than nestedness in maintaining ß-diversity. Our results, therefore, indicate that maintaining natural flow dynamics and habitat continuity are of vital importance for conservation of fish assemblages and diversity in headwater streams and that the dam construction in such areas, therefore, induces conservation problems.

Author Contributions

This paper received individual contributions from each author as specified: L.H., J.H. and Q.Z. developed the basic idea and designed the investigation. J.H., L.H., Y.M., Z.W. and Q.Z. performed the sampling collection and processing. L.H., Y.M., N.W. performed the data analyses. L.H. and J.H. led the original writing. E.J., Z.W. and N.W. reviewed and made corrections to improve the paper. All authors read the full paper and agreed to its publication.

Acknowledgments

This study was supported by National Natural Science Foundation of China (Grant/Award Number: 51509042), Guangxi Natural Science Foundation (Grant/Award Numbers: 2016GXNSFAA380104; 2018GXNSFAA281022), Guangxi ‘Bagui Scholar’ Construction Project (Grant/Award Number: 2016A10), and Guangxi Science and Technology Planning Project (GuiKe-AD18126018). Erik Jeppesen was supported by WATEC-Centre for Water Technology, Aarhus University. We also thank Anne Mette Poulsen for English revision, and two anonymous reviewers’ constructive comments which greatly improved our manuscript.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding this publication.

References

  1. Jowett, I.G.; Richardson, J. Fish communities in New Zealand rivers and their relationship to environmental variables. N. Z. J. Mar. Freshw. Res. 2003, 37, 347–366. [Google Scholar] [CrossRef] [Green Version]
  2. Zeni, J.O.; Casatti, L. The influence of habitat homogenization on the trophic structure of fish fauna in tropical streams. Hydrobiologia 2014, 726, 259–270. [Google Scholar] [CrossRef]
  3. Silvano, R.A.M.; do Amaral, B.D.; Oyakawa, O.T. Spatial and temporal patterns of diversity and distribution of the Upper Jurua River fish community (Brazilian Amazon). Environ. Biol. Fishes 2000, 57, 25–35. [Google Scholar] [CrossRef]
  4. Chick, J.H.; Pegg, M.A.; Koel, T.M. Spatial patterns of fish communities in the Upper Mississippi River System: Assessing fragmentation by low-head dams. River Res. Appl. 2006, 22, 413–427. [Google Scholar] [CrossRef]
  5. Li, J.; Huang, L.; Zou, L.; Kano, Y.; Sato, T.; Yahara, T. Spatial and temporal variation of fish assemblages and their associations to habitat variables in a mountain stream of north Tiaoxi River, China. Environ. Biol. Fishes 2012, 93, 403–417. [Google Scholar] [CrossRef]
  6. Bennett, M.G.; Kozak, J.P. Spatial and temporal patterns in fish community structure and abundance in the largest US river swamp, the Atchafalaya River floodplain, Louisiana. Ecol. Freshw. Fish 2016, 25, 577–589. [Google Scholar] [CrossRef]
  7. Espirito-Santo, H.M.V.; Zuanon, J. Temporary pools provide stability to fish assemblages in Amazon headwater streams. Ecol. Freshw. Fish 2017, 26, 475–483. [Google Scholar] [CrossRef]
  8. Pegg, M.A.; McClelland, M.A. Spatial and temporal patterns in fish communities along the Illinois River. Ecol. Freshw. Fish 2004, 13, 125–135. [Google Scholar] [CrossRef]
  9. Senay, C.; Taranu, Z.E.; Bourque, G.; Macnaughton, C.J.; Lanthier, G.; Harvey-Lavoie, S.; Boisclair, D. Effects of river scale flow regimes and local scale habitat properties on fish community attributes. Aquat. Sci. 2017, 79, 13–26. [Google Scholar] [CrossRef]
  10. Daufresne, M.; Boet, P. Climate change impacts on structure and diversity of fish communities in rivers. Glob. Chang. Biol. 2007, 13, 2467–2478. [Google Scholar] [CrossRef]
  11. Wyzga, B.; Amirowicz, A.; Radecki-Pawlik, A.; Zawiejska, J. Hydromorphological Conditions, Potential Fish Habitats and the Fish Community in a Mountain River Subjected to Variable Human Impacts, the Czarny Dunajec, Polish Carpathians. River Res. Appl. 2009, 25, 517–536. [Google Scholar] [CrossRef]
  12. Fischer, J.; Paukert, C.; Daniels, M. Fish Community Response to Habitat Alteration: Impacts of Sand Dredging in the Kansas River. Trans. Am. Fish. Soc. 2012, 141, 1532–1544. [Google Scholar] [CrossRef] [Green Version]
  13. Shukla, R.; Bhat, A. Beta-diversity partitioning and drivers of variations in tropical fish community structure in central India. Aquat. Sci. 2018, 80, 18. [Google Scholar] [CrossRef]
  14. Zbinden, Z.D.; Matthews, W.J. Beta diversity of stream fish assemblages: Partitioning variation between spatial and environmental factors. Freshw. Biol. 2017, 62, 1460–1471. [Google Scholar] [CrossRef]
  15. Fausch, K.D.; Torgersen, C.E.; Baxter, C.V.; Li, H.W. Landscapes to riverscapes: Bridging the gap between research and conservation of stream fishes. Bioscience 2002, 52, 483–498. [Google Scholar] [CrossRef]
  16. Leite, G.F.M.; Silva, F.T.C.; Goncalves, J.F.J.; Salles, P. Effects of conservation status of the riparian vegetation on fish assemblage structure in neotropical headwater streams. Hydrobiologia 2015, 762, 223–238. [Google Scholar] [CrossRef]
  17. Legendre, P. Interpreting the replacement and richness difference components of beta diversity. Glob. Ecol. Biogeogr. 2014, 23, 1324–1334. [Google Scholar] [CrossRef] [Green Version]
  18. Koleff, P.; Gaston, K.J.; Lennon, J.J. Measuring beta diversity for presence–absence data. J. Anim. Ecol. 2003, 72, 367–382. [Google Scholar] [CrossRef]
  19. Baselga, A. Partitioning the turnover and nestedness components of beta diversity. Glob. Ecol. Biogeogr. 2010, 19, 134–143. [Google Scholar] [CrossRef]
  20. Baselga, A.; Leprieur, F. Comparing methods to separate components of beta diversity. Methods Ecol. Evol. 2015, 6, 1069–1079. [Google Scholar] [CrossRef]
  21. Griffiths, D. Connectivity and vagility determine beta diversity and nestedness in North American and European freshwater fish. J. Biogeogr. 2017, 44, 1723–1733. [Google Scholar] [CrossRef]
  22. Vitorino, O.B.; Fernandes, R.; Agostinho, C.S.; Pelicice, F.M. Riverine networks constrain beta-diversity patterns among fish assemblages in a large Neotropical river. Freshw. Biol. 2016, 61, 1733–1745. [Google Scholar] [CrossRef]
  23. Leprieur, F.; Tedesco, P.A.; Hugueny, B.; Beauchard, O.; Durr, H.H.; Brosse, S.; Oberdorff, T. Partitioning global patterns of freshwater fish beta diversity reveals contrasting signatures of past climate changes. Ecol. Lett. 2011, 14, 325–334. [Google Scholar] [CrossRef] [PubMed]
  24. Kang, B.; He, D.M.; Perrett, L.; Wang, H.Y.; Hu, W.X.; Deng, W.D.; Wu, Y.F. Fish and fisheries in the Upper Mekong: Current assessment of the fish community, threats and conservation. Rev. Fish Biol. Fish. 2009, 19, 465–480. [Google Scholar] [CrossRef]
  25. Xie, J.Y.; Tang, W.J.; Yang, Y.H. Fish assemblage changes over half a century in the Yellow River, China. Ecol. Evol. 2018, 8, 4173–4182. [Google Scholar] [CrossRef] [PubMed]
  26. Xing, Y.C.; Zhang, C.G.; Fan, E.Y.; Zhao, Y.H. Freshwater fishes of China: Species richness, endemism, threatened species and conservation. Divers. Distrib. 2016, 22, 358–370. [Google Scholar] [CrossRef]
  27. Shuai, F.M.; Li, X.H.; Chen, F.C.; Li, Y.F.; Lek, S. Spatial patterns of fish assemblages in the Pearl River, China: Environmental correlates. Fundam. Appl. Limnol. 2017, 189, 329–340. [Google Scholar] [CrossRef]
  28. Yan, Y.Z.; Xiang, X.Y.; Chu, L.; Zhan, Y.J.; Fu, C.Z. Influences of local habitat and stream spatial position on fish assemblages in a dammed watershed, the Qingyi Stream, China. Ecol. Freshw. Fish 2011, 20, 199–208. [Google Scholar] [CrossRef]
  29. Kang, B.; Deng, J.M.; Wu, Y.F.; Chen, L.Q.; Zhang, J.; Qiu, H.Y.; Lu, Y.; He, D.M. Mapping China’s freshwater fishes: Diversity and biogeography. Fish Fish. 2014, 15, 209–230. [Google Scholar] [CrossRef]
  30. Kang, B.; Huang, X.X.; Yan, Y.Z.; Yan, Y.R.; Lin, H.D. Continental-scale analysis of taxonomic and functional fish diversity in the Yangtze river. Glob. Ecol. Conserv. 2018, 15, e00442. [Google Scholar] [CrossRef]
  31. Li, R.; Chen, Q.; Tonina, D.; Cai, D. Effects of upstream reservoir regulation on the hydrological regime and fish habitats of the Lijiang River, China. Ecol. Eng. 2015, 76, 75–83. [Google Scholar] [CrossRef]
  32. Wu, Z.; Huang, L. The Relationship between Fish Diversity and Environmental Variables; Science Press: Beijing, China, 2017. [Google Scholar]
  33. Zhu, Z.; Wu, Z.; Huang, L.; Feng, W.; Ding, Y. Species composition and diversity of fish in the upper reach of Lijiang River. Sichuan J. Zool. 2015, 34, 126–132. [Google Scholar]
  34. Anonymous. Freshwater Fishes of Guangxi, China, 2nd ed.; Guangxi People’s Publishing House: Nanning, China, 2006. [Google Scholar]
  35. Barbour, M.T.; Gerritsen, J.; Snyder, B.D.; Stribling, J.B. Rapid Bioassessment Protocols for Use in Streams and Wadable Rivers: Periphyton, Benthic Invertebrates and Fish, 2nd ed.; U.S. Environmental Protection Agency, Office of Water: Washington, DC, USA, 1999.
  36. Harrelson, C.C.; Rawlins, C.L.; Potyondy, J.P. Stream Channel Reference Sites: An Illustrated Guide to Field Technique; General Technical Report RM-245; U.S. Department of Agriculture, Forest Service, Rocky Mountain Forest and Range Experiment Station: Fort Collins, CO, USA, 1994.
  37. Baselga, A.; Orme, C.D.L. betapart: An R package for the study of beta diversity. Methods Ecol. Evol. 2012, 3, 808–812. [Google Scholar] [CrossRef]
  38. Svenning, J.C.; Baktoft, K.H.; Balslev, H. Land-use history affects understorey plant species distributions in a large temperate-forest complex, Denmark. Plant Ecolgy 2009, 201, 221–234. [Google Scholar] [CrossRef]
  39. Lepš, J.; Šmilauer, P. Multivariate Analysis of Ecological Data Using CANOCO; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  40. Wu, N.; Qu, Y.; Guse, B.; Makarevičiūtė, K.; To, S.; Riis, T.; Fohrer, N. Hydrological and environmental variables outperform spatial factors in structuring species, trait composition, and beta diversity of pelagic algae. Ecol. Evol. 2018, 8, 2947–2961. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Peres-Neto, P.R.; Legendre, P.; Dray, S.; Borcard, D. Varation partitioning of species data matrices: Estimation and comparison of fractions. Ecology 2006, 87, 2614–2625. [Google Scholar] [CrossRef]
  42. Legendre, P. Studying beta diversity: Ecological variation partitioning by multiple regression and canonical analysis. J. Plant Ecol. 2008, 1, 3–8. [Google Scholar] [CrossRef]
  43. Baselga, A. The relationship between species replacement, dissimilarity derived from nestedness, and nestedness. Glob. Ecol. Biogeogr. 2012, 21, 1223–1232. [Google Scholar] [CrossRef] [Green Version]
  44. Sály, P.; Takács, P.; Kiss, I.; Bíró, P.; Erős, T. The relative influence of spatial context and catchment- and site-scale environmental factors on stream fish assemblages in a human-modified landscape. Ecol. Freshw. Fish 2011, 20, 251–262. [Google Scholar] [CrossRef]
  45. Godinho, F.N.; Ferreira, M.T.; Santos, J.M. Variation in fish community composition along an Iberian river basin from low to high discharge: Relative contributions of environmental and temporal variables. Ecol. Freshw. Fish 2000, 9, 22–29. [Google Scholar] [CrossRef]
  46. Magalhães, M.F.; Batalha, D.C.; Collares-Pereira, M.J. Gradients in stream fish assemblages across a Mediterranean landscape: Contributions of environmental factors and spatial structure. Freshw. Biol. 2002, 47, 1015–1031. [Google Scholar] [CrossRef]
  47. Matthews, W.J. Patterns in Freshwater Fish Ecology; Chapman & Hall: New York, NY, USA, 1998. [Google Scholar]
  48. Jackson, D.A.; Peres-Neto, P.R.; Olden, J.D. What controls who is where in freshwater fish communities—The roles of biotic, abiotic, and spatial factors. Can. J. Fish. Aquat. Sci. 2001, 58, 157–170. [Google Scholar]
  49. Fu, C.Z.; Wu, J.H.; Chen, J.K.; Qu, Q.H.; Lei, G.C. Freshwater fish biodiversity in the Yangtze River basin of China: Patterns, threats and conservation. Biodivers. Conserv. 2003, 12, 1649–1685. [Google Scholar] [CrossRef]
  50. Li, Y.R.; Tao, J.; Chu, L.; Yan, Y.Z. Effects of anthropogenic disturbances on alpha and beta diversity of fish assemblages and their longitudinal patterns in subtropical streams, China. Ecol. Freshw. Fish 2018, 27, 433–441. [Google Scholar] [CrossRef]
  51. 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]
  52. Colwell, R.K.; Chao, A.; Gotelli, N.J.; Lin, S.-Y.; Mao, C.X.; Chazdon, R.L.; Longino, J.T. Models and estimators linking individual-based and sample-based rarefaction, extrapolation and comparison of assemblages. J. Plant Ecol. 2012, 5, 3–21. [Google Scholar] [CrossRef] [Green Version]
  53. Chao, A. Non-parametric estimation of the number of classes in a population. Scand. J. Stat. 1984, 11, 265–270. [Google Scholar]
  54. Myers, N.; Mittermeier, R.A.; Mittermeier, C.G.; da Fonseca, G.A.B.; Kent, J. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853. [Google Scholar] [CrossRef]
  55. Araújo, E.S.; Marques, E.E.; Freitas, I.S.; Neuberger, A.L.; Fernandes, R.; Pelicice, F.M. Changes in distance decay relationships after river regulation: Similarity among fish assemblages in a large Amazonian river. Ecol. Freshw. Fish 2013, 22, 543–552. [Google Scholar] [CrossRef]
  56. Pelicice, F.M.; Pompeu, P.S.; Agostinho, A.A. Large reservoirs as ecological barriers to downstream movements of Neotropical migratory fish. Fish Fish. 2015, 16, 697–715. [Google Scholar] [CrossRef]
  57. Cooper, A.R.; Infante, D.M.; Daniel, W.M.; Wehrly, K.E.; Wang, L.; Brenden, T.O. Assessment of dam effects on streams and fish assemblages of the conterminous USA. Sci. Total Environ. 2017, 586, 879–889. [Google Scholar] [CrossRef] [PubMed]
  58. Socolar, J.B.; Gilroy, J.J.; Kunin, W.E.; Edwards, D.P. How Should Beta-Diversity Inform Biodiversity Conservation? Trends Ecol. Evol. 2016, 31, 67–80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Study area and sampling sites in upstream of Lijiang River.
Figure 1. Study area and sampling sites in upstream of Lijiang River.
Water 11 00680 g001
Figure 2. Species rarefaction curves in upstream Lijiang River. Left: Individual-based rarefaction curve; right: sample-based rarefaction curve.
Figure 2. Species rarefaction curves in upstream Lijiang River. Left: Individual-based rarefaction curve; right: sample-based rarefaction curve.
Water 11 00680 g002
Figure 3. Contributions of the environmental, hydrological, and spatial variables to the variance in species composition in upstream of Lijiang River.
Figure 3. Contributions of the environmental, hydrological, and spatial variables to the variance in species composition in upstream of Lijiang River.
Water 11 00680 g003
Figure 4. Relationship between Sørensen indices (ßSIM, ßNES, and ßSOR) and dissimilarity matrices based on Euclidean distances for spatial (Svdis), hydrological (Hvdis), and environmental (Evdis) variables in upstream Lijiang River. For other abbreviations see Table 1.
Figure 4. Relationship between Sørensen indices (ßSIM, ßNES, and ßSOR) and dissimilarity matrices based on Euclidean distances for spatial (Svdis), hydrological (Hvdis), and environmental (Evdis) variables in upstream Lijiang River. For other abbreviations see Table 1.
Water 11 00680 g004
Figure 5. Relationship between Jaccard indices (ßJTU, ßJNE, and ßJAC) and dissimilarity matrices based on Euclidean distances for spatial (Svdis), hydrological (Hvdis), and environmental (Evdis) variables in upstream Lijiang River. For other abbreviations see Table 1.
Figure 5. Relationship between Jaccard indices (ßJTU, ßJNE, and ßJAC) and dissimilarity matrices based on Euclidean distances for spatial (Svdis), hydrological (Hvdis), and environmental (Evdis) variables in upstream Lijiang River. For other abbreviations see Table 1.
Water 11 00680 g005
Figure 6. Species richness, fish abundance, and Shannon index in different stream orders in upstream Lijiang River.
Figure 6. Species richness, fish abundance, and Shannon index in different stream orders in upstream Lijiang River.
Water 11 00680 g006
Table 1. Abbreviations used in this study.
Table 1. Abbreviations used in this study.
ElementsAbbreviation
Environmental variablesEv
Hydrological variablesHv
Spatial variablesSv
Dissimilarity matrices based on Euclidean distances for environmental variablesEvdis
Dissimilarity matrices based on Euclidean distances for hydrological variablesHvdis
Dissimilarity matrices based on Euclidean distances for spatial variablesSvdis
Sørensen dissimilarityßSOR
Simpson dissimilarity (=turnover component of Sørensen dissimilarity)ßSIM
Nestedness-resultant component of Sørensen dissimilarityßNES
Jaccard dissimilarityßJAC
Turnover component of Jaccard dissimilarityßJTU
Nestedness-resultant component of Jaccard dissimilarityßJNE
Table 2. Description of hydrological (Hv), environmental (Ev), and spatial (Sv) variables in upstream Lijiang River.
Table 2. Description of hydrological (Hv), environmental (Ev), and spatial (Sv) variables in upstream Lijiang River.
Abiotic VariablesMeanMin.Max.
Environmental variables (Ev)
Water temperature (°C)22.118.825.1
pH7.86.99.7
Dissolved oxygen (mg/L)8.46.19.8
Conductivity (μS/cm)62.416.1157.4
Turbidity (NTU)3.10.218.4
Substrate size (cm)10.27.215.9
% Hydrostatic habitat (pool, plunge)28.40.076.9
% Fastflow habitat (riffle, cascade)36.37.567.0
% Slowflow habitat (run, glide) 35.40.075.0
Hydrological variables (Hv)
Wetted width (m)14.22.463.5
Depth (cm)34.316.153.5
Flow velocity (cm/s)31.214.259.4
Distance from the dam of Fuzikou Reservoir (km)20.42.235.6
Altitude (m)324.4226.6469.8
Spatial variables (Sv)
Latitude (°N)25.7825.7025.88
Longitude (°E)110.46110.50110.57
PCNM1: Principal coordinates of neighborhood matrix 10.00−0.250.32
PCNM2: Principal coordinates of neighborhood matrix 20.00−0.290.34
PCNM3: Principal coordinates of neighborhood matrix 30.00−0.380.37
PCNM4: Principal coordinates of neighborhood matrix 40.00−0.420.37
PCNM5: Principal coordinates of neighborhood matrix 50.00−0.250.48
PCNM6: Principal coordinates of neighborhood matrix 60.00−0.420.39
PCNM7: Principal coordinates of neighborhood matrix 70.00−0.620.42
PCNM8: Principal coordinates of neighborhood matrix 80.00−0.370.40
PCNM9: Principal coordinates of neighborhood matrix 90.00−0.50.36
PCNM10: Principal coordinates of neighborhood matrix 100.00−0.430.39
PCNM11: Principal coordinates of neighborhood matrix 110.00−0.610.45
PCNM12: Principal coordinates of neighborhood matrix 120.00−0.710.70
PCNM13: Principal coordinates of neighborhood matrix 130.00−0.300.56
PCNM14: Principal coordinates of neighborhood matrix 140.00−0.500.76
PCNM15: Principal coordinates of neighborhood matrix 150.00−0.400.50
PCNM16: Principal coordinates of neighborhood matrix 160.00−0.610.43
PCNM17: Principal coordinates of neighborhood matrix 170.00−0.420.82
Table 3. Fish species and abundance collected in upstream Lijiang River.
Table 3. Fish species and abundance collected in upstream Lijiang River.
Fish SpeciesFish Abundance (N)
Bagridae
1. Tachysurus adiposalis (Oshima, 1919)7
2. Tachysurus albomarginatus (Rendahl, 1928)34
Balitoridae
3. Erromyzon sinensis (Chen, 1980)312
4. Pseudogastromyzon fangi (Nichols, 1931)736
5. Vanmanenia lineata (Fang, 1935)26
6. Vanmanenia pingchowensis (Fang, 1935)86
Cichlidae
7. Oreochromis niloticus (Linnaeus, 1758) *2
Cobitidae
8. Misgurnus anguillicaudatus (Cantor, 1842)36
Cyprinidae
9. Acrossocheilus kreyenbergii (Regan, 1908)19
10. Acrossocheilus parallens (Nichols, 1931)381
11. Carassius auratus (Linnaeus, 1758)1
12. Cyprinus carpio Linnaeus, 17585
13. Microphysogobio chenhsienensis (Fang, 1938)29
14. Microphysogobio kiatingensis (Wu, 1930)3
15. Onychostoma barbatulum (Pellegrin, 1908)1
16. Onychostoma gerlachi (Peters, 1881)1
17. Opsariichthys bidens Günther, 187388
18. Parasinilabeo assimilis Wu & Yao, 197716
19. Pseudorasbora parva (Temminck & Schlegel, 1846)12
20. Rhodeus ocellatus (Kner, 1866)1
21. Squalidus atromaculatus (Nichols & Pope, 1927)27
22. Zacco platypus (Temminck & Schlegel, 1846)623
Gobiidae
23. Rhinogobius duospilus (Herre, 1935)110
24. Rhinogobius filamentosus (Wu, 1939)4
25. Rhinogobius similis Gill, 1859143
26. Rhinogobius leavelli (Herre, 1935)183
Mastacembelidae
27. Macrognathus aculeatus (Bloch, 1786) 2
Nemacheilidae
28. Oreonectes platycephalus Günther, 18682
29. Schistura fasciolata (Nichols & Pope, 1927)22
30. Schistura incerta (Nichols, 1931)9
31. Traccatichthys pulcher (Nichols & Pope, 1927)1
Percichthyidae
32. Coreoperca whiteheadi Boulenger, 190025
33. Siniperca scherzeri Steindachner, 18921
34. Siniperca undulata Fang & Chong, 19321
Siluridae
35. Pterocryptis anomala (Herre, 1933)2
36. Pterocryptis cochinchinensis (Valenciennes, 1840)1
Sisoridae
37. Glyptothorax fokiensis (Rendahl, 1925)16
Notes: * exotic species.
Table 4. Results of forward selection of environmental (Ev), hydrological (Hv), and spatial (Sv) variables for species composition in upstream Lijiang River.
Table 4. Results of forward selection of environmental (Ev), hydrological (Hv), and spatial (Sv) variables for species composition in upstream Lijiang River.
Abiotic VariablesFPAbiotic VariablesFP
Environmental variables (Ev)4.52<0.001Spatial variables (Sv)4.39<0.001
Turbidity4.550.005PCNM19.610.005
Conductivity3.930.005PCNM23.810.01
Hydrological variables (Hv)4.89<0.001PCNM32.770.025
Altitude7.930.005PCNM83.160.015
Depth1.950.03
Table 5. Results of Mantel and partial Mantel test for the correlation between ß-diversity and dissimilarity matrices based on Euclidean distances for spatial (Svdis), hydrological (Hvdis), and environmental variables (Evdis) in upstream Lijiang River. For other abbreviations: see Table 1.
Table 5. Results of Mantel and partial Mantel test for the correlation between ß-diversity and dissimilarity matrices based on Euclidean distances for spatial (Svdis), hydrological (Hvdis), and environmental variables (Evdis) in upstream Lijiang River. For other abbreviations: see Table 1.
IndexEvdisHvdisSvdisEvdis aHvdis aSvdis a
ßSIM0.269 **0.155 *0.232 ***0.188 *0.0340.128
ßSNE0.0330.150−0.0580.0470.162−0.085
ßSOR0.356 ***0.366 ***0.200 **0.281 *0.255 *0.042
ßJTU0.226 **0.1350.225 ***0.154 *0.0280.139
ßJNE−0.0100.097−0.0910.0180.126−0.099
ßJAC0.346 ***0.373 ***0.213 **0.264 *0.264 *0.061
a The pure effect while controlling for the other two distances. * P < 0.05, ** P < 0.01, *** P < 0.001

Share and Cite

MDPI and ACS Style

Huang, L.; Huang, J.; Wu, Z.; Mo, Y.; Zou, Q.; Jeppesen, E.; Wu, N. Beta Diversity Partitioning and Drivers of Variations in Fish Assemblages in a Headwater Stream: Lijiang River, China. Water 2019, 11, 680. https://doi.org/10.3390/w11040680

AMA Style

Huang L, Huang J, Wu Z, Mo Y, Zou Q, Jeppesen E, Wu N. Beta Diversity Partitioning and Drivers of Variations in Fish Assemblages in a Headwater Stream: Lijiang River, China. Water. 2019; 11(4):680. https://doi.org/10.3390/w11040680

Chicago/Turabian Style

Huang, Liangliang, Jian Huang, Zhiqiang Wu, Yuanmin Mo, Qi Zou, Erik Jeppesen, and Naicheng Wu. 2019. "Beta Diversity Partitioning and Drivers of Variations in Fish Assemblages in a Headwater Stream: Lijiang River, China" Water 11, no. 4: 680. https://doi.org/10.3390/w11040680

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