Table of Contents
ISRN Civil Engineering
Volume 2012 (2012), Article ID 345979, 11 pages
Research Article

A Barotropic Model of the Red Sea Circulation

Faculty of Engineering, Cairo University, Giza 12613, Egypt

Received 11 January 2012; Accepted 2 February 2012

Academic Editors: I. Raftoyiannis and N. Xie

Copyright © 2012 Mohamed M. El-Shabrawy et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


A linear depth-averaged numerical model of the horizontal barotropic circulation of the Red Sea is developed. The flow considered is incompressible, quasi-steady and free from vertical buoyancy currents and horizontal density currents. The model predicts general water circulation as affected by friction stresses at the irregular Red Sea bottom and coastal features, variable Coriolis force, vertical, and lateral turbulent friction, and seasonal non-uniformity of wind stress at the sea surface. The implicit finite-difference solution of the boundary value problem is verified with previous solutions for a rectangular constant-depth sea and for an elliptical lake with uniformly-sloping bottom topography. The model output is shown as plots of seasonally-averaged mass-transport vectors and circulation streamlines in the Red Sea basin.

1. Introduction

The Red Sea is 2040 km length, 280 km averaged width, and 491 m averaged depth. Also, its latitude extent between 12°N and 28°N warrants an effective and variable Coriolis force, the β-effect [1]. The two-dimensional models for this effect on ocean currents were devised for constant-depth rectangular basins [2, 5]. The combined action of depth variation and the β-effect on the two-dimensional currents were also modeled [1, 3, 4]. The driving wind was assumed constant [3, 4] although it is variable [1]. The bottom friction was modeled by linear mass-transport dependence [3, 4] while lateral turbulent friction was neglected. The lateral turbulent friction was modeled, in [1], by linear dependence on mass- transport dependence. This assumption had no effect the order of the elliptic partial differential equation governing the depth-integrated circulation streamline. Analytical solutions were obtained, in [1], for constant depth basins with variable Coriolis parameter and spatially-variable wind stress distribution. All of these two-dimensional models [1, 35], have neglected lateral turbulent friction.

This two-dimensional model considers the manner by which barotropic circulation in the Red Sea is affected by its irregular bottom bathymetry, latitude variation of the Coriolis force, vertical and lateral turbulent friction stresses, and seasonal variation of the surface wind stress. The mass and momentum conservation equations are vertically integrated and combined to form a single partial differential equation for the mass transport stream function. Bottom friction is modeled by linear mass- transport dependence while lateral turbulent friction has a third order mass- transport dependence. The external driving force acts at the sea surface as a non-uniform seasonally-averaged wind stress. The governing equations of the three-dimensional flow are modified by the model assumptions of incompressible and hydrostatic flow. This allows the derivation of the governing equations of the two-dimensional baro-tropic (homogeneous sea water density) and quasi-steady circulation. The depth-integrated equations are hence manipulated to result in a bi-harmonic equation of the mass-transport stream function.

A finite difference scheme formulation of the bi-harmonic equation over a Cartesian mesh with variable grid size is derived. The model is verified versus simple cases of known analytical solutions. The test cases include both constant and variable depth basins with uniform and variable wind stress distributions. The model numerical results of the circulations of the Red Sea (with its Spatial wind stress distribution, irregular coastal configuration and bottom topography) are presented in plots of mass-transport vectors and circulation streamlines.

2. Governing Equations

The continuity equation for the compressible flow, first derived by Leonhard Euler, [8] is 𝜕𝜌+𝜕𝑡𝜕(𝜌𝑢)+𝜕𝑥𝜕(𝜌𝑣)+𝜕𝑦𝜕(𝜌𝑣)𝜕𝑧=0,(1) where 𝜌 is the fluid density, 𝑢,𝑣,𝑤 are velocity components in 𝑥,𝑦,𝑧 directions, and 𝑡 is the time. Also, the momentum equation [6] in the Cartesian form is𝜕𝑢𝜕𝑡+𝑢𝜕𝑢𝜕𝑥+𝑣𝜕𝑢𝜕𝑦+𝑤𝜕𝑢1𝜕𝑧=𝜌𝜕𝑝𝜕𝑥+𝑓𝑣2𝜔cos𝜑𝑤+𝜇𝑣𝜕2𝑢𝜕𝑧2+𝜇𝜕2𝑢𝜕𝑥2+𝜕2𝑢𝜕𝑦2,(2a)𝜕𝑣𝜕𝑡+𝑢𝜕𝑣𝜕𝑥+𝑣𝜕𝑣𝜕𝑦+𝑤𝜕𝑣1𝜕𝑧=𝜌𝜕𝑝𝜕𝑦+𝑓𝑢+𝜇𝑣𝜕2𝑣𝜕𝑧2+𝜇𝜕2𝑣𝜕𝑥2+𝜕2𝑣𝜕𝑦2,(2b)𝜕𝑤𝜕𝑡+𝑢𝜕𝑤𝜕𝑥+𝑣𝜕𝑤𝜕𝑦+𝑤𝜕𝑤1𝜕𝑧=𝜌𝜕𝑝𝜕𝑧𝑔2𝜔cos𝜑𝑢+𝜇𝑣𝜕2𝑤𝜕𝑧2+𝜇𝜕2𝑤𝜕𝑥2+𝜕2𝑤𝜕𝑦2,(2c)where 𝑝 is the pressure field, 𝑥,𝑦,𝑧 are Cartesian coordinate system fixed with respect to the rotating earth, with 𝑧 upward, and 𝑓 is the Coriolis parameter. 𝜇,𝜇𝑣 are depth averaged coefficients that model the momentum exchange due to eddies generated by the lateral turbulent and vertical mixing respectively, and 𝜑 is latitudinal angle.

