Water quality and restoration in a coastal subdivision stormwater pond

Stormwater ponds are commonly used in residential and commercial areas to control flooding. The accumulation of urban contaminants in stormwater ponds can lead to a number of water quality problems including high nutrient, chemical contaminant, and bacterial levels. This study examined the interaction between land use and coastal pond water quality in a South Carolina residential subdivision pond. Eutrophic levels of chlorophyll and phosphorus were present in all seasons. Harmful cyanobacterial blooms were prevalent during the summer months. Microcystin toxin and fecal coliform bacteria levels were measured that exceeded health and safety standards. Low concentrations of herbicides (atrazine and 2,4-D) were also detected during summer months. Drainage from the stormwater pond may transport contaminants into the adjacent tidal creek and estuary. A survey of residents within the pond's watershed indicated poor pet waste management and frequent use of fertilizers and pesticides as possible contamination sources. Educational and outreach activities were provided to community members to create an awareness of the water quality conditions in the pond. Pond management strategies were then recommended, and selected mitigation actions were implemented. Water quality problems identified in this study have been observed in other coastal stormwater ponds of varying size and salinity, leading this project to serve as a potential model for coastal stormwater pond management.     

IMPACTS OF URBANIZATION ON NUTRIENT CONCENTRATIONS IN SMALL SOUTHEASTERN COASTAL STREAMS1

ABSTRACT: Coastal watersheds in the southeastern United States are rapidly changing due to population growth and attendant increases in residential development, industry, and tourism related commerce. This research examined spatial and temporal patterns of nutrient concentrations in streams from 10 small watersheds (< 4 km²) that drain into Murrells Inlet (impacted) and North Inlet (pristine), two high salinity estuaries along the South Carolina coast. Monthly grab samples were collected during baseflow during 1999 and analyzed for total and dissolved inorganic and organic forms of nitrogen and phosphorus. Data were grouped into forested wetland creeks (representing predevelopment reference sites), urban creeks, and urban ponds. DON and NH4 concentrations were greater in forested streams than in urban streams. NO3 and TP concentrations were greatest in urban streams. Seasonally, concentrations were highest during summer for TN, NH4, DON, and TP, while NO3 concentrations were greatest during winter. Nutrient ratios clearly highlighted the reduction in organic nitrogen due to coastal development. Multiple regression models to predict instream nutrient concentrations from land use in Murrells Inlet suggest that effects are not significant (small r²). The findings indicate that broad land use/land cover classes cannot be used to predict nutrient concentrations in streams in the very small watersheds in our study areas.     

GROUPWISE MODELING STUDY OF BACTERIALLY IMPAIRED WATERSHEDS IN TEXAS: CLUSTERING ANALYSIS1

ABSTRACT: Under the Clean Water Act (CWA) program, the Texas Commission on Environmental Quality (TCEQ) listed 110 stream segments in the year 2000 with pathogenic bacteria impairment. A study was conducted to evaluate the probable sources of pollution and characterize the watersheds associated with these impaired water bodies. The primary aim of the study was to group the water bodies into clusters having similar watershed characteristics and to examine the possibility of studying them as a group by choosing models for total maximum daily load (TMDL) development based on their characteristics. This approach will help to identify possible sources and determine appropriate models and hence reduce the number of required TMDL studies. This in turn will help in reducing the effort required to restore the health of the impaired water bodies in Texas. The main characteristics considered for the classification of water bodies were land use distribution within the watershed, density of stream network, average distance of land of a particular use to the closest stream, household population, density of on‐site sewage facilities (OSSFs), bacterial loading from different types of farm animals and wildlife, and average climatic conditions. The climatic data and observed instream fecal coliform bacteria concentrations were analyzed to evaluate seasonal variability of instream water quality. The grouping of water bodies was carried out using the multivariate statistical techniques of factor analysis/principal component analysis, cluster analysis, and discriminant analysis. The multivariate statistical analysis resulted in six clusters of water bodies. The main factors that differentiated the clusters were found to be bacterial contribution from farm animals and wildlife, density of OSSFs, density of households connected to public sewers, and land use distribution.     

