Agriculture, urbanization, and human activities, if not managed carefully, can expose a water body to environmental degradation, decreased water quality, and ultimately impaired benthic macroinvertebrate assemblage conditions. In streams where the benthic macroinvertebrates are impaired, the stream itself will not be meeting the water quality standards set forth in the Clean Water Act. As a result, the goal of this study was to establish relationships between benthic macroinvertebrates and their stressors so that stressor levels that would not adversely impact the benthic macroinvertebrates could be determined. Stressors such as sediment, habitat, water quality, landuse, watershed characteristics, and livestock numbers impact the benthic macroinvertebrate assemblage conditions. Since sediment is recognized as the Nation's leading pollutant and since the benthic macroinvertebrates live in the sediment on a stream bottom, this study placed emphasis upon the investigation of sediment as a primary stressor to the benthic macroinvertebrates. The specific objectives of this study were to develop relationships between the benthic macroinvertebrates and sediment and other stressors for Virginia streams, to evaluate the accuracy of the stressor/benthos relationships, and to discuss the implications of the study results for development of benthic TMDLs. A procedure to determine the relationships between stressors and benthic macroinvertebrate assemblage conditions was developed. Existing data on sediment, habitat, water quality, landuse, watershed characteristics, livestock numbers, and benthic macroinvertebrate assemblage conditions were compiled for 34 stations with 105 samples collected from the fall of 1996 to the fall of 1998. The 34 stations were located within 13 counties in Virginia (Rockbridge, Rockingham, Augusta, Frederick, Shenandoah, Page, Loudoun, Fairfax, Prince William, Fauquier, Culpeper, Rappahannock, and Madison) and in watersheds dominated by agricultural, urban, and forested landuses. Virginia currently uses the Rapid Bioassessment Protocol (RBP) method in its Biological Assessment Program. The RBP compares habitat and biological measures of the benthic macroinvertebrates to reference conditions using individual metrics. VADEQ's Biomonitoring Database, together with Ambient Water Quality Monitoring reports, GIS data layers, and VADCR's Hydrologic Unit Animal Census Database provided all of the necessary information for the stressor variables and benthic macroinvertebrate conditions. Accordingly, the stressor/benthos relationships were evaluated using statistical analyses procedures such as forward, backward, and stepwise multiple regression techniques; correlation analysis; principal component analysis; and r-square analysis. The statistical results indicated that sediment alone cannot be used to assess the benthic macroinvertebrate assemblage conditions. Other stressors such as dissolved oxygen, flow, % urban land, total suspended solids, temperature, stream velocity, substrate, hardness and alkalinity greatly impact the benthic macroinvertebrate assemblage conditions. The study results also indicate that the individual metrics within the RBP procedure are just as critical as the final RBP values in describing the benthic macroinvertebrate assemblage conditions. Upon completing the steps needed to develop stressor/benthos relationships, the validity of the relationships were verified for their application to other streams in Virginia. Validation was completed using 10 stations with 29 samples from the fall of 1996 to the fall of 1998. The 10 stations were located within 8 counties in Virginia (Bedford, Montgomery, Pulaski, Giles, Botetourt, Albemarle, Orange, and Culpeper) and in watersheds dominated by agricultural, urban, and forested landuses to correspond with the stations used to develop the stressor/benthos relationships. The implications of the relationships with regard to TMDLs were also studied using total suspended solids (TSS) loadings, turbidity levels, and embeddedness levels as the stressors of concern. The results for all benthic stations within Virginia showed that moderately impaired streams generally need to reduce the amount of embeddedness by 11, reduce turbidity levels by 5 FTU (57%), and reduce TSS values by 7 mg/L (68%) to meet a threshold value that would no longer adversely impact the benthic macroinvertebrates. Similarly, for the severely impaired stations throughout Virginia to meet threshold values, embeddedness amounts need to be reduced by 22, turbidity reduced by 57 FTU (93%), and TSS reduced by 74 mg/L (96%). This study was important since the proposed stressor/benthos relationships can provide policymakers with a useful tool to determine stressor thresholds that will not adversely impact the benthic macroinvertebrate assemblage conditions for use in developing benthic TMDLs in Virginia. The stressor/benthos relationships could also be used to determine the impact of certain activities or stressors on the benthic macroinvertebrates assemblage conditions in a given stream. / Master of Science

