Temporal and Spatial Trends in Nutrient and Sediment Loading to Lake Tahoe,California‐Nevada,USA

Since 1980, the Lake Tahoe Interagency Monitoring Program (LTIMP) has provided stream‐discharge and water quality data—nitrogen (N), phosphorus (P), and suspended sediment—at more than 20 stations in Lake Tahoe Basin streams. To characterize the temporal and spatial patterns in nutrient and sediment loading to the lake, and improve the usefulness of the program and the existing database, we have (1) identified and corrected for sources of bias in the water quality database; (2) generated synthetic datasets for sediments and nutrients, and resampled to compare the accuracy and precision of different load calculation models; (3) using the best models, recalculated total annual loads over the period of record; (4) regressed total loads against total annual and annual maximum daily discharge, and tested for time trends in the residuals; (5) compared loads for different forms of N and P; and (6) tested constituent loads against land use‐land cover (LULC) variables using multiple regression. The results show (1) N and P loads are dominated by organic N and particulate P; (2) there are significant long‐term downward trends in some constituent loads of some streams; and (3) anthropogenic impervious surface is the most important LULC variable influencing water quality in basin streams. Many of our recommendations for changes in water quality monitoring and load calculation methods have been adopted by the LTIMP.     

Toward Explaining Nitrogen and Phosphorus Trends in Chesapeake Bay Tributaries, 1992–2012

Understanding trends in stream chemistry is critical to watershed management, and often complicated by multiple contaminant sources and landscape conditions changing over varying time scales. We adapted spatially referenced regression (SPARROW) to infer causes of recent nutrient trends in Chesapeake Bay tributaries by relating observed fluxes during 1992, 2002, and 2012 to contemporary inputs and watershed conditions. The annual flow‐normalized nitrogen flux to the bay from its watershed declined by 14% to 127,000 Mg (metric tons) between 1992 and 2012, due primarily (more than 80% of the decline) to reduced point sources. The remainder of the decline was due to reduced atmospheric deposition (13%) and urban nonpoint sources. Agricultural inputs, which contribute most nitrogen to the bay, changed little, although trends in the average nitrogen yield (flux per unit area) from cropland and pasture to streams in some settings suggest possible effects of evolving nutrient applications or other land management practices. Point sources of phosphorus to local streams declined by half between 1992 and 2012, while nonpoint inputs were relatively unchanged. Annual phosphorus delivery to the bay increased by 9% to 9,570 Mg between 1992 and 2012, however, due mainly to reduced retention in the Susquehanna River at Conowingo Reservoir.     

Estimating Relevance of Organic Carbon,Nitrogen, and Phosphorus Loads to a Blackwater River Estuary1

Abstract: In blackwater river estuaries, a large portion of external carbon, nitrogen, and phosphorus load are combined in complex organic molecules of varying recalcitrance. Determining their lability is essential to establishing the relationship between anthropogenic loads and eutrophication. A method is proposed in which organic C, N, and P are partitioned into labile and refractory forms, based upon first‐order decay estimated by biochemical oxygen demand relative to total organic carbon, and C:N and C:P ratios as a function of organic carbon lability. The technique was applied in developing total maximum daily loads for the lower St. Johns, a blackwater Atlantic coastal plain river estuary in Northeast Florida. Point source organic nutrients were determined to be largely labile. Urban runoff was found to have the highest relative labile organic N and P content, followed by agricultural runoff. Natural forest and silviculture runoff were high in refractory organic N and P. Upstream labile C, N, and P loads were controlled by autochthonous production, with 34‐50% of summer total labile carbon imported as algal biomass. Differentiation of labile and refractory organic forms suggests that while anthropogenic nutrient enrichment has tripled the total nitrogen load, it has resulted in a 6.7‐fold increase in total labile nitrogen load.     

Development and Application of the 2010 Chesapeake Bay Watershed Total Maximum Daily Load Model

