P-wave tomography and dynamics of the crust and upper mantle beneath western Tibet

The continental collision between the Indian and Asian plates plays a key role in the geologic and tectonic evolution of the Tibetan plateau. In this article we present high-resolution tomographic images of the crust and upper mantle derived from a large number of high-quality seismic data from the ANTILOPE project in western Tibet. Both local and distant earthquakes were used in this study and 35,115 P-wave arrival times were manually picked from the original seismograms. Geological and geochemical results suggested that the subducting Indian plate has reached northward to the Lhasa terrane, whereas our new tomography shows that the Indian plate is currently sub-horizontal and underthrusting to the Jinsha river suture at depths of ~ 100 to ~ 250 km, suggesting that the subduction process has evolved over time. The Asian plate is also imaged clearly from the surface to a depth of ~ 100 km by our tomography, and it is located under the Tarim Basin north of the Altyn Tagh Fault. There is no obvious evidence to show that the Asian plate has subducted beneath western Tibet. The Indian and Asian plates are separated by a prominent low-velocity zone under northern Tibet. We attribute the low-velocity zone to mantle upwelling, which may account for the warm crust and upper mantle beneath that region, and thus explain the different features of magmatism between southern and northern Tibet. But the upwelling may not penetrate through the whole crust. We propose a revised geodynamic model and suggest that the high-velocity zones under Lhasa terrane may reflect a cold crust which has interrupted the crustal flow under the westernmost Tibetan plateau.     

Mesoproterozoic granulites of the Shillong–Meghalaya Plateau: Evidence of westward continuation of the Prydz Bay Pan-African suture into Northeastern India

Thermo-chronologic considerations in the Australo-Antarctic domains suggest that the Prydz Bay Pan-African suture of East Antarctica continues westward into India. However, the location of the suture within Eastern India has so far been uncertain because of a lack of adequate thermo-chronological information. In this study, electron microprobe (EPMA) monazite dates and mineral paragenesis of granulite facies metapelites are reported from two areas of the Shillong–Meghalaya gneissic complex (SMGC), a crustal block located in the extreme northeast of the Indian shield close to the Australo-Antarctic block in Neoproterozoic-Cambrian paleomagnetic reconstructions of the Rodinia supercontinent. In the Garo-Goalpara Hills region, a well constrained Mesoproterozoic age of 1596 ± 15 Ma (n = 103) is correlated with a counterclockwise pressure–temperature path with near peak conditions of 7–8 kbar and 850 °C. Rare matrix monazite rims record younger ages (1032–1273 Ma). At Sonapahar region, 50 km ESE of Garo-Goalpara Hills, homogenous monazite grains in granulite facies metapelites yield EPMA dates tightly clustered at 500 ± 14 Ma (n = 36) irrespective of their textural setting in a well-annealed mineral matrix. In a few zoned monazite grains, the cores yield older ages of 1078 ± 31 Ma (n = 10) and 1472 ± 38 Ma (n = 13). The 500 Ma date corresponds with the ca. 880–480 Ma Rb–Sr dates of porphyritic granites that predominantly intruded the east-central part of the SMGC. We propose that the progressively eastward dominance of Cambro-Ordovician ages in the SMGC indicates a Pan-African final amalgamation of the Indian plate with the Australo-Antarctic plate and a northward extension of the Prydz Bay suture through the SMGC, with the western boundary of the suture possibly located between the Garo-Goalpara Hills and Sonapahar areas.     

Co and postseismic characteristics of Indian sub-continent in response to the 2004 Sumatra earthquake

The M_w 9.3 Sumatra earthquake of December 26, 2004 caused extensive coseismic displacements globally, measurements of which were made essentially using modern geodetic techniques. This earthquake induced considerable perturbation in stress distribution as far as ∼8000 km away from the epicenteral region, which is tending to relax to its normal rates as seen from postseismic transient deformation. The monitoring of crustal displacements from strategically located sites using GPS provides coseismic as well as postseismic deformation that facilitates the understanding of the fault geometry, elastic thickness, postseismic relaxation mechanisms, rheology and earthquake recurrence time interval.We investigated coseismic and postseismic GPS derived displacements in Indian region together with the GPS data collected from Andaman and Sumatra region. It is found that while EW displacements are significantly large in peninsular India, those in the region to the north of Central India Tectonic Zone (CITZ) are relatively small. We could delineate the postseismic transients from position time series and interpreted them in terms of viscoelastic relaxation. It is inferred that the postseismic deformation is characterized by a power-law viscoelastic flow in the mantle. In Indian peninsula region, the timescale parameter of the exponential decay (τ = 250 days) would require an extremely low viscosity for the upper mantle. Relying on the prevailing coseismic and postseismic displacement fields, the present study also reflects upon the contemporary litho-tectonics of the Indian sub-continent.     

