Relating Land Use/Cover and Landscape Pattern to the Water Quality under the Simulation of SWAT in a Reservoir Basin, Southeast China

: Understanding the relationship between land use/cover pattern and water quality could provide guidelines for non-point source pollution and facilitate sustainable development. The previous studies mainly relate the land use/cover of the entire region to the water quality at the monitoring sites, but the water quality at monitoring sites did not totally reﬂect the water environment of the entire basin. In this study, the land use/cover was monitored on Google Earth Engine in Tang-Pu Reservoir basin, China. In order to reﬂect the water quality of the whole study area, the spatial distribution of the determinants for water quality there, i.e., the total nitrogen and total phosphorus (TN&TP), were simulated by the Soil and Water Assessment Tool (SWAT). The redundancy analysis explored the correlations between land use/cover pattern and simulated TN&TP. The results showed that: (1) From 2009 to 2019, forest was the dominant land cover, and there was little land use/cover change. The landscape fragmentation increased, and the connectivity decreased. (2) About 25% TP concentrations and nearly all the TN concentrations at the monitoring points did not reach drinking water standard, which means nitrogen and phosphorus pollution were the most serious problems. The highest output per unit TN&TP simulated by SWAT were 44.50 kg/hm 2 and 9.51 kg/hm 2 and occurred in areas with highly fragile landscape patterns. (3) TN&TP correlated positively with cultivated and construction land but negatively with forest. The correlation between forest and TN&TP summited at 500–700-m buffer and construction land at 100-m buffer. As the buffer size increased, the correlation between the cultivated land, and the TN weakened, while the correlation with the TP increased. TN&TP correlated positively with the Shannon’s Diversity Index and negatively with the Contagion Index. This study provides a new perspective for exporting the impact of land use/cover pattern on water quality.


Introduction
Maintaining clean water exerts a pivotal function on habitat protection, agricultural production, and public health. However, it is challenging, especially due to the non-point source (NPS) pollution [1][2][3][4]. NPS pollution is mainly composed of soil sediment particles, nitrogen and phosphorus, and other nutrients, pesticides, etc., which enter the water, soil, or atmospheric environment through surface runoff, soil erosion, farmland drainage, etc. [5]. The rapid spread of urbanization, severe soil erosion, and excessive application of pesticides and fertilizers all have significant contributions to the NPS pollution of the regional water environment-an important environmental problem that must be resolved urgently [6,7].
The influence of land use/cover type on water quality has been a concern since the 1970s [8]. Previous studies have shown that the deterioration of water quality is driven by improper land use/cover due to anthropogenic activities [9,10]. Inappropriate land use/cover has an impact on the hydrological process of the basin by causing changes in the properties of the underlying surface and soil composition [11,12]. It can affect the process of runoff generation in the basin, further affecting the migration and transformation of NPS pollutants, especially nitrogen and phosphorus [5,13]. Castillo et al. [9] and Rodrigues et al. [14] noted that land use/cover type can increase nutrient loading and reduce nutrient retention, resulting in the eutrophication and hypoxia of the ecosystem. Other studies have found that the artificial regulation of land use/cover can control NPS pollution, thereby effectively supporting ecological functions [15]. Generally, studies have mainly focused on the impact of land use/cover type on water chemistry at different temporal and spatial scales. Mello et al. [16] collected samples to evaluate the land use/cover effects on water quality using Pearson's correlation and redundancy analysis (RDA). Their results indicated that forest cover played a significant role in keeping the water clean, whereas agriculture and urban areas led to water quality degradation and placed immense pressure on the ecosystem. Based on Bayesian Networks, Xu et al. [17] established 20 monitoring stations and studied the response of water quality to land use/cover. They found that controlling water contaminant sources at a less than 20-km radii around monitoring stations could improve water quality efficiently, with farmlands acting as a "source" in wet seasons and as a "sink" in dry seasons.
Land use/cover composition is generally related to water quality, and landscape configuration controlling various biogeochemical and physical processes of the basin may be a more sensitive predictor of water quality [18][19][20]. There are many studies revealing that landscape metrics were significantly associated with water quality. Ding et al. [21] established 56 sampling sites and chose nine water parameters to reflect the water quality. They developed empirical models relating landscape metrics to the water quality using multivariate statistical analyses. The results illustrated that clean water quality was linked to the lower largest patch index of constructed land in plain catchments, whereas it was linked to a lower value for the patch densities of farmland, grassland, and orchard mountain watershed. Based on spatial filter regression of Moran's Eigenvector, Janardan et al. [22] revealed that different land use/cover metrics had different scale effects, and water quality displayed highly temporal variations. At present, many studies have established sampling sites to measure water quality parameters so as to explore the impact of land use/cover and landscape pattern on water quality. However, the water quality at several sampling sites does not necessarily reflect the water quality of the entire basin. Thus, to better evaluate the complex impacts of land use/cover type and landscape pattern on water quality, it is necessary to quantify the correlation between land use/cover pattern and the water quality of the entire basin.
The Soil and Water Assessment Tool (SWAT) is a continuous, long-term, physically based distributed model developed to assess the impacts of climate and land management on hydrological components, sediment loading, and pollution transport in basins [23]. It has a wide range of applications and can be used for NPS pollution analysis [24,25]. When there is a lack of data, the SWAT model can be established by applying model parameter transplantation parameter transfer [26,27]. Zeiger et al. [28] used SWAT to simulate land use/cover impacts on streamflow, sediment, and nutrient loading in a karst catchment. Sertel et al. [25] evaluated the impact of land cover/use changes on watershed responses and hydrological processes in a rapidly urbanizing mid-latitude water supply catchment by applying SWAT. At this stage, many studies have focused on setting up a SWAT model and simulating runoff, sediment, and nutrients by land use/cover data. In addition, the NPS pollution scenarios based on the change of land use/cover have been studied extensively [25,29,30]. Rarely, studies established a correlation between land use/cover and water quality simulated by SWAT model.
It is crucial to investigate the relationship between land use/cover pattern and water quality, as this is conducive to optimizing land use/cover patterns and improving waterbased ecological environment [31,32]. Previous studies have mainly focused on measuring water quality indicators at sampling sites to reflect the regional water environment, but this does not represent the pollution in its entirety. This study used the SWAT model to simulate the NPS pollution of the entire basin of the Tang-Pu Reservoir to explore the influence of land use/cover and landscape pattern on water quality. As a major national reservoir in southeast China, the Tang-Pu Reservoir has played a critical role in supplying water to the local community for drinking and agricultural purposes [33]. The maintenance of good water quality in the basin of the Tang-Pu Reservoir is of great significance to the stability of the local society and the lives of residents. Thus, understanding the relationship between land use/cover patterns and water quality can inform water pollution control and land use/cover management and help the government to make accurate decisions in the protection of water resources in the Tang-Pu Reservoir basin and similar basins. The main objectives of this study were to: (1) monitor the land use/cover and landscape pattern, (2) characterize the water quality and simulate the NPS pollution of water quality in the Tang-Pu Reservoir basin, and (3) reflect the relationships between TN&TP and land use/cover types and landscape pattern.

