Heterogeneity in crustal structure across the Southern Granulite Terrain (SGT): Inferences from an analysis of gravity and magnetic fields in the Periyar plateau and adjoining areas

The study region forms the western part of the Madurai block (southern block) and shares several lithological characteristics of the Proterozoic exhumed South Indian Granulite Terrain (SGT). The crustal structure of the area has been derived from gravity data, constrained partly by aeromagnetic data. The Bouguer anomaly map of the region prepared based on detailed gravity observations shows a number of features (i) the Periyar lineament separates two distinctly different gravity fields, one, a high gravity gradient tending to be positive towards the coast in south west and significant gravity lows ranging from − 85 to as low as − 150 mGal in the NE covering a large part of the Periyar plateau (ii) within the broad gravity low, three localised circular anomalies of considerable amplitude occur in the region of Munnar granite. A magnetic low region in the central part coincides with the area of retrogressed charnockites and the major lineaments suggestive of a genetic link and considerable downward extent. The crustal models indicate that the upper layer containing exhumed lower crustal rocks (2.76 gm/cc) is almost homogeneous, most part of the gravity field resulting from variations in intracrustal layers of decharnockitised hornblendic gneisses and granite bodies. Below it, a denser layer (2.85 gm/cc) of unknown composition exists with Moho depth ranging from 36 to 41 km. The structure below the region is compared with that of two other segments of the SGT from which it differs markedly. The Wynad plateau forming the western part of the Northern Block of the SGT is characterised by a heterogeneity due to the presence of contrasting crustal blocks on either side of the Bavali shear zone, possibly a westward extension of the Moyar shear zone and presence of high density material in the mid-to-lower crustal portions. The crust below the Kuppam–Palani transect has a distinctive four-layer structure with a mid-crustal low density layer. The differences in crustal structure are consistent with the different tectonic settings of the three regions discussed in the paper. It is suggested that the crustal structure below the Kuppam–Palani transect corridor is not representative of the SGT as a whole, an aspect of great relevance to intra-continental comparisons and trans-continental reconstructions of continent configurations of the Gondwanaland.     

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.     

Long wavelength gravity anomalies over India: Crustal and lithospheric structures and its flexure

Long wavelength gravity anomalies over India were obtained from terrestrial gravity data through two independent methods: (i) wavelength filtering and (ii) removing crustal effects. The gravity fields due to the lithospheric mantle obtained from two methods were quite comparable. The long wavelength gravity anomalies were interpreted in terms of variations in the depth of the lithosphere–asthenosphere boundary (LAB) and the Moho with appropriate densities, that are constrained from seismic results at certain points. Modeling of the long wavelength gravity anomaly along a N–S profile (77°E) suggest that the thickness of the lithosphere for a density contrast of 0.05 g/cm³ with the asthenosphere is maximum of ∼190 km along the Himalayan front that reduces to ∼155 km under the southern part of the Ganga and the Vindhyan basins increasing to ∼175 km south of the Satpura Mobile belt, reducing to ∼155–140 km under the Eastern Dharwar craton (EDC) and from there consistently decreasing south wards to ∼120 km under the southernmost part of India, known as Southern Granulite Terrain (SGT).The crustal model clearly shows three distinct terrains of different bulk densities, and thicknesses, north of the SMB under the Ganga and the Vindhyan basins, and south of it the Eastern Dharwar Craton (EDC) and the Southern Granulite Terrain (SGT) of bulk densities 2.87, 2.90 and 2.96 g/cm³, respectively. It is confirmed from the exposed rock types as the SGT is composed of high bulk density lower crustal rocks and mafic/ultramafic intrusives while the EDC represent typical granite/gneisses rocks and the basement under the Vindhyan and Ganga basins towards the north are composed of Bundelkhand granite massif of the lower density. The crustal thickness along this profile varies from ∼37–38 km under the EDC, increasing to ∼40–45 km under the SGT and ∼40–42 km under the northern part of the Ganga basin with a bulge up to ∼36 km under its southern part. Reduced lithospheric and crustal thicknesses under the Vindhyan and the Ganga basins are attributed to the lithospheric flexure of the Indian plate due to Himalaya. Crustal bulge due to lithospheric flexure is well reflected in isostatic Moho based on flexural model of average effective elastic thickness of ∼40 km. Lithospheric flexure causes high heat flow that is aided by large crustal scale fault system of mobile belts and their extensions northwards in this section, which may be responsible for lower crustal bulk density in the northern part. A low density and high thermal regime in north India north of the SMB compared to south India, however does not conform to the high S-wave velocity in the northern part and thus it is attributed to changes in composition between the northern and the southern parts indicating a reworked lithosphere. Some of the long wavelength gravity anomalies along the east and the west coasts of India are attributed to the intrusives that caused the breakup of India from Antarctica, and Africa, Madagascar and Seychelles along the east and the west coasts of India, respectively.     