The longest voyage: Tectonic,magmatic, and paleoclimatic evolution of the Indian plate during its northward flight from Gondwana to Asia

The tectonic evolution of the Indian plate, which started in Late Jurassic about 167 million years ago (~ 167 Ma) with the breakup of Gondwana, presents an exceptional and intricate case history against which a variety of plate tectonic events such as: continental breakup, sea-floor spreading, birth of new oceans, flood basalt volcanism, hotspot tracks, transform faults, subduction, obduction, continental collision, accretion, and mountain building can be investigated. Plate tectonic maps are presented here illustrating the repeated rifting of the Indian plate from surrounding Gondwana continents, its northward migration, and its collision first with the Kohistan–Ladakh Arc at the Indus Suture Zone, and then with Tibet at the Shyok–Tsangpo Suture. The associations between flood basalts and the recurrent separation of the Indian plate from Gondwana are assessed. The breakup of India from Gondwana and the opening of the Indian Ocean is thought to have been caused by plate tectonic forces (i.e., slab pull emanating from the subduction of the Tethyan ocean floor beneath Eurasia) which were localized along zones of weakness caused by mantle plumes (Bouvet, Marion, Kerguelen, and Reunion plumes). The sequential spreading of the Southwest Indian Ridge/Davie Ridge, Southeast Indian Ridge, Central Indian Ridge, Palitana Ridge, and Carlsberg Ridge in the Indian Ocean were responsible for the fragmentation of the Indian plate during the Late Jurassic and Cretaceous times. The Réunion and the Kerguelen plumes left two spectacular hotspot tracks on either side of the Indian plate. With the breakup of Gondwana, India remained isolated as an island continent, but reestablished its biotic links with Africa during the Late Cretaceous during its collision with the Kohistan–Ladakh Arc (~ 85 Ma) along the Indus Suture. Soon after the Deccan eruption, India drifted northward as an island continent by rapid motion carrying Gondwana biota, about 20 cm/year, between 67 Ma to 50 Ma; it slowed down dramatically to 5 cm/year during its collision with Asia in Early Eocene (~ 50 Ma). A northern corridor was established between India and Asia soon after the collision allowing faunal interchange. This is reflected by mixed Gondwana and Eurasian elements in the fossil record preserved in several continental Eocene formations of India. A revised India–Asia collision model suggests that the Indus Suture represents the obduction zone between India and the Kohistan–Ladakh Arc, whereas the Shyok-Suture represents the collision between the Kohistan–Ladakh arc and Tibet. Eventually, the Indus–Tsangpo Zone became the locus of the final India–Asia collision, which probably began in Early Eocene (~ 50 Ma) with the closure of Neotethys Ocean. The post-collisional tectonics for the last 50 million years is best expressed in the evolution of the Himalaya–Tibetan orogen. The great thickness of crust beneath Tibet and Himalaya and a series of north vergent thrust zones in the Himalaya and the south-vergent subduction zones in Tibetan Plateau suggest the progressive convergence between India and Asia of about 2500 km since the time of collision. In the early Eohimalayan phase (~ 50 to 25 Ma) of Himalayan orogeny (Middle Eocene–Late Oligocene), thick sediments on the leading edge of the Indian plate were squeezed, folded, and faulted to form the Tethyan Himalaya. With continuing convergence of India, the architecture of the Himalayan–Tibetan orogen is dominated by deformational structures developed in the Neogene Period during the Neohimalayan phase (~ 21 Ma to present), creating a series of north-vergent thrust belt systems such as the Main Central Thrust, the Main Boundary Thrust, and the Main Frontal Thrust to accommodate crustal shortening. Neogene molassic sediment shed from the rise of the Himalaya was deposited in a nearly continuous foreland trough in the Siwalik Group containing rich vertebrate assemblages. Tomographic imaging of the India–Asia orogen reveals that Indian lithospheric slab has been subducted subhorizontally beneath the entire Tibetan Plateau that has played a key role in the uplift of the Tibetan Plateau. The low-viscosity channel flow in response to topographic loading of Tibet provides a mechanism to explain the Himalayan–Tibetan orogen. From the start of its voyage in Southern Hemisphere, to its final impact with the Asia, the Indian plate has experienced changes in climatic conditions both short-term and long-term. We present a series of paleoclimatic maps illustrating the temperature and precipitation conditions based on estimates of Fast Ocean Atmospheric Model (FOAM), a coupled global climate model. The uplift of the Himalaya–Tibetan Plateau above the snow line created two most important global climate phenomena—the birth of the Asian monsoon and the onset of Pleistocene glaciation. As the mountains rose, and the monsoon rains intensified, increasing erosional sediments from the Himalaya were carried down by the Ganga River in the east and the Indus River in the west, and were deposited in two great deep-sea fans, the Bengal and the Indus. Vertebrate fossils provide additional resolution for the timing of three crucial tectonic events: India–KL Arc collision during the Late Cretaceous, India–Asia collision during the Early Eocene, and the rise of the Himalaya during the Early Miocene.     

