MINIMIZING THE IMPACT OF URBANIZATION ON LONG TERM RUNOFF1

Increasing concern about the problems caused by urban sprawl has encouraged development and implementation of smart growth approaches to land use management. One of the goals of smart growth is water resources protection, in particular minimizing the runoff impact of urbanization. To investigate the magnitude of the potential benefits of land use planning for water resources protection, possible runoff impacts of historical and projected urbanization were estimated for two watersheds in Indiana and Michigan using a long term hydrological impact analysis model. An optimization component allowed selection of land use change placements that minimize runoff increase. Optimizing land use change placement would have reduced runoff increase by as much as 4.9 percent from 1973 to 1997 in the Indiana study watershed. For nonsprawl and sprawl scenarios in the Michigan watershed for 1978 to 2040, optimizing land use change placement would have reduced runoff increase by 12.3 percent and 20.5 percent, respectively. The work presented here illustrates both an approach to assessing the magnitude of the impact of smart growth and the significant potential scale of smart growth in moderating runoff changes that result from urbanization. The results of this study have significant implications for urban planning.     

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.     

Assessing Watershed-Scale, Long-Term Hydrologic Impacts of Land-Use Change Using a GIS-NPS Model

Land-use change, dominated by an increase in urban/impervious areas, has a significant impact on water resources. This includes impacts on nonpoint source (NPS) pollution, which is the leading cause of degraded water quality in the United States. Traditional hydrologic models focus on estimating peak discharges and NPS pollution from high-magnitude, episodic storms and successfully address short-term, local-scale surface water management issues. However, runoff from small, low-frequency storms dominates long-term hydrologic impacts, and existing hydrologic models are usually of limited use in assessing the long-term impacts of land-use change. A long-term hydrologic impact assessment (L-THIA) model has been developed using the curve number (CN) method. Long-term climatic records are used in combination with soils and land-use information to calculate average annual runoff and NPS pollution at a watershed scale. The model is linked to a geographic information system (GIS) for convenient generation and management of model input and output data, and advanced visualization of model results. The L-THIA/NPS GIS model was applied to the Little Eagle Creek (LEC) watershed near Indianapolis, Indiana, USA. Historical land-use scenarios for 1973, 1984, and 1991 were analyzed to track land-use change in the watershed and to assess impacts on annual average runoff and NPS pollution from the watershed and its five subbasins. For the entire watershed between 1973 and 1991, an 18% increase in urban or impervious areas resulted in an estimated 80% increase in annual average runoff volume and estimated increases of more than 50% in annual average loads for lead, copper, and zinc. Estimated nutrient (nitrogen and phosphorus) loads decreased by 15% mainly because of loss of agricultural areas. The L-THIA/NPS GIS model is a powerful tool for identifying environmentally sensitive areas in terms of NPS pollution potential and for evaluating alternative land use scenarios for NPS pollution management.     

MOBILITY AND AQUATIC TOXICITY OF COPPER IN AN URBAN WATERSHED1

ABSTRACT: Abundant use of copper based products has resulted in increased violation of copper water quality criteria in runoff from urban storm water systems. The objectives of this work were to understand the mobility and toxicity of copper in an urban watershed and to apportion the amount of copper entering the freshwater receiving stream from different urban land covers using a mass balance approach. Sixteen rainfall events collected from the University of Connecticut study watershed between August 1998 and September 2000 were analyzed to assess copper flux in an urban storm water system. Mean flow weighted dissolved copper concentrations observed in the study for copper based architectural material runoff, pervious area runoff, impervious area runoff, and in the receiving stream were 1210 ± 840, 9 ± 3, 8 ± 2, and 14 ± 7 μg/L, respectively. Mean dissolved copper concentrations in the receiving stream exceeded Connecticut's water quality criteria. Despite exceeding the dissolved concentration based criteria, cupric ion concentrations at the system outlet remained below 0.05 μg/L for all storms analyzed, and no acute toxicity (using Daphnia pulex as the test organism) was measured in samples collected from the stream.     

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.     