3. Model Assumptions

To build a model of the Red Sea circulation, we assume the following.(1)The time variations on the daily, weekly, and monthly scales are averaged. Consequently, all physical quantities are seasonally averaged and the flow is considered “quasi-steady” 𝜕𝜕𝑡=0.(3)(2)Turbulent shear forces in the momentum equation ((2a), (2b), and (2c)) are simulated by depth averaged eddy (momentum exchange) coefficients 𝜇,𝜇𝑣 for lateral turbulent and vertical mixing, respectively. These coefficients have constant values.(3)All types of wave –like and small-scale motions are neglected. Only current motion and large scale circulation are considered.(4)The sea water density is assumed homogenous and hence density currents are neglected. (5)The atmospheric pressure is assumed uniform over the Red Sea.(6)The flow is assumed incompressible the vertical profile of the pressure field is assumed to be hydrostatic. Vertical acceleration is neglected.(7)The Red Sea basin latitudinal extent is large enough to allow the Coriolis force and its latitude variation (the beta effect) to have a significant effect.(8)Linearization is achieved by assuming the convective (non-linear) acceleration much smaller than the Coriolis (linear) acceleration [8].

4. Barotropic Model Equations

According to the model assumptions, the governing equations take the form [4, 6],𝜕𝑢+𝜕𝑥𝜕𝑣+𝜕𝑦𝜕𝑤𝜕𝑧=0,(4)1𝑓𝑣=𝜌𝜕𝑝𝜕𝑥+𝜇𝑣𝜕2𝑢𝜕𝑧2+𝜇𝜕2𝑢𝜕𝑥2+𝜕2𝑢𝜕𝑦2,(5)1𝑓𝑢=𝜌𝜕𝑝𝜕𝑦+𝜇𝑣𝜕2𝑣𝜕𝑧2+𝜇𝜕2𝑣𝜕𝑥2+𝜕2𝑣𝜕𝑦2,(6)𝜕𝑝𝜕𝑧=𝜌𝑔.(7) This set of linear partial equations in 𝑢, 𝑣, 𝑤 and 𝑝 are subject to boundary conditions; a wind stress at the mean sea surface, 𝑧=0, and a frictional stress at the sea bottom, 𝑧=𝐻𝜏𝑤𝑥,𝜏𝑤𝑦=𝜇𝑣𝜕𝑢,𝜕𝑧𝜕𝑣𝜕𝑧𝑧=0,𝜏𝑏𝑥,𝜏𝑏𝑦=𝜇𝑣𝜕𝑢,𝜕𝑧𝜕𝑣𝜕𝑧𝑧=𝐻.(8) Define the components of mass transport between the sea surface,𝜁, and the bottom, 𝐻, in the x- and y-directions as𝑀𝑥=𝜁𝐻𝑢𝑑𝑧,𝑀𝑦=𝜁𝐻𝑣𝑑𝑧.(9) Integrating equation (4) over 𝑧 between the boundaries 𝐻(𝑥,𝑦) and 𝜁(𝑥,𝑦), and using the conditions of integrity of the free surface and impermeability of the bottom [5, 6],𝜕𝑀𝑥+𝜕𝑥𝜕𝑀𝑦𝜕𝑦=0.(10) A mass transport stream function 𝜓(𝑥,𝑦)can then be defined𝑀𝑥=𝜕𝜓𝜕𝑦,𝑀𝑦=𝜕𝜓𝜕𝑥.(11) Integration of (7) over 𝑧 between the boundaries 𝐻(𝑥,𝑦) and 𝜁(𝑥,𝑦) gives𝑝=𝑝𝑎+𝜌𝑔(𝜁+𝐻),(12) where 𝑝𝑎is the atmospheric pressure.