Location of the pilot borehole for investigations of reservoir triggered seismicity at Koyna,India

Artificial water reservoir triggered earthquakes are now known to have occurred at over 120 sites globally. The part played by the reservoirs in triggering is not exactly known due to lack of near field observations of triggered earthquakes. Koyna, located near the west coast of India, where triggered earthquakes have been occurring since 1962 provides an excellent site for near field observations of the target M ≥ 2 earthquakes. A 6 borehole seismic network has been deployed recently in the Koyna region at depths of 981–1522 m to improve the hypocenter locations. During May–December 2015, a total of 1039 earthquakes of M_L ≥ 0.5 were located using the borehole seismic network. The region is also monitored through a dense network of 23 surface broad-band stations. Our analysis indicates a significant improvement in the estimation of absolute locations of earthquakes with errors of the order of ± 300 m, combining both the networks. Based on seismicity, and logistics, a block of 2 × 2 km² area has been chosen for drilling the first pilot borehole of ~ 3 km depth, where M ≥ 2 earthquakes have been occurring frequently since 2005.     

Deep electrical structure over the igneous arc of the Indo Burman Orogen in Sagaing province,Myanmar from magnetotelluric studies

Magnetotelluric studies over the igneous arc of the Indo Burman range in the Sagaing province of Myanmar have delineated the high resistivity Indian plate subducting westwards beneath the Burmese block to depths of 30 km and beyond. The thick moderately resistive (20–100 Ω m) layer overlying the subducting Indian plate may be due to the low resistivity sediments. The entire region is covered with prominent sedimentary layer with a conductance varying between 20 and 3000 S showing a general increase from the east to west, suggesting that their thickness increases toward the west. The large unsystematic variations in the conductance are indicative of the widely varying depositional environments and also possible vertical block movements during the course of their deposition. A west dipping low resistivity zone to the east of Burmese block seems to demarcate its eastern limit, suggesting the possibility of a hitherto unknown deep seated fault, which is also supported by the several earthquake foci located over this zone. The nature of the crustal movements over this fault is not immediately apparent. Possibility exists that the Sagaing fault is an en echelon fault and the present feature observed here is a part of this en echelon fault. The possibility of channel flows of the weakened rocks in the deep crust observed in the vicinity of the eastern Himalayan syntaxis may also cause such low resistivity zones.     

Present-day crustal motion around the Pamir Plateau from GPS measurements

We have used geodetic techniques to improve constraints on the crustal motion of the Pamir Plateau. Three campaigns of Global Position System data acquisition between 2011 and 2015 demonstrate that, in association with the India–Asia collision, a complex pattern of crustal motion exists in the Pamir Plateau. In a north–south direction from the Indian Plate to the Hazak Block, the crust has absorbed ~ 35 mm/yr of shortening, of which ~ 35% is distributed around the Hindu Kush region (~ 12 mm/yr), and another ~ 35% is taken up around the Alai Valley (also ~ 12 mm/yr). Global Position System measurements also show ~ 5 mm/yr of shortening between the Pamir Plateau and the Tajik Basin, whereas between the Pamir and the Tarim Basin, an ~ 10 mm/yr extension rate is observed. With respect to the stable Eurasian Plate, the Pamir rotates counterclockwise at a rate of ~ 1.822°Myr^− 1, with an Euler pole positioned about the west end of the Tajik Basin (37.03 ± 0.74°N, 65.89 ± 0.12°E). The strain rate field calculated from Global Position System velocities reveals that the crustal motion is consistent with localized deformation around the Hindu Kush and the Alai Valley, the latter representing a zone with strong shallow seismic activity.     

The plate coupling in the Tokai District,the Central Japan,inferred from the different data using triangular dislocation elements

We estimate interseismic coupling on the subducting plate interface in the Tokai area, central Japan, by inverting two geodetic data sets. The data record surface motion between March 1996 to May 2000; one represents vertical motion deduced from the leveling observations and the other is the horizontal velocity field deduced from GPS observations. In the inversion, we employed the analytical solutions of surface displacement due to a triangular dislocation element embedded in a homogeneous elastic half space in order to represent the curved plate interface. The vertical data show that the most strongly coupled portion of the subduction interface is concentrated beneath Omaezaki Cape, while the horizontal data show strongest coupling in the shallower region of the subducting plate interface. The estimated maximum value of coupling from the horizontal data is 40 mm/year, while that from vertical data is 25 mm/year.     