Crustal structure across southern Mexico inferred from gravity data

We present a gravity model of the crustal structure in southern Mexico based on interpretation of a detailed marine gravity profile perpendicularly across the Middle America Trench offshore from Acapulco, and a regional gravity transect extending into continental Mexico across the Sierra Madre del Sur, the central sector of the Trans-Mexican Volcanic Belt, the Sierra Madre Oriental, the Coastal Plain, and into the Gulf of Mexico. The elastic thickness of the Cocos lithospheric plate was found to be 30 km. In agreement with a previous seismic refraction study, no major differences in crustal structure were observed on both sides of the O’Gorman Fracture Zone. The gravity high seaward of the trench is interpreted as due to the incipient flexure and crustal thinning. The gravity low at the axis of the trench is explained by the increase in water depth and the existence of low-density accreted or continental-derived sediments (2.25 and 2.40 g/cm³). A gravity high of 50 mGal extending about 100 km landward is interpreted as caused by local shoaling of the Moho. The crust attains a thickness of 42 km under the Trans-Mexican Volcanic Belt but thins beneath the Coastal Plain and the continental slope of the Gulf of Mexico. Gravity highs around the Sierra de Tamaulipas are interpreted in terms of relief of the lower–upper crustal interface, implying a shallow basement.     

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.     

Crustal character and thickness over the Dangerous Grounds and beneath the Northwest Borneo Trough

New deep reflection seismic, bathymetry, gravity and magnetic data have been acquired in a marine geophysical survey of the southern South China Sea, including the Dangerous Grounds, Northwest Borneo Trough and the Central Luconia Platform. The seismic and bathymetry data map the topography of shallow density interfaces, allowing the application of gravity modeling to delineate the thickness and composition of the deeper crustal layers. Many of the strongest gravity anomalies across the area are accounted for by the basement topography mapped in the seismic data, with substantial basement relief associated with major rift development. The total crustal thickness is however quite constant, with variations only between 25 and 30 km across the Central Luconia Platform and Dangerous Grounds. The Northwest Borneo Trough is underlain by thinned crust (25–20 km total crustal thickness) consistent with the substantial water depths. There is no evidence of any crustal suture associated with the trough, nor any evidence of relict oceanic crust beneath the trough. The crustal thinning also does not extend along the complete length of the trough, with crustal thicknesses of 25 km and more modeled on the most easterly lines to cross the trough. Modeled magnetic field variations are also consistent with the study area being underlain by continental crust, with the magnetic field variations well explained by irregular magnetisations consistent with inhomogeneous continental crust, terminating at the basement unconformity as mapped from the seismic data.     

Gravity field and isostatic state of Ethiopia and adjacent areas

Over 35,000 onshore and offshore gravity stations have been compiled in order to test isostatic models against geologic structures over a part of the Afro–Arabian shield. The area of Ethiopia covers an important part of this system because it contains the major section of the ≈5000 km Afro–Arabian rift and includes the transition between the Arabo-Nubian-Shield (ANS) and the Mozambique Belt (MB).Isostatic residual anomalies have been calculated using both Airy and Vening-Meinesz (flexural rigidity D = 10²² Nm) models. The isostatic residual anomalies outline the major Precambrian belts, the Cenozoic rifts and associated major structures. Positive residual anomalies associated with the main Ethiopian Rift (MER) and Kenyan rift systems could be the expressions of an axial intrusive body and swarms of local faults and fractures. The residual anomalies indicate relative stability in the MER and increased tectonic activity in the areas of the Red Sea, Gulf of Aden and Afar. Near-zero isostatic residuals flank the MER and Kenya rifts and are found within the Danakil Alps and some plateau regions.The small mean isostatic residual anomaly (about 8 mGal) and the isostatic analysis show a slight positive bias indicating under compensation. The undercompensation may imply that there are upper crustal features that are not compensated regionally (probably supported by the rigidity of the lithosphere) and isostatic disequilibrium in the region. Therefore, the high topography of Ethiopia and East African plateau is partly compensated by thicker crust (broad negative isostatic regional anomaly) and partly by dynamic forces.The results of the qualitative interpretation form the basis of continuing three-dimensional gravity modelling and quantitative analysis that also integrates data from eastern Sudan.     