Integration of (5) and (6) over 𝑧 between the boundaries 𝐻(𝑥,𝑦) and 𝜁(𝑥,𝑦) gives𝑓𝑀𝑦1=𝜌𝜕𝜕𝑥𝜁𝐻𝑝𝑑𝑧+𝜇𝑣𝜕𝑢𝜕𝑧𝜁𝜕𝑢𝜕𝑧𝐻+𝜇𝜕2𝑀𝑥𝜕𝑥2+𝜕2𝑀𝑥𝜕𝑦2,(13)𝑓𝑀𝑥1=𝜌𝜕𝜕𝑦𝜁𝐻𝑝𝑑𝑧+𝜇𝑣𝜕𝑣𝜕𝑧𝜁𝜕𝑣𝜕𝑧𝐻+𝜇𝜕2𝑀𝑦𝜕𝑥2+𝜕2𝑀𝑦𝜕𝑦2.(14) Dividing both equations by 𝐻 differentiating (13) for 𝑦 and (14) for 𝑥, and subtracting to eliminate pressure from the governing equations,𝜕𝑓𝜕𝑥𝐻𝑀𝑥+𝜕𝑓𝜕𝑦𝐻𝑀𝑦=𝜕𝜏𝜕𝑦𝑤𝑥𝐻+𝜕𝜏𝜕𝑦𝑏𝑥𝐻+𝜕𝜏𝜕𝑥𝑤𝑦𝐻𝜕𝜏𝜕𝑥𝑏𝑦𝐻+𝜇𝜕1𝜕𝑥𝐻𝜕2𝑀𝑦𝜕𝑥2+1𝐻𝜕2𝑀𝑦𝜕𝑦2𝜕1𝜕𝑦𝐻𝜕2𝑀𝑥𝜕𝑥2+1𝐻𝜕2𝑀𝑥𝜕𝑦2.(15) It can be easily simplified to be𝜇𝐻2𝐻𝜕4𝜓𝜕𝑥4𝜕+24𝜓𝜕𝑥2𝜕𝑦2+𝜕4𝜓𝜕𝑦4𝜕𝐻𝜕𝜕𝑥3𝜓𝜕𝑥3+𝜕3𝜓𝜕𝑥𝜕𝑦2𝜕𝐻𝜕𝜕𝑦3𝜓𝜕𝑦3+𝜕3𝜓𝜕𝑦𝜕𝑥2+𝑓𝐻2𝜕𝐻𝜕𝑥𝜕𝜓𝛽𝜕𝑦𝐻𝑓𝐻2𝜕𝐻𝜕𝑦𝜕𝜓=1𝜕𝑥𝐻𝜕𝜏𝑤𝑦𝜏𝜕𝑥𝑤𝑦𝐻2𝜕𝐻1𝜕𝑥𝐻𝜕𝜏𝑤𝑥𝜏𝜕𝑦𝑤𝑥𝐻2𝜕𝐻+1𝜕𝑦𝐻𝜕𝜏𝑏𝑦𝜏𝜕𝑥𝑏𝑦𝐻2𝜕𝐻1𝜕𝑥𝐻𝜕𝜏𝑏𝑥𝜏𝜕𝑦𝑏𝑥𝐻2𝜕𝐻.𝜕𝑦(16) Following the common assumption of the turbulent shear stress as linear function of horizontal mass transport [6, 7], using𝜏𝑏𝑥=𝑐𝑀𝑥and𝜏𝑏𝑦=𝑐𝑀𝑦 in (16) leads to𝜇𝐻2𝐻𝜕4𝜓𝜕𝑥4𝜕+24𝜓𝜕𝑥2𝜕𝑦2+𝜕4𝜓𝜕𝑦4𝜕𝐻𝜕𝜕𝑥3𝜓𝜕𝑥3+𝜕3𝜓𝜕𝑥𝜕𝑦2𝜕𝐻𝜕𝜕𝑦3𝜓𝜕𝑦3+𝜕3𝜓𝜕𝑦𝜕𝑥2+𝑐𝐻𝜕2𝜓𝜕𝑥2+𝜕2𝜓𝜕𝑦2+𝑓𝐻2𝜕𝐻𝑐𝜕𝑥𝐻2𝜕𝐻𝜕𝑦𝜕𝜓𝛽𝜕𝑦𝐻𝑓𝐻2𝜕𝐻𝑐𝜕𝑦𝐻2𝜕𝐻𝜕𝑥𝜕𝜓=1𝜕𝑥𝐻𝜕𝜏𝑤𝑦𝜏𝜕𝑥𝑤𝑦𝐻2𝜕𝐻1𝜕𝑥𝐻𝜕𝜏𝑤𝑥𝜏𝜕𝑦𝑤𝑥𝐻2𝜕𝐻,𝜕𝑦(17) where c is the bottom friction coefficient with the value taken to be either  .0025 or  .005, depending on the author. Studies [8, 9] have used 0.0025, whereas [10] has used a value of 0.005. Rao and Murty [11] have used values 0.0025, 0.0125, and 0.000625 in the modeling of the Lake Ontario. They found that as the bottom friction coefficient decreased the maximum mass transport increased.

