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.

From Plinthic Acrisols to Plinthosols and Gleysols: iron and groundwater dynamics in the tertiary sediments of the upper Amazon basin

Topography has been reported to be the major factor ruling the spatial distribution of Acrisols, Plinthosols and Gleysols on the seasonally flooded, low elevation plateaux of the upper Amazon basin occupied by Tertiary (Ica & Solimoes) sediments. In this study, detailed morphological and mineralogical investigations conducted in a representative 25-ha site were combined with hydro-geochemical data to relate the vertical and lateral soil differentiations observed to the hydro-geological history of that part of the basin. As a result of the uplift of the Andes, several cuts in the extensive Tertiary marshlands have formed, at first, slightly incised plateaux of low elevation. There, weathering under hot and humid climates would have generated a reddish, freely drained and bioturbated topsoil layer and the vertical differentiation in subsoil sediments of a plinthite over an iron-depleted mottled clay. The second episode of soil differentiation is linked to the replacement of the forest by a savannah under the drier climates of the late Pleistocene, which favours surface runoff and the infill of the incisions by fine particles. This infill, combined with the return to the present humid climate, has then enabled the local groundwater to rise on the plateaux and to generate episaturation at the topsoil/subsoil transition close to the depressions. Nowadays, ferrous iron is released from the partly iron-depleted topsoil weathering front at high water levels during the rainy seasons. It moves from footslope to low-lying positions and from top to bottom in the soil profile according to the groundwater dynamics. The present general trend is thus to the lateral export of iron at high water levels due to subsurface and overland flows, its vertical transfer during the recession of the groundwater and accumulation in a nodular plinthite. In the latter, ferrous iron is adsorbed onto its softest iron masses where it feeds the neoformation of ferrihydrite that rapidly dehydrates into haematite.

Spatial organisation of groundwater dynamics and streamflow response from different hydropedological units in a montane catchment

Funding was provided by the Leibniz Association (SAW-2012-IGB 4167) within the International Leibniz Graduate School: Aquatic boundaries and linkages- Aqualink. We would like to thank the NRI staff for their help during field work.

The Distribution and Quality of Hine's Emerald Dragonfly Habitat in Relation to Surface and Groundwater Dynamics in the Lockport Prairie Preserve Preliminary Report (February 1, 2001-January 31, 2002)

INHS Technical Report prepared for Corporation for Openlands, US Fish and Wildlife Service

Numerical simulation of a new two-dimensional variable-order fractional percolation equation in non-homogeneous porous media

Percolation flow problems are discussed in many research fields, such as seepage hydraulics, groundwater hydraulics, groundwater dynamics and fluid dynamics in porous media. Many physical processes appear to exhibit fractional-order behavior that may vary with time, or space, or space and time. The theory of pseudodifferential operators and equations has been used to deal with this situation. In this paper we use a fractional Darcys law with variable order Riemann-Liouville fractional derivatives, this leads to a new variable-order fractional percolation equation. In this paper, a new two-dimensional variable-order fractional percolation equation is considered. A new implicit numerical method and an alternating direct method for the two-dimensional variable-order fractional model is proposed. Consistency, stability and convergence of the implicit finite difference method are established. Finally, some numerical examples are given. The numerical results demonstrate the effectiveness of the methods. This technique can be used to simulate a three-dimensional variable-order fractional percolation equation.

Chloride as a signature indicator of soil textural and hydrologic stratigraphies in variable charge deep profiles

