High P-T granulite relicts from the Imjingang belt, South Korea:Tectonic significance

This study presents a detailed description on crustal metamorphic signatures of garnet-clinopyroxene-quartz-rutile-bearing high P-T granulites, Samgot unit, Imajingang belt, northwestern Korean Peninsula that formed during Permo-Triassic regional metamorphism related to the amalgamation of East Asian continental fragments. Lenses and blocks of high P-T granulites and garnet-bearing leucosomes occur within mafic metamorphic rocks (mainly amphibolites). The mafic blocks comprise relicts of granoblastic garnet and clinopyroxene with medium-grained quartz and rutile. These relict mineral assemblages are confined to local micro-domains and constitute remnants of peak metamorphism. Plagioclase and amphibole form only as retrograde phases in medium ton coarse-grained moats that rim grain boundaries between relict peak mineral assemblages. This microstructure represents the reaction between garnet, clinopyroxene, quartz and rutile in the presence of melt to form amphibole, plagioclase and titanite with minor biotite. The leucosome domains consist of euhedral garnets within the quartz-K feldspar-plagioclase (granitic) matrix, probably representing peritectic garnet growth along with melting. The rare earth element (REE) composition of minerals also support the peritectic garnet growth with a positive Eu/Eu* (positive Eu anomaly), while the relict garnet shows a slight negative anomaly typical for high-grade granulites. The peak-metamorphic conditions calculated from thermodynamic modeling and compositional isopleths indicate a temperature around c. 900 degrees C at a pressure around c. 20 kbar. The present P-T path indicates a clear multi-stage decompression history with initial decompression and cooling followed by a stage of decompression during hydration possibly during Late Triassic exhumation. The results from this study together with the presence of eclogites from the Hongsung area suggest that the Imjingang area and the western Gyeonggi massif likely resided at crustal levels deeper than those of the eastern and southern part of the Gyeonggi massif. (C) 2009 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

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.

Spatial variations of effective elastic thickness over the Ninetyeast Ridge and implications for its structure and tectonic evolution

We present new data on the strength of oceanic lithosphere along the Ninetyeast Ridge (NER) from two independent methods: spectral analysis (Bouguer coherence) using the fan wavelet transform technique, and spatial analysis (flexure inversion) with the convolution method. The two methods provide effective elastic thickness (T-e) patterns that broadly complement each other, and correlate well with known surface structures and regional-scale features. Furthermore, our study presents a new high resolution database on the Moho configuration, which obeys flexural isostasy, and exhibit regional correlations with the T-e variations. A continuous ridge structure with a much lower T-e value than that of normal oceanic lithosphere provides strong support for the hotspot theory. The derived T-e values vary over the northern (higher T-e similar to 10-20 km), central (anomalously low T-e similar to 0-5 km), and southern (low T-e similar to 5 km) segments of the NER. The lack of correlation of the T-e value with the progressive aging of the lithosphere implies differences in thermo-mechanical setting of the crust and underlying mantle in different parts of the NER, again indicating diversity in their evolution. The anomalously low T-e and deeper Moho (similar to 22 km) estimates of the central NER (between 0.5 degrees N and 17 degrees S) are attributed to the interaction of a hotspot with the Wharton spreading ridge that caused significant thermal rejuvenation and hence weakening of the lithosphere. The higher mechanical strength values in the northern NER (north of 0.5 degrees N) may support the idea of off-ridge emplacement and a relatively large plate motion at the time of volcanism. The low T-e and deeper Moho (similar to 22 km) estimates in the southern part (south of 17 degrees S) suggest that the lithosphere was weak and therefore younger at the time of volcanism, and this supports the idea that the southern NER was emplaced on the edge of the Indian plate. (C) 2013 Elsevier B.V. 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.

Geomorphology reveals active decollement geometry in the central Himalayan seismic gap

The similar to 700-km-long ``central seismic gap'' is the most prominent segment of the Himalayan front not to have ruptured in a major earthquake during the last 200-500 yr. This prolonged seismic quiescence has led to the proposition that this region, with a population >10 million, is overdue for a great earthquake. Despite the region's recognized seismic risk, the geometry of faults likely to host large earthquakes remains poorly understood. Here, we place new constraints on the spatial distribution of rock uplift within the western similar to 400 km of the central seismic gap using topographic and river profile analyses together with basinwide erosion rate estimates from cosmogenic Be-10. The data sets show a distinctive physiographic transition at the base of the high Himalaya in the state of Uttarakhand, India, characterized by abrupt strike-normal increases in channel steepness and a tenfold increase in erosion rates. When combined with previously published geophysical imaging and seismicity data sets, we interpret the observed spatial distribution of erosion rates and channel steepness to reflect the landscape response to spatially variable rock uplift due to a structurally coherent ramp-flat system of the Main Himalayan Thrust. Although it remains unresolved whether the kinematics of the Main Himalayan Thrust ramp involve an emergent fault or duplex, the landscape and erosion rate patterns suggest that the decollement beneath the state of Uttarakhand provides a sufficiently large and coherent fault segment capable of hosting a great earthquake.

