Causality analysis of groundwater dynamics based on a Vector Autoregressive model in the semi-arid basin of Gundal (South India)

Causal relationships existing between observed levels of groundwater in a semi-arid sub-basin of the Kabini River basin (Karnataka state, India) are investigated in this study. A Vector Auto Regressive model is used for this purpose. Its structure is built on an upstream/downstream interaction network based on observed hydro-physical properties. Exogenous climatic forcing is used as an input based on cumulated rainfall departure. Optimal models are obtained thanks to a trial approach and are used as a proxy of the dynamics to derive causal networks. It appears to be an interesting tool for analysing the causal relationships existing inside the basin. The causal network reveals 3 main regions: the Northeastern part of the Gundal basin is closely coupled to the outlet dynamics. The Northwestern part is mainly controlled by the climatic forcing and only marginally linked to the outlet dynamic. Finally, the upper part of the basin plays as a forcing rather than a coupling with the lower part of the basin allowing for a separate analysis of this local behaviour. The analysis also reveals differential time scales at work inside the basin when comparing upstream oriented with downstream oriented causalities. In the upper part of the basin, time delays are close to 2 months in the upward direction and lower than 1 month in the downward direction. These time scales are likely to be good indicators of the hydraulic response time of the basin which is a parameter usually difficult to estimate practically. This suggests that, at the sub-basin scale, intra-annual time scales would be more relevant scales for analysing or modelling tropical basin dynamics in hard rock (granitic and gneissic) aquifers ubiquitous in south India. (c) 2012 Elsevier B.V. All rights reserved.

Ganga-Brahmaputra river discharge from Jason-2 radar altimetry: An update to the long-term satellite-derived estimates of continental freshwater forcing flux into the Bay of Bengal

This paper discusses the use of Jason-2 radar altimeter measurements to estimate the Ganga-Brahmaputra surface freshwater flux into the Bay of Bengal for the period mid-2008 to December 2011. A previous estimate was generated for 1993-2008 using TOPEX-Poseidon, ERS-2 and ENVISAT, and is now extended using Jason-2. To take full advantages of the new availability of in situ rating curves, the processing scheme is adapted and the adjustments of the methodology are discussed here. First, using a large sample of in situ river height measurements, we estimate the standard error of Jason-2-derived water levels over the Ganga and the Brahmaputra to be respectively of 0.28 m and 0.19 m, or less than similar to 4% of the annual peak-to-peak variations of these two rivers. Using the in situ rating curves between water levels and river discharges, we show that Jason-2 accurately infers Ganga and Brahmaputra instantaneous discharges for 2008-2011 with mean errors ranging from similar to 2180 m(3)/s (6.5%) over the Brahmaputra to similar to 1458 m(3)/s (13%) over the Ganga. The combined Ganga-Brahmaputra monthly discharges meet the requirements of acceptable accuracy (15-20%) with a mean error of similar to 16% for 2009-2011 and similar to 17% for 1993-2011. The Ganga-Brahmaputra monthly discharge at the river mouths is then presented, showing a marked interannual variability with a standard deviation of similar to 12500 m(3)/s, much larger than the data set uncertainty. Finally, using in situ sea surface salinity observations, we illustrate the possible impact of extreme continental freshwater discharge event on the northern Bay of Bengal as observed in 2008.

Impact of pit-toilet leachate on groundwater chemistry and role of vadose zone in removal of nitrate and E-coli pollutants in Kolar District, Karnataka, India

