C. Polprasert, M.G. Dissanayake, N.C. Thanh
Waste stabilization ponds (WSP) are becoming pop-ular for treating wastewater, particulary in tropical and subtropical regions where there is an abundance of sun-light, and the ambient temperature is normally high.The ability of WSP system to reduce the biochemical oxygen demand(BOD) of wastewater is well established. Mathematical models have been developed to describe the kinetis of organic degradation in these ponds. How-ever, equally. Important is the effectiveness of WSP systems in reducing pathogenic microorganisms. Because of a lack of sound desing criteria, there are still some doubts as to whether WSP can meet present effluent standards set by many authorities without disinfection.
Kinetic models of bacteria die-off can be expanded to include the effects of algae concentration and organic
loading
Development of mathematical model for bacterial die-off in WSP
Experimental investigation
Result and discussion
Verification of model responses and application
Summary and conclusions
Appendix
Acknowledgments
References
Based on first order kinetics and assuming completely mixed conditions, Marais and Shaw9 proposed a model for the die - off of indicator bacteria in WSP. Because
temperature was found to affect the bacterial removal efficiency substantially, Marais6 altered the model and derived a first-order equation in which the first-order rate constant was assumed to be temperature dependent. Other coliform decay models in WSP, developed by Klock,10 Bowles et al.,11 and Ferrara and Harleman,12 were of the first-order reactions in which the decay rate is temperature dependent.
In fact, the WSP should be considered as a complex system encompassing the existence of several living species, especially the interrelationship of algae and bacteria, which bring about an ecological pattern different from pure culture behavior. Numerous authors11,13 have pointed to a need to improve existing models of coliform decay. The comprehensive model should include the relationship of coliform die-off to other major parameters: algal biomass concentration (Cs), temperature(T), organic loading (OL), sunlight intensity (I), sunlight duration (L), hydraulic detention time (0), substrate degradation rate (Ks), and pond dispersion (d). A research program was undertaken to develop mathematical relationships of the bacterial die-off in WSP14 incorporating two proposed models, one for the bacterial die-off rate coefficient, k, and the other for the algal concentration, Cs. Verification of the results obtained was made with experimental data from the full-scale WSP and some published data for existing ponds in northeast Brazil.15,16 The portion of the work pertaining to mathematical relationships of bacterial die-off (the k model) and data verifications is reported herein.
Ne = No
(1 + kT0)n

where Ne and No are effluent and influent fecal coliform numbers per 100 ml wastewater, respectively; q is in days; and kT, the first-order rate constant for fecal coliform removal at T°C, in days-1, has an approximate value in the temperature range 2 to 212°C, or:
2.6(1.19)T-20 (2)kT =
A later study by Mara et al.16 found the above equations valid for temperatures up to 30°C. Although q might appear to affect the bacterial die-off directly, it actually induces changes in the pond environment, such as the variation of d, Cs, pH, and nutrients that consequently influence the bacterial die-off process. To improve existing models, Thirumurthi17 considered the non-ideal flow in WSP for the model developed for BOD reduction, which included the dispersion number (d) to evaluate the liquid flow and mixing characteristics. The value of d incorporates physical flow characteristics, such as the pond shape, presence of inert zones, flow velocity, and mixing conditions such as wind curretns, thermal stratification, and turbulence.
For ideal plug flow conditions, d is zero, whereas for ideal completely mixed conditions, d reaches infinity. As the value of d becomes greater, the flow deviates from the plug flow towards the completely mixed conditions.
In view of the above mentioned factors, the bacterial die-off in WSP could be related as follows:
Die-off = f(k,d,0) (3)
where
k = bacterial die-off rate coefficient, day-1.
Since the movement of bacteria with seepage water to the outside pond environment could be regarded as minimal when compared to the amount present in the pond,18 the downward diffusion of coliform bacteria was disregarded in the model development. Under steady state conditions, Levenspiels19 flow model for chemical reaction and dispersion was modified in this study to account for the bacterial die-off and dispersion. A steady flow WSP of length L0 through which fluid is flowing with a constant velocity, and in which bacteria is mixed axially with a dispersion coefficient, d was considered as shown in Figure 1. By referring to the elementary section of the pond, the bacterial density mass balance was written as:
input = output + disappearance by pond action
+ accumulation (4)
Equation 4 was modified as in Equation 5:
(output input)bulk flow + (output input)axial dispersion
+ disappearance by pond action
+ accumulation = 0 (5)
in which the final elementary equation is shown in Equation 6.
d . d ²N - d N - k0N = 0 (6)
d Z² d Z
Wehner and Wilhelm20 solved this type of equation for any kind of entrance and exist conditions. The solution of Equation 6 is given below:
Ne = 4ae(1/2d) (7)
No (1 + a)²e(a/2d) (1 a)²e(-a/2d)
where
a = V 1 + 4k0d;
and Ne, No, k, and
d are as defined previously.
Equation 7 can be valid for WSP in which reactions are occurring uniformly throughout the
pond depth at a rate coefficient k. However, in field conditions, owing to uneven
light penetration and possible stratification, the k value can be slightly
different at various depths. In this study the evaluation of the k value was made
for the overall pond depth to compensate for uneven bacterial die-off rates.
The proposed k model. Although k is dependent on tmeperature (Equation 2), other investigators observed k to vary with pH, dissolved oxygen (DO), and nutrient content in the pond.2,4,8 Because pH and DO vary as the algal concentration changes2 and mixed algal species enhance greater bacterial die-off rates,5 the dependence of k on the external factors can be related as:
(8)k = f(T,Cs,OL)
Because the influence of T, OL, and Cs on k are complex, the k model was postulated to be non-linear as follows:
ek = R1w1Tw2CsW3OL (9)
in which R1, w1, w2, and w3 are constants . In order to test the potential lethal effect of sunlight on bacterial die-off in WSP, as reported by Moeller and Calkins,21 the k model could be presented in an alternative form as in Equations 10 and 11.
k = f(T, Cs, OL, UVI) (10)

