Assessment of Wastewater Reuse Effects on Nutrient Loads from Paddy Field Using Field-Scale Water Quality Model

CREAMS-PADDY, a modified version of the field-scale CREAMS model, simulates the hydrologic, sediment, and nutrient cycles in paddy fields. The CREAMS-PADDY model was applied to estimate the effects of using wastewater for irrigation on nutrient loads from paddy fields in Republic of Korea. The model was calibrated and validated using data from two rice paddy fields. The coefficient of determination between observed and simulated total nitrogen and total phosphorus were 0.92 and 0.57, respectively, for the calibration period and 0.84 and 0.73 for the validation period. Simulations showed that when using wastewater for irrigation, the total nitrogen loads increased by 210% and total phosphorus by 1,270% when compared with conventional water irrigation. The total nitrogen and total phosphorus concentration in the ponded water increased by 254 and 534%, respectively, when compared with conventional water irrigation. The effect of reducing N and P fertilizer application rates by 10, 30, and 50% on nutrient loads exiting a paddy field were also simulated using the validated CREAMS-PADDY model. These simulations indicated that total phosphorus loads from the paddy were reduced only slightly by reducing the fertilizer, while total nitrogen loads were reduced by as much as 8.8, 16.6, and 24.4% when N ferlitizer rates were reduced by 10, 30, and 50%, respectively. An erratum to this article can be found at     

Long-term trends and short-term variability of water quality in Skive Fjord, Denmark – nutrient load and mussels are the primary pressures and drivers that influence water quality

Nineteen years of monitoring data from the eutrophic Skive Fjord, Denmark were examined for linkages to external pressures and drivers, including nutrient inputs, meteorology and stocks of blue mussels. Linkages were examined by: 1) time-series analysis to document effects of nutrient reduction programs, 2) Pearson Rank correlations, 3) multivariate statistical analysis (PLS) to identify water quality variables with high predictability and their linkages to pressures, and 4) regression analysis to quantify relationships between pressures and water quality. Freshwater input, nitrogen load and phosphorus load showed decreasing trends through the period 1984–2002. The load reductions were only partially translated into trends in water quality: phosphorus decreased in most seasons, while total nitrogen decreased during winter and spring only. Phosphorus concentration had the highest predictability (explained by seasonal temperature variation) followed by transparency, silicate, tot-N, chlorophyll-a, primary productivity, phytoplankton diversity and phytoplankton turnover. The variation in pressures other than nutrient input confounded the relations between loads and water quality. High biomass of mussels led to reduced chlorophyll-a and increased transparency, while short-term variability in water column mixing led to changes in chlorophyll-a due to nutrient entrainment and coupling to benthic mussels.     

Implication of two in-stream processes in the fate of nutrients discharged by sewage system into a temporary river

The aim of this study was to better understand the fate of nutrients discharged by sewage treatment plants into an intermittent Mediterranean river, during a low-flow period. Many pollutants stored in the riverbed during the low-flow period can be transferred to the downstream environments during flood events. The study focused on two processes that affect the fate and the transport of nutrients, a physical process (retention in the riverbed sediments) and a biological process (denitrification). A spatial campaign was carried out during a low-flow period to characterize the nutrient contents of both water and sediments in the Vène River. The results showed high nutrient concentrations in the water column downstream of the treated wastewater disposal (up to 13,315 μg N/L for ammonium and 2,901 μg P/L for total phosphorus). Nutrient concentrations decreased rapidly downstream of the disposal whereas nutrient contents in the sediments increased (up to 1,898 and 784 μg/g for total phosphorus and Kjeldahl nitrogen, respectively). According to an in situ experiment using sediment boxes placed in the riverbed for 85 days, we estimated that the proportion of nutrients trapped in the sediments represents 25% (respectively 10%) of phosphorus (respectively nitrogen) loads lost from the water column. In parallel, laboratory tests indicated that denitrification occurred in the Vène River, and we estimated that denitrification likely coupled to nitrification processes during the 85 days of the experiment was significantly involved in the removal of nitrogen loads (up to 38%) from the water column and was greater than accumulation processes.     

Assessing the influence of nutrient reduction on water quality using a three-dimensional model: case study in a tidal estuarine system

