The Characteristics of Net Anthropogenic Nitrogen and Phosphorus Inputs (NANI/NAPI) and TN/TP Export Fluxes in the Guangdong Section of the Pearl River (Zhujiang) Basin

: Human activities have greatly inﬂuenced the inputs and cycling pathways of nitrogen (N) and phosphorus (P), causing dramatic environmental problems in the Pearl River Basin. In this study, the characteristics of net anthropogenic nitrogen and phosphorus inputs (NANI/NAPI) were analyzed in the Guangdong section of the Pearl River Basin from 2016 to 2020. NANI showed a very slight decrease trend from (1.51 ± 0.09) × 10 4 to (1.36 ± 0.08) × 10 4 kg · N · km − 2 · yr − 1 , while the average intensity of NAPI was 3.8 × 10 3 kg · P · km − 2 · yr − 1 . Both NANI and NAPI intensities were at high levels, resulting in the serious deterioration of water quality in the Pearl River Basin. Fertilizer input was the most important component for the intensities of NANI and NAPI, accounting for 38–42% and 53–56%. However, in the Pearl River Delta, the major components of NANI and NAPI were the human and animal consumption (food/feed) inputs and non-food net phosphorus input. The input of NANI and NAPI should be controlled for different areas, based on the differing driving forces, to alleviate the deterioration of water quality. This study of NANI and NAPI in the Pearl River Basin is one of the important prerequisites for clarifying the input and water quality, providing support for further effective control of nitrogen and phosphorus pollution in the Pearl River.


Introduction
Biochemical flows (nitrogen and phosphorus), together with climate change, landsystem change, and biodiversity loss, have become the most noteworthy ecological problems in the world, exceeding the "Planetary Boundaries" [1].These extensive biochemical flows lead to severe environmental problems in watersheds or oceans, such as eutrophication, algal blooms, and ocean acidification [2].However, with economic development, population growth, and food demand, the use of chemicals containing nitrogen (N) and phosphorus (P) are rapidly increasing, leading to excessive N and P input in the watershed system [3].Therefore, it is essential to identify the main sources of anthropogenic N and P inputs and take abatement actions to prevent the deterioration of water quality.
Reasonable and effective control measures need to make an inventory of the input intensity of each source.To accurately estimate the nutrient inputs and pathways, various models were used for cycling pathways, such as SWAT [4], MARINA [5], and HSPF [6].In addition, the establishment of a pollutant input inventory, based on the mass balance method to calculate the input rates, was applied to large watersheds [7].To assess the amount of nitrogen input caused by human activities in the watershed, the concept of net anthropogenic nitrogen input (NANI) was proposed by Howarth et al. in 1996 [8].Compared with other models, the NANI model can more accurately predict nitrogen flux as a linear or index function of NANI, especially in less agricultural and more forested watersheds [9][10][11][12].Based on the method of the NANI estimation model, net anthropogenic phosphorus input (NAPI) was expanded by Russell et al., adding the net flow input module of human activities in the basin, including the amount of fertilizer consumption, net phosphorus input of food/feed, and non-food phosphorus (mainly from detergents used by humans) [13].In addition, the spatiotemporal distributions of NANI and NAPI have been calculated by studies, providing precise information for pollution control management for policymakers [14][15][16][17].These NANI and NAPI intensity estimations, applied to all river basins, can assess the nitrogen and phosphorus inputs caused by human activities in the region and analyze their input sources and spatial and temporal distributions; this can effectively identify the regions and influencing factors with more nitrogen and phosphorus accumulation, and put forward management suggestions according to local conditions.
Guangdong province has the most developed economy and the highest population density in the Pearl River Basin of China, with high-density industrial and agricultural activities.A large number of industrial, agricultural, and other human activities has led to excessive N and P inputs, which are mainly responsible for environmental risk in the Pearl River Basin.Especially in the Pearl River Delta, nitrogen and phosphorus inputs are extremely high [7].However, there are obvious differences between the Pearl River Delta and other regions, because agricultural activities are frequent in non-Pearl River Delta regions.The agricultural and industrial activities in the rapidly urbanizing basin have dramatically altered the nitrogen and phosphorus cycles compared to other areas.In this study, based on the NANI/NAPI estimation approach, we evaluated the NANI and NAPI intensities in the Guangdong section of the Pearl River Basin.Specifically, this study (a) calculates the NANI and NAPI in the Guangdong section of the Pearl River Basin from 2016 to 2020, (b) analyzes the temporal and spatial intensity of NANI and NAPI, especially between the Pearl River Delta and other regions, (c) analyzes the variation and driving force of NANI and NAPI, and (d) to explores the relationship between TN and TP flux, and NANI and NAPI.

