Nonpoint Source Pollution Assessment of Wujiang RiverWatershed in Guizhou Province,SW China

The amount of pollution from nonpoint sources flowing in the streams of the Wujiang River watershed in Guizhou Province, SW China, is estimated by a geographic information system (GIS)-based method using rainfall, surface runoff and land use data. A grid of cells of 100 m in size is laid over the landscape. For each cell, mean annual surface runoff is estimated from rainfall and percent land use, and expected pollutant concentration is estimated from land use. The product of surface runoff and concentration gives expected pollutant loading from that cell. These loadings are accumulated going downstream to give the expected annual pollutant loadings in streams and rivers. By dividing these accumulated loadings by the similarly accumulated mean annual surface runoff, the expected pollutant concentration from nonpoint sources is determined for each location in a stream or river. Observed pollutant concentrations in the watershed are averaged at each sample point and compared to the expected concentrations at the same locations determined from the grid cell model. In general, annual nonpoint source nutrient loadings in the Wujiang River watershed are seen to be predominantly from the agricultural and meadow areas. The total annual loadings through the outlet of the watershed are 40,309 and 2,607 tons for total nitrogen (TN) and total phosphorus (TP), respectively.     

The importance of small urbanized watersheds to pollutant loading in a large oligotrophic subalpine lake of the western USA

Urban land use has been implicated as a major contributor of nonpoint source pollution in aquatic systems. Through increased nonpoint delivery of pollutants, including constituents found in stormwater, Lake Tahoe is undergoing a marked decline in its transparency, primarily due to increasing production of algae from enhanced nutrient loading and delivery of fine particles to the lake from the watershed. In response to these findings, a regional restoration effort is underway to improve basin watersheds and the water quality in Lake Tahoe. In this study, stormwater autosamplers were used to collect flow-weighted composite samples that characterized event mean concentrations for event and nonevent conditions within a small, urbanized watershed in the Tahoe basin. An event-specified constant-concentration water quality model was then applied to the event mean concentration and continuous streamflow data to estimate pollutant loads for nitrate, nitrite, ammonia, orthophosphate, and suspended sediment. These data were compared with previously reported load estimates from 10 primary monitored streams in larger watersheds of the Tahoe basin. Results from a linear regression analysis demonstrate strong and significant relationships between watershed impervious area and pollutant loadings from Lake Tahoe watersheds. These small, urbanized watersheds and intervening zones, which only comprise 10 % of the total Lake Tahoe drainage area, include a significant portion of the total Lake Tahoe impervious area. The findings of this study suggest that small, urbanized watersheds and intervening zones are disproportionately important contributors of nonpoint source pollution, including nutrients and suspended particles.     

Phosphorus Load Reduction Goals for Feitsui Reservoir Watershed, Taiwan

The present paper describes an effort for developing the total maximum daily load (TMDL) for phosphorus and a load reduction strategy for the Feitsui Reservoir in Northern Taiwan. BASINS model was employed to estimate watershed pollutant loads from nonpoint sources (NPS) in the Feitsui Reservoir watershed. The BASINS model was calibrated using field data collected during a 2-year sampling period and then used to compute watershed pollutant loadings into the Feitsui Reservoir. The simulated results indicate that the average annual total phosphorus (TP) loading into the reservoir is 18,910 kg/year, which consists of non-point source loading of 16,003 kg/year, and point source loading of 2,907 kg/year. The Vollenweider mass balance model was used next to determine the degree of eutrophication under current pollutant loading and the load reduction needed to keep the reservoir from being eutrophic. It was estimated that Feitsui Reservoir can becoming of the oligotrophic state if the average annual TP loading is reduced by 37% or more. The results provide the basis on which an integrated control action plan for both point and nonpoint sources of pollution in the watershed can be developed.     

Influence of different nitrate–N monitoring strategies on load estimation as a base for model calibration and evaluation

