Next Article in Journal
Water, Energy, and Emissions Nexus: Effect of Inflows in Urban Drainage Systems
Previous Article in Journal
Spatio-Temporal Characteristics of Drought and Its Relationship with El Niño-Southern Oscillation in the Songhua River Basin from 1960 to 2019
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Understanding the Combined Effects of Land Cover, Precipitation and Catchment Size on Nitrogen and Discharge—A Case Study of the Mississippi River Basin

Faculty of Geography, Philipps-University of Marburg, Deutschhausstr. 10, 35037 Marburg, Germany
*
Author to whom correspondence should be addressed.
Water 2022, 14(6), 865; https://doi.org/10.3390/w14060865
Submission received: 26 January 2022 / Revised: 27 February 2022 / Accepted: 8 March 2022 / Published: 10 March 2022
(This article belongs to the Section Water Quality and Contamination)

Abstract

:
Biological processes of rivers are strongly influenced by concentration and fluxes of nitrogen (N) levels. In order to restrain eutrophication, which is typically caused by urbanisation and agricultural expansion, nitrogen levels must be carefully controlled. Data from 2013 to 2017 were gathered from 26 sub-catchments in the Mississippi River basin to assess the effects that catchment size, land cover, and precipitation can have on the discharge and total nitrogen (TN) and how TN yields deviate from a generalised local trend. The findings indicated that land cover and precipitation had a determinative effect on area-weighted discharge (Qarea). More specifically, Qarea had significant positive (directly proportional) relationships with precipitation, forest, and urbanised land cover, and significant negative (inversely proportional) relationships with grassland/pasture and scrub/shrub land covers. Concurrently, the TN concentration significantly increased in the presence of agricultural land cover, but significantly decreased in forest land cover. The TN yield (TN concentration × Qarea) was largely determined by Qarea because the latter was observed to fluctuate more dramatically than concentration levels. Consequently, the TN yield exhibited the same relationships that Qarea had with precipitation and land covers. The TN yield changed significantly (p < 0.05) and positively with instantaneous discharge across all sites. Nevertheless, the rate of TN yield variations with discharge displayed a significant (p < 0.0001) negative (r2 = 0.80) relation with the catchment size. Ultimately, this study used discharge readings to facilitate the prediction of TN concentrations and yields across various catchment areas in the Mississippi River basin and provided a robust model for future research in this area.

Graphical Abstract

1. Introduction

Human nutrition has considerably improved owing to the use of nitrogen-based fertilisers in agriculture [1,2]. Yet, the overuse of these fertilisers has resulted in unintended health-related and environmental issues [3,4]. For example, the overuse of fertilisers has engendered eutrophication, which can cause more frequent algal blooms, oxygen depletion, water turbidity, and biodiversity loss [5,6]. These consequences have progressively negative effects, for example, oxygen depletion can disband heavy metals that were once precipitated or bound to river sediment particles [7]. These metals are generally associated with harmful effects on water quality and ecosystem wellness [8,9].
Since the Industrial Revolution, nitrogen levels have progressively increased across the globe as a result of human activities [10]. For example, the annual use of commercial nitrogen-based fertilisers increased from 594,000 to 11.5 million tons between 1945 and 1985 in the United States of America [11]. Simultaneously, the concentration of nitrogen in America’s water has increased to critical levels [12,13]. Global mobilisable nitrogen levels (Nmob) have increased dramatically since industrialisation. Prior to industrialisation, the level was 111 Tg Nmob/year, and at present it is 223 Tg Nmob/year (Table 1). The pre-industrial nitrogen levels are consistent with estimates made in existing research by [14,15] and cohere with those generated by [16] despite being slightly lower. As a result of increased anthropogenic emissions, such as those produced through industry, the cultivation of crops, and the use of industrial fertilisers, the present Nmob levels are high [17].
The Mississippi River basin is the third biggest basin in the world, covering approximately 41% of the contiguous United States [19], and is the major source of freshwater and nutrients to the Gulf of Mexico [20]. Ref. [21] approximated the amount of nitrogen delivered to the Gulf of Mexico from the Mississippi River basin. They identified six sources of nitrogen, viz., chemical fertilisers, atmospheric deposition, manure, fixation and other legume sources, urban effluent, and wastewater treatment facilities, which account for 41, 26, 10, 9, 7, and 7%, respectively. The Mississippi River has undergone significant engineering changes throughout the years as a result of human activities, such as urbanization and agricultural growth. These activities have contributed remarkably to the deterioration of water quality in the northern Gulf of Mexico and to the occurrence of seasonal hypoxia (dissolved oxygen less than 2 mg L−1) each summer since the 1980s [22]. This is the world’s second biggest human-induced hypoxic zone, ranging in size from 40 km2 in 1988 owing to drought to 22,730 km2 in 2017 [23,24]. The Mississippi River/Gulf of Mexico Watershed Nutrient Task Force established a target of lowering the hypoxic zone’s size to 5000 km2 by the year 2035 [25]. However, long-term research has shown that hypoxia has tended to extend over time, posing negative ecological and economic impacts [26,27].
Land cover change is a significant factor influencing N export from the land. The Mississippi River basin is one of the most productive agricultural areas in the United States and has been notably influenced by land cover changes [28,29]. Conversion of natural vegetation to croplands has resulted in a large increase in nitrogen loading beyond baseline levels in the Mississippi River basin prior to the heavy use of synthetic nitrogen fertilisers [30]. On the other hand, climatic change is predicted to also alter nitrogen loadings. Various climate studies predicted certain changes in future climatological conditions, especially for the primary hydrological controller, precipitation [31]. If precipitation alters, the discharge and accumulation of nitrogen will change considerably [32,33,34]. Existing research contends that the accumulation and discharge of nitrogen will increase in tandem with increased precipitation [35,36,37,38,39,40,41]. Similarly, [42] found that the degree of interannual variability in nitrogen loading from the Mississippi basin may vary by a factor of 2.3, with 76% of the variability related to annual precipitation fluctuations. However, the degree to which precipitation fluctuation has an effect on N yield and transport to the Gulf remains unknown. This unpredictability makes identifying and implementing effective N reduction techniques in the case of catastrophic occurrences more difficult [43]. Variation in freshwater nitrogen cycles depends on multiple factors, rather than a single human activity or climate change variable [44]. This is well demonstrated by [45], who observed that N concentrations decrease with forest area ratios and decrease further as levels of precipitation increase. On the other hand, they found that the slope of the concentration and agricultural lands was positive and increased further with precipitation. Likewise [46], concluded that the nitrogen flux from agricultural land is high and increases during periods of heavy rainfall. Consequently, future N fluctuations will result from variations in precipitation and land cover [31]. However, it remains impossible to state whether a trend exists between nitrogen concentrations and yields, with respect to precipitation and/or land cover variations. Existing research typically considers the impacts of land cover and precipitation individually. Limited work has been done to investigate the combined effects of precipitation and land cover on the nitrogen delivery within freshwater, and this research often analyses a particular spatial scale e.g., [45,47,48]. Therefore, further research is needed to understand the behaviour of nitrogen in relation to spatial scale.
Parsons et al. [49,50] detected a relationship between catchment size and erosion rates and hypothesised that a similar relationship may be found between the behaviour of nitrogen and catchment size. Therefore, this hypothesis was investigated in the current study, i.e., catchment area that controls erosion rate could also control nitrogen yield in the Mississippi basin. The present study investigated the combined effects of precipitation and land cover on discharge and nitrogen concentrations within variously sized sub-catchments of the Mississippi basin. Such investigation is important to characterize the relative importance of these drivers (precipitation and land cover) on future changes of the hydrological regime. From the viewpoint of stakeholders or water resources managers, the findings of this research may help watershed stewardship programs within the context of climate and land cover change adaptation. For instance, findings from such studies could be used to help sustain the protection of natural habitat (e.g., wooded regions and wetlands), define long-term impacts on ecological objectives, such as algal blooms, in surface waterways, and identify long-term land cover zoning plans that protect water quality.