Assessment of chemistry of groundwater infiltrated by pit-toilet leachate and contaminant removal by vadose zone form the focus of this study. The study area is Mulbagal Town in Karnataka State, India. Groundwater level measurements and estimation of unsaturated permeability indicated that the leachate recharged the groundwater inside the town at the rate of 1 m/day. The average nitrate concentration of groundwater inside the town (148 mg/L) was three times larger than the permissible limit (45 mg/L), while the average nitrate concentration of groundwater outside the town (30 mg/L) was below the permissible limit. The groundwater inside the town exhibited E. coli contamination, while groundwater outside the town was free of pathogen contamination. Infiltration of alkalis (Na+, K+) and strong acids (Cl-, SO4 (2-)) caused the mixed Ca-Mg-Cl type (60 %) and Na-Cl type (28 %) facies to predominate groundwater inside the town, while, Ca-HCO3 (35 %), mixed Ca-Mg-Cl type (35 %) and mixed Ca-Na-HCO3 type (28 %) facies predominated groundwater outside/periphery of town. Reductions in E. coli and nitrate concentrations with vadose zone thickness indicated its participation in contaminant removal. A 4-m thickness of unsaturated sand + soft, disintegrated weathered rock deposit facilitates the removal of 1 log of E. coli pathogen. The anoxic conditions prevailing in the deeper layers of the vadose zone (> 19 m thickness) favor denitrification resulting in lower nitrate concentrations (28-96 mg/L) in deeper water tables (located at depths of -29 to -39 m).

Prediction of Ground Water Levels in the Uplands of a Tropical Coastal Riparian Wetland using Artificial Neural Networks

Artificial Neural Networks (ANNs) have been found to be a robust tool to model many non-linear hydrological processes. The present study aims at evaluating the performance of ANN in simulating and predicting ground water levels in the uplands of a tropical coastal riparian wetland. The study involves comparison of two network architectures, Feed Forward Neural Network (FFNN) and Recurrent Neural Network (RNN) trained under five algorithms namely Levenberg Marquardt algorithm, Resilient Back propagation algorithm, BFGS Quasi Newton algorithm, Scaled Conjugate Gradient algorithm, and Fletcher Reeves Conjugate Gradient algorithm by simulating the water levels in a well in the study area. The study is analyzed in two cases-one with four inputs to the networks and two with eight inputs to the networks. The two networks-five algorithms in both the cases are compared to determine the best performing combination that could simulate and predict the process satisfactorily. Ad Hoc (Trial and Error) method is followed in optimizing network structure in all cases. On the whole, it is noticed from the results that the Artificial Neural Networks have simulated and predicted the water levels in the well with fair accuracy. This is evident from low values of Normalized Root Mean Square Error and Relative Root Mean Square Error and high values of Nash-Sutcliffe Efficiency Index and Correlation Coefficient (which are taken as the performance measures to calibrate the networks) calculated after the analysis. On comparison of ground water levels predicted with those at the observation well, FFNN trained with Fletcher Reeves Conjugate Gradient algorithm taken four inputs has outperformed all other combinations.

Short- and long-term temporal variations in salinity and the oxygen, carbon and hydrogen isotopic compositions of the Hooghly Estuary water, India

The Hooghly River estuary provides a unique experimental site to understand the effect of monsoonal river discharge on freshwater and seawater mixing. Water samples collected bi-weekly for a duration of 17 months were analyzed for salinity, delta O-18,delta C-13(DIC), as well as delta D to investigate the differential mixing of freshwater and seawater. The differences in salinity and delta O-18 of samples collected during low and high tides on the same day are strongly correlated suggesting a well mixed water column at our sampling site. Low salinity and depleted delta O-18 during monsoon is consistent with increased river discharge as well as high rainfall. We identified different slopes in a delta O-18 versus salinity plot for the estuary water samples collected during monsoon and non-monsoon seasons. This is driven by composition of the freshwater source which is dominated by rainwater during monsoon and rivers during non-monsoon months. Selected delta D analyses of samples indicate that groundwater contributes significantly to the Hooghly Estuary during low rainfall times of the year. delta C-13(DIC) measured in the water recorded low values towards the end of monsoon indicating low productivity (i.e. increased organic respiration) while progressively increasing delta C-13(DIC) values from October till January as well as during some of the pre-monsoon months can be explained by increasing productivity. Very low delta C-13(DIC) (similar to-20%0) suggests involvement of carbon derived from anaerobic oxidation of organics and/or methane with potential contribution from increased anthropogenic water supply. An estimate of seawater incursion into the Hooghly Estuary at different times of the year is obtained by using salinity data in a two-component mixing model. Presence of seawater was found maximum (31-37%) during February till July and lowest (less than or equal to 6%) from September till November. We notice a temporal offset between Ganges River discharge farther upstream at Farakka and salinity variation at the Hooghly Estuary. We believe that this time lag is a result of the physical distance between Farakka and Kakdweep (our sampling location) and put constraints on the travel time of river water during early monsoon. (c) 2012 Published by Elsevier B.V.

