Characterization of groundwater chemistry under the influence of lithologic and anthropogenic factors along a climatic gradient in Upper Cauvery basin, South India

Hydrogeological and climatic effect on chemical behavior of groundwater along a climatic gradient is studied along a river basin. `Semi-arid' (500-800 mm of mean annual rainfall), `sub-humid' (800-1,200 mm/year) and `humid' (1,200-1,500 mm/year) are the climatic zones chosen along the granito-gneissic plains of Kabini basin in South India for the present analysis. Data on groundwater chemistry is initially checked for its quality using NICB ratio (<+/- 5 %), EC versus TZ+ (similar to 0.85 correlation), EC versus TDS and EC versus TH analysis. Groundwater in the three climatic zones is `hard' to `very hard' in terms of Ca-Mg hardness. Polluted wells are identified (> 40 % of pollution) and eliminated for the characterization. Piper's diagram with mean concentrations indicates the evolution of CaNaHCO3 (semi-arid) from CaHCO3 (humid zone) along the climatic gradient. Carbonates dominate other anions and strong acids exceeded weak acids in the region. Mule Hole SEW, an experimental watershed in sub-humid zone, is characterized initially using hydrogeochemistry and is observed to be a replica of entire sub-humid zone (with 25 wells). Extension of the studies for the entire basin (120 wells) showed a chemical gradient along the climatic gradient with sub-humid zone bridging semi-arid and humid zones. Ca/Na molar ratio varies by more than 100 times from semi-arid to humid zones. Semi-arid zone is more silicaceous than sub-humid while humid zone is more carbonaceous (Ca/Cl similar to 14). Along the climatic gradient, groundwater is undersaturated (humid), saturated (sub-humid) and slightly supersaturated (semi-arid) with calcite and dolomite. Concentration-depth profiles are in support of the geological stratification i.e., not approximate to 18 m of saprolite and similar to 25 m of fracture rock with parent gneiss beneath. All the wells are classified into four groups based on groundwater fluctuations and further into `deep' and `shallow' based on the depth to groundwater. Higher the fluctuations, larger is its impact on groundwater chemistry. Actual seasonal patterns are identified using `recharge-discharge' concept based on rainfall intensity instead of traditional monsoon-non-monsoon concept. Non-pumped wells have low Na/Cl and Ca/Cl ratios in recharge period than in discharge period (Dilution). Few other wells, which are subjected to pumping, still exhibit dilution chemistry though water level fluctuations are high due to annual recharge. Other wells which do not receive sufficient rainfall and are constantly pumped showed high concentrations in recharge period rather than in discharge period (Anti-dilution). In summary, recharge-discharge concept demarcates the pumped wells from natural deep wells thus, characterizing the basin.

Comparison of CloudSat and TRMM radar reflectivities

Comparison of reflectivity data of radars onboard CloudSat and TRMM is performed using coincident overpasses. The contoured frequency by altitude diagrams (CFADs) are constructed for two cases: (a) only include collocated vertical profiles that are most likely to be raining and (b) include all collocated profiles along with cloudy pixels falling within a distance of about 50 km from the centre point of coincidence. Our analysis shows that for both cases, CloudSat underestimates the radar reflectivity by about 10 dBZ compared to that of TRMM radar below 15 km altitude. The difference is well outside the uncertainty value of similar to 2 dBZ of each radar. Further, CloudSat reflectivity shows a decreasing trend while that of TRMM radar an increasing trend below 4 km height. Basically W-band radar that CloudSat flies suffers strong attenuation in precipitating clouds and its reflectivity value rarely exceeds 20 dBZ though its technical specification indicates the upper measurement limit to be 40 dBZ. TRMM radar, on the other hand, cannot measure values below 17 dBZ. In fact combining data from these two radars seems to give a better overall spatial structure of convective clouds.

Numerical simulations of CO2 migration during charnockite genesis

Charnockite is considered to be generated either through the dehydration of granitic magma by CO2 purging or by solid-state dehydration through CO2 metasomatism during granulite facies metamorphism. To understand the extent of dehydration, CO2 migration is quantitatively modeled in silicate melt and metasomatic fluid as a function of temperature, H2O wt%, pressure, basal CO2 flux and dynamic viscosity. Numerical simulations show that CO2 advection through porous and permeable high-grade metamorphic rocks can generate dehydrated patches close to the CO2 flow path, as illustrated by the occurrences of ``incipient charnockites.'' CO2 reaction-front velocity constrained by field observations is 0.69 km/m.y., a reasonable value, which matches well with other studies. On the other hand, temperature, rate of cooling, and basal CO2 flux are the critical parameters affecting CO2 diffusion through a silicate melt. CO2 diffusion through silicate melt can only occur at temperature greater than 840 degrees C and during slow cooling (<= 3.7 x 10(-5) degrees C/yr), features that are typical of magma emplacement in the lower crust. Stalling of CO2 fluxing at similar to 840 degrees C explains why some deep-level plutons contain both hydrous and anhydrous (charnockitic) mineral assemblages. CO2 diffusion through silicate melt is virtually insensitive to pressure. Addition of CO2 basal flux facilitates episodic dehydrated melt migration by generating fracture pathways.

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

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.

Assessing GCM Convergence for India Using the Variable Convergence Score

