# A Barotropic Model of the Red Sea Circulation

Authors: Mohamed M. El-Shabrawy; Khaled M. Fasseih; Mamdouh A. Zaki
Journal: ISRN Civil Engineering
 (2012)

---

## Abstract

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.


---

## Body

## 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, 3–5], 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
(1)∂ρ∂t+∂(ρu)∂x+∂(ρv)∂y+∂(ρv)∂z=0,
where ρ is the fluid density, u,v,w are velocity components in x,y,z directions, and t is the time. Also, the momentum equation [6] in the Cartesian form is(2a)∂u∂t+u∂u∂x+v∂u∂y+w∂u∂z=-1ρ∂p∂x+fv-2ωcos⁡⁡φw+μv∂2u∂z2+μh(∂2u∂x2+∂2u∂y2),(2b)∂v∂t+u∂v∂x+v∂v∂y+w∂v∂z=-1ρ∂p∂y+fu+μv∂2v∂z2+μh(∂2v∂x2+∂2v∂y2),(2c)∂w∂t+u∂w∂x+v∂w∂y+w∂w∂z=-1ρ∂p∂z-g-2ωcos⁡⁡φu+μv∂2w∂z2+μh(∂2w∂x2+∂2w∂y2),where p is the pressure field, x,y,z are Cartesian coordinate system fixed with respect to the rotating earth, with z upward, and f is the Coriolis parameter. μh,μv 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”(3)∂∂t=0.(2)
Turbulent shear forces in the momentum equation ((2a), (2b), and (2c)) are simulated by depth averaged eddy (momentum exchange) coefficients μh,μv 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],(4)∂u∂x+∂v∂y+∂w∂z=0,(5)-fv=-1ρ∂p∂x+μv∂2u∂z2+μh(∂2u∂x2+∂2u∂y2),(6)fu=-1ρ∂p∂y+μv∂2v∂z2+μh(∂2v∂x2+∂2v∂y2),(7)∂p∂z=-ρg.
This set of linear partial equations in u, v, w and p are subject to boundary conditions; a wind stress at the mean sea surface, z=0, and a frictional stress at the sea bottom, z=-H(8)(τxw,τyw)=μv(∂u∂z,∂v∂z)z=0,(τxb,τyb)=μv(∂u∂z,∂v∂z)z=-H.
Define the components of mass transport between the sea surface,ζ, and the bottom, -H, in the x- and y-directions as(9)Mx=∫-Hζudz,My=∫-Hζvdz.
Integrating equation (4) over z between the boundaries -H(x,y) and ζ(x,y), and using the conditions of integrity of the free surface and impermeability of the bottom [5, 6],(10)∂Mx∂x+∂My∂y=0.
A mass transport stream function ψ(x,y)can then be defined(11)Mx=∂ψ∂y,My=-∂ψ∂x.
Integration of (7) over z between the boundaries -H(x,y) and ζ(x,y) gives(12)p=pa+ρg(ζ+H),
where pais the atmospheric pressure.Integration of (5) and (6) over z between the boundaries -H(x,y) and ζ(x,y) gives(13)-fMy=-1ρ∂∂x∫-Hζpdz+μv((∂u∂z)ζ-(∂u∂z)-H)+μh(∂2Mx∂x2+∂2Mx∂y2),(14)fMx=-1ρ∂∂y∫-Hζpdz+μv((∂v∂z)ζ-(∂v∂z)-H)+μh(∂2My∂x2+∂2My∂y2).
Dividing both equations by H  differentiating (13) for y and (14) for x, and subtracting to eliminate pressure from the governing equations,(15)∂∂x(fH)Mx+∂∂y(fH)My=(-∂∂y(τxwH)+∂∂y(τxbH)+∂∂x(τywH)-∂∂x(τybH))+μh(∂∂x(1H∂2My∂x2+1H∂2My∂y2)-∂∂y(1H∂2Mx∂x2+1H∂2Mx∂y2)).
It can be easily simplified to be(16)-μhH2[H(∂4ψ∂x4+2∂4ψ∂x2∂y2+∂4ψ∂y4)-∂H∂x(∂3ψ∂x3+∂3ψ∂x∂y2)-∂H∂y(∂3ψ∂y3+∂3ψ∂y∂x2)]+fH2∂H∂x∂ψ∂y-(βH-fH2∂H∂y)∂ψ∂x=[(1H∂τyw∂x-τywH2∂H∂x)-(1H∂τxw∂y-τxwH2∂H∂y)+(1H∂τyb∂x-τybH2∂H∂x)-(1H∂τxb∂y-τxbH2∂H∂y)].
Following the common assumption of the turbulent shear stress as linear function of horizontal mass transport [6, 7], usingτxb=cMxandτyb=cMy in (16) leads to(17)-μhH2[H(∂4ψ∂x4+2∂4ψ∂x2∂y2+∂4ψ∂y4)-∂H∂x(∂3ψ∂x3+∂3ψ∂x∂y2)-∂H∂y(∂3ψ∂y3+∂3ψ∂y∂x2)]+cH(∂2ψ∂x2+∂2ψ∂y2)+(fH2∂H∂x-cH2∂H∂y)∂ψ∂y-(βH-fH2∂H∂y-cH2∂H∂x)∂ψ∂x=(1H∂τyw∂x-τywH2∂H∂x)-(1H∂τxw∂y-τxwH2∂H∂y),
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 Figure1. Expanding the spatial derivatives using the methods of [12, 14] of (17) will lead to a 13-point finite difference equation ofψ: (18)C1ψi-2,j+C2ψi-1,j+C3ψi,j+C4ψi+1,j+C5ψi+2,j+C6ψi-1,j+1+C7ψi-1,j-1+C8ψi,j+1+C9ψi,j-1+C10ψi,j+2+C11ψi,j-2+C12ψi+1,j+1+C13ψi+1,j-1=RHS,
where(19)C1=-μhHi-2,j2((-6(s3-2s1-s5)h3s7(s7+s3)H1H2)(∂H∂x)i-2,j+24Hi-2,jh4s7(s7+s3)H1H2),C2=-μhHi-1,j2((6(s7+s3-2s1-s5)h3s7s3(s7+s3)(H2-s7)+2(b2+b4)h(s1+s3))(∂H∂x)i-1,j-Hi-1,j(24h4s7s3(s7+s3)(H2-s7)+4b3(b2+b4))+cb3Hi-1,j+1h(s1+s3)(βHi-1,j-fHi,j2(∂H∂y)i-1,j-cHi-1,j2(∂H∂x)i-1,j)),C3=-μhHi,j2(Hi,j(24h4s2s4(s8+s4)(s2+s6)+2h4s1s3(s7+s3)(s1+s5)+(b2+b4)(b1+b3))-(∂H∂x)i,j×(6(s7+2s3-2s1-s5)h3s1s3(s7+s3)(s1+s5))-(∂H∂y)i,j(6(s8+2s4-2s2-s6)h3s2s4(s8+s4)(s2+s6)))-2cHi,j(b1+b3+b2+b4),C4=-μhHi+1,j2(∂H∂x(∂H∂x)i+1,j(-6(s7+2s3-s1-s5)h3s1s5(s1+s3)H1-2(b2+b4)h(s1+s3))(∂H∂x)i+1,j-Hi+1,j(24h4s1s5(s1+s3)H1+4b1(b2+b4))+cb1Hi+1,j-1h(s1+s3)(βHi+1,j-fHi,j2(∂H∂y)i+1,j-cHi+1,j2(∂H∂x)i+1,j)),C5=-μhHi+2,j2((6(s7+2s3-s1)h3H2(H2-s7)(s1+s5)s5)(∂H∂x)i+2,j+24Hi+2,jh4H2(H2-s7)(s1+s5)s5),C6=-μhHi-1,j+12(2b2b3Hi-1,j+1-(b2h(s1+s3))(∂H∂x)i-1,j+1-(∂H∂y)i-1,j+1(b3h(s2+s4))),C7=-μhHi-1,j-12(2b4b3Hi-1,j-1-(b4h(s1+s3))(∂H∂x)i-1,j-1-(∂H∂y)i-1,j-1(b3h(s2+s4))),C8=-μhHi,j+12(Hi,j+1(-24h4s2s6(s2+s4)H1-4b2(b1+b3))+(∂H∂y)i,j+1(6(s8+2s4-s2-s6)h3s2s6(s2+s4)H3+2(b1+b3)h(s2+s4)))+2cb2Hi,j+1+(fHi,j+12(∂H∂x)i,j+1-cHi,j+12(∂H∂y)i,j+1)(1h(s2+s4)),C9=-μhHi,j-12(Hi,j-1(-24h4s8s4(s8+s4)(H4-s8)-4b4(b1+b3))+(∂H∂y)i,j-1(6(s8+s4-2s2-s6)h3s8s4(s8+s4)(H4-s8)-2(b1+b3)h(s2+s4)))+2cb4Hi,j-1-(fHi,j-12(∂H∂x)i,j-1-cHi,j-12(∂H∂y)i,j-1)(1h(s2+s4)),C10=-μhHi,j+22((-6(s8+2s4-s2)h3H4(H4-s8)(s2+s6)s6)(∂H∂x)i,j+2+24Hi,j+2h4H2(H2-s7)(s1+s5)s5),C11=-μhHi,j-22((-6(s4-2s2-s6)h3s8(s8+s4)H3H4)(∂H∂x)i,j-2+24Hi,j-2h4s8(s8+s4)H3H4),C12=-μhHi+1,j+12(2b2b1Hi+1,j+1-(b2h(s1+s3))(∂H∂x)i+1,j+1-(∂H∂y)i+1,j+1(b1h(s2+s4))),C13=-μhHi+1,j-12(2b4b1Hi+1,j-1-(b4h(s1+s3))(∂H∂x)i+1,j-1-(∂H∂y)i+1,j-1(b1h(s2+s4))),RHS=(1H∂τyw∂x-τywH2∂H∂x)i,j-(1H∂τxw∂y-τxwH2∂H∂y)i,j.Figure 1
Schematic grid diagram.Where(20)b1=2s1(s1+s3)b2=2s2(s2+s4)b3=2s3(s1+s3)b4=2s4(s2+s4)b0=b1+b2+b3+b4.


