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.     

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.     

Structure and evolution of the Cauvery Shear Zone system, Southern Granulite Terrain, India: Evidence from deep seismic and other geophysical studies

The Southern Granulite Terrain with exposed Archean lower crustal rocks is studied using various geophysical tools. The crustal structure derived from seismic reflection and refraction/wide-angle reflection studies is used to understand the tectonic evolution of the region. Deep seismic reflection section along the Kolattur–Palani segment shows an oppositely dipping reflection fabric near the Moyar–Bhavani shear zone, which is interpreted as a signature of collision between the Dharwar craton and another crustal block in the south. The thickened crust due to collision was delaminated during the orogenic collapse and modified the central part, covering the Cauvery Shear Zone system, located between the Moyar–Bhavani and Karur–Oddanchatram shear zones. The delaminated lower crust is altered by magmatic underplating as evidenced by the high velocity layer just above the Moho. The velocity model of the region indicates crustal thickening at the boundary of the Dharwar craton and Moyar–Bhavani shear zone and thinning further south. Back-scattered seismic wave field with negative moveout and the Moho-offset indicate the spatial location and strike-slip nature of the shear zones. Present study suggests that the late Archean collision and suturing of the Dharwar craton with the southern crustal block at the Moyar–Bhavani shear zone may be responsible for the evolution of late Archean granulites. Late Neoproterozoic rifting is observed along the paleo-fault zones. The seismic studies constrained by gravity, magnetic and magnetotelluric data suggest that the Moyar–Bhavani and Karur–Oddanchatram shear zones of the Cauvery Shear Zone system mark terrane boundaries/suture zones.     

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.     

Regional geoelectric structure beneath Deccan Volcanic Province of the Indian subcontinent using magnetotellurics

Understanding deep continental structure and the seismotectonics of Deccan trap covered region has attained greater importance in recent years. For imaging the deep crustal structure, magnetotelluric (MT) investigations have been carried out along three long profiles viz. Guhagarh–Sangole (GS), Sangole–Partur (SP), Edlabad–Khandwa (EK) and one short profile along Nanasi–Mokhad (NM). The results of GS, SP and NM profiles show that the traps lie directly over high resistive basement with thin inter-trappean sediments, where large thickness of sediments, of the order of 1.5–2.0 km, has been delineated along EK profile across Narmada–Son–Lineament zone. The basement is intersected by faults/fractures, which are clearly delineated as narrow steep conducting features at a few locations. The conducting features delineated along SP profile are also seen from the results of aeromagnetic anomalies. Towards the southern part of the profile, these features are spatially correlated with Kurduwadi rift proposed earlier from gravity studies. Apart from the Kurduwadi rift extending to deep crustal levels, the present study indicates additional conductive features in the basement. The variation in the resistivity along GS profile can be attributed to crustal block structure in Koyna region. Similar block structure is also seen along NM profile.Deccan trap thickness, based on various geophysical methods, varies gradually from 1.8 km towards west to 0.3 km towards the east. While this is the general trend, a sharp variation in the thickness of trap is observed near Koyna. The resistivity of the trap is more (150–200 Ω m) towards the west as compared to the east (50–60 Ω m) indicating more compact or denser nature for the basalt towards west. The upper crust is highly resistive (5000–10,000 Ω m), and the lower crust is moderately resistive (500–1000 Ω m). In the present study, seismotectonics of the region is discussed based on the regional geoelectrical structure with lateral variation in the resistivity of the basement and presence of anomalous conductors in the crust.     

Comparison of gravity anomaly between mature and immature intra-oceanic subduction zones in the western Pacific

