Crustal structure of the western part of the Southern Granulite Terrain of Indian Peninsular Shield derived from gravity data

The Southern Granulite Terrain (SGT) is composed of high-grade granulite domain occurring to the south of Dharwar Craton (DC). The structural units of SGT show a marked change in the structural trend from the dominant north–south in DC to east–west trend in SGT and primarily consist of different crustal blocks divided by major shear zones. The Bouguer anomaly map prepared based on nearly 3900 gravity observations shows that the anomalies are predominantly negative and vary between −125 mGal and +22 mGal. The trends of the anomalies follow structural grain of the terrain and exhibit considerable variations within the charnockite bodies. Two-dimensional wavelength filtering as well as Zero Free-air based (ZFb) analysis of the Geoid-Corrected Bouguer Anomaly map of the region is found to be very useful in preparing regional gravity anomaly map and inversion of this map gave rise to crustal thicknesses of 37–44 km in the SGT. Crustal density structure along four regional gravity profiles cutting across major shear zones, lineaments, plateaus and other important geological structures bring out the following structural information. The Bavali Shear Zone extending at least up to 10 km depth is manifested as a plane separating two contrasting upper crustal blocks on both sides and the gravity high north of it reveals the presence of a high density mass at the base of the crust below Coorg. The steepness of the Moyar and Bhavani shears on either side of Nilgiri plateau indicates uplift of the plateau due to block faulting with a high density mass at the crustal base. The Bhavani Shear Zone is manifested as a steep southerly dipping plane extending to deeper levels along which alkaline and granite rocks intruded into the top crustal layer. The gravity high over Palghat gap is due to the upwarping of Moho by 1–2 km with the presence of a high density mass at intermediate crustal levels. The gravity low in Periyar plateau is due to the granite emplacement, mid-crustal interface and the thicker crust. The feeble gravity signature across the Achankovil shear characterized by sharp velocity contrast indicates that the shear is not a superficial structure but a crustal scale zone of deformation reaching up to mid-crustal level.     

Multiple episodes of rifting in Central and East Africa: A re-evaluation of gravity data

A compilation of new and existing gravity data, as well as geophysical and geological data, is used to assess the cumulative effects of multiple rifting episodes on crustal and upper mantle density structures beneath the Uganda-Kenya-Ethiopia-Sudan border region. This compilation includes new gravity and geological data collected in 1990 in south-western Ethiopia. Variations in the trends and amplitudes of Bouguer gravity anomalies reveal three overlapping rift systems: Mesozoic, Paleogene and Miocene-Recent. Each of these rift systems is a number of 40–100 km long sedimentary basins, and each system is approximately 1000 km long. The Bouguer anomaly patterns indicate that the Ethiopian and East African plateaux and corresponding gravity anomalies are discrete tectonic features. Models of structural and gravity profiles of two basins (Omo and Chew Bahir basins) suggest that pre-Oligocene (Cretaceous?) strata underlie 3 km or more of Neogene-Recent strata within the northern Kenya rift, and that more than 2 km of Neogene-Recent strata underlie parts of the southern Main Ethiopian rift. The superposition of perhaps three rifting episodes in the Lake Turkana (Omo) region has led to 90% crustal thinning (β ≈ 2).     

Crustal structure and rift tectonics across the Cauvery–Palar basin,Eastern Continental Margin of India based on seismic and potential field modelling

The Cauvery–Palar basin is a major peri-cratonic rift basin located along the Eastern Continental Margin of India (ECMI) that had formed during the rift-drift events associated with the breakup of eastern Gondwanaland (mainly India–Sri Lanka–East Antarctica). In the present study, we carry out an integrated analysis of the potential field data across the basin to understand the crustal structure and the associated rift tectonics. The composite-magnetic anomaly map of the basin clearly shows the onshore-to-offshore structural continuity, and presence of several high-low trends related to either intrusive rocks or the faults. The Curie depth estimated from the spectral analysis of offshore magnetic anomaly data gave rise to 23 km in the offshore Cauvery–Palar basin. The 2D gravity and magnetic crustal models indicate several crustal blocks separated by major structures or faults, and the rift-related volcanic intrusive rocks that characterize the basin. The crustal models further reveal that the crust below southeast Indian shield margin is ～36 km thick and thins down to as much as 13–16 km in the Ocean Continent Transition (OCT) region and increases to around 19–21 km towards deep oceanic areas of the basin. The faulted Moho geometry with maximum stretching in the Cauvery basin indicates shearing or low angle rifting at the time of breakup between India–Sri Lanka and the East Antarctica. However, the additional stretching observed in the Cauvery basin region could be ascribed to the subsequent rifting of Sri Lanka from India. The abnormal thinning of crust at the OCT is interpreted as the probable zone of emplaced Proto-Oceanic Crust (POC) rocks during the breakup. The derived crustal structure along with other geophysical data further reiterates sheared nature of the southern part of the ECMI.     

