Seismotectonics of the April-May 2015 Nepal earthquakes: An assessment based on the aftershock patterns, surface effects and deformational characteristics

Occurrence of the April 25, 2015 (Mw 7.8) earthquake near Gorkha, central Nepal, and another one that followed on May 12 (Mw 7.3), located similar to 140 km to its east, provides an exceptional opportunity to understand some new facets of Himalayan earthquakes. Here we attempt to assess the seismotectonics of these earthquakes based on the deformational field generated by these events, along with the spatial and temporal characteristics of their aftershocks. When integrated with some of the post-earthquake field observations, including the localization of damage and surface deformation, it became obvious that although the mainshock slip was mostly limited to the Main Himalayan Thrust (MHT), the rupture did not propagate to the Main Frontal Thrust (MFT). Field evidence, supported by the available InSAR imagery of the deformation field, suggests that a component of slip could have emerged through a previously identified out-of-sequence thrust/active thrust in the region that parallels the Main Central Thrust (MCT), known in the literature as a co-linear physiographic transitional zone called PT2. Termination of the first rupture, triggering of the second large earthquake, and distribution of aftershocks are also spatially constrained by the eastern extremity of PT2. Mechanism of the 2015 sequence demonstrates that the out-of-sequence thrusts may accommodate part of the slip, an aspect that needs to be considered in the current understanding of the mechanism of earthquakes originating on the MHT. (c) 2015 Elsevier Ltd. All rights reserved.

On the control of rapidly rotating convection by an axially varying magnetic field

The magnetic field in rapidly rotating dynamos is spatially inhomogeneous. The axial variation of the magnetic field is of particular importance because tall columnar vortices aligned with the rotation axis form at the onset of convection. The classical picture of magnetoconvection with constant or axially varying magnetic fields is that the Rayleigh number and wavenumber at onset decrease appreciably from their non-magnetic values. Nonlinear dynamo simulations show that the axial lengthscale of the self-generated azimuthal magnetic field becomes progressively smaller as we move towards a rapidly rotating regime. With a small-scale field, however, the magnetic control of convection is different from that in previous studies with a uniform or large-scale field. This study looks at the competing viscous and magnetic mode instabilities when the Ekman number E (ratio of viscous to Coriolis forces) is small. As the applied magnetic field strength (measured by the Elsasser number Lambda) increases, the critical Rayleigh number for onset of convection initially increases in a viscous branch, reaches an apex where both viscous and magnetic instabilities co-exist, and then falls in the magnetic branch. The magnetic mode of onset is notable for its dramatic suppression of convection in the bulk of the fluid layer where the field is weak. The viscous-magnetic mode transition occurs at Lambda similar to 1, which implies that small-scale convection can exist at field strengths higher than previously thought. In spherical shell dynamos with basal heating, convection near the tangent cylinder is likely to be in the magnetic mode. The wavenumber of convection is only slightly reduced by the self-generated magnetic field at Lambda similar to 1, in agreement with previous planetary dynamo models. The back reaction of the magnetic field on the flow is, however, visible in the difference in kinetic helicity between cyclonic and anticyclonic vortices.

Isotopic and geochemical characterization of invader tilapia fishes from water bodies of West Bengal and Karnataka, India