Model-based predictions of the impact of land management practices on nutrient loading require measured nutrient flux data for model calibration and evaluation. Consequently, uncertainties in the monitoring data resulting from sample collection and load estimation methods influence the calibration, and thus, the parameter settings that affect the modeling results. To investigate this influence, we compared three different time-based sampling strategies and four different load estimation methods for model calibration and compared the results. For our study, we used the river basin model Soil and Water Assessment Tool on the intensively managed loess-dominated Parthe watershed (315 km²) in Central Germany. The results show that nitrate–N load estimations differ considerably depending on sampling strategy, load estimation method, and period of interest. Within our study period, the annual nitrate–N load estimation values for the daily composite data set have the lowest ranges (between 9.8% and 15.7% maximum deviations related to the mean value of all applied methods). By contrast, annual estimation results for the submonthly and the monthly data set vary in greater ranges (between 24.9% and 67.7%). To show differences between the sampling strategies, we calculated the percentage deviation of mean load estimations of submonthly and monthly data sets as related to the mean estimation value of the composite data set. For nitrate–N, the maximum deviation is 64.5% for the submonthly data set in the year 2000. We used average monthly nitrate–N loads of the daily composite data set to calibrate the model to achieve satisfactory simulation results [Nash–Sutcliffe efficiency (NSE) 0.52]. Using the same parameter settings with submonthly and monthly data set, the NSE dropped to 0.42 and 0.31, respectively. Considering the different results from the monitoring strategy and the load estimation method, we recommend both the implementation of optimized monitoring programs and the use of multiple load estimation methods to improve water quality characterization and provide appropriate model calibration and evaluation data.     

Screening Assessment of Dissolved Zinc from Inactive Mines in a Mountain Watershed

A screening-level assessment of dissolved zinc from inactive and abandoned metal mines in the Cement Creek Basin was performed. The basin is part of the Upper Animas River Basin in the San Juan Mountains of southwestern Colorado. Stream discharge and dissolved zinc concentrations were measured at 49 stations below nonpoint sources including tailings and waste rock, point sources such as adits, and background areas. One measurement was made at a station during three flow events: storm runoff, peak snowmelt runoff, and baseflow. The highest concentrations occurred in the upper part of the basin immediately downstream from nonpoint and point source discharges, especially in the North Fork of Cement Creek. The mean concentration in Cement Creek was highest during baseflow (1350 g l-1) and lowest during snowmelt (796 g l-1). Most exceedances of national acute and chronic criteria for brown trout were chronic criteria exceedances in the upper part of the basin. Subareas with the greatest extent of nonpoint source areas in the upper part of the watershed, especially those contributing to Upper Cement Creek and the North Fork of Cement Creek, generally had the highest loadings and unit area loadings. The greatest loadings from all subareas to Cement Creek occurred during snowmelt (219 000 g day-1 and 17 400 kg for the snowmelt season). The highest unit area loadings from all subareas also occurred during snowmelt (190 g ha-1 day-1 and 15 147 g ha-1 for the snowmelt season). Loadings from subareas with extensive nonpoint source areas were always much greater than those from point sources and background areas.     

A Spatially Explicit Model for Estimating Annual Average Loads of Nonpoint Source Nutrient at the Watershed Scale

The overloaded nonpoint source (NPS) nutrients in upper streams always result in the nutrient enrichment at lakes and estuaries downstream. As NPS pollution has become a serious environmental concern in watershed management, the information about nutrient output distribution across a watershed has been critical in the designing of regional development policies. But existing watershed evaluation models often encounter difficulties in application because of their complicated structures and strict requirements for the input data. In this paper, a spatially explicit and process-based model, Integrated Grid’s Exporting and Delivery model, was introduced to estimate annual in-stream nutrient levels. Each grid cell in this model was regarded as having potentials of both exporting new nutrients and trapping nutrients passing by. The combined nutrient dynamics of a grid is mainly determined by the grid’s features in land use/land cover, soil drainage, and geomorphology. This simple-concept model was tested at some basins in north Georgia in the USA. Stations in one basin were used to calibrate the model. Then an external validation was employed by applying the calibrated model to stations in the other neighbor basins. Model evaluation statistics implied the model’s validity and good performance in estimating the annual NPS nutrients’ fluxes at the watershed scale. This study also provides a promising prospect that in-stream annual nutrient loads can be accurately estimated from a few public available datasets.     

Quantitative assessment of eutrophication: A scoring system for characterising water quality in coastal marine ecosystems

