^{1}

^{*}

^{1}

^{1}

^{1}

The provision of economic resources to countries at the shore of Gulf of Guinea is a very important effort that has helped West Africa to develop to some extent. Taking the study area as a system it will be important to know the dynamics that occur in it to help neighouring countries predict its physical and thermodynamics states at all times. Ghana is located in West Africa and it is bordered in the south by the Gulf of Guinea or Atlantic Ocean. The objective of this research is to characterize the hydrodynamic circulation of the Gulf of Guinea areas neighouring Ghana. A 3-Dimensional hydrodynamic model was implemented in the territorial waters of Ghana using the Modelo Hidrodinamico, MOHID model to represent the dynamics and study the complex circulation pattern of the sea. To achieve this in an affordable computational time, nested domain approach was used to implement the hydrodynamic model in both 2 and 3-Dimensional gridded levels. The first level is a barotropic model with only tide. The nested domains of the rest of the levels are baroclinic forced with atmospheric and oceanic elements. To quantify its accuracy, the model was validated and calibrated in three stages; first, the frequency of the water level, followed by the circulation pattern and last, analysing the nature and profiles of the atmospheric and oceanic elements. The implemented model showed good agreement with the measured water surface level in the domain with mean error values not exceeding 14.00% of the measured data and with correlation factors higher than 0.80. Also, the intensity and direction of velocity observed in the current data are well represented by the model at the water surface levels with mean errors lower than 20.00% of the measured data components. The profiles obtained for both the temperature and salinity at shore show completely a straight line for salinity and also a straight line which is slightly curved at the top for the temperature profile. The vertical straight line for the salinity profile shows that at the shore the waters are not stratified vertically or they are well mixed. The slight curve at the top of the temperature graph accounts for the mixing dynamics that occur close to land or the effect of heat flux at the surface. It can be concluded that the hydrodynamic model obtained by this study is the true reflection of the territorial waters of Ghana.

