Hydrogeology of the Brunswick (Passaic) Formation and Implications for Ground Water Monitoring Practice

Fractured shales of the Brunswick Formation provide a major aquifer in the most industrialized region of New Jersey. Numerous cases of ground water contamination have been documented in this formation. However, effectiveness of monitoring and remediation efforts is often hampered by the use of inappropriate concepts regarding ground water flow controls in this complex aquifer system. One such concept presumes that near-vertical fractures parallel to the strike of beds provide principal passages for the flow and produce an anisotropic response to pumping stress. Field evidence presented in this paper confirms that the Brunswick Formation hosts a gently dipping, multiunit, leaky aquifer system that consists of thin water-bearing units and thick intervening aquitards. The water-bearing units are associated with major bedding partings and/or intensely fractured seams. Layered heterogeneity of such a dipping multiunit aquifer system produces an anisotropic flow pattern with preferential flow along the strike of beds. Within the weathered zone, the permeability of the water-bearing units can be greatly reduced. The commonly used hydrogeologic model of the Brunswick as a one-aquifer system, sometimes with vaguely defined "shallow" and "deep" zones, often leads to the development of inadvertent cross-flows within monitoring wells. If undetected, cross-flows may promote contaminant spread into deeper units and impair the quality of hydrogeologic data. Hydrogeologic characterization of the Brunswick shales at any given site should be aimed primarily at identification of the major water-bearing and aquitard units. Recommended techniques for this characterization include fluid logging and other in-well tests.     

Resolution of MRS applied to the characterization of hard-rock aquifers

The performance of the Magnetic Resonance Sounding (MRS) method applied to the investigation of heterogeneous hard-rock aquifers was studied. It was shown using both numerical modeling and field measurements that MRS could be applied to the investigation of the weathered part of hard-rock aquifers when the product of the free water content multiplied by the thickness of the aquifer is >0.2 (for example, 10-m-thick layer with a 2% water content). Using a currently available one-dimensional MRS system, the method allows the characterization of two-dimensional subsurface structures with acceptable accuracy when the size of the subsurface anomaly is equal to or greater than the MRS loop. However, the fractured part of hard-rock aquifers characterized by low effective porosity (<0.5%) cannot be resolved using currently available MRS equipment. It was found that shallow water in the weathered part of the aquifer may screen MRS signals from deeper water-saturated layers, thus further reducing the possibility of investigating deeper fractured aquifers. A field study using the NUMIS(plus) MRS system developed by IRIS Instruments was carried out on an experimental watershed in southern India. A heterogeneous unconfined aquifer in a gneissic formation was successfully localized, and MRS results were confirmed by drilling shortly after the geophysical study. The top of the aquifer revealed by MRS was found to be in a good agreement with observed static water level measurements in boreholes.     

Hydraulics of a partially penetrating well with skin zone in a confined aquifer

A steady/quasi-steady model is developed for predicting flow into a partially penetrating well with skin zone in a confined aquifer overlying an impervious layer. The model takes into account flow through the bottom of the wellbore, finite skin thickness and finite horizontal and vertical extent of the aquifer. Moreover, the solution can be easily extended to include the mixed-type boundary condition at the well face, where a Dirichlet in the form of a specified hydraulic head and a Neumann in the form of zero flux coexist at the same time at different portions of the well face. The validity of the proposed solution is tested by comparing a few results obtained from the developed model with corresponding results obtained by analytical and numerical means. The study shows that, among other factors remaining constant, both the horizontal and vertical extent of an artesian aquifer, thickness of the skin zone, bottom flow and conductivity contrast of the skin and formation zones, play an important part in deciding flow to a well dug in the aquifer, and hence these factors must be considered while analyzing the problem. The model proposed here can be used to estimate skin thickness as well as hydraulic conductivities of the skin and formation zones of a well with skin zone in an artesian aquifer underlain by an impervious layer by utilizing pumping test data falling in the steady or quasi-steady state of a typical pumping test. As the proposed solution is of a general nature in the sense that it can handle, apart from partial penetration and bottom flow, the finite size skin zone and finite horizontal and vertical extent of an artesian aquifer together with the mixed-type boundary condition at the well face, it is hoped that the predictions coming out of the model will be more realistic than those obtained using solutions developed with more stringent assumptions.     

A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

Using temperature to test models of flow near Yucca Mountain,Nevada