Assessing Surface Water Quality and Its Relation with Urban Land Cover Changes in the Lake Calumet Area, Greater Chicago

Urban land use and land cover change significantly affect spatial and temporal patterns of runoff, which in turn impacts surface water quality. With the exponential growth in urban areas over the past three decades, changes in land use and land cover to cater for the growth of cities has been a conspicuous spectacle in urban spaces. The main goal of this study was to assess the impacts of land cover change on runoff and surface water quality using a partial area hydrology framework. The study employed ArcHydro GIS extension and a modified version of Long-Term Hydrologic and Nonpoint Source Pollution model (L-THIA-NPS) in estimating runoff and nonpoint source pollutant concentration around Lake Calumet between 1992 and 2001. Data employed include National Land Cover Data set, rainfall data, digital elevation model (DEM), Soil Survey Geographic (SSURGO) data, and The United States Environmental Protection Agency’s STORET (storage and retrieval) water quality data. The model was able to predict surface water quality reasonably well over the study period. Sensitivity analysis facilitated a manual calibration of the model. Model validation was executed by comparing simulated results following calibration and observed water quality data for the study area. The study demonstrates that the level of concentration of nonpoint source pollutants in surface water within an urban watershed heavily depends on the spatiotemporal variations in areas that contribute towards runoff compared to the spatial extent of change in major land use/land cover.     

EFFECTS OF PRECIPITATION AND LAND USE ON STORM RUNOFF1

ABSTRACT: Storm-runoff quantity and quality were studied in three watersheds located near St. Paul in Ramsey County, Minnesota, from April 15 through September 15 of 1984, 1985, and 1986 to qualitatively determine the effects of precipitation and selected land uses on storm runoff. In respect to precipitation effects, differences in stormrunoff quantity between years in an urban watershed that lacks wetlands appear to be related to the average storm size (amount of precipitation) during the study period of each year. In contrast, the differences in storm-runoff quantity from watersheds that contain wetlands appear to be related to total precipitation during study period of each year. In respect to land use, the differences in storm-runoff quantity appear to be related to the amounts of impervious and wetland area. The watershed that contains the largest amount of impervious area and smallest amount of wetland area has the largest amount of storm runoff. Differences in storm-runoff quality appear to be related to the amounts of wetland and lake area. The watershed that contains the largest amounts of wetland and lake area has the smallest storm-runoff loading of suspended solids, phosphorus, and nitrogen. The wetland and lake areas likely retain the loading and, subsequently, lower the amount of storm-runoff loading exported from a watershed.     

EXPORT COEFFICIENT MODELING TO ASSESS PHOSPHORUS LOADING IN AN URBAN WATERSHED1

ABSTRACT: An export coefficient modeling approach was used to assess the influence of land use on phosphorus loading to a Southern Ontario stream. A model was constructed for the 1995–1996 water year and calibrated within ± 3 percent of the observed mean concentration of total phosphorus. It was found that runoff from urban areas contributed most to the loading of phosphorus to the stream. When the model was assessed by running it for the 1977–1978 water year, using water quality and land use data collected independently, agreement within ± 7 percent was obtained. The model was then used to forecast the impact of future urban development proposed for the watershed, in terms of phosphorus loading, and to evaluate the reduction in loading resulting from several urban best management practices (BMP). It was determined that phosphorus removal will have to be applied to all the urban runoff from the watershed to appreciably reduce stream phosphorus concentration. Of the BMP designs assessed, an infiltration pond system resulted in the greatest phosphorus load reduction, 50 percent from the 1995–1996 baseline.     

STORM RUNOFF PREDICTION BASED ON A SPATIALLY DISTRIBUTED TRAVEL TIME METHOD UTILIZING REMOTE SENSING AND GIS1