The western Pacific hosts major subduction systems such as Izu–Bonin–Mariana and Tonga–Kermadec, but also less conspicuous systems such as Yap, Mussau and Hjort trenches which constitute the young, incomplete, or ultraslow-member in the evolutionary spectrum of subduction zones. We used satellite-derived gravity data to compare well-developed and immature subduction systems. It is shown that at spatial resolution > 10–20 km or so, the satellite data have accuracy comparable to ship-board gravity measurements over intra-oceanic subduction zones. In the isostatic residual gravity anomaly map, the width of non-isostatically-compensated region of the mature subduction zones is much wider than that of immature ones. More importantly, when the gravitational attraction due to seafloor is removed, a large difference exists between the mature and immature subduction zones in the overriding plate side. Mature subduction zones exhibit broad low gravity anomalies of ~ 200–250 mGal centered at distances of 150–200 km from the trench which are not found over immature subduction zones. The cause of the broad low gravity anomalies over mature subduction zones is debatable due to lack of information on the deep crust and upper mantle structure and property. We discuss the following four causes: (1) serpentinization of the upper mantle beneath the forearc; (2) presence of partial melt in the mantle wedge caused by release of volatiles from the slab, frictional heating and distributed by mantle circulation; (3) difference in density structure between the overriding and subducting plates caused by difference in age and thermal structures with and without compositional stratification between crust and mantle; and (4) anomalous thickness of the arc not explained by isostasy. Our analysis suggests that serpentinization cannot explain the observed gravity anomaly which appears ~ 150–200 km from the trench. Although the extent and distribution of partial melt within the mantle wedge remain in question, to our best estimate, partial melting contributes little (< 50 mGal) to the total negative gravity anomaly. The difference in density structure reflecting temperature difference can only explain less than half of the low gravity anomaly. The sinking of lighter crustal material produces a large negative anomaly in the forearc but its location does not match the observed gravity anomaly. It appears that one cannot explain the total difference in gravity anomaly without invoking anomalous thickness of the arc. Although we could not identify the sole or combination of factors that give rise to the low gravity anomaly in mature subduction zones, the comparison of gravity anomalies between mature and immature subduction zones is likely to provide an important constraint for understanding the evolution and structure of subduction zones as more complementary evidences become available.     

Low-temperature thermochronological record of exhumation of the Bitterroot metamorphic core complex, northern Cordilleran Orogen

The Bitterroot metamorphic core complex is an exhumed, mid-crustal, plutonic–metamorphic complex that formed during crustal thickening and subsequent extension in the hinterland of the North American Cordilleran Orogen, in the northern Idaho batholith region. Extension was accommodated mainly on the Bitterroot mylonite zone, a 500–1500-m-thick shear zone that deforms granitic intrusive rocks as young as 53–52 Ma, as well as older high-grade metamorphic rocks and plutons. Exhumation of the core complex, in Eocene time, is marked in the shear zone by the transition from amphibolite-facies mylonitization, to greenschist-facies mylonitization, chloritic brecciation, to brittle faulting that progressed from shallower crustal levels in the west to deeper crustal levels in the east from ca. 53 –30 Ma based on U–Pb, Ar–Ar, and fission-track data. Apatite and zircon fission-track data record the lower-temperature part of the exhumation history and help define when the shear zone became inactive, as well as the transition from rapid, core complex-style extension to slower basin-and-range-style extension. They indicate that the western part of the complex was exhumed to within 1–2 km of the surface by 48–45 Ma, while the eastern part of the complex was still at amphibolite-facies conditions and that the eastern part of the complex was not exhumed below 60 °C until after 30 Ma. Younger apatite fission-track ages (≤26 Ma) on the eastern range front of the Bitterroot Mountains suggest that the present topographic expression of the mylonite front was due to Miocene high-angle faulting and widening of the Bitterroot Valley.     

Evidence for Proterozoic Collision from Airborne Magnetic and Gravity Studies in Southern Granulite Terrain, India and Signatures of Recent Tectonic Activity in the Palghat Gap