Ground water temperatures in the fractured volcanic aquifer near Yucca Mountain, Nevada, have previously been shown to have significant spatial variability with regions of elevated temperatures coinciding roughly with near-vertical north-south trending faults. Using insights gained from one-dimensional models, previous investigators have suggested upwelling along faults from an underlying aquifer as a likely explanation for this ground water temperature pattern. Using a three-dimensional coupled flow and heat-transport model, we show that the thermal high coinciding with the Paintbrush fault zone can be explained without significant upwelling from the underlying aquifer. Instead, the thermal anomaly is consistent with thermal conduction enhanced slightly by vertical ground water movement within the volcanic aquifer sequence. If more than approximately 400 m3/day of water enters the volcanic aquifer from below along a 10 km fault zone, the calculated temperatures at the water table are significantly greater than the measured temperatures. These results illustrate the potential limitations in using one-dimensional models to interpret ground water temperature data, and underscore the value in combining temperature data with fully coupled three-dimensional simulations.     

Evaluation of an over-used coastal aquifer (Orissa,India) using statistical approaches

Abstract

This study evaluates an over-exploited aquifer (Balasor, India) and also explores the possibilities of sustainable management using several statistical approaches. First, bootstrap analysis indicates that groundwater pumping has resulted in the reduction of mean cultivated area as the average irrigation capacity per bore well dropped from 3.74 ha to 1.5 ha within a period of 10 years of operation. However, modelling the groundwater levels using the seasonal autoregressive integrated moving average (SARIMA) procedure showed no evidence of large-scale groundwater withdrawals. The derived models can be used for water table forecasting and also for infilling the missing observations. The empirical relationship among pumping test results indicated that well depth and aquifer thickness significantly influence the discharge from the aquifer. This discharge may have encouraged the resource-rich farmers to exploit the lower aquifer. Based on a contour plot, the zone of groundwater exploitation was estimated to vary from 25 to 60 m below the surface. Therefore, a uniform aquifer exploitation policy needs to be implemented to curb the vertical competition in exploiting the aquifer and to develop sustainable management of the groundwater.

Citation Panda, D. K. & Kumar, A. (2011) Evaluation of an over-used coastal aquifer (Orissa, India) using statistical approaches. Hydrol. Sci. J. 56(3), 486–497.     

Heat Tracing in a Fractured Aquifer with Injection of Hot and Cold Water

Heat as a tracer in fractured porous aquifers is more sensitive to fracture-matrix processes than a solute tracer. Temperature evolution as a function of time can be used to differentiate fracture and matrix characteristics. Experimental hot (50 °C) and cold (10 °C) water injections were performed in a weathered and fractured granite aquifer where the natural background temperature is 30 °C. The tailing of the hot and cold breakthrough curves, observed under different hydraulic conditions, was characterized in a log–log plot of time vs. normalized temperature difference, also converted to a residence time distribution (normalized). Dimensionless tail slopes close to 1.5 were observed for hot and cold breakthrough curves, compared to solute tracer tests showing slopes between 2 and 3. This stronger thermal diffusive behavior is explained by heat conduction. Using a process-based numerical model, the impact of heat conduction toward and from the porous rock matrix on groundwater heat transport was explored. Fracture aperture was adjusted depending on the actual hydraulic conditions. Water density and viscosity were considered temperature dependent. The model simulated the increase or reduction of the energy level in the fracture-matrix system and satisfactorily reproduced breakthrough curves tail slopes. This study shows the feasibility and utility of cold water tracer tests in hot fractured aquifers to boost and characterize the thermal matrix diffusion from the matrix toward the flowing groundwater in the fractures. This can be used as complementary information to solute tracer tests that are largely influenced by strong advection in the fractures.     

ELECTROMAGNETIC PROSPECTING FOR GROUNDWATER IN PRECAMBRIAN TERRAINS IN THE REPUBLIC OF UPPER VOLTA*

The “Autorité des amenegements des valées des Voltas (AVV)” is establishing new rural settlements in the Volta valleys. First, a survey of available water supplies is performed. Economic aquifers in Precambrian terrains are deep (15–50 m) and usually occur in fractured zones accompanying faults. Such zones can be identified on aerial photographs, but their precise location on the ground is virtually impossible by visual means. Because of the small size of the aquifers, a location error of 5 m can make the difference between a productive well and a dry hole. Traditionally, resistivity profiling has been used as the means of locating the fractured zones in the field. Our studies suggest that the task can be performed faster, cheaper and more accurately by VLF and EM methods. Because of the limited choice of transmitting stations reccivable in Upper Volta, the VLF method is not sufficiently sensitive to detect conductors with a strike between 45° and 105°. The results obtained with a multifrequency, horizontal-loop EM (HLEM) system were satisfactory in all investigated areas. During the 1980 field season, 35 target areas were surveyed. Of the 24 holes drilled so far, 23 are productive. The weathered layer is a source of distinctive HLEM anomalies, which are characteristic of the underlying rocks. Therefore, different interpretational procedures had to be developed for granitic and volcano-sedimentary areas. Despite the high background level of in-phase and quadrature components, which varied with thickness and conductivity of the weathered layer, aquifers could be detected at a depth greater than 30 m. Attempts were made to interpret the HLEM results quantitatively using two models: a three-layer medium and a valley discontinuity. The latter model is more realistic, but more scale modelling will have to be performed to permit development of viable interpretational procedures. Meanwhile, phasor diagrams based on drilling and resistivity sounding data can be used to estimate the aquifer depth.     

