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.

ON POLYTOPAL UPPER BOUND SPHERES

Generalizing a result (the case k = 1) due to M. A. Perles, we show that any polytopal upper bound sphere of odd dimension 2k + 1 belongs to the generalized Walkup class K-k(2k + 1), i.e., all its vertex links are k-stacked spheres. This is surprising since it is far from obvious that the vertex links of polytopal upper bound spheres should have any special combinatorial structure. It has been conjectured that for d not equal 2k + 1, all (k + 1)-neighborly members of the class K-k(d) are tight. The result of this paper shows that the hypothesis d not equal 2k + 1 is essential for every value of k >= 1.

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.

Stability of long unsupported twin circular tunnels in soils

The stability of two long unsupported circular parallel tunnels aligned horizontally in fully cohesive and cohesive-frictional soils has been determined. An upper bound limit analysis in combination with finite elements and linear programming is employed to perform the analysis. For different clear spacing (S) between the tunnels, the stability of tunnels is expressed in terms of a non-dimensional stability number (gamma H-max/c); where H is tunnel cover, c refers to soil cohesion, and gamma(max) is maximum unit weight of soil mass which the tunnels can bear without any collapse. The variation of the stability number with tunnels' spacing has been established for different combinations of H/D, m and phi; where D refers to diameter of each tunnel, phi is the internal friction angle of soil and m accounts for the rate at which the cohesion increases linearly with depth. The stability number reduces continuously with a decrease in the spacing between the tunnels. The optimum spacing (S-opt) between the two tunnels required to eliminate the interference effect increases with (i) an increase in H/D and (ii) a decrease in the values of both m and phi. The value of S-opt lies approximately in a range of 1.5D-3.5D with H/D = 1 and 7D-12D with H/D = 7. The results from the analysis compare reasonably well with the different solutions reported in literature. (C) 2013 Elsevier Ltd. All rights reserved.

The influence of debris cover and glacial lakes on the recession of glaciers in Sikkim Himalaya, India

We investigated area changes in glaciers covering an area of similar to 200 km(2) in the Tista basin, Sikkim, Eastern Indian Himalaya, between similar to 1990 and 2010 using Landsat Thematic Mapper (TM) and Indian Remote-sensing Satellite (IRS) images and related the changes to debris cover, supraglacial lakes and moraine-dam lakes. The glaciers lost an area of 3.3 +/- 0.8% between 1989/90 and 2010. More detailed analysis revealed an area loss of 2.00 +/- 0.82, 2.56 +/- 0.61 and 2.28 +/- 2.01 km(2) for the periods 1989-97, 1997-2004/05 and 2004-2009/10, respectively. This indicates an accelerated retreat of glaciers after 1997. On further analysis, we observed (1) the formation and expansion of supraglacial lakes on many debris-covered glaciers and (2) the merging of these lakes over time, leading to the development of large moraine-dam lakes. We also observed that debris-covered glaciers with lakes lose a greater area than debris-covered glaciers without lakes and debris-free glaciers. The climatic data for 24 years (1987-2011), measured at the Gangtok meteorological station (1812 m a.s.l.), showed that the region experienced a 1.0 degrees C rise in the summer minimum temperature and a 2.0 degrees C rise in the winter minimum temperature, indicating hotter summers and warmer winters. There was no significant trend in the total annual precipitation. We find that glacier retreat is caused mainly by a temperature increase and that debris-covered glaciers can retreat at a faster rate than debris-free glaciers, if associated with lakes.

Vertical uplift resistance of two horizontal strip anchors with common vertical axis

With an application of the upper bound finite element limit analysis, the vertical pullout capacity of a group of two horizontal strip plate anchors, with the common vertical axis and placed in a cohesive-frictional soil, has been computed. The variation of the uplift factors Fc, Fq and Fy, due to the contributions of soil cohesion, surcharge pressure and unit weight, respectively, has been evaluated for different combinations of S/B and H/B. As compared to single isolated anchor, the group of two anchors generates significantly greater magnitude of Fc for Φ ≤ 20° especially with greater values of H/B and under fully bonded anchor-soil interface condition. The factor Fc attains almost the maximum value when the upper anchor plate is placed midway between ground surface and the lower anchor plate. The factors Fq and Fy, on the other hand, for a group of two anchors are found to remain almost equal to that of a single isolated anchor as long as the levels of the lower plate in the group and the single isolated anchor are kept the same.