Two independent methods have been used for the solution of the simultaneous set of inhomogeneous equations presented by (17). The first one is an exact method called the LU-decomposition method [12]. The second is an iterative method called the Liebmann accelerated point over-relaxation method [13]. The LU-decomposition method does not involve any iteration procedures and is consequently more accurate than the relaxation method; it may not be suitable for multilayer situations and for more sophisticated models than the present one simply because of computer storage limitations in dealing with large matrices. For this reason we used the LU-decomposition method to be more accurate.

5. Numerical Solution

When the basin is rectangular with sides parallel to grid lines, the standard difference replacements can be applied at all internal grid points. When the basin is non-rectangular, this is the default of all real basins; internal grid points adjacent to the boundary require special treatment. Such a grid point is depicted in Figure 1. Expanding the spatial derivatives using the methods of [12, 14] of (17) will lead to a 13-point finite difference equation of𝜓: 𝐶1𝜓𝑖2,𝑗+𝐶2𝜓𝑖1,𝑗+𝐶3𝜓𝑖,𝑗+𝐶4𝜓𝑖+1,𝑗+𝐶5𝜓𝑖+2,𝑗+𝐶6𝜓𝑖1,𝑗+1+𝐶7𝜓𝑖1,𝑗1+𝐶8𝜓𝑖,𝑗+1+𝐶9𝜓𝑖,𝑗1+𝐶10𝜓𝑖,𝑗+2+𝐶11𝜓𝑖,𝑗2+𝐶12𝜓𝑖+1,𝑗+1+𝐶13𝜓𝑖+1,𝑗1=RHS,(18) where𝐶1=𝜇𝐻2𝑖2,𝑗𝑠632𝑠1𝑠53𝑠7𝑠7+𝑠3𝐻1𝐻2𝜕𝐻𝜕𝑥𝑖2,𝑗+24𝐻𝑖2,𝑗4𝑠7𝑠7+𝑠3𝐻1𝐻2,𝐶2=𝜇𝐻2𝑖1,𝑗6𝑠7+𝑠32𝑠1𝑠53𝑠7𝑠3𝑠7+𝑠3𝐻2𝑠7+2𝑏2+𝑏4𝑠1+𝑠3𝜕𝐻𝜕𝑥𝑖1,𝑗𝐻𝑖1,𝑗244𝑠7𝑠3𝑠7+𝑠3𝐻2𝑠7+4𝑏3𝑏2+𝑏4+𝑐𝑏3𝐻𝑖1,𝑗+1𝑠1+𝑠3𝛽𝐻𝑖1,𝑗𝑓𝐻2𝑖,𝑗𝜕𝐻𝜕𝑦𝑖1,𝑗𝑐𝐻2𝑖1,𝑗𝜕𝐻𝜕𝑥𝑖1,𝑗,𝐶3=𝜇𝐻2𝑖,𝑗𝐻𝑖,𝑗244𝑠2𝑠4𝑠8+𝑠4𝑠2+𝑠6+24𝑠1𝑠3𝑠7+𝑠3𝑠1+𝑠5+𝑏2+𝑏4𝑏1+𝑏3𝜕𝐻𝜕𝑥𝑖,𝑗×6𝑠7+2𝑠32𝑠1𝑠53𝑠1𝑠3𝑠7+𝑠3𝑠1+𝑠5𝜕𝐻𝜕𝑦𝑖,𝑗6𝑠8+2𝑠42𝑠2𝑠63𝑠2𝑠4𝑠8+𝑠4𝑠2+𝑠62𝑐𝐻𝑖,𝑗𝑏1+𝑏3+𝑏2+𝑏4,𝐶4=𝜇𝐻2𝑖+1,𝑗6𝑠7+2𝑠3𝑠1𝑠53𝑠1𝑠5𝑠1+𝑠3𝐻12𝑏2+𝑏4𝑠1+𝑠3𝜕𝐻𝜕𝑥𝑖+1,𝑗𝐻𝑖+1,𝑗244𝑠1𝑠5𝑠1+𝑠3𝐻1+4𝑏1𝑏2+𝑏4+𝑐𝑏1𝐻𝑖+1,𝑗1𝑠1+𝑠3𝛽𝐻𝑖+1,𝑗𝑓𝐻2𝑖,𝑗𝜕𝐻𝜕𝑦𝑖+1,𝑗𝑐𝐻2𝑖+1,𝑗𝜕𝐻𝜕𝑥𝑖+1,𝑗,𝐶5=𝜇𝐻2𝑖+2,𝑗6𝑠7+2𝑠3𝑠13𝐻2𝐻2𝑠7𝑠1+𝑠5𝑠5𝜕𝐻𝜕𝑥𝑖+2,𝑗+24𝐻𝑖+2,𝑗4𝐻2𝐻2𝑠7𝑠1+𝑠5𝑠5,𝐶6=𝜇𝐻2𝑖1,𝑗+12𝑏2𝑏3𝐻𝑖1,𝑗+1𝑏2𝑠1+𝑠3𝜕𝐻𝜕𝑥𝑖1,𝑗+1𝜕𝐻𝜕𝑦𝑖1,𝑗+1𝑏3𝑠2+𝑠4,𝐶7=𝜇𝐻2𝑖1,𝑗12𝑏4𝑏3𝐻𝑖1,𝑗1𝑏4𝑠1+𝑠3𝜕𝐻𝜕𝑥𝑖1,𝑗1𝜕𝐻𝜕𝑦𝑖1,𝑗1𝑏3𝑠2+𝑠4,𝐶8=𝜇𝐻2𝑖,𝑗+1𝐻𝑖,𝑗+1244𝑠2𝑠6𝑠2+𝑠4𝐻14𝑏2𝑏1+𝑏3+𝜕𝐻𝜕𝑦𝑖,𝑗+16𝑠8+2𝑠4𝑠2𝑠63𝑠2𝑠6𝑠2+𝑠4𝐻3+2𝑏1+𝑏3𝑠2+𝑠4+2𝑐𝑏2𝐻𝑖,𝑗+1+𝑓𝐻2𝑖,𝑗+1𝜕𝐻𝜕𝑥𝑖,𝑗+1𝑐𝐻2𝑖,𝑗+1𝜕𝐻𝜕𝑦𝑖,𝑗+11𝑠2+𝑠4,𝐶9=𝜇𝐻2𝑖,𝑗1𝐻𝑖,𝑗1244𝑠8𝑠4𝑠8+𝑠4𝐻4𝑠84𝑏4𝑏1+𝑏3+𝜕𝐻𝜕𝑦𝑖,𝑗16𝑠8+𝑠42𝑠2𝑠63𝑠8𝑠4𝑠8+𝑠4𝐻4𝑠82𝑏1+𝑏3𝑠2+𝑠4+2𝑐𝑏4𝐻𝑖,𝑗1𝑓𝐻2𝑖,𝑗1𝜕𝐻𝜕𝑥𝑖,𝑗1𝑐𝐻2𝑖,𝑗1𝜕𝐻𝜕𝑦𝑖,𝑗11𝑠2+𝑠4,𝐶10=𝜇𝐻2𝑖,𝑗+2𝑠68+2𝑠4𝑠23𝐻4𝐻4𝑠8𝑠2+𝑠6𝑠6𝜕𝐻𝜕𝑥𝑖,𝑗+2+24𝐻𝑖,𝑗+24𝐻2𝐻2𝑠7𝑠1+𝑠5𝑠5,𝐶11=𝜇𝐻2𝑖,𝑗2𝑠642𝑠2𝑠63𝑠8𝑠8+𝑠4𝐻3𝐻4𝜕𝐻𝜕𝑥𝑖,𝑗2+24𝐻𝑖,𝑗24𝑠8𝑠8+𝑠4𝐻3𝐻4,𝐶12=𝜇𝐻2𝑖+1,𝑗+12𝑏2𝑏1𝐻𝑖+1,𝑗+1𝑏2𝑠1+𝑠3𝜕𝐻𝜕𝑥𝑖+1,𝑗+1𝜕𝐻𝜕𝑦𝑖+1,𝑗+1𝑏1𝑠2+𝑠4,𝐶13=𝜇𝐻2𝑖+1,𝑗12𝑏4𝑏1𝐻𝑖+1,𝑗1𝑏4𝑠1+𝑠3𝜕𝐻𝜕𝑥𝑖+1,𝑗1𝜕𝐻𝜕𝑦𝑖+1,𝑗1𝑏1𝑠2+𝑠4,1RHS=𝐻𝜕𝜏𝑤𝑦𝜏𝜕𝑥𝑤𝑦𝐻2𝜕𝐻𝜕𝑥𝑖,𝑗1𝐻𝜕𝜏𝑤𝑥𝜏𝜕𝑦𝑤𝑥𝐻2𝜕𝐻𝜕𝑦𝑖,𝑗.(19)