Lithospheric structure of the North China Craton: Integrated gravity,geoid and topography data

The lithospheric structure of ancient cratons provides important constraints on models relating to tectonic evolution and mantle dynamics. Here we present the 3D lithospheric structure of the North China Craton (NCC) from a joint inversion of gravity, geoid and topography data. The NCC records a prolonged history of Archean and Paleoproterozoic accretion of crustal blocks through subduction and collision building the cratonic architecture, which was subsequently differentially destroyed during Mesozoic through extensive magmatism. The thermal structure obtained in our study is considered to define the lithosphere-asthenosphere boundary (LAB) of the NCC, and reflects the density variations within the mantle lithosphere. Employing the Moho depths from deep seismic sounding profiles for the inversion, and based on repeated computations using different parameters, we estimate the Moho depth, LAB depth and average crustal density of the craton. The Moho depth varies from 28 to 50 km and the LAB depth varies from 105 to 205 km. The LAB and Moho show concordant thinning from West to East of the NCC. The average crustal density is 2870 kg m^− 3 in the western part of the NCC, higher than that in the eastern part (2750 kg m^− 3). The results of joint inversion in our study yielded LAB depth and lithospheric thinning features similar to those estimated from thermal and seismic studies, although our results show different depth and variations in the thickness. The lithosphere gently thins from 145 to 105 km in the eastern NCC, where as the thinning is much less pronounced in the western NCC with average depth of about 175 km. The joint inversion results in this study provide another perspective on the lithospheric structure from the density properties and corresponding geophysical responses in an ancient craton.     

Variscan Sn–W–Mo metallogeny in the gravity picture of the Krušné hory/Erzgebirge granite batholith (Central Europe)

Late Variscan wolframite (± molybdenite) and cassiterite–wolframite greisen, skarn and vein deposits occur in a close spatial association with the granites of the Krušné hory/Erzgebirge batholith (KHEB) in Central Europe. We examined the distribution of the deposits in relation to the gravity field affected by Late Variscan granites using the data from previous gravity and metallogenic studies. Late Variscan granites are differentiated into earlier biotite monzogranites (low-F granites) and later biotite or lithium mica syenogranites (high-F granites) in accordance with the previous classifications. All the outcrops of granites in the KHEB region and their hidden continuation are confined to the Bouguer anomaly contour of − 20 mGal. The Sn–W–Mo (rare metal) deposits and occurrences are within the gravity contour of − 30 mGal with the exception of the Grossschirma stratiform tin deposit in the Freiberg polymetallic ore district. We constructed a geological model based on the gravity data along two profiles across the KHEB showing the position of some rare metal deposits and of outcropping and hidden granite bodies. The models show that the overlapping of earlier and later granites is in the areas of the most intense regional gravity minima. These coincide with the Eastern Volcano-Plutonic Complex (Altenberg minimum), which encloses large volumes of felsic extrusives, microgranite dikes and granites, and the Western Plutonic Complex (Eibenstock minimum), with small volumes of felsic dikes and predominance of earlier and later granites, with no extrusives preserved. There is no distinct relationship between the masses of Late Variscan granites and the distribution and the sizes of associated W ± Mo and Sn–W deposits. We prefer the idea that rare metal mineralization was formed by hydrothermal fluids derived from outside of presently outcropping granites. It originated in two cycles: one connected with the formation of earlier granites producing W ± Mo associations and the other one associated with later granites connected with Sn–W mineralization. Mineralizing fluids were probably generated by mantle–crustal interaction in the crust near the mantle–crust boundary as also indicated by lamprophyric intrusions coeval with the Late Variscan granitic magmatism.     

Crustal shear-wave velocity structure beneath northeast India from teleseismic receiver function analysis

We investigated the seismic shear-wave velocity structure of the crust beneath nine broadband seismological stations of the Shillong–Mikir plateau and its adjoining region using teleseismic P-wave receiver function analysis. The inverted shear wave velocity models show ∼34–38 km thick crust beneath the Shillong Plateau which increases to ∼37–38 km beneath the Brahmaputra valley and ∼46–48 km beneath the Himalayan foredeep region. The gradual increase of crustal thickness from the Shillong Plateau to Himalayan foredeep region is consistent with the underthrusting of Indian Plate beyond the surface collision boundary. A strong azimuthal variation is observed beneath SHL station. The modeling of receiver functions of teleseismic earthquakes arriving the SHL station from NE backazimuth (BAZ) shows a high velocity zone within depth range 2–8 km along with a low velocity zone within ∼8–13 km. In contrast, inversion of receiver functions from SE BAZ shows high velocity zone in the upper crust within depth range ∼10–18 km and low velocity zone within ∼18–36 km. The critical examination of ray piercing points at the depth of Moho shows that the rays from SE BAZ pierce mostly the southeast part of the plateau near Dauki fault zone. This observation suggests the effect of underthrusting Bengal sediments and the underlying oceanic crust in the south of the plateau facilitated by the EW-NE striking Dauki fault dipping 30⁰ toward northwest.     

