JKAU: Mar. Sci., Vol. 18, pp: 149-167 (2007 A.D. /1428 A.H.)

Application of a Mathematical Model for Estimating the Pattern of Pollution Dispersion at the Coast of

Ebtessam E.E. Mohamed, El-Sayed M. El-sayed*, Nabil N. Saad, and Hassan A. Abou Tahoun Marine Physics Lab., National Institute of Oceanography and Fisheries, Alexandria *Department of physics, Ain Shams University, ,

Abstract. During December 1996, 27 sampling stations were explored along the coastal zone in front of Alexandria situated between El- Agamy beach in the West and Abu-Qir headland in the East. Also, current measurements have been obtained by using current meter and Drougs. The dispersion coefficient Kx and Ky of water during the time of measurements in front of Alexandria were calculated. A predictive model based on the equations described by Bonazountas (1987) was applied to estimate initial dilution and wastewater field in front of Alexandria coastal zone.

Introduction The state of pollutants in inshore water of the Mediterranean has reached critical levels. This is mainly due to high quantities of domestic sewage discharged either partially treated or even untreated into the sea through rivers, outlets and pipelines. Alexandria is one of the largest coastal cities on the Mediterranean; the waste disposal of the city is mainly discharged into the through many short outlets. Some of these disposal sites, e.g., the pumping station at Abu Qir and El Max bay, discharge municipal wastewater mixed with drainage water. El Max bay is mainly polluted by El Umum drainage water, which 6 pours about 6x10 (m3/day) of untreated industrial, agricultural and domestic’s effluents. On the other hand, a major outfall at Kayet Bay and minor at different coastal sites extending from the Eastern Harbour to 149 150 Ebtessam E.E. Mohamed et al.,

El Maamora dispose sewage wastes into the beaches. The bulk volume of domestic wastewater, which exceeds 400,000 m3/day, is discharged into the sea through the principal pipe at Kayet Bay outfall. Much intensive attention was paid to study the environmental effect on El-max and Eastern Harbour. The major hydrographic, current measurements, the physico-chemical characteristics and sediment of the area under investigation had been studied by; El-Wakeel and El-Sayed (1978); Abdallah (1979); Farag (1982); Rifaat (1982); Abul-Dahab (1986), Halim (1988, 2004); El-Nagar (1994); Said et al (1994); Farag. et al (2000). This study is carried out to evaluate the effect of sewage and drainage wastewater on the coast of Alexandria particularly in front of Kayet Bay outfall and El Umum drain. Also, to illustrate the final patterns resulting from mixing between the ambient seawater and the discharge pollutant.

Data and Methods of Analysis During Dec. 1996, the coastal zone in front of Alexandria between El Agami and Abu Qir headline was surveyed. Through this period 27 sampling stations were explored along five sections almost normal to the coast. Also, some stations were taken scattered around Kayet Bay outfall and inside the eastern harbour. Current measurements have been obtained by two different ways:

1. Current meter model Anderaa RCM-7 had been used to gather current data at constant depth 3 m under surface at station 1-22 after anchoring the boat in every station. The current meter was adjusted to average and record one reading every five minutes.

2. Detected the paths of wind-induced currents generated in the upper surface layer by using two drogues submerged to about 40 cm under sea surface. A motorboat followed up the movement of the two drogues and their locations were fixed approximately every thirty minutes.

The locations of sampling stations taken during December 1996 are shown in Fig. 1. Application of a Mathematical Model for … 151

Kayet Bay Latitude

El Max Bay

Longitude Fig. 1. The locations of sampling stations taken during December 1996. Current Regime and its Dispersing Ability It is clear from the current regime in the area of investigation (Fig. 2) the predominant direction of flow vectors was to the southwest, except El Max sector, with their speed fluctuating between 1.007 and 7.5 cm/sec. The current was coastward in front of the area between Kayet Bay and Ras El teen and inside the Eastern Harbour.

Longitude Fig. 2. Current regime in the area of study during December 1996, (current meter data).

152 Ebtessam E.E. Mohamed et al.,

In El Max sector there was no predominant direction. The near shore current was flowing seawards, while the offshore one moved to the west and it was almost parallel to the coastline. Generally, the tide is ignored in this study because its components was very small and contribute nothing in the current intensity. The tracks of the drogues are illustrated in Fig. 3. It is clear that the predominant direction of the surface currents was the southwest. The current value calculated from drogues experiments are used in the calculating of diffusion coefficient (Kx and Ky) using the method described by List et al (1990) and can be summarized as follow: Consider the location coordinates (x,y) of a drogue J at a time I are defined by (Xij, Yij).From these coordinates it is possible to develop estimates of the Lagrangian velocity of the drogue and the position of centroid of the drogue set may be computed. The centroid of set N drogues at a time I is given by:

