Critical Sampling Points Methodology: Case Studies of Geographically Diverse Watersheds

Only with a properly designed water quality monitoring network can data be collected that can lead to accurate information extraction. One of the main components of water quality monitoring network design is the allocation of sampling locations. For this purpose, a design methodology, called critical sampling points (CSP), has been developed for the determination of the critical sampling locations in small, rural watersheds with regard to total phosphorus (TP) load pollution. It considers hydrologic, topographic, soil, vegetative, and land use factors. The objective of the monitoring network design in this methodology is to identify the stream locations which receive the greatest TP loads from the upstream portions of a watershed. The CSP methodology has been translated into a model, called water quality monitoring station analysis (WQMSA), which integrates a geographic information system (GIS) for the handling of the spatial aspect of the data, a hydrologic/water quality simulation model for TP load estimation, and fuzzy logic for improved input data representation. In addition, the methodology was purposely designed to be useful in diverse rural watersheds, independent of geographic location. Three watershed case studies in Pennsylvania, Amazonian Ecuador, and central Chile were examined. Each case study offered a different degree of data availability. It was demonstrated that the developed methodology could be successfully used in all three case studies. The case studies suggest that the CSP methodology, in form of the WQMSA model, has potential in applications world-wide.     

A Water Quality Monitoring Network Design Methodology for the Selection of Critical Sampling Points: Part I

The principal instrument to temporally and spatially manage water resources is a water quality monitoring network. However, to date in most cases, there is a clear absence of a concise strategy or methodology for designing monitoring networks, especially when deciding upon the placement of sampling stations. Since water quality monitoring networks can be quite costly, it is very important to properly design the monitoring network so that maximum information extraction can be accomplished, which in turn is vital when informing decision-makers. This paper presents the development of a methodology for identifying the critical sampling locations within a watershed. Hence, it embodies the spatial component in the design of a water quality monitoring network by designating the critical stream locations that should ideally be sampled. For illustration purposes, the methodology focuses on a single contaminant, namely total phosphorus, and is applicable to small, upland, predominantly agricultural-forested watersheds. It takes a number of hydrologic, topographic, soils, vegetative, and land use factors into account. In addition, it includes an economic as well as logistical component in order to approximate the number of sampling points required for a given budget and to only consider the logistically accessible stream reaches in the analysis, respectively. The methodology utilizes a geographic information system (GIS), hydrologic simulation model, and fuzzy logic.     

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.     

基于GWLF模型的于桥水库入库河流水文模拟及氮磷负荷估算

流域模型技术应用是当前开展面源污染防治的重要工具,而水文过程的准确模拟是进行污染负荷估算的首要环节和关键步骤。为了弄清近年来于桥水库入库河流氮、磷输入负荷,选取GWLF模型对水平口子流域的水文过程进行模拟,首先利用2006—2018年气象、水文资料率定模型水文参数,然后将参数推广到整个流域,对2019—2020年3条主要入库河流流量进行模拟,最后乘以相应河流断面的总氮、总磷浓度估算氮磷输入负荷。结果显示：GWLF模型适用于研究区的水文过程模拟,校准期和验证期的纳氏系数分别为0.89和0.91,平均相对误差分别为12.2%和13.1%;2020年总氮入库负荷为3 977.0 t,其中引滦调水贡献占57.0%,3条入库河流共贡献43.0%;总磷入库负荷为48.8 t,其中引滦调水贡献占68.6%,3条入库河流共贡献31.4%。GWLF模型输入数据需求量较少,模型参数较少,模拟效果较好,适用于中小型流域的水资源和水环境管理,具有一定的推广应用前景。     

Nonpoint source pollution loading from an undistributed tropic forest area

