Optimization of the N footprint model and analysis of nitrogen pollution in irrigation areas: A case study of Ningxia Hui Autonomous Region, China.

Water diverted from rivers for irrigation areas often contains large amounts of nitrogen (N), which is frequently overlooked and its role in contributing to N pollution is unknown. To investigate the influence of water diversion on N in different systems within irrigation areas, we developed and optimized the N footprint model, taking into account the N carried by irrigation water diversion and drainage in irrigated areas. This optimized model can serve as a reference for evaluating N pollution in other irrigated areas. By analyzing 29 years (1991-2019) of statistical data from a diverted irrigation area in Ningxia Hui Autonomous Region (Ningxia), China, the study assessed the contribution of water diversion to N in agriculture, animal husbandry, and human domestic activities. The results demonstrated that water diversion and drainage accounted for 10.3% and 13.8% in whole system, of the total N input and output in Ningxia, highlighting the potential N pollution risks associated with these activities. Additionally, the use of fertilizers in the plant subsystem, feed in the animal subsystem, and sanitary sewage in the human subsystem represented the main sources of N pollution in each subsystem. On a temporal scale, the study found that N loss increased year by year before reaching a stable level, indicating that N loss had reached its peak in Ningxia. The correlation analysis suggested that rainfall could regulate N input and output in irrigated areas by showing a negative correlation with water diversion, agricultural water consumption, and N from irrigated areas. Moreover, the study revealed that the amount of N brought by water diverted from rivers for irrigation should be taken into account when calculating the amount of fertilizer N required in the irrigation area.

Comparing critical source areas for the sediment and nutrients of calibrated and uncalibrated models in a plateau watershed in southwest China.

Controlling non-point source pollution is often difficult and costly. Therefore, focusing on areas that contribute the most, so-called critical source areas (CSAs), can have economic and ecological benefits. CSAs are often determined using a modelling approach, yet it has proved difficult to calibrate the models in regions with limited data availability. Since identifying CSAs is based on the relative contributions of sub-basins to the total load, it has been suggested that uncalibrated models could be used to identify CSAs to overcome data scarcity issues. Here, we use the SWAT model to study the extent to which an uncalibrated model can be applied to determine CSAs. We classify and rank sub-basins to identify CSAs for sediment, total nitrogen (TN), and total phosphorus (TP) in the Fengyu River Watershed (China) with and without model calibration. The results show high similarity (81%-93%) between the identified sediment and TP CSA number and locations before and after calibration both on the yearly and seasonal scale. For TN alone, the results show moderate similarity on the yearly scale (73%). This may be because, in our study area, TN is determined more by groundwater flow after calibration than by surface water flow. We conclude that CSA identification with the uncalibrated model for TP is always good because its CSA number and locations changed least, and for sediment, it is generally satisfactory. The use of the uncalibrated model for TN is acceptable, as its CSA locations did not change after calibration; however, the TN CSA number changed by over 60% compared to the figures before calibration on both yearly and seasonal scales. Therefore, we advise using an uncalibrated model to identify CSAs for TN only if water yield composition changes are expected to be limited. This study shows that CSAs can be identified based on relative loading estimates with uncalibrated models in data-deficient regions.

Spatio-temporal variation of net anthropogenic nitrogen inputs (NANI) from 1991 to 2019 and its impacts analysis from parameters in Northwest China.

At present, excessive nutrient inputs caused by human activities have resulted in environmental problems such as agricultural non-point source pollution and water eutrophication. The Net Anthropogenic Nitrogen Inputs (NANI) model can be used to estimate the nitrogen (N) inputs to a region that are related to human activities. To explore the net nitrogen input of human activities in the main grain-producing areas of Northwestern China, the county-level statistical data for the Ningxia province and NANI model parameters were collected, the spatio-temporal distribution characteristics of NANI were analyzed and the uncertainty and sensitivity of the parameters for each component of NANI were quantitatively studied. The results showed that: (1) The average value of NANI in Ningxia from 1991 to 2019 was 7752 kg N km-2 yr-1. Over the study period, the inputs first showed an overall increase, followed by a decrease, and then tended to stabilize. Fertilizer N application was the main contributing factor, accounting for 55.6%. The high value of NANI in Ningxia was mainly concentrated in the Yellow River Diversion Irrigation Area. (2) The 95% confidence interval of NANI obtained by the Monte Carlo approach was compared with the results from common parameters in existing literature. The simulation results varied from -6.4% to 27.4% under the influence of the changing parameters. Net food and animal feed imports were the most uncertain input components affected by parameters, the variation range was -20.7%-77%. (3) The parameters of inputs that accounted for higher proportions of the NANI were more sensitive than the inputs with lower contributions. The sensitivity indexes of the parameters contained in the fertilizer N applications were higher than those of net food and animal feed imports and agricultural N-fixation. This study quantified the uncertainty and sensitivity of parameters in the process of NANI simulation and provides a reference for global peers in the application and selection of parameters to obtain more accurate simulation results.

