Hydrodynamic Model for Operational Forecasting in Coastal Waters of Ghana ()
1. Introduction
The provision of economic resources to countries at the shore of Gulf of Guinea makes it very important to countries in West Africa [1]. Taking it as a system it will be important to know the dynamics (sea surface temperature, changes in density, current pattern, stratification of the sea, etc.) that occur in it to help neighouring countries predict its state at all times [2] [3] [4]. The use of numerical tools has made this possible. Numerical tools have the ability to help stakeholders link cause to effect, estimate and predict changes; and also as a tool to manage environmental phenomena with the best approach i.e. phytoplankton growth, water quality disturbance, etc. [5] [6] [7]. Oceanic circulation is caused by two processes namely, wind-driven or density-driven circulation and/or both and this is the key cause of bio-physical processes taking place in the Gulf of Guinea. Oceanic circulation influences all the activities in the oceanic ecosystem, including sediment transport, turbidity pattern, particle dispersal, etc. One other notable aspect of the environment ecosystem is the distribution, health and economic impact of bio-chemical particles (e.g. pollutant) in the marine world [5] [8] [9]. Circulation current patterns are of considerable interest to the Gulf of Guinea indigenous community who depend on the natural resources for their livelihood, as well as to industries that operate the region such as fisheries and shipping [1].
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 [1] [10].
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.
2. Methodology
In this research, Modelo Hidrodinâmico, MOHID model, a 3-Dimensional water modeling system was used [2]. It is designed in such a way that the domain is divided into various medium; land, air and water. The architecture is such that communication is done between media by interface modules (i.e. InterfaceSedimentWater Module and InterfaceWaterAir Module). These interfacial modules handle the fluxes between two media, either air and water or water and sediment.
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 (2nd and 3rd), 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.
2.1. Implementation
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;
(1)
where
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:
(2)
(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;
(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 xi can be divided into three contributions;
(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 [7] by considering salinity and temperature;
(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 [1] [8] [11] [12] to algebraically represent and solve the equations. The algebraic representation of the equation is applied to a cell control volume.
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 [3] manner. Arakawa C arrangement allows velocities to be located at the center of the cell under consideration; horizontal velocities are located in the center of the west (u-velocity) and south (v-velocities) faces, while elevation, turbulent magnitudes and tracers are placed in the center. Also, a staggering in the vertical direction is used, with vertical velocity w, tracers and turbulent magnitudes vertically placed in the top and bottom faces and horizontal velocities and elevation in the center of the element (in vertical). The temporal discretization described in [13] using semi-implicit algorithm: Alternate Direction Implicit was used. In all simulations, the 6 equations algorithm discretization schemes by [14] are used.
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 [11].
2.2. Boundary Conditions
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:
(7)
By means of the wind surface stress, τw, the Diffusive flux of momentum is explicitly imposed:
(8)
Wind stress is estimated from a quadratic friction law:
(9)
where drag coefficient, CD, 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:
(10)
The drag coefficient, CD, is estimated from the equation:
(11)
where κ and
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 [15]. Also, at the bottom, a zero flux is assumed for both salinity and temperature.
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;
(12)
(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 [4] [7] [16].
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 [12].
3. Results and Discussion
Figure 1 shows the levels considered during the hydrodynamic implementation to reduce computational time and memory allocation. The last level containing the 3rd domain is the smallest with higher resolution compared to the first and second levels containing the 1st and 2nd 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. Figure 2 shows the locations of the tidal gauges for the four places.
Figure 3 shows graphs of water level obtained from the field and also from the hydrodynamic model and Figure 4, Figure 5 show temperature and salinity profiles of the same hydrodynamic model.
The hydrodynamic pattern for the region is also shown in Figures 6(a)-(d) for particular time during the simulation period.
Figure 1. Resolutions of the various domain LEVELS.
Figure 2. Tide Gauges location for validation.
(a)(b)
Figure 3. Validation of Hydrodynamic Model ((a) Takoradi Tide Gauge Series, and (b) Tema Tide Gauge Series).
(a)(b)
Figure 4. Salinity profiles for Hydrodynamic Model.
(a)(b)
Figure 5. Temperature profiles for Hydrodynamic Model.
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 (Figure 3) at an accuracy margin of 14.00%. The difference in tide measurements especially the simulated results deviating slightly from the measured data can be attributed to; the proximity of the installed gauges nearer to the shores. When tidal gauges are installed near the shore their readings are affected for example, by human activities around the tidal gauges; docking ships at the habour (which may agitate the sea water causing ripples) and obsolete measuring tools like dialed instrument can give some margin of errors. Also due to sedimentation from shoreline morphology or silt transportation some errors are introduced. Another contributing factor is the bathymetry resolution or poor bathymetry accuracy for lower areas. Hence, the differences in results.
Figure 6. Hydrodynamic pattern for sellected dates.
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.
4. Conclusions
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.
Acknowledgements
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.