Water quality and unit nonpoint sources (NPS) pollution load from a forest area were studied in a mountainous watershed in Taiwan. The flow rates were measured with rectangular weirs and samples taken for water quality analysis in both non-rainy and rainy days for 2 years. The subroutine of the Hydrological Simulation Program--FORTRAN was used to simulate runoff for additional 3 years. Total annual loads of various water quality parameters were then estimated by a regression model. Most of the parameter concentrations are higher during the rainy days; their values are typically higher as compared to data from other undisturbed forest areas. Nevertheless, the concentration ratio of dissolved inorganic nitrogen to TN or PO4(3-) -P to TP shows TN or TP no correlations with the flow rates, whereas the concentrations of SS and TP are positively correlated with the flow rate. The fluctuation of annual load from this watershed is significant. For example, six major events of the entire year, for which the total duration is merely 6.4 days, contribute 42% of the annual precipitation and at least 40% of the annual NPS loads. The management for controlling the NPS pollution from this forest watershed is discussed.     

Uncertainty in watershed response predictions induced by spatial variability of precipitation

Negligence to consider the spatial variability of rainfall could result in serious errors in model outputs. The objective of this study was to examine the uncertainty of both runoff and pollutant transport predictions due to the input errors of rainfall. This study used synthetic data to represent the “true” rainfall pattern, instead of interpolated precipitation. It was conducted on a synthetic case area having a total area of 20 km² with ten subbasins. Each subbasin has one rainfall gauge with synthetic precipitation records. Six rainfall storms with varied spatial distribution were generated. The average rainfall was obtained from all of the ten gauges by the arithmetic average method. The input errors of rainfall were induced by the difference between the actual rainfall pattern and estimated average rainfall. The results show that spatial variability of rainfall can cause uncertainty in modeling outputs of hydrologic, which would be transport to pollutant export predictions, when uniformity of rainfall is assumed. Since rainfall is essential information for predicting watershed responses, it is important to consider the properties of rainfall, particularly spatial rainfall variability, in the application of hydrologic and water quality models.     

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.     

Modeling Runoff Response to Land Cover and Rainfall Spatial Variability in Semi-Arid Watersheds

Hydrologic response is an integrated indicator of watershed condition, and significant changes in land cover may affect the overall health and function of a watershed. This paper describes a procedure for evaluating the effects of land cover change and rainfall spatial variability on watershed response. Two hydrologic models were applied on a small semi-arid watershed; one model is event-based with a one-minute time step (KINEROS), and the second is a continuous model with a daily time step (SWAT). The inputs to the models were derived from Geographic Information System (GIS) theme layers of USGS digital elevation models, the State Soil Geographic Database (STATSGO) and the Landsat-based North American Landscape Characterization classification (NALC) in conjunction with available literature and look up tables. Rainfall data from a network of 10 raingauges and historical stream flow data were used to calibrate runoff depth using the continuous hydrologic model from 1966 to 1974. No calibration was carried out for the event-based model, in which six storms from the same period were used in the calculation of runoff depth and peak runoff. The assumption on which much of this study is based is that land cover change and rainfall spatial variability affect the rainfall-runoff relationships on the watershed. To validate this assumption, simulations were carried out wherein the entire watershed was transformed from the 1972 NALC land cover, which consisted of a mixture of desertscrub and grassland, to a single uniform land cover type such as riparian, forest, oak woodland, mesquite woodland, desertscrub, grassland, urban, agriculture, and barren. This study demonstrates the feasibility of using widely available data sets for parameterizing hydrologic simulation models. The simulation results show that both models were able to characterize the runoff response of the watershed due to changes of land cover.     

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.     

Watershed-Based Survey Designs

Watershed-based sampling design and assessment tools help serve the multiple goals for water quality monitoring required under the Clean Water Act, including assessment of regional conditions to meet Section 305(b), identification of impaired water bodies or watersheds to meet Section 303(d), and development of empirical relationships between causes or sources of impairment and biological responses. Creation of GIS databases for hydrography, hydrologically corrected digital elevation models, and hydrologic derivatives such as watershed boundaries and upstream–downstream topology of subcatchments would provide a consistent seamless nationwide framework for these designs. The elements of a watershed-based sample framework can be represented either as a continuous infinite set defined by points along a linear stream network, or as a discrete set of watershed polygons. Watershed-based designs can be developed with existing probabilistic survey methods, including the use of unequal probability weighting, stratification, and two-stage frames for sampling. Case studies for monitoring of Atlantic Coastal Plain streams, West Virginia wadeable streams, and coastal Oregon streams illustrate three different approaches for selecting sites for watershed-based survey designs.     