ABSTRACT: In this study, remotely sensed data and geographic information system (GIS) tools were used to estimate storm runoff response for Simms Creek watershed in the Etonia basin in northeast Florida. Land cover information from digital orthophoto quarter quadrangles (DOQQ), and enhanced thematic mapper plus (ETM+) were analyzed for the years 1990, 1995, and 2000. The corresponding infiltration excess runoff response of the study area was estimated using the U.S. Department of Agriculture (USDA), Natural Resources Conservation Service Curve Number (NRCS‐CN) method. A digital elevation model (DEM)/GIS technique was developed to predict stream response to runoff events based on the travel time from each grid cell to the watershed outlet. A comparison of predicted to observed stream response shows that the model predicts the total runoff volume with an efficiency of 0.98, the peak flow rate at an efficiency of 0.85, and the full direct runoff hydrograph with an average efficiency of 0.65. The DEM/GIS travel time model can be used to predict the runoff response of ungaged watersheds and is useful for predicting runoff hydrographs resulting from proposed large scale changes in the land use.     

An analysis of urban development and its environmental impact on the Tampa Bay watershed

Urbanization has transformed natural landscapes into anthropogenic impervious surfaces. Urban land use has become a major driving force for land cover and land use change in the Tampa Bay watershed of west-central Florida. This study investigates urban land use change and its impact on the watershed. The spatial and temporal changes, as well as the development density of urban land use are determined by analyzing the impervious surface distribution using Landsat satellite imagery. Population distribution and density are extracted from the 2000 census data. Non-point source pollution parameters used for measuring water quality are analyzed for the sub-drainage basins of Hillsborough County. The relationships between 2002 urban land use, population distribution and their environmental influences are explored using regression analysis against various non-point source pollutant loadings in these sub-drainage basins. The results suggest that strong associations existed between most pollutant loadings and the extent of impervious surface within each sub-drainage basin in 2002. Population density also exhibits apparent correlations with loading rates of several pollutants. Spatial variations of selected non-point source pollutant loadings are also assessed.     

Modeling Agricultural Nonpoint Source Pollution Using a Geographic Information System Approach

Agricultural non-point source (NPS) pollution, primarily sediment and nutrients, is the leading source of water-quality impacts to surface waters in North America. The overall goal of this study was to develop geographic information system (GIS) protocols to facilitate the spatial and temporal modeling of changes in soils, hydrology, and land-cover change at the watershed scale. In the first part of this article, we describe the use of GIS to spatially integrate watershed scale data on soil erodibility, land use, and runoff for the assessment of potential source areas within an intensively agricultural watershed. The agricultural non-point source pollution (AGNPS) model was used in the Muddy Creek, Ontario, watershed to evaluate the effectiveness of management strategies in decreasing sediment and nutrient [phosphorus (P)] pollution. This analysis was accompanied by the measurement of water-quality parameters (dissolved oxygen, pH, hardness, alkalinity, and turbidity) as well as sediment and P loadings to the creek. Practices aimed at increasing year-round soil cover would be most effective in decreasing sediment and P losses in this watershed. In the second part of this article, we describe a method for characterizing land-cover change in a dynamic urban fringe watershed. The GIS method we developed for the Blackberry Creek, Illinois, watershed will allow us to better account for temporal changes in land use, specifically corn and soybean cover, on an annual basis and to improve on the modeling of watershed processes shown for the Muddy Creek watershed. Our model can be used at different levels of planning with minimal data preprocessing, easily accessible data, and adjustable output scales.     

Estimating the Variability and Confidence of Land Use and Imperviousness Relationships at a Regional Scale1