ek = Ra1wa1Twa2Cswa2wa3OLwa4UVI (11)
where UVI = ultraviolet rays index, dimensionaless, ad Ra1, wa1, wa2, wa3, and wa4 are constants. The UVI index was weighed in accordance with the sunlight intensity (I, cal/day/cm²).
UVI = 1 + I (12)
100
In the absence of sunlight, UVI was
considered to be unity. The constants in Equations 9 and 11 were deter-mined by entering a
series of experimental values obtained from laboratory and pilot scale WSP for different
combinations of T, OL, Cs, and UVI into the multiple regression
equations. The accuracy of these parameters depends upon arriving at a low standard error,
a high multiple correlation, and a high F-value. The total coliform species was
chosen in the experiments as the indicator bacteria for the k-model, but it can be
extended to predict the k-values of other enteric bacterial species.
To accommodate for any species of bacteria, a species constant, l , was included in the
model (Equation 9).
EkF = R1l w1Tw2Csw3OL (13)
Because total coliforms were chosen as the indicator bacteria in this model, the species constant, l , is unity. A series of kF values for fecal coliforms were obtained experimentally for a combination of T, OL, and Cs values. By entering these values into Equation 13, the average R1 value for fecal coliforms (Rf) could be determined, in

Parameter
Range
Temperature, °C
20.0 33.0
Light intensity, cal/cm² d
80.0 3 30.0
Light duration, h/d
8.0 14.0
Influent COD, mg/l
150.0 170.0
Dispersion number (d)
0.100 0.200
Detention time, days
3.0 10.0
______________________________________________________________________________
which:
RF = l 1R1 (14)
where
RF = R value for fecal coliforms,
R1 = R value for total coliforms, and
l 1 = Species constant for fecal coliforms.
Laboratory-scale WSP (LWSP). Five rectangular LWSP units, made of 6-mm thick PVC sheet, each with different dimensions, were used in the experiments (Figure 2). The LWSP inlet was connected to a flow inducer to obtain a constant influent flow, and two baffles were fixed ner the pond inlet and outlet to minimize short-circuiting and free-flow conditions. A 500-1 polyethylene vessel with a stirrer was used as the feed tank to which feedlines and pumps were connected to facilitate continuous operation of the system. The experiments were conducted in a controlled temperature room and pond illumination was accomplished by providing a set of fluorescent bulbs fitted to a wooden stand. To overcome the red and orange wavelenght deficiency in fluorescent bulbs, a number
Table 2 Ranges of operating conditions for pilot-scale WSPPWSP 2 &
Temperature, °C
28.6 32.6
28.6 32.6
Light intensity, cal/cm² . d
321.6 643.5
321.6 643.5
Light duration, h/d
11.5 12.5
11.5 12.5
Influent COD, mg/l
160.0 180.0
40.0 95.0
Dispersion number (d)
0.115 0.195
0.125 0.215
Detention time, days
5.8 17.9
3.6 11.1
Dimension, I X w X d, m3
4.5 X 20 X 0.86
4.0 X 2.0 X 06