Density structure of the cratonic mantle in southern Africa: 1. Implications for dynamic topography

The origin of high topography in southern Africa is enigmatic. By comparing topography in different cratons, we demonstrate that in southern Africa both the Archean and Proterozoic blocks have surface elevation 500–700 m higher than in any other craton worldwide, except for the Tanzanian Craton. An unusually high topography may be caused by a low density (high depletion) of the cratonic lithospheric mantle and/or by the dynamic support of the mantle with origin below the depth of isostatic compensation (assumed here to be at the lithosphere base). We use free-board constraints to examine the relative contributions of the both factors to surface topography in the cratons of southern Africa. Our analysis takes advantage of the SASE seismic experiment which provided high resolution regional models of the crustal thickness.We calculate the model of density structure of the lithospheric mantle in southern Africa and show that it has an overall agreement with xenolith-based data for lithospheric terranes of different ages. Density of lithospheric mantle has significant short-wavelength variations in all tectonic blocks of southern Africa and has typical SPT values of ca. 3.37–3.41 g/cm³ in the Cape Fold and Namaqua–Natal fold belts, ca. 3.34–3.35 g/cm³ in the Proterozoic Okwa block and the Bushveld Intrusion Complex, ca. 3.34–3.37 g/cm³ in the Limpopo Belt, and ca. 3.32–3.33 g/cm³ in the Kaapvaal and southern Zimbabwe cratons.The results indicate that 0.5–1.0 km of surface topography, with the most likely value of ca. 0.5 km, cannot be explained by the lithosphere structure within the petrologically permitted range of mantle densities and requires the dynamic (or static) contribution from the sublithospheric mantle. Given a low amplitude of regional free air gravity anomalies (ca. + 20 mGal on average), we propose that mantle residual (dynamic) topography may be associated with the low-density region below the depth of isostatic compensation. A possible candidate is the low velocity layer between the lithospheric base and the mantle transition zone, where a temperature anomaly of 100–200 °C in a ca. 100–150 km thick layer may explain the observed reduction in Vs velocity and may produce ca. 0.5–1.0 km to the regional topographic uplift.     

40Ar/39Ar age pattern associated with differential uplift along the Eastern Highlands shear zone,Cape Breton Island,Canadian Appalachians

The Eastern Highlands shear zone in Cape Breton Island is a crustal scale thrust. It is characterized by an amphibolite-facies deformation zone ∼5 km wide formed deep in the crust that is overprinted by a greenschist-facies mylonite zone ∼1 km wide that formed at a more shallow level. Hornblende ⁴⁰Ar/³⁹Ar plateau ages on the hanging wall decrease towards the centre of the shear zone. In the older zone (over 7.8 km from the centre), the ages are between ∼565 and ∼545 Ma; in the younger zone (within 4.5 km of the centre), they are between ∼425 and ∼415 Ma; and in the transitional zone in between, they decrease abruptly from ∼545 to ∼425 Ma. Pressures of crystallization of plutons in the hanging wall, based on the Al-in-hornblende barometer and corresponding to depth of emplacement, increase towards the centre of the shear zone and indicate a differential uplift of up to ∼28 km associated with movement along the shear zone. The age pattern is interpreted to have resulted from the differential uplift. The pressure data show that rocks exposed in the younger zone were buried deep in the crust and did not cool through the hornblende Ar blocking temperature (∼500°C) until differential uplift occurred. The ⁴⁰Ar/³⁹Ar ages in the zone (∼425–415 Ma) thus date shear zone movement or the last stage of it. In contrast, rocks in the older zone were more shallowly buried before differential uplift and cooled through the blocking temperature soon after the emplacement of ∼565–555 Ma plutons in the area, long before shear zone movement. The transitional zone corresponds to the Ar partial retention zone before differential uplift. The ⁴⁰Ar/³⁹Ar age pattern thus reflects a Neoproterozoic to Silurian cooling profile that was exposed as a result of differential uplift related to movement along the shear zone. A similar K–Ar age pattern has been reported for the Alpine fault in New Zealand. It is suggested that such isotopic age patterns can be used to help constrain the ages, kinematics, displacements and depth of penetration of shear zones.     