The Phase 5.3 Watershed Model simulates the Chesapeake watershed land use, river flows, and the associated transport and fate of nutrient and sediment loads to the Chesapeake Bay. The Phase 5.3 Model is the most recent of a series of increasingly refined versions of a model that have been operational for more than two decades. The Phase 5.3 Model, in conjunction with models of the Chesapeake airshed and estuary, provides estimates of management actions needed to protect water quality, achieve Chesapeake water quality standards, and restore living resources. The Phase 5.3 Watershed Model tracks nutrient and sediment load estimates of the entire 166,000 km² watershed, including loads from all six watershed states. The creation of software systems, input datasets, and calibration methods were important aspects of the model development process. A community model approach was taken with model development and application, and the model was developed by a broad coalition of model practitioners including environmental engineers, scientists, and environmental managers. Among the users of the Phase 5.3 Model are the Chesapeake watershed states and local governments, consultants, river basin commissions, and universities. Development and application of the model are described, as well as key scenarios ranging from high nutrient and sediment load conditions if no management actions were taken in the watershed, to low load estimates of an all‐forested condition.     

Development of the Chesapeake Bay Watershed Total Maximum Daily Load Allocation

Nutrient load allocations and subsequent reductions in total nitrogen and phosphorus have been applied in the Chesapeake watershed since 1992 to reduce hypoxia and to restore living resources. In 2010, sediment allocations were established to augment nutrient allocations supporting the submerged aquatic vegetation resource. From the initial introduction of nutrient allocations in 1992 to the present, the allocations have become more completely applied to all areas and loads in the watershed and have also become more rigorously assessed and tracked. The latest 2010 application of nutrient and sediment allocations were made as part of the Chesapeake Bay total maximum daily load and covered all six states of the Chesapeake watershed. A quantitative allocation process was developed that applied principles of equity and efficiency in the watershed, while achieving all tidal water quality standards through an assessment of equitable levels of effort in reducing nutrients and sediments. The level of effort was determined through application of two key watershed scenarios: one where no action was taken in nutrient control and one where maximum nutrient control efforts were applied. Once the level of effort was determined for different jurisdictions, the overall load reduction was set watershed‐wide to achieve dissolved oxygen water quality standards. Further adjustments were made to the allocation to achieve the James River chlorophyll‐a standard.     

Evaluation of regression methodology with low-frequency water quality sampling to estimate constituent loads for ephemeral watersheds in Texas

Water quality regulation and litigation have elevated the awareness and need for quantifying water quality and source contributions in watersheds across the USA. In the present study, the regression method, which is typically applied to large (perennial) rivers, was evaluated in its ability to estimate constituent loads (NO(3)-N, total N, PO(4)-P, total P, sediment) on three small (ephemeral) watersheds with different land uses in Texas. Specifically, regression methodology was applied with daily flow data collected with bubbler stage recorders in hydraulic structures and with water quality data collected with four low-frequency sampling strategies: random, rise and fall, peak, and single stage. Estimated loads were compared with measured loads determined in 2001-2004 with an autosampler and high-frequency sampling strategies. Although annual rainfall and runoff volumes were relatively consistent within watersheds during the study period, measured annual nutrient and sediment concentrations and loads varied considerably for the cultivated and mixed watersheds but not for the pasture watershed. Likewise, estimated loads were much better for the pasture watershed than the cultivated and mixed landuse watersheds because of more consistent land management and vegetation type in the pasture watershed, which produced stronger correlations between constituent loads and mean daily flow rates. Load estimates for PO(4)-P were better than for other constituents possibly because PO(4)-P concentrations were less variable within storm events. Correlations between constituent concentrations and mean daily flow rate were poor and not significant for all watersheds, which is different than typically observed in large rivers. The regression method was quite variable in its ability to accurately estimate annual nutrient loads from the study watersheds; however, constituent load estimates were much more accurate for the combined 3-yr period. Thus, it is suggested that for small watersheds, regression-based annual load estimates should be used with caution, whereas long-term estimates can be much more accurate when multiple years of concentration data are available. The predictive ability of the regression method was similar for all of the low-frequency sampling strategies studied; therefore, single-stage or random strategies are recommended for low-frequency storm sampling on small watersheds because of their simplicity.     

Groundwater Nutrient Concentrations and Mass Loading Rates at Iowa Golf Courses