A coupled three-dimensional hydrodynamic and water quality model has been developed and applied to the Danshuei River estuarine system and adjacent coastal sea. The water quality model considers various species of nitrogen, phosphorus, organic carbon, and phytoplankton as well as dissolved oxygen and is driven by a three-dimensional hydrodynamic model. The hydrodynamic and water quality models were validated with observations of water surface elevation, velocity, salinity distribution, and water quality parameters. Statistical error analysis shows that predictions of hydrodynamics, salinity, dissolved oxygen, and nutrients from the model simulation quantitatively agreed with the observed data. The validated model was then applied to predict water quality conditions as a result of a reduction in nutrient loadings based on different engineering strategies. The simulated results revealed that the dissolved oxygen concentration would increase significantly and would be higher than 2 mg/L in the main stream and in three tributaries to meet the minimum statutory requirement for dissolved oxygen. Active estuarine management focused on the reduction of anthropogenic nutrient loads is needed for improvement in water quality.     

A hydrologic network supporting spatially referenced regression modeling in the Chesapeake Bay Watershed

The U.S. Geological Survey has developed a methodology for statistically relating nutrient sources and land-surface characteristics to nutrient loads of streams. The methodology is referred to as SPAtially Referenced Regressions On Watershed attributes (SPARROW), and relates measured stream nutrient loads to nutrient sources using nonlinear statistical regression models. A spatially detailed digital hydrologic network of stream reaches, stream-reach characteristics such as mean streamflow, water velocity, reach length, and travel time, and their associated watersheds supports the regression models. This network serves as the primary framework for spatially referencing potential nutrient source information such as atmospheric deposition, septic systems, point-sources, land use, land cover, and agricultural sources and land-surface characteristics such as land use, land cover, average-annual precipitation and temperature, slope, and soil permeability. In the Chesapeake Bay watershed that covers parts of Delaware, Maryland, Pennsylvania, New York, Virginia, West Virginia, and Washington D.C., SPARROW was used to generate models estimating loads of total nitrogen and total phosphorus representing 1987 and 1992 land-surface conditions. The 1987 models used a hydrologic network derived from an enhanced version of the U.S. Environmental Protection Agency's digital River Reach File, and course resolution Digital Elevation Models (DEMs). A new hydrologic network was created to support the 1992 models by generating stream reaches representing surface-water pathways defined by flow direction and flow accumulation algorithms from higher resolution DEMs. On a reach-by-reach basis, stream reach characteristics essential to the modeling were transferred to the newly generated pathways or reaches from the enhanced River Reach File used to support the 1987 models. To complete the new network, watersheds for each reach were generated using the direction of surface-water flow derived from the DEMs. This network improves upon existing digital stream data by increasing the level of spatial detail and providing consistency between the reach locations and topography. The hydrologic network also aids in illustrating the spatial patterns of predicted nutrient loads and sources contributed locally to each stream, and the percentages of nutrient load that reach Chesapeake Bay.     

Water quality status and trends in agriculture-dominated headwaters; a national monitoring network for assessing the effectiveness of national and European manure legislation in The Netherlands

Large nutrient losses to groundwater and surface waters are a major drawback of the highly productive agricultural sector in The Netherlands. The resulting high nutrient concentrations in water resources threaten their ecological, industrial, and recreational functions. To mitigate eutrophication problems, legislation on nutrient application in agriculture was enforced in 1986 in The Netherlands. The objective of this study was to evaluate this manure policy by assessing the water quality status and trends in agriculture-dominated headwaters. We used datasets from 5 agricultural test catchments and from 167 existing monitoring locations in agricultural headwaters. Trend analysis for these locations showed a fast reduction of nutrient concentrations after the enforcement of the manure legislation (median slopes of ?0.55 mg/l per decade for total nitrogen (N-tot) and ?0.020 mg/l per decade for total phosphorus (P-tot)). Still, up to 76 % of the selected locations currently do not comply with either the environmental quality standards (EQSs) for nitrogen (N-tot) or phosphorus (P-tot). This indicates that further improvement of agricultural water quality is needed. We observed that weather-related variations in nutrient concentrations strongly influence the compliance testing results, both for individual locations and for the aggregated results at the national scale. Another important finding is that testing compliance for nutrients based on summer average concentrations may underestimate the agricultural impact on ecosystem health. The focus on summer concentrations does not account for the environmental impact of high winter loads from agricultural headwaters towards downstream water bodies.     

Modelling the Effects of Inflow Parameters on Lake Water Quality