The otoliths (N = 12) of freshwater invasive species tilapia (Tilapia mossambicus) collected from two water bodies located at Kolkata and Bangalore, India, were analyzed for stable isotopes (delta 18O, delta 14C) and major and trace elements in order to assess the suitability of using otoliths as a tracer of aquatic environmental changes. The stable isotope analysis was done using the dual inlet system of a Finnigan-MAT 253 isotope ratio mass spectrometer (Thermo-Fisher, Bremen, Germany). Concentrations of major and trace elements were determined using a Thermo X-Series II quadrupole mass spectrometer. The stable isotope composition in tilapia otolith samples from Bangalore and Kolkata water bodies are quite good agreeing with that of the respective lake/pond and rain water. Elemental composition revealed in a pattern of Ca > Fe > Na > Sr > K > Ba > Cr > Mg > As > Mn > Zn > Co > Cu > Cd > Pb. The otoliths from Kolkata pond water are more enriched in Ba, Zn, Pb, Mn, Se, Cu, Zn, Cd, and Ni whereas Cr and As were found to be higher in otolith samples from Bangalore lake. The enrichment factor (EF) values of Cr were higher for both the sampling location in comparison with other metals, although all the studied metals exhibited EF values >1. The PCA shows clustering of metals in the otolith which are related either with the metabolic and physiological attributes or waterborne source. The study demonstrated the potential of stable isotope techniques to distinguish otolith specimens from varied climatic zone, while elemental composition recorded the quality of water at both the locations. The role of climate driving the quality of water can be understood by detailed and continuous monitoring of otolith specimens in the future. Future method allows reconstruction of climate and water quality from old specimens from field exposures or museum collection.

Onset of plane layer magnetoconvection at low Ekman number

We study the onset of magnetoconvection between two infinite horizontal planes subject to a vertical magnetic field aligned with background rotation. In order to gain insight into the convection taking place in the Earth's tangent cylinder, we target regimes of asymptotically strong rotation. The critical Rayleigh number Ra-c and critical wavenumber a(c) are computed numerically by solving the linear stability problem in a systematic way, with either stress-free or no-slip kinematic boundary conditions. A parametric study is conducted, varying the Ekman number E (ratio of viscous to Coriolis forces) and the Elsasser number. (ratio of the Lorentz force to the Coriolis force). E is varied from 10(-9) to 10(-2) and. from 10(-3) to 1. For a wide range of thermal and magnetic Prandtl numbers, our results verify and confirm previous experimental and theoretical results showing the existence of two distinct unstable modes at low values of E-one being controlled by the magnetic field, the other being controlled by viscosity (often called the viscous mode). It is shown that oscillatory onset does not occur in the range of parameters we are interested in. Asymptotic scalings for the onset of these modes are numerically confirmed and their domain of validity is precisely quantified. We show that with no-slip boundary conditions, the asymptotic behavior is reached for E < 10(-6) and establish a map in the (E, Lambda) plane. We distinguish regions where convection sets in either through the magnetic mode or through the viscous mode. Our analysis gives the regime in which the transition between magnetic and viscous modes may be observed. We also show that within the asymptotic regime, the role played by the kinematic boundary conditions is minimal. (C) 2015 AIP Publishing LLC.

Interlaboratory comparison of magnesium isotopic compositions of 12 felsic to ultramafic igneous rock standards analyzed by MC-ICPMS

To evaluate the interlaboratory mass bias for high-precision stable Mg isotopic analysis of natural materials, a suite of silicate standards ranging in composition from felsic to ultramafic were analyzed in five laboratories by using three types of multicollector inductively coupled plasma mass spectrometer (MC-ICPMS). Magnesium isotopic compositions from all labs are in agreement for most rocks within quoted uncertainties but are significantly (up to 0.3 parts per thousand in Mg-26/Mg-24, > 4 times of uncertainties) different for some mafic samples. The interlaboratory mass bias does not correlate with matrix element/Mg ratios, and the mechanism for producing it is uncertain but very likely arises from column chemistry. Our results suggest that standards with different matrices are needed to calibrate the efficiency of column chemistry and caution should be taken when dealing with samples with complicated matrices. Well-calibrated standards with matrix elements matching samples should be used to reduce the interlaboratory mass bias.

Linking monazite geochronology with fluid infiltration and metamorphic histories: Nature and experiment