Estimation of the nonpoint source nitrogen load in a strongly disturbed watershed of the North China Plain.

Identification of nonpoint source (NPS) pollution is a great challenge in the North China Plain, which has modified rivers and insufficient data. In this study, a simple and reasonable method was developed to estimate the total nitrogen (TN) load in rural areas of the North China Plain. The method was found to work well and produce results consistent with monitoring data when considering various TN sources and transfer mechanisms. The annual TN loads from rural living, livestock and the farmlands were 121.9 × 10(3), 45.6 × 10(3) and 78.5 × 10(3) kg/yr, respectively. The TN load in the region along the river contributed much more to the NPS pollution than that in areas far from the river, with average TN loads of approximately 3394 and 602 kg km(-2) yr(-1), respectively. Overall, the results indicate that this method is suitable for NPS load estimates in severely disturbed watersheds with insufficient data.

Contrasting patterns of nutrient dynamics during different storm events in a semi-arid catchment of northern China.

Nutrient discharge during storm events is a critical pathway for nutrient export in semi-arid catchments. We investigated nutrient dynamics during three summer storms characterized by different rainfall magnitude in 2012 in a semi-arid catchment of northern China. The results showed that, in response to storm events, nutrient dynamics displayed big variation in temporal trends of nutrient concentration and in nutrient concentration-flow discharge relationships. Nutrient concentrations had broader fluctuations during an extreme storm than during lesser storms, whereas the concentration ranges of the a moderate storm were no broader than those of a smaller one. The different concentration fluctuations were caused by storm magnitude and intensity coupled with the antecedent rainfall amount and cumulative nutrients. Correlation coefficients between nutrient concentrations and flow discharge varied from positive to negative for the three different events. There were no consistent hysteresis effects for the three different events, and no hysteresis effects were observed for any of the variables during the moderate storm (E2). Our findings provide useful information for better understanding nutrient loss mechanisms during storm events in semi-arid areas of a monsoon climate region.

Investigating hydrological transport pathways of dissolved organic carbon in cold region watershed based on a watershed biogeochemical model.

Dissolved organic carbon (DOC) is a significant component of regional and global carbon cycles and an important surface water quality indicator. DOC affects the processes of solubility, bioavailability and transport for a number of contaminants, such as heavy metals. Therefore, it is crucial to understand DOC fate and transport in the watershed and the transport pathways of DOC load. We modified a previously developed watershed-scale organic carbon model by incorporating the DOC load from glacier melt runoff and used the modified model to simulate periodic daily DOC load in the upper Athabasca River Basin (ARB) in the cold region of western Canada. The calibrated model achieved an overall acceptable performance for simulating daily DOC load with model uncertainties mainly from the underestimation of peak loads. Parameter sensitivity analysis indicates that the fate and transport of DOC load in upper ARB are mainly controlled by DOC production in the soil layers, DOC transport at the soil surface, and reactions in the stream system. The modeling results indicated that the DOC load is mainly from the terrestrial sources and the stream system was a negligible sink in the upper ARB. It also indicated that rainfall-induced surface runoff was the major transport pathway of DOC load in the upper ARB. However, the DOC loads transported by glacier melt runoff were negligible and only accounted for 0.02% of the total DOC loads. In addition, snowmelt-induced surface runoff and lateral flow contributed 18.7% of total DOC load, which is comparable to the contribution from the groundwater flow. Our study investigated the DOC dynamics and sources in the cold region watershed in western Canada and quantified the contribution of different hydrological pathways to DOC load, which could provide a useful reference and insight for understanding watershed-scale carbon cycle processes.

Long-term variations of water quality and nutrient load inputs in a large shallow lake of Yellow River Basin: Implications for lake water quality improvements.