Implications of gravity data from East Kalimantan and the Makassar Straits: a solution to the origin of the Makassar Straits?

Recent free-air gravity data covering the Makassar Straits is integrated with Bouguer gravity data from onshore East Kalimantan to provide new insights into the basement structure of the region. Onshore Kalimantan, gravity highs on the northern margin of the Kutai Basin trend NNE–SSW and N–S and correspond with the axes of inverted Eocene half-grabens. NW–SE trending lows correspond to deep seated basement weaknesses reactivated as normal faults during the Tertiary. An intra-basin gravity high trending NNE–SSW, the Kutai Lakes Gravity High, is modelled as folded high density Paleogene sediments flanked by syn-inversion synclines infilled with low density sediments. Offshore Kalimantan, the Makassar Straits include two basins offset by an en-echelon fault zone, suggestive of an extensional origin. The regional signature of the free-air anomaly data mirrors the bathymetry, but this effect can be reduced by the use of filters in order to examine the basin architecture. The free-air gravity minimum in the Makassar Strait is only −20 mGal, much smaller than that appropriate for a foreland basin, and more indicative of an extensional basin. The steepness of the gradients on the flanks of the basins indicates fault control of their margins. A regional 2D profile across the North Makassar Basin suggests the presence of attenuated crust (<14 km) in the basin axis at the present day, whereas flexural backstripping implies the presence of oceanic crust of middle Eocene age. The presence of oceanic crust in the North Makassar Straits Basin has implications for regional plate tectonic models.     

利用重力场研究东北地区断裂分布及构造分区

为了深入研究东北地区断裂分布及构造分区,利用全国1∶2 500 000布格重力异常图,分析了泛东北地区（包括东北3省及内蒙古自治区的东北部分地区）的区域重力场特征。利用布格异常变化形态,并参考地质资料,在研究区内划分出岩石圈断裂14条,壳内断裂11条。根据大兴安岭和依兰—伊通两条明显的重力异常梯级带把研究区划分出兴安、松嫩以及张广才岭等三大重力异常区,在划分断裂和分析重力场特征的基础上进行了构造分区。其研究成果为认识和研究泛东北地区油气与矿产分布提供了重力场方面的依据。     

Integrated 3D density modelling and segmentation of the Dead Sea Transform

A 3D interpretation of the newly compiled Bouguer anomaly in the area of the “Dead Sea Rift” is presented. A high-resolution 3D model constrained with the seismic results reveals the crustal thickness and density distribution beneath the Arava/Araba Valley (AV), the region between the Dead Sea and the Gulf of Aqaba/Elat. The Bouguer anomalies along the axial portion of the AV, as deduced from the modelling results, are mainly caused by deep-seated sedimentary basins (D > 10 km). An inferred zone of intrusion coincides with the maximum gravity anomaly on the eastern flank of the AV. The intrusion is displaced at different sectors along the NNW–SSE direction. The zone of maximum crustal thinning (depth 30 km) is attained in the western sector at the Mediterranean. The southeastern plateau, on the other hand, shows by far the largest crustal thickness of the region (38–42 km). Linked to the left lateral movement of approx. 105 km at the boundary between the African and Arabian plate, and constrained with recent seismic data, a small asymmetric topography of the Moho beneath the Dead Sea Transform (DST) was modelled. The thickness and density of the crust suggest that the AV is underlain by continental crust. The deep basins, the relatively large intrusion and the asymmetric topography of the Moho lead to the conclusion that a small-scale asthenospheric upwelling could be responsible for the thinning of the crust and subsequent creation of the Dead Sea basin during the left lateral movement. A clear segmentation along the strike of the DST was obtained by curvature analysis: the northern part in the neighbourhood of the Dead Sea is characterised by high curvature of the residual gravity field. Flexural rigidity calculations result in very low values of effective elastic lithospheric thickness (t _e < 5 km). This points to decoupling of crust in the Dead Sea area. In the central, AV the curvature is less pronounced and t _e increases to approximately 10 km. Curvature is high again in the southernmost part near the Aqaba region. Solutions of Euler deconvolution were visualised together with modelled density bodies and fit very well into the density model structures. An erratum to this article can be found at     