Upper bound solution for pullout capacity of vertical anchors in sand using finite elements and limit analysis

The horizontal pullout capacity of vertical anchors embedded in sand has been determined by using an upper bound theorem of the limit analysis in combination with finite elements. The numerical results are presented in nondimensional form to determine the pullout resistance for various combinations of embedment ratio of the anchor (H/B), internal friction angle (ϕ) of sand, and the anchor-soil interface friction angle (δ). The pullout resistance increases with increases in the values of embedment ratio, friction angle of sand and anchor-soil interface friction angle. As compared to earlier reported solutions in literature, the present solution provides a better upper bound on the ultimate collapse load.

Development of a regional model for the North Indian Ocean

We have developed a one-way nested Indian Ocean regional model. The model combines the National Oceanic and Atmospheric Administration (NOAA) Geophysical Fluid Dynamics Laboratory's (GFDL) Modular Ocean Model (MOM4p1) at global climate model resolution (nominally one degree), and a regional Indian Ocean MOM4p1 configuration with 25 km horizontal resolution and 1 m vertical resolution near the surface. Inter-annual global simulations with Coordinated Ocean-Ice Reference Experiments (CORE-II) surface forcing over years 1992-2005 provide surface boundary conditions. We show that relative to the global simulation, (i) biases in upper ocean temperature, salinity and mixed layer depth are reduced, (ii) sea surface height and upper ocean circulation are closer to observations, and (iii) improvements in model simulation can be attributed to refined resolution, more realistic topography and inclusion of seasonal river runoff. Notably, the surface salinity bias is reduced to less than 0.1 psu over the Bay of Bengal using relatively weak restoring to observations, and the model simulates the strong, shallow halocline often observed in the North Bay of Bengal. There is marked improvement in subsurface salinity and temperature, as well as mixed layer depth in the Bay of Bengal. Major seasonal signatures in observed sea surface height anomaly in the tropical Indian Ocean, including the coastal waveguide around the Indian peninsula, are simulated with great fidelity. The use of realistic topography and seasonal river runoff brings the three dimensional structure of the East India Coastal Current and West India Coastal Current much closer to observations. As a result, the incursion of low salinity Bay of Bengal water into the southeastern Arabian Sea is more realistic. (C) 2013 Elsevier Ltd. All rights reserved.

Evolution of superconductivity in PrFe1-xCoxAsO (x=0.0-1.0)

We report the synthesis and physical property characterization of Prfe(1-x)Co(x)AsO (x=0.0-1.0). The studied samples are synthesized by through the solid state reaction route via the vacuum encapsulation method. The pristine compound PrFeAsO does not show superconductivity, but rather exhibits a metallic step like transition due to spin density wave (SOW) ordering of Fe moments (Fe-SDW) below 150 K, Followed by another upward step due to anomalous ordering of Pr moments (Pr-TN) at 12 K. Both the Fe-SDW and Pr-TN temperatures decrease monotonically with Co substitution at Fe site Superconductivity appears in a narrow range of x from 0.07 to 0.25 with maximum T-c at 11.12 K for x=0.15. Samples with x >= 0.25 exhibit metallic behavior right from 300 K down to 2 K, without any Fe-SDW or Pr-TN steps in resistivity. In fact, though Fe-SDW decreases monotonically, the pr(TN) disappeared even with x=0.02. The magneto transport measurements below 14 Ton superconducting polycrystalline Co doped Pi FeAs0 lead to extrapolated values of the upper critical fields H-c2(0)] of up to 60 T. (C) 2014 Elsevier Ltd. All rights reserved.

Solving Axisymmetric Stability Problems by Using Upper Bound Finite Elements, Limit Analysis, and Linear Optimization

A numerical formulation has been proposed for solving an axisymmetric stability problem in geomechanics with upper bound limit analysis, finite elements, and linear optimization. The Drucker-Prager yield criterion is linearized by simulating a sphere with a circumscribed truncated icosahedron. The analysis considers only the velocities and plastic multiplier rates, not the stresses, as the basic unknowns. The formulation is simple to implement, and it has been employed for finding the collapse loads of a circular footing placed over the surface of a cohesive-frictional material. The formulation can be used to solve any general axisymmetric geomechanics stability problem.

Effect of the Flow Rule on the Bearing Capacity of Strip Foundations on Sand by the Upper-Bound Limit Analysis and Slip Lines

