A survey using the electrical resistivity and test pit methods is carried out in the Douala subbasin in order to investigate the sub-surface layering and evaluate the thickness of each layer through the variation of the resistivity with the depth. Applying the Schumberger array, a total of 48 vertical electrical sounding is conducted along 8 profiles directed W-E and using traditional method, 25 test pits with depth between 15 and 19 meters are done. Qualitative interpretation is carried out through the resistivity and test pit maps. The test pit maps reveal the presence of the following layers: The first layer of thickness varies from 0.2 to 2.0 m (average of 1.9 m). The thickness of the second layer varies between 0.2 and 4.0 m with an average of 2.8 m. The thickness of the third layer is between 2.0 and 4.0 m with an average of 3.8 m. That of the fourth and fifth layers ranges respectively between 4.0 and 6.0 m (average of 5.4 m) and between 2.0 and 6.0 m with an average of 4.4 m. Isoresistivities maps from resistivities computed for different AB /2 values (1.5; 4.4; 13.5; 19; 27.5; 58 m) are generally characterized by the apparent resistivities which define many domains in function of the AB/2 values. The superficial layer is characterized by electric resistivity values ranging from 720 to 2073 Ωm for AB/2 = 1.5 m and probably represents pedologic materials. For other values of AB/2, resistivity varies from 253 to 32,381 Ωm ( AB/2 = 4.4 m), 0.1 to 1838 ( AB/2 = 13.2 m) with an average of 341 Ωm, 0.3 to 1608 Ωm ( AB/2 = 19 m) with an average of 279 Ωm, 6.4 to 1427 Ωm ( AB/2 = 27.5 m) with an average 291 Ωm, and 19.7 to 3555 Ωm ( AB/2 = 58 m) with an average of 850 Ωm. These different values of resistivity represent the sedimentary materials.
The development of a country or even a region depends among other on the wealth of the subsoil. The superficial layer of soil being this interesting part of agriculture, the anchoring of the foundations of engineering structures, the search for certain natural resources, etc. sedimentary basins generally have a more or less complex structure and are sometimes interesting areas for the exploration of deposits. In the Douala sub-basin (south Cameroon, central Africa), several stratigraphic studies, such as Von Koenen (1896), Solger (1904), Riedel (1932), Dartevelle and Brebion (1956), Dartevelle et al. (1957), Reyre (1966), Dumort (1968), Njike Ngaha (1984), Regnoult (1986), Brownfield and Charpentier (2006), showed clay deposits in various formation, the main being the formations of Mundeck, Logbadjeck, Logbaba, N’kapa, Matanda, Souellaba and Wouri. However, the Miocene Matanda formation ( Dumort , 1968; Njike Ngaha, 1984) is one of the most interesting because of its stratigraphic position (it outcrops and is the least eroded, with a high potential for clays) and its geographical location (in and economic city with the advantage of considering and establishment in industries). In the Bomkoul area, Njopwouo and Wandji (1985) and Njopwouo and Kong (1986) used to show the kaolinitic nature of two identified raw clayey materials (gray and spotted clays) by chemical and mineralogical analysis. Elimbi & Njopwouo (2002) have also shown that these raw clays materials have a good characteristic for ceramics. Njopwouo (1984) studied the use of these clays materials in the polymerization of styrene and in the reinforcement of natural rubber. The preliminary work done by Mbog (2010) and Maliengoue (2010) describe the clay of materials of the Bomkoul area. In addition, the works of Ngon Ngon et al. (2012) have shown that these clays can be used in the manufacture of bricks, stoneware, and tiles, with the exception of pyrite rich clays. Ngon Ngon et al. (2012) reveal three facies: gray, dark gray and variegated. The studies conducted by Mbog (2016) focused exclusively on the differentiation of soil clay from the alteration of sedimentary materials, sedimentary clays which are deposits of fine sediments. However, despite these works, many shadows remain and require in-depth studies to better understand the mapping of the local and regional structures of this locality. In addition, the drillings and pits made so far have not completely crossed the different types of materials (sand, clay and sand stone) to ease the differentiation in depth their limits and geometry of deposit.
The study area is located in Douala sub-basin which forms the northern part of the Cameroon Douala/Kribi-Campo Basin. The Douala sub-basin is located in the gulf of guinea between the Cameroun Volcanic Line in the north and the Kribi-Campo sub basin in the south. It is characterized by low surface of erosion (<200 m of altitude) which corresponds to interfluves with convex and convexo-concave hillsides joined to narrow and either more or less deep valleys. The lithostratigraphy of Douala sub-basin is made of seven major Formations related to its geodynamic and sedimentary evolution (Nguene et al., 1992; Regnoult, 1986; SNH/UD, 2005). 1) The syn-rift period represented by the Mundeck Formation (Aptian Cenomanian) is discordant onto the Precambrian basement and consists of continental and fluvio-deltaic deposits, i.e. clays, coarse-grained sandstones, conglomerates. The post rift sequence includes; 2) The Logbadjeck Formation (Cenomanian-Campanian), discordant onto the Mundeck Formation and composed of microconglomarates, sand, sandstone, limestone, and clays; 3) The Logbaba Formation (Maastrichtian), maintly composed of sandstone, sand and fossiliferous clays; 4) The N’kapa Formation (Paleocene Eocene), rich in marl and clay with lenses of sand and fine to coarse-grained sandstone; 5) The Souellaba Formation (Oligocene) lying unconformably on N’kapa deposits and characterized by marl deposits with some interstratified lenses and sand channels; 6) The Matanda Formation (Miocene), dominated by deltaic facies interstratified with volcano-clasties layers; and 7) The Wouri Formation (Plio-Pleistocene) which consists of gravelly and sandy deposits with a clayey or kaolinic matrix.
The Bomkoul area is located to the north-east of Douala region between latitudes 48,050 to 48,060N and longitudes 98,470 to 98,490E in Douala sub-basin, within a humid equatorial climate. The geomorphology of the study area is a coastal plain with low altitudes (14 - 120 m). This area shows hills with flat and sharp summits ant is deeply incised by V and U shape valleys. According to the geological map of SNH/UD (2005) report, the relative age of the deposit of Bomkoul area is Upper Tertiary Miocene-Pliocene corresponding to the Matanda-Wouri-Formation (
In order to evaluate the geoelectrical resistivity map of the area, resistivity data was collected using an ABEM 1000 Terrameter following the step of 1 km. Fourty eight Vertical Electrical Sounding (VES) numbered S1… S48 and twenty four pits were done (
Electrical resistivity (ρ) is an intrinsic property of all earth materials. Resistivity is the reciprocal electrical current. In most porous rock systems, ionic conduction by interstitial fluids and surface conduction at the interface between solid rock matrices and electrolytic solution are responsible for a major part of the electric current flowing through a formation (Keleko et al., 2013). Electrical
resistivity techniques measure earth resistivity by passing and electrical current into the ground and measuring the resulting potentials created in the earth.
This method involves the supply of direct current into the grown through a pair of current electrodes and the measurement of the resulting potential through another pair of electrodes called potential electrodes. Since the current is known and the potential can be measured, an apparent resistivity can by calculated using Equation (1)
ρ a = π Δ V I [ ( A B 2 ) 2 − ( M N 2 ) 2 M N ] (1)
AB: current electrode spacing in meter, ΔV: potential differences in volts, MN: potential electrode spacing in meter, I: electric current in amperes, π = 22/7, ρa: apparent resistivity.
For Schlumberger soundings, the apparent resistivity values ρa were plotted against half current electrode spacing (AB/2) on a log-log graph and a smooth curve was drawn for each of the soundings. Then, the sounding curve was interpreted to determine the true resistivities and thicknesses of the sub surface layer. The resistivity of the subsurface material observed is a function of the magnitude of the current, the recorded potential difference and the geometry of the electrode array used. The depth of penetration is proportional to the Schlumberger array which uses closely spaced potential electrodes and widely spaced current electrodes. Generally, the depth of infiltration is small in this method and only shallow subsurface layers have been surveyed (Danielsen et al., 2007).
A total number of 48 Vertical Sounding (VES) with a step of 1 km were conducted along 8 geometric profiles. The soundings are arranged along profiles ranging approximately East-West. The Schlumberger electrode configuration having a maximum current electrode spread of 166 m was used. The current electrode spacin begins with a distance equal to 3 m and extends up to 166 m.
The results of the interpretation are represented in the form of the resistivity values that can be used for preparing the isoapparent electric resistivity map which reflect both lateral and vertical variations in resistivity. The surveys were carried out in dry weather facilitating implementation. The investigation area is very accessible in some places and the device could extend to over 240 m at the surface. Theoretical investigation depth is given by Barker’s formula (Banton & Bangoy, 1997) as Depth = 0.19 × L where L = AB for each measuring point.
The field results of the study are presented in qualitative interpretations of resistivity maps and quantitative interpretation of test pits. In the qualitative interpretations, the shape of the field curve is observed to get an idea qualitatively about the number of layers. The results involved isoapparent electrical resistivity maps and help to study the variation of the resistivity in the precise depth. This isoapparent resistivity maps reflect the lateral variation of apparent resistivity over a horizontal plane at a certain depth. In other words, theses maps indicate distribution of apparent resistivity in the area against distance of current electrodes. The maximum depth penetration of the AMNB method is 1/3 to 1/4 of the maximum distance of current electrodes (Frohlich et al., 1996). In the quantitative interpretation of test pits, layer thickness is obtained and helps to evaluate the real thickness of subsurface layer and better understand the variation of resistivity with the depth. Since current is conducted electrolytically by the groundwater (in the saturated layers) and surfacely by effective surface of minerals (in the dry layers) (Zohdy et al., 1974), the low values of apparent resistivity are attributed to the presence of saturated layer and the high values in other areas to the presence of unconsolidated and dry layer.
In the case study, isoapparent electric resistivity maps were constructed at AB = 3, 8.8, 26.4, 55 and 116 m. The choice of AB values depends on the variability of resistivity values observed on maps. These maps reflect the lateral variations of the electric resistivity at depths of 0.75 - 1, 2.20 - 2.93, 6.6 - 8.8, 13.75 - 18.33 and 29 - 38.67 m, respectively.
• Isoapparent electric resistivity map for AB = 3 m
Isoapparent resistivity map for AB = 3 m reflects the lateral variation over a horizontal plane at a depth of about 0.75 - 1.00 m. The resistivity varies enough increasing up to 2100 Ωm at some point and decreasing up to 700 Ωm at other points (
• Isoapparent electric resistivity map for AB = 8.8 m
Resistivities values are almost constant around 2000 Ωm in the entire area except the area around station S9 where it increases up to 36,000 Ωm indicating the intrusive formation.
• Isoapparent electric resistivity map for AB = 26.4 m
It shows a relatively constant resistivity (0 - 700 Ωm) but for some points there is a rise in resistivity. Exception is that of stations S1, S2 and S48 where the values of resistivities vary between 700 - 2300 Ωm. Low resistivity indicates the presence of aquifer in the area.
• Isoapparent electric resistivity map for AB = 38 m
It shows a relatively constant resistivity and varying between 0 and 700 Ωm except stations S1, S2 and S48 where the values of resistivities varies between 700 - 1800 Ωm. This slight variation of the resistivity compared to that obtained with AB = 26.4 m could reflect the fact that these two resistivity maps are obtained in the same layer. It indicates the presence of aquifer in the area at this depth except the zones of high resistivities.
• Isoapparent electric resistivity map for AB = 55 m
• Isoapparent electric resistivity map for AB = 116 m
The superposition of maps (
The description of different layers of clay materials (
In Bomkoul clays, silts, sandstones, nodular clays constitute the major field macroscopic differentiated facies. Several traditional sections located at different topographical levels have revealed thick clay deposits (19 m) and bottom to top, topped by materials representing almost the entire surface of the area with local
| Pits | Layer 1 | Layer 2 | Layer 3 | Layer 4 | Layer 5 | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| Thickness (m) | Feature | Thickness (m) | Feature | Thickness (m) | Feature | Thickness (m) | Feature | Thickness (m) | Feature | |
| P1 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P2 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P3 | 2m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P4 | 2 m | Sandy clay | 2 m | Sandstone | 4 m | Clay | 6 m | Compact clayrich | ||
| P5 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P6 | / | Sand | / | Sand | / | Sand | / | Sand | / | Sand |
| P7 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P8 | 2 m | Sand clay | 4 m | Sandstone | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Clay |
| P9 | 2 m | Sand clay | 4 m | Sandstone | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Clay |
| P10 | 2 m | sandy clay | 2 m | Sandstone | 4 m | Mottled sandy clay | 4 m | Clay | 6 m | Compact clayrich |
| P11 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P12 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P13 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P14 | 2 m | Sandy clay | 2 m | Sandstone | 4 m | Mottled sandy clay | 4 m | Clay | 6 m | Compact clayrich |
| P15 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P16 | 20 cm | Clay gravel and sandstone | 3 m | Clay nodular | 4 m | Mottled sandy clay | 4 m | Silts | 4 m | Sandy clay |
| P17 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P18 | 20 cm | Clay gravel and sandstone | 3 m | Clay nodular | 4 m | Silts | 6 m | Sandy clay | 4 m | Sandy clay |
| P19 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P20 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P21 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P22 | 2 m | Sandy clay | 2 m | Sandstone | 4 m | Mottled sandy clay | 4 m | Clay | 6 m | Compact clayrich |
| P23 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
| P24 | 2 m | Nodular clay | 20 cm | Clay | 2 m | Nodular clay | 4 m | Mottled clay | 2 m | Sandy clay |
| P25 | 2 m | Sand clay | 3 m | Clay nodular | 4 m | Mottled clay | 6 m | Sandy clay | 4 m | Sandy clay |
variations. Fifteen lithofacies were identified and described in
The representative facies are dominated by F1 to F9 and S1 to S3 codes.
A geological survey has permitted the description of the four clayey profiles in terms of texture, structure, distribution and color. The color was obtained using a Munsell soil color chart, and the terminology adopted for the description was that of Miall (1996) and Postma (1990). Fifteen lithofacies were observed in the four profiles and the traditional pits studied (
Twenty-five pits of 8 to 19 m were made to study the thickness of the different layers (
| Facies codes | Descriptions | Colour | Sedimentary Structures | Deposit process |
|---|---|---|---|---|
| S1 | Sandstone medium to coarse grained (30 cm to 2 m thickness) | 5Y7/6 to 5R4/6 to 8/1 | Horizontal laminations, oblique ferruginous massive red | Cohesive debris flows |
| S2 | Sandstone horizon end alternating way (1 - 4 m thickness) | 5R4/6 to 5R5/4 | ||
| S3 | Sandy clay (2 m) | 5R5/4 | ||
| F12 | Horizon clay gravel and sandstone fragments (10 - 20 cm) | 5R4/6 | ||
| F11 | Horizon nodular clay, gastropods and bivalves matrix (1 - 2 m) | 5R4/6 | ||
| F10 | Silts with gastropods and bivalves bioclast (4 m) | 5GY4/1 | Subhorizontal laminations, low bioturbations, 10 cm thickness | Suspension deposit in a quiet environment |
| F9 | Sand micaceous clay, gastropods and bivalves mould (2 m) | 5GY4/1 | Horizontal laminations | Planar bed flows with an average hydrodynamic |
| F8 | Horizon clay nodular (2 à 3 m) | 5R5/4 to 5Y7/6 | ||
| F7 | Micaceous sandy clay, coal, gastropods and bivalves mould (6 m) | 5G4/1 | Fine horizontal laminations, 60 cm thickness | Bed flows (Miall, 1978) with an average planar hydrodynamics |
| F6 | Micaceous sandy clay, coal, gastropods and bivalves mould (4 m) | 5G6/1 | Horizontal laminations | |
| F5 | Mottled clay, plates and some ferruginous nodules (2 - 4 m) | 5R4/6 | Fine horizontal laminations, 60 cm thickness | |
| F4 | Micaceous clay studded with gastropods moulds and crusts of hardened plates (6 m) | 5G6/1 | Horizontal laminations | Suspension deposits in a weak current environment with low energy (Miall, 1996) |
| F3 | Mottled sandy clay, presence of some ferruginous nodules (4 m) | 5Y6/1 | Horizontal laminations | Deposit in the lower part of a lower flow regime (Miall, 1978; Postma, 1990) |
| F2 | Clay with gastropods moulds and coal (4 m) | 5G6/19+ | Horizontal laminations | |
| F1 | Compact clay rich muscovite and charcoal fragments (4 - 6 m) | 5G6/1 | Horizontal laminations |
and 2.3 m with small thicknesses (domain 1) observed in the center and east of the study area and varying between 0.2 and 1 m. The second layer presents rather small thicknesses (domain 1) varying between 0.2 and 1.8 m in the North sector while the large thicknesses (domain 2) varying between 1.8 and 4.2 m are observed in the rest of the study area with a peak around wells 8 and 9. Layers 3 and 4 have almost the same characteristics: low thickness in the north and center of the map varies between 2 and 3 m for layer 3 and between 3.8 and 4.8 m for layer 4. The fifth layer has small thicknesses (domain 1) varying between 2 and 3.8 m in the northern sector while the values (domain 2) varying between 3.8 and 6.2 m are observed in the rest of the study area with peaks around pits 8, 9, 10 and 14.
Manual pits generally show layer thicknesses ranging from 2 m to 6 m. The perspective representation (
The analysis of resistivity maps coupled with the maps of the thicknesses of the soil layers obtained from the manual pits reveals information both on the variation of the resistivity of the soil and the geometric characteristics of the different layers of ground situated between the surface and the depth of investigation.
The first layer with an average thickness of 2 m is located between the surface and approximately 3 m of depth. The exceptions are P16 and P18 stations where there is a thickness of 20 cm. It has an apparent resistivity varying between 700 and 2000 Ωm at the first meter, reaching 35,000 Ωm between 1 to 2 m deep around the S9 station and a constant value of the resistivity (2000 Ωm) at this depth for the remaining of stations. The variation of the resistivity at the first layer reflects the distribution of materials of different natures such as sandy clays and lateritic clays while those observed in other zones reflect only the presence of sandy clays.
The high resistivity values (≥2000 Ωm) observed in this layer also reflect the presence of sandy clays, but with a cohesive debris flows like deposit process and a sedimentary structure which are oblique ferruginous massive red and horizontal/subhorizontal laminations. Those of weak resistivities have planar bed flows with an average hydrodynamic deposit process. Planar bed flows with an average hydrodynamic deposit process had been identified by (Mbog et al., 2016). The almost constant resistivity values between 1 and 2 m depth could indicate homogenization of the soil with an intrusion around the S9 station (
The thickness of the second layer varies between 2 and 4 m with the exception of the pit P24 where is 20 cm. It has high resistivity values up to 2300 Ωm in the northeast and southwest (
The third layer, situated between 5 and 9 m deep, has an average thickness of 4 m and the same variation of resistivity as the previous layer, but a slight difference is observed on the nature of the soil which varies between clay, nodular clay, mottled clay and silts.
The fourth layer is located between 9 and 15 m deep, with a thickness of 4 to 6 m and the same variation of the resistivity (
The fifth studied layer located between 15 and 19 m depth shows a significant variation of the resistivity (
Odunaike et al. (2013) reveal five to six layers of soil in Ogun State of Nigeria with the topsoil comprising of clay/sand/laterite, sandy clay, clayey sand, coarse sand and hardpan. The information reveals that the subsurface is representing a typical example of a sedimentary environment and the aquifer units.
According to the results of the resistivity survey carried out in the Douala sub-basin, it is now possible to reach the following conditions: the study of isoapparent resistivity maps showed that, generally the direction of resistivity contours is NW-SE approximately. This direction corresponds to the direction of a temporary river that drains the area. The distribution of apparent resistivity values in all maps showed that high values are due to coarse grain talus. Also, in other parts of the study area, the apparent resistivity values are the lowest due to fine grain alluvium. Five major layers were identified from the results of qualitative and quantitative interpretation of test pits in the sub basin. The first layer is interpreted as near-surface layer and has highly variable resistivity ranging from about 700 to 35,000 Ωm and a thickness of less than 1.9 m. The second one is dry alluvium layer with a resistivity of less than 2300 Ωm and a thickness of more than 4 m. The third layer that corresponds to the saturated layer has a resistivity and thickness of less than 2300 Ωm and 4 m respectively. The fourth layer is interpreted as a coarse-grained sand. This layer in most parts of the area is probably slate with a resistivity varying from 700 Ωm to 1600 Ωm and the thickness of 2 m.
The authors are grateful to all participants of the field works specially Pr Jacques ETAME.
The authors declare no conflicts of interest regarding the publication of this paper.
Mbog, M. B., Kenfack, J. V., Ngon, G. F. N., Tassongwa, B., Bayiga, E. C., & Etame, J. (2019). Morpho-Structural Mapping Constraints from Geophysical and Test Pit Investigations: Case Study of the Bomkoul Locality in Douala Sedimentary Basin, Cameroon, Central Africa. Journal of Geoscience and Environment Protection, 7, 136-153. https://doi.org/10.4236/gep.2019.710011