Interaction and influence of two creeks on Escherichia coli concentrations of nearby beaches: exploration of predictability and mechanisms

The impact of river outfalls on beach water quality depends on numerous interacting factors. The delivery of contaminants by multiple creeks greatly complicates understanding of the source contributions, especially when pollution might originate up- or down-coast of beaches. We studied two beaches along Lake Michigan that are located between two creek outfalls to determine the hydrometeorologic factors influencing near-shore microbiologic water quality and the relative impact of the creeks. The creeks continuously delivered water with high concentrations of Escherichia coli to Lake Michigan, and the direction of transport of these bacteria was affected by current direction. Current direction reversals were associated with elevated E. coli concentrations at Central Avenue beach. Rainfall, barometric pressure, wave height, wave period, and creek specific conductance were significantly related to E. coli concentration at the beaches and were the parameters used in predictive models that best described E. coli variation at the two beaches. Multiple inputs to numerous beaches complicates the analysis and understanding of the relative relationship of sources but affords opportunities for showing how these complex creek inputs might interact to yield collective or individual effects on beach water quality.     

A Modeling System to Assess Land Cover Land Use Change Effects on SAV Habitat in the Mobile Bay Estuary

Estuarine ecosystems are largely influenced by watersheds directly connected to them. In the Mobile Bay, Alabama watersheds we examined the effect of land cover and land use (LCLU) changes on discharge rate, water properties, and submerged aquatic vegetation, including freshwater macrophytes and seagrasses, throughout the estuary. LCLU scenarios from 1948, 1992, 2001, and 2030 were used to influence watershed and hydrodynamic models and evaluate the impact of LCLU change on shallow aquatic ecosystems. Overall, our modeling results found that LCLU changes increased freshwater flows into Mobile Bay altering temperature, salinity, and total suspended sediments (TSS). Increased urban land uses coupled with decreased agricultural/pasture lands reduced TSS in the water column. However, increased urbanization or agricultural/pasture land coupled with decreased forest land resulted in higher TSS concentrations. Higher sediment loads were usually strongly correlated with higher TSS levels, except in areas where a large extent of wetlands retained sediment discharged during rainfall events. The modeling results indicated improved water clarity in the shallow aquatic regions of Mississippi Sound and degraded water clarity in the Wolf Bay estuary. This integrated modeling approach will provide new knowledge and tools for coastal resource managers to manage shallow aquatic habitats that provide critical ecosystem services.     

FORECASTING FUTURE LAND USE FOR WATERSHED ASSESSMENT1

ABSTRACT: Protecting surface water quality in watersheds undergoing demographic change requires both the management of existing threats and planning to address potential future stresses arising from changing land use. Many reservoirs and threatened waterbodies are located in areas undergoing rapid population growth, and increases in density of residential and commercial land use, accompanied by increased amount of impervious surface area, can result in increased pollutant loading and degradation of water quality. Effective planning to address potential threats, including zoning and growth management, requires analytical tools to predict and compare the impacts of different management options. The focus of this paper is not on developing demographic projections, but rather the translation of such projections into changes in land use which form the basis for assessment of future watershed loads. Land use change can be forecast at a variety of spatial and temporal scales. A semi-lumped, GIS-based, transition matrix approach is recommended as consistent with the level of complexity achievable in most watershed models. Practical aspects of forecasting future land use for watershed assessment are discussed. Several recent reservoir water supply projection studies are used to demonstrate a general framework for simulating changes in land use and resulting impacts on water quality. In addition to providing a technical basis for selecting optimal management alternatives, such a tool is invaluable for demonstrating to different stakeholder groups the trade-offs among management alternatives, both in terms of water quality and future land use patterns within the watershed.     