A one-dimensional lake water quality model which includes water temperature, phytoplankton, phosphorus as phosphate, nitrogen as ammonia, nitrogen as nitrate and dissolved oxygen concentrations, previously calibrated for Lake Calhoun (USA) is applied to Uokiri Lake (Japan) for the year 1994. The model simulated phytoplankton and nutrient concentrations in the lake from July to November. Most of the water quality parameters are found to be the same as for Lake Calhoun. To predict probable lake water quality deterioration from algal blooming due to increased nutrient influx from river inflow, the model was run for several inflow water conditions. Effects of inflow nutrient concentration, inflow volume, inflow water temperatures are presented separately. The effect of each factor is considered in isolation although in reality more than one factor can change simultaneously. From the results it is clear that inflow nutrient concentration, inflow volume and inflow water temperature show very regular and reasonable impacts on lake water quality.     

Estimation of riverine loads of nitrogen and phosphorus to the Baltic Sea, 1970–1993

This article presents the results of the first critical examination of time series of riverine nutrient-load data for the entire Baltic Sea drainage area. Water quality data collected by or for the different national environmental agencies were compiled and analysed statistically to identify and remove inconsistent or obviously incorrect observations. Moreover, sampling tours were undertaken to acquire additional information about the present nutrient concentrations in the largest rivers in the study area. Gaps in the time series of approved data were then filled in by employing statistical interpolation and extrapolation methods. Thereafter, the concentration and runoff data were combined to obtain estimates of monthly nutrient loads for the time period 1970–93. The results of the calculations showed that although there had been substantial changes in land use, atmospheric deposition and wastewater treatment in many parts of the study area, the total riverine loads of nitrogen (N) and phosphorus (P) to the Baltic Sea have been fairly constant since 1980, and most likely also since 1970. Moreover, the interannual variation was clearly correlated to the runoff. The mean annual loads for the time period 1980–93 were found to be about 825 000 tonnes N and 41 000 tonnes P, respectively. This implies that (i) several other investigators have strongly underestimated the riverine loads of nutrients, especially the nitrogen, and that (ii) the riverine loads by far exceed the input to the Baltic Sea from other sources, {i.e.} atmospheric deposition, direct emissions from cities and industries along the Baltic Sea coast and nitrogen fixation by marine algae.     

The Pollution Load by Nitrogen and Phosphorus IN the Jadro River

The objective of the investigations of the Jadro River, located in Croatia, was to estimate the nitrogen and phosphorus loads in the Jadro River spring and its streamflow by calculating the load in kg/day or tons/year and to compare this with the load for the maximum allowed concentrations (MAC) for drinking water (Official Bulletin, No 46/94) expressed in kg/day or tons/year.Daily pollution loads at the Jadro River spring for total N ranged from 0 to 304 kg, for NH₃-N from 0 to 38 kg, for NO₃-N from 0–1321 kg and for PO₄-P from 0–92 kg in the period from September 1993 to September 2003. When compared with MAC loads the results prove that the Jadro River spring is not polluted by nitrogen compounds and phosphorus.The average annual load for total N ranged from 10 to 33 t, for NH₃-N from 0.25 to 5.15 t, for NO₃-N from 40 to 190 t, and for PO₄-P from 0.3 to 11.5 t. The nitrogen compounds and phosphorus loads vary from one year to another without any constant decreasing or increasing trends. The annual average loads compared with the average annual MAC loads (especially for NH₃-N and PO₄-P) show that there were no threats of constant pollution of the spring.The loads for total N and PO₄-P along the Jadro River flow from the spring to the fishpond entrance were monitored over a five year period from September 1999 to September 2004. The results show no regularities. The highest annual total nitrogen load of 45 t was recorded at the entrance to the fishpond during the 2002/2003 period. The highest annual PO₄-P load of 10 t was measured at the Vidović Bridge during the 2003/2004 period; however, the concentrations of N and P did not exceed the MAC concentrations which are prescribed for drinking water.According to the investigation results of the daily and average annual loads compared with MAC loads for drinking water, it can be concluded that the Jadro River spring and its streamflow are not polluted by nitrogen and phosphorus.     

望虞河西岸河网区沉积物营养盐和重金属污染时空分布特征

