A Case Study of Contaminant Transport Modelling Tuzla Oil Spill
A Case Study of Contaminant Transport Modelling Tuzla Oil Spill
A Case Study of Contaminant Transport Modelling Tuzla Oil Spill
ABSTRACT
On February 13th, 1997, a tanker named TPAO exploded in Tuzla shipyards located on the
northeastern coast of the Sea of Marmara. During the fire, an estimated amount of 215 tons of oil was
VSLOOHG LQWR WKH $\GQON %D\ %DVHG RQ ILHOG REVHUYDWLRQV DQG FRPSXWHU VLPXODWLRQV WKH SUHVHQW
paper investigates the water circulation and the distribution of contaminant concentrations in time
steps following the accident. Two separate model grids are developed for the Sea of Marmara and the
$\GQON %D\ WR DFFRXQW IRU ERWK WKH ODUJH VFDOH HIIHFWV RI WKH ZDWHU FLUFXODWLRQ DQG WKH ORFDO
conditions determined by the physical boundaries. The simulation results for the distribution of
surface currents and the fate of oil spill showed that the contamination mainly affected the immediate
YLFLQLW\ RI WKH $\GQON %D\ 7KH VKDSH RI WKH ED\ DQG WKH ORFDO ZLQG DQG FXUUHQW FOLPDWH DUH IRXQG WR
be the key factors to limit contamination within the bay waters. Model runs for possible oil spill
scenarios in the Sea of Marmara showed that the most critical point is located at the southern entrance
of the Strait of Istanbul. At this location, strong coastal currents generated by the natural channel
hydrodynamics, force the contaminant to take a shore parallel course either westerly or easterly
depending on the wind induced circulation in the Sea of Marmara.
INTRODUCTION
Marmara Region (Figure 1) is an important coastal settlement in Turkey with rapidly increasing
population and industrial activities. Being enclosed by this most industrialized region in Turkey, the
Sea of Marmara and the Turkish Straits are subject to intensive navigation activity. With the recent
increases in sea traffic, these waterways have become a prime site for oil spill pollution. In the
Marmara Region, nearly 450 sea accidents have been reported within the last forty years. Some of
these accidents resulted in historic oil spills with severe impacts on the environment. The present
study investigates a recent accident in Tuzla, which is located off the Northeast Coast of the Sea of
Marmara.
In this study, two separate model grids are developed, one for the Sea of Marmara and another one
IRU WKH $\GQON %D\ WR VLPXODWH WKH ZDWHU FLUFXODWLRQ 8VLQJ WKH UHVXOWV RI WKH K\GURG\QDPLF PRGHO D
contaminant transport model predicts the location and the concentration of the contamination for
different time projections following the accident.
The main objective of the computer modelling is the prediction of the contaminant behavior, which
is the key element for oil pollution control in case of an accident. This provides the opportunity to
take the necessary immediate precautions by means of containment and removal. Once the model
grids are developed, this study enables to develop scenarios for different physical conditions in the
Sea of Marmara and the Turkish Straits (wind speed and direction, boundary conditions, type and
amount of the spill etc.). The results of such a study can be used as a basis for a Regional Oil Spill
Contingency Plan.
1
Grad. Asst., Dept. of Civil Engg, %RD]LoL 8QLY ,VWDQEXO 7XUNH\
2
*UDG $VVW ,QVW RI (QY 6FL %RD]LoL 8QLY ,VWDQEXO 7XUNH\
3
$VVW 3URI 'HSW RI &LYLO (QJJ %RD]LoL 8QLY ,VWDQEXO 7XUNH\
Figure 1. Location Map of the Sea of Marmara and the Turkish Straits
The first of these models, RMA-2, is a two-dimensional, depth averaged, free surface finite
element model, which can simulate the current circulation. After a finite element mesh has been
constructed and boundary conditions and material properties have been defined, the water surface
elevation and flow velocity at each grid point can be computed.
Based on the hydrodynamic solution obtained by RMA-2, a second numerical model, RMA-4 is
used to simulate the contaminant transport. The contaminant transport model requires as input the
initial spill conditions as a set of point loads in addition to the physical parameters used in the
hydrodynamic model.
h ( uh ) ( vh )
Continuity Equation: + + =0
t x y
u u u h xx 2u xy 2u sx bx
+u +v = g +
t x y x x 2 y 2 h h
v v v h yx 2 v yy 2 v sy by
+u +v = g +
t x y y x 2 y 2 h h
where,
t is the time
x, y are the horizontal directions
u, v are the water particle velocities
h is the water depth
g is the gravitational acceleration
xx, xy, yx, yy are the turbulent exchange coefficients
sx, sy are the surface shear stresses due to wind
bx, by are the bottom shear stresses
is the water density
For the Marmara Model, the boundary conditions are given in terms of head difference between
the Strait of Istanbul (Bosphorus) and the Strait of anakkale (Dardanelles) and the outflow rate at
Dardanelles. These parameters are assumed to have relatively small effects on the final distribution of
the contaminant concentrations in the fine grid model with respect to the main forcing caused by the
wind shear. Wind conditions are taken for the accident day and the following two days. As a result,
the distribution of the water velocities in the Sea of Marmara is obtained for 72 hours (Figure 5.).
For the Tuzla Model, the hydrodynamic conditions at the open boundary are obtained from the
results of the Marmara Model. The same wind conditions are used previously in the Marmara Model.
As a result, the current circulation in AyGQON %D\ LV obtained for 72 hours of wind data (Figure 3).
c c c 2c 2c
h( + u + v Dx 2 Dy 2 + kc ) = 0
t x y x y
where,
MODEL RESULTS
:KHQ WKH RLO VSLOO VLPXODWLRQ LV FRQGXFWHG IRU WKH $\GQON %D\ ZLWK WKH H[LVWLQJ FRQGLWLRQV RI WKH
DFFLGHQW LW LV REVHUYHG WKDW WKH VSLOOHG RLO KDV QRW EHHQ VSUHDG RXWVLGH WKH $\GQON %D\ )LJXUH
This is due to the wind conditions, the amount of spill and the closed current patterns (Figure 3). Field
observations confirm same type of behavior for the fate of the contamination indicating that the
occurred spill had no large scale effect on other regions of the Sea RI 0DUPDUD 2NX HW DO
+RZHYHU LW ZDV QRW SRVVLEOH WR IRUHFDVW WKH VHYHULW\ RI LWV ORFDO HIIHFWV RQ WKH $\GQON %D\
As it is stated before, the oil spill simulation is conducted for different locations (Silivri, Tekirda,
DUN|\, Karabiga, Gemlik, Bosphorus) in the Sea of Marmara with the existing conditions after the
modelling of current circulation (Figure 5). The evaluation of the findings of the Marmara Model put
forth that the most critical point in terms of pollution is the entrance of the Bosphorus to the Sea of
Marmara (Figure 6a, 6b). This result was expected because of the strong surface currents at the
entrance point.
Figure 4. Contaminant Concentration Distribution in AydQON %D\ 33 hours after the spill)
Figure 5. Current Circulation in the Sea of Marmara (Feb. 13th,1997)
Figure 6a. Contaminant Distribution in the Sea of Marmara (12 hours after the spill)
Figure 6b. Contaminant Distribution in the Sea of Marmara (33 hours after the spill)
RESPONSE AND RECOVERY
Ideally there are mainly six response alternatives to combat an oil spill (OSRL, 1995): (1) Monitor
and evaluate, (2) disperse the oil with chemicals, (3) contain and recover the oil at sea, (4) protect
vulnerable resources, (5) burning and (6) clean-up the shoreline.
In the case of the TPAO oil spill, it is seen that an ideal containment couldnt be performed,
LQVWHDG WKH RLO ZDV FROOHFWHG E\ VNLPPHUV RQO\ ZKHUH HYHU LW ZDV REVHUYHG DW WKH VHD VXUIDFH 6DUND\D
et al., 1997). However, this operation was not fully effective to control the spill and the oil came to
the shore. The shore portion of the oil was collected by mechanical pick-up equipment and through
hydraulic cleaning. In Tuzla, the effectiveness of these clean-up operations is questionable due to the
continuous contamination coming from the shipyards. The lack of the prediction of oil movement
made it more difficult to control the spill.
During an oil spill accident, the following initial actions should be taken (OSRL, 1995); source
isolation and containment, data collection, prediction of the trajectory of the oil spill, spill
surveillance, consideration of strategies and equipment requirements and finally the response. These
steps couldn't be carried out for the TPAO Spill, therefore it cant be called a planned response. The
present study is an example of modelling efforts to help to predict the movement of the spill
immediately after the accident following the data collection. This may be used for determining the
strategy alone in the case where detailed surveillance is not possible and the response (containment,
equipment planning etc.) can be planned according to the modelling results.