Soil properties that influence water movement through profiles are important for determining flow paths, reactions between soil and solute, and the ultimate destination of solutes. This is particularly important in high rainfall environments. For highly weathered deep profiles, we hypothesize that abrupt changes in the distribution of the quotient [QT = (silt + sand)/clay] reflect the boundaries between textural units or textural (TS) and hydrologic (HS) stratigraphies. As a result, QT can be used as a parameter to characterize TS and as a surrogate for HS. Secondly, we propose that if chloride distributions were correlated with QT, under non-limiting anion exchange, then chloride distributions can be used as a signature indicator of TS and HS. Soil cores to a depth of 12.5 in were taken from 16 locations in the wet tropical Johnstone River catchment of northeast Queensland, Australia. The cores belong to nine variable charge soil types and were under sugarcane (Saccharun officinarum-S) production, which included the use of potassium chloride, for several decades. The cores were segmented at I m depth increments and subsamples were analysed for chloride, pH, soil water content (theta), clay, silt and sand contents. Selected bores were capped to serve as piezometers to monitor groundwater dynamics. Depth incremented QT, theta and chloride correlated, each individually, significantly with the corresponding profile depth increments, indicating the presence of textural, hydrologic and chloride gradients in profiles. However, rapid increases in QT down the profile indicated abrupt changes in TS, suggesting that QT can be used as a parameter to characterize TS and as a surrogate for HS. Abrupt changes in chloride distributions were similar to QT, suggesting that chloride distributions can be used as a signature indicator of QT (TS) and HS. Groundwater data indicated that chloride distributions depended, at least partially, on groundwater dynamics, providing further support to our hypothesis that chloride distribution can be used as a signature indicator of HS. Copyright (c) 2005 John Wiley & Sons, Ltd.

Characterisation of sea-water intrusion in the Pioneer Valley, Australia using hydrochemistry and three-dimensional numerical modelling

Sea-water intrusion is actively contaminating fresh groundwater reserves in the coastal aquifers of the Pioneer Valley,north-eastern Australia. A three-dimensional sea-water intrusion model has been developed using the MODHMS code to explore regional-scale processes and to aid assessment of management strategies for the system. A sea-water intrusion potential map, produced through analyses of the hydrochemistry, hydrology and hydrogeology, offsets model limitations by providing an alternative appraisal of susceptibility. Sea-water intrusion in the Pioneer Valley is not in equilibrium, and a potential exists for further landward shifts in the extent of saline groundwater. The model required consideration of tidal over-height (the additional hydraulic head at the coast produced by the action of tides), with over-height values in the range 0.5-0.9 m giving improved water-table predictions. The effect of the initial water-table condition dominated the sensitivity of the model to changes in the coastal hydraulic boundary condition. Several salination processes are probably occurring in the Pioneer Valley, rather than just simple landward sea-water advancement from modern sources of marine salts. The method of vertical discretisation (i.e. model-layer subdivision) was shown to introduce some errors in the prediction of watertable behaviour.

Geophysical characterization of the complex dynamics of groundwater and seawater exchange in a highly stressed aquifer system linked to a coastal lagoon (SE Spain)

Anthropogenic pressure influences the two-way interactions between shallow aquifers and coastal lagoons. Aquifer overexploitation may lead to seawater intrusion, and aquifer recharge from rainfall plus irrigation may, in turn, increase the groundwater discharge into the lagoon. We analyse the evolution, since the 1950s up to the present, of the interactions between the Campo de Cartagena Quaternary aquifer and the Mar Menor coastal lagoon (SE Spain). This is a very heterogeneous and anisotropic detrital aquifer, where aquifer–lagoon interface has a very irregular geometry. Using electrical resistivity tomography, we clearly identified the freshwater–saltwater transition zone and detected areas affected by seawater intrusion. Severity of the intrusion was spatially variable and significantly related to the density of irrigation wells in 1950s–1960s, suggesting the role of groundwater overexploitation. We distinguish two different mechanisms by which water from the sea invades the land: (a) horizontal advance of the interface due to a wide exploitation area and (b) vertical rise (upconing) caused by local intensive pumping. In general, shallow parts of the geophysical profiles show higher electrical resistivity associated with freshwater mainly coming from irrigation return flows, with water resources mostly from deep confined aquifers and imported from Tagus river, 400 km north. This indicates a likely reversal of the former seawater intrusion process.

Reading Our Land Use Legacy in Groundwater: Nitrate Dynamics in the Judith River Watershed

Elevated nitrate in groundwater is common is agricultural areas where fertilizer has been added at high rates for decades. Within the Judith River Wastershed, high native soil fertility allowed for dryland wheat production without N fertilization until the 1980s, yet elevated nitrate levels were frequently observed in shallow aquifers. Dr. Stephanie Ewing presents results for soil, groundwater and surface water analyses from a hydrologically isolated strath terrace near Moccasin, MT. In context of this uniquely well constrained field setting, these observed data, along with land use history and a simple mass balance model, revel the long term development and perturbation of native soil fertility with cultivation.