2. Materials and Methods

The TN budgets of 26 sub-catchments (hereafter, catchments) in the Mississippi River basin were quantified (Figure 1). Using the United States Geological Survey [51], data pertaining to the monthly nitrogen concentrations and discharges within each catchment between 2013 and 2017 were gathered. The National Water Quality Network (NWQN) contains 110 river sites across the United States monitored by the USGS National Water Quality Program. This network includes monthly TN concentrations and discharge data. The selection of monthly scales was contingent on data availability. The National Land Cover Database was used to collect data about the land cover with a 30 m resolution. Precipitation data for the stations dispersed over each catchment were gathered from the National Oceanic and Atmospheric Administration [52] from 2013 to 2017. The Thiessen polygon technique (Equations (1) and (2)) was used to interpolate the precipitation map for each catchment. In this technique, the precipitation reported at each station is weighted according to the region nearest to the station. The following approach is used to determine the weighing area: The precipitation stations are connected to create a triangular network. To form a polygon around each station, perpendicular bisectors for each of the triangle’s sides are created. Thiessen polygons is the name given to these bounding polygons. If P1, P2,…, and PM denote the precipitation magnitudes measured at stations 1, 2,…, and M; and A1, A2,… AM denote the corresponding areas of the Thiessen polygons, then the average precipitation across the catchment P for a catchment area A is given by:
P ¯   = P 1 A 1 + P 2 A 2 + + P M A M A 1 +   A 2 + + A M
Thus in general for M stations:
P ¯   = i = 1 M P i A i A = i = 1 M P i A i A
The ratio Ai A is referred to as the station’s weighting factor. The Thiessen polygon approach is preferable to the arithmetic-average method for computing the average precipitation across an area because weighting is applied to the individual stations on a rational basis. Once the weighting parameters are defined, calculating the average P for a set network of stations is quite straightforward [53]. To construct the Thiessen polygons in ArcGIS, we selected Arc Toolbox > Analysis Tools > Proximity > Create Thiessen Polygons. The weighting factor assigned to each station was multiplied by the station’s P. The sum of the weightage factor and p values for all stations was the catchment’s average p value. Catchments were selected on the basis of their areas, precipitation levels, and the availability of data relating to the TN concentrations and nitrogen discharge. The catchments under study are geographically widespread, demonstrating diverse levels of areas, land cover, and precipitation (Figure 1). The catchment areas range from approximately 3210 km2 to 1,848,000 km2 (Table 2).
The various forms of land cover included pasture/grassland, forest, agricultural land, wetland/water, urbanised land, scrub/shrub lands, and barren lands.
In order to limit the effects of anomalous results among the TN concentration readings, area-weighted discharge values (Qarea), and TN yields between catchments, the median measurements of each variable were employed in the regression analysis (Table 3). This meant that the median TN concentration, median Qarea, and median TN yield were assessed to determine whether a relationship existed between these variables, and the median monthly precipitation, catchment size, and land cover under study. The TN levels were determined using the following formula:
TN   yield   =   TN ×   Q / A
TN yield describes the total nitrogen yields at the catchment outlet (mg TN m−2 min−1); [TN] refers to the total nitrogen concentrations (mg TN L−1); Q indicates the level of discharge at the catchment outlet (L min−1); and A refers to the catchment area (m2).
The relationship between the TN yield and instantaneous discharge were quantified at each catchment. In order to understand the effects of nitrogen discharge on the TN yield, the latter was plotted against discharge for each catchment. A positive correlation between the TN yield and the levels of discharge was expected; however, this was initially tested in accordance with the relationship between nitrogen concentration and discharge. The presence of chemostatic, flushing, or dilution behaviour (i.e., instances in which chemical concentration remains constant, increases, or decreases with discharge respectively) is often assessed by examining the relationship between element concentration and discharge [54]. The TN yield had positive relations with both nitrogen concentrations and discharge levels according to Equation (3). Therefore, if the TN concentrations did not strongly decline with discharge (dilution), constant or ascending concentrations with discharges (chemostatic or flushing) resulted in a positive relationship between yield and discharge. Subsequently, the slope of the relationship between the TN yield and discharge (hereafter, TN yield rate) was quantified for each catchment. The TN yield rate for each catchment was then projected versus its area. A statistical analysis was undertaken using JMP14.2 (SAS Institute, Cary, NC, USA). Within our analysis, the significance level was 0.05. If the p value was found to be less than or equivalent to 0.05, the result was determined to be statistically significant.

3. Results

The analysis found that no significant trend existed between catchment area and precipitation (Table 4), and this indicated that our results were not influenced by the catchment area or climate conditions. The collinearity (r2) between precipitation and types of land cover reached 78%. Applying variance inflation factors (VIF) to investigate the multicollinearity among catchment area, precipitation, and land cover types showed a significant multicollinearity. Specifically, agriculture land cover exhibited high VIFs (more than 10 *) indicating that this land cover category was highly correlated with at least one of the other parameters in the regression. Excluding this land cover type from the analysis to reduce multicollinearity resulted in low VIFs (less than 10). The ratio of forest and urbanised land cover displayed positive significant (p < 0.05) correlations with precipitation, and the percentages of barren, shrub/scrub land, and grassland/pasture had negative correlations with the levels of precipitation (Table 4). Our study found that barren land cover exhibited a significant positive correlation with scrub/shrub land cover. Results also showed a significant positive correlation between forest and urbanised land covers, a significant negative correlation between wetland/water and grassland/pasture cover, a significant positive correlation between grassland/pasture and scrub/shrub land cover, a significant negative correlation between grassland/pasture and urbanised land covers, and a significant negative correlation between scrub/shrub and urbanised land covers (Table 4).

3.1. Total Nitrogen, Precipitation, Land Cover, and Catchment Size

This analysis found that the percentage of agricultural land had a significant positive relationship with the median TN concentration (p < 0.05), while the percentage of forest land cover had a negative correlation with the median TN concentration (Table 5).
*A rule of thumb is that if the VIF is more than 10 then multicollinearity exists [55]. In response to precipitation, and ratios of forest and urbanised land cover, the median Qarea exhibited a significant positively correlated variation (Table 5). Conversely, the median Qarea varied in negative relations with the ratio of grassland/pasture and scrub/shrub lands (Table 5). In addition, a positive correlation existed between the median TN yield and precipitation, agricultural and urbanised land covers, while a negative correlation existed between the median TN yield and grassland/pasture and scrub/shrub land covers (Table 5). The slopes of regression between the TN yield and precipitation and land cover categories were analysed. The slopes of such regression between yield and precipitation, urban, pasture/grassland, shrub/scrub, and agriculture land covers were 0.002324, 0.000235, −0.000063, −0.000059, and 0.000041, respectively. Such regression indicated that the TN yield was more sensitive to precipitation variations followed by urban, pasture/grassland, shrub/scrub, and agriculture land cover variations, respectively. Across the investigated catchments, the relationship between nitrogen concentration and discharge showed a positive correlation in 23 catchments; a negative correlation in 2 catchments, and an insignificant correlation in 1 catchment (Table 2; Figure 2). Finally, across every catchment, a strong positive relation was observed between the TN yield and the discharge (Table 2; Figure 2).

3.2. TN Yield Rates

