Performance Assessments of a Novel Well Design for Reducing Exposure to Bedrock‐Derived Arsenic

Arsenic in groundwater is a serious problem in New England, particularly for domestic well owners drawing water from bedrock aquifers. The overlying glacial aquifer generally has waters with low arsenic concentrations but is less used because of frequent loss of well water during dry periods and the vulnerability to surface‐sourced bacterial contamination. An alternative, novel design for shallow wells in glacial aquifers is intended to draw water primarily from unconsolidated glacial deposits, while being resistant to drought conditions and surface contamination. Its use could greatly reduce exposure to arsenic through drinking water for domestic use. Hypothetical numerical models were used to investigate the potential hydraulic performance of the new well design in reducing arsenic exposure. The aquifer system was divided into two parts, an upper section representing the glacial sediments and a lower section representing the bedrock. The location of the well, recharge conditions, and hydraulic properties were systematically varied in a series of simulations and the potential for arsenic contamination was quantified by analyzing groundwater flow paths to the well. The greatest risk of arsenic contamination occurred when the hydraulic conductivity of the bedrock aquifer was high, or where there was upward flow from the bedrock aquifer because of the position of the well in the flow system.     

The Effect of Modeled Recharge Distribution on Simulated Groundwater Availability and Capture

Simulating groundwater flow in basin‐fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin‐fill aquifers by direct infiltration and transport through faults and fractures in the high‐elevation areas, by flowing overland through high‐elevation areas to infiltrate at basin‐fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin‐fill aquifers by calibrating a groundwater‐flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady‐state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.     

Simulation of groundwater flow in a crystalline rock aquifer system in Southern Ghana – An evaluation of the effects of increased groundwater abstraction on the aquifers using a transient groundwater flow model

Monitored groundwater level data, well logs, and aquifer data as well as the relevant surface hydrological data were used to conceptualise the hydrogeological system of the Densu Basin in Southern Ghana. The objective was to numerically derive the hydraulic conductivity field for better characterization of the aquifer system and for simulating the effects of increasing groundwater abstraction on the aquifer system in the basin. The hydraulic conductivity field has been generated in this study through model calibration. This study finds that hydraulic conductivity ranges between a low of 2 m/d in the middle sections of the basin and about 40 m/d in the south. Clear differences in the underlying geology have been indicated in the distribution of aquifer hydraulic conductivities. This is in consonance with the general assertion that the hydrogeological properties of the aquifers in the crystalline basement terrains are controlled by the degree of fracturing and/or weathering of the country rock. The transient model suggest aquifer specific storage values to range between 6.0 × 10^?5 m^?1 and 2.1 × 10^?4 m^?1 which are within acceptable range of values normally quoted for similar lithologies in the literature. There is an apparent subtle decrease in groundwater recharge from about 13% of the annual precipitation in 2005 to about 10.3% of the precipitation in 2008. The transient model was used to simulate responses of the system to annual increment of groundwater abstraction by 20% at the 2008 recharge rates for the period 2009 – 2024. The results suggest that the system will not be able to sustain this level of abstraction as it would lead to a basin wide drawdown in the hydraulic head by 4 m by the end of the prediction period. It further suggests a safe annual increment in groundwater abstraction by 5% under business as usual recharge conditions. Identification and protection of groundwater recharge areas in the basin are recommended in order to safeguard the integrity of the resource under the scenario of increased abstraction for commercial activities in the basin. Copyright © 2012 John Wiley & Sons, Ltd.     

A numerical investigation of bedrock groundwater recharge and exfiltration on soil mantled hillslopes

Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first-order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated-soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.     

Using Domestic Well Records to Determine Fractured Bedrock Watersheds and Recharge Rates

This study presents an approach for delineating groundwater basins and estimating rates of recharge to fractured crystalline bedrock. It entailed the use of completion report data (boring logs) from 2500 domestic wells in bedrock from the Coventry Quadrangle, which is located in northeastern Connecticut and characterized by metamorphic gneiss and schist. Completion report data were digitized and imported into ArcGIS^® for data analysis. The data were processed to delineate groundwater drainage basins for the fractured rock based on flow conditions and to estimate groundwater recharge to the bedrock. Results indicate that drainage basins derived from surface topography, in general, may not correspond with bedrock drainage basins due to scale. Estimates of recharge to the bedrock for the study area indicate that only a small fraction of the precipitation or the amount of water that enters the overburden recharges the rock. The approach presented here can be a useful method for water resource‐related assessments that involve fractured rock aquifers.     