The composite airborne total intensity map of the Southern Granulite Terrain (SGT) at an average elevation of 7000' (≈ 2100 m) shows bands of bipolar regional magnetic anomalies parallel to the structural trends suggesting the distribution of mafic/ultramafic rocks that are controlled by regional structures/shear zones and thrusts in this region. The spectrum and the apparent susceptibility map computed from the observed airborne magnetic anomalies provide bands of high susceptibility zones in the upper crust associated with known shear zones/thrusts such as Transition Zone, Moyar-Bhavani and Palghat-Cauvery Shear Zones (MBSZ and PCSZ). The quantitative modelling of magnetic anomalies across Transition Zone, MBSZ and PCSZ suggest the presence of mafic rocks of susceptibility (1.5-4.0 × 10⁻³ CGS units) in upper crust from 8-10 km extending up to about 21-22 km, which may represent the level of Curie point geotherm as indicated by high upper mantle heat flow in this section.Two sets of paired gravity anomalies in SGT and their modelling with seismic constraints suggest gravity highs and lows to be caused by high density mafic rocks along Transition Zone and Cauvery Shear Zone (CSZ) in the upper crust at depth of 6-8 km and crustal thickening of 45-46 km south of them, respectively. High susceptibility and high density rocks (2.8 g/cm³) along these shear zones supported by high velocity, high conductivity and tectonic settings suggest lower crustal mafic/ultramafic granulite rocks thrusted along them. These signatures with lower crustal rocks of metamorphic ages of 2.6-2.5 Ga north of PCSZ and Neoproterozoic period (0.6-0.5 Ga) south of it suggest that the SGT represents mosaic of accreted crust due to compression and thrusting. These observations along with N-verging thrusts and dipping reflectors from Dharwar Craton to SGT suggest two stages of N-S directed compression: (i) between Dharwar Craton and northern block of SGT during 2.6-2.5 Ga with Transition Zone and Moyar Shear towards the west as thrust, and (ii) between northern and southern blocks of SGT with CSZ as collision zone and PCSZ as thrust during Neoproterozoic period (0.6-0.5 Ga). The latter event may even represent just a compressive phase without any collision related to Pan-African event. The proposed sutures in both these cases separate gravity highs and lows of paired gravity anomalies towards north and south, respectively. The magnetic anomalies and causative sources related to Moyar Shear, MBSZ and PCSZ join with those due to Transition Zone, Mettur and Gangavalli Shears in their eastern parts, respectively to form an arcuate-shaped diffused collision zone during 2.6-2.5 Ga.Most of the Proterozoic collision zones are highlands/plateaus but the CSZ also known as the Palghat Gap represents a low lying strip of 80-100 km width, which however, appears to be related to recent tectonic activities as indicated by high upper mantle heat flow and thin crust in this section. It is supported by low density, low velocity and high conductive layer under CSZ and seismic activity in this region as observed in case of passive rift valleys. They may be caused by asthenospheric upwarping along pre-existing faults/thrusts (MBSZ and PCSZ) due to plate tectonic forces after the collision of Indian and Eurasian plates since Miocene time.     

Seismic and geological structure of the crust in the transition from Baltica to Palaeozoic Europe in SE Poland—CELEBRATION 2000 experiment, profile CEL02

The CELEBRATION 2000 together with the earlier POLONAISE'97 deep seismic sounding experiments was aimed at the recognition of crustal structure in the border zone between the Precambrian East European Craton (Baltica) and Palaeozoic Europe. The CEL02 profile of the CELEBRATION family is a 400-km long SW–NE transect, running in Poland from the Upper Silesia Block (USB), across the Małopolska Block (MB) and the Trans-European Suture Zone (TESZ) to the East European Craton (EEC). The structure along CEL02 was interpreted using both 2D tomography and forward ray-tracing techniques as well as 2D gravity modelling.The crustal thickness along CEL02 varies from 32–35 km in the USB to 45–47 km beneath the TESZ and the EEC. The USB is a clearly distinctive crustal block with the characteristic high velocity lower crust (7.1–7.2 km/s), interpreted as a fragment of Gondwana. The Kraków–Lubliniec Fault is a terrane boundary produced by soft docking of the USB with the MB. The Małopolska crust fundamentally differs from the USB and has a strong connection with Baltica. It is a transitional, 150- to 200-km wide unit composed of the extended Baltican lower crust and the overlying low velocity (5.15–5.9 km/s) Neoproterozoic metasediments in the up to 18-km thick upper crust. The Łysogóry Unit has its crustal structure identical with that of Małopolska, thus it is connected with Baltica and cannot be interpreted as a Gondwana-derived terrane. Higher velocity and density bodies found below the Mazovia–Lublin Graben at a depth of 12 km and at the base of the lower crust, might be a result of mantle-derived mafic intrusions accompanying the extension of Baltica. By the preliminary 2D gravity modelling, we have reconfirmed the need for considering the increased TESZ mantle density in comparison to the EEC and USB mantle.     

Depths to density discontinuities beneath the Adamawa Plateau region, Central Africa, from spectral analyses of new and existing gravity data