Abstract: Impervious cover is a commonly used metric to help explain or predict anthropogenic impacts on aquatic resources; often it is used as a surrogate for intensity of human impacts when evaluating effects on aquatic resources. The most common way to estimate imperviousness is based on relationships with land use. Few studies have evaluated how the relationship between impervious surface and land use varies among geographies with different levels of development and between types of imagery used to assign land use type. In this study, we assess variability in estimates of imperviousness based on two locally available land use datasets: one based on aerial imagery (2‐m resolution) and another based on satellite imagery (30‐m resolution). The ranges and variability in imperviousness within land use categories were assessed at several spatial scales, including within counties, between counties, and between watersheds. Results indicate that there was considerable variability for all developed land use types. Estimated impervious cover often varied over a range of 20‐40% points within a land use category. Furthermore, there were clear spatial patterns both between and within counties, with impervious cover for a given land use type being higher near the urban centers and lower at the margins of development. Estimates of imperviousness for 12 study watersheds indicated that variability increased with increasing watershed development, making it difficult to confidently set management or regulatory targets based on impervious cover. This study suggests that locally derived, high resolution satellite or aerial imagery should be used to estimate imperviousness when a high level of accuracy and precision is required for regulatory or management decisions. Furthermore, the error associated with impervious land use relationships should be accounted for when using impervious cover in runoff or water quality models, or when making management decisions regarding stream health.     

MODELING FRAMEWORK FOR MANAGING COPPER RUNOFF IN URBAN WATERSHEDS1

ABSTRACT: A modeling framework was developed for managing copper runoff in urban watersheds that incorporates water quality characterization, watershed land use areas, hydrologic data, a statistical simulator, a biotic ligand binding model to characterize acute toxicity, and a statistical method for setting a watershed specific copper loading. The modeling framework is driven by export coefficients derived from water quality parameters and hydrologic inputs measured in an urban watershed's storm water system. This framework was applied to a watershed containing a copper roof built in 1992. A series of simulations was run to predict the change in receiving stream water chemistry caused by roof aging and to determine the maximum copper loading (at the 99 percent confidence level) a watershed could accept without causing acute toxicity in the receiving stream. Forecasting the amount of copper flux responsible for exceeding the assimilation capacity of a watershed can be directly related to maximum copper loadings responsible for causing toxicity in the receiving streams. The framework developed in this study can be used to evaluate copper utilization in urban watersheds.     

SIMULATING HYDROLOGIC AND WATER QUALITY IMPACTS IN AN URBANIZING WATERSHED1

ABSTRACT: The Hydrologic Simulation Program‐Fortran (HSPF) was calibrated and used to assess the future effects of various land development scenarios on water quality in the Polecat Creek watershed in Caroline County, Virginia. Model parameters related to hydrology and water quality were calibrated and validated using observed stream flow and water quality data collected at the watershed outlet and the outlets of two sub water sheds. Using the county's Comprehensive Plan, land use scenarios were developed by taking into account the trends and spatial distributions of future development. The simulation results for the various land use scenarios indicate that runoff volume and peak rate increased as urban areas increased. Urbanization also increased sediment loads mainly due to increases in channel erosion. Constituent loads of total Kjeldal nitrogen, orthophosphorus, and total phosphorous for Polecat Creek watershed slightly decreased under future development scenarios. These reductions are due to increases in urban areas that typically contribute smaller quantities of nitrogen and phosphorous, as compared to agricultural areas. However, nitrate loads increased for the future land use scenarios, as compared to the existing land use. The increases in nitrate loads may result from increases in residential land and associated fertilizer use and concurrent decreases in forested land. The procedures used in this paper could assist local, state, and regional policy makers in developing land management strategies that minimize environmental impacts while allowing for future development.     

STREAM HEALTH RANKINGS PREDICTED BY SATELLITE DERIVED LAND COVER METRICS1

