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.

Hydraulic Subsurface Measurements and Hydrodynamic Modeling as Indicators for Groundwater Flow Systems in the Rotondo Granite, Central Alps (Switzerland)

Regional groundwater flow in high mountainous terrain is governed by a multitude of factors such as geology, topography, recharge conditions, structural elements such as fracturation and regional fault zones as well as man-made underground structures. By means of a numerical groundwater flow model, we consider the impact of deep underground tunnels and of an idealized major fault zone on the groundwater flow systems within the fractured Rotondo granite. The position of the free groundwater table as response to the above subsurface structures and, in particular, with regard to the influence of spatial distributed groundwater recharge rates is addressed. The model results show significant unsaturated zones below the mountain ridges in the study area with a thickness of up to several hundred metres. The subsurface galleries are shown to have a strong effect on the head distribution in the model domain, causing locally a reversal of natural head gradients. With respect to the position of the catchment areas to the tunnel and the corresponding type of recharge source for the tunnel inflows (i.e. glaciers or recent precipitation), as well as water table elevation, the influence of spatial distributed recharge rates is compared to uniform recharge rates. Water table elevations below the well exposed high-relief mountain ridges are observed to be more sensitive to changes in groundwater recharge rates and permeability than below ridges with less topographic relief. In the conceptual framework of the numerical simulations, the model fault zone has less influence on the groundwater table position, but more importantly acts as fast flow path for recharge from glaciated areas towards the subsurface galleries. This is in agreement with a previous study, where the imprint of glacial recharge was observed in the environmental isotope composition of groundwater sampled in the subsurface galleries. Copyright © 2012 John Wiley & Sons, Ltd.

Detection of Groundwater flow zones in low productivity boreholes (Submitted)

Conventional methods of detecting groundwater flow zones in open boreholes installed in fractured bedrock aquifers rely on either contrasts in water quality or on significant rates of vertical flow. In low productivity boreholes these methods have proced of limited value. Tracer tests completed in six low productivity bedrock boreholes installed into fresh Precambrian crystalline bedrock revealed measurable dilution, despite very low tranmissivities. Fluorescent tracer profiles generated during discharge pumping permitted identification of the principle zones contributing to flow. High resolution acoustic televiewer strike and dip measurements of fractures corresponding to these zones revealed a consistent pattern of regional lineament trends and suggested a strongly anisotropic flow pattern in bedrock.

Radon emissions related to the granitic Precambrian shield in southern Brazil

The equivalent uranium (eU) activity concentration was analysed in selected granite samples at several sites in Porto Alegre, Southern Brazil, to obtain information on the radon (222Rn) generation by the aquifer rock matrices. Radon analyses of ground water and soil samples were also performed. Several samples exhibited a dissolved 222Rn activity concentration exceeding the World Health Organization maximum limit of 100 Bq l-1. The dissolved radon content in ground waters from the Fractured Precambrian Aquifer System exhibited a direct significant correlation with the eU in the rock matrices, which is a typical result of water-rock interactions. Variation in the soil's porosity was confirmed as an important factor for 222Rn release, as expected, due to its gaseous nature. Thus, although the calcic-alkaline to alkaline Precambrian granitoid rocks of the study area are important reservoirs for underground resources, they can release high amounts of radon gas into the liquid phase. © 2013 Copyright Taylor and Francis Group, LLC.

Images of a groundwater flow system in Nicaragua

Conceptualization of groundwater flow systems is necessary for water resources planning. Geophysical, hydrochemical and isotopic characterization methods were used to investigate the groundwater flow system of a multi-layer fractured sedimentary aquifer along the coastline in Southwestern Nicaragua. A geologic survey was performed along the 46 km2 catchment. Electrical resistivity tomography (ERT) was applied along a 4.4 km transect parallel to the main river channel to identify fractures and determine aquifer geometry. Additionally, three cross sections in the lower catchment and two in hillslopes of the upper part of the catchment were surveyed using ERT. Stable water isotopes, chloride and silica were analyzed for springs, river, wells and piezometers samples during the dry and wet season of 2012. Indication of moisture recycling was found although the identification of the source areas needs further investigation. The upper-middle catchment area is formed by fractured shale/limestone on top of compact sandstone. The lower catchment area is comprised of an alluvial unit of about 15 m thickness overlaying a fractured shale unit. Two major groundwater flow systems were identified: one deep in the shale unit, recharged in the upper-middle catchment area; and one shallow, flowing in the alluvium unit and recharged locally in the lower catchment area. Recharged precipitation displaces older groundwater along the catchment, in a piston flow mechanism. Geophysical methods in combination with hydrochemical and isotopic tracers provide information over different scales and resolutions, which allow an integrated analysis of groundwater flow systems. This approach provides integrated surface and subsurface information where remoteness, accessibility, and costs prohibit installation of groundwater monitoring networks.