The Relative Importance of Road Density and Physical Watershed Features in Determining Coastal Marsh Water Quality in Georgian Bay

We used a GIS-based approach to examine the influence of road density and physical watershed features (watershed size, wetland cover, and bedrock type) on water quality in coastal marshes of Georgian Bay, Ontario. We created a GIS that included landscape information and water-quality data from a 9-year synoptic survey of 105 coastal marshes covering 28 quaternary watersheds. Multiple regressions and partial correlations were used to discern confounding effects of human-induced (road density) versus natural physical watershed determinants of water quality. Road density was the dominant factor influencing many water quality variables, showing positive correlations with specific conductivity (COND), total suspended solids (TSS), and inorganic suspended solids (ISS) and a negative correlation with overall Water Quality Index scores. Road density also showed positive correlations with total nitrate nitrogen (TNN) and total phosphorus (TP). By comparison, larger watershed area was the main factor leading to elevated TP concentrations. The proportion of the watershed occupied by wetlands explained the largest amount of variation in TNN concentrations (negative correlation) and was also negatively correlated with COND and positively correlated with TSS and ISS when we controlled for road density. Bedrock type did not have a significant effect in any of the models. Our findings suggest that road density is currently the overriding factor governing water quality of coastal marshes in Georgian Bay during the summer low-flow period. We recommend that natural variation in physical watershed characteristics be considered when developing water quality standards and management practices for freshwater coastal areas.     

Impacts of Climate Variability and Land Use Alterations on Frequency Distributions of Terrestrial Runoff Loading to Coastal Waters in Southern California1

Abstract: The transport of water, sediment, dissolved and particulate chemicals, and bacteria from coastal watersheds affects the nearshore marine and estuarine waters. In southern California, coastal watersheds deliver water and associated constituents to the nearshore system in discrete pulses. To better understand the pulsed nature of these watersheds, frequency distributions of simulated runoff events are presented for: (1) three land use conditions (1929, 1998, 2050); (2) three time periods (all water years 1989‐2002), only El Nino years (1992, 1993, 1995, 1998); and only non‐El Nino years; and (3) three regions (watershed, uplands, and lowlands). At the watershed scale, there was a significant increase (>200%) in mean event runoff from 1929 to 2050 (0.4‐1.3 cm) due to localized urbanization, which shifted the dominant sources of runoff from the mountains in 1929 (78% of watershed runoff) to the coastal plane for 2050 conditions (51% of watershed runoff). Inter‐annual climate variability was strong in the rainfall and runoff frequency distributions, with mean event rainfall and runoff 66 and 60% larger in El Nino relative to non‐El Nino years. Combining urbanization and climate variability, 2050 land conditions resulted in El Nino years being five times more likely to produce large (>3.0 cm) runoff events relative to non‐El Nino years. Combining frequency distributions of event runoff with regional nutrient export relationships, we show that in El Nino years, one in five events produced runoff ≥2.5 cm and temporary nearshore nitrate and phosphate concentrations of 12 and 1.4 μM, respectively, or approximately 5‐10 times above ambient conditions.     

Effects of Changing Environments of Mangrove Creeks on Fish Communities at Trat Bay, Thailand

Effects of changing environments of riverside mangroves, coastal land uses, and water quality on fish communities were studied in Bangphra and Thaprik creeks, Trat Bay, Thailand. Regression analysis revealed that fish species richness in the wet season had a negative relationship with water transparency, nitrate, and phosphate and a positive relationship with zooplankton. In the dry season, species richness had a negative relationship with nitrate and phosphate and a positive relationship with salinity, pH, and zooplankton. Abundances and species richness of fish declined over distance from downstream to upstream in both creeks. Riparian mangroves and water quality also declined with distance upstream in both creeks. Results from one-way analysis of variance and Tukey’s HSD test revealed that the highest zooplankton volume with the lowest amounts of nitrate and phosphate were observed at the downstream station in both creeks in each season. Low zooplankton volume with high amounts of nitrate and phosphate were found at the midstream and upstream stations of the creeks. The midstream and upstream stations of Bangphra Creek had low to moderate abundance of mangroves along the riversides, whereas shrimp farms were mainly found along the riversides at the midstream and upstream stations of Thaprik Creek. Correlation analysis results of land-use types and the significant habitat factors were discussed. This study found that mangrove degradation, shrimp farming, and residential and agricultural areas altered water quality and the health of fish habitats, causing the decreases in fish abundance and species richness.     