ABSTRACT: Land cover and land use change have long been known to influence the chemical, physical, and biological characteristics of streams. This study makes use of land cover maps derived from fine resolution satellite imagery and an extensive stream quality dataset to determine the relationship between small watershed health rankings and land cover composition and configuration. Landscape metrics were derived from digital impervious surface area (ISA), tree cover (percent), and agricultural crop maps within Montgomery County, Maryland. Watershed rankings were developed by state and county collaborators (MD‐DNR and MCDEP) using extensive biological and chemical measurements. In stepwise logistic regression models the factors accounting for the most variation in stream health ranking were the percent ISA, followed by the percent of tree cover. Riparian buffer zone tree cover was also a significant predictor. Of the metrics that considered the spatial configuration of the landscape, a contagion index and the percent of ISA in the flow path from the ISA to the stream were also found to be significant predictors of stream health. Despite limited ability to characterize landscape configuration or narrow riparian buffer zone vegetation with coarser resolution imagery (from Landsat), model results were not significantly different from those based on the use of fine‐resolution ISA information, suggesting that broader area applications of the approach are possible. The results indicate that management practices designed to improve stream water quality should focus on the amount of ISA and tree cover in both the watershed and within the buffer zone.     

Measuring the Impact of Urbanization on Channel Widths Using Historic Aerial Photographs and Modern Surveys1

Abstract: Land use in a watershed is commonly held to exert a strong influence on trunk channel form and process. Land use changes act over human time‐scales, which are short enough to measure effects on channels directly using historic aerial photographs. We show that high‐resolution topographic surveys for the channels of paired watersheds in the Lehigh Valley, Pennsylvania, are comparable, but have channel widths that have changed dramatically in the past five decades. The two watersheds, Little Lehigh Creek and Sacony Creek, are similar in most aspects except in their respective amount of urban land use. Aerial photographs of the urbanized Little Lehigh Creek show that a majority of the measured widths (67 of 85) were statistically wider in 1999 than in 1947. In contrast, the measured widths from the agricultural Sacony Creek are more evenly distributed among those that widened (18), narrowed (28), and those that were statistically unchanged (6) from 1946 to 1999. From 1946 to 1999 the only section of Sacony Creek that widened was that reach downstream of the only sizable urban area in the watershed. The current land use in Sacony Creek watershed resembles that of 1946, while the Little Lehigh Creek watershed has more than tripled its urban area. These data, in concert with other recent hydrologic data from the watersheds suggest that the increase in urban area‐generated peak discharges is the mechanism behind the widening that occurred in the Little Lehigh Creek. These wider channels can affect water quality, aquatic habitat, suspended sediment loads, and river esthetics.     

Weighting Nitrogen and Phosphorus Pixel Pollutant Loads to Represent Runoff and Buffering Likelihoods

Watershed models often estimate annual nitrogen (N) or phosphorus (P) pollutant loads in rural areas with export coefficient (EC) (kg/ha/yr) values based on land cover, and in urban areas as the product of spatially uniform event mean concentration (EMC) (mg/L) values and runoff volume. Actual N and P nonpoint source (NPS) pollutant loading has more spatial complexity due to watershed variation in runoff likelihood and buffering likelihood along surface and subsurface pathways, which can be represented in a contributing area dispersal area (CADA) NPS model. This research develops a CADA NPS model to simulate how watershed properties of elevation, land cover, and soils upslope and downslope of each watershed pixel influence nutrient loading. The model uses both surface and subsurface runoff indices (RI), and surface and subsurface buffer indices (BI), to quantify the runoff and buffering likelihood for each watershed pixel, and generate maps of weighted EC and EMC values that identify NPS pollutant loading hotspots. The research illustrates how CADA NPS model maps and pixel loading values are sensitive to the spatial resolution and accuracy of elevation and land cover data, and model predictions can represent the lower and upper bounds of NPS loading. The model provides managers with a tool to rapidly visualize, rank, and investigate likely areas of high nutrient export.     

CHANGING RAINFALL-RUNOFF RELATIONSHIPS IN THE URBANIZING PEACHTREE CREEK WATERSHED,ATLANTA, GEORGIA1