Study Region
The Tang-Pu Reservoir is located in the Shaoxing City, China (29 • 30 -29 • 57 N, 120 • 24 -120 • 55 E). It belongs to the Xiaoshun River Basin, a tributary of the Cao'e River, and encompasses an area of approximately 550 km 2 with a total length of 70 km. The total storage capacity of the reservoir is 235 million m 3 , and the average water supply is 278.2 million m 3 . The Tang-Pu Reservoir supplies water to the local community of 5 million people for drinking and agricultural purposes. The topography in the basin is characterized by low mountain ranges and hills, and the terrain rises from northeast to southwest ( Figure 1). The basin has a typical subtropical monsoonal climate, with an annual cycle of dry and wet seasons. The mean annual air temperature is 16.5 • C, and the mean annual amount of sunshine is 1903 h. The mean annual precipitation is 1564.4 mm, occurring during March, June, and September. In recent years, the average runoff of the basin is approximately 366 million m 3 , most of which is from precipitation. The soil type in this area is dominated by red soil and Acrisol. The predominant land use/cover type is forest, and the natural vegetation in the area is dominated by tree layers, most of which are masson pine. The main crops cultivated in the basin are cash crops, such as Torreya and tea, which are the main source of economy for the inhabitants of the area. The water environment of the river basin has been polluted by nitrogen and phosphorus recently, which has caused slight eutrophication of the reservoir environment, and the frequency of occurrence of eutrophication has increased [33].