Satellite Gravity Anomalies and Their Correlation with the Major Tectonic Features in the South China Sea

The South China Sea (SCS) is a region of interaction among three major plates: the Pacific, Indo-Australian and Eurasian. The collision of the Indian subcontinent with the Eurasian plate in the northwest, back-arc spreading at the center, and subduction beneath the Philippine plate along Manila trench in the east and the collision along Palawan trough in the south have produced complex tectonic features within and along the SCS. This investigation examines the satellite-derived gravity anomalies of the SCS and compares them with major tectonic features of the area. A map of Bouguer gravity anomaly is derived in conjunction with available seafloor topography to investigate the crustal structure. The residual isostatic gravity anomaly is calculated assuming that the Cenozoic sedimentary load is isostatically compensated. The features in the gravity anomalies in general correlate remarkably well with the major geological features, including offsets in the seafloor spreading segments, major faults, basins, seamounts and other manifestations of magmatism and volcanism on the seafloor. They also correlate with the presumed location of continental-oceanic crust boundary. The region underlain by oceanic crust in the central part of the SCS is characterized by a large positive Bouguer gravity anomaly (220–330 mgal) as well as large free-air and residual isostatic anomalies. There are, however, important differences among spreading segments. For example, in terms of free-air gravity anomaly, the southwest section of mid-ocean has an approximately 50 km wide belt of gravity low superimposed on a broad high of 45 mgal running NW–SE, whereas there are no similar features in other spreading segments. There are indications that gravity anomalies may represent lateral variation in upper crustal density structure. For instance, free air and isostatic anomalies show large positive anomalies in the east of the Namconson basin, which coincide with areas of dense volcanic material known from seismic surveys. The Red River Fault system are clearly identified in the satellite gravity anomalies, including three major faults, Songchay Fault in the southwest, Songlo Fault in the Northeast and Central Fault in the center of the basin. They are elongated in NW–SE direction between 20±30'N and 17°N and reach to Vietnam Scarp Fault around 16°30'N. It is also defined that the crustal density in the south side of the Central Basin is denser than that in the north side of the Central Basin.     

The seismic b-value and its correlation with Bouguer gravity anomaly over the Shillong Plateau area: Tectonic implications

Clues to the understanding of intra- and inter-plate variations in strength or stress state of the crust can be achieved through different lines of evidence and their mutual relationships. Among these parameters Bouguer gravity anomalies and seismic b-values have been widely accepted over several decades for evaluating the crustal character and stress regime. The present study attempts a multivariate analysis for the Shillong Plateau using the Bouguer gravity anomaly and the earthquake database, and establishes a causal relationship between these parameters. Four seismic zones (Zones I–IV), with widely varying b-values, are delineated and an excellent correlation between the seismic b-value and the Bouguer gravity anomaly has been established for the plateau. Low b-values characterize the southwestern part (Zone IV) and a zone (Zone III) of intermediate b-values separates the eastern and western parts of the plateau (Zones I and II) which have high b-values. Positive Bouguer anomaly values as high as +40 mgal, a steep gradient in the Bouguer anomaly map and low b-values in the southwestern part of the plateau are interpreted as indicating a thinner crustal root, uplifted Moho and higher concentration of stress. In comparison, the negative Bouguer anomaly values, flat regional gradient in the Bouguer anomaly map and intermediate to high b-values in the northern part of the plateau are consistent with a comparatively thicker crustal root and lower concentration of stress, with intermittent dissipation of energy through earthquake shocks. Further, depth wise variation in the b-value for different seismic zones, delineated under this study, allowed an appreciation of intra-plateau variation in crustal thickness from ∼30 km in its southern part to ∼38 km in the northern part. The high b-values associated with the depth, coinciding with lower crust, indicate that the Shillong Plateau is supported by a strong lithosphere.     