Golf courses are often considered by the public to be significant nitrogen (N) and phosphorus (P) nonpoint sources but only limited information exists on nutrient concentrations and loads in golf course groundwater. In this study, we measured N and P concentrations in groundwater and available surface water at six randomly selected Iowa golf courses to assess the loading risk posed by these facilities to groundwater and local rivers. At each course, three shallow monitoring wells were installed, one each on representative tee, fairway, and rough locations. Wells and available surface water were sampled on eight occasions during 2015 and 2016. NO₃‐N concentrations were not detected above 1 mg/L at three of the six courses monitored in this study and the overall mean NO₃‐N concentration in Iowa golf courses was 2.2 mg/L. The mass of NO₃‐N recharged to groundwater averaged 3.3 kg/ha at the six courses, which represents approximately one‐tenth of the NO₃‐N load exported by the watershed that contains the course and represented approximately 0.1 to 8% of the fertilizer N applied. Groundwater orthophosphorus concentrations averaged 0.13 mg/L and were similar to those measured in a variety of settings across Iowa. Study results should prove useful in evaluating nutrient contributions from golf courses in Midwestern states where nutrient reduction strategies are being pursued.     

Integration of SWAT and HSPF for Simulation of Sediment Sources in Legacy Sediment‐Impacted Agricultural Watersheds

A total maximum daily load for the Chesapeake Bay requires reduction in pollutant load from sources within the Bay watersheds. The Conestoga River watershed has been identified as a major source of sediment load to the Bay. Upland loads of sediment from agriculture are a concern; however, a large proportion of the sediment load in the Conestoga River has been linked to scour of legacy sediment associated with historic millpond sites. Clarifying this distinction and identifying specific segments associated with upland vs. channel sources has important implications for future management. In order to address this important question, we combined the strengths of two widely accepted watershed management models — Soil and Water Assessment Tool (SWAT) for upland agricultural processes, and Hydrologic Simulation Program FORTRAN (HSPF) for instream fate and transport — to create a novel linked modeling system to predict sediment loading from critical sources in the watershed including upland and channel sources, and to aid in targeted implementation of management practices. The model indicates approximately 66% of the total sediment load is derived from instream sources, in agreement with other studies in the region and can be used to support identification of these channel source segments vs. upland source segments, further improving targeted management. The innovated linked SWAT‐HSPF model implemented in this study is useful for other watersheds where both upland agriculture and instream processes are important sources of sediment load.     

Monitored and Modeled Correlations of Sediment and Nutrients with Chesapeake Bay Water Clarity

This article analyzes the correlations of the observed and modeled light attenuation coefficient, Kd, with in situ total suspended solids (TSS) and chlorophyll‐a concentrations in Chesapeake Bay (CB) tidal waters, and with sediment and nutrient loads from the Chesapeake watershed. Light attenuation is closely related to in situ TSS and chlorophyll‐a concentrations, however, the strength of the correlation differs among the CB segments. There are distinct differences between saline and tidal fresh segments in the main Bay, but less distinction among saline and tidal fresh segments in the tidal tributaries. The correlation between Kd with sediment and nutrient loads is complicated by the lag times of TSS and the chlorophyll‐a responses to reductions in nutrient and sediment loads from the watershed, and also due to the diverse load sources. Three sets of model sensitivity scenarios were performed with: (1) differential sediment and nutrient loads; (2) selective sediment source types; and (3) geographically isolated inputs. The model results yield similar findings as those based on observed data and provide information regarding the effect of sediment on specific water bodies. Based on the model results a method was developed to determine sediment and nutrient load reductions needed to achieve the water clarity standards of the CB segments.     

MODELING THE HYDROCHEMISTRY OF THE CANNONSVILLE WATERSHED WITH GENERALIZED WATERSHED LOADING FUNCTIONS (GWLF)1

ABSTRACT: A previous modeling study used the Generalized Watershed Loading Functions (GWLF) model to simulate stream‐flow, and nutrient and sediment loads to Cannonsville Reservoir from the West Branch Delaware River (WBDR). We made several model revisions, calibrated key parameters, and tested the original GWLF model and a revised GWLF model using more recent data. Model revisions included: addition of unsaturated leakage between unsaturated and saturated subsurface reservoirs; revised timing of sediment export; inclusion of urban sediments and dissolved nutrients; tracking of particulate nutrients from point sources; and revised timing of septic system loads. The revision of sediment yield timing resulted in significant improvements in monthly sediment and particulate phosphorus predictions as compared to the original model. Addition of unsaturated leakage improved hydrologic predictions during low flow months. The other model changes improve realism without adding significant model complexity or data requirements. Goodness of fit of revised model predictions versus stream measurements, as measured by the Nash‐Sutcliff coefficient of model efficiency, exceeded 0.8 for streamflow‐0.7 for sediment yield and dissolved nitrogen (N) and 0.6 for particulate and dissolved phosphorus (P). The revised GWLF model, with limited calibration, provides reasonable estimates of monthly streamflow, and nutrient and sediment loads in the Cannonsville watershed.     