Proterozoic orogens in southern Peninsular India: Contiguities and complexities

The Precambrian terranes of southern Peninsular India have been central to discussions on the history of formation and breakup of supercontinents. Of particular interest are the Proterozoic high grade metamorphic orogens at the southern and eastern margins of the Indian shield, skirting the 3.4 Ga Dharwar craton which not only preserve important records of lower crustal processes and lithospheric geodynamics, but also carry imprints of the tectonic framework related to the assembly of the major Neoproterozoic supercontinents – Rodinia and Gondwana. These Proterozoic orogens are described as Southern Granulite Terrane (SGT) in the southern tip and the Eastern Ghats Mobile Belt (EGMB) in the eastern domains of the peninsula. The contiguity of these orogens is broken for a distance of ∼400 km and disappears in the Bay of Bengal. These orogens expose windows of middle to lower crust with well-preserved rock records displaying multiple tectonothermal events and multiphase exhumation paths.Recent studies in these orogens have led to the recognition of discrete crustal blocks or terranes separated by major shear zone systems, some of which represent collisional sutures. The SGT and EGMB carry several important features such as fold-thrust tectonics, regional granulite facies metamorphism of up to ultrahigh-temperature conditions in some cases, multiple P–T paths, development of lithospheric shear zones, emplacement of ophiolites, presence of alkaline and anorthositic complexes, development of crustal-scale “flower structures”, transpressional strains, and reactivation tectonics. A heterogeneous distribution of different metamorphic and magmatic assemblages with distinct spatial and temporal strain variations in shaping the fabric elements in different blocks is identified. Both EGMB and SGT share a common transpressional deformation history during the latest Neoproterozoic characterized by the steepening of the initial low angle crustal scale structures leading to a subvertical grain conducive to reactivation tectonics. Our synthesis of the spatial distribution, geometry, kinematics and the transpressional strain of the shear zone systems provides insights into the tectono-metamorphic history of the Proterozoic orogens of southern India and their contiguity and complexities. Recent understanding of subduction, accretion and collisional history along these zones together with a long lived transpressional tectonic regime imply that these orogens witnessed identical tectonic regimes at different times in Earth history, although the major and common structural architecture was built during the final assembly of the Gondwana supercontinent.     

Anomalous crustal and lithospheric mantle structure of southern part of the Vindhyan Basin and its geodynamic implications

Tectonically active Vindhyan intracratonic basin situated in central India, forms one of the largest Proterozoic sedimentary basins of the world. Possibility of hydrocarbon occurrences in thick sediments of the southern part of this basin, has led to surge in geological and geophysical investigations by various agencies. An attempt to synthesize such multiparametric data in an integrated manner, has provided a new understanding to the prevailing crustal configuration, thermal regime and nature of its geodynamic evolution. Apparently, this region has been subjected to sustained uplift, erosion and magmatism followed by crustal extension, rifting and subsidence due to episodic thermal interaction of the crust with the hot underlying mantle. Almost 5–6 km thick sedimentation took place in the deep faulted Jabera Basin, either directly over the Bijawar/Mahakoshal group of mafic rocks or high velocity-high density exhumed middle part of the crust. Detailed gravity observations indicate further extension of the basin probably beyond NSL rift in the south. A high heat flow of about 78 mW/m² has also been estimated for this basin, which is characterized by extremely high Moho temperatures (exceeding 1000 °C) and mantle heat flow (56 mW/m²) besides a very thin lithospheric lid of only about 50 km. Many areas of this terrain are thickly underplated by infused magmas and from some segments, granitic–gneissic upper crust has either been completely eroded or now only a thin veneer of such rocks exists due to sustained exhumation of deep seated rocks. A 5–8 km thick retrogressed metasomatized zone, with significantly reduced velocities, has also been identified around mid to lower crustal transition.     

Lithosphere,crust and basement ridges across Ganga and Indus basins and seismicity along the Himalayan front,India and Western Fold Belt,Pakistan