The influence of the flow rule on the bearing capacity of strip foundations placed on sand was investigated using a new kinematic approach of upper-bound limit analysis. The method of stress characteristics was first used to find the mechanism of the failure and to compute the stress field by using the Mohr-Coulomb yield criterion. Once the failure mechanism had been established, the kinematics of the plastic deformation was established, based on the requirements of the upper-bound limit theorem. Both associated and nonassociated plastic flows were considered, and the bearing capacity was obtained by equating the rate of external plastic work to the rate of the internal energy dissipation for both smooth and rough base foundations. The results obtained from the analysis were compared with those available from the literature. (C) 2014 American Society of Civil Engineers.

Stability of a circular tunnel in presence of pseudostatic seismic body forces

The stability of a long circular tunnel in a cohesive frictional soil medium has been determined in the presence of horizontal pseudo-static seismic body forces. The tunnel is supported by means of lining and anchorage system which is assumed to exert uniform internal compressive normal pressure on its periphery. The upper bound finite element limit analysis has been performed to compute the magnitude of the internal compressive pressure required to support the tunnel. The results have been presented in terms of normalized compressive normal stress, defined in terms of sigma(i)/c; where sigma(i) is the magnitude of the compressive normal pressure on the periphery of the tunnel and c refers to soil cohesion. The variation of sigma(i)/c with horizontal earthquake acceleration coefficient (alpha(h)) has been established for different combinations of H/D, gamma D/c and phi where (i) H and D refers to tunnel cover and diameter, respectively, and (ii) gamma and phi correspond to unit weight and internal friction angle of soil mass, respectively. Nodal velocity patterns have also been plotted for assessing the zones of significant plastic deformation. The analysis clearly reveals that an increase in the magnitude of the earthquake acceleration leads to a significant increment in the magnitude of internal compressive pressure. (C) 2014 Elsevier Ltd. All rights reserved.

Performance Analysis of a Cooperative MAC Protocol of Wireless Ad Hoc Networks

Different medium access control (MAC) layer protocols, for example, IEEE 802.11 series and others are used in wireless local area networks. They have limitation in handling bulk data transfer applications, like video-on-demand, videoconference, etc. To avoid this problem a cooperative MAC protocol environment has been introduced, which enables the MAC protocol of a node to use its nearby nodes MAC protocol as and when required. We have found on various occasions that specified cooperative MAC establishes cooperative transmissions to send the specified data to the destination. In this paper we propose cooperative MAC priority (CoopMACPri) protocol which exploits the advantages of priority value given by the upper layers for selection of different paths to nodes running heterogeneous applications in a wireless ad hoc network environment. The CoopMACPri protocol improves the system throughput and minimizes energy consumption. Using a Markov chain model, we developed a model to analyse the performance of CoopMACPri protocol; and also derived closed-form expression of saturated system throughput and energy consumption. Performance evaluations validate the accuracy of the theoretical analysis, and also show that the performance of CoopMACPri protocol varies with the number of nodes. We observed that the simulation results and analysis reflects the effectiveness of the proposed protocol as per the specifications.

Upper Bounds on the Size of Grain-Correcting Codes

In this paper, we revisit the combinatorial error model of Mazumdar et al. that models errors in high-density magnetic recording caused by lack of knowledge of grain boundaries in the recording medium. We present new upper bounds on the cardinality/rate of binary block codes that correct errors within this model. All our bounds, except for one, are obtained using combinatorial arguments based on hypergraph fractional coverings. The exception is a bound derived via an information-theoretic argument. Our bounds significantly improve upon existing bounds from the prior literature.

Regimes of nonlinear depletion and regularity in the 3D Navier-Stokes equations

The periodic 3D Navier-Stokes equations are analyzed in terms of dimensionless, scaled, L-2m-norms of vorticity D-m (1 <= m <= infinity). The first in this hierarchy, D-1, is the global enstrophy. Three regimes naturally occur in the D-1-D-m plane. Solutions in the first regime, which lie between two concave curves, are shown to be regular, owing to strong nonlinear depletion. Moreover, numerical experiments have suggested, so far, that all dynamics lie in this heavily depleted regime 1]; new numerical evidence for this is presented. Estimates for the dimension of a global attractor and a corresponding inertial range are given for this regime. However, two more regimes can theoretically exist. In the second, which lies between the upper concave curve and a line, the depletion is insufficient to regularize solutions, so no more than Leray's weak solutions exist. In the third, which lies above this line, solutions are regular, but correspond to extreme initial conditions. The paper ends with a discussion on the possibility of transition between these regimes.