Data
The data used in the study are shown in Table 1 The Shaoxing Environmental Protection Section has established two sampling sites in the Tang-Pu Reservoir basin. One monitoring site, named Spillway, is located in the upstream with little anthropogenic activity, while the other monitoring site, named Shuang-Jiang Stream, is near the reservoir estuary ( Figure 1). Water samples were collected once a month during 2017-2019. Six water quality parameters, including water temperature (TEM), pH, dissolved oxygen (DO), electrical conductivity (EC), TN, TP, and permanganate index (COD Mn), were analyzed according to the environmental quality standards for surface water (Chinese State Environmental Protection Bureau 2002). TEM ( • C), and pH and EC (µs·cm −1 ) were measured on site using a portable conductivity meter (METTLER-SG3, Mettler Toledo, Columbus, OH, USA). DO (mg·L −1 ) was measured on site using a portable oxygen meter (YSI DO200, YSI Incorporated Company, Yellow Springs, OH, USA). TN (mg·L −1 ) was detected using the digestion-UV spectrophotometric method (UV2800, UNICO Company, South Brunswick, NJ, USA). COD Mn (mg·L −1 ) was determined by the potassium permanganate method. TP (mg·L −1 ) was measured using a spectrophotometer (UNICO2100, UNICO Company, South Brunswick, NJ, USA). The pretreatment and the determination of the parameters followed the national standard methods of examining water and wastewater, which were provided by the Ministry of Ecology and Environment, China.  The SWAT was applied in the Tang-Pu Reservoir to assess the impacts of land use/cover on hydrological components. The input data used in the SWAT model included a digital elevation model (DEM), soil data, digital land use/cover maps, meteorological data, hydrological data, and agricultural management information. The quality of the calculations performed is greatly influenced by the detail of the input data. We used a 12.5-m resolution DEM, from which the delineation of the sub-basins was determined as well as some topographic attributes (area, slope, slope length) and characteristics of their channel network [34]. The soil map was constructed from the Harmonized World Soil Database acquired from the Institute of Soil Science, Chinese Academy of Sciences. The soil input file defines the physical properties used by SWAT for all layers in the soil. There were many variables in the soil input file, such as soil type, depth, saturated hydraulic conductivity, and content of clay, silt, and sand. The local meteorological data, including daily precipitation, minimum and maximum temperature, average wind speed, average relative humidity, and solar radiation in the period from 2015 to 2019, were derived from Shaoxing meteorological stations. The hydrological data included the monthly flow data of the Shuang-Jiang Stream monitoring station from 2015 to 2019, and the water quality data included the total nitrogen and total phosphorus (TN&TP) data of the Shuang-Jiang Stream monitoring station from 2017 to 2019. These data were derived from the Shaoxing Ecological Environment Bureau. Agricultural management information was obtained from the local statistics yearbook.

Land Use/Cover Classification
Random forest classification has become an effective method in the remote sensing classification of machine learning [35]. It is an ensemble classifier that based on trees, and it can create new training sets by bagging techniques. When running the random forest method, random feature subspace enables faster tree building, and out-of-bag estimates evaluate the relative importance of each input feature [36].
Referring to the classification system of the Resource and Environment Science and Data Center, the land use/cover classes were categorized into five groups: (1) cultivated land, including paddy field and irrigated land; (2) garden land, including the orchard, tea garden, and other garden; (3) water, including rivers, reservoirs, and ponds; (4) forest, including wooded areas and mixed forest areas; and (5) constructed land, including residential land, road, and industrial land ( Figure 1). The five land use/cover types were classified based on GEE using the random forest algorithm. About 105 samples were selected for each category, of which 80% of the samples were used for random classification, and 20% were used for accuracy validation.

Landscape Metrics for Landscape Pattern
Based on the classification map of land use/cover, landscape metrics were analyzed to characterize land use/cover patterns and investigate the effects of landscape configuration on water quality with the computer program FRAGSTATS 3.3. Because of the correlation among landscape metrics, a total of nine metrics were extracted to reduce the dependency and redundancy of the indices (Table 2), including Patch Density (PD), Largest Patch Index (LPI), Perimeter-Area Fractal Dimension (PAFRAC), Landscape Shape Index (LSI), Landscape Division Index (DIVISION), Splitting Index (SPLIT), Contagion Index (CONTAG), Shannon's Diversity Index (SHDI), and Shannon's Evenness Index (SHEI). Table 2. Descriptions of land use/cover types and landscape pattern metrics.

Metric Definition Description
PD Patch Density is the number of corresponding patches divided by total landscape area. Reflects the fragmentation of landscape.

LPI
The area of the largest patch of the corresponding patch type divided by total landscape area. A measure of dominance. PAFRAC 2 divided by the slope of regression line obtained by regressing the logarithm of patch area against the logarithm of patch perimeter.
Reflects shape complexity across a range of patch sizes.

LSI
The sum of the entire landscape boundary and all edges within the landscape boundary divided by the total landscape area.
Reflects the complex shape of the patches that make up the landscape.

DIVISION
1 minus the sum of patch area divided by total landscape area, quantity squared, summed across all patches of the corresponding patch type.
Reflects the contagion/interspersion of the landscape SPLIT SPLIT equals the total landscape area squared divided by the sum of patch area squared, summed across all patches of the corresponding patch type.
SPLIT increases as the focal patch type is increasingly reduced in area and subdivided into smaller patches Table 2. Cont.