Paleomagnetic versus GPS determined tectonic rotation around eastern Himalayan syntaxis in East Asia

Northward indentation of the Indian Plate has brought about significant tectonic deformation into East Asia. A record of long-term tectonic deformation in this area for the past 50 M yr, particularly the vertical axis rotation, is available through paleomagnetic data. In order to depict rotational deformation in this area with respect to Eurasia, we compiled reliable paleomagnetic data sets from 79 localities distributed around eastern Himalayan syntaxis in East Asia. This record delineates that a zone affected by clockwise rotational deformation extends from the southern tip of the Chuan Dian Fragment to as far as the northwestern part of the Indochina Peninsula. A limited zone that experienced a significant amount of clockwise rotation after an initial India–Asia collision is now located at 23.5°N, 101°E, far away from an area (27.5°N, 95.5°E) where an intense rotational motion has been viewed by a snapshot of GPS measurements. This discrepancy in clockwise rotated positions is attributed to southeastward extrusion of the tectonic blocks within East Asia as a result of ongoing indentation of the Indian Plate. A quantitative comparison between the GPS and paleomagnetically determined clockwise rotation further suggests that following an initial India–Asia collision the crust at 30°N, 94°E paleoposition was subjected to southeastward displacement together with clockwise rotation, which eventually reached to present-day position of 23.5°N, 101°E, implying a crustal displacement of about 1000 km during the past 50 M yr.     

Strike-slip tectonics during the Neoproterozoic–Cambrian assembly of East Gondwana: Evidence from a newly discovered microcontinent in the Indian Ocean (Batavia Knoll)

The first dredge samples recovered from the Batavia Knoll in the eastern Indian Ocean demonstrate it is a rifted continental fragment that now lies submerged in 2 km of water between the Perth Abyssal Plain and Wharton Basin, ~ 1600 km offshore Western Australia. This knoll is the second microcontinent to have been discovered in this region, lying ~ 180 km to the NE of Gulden Draak Knoll; both knolls rifted from the Indian margin at c. 104–101 Ma during the breakup of East Gondwana. Dredged rocks include granite, granodiorite, charnockite and quartz-rich pegmatite, preserving igneous microstructure and biotite-schlieren. The composite nature of many samples suggests derivation from a felsic igneous complex. This basement is at least partly covered by Late Cretaceous sedimentary strata. Zircon grains from five representative granitoid samples primarily display modified oscillatory zoning and U-Pb data that spread along concordia over ~ 50–70 m.y., suggesting protracted magmatism and Pb-loss during the latest Neoproterozoic–early Cambrian (c. 580–530 Ma). Despite the span in U-Pb zircon ages, the primary igneous populations exhibit a relatively narrow range of initial Hf isotopic ratios (Hf_i = 0.281290–0.281814; εHf_i = − 41.40 to − 23.32). This evolved nature of the Hf signature suggests derivation of melts from ancient crust (T_DM2 = 4.05–2.92 Ga) with very limited juvenile input. Plate models based on sparse paleomagnetic data imply significant relative motion between India and Australia during the late Neoproterozoic–Cambrian. We suggest the new isotopic data from the Batavia Knoll granitic rocks, coupled with other regional geological constraints, are consistent with magmatic intrusion and subsequent disturbance of U-Pb systems along a late Neoproterozoic–Cambrian sinistral strike-slip margin related to Gondwana assembly.     

A Late Cretaceous and Cenozoic reconstruction of the Southwest Pacific region: Tectonics controlled by subduction and slab rollback processes

A Cenozoic tectonic reconstruction is presented for the Southwest Pacific region located east of Australia. The reconstruction is constrained by large geological and geophysical datasets and recalculated rotation parameters for Pacific–Australia and Lord Howe Rise–Pacific relative plate motion. The reconstruction is based on a conceptual tectonic model in which the large-scale structures of the region are manifestations of slab rollback and backarc extension processes. The current paradigm proclaims that the southwestern Pacific plate boundary was a west-dipping subduction boundary only since the Middle Eocene. The new reconstruction provides kinematic evidence that this configuration was already established in the Late Cretaceous and Early Paleogene. From ∼ 82 to ∼ 52 Ma, subduction was primarily accomplished by east and northeast-directed rollback of the Pacific slab, accommodating opening of the New Caledonia, South Loyalty, Coral Sea and Pocklington backarc basins and partly accommodating spreading in the Tasman Sea. The total amount of east-directed rollback of the Pacific slab that took place from ∼ 82 Ma to ∼ 52 Ma is estimated to be at least 1200 km. A large percentage of this rollback accommodated opening of the South Loyalty Basin, a north–south trending backarc basin. It is estimated from kinematic and geological constraints that the east–west width of the basin was at least ∼ 750 km. The South Loyalty and Pocklington backarc basins were subducted in the Eocene to earliest Miocene along the newly formed New Caledonia and Pocklington subduction zones. This culminated in southwestward and southward obduction of ophiolites in New Caledonia, Northland and New Guinea in the latest Eocene to earliest Miocene. It is suggested that the formation of these new subduction zones was triggered by a change in Pacific–Australia relative motion at ∼ 50 Ma. Two additional phases of eastward rollback of the Pacific slab followed, one during opening of the South Fiji Basin and Norfolk Basin in the Oligocene to Early Miocene (up to ∼ 650 km of rollback), and one during opening of the Lau Basin in the latest Miocene to Present (up to ∼ 400 km of rollback). Two new subduction zones formed in the Miocene, the south-dipping Trobriand subduction zone along which the Solomon Sea backarc Basin subducted and the north-dipping New Britain–San Cristobal–New Hebrides subduction zone, along which the Solomon Sea backarc Basin subducted in the west and the North Loyalty–South Fiji backarc Basin and remnants of the South Loyalty–Santa Cruz backarc Basin subducted in the east. Clockwise rollback of the New Hebrides section resulted in formation of the North Fiji Basin. The reconstruction provides explanations for the formation of new subduction zones and for the initiation and termination of opening of the marginal basins by either initiation of subduction of buoyant lithosphere, a change in plate kinematics or slab–mantle interaction.     

Geodetic and seismological investigation of crustal deformation near Izmir (Western Anatolia)

The Aegean region including western Turkey, mainland Greece, and the Hellenic Arc is the most seismological and geodynamical active domain in the Alpine Himalayan Belt. In this study, we processed 3 years of survey-mode GPS data and present the analysis of a combination of geodetic and seismological data around Izmir, which is the third most populated city in Turkey. The velocities obtained from 15 sites vary between 25 mm/yr and 28 mm/yr relative to the Eurasian plate. The power law exponent of earthquake size distribution (b-value) ranges from 0.8 to 2.8 in the Izmir region between 26.2°E and 27.2°E. The lowest b-value zones are found along Karaburun Fault (b = 0.8) and, between Seferihisar and Tuzla Faults (b = 0.8). A localized stress concentration is expected from numerical models of seismicity along geometrical locked fault patches. Therefore, areas with lowest b-values are considered to be the most likely location for a strong earthquake, a prediction that is confirmed by the 2005 M_w = 5.9 Seferihisar earthquake sequences, with epicentres located to the south of the Karaburun Fault. The north–south extension of the Izmir area is corroborated by extension rates up to 140 nanostrain/yr as obtained from our GPS data. We combined the 3-year GPS velocity field with the published velocity field to determine the strain rate pattern in the area. The spatial distribution of b-value reflects the normal background due to the tectonic framework and is corroborated by the geodetic data. b-Values correlate with strain pattern. This relationship suggests that decrease of b-values signifies accumulating strain.     

Tectonic drivers and the influence of the Kerguelen plume on seafloor spreading during formation of the early Indian Ocean

The Perth Abyssal Plain (PAP), located offshore southwest Australia, formed at the centre of Mesozoic East Gondwana breakup and Kerguelen plume activity. Despite its importance as a direct and relatively undisturbed recorder of this early spreading history, sparse geophysical data sets and lack of geological sampling hamper our understanding of the evolution of the PAP. This study combines new bathymetric profiles across the PAP with petrographic and geochemical data from the first samples ever to be dredged from the flanks of the Dirck Hartog Ridge (DHR), a prominent linear bathymetric feature in the central PAP, to better constrain the formation of the early Indian Ocean floor and the influence of the Kerguelen plume. Seafloor spreading in the PAP initiated at ~ 136 Ma with spreading observed to occur at (half) rates of ~ 35 mm/yr. Changes in spreading rate are difficult to discern after the onset of the Cretaceous Quiet Zone at ~ 120 Ma, but an increase in seafloor roughness towards the centre of the PAP likely resulted from a half-spreading rate decrease from 35 mm/yr (based on magnetic reversals) to ~ 24 mm/yr at ~ 114 Ma. Exhumed gabbro dredged from the southernmost dredge site of the DHR supports a further slowdown of spreading immediately prior to full cessation at ~ 102 Ma. The DHR exhibits a high relief ridge axis and distinctive asymmetry that is unusual compared to most active or extinct spreading centres. The composition of mafic volcanic samples varies along the DHR, from sub-alkaline dolerites with incompatible element concentrations most similar to depleted-to-normal mid-ocean ridge basalts in the south, to alkali basalts similar to ocean island basalts in the north. Therefore, magma sources and degrees of partial melting varied in space and time. It is likely that the alkali basalts are a manifestation of later excess volcanism, subsequent to or during the cessation of spreading. In this case, enriched signatures may be attributed tectonic drivers and melting of a heterogeneous mantle, or to an episodic influence of the Kerguelen plume over distances greater than 1000 km. We also investigate possible scenarios to explain how lower crustal rocks were emplaced at the crest of the southern DHR. Our results demonstrate the significance of regional tectonic plate motions on the formation and deformation of young ocean crust, and provide insight into the unique DHR morphology.     