Recharge and Groundwater Flow Within an Intracratonic Basin,Midwestern United States

The conservative nature of chloride (Cl^?) in groundwater and the abundance of geochemical data from various sources (both published and unpublished) provided a means of developing, for the first time, a representation of the hydrogeology of the Illinois Basin on a basin‐wide scale. The creation of Cl^? isocons superimposed on plan view maps of selected formations and on cross sections across the Illinois Basin yielded a conceptual model on a basin‐wide scale of recharge into, groundwater flow within and through the Illinois Basin. The maps and cross sections reveal the infiltration and movement of freshwater into the basin and dilution of brines within various geologic strata occurring at basin margins and along geologic structures. Cross‐formational movement of brines is also seen in the northern part of the basin. The maps and cross sections also show barriers to groundwater movement created by aquitards resulting in areas of apparent isolation/stagnation of concentrated brines within the basin. The distribution of Cl^? within the Illinois Basin suggests that the current chemical composition of groundwater and distribution of brines within the basin is dependent on five parameters: (1) presence of bedrock exposures along basin margins; (2) permeability of geologic strata and their distribution relative to one another; (3) presence or absence of major geologic structures; (4) intersection of major waterways with geologic structures, basin margins, and permeable bedrock exposures; and (5) isolation of brines within the basin due to aquitards, inhomogeneous permeability, and, in the case of the deepest part of the basin, brine density effects.     

Distinguishing Mountain Front and Mountain Block Recharge in an Intermontane Basin of the Himalayan Region

Intermontane basin aquifers worldwide, particularly in the Himalayan region, are recharged largely by the adjoining mountains. Recharge in these basins can occur either by water infiltrating from streams near mountain fronts (MFs) as mountain front recharge (MFR) or by sub-surface mountain block infiltration as mountain block recharge (MBR). MFR and MBR recharge are challenging to distinguish and are least quantified, considering the lack of extensive understanding of the hydrological processes in the mountains. This study used oxygen and hydrogen isotopes (δ¹⁸O and δ²H), electrical conductivity (EC) data, hydraulic head, and groundwater level data to differentiate MFR and MBR. Groundwater level data provide information about the groundwater-surface water interactions and groundwater flow directions, whereas isotopes and EC data are used to distinguish and quantify different recharge sources. The present methodology is tested in an intermontane basin of the Himalayan region. The results suggest that karst springs (KS) and deep groundwater (DGW) recharge are dominated by snowmelt (47% ± 10% and 46% ± 9%) as MBR from adjacent mountains, insignificantly affected by evaporation. The hydraulic head data and isotopes indicate Quaternary shallow groundwater (SGW) aquifer system recharge as MFR of local meteoric water with significant evaporation. The results indicate several flow paths in the aquifer system, a local flow for KS, intermediate flow for SGW, and regional flow for DGW. The findings will significantly impact water resource management in the area and provide vital baseline knowledge for sustainable groundwater management in other Himalayan intermontane basins.     

Temporal variation of stable isotopes in a precipitation–groundwater system: implications for determining the mechanism of groundwater recharge in high mountain–hills of the Loess Plateau,China

The groundwater in shallow loess aquifers in high mountain–hills in the western Loess Plateau in China is almost the sole water resource for local residents. However, the question about how the loess groundwater naturally circulates in these high mountain–hills, characterized by low precipitation and high potential evaporation, remains unclear. The objectives of this study are to evaluate the application of hydrogen and oxygen isotopes to (1) examine temporal variations of the isotopic composition of precipitation and shallow groundwater and (2) uncover the mechanism of groundwater recharge in high mountain–hills. Results from 2 years of monitoring data show a difference in the stable isotopes for groundwater and local precipitation between the winter and summer periods. Similar to precipitation, stable isotopes in groundwater are observed to be depleted in winter and enriched in summer, particularly in oxygen isotope. A prominent characteristic is that H and O isotopes of groundwater show a very clear response to strong precipitation in the rainy season in 2013. The results highlight that local precipitation is the likely recharge source for groundwater in shallow loess aquifers. Annual recharge from local precipitation maintains the groundwater resource in the shallower loess aquifer. The mechanisms governing shallow loess groundwater recharge in high mountain–hills were evaluated. In addition to possible vertical slow percolation of soil water through the unsaturated zone, rapid groundwater recharge mechanisms have been identified as temporal preferential infiltration through sinkholes, slip surface or landslide surface and through the interface of loess layer and palaeo‐soils. Most groundwater can be recharged after a heavy rainy season. Copyright © 2015 John Wiley & Sons, Ltd.     