Many studies have investigated water quality changes in freshwater lakes, however, studies examining long-term relationships between lake water quality and total nitrogen (TN) and total phosphorus (TP) load inputs and investigating the causes that indicate improvements in water quality are limited. In this study, we utilized the LOADEST model to estimate TN and TP load inputs, assessed lake trophic status using the integrated nutrient index method, and explored trends and relationships between nutrient load inputs and water quality in Wuliangsuhai Lake, a large shallow lake of Yellow River Basin in China. Additionally, we identified the causes for recent water quality improvements and proposed future management strategies to further improve the water quality. Our findings revealed that water quality in Wuliangsuhai Lake of Yellow River basin has been improved mainly due to the abatements of nutrient load inputs from the watershed. Between 2010 and 2020, TN and TP loads from the watershed decreased significantly by 65.12 % and 89.4 %, respectively. TN and TP concentrations also notably decreased across the lake areas, including the inlet (91.21 % and 95.59 %), central (73.49 % and 87.12 %), and outlet (40.68 % and 40.54 %) areas. Correlation analysis confirms a strong positive relationship between lake water quality and nutrient load inputs (excluding the outlet area), highlighting the impact of nutrient inputs on lake water quality. The results indicated that the recent water quality improvements in the lake was mainly because of effective control for point source pollution from industrial wastewater and domestic sewage and the non-point source pollution control holds the potential to further improve the water quality.

A numerical modeling framework for simulating the key in-stream fate processes of PAH decay in Muskeg River Watershed, Alberta, Canada.

Most previous water quality studies oversimplified in-stream processes for modeling the fate and transport of critical organic contaminants, such as Polycyclic Aromatic Hydrocarbons (PAHs). Taking four selected PAHs as representative organic contaminants, we developed a numerical modeling framework using a Water Quality Analysis Simulation Program 8 (WASP8) and a well-established watershed model, i.e., Soil and Water Assessment Tool (SWAT) to: (1) address the influence of in-stream processes, including direct photolysis, volatilization, partitioning of PAHs to suspended solids, and DOC complexation processes on PAH concentrations; and (2) establish relationships between spatiotemporal distribution of environmental factors (e.g., ice coverage, water temperature, wind, and light attenuation), in-stream processes, and PAH concentrations at a watershed scale. Using calibrated SWAT and WASP8 models, we evaluated the impacts of seasonal changes in environmental factors on in-stream processes in the Muskeg River watershed, which is part of the Athabasca Oil Sands Region (AOSR), the third-largest crude oil reserves of the world in western Canada. Among four selected PAHs, simulation results suggest that Naphthalene primarily decay in the water through volatilization or direct photolysis. For Phenanthrene, Pyrene, and Chrysene, DOC complexation, volatilization, and direct photolysis all contribute to their decay in the water, with a strong dependence on seasonality. Model simulations indicated that direct photolysis and volatilization rates are meager in cold seasons, mainly due to low river temperature and ice coverage. However, these processes gradually resume when entering the warm season. In summary, the model simulation results suggest that critical in-stream processes such as direct photolysis, volatilization, and partitioning and their relationship with environmental factors should be considered when simulating the fate and transport of organic contaminants in the river systems. Our results also reveal that the relationship between environmental factors and fate processes affecting PAH concentrations can vary across a watershed and in different seasons.

Hydro-climate and biogeochemical processes control watershed organic carbon inflows: Development of an in-stream organic carbon module coupled with a process-based hydrologic model.

Dissolved organic carbon (DOC) in surface waters directly influences the speciation, transport, and fate of heavy metals, as well as the partitioning of organic contaminants. However, the lack of process-based watershed-scale models for simulating carbon cycling and transport has limited the effective watershed management to control organic carbon fluxes to source waters and throughout the river systems. Here, a process-based in-stream organic carbon (OC) module was developed, coupled with the physically process-based Soil and Water Assessment Tool (SWAT), and linked with its existing soil carbon module to simulate dynamics of both particulate organic carbon (POC) and DOC. The advanced model simulates a large spectrum of OC processes from landscapes to stream networks throughout the watersheds. In-stream organic carbon processes related to POC and DOC as state variables are modeled in the water column, and the transformations between different carbon species and interactions between OC with algae are considered. The module's ability to simulate total organic carbon (TOC) loads was assessed, and the monthly and seasonal variations were captured over 14 years. Simulations for TOC loads suggested that spring snowmelt and summer rainfall runoff events are the main driving forces behind TOC export in the watershed. The parameter sensitivity analysis and dynamic reaction rate simulated in the streams suggested that TOC dynamics in the study area are controlled by both landscape and in-stream processes. The spatiotemporal analysis of the simulated TOC load showed that 55.8% of total terrestrial OC exports into the streams are removed due to in-stream POC settling and DOC mineralization, confirming the necessity of integrating terrestrial and aquatic OC processes for process understanding and for modelling and management of water quality at the watershed scale. The developed OC module is a potentially effective tool for simulating the OC cycle at the watershed scale and can be applied further to water treatment plans and watershed management.