The contamination of inland surface water by point and nonpoint sources is a widespread human health concern. To address this problem, a thorough understanding of the sources and persistence of this pollution is necessary. Using fecal indicator bacteria (FIB) as a surrogate for enteric pathogens, the Environmental Protection Agency (EPA) and state agencies have developed surface water quality standards. If a segment does not meet these standards, a Total Maximum Daily Load (TMDL) must be developed for the watershed to identify sources of bacterial pollution. Currently, FIB pollution is the leading pollution type addressed by TMDLs in the country. One source of FIB identified in TMDLs is colonies of birds roosting under bridges. It has been proposed that the birds’ feces can augment the FIB concentrations downstream of bridges. In this year-long study of Bull Creek in Austin, Texas, it was determined that the concentrations of Escherichia coli and fecal coliform downstream of a bridge were significantly greater when migratory cliff swallows were nesting under the bridge. The downstream concentrations of both FIB exceeded contact recreation standards. Data also suggest that FIB from the feces could be swept into the stream by runoff from storm events. No enterococci loading was observed in conjunction with the swallows. This study also investigated the affects of sediment resuspension on surface water FIB concentrations and the persistence of FIB in shady, inland streams. The resuspension of sediments with attached FIB could also increase the FIB concentration in the water column and increase its persistence. The results of reactor-based experiments demonstrated that the concentration of E. coli in water from an effluent-dominated stream increased by a factor of 3 when riverine sediments were resuspended and exceeded single sample standards for contact recreation, suggesting sediments as a reservoir of E. coli. Additionally, concentrations of E. coli decreased by approximately 90% and 70% over 2 days in reactors containing stream water and sediment-laden stream water, respectively. / text

This study examined nonpoint source pollution via tributaries to Pineview Reservoir. Since few literature values of export coefficients are available for snowfall-dominated watersheds such as Pineview, locally scaled rates were quantified using an upstream-downstream bracketing technique. Nitrogen and phosphorus grab samples were manually collected and discharge measurements were conducted during the annual study period. Additionally, high-frequency monitoring sensors that measured EC, temperature, turbidity, and water level were deployed at the up- and downstream sites to represent short-duration transport events and to examine watershed processes on a more representative time scale. Daily nutrient loads were estimated from grab samples and flow rates using the Rank- Data (RD) distribution method and, using surrogate relationships for discharge and total phosphorus (TP) concentration from high-frequency sensors, half-hour TP loads were calculated. Short-duration snow melt events were identified by turbidity spikes and increased air temperatures as well as, in some cases, hydrograph peaks. During these events, export coefficients from 0.31 to 0.54 g TP/ha/hr were observed for low and high elevation snow melt events. While losses to ground water were observed for one study reach, the study reach with positive load gains had annual export coefficients of 0.018 g TP/ha/hr from high-frequency loads and 7.5e-5 g NOx-N/ha/hr, 3.5e-6 g SRP/ha/hr, and 9.1e-6 g TP/ha/hr from RD loads. These rates were 1,000 to 2,000 times greater than available literature values typical of rainfall-dominated watersheds. This study showed the importance of erosive, snow melt events on nutrient transport and the need for high-frequency monitoring representing short-duration events for accurate estimation of export coefficients. Management strategies to reduce nutrients from tributaries should focus on erosion control in the Pineview Reservoir watershed.

