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.

A method for predicting the consolidated undrained bearing capacity of shallow foundations

The effect of consolidation on the undrained bearing capacity of both rough and smooth strip and circular surface foundations is investigated, examining the influence of the magnitude and duration of an applied preload and the initial over-consolidation ratio of the deposit. The investigation comprised small strain finite-element analysis, with the soil response represented by Modified Cam Clay. The results are distilled into dimensionless and generalised forms, from which simple trends emerge. Based on these results, a simple method for predicting the consolidated undrained bearing capacity is proposed.

Do stable isotopes in carbonate cement of Mio-Pleistocene Himalayan sediments record paleoecological and paleoclimatic changes?

Non-pedogenic carbonates, such as carbonate cement and nodules in the sandstones, are quite common in the terrestrial geological record. Unlike pedogenic carbonates, their stable isotope ratios lack investigations for paleo-climatic reconstructions. The present investigation therefore, explores the possibility of use of stable isotope studies of non-pedogenic carbonates from the Mb-Pleistocene Siwalik Group of sediments exposed in the Ramnagar sub-basin of the NW Himalaya. Petrographic studies reveal the dominance of micrite fabric in carbonate nodules both of pedogenic and non-pedogenic samples irrespective of specific stratigraphic unit However, calcite as cement in the sandstones shows the dominance of micrite fabric in the younger in age sediments. Seventy-two non-pedogenic carbonate samples from the carbonate nodules and cement in the Siwalik sandstones, ranging in age between similar to 1 Ma and 12.2 Ma, were analyzed for delta C-13 and delta O-18 values. The delta C-13 values vary from -24.77 parts per thousand to -1.1 parts per thousand and delta O-18 values vary from -15.34 parts per thousand to -7.81 parts per thousand. Pedogenic and non-pedogenic carbonates ranging in age between similar to 1 Ma and 6 Ma have largely similar delta C-13 values and the range of delta C-13 values indicate the dominance of C-4 type of vegetation. However, unlike pedogenic carbonates which showed the dominance of C-3 type of vegetation pre- 7 Ma on the basis of delta C-13 -depleted isotopic values (Singh et al., 2011), delta C-13 values are largely enriched in the corresponding aged non-pedogenic carbonates revealing no information on specific type of vegetation. Likewise, paleoprecipitational reconstructions from delta O-18 values in pedogenic carbonates showed a progressive increase in aridity from similar to 12 Ma to recent excluding short term increases in rainfall/monsoon intensity at around 10 Ma, 5 Ma, and 1.8 Ma (Singh et al., 2012). On the contrary, such reconstructions are not possible from the delta O-18 values of non-pedogenic carbonates and indeed the delta O-18 values of non-pedogenic carbonates are largely depleted to as much as 6 parts per thousand from the corresponding pedogenic carbonates. Such differences in delta C-13 and delta O-18 values of non-pedogenic carbonates from pedogenic carbonates are primarily due to the dependence of the former on groundwater conditions responsible for precipitating carbonate. Further, a comparison of isotopic values between non-pedogenic and pedogenic carbonates can be interpreted that post-6 Ma and pre-6 Ma non-pedogenic carbonates were largely formed by shallow and deep groundwater conditions respectively. The result of these investigative studies therefore, suggests that the stable delta C-13 and delta O-18 values of non-pedogenic carbonates, unlike the pedogenic carbonates and irrespective of nature of calcite fabric, showed their little importance in paleoclimatic and paleoecological reconstructions. (C) 2014 Elsevier B.V. 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.

Three-dimensional localization of multiple acoustic sources in shallow ocean with non-Gaussian noise

In this paper, a low-complexity algorithm SAGE-USL is presented for 3-dimensional (3-D) localization of multiple acoustic sources in a shallow ocean with non-Gaussian ambient noise, using a vertical and a horizontal linear array of sensors. In the proposed method, noise is modeled as a Gaussian mixture. Initial estimates of the unknown parameters (source coordinates, signal waveforms and noise parameters) are obtained by known/conventional methods, and a generalized expectation maximization algorithm is used to update the initial estimates iteratively. Simulation results indicate that convergence is reached in a small number of (<= 10) iterations. Initialization requires one 2-D search and one 1-D search, and the iterative updates require a sequence of 1-D searches. Therefore the computational complexity of the SAGE-USL algorithm is lower than that of conventional techniques such as 3-D MUSIC by several orders of magnitude. We also derive the Cramer-Rao Bound (CRB) for 3-D localization of multiple sources in a range-independent ocean. Simulation results are presented to show that the root-mean-square localization errors of SAGE-USL are close to the corresponding CRBs and significantly lower than those of 3-D MUSIC. (C) 2014 Elsevier Inc. All rights reserved.