THE EFFECTS OF CLIMATE CHANGE ON STREAM FLOW AND NUTRIENT LOADING'

ABSTRACT: This study assesses the potential impact of climate change on stream flow and nutrient loading in six watersheds of the Susquehanna River Basin using the Generalized Watershed Loading Function (GWLF). The model was used to simulate changes in stream flow and nutrient loads under a transient climate change scenario for each watershed. Under an assumption of no change in land cover and land management, the model was used to predict monthly changes in stream flow and nutrient loads for future climate conditions. Mean annual stream flow and nutrient loads increased for most watersheds, but decreased in one watershed that was intensively cultivated. Nutrient loading slightly decreased in April and late summer for several watersheds as a result of early snowmelt and increasing evapotranspiration. Spatial and temporal variability of stream flow and nutrient loads under the transient climate scenario indicates that different approaches for future water resource management may be useful.     

Nitrogen loads through baseflow, stormflow, and underflow to Rehoboth Bay, Delaware

A detailed study of water and nitrogen (N) discharge from a small, representative subwatershed of Rehoboth Bay, Delaware, was conducted to determine total N loads to the bay. The concentrations of ammonium (NH4(+)), nitrate + nitrite (NO3(-) + NO2(-)), and dissolved and particulate organic N were determined in baseflow and storm waters discharging from Bundicks Branch from October 1998 to April 2002. A novel hydrographic separation model that accounts for significant decreases in baseflow during storm events was developed to estimate N loads during unsampled storms. Nitrogen loads based on gauged flows alone (7100-19,100 kg/yr) significantly underestimated those based on land use-land cover (LULC) and estimated N export factors from different classes of LULC (32,000-40,600 kg/yr). However, when ungauged underflow and associated N loads were included in the total loads (25,500-33,800 kg/yr), there was much better agreement with LULC export models. This suggests that in permeable coastal plain sediments, underflow contributes significantly to N fluxes to estuarine receiving waters, particularly in drier years. Based on the similarity in LULC, N loads from the Bundicks Branch subwatershed were used to estimate upland loads to the entire Rehoboth Bay Watershed (259,000-316,000 kg/yr). These N loads from the watershed were much greater than those from direct atmospheric deposition (49,000-64,500 kg/yr) and from a local wastewater treatment plant (9700-13,700 kg/yr). While the watershed was the principal source of N at all times during the year, the relative contributions from the watershed, wastewater, and direct atmospheric deposition varied predictably with season.     

The Use of Carbon and Nitrogen Isotopes to Study Watershed Erosion Processes1

Abstract:  Tracer studies are needed to better understand watershed soil erosion and calibrate watershed erosion models. For the first time, stable nitrogen and carbon isotopes (δ¹⁵N and δ¹³C) and the carbon to nitrogen atomic ratio (C/N) natural tracers are used to investigate temporal and spatial variability of erosion processes within a sub‐watershed. Temporal variability was assessed by comparing δ¹⁵N, δ¹³C, and C/N of eroded‐soils from a non‐equilibrium erosion event immediately following freezing and thawing of surface soils with two erosion events characterized by equilibrium conditions with erosion downcutting. Spatial variability was assessed for the equilibrium events by using the δ¹⁵N and δ¹³C signatures of eroded‐soils to measure the fraction of eroded‐soil derived from rill/interrill erosion on upland hillslopes as compared to headcut erosion on floodplains. In order to perform this study, a number of tasks were carried out including: (1) sampling source‐soils from upland hillslopes and floodplains, (2) sampling eroded‐soils with an in situ trap in the stream of the sub‐watershed, (3) isotopic and elemental analysis of the samples using isotope ratio mass spectrometry, (4) fractioning eroded‐soil to its upland rill/interrill and floodplain headcut end‐members using an unmixing model within a Bayesian Markov Chain Monte Carlo framework, and (5) evaluating tracer unmixing model results by comparison with process‐based erosion prediction models for rill/interrill and headcut erosion processes. Results showed that finer soil particles eroded during the non‐equilibrium event were enriched in δ¹⁵N and δ¹³C tracers and depleted in C/N tracer relative to coarser soil particles eroded during the equilibrium events. Correlation of tracer signature with soil particle size was explainable based on known biogeochemical processes. δ¹⁵N and δ¹³C were also able to distinguish between upland rill/interrill erosion and floodplain headcut erosion, which was due to different plant cover at the erosion sources. Results from the tracer unmixing model highlighted future needs for coupling rill/interrill and headcut erosion prediction models.     