Study Area
The Pearl River, the second largest river and the third longest river in China with an annual runoff of more than 300 billion m 3 , is the most economically developed and densely populated region of China, and contributed 9.2% of the gross domestic product in 2016.In particular, the Guangdong section of the Pearl River Basin (21 • 13 -25 • 31 N, 109 • 39 -117 • 19 E) has the most developed economy and the highest population density.The downstream of the Pearl River overflows into a river network in Guangdong Province, converging in the Pearl River Delta network river area and finally flowing into the South China Sea.However, due to high-density industrial and agricultural activities, a large amount of N and P pollution has input in the Pearl River Basin, leading to the deterioration of water quality [18][19][20].The poor-quality water has serious impacts on the acidification or eutrophication of the ecological environment of the basin, estuary, and oceans [21][22][23].Therefore, Guangdong province was selected to calculate NAPI and NAPI intensities.The statistics and analysis of nitrogen and phosphorus input intensities were conducted in Guangdong Province and its 21 major cities, and the monitoring sections were selected in the drainage system of the Pearl River Basin, including West River, East River, North River, and various rivers on the Pearl River Delta (Figure 1).
sities.The statistics and analysis of nitrogen and phosphorus input intensities w ducted in Guangdong Province and its 21 major cities, and the monitoring sectio selected in the drainage system of the Pearl River Basin, including West River, Ea North River, and various rivers on the Pearl River Delta (Figure 1).