Evidence of Mid-ocean ridge and shallow subduction forearc magmatism in the Nagaland-Manipur ophiolites,northeast India: constraints from mineralogy and geochemistry of gabbros and associated mafic dykes

We discuss here the mineralogical and geochemical characteristics of mafic intrusive rocks from the Nagaland-Manipur Ophiolites (NMO) of Indo-Myanmar Orogenic Belt, northeast India to define their mantle source and tectonic environment. Mafic intrusive sequence in the NMO is characterized by hornblende-free (type-I) and hornblende-bearing (type-II) rocks. The type-I is further categorized as mafic dykes (type-Ia) of tholeiitic N-MORB composition, having TiO₂ (0.72–1.93 wt.%) and flat REE patterns (La_N/Yb_N = 0.76–1.51) and as massive gabbros (type-Ib) that show alkaline E-MORB affinity, having moderate to high Ti content (TiO₂ = 1.18 to 1.45 wt.%) with strong LREE-HREE fractionations (La_N/Yb_N = 4.54–7.47). Such geochemical enrichment from N-MORB to E-MORB composition indicates mixing of melts derived from a depleted mantle and a fertile mantle/plume source at the spreading center. On the other hand, type-II mafic intrusives are hornblende bearing gabbros of SSZ-type tholeiitic composition with low Ti content (TiO₂ = 0.54 wt.%–0.86 wt.%) and depleted LREE pattern with respect to HREE (La_N/Yb_N = 0.37–0.49). They also have high Ba/Zr (1.13–2.82), Ba/Nb (45.56–151.66) and Ba/Th (84.58–744.19) and U/Th ratios (0.37–0.67) relative to the primitive mantle, which strongly represents the melt composition generated by partial melting of depleted lithospheric mantle wedge contaminated by hydrous fluids derived from subducting oceanic lithosphere in a forearc setting. Their subduction related origin is also supported by presence of calcium-rich plagioclase (An_16.6–32.3). Geothermometry calculation shows that the hornblende bearing (type-II) mafic rocks crystallized at temperature in range of 565°–625 °C ± 50 (at 10 kbar). Based on these available mineralogical and geochemical evidences, we conclude that mid ocean ridge (MOR) type mafic intrusive rocks from the NMO represent the section of older oceanic crust which was generated during the divergent process of the Indian plate from the Australian plate during Cretaceous period. Conversely, the hornblende-bearing gabbros (type-II) represent the younger oceanic crust which was formed at the forearc region by partial melting of the depleted mantle wedge slightly modified by the hydrous fluids released from the subducting oceanic slab during the initial stage of subduction of Indian plate beneath the Myanmar plate.     

Thermo-mechanical modeling of the obduction process based on the Oman Ophiolite case

Obduction emplaces regional-scale fragments of oceanic lithosphere (ophiolites) over continental lithosphere margins of much lower density. For this reason, the mechanisms responsible for obduction remain enigmatic in the framework of plate tectonics. We present two-dimensional (2D) thermo-mechanical models of obduction and investigate the possible dynamics and physical controls of this process. Model geometry and boundary conditions are based on available geological and geochronological data and numerical modeling results are validated against petrological and structural observations of the Oman (Semail) Ophiolite. Our model reproduces the stages of oceanic subduction initiation away from the Arabian margin, the emplacement of the Oman Ophiolite on top of it, and the domal exhumation of the metamorphosed margin through the ophiolitic nappe. A systematic study indicates that 350–400 km of bulk shortening provides the best fit for both maximum pressure–temperature conditions of the metamorphosed margin (1.5–2.5 GPa/450–600 °C) and the dimension of the ophiolitic nappe (~ 170 km width). Our results confirm that a thermal anomaly located close to the Arabian margin (~ 100 km) is needed to initiate obduction. We further suggest that a strong continental basement rheology is a prerequisite for ophiolite emplacement.     

Variable behavior of the Dead Sea Fault along the southern Arava segment from GPS measurements