Subglacial Meltwater Recharge in the Dongkemadi River Basin,Yangtze River Source Region

Glaciers on the Tibetan Plateau play an important role in the local hydrological cycle. However, there are only few studies on groundwater in the alpine basins in the Tibetan Plateau which considered the effects of glaciers. Glaciers are extensively distributed in the Dongkemadi River Basin, which is a representative alpine basin in the Yangtze River source region. This study focuses on building a numerical groundwater flow model with glaciations using HydroGeoSphere (HGS) to simulate subglacial meltwater recharge to groundwater in the Dongkemadi River Basin in response to future climate changes. Effects of hydraulic conductivity, precipitation, and temperature on subglacial meltwater recharge to groundwater were discussed. Glacier changes in the future 50 years were predicted under different climate change scenarios. Results show that: (1) the average thickness of the glacier will change significantly; (2) the simulated rate of annual mean subglacial meltwater recharge to groundwater is 4.58 mm, which accounts for 6.33% of total groundwater recharge; and (3) hydraulic conductivity has the largest influence on subglacial meltwater recharge to groundwater, followed by temperature and precipitation. Results of this study are also important to sustainable water resource usage in the Yangtze River source region.     

Conceptualizing leakage and storage contributions from long open interval wells in regional deep basin flow models

Deep basin aquifers are increasingly used in water‐stressed areas, though their potential for sustainable development is inhibited by overlying aquitards and limited recharge rates. Long open interval wells (LOIWs)—wells uncased through multiple hydrostratigraphic units—are present in many confined aquifer systems and can be an important mechanism for deep basin aquifers to receive flow across aquitards. LOIWs are a major control on flow in the deep Cambrian–Ordovician sandstone aquifers of the upper Midwest, USA, providing a source of artificial leakage from shallow bedrock aquifers and equilibrating head within the sandstone aquifers despite differential pumpage. Conceptualizing and quantifying this anthropogenic flow has long been a challenge for groundwater flow modellers, particularly on a regional scale. Synoptic measurements of active production wells and well completion data for northeast Illinois form the basis for a transient, head‐specified MODFLOW model that determines mass balance contributions to the region and estimates LOIW leakage to the aquifers. Using this insight, transient LOIW leakage was simulated using transiently changing K_V zones in a traditional, Q‐specified MODFLOW‐USG model, a novel approach that allows the K_V in a cell containing a LOIW to change transiently by use of the time‐variant materials (TVM) package. With this modification, we achieved a consistent calibration through time, averaging 19.9 m root mean squared error. This model indicates that artificial leakage via LOIWs contributed a minimum of 10–13% of total flow to the sandstone aquifers through the entire history of pumping, up to 50% of flow around 1930. Removal from storage exceeds 40% of flow during peak withdrawals, much of this flow sourced from units other than the primary sandstone aquifers via LOIWs. As such, understanding the timing and magnitude of LOIW leakage is essential for predicting future water availability in deep basin aquifers.     

A low‐dimensional hillslope‐based catchment model for layered groundwater flow

Despite the strong interaction between surface and subsurface waters, groundwater flow representation is often oversimplified in hydrological models. For instance, the interplay between local or shallow aquifers and deeper regional‐scale aquifers is typically neglected. In this work, a novel hillslope‐based catchment model for the simulation of combined shallow and deep groundwater flow is presented. The model consists of the hillslope‐storage Boussinesq (hsB) model representing shallow groundwater flow and an analytic element (AE) model representing deep regional groundwater flow. The component models are iteratively coupled via a leakage term based on Darcy's law, representing delayed recharge to the regional aquifer through a low conductivity layer. Simulations on synthetic single hillslopes and on a two‐hillslope open‐book catchment are presented, and the results are compared against a benchmark three‐dimensional Richards equation model. The impact of hydraulic conductivity, hillslope plan geometry (uniform, convergent, divergent), and hillslope inclination (0.2%, 5%, and 30%) under drainage and recharge conditions are examined. On the single hillslopes, good matches for heads, hydrographs, and exchange fluxes are generally obtained, with the most significant differences in outflows and heads observed for the 30% slope and for hillslopes with convergent geometry. On the open‐book catchment, cumulative outflows are overestimated by 1–4%. Heads in the confined and unconfined aquifers are adequately reproduced throughout the catchment, whereas exchange fluxes are found to be very sensitive to the hillslope drainable porosity. The new model is highly efficient computationally compared to the benchmark model. The coupled hsB/AE model represents an alternative to commonly used groundwater flow representations in hydrological models, of particular appeal when surface–subsurface exchanges, local aquifer–regional aquifer interactions, and low flows play a key role in a watershed's dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

