Next Article in Journal
Adsorption Studies on Magnetic Nanoparticles Functionalized with Silver to Remove Nitrates from Waters
Previous Article in Journal
Contributions of Vegetation Greening and Climate Change to Evapotranspiration Trend after Large-Scale Vegetation Restoration on the Loess Plateau, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Catchment versus Riparian Buffers: Which Land Use Spatial Scales Have the Greatest Ability to Explain Water Quality Changes in a Typical Temperate Watershed?

1
Beijing Key Laboratory of Traditional Chinese Medicine Protection and Utilization, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
2
State Key Laboratory of Earth Surface Processes and Resource Ecology, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
3
School of Environment, Beijing Normal University, Beijing 100875, China
4
College of Urban and Environmental Sciences, Northwest University, Xi’an 710127, China
*
Author to whom correspondence should be addressed.
Water 2021, 13(13), 1758; https://doi.org/10.3390/w13131758
Submission received: 25 May 2021 / Revised: 19 June 2021 / Accepted: 22 June 2021 / Published: 25 June 2021
(This article belongs to the Section Water Quality and Contamination)

Abstract

:
Identifying the multi-scale spatial relationship between land use and water quality is critical for determining the priorities and key areas of river management. To more accurately identify the scale effect of land-use patterns on water quality and quantitatively distinguish the difference in the impacts of land-use composition and configuration on water quality, we used 94 sites to extract the upstream catchment and riparian buffer zone with different widths. The results showed that the ability of land use variables with different buffer widths to explain water quality differed slightly from the ability of these variables at the catchment scale, and the joint explanatory ability of land use composition and configuration was greater than that of each individually. The patch density and landscape shape index of cultivated land, shrubland, and built-up land in the buffer area close to the water bodies were the main factors for the increase in the concentration of total nitrogen, nitrate nitrogen, total phosphorus, and suspended solids. As the width of the buffer increased, the role of the percent of land use increased. Our research indicates that water quality management needs to adopt a multi-scale perspective and focus on key local areas while coordinating at a broader scale.

Graphical Abstract

1. Introduction

Due to the dual threats of climate change and human activities, most river ecosystems in the world are currently being damaged to varying degrees [1,2,3]. Land use has been proven to be one of the most significant human activities affecting river ecosystems [4,5]. A large amount of industrial and urban domestic wastewater discharge in basins, agricultural production, grassland destruction, deforestation, and other activities have led to ecological and environmental problems such as destruction of river habitats, deterioration of water quality, and degradation of aquatic community structure [6,7,8]. The Global Risks Report 2020 also pointed out that water resources and related environmental problems were one of the most serious current crises and that water quality problems were particularly significant (World Economic Forum, WEF) [9]. Therefore, identifying the relationship between land use and water quality is of great significance to improving the water quality of rivers.
The relationship between land use and water quality is very sensitive to scale changes [10,11,12]. Many studies have shown that the impact of land use on water quality varies from the riparian zone to the catchment or other spatial scales [13,14]. As a land-water ecotone, the riparian zone is the last line of defense for river water quality protection [15,16] and plays an important role in intercepting pollutants, improving water quality, and maintaining habitat diversity for river ecosystems [17,18]. The difference in the relationship between land use and water quality caused by various spatial scales is a challenge for current research [19]. Existing studies disagree over the spatial scales that have the strongest correlation between land use and water quality in a specific watershed [11,20,21], and this disagreement is mainly caused by the uniqueness of each basin [7]. At the same time, confusion among related concepts and irregular definition methods at different scales are also important reasons that it is difficult to compare similar studies [22]. Therefore, the accurate identification of the contribution of land use to multi-scale impacts on water quality still requires further scientific examination.
Previous studies have shown that land-use patterns (mainly in terms of the composition and spatial configuration of land use/cover types) play an important role in the formation, release, absorption, and interception of pollutants, which in turn affect the water quality of rivers [23,24]. Land use configuration refers to the spatial distribution of land use patches and their relationships to one another [25], which affects the exchange of water flows on land surfaces among different land use patches [26]. In recent years, research on the impact of land-use configuration on water quality has gradually attracted attention, and a landscape index has become an important analytical indicator [27]. Studies have shown that utilizing the landscape indices used to describe land use configuration may be better than using the proportion of each land use patch type (land use composition) to predict water quality changes [13,19,28,29]. However, due to the randomness and collinearity of landscape indices selection [30], the difference in the explanatory power of land use composition and configuration in terms of river water quality is uncertain. A lack of understanding of this relationship will restrict the formulation and implementation of specific control measures by watershed managers [31]. What is gratifying is that with the development of multivariate statistical techniques, gradient analysis methods have been introduced to quantitatively distinguish the contribution rates of different land-use variables to the spatial differentiation of water quality [12]. Therefore, it is urgent to further explore how to quantify the impact of land use composition and configuration on water quality and determine the contribution of major land-use variables to the water quality of rivers [11].
The Luanhe River basin in China plays an important role in the “Diversion Project from Luan-he to Tianjin”, several major tributaries of the upper and middle streams of the Luanhe River basin flow into the mainstream, eventually draining into the Daheiting and Panjiakou reservoirs, and the starting point of the Project is Daheiting reservoir, which provides valuable water resources for Tianjin and Tangshan. Great changes in land use have occurred due to urban expansion and agricultural development [32,33], which result in latent dangers to the water supply. Therefore, it is of great significance to explore the multi-scale relationship between land use and water quality in this basin to improve water quality and the healthy development of aquatic ecology. The objectives of this study were threefold: (1) to reveal the relationship between land use patterns and river water quality at different spatial scales; (2) to identify the main land-use variables that affect water quality at different spatial scales; and (3) to distinguish the differential impact of land use composition and spatial configuration on water quality changes.

2. Methods and Materials

2.1. Study Area

