The Geology of Senegal shows a sedimentary basin covering more than 90% of surface of the Country. The formations of the sedimentary basin have ages ranging from Maastrichtian to Early Miocene upon which, we have recent Quaternary sedimentary formations [1]. Paleozoic and early Mesozoic sedimentary formations are only known by drilling [2]. These sedimentary formations lay down on a Cambrian to Precambrian basement shallowing up from West to East. Therefore, the formation of the basement outcrops on the Eastern part of Senegal, mainly from Bakel, to the Guinean and Malian frontiers [3].

This repartition of the geological formations describes also a varying topography decreasing mainly from East to West. This topography allows to the Senegal River to flow in the same direction, from its source in Guinea to the North of Saint Louis, in the vicinity of the Atlantic Ocean, Western Senegal (Figure 1). The Senegal River has a course of 1790 km. It started in Guinea and delineates the frontiers between Senegal and Mali on its Eastern section, and between Senegal and Mauritania in its Northern part [4].

Due to the international course of Senegal River (Figure 1), the crossed States (Guinea, Senegal, Mali and Mauritania) set up an intergovernmental organization to manage the River. The organization is responsible of managing the watershed of the Senegal River. This includes several damns built on the River for

hydroelectricity production and for sharing water for drinking and agriculture between the member States [5]. Another responsibility of the organization is to enhance and maintain navigability on the Senegal River [6] [7].

To enhance navigability on the Senegal River and to enhance regional integration and cooperation, it has been planned since a while, to build a port at the locality of Bakel. Bakel is located in Senegal. Nevertheless, it forms a triple point between Mali, Senegal and Mauritania. Its geographical position makes of Bakel a good candidate for building an international port, benefiting to as many member States as possible. Therefore, investigating the lithology and the geometry of the basement will provide, additionally to the geological knowledge of the Panafrican basement, valuable information in the preliminary design phase of the Port.

It is well known and established that the basement is outcropping when we move towards Bakel in the East. Subsequently the basement will be shallowing up toward this direction. Several studies have been performed in the aim to map the depth of the basement [8] [9] [10] [11], nevertheless, the scale and the resolution of these studies were not high enough to provide sufficient information on the thickness of the sedimentary layer along the Senegal river course.

The aims of this study is to map at a high resolution, in the locality of Bakel, the geometric relationship between alluviums, under the Senegal river bed, and substratum, to assess the hardness of the rocks present underground using a geophysical investigation method, more specifically seismic refraction tomography technique which gives the underground repartition of seismic P-wave velocities.

The site of this study is located in Bakel, Eastern Senegal. Bakel is on the frontier between the Meso-cenozoic Senegalo-mauritanian sedimentary basin. It corresponds to the outcropping of the Proterozoic basement, locally represented by the Panafrican Mountain chain [3]. The Panafrican chain is also referred to as Mauritanides. The Mauritanides have been successively affected by the Panafrican orogeny (1000 to 595 Ma) and the Hercynian orogeny dated between 300 and 250 Ma [8] [12].

The Mauritanides is a thrust belt long of 2000 km surrounding the West African Craton [13]. It includes three main parts: the Adrar Soutouf Massif in Southern Morocco [14], the Central Mauritanides in Mauritania, starting from Bakel in Senegal to Southern Morocco referred to as Anti-Atlas [15], the Bassarides Belt in the Northern Guinea and Southern Senegal [10] and the Rokelides Belt from Guinea to Liberia [16] [17].

At Bakel, the Mauritanides is locally represented by formations referred to as the “Group of Bakel” belonging to the “Complex of Bakel”. The Complex of Bakel is a thrust zone formed of several groups of meta-volcanic and meta-sedimentary allochthones units. In fact, the Group of Bakel overlaps allochthone units in the East and is overlapped from the West [18].

The “Complex of Bakel” shows several hills around Bakel. The hills are formed of schists showing inclusions of sericite or hematite and reinforced with quarzitic bars.

The thrust belt formations described above present local plains occupied by Holocene alluvial deposits linked to the regime of the Senegal River and/or its tributaries [19]. More details about the structure of the Mauritanides thrust belt in the vicinity of Bakel and the studied area can be found on the geological map (Figure 2).

The geometric relationship of the rock types around the Senegal River bed in the site of study, the alluvium layer, the weathered bedrock and the depth of the fresh hard bedrock has been investigated using 4 seismic refraction profiles (Figure 3). In fact, seismic investigation techniques are widely used in soil mechanical properties studies, particularly in rippability or ease of excavation assessment of rocks which is a significant aspect of the preliminary work of any civil engineering project [20] [21] [22] [23]. Seismic Refraction is a cost effective powerful investigation tool for shallow survey, and is increasingly used for small scale mapping of near surface, particularly in site investigations for civil and geotechnical engineering [24]. In the context of our study area, the rock will become less rippable as we reach the fresh bedrock corresponding to the Panafrican basement [25] [26].

