|
International Journal of Environmental Research, Vol. 4, No. 1, 2010, pp. 153-160 Prediction of Dissolved Oxygen through Mathematical Modeling Naik, V. K. 1 and Manjapp, S.2 1Department
of Built & Natural Environment, Caledonian College of ngineering, Seeb,
Muscat, Oman
Received 22 Dec. 2008; Revised 13 Aug. 2009; Accepted 25 Aug. 2009 Code Number: er10016 ABSTRACT Water quality aspect with regard to dissolved oxygen was studied for a 24 kilometer stretch of Malaprabha River in Karnataka State, India. Main objective of the research was to simulate the predicted dissolved oxygen depletion due to the waste load allocation in the river with the ambient observed values of dissolved oxygen, and to ascertain the application of mathematical modeling for predicting the dissolved oxygen in a mixing zone. The entire river stretch selected for the study was divided in to four stations and water samples were collected from each station, and analyzed for different parameters. Station S1 was purely the upstream of the sewage discharge point and hence reflected the water quality of the river without pollution. Station S2 and downstream stations were the part of mixing zone where the discharged sewage at S2 kept spreading. Dissolved oxygen levels at these four stations were calculated through mathematical modeling, considering all variables that affect dissolved oxygen variation. Further, this predicted dissolved oxygen found with the application of mathematical modeling was then simulated with the actual observed values of dissolved oxygen in these four stations. The results were highly encouraging, with the predicted values almost agreeing with the observed values. Thus the application of such mathematical modeling may be useful for a river stretch where water pollution due to the discharge of effluents is expected, and will help the engineers in arriving at a proper waste-load allocation for the river. Key word: Mixing zone, In-stream, Transects, Sources and sinks,Diffusion factor INTRODUCTION Analytical solutions for steady state dissolved oxygen (DO) prediction in mixing zones of rivers are constructed by employing the product law concept. Dissolved Oxygen (DO) is an important water quality parameter affecting the health of a river and hence, great importance is attached to maintain the DO at desirable level. Dissolved oxygen is important to aquatic life, as detrimental effects can occur when DO levels drop below 4 mg/ lt. or 5 mg/ lt. depending on the aquatic species (Chapra, 1997; Thomann and Mueller, 1987). Oxygen level in Malaprabha River is adversely affected by the discharge of untreated sewage from a town on its bank in to the river. Upstream part of the river water quality before the discharge of sewage is fairly good, and the oxygen level is mainly controlled by the natural forces in this stretch. An effort is made in this paper to predict dissolved oxygen in the mixing zone due to the impact of sewage disposal from the nearby town. Dissolved oxygen of a stream plays an important role in the river eco-system affecting the flora and fauna of the stream (Wang et al., 2003). Generally, a single point measurement of dissolved oxygen by the traditional method is not reliable, and sometimes may be misleading in assessing the state of the river water body. This is due to the fact that it does not reflect the effect of biological activity and also do not reveal extremes in dissolved oxygen concentrations. The aim of the study was to arrive at a mathematical solution through modeling which would be able to predict the dissolved oxygen variation with distance. Attention was also focused on the lateral distribution of dissolved oxygen across the river width. The cross section of the river channel within the mixing zone was divided in to two parts, viz., Limited Use Zone or Zone of non-compliance –LUZ and Zone of Passage –ZOP, requiring water quality compliance along designated LUZ boundary, as shown in Fig.2 (Thomann, 1987; Chapra, 1997). The location of the Malaprabha River is as in (Fig.1). MATERIALS & METHODSField tests were conducted from June, 2005 to May, 2006 to monitor the water quality parameters at 4 Stations – S1, S2, S3 and S4 as shown in Fig.2. The study stretch selected was 24 kilometer, starting from its birth place. A part of the stretch which formed the upstream of the sewage disposal point was 5 km, and the downstream stretch from the discharge of sewage was 19 kilometer and sewage from the town is discharged into the stream through a piped outfall located at the left bank (looking d/s), as shown in Fig.2. The channel alignment is fairly straight, but significant variations occur in widths and depths. The upper portion of the stretch is shallow with depths up to 1m, with rocky bottom, whereas the lower portion is deep, up to 2 to 2.5 meter, with muddy bottom. The in-stream and sewage characteristics were analyzed during this period. The method of sample collection, transport to the laboratory, preservation and analyses were all done as per the standard procedures(American Public Health Administration, 1998). The in-stream monitoring was carried out for four to five lateral points at each station. Water temperature and DO were measured on site only during the sampling period. The analyzed result of sewage for important parameters, re-oxygenation and deoxygenation rate constants, respiration and photosynthesis rates and other hydro-geometric characteristics of the river stretch are presented in Table 1. The decay and re-aeration rates for CBOD and DO, respiration and photosynthesis rates along with sediment oxygen demand – were all found by applying standard methods and procedures (Thomann, 1987; Chapra, 1997). The measured DO profile just at the u/s of sewage disposal and the downstream for S1, S2 and S3 are presented in Fig.3. which exhibit both lateral and longitudinal variations. The DO values along the outfall (left) bank at S1, S2, S3 were found to be lower than those at the opposite bank, displaying increase in the lateral direction, clearly indicating the plume behavior of the waste dispersion. The other stations show relatively uniform distribution of DO laterally. The analytical formulation for 2-D DO model assumes that, steady state conditions prevail; cross sectional distribution of background water quality constituents are uniform, densities of sewage and receiving water are same, far-field concentrations are not affected by near-field mixing processes, longitudinal dispersion effects are negligible. The DO distribution within the mixing zone is affected by pollutant loadings from two principle inputs, i.e. the upstream flow and the sewage effluent discharge. If the initial DO is lower than DO saturation value, there would be initial DO deficit in each input (Morgan et al., 2006). The oxygen sinks considered here, include the two DO deficits, carbonaceous biochemical oxygen demand (CBOD) and nitrogenous oxygen demand (NOD) that are present in the upstream water and outfall discharge, while the DO source is due to re-aeration. The formulation for DO distribution in mixing zone consists of three steps; I) The 1-D relationships for DO deficits due to upstream loadings and distributed internal sources and sinks are considered assuming cross sectional uniformity. II) 2-D relationships describe the effect of sewage discharge which include:
iii) The principle of superposition is invoked to obtain the total DO deficit as the sum of the 1-D and the 2-D deficits (Gowda, 1980). The DO deficit, Db due to the combined effects of the decay of background CBOD, NOD and initial DO deficits, and of re-aeration are described by the 1-D transport equation (Chapra, 1997). where: Db = (Cs - Cb), inmg/ lt. The solution for eqn.(1) is: Where, Lbo& Nbo are the initial values of background CBOD & NOD. ) Db= (Cs-Cb) in mg/ lt. Dbo = Cs -Cbo = Initial background DO deficit at x =0 Cbo = Initial background DO value The diffusion factor Fc is given by: And the factor Fn is given by: The concentration of NOD is calculated as NOD = 4.57×TKN +1.14× NO2 (5) TKN = Total Kjeldahl Nitrogen; and NO2 = Nitrogen Dioxide in the water sample. For the convenience, the effects of effluent DO and of CBOD & NOD on the in stream DO deficit are considered separately. The readily available 1-D &2-D equqtions3 are utilized to construct analytical solution for the 2-D deficit distribution in mixing zones through the application of product law concept, as outlined below: i. Effect of Effluent DO: The in-stream DO deficit due to the effluent DO deficit may be described by the 2-D transport equation: where: De = DO deficit at a point (x,q) in the mixing zone due to the effluent DO deficit Deo Eq = Transverse diffusion factor q =Cumulative discharge denoting the lateral co-ordinate mx = Scaling factor The solution to eqn. (4) is given by: Qe = Effluent flow rate And G (Φ, p) = A factor depending on the type of outfall and location and equal to 1 for steady state condition ii) Effect of Effluent CBOD: The 2-D distribution of CBOD in the stream due to the effluent CBOD is given by: wherein, L = CBOD concentration at a point (x, q) Solution for the above equation for an effluent CBOD value Le is The 2-d distribution of DO deficit in the stream channel due to the combined effect of decay and lateral dispersion of CBOD, and re-aeration, is described by the transport equation in which, Dc =DO deficit at a point (x, q) due to CBOD.Solution to eqn. (10) is: By following the similar procedure as above for CBOD, we can obtain the expressions the distributions of NOD, N and the corresponding DO deficit, Dn can be obtained, by replacing N, Ne, Kn and Fn in place of L, Le, Kd and Fc respectively in above equations. Total DO Deficit: The distribution of total DO deficit, Dt due to the combined effects of upstream inputs, internal source and sinks, re-aeration and effluent loading is described by the transport equation: where: Sp =(S + R - P) denoting the algebraic sum of the distributed sinks and source. Dt = (Db + Dp + Ds) Ds = De + Dc + Dn Dp is the DO deficit due to Sp given by equation The analytical solution for the equation (12) is given by the sum of the DO deficits of the corresponding solutions, which can be written as: In the above equation, the terms within the three sets of the square brackets on the right hand side denote, in order, the DO deficits due to the background, distributed source & sinks and effluent. However, the effect of Sp is not considered in this article for sake of simplicity.The predicted values for the stream under the study, is estimated with the use of the above equations and the values are as shown in Table 2. These predicted and the observed values of Dissolved Oxygen are presented in Table 3. The comparison between the observed and the predicted DO profile for the selected study area of Malaprabha River is simulated as shown in Fig.4. RESULTS & DISCUSSION Analyzed results of sewage sample for important parameters, re-oxygenation and deoxygenation rate constants, respiration and photosynthesis rates and other hydro-geometrical characteristics for the study area are presented in Table 1. Measured dissolved oxygen profile just at the upstream of sewage disposal (S1) and the downstream stations S2 and S3 is presented in Fig.3, which exhibits both lateral and longitudinal variations (Tiseer et al., 2008). Due to gradual spreading of the effluent, lateral concentration gradients were established in the mixing zone. This had resulted in spatial distribution of DO as shown pictorially in Fig.3, with minimum DO distribution occurring along the outfall bank, while it was near to saturation along the opposite bank. The distribution factor, however, would depend on the type of outfall, i.e. piped or diffused –and its location in the river. Further, at first two stations, i.e. S1 and S2, the DO concentrations near the opposite bank was close to those at the upstream station S1, but lower DO values at the other two stations –S3 and S4. This clearly indicated the plume spreading effect as in Fig.2. where extent and the distance depending on the velocity and nature of flow of the river water (JCE, 2005; Vellidas et al., 2006). Values of predicted dissolved oxygen deficit are presented in Table 2. Steady states analytical DO model for the mixing zone of Malaprabha River was applied and these DO deficit results were obtained. Actual field DO values were also observed at these stations, and both the values are presented in Table 3. It can be observed from Fig. 4 that the observed and predicted values almost agreed, and the predicted DO values showing slightly higher results. As the DO variations occur due to many factors, such as biological activities in the river, re-aeration and deoxygenation rates etc., -a 2 to 3 percent variation in DO variation is acceptable (Wang et al., 2003). The results indicated that the concentration of DO was less at the effluent disposal point and also along the bank, whereas it was more near the opposite bank. Dispersion coefficient for biochemical oxygen demand varied from 0.151 to 0.07, with 0.151being highest at station S2, indicating the maximum spread of the effluent at the disposal location. Reduced values of dispersion coefficients were observed at other stations, which showed that the effluent dispersion was on lower side due to the effect of dilution (Steele, 1989). Similar observations can be made with other parameters presented in Table 2. Sources and sinks, determined by standard procedures (Chapra, 1997; Thomann, 1987) for the dissolved oxygen indicated less biological activities at station S1 and more at S2. This clearly showed that contamination of river water is very less at upstream of the effluent disposal point and more on the downstream side (Dogan et al., 2009). CONCLUSIONSteady state 2-D analytical solutions are presented for predicting DO distributing in mixing zones of Malaprabha River, by the use of the concept of product law. The effect of background conditions and the effect of sewage quality on DO deficit are considered separately for convenience (Gowda, 1984). For each Station, the deficit due to initial DO level and those due to CBOD and NOD are independently considered. The product law is applied for the analytical solution for the DO deficit prediction.For the Malaprabha River, for which this analytical modeling approach is adopted, the observed and predicted vales of DO deficit are presented in Fig.4. From observations, it can be seen that both the values of the DO deficits agree on the upstream and downstream stretch of the sewage disposal except between the transect 5 and 10, which may be attributed to local conditions (Araujo and Garutti, 2005). Hence, this simulation with respect to DO deficit, confirms the application of this method for predicting the DO deficit values for any changes in the sewage characteristics and stream water, and also for the changes in hydraulic characteristics of both, sewage and river water. The analytical approach can be further extended for other parameters of the river water, such as pH-value, temperature, nitrates, phosphates etc., with modifications in the above equation. ACKNOWLEDGEMENTAuthors thank the Executives and the Dean of Caledonian College of Engineering and the Principals of Gogte Institute of Technology, Belgaum and Bapuji Institute of Engineering & Technology, Davanagere, Karnataka, India for their encouragement and guidance. REFERENCES
Copyright 2010 - International Journal of Environmental Research The following images related to this document are available:Photo images[er10016t1.jpg] [er10016f2.jpg] [er10016t3.jpg] [er10016t2.jpg] [er10016f1.jpg] [er10016f4.jpg] [er10016f3.jpg] |
|