Modeling Nutrient Release in the Tai Lake Basin of China: Source Identification and Policy Implications

Because nutrient enrichment has become increasingly severe in the Tai Lake Basin of China, identifying sources and loads is crucial for watershed nutrient management. This paper develops an empirical framework to estimate nutrient release from five major sectors, which requires fewer input parameters and produces acceptable accuracy. Sectors included are industrial manufacturing, livestock breeding (industrial and family scale), crop agriculture, household consumption (urban and rural), and atmospheric deposition. Results show that in the basin (only the five sectors above), total nutrient loads of nitrogen (N) and phosphorus (P) into aquatic systems in 2008 were 33043.2 tons N a^?1 and 5254.4 tons P a^?1, and annual area-specific nutrient loads were 1.94 tons N km^?2 and 0.31 tons P km^?2. Household consumption was the major sector having the greatest impact (46 % in N load, 47 % in P load), whereas atmospheric deposition (18 %) and crop agriculture (15 %) sectors represented other significant proportions of N load. The load estimates also indicate that 32 % of total P came from the livestock breeding sector, making it the second largest phosphorus contributor. According to the nutrient pollution sectors, six best management practices are selected for cost-effectiveness analysis, and feasible options are recommended. Overall, biogas digester construction on industrial-scale farms is proven the most cost-effective, whereas the building of rural decentralized facilities is the best alternative under extreme financial constraint. However, the reduction potential, average monetary cost, and other factors such as risk tolerance of policy makers should all be considered in the actual decision-making process.     

Herbicide Transport in Goodwater Creek Experimental Watershed: II. Long‐Term Research on Acetochlor,Alachlor, Metolachlor,and Metribuzin1

Lerch, R.N., E.J. Sadler, C. Baffaut, N.R. Kitchen, and K.A. Sudduth, 2010. Herbicide Transport in Goodwater Creek Experimental Watershed: II. Long‐Term Research on Acetochlor, Alachlor, Metolachlor, and Metribuzin. Journal of the American Water Resources Association (JAWRA) 1‐15. DOI: 10.1111/j.1752‐1688.2010.00504.x Abstract: Farmers in the Midwestern United States continue to be reliant on soil‐applied herbicides for weed control in crop production, and herbicide contamination of streams remains an environmental problem. The main objective of this study was to analyze trends in concentration and load of acetochlor, alachlor, metolachlor, and metribuzin in Goodwater Creek Experimental Watershed (GCEW) from 1992 to 2006. A secondary objective was to document the effects of best management practices (BMPs) implemented within GCEW on herbicide transport trends. Median relative herbicide loads, as a percent of applied, were 3.7% for metolachlor, 1.3% for metribuzin, 0.36% for acetochlor, and 0.18% for alachlor. The major decrease in alachlor use and increase in acetochlor use caused shifts in flow‐weighted concentrations that were observed over the entire concentration range. The smaller decrease in metolachlor use led to a consistent decreasing time trend only for the upper end of the concentration distribution. Metribuzin also showed moderate decreases in concentration with time since 1998. Annual loads were generally correlated to second quarter discharge. Despite extensive education efforts in the watershed, conservation BMPs within GCEW were mainly implemented to control erosion, and therefore had no discernable impact on reducing herbicide transport. Overall, changes in herbicide use and second quarter discharge had the greatest effect on trends in flow‐weighted concentration and annual load.     

Hot Spots and Hot Moments in Riparian Zones: Potential for Improved Water Quality Management1