Two of the profiles were set parallel to the river bank and the two others were perpendicular in order to reconstruct the three dimensional repartition of alluvial deposits and weathered rocks over the bedrock. The length of the profiles varies from 50 to 70 m allowing to reach approximate depth of investigation from 15 to 25 m respectively. The length of the profiles was driven by the landcover constraints. In fact, the study site is bordered in the West by the City of Bakel, and in the East by the river. Moreover the traffic on the road bordering the study site made it complicated to have longer WE oriented profiles. This situation resulted in river parallel profiles longer than perpendicular ones limited by the extend and the human activities of Bakel City. Nevertheless, the obtained profiles length, which controls the depth of investigation was enough to achieve the targeted depth of investigation that was 15 to 20 meters.

For the different profiles, the receiver spacing varies between 2.5 and 5 m to achieve a good resolution.

We used a 24 channels Daqlink III seismograph with a set of 10 Hz geophones. The data acquisition parameters for gain, sampling rate and time window were set respectively to 1, 0.125 millisecond and 0.5 second. The seismic source, was an 80 kg mass dropped by an electric motor from 2 meters height on a plate fixed on the soil. For each profile, 11 shots were performed consisting of 4 offset shots, 2 end shots, 4 internal shots (between geophones 3 - 4, 6 - 7, 18 - 19, 21 - 22) and 1 center shot between geophones 12 and 13. For more accurate results, the coordinates for the sources and receivers were collected using a Leica differential GPS.

The first arrivals have been picked on each shot of the acquired seismic refraction data under Vibrascope which was also used as acquisition software (Figure 4). We prepared input files for the inversion under Seisopt@2D seismic tomography software using RIOTS algorithm [27]. The three input files for the inversion software consist of:

• The source file, containing the coordinates of all the shots for a given profile. For a given profile and for each shot the position (x), the altitude (z) and the number of receivers (N) are specified.

• The receiver file, containing for each shot, the position (x) and the altitude (z) of each of the 24 geophones.

• The picking file, containing for each shot, the time of first arrival for each receiver.

The inversion parameters and the quality control of the output model were tuned mainly by comparing observed and predicted hodochrons (Figure 5).

The inversion produced four seismic 2D profiles showing P-wave velocity variation with depth (Figure 8). Even if the 2D profiles show spatial change of the seismic P wave velocity with depth, they do not allow understanding easily their geometrical relationship and the tridimensional geometry of the bedrock and the overlaying alluvium. To integrate the 2D profiles into a tridimensional model, we computed a 3D block model using the Rockware Rockworks® software.

The obtained 3D block model allows a better comprehension of the geometry of the existing geological formations. And to achieve a major goal of the study which consists of identifying the different types of lithology, we need to interpret the seismic P wave velocity to rock type. The interpretation of seismic velocity to rock type was done based on the outcropping rocks along the seismic profiles and confirmed using Caterpillar rippability chart [28] presented below (Figure 6).

The main steps of our approach are summarized in the flow diagram below (Figure 7).

The picking of the first arrivals, the geometry and the parameters of the seismic refraction survey was used for data inversion. Inversion allowed to reconstruct the seismic velocity model from the observed data [29]. The inversion of the seismic refraction data resulted in 4 seismic profiles showing the repartition of P-wave velocity with depth (Figure 8).

The obtained results show that the measured seismic P-wave velocities vary from 290 to 5200 m/s. All the four profiles show increasing velocity with depth. The comparison of seismic velocity repartition with outcropping geological formations allows to deduce that alluvium, in the upper layer, shows lower velocities and the fresh and hard rock shows higher velocities. The transition between the velocities of alluvium and hard rock is formed by the weathered rock.

To determine the lithology from the obtained seismic tomography profiles, we used information from the outcropping rocks in the study area to convert the P-wave velocity to rock type with the following approach: we consider the three main existing facies: the alluvium, the fresh quartzite and the mixture. The mixture is formed of alluvium containing big blocks or weathered quartzite (Figure 9).

The location A shows a massif of weathered rock passing into its center to hard rock and corresponding to the Southern end (left part) of profile 4. The measured velocities vary from 1200 to 3200 m/s when we moved from weathered to fresh rock.

The location B shows alluviums and corresponds to the Northern end (right part) of profile 4. The measured velocities vary from 290 to 800 m/s.

We can deduct from these two observations that the transition between these two intervals, is geologically formed by weathered rocks with velocities varying from 800 to 1200 m/s. The correlation between seismic P-wave velocities and lithology is summarized in Table 1.