Modelling solute and microorganism transport through a peri-alpine gravel aquifer by considering the aquifer as a dual porosity fractured medium.

Approach:
In-situ passive gradient comparative artificial tracer testing, undertaken using solutes (Uranine and Iodide), Bacteria (E.coli and P.putida) and bacteriophage (H40/1), permitted comparison of the mobility of different sized microorganisms relative to solutes in the sand and gravel aquifer underlying Dornach, Germany.
Tracer breakthrough curves reveal that even though uranine initially arrived at observation wells at the same time as microbiological tracers, maximum relative concentrations were sometimes less than those of microbiological tracers, while solute breakthrough curves proved more disperse.
Monitoring uranine breakthrough with depth suggested tracers arrived in observation wells in discrete 0.5m-1m thick intervals, over the aquifer’s 12m saturated thickness. Nearby exposures of aquifer material suggested that the aquifer consisted of sandy gravels enveloping sequences of open framework (OW) gravel up to 1m thick. Detailed examination of OW units revealed that they contained lenses of silty sand up to 1m long x 30cm thick., while granulometric data suggested that the gravel was two to three orders of magnitude more permeable than the enveloping sandy gravel.
Solute and microorganism tracer responses could not be simulated using conventional advective-dispersive equation solutions employing the same velocity and dispersion terms. By contrast solute tracer responses, modelled using a dual porosity approach for fractured media (DP-1D) corresponded well to observed field data. Simulating microorganism responses using the same transport terms, but no dual porosity term, generated good model fits and explained the higher relative concentration of the bacteria, compared to the non-reactive solute, even with first order removal to account for lower RR. Geologically, model results indicate that the silty units within open framework gravels are accessible to solute tracers, but not to microorganisms.
Importance:
Results highlight the benefits of geological observations developing appropriate conceptual models of solute and micro organism transport and in developing suitable numerical approaches to quantifying microorganism mobility at scales appropriate for the development of groundwater supply (wellhead) protection zones.

Numerical assessment of potential impacts of hydraulically fractured Bowland Shale on overlying aquifers

Natural gas extracted from hydraulically fractured shale formations potentially has a big impact on the global energy landscape. However, there are concerns of potential environmental impacts of hydraulic fracturing of the shale formations, particularly those related to water quality. To evaluate the potential impact of hydraulically fractured shale on overlying aquifers, we conduct realizations of numerical modeling simulations to assess fluid flow and chloride transport from a synthetic Bowland Shale over a period of 11,000 years. The synthetic fractured shale was represented by a three-dimensional discrete fracture model that was developed by using the data from a Bowland Shale gas exploration in Lancashire, UK. Chloride mass exchange between fractures and the rock matrix was fully accounted for in the model. The assessment was carried out to investigate fluid and chloride mass fluxes before, during, and after hydraulic fracturing of the Bowland Shale. Impacts of the upward fracture height and aperture, as well as hydraulic conductivity of the multilayered bedrock system, are also included this assessment. This modeling revealed that the hydraulically fractured Bowland Shale is unlikely to pose a risk to its overlying groundwater quality when the induced fracture aperture is ≤200 µm. With the fracture aperture ≥1000 µm, the upward chloride flux becomes very sensitive to the upward fracture height growth and hydraulic conductivity of the multilayered bedrock system. In the extremely unlikely event of the upward fracture growth directly connecting the shale formation to the overlying Sherwood Sandstone aquifer with the fracture aperture ≥1000 µm, the upward chloride mass flux could potentially pose risks to the overlying aquifer in 100 years. The model study also revealed that the upward mass flux is significantly intercepted by the horizontal mass flux within a high permeable layer between the Bowland Shale and its overlying aquifers, reducing further upward flux toward the overlying aquifers.