A moving mesh algorithm for 3-D regional groundwater flow with water table and seepage face

A numerical algorithm is described for solving the free-surface groundwater flow equations in 3-D large-scale unconfined aquifers with strongly heterogeneous conductivity and surface recharge. The algorithm uses a moving mesh to track the water-table as it evolves according to kinematic and seepage face boundary conditions. Both steady-state and transient algorithms are implemented in the SECO-Flow 3-D code and demonstrated on stratigraphy based on the Delaware Basin of south-eastern New Mexico.     

Electric-hydraulic conductivity correlation in fractured crystalline bedrock: Central Landfill, Rhode Island, USA

Remote sensing and geoelectrical methods were used to find water-bearing fractures in the Scituate granite under the Central Landfill of Rhode Island. These studies were necessary to evaluate the integrity of the sanitary landfill and for planning safe landfill extensions. The most useful results were obtained with fracture trace analysis using Landsat and SLAR imagery in combination with ground-based resistivity measurements using Schlumberger vertical electrical soundings based on the assumption of horizontally layered strata. Test borings and packer tests confirmed, in the presence of a lineament and low bedrock resistivity, the probable existence of high bedrock fracture density and high average hydraulic conductivity. However, not every lineament was found to be associated with high fracture density and high hydraulic conductivity. Lineaments alone are not a reliable basis for characterising a landfill site as being affected by fractured bedrock. Horizontal fractures were found in borings located away from lineaments. High values of hydraulic conductivity were correlated with low bedrock resistivities. Bedrock resistivities between 60 and 700 Ω m were associated with average hydraulic conductivities between 4 and 60 cm/day. In some cases very low resistivities were confined to the upper part of the bedrock where the hydraulic conductivity was very large. These types of fractures apparently become narrower in aperture with depth. Bedrock zones having resistivities greater than 1000 Ω m showed, without exception, no flow to the test wells. Plots of bedrock resistivity versus the average hydraulic conductivity indicate that the resistivity decreases with increasing hydraulic conductivity. This relationship is inverse to that found in most unconsolidated sediments and is useful for estimating the hydraulic conductivity in groundwater surveys in fractured bedrock. In appropriate settings such as the Central Landfill site in New England, this electric-hydraulic correlation relationship, supplemented by lineament trace analysis, can be used effectively to estimate the hydraulic conductivity in bedrock from only a limited number of resistivity depth soundings and test wells.     

Computer analysis of regional groundwater flow and boundary conditions in the basin of Mexico

Mexico City is situated in the Valley of Mexico on the extensive lacustrine clays that overlay highly productive aquifers of both volcanic and sedimentary origin. The Valley is closed by volcanic mountains. The natural hydraulic boundary conditions associated withe mountain ranges and their relationship to the important aquifers were studied using a two-dimensional, steady-state finite-element model in cross section. Four cross sections were analysed under hydrologic conditions existing prior to the large scale pumping of the aquifers. Factors such as bulk hydraulic conductivities and regional infiltration rates were obtained from field observations and the literature to assess location of the associated groundwater divides, and the water-table in the mountains. The modeled flow patterns are consistent with the historical hydrologic records piezometric characteristics and observed surface features of the groundwater in the Basin of Mexico. From the modeling results, the groundwater recharge in the mountains is 30–50% of the mean average precipitation. Higher and lower rates result in a flow regime that is not compatible with field observations. In general the location of the divides in the mountains is displaced towards the Valley of Mexico, which influences the groundwater budget of the Valley. The water table in places is several hundred metres below ground surface, in accordance with field observations of a very thick unsaturated zone. Before major aquifer exploitation began about 50 years ago, 40–50% of the total discharge into the Valley was by upward flow through the lacustrine deposits. The best results were obtained using a subsurface distribution of hydrostratigraphic units based on recently published geological interpretations.     

The role of recharge and evapotranspiration as hydraulic drivers of ion concentrations in shallow groundwater on Everglades tree islands,Florida (USA)