Gravity anomalies and subsurface geology in the Westerwald volcanic area, Germany

The area of the present study constitutes an alkaline volcanic province in the eastern sector of the Rhenish massif. A series of gravity measurements were carried out on the volcanic fields of the Westerwald. Three-dimensional modelling and wavelength filtering processing techniques were used to analyze the gravity data. The filtered Bouguer anomaly maps show two major regional gravity features: (a) Increasing Bouguer values towards the northeastern part of the study area could be caused by lateral lithological variations within the upper crust. (2) Local negative Bouguer values in the southwest correlate with magmatic materials of intermediate type. The modelling results indicate that the volcanics of the Westerwald are underlain by two different magmatic complexes at a depth in the range 3.3–10 km with density values of 2680 and 2750 kg/m³. The densities assigned to the local igneous intrusions are in the range of 2314–2948 kg/m³ and at depths between 0.4 and 1.3 km. In the NE a diabase bed was modelled to a maximum depth of approximately 1.6 km using the assigned density of 2800 kg/m³.     

最佳向上延拓高度的估计

提出根据两个相邻高度重力异常向上延拓值相关系数与高度的关系 ,以估计应用向上延拓分离区域及剩余重力异常的最佳向上延拓高度的方法。二维模型计算表明 ,不同高度的观测重力异常向上延拓值和观测面上区域重力异常值的互相关系数与高度的关系曲线 ,存在一个明显的极大值 ;这个极大值对应的高度 ,就是从观测异常中分离出这一区域重力异常所需要的最佳向上延拓高度。两个相邻高度重力异常向上延拓值之间的互相关系数与高度的关系曲线 ,存在一个明显的转折点 ,这个转折点对应的高度 ,就是所求的最佳向上延拓高度。应用本方法处理华南北部地区布格重力异常的结果表明 ,由于引起本区区域重力异常的地质因素 ,除了莫霍面及上地壳底面外 ,还受到本区广泛分布甚至出露的花岗岩的影响 ;所以为了从观测异常中分离这一区域异常所需要的最佳向上延拓高度为 2 0 0km ,小于莫霍面及上地壳底界面的平均深度。为了从观测异常中分离出由莫霍面引起的重力异常所需要的向上延拓高度 ,达到 15 0km。因此 ,应用本方法处理实测重力资料 ,必须首先了解引起区域重力异常的场源情况。     

Application of gravity and magnetic methods for the crustal study and delineating associated ores in the western limb of Hazara Kashmir Syntaxis,Northwest Himalayas,Pakistan

The present geophysical study deals with the ores and crustal demonstration of southeastern Hazara and its adjoining areas of Azad Jammu and Kashmir, Pakistan, on the basis of terrestrial gravity and magnetic data. Tectonically, the study area lies in the Lesser Himalayas as well as to an extent in the sub-Himalaya, more specifically in the western limb of Hazara Kashmir Syntaxis. In this study, 567 gravity and 508 magnetic stations have been measured with CG-5 gravimeter and proton precession magnetometer, respectively. The collected data have been processed by applying standard corrections and then different types of maps were prepared. The ores in the area have been delineated by the qualitative interpretation of residual Bouguer anomaly and reduction to pole total magnetic intensity maps, whereas regional structures are demarcated by the Bouguer anomaly and regional Bouguer anomaly maps. The positive contour closures on the residual Bouguer anomaly map indicate the iron ore and phosphate, whereas negative contour closures are the effects of low-density material which consists of gypsum and soapstone. The pole-reduced total intensity map also shows the negative and positive contour closures almost in the same localities and confirms the residual Bouguer anomaly map. The geological model computed on the basis of Bouguer anomaly demarcated a series of faults between different rock units in the study area. The Kashmir Boundary Thrust cuts the western limb of Hazara Kashmir Syntaxis near the apex in the north of Muzaffarabad and marks the boundary between Murree Formation and carbonates of Abbottabad Formation. The gravity model also suggests that the thickness of the crust increases towards the northeast.     

Crustal structure and nature of emplacement of the 85°E Ridge in the Mahanadi offshore based on constrained potential field modeling: Implications for intraplate plume emplaced volcanism

