Interaction among Controlling Factors on Riverine DIN Export in Small Mountainous Rivers of Taiwan: Inseparable Human-Landscape System

: Increasing anthropogenic nitrogen (N) emission via di ﬀ erent pathways has shown prominent impact on aquatic ecosystems for decades, but the e ﬀ ects of interaction among climate-, landscape- and human-associated variables on riverine DIN (dissolved inorganic nitrogen, mainly NO 3 − and NH 4 + ) export are unclear. In this study, the data of 43 watersheds with a wide range of climate-, landscape-and human-associated gradients across Taiwan were evaluated with partial redundancy analysis (pRDA) to examine their interactive controls on riverine DIN export. Results show that the annual riverine DIN export in Taiwan is approximately 3100 kg-N km − 2 yr − 1 , spanning from 230 kg-N km − 2 yr − 1 in less disturbed watersheds (eastern and central Taiwan) to 10,000 kg-N km − 2 yr − 1 in watersheds with intensive human intervention (southwestern and northern Taiwan). NO 3 − is generally the single dominant form of DIN, while NH 4 + renders signiﬁcance in disturbed watersheds. Nearly all environmental variables display a positive correlation with DIN export, except for landscape setting variables (e.g., slope, area, channel length), which show a negative relationship. In terms of seasonal pattern, climate and human-landscape variables are related to NO 3 − export independently in the wet season, yet in the dry season climate-human variables jointly dominate NO 3 − export. Meanwhile, human-landscape (LH) variables ( λ 1 of LH > 0.60) control NH 4 + exports in both seasons, and human-associated (H) variables ( λ 1 of H = 0.13) have a minor e ﬀ ect on NH 4 + exports in dry season. Precisely, the contribution of controlling variables on DIN export vary with species and seasons, indicating water quality management could be time-dependent, which should be taken into consideration for designing mitigation strategies.


Introduction
Reactive nitrogen, a vital and essential nutrient for organisms and ecosystems, plays a key role in maintaining biodiversity and functions of ecosystems [1,2]. Over the past half century, the rapid increasing anthropogenic N emissions inevitably accelerated N deposition into the biosphere [3], and consequently exceeded the N-requirement for terrestrial ecosystems [4]. Studies show that N emissions and depositions have been declining in Europe and the U.S. since 2000, East and South Asia, in contrast, have become the hot spots of pollutant emissions due to population growth and intense agricultural activities [5,6]. Regions located near the emission source of East Asia receive excessive N deposition under favorable climatic conditions, such as the East Asian monsoon that blows directly to Taiwan with abundant rainfall [7]. The overloaded dissolved inorganic nitrogen (DIN) (majorly NO 3 − , and NH 4 + ) led to eutrophication and harmful algal blooms which deteriorated water quality and caused damages in aquatic communities [8]. However, the DIN exports, particularly for seasonal changes, in subtropical mountainous watersheds are still unclear. Many studies have demonstrated that even a mild replacement of natural vegetation with agricultural land use within watersheds would have significant impacts on hydrochemical processes [9,10], especially for DIN export from non-point pollution sources [11][12][13]. Landscape features such as slope, soil type/moisture, channel length, watershed area and relief may also regulate water quality [14]. In addition, climatic factors play a principal role in nutrient cycling in the era of warming climate and increasing extreme events. A warmer temperature will accelerate biogeochemical processes that would alter enzyme reactions, e.g., nitrification and denitrification, via microbial activities [15]. The torrential rainfall caused by synoptic weather, such as thunderstorms and tropical cyclones during summer, can bring approximately 30% to 50% of annual precipitation, such that a considerable amount of DIN will be flushed out, as evidenced in Taiwan [11,12]. However, the effects of the independent and interactive relationship among landscape patterns, climatic factors and anthropogenic disturbance are still not clear on riverine DIN export [16].
Previous analyses widely used linear or non-linear regression models to estimate the DIN export, but failed to separate the collinearity among variables, which might violate the statistical independent assumption [17]. To tackle the issues of collinearity, a series of well-developed methods, namely, principal components analysis (PCA), redundancy analysis (RDA) and partial RDA (pRDA), have been utilized to clarify the relative importance of independent variables and their interactive effects on dependent variables [18,19]. In this study, we apply PCA and pRDA to evaluate the contribution of climatic-, landscape-and human-associated variables and their interactions on DIN export based on 43 island-wide watersheds in Taiwan. This synthesis can help to disentangle the intertwined effects of these fundamental environmental factors on the behaviors of nutrient fluxes. Specifically, the objectives of this study are to (1) quantify riverine DIN export in subtropical watersheds in Taiwan, (2) explore the interplay of human disturbance, climatic factors and landscape characteristics on riverine NO 3 − and NH 4 + exports, respectively, and (3) figure out the spatial and seasonal variation of the controlling factors' influences on DIN exports.