## 6. Boundary Conditions

Equation (18) is of fourth order in x,y, and one must satisfy eight lateral boundary conditions. Thus,  Σ(x,y)  is the closed coastal boundary, we can set [15](21)ψ=0on⁡Σ.
Additionally, we can stipulate no slip condition(22)∂ψ∂n=0on⁡Σ.


## 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(23)τxw=-τocos⁡⁡πyL,τy=0w,
wind stress curl is zero at y=0 and L, the southern and northern extents of a closed basin of extent L in the latitudinal direction and extent B in the longitudinal direction.Define the small perturbation parameterε=δ/Wto 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](24)ψ(x,y)=πBτoβL(xB-1+e-ax)sin⁡(πyL),My=-πBτoβL(1B-ae-ax)sin⁡(πyL).
where a=β/c. 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 My for the value0.05. The results are for typical values ofτo=0.1Nm-2,W=4000km,L=4000km,β=2*10-11m-1s-1, and c=2*10-6s-1.Numerical and analytical streamlines in a flat-bottom rectangular ocean for the Stommel problemε=0.05.
(a)(b)Figure 3
Numerical and analyticalMy in a flat-bottom rectangular ocean for the Stommel problem ε=0.05.It is seen that the values ofMy 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 thex-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.


## 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(23)τxw=-τocos⁡⁡πyL,τy=0w,
wind stress curl is zero at y=0 and L, the southern and northern extents of a closed basin of extent L in the latitudinal direction and extent B in the longitudinal direction.Define the small perturbation parameterε=δ/Wto 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](24)ψ(x,y)=πBτoβL(xB-1+e-ax)sin⁡(πyL),My=-πBτoβL(1B-ae-ax)sin⁡(πyL).
where a=β/c. 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 My for the value0.05. The results are for typical values ofτo=0.1Nm-2,W=4000km,L=4000km,β=2*10-11m-1s-1, and c=2*10-6s-1.Numerical and analytical streamlines in a flat-bottom rectangular ocean for the Stommel problemε=0.05.
(a)(b)Figure 3
Numerical and analyticalMy in a flat-bottom rectangular ocean for the Stommel problem ε=0.05.It is seen that the values ofMy 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 thex-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. Table1 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).
Jan.Feb.Mar.Apr.MayJun.Jul.Aug.Sept.oct.Nov.Dec.14.44.42.42.42.4−0.1−2.4−2.4−0.94.44.44.424.44.44.42.42.4−0.1−2.4−2.4−0.94.44.46.736.76.76.74.42.4−0.9−2.4−2.4−0.94.49.39.346.74.44.42.40.9−0.9−0.9−0.9−0.94.46.76.754.42.42.40.90−2.4−0.9−0.9−0.92.44.44.462.40.90.90.10.1−2.4−2.4−0.9−0.90.94.42.470.9000.10.9−2.4−2.4−0.4−0.902.42.48−0.9−0.9−0.9−0.9−0.9−4.4−2.4−2.4−2.40−0.9−0.99−0.9−2.4−0.9−0.9−2.4−4.4−2.4−2.4−2.40−0.9−0.910−2.4−2.4−2.4−2.4−2.4−4.4−2.4−2.4−2.4−2.4−0.9−2.411−2.4−4.4−2.4−2.4−2.4−4.4−2.4−2.4−2.42.4−2.4−2.412−2.4−4.4−2.4−2.4−2.4−4.4−2.4−2.4−4.4−2.4−2.4−2.413−4.4−4.4−4−2.4−2.4−4.4−2.4−2.4−4.4−2.4−2.4−2.414−4.4−4.4−4.4−2.4−4.4−4.4−2.4−2.4−4.4−4.4−4.4−2.415−4.4−4.4−4.4−4.4−4.4−4.4−4.4−4.4−6.7−4.4−4.4−4.416−4.4−4.4−4.4−4.4−4.4−6.7−6.7−6.7−9.3−6.7−4.4−4.417−4.4−4.4−4.4−4.4−6.7−6.7−6.7−6.7−6.7−6.7−4.4−4.4Table 2
Monthly mean values of the latitudinal surface wind velocity (m/s).
Jan.Feb.Mar.Apr.MayJun.Jul.Aug.Sept.oct.Nov.Dec.1−2.9−2.9−1.0310.600.60.4−0.15−1.76−4.4−4.42−2.9−2.9−1.882.10.600.60.4−0.15−1.76−4.4−6.73−4.4−4.4−2.882.10.60.20.60.4−0.15−1.76−9.3−9.34−4.4−2.9−1.910.20.20.20.150.15−1.76−6.7−6.75−2.9−1.6−1.030.400.60.150.150.150.96−4.4−4.46−1.60.6−0.4000.60.40.150.150.36−4.4−2.47−0.60000.20.60.40.40.150−2.4−0.980.60.60.40.40.21.10.40.40.40−0.9−2.490.61.60.40.40.61.10.40.40.40−0.9−0.9101.61.6110.61.10.40.40.40.960.92.4111.61.6110.61.10.60.40.70.962.42.4121.62.9110.61.10.60.60.70.962.42.4132.92.91.910.61.10.60.40.70.962.42.4142.92.91.911.11.10.60.40.71.764.42.4152.92.91.91.91.11.70.70.71.11.764.44.4162.92.91.91.91.11.71.11.11.552.74.44.4172.92.92.12.11.71.71.11.11.12.64.44.4