Incorporating a non-reactive heavy metal simulation module into SWAT model and its application in the Athabasca oil sands region.

Heavy metal contaminations in an aquatic environment is a serious issue since the exposure to toxic metals can cause a variety of public health problems. A watershed-scale model is a useful tool for predicting and assessing heavy metal fate and transport in both terrestrial and aquatic environments. In this study, we developed a simulation module for non-reactive heavy metals and incorporated it into the widely used Soil Water Assessment Tool (SWAT) model. The simulated processes in the developed model include heavy metal deposition, partitioning in soil and water, and transport by different pathways in both terrestrial and aquatic environments. Three-phase partitioning processes were considered in the module by simulating heavy metals portioning to dissolved organic carbon in the soil and stream. This developed module was used for watershed-scale simulation of heavy metal processes in the Muskeg River watershed (MRW) of the Athabasca oil sands region in western Canada for the first time. The daily streamflow and sediment load from 2015 to 2017 were first calibrated and validated. Subsequently, the daily Lead and Copper loads at the outlet station were used for heavy metal calibration and validation. The performances for the daily heavy metal loads simulation during the whole simulation period can be considered as "satisfactory" based on the recommended model performance criteria with the Nash-Sutcliffe efficiency as 0.41 and 0.71 for Pb and Cu loads, respectively. The simulation results indicate that the spring and summer are hot moments for heavy metal transport and the snowmelt in spring and rainfall runoff events in summer are the main driving forces for the metal transport in the MRW. We believe the developed model can be a useful tool for simulating the fate and transport of non-reative heavy metals at watershed scale and further used to assess management scenarios for mitigating heavy metal pollution in the Athabasca oil sands region.

Assessing climate change impacts on stream temperature in the Athabasca River Basin using SWAT equilibrium temperature model and its potential impacts on stream ecosystem.

Stream temperatures, which influence dynamics and distributions of the aquatic species and kinetics of biochemical reactions, are expected to be altered by the climate change. Therefore, predicting the impacts of climate change on stream temperature is helpful for integrated water resources management. In this study, our previously developed Soil and Water Assessment Tool (SWAT) equilibrium temperature model, which considers both the impacts of meteorological condition and hydrological processes, was used to assess the climate change impact on the stream temperature regimes in the Athabasca River Basin (ARB), a cold climate region watershed of western Canada. The streamflow and stream temperatures were calibrated and validated first in the baseline period, using multi-site observed data in the ARB. Then, climate change impact assessments were conducted based on three climate models under the Representative Concentration Pathways 4.6 and 8.5 scenarios. Results showed that warmer and wetter future condition would prevail in the ARB. As a result, streamflow in the basin would increase despite the projected increases in evapotranspiration due to warmer condition. On the basin scale, annual stream temperatures are expected to increase by 0.8 to 1.1 °C in mid-century and by 1.6 to 3.1 °C in late century. Moreover, the stream temperature changes showed a marked temporal pattern with the highest increases (2.0 to 7.4 °C) in summer. The increasing stream temperatures would affect water quality dynamics in the ARB by decreasing dissolved oxygen concentrations and increasing biochemical reaction rates in the streams. Such spatial-temporal changes in stream temperature regimes in future period would also affect aquatic species, thus require appropriate management measures to attenuate the impacts.

Integrating terrestrial and aquatic processes toward watershed scale modeling of dissolved organic carbon fluxes.