The Dead Sea Fault is a major strike-slip fault bounding the Arabia plate and the Sinai subplate. On the basis of three GPS campaign measurements, 12 years apart, at 19 sites distributed in Israel and Jordan, complemented by Israeli permanent stations, we compute the present-day deformation across the Wadi Arava fault, the southern segment of the Dead Sea Fault. Elastic locked-fault modelling of fault-parallel velocities provides a slip rate of 4.7 ± 0.7 mm/yr and a locking depth of 11.6 ± 5.3 km in its central part. Along its northern part, south of the Dead Sea, the simple model proposed for the central profile does not fit the velocity field well. To fit the data, two faults have to be taken into account, on both sides of the sedimentary basin of the Dead Sea, each fault accommodating ∼ 2 mm/yr. Locking depths are small (less than 2 km on the western branch, ∼ 6 km on the eastern branch). Along the southern profile, we are once again unable to fit the data using the simple model, similar to the central profile. It is very difficult to propose a velocity greater than 4 mm/yr, i.e. smaller than that along the central profile. This leads us to propose that a part of the relative movement from Sinai to Arabia is accommodated along faults located west of our profiles.     

Pliocene volcano-tectonics and paleogeography of the Turkana Basin,Kenya and Ethiopia

The distribution of hominin fossil sites in the Turkana Basin, Kenya is intimately linked to the history of the Omo River, which affected the paleogeography and ecology of the basin since the dawn of the Pliocene. We report new geological data concerning the outlet channel of the Omo River between earliest Pliocene and final closure of the Turkana Basin drainage system in the latest Pliocene to earliest Quaternary. Throughout most of the Pliocene the Omo River entered the Turkana Basin from its source in the highlands of Ethiopia and exited the eastern margin of the basin to discharge into the Lamu embayment along the coast of the Indian Ocean. During the earliest Pliocene the river’s outlet was located in the northern part of the basin, where a remnant outlet channel is preserved in basalts that pre-date eruption of the Gombe flood basalt between 4.05 and 3.95 Ma. The outlet channel was faulted down to the west prior to 4.05 Ma, forming a natural dam behind which Lake Lonyumun developed. Lake Lonyumun was drained between 3.95 and 3.9 Ma when a new outlet channel formed north of Loiyangalani in the southeastern margin of the Turkana Basin. That outlet was blocked by Lenderit Basalt lava flows between 2.2 and 2.0 Ma. Faulting that initiated either during or shortly after eruption of the Lenderit Basalt closed the depression that is occupied by modern Lake Turkana to sediment and water.Several large shield volcanoes formed east of the Turkana Basin beginning by 2.5–3.0 Ma, volcanism overlapping in time, but probably migrating eastward from Mount Kulal on the eastern edge of the basin to Mount Marsabit located at the eastern edge of the Chalbi Desert. The mass of the volcanic rocks loaded and depressed the lithosphere, enhancing subsidence in a shallow southeast trending depression that overlay the Cretaceous and Paleogene (?) Anza Rift. Subsidence in this flexural depression guided the course of the Omo River towards the Indian Ocean, and also localized accumulations of lava along the margins of the shield volcanoes. Lava flows at Mount Marsabit extended across the Omo River Valley after 1.8–2.0 Ma based on estimated ages of fossils in lacustrine and terrestrial deposits, and possibly by as early as 2.5 ± 0.3 Ma based on dating of a lava flow. During the enhanced precipitation in latest Pleistocene and earliest Holocene (11–9.5 ka) this flexural depression became the site of Lake Chalbi, which was separated from Lake Turkana by a tectonically controlled drainage divide.     

Detrital zircon U–Pb geochronology,trace-element and Hf isotope geochemistry of the metasedimentary rocks in the Eastern Himalayan syntaxis: Tectonic and paleogeographic implications

The origin of the Greater Himalayan Sequence in the Himalaya and the paleogeographic position of the Lhasa terrane within Gondwanaland remain controversial. In the Eastern Himalayan syntaxis, the basement complexes of the northeastern Indian plate (Namche Barwa Complex) and the South Lhasa terrane (Nyingchi Complex) can be studied to explore these issues. Detrital zircons from the metasedimentary rocks in the Namche Barwa Complex and Nyingchi Complex yield similar U–Pb age spectra, with major age populations of 1.00–1.20 Ga, 1.30–1.45 Ga, 1.50–1.65 Ga and 1.70–1.80 Ga. The maximum depositional ages for their sedimentary protoliths are ~ 1.0 Ga based on the mean ages of the youngest three detrital zircons. Their minimum depositional ages are ~ 477 Ma for the Namche Barwa Complex and ~ 499 Ma for the Nyingchi Complex. Detrital zircons from the Namche Barwa Complex and Nyingchi Complex also display similar trace-element signatures and Hf isotopic composition, indicating that they were derived from common provenance. The trace-element signatures of 1.30–1.45 Ga detrital zircons indicate that the 1.3–1.5 Ga alkalic and mafic rocks belt in the southeastern India is a potential provenance. Most 1.50–1.65 Ga zircons have positive ε_Hf(t) values (+ 1.2 to + 9.0), and most 1.70–1.80 Ga zircons have negative ε_Hf(t) values (− 7.1 to − 1.9), which are compatible with those of the Paleo- to Mesoproterozoic orthogneisses in the Namche Barwa Complex. Provenance analysis indicates that the southern Indian Shield, South Lhasa terrane and probably Eastern Antarctica were the potential detrital sources. Combined with previous studies, our results suggest that: (1) the Namche Barwa Complex is the northeastern extension of the Greater Himalaya Sequence; (2) the metasedimentary rocks in the Namche Barwa Complex represent distal deposits of the northern Indian margin relative to the Lesser Himalaya; (3) the South Lhasa terrane was tectonically linked to northern India before the Cambrian.     

