“I“ w-' ' .‘ ‘ ‘ >00: . STREAM WATER QUALITY AS RELATED TO UBIB‘ANEZANON OF ”9 WATERSHED This}: he fire Dogma of M. 5. MICHEGAN STATE UNWEBBITY Alvin Lee Jensen 1966 Illllllml 31 WW Lu! 1m in fl 1111 “W l ‘ LIBRARY Michigan State University firal‘rflq EJHF; _&j'.‘ll\!] Q i 1““ ’ I ‘iiv thi Uu:;_. ultra-I W“; ZZ/s/y' 3;: gig?” ' '1 . C! ”3‘1." {”41“ ".5 I,’ I. 4: wivdw 3?" LP“ IE 1 6 1899 JUI: 51:31 7-333“ 037433.319 23% STREAM WATER QUALITY AS RELATED TO URBANIZATION OF ITS WATERSHED BY Alvin Lee Jensen A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE Department of Fisheries and Wildlife 1966 ACKNOWLEDGEMENTS I thank Dr. Robert C. Ball for continuous financial support, and for giving me freedom to pursue my own goals. His dedication to the task of improving the quality of American life through preserving the quality of our natural bodies of water, and the energy with which he pursues this task will be a continued inspiration to me. I thank Dr. Frank Peabody for doing the bacteriological analyses appearing in this thesis. I thank my wife Joan for allow- ing me to spend most of my time in pursuit of my academic career. I thank my committee members and I thank all others who have in some way contributed to this work, for it could not have been completed without many of their suggestions. This work was supported by the National Institutes of Health through research grant 711111. ii "America today stands poised on a pinnacle of wealth and power, yet we live in a land of vanishing beauty, of increasing ugliness, of shrinking open space, and of an over-all environment that is diminished daily by pollution and noise and blight." Udall, 1963 iii TABLE OF CONTENTS Page INTRODUCTION . . . . . . . . . . . . . . . . . . . . . . 1 DESCRIPTION OF STUDY AREA. . . . . . . . . . . . . . . . 5 METHODS. . . . . . . . . . . . . . . . . . . . . . . . . 10 Sampling Procedures . . . . . . . . . . . . . . . . 10 Bacteriological Methods . . . . . . . . . . . . . . 12 Bioassay Procedures . . . . . . . . . . . . . . . . 15 Bottom Sampling Methods . . . . . . . . . . . . . . 15 Chemical Methods. . . . . . . . . . . . . . . . 14 Total Alkalinity as CaCO3. . . . . . . . . . . 14 Ammonia Nitrogen . . . . . . . . . . . . . . . 14 Biochemical OxygenODemand. . . . . . . . . . 15 Conductivity at 18 Centigrade . . . . . . . . 15 Nitrate Nitrogen . . . . . . . . . . . . . . . 16 Nitrite Nitrogen . . . . . . . . . . . . . . . 16 Organic Nitrogen . . . . . . . . . . . . . . . 16 Dissolved Oxygen . . . . . . . . . . . . . . . 16 Total and Ortho—phosphorous. . . . . . . . . . 17 Physical Methods. . . . . . . . . . . . . . . . . . 17 Total Residue. . . . . . . . . . . . . . . . . 17 Temperature. . . . . . . . . . . . . . . . . . 17 pH . . . . . . . . . . . . . . . . . . . . . . 18 Turbidity. . . . . . . . . . . . . . . . . . . 18 Velocity of River Waters . . . . . . . . . . . 18 Distance Between Stations. . . . . . . . . . . 18 Statistical Methods. . . . . . . . . . . . . . 19 RESULTS AND DISCUSSION . . . . . . . . . . . . . . . . . 22 Alkalinity. . . . . . . . . . . . . . . . . . . . . 20 Dissolved Oxygen. . . . . . . . . . . . . . . . . . 25 Biochemical Oxygen Demand . . . . . . . . . . . . . 56 Phosphorous . . . . , . . . . . . . . . . . . . . . 49 Summary of Nitrogen Determinations. . . . - . . . 54 Results and Discussion of Physical Tests. . . . . . 58 pH . . . . . . . . . . . . . . . . . . . . . . 58 Temperature. . . . . . . . . . . . . . . . . . 59 Total Residue. . . . . . . . . . . . . . . . . 59 Turbidity. . . . . . . . . . . . . . . . . . . 60 iv TABLE OF CONTENTS - Continued Page Conductivity. . . . . . . . . . . . . . . . . 69 Bacteriology. . . . . . . . . . . . . . . . . 70 Bottom Sampling Data. . . . . . . . . . . . . 79 RESULTS AND DISCUSSION OF DRAIN EFFLUENT ANALYSES . . . 86 SUMMARY . . . . . . . . . . . . . . . . . . . . . . . . 9O LITERATURE CITED. . . . . . . . . . . . . . . . . . . . 95 APPENDICES. . . . . . . . . . . . . . . . . . . . . . . 97 LIST OF TABLES TABLE Page 1. Summary of randomized block design analysis of variance and test for linearity for river total alkalinity determinations. . . . . . . . . . . . . 22 2. Data for diurnal oxygen determinations . . . . . . 27 5. Summary of one-way analysis of variance and test for linearity for river dissolved oxygen determin- ations . . . . . . . . . . . . . . . . . . . . . . 55 4. Data for determination of biochemical oxygen demand rate constant . . . . . . . . . . . . . . . 59 5. Summary of one—way analysis of variance and test for linearity for river station biochemical oxygen demand determinations. . . . . . . . . . . . . . . 42 6. Summary of one-way analysis of variance and test for linearity for river total—phosphorous determin- ations . . . . . . . . . . . . . . . . . . . . . . 55 7. Summary of one-way analysis of variance and test for linearity for river ortho-phosphorous determin- ations O O O O C O O O O O O O O O O O O O C O O O 55 8. Summary of one-way analysis of variance and test for linearity for river station nitrate nitrogen . 56 9. Summary of one-way analysis of variance and test for linearity for river ammonia nitrogen determin- ations O O O O O O O O O O O C O O O O O O O O O O 56 10. Summary of randomized block design analysis of variance for ammonia nitrogen determinations . . . 57 11. Summary of one-way analysis of variance and test for linearity for river organic nitrogen . . . . . 57 12. Summary of one-way analysis of variance and test for linearity for river pH determinations. . . . . 67 vi LIST OF TABLES - Continued TABLE 15. 14. 15. 16. 17. 18. 19. Page Summary of one-way analysis of variance and test for linearity for river turbidity measurements. . . 67 Summary of randomized blocks design analysis of variance for river temperatures . . . . . . . . . . 68 Summary of one—way analysis of variance for river conductivity determinations . . . . . . . . . . . . 71 Most probable number of gram negative bacteria per 100 ml. at river sampling stations. . . . . . . . . 76 Summary of bottom organisms . . . . . . . . . . . . 82 Summary of one-way analysis of variance and test for linearity for tubificidae . . . . . . . . . . . 85 Summary of one—way analysis of variance for Chironominae. . . . . . . . . . . . . . . . . . . . 85 vii LIST OF FIGURES FIGURE 1. 2. 5. 4. 5. 6. 10. 11. 12. 15. Map of study area. . . . . . . . . . . . . . . . . Average yearly discharge of the Red Cedar River at Farm Lane. . O C O C O O O O O O O O O O O O O O 0 Mean values for total alkalinity at river sampling Stations 0 O O O O C O O O O O O O O O O O O C O O Diurnal oxygen curves for 9-17—64. . . . . . . . . Diurnal oxygen curves for 8-15-64. . . . . . . . . BOD and dissolved oxygen curves. . . . . . . . . . Mean values for total—phosPhorous at sampling stations . . . . . . . . . . . . . . . . . . . . . Mean values for pH determinations at river Sampling stations. . . . . . . . . . . . . . . . . Mean values for temperature determinations at river sampling stations. . . . . . . . . . . . . . Mean values for turbidity measurements at river sampling stations. . . . . . . . . . . . . . . . . Mean values for conductivity determinations at river sampling stations. . . . . . . . . . . . . . Most probable number of gram—negative bacteria per 100 ml. at river sampling stations . . . . . . . . Mean number of Tubificidae per square foot at river sampling stations. . . . . . . . . . . . . . viii 24 29 51 48 52 62 64 66 75 78 84 LIST OF APPENDICES APPENDIX Page A. Drain effluent sampling data. . . . . . . . . . . 97 B. Description of drains entering river. . . . . . . 107 C. River station sampling data . . . . . . . . . . . 115 ix INTRODUCTION This study of a 2.65 mile reach of the Red Cedar River is part of an extensive program of studies on the river which began in 1958. Brehmer (1958) reported on a study of nutrient accrual, uptake, and regeneration. He indicated 100 ug. phosphorus per liter and 0.5 mg. inorganic nitrogen per liter are introduced into the river from Williamston's Sewage Treatment Plant during normal river flow. But these nutrients appeared to have been removed within a distance of 0.6 miles downstream. Vannote (1961) investigated the chemistry and hydrology of the river and its tributaries. Vannote reported a significant gradient of phosphorous enrich- ment with distance downstream. Grezenda (1960) stated run- off is the major source of inorganic nitrogen, while sanitary drains are the main source of phosphorous. Data of these workers indicates there is a significant amount of pollution entering the river. Determinations made at Farm Lane indi- cated 9.46 metric tons of nitrogen are flushed downstream annually. This does not take into consideration the contribu- tion from East Lansing's Kalamazoo Street Sewage Treatment Plant. King (1964) states Fowlerville plating plant wastes have eliminated all aquatic macrofauna, except Tubificidae, for a distance of 15 miles downstream from the plating plant. King also reports inorganic sediments from highway construc- tion have reduced Red Cedar aufwuchs production 68 percent. The study described here was motivated by visible signs of stream deterioration between Hagadorn Road and East Kalamazoo Street. Objectives of this study were to locate pollutant sources, determine the amounts entering the river and to evaluate effects of these pollutants on water quality and on composition of river bottom fauna. For clarification, pollution as used in this study, is defined as every con- tamination or alteration of physical, chemical, or biological properties of a lake or stream which will or may lower its value to man. Water pollutants can be separated into three categories: industrial pollutants including pesticides, erosional products, and domestic sewage. Determinations made during this study were directed towards establishing the level of domestic and erosional pollution. These were the most obvious sources of pollution. This is not to underestimate the importance of other pollutants, of which there are many. In this study an effort has been made to integrate biology, chemistry, and sanitary engineering, for the problem of water pollution is not completely encompassed by any one of these disciplines. Perhaps water pollution studies are best interpreted from the perspective of human ecology--the study of the relationships between man and his environment. DESCRIPTION OF STUDY AREA The Red Cedar River is a central, lower Michigan warm- water stream arising from a marshy area near Cedar Lake in Marion Township. The Red Cedar River flows 40 miles through Michigan farmland and woodlots, flowing by several small towns, through city parks, and through the campus of a large university before its confluence with the Grand River at Lansing, Michigan. King (1964) describes aspects of the river biology, Vannote (1961) describes the river hydrology and phosphate chemistry, and Meehan (1958) describes the river climatology. This study was completed on the river reach located be- tween Hagadorn Road and East Kalamazoo Street. The study section is 2.65 miles in length. Stream width varies from 25 to 80 feet. Water depths vary from 15 inches over silt— covered sand flats to approximately 6 feet in the impoundment formed by a dam located at the study section mid-point. The dam was constructed primarily to provide a reservoir to be used as a source of cooling water for Michigan State University power plants. A secondary use of the impoundment is for recreation, but poor water quality presently limits this use. During this study period from August 1964 to November, 1964, volume flow of the Red Cedar River was considerably less than the average summer flow from 1947 to 1965. 5 .mCOHumum mQHHmEmm Hm>HH mo mGOHumooH mCH30£m .mwum wvsum mo mm: .H mudmflm A! - I , o to 2.. I. . \ .1 _ 4 .5: 4o cow. 8.. ave n 2255 I 2.23.. U «W 2.5.2.22: 3525 035.2322 M A. o @2523 55 T mm 3 2.33 o o 0 335 Foo... __ u m a 20.25 L A K T m «é E « a. «9 .w 6o o m 29.25 _ 20.25 t 22.25 02429 ¢u7r¢ rhea: Figure 2 illustrates a dramatic decrease in volume flow of the Red Cedar during the past several years. A correlation study of ground water, precipitation, and cfs flow of the Red Cedar indicated decreased flow is poorly correlated with both precipitation and ground water. Regression coefficients obtained were 0.55 and -0.16 respectively. A more complete analysis of available data is necessary. Bottom sampling showed the river bed at Hagadorn Road to be coarse gravel covered with two or more inches of silt. Silt accumulation is thought to be caused by decreased carry- ing capacity of waters entering the impoundment where velocity decreases, and ability to carry particles in suspension diminishes.r Decaying leaves and other detritus blanket the stream bed of the large pond-like area upstream from Farm Lane Bridge. Above East Kalamazoo Street, below drain 58, a large sludge bed has formed, and extends downstream for nearly a quarter mile. Below the Kalamazoo Street Sewage Treatment Plant sludge beds several feet thick have formed, and in some locations nearly protrude from the water surface. Recent years have brought a large increase in construc- tion along the banks of the Red Cedar River between Hagadorn Road and East Kalamazoo Street. New apartment buildings and university facilities are numerous. Construction of these buildings, roads, and bridges has resulted in the flushing of a huge sediment load into the river. Rapid growth of Figure 2. Average yearly discharge for the Red Cedar River for 1947-1964. (Data from U. S. Geological Survey.) Regression equation calculated for data is, § = 550.81 - 13.95 x Y = Estimated discharge in cubic feet per second. X = Time in years from 1947 = 0. mom. ”mm. .mm. mom. hum. mmm. ”mm. .09 me. NVm. 0 0 0 0 0 ozoowm mun. Puma 0530 z. mom._._Z_.._x0 ou>..0wm.o amt... cum .0! STATION 5 3.00 200 LOO 3(33 840nm 2xo om>.._omw_o amt]. «ma .0: Ida“? 8(Oflm 3(09N. 8a0fl0. Suofin Snoflo Snofl¢ Inoflu Incas— 8(Oflo. 3 K> n Miles downstream from Hagadorn Road, 54 indicates a large decrease in oxygen percent saturation as the waters flow downstream. A regression analysis, Table 5, indicated the existence of significant differences between means at sampling stations. Scheffe's test for significance of differences between individual means indicated the mean for Station 5 and the means for Stations 5, 2, and 1 were significantly different. Also the difference between means for Stations 1 and 4 was significant. Oxygen concentrations decreased significantly from one downstream station to the next, with the average value at Hagadorn Road being only 75 percent saturation. Values at all downstream points were below levels considered permissible for maintaining water quality. Large oxygen deficits are caused by the increased depth and low turbulance of the impoundment, by sludge accumulation, by decomposition of leaves and litter, and to a major extent from the stabilization of large quanti- ties of sewage entering the river. Addition of municipal pollution to a stream has a dramatic effect on dissolved oxygen concentrations, as bacteria consume oxygen in the process of stabilizing decomposable organic matter abundant in municipal sewage. Municipal sewage also has a high nitrogen and phOSphorous concentration resulting from the breakdown of proteins, amino acids, and household detergents. To illustrate quantitatively the relationship existing between dissolved oxygen and concentrations of components of 55 Table 5. Summary of one-way analysis of variance and test for linearity for river dissolved oxygen determinations. Source SS df MS Fexp. F0.95 Treatment. 19746 4 4956.50 18.46 2.58 Regression 18959 1 18959.00 71.61 4.04 Deviation 787 5 262.55 0.99 2.80 from linearity Error 15501 50 264.72 Total 55247 54 56 domestic sewage effluents, phosphorous, nitrogen, and BOD, a multiple linear regression analysis was made. The following equation was fitted to the data by the method of least squares (Li, 1964), § = 81.50 - 0.47xl - 25.11x2 - 1.78x3 where, Y = Estimated mg. per liter of dissolved oxygen for a given set of (X1, X2, X3). X1 = Mg. per liter Biochemical oxygen demand. X2 = Mg. per liter total phosPhorous. X3 = Mg per liter nitrate + ammonia nitrogen. To test how well this regression equation predicts dissolved oxygen concentrations given the BOD, total phOSphorous, and nitrogen values, a correlation analysis was made between values predicted by this equation and values observed in this study. A correlation coefficient R, of 0.85 was obtained. This high value indicates phOSphorous, nitrogen, and BOD concentrations are highly correlated with decreasing oxygen concentrations. And it can be concluded sewage entering the river is the factor most reSponsible for observed oxygen concentrations. This brings us to another important topic, biochemical oxygen demand and stream self purification. Biochemical Oxygen Demand Biochemical oxygen demand, BOD, is defined as the amount of dissolved oxygen required by bacteria to stabilize 57 decomposable organic matter in a volume of water under aerobic conditions (Sawyer, 1960). An Englishman, Sir Edward Frankland, formulated the BOD test in 1870. Since then an impressive theory has been established describing the kenetics of the oxidative processes which occur in streams receiving sewage. Oxidation of organic matter in water takes place in two stages, requiring about 20 days to complete. The first stage involves oxidation of carboniferous material. The second stage is oxidation of nitrogenous material by the bacteria Nitro- somonas and Nitrobacter as follows, Nitrosomonas 2NH3+502 2NO2-+2H++2H20 Nitrobacter + Z‘Nog:+o,2 )2NO2_+2H To avoid the error which would be introduced by the second stage of this process, the standard BOD determination is run for only 5 days, and measures only 70 to 80 percent of the total BOD which is termed ultimate BOD (Phelps, 1944). The BOD analysis is essentially a bioassay procedure, and as such the presence of trace amounts of toxic materials like chlorine, copper and arsenic lead to underestimation of BOD values. Chlorine concentrations in the river were high. If the oxygen concentration is not limiting, decomposition of organic matter in streams can be considered to be a first order chemical reaction, with the reaction rate depending only on concentrations of decomposable organic material. The first stage of the BOD reaction can be formulated in the differential 58 d(BOD) = dt k is the BOD reaction rate constant (Phelps, 1944). Solution equation for a first order reaction, - k(BOD), where of this equation gives, BODt = BODO e—kt, where BODt is the BOD remaining at time t, BODO is the initial BOD, and t is again time. Engineers often use the retarded exponential, BODt = L(1 - e'kt the ultimate BOD. At the present time there is no satisfactory ), but this formula requires evaluation of L, way to obtain L (Fair and Geyer, 1961). Using the simple exponential, determination of the BOD reaction rate constant k allows prediction of the residual BOD at points a known time t of flow downstream from the point where BODO was determined. To estimate the number of persons polluting the river between Station 1 and Station 5, the BOD reaction rate con- stant was determined, and an estimate of river velocity was made. Results are reported in Table 4. For each station two sets of 5 samples each were started. Two samples were termi- nated after 1, 2, 5, 4, and 5 days. If the BOD reaction can be treated as being linear, the regression of ln %§%%§%- on time should not deviate significantly from linearity. The data for each station were tested for linearity, and the linear regression equation was derived. Only at Station 1 did the data deviate significantly from linearity. The following equations were derived. Station 1: § = 0.64 + 0.06 t Station 2: Y = 0.61 + 0.05 t Station 5: Y = 0.49 + 0.12 t Table 4. Data for determination of BOD reaction rate constant. Data expressed as ln(BOD). Station Hours before termination T 12 24 56 48 96 'J 1 1.000 0.587 0.875 0.215 1.197 1.000 1.000 1.000 0.651 1.787 Ti' 2.000 1.578 1.875 0.864 2.984 9.285 2 0.059 0.695 1.001 0.770 0.524 -0.128 0.715 0.565 2.007 1.571 Ti' -0.089 1.206 1.566 2.777 1.895 7.555 5 0.695 0.519 0.859 0.015 1.165 1.058 0.940 1.482 1.517 1.778 Ti' 1.751 1.459 2.541 2.552 2.941 11.024 4 -0.255 0.589 0.806 0.588 1.258 -0.529 0.675 1.151 0.610 1.609 Ti' -0.562 1.262 1.957 1.198 2.847 6.702 5 2.079 2.216 2.407 2.518 2.791 1.717 2.255 2.965 2.515 2.708 5.796 4.469 5.572 4.851 5.499 25.967 T.. 1 40 Station 4: Y = -1.16 + 0.57 t Station 5: § = 1.48 + 0.19 t where, Y = Estimated ln (BOD ) . (BODO) t = Time in hours. The average value of the slopes of these lines was calcu- lated to be 0.15. This compares favorably with the value of 0.17 given in the literature (Sawyer, 1960). The average expected BOD value at station 5 was calculated to be 5.57 using -(0.15)(O.97) the formula BOD = (5.69) e . The value for 0.97 mg. per liter BOD added between Stations 1 and 5 was determined to be 7.05. During the period of this study the U.S.G.S. reports 16.7 cfs as the average flow of the Red Cedar River at the Farm Lane gauging station. River velocity was estimated as 0.05 feet per second. It has been shown that the mean BOD of sewage produced per person per day is 54,000 mg. (Fair and Geyer, 1961). It must then be concluded from the following calculations, 16.7 cfs = 56,072,000 liters per day. Mg. per liter entering Red Cedar per day. (56,072,000)(7.05) = 252,504,000 Persons polluting Red Cedar in study area. 252,504,000 54,000 = 4'674 that 4,674 persons are polluting the Red Cedar River between Hagadorn Road and West Kalamazoo Street. Because of the toxic 41 effect of chlorine introduced at the Kalamazoo Street Sewage Treatment Plant, this value is a minimum. This method of determining the extent to which a river reach is polluted is in some ways unsatisfactory. First, the BOD reaction is not a first order reaction, but rather a mix- ture of first, second, and third order reactions. Secondly, the deoxygenation rate constant as determined in a BOD bottle is probably different from the stream deoxygenation rate constant. And finally, the time of flow between stations is difficult if not impossible to determine. An empirical method can surmount these difficulties. Analysis of variance, Table 5, indicated data had significant curvature, and that significant differences existed between sampling station means. Under any set of environmental con- ditions BOD is some function of distance downstream from our first sampling station. By Taylor's theorm, this function can be approximated by a polynomial. The following least squares regression equation was calcu- lated (Li, 1965). Y=5.8+5.4x-7.8x2+2.5x3 Y = Estimated mg. per liter BOD. X = Distance downstream from Station 1. A correlation coefficient R = 0.66 was obtained between ob- served values and values predicted. This cubic regression equation can be used to estimate mg. per liter BOD added to 42 Table 5. Summary of onedway analysis of variance and test for linearity for river BOD determinations. Source SS df MS Fexp. F0.95 Treatment 251.69 4 62.92 10.81 2.67 Regression 66.99 1 66.99 11.51 4.15 Deviation from 184.70 5 61.56 10.57 2.90 linearity Error 186.28 52 5.82 Total 457.97 59 45 the river in the study section. Differentiating gives the differential equation, 93:. = _ 2 dX 5.4 15.0 X + 7.5 X which expresses instantaneous change of BOD with respect to distance downstream as a polynomial function. Solving this equation gives, 0 2. Y = f 65(5.4 - 15.0 x + 7.5 x2) dX = 7.05 0 which is nearly the same value obtained by the Streeter—Phelps method. In larger stream reaches, and more heterogeneous stream reaches the values would probably not be so close. The method presented here is thought to be superior to the standard engineering method in that it does not involve de- termination of the stream deoxygenation constant and time of flow between sampling stations. This same procedure can be used for studying other stream parameters, for example, phOSphorous or nitrogen. Small increments of organic material added to a stream are broken down, by bacterial action, into simple inorganic molecules. Oxygen used by bacteria in stabilizing these organic materials is replaced through physical and photo- synthetic reaeration. This process is termed natural stream purification. Engineers have long sought a means for de- termining what level of BOD can be loaded into a stream with- out depressing the oxygen to undesired levels. Streeter and 44 Phelps proposed the following differential equation as a model for the process of stream reaeration and deoxygenation (Streeter and Phelps, 1925). d(DO) dt = K1 (BOD) - K2 (D0) Solution of this equation gives, K1(BODQ) (e-Klt _ e-th K2t + _ K2 _ K1 ) DOOe DO= which can be used to predict DO deficits which will be in- curred by a BODO loading a known time of flow downstream. Application of this equation requires determination of K1, the deoxygenation constant; and K2, the reaeration constant. Both of these constants are difficult and in some cases perhaps impossible to obtain accurately. K1 is the rate of chemical and biological removal of oxygen from the water. An estimate of K1 can be obtained from the BOD reaction rate constant, but this does not take into consideration the oxygen demand of bottom sludges, or the respiratory needs of stream dwelling plants and animals. Velz (1958), has shown that bottom sludges have a considerable oxygen demand. K2, the reaeration co- efficient is equally difficult to estimate. Streeter and n Phelps (1925) have empirically estimated that K2 =-%¥— , where, C = Surface slope and channel irregularities factor. V = Velocity of water movement. n = Factor derived from relationship of river stage to velocity. H = Average water depth. 45 O'Conner and Dobbins (1956), have derived theoretically the equation, 480 D5 S5 H K2: where, D = Coefficient of molar diffusion. S Stream slope. H Average depth. From both of these equations it can be deduced that increased depth and decreased slope have a tremendous effect on stream oxygen dynamics. Moore 2; gl. (1950) have presented a method for obtaining these constants using nomograms. Nemerow (1965) gives the constants simple formulas which, however, yield results of questionable accuracy. But use of this simple linear model is also a question— able procedure. Dobbins (1964) has found this method fails to consider the following important factors. 1. Removal of BOD by sedimentation or adsorption. 2. Addition of BOD by scour of bottom deposits or by diffusion from bottom deposits. 5. Addition of new increments of BOD from other sources. 4. Removal of DO by diffusion to benthal deposits to satisfy benthal oxygen demand. 5. Addition of oxygen by plants. 6. Removal of DO by respiration of aquatic organisms. 7. Changes in channel configuration. 8. Diurnal fluctuations in BOD, DO, temperature, etc. 46 It must be concluded that although the Streeter—Phelps model is used in nearly all engineering surveys, it is far too simple to accurately describe reaeration and deoxygenation of a natural stream. Churchill and Buckingham (1956) have presented a statis- tical method for evaluation of stream loading capacity. Using the multiple-linear regression equation, Q = bO + blxl + b2x2 + b3X3 where, Y = D0 in ppm. at sag point in stream DO curve X1 = 5 Day BOD in ppm. X2 = Water temperature X3 = Stream discharge, it was possible to accurately predict the maximum downstream DO sag caused by a load of organic pollution. But this method is not applicable if more than one source of pollution exists in the stream reach being studied, or if the stream reach is not homogeneous. In some cases the natural purification capacity of a stream may be estimated using the following procedure. BOD and DO can be estimated as functions of distance downstream from the first sampling station. This gives two parametric polynomial equations, the graphs of which are shown in Figure 6. BOD = fl (distance) D0 = f2 (distance) Figure 6. 47 Means and regression equations for percent saturation dissolved oxygen and mg. per liter biochemical oxygen demand. Distances on abscissia are proportional to distances between stations. Observations were taken at approxi— mately 8:00 A.M. Each oxygen value is the mean of 55 observations. Each BOD value is the mean of 16 observations. 6 O BIOCHEMICAL OXYGEN DEMAND 90 O PERCENT 8 ATURATION OISSOLVED OXYGEN oz<2wo zmo>xo 129291005 EMF... «$1.02 wmwmnm 9876543 0 .. O o o o o o o o 8 7 6 5 4 3 2 zwm>x0 om>40wm5 ZO_._.._._o_mm3._. zomxos. SAMPLING STATION Table 12. Summary of one-way analysis of variance linearity for river pH determinations. 67 and test for Source SS df MS Fexp. F0.95 Treatments 0.87 4 0.21 5.00 2.65 Regression 0.04 1 0.04 0.57 4.12 Deviation from 0.85 5 0.27 5.85 2.88 linearity Error 2.67 55 0.07 Total 5.65 59 Table 15. Summary of one-way analysis of variance and test linearity for river turbidity determinations. for Source SS df MS Fexp. F0.95 Treatments 2885.26 4 720.81 17.57 2.69 Regression 426.17 1 426.17 10.59 4.17 Deviation from 457.09 5 152.56 5.71 2.92 linearity Error 1260.45 50 41.01 Total 4145.69 54 f““‘-_f?1 u-‘.-‘4._ -a—. _ ) 68 Table 14. Summary of randomized block design analysis of variance for river temperatures. Source SS df MS Fexp. F0.95 Treatment 266.25 4 66.55 4.88 2.47 Blocks 2751.85 25 110.07 8.08 1.60 Non-additivity 5.91 1 5.91 0.28 5.95 Remainder 1558.52 99 15.72 Error 1562.25 100 Total 4580.29 129 69 Conductivity» Ohm's law for metallic conductors states that I ='%, where R ='EL , is termed the resistance. p Is the Specific A resistance or resistivity of the conducting medium, L the length of the resistor, and A the cross sectional area of the resistor. As soon as sufficiently sensitive measuring devices became available it became apparent that solutions of electro- lytes as well as metallic conductors followed Ohm's law. According to the Debye-Huckel theory, ions in solution migrate towards the pole with electrical charge opposite to their own when placed into an electric field. This migration creates a measurable electric current proportional to concentrations of the dissolved ions. Conductance, conductivity, and molar equivalents concentration of solute are related by the expres- sion (Moore, 1962), Conductivity = the reciprocal of resistivity. Measured as (ohm-l cm'l). where, F -/\-= Conductance = the reciprocal of resistance. Measured as (ohm-1). C eq Molar equivalents solute. It can be seen from this equation that as concentrations of solutes increase conductivity of the solution increases. An approximation of the mg. per liter of cations or anions —1 in a solution can be obtained by multiplying u-ohms conduct- ance by 0.01 (Standard Methods, 1960). 70 Analysis of conductivity data for stream sampling sta- tions indicated significant differences existed between sampling station means, and that data do not deviate signifi— cantly from linearity, Table 15. The linear regression equation, §= 59.1 + 8.0x Y = Estimated conductance in ohm"1 x 10‘3 X = Miles downstream from Station 1, indicates a consistant large increase in conductance at down- stream sampling stations. Results of Scheffe's test indi- cated means for conductance measurements at Stations 4 and 5 were significantly different from that at Station 1. The mean for Station 5 is also significantly different from the mean for Station 2. As every box of salt purchased in the East Lansing com— munity eventually finds its way down the drain, it is not dif- ficult to account for the pattern of conductivity measurements. illustrated in Figure 20. Run-off may also account for some increase in conductivity. But Ellis (1956) states that erosion silt does not significantly alter the salt complex, or the amount of electrolytes in river waters. Bacteriology Water has long been known to be a possible transmitter of disease. Hippocrates recommended drinking water be filtered and boiled (Frank, 1955). Typhoid fever, paratyphoid fever, 71 Table 15. Summary of one-way analysis of variance and test for linearity for river conductance determinations. Source SS df MS Fexp. F0.95 Treatment 0.227 4 0.056 8.00 2.69 Regression 0.190 1 0.190 27.14 4.17 Deviation from 0.057 5 0.012 1.71 2.92 linearity Error 0.215 50 0.007 Total 0.442 54 72 Figure 11. Conductivity determinations for river sampling stations expressed as mho's. Distances on abscissia are proportional to distances between sampling stations. Each value is the mean of 5 observations. 800 700 wOI 2 600 z. 0 0 0 0 5 4 >._._>_._.ODQ 20 O O 3 o 200 SAMPLING STATION 74 bacillary dysentery, amebic dysentery, asiatic cholera, diarrhea, infectious hepatitis, infectious jaundice, polio- myelitis, and anthrax are some diseases commonly associated with water. Unpolluted streams contain an abundant bacterial popu- lation averaging, for example, 1900 per 100 ml. in a Great Britain survey of twenty English rivers (Gainy and Lord, 1952). But the number of free living bacteria is usually only a dozen or so per 100 ml., composed primarily of Species of Micrococcus, Flavobacterium, Achromobacterium, Bacillus, Proteus, and Leptospira (Frobisher, 1957). Surface run-off contributes considerable numbers of bacteria to streams, causing bacteria counts to be significantly higher after heavy showers (Frobisher, 1957). Cultivated soils contains up to several billion bacteria per gram. It has not been possible to devise specific tests for water-borne disease producing micro-organisms. It is there- fore common procedure to test for relatively high numbers of coliform bacteria rather than to attempt tests Specifically for Salmonella typhi, the typhoid producing bacilli, or other disease producing organisms. The coliform group of bacteria includes all aerobic or facultative anerobic Gram-negative non-Spore forming bacilli capable of fermenting lactose broth with gas formation (Gainy and Lord, 1952). Although most coliforms are harmless, high populations are associated with the presence of pathogenic 75 bacteria and water containing large numbers are regarded as dangerous (Gainy and Lord, 1952). The United States Treasury Department prohibits inter- state common carriers to serve drinking water with more than one Escherichia coli (a common coliform bacteria) per 100 m1. Sewage plant effluents usually contain between 10,000,000 and 1,000,000,000 Escherichia coli bacteria per 100 ml. prior to dilution (Gainy and Lord, 1952). The most probable number (MPN) is an index of the number of coliform bacteria which, more probably than any other number, would give the results shown by the laboratory examination (Standard Methods, 1960). As can be seen from Figure 12, the most probable number of coliform bacteria in the Red Cedar River is considerably higher than is desirable. At one point the bacteria count exceeds a value commonly found in sewage plant effluents. Bacteria counts upstream at Hagadorn Road are excessively high, indicating, as does BOD and dissolved oxygen data, that con— siderable pollution is entering the river above this point. The source of these upstream pollutants may be saturated septic tank leaching fields. As the river flows through the Michigan State University campus the MPN of coliform bacteria increases exponentially on logarithmic graph paper, indicating several sources of domestic sewage enter the river in the study area. 76 Table 16. Most probable number of gram negative bacteria per 100 ml. of river water for five sampling stations. . Dilution Statlon 10 10'1 10"2 10‘3 10'4 10‘5 10-: 10"7 MPN 8-24-84 1 5/5 2/5 1/5 0/5 0/5 15000 2 5/5 2/5 1/5 0/5 0/5 15000 5 5/5 5/5 2/5 0/5 0/5 15000 4 5/5 5/5 2/5‘ 0/5 0/5 110000 5 5/5' 5/5 5/5 2/5 0/5 1100000 9-15-64 1 2/5 0/5 0/5 0/5 0/5 0/5 9100 1/5 0/5 0/5 0/5 0/5 0/5 5800 5 5/5 1/5 1/5 0/5 0/5 0/5 45000 4 2/5 0/5 0/5 0/5 0/5 0/5 9100 5 5/5 0/5 0/5 0/5 0/5 0/5 25000 11-18-84 1 5/5 5/5 5/5 0/5 0/5 0/5 2400 5/5 5/5 2/5 0/5 0/5 0/5 9500 5 5/5 5/5 5/5 5/5 2/5 0/5 450000 4 5/5 5/5 5/5 0/5 0/5 0/5 24000 5 5/5 5/5 5/5 5/5 5/5 0/5 2400000 77 Figure 12. Most probable number (MPN) for coliform bacteria per 100 ml. of river water at sampling stations. Distances on abscissia are propor- tional to distances between sampling stations. 0 0 O 2 n 0 00 o 432 32 I I000 300 200 0 O. x 42 00. «in. <_mm._.0_._.<0wz 2410.05 mmiaz m4m