Figure 1: Schematic grid diagram.


6. Boundary Conditions

Equation (18) is of fourth order in 𝑥,𝑦, and one must satisfy eight lateral boundary conditions. Thus,  Σ(𝑥,𝑦) is the closed coastal boundary, we can set [15]𝜓=0onΣ.(21) Additionally, we can stipulate no slip condition𝜕𝜓𝜕𝑛=0onΣ.(22)

7. Analytical Verifications

In this section the model is verified versus the analytical solution of Stommel over a rectangular domain with constant depth. The verification for variable depth will be carried against results of Alan and John [3].

7.1. Verifications against Stommel

Stommel retained the bottom friction terms to derive his model. The bottom friction was indeed the dissipative ingredient in the famous paper by Stommel [2] who first recognized the need to introduce dissipation into the system to obtain closed streamlines [15].

Stommel [2] assumed a distribution of wind stress of the form𝜏𝑤𝑥=𝜏𝑜cos𝜋𝑦𝐿,𝜏𝑤𝑦=0,(23) wind stress curl is zero at 𝑦=0 and 𝐿, the southern and northern extents of a closed basin of extent 𝐿 in the latitudinal direction and extent 𝐵 in the longitudinal direction.

Define the small perturbation parameter 𝜀=𝛿/𝑊to represent the ratio of boundary layer thickness to the width of the basin. The solution of the Stommel problem is obtained by asymptotic expansion as [16]𝜓(𝑥,𝑦)=𝜋𝐵𝜏𝑜𝑥𝛽𝐿𝐵1+𝑒𝑎𝑥sin𝜋𝑦𝐿,𝑀𝑦=𝜋𝐵𝜏𝑜1𝛽𝐿𝐵𝑎𝑒𝑎𝑥sin𝜋𝑦𝐿.(24) where 𝑎=𝛽/𝑐. The analytical and numerical solutions are compared for the two values0.025,0.05 for parameter𝜀. Figure 2 shows analytical and numerical solutions of 𝜓 for the value 0.05. Also, Figure 3 compares between the analytical and numerical solutions 𝑀𝑦 for the value0.05. The results are for typical values of𝜏𝑜=0.1Nm2,𝑊=4000km,𝐿=4000km,𝛽=21011m1s1, and 𝑐=2106s1.