Metric Definition Description
CONTAG Extent to which patch types are aggregated or clumped as a percentage of the maximum possible.
Reflects the degree of agglomeration or extension trend of different patch types.

SHDI
The number of different patch types and the proportional area distribution among patch types.
Compares and analyzes landscape diversity and heterogeneity.
SHEI SHEI equals minus the sum, across all patch types, of the proportional abundance of each patch type multiplied by that proportion, divided by the logarithm of the number of patch types.
Analyzes the diversity of the landscape and reflects a landscape dominated by one or a few dominant patch types.

SWAT
SWAT was applied in the Tang-Pu Reservoir to assess the impacts of land use/cover on hydrological components [23]. In the SWAT model, a basin is divided into multiple sub-basins based on a DEM. The sub-basins are further divided into a series of uniform hydrological response units (HRUs) based on soil and land use/cover. After simulating the hydrological components, sediment yield, and nutrient cycles in each HRU, the results are aggregated for the sub-basins [37]. After adding all kinds of data, i.e., the soil chemical attribute data, local meteorological data, hydrological data, and agricultural management information, to the model, the SWAT model can be operated normally [38]. There are many parameters that have a great effect on the simulation results. To improve the reliability and applicability of the model, it is necessary to continuously modify the parameters of the model to verify the results [29,31]. The SWAT-CUP software provides a variety of uncertainty analysis algorithms. Among them, Sequential Uncertainty Fitting (SUFI2) was selected for sensitivity analysis, calibration, and verification of the model because of its high efficiency and accuracy according to previous studies [39]. In this paper, we regarded the sub-basin where the monitoring site of Shuang-Jiang Stream located as the area for calibration and verification. We selected January 2017 to December 2018 as the nutrient calibration period and January 2019 to December 2019 as the nutrient verification period to check their performances.
During the calibration and validation, the correlation coefficient (R 2 ) and the Nash-Sutcliffe efficiency coefficient (NS) were selected to evaluate the applicability of the model and quantify the degree of fit between the measured value and simulation value. The validated model was applicable as the nutrient indices met the criteria recommended by Moriasi et al. [40], i.e., R 2 > 0.0 and NS > 0.35 for phosphorus; R 2 > 0.30 and NS > 0.35 for nitrogen. The mathematical formulas for the R 2 and NS calculation are as follows: where O i is the ith observed data, O is the mean value of observed data, P i is the ith simulation data, P is the mean value of simulation data, and the n is the total number of the observations.

Statistical analysis
RDA is an ordination method of the regression analysis combined with principal component analysis, which can reflect correlation between response variables and explanatory variables [41]. The explanatory variable that most influential for response can be found among the known environmental factors. RDA concentrates environmental factors in several ranking axes to indicate their impact on response variables, and the ranking axis is linear combination of the environmental factors [42]. RDA is widely used to determine the relationship between environmental factors and land use/cover in recent years. [21,41]. RDA allows us to examine the impact of multiple land use/cover variables on TN&TP simultaneously. There are two important results that can be output from RDA. One of which is ratio (%) of the total variance of TN&TP explained by ranking axes, and the other is the ordination diagrams, which definitely reveal the correlation between land use/cover pattern (environmental factor) and water quality parameters (species). In the species-environment ranking diagram, the direction and length of the arrow represent the correlation between species and environmental factors. The angle between the environmental factor and the species arrow indicate the correlation between the environmental factor and the species. If the angle between the two arrows is less than 90 • , it means that there is a positive correlation between them; if it is greater than 90 • , it indicates a negative correlation; and if it close to 90 • , the correlation is weak. The length of the environmental factor arrow indicates the degree of its influence on the species data. The longer the arrow, the greater the explanation of the environmental factor to the response variable [43]. The RDAs were performed using the CANOCO 5.0 program (Microcomputer Power Company, Ithaca, New York, NY, USA).

Land Use/Cover and Landscape Pattern Change Detection
The results of confusion matrix in 2009 and 2019 are shown in Tables 3 and 4. They show that the overall accuracy of the land use/cover classification in the study area in 2009 and 2019 were 89% and 91%, respectively, and the Kappa coefficients land use/cover were 0.86 and 0.88, respectively, indicating that the classification accuracy met the requirements of land use/cover and cover classification. The land use/cover classification map obtained by the random forest algorithm under the GEE platform is shown in Figure 2. Forest was the dominant land cover in the basin, and its percentage coverage accounted for more than 60%. Forest displayed a decreasing tendency during the years from 2009 to 2019 and was mainly converted to cultivated land and garden land. Garden land was also an important land use/cover in the basin, accounting for a percentage land cover of 16%. There was a slight increasing trend for the area of garden land from the period 2009 to 2019. Cultivated land was widely distributed across the basin, occupying > 11%. Cultivated land in the basin decreased between 2009 and 2019. Constructed land was mainly dispersed around the constructed land. Constructed land exhibited the most obvious change in the area, increasing by 3.53 km 2 in the 10-year study period, and it occupied the area of surrounding cultivated land and gradually expanded outward. These changes emphasize that the Tang-Pu Reservoir is a mountainous area that is mainly covered by forest and is thus less affected by human disturbance. The increase in constructed land was mainly reflected in the increase in rural residential land and road land.