An integrated interpretation of multi-channel seismic reflection, gravity and magnetic datasets belonging to northern most part of the 85°E Ridge in the Mahanadi offshore is carried out to study the crustal structure and mode of its emplacement. The basement structure map of the ridge reveals that it is 130–150 km wide and is composed of an eastern high which appears as a continuous, broad and smooth topographyand the western high characterized by several steep isolated highs. The seismic velocities reported for the first time over the ridge indicate several sedimentary sequences ranging in velocities between 1.6 and 4.0 km/s above the acoustic basement top. The salient aspects of the sedimentary velocities are; a low velocity layer (2.6–3.2 km/s) within the Cretaceous sequence in the intervening depressions encompassing the flank region, and a regionally widespread higher velocity layer (3.5–3.8 km/s) belonging to the Eocene–Oligocene section overlying the ridge. A layer having a velocity of 4.2–4.7 km/s probably made of volcanoclastic rocks is observed immediately below the acoustic basement. The sediment isopach maps presented here for three major horizons are used to compute the 3-D sediment gravity effect to obtain a crustal Bouguer anomaly map of the region. Detailed analysis of the gravity and magnetic anomaly maps clearly demonstrates the continuity of ridge up to the Mahanadi coast at Chilka Lake. Seismically constrained gravity and magnetic models indicate that the ridge is composed of volcanic material that was emplaced on continental crust in the shelf-slope areas and over the oceanic crust in the deep offshore areas. The modeled crustal structure below the ridge further indicates volcanic emplacement of the ridge on a relatively younger lithosphere. We propose two alternative models for the emplacement of the ridge.     

Regional gravity analysis of the crustal structure of Tunisia

Gravity data were integrated with seismic refraction/reflection data, well data and geological investigations to determine a general crustal structure of Tunisia. The gravity data analysis included the construction of a complete Bouguer gravity anomaly map, residual gravity anomaly maps, horizontal gravity gradient maps and a 2.5-D gravity model. Residual gravity anomaly maps illustrate crustal anomalies associated with various structural domains within Tunisia including the Sahel Block, Saharian Flexure, Erg Oriental Basin, Algerian Anticlinorium, Gafsa Trough, Tunisian Trough, Kasserine Platform and the Tell Mountains. Gravity anomalies associated with these features are interpreted to be caused either by thickening or thinning of Palæozoic and younger sediments or by crustal thinning. Analysis of the residual gravity anomaly and horizontal gravity gradient maps also determined a number of anomalies that may be associated with previously unknown structures. A north-south trending gravity model in general indicated similar subsurface bodies as a coincident seismic model. However, thinner Mesozoic sediments within the Tunisian Trough, thinner Palæozoic sediments in the Gafsa Trough, and a greater offset on the Saharian Flexure were required by the gravity data. Additionally, basement uplifts under the Kasserine Platform and Gafsa Trough, not imaged by seismic data, were required by the gravity data. The gravity model revealed two previously unknown basins north and south of the Algerian Anticlinorium (5 km), while the Erg Oriental Basin is composed of at least two sub-basins, each with a depth of 5 km.     

Crustal density structure in the Spanish Central System derived from gravity data analysis (Central Spain)

Shallow and deep sources generate a gravity low in the central Iberian Peninsula. Long-wavelength shallow sources are two continental sedimentary basins, the Duero and the Tajo Basins, separated by a narrow mountainous chain called the Spanish Central System. To investigate the crustal density structure, a multitaper spectral analysis of gravity data was applied. To minimise biases due to misleading shallow and deep anomaly sources of similar wavelength, first an estimation of gravity anomaly due to Cenozoic sedimentary infill was made. Power spectral analysis indicates two crustal discontinuities at mean depths of 31.1 ± 3.6 and 11.6 ± 0.2 km, respectively. Comparisons with seismic data reveal that the shallow density discontinuity is related to the upper crust lower limit and the deeper source corresponds to the Moho discontinuity. A 3D-depth model for the Moho was obtained by inverse modelling of regional gravity anomalies in the Fourier domain. The Moho depth varies between a mean depth of 31 km and 34 km. Maximum depth is located in a NW–SE trough. Gravity modelling points to lateral density variations in the upper crust. The Central System structure is described as a crustal block uplifted by NE–SW reverse faults. The formation of the system involves displacement along an intracrustal detachment in the middle crust. This detachment would split into several high-angle reverse faults verging both NW and SE. The direction of transport is northwards, the detachment probably being rooted at the Moho.     