A scoring system based on nutrient concentrations was developed to assess coastal water quality according to the trophic level. Three nutrient data sets from eutrophic, mesotrophic and oligotrophic waters were used as the reference information for setting up a semi quantitative water quality scale (from 0 to 5) to express different nutrient loadings. The validity and sensitivity of the method was applied to a number of stations spaced out along the coastal area of Rhodes (Greece). A score for each nutrient/sampling site was calculated and the scorecard formed, was the data matrix used for numerical classification of the stations. The results showed (a) good discrimination between cutrophic, mesotrophic and oligotrophic waters (b) nitrate among the nutrients showed the maximum sensitivity in characterising pollution levels. The reference data sets used for assessing eutrophication levels ensured the objectivity of the method. The proposed method is described step-by-step and it is suggested that the method can be further adapted to describe other forms of pollution becoming a useful quantitative technique in coastal management practices.     

Non-intrusive characterization methods for wastewater-affected groundwater plumes discharging to an alpine lake

Streams and lakes in rocky environments are especially susceptible to nutrient loading from wastewater-affected groundwater plumes. However, the use of invasive techniques such as drilling wells, installing piezometers or seepage meters, to detect and characterize these plumes can be prohibitive. In this work, we report on the use of four non-intrusive methods for this purpose at a site in the Rocky Mountains. The methods included non-invasive geophysical surveys of subsurface electrical conductivity (EC), in-situ EC measurement of discharging groundwater at the lake-sediment interface, shoreline water sampling and nutrient analysis, and shoreline periphyton sampling and analysis of biomass and taxa relative abundance. The geophysical surveys were able to detect and delineate two high-EC plumes, with capacitively coupled ERI (OhmMapper) providing detailed two-dimensional images. In situ measurements at the suspected discharge locations confirmed the presence of high-EC water in the two plumes and corroborated their spatial extent. The nutrient and periphyton results showed that only one of the two high-EC plumes posed a current eutrophication threat, with elevated nitrogen and phosphorus levels, high localized periphyton biomass and major shifts in taxonomic composition to taxa that are commonly associated with anthropogenic nutrient loading. This study highlights the need to use non-intrusive methods in combination, with geophysical and water EC-based methods used for initial detection of wastewater-affected groundwater plumes, and nutrient or periphyton sampling used to characterize their ecological effects.     

Spatial properties of nutrient concentrations in the Baltic Sea

Using existing monitoring data, the present study attempts to characterize spatial patterns of surface water nutrients in the Baltic Sea. The analysis was made for two different spatial scales, differences between and within sub-basins.Non-parametric methods were chosen to reduce problems with the distributional properties of the data. By dividing the data into four seasons care was taken as to seasonality. To avoid bias due to different laboratory analysis, only data from Nordic countries were used, whose results did not apparently depart from each other. Bias due to different sampling frequency was another problem. This was reduced by using only the last observation in each season for every station and year.The results suggested differences in the nutrient concentrations between basins. The two northernmost basins (Bothnian Bay and Bothnian Sea) had lower phosphate concentrations and higher silicate concentrations compared to the rest of the Sea. Bothnian Bay and the Gulf of Finland had higher nitrate concentrations.The concentration structure within basins was studied using transects in both latitudinal and longitudinal direction. A gradient for phosphate and nitrate was found in the Gulf of Finland, with lower concentrations at the mouth. The Bothnian Sea showed lower concentrations in the middle of the basin compared to the coasts. The highest concentrations of phosphate were found close to the Finnish coast and for silicate the highest concentrations were located near the Swedish coast. It was not possible to study variation in the west-east direction within the Baltic proper, due to data shortage. For the transects in the north-south direction no differences were detected for nitrate and silicate. Phosphate gave one significant result during autumn for the transect in the eastern part of the Baltic proper.This study revealed several problems associated with the data available. Uneven sampling in space and time put severe constraints on the study. A better design of the monitoring program is suggested, where spatial properties are considered explicitly. Revision of the program in this direction is needed for reasonable calculation of total amounts and concentrations representative for a basin.     

Performance evaluation of a full-scale natural treatment system for nonpoint source and point source pollution removal