New gravity data from the Adamawa Uplift region of Cameroon have been integrated with existing gravity data from central and western Africa to examine variations in crustal structure throughout the region. The new data reveal steep northeast-trending gradients in the Bouguer gravity anomalies that coincide with the Sanaga Fault Zone and the Foumban Shear Zone, both part of the Central African Shear Zone lying between the Adamawa Plateau and the Congo Craton. Four major density discontinuities in the lithosphere have been determined within the lithosphere beneath the Adamawa Uplift in central Cameroon using spectral analysis of gravity data: (1) 7–13 km; (2) 19–25 km; (3) 30–37 km; and (4) 75–149 km. The deepest density discontinuities determined at 75–149 km depth range agree with the presence of an anomalous low velocity upper mantle structure at these depths deduced from earlier teleseismic delay time studies and gravity forward modelling. The 30–37 km depths agree with the Moho depth of 33 km obtained from a seismic refraction experiment in the region. The intermediate depth of 20 km obtained within region D may correspond to shallower Moho depth beneath parts of the Benue and Yola Rifts where seismic refraction data indicate a crustal thickness of 23 km. The 19–20 km depths and 8–12 km depths estimated in boxes encompassing the Adamawa Plateau and Cameroon Volcanic Line may may correspond to mid-crustal density contrasts associated with volcanic intrusions, as these depths are less than depths of 25 and 13 km, respectively, in the stable Congo Craton to the south.     

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.     

Gravity Field and Subsurface Geometry of the Kalpatta Granite, South India and the Tectonic Significance

The Kalpatta granite, of Pan-African age occurs in the southern granulite terrain of Peninsular India. Bouguer anomaly map of the Kalpatta and adjoining areas reveals a gravity low of 8–10 mGal centered over the Kalpatta granite and a minor low of 4–6 mGal over the adjacent Ambalavayal granite pluton. The residual anomaly map prepared for the Kalpatta granite has been utilized to obtain depth extent and 3-D geometry of the pluton. The analysis suggests that the Kalpatta granite is an elliptical and somewhat pear shaped body with horizontal dimensions of 6–11 km and extending to a depth of 6.5km, has steeply inward dipping contacts and the shape seen on the surface continues throughout its depth. The smooth oval shape could indicate low ductility contrast, deeper level of emplacement and permissive nature of the pluton. The 3-D depth model indicates an oblique section with much deeper levels exposed in the south, in other words, the crustal block encompassing the pluton has suffered a NNW tilt during uplift after the emplacement. It is further inferred that there was no post intrusive shape modification, the NW tilting of the region and denudation gave rise to the present outcrop pattern of the body.     

Gravity signals from the lithosphere in the Central European Basin System

We study the gravity signals from different depth levels in the lithosphere of the Central European Basin System (CEBS). The major elements of the CEBS are the Northern and Southern Permian Basins which include the Norwegian–Danish Basin (NDB), the North-German Basin (NGB) and the Polish Trough (PT). An up to 10 km thick sedimentary cover of Mesozoic–Cenozoic sediments, hides the gravity signal from below the basin and masks the heterogeneous structure of the consolidated crust, which is assumed to be composed of domains that were accreted during the Paleozoic amalgamation of Europe. We performed a three-dimensional (3D) gravity backstripping to investigate the structure of the lithosphere below the CEBS.Residual anomalies are derived by removing the effect of sediments down to the base of Permian from the observed field. In order to correct for the influence of large salt structures, lateral density variations are incorporated. These sediment-free anomalies are interpreted to reflect Moho relief and density heterogeneities in the crystalline crust and uppermost mantle. The gravity effect of the Moho relief compensates to a large extent the effect of the sediments in the CEBS and in the North Sea. Removal of the effects of large-scale crustal inhomogeneities shows a clear expression of the Variscan arc system at the southern part of the study area and the old crust of Baltica further north–east. The remaining residual anomalies (after stripping off the effects of sediments, Moho topography and large-scale crustal heterogeneities) reveal long wavelength anomalies, which are caused mainly by density variations in the upper mantle, though gravity influence from the lower crust cannot be ruled out. They indicate that the three main subbasins of the CEBS originated on different lithospheric domains. The PT originated on a thick, strong and dense lithosphere of the Baltica type. The NDB was formed on a weakened Baltica low-density lithosphere formed during the Sveco-Norwegian orogeny. The major part of the NGB is characterized by high-density lithosphere, which includes a high-velocity lower crust (relict of Baltica passive margin) overthrusted by the Avalonian terrane. The short wavelength pattern of the final residuals shows several north–west trending gravity highs between the Tornquist Zone and the Elbe Fault System. The NDB is separated by a gravity low at the Ringkøbing–Fyn high from a chain of positive anomalies in the NGB and the PT. In the NGB these anomalies correspond to the Prignitz (Rheinsberg anomaly), the Glueckstadt and Horn Graben, and they continue further west into the Central Graben, to join with the gravity high of the Central North Sea.     