Recently, evapotranspiration has been hypothesized to promote the secondary formation of calcium carbonate year‐round on tree islands in the Everglades by influencing groundwater ions concentrations. However, the role of recharge and evapotranspiration as drivers of shallow groundwater ion accumulation has not been investigated. The goal of this study is to develop a hydrologic model that predicts the chloride concentrations of shallow tree island groundwater and to determine the influence of overlying biomass and underlying geologic material on these concentrations. Groundwater and surface water levels and chloride concentrations were monitored on eight constructed tree islands at the Loxahatchee Impoundment Landscape Assessment (LILA) from 2007 to 2010. The tree islands at LILA were constructed predominately of peat, or of peat and limestone, and were planted with saplings of native tree species in 2006 and 2007. The model predicted low shallow groundwater chloride concentrations when inputs of regional groundwater and evapotranspiration‐to‐recharge rates were elevated, while low evapotranspiration‐to‐recharge rates resulted in a substantial increase of the chloride concentrations of the shallow groundwater. Modeling results indicated that evapotranspiration typically exceeded recharge on the older tree islands and those with a limestone lithology, which resulted in greater inputs of regional groundwater. A sensitivity analysis indicated the shallow groundwater chloride concentrations were most sensitive to alterations in specific yield during the wet season and hydraulic conductivity in the dry season. In conclusion, the inputs of rainfall, underlying hydrologic properties of tree islands sediments and forest structure may explain the variation in ion concentration seen across Everglades tree islands. Copyright © 2012 John Wiley & Sons, Ltd.     

The chemical and hydrologic structure of Poa´s Volcano, Costa Rica

Comparison of the chemical characteristics of spring and river water draining the flanks of Poa´s Volcano, Costa Rica indicates that acid chloride sulfate springs of the northwestern flank of the volcano are derived by leakage and mixing of acid brines formed in the summit hydrothermal system with dilute flank groundwater. Acid chloride sulfate waters of the Rio Agrio drainage basin on the northwestern flank are the only waters on Poa´s that are affected by leakage of acid brines from the summit hydrothermal system. Acid sulfate waters found on the northwestern flank are produced by the interaction of surface and shallow groundwater with dry and wet acid deposition of SO₂ and H₂SO₄ aerosols, respectively. The acid deposition is caused by a plume of acid gases that is released by a shallow magma body located beneath the active crater of Poa´s.No evidence for a deep reservoir of neutral pH sodium chloride brine is found at Poa´s. The lack of discharge of sodium chloride waters at Poa´s is attributed to two factors: (1) the presence of a relatively volatile-rich magma body degassing at shallow depths (< 1 km) into a high level summit groundwater system; and (2) the hydrologic structure of the volcano in which high rates of recharge combine with rapid lateral flow of shallow groundwater to prevent deep-seated sodium chloride fluids from ascending to the surface. The shallow depth of the volatile-rich magma results in the degassing of large quantities of SO₂ and HCl. These gases are readily hydrolyzed and quickly mix with meteoric water to form a reservoir of acid chloride-sulfate brine in the summit hydrothermal system. High recharge rates and steep hydraulic gradients associated with elevated topographic features of the summit region promote lateral flow of acid brines generated in the summit hydrothermal system. However, the same high recharge rates and steep hydraulic gradients prevent lateral flow of deep-seated fluids, thereby masking the presence of any sodium chloride brines that may exist in deeper parts of the volcanic edifice.Structural, stratigraphic, and topographic features of Poa´s Volcano are critical in restricting flow of acid brines to the northwestern flank of the volcano. A permeable lava-lahar sequence that outcrops in the Rio Agrio drainage basin forms a hydraulic conduit between the crater lake and acid chloride sulfate springs. Spring water residence times are estimated from tritium data and indicate that flow of acid brines from the active crater to the Rio Agrio source springs is relatively rapid (3 to 17 years). Hydraulic conductivity values of the lava-lahar sequence calculated from residence time estimates range from 10⁻⁵ to 10⁻⁷ m/s. These values are consistent with hydraulic conductivity values determined by aquifer tests of fractured and porous lava/pyroclastic sequences at the base of the northwestern flank of the volcano.Fluxes of dissolved rock-forming elements in Rio Agrio indicate that approximately 4300 and 1650 m³ of rock are removed annually from the northwest flank aquifer and the active crater hydrothermal system, respectively. Over the lifetime of the hydrothermal system (100's to 1000's of years), significant increases in aquifer porosity and permeability should occur, in marked contrast to the reduction in permeability that often accompanies hydrothermal alteration in less acidic systems. Average fluxes of fluoride, chloride and sulfur calculated from discharge and compositional data collected in the Rio Agrio drainage basin over the period 1988–1990 are approximately 2, 38 and 30 metric tons/day. These fluxes should be representative of minimum volatile release rates at Poa´s in the last 10 to 20 years.     