以无锡市望虞河西岸河网区为研究区域,于2018年12月—2019年9月分冬、春、夏、秋4个季节采集了上覆水和沉积物样本,研究氮磷营养因子在沉积物-上覆水界面的释放规律,探究主要水质指标、沉积物氮磷和重金属含量的时空分布特征。结果表明:研究区域内的水质状况不佳,主要是氮含量超标,污染严重的点位集中在人类活动较为密集的市区和工业区;受上覆水氮磷浓度、溶解氧浓度等变化的影响,沉积物营养盐含量季节变化明显,污染状况较为严重,75%的点位属于中度和重度污染;沉积物重金属含量季节变化不显著,部分重金属的空间分布特征相似,污染程度排序为Zn>Cd>Cu>Pb>Ni>Cr>As。     

Sensitivity of Simulated Conservation Practice Effectiveness to Representation of Field and In-Stream Processes in the Little River Watershed

Evaluating the effectiveness of conservation practices (CPs) is an important step to achieving efficient and successful water quality management. Watershed-scale simulation models can provide useful and convenient tools for this evaluation, but simulated conservation practice effectiveness should be responsive to parameter values used to represent the practices in the modeling. The objectives of this study were to (1) assess the impacts of a set of conservation practices on hydrology and water quality of a watershed and (2) evaluate the sensitivity of Soil and Water Assessment Tool (SWAT) modeling outputs and simulated conservation practice effectiveness to parameters. The modeling study was conducted in an agricultural watershed, the subwatershed K (16.9 km²) of the Little River Experimental watershed located in the South Atlantic Coastal Plain of the USA. Sensitivity analysis showed that hydrologic response unit (HRU) and watershed-scale simulations for water quality were most sensitive to CN and FILTERW parameters. Load reduction rates as a function of increased aerial coverage of the conservation practices were greatest for total phosphorus (TP), followed by sediment and total nitrogen (TN). The results indicated that conservation practices would have a limited impact on stream flow volume but could have a significant impact on sediment and TP loads within this region. Watershed-scale TN and TP loads were also sensitive to an in-stream nutrient transformation process represented using the QUAL2E algorithm in SWAT. The study clearly demonstrated the most sensitive model parameters and the optimal conservation practices for this watershed.     

An appraisal of changes in seasonal water quality during passage through a shallow reservoir in Western Poland

Following restoration changes in Antoninek Reservoir physico-chemical and biological processes in the water column and bottom sediments were measured to outline mechanisms of changes in nitrogen, phosphorus and organic matter concentrations during water flow through this reservoir. Intensive mineralisation of organic matter in the shallow sediments stimulated primary production and influenced increasing ammonia and nitrite nitrogen concentrations. Two main factors affected concentrations of phosphorus: (1) its presence in the external loads of river waters entering the reservoir, more important in the colder seasons as the water discharge was higher and (2) from the internal loads coming from bottom sediments. The quality of the river water during its flow through this reservoir improved for most parameters and seasons. However, concentrations of nutrients were still high in waters flowing out from the reservoir and in some months they were higher in the outflow than in waters entering the reservoir.     

Using process-based models to filter out natural variability in observed concentrations of nitrogen and phosphorus in river water

Advances in process-based modelling of loads of nitrogen and phosphorus carried by rivers have created new possibilities to interpret time series of water quality data. We examined how model runs with constant anthropogenic forcing can be used to estimate and filter out weather-driven variation in observational data and, thereby, draw attention to other features of such data. An assessment of measured and modelled nutrient concentrations at the outlets of 45 Swedish rivers provided promising results for total nitrogen. In particular, joint analyses of observational data and outputs from the catchment model S-HYPE strengthened the evidence that downward trends in nitrogen were due to mitigation measures in agriculture. Evaluation of modelled and observed total phosphorus concentrations revealed considerable bias in the collection or chemical analysis of water samples and also identified weaknesses in the model outputs. Together, our results highlight the need for more efficient two-way communication between environmental modelling and monitoring.     

Determining storm sampling requirements for improving precision of annual load estimates of nutrients from a small forested watershed

This study sought to determine the lowest number of storm events required for adequate estimation of annual nutrient loads from a forested watershed using the regression equation between cumulative load (∑L) and cumulative stream discharge (∑Q). Hydrological surveys were conducted for 4 years, and stream water was sampled sequentially at 15-60-min intervals during 24 h in 20 events, as well as weekly in a small forested watershed. The bootstrap sampling technique was used to determine the regression (∑L-∑Q) equations of dissolved nitrogen (DN) and phosphorus (DP), particulate nitrogen (PN) and phosphorus (PP), dissolved inorganic nitrogen (DIN), and suspended solid (SS) for each dataset of ∑L and ∑Q. For dissolved nutrients (DN, DP, DIN), the coefficient of variance (CV) in 100 replicates of 4-year average annual load estimates was below 20% with datasets composed of five storm events. For particulate nutrients (PN, PP, SS), the CV exceeded 20%, even with datasets composed of more than ten storm events. The differences in the number of storm events required for precise load estimates between dissolved and particulate nutrients were attributed to the goodness of fit of the ∑L-∑Q equations. Bootstrap simulation based on flow-stratified sampling resulted in fewer storm events than the simulation based on random sampling and showed that only three storm events were required to give a CV below 20% for dissolved nutrients. These results indicate that a sampling design considering discharge levels reduces the frequency of laborious chemical analyses of water samples required throughout the year.     

