Estimating the nutrient thresholds of a typical tributary in the Liao River basin, Northeast China
Estimating regional nutrient criteria for streams and rivers is a key step toward protecting river water quality and restoring the health of aquatic ecosystems. Using a multivariable statistical analysis technique, nutrients were identified as the main factor influencing the degradation of the benthic macroinvertebrate community. Three chemical methods (the reference stream distribution approach, all-streams distribution approach and Y-intercept approach) and one biological method (the stress-response approach) were applied to evaluate the nutrient thresholds in the Qing River basin. The reference stream distribution approach and all-streams distribution approach were based on calculating a predetermined percentile of reference streams and all-streams water quality data set, respectively. The Y-intercept approach was based on determining the influence of human activity on water quality by linear regression models. The biological method was based on the response of the benthic macroinvertebrate community structure to changes in water quality. The chemical thresholds were 0.750–1.288 mg/L for total nitrogen (TN) and 0.035–0.046 mg/L for total phosphorus (TP); the biological thresholds were 1.050–1.655 for TN and 0.052–0.101 for TP. The results from the chemical approaches were verified using the biological method, resulting in preliminarily recommended thresholds of 1.000 mg/L TN and 0.040 mg/L TP in the Qing River system.
Eutrophication is a common environmental problem throughout the world, primarily occurring in enclosed or semi-enclosed water bodies such as lakes, reservoirs and bays1,2,3,4. With the excessive exploitation of natural resources and the rapid development of industry and agriculture, the nutrient concentrations in rivers have increased in recent years4,5,6,7. Excessive nutrient levels can negatively impact river water and cause various problems, including affecting the health of humans and livestock, changing the structures of river ecosystems, and reducing the aesthetic and recreational features of the rivers2,8,9,10. Furthermore, because rivers are the source of water for many lakes, coastal waters and wetlands, increased nutrients may increase eutrophication in these downstream waters as well4,8,11,12.
Establishment of water quality criteria and standards is a key step toward protecting water quality8. A number of developed countries have been or are formulating nutrient criteria for surface waters. The United States Environmental Protection Agency (US EPA) issued a national strategy for developing a regional set of nutrient criteria in 19989, and a nutrient criteria technical guidance manual for rivers and streams was published two years later8. Dodds et al.3 estimated the boundaries of oligotrophic, mesotrophic and eutrophic states of rivers using published data on the water quality of large- and medium-scale rivers in the United States with a frequency distribution method in 1998. In 2000, the US EPA8 published methods for evaluating reference conditions, including a reference stream distribution approach based on a predetermined percentile of the water quality data set of reference streams and an all-streams distribution approach based on a predetermined percentile of the water quality data set of the general stream population. Smith et al.13 evaluated the natural background concentrations of nutrients in streams and rivers of the contiguous United States by integrating nutrient yield models with the SPARROW transport model from existing reference streams. Dodds and Oakes14 employed a Y-intercept approach for establishing reference conditions across watersheds based on estimating the influence of land use on river nutrient concentrations by means of multiple linear regression analysis. Some researchers15,16,17 have evaluated the environmental thresholds for nutrients in streams using a stress-response approach based on the relationships between biological response parameters and river water quality parameters.
The research on water quality criteria and standards is sparse in China compared with developed countries. Water quality criteria are mainly derived from the relevant works in foreign countries18, which reduces the applicability of these criteria due to the regional variability of aquatic biota19. Nutrient levels and reactions vary greatly across a country, and the US EPA has suggested that recommended nutrient criteria must reflect the characteristics of different geographical regions and types of water bodies9.
The existing Chinese surface water environmental quality standards (GB3838-2002) have wide coverage and poor regional specificity, which may result in problems of under- or over-protection for river water management. In addition, controls on total nutrient quantities consider mainly point source parameters in China; few consider non-point source parameters. The water quality standards lack total nitrogen (TN) control standards for streams and rivers, and as such, they cannot address the main causes of nutrient pollution in water bodies via non-point source pollution. Establishing regional nutrient criteria and standards is a fundamental step toward reducing nutrient pollution in aquatic ecosystems in China.
The study of nutrient criteria for streams and rivers is still in its infancy in China. Wu et al.20 (2010) have proposed Chinese water environment quality criteria. Chen and Lu21 have established nutrient criteria for streams and rivers with three chemical methods in the Cao-E River basin. Wu et al.22 have estimated river nutrient thresholds based on benthic macroinvertebrate assemblages in the upper reaches of Xitiao Stream in Zhejiang, China. However, there are few comparative studies of nutrient criteria for streams and rivers based on chemical and biological approaches in China.
Exploring the influence of river basin characteristics on river water quality can help identify the main pollutants that affect benthic macroinvertebrate assemblages and evaluate the nutrient thresholds. According to the geographical features of Northeast China, nutrient criteria of streams and rivers were evaluated in the Qing River system using a combination of chemically derived threshold approaches (the reference stream distribution approach, all-streams distribution approach, Y-intercept approach) and a biologically derived threshold approach (stress-response approach). Each of these approaches has its own application scope and shortcoming16,21. The chemically derived thresholds were cross-calibrated with the biologically derived thresholds. The use of these three chemical methods and one biological method exploits the advantages and avoids the disadvantages of each method. These four complementary methods can provide strong support for the establishment of water quality criteria and standards in China and scientific and technological support for the implementation of water pollution controls.
Results and Discussion
Influence of river basin characteristics on river water quality
The river basin characteristics, the river water quality parameters and the potential relationships between them at the sampling sites are discussed to identify the main pollutants and further assist in evaluating the nutrient thresholds in the Qing River basin. Table 1 gives the descriptive statistics for the physical and chemical parameters of the water, the river habitat assessment indicators and land use. Sulfate and pH exhibited less spatial variability than the other physical and chemical parameters, with coefficients of variation of 22.27% and 8.06%, respectively. The fecal coliform count (FCC), total bacterial count (TBC) and total phosphorus (TP) exhibited more spatial variability, with coefficients of variation of 134.08%, 100.21% and 100.00%, respectively. Sites with higher TP concentrations were located in the middle and lower reaches of the Qing River (monitoring sites 3, 13, 16 and 34; Fig. 1).
Table 1 Characteristics of the Qing River.
Full size table
Locations of sampling stations. Latitude and longitude were measured using a hand-held GPS unit, and ArcGIS 10.0 Desktop GIS software (http://www.esri.com/arcgis/about-arcgis) was used to plot the sampling stations.
Full size image
The average values of the NH4+-N concentration, CODCr, BOD5 and FCC (Table 1) were 0.5 mg/L, 16.2 mg/L, 3.2 mg/L and 723.4 num/L, representing grade II, grade III, grade III and grade II, respectively, based on the environmental quality standards for surface water in China (GB3838-2002). The sub-watersheds of the 35 monitoring sites in the Qing River system had high land use variability. The ion concentration, represented by the EC, was very significantly associated with the chloride concentration, sulfate concentration and land use area percentage and significantly associated with the TP concentration (Table 2). Though the EC variations were mainly determined by geographical differences, increased EC values are usually an indicator of urbanization23,24, which was confirmed by the strong relationship between the chloride concentration and EC in the Qing River basin. The sampling sites with higher nutrient concentrations usually had higher EC values because poor runoff management led to higher nutrient and ion concentrations. It has been shown previously that EC usually increases with the urban land area percentage, even in limestone regions with high background EC levels17.
Table 2 Pearson correlation coefficients among the water chemistry parameters, habitat parameters and land uses at the sampling stations in the Qing River basin.
Full size table
The average river habitat score of the sub-watersheds was 131 (total score of 200); the lowest average scores (12) were for channel alteration, vegetative protection and riparian vegetative zone width (Table 1), which showed that these parameters were the most limiting river habitat parameters in the Qing River basin. The highest average scores (14) were for embeddedness and sediment deposition, showing that these parameters were at acceptable levels. The river habitat parameters were very significantly negatively associated with CODCr and significantly negatively associated with the nutrient concentrations (TN or/and TP), and the TN concentration was very significantly associated with the NO3−N concentration (Table 2), which indicated that the main pollution type was non-point source. The TN concentration was higher than the TP concentration by multiple orders of magnitude, and the molar N:P ratios at most of the sampling sites were greater than 17 (Table 1).
The TN concentration, EC, CODCr and BOD5 were very significantly associated with the main land use percentage in the Qing River system (Table 2), which indicated that land use made the greatest contribution to the nutrient, ion, and oxygen consumption in the River. The total habitat scores decreased as the urban land and cropland percentages increased (Table 2). The nutrient concentrations in the river water increased with increasing human activities. The correlation analysis (Table 2) showed that the TP level was closely associated with the urban land, cropland and woodland percentages (r values of 0.779, 0.497 and −0.549, respectively). Higher ion and TP concentrations (monitoring sites 3, 13, 16 and 34; Fig. 1) were found near urban land (wastewater treatment plant outflow, domestic wastewater and livestock wastewater). The TN concentration was very significantly associated with the cropland and woodland percentages (r values of 0.650 and −0.615, respectively) because the research area is mainly affected by agricultural non-point source pollution21,24,25,26,27. Poor agricultural management practices and mining activities can input suspended solids and sediments, including high levels of nutrients into rivers21,25,28. pH had a low correlation with the urban land percentage (r value of 0.121) mainly because the area was located in the agricultural region of the Qing River basin (the research area contained little urban land, and most of the urban land was located downstream of the Qing River basin). The relationship between land use and nutrient concentrations was explored by multiple regression analysis, and the results (Table 3) showed that cropland and population density could explain 42.3% and 66.1% of the variance of TN and TP, respectively. The cropland and residential domestic pollution contributed substantially to the increased river nutrient concentrations.
Table 3 Results of the regression analyses between nutrient (TN and TP) concentrations and land uses in the Qing River basin.
Full size table
Relationship between benthic macroinvertebrate community characteristics and river nutrient concentrations
Benthic macroinvertebrates, the most commonly used environmental indicator organisms in biological water quality monitoring, can represent the cumulative effect of water body disturbance and pollution stress in a past period of time; many studies16,17,22 indicated that nutrient pollution can negatively affects the benthic macroinvertebrate community structure. A total of 70 benthic macroinvertebrate taxa belonging to 11 orders (Diptera, Ephemeroptera, Trichoptera, Plecoptera, Coleoptera, Hemiptera, Decapoda, Basommatophora, Gnathobdellida, Rhynchobdellida, Oligochaeta: Plesiopora) and three phyla (Annelida, Mollusca, Arthropoda) were collected over the study period. The relative abundance of taxonomic groups at the order level revealed that the most abundant taxa in the study area, especially downstream, were in the order Diptera (37.52%). Ephemeroptera (26.16%) and Trichoptera (19.23%) were the next most abundant orders in the study area and were the most represented groups in the forest areas. Four biological metrics, which are based on the evaluation of macroinvertebrate measures and are commonly used for assessing water quality and river ecosystem health, were significantly related to the TN and TP concentrations (Table 4) and were calculated for all the sampling stations (Fig. 1): the Ephemeroptera, Plecoptera, and Trichoptera (EPT) relative abundance, EPT taxonomic richness, modified family biotic index (FBI), and Diptera and non-insect relative abundance. The EPT relative abundance percentages (range of 8–74%) and EPT taxonomic richness (range of 1–18) were higher at monitoring sites 27–30 (Fig. 1). The modified FBI values (range of 3.79–6.79) were the highest at monitoring sites 2–5. The Diptera and non-insect relative abundance percentages ranged from 16% to 86%, were the lowest at monitoring sites 27–30 and 32 (16–19%) and the highest at monitoring sites 2–5 (76–86%).
Table 4 Pearson correlation coefficients between the environmental variables and benthic macroinvertebrate metrics.
Full size table
Unbalanced levels of nutrients, ions and DO are the possible causes of biological degradation in rivers29,30. The results (Table 4) showed that DO was not significantly related with the 4 benthic macroinvertebrate indices. Although the indices were significantly associated with EC, these correlations were weaker than their correlations with the nutrient concentrations (TN and TP). TP and TN were significantly associated with the EPT relative abundance % (r = −0.430, p < 0.01, and r = −0.409, p < 0.05, respectively), EPT taxonomic richness (r = −0.433, p < 0.01, and r = −0.398, p < 0.05, respectively), modified FBI (r = 0.383, p < 0.05, and r = 0.404, p < 0.05, respectively) and Diptera and non-insect relative abundance (r = 0.406, p < 0.05, and r = 0.424, p < 0.05, respectively), and TP also had very significant relationships with the EPT relative abundance % and EPT taxonomic richness (Table 4). The EPT relative abundance % (r = −0.459, p < 0.05, and r = 0.771, p < 0.01, respectively), EPT taxonomic richness (r = −0.704, p < 0.01, and r = 0.710, p < 0.01, respectively), modified FBI (r = 0.741, p < 0.01, and r = −0.743, p < 0.01, respectively) and Diptera and non-insect relative abundance (r = 0.799, p < 0.01, and r = −0.801, p < 0.01, respectively) had significant relationships with the percentages of cropland and woodland (Table 4), which showed that the substantial nutrient loss caused by agricultural cultivation has altered the community structure of benthic macroinvertebrates. As shown in Fig. 2, the biotic indices varied considerably at low nutrient concentrations but varied only slightly at high concentrations. The EPT relative abundance % and EPT taxonomic richness increased first and then decreased as the TP concentration increased. The EPT relative abundance % and EPT taxonomic richness decreased as the TN concentration increased. The modified FBI and Diptera and non-insect relative abundance increased as the TP and TN concentrations increased. Nutrients were not the main factors affecting the aquatic species distributions at low concentrations, and at high concentrations, they obviously hindered the growth of certain species; these results are similar to those of previous studies16,22.
Relationships between the benthic macroinvertebrate metrics and nutrient (TN and TP) concentrations in the Qing River basin. The smooth lines show the locally weighted scatterplot smoothing (LOWESS).
Full size image
Nutrient criteria evaluation
Based on the Y-intercept approach14, the TP criterion was 0.035 mg/L with a 95% confidence interval of 0.019–0.050 mg/L in the Qing River system, and the TN criterion was 0.749 mg/L with a 95% confidence interval of 0.477–1.021 mg/L (Table 5). The TN and TP thresholds evaluated according to the reference stream distribution approach were 1.288 mg/L and 0.046 mg/L, respectively; the TN and TP thresholds evaluated using the all-streams distribution approach were 0.724 mg/L and 0.024 mg/L, respectively (Table 6). Nutrient criteria can be evaluated with a stress-response approach based on benthic macroinvertebrate community characteristics. Because the EPT relative abundance %, EPT taxonomic richness, modified FBI and Diptera and non-insect relative abundance play important roles in river health assessments16,31,32, these four biotic indices that were significantly correlated with both TN and TP were selected to establish the river nutrient criteria. The number of pollution-sensitive species decreased rapidly, and the stability of the community structure deteriorated when the nutrient (TN and TP) concentrations exceeded their respective thresholds (Fig. 2). Regression tree analysis with the four biotic indices (EPT relative abundance %, EPT taxonomic richness, modified FBI and Diptera and non-insect relative abundance) showed that the evaluated TN thresholds were 1.655 mg/L, 1.050 mg/L, 1.574 mg/L and 1.610 mg/L, respectively, and the evaluated TP thresholds were 0.084 mg/L, 0.052 mg/L, 0.101 mg/L and 0.084 mg/L, respectively (Fig. 2, Table 6). Among the evaluated nutrient thresholds obtained using the four approaches, the values determined with the stress-response approach were highest, followed by those determined with the reference stream distribution approach, the Y-intercept approach and the all-streams distribution approach, in decreasing order (Table 6).
Table 5 Predicted nutrient (TN and TP) thresholds and their confidence intervals determined using the Y-intercept approach for the Qing River basin.
Full size table
Table 6 Evaluated nutrient (TN and TP) thresholds for the Qing River basin in Northeast China.
Full size table
There are few reference streams or reaches in the developed region, and the nutrient loads of reference streams or reaches can vary significantly in the same ecological region due to differences in runoff13,14,33. The 5–25% frequency distribution for all streams was selected to evaluate the nutrient thresholds9. This 5–25% frequency distribution was a deduced value, and the selected frequency distribution can generate nutrient thresholds that result in under- or over-protection of the water quality34. The Y-intercept approach is simple and widely applicable and is primarily used in studies considering several major influencing factors, including land use and population density. However, river water quality is affected by many factors, such as livestock and poultry breeding35,36,37, point source pollution24, agricultural management practices38,39, atmospheric deposition39,40, topography41,42, and the riparian buffer26,43,44. In addition, the agricultural watershed is a region of internal change affected by human activities; even all human activities, including agricultural activities, were ceased, it is unclear whether nutrient outputs would be restored to their natural levels16. Some studies16,45 have indicated that the other factors might be equally or more important than nutrient levels in affecting the aquatic community structure of a certain area. In contrast, the stress-response approach based on biological attributes allows identification of nutrient thresholds consistent with good ecological condition; however, such models often have low explanatory power, especially for streams14,16,45.
Each nutrient criteria evaluation approach has disadvantages, and the comparison of the chemically derived thresholds (reference stream distribution approach, all-streams distribution approach and Y-intercept approach) and the biologically derived threshold (the stress-response approach based on benthic macroinvertebrate community structure) strengthens the validity of the nutrient criteria for the protection and rehabilitation of aquatic life. The nutrient thresholds estimated using the chemically derived approaches were 0.724–1.288 mg/L for TN and 0.024–0.046 mg/L for TP (Table 6). The chemically derived TN and TP thresholds could be used as the nutrient thresholds for the Qing River basin because they mainly lie within the bounds of all biological thresholds16,17,32. In addition, because the pollution in the study area can be categorized as have an agricultural non-point source at most sampling sites, the recommended preliminary nutrient thresholds are 1.000 mg/L for TN and 0.040 mg/L for TP (approximately equal to the averages of the chemical thresholds and less than the biological thresholds). Further work is required to evaluate other biological responses (e.g., algal and fish biological metrics) to nutrients and perform field evaluations of the applicability of the preliminary recommended nutrient thresholds.
Materials and Methods
The Qing River, a main tributary of the Liao River, is located in Northeast China. The Qing River basin is between N42.30°–42.70° and E123.80°–124.8° and has a catchment area of 4785 km2 (Fig. 1). The Qing River has a length of 217 km and 8 tributaries, including the Kou and Mazhong Rivers. The Qing River basin, with an average annual precipitation and temperature of 692 mm and 6.5 °C, respectively, has a temperate monsoon climate, with precipitation mainly occurring in the summer (July and August). The terrain is quite varied, from hills and mountains to alluvial plains. The vegetation coverage decreases from east to west. Land uses mainly include woodland, cropland, water and urban land. The eastern region of the Qing River basin, with little human activity, is predominantly woodland, and the rivers in this area have higher gradients, faster flows and a larger average grain size of riverbed sediments than those in the western region. Residential density, cropland, urban land percentage and human activities are all higher in the western region of the Qing River basin46. The riparian zone is more consolidated, and the river pollution load is higher in the urban areas. The Qing River basin includes Kaiyuan City, the Qinghe District and Xifeng County, with a population of over one million and an annual GDP of more than 20 billion Chinese yuan; the agricultural activities are dominated by cropland cultivation and animal husbandry in the Qinghe District47.
Sample collection and river habitat investigation
The riparian characteristics, physical habitat parameters, hydrochemistry and biological integrity indices were evaluated in different periods in 2009 and 2011. A habitat survey was conducted utilizing the evaluation indices and methods for river habitats48,49; 10 indices were rated on a scale of 0–20 at the sampling sites, with high scores indicating a high-quality habitat. The monitored parameters included the TN, TP, ammonia nitrogen (NH4+-N), dissolved oxygen (DO), chloride, and sulfate concentrations; 5-day biochemical oxygen demand (BOD5); potassium bichromate index (CODCr); pH; temperature (T); electrical conductivity (EC); FCC; and TBC. The sampling, preservation and analytical procedures were performed according to the national standard methods of China50; the analytical methods are listed in Table 7.
Table 7 Analytical methods used to assess water quality parameters of the Qing River basin.
Full size table
Macroinvertebrate samples were collected at four locations within each of the 35 sites using a surber sampler with a size of 0.3 × 0.4 m and a pore diameter of 425 µm. The net was placed on rocky river substrates, and an area of 0.12 m2 upstream of the net was disturbed to dislodge and wash the macroinvertebrates into the net. Individual rocks in the sampled area were picked up, and the attached organisms were removed. The four replicate samples were combined, and all of the obtained aquatic organisms were preserved in 70% ethanol51,52.
Derived nutrient thresholdsReference stream distribution approach
The 75% frequency distribution of the nutrient data obtained for reference streams was applied to estimate the thresholds8.
All-streams distribution approach
The 5th to 25th percentile of the frequency distribution of the nutrient data obtained for all streams was selected to estimate the nutrient thresholds in the study area8.
Multiple linear regression analysis was used to evaluate the thresholds. The logarithms of the in-stream nutrient concentrations are the dependent parameters, the land use percentages are the independent parameters, and the intercept of the regression model represents the nutrient thresholds14.
Based on the macroinvertebrate data, nine macroinvertebrate community indices were calculated: species richness, macroinvertebrate density, Shannon’s diversity index, the modified Shannon-Wiener index, evenness, EPT taxonomic richness, EPT relative abundance, the modified FBI, and the Diptera and non-insect relative abundance. Four community indices (EPT relative abundance %, EPT taxonomic richness, modified FBI and Diptera and non-insect relative abundance) were significantly correlated (p < 0.05) with both TN and TP and were therefore retained for further evaluation. EPT taxonomic richness is determined from the total number of taxa in the orders EPT53. EPT relative abundance is calculated as the sum of the number of taxa in the orders EPT divided by the total number of organisms in the sample54. The modified FBI considers organisms other than arthropods using genus- and species-level tolerance values and is calculated as follows in Equation (1)55:
where n i = individual number within a taxon, t i = stressor tolerance of a taxon, and N=total number of organisms in the sample. The Diptera and non-insect relative abundance is calculated as the sum of the Diptera taxa and non-insect taxa divided by the total number of organisms in the sample. The relationship between a given biotic index and nutrient was explored via LOWESS, without predefined limits of acceptable ecological conditions (i.e., metrics of invertebrate community composition), and regression tree analysis was used to evaluate the environmental thresholds for TP or TN, where the biotic metric revealed the break point (the greatest degree of change)56. Such analysis was only applied to the taxonomic indices possessing significant relationships with nutrients (TP or TN)16. The basic data statistics were analyzed in Microsoft Excel 2010, and the multivariate statistical analyses were conducted in SPSS 19.0.
No specific permits were required for the described field studies; the sampling did not cause any disturbance to the environment or to the protected species at the sampling sites.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Turner, R. E. & Rabalais, N. N. Coastal eutrophication near the Mississippi river delta. Nature. 368, 619–621 (1994).
ADSArticle Google Scholar2.
Emelko, M. B., Silins, U., Bladon, K. D. & Stone, M. Implications of land disturbance on drinking water treatability in a changing climate: Demonstrating the need for “Source Water Supply and Protection” Strategies. Water. Res. 45, 461–472 (2011).
CASArticlePubMed Google Scholar3.
Dodds, W. K., Jones, J. R. & Welch, E. B. Suggested classification of stream trophic state: distributions of temperate stream types by chlorophyll, total nitrogen, and phosphorus. Water. Res. 32, 1455–1462 (1998).
CASArticle Google Scholar4.
Huang, L. et al. Comparison of phosphorus fractions and phosphatase activities in coastal wetland soils along vegetation zones of Yancheng National Nature Reserve, China. Estuar. Coast. Shelf. S. 157, 93–98 (2015).
ADSCASArticle Google Scholar5.
Hou, L. Y., Hu, A. Y. & Yu, C. P. Evaluation on spatial-temporal distribution of nutrients and potential eutrophication in surface waters of Jiulongjiang River and its Estuaty. J. Appl. Oceanography. (in Chinese). 33, 369–378 (2014).
Zhang, S. et al. Trophic States and nutrient output of tributaries in Three Gorges Reservoir area (in Chinese). Environ. Sci. 28, 500–505 (2007).
Ji, X. L., Dahlgren, R. A. & Zhang, M. H. Comparison of seven water quality assessment methods for the characterization and management of highly impaired river systems. Environ. Monit. Assess. 188, 15 (2016).
ArticlePubMed Google Scholar8.
US EPA. Nutrient Criteria Technical Guidance Manual: Rivers and Streams. (Office of Water & Office of Science and Technology, Washington, DC, USA, 2000).
US EPA. National Strategy for the Development of Regional Nutrient Criteria. (Office ofWater, Washington, DC, USA, 1998).
Harding, L. W. et al. Variable climatic conditions dominate recent phytoplankton dynamics in Chesapeake bay. Sci. Rep-UK. 6, 23773 (2016).
ADSCASArticle Google Scholar11.
Grizzetti, B. et al. Human pressures and ecological status of European rivers. Sci. Rep-UK. 7, 205 (2017).
ADSCASArticle Google Scholar12.
Huang, L., Li, Y., Zhang, Y. & Guan, Y. A simple method to separate phosphorus sorption stages onto solid mediums. Ecol. Eng. 69, 63–69 (2014).
Article Google Scholar13.
Smith, R. A., Alexander, R. B. & Schwarz, G. E. Natural background concentrations of nutrients in streams and rivers of the conterminous United States. Environ. Sci. Technol. 37, 3039–3047 (2003).
ADSCASArticlePubMed Google Scholar14.
Dodds, W. K. & Oakes, R. M. A technique for establishing reference nutrient concentrations across watersheds affected by humans. Limnol. Oceanogr-Meth. 2, 333–341 (2004).
Article Google Scholar15.
Evans-White, M. A., Haggard, B. E. & Scott, J. T. A review of stream nutrient criteria development in the United States. J. Environ. Qual. 42, 1002 (2013).
CASArticlePubMed Google Scholar16.
Chambers, P. A. et al. Development of environmental thresholds for nitrogen and phosphorus in streams. J. Environ. Qual. 41, 7–20 (2012).
CASArticlePubMed Google Scholar17.
Zheng, L., Gerritsen, J., Beckman, J., Ludwig, J. & Wilkes, S. Land use, geology, enrichment, and stream biota in the eastern ridge and valley ecoregion: Implications for nutrient criteria development. J. Am. Water. Resour. As. 44, 1521–1536 (2008).
CASArticle Google Scholar18.
Xia, Q., Chen, Y. Q. & Liu, X. B. Water Quality Criteria and Standard (in Chinese). (Standards Press of China, 2004).
Meng, W., Zhang, Y. & Zheng, B. H. The quality criteria, standards of water environment and the water pollutant control strategy on watershed (in Chinese). Res. Environ. Sci 19, 1–6 (2006).
Wu, F. et al. China embarking on development of its own national water quality criteria system. Environ. Sci. Technol. 44, 7992–7993 (2010).
ADSCASArticlePubMed Google Scholar21.
Chen, J. B. & Lu, J. Establishment of reference conditions for nutrients in an intensive agricultural watershed, Eastern China. Environ. Sci. pollut. R. 21, 2496–2505 (2013).
Article Google Scholar22.
Wu, D. H., Yu, H. Y., Wu, H. Y., Zhou, B. & Wang, B. X. Estimation of river nutrients thresholds based on benthic macroinvertebrate assemblages: a case study in the upper reaches of Xitiao stream in Zhejiang China (in Chinese). Chin. J. Appl. Ecol. 21, 483–488 (2010).
CAS Google Scholar23.
Paul, M. J. & Meyer, J. L. Streams in the urban landscape. Annu. Revecol. Syst. 32, 333–365 (2001).
Article Google Scholar24.
Chen, J., Li, F., Fan, Z. & Wang, Y. Integrated application of multivariate statistical methods to source apportionment of watercourses in the Liao River basin, Northeast China. Int. J. Env. Res. Pub. He. 13, 1035 (2016).
Article Google Scholar25.
Chen, J. & Lu, J. Effects of land use, topography and socio-economic factors on river water quality in a mountainous watershed with intensive agricultural production in East China. PLoS One. 9, e102714 (2014).
ArticlePubMedPubMed Central Google Scholar26.
Chang, H. Spatial analysis of water quality trends in the Han River basin, South Korea. Water. Res. 42, 3285–3304 (2008).
CASArticlePubMed Google Scholar27.
Coulter, C. B., Kolka, R. K. & Thompson, J. A. Water quality in agricultural, urban, and mixed land use watersheds. J. Am. Water. Resour. As. 40, 1593–1601 (2004).
CASArticle Google Scholar28.
Chen, X. et al. Isomeric specific partitioning behaviors of perfluoroalkyl substances in water dissolved phase, suspended particulate matters and sediments in Liao River basin and Taihu Lake, China. Water. Res. 80, 235–244 (2015).
CASArticlePubMed Google Scholar29.
Li, F. Y., Fan, Z. P. & Chen, J. B. Map of Aquatic Ecoregions in the Liao River Basin (in Chinese). (Liaoning Science and Technology Press, 2015).
Zhang, Y., Guo, F., Meng, W. & Wang, X. Water quality assessment and source identification of Daliao River basin using multivariate statistical methods. Environ. Monit. Assess. 152, 105–121 (2009).
CASArticlePubMed Google Scholar31.
Weigel, B. M. & Robertson, D. M. Identifying biotic integrity and water chemistry relations in nonwadeable rivers of Wisconsin: toward the development of nutrient criteria. Environ. Manage. 40, 691–708 (2007).
ADSArticlePubMed Google Scholar32.
Wang, L., Robertson, D. M. & Garrison, P. J. Linkages between nutrients and assemblages of macroinvertebrates and fish in wadeable streams: implication to nutrient criteria development. Environ. Manage. 39, 194–212 (2007).
ADSArticlePubMed Google Scholar33.
Hughes, R. M., Larsen, D. P. & Omernik, J. M. Regional reference sites: a method for assessing stream potentials. Environ. Manage. 10, 629–635 (1986).
ADSCASArticle Google Scholar34.
Suplee, M. W., Varghese, A. & Cleland, J. Developing nutrient criteria for streams: An evaluation of the frequency distribution method. J. Am. Water. Resour. As. 43, 453–472 (2007).
CASArticle Google Scholar35.
Buck, O., Niyogi, D. K. & Townsend, C. R. Scale-dependence of land use effects on water quality of streams in agricultural catchments. Environ. Pollut. 130, 287–299 (2004).
CASArticlePubMed Google Scholar36.
Wu, L., Long, T., Liu, X. & Mmereki, D. Simulation of soil loss processes based on rainfall runoff and the time factor of governance in the Jialing River Watershed, China. Environ. Monit. Assess. 184, 3731–3748 (2012).
ArticlePubMed Google Scholar37.
Yang, Y., Chen, Y., Zhang, X., Ongley, E. & Zhao, L. Methodology for agricultural and rural NPS pollution in a typical county of the North China Plain. Environ. Pollut. 168, 170–176 (2012).
CASArticlePubMed Google Scholar38.
Jacob, J. S. & Lopez, R. Is Denser Greener? An evaluation of higher density development as an urban stormwater-quality best management practice. J. Am. Water. Resour. As. 45, 687–701 (2009).
Article Google Scholar39.
Zhai, S., Yang, L. & Hu, W. Observations of atmospheric nitrogen and phosphorus deposition during the period of algal bloom formation in Northern Lake Taihu, China. Environ. Manage. 44, 542–551 (2009).
ADSArticlePubMed Google Scholar40.
Huang, P. et al. Atmospheric deposition as an important nitrogen load to a typical agroecosystem in the Huang-Huai-Hai Plain. 1. Measurement and preliminary results. Atmos. Environ. 45, 3400–3405 (2011).
ADSCASArticle Google Scholar41.
Mehaffey, M. H. et al. Linking land cover and water quality in New York City’s water supply watersheds. Environ. Monit. Assess. 107, 29–44 (2005).
CASArticlePubMed Google Scholar42.
Ye, L., Cai, Q., Liu, R. & Cao, M. The influence of topography and land use on water quality of Xiangxi River in Three Gorges Reservoir region. Environmen Geol. 58, 937–942 (2009).
ADSCASArticle Google Scholar43.
Anbumozhi, V., Radhakrishnan, J. & Yamaji, E. Impact of riparian buffer zones on water quality and associated management considerations. Ecol. Eng. 24, 517–523 (2005).
Article Google Scholar44.
Maillard, P. & Pinheiro Santos, N. A. A spatial-statistical approach for modeling the effect of non-point source pollution on different water quality parameters in the Velhas river watershed-Brazil. J. Environ. Manage. 86, 158–170 (2008).
CASArticlePubMed Google Scholar45.
Dodds, W. K., Smith, V. H. & Zander, B. Developing nutrient targets to control benthic chlorophyll levels in streams: a case study of the Clark Fork River. Water. Res. 31, 1738–1750 (1997).
CASArticle Google Scholar46.
Li, F. Y. et al. Impacts of land use on the spatial and temporal distribution characteristics of the macroinvertebrates in the Liao River basin (in Chinese). Acta. Sci. Circum. 36, 767–777 (2016).
CAS Google Scholar47.
Li, Y., Hu, C., Li, F. Y. & Huang, Y. G. The application demonstration of technique systems for Qinghe Watershed water quality target management (in Chinese). Environ. Monit. Forew. 5, 1–6 (2013).
MathSciNetCAS Google Scholar48.
Zheng, B. H., Zhang, Y. & Li, Y. B. Study of indicators and methods for river habitat assessment of Liao River basin (in Chinese). Acta. Sci. Circum. 27, 928–936 (2007).
ADSCAS Google Scholar49.
US EPA Office of Water. Rapid bioassessment Protocols for Use in Wadeable Streams and Rivers: Periphyton, Benthic Macroinvertebrates, and Fish. (United States Environmental Protection Agency, 1999).
State Environment Protection Bureau of China. Water and Wastewater Analysis Method (in Chinese). 4th ed. (China Environmental Science Press, 2002).
Morse, J. C., Yang, L. & Tian, L. Aquatic Insects of China Useful for Monitoring Water Quality (in Chinese). (Hohai University Press,1994).
Dalian Fisheries College. Freshwater Biology1: Taxonomic Section (in Chinese). (Agricultural Press, 1982).
Morse, C. C., Huryn, A. D. & Cronan, C. Impervious Surface Area as a Predictor of the Effects of Urbanization On Stream Insect Communities in Maine, USA. Environ. Monit. Assess. 89, 95–127 (2003).
ArticlePubMed Google Scholar54.
Silveira, M. P., Baptista, D. F., Buss, D. F., Nessimian, J. L. & Egler, M. Application of Biological Measures for Stream Integrity Assessment in South-East Brazil. Environ. Monit. Assess. 101(1–3), 117 (2005).
CASPubMed Google Scholar55.
Mandaville, S. M. Benthic Macroinvertebrates in Freshwaters- Taxa Tolerance Values, Metrics, and Protocols (2002).
Moisen, G. G. CART—Classification and Regression Trees. (Wadsworth International Group, 1984).
This work was supported by the National Natural Science Foundation of China (No. 41501551), the Major Science and Technology Program for Water Pollution Control and Treatment (No. 2012ZX07505-001-01), the Research Fund for Scientific Talent of LSHU (No. 00100850) and the Natural Science Foundation of Jiangsu Province (No. BK20150165).
Jiabo Chen and Fayun Li contributed equally to this work.
AffiliationsNational & Local United Engineering Laboratory of Petroleum Chemical Process Operation Optimization and Energy Conservation Technology, Liaoning Shihua University, Fushun, 113001, ChinaJiabo Chen, Fayun Li & Yanjie WangInstitute of Eco-environmental Sciences, Liaoning Shihua University, Fushun, 113001, ChinaJiabo Chen, Fayun Li & Yanjie WangYixing Urban Supervision & Inspection Administration of Product Quality, China National Supervision & Inspection Center of Environmental Protection Equipment Quality (Jiangsu, preparation), Yixing, 214205, ChinaYun Kong
AuthorsSearch for Jiabo Chen in:Nature Research journals • PubMed • Google Scholar Search for Fayun Li in:Nature Research journals • PubMed • Google Scholar Search for Yanjie Wang in:Nature Research journals • PubMed • Google Scholar Search for Yun Kong in:Nature Research journals • PubMed • Google Scholar
F.L., J.C. and Y.W. designed the experiments. F.L., J.C. and Y.W. performed the experiments. Y.K. and Y.W. analyzed the data. J.C. wrote the paper.
The authors declare no competing interests.
Correspondence to Jiabo Chen or Fayun Li.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
26 July 2017
16 February 2018
28 February 2018