B B

Longitude Fig. 3. Tracks of Drogues experiment in front of Alexandria coast during December 1996. Application of a Mathematical Model for … 153

∑ X ij (1) j X i = N And similar for Y coordinate

∑Yij (2) j Y i = N

The mean velocity of the drogue set is simply the time derivative of the position coordinates. The variance in the drogue position defined by:

2 (3) ∑()− XX iij j σ 2 = xi N −1

2 (4) ∑()− YY iij j σ 2 = yi N −1

Following Okubo (1974), the dispersion of the drogue distribution can be estimated as:

222 i ( += σσσ yixi ) 2/ (5)

And the relative dispersion coefficient K is given by:

1 ∂ 2 1 Δσσ 2 (6) tiK )( = i ≅ i 2 ∂ 2 Δtt

The spatially dependent relative dispersion coefficient is given by:

1 ∂ 2 1 Δσσ 2 (7) tiK )( = xi ≅ xi x 2 ∂ 2 Δtt

2 2 1 ∂ yi 1 Δσσ yi (8) tiK )( = ≅ y 2 ∂ 2 Δtt

154 Ebtessam E.E. Mohamed et al.,

Then

( += KKK yx ) 2/ (9)

The dispersion coefficients of the coastal water in front of Alexandria were found to be:

Kx =45049 cm2/sec Ky = 870609cm2/sec

The obtained value of Kx and Ky means that the dispersion parallel to the coastline is much higher than that normal to the coast. Theoretical Approach for the Model

There is a great diversity in the flow patterns that may develop as the discharged wastes stream mixes with the ambient water. These flow patterns will determine the configuration, size and intensity of the mixing process as well as any impact of the discharged pollutants on water body surface, bottom and shoreline of the critical area. In the early decades of th the 20 century, many scientists developed theories and equations, which deal with the physical and hydrodynamic features of mixing process between discharge effluents and their receiving ambient water bodies. The behavior of an outfall discharge wastewater at a certain depth beneath sea surface can be divided onto three successive steps:

1. Near the orifice of the outfall, the discharges are influenced by jet momentum.

2. Near the sea surface, the discharges are governed by buoyancy.

3. The dispersion and disposition of discharges by current action from their final level away to the surrounding areas. The construction of the used constructive model was based on the equations described by Bonazountas (1987) to estimate initial dilution and waste water field. Bonazountas (1987) calculated the thickness of the wastewater field by assuming that, the average dilution at elevation Y is proportional to Y, this is nearly true in all buoyant line plume and jet cases for Y sufficient Application of a Mathematical Model for … 155 large, thus assuming Saw is the same as the average dilution in the plume and Y=Yb then Saw/Sa = Yb/Ymax , where:

Saw : the average dilution in the waste field Sa : the calculated value at the top of the plume. Yb : the Y coordinates at the bottom of the waste field Ymax: the Y coordinates at the top of the waste field. Qo : the source discharge (volume flux). U : the current speed. b : the width of the waste field (normal to the current). g : gravitational acceleration. q : outfall discharge. Δρ :residual density difference at top of plume. ρ :density of seawater.

The thickness of the waste field is

H =Ymax – Yb =Ymax (1-Yb /Ymax) or

H = Ymax((QoSa / ubYmax)/(1 +QoSa /ubYmax))

To apply this equation, it is necessary to determine u, b and Sa, also, it should be noted that:

If the product (ub) is very small, then the analysis fails because the assumption Saw /Sa =Yb/Ymax break down for a current which is perpendicular to the diffuser, (b) may simply taken as the length of the diffuser.

For a current parallel to the diffuser, the width (b) must be estimated differently.

For uniform ambient density, Bonazountas (1987) demonstrated that: The actual average dilution including the effect of blocking is given by:

1/3 Saw = 0.38 (gq) Ymax / q

G = g (Δρ / ρ)

And the average dilution Sa in the absence of blocking is given by:

1/3 Sa =0.54 (Gq) Ymax/q

156 Ebtessam E.E. Mohamed et al.,

Model Application

The different equations of Bonazountas (1987) were applied to obtain the pass of the pollutant and its decay on the sea surface under two conditions: - Considering the ambient water moving with constant speed and direction (u=5cm/sec & v=1 cm/sec). - Using two sets of current data collected in September 29 and 30 1994 in front of Alexandria. Applying the first condition from which one can obtain that pollutant layer is spread in ellipse form as shown in Fig. 4 and 5. Applying the second condition using the two different sets of current vectors measured in front of Alexandria.

The first set of current data was collected in sep.29, 1994 in front of Alexandria. In the region situated between El-Agamy and Abu Qir. The dominant current direction for this set of data was NWW, N and NNE their speed in the range from 4 to 49 cm/sec.

The outcoming pattern at the sea surface for zinc and suspended particles discharged from Kayet Bay outfall at 8m depths (Fig. 6-a and b) Showed disposition and dispersion of waste discharges northwards and eastwards with a bulk transport in the SEE direction, i.e., the coastline of Mostafa Kamel was a target of pollutant accumulations. Also the out coming patterns of zinc and suspended particles discharged directly on the sea surface from El-Umum drain, (Fig. 7-a and b), showed disposition and dispersion of waste discharges northward and eastwards, also, with a bulk transport in the SEE direction, thus pollutants would be carried into the western Harbour and Kayet Bay near shore water as well as the open sea. The analysis and current data for the second current set, which was collected in September 30, 1994, are shown in the following table:

Number of readings Direction Intensity (cm/sec) 25 NNE, NE, NEE 2-5 19 E 3-4 14 SEE, SE, SSE 4 9 NW, NNW 13-37 2 SSW, SW 2-9 1 W 3 Application of a Mathematical Model for … 157

) m ( Distance

Distance (m)

Fig. 4. Distribution pattern of suspended particles plume of Kayet Bay outfall at the surface at u=5 cm/sec and v=1 cm/sec.

Distance (m)

Distance (m)

Fig. 5. Distribution pattern of Zinc plume of El Umum Drain at the sea surface at u= 5cm/sec and v=1 cm/sec.

158 Ebtessam E.E. Mohamed et al., ) m ( Distance

Distance (m)

Distance (m) Fig. 6. Distribution pattern of a) Zinc and b) Suspended particles on the sea surface at Kayet Bay using current data measured in front of Alexandria on September 29, 1994. Application of a Mathematical Model for … 159 ) m ( Distance

) m ( Distance

Distance (m) Fig. 7. Distribution pattern of a) Zinc and b) Suspended particles on the sea surface in front of El-Umum Drain using current data measured in front of Alexandria on September 29, 1994.