Vidon, Philippe, Craig Allan, Douglas Burns, Tim P. Duval, Noel Gurwick, Shreeram Inamdar, Richard Lowrance, Judy Okay, Durelle Scott, and Steve Sebestyen, 2010. Hot Spots and Hot Moments in Riparian Zones: Potential for Improved Water Quality Management. Journal of the American Water Resources Association (JAWRA) 46(2):278-298. DOI: 10.1111/j.1752-1688.2010.00420.x Abstract: Biogeochemical and hydrological processes in riparian zones regulate contaminant movement to receiving waters and often mitigate the impact of upland sources of contaminants on water quality. These heterogeneous processes have recently been conceptualized as “hot spots and moments” of retention, degradation, or production. Nevertheless, studies investigating the importance of hot phenomena (spots and moments) in riparian zones have thus far largely focused on nitrogen (N) despite compelling evidence that a variety of elements, chemicals, and particulate contaminant cycles are subject to the influence of both biogeochemical and transport hot spots and moments. In addition to N, this review summarizes current knowledge for phosphorus, organic matter, pesticides, and mercury across riparian zones, identifies variables controlling the occurrence and magnitude of hot phenomena in riparian zones for these contaminants, and discusses the implications for riparian zone management of recognizing the importance of hot phenomena in annual solute budgets at the watershed scale. Examples are presented to show that biogeochemical process-driven hot spots and moments occur along the stream/riparian zone/upland interface for a wide variety of constituents. A basic understanding of the possible co-occurrence of hot spots and moments for a variety of contaminants in riparian systems will increase our understanding of the influence of riparian zones on water quality and guide management strategies to enhance nutrient or pollutant removal at the landscape scale.     

Assessment of chlorophyll-a as a criterion for establishing nutrient standards in the streams and rivers of Illinois

Nutrient enrichment is a frequently cited cause for biotic impairment of streams and rivers in the USA. Efforts are underway to develop nutrient standards in many states, but defensible nutrient standards require an empirical relationship between nitrogen (N) or phosphorus (P) concentrations and some criterion that relates nutrient levels to the attainment of designated uses. Algal biomass, measured as chlorophyll-a (chl-a), is a commonly proposed criterion, yet nutrient-chl-a relationships have not been well documented in Illinois at a state-wide scale. We used state-wide surveys of >100 stream and river sites to assess the applicability of chl-a as a criterion for establishing nutrient standards for Illinois. Among all sites, the median total P and total N concentrations were 0.185 and 5.6 mg L(-1), respectively, during high-discharge conditions. During low-discharge conditions, median total P concentration was 0.168 mg L(-1), with 25% of sites having a total P of > or =0.326 mg L(-1). Across the state, 90% of the sites had sestonic chl-a values of < or =35 microg L(-1), and watershed area was the best predictor of sestonic chl-a. During low discharge there was a significant correlation between sestonic chl-a and total P for those sites that had canopy cover < or =25% and total P of < or =0.2 mg L(-1). Results suggest sestonic chl-a may be an appropriate criterion for the larger rivers in Illinois but is inappropriate for small rivers and streams. Coarse substrate to support benthic chl-a occurred in <50% of the sites we examined; a study using artificial substrates did not reveal a relationship between chl-a accrual and N or P concentrations. For many streams and rivers in Illinois, nutrients may not be the limiting factor for algal biomass due to the generally high nutrient concentrations and the effects of other factors, such as substrate conditions and turbidity.     

Influence of Sampling Frequency on Estimation of Annual Total Phosphorus and Total Suspended Solids Loads1

Abstract: The determination of sediment and nutrient loads is typically based on the collection and analysis of grab samples. The frequency and regularity of traditional sampling may not provide representation of constituent loading, particularly in systems with flashy hydrology. At two sites in the Little Bear River, Utah, continuous, high‐frequency turbidity was used with surrogate relationships to generate estimates of total phosphorus and total suspended solids concentrations, which were paired with discharge to estimate annual loads. The high frequency records were randomly subsampled to represent hourly, daily, weekly, and monthly sampling frequencies and to examine the effects of timing, and resulting annual load estimates were compared to the reference loads. Higher frequency sampling resulted in load estimates that better approximated the reference loads. The degree of bias was greater at the more hydrologically responsive site in the upper watershed, which required a higher sampling frequency than the lower watershed site to achieve the same level of accuracy in estimating the reference load. The hour of day and day of week of sampling impacted load estimation, depending on site and hydrologic conditions. The effects of sampling frequency on the determination of compliance with a water quality criterion were also examined. These techniques can be helpful in determining necessary sampling frequency to meet the objectives of a water quality monitoring program.