Three-Dimensional Simulation of Qeshm Channel Currents
Total Page:16
File Type:pdf, Size:1020Kb
Journal of the Persian Gulf (Marine Science)/Vol. 2/No. 3/March 2011/8/9-16 Three-Dimensional Simulation of Qeshm Channel Currents Mahmoudov, Masoud1*; Chegini, Vahid1; Montazeri Namin, Masoud2 1- Iranian National Institute for Oceanography, Tehran, IR Iran 2- Tehran University, Tehran, IR Iran Received: April 2010 Accepted: December 2010 © 2011 Journal of the Persian Gulf. All rights reserved. Abstract Qeshm Channel is a shallow and narrow waterway located between Qeshm Island and the mainland in the vicinity of Hormuzgan Province, Islamic Republic of Iran. This channel is important because of its economic, industrial, fisheries and navigation role it plays as well as environmental issues it presents in the region. A prognostic study was performed to simulate currents in this channel, using COHERENS model. This model is a three-dimensional hydrodynamic model. Simulation of currents was carried out in 20 sigma levels from the seabed to the water surface during one month. It was assumed that the variations of seawater temperature and salinity, four main tidal constituents and the regional wind were the most effective factors in the numerical simulation. Finally, sensitivity analysis was carried out for each factor and the outputs of simulation were verified using the field data recorded by the experts from Iranian National Center for Oceanography. Very good agreements were found between the numerical results and the field data. Keywords: Tide, Physical parameters, Currents, Qeshm channel Downloaded from jpg.inio.ac.ir at 5:28 IRST on Monday September 27th 2021 1. Introduction example, modeling of North Sea (Patrick and Luyten, et.al, 1999), modeling the impact of the In the past decades, a variety of mathematical and Scheldt and Rhine/Meuse plumes on the salinity hydrodynamic models were developed as demands distribution in Belgian waters (Lacroix and Ruddick, from scientists and researchers increased. Some of 2004), hydrodynamic investigation of Sacca di Goro these models were developed based on fundamental coastal lagoon (Marinov and Norro, 2005), equations of motion such as Navier-Stocks circulation of the Persian Gulf (Kämpf and equations. COHERENS (COupled Hydrodynamic Sadrinasab, 2005), dispersion of contaminant from Ecological REgioNal Shelf Seas), however, was a Karoon estuary in the Persian Gulf (Poorkiani and three-dimensional, N-S based hydrodynamic model Sadrinasab, 2007) and three dimensional dispersion which also included biology, contaminant and of pollutants in the Persian Gulf (Mahmoudi and sediment modules. Sadrinasab, 2008). Furthermore, this model was also COHERENS has been used by many investigators used in a study that investigated the two- in the Persian Gulf region and other seas; for dimensional modeling of Khuran Strait employing CECAD in Cartesian coordinate (Ghiyasi, 2006). * E-mail: [email protected] Sabagh Yazdi and Sharbati(2008) applied NASIR Mahmoudov et al. / Three-Dimensional Simulation of Qeshm Channel Currents software in modeling of tidal currents in Qeshm Records of weather stations at Qeshm airport (N: Channel, In this study, however, the hydrodynamic 26º45’- E: 55º54’) and synoptic stations located module of COHERENS was applied to study the south of Qeshm city (N: 26º57’– E: 56º16’) were current hydraulics of the Qeshm Channel. used as meteorological data in the simulation. Since the direction of prevailing wind was north-south, 2. Materials and Methods perpendicular to the channel, the effect of wind on currents was considered not to be intensive. Climate 2.1. Description of the Study Area of this region is usually warm and humid. In summer, heat reaches up to 45 ºC and in winter, Hara Biosphere Reserve and Mangrove forests in temperature rarely falls below 12 ºC. The average the Qeshm area are important ecosystems of national number of rainy days per annum is about 12 days and international importance in the region; however, and the annual precipitation is less than 175 mm. Qeshm Channel is located in a region ecologically Salinity of water in this region is influenced by vulnerable to damaging activities such as, oil spill, currents from the Persian Gulf and the Sea of Oman inland tank leakage and ships ballast water flushing (Johns, 1998). Local variation of salinity is poor due and development occuring around it 0 Because of to narrow width and shallow depth of the channel. construction of several ports in this region, gaining Average salinity is 38 PSU and seasonal fluctuation knowledge of current induced sediment transport in is about 4 PSU. Average seasonal water temperature this region is quite necessary. The study area fluctuates between 22 ºC in winter up to 32 ºC in extended from 55º16’ to 56º16’ longitude and 26º39’ summer with a diurnal fluctution of about 0.2 ºC to to 27º10’ latitude (Fig. 1). The east-west bound 1.2 ºC (Hajizade Zaker, 2005). channel is 110 Km long and is extended between Tide is mixed type because of the influence of the Bandar Abbas and Qeshm City towards Bandar- Persian Gulf (Najafi, 1997). The main diurnal and Lengeh and Bandar-Khamir and Basaeedu area. semidiurnal constituents are M2, S2, O1 and K1 Maximum and minimum widths of the channel are (Foreman, 1977). Each constituent has its specific 25 and 3.5 Km, between Bandar-abbas and Qeshm amplitude, period and phase. However, amplitudes Downloaded from jpg.inio.ac.ir at 5:28 IRST on Monday September 27th 2021 city and between Pohl and Laft ports in Khuran and phases differ per location. Tidal characteristics at Strait, respectively. Maximum depth of the channel Bandar-Abbas (east of channel) and Basaidu (west of is about 33 m near eastern throat and overall it is channel) are presented in Table 1. considered shallow. Table 1- Characteristics of Main Tidal Constituents at Bandar- Abbas and Basaidu m2 s2 o1 k1 (deg) (deg) (deg) (deg) (deg) (deg) (deg) (deg) A(m) A(m) A(m) A(m) T(h) T(h) T(h) lOcation lat log T(h) 12.42 12.42 293.8 12.00 349.8 25.82 23.93 Bandar- 26º 55º 1.01 0.43 0.19 54.8 0.37 77.3 Abbas 39’N 16’E 12.42 12.42 339.0 12.00 25.82 23.93 125.6 0.85 0.85 0.37 0.25 88.3 0.30 27º 56º 0.4 Basaidu 10’N 16’E Fig. 1. Geographic Situation of the Qeshm Channel As for biological aspects, this area has the greatest 10 Journal of the Persian Gulf (Marine Science)/Vol. 2/No. 3/March 2011/8/9-16 Hara forest in the Persian Gulf which is the breeding S S S S u v w ground and habitat for seasonal migrant sea birds t x x x 1 2 3 (6) C S C S C S such as herons and flamingoes. More than 93 bird D S T D S T D S T D T T D H T D H T x x x x x x species and 32 species of fish, marine turtles, marine 3 E 3 U 1 E 1 U 2 E 2 U snakes, jackals and hyenas inhabit the region. Also, where u, v and w are velocity components of the several species of crustaceans, lobsters and mollusks current, T stands for temperature, S the salinity, live in the area. f 2sin the Coriolis frequency, 2 / 86164 rad/s the earth rotation frequency, g the gravity 3. The Model acceleration, p the pressure, T and T the vertical eddy viscosity and diffusion coefficients, H the 3.1. Governing Equations horizontal diffusion coefficient for salinity and temperature, the density, 0 a reference density, In this study, simulation was performed by the c p the seawater specific heat at constant pressure and three-dimensional version of COHERENS. This I (x1 , x2 , x3 , t) solar irradiance. The horizontal model is based on vertical sigma coordinate. components of the stress tensor are defined by: The model is established on the hydrostatic version of Naveir-Stokes equations, continuity and u 2 11 H (7) standard equation for scalar parameters. The x1 Boussinesq approximation is included in horizontal C S D u v T momentum equations. The governing equations are 21 12 H D T (8) E x x U as follows: 2 1 v 22 2 H (9) u u u u x u v w fv 2 t x1 x2 x3 (1) C S 1 p D u T where H is the horizontal diffusion coefficient D T T 11 21 x x x x x 0 1 3 E 3 U 1 2 for momentum. Downloaded from jpg.inio.ac.ir at 5:28 IRST on Monday September 27th 2021 v v v v 3.2. Mesh Generation and Bathymetry u v w fu t x1 x2 x3 C S (2) 1 p D v T Considering extension and water mass of Qeshm D T T 12 22 x x x x x 0 2 3 E 3 U 1 2 Channel, computational domain was generated with 100 x 60 uniform structural square meshes with a p length of 1000m in both x and y directions (Fig. 2). g (3) x3 Also, 20 sigma levels were assumed from seabed to surface in vertical to correspond the location of current metering instruments at different depth to u v w 0 (4) computational cells. Coordinates of computational x1 x2 x3 cells vary from (i=1, j=1, k=1) to (i=100, j=60, k=20), in which i, j and k are standard Cartesian T T T T u v w coordinate indices. t x x x 1 2 3 (5) C S C S C S Hydrographic data was prepared by transforming 1 I D T T D T T D T T D T T D H T D H T 0cp x3 x3 E x3 U x1 E x1 U x2 E x2 U Iranian National Geographical Organization 11 Mahmoudov et al.