Dense shallow granular flows

Simplified equations are derived for a granular flow in the `dense' limit where the volume fraction is close to that for dynamical arrest, and the `shallow' limit where the stream-wise length for flow development (L) is large compared with the cross-stream height (h). The mass and diameter of the particles are set equal to 1 in the analysis without loss of generality. In the dense limit, the equations are simplified by taking advantage of the power-law divergence of the pair distribution function chi proportional to (phi(ad) - phi)(-alpha), and a faster divergence of the derivativ rho(d chi/d rho) similar to (d chi/d phi), where rho and phi are the density and volume fraction, and phi(ad) is the volume fraction for arrested dynamics. When the height h is much larger than the conduction length, the energy equation reduces to an algebraic balance between the rates of production and dissipation of energy, and the stress is proportional to the square of the strain rate (Bagnold law). In the shallow limit, the stress reduces to a simplified Bagnold stress, where all components of the stress are proportional to (partial derivative u(x)/partial derivative y)(2), which is the cross-stream (y) derivative of the stream-wise (x) velocity. In the simplified equations for dense shallow flows, the inertial terms are neglected in the y momentum equation in the shallow limit because the are O(h/L) smaller than the divergence of the stress. The resulting model contains two equations, a mass conservation equations which reduces to a solenoidal condition on the velocity in the incompressible limit, and a stream-wise momentum equation which contains just one parameter B which is a combination of the Bagnold coefficients and their derivatives with respect to volume fraction. The leading-order dense shallow flow equations, as well as the first correction due to density variations, are analysed for two representative flows. The first is the development from a plug flow to a fully developed Bagnold profile for the flow down an inclined plane. The analysis shows that the flow development length is ((rho) over barh(3)/B) , where (rho) over bar is the mean density, and this length is numerically estimated from previous simulation results. The second example is the development of the boundary layer at the base of the flow when a plug flow (with a slip condition at the base) encounters a rough base, in the limit where the momentum boundary layer thickness is small compared with the flow height. Analytical solutions can be found only when the stream-wise velocity far from the surface varies as x(F), where x is the stream-wise distance from the start of the rough base and F is an exponent. The boundary layer thickness increases as (l(2)x)(1/3) for all values of F, where the length scale l = root 2B/(rho) over bar. The analysis reveals important differences between granular flows and the flows of Newtonian fluids. The Reynolds number (ratio of inertial and viscous terms) turns out to depend only on the layer height and Bagnold coefficients, and is independent of the flow velocity, because both the inertial terms in the conservation equations and the divergence of the stress depend on the square of the velocity/velocity gradients. The compressibility number (ratio of the variation in volume fraction and mean volume fraction) is independent of the flow velocity and layer height, and depends only on the volume fraction and Bagnold coefficients.

Low-frequency modes in an equatorial shallow-water model with moisture gradients

Motivated by observations of the mean state of tropical precipitable water (PW), a moist, first baroclinic mode, shallow-water system on an equatorial beta-plane with a background saturation profile that depends on latitude and longitude is studied. In the presence of a latitudinal moisture gradient, linear analysis of the non-rotating problem reveals large-scale, symmetric, eastward and westward propagating unstable modes. The introduction of a zonal moisture gradient breaks the east-west symmetry of the unstable modes. The effects of rotation are then included by numerically solving the resulting eigenvalue problem on an equatorial beta-plane. With a purely meridional moisture gradient, the system supports large-scale, low-frequency, eastward and westward moving neutral modes. Some of the similarities, and some of the discrepancies of these modes with intraseasonal tropical waves are pointed out. Finally, a zonal moisture gradient in the presence of rotation renders some of the aforementioned neutral modes unstable. In particular, according to observations of large-scale, low-frequency tropical variability, it is seen that regions where the background saturation profile increases (decreases) to the east favour eastward (westward) moving moist modes.

