Next Article in Journal
Detection of SARS-CoV-2 and Variants in Hospital Wastewater in a Developing Country
Next Article in Special Issue
Niche Analysis of the Main Fish in the Lhasa River Basin
Previous Article in Journal
Unraveling Physical and Chemical Effects of Textile Microfibers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Different Land Use Types on the Taxonomic and Functional Diversity of Macroinvertebrates in an Urban Area of Northern China

Liaoning Provincial Key Laboratory for Hydrobiology, College of Fisheries and Life Science, Dalian Ocean University, Dalian 116000, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(23), 3793; https://doi.org/10.3390/w14233793
Submission received: 19 October 2022 / Revised: 15 November 2022 / Accepted: 15 November 2022 / Published: 22 November 2022
(This article belongs to the Special Issue Ecology of Freshwater Fishes)

Abstract

:
The urbanization of riverine landscapes is an increasing threat to river ecosystems. However, it is unclear which metrics can best assess the response of macroinvertebrates to the conversion of forested lands to urban and agricultural lands. The main goal of this study is to examine whether trait-based approaches are more sensitive than taxonomic approaches in distinguishing macroinvertebrate responses to different land use types in a highly urbanized area of northern China. Results based on 14 environmental variables showed a significant difference across a human-induced environmental gradient. The results showed that no significant differences were observed in terms of taxonomic diversity indices between the different land use types. Functional evenness (FEve) and Rao’s quadratic entropy decreased with the increase in urban area caused by the intensification of human activity, demonstrating that functional diversity is more sensitive than taxonomic diversity in discriminating between different land use types. In addition, the results based on RLQ (physical–chemical variables (R), macroinvertebrate taxa (L), and species traits (Q)) and fourth-corner analyses indicated that the trait states of bi- or multivoltine, high dispersal capacity, and not-streamlined body shape were much higher in the agricultural area and positively related to farmland percentage. Taxa with large body size were dominant in urban areas and were positively correlated with EC. Overall, the observed responses of traits to environmental variables suggest that trait-based approaches should be incorporated into land use management for river restoration.

1. Introduction

The ecological impact of land use change on river and stream health is an important environmental issue worldwide Allan [1,2]. Increasing evidence suggests that changes in community composition caused by land use have been associated with the homogenization of taxa and the loss of taxonomic diversity [3]. For example, increased agricultural and urban land use in watersheds leads to the discharge of large amounts of domestic sewage and some industrial wastewater into urban rivers [4,5]. The abundance of pollution-tolerant taxa increases dramatically, and sensitive taxa are significantly reduced or disappear [6,7]. In addition, the expansion of agricultural land use within a watershed alters the nutrient dynamics due to the use of chemical fertilizers, and excess nitrogen and phosphorus inputs have modified community structure and biodiversity [8]. In this context, ecologists aim to understand the effects of land use change on biodiversity and explore which assessment methods can accurately examine the responses.
Macroinvertebrates are valuable and excellent indicators for biomonitoring and bioassessment [9]; compared to other aquatic indicators, macroinvertebrates have many advantages, such as sensitivity to stream disturbance, short life cycles, limited mobility, and ease of collection and identification [10]. Macroinvertebrates can also indicate habitat degradation and fragmentation caused by land use change in the context of increasing concentrations of agricultural and urban pollutants [1,11,12].
The response of macroinvertebrates to the conversion of forest land to farmland and urban land is usually examined from the perspectives of taxonomic diversity [13]. For instance, changes in land use types can significantly alter community composition within a region by altering the physicochemical characteristics of the aquatic ecosystem [7]. The taxonomic diversity of macroinvertebrates significantly decreases along land use intensity and urbanization gradients [14]. However, taxonomy-based tools typically focus only on biodiversity measurements, ignoring the relationships between the environmental characteristics and the traits of macroinvertebrates in the ecosystem [15,16]. In aquatic ecosystems, trait-based methods have been confirmed to have certain advantages in assessing the impact of different land uses on macroinvertebrates [17,18,19]. The biological trait profile of communities offers an alternative approach to traditional measures of macroinvertebrate taxonomic identity and is less constrained by biogeographic influences [16]. Across broader regional extents and natural environmental gradients, the responses of certain functional traits to land use changes at various spatial scales were consistent because most stressors affect only certain trait categories [20]. Hence, in this context, assessing the response of functional traits to different land use types can not only distinguish multiple stressors but also provide early warning signals of ecosystem processes responses to such impacts before actual species loss [21,22,23].
This study aimed to assess the differences in taxonomic and functional diversity of macroinvertebrates in three land use types—forest land, farmland, and urban land—and the reasons for such differences during urbanization. Jinan, as a highly urbanized city, has gradually transformed most of its forest land into agricultural and urban land through disturbance from frequent human activities. The biodiversity of river ecosystems has decreased, and the functional composition has changed. Therefore, Jinan can serve as a typical city for studying the effects of land use change on urban river biodiversity and provide a theoretical basis for biodiversity conservation in urban rivers.
Three hypotheses were tested in our study. First, increasingly frequent human activities have significantly altered species composition in aquatic ecosystems. We hypothesized that species composition varies by land use type, with forested lands having a richer species composition and higher taxonomic diversity [24,25]. Second, some traits are more sensitive to land use change, and stressors associated with the shift from forest to urbanization and agricultural land can act as environmental filters that allow for differences in macroinvertebrate trait composition across land use types [26,27]. Third, trait-based approaches are more sensitive to environmental changes and can provide early warning signals of land use change before actual species loss [19,28]; therefore, rather than taxonomic diversity approaches, functional trait-based approaches can better distinguish differences in macroinvertebrates between land use types.

2. Materials and Methods

2.1. Study Area

This study was conducted in the regions of Jinan. Jinan is the capital of Shandong Province, located in the northern region of China. It has a short springs and autumns; hot and rainy summers; and long, cold winters. Based on the characteristics of land use types, the city of Jinan was categorized into mountain–hilly (Table S2), urban (Table S3), and agricultural land (Table S4) [29]. The southern area has a high degree of forest cover and is classified as a mountain–hilly area, and the primary flowing water system is the Yellow River. The middle region is more urbanized than the other areas and is therefore defined as an urban area, with the Xiaoqing River flowing through it. The northern region is a traditional agricultural area, and the Tuhai River, which is part of the Haihe River system, flows through the region as a large drainage river.
In 2015, the population of permanent residents in Jinan was 6.8 million, which consisted of 4.1 million individuals living in urban areas and 2.7 million individuals living in rural areas. With the increase in population, land use and cover have dramatically changed [30] in this region over the past decades, reflecting the main types and degrees of anthropogenic disturbance in urban ecosystems. Repeat sampling was conducted at 44 sites during the spring, summer, and fall of each year from 2014–2016. These sites were divided into 15 mountain–hilly, 15 urban, and 14 agricultural sites (Figure 1).

2.2. Data Collection

