Azerbaijan is located in the Alpine Himalayan collisional zone and is characterized by its complex and variable geological structure. To study Azerbaijan’s deep structure, twelve main regional interpreting profiles were selected for comprehensive 3D combined gravity-magnetic modeling. The development of the initial physical-geological models (PGMs) was based on known surface geology, drilling data, previous seismic, magnetotelluric and thermal data analysis, examination of the richest petrophysical data, as well as quantitative and qualitative gravity/magnetic data examination. The PGMs thus reflect the key structural-formational specifics of Azerbaijan’s geological structure, beginning from the subsurface (hundreds of meters) up to the Moho discontinuity (40 - 60 km). The PGMs revealed common factors controlling ore- and hydrocarbon bearing formations, primarily the boundaries of tectonic blocks, masked faults, and buried uplifts of magmatic rocks. Many of these factors can be used to investigate long-term geodynamic activity at a depth. The article summarizes many years of investigation by exemplifying the most typical PGMs for the Greater and Lesser Caucasus, the Kura depression, and central and northern Azerbaijan.
Investigating the South Caspian basin to localize commercial oil-and-gas deposits is of great interest to geoscientists [e.g., 1-6]. At the same time, Azerbaijan is known to be rich in economic minerals. Besides oil and gas there are iron, copper, gold, mercury, zinc, molybdenum, cobalt, aluminum and other deposits [
The complexity of Azerbaijan’s geological structure stems from its location in the Alpine Himalayan collisional zone. The NE part of Azerbaijan is part of the Pre-Caucasian foreland filled by Cenozoic terrigenous sediments. A heterogenic Nakhichevan folding system is located in the SW part, where carbonate Paleozoic strata and Cenozoic magmatic formations are mixed (
According to Khain [
The Alpine tectono-magmatic cycle is characterized by more complete geological data. As a whole, the Azerbaijan territory is typical of frequently changing geological associations on the vertical and lateral axes,
the presence of multifarious fold and fault structures of different orders, and regional and local metamorphism. All these factors make the development of reliable models of these media highly complex.
Gravity-magnetic data processing is generally intended to reduce and eliminate noise factors of different origins and intensities. The main problem faced by qualitative interpretation is to single out a desired target, whereas quantitative interpretation needs to determine and refine the target parameters. Thus, geological problems need to be resolved in terms of: (1) the capabilities of the geophysical method selected for measurements of the field containing the information required, (2) the properties of the medium under study, its capability to generate detectable signals (anomalies), (3) the methods for data processing and interpretation; namely, their ability to extract information from the field revealing the effects of the geological objects.
The complex geological structure of Azerbaijan determines the highly intricate nature of the gravity and magnetic fields, which reflect the superimposed effects of outcropped bodies and structures and those buried at various depths. Therefore to identify anomalies in the different classes of geological sources, not only observed geophysical fields need to be used but also their transformations. Such a classification can be carried out by the use of spatial features (regional and local anomalies), morphological features (isometric anomalies, elongated anomalies or ledges), or the sign and intensity of geophysical fields.
To study the gravity anomalies in Azerbaijan, the results of Bouguer gravity analytical continuation to levels of 4, 8 and 20 km were examined. on the basis of these transformations, difference DgB(0-4), DgB(4-10) and DgB(8-20) maps were computed. It was found that the residual field DgB(0-4) in orogens reflects the influence of the near-surface Alpine structures composed of volcanogenic and volcanogenic-sedimentary associations and intrusions; in the Kura depression this field reflects the distribution of dense heterogeneities in the upper part of the
sedimentary cover. The difference field DgB(8-20) in orogens primarily reflects the influence of the structure and composition of the metamorphic basement (Baikalian folded foundation) and deep parts of some intrusions; in the Kura depression it mainly reflects the influence of the Mesozoic complexes. The pattern of the difference field DgB(4-10) is close to DgB(8-20), but is more differentiated.
To examine the regional magnetic anomalies the following magnetic data were used: field DT (in the Middle Kura Depression – field Za) at the levels of 6 and 10 km, and a field of DT (Za) horizontal gradient at the level of 2 km. Attempts to reveal local magnetic anomalies as a residual field of DT(0-6) or DT(0-10) were unsuccessful because of sharp differentiation of the observed magnetic field in the Lesser Caucasus and in some other areas.
The potential field characterization procedure was based on the richest data from petrophysical studies (more than 80 000 determinations of magnetic susceptibility and more than 25 000 determinations of density) conducted in different areas of Azerbaijan [
The geomagnetic section of the Greater Caucasus mega-anticlinorium is characterized on the whole by low values of magnetic susceptibility (= 80×10−6 SI), which sharply increase for volcanogenic and volcanogenic-sedimentary associations (
Quaternary (coarse molassa) with = 100×10-6 SI, (2) the Eocene-Miocene (thin molassa) with = 100 ×10−6 SI, and (3) the Jurassic-Cretaceous with = 30×10−6 SI.
In the density section, four floors could be identified (
The petromagnetic section is characterized by high differentiation: magnetic and strongly magnetic rocks predominate here (
The petrodensity section has a weak differentiation that complicates singling out accurate density boundaries (
Three petromagnetic floors were found: (1) the Pliocene-Quaternary (=300·10−6 SI), the Eocene-Miocene
(= 100·10−6 SI), and the Cretaceous (volcanogenic and volcanogenic-sedimentary associations) with= 1500·10−6 SI.
Low-density deposits compose most of the section of the Kura depression; only the Jurassic-Cretaceous associations discovered in the lower part of the section relate to the middle-density associations. The three petrodensity floors have the following characteristics: (1) the Quaternary (= 2.00 g/cm3), (2) the Paleogene – Neogene (= 2.18 g/cm3), and (3) the Jurassic – Cretaceous (= 2.62 g/cm3).
Two petromagnetic floors were found in the Talysh anticlinorium: the Eocene-Neogene (= 300×10−6 SI), and (2) the Paleocene - Middle Eocene (= 700×10−6 SI). In the second floor, inverse high-magnetized zones were identified.
The density values range considerably: (1) Pliocene (= 2.00 g/cm3), (2) Oligocene-Miocene (= 2.25 g/cm3), (3) Eocene (= 2.53 g/cm3), and (4) Paleocene (= 2.45 g/cm3).
The specificity shared by all these sequences is the availability of two geostructures petrodensity boundaries typical of all (
associations (the contrast Dsaver consists of 0.14 - 0.44 g/cm3). However, in the latter, the boundary is not differentiated clearly in the sequence of the Lesser Caucasus: Dsaver reaches only 0.05 g/cm3. Here the reference horizon of lower density and magnetization are liparites and plagio-liparites of the Upper Bajocian.
As a whole, the densest rocks are typical of intrusive bodies. The density of intrusive rocks of the Early Alpine sub-stage ranges approximately from 2.80 to 2.95 g/cm3(excluding Bajocian plagiogranites with a density of 2.50 - 2.60 g/cm3). The density of the Late Alpine intrusives ranges on average from 2.70 to 2.85 g/cm3; the density of pre-orogenic intrusives exhibits a clear tendency to decrease to the class of intermediate density (saver = 2.65 - 2.75 g/cm3). The intrusives of the orogenic sub-stage have the lowest density (saver = 2.55 - 2.60 g/cm3).
The rocks composing effusive associations mainly belong to the class of increased (saver = 2.70 - 2.85 g/cm3) and intermediate (saver = 2.50 - 2.70 g/cm3) density. This increased density is typical basically of effusive rocks of the Early Alpine sub-stage. The density of the later sub-stages decreases gradually, and rocks of the orogenic sub-stage almost completely belong to the class of low density (saver = 2.45 - 2.60 g/cm3).
Among the sedimentary deposits, the sandy-shales of the Greater Caucasus, as well as carbonates of the Late Jurassic-Early Cretaceous are the densest (saver = 2.60 - 2.75 g/cm3). The tendency for density to decrease from the early to late stages is also found here. The main petrodensity discontinuities show similar peculiarities in the central and western regions of the Caucasus, namely in Armenia [
As a first approximation, the geomagnetic sequence can be divided into three floors corresponding to the main tectonic stages. The lower floor (Pre-Jurassic), represented mainly by terrigeneous-carbonate formations penetrated by acid intrusions, is practically non-magnetic. Some exclusions can be found in the volcanogenic associations revealed in the Shamkhor anticlinorium of the Lesser Caucasus. The middle floor (Mesozoic-Eocene), within which thick volcanites of basic and intermediate composition have developed, is differentiated by its increased magnetization. The upper floor (Post-Eocene) represented mainly by the terrigeneous deposits of the sedimentary cover is characterized by a lower magnetization (although in some areas Neogene-Quaternary lavas with high magnetization are known).
In the intermediate floor volcanogenic associations of the Middle Jurassic have the highest average magnetic susceptibility – in the Lesser Caucasus they are approximately 4 km thick. The remanent magnetization of these rocks may be parallel or antiparallel to the induced magnetization; relatively low-magnetic associations are usually typical of the zones of inverse magnetization.
The Cenozoic volcanogenic associations of the Talysh as well as lava covering the central part of the Lesser Caucasus have other characteristics. Here paleomagnetic zones of different polarities alternate. Analyses indicate that upward continued magnetic anomalies (associated with the gently sloping Cenozoic associations) are rapidly attenuated. The lavas are characterized by significant magnetization values compared to pyroclastic associations of a similar composition. The magnetization of ultrabasites is approximately proportional to their degree of serpentinization (at the same time the density of the ultrabasites decreases as a function of their serpentinization).
By analyzing the sequence of the main Caucasian petrodensity and petromagnetic floors (
To conduct a quantitative interpretation of magnetic and gravity anomalies, advanced methods have specially been developed for complex environments [
One example of such an interpretation is the Guton magnetic anomaly situated in NW Azerbaijan, near the border with Russia (
abovementioned characteristics point to the basic-acid composition of the anomalous body (intrusion). The conclusion as to the significant vertical thickness of this body concurs with the geothermic data for the depth of the Curie discontinuity in this area (about 30 km) [
The development of 3-D PGM was performed in three
basic steps:
The first step is as follows:
(A) A geological Section is constructed where all the intrusive, effusive and other associations, as well as faults and the surface of folded foundation are compiled on the basis of geological data within a strip 15 - 20 km wide. The interpreting Section is located in the middle of this strip.
Such a Section characterizes the upper portion of the Earth’s crust, which is 2 - 3 to 5 - 8 km in thickness from the Earth’s surface to the Baikalian basement. Deeper parts of the intrusive bodies and certain faults are formed by extrapolation of the available constructions, general geological considerations and the results of previous geophysical analyses.
(B) A preliminary petrophysical model of the Section is developed. Here all the geological bodies acquire density and magnetization values according to the preceding petrophysical data and results of geophysical field interpretations. When no data are available on the magnetization direction, it is assumed to be parallel to the normal geomagnetic field. Further, the magnetization direction is refined in the course of physical-geological modeling. The petrophysical model includes deep-seated layers of the Earth’s crust: (1) the “basaltic”, (2) the intermediate between the crust and the upper mantle, and (3) the upper mantle. Their surfaces are constructed and physical properties are associated with them according to the data from previous seismic, magnetotelluric and other deep geophysical studies.
(C) The initial (preliminary) petrophysical model includes hidden bodies as well. Their location, thickness, depth, density and magnetization are obtained from a quantitative analysis of magnetic and gravity anomalies as well as seismic data.
The second step selects the gravity and magnetic fields along the interpreting profiles using the Geological Space Field Computing (GSFC) program [19,35]. Each time the fields are from different bodies, groups of bodies and the total computed model are displayed and compared to the observed gravity and magnetic fields. Using the results of this comparison, the changes that match the gravity and magnetic effects are introduced into the model of the medium. Computations, comparisons of fields and model modifications are repeated until the desired fit between the computed and observed fields is obtained.
Then, a regional gravity (and sometimes magnetic) field is roughly selected. As a rule, the densities of deep-seated complexes are not changed; the modi-fications only affect the shape of their roof. Next, fields of local bodies are selected. If necessary, this is followed by a verification of the regional field and the field of the local bodies.
At each computational step, a separate analysis of gravity and magnetic fields is carried out. Coordinated variations are verified in the subsequent steps, and then introduced into the model. This procedure leads to an integrated quantitative interpretation for anomalous gravity and magnetic fields. The selection ends when the computed gravity and magnetic fields coincide accurately with the observed fields.
The third step involves a detailed geological interpretation of these models. A 3-D PGM of the area under investigation is developed based on the qualitative and quantitative data. This yields the final geological maps, and the models are characterized by a more complete rendering of the geological targets, mineral controls and mineral deposits, including deep-seated ones.
The geological interpretation of the complexes and local bodies of the selected (final) petrophysical model does not usually present any difficulties, since in the implementation of the interactive selection system almost all the bodies in the model acquire some specific geological content. The geological nature of new sources introduced into the model during the selection and reflected either in the initial geological section, or in the initial PGM, is determined according to the similarity of their physical properties, dimensions, and depth of occurrence with respect to the known targets. The age of the bodies is determined according to their interrelations with the host rocks.
Profiles 1 and 2 crossing the Lesser Caucasus (Figures 8 and 9) illustrate the very complex geological structure of this region. The Late-Alpine effusives in these PGMs compose an ophiolite zone (which is a relic of the ocean crust). It is thought that the same rocks occur in the NE immersion of the Lesser Caucasus. Pre-orogenic and orogenic intrusive and effusive rocks are fixed in the southern parts of these PGMs. Thick sedimentary deposits are found in northern parts of these profiles. A smooth high of the Moho discontinuity is observed from south to north from a depth of 52 - 54 km up to 42 km. Classes of disturbing objects were revealed here such as acid intrusions of lower density and magnetization, basic magmatic rocks of increased density and magnetization and fault zones. It was determined that the clearest density boundaries were associated with the base of the Cenozoic sedimentary strata and to a lesser degree with the base of the Alpine complexes. According to the modeling, geomagnetic boundaries were associated mainly with the roof and bottom of the Mesozoic floor of heightened magnetization.
In the southern part of profile 2, a zone of strong decompaction is observed in the root parts of the Dalidag intrusive at depths from 4 to 17 km, and by a horizontal thickness of 30 km. This rock decompaction zone produces the greatest Kelbadzhar-Dalidag gravity minimum. This large “granite room” apparently was a source of pre-orogenic and orogenic granitoid magmatism.
Three profiles (No. 4, 5 and 6) were studied across the Greater Caucasus mega-anticlinorium. Profile 6, presented in
depth of the Moho discontinuity from the south to the middle part of this profile (from 47 to 60 km). The behavior of the gravity field depends primarily on the Moho depth.
In this PGM the Early-Alpine cores of the Greater Caucasus are clearly linked to the development of Jurassic associations present in the Lower-Middle-Jurassic sandy-shale deposits. A surface of the Pre-Alpine foundation occupies the highest position in this core. In the PGM (the Greater Caucasus) their wide distribution was predicted given the results of 3D combined gravitymagnetic modeling of the Early-Alpine LowerMiddle-Jurassic sandy-shale deposits of the Greater Caucasus (occupying this position as a result of overthrust), as well as the sedimentary cover of the Kura Depression. In the northern parts of this profile, rocks of significant magnetization were identified at a depth of 9 - 10 km. These rocks are presumably related to tuff-terrigeneous deposits of the Middle-Late-Paleozoic age.
For many years the prevailing view was that in the Kura depression separating mega-anticlinoria of the Greater and the Lesser Caucasus, thick sedimentary deposits occur in the crystalline Pre-Alpine basement, and subvertical deep faults divide these structures. On a projection of the buried uplift of the basement to the Earth’s surface, which is assumed to exist on the basis of high densities and velocities of elastic waves, the Saatly Super-Deep borehole (SD-1) was drilled in 1965. However, analysis of the magnetic properties of rocks and magnetic survey results showed that the basement was not magnetized, and the Mesozoic magmatic associations of basic and intermediate composition (of high magnetization) occupied most of the geological section in the Middle Kura
depression [
The validity of this interpretation was confirmed by the results of SD-1 drilling: the borehole exposed Mesozoic volcanogenic rocks at a depth of 3.6 km and did not come from the rocks at its bottom at 8.2 km [
An integrated PGM (gravity-magnetic-seismic with elements of thermal field analysis) of the Earth’s crust along profile No. 9 is shown in
The analysis indicates that the area of the TalyshVandam gravity maximum (TVGM) has a highly inhomogeneous geological structure. Separate elements of the TVGM – anomalies of the second order – reflect differently, showing that they originated independently in different areas of the Earth’s crust. Taking into account the location of the lower edge of the magnetized masses and similarities between this geological section and the Lesser Caucasian, it was presumed that magmatic rocks should be found down to 10 km in depth (comparable to the Lower Bajocian rocks of the Lesser Caucasus).
Obviously, the main source of the TVGM is associated with the underlying highly dense, strongly metamorphosized (initially chiefly sedimentary associations) nonmagnetic or low-magnetic Pre-Baikalian floor (this type of floor has a submeridional strike in Russian and African platforms [29,36]). The depth of the upper edge of these highly dense rocks was estimated at 9.5 km.
Thus, SD-1 has yet to discover (it is unlikely that these drilling operations will be pursued in the coming years) the source of the TVGM, but clarified the origin of the Gandja magnetic maximum. This is crucial not only for analysis of the tectonic-magmatic evolution of the Caucasus region, but also for prospecting of its oreand oil and gasbearing deposits. For example, the geomag-
netic models makes it possible to extend prospecting of oil deposits in the Middle Kura depression, which are associated with the zones of protrusions of the Mesozoic magnetoiactive associations. Many geophysical methods have mapped the dense Mesozoic associations, but only magnetic prospecting has been able to reveal the geological associations of basic and middle consistency [
The complex geological structure of Azerbaijan determines the highly intricate nature of its gravity and magnetic fields, which reflect the superimposed effects of outcropped bodies and structures and those buried at various depths. It was shown that a combined gravity-magnetic investigation (including the field behaviors at different levels, qualitative and modern quantitative analyses, and 3D combined modeling of potential fields) is powerful tool for studying such complex geological regions as the land of Azerbaijan. A significant role is played by the analysis of petrophysical characteristics, as well as an integrated examination of seismic, magnetotelluric and thermal data. On the basis of a comprehend-sive combined analysis of gravity and magnetic fields, a series of 3-D regional PGMs (including the area of the Saatly Super-Deep borehole), covering Azerbaijan and some adjacent land areas, were developed. The key examples of PGMs presented here highlight the significant formational and tectonic peculiarities of Azerbaijan’s deep structure and could be used for various types of geological-geophysical zonation. The next steps should include a more detailed physical-geological examination of these PGMs.
The author would like to thank the anonymous reviewers who thoroughly reviewed the manuscript, and whose critical comments and valuable suggestions were so helpful in preparing this paper.
This article is dedicated to the memory of Prof. Boris Khesin who passed away on November 8, 2010 in his 79th year.
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[