Characteristics of Water Quality
The descriptive statistics of the original data for the seven water quality parameters are shown in Table 5. Each parameter can be identified the classification according to the environmental quality standards for surface water in China (GB3838-2002) [44]. Among them, class III is the minimum standard for human living drinking water, class IV is applicable to industrial water areas, and class V is applicable to agricultural water areas. DO varied greatly from 6.52 to 12.18 mg/L, with 78% of the samples exceeding water quality standards class I. The average concentration of COD Mn was 1.90 mg/L. More than half of the samples (65%) exceeded the class I water quality standard, with 32% exceeding class II water quality standards. The mean value of TP was 0.04 mg/L, with 7, 36, 32, and 22% of the samples exceeding water the quality standards class I, class II, class III, and class IV, respectively. In other words, 25% of the sampling sites for TP value did not meet class III, which is drinking water standard. Nitrogen pollution was the most serious pollution problem in the Tang-Pu Reservoir. The ratios for the mean values of TN compared to the class V standard of 2.0 mg N/L were 1.19. Nearly all the TN concentrations did not reach class III for the drinking water standard. Only about 35% of TN concentrations exceeded the class V water quality standard. The highest concentrations of TN were 2.88 times higher than the class V water quality standard. In summary, all of the water quality indicators in the Tang-Pu Reservoir basin had varying degrees of excessive pollution, especially TN&TP. Nitrogen and phosphorus pollution were the most serious pollution, causing eutrophication of the reservoir. Note: N represents the number of samples, SD represents the standard deviation, and CV represents the coefficient of variation.
From the above study, it is clear that TN&TP were the most influencing factors for water quality in the study area, causing eutrophication of the reservoir. The following study focuses on nitrogen and phosphorus pollution. SWAT model simulated the output and their spatial distribution of TN&TP, and the relationship between land use/cover pattern and TN&TP in the basin was analyzed.

Model Calibration and Validation
Five sensitive parameters were identified: ERORGN, NPERCO, SDNCO, PHOSKD, and PPERCO. Each sensitive parameter was manually calibrated, and the optimal values are listed in Table 6. The comparison of measured and simulation values in calibration and verification period is provided in Figure 4. During the calibration period, the R 2 values of TN and TP were 0.81 and 0.72, respectively, and the NS values of TN and TP were 0.75 and 0.65, showing a close agreement between the measured and simulated values. During the verification period, the R 2 values of TN and TP both reached 0.4 and met the criteria (for TN, R 2 value was 0.7, and for TP, R 2 value was 0.75), and the NS values of TN and TP both exceeded 0.35 and met the criteria (for TN, NS value was 0.63, and for TP, NS value was 0.71), which was considered acceptable and applicable for future application in this study.  For the average monthly TN value, the simulated value of 30.73 t is lower than the measured value of 37.21 t, and for the average monthly TP value, the simulated value of 1.64 t is lower than the measured value of 1.47 t. The reason for the underestimation is precipitation. Since these were wet years in the simulation period, strong rainfall runoff washed NPS pollutants into the river, making a small part of the pollution output remained in the soil and seepage pits. The coefficient of pollutants entering the river was relatively large.

Spatial Variation Characteristics of TN&TP Output Simulated by SWAT
Under the simulated results of SWAT model, the annual TN&TP output in the whole study area is 300.87 t and 24.23 t, respectively, and the unit TN&TP output is 6.61 kg/hm 2 and 0.53 kg/hm 2 , respectively.
Using the natural break point method, the output and output per unit area of TN&TP in the river section were divided into five categories ( Figure 5). The TN&TP displayed similar spatial distribution characteristics. In terms of output and output per unit area, TN was always higher than TP. The TN&TP output were the highest at the junction of Southern Stream and Northern Stream, located in sub-basins 2, 3, and 6, and the TN&TP in the entire basin reached the maximum as the rivers merged. The second and third grades were mainly distributed near the river channels of Southern Stream (located in the sub-basins 10, 12, 16, 17, 21, and 22) and Northern Stream (located in the sub-basins 7, 9, and 11), respectively. The fourth grade was mainly located at the source of Southern Stream and Northern Stream, located in sub-basins 12, 18, 22, 23, 24, and 25. The slope of this area was larger, and with the surface runoff formed by rainfall washing into the river, the TN&TP attached to the sediment, causing the TN&TP to increase. The areas with the lowest TN&TP outputs were scattered in the basin. The areas with large TN&TP output per unit area were located in sub-basins 4, 10, 12, 15, and 22, where the cultivated land and construction land were concentrated and the landscape pattern was relatively complex.