Dynamic Hydrologic Simulation of the Bear Brook Watershed in Maine (BBWM)

Bear Brook Watershed in Maine (BBWM) consists of a pair of research watersheds, East Bear Brook (EBB) and West Bear Brook (WBB). Years of research and observations have shown both watersheds have high similarity in geographic and hydrologic characteristics; a simple comparison of hydrographs from these two watersheds further substantiates this similarity. The Object Watershed Link Simulation (OWLS) model was developed and used to simulate the hydrological processes within the BBWM. The OWLS model is a 3-dimensional, vector-based, visualized, physically-based, distributed watershed hydrologic model. Simulation results not only provide a close examination of hydrologic processes within a watershed, but also dynamically visualize the processes of flow separations and Variable Source Areas (VSA). Results from flow separations suggest that surface flow from riparian area is the predominate component for the flood rising limb and that macropore flow from riparian area dominates during the falling limb. Soil matrix flow has little effect flood period but is a persistent contributor to base flow. Results from VSA visualization demonstrate 3-D dynamic changes in surface flow distribution and suggest that downstream riparian areas are the major contributing area for peak flow. As water chemistry is highly relevant to the flow paths within a watershed, simulations have provided valuable information about source of stream flow and the water migration dynamics to support the study of watershed chemistry in the BBWM. More specific linkages between the chemistry behavior and the dynamic hydrologic processes should become the next simulation effort in the watershed study. There are many questions that are critical to watershed chemistry studies like: which flow component (surface flow, macropore flow, soil matrix flow) predominates during peak flows? How do the flow components distribute during a flood event? How do flow contributions differ between these two watersheds? Which portion of the watershed contributes the most to the peak flows? These questions remain unknown from previous observations and only can be addressed with a physically-based distributed model.     

Nutrient concentrations in Maryland non-tidal streams

Using a spatially extensive database from the Maryland Biological Stream Survey (MBSS), we describe nutrient relationships of small-order, non-tidal streams to Maryland watershed basins, Maryland Tributary Strategy basins, and stream order. In addition, we estimate the number of stream km affected by nutrient loading, using derived nutrient criteria. Based on the MBSS spring water quality sampling, we determined several important factors relating to nutrient levels in non-tidal streams. There are strong linear relationships of nutrients to the percentage of agriculture and forested land present within MBSS sampling strata. Both mean total nitrogen (TN) and mean total phosphorus (TP) levels for watershed basins by stream order show exceedances of derived nutrient reference criteria for Maryland. Four Maryland basins have over 85% of their stream kilometers exceeding the TN criterion, with three basins over 90% of the TP criterion. To protect small stream integrity in Maryland, we recommend an upper stream TN criterion between 1.34 and 1.68 mg/L and an upper stream TP criterion between 0.025 and 0.037 mg/L, based on quantile analyses. Elevated levels of both TN and TP are present in non-tidal streams, with subsequent nutrient inputs into the upper freshwater tidal reaches of the Chesapeake Bay.     

基于DPeRS模型的渭河典型断面汇水区面源污染评估及污染成因分析

近年来甘肃渭河桦林断面月度水质不稳定达标的问题引起了管理部门的广泛关注,掌握桦林断面汇水范围面源污染现状,对控制流域面源污染和促进水质稳定达标具有重要意义。采用遥感分布式污染估算（DPeRS）面源污染评估模型,对2018年黄河流域甘肃桦林断面汇水区面源污染空间分布特征进行分析,开展多类型污染量产排特征解析。结果表明：农业面源污染量方面,2018年甘肃桦林断面汇水区总氮（TN）、总磷（TP）、氨氮（NH⁺₄-N）、重铬酸盐指数（COD_Cr）面源污染排放量分别为11 591,2 697,7 141和1 458 t,入河量分别为2 184,512,1347,263 t;空间分布上,氮型（TN和NH⁺₄-N）排放负荷高值区主要分布在陇西县、武山县县段和岷县县段;武山县县段TP排放负荷较为突出;COD_Cr型面源污染高负荷区主要分布在陇西县、渭源县县段和武山县县段。农业面源污染物入河排放负荷空间分布差异明显,氮磷型（TN、NH⁺₄-N和TP）入河高负荷区主要分布在武山县县段、陇西县、临洮县县段;COD_Cr型面源污染入河高负荷区呈分散分布。漳县西部地区水土流失量较高,漳县西部、陇西县和渭源县县段北部局部地区泥沙负荷量较高。枯水期污染治理仍是保障水质稳定达标的关键期,农田径流是渭河桦林断面所在汇水区氮磷型面源污染的首要污染类型,畜禽养殖是COD_Cr型面源污染的首要污染类型。     