160 Ebtessam E.E. Mohamed et al., Distance (m) Distance (m)

Distance (m) Fig. 8. Distribution pattern of a) Zinc and b) Suspended particles on the sea surface at Kayet Bay using current data measured in front of Alexandria on September 30, 1994. Application of a Mathematical Model for … 161 Distance (m) Distance (m)

Distance (m) Fig. 9. Distribution pattern of a) Zinc and b) Suspended particles on the sea surface in front of El-Umum Drain using current data measured in front of Alexandria on September 30, 1994.

162 Ebtessam E.E. Mohamed et al.,

The out coming patterns at the sea surface for zinc and suspended particles discharged from Kayet Bay outfall at 8m depth are illustrated in Fig. 8. The transport of this discharged wastewater was the major in the NE, E and SE directions and the minor in S direction, while it was minimal in the N and W directions. Thus pollutant would be carried towards the coastal zone of Kayet Bay and into the eastern Harbour. The patterns of waste discharges for zinc and suspended particles from El-Umum drain came in complete accordance with the current vectors illustrated in the table above. The disposition of discharges was the largest in the NNE, NE, NEE, E, SEE, SE and SSE directions and the smallest in the NNW and NW directions, (Fig. 9). While it was the minimal in the W direction. The discharged waste from El-Umum drain would be dispersed in the open sea and carried into the western Harbour. Model Verification

Regarding the observed current data obtained during the cruise of December 1996, it is clear that, in front of Kayet Bay, the prevailing directions of currents were, SSW, SW and SWW with their intensities in the range from 3 to 7.2cm/sec, while in El-Max sector, the prevailing directions were N, S and SW with their intensities in the range from 7.2 to 11.7cm/sec. Applying the model when the exponential part is equal to 1. A constant concentration was calculated at certain times with difference in traveling space which depends on the current speed and direction, as shown in Fig. 10 and 11. Considering the exponential part with its real calculated values since it depends on the current speed and the effluents discharged through Kayet Bay outfall. The patterns of dispersion of zinc and suspended particles on the sea surface are illustrated in Fig. 12. The discharge water (Fig. 12) was disposed and dispersed towards the south and the west with bulk transport in the SSW, SW and SWW directions, i.e. the waste discharges were carried towards the coastal area from Kayet Bay to Ras El-Teen, following the current direction. Figure 13 illustrates the obtained patterns of dispersion of zinc and suspended particles discharged directly on the sea surface from El Umum drain. It was obvious that the waste discharges from El-Umum drain during December 1996 were disposed and dispersed northwards and Application of a Mathematical Model for … 163 westwards with bulk transport in SW direction i.e. the waste discharges were carried towards the coastal area from El- Dikhila to El- Agamy beach.