Migmatised metapelites from the Kodaikanal region, central Madurai Block, southern India have undergone ultrahigh-temperature metamorphism (950-1000 degrees C; 7-8 kbar). In-situ electron microprobe Th-U-Pb isochron (CHIME) dating of monazites in a leucosome and surrounding silica-saturated and silica-poor restites from the same outcrop indicates three principal ages that can be linked to the evolutionary history of these rocks. Monazite grains from the silica-saturated restite have well-defined, inherited cores with thick rims that yield an age of ca. 1684 Ma. This either dates the metamorphism of the original metapelite or is a detrital age of inherited monazite. Monazite grains from the silica-poor restite, thick rims from the silica-saturated restite, and monazite cores from the leucosome have ages ranging from 520 to 540 Ma suggesting a mean age of 530 Ma within the error bars. In the leucosome the altered rim of the monazite gives an age of ca. 502 Ma. Alteration takes the form of Th-depleted lobes of monazite with sharp curvilinear boundaries extending from the monazite grain rim into the core. We have replicated experimentally these altered rims in a monazite-leucosome experiment at 800 degrees C and 2 kbar. This experiment, coupled with earlier published monazite-fluid experiments involving high pH alkali-bearing fluids at high P-T, helps to confirm the idea that alkali-bearing fluids, in the melt and along grain boundaries during crystallization, were responsible for the formation of the altered monazite grain rims via the process of coupled dissolution-reprecipitation. (C) 2015 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.

Fingerprinting environmental conditions and related stress using stable isotopic composition of rice (Oryza sativa L.) grain organic matter

Sea level rise (SLR) is a primary factor responsible for inundation of low-lying coastal regions across the world, which in turn governs the agricultural productivity. In this study, rice (Oryza sativa L.) cultivated seasonally in the Kuttanad Wetland, a SLR prone region on the southwest coast of India, were analysed for oxygen, hydrogen and carbon isotopic ratios (delta O-18, delta H-2 and delta C-13) to distinguish the seasonal environmental conditions prevalent during rice cultivation. The region receives high rainfall during the wet season which promotes large supply of fresh water to the local water bodies via the rivers. In contrast, during the dry season reduced river discharge favours sea water incursion which adversely affects the rice cultivation. The water for rice cultivation is derived from regional water bodies that are characterised by seasonal salinity variation which co-varies with the delta O-18 and delta H-2 values. Rice cultivated during the wet and the dry season bears the isotopic imprints of this water. We explored the utility of a mechanistic model to quantify the contribution of two prominent factors, namely relative humidity and source water composition in governing the seasonal variation in oxygen isotopic composition of rice grain OM. delta C-13 values of rice grain OM were used to deduce the stress level by estimating the intrinsic water use efficiency (WUEi) of the crop during the two seasons. 1.3 times higher WUE, was exhibited by the same genotype during the dry season. The approach can be extended to other low lying coastal agro-ecosystems to infer the growth conditions of cultivated crops and can further be utilised for retrieving paleo-environmental information from well preserved archaeological plant remains. (c) 2015 Elsevier Ltd. All rights reserved.

Dynamos driven by weak thermal convection and heterogeneous outer boundary heat flux

We use numerical dynamo models with heterogeneous core-mantle boundary (CMB) heat flux to show that lower mantle lateral thermal variability may help support a dynamo under weak thermal convection. In our reference models with homogeneous CMB heat flux, convection is either marginally supercritical or absent, always below the threshold for dynamo onset. We find that lateral CMB heat flux variations organize the flow in the core into patterns that favour the growth of an early magnetic field. Heat flux patterns symmetric about the equator produce non-reversing magnetic fields, whereas anti-symmetric patterns produce polarity reversals. Our results may explain the existence of the geodynamo prior to inner core nucleation under a tight energy budget. Furthermore, in order to sustain a strong geomagnetic field, the lower mantle thermal distribution was likely dominantly symmetric about the equator. (C) 2015 Elsevier B.V. All rights reserved.

Isotopic disequilibrium in Globigerina bulloides and carbon isotope response to productivity increase in Southern Ocean