Study Site
Taiwan, located in the northwest Pacific Ocean, is a tropical/subtropical mountainous island in the East Asian monsoon climate zone. Elevation ranges from sea level to approximately 4000 m in a short horizontal distance (<75 km). The Tropic of Cancer crosses the central part of Taiwan, which divides the island into tropical monsoon climate in the south and subtropical monsoon climate in the north. The mean annual temperature (MAT) is 22 • C across the island, ranging from 15.0 • C/18.2 • C in the north/south in January to 28.7 • C/28.4 • C in July, and the MAT decreases with increasing altitude [20]. The mean annual precipitation (MAP) is 2500 mm for the entire island but shows high spatial variability ranging from less than 1500 mm in southwest Taiwan to over 4000 mm in the mountains of northeastern Taiwan [20]. There are more than 75% of MAP falls in the humid summer (May to October), while winter-spring (November to April) is a relatively dry period.
Natural forests, plantation and bamboo forests cover 60% of the land, while farmlands and urbanized areas occupy 29% and 11%, respectively ( Figure 1). The dominant vegetation types and land cover of Taiwan change from mixed and conifer forests at the mid and high elevation (>2000 m a.s.l.), evergreen broadleaf forests at the low and mid elevation (200−2000 m a.s.l.) to urban, buildup area and farmland on coastal plains (<800 m a.s.l.) (Figure 1). In this study, in order to estimate island-wide watershed DIN export, our sampling scheme covers a total of 43 sites. Individual watershed represented by each site can be as small as 82 km 2 and as large as 2969 km 2 , and as a whole, they occupy more than 70% of the island and distribute evenly across Taiwan (Figure 1). Twenty-nine of the 43 watersheds are majorly covered by natural forest (>70% of forest cover), while 14 of the 43 are situated in the transition between mountain and plain regions where agricultural land cover and buildup area account for 10-65% and 0-18%, respectively ( Figure 1 and Table S1). The average slope of the watersheds varies from 8% to 76%.

DIN Sampling and Streamflow Simulation
A monthly sampling scheme was conducted at all sampling sites (watersheds) during 2015-2016. Each route would be finished within 2 days. Stream water samples were collected by plunging a 1-L PE (polyethylene) bucket into stream, and the water was immediately filtered through 0.7 µm filters. A 15 mL subsample of filtrate was frozen on-site in liquid nitrogen and kept frozen until laboratory analysis at National Taiwan University. Nitrate, nitrite and ammonium content were analyzed using ion chromatography (IC) using a Dionex ICS−1500 (Thermo Fisher Scientific Inc. ® Sunnyvale, CA, USA) with a detection limit of 0.2, 0.2, and 0.4 µM, respectively. Our DIN calculation includes NO 2 − , NO 3 − and NH 4 + . Among them, nitrite is easily oxidized to nitrate and accounts for a small fraction (<5.0%). Therefore, we mainly analyzed and discussed NO 3 − and NH 4 + .
Streamflow was acquired from the water level stations maintained by the Water Resource Agency (www.wra.gov.tw). The water level records were converted to streamflow via an individual rating curve and cross section approach. For some periods of missing records and ungauged sites, a hydrologic model (HBV, the Hydrologiska Byråns Vattenbalansavdelning model [21]) was used to fill the data gaps [22]. The historical observed daily streamflow was utilized to train the parameter set to fit low, normal and extreme values of simulated streamflow using the performance measure of NSE (Nash-Sutcliffe efficiency coefficient) [23]. The calibrated parameter set was then applied to the watersheds using their own climatic inputs and terrain information to simulate their daily streamflow during 2015-2016.