Differences in Land Use/Cover on TN&TP Output
The impact of different land use/cover types on the TN&TP output is explored in Table 7. The TN annual output from garden land contributed the most, accounting for 45.98%, followed by forest and cultivated land. Construction land had the lowest annual output, being only 0.94%. The TP annual output from cultivated land contributed the most, accounting for 44.55%, followed by forest garden land and construction land. From the perspective of the TN&TP output per unit area, the land use/cover types were ranked as follows in term of TN: garden land (18.83 kg/km 2 ) > cultivated land (15.49 kg/km 2 ) > forest (2.76 kg/km 2 ) > construction land (1.92 kg/km 2 ) and as follows in terms of TP: cultivated (2.15 kg/km 2 ) > garden land (0.60 kg/km 2 ) > construction land (0.32 kg/km 2 ) > forest (0.28 kg/km 2 ).

Differences in the Influences of Land Use/Cover on TN&TP among the Scales
The rivers in each sub-basin were divided into six spatial buffer scales of 100 m, 300 m, 500 m, 700 m, 1000 m, and the whole catchment. Explore the impact of land use/cover types on TN&TP under different buffer zones using RDA. From the species-environment ranking of RDA (Figure 6), the distribution of species, the relationship between species and environmental factors, and the relationship between different environmental factors could be derived.
Regardless of the buffer zone, cultivated land and construction land were positively correlated with the TN&TP outputs. On the contrary, forest was negatively correlated with the TN&TP outputs. No matter under which buffer zone, the axes of the garden were almost perpendicular to the axes of TN&TP. Therefore, there was no significant relationship between the proportion of garden land and water area and the output of TN&TP. The relationship between the output per unit area of TN&TP and land use/cover varied with the spatial scale. As the buffer size increased, the proportion of cultivated land area decreased. The angle between the axis of cultivated land and TN became larger, while that between cultivated land and TP become smaller from 100-m buffer to the whole catchment. That is to say, the correlation between the cultivated land and the TN output per unit area weakened, while the correlation with the TP output per unit area increased. For construction land, the proportion of its area gradually decreased as the buffer zone increased. The correlation between construction land and TN&TP per unit area reached the maximum in the buffer range of 100 m buffer scale, and the correlation was weaker at other buffer scales. The proportion of forest area increased with the increase of the buffer zone. The forest had an increased negative correlation with the TN&TP output per unit area at the buffer scale from 100 m to 500 m. The correlation between forest and TN&TP summited at the buffer scales of 500-700 m, after which the negative correlation began to weaken as it became larger than the 700-m buffer scale.

The Impact of Landscape Pattern on TN&TP
The water quality change of the river basin was closely related to the composition of the landscape pattern. In this study, the landscape pattern index and the TN&TP inflow into the river in the sub-basin were used as environmental factors and species variables to explore the impact of the landscape pattern index on the output of TN&TP. Note from Table  8 that the characteristic value of the A1 principal axis of the RDA ranking was 0.731, and the interpretation amount was 73.1%. The correlation coefficients of the A1 principal axis and the A2 principal axis were 0.857 and 0.839, respectively, indicating that the landscape pattern of the Tang-Pu Reservoir basin was related to the output of TN&TP. To further verify the relevance of each landscape pattern factor to the TN&TP output, a t-value double sequence diagram was made for the results of the RDA (Figure 7). The results showed that SHDI and CONTAG had a greater impact on TN&TP. IJI, SPLIT, SHDI, and DIVISION had a strong positive correlation with the output of TN&TP. The larger the SPLIT, IJI, DIVISION, and SHDI, the higher the fragmentation and separation of the landscape and the richer the basin landscape, which may increase the output of TN&TP. There was a certain negative correlation between LPI, CONTAG, and AI and the output of TN&TP. Larger LPI, CONTAG, and AI indicated better landscape integrity and connectivity among landscape patches. It showed that there was less interference from human activities, resulting in a lower degree of landscape fragmentation, so the TN&TP output decreased. In addition, the correlation between the other indexes and the TN&TP output was not very strong. (b) correlation between landscape metrics and TP; if the arrow of a certain environment completely falls on a blue circle, it indicates that it is significantly positively correlated with the species. Oppositely, if the arrow of a certain environmental factor falls within the red circle, it indicates that the species and the environmental factor are significantly negatively correlated. The more the arrow falls within a circle, the more significant the correlation between the species and the environmental factor.