Modeling Nutrient Export From Coastal California Watersheds

Nitrate and phosphate export coefficient models were developed for coastal watersheds along the Santa Barbara Channel in central California. One approach was based on measurements of nutrient fluxes in streams from specific land use classes and included a watershed response function that scaled export up or down depending on antecedent moisture conditions. The second approach for nutrient export coefficient modeling used anthropogenic nutrient loading for land use classes and atmospheric nutrient deposition to model export. In an application of the first approach to one watershed, the nitrate and phosphate models were within 20% of measured values for most storms. When applied to another year, both nitrate and phosphate models generally performed adequately with annual, storm‐flow, and base‐flow values within 20% of measured nutrient loadings. Less satisfactory results were found when applied to neighboring watersheds with difference percentages of land use and hydrologic conditions. Application of the second approach was less successful than the first approach.     

Assessment of Economic and Water Quality Impacts of Land Use Change Using a Simple Bioeconomic Model

The objective of this study is to assess the economic and water quality impact of land use change in a small watershed in the Wiregrass region of Alabama. The study compares changes in water quality and revenue from agricultural and timber production due to changes in land use between years 1992 and 2001. The study was completed in two stages. In the first stage, a biophysical model was used to estimate the effect of land use change on nitrogen and phosphorus runoff and sediment deposition in the main channel; in the second stage, farm enterprise budgeting tools were used to estimate the economic returns for the changes in land use condition. Both biophysical and economic results are discussed, and a case for complex optimization to develop a decision support system is presented.     

TMDL Balance: A Model for Coastal Water Pollutant Loadings

Bacterial contamination accounts for more than 60% of the impairments included on the 2008 Texas 303(d) List. Many of these bacterial impairments are along the Texas Gulf Coast because coastal waters often are regulated for oyster harvesting, which have strict water quality standards. Under the Clean Water Act, each one of these impaired waterbodies requires a total maximum daily load (TMDL) study to be performed. A recent, statewide study recommended the development and application of simple modeling approaches to address the majority of Texas's bacteria TMDLs, including “… simple load duration curve, GIS [geographic information systems], and/or mass balance models.” We developed the TMDL Balance model in response to this recommendation. TMDL Balance is a steady state, mass balance, GIS‐based model for simulating pollutant loads and concentrations in coastal systems. The model uses plug‐flow reactor and continuously‐stirred tank reactor equations to route spatially distributed point and nonpoint source loads through a watershed via overland flow, non‐tidal flow, and tidal flow, decaying the loads via first‐order kinetics. In this paper, we explain the development of the watershed loading portion of the TMDL Balance model, demonstrating the methodology through a case study: computing bacterial loads in the Copano Bay watershed of southeast Texas. The application highlights an example of distributing bacterial sources spatially based on land use data.     

MODELING DEVELOPED COASTAL WATERSHEDS WITH THE AGRICULTURAL NON-POINT SOURCE MODEL1