Dissolved organic carbon (DOC) is not only a critical component of global and regional carbon budgets, but also an important precursor for carcinogenic disinfection byproducts (DBP) generated during drinking water disinfection process. The lack of process based watershed scale model for carbon cycling has been a limiting factor impeding effective watershed management to control DOC fluxes to source waters. Here, we integrated terrestrial and aquatic carbon processes into the widely tested Soil and Water Assessment Tool (SWAT) watershed model to enable watershed-scale DOC modeling (referred to as SWAT-DOC hereafter). The modifications to SWAT mainly fall into two groups: (1) DOC production in soils and its transport to aquatic environment by different hydrologic processes, and (2) riverine transformation of DOC and their interactions with particular organic carbon (POC), inorganic carbon and algae (floating and bottom). We tested the new SWAT-DOC model in the Cannonsville watershed, which is part of the New York City (NYC) water supply system, using long-term DOC load data (from 1998 to 2012) derived from 1399 DOC samplings. The calibration and verification results indicate that SWAT-DOC achieved satisfactory performance for both streamflow and DOC at daily and monthly temporal scales. The parameter sensitivity analysis indicates that DOC loads in the Cannonsville watershed are controlled by the DOC production in soils and its transport in both terrestrial and aquatic environments. Further model uncertainty analysis indicates high uncertainties associated with peak DOC loads, which are attributed to underestimation of high streamflows. Therefore, future efforts to enhance SWAT-DOC to better represent runoff generation processes hold promise to further improve DOC load simulation. Overall, the wide use of SWAT and the satisfactory performance of SWAT-DOC make it a useful tool for DOC modeling and mitigation at the watershed scale.

Developing a non-point source P loss indicator in R and its parameter uncertainty assessment using GLUE: a case study in northern China.

Uncertainty analysis is an important prerequisite for model application. However, the existing phosphorus (P) loss indexes or indicators were rarely evaluated. This study applied generalized likelihood uncertainty estimation (GLUE) method to assess the uncertainty of parameters and modeling outputs of a non-point source (NPS) P indicator constructed in R language. And the influences of subjective choices of likelihood formulation and acceptability threshold of GLUE on model outputs were also detected. The results indicated the following. (1) Parameters RegR2, RegSDR2, PlossDPfer, PlossDPman, DPDR, and DPR were highly sensitive to overall TP simulation and their value ranges could be reduced by GLUE. (2) Nash efficiency likelihood (L1) seemed to present better ability in accentuating high likelihood value simulations than the exponential function (L2) did. (3) The combined likelihood integrating the criteria of multiple outputs acted better than single likelihood in model uncertainty assessment in terms of reducing the uncertainty band widths and assuring the fitting goodness of whole model outputs. (4) A value of 0.55 appeared to be a modest choice of threshold value to balance the interests between high modeling efficiency and high bracketing efficiency. Results of this study could provide (1) an option to conduct NPS modeling under one single computer platform, (2) important references to the parameter setting for NPS model development in similar regions, (3) useful suggestions for the application of GLUE method in studies with different emphases according to research interests, and (4) important insights into the watershed P management in similar regions.

Modeling nitrous oxide emissions from rough fescue grassland soils subjected to long-term grazing of different intensities using the Soil and Water Assessment Tool (SWAT).

Given the rising nitrous oxide (N2O) concentration in the atmosphere, it has become increasingly important to identify hot spots and hot moments of N2O emissions. With field measurements often failing to capture the spatiotemporal dynamics of N2O emissions, estimating them with modeling tools has become an attractive alternative. Therefore, we incorporated several semi-empirical equations to estimate N2O emissions with the Soil and Water Assessment Tool from nitrification and denitrification processes in soil. We then used the model to simulate soil moisture and the N2O flux from grassland soils subjected to long-term grazing (> 60 years) at different intensities in Alberta, Canada. Sensitivity analysis showed that parameters controlling the N2O flux from nitrification were most sensitive. On average, the accuracy of N2O emission simulations were found to be satisfactory, as indicated by the selected goodness-of-fit statistics and predictive uncertainty band, while the model simulated the soil moisture with slightly higher accuracy. As expected, emissions were higher from the plots with greater grazing intensity. Scenario analysis showed that the N2O emissions with the recommended fertilizer rate would dominate the emissions from the projected wetter and warmer future. The combined effects of fertilization and wetter and warmer climate scenarios would increase the current N2O emission levels by more than sixfold, which would be comparable to current emission levels from agricultural soils in similar regions.

Assessing climate change impacts on fresh water resources of the Athabasca River Basin, Canada.