Overall, the TN yields increased with the levels of discharge, though the rate of such increase became smaller with the catchment size. The TN yield rate versus catchment area had a significant negative correlation (p < 0.0001; r2 = 0.80) (Figure 3). This can be seen in the following equation:
Log   Y = 12.37 1.02 × Log   X
where Log indicates the logarithm based 10, Y represents the TN yield rate and X refers to the catchment area.

4. Discussion

Precipitation is documented to be a controlling parameter for the growth and distribution of vegetation e.g., [56,57,58]. A substantial positive relation was reported between precipitation levels and vegetation density [59]. Likewise, [60] studied the effects of precipitation levels on land cover variation. According to their research, precipitation displayed a gradual increase from barren lands to grasslands to agricultural and forest lands. Simultaneously, land cover itself can affect precipitation levels. Ref. [61] found that highly urbanised environments typically experienced more frequent heavy precipitation. In their studies of Amazonia [62,63], found that precipitation levels declined as a result of deforestation, depending on the variant of land cover that substituted the forested area. They found that replacing forested areas with agricultural land resulted in greater precipitation than replacement with pasture/grassland. Contrastingly [64], found that the conversion of afforestation in the United States increased precipitation levels. Other research found that precipitation significantly enhances with vegetation density [65,66]. This is due to the fact that greater vegetation density quickens the rate of evapotranspiration, and thereby results in a lower vapour pressure deficit, which, in turn, increases the formation of clouds and rainfall [67,68]. The abovementioned explanation can justify the existence of the positive relations between urbanised and forested lands and precipitation levels, and the negative relations between bare lands, grassland/pasture, and scrub/shrub lands with precipitation (Table 4). Qarea increased with the level of precipitation (r = 0.87; Table 5), and this could be attributed to the enhanced water availability for runoff [69,70]. Since precipitation drives Qarea, the latter followed the same correlations between the former land cover types (i.e., positive relations with urban and forest land cover, and negative relations with grassland/pasture and scrub/shrub land cover) (Table 4 and Table 5).

4.1. Total Nitrogen, Precipitation, Land Cover, and Area

Across the investigated catchments, the TN concentration was found to increase with the ratio of crop land cover (Table 5). The use of agricultural fertiliser generates higher levels of nitrogen in the soil and water runoff [71,72]. Moreover, tillage and other agricultural activities can cause erosion, which also increases the loss of nitrogen. Many studies have found that tillage causes markedly high levels of soil erosion and nitrogen loss [73,74,75]. On the other hand, TN concentrations are quite low in forested areas due to the limited input of nitrogen, ongoing microbial activities, and dense vegetation within that environment (Table 5). The accumulation of atmospheric nitrogen is the most significant cause of nitrogen build-up in watersheds within areas of high natural vegetation (such as forests). Nevertheless, these levels remain significantly lower than that caused by human activity (such as the use of fertilisers and the production of wastewater) within agricultural and urbanised areas [76]. Without the influence of anthropogenic inputs, the transfer of nitrogen from upland forests into streams can be eliminated by soil-stream interfaces [77,78,79], microbial absorption assimilation, or the denitrification process [80,81,82]. In comparison with cropland or grassland, forested land has the lowest levels of microbial nitrogen fixation [83]. In turn, dense vegetation within forested areas limits the incidence of soil erosion and runoff because soils are more stable and the vegetation can absorb precipitation. Therefore, nitrogen loss is limited [84]. Surprisingly, the results of this study did not show a correlation between the ratio of urbanised land cover and the TN concentration (Table 5). Within a previous study that involved mixed land covers such as the watersheds of the Menominee River, Altamaha River, Connecticut River, and Upper Snake River, low TN concentrations were found. Within these large watersheds, the runoff from urbanised and agricultural environments can become diluted in forested and other relatively undeveloped lands [85]. The results of the current analysis cohere with this study, as the nitrogen which had originated from urbanised areas was likely to have been diluted upon reaching undeveloped lands.
A positive correlation was found between TN yields and discharge, which can be explained based on the positive relations between TN concentrations and discharge. The presence of positive relations between the TN concentration and discharge level across the majority of catchments (Table 2) can rationalise the positive relations between TN yields and discharge levels (Table 2; Figure 2). Significantly, within the two catchments that showed a negative correlation between the TN concentrations and discharge, there were higher levels of urbanised land than in other catchments. More specifically, these two catchments had a significant amount of high and medium intensity urbanised land compared to other catchments. This finding was affirmed with reference to data from SPARROW (Spatially Referenced Regression on Watershed attributes) [86], which shows that nitrogen typically emerges as a result of wastewater treatment processes, urbanised land, and the use of agricultural fertiliser and manure. With regard to the two catchments mentioned previously, SPARROW data showed that these catchments had a higher level of nitrogen discharge from wastewater treatment than other catchments. Within these catchments, the presence of a negative correlation between TN concentration and discharge meant a high TN concentration at a low discharge level which could be ascribed to a constant source of nitrogen that is diluted with high discharges (i.e., urban wastewater in the current study). In spite of the negative correlation between the TN concentration and the Qarea in these two catchments, a positive correlation still existed between the TN yield and the Qarea. This was the result of the rate of variation between the TN concentration levels and the Qarea. For example, the South Platte River (near Kersey, Colorado) showed the highest negative correlation between the TN concentration and the Qarea. However, within this river, the TN yield still increased in tandem with the levels of Qarea. These trends are explained by the fact that the TN concentration declined by a factor of 4.3, whereas the Qarea levels increased by a factor of 39.8. Therefore, the results of the present study cohered with the work of [87], which illustrated that streamflow affects nutrient yield variation more than concentration levels. Therefore, this study showed that the N yield follows Qarea, in that it shares a positive relation with precipitation and urbanised land and has a negative correlation with grassland/pasture and scrub/shrub land cover (Table 5). Notably, no correlation was found between the TN yield and forested land, even though forested cover maintains reasonably low TN concentration (Table 5). In fact, the TN yield remained consistent within forested areas, as the low TN concentrations were counterbalanced by a heightened Qarea (Table 5).

4.2. TN Yield Rate

As mentioned previously, the TN yield increased in tandem with the levels of discharge. However, this rate of increase was contingent on the size of the catchment. Previous research has asserted that the erosion of soil by water contributes significant amounts of nitrogen to various ecosystems [88,89,90,91]. In catchments greater than 10 km2, the level of erosion and, consequently, yield was documented to decrease with catchment size [92,93]. Within these catchments (more than 10 km2), the sediment sinking potential is typically higher than the sediment sourcing potential, leading to a reduction in sediment yield [94]. This can be ascribed to the fact that large catchments encompass greater floodplain development and more foot slope terrains, in which sediment can be stored [92]. It is also more likely that sediment will be deposited before a catchment’s outlet point, as it must travel a greater distance [95,96]. Many research works have found a negative relation between sediment yields and catchment size [97,98,99], which explains the negative relation between the TN yields and catchment size within this study (Figure 3). This means that smaller catchments experience higher erosion levels and as erosion is a primary controller for nitrogen sourcing, such catchments have higher nitrogen yield rates than larger catchments [100].
Within the present study, a geographically diverse array of catchment areas, various forms of land cover, and differing levels of precipitation were analysed. This study illustrated that the findings presented within Figure 3 can be utilised as a means to predict the TN yield rate throughout the Mississippi basin. By investigating the TN yield rate within a certain catchment area, the TN yields and concentrations levels can be deduced in each instance of discharge.

5. Conclusions