Effect of Groundwater Seepage on Uplift Resistance of Buried Pipelines

The uplift resistance of pipelines buried in sands, in the presence of inclined groundwater flow, considering both upward and downward flow directions, has been determined by using the lower bound finite elements limit analysis in conjunction with nonlinear optimization. A correction factor (f (gamma) ), which needs to be multiplied with the uplift factor (F (gamma) ), has been computed to account for groundwater seepage. The variation of f (gamma) has been obtained as a function of i(gamma (w) /gamma (sub) ) for different horizontal inclinations (theta) of groundwater flow; where i = absolute magnitude of hydraulic gradient along the direction of flow, gamma (w) is the unit weight of water and gamma (sub) is the submerged unit weight of soil mass. For a given magnitude of i, there exists a certain critical value of theta for which the magnitude of f (gamma) becomes the minimum. An example has also been presented to illustrate the application of the results obtained for designing pipelines in presence of groundwater seepage.

Seismic Bearing Capacity of Shallow Embedded Foundations on a Sloping Ground Surface

By using the lower-bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization, bearing-capacity factors, N-c and N-gamma q, with an inclusion of pseudostatic horizontal seismic body forces, have been determined for a shallow embedded horizontal strip footing placed on sloping ground surface. The variation of N-c and N-gamma q with changes in slope angle (beta) for different values of seismic acceleration coefficient (k(h)) has been obtained. The analysis reveals that irrespective of ground inclination and the embedment depth of the footing, the factors N-c and N-gamma q decrease quite considerably with an increase in k(h). As compared with N-c, the factor N-gamma q is affected more extensively with changes in k(h) and beta. Unlike most of the results reported in literature for the seismic case, the present computational results take into account the shear resistance of soil mass above the footing level. An increase in the depth of the embedment leads to an increase in the magnitudes of both N-c and N-gamma q. (C) 2014 American Society of Civil Engineers.

Use of permeable reactive barrier to mitigate groundwater nitrate contamination from on-site sanitation

Nearly 50% of India's population depends on variants of pit-toilet systems for human waste disposal. Nitrate contamination of groundwater by pit-toilet leachate is a major environmental concern in the country as it sources a major proportion (50-80%) of potable water from aquifers. Therefore, minimizing nitrate contamination of groundwater due to leachate infiltration from pit-toilet systems is essential. Batch and column experiments demonstrated the capability of bentonite-enhanced sand (BES) specimens to reduce nitrate concentrations in synthetic solutions (initial NO3-N concentration = 22.7 mg/L, C/N = 3) by about 85-90% in 10 to 24 hour by a heterotrophic denitrification process. Based on the laboratory results, it is recommended that use of a BES-permeable reactive barrier layer at the base of pit-toilets will facilitate heterotrophic denitrification and mitigate nitrate contamination of the underlying aquifer.

Multiobjective fuzzy optimization for sustainable groundwater management using particle swarm optimization and analytic element method

Groundwater management involves conflicting objectives as maximization of discharge contradicts the criteria of minimum pumping cost and minimum piping cost. In addition, available data contains uncertainties such as market fluctuations, variations in water levels of wells and variations of ground water policies. A fuzzy model is to be evolved to tackle the uncertainties, and a multiobjective optimization is to be conducted to simultaneously satisfy the contradicting objectives. Towards this end, a multiobjective fuzzy optimization model is evolved. To get at the upper and lower bounds of the individual objectives, particle Swarm optimization (PSO) is adopted. The analytic element method (AEM) is employed to obtain the operating potentio metric head. In this study, a multiobjective fuzzy optimization model considering three conflicting objectives is developed using PSO and AEM methods for obtaining a sustainable groundwater management policy. The developed model is applied to a case study, and it is demonstrated that the compromise solution satisfies all the objectives with adequate levels of satisfaction. Sensitivity analysis is carried out by varying the parameters, and it is shown that the effect of any such variation is quite significant. Copyright (c) 2015 John Wiley & Sons, Ltd.