Coupling geomorphology parameters with lithology for micro-level groundwater resource assessment: a case study from semi-arid hard rock terrain in Tamil Nadu, India

The standard procedure of groundwater resource estimation in India till date is based on the specific yield parameters of each rock type (lithology) derived through pumping test analysis. Using the change in groundwater level, specific yield, and area of influence, groundwater storage change could be estimated. However, terrain conditions in the form of geomorphological variations have an important bearing on the net groundwater recharge. In this study, an attempt was made to use both lithology and geomorphology as input variables to estimate the recharge from different sources in each lithology unit influenced by the geomorphic conditions (lith-geom), season wise separately. The study provided a methodological approach for an evaluation of groundwater in a semi-arid hard rock terrain in Tirunelveli, Tamil Nadu, India. While characterizing the gneissic rock, it was found that the geomorphologic variations in the gneissic rock due to weathering and deposition behaved differently with respect to aquifer recharge. The three different geomorphic units identified in gneissic rock (pediplain shallow weathered (PPS), pediplain moderate weathered (PPM), and buried pediplain moderate (BPM)) showed a significant variation in recharge conditions among themselves. It was found from the study that Peninsular gneiss gives a net recharge value of 0.13 m/year/unit area when considered as a single unit w.r.t. lithology, whereas the same area considered with lith-geom classes gives recharge values between 0.1 and 0.41 m/year presenting a different assessment. It is also found from this study that the stage of development (SOD) for each lith-geom unit in Peninsular gneiss varies from 168 to 230 %, whereas the SOD is 223 % for the lithology as a single unit.

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.

Evidence for Permo-Triassic collision in Far East Asia: The Korean collisional orogen

Eclogites from paragneiss in the Korean Peninsula are characterized by a peak pressure assemblage of garnet + omphacite + quartz + rutile, that is overprinted by multiphase symplectites involving augite, amphibole, orthopyroxene, ilmenite and plagioclase and by a similar high-pressure assemblage with a pronounced absence of the omphacite component in clinopyroxene formed during the peak and orthopyroxene in the retrograde stage. Eclogites were metamorphosed at a minimum pressures of not, vert, similar 20–23 kbar at temperatures of not, vert, similar 840–1000 °C, equivalent to a crustal depth of not, vert, similar 70–75 km, whereas high-pressure granulite in Late Paleozoic rocks underwent metamorphic conditions of not, vert, similar 18–19 kbar at not, vert, similar 950 °C with a minimum crustal depth of not, vert, similar 60–65 km. The presence of the eclogites and high-pressure granulite suggests deep-seated subduction of crustal complexes with metamorphism at different crustal levels. The eclogites were exhumed quickly resulting in near- isothermal decompression. On the other hand, the multistage exhumation of the high-pressure granulites suggests retrograde overprinting after initial decompression. The similarity of these petrological characteristics, metamorphic conditions and also the regional structural styles with those of the Sulu belt (China) strongly suggests the existence of a Permo-Triassic Alpine-type “Korean collision belt” in Far East Asia. This model provides a better understanding of the paleogeograpic evolution of Permo-Triassic East Asia, including a robust tectonic correlation of the Korean collision belt with the Qinling–Dabie–Sulu collision belt.

Counterclockwise exhumation of a hot orogen: The Paleoproterozoic ultrahigh-temperature granulites in the North China Craton

Regional metamorphic belts provide important constraints on the plate tectonic architecture of orogens. We report here a detailed petrologic examination of the sapphirine-bearing ultra-high temperature (UHT) granulites from the Jining Complex within the Khondalite Belt of the North China Craton (NCC). These granulites carry diagnostic UHT assemblages and their microstructures provide robust evidence to trace the prograde, peak and retrograde metamorphic evolution. The P–T conditions of the granulites estimated from XMgGrt(Mg/Fe + Mg) − XMgSpr isopleth calculations indicate temperature above 970 °C and pressures close to 7 kbar. We present phase diagrams based on thermodynamic computations to evaluate the mineral assemblages and microstructures and trace the metamorphic trajectory of the rocks. The evolution from Spl–Qtz–Ilm–Crd–Grt–Sil to Spr–Qtz–Crd–Opx–Ilm marks the prograde stage. The Spl–Qtz assemblage appears on the low-pressure side of the P–T space with Spr–Qtz stable at the high-pressure side, possibly representing an increase in pressure corresponding to compression. The spectacular development of sapphirine rims around spinel enclosed in quartz supports this inference. An evaluation of the key UHT assemblages based on model proportion calculation suggests a counterclockwise P–T path. With few exceptions, granulite-facies rocks developed along collisional metamorphic zones have generally been characterized by clockwise exhumation trajectories. Recent evaluation of the P–T paths of metamorphic rocks developed within collisional orogens indicates that in many cases the exhumation trajectories follow the model subduction geotherm, in accordance with a tectonic model in which the metamorphic rocks are subducted and exhumed along a plate boundary. The timing of UHT metamorphism in the NCC (c. 1.92 Ga) coincides with the assembly of the NCC within the Paleoproterozoic Columbia supercontinent, a process that would have involved subduction of passive margins sediments and closure of the intervening ocean. Thus, the counterclockwise P–T path obtained in this study correlates well with a tectonic model involving subduction and final collisional suturing, with the UHT granulites representing the core of the hot or ultra-hot orogen developed during Columbia amalgamation.