Relationship between Land Use/Cover and TN&TP
Land use/cover types and their spatial distribution affect the water pollution status of water bodies at the basin [4,10,16]. Owing to differences in land use/cover patterns and human activity intensity, they can be significantly correlated, i.e., different land use/cover types have different pollutant outputs [45]. During rainfall, surface pollutants flow into the water with surface runoff, and pollutants may be deposited, absorbed, and reprecipitated when flowing through different land use/cover types [46].
This study showed that construction land and cultivated land had negative effects on water quality, whereas forest had positive effects on water quality. These results are consistent with most previous studies [6,47,48]. Regardless of the buffer zone, the construction land was positively correlated with the TN&TP outputs. The pollution in construction land is transferred to streams via runoff from impervious surfaces, which adversely affects stream water quality after rainfall runoff is washed away [49,50]. The limited construction land area in the Tang-Pu Reservoir basin is the main factor contributing to the poor TN&TP output. Cultivated land strongly influences the nutrient concentrations of stream water [51]. Owing to fertilization, pesticides, and agricultural management measures such as irrigation, the TN&TP output in cultivated land is relatively high, which causes eutrophication of the water during and following precipitation events [52,53].
The large area coverage of forest did not produce a large TN&TP output. Vegetation plays a positive role in soil and water conservation so that the N and P elements in the soil are lost less with surface runoff. That is to say, forest is considered as a sink for potential pollutants in water bodies [52,54]. This is also in line with some research results that the water quality of sub-basins dominated by forest is better than those dominated by cultivated land and garden land [16,55].

Scale Effect on the Influences of Land Use/Cover on TN&TP
The influence of land use/cover pattern on water quality is scale dependent [53,56,57]. Herein, different land use/cover patterns had various scale effects. For example, construction land had greater impacts on TN at smaller scales. This can be explained from two aspects. From the perspective of the nature of N, it may be due to the denitrification of part of the N under the action of microorganisms during the long-distance runoff transport of the N fertilizer applied to the cultivated land [58] or that the positively charged NH + 4 is easily absorbed by soil particles and used by plants. From the perspective of land use/cover structure, as the proportion of forest increases, forest intercepts some of the TN generated by agricultural activities. As a result, the correlation between TN and cultivated land becomes worse as the buffer size increases [54,55]. Construction land was more strongly correlated with TN&TP at smaller scales. This may be related to the distribution of rural residential land along the river and the artificial disturbances, which greatly damage river ecosystems. Forest was highly correlated with improved water quality at the buffer scale of 500-700 m. The wider the riparian zone, the greater the effect of intercepting pollutants [59]. Riparian vegetation plays an important role in maintaining and improving water quality [60]. These findings showed that the area of cultivated land and construction land should be reduced in the buffer zone closer to the river, while the application of fertilizers in cultivated land should be reduced. And the area of forest land should be increased in the buffer zone of 500-700 m to intercept N and P elements.
The complexity of the scale effects indicates how challenging it is to select a single optimal scale to control a large number of water quality parameters because of the different biogeochemical processes occurring at different scales [61]. Therefore, a multi-scale perspective should be provided for sustainable utilization of water resources and land use/cover planning. In exploring the relevance of land use/cover and water quality, the result is not only affected by the proportion of land use/cover area but also by the spatial pattern of land use/cover to a certain extent; that is to say, the patch fragmentation and spatial interconnection of land use/covers played an important role in modulating land use/cover effects on water quality in this study, as supported by several studies [62,63]. Sun et al. [64] found that an aggregated forest landscape had a greater ability to fix and absorb pollutants in streams than a dispersed forest landscape. This suggested that it might be possible to engineer a land use/cover configuration to improve water quality without appreciably altering the land use/cover composition in the study area.