ABSTRACT: Many coastal states are facing increasing urban growth along their coast lines. The growth has caused urban non-point source nitrogen runoff to be a major contributor to coastal and estuarine enrichment. Water resource managers are responsible for evaluating the impacts from point and non-point sources in developed watersheds and developing strategies to manage future growth. Non-point source models provide an effective approach to these management challenges. The Agricultural Non-Point Source Model (AGNPS) permits the incorporation of important spatial information (soils, landuse, topography, hydrology) in simulating surface hydrology and nitrogen non-point source runoff. The AGNPS model was adapted for developed coastal watersheds by deriving urban coefficients that reflect urban landuse classes and the amount of impervious surface area. Popperdam Creek watershed was used for model parameter development and model calibration. Four additional watersheds were simulated to validate the model. The model predictions of the peak flow and total nitrogen concentrations were close to the field measurements for the five sub-basins simulated. Measured peak flow varied by 30 fold among the sub-basins. The average simulated peak flow was within 14 percent of the average measured peak flow. Measured total nitrogen loads varied over an order of magnitude among the sub-basins yet error between the measured and simulated loads for a given sub-basin averaged 5 percent. The AGNPS model provided better estimates of nitrogen loads than widely used regression methods. The spatial distribution of important watershed characteristics influenced the impacts of urban landuse and projecting future residential expansion on runoff, sediment and nitrogen yields. The AGNPS model provides a useful tool to incorporate these characteristics, evaluate their importance, and evaluate fieldscale to watershed-scale urban impacts.     

Modeling the relationship between land use and surface water quality

It is widely known that watershed hydrology is dependent on many factors, including land use, climate, and soil conditions. But the relative impacts of different types of land use on the surface water are yet to be ascertained and quantified. This research attempted to use a comprehensive approach to examine the hydrologic effects of land use at both a regional and a local scale. Statistical and spatial analyses were employed to examine the statistical and spatial relationships of land use and the flow and water quality in receiving waters on a regional scale in the State of Ohio. Besides, a widely accepted watershed-based water quality assessment tool, the Better Assessment Science Integrating Point and Nonpoint Sources (BASINS), was adopted to model the plausible effects of land use on water quality in a local watershed in the East Fork Little Miami River Basin. The results from the statistical analyses revealed that there was a significant relationship between land use and in-stream water quality, especially for nitrogen, phosphorus and Fecal coliform. The geographic information systems (GIS) spatial analyses identified the watersheds that have high levels of contaminants and percentages of agricultural and urban lands. Furthermore, the hydrologic and water quality modeling showed that agricultural and impervious urban lands produced a much higher level of nitrogen and phosphorus than other land surfaces. From this research, it seems that the approach adopted in this study is comprehensive, covering both the regional and local scales. It also reveals that BASINS is a very useful and reliable tool, capable of characterizing the flow and water quality conditions for the study area under different watershed scales. With little modification, these models should be able to adapt to other watersheds or to simulate other contaminants. They also can be used to study the plausible impacts of global environmental change. In addition, the information on the hydrologic effects of land use is very useful. It can provide guidelines not only for resource managers in restoring our aquatic ecosystems, but also for local planners in devising viable and ecologically-sound watershed development plans, as well as for policy makers in evaluating alternate land management decisions.     

Predicting estuarine sediment metal concentrations and inferred ecological conditions: an information theoretic approach

Empirically derived relationships associating sediment metal concentrations with degraded ecological conditions provide important information to assess estuarine condition. Resources limit the number, magnitude, and frequency of monitoring activities to acquire these data. Models that use available information and simple statistical relationships to predict sediment metal concentrations could provide an important tool for environmental assessment. We developed 45 predictive models for the total concentrations of copper, lead, mercury, and cadmium in estuarine sediments along the Southern New England and Mid-Atlantic regions of the United States. Using information theoretic model-averaging approaches, we found total developed land and percent silt/clay of estuarine sediment were the most important variables for predicting the presence of all four metals. Estuary area, river flow, tidal range, and total agricultural land varied in their importance. The model-averaged predictions explained 78.4, 70.5, 56.4, and 50.3% of the variation for copper, lead, mercury, and cadmium, respectively. Overall prediction accuracies of selected sediment benchmark values (i.e., effects ranges) were 83.9, 84.8, 78.6, and 92.0% for copper, lead, mercury, and cadmium, respectively. Our results further support the generally accepted conclusion that sediment metal concentrations are best described by the physical characteristics of the estuarine sediment and the total amount of urban land in the contributing watershed. We demonstrated that broad-scale predictive models built from existing monitoring data with information theoretic model-averaging approaches provide valuable predictions of estuarine sediment metal concentrations and show promise for future environmental modeling efforts in other regions.     