This study found that the TN yield (TN concentration × Qarea) was largely determined by Qarea because the latter fluctuated more dramatically than nitrogen concentration levels. In addition, the TN concentration and Qarea differed according to the level of precipitation and land cover conditions. The TN concentration was found to increase in positive correlations with discharge across most catchments under study except in four catchments. In those four catchments, urban effluent was an important N source, which represented a steady source that was diluted at higher discharges. Variations in precipitation and/or land cover could affect discharge and/or TN concentration and consequently the TN yield. The slopes of regression between the TN yield and precipitation and land cover categories indicated that the N yield was more sensitive to precipitation variations followed by urban, pasture/grassland, shrub/scrub, and agriculture land cover variations, respectively. The TN yields increased in tandem with the levels of discharge, even though the rate of such an increase declined with catchment size. Overall, this study involved the analysis of broad spatial scales, land cover, and precipitation, across several catchments within the Mississippi River basin. Ultimately, the study’s findings supported the use of discharge measurements and catchment size as a means to predict TN concentrations and yields.

Author Contributions

Conceptualization, data curation, formal analysis, methodology, software, and writing—original draft by H.A. and C.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. The APC was paid by The University of Marburg.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stewart, W.M.; Roberts, T.L. Food Security and the Role of Fertilizer in Supporting it. Procedia Eng. 2012, 46, 76–82. [Google Scholar] [CrossRef] [Green Version]
  2. Davidson, E.A.; Nifong, R.L.; Ferguson, R.B.; Palm, C.; Osmond, D.L.; Baron, J.S. Nutrients in the nexus. J. Environ. Stud. Sci. 2016, 6, 25–38. [Google Scholar] [CrossRef] [Green Version]
  3. Caccia, V.G.; Boyer, J.N. Spatial patterning of water quality in Biscayne Bay Florida as a function of land use and water management. Mar. Pollut. Bull. 2005, 50, 1416–1429. [Google Scholar] [CrossRef] [PubMed]
  4. Howarth, R.; Swaney, D.; Billen, G.; Garnier, J.; Hong, B.; Humborg, C.; Johnes, P.; Mörth, C.-M.; Marino, R. Nitrogen fluxes from the landscape are controlled by net anthropogenic nitrogen inputs and by climate. Front. Ecol. Environ. 2012, 10, 37–43. [Google Scholar] [CrossRef] [Green Version]
  5. World Health Organisation (WHO); European Commission. Eutrophication and Health; WHO (Regional Office for Europe) and the European Commission: Luxembourg, 2002; Available online: https://ec.europa.eu/environment/water/water-nitrates/pdf/eutrophication.pdf (accessed on 17 February 2022).
  6. Russell, M.J.; Weller, D.E.; Jordan, T.E.; Sigwart, K.J.; Sullivan, K.J. Net anthropogenic phosphorus inputs: Spatial and temporal variability in the Chesapeake Bay region. Biogeochemistry 2008, 88, 285–304. [Google Scholar] [CrossRef]
  7. Antweiler, R.C.; Goolsby, D.A.; Taylor, H.E. Nutrients in the Mississippi River. In Contaminants in the Mississippi River 1987–1992; Meade, R.H., Ed.; U.S. Geological Survey Circular: Denver, CO, USA, 1995; Volume 1133, pp. 73–86. [Google Scholar]
  8. Jonsson, K.; Johansson, H.; Worman, A. Hyporheic exchange of reactive and conservative solutes in streams tracer methodology and model interpretation. J. Hydrol. 2003, 278, 153–171. [Google Scholar] [CrossRef]
  9. Ip, C.C.; Li, X.-D.; Zhang, G.; Wai, O.W.; Li, Y.-S. Trace metal distribution in sediments of the Pearl River Estuary and the surrounding coastal area, South China. Environ. Pollut. 2007, 147, 311–323. [Google Scholar] [CrossRef] [Green Version]
  10. Schlesinger, W.H.; Bernhardt, E.S. Biogeochemistry: An Analysis of Global Change, 3rd ed.; Academic Press: Waltham, MA, USA, 2013. [Google Scholar]
  11. Alexander, R.B.; Smith, R.A. County-Level Estimates of Nitrogen and Phosphorus Fertilizer Use in the United States, 1945 to 1985; Open-File Report 90-130; U.S. Geological Survey: Reston, VA, USA, 1990; pp. 15.
  12. Howarth, R.; Sharpley, A.; Walker, D. Sources of Nutrient Pollution to Coastal Waters in the United States: Implications for Achieving Coastal Water Quality Goals. Estuaries 2002, 25, 656–676. [Google Scholar] [CrossRef]
  13. Turner, R.E.; Rabalais, N.N.; Justic, D. Gulf of Mexico hypoxia: Altered states and a legacy. Environ. Sci. Technol. 2008, 42, 2323–2327. [Google Scholar] [CrossRef]
  14. Dentener, F.J.; Crutzen, P.J. A three dimensional model of the global ammonia cycle. J. Atmos. Chem. 1994, 19, 331–369. [Google Scholar] [CrossRef]
  15. Smil, V. Nitrogen in crop production: An account of global flows. Glob. Biogeochem. Cycles 1999, 13, 647–662. [Google Scholar] [CrossRef] [Green Version]
  16. Galloway, J.N.; Dentener, F.J.; Capone, D.G.; Boyer, E.W.; Howarth, R.W.; Seitzinger, S.P.; Asner, G.P.; Cleveland, C.C.; Green, P.A.; Holland, E.A.; et al. Nitrogen cycles: Past, present, and future. Biogeochemistry 2004, 70, 153–226. [Google Scholar] [CrossRef]
  17. Green, P.A.; Vörösmarty, C.J.; Meybeck, M.; Galloway, J.; Peterson, B.J.; Boyer, E. Pre-industrial and contemporary fluxes of nitrogen through rivers: A global assessment based on typology. Biogeochemistry 2004, 68, 71–105. [Google Scholar] [CrossRef]
  18. Worldometers.info. 7 Continents. 2021. Available online: https://www.worldometers.info/geography/7-continents/ (accessed on 17 February 2022).
  19. Battaglin, W.A.; Aulenbach, B.T.; Aldo, V.; Buxton, H.T. Changes in Streamflow and the Flux of Nutrients in the Mississippi-Atchafalaya River Basin, USA, 1980–2007; U.S. Geological Survey: Reston, VA, USA, 2010.
  20. Dale, V.H.; Armitage, T.; Bianchi, T.; Blumberg, A.; Boynton, W.; Conley, D.J.; Crumpton, W.; David, M.; Gilbert, D.; Howarth, R.W.; et al. Hypoxia in the Northern Gulf of Mexico; Springer: New York, NY, USA, 2010; p. 284. [Google Scholar]
  21. Robertson, D.M.; Saad, D.A. SPARROW Models Used to Understand Nutrient Sources in the Mississippi/Atchafalaya River Basin. J. Environ. Qual. 2014, 42, 1422. [Google Scholar] [CrossRef] [PubMed]
  22. Tian, H.; Xu, R.; Pan, S.; Yao, Y.; Bian, Z.; Cai, W.; Hopkinson, C.S.; Justic, D.; Lohrenz, S.; Lu, C.; et al. Long-Term Trajectory of Nitrogen Loading and Delivery From Mississippi River Basin to the Gulf of Mexico. Glob. Biogeochem. Cycles 2020, 34, e2019GB006475. [Google Scholar] [CrossRef]
  23. Scavia, D.; Rabalais, N.N.; Turner, R.E.; Justić, D.; Wiseman, W.J. Predicting the response of Gulf of Mexico hypoxia to variations in Mississippi River nitrogen load. Limnol. Oceanogr. 2003, 48, 951–956. [Google Scholar] [CrossRef] [Green Version]
  24. Turner, R.E.; Rabalais, N.N.; Justic, D. Predicting summer hypoxia in thenorthern Gulf of Mexico: Riverine N, P, and Si loading. Mar. Pollut. Bull. 2006, 52, 139–148. [Google Scholar] [CrossRef]
  25. Mississippi River/Gulf of Mexico Watershed Nutrient Task Force. Action plan for reducing, mitigating, and controlling hypoxia in the northern Gulf of Mexico. US Environmental Protection Agency, Mississippi River/Gulf of Mexico Watershed Nutrient Task Force: Washington, DC, USA, 2001. Available online: https://permanent.fdlp.gov/lps119826/LPS119826.pdf (accessed on 17 February 2022).
  26. Rabalais, N.N.; Turner, T.; Gupta, B.S.; Boesch, D.; Chapman, P.; Murrell, M. Hypoxia in the northern Gulf of Mexico: Does the science support the plan to reduce, mitigate, and control hypoxia? Estuaries Coasts 2007, 30, 753–772. [Google Scholar] [CrossRef]
  27. Sinha, E.; Michalak, A.; Balaji, V. Eutrophication will increase during the 21st century as a result of precipitation changes. Science 2017, 357, 405–408. [Google Scholar] [CrossRef] [Green Version]
  28. Foley, J.A.; Kucharik, C.J.; Twine, T.E.; Coe, M.T.; Donner, S.D. Land use, land cover, and climate change across the Mississippi basin: Impacts on selected land and water resources. Ecosyst. Land Use Chang. 2004, 15, 249–261. [Google Scholar]
  29. Ren, W.; Tian, H.; Cai, W.J.; Lohrenz, S.E.; Hopkinson, C.S.; Huang, W.J.; Yang, J.; Tao, B.; Pan, S.; He, R. Century-long increasing trend and variability of dissolved organic carbon export from the Mississippi River basin driven by natural and anthropogenic forcing. Glob. Biogeochem. Cycles 2016, 30, 1288–1299. [Google Scholar] [CrossRef]
  30. Van Meter, K.J.; Basu, N.B.; Van Cappellen, P. Two centuries of nitrogen dynamics: Legacy sources and sinks in the Mississippi and Susquehanna River Basins. Glob. Biogeochem. Cycles 2017, 31, 2–23. [Google Scholar] [CrossRef] [Green Version]
  31. Robertson, D.M.; Saad, D.A.; Christiansen, D.E.; Lorenz, D.J. Simulated impacts of climate change on phosphorus loading to Lake Michigan. J. Great Lakes Res. 2016, 42, 536–548. [Google Scholar] [CrossRef] [Green Version]
  32. Hungate, B.A.; Dukes, J.S.; Shaw, M.R.; Luo, Y.; Field, C.B. Nitrogen and climate change. Science 2003, 302, 1512–1513. [Google Scholar] [CrossRef] [Green Version]
  33. Wiley, M.J.; Hyndman, D.W.; Pijanowski, B.C.; Kendall, A.D.; Riseng, C.; Rutherford, E.S.; Cheng, S.T.; Carlson, M.L.; Tyler, J.A.; Stevenson, R.J.; et al. A multi-modeling approach to evaluating climate and land use change impacts in a Great Lakes River Basin. Hydrobiologia 2010, 657, 243–262. [Google Scholar] [CrossRef]
  34. USEPA United States Environmental Protection Agency. Watershed Modeling to Assess the Sensitivity of Streamflow, Nutrient, and Sediment Loads to Potential Climate Change and Urban Development in 20 U.S. Watersheds; National Center for Environmental Assessment: Washington, DC, USA, 2013. Available online: http://www.epa.gov/ncea (accessed on 17 February 2022).
  35. Chang, H. Water Quality Impacts of Climate and Land Use Changes in Southeastern Pennsylvania. Prof. Geogr. 2004, 56, 240–257. [Google Scholar]
  36. Ulen, B.; Johansson, G. Long-term nutrient leaching from a Swedish arable field with intensified crop production against a background of climate change. Acta Agric. Scand. Sect. B Soil Plant Sci. 2009, 59, 157–169. [Google Scholar] [CrossRef]
  37. Stuart, M.E.; Gooddy, D.C.; Bloomfield, J.P.; Williams, A.T. A review of the impact of climate change on future nitrate concentrations in groundwater of the UK. Sci. Total Environ. 2011, 409, 2859–2873. [Google Scholar] [CrossRef] [Green Version]
  38. Suddick, E.C.; Davidson, E.A. The Role of Nitrogen in Climate Change and the Impacts of Nitrogen-Climate Interactions on Terrestrial and Aquatic Ecosystems, Agriculture, and Human Health in the United States: A Technical Report Submitted to the US National Climate Assessment; North American Nitrogen Center of the International Nitrogen Initiative (NANC-INI); Woods Hole Research Center: Falmouth, ME, USA, 2012. [Google Scholar]
  39. Loecke, T.D.; Burgin, A.J.; Riveros-Iregui, D.A.; Ward, A.S.; Thomas, S.A.; Davis, C.A.; St Clair, M.A. Weather whiplash in agricultural regions drives deterioration of water quality. Biogeochemistry 2017, 133, 7–15. [Google Scholar] [CrossRef] [Green Version]
  40. Strickling, H.L.; Obenour, D.R. Leveraging Spatial and Temporal Variability to Probabilistically Characterize nutrient sources and export rates in a developing watershed. Water Resour. Res. 2018, 54, 5143–5162. [Google Scholar] [CrossRef]
  41. Iqbal, J.; Necpalova, M.; Archontoulis, S.V.; Anex, R.P.; Bourguignon, M.; Herzmann, D.; Mitchell, D.C.; Sawyer, J.E.; Zhu, Q.; Castellano, M.J. Extreme weather-year sequences have nonadditive effects on environmental nitrogen losses. Glob. Chang. Biol. 2018, 24, e303–e317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Sinha, E.; Michalak, A.M. Precipitation dominates interannual variability ofriverine nitrogen loading across the continental United States. Environ. Sci. Technol. 2016, 50, 12874–12884. [Google Scholar] [CrossRef] [PubMed]
  43. Lu, C.; Zhang, J.; Tian, H.; Crumpton, W.G.; Helmers, M.J.; Cai, W.-J.; Hopkinson, C.S.; Lohrenz, S.E. Increased extreme precipitation challenges nitrogen load management to the Gulf of Mexico. Commun. Earth Environ. 2020, 1, 1–10. [Google Scholar] [CrossRef]
  44. Xia, X.H.; Zhang, S.B.; Li, S.L.; Zhang, L.W.; Wang, G.Q.; Zhang, L. The cycle of nitrogen in river systems: Sources, transformation, and flux. Environ. Sci. Process. Impacts 2018, 20, 863–891. [Google Scholar] [CrossRef]
  45. Ide, J.; Takeda, I.; Somura, H.; Mori, Y.; Sakuno, Y.; Yone, Y.; Takahashi, E. Impacts of hydrological changes on nutrient transport from diffuse sources in a rural river basin, western Japan. J. Geophys. Res. Biogeosci. 2019, 124, 2565–2581. [Google Scholar] [CrossRef]
  46. Thompson, E. How land use affects nutrient pollution in a changing climate. EOS 2019, 100. [Google Scholar] [CrossRef]
  47. Wu, L.; Long, T.; Liu, X.; Guo, J. Impacts of climate and land-use changes on the migration of non-point source nitrogen and phosphorus during rainfall-runoff in the Jialing River Watershed, China. J. Hydrol. 2012, 475, 26–41. [Google Scholar] [CrossRef]
  48. El-Khoury, A.; Seidou, O.; Lapen, D.; Que, Z.; Mohammadian, M.; Sunohara, M.; Bahram, D. Combined impacts of future climate and land use changes on discharge, nitrogen and phosphorus loads for a Canadian river basin. J. Environ. Manag. 2015, 151, 76–86. [Google Scholar] [CrossRef] [Green Version]
  49. Parsons, A.J.; Wainwright, J.; Powell, D.M.; Kaduk, J.; Brazier, R.E. A conceptual model for determining soil erosion by water. Earth Surf. Process. Landf. 2004, 29, 1293–1302. [Google Scholar] [CrossRef]
  50. Parsons, A.J.; Brazier, R.E.; Wainwright, J.; Powell, D.M. Scale relationships in hillslope runoff and erosion. Earth Surf. Process. Landf. 2006, 31, 1384–1393. [Google Scholar] [CrossRef]
  51. USGS United States Geological Survey. 2022. Available online: https://www.sciencebase.gov/catalog/item/617987dbd34ea58c3c6fa16e (accessed on 17 February 2022).
  52. National Oceanic and Atmospheric Administration (NOAA). 2022. Available online: https://www.psl.noaa.gov/data/gridded/ (accessed on 17 February 2022).
  53. Subramanya, K. Engineering Hydrology, 3rd ed.; McGraw-Hill Education: New Delhi, India, 2009; p. 450. [Google Scholar]
  54. Godsey, S.E.; Kirchner, J.W.; Clow, D.W. Concentration–discharge relationships reflect chemostatic characteristics of US catchments. Hydrol. Process. 2009, 23, 1844–1864. [Google Scholar] [CrossRef]
  55. PennState. 2018. Available online: https://online.stat.psu.edu/stat462/node/180/ (accessed on 25 February 2022).
  56. Mi, X.C.; Zhang, J.-T.; Zhang, F.; Shangguan, T.L. Analysis of relationships between vegetation and climate in Shanxi Plateau. Acta Phytoecol. Sin. 1996, 20, 549–560, (In Chinese with English abstract). [Google Scholar]
  57. Zhang, J.-T. A study on relations of vegetation, climate and soils in Shanxi province, China. Plant Ecol. 2002, 162, 23–31. [Google Scholar]
  58. Zhang, J.-T.; Ru, W.; Li, B. Relationships between vegetation and climate on The loess plateau in china. Folia Geobot. 2006, 41, 151–163. [Google Scholar] [CrossRef]
  59. Goward, S.N.; Prince, S.D. Transient effects of climate on vegetation dynamics: Satellite observations. J. Biogeogr. 1995, 22, 549–563. [Google Scholar] [CrossRef]
  60. Fan, X.; Ma, Z.; Yang, Q.; Han, Y.; Mahmood, R. Land use/land cover changes and regional climate over the Loess Plateau during 2001–2009. Part II: Interrelationship from observations. Clim. Chang. 2015, 129, 441–455. [Google Scholar] [CrossRef] [Green Version]
  61. Kishtawal, C.M.; Niyogi, D.; Tewari, M.; PielkeSr, R.A.; Shepherd, J.M. Urbanization signature in the observed heavy rainfall climatology over India. Int. J. Climatol. 2010, 30, 1908–1916. [Google Scholar] [CrossRef] [Green Version]
  62. Costa, M.H.; Yanagi, S.N.; Oliveira, P.J.; Ribeiro, A.; Rocha, E.J. Climate change in Amazonia caused by soybean cropland expansion, as compared to caused by pastureland expansion. Geophys. Res. Lett. 2007, 34, L07706. [Google Scholar] [CrossRef]
  63. Sampaio, G.; Nobre, C.; Costa, M.; Satyamurty, P.; Soares-Filho, B.; Cardoso, M. Regional climate change over eastern Amazonia caused by pasture and soybean cropland expansion. Geophys. Res. Lett. 2007, 34, L17709. [Google Scholar] [CrossRef] [Green Version]
  64. Notaro, M.; Liu, Z. Observed vegetation–climate feedbacks in the United States. J. Clim. 2006, 19, 763–786. [Google Scholar] [CrossRef] [Green Version]
  65. Clark, C.A.; Arritt, R.W. Numerical Simulations of the Effect of Soil Moisture and Vegetation Cover on the Development of Deep Convection. J. Appl. Meteorol. Climatol. 1995, 34, 2029–2045. [Google Scholar] [CrossRef] [Green Version]
  66. Sud, Y.C.; Mocko, D.M.; Walker, G.K. Influence of land surface fluxes on precipitation: Inferences from simulations forced with four ARM-CART SCM datasets. J. Clim. 2001, 14, 3666–3691. [Google Scholar] [CrossRef]
  67. Freedman, J.M.; Fitzjarrald, D.R.; Moore, K.E.; Sakai, R.K. Boundary layer clouds and vegetation–atmosphere feedbacks. J. Clim. 2001, 14, 180–197. [Google Scholar] [CrossRef]
  68. McPherson, R.A. A review of vegetation–atmosphere interactions and their influences on mesoscale phenomena. Prog. Phys. Geogr. 2007, 31, 261–285. [Google Scholar] [CrossRef]
  69. Zhou, Y.; Yang, Z.; Zhang, D.; Jin, X.; Zhang, J. Inter-catchment comparison of flow regime between the Hailiutu and Huangfuchuan rivers in the semi-arid Erdos Plateau, Northwest China. Hydrol. Sci. J. 2015, 60, 688–705. [Google Scholar] [CrossRef]
  70. Chai, Y.; Zhu, B.; Yue, Y.; Yang, Y.; Li, S.; Ren, J.; Xiong, H.; Cui, X.; Yan, X.; Li, Y. Reasons for the homogenization of the seasonal discharges in the Yangtze River. Hydrol. Res. 2020, 51, 470–483. [Google Scholar] [CrossRef] [Green Version]
  71. Jaynes, D.B.; Colvin, T.S.; Karlen, D.L.; Cambardella, C.A.; Meek, D.W. Nitrate Loss in Subsurface Drainage as Affected by Nitrogen Fertilizer Rate. J. Environ. Qual. 2001, 30, 1305–1314. [Google Scholar] [CrossRef]
  72. Xu, Z.; Zhang, X.; Xie, J.; Yuan, G.; Tang, X.; Sun, X.; Yu, G. Total Nitrogen Concentrations in Surface Water of Typical Agro- and Forest Ecosystems in China, 2004–2009. PLoS ONE 2014, 9, e92850. [Google Scholar] [CrossRef]
  73. Chichester, F.W.; Richardson, C.W. Sediment and nutrient loss from clay soils as affected by tillage. J. Environ. Qual. 1992, 21, 587–590. [Google Scholar] [CrossRef]
  74. Franklin, D.H.; Truman, C.C.; Potter, T.L.; Bosch, D.D.; Strickland, T.C.; Jenkins, M.B.; Nuti, R.C. Nutrient losses in runoff from conventional and no-till pearl millet on pre-wetted Ultisols fertilized with broiler litter. Agric. Water Manag. 2012, 113, 38–44. [Google Scholar] [CrossRef]
  75. Issaka, F.; Zhang, Z.; Zhao, Z.Q.; Asenso, E.; Li, J.H.; Li, Y.T.; Wang, J.J. Sustainable conservation tillage improves soil nutrients and reduces nitrogen and phosphorous losses in maize farmland in southern China. Sustainability 2019, 11, 2397. [Google Scholar] [CrossRef] [Green Version]
  76. Valiela, I.; Bowen, J.L. Nitrogen sources to watersheds and estuaries: Role of land cover mosaics and losses within watersheds. Environ. Pollut. 2002, 118, 239–248. [Google Scholar] [CrossRef]
  77. McClain, M.E.; Richey, J.E.; Pimentel, T.P. Groundwater nitrogen dynamics at the terrestriallotic interface of a small catchment in the Central Amazon Basin. Biogeochemistry 1994, 27, 113–127. [Google Scholar] [CrossRef]
  78. Hedin, L.O.; von Fischer, J.C.; Ostrom, N.E.; Kennedy, B.P.; Brown MGRobertson, G.P. Thermodynamic constraints on nitrogen transformations and other biogeochemical processes at soil stream interfaces. Ecology 1998, 79, 684–703. [Google Scholar] [CrossRef] [Green Version]
  79. Hill, A.R.; Devito, K.J.; Campagnolo, S.; Sanmugadas, K. Subsurface denitrification in a forest riparian zone: Interactions between hydrology and supplies of nitrate and organic carbon. Biogeochemistry 2000, 51, 193–223. [Google Scholar] [CrossRef]
  80. Groffman, P.M.; Gold, A.J.; Simmons, R.C. Nitrate dynamics in riparian forests: Microbial studies. J. Environ. Qual. 1992, 21, 666–671. [Google Scholar] [CrossRef] [Green Version]
  81. Groffman, P.M.; Holland, E.; Myrold, D.D.; Robertson, G.P.; Zou, X. Denitrification. In Standard Soil Methods for Long Term Ecological Research; Robertsonn, G.P., Bledsoe, C.S., Coleman, D.C., Sollins, P., Eds.; Oxford University Press: New York, NY, USA, 1999; pp. 272–288. [Google Scholar]
  82. Lowrance, R. Groundwater nitrate and denitrification in a coastal plain riparian forest. J. Environ. Qual. 1992, 21, 401–405. [Google Scholar] [CrossRef]
  83. Horel, A.; Potyó, I.; Szili-Kovács, T.; Molnár, S. Potential nitrogen fixation changes under different land uses as influenced by seasons and biochar amendments. Arab. J. Geosci. 2018, 11, 559. [Google Scholar] [CrossRef] [Green Version]
  84. White, M.D.; Greer, K.A. The effects of watershed urbanization on the stream hydrology and riparian vegetation of Los Peñasquitos Creek, California. Landsc. Urban Plan. 2006, 74, 125–138. [Google Scholar] [CrossRef]
  85. USGS United States Geological Survey. The Quality of Our Nation’s Waters—Nutrients and Pesticides; U.S. Geological Survey Circular: Reston, VA, USA, 1999; Volume 1225, p. 82. Available online: https://pubs.usgs.gov/circ/1999/1225/report.pdf (accessed on 17 February 2022).
  86. USGS United States Geological Survey. 2012. Available online: https://sparrow.wim.usgs.gov/sparrow-midwest-2012/ (accessed on 17 February 2022).
  87. Schlesinger, W.H.; Ward, T.J.; Anderson, J. Nutrient losses in runoff from grassland and shrubland habitats in Southern New Mexico: II. Field plots. Biogeochemistry 2000, 49, 69–86. [Google Scholar] [CrossRef]
  88. Bertol, I.; Mello, E.L.; Guadagnin, J.C.; Zaparolli, A.L.V.; Carrafa, M.R. Nutrients losses by water erosion. Sci. Agric. 2003, 60, 581–586. [Google Scholar] [CrossRef]
  89. Guadagnin, J.C.; Bertol, I.; Cassol, P.C.; Amaral, A.J. Soil, water and nitrogen losses through erosion under different tillage systems. Rev. Bras. Cienc. Solo 2005, 29, 277–286. [Google Scholar] [CrossRef]
  90. Berhe, A.A.; Arnold, C.; Stacy, E.; Lever, R.; McCorkle EAraya, S.N. Soil erosion controls on biogeochemical cycling of carbon and nitrogen. Nat. Educ. Knowl. 2014, 5, 2. [Google Scholar]
  91. Bramorski, J.; Trivelin, P.C.O.; Crestana, S. Nitrogen loss by erosion from mechanically tilled and untilled soil under successive simulat-ed rainfalls. Rev. Bras. Cienc. Solo. 2015, 39, 1204–1211. [Google Scholar] [CrossRef] [Green Version]
  92. de Vente, J.; Poesen, J. Predicting soil erosion and sediment yield at the basin scale: Scale issues and semi-quantitative models. Earth Sci. Rev. 2005, 71, 95–125. [Google Scholar] [CrossRef]
  93. de Vente, J.; Poesen, J.; Arabkhedri, M.; Verstraeten, G. The sediment delivery problem revisited. Prog. Phys. Geogr. 2007, 31, 155–178. [Google Scholar] [CrossRef] [Green Version]
  94. Osterkamp, W.R.; Toy, T.J. Geomorphic considerations for erosion prediction. Environ. Geol. 1997, 29, 152–157. [Google Scholar] [CrossRef]
  95. Walling, D.E. The sediment delivery problem. J. Hydrol. 1983, 65, 209–237. [Google Scholar] [CrossRef]
  96. Syvitski, J.P. Supply and flux of sediment along hydrological pathways: Research for the 21st century. Glob. Planet. Chang. 2003, 39, 1–11. [Google Scholar] [CrossRef]
  97. Milliman, J.D.; Meade, R.H. World-wide delivery of river sediment to the oceans. J. Geol. 1983, 91, 1–21. [Google Scholar] [CrossRef]
  98. Milliman, J.D.; Syvitski, J.P. Geomorphic/tectonic control of sediment discharge to the ocean: The importance of small mountainous rivers. J. Geol. 1992, 100, 525–544. [Google Scholar] [CrossRef]
  99. Lane, L.J.; Hernandez, M.; Nichols, M. Processes controlling sediment yield from watersheds as function of spatial scale. Environ. Model. Softw. 1997, 12, 355–369. [Google Scholar] [CrossRef]
  100. Mutema, M.; Chaplot, V.; Jewitt, G.; Chivenge, P.; Bloschl, G. Annual water, sediment, nutrient, and organic carbon fluxes in river basins: A global meta-analysis as a function of scale. Water Resour. Res. 2015, 51, 8949–8972. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Study sites.
