Spatiotemporal Dynamics of Nitrogen Budgets under Anthropogenic Activities in Metropolitan Areas

: The rapid growth of metropolitan regions is closely associated with high nitrogen (N) ﬂows, which is known as the most important reason for widespread water pollution. It is, therefore, crucial to explore the spatiotemporal patterns of N budgets under intensive human activity. In this study, we estimated the long-term (2000–2015) N budgets by integrating the net anthropogenic nitrogen input (NANI) and the export coefﬁcient model (ECM) in the Yangtze River Delta Urban Agglomeration (YRDUA), a typical metropolitan area with strong human disturbances. The results revealed that the NANI decreased by 10% from 2000 to 2015, while N exports showed a 6% increase. Hotspots for N budgets were found in the northeastern areas, where cropland and construction land were dominant. The linear regression showed a close relationship between the NANI and N export, and about 18% of the NANI was exported into the river system. By revealing the critical sources and drivers of N budgets over time, our work aimed to provide effective information for regional policy on nitrogen management. Future strategies, such as improving the fertilizer efﬁciency, optimizing the land use pattern, and controlling the population density, are necessary in order to address the environmental challenge concerns of excessive N.


Introduction
Over the past several decades, excessive nitrogen (N) enrichment in terrestrial and aquatic ecosystems has become a global concern [1][2][3]. The imbalance of the N cycle has led to a range of ecological damage, such as eutrophication, biodiversity loss, acid rain, and soil acidification, which is mainly due to the increasing human interference related to overfertilization, dense population, massive food demand, high livestock breeding rate, and increased combustion of fossil fuels [4][5][6][7]. Some researchers have even documented that the safe operating boundaries have already been exceeded in the systems of the nitrogen cycle that are under human interference, as well as in the rate of biodiversity loss and climate change; these have been identified as the three largest threats to our planet in the Anthropocene [8]. Such conditions could generate unacceptable environmental changes if effective strategies are not taken. Given that the global N overload continues to be a critical problem [1,[9][10][11], it is significant to understand and characterize the spatiotemporal patterns of the N budget, the main sources, and the impacts of human-induced N in high-urbanization areas.
While N inputs and exports are the two vital parts of the N cycle, research has proved that riverine N flux is greatly related to anthropogenic N inputs in many regions across the world [12][13][14][15][16]. The ratio of N exports to N inputs for rivers can also reflect the capacity for N retention in terrestrial ecosystems and the potential risks of N pollution in aquatic ecosystems [17][18][19][20][21]. A previous review showed that the proportion of N exports ranged from 3% to 118%, with an average ratio of 24% [20], and this was influenced by various factors, such as hydroclimate, land use type, the degree of N saturation, and other human activities [13]. For example, the N export ratio was lower in Baiyangdian Basin (3%), a semiarid watershed, than in Taihu Basin (38%), a humid watershed [22,23]. However, there have been few systematic reports on the dynamic spatiotemporal patterns of N budgets, as well as on the dynamic response of riverine N exports to anthropogenic N inputs or how N budgets have varied over time due to anthropogenic activities in metropolitan areas that have experienced rapid urbanization and socioeconomic development.
Revealing the spatiotemporal patterns of N budgets can provide effective information for improving nitrogen pollution management and can point to the dominant sources of pollution for potential reduction strategies [24]. The net anthropogenic nitrogen input (NANI) model, first introduced by Howarth et al. [25], has been widely used to estimate the anthropogenic N inputs in various regions across the world, and has been shown to be a simple and effective mass balance approach that is detailed enough to incorporate available data on individual crops, livestock, and people [6,12,13,[22][23][24][26][27][28][29][30][31][32]. This method is typically calculated as the sum of four major components: fertilizer N application, N in net food and feed imports, atmospheric N deposition, and biological N fixation. Considering the modes of delivery of N, Zhang et al. [28] divided the N inputs into non-point-source and point-source, providing a more accurate estimate of nutrient inputs, especially in areas with heavy point-source N pollution.
Given that increases in nutrient inputs are generally accompanied by corresponding increases in nutrient exports [12,13,22,23], it is necessary to assess N exports when evaluating the interactions between human activities and the N cycle. There is a series of models, including both mechanistic and empirical models, that can be applied to assess nutrient exports. Among them, the export coefficient model (ECM) is a strong empirical model that has been proved to be suitable for simulating the flux of nutrient pollutants entering rivers in plain areas [33][34][35], particularly when there are insufficient monitoring data for studies. This model, which was developed in North America in the 1970s, was confirmed to be able to resolve the issues related to the flat terrain and intricate river network, and can be used to estimate N exports into a river system when evaluating the nutrient contaminant export load [18,22,23,[36][37][38]. Recently, researchers studied the relationship between N inputs and riverine N exports by combining the NANI and ECM models in plain basins, such as the Taihu Basin [23] and Baiyangdian Basin [22]. Overall, the simulations for estimating nitrogen inputs and exports can be used to explore not only the exact N flow budgets, but also the corresponding risks of N load.
Urban ecosystems are complex systems that couple humans and nature. Exploring the effects of human activities in cities is of great significance for sustainable development, especially in the context of rapid urbanization. As one of the most developed urban areas in China, approximately 16% of the total population of China is concentrated in the Yangtze River Delta Urban Agglomeration (YRDUA), though it only accounts for 2% of the land area of China. Due to the humid climate and the flat terrain, the human activities for living and producing are very intense in this region. For instance, there are a large number of agricultural production activities in the central and northern parts, as well as dense population in the eastern and central parts. The Taihu Lake Basin in the center is widely regarded as an important area with excessive nutrient enrichment and a fragile water environment. Therefore, the YRDUA is a representative area with enhanced human disturbance in the context of rapid urbanization. Results from this region would also have implications for other urban areas, as most regions around the world are experiencing continuous socioeconomic development.
In this study, we integrated the NANI and ECM models to quantitatively evaluate long-term (2000-2015) N inputs and exports at the city scale in the YRDUA, a typical urban agglomeration in China, and identified the spatial and temporal dynamics, as well as the potential driving effects of socioeconomic factors on N budgets. The specific goals were: (1) to quantify the N inputs and exports and identify their main sources; (2) to reveal the dynamic spatial-temporal characteristics of N budgets and identify the impact of anthropogenic driving factors; and (3) to explore the N export ratio and to illustrate the response of riverine N exports to anthropogenic N inputs.

Study Area
The Yangtze River Delta Urban Agglomeration (YRDUA), located on the largest alluvial plain in Eastern China (29 •   As one of the six world-class urban agglomerations, the YRDUA has experienced fast urbanization and industrialization in the past few decades. Population settlements in this area have reached 1063 people km 2 , which is approximately seven times the nation's average. The land cover mainly consists of cropland and forest land, accounting for 51% and 25% of the total area, respectively. Overall, cropland is mainly distributed in the northeast, while forest land is distributed in the southwest. In addition, construction land is densely concentrated in the central and eastern areas (especially in the Taihu Lake Basin) ( Figure S1). The large proportion of agricultural land, the high-density population, and the development of industry and transportation have caused a significant influence on the ecological environment of this region. In addition, there is a dense river system, with more than 200 lakes [39]. Taihu Lake, the third-largest freshwater lake in China, is located in the center of this area. Large-scale eutrophication is a serious problem in the surface waters of the region, and especially in Taihu Lake [40]. Currently, the contradiction between the human interference and ecological security has yet to be effectively resolved.

Accounting for Anthropogenic N Budgets
To estimate the human-controlled nitrogen budgets entering the territory, we used the NANI (net anthropogenic nitrogen input) approach to calculate the anthropogenic N inputs and the ECM (export coefficient model) to count the N exports. An overview of the model framework and data requirements is presented in Figure 2. Additionally, the methodology for evaluating N budgets is described below.

Non-Point-Source Inputs
Non-point-source N inputs (NANI n ) are mainly composed of atmospheric N deposition, fertilizer N application, net food and feed import N, and biological N fixation. The basic calculation (Equation (1)) is presented below: where NANI n is the N from non-point-source inputs, N fer is the N fertilizers applied, N fix is the biological N fixation, N dep is the atmospheric N deposition, and N im is the net food and feed N imports in urban and rural regions. The units are kg N km −2 a −1 . Among these components, only N im is allowed to be negative because its calculation is based on the balance of local production and consumption [41]. It is negative when the local food and feed supply exceed the demand. The net food and feed imports can be calculated with Equation (2): (2) where N selfo and N selfe are the N consumption by humans and livestock, respectively. N harv is the N in crops, and N liv is the N in animal products. The units of the above are kg N km −2 a −1 .

Point-Source Inputs
Point-source N inputs (NANI p ) are estimated as the sum of industrial and urban point discharges that could directly flow into the river systems [28]. The calculation (Equation (3)) is where N ind is the N discharged by industrial production, N urban is the N discharged by urban inhabitants, and the units are kg N km −2 a −1 . I rem−tn refers to the average N removal rate by a sewage plant. I sew is the proportion of sewage effluent connected with sewage plants. The above are expressed as percentages.
Based on the assumption that N emitted in wastewater by households and industries is generally connected to the same sewerage system [24], the calculation of I sew can be obtained from the following equation (Equation (4)): where W ind and W urban are the volume of wastewater produced by industrial production and domestic households, respectively; W sew is the actual treatment volume by sewage plants.

Estimation of N Exports
To understand the nitrogen exports in different subsystems, we applied the export coefficient model (ECM) to estimate N exports. The ECM model is a mathematical weighted equation that can be applied to assess the pollution load exported from different sources to water bodies [38]. Here, we estimated the nitrogen exported from agricultural land, forest land, grass land, and other land-use types, as well as livestock breeding, urban life, and rural life. The expression (Equation (5)) is where L is the amount of nitrogen exported (kgN a −1 ), and E i is the export coefficient for each source (kg N km −2 a −1 ). A i is the area of the land-use type (km −2 ) or the number of livestock (capita) or population (people). The export coefficients used here (Table S4) refer to previous studies [22,23].

Data Sources
City-level datasets, such as population, crop yields, amounts of fertilizer applied, livestock number, and industrial production, were adopted for counting anthropogenic nitrogen inputs. Most of the data used to calculate N budgets or represent the human activity index were obtained from local statistical yearbooks, national agriculture yearbooks, and national economic and social statistical bulletins. The raster data of wet and dry atmospheric N deposition from 2000 to 2015 were derived from Jia et al. [42] and can be accessed at http://www.sciencedb.cn/dataSet/handle/607 (accessed on 12 December 2020). The land-use data set (2000-2015) originated from the Data Center for Resource and Environmental Science, Chinese Academy of Sciences (RESDC) (http://www.resdc.cn (accessed on 12 December 2020)). Additionally, the calculation parameters applied in this study refer to previous research and are shown in the Supplementary Materials.

Data Analysis 2.4.1. The Relationship between N Budgets and Anthropogenic Activities
To assess the effects of human activities on nitrogen inputs and exports, we used the method of linear regression. Three categories of human activity indexes, including nine sub-indexes, were chosen to establish the human activity index system (Figure 3), which consisted of the population density (PD), the urbanization rate (UR), the gross domestic product (GDP), the percentage of agricultural land (PAL), the percentage of forest land (PFL), the percentage of construction land (PCL), the proportion of primary industry (PPI), the proportion of secondary industry (PSI), and the proportion of tertiary industry (PTI).

Hotspot Analysis
Hotspot analysis was conducted via the "Hot Spot Analysis (Getis-Ord G * i )" tool in the Spatial Statistics Tools, which is a part of the ArcGIS software 10.2, and was used to identify the areas of the greatest concern [24]. The formula of Getis-Ord G * i calculated for each feature in a dataset is given as (Equation (6)): where x j is the value for feature j, n is equal to the total number of features, W i,j is the spatial weight between feature i and j, and For statistically significant (p < 0.01 or p < 0.05 or p < 0.1) positive G * i , a larger G * i value indicates that the clustering of high values (hotspot) is more intense. For statistically significant (p < 0.01 or p < 0.05 or p < 0.1) negative G * i , smaller G * i means more intense the clustering of low values (cold spot).

Sensitivity Analysis
Sensitivity analysis is an effective method for assessing the contributions of the main components of N budgets [12]. The sensitivity (S) of a variable y to a parameter x is calculated by examining the effect of a change in x on the response of y relative to a baseline value. The sensitivity is defined here as the proportional change of variable y, relative to baseline y b , divided by the proportional change in parameter x, relative to baseline value x b . In this study, we applied a ±10% change for the main parameters. The expression (Equation (9)) is where x is the parameter, and y is the variable; x b and y b are the baselines of x and y, respectively.

Temporal Dynamics and Sources of Anthropogenic N Inputs
On average, the area-weighted NANI in the YRDUA decreased from 18,766.79 kg N km −2 a −1 in 2000 to 16,798.08 kg N km −2 a −1 in 2015, showing a downward trend ( Table 1). The NANI values of the 26 cities in this region ranged from −10,563.83 to 54,168.68 kg N km −2 a −1 , with an average of 19,331.37 kg N km −2 a −1 in 2000 and, in 2015, ranged from −8612.49 to 39,208.17 kg N km −2 a −1 , with an average of 17,392.74 kg N km −2 a −1 . This is due to the distinctions in the urban development characteristics, which caused such a large difference in the NANI intensity between these cities. Non-point source N accounted for approximately 90% of the total N inputs, indicating that the NANI in the YRDUA mainly came from non-point-source inputs. Over the study period, the change in total anthropogenic N inputs was consistent with the non-point-source N inputs, while the point-source N inputs showed an upward trend first, and then decreased. Among the components of the NANI, the fertilizer N application was proved to be the dominant source for N inputs across the 26 cities studied (Figure 4), contributing about 58%, which implied that agricultural production was one of the main drivers of N input intensity. Accompanied by intensive fertilizer application, the atmosphere N deposition, which became the second major source of N inputs, presented an increasing trend, and its amount increased by 29% from 2000 (2518 kg N km −2 a −1 ) to 2015 (3238 kg N km −2 a −1 ). Domestic sewage N discharge was the third major source of the NANI in the YRDUA, with an average proportion of about 10%. In contrast, the N contents from biological N fixation and industrial N production were relatively stable. The variation rates of N input sources over the study period are shown in Figure 5. During 2000-2005, four out of six components of the NANI increased. However, in 2005-2015, there was only one increase; in other words, except for the atmosphere N deposition, the others showed a downward trend. Thus, we concluded that, in addition to the fertilizer N application, which is generally the main source of N inputs, atmospheric N deposition, which mainly results from fossil fuel combustion, traffic emissions, industrial exhaust emissions, etc., also requires sufficient attention in metropolitan areas.

Temporal Dynamics and Sources of Anthropogenic N Exports
The dynamic characteristics of N exports to rivers in the YRDUA during 2000-2015 are shown in Figure 6. In this period, area-weighted N exports grew from 3047.52 kg N km −2 a −1 in 2000 to 3233.46 kg N km −2 a −1 in 2015, with an increase of 6%. Various sources were calculated, including farming, livestock breeding, rural life, urban life, and other land-use types. We observed that the contributions from the main pollution sources changed considerably between 2000 and 2015. Agricultural production was the predominant source of N exports in the YRDUA, contributing to 35% of total N exports, and this amount showed a downward trend in 2000-2015, which was consistent with the reduction of cultivated land, which fell by 7%.
On the contrary, the amount of N exported from urban life, the second main source, showed a significant increasing trend, rising about 106% from 2000 to 2015. This is due to the rapid increase in urban population-more specifically, the number of urban residents in 2015 was about twice as many as in 2000. The N exported from both rural life and livestock breeding showed a decreasing trend from 2000 to 2015, in which the trend of rural life was more obvious. Forest land, grass land, and other land-use types remained relatively stable during the study period. Overall, changes in N exports were inextricably linked to urbanization, as well as the intensification of agriculture and aquaculture.

Overview of N Flows in the Metropolitan Area
An overview of the mean N flows in the YRDUA during 2000-2015 is presented in Figure 7. According to the traces of N flows, we found that the agricultural system was the most crucial part, and received the highest N mass compared with other nodes (e.g., urban and rural inhabitants and livestock). However, the utilization efficiency of N showed a considerable mismatch with the intensity of fertilizer N inputs into this system, which resulted in almost 70% of N inputs not being effectively used. The surplus N was either exported to the water bodies through a non-point pathway or retained in the system. As for the human and livestock systems, the N derived from urban and rural residents and livestock was apparently also the major component of non-point-source pollution. Even if we assume that the N production of local crops and animals was entirely used for domestic demand, this metropolitan region still needed food and feed N imports at the overall level. Notably, several assumptions were made to generate this diagram (Figure 7), including that (1) we assumed that each node (i.e., crops, humans, livestock, etc.) was subjected to the conservation of mass; (2) the N produced by local crops and animals were assumed to be completely used to meet the demands of local human consumption and livestock breeding; (3) according to the difference in N demand between people and animals, we assumed that 70% of the crop production N flowed to livestock, while 30% was supplied to humans; (4) the percentage of sewage effluent that was treated by sewage plants was assumed to be 80% based on a multi-year average.

Spatial Patterns and Hotspot Analysis of N Budgets
The dynamic spatial distribution pattern of N budgets is shown in Figure 8. The overall distribution of the NANI presented a similar spatial pattern to that of non-pointsource N inputs (NANI n ) during 2000-2015, showing that the N inputs were higher in the northeastern plain area than in the southwestern hilly region. It appears that the N inputs were generally high in areas where the population density was high and were low in areas where woodland was the main land-use type. We found that the NANI intensities in Shanghai and Jiaxing were invariably high throughout this period; however, the reasons for their high values are not the same. Specifically, net food and feed N import contributed largely to the N inputs of Shanghai, while fertilizer N application was the most important source for Jiaxing. This mostly resulted from the different development patterns. Shanghai is the city with the highest population density in the YRDUA, which has led to a great demand for food N imports, and this is different from Jiaxing, which has relatively developed agricultural production. The N exports also showed a similar spatial trend of N inputs, in which Shanghai and Jiaxing were also the most prominent. Considering the distribution of the land-use patterns ( Figure S1) and population residences, we concluded that high-intensity N budgets mainly occurred in areas with intensive agricultural activities and dense population.
The cities with high N budgets in the YRDUA were mainly concentrated in the Taihu Lake Basin and its north, including Shanghai, Jiaxing, Yancheng, Taizhou (JS), Nantong, Changzhou, Wuxi, and Nanjing. These cities were also identified as the hotspots (p < 0.01) of N budgets during the study period ( Figure S2). Through our research, the hotspots covered 19% of the total land, but contributed to more than 30% of the total N budgets. In contrast, the cold spots occupied 14% of the total area and contributed less than 7% of the total N budgets. Overall, the spatial distribution of hotspots remained relatively stable during 2000-2015, implying that the relatively high values of N inputs found in 2000 were likely to remain in 2015.

Potential Driving Factors of N Budgets
To explore how human activities affected the N budgets in the YRDUA, linear regression analysis was performed to explore the relationship between the human activity index and N budgets ( Figure 9). All data have been standardized here. Three categories of indicators, including the levels of population and economic development (PED), land-use pattern (LUP), and industrial structure (IS), were selected in this analysis (Figure 3). The results showed that, except for the urbanization rate (UR), the other PED indexes were significantly positively correlated with the N inputs and exports (Table S5). This indicates that cities with higher population densities or higher levels of economic development are more likely to have higher N input and export intensities.
Among the LUP indexes, the percentage of agricultural land (PAL) and the percentage of construction land (PCL) had clear positive effects on the intensity of nitrogen flux, while the N inputs and exports were inversely proportional to the proportion of forest land (PFL). Consistently with the spatial distribution pattern (Figure 8), the areas with high N budgets were located in the northeast of the YRDUA, where there was a high proportion of cropland or construction land. The low N budgets were distributed in the southwestern cities, where the land-use patterns were dominated by forest areas. Regarding the IS indexes, the three sub-indexes (PPI, PSI, and PTI) were observed to have no obvious correlation with the N inputs and exports, indicating that the characteristics of the industrial structure had no direct driving effects on the nitrogen flux in urban areas.

Dynamic Response of N Exports to N Inputs
The riverine N exports from the 26 cities were significantly and positively correlated with the NANI during 2000-2015 in the YRDUA (Figure 10). The slopes between the NANI and riverine N exports were 0.14, 0.16, 0.19, and 0.21 in 2000, 2005, 2010, and 2015, respectively. The slope of the fitted equation showed a significant increasing trend from 2000 to 2015. This is consistent with a previous study on the Taihu Lake Basin, in which the export ratio also showed remarkable growth during 1980-2010. Overall, 18% of the N inputs were exported to the water bodies across the 26 cities in the YRDUA, which was evidently higher (about 3%) than in the Baiyangdian Basin [22]. Considering the case that they are both plain areas, hydroclimate might be a major cause of this difference. Therefore, it can be concluded that, compared with arid or semi-arid areas, the NANI could have a greater influence on the riverine N in humid and rainy regions.
As for the NANI n and NANI p , the N export ratio also had a distinct rising trend during the study period. On average, 19% of the NANI n and 86% of the NANI p were exported into water bodies. Even though the non-point-source inputs were fairly high, they are usually relevant to strong landscape retention processes, which can dramatically reduce the masses flowing into river systems.

Anthropogenic N Inputs in the YRDUA
In this study, long-term N inputs in the YRDUA were quantitatively evaluated by using the NANI method. For further understanding, we compared the characteristics of the NANI n in major basins and regions around the world [12,16,[22][23][24]26,31,[43][44][45][46][47][48][49] (Figure 11). We only compared the characteristics of non-point-source N here, as it was the dominant source in this study and there have been relatively few studies involving point-source nitrogen so far. In general, the multi-year averaged NANI n in the YRDUA was around three times higher than that reported in mainland China [44,45].
Except for the Huai River Basin (HRB) and Baiyangdian Basin (BYD), the NANI n of the YRDUA was higher than those of the other basins or regions shown in the figure. Fertilizer N application and atmospheric N deposition were the main causes for the relatively high values of the NANI n in these two basins, and this is consistent with the two main sources of the NANI in the YRDUA. Unlike other areas where fertilizer N was the most vital source of N inputs, the atmospheric N deposition was the dominant source in Beijing, which was closely related to the urban characteristics there. This also revealed that fossil fuel combustion could be the main source of N pollution in metropolises.
Generally, land-use types dominated by cultivated land were largely responsible for large N inputs [50]. For instance, the percentage of agricultural land (51%) in the YRDUA was smaller than that in the Huai River Basin [24], and the NANI n was comparatively lower due to the relatively less intensive agricultural production and fertilizer application (Table S6). In contrast, the NANI n in the YRDUA was higher than that in the Erhai Basin [48], where forest was the main land-use type, covering about 69%, which is nearly three times the proportion of forest land in the YRDUA. In addition, the high population density also contributed a great deal to the N input intensity, which was demonstrated by the examples of the Taihu Basin, Chaohu Basin, Beijing Municipality, and YRDUA ( Figure 11). The uncertainties in our study are mainly relevant to the data sources and parameters because the N budgets were calculated based on two empirical models. First, most of the data were originally derived from statistical yearbooks, which are known to be the most reliable data source in China [51]. Many parameters used in this study for assessing N budgets through individual crops, livestock, population, etc. can also cause uncertainties. The most influential parameters of N budgets were identified through sensitivity analysis ( Figure S3). The results showed that the most sensitive input term of the non-point-source components of the NANI was fertilizer application, followed by food and feed inputs. Regarding the point sources, the most sensitive component was the urban domestic N discharge. Ensuring the accuracy of sensitivity parameters is essential for reducing the biases and uncertainties of estimation. In addition, the dynamics of the coefficients might also result in potential uncertainties. Thus, it is essential to be careful when selecting the associated parameters.

Influence of Anthropogenic Inputs on Riverine N Exports
Studying the response of riverine N exports to human-caused N inputs can promote understanding of the influence of anthropogenic N inputs on regional water pollution at the macroscopic scale. We examined the N export ratio and found that a mean of about 18% of the N inputs were exported to water bodies during 2000-2015. The N export ratio showed an increasing trend from 2000 to 2015, with a growth rate of 7%. Compared with the global N export ratio (25%) [52], that of the YRDUA was significantly lower, which indicates that there might be a considerable amount of N retained in this urban agglomeration, with potential harm to the terrestrial ecosystem.
For non-point-source and point-source pollution, although non-point-source N inputs were dominant in the YRDUA (contributing about 90%), we found that only about 19% of non-point-source N was exported to water bodies. This is because non-point-source N inputs are generally related to landscape retention processes, which can significantly reduce the N masses [24], and the denitrification may also be a significant loss pathway in agricultural watersheds. However, attention should also be paid to the issue of nitrogen saturation. As mentioned by Lian et al. [23], the excess N could not be retained in the Taihu Basin because the N surplus there was so large. This indicates that the ecosystem has a certain storage capacity for nitrogen; once it oversteps its capacity, the N retention rate of a landscape will be greatly reduced. Regarding the limited amount of green space (forest land and grass land), which typically has a higher N retention rate than other land-use types, in the YRDUA (accounting for only 29% of total land in 2015), it is vital to investigate the storage capacity and nitrogen saturation limit in rapidly urbanizing areas in the future.

Possible Implications for Reducing N Loading
The spatial-temporal characteristics and main sources of N budgets in the YRDUA, a typical developed urban agglomeration, were explored in our study. It is crucial to take effective measures and strategies for mitigating N pollution. Although the overall anthropogenic N inputs showed a downward trend, the contribution of the increasing atmospheric N deposition cannot be ignored. In addition, areas identified as hotspots ( Figure S2) were always areas of high N inputs and exports during the study period, which should be paid more attention in order to control excessive N. Regarding the agricultural system, which was identified as the most important node of N flows, it is necessary to improve the production efficiency of limited cropland-for instance, by carrying out different crop rotations and increasing the effective yield of crops-as production is intensified on land that is not converted for urban uses.
In addition, the scientific application of chemical fertilizers, improvement of fertilizer use efficiency, or development of biological fertilizers (e.g., green manure and manure piles) can also contribute to reducing the N load of the farming system. Regarding the livestock and human systems, effective use of their N outputs, such as human and animal feces or animal remains, is a useful way to promote N recycling. For industry, excessive dependence on fossil fuels greatly promotes the increase in atmospheric N deposition; therefore, it is necessary to use strategies such as alternative and renewable energy sources, including solar, wind, and energy derived from biomass.
Since sewage treatment plants play an important role in reducing the point-source pollution, they should also receive attention when attempting to control N. According to the sensitivity analysis of N inputs (Figure S3), if the sewage treatment capacity or N removal rate can increase by 10%, the point-source inputs of the NANI can be reduced by 9%. Therefore, improving the sewage treatment could be significant for the reduction of point-source inputs in the future.
Furthermore, optimizing the land-use patterns will also make a great contribution to the reduction of the nutrient load. As mentioned above, the proportion of forest land has a significant negative impact on the intensity of nitrogen budgets. If land-use patterns can be managed and optimized appropriately, with urban planning guiding sustainable development and natural environments that serve as a nutrient sink where nitrogen can be consumed by biological processes, the ecological damage from anthropogenic nitrogen sources can be minimized.

Conclusions
In this study, we evaluated the dynamic spatial-temporal characteristics of N budgets, the N export ratio, and the main anthropogenic driving factors of N inputs and exports in the YRDUA, a developed urban agglomeration area with intensive human activities. On average, the value of the NANI was about three times higher than the counterpart reported in mainland China over the study period. While non-point-source N contributed a large portion of the total inputs, fertilizer N application was the most important component of anthropogenic N inputs.
Although fertilizer nitrogen input was controlled to a certain extent from 2000 to 2015, the significant increase in the atmospheric nitrogen deposition revealed that controlling the use of fossil fuels and the emission of exhaust gas had great significance for the nitrogen management in rapidly urbanizing regions. For the point-source inputs, although they could be effectively controlled with the improvement of sewage treatment capacity, when this ability develops to a certain level, the continuous growth of the population will bring huge pressure on the regional nitrogen control. As was shown in the analysis of potential driving factors, the population density had a prominent positive correlation with the nitrogen budgets, especially the point-source nitrogen input. Through the analysis of the riverine nitrogen export ratio, it was found that around 82% of the nitrogen remained in the ecosystem, which may cause nitrogen saturation that will need to be addressed in the future.