Effect of long-term application of biosolids for land reclamation on surface water chemistry

Biosolids are known to have a potential to restore degraded land, but the long-term impacts of this practice on the environment, including water quality, still need to be evaluated. The surface water chemistry (NO3-, NH4+, and total P, Cd, Cu, and Hg) was monitored for 31 yr from 1972 to 2002 in a 6000-ha watershed at Fulton County, Illinois, where the Metropolitan Water Reclamation District of Greater Chicago was restoring the productivity of strip-mined land using biosolids. The mean cumulative loading rates during the past 31 yr were 875 dry Mg ha(-1) for 1120-ha fields in the biosolids-amended watershed and 4.3 dry Mg ha(-1) for the 670-ha fields in the control watershed. Biosolids were injected into mine spoil fields as liquid fertilizer from 1972 to 1985, and incorporated as dewatered cake from 1980 to 1996 and air-dried solids from 1987 to 2002. The mean annual loadings of nutrients and trace elements from biosolids in 1 ha were 735 kg N, 530 kg P, 4.5 kg Cd, 30.7 kg Cu, and 0.11 kg Hg in the fields of the biosolids-amended watershed, and negligible in the fields of the control watershed. Sampling of surface water was conducted monthly in the 1970s, and three times per year in the 1980s and 1990s. The water samples were collected from 12 reservoirs and 2 creeks receiving drainage from the fields in the control watershed, and 8 reservoirs and 4 creeks associated with the fields in the biosolids-amended watershed for the analysis of NO3- -N (including NO2- N), NH4+-N, and total P, Cd, Cu, and Hg. Compared to the control (0.18 mg L(-1)), surface water NO3- -N in the biosolids-amended watershed (2.23 mg L(-1)) was consistently higher; however, it was still below the Illinois limit of 10 mg L(-1) for public and food-processing water supplies. Biosolids applications had a significant effect on mean concentrations of ammonium N (0.11 mg L(-1) for control and 0.24 mg L(-1) for biosolids) and total P (0.10 mg L(-1) for control and 0.16 mg L(-1) for biosolids) in surface water. Application of biosolids did not increase the concentrations of Cd and Hg in surface water. The elevation of Cu in surface water with biosolids application only occurred in some years of the first decade, when land-applied sludges contained high concentrations of trace metals, including Cu. In fact, following the promulgation of 40 CFR Part 503, the concentrations of all three metals fell below the method detection level (MDL) in surface water for nearly all samplings. Nitrate in the surface water tends to be higher in spring, and ammonium, total P, and total Hg in summer and fall. Mean nitrate, ammonium, and total phosphorus concentrations were found to be greater in creeks than reservoirs. The results indicate that application of biosolids for land reclamation at high loading rates from 1972 to 2002, with adequate runoff and soil erosion control, had only a minor impact on surface water quality.     

LAND TREATMENT IN WATERSHED PROGRAM AREAS1

Several types of watershed program areas, principally but not solely those operated under authority of the Watershed Protection and Flood Prevention Act, are greatly altering the American landscape for good. While structural measures tend to be the most spectacular feature of these watersheds, land treatment measures (soil, water, woodland, and wildlife conservation practices) actually are basic to the long term stability and maximum contributions of these watershed programs.     

WATERSHED SCALING EFFECT ON BASE FLOW NITRATE,VALLEY AND RIDGE PHYSIOGRAPHIC PROVINCE1