Environmental tracer transport (H and SF6) in the saturated and unsaturated zones and its use in nitrate pollution management

An important quantity in groundwater protection is the residence time of water in an aquifer. It relates to both the travel time of a pollutant to arrive at a well and the time span required for self-purification of a polluted aquifer after removal of pollutant inputs. Time scales for aquifers can be gained from artificial tracer experiments or from environmental tracer data, the latter offering the only realistic alternative if time scales of years or decades have to be taken into account.

Different tracers show different time scales due to their different transport mechanisms especially in the unsaturated zone. While solute tracers are moved advectively with the seepage water, gas tracers pass the unsaturated zone diffusively through the air phase. Depending on the properties of the unsaturated zone (hydraulic properties, thickness) this difference in behavior can be used to separate the subsurface transport process into the unsaturated and the saturated parts.

In a field study in Germany, SF₆ and ³H were used as environmental tracers. Both have a relatively well-known input function. Interpretation of data from observation wells by a box model approach led to spatially and temporally varying residence times. This was an indication that the influence of the unsaturated zone could not be neglected. While the gas tracer SF₆ shows only residence times in the saturated zone, the tracer ³H reflects the whole travel time of water including both the unsaturated and saturated zones. Using a one-dimensional plug-flow model for the unsaturated zone combined with a detailed two-dimensional flow and transport model for the saturated zone leads to a holistic and consistent interpretation of the measured tracer concentrations. The observed pattern of old water under thick loess cover and younger water under areas where the fractured basalt aquifer crops out is reproduced after adjusting only two parameters: the effective porosity of the saturated aquifer and the product of field capacity and thickness of the unsaturated zone. While the effective porosity of the saturated zone is adjusted by means of the SF₆ data, the field capacity of the loess layer is adjusted by means of the ³H observations. The thickness of the unsaturated zone is deduced from geological and pedological maps. All flow data are obtained from a calibrated flow model, which is based on geological data, observed heads and pumping tests only.

The transport model for the saturated zone was calibrated by fitting the porosity by means of gaseous tracer concentrations (SF₆). The combined saturated–unsaturated zone model was then calibrated by fitting the field capacity of the unsaturated zone by means of ³H concentrations. With this model it was possible to verify the observed NO₃ concentrations at the drinking water wells and to develop predictions for their future development under various scenarios of fertilizer input reduction in specific areas.     

A Rapid Drawdown Event Related to Quarry Operations in a Dolomite Aquifer

Two active quarries are mining stone from the Silurian dolomite aquifer in Waukesha County in southeastern Wisconsin. The village in which the quarries are located uses local zoning to control the depth of mining and to institute a long-term water level monitoring program and well guarantee/one with the quarry owners. Water levels dropped as much as 40 feet in at least 24 residential wells surrounding the quarries over a period of a few hours to days. The rapid decline in head was caused by a single boring drilled lo a depth of 75 feel he low the floor of one quarry. The borehole penetrated a localized fracture zone under confined artesian head. Water levels recovered to previous static levels within nine days after grouting the borehole. The rapid drawdown event demonstrates the potential impact of mining in fractured aquifers. The apparent complete recovery of the aquifer demonstrates that quick response can sometimes restore an aquifer. However, the potential for blasting into a similar zone illustrates the need for a well-thought-out aquifer monitoring program and emergency response plan. The experience of the village is a good example of managing conflicting uses of a finite resource and collecting baseline data needed to make informed decisions.     

Field data and flow system response in clay (vertisol) shale terrain,north central Texas,USA

