Effects of Urbanization on Flow Duration and Stream Flashiness: A Case Study of Puget Sound Streams,Western Washington,USA

The overall influence of urbanization on how flows of different frequency might change over time, while important in hydrologic design, remains imprecisely known. In this study, we investigate the effects of urbanization on flow duration curves (FDCs) and flow variability through a case study of eight watersheds that underwent different amounts of growth, in the Puget Sound region in Western Washington State, United States. We computed annual FDCs from flow records spanning 1960‐2010 and, after accounting for the effects of precipitation, we conducted statistical trend analyses on flow metrics to quantify how key FDC percentiles changed with time in response to urbanization. In the urban watersheds, the entire FDC tended to increase in magnitude of flow, especially the 95th‐99th percentile of the daily mean flow series, which increased by an average of 43%. Stream flashiness in urban watersheds was found to increase by an average of 70%. The increases in FDC magnitude and flashiness in urbanizing watersheds are most likely a result of increasing watershed imperviousness and altered hydrologic routing. Rural watersheds were found to have decreasing FDC magnitude over the same time period, which is possibly due to anthropogenic extractions of groundwater, and increasing stream flashiness, which is likely a result of reductions in base flow and increasing precipitation intensity and variability.     

Flow,Organic, and Inorganic Sediment Yields from a Channelized Watershed in the South Carolina Lower Coastal Plain

Many small streams in coastal watersheds in the southeastern United States are modified for agricultural, residential, and commercial development. In the South Carolina Lower Coastal Plain, low‐relief topography and a shallow water table make stream channelization ubiquitous. To quantify the impacts of urbanization and stream channelization, we measured flow and sediment from an urbanizing watershed and a small forested watershed. Flow and sediment export rates were used to infer specific yields from forested and nonforested regions of the urbanizing watershed. Study objectives were to: (1) quantify the range of runoff‐to‐rainfall ratios; (2) quantify the range of specific sediment yields; (3) characterize the quantity and quality of particulate matter exported; and (4) estimate sediment yield attributable to agriculture, development, and channelization activities in the urbanizing watershed. Our results showed that the urban watershed exported over five times more sediment per unit area compared with the forested watershed. Sediment concentration was related to flow flashiness in the urban watershed and to flow magnitude in the forested watershed. Sediments from the forested watershed were dominated by organic matter, whereas mineral matter dominated sediment from the urban stream. Our results indicated that a significant shift in sediment quality and quantity are likely to occur as forested watersheds are transformed by urbanization in coastal South Carolina.     

Effects of Watershed Impervious Cover on Dissolved Silica Loading in Storm Flow1

Abstract: Dissolved silica (DSi) availability is a factor that affects the composition of algal populations in aquatic ecosystems. DSi cycling is tightly linked to the hydrological cycle, which is affected by human alterations of the landscape. Development activities that increase impervious cover change watershed hydrology and may increase the discharge of DSi‐poor rainwater and decrease the discharge of DSi‐rich ground water into aquatic ecosystems, possibly shifting algal community composition toward less desirable assemblages. In this study, DSi loadings from two adjacent coastal watersheds with different percent impervious cover were compared during four rain and five nonrain events. Loadings in the more impervious watershed contained a significantly larger proportion of surface runoff than base flow (ground‐water discharge) and had lower [DSi] water during rain events than the less impervious watershed. Application of the Soil Conservation Service Curve Number (CN) method showed that the minimum rainfall height necessary to yield runoff was significantly lower for the more impervious watershed, implying that runoff volumes increase with impervious cover as well as the frequency of runoff‐yielding events. Empirical data collected during this study and estimates derived from the CN method suggest that impervious cover may be responsible for both short‐term DSi limitation during rain events as well as long‐term reduction of DSi inputs into aquatic ecosystems.     

Urbanization and Its Effect On Runoff in the Whiteoak Bayou Watershed,Texas1

Abstract: The capacity of a watershed to urbanize without changing its hydrologic response and the relationship between that response and the spatial configuration of the developed areas was studied. The study was conducted in the Whiteoak Bayou watershed (223 km²), located northwest of Houston, Texas, over an analysis period from 1949 to 2000. Annual development data were derived from parcel data collected by the Harris County Appraisal District. Using these data, measures of the spatial configuration of the watershed urban areas were calculated for each year. Based on regression models, it was determined that the annual runoff depths and annual peak flows depended on the annual precipitation depth, the developed area and the maximum 12‐h precipitation depth on the day and day before the peak flow took place. It was found that, since the early 1970s, when the watershed reached a 10% impervious area, annual runoff depths and peak flows have increased by 146% and 159%, respectively. However, urbanization is responsible for only 77% and 32% of the increase, respectively, while precipitation changes are responsible for the remaining 39% and 96%, respectively. Likewise, an analysis of the development data showed that, starting in the early 1970s, urbanization in the watershed consisted more of connecting already developed areas than of creating new ones, which increases the watershed’s conveyance capacity and explains the change in its response. Before generalizing conclusions, though, further research on other urban watersheds with different urbanization models appears to be necessary.     

Spatial and Temporal Patterns of Low Streamflow and Precipitation Changes in the Chesapeake Bay Watershed

Spatial and temporal patterns in low streamflows were investigated for 183 streamgages located in the Chesapeake Bay Watershed for the period 1939–2013. Metrics that represent different aspects of the frequency and magnitude of low streamflows were examined for trends: (1) the annual time series of seven‐day average minimum streamflow, (2) the scaled average deficit at or below the 2% mean daily streamflow value relative to a base period between 1939 and 1970, and (3) the annual number of days below the 2% threshold. Trends in these statistics showed spatial cohesion, with increasing low streamflow volume at streamgages located in the northern uplands of the Chesapeake Bay Watershed and decreasing low streamflow volume at streamgages in the southern part of the watershed. For a small subset of streamgages (12%), conflicting trend patterns were observed between the seven‐day average minimum streamflow and the below‐threshold time series and these appear to be related to upstream diversions or the influence of reservoir‐influenced streamflows in their contributing watersheds. Using multivariate classification techniques, mean annual precipitation and fraction of precipitation falling as snow appear to be broad controls of increasing and decreasing low‐flow trends. Further investigation of seasonal precipitation patterns shows summer rainfall patterns, driven by the Atlantic Multidecadal Oscillation, as the main driver of low streamflows in the Chesapeake Bay Watershed.     

FOREST COVER,IMPERVIOUS‐SURFACE AREA,AND THE MITIGATION OF STORMWATER IMPACTS1

ABSTRACT: For 20 years, King County, Washington, has implemented progressively more demanding structural and nonstructural strategies in an attempt to protect aquatic resources and declining salmon populations from the cumulative effects of urbanization. This history holds lessons for planners, engineers, and resource managers throughout other urbanizing regions. Detention ponds, even with increasingly restrictive designs, have still proven inadequate to prevent channel erosion. Costly structural retrofits of urbanized watersheds can mitigate certain problems, such as flooding or erosion, but cannot restore the predevelopment flow regime or habitat conditions. Widespread conversion of forest to pasture or grass in rural areas, generally unregulated by most jurisdictions, degrades aquatic systems even when watershed imperviousness remains low. Preservation of aquatic resources in developing areas will require integrated mitigation, which must including impervious‐surface limits, forest‐retention policies, stormwater detention, riparian‐buffer maintenance, and protection of wetlands and unstable slopes. New management goals are needed for those watersheds whose existing development precludes significant ecosystem recovery; the same goals cannot be achieved in both developed and undeveloped watersheds.     

INCREASED BASEFLOW IN IOWA OVER THE SECOND HALF OF THE 20TH CENTURY1

ABSTRACT: Historical trends in annual discharge characteristics were evaluated for 11 gauging stations located throughout Iowa. Discharge records from nine eight‐digit hydrologic unit code (HUC‐8) watersheds were examined for the period 1940 to 2000, whereas data for two larger river systems (Cedar and Des Moines Rivers) were examined for a longer period of record (1903 to 2000). In nearly all watersheds evaluated, annual base flow, annual minimum flow, and the annual base flow percentage significantly increased over time. Some rivers also exhibited increasing trends in total annual discharge, whereas only the Maquoketa River had significantly decreased annual maximum flows. Regression of stream discharge versus precipitation indicated that more precipitation is being routed into streams as base flow than as storm flow in the second half of the 20th Century. Reasons for the observed stream flow trends are hypothesized to include improved conservation practices, greater artificial drainage, increasing row crop production, and channel incision. Each of these reasons is consistent with the observed trends, and all are likely responsible to some degree in most watersheds.     

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.     

Effects of Impervious Cover at Multiple Spatial Scales on Coastal Watershed Streams1

Abstract: The spatial scale and location of land whose development has the strongest influence on aquatic ecosystems must be known to support land use decisions that protect water resources in urbanizing watersheds. We explored impacts of urbanization on streams in the West River watershed, New Haven, Connecticut, to identify the spatial scale of watershed imperviousness that was most strongly related to water chemistry, macroinvertebrates, and physical habitat. A multiparameter water quality index was used to characterize regional urban nonpoint source pollution levels. We identified a critical level of 5% impervious cover, above which stream health declined. Conditions declined with increasing imperviousness and leveled off in a constant state of impairment at 10%. Instream variables were most correlated (0.77 ≤ |r| ≤ 0.92, p < 0.0125) to total impervious area (TIA) in the 100‐m buffer of local contributing areas (～5‐km² drainage area immediately upstream of each study site). Water and habitat quality had a relatively consistent strong relationship with TIA across each of the spatial scales of investigation, whereas macroinvertebrate metrics produced noticeably weaker relationships at the larger scales. Our findings illustrate the need for multiscale watershed management of aquatic ecosystems in small streams flowing through the spatial hierarchies that comprise watersheds with forest‐urban land use gradients.     

NATURAL RESTABILIZATION OF STREAM CHANNELS IN URBAN WATERSHEDS1

ABSTRACT: Stream channels are known to change their form as a result of watershed urbanization, but do they restabilize under subsequent conditions of constant urban land use? Streams in seven developed and developing watersheds (drainage areas 5–35 km²) in the Puget Sound lowlands were evaluated for their channel stability and degree of urbanization, using field and historical data. Protocols for determining channel stability by visual assessment, calculated bed mobility at bankfull flows, and resurveyed cross‐sections were compared and yielded nearly identical results. We found that channel restabilization generally does occur within one or two decades of constant watershed land use, but it is not universal. When (or if) an individual stream will restabilize depends on specific hydrologic and geomorphic characteristics of the channel and its watershed; observed stability is not well predicted by simply the magnitude of urban development or the rate of ongoing land‐use change. The tendency for channel restabilization suggests that management efforts focused primarily on maintaining stability, particularly in a still‐urbanizing watershed, may not always be necessary. Yet physical stability alone is not a sufficient condition for a biologically healthy stream, and additional rehabilitation measures will almost certainly be required to restore biological conditions in urban systems.     

Effects of Changing Forest and Impervious Land Covers on Discharge Characteristics of Watersheds

Effects of changing patterns of forest and impervious land covers on hydrologic regimes of watersheds were evaluated for urban and rural areas of the lower Cedar River drainage near Seattle, Washington. Land cover characterizations were used in a spatially explicit hydrology model to assess effects of land covers on watershed hydrology during presettlement conditions (full forest cover), 1991 and 1998. For the presettlement to 1991 period, urban watersheds showed decreases in forest covers (range 63% to 83%) and increases in impervious surfaces (range 43% to 71%). Rural watersheds showed similar patterns but smaller changes, with forest covers decreasing (range 28% to 34%) and impervious surfaces increasing (range 8% to 15%). For the 1991 and 1998 period, changes in forest covers for urban and rural watershed were <24%, with losses in some watersheds and regeneration in others. Impervious surfaces continued to increase, but increases were larger in rural (range 38% to 60%) than in urban watersheds (range 4% to 27%). Flood-frequency curves indicated that discharge rates (m sec^–1) for all watersheds were higher in 1991 and 1998 than historical and suggested that chances for floods increase because of changing land covers. The largest increases in discharge rates were in urban watersheds, with rates for 2-year, 10-year, and 25-year recurrence intervals being more than two times greater than the rate during historical conditions. Changes in flow regimes were indicated by presettlement discharge levels of less frequent recurrence intervals (10-year and 25-year) occurring in posturbanization times (1991 and 1998) during more frequent intervals (2-year and 10-year). Normalized flows (m yr^–1) of watersheds for 2-year, 10-year, and 25-year recurrence intervals indicate how flow regimes in 1991 and 1998 can change as functions of different areas of land covers. During 1991 and 1998, abrupt increases in flows occurred when forest covers were low (range 17% to 37%) and impervious surfaces were >46%. In contrast, the lowest flows occurred when forest covers were most extensive (range 59% to 81%) and impervious surfaces were <23%. We conclude that our use of spatial characterizations of impervious surfaces and forested covers in a spatially explicit hydrology model provides a robust approach for revealing how variations in different types and spatial distributions of land covers can affect flood regimes and flows of different watersheds.     

BASE FLOW RECESSION RATES,LOW FLOWS,AND HYDROLOGIC FEATURES OF SMALL WATERSHEDS IN PENNSYLVANIA,USA1

This paper examines the relationships between measurable watershed hydrologic features, base flow recession rates, and the Q_7,10 low flow statistic (the annual minimum seven‐day average streamflow occurring once every 10 years on average). Base flow recession constants were determined by analyzing hydrograph recession data from 24 small (>130 km²), unregulated watersheds across five major physiographic provinces of Pennsylvania, providing a highly variable dataset. Geomorphic, hydrogeologic, and land use parameters were determined for each watershed. The base flow recession constant was found to be most strongly correlated to drainage density, geologic index, and ruggedness number (watershed slope); however, these three parameters are intercorrelated. Multiple regression models were developed for predicting the recession rate, and it was found that only two parameters, drainage density and hydrologic soil group, were required to obtain good estimates of the recession constant. Equations were also developed to relate the recession rates to Q_7,10 per unit area, and to the Q_7,10/Q₅₀ ratio. Using these equations, estimates of base flow recession rates, Q_7,10, and streamflow reduction under drought conditions can be made for small, ungaged basins across a wide range of physiography.     

STREAM HEALTH AFTER URBANIZATION1

ABSTRACT: Urban development has compromised the quality of physical elements offish habitat in low‐order spawning and rearing streams. In order to identify where priorities should lie in stream rehabilitation, field surveys of a number of streams were conducted near Vancouver, British Columbia. All of the streams were located in watersheds which were urbanized approximately 20 years earlier. The study watersheds ranged from 5 to 77 percent total impervious area (percent TIA). The urban streambeds were found to have less fine material and slightly higher values of intragravel dissolved oxygen than in rural streams. This improved gravel quality is attributed to the higher peak flows generated by impervious areas, and the reduced recruitment of fine material in the urban watersheds. Summer base flow was uniformly low when imperviousness was above 40 percent, evidenced by a decrease in velocity rather than water depth. Large woody debris (LWD) was scarce in all streams with > 20 percent TIA. A healthy buffer zone and abundant LWD were found to stabilize stream banks. The introduction of LWD is considered the most important strategy for stream rehabilitation. Stormwater detention ponds, in contrast, are concluded to have few hydrological benefits if constructed after a stream has reached its urban equilibrium.     

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.     

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.     

Use of Streamflow Data to Estimate Base Flow/Ground‐Water Recharge For Wisconsin1

Abstract: The average annual base flow/recharge was determined for streamflow‐gaging stations throughout Wisconsin by base‐flow separation. A map of the State was prepared that shows the average annual base flow for the period 1970‐99 for watersheds at 118 gaging stations. Trend analysis was performed on 22 of the 118 streamflow‐gaging stations that had long‐term records, unregulated flow, and provided aerial coverage of the State. The analysis found that a statistically significant increasing trend was occurring for watersheds where the primary land use was agriculture. Most gaging stations where the land cover was forest had no significant trend. A method to estimate the average annual base flow at ungaged sites was developed by multiple‐regression analysis using basin characteristics. The equation with the lowest standard error of estimate, 9.5%, has drainage area, soil infiltration and base flow factor as independent variables. To determine the average annual base flow for smaller watersheds, estimates were made at low‐flow partial‐record stations in 3 of the 12 major river basins in Wisconsin. Regression equations were developed for each of the three major river basins using basin characteristics. Drainage area, soil infiltration, basin storage and base‐flow factor were the independent variables in the regression equations with the lowest standard error of estimate. The standard error of estimate ranged from 17% to 52% for the three river basins.     