Influence of anthropogenic contamination on groundwater chemistry in Mulbagal town, Kolar District, India

Geochemical studies are performed to examine the impact of leachate infiltration from on-site sewage disposal systems on the groundwater chemistry in Mulbagal town, Kolar District, Karnataka State, India. The leachate infiltration imposed nitrate concentrations ranging from 4 mg/L to 388 mg/L in the groundwater samples; it was observed that 79% of the samples exhibited nitrate concentrations in excess of drinking water permissible limit (45 mg/L). The average (of 43 measurements) E. coli levels in the groundwater samples corresponded to 189 MPN/100 mL and 55% of the samples tested exhibit pathogen contamination. Results also showed that the groundwater in the study area is characterized by acidic pH, large calcium + magnesium ion and Na/Cl ratios of < unity causing majority of the ground water samples to classify as Ca-Mg-Cl type and Na-Cl type. Saturation index (SI) computation using Visual MINTEQ program showed that the groundwater samples are under-saturated with respect to calcite. The theoretical SI values (of calcite) however suggested that the groundwater samples ought to be over-saturated with calcite. Under-saturation of the groundwater samples with calcite is attributed to increased dissolution of the mineral in the acidic environment of the groundwater.

Sea-level changes and its impact on coastal archaeological monuments: seven pagodas of Mahabalipuram, a case study