Spectral analysis of the digital data of the Bouguer anomaly of North India including Ganga basin suggest a four layer model with approximate depths of 140, 38, 16 and 7 km. They apparently represent lithosphere–asthenosphere boundary (LAB), Moho, lower crust, and maximum depth to the basement in foredeeps, respectively. The Airy’s root model of Moho from the topographic data and modeling of Bouguer anomaly constrained from the available seismic information suggest changes in the lithospheric and crustal thicknesses from ∼126–134 and ∼32–35 km under the Central Ganga basin to ∼132 and ∼38 km towards the south and 163 and ∼40 km towards the north, respectively. It has clearly brought out the lithospheric flexure and related crustal bulge under the Ganga basin due to the Himalaya. Airy’s root model and modeling along a profile (SE–NW) across the Indus basin and the Western Fold Belt (WFB), (Sibi Syntaxis, Pakistan) also suggest similar crustal bulge related to lithospheric flexure due to the WFB with crustal thickness of 33 km in the central part and 38 and 56 km towards the SE and the NW, respectively. It has also shown the high density lower crust and Bela ophiolite along the Chamman fault. The two flexures interact along the Western Syntaxis and Hazara seismic zone where several large/great earthquakes including 2005 Kashmir earthquake was reported.The residual Bouguer anomaly maps of the Indus and the Ganga basins have delineated several basement ridges whose interaction with the Himalaya and the WFB, respectively have caused seismic activity including some large/great earthquakes. Some significant ridges across the Indus basin are (i) Delhi–Lahore–Sargodha, (ii) Jaisalmer–Sibi Syntaxis which is highly seismogenic. and (iii) Kachchh–Karachi arc–Kirthar thrust leading to Sibi Syntaxis. Most of the basement ridges of the Ganga basin are oriented NE–SW that are as follows (i) Jaisalmer–Ganganagar and Jodhpur–Chandigarh ridges across the Ganga basin intersect Himalaya in the Kangra reentrant where the great Kangra earthquake of 1905 was located. (ii) The Aravalli Delhi Mobile Belt (ADMB) and its margin faults extend to the Western Himalayan front via Delhi where it interacts with the Delhi–Lahore ridge and further north with the Himalayan front causing seismic activity. (iii) The Shahjahanpur and Faizabad ridges strike the Himalayan front in Central Nepal that do not show any enhanced seismicity which may be due to their being parts of the Bundelkhand craton as simple basement highs. (iv) The west and the east Patna faults are parts of transcontinental lineaments, such as Narmada–Son lineament. (v) The Munghyr–Saharsa ridge is fault controlled and interacts with the Himalayan front in the Eastern Nepal where Bihar–Nepal earthquakes of 1934 has been reported. Some of these faults/lineaments of the Indian continent find reflection in seismogenic lineaments of Himalaya like Everest, Arun, Kanchenjunga lineaments. A set of NW–SE oriented gravity highs along the Himalayan front and the Ganga and the Indus basins represents the folding of the basement due to compression as anticlines caused by collision of the Indian and the Asian plates. This study has also delineated several depressions like Saharanpur, Patna, and Purnia depressions.     

Distinct variations of crustal shear wave velocity structure and radial anisotropy beneath the North China Craton and tectonic implications

This study presents the crustal shear wave velocity structure and radial anisotropy along two linear seismic arrays across the North China Craton (NCC) from ambient noise tomography. About a half to one year long ambient noise data from 87 stations were used for obtaining the inter-station surface wave empirical Green's functions (EGFs) from cross-correlation. Rayleigh and Love dispersion curves within the period band 5–30 s were measured from the EGFs of the vertical and transverse components, respectively. These dispersion data were then used to determine the crustal shear wave velocity structure (V_SV and V_SH) and radial anisotropy (2(V_SH ? V_SV) / (V_SH + V_SV)) from point-wise linear inversion with constraints from receiver function analysis. Our results reveal substantial structural variations among different parts of the NCC. The Bohai Bay Basin in the eastern NCC is underlain by a thin crust (~ 30 km) with relatively low velocities (particularly V_SV) and large positive radial anisotropy in the middle to lower crust. Such a crustal structure is no longer of a cratonic type and may have resulted from the widespread tectonic extension and intensive magmatism in this region since late Mesozoic. Beneath the Ordos Basin in the western NCC, the crust is relatively thicker (≥ 40 km) and well stratified, and presents a large-scale low velocity zone in the middle to lower crust and overall weak radial anisotropy except for a localized lower crust anomaly. The overall structural features of this region resemble those of typical Precambrian shields, in agreement with the long-term stability of the region. The crustal structure under the Trans North China Orogen (TNCO, central NCC) is more complicated and characterized by smaller scale velocity variations, strong positive radial anisotropy in the middle crust and rapid change to weak-to-negative anisotropy in the lower crust. These features may reflect complex deformations and crust–mantle interactions, probably associated with tectonic extension and magmatic underplating during the Mesozoic to Cenozoic evolution of the region. Our structural images in combination with previous seismic, geological and geochemical observations suggest that the Phanerozoic lithospheric reactivation and destruction processes may have affected the crust (especially the middle and lower crust) of the eastern NCC, and the effect probably extended to the TNCO, but may have minor influence on the crust of the western part of the craton.     