Export Estimation
Based on the discrete NO 3 − and NH 4 + concentration and continuous streamflow rates, individual NO 3 − and NH 4 + export of the 43 sites were then estimated using an R software package, loadflex, which provides several common methods (e.g., interpolations, regressions and composite method) for export estimation [24]. The composite method synthesizing rectangular interpolation and regression models is applied for export estimation. The mean values of two export methods in 2015 and 2016 were used as export results for further analysis. The rectangular interpolation has usually been used for studies on solute and sediment exports, where horizontal lines are drawn through observations in a plot of concentrations against time, and each horizontal line is connected to the next by a vertical line midway between successive observations [25]. The regression approach is a longstanding interpolation in estimating watershed solute exports ( Figure 2). It often requires less data than other models if the data can span over the range of predictors instead of the full time period of interest [26]. A simple regression equation based on observed nutrient concentrations and concurrent streamflow (Q) with an exponent function is used to represent the hydrological influence on transport as Equation (1) below: where Q j [mm d −1 ] is the daily streamflow rate on j-th day; C j [mg-N L −1 ] is an estimated concentration of NO 3 − and NH 4 + on the j-th day, m is the conversion factor to convert the calculated values into a specific unit [kg-N km −2 yr −1 ], and a 0 and a 1 are regressive coefficients. Coefficient a 0 is generally highly associated with the mean of observed nutrient concentration, and a 1 indicates the hydrological influence. A larger coefficient, a 1 (>0), indicates enhanced concentration with increasing streamflow, whereas a smaller value reflects the dilution effect because concentration decreases with the increase of discharge percentage. From Equation (1), we can estimate the concentration and export for non-measured days by introducing continuous daily streamflow [27]. According to the hydrologic seasonality, we summarized the daily export from May to October as wet season export and the summation of other daily exports as dry season export.

Variation Partitioning: PCA, RDA and pRDA
To tackle collinearity issues, many researchers have used principal components analysis (PCA) to reduce a number of correlated variables into a set of uncorrelated variables, which reserves its total variance and uncovers its hidden patterns [28]. In addition, in order to realize the relative importance of different explanatory variables and their interactive effects, direct gradient analysis such as redundancy analysis (RDA) and its successive partial constrained ordinations, i.e., partial RDA (pRDA) have been commonly proposed [18,19]. This allows researchers to explore the relationships between predictor variables and dependent variables by removing the intertwined effects among them [29]. However, most previous studies utilized PCA or RDA methods focusing on the relationship between biological phenomena and environmental influence, and there are only a few studies on water quality (e.g., [30]). Therefore, exploring the likely collinear controlling factors and their internal relationships to riverine DIN transport based on PCA and RDA will be valuable.
In order to interpret the spatial and temporal patterns of riverine NO 3 − and NH 4 + export and unravel the dependency among controlling factors, i.e., the human disturbance, climatic factors and landscape settings (Table 1), the whole analysis was carried out in three steps [29]: (i) the PCA was applied to find out a set of uncorrelated variables, (ii) detailed relations between export and each one of the controlling factors were displayed using a scatterplot matrix, and (iii) the RDA and pRDA were conducted to disentangle the contribution of the major variables.