The name `Seven Pagodas' has served as a nickname for the south Indian port of Mahabalipuram since the early European explorers used it as landmark for navigation as they could see summits of seven temples from the sea. There are many theories concerning the name Seven Pagodas. The present study has compared coastline and adjacent seven monuments illustrated in a 17th century Portolan Chart (maritime map) with recent remote sensing data. This analysis throws new light on the name ``Seven Pagodas'' for the city. This study has used DEM of the site to simulate the coastline which is similar to the one depicted in the old portolan chart. Through this, the then sea level and corresponding flooding extent according to topography of the area and their effect on monuments could be analyzed. Most importantly this work has in the process identified possibly the seven monuments that constituted the name Seven Pagodas and this provides an alternative explanation to one of the mysteries of history. This work has demonstrated unique method of studying coastal archaeological sites. As large numbers of heritage sites around the world are on coastlines, this methodology has potential to be very useful for coastal heritage preservation and management.

Improved modeling of groundwater recharge in agricultural watersheds using a combination of crop model and remote sensing

For improved water management and efficiency of use in agriculture, studies dealing with coupled crop-surface water-groundwater models are needed. Such integrated models of crop and hydrology can provide accurate quantification of spatio-temporal variations of water balance parameters such as soil moisture store, evapotranspiration and recharge in a catchment. Performance of a coupled crop-hydrology model would depend on the availability of a calibrated crop model for various irrigated/rainfed crops and also on an accurate knowledge of soil hydraulic parameters in the catchment at relevant scale. Moreover, such a coupled model should be designed so as to enable the use/assimilation of recent satellite remote sensing products (optical and microwave) in order to model the processes at catchment scales. In this study we present a framework to couple a crop model with a groundwater model for applications to irrigated groundwater agricultural systems. We discuss the calibration of the STICS crop model and present a methodology to estimate the soil hydraulic parameters by inversion of crop model using both ground and satellite based data. Using this methodology we demonstrate the feasibility of estimation of potential recharge due to spatially varying soil/crop matrix.

Bearing Capacity of Foundations with Inclined Groundwater Seepage

The ultimate bearing capacity of strip foundations in the presence of inclined groundwater flow, considering both upward and downward flow directions, has been determined by using the lower bound finite-element limit analysis. A numerical solution has been generated for both smooth and rough footings placed on frictional soils. A correction factor (f gamma), which needs to be multiplied with the N gamma-term, has been computed to account for groundwater seepage. The variation of f gamma has been obtained as a function of the hydraulic gradient (i) for various inclinations of groundwater flow. For a given magnitude of i, there exists a certain critical inclination of the flow for which the value of f gamma is minimized. With an upward flow, for all flow inclinations, the magnitude of f gamma always reduces with an increase in the value of i. An example has also been provided to illustrate the application of the obtained results when designing foundations in the presence of groundwater seepage.

Constraining complex aquifer geometry with geophysics (2-D ERT and MRS measurements) for stochastic modelling of groundwater flow

Stochastic modelling is a useful way of simulating complex hard-rock aquifers as hydrological properties (permeability, porosity etc.) can be described using random variables with known statistics. However, very few studies have assessed the influence of topological uncertainty (i.e. the variability of thickness of conductive zones in the aquifer), probably because it is not easy to retrieve accurate statistics of the aquifer geometry, especially in hard rock context. In this paper, we assessed the potential of using geophysical surveys to describe the geometry of a hard rock-aquifer in a stochastic modelling framework. The study site was a small experimental watershed in South India, where the aquifer consisted of a clayey to loamy-sandy zone (regolith) underlain by a conductive fissured rock layer (protolith) and the unweathered gneiss (bedrock) at the bottom. The spatial variability of the thickness of the regolith and fissured layers was estimated by electrical resistivity tomography (ERT) profiles, which were performed along a few cross sections in the watershed. For stochastic analysis using Monte Carlo simulation, the generated random layer thickness was made conditional to the available data from the geophysics. In order to simulate steady state flow in the irregular domain with variable geometry, we used an isoparametric finite element method to discretize the flow equation over an unstructured grid with irregular hexahedral elements. The results indicated that the spatial variability of the layer thickness had a significant effect on reducing the simulated effective steady seepage flux and that using the conditional simulations reduced the uncertainty of the simulated seepage flux. As a conclusion, combining information on the aquifer geometry obtained from geophysical surveys with stochastic modelling is a promising methodology to improve the simulation of groundwater flow in complex hard-rock aquifers. (C) 2013 Elsevier B.V. All rights reserved.

Modelling the impact of extensive irrigation on the groundwater resources

Drastic groundwater resource depletion due to excessive extraction for irrigation is a major concern in many parts of India. In this study, an attempt was made to simulate the groundwater scenario of the catchment using ArcSWAT. Due to the restriction on the maximum initial storage, the deep aquifer component in ArcSWAT was found to be insufficient to represent the excessive groundwater depletion scenario. Hence, a separate water balance model was used for simulating the deep aquifer water table. This approach is demonstrated through a case study for the Malaprabha catchment in India. Multi-site rainfall data was used to represent the spatial variation in the catchment climatology. Model parameters were calibrated using observed monthly stream flow data. Groundwater table simulation was validated using the qualitative information available from the field. The stream flow was found to be well simulated in the model. The simulated groundwater table fluctuation is also matching reasonably well with the field observations. From the model simulations, deep aquifer water table fluctuation was found very severe in the semi-arid lower parts of the catchment, with some areas showing around 60m depletion over a period of eight years. Copyright (c) 2012 John Wiley & Sons, Ltd.

Bearing capacity of foundations subjected to groundwater flow

The ultimate bearing capacity of strip foundations subjected to horizontal groundwater flow has been computed by making use of the stress characteristics method which is well known for its capability in solving quite accurately different stability problems in geotechnical engineering. The numerical solution has been generated both for smooth and rough footings placed on frictional soils. A correction factor (fγ) associated with Nγ term, to account for the existence of ground water flow, has been introduced. The variation of fγ has been obtained as a function of hydraulic gradient (i) for different values of soil frictional angle. The magnitude of fγ reduces continuously with an increase in the value of i.

Influence of exhaust gas heating and L/D ratios on the discharge efficiencies for an activated carbon natural gas storage system

A transient 2D axi-symmetric and lumped parameter (LP) model with constant outflow conditions have been developed to study the discharge capacity of an activated carbon bed. The predicted discharge times and variations in bed pressure and temperature are in good agreement with experimental results obtained from a 1.82 l adsorbed natural gas (ANG) storage system. Under ambient air conditions, a maximum temperature drop of 29.5 K and 45.5 K are predicted at the bed center for discharge rates of 1.0 l min(-1) and 5.0 l min(-1) respectively. The corresponding discharge efficiencies are 77% and 71.5% respectively with discharge efficiencies improving with decreasing outflow rates. Increasing the LID ratio from 1.9 to 7.8 had only a marginal increase in the discharge efficiency. Forced convection (exhaust gas) heating had a significant effect on the discharge efficiency, leading to efficiencies as high as 92.8% at a discharge of 1.0 l min(-1) and 88.7% at 5 l min(-1). Our study shows that the LP model can be reliably used to obtain discharge times due to the uniform pressure distributions in the bed. Temperature predictions with the LP model were more accurate at ambient conditions and higher discharge rates, due to greater uniformity in bed temperatures. For the low thermal conductivity carbon porous beds, our study shows that exhaust gas heating can be used as an effective and convenient strategy to improve the discharge characteristics and performance of an ANG system. (C) 2013 Elsevier Ltd. All rights reserved.

Influence of anthropogenic contamination on fluoride concentration in groundwater: A study of Mulbagal town, Kolar district, Karnataka

Groundwater contamination is a serious concern in India. Major geogenic contaminants include fluoride, arsenic and iron, while common anthropogenic contaminants include nitrate, metals, organics and microbial contamination. Besides, known point and diffuse sources, groundwater c ontamination from inf iltration of pit to ilet leachate is an emerging concern. The study area of this paper is Kolar district in Karnataka that is hot spot of fluoride contamination. The absence of fluoride contamination in Mulbagal town and the alterations in groundwater chemistry from infiltration of pit toilet leachate motivated the author to examine the possible linkages between anthropogenic contamination and fluoride concentration in groundwater of Mulbagal town. Analysis of the groundwater chemistry revealed that the groundwater in Mulbagal town is under saturated with respect to calcite that suppresses the dissolution of fluorite and the fluoride concentration in the groundwater. The slightly acidic pH of the groundwater is considered responsible to facilitate calcite dissolution under saturation.

Hydrographical characteristics and oxygen isotopic signatures of water in a coastal environment (Mangalore) along the southeastern Arabian Sea

Coastal marine environments are important links between the continents and the open ocean. The coast off Mangalore forms part of the upwelling zone along the southeastern Arabian Sea. The temperature, salinity, density, dissolved oxygen and stable oxygen isotope ratio (delta O-18) of surface waters as well as those of bottom waters off coastal Mangalore were studied every month from October 2010 to May 2011. The coastal waters were stratified in October and November due to precipitation and runoff. The region was characterised by upwelled bottom waters in October, whereas the region exhibited a temperature inversion in November. The surface and bottom waters presented almost uniform properties from December until April. The coastal waters were observed to be most dense in January and May. Comparatively cold and poorly oxygenated bottom waters during the May sampling indicated the onset of upwelling along the region. delta O-18 of the coastal waters successfully documented the observed variations in the hydrographical characteristics of the Mangalore coast during the monthly sampling period. We also noted that the monthly variability in the properties of the coastal waters of Mangalore was related to the hydrographical characteristics of the adjacent open ocean inferred from satellite-derived surface winds, sea surface height anomaly data and sea surface temperatures.