Some relationships between lithology,basin form and hydrology: a case study from the Thames basin,UK

The role of lithology in influencing basin form and function is explored empirically by investigating correlations between a range of catchment variables, where the spatial unit of analysis is not surface catchments but lithologically coherent groundwater units. Using the Thames basin, UK, as a case study, nine groundwater units have been identified. Values for 11 hydrological and geomorphological variables, including rainfall, drainage density, Baseflow Index, aquifer porosity, storage coefficient and log‐hydraulic conductivity, aquifer and drainage elevation, river incision, and hypsometric integral have been estimated for each of the groundwater units in the basin, and Pearson correlation coefficients calculated for all pairs of variables. Seven of the correlation coefficients are found to be significant at a confidence level of > 99%. Negative correlations between drainage density and log aquifer hydraulic conductivity, and between drainage density and river incision, and positive correlations between log‐hydraulic conductivity and river incision, log‐hydraulic conductivity and Baseflow Index, and between Baseflow Index and river incision are inferred to have consistent causal explanations. For example, incision of rivers into aquifers leads to relative increases in hydraulic gradients in the vicinity of rivers which, in turn, promotes the development of secondary porosity increasing both aquifer hydraulic conductivity and, hence, Baseflow Index. The implication of this interpretation is that the geomorphological evolution of basins is intimately linked to the evolution of hydraulic conductivity of the underlying aquifers. This is consistent with, and supports the notion of a coupled complexly evolving surface water‐groundwater system. Copyright © 2011 John Wiley & Sons, Ltd.     

Correlating Bedrock Folds to Higher Rates of Arsenic Detection in Groundwater,Southeast Wisconsin,USA

Arsenic in private drinking water wells is a significant problem across much of eastern Wisconsin, USA. The release mechanism and stratigraphic distribution of sulfide and iron (hydr)oxide sources of arsenic in bedrock aquifers are well understood for northeastern Wisconsin. However, recent geologic mapping has identified numerous small bedrock folds to the south, and the impact of these geologic structures on local groundwater flow and well contamination has been little studied. This paper examines the hydrologic and structural effects of the Beaver Dam anticline, southeast Wisconsin, on arsenic in groundwater in the region. Multivariate logistic regression shows wells near the Beaver Dam anticline are statistically more likely to detect arsenic in groundwater compared to wells farther away. Structural and hydrologic changes related to folding are interpreted to be the cause. Core drilled near the fold axis is heavily fractured, and many fractures are filled with sulfides. Elevated hydraulic conductivity estimates are also recorded near the fold axis, which may reflect a higher concentration of vertical fractures. These structural and hydrologic changes may have led to systematic changes in the distribution and concentration of arsenic-bearing mineral hosts, resulting in the observed detection pattern. For areas with similar underlying geology, this approach may improve prediction of arsenic risk down to the local level.     

Contributing Areas to Domestic Wells in Dipping Sedimentary Rocks under Extreme Recharge Events

We use particle tracking to determine contributing areas (CAs) to wells for transient flow models that simulate cyclic domestic pumping and extreme recharge events in a small synthetic watershed underlain by dipping sedimentary rocks. The CAs consist of strike-oriented bands at locations where the water table intersects high-hydraulic conductivity beds, and from which groundwater flows to the pumping well. Factors that affect the size and location of the CAs include topographic flow directions, rock dip direction, cross-bed fracture density, and position of the well relative to streams. For an effective fracture porosity (n_e) of 10⁻⁴, the fastest advective travel times from CAs to wells are only a few hours. These results indicate that wells in this type of geologic setting can be highly vulnerable to contaminants or pathogens flushed into the subsurface during extreme recharge events. Increasing n_e to 10⁻³ results in modestly smaller CAs and delayed well vulnerability due to slower travel times. CAs determined for steady-state models of the same setting, but with long-term average recharge and pumping rates, are smaller than CAs in the models with extreme recharge. Also, the earliest-arriving particles arrive at the wells later in the steady-state models than in the extreme-recharge models. The results highlight the importance of characterizing geologic structure, simulating plausible effective porosities, and simulating pumping and recharge transience when determining CAs in fractured rock aquifers to assess well vulnerability under extreme precipitation events.