Figure 2: Numerical and analytical streamlines in a flat-bottom rectangular ocean for the Stommel problem 𝜀=0.05.
Figure 3: Numerical and analytical 𝑀𝑦 in a flat-bottom rectangular ocean for the Stommel problem 𝜀=0.05.

It is seen that the values of 𝑀𝑦 are not zeros at western and eastern boundaries because the horizontal friction is neglected.

7.2. Verifications against Alan and John [3]

To verify the model over irregular bathymetry we shall considered the simplified model of the Lake Ontario basin used by Lockner [17]. It is an elliptic shoreline with a major axis of 300 km and a minor axis of 87 km, and sloping uniformly everywhere toward a maximum depth of 180 m in the center, Figure 4.

Figure 4: Depth contours for the simple elliptic basin.

The basin is subjected to uniform wind stress opposite to the x-direction. The streamlines from the numerical solution (Figure 5) are perfectly as the results obtained by Alan and John [3]. Although the wind has zero curl, a double-gyre pattern appeared mainly due to depth variation. Depth variation is a main characteristic of real seas.

Figure 5: Streamline contours for the simple elliptic basin.

8. Circulations of the Red Sea

8.1. The Red Sea Wind Distribution

Wind velocity data were obtained for each region of the Red Sea for every month of the year. That is, monthly averaged regional mean wind velocity vectors were used. A monthly mean wind is assumed to flow persistently through the entire month. A regional mean wind is assumed to blow uniformly over the areal extent of the associated region. The Red Sea and its adjacent coasts were simulated by a rectangle 300 km wide and 2040 km long. This rectangle is subdivided into 17 Regions, each has a rectangular form with 300 km × 120 km. regions were numbered from south to north. Table 1 represents monthly mean values of the longitudinal surface wind velocity, while Table 2 represents monthly mean values of the latitudinal surface wind velocity [18].

Table 1: Monthly mean values of the longitudinal surface wind velocity (m/s).
Table 2: Monthly mean values of the latitudinal surface wind velocity (m/s).
8.2. Momentum Transfer at Water Surface

