Multiple garnet growth in garnet-kyanite-staurolite gneiss, Pangong metamorphic complex, Ladakh Himalaya: New constraints on tectonic setting

Garnet-kyanite-staurolite gneiss in the Pangong complex, Ladakh Himalaya, contains porphyroblastic euhedral garnets, blades of kyanite and resorbed staurolite surrounded by a fine-grained muscovite-biotite matrix associated with a leucogranite layer. Sillimanite is absent. The gneiss contains two generations of garnet in cores and rims that represent two stages of metamorphism. Garnet cores are extremely rich in Mn (X(Sps) = 0.35-038) and poor in Fe (X(Alm) = 0.40-0.45), whereas rims are relatively Mn-poor (X(Sps) =0.07-0.08), and rich in Fe (X(Alm), = 0.75-0.77). We suggest that garnet cores formed during prograde metamorphism in a subduction zone followed by abrupt exhumation, during early collision of the Ladakh arc and Karakoram block. The subsequent India-Asia continental collision subducted the metamorphic rocks to a mid-crustal level, where the garnet rims overgrew the Mn-rich cores at ca. 680 degrees C and ca. 8.5 kbar. PT calculations were estimated from phase diagrams calculated using a calculated bulk chemical composition in the Mn-NCKFMASHT system for the garnet-kyanite-staurolite-bearing assemblage. Muscovites from the metamorphic rocks and associated leucogranites have consistent K-Ar ages (ca. 10 Ma), closely related to activation of the Karakoram fault in the Pangong metamorphic complex. These ages indicate the contemporaneity of the exhumation of the metamorphic rocks and the cooling of the leucogranites. (C) 2011 Elsevier B.V. All rights reserved.

The 18 September 2011, North Sikkim earthquake

The 18 September 2011, magnitude Mw 6.9 earthquake close to the Nepal-Sikkim border caused significant damage due to ground shaking and caused several landslides. Observations from the post-earthquake surveys in the affected areas within Sikkim suggest that the poorly engineered, multistoried structures were relatively more impacted. Those located on alluvial terraces were also affected. The morphology of the region is prone to landslides and the possibility for their increased intensity during the forthcoming monsoon need to be considered seriously. From the seismotectonic perspective, the mid-crustal focal depth of the North Sikkim earthquake reflects the ongoing deformation of the subducting Indian plate.

Dynamics of the North American continent

The forces that cause deformation of western North America have been debated for decades. Recent studies, primarily based on analysis of crustal stresses in the western United States, have suggested that the deformation of the region is mainly controlled by gravitational potential energy (GPE) variations and boundary loads, with basal tractions due to mantle flow playing a relatively minor role. We address these issues by modelling the deviatoric stress field over western North America from a 3-D finite element mantle circulation model with lateral viscosity variations. Our approach takes into account the contribution from both topography and shallow lithosphere structure (GPE) as well as that from deeper mantle flow in one single model, as opposed to separate lithosphere and circulation models, as has been done so far. In addition to predicting the deviatoric stresses we also jointly fit the constraints of geoid, dynamic topography and plate motion both globally and over North America, in order to ensure that the forces that arise in our models are dynamically consistent. We examine the sensitivity of the dynamic models to different lateral viscosity variations. We find that circulation models that include upper mantle slabs yield a better fit to observed plate velocities. Our results indicate that a model of GPE variations coupled with mantle convection gives the best fit to the observational constraints. We argue that although GPE variations control a large part of the deformation of the western United States, deeper mantle tractions also play a significant role. The average deviatoric stress magnitudes in the western United States range 30-40 MPa. The cratonic region exhibits higher coupling to mantle flow than the rest of the continent. We find that a relatively strong San Andreas fault gives a better fit to the observational constraints, especially that of plate velocity in western North America.

Chromium isotope variations (delta Cr-53/52) in mantle-derived sources and their weathering products: Implications for environmental studies and the evolution of delta Cr-53/52 in the Earth's mantle over geologic time