Under the Clean Water Act (CWA) program the Texas Commission on Environmental Quality (TCEQ) listed 110 stream segments with pathogenic bacteria impairment in 2000. The current study was conducted to characterize the watersheds associated with the impaired waterbodies. The main characteristics considered for the classification of waterbodies were designated use of the waterbody, land use distribution, density of stream network, average distance of a land of a particular use to the closest stream, household population, density of on-site sewage facilities (OSSF), bacterial loading due to the presence of different types of farm animals and wildlife, and average climatic conditions. The availability of observed in-stream fecal coliform bacteria concentration data was evaluated to obtain subgroups of data-rich and data-poor watersheds within a group. The climatic data and observed in-stream fecal coliform bacteria concentrations were analyzed to find out seasonal variability of the water quality. The watershed characteristics were analyzed using the multivariate statistical analysis techniques such as factor analysis/principal component analysis, cluster analysis, and discriminant analysis. Six groups of watersheds were formed as result of the statistical analysis. The main factors that differentiate the clusters were found to be bacterial contribution from farm animals and wildlife, density of OSSF, density of households connected to public sewers, and the land use distribution. Two watersheds were selected each from two groups of watersheds. Hydrological Simulation Program-FORTRAN (HSPF) model was calibrated for one watershed within each group and tested for the other watershed in the same group to study the similarity in the parameter sets due to the similarity in watershed characteristics. The study showed that the watersheds within a given cluster formed during the multivariate statistical analysis showed similar watershed characteristics and yielded similar model results for similar model input parameters. The effect of parameter uncertainty on the in-stream bacterial concentration predictions by HSPF was evaluated for the watershed of Salado Creek, in Bexar County. The parameters that control the HSPF model hydrology contributed the most variance in the in-stream fecal coliform bacterial concentrations corresponding to a simulation period between 1 January 1995 and 31 December 2000.

Under the Clean Water Act (CWA) program the Texas Commission on Environmental Quality (TCEQ) listed 110 stream segments with pathogenic bacteria impairment in 2000. The current study was conducted to characterize the watersheds associated with the impaired waterbodies. The main characteristics considered for the classification of waterbodies were designated use of the waterbody, land use distribution, density of stream network, average distance of a land of a particular use to the closest stream, household population, density of on-site sewage facilities (OSSF), bacterial loading due to the presence of different types of farm animals and wildlife, and average climatic conditions. The availability of observed in-stream fecal coliform bacteria concentration data was evaluated to obtain subgroups of data-rich and data-poor watersheds within a group. The climatic data and observed in-stream fecal coliform bacteria concentrations were analyzed to find out seasonal variability of the water quality. The watershed characteristics were analyzed using the multivariate statistical analysis techniques such as factor analysis/principal component analysis, cluster analysis, and discriminant analysis. Six groups of watersheds were formed as result of the statistical analysis. The main factors that differentiate the clusters were found to be bacterial contribution from farm animals and wildlife, density of OSSF, density of households connected to public sewers, and the land use distribution. Two watersheds were selected each from two groups of watersheds. Hydrological Simulation Program-FORTRAN (HSPF) model was calibrated for one watershed within each group and tested for the other watershed in the same group to study the similarity in the parameter sets due to the similarity in watershed characteristics. The study showed that the watersheds within a given cluster formed during the multivariate statistical analysis showed similar watershed characteristics and yielded similar model results for similar model input parameters. The effect of parameter uncertainty on the in-stream bacterial concentration predictions by HSPF was evaluated for the watershed of Salado Creek, in Bexar County. The parameters that control the HSPF model hydrology contributed the most variance in the in-stream fecal coliform bacterial concentrations corresponding to a simulation period between 1 January 1995 and 31 December 2000.

Land uses dominated by human activity can have a significant effect on ecological processes. In Florida, oxygen depletion is the most common impairment in lake, stream, and coastal water bodies. The continual growth and development in Florida, along with a conversion to more human intense land uses warrants study and discussion on impacts to dissolved oxygen regimes along a gradient of human disturbance. This research study is designed to identify observable trends in dissolved oxygen regimes along a gradient of increasing human intensity. Twenty-six stations in the Tampa Bay area were selected to represent lotic systems in west-central Florida. Data was collected quarterly, during four-day deployments, using a deployable data sonde. Grab samples for nutrients and chlorophyll-a provided antecedent data to explain observed trends. Physical components of streams, such as channelization were also taken into account. Biological integrity of streams was assessed to identify if altered dissolved oxygen regimes as a result of human land use significantly affect the health of the systems. Analysis included the use of Spearman rank order correlations to identify patterns. Dissolved oxygen regimes were correlated with the Landscape Development Intensity Index (LDI). Nutrients, primary productivity, and physical alteration to the streambed play a significant role in understanding how land use affects dissolved oxygen regimes. Results indicate the intensity of human land use has a significant effect on dissolved oxygen regimes and has significant policy implications for Florida's Total Maximum Daily Load (TMDL) program. Diel variation in oxygen measurements may be a more appropriate indicator of impairment and stream biological integrity.