ABSTRACT: A study of stream base flow and NO₃‐N concentration was conducted simultaneously in 51 subwatersheds within the 116‐square‐kilometer watershed of East Mahantango Creek near Klingerstown, Pennsylvania. The study was designed to test whether measurable results of processes and observations within the smaller watersheds were similar to or transferable to a larger scale. Ancillary data on land use were available for the small and large watersheds. Although the source of land‐use data was different for the small and large watersheds, comparisons showed that the differences in the two land‐use data sources were minimal. A land use‐based water‐quality model developed for the small‐scale 7.3‐square‐kilometer watershed for a previous study accurately predicted NO₃‐N concentrations from sampling in the same watershed. The water‐quality model was modified and, using the imagery‐based land use, was found to accurately predict NO₃‐N concentrations in the subwatersheds of the large‐scale 116‐square‐kilometer watershed as well. Because the model accurately predicts NO₃‐N concentrations at small and large scales, it is likely that in second‐order streams and higher, discharge of water and NO₃‐N is dominated by flow from smaller first‐order streams, and the contribution of ground‐water discharge to higher order streams is minimal at the large scale.     

Modelling hydroenvironmental and health risk assessment parameters along the South Wales Coast

This paper highlights the increasing concerns relating to hydroenvironmetal issues and cites recent examples of the challenges now being regularly faced by hydroenvironmetal scientists and engineers. The limitations and restrictions of both physical (or laboratory) and numerical (or computer based) hydraulic models used in the planning and management of aquatic basins are discussed. General details are given of numerical models used for flow and water quality concentration predictions in estuarine waters, with particular application to the challenges occurring along the South Wales coast. A highly accurate and non-diffusive finite difference scheme that solves the transport equation for predicting water quality indicators and suspended sediment concentration distributions is also discussed. In particular, details are outlined of the extension of the water quality indicators of faecal coliforms, as required to comply with the EU Bathing Water Directive, to predict health risk assessment, in the form of predicting the risk of gastroenteritis. Three example research projects along the South Wales coast are described; the projects involve the application of two-dimensional and three-dimensional hydroenvironmetal models to predict flow patterns and water quality indicator organism distributions in the coastal receiving waters. These studies include: (i) a curvilinear finite difference approach to modelling flows in the Bristol Channel, (ii) coastal health risk predictions in Swansea Bay using combined water quality and epidemiological models, and (iii) combined sewer overflow discharges into Cardiff Bay.     

Development of bacteria and benthic total maximum daily loads: a case study, Linville Creek, Virginia

Two total maximum daily load (TMDL) studies were performed for Linville Creek in Rockingham County, Virginia, to address bacterial and benthic impairments. The TMDL program is an integrated watershed management approach required by the Clean Water Act. This paper describes the procedures used by the Center for TMDL and Watershed Studies at Virginia Tech to develop the Linville Creek TMDLs and discusses the key lessons learned from and the ramifications of the procedures used in these and other similar TMDL studies. The bacterial impairment TMDL was developed using the Hydrological Simulation Program-Fortran (HSPF). Fecal coliform loads were estimated through an intensive source characterization process. The benthic impairment TMDL was developed using the Generalized Watershed Loading Function (GWLF) model and the reference watershed approach. The bacterial TMDL allocation scenario requires a 100% reduction in cattle manure direct-deposits to the stream, a 96% reduction in nonpoint-source loadings to the land surface, and a 95% reduction in wildlife direct-deposits to the stream. Sediment was identified as the primary benthic stressor. The TMDL allocation scenario for the benthic impairment requires an overall reduction of 12.3% of the existing sediment loads. Despite the many drawbacks associated with using watershed-scale models like HSPF and GWLF to develop TMDLs, the detailed watershed and pollutant-source characterization required to use these and similar models creates information that stakeholders need to select appropriate corrective measures to address the cause of the water quality impairment when implementing the TMDL.