Figure 1. Study sites.
Water 14 00865 g001
Figure 2. Instantaneous TN yields (mg TN m−2 min−1) versus instantaneous discharge (L min−1) for the studied rivers. (a) is for the first panel, and (b) for the second panel. * Significant correlations (p < 0.05).
Figure 2. Instantaneous TN yields (mg TN m−2 min−1) versus instantaneous discharge (L min−1) for the studied rivers. (a) is for the first panel, and (b) for the second panel. * Significant correlations (p < 0.05).
Water 14 00865 g002aWater 14 00865 g002b
Figure 3. TN yield rates (i.e., the slopes of TN yield—Q regression) versus catchment size.
Figure 3. TN yield rates (i.e., the slopes of TN yield—Q regression) versus catchment size.
Water 14 00865 g003
Table 1. Pre-industrial and contemporary mobilisable nitrogen loadings [17].
Table 1. Pre-industrial and contemporary mobilisable nitrogen loadings [17].
Pre-Industrial Nmob (Tg/Year)Contemporary Nmob (Tg/Year)
DepositionFixationTotalDepositionFixationFertiliserLivestock LoadPeople LoadTotalWorld Population Share (%) *
Africa3.6331.9935.616.5825.020.946.432.2541.2217.20
Asia3.2925.4528.7311.2122.6220.2122.4112.789.1559.54
Australia0.466.997.450.465.70.191.480.097.910.33
Europe0.623.924.544.43.065.4810.133.0926.169.59
North America1.279.8111.076.168.765.485.851.9528.217.60
Oceania0.020.340.350.030.170.070.580.020.870.21
South America2.7520.1622.913.5116.121.596.631.2129.065.53
Totals12991113281345421223100
* [18].
Table 2. Study rivers, catchment area (km2), precipitation (mm), percentage of land cover (%), and the correlations (r2) of discharge level (L min−1) with both TN concentrations (mg L−1) and TN yields (mg m−2 min−1).
Table 2. Study rivers, catchment area (km2), precipitation (mm), percentage of land cover (%), and the correlations (r2) of discharge level (L min−1) with both TN concentrations (mg L−1) and TN yields (mg m−2 min−1).
No.RiverArea (km2)Mean Annual Precipitation (m)Pasture/
Grassland
ForestBarrenAgricultureUrbanWater/
Wetland
Shrub/
Scrub
Correlation (r2) between Q-CCorrelation (r2) between Q-Y
1Little Arkansas River near Sedgwick, KS32100.7423.73.70.161.59.41.40.20.02+0.96+
2Grand River near Sumner, MO17,8200.9147.217.70.226.44.73.40.40.18+0.69+
3Elkhorn River at Waterloo, NE17,8690.6431.21.60.257.84.24.80.20.71+0.97+
4South Platte River near Kersey, CO25,0190.3730.530.80.64.48.94.220.60.59−0.90+
5White River at Hazleton, IN29,2790.969.531.90.344.111.72.20.30.28+0.92+
6Iowa River at Wapello, IA32,3690.878.64.20.275.97.23.70.20.54+0.94+
7Yazoo River below Steele Bayou near Long Lake, MS34,5901.2210.526.80.239.15.415.92.10.09+0.84+
8North Canadian River near Harrah, OK35,6800.6748.610.40.2245.3110.50.31−0.99+
9Des Moines River at Keosauqua, IA36,3600.8513.49.40.265.97.63.20.30.61+0.92+
10Osage River near St. Thomas, MO37,7690.9745.231.40.212.65.84.30.50.14+0.94+
11Illinois River at Valley City, IL69,2590.924.911.70.364.814.23.80.30.41+0.95+
12Wabash River at New Harmony, IN75,7200.976.320.50.2619.22.50.30.24+0.91+
13Mississippi River at Hastings, MN96,0900.699.1160.3467.120.80.70.51+0.94+
14Tennessee River at Highway 60 near Paducah, KY104,4491.1821.958.50.33.49.74.12.10.45+0.93+
15Kansas River at DeSoto, KS154,7700.6439.520.25340.90.40.58+0.97+
16Red River at Alexandria, LA174,8190.8636.621.40.514.54.7616.30.26+0.99+
17Yellowstone River near Sidney, MT178,9190.363312.513.71.32.146.40.31+0.89+
18Platte River at Louisville, NE221,1100.4150.88.90.315.13.53.4180.52+0.96+
19Mississippi River at Clinton, IA221,7100.8410.726.30.236.56.319.20.80.13+0.83+
20Ohio River at Cannelton Dam at Cannelton, IN251,2300.9916.959.60.510.69.91.60.90.58+0.98+
21Arkansas River at David D Terry Lock and Dam below Little Rock, AR410,3300.6544.115.20.320.94.61.713.20.29+0.97+
22Mississippi River Below Grafton, IL443,6700.8710.618.50.250.48.111.70.50.56+0.93+
23Ohio River at Olmsted, IL525,7701.0517.153.70.417.87.52.510.42+0.96+
24Missouri River at Omaha, NE836,0500.46428.50.622.32.33.520.80.16+0.37+
25Missouri River at Hermann, MO1,353,3700.50439.40.525.83.23.214.90.39+0.88+
26Mississippi River at Thebes, IL1,847,1790.5934.411.50.432.74.55.5110.62+0.94+
Note: Bold numbers denote significant correlations (p < 0.05). + Positive correlations. − Negative correlations.
Table 3. Catchment area, Qarea, TN concentrations, and TN yields of the study sites.
Table 3. Catchment area, Qarea, TN concentrations, and TN yields of the study sites.
Catchment No.Area (km2)Median Qarea
(mm min−1)
Median Concentrations (mg TN L−1)Median Yields
(mg TN m−2 min−1)
132103.97 × 1052.6368.86 × 105
217,8202.78 × 1042.5519.28 × 104
317,8691.62 × 1047.1521.30 × 103
425,0196.29 × 1056.0733.90 × 104
529,2798.37 × 1042.8342.48 × 103
632,3696.01 × 1048.0764.84 × 103
734,5907.76 × 1041.4201.07 × 103
835,6801.05 × 1055.3945.27 × 105
936,3604.24 × 1049.2683.66 × 103
1037,7693.12 × 1040.7822.48 × 104
1169,2596.05 × 1045.1683.16 × 103
1275,7207.49 × 1044.1943.13 × 103
1396,0903.19 × 1044.5001.40 × 103
14104,4497.51 × 1040.6634.97 × 104
15154,7703.11 × 1052.2951.62 × 104
16174,8191.80 × 1041.0521.89 × 104
17178,9197.37 × 1050.8456.03 × 105
18221,1105.79 × 1053.8042.17 × 104
19221,7104.37 × 1042.9291.39 × 103
20251,2308.51 × 1041.9061.54 × 103
21410,3301.12 × 1040.9431.04 × 104
22443,6704.93 × 1044.0802.05 × 103
23525,7708.75 × 1041.8221.76 × 103
24836,0506.76 × 1052.5361.71 × 104
251,353,3709.02 × 1052.9593.01 × 104
261,847,1792.27 × 1043.4938.74 × 104
Table 4. Correlation (r) between catchment area (km2), precipitation (m), and percentage of land cover (%).
Table 4. Correlation (r) between catchment area (km2), precipitation (m), and percentage of land cover (%).
Area (km2)Mean Annual Precipitation (m)Bare LandForestWetland/WaterGrassland/PastureScrub/ShrubUrban
Area (km2)1.00
Mean annual precipitation (m)−0.331.00
Bare land0.30−0.501.00
Forest−0.100.580.151.00
Wetland/Water−0.040.19−0.180.051.00
Grassland/Pasture0.26−0.550.22−0.27−0.421.00
Scrub/Shrub0.25−0.710.86−0.19−0.210.481.00
Urban−0.360.55−0.330.41−0.01−0.70−0.571.00
Note: Bold numbers indicate significant correlations (p < 0.05).
Table 5. Correlation (r) of Qarea, TN concentrations, TN yields with catchment size, precipitation, and land cover types.
Table 5. Correlation (r) of Qarea, TN concentrations, TN yields with catchment size, precipitation, and land cover types.
Median Qarea
(mm min−1)
Median TN Concentration (mg L−1)Median TN Yield
(mg m−2 min−1)
Area (km2)−0.18−0.16−0.20
Mean annual precipitation (m)0.83−0.180.43
Barren −0.21−0.29−0.32
Agriculture0.180.620.70
Forest0.67−0.44−0.01
Wetland/Water0.160.000.07
Grassland/Pasture−0.75−0.29−0.76
Scrub/Shrub−0.54−0.23−0.49
Urban0.640.210.55
Note: Bold numbers indicate significant correlations (p < 0.05).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Allafta, H.; Opp, C. Understanding the Combined Effects of Land Cover, Precipitation and Catchment Size on Nitrogen and Discharge—A Case Study of the Mississippi River Basin. Water 2022, 14, 865. https://doi.org/10.3390/w14060865

AMA Style

Allafta H, Opp C. Understanding the Combined Effects of Land Cover, Precipitation and Catchment Size on Nitrogen and Discharge—A Case Study of the Mississippi River Basin. Water. 2022; 14(6):865. https://doi.org/10.3390/w14060865

Chicago/Turabian Style

Allafta, Hadi, and Christian Opp. 2022. "Understanding the Combined Effects of Land Cover, Precipitation and Catchment Size on Nitrogen and Discharge—A Case Study of the Mississippi River Basin" Water 14, no. 6: 865. https://doi.org/10.3390/w14060865

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop