Search

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.

This study investigated collaborative groundwater-flow modeling and scenario analysis in the Little Plover River basin, Wisconsin, USA where an unconfined aquifer supplies groundwater for agricultural irrigation, industrial processing, municipal water supply, and stream baseflow. We recruited stakeholders with diverse interests to identify, prioritize, and evaluate scenarios defined as management changes to the landscape. Using a groundwater flow model, we simulated the top 10 stakeholder-ranked scenarios under historically informed dry, average, and wet weather conditions and evaluated the ability of scenarios to meet government-defined stream flow performance measures. Results show that multiple changes to the landscape are necessary to maintain optimum stream flow, particularly during dry years. Yet, when landscape changes from three scenarios-transferring water from the local waste water treatment plant to basin headwaters, moving municipal wells further from the river and downstream, and converting 240 acre (97 ha) of irrigated land to unirrigated land-were simulated in combination, the probability of meeting or exceeding optimum flows rose to 75, 65, and 34% at upper, mid, and lower stream gages, respectively, in dry climate conditions. Discussions with stakeholders reveal that the collaborative model and scenario analysis process resulted in social learning that built upon the existing complex and dynamic institutional landscape. The approach provided a forum for solution-based discussions, and the model served as an important mediation tool for the development and evaluation of community-defined scenarios in a high conflict environment. Today, stakeholders continue to work collaboratively to overcome challenges and implement voluntary solutions in the basin.

A rural subdivision in south central Wisconsin was instrumented with monitoring wells and lysimeters before, during, and after its construction to examine the impacts of the unsewered subdivision on groundwater quality and quantity. Prior to construction, the 78-acre (32 ha) site was farmland. Sixteen homes were constructed beginning in 2003. Initial monitoring from 2002 to 2005 showed that groundwater beneath the site had been impacted by previous agricultural use, with nitrate-N values as high as 30 mg/L and some detections of the herbicide atrazine. Our 12-year study shows that the transition from agricultural to residential land use has changed groundwater quality in both negative and positive ways. Although groundwater elevations showed typical seasonal fluctuations each year, there were no measurable changes in groundwater levels or general flow directions during the 12-year study period. Chloride values increased in many wells, possibly as a result of road salting or water softener discharge. Nitrate concentrations varied spatially and temporally over the study period, with some initial concentrations substantially above the drinking water standard. In some wells, nitrate and atrazine levels have declined substantially since agriculture ceased. However, atrazine was still present at trace concentrations throughout the site in 2014. Wastewater tracers show there are small but detectable impacts from septic effluent on groundwater quality. Particle traces based on a groundwater flow model are consistent with the hypothesis that septic leachate has impacted groundwater quality.

Pathogen contamination from leaky sanitary sewers poses a threat to groundwater quality in urban areas, yet the spatial and temporal dimensions of this contamination are not well understood. In this study, 16 monitoring wells and six municipal wells were repeatedly sampled for human enteric viruses. Viruses were detected infrequently, in 17 of 455 samples, compared to previous sampling at these wells. Thirteen of the 22 wells sampled were virus-positive at least once. While the highest virus concentrations occurred in shallower wells, shallow and deep wells were virus-positive at similar rates. Virus presence in groundwater was temporally coincident, with 16 of 17 virus-positive samples collected in a six-month period. Detections were associated with precipitation and occurred infrequently during a prolonged drought. The study purposely included sites with sewers of differing age and material. The rates of virus detections in groundwater were similar at all study sites during this study. However, a relationship between sewer age and virus detections emerged when compared to data from an earlier study, conducted during high precipitation conditions. Taken together, these data indicate that sewer condition and climate affect urban groundwater contamination by human enteric viruses.

A rural subdivision in south central Wisconsin was instrumented with monitoring wells and lysimeters before, during, and after its construction to examine the impacts of the unsewered subdivision on groundwater quality and quantity. Prior to construction, the 78-acre (32 ha) site was farmland. Sixteen homes were constructed beginning in 2003. Initial monitoring from 2002 to 2005 showed that groundwater beneath the site had been impacted by previous agricultural use, with nitrate-N values as high as 30 mg/L and some detections of the herbicide atrazine. Our 12-year study shows that the transition from agricultural to residential land use has changed groundwater quality in both negative and positive ways. Although groundwater elevations showed typical seasonal fluctuations each year, there were no measurable changes in groundwater levels or general flow directions during the 12-year study period. Chloride values increased in many wells, possibly as a result of road salting or water softener discharge. Nitrate concentrations varied spatially and temporally over the study period, with some initial concentrations substantially above the drinking water standard. In some wells, nitrate and atrazine levels have declined substantially since agriculture ceased. However, atrazine was still present at trace concentrations throughout the site in 2014. Wastewater tracers show there are small but detectable impacts from septic effluent on groundwater quality. Particle traces based on a groundwater flow model are consistent with the hypothesis that septic leachate has impacted groundwater quality.

Viruses are attractive tracers of short (

Until recently, few water utilities or researchers were aware of possible virus presence in deep aquifers and wells. During 2008 and 2009 we collected a time series of virus samples from six deep municipal water-supply wells. The wells range in depth from approximately 220 to 300 m and draw water from a sandstone aquifer. Three of these wells draw water from beneath a regional aquitard, and three draw water from both above and below the aquitard. We also sampled a local lake and untreated sewage as potential virus sources. Viruses were detected up to 61% of the time in each well sampled, and many groundwater samples were positive for virus infectivity. Lake samples contained viruses over 75% of the time. Virus concentrations and serotypes observed varied markedly with time in all samples. Sewage samples were all extremely high in virus concentration. Virus serotypes detected in sewage and groundwater were temporally correlated, suggesting very rapid virus transport, on the order of weeks, from the source(s) to wells. Adenovirus and enterovirus levels in the wells were associated with precipitation events. The most likely source of the viruses in the wells was leakage of untreated sewage from sanitary sewer pipes.

Groundwater models often serve as management tools to evaluate competing water uses including ecosystems, irrigated agriculture, industry, municipal supply, and others. Depletion potential mapping-showing the model-calculated potential impacts that wells have on stream baseflow-can form the basis for multiple potential management approaches in an oversubscribed basin. Specific management approaches can include scenarios proposed by stakeholders, systematic changes in well pumping based on depletion potential, and formal constrained optimization, which can be used to quantify the tradeoff between water use and stream baseflow. Variables such as the maximum amount of reduction allowed in each well and various groupings of wells using, for example, K-means clustering considering spatial proximity and depletion potential are considered. These approaches provide a potential starting point and guidance for resource managers and stakeholders to make decisions about groundwater management in a basin, spreading responsibility in different ways. We illustrate these approaches in the Little Plover River basin in central Wisconsin, United States-home to a rich agricultural tradition, with farmland and urban areas both in close proximity to a groundwater-dependent trout stream. Groundwater withdrawals have reduced baseflow supplying the Little Plover River below a legally established minimum. The techniques in this work were developed in response to engaged stakeholders with various interests and goals for the basin. They sought to develop a collaborative management plan at a watershed scale that restores the flow rate in the river in a manner that incorporates principles of shared governance and results in effective and minimally disruptive changes in groundwater extraction practices.

Using simple numerical groundwater flow models, we tested the impacts of suburban developments on groundwater levels and discharge to streams. We used lot sizes of 1, 3 and 5 acres (4000, 12,000 and 20,000 m2) with one domestic well per lot that pumped water from shallow aquifers. Our modelling showed that pumping had little impact on water levels and groundwater discharge to streams if the developed area is of a moderate size. However, domestic wells had the potential to impact local groundwater levels and baseflows in large developments. In township-wide development scenarios of 1-acre (4000 m2) lots, simulated drawdowns beneath developed areas ranged from 1 to 18 ft (0.3 to 5.5 m), and baseflow reductions ranged from 20 to 40%. Impacts generally were inversely proportional to lot size, recharge rate and hydraulic conductivity of the aquifer materials. Developments using individual domestic wells have the potential to impact local groundwater levels and surface water features. The impacts can range from...

Septic systems that are built in compliance with regulations are generally not expected to be the cause of groundwater borne disease outbreaks, especially in areas with thick vadose zones. However, this case study demonstrates that a disease outbreak can occur in such a setting and outlines the combination of epidemiological, microbiological, and hydrogeological methods used to confirm the source of the outbreak. In early June 2007, 229 patrons and employees of a new restaurant in northeastern Wisconsin were affected by acute gastroenteritis; 6 people were hospitalized. Epidemiological case-control analysis indicated that drinking the restaurant's well water was associated with illness (odds ratio = 3.2, 95% confidence interval = 0.9 to 11.4, P = 0.06). Microbiological analysis (quantitative reverse transcription-polymerase chain reaction) measured 50 genomic copies per liter of norovirus genogroup I in the well water. Nucleotide sequencing determined the genotype as GI.2 and further showed the identical virus was present in patrons' stool specimens and in the septic tank. Tracer tests using dyes injected at two points in the septic system showed that effluent was traveling from the tanks (through a leaking fitting) and infiltration field to the well in 6 and 15 d, respectively. The restaurant septic system and well (85-m deep, in a fractured dolomite aquifer) both conformed to state building codes. The early arrival of dye in the well, which was 188 m from the septic field and located beneath a 35-m thick vadose zone, demonstrates that in highly vulnerable hydrogeological settings, compliance with regulations may not provide adequate protection from fecal pathogens.

Problem: Groundwater contamination is a concern in rural residential subdivisions where numerous septic systems and private wells are sited in close proximity. Although most state codes regulate the construction and location of private wells, these regulations do not usually account for site-specific conditions that may impact drinking-water quality. Purpose: Groundwater models provide a technical basis for delineating groundwater flow to domestic wells. Despite their widespread use in the hydrologic sciences, planners and developers rarely have access to such models. We aimed to assess existing regulations for domestic wells and septic systems and illustrate how groundwater models can be used to evaluate strategies for additional drinking-water protection in unsewered residential subdivisions. Methods: We developed groundwater flow models for two subdivisions in southern Wisconsin, analyzed the results, and disseminated them to local officials, developers, and residents. Results and conclusions: Models o...

There are few studies on the hydrogeology of sedimentary rock aquitards although they are important controls in regional ground water flow systems. We formulate and test a three-dimensional (3D) conceptual model of ground water flow and hydrochemistry in a fractured sedimentary rock aquitard to show that flow dynamics within the aquitard are more complex than previously believed. Similar conceptual models, based on regional observations and recently emerging principles of mechanical stratigraphy in heterogeneous sedimentary rocks, have previously been applied only to aquifers, but we show that they are potentially applicable to aquitards. The major elements of this conceptual model, which is based on detailed information from two sites in the Maquoketa Formation in southeastern Wisconsin, include orders of magnitude contrast between hydraulic diffusivity (K/S(s)) of fractured zones and relatively intact aquitard rock matrix, laterally extensive bedding-plane fracture zones extending over distances of over 10 km, very low vertical hydraulic conductivity of thick shale-rich intervals of the aquitard, and a vertical hydraulic head profile controlled by a lateral boundary at the aquitard subcrop, where numerous surface water bodies dominate the shallow aquifer system. Results from a 3D numerical flow model based on this conceptual model are consistent with field observations, which did not fit the typical conceptual model of strictly vertical flow through an aquitard. The 3D flow through an aquitard has implications for predicting ground water flow and for planning and protecting water supplies.

Confined aquifers are overlain by low-permeability aquitards that are commonly assumed to protect underlying aquifers from microbial contaminants. However, empirical data on microbial contamination beneath aquitards is limited. This study determined the occurrence of human pathogenic viruses in well water from a deep sandstone aquifer confined by a regionally extensive shale aquitard. Three public water-supply wells were each sampled 10 times over 15 months. Samples were analyzed by reverse transcription-polymerase chain reaction (RT-PCR) for several virus groups and by cell culture for infectious enteroviruses. Seven of 30 samples were positive by RT-PCR for enteroviruses; one of these was positive for infectious echovirus 18. The virus-positive samples were collected from two wells cased through the aquitard, indicating the viruses were present in the confined aquifer. Samples from the same wells showed atmospheric tritium, indicating water recharged within the pastfew decades. Hydrogeologic conditions support rapid porous media transport of viruses through the upper sandstone aquifer to the top of the aquitard 61 m below ground surface. Natural fractures in the shale aquitard are one possible virus transport pathway through the aquitard; however, windows, cross-connecting well bores, or imperfect grout seals along well casings also may be involved. Deep confined aquifers can be more vulnerable to contamination by human viruses than commonly believed.

Sandstones often escape extensive hydrogeologic characterization due to their high primary porosity and perceived homogeneity of permeability. This study provides evidence for laterally extensive, high permeability zones in the Tunnel City Group, an undeformed, Cambrian-aged sandstone unit that exists in the subsurface throughout much of central and southern Wisconsin, USA. Several discrete high-permeability zones were identified in boreholes using flow logging and slug tests, and the interconnectedness of the features was tested using a site-specific numerical model for springs in the region. Explicit incorporation of a high-permeability layer leads to improvements in the flux calibration over simulations that lack the features, thus supporting the hydraulic continuity of high-permeability zones in the sandstone aquifer over tens of kilometers. The results suggest that stratigraphically controlled heterogeneities like contrasts in lithology or bedding-plane fractures, which have been shown to strongly influence the flow of groundwater in more heterogeneous sedimentary rocks, may also deserve close examination in sandstone aquifers.

Hydrogeologists commonly use particle tracking to determine groundwater flow paths by tracing the movement of imaginary particles placed in the flow field of a numerical model. Particle tracking codes such as MODPATH (Pollock 2012) calculate the velocity field from the simulated head solution and use it to track the movement of the particles. The method is commonly used to evaluate model boundary conditions, determine the locations of recharge and discharge areas in a modeled aquifer, and approximate the movement of contaminants (Anderson and Woessner 1992). A particularly common use of particle tracking is to delineate the capture zone of a pumping well, which is the aquifer volume that contributes water to it. Capture zone delineation used for wellhead protection studies is often based on a time-of-travel criterion, such as a 1- or 5-year capture zone (U.S. Environmental Protection Agency 1987). Capture zones in steady-state models can be delineated using particles moving in forward or reverse directions. In forward tracking, particles are placed in the model domain and are traced as they move downgradient toward a sink (a well or discharge area). The capture zone is delineated by the pathlines of those particles that arrive at the sink in a specified time. Reverse tracking involves placing particles at a sink and tracking them in the reverse (upgradient) flow direction to their sources. Because all the particles start at a common location, the ensemble of the pathlines defines the capture zone of the sink. Steady-state flow models are often used with particle-tracking codes to delineate capture zones because of simplicity and the desirability of representing long-term average conditions.

This study integrated surface and subsurface stratigraphic data with geophysical logs and hydrogeologic data in order to characterize the hydraulic properties of the Silurian dolomite in northeastern Wisconsin. Silurian stratigraphy consists of predictable alternations of characteristic facies associations. A vertical profile of hydraulic conductivity, obtained from short-interval packer tests in a core hole that penetrates a majority of the Silurian section, indicates that hydraulic conductivity ranges over five orders of magnitude (10–1 to 10–6 cm/s). Matrix conductivity is generally low and varies with texture; the finer-grained restricted-marine and transitional facies being less permeable than the coarser-grained open-marine facies. High-conductivity values are generally associated with bedding-plane fractures, and fracture frequency is greater in the restricted-marine facies.

A wellhead protection study for the city of Sturgeon Bay, Wisconsin, USA, demonstrates the necessity of combining detailed hydrostratigraphic analysis with groundwater modeling to delineate zones of contribution for municipal wells in a fractured dolomite aquifer. A numerical model (MODFLOW) was combined with a particle tracking code (MODPATH) to simulate the regional groundwater system and to delineate capture zones for municipal wells. The hydrostratigraphic model included vertical and horizontal fractures and high-permeability zones. Correlating stratigraphic interpretations with field data such as geophysical logs, packer tests, and fracture mapping resulted in the construction of a numerical model with five high-permeability zones related to bedding planes or facies changes. These zones serve as major conduits for horizontal groundwater flow. Dipping fracture zones were simulated as thin high-permeability layers. The locations of exposed bedrock and surficial karst features were used to identify areas of enhanced recharge. Model results show the vulnerability of the municipal wells to pollution. Capture zones for the wells extend several kilometers north and south from the city. Travel times from recharge areas to all wells were generally less than one year. The high seasonal variability of recharge in the study area made the use of a transient model necessary.

Viruses are attractive tracers of short (3 year) travel times in aquifers because they have unique genetic signatures, are detectable in trace quantities, and are mobile in groundwater. Virus "snaphots" result from infection and disappearance in a population over time; therefore, the virus snapshot shed in the fecal wastes of an infected population at a specific point in time can serve as a marker for tracking virus and groundwater movement. The virus tracing approach and an example application are described to illustrate their ability to characterize travel times in high-groundwater velocity settings, and provide insight unavailable from standard hydrogeologic approaches. Although characterization of preferential flowpaths does not usually characterize the majority of other travel times occurring in the groundwater system (e.g., center of plume mass; tail of the breakthrough curve), virus approaches can trace very short times of transport, and thus can fill an important gap in our current hydrogeology toolbox.

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.

Most wellhead protection studies in fractured-rock aquifers rely on the assumption that the aquifer approximates a porous medium at the scale of the wellhead protection area. Significant errors can result if the assumption is incorrectly applied. Some authors have developed theoretical and subjective criteria for determining when the porous-media approximation is appropriate. Most of these criteria, however, require detailed field work to test the validity of the porous-media approximation. Experiments have been carried out with Roucleau's two-dimensional discrete fracture flow model coupled with a particle-tracking code to determine when the porous-media approximation is appropriate for delineating the capture zone of a well drilled in fractured rock. Specifically, the effects of anisotropy and fracture density on the capture-zone determination have been examined. It has been found that, even in densely fractured aquifers, the zone of contribution determined by the fracture-flow model is significantly larger than the capture zone determined by porous-media-based models.

Piezometer nests were installed at study sites in each of five north-south-trending end moraines of the late Pleistocene Oak Creek Formation in southeastern Wisconsin. The formation is composed primarily of a fine-grained glacial diamicton (till) and laterally continuous and discontinuous, coarse-grained lake and meltwater stream sediment. It overlies the Silurian dolomite aquifer, which is a source of drinking water to rural areas. The average vertical linear velocity and age of ground water in the Oak Creek Formation were estimated by three methods: Darcy's Law, environmental isotopes including 3 H, δ 2 H, δ 18 O, and 14 C (dissolved inorganic carbon), and solute transport modeling of 18 O. The F-1 and Metro sites in the Tinley moraine showed excellent agreement among the three estimates of vertical velocity and showed the lowest velocity values (0.3–0.5 cm year −1 downward), which suggests that diffusion controls vertical mass transport at these sites. Although the extrapolated maximum age of ground water is 35 000 years, ground water below 75 m at these sites is probably not older than 15 000 years, which is the maximum age of the formation. Estimates of velocity showed less agreement at study sites in the Lake Border moraine system to the east and ranged from about 0.2 to 20.7 cm year −1 ; maximum groundwater age could range from 213 to 6000 years. Higher and more variable velocities, perhaps owing to thinner and more heterogeneous sediment in these areas, suggest that diffusion may not dominate vertical mass transport. Heterogeneity and fractures may also promote the development of groundwater flow systems dominated by lateral flow. Because of the uncertainty about the nature of groundwater flow, velocity, and age in the formation east of the Tinley moraine, future waste-disposal activity in southeastern Wisconsin should be confined to the thickest parts of the Tinley moraine near the present F-1 and Metro sites.

The detection of pharmaceuticals and other organic wastewater contaminants (OWCs) in ground water and surface-water bodies has raised concerns about the possible ecological impacts of these compounds on nontarget organisms. On-site wastewater treatment systems represent a potentially significant route of entry for organic contaminants to the environment. In this study, effluent samples were collected and analyzed from conventional septic systems and from systems using advanced treatment technologies. Six of 13 target compounds were detected in effluent from at least one septic system. Caffeine, paraxanthine, and acetaminophen were the most frequently detected compounds, and estrogenic activity was detected in 14 of 15 systems. The OWC concentrations were significantly lower in effluent after sand filtration (p < 0.01) or aerobic treatment (p < 0.05) as compared with effluent that had not undergone advanced treatment. In general, concentrations in conventional systems were comparable to those measured in previous studies of municipal wastewater treatment plant (WWTP) influent, and concentrations in systems after advanced treatment were comparable to previously measured concentrations in WWTP effluent. These data indicate that septic systems using advanced treatment can reduce OWCs in treated effluent to similar concentrations as municipal WWTPs.

Recharge varies spatially and temporally as it depends on a wide variety of factors (e.g. vegetation, precipitation, climate, topography, geology, and soil type), making it one of the most difficult, complex, and uncertain hydrologic parameters to quantify. Despite its inherent variability, groundwater modellers, planners, and policy makers often ignore recharge variability and assume a single average recharge value for an entire watershed. Relatively few attempts have been made to quantify or incorporate spatial and temporal recharge variability into water resource planning or groundwater modelling efforts. In this study, a simple, daily soil–water balance model was developed and used to estimate the spatial and temporal distribution of groundwater recharge of the Trout Lake basin of northern Wisconsin for 1996–2000 as a means to quantify recharge variability. For the 5 years of study, annual recharge varied spatially by as much as 18 cm across the basin; vegetation was the predominant control on this variability. Recharge also varied temporally with a threefold annual difference over the 5-year period. Intra-annually, recharge was limited to a few isolated events each year and exhibited a distinct seasonal pattern. The results suggest that ignoring recharge variability may not only be inappropriate, but also, depending on the application, may invalidate model results and predictions for regional and local water budget calculations, water resource management, nutrient cycling, and contaminant transport studies. Recharge is spatially and temporally variable, and should be modelled as such. Copyright © 2009 John Wiley & Sons, Ltd.

In regions where ground water is generally younger than about 30 years, developing the tritium input history of an area for comparison with the current tritium content of ground water allows quantitative estimates of minimum ground-water age. The tritium input history for central Wisconsin has been constructed using precipitation tritium measured at Madison, Wisconsin and elsewhere. Weighted tritium inputs to ground water reached a peak of over 2,000 TU in 1964, and have declined since that time to about 20-30 TU at present. In the Buena Vista basin in central Wisconsin, most ground-water samples contained elevated levels of tritium, and estimated minimum ground-water ages in the basin ranged from less than one year to over 33 years. Ground water in mapped recharge areas was generally younger than ground water in discharge areas, and estimated ground-water ages were consistent with flow system interpretations based on other data. Estimated minimum ground-water ages increased with depth in areas of downward ground-water movement. However, water recharging through thick moraine sediments was older than water in other recharge areas, reflecting slower infiltration through the sandy till of the moraine.

Confined aquifers are overlain by low-permeability aquitards that are commonly assumed to protect underlying aquifers from microbial contaminants. However, empirical data on microbial contamination beneath aquitards is limited. This study determined the occurrence of human pathogenic viruses in well water from a deep sandstone aquifer confined by a regionally extensive shale aquitard. Three public water-supply wells were each sampled 10 times over 15 months. Samples were analyzed by reverse transcription-polymerase chain reaction (RT-PCR) for several virus groups and by cell culture for infectious enteroviruses. Seven of 30 samples were positive by RT-PCR for enteroviruses; one of these was positive for infectious echovirus 18. The virus-positive samples were collected from two wells cased through the aquitard, indicating the viruses were present in the confined aquifer. Samples from the same wells showed atmospheric tritium, indicating water recharged within the pastfew decades. Hydrogeologic conditions support rapid porous media transport of viruses through the upper sandstone aquifer to the top of the aquitard 61 m below ground surface. Natural fractures in the shale aquitard are one possible virus transport pathway through the aquitard; however, windows, cross-connecting well bores, or imperfect grout seals along well casings also may be involved. Deep confined aquifers can be more vulnerable to contamination by human viruses than commonly believed.

Aquitards protect underlying aquifers from contaminants and limit recharge to those aquifers. Understanding the mechanisms and quantity of ground water flow across aquitards to underlying aquifers is essential for ground water planning and assessment. We present results of laboratory testing for shale hydraulic conductivities, a methodology for determining the vertical hydraulic conductivity (K(v)) of aquitards at regional scales and demonstrate the importance of discrete flow pathways across aquitards. A regional shale aquitard in southeastern Wisconsin, the Maquoketa Formation, was studied to define the role that an aquitard plays in a regional ground water flow system. Calibration of a regional ground water flow model for southeastern Wisconsin using both predevelopment steady-state and transient targets suggested that the regional K(v) of the Maquoketa Formation is 1.8 x 10(-11) m/s. The core-scale measurements of the K(v) of the Maquoketa Formation range from 1.8 x 10(-14) to 4.1 x 10(-12) m/s. Flow through some additional pathways in the shale, potential fractures or open boreholes, can explain the apparent increase of the regional-scale K(v). Based on well logs, erosional windows or high-conductivity zones seem unlikely pathways. Fractures cutting through the entire thickness of the shale spaced 5 km apart with an aperture of 50 microns could provide enough flow across the aquitard to match that provided by an equivalent bulk K(v) of 1.8 x 10(-11) m/s. In a similar fashion, only 50 wells of 0.1 m radius open to aquifers above and below the shale and evenly spaced 10 km apart across southeastern Wisconsin can match the model K(v).

Previous site-specific studies designed to assess the impacts of unsewered subdivisions on ground water quality have relied on upgradient monitoring wells or very limited background data to characterize conditions prior to development. In this study, an extensive monitoring program was designed to document ground water conditions prior to construction of a rural subdivision in south-central Wisconsin. Previous agricultural land use has impacted ground water quality; concentrations of chloride, nitrate-nitrogen, and atrazine ranged from below the level of detection to 296 mg/L, 36 mg/L, and 0.8 microg/L, respectively, and were highly variable from well to well and through time. Seasonal variations in recharge, surface topography, aquifer heterogeneities, surficial loading patterns, and well casing depth explain observed variations in ground water chemistry. This variability would not have been detected if background conditions were determined from only a few monitoring wells or inferred from wells located upgradient of the subdivision site. This project demonstrates the importance of characterizing both ground water quality and chemical variability prior to land-use change to detect any changes once homes are constructed.

[1] Hydraulic heads in a dolomitic shale bedrock aquitard in Wisconsin, USA, are apparently not in equilibrium with drawdown in the underlying aquifer system caused by pumping for municipal supply over the last century. Measurements of head with depth, downhole geophysical logs, and estimates of hydraulic conductivity indicate very low vertical hydraulic diffusivity, and show that high horizontal permeability within the aquitard due to bedding-plane fracture zones can allow lateral groundwater flow. Unlike the hydrogeological conceptual models used in many investigations, flow in aquitards cannot always be inferred to be primarily one-dimensional and vertical. Failure to account for transient conditions and lateral flow in similar settings of intensive groundwater pumping could lead to error in estimates of aquitard leakage and underlying aquifer properties.

There is an acute need for assessment of the groundwater resources of the United States at regional and national scales. This is the principal finding of a National Research Council (NRC) report titled “Investigating Groundwater Systems on Regional and National Scales” published in December 2000. While groundwater knows no political boundaries, groundwater studies, unfortunately do. In the United States, most regional-scale groundwater studies are conducted either by state geological surveys or by the U.S. Geological Survey (USGS). Most groundwater studies conducted by the USGS, with the exception of those of the National Water Quality Assessment (NAWQA) Program, are done under the Federal-State Cooperative Water (Co-op) Program. Co-op studies only occasionally cross county lines, much less state lines. This is because they must be funded at least 50% by a local partner.

The concept of the dividing surface and streamline is introduced to delineate groundwater flow which is intercepted by partially penetrating ditches. The depth of the dividing surface or dividing streamline is used to quantify the degree to which the ditches act as hydraulic barriers to limit the spread of contaminated groundwater in shallow water table aquifers under steady-state conditions. An analytical solution is derived which relates this depth to the water table gradient, the head difference between the ditch stage and the aquifer, the ditch width and the ratio of the vertical to horizontal hydraulic conductivity of the aquifer. The roles of recharge, aquifer heterogeneity and ditch depth in controlling the depth of the dividing surface or streamline, and thereby the effectiveness of interceptor ditches, are examined by using a numerical model.

Studies of the relationships between ground-water systems and surface-water systems (lakes) generally require knowledge of the hydrogeologic properties of lakebed materials. Direct measurement of these properties may be prohibitively expensive or difficult in an offshore environment. A correlation between longitudinal conductance, a geoelectric parameter, and lakebed leakance (vertical hydraulic conductivity divided by thickness) may provide a rapid and inexpensive method for estimating lakebed hydrogeological properties using data collected by offshore seismic and electrical surveys. A test of the method at three study sites in Lake Michigan shows a positive linear relationship between the logarithm of longitudinal conductance and the logarithm of lakebed leakance at the three sites.

We have developed a method to calculate ground- water recharge rates using the mass-balance equation, water- table elevation data, estimates of hydraulic conductivity, and aquifer thickness data, and have applied this method to produce a map of the recharge and discharge patterns for a ground-water basin in central Wisconsin. This recharge mapping method is simplified using a modified computer program, the USGS Modular Groundwater Flow Model (McDonald and Harbaugh, 1984). The modeled recharge pattern compares favorably with a recharge map based on field observations. Because recharge rates are extremely sensitive to hydraulic conductivity, the magnitudes of the calculated rates are less reliable than the patterns of recharge and discharge areas. However, introducing stream discharge data constrains the model to produce net recharge rates averaged over the basin which agree with estimates of the basin yield. Because the method is insensitive to the position of lateral boundaries, it can be used to map recharge over areas within basins that are not physically bounded. Recharge maps made with this method can be used to design ground-water monitoring networks and as frameworks for interpreting geochemical or potentiometric data.

An inverse relationship exists between the longitudinal electrical conductance (S) and the leakance (I.) of sediments between a bedrock aquifer and a surface-water body. The empirical relationship for western Lake Michigan, first reported by Bradbury and Taylor (1984), has been refined after additional measurements. Leakance can now be predicted from longitudinal conductance with 95% confidence to within ±½ order of magnitude. Evidence is also presented that sediment samples collected by coring give erroneously low values for hydraulic conductivity and leakance.

The U.S. Geological Survey's modular code (McDonald and Harbaugh, 1984) was used to reconstruct past ground-water flow paths and make predictions of future flow patterns to understand movement of sulfate-contaminated ground water through a carbonate aquifer near an underground mine. The abandoned mine is in the zinc-lead district near the town of Shullsburg in southwestern Wisconsin. The model predicts observed hydraulic head distributions and ground-water discharge rates with fair accuracy. Instantaneous dewatering rates are underestimated, but long-term rates can be estimated with the model. The persistence of a cone of depression at present and in the future allows uncontaminated water from surrounding areas to dilute the sulfate, as confirmed by chemical sampling. The model predicts radial flow will persist near the Shullsburg mines until approximately 1988. The model represents a significant test of the ability of a ground-water flow model to predict the impacts of underground mining on ground-water systems. The method used required field data collected before, during, and after mining. A model calibrated to just one or two of these data sets is likely to have larger associated uncertainties.

This study compares the accuracy of two types of water table maps both of which were constructed with the object of optimizing future mapping efforts in similar environments. The. first type of map is based solely on office information, with no field verification. The second type of map is based on careful field mapping using numerous measurement points. The office-derived maps were based on topography, surface water features, existing reports, maps and data in the files of the Wisconsin Geological and Natural History Survey; the data were not field-verified. The field-derived maps used a dense network of 236 piezometers at 176 sites in an area of approximately 170 square miles. The field project was much more expensive and labor-intensive than was the construction of office-derived maps for the same area. The two methods produce water table maps which agree to an appreciable extent, the greatest agreement being in areas having ground water-fed streams. Differences in water table elevations indicated by the two methods range from negligible to approximately 5 feet. Thus, depending upon the availability of existing information, relatively accurate water table elevations can be delineated in similar sandy unconfined aquifers without time-consuming and expensive field work that drilling and piezometer installation entails. Preliminary construction of office-derived water table maps enables researchers to use their resources efficiently. In some situations, expensive installation of wells and piezometers for a regional monitoring network may add little accuracy to the regional map. For localized problems, collection of additional field data will always be necessary, but can be guided by the office-derived maps. The authors caution that this technique may only be applicable to sandy, unconfined aquifers in humid climates.

The ground-water pattern in the vicinity of an interceptor ditch used for agricultural drainage is simulated using flow and stream function models to explore the possibility of using interceptor ditches to control the spread of agricultural chemicals in shallow aquifers. The concept of the dividing streamline is used to quantify the effectiveness of the interceptor ditch for limiting the spread of potential contaminants. The simulation results suggest that the interceptor ditch may create an effective hydraulic barrier to the movement of agricultural contaminants in ground water, but its effectiveness is subject to seasonal changes in hydrological conditions, including the regional water-table gradient, the distribution of recharge and discharge, and the ditch stage. A quantitative relationship between effectiveness of the ditch and these hydrological conditions is established.