When wind blows over a sea, it exerts a shearing stress on the water surface. Momentum is transferred from the wind to the surface and due to the turbulence transported further down into the water. Accurate knowledge of the force exerted by the wind on the water surface is a prerequisite for every attempt to analyze the movements of water. The formula used for determining the wind shear stress from a given wind speed is [3, 19]𝜏𝑤𝑥=𝜌𝑎𝑐𝑠𝑊𝑥𝑊2𝑥+𝑊2𝑦,𝜏𝑤𝑦=𝜌𝑎𝑐𝑠𝑊𝑦𝑊2𝑥+𝑊2𝑦,(25) where 𝑐𝑠  is wind stress coefficient, 𝑊𝑥,𝑊𝑦, are the component wind speeds in x- and y-directions, respectively, and 𝜌𝑎 is the air density. For ground surfaces the characteristic roughness is truly a characteristic of the surface. The water surface, however, changes character with increasing wind speed and becomes more corrugated. Also, the dimensions of a sea influence the characteristic roughness of the water. Although it may not be possible to truly define a roughness parameter of a sea, this parameter or rather the stress coefficient has been determined from wind measurements and from observations of the wind set-up. Wilson [7] has analyzed the literature on the stress coefficient over oceans and computed average values of the stress coefficient referred to as an elevation of 10 meters:𝑐𝑠=1.5103±0.8103𝑐for𝑊<6m/s,𝑠=2.4103±0.6103for𝑊>6m/s.(26)

From the Great Lakes of North America much lower values have been reported [20]. Observations indicated that an appropriate value for the Great Lakes is 1.2103. Donelan et al. [21] estimated the wind stress from the steady-state water setup of Lake Ontario and found the wind stress coefficient to be 1.3103 for stable and neutral conditions and 1.6103 for unstable conditions.

For enclosed water bodies and fetches of the order 10 km the reported values on the wind stress coefficient are scarce. Bengtsson [22] studied Lake Vomb in southern Sweden. It was found to use 1.2103 for speeds higher than 5.5 m/s and 0.9103 for speeds lower than 4.5 m/s. Alan and John [3] used a value of 2.7103 in their simplified elliptic representation of the Lake Ontario. Shankar et al. [19] used a value of2.55103. In our simulation of the Red Sea a value of 2.0103 is used.

8.3. The Red Sea Bottom Depth Distribution

The Red Sea long narrow basin separates Africa from Asia and extends from NNW to SSE between latitudes 12°30'N to 30° N, in an almost straight line. Its total length is 1932 km and average width is 280 km. In width it decreases from 306 km, near Massawa, to 26 km in the Straits of Bab Al Mandab.

The Red Sea latitude extent from 12° 30'N to 30° N warrants an effective and variable Coriolis force, the 𝛽eect [1]. The central channel reaches great depths of more than 2000 m. However, the average Red Sea depth is 450 m only. This central channel is fringed in the southern Red Sea by shallow areas less than 50 m deep. See Figure 6 that describes the Red Sea basin characteristics.

Figure 6: The Red Sea Basin.

The area of the Red Sea is 4.6×105 km2 and its mean depth is 500 m. The maximum recorded depth is 3039 m in the axial trough at 19° 35'N, 38° 40'E. The Red Sea from the Gulf of Aden lies to the north of Bab Al Mandab near Hannish Island, where the channel is only about 130 m at its deepest.

9. The Red Sea Circulation and Mass Transport

The model can predict the Red Sea circulation driven by different wind distributions. The coefficient of bottom friction is taken as 0.0025. The output presented below included both the stream function and mass transport vectors. Simulations are carried for different cases, averaged values of summer, winter and also for of the entire year. Figures 7, 8, and 9 show the seasonally-averaged stream lines and mass transport vectors for winter, summer, and the whole year respectively. Each figure is composed of three distributions: (a) The driving wind stress distribution, (b) The resulting field of mass-transport vectors, and (c) The resulting pattern of circulation streamlines.

Figure 7: Wind stress vectors (a) mass transport vectors (b) and Streamlines (c) Values in the Red Sea are winter-averaged.
Figure 8: Wind stress vectors (a) mass transport vectors (b) and Streamlines (c) Values in the Red Sea are summer-averaged.
Figure 9: Wind stress vectors (a) mass transport vectors (b) and Streamlines (c) Values in the Red Sea are year-averaged.