This study presents a full-scale performance of a natural treatment system (NTS) facility in Taiwan with nearly 2 years of observations. The study site, composed of several treatment ponds in series, was designed primarily to reduce polluted stormwater runoff from tea gardens and partially to untreated domestic wastewater from nearby villages. Thus, both nonpoint source and point source pollution are treated in this system. From 28 field samplings in 2006-2007, the NTS site shows satisfactory treatment performance and the effluent water quality is significantly improved. Seven of the 28 sampling events are storm events (nonpoint source pollution) and the remainder are from regular monitoring (point source pollution). The average volume of influent and effluent is 533 CMD and 196 CMD, respectively. In order to determine the removal efficiency, several assessment measures are employed in an attempt to obtain unbiased conclusions. They are removal rate (RR), efficiency rate (ER), summation of loads (SOL), flux rate (FR), and effluent probability method (EPM). The average percent removal efficiency of NH(3)-N is 53.5-75.2% and of TP is 59.0-84.7%, in which the highest result is calculated by SOL method and the lowest rate is obtained from RR. In FR evaluation, larger treatment capacity for NH(3)-N than for TP is provided in the site and the average FR is respectively 0.230 g/m(2) day and 0.017 g/m(2) day. Of the methods examined, EPM is the only method capable of illustrating data distribution. Finally, recommendations on the usefulness of these measures are summarized to facilitate the understandings of NTS performance evaluations.     

Application of multivariate statistical techniques in the assessment of water quality in the Southwest New Territories and Kowloon, Hong Kong

The application of different multivariate statistical techniques for the interpretation of a complex data matrix obtained during 2000?C2007 from the watercourses in the Southwest New Territories and Kowloon, Hong Kong was presented in this study. The data set consisted of the analytical results of 23 parameters measured monthly at 16 different sampling sites. Hierarchical cluster analysis grouped the 12 months into two periods and the 16 sampling sites into three groups based on similarity in water quality characteristics. Discriminant analysis (DA) provided better results both temporally and spatially. DA also offered an important data reduction as it only used four parameters for temporal analysis, affording 84.2% correct assignations, and eight parameters for spatial analysis, affording 96.1% correct assignations. Principal component analysis/factor analysis identified four latent factors standing for organic pollution, industrial pollution, nonpoint pollution, and fecal pollution, respectively. KN1, KN4, KN5, and KN7 were greatly affected by organic pollution, industrial pollution, and nonpoint pollution. The main pollution sources of TN1 and TN2 were organic pollution and nonpoint pollution, respectively. Industrial pollution had high effect on TN3, TN4, TN5, and TN6.     

Application of integrated GIS and multimedia modeling on NPS pollution evaluation

In Taiwan, nonpoint source (NPS) pollution is one of the major causes of the impairment of surface waters. I-Liao Creek, located in southern Taiwan, flows approximately 90 km and drains toward the Kaoping River. Field investigation results indicate that NPS pollution from agricultural activities is one of the main water pollution sources in the I-Liao Creek Basin. Assessing the potential of NPS pollution to assist in the planning of best management practice (BMP) is significant for improving pollution prevention and control in the I-Liao Creek Basin. In this study, land use identification in the I-Liao Creek Basin was performed by properly integrating the skills of geographic information system (GIS) and global positioning system (GPS). In this analysis, 35 types of land use patterns in the watershed area of the basin are classified with the aid of Erdas Imagine process system and ArcView GIS system. Results indicate that betel palm farms, orchard farms, and tea gardens dominate the farmland areas in the basin, and are scattered around on both sides of the river corridor. An integrated watershed management model (IWMM) was applied for simulating the water quality and evaluating NPS pollutant loads to the I-Liao Creek. The model was calibrated and verified with collected water quality and soil data, and was used to investigate potential NPS pollution management plans. Simulated results indicate that NPS pollution has significant contributions to the nutrient loads to the I-Liao Creek during the wet season. Results also reveal that NPS pollution plays an important role in the deterioration of downstream water quality and caused significant increase in nutrient loads into the basin's water bodies. Simulated results show that source control, land use management, and grassy buffer strip are applicable and feasible BMPs for NPS nutrient loads reduction. GIS system is an important method for land use identification and waste load estimation in the basin. Linking the information of land utilization with the NPS pollution simulation model may further provide essential information of potential NPS pollution for all subregions in the river basin. Results and experience obtained from this study will be helpful in designing the watershed management and NPS pollution control strategies for other similar river basins.     

A comparison of X-ray fluorescence and wet chemical analysis of air filter samples from a scrap lead smelting operation