Figure 4 -- Layout of AIT waste stabilization ponds (FWSP).
Pilot-scale WSP (PWSP). The PWSP exeperiments were carried out under ambient conditions, with three ponds connected in series (Figure 3). A constant head tank was used to feed wastewater to PWSP at a constant rate, but a pump with a timer was used for sewage feeding periodically to avoid pipe blockage. The ranges of operating conditions of the PWSP experiments are shown in Table 2. The same measurements taken for the LWSP, plus the diurnal variations of pH and DO, were taken for each PWSP.
Full-scale WSP (FWSP). The wastewater treatment facility (Figure 4) at AIT comprises two parallel sets of ponds in series (FWSP). The FWSP1 serves a facultative function, while the FWSP2 is a maturation pond. The final effluent is discharged into a nearby canal. The EWSP operation is under the control of the AIT physical plant, but cooperation was provided in the measurements of flows and other important parameters similar to the LWSP and PWSP.
Tracer study. To determine the dispersion characteristics of the LWSP2, PWSP, and FWSP, tracer studies were carried out for different q of each pond. Sodium chloride (NaCl) was the impulse tracer material for all the experiments, as in the experiments of Thirumurthi17. According to Levenspiel19, the response tracer concentration should be monitored at the exit stream at fixed time intervals. The amount of input impulse tracer concentration varied from 30 g in the LWSP to about 700 kg in the FWSP. The base amount of NaCl present in the wastewater was taken into consideration when the exit responses were monitored. The calculation of the d values was made with theroposed by Levenspiel and Smith, which is described below:
Mean detention time (actual), q = S q
i Ci (15)
S Ci
Standard deviation, s ² = S Ciq i - (q )² (16)
S Ci
if Æ = q i
q
then s ² = 2d 2d² (1-e-1/d) (17)
The term, d, can be calculated by trial and error from Equation 17 in which.
q i = time after impulse injection, days; and
Ci = tracer response concentration at the exist stream, mg/l.
Though tracer studies require careful
attention and much time, they are one of the most reliable techniques for evaluating the
actual flow pattern in a WSP. Previous investigators, such as Thirumurthi17,
Bowles et al., and Uhlmann, based their mathematical formulations of WSP
performance on partial mixing conditions and flow pattern.
Analytical methods. All analyses were undertaken according to the methods described
in "Standard Methods"24. DO measurements were conducted with a DO
meter. PH of the water samples was measured with a fluent flows of the LWSP. Intensity of
artificial light and sunlight was measured with a bimetallic actinograph instrument. The
samples for algal analysis were initially filtered through a 150-m m mesh to separate out
zooplanktons, measured for the percent absorbance with a spectrophotometer at a wavelength
of 650 mm, and weighed for dry algal concentration with the specimen curve prepared in
advance.
The total coliform counts were enumerated by the standard MPN method using the five-tube dilution technique. Appropriate dilutions, prepared in buffered dilution water, were inoculated into lactose broth, presumptively tested by incubation at 35°C for 24 hours. All the positive tubes were the confirmatively tested by subculturing into brilliant green lactose bile broth at 35°C for 48 hours. The fecal coliform MPN counts were determined by subculturing all the positive presumptive tubes into an EC medium at 44.5°C for 24 hours.
Performance of LWSP. The data of LWSP performance, presented in Table 4, show the variation in percent reduction of total coliforms (PAT) and fecal coliforms (PAF) ranging from 73 to 96 and 78 to 97, respectively, with changes in the magnitude of other environmental factors. The COD removal as reported by Dissanayake14 was found to range between 40 and 80%. Because the LWSP were operated as single-cell ponds, it is unlikely they could achieve better results than those reported above. Both Marais6 and Parker8 noted the superior performance of a series of ponds as compared with a single cell pond. The concentrations of Cs in the ponds varied from 28 to 260 mg/l, depending on the values of d, 0, OL, T, L, and I. No attempt was made to investigate
Table 3 AIT campus wastewater characteristics and performance of the FWSP ____________________________________________________________________________________________ Raw wastewaterEffluent of
Effluent of
Parameter
Range
Mean
FWSP 1
FWSP 2
____________________________________________________________________________________________
COD, mg/l (filtered)
100-440
180
50
30
BOD, mg/l (filtered)
60-150
120
35
15
SS, mg/l
20-120
110
80
115
VSS, mg/l
10-90
45
--
--
Org-N, mg/l
2-5
--
--
--
NH3-N, mg/l
4-10
5.5
3.0
Trace
P04-3, mg/l
0.5-1.5
--
0.4
0.3
PH
6.2-7.5
7.3
7.5
8.0
Temperature, °C
28-33
_
28-33
2 8-33
Total coliforms, MPN/100 ml 22 x 105-16 x 106
11x10 6
15 X 105
24 X 103
Fecal coliforms, MPN/100 ml 16 x 103-92 x 104
46 x 10 3
93 X 102
11 X 102
Detention time (0), days
&a mp;a mp;n
bsp;
8
20
Dimension, I X w X d, m3
&a mp;a mp;n bsp; 50 X 23 X 2.4
120 X 43 X1.3
____________________________________________________________________________
________________________________________________________________________________________________________________________
Percentage
reduction
_____________
Total
Fecal
Dispersion
Detention
Algal
Organic
Light Light
coliform coliform
number
time
concentra
loading Temperature
duration intensity
PAT
PAF
d
0
(mg/1) (kg COD
(°C)
(h/d) (cal/cm2.d)
&a mp;a mp;n
bsp; C
OL
T
L
I
______________________________________________________________________________________________________________
74.71
79.41
0.160
3.0 215.0
118.0
22.5
14.0 326.1
87.71 89.54
0.110
3.0 140.2
105.6
30.0
14.0
81.7
86.14 88.00
0.160
3.0 129.4
101.7
30.0
14.0
81.7
77.72 81.76
0.170
3.0 31.8
113.4
30.0
14.0
81.7
76.17 80.00
0.190
3.0
31.8 106.8
30.0
14.0 81.7
77.44 80.00
0.200
3.0 28.0
52.7
30.0
14.0 81.7
73.00 79.41
0.165
4.0 182.0
89.1
22.5
14.0 326.1
82.60 96.07
0.120
6.0 89.5
54.7
22.5
14.0 326.1
93.51 97.39
0.120
6.0 158.6
58.9
33.0
14.0
326.1
96.16 97.39
0.120
6.0 136.0
57.6
33.0
16.0
372.7
77.50 80.93
0.120
6.0 260.0
57.0
22.5
12.0 279.5
77.50 78.12
0.120
6.0
50.0 57.0
22.5
10.0
232.9
78.80 80.00
0.120
6.0 150.0
55.9
22.5
14.0
326.1
73.71 79.41
0.120
6.0
62.0 55.9
22.5
8.0
186.3
83.07 90.00
0.110
8.0 177.0
42.5
22.5
12.0
279.5
86.15 91.53
0.100
10.0 148.0
34.0
22.5
14.0
326.1
____________________________________________________________________________________________________________________________
Percentage |
|||||||||
Total |
Fecal |
Dispersion |
Detention |
Algal |
Organic |
Temperature |
Light |
Light |
Remarks |
| 84.09 80.00 77.14 90.00 94.25 85.86 95.42 94.25 94.13 97.40 94.55 77.21 96.19 97.31 99.00 98.24 97.57 97.39 99.78 98.82 95.50 98.67 95.72 87.14 99.00 95.58 90.00 99.86 98.18 98.77 98.54 96.00 85.00 98.96 95.47 97.20 92.00 97.50 95.57 86.15 98.77 95.00 97.60 92.27 83.52 99.45 97.38 99.41 |
87.85 84.11 80.00 90.00 96.19 89.42 96.82 94.81 93.92 97.54 90.00 81.85 97.60 98.71 93.33 99.23 96.15 96.00 98.70 99.09 90.00 98.62 95.72 90.00 99.68 97.64 87.50 99.85 99.00 96.53 99.70 97.03 85.00 98.25 91.22 98.60 93.92 98.98 95.00 89.23 99.00 93.57 98.35 98.91 86.00 97.17 97.21 99.09 |
0.195 0.215 0.215 0.175 0.190 0.190 0.145 0.185 0.185 0.130 0.175 0.175 0.130 0.175 0.175 0.125 0.155 0.155 0.125 0.155 0.155 0.125 0.155 0.155 0.125 0.155 0.155 0.125 0.155 0.155 0.125 0.155 0.155 0.125 0.150 0.150 0.145 0.180 0.180 0.125 0.150 0.150 0.125 0.155 0.155 0.115 0.125 0.125 |
5.8 3.6 3.6 7.0 4.3 4.3 8.9 5.6 5.6 9.4 5.8 5.8 9.9 6.1 6.1 11.5 7.1 7.1 11.7 7.2 7.2 12.4 7.7 7.7 12.3 7.7 7.7 12.8 7.9 7.9 13.4 8.3 8.3 13.4 8.3 8.3 13.4 8.3 8.3 14.1 8.8 8.8 14.9 7.2 7.2 17.9 11.1 11.1 |
120.0 128.0 113. 121.0 210.0 160.0 36.5 172.0 162.0 154.0 178.0 85.0 108.0 286.0 166.0 178.0 220.0 122.0 230.0 375.0 290.0 180.0 200.0 80.0 151.0 210.0 140.0 270.0 265.0 260.0 163.0 107.0 80.0 150.0 140.0 170.0 155.0 153.5 130.0 104.0 275.0 160.0 110.0 157.0 70.0 36.5 70.0 162.0 |
263.3 70.6 70.6 225.7 94.3 98.3 161.3 90.9 72.0 157.9 61.0 80.3 152.9 63.4 87.0 132.5 30.4 65.6 130.3 69.7 84.9 19.2 40.1 52.2 123.2 48.0 54.9 121.6 44.8 58.3 110.0 39.5 36.4 113.3 44.8 39.7 161.5 70.4 01.7 103.5 60.2 76.9 120.8 51.6 45.7 85.5 20.8 29.8 |
29.2 29.2 29.2 28.6 28.6 28.6 30.9 30.9 30.9 29.3 29.3 29.3 31.4 31.4 31.4 30.3 30.3 30.3 29.1 29.1 29.1 30.7 30.7 30.7 30.2 30.2 30.2 32.2 32.2 32.2 29.0 29.0 29.0 30.5 30.5 30.5 31.2 31.2 31.2 31.6 31.6 31.6 29.5 29.5 29.5 32.6 32.6 32.6 |
11.5 11.5 11.5 11.5 11.5 11.5 12.0 12.0 12.0 12.5 12.5 12.5 12.5 12.5 12.5 11.5 11.5 11.5 11.5 11.5 11.5 12.0 12.0 12.0 12.0 12.0 12.0 12.5 12.5 12.5 12.0 12.0 12.0 12.0 12.0 12.0 12.0 12.0 12.0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 |
459.7 459.7 459.7 474.9 474.9 474.9 518.4 518.4 518.4 475.0 475.0 475.0 582.2 582.2 582.2 509.4 509.4 509.4 526.0 526.0 526.0 566.9 566.9 566.9 643.5 643.5 643.5 643.5 643.5 643.5 462.2 462.2 462.2 505.6 505.6 505.6 459.7 459.7 459.7 516.9 516.9 516.9 321.6 321.6 321.6 505.6 505.6 505.6 |
PWSP1 PWSP2 PWSP3 PWSP1 PWSP2 PWSP3 PWSP1 PWSP2 PWSP3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 PWSP 1 PWSP 2 PWSP 3 |


level), and Cs in the pond water. In this
study, the DO at dawn did not reach zero because the applied loadings of the PWSP were les
than 280 kg COD/ha d, the level considered as facultative.3,29 The
decrease of Cs in the later period of a day could be caused by light and
nutrient limitations, competition for food, algal cell sedimen-tation, and unfavourable
conditions existing at high pH conditions.
Performanceof FWSP. The data in Table 3 indicated a considerable reduction of COD
and BOD int the FWSP operated in series as their concentrations in the FWSP2 effluent were
only about 30 and 15 mg/l, respectively. However, the SS concentration of 155 mg/l was
con-sidered high, which was caused by algal cell suspension in the effluent. There were
decreases in the total and fecal coliform densities of about 80% in the FWSP1 (0 = 8 days),
and about 90 to 98% in the FWSP2 (0 = 20 DAYS). The removal of organics and
bacteria in the FWSP are comparable with other ponds reported in the literature.29