Fig. 10. Distribution pattern of waste water at sea surface in front of Kayet Bay outfall and its decay with time.

Fig. 11. Distribution pattern of waste water at sea surface in front of El Umum Drain and its decay with time.

164 Ebtessam E.E. Mohamed et al.,

) m ( Distance Distance (m)

Distance (m)

Fig. 12. Distribution pattern of a) Zinc and b) Suspended particles on the sea surface at Kayet Bay using current data measured in front of Alexandria on December 30, 1996. Application of a Mathematical Model for … 165 ) m ( Distance

Distance (m)

Distance (m) Fig. 13. Distribution pattern of a) Zinc and b) Suspended particles on the sea surface in front of El-Umum Drain using current data measured in front of Alexandria on September, 1996.

166 Ebtessam E.E. Mohamed et al.,

Conclusion

In the case of Kayet Bay outfall: 1. Pollutants would be carried to the coastal area between Kayet Bay and Ras El-Teen if currents in front of Kayet Bay were SSW or SW or SWW. 2. Pollutants would be carried to Moustafa Kamel coastline if currents were SEE. In the case of El-Umum drain: 1. Waste discharges would be carried to El-Dikhila area if currents were SW. 2. Waste discharges would be carried to the western Harbour if currents were SSE or SE or SEE. 3. Waste discharges would be carried towards Kayet Bay region if currents were NE.

References

Abdel-Halim, A.M. (2004) Chemical equilibria studies of some phosphorus compounds under different conditions, Ph. D. Thesis, Faculty of Science, Alexandria University, Egypt. Abdullh, A.M. (1979) Study of the currents and hydrographic structure of the water mass in front of Alexandria coast. M. Sc. Thesis, Faculty of Science, Alexandria University, Egypt. Abuoul-Dahab, O. and Halim, Y. (1986) Impact of land – Based sources on chromium species and concentrations in coastal waters West of Alexandria, Rapp. Comm., Int. Mer Medit., 30 (2): 108. Bonazountas, M. (1987) Training Course on Modeling Outfall & Coastal Water Quality, Athens, Greece. El-Naggar, M.F. (1994) Evaluation of chemical effects on the new sewage outfall in Kayet – Bay area on seawater quality, M. Sc. Thesis, Alexandria University. Farag, M.M. (1982) Circulation pattern and hydrographic structure of El-Max and Western Harbour areas, M. Sc. Thesis, Faculty of Science, Alexandria University, Egypt. Farag, M.M., Aboul-Naga, W.M., Shalaby, E.A., Baghdadi, H.H. and Kamel, M.S. (2000) Tracing some water quality after the construction of El- New harbour using Remote Sensing, Bull. Nat. Inst. of Oceanogr. & Fish. A. R. E., 26: 43-70. List, E.J., Gartrell, G. and Winnant, C.D. (1990) Diffusion and Dispersion in Coastal Waters, J. Hydr. Engrg., 116 (10): 1158-1179. Okubo, A. (1974) Some speculation on oceanic diffusion diagrams, Repport et process-Verbaux des Reuniions, Commission Internationale Pour L,exploration scientifique de la Mer., 167: 77-85. Rifaat, A.E. (1982) Bottom sediments of the Western Harbour of Alexandria and its neighborhood, M. Sc. Thesis, Faculty of Science, Alexandria University, Egypt. Said, M.A., El-Deek, M.T.H. and Shriadah, M.A. (1994) Effect of pollution on the hydrochemical characteristics of different water types in El-Max area West of Alexandria, Egypt, Acta Adriat., 34(1/2): 9 – 19. Application of a Mathematical Model for … 167

* – – – $% – !" # – *

. # $ ! " " " * +, -# ) %& ' '( 0( ! - ./ !% / 3! .1& -# " " 6 $ / , 4 5' '<6 *, 4' 9:;; ' 47, 8 " < <) "(< >?) "4 5' ' *6 !,/ ./ .6/ +, */ - #< > ' , " *6 ( < . < / - $ " !% 4 47, 4 -. ! < " . 4 / (