Personal and area air samples were taken at a scrap lead smelter operation in a bullet manufacturing facility. Samples were taken using the 37-mm styrene-acrylonitrile closed-face filter cassette (CFC, the current US standard device for lead sampling), the 37-mm GSP or "cone" sampler, the 25-mm Institute of Occupational Medicine (IOM) inhalable sampler, and the 25-mm Button sampler (developed at the University of Cincinnati). Polyvinylchloride filters were used for sampling. The filters were pre- and post-weighed, and analyzed for lead content using a field-portable X-ray fluorescence (XRF) analyzer. The filters were then extracted with dilute nitric acid in an ultrasonic extraction bath and the solutions were analyzed by inductively coupled plasma optical emission spectroscopy. The 25-mm filters were analyzed using a single XRF reading, while three readings on different parts of the filter were taken from the 37-mm filters. The single reading from the 25-mm filters was adjusted for the nominal area of the filter to obtain the mass loading, while the three readings from the 37-mm filters were inserted into two different algorithms for calculating the mass loadings, and the algorithms were compared. The IOM sampler was designed for material collected in the body of the sampler to be part of the collected sample as well as that on the filter. Therefore, the IOM sampler cassettes were rinsed separately to determine if wall-loss corrections were necessary. All four samplers gave very good correlations between the two analytical methods above the limit of detection of the XRF procedure. The limit of detection for the 25-mm filters (5 microg) was lower than for the 37-mm filters (10 microg). The percentage of XRF results that were within 25% of the corresponding ICP results was evaluated. In addition, the bias from linear regression was estimated. Linear regression for the Button sampler and the IOM sampler using single readings and the GSP using all tested techniques for total filter loading gave acceptable XRF readings at loadings equivalent to sampling at the OSHA 8-hour Action Level and Permissible Exposure Limit. However, the CFC only had acceptable results when the center reading corrected for filter area was used, which was surprising, and may be a result of a limited data set. In addition to linear regression, simple estimation of bias indicated reasonable agreements between XRF and ICP results for single XRF readings on the Button sampler filters, (82% of the individual results within criterion), and on the IOM sampler filters (77% or 61%--see text), and on the GSP sampler filters using the OSHA algorithm (78%). As a result of this pilot project, all three samplers were considered suitable for inclusion in further field research studies.     

Evaluating performance of stormwater sampling approaches using a dynamic watershed model

Accurate quantification of stormwater pollutant levels is essential for estimating overall contaminant discharge to receiving waters. Numerous sampling approaches exist that attempt to balance accuracy against the costs associated with the sampling method. This study employs a novel and practical approach of evaluating the accuracy of different stormwater monitoring methodologies using stormflows and constituent concentrations produced by a fully validated continuous simulation watershed model. A major advantage of using a watershed model to simulate pollutant concentrations is that a large number of storms representing a broad range of conditions can be applied in testing the various sampling approaches. Seventy-eight distinct methodologies were evaluated by "virtual samplings" of 166 simulated storms of varying size, intensity and duration, representing 14 years of storms in Ballona Creek near Los Angeles, California. The 78 methods can be grouped into four general strategies: volume-paced compositing, time-paced compositing, pollutograph sampling, and microsampling. The performances of each sampling strategy was evaluated by comparing the (1) median relative error between the virtually sampled and the true modeled event mean concentration (EMC) of each storm (accuracy), (2) median absolute deviation about the median or "MAD" of the relative error or (precision), and (3) the percentage of storms where sampling methods were within 10% of the true EMC (combined measures of accuracy and precision). Finally, costs associated with site setup, sampling, and laboratory analysis were estimated for each method. Pollutograph sampling consistently outperformed the other three methods both in terms of accuracy and precision, but was the most costly method evaluated. Time-paced sampling consistently underestimated while volume-paced sampling over estimated the storm EMCs. Microsampling performance approached that of pollutograph sampling at a substantial cost savings. The most efficient method for routine stormwater monitoring in terms of a balance between performance and cost was volume-paced microsampling, with variable sample pacing to ensure that the entirety of the storm was captured. Pollutograph sampling is recommended if the data are to be used for detailed analysis of runoff dynamics.     

Assessing temporal representativeness of water quality monitoring data