AgInput: An Agricultural Nutrient and Pesticide Source Model

Agriculture can be a major nonpoint source (NPS) of nutrient and pesticide contamination in the environment. Available databases do not provide accurate and dynamic data on fertilizer and pesticide application, which limits the ability of complex watershed models to simulate contaminant loads into impaired water bodies. A model for estimating agricultural nutrient and pesticide input for watershed modeling has been developed. Climate, soils, and major agricultural operations are considered within the model, so that it can be adapted to any watershed or subregion within a watershed. The timing of the agricultural operations is a function of the weather data, providing realistic results at daily, monthly, or annual application rates. The model also predicts irrigation demand and biomass production, which can be used to calibrate the model. Model output can be used in any watershed model that considers agricultural land uses. Two case studies were evaluated, using grape vineyards in the Napa River and strawberry production in Newport Bay as examples. The predicted time to maturity corresponded well with actual data. Irrigation and fertilizer needs were very sensitive to weather input. Although the model can generate weather from long-term averages, the simulated results are best when at least observed precipitation and temperature are provided, to capture extreme events. The model has data for 98 crops and 126 pesticides, based on the California Department of Pesticide Regulation database. The databases are easily modifiable by the user to adapt them to local conditions. The output from AgInput is much needed for watershed modeling and for development of total maximum daily loads (TMDLs), based on realistic targets of irrigation, nutrient, and pesticide inputs. The model is available for free download at .     

Water Quality Assessment of an Untreated Effluent Impacted Urban Stream: The Bharalu Tributary of the Brahmaputra River, India

Guwahati, the lone city on the bank of the entire midstream of the Brahmaputra River, is facing acute civic problem due to severe depletion of water quality of its natural water bodies. This work is an attempt towards water quality assessment of a relatively small tributary of the Brahmaputra called the Bharalu River flowing through the city that has been transformed today into a city drainage channel. By analyzing the key physical, chemical and biological parameters for samples drawn from different locations, an assessment of the dissolved load and pollution levels at different segments in the river was made. Locations where the contaminants exceeded the permissible limits during different seasons were identified by examining spatial and temporal variations. A GIS developed for the watershed with four layers of data was used for evaluating the influence of catchment land use characteristics. BOD, DO and total phosphorus were found to be the sensitive parameters that adversely affected the water quality of Bharalu. Relationship among different parameters revealed that the causes and sources of water quality degradation in the study area were due to catchments input, anthropogenic activities and poor waste management. Elevated levels of total phosphorus, BOD and depleted DO level in the downstream were used to develop an ANN model by taking total phosphorus and BOD as inputs and dissolved oxygen as output, which indicated that an ANN based predictive tool can be utilized for monitoring water quality in the future.     

湖州合溪新港污染物入湖通量估算与分析

根据2021年5月—2022年4月合溪新港河流水量、水质(TN和TP)的同步监测数据，利用通量模型核算了合溪新港污染物(TN和TP)通量。通过测算合溪新港TN、TP通量与断面降雨强度、水质的响应关系，分析了该区域的污染类型及特点，为后期水质污染调查及通量研究提供了新思路。结果表明：合溪新港流量与降雨量存在明显相关关系，在强降雨期(7—8月)水体流量最高，占监测周期总流量的57.77%;少雨期则流量最低，且会出现湖水倒灌现象(11—12月)。通过分析合溪新港TN、TP通量与流量、水质的相关关系，确定了该流域污染类型为点源污染及农业面源污染共存的混合型污染，且在高强度降雨时污染物负荷量较大。综上，可针对农业面源污染对该流域治理提出相关对策，建立农业面源污染防治体系，以有效降低TN和TP污染物的入湖通量，减少太湖TN和TP污染物负荷量。     