Data Collection
To estimate NANI and NAPI intensities, the basic data used to calculate i included three categories: social data (population, land-use type), agricultural a data (fertilizer application amount, agricultural product planting area and yie poultry stock), economic data (Gross Domestic Product-GDP, Gross Output V Agriculture-GOVA, and Gross Output Value of Industrial-GOVI.To evaluate the quality and flux of nitrogen and phosphorus in the Pearl Rive the water quality data for 2016-2020 from 96 monitoring sections of the Pearl Riv in Guangdong were collected from the automatic water quality monitoring plat the China General Environmental Monitoring Station (www.cnemc.cn(accesse July 2022)) (Figure 1).The data were the monthly average data of monitoring sectio data were filtered to eliminate the extreme values and then the missing points we  To evaluate the quality and flux of nitrogen and phosphorus in the Pearl River Basin, the water quality data for 2016-2020 from 96 monitoring sections of the Pearl River Basin in Guangdong were collected from the automatic water quality monitoring platform of the China General Environmental Monitoring Station (www.cnemc.cn(accessed on 12 July 2022)) (Figure 1).The data were the monthly average data of monitoring sections.The data were filtered to eliminate the extreme values and then the missing points were completed by using the spatial interpolation of inverse distance weighting (IDW) and the Kriging method in GIS (ArcGIS 10.5) [24,25].

The Flux Calculation
The fluxes of nitrogen and phosphorus are obtained via the product of average period concentration and average period flow, showing the period averaging method.
where W is the flux of nutrient elements in the Pearl River; Q(t) is the function of the period, and C(t) is the function of the period; Qi is the instantaneous runoff flux; Ci is the instantaneous concentration of the nutrient element; K is adjustable constant; Qr is the monthly average runoff flux.

Calculation of NANI and NAPI
Net anthropogenic nitrogen and phosphorus inputs are accounting methods for controlling the input intensities of nitrogen and phosphorus sources in watershed human activities, generally including 6 parts: (a) atmospheric deposition, (b) fertilizer input, (c) seeding input, (d) agricultural fixation, (e) human and animal consumption (food/feed) inputs, and (f) other input (the non-food inputs of N and P by rural residents), which represents the input of exogenous nitrogen and phosphorus in the basin caused by human activities.The calculations are as follows: NAPI = P dep + P fer + P seed + P non + P im (4) where N dep and P dep are atmospheric depositions; N fer and P fer are chemical fertilizer inputs; N seed and P seed are the seeds input intensity of each crop; N im and P im are the net food and feed inputs of N and P; N fix is the agricultural fixation by crops, while P non is the non-food inputs of P by rural residents (such as the P-based detergents used for household cleaning).All the units are kg•km −2 •yr −1 .
(1) Calculation of atmospheric deposition To accurately estimate atmospheric nitrogen deposition in the Pearl River Basin, the literature on nitrogen dry and wet deposition from 2016 to 2020 was comprehensively gathered.The data were from scientific databases, including the CNKI database, WanFang database, ISI Web of Science, Google Scholar, and Science Direct databases, listed in Table 1.In the Pearl River and nationwide Pearl River Basin, the atmospheric nitrogen deposition presents a stable trend, with an approximate decrease from 2016 to 2020 [26].The average atmospheric total nitrogen deposition was 3567 kg•km −2 •yr −1 , mainly including NO x , NH y , and organic nitrogen (ON).The rate of NH y and NO x were approximately 1.1-1.3[27,28].NO x and ON were derived from fuel combustion sources, dust, and vehicle input [29][30][31].However, NH y was derived from agricultural activities, such as animal husbandry and agricultural fertilization, accounting for more than 80% [32][33][34][35], and just a little from vehicle inputs [36,37].Therefore, to avoid duplicate calculations with agricultural sources, 80% of NH y deposition was subtracted from the atmospheric nitrogen deposition calculation.Based on the comprehensive literature, the atmospheric nitrogen deposition rate (N dep ) in the Pearl River Basin was determined as 2.4 For atmospheric phosphorus deposition, it can be almost limited.Based on Pan et al., the global P deposition was ~100 kg•P•km −2 •yr −1 covering the years 1954 to 2020, and P atmospheric deposition in Asia was slightly higher than that in Europe [38].Therefore, according to Ma et al. and Zhang et al., the phosphorus deposition rates (P dep ) in the Pearl River Basin was 92 kg•P•km −2 •yr −1 from 2016 to 2018 and 108 kg•P•km −2 •yr −1 from 2019 to 2020 [39,40].(2) Calculation of fertilizer input The fertilizers required for agricultural activities included chemical fertilizers and organic fertilizers.Since organic fertilizer mainly comes from the interior of the region and is not an external nitrogen source input, only chemical fertilizer was considered in fertilizer input.The input amount of nitrogen fertilizer (N fer ) was derived from chemical N (ammonium nitrate, ammonium bicarbonate, urea-N, etc.) and combined fertilizers with N content of 12.8% [12].Similarly, the input amount of phosphorus fertilizer (P fer ) was derived from chemical P and combined fertilizers with P content of 15.0% [14].
(3) Calculation of seeds input intensity of each crop The content of seed was also an input for agricultural budgets.Wheat, rice, yam, soybeans, vegetables, and peanuts were considered to estimate seed nitrogen and phosphorus inputs, which were the major agricultural plants in the Pearl River Basin [7].The estimation of N seed and P seed was based on the area of cultivated land in the Pearl River Basin, multiplied by the input of nitrogen and phosphorus per unit area of the seed of each crop, with the different correlation coefficients shown in Table 2. (4) Calculation of agricultural fixation Some crops have the quality of nitrogen fixation, which is an important source of nitrogen input in the Pearl River Basin.On the contrary, there is no phosphorus fixation.Therefore, N fix was calculated via multiplying each N-fixing crop area by its N fixation rate, as shown in Table 3. (5) Calculation of net food and feed inputs of N and P Net food and feed nitrogen and phosphorus inputs refer to the quality balance between the N and P outputs of livestock and crops and the N and P consumptions of humans and livestock, calculated by subtracting the content of nitrogen and phosphorus in livestock products and crop products from the consumption of nitrogen and phosphorus in human food and livestock feed: where HC and LC were the consumption of nitrogen and phosphorus by humans and livestock, respectively.CP was the nitrogen output of crop products, while CPL was the nitrogen loss of crop products due to corruption; similarly, LP was the nitrogen output of livestock products, and LPL was the nitrogen loss due to corruption and inedibility of livestock products.The different correlation coefficients were shown in Table 4, and the loss due to corruption and inedibility was calculated as 10% of the total.When the production of food and feed in the Pearl River Basin was not enough to meet the needs of humans and animals, the result is positive and needs to be input, i.e., net food and feed inputs of N im and P im , and vice versa.(6) Calculation of the non-food input Due to the economic development and population density increase in the Pearl River Basin, non-food net phosphorus input (P non ) was also one of the important sources of NAPI.P non mainly comes from detergents used in people's daily life.P non was calculated by multiplying the total number of study units by the per capita non-food phosphorus input coefficient (0.63 kg

Uncertainty Analysis
To gain insight into the uncertainty of NANI and NAPI estimates, an uncertainty analysis was performed using the Monte Carlo simulation.Due to the regional characteristics of N and content input intensity of nitrogen and phosphorus sources, it was assumed that all the coefficients of NANI and NAPI calculation models approximately follow a normal distribution with a coefficient of variation of 30% [18].A total of 10,000 Monte Carlo simulations were performed to obtain the mean and 95% confidence interval.

Grey Relational Analysis
Grey relational analysis is one of the grey system theories of the geometric proximity between different discrete sequences, to solve certain problems for solving the complicated interrelationships among the multiple performance characteristics [48,49].The grey relational degree, which is viewed as a measure of the similarities of discrete data that may be ordered in a sequential sequence, is used to characterize closeness, i.e., the higher grey relational degree has a higher relationship [50].Furthermore, grey relational analysis has also been used in amounts of studies evaluating environmental water quality and the impact of socioeconomic factors [7,[50][51][52].To analyze the impact of socioeconomic factors, grey relational grade analysis was conducted between NANI and NAPI intensities and the socioeconomic factors, such as Gross Domestic Product (GDP), Gross Agricultural Output Value (GAOV), N Fertilizer (NF), P Fertilizer (PF), Population Density (PD), and Urbanization Rate (UR).

Temporal and Spatial of NANI and NAPI
The temporal and spatial changes of NANI in the Guangdong section of the Pearl River Basin from 2016 to 2020 are shown in Figures 2 and 3. On the spatial scale, it exhibited obvious regional distribution, the areas being the southwest coast (Zhanjiang, Maoming), eastern coast (Shantou, Jieyang), northern region (Shaoguan, Qingyuan, Heyuan), and the Pearl River Delta.The intensity of NANI in the Pearl River Basin was higher in the Pearl River Delta and coastal (southwest and eastern) regions than that in the northern region (Figure 2).The highest area was the Pearl River Delta, with an average NANI intensity of (2.2-2.3)× 10 4 kg•N•km −2 •yr −1 from 2016 to 2020.In particular, Shenzhen, the most developed city with the highest population density, had a NANI intensity of (4.1-4.9)× 10 4 kg•N•km −2 •yr −1 .The intensity of NANI in the northern region was the lowest, with an average intensity NANI of (8.1-9.4) × 10 3 kg•N•km −2 •yr −1 from 2016 to 2020.On the temporal scale, the NANI intensity of Shenzhen continued to increase, while all other cities showed a slight trend of decreasing or maintaining NANI intensities from 2016 to 2020.This decreasing or maintaining trend began in 2012, with a decrease in chemical product consumption [7].The temporal and spatial changes of NAPI were similar to NANI, which also presents higher in the Pearl River Delta and coastal regions than that in the northern region (Figure 3).Unlike NANI, on the temporal scale, NAPI remained stable from 2016 to 2020.The unbalanced population distribution and uneven economic growth across the sub-basins were connected to the temporal and geographical disparities in the distribution of NANI and NAPI intensities across the Guangdong region of the Pearl River Basin [53,54].This spatial pattern of NANI and NAPI was similar to other basins or watersheds in the world [10].
The temporal and spatial changes of nutrient element concentrations had similar patterns compared with NANI and NAPI intensities, which were higher in the Pearl River Delta and coastal regions than that in the northern region (Figures 2 and 3).It indicated that NANI and NAPI had a positively correlated influence on the N and P concentration in the Pearl River.The concentration of TN and TP had different pollution degrees and patterns; these showed that that TN was exceeding the standard across of basin, while the water quality of the Pearl River was fine for phosphorus.High intensities of NANI and NAPI had a very serious impact on the water quality of the Pearl River.The average concentration of TN in the Guangdong section of The Pearl River were 2.65, 2.96, 2.86, 2.40, and 2.31 mg•L −1 from 2016 to 2020, respectively, and the average concentration of TP was 0.15, 0.17, 0.16, 0.15, and 0.08 mg•L −1 from 2016 to 2020, respectively.As shown in Figure 4, poor quality water sections for TN were 52%, 55%, 50%, 39%, and 47% from 2016 to 2020, respectively, based on the Quality Standard and Limit of Surface Water Environment of China (GB3838-2002).As for phosphorus, the poor quality water sections at the source of drinking water were 83%, 86%, 87%, and 94% from 2016 to 2020, respectively, which had been alleviated in 2020.In addition, there was a difference between nutrient element concentrations and NANI and NAPI intensities on the spatial scale, which was evident in the two cities on the west bank of the Pearl River Estuary, i.e., Zhuhai and Zhongshan.Since the population density and agricultural activities of these two cities are not very high, the calculation of NANI and NAPI showed low intensity.The TN and TP in the Pearl River Estuary converged the pollution from upstream, showing a higher TN and TP concentration distribution [18,20].Therefore, the high intensity of NANI led to serious TN pollution in the Pearl River Basin of Guangdong Province, carrying great risks, especially in the Pearl River Delta region.Similarly, the Pearl River Delta region exhibited TP pollution by NAPI.The water quality deterioration had a very adverse impact on the ecological environment and domestic water safety of the Pearl River Basin.The temporal and spatial changes of nutrient element concentrations had similar patterns compared with NANI and NAPI intensities, which were higher in the Pearl River Delta and coastal regions than that in the northern region (Figures 2 and 3).It indicated that NANI and NAPI had a positively correlated influence on the N and P concentration in the Pearl River.The concentration of TN and TP had different pollution degrees and patterns; these showed that that TN was exceeding the standard across of basin, while the   The temporal and spatial changes of nutrient element concentrations had similar patterns compared with NANI and NAPI intensities, which were higher in the Pearl River Delta and coastal regions than that in the northern region (Figures 2 and 3).It indicated that NANI and NAPI had a positively correlated influence on the N and P concentration in the Pearl River.The concentration of TN and TP had different pollution degrees and in the Pearl River Estuary converged the pollution from upstream, showing a higher TN and TP concentration distribution [20,18].Therefore, the high intensity of NANI led to serious TN pollution in the Pearl River Basin of Guangdong Province, carrying great risks, especially in the Pearl River Delta region.Similarly, the Pearl River Delta region exhibited TP pollution by NAPI.The water quality deterioration had a very adverse impact on the ecological environment and domestic water safety of the Pearl River Basin.The human and animal consumption (food/feed) inputs were the second dominating contributor of NANI, accounting for 31.4-33.2%and the non-food net phosphorus input was the second dominating contributor of NAPI, accounting for 27.8-30.2%.Similarly, there were regional differences in the driving force of NANI and NAPI.In the Pearl River Delta, economically developed and densely populated, human food and animal feed are the main components [7].However, the components of NANI and NAPI had strong regional deviations.For instance, the major component of NANI in Shenzhen, the core city of the Pearl River Delta region, was the human and animal consumption (food/feed) inputs, accounting for 83-92%.The major components of NAPI in Shenzhen were the net food and feed inputs, and the non-food inputs, accounting for 35-43% and 45-53%, respectively.Shenzhen requires many food imports due to extremely high population density (8828 person•km −2 in 2020) and almost no agricultural activities.On the contrary, in the northern region of the Pearl River Basin, fertilizer input was the most important component for the intensities of NANI and NAPI, due to intensive agricultural activities.For example, fertilizer input accounted for 43-54% of NANI intensity, and 66-75% of NAPI intensity in Heyuan.

NANI and NAPI and the Driving Force
imports due to extremely high population density (8828 person•km −2 in 2020) and almost no agricultural activities.On the contrary, in the northern region of the Pearl River Basin, fertilizer input was the most important component for the intensities of NANI and NAPI, due to intensive agricultural activities.For example, fertilizer input accounted for 43-54% of NANI intensity, and 66-75% of NAPI intensity in Heyuan.To analyze regional differences, grey relational grade analysis was conducted between NANI and NAPI intensities and the socioeconomic factors [7,52].From the perspective of the whole Guangdong Province, the grey relational grade analysis showed that the primary factor in changing NANI and NAPI intensities was the Gross Agricultural Output Value in the Guangdong section of the Pearl River Basin, with relational grades of 0.914 and 0.923 for NANI and NAPI, respectively (Table 5).Fertilizers N and P also had a high grey correlation with the inputs of NANI and NAPI.Although Guangdong Province has a developed economy, high population density, and high urbanization rate, its NANI and NAPI mainly come from agricultural non-point source emissions.Therefore, agricultural activities were still the main factor for NANI and NAPI intensities in Guangdong Province.However, there were regional differences within Guangdong Province, i.e., the Pearl River Delta and non-Pearl River Delta region.In urbanized areas of the Pearl River Delta region, the industrial and agricultural structure changed rapidly, and intensive fertilizer application was no longer the primary cause of NANI and NAPI.In the Pearl River Delta, the population density and urbanization rate can explain the sources of NANI and NAPI, with higher grey relational grades.However, this tight correlation was lost in the non-Pearl River Delta region and had been replaced by the Fertilizers N and P, directly To analyze regional differences, grey relational grade analysis was conducted between NANI and NAPI intensities and the socioeconomic factors [7,52].From the perspective of the whole Guangdong Province, the grey relational grade analysis showed that the primary factor in changing NANI and NAPI intensities was the Gross Agricultural Output Value in the Guangdong section of the Pearl River Basin, with relational grades of 0.914 and 0.923 for NANI and NAPI, respectively (Table 5).Fertilizers N and P also had a high grey correlation with the inputs of NANI and NAPI.Although Guangdong Province has a developed economy, high population density, and high urbanization rate, its NANI and NAPI mainly come from agricultural non-point source emissions.Therefore, agricultural activities were still the main factor for NANI and NAPI intensities in Guangdong Province.However, there were regional differences within Guangdong Province, i.e., the Pearl River Delta and non-Pearl River Delta region.In urbanized areas of the Pearl River Delta region, the industrial and agricultural structure changed rapidly, and intensive fertilizer application was no longer the primary cause of NANI and NAPI.In the Pearl River Delta, the population density and urbanization rate can explain the sources of NANI and NAPI, with higher grey relational grades.However, this tight correlation was lost in the non-Pearl River Delta region and had been replaced by the Fertilizers N and P, directly reflecting high-density agricultural activities rather than GAOV.The primary input components of NANI and NAPI upstream of Pearl River were mostly caused by agricultural systems as fertilizer (>40% NANI, >60% NAPI).The Pearl River Delta Basin's fertilizer input intensities for N and P were not high, because of a lower amount of agricultural land (only 15%) [7].This discrepancy results from the different rates of urbanization and growth of the population.To meet the demand for nitrogen and phosphorus consumption, due to the high urbanization rate of cities, there might be an increase in agricultural activities in surrounding cities and increase the budget of N and P [55,56].

Riverine Fluxes and Linking Basin NAPI and NAPI Budget
To explore the relationship between the transportation of the nutrient elements and NANI and NAPI in the Guangdong section of the Pearl River Basin, the temporal and spatial changes of the TN and TP fluxes were plotted on a thematic map, as shown in Figures 6 and 7.Both the TN and TP fluxes were significantly larger in the Pearl River Estuary than those in the upstream monitoring sections, and obviously, it was clear that the fluxes of TN and TP increased from upstream to downstream of the Pearl River Basin.Similar to the temporal and spatial changes of the NANI and NAPI intensities, the TN and TP fluxes were higher in the Pearl River Delta and coastal regions than that in the northern region.It might indicate that the fluxes of nutrient elements in the Pearl River Basin were accumulating throughout the basin.The average TN flux of each monitoring section was 6.4 × 10 5 , 7.4 × 10 5 , 7.0 × 10 5 , 6.6 × 10 5 , and 3.2 × 10 5 t•yr −1 from 2016 to 2020, respectively.As for TP, the average TN flux was 3.7 × 10 3 , 4.2 × 10 3 , 4.4 × 10 3 , 3.8 × 10 3 , and 1.3 × 10 3 t•yr −1 , from 2016 to 2020, respectively.The fluxes of TN and TP were decreased compared with previous studies, but they were still at a high level and had adverse effects on the marine ecological environment [17,57,58].The above models accounted for 92% and 78% of the variation in annual TN and TP fluxes.This model's output showed that, in comparison to individual variables, the combined effect of NANI and water discharge may better explain the interannual variations in riverine TN and TP export from the Pearl River Basin.The TN and TP fluxed accounted for 0.49-45.3%and 0.62-24.8% of NANI and NAPI in all cities, respectively, except Zhuhai and Zhongshan.Howarth et al. also found that 20-25% of the total nitrogen input from human activities in the watershed finally enters the water body, and 75-80% of the nitrogen was stored in the terrestrial ecosystem or enters the atmosphere [8].By comparing the relevant research results of NANI in various regions, Swaney et al. proved the relationship between the input of human nitrogen into rivers and the output of nitrogen to coastal waters, indicating that rivers usually export 15-30% of NANI to coastal waters.This indicates that the nitrogen input into the basin by human activities is the main nitrogen source in most of the world's watersheds and coastal waters [11].Therefore, the high intensities of NANI and NAPI might cause N and P to remain in the sediment, which will increase the N and P input load into the ocean.This N and P might accumulate long after inputs and continue to be released, by balance calculation of input and output flux, over the past 30-70 years [60].Strengthening the collaborative management of NANI and NAPI and residual N and P is the key to effectively controlling nitrogen and phosphorus pollution in the Pearl River.

Ecological and Management Implications
The riverine TN and TP fluxes were significantly correlated with NANI and NAPI via nonlinear correlation.Therefore, the input of NANI and NAPI should be controlled based on the driving force, to alleviate the deterioration of water quality and its impact on the oceans.The majority of environmental management studies on regulating N and P in intensive agriculture systems were geared toward reducing pollution in river basins [7].However, it is necessary to consider the NANI and NAPI input of regional heterogeneity.Human food and animal feed were the most important contributors of NANI and NAPI in the Pearl River Delta, while fertilizers N and P were the most important contributors in the non-Pearl River Delta region.Therefore, the investigation of pollutant sources and the implementation of source-control policies is indispensable.According to the contributors of NANI and NAPI of regions in the Guangdong Province, the management and control measures of N and P can be established in agricultural activities and personal consumption habits.In the Pearl River Delta, the reduction in the consumption of food N and P would benefit from actions taken to build a more sustainable social and environmental food system.In the non-Pearl River Delta region, measures such as precision fertilization should be taken to improve the fertilizer utilization rate and reduce fertilizer N and P usage.In addition, intensive water treatment facilities are required to improve discharge efficiency, to improve the water quality of the Pearl River Basin.

Conclusions
The data of average riverine nitrogen and phosphorus fluxes and net anthropogenic inputs of nitrogen and phosphorus from 2016 to 2020 in the Guangdong section of the Pearl River Basin were summarized and analyzed, in order to further effectively control nitrogen and phosphorus problems in the Pearl River.It provided a better understanding of anthropogenic inputs of N and P by the calculation of NANI and NAPI in the Guangdong section of Pearl River, where the intensities of NANI and NAPI were still higher than levels in China.The driving force of the intensity of NANI and NAPI in the Pearl River Basin was different from Pearl River Delta, coastal and northern regions.Agricultural activities were still the primary factor in changing NANI and NAPI intensities, how- The estimation methods of NANI and NAPI intensities were calculated by the coefficient of each pollution source from scientific databases, references, field scale measurement, small watershed monitoring, and model inversion, which existed for regional deviation.To better estimate the NANI and NAPI intensities in different regions, water quality monitoring, and sewage discharge were used for comprehensive evaluation.The riverine TN and TP fluxes were significantly correlated with NANI and NAPI via nonlinear correlation, which was tightly linked to human activities [16,59].A multivariate regression model for predicting the TN and TP output fluxes ([F N , kg•N•km −2 •yr −1 ] and ([F P , kg•N•km −2 •yr −1 ]) with flow Q (10 10 m 3 /yr), and NANI or NAPI in the Guangdong section of the Pearl River Basin, was constructed.According to the simulation accuracy index (R 2 ), the model formed of various combinations of the variables was compared and selected, and the following optimal multiple regression model forms were selected: The above models accounted for 92% and 78% of the variation in annual TN and TP fluxes.This model's output showed that, in comparison to individual variables, the combined effect of NANI and water discharge may better explain the interannual variations in riverine TN and TP export from the Pearl River Basin.The TN and TP fluxed accounted for 0.49-45.3%and 0.62-24.8% of NANI and NAPI in all cities, respectively, except Zhuhai and Zhongshan.Howarth et al. also found that 20-25% of the total nitrogen input from human activities in the watershed finally enters the water body, and 75-80% of the nitrogen was stored in the terrestrial ecosystem or enters the atmosphere [8].By comparing the relevant research results of NANI in various regions, Swaney et al. proved the relationship between the input of human nitrogen into rivers and the output of nitrogen to coastal waters, indicating that rivers usually export 15-30% of NANI to coastal waters.This indicates that the nitrogen input into the basin by human activities is the main nitrogen source in most of the world's watersheds and coastal waters [11].Therefore, the high intensities of NANI and NAPI might cause N and P to remain in the sediment, which will increase the N and P input load into the ocean.This N and P might accumulate long after inputs and continue to be released, by balance calculation of input and output flux, over the past 30-70 years [60].Strengthening the collaborative management of NANI and NAPI and residual N and P is the key to effectively controlling nitrogen and phosphorus pollution in the Pearl River.

Ecological and Management Implications
The riverine TN and TP fluxes were significantly correlated with NANI and NAPI via nonlinear correlation.Therefore, the input of NANI and NAPI should be controlled based on the driving force, to alleviate the deterioration of water quality and its impact on the oceans.The majority of environmental management studies on regulating N and P in intensive agriculture systems were geared toward reducing pollution in river basins [7].However, it is necessary to consider the NANI and NAPI input of regional heterogeneity.Human food and animal feed were the most important contributors of NANI and NAPI in the Pearl River Delta, while fertilizers N and P were the most important contributors in the non-Pearl River Delta region.Therefore, the investigation of pollutant sources and the implementation of source-control policies is indispensable.According to the contributors of NANI and NAPI of regions in the Guangdong Province, the management and control measures of N and P can be established in agricultural activities and personal consumption habits.In the Pearl River Delta, the reduction in the consumption of food N and P would benefit from actions taken to build a more sustainable social and environmental food system.In the non-Pearl River Delta region, measures such as precision fertilization should be taken to improve the fertilizer utilization rate and reduce fertilizer N and P usage.In addition, intensive water treatment facilities are required to improve discharge efficiency, to improve the water quality of the Pearl River Basin.

Conclusions
The data of average riverine nitrogen and phosphorus fluxes and net anthropogenic inputs of nitrogen and phosphorus from 2016 to 2020 in the Guangdong section of the Pearl River Basin were summarized and analyzed, in order to further effectively control nitrogen and phosphorus problems in the Pearl River.It provided a better understanding of anthropogenic inputs of N and P by the calculation of NANI and NAPI in the Guangdong section of Pearl River, where the intensities of NANI and NAPI were still higher than levels in China.The driving force of the intensity of NANI and NAPI in the Pearl River Basin was different from Pearl River Delta, coastal and northern regions.Agricultural activities were still the primary factor in changing NANI and NAPI intensities, however, the consumption of food/feed N and P had changed the main contributor in the Pearl River Delta due to rapid population growth and urban development.The high intensities of NANI and NAPI caused the Pearl River to deteriorate in water quality, which made ~50% of the Pearl River seriously substandard.Strengthening the collaborative management of NANI and NAPI is, therefore, the key to effectively controlling nitrogen and phosphorus pollution in the Pearl River.

Figure 1 .
Figure 1.The geographic location of the Pearl River in Guangdong Province.

Figure 1 .
Figure 1.The geographic location of the Pearl River in Guangdong Province.
urban 1795 (a, b) natural site 1140 (a, b) urban 1509 (a, b) natural site 912 (a, b) / urban 3304 ± 952 (a, b) natural site 2052 ± 1022 (a, b) is obtained by the chemical composition of precipitation; b by mixed deposition; d by dry deposition; w by wet deposition.

Figure 4 .
Figure 4.The distribution of Pearl River water quality assessment from 2016 to 2020.(a) TN; (b) TP The annual changes in NANI and NAPI intensities in the Guangdong section of the Pearl River Basin from 2016 to 2020 are shown in Figure 5.The intensities of NANI and NAPI tend to be stable.NANI showed a very slight trend of decrease from (1.51 ± 0.09) × 10 4 kg•N•km −2 •yr −1 to (1.36 ± 0.08) × 10 4 kg•N•km −2 •yr −1 , however, that was much higher than the Chinese average (3791 kg•N•km −2 •yr −1 ) [15].Compared with N, the intensity of NAPI remained stable, and the average intensity of NAPI in the Guangdong section of the Pearl River Basin was 3.8 × 10 3 kg•P•km −2 •yr −1 , but still at a moderate level compared with other watersheds (Figure 5b) [14,17].In the Guangdong section of the Pearl River Basin, fertilizer input was the most important component for the intensities of NANI and NAPI, accounting for 38.0-42.3% and 53.1-56.1%,respectively.The human and animal consumption (food/feed) inputs were the second dominating contributor of NANI, accounting for

Figure 4 .
Figure 4.The distribution of Pearl River water quality assessment from 2016 to 2020.(a) TN; (b) TP.
The annual changes in NANI and NAPI intensities in the Guangdong section of the Pearl River Basin from 2016 to 2020 are shown in Figure 5.The intensities of NANI and NAPI tend to be stable.NANI showed a very slight trend of decrease from (1.51 ± 0.09) × 10 4 kg•N•km −2 •yr −1 to (1.36 ± 0.08) × 10 4 kg•N•km −2 •yr −1 , however, that was much higher than the Chinese average (3791 kg•N•km −2 •yr −1 ) [15].Compared with N, the intensity of NAPI remained stable, and the average intensity of NAPI in the Guangdong section of the Pearl River Basin was 3.8 × 10 3 kg•P•km −2 •yr −1 , but still at a moderate level compared with other watersheds (Figure 5b) [14,17].In the Guangdong section of the Pearl River Basin, fertilizer input was the most important component for the intensities of NANI and NAPI, accounting for 38.0-42.3% and 53.1-56.1%,respectively.

Figure 5 .
Figure 5.The changes in each component of NANI and NAPI intensities from 2016 to 2020.(a) NANI; (b) NAPI.

Figure 5 .
Figure 5.The changes in each component of NANI and NAPI intensities from 2016 to 2020.(a) NANI; (b) NAPI.

Figure 6 .
Figure 6.Spatial and temporal changes of TN flux from 2016 to 2020.

Figure 6 .
Figure 6.Spatial and temporal changes of TN flux from 2016 to 2020.

Figure 7 .
Figure 7. Spatial and temporal changes of TP flux from 2016 to 2020.

Figure 7 .
Figure 7. Spatial and temporal changes of TP flux from 2016 to 2020.

Table 1 .
Atmospheric nitrogen deposition observation in the Pearl River Basin.

Table 5 .
Relational grade by grey relational analysis.