Proper management of blue and green water resources is important for the sustainability of ecosystems and for the socio-economic development of river basins such as the Athabasca River Basin (ARB) in Canada. For this reason, quantifying climate change impacts on these water resources at a finer temporal and spatial scale is often necessary. In this study, we used a Soil and Water Assessment Tool (SWAT) to assess climate change impacts on fresh water resources, focusing explicitly on the impacts to both blue and green water. We used future climate data generated by the Canadian Center for Climate Modelling and Analysis Regional Climate Model (CanRCM4) with a spatial resolution of 0.22°×0.22° (~25km) for two emission scenarios (RCP 4.5 and 8.5). Results projected the climate of the ARB to be wetter by 21-34% and warmer by 2-5.4°C on an annual time scale. Consequently, the annual average blue and green water flow was projected to increase by 16-54% and 11-34%, respectively, depending on the region, future period, and emission scenario. Furthermore, the annual average green water storage at the boreal region was expected to increase by 30%, while the storage was projected to remain fairly stable or decrease in other regions, especially during the summer season. On average, the fresh water resources in the ARB are likely to increase in the future. However, evidence of temporal and spatial heterogeneity could pose many future challenges to water resource planners and managers.

Estimating internal P loading in a deep water reservoir of northern China using three different methods.

Much attention had been paid to reducing external loading of nutrients to improve water quality, while internal loading from sediment, which has been largely neglected, is also an important source for water eutrophication. The internal load in deep lakes or reservoirs is not easy to be detected and be quantified. In this study, three different methods (mass balance method, Fick's law, and regression equation) were combined to calculate the gross or/and net P release from sediment using limited data. Our results indicated that (1) the methods of mass balance and regression equation give similar results of sediment P release rate, with values of 0.889 and 0.902 mg m(2) d(-1), respectively, while the result of Fick's law was much lower (0.400 mg m(2) d(-1)); (2) Hot periods of sediment releasing were suggested to occur from March to April and from August to September, which correspond to periods of high risks of algae blooms. The remaining months of the year were shown as net nutrient retention; (3) for the whole region, Baihedam and Chaohekuqu were identified as zones with a higher possibility to release P from sediment. (4) P loading to the Miyun Reservoir was greater in the inflow than in the outflow, suggesting a portion of the inflow P load was retained in the water or sediment; hence, release of sediment P may continue to be a major source of phosphorus in the future.

[Net anthropogenic nitrogen input to Huaihe River Basin, China during 1990-2010].

Social economy in Huaihe River Basin had undergone enormous changes during 1990-2010. The grain yield had increased by 58%, from 64.14 million tons to 101.21 million tons, and the urbanization rate had increased by 22%, from 13% to 35%. Assessing the negative impacts of these high intensive human activities caused by rapid social development on terrestrial ecosystem would serve as a scientific basis for quantitative management of regional ecology. This paper estimated the spatial and temporal distribution of net anthropogenic nitrogen input (NANI) in Huaihe River Basin during 1990-2010. The results showed that there was an increasing trend in NANI in the period of 1990-2001, and after that this trend was slower. The NANI increased from approximately 17232 kg N · km(-2) · a(-1) in 1990 to a peak of 28771 kg N · km(-2) · a(-1) in 2003, and then declined to 26415 kg N · km(-2) · a(-1) in 2010. Chemical fertilizer and atmospheric deposition were the largest two sources of NANI, followed by food & feed import and biological nitrogen. Contributions from both chemical fertilizer and atmospheric deposition had been increasing continuously, respectively from 64% and 16% in 1990 to 77% and 19%. Our findings implied that the shift from fertilizer-supported agriculture and fossil fuel-supported industry to sci-tech lead economic development is urgently needed.

Variations in source apportionments of nutrient load among seasons and hydrological years in a semi-arid watershed: GWLF model results.

Quantifying source apportionments of nutrient load and their variations among seasons and hydrological years can provide useful information for watershed nutrient load reduction programs. There are large seasonal and inter-annual variations in nutrient loads and their sources in semi-arid watersheds that have a monsoon climate. The Generalized Watershed Loading Function model was used to simulate monthly nutrient loads from 2004 to 2011 in the Liu River watershed, Northern China. Model results were used to investigate nutrient load contributions from different sources, temporal variations of source apportionments and the differences in the behavior of total nitrogen (TN) and total phosphorus (TP). Examination of source apportionments for different seasons showed that point sources were the main source of TN and TP in the non-flood season, whereas contributions from diffuse sources, such as rural runoff, soil erosion, and urban areas, were much higher in the flood season. Furthermore, results for three typical hydrological years showed that the contribution ratios of nutrient loads from point sources increased as streamflow decreased, while contribution ratios from rural runoff and urban area increased as streamflow increased. Further, there were significant differences between TN and TP sources on different time scales. Our findings suggest that priority actions and management measures should be changed for different time periods and hydrological conditions, and that different strategies should be used to reduce loads of nitrogen and phosphorus effectively.