Assessment of surface water quality using multivariate statistical techniques in red soil hilly region: a case study of Xiangjiang watershed, China

In the study, multivariate statistical methods including factor, principal component and cluster analysis were applied to analyze surface water quality data sets obtained from Xiangjiang watershed, and generated during 7 years (1994-2000) monitoring of 12 parameters at 34 different profiles. Hierarchical cluster analysis grouped 34 sampling sites into three clusters, including relatively less polluted (LP), medium polluted (MP) and highly polluted (HP) sites, and based on the similarity of water quality characteristics, the watershed was divided into three zones. Factor analysis/principal component analysis, applied to analyze the data sets of the three different groups obtained from cluster analysis, resulted in four latent factors accounting for 71.62%, 71.77% and 72.01% of the total variance in water quality data sets of LP, MP and HP areas, respectively. The PCs obtained from factor analysis indicate that the parameters for water quality variations are mainly related to dissolve heavy metals. Thus, these methods are believed to be valuable to help water resources managers understand complex nature of water quality issues and determine the priorities to improve water quality.     

Use of neural networks for monitoring surface water quality changes in a neotropical urban stream

This paper reports the using of neural networks for water quality analysis in a tropical urban stream before (2002) and after sewerage building and the completion of point-source control-based sanitation program (2003). Mathematical modeling divided water quality data in two categories: (a) input of some in situ water quality variables (temperature, pH, O2 concentration, O2 saturation and electrical conductivity) and (b) water chemical composition (N-NO2(-); N-NO3(-); N-NH4(+) Total-N; P-PO4(3-); K+; Ca2+; Mg+2; Cu2+; Zn2+ and Fe+3) as the output from tested models. Stream water data come from fortnightly sampling in five points along the Ipanema stream (Southeast Brazil, Minas Gerais state) plus two points downstream and upstream Ipanema discharge into Doce River. Once the best models are consistent with variables behavior we suggest that neural networking shows potential as a methodology to enhance guidelines for urban streams restoration, conservation and management.     

Predicting Water Quality Impaired Stream Segments using Landscape-Scale Data and a Regional Geostatistical Model: A Case Study in Maryland

In the United States, probability-based water quality surveys are typically used to meet the requirements of Section 305(b) of the Clean Water Act. The survey design allows an inference to be generated concerning regional stream condition, but it cannot be used to identify water quality impaired stream segments. Therefore, a rapid and cost-efficient method is needed to locate potentially impaired stream segments throughout large areas. We fit a set of geostatistical models to 312 samples of dissolved organic carbon (DOC) collected in 1996 for the Maryland Biological Stream Survey using coarse-scale watershed characteristics. The models were developed using two distance measures, straight-line distance (SLD) and weighted asymmetric hydrologic distance (WAHD). We used the Corrected Spatial Akaike Information Criterion and the mean square prediction error to compare models. The SLD models predicted more variability in DOC than models based on WAHD for every autocovariance model except the spherical model. The SLD model based on the Mariah autocovariance model showed the best fit (r² = 0.72). DOC demonstrated a positive relationship with the watershed attributes percent water, percent wetlands, and mean minimum temperature, but was negatively correlated to percent felsic rock type. We used universal kriging to generate predictions and prediction variances for 3083 stream segments throughout Maryland. The model predicted that 90.2% of stream kilometers had DOC values less than 5 mg/l, 6.7% were between 5 and 8 mg/l, and 3.1% of streams produced values greater than 8 mg/l. The geostatistical model generated more accurate DOC predictions than previous models, but did not fit the data equally well throughout the state. Consequently, it may be necessary to develop more than one geostatistical model to predict stream DOC throughout Maryland. Our methodology is an improvement over previous methods because additional field sampling is not necessary, inferences about regional stream condition can be made, and it can be used to locate potentially impaired stream segments. Further, the model results can be displayed visually, which allows results to be presented to a wide variety of audiences easily.