Two-dimensional spatial analysis of the seismic b-value and the Bouguer gravity anomaly in the southeastern part of the Zagros Fold-and-Thrust Belt,Iran: Tectonic implications

Environmental managers and protection agencies try to assess the magnitudes of earthquakes in regions of seismic activity. For several decades they have used the seismic b-values and Bouguer anomalies for evaluating the crustal character and stress regimes. We have analyzed geostatistically data on both variables to map their spatial distributions in the southeast of the Zagros of Iran. We found a strong correlation between the distribution of the b-value and the Bouguer gravity anomaly in the region. The large Bouguer gravity anomaly values and small b-values all accord with there being a thinner crustal root and a larger concentration of stress in the center. The small to moderate Bouguer gravity anomaly values and intermediate to large b-values accord with the thicker crustal root and the smaller concentration of stress in the northeast. We conclude the southeast of the Zagros, consists of heterogeneous crust, such that accounts for its varied tectonics.     

The structure of the lithosphere beneath the Eastern rift, East Africa, deduced from gravity studies

A compilation of all published and unpublished gravity data for the Eastern rift between latitudes 1°N and 5°S is presented. The Bouguer anomaly map reveals that the shape of the negative regional anomaly associated with the rift is approximately two-dimensional, striking east of north, of width 350 ± 50 km and amplitude500 ± 100 g.u. relative to the background value of−1300 ± 100 g.u. to the west. The regional anomaly is interpreted in terms of an upward thinning of the lithosphere and replacement by low-density asthenosphere. This model is different from previous interpretations in that major lithospheric thinning is restricted to the region of the Eastern rift affected by the domal uplift and does not extend beneath the Lake Victoria region to the west. The gravity and seismic models are compatible if the anomalous upper mantle (asthenospheric part), beneath the rift, is in a state of partial melt. A consequence of the revised regional anomaly is that it reduces previous amplitude estimates of the axial positive residual anomaly within the rift by at least 50% and generates negative anomalies over the rift shoulders in areas covered by Cenozoic volcanics. These negative anomalies are considered to be caused by the low density of the surface volcanics. Within the rift, elongated negative anomalies of amplitude 100–350 g.u. are associated with sedimentary basins and are attributed to low-density sediments up to 3 km thick. The positive residual anomaly along the axis of the rift can be interpreted in terms of either a dyke injection zone less than 15 km wide or by a dense infill body about 2.5 km thick. The positive anomaly is shown to be confined to the volcanic province of the Eastern rift and has its southern termination in the Magadi—Natron area, just north of where the Kenya rift valley changes to block faulting in N. Tanzania. This termination coincides with a change in the spatial distribution of the seismic and geothermal activity.     

闽南区域地壳稳定性分区及其特征

通过一系列区域地壳稳定性背景资料的分析，以地壳结构、深断裂、活动断裂、第四纪升降速率、大地热流值、布格异常梯度、地壳应变能量、地震最大震级、基本烈度等综合指标，将闽南区域的地壳稳定性，自西而东划分为稳定区、基本稳定区和次稳定区。     

Thickness and geothermal gradient of Neoarchean continental crust: Inference from the southeastern North China Craton

The thickness and geothermal gradient of Archean continental crust are critical factors for understanding the geodynamic processes in Earth's early continental crust. Archean tonalite-trondhjemite-granodiorite (TTG) gneisses provide one of the potential indicators of paleo-crustal thickness and geothermal gradient because crust-derived TTG melts are generally thought to originate from partial melting of mafic rocks at the crustal root. In the Western Shandong Province (WSP) of the North China Craton (NCC), two episodes of Neoarchean TTG magmatism are recognized at ~2.70 Ga and ~2.55 Ga which were sourced from partial melting of juvenile crustal components. The ~2.70 Ga TTG gneisses show highly fractionated rare earth element (REE) patterns (average (La/Yb)_N = 39), whereas the ~2.55 Ga TTG gneisses have relatively less fractionated REE patterns (average (La/Yb)_N = 18). Petrogenetic evaluation suggest that the magmatic precursors of the TTG gneisses of both episodes originated from partial melting of juvenile crustal materials at different crustal depths with residual mineral phases of Grt, Cpx, Amp, Pl and Ilm. Together with the garnet proportion in the residue, the P–T pseudosections of equilibrium mineral assemblages, and the temperature calculated from Titanium-in-zircon thermometer, we estimate the Neoarchean crustal thicknesses as 44–51 km with geothermal gradients of 17 to 20 °C/km for the ~2.70 Ga TTG gneisses whereas the ~2.55 Ga TTG gneisses show lesser crustal thicknesses of 35–43 km with geothermal gradients of 19 to 26 °C/km, with an approximately 10 km difference in crustal thickness. Our estimates on the thicknesses and geothermal gradients of the Neoarchean crust are similar to those (~41 km, ~20 °C/km) of the modern average continental crust, indicating that a modern-style plate tectonic regime may have played an important role in the formation and evolution of the Neoarchean continental crust in the NCC.     