Preventing impairment of waterbodies requires control, reduction and interception of contaminant losses at the field and subwatershed level. Three specific research objectives were accomplished in this study: 1) compare the HSPF-predicted flow, sediment, total nitrogen (TN) and total phosphorus (TP) loads resulting from simulation of spatially distributed site-specific and county-level disaggregated land use data at subwatershed and watershed levels, 2) evaluate the effects of site-specific and county-level disaggregated BMP data on modeled BMP responses in HSPF-predicted flow, sediment, TN and TP loads at the watershed level, and 3) analyze the long-term effects of the two spatial BMP datasets on achieving the Chesapeake Bay Total Maximum Daily Load (TMDL) goals for sediment, TN and TP. Site-specific data are derived from the local watershed inventory while disaggregated data are based on county-level aggregated data that are distributed to portions of river segments that intersect each county. The study site was the Upper Opequon Watershed in northern Virginia. Results for each research objective are as follows: (1) HSPF-predicted flow, sediment, TN and TP were higher using disaggregated land use data in subwatersheds at monthly and annual time-steps. (2) Predicted load reductions were higher with site-specific BMP data than with disaggregated data for the study watershed. (3) Current levels of cost-shared BMP implementation in the Upper Opequon Watershed using either site-specific or county-level disaggregated BMP datasets do not meet the Chesapeake Bay TMDL goals. Increasing BMP implementation level to 100% of the available land also failed to meet TMDL target goals. Generally, use of disaggregated land use data does not accurately represent the existing watershed conditions. Further, for the study watershed, use of disaggregated county-level BMP data poorly represented actual watershed conditions, which resulted in higher pollutant yields and higher levels of BMPs needed to meet water quality goals. The study suggests that site-specific land use and BMP data must be used during TMDL implementation planning to maintain credibility with local stakeholders and improve the accuracy of the developed implementation plans. / Ph. D.

High levels of fecal indicator bacteria (e.g. E. coli) are the leading cause of identified surface water impairments in the United States. The US Clean Water Act of 1972 requires that jurisdictions establish priority rankings for impaired waterways and develop a Total Maximum Daily Load (TMDL) plan for each. Although past research indicates that the risk of illness to humans varies by source of fecal contamination, current watershed assessments are developed according to total concentration of indicator bacteria, with all sources weighed equally. A stochastic model using Quantitative Microbial Risk assessment (QMRA) principles to translate source-specific (e.g. human, livestock) daily average concentrations of E.coli into a daily average risk of gastroenteritis infection was developed and applied to Pigg River, an impaired watershed in southern Virginia. Exposure was calculated by multiplying a ratio of source related reference pathogens to predicted concentrations of E.coli and a series of qualifying scalars. Risk of infection was then determined using appropriate dose response relationships. Overall, human and goose sources resulted in the greatest human health risk, despite larger overall E.coli loading associated with cattle. Bacterial load reductions specified in the Pigg River TMDL were applied using Hydrological Simulation Program- FORTRAN (HSPF) to assess the effect these reductions would have on the risk of infection attributed to each modeled bacterial source. Although individual risk sources (neglecting geese) were reduced below the EPA limit of 8 illnesses per 1000 exposures, the combined risk of illness varied between 0.006 and 64 illnesses per 1000 exposures. / Master of Science