These plots of seasonal Red Sea mass-transport vectors and streamlines demonstrate the manner by which they are influenced by the interaction of two patterns: the wind-stress pattern and the bottom-bathymetry pattern. The wind-driven circulation show some features which can be summarized bellow: (1)The largest mass transport in all cases is in the deepest area in the Red sea. The central area is in the northern part. (2)Mass transport direction is the same as wind direction along the central part and then circulates in the opposite direction beside the shoreline. (3)The mass transport values in the northern part are much greater than the values in the southern part. Complex depth variations in the northern part cause this difference. (4)The wind in the summer blows in the south-west direction over the total area but in the winter it blows in the south-west direction in the northern part and in the north-east direction in the southern part. As a result of this variation flow patterns are more complex in the winter than in the summer. (5)The summer has a maximum value of streamline of 23 Sverdrup while the winter has a maximum value of 30Sverdrup while the maximum value for the averaged-year simulation is24Sverdrup.


  1. M. A. F. M. Zaki, “Combined wind and density-driven circulation in enclosed seas with variable coriolis parameter,” Journal of Engineering and Applied Sciences, vol. 42, pp. 433–446, 1995. View at Google Scholar
  2. H. Stommel, “The westward intensification of wind-diven ocean currents,” American Geophysical Union, vol. 29, no. 2, 1948. View at Google Scholar
  3. J. W. Alan and H. T. John, “Steady wind driven currents in a large lake with depth-dependent eddy viscosity,” Journal of Physical Oceanography, vol. 6, no. 1, 1976. View at Google Scholar · View at Scopus
  4. W. Lick, “Numerical modelling of lake currents,” Annual Review of Fluid Mechanics, vol. 4, pp. 49–74, 1976. View at Google Scholar · View at Scopus
  5. L. H. Kantha and C. A. Clayson, Numerical Models of Oceans and Oceanic Processes, Academic Press, 2000.
  6. G. Neumann and W. J. Pierson, Principles of Physical Oceanography, Prentice-Hall, New Jersey, NJ, USA, 1966.
  7. B. W. Wilson, “Note on surface wind stress over water at low and high speeds,” Journal of Geophysical Research, vol. 65, no. 10, pp. 3377–3382, 1960. View at Google Scholar
  8. P. D. Thompson, Numerical Weather Analysis and Prediction, Macmillan, New York, NY, USA, 1961.
  9. J. O. Backhaus and E. Maier-Reimer, “On seasonal circulation in the North Sea,” in North Sea Dynamics, J. Sundermann and W. Lenz, Eds., p. 693, Springer, Heidelberg, Germany, 1983. View at Google Scholar
  10. G. K. Furnes, “Wind effects in the North Sea,” Journal of Physical Oceanography, vol. 10, no. 6, pp. 978–984, 1980. View at Google Scholar · View at Scopus
  11. D. B. Rao and T. S. Murty, “Calculation of the steady state wind-driven circulations in Lake Ontario,” Archiv für Meteorologie, Geophysik und Bioklimatologie Serie A, vol. 19, no. 2, pp. 195–210, 1970. View at Publisher · View at Google Scholar · View at Scopus
  12. A. R. Mitchell and D. F. Griffths, Computational Methods in Partial Differential Equations, John Wiley & Sons, New York, NY, USA, 1980.
  13. S. C. Chapra and R. P. Canale, Numerical Methods for Engineers, McGraw-Hill, 5th edition, 2006.
  14. A. K. Singh and B. S. Bhadauria, “Finite difference formulae for unequal sub-intervals using lagrange’s interpolation formula,” International Journal of Mathematics and Analysis, vol. 3, no. 17, pp. 815–827, 2009. View at Google Scholar
  15. G. L. Mellor, Introduction to Physical Oceanography, Springer, New York, NY, USA, 1996.
  16. M. Van Dyke, Perturbation Methods in Fluid Mechanics, Academic Press, New York, NY, USA, 1964.
  17. D. Lockner, “Sensitivity of a numerical circulation model of Lake Ontario to change in Lake Symmetry and friction depth and to variable wind stress,” IFYGL. Rochester Embayment Project Report No. 2, Department of Geological Sciences, University of Rochester, 1973. View at Google Scholar
  18. M. A. F. M. Zaki and I. E. Mobarek, “A model of the Red Sea air pollution II- the computational model,” Scientific Engineering Bulletin, Faculty of Engineering, Cairo University, vol. 2, 1981. View at Google Scholar
  19. N. J. Shankar, X. B. Chao, H. F. Cheong, and Y. Q. Li, “Numerical simulation of oil spills in the Singapore coastal waters,” in Recent Advances in Marine Science and Technology, N. Saxena, Ed., pp. 177–188, PACON International, 2001. View at Google Scholar
  20. J. P. Bruce, D. V. Andersson, and G. K. Rodgers, “Temperature, humidity and wind profiles over the great lakes,” in Proceedings of the 4th Conference on Great Lakes Research, pp. 65–70, 1961.
  21. M. A. Donelan, F. C. Elder, and P. F. Hamblin, “Determination of the aerodynamic drag coefficient from wind set-up,” in Proceedings of the 17th Conference on Great Lakes Research, pp. 778–788, 1974.
  22. L. Bengtsson, Wind Stress on Small Lakes, Tekniska Högskolan, Lund, Sweden, 1973.