Relationship between Landscape Pattern and TN&TP
The pattern and process theory of landscape ecology explains the spatial differences between water quality. One of the reasons is that water quality is a comprehensive reflection of multiple land use/cover composition and land use/cover spatial distribution [65]. Thus, when there is heterogeneity in different sub-basins, the degree of impact on the water quality of river also differs. Our landscape pattern conclusions were consistent with some previous studies, but there were also inconsistencies. For example, similar to the results of Lee et al. [66] and Shi et al. [53], water quality was positively correlated with SHDI and negatively correlated with CONTAG. The reason is that the height of SHDI represents the richness of diversity and heterogeneity of the landscape, which is related to the increase in the types of landscape patches, leading to a balanced distribution of various patch types [67]. It may also be related to the richness of patches and the complexity of the landscape structure [68]. CONTAG reflects the degree of separation and spread of the landscape. A high spread value indicated that a certain dominant patch type in the landscape has formed good connectivity, and the degree of fragmentation of the landscape is low. Obviously, for this study area where forest is the dominant landscape type, a high SHDI and low LPI indicate intense human interference, and the dominance of pollution source landscape types (i.e., construction land and cultivated land) will increase. The degree of fragmentation of forest will increase simultaneously, resulting in the increasing risk of water quality degradation. Contrary to our result, Sun et al. [69] and Zhang et al. [70] found that with the increase of SHDI, river water quality tended to be better, and at the same time, CONTAG was positively correlated with river water quality. The reason for this result may be related to the dominant species in the study area. The dominant landscape type of the study areas of Sun et al. [69] and Zhang et al. [70] were pollution source landscapes, that is, basins dominated by construction land or cultivated land. This is contrary to our study area, and therefore the higher the CONTAG, the more obvious the plaque aggregation and connectivity. It led to an increase in the concentration of pollutants in the water. The same applies for SHDI. In addition, the reason for the inconsistency of some results may be related to the research scale. Li et al. [71] studied the correlation between the landscape pattern and water quality from a series of spatial scales, and the results showed that it is at the sub-basin scale that CONTAG and the concentration index of various pollutants in water have a negative correlation, while SHDI is associated with a positive correlation. At other spatial scales, the result is the opposite.
In general, owing to differences in regional characteristics and spatial and temporal scales, the landscape index and water quality did not reach a unified conclusion and demonstrated many uncertainties, indicating that further landscape-water quality correlation research is needed.

Limitations
Some limitations should be noted in this study. First, due to the lack of long-term data in the basin, only one year of water quality was simulated by SWAT model. Thus, only the effect of land use/cover on the NPS pollution in 2019 was analyzed. Although the results of the relationship between land use/cover and TN&TP in 2019 remain consistent with most studies, it is recommended that more years need to be investigated to improve the confidence of the result. Second, the SWAT model can only simulate the spatial distribution of nutrient contamination in the basin but not the other water quality parameters, i.e., DO, EC, and COD Mn . TN&TP do not fully represent the water quality of the basin. Based on the study results, more spatial distribution of other water quality parameters can be simulated to explore the impact of land use/cover on water quality indicators.

Conclusions
A better understanding of the impact of land use/cover pattern on water quality could provide guidelines for NPS pollution caused by improper land use/cover to achieve sustainable water resources utilization. Although SWAT model can simulate the TN&TP output to display the overall NPS pollution situation in entire basin, its results were rarely used to explore the influence of land use/cover and landscape pattern on water quality. Based on the outcomes of monitoring the land use/cover pattern and assessing the water quality in the Tang-Pu Reservoir basin, our study simulated the output values and their spatial distribution of the most influencing factors for water quality, that is, the TN&TP, and explored the correlations between land use/cover pattern and TN&TP output by the redundancy analysis.
We found that the land use/cover type was little affected between the period 2009 and 2019. Forest was the dominant land cover, displaying a decreasing tendency. Constructed land exhibited the most obvious change, occupying the surrounding cultivated land, and cultivated land was decreasing in the basin. The landscape fragmentation increased, and the connectivity decreased. The assessment characterized the water quality showed that TN&TP was the most influencing factor, and the high N and P pollution in the reservoir contributed a great deal to the eutrophication problems. Under the simulation of SWAT, the TN&TP outputs were highest at the junction of rivers, and the output per unit TN&TP was at the areas where the landscape pattern was highly fragile. We also found that owing to fertilization, pesticides, and agricultural management measures, i.e., irrigation, the output per unit area TN&TP of cultivated land and garden land was higher than that of construction land and forest. There existed a positive correlation between the TN&TP and cultivated land and construction land, which means that the expansion of cultivated land and construction land inevitably increased the NPS pollution in the basin. Vegetation plays a positive role in soil and water conservation so that the N and P elements in the soil are lost less with surface runoff. Meanwhile, different land use/cover types had different scale effects. The correlation between forest and TN&TP summited at 500-700-m buffer and construction land at 100-m buffer. As the buffer size increased, the correlation between the cultivated land and the TN weakened, while the correlation with the TP increased. These findings indicated that land use/cover planning in water quality management needs to adopt a multi-scale perspective. The TN&TP outputs in the sub-basin were positively correlated with SHDI and negatively correlated with CONTAG in the area where forest is the dominant land cover. Therefore, the area of cultivated land and construction land should be reduced in the buffer zone closer to the river, and the application of fertilizers in cultivated land should be reduced. The more connected forest land should be increased in the 500-700-m buffer zone to intercept N and P elements. These will be necessary to promote sustainable water resources utilization.  Data Availability Statement: All data generated or analyzed during this study are included in this published article.

Conflicts of Interest:
The authors declare no competing interest.