Using Geographically Isolated Loading Scenarios to Analyze Nitrogen and Phosphorus Exchanges and Explore Tailored Nutrient Control Strategies for Efficient Management

A set of geographically isolated differential nitrogen (N) and phosphorus (P) load model scenarios from major Chesapeake basins provides information on the relative impact of nutrient loads on primary production and dissolved oxygen in the Chesapeake Bay. Model results show the relationships of deep water dissolved oxygen with nutrient limitation-related algal blooms, organic carbon loads from the watershed, estuarine circulation, nutrient cycling, and nutrient diagenesis. The combined effect of changes in load from multiple basins is additive for changes in both chlorophyll-a and deep water dissolved oxygen concentrations. Management of both N and P are required in the Chesapeake watershed and tidal waters to achieve water quality standards, but overall efficiencies could be gained with strategies that place greater emphasis on P control in the upper Bay and greater emphasis on N control in the lower Bay. The areas of the Bay with the highest degree of dissolved oxygen degradation that generally drive management decisions are mostly P-limited and are significantly influenced by the load from the upper Bay’s basins. Reducing P from the upper Bay’s basins will intensify P limitation and would allow an increase in N of about six times the weight of P reduction. Combining the relative nutrient reduction effectiveness with the relative control cost information could improve management efficiency and provide benefits at a lower cost. This article describes initial steps that can be taken to examine the benefits from N-P exchanges.     

Load estimation and source apportionment of nonpoint source nitrogen and phosphorus based on integrated application of SLURP model, ECM, and RUSLE: a case study in the Jinjiang River, China

The nonpoint source (NPS) pollution is difficult to manage and control due to its complicated generation and formation. Load estimation and source apportionment are an important and necessary process for efficient NPS control. Here, an integrated application of semi-distributed land use-based runoff process (SLURP) model, export coefficients model (ECM), and revise universal soil loss equation (RUSLE) for the load estimation and source apportionment of nitrogen and phosphorus was proposed. The Jinjiang River (China) was chosen for the evaluation of the method proposed here. The chosen watershed was divided into 27 subbasins. After which, the SLURP model was used to calculate land use runoff and to estimate loads of dissolved nitrogen and phosphorus, and ECM was applied to estimate dissolved loads from livestock and rural domestic sewage. Next, the RUSLE was employed for load estimation of adsorbed nitrogen and phosphorus. The results showed that the 12,029.06 t?a^?1 pollution loads of total NPS nitrogen (TN) mainly originated from dissolved nitrogen (96.24 %). The major sources of TN were land use runoff, which accounted for 45.97 % of the total, followed by livestock (32.43 %) and rural domestic sewage (17.83 %). For total NPS phosphorous (TP), its pollution loads were 570.82 t?a^?1 and made up of dissolved and adsorbed phosphorous with 66.29 and 33.71 % respectively. Soil erosion, land use runoff, rural domestic sewage, and livestock were the main sources of phosphorus with contribution ratios of 33.71, 45.73, 14.32, and 6.24 % respectively. Therefore, land use runoff, livestock, and soil erosion were identified as the main pollution sources to influence loads of NPS nitrogen and phosphorus in the Jinjiang River and should be controlled first. The method developed here provided a helpful guideline for conducting NPS pollution management in similar watershed.     

Nutrient spatial pattern of the upstream,mainstream and tributaries of the Three Gorges Reservoir in China