Drinking-water quality and variations in water levels in the fractured crystalline-rock aquifer, west-central Jefferson County, Colorado /

Mode of access: Internet.

Visualisation modelling in 3D of an alluvial aquifer system at a valley-wide scale to understand groundwater status : Lockyer, QLD

Visualisation provides a method to efficiently convey and understand the complex nature and processes of groundwater systems. This technique has been applied to the Lockyer Valley to aid in comprehending the current condition of the system. The Lockyer Valley in southeast Queensland hosts intensive irrigated agriculture sourcing groundwater from alluvial aquifers. The valley is around 3000 km2 in area and the alluvial deposits are typically 1-3 km wide and to 20-35 m deep in the main channels, reducing in size in subcatchments. The configuration of the alluvium is of a series of elongate “fingers”. In this roughly circular valley recharge to the alluvial aquifers is largely from seasonal storm events, on the surrounding ranges. The ranges are overlain by basaltic aquifers of Tertiary age, which overall are quite transmissive. Both runoff from these ranges and infiltration into the basalts provided ephemeral flow to the streams of the valley. Throughout the valley there are over 5,000 bores extracting alluvial groundwater, plus lesser numbers extracting from underlying sandstone bedrock. Although there are approximately 2500 monitoring bores, the only regularly monitored area is the formally declared management zone in the lower one third. This zone has a calibrated Modflow model (Durick and Bleakly, 2000); a broader valley Modflow model was developed in 2002 (KBR), but did not have extensive extraction data for detailed calibration. Another Modflow model focused on a central area river confluence (Wilson, 2005) with some local production data and pumping test results. A recent subcatchment simulation model incorporates a network of bores with short-period automated hydrographic measurements (Dvoracek and Cox, 2008). The above simulation models were all based on conceptual hydrogeological models of differing scale and detail.

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.

Geophysical and hydrogeological characterisation of the impacts of on-site wastewater treatment discharge to groundwater in a poorly productive bedrock acquifer

Contaminants discharging from on-site wastewater treatment systems (OSWTSs) can impact groundwater quality, threatening human health and surface water ecosystems. Risk of negative impacts becomes elevated in areas of extreme vulnerability with high water tables, where thin unsaturated intervals limit vadose zone attenuation. A combined geophysical/hydrogeological investigation into the effects of an OSWTS, located over a poorly productive aquifer (PPA) with thin subsoil cover, aimed to characterise effluent impacts on groundwater. Groundwater, sampled from piezometers down-gradient of the OSWTS percolation area displayed spatially erratic, yet temporally consistent, contaminant distributions. Electrical resistivity tomography identified an area of gross groundwater contamination close to the percolation area and, when combined with seismic refraction and water quality data, indicated that infiltrating effluent reaching the water table discharged to a deeper more permeable zone of weathered shale resting on more competent bedrock. Subsurface structure, defined by geophysics, indicated that elevated chemical and microbiological contaminant levels encountered in groundwater samples collected from piezometers, down-gradient of sampling points with lower contaminant levels, corresponded to those locations where piezometers were screened close to the weathered shale/competent rock interface; those immediately up-gradient were too shallow to intercept this interval, and thus the more impacted zone of the contaminant plume. Intermittent occurrence of faecal indicator bacteria more than 100 m down gradient of the percolation area suggested relatively short travel times. Study findings highlight the utility of geophysics as part of multidisciplinary investigations for OSWTS contaminant plume characterisation, while also demonstrating the capacity of effluent discharging to PPAs to impact groundwater quality at distance. Comparable geophysical responses observed in similar settings across Ireland suggest the phenomena observed in this study are more widespread than previously suspected.

Analysis of groundwater-level response to rainfall and recharge estimates in fractured hard rock aquifers, NW Ireland