During periods of reduced precipitation, flow in low-order, upland streams may be reduced and may stop completely. Under these "low flow" conditions, fecal bacteria directly deposited in the stream dominate in-stream bacteria loads. When developing a Total Maximum Daily Load (TMDL) to address a bacterial impairment in an upland, rural watershed, direct deposit (DD) fecal bacteria sources (livestock and wildlife defecating directly in the stream) often drive the source-load reductions required to meet water quality criteria. Due to limitations in the application of existing watershed-scale water quality models, under low-flow conditions the models can predict unrealistically high in-stream fecal bacteria concentrations. These unrealistically high simulated concentrations result in TMDL bacteria source reductions that are much more severe than what actually may be needed to meet applicable water quality criteria. This study used the Hydrological Simulation Program-FORTRAN (HSPF) to compare three low-flow DD simulation approaches and combinations (treatments) on two Virginia watersheds where bacterial impairment TMDLs had been previously developed and where low-flow conditions had been encountered. The three methods; Flow Stagnation (FS), DD Stage Cut-off (SC), and Stream Reach Surface Area (SA), have all been used previously to develop TMDLs. A modified version of the Climate Generation (CLIGEN) program was used to stochastically generate climate inputs for multiple model simulations. Violations of Virginia's interim fecal coliform criteria and the maximum simulated in-stream fecal coliform concentration were used to compare each treatment using ANOVA and Kruskal Wallis rank sum procedures. Livestock DD bacteria sources were incrementally reduced (100%, 50%, 15%, 10%, 5%) to represent TMDL load reduction allocation scenarios (allocation levels). Results from the first watershed indicate that the FS method simulated significantly lower instantaneous criterion violation rates at all allocation levels than the Control. The SC method reduced the livestock DD load compared to the Control, but produced significantly lower instantaneous criterion violation rates only at the 100% allocation level. The SA method did not produce significantly different instantaneous criterion violation rates compared to the Control. Geometric mean criterion violation rates were not significantly different from the Control at any allocation level. The distributions of maximum in-stream fecal coliform concentrations simulated by the combinations SC + FS and SC + SA + FS were both significantly different from the Control at the 100% allocation level. The second watershed did not produce low-flow conditions sufficient to engage the FS or SC methods. However, the SA method produced significantly different instantaneous violation rates than the Control at all allocation levels, which suggests that the SA method continues to affect livestock DD loads when low-flow conditions are not simulated in the watershed. No significant differences were found in the geometric mean violation rate or distribution of maximum simulated in-stream fecal coliform concentrations compared to the Control at any allocation level. This research suggests that a combination of the SC and FS methods may be the most appropriate treatment for addressing unrealistically high concentrations simulated during low-flow conditions. However, this combination must be used with caution as the FS method may increase the maximum simulated in-stream fecal coliform concentration if HSPF simulates zero volume within the reach. / Master of Science

For surface-water resource management studies, it is important to quantify all of the mechanisms that contribute to water quantity and influence water quality. In this regard, various methods have been used to ground-water fluxes in lake systems. These have included physical measurements (e.g., seepage meters), flow-net analyses, water budgets, chemical tracers, ground-water flow models, and statistical analyses. The method developed for this study for calculating ground-water inflow uses a simplified, 1-layer (surficial aquifer) ground-water flow model. The test area was on a set of lakes known as the Winter Haven Chain of Lakes in Polk County, Florida. The technique combines the use of a numerical model (MODFLOW) with an inverse prediction technique (PEST) to determine net surficial recharge rates. Within the model, the lakes were represented as constant-head boundaries. A general, surficial ground water no-flow boundary was delineated around the entire lake system based on the topographic boundaries. The model used annual average lake elevations to create a constant-head boundary for each lake for each year. Annual average elevations of surficial well heads were used as target well data. Model results generally support previous studies in the region, concluding that the lake chain receives significant inflow from the surficial aquifer and leaks to the Floridan aquifer. As a consequence, ground-water quality constituency was found to be of critical importance. One of the most important observations from this study is the need for accurate ground-water concentrations for ridge lake water quality management. The initial measured values used in this study were highly variable, uncertain, and likely underestimated the effect that ground water has on nutrient loading to the Winter Haven Chain of Lakes.