A geostatistical re-interpretation of gravity surveys in the Yagoua, Cameroon region
Philippe Njandjock Nouck1*, Chamberlin Kenfack1, Ahmad Diab Diab2, Kasi Njeudjang3, Larissa Jorelle Meli'I4 and Rodrigue Kamseu5
1 Department of Physics Faculty of Science University of Yaoundé I Cameroon. *Corresponding author:
2 Department of Physics Faculty of Science University of Abeché. Chad.
3 Pecten oil and gas Company Douala, Cameroon.
]]> 4 Departments of Physics Advanced Teachers Training College of Nkolnda-Nsimalen, Cameroon.5 Centre de Formation spécialisée de Géostatistique MINES-Paris France.
Received: September 28, 2012.
Accepted: March 13, 2013.
Published on line: September 30, 2013.
Resumen
Desde 1960 se han realizado diversos estudios de gravedad en la región Yagoua del norte de Camerún. Se recabaron datos de gravedad en una área amplia que abarca aproximadamente 11.628 km2. Estos datos son insuficientes, irregulares, dispersos y no permiten eficientemente continuaciones ascendentes y descendentes del campo gravitatorio, derivadas y otras operaciones que requieren datos reticulados regulares. Algunas anomalías en el mapa Collignon (1968) pueden correlacionarse con la estructura geológica que se conoce, pero no aparecen en los mapas de Louis (1970) y Poudjom et al. (1996). Para producir los datos de gravedad reticulares regulares y mejor control de las anomalías, derivadas de estructuras geológicas, se aplicó el método de Kriging a una línea de base de datos-188. Se ensayaron para este propósito varios modelos de variograma. Se encontró que un modelo esférico era la mejor opción; se ha elaborado un nuevo conjunto de datos Kriging con unos 10.100 resultados y un nuevo mapa con los datos Kriged Bouguer. Este mapa contiene anomalías positivas en las zonas Maroua-Mindif y Maga (1968) en el mapa Collignon, que no estaban presentes en los mapas de Louis (1970) y Poudjom et al. (1996). Las anomalías positivas de Guibi-Doukoula y Yagoua, que no se encuentran separados en los mapas de Louis (1970) y Poudjom et al. (1996), aparecen claramente distintas a como fueron previstas por Collignon (1968). Los nuevos resultados pueden ser utilizados para los estudios gravimétricos posteriores.
Palabras clave: Anomalía Bouguer, datos de gravedad, geoestadística, variograma, Kriging, Yagoua.
]]>Abstract
Since 1960, many gravity studies have been carried out in the Yagoua region of northern Cameroon. Gravity data was collected over a wide area of approximately 11628 km2. These data are insufficient, irregular, scattered and do not efficiently permit gravity field downward and upward continuations, derivatives and other operations that might require regular gridded data. Some anomalies on the Collignon map (1968), may correlate with known geological structure but do not appear on maps by Louis (1970) and Poudjom et al. (1996). To produce regular gridded gravity data and better control anomalies due to geological structures, the kriging method was applied to a 188-data baseline. Several variogram models were tested for this purpose. It was found that a spherical variogram model is the best; it has produced a new kriging dataset of about 10,100 data and a new map of kriged Bouguer data. This map contains positive anomalies in the Maroua-Mindif and Maga areas on the Collignon (1968) map, which were not present on Louis (1970) and Poudjom et al. (1996) maps. The positive anomalies of Guibi-Doukoula and Yagoua, not separated on the Louis (1970) and Poudjom et al. (1996) maps, show up as clearly distinct as previewed by Collignon (1968).The new results can be used for subsequent gravimetric studies.
Key words: Bouguer anomaly, gravity data, geostatistics, variogram, kriging, Yagoua.
Introduction
Geostatistics is applied in the Earth sciences as an interpolation procedure that uses an available dataset to obtain an optimal, linear and unbiased estimation of a property whose estimation error is minimized (Matheron, 1973). It operates on a random variable for which a set of possible values is known but whose final result requires a measurement. The objective of this study is to re-examine a gravity data set from irregular gridded data. For this operation known as kriging, we make a crucial choice of a variogram. This is a better choice, because more kriged values are closer to reality. The principal difficulty is to find a variogram model that fits the data to be interpolated. In this work, we calculate an experimental variogram from existing gravity data, we evaluate the RMS between the variogram and we use various cross-data validation criteria to ensure that the variogram model chosen is optimal in order to correlate with the experimental variogram. We also offer a new gravity dataset and we propose a new gravity map of the Yagoua region in northern Cameroon.
Study area and Gravity data
Presentation of the study area
]]> Yagoua region (Figure 1), as others in Central-Africa, is the product of a complex period of continental disruption associated to plate tectonic fragmentation of Gondwana (Genik, 1992; Njandjock, 2004). It is characterised by polyphase rifting, separated by tectonic events that can be linked to regional deformation, hiatus of sedimentation and unconformities in seismic sections and outcrops. The region belonging to the Panafrican belt, is bounded to the South by the Doba basin, to the North by Lake Chad basin, to the East by Doseo and Salamat basins, and to the North West by the Mandara Mountains and the Southeast by the Kaele dome. The study area in Cameroon, covers an area of about 11628 km2 and stretches from latitude 9°45 'N to 10°47'N to longitude 14°20' E to 15°30'E.
Geological context
The Yagoua region is located in the southern part of the Logone Birni Basin (LBB) characterised by Quaternary sediments and belonging to the West and Central African Rift System. The geology of the region (Figure 1) is underlain by a large sedimentary formation and a Precambrian basement. The basement consists of acid and metamorphic formations. gneisses, migmatites, diorites, anatexites, syenites, syn-tectonic to post-tectonic granites, basalts and shales. Rare basalts are observed in the Kaélé region (Maurin, 2002). Among the bedrock outcrops, the Precambrian is dominant. It is characterized by migmatites and anatexites located southwest of the region in Maroua and Kaélé. Some remarkable features include the Mindif intrusion composed of syenites (Eno Belinga, 1984), the North Maroua gabbro hills and the Kaélé hills. The sedimentation in this area commenced during the Neocomian-Albian rifting period and consists of sandstones, clays and shales. The sedimentary formation is everywhere overlain by sandy alluvial formations and dunes attributable to the Pre-Bima formation (Louis, 1970; Genik, 1992).
Gravity data
Gravity data used in this study were acquired during surveys conducted by various organizations and researchers: Collignon (1968); Louis (1970); Elf-Serepca, 1980; Poudjoum et al., 1995-1996; Njandjock, 2002 to 2004 and Njandjock et al., 2006. Measurements were carried out along profiles to detect significant variations of geological facies. The station locations were obtained from topographic maps and compass traverses. The elevation of the stations was obtained from barometric readings, using Wallace and Tierman or Thomnen altimeters. Variations in the gravity field were measured using Worden, Lacoste & Romberg and North American gravimeters. The gravity data were converted to Cartesian coordinates and need to be interpolated. The conversion to kilometric distances was by UTM (Universal Transverse Mercator) on the Clarke's ellipsoid (1880) with the Prime Meridian. Bouguer anomaly maps by Collignon(1968), Louis (1970) and Poudjom (1996) are shwon on figures 2, 3, 4 and 5.
]]> Method
The main purpose of a geostatistical study is to construct a mathematical model of the random function Z (X), based on an experimental data set Zexp (Xi), a single realization of Z (X), where x and x + h are two points separated by a distance h, and z (x) and z(x + h) are associated random variables. The variogram is known and the semi-variance of the difference [z (x)-z (x + h)] is (Chiles and Delfiner, 1999; Kumar and Remadevi 2006):
Experimentally, it is calculated by the equation
Supposing that we want to estimate a block V centered at X0. ZV denote the true value (unknown) of this block and ZV* the estimate that is obtained:
Zi being the random variables corresponding to sample points. We want to minimize
Substituting the expression of the estimated in this equation, we obtained
]]>
Since the estimate is unbiased, we set
To minimize σ2e, we use Lagrange's method and we form the Lagrangian
Where µ is the Lagrange multiplier. The minimum is reached when all the λi and µ partial derivatives are canceled out. This leads to the following ordinary kriging system:
Minimum variance estimation called kriging variance is obtained by substituting the kriging equations in the general expression of the variance estimation.
Results and Discussion
Experimental variogram of the Yagoua region
The experimental variogram is now calculated from equation (2). Results are ranked by increasing distance, and then grouped in intervals centered on increasing multiples of h= 5. The results shown in Table 1, are used to plot the variogram (Figure 6).
The variogram characteristics are: range a = 25; level C + C0 = 81; nugget; C0 = 1.
]]> The known mathematical variogram models in the literature were used to calculate different models of variograms and compare them to the experimental one (Figures 6 and 7). The results are summarized in table 2. In these calculations, several models (for example the Gaussian model) not shown here have been eliminated. The superposition of curves giving the variogram as a function of h is given in Figure 7. The dark blue curve represents the experimental variogram of the Yagoua region. Our aim is to find among the other curves of this graph the one that best approximates the dark blue curve. The outer gray curve represents the gravity model whose elimination is visible; it is difficult to comment on the other curves. Indeed within this paragraph, it is difficult to say exactly which curve (red, green, purple, blue or orange) is closest to the dark blue curve. Qualitatively all these models are eligible. For this reason we will make a quantitative analysis using the equation
to calculate the relative changes, and standard deviation of these relative changes between each variogram model and the experimental variogram.
Changes between experimental variogram values and those of different models are reported in table 3. The standard deviation values between the experimental variogram and the different model is 2.0 for spherical, 3.3 for quadratic, 5.2 for pentaspherical, 7.5 for exponential and 13.2 for Gaussian (table 3). It is clear that the model giving the best variogram model are the spherical, tracking quadratic and pentaspherical models. In order to confirm or invalidate this result, we made kriged maps using several models. The goal is to determine which variogram model reproduces the experimental values. For this purpose, we have omitted 15 of the 188 baseline. Then we realized the kriged map of spherical, square, pentaspherical, exponential, cubic, Gaussian, rational quadratic and power models with the 173 remaining data. On these kriged maps, we reported values of anomalies the hidden points. Abnormal values obtained by kriging with different variogram models are shown in table 4. The relative differences between the anomalies from various models and the current anomalies are contained in table 5.
The lowest standard deviation 2.69 is obtained with the spherical model, followed by 2.86 for the exponential model and 2.87 for the pentaspherical model. The spherical model is the one that is best reproduced by kriged data or field values. In both approaches, the first four eligible models are the same (spherical, pentaspherical, quadratic and exponential) but in a different ranking than the second, the first being the spherical model. The differences between the values of standard deviation from one model to another are not great. This was predictable because curves in Figure 7 are all very close to the experimental variogram curve.
A new Bouguer map of Yagoua region
The new Bouguer anomaly map (Figure 8) obtained by the kriging method reveals several anomalies which can be grouped into two major families: a family of negative anomalies and a family of positive anomalies. The positive anomalies are located south of Kaele, at Guibi-Doukoula, Yagoua, Guidiguis and at Maroua-Mindif. The positive anomalies of Guibi-Doukoula and Yagoua, observed on the new map, also appear in other documents. They appeared separated on Collignon (1968) but are absent on Louis (1970) and Poudjom et al. (1996) maps. The reappear clearly separated on the new map proposed by the current study. In addition, this map brings out the positive anomalies of Maga and Maroua-Mindif formerly reported by Collignon (1968) and absent on Louis (1970) and Poudjom et al. (1996) maps. By comparing these results (Figure 9) with the geological map (Figure 1) and, according to Manga et al. (2001), Njandjock (2004) and Njandjock et al. (2006), the positive anomalies of Guibi-Doukoula, Yagoua and Guididuis correlate with an uplift of gneissic-granite type dense materials. The positive anomaly of Maroua-Mindif correlates with Mindif's "tooth" syenite block. The negative anomaly directed NE-SW extends from west Kaele to north Moulvouday and occupies the center of the study area. This negative anomaly is surrounded by positive anomalies and may coincide with a Quaternary sedimentary cover at Moulvouday. Compared to Manga et al. (2001) and Njandjock et al. (2006), this large anomaly may be the gravity signature of the Moulvouday-Yagoua sedimentary basin.
Conclusion
]]> In this study, the experimental variogram was calculated from the existing gravity data. The RMS and the hindered points technique showed that the best theoretical variogram for these data were spherical. This variogram was therefore used to produce a new regular gravity data base and a new gravity map of the Yagoua region which can be used in further researches. This study includes mapping of the Moulvouday negative anomaly which may correlate with sedimentary thickening, while positive anomalies of Guibi-Doukoula and Yagoua may correspond to granite-gneiss dense basement material uplift. In addition, we confirm the positive anomalies of Maga and Maroua-Mindif which no longer appeared on Louis (1970) and Poudjom et al. (1996) but were reported by Collignon (1968). In addition, positive anomalies of Guibi-Doukoula and Yagoua that appeared separated on the Collignon map (1968), and which were present on the Louis (1970) and Poudjom et al. (1996) maps, reappear clearly separated.
Acknowledgements
We would like to thank Dr. Cinna Lomnitz and Mister ACHAKENG NKEMKA John of Ministry of Scientific Research and Innovation-Cameroon, for corrections and suggestions on the manuscript.
Bibliography
Bessoles B., Trompette R., 1980, The Pan-African chain, mobile zone of Central Africa mobile zone and sudden, BRGM. Memory 92. [ Links ]
Chiles and Delfiner, 1999, Geostatistics, Modeling spatial uncertainty. Wiley series on probability and statistics. [ Links ]
]]>Collignon F., 1968, Gravimétrie de reconnaissance du Cameroun. Orstom, 37p. [ Links ]
Elf-Serepca, 1981, Geological Map of Garoua basin's, Operating service. [ Links ]
Eno Belinga S.M, 1984, Geology of Cameroon. University Bookstore in Yaoundé, p.307. [ Links ]
Genik G.J, 1992, Regional structural gravimetry and petroleum aspects of basins in Niger, Chad and Central African Republic (CAR), Tectonophysic, 313, 169-185. [ Links ]
Kumar V., Remadevi, 2006, Kriging of Ground-water Levels – A case Study. Journal of Spatial Hydrology, 6, 1, 81-94. [ Links ]
]]>Louis P., 1970, Contribution géophysique à la connaissance géologique du bassin du lac Tchad. Orstom, 42, Paris 311p. [ Links ]
Manga S.C., Loule J.P., Koum J.J., 2001, Tectonostratigraphic evolution and prospectivity of Logone Birni Basin, North-Cameroon-Central Africa. AAPG Extended Abstracts, 124. [ Links ]
Matheron G., 1973, The intrinsic random functions and their applications. Advances in applied probability, 5, 439-468p. [ Links ]
Maurin S., 2002, Geomorphology, atlas of the province of the far North Cameroon. Plate 04, MINRESI-INC-Cameroon. [ Links ]
Njandjock Nouck P., 2004, Gravimetric investigation and mapping of geological structures in the superficial and deep-Cameroon region Yagoua Ph.D Thesis in Physics. University of Yaounde I. [ Links ]
]]>Njandjock Nouck P., 2007, Note on geostatistical courses, University of Yaounde I. [ Links ]
Poudjom-Djomani Y.H., Legeley-Padovani A., Boukele D.B., Nnange J.M., Ateba-Bekoa, Alboul Y., Fairhead J.D., 1996, Levés gravimétriques de reconnaissance-Cameroun. Orstom, France, 30p. [ Links ]
Vicat J.P., Bilong P., 1998, Geology and environment in Cameroon. University Press of Yaounde. [ Links ]
]]>