The water budget in clay shale terrain is controlled by a complex interaction between the vertisol soil layer, the underlying fractured rock, land use, topography, and seasonal trends in rainfall and evapotranspiration. Rainfall, runoff, lateral flow, soil moisture, and groundwater levels were monitored over an annual recharge cycle. Four phases of soil–aquifer response were noted over the study period: (1) dry‐season cracking of soils; (2) runoff initiation, lateral flow and aquifer recharge; (3) crack closure and down‐slope movement of subsurface water, with surface seepage; (4) a drying phase. Surface flow predominated within the watershed (25% of rainfall), but lateral flow through the soil zone continued for most of the year and contributed 11% of stream flow through surface seepage. Actual flow through the fractured shale makes up a small fraction of the water budget but does appear to influence surface seepage by its effect on valley‐bottom storage. When the valley soil storage is full, lateral flow exits onto the valley‐bottom surface as seasonal seeps. Well response varied with depth and hillslope position. FLOWTUBE model results and regional recharge estimates are consistent with an aquifer recharge of 1·6% of annual precipitation calculated from well heights and specific yield of the shale aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

Determining anisotropic transmissivity using a simplified Papadopulos method

The straight-line method presented by Papadopulos requires a minimum of three observation wells for determining the transmissivity tensor of a homogeneous and anisotropic aquifer. A simplification of this method was developed for fractured aquifers where the principal directions of the transmissivity tensor are known prior to implementation, such as when fracture patterns on outcropping portions of the aquifer may be used to infer the principal directions. This new method assumes that observation wells are drilled along the two principal directions from the pumped well, thus reducing the required number of observation wells to two. This method was applied for an aquifer test in the fractured Navajo Sandstone of southwestern Utah and yielded minimum and maximum principal transmissivity values of 70 and 1800 m(2)/d, respectively, indicating an anisotropy ratio of approximately 24 to 1.     

A Tracer Test for Detecting Cross Contamination Along a Monitoring Well Column

A tracer test was used to evaluate whether cross contamination exists along a monitoring well completed through a shallow ground water system in fractured clay and screened in a sand and gravel aquifer. The fractured clay is separated from the sand and gravel deposit by a layer of highly plastic unfractured clay. A natural vertical downward hydraulic gradient of approximately 0.5 exists between the shallow system and the sand and gravel aquifer. Ground water contamination was detected in an adjacent monitoring well screened in the fractured clay and in the monitoring well screened in the sand and gravel deposit. No ground water contamination was apparent in an intermediate well screened in the unfractured clay layer. A tracer of sodium bromide was injected into a shallow boring near the monitoring wells. The tracer was detected in the monitoring well in the sand and gravel aquifer after three to seven days. The bromide concentration continued to increase in this well with time while the concentration in the shallow boring declined. This trend of tracer concentration indicates the tracer has in fact migrated downward and possibly traveled along the well column.     

A New Randomized Binary Prior Model for Hydraulic Tomography in Fractured Aquifers

We present a novel pilot-point-based hydraulic tomography (HT) inversion procedure to delineate preferential flow paths and estimate hydraulic properties in a fractured aquifer. Our procedure considers a binary prior model developed using a randomized algorithm. The randomized algorithm involves discretizing the domain into grid cells, assigning a binary label to each cell, traversing the grid randomly, and choosing the optimal grid configuration cell-by-cell. This binary prior model is used to guide the placement of pilot points and to constrain aquifer parameters during pilot-point-based HT inversion. A two-dimensional fractured granite rock block was considered to test our methodology under controlled laboratory conditions. Multiple pumping tests were conducted at selected ports and the pressure responses were monitored. The pumping datasets thus obtained were preprocessed using median filters to remove random noise, and then analyzed using the proposed procedure. The proposed binary prior algorithm was implemented in C++ by supplying the forward groundwater model, HydroGeoSphere (HGS). Pilot-point-assisted HT inversion was performed using the parameter-estimation tool, coupled to HGS. The resulting parameter distributions were assessed by: (1) a visual comparison of the K- and S_s-tomograms with the known topology of the fractures and (2) comparing model predictions with measurements made at two validation ports that were not used in calibration. The performance assessment revealed that HT with the proposed randomized binary prior could be used to recover fracture-connectivity and to predict drawdowns in fractured aquifers with reasonable accuracy, when compared to a conventional pilot-point inversion scheme.     

Reverse Water‐Level Fluctuations Associated with Fracture Connectivity

Reverse water‐level fluctuations (RWFs), a phenomenon in which water levels rise briefly in response to pumping, were detected in monitoring wells in a fractured siliciclastic aquifer system near a deep public supply well. The magnitude and timing of RWFs provide important information that can help interpret aquifer hydraulics near pumping wells. A RWF in a well is normally attributed to poroelastic coupling between the solid and fluid components in an aquifer system. In addition to revealing classical pumping‐induced poroelastic RWFs, data from pressure transducers located at varying depths and distances from the public supply well suggest that the RWFs propagate rapidly through fractures to influence wells hundreds of meters from the pumping well. The rate and cycling frequency of pumping is an important factor in the magnitude of RWFs. The pattern of RWF propagation can be used to better define fracture connectivity in an aquifer system. Rapid, cyclic head changes due to RWFs may also serve as a mechanism for contaminant transport.     