Structure of the crust and uppermost mantle of the Yanshan Belt and adjacent regions at the northeastern boundary of the North China Craton from Rayleigh Wave Dispersion Analysis

We have used Rayleigh wave dispersion analysis and inversion to produce a high resolution S-wave velocity imaging profile of the crust and uppermost mantle structure beneath the northeastern boundary regions of the North China Craton (NCC). Using waveform data from 45 broadband NCISP stations, Rayleigh wave phase velocities were measured at periods from 10 to 48 s and utilized in subsequent inversions to solve for the S-wave velocity structure from 15 km down to 120 km depth. The inverted lower crust and uppermost mantle velocities, about 3.75 km/s and 4.3 km/s on average, are low compared with the global average. The Moho was constrained in the depth range of 30–40 km, indicating a typical crustal thickness along the profile. However, a thin lithosphere of no more than 100 km was imaged under a large part of the profile, decreasing to only ~ 60 km under the Inner Mongolian Axis (IMA) where an abnormally slow anomaly was observed below 60 km depth. The overall structural features of the study region resemble those of typical continental rift zones and are probably associated with the lithospheric reactivation and tectonic extension widespread in the eastern NCC during Mesozoic–Cenozoic time. Distinctly high velocities, up to ~ 4.6 km/s, were found immediately to the south of the IMA beneath the northern Yanshan Belt (YSB), extending down to > 100-km depth. The anomalous velocities are interpreted as the cratonic lithospheric lid of the region, which may have not been affected by the Mesozoic–Cenozoic deformation process as strongly as other regions in the eastern NCC. Based on our S-wave velocity structural image and other geophysical observations, we propose a possible lithosphere–asthenosphere interaction scenario at the northeastern boundary of the NCC. We speculate that significant undulations of the base of the lithosphere, which might have resulted from the uneven Mesozoic–Cenozoic lithospheric thinning, may induce mantle flows concentrating beneath the weak IMA zone. The relatively thick lithospheric lid in the northern YSB may serve as a tectonic barrier separating the on-craton and off-craton regions into different upper mantle convection systems at the present time.     

Local thickening of the Cascadia forearc crust and the origin of seismic reflectors in the uppermost mantle

Seismic reflection profiles from three different surveys of the Cascadia forearc are interpreted using P wave velocities and relocated hypocentres, which were both derived from the first arrival travel time inversion of wide-angle seismic data and local earthquakes. The subduction decollement, which is characterized beneath the continental shelf by a reflection of 0.5 s duration, can be traced landward into a large duplex structure in the lower forearc crust near southern Vancouver Island. Beneath Vancouver Island, the roof thrust of the duplex is revealed by a 5–12 km thick zone, identified previously as the E reflectors, and the floor thrust is defined by a short duration reflection from a < 2-km-thick interface at the top of the subducting plate. We show that another zone of reflectors exists east of Vancouver Island that is approximately 8 km thick, and identified as the D reflectors. These overlie the E reflectors; together the two zones define the landward part of the duplex. The combined zones reach depths as great as 50 km. The duplex structure extends for more than 120 km perpendicular to the margin, has an along-strike extent of 80 km, and at depths between 30 km and 50 km the duplex structure correlates with a region of anomalously deep seismicity, where velocities are less than 7000 m s^− 1. We suggest that these relatively low velocities indicate the presence of either crustal rocks from the oceanic plate that have been underplated to the continent or crustal rocks from the forearc that have been transported downward by subduction erosion. The absence of seismicity from within the E reflectors implies that they are significantly weaker than the overlying crust, and the reflectors may be a zone of active ductile shear. In contrast, seismicity in parts of the D reflectors can be interpreted to mean that ductile shearing no longer occurs in the landward part of the duplex. Merging of the D and E reflectors at 42–46 km depth creates reflectivity in the uppermost mantle with a vertical thickness of at least 15 km. We suggest that pervasive reflectivity in the upper mantle elsewhere beneath Puget Sound and the Strait of Georgia arises from similar shear zones.     