Geodynamics of NW India: Subduction, lithospheric flexure, ridges and seismicity

Spectral analysis of digital data of the Bouguer anomaly map of NW India suggests maximum depth of causative sources as 134 km that represents the regional field and coincides with the upwarped lithosphere — asthenosphere boundary as inferred from seismic tomography. This upwarping of the Indian plate in this section is related to the lithospheric flexure due to its down thrusting along the Himalayan front. The other causative layers are located at depths of 33, 17, and 6 km indicating depth to the sources along the Moho, lower crust and the basement under Ganga foredeep, the former two also appear to be upwarped as crustal bulge with respect to their depths in adjoining sections. The gravity and the geoid anomaly maps of the NW India provide two specific trends, NW-SE and NE-SW oriented highs due to the lithospheric flexure along the NW Himalayan fold belt in the north and the Western fold belt (Kirthar -Sulaiman ranges, Pakistan) and the Aravalli Delhi Fold Belt (ADFB) in the west, respectively. The lithospheric flexures also manifest them self as crustal bulge and shallow basement ridges such as Delhi — Lahore — Sagodha ridge and Jaisalmer — Ganganagar ridge. There are other NE-SW oriented gravity and geoid highs that may be related to thermal events such as plumes that affected this region. The ADFB and its margin faults extend through Ganga basin and intersect the NW Himalayan front in the Nahan salient and the Dehradun reentrant that are more seismogenic. Similarly, the extension of NE-SW oriented gravity highs associated with Jaisalmer — Ganganagar flexure and ridge towards the Himalayan front meets the gravity highs of the Kangra reentrant that is also seismogenic and experienced a 7.8 magnitude earthquake in 1905. Even parts of the lithospheric flexure and related basement ridge of Delhi — Lahore — Sargodha show more seismic activity in its western part and around Delhi as compared to other parts. The geoid highs over the Jaisalmer — Ganganagar ridge passes through Kachchh rift and connects it to plate boundaries towards the SW (Murray ridge) and NW (Kirthar range) that makes the Kachchh as a part of a diffused plate boundary, which, is one of the most seismogenic regions with large scale mafic intrusive that is supported from 3-D seismic tomography. The modeling of regional gravity field along a profile, Ganganagar — Chandigarh extended beyond the Main Central Thrust (MCT) constrained from the various seismic studies across different parts of the Himalaya suggests crustal thickening from 35-36 km under plains up to ～56 km under the MCT for a density of 3.1 g/cm³ and 3.25 g/cm³ of the lower most crust and the upper mantle, respectively. An upwarping of ～3 km in the Moho, crust and basement south of the Himalayan frontal thrusts is noticed due to the lithospheric flexure. High density for the lower most crust indicates partial eclogitization that releases copious fluid that may cause reduction of density in the upper mantle due to sepentinization (3.25 g/cm³). It has also been reported from some other sections of Himalaya. Modeling of the residual gravity and magnetic fields along the same profile suggest gravity highs and lows of NW India to be caused by basement ridges and depressions, respectively. Basement also shows high susceptibility indicating their association with mafic rocks. High density and high magnetization rocks in the basement north of Chandigarh may represent part of the ADFB extending to the Himalayan front primarily in the Nahan salient. The Nahan salient shows a basement uplift of ～ 2 km that appears to have diverted courses of major rivers on either sides of it. The shallow crustal model has also delineated major Himalayan thrusts that merge subsurface into the Main Himalayan Thrust (MHT), which, is a decollment plane.