Baseflow dynamics : multi-tracer surveys to assess variable groundwater contributions to montane streams under low flows

Peer reviewed

Habitat dynamics in the bed sediments of an intermittent upland stream

The effects of the hydrological regime on temporal changes to physical characteristics of substratum habitat, sediment texture of surface sediments (<10 cm), were investigated in a sub-tropical headwater stream over four years. Surface discharge was measured together with vertical hydraulic gradient and groundwater depth in order to explore features of sediment habitat that extend beyond the streambed surface. Whilst the typical discharge pattern was one of intermittent base flows and infrequent flow events associated with monsoonal rain patterns, the study period also encompassed a drought and a one-in-a-hundred-year flood. Rainfall and discharge did not necessarily reflect the actual conditions in the stream. Although surface waters were persistent long after discharge ceased, the streambed was completely dry on several occasions. Shallow groundwater was present at variable depths throughout the study period, being absent only at the height of the drought. The streambed sediments were mainly gravels, sand and clay. Finer sediment fractions showed a marked change in grain size over time, although bedload movement was limited to a single high discharge event. In response to a low discharge regimen (drought), sediments characteristically showed non-normal distributions and were dominated by finer materials. A high-energy discharge event produced a coarsening of sands and a diminished clay fraction in the streambed. Particulate organic matter from sediments showed trends of build-up and decline with the high and low discharge regimes, respectively. Within the surface sediment intersticies three potential categories of invertebrate habitat were recognised, each with dynamic spatial and temporal boundaries.

Studies of Earth Dynamics with the Superconducting Gravimeter

The superconducting (or cryogenic) gravimeter (SG) is based on the levitation of a super­conducting sphere in a stable magnetic field created by current in superconducting coils. Depending on frequency, it is capable of detecting gravity variations as small as 10-11ms-2. For a single event, the detection threshold is higher, conservatively about 10-9 ms-2. Due to its high sensitivity and low drift rate, the SG is eminently suitable for the study of geodynamical phenomena through their gravity signatures. I present investigations of Earth dynamics with the superconducting gravimeter GWR T020 at Metsähovi from 1994 to 2005. The history and key technical details of the installation are given. The data processing methods and the development of the local tidal model at Metsähovi are presented. The T020 is a part of the worldwide GGP (Global Geodynamics Project) network, which consist of 20 working station. The data of the T020 and of other participating SGs are available to the scientific community. The SG T020 have used as a long-period seismometer to study microseismicity and the Earth s free oscillation. The annual variation, spectral distribution, amplitude and the sources of microseism at Metsähovi were presented. Free oscillations excited by three large earthquakes were analyzed: the spectra, attenuation and rotational splitting of the modes. The lowest modes of all different oscillation types are studied, i.e. the radial mode 0S0, the "football mode" 0S2, and the toroidal mode 0T2. The very low level (0.01 nms-1) incessant excitation of the Earth s free oscillation was detected with the T020. The recovery of global and regional variations in gravity with the SG requires the modelling of local gravity effects. The most important of them is hydrology. The variation in the groundwater level at Metsähovi as measured in a borehole in the fractured bedrock correlates significantly (0.79) with gravity. The influence of local precipitation, soil moisture and snow cover are detectable in the gravity record. The gravity effect of the variation in atmospheric mass and that of the non-tidal loading by the Baltic Sea were investigated together, as sea level and air pressure are correlated. Using Green s functions it was calculated that a 1 metre uniform layer of water in the Baltic Sea increases the gravity at Metsähovi by 31 nms-2 and the vertical deformation is -11 mm. The regression coefficient for sea level is 27 nms-2m-1, which is 87% of the uniform model. These studies are associated with temporal height variations using the GPS data of Metsähovi permanent station. Results of long time series at Metsähovi demonstrated high quality of data and correctly carried out offsets and drift corrections. The superconducting gravimeter T020 has been proved to be an eminent and versatile tool in studies of the Earth dynamics.