Crustal structure of the northeastern margin of the Tibetan plateau from the Songpan-Ganzi terrane to the Ordos basin

The 1000-km-long Darlag–Lanzhou–Jingbian seismic refraction profile is located in the NE margin of the Tibetan plateau. This profile crosses the northern Songpan-Ganzi terrane, the Qinling-Qilian fold system, the Haiyuan arcuate tectonic region, and the stable Ordos basin. The P-wave and S-wave velocity structure and Poisson's ratios reveal many significant characteristics in the profile. The crustal thickness increases from northeast to southwest. The average crustal thickness observed increases from 42 km in the Ordos basin to 63 km in the Songpan-Ganzi terrane. The crust becomes obviously thicker south of the Haiyuan fault and beneath the West-Qinlin Shan. The crustal velocities have significant variations along the profile. The average P-wave velocities for the crystalline crust vary between 6.3 and 6.4 km/s. Beneath the Songpan-Ganzi terrane, West-Qinling Shan, and Haiyuan arcuate tectonic region P-wave velocities of 6.3 km/s are 0.15 km/s lower than the worldwide average of 6.45 km/s. North of the Kunlun fault, with exclusion of the Haiyuan arcuate tectonic region, the average P-wave velocity is 6.4 km/s and only 0.5 km/s lower than the worldwide average. A combination of the P-wave velocity and Poisson's ratio suggests that the crust is dominantly felsic in composition with an intermediate composition at the base. A mafic lower crust is absent in the NE margin of the Tibetan plateau from the Songpan-Ganzi terrane to the Ordos basin. There are low velocity zones in the West-Qinling Shan and the Haiyuan arcuate tectonic region. The low velocity zones have low S-wave velocities and high Poisson's ratios, so it is possible these zones are due to partial melting. The crust is divided into two layers, the upper and the lower crust, with crustal thickening mainly in the lower crust as the NE Tibetan plateau is approached. The results in the study show that the thickness of the lower crust increases from 22 to 38 km as the crustal thickness increases from 42 km in the Ordos basin to 63 km in the Songpan-Ganzi terrane south of the Kunlun fault. Both the Conrad discontinuity and Moho in the West-Qinling Shan and in the Haiyuan arcuate tectonic region are laminated interfaces, implying intense tectonic activity. The arcuate faults and large earthquakes in the Haiyuan arcuate tectonic region are the result of interaction between the Tibetan plateau and the Sino–Korean and Gobi Ala Shan platforms.     

Gravity anomalies and associated tectonic features over the Indian Peninsular Shield and adjoining ocean basins

A combined gravity map over the Indian Peninsular Shield (IPS) and adjoining oceans brings out well the inter-relationships between the older tectonic features of the continent and the adjoining younger oceanic features. The NW–SE, NE–SW and N–S Precambrian trends of the IPS are reflected in the structural trends of the Arabian Sea and the Bay of Bengal suggesting their probable reactivation. The Simple Bouguer anomaly map shows consistent increase in gravity value from the continent to the deep ocean basins, which is attributed to isostatic compensation due to variations in the crustal thickness. A crustal density model computed along a profile across this region suggests a thick crust of 35–40 km under the continent, which reduces to 22/20–24 km under the Bay of Bengal with thick sediments of 8–10 km underlain by crustal layers of density 2720 and 2900/2840 kg/m³. Large crustal thickness and trends of the gravity anomalies may suggest a transitional crust in the Bay of Bengal up to 150–200 km from the east coast. The crustal thickness under the Laxmi ridge and east of it in the Arabian Sea is 20 and 14 km, respectively, with 5–6 km thick Tertiary and Mesozoic sediments separated by a thin layer of Deccan Trap. Crustal layers of densities 2750 and 2950 kg/m³ underlie sediments. The crustal density model in this part of the Arabian Sea (east of Laxmi ridge) and the structural trends similar to the Indian Peninsular Shield suggest a continent–ocean transitional crust (COTC). The COTC may represent down dropped and submerged parts of the Indian crust evolved at the time of break-up along the west coast of India and passage of Reunion hotspot over India during late Cretaceous. The crustal model under this part also shows an underplated lower crust and a low density upper mantle, extending over the continent across the west coast of India, which appears to be related to the Deccan volcanism. The crustal thickness under the western Arabian Sea (west of the Laxmi ridge) reduces to 8–9 km with crustal layers of densities 2650 and 2870 kg/m³ representing an oceanic crust.     