The effectiveness of different monitoring methods in detecting temporal changes in water quality depends on the achievable sampling intervals, and how these relate to the extent of temporal variation. However, water quality sampling frequencies are rarely adjusted to the actual variation of the monitoring area. Manual sampling, for example, is often limited by the level of funding and not by the optimal timing to take samples. Restrictions in monitoring methods therefore often determine their ability to estimate the true mean and variance values for a certain time period or season. Consequently, we estimated how different sampling intervals determine the mean and standard deviation in a specific monitoring area by using high frequency data from in situ automated monitoring stations. Raw fluorescence measurements of chlorophyll a for three automated monitoring stations were calibrated by using phycocyanin fluorescence measurements and chlorophyll a analyzed from manual water samples in a laboratory. A moving block bootstrap simulation was then used to estimate the standard errors of the mean and standard deviations for different sample sizes. Our results showed that in a temperate, meso-eutrophic lake, relatively high errors in seasonal statistics can be expected from monthly sampling. Moreover, weekly sampling yielded relatively small accuracy benefits compared to a fortnightly sampling. The presented method for temporal representation analysis can be used as a tool in sampling design by adjusting the sampling interval to suit the actual temporal variation in the monitoring area, in addition to being used for estimating the usefulness of previously collected data.     

Estimation of Land Use Specific Runoff and Pollutant Concentration for Tapi River Basin in India

Non-point source (NPS) pollution is the result of various land use practices such as agriculture, sites of construction and waste disposal, urban development and so on. The control of NPS pollution is possible by regular monitoring and assessment on watershed basis to educate people for implementing well-known structural and non-structural measures. Recent trend is to use GIS based modelling tool for assessment of rainfall-runoff and non-point loading. The approach requires generation and analysis of basin wide data on various features of land and estimates of Event Mean Concentrations (EMCs) of pollutants in the runoff. In the present paper, basin wide data in different districts of Tapi basin has been analysed for land use distribution; fertilizer application; low, medium and high-density habitation; and annual rainfall. Coefficients of runoff have been estimated considering pervious and impervious area for different land use types, and compared with the reported values for Indian conditions. The estimated mean annual runoff flow indicated that two districts Jalgaon and Dhule contribute maximum runoff to the Tapi River. Estimates of EMCs for BOD and nutrients (N and P) in the runoff from various districts are useful in GIS-based modelling study for NPS pollution assessment.     

The magnitude of variability produced by methods used to estimate annual stormwater contaminant loads for highly urbanised catchments

Stormwater contaminant loading estimates using event mean concentration (EMC), rainfall/runoff relationship calculations and computer modelling (Model of Urban Stormwater Infrastructure Conceptualisation—MUSIC) demonstrated high variability in common methods of water quality assessment. Predictions of metal, nutrient and total suspended solid loadings for three highly urbanised catchments in Sydney estuary, Australia, varied greatly within and amongst methods tested. EMC and rainfall/runoff relationship calculations produced similar estimates (within 1 SD) in a statistically significant number of trials; however, considerable variability within estimates (～50 and ～25 % relative standard deviation, respectively) questions the reliability of these methods. Likewise, upper and lower default inputs in a commonly used loading model (MUSIC) produced an extensive range of loading estimates (3.8–8.3 times above and 2.6–4.1 times below typical default inputs, respectively). Default and calibrated MUSIC simulations produced loading estimates that agreed with EMC and rainfall/runoff calculations in some trials (4–10 from 18); however, they were not frequent enough to statistically infer that these methods produced the same results. Great variance within and amongst mean annual loads estimated by common methods of water quality assessment has important ramifications for water quality managers requiring accurate estimates of the quantities and nature of contaminants requiring treatment.     

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.     

Linking phosphorus export and hydrologic modeling: a case study in Central Italy

Soil erosion is an open topic, not only because soil fertility is lost, but also because nutrients are spilled into water bodies, thereby causing pollution. Research carried out in this field has amply described this process, but the interaction between these factors is complex and experimental research is needed to understand the production of loads of nutrients for different land uses. This paper describes a long-term monitoring case study using high-resolution rainfall data and runoff samples, carried out in the Lake Vico basin (Central Italy) to determine the phosphorus (P) export during erosive rainfall events. State of the art GIS-based basin characterization and advanced rainfall-runoff models are employed in order to describe the relationship between nutrient export and rainfall or runoff time distribution. Results show that the phosphorus export is strongly related to such time distributions, and less to the cumulative amount of rainfall or runoff.     

Status of the Niagara River point source discharge information: Sampling design and estimation of loading

Two basic requirements, low bias and high precision, are necessary for generating reliable estimates for the load from point and nonpoint sources of pollution. Biases and low precision can be the result of using a bad sampling design and/or inadequate method of estimation. The effects of biases can be reduced at the design stage prior to the data collection or at the data anlysis stage. This paper discusses the statistical issues involved in generating adequate load estimations using recently published point source discharge data from the Niagara River to illustrate these issues.