The provision of economic resources to countries at the shore of Gulf of Guinea makes it very important to countries in West Africa [

Ghana is located in West Africa. It is bordered with Ivory Coast at the west, north by Burkina Faso, east by Togo, and south by the Gulf of Guinea or Atlantic Ocean. It is located between longitude 3.25˚W to 1.19˚E and latitude 4.74˚N to 11.18˚N. Accra is the capital city of Ghana situated in the Greater Accra Region and it stretches along the coast. Apart from Accra, there are other major cities like Takoradi, Sekondi, Cape Coast, Elmina, Keta, Tema, etc. also along the coast and the main occupation of the inhabitants is fishing. The coast is also filled with natural, sensitive and strategic resources like forts, castles, industries, water bodies (i.e. Volta river estuary and lagoon) and tourist sites. The discovery of commercial quantities of oil has placed all these resources under threat. The economic importance of the coast including ports and ports activities, high traffic of vessel activities, exploration exercises, local industries, fisheries and tourism justify the need for studying the hydrodynamics of the area [

The objective of this research is to study the hydrodynamic circulation of the Gulf of Guinea areas neighouring Ghana. Another goal is to validate the MOHID water system model using data available to proof it capability to hydrodynamic modeling.

In this research, Modelo Hidrodinâmico, MOHID model, a 3-Dimensional water modeling system was used [

In order to simulate all possible dynamics in an affordable computer processing unit (CPU) time, nested domain approach was used to implement the hydrodynamic model in both 2 and 3-Dimensional gridded levels. It is made up of three domains; Domains 1, 2, and 3 are contained in levels 1, 2, and 3, respectively. The first level is a barotropic model with only tide. The nested domains of the rest of the levels are baroclinic.

At the open boundaries, the main domain was forced with tide from the FES (Finite Element Solution), 2012 tide model. The domains of the second and third of the second and third levels are subdomains of the first level domain. Open boundaries come in the sense of surface and lateral boundaries with cell volume not sided by land. At these levels (2^{nd} and 3^{rd}), they are purely 3-Dimensional in grid. The free surface was forced with atmospheric conditions; Pressure from the Atmospheric, Wind velocity, Solar radiation, Air temperature, Relative Humidity, downward Long Wave Radiation and Cloud cover. The atmospheric data was taken from The Global Forecast System (GFS) which is a weather model.

A database containing the following average daily distribution information parameters was forced at the open boundary or lateral boundary (i.e. Temperature, Salinity, Velocity, and Water level). The forcing was done for both domains in the second and third levels for the horizontal grids with a vertical resolution of 43 layers. Closed boundaries were taken to be land cell where computation will not be performed.

The results were analysed statistically to investigate the margin of deviation from measured data.

This section introduces the equations and assumptions involved in the implementation of the hydrodynamic model of the coastal waters of Ghana.

Hydrodynamics; the model solves the 3-Dimensional Navier Stokes incompressible primitive equations. Hydrostatic equilibrium is accounted for and also Boussinesq and Reynolds approximations as;

∂ ∂ t ∫ V v H d V = A + B + C + D (1)

where

A = − ∮ A v H ( v ⋅ n ) d A

B = ∮ A v H ( v ⋅ n ) d A

C = − 1 ρ ∮ A p ⋅ n H d A

D = ∮ V 2 Ω ∗ v H d V + F

The model is made to allow all these assumed approximations in its solution.

Considering momentum balance equations, a Cartesian scale for horizontal velocities is given as:

∂ t u = − ∂ y ( u v ) − ∂ z ( u w ) − ∂ x ( u u ) + f v − 1 / ρ o ( ∂ x P ) + ( ∂ x ( v H + v ) ∂ x u ) + ( ∂ y ( v H + v ) ∂ y u ) + ( ∂ z ( v t + v ) ∂ z u ) (2)

∂ t v = − ∂ y ( v v ) − ∂ z ( u w ) − ∂ x ( v u ) − f u − 1 / ρ o ( ∂ y P ) + ( ∂ x ( v + v H ) ∂ x v ) + ( ∂ y ( v + v H ) ∂ y v ) + ( ∂ z ( v t + v ) ∂ z v ) (3)

where, u, v and w are the velocity vector components in the Cartesian coordinate that is the x, y and z directions respectively, f the Coriolis parameter, ν_{H} and ν_{t} the viscosities for turbulent in the various directions that is horizontal and vertical directions, pressure is designated by p. The first three terms on the right stand for the changing velocities that is advective transport, the fourth term is Coriolis force, and the three terms following are the pressure gradient and turbulent diffusion (last three terms).

The mass balance equation (that is, the incompressible continuity equation) was used to compute the vertical velocity, w by integrating the depth, z and bottom, -h;

∂ x u + ∂ y v + ∂ z w = 0 (4)

The Reynolds and Boussinesq approximations were assumed in the model. Also, Hydrostatic equilibrium was assumed.

By using this expression and the Boussinesq approximation, the horizontal pressure gradient in the direction x_{i} can be divided into three contributions;

∂ x i P = ∂ x i P a t m − g ρ 0 ∂ x i η − g ∫ z η ∂ x i ρ d z (5)

The density is estimated from the transported salinity and temperature within the domain by the water properties module. This estimation is done from the equation of state [

ρ = ( 5890 − 0.375 T 2 + 38 T + 3 S ) / ( ( 1779.5 − 0.00745 T 2 + 11.25 T ) − ( 3.8 + 0.01 T ) S 2 + 0.698 ( 5890 − 0.375 T 2 + 38 T + 3 S ) ) (6)

Discretization; the Module for Geometry computes the boundary areas as well as the volume of the domain considering the bathymetry data and the free surface elevation. The model uses a finite volume approach [

In the model, volume of an element is represented in a way that allows a degree of freedom in the vertical direction and in the horizontal grid, a Cartesian orthogonal is allowed which is staggered in an Arakawa C [

The ADI scheme calculates the horizontal velocities alternatively both explicitly and implicitly; that calculates one component of the horizontal velocity explicitly while the other is calculated implicitly, making sure that the internal and external modes have the same time steps [

The model has the features of allowing five different boundary conditions; free surface, moving boundary, lateral closed or opened boundary and bottom boundary.

Free surface; a zero vertical advective fluxes over the surface is assumed. A flux at the surface (i.e. W as flux) is:

W flux | surface = 0 (7)

By means of the wind surface stress, τ_{w}, the Diffusive flux of momentum is explicitly imposed:

∂ v ¯ H / ∂ z | surface = τ w (8)

Wind stress is estimated from a quadratic friction law:

τ w = C D ρ a W | W | (9)

where drag coefficient, C_{D}, depends on the wind speed. W and ρ_{a} are the measurements at a height of 10.00 m over the free sea surface for wind speed and air density respectively.

Bottom boundary; the bottom advective fluxes are imposed as null and diffusive flux of momentum is estimated by means of a bottom stress that is calculated by a non-slip method with a quadratic law that depends on the near-bottom velocity. So, the diffusive term at the bottom is written as:

v ( ∂ v ¯ H / ∂ z ) | bottom = C D v H | v H | (10)

The drag coefficient, C_{D}, is estimated from the equation:

C D = ( K / log ( ( z + z 0 b ) / z 0 b ) ) 2 (11)

where κ and z 0 b are the von Karman constant and the bottom roughness length respectively. The equation shows the characteristic of the law of the wall near boundaries especially at the bottom [

Lateral closed boundaries; computation is limited by land. At this boundary the cells are impermeable. At the cell faces, surfacing land or in contact with land, a zero water and momentum diffusion fluxes are specified. During the implementation of this closed boundary a free slip condition is assumed;

∂ v ¯ H / ∂ η = 0 (12)

v ⋅ n = 0 (13)

Open boundaries; with the concept of nesting domains, localized areas like lagoons, lakes, estuaries, etc. can be studied. When a space within a region is extracted for a study, it should be implemented in such a way that guarantees values of the properties for both the region and the environment are allowed to enter and leave each other. The opening through which these properties are allowed between these two bodies or the surface/boundary through which they communicate is called open boundary [

Moving boundaries; when the position of the boundary changes with time it is referred to as moving boundary. Moving boundaries are closed boundaries in nature and change its position with time. They often exist at inter tidal zones where the points at the boundary are covered or uncovered by water with time within the domain. A stable algorithm is implemented within MOHID for these regions not to affect the hydrodynamics estimations, for example, at estuaries [

^{rd} domain is the smallest with higher resolution compared to the first and second levels containing the 1^{st} and 2^{nd} domains respectively.

The minimum and maximum depths for the first domain are 10 and 433 m respectively, stretching from longitude 0.2˚N to 6.0˚N between latitude 4.95˚E and 4.75˚W. It also has a 2-Dimensional geometry with one sigma vertical layer grid extending from the surface to the maximum depth of the bathymetry data (that is, 433 km). Validation of simulated results was carried out at the following locations; Tema, Accra, Cape Coast and Takoradi at the third level.

The hydrodynamic pattern for the region is also shown in Figures 6(a)-(d) for particular time during the simulation period.

At each stage in the implementation, the model was validated. An overall validation was carried out at the end of the hydrodynamic model (i.e. the third level) looking at the water level, water currents and the pattern at the surface, salinity and temperature profiles.

From the graphs (Figures 3-6) the model simulations predicted very well and were very close to the measured data especially the tide frequency (

The water current was also found to correlate with the measured data at a factor of 20.00%. The pattern was in the same sense as the measured.

The temperature and salinity were analysed looking at their nature as the country lack weather station that keep such ocean elements. During the simulation also there was no algo float located within the domain. Both the temperature and salinity showed similar graph (i.e. almost straight line). The kind of line that was obtained by the salinity can be explained by the fact that the waters at the shore are well mixed. A slight difference was observed for the temperature at the top of the curve (i.e. slight curving at the top). The slight curve can be attributed to the fact that there is a mixing of hot fresh water from land or by the presence of heat flux at the surface of the water.

The 3-Dimensional Hydrodynamic model was used to simulate the dynamics of the Gulf of Guinea areas bordering Ghana. The model results were compared with measure data. The model results shown similar trend for most of the quantity validated (i.e. water circulation pattern; water levels; and profile for temperature and salinity). Maximum deviations of 14.00% and 20.00% for water level and water velocity respectively make them appropriate to be used for hydrodynamic forecasting in the coastal waters of Ghana.

The text edit has been completed, the paper is ready for the template. Duplicate the template file by using the Save As command, and use the naming convention prescribed by your journal for the name of your paper. In this newly created file, highlight all of the contents and import your prepared text file. You are now ready to style your paper.

The authors of the paper gratefully acknowledge the contribution of MARETEC, of Instituto Superior Tecnico, Universidade de Lisboa, Portugal for providing the research team with the required measuring equipment and support. The authors also wish to thank Prof. Ramiro Neves and Dr. Ligia Pinto from the same institution.

The authors declare no conflicts of interest regarding the publication of this paper.

Uba, F., Essandoh, E.O., Nyantakyi, E.K. and Anumah, P. (2020) Hydrodynamic Model for Operational Forecasting in Coastal Waters of Ghana. Open Journal of Modelling and Simulation, 8, 48-59. https://doi.org/10.4236/ojmsi.2020.82004