Geometry of sand-bed channels with seepage

Literature suggests that apart from a significant percentage of water loss from seepage, downward seepage from the channel causes an increase in the mobility of channel bed materials and thus changes the channel stability. Consequently, regime conditions (which provide the relationship among hydraulic parameters) should also be revised by incorporating downward seepage as an additional parameter. In the present work, regime conditions for the prediction of channel geometry in alluvial channels affected by downward seepage have been formulated on the basis of experimental observations. (C) 2011 Elsevier B.V. All rights reserved.

Profiling of late Tertiary-early Quaternary surface in the lower reaches of Narmada valley using microtremors

In this paper, we propose the first approximation for thickness of Quaternary sediment and late Quaternary early Tertiary topography for the part of lower reaches of Narmada valley in a systematic way using the shallow seismic method, that records both horizontal and vertical components of the microtremor (ambient noise) caused by natural processes. The measurements of microtremors were carried out at 31 sites spaced at a grid interval of 5 km s using Lennartz seismometer (5 s period) and City shark-II data acquisition system. The signals recorded were analysed for horizontal to the vertical (H/V) spectral ratio using GEOPSY software. For the present study, we concentrate on frequency range between 0.2 Hz and 10 Hz. The thickness of unconsolidated sediments at various sites is calculated based on non-linear regression equations proposed by Ibs-von Seht and Wohlenberg (1999) and Parolai et al. (2002). The estimated thickness is used to plot digital elevation model and cross profiles correlating with geomorphology and geology of the study area. (C) 2011 Elsevier Ltd. All rights reserved.

Deterministic seismic hazard macrozonation of India

Earthquakes are known to have occurred in Indian subcontinent from ancient times. This paper presents the results of seismic hazard analysis of India (6 degrees-38 degrees N and 68 degrees-98 degrees E) based on the deterministic approach using latest seismicity data (up to 2010). The hazard analysis was done using two different source models (linear sources and point sources) and 12 well recognized attenuation relations considering varied tectonic provinces in the region. The earthquake data obtained from different sources were homogenized and declustered and a total of 27,146 earthquakes of moment magnitude 4 and above were listed in the study area. The sesismotectonic map of the study area was prepared by considering the faults, lineaments and the shear zones which are associated with earthquakes of magnitude 4 and above. A new program was developed in MATLAB for smoothing of the point sources. For assessing the seismic hazard, the study area was divided into small grids of size 0.1 degrees x 0.1 degrees (approximately 10 x 10 km), and the hazard parameters were calculated at the center of each of these grid cells by considering all the seismic sources within a radius of 300 to 400 km. Rock level peak horizontal acceleration (PHA) and spectral accelerations for periods 0.1 and 1 s have been calculated for all the grid points with a deterministic approach using a code written in MATLAB. Epistemic uncertainty in hazard definition has been tackled within a logic-tree framework considering two types of sources and three attenuation models for each grid point. The hazard evaluation without logic tree approach also has been done for comparison of the results. The contour maps showing the spatial variation of hazard values are presented in the paper.

The hazard potential of the western segment of the Makran subduction zone, northern Arabian Sea

Evaluating the hazard potential of the Makran subduction zone requires understanding the previous records of the large earthquakes and tsunamis. We address this problem by searching for earthquake and tectonic proxies along the Makran Coast and linking those observations with the available constraints on historical seismicity and the tell-tale characteristics of sea floor morphology. The earthquake of Mw 8.1 of 1945 and the consequent tsunami that originated on the eastern part of the Makran are the only historically known hazardous events in this region. The seismic status of the western part of the subduction zone outside the rupture area of the 1945 earthquake remains an enigma. The near-shore shallow stratigraphy of the central part of Makran near Chabahar shows evidence of seismically induced liquefaction that we attribute to the distant effects of the 1945 earthquake. The coastal sites further westward around Jask are remarkable for the absence of liquefaction features, at least at the shallow level. Although a negative evidence, this possibly implies that the western part of Makran Coast region may not have been impacted by near-field large earthquakes in the recent past-a fact also supported by the analysis of historical data. On the other hand, the elevated marine terraces on the western Makran and their uplift rates are indicative of comparable degree of long-term tectonic activity, at least around Chabahar. The offshore data suggest occurrences of recently active submarine slumps on the eastern part of the Makran, reflective of shaking events, owing to the great 1945 earthquake. The ocean floor morphologic features on the western segment, on the contrary, are much subdued and the prograding delta lobes on the shelf edge also remain intact. The coast on the western Makran, in general, shows indications of progradation and uplift. The various lines of evidence thus suggest that although the western segment is potentially seismogenic, large earthquakes have not occurred there in the recent past, at least during the last 600 years. The recurrence period of earthquakes may range up to 1,000 years or more, an assessment based on the age of the youngest dated coastal ridge. The long elapsed time points to the fact that the western segment may have accumulated sufficient slip to produce a major earthquake.