Land use analysis was conducted using digital elevation models (DEM) with 30 m resolution in 2015, provided by the Chinese Academy of Sciences (http://www.cnic.cn/, accessed on 27 February 2022), to evaluate land use in the riparian zone around these sites. The Tabulate Area function in ArcGIS 10.5 (Esri, Inc, Redlands, CA, USA) was used to calculate the percentage of farmland, forestland, grassland, water area, and urban area within a 300 m buffer upstream of each sampling site to illustrate changes in the natural habitat structure of the river under different land use practices [31].
Eleven traits with thirty-six different trait categories were considered (Table 1 and Table S5). Traits included those related to life history (voltinism and development), mobility (dispersal and drift), morphology (e.g., attachment, armoring, shape, respiration, and size), and ecology (rheophily, habit, and trophic habit). Trait categories were selected based on their mechanistic relationships with urban pollution. Trait availability and ease of measurement were taken into account during selection. There is considerable evidence in the literature that selected traits markedly respond to anthropogenic stress [32,33,34,35,36,37]. Information on trait categories identified was obtained mainly from the literature [38,39,40]. When trait information is not available at the genus or species taxonomic levels, we typically retrieve it at the family level instead. Affinity scores were used in a fuzzy coding approach that considers the variability of macroinvertebrate life histories and better accounts for differences between species within a family and between species at different growth stages [41,42]. We use an affinity score range of 0–5 to describe the affinity of an organism for a trait attribute, with a score of 0 indicating that the organism has no affinity for a particular trait attribute and a score of 5 indicating a maximum affinity for that trait attribute. A score of 1 was used to indicate low affinity and a score of 3 to indicate moderate affinity [43]. Potential responses of macroinvertebrate community traits to land use type changes and associated stressors were predicted a priori.
The measured physicochemical variables included water temperature, pH, electrical conductivity (EC), dissolved oxygen (DO) concentration, total nitrogen (TN), ammonium nitrogen (NH3-N), nitrate nitrogen, chemical oxygen demand (COD), biochemical oxygen demand (BOD), and total phosphorus (TP). Water temperature (Temp), DO, and EC were determined using a water quality analyzer (YSI85). After sampling, all water samples were placed in low-temperature incubators for storage and brought back to the laboratory. Physicochemical variables, such as COD, TN, NH3-N, and TP, were measured following the Chinese surface water quality standards [45]. Macroinvertebrates were collected from 2014 to 2016 in the spring, autumn, and summer using a D-frame kick net (30 × 30 cm aperture) of mesh size 500 μm. Sampling was repeated at three locations within the same riffle. All samples were transferred through a 60-mesh sieve into 300 mL plastic bottles and stored in 90% alcohol solution. After the samples were brought back to the laboratory, the macroinvertebrates were separated from the impurities, such as rocks and debris. The benthic samples were transferred into 100 mL plastic bottles using the manual picking method, stored in 95% alcohol solution, and sorted and counted under a microscope or dissecting microscope. Macroinvertebrates were identified at the genus level or at the lowest practical taxonomic level in the laboratory using a compound light microscope and local guidelines for freshwater invertebrate identification [46,47,48,49,50,51].

2.3. Data Analysis

Taxonomic diversity across different land use types was compared using the Shannon diversity index, Pielou’s diversity index, and species richness. Functional diversity analyses were conducted using the relative abundance of each trait category in the FD package to calculate the functional evenness (FEve), functional richness (FRic), functional divergence (FDiv), functional dispersion (FDis), and Rao’s quadratic entropy (RaoQ) [52]. Functional richness (FRic) measures the size of the trait space occupied by the community [53]. Functional evenness (FEve) measures the homogeneity of the distribution of functional indices in the trait space, and functional divergence (FDiv) measures the proportion of species abundance at the edge of the trait space [54]. Functional dispersion (FDis) uses the relative abundance of species traits to calculate the dispersion of traits. Rao’s quadratic entropy quantifies the divergent aspects of functional diversity as it not only takes into account differences in traits but also involves differences in abundance between species [55,56]. These indices can only reflect one aspect of functional diversity each, and combining them would provide a more accurate description [57].
R-environmental characteristics of sampling sites, L-species abundance across samples, Q-species traits (RLQ), and fourth-corner methods were used to identify the bivariate relationships between environmental variables and functional traits [37]. RLQ analysis maximizes the covariance between traits and environmental variables mediated by species abundance [16,58]. The fourth-corner approach quantifies and testes all correlations between trait categories and environmental variables. Three matrices were developed for the RLQ analysis, including a matrix of environmental variables, a matrix of functional traits, and a matrix of taxa abundance. Correspondence ordination analysis (CA) was performed on the taxa abundance matrices, principal component analysis (PCA) was applied to the functional trait matrices, and the Hill–Smith function was performed on the environmental variables [58]. The RLQ function was designed to integrate the results of the separate analyses performed on the R, L, and Q matrices. Monte Carlo tests were used to test the overall significance of Model 2 (H1: assumes no relationship between R and L) and Model 4 (H2: assumes no relationship between L and Q) [59]. According to the results of the RLQ ordination, the fourth-corner test was used to elucidate the relationship between functional traits and environmental variables to further confirm the ecological preferences of different land use types. The p values for the fourth-corner analysis were adjusted using the false discovery rate method.
The FD package was used to calculate the functional diversity indices [60], and the ade4 package was used for the RLQ and fourth-corner methods [61,62]. All statistical analyses were performed in R version 4.2.0 (R Development Core Team 2022). The Kruskal-Wallis test was performed using Origin 2022 (Origin Lab, Hampton, CA, USA).

3. Results

3.1. Environmental Variables

The results showed highly significant differences in physicochemical variables between the different land use types (Kruskal-Wallis test: p < 0.05, Table 2). Among these variables, the pH value, DO concentration, and percentage of forestland in mountain–hilly habitat were the highest. In contrast, the EC, ammonium nitrogen, TN, TP, COD, % farmland, and % urban land were markedly higher in the urban and agricultural areas than in mountain–hilly areas.

3.2. Taxonomic and Functional Diversity

In total, 73 taxa were collected during the 2014–2016 sampling period, belonging to 12 orders, 25 families, and 30 genera. All taxonomic diversity indices used, namely species richness, the Shannon-Wiener index, Simpson diversity, and species evenness, did not differ markedly between different site groups (Figure 2).
Functional evenness (FEve) and Rao’s quadratic entropy (RaoQ) were significantly higher in mountainous hilly areas than in agricultural areas (Figure 3). The functional diversity indices FRic, FDis, and FDiv exhibited no marked difference between the different land use types.
In terms of the functional traits, 23 trait categories showed significant differences across land use types (Kruskal-Wallis test: p < 0.05, Table 1, Figure 4, Table S1). Non-seasonal taxa dominated at agricultural sites, while fast-seasonal taxa showed high abundance at mountain–hilly sites. Regarding “dispersal”, urban and agricultural areas had a much higher abundance of individuals with high dispersal ability. In terms of “drift”, taxa with drifting ability were abundant at mountain–hilly sites. In addition, mountain–hilly areas exhibited a significantly higher abundance of taxa characterized by the traits “free-ranging” and “soft-bodied”. In the case of “size” traits, small individuals predominated in mountain–hilly areas, whereas medium-sized taxa were most abundant at agricultural sites. In addition, the mountain–hilly area was dominated by burrowers, and the clingers were most abundant in the agricultural area. Predators were most prevalent at the agricultural sites, while mountain–hilly areas had a higher abundance of collector–gatherers. No significant difference was observed between the three land use types for all trait categories of body shape, respiration, and rheophily.

3.3. Relationship between Functional Diversity and Environmental Variables

The first two axes accounted for 91.9% (RLQ1 = 78.5%, RLQ2 = 13.4%) of the total variance between functional traits and environmental variables across the whole ecoregion (Table 3). The results of the fourth-corner analysis highlighted the bivariate association between environmental variables and functional traits (Figure 5). Among these variables, five trait categories showed predictive responses to environmental land use types, such as bi- or multivoltine, high dispersal capability, a rare occurrence of drift, not-streamlined body shape, and depositional habitat, which showed considerably positive correlations with the percent of farmland. The following traits of preference for free-ranging, soft-bodied forms and collector-gathering were positively correlated with total nitrogen. The composition of the five trait categories responded predictively with the TP concentration. For instance, organisms with a fast seasonal life cycle, a low dispersal capability, an abundant occurrence in drift, a small body size, and a burrowing habit were dominated by a high disturbance gradient with high phosphorus concentrations. Meanwhile, these traits were also negatively correlated with DO.

4. Discussion

Our results indicated that there were significant differences in the composition of macroinvertebrate traits between land use types in urban ecoregions. An interesting result showed that individuals characterized by traits related to resistance and resilience (i.e., small body size) occurred more frequently in mountain-hilly sites. The taxonomic diversity index in the present study was far less sensitive than the functional diversity index in distinguishing different land use types.

4.1. Response of Taxonomic Diversity to Different Types of Land Use

Inconsistent with previous studies [63,64,65], our results showed no significant differences in species composition and taxonomic diversity between land use types. The results indicated that taxonomic methods were not powerful tools for monitoring and evaluating the environmental impacts of land use changes on macroinvertebrates in highly urbanized ecoregions. This finding was in accordance with our third hypothesis that the taxonomic diversity index demonstrated low predictive capability for different land use types compared to the functional diversity index.

4.2. Response of Functional Diversity to Different Types of Land Use

As hypothesized, functional diversity showed a significant difference along the disturbance gradient, which could be a valuable indicator for distinguishing different land use types. Mountain–hilly areas exhibited higher FEve values than agricultural areas. Higher FEve values indicated that the distribution of trait abundance was relatively evenly disturbed in the functional space. Some researchers have reported relevant studies in which functional evenness (FEve) decreased with increasing disturbance, which was consistent with our findings [66,67]. RaoQ represents a mix between FRic and FDiv, which was demonstrated to be a promising indicator of disturbance [68]. The RaoQ values were significantly higher in mountain-hilly areas than in agricultural areas, indicating that land use changes caused by human activities, especially agricultural activities, lead to a decrease in water quality and habitat degradation in ecosystems [69]. Frequent anthropogenic stressors cause community variation and screening of biological trait composition, ultimately leading to changes in the functional diversity of aquatic insect communities.
In terms of the trait of large body size, we predicted that large body size would be less associated with highly disturbed sites and that small body size would be more prevalent as urban pollution increased. However, our results indicated the opposite trend [37]. The present study shows that taxa with large body size, such as gastropods, were dominant in urban areas. In addition, small-bodied species, such as chironomid species, were the most dominant taxa in mountain–hill areas. Our results support that those who reported large body size were associated with heavily impacted urban rivers [70,71,72]. In our study, electrical conductivity is an important variable influencing the pattern of body size distribution because it affects the osmotic balance. Macroinvertebrates are exposed to contaminants through external contact, so individuals with larger body size can minimize external exposure to dissolved salts by reducing the body surface area to volume ratio. Thus, high electrical conductivity in agricultural and urban areas contributed to the high relative abundance of the large-bodied species [19]. Small-bodied chironomids were associated with traits such as a fast seasonal life cycle, low dispersal ability, high drift occurrence, and burrowing. Previous studies have shown that these traits were highly resilient and resistant in the face of environmental disturbances, reducing the time to recover from disturbances [21]. This conclusion is supported by most studies on natural river disturbances [36,37,73]. However, in our study, rivers under different land use types were subjected to different levels of pollution from human activities, resulting in low species richness and mostly pollution-tolerant species. Therefore, our results do not support using these traits as indicators of urban river pollution levels. The trait states of lack of free-ranging, no body armor, and collector–gatherer feeding mode were associated with mountain–hilly sites and positively correlated with increasing TN. Burrowing species preferred soft-bottom sediment habitats, thus coping with environmental pollution by using sediment as a buffer mechanism in water bodies [74]. The dominance of burrowers could be explained by the increased deposition of fine sediments in streambeds due to anthropogenic stressors in mountain–hilly areas, which reduces sediment heterogeneity [75]. The dispersal ability of aquatic dispersers in mountain–hilly areas was weak, which may be caused by the poor connectivity of aquatic systems and habitat fragmentation [76,77]. A higher proportion of gatherers was observed in the mountain–hilly areas and was significantly positively correlated with the TP and COD. Gatherers dominated the habitat of the mountain–hilly area, reflecting the fairly high levels of fine particulate organic matter in both water bodies and sediment [78]. The absolute predominance of gatherers may be partial because the main filterers are chironomids, which is inconsistent with the research conclusions of other scholars [79].
The percentage of farmland was positively correlated with a preference for bivoltine or multivoltine, high dispersal ability, a rare occurrence in drift, non-streamlined body shape, and depositional habitats. Previous studies have shown that taxa with traits of bi- or multivoltine mainly occur in disturbed sites because they have strong recolonization and refuge exploitation abilities in perturbed environments. For instance, studies [36] reported that the presence of organisms that exhibit bivoltine to pesticides impacted rivers, which was consistent with our research. The number of organisms with rare occurrences in drift increased with ongoing habitat destruction and degradation, mainly because the loss of boulders reduced habitat availability for taxa. The accumulation of sediments and the creation of depositional environments [80] for taxa prevented them from leaving stressed environments for more selective conditions. The trait of a non-streamlined body shape was also closely related to environmental stressors [81]. It has been reported previously that organisms with a streamlined body shape can reduce the drag effects caused by water currents. In the case of trophic habit, herbivores were highly abundant in agricultural areas, showing a predictable response to the anthropogenic disturbance gradient. The suspended particulate organic matter content was relatively high because of the large-scale use of fertilizer, herbicides, and pesticides, which contributed to the enhancement of aquatic plant growth, such as algae, and the increasing proportion of herbivores [70].
Respiratory traits indicated the response of the taxa to DO in the water column, which reflected DO concentration changes. Taxa breathing with gills were dominant under high DO conditions [82], while species with aerial respiration were dominant at sites with low DO and degraded ecosystem habitats. Our study found that the DO concentration was higher across the three land use types and that there was a higher abundance of taxa with gill respiration. The finding also supports the general idea that in high-DO environments, organisms with gill respiration are more common [83].

5. Conclusions

In conclusion, a comprehensive understanding of the differences in community characteristics and trait combinations under different land use types can contribute to the prediction of potential changes in aquatic community structure and ecosystem functioning in urban rivers. The results showed that the functional approaches of macroinvertebrates were better than taxonomic approaches for evaluating the ecological quality of urban river ecosystems. In addition, environmental variables played a dominant role in determining community trait composition. Therefore, this study has important implications for freshwater bioassessment in urban areas, especially in other areas with similar land use types.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w14233793/s1, Table S1. Kruskal-Wallis test results of the community-weighted means of trait values (CWM) of each trait category (mean ± SD) distributed across different types of land-use. Significant differences of environmental variables across different types of land-use sites are indicated in bold (Kruskal-Wallis Test: p < 0.05). Table S2. The proportions of land cover surrounding each site in mountain-hilly area. Table S3. The proportions of land cover surrounding each site in urban area. Table S4. The proportions of land cover surrounding each site in agricultural area. Table S5. Species distribution under the three analyzed land use types. Trait abbreviations and descriptions of the are provided in Table 1.

Author Contributions

Conceptualization, X.Y. and A.L.; methodology, A.L.; software, A.L., J.L. and B.Z.; formal analysis, A.L. and J.L.; data curation, A.L.; writing—original draft preparation, A.L.; writing—review and editing, X.Y.; project administration, X.Y.; funding acquisition, X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (41977193), and the National Science and Technology Basic Resources Survey Program of China (2019FY101700).

Data Availability Statement

Data are contained within the article or Supplementary Material.

Acknowledgments

Thanks are due to Jiaxin Li for assistance with the experiments and to Biao Zheng for valuable discussion.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Allan, J.D. Landscapes and riverscapes: The influence of land use on stream ecosystems. Annu. Rev. Ecol. Evol. Syst. 2004, 35, 257–284. [Google Scholar] [CrossRef] [Green Version]
  2. Gyawali, S.; Techato, K.; Monprapussorn, S.; Yuangyai, C. Integrating Land use and Water quality for Environmental based land use planning for U-tapao River Basin, Thailand. In Proceedings of the PSU-USM International Conference on Humanities and Social Sciences, Hat Yai, Thailand, 3–4 October 2011; pp. 556–563. [Google Scholar]
  3. Barnum, T.R.; Weller, D.E.; Williams, M. Urbanization reduces and homogenizes trait diversity in stream macroinvertebrate communities. Ecol. Appl. 2017, 27, 2428–2442. [Google Scholar] [CrossRef] [PubMed]
  4. Rueda, J.; Camacho, A.; Mezquita, F.; Hernandez, R.; Roca, J.R. Effect of episodic and regular sewage discharges on the water chemistry and macroinvertebrate fauna of a Mediterranean stream. Water Air Soil Pollut. 2002, 140, 425–444. [Google Scholar] [CrossRef]
  5. Covich, A.P.; Palmer, M.A.; Crowl, T.A. The role of benthic invertebrate species in freshwater ecosystems—Zoobenthic species influence energy flows and nutrient cycling. Bioscience 1999, 49, 119–127. [Google Scholar] [CrossRef] [Green Version]
  6. Piano, E.; Doretto, A.; Mammola, S.; Falasco, E.; Fenoglio, S.; Bona, F. Taxonomic and functional homogenisation of macroinvertebrate communities in recently intermittent Alpine watercourses. Freshw. Biol. 2020, 65, 2096–2107. [Google Scholar] [CrossRef]
  7. Fierro, P.; Bertrán, C.; Tapia, J.; Hauenstein, E.; Peña-Cortés, F.; Vergara, C.; Cerna, C.; Vargas-Chacoff, L. Effects of local land-use on riparian vegetation, water quality, and the functional organization of macroinvertebrate assemblages. Sci. Total Environ. 2017, 609, 724–734. [Google Scholar] [CrossRef]
  8. Townsend, C.R.; Uhlmann, S.S.; Matthaei, C.D. Individual and combined responses of stream ecosystems to multiple stressors. J. Appl. Ecol. 2008, 45, 1810–1819. [Google Scholar] [CrossRef]
  9. Bonada, N.; Prat, N.; Resh, V.H.; Statzner, B. Developments in aquatic insect biomonitoring: A comparative analysis of recent approaches. Annu. Rev. Entomol. 2006, 51, 495–523. [Google Scholar] [CrossRef] [Green Version]
  10. Serrano Balderas, E.C.; Grac, C.; Berti-Equille, L.; Armienta Hernandez, M.A. Potential application of macroinvertebrates indices in bioassessment of Mexican streams. Ecol. Indic. 2016, 61, 558–567. [Google Scholar] [CrossRef]
  11. Cortelezzi, A.; Sierra, M.V.; Gomez, N.; Marinelli, C.; Capitulo, A.R. Macrophytes, epipelic biofilm, and invertebrates as biotic indicators of physical habitat degradation of lowland streams (Argentina). Environ. Monit. Assess. 2013, 185, 5801–5815. [Google Scholar] [CrossRef]
  12. Blann, K.L.; Anderson, J.L.; Sands, G.R.; Vondracek, B. Effects of Agricultural Drainage on Aquatic Ecosystems: A Review. Crit. Rev. Environ. Sci. Technol. 2009, 39, 909–1001. [Google Scholar] [CrossRef]
  13. Matomela, N.; Chakona, A.; Kadye, W. Comparative assessment of macroinvertebrate communities within three Afromontane headwater streams influenced by different land use patterns. Ecol. Indic. 2021, 129, 107972. [Google Scholar] [CrossRef]
  14. Liu, Z.; Li, Z.; Castro, D.M.P.; Tan, X.; Jiang, X.; Meng, X.; Ge, Y.; Xie, Z. Effects of different types of land-use on taxonomic and functional diversity of benthic macroinvertebrates in a subtropical river network. Environ. Sci. Pollut. Res. 2021, 28, 44339–44353. [Google Scholar] [CrossRef] [PubMed]
  15. Carter, J.L.; Resh, V.H. After site selection and before data analysis: Sampling, sorting, and laboratory procedures used in stream benthic macroinvertebrate monitoring programs by USA state agencies. J. N. Am. Benthol. Soc. 2001, 20, 658–682. [Google Scholar] [CrossRef] [Green Version]
  16. Menezes, S.; Baird, D.J.; Soares, A.M.V.M. Beyond taxonomy: A review of macroinvertebrate trait-based community descriptors as tools for freshwater biomonitoring. J. Appl. Ecol. 2010, 47, 711–719. [Google Scholar] [CrossRef]
  17. Gayraud, S.; Statzner, B.; Bady, P.; Haybachp, A.; Schöll, F.; Usseglio-Polatera, P.; Bacchi, M. Invertebrate traits for the biomonitoring of large European rivers: An initial assessment of alternative metrics. Freshw. Biol. 2003, 48, 2045–2064. [Google Scholar] [CrossRef] [Green Version]
  18. Bady, P.; Dolédec, S.; Fesl, C.; Gayraud, S.; Bacchi, M.; Schöll, F. Use of invertebrate traits for the biomonitoring of European large rivers: The effects of sampling effort on genus richness and functional diversity. Freshw. Biol. 2005, 50, 159–173. [Google Scholar] [CrossRef]
  19. Doledec, S.; Statzner, B. Invertebrate traits for the biomonitoring of large European rivers: An assessment of specific types of human impact. Freshw. Biol. 2008, 53, 617–634. [Google Scholar] [CrossRef]
  20. Dolédec, S.; Phillips, N.; Townsend, C. Invertebrate community responses to land use at a broad spatial scale: Trait and taxonomic measures compared in New Zealand rivers. Freshw. Biol. 2011, 56, 1670–1688. [Google Scholar] [CrossRef]
  21. Townsend, C.R.; Hildrew, A.G. Species traits in relation to a habitat templet for river systems. Freshw. Biol. 1994, 31, 265–275. [Google Scholar] [CrossRef]
  22. Zuellig, R.E.; Schmidt, T.S. Characterizing invertebrate traits in wadeable streams of the contiguous US: Differences among ecoregions and land uses. Freshw. Sci. 2012, 31, 1042–1056. [Google Scholar] [CrossRef] [Green Version]
  23. Alahuhta, J.; Erős, T.; Kärnä, O.-M.; Soininen, J.; Wang, J.; Heino, J. Understanding environmental change through the lens of trait-based, functional, and phylogenetic biodiversity in freshwater ecosystems. Environ. Rev. 2018, 27, 263–273. [Google Scholar] [CrossRef] [Green Version]
  24. Du, X.; Song, D.; Ming, K.; Jin, X.; Wang, H.B.; Wang, L.; Liu, H.; Zhao, C.; Huo, T.B. Response of macroinvertebrate communities to land use and water quality in Wudalianchi Lake. Ecol. Evol. 2021, 11, 1368–1377. [Google Scholar] [CrossRef] [PubMed]
  25. Castro-Lopez, D.; Rodriguez-Lozano, P.; Arias-Real, R.; Guerra-Cobian, V.; Prat, N. The Influence of Riparian Corridor Land Use on the Pesqueria River’s Macroinvertebrate Community (NE Mexico). Water 2019, 11, 1930. [Google Scholar] [CrossRef] [Green Version]
  26. de Castro, D.M.P.; Doledec, S.; Callisto, M. Land cover disturbance homogenizes aquatic insect functional structure in neotropical savanna streams. Ecol. Indic. 2018, 84, 573–582. [Google Scholar] [CrossRef]
  27. Firmiano, K.R.; Castro, D.M.P.; Linares, M.S.; Callisto, M. Functional responses of aquatic invertebrates to anthropogenic stressors in riparian zones of Neotropical savanna streams. Sci. Total Environ. 2021, 753, 141865. [Google Scholar] [CrossRef]
  28. Leitao, R.P.; Zuanon, J.; Mouillot, D.; Leal, C.G.; Hughes, R.M.; Kaufmann, P.R.; Villeger, S.; Pompeu, P.S.; Kasper, D.; de Paula, F.R.; et al. Disentangling the pathways of land use impacts on the functional structure of fish assemblages in Amazon streams. Ecography 2018, 41, 219–232. [Google Scholar] [CrossRef]
  29. Liu, S.; Dong, X.; Zhou, X.; Chen, S. Characteristic analysis of weather and climate memory in temperature: Taking Jinan and Qingdao as examples. J. Mar. Meteorol. 2022, 42, 107–114. [Google Scholar] [CrossRef]
  30. Li, H.; Wang, S.; Ji, G.; Zhang, L. Changes in land use and ecosystem service values in Jinan, China. Energy Procedia 2011, 5, 1109–1115. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, J.; Jiang, X.; Li, Z.; Meng, X.; Heino, J.; Xie, Z.; Wang, X.; Yu, J. Changes in multiple facets of macroinvertebrate alpha diversity are linked to afforestation in a subtropical riverine natural reserve. Environ. Sci. Pollut. Res. 2018, 25, 36124–36135. [Google Scholar] [CrossRef]
  32. Poff, N.L.; Olden, J.D.; Vieira, N.K.M.; Finn, D.S.; Simmons, M.P.; Kondratieff, B.C. Functional trait niches of North American lotic insects: Traits-based ecological applications in light of phylogenetic relationships. J. N. Am. Benthol. Soc. 2006, 25, 730–755. [Google Scholar] [CrossRef] [Green Version]
  33. Edegbene, A.O.; Arimoro, F.O.; Odume, O.N. How does urban pollution influence macroinvertebrate traits in forested riverine systems? Water 2020, 12, 3111. [Google Scholar] [CrossRef]
  34. Erasmus, J.H.; Lorenz, A.W.; Zimmermann, S.; Wepener, V.; Sures, B.; Smit, N.J.; Malherbe, W. A diversity and functional approach to evaluate the macroinvertebrate responses to multiple stressors in a small subtropical austral river. Ecol. Indic. 2021, 131, 108206. [Google Scholar] [CrossRef]
  35. Liu, Y.; Tian, Y.; Gao, Y.; Cui, D.; Zhang, W.; Jiao, Z.; Yao, F.; Zhang, Z.; Yang, H. The Impacts of Different Anthropogenic Disturbances on Macroinvertebrate Community Structure and Functional Traits of Glacier-Fed Streams in the Tianshan Mountains. Water 2022, 14, 1298. [Google Scholar] [CrossRef]
  36. Kuzmanovic, M.; Dolédec, S.; de Castro-Catala, N.; Ginebreda, A.; Sabater, S.; Muñoz, I.; Barceló, D. Environmental stressors as a driver of the trait composition of benthic macroinvertebrate assemblages in polluted Iberian rivers. Environ. Res. 2017, 156, 485–493. [Google Scholar] [CrossRef] [PubMed]
  37. Diaz, A.M.; Alonso, M.L.S.; Gutierrez, M. Biological traits of stream macroinvertebrates from a semi-arid catchment: Patterns along complex environmental gradients. Freshw. Biol. 2008, 53, 1–21. [Google Scholar] [CrossRef]
  38. Edegbene, A.O.; Arimoro, F.O.; Odume, O.N. Exploring the distribution patterns of macroinvertebrate signature traits and ecological preferences and their responses to urban and agricultural pollution in selected rivers in the Niger Delta ecoregion, Nigeria. Aquat. Ecol. 2020, 54, 553–573. [Google Scholar] [CrossRef]
  39. Palmer, C.G.; Maart, B.; Palmer, A.R.; O’Keeffe, J.H. An assessment of macroinvertebrate functional feeding groups as water quality indicators in the Buffalo River, eastern Cape Province, South Africa. Hydrobiologia 1996, 318, 153–164. [Google Scholar] [CrossRef]
  40. Usseglio-Polatera, P.; Bournaud, M.; Richoux, P.; Tachet, H. Biological and ecological traits of benthic freshwater macroinvertebrates: Relationships and definition of groups with similar traits. Freshw. Biol. 2000, 43, 175–205. [Google Scholar] [CrossRef]
  41. Mondy, C.P.; Usseglio-Polatera, P. Using fuzzy-coded traits to elucidate the non-random role of anthropogenic stress in the functional homogenisation of invertebrate assemblages. Freshw. Biol. 2014, 59, 584–600. [Google Scholar] [CrossRef]
  42. Chevene, F.; DolÉAdec, S.; Chessel, D. A fuzzy coding approach for the analysis of long-term ecological data. Freshw. Biol. 1994, 31, 295–309. [Google Scholar] [CrossRef]
  43. Frank, C.A. Exploring the community structure of Afrotropical macroinvertebrate traits and ecological preferences along an agricultural pollution gradient in the Kat River, Eastern Cape, South Africa. Ecol. Indic. 2022, 135, 108570. [Google Scholar] [CrossRef]
  44. Dolédec, S.; Statzner, B.; Bournard, M. Species traits for future biomonitoring across ecoregions: Patterns along a human-impacted river. Freshw. Biol. 1999, 42, 737–758. [Google Scholar] [CrossRef]
  45. Huang, X.; Chen, W.; Cai, Q. Survey, observation and analysis of lake ecology. In Standard Methods for Observation and Analysis in Chinese Ecosystem Research Network, Series V; Standards Press of China: Beijing, China, 1999. [Google Scholar]
  46. Morse, J.C.; Yang, L.; Tian, L. Aquatic Insects of China Useful for Monitoring Water Quality; Hohai University Press: Nanjing, China, 1984. [Google Scholar]
  47. Kawai, S. Aquatic Insects of Japan: Manual with Keys and Illustrations; Tokai University Press: Tokyo, Japan, 1925. [Google Scholar]
  48. Brinkhurst, R.O. Guide to the Freshwater Aquatic Microdrile Oligochaetes of North America; Department of Fisheries and Oceans: Vancouver, BC, Canada, 1986. [Google Scholar]
  49. Epler, J.H. Identification Manual for the Larval Chironomidae (Diptera) of Florida; Bureau of Surface Water Management, Florida Department of Environmental Protection: Tallahassee, FL, USA, 1995. [Google Scholar]
  50. Merritt, R.W.; Cummins, K.W. An Introduction to the Aquatic Insects of North America; Kendall Hunt: Dubuque, IA, USA, 1996. [Google Scholar]
  51. Wiggins, G.B. Key to Larvae of North American Families of Trichoptera; University of Toronto Press: Toronto, ON, Canada, 1996; pp. 43–48. [Google Scholar]
  52. Laliberte, E.; Legendre, P. A distance-based framework for measuring functional diversity from multiple traits. Ecology 2010, 91, 299–305. [Google Scholar] [CrossRef]
  53. Cornwell, W.; Schwilk, D.; Ackerly, D. A trait-based test for habitat filtering: Convex hull volume. Ecology 2006, 87, 1465–1471. [Google Scholar] [CrossRef]
  54. Villéger, S.; Mason, N.W.H.; Mouillot, D. New multidimensional functional diversity indices for a multifaceted framework in functional ecology. Ecology 2008, 89, 2290–2301. [Google Scholar] [CrossRef] [Green Version]
  55. Brown, L.E.; Milner, A.M. Rapid loss of glacial ice reveals stream community assembly processes. Glob. Chang. Biol. 2012, 18, 2195–2204. [Google Scholar] [CrossRef] [Green Version]
  56. Pavoine, S.; Dolédec, S. The apportionment of quadratic entropy: A useful alternative for partitioning diversity in ecological data. Environ. Ecol. Stat. 2005, 12, 125–138. [Google Scholar] [CrossRef]
  57. Mason, N.; Mouillot, D.; Lee, W.G.; Setl, W.H. Functional Richness, Functional Evenness and Functional Divergence: The Primary Components of Functional Diversity. Oikos 2005, 111, 112–118. [Google Scholar] [CrossRef]
  58. Dolédec, S.; Chessel, D.T.; ter Braak, C.; ter Braak, C.J.F.; Champely, S. Matching species traits to environmental variables: A new three-table ordination method. Environ. Ecol. Stat. 1996, 3, 143–166. [Google Scholar] [CrossRef]
  59. Pallottini, M.; Cappelletti, D.; Fabrizi, A.; Gaino, E.; Goretti, E.; Selvaggi, R.; Céréghino, R. Macroinvertebrate functional trait responses to chemical pollution in agricultural–industrial landscapes. River Res. Appl. 2017, 33, 505–513. [Google Scholar] [CrossRef]
  60. Laliberté, E.; Legendre, P.; Shipley, B. FD: Measuring functional diversity from multiple traits, and other tools for functional ecology. R Package Version 2014, 1, 1–12. [Google Scholar]
  61. Dray, S.; Dufour, A.-B. The ade4 Package: Implementing the Duality Diagram for Ecologists. J. Stat. Softw. 2007, 22, 1–20. [Google Scholar] [CrossRef] [Green Version]
  62. Dray, S.; Choler, P.; Dolédec, S.; Peres-Neto, P.R.; Thuiller, W.; Pavoine, S.; ter Braak, C.J.F. Combining the fourth-corner and the RLQ methods for assessing trait responses to environmental variation. Ecology 2014, 95, 14–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Jonsson, M.; Burrows, R.M.; Lidman, J.; Fältström, E.; Laudon, H.; Sponseller, R.A. Land use influences macroinvertebrate community composition in boreal headwaters through altered stream conditions. Ambio 2017, 46, 311–323. [Google Scholar] [CrossRef] [Green Version]
  64. Luo, K.; Hu, X.; He, Q.; Wu, Z.; Cheng, H.; Hu, Z.; Mazumder, A. Impacts of rapid urbanization on the water quality and macroinvertebrate communities of streams: A case study in Liangjiang New Area, China. Sci. Total Environ. 2018, 621, 1601–1614. [Google Scholar] [CrossRef]
  65. Martins, R.; Couceiro, S.; Melo, A.; Moreira, M.; Hamada, N. Effects of urbanization on stream benthic invertebrate communities in Central Amazon. Ecol. Indic. 2017, 73, 480–491. [Google Scholar] [CrossRef]
  66. Schriever, T.A.; Bogan, M.T.; Boersma, K.S.; Cañedo-Argüelles, M.; Jaeger, K.L.; Olden, J.D.; Lytle, D.A. Hydrology shapes taxonomic and functional structure of desert stream invertebrate communities. Freshw. Sci. 2015, 34, 399–409. [Google Scholar] [CrossRef] [Green Version]
  67. Gerisch, M.; Agostinelli, V.; Henle, K.; Dziock, F. More species, but all do the same: Contrasting effects of flood disturbance on ground beetle functional and species diversity. Oikos 2012, 121, 508–515. [Google Scholar] [CrossRef]
  68. Péru, N.; Dolédec, S. From compositional to functional biodiversity metrics in bioassessment: A case study using stream macroinvertebrate communities. Ecol. Indic. 2010, 10, 1025–1036. [Google Scholar] [CrossRef]
  69. Fierro, P.; Valdovinos, C.; Lara, C.; Saldias, G.S. Influence of Intensive Agriculture on Benthic Macroinvertebrate Assemblages and Water Quality in the Aconcagua River Basin (Central Chile). Water 2021, 13, 492. [Google Scholar] [CrossRef]
  70. Odume, O.N. Searching for urban pollution signature and sensitive macroinvertebrate traits and ecological preferences in a river in the Eastern Cape of South Africa. Ecol. Indic. 2020, 108, 105759. [Google Scholar] [CrossRef]
  71. Odume, O.N.; Palmer, C.G.; Arimoro, F.O.; Mensah, P.K. Patterns of chironomid body-size distribution in an effluent-impacted river in the Eastern Cape, South Africa. Afr. J. Aquat. Sci. 2014, 39, 377–386. [Google Scholar] [CrossRef]
  72. Statzner, B.; BÊChe, L.A. Can biological invertebrate traits resolve effects of multiple stressors on running water ecosystems? Freshw. Biol. 2010, 55, 80–119. [Google Scholar] [CrossRef]
  73. Feio, M.J.; Dolédec, S. Integration of invertebrate traits into predictive models for indirect assessment of stream functional integrity: A case study in Portugal. Ecol. Indic. 2012, 15, 236–247. [Google Scholar] [CrossRef]
  74. Tomanova, S.; Usseglio-Polatera, P. Patterns of benthic community traits in neotropical streams: Relationship to mesoscale spatial variability. Fundam. Appl. Limnol. 2007, 170, 243–255. [Google Scholar] [CrossRef]
  75. Bryce, S.A.; Lomnicky, G.A.; Kaufmann, P.R. Protecting sediment-sensitive aquatic species in mountain streams through the application of biologically based streambed sediment criteria. J. N. Am. Benthol. Soc. 2010, 29, 657–672. [Google Scholar] [CrossRef] [Green Version]
  76. Brown, B.L.; Swan, C.M. Dendritic network structure constrains metacommunity properties in riverine ecosystems. J. Anim. Ecol. 2010, 79, 571–580. [Google Scholar] [CrossRef]
  77. Heino, J.; Melo, A.S.; Siqueira, T.; Soininen, J.; Valanko, S.; Bini, L.M. Metacommunity organisation, spatial extent and dispersal in aquatic systems: Patterns, processes and prospects. Freshw. Biol. 2015, 60, 845–869. [Google Scholar] [CrossRef]
  78. Miserendino, M.L.; Masi, C.I. The effects of land use on environmental features and functional organization of macroinvertebrate communities in Patagonian low order streams. Ecol. Indic. 2010, 10, 311–319. [Google Scholar] [CrossRef]
  79. Statzner, B. Mineral grains in caddisfly pupal cases and streambed sediments: Assessing resource use and its limitation across various river types. Ann. Limnol.-Int. J. Limnol. 2011, 47, 103–118. [Google Scholar] [CrossRef]
  80. Burdon, F.J.; Harding, M. Habitat loss drives threshold response of benthic invertebrate communities to deposited sediment in agricultural streams. Ecol. Appl. 2013, 23, 1036–1047. [Google Scholar] [CrossRef] [PubMed]
  81. Edegbene, A.O.; Adam, M.B.; Gambo, J.; Osimen, E.C.; Ikomi, R.B.; Ogidiaka, E.; Omovoh, G.O.; Akamagwuna, F.C. Searching for indicator macroinvertebrate traits in an Afrotropical riverine system: Implication for ecosystem biomonitoring and sustainability. Environ. Monit. Assess. 2021, 193, 711. [Google Scholar] [CrossRef]
  82. Chapman, L.J.; Schneider, K.R.; Apodaca, C.; Chapman, C.A. Respiratory ecology of macroinvertebrates in a swamp-river system of east Africa. Biotropica 2004, 36, 572–585. [Google Scholar] [CrossRef]
  83. Chadwick, M.A. Stream Ecology: Structure and Function of Running Waters. Freshw. Biol. 2008, 53, 1914. [Google Scholar] [CrossRef]
Figure 1. The distribution of 44 sampling sites in Jinan, Shandong Province, China.
Figure 1. The distribution of 44 sampling sites in Jinan, Shandong Province, China.
Water 14 03793 g001
Figure 2. Box plots of taxonomic diversity indexes of richness, Simpson’s diversity, Shannon–Wiener index, and evenness of macroinvertebrates distributed across different land use types. The lines in the boxes represent the median, the boxes represent the 1st and 3rd quartiles, and the lines at the top and bottom of the boxes represent the 5th and 95th percentiles, respectively. (a) Species richness in different land use types; (b) species diversity in different land use types; (c) species Shannon–Wiener index in different land use types; (d) species evenness in different land use types.
Figure 2. Box plots of taxonomic diversity indexes of richness, Simpson’s diversity, Shannon–Wiener index, and evenness of macroinvertebrates distributed across different land use types. The lines in the boxes represent the median, the boxes represent the 1st and 3rd quartiles, and the lines at the top and bottom of the boxes represent the 5th and 95th percentiles, respectively. (a) Species richness in different land use types; (b) species diversity in different land use types; (c) species Shannon–Wiener index in different land use types; (d) species evenness in different land use types.
Water 14 03793 g002
Figure 3. Box plots of functional diversity indexes of richness functional richness, functional evenness, functional divergence, functional dispersion, and Rao’s quadratic entropy of macroinvertebrates distributed across different land use types. * significance level at 0.05 for non-parametric tests. ** significance level at 0.01 for non-parametric tests. (a) Functional richness (FRic) in different land use types; (b) functional evenness (FEve) in different land use types; (c) functional divergence (FDiv) in different land use types; (d) functional dispersion (FDis) in different land use types; (e) Rao’s quadratic entropy (RaoQ) in different land use types.
Figure 3. Box plots of functional diversity indexes of richness functional richness, functional evenness, functional divergence, functional dispersion, and Rao’s quadratic entropy of macroinvertebrates distributed across different land use types. * significance level at 0.05 for non-parametric tests. ** significance level at 0.01 for non-parametric tests. (a) Functional richness (FRic) in different land use types; (b) functional evenness (FEve) in different land use types; (c) functional divergence (FDiv) in different land use types; (d) functional dispersion (FDis) in different land use types; (e) Rao’s quadratic entropy (RaoQ) in different land use types.
Water 14 03793 g003aWater 14 03793 g003b
Figure 4. The structure of 36 categories of eleven traits at three land use types. Each bar represents the mean percentage with standard deviation of each trait category. (a) The trait of voltinism in different land use types; (b) the trait of development in different land use types; (c) the trait of dispersal in different land use types; (d) the trait of drift in different land use types; (e) the trait of attachment in different land use types; (f) the trait of armoring in different land use types; (g) the trait of shape in different land use types; (h) the trait of respiration in different land use types; (i) the trait of size in different land use types; (j) the trait of rheophily in different land use types; (k) the trait of habit in different land use types; (l) the trait of trophic habit in different land use types. Various letters denote the significant differences according to Kruskal-Wallis test at a 0.05 probability level.
Figure 4. The structure of 36 categories of eleven traits at three land use types. Each bar represents the mean percentage with standard deviation of each trait category. (a) The trait of voltinism in different land use types; (b) the trait of development in different land use types; (c) the trait of dispersal in different land use types; (d) the trait of drift in different land use types; (e) the trait of attachment in different land use types; (f) the trait of armoring in different land use types; (g) the trait of shape in different land use types; (h) the trait of respiration in different land use types; (i) the trait of size in different land use types; (j) the trait of rheophily in different land use types; (k) the trait of habit in different land use types; (l) the trait of trophic habit in different land use types. Various letters denote the significant differences according to Kruskal-Wallis test at a 0.05 probability level.
Water 14 03793 g004
Figure 5. The summary of combining RLQ and fourth-corner analysis for macroinvertebrates traits and environmental variables. Red/blue cells indicate significant positive/negative relationships, respectively; grey cells indicate no significant relationships. The codes for macroinvertebrate traits and environmental variables are explained by Table 1 and Table 2, respectively.
Figure 5. The summary of combining RLQ and fourth-corner analysis for macroinvertebrates traits and environmental variables. Red/blue cells indicate significant positive/negative relationships, respectively; grey cells indicate no significant relationships. The codes for macroinvertebrate traits and environmental variables are explained by Table 1 and Table 2, respectively.
Water 14 03793 g005
Table 1. Macroinvertebrate traits and categories of used in the assessment of the trait–environment relationships across three different land use types in the present study.
Table 1. Macroinvertebrate traits and categories of used in the assessment of the trait–environment relationships across three different land use types in the present study.
TraitsCategoryCodeReference
VoltinismSemivoltineVolt1[3,32,44]
UnivoltineVolt2
Bi- or multivoltine Volt3
DevelopmentFast seasonalDeve1[32]
Slow seasonalDeve2
NonseasonalDeve3
DispersalNoDisp1[3,32]
LowDisp2
HighDisp3
DriftRareDrft1[3,32]
CommonDrft2
Very abundantDrft3
AttachmentFree-rangingAtch1[32]
Sessile, sedentaryAtch2
BothAtch3
ArmoringNone (soft-bodied forms)Armr1[32]
Poor (heavily sclerotized)Armr2
Good (e.g., some cased caddisflies)Armr3
ShapeStreamlined (flat, fusiform)Shap1[32,44]
Not streamlined (cylindrical, round, or bluff)Shap2
RespirationTegumentResp1[3,32,44]
GillsResp2
Air (Plastron or spiracle)Resp3
SizeSmall (<9 mm)Size1[3,32,44]
Medium (9–16 mm)Size2
Large (>16 mm)Size3
RheophilyDepositional onlyRheo1[3,32]
Depositional and erosionalRheo2
ErosionalRheo3
HabitBurrowHabit1[3,32,44]
ClimbHabit2
SprawlHabit3
ClingHabit4
SwimHabit5
Trophic habitCollector–gathererTrop1[3,32,44]
Collector–filtererTrop2
Herbivore (scraper, piercer, and shedder)Trop3
Predator (piercer and engulfer)Trop4
Shredder (detritivore)Trop5
Table 2. Comparison of environmental variables (mean ± SD) in different land use type sites. Significant differences in environmental variables across different types of land use are indicated in bold (Kruskal-Wallis Test: p < 0.05). Mountain–hilly areas (15 sites); urban areas (15 sites); agricultural areas (14 sites).
Table 2. Comparison of environmental variables (mean ± SD) in different land use type sites. Significant differences in environmental variables across different types of land use are indicated in bold (Kruskal-Wallis Test: p < 0.05). Mountain–hilly areas (15 sites); urban areas (15 sites); agricultural areas (14 sites).
VariableMountain-Hilly AreaUrban AreaAgricultural Areapχ2
Water temperature (°C)21.57 ± 4.8921.48 ± 4.4220.94 ± 4.880.9300.145
pH8.24 ± 0.377.91 ± 0.358.17 ± 0.33<0.00133.132
Electrical conductivity (EC) (s/cm)634.65 ± 393.511145.73 ± 700.551723.13 ± 863.66<0.001118.534
Dissolved oxygen concentration (DO) (mg/L)8.15 ± 1.527.44 ± 2.357.49 ± 2.450.1993.222
Total nitrogen (TN) (mg/L)3.94 ± 4.165.86 ± 4.494.69 ± 9.32<0.00115.651
Ammonium nitrogen (NH4-N) (mg/L)0.96 ± 3.111.19 ± 1.642.05 ± 8.78<0.00111.931
Nitrate nitrogen (NO3-N) (mg/L)2.2 ± 1.953.81 ± 3.731.35 ± 1.72<0.00127.563
Chemical oxygen demand (COD) (mg/L)20.56 ± 17.8921.88 ± 14.2832.18 ± 31.45<0.00127.522
Biochemical oxygen demand (BOD) (mg/L)2.97 ± 4.253.58 ± 3.765.48 ± 7.53<0.00114.85
Total phosphorus (TP) (mg/L)0.28 ± 0.690.33 ± 0.580.39 ± 1.02<0.00124.121
Farmland (%)34.31 ± 33.6724.52 ± 24.3359.90 ± 28.38<0.0019.293
Forest land (%)4.86 ± 16.36 6.074
Grassland (%)2.91 ± 8.82 6.074
Water area (%)24.50 ± 24.5717.39 ± 28.7612.35 ± 23.550.1903.319
Urban land (%)30.40 ± 28.0558.03 ± 41.0926.52 ± 26.750.1204.233
Table 3. Results of the RLQ analysis. Eigenvalues and the first two RLQ axes account for the percentage of total inertia; covariation is the covariance between the two new sets of factor scores projected onto the first two RLQ axes; correlation is the correlation between the two new sets of factor scores projected onto the first two RLQ axes. Cumulative inertia is the variance of each set of factor scores; ratio is the variance of the first RLQ axis as a percentage of the first axis analyzed separately.
Table 3. Results of the RLQ analysis. Eigenvalues and the first two RLQ axes account for the percentage of total inertia; covariation is the covariance between the two new sets of factor scores projected onto the first two RLQ axes; correlation is the correlation between the two new sets of factor scores projected onto the first two RLQ axes. Cumulative inertia is the variance of each set of factor scores; ratio is the variance of the first RLQ axis as a percentage of the first axis analyzed separately.
Axes1Axes2
Eigenvalue18.123.10
% of total co-inertia78.47%13.43%
Covariation4.261.76
Correlation0.660.46
Cumulative inertia (environment)3.706.53
Ratio (environment)84.61%91.26%
Cumulative inertia (traits)11.3516.47
Ratio (traits)99.08%83.98%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lu, A.; Li, J.; Zheng, B.; Yin, X. Effect of Different Land Use Types on the Taxonomic and Functional Diversity of Macroinvertebrates in an Urban Area of Northern China. Water 2022, 14, 3793. https://doi.org/10.3390/w14233793

AMA Style

Lu A, Li J, Zheng B, Yin X. Effect of Different Land Use Types on the Taxonomic and Functional Diversity of Macroinvertebrates in an Urban Area of Northern China. Water. 2022; 14(23):3793. https://doi.org/10.3390/w14233793

Chicago/Turabian Style

Lu, Aoran, Jiaxin Li, Biao Zheng, and Xuwang Yin. 2022. "Effect of Different Land Use Types on the Taxonomic and Functional Diversity of Macroinvertebrates in an Urban Area of Northern China" Water 14, no. 23: 3793. https://doi.org/10.3390/w14233793

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