First tectonic-geomorphology study along the Longmu–Gozha Co fault system,Western Tibet

The Longmu–Gozha Co left-lateral strike-slip fault system (LGCF) is located in remote western Tibet, forming a triple junction with both the Altyn Tagh fault (ATF) and the Karakorum fault (KF), the two major strike-slip faults in the region. The Ashikule, Gozha Co and Longmu Co faults are clear and distinct left-stepping en-echelon faults, together forming the LGCF system. Although poorly documented, quantifying its activity remains a key problem to understand the kinematics and the tectonic history of the westernmost Tibetan Plateau. Indeed, the Karakax fault (NW segment of the ATF), LGCF and KF together control the tectonics of western Tibet which itself controls the extrusion of Tibet towards the east, with the LGCF acting as a natural boundary for eastward motion of the Tibetan Plateau due to India's northward impingement. The LGCF system shows clear and impressive morphological indications of left-lateral active shear, that we quantify using field measurements (terrestrial LIDAR) along with ¹⁰Be surface-exposure dating. Our data suggest a slip-rate < 3 mm/yr, consistent with geodetic and block model studies. While it is on the order of the Karakax fault slip-rate (~ 2 mm/yr), it is smaller than those along the ATF and KF (> 9 and > 8 mm/yr, respectively), yielding a few mm/yr of extension accommodated most likely in the Ashikule graben and surroundings, located between the ATF and Karakax faults. Numerous evidences of recent tectonic-related events are present in the vicinity, such as the 1951 volcanic eruption as well as the 2008 and 2014 Ms 7.3 Yutian earthquakes, attesting of its high activity. In addition, the LGCF's en-echelon geometry and identical direction with the ATF, as well as smaller geological offsets and lower slip-rate compared to those on the surrounding faults, suggest that this segment of the ATF may be the most recent.     

Arc magmatism in the Delhi Fold Belt: SHRIMP U–Pb zircon ages of granitoids and implications for Neoproterozoic convergent margin tectonics in NW India

The northwestern region of Peninsular India preserves important records of Precambrian plate tectonics and the role of Indian continent within Proterozoic supercontinents. In this study, we report precise SHRIMP zircon U–Pb ages from granitoids from the Sirohi terrane located along the western fringe of the Delhi Fold Belt in Rajasthan, NW India. The data reveal a range of Neoproterozoic ages from plagiogranite of Peshua, foliated granite of Devala, and porphyritic granite of Sai with zircon crystallization from magmas at 1015 ± 4.4 Ma, 966.5 ± 3.5 and 808 ± 3.1 respectively. The plagiogranite shows high SiO₂, Na₂O and extremely low K₂O, Rb, Ba, comparable with typical oceanic plagiogranites. These rocks possess low LREE and HREE concentrations and a relatively flat LREE–HREE slope, a well-developed negative Eu-anomaly and conspicuous Nb and Ti anomalies. Compared to the plagiogranite, the foliated Devala granite shows higher SiO₂ and moderate Na₂O, together with high K₂O and comparatively higher Rb, Ba, Sr and REE, with steep REE profiles and a weak positive Eu anomaly. In contrast to the plagiogranite and foliated granite, the porphrytic Sai granite has comparatively lower SiO₂ moderately higher Na₂O, extremely high Y, Zr, Nb and elevated REE. The geochemical features of the granitoids [HFSE depletion and LILE enrichment, Nb- and Ta-negative anomalies], and their plots in the fields of Volcanic Arc Granites and those from active continental margins in tectonic discrimination diagrams suggest widespread Neoproterozoic arc magmatism with changing magma chemistry in a protracted subduction realm. Our results offer important insights into a long-lived active continental margin in NW India during early and mid Neoproterozoic, consistent with recent similar observations on Cryogenian magmatic arcs widely distributed along the margins of the East African Orogen, and challenge some of the alternate models which link the magmatism to extensional tectonics associated with Rodinia supercontinent breakup.