### 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](25)τxw=ρacsWxWx2+Wy2,τyw=ρacsWyWx2+Wy2,
where cs  is wind stress coefficient, Wx,Wy, are the component wind speeds in x- and y-directions, respectively, and ρa 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:(26)cs=1.5*10-3±0.8*10-3forW<6m/s,cs=2.4*10-3±0.6*10-3forW>6m/s.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.2*10-3. 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.3*10-3 for stable and neutral conditions and 1.6*10-3 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.2*10-3 for speeds higher than 5.5 m/s and 0.9*10-3 for speeds lower than 4.5 m/s. Alan and John [3] used a value of 2.7*10-3 in their simplified elliptic representation of the Lake Ontario. Shankar et al. [19] used a value of2.55*10-3. In our simulation of the Red Sea a value of 2.0*10-3 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 β-effect  [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.


## 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. Table1 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).
Jan.Feb.Mar.Apr.MayJun.Jul.Aug.Sept.oct.Nov.Dec.14.44.42.42.42.4−0.1−2.4−2.4−0.94.44.44.424.44.44.42.42.4−0.1−2.4−2.4−0.94.44.46.736.76.76.74.42.4−0.9−2.4−2.4−0.94.49.39.346.74.44.42.40.9−0.9−0.9−0.9−0.94.46.76.754.42.42.40.90−2.4−0.9−0.9−0.92.44.44.462.40.90.90.10.1−2.4−2.4−0.9−0.90.94.42.470.9000.10.9−2.4−2.4−0.4−0.902.42.48−0.9−0.9−0.9−0.9−0.9−4.4−2.4−2.4−2.40−0.9−0.99−0.9−2.4−0.9−0.9−2.4−4.4−2.4−2.4−2.40−0.9−0.910−2.4−2.4−2.4−2.4−2.4−4.4−2.4−2.4−2.4−2.4−0.9−2.411−2.4−4.4−2.4−2.4−2.4−4.4−2.4−2.4−2.42.4−2.4−2.412−2.4−4.4−2.4−2.4−2.4−4.4−2.4−2.4−4.4−2.4−2.4−2.413−4.4−4.4−4−2.4−2.4−4.4−2.4−2.4−4.4−2.4−2.4−2.414−4.4−4.4−4.4−2.4−4.4−4.4−2.4−2.4−4.4−4.4−4.4−2.415−4.4−4.4−4.4−4.4−4.4−4.4−4.4−4.4−6.7−4.4−4.4−4.416−4.4−4.4−4.4−4.4−4.4−6.7−6.7−6.7−9.3−6.7−4.4−4.417−4.4−4.4−4.4−4.4−6.7−6.7−6.7−6.7−6.7−6.7−4.4−4.4Table 2
Monthly mean values of the latitudinal surface wind velocity (m/s).
Jan.Feb.Mar.Apr.MayJun.Jul.Aug.Sept.oct.Nov.Dec.1−2.9−2.9−1.0310.600.60.4−0.15−1.76−4.4−4.42−2.9−2.9−1.882.10.600.60.4−0.15−1.76−4.4−6.73−4.4−4.4−2.882.10.60.20.60.4−0.15−1.76−9.3−9.34−4.4−2.9−1.910.20.20.20.150.15−1.76−6.7−6.75−2.9−1.6−1.030.400.60.150.150.150.96−4.4−4.46−1.60.6−0.4000.60.40.150.150.36−4.4−2.47−0.60000.20.60.40.40.150−2.4−0.980.60.60.40.40.21.10.40.40.40−0.9−2.490.61.60.40.40.61.10.40.40.40−0.9−0.9101.61.6110.61.10.40.40.40.960.92.4111.61.6110.61.10.60.40.70.962.42.4121.62.9110.61.10.60.60.70.962.42.4132.92.91.910.61.10.60.40.70.962.42.4142.92.91.911.11.10.60.40.71.764.42.4152.92.91.91.91.11.70.70.71.11.764.44.4162.92.91.91.91.11.71.11.11.552.74.44.4172.92.92.12.11.71.71.11.11.12.64.44.4


## 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](25)τxw=ρacsWxWx2+Wy2,τyw=ρacsWyWx2+Wy2,
where cs  is wind stress coefficient, Wx,Wy, are the component wind speeds in x- and y-directions, respectively, and ρa 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:(26)cs=1.5*10-3±0.8*10-3forW<6m/s,cs=2.4*10-3±0.6*10-3forW>6m/s.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.2*10-3. 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.3*10-3 for stable and neutral conditions and 1.6*10-3 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.2*10-3 for speeds higher than 5.5 m/s and 0.9*10-3 for speeds lower than 4.5 m/s. Alan and John [3] used a value of 2.7*10-3 in their simplified elliptic representation of the Lake Ontario. Shankar et al. [19] used a value of2.55*10-3. In our simulation of the Red Sea a value of 2.0*10-3 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 β-effect  [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. Figures7, 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.Wind stress vectors (a) mass transport vectors (b) and Streamlines (c) Values in the Red Sea are winter-averaged.
(a)(b)(c)Wind stress vectors (a) mass transport vectors (b) and Streamlines (c) Values in the Red Sea are summer-averaged.
(a)(b)(c)Wind stress vectors (a) mass transport vectors (b) and Streamlines (c) Values in the Red Sea are year-averaged.
(a)(b)(c)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 23Sverdrup while the winter has a maximum value of 30Sverdrup while the maximum value for the averaged-year simulation is24Sverdrup.



---

*Source: 345979-2012-06-14.xml*