Cumberland batholith,Trans-Hudson Orogen,Canada: Petrogenesis and implications for Paleoproterozoic crustal and orogenic processes

Large volume, plutonic belts, such as the ～ 221,000 km², ca. 1.865–1.845 Ga Cumberland batholith (CB) of the Trans-Hudson Orogen in Canada, are major components of Paleoproterozoic orogenic belts. In many cases, they have been interpreted as continental arc batholiths. The petrogenesis and tectonic context of the CB and implications for crustal growth and recycling are interpreted herein based on a 900 km geochemical-isotopic (Nd–O) transect across it and into granitoid plutons within bounding Archean cratons in central and southern Baffin Island.The mainly granulite grade CB, emplaced over an age span of between 14 and 24 Ma, consists mainly of high-K to shoshonitic monzogranite and granodiorite, but also includes low- and medium-K granitoid rocks. Metaluminous to slightly peraluminous compositions and δ¹⁸O (VSMOW) values (+ 6 to + 10‰) indicate derivation from infracrustal (I-type) sources. ε_Nd 1.85 Ga signatures (? 12 to ? 2) of both mafic and felsic units suggest a dominance of evolved sources. Isotopic signatures in the interior of the CB (? 2 to ? 7) are more radiogenic than those within Archean domains in central (? 8 to ? 15) and southern (? 5 to ? 19) Baffin Island. The isotopic transect is interpreted as ‘imaging’ an accreted microcontinental block (Meta Incognita) and bounding Archean cratons. The CB includes granites of arc, within-plate (A-type) and post-collisional affinity and volumetrically minor mafic rocks with both arc and non-arc features. (La/Yb)_CN and Sr/Y values range from < 1 to 225 and < 1 to 611, respectively. In these respects, some CB granitoid rocks resemble Paleozoic adakitic granites, interpreted as partial melts of greatly thickened crust within post-collisional settings, such as Tibet. Thus, the CB likely encompasses various non-consanguineous magmatic suites generated at deep- to mid-crustal depths. Although CB granitoid rocks undoubtedly had important crustal sources, it is hard to assess the relative contribution of mantle-derived magmas.The CB is best interpreted as a post-accretion batholith resulting from large-scale lithospheric mantle delamination followed by the upwelling of hot asthenospheric mantle leading to voluminous crustal partial melting. Contributors to crustal instability which may have facilitated such delamination included: (a) a collage of recently assembled small cratons underlain by hot, weak lithosphere with mantle-depth structural breaks within this segment of the Trans-Hudson Orogen; (b) the gabbro-eclogite phase transformation, and (c) a greatly thickened crustal section (> 60 km), as evidenced by adakitic granites.     

Analysis of airborne magnetic and gravity anomalies of peninsular shield, India integrated with seismic and magnetotelluric results and gravity anomalies of Madagascar, Sri Lanka and East Antarctica

Modelling of gravity and airborne magnetic data integrated with seismic studies suggest that the linear gravity and magnetic anomalies associated with Moyar Bhavani Shear Zone (MBSZ) and Palghat Cauvery Shear Zone (PCSZ) are caused by high density and high susceptibility rocks in upper crust which may represent mafic lower crustal rocks. This along with thick crust (44–45 km) under the Southern Granulite Terrain (SGT) indicates collision of Dharwar craton towards north and SGT towards south with N–S directed compression during 2.6–2.5 Ga. This collision may be related to contemporary collision northwards between Eastern Madagascar–Western Dharwar Craton (WDC) and Eastern Dharwar Craton (EDC). Arcuate shaped N and S-verging thrusts, MBSZ-Mettur Shear and PCSZ-Gangavalli Shear, respectively across Cauvery Shear zone system (CSZ) in SGT also suggest that the WDC, EDC and SGT might have collided almost simultaneously during 2.6–2.5 Ga due to NW–SE directed compressional forces with CSZ as central core complex in plate tectonics paradigm preserving rocks of oceanic affinity. Gravity anomalies of schist belts of WDC suggest marginal and intra arc basin setting.The gravity highs of EGFB along east coast of India and regional gravity low over East Antarctica are attributed to thrusted high-density lower crustal/upper mantle rocks at a depth of 5–6 km along W-verging thrust, which is supported by high seismic velocity and crustal thickening, respectively. It may represent a collision zone at about 1.0 Ga between India and East Antarctica. Paired gravity anomalies in the central part of Sri Lanka related to high density intrusives under western margin of Highland Complex and crustal thickening (40 km) along eastern margin of Highland Complex with several arc type magmatic rocks of about 1.0 Ga in Vijayan Complex towards the east may represent collision between them with W-verging thrust as in case of EGFB. The gravity high of Sri Lanka in the central part falls in line with that of EGFB, in case it is fitted in Gulf of Mannar and may represent the extension of this orogeny in Sri Lanka.     