Despite fractured hard rock aquifers underlying over 65% of Ireland, knowledge of key processes controlling groundwater recharge in these bedrock systems is inadequately constrained. In this study, we examined 19 groundwater-level hydrographs from two Irish hillslope sites underlain by hard rock aquifers. Water-level time-series in clustered monitoring wells completed at the subsoil, soil/bedrock interface, shallow and deep bedrocks were continuously monitored hourly over two hydrological years. Correlation methods were applied to investigate groundwater-level response to rainfall, as well as its seasonal variations. The results reveal that the direct groundwater recharge to the shallow and deep bedrocks on hillslope is very limited. Water-level variations within these geological units are likely dominated by slow flow rock matrix storage. The rapid responses to rainfall (⩽2 h) with little seasonal variations were observed to the monitoring wells installed at the subsoil and soil/bedrock interface, as well as those in the shallow or deep bedrocks at the base of the hillslope. This suggests that the direct recharge takes place within these units. An automated time-series procedure using the water-table fluctuation method was developed to estimate groundwater recharge from the water-level and rainfall data. Results show the annual recharge rates of 42–197 mm/yr in the subsoil and soil/bedrock interface, which represent 4–19% of the annual rainfall. Statistical analysis of the relationship between the rainfall intensity and water-table rise reveal that the low rainfall intensity group (⩽1 mm/h) has greater impact on the groundwater recharge rate than other groups (>1 mm/h). This study shows that the combination of the time-series analysis and the water-table fluctuation method could be an useful approach to investigate groundwater recharge in fractured hard rock aquifers in Ireland.

Delineating environmental groundwater vulnerability and protection zones mapping in fractured rock masses

Hard‐rock watersheds commonly exhibit complex geological bedrock and morphological features. Hydromineral resources have relevant economic value for the thermal spas industry. The present study aims to develop a groundwater vulnerability approach in Caldas da Cavaca hydromineral system (Aguiar da Beira, Central Portugal) which has a thermal tradition that dates back to the late 19th century, and contribute to a better understanding of the hydrogeological conceptual site model. In this work different layers were overlaid, generating several thematic maps to arrive at an integrated framework of several key‐sectors in Caldas da Cavaca site. Thus, to accomplish a comprehensive analysis and conceptualization of the site, a multi‐technical approach was used, such as, field and laboratory techniques, where several data was collected, like geotectonics, hydrology and hydrogeology, hydrogeomorphology, hydrogeophysical and hydrogeomechanical zoning aiming the application of the so‐called DISCO method. All these techniques were successfully performed and a groundwater vulnerability to contamination assessment, based on GOD‐S, DRASTIC‐Fm, SINTACS, SI and DISCO indexes methodology, was delineated. Geographical Information Systems (GIS) technology was on the basis to organise and integrate the geodatabases and to produce all the thematic maps. This multi‐technical approach highlights the importance of groundwater vulnerability to contamination mapping as a tool to support hydrogeological conceptualisation, contributing to better decision‐making of water resources management and sustainability.

An Integrated Hydro geological Study of the Muvattupuzha River Basin, Kerala, India

The present investigation on the Muvattupuzha river basin is an integrated approach based on hydrogeological, geophysical, hydrogeochemical parameters and the results are interpreted using satellite data. GIS also been used to combine the various spatial and non-spatial data. The salient finding of the present study are accounted below to provide a holistic picture on the groundwaters of the Muvattupuzha river basin. In the Muvattupuzha river basin the groundwaters are drawn from the weathered and fractured zones. The groundwater level fluctuations of the basin from 1992 to 2001 reveal that the water level varies between a minimum of 0.003 m and a maximum of 3.45 m. The groundwater fluctuation is affected by rainfall. Various aquifer parameters like transmissivity, storage coefficient, optimum yield, time for full recovery and specific capacity indices are analyzed. The depth to the bedrock of the basin varies widely from 1.5 to 17 mbgl. A ground water prospective map of phreatic aquifer has been prepared based on thickness of the weathered zone and low resistivity values (<500 ohm-m) and accordingly the basin is classified in three phreatic potential zones as good, moderate and poor. The groundwater of the Muvattupuzha river basin, the pH value ranges from 5.5 to 8.1, in acidic nature. Hydrochemical facies diagram reveals that most of the samples in both the seasons fall in mixing and dissolution facies and a few in static and dynamic natures. Further study is needed on impact of dykes on the occurrence and movement of groundwater, impact of seapages from irrigation canals on the groundwater quality and resources of this basin, and influence of inter-basin transfer of surface water on groundwater.