The Luanhe River basin originates in the Bayantugur Mountains in Fengning County, Hebei Province, in North China between 115°32′–118°53′ E and 40°11′–42°45′ N and covers an area of 3.6 × 104 km2 (Figure 1). This area is influenced by a typical temperate monsoon climate and has had an average annual atmospheric temperature of 6.1 °C and a mean annual precipitation of 444.1 mm over the last 30 years (http://data.cma.cn, accessed on 16 August 2017). The mean annual average runoff from the study area is 39.3 × 108 m3, approximately 66% of which is from precipitation from June to September [33]. Brown earth, cinnamon soils, and Castanozems account for more than 80% of the total area of the basin, with poor corrosion resistance and high nitrogen content. The predominant land use type is shrubland, which occupied 41.75% of the basin in 2016. Furthermore, its average slope is the largest among all land use types (Table 1). The crops were mainly corn, rice, and soybean, and approximately 3.8 × 106 inhabitants lived in the study area in 2016 [34,35].

2.2. Water Sampling and Parameter Measurements

Considering that the precipitation in the study area is concentrated from June to August and that nonpoint source pollution also mainly occurs during this period, water samples were collected in July 2017 [33]. Based on the principle of stratified sampling and the uniformity of spatial distribution, traffic accessibility, and others, a total of 94 sampling sites were established in the study area, covering various types of land uses and major tributaries in the basin. Water samples were collected based on the “Standard Methods for the examination of Water and Wastewater (version 4), 2002” [36], and 1 L plexiglass water collectors were used to collect water samples at the surface of the river (water depth at the center of the sampling point or 5 m from the bank at approximately 0.5 m). Water samples were collected three times at each sampling site and mixed thoroughly. The mixed water samples were divided into two groups: one was the raw water sample without any treatment; the other water sample was filtered through precombusted 0.45 μm glass fiber filters (WCN grade: Whatman GF/F filters, Maidstone, UK) for subsequent soluble substance analysis. A small amount of concentrated sulfuric acid (H2SO4) was added to the two water samples for acid preservation. The filter was stored in a dark place protected from light and was used to determine chlorophyll a (Chla) and suspended solids (SS). The water sample bottle was immersed in a 1:10 hydrochloric acid solution for 24 h, washed with distilled water, and dried; the glass cellulose film was pretreated at 450 °C for 6 h. All samples were stored frozen at a low temperature and returned to the laboratory for indoor analysis in the shortest possible amount of time.
We analyzed 11 water quality parameters, including pH, electrical conductivity (EC, μs/cm), dissolved oxygen (DO, mg/L), permanganate index (CODMn, mg/L), 5-day biochemical oxygen demand (BOD5, mg/L), total phosphorus (TP, mg/L), total nitrogen (TN, mg/L), ammonia nitrogen (NH3-N, mg/L), nitrate-nitrogen (NO3-N, mg/L), chlorophyll a (Chla, μg/L), and suspended solids (SS, mg/L). The water chemical analysis followed the national standard methods for examining water and wastewater in China [36].

2.3. Quantification of Land Use Indicators

The mainstream and catchment boundaries were extracted from a digital elevation model (DEM, ASTER GDEM V2) with a spatial resolution of 30 m using the ArcSWAT extension module, which is the graphical user interface of the SWAT (Soil and Water Assessment Tool) model on the ArcGIS platform (ESRI Company, Redlands, CA, USA) [37]. The land use data were derived from 2016 Landsat 8 satellite images, downloaded from https://glovis.usgs.gov/ (accessed on 1 August 2018). Images were preprocessed by ENVI 5.3 software (Exelis Visual Information Solutions, Boulder, CO, USA) for radiometric calibration, atmospheric correction, and image sharpening. The processed images were used to create a land-use map through a supervised maximum likelihood classification, and the land uses were classified into 22 types. The overall accuracy of each land use type was above 90%, and the kappa coefficient was above 85%. Land uses were reclassified into seven categories according to the needs of the study: (1) cultivated land (cul), including paddy fields and dry land; (2) woodland (woo), including sparse wood and other forests; (3) shrubland (shr), including bushes and orchards; (4) grassland (gra), including pastures and grass communities; (5) water bodies (wat), including rivers, wetlands, and reservoirs; (6) built-up land (bui), including urban, rural residential areas and other infrastructure; and (7) other land (oth). The proportions of cultivated land, woodland, shrubland, grassland, water bodies, built-up land, and other land in the study area were 17.97%, 16.04%, 41.75%, 17.28%, 0.51%, 3.94%, and 2.50%, respectively (Figure 1).
The land-use variables that were considered to have a great impact on water quality changes in previous studies were selected [38]. We mainly focused on the relationship between land-use variables and water quality parameters at the class level, including the percent of landscape (PLAND), patch density (PD), landscape shape index (LSI), aggregation index (AI), and large patch index (LPI) [24]. In addition, patch density (PD), landscape shape index (LSI), contagion (CONTAG), and the Shannon’s diversity index (SHDI) were calculated at the landscape level to understand the overall land-use characteristics at different scales [12]. It is worth noting that the purpose of our analysis of land-use variables was to link the size, shape, and connectivity of land use patches with specific impacts on ecological water processes rather than just to evaluate and quantify the index in isolation [39]. These metrics were calculated using FRAGSTATS 4.2 software expressed not only the proportion of land use but also the fragmentation of the land use, shape complexity, connectivity, and diversity [40].
Land-use variables were identified at multiple spatial scales (Figure 1), which were the catchment scale (entire upstream catchment above the sampling site) [13] and the scale of a riparian buffer with different widths (100 m, 200 m, 500 m, 1000 m, 1500 m, and 2000 m). The riparian buffer widths were extracted and calculated in batches by Python scripts with the ArcPy site package (ArcGIS 10.3 License, ESRI Company, USA).

2.4. Statistical Analysis

First, the Kolmogorov-Smirnov test was used to examine the normality of the water quality variables, and the log (x + 1) transformation of the water quality variables were performed to meet the assumptions of normality before redundancy analysis (RDA). The land-use variables that greatly contributed to water quality (p < 0.05) were screened through stepwise regression analysis, RDA was performed for the analysis of the water quality and land-use patterns relationships, and Monte Carlo permutation (n = 499) was adopted to test the significance of the RDA results [41]. The cosine value of the angle between the arrow of the land-use variable and the arrow of the water quality variable in the RDA sequence diagram was positive, indicating that the relationship between the two was positively correlated, and the longer the arrow was, the stronger the correlation between the two variables [42]. Furthermore, variation partitioning analysis (VPA) was used to perform partial redundancy analysis on land-use composition and configuration indicators, decomposing the respective individual explanatory ability and joint explanatory ability of the land use composition and configuration for water quality changes [43]. Stepwise regression analysis was implemented in SPSS 20.0 software (IBM Company, Chicago, IL, USA), and both RDA and VPA were implemented in the CANOCO 5.0 program (Microcomputer Power Company, El Segundo, CA, USA).

3. Results

3.1. Differences in Land Use Variables and Water Quality Characteristics

For the land-use composition of the riparian buffer widths and catchment (Figure 2), the proportion of cultivated land and built-up land in the riparian zone showed a clear downward trend as the width of the buffer increased, while the changing trend for the area of woodland and shrubland had the opposite trend. Additionally, among the land use types, shrubland area accounted for the largest proportion and was the dominant land use type in the basin. The proportion of grassland, water bodies, and other land changed relatively little. For the land-use spatial configuration (Figure 3), the PD in the 100 m buffer was the largest. As the width of the buffer increased, the PD gradually decreased, and the PD value of the 2000 m buffer was generally equal to the PD value of the catchment. The LSI and CONTAG both decreased slightly and then gradually increased as the width of the buffer increased, while the SHDI showed the opposite trend. From the overall change trend in land use, it can be seen that the land-use composition and configuration characteristics of the 2000 m buffer and catchment were relatively similar.
There were significant spatial differences in water quality in the study area (Figure 4). The spatial distributions of CODMn, TP, EC, TN, NH3-N, and NO3-N were quite different. Among them, the larger areas of CODMn and NH3-N appeared in the farming-pastoral ecotone in the upper reaches of the basin. The larger areas of TP, EC, TN, and NO3-N were mainly distributed in the hilly areas in the middle and lower reaches of the basin, and the spatial distribution of the other water quality parameters was not significantly different.

3.2. Relationship between Land Use and Water Quality at Multiple Scales

Based on stepwise regression to screen significant explanatory variables (Table 2), the forward selection method in the RDA process was used to filter out the dominant land use variable groups (p < 0.05), and the relationship between land-use variables and water quality was visually shown through an ordination diagram (Figure 5). Table 3 was the explanatory ability of different spatial scales, which showed the land-use variables selected for the catchment and each buffer width explained 29.8% to 35.7% of the changes in water quality. Although the explanatory ability of the 1000 m buffer was the largest (35.7%), the difference in the explanatory ability of the catchment and the 200 m buffer was only 1.9% and 2.5%, respectively.
Figure 5 showed that TN, NO3-N, CODMn, TP, and SS had significant relationships with the land-use variables. Among them, TN and NO3-N were positively correlated with PDcul and PDshr in the 100 m and 200 m buffers, respectively, and were negatively correlated with AIgra. At the 500–2000 m buffer and the catchment scale, TN and NO3-N were mainly positively correlated with PLANDshr and PDwoo. CODMn had a significant positive correlation with AIgra and a negative correlation with PDwoo at the catchment scale and for almost all buffers. TP and SS were mainly positively correlated with LSIshr and PLANDbui at the catchment scale and for almost all buffers. The relationships between other water quality parameters and land-use variables were not obvious.

3.3. Differences in the Capability of Land Use Composition and Configuration to Explain Water Quality at Different Scales

Based on the RDA of the water quality differences explained by the land-use variables, we further used VPA to distinguish the respective explanatory ability of land use composition and configuration. The results (Figure 6) showed that at the catchment scale and for all buffer widths, the shared ability of land use composition and configuration to interpret water quality was greater than that of each individually, and the explanatory ability of land use configuration was greater than that of land use composition.
In terms of the change trends at spatial scales, the individual explanatory ability of the land-use composition and configuration in the 100 m buffer was the largest. With increasing buffer width, the individual explanatory ability of the land-use composition ranged from 0.1% to 1.9%. The range of individual explanatory abilities was 5.8–10.7%, and the overall trend decreased. Additionally, the joint explanatory ability of land use composition and configuration in the 100–1000 m buffer widths showed a gradual increase and remained stable when the buffer exceeded 1000 m.

4. Discussion

4.1. Effective Spatial Scale Identification of Land Use Patterns on Water Quality

Different from conclusions in previous studies [5,7,19], the RDA results showed that there is little difference in the ability to explain water quality changes between the catchment and the buffer zone with different widths (Figure 5). However, it is worth noting that the explanatory rate of the RDA results has its own complex rules, and its ability to analyze the impact of land use and water quality is limited. Specifically, in our study area, in the near-distance riparian buffer widths (e.g., within 100 and 200 m), grassland, built-up land, and cultivated land had a strong ability to explain water quality changes, and shrubland and woodland became the main land use variables to explain water quality changes as the width of the buffer increased. This result was largely due to fertilizer, livestock manure, and domestic sewage from farmlands or urban areas more easily entering the river within these areas [44,45]. The impact of buffer areas far away from water bodies on water quality may be more complicated [13]. The land use in the basin is dominated by shrubland and woodland, their slopes are relatively large, and the corresponding soils are mostly brown earth and cinnamon soils (Table 1). This soil has poor anti-corrosion ability and is more prone to soil erosion, which makes it easier for nitrogen, phosphorus, and other nutrients to enter water bodies through surface runoff [13]. Given this scenario, we should accurately understand the multi-spatial scale relationship between land use and water quality and should not take the explanatory rate as the only basis for the impact of land use on water quality, and we should also pay attention to identifying which land-use variables contribute to the explanatory rate.
Furthermore, the VPA results showed that the ability of land use composition and configuration to jointly explain water quality changes at different scales was greater than the individual ability of these two factors (Figure 6); this result indicates that there was collinearity between the land use composition and configuration and that the joint influence of the two was an important reason for the variation in the water quality, which is consistent with existing views [13]. Therefore, we need to understand the impact of land-use composition and configuration at different spatial scales on water quality changes from a holistic perspective instead of focusing at a certain scale in isolation.

4.2. Main Land-Use Variables That Affect Water Quality at Different Spatial Scales

The RDA results show that TN and NO3-N were positively correlated with PDcul and PDshr in the riparian buffer zone within 200 m of the water bodies, and this result may have been related to agricultural nonpoint source pollution [46]. We know that summer is a critical period for crop growth, when the frequency of fertilization increases and the surface runoff is relatively high, which makes it easier for nutrients such as nitrogen and phosphorus to migrate into water bodies [47]. In addition, the fragmentation of forests is not conducive to the filtration and retention of pollutants [48], which is also an important reason for this correlation. It is worth noting that there was a very significant positive correlation between TN, NO3-N, and PLANDshr in the buffer beyond 200 m and at the catchment scale. This result is different from the conclusion that forestland is generally regarded as a medium that can intercept pollutants and improve water quality [49]. The reason for this result is that on the one hand, it may be related to the classification of land use types [14]. When we interpreted land use types, the orchards and traditional Chinese medicine plantations may have been classified into the shrubland category due to the constraints on satellite image resolution. Fertilization during the planting of orchards and Chinese medicinal materials increases the nitrogen content of the soil. On the other hand, the average slope of the shrubland was the largest among all land use types (Table 1), which was conducive to soil erosion after rainfall and accelerated the flow of fertilizer into the river [50].
In addition, TN and NO3-N were negatively correlated with AIgra in the buffer close to the water bodies, which was generally confirmed by existing studies that show that grassland has a purifying effect on water quality. In addition, CODMn was positively correlated with AIgra at all scales, which may have been caused by the high humus in the grassland on the Bashang Plateau [33] and may also have been related to the enhanced control of dominant land-use types (grassland) in the overall landscape [40].
TP and SS were positively correlated with PLANDbui at almost all scales, which is consistent with the positive correlation between built-up land and river water quality degradation [47,51,52]. The built-up land in the study area is mainly distributed near water, and the lack of effective treatment of the large amount of surface debris produced by resident life increases the concentration of TP and SS [53]. Additionally, the increase in impervious surface area (ISA) in the process of urbanization made it impossible for natural vegetation and soil to filter and intercept pollutants [54].

4.3. Implications for Management

Improving river water quality is a complex and long process that requires continuous control and intervention by managers [31,55]. In contrast to the results of some studies that focus on land use management, our results mainly focus on improving water quality by rationally optimizing the composition and configuration of land use at different spatial scales [56]. This approach calls for managers to focus on planning the configuration of natural forests and grasslands in the riparian zone in a scientific way to reduce vegetation fragmentation [57,58]. Furthermore, attention must be paid to the control of agricultural nonpoint source pollution [59], and cultivated land should be moved away from the buffer zones close to water bodies under the premise of ensuring that the total area of the cultivated land in the region remains unchanged [60,61].
At the catchment scale and for the buffer areas far away from the water bodies, especially in hilly areas with large slopes, more attention should be given to protecting natural woodlands and shrublands and increasing the ratio of woodland and grassland area to ensure the healthy development of the aquatic environment.

5. Conclusions

Based on the results of the relationship between river water quality and land-use variables, we found that there was little difference in the ability of the riparian zone and the catchment scale to explain water quality changes. Cultivated land, built-up land, and shrubland were the main factors that influenced the increase in the concentrations of TN, NO3-N, TP, and SS in the river. The good connectivity of the grasslands accounted for the high CODMn. The joint explanatory ability of land use composition and configuration was greater than their separate explanatory ability, and the closer to the water bodies, the more obvious the land-use configuration role was. With the increase in the width of the buffer, the impact of land-use composition on water quality gradually increased.
Our research results support that water quality management is a comprehensive problem that requires overall consideration and attention in local areas. The results provide a scientific basis for river water ecological management. Although nonpoint source pollution is generally more serious in summer than in other seasons, only sampling water quality in summer may have concealed some information caused by seasonal differences. In future work, we will collect water quality data in different seasons to further explore whether seasonal changes in the study area will affect the relationship between land use and water quality.

Author Contributions

Conceptualization, M.S. and Q.L.; Methodology, M.S. and M.K.; Software, M.S.; Validation, M.S. and Y.J.; Formal analysis, M.S.; Investigation, M.S., Y.T., Y.L. and X.X.; Data Curation, M.S. and Y.J.; Writing—Original Draft Preparation, M.S.; Writing—Review & Editing, M.S.; Visualization, M.S.; Supervision, Y.J.; Project Administration, Y.J.; Funding Acquisition, Y.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science and Technology Major Project: Water Pollution Control and Management Technology of China, grant number No. 2017ZX07301-001-03 and the APC was funded by the National Science and Technology Major Project: Water Pollution Control and Management Technology of China, grant number No. 2017ZX07301-001-03.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

This work was supported by the National Science and Technology Major Project: Water Pollution Control and Management Technology of China (No. 2017ZX07301-001-03). We would like to thanks the anonymous reviewers for their contributions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vörösmarty, C.J.; McIntyre, P.B.; Gessner, M.O.; Dudgeon, D.; Prusevich, A.; Green, P.; Davies, P.M. Global threats to human water security and river biodiversity. Nature 2010, 467, 555–561. [Google Scholar] [CrossRef]
  2. Feld, C.K.; Fernandes, M.R.; Ferreira, M.T.; Hering, D.; Ormerod, S.J.; Venohr, M.; Gutiérrez-Cánovas, C. Evaluating riparian solutions to multiple stressor problems in river ecosystems—a conceptual study. Water Res. 2018, 139, 381–394. [Google Scholar] [CrossRef] [PubMed]
  3. Liu, Q.; Tian, Y.; Liu, Y.; Yu, M.; Hou, Z.; He, K.; Jiang, Y. Relationship between dissolved organic matter and phytoplankton community dynamics in a human-impacted subtropical river. J. Clean. Prod. 2021, 289, 125–144. [Google Scholar] [CrossRef]
  4. Tu, J. Spatially varying relationships between land use and water quality across an urbanization gradient explored by geographically weighted regression. Appl. Geogr. 2011, 31, 376–392. [Google Scholar] [CrossRef]
  5. Mainali, J.; Chang, H. Landscape and anthropogenic factors affecting spatial patterns of water quality trends in a large river basin, South Korea. J. Hydrol. 2018, 564, 26–40. [Google Scholar] [CrossRef]
  6. Barbosa, A.E.; Fernandes, J.N.; David, L.M. Key issues for sustainable urban stormwater management. Water Res. 2012, 46, 6787–6798. [Google Scholar] [CrossRef]
  7. Sun, Y.; Guo, Q.; Liu, J.; Wang, R. Scale effects on spatially varying relationships between urban landscape patterns and water quality. Environ. Manag. 2014, 54, 272–287. [Google Scholar] [CrossRef]
  8. Nafi’Shehab, Z.; Jamil, N.R.; Aris, A.Z.; Shafie, N.S. Spatial variation impact of landscape patterns and land use on water quality across an urbanized watershed in Bentong, Malaysia. Ecol. Indic. 2021, 122, 107254. [Google Scholar] [CrossRef]
  9. Franco, E.G. The Global Risks Report 2020. In World Economic Forum. 2020. Available online: https://www.weforum.org/reports/the-global-risks-report-2020 (accessed on 28 January 2020).
  10. Sliva, L.; Williams, D.D. Buffer zone versus whole catchment approaches to studying land use impact on river water quality. Water Res. 2001, 35, 3462–3472. [Google Scholar] [CrossRef]
  11. Shi, P.; Zhang, Y.; Li, Z.; Li, P.; Xu, G. Influence of land use and land cover patterns on seasonal water quality at multi-spatial scales. Catena 2017, 151, 182–190. [Google Scholar] [CrossRef]
  12. Hille, S.; Andersen, D.K.; Kronvang, B.; Baattrup-Pedersen, A. Structural and functional characteristics of buffer strip vegetation in an agricultural landscape–high potential for nutrient removal but low potential for plant biodiversity. Sci. Total Environ. 2018, 628, 805–814. [Google Scholar] [CrossRef]
  13. Ding, J.; Jiang, Y.; Liu, Q.; Hou, Z.; Liao, J.; Fu, L.; Peng, Q. Influences of the land use pattern on water quality in low-order streams of the Dongjiang River basin, China: A multi-scale analysis. Sci. Total Environ. 2016, 551, 205–216. [Google Scholar] [CrossRef] [PubMed]
  14. Xu, Q.; Wang, P.; Shu, W.; Ding, M.; Zhang, H. Influence of landscape structures on river water quality at multiple spatial scales: A case study of the Yuan river watershed, China. Ecol. Indic. 2021, 121, 107226. [Google Scholar] [CrossRef]
  15. Amuchástegui, G.; Di Franco, L.; Feijoó, C. Catchment morphometric characteristics, land use and water chemistry in Pampean streams: A regional approach. Hydrobiologia 2016, 767, 65–79. [Google Scholar] [CrossRef]
  16. Leal, C.G.; Barlow, J.; Gardner, T.A.; Hughes, R.M.; Leitão, R.P.; Mac Nally, R.; Pompeu, P.S. Is environmental legislation conserving tropical stream faunas? A large-scale assessment of local, riparian and catchment-scale influences on Amazonian fish. J. Appl. Ecol. 2018, 55, 1312–1326. [Google Scholar] [CrossRef] [PubMed]
  17. Meixler, M.S.; Bain, M.B. Landscape scale assessment of stream channel and riparian habitat restoration needs. Landsc. Ecol. Eng. 2010, 6, 235–245. [Google Scholar] [CrossRef]
  18. Dosskey, M.G.; Vidon, P.; Gurwick, N.P.; Allan, C.J.; Duval, T.P.; Lowrance, R. The role of riparian vegetation in protecting and improving chemical water quality in streams 1. J. Am. Water Resour. Assoc. 2010, 46, 261–277. [Google Scholar] [CrossRef]
  19. Li, K.; Chi, G.; Wang, L.; Xie, Y.; Wang, X.; Fan, Z. Identifying the critical riparian buffer zone with the strongest linkage between landscape characteristics and surface water quality. Ecol. Indic. 2018, 93, 741–752. [Google Scholar] [CrossRef]
  20. Xiao, R.; Wang, G.; Zhang, Q.; Zhang, Z. Multi-scale analysis of relationship between landscape pattern and urban river water quality in different seasons. Sci. Rep. 2016, 6, 25250. [Google Scholar] [CrossRef] [Green Version]
  21. Zhang, J.; Li, S.; Dong, R.; Jiang, C.; Ni, M. Influences of land use metrics at multi-spatial scales on seasonal water quality: A case study of river systems in the Three Gorges Reservoir Area, China. J. Clean. Prod. 2019, 206, 76–85. [Google Scholar] [CrossRef]
  22. Gove, N.E.; Edwards, R.T.; Conquest, L.L. Effects of scale on land use and water quality relationships: A longitudinal basin-wide perspective 1. J. Am. Water Resour. Assoc. 2001, 37, 1721–1734. [Google Scholar] [CrossRef]
  23. Mitchell, M.G.; Bennett, E.M.; Gonzalez, A. Linking landscape connectivity and ecosystem service provision: Current knowledge and research gaps. Ecosystems 2013, 16, 894–908. [Google Scholar] [CrossRef]
  24. Clément, F.; Ruiz, J.; Rodríguez, M.A.; Blais, D.; Campeau, S. Landscape diversity and forest edge density regulate stream water quality in agricultural catchments. Ecol. Indic. 2017, 72, 627–639. [Google Scholar] [CrossRef]
  25. McGarigal, K.; Marks, B.J. Spatial Pattern Analysis Program for Quantifying Landscape Structure; Gen. Tech. Rep. PNW-GTR-351; US Department of Agriculture, Forest Service, Pacific Northwest Research Station: Portland, OR, USA, 1995; pp. 1–122.
  26. Ricart, A.M.; Dalmau, A.; Pérez, M.; Romero, J. Effects of landscape configuration on the exchange of materials in seagrass ecosystems. Mar. Ecol. Prog. Ser. 2015, 532, 89–100. [Google Scholar] [CrossRef]
  27. Shen, Z.; Hou, X.; Li, W.; Aini, G. Relating landscape characteristics to non-point source pollution in a typical urbanized watershed in the municipality of Beijing. Landsc. Urban Plan. 2014, 123, 96–107. [Google Scholar] [CrossRef]
  28. Uuemaa, E.; Antrop, M.; Roosaare, J.; Marja, R.; Mander, Ü. Landscape metrics and indices: An overview of their use in landscape research. Living Rev. Landsc. Res. 2009, 3, 1–28. Available online: http://hdl.handle.net/1854/LU-695518 (accessed on 9 June 2018). [CrossRef]
  29. Uuemaa, E.; Mander, Ü.; Marja, R. Trends in the use of landscape spatial metrics as landscape indicators: A review. Ecol. Indic. 2013, 28, 100–106. [Google Scholar] [CrossRef]
  30. Zhang, G.; Guhathakurta, S.; Dai, G.; Wu, L.; Yan, L. The control of land-use patterns for stormwater management at multiple spatial scales. Environ. Manag. 2013, 51, 555–570. [Google Scholar] [CrossRef]
  31. Liu, J.; Liu, X.; Wang, Y.; Li, Y.; Jiang, Y.; Fu, Y.; Wu, J. Landscape composition or configuration: Which contributes more to catchment hydrological flows and variations? Landsc. Ecol. 2020, 35, 1531–1551. [Google Scholar] [CrossRef]
  32. Cheng, X.; Chen, L.; Sun, R. Modeling the non-point source pollution risks by combing pollutant sources, precipitation, and landscape structure. Environ. Sci. Pollut. R. 2019, 26, 11856–11863. [Google Scholar] [CrossRef]
  33. Tian, Y.; Jiang, Y.; Liu, Q.; Dong, M.; Xu, D.; Liu, Y.; Xu, X. Using a water quality index to assess the water quality of the upper and middle streams of the Luanhe River, northern China. Sci. Total Environ. 2019, 667, 142–151. [Google Scholar] [CrossRef]
  34. Hebei Province Bureau of Statistics. Hebei Province Statistical Yearbook 2017; China Statistics Press: Beijing, China, 2017. [Google Scholar]
  35. Inner Mongolia Autonomous Region Bureau of Statistics. Inner Mongolia Autonomous Region Statistical Yearbook 2017; China Statistics Press: Beijing, China, 2017. [Google Scholar]
  36. State Environmental Protection Administration. The National Standards of the People’s Republic of China: Environmental Quality Standards for Surface Water (GB 3838-2002); Chinese Environmental Sciences Press: Beijing, China, 2002.
  37. Arnold, J.G.; Srinivasan, R.; Muttiah, R.S.; Williams, J.R. Large area hydrologic modeling and assessment part I: Model development 1. J. Am. Water Resour. Assoc. 1998, 34, 73–89. [Google Scholar] [CrossRef]
  38. Alberti, M.; Booth, D.; Hill, K.; Coburn, B.; Avolio, C.; Coe, S.; Spirandelli, D. The impact of urban patterns on aquatic ecosystems: An empirical analysis in Puget lowland sub-basins. Landsc. Urban Plan. 2007, 80, 345–361. [Google Scholar] [CrossRef]
  39. Kupfer, J.A. Landscape ecology and biogeography: Rethinking landscape metrics in a post-FRAGSTATS landscape. Prog. Phys. Geogr. 2012, 36, 400–420. [Google Scholar] [CrossRef]
  40. McGarigal, K.; Cushman, S.A.; Ene, E. FRAGSTATS v4: Spatial Pattern Analysis Program for Categorical and Continuous Maps. Computer Software Program Produced by the Authors at the University of Massachusetts, Amherst. 2012. Available online: http://www. umass. edu/landeco/research/fragstats/fragstats.html (accessed on 8 June 2018).
  41. Braak, C.J.; Smilauer, P. Canoco Reference Manual and User’s Guide: Software for Ordination, Version 5.0. Ithaca USA: Microcomputer Power; Wageningen University: Wageningen, The Netherlands, 2012; p. 496. [Google Scholar]
  42. Jongman, E.; Jongman, S.R.R. Data Analysis in Community and Landscape Ecology; Cambridge University Press: Cambridge, UK, 1995; Available online: www.cambridge.org/9780521475747 (accessed on 8 June 2018).
  43. Økland, R.H.; Eilertsen, O. Canonical correspondence analysis with variation partitioning: Some comments and an application. J. Veg. Sci. 1994, 5, 117–126. [Google Scholar] [CrossRef]
  44. de Mello, K.; Valente, R.A.; Randhir, T.O.; dos Santos, A.C.A.; Vettorazzi, C.A. Effects of land use and land cover on water quality of low-order streams in Southeastern Brazil: Watershed versus riparian zone. Catena 2018, 167, 130–138. [Google Scholar] [CrossRef]
  45. Bawa, R.; Dwivedi, P. Impact of land cover on groundwater quality in the Upper Floridan Aquifer in Florida, United States. Environ. Pollut. 2019, 252, 1828–1840. [Google Scholar] [CrossRef] [PubMed]
  46. Guo, Q.; Ma, K.; Yang, L.; He, K. Testing a Dynamic Complex Hypothesis in the Analysis of Land Use Impact on Lake Water Quality. Water Resour. Manag. 2010, 24, 1313–1332. [Google Scholar] [CrossRef] [Green Version]
  47. Ouyang, W.; Skidmore, A.K.; Toxopeus, A.G.; Hao, F. Long-term vegetation landscape pattern with non-point source nutrient pollution in upper stream of Yellow River basin. J. Hydrol. 2010, 389, 373–380. [Google Scholar] [CrossRef]
  48. Haidary, A.; Amiri, B.J.; Adamowski, J.; Fohrer, N.; Nakane, K. Assessing the impacts of four land use types on the water quality of wetlands in Japan. Water Resour. Manag. 2013, 27, 2217–2229. [Google Scholar] [CrossRef]
  49. Pacheco, F.A.L.; Santos, R.M.B.; Fernandes, L.S.; Pereira, M.G.; Cortes, R.M.V. Controls and forecasts of nitrate yields in forested watersheds: A view over mainland Portugal. Sci. Total Environ. 2015, 537, 421–440. [Google Scholar] [CrossRef]
  50. Ahearn, D.S.; Sheibley, R.W.; Dahlgren, R.A.; Anderson, M.; Johnson, J.; Tate, K.W. Land use and land cover influence on water quality in the last free-flowing river draining the western Sierra Nevada, California. J. Hydrol. 2005, 313, 234–247. [Google Scholar] [CrossRef]
  51. Lee, S.W.; Hwang, S.J.; Lee, S.B.; Hwang, H.S.; Sung, H.C. Landscape ecological approach to the relationships of land use patterns in watersheds to water quality characteristics. Landsc. Urban Plan. 2009, 92, 80–89. [Google Scholar] [CrossRef]
  52. Kennen, J.G.; Riva-Murray, K.; Beaulieu, K.M. Determining hydrologic factors that influence stream macroinvertebrate assemblages in the northeastern US. Ecohydrol. Ecosyst. Land Water Process Interact. Ecohydrogeomorphology 2010, 3, 88–106. [Google Scholar] [CrossRef]
  53. Namugize, J.N.; Jewitt, G.; Graham, M. Effects of land use and land cover changes on water quality in the uMngeni river catchment, South Africa. Phys. Chem. Earth Parts A B C 2018, 105, 247–264. [Google Scholar] [CrossRef]
  54. Jackson, C.R.; Bahn, R.A.; Webster, J.R. Water quality signals from rural land use and exurbanization in a mountain landscape: What’s clear and what’s confounded? J. Am. Water Resour. Assoc. 2017, 53, 1212–1228. [Google Scholar] [CrossRef]
  55. Ding, J.; Jiang, Y.; Fu, L.; Liu, Q.; Peng, Q.; Kang, M. Impacts of land use on surface water quality in a subtropical River Basin: A case study of the Dongjiang River Basin, Southeastern China. Water 2015, 7, 4427–4445. [Google Scholar] [CrossRef] [Green Version]
  56. Tudesque, L.; Tisseuil, C.; Lek, S. Scale-dependent effects of land cover on water physico-chemistry and diatom-based metrics in a major river system, the Adour-Garonne basin (South Western France). Sci. Total Environ. 2014, 466, 47–55. [Google Scholar] [CrossRef]
  57. Amiri, B.J.; Nakane, K. Modeling the linkage between river water quality and landscape metrics in the Chugoku District of Japan. Water Resour. Manag. 2009, 23, 931–956. [Google Scholar] [CrossRef]
  58. Ai, L.; Shi, Z.H.; Yin, W.; Huang, X. Spatial and seasonal patterns in stream water contamination across mountainous watersheds: Linkage with landscape characteristics. J. Hydrol. 2015, 523, 398–408. [Google Scholar] [CrossRef]
  59. Mokarram, M.; Saber, A.; Sheykhi, V. Effects of heavy metal contamination on river water quality due to release of industrial effluents. J. Clean. Prod. 2020, 277, 123380. [Google Scholar] [CrossRef]
  60. Fernandes, J.D.F.; de Souza, A.L.; Tanaka, M.O. Can the structure of a riparian forest remnant influence stream water quality? A tropical case study. Hydrobiologia 2014, 724, 175–185. [Google Scholar] [CrossRef]
  61. Gu, Q.; Hu, H.; Ma, L.; Sheng, L.; Yang, S.; Zhang, X.; Chen, L. Characterizing the spatial variations of the relationship between land use and surface water quality using self-organizing map approach. Ecol. Indic. 2019, 102, 633–643. [Google Scholar] [CrossRef]
Figure 1. Location and characteristics of the study area: (a) the location of the 94 sampling sites in the Luanhe River basin and the land use types; (b) catchment (entire upstream catchment above the sampling site); and (c) riparian buffer with different widths.
Figure 1. Location and characteristics of the study area: (a) the location of the 94 sampling sites in the Luanhe River basin and the land use types; (b) catchment (entire upstream catchment above the sampling site); and (c) riparian buffer with different widths.
Water 13 01758 g001
Figure 2. The variations in the land use composition of the riparian buffer widths and catchment. Land use types include the cultivated land (cul), woodland (woo), shrubland (shr), grassland (gra), water bodies (wat), built-up land (bui) and other land (oth).
Figure 2. The variations in the land use composition of the riparian buffer widths and catchment. Land use types include the cultivated land (cul), woodland (woo), shrubland (shr), grassland (gra), water bodies (wat), built-up land (bui) and other land (oth).
Water 13 01758 g002
Figure 3. Land use metrics representing the spatial configurations of the riparian buffer widths and catchment: (a) PD means the number of patches per unit area (number per 100 ha); (b) LSI provides a standardized measure of total edge or edge density that adjusts for the size of the landscape; (c) CONTAG describes the degree of agglomeration or spreading trend of different patch types in the landscape; (d) SHDI equals minus the sum, across all patch types, of the proportional abundance of each patch type multiplied by that proportion.
Figure 3. Land use metrics representing the spatial configurations of the riparian buffer widths and catchment: (a) PD means the number of patches per unit area (number per 100 ha); (b) LSI provides a standardized measure of total edge or edge density that adjusts for the size of the landscape; (c) CONTAG describes the degree of agglomeration or spreading trend of different patch types in the landscape; (d) SHDI equals minus the sum, across all patch types, of the proportional abundance of each patch type multiplied by that proportion.
Water 13 01758 g003aWater 13 01758 g003b
Figure 4. The spatial variations in the water quality parameters in the Luanhe River basin.
Figure 4. The spatial variations in the water quality parameters in the Luanhe River basin.
Water 13 01758 g004
Figure 5. Redundancy analysis (RDA) biplots showing the correlation between water quality variables (black arrows) and land use variables (red arrows). (af) respectively represent the scale of a riparian buffer with different widths (100 m, 200 m, 500 m, 1000 m, 1500 m, and 2000 m), and (g) represents catchment scale.
Figure 5. Redundancy analysis (RDA) biplots showing the correlation between water quality variables (black arrows) and land use variables (red arrows). (af) respectively represent the scale of a riparian buffer with different widths (100 m, 200 m, 500 m, 1000 m, 1500 m, and 2000 m), and (g) represents catchment scale.
Water 13 01758 g005aWater 13 01758 g005b
Figure 6. Fractions of the variation explained by individual and joint effects of the land use composition and configuration.
Figure 6. Fractions of the variation explained by individual and joint effects of the land use composition and configuration.
Water 13 01758 g006
Table 1. Average slope of main land use types.
Table 1. Average slope of main land use types.
Land Use TypesCultivated LandWoodLandShrubLandGrassLandBuilt-Up Land
Average slope (°)9.2019.0719.598.669.55
Table 2. Significant explanatory variables screened out by stepwise regression analysis.
Table 2. Significant explanatory variables screened out by stepwise regression analysis.
ScalesExplanatory Variables
Riparian buffer width100 mPLANDcul, PDcul, LSIcul, PDwoo, PLANDshr, LSIshr, PLANDgra, AIgra, LSIgra, PDgra, PLANDbui
200 mPLANDcul, LPIcul, PDwoo, AIwoo, AIshr, PDshr, AIgra, LSIgra, LSIbui, PLANDbui, PDbui
500 mLSIcul, PDwoo, PLANDshr, LSIshr, LPIgra, LSIgra, AIgra, AIbui, PLANDbui
1000 mPLANDcul, LSIcul, LSIwoo, PLANDwoo, PDwoo, PLANDshr, AIshr, LPIgra, AIgra, AIbui, PLANDbui
1500 mLSIcul, LSIwoo, PDwoo, PLANDshr, LSIshr, AIgra, LPIgra, PLANDbui, AIbui, PDbui
2000 mLSIcul, PLANDwoo, PDwoo, PLANDshr, LSIgra, AIgra, PDgra, PLANDbui
CatchmentLSIcul, PLANDwoo, LSIwoo, PLANDshr, AIgra, PDgra, LSIgra, PDbui, PLANDbui
Table 3. RDA ranking axes and the total variance explanatory rate at each buffer width and catchment scale.
Table 3. RDA ranking axes and the total variance explanatory rate at each buffer width and catchment scale.
ScalesAxis 1Axis 2Axis 3Axis 4Explained Variance (%)
Riparian buffer width100 m
EV0.2050.0760.0130.00229.8
CPC (%)68.8794.4798.8799.55
200 m
EV0.2120.0860.0110.00333.8
CPC (%)67.5695.1798.7599.58
500 m
EV0.2220.0970.00090.00533.5
CPC (%)66.1895.2698.0299.54
1000 m
EV0.2340.1090.0090.00435.7
CPC (%)65.4595.8998.4699.56
1500 m
EV0.2330.0980.0110.01135.6
CPC (%)65.5392.9096.0999.19
2000 m
EV0.2200.0900.0190.00333.4
CPC (%)65.8292.8198.5299.54
CatchmentEV0.2120.0970.0170.00533.2
CPC (%)63.8092.8598.0299.57
Notes: “EV” is the abbreviation for “Eigen values”, “CPC (%)” is the abbreviation for “Cumulative percentage correlation of land use-water quality variable (%)”.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Song, M.; Jiang, Y.; Liu, Q.; Tian, Y.; Liu, Y.; Xu, X.; Kang, M. Catchment versus Riparian Buffers: Which Land Use Spatial Scales Have the Greatest Ability to Explain Water Quality Changes in a Typical Temperate Watershed? Water 2021, 13, 1758. https://doi.org/10.3390/w13131758

AMA Style

Song M, Jiang Y, Liu Q, Tian Y, Liu Y, Xu X, Kang M. Catchment versus Riparian Buffers: Which Land Use Spatial Scales Have the Greatest Ability to Explain Water Quality Changes in a Typical Temperate Watershed? Water. 2021; 13(13):1758. https://doi.org/10.3390/w13131758

Chicago/Turabian Style

Song, Minmin, Yuan Jiang, Qi Liu, Yulu Tian, Yang Liu, Xia Xu, and Muyi Kang. 2021. "Catchment versus Riparian Buffers: Which Land Use Spatial Scales Have the Greatest Ability to Explain Water Quality Changes in a Typical Temperate Watershed?" Water 13, no. 13: 1758. https://doi.org/10.3390/w13131758

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