General circulation models (GCMs) use transient climate simulations to predict climate conditions in the future. Coarse-grid resolutions and process uncertainties necessitate the use of downscaling models to simulate precipitation. However, in the downscaling models, with multiple GCMs now available, selecting an atmospheric variable from a particular model which is representative of the ensemble mean becomes an important consideration. The variable convergence score (VCS) provides a simple yet meaningful approach to address this issue, providing a mechanism to evaluate variables against each other with respect to the stability they exhibit in future climate simulations. In this study, VCS methodology is applied to 10 atmospheric variables of particular interest in downscaling precipitation over India and also on a regional basis. The nested bias-correction methodology is used to remove the systematic biases in the GCMs simulations, and a single VCS curve is developed for the entire country. The generated VCS curve is expected to assist in quantifying the variable performance across different GCMs, thus reducing the uncertainty in climate impact-assessment studies. The results indicate higher consistency across GCMs for pressure and temperature, and lower consistency for precipitation and related variables. Regional assessments, while broadly consistent with the overall results, indicate low convergence in atmospheric attributes for the Northeastern parts of India.

The role of buoyancy in polarity reversals of the geodynamo

We investigate polarity reversals in the geodynamo using a rotating, convection-driven dynamo model. As the flow in rapidly rotating convection is dominated by columns aligned with the axis of rotation, the focus is on the dynamics of columnar vortices. By studying the growth of a seed magnetic field to a stable axial dipole field, we show that the magnetic field acts in ways that significantly enhance the relative helicity between cyclonic and anticyclonic vortices. This flow asymmetry is the hallmark of a dipolar dynamo. Strong buoyancy, on the other hand, offsets the effect of the magnetic field, establishing parity between positive and negative vortices. As the dipole field is deprived of the helicity required to support itself, the dynamo is pushed into a reversing state. This is a likely regime for polarity reversals in the Earth's core. The integral lengthscale at which buoyancy injects energy is not significantly different from the convective flow lengthscale, which implies that buoyancy does not feed vortices at the small scales where non-linear inertia is present. The lengthscale at which the Lorentz force acts in the reversing dynamo is small, which may allow the passive presence of non-linear inertia in the small scales.

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.

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.

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.

Joint modeling of lithosphere and mantle dynamics: Evaluation of constraints from global tomography models

With the advances in technology, seismological theory, and data acquisition, a number of high-resolution seismic tomography models have been published. However, discrepancies between tomography models often arise from different theoretical treatments of seismic wave propagation, different inversion strategies, and different data sets. Using a fixed velocity-to-density scaling and a fixed radial viscosity profile, we compute global mantle flow models associated with the different tomography models and test the impact of these for explaining surface geophysical observations (geoid, dynamic topography, stress, and strain rates). We use the joint modeling of lithosphere and mantle dynamics approach of Ghosh and Holt (2012) to compute the full lithosphere stresses, except that we use HC for the mantle circulation model, which accounts for the primary flow-coupling features associated with density-driven mantle flow. Our results show that the seismic tomography models of S40RTS and SAW642AN provide a better match with surface observables on a global scale than other models tested. Both of these tomography models have important similarities, including upwellings located in Pacific, Eastern Africa, Iceland, and mid-ocean ridges in the Atlantic and Indian Ocean and downwelling flows mainly located beneath the Andes, the Middle East, and central and Southeast Asia.

Mechanical and thermal structure of a cylindrical plume in the earth's mantle

A mantle plume is understood as a hot, narrow, upwelling flow in the earth's mantle and accompanied by an efficient transfer of mass and energy from deep to upper layer of the earth. The cylindrical plume in earth's mantle plays an important role in explaining the origin of the surface hot spots and linear island chains. From the basic hydrodynamical equations, the detailed mechanical and thermal structure of a cylindrical plume of Newtouian fluids with temperature and pressure-dependent viscosity are given in the present paper. For two sets of rheological parameters the radial profiles of upward velocity, temperature and viscosity in the plume and radiuses of the plume at various depths have been calculated.

Distributions of temperature, velocity and viscosity at symmetrical axis of cylindrical plume in earth's mantle with variable viscosity

The parameters at the symmetrical axis of a cylindrical plume characterize the strength of this plume and provide a boundary condition which must be given to investigate the structure of a plume. For Newtonian fluid with a temperature-and pressure-dependence viscosity, an asymptotical solution of hydrodynamic equations at the symmetrical axis of the plume is found in the present paper. The temperature, upward velocity and viscosity at the symmetrical axis have been obtained as functions of depth, The calculated results have been given for two typical sets of Newtonian rheological parameters. The results obtained show that the temperature distribution along the symmetrical axis is nearly independent of the theological parameters. The upward velocity at the symmetrical axis, however, is strongly dependent on the rheological parameters.

Contributions to the geology and paleontology of the Canal Zone, Panama, and geologically related areas in Central America and the West Indies. Decapod crustaceans from the Panama region

Fifty-eight species of Decapods are enumerated from the collections examined by the author. Three species described by other authors are inserted in systematic order, thus making the list complete for the Panama region. All available material in the United States National Museum from Panama and Costa Rica is included; it ranges in age from the Oligocene (Culebra formation) to the Pleistocene.