Oxygen and carbon isotope ratios in planktonic foraminifera Globigerina bulloides collected from tow samples along a transect from the equatorial Indian ocean to the Southern Ocean (45 degrees E and 80 degrees E and 10 degrees N to 53 degrees S) were analysed and compared with the equilibrium delta O-18 and delta C-13 values of calcite calculated using the temperature and isotopic composition of the water column. The results agree within similar to 0.25% for the region between 10 degrees N and 40 degrees S and 75-200 m water depth which is considered to be the habitat of Globigerina bulloides. Further south (from 40 degrees S to 55 degrees S), however, the measured delta O-18 and delta C-13 values are higher than the expected values by similar to 2% and similar to 1% respectively. These enrichments can be attributed to either a `vital effect' or a higher calcification rate. An interesting pattern of increase in the delta C-13(DIC) value of the surface water with latitude is observed between 35 degrees S and similar to 60 degrees S, with a peak at similar to 42 degrees S. This can be caused by increased organic matter production and associated removal. A simple model accounting for the increase in the delta C-13(DIC) values is proposed which fits well with the observed chlorophyll abundance as a function of latitude.

Influence of the upwelling events on the delta C-13 and delta O-18 of the benthic bivalve shells of the South Western Continental Margin of India

In this study, the influence of the spatial and temporal variability of upwelling intensity and the associated biological productivity observed during different phases of summer monsoon along the southwestern continental margin of India (SWCMI) on the delta C-13 and delta O-18 of the inorganic biogenic carbonate shells was investigated. Multispecies benthic bivalve shells (1-5 mm) separated from ten surface sediment samples of SWCMI (off 12 degrees N, 10 degrees N and 9 degrees N) collected during the onset (OSM) and peak (PSM) phase of the summer monsoon of 2009 were analysed for delta C-13 and delta O-18. Sea surface temperature along the study region indicates prominent upwelling in PSM than in OSM. A comparison of analytical and predicted values for delta O-18 in the bivalve shells confirmed their in situ origin during both the sampling periods. During PSM, the delta C-13 values in the benthic bivalve shells were more depleted in C-13 than during OSM which recorded lower values of delta C-13 in dissolved inorganic carbon of bottom waters expected in the study region in PSM due to the upwelled waters, high surface productivity and the associated high degradation of the organic matter in the subsurface and bottom waters. However, this depletion of delta C-13 was not observed in benthic bivalve shells obtained from 10 degrees N, since it is influenced by high export fluxes of carbon from the Cochin estuary since early monsoon months.

Tracking the migration of the Indian continent using the carbonate clumped isotope technique on Phanerozoic soil carbonates

Approximately 140 million years ago, the Indian plate separated from Gondwana and migrated by almost 90 degrees latitude to its current location, forming the Himalayan-Tibetan system. Large discrepancies exist in the rate of migration of Indian plate during Phanerozoic. Here we describe a new approach to paleo-latitudinal reconstruction based on simultaneous determination of carbonate formation temperature and delta O-18 of soil carbonates, constrained by the abundances of C-13-O-18 bonds in palaeosol carbonates. Assuming that the palaeosol carbonates have a strong relationship with the composition of the meteoric water, delta O-18 carbonate of palaeosol can constrain paleo-latitudinal position. Weighted mean annual rainfall delta O-18 water values measured at several stations across the southern latitudes are used to derive a polynomial equation: delta(18)Ow = -0.006 x (LAT)(2) - 0.294 x (LAT) - 5.29 which is used for latitudinal reconstruction. We use this approach to show the northward migration of the Indian plate from 46.8 +/- 5.8 degrees S during the Permian (269 M. y.) to 30 +/- 11 degrees S during the Triassic (248 M. y.), 14.7 +/- 8.7 degrees S during the early Cretaceous (135 M. y.), and 28 +/- 8.8 degrees S during the late Cretaceous ( 68 M. y.). Soil carbonate delta O-18 provides an alternative method for tracing the latitudinal position of Indian plate in the past and the estimates are consistent with the paleo-magnetic records which document the position of Indian plate prior to 135 +/- 3 M. y.

Archean tectonics and crustal evolution of the Biligiri Rangan Block, southern India

The Southern Granulite Terrain in India is a collage of crustal blocks ranging in age from Archean to Neoproterozoic. This study investigate the tectonic evolution of one of the northernmost block- the Biligiri Block (BRB) through a multidisciplinary approach involving field investigation, petrographic studies, LA-ICPMS zircon U-Pb geochronology, Hf isotopic analyses, metamorphic P-T phase diagram computations, and crustal thickness modeling. The garnet bearing quartzofeldspathic gneiss from the central BRB preserve Mesoarchean magmatic zircons with ages between 3207 and 2806 Ma and positive epsilon Hf value (+2.7) which possibly indicates vestiges of a Mesoarchean primitive continental crust. The occurrence of quartzite-iron formation intercalation as well as ultramafic lenses along the western boundary of the BRB is interpreted to indicate that the Kollegal structural lineament is a possible paleo-suture. Phase diagram computation of a metagabbro from the southwestern periphery of the Kollegal suture zone reveals high-pressure (similar to 18.5 kbar) and medium-temperature (similar to 840 degrees C) metamorphism, likely during eastward subduction of the Western Dharwar oceanic crust beneath the Mesoarchean BRB. In the model presented here, slab subduction, melting and underplating processes generated arc magmatism and subsequent charnockitization within the BRB between ca. 2650 Ma and ca. 2498 Ma. These results thus reveal Meso- to Neoarchean tectonic evolution of the BRB. The spatial variation of crustal thickness, derived from flexure inversion technique, provides additional constraints on the tectonic linkage of the BRB with its surrounding terrains. In conjunction with published data, the Moyar and the Kollegal suture zones are considered to mark the trace of ocean closure along which the Nilgiri and Biligiri Rangan Blocks accreted on to the Western Dharwar Craton. (C) 2016 Elsevier B.V. All rights reserved.

Stalagmite growth perturbations from the Kumaun Himalaya as potential earthquake recorders

The central part of the Himalaya (Kumaun and Garhwal Provinces of India) is noted for its prolonged seismic quiescence, and therefore, developing a longer-term time series of past earthquakes to understand their recurrence pattern in this segment assumes importance. In addition to direct observations of offsets in stratigraphic exposures or other proxies like paleoliquefaction, deformation preserved within stalagmites (speleothems) in karst system can be analyzed to obtain continuous millennial scale time series of earthquakes. The Central Indian Himalaya hosts natural caves between major active thrusts forming potential storehouses for paleoseismological records. Here, we present results from the limestone caves in the Kumaun Himalaya and discuss the implications of growth perturbations identified in the stalagmites as possible earthquake recorders. This article focuses on three stalagmites from the Dharamjali Cave located in the eastern Kumaun Himalaya, although two other caves, one of them located in the foothills, were also examined for their suitability. The growth anomalies in stalagmites include abrupt tilting or rotation of growth axes, growth termination, and breakage followed by regrowth. The U-Th age data from three specimens allow us to constrain the intervals of growth anomalies, and these were dated at 4273 +/- 410 years BP (2673-1853 BC), 2782 +/- 79 years BP (851-693 BC), 2498 +/- 117 years BP (605-371 BC), 1503 +/- 245 years BP (262-752 AD), 1346 +/- 101 years BP (563-765 AD), and 687 +/- 147 years BP (1176-1470 AD). The dates may correspond to the timings of major/great earthquakes in the region and the youngest event (1176-1470 AD) shows chronological correspondence with either one of the great medieval earthquakes (1050-1250 and 1259-1433 AD) evident from trench excavations across the Himalayan Frontal Thrust.