Atoll island hydrogeology: flow and freshwater occurrence in a tidally dominated system

A layered-aquifer model of groundwater occurrence in an atoll island was tested with a solute-transport numerical model. The computer model used, SUTRA, incorporates density-dependent flow. This can be significant in freshwater-saltwater interactions associated with the freshwater lens of an atoll island. Boundary conditions for the model included ocean and lagoon tidal variations. The model was calibrated to field data from Enjebi Island, Enewetak Atoll, and tested for sensitivity to a variety of parameters. This resulted in a hydraulic conductivity of 10 m day⁻¹ for the surficial aquifer and 1000 m day⁻¹ for the deeper aquifer; this combination of values gave an excellent reproduction of the tidal response data from test wells. The average salinity distribution was closely reproduced using a dispersivity of 0.02m. The computer simulation quantitatively supports the layered-aquifer model, including under conditions of density-dependent flow, and shows that tidal variations are the predominant driving force for flow beneath the island. The oscillating, vertical flow produced by the tidal variations creates an extensive mixing zone of brackish water. The layered-aquifer model with tidally driven flow is a significant improvement over the Ghyben-Herzberg-Dupuit model as it is conventionally applied to groundwater studies for many Pacific reef islands.     

Poroelastic Effects on Fracture Characterization

Reverse water‐level fluctuations have been widely observed in aquitards or aquifers separated from a pumped confined aquifer (Noordbergum effect) immediately after the initiation of pumping. This same reverse fluctuation has been observed in a fractured crystalline‐rock aquifer at the Coles Hill uranium site in Virginia in which the reverse water‐level response occurs within a pumped fracture and results from an instantaneous strain response to pumping that precedes the pore‐pressure response in observation wells of sufficient distance from the pumped well. This response is referred to as the Mandel‐Cryer effect. The unique aspect of this water level rise during a controlled 24 h pumping test was that the reverse water levels lasted for approximately 100 min and reached a magnitude of nearly 1 cm prior to a typical drawdown response. The duration and magnitude of the response reflects the poromechanical properties of the fractured host rock and hydraulic properties of the pumped fracture. An axisymmetric flow and deformation model were developed using Biot2 in an effort to simulate the observed water‐level response along an assumed 0.5 to 1.0 cm aperture horizontal fracture 176 m from the pumping well and to identify the importance of the poroelastic effect. Results indicate that traditional aquifer‐testing methods that ignore the poromechanical response are not significantly different than results that include the response. However, the poroelastic effect allows for more accurate and efficient parameter calibration.     

Semi‐analytical solution for a slug test in partially penetrating wells including the effect of finite‐thickness skin

This paper presents a new semi‐analytical solution for a slug test in a well partially penetrating a confined aquifer, accounting for the skin effect. This solution is developed based on the solution for a constant‐flux pumping test and a formula given by Peres and co‐workers in 1989. The solution agrees with that of Cooper and co‐workers and the KGS model when the well is fully penetrating. The present solution can be applied to simulate the temporal and spatial head distributions in both the skin and formation zones. It can also be used to demonstrate the influences of skin type or skin thickness on the well water level and to estimate the hydraulic parameters of the skin and formation zones using a least‐squares approach. The results of this study indicate that the determination of hydraulic conductivity using a conventional slug‐test data analysis that neglects the presence of a skin zone will give an incorrect result if the aquifer has a skin zone. Copyright © 2008 John Wiley & Sons, Ltd.     

Simultaneous determination of aquifer parameters and zone structures with fuzzy c-means clustering and meta-heuristic harmony search algorithm

This study proposes an inverse solution algorithm through which both the aquifer parameters and the zone structure of these parameters can be determined based on a given set of observations on piezometric heads. In the zone structure identification problem fuzzy c-means (FCM) clustering method is used. The association of the zone structure with the transmissivity distribution is accomplished through an optimization model. The meta-heuristic harmony search (HS) algorithm, which is conceptualized using the musical process of searching for a perfect state of harmony, is used as an optimization technique. The optimum parameter zone structure is identified based on three criteria which are the residual error, parameter uncertainty, and structure discrimination. A numerical example given in the literature is solved to demonstrate the performance of the proposed algorithm. Also, a sensitivity analysis is performed to test the performance of the HS algorithm for different sets of solution parameters. Results indicate that the proposed solution algorithm is an effective way in the simultaneous identification of aquifer parameters and their corresponding zone structures.