ABSTRACT: Peachtree Creek is a gaged watershed that has experienced a substantial increase in urbanization. The relationships of runoff to rainfall were studied for total annual flows, low flows, and peak flows. For each type of flow the relationship in the later, more urbanized period was compared to that in the earlier, less urbanized period. An increase in total runoff in wet years was observed as urbanization increased, but a decrease occurred during dry years. For low flows a similar decrease of runoff in dry years was found. An increase in peak runoff was observed over most of the range of precipitation. Increasing peak flows and declining low flows can be adequately explained by urban hydrologic theoryshed. which focuses on the effects of urban impervious surfaces upon direct runoff and infiltration. However, a decline of total runoff in dry years can be explained only by taking into account evapotranspiration as well. The concept of advectively assisted urban evapotranspiration, previously discovered by climatologists, is needed to explain such a loss of total runoff. Urban hydrologic theory must take into account vegetation and evapotranspiration, as well as impervious surfaces and their direct runoff, to explain the magnitude of total annual flows and low flows. Urban stormwater management should address the restoration of low flows, as well as the control of floods.     

A conceptual framework of agricultural land use planning with BMP for integrated watershed management

Land use planning is an important element of the integrated watershed management approach. It not only influences the environmental processes such as soil and stream bed erosion, sediment and nutrient concentrations in streams, quality of surface and ground waters in a watershed, but also affects social and economic development in that region. Although its importance in achieving sustainable development has long been recognized, a land use planning methodology based on a systems approach involving realistic computational modeling and meta-heuristic optimization is still lacking in the current practice of integrated watershed management. The present study proposes a new approach which attempts to combine computational modeling of upland watershed processes, fluvial processes and modern heuristic optimization techniques to address the water-land use interrelationship in its full complexity. The best land use allocation is decided by a multi-objective function that minimizes sediment yields and nutrient concentrations as well as the total operation/implementation cost, while the water quality and the production benefits from agricultural exploitation are maximized. The proposed optimization strategy considers also the preferences of land owners. The runoff model AnnAGNPS (developed by USDA), and the channel network model CCHE1D (developed by NCCHE), are linked together to simulate sediment/pollutant transport process at watershed scale based on any assigned land use combination. The greedy randomized adaptive Tabu search heuristic is used to flip the land use options for finding an optimum combination of land use allocations. The approach is demonstrated by applying it to a demonstrative case study involving USDA Goodwin Creek experimental watershed located in northern Mississippi. The results show the improvement of the tradeoff between benefits and costs for the watershed, after implementing the proposed optimal land use planning.     

INFLUENCES OF WATERSHED URBANIZATION AND INSTREAM HABITAT ON MACROINVERTEBRATES IN COLD WATER STREAMS1

ABSTRACT: We analyzed data from riffle and snag habitats for 39 small cold water streams with different levels of watershed urbanization in Wisconsin and Minnesota to evaluate the influences of urban land use and instream habitat on macroinvertebrate communities. Multivariate analysis indicated that stream temperature and amount of urban land use in the watersheds were the most influential factors determining macroinvertebrate assemblages. The amount of watershed urbanization was nonlinearly and negatively correlated with percentages of Ephemeroptera‐Plecoptera‐Trichoptera (EPT) abundance, EPT taxa, filterers, and scrapers and positively correlated with Hilsenhoff biotic index. High quality macroinvertebrate index values were possible if effective imperviousness was less than 7 percent of the watershed area. Beyond this level of imperviousness, index values tended to be consistently poor. Land uses in the riparian area were equal or more influential relative to land use elsewhere in the watershed, although riparian area consisted of only a small portion of the entire watershed area. Our study implies that it is extremely important to restrict watershed impervious land use and protect stream riparian areas for reducing human degradation on stream quality in low level urbanizing watersheds. Stream temperature may be one of the major factors through which human activities degrade cold‐water streams, and management efforts that can maintain a natural thermal regime will help preserve stream quality.