Here we report chromium isotope compositions, expressed as delta Cr-53/ 52 in per mil (&) relative to NIST 979, measured in selected Cr-rich minerals and rocks formed by the primary magmatic as well as the secondary metamorphic and weathering processes. The main objectives of this study were: (i) to further constrain the isotope composition of the Earth's mantle Cr inventory and its possible variation during geological history, based on the analysis of globally distributed and stratigraphically constrained mantle-derived chromites; and (ii) to investigate the magnitude and systematics of Cr isotope fractionation during oxidative weathering and secondary alteration (i. e., hydration, serpentinization) of the magmatic Cr sources. Specifically, we analyzed delta Cr-53/ 52 in a set of globally distributed mantle-derived chromites (FeMgCr2O4, n = 30) collected from various locations in Europe, Asia, Africa and South America, and our results confirm that a chromite-hosted Earth's mantle Cr inventory is uniform at - 0.079 +/- 0.129& (2SD), which we named here as a ` canonical' mantle d 53/ 52 Cr signature. Furthermore our dataset of stratigraphically constrained chromites, whose crystallization ages cover most of the Earth's geological history, indicate that the bulk Cr isotope composition of the chromite-hosted mantle inventory has remained uniform, within about +/- 0.100&, since at least the Early Archean times (similar to 3500 million years ago, Ma). To investigate the systematics of Cr isotope fractionation associated with alteration processes we analyzed a number of secondary Cr-rich minerals and variably altered ultramafic rocks (i. e., serpentinized harzburgites, lherzolites) that revealed large positive delta Cr-53/ 52 anomalies that are systematically shifted to higher values with an increasing degree of alteration and serpentinization. The degree of aqueous alteration and serpentinization was quantified by the abundances of fluid-mobile (Rb, K) elements, and by the Loss On Ignition (LOI) parameter, which determines the amount of structurally bound water (OH/ H2O) present in secondary hydrated minerals like serpentine. Overall, we observed that altered ultramafic rocks that yielded the highest LOI values, and the lowest amounts of fluid mobile elements, also yielded the heaviest delta Cr-53/ 52 signatures. Therefore, we conclude that secondary alteration (i.e., hydration, serpentinization) of ultramafic rocks in near-surface oxidative environments tend to shift the bulk Cr isotope composition of the weathered products to isotopically heavier values, pointing to a dynamic redox cycling of Cr in the Earth's crustal and near-surface environments. Hence, if validated by future

Elastic thickness structure of the Andaman subduction zone: Implications for convergence of the Ninetyeast Ridge

We use the Bouguer coherence (Morlet isostatic response function) technique to compute the spatial variation of effective elastic thickness (T-e) of the Andaman subduction zone. The recovered T-e map resolves regional-scale features that correlate well with known surface structures of the subducting Indian plate and the overriding Burma plate. The major structure on the India plate, the Ninetyeast Ridge (NER), exhibits a weak mechanical strength, which is consistent with the expected signature of an oceanic ridge of hotspot origin. However, a markedly low strength (0< T-e <3 km) in that region, where the NER is close to the Andaman trench (north of 10 N), receives our main attention in this study. The subduction geometry derived from the Bouguer gravity forward modeling suggests that the NER has indented beneath the Andaman arc. We infer that the bending stresses of the viscous plate, which were reinforced within the subducting oceanic plate as a result of the partial subduction of the NER buoyant load, have reduced the lithospheric strength. The correlation, T-e < T-s (seismogenic thickness) reveals that the upper crust is actively deforming beneath the frontal arc Andaman region. The occurrence of normal-fault earthquakes in the frontal arc, low Te zone, is indicative of structural heterogeneities within the subducting plate. The fact that the NER along with its buoyant root is subducting under the Andaman region is inhibiting the subduction processes, as suggested by the changes in trench line, interrupted back-arc volcanism, variation in seismicity mechanism, slow subduction, etc. The low T-e and thinned crustal structure of the Andaman back-arc basin are attributed to a thermomechanically weakened lithosphere. The present study reveals that the ongoing back-arc spreading and strike-slip motion along the West Andaman Fault coupled with the ridge subduction exerts an important control on the frequency and magnitude of seismicity in the Andaman region. (C) 2013 Elsevier Ltd. All rights reserved.

Tectonic inheritance of the Indian Shield: New insights from its elastic thickness structure

A new evaluation of the elastic thickness (Te) structure of the Indian Shield, derived from isotropic fan wavelet methodology, documents spatial variations of lithospheric deformation in different tectonic provinces correlated with episodic tectono-thermal events. The Te variations corroborated by shear velocity, crustal thickness, and seismogenic thickness reveal the heterogeneous rheology of the Indian lithosphere. The thinned, attenuated lithosphere beneath Peninsular India is considered to be the reason for its mechanically weak strength (<30 km), where a decoupled crust-mantle rheology under different surface/subsurface loading structures may explain the prominent low Te patterns. The arcuate Te structure of the Western Dharwar province and a NNE-trending band of low Te anomaly in the Southern Granulite Terrane are intriguing patterns. The average Te values (40-50 km) of the Central Indian Tectonic Zone, the Bastar Craton, and the northern Eastern Ghats Mobile Belt are suggestive of old, stable, Indian lithosphere, which was not affected by any major tectono-thermal events after cratonic stabilization. We propose that the anomalously high Te (60-85 km) and high S-wave velocity zone to the north of the Narmada-Son Lineament, mainly in NW Himalaya, and the northern Aravalli and Bundelkhand Cratons, suggest that Archean lithosphere characterized by a high velocity mantle keel supports the orogenic topographic loads in/near the Himalaya. The Te map clearly segments the volcanic provinces of the Indian Shield, where the signatures of the Reunion, Marion, and Kerguelen hotspots are indicated by significantly low Te patterns that correlate with plume- and rift-related thermal and mechanical rejuvenation, magmatic underplating, and crustal necking. The correlations between Te variations and the occurrence of seismicity over seismically active zones reveal different causal relationships, which led to the current seismogenic zonation of the Indian Shield. (C) 2013 Elsevier B.V. All rights reserved.

India-Madagascar paleo-fit based on flexural isostasy of their rifted margins

The present study contributes new constraints on, and definitions of, the reconstructed plate margins of India and Madagascar based on flexural isostasy along the Western Continental Margin of India (WCMI) and the Eastern Continental Margin of Madagascar (ECMM). We have estimated the nature of isostasy and crustal geometry along the two margins, and have examined their possible conjugate structure. Here we utilize elastic thickness (Te) and Moho depth data as the primary basis for the correlation of these passive margins. We employ the flexure inversion technique that operates in spatial domain in order to estimate the spatial variation of effective elastic thickness. Gravity inversion and flexure inversion techniques are used to estimate the configuration of the Moho/Crust-Mantle Interface that reveals regional correlations with the elastic thickness variations. These results correlate well with the continental and oceanic segments of the Indian and African plates. The present study has found a linear zone of anomalously low-Te (1-5 km) along the WCMI (similar to 1680 km), which correlates well with the low-Te patterns obtained all along the ECMM. We suggest that the low-Te zones along the WCMI and ECMM represent paleo-rift inception points of lithosphere thermally and mechanically weakened by the combined effects of the Marion hotspot and lithospheric extension due to rifting. We have produced an India-Madagascar paleo-fit representing the initial phase of separation based on the Te estimates of the rifted conjugate margins, which is confirmed by a close-fit correlation of Moho geometry and bathymetry of the shelf margins. The matching of tectonic lineaments, lithologies and geochronological belts between India and Madagascar provide an additional support for the present plate reconstruction. (C) 2014 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

Neoarchean arc magmatism followed by high-temperature, high-pressure metamorphism in the Nilgiri Block, southern India

The Nilgiri Block, southern India is an exhumed lower crust formed through arc magmatic processes in the Neoarchean. The main lithologies in this terrane include charnockites, gneisses, volcanic tuff, metasediments, banded iron formation and mafic-ultramafic bodies. Mafic-ultramafic rocks are present towards the northern and central part of the Nilgiri Block. We examine the evolution of these mafic granulites/metagabbros by phase diagram modeling and U-Pb sensitive high resolution ion microprobe (SHRIMP) dating. They consist of a garnet-clinopyroxene-plagioclase-hornblende-ilmenite +/- orthopyroxene +/- rutile assemblage. Garnet and clinopyroxene form major constituents with labradorite and orthopyroxene as the main mineral inclusions. Labradorite, identified using Raman analysis, shows typical peaks at 508 cm(-1), 479 cm(-1), 287 cm(-1) and 177 cm(-1). It is stable along with orthopyroxene towards the low-pressure high-temperature region of the granulite fades (M1 stage). Subsequently, orthopyroxene reacted with plagioclase to form the peak garnet + clinopyroxene + rutile assemblage (M2 stage). The final stage is represented by amphibolite facies-hornblende and plagioclase-rim around the garnet-clinopyroxene assemblage (M3 stage). Phase diagram modeling shows that these mafic granulites followed an anticlockwise P-T-t path during their evolution. The initial high-temperature metamorphism (M1 stage) was at 850-900 degrees C and similar to 9 kbar followed by high-pressure granulite fades metamorphism (M2 stage) at 850-900 degrees C and 14-15 kbar. U-Pb isotope studies of zircons using SHRIMP revealed late Neoarchean to early paleoproterozoic ages of crystallization and metamorphism respectively. The age data shows that these mafic granulites have undergone arc magmatism at ca. 25392 +/- 3 Ma and high-temperature, high-pressure metamorphism at ca. 2458.9 +/- 8.6 Ma. Thus our results suggests a late Neoarchean arc magmatism followed by early paleoproterozoic high-temperature, high-pressure granulite facies metamorphism due to the crustal thickening and suturing of the Nilgiri Block onto the Dharwar Craton. (C) 2015 Elsevier B.V. All rights reserved.

Evolution of high-pressure mafic granulites and pelitic gneisses from NE Madagascar: Tectonic implications

The occurrence of high-pressure mafic-ultramafic bodies within major shear zones is one of the indicators of paleo-subduction. In mafic granulites of the Andriamena complex (north-eastern Madagascar) we document unusual textures including garnet-clinopyroxene-quartz coronas that formed after the breakdown of orthopyroxene-plagioclase-ilmenite. Textural evidence and isochemical phase diagram calculations in the Na2O-CaO-K2O-FeO-MgO-Al2O3-SiO2-H2O-TiO2 system indicate a pressure-temperature (P-T) evolution from an isothermal (780 degrees C) pressure up to c. 24 kbar to decompression and cooling. Such a P-T trajectory is typically attained in a subduction zone setting where a gabbroic/ultramafic complex is subducted and later exhumed to the present crustal level during oceanic closure and final continental collision. The present results suggest that the presence of such deeply subducted rocks of the Andriamena complex is related to formation of the Betsimisaraka suture. LA-ICPMS U-Pb zircon dating of pelitic gneisses from the Betsimisaraka suture yields low Th/U ratios and protolith ages ranging from 2535 to 2625 Ma. A granitic gneiss from the Alaotra complex yields a zircon crystallization age of ca. 818 Ma and Th/U ratios vary from 1.08 to 2.09. K-Ar dating of muscovite and biotite from biotite-kyanite-sillimanite gneiss and garnet-biotite gneiss yields age of 486 +/- 9 Ma and 459 +/- 9 Ma respectively. We have estimated regional crustal thicknesses in NE Madagascar using a flexural inversion technique, which indicates the presence of an anomalously thick crust (c. 43 km) beneath the Antananarivo block. This result is consistent with the present concept that subduction beneath the Antananarivo block resulted in a more competent and thicker crust. The textural data, thermodynamic model, and geophysical evidence together provide a new insight to the subduction history, crustal thickening and evolution of the high-pressure Andriamena complex and its link to the terminal formation of the Betsimisaraka suture in north-eastern Madagascar. (C) 2015 Elsevier B.V. All rights reserved.