Use of the Polynomial Separation and the Gravity Spectral Analysis to Estimate the Depth of the Northern Logone Birni Sedimentary Basin (CAMEROON) ()
1. Introduction
The Logone Birni sedimentary Basin is located in the far north region of Cameroon, south of Lake Chad (Figure 1(a)). It covers an area of approximately 27.000 km2 and it is part of the West and Central African Rift system (WCARS). The West and Central African Rift system is divided into two Cretaceous subsystems genetically linked, but physically separated [1] , named the West African Rift Subsystem (WARS) and the Central African Rift Subsystem (CARS). [2] indicated that these two rift subsystems, although genetically related, present some structural differences. To study the structure and understand the tectonic setting of the Logone Birni Basin, geophysical analysis and interpretation based on seismic reflection data coupled with gravity, magnetic and remote sensing data were used. The structure basement conditions of basin infill and distribution of buried volcanic bodies’ were elucidated. The results confirmed the basic knowledge of the ancient tectonic and geophysical models of WCARS [2] [3] .
Figure 1. (a) Simplified geological map of Cameroon and adjacent areas (modified from [2] ). The map shows distribution of different sedimentary basins: localization of the Logone Birni Basin in far north of Cameroon; (b) Simplified geological map of the study area (modified from [5] ).
Significant hydrocarbon deposits were discovered in the sedimentary basins which are part of the WARS and CARS (e.g., the Doba basin in Chad, Muglad and Melut basins in Sudan). Despite the scientific and economic importance of these zones, recent surveys date from the 1960s to the 1990s [1] [4] , and most of the data are owned by oil companies, and not always available to universities.
In this study, we calculated the average depths of the basin using spectral analysis of the residual anomaly maps obtained after filtering the Bouguer anomaly map computed by [4] [5] the northern Logone Birni Basin (Figure 1(a)).
2. Regional Geology
The study area (Figure 1(a), Figure 1(b)) is located in the transition zone between the Central and the West African Rift System (WCARS) [6] [7] . This area was affected by the Pan-African orogeny (750 - 550 million years) which generated the main part of lineaments and faults in the Palaeozoic (550 - 130 million years), the Cretaceous (130 - 75 million years), the Maastrichtian-Palaeogene (75 - 30 million years) and the Neogene-Recent (30 - 0 million years).
The region is under a tropical climate, which favours the growth of a shrubby savannah and sparse forest vegetation. The Lake Chad region has the lowest altitude ( 280 m ). The main rivers that drain the region are the Chari and the Logone; the Chari is tributary to the Lake Chad while the Logone is an effluent of the Chari. The geological history of the region is related to the formation of the Chad basin.
Throughout the Lower Cretaceous, the Chad basin was an extension zone between the western and central Africa plates. It was subject to various tectonic activities from the primary era to the Quaternary era marked by relief inversions [8] [9] . The present sedimentary cover which has an average depth of 600 m or more is the result of an interaction between two main factors: the Lake Chad transgression in the humid period and the lacustrine transgressions during the more arid periods marked by Aeolian erosion that strongly removes deposited materials [10] [11] [12] . The lacustrine transgressions can be delimitated over three main geological periods. This delimitation was made possible thanks to the discovery of sedimentary basins that correspond to palaeo-shorelines of the lake [13] [14] [15] [16] :
- The first transgression might have started at the end of the Tertiary, resulting to basin infill in the Middle Quaternary from material derived from partial ablation of ancient deposits and border massifs. The process was amplified by subsidence of the basin and, inversely, by massifs uplifting.
- The second transgression dates from about 21.000 years.
- The third transgression linked to the third humid period, dated from 12.000 years ago and marks the constitution of the first trace of the actual drainage pattern.
The northern Logone Birni Basin is a sedimentary plain covered by sandy clayey alluvial deposits of quaternary age. Its lowest altitude makes it vulnerable to regular flooding by the River Chari and the Lake Chad. The sedimentary cover date from the Tertiary to the Quaternary throughout the region. This cover is made by river alluvia, lacustrine, wind sediments, and shows three main series [10] [14] (Figure 1(b)):
- The Bodele series constituted by lightly differentiated sediments of the Upper Tertiary composed of sand and sandstone with some clayey intercalations. This series which date from the Pliocene is mainly fluvial.
- The Soulias series (Middle and upper Pleistocene) constituted by Aeolian sand and lacustrine limestone.
- The Labde series, which includes thin lacustrine deposits date from 2400 years to recent.
The sequence is underlained by granitic, gneissic and migmatitic basement rocks which appear at variable depths.
3. Methods
3.1. The Polynomial Separation Method
The polynomial separation method was used to produce the first, second and third degree residual maps. The algorithm by [17] [18] was used to adjust the polynomial surfaces to the Bouguer anomaly map [5] [19] [20] . This method is based on the analytical least square method and the polynomial decomposition series.
The least-square method was used to compute the mathematical surface which gave the best fits to the gravity field within specific limits [5] [19] [20] . This surface is considered to be the regional gravity anomaly. The residual was obtained by subtracting the regional field from the observed gravity field. In practice, the regional surface is considered as a two-dimensional polynomial. The order of this polynomial depends on the complexity of the geology in the study area. The first, second and third-order polynomial surfaces of the regional anomaly obtained in this work is not presented, but the corresponding residual anomaly is presented in Figure 3.
Mathematic Formulation of the Method
The Bouguer anomaly B(x, y) in the given point M(x, y) of the earth in cartesiane coordinates is governed by the relation:
(3.1)
is the sum of the residual anomaly
and the regional anomaly
.
The surface
which is adapted to the gravity field data
is given by the following relation (Radhakrishhna and Krishnamacharyulu, 1990):
(3.2)
where N is the order of the polynomial,
the number of
terms of the polynomial, and
the coefficients to be determined:
The first order Polynomial is:
(3.3)
The second order polynomial is:
(3.4)
The third order polynomial is:
(3.5)
We denote by
the difference between the homologous points of the experimental and analytical surfaces respectively and by N0 the number of stations Pi in which the Bouguer anomaly is known. The adjustment of the surfaces which consists in making the quadratic deviation minimal is expressed by:
then
with
(3.6)
and
We then obtain a system of (M) equations with (M) unknowns. The unknowns are the coefficients
of the polynomial
of order N. Once the coefficients are determined we determine the analytic regional anomaly
and the residual by:
(3.7)
The polynomial method is particularly used when the amplitude of the residual anomalies is negligible compared to the regional one. Apart from the polynomial method there are other methods such as the upward continuation method.
3.2. The Spectral Analysis Method
The spectral analysis method was used to investigate the wave numbers of gravity residual anomalies and to estimate the depths of the bottom and top of the residual anomalies source bodies.
Many scientists used the calculation of the power spectrum from the Fourier coefficients to obtain the average depth of the disturbing surface or equivalently the average depth to the top of the disturbing body [21] [22] [23] . In order to estimate the depths of the bodies responsible for the observed gravity anomalies on the residual maps (Figure 3), it is necessary to define the power spectrum of a gravity anomaly in relation to the average depths of the disturbing interface. It is also important to note that the final equations are dependent on the definition of the wave numbers during Fourier transform. The fast Fourier transform (FFT) has been applied to an anomaly profile with n data.
Mathematic Formulation of the Method
The approach is based on the following procedure that allows the determination of different depths [5] [24] . For the observed residual anomalies:
(1)
In this function the wavenumber k is defined as
and
is the amplitude coefficients of the spectrum is defined as,
(2)
For an altitude of z = 0, the amplitude coefficients of the spectrum can be written as,
(3)
Then the amplitude coefficients of the spectrum can be rewritten in terms of (3) as,
(4)
The power spectrum
is defined as,
(5)
Taking the logarithm of both sides,
(6)
We can plot the wavenumber, k, against
to attain the average depth to the disturbing interface. The interpretation of this plot requires an estimate of the best line fit of the lowest wavenumbers where a change in gradient is observed. The average depth can be estimated from Equation (6) as,
(7)
where
is the average depth,
and
are derivative of P and k respectively. The observed residual anomalies depths using power spectral approach as shown in Figures 4-6.
4. Results and Discussion
4.1. Interpretation of Bouguer Anomaly Maps
A map of the Bouguer anomalies represented with dotted lines is shown in Figure 2(a). This map contains about 700 measurement points in a zone limited between latitudes 12˚11' and 13˚16' North, and longitudes 14˚16' and 15˚26' East and covers an area of approximately 16,200 km2. This map with a 5 km grid size was obtained by kriging, with a variogram model that has an angle of 0 degree and an anisotropy coefficient of [4] [17] . The Bouguer anomaly map of the Northern Logone-Birni area (Figure 2(b)) is deduced from the map of Figure 2(a)
Figure 2. (a) Map of grid gravity stations distribution of Bouguer anomalies from [5] modified; (b) Bouguer anomaly map of the region from [5] modified.
with data uniformly sampled. This Bouguer gravity map is shown in Figure 2(b) with contour interval of 5 mGal. The observed data were separated into regional and residual anomalies before interpretation.
4.2. Interpretation of Residual Anomaly Maps
Qualitative interpretation is based on the geophysical information that can be extracted from the residual anomaly maps and their relationship with the geology of the region. The main features of the residual anomaly map (Figure 3) are very similar to those of the Bouguer anomaly map (Figure 2(b)).
The Bouguer anomaly map (Figure 2(b)) represents the resulting effects of the near surface and deep geological units within the survey area. This Bouguer anomaly map is characterized by the relative strong negative anomalies around Goulfey (−49 mGal), Tom Merifine (−44 mGal) and Tourba (−44 mGal). Those anomalies correspond to the northern Logone Birni basin infilled with sediments ( [5] [10] [13] [15] ). This map is also characterized by relative anomalies highs at Hadjer El Hamis and Djermaya which are in agreement with the uplifting of the basement as a result of magmatic intrusion in the corner of the basin by denser rocks such as rhyolites [5] [12] [15] [24] [25] .
The residual anomaly map (Figure 3) obtained after polynomial separation portrays positive and negative values which may indicate an intrusion of the basement into sediments, and the northern Logone Birni sub-basin respectively. The residual anomalies found in the area are almost elliptical in shape and trend N-S, NW-SE, NE-SW; they are accompanied by a steep dipping gradient, which can be associated with structural geological boundaries.
The results of the structural trends of different anomalies are NW-SE and NE-SW. Tables 1-3, are in agreement with the results of [2] [3] and [7] which show that structural style in the Logone Birni Basin are mainly NE-SW. These styles correspond to tilted fault blocks associated with NNE-SSW senestral extension movements. While the NW-SE directions correspond to the directions WNW-ESE with normal faults and fault blocks associated with ENE-WSW dextral transtensional movements.
The structural directions observed here are similar to those of several major oil fields established in Sudan, for example, Muglad and Melut basins [25] [26] .
4.3. Interpretation of Spectral Analysis Profiles (Segments)
The spectral analysis curves from profiles P1, P2 and P3 of Figures 4-6 shows the variation of the logarithms of the power spectrum with respect to the wavenumbers, and the slopes (H) selected to estimate the depths. The shape of the power spectrum depends mainly on the values of the wavenumber along the profile. These values can be qualitatively estimated by the shape of the profile which depends on the geology of the area. The high wavenumbers (0.3; 0.6) are related to the shallow anomalies, and the low frequencies (0.0; 0.3) are related to the deepest anomalies ( [2] [4] ).
The profiles (P1, P2 and P3) have been selected along the large residual anomaly observed on the three residual maps, which correlate with the northern Logone Birni Basin (Figure 3). P1 is approximately orientated NW-SE and crosses the town of Goulfey. P2 and P3 are orientated NW-SE and cut-across the towns of Tom Merifine and Tourba.
Table 1 and Figure 4 show along profiles P1, the deepest anomalies with a
Figure 3. Residual maps of the study area showing gravity signature of the northern Logone Birni basin and profiles P1, P2 and P3: (a) The first order residual anomaly map of the region showing the profiles; (b) The second order residual anomaly map of the region showing the profiles; (c) The third order residual anomaly map of the region showing the profiles.
Figure 4. Power spectrum of profiles P1, P2 and P3 associated to the first degree trend gravity map.
Table 1. Depths obtained by spectral analysis and characteristics of residual maps and profiles, for first degree polynomial fitting.
Table 2. Depths obtained by spectral analysis and characteristics of residual maps and profiles, for second degree polynomial fitting.
Table 3. Depths obtained by spectral analysis and characteristics of residual maps and profiles, for third degree polynomial fitting.
depth of about 4.55, 4.00 and 1.63 km , while the depth of the shallowest anomalies is at 0.81, 0.47, and 0.25 km .
Table 2 and Figure 5 show along profiles P2, the deepest anomalies with a depth of about 2.10, 2.09 and 1.67 km and the shallowest anomalies with a depth of 0.34, 0.34 and 0.24 km .
Table 3 and Figure 6 show along profiles P3, the deepest anomalies with a depth of about 3.26, 2.96 and 3.84 km and the shallowest anomalies with a depth of 0.73, 0.49 and 0.49 km .
The estimated depths are quite close to those reported by [9] and by [24] (0.325 - 0.600 km ). Similar estimated depth values were also obtained using a combination of geomorphological data, Lansat images and resistivity sounding and, by [1] ( 5 km ) using a gravity anomaly grid in some parts of the Chad basin (Termit and Bornu basin).
5. Conclusions
We applied gravity polynomial fitting and spectral analysis to observed residual anomalies from the northern Logone Birni sedimentary basin of the far north region of Cameroon, to estimate the structure of the basin as well the depth of residual anomalies. The residual anomaly maps were obtained by fitting the Bouguer anomalies using polynomial surfaces of first, second and third orders
Figure 5. Power spectrum of profiles P1, P2 and P3 associated to the second degree trend gravity map.
Figure 6. Power spectrum of profiles P1, P2 and P3 associated to the third degree trend gravity map.
and 3 profiles were extracted to estimate the residual anomaly depths. The shape of the residual anomalies were found to be almost elliptical and trend N-S, NW-SE, NE-SW and are accompanied by a steep dipping gradient. Interpretations of the results show that the structural style of the Logone Birni Basin is mainly NE-SW.
Three distinct residual anomalies were identified in the basin. Over the Goulfey, the deepest anomaly is located at a depth of about 4.55, 4.00 and 1.63 km , and the shallowest anomaly at 0.81, 0.47, 0.25 km . Over the Tom Merifine, the deepest anomaly is located at a depth of about 2.10, 2.09 and 1.67 km , and the shallowest anomaly at 0.34, 0.34, 0.24 km . Over the Tourba, the deepest anomaly is at a depth of about 3.26, 2.96 and 3.84 km , and the shallowest anomaly at 0.73, 0.49, 0.49 km . The average residual depth along the basin varies between 0.24 km and 4.55 km with the minimum at Tom Merifine and the maximum at Goulfey.