The repartition of the P-wave velocity based on outcropping rocks has been confirmed based on the Caterpillar D9R rippability chart (Figure 6) that slightly consider the same intervals for alluviums, weathered and fresh bedrock based on their rippability.

The obtained seismic velocity data have been compiled to a database containing the UTM WGS84 x and y coordinates, the z level and the velocity. The obtained database has been used to build a block model of the seismic velocity in

the investigated area. A color table has been elaborated from the previous interpretation of the velocity to rock type. The results of the tridimensional modeling of the data are presented below.

The tridimensional model of the investigated area (Figure 10) shows that the small gulf observed in the Southern part of the investigated area is formed of hard fresh rock passing laterally to weathered rock (Figure 10(a) and Figure 10(c)). The gulf is followed by a bay covered with alluvium (Figure 10(b) and Figure 10(c)). The thickness of the alluvium deposits increases toward the North corresponding to the flowing direction of the Senegal river (Figure 10(a) and Figure 10(d)).

Partly weathered bed rock is obviously present in the transition between the alluviums and the fresh bedrock. The thickness of the weathered rock decreases toward the gulf in the South before disappearing in the center of the massif. Toward the West, the thickness of weathered rock increases when we move far from the river banks.

The fresh bedrock shows an irregular surface, deepening toward the East and the Nord. We can also notice that the maximum depth corresponding to maximum alluvium thickness is obtained in the bay. The deepest point for the bedrock, located in the center of the bay, is at 0 m above sea level, while the highest point is at 18 m and is observed on the hills boarding the river. The reference topographic level corresponding to the outcrops of fresh bedrock observed in the massif of the gulf is 14 meters.

The tridimensional model of the investigated area shows that course of the Senegal River is linked to the lithology. In fact, the sharp bending of the river is due to the massif of fresh and partially weathered rock forming a gulf that seems to block the sediments upstream. Right after the gulf, alluviums reappear and let to the river the possibility to shape a large bay downstream. The lithology is in

its turn linked to the morphology as, the thickness of the alluvium is very important in the bay and less important in the gulf which corresponds respectively to low and high topography. In fact the sedimentation and erosion rates are controlled by natural and/or anthropic factors [30] [31]. Both erosion and deposition are continuous processes in the river in an attempt to reach a new equilibrium in channel geometry [32]. In the case of the study area, the massif of rock, responsible of the sharp bend of the river dissipates considerably the energy upstream the bay. The barrier formed by the massif seems to be responsible of the decrease of energy that allows more important sedimentation in the bay, while the sharp bending of the river course and the decrease of the depth, induce breaking waves and subsequently important water agitation which cause erosion in the cap.

Moreover, he obtained model allowed estimating the volume of alluvium that could be removed by dredging, in the investigated area which is around 13,000 m³, while the volume of weathered rock is 19,000 m³. The total surface occupied by the model is approximately 7500 m².

The seismic investigation of the Senegal River banks in the City Bakel gives valuable information on the geometric relationship between the Proterozoic Panafrican basement formed of hard rocks and the overlaying alluvial deposits. Also, to enhance the navigability and insure draught requirements for future boats and vessels navigating in the port, it is necessary to dredge the riverbed. The dredging can involve alluvium, very weathered bed rock that the obtained model can allow to assess and estimate.

The seismic velocity model allows mapping the spatial repartition of the rock types that seem to have a huge influence on the course of the river, the erosion and deposition of the sediments and therefore, their geometry. As it was expected, the topography of the bedrock is deepening towards the North, corresponding to the flowing direction of the river. Finally, the obtained tridimensional seismic P-wave velocity model allows to know that the dredging will be easier toward the North.

The sharp bending and the rapid change in lithology lead to suspect the presence of an EW fault or discontinuity between the gulf and the bay. Further investigations are needed to confirm its presence, its relationship to the geological framework, and eventually its regime, as that can interact with the projected infrastructures. Also, decreasing the spacing between the profiles can increase the resolution and give more accurate repartition of seismic velocities. It is worth to notice that all the profiles were localized onshore; acquiring offshore data can widen the investigated area and give valuable information to Bakel port construction planners. Finally, performing some drills will allow to tune the interpreted model and be more accurate on the lithology corresponding to each measured seismic P-wave velocity.

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

Ndiaye, M., Sall, O.A., Thiam, A., Sarr, D., Badji, M. and Ndoye, I. (2020) Investigating the Depth and the Geometry of the Quarzitic Panafrican Basement Using Near-Surface 3D Seismic Refraction Tomography: Case Study of the Locality of Bakel (Senegal). International Journal of Geosciences, 11, 345-359. https://doi.org/10.4236/ijg.2020.116018