Geophysical constraints on geodynamic processes at convergent margins: A global perspective

Convergent margins, being the boundaries between colliding lithospheric plates, form the most disastrous areas in the world due to intensive, strong seismicity and volcanism. We review global geophysical data in order to illustrate the effects of the plate tectonic processes at convergent margins on the crustal and upper mantle structure, seismicity, and geometry of subducting slab. We present global maps of free-air and Bouguer gravity anomalies, heat flow, seismicity, seismic Vs anomalies in the upper mantle, and plate convergence rate, as well as 20 profiles across different convergent margins. A global analysis of these data for three types of convergent margins, formed by ocean–ocean, ocean–continent, and continent–continent collisions, allows us to recognize the following patterns. (1) Plate convergence rate depends on the type of convergent margins and it is significantly larger when, at least, one of the plates is oceanic. However, the oldest oceanic plate in the Pacific ocean has the smallest convergence rate. (2) The presence of an oceanic plate is, in general, required for generation of high-magnitude (M > 8.0) earthquakes and for generating intermediate and deep seismicity along the convergent margins. When oceanic slabs subduct beneath a continent, a gap in the seismogenic zone exists at depths between ca. 250 km and 500 km. Given that the seismogenic zone terminates at ca. 200 km depth in case of continent–continent collision, we propose oceanic origin of subducting slabs beneath the Zagros, the Pamir, and the Vrancea zone. (3) Dip angle of the subducting slab in continent–ocean collision does not correlate neither with the age of subducting oceanic slab, nor with the convergence rate. For ocean–ocean subduction, clear trends are recognized: steeply dipping slabs are characteristic of young subducting plates and of oceanic plates with high convergence rate, with slab rotation towards a near-vertical dip angle at depths below ca. 500 km at very high convergence rate. (4) Local isostasy is not satisfied at the convergent margins as evidenced by strong free air gravity anomalies of positive and negative signs. However, near-isostatic equilibrium may exist in broad zones of distributed deformation such as Tibet. (5) No systematic patterns are recognized in heat flow data due to strong heterogeneity of measured values which are strongly affected by hydrothermal circulation, magmatic activity, crustal faulting, horizontal heat transfer, and also due to low number of heat flow measurements across many margins. (6) Low upper mantle Vs seismic velocities beneath the convergent margins are restricted to the upper 150 km and may be related to mantle wedge melting which is confined to shallow mantle levels.     

Seismic structure of the crust and uppermost mantle of South America and surrounding oceanic basins

We present a new set of contour maps of the seismic structure of South America and the surrounding ocean basins. These maps include new data, helping to constrain crustal thickness, whole-crustal average P-wave and S-wave velocity, and the seismic velocity of the uppermost mantle (P_n and S_n). We find that: (1) The weighted average thickness of the crust under South America is 38.17 km (standard deviation, s.d. ±8.7 km), which is ∼1 km thinner than the global average of 39.2 km (s.d. ±8.5 km) for continental crust. (2) Histograms of whole-crustal P-wave velocities for the South American crust are bi-modal, with the lower peak occurring for crust that appears to be missing a high-velocity (6.9–7.3 km/s) lower crustal layer. (3) The average P-wave velocity of the crystalline crust (P_cc) is 6.47 km/s (s.d. ±0.25 km/s). This is essentially identical to the global average of 6.45 km/s. (4) The average P_n velocity beneath South America is 8.00 km/s (s.d. ±0.23 km/s), slightly lower than the global average of 8.07 km/s. (5) A region across northern Chile and northeast Argentina has anomalously low P- and S-wave velocities in the crust. Geographically, this corresponds to the shallowly-subducted portion of the Nazca plate (the Pampean flat slab first described by Isacks et al., 1968), which is also a region of crustal extension. (6) The thick crust of the Brazilian craton appears to extend into Venezuela and Colombia. (7) The crust in the Amazon basin and along the western edge of the Brazilian craton may be thinned by extension. (8) The average crustal P-wave velocity under the eastern Pacific seafloor is higher than under the western Atlantic seafloor, most likely due to the thicker sediment layer on the older Atlantic seafloor.