SLP100
The average slope in the 100 m buffer zone SLP200 The average slope in the 200 m buffer zone SLP500 The average slope in the 500 m buffer zone SLP1000 The average slope in the 1000 m buffer zone SLP2000 The The percentage of buildup area in the 100 m buffer zone BD200 The percentage of buildup area in the 200 m buffer zone BD500 The percentage of buildup area in the 500 m buffer zone BD1000 The percentage of buildup area in the 1000 m buffer zone BD2000 The percentage of buildup area in the 2000 m buffer zone BD The percentage of the buildup area in watershed The percentage of agriculture in the 100 m buffer zone AGR200 The percentage of agriculture in the 200 m buffer zone AGR500 The percentage of agriculture in the 500 m buffer zone AGR1000 The percentage of agriculture in the 1000 m buffer zone AGR2000 The percentage of agriculture in the 2000 m buffer zone AGR The percentage of the agriculture in the watershed PCA was applied to reduce redundant information and to transform the original correlated data into another set of uncorrelated variables. The PCA keeps only a few independent sets (patterns) of environmental data that are distinct from each other, which will help to realize the effects of various characteristics of watersheds on NO 3 − and NH 4 + exports in our study [29,31]. The varimax rotation was selected to better separate divergent groups of variables, as suggested [32]. The environmental variables were centered and standardized in order to approximate normally distributed random errors and then were derived from the PCs via a standardized linear projection which maximizes the variance in the projected space [33]. For a set of observed-dimensional data vectors, { t n }, n ∈ {1, · · · , N}, the q principal axes { w j , j ∈ 1, · · · , q , could be derived as the orthonormal axes onto which the retained variance under projection is maximal. It can be shown that the vectors w j are given by the q dominant eigenvectors (i.e., those with the largest associated eigenvalues λ j ) of the sample covariance matrix. The outcomes of PCA help us to identify relationships between these variables and determine which variables require further investigation. The variables with loading higher than 0.1 in the first and second PCs were kept for the following RDA and pRDA analysis to constrain the ordination of environmental variables and to avoid the collinearity problem [31,34]. Moreover, we know water quality is regulated by riparian zones along the river and stream networks, but what needs to be clarified is spatially to what extent their individual effect is [35]. Here, we delineated the buffer zones of 100, 200, 500, 1000 and 2000 m along the stream network using the buffer tool in ArcGIS v.10.7. (ESRI Inc., Redlands, CA, USA) The environmental variables within the entire watershed and five buffer zones were also retrieved as previous studies suggested [30,36,37]. The land cover/land use data were acquired from the Ministry of the Interior of Taiwan (Figure 1b; https://www.moi.gov.tw), and the digital elevation model (DEM) data were derived from the open data platform in Taiwan (https://data.gov.tw/), which were provided as input for calculations of landscape settings and human disturbance variables ( Table 1).
The RDA and pRDA were further applied to quantify the individual effect and integrative contribution among human disturbance, climatic factors and landscape setting on riverine NO 3 − and NH 4 + exports [18]. RDA extends the algorithm of PCA with a response matrix Y (with n objects and p variables) by an explanatory matrix X (with n objects and m variables). First, RDA produces a matrix of fitted valuesŶ through Equation (2) and second, runs a PCA based onŶ [38]. For pRDA, the additional explanatory variables, called covariables, are assembled in matrix W; the linear effects of the explanatory variables in X on the response variables in Y are adjusted for the effects of the covariables in W [39]. In our study, the total variance of riverine NO 3 − and NH 4 + exports could be explained by the variables derived from human disturbance, climatic factors and landscape setting, and their individual contribution of NO 3 − and NH 4 + export can be finally figured out. We further partitioned the total variation of the riverine NO 3 − and NH 4 + response variables using three steps (Table 2). First, canonical ordination with no covariables was used to estimate the total amount of variance explained (as sum of canonical eigenvalues) in the NO 3 − and NH 4 + export attributed to all explanatory variables, human disturbances (H), climatic factors (C) and landscape setting (L), and the total unexplained variance (1-HCL). Second, the combinations of various covariables were considered to calculate the separate effect of each variable (H, C or L), in which an individual predictor variable was run (e.g., H) with the remaining other two as covariables (e.g., C&L). Third, a series of partial canonical ordinations were used to calculate the unique and interactive effects for each set of predictors (e.g., C&L-H) by considering the interaction term of interest as explanatory (C&L) and excluding the effect of not interest (e.g., H). For more details of calculations, please refer to related studies [19,29].

Scatterplot Matrix
The correlations between nutrient exports and variables at different buffer zones show that significantly higher coefficients are generally found between annual and seasonal nutrient exports and variables at watershed scale ( Figures S1-S3), whereas the fraction of buildup area within a 100 m buffer has higher correlations to DIN (r = 0.82, p < 0.01) in the dry season than other variables among different buffer zones ( Figure S3). Streamflow, the fraction of agricultural land cover of the entire watershed, buildup area within a 100 m buffer zone and buildup area of the entire watershed show significant positive relationships to DIN export (r = 0.52-0.64, p < 0.01; Figure 5). However, there is a negative relationship found between DIN and slope of entire watershed (r = −0.64, p < 0.01; Figure 5). Streamflow positively correlates to NO 3 − (r = 0.72, p < 0.01) and NH 4 + (r = 0.22, p < 0.01; Figure 5).

PCA of Environmental Variables
The first two principal components with eigenvalues > 5.0 are retained, accounting for 70.6% and 72.5% of variances on NO 3 − and NH 4 + exports, respectively ( Figure 6). The combinations of variables within the watershed scale (SF, SLP, AGR, and BD) and buildup area in a 100 m buffer (BD100) can explain much more variances of NO 3 − and NH 4 + exports than the performance using variables derived from the separate buffer zones (≤60%; Figure S4). Generally, the spatial variability of buildup (BD) and slope (SLP) are the main ingredients of the first PCs for NO 3 − and NH 4 + (Figure 6). The second PC is associated with seasonal variables, i.e., streamflow (SFWet and SFDry) and rainfall (RWet and RDry) during wet and dry periods. Most environmental variables displayed positive correlations with these two response variables (NO 3 − and NH 4 + ), except for landscape setting variables such as area (A), channel length (CL), longest channel length (LCL), relief (R) and slope (SLP), i.e., the opposite direction to NO 3 − and NH 4 + export (blue lines in Figure 6). During the wet season, the smaller projected angle is between the fraction of BD and NO 3 − export (NiWet), indicating high relevance between human impact and the wet season NO 3 − export, while during the dry season, NO 3 − export relates to the vectors dominated by streamflow and rainfall. However, unlike NO 3 − export, NH 4 + export is dominated by BD regardless of the different seasons. According to the results derived from PCA, five environmental variables including streamflow, slope, the fraction of agricultural land cover, buildup area within a 100 m buffer zone and buildup area of the entire watershed with higher loading were selected for further analysis.

Variance Partitioning-RDA and pRDA
The RDA model shows that the total variance explained (total inertia) by five selected environmental variables is always higher on NO 3 − exports than on NH 4 + exports. These predictor variables can explain 86% and 79% of the total variance of NO 3 − and NH 4 + exports, respectively, in the dry season, but only 73% and 69% in the wet season (Table 4). increase from 0.28 in the wet season, 0.52 in annual export, to 0.59 in the dry season. However, the effects of streamflow are low (λ 1 < 0.05) for NH 4 + annual and seasonal exports (Table 4). In contrast, slope has lower marginal effects on annual and seasonal NO 3 − exports (λ 1 = 0.12-0.32) than on NH 4 + exports (λ 1 = 0.56-0.64), and the marginal effects of slope are higher in the wet season (λ 1 = 0.32-0.64) than in the dry season (λ 1 = 0.12-0.61) for nutrient export ( Table 4)

Characteristics of DIN Concentrations and Exports in Taiwan
On average, DIN concentration from 43 island-wide watersheds across Taiwan is 1.66 mg-N L −1 , and the two main DIN species present a mediated dilution effect with streamflow (Table 3 and Figure  S5). The average riverine DIN export reaches 3100 kg-N km −2 yr −1 , which is much greater than the global mean (208 kg-N km −2 yr −1 ) [12]. High rainfall and streamflow, N deposition and N fertilizer application for agricultural production at upstream regions can account for the significant nutrients streamflow [11]. The results reveal that DIN export and concentration vary spatially, ranging from 200 kg-N km −2 yr −1 , 0.3 mg-N L −1 in less disturbed watersheds (site 39; Figure 3), to over 10,000 kg-N km −2 yr −1 , 8.8 mg-N L −1 in highly disturbed watersheds (site 16; Figure 3). The spatial DIN surge also indicates that the environmental background actually exports DIN and consequently induces the risk of eutrophication in the downstream. Previous studies suggested that intact forested watersheds demonstrated high N retention capacity but the capacity would collapse with significant land cover conversion [16,40]. In our study, higher NH 4 + exports in highly disturbed watersheds with smaller drainage area (<1000 km 2 ) and steep slope (>30%) show that these environmental backgrounds are unfavorable to ammonia oxidation or assimilation due to rapid transport [41]. In addition, most of the watersheds with higher NH 4 + exports than NO 3 − appear in plain areas where human disturbances are high and sewage systems are deficient [42].
The island-wide DIN estimation reveals that the wet season (May to October) contributes 78% of the annual DIN export, 3100 kg-N km −2 yr −1 , which is consistent with previous findings that hydrological processes control DIN export [43,44]. Because nutrient export and transport is regulated by precipitation and streamflow, this controlling factor is more dominant in the tropics and subtropics [12,45]. The torrential rainfall brought by typhoons during the humid growing summer in Taiwan accounts for 30−50% of the annual precipitation, 2500 mm yr −1 [13,46], and consequently causes vital effects on biogeochemical processes, i.e., a huge amount of nutrients being flushed out from terrestrial ecosystems to aquatic ecosystems. It also explains why the nutrient streamflow in Taiwan is in the leading place worldwide [12].

Climatic Control
In tropical/subtropical mountainous Taiwan, abundant rainfall usually leads to a great amount of net nutrient exports even during the growing summer, which is distinct from the findings in temperate forest ecosystems [11,45,47]. Water is the conveyor of ion movement. In spite of the dilution of NO 3 − concentration during flood periods (wet season), the extensive runoff by typhoons flushes over surface and near-surface and leads to greater DIN exports [13].
The relation between export and streamflow (F-Q relation), which definitely shows export change with streamflow, is particularly crucial as regarding nutrient balance and transport, although the F-Q relation, which inevitably incorporates streamflow in the calculation, could likely lead to the "spurious relation". However, the F-Q relation, which can present the dominance of supply-limited or kinetically-limited under different hydrologic conditions, helps to indicate the nutrient budget balance, and the transported amount is also important. Notably, there are totally different hydrological controls on NO 3 − and NH 4 + transport. One is that streamflow plays a strong role on NO 3 − export (r = 0.72; λ 1 = 0.28-0.59), which is in contrast to a relatively weak relation with NH 4 + export (r = 0.22; λ 1 < 0.05). The positive NO 3 − streamflow relation indicates the sources are relatively sufficient, as compared to kinetic transport. One possible interpretation is that warmer and more humid conditions (e.g., higher soil moisture) during the wet summer season are favorable for nitrification and promotes NO 3 − accumulation and then transports it from the soil to the aquatic system [48,49]. On the contrary, NH 4 + is easily converted and emitted to the atmosphere via microbial activities in warm and humid conditions [15,50], which has lower retention in soil capacity compared to NO 3 − . Whether NH 4 + is source-limited casts a shadow of doubt on hydrologic control on NH 4 + transport.

The Consideration of Landscape and Buffer Zone
Streams receive nutrients from aquatic ecosystems and adjacent terrestrial landscapes, such that landscape configuration at watershed-or buffer zone-scale has important influences on stream water quality, ecological process and biodiversity. Watersheds with steeper slopes usually export more nutrients [51], yet our results surprisingly show that slope is negatively related to NO 3 − and NH 4 + exports ( Figure 5). The possible reason is that most upstream watersheds with steep slopes are covered by pristine forests where anthropogenic sources are scarce [10,52]. Studies suggested that elevation might be a suitable parameter to predict water quality [30], but the collinearity between landscape (e.g., elevation or SLP) and human-made land cover (e.g., AGR and BD) from the high, steep montane region to the low, flat plain in Taiwan keeps us from using both simultaneously, even they might have significant contribution on nutrient exports. Many studies demonstrate the control of land use and landscape on water quality at watershed scale [28,30], but few studies consider the distance from source area to river, i.e., the regulation of the riparian or buffer zone on nutrient movement. Plants and microbial activities within the riparian zone can help to uptake a great quantity of water, nutrients and sediments, mitigating the nutrients export to aquatic ecosystems within the watershed [53,54]. However, the effective buffer distance is uncertain, and it could vary with the elements which are concerned [36,55]. A study conducted in Puerto Rico examined how the landscape pattern changes affected water quality in-stream and found that turbidity and dissolved oxygen responded to land use and land cover (LUC) at watershed scale [35], phosphorus concentration and fecal matter responded to LUC at sub-watershed scale, whereas nitrogen concentrations linked to LUC in riparian buffers of larger watersheds [35]. Another study conducted in the Saginaw Bay of central Michigan investigating 62 catchments with a gradient of disturbed land cover showed that the land use factors within a 100 m buffer zone adjacent to the river could explain much higher (or equal) variance of NO 3 − and NH 4 + concentrations than those derived from the watershed scale [52]. Conversely, the relationships between total nitrogen exports and land use at the watershed scale were better than riparian buffer zones of 200 and 400 m in highly disturbed rivers in Illinois and Texas [56]. Previous assessments, based upon mountainous background, demonstrated that the buffer zone within a 100 m riparian zone, the buildup area particularly, plays an important role in regulating DIN exports [56][57][58]. In Taiwan, the deficiency of sewage systems in buildup areas located near river networks contributes a significant amount to the DIN exports [42]. The percentage of buildup at the watershed scale highlighted in our study indicates that dispersal non-point sewage sources in the buildup area and scattered agricultural activities would be critical in assessing NH 4 + export.

Human Disturbance
The contribution of human activities to DIN exports (including NO 3 − and NH 4 + ) from land to water has been underscored due to agricultural activities (e.g., fertilizer application) and urbanization (e.g., domestic wastewater) [52,59]. To meet the ever-increasing demand, agricultural activities, such as high value fruit, montane cabbage and tea plantations, are pervasive in mountainous Taiwan. Consequently, the excessive addition of inorganic N fertilizer and organic manure on the thin soil layer is readily flushed out to aquatic ecosystems downstream during wet season [12,60]. The N retention or removal capacity within watersheds will dramatically decrease if forests are transformed to agricultural areas [61]. However, though the effects of agricultural land on NO 3 − and NH 4 + exports are statistically significant, their contributions are weaker than previous studies suggested [57,58]. A Canadian study suggested that urban land use has a stronger effect on water quality than agriculture has [14]. Obviously, the importance of agricultural and urban land use on DIN exports is not easy to identify and separate, not only due to the area, but also the "intensity" (e.g., intensive agriculture or dense population) and spatial configuration. For example, the effect of urban areas is minor in Finland, because the population and settlements are mostly scattered so that wastewater can be purified before it flows into main streams [28]. In this study, both NO 3 − and NH 4 + exports are significantly correlated to AGR and BD, which shows that AGR and BD are important for DIN exports ( Figure 5). Notably, higher correlation coefficients of AGR and BD to NH 4 + export (r = 0.67-0.75) than those to NO 3 − export (r = 0. 32-0.40) indicates that land use pattern is more effective to explain the variance of NH 4 + export. Conceptually, NH 4 + is easier to be taken up by plants and to be oxidized through nitrification and thus NO 3 − is the main species of DIN within agriculture-dominated catchments. The intensive urban developments located near stream riparian in Taiwan deteriorates water quality directly and contributes a significant amount to the DIN exports. This phenomenon is reflected by the higher effect of BD100 (buildup area in a 100 m buffer zone) than the BD of the entire watershed on NO 3 − export during the dry season (Table 4). Therefore, domestic wastewater inevitably elevates NH 4 + export, which indicates a strong intrinsic collinearity with agriculture and buildup area. Such intrinsic collinearity presents an inseparable human-landscape system, and their interplay could not be distinguished perfectly.

Interactive Effects among Variables
It was noticeable that interactive effects between landscape and human variables can explain most variabilities of NO 3 − export in the wet season (51.59% of the total variance) and seasonal NH 4 + export (both >80% of the total variance). Meanwhile, interactive effects between climatic and human variables explain 44.15% of the total NO 3 − export variance in the dry season. Thus, the interactive effect between landscape setting and human disturbance and climatic and human disturbance will result in a high efficiency of prediction regarding NO 3 − and NH 4 + export in different conditions. One explanatory variable might be partially linked to other variables and altogether would improve or reduce the predictive power in evaluation of NO 3 − and NH 4 + export [62]. For example, the opposite influences of streamflow and slope might suppress their capability to assess NO 3 − and NH 4 + exports without eliminating the collinearity between them. The mixed effects of environmental variables also suggest that it is necessary and will be more effective to apply an integrative management strategy [63].
There is a significant change between the controlling interactive effects of NO 3 − export in the wet and dry seasons. In the wet season, streamflow and interactive effects between landscape setting and human disturbance might reveal that they are two unique groups controlling the variance of NO 3 − export. This might be reasonable because the abundant rainfall during the humid summer will decrease the spatial variability of streamflow across Taiwan, which leads to a weaker relationship between streamflow and the interactive effects of landscape setting and human disturbance. On the contrary, the increase of the interactive effects between climatic and human disturbance might indicate the increasing spatial variability of streamflow in the dry season. Therefore, when we predict NO 3 − export in the wet season, both streamflow and the interactive effects between landscape setting and human disturbance cannot be ignored. Furthermore, the strong interactive effects between landscape setting and human disturbance on NH 4 + export, indicates that these two groups are the primary control of NH 4 + export. Because landscape setting and human disturbance are highly related ( Figure 5), the variance explained by them might be similar and it might be difficult to separate their individual effect from both combined. Some critical variables are not included in this study but might be significant in other regions. Studies found a strong negative relationship between nitrogen content and soil type [14]. However, most watersheds in Taiwan cross various soil substrates and geological units, which challenges us to unambiguously define a specific soil type for each watershed. In addition, long-lasting atmospheric nitrogen deposition could lead to N saturation in temperate forest ecosystems, resulting in net loss of nitrate to streams and consequent acidification of stream water [1,16,64]. A large-scale study indicated that riverine N export could be predicted by atmospheric N deposition rates [65]. Therefore, the influence of interaction between atmospheric deposition, land use and hydroclimate should be considered in following syntheses of DIN responses or developing models for riverine DIN export [66]. for the wet and dry seasons, respectively, and climatic variables (e.g., streamflow) are not effective variables for NH 4 + transport. The human-landscape variables are the major factors to explain the total variance of NH 4 + export (over 80%), in both the wet and dry seasons. The shift of interactive effects of variables on nutrient export is important for water quality management at watershed scale and designing mitigation strategies. Inevitably, the effects of intrinsic collinearity in the human-landscape system cannot be clearly separated due to spurious correlation, though the statistical approach provides some cues. For example, paired AGR and BD or SLP and BD are highly collinear but difficult to single out for estimating nutrient export and for interpretation. Nevertheless, with the accumulation of these studies, it is more possible to clarify the interactive effects, which could be of great help in advancing the understanding of DIN export mechanisms and global synthesized assessment. yr −1 ). Figure S1: