

Section 1
Section 2
Section 3
Section 4
Section 5
Section 6
Section 7
Section 8
Section 9
Section 10
Appendices
Appendix A
Appendix B
Appendix C
Appendix D
Appendix E
Appendix F
Appendix G
This summary was developed by USEPA Region 5 - Chicago.
This chapter provides an overview of practical methods for estimating short-term and long-term surface water quality impacts related to urban and suburban development sites. These methods may be used by the urban planner or developer to estimate the impacts on receiving waters from development which may result based on various planning conditions and assumptions. Using the approaches presented here, mitigation of water quality impacts may be tested based upon incorporation of various types of management practices.
There are two principal types of water quality impacts typically associated with urban and suburban development. The first includes the impacts related to the construction phase of development as soils which are destabilized due to clearing grading and excavation are subject to increased erosion by wind and water. Eroded soils associated with construction activity can be carried offsite and deposited in receiving waters such as lakes, rivers and wetlands. Adverse impacts related to these sediments include increased turbidity and habitat modification, including smothering of spawning beds. While the construction phase itself may be relatively short-lived, the impacts to receiving waters from poorly managed construction activities may be extremely severe and long-lasting, particularly to sensitive areas such as wetlands and inland lakes.
Once the construction phase is over, other receiving water quality impacts may become more pronounced due to potentially dramatic changes to the area's hydrology (reduced baseflow and exaggerated peak flow volumes), and the change in land use compared to predevelopment conditions. The increase in impervious areas causes a resultant increase in runoff rates and volumes. This can result in increased streambank erosion and associated water quality problems.
The increased runoff also accelerates the transport of land-borne pollutants into receiving waters. Typical pollutants which may be found in urban storm water at significant levels include heavy metals, oil and grease, pesticides, fertilizers and other nutrients, and toxic organic contaminants. Runoff from roadways and parking lots may cause significant elevations in receiving water temperatures during summer months. Winter road deicing activities can contribute high levels of chlorides or sediment.
In order to properly manage and maintain urban water resources, the impacts associated with new development must be carefully evaluated. Post-development impacts may be evaluated in terms of short-term (acute) impacts, and long-term (chronic) impacts. Short term impacts include the changes to a receiving water's chemistry, hydrology, temperature, etc, caused by individual runoff events, and are typically on a timescale of hours to days. Long-term impacts are those which are manifested in the weeks-to-years timescale, and include changes to the dry and wet weather hydrology, streambank morphology, and water chemistry of the receiving water. Long-term chemical impacts are most critical for receiving waters with longer residence times such as lakes and wetlands, and for slower moving stream segments.
In terms of the changes to a receiving water's chemistry due to urban runoff, pollutant concentrations are best used to evaluate short-term effects, while pollutant loadings are appropriate for assessing long-term impacts. Land use planners and developers need to understand these impacts and carefully plan in order to mitigate the negative water quality impacts of development. Part of the analysis should be to evaluate changes in both the annual mass of pollutants exported from a developing area (pollutant loading), and instream pollutant concentration related to runoff from new development or redevelopment.
Loading estimates may focus on nutrients such as phosphorus and nitrogen which contribute to algal blooms in lakes and ponds when the assimilative capacity is exceeded. Estimated loadings can be compared with any existing load allocation limitations for a given receiving water. Even when load allocations do not currently exist, loading estimates are very useful for predicting gross changes in the export of various parameters (sediment, oxygen demanding substances, toxic metals and organics, nutrients), and allow for the analysis of various best management practice (BMP) alternatives to modulate any increased loading of pollutants.
Concentration estimates can be compared with applicable State water quality standards to provide an indication of the likelihood that those standards will be exceeded as a result of storm water discharges. This analysis will help in the planning of BMPs to reduce short term impacts such as acute aquatic toxicity, biochemical oxygen demand and bacteria.
The following is a summary of three methods which may be used to estimate water quality impacts of new development with respect to increased pollutant loading and pollutant concentration.
Pollutant export estimates for a wide variety of pollutants under various planning assumptions can be estimated using the Simple Method (Schueler, 1987). The method is very easy to use as it requires only information which is readily available and does not involve the use of computer models to calculate load estimates. It is recommended that the method be limited in application to sites less than 1 square mile in area.
The annual mass export of a given pollutant in urban runoff may estimated by the following basic form of the Simple Method:
(EQ 1) L = (P)(Pj)(Rv)(C)(A)(0.227) - (where concentration is in mg/l), or
(EQ 1a) L = (P)(Pj)(Rv)(C)(A)(0.000227) - (where concentration is in ug/l)
where:
L = annual mass of pollutant export (lbs/yr)
P = annual precipitation (inches)
Pj = correction factor for smaller storms which do not produce runoff (dimensionless)
Rv = runoff coefficient (dimensionless)
C = average concentration of pollutant
A = site area (acres)
Annual precipitation Where site specific values for P are not available these can be estimated from Figure 1 and Table 1. In Illinois, reasonable estimates fall between 30 inches per year in the northern and central parts of the State to 42 inches per year in the extreme southern section.
FIGURE 1. Rain Zones for the United States (EPA, 1989)

| Rain Zone |
Number of Storms |
COV |
Precip (in) |
COV |
|---|---|---|---|---|
Northeast |
70 |
0.13 |
34.6 |
0.18 |
Northeast- Coastal |
62 |
0.12 |
41.4 |
0.21 |
Mid-Atlantic |
62 |
0.13 |
39.5 |
0.18 |
Central |
68 |
0.14 |
41.9 |
0.19 |
North Central |
55 |
0.16 |
29.8 |
0.22 |
Southeast |
65 |
0.15 |
49.0 |
0.20 |
East Gulf |
68 |
0.17 |
53.7 |
0.23 |
East Texas |
41 |
0.22 |
31.2 |
0.29 |
West Texas |
30 |
0.27 |
17.3 |
0.33 |
Southwest |
20 |
0.30 |
7.4 |
0.37 |
West Inland |
14 |
0.38 |
4.9 |
0.43 |
Pacific South |
19 |
0.36 |
10.2 |
0.42 |
Northwest Inland |
31 |
0.23 |
11.5 |
0.29 |
Pacific Central |
32 |
0.26 |
18.4 |
0.33 |
Pacific Northwest |
71 |
0.15 |
35.7 |
0.19 |
COV = Coefficient of Variation = Standard Deviation/Mean
Correction Factor This factor is used to account for smaller storms which produce no runoff. The value of Pj may be estimated to be 0.9 where more precise data are unavailable.
Runoff Coefficient Rv represents that fraction of precipitation which appears as runoff. This may be estimated from the following:
(EQ 2)Rv = 0.05 + 0.009(I) (Schueler, 1987)
where I is the impervious area for the site expressed as percent. I may be estimated by summing the area of impervious surfaces dividing by the total area.
Alternatively, I may be estimated for residential areas by:
(EQ 3) I = 9(PD)1/2 (Shelly, 1988)
where PD is the population density in persons/acre.
Pollutant Concentration The concentration of pollutant C, may be determined from flow-weighted composite samples representative of annual average values in urban runoff from a given area. Where such data are not available, estimates may be based on data from the NURP database or other reliable sources. A table of C values compiled from NURP data is provided in Table 2. Other data on pollutant concentrations (Schueler, 1987), are presented in Tables 3 and 4.
| Pollutant |
For Media Urban Site |
For 90th Percentile Urban Site |
Coefficient of Variation |
|---|---|---|---|
TSS (mg/l) |
100 |
300 |
1-2 |
BOD (mg/l) |
9 |
15 |
0.5-1.0 |
COD (mg/l) |
65 |
140 |
0.5-1.0 |
Tot. P (mg/l) |
0.33 |
0.70 |
0.5-1.0 |
Sol. P (mg/l) |
0.12 |
0.21 |
0.5-1.0 |
TKN (mg/l) |
1.5 |
3.3 |
0.5-1.0 |
NO2 + 3 -N (mg/l) |
0.68 |
1.75 |
0.5-1.0 |
Copper (ug/l) |
34 |
93 |
0.5-1.0 |
Lead (ug/l) |
144 |
350 |
0.5-1.0 |
Zinc (ug/l) |
160 |
500 |
0.5-1.0 |
| Pollutant |
Residential |
Mixed |
Commercial |
Open/ Nonurban |
||||
|---|---|---|---|---|---|---|---|---|
| Med |
COV |
Med |
COV |
Med |
COV |
Med |
COV |
|
BOD mg/l |
10.0 |
0.41 |
7.8 |
0.52 |
9.3 |
0.31 |
-- |
-- |
COD mg/l |
73 |
0.55 |
65 |
0.58 |
57 |
0.39 |
40 |
0.78 |
TSS mg/l |
101 |
0.96 |
67 |
1.14 |
69 |
0.85 |
70 |
2.92 |
Total P |
383 |
0.69 |
263 |
0.75 |
201 |
0.67 |
121 |
1.66 |
Soluble P |
143 |
0.46 |
56 |
0.75 |
80 |
0.71 |
26 |
2.11 |
TKN |
1900 |
0.73 |
1288 |
0.50 |
1179 |
0.43 |
965 |
1.00 |
Nit.NO2+NO3 |
736 |
0.83 |
558 |
0.67 |
572 |
0.48 |
543 |
0.91 |
Copper ug/l |
144 |
0.75 |
114 |
1.35 |
104 |
0.68 |
30 |
1.52 |
Lead ug/l |
33 |
0.99 |
27 |
1.32 |
29 |
0.81 |
-- |
-- |
Zinc ug/l |
135 |
0.84 |
154 |
0.78 |
226 |
1.07 |
195 |
0.66 |
Source: NURP (EPA 1983)
| Pollutant |
Concentration |
|---|---|
COD (mg/l) |
>40 |
Tot. P (mg/l) |
0.15 |
Sol. P (mg/l) |
0.04 |
TKN (mg/l) |
0.61 |
NO2 -N (mg/l) |
0.17 |
A proposed 25 acre development in Northeastern Illinois would convert a woodland area (I = 2%) to single family homes and townhouses. The total imperviousness would be 40%. Estimate the post-development increases in phosphorus and total Kjeldahl nitrogen (TKN) loadings.
The annual precipitation is assumed to be 30 inches/year (Table 1). The runoff coefficient is calculated from EQ 2.
Prior to development, I = 2%:
Rv = 0.05 + 0.009(2) = 0.068
After development:
Rv = 0.05 + 0.009(40) = 0.41
The concentration values C, are taken from Tables 3 and 4. (Mean NURP concentration values are assumed)
| Parameter | Pre-development | Post-development |
|---|---|---|
| P | 30 inches/year | 30 inches/year |
| Pj | 0.9 | 0.9 |
| Rv | 0.068 | 0.41 |
| C (TKN) | 0.61 mg/l | 1.5 mg/l |
| C (total P) | 0.15 mg/l | 0.33 mg/l |
| A | 25 acres | 25 acres |
Annual loadings are computed from EQ. 1:
Pre-development:
TKN = [(30 in/yr)(0.9)(0.068)(0.61 mg/l)(25 acres)(0.227)] = 6.4 lbs/yr
P-total = [(30 in/yr)(0.9)(0.068)(0.15 mg/l)((25 acres)(0.227)] = 1.6 lbs/yr
Post-development:
TKN = [(30 in/yr)(0.9)(0.41)(1.5 mg/l)(25 acres)(0.227)] = 94.2 lbs/yr
P-total = [(30 in/yr)(0.9)(0.41)(0.33 mg/l)((25 acres)(0.227)] = 20.7 lbs/yr
Results:
| Parameter | Pre-devel. | Post-devel. | %Increase |
|---|---|---|---|
| TKN | 6.4 lbs/yr | 94.2 lbs/yr | 1472% |
| P-total | 1.6 lbs/yr | 20.7 lbs/yr | 1294% |
For the above example, what would be the post development nutrient increase if the total imperviousness were limited to 25%?
From Equation 2, Rv = 0.28. From Equation 1, the post development increase would be:
| Parameter | Pre-devel. | Post-devel. | % Increase |
|---|---|---|---|
| TKN | 6.4 lbs/yr | 64.3 lbs/yr | 1000 |
| P-total | 1.6 lbs/yr | 14.1 lbs/yr | 884 |
Suggest a BMP or group of BMPs which could potentially limit the export of total phosphorus to within 50 % of pre-development levels.
Table 15 provides a summary of pollutant removal efficiencies for
various storm water runoff control practices. Wet ponds or multiple pond
systems are the most reliable practices for controlling nutrients in
runoff, and are also generally effective in removing other pollutants of
concern.
Loading Functions The simple method may be used to convert typical concentration values to estimates of annual mass loadings. Also known as loading functions, these estimates can be based upon unit area for various land use types and per cent site imperviousness, or other similar constants. Loading functions allow for a direct estimation of pollutant loading for various land use types. Table 4 presents calculated loading functions from various land use types.
| Land Use |
I |
Total Phos. |
TKN |
BOD 5-day |
Zinc |
Lead |
|---|---|---|---|---|---|---|
| Rural Residential |
0 5 10 |
0.10 0.19 0.28 |
0.45 0.86 1.27 |
2.70 5.16 7.62 |
0.05 0.09 0.14 |
0.04 0.08 0.12 |
| Large Lot Single Family |
10 15 20 |
0.28 0.38 0.46 |
1.27 1.73 2.09 |
7.62 10.4 12.5 |
0.04 0.18 0.22 |
0.12 0.17 0.20 |
| Medium Density |
20 25 30 35 |
0.46 0.55 0.64 0.74 |
2.09 2.50 2.92 3.36 |
12.5 15.0 17.5 20.2 |
0.22 0.27 0.31 0.36 |
0.20 0.24 0.28 0.32 |
| Town-house |
35 40 45 50 |
0.74 0.83 0.92 1.01 |
3.36 3.77 4.18 4.59 |
20.2 22.6 25.1 27.5 |
0.36 0.40 0.45 0.49 |
0.32 0.36 0.40 0.44 |
| Garden Apartment |
50 55 60 |
1.01 1.10 1.19 |
4.59 5.00 5.41 |
27.5 30.0 32.5 |
0.49 0.53 0.58 |
0.44 0.48 0.52 |
| High Rise, Light Commercial/ Industrial |
60 65 70 75 80 |
1.19 1.28 1.37 1.46 1.55 |
5.41 5.82 6.23 6.64 7.05 |
32.5 34.9 37.4 39.8 42.3 |
0.58 0.62 0.66 0.70 0.75 |
0.52 0.56 0.60 0.64 0.68 |
| Heavy Commercial, Shopping Center |
80 85 90 95 100 |
1.55 1.64 1.72 1.82 1.91 |
7.05 7.45 7.82 8.27 8.68 |
42.3 44.7 46.9 49.6 52.1 |
0.75 0.80 0.83 0.88 0.93 |
0.68 0.72 0.75 0.79 0.83 |
1 P=30 inches, Pj=0.9, Rv=0.05+0.009(I), A= 1 acre,
C= mean NURP values from Table 2
Similar loading functions can be made using the data in Table 3 or other appropriate data.
Total annual loads are estimated by multiplying the area associated with each given land use type by the loading function for that land use:
(EQ. 4) L = SUM (Lx)(Ax)
where:
L = total loading (lbs/yr)
Lx = loading function for land use x (lbs/acre/yr)
Ax = area of land use x (acre)
In another variation of the Simple Method, Heaney (Mills et al, 1985) has developed a loading function based on population density and street cleaning frequency:
(EQ.5) Lx = (ax)(Fx)(Yx)(P)
where:
Lx = loading function for land use x (lbs/acre)
ax = pollutant concentration factor (lbs/acre/in)
Fx = population density function
Yx = street cleaning factor
P = annual precipitation (inches)
Total loading is calculated using EQ. 5. Typical ax values are given in Table 5.
| Land Use |
BOD |
TSS |
PO4 |
Nit. |
|---|---|---|---|---|
Residential |
0.78 |
16 |
0.033 |
0.13 |
Commercial |
3.13 |
22 |
0.073 |
0.29 |
Industrial |
1.18 |
29 |
0.069 |
0.27 |
Other Developed |
0.11 |
2.7 |
0.009 |
0.06 |
The population density function, Fx is a dimensionless parameter. Typical empirical values for Fx adapted from Heaney et al. are:
1.0 for commercial and industrial development, and
(EQ. 6) 0.142 + 0.134 [0.405(PD)]0.54 for residential
where:
PD = population density (persons/acre)
The street cleaning factor is based upon the street sweeping interval in days (Ns):
Yx = Ns/20 for Ns < 20 days
Yx = 1.0 for Ns > 20 days
EXAMPLE 4
Referring to Example 1, assume that the population density is 25 persons/acre and that street sweeping is performed 1/month, what is the annual nitrogen and phosphorus loading from the development as predicted from EQ. 4 ? How does this compare with the prediction from Example 1?
ANSWER
From TABLE 4, ax is 0.033 (lbs/acre/in).
From EQ. 6, the population density function, Fx = 0.87
Yx is set to 1.0, since street sweeping frequency is less than 1/20 days.
P = 30 in/yr.
From EQ 5, the phosphorus loading function is:Lx = (0.033 lbs/acre-inch)(0.87)(1.0)(30 in/yr) = 0.86 lbs/acre-yr
The total annual load from EQ 3 is:
L = (25 acres)(0.86 lbs/acre-yr) = 22 lbs/yr (Using the simple method in Example 1 predicted 20.7 lbs/yr)
2. PHOSPHORUS LOAD ALLOCATIONS - THE MAINE DEP METHOD
The Maine Department of Environmental Protection has developed a detailed application of the loading functions method for determining changes in phosphorus loadings which may be expected as a result of different urban and suburban development scenarios (Dennis et al. 1989). Estimated phosphorus loadings can be compared with specified phosphorus loading allocations for Maine lakes. In addition, the procedure allows for the estimation of phosphorus loading mitigation, based on the use of various combinations of BMPS. By use of simple desk-top calculations, planners and developers are able to estimate in advance whether proposed development areas will comply with the State's phosphorus loading allocations.
The acceptable increase in phosphorus export is determined by:
(EQ. 7) Lp = (FC)/D
where:
Lp = acceptable increase in the phosphorus loading function (lbs/acre/yr)
1F = phosphorus coefficient for the lake watershed (lbs/ppb/yr)
C = acceptable increase in lake phosphorus concentration (ppb)
D = future area to be developed over next 50 years in the watershed (acres)
1 F factors have been determined for specific lakes in Maine. Similar targets may be established for waterbodies in other areas. In the absence of specific loading limitations, the process may be used to estimate the increase in phosphorus loading resulting from a proposed development.
C, the acceptable increase in phosphorus concentration, is a function of
existing water quality and the level of desired protection. C values are
given in Table 6.
| Water Quality Category | Lake Protection Level |
||
|---|---|---|---|
| High | Medium | Low | |
Outstanding |
0.5 |
1.0 |
1.0 |
Good |
1.0 |
1.5 |
2.0 |
Moderate/Stable |
1.0 |
1.25 |
1.5 |
Moderate/Sensitive |
0.75 |
1.0 |
1.25 |
Poor/Restorable |
0.1 |
0.5 |
NA |
Poor/Low Priority |
2.0 |
4.0 |
6.0 |
D is determined as the total area minus already developed and undevelopable land (steep slopes, wetlands, parks, etc.) and multiplying by a development factor which estimates the portion of undeveloped land which is likely to be developed. In Maine these development factors are:
0.20 - 0.35 for lake areas near growth centers
0.15 - 0.25 for lake areas subject to seasonal development
0.10 - 0.20 for lakes for which development is shoreline dependent
0.10 - 0.15 for lake areas not subject to development pressure
It is recommended that conservative upper estimates be used for development factors.
The permitted phosphorus export (PPE) for a site is simply:
(EQ. 8) PPE = (Lp)(A)
where:
PPE = permitted phosphorus export for proposed development (lbs/acre)
A = the proposed area of the site (acres)
The proposed area of the development (A) should include all areas except those which are undevelopable such as wetlands > 1 acre, and steep slopes.
The total predicted phosphorus export (TE) for a development site is the summation of export values from roadways, individual houselots, multi-unit housing, commercial and industrial development. Credit is given for phosphorus control measures which are employed. The general equation for Phosphorus export is:
(EQ. 9) (TE) = summation (RE) + (HE) + (CE)
where:
(TE) = total predicted phosphorus export
(RE) = phosphorus export from roadways
(HE) = phosphorus export from individual house lots
(CE) = phosphorus export from multi-unit housing, commercial, and industrial development,
Phosphorus Export from Residential Area Roadways (RE)
Road surface phosphorus export is determined as follows:
(EQ. 9) (RE) = [(FT)(LBS)(TFb)(TFwp)(TFi)(TFo)]/100
where:
(FT) = length of roadway being evaluated (feet)
(LBS) = annual export of phosphorus from 100 feet of roadway, before treatment
(TFb) = treatment factor for buffer strips
(TFwp) = treatment factor for wet ponds
(TFi) = treatment factor for infiltration practice
(TFo) = treatment factor for other treatment factor
The annual export per 100 feet of roadway is calculated as:
(EQ. 10) (LBS) = (road surface width)(0.012) + (road ditch width)(0.004)
Treatment factors (TF) for all the above calculations and those that follow must be numbers between 0 and 1.0 which reflect the long term phosphorus removal efficiency of the treatment practice or practices employed. Tables in Appendix C-1 present some recommended values. Note that lower numbers reflect higher removal efficiencies. It is also evident that the calculation gives greater credit where redundant treatment practices are employed.
Phosphorus export from individual house lots (HE)
The annual phosphorus export from an individual houselot is calculated as:
(EQ. 11) (HE) = (BP)(TFb)(TFwp)(TFi)(TFo)
where:
(BP) = phosphorus export before treatment
(TFb) = treatment factor for buffer strips
(TFwp) = treatment factor for wet ponds
(TFi) = treatment factor for infiltration practice
(TFo) = treatment factor for other treatment factor
| Hydrologic Group |
Area Cleared per Lot |
||
|---|---|---|---|
| <10,000 ft |
>10,000 ft |
>15,000 ft |
|
A |
.27 (.2) |
.30 (.30) |
.35 (.35) |
B |
.32 (.40) |
.39 (.46) |
.49 (.54) |
C |
.34 (.48) |
.44 (.56) |
.58 (.67) |
D |
.36 (.62) |
.47 (.62) |
.62 (.74) |
Note: Values in parentheses are appropriate for sites where more than 40% of timber volume has been harvested within the last 5 years.
Phosphorus Export from Multi-unit Housing, Commercial, and Industrial Development (CE)
Phosphorus Export from multi-unit housing, commercial, and industrial development is calculated as:
(EQ. 12) (CEx) = (Lx)((BLx)(TFb)(TFwp)(TFi)(TFo)
where:
(Lx) = altered land surface area (acres)
(BLx) = additional phosphorus export per acre of altered land surface (lbs/acre)
(TFb) = treatment factor for buffer strips
(TFwp) = treatment factor for wet ponds
(TFi) = treatment factor for infiltration practice
(TFo) = treatment factor for other treatment factor
Values for additional phosphorus export associated with altered land uses are
found in Table 12.
| Land Use Category |
Phosphorus Export Before Treatment |
|---|---|
Lawn A |
.30 lbs/acre |
Lawn B |
.65 lbs/acre |
Lawn C |
.97 lbs/acre |
Lawn D |
1.1 lbs/acre |
Road Ditch |
1.8 lbs/acre |
Road Surface |
5.3 lbs/acre |
Impervious Surfaces |
3.5 lbs/acre |
The total loading from multi-unit housing, commercial and industrial areas is the summation of all areas for various land use categories.
Referring to the proposed development in Example 1, consists of 40 single family units with an average lot size of 0.31 acres, and 29 townhouses with an average of 0.43 acres. The resident soil is Type C. All the runoff will be treated by a wet detention basin. Runoff from the townhouses and roads will also be treated by a 150 foot buffer strip with a slope of 12%. The wet detention basin will be designed with a length to width ratio of 3:1, and mean depth of 5 ft.
The total road length to be added as part of the subdivision is 1,100 feet. The road width is 38 feet with 4 foot shoulders. The main access road which is included in the total road length is 700 feet and has 5 foot ditches on each roadside.
Calculate the additional phosphorus export associated with the proposed development. Also calculate without treatment and compare with Examples 1 and 4.
Referring to Appendix C-1, assume a treatment factor of 0.7 for buffer strips 0.5 for the wet detention pond. The analysis must consider the runoff from the houses and townhouses separately. Assuming that < 10,000 ft2 cleared, phosphorus loading from the single-family dwellings is:
(HE) = (40 lots)(0.48lbs/lot/yr)(0.5) = 9.6 lbs/yr
From the townhouses:
(CE) = (29 lots)(0.43 acres/lot)(0.97 lbs/acre/yr)(0.5)(0.7) = 4.2 lbs/yr
For the roadways, use EQ. 9:
(LBS) = [(38 + 4 + 4)(0.012)] + [(5 + 5)(0.004)] = 0.59 lbs/100 ft
(RE) = (11)(0.59)(0.7)(0.5) = 2.3 lbs/yr
The total additional phosphorus export is then:
HE + CE + RE = 16.1 lbs/yr
Without treatment:
(HE) = (40 lots)(0.48lbs/lot/yr) = 19.2 lbs/yr
From the townhouses:
(CE) = (29 lots)(0.43 acres/lot)(0.97 lbs/acre/yr) = 12.1 lbs/yr
For the roadways, use EQ. 9:
(LBS) = [(38 + 4 + 4)(0.012)] + [(5 + 5)(0.004)] = 0.59 lbs/100 ft
(RE) = (11)(0.59) = 6.5 lbs/yr
The total additional phosphorus export is then:
HE + CE + RE = 37.8 lbs/yr
This is somewhat higher than the solutions to Examples 1 and 4 (20.7 lbs/year and 22 lbs/year, respectively).
A 7.4 acre office complex in an area with type B soils is proposed. The site includes 4.9 acres of lawn area. Rooftop accounts for 0.9 acres. and 1.1 acre for parking, and 0.3 acres for road surface, and 0.2 acres for road ditch. Calculate the additional phosphorus export.
All flows are to be treated by a 100 ft buffer strip with 10% slope and a wet detention pond with a 4:1 length to width ratio and a mean depth of 4 feet.
From Appendix C - 1, TABLES C - 1.1 and C - 1.2, the Treatment factors are 0.6 for the buffer strip and 0.48 for the wet pond. Phosphorus loadings from the various areas are from Table 12 and EQ 9. The total loading is:
(CE) = [(4.9)(0.65) + (0.9)(3.5) + (1.1)((3.5) + (0.3)(5.3) + (0.2)(3.5)] x [(0.6)(0.48)] = 2.77 lbs phosphorus/yr
While the Maine procedure was conceived for use in loading estimates of a particular pollutant (phosphorus), and is specific to the State of Maine, the basic concept can be expanded for use with other pollutants of concern in any type of receiving water anywhere in the country. In order to adapt this procedure, the following types of information are necessary:
Where available, these may be compared with target loading ceilings for pollutants of concern.
Table 13 summarizes data on concentrations of various pollutants in runoff from to urban catchments in Wisconsin (Bannerman et al, 1992). In Table 14 these are converted to annual pollutant load per acre, based on 30 inches of precipitation annually. These can be converted to lbs per square foot by dividing by 43,560.
| Source1 |
Mean Pollutant Concentration2 |
||||||
|---|---|---|---|---|---|---|---|
| TSS mg/l | Total Phos mg/l |
Cd ug/l | Cr ug/l | Cu ug/l | Pb ug/l | Zn ug/l | |
IndustRoof 1 |
54 |
.13 |
.3 |
-- |
7 |
8 |
1348 |
Arterial ST 1 |
875 |
1.01 |
2.8 |
26 |
85 |
85 |
629 |
Arterial ST 2 |
241 |
.53 |
2.6 |
18 |
50 |
55 |
554 |
Feeder ST 1 |
969 |
1.57 |
3.7 |
17 |
97 |
107 |
574 |
Feeder ST 2 |
1085 |
1.77 |
.8 |
7 |
25 |
38 |
245 |
Parking Lot 1 |
475 |
.48 |
1.2 |
16 |
47 |
62 |
361 |
Parking Lot 2 |
91 |
.26 |
.08 |
7 |
21 |
30 |
249 |
Outfall 1 |
174 |
.38 |
1.1 |
7 |
31 |
26 |
295 |
Outfall 2 |
374 |
.86 |
.6 |
5 |
20 |
40 |
254 |
ResiDriveway 2 |
193 |
1.5 |
.5 |
2 |
20 |
20 |
113 |
FlatRoof 2 |
19 |
.24 |
.4 |
-- |
10 |
10 |
363 |
Collector ST 2 |
386 |
1.22 |
1.7 |
13 |
61 |
62 |
357 |
ResiLawn 2 |
457 |
3.47 |
-- |
-- |
13 |
-- |
60 |
ResiRoof 2 |
36 |
.19 |
.2 |
-- |
5 |
10 |
153 |
Source: Bannerman et al, 1992
1Study area 1 described as mainly industrial; Study area 2 described as medium density residential
2Area = 1 acre P = 30 inches/yr
| Source1 |
Pollutant Loading (lbs/acre/year)2 |
||||||
|---|---|---|---|---|---|---|---|
| TSS | Total Phos |
Cd | Cr | Cu | Pb | Zn | |
IndustRoof 1 |
367 |
.88 |
.002 |
-- |
.05 |
.05 |
9.2 |
Arterial ST 1 |
5950 |
6.9 |
0.2 |
.18 |
.58 |
.58 |
4.3 |
Arterial ST 2 |
1639 |
3.6 |
.018 |
.12 |
.34 |
.37 |
3.8 |
Feeder ST 1 |
6589 |
10.7 |
.025 |
.12 |
.66 |
.73 |
3.9 |
Feeder ST 2 |
7378 |
12.0 |
.005 |
.05 |
.17 |
.26 |
1.67 |
Parking Lot 1 |
3230 |
3.26 |
.008 |
.11 |
.32 |
.42 |
2.45 |
Parking Lot 2 |
619 |
1.77 |
.0005 |
.05 |
.14 |
.20 |
1.69 |
Outfall 1 |
1183 |
2.58 |
.007 |
.05 |
.21 |
.18 |
2.01 |
Outfall 2 |
2543 |
5.84 |
.004 |
.03 |
.136 |
.27 |
1.73 |
ResiDriveway 2 |
1312 |
10.2 |
.003 |
.014 |
.136 |
.136 |
.768 |
FlatRoof 2 |
129 |
1.63 |
.003 |
-- |
.068 |
.068 |
2.47 |
Collector ST 2 |
2625 |
8.30 |
.012 |
.088 |
.415 |
.422 |
2.43 |
ResiLawn 2 |
3108 |
23.6 |
-- |
-- |
.088 |
-- |
0 |
ResiRoof 2 |
245 |
1.29 |
.0014 |
-- |
.034 |
.068 |
1.04 |
Source: Bannerman et al, 1992
1Study area 1 described as mainly industrial; Study area 2 described as medium density residential
2Area = 1 acre P = 30 inches/year
Schueler has provided an assessment of effectiveness of various control practices in removing pollutants (Schueler, 1992). A summary is provided in Table 15.
| Storm Water Management |
Pollutant Removal Efficiency |
|||
|---|---|---|---|---|
| TSS |
Nutrients |
Organic |
Metals |
|
Extended Detention |
30-70% |
low to neg. for soluble nutrients |
15-40% for COD |
|
Wet Ponds |
50-90% |
30-90% for total P 40-80% for sol. nutr. |
mod-high removal |
mod-high removal |
Stormwater Wetlands |
slightly higher than wet ponds |
somewhat lower than wet ponds |
||
Multiple Pond System |
Varies with design, but typically enhanced |
|||
Infiltration Trenches |
+90% |
60% |
+90% |
|
Infiltration Basins |
||||
Porous Pavement |
up to 80% |
up to 60% for Phos up to 80% for Nit |
High |
High |
Sand Filters |
85% |
40% for dissolved Phosphorus 35% for Nit |
50-70% |
|
Peat Sand Filters |
85% |
70% for Phosphorus 50% for Nit |
90% for BOD |
|
Grassed Swales |
up to 70% |
30% for Phos 25% for Nit |
50-90% |
|
Filter Strip |
28% |
|||
Source: Schueler et al, 1992
EXAMPLE 7
Returning to Example 1 and Example 5. Assume that the single family units have roofs which are 2,400 square feet (0.055 acre), and the townhouse roofs are 3,000 square feet (0.069 acre). Residential driveways are assumed to average 610 square feet (0.014 acres). Parking for the townhouses consists of 11 lots at an average of 6,500 square feet (0.15 acres) each.
Calculate the zinc loading for the 25 acre single family/townhouse development, using the Tables 14 and 15, and alternatively, EQ 1, before treatment. How do the results compare?
Discussion
The total roof area is:
(40)(0.055) + (29)(0.069) = 4.2 acres rooftops
Total driveway area is:
(40)(0.014 acres) = 0.56 acres driveway
Total parking lot area is:
(11)(0.15 acres) = 1.65 acres parking lot.
The area associated with the roadways (feeder street) is:
(38 + 4 + 4)(1,100) = 50,600 square feet or 1.2 acres
The remaining area is considered to be residential lawn:
25 - [4.2 + 0.56 + 1.65 + 1.2] = 17.4 acres residential lawn
Loading estimates:
Referring to Table 14, the estimated loading rate from rooftops is assumed to be 1.04 lbs zinc/yr, and 0.68 lbs lead/yr. Total loading for the development is:
Zinc - (4.2 acres)(1.04 lbs/acre/yr) = 4.24 lbs/yr
Lead - (4.2 acres)(0.68 lbs/acre/yr) = 2.86 lbs/yr
The estimated loading rate from residential driveways is assumed to be 0.768 lbs zinc/acre/yr, and 0.136 lbs lead/acre/yr. Total loading from residential driveways is:
Zinc - (0.56 acres)(0.768 lbs/acre/yr = 0.43 lbs/yr
Lead - (0.56 acres)(0.136 lbs/acre/yr = 0.076 lbs/yr
The estimated loading rate from parking lots is 1.69 lbs zinc/acre/yr and 0.20 lbs lead/acre/yr (assume study area 2 - medium density residential). Total loading from parking lots is:
Zinc - (1.65 acres)(1.69 lbs/acre/yr) = 2.60 lbs/yr
Lead - (1.65 acres)(0.20 lbs/acre/yr) = 0.33 lbs/yr
The estimated loading rate from the feeder streets is 1.67 lbs zinc/acre/yr, and 0.26 lbs lead/acre/yr. Total loading is:
Zinc - (1.2 acres)(1.67 lbs/acre/yr) = 2.00 lbs/yr
Lead - (1.2 acres)(0.26 lbs/acre/yr) = 0.31 lbs/yr
The estimated loading rate from residential lawns is 0 for zinc and unknown for lead. Therefore, assume no significant increase in metals loading from these areas.
Summing the above the total loading from all areas is gives:
Zinc - 4.24 + 0.43 + 2.60 + 2.00 = 9.3 lbs/yr
Lead - 2.86 + 0.076 + 0.33 + 0.31 = 3.6 lbs/yr
Using EQ 1a, and referring back to Table 2 and Example 1:
Zinc = (30 in/yr)(0.9)(0.068)(0.160 mg/l)(25 acres) = 7.3 lbs/yr
Lead = (30 in/yr)(0.9)(0.068)(0.144 mg/l)(25 acres) = 6.6 lbs/yr
The methods provide reasonable agreement for the proposed development.
The approaches outlined above offer tools for predicting changes in long-term loading rates of pollutants to surface waters as an aid to planning activities. These methods do not provide for estimating short term impacts of urban runoff. Such impacts are more properly viewed as the result of instream pollutant concentrations rather than average loading rates. Predicted instream concentrations can be compared with state water quality standards as a means of predicting water quality standards violations due to urban runoff.
Typically much more complex computer models are employed to predict short term wet weather impacts to receiving waters brought about by urbanization. These models integrate hydrological and instream chemical processes in order to estimate instream pollutant concentrations. Models such as STORM and SWMM require significant data input and site specific verification.
Analysis of data collected as part of the National Urban Runoff Program (NURP) indicates that event mean pollutant concentrations may adequately be specified as a lognormal distribution (EPA, 1986). Because of this, the expected concentration for a given probability for a given pollutant in urban runoff can be determined for a particular data set if the central tendency (median or mean value) and the variability (coefficient of variation or standard deviation) are known. This concentration can be compared to some reference concentration such as a water quality standard to indicate the likelihood that an acute water quality impact will occur in the receiving water. Alternatively, the probability that a given concentration level (such as a water quality standard) will be exceeded can be estimated.
The expected runoff concentration for pollutant x is:
(EQ. 13) Cx = Cm (exp [Z (ln{1+COV}2)1/2]
where:
* Cx = expected concentration of pollutant x
* Z = standard normal probability (for specified probability of occurrence)
* Cm = median pollutant concentration
* COV = coefficient of variation
* For log-transormed data
The probability that a specified concentration will be exceeded can be determined by substituting the concentration level of interest for Cx in EQ. 13, solving the equation for Z, and locating the associated probability for the calculated Z value:
(EQ. 14) Z = (ln[Cx/Cm])/[(ln(1+COV2)1/2]
Z values for various probabilities of occurrence is presented in Table 16.
Median event mean concentrations and coefficients of variation for NURP data for all land use types are presented in Table 3. If sufficient local data are available these may also be used provided they are transformed into logrhythmic form. Illinois Water quality standards for various pollutants is presented in Table 17. For certain metals these are based on hardness.
| Z Value |
Probability of Exceedance |
|---|---|
3.090 |
0.1% |
2.326 |
1% |
2.054 |
2% |
1.881 |
3% |
1.751 |
4% |
1.645 |
5% |
1.476 |
7% |
1.282 |
10% |
1.036 |
15% |
0.842 |
20% |
0.674 |
25% |
0.524 |
30% |
0.385 |
35% |
0.253 |
40% |
0.000 |
50% |
-0.253 |
60% |
-0.524 |
70% |
-0.842 |
80% |
-1.282 |
90% |
| Pollutant |
Acute Standard |
Chronic Standard |
|---|---|---|
Arsenic (ug/l) |
360 |
|
Cadmium (ug/l) |
exp[A + B ln(H)] but not > 50ug/l A=-2.98, B=1.128 |
exp[A + B ln(H)] A=-3.49, B=0.785 |
Hexavalent Chromium (ug/l) |
16 |
11 |
Trivalent Chromium (ug/l) |
exp[A + B ln(H)] A=3.688, B=0.819 |
exp[A + B ln(H)] A=1.561, B=0.819 |
Copper (ug/l) |
exp[A + B ln(H)] A=-1.46, B=0.942 |
exp[A + B ln(H)] A=-1.47, B=0.855 |
Cyanide (ug/l) |
22 |
5.2 |
Lead (ug/l) |
exp[A + B ln(H)] but not >100ug/l A=-1.46, B=1.273 |
NA |
Mercury (ug/l) |
0.5 |
NA |
Barium (mg/l) |
5.0 |
NA |
Boron (mg/l) |
1.0 |
NA |
Chloride (mg/l) |
500 |
NA |
Fluoride (mg/l) |
1.4 |
NA |
Iron Mg/l) |
1.0 |
NA |
Manganese (mg/l) |
1.0 |
NA |
Nickel (mg/l) |
1.0 |
NA |
Phenols (mg/l) |
0.1 |
NA |
Selenium (mg/l) |
1.0 |
NA |
Silver (ug/l) |
5.0 |
NA |
Sulfate (mg/l) |
500 |
NA |
Total Dissolved Solids (mg/l) |
1000 |
NA |
Zinc (mg/l) |
1.0 |
NA |
H = Hardness
EXAMPLE 8
For the development described in Example 1, what is the probability that water quality standards will be violated for zinc and lead, assuming no treatment.
Answer
From Table 15, the acute water quality standard for lead is:
WQSacute = exp[-1.46 + 1.273(ln(H))]
Assuming the Hardness = 100 mg/l, then:
WQSacute = 82 ug/l lead
From Table 15, the acute water quality standard for zinc is:
WQSacute = 1.0 mg/l zinc
From Table 3, the median concentrations and Coefficients of variation for lead and zinc are
Lead: Cm = 33 ug/l, COV = 0.99
Zinc: Cm = 135 ug/l, COV = 0.84
Applying EQ 14, the probability that lead and zinc water quality standards would be exceeded for any given storm (assuming no treatment or dilution) would be estimated to be:
Lead: Z = (ln[82/33])/[(ln(1+ 0.99)2)]1/2
Z = 0.77, which corresponds to an exceedance probability of 20 - 25%
Zinc: Z = (ln[1000/135])/[(ln(1+ 0.84)2)]1/2
Z = 1.81, which corresponds to an exceedance probability of less than 5%
Bannerman, R.T., Dodds, R., Owens, D., Hughes, P. 1992. Sources of Pollutants in Wisconsin Stormwater. Wisconsin Department of Natural Resources Grant Report.
Dennis, J., Noel, J., Miller, D., Eliot, C. 1989. Phosphorus Control in Lake Watersheds: A Technical Guide to Evaluating New Development. Maine Department of Environmental Protection.
Mills, W.B., Porcella, D.B., Ungs, M.J., Gherini, S.A., Summers, K.V., Mok, L, Rupp, G.L., Bowie, G.L., Haith, D.A. Water Quality Assessment: A Screening Procedure for Toxic Pollutants. EPA - 600/6-85-002A
Scheuler, T.R. 1987. Controlling Urban Runoff: A Practical Manual for Planning and Designing Urban BMPs. Metropolitan Washington Council of Governments.
Scheuler, T.R., Kumble, P.A., Heraty, M.A. 1992. A Current Assessment of Urban Best Management Practices. Techniques for Reducing Non-Point Source Pollution in the Coastal Zones. Metropoltan Washington Council of Governments.
Shelly, P. Technical Memorandum to S.A.I.C, 11/03/88.
U.S. Environmental Protection Agency (EPA). 1983. Results of the Nationwide Urban Runoff Program. Volume I. Final Report. Water Planning Division. Washington, D.C.
U.S. Environmental Protection Agency (EPA). Urban Targeting and BMP Selection. Terrene Institute, Washington, D.C.
| Treatment Factor- Wooded (Non-Wooded) |
|||||
|---|---|---|---|---|---|
| Slope |
25 ft |
50 ft |
100 ft |
150 ft |
200 ft |
0-10% |
.75 (.95) |
.4 (.6) |
.2 (.4) |
.1 (.3) |
0 (.2) |
11-15% |
.8 (1.0) |
5. (.7) |
.25 (.45) |
.1 (.3) |
0 (.2) |
16-20% |
.8 (1.0) |
.7 (.9) |
.5 (.7) |
.25 (.45) |
.1 (.3) |
21-30% |
.8 (1.0) |
.75 (.95) |
.7 (.9) |
.6 (.8) |
.3 (.6) |
| Treatment Factor- Wooded (Non-Wooded) |
|||||
|---|---|---|---|---|---|
| Slope |
25 ft |
50 ft |
100 ft |
150 ft |
200 ft |
0-10% |
.75 (.95) |
.6 (.8) |
.4 (.6) |
.2 (.4) |
.1 (.2) |
11-15% |
.8 (1.0) |
.75 (.95) |
.5 (.7) |
.2 (.4) |
.1 (.2) |
16-20% |
.8 (1.0) |
.8 (.1.0) |
.65 (.85) |
.4 (.6) |
.2 (.4) |
21-30% |
.8 (1.0) |
.8 (.1.0) |
.7 (.9) |
.5 (.7) |
.3 (.6) |
| Treatment Factor- Wooded (Non-Wooded) |
|||||
|---|---|---|---|---|---|
| Slope |
25 ft |
50 ft |
100 ft |
150 ft |
200 ft |
0-10% |
.8 (1.0) |
.7 (.9) |
.55 (.75) |
.45 (.65) |
.35 (.55) |
11-15% |
.8 (1.0) |
.75 (.95) |
.6 (.8) |
.5 (.7) |
.4 (.65) |
16-20% |
.8 (1.0) |
.8 (1.0) |
.7 (.9) |
.6 (.8) |
.5 (.65) |
21-30% |
.8 (1.0) |
.8 (1.0) |
.75 (.95) |
.65 (.85) |
.5 (.75) |
| Treatment Factor- Wooded (Non-Wooded) |
|||||
|---|---|---|---|---|---|
| Slope |
25 ft |
50 ft |
100 ft |
150 ft |
200 ft |
0-10% |
.9 (1.0) |
.8 (.65) |
.75 (.8) |
.7 (.8) |
.6 (.75) |
11-15% |
.9 (1.0) |
.85 (1.0) |
.8 (.9) |
.75 (.9) |
.65 (.8) |
16-20% |
.9 (1.0) |
.9 (1.0) |
.85 (1.0) |
.8 (1.0) |
.7 (.85) |
21-30% |
.9 (1.0) |
.9 (1.0) |
.9 (1.0) |
.8 (1.0) |
.75 (.9) |
Source: Maine DEP, 1989
| Number of Storm Volumes |
|||||
|---|---|---|---|---|---|
| 1/2 V |
1V |
2V |
3V |
||
| MEAN DEPTH |
3 ft |
.50 |
.4 |
.33 |
.31 |
5 ft |
.47 |
.34 |
.27 |
.24 |
|
7 ft |
.44 |
.32 |
.24 |
.20 |
|
| Number of Storm Volumes |
|||||
|---|---|---|---|---|---|
| 1/2 V |
1V |
2V |
3V |
||
| MEAN DEPTH |
3 ft |
.56 |
.48 |
.42 |
.40 |
5 ft |
.52 |
.43 |
.36 |
.33 |
|
7 ft |
.51 |
.41 |
.33 |
.29 |
|
| Number of Storm Volumes |
|||||
|---|---|---|---|---|---|
| 1/2 V |
1V |
2V |
3V |
||
| MEAN DEPTH |
3 ft |
.61 |
.55 |
.51 |
.49 |
5 ft |
.59 |
.51 |
.45 |
.43 |
|
7 ft |
.58 |
.49 |
.42 |
.39 |
|
| Number of Storm Volumes |
|||||
|---|---|---|---|---|---|
| 1/2 V |
1V |
2V |
3V |
||
| MEAN DEPTH |
3 ft |
.46 |
.34 |
.26 |
.23 |
5 ft |
.43 |
.31 |
.22 |
.18 |
|
7 ft |
.43 |
.30 |
.19 |
.16 |
|
| Number of Storm Volumes |
|||||
|---|---|---|---|---|---|
| 1/2 V |
1V |
2V |
3V |
||
| MEAN DEPTH |
3 ft |
.50 |
.39 |
.31 |
.28 |
5 ft |
.48 |
.35 |
.26 |
.23 |
|
7 ft |
.47 |
.34 |
.24 |
.22 |
|
| Number of Storm Volumes |
|||||
|---|---|---|---|---|---|
| 1/2 V |
1V |
2V |
3V |
||
| MEAN DEPTH |
3 ft |
.53 |
.43 |
.36 |
.34 |
5 ft |
.52 |
.41 |
.33 |
.30 |
|
7 ft |
.51 |
.40 |
.31 |
.26 |
|
| Number of Storm Volumes |
|||||
|---|---|---|---|---|---|
| 1/2 V |
1V |
2V |
3V |
||
| MEAN DEPTH |
3 ft |
.44 |
.33 |
.23 |
.19 |
5 ft |
.43 |
.30 |
.19 |
.16 |
|
7 ft |
.42 |
.27 |
.18 |
.15 |
|
| Number of Storm Volumes |
|||||
|---|---|---|---|---|---|
| 1/2 V |
1V |
2V |
3V |
||
| MEAN DEPTH |
3 ft |
.47 |
.35 |
.26 |
.23 |
5 ft |
.46 |
.33 |
.23 |
.19 |
|
7 ft |
.46 |
.32 |
.22 |
.18 |
|
| Number of Storm Volumes |
|||||
|---|---|---|---|---|---|
| 1/2 V |
1V |
2V |
3V |
||
| MEAN DEPTH |
3 ft |
.50 |
.39 |
.31 |
.27 |
5 ft |
.49 |
.36 |
.27 |
.24 |
|
7 ft |
.48 |
.35 |
.26 |
.22 |
|