Estimation of Nonpoint Source Nitrate Concentrations in Indiana Rivers Based on Agricultural Drainage in the Watershed

Subsurface tile‐drained agricultural fields are known to be important contributors to nitrate in surface water in the Midwest, but the effect of these fields on nitrate at the watershed scale is difficult to quantify. Data for 25 watersheds monitored by the Indiana Department of Environmental Management and located near a U.S. Geological Survey stream gage were used to investigate the relationship between flow‐weighted mean concentration (FWMC) of nitrate‐N and the subsurface tile‐drained area (DA) of the watershed. The tile DA was estimated from soil drainage class, land use, and slope. Nitrate loads from point sources were estimated based on reported flows of major permitted facilities with mean nitrate‐N concentrations from published sources. Linear regression models exhibited a statistically significant relationship between annual/monthly nonpoint source (NPS) nitrate‐N and DA percentage. The annual model explained 71% of the variation in FWMC of nitrate‐N. The annual and monthly models were tested in 10 additional watersheds, most with absolute errors within 1 mg/l in the predicted FWMC. These models can be used to estimate NPS nitrate for unmonitored watersheds in similar areas, especially for drained agricultural areas where model performance was strongest, and to predict the nitrate reduction when various tile drainage management techniques are employed.     

Impacts of Urbanization on Stream Habitat and Fish Across Multiple Spatial Scales

We analyzed the relation of the amount and spatial pattern of land cover with stream fish communities, in-stream habitat, and baseflow in 47 small southeastern Wisconsin, USA, watersheds encompassing a gradient of predominantly agricultural to predominantly urban land uses. The amount of connected impervious surface in the watershed was the best measure of urbanization for predicting fish density, species richness, diversity, and index of biotic integrity (IBI) score; bank erosion; and base flow. However, connected imperviousness was not significantly correlated with overall habitat quality for fish. Nonlinear models were developed using quantile regression to predict the maximum possible number of fish species, IBI score, and base flow for a given level of imperviousness. At watershed connected imperviousness levels less than about 8%, all three variables could have high values, whereas at connected imperviousness levels greater than 12% their values were inevitably low. Connected imperviousness levels between 8 and 12% represented a threshold region where minor changes in urbanization could result in major changes in stream condition. In a spatial analysis, connected imperviousness within a 50-m buffer along the stream or within a 1.6-km radius upstream of the sampling site had more influence on stream fish and base flow than did comparable amounts of imperviousness further away. Our results suggest that urban development that minimizes amount of connected impervious surface and establishes undeveloped buffer areas along streams should have less impact than conventional types of development.     

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

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

Nutrient and Sediment Production,Watershed Characteristics,and Land Use in the Tahoe Basin,California‐Nevada1

Abstract: In efforts to control the degradation of water quality in Lake Tahoe, public agencies have monitored surface water discharge and concentrations of nitrogen, phosphorus, and suspended sediment in two separate sampling programs. The first program focuses on 20 watersheds varying in size from 162 to 14,000 ha, with continuous stream gaging and periodic sampling; the second focuses on small urbanized catchments, with automated sampling during runoff events. Using data from both programs, we addressed the questions (1) what are the fluxes and concentrations of nitrogen and phosphorus entering the lake from surface runoff; (2) how do the fluxes and concentrations vary in space and time; and (3) how are they related to land use and watershed characteristics? To answer these questions, we calculated discharge‐weighted average concentrations and annual fluxes and used multiple regression to relate those variable to a suite of GIS‐derived explanatory variables. The final selected regression models explain 47‐62% of the variance in constituent concentrations in the stormwater monitoring catchments, and 45‐72% of the variance in mean annual yields in the larger watersheds. The results emphasize the importance of impervious surface and residential density as factors in water quality degradation, and well‐developed soil as a factor in water quality maintenance.