P-wave velocity structure under the Changbaishan volcanic region, NE China, from wide-angle reflection and refraction data

We present velocity models determined by inverting refracted and reflected arrivals along two active source lines in the Changbaishan volcanic region, NE China. We resolve a prominent low-velocity zone (LVZ) in the crust, with velocities as low as 5.4 km/s. Away from the LVZ, the velocity gradients in the crust are relatively smooth, with average P-wave velocities of about 6.0–6.5 km/s. The Moho is at about 35 km depth, thickening to about 40 km under the Tianchi volcano, and thinning to about 30 km under the LVZ. The LVZ is located about 30–60 km to the north of the summit of the Tianchi volcano (the most recently active volcano in the region), is about 30–75 km in north–south extent, is at most 35 km in east–west extent, and is in the depth range of about 10–25 km below the surface. We use these results to constrain receiver function inversions, and show that the receiver functions in the region are compatible with our findings. With these data alone, the significance of the LVZ in non-unique, although we do not see any evidence to support the presence of partial melt in the crust, and favor the interpretation that the LVZ indicates a residual crustal magma chamber.     

A new tectonic model for the Laurentia−Avalonia−Baltica sutures in the North Sea: A case study along MONA LISA profile 3

We present a new model for the lithospheric structure of the transitions between Laurentia, Avalonia and Baltica in the North Sea, northwestern Europe based on 2¾D potential field modelling of MONA LISA profile 3 across the Central Graben, with constraints from seismic P-wave velocity models and the crustal normal incidence reflection section along the profile. The model shows evidence for the presence of upper-and lower Palaeozoic sedimentary rocks as well as differences in crustal structure between the palaeo-continents Laurentia, Avalonia and Baltica. Our new model, together with previous results from transformations of the gravity and magnetic fields, demonstrates correlation between crustal magnetic domains along the profile and the terrane affinity of the crust. This integrated interpretation indicates that a 150 km wide zone, characterized by low-grade metamorphosis and oblique thrusting of Avalonia crust over Baltica lower crust, is characteristic for the central North Sea area. The magnetic susceptibility and the density across the Coffee Soil Fault range from almost zero and 2715 kg/m³ in Avalonia crust to 0.05 SI and 2775 kg/m³ in Baltica crust. The model of MONA LISA profile 3 indicates that the transition between Avalonia and Baltica is located beneath the Central Graben with a ramp–flat–ramp geometry. Our results indicate that the initial rifting of the Central Graben and the Viking Graben was controlled by the location of the Caledonian collisional suture, located at the Coffee Soil Fault, and that the deep crustal part of Baltica extends further to the west than hitherto believed.     

Crustal structure of mainland China from deep seismic sounding data

Since 1958, about ninety seismic refraction/wide angle reflection profiles, with a cumulative length of more than sixty thousand kilometers, have been completed in mainland China. We summarize the results in the form of (1) a new contour map of crustal thickness, (2) fourteen representative crustal seismic velocity–depth columns for various tectonic units, and, (3) a Pn velocity map. We found a north–south-trending belt with a strong lateral gradient in crustal thickness in central China. This belt divides China into an eastern region, with a crustal thickness of 30–45 km, and a western region, with a thickness of 45–75 km. The crust in these two regions has experienced different evolutionary processes, and currently lies within distinct tectonic stress fields. Our compilation finds that there is a high-velocity (7.1–7.4 km/s) layer in the lower crust of the stable Tarim basin and Ordos plateau. However, in young orogenic belts, including parts of eastern China, the Tianshan and the Tibetan plateau, this layer is often absent. One exception is southern Tibet, where the presence of a high-velocity layer is related to the northward injection of the cold Indian plate. This high-velocity layer is absent in northern Tibet. In orogenic belts, there usually is a low-velocity layer (LVL) in the crust, but in stable regions this layer seldom exists. The Pn velocities in eastern China generally range from 7.9 to 8.1 km/s and tend to be isotropic. Pn velocities in western China are more variable, ranging from 7.7 to 8.2 km/s, and may display azimuthal anisotropy.