A comprehensive monitoring program was conducted to investigate the nutrient spatial pattern in the mainstream of the Yangtze River from the Baihetan Dam down to the Three Gorges Dam located at the upper region of the Yangtze River in China. Samples were taken from 33 different sites from July 30 to August 19, 2011. The nutrient patterns of the three representative tributaries of the Three Gorges Reservoir (TGR)—the Modao, the Daning, and the Xiangxi Rivers—were also investigated. The results show that the mainstream of the TGR has a higher concentration of nitrogen and a lower concentration of phosphorus than that of the upper mainstream before the TGR. Moreover, it was found that nitrate-nitrogen (NO₃-N) is the main nitrogen component, while particulate phosphorus predominates the total phosphorus (TP). It was found that the three representative tributaries of the TGR have lower total nitrogen (TN) concentrations compared to the corresponding sections of the mainstream TGR. Based on the nutrient spatial pattern, the nutrient flux was calculated. The total fluxes of TN, NO₃-N, TP, and orthophosphate (PO₄-P) from the upstream reach into the TGR are 2,155.06, 1,674.97, 212.98, and 83.42 t day^?1, respectively. The amount of nutrients imported from the TGR into its tributaries is more than the amount exported. It was determined that the Xiangxi River has the largest net rate of imported nitrogen at 7.66 t day^?1, whereas the Daning River has the largest net rate of imported phosphorus at 1.75 t day^?1. In addition, compared with the nutrients imported from the TGR into its tributaries, the nutrient flux from the upstream reach into the TGR contributes approximately less than 3 %.     

Evaluation of Point and Diffuse Sources of Nutrients in a River Basin on Base of Monitoring Data

The methodology of materials accounting is presented and applied to developing nutrient balance (nitrogen and phosphorus) in a river basin. The method is based on the balance principle: inputs and outputs of each nitrogen and phosphorus related sub-systems were balanced. The application of the methodology strategies was illustrated by means of a case study of the Krka river, Slovenia. Different pathways of emission to surface waters were taken into account: WWTP discharges, direct discharges, erosion/runoff and baseflow. Total annual emission into the river Krka was estimated to be 362 tonnes N/year and 73.3 tonnes P/year. The main sources of nitrogen are diffuse sources, emitted via baseflow (52%). Other important sources are effluents from WWTP, which account for 36% of total emissions. Other sources like erosion and direct discharges to surface water (animal manure, industry, households) are of lower magnitude. Erosion is main source of phosphorus emission (55% of total emission), WWTP effluents account for 37% of total emission, while other sources are less important. Besides reduction of point sources by means of wastewater collection and implementation of nutrient removal technology, managing agricultural nitrogen and phosphorus to protect water quality should become a major challenge in the Krka river basin.     

洮滆水系湖库富营养化生态风险的特点与比较

茅东水库、长荡湖、涌湖、太湖竺山湾是洮滆水系从上游到下游排列的4大典型湖库,2008年的监测分析表明,氮、磷是该水系湖库富营养化的主要污染因子,并沿流域呈加剧趋势,上下游TP质量浓度为0.081～0.296 mg/L,差异小,而TN质量浓度为0.314～5.67 mg/L,差异大,长荡湖到涌湖是洮滆水系首要污染物TN快...     

Hydrology, suspended sediment dynamics and nutrient loading in Lake Takkobu, a degrading lake ecosystem in Kushiro Mire, northern Japan

Suspended sediment and nutrient loadings from agricultural watersheds have lead to habitat degradation in Lake Takkobu. To examine their relationships with land-use activities, we monitored sediment, nutrient and water discharges into the lake for a 1-year sampling period. The Takkobu River contributed the largest portion of the annual water discharge into the lake, compared with the other tributaries. During dry conditions, lake water flowed into the Kushiro River, and conversely during flooding, Kushiro River water flowed into the lake. Inflows from the Kushiro River had a high proportion of inorganic matter, with high concentrations of total nitrogen and total phosphorus, attributed to agricultural land-use development and stream channelization practiced since the 1960s in the Kushiro Mire. Nutrient loadings from these two rivers were significantly higher during flooding than in dry conditions. However, there was no clear correlation between river discharge and nutrient concentrations. Since land-use activities in the Kushiro River and Takkobu River watersheds were concentrated near rivers, nutrients easily entered the drainage system under low flow conditions. In contrast, water discharged from small, forest-dominated watersheds contained a low proportion of inorganic matter, and low nutrient concentrations. The suspended sediment delivered to the lake during the sample period was estimated as approximately 607 tons, while the total nitrogen and total phosphorus inflows were about 10,466 and 1,433 kg, respectively. Suspended sediment input into the lake was 65%, and total nitrogen and total phosphorus were 40% and 48%, respectively, being delivered by the Kushiro River.