It seems that another polishing pond, connected in series to the FWSP2, should be provided to achieve the stringent standards of SS and bacteria required by some regulatory agencies.26
It can be seen from the results of the
LWSP, PWSP, and FWSP that bacterial die-off is a complex phenomenon involving various
factors and interactions in the WSP. Because in actuality, the WSP are not in either
completely mixed or plug-flow conditions, and not just a single parameter, such as T,
controls the bacterial die-off rate, the proposed Equations 7 and 9 or 11, as models for
bacterial die-off, would better represent the pond performance.
The k model. The experimental k values were determined by entering
each data ser of PAT, d, and q in Tables 4 and 5 as input variables into Equation
7. Then, together with the respective values of T, OL, Cs, and UVI, the
experimental values as determined were entered into the multiple regression equations
(Equations 9 and 11) to evaluate the regresion coefficients (or constants) using a
computer program form a scientific subroutine package. The multiple regression results
including the multiple correlation coefficients, standard errors of estimate, and
F-statistic values obtained for Equations9 and 11 were, respectively, 0.695 and 0.719,
0.152 and 0.294, and 19.039 and 16.109. Although Equation 9 had a slightly lower multiple
correlation coefficient, its standard error of estimate and F-statistic value were lower
and higher, respectively, than Equation 11, indicating that both equations represented a
certain degree of correlation withing the same range btween k and the parameters T,
Cs, OL, and I. In a WSP, k would depend on factors
other than just those above mentioned parameters. But for practical reasons, Equation 9
was chosen to represent the k model, which can be written as in Equation 18:
ek = 0.6351(1.0281)T(1.0016)Cs(0.9994)OL (18)
The empirical relationship obtained for k in Equation 18 is essentially an extension of Equation 2, which quantifies the dependence of k on temperature only, and it eliminates the necessity for a temperature correction. However, a comparison of the k values based on Equations 18 and 2 is not relevant because the equations governing bacterial die-off (Equations 7 and 1, respectively) are of different order. Although Equation 18 does not include the parameter I or UVI, the direct relationship between I and Cs has been substantially reported,2 thus the term Cs indirectly represents the influence of I on k.
The regression coefficients represent the
factor by which each parameter (T, Cs, and OL) contributes to the
independent parameter determining whether the contribution is negative or positive. For
Equation 18, an increment in either T or Cs results in an increase in the k
vaue while the opposite is true for OL, these phenomena could be explained as
follws. Although Marais6 observed that beyond a temperature of 21°C, k
decreased as temperature increased, a later study by Mara and Silva16 showed
that k increased with temperature up to at least 30°C. At a temperature range of
21 to 30°C, bacterial activity is stimulated and results in a higher probability of
bacterial survival, but other factors, such as increased algal activity and pH, cause
higher die-off rates by offsetting the stimulated bacterial activity. The algal
concentration in the pond was found to increase with the onset of the light period. With
the consequent increae in pH (Oswald and Gotaas30 and Figures 5 and 7), the pH
variation from about 8 to 10 canyield a hostile environment for bacterial survival4.
The regression coefficient for OL indicates that an increase in organic loading
causes a slight decrease in k because of more nutrient availability for the
bacterial activity.28 Therefore, the regression coefficients and their signs,
as shown in Equation 18, are considered to be in accordance with actual pond performance
and the literature, when operated within the ranges employed for the LWSP and PWSP in this
study.
Species coefficient for fecal coliforms (l 1). The fecal coliform
die-off coefficients (kF) were calculated by entering the
respective values of PAF, d, and 0 in Tables 4 and 5 into Equation 7. The
values of R1l 1 or RF could be determined by entering the values of
the calculated kF, T, Cs, and OL into Equation 13, in which the
values of w1, w2, and w3 were the same as in Equations 9
or 18. By entering the mean RF value into Equation 14, the species constant, l 1
was found to be 1.1274, and the kF is given as follows:
ekF = 1.1274(0.6351)(1,0281)T(1.0016)Cs(0.994)OL (19)
Thoug the k to kF ratio is constant, the resulting PAT to PAF ratio is variable because Equations 18 and 19 have to be used with Equation 7, which has a higher order to predict the bacterial survival ratio in WSP. From Equations 18 and 19, the value of kF is marginally higher than k, which is in accordance with the results of Fransmathes, who reported better survival of the total coliforms than the fecal coliforms in WSP.
Comparison with predicted results |
Comparasion with |
||||
CC |
MPR |
SDPR |
MPR |
SDPR |
|
Total coliform percentage reduction (PAT) based on 22 FWSP experimental dataa Fecal coliform percentage reduction (PAF), based on 22 FWSP experimental dataa Fecal coliform percentage reduction (PAF), based on 14 experimental data on ponds in NE Brazil |
0.986b 0.988c 0.994c |
0.302 0.198 0.352 |
1.090 0.675 1.061 |
-2.730 -2.451 -2.705 |
3.934 3.302 2.961 |
Note CC = correction coefficient.
MPR = mean of percentage errors.
SDPR = standard deviation of percentage errors (excluding errors over 10%)
a = the ranges of I, L and d were similar to those reported in Table 2.
b = the values of k were determined from Equation 18.
c = the values of kF were deterined from Equation 19.
Comparing the 22 experimental data observed PAT in the FWSP with those predicted by Equations 7 and 18 gave the correlation coefficient of 0.986, indicating that the models developed in this study performed with a high degree of accuracy in the prediction of the total coliform die-off. The mean and standard deviation of percentage errors (MPR and SDPR, respectively) as resulted form the PAT prediction by Equation 7, were lower than that of Equation 7 yielded a correlation coefficient of 0.988, and there were lower mean and standard deviation of percentage errors in Equation 7 than that of Equation 1. By assuming the values of I, L, and d to be 390 cal/cm² d, 11.5 h/d, and 0.150, respectively, for the pond data in northeast Brazil, the comparison between the fecal coliform percentage reduction and those predicted by Equations 7 and 1 could be made. The observed PAF and those predicted by Equation 7 had a correlation coefficient of 0.994, further asserting te accuracy of the proposed models. There were lower mean and standard deviation of percentage errors of the predicted PAF by Equation 7 than those by Equation 1. The main reasons for this difference are probably because Equation 1 assumes complete-mixing conditions in the pond and it includes only T and 0, whereas the proposed models cater for actual dispersion characteristics of the pond and encompass other important parameters, such as Cs and OL.
Because the values of the multiple regression coefficients in Equation 18 appear to be close to unity, it could be interpreted that the model was not sensitive. An attempt was made to observe the effects that variation of the T, Cs, and OL have upon the value of k. It was found that as T increases from 10 to 30°C, k changes by as much as 0.55; as Cs increases from 25 to 300 mg/l/k changes by as much as 0.44; and when OL increases from 50 to 300 kg. COD/had, k decreases by as much as 0.15. For the parameter ranges mentioned above, k changes from a minimun value of 0.17 to a maximum of 0.68. Equation 7 was tested to observe the effects that variation in k and d have on the coliform survival ratio. By assuming the following: T = 25°C, OL = 100 kg COD/ha d, Cs = 250 mg/l, 0 = 10 days, and d = 0.4; it was found that if k changes from 0.5 to 1.0, the Ne/No ratio decreases from 0.061 to 0.012, a decrease in coliform density of about 80%. However, with the same values of T, OL, Cs, 0, and the k value of 0.5, it was observed that as d changes from 0.4 to 0.8, the Ne/No ratio increases from 0.061 to 0.089, an increase in coliform density of about 46%. Hence it is apparent that the changes in the values of T, Cs, and OL can cause significant changes in the k value, which in turn, together with the d value, can cause significant variations in effluent coliform quality.
The Wehner and Wilhelm equation (Equation 8) may seem complicated, but charts may be prepared similar to that of Thirumuthi17 (for BOD reduction) to read off directly the coliform survival ratio as a function of the pond mixing pattern (d) and the value of k0. Because the second term in the denominator in Equation 7 is small, it may be neglected, in which case the equation can be simplified as:
Ne = 4ae(1-a/2d)
(20)
No (1 + a)²
The error of Equation 20 may be significant when the value of d exceeds 2.0. However, based on the data in Tables 4 and 5 of this study and of Nashashibi,32 d seldom exceeds 1.0 because of low hydraulic loads. In practice, the values of T, OL, Cs, and d need to be known if the coliform survival ratio in WSP is to be predicted from Equations 7 and 18 or 19. T can be measured onsite directly and OL can be calculated or taken from the literature,25,29 thus only Cs and d need to be determined. Oswald and Gotaas,30 reported the rate of algal yield to vary from 2.5 tons/ha moth in the winter to 12.5 tons/ha month in the summer. From continuous cultures, they found the dry weight algal cell material to be a logarithmic function of the substrate (BOD), up to approximately 400 mg/l (for the wastewater BOD ranging from 75 to 375 mg/l). The corresponding algal cell concentrations were from 125 to 325 mg/l, respectively. A model was recently developed14 that relates the Cs in WSP to the substrate degradation rate, hydraulic detention time, and other factors, such as T, L¸and I. It was found to be satisfactory in predicting the algal concentration. The value of d can be obtained by conducting tracer studies in existing, similarly shaped ponds. It should be noted that the proposed equations
(Equations 7, 18, and 19) are applicable
to the types of single-cell ponds experimented with (facultative and maturation), and the
ranges of operating conditions conducted in this study. Additional research is required to
modify these equations so that they can be used to predict bacterial die-off in a series
of WSP.
Design example. Determine the hydraulic detention time, 0, of WSP treating
domestic wastewater so that the final effluent will contain less than 100 fecal
coliforms/100 ml, which is the standard for reuse in restricted irrigation. The initial
COD and fecal coliform concentrations are 300 mg/l and 10/100 ml, respectively. The
expected monthly temperature in the pond in tropical areas is about 20°C in winter.
Tracer studies and water analyses in the existing, similarly shaped pond indicate that the
values of d and Cs are approximately 0.200 and 200 mg/l, respectively.
Solution 1. Select the OL value to be 200 kg. COD/ha d. The value of kF
can be determined according to Equation 19:
ekF = 1,1274(0.6435)(1.0281)20(1.0016)200(0.9994)200
in which kF was found to
be equal to 0.433 day -1. Take a = Ö 1+4kF0 d and put it
into Equation 20. By trial and error, the WSP hydraulic detention time (0) is found
to be 34 days to produce effluent with fecal coliform concentration less than 100/100 ml.
Solution 2. If Equations 1 and 2 are to be used in the design, take from Equation
2, kT = 2.6(1.19)20-20 = 2.6 day 1.
From Equation 1, the hydraulic detention time, 0, can be determined as follows:
if n = 1, 0 = 384 days,
n = 2, 0 = 12 days, and
n = 3, 0 = 3.5 days.
It is difficult to say, based on the above results, which method is superior in terms of WSP design for bacterial reduction. Although Solution 2 gave the anomalous value of 0 for the case n = 1, for the case n = 1, for the case n = 2, the total 0 is 12 x 2 = 24 dys, which is less than the Solution 1 result. At this stage, it could be generally stated that the proposed equations (Equations 7, 18, and 19) should be used in the design of single-cell facultative or maturation ponds, and that bacterial die-off in the ponds can be predicted by these equations with high accuracy. Equations 1 and 2 should be used in the design of WSP in series, as they gave more choices in selecting the n-series of ponds.
(1) A multiple linear regression equation was developed to relate the bacterial die-off rate coefficient (k) to other parameters, such as T, Cs, and OL. This equation can be extended to predict the k values of other bacterial species by the introduction of the species coefficients into the regression coefficient.
(2) The non-ideal flow equation of Wehner and Wilhelm was proposed for the prediction of the bacterial survival ratio in WSP. Besides the k value, this equation accounts for the pond dispersion factor, d, and the actual hydraylic detention time, 0.
(3) The Wehner and Wilhelm equation was found to yield better results in predicting bacterial die-off in WSP than the first-order rate equation. This is probably because the former includes other important parameters influencing pond performance and bacterial die-off. Verification of the model responses was made with data of the full scale WSP and some ponds in northeast Brazil.
(4) Sensitivity analyses were carried out to show the effects that changes in magnitude of these influencial parameters had upon the k value and the bacterial survival ratio.
Cs = algal concentration, mg/l
d = dispersion number, dimensionaless
I = light intensity, calories/cm² d
k = bacterial die-off rate coefficient, day-1
Kf = fecal coliform die-off rate coefficient, day -1
kT = first-order rate constant for bacterial die- off, day -1
L = light duration, h/d
n = number of waste stabilization ponds in series
Ne,No = effluent and influent bacterial densities, respectively, MPN/100 ml
OL = influent COD loading rate, kg COD / ha d
PAF = percentage fecal coliform reduction
PAT = percentage total coliform reduction
R1,w1 = multiple regression coefficients for
w2,w3 the k-model
Rai,wa1, = multiple regression coefficient for the
wa2, wa3 k-model with UVI as a parameter
wa4
RF = a multiple regression coefficient for the k-model of fecal coliforms
T = pond wter temperature, °C
UVI = ultra-violet ray index, dimensionaless
0 = pond hydraulic detention time, days
l = bacterial species constant, dimensionaless
l 1 = species constant for fecal coliforms, dimensionaless
Hermann, E.R., and Gloyna, E.F., "Waste Stabilization Ponds Experimental Investigation." Sew, and Ind. Wastes, 30, 511 (1958).
Oswald, W.J., "Light Conversion Efficiency of Algae Grown in Sewage". J. San. Eng. Div., Amer. Soc. Civil Eng. 86, SA4, 71 (1960).
Porges, R., and Mackenthun, K.M. "Waste Stabilization Ponds, Use, Function and Biota". Biotech. Bioeng., 5, 255 (1963).
Parhad, N.M., and Rao, N.U., "Effect of pH on survival of escherichia coli". J. Water Pollut. Control Fed., 46, 980 (1974).
Davis, E. M., and Gloyna, E.F., "Bacterial Dieoff in Ponds". J. San. Eng. Div., Amer. Soc. Civil Eng., 98, SA1, 59 (1972).
Marais, G. V. R., "Faecal Bacterial Kinetics in Stabilization Ponds". J. Env. Eng. Div., amer. Soc. civil Eng., 100, EE1, 119 (1974).
Mancini, J.L., "Numerical estimates of coliform mortality rates under various conditions". J. Water Pollut, Control Fed., 50, 2477 (1978).
Parker, C.D., "Microbiological aspects of lagoon treatment". J. Water Pollut. Control Fed. 34, 149(1962).
Marais, G. v. R. , and Shaw, V.A., "A Rational Theory for the Design of Sewage Stabilization Ponds in Central and South Africa". Trans. S. Afr. Inst. Civil Eng., 13, 205 (1961).
Klock, J.W., "Survival of coliform bacteria in wastewater treatment lagoons". J. Water Pollut. Control Fed., 43, 2071(1971).
Bowles, D.S., et al., "Coliform decay rates in waste stabilization ponds". J. Water Pollut. Control Fed., 51, 87 (1979).
Ferrara, R..A., and Harleman, D.R.F., "Dynamic Nutrient Cycle Model for Waste Stabilization Ponds". J. Env. Eng. Div., Amer. Soc. Civil Eng., 106, EEI, 37 (1980).
Finney, B.A., and Middlebrooks, E.J., "Facultative waste stabilization pond desing". J. Water Pollut. Control Fed., 52, 134 (1980).
Dissanayake, M.G., "Kinetics of Bacterial Die-off in Waste Stabilization Ponds". Doctoral dissertation. No.EV-80-2, Asian Institute of Technology, Bangkok, Thailand (1980).
Mara, D.D. et al., "Design Verification for Tropical Oxidation Ponds". J. Env. Div. Amer. Soc. Civil Eng., 105, EE1, 151 (1979).
Mara, D.D., and Silva, S.A., "Sewage Treatment in Waste Stabilization Ponds: Recent Research in Northeast Brazil". Prog. Water Tech. 11, 341(1979).
Thirumurthi, D., "Design Principles of Waste Stabilization Ponds". J. San. Eng. Div., Amer. Soc. Civil Eng. 95, SA2, 311 (1969).
Krone, R.B., et al., "Movement of Coliform Bacteria Through Porous Media". Sew, and Ind. Wastes, 30, 1(1958).
Levenspiel, O., "Chemical Reaction Engineering". 2nd Ed., John Wiley and Sons, Inc., New York, N.Y. (1972).
Wehner, J.F. and Wilhelm, R.H., "Boundary Conditions of Flow Reactor". Chem. Eng. Sci., 6, 89 (1956).
Moeller, J.R., and Calkins, J., J., "Bactericidal agents in wastewater lagoons and lagoon design". J. Water Pollut. Control Fed., 52, 2442, (1980).
Levenspiel, O., and Smith, W.K., "Notes on the Diffusion-type Model for the Longitudinal Mixing of fluids in flow". Chem. Eng. Sci., 6, 227 (1957).
Uhlmann, D., "BOD Removal Rates of Waste Stabilization Ponds as a Function of Loading, Retention Time, Temperature and Hydraulic Flow Pattern". Water Res., 13, 193 (1979).
"Standard Methods for the Examination of Water and Wastewater". 14th Ed., Amer. Pub. Health Assoc., New York, N.Y. (1976).
"Wastewater Engineering Collection, Treatment Disposal". Metcalf and Eddy, Inc., McGraw - Hill, New York, N.Y. (1972).
World Health Organization, "Reuse of Effluents: Methods of Wastewater Treatment and Health Safegurds". WHO Technical Report Ser. 517, Geneva (1973).
Gann, J.D., et. Al., "Aerobic bacteriology of waste stabilization ponds". J. Water Pollut. Control. Fed. 40, 185, (1986).
McGarry, M.G., and Bouthiliar, P.H. "Survival of S. typhi in Sewage Oxidation Ponds". J. San Eng. Div., Amer. Soc. Civil. Eng. 92, SA4, 33 (1966).
McGarry, M.G., and Pescod, M.B., "Waste Stabilization Pond Criteria for Tropical Asia". Proc. 2nd. Intl. Symp. Waste Treatment Lagoons, Kansas City, Mo. (1970).
Oswald, W.J., and Gotaas, H.B., "Photosynthesis in Sewage Treatment". Trans. Amer. Soc. Civil Eng., 81, 73 (1955).
Fransmathes, J.R. "Bacteria and Lagoons." Water and Sew. Works, 117, 90 (1970).
Nashashibi, O.I., "A Process Approach to the Design of Waste Stabilization Ponds." M.S. thesis, Univ. Of Miami, Coral Gables, Florida (1976).
Comentarios al Webmaster |
[ Homepage CEPIS ] |