Physico-Chemical Characteristics of Gushing Water Aquifers in the Coastal Sedimentary Basin of Benin (West Africa) ()
1. Introduction
Groundwater is a valuable asset for human consumption and other multiple uses; this is especially true when groundwater is thermal (IEA, 2011). Also, they are easily accessible when they are gushing (Adissin et al., 2019). Indeed, gushing water points are particularly appealing because households access them without spending energy in pumping. Moreover, the permanent flow of these waters promotes their multiple uses for the development of small-scale agriculture (market gardening, irrigation) and livestock (Ahossi, 2020). As thermal water, they also serve other purposes such as energy supply and sanitary treatment (Kapasa, 2014).
In Benin, in the coastal sedimentary basin, boreholes have captured gushing water aquifers, sometimes thermal water in some communities (Boukari & Alassane, 2007). Besides the aforementioned uses, thermal waters are used for palm oil processing activities. Among the multitude of gushing water boreholes, one can notice an unexplained distribution of gushing and thermal water points to this day.
This study aims to characterize gushing water aquifers with an emphasis on the origin of the distribution of thermal abnormalities.
2. Presentation of the Study Area
2.1. Localization
The territorial scope of the coastal sedimentary basin of Benin encompasses the localities situated in the latitudes 7˚36' and 7˚12' North and longitudes 1˚36' and 2˚48' East (Figure 1). The total area is about 11,200 km2 (Boukari, 2002). From administrative standpoint, it includes the entire Atlantic, Coastal, Mono, Ouémé, Plateau and part of the Zou and Couffo provinces.
Figure 1. Location of the study area in Dahomey Bay.
2.2. Climate
In Benin, the annual average temperatures range from 26.5 to 27.5 degrees Celsius depending on the region (Boukari, 2002). The relative humidity is very high everywhere during the rainy season (80% on average in July and September). This rate drops during the warmer season (mainly between January and February) to about 60% in the South (Le Barbé et al., 1993). The annual rainfall system to which the coastal sedimentary basin belongs is equatorial (or sub-equatorial) and consists of four seasons (a large warm season from November to April, a large rainy season from April to July, a small warm season from July to September, and a small rainy season from September to November).
2.3. Geomorphology
The coastal sedimentary basin is located in the southern part of Benin and consists of sediments supported by a crystallophyllian base (Boukari, 2002). The morphology of the coastal sedimentary basin is characterized by two sets of plateaux arranged on either side of an ENE-WSW-oriented median depression, the Lama depression (Figure 2). The main rivers (Mono, Couffo and Ouémé) flow southward and form respectively the valleys of Mono, Couffo and Ouémé (Le Barbé et al., 1993).
Figure 2. Geomorphological map of the coastal sedimentary basin of Benin (Slansky (1962)
2.4. Geology
The geology of the Coastal Sedimentary Basin of Benin is the same as the Bay of Dahomey (Figure 3), which extends from Ghana on the west, and to Nigeria on the east (Slansky, 1962). These deposits have a monoclinal structure characterized by differential subsidence increasing toward the SSE (Dray et al., 1989). Based on lithological and sedimentological landmarks that support successive variations in sea level, eight stratigraphic units have been enumerated (Figure 3) from the Turonian-Coniacien (Unit I) to the Quaternary (Unit VIII), (Slansky, 1962; Lang & Houessou, 1978; Zevounou, 1992; Oyede et al., 2006; Zevounou et al., 2012): Unit I (Turonian-Coniacien): quartz sands with pebbles, intercalations of lignitous and marno-limestone levels; Unit II (Maastrichtian, Lower and Middle Paleocene): sands, argillites, marls, limestone; Units III to IV (Upper Paleocene, Middle Eocene): Kaolinitic clay, argillites, sometimes marl; Units V to VII (Mio-Pliocene, “Continental Terminal”): quartz sands, gravel; Unit VIII (Quaternary): sands and argillites.
Figure 3. Geological map of the Benin County Sedimentary Basin (extracted from geological sheets at 1/200,000 Lokossa-Porto Novo and Abomey-Zangnanado (IRB, 1989)).
2.5. Hydrogeology
By their lithological nature (Figure 3), the geological formations of the coastal sedimentary basin of Benin are grouped into four hydrogeological reservoirs. From north to south, there are listed as follow:
- The Upper Cretaceous sand aquifer or Unit I (Figure 3). It can reach up to 150 meters of thickness but remains mostly less accessible to the south of the basin due to its rapid sinking (Boukari, 2005). While its nappe is free in the north of the depression of Lama, where it emerges, it is captivated, in the South, under the argillites and the marls of both Maastrichtian and Eocene, where it sinks rapidly.
- The aquifer of the Paleocene limestone or IIb unit is essentially captivated under the argillites and marls of the Upper Paleocene and the Lower and Middle Eocene (Figure 3). Its thickness varies from a few meters to a maximum of thirty meters (Boukari, 2005). It is sometimes outcropping or sub-outcropping inside the Lama depression, then rapidly sinks southward and gets captivated under the clays in the valleys of the main BSC streams (Boukari & Alassane, 2007).
- The sand aquifer of the Upper Miocene-Pliocene or Units VI and VII (Figure 3). Its maximum thickness is generally around 120 meters and its nappe is usually free. Moreover, it is superficially subdivided into three hydrogeological units morphologically represented by the southern plateaus of the basin (Boukari, 2005).
- The Quaternary aquifers or Unit VIII are represented by the aquifers of the coastal sands, on one hand, and those of the alluvial sands on the other (Figure 3). Their maximum thicknesses are around 20 meters (Boukari & Alassane, 2007). In the alluvial valleys of the Mono River and especially the complex Oueme-Lac Nokoué, aquifers are often loaded under more or less superficial argillites with frequent phenomena of gushing artesian.
3. Material and Methods
The methodological structuring of this study is described in three phases: the collection of existing data; on-site work and processing and analysis of the data obtained.
3.1. Material
The collected data relate to the lithology of boreholes that capture the gushing water aquifers of the Upper Cretaceous and Paleocene, especially those with thermal water. The topographic sheets, the geological sheets (Abomey-Zangnanado, Lokossa-Porto Novo) made at a scale of 1/200,000 were also collated. The temperature of the borehole waters that capture the Cretaceous, Paleocene and Miocene-Pliocene aquifers has been collated. These data are obtained from the structures responsible for water management (General Directorate of Water, Department of Water) and the Benin Office of Geological and Mineral Research (OBRGM).
The on-site work is performed with the aim of mapping the geothermal abnormal points in the coastal sedimentary basin; therefore, a GPS has been used. The measurement of the physico-chemical parameters and the water sampling of the thermal boreholes allowed the geochemistry of the thermal anomaly zones to be highlighted. At each thermal water point encountered, the geographical coordinates are recorded along with pH and temperature values, conductivity and water collection.
3.2. Data Analysis and Processing
In order to better describe the geometry and structure of the gushing aquifers in the study area, two hydrogeological sections were performed along the north-south oriented lines AA' and BB'. Chemical analyses were performed on eight samples and were focused on some major water elements (
, Ca2+, Mg2+,
).
Temperature and depth data collected from the boreholes are entered into Surfer and QGIS to obtain maps of iso-depht and isotherms after validation with the coefficient of some statistical parameters such as the NASH and RSR coefficient of the Cretaceous, Paleocene and Neogene aquifers. Table 1 shows the characteristics of the fitted variograms.
Table 1. Characteristics of variograms.
Figure 4 shows respectively the variograms of customized models of the depth and temperature of the Neogene (aa'), Paleocene (bb') and Cretaceous (cc') of the coastal sedimentary basin. One can notice upon observation of these figures (Figure 4) that the model passes through the majority of points. The interpolation method used is Kriging and the power model was chosen for the kriging. This model adjusted for the raw variogram allowed us to obtain depth and temperature values in areas where there are no measurements (Lawin, 2007). The expression of the variogram of the power model is:
(1)
with
; C0 petite; h the distance
After the visual assessment of the overlay, the criteria of Nash (2) and RSR (3) allowed us to control the adjustment made so that the theoretical variogram better approaches the interesting part of the raw variogram (Figure 4). These criteria are defined as follows:
➢ Nash (coefficient of Nash, NSE)
(2)
The adjustment is optimal when the Nash coefficient efficiency is maximum (close to 1).
➢ RSR (standard observations deviation ratio)
(3)
The model is optimal when RSR is less or equal to 0.70 (Moriasi et al., 2007). This method allows us to identify the best set of parameters (C0, C and a) that characterize the modeled variogram (Lawin, 2007).
4. Results
4.1. Distribution of Captured Aquifers and Depth Variation
4.1.1. Analysis of AA' and BB' Lithostratigraphic Sections
The N-S lithostratigraphic correlations of section AA' in the Ouéme river watershed reveal six facies (Figure 5). From top to bottom, these are clay or laterite, gravel, clay sand, limestone and sand. The summit formations are represented by clays or laterites in the works. In the northern part of the section, the base is represented by the Gneiss and is located at a depth of about 100 meters. The sandy layer is located at a depth of about 90 meters and gradually sinks toward the southern part of the section. The thermal water points captured the sands at an average depth of 160 meters to about 200 meters. The temperature of these waters is between 39 and 50 degrees Celsius (Sagon and Gohissanou, respectively). In the southern part of the section, at the localities of Ahitonou, the formations captured are limestones. The Ahitonou borehole captures these limestones at a depth of about 60 meters with a temperature of 41˚C while at a few hundred meters away (Atchabita), the same limestones are captured at 200 meters with a temperature of 41˚C. The two boreholes are separated by the Zounvo river, which is a tributary stream of Ouémé river. This offset observed in the formations of these boreholes and the presence of the stream suggests that the area is crossed by faults. In the Southern part of the section, the limestones are located at a depth of more than 400 meters. Also, the thickness of the clays becomes consistent in this part of the section. Indeed, in the localities of Hetin-Sota, the limestones sink abruptly to a depth of about 417 meters. The water temperature of this locality is 69˚C. This sudden variation in the depth of the limestone shows that the area is crossed by faults.
Figure 5. North-South hydrogeological cut following the AA' cut line.
The N-S lithostratigraphic correlations of section BB' of the Couffo basin reveal seven types of facies in the studied area (Figure 6). They are from top to bottom clay, laterite or bar soil, gravel, clay sand, sand and gneiss. The summit formations are mostly represented by clays. There are a few levels of laterite sometimes. In the Northern part of the section, the localities of Affomai and Banigbe the clays have an average thickness of about 20 meters. Conversely, in the southern part of the section, from Adjame to the border of Possotomé, their thickness goes up to 140 meters with some layers of calcareous and marl levels. At the North of the section, the base is represented by the gneiss and located at about 70 meters depth. At the South of the section, in the localities of Agbodji, the sands subside to a depth of 300 meters while at five kilometers from these communities to the north of the section, the sands are located at an average depth of about 60 meters. The presence of Couffo and Ahémé lake streams in the proximity of the Possotomé and Eden thermal water boreholes (55˚C and 50˚C respectively) and the lags observed along the section suggest the presence of a fault system that was used to set them up.
Figure 6. North-South hydrogeological cut following BB' cut line.
4.1.2. Depth of Aquifers Captured in the BSC
Figure 7 shows the depth of the major aquifers in the coastal sedimentary basin. It is noticeable that the depth of aquifers increases from north to south. Indeed, the depth of the Cretaceous aquifer, which is about 40 meters in the northern zone where it outcrops, increases southwards to 280 meters under the outcrops of Paleogene marl and limestone of Lama Depression.
Figure 7. Aquifer depth map of the coastal sedimentary basin.
The Paleogene aquifer is 80 meters deep on average in the Lama depression, where it rises to 380 meters under the clays and marls of the valleys of the main BSC Rivers. The Neogene aquifer is captured on the southern plateaus of the BSC. The depth of this aquifer is 100 meters at the north and decreases to about 45 meters at the south of the plateaus.
4.2. Thermal Variation of Aquifer Water in the BSC
Figure 8 shows isotherm of water in the coastal sedimentary basin aquifers. The water temperatures of the major aquifers are noticed to be around 28 degrees Celsius to 34 degrees Celsius. Indeed, the temperature of the Cretaceous aquifer is about 29 degrees Celsius in the zone where it rises and 30 - 34 degrees Celsius under the Paleogene clays and marls. The temperature of the Paleogene aquifer varies between 27 and 30 degrees Celsius. In the Neogene aquifer, trapped under the southern plateaus of the BSC, the water temperature is about 28 to 30 degrees Celsius.
Although the isotherms of the main aquifers are between 28 degrees Celsius and 34 degrees Celsius, some points in these aquifers have temperatures ranging from 38 to 69 degrees Celsius. These thermal anomalies are found in the valleys of the main rivers (Figure 8).
Figure 8. Aquifer isotherm map of the coastal sedimentary basin.
4.3. Distribution of pH, Conductivity and Chemical Elements of BSC Thermal Waters
4.3.1. pH Variation of Thermal Waters
Figure 9 illustrates the pH distribution of thermal waters in the Coastal Sedimentary of Benin. The thermal waters are found to be acidic in the north (pH between 4.8 and 5.5) and neutral to basic in the south (pH between 7.3 and 8.5).
Figure 9. Distribution of thermal water boreholes and their pH in the BSC.
4.3.2. Composition in Bicarbonate, Calcium, Magnesium, Nitrate and Ammonium Ions in the Thermal Waters of the BSC
The table below (Table 2) presents the physical and chemical characteristics of thermal waters of the coastal sedimentary basin of Benin. It is noted that bicarbonate ions have a high concentration in waters collected on the southern part of the BSC (48.2 mg/L at 85.84 mg/L). Bicarbonate ion concentration is high from 48 mg/L to 85 mg/L in the south and low from 8 mg/L to 39 mg/L in the North. Conversely, in the north of the BSC, calcium and magnesium have a high concentration in the collected thermal waters. The concentration in ion Ca2+ range from 6 mg/L to 13 mg/L with exceptionally high values at Gohissanou (47.17 mg/L), Kpokissa (57 mg/L) and Hetin-Sota (64 mg/L). On the other hand, the composition of Mg2+ changes from 5 to 22 mg/L with an exceptionally high value in Hètin-Sota (78.86 mg/L). The conductivity values of the captured thermal waters are raised from 940 μs/cm to 1850 μs/cm in the south and low from 154 μ/cm to 170 μs/cm to the north.
Table 2. The Physico-chemical elements of thermal waters.
5. Discussion
From a geothermal and hydrogeological standpoint, the coastal sedimentary basin of Benin is characterized by thermal anomalies and artesian phenomena. Previous studies (Boukari & Alassane, 2007) have shown that aquifers that are the object of artesian phenomena are those of Cretaceous and Paleocene. The mapping of thermal anomaly zones shows that thermal waters are particularly observed in certain points in the valleys of the main rivers on either side of the Lama depression. The origin of thermal variations seems more in relationship with tectonics.
According to Taillefer (2017), the process of groundwater warming occurs based on three origins. We can distinguish the magmatic hydrothermal systems (Faulds et al., 2010), derived from the heat released by magmatic bodies; hydrothermal systems of geothermal gradient type (Taillefer, 2017), (deep aquifers from 2000 meters to 3000 meters) and hydrothermal systems resulting from the rise of heat from the lithosphere through the faults (rise of the isotherms) (Trabelsi et al. (2007)
The cross-sections made following the N-S directions reveal offsets in the geological layers which result in the thermal variation of the waters collected at Ahitonou (66 meters deep with a temperature of 41˚C) and Atchabita (deep 200 meters with a temperature of 27˚C). The continuous deformation zones are areas of weakness that serve as a channel for heat from the lithosphere and its distribution in nearby aquifers. Thus, the boreholes positioned in the areas of such phenomena will have a higher temperature compared to those in the vicinity. The random distribution of thermal water boreholes in the valleys of the main BSC streams is believed to be linked to tectonic events.
Since the chemical content of the waters is the result of the rocks crossed (Vernoux et al., 2014), the acidic behavior of the waters collected at the north of the BSC show that the northern part is under the influence of the basement. The aquifer captured to the north is from the sands of the Cretaceous which is made of major discordance at the base. On the contrary, in the South, the thermal waters have a neutral to basic pH. The high concentration of bicarbonate ion in these waters shows that they originate in carbonate reservoirs. This higher mineralization is mainly related to a solution of Sulfate chlorides originating from attapulgite-type clay levels and sublittoral deposits containing carbonates, phosphates and sulphates (Dray et al., 1989).
Tectonic manifestations are believed to be the cause of the warming of the gushing aquifers waters.
6. Conclusion
In the coastal sedimentary basin in southern Benin, the water of gushing aquifers captured in the valleys of major rivers has high temperatures, especially in some localities. These aquifers are made up of sands with an acid pH influenced by the crystalline base underlying the north, and limestone, then marl with a neutral or basic pH to the south. Clearly, the artesian phenomena of the aquifers of the valleys of the major streams are closely related to the topography which is a factor controlling their piezometry. The thermal water boreholes irregularly distributed among all the water boreholes gushing from the river valleys are said to be related to fault crossings. However, this hypothesis remains to be elucidated by our forthcoming linear studies, sedimentological and geochemical.