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.     

Using hydrogeochemical data to trace groundwater flow paths in a cold alpine catchment

Significant uncertainty remains in understanding the groundwater flow pathways in the northeastern Qinghai–Tibet Plateau. Hydrogeochemical and isotopic data as well as hydrogeological data were combined to explore the groundwater flow path in a representative cold alpine catchment in the headwater region of the Heihe River. The results indicate that the suprapermafrost groundwater chemical components were mainly affected by calcite dissolution and evaporation, whereas the geochemistry of subpermafrost groundwater was controlled by dolomite and gypsum dissolution, calcite precipitation, and albite and halite dissolution. Distinct hydrogeochemical characteristics and controlling processes suggest a poor hydraulic connectivity between the suprapermafrost and subpermafrost groundwater. The hydraulic connectivity between permafrost groundwater and groundwater in the seasonally frozen area was confirmed by their similar hydrogeochemical features. In the seasonally frozen area, a silty clay layer with low permeability separates the aquifer into the deep (depth >20 m) and shallow (depth <20 m) flow paths. The deep groundwater was characterized by the enhanced dedolomitization and enhanced cation exchange processes compared with the shallow groundwater. Groundwater in the seasonally frozen area finally discharges as base flow into the stream. These results provide useful information about the groundwater flow systems in the unique alpine gorge catchments in Qinghai–Tibet Plateau. The above findings suggest that the permafrost distribution and the aquifer structures within the seasonally frozen area have significant impact on groundwater flow paths. Cross‐validation by drilling work and hydrograph data confirms that the hydrogeochemical and isotopic tracers combined with field investigations can be relatively low‐cost tools in interpreting the groundwater flow paths in similar alpine catchments.     

Imaging Hydrological Processes in Headwater Riparian Seeps with Time‐Lapse Electrical Resistivity

Delineating hydrologic and pedogenic factors influencing groundwater flow in riparian zones is central in understanding pathways of water and nutrient transport. In this study, we combined two‐dimensional time‐lapse electrical resistivity imaging (ERI) (depth of investigation approximately 2 m) with hydrometric monitoring to examine hydrological processes in the riparian area of FD‐36, a small (0.4 km²) agricultural headwater basin in the Valley and Ridge region of east‐central Pennsylvania. We selected two contrasting study sites, including a seep with groundwater discharge and an adjacent area lacking such seepage. Both sites were underlain by a fragipan at 0.6 m. We then monitored changes in electrical resistivity, shallow groundwater, and nitrate‐N concentrations as a series of storms transitioned the landscape from dry to wet conditions. Time‐lapse ERI revealed different resistivity patterns between seep and non‐seep areas during the study period. Notably, the seep displayed strong resistivity reductions (～60%) along a vertically aligned region of the soil profile, which coincided with strong upward hydraulic gradients recorded in a grid of nested piezometers (0.2‐ and 0.6‐m depth). These patterns suggested a hydraulic connection between the seep and the nitrate‐rich shallow groundwater system below the fragipan, which enabled groundwater and associated nitrate‐N to discharge through the fragipan to the surface. In contrast, time‐lapse ERI indicated no such connections in the non‐seep area, with infiltrated rainwater presumably perched above the fragipan. Results highlight the value of pairing time‐lapse ERI with hydrometric and water quality monitoring to illuminate possible groundwater and nutrient flow pathways to seeps in headwater riparian areas.     

Groundwater contributions of flow and nitrogen in a headwater agricultural watershed

Nonpoint sources of nitrogen (N) and other nutrients are a major source of water pollution within the Chesapeake Bay watershed and other basins around the world. Human activities associated with agricultural practices can account for a large percentage of N loadings delivered to streams and rivers. This work aims to improve understanding of N transport from groundwater to surface waters, quantifying the principal hydrological processes driving water and N fluxes into and out of a headwater agricultural stream reach. The study site is a 175-m stream reach in a heavily cultivated 40-ha watershed in east-central Pennsylvania. This subwatershed is underlain by fractured shale bedrock, and receives most of its baseflow from groundwater, either by diffuse matrix discharge through the streambed or by localized discharge through riparian seeps. Samples of stream, seep, and shallow groundwater were collected approximately monthly under steady hydrologic conditions in 2017. Calculated matrix flow from hydraulic head and conductivity measurements paired with differential stream gauging was used to solve for the riparian seep flux using a mass balance approach. Riparian seep fluxes ranged from 45 to 217 m³/d, transporting 0.6–4.2 kg N d⁻¹ of nitrate-N from the fractured bedrock aquifer to the stream. Hydrochemical data suggest that the stream is mainly disconnected from the underlying aquifer and that seeps supply essentially all water and N to the system. Seeps are likely sourced with N in nearby agricultural fields and accelerated through the system with shorter residence times than shallow groundwater. Water isotope data reinforced this notion. This study underscores the importance of agriculture as a source of N to ground and surface waters. Identifying source areas that are causing groundwater enrichment of N and seep areas where N discharges to streams is beneficial for developing N pollution mitigation strategies and implementing management practices that aim to reduce nutrient loads to the Chesapeake Bay.     

Modeling a Large-Scale Historic Aquifer Test: Insight into the Hydrogeology of a Regional Fault Zone

Faults can act as flow barriers or conduits to groundwater flow by introducing heterogeneity in permeability. We examine the hydrogeology of the Sandwich Fault Zone, a 137 km long zone of high-angle faults in northern Illinois, using a large-scale historic aquifer test. The fault zone is poorly understood at depth due to the majority of the faults being buried by glacial deposits and its near-vertical orientation which limits geologic sampling across faults. The aquifer test—perhaps one of the largest in terms of overall withdrawal in North American history—was conducted in 1942 at a facility adjacent to the fault zone. More than 34,000 m³/day was pumped for 37 days from nine multiaquifer wells open to the stratified Cambrian-Ordovician sandstone aquifer system. We modeled the aquifer test using a transient MODFLOW-USG model and simulated pumping wells with the CLN package. We tested numerous fault core/damage zone conceptualizations and calibrated to drawdown values recorded at production and observation wells. Our analysis indicates that the fault zone is a low-permeability feature that inhibits lateral movement of groundwater and that there is at least an order of magnitude decrease in horizontal hydraulic conductivity in the fault core compared to the undeformed sandstone. Large head declines have occurred north of the fault zone (over 300 m since predevelopment conditions) and modifying fault zone parameters significantly affects calibration to regional drawdown on a decadal scale. The flow-barrier behavior of the fault zone has important implications for future groundwater availability in this highly stressed region.     

Controls on Ground Water Chemistry in the Central Couloir Sud Rifain,Morocco

Irrigation, urbanization, and drought pose challenges for the sustainable use of ground water in the central Couloir sud rifain, a major agricultural region in north-central Morocco, which includes the cities of Fès and Meknès. The central Couloir is underlain by unconfined and confined carbonate aquifers that have suffered declines in hydraulic head and reductions in spring flow in recent decades. Previous studies have surveyed ground water flow and water quality in wells and springs but have not comprehensively addressed the chemistry of the regional aquifer system. Using graphical techniques and saturation index calculations, we infer that major ion chemistry is controlled (1) in the surficial aquifer by cation exchange, calcite dissolution, mixing with deep ground water, and possibly calcite precipitation and (2) in the confined aquifer and warm springs by calcite dissolution, dolomite dissolution, mixing with water that has dissolved gypsum and halite, and calcite precipitation. Analyses of ²H and ¹⁸O indicate that shallow ground water is affected by evaporation during recharge (either of infiltrating precipitation or return flow), whereas deep ground water is sustained by meteoric recharge with little evaporation. Mechanisms of recharge and hydrochemical evolution are broadly consistent with those delineated for similar regional aquifer systems elsewhere in Morocco and in southern Spain.     

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.     

Toward a new conceptual model for groundwater flow in merokarst systems: Insights from multiple geophysical approaches

Merokarst aquifers — relatively thin (<1–2 m) karstified carbonate units interbedded between mudstone, shale, or sandstone — constitute a significant proportion of carbonate terrain and underlie a large portion of the west- and south-central USA, yet few advances have been made in our understanding of porosity development and flow-path generation in these complex systems in decades. Toward this end, we used a multi-geophysical approach at the well-studied Konza Prairie Biological Station (KPBS), a part of the larger Flint Hills (25,734 km²), underlain by thin limestone units (1–2 m thick) interbedded with mudstone/shale units (2–4 m thick), to elucidate hydrologic connectivity and potential controls on known groundwater flow directions. We combined electrical resistivity tomography (ERT), surface and borehole nuclear magnetic resonance (NMR), and ground penetrating radar (GPR) measurements across a low order catchment where over 25 boreholes and groundwater wells sampling perched aquifers could be used to constrain interpretation of lithology, potential flow paths, and permeability. Data revealed that groundwater export may be an unappreciated component of lateral-flow-dominated models used to represent merokarst in that: (a) potentiometric surfaces from two limestone units showed groundwater flows toward a hydrologic depression, opposite the direction of stream flow, in the upstream portion of the catchment, (b) long term measures of groundwater levels revealed a greater variance and overall water storage in this same upstream area compared to wells near the outlet, and (c) ERT and NMR results indicate pronounced lateral heterogeneity within a given unit, suggestive of a greater degree of vertical hydrological connectivity than usually considered for horizontally-layered merokarst. Our data suggest vertical connectivity can shunt water to depth in these “sandwiched” merokarst aquifers, yielding atypical groundwater flow directions and unrealized deep export of weathering solutes and carbon.     

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.     

Interactions between groundwater and surface water in a virginia coastal plain watershed. 1. Hydrological flowpaths

A field study of surface water and groundwater interactions during baseflow and stormflow conditions was performed at the Reedy Creek watershed in the Virginia Coastal Plain. Three estimates of the average saturated hydraulic conductivity (K_s) of the unconfined aquifer were in reasonable agreement (ranging from 0.0033 to 0.010 cm/s), indicating that baseflow in the creek is entirely from the drainage of shallow groundwater from the relatively thin (1–6 m thick) unconfined aquifer. This relatively permeable surficial aquifer was found to be underlain by dark, olive grey, clay-silt and diatomaceous Miocene deposits of low permeability known as the Calvert Formation, which is believed to function as a confining bed in the area. A chemical hydrograph separation technique was used to resolve the contributions of [old] (pre-event) and [new] (event) water to stormflow. Results from a major rainstorm indicated that old water dominated the stormflow response of the watershed, although the new water contribution approached 40% at the hydrograph peak. Stormflow at Reedy Creek appears to result from saturation overland flow from variable source areas which include the stream channels and a significant part of the riparian wetland area. This response appears to be attributable to the transient dynamics of the shallow groundwater flow system and to the formation of localized groundwater mounds which raise the water-table to the wetland surface.     

Coastal aquifer response to extreme storm events in Emilia‐Romagna,Italy

With global warming and sea level rise, many coastal systems will experience increased levels of inundation and storm flooding, especially along sandy lowland coastal areas, such as the Northern Adriatic coast (Italy). Understanding how extreme events may directly affect groundwater hydrology in shallow unconfined coastal aquifers is important to assess coastal vulnerability and quantify freshwater resources. This study investigates shallow coastal aquifer response to storm events. The transitory and permanent effects of storm waves are evaluated through the real time monitoring of groundwater and soil parameters, in order to characterize both the saturated and unsaturated portions of the coastal aquifer of Ravenna and Ferrara (southern Po Delta, Italy). Results highlight a general increase in hydraulic head and soil moisture, along with a decrease in groundwater salinity and pore water salinity due to rainfall infiltration during the 2 days storm event. The only exceptions are represented by the observation wells in proximity to the coastline (within 100 m), which recorded a temporary increase in soil and water salinity caused by the exceptional high waves, which persist on top of the dune crest during the storm event. This generates a saline plume that infiltrates through the vadose zone down to the saturated portion of the aquifer causing a temporary disappearance of the freshwater lens generally present, although limited in size, below the coastal dunes. Despite the high hydraulic conductivity, the aquifer system does not quickly recover the pre‐storm equilibrium and the storm effects are evident in groundwater and soil parameters after 10 days past the storm overwash recess.     

Groundwater flow and implications for groundwater contamination north of Prewitt,New Mexico,U.S.A.

In the southern San Juan Basin, New Mexico, strata of Permian and younger age dip gently toward the center of the basin. Most previous investigators believed that recharge to these strata occurred by precipitation on the outcrops and groundwater flowed downdip to the north and northeast. Recent water-level measurements in an undeveloped part of the basin near Prewitt, New Mexico, show that groundwater at shallow depths in alluvium and bedrock flows southward, opposite to the dip direction, and toward a major ephemeral drainage in a strike valley. North of this area, groundwater in deep bedrock aquifers does appear to flow northward. This information suggests that there are two groundwater circulation patterns; a shallow one controlled by topography and a deeper one controlled by geologic structure.Significant amounts of recharge to sandstone aquifers by infiltration through outcrops is unlikely due to the near-vertical exposures on cliffs, the gentle dip of the strata, and small annual precipitation. Numerical model results suggest that recharge to bedrock aquifers may be from downward leakage via aquitards over large areas and leakage from narrow alluvial aquifers in the subcrop area. The recharge mechanism is controlled by the hydraulic conductivity of the strata.As the flow path is controlled by hydraulic conductivity contrasts, geologic structure, and topography, contamination movement from surface impoundments is likely to be difficult to predict without a thorough hydrogeological site investigation.     

Assessing hydrogeologic controls on dynamic groundwater storage using long‐term instrumental records of water table levels

This study analyzes a long‐term regional compilation of water table response to climate variability based on 124 long‐term groundwater wells distributed across New England, USA, screened in a variety of geologic materials. The New England region of the USA is located in a humid‐temperature climate underlain by low‐storage‐fractured metamorphic and crystalline bedrock dissected by north–south trending valleys filled with glacial and post‐glacial valley fill sediments. Uplands are covered by thin glacial till that comprises more than 60% of the total area. Annual and multi‐annual responses of the water table to climate variability are assessed to understand how local hydraulic properties and hydrogeologic setting (located in recharge/discharge region) of the aquifer influence the hydrologic sensitivity of the aquifer system to climate variability. This study documents that upland aquifer systems dominated by thin deposits of surface till comprise ~70% of the active and dynamic storage of the region. Total aquifer storage changes of +5 to ?7 km³ occur over the region during the study interval. The storage response is dominated by thin and low permeability surficial till aquifer that fills and drains on a multi‐annual basis and serves as the main mechanism to deliver water to valley fill aquifers and underlying bedrock aquifers. Whereas the till aquifer system is traditionally neglected as an important storage reservoir, this study highlights the importance of a process‐based understanding of how different landscape hydrogeologic units contribute to the overall hydrologic response of a region.     

Patterns and dynamics of river–aquifer exchange with variably-saturated flow using a fully-coupled model

The parallel physically-based surface–subsurface model PARFLOW was used to investigate the spatial patterns and temporal dynamics of river–aquifer exchange in a heterogeneous alluvial river–aquifer system with deep water table. Aquifer heterogeneity at two scales was incorporated into the model. The architecture of the alluvial hydrofacies was represented based on conditioned geostatistical indicator simulations. Subscale variability of hydraulic conductivities (K) within hydrofacies bodies was created with a parallel Gaussian simulation. The effects of subscale heterogeneity were investigated in a Monte Carlo framework. Dynamics and patterns of river–aquifer exchange were simulated for a 30-day flow event. Simulation results show the rapid formation of saturated connections between the river channel and the deep water table at preferential flow zones that are characterized by high conductivity hydrofacies. Where the river intersects low conductivity hydrofacies shallow perched saturated zones immediately below the river form, but seepage to the deep water table remains unsaturated and seepage rates are low. Preferential flow zones, although only taking up around 50% of the river channel, account for more than 98% of total seepage. Groundwater recharge is most efficiently realized through these zones. Subscale variability of K_sat slightly increased seepage volumes, but did not change the general seepage patterns (preferential flow zones versus perched zones). Overall it is concluded that typical alluvial heterogeneity (hydrofacies architecture) is an important control of river–aquifer exchange in rivers overlying deep water tables. Simulated patterns and dynamics are in line with field observations and results from previous modeling studies using simpler models. Alluvial heterogeneity results in distinct patterns and dynamics of river–aquifer exchange with implications for groundwater recharge and the management of riparian zones (e.g. river channel-floodplain connectivity via saturated zones).     

Impact of complex aquifer geometry on groundwater storage in high‐elevation meadows of the Sierra Nevada Mountains,CA

Recent research has indicated that Sierra Nevada meadows are hydrologically more complex than previously considered. Improved understanding of the effects of aquifer parameters and climate change on water resources in and downstream of meadows is critically needed to effectively manage mountain meadows for ecosystem services and watershed contributions. This research investigates the roles of bedrock geometry, saturated hydraulic conductivity, and meadow gradient in affecting groundwater storage dynamics and surface‐water outflows in site‐scale high‐elevation meadows. Under current and projected lower snowpack conditions, we modeled groundwater flow in representative high‐elevation meadows considering 2 conceptual aquifer thickness models: uniform and variable thickness. Spatially, variable aquifer thicknesses interpreted from bedrock depths (0–28 m) were identified from a high‐resolution ground‐penetrating radar survey conducted at Tuolumne Meadows, CA. Our interpreted bedrock surface indicated several buried U‐shaped valleys including a buried ridge that separates 2 U‐shaped valleys. Groundwater flow simulations show that an increase in meadow gradient and hydraulic conductivity led to a decrease in seasonal storage and an increase in surface‐water outflow. However, models with varying bedrock geometries change the magnitude and timing of these processes. Uniform thickness models overestimated storage at the model edges and resulted in higher projected volumes of water being released to streams earlier than previously observed.     

A Single Packer Method for Characterizing Water Contributing Fractures in Crystalline Bedrock Wells

Water levels and water quality of open borehole wells in fractured bedrock are flow-weighted averages that are a function of the hydraulic heads and transmissivities of water contributing fractures, properties that are rarely known. Without such knowledge using water levels and water quality data from fractured bedrock wells to assess groundwater flow and contaminant conditions can be highly misleading. This study demonstrates a cost-effective single packer method to determine the hydraulic heads and transmissivities of water contributing fracture zones in crystalline bedrock wells. The method entails inflating a pipe plug to isolate sections of an open borehole at different depths and monitoring changes in the water level with time. At each depth, the change in water level with time was used to determine the sum of fracture transmissivities above the packer and then to solve for individual fracture transmissivity. Steady-state wellbore heads along with the transmissivities were used to determine individual fracture heads using the weighted average head equation. The method was tested in five wells in crystalline bedrock located at the University of Connecticut in Storrs. The single packer head and transmissivity results were found to agree closely with those determined using conventional logging methods and the dissolved oxygen alteration method. The method appears to be a simple and cost-effective alternative in obtaining important information on flow conditions in fractured crystalline bedrock wells.     

Groundwater occurrence and flow pattern in the Enugu coal-mine area,Anambra State,Nigeria

Abstract

The distribution of major geological units, static water level data, water chemistry data, and observations of surface features influenced by groundwater seepage were used to ascertain the nature of groundwater occurrence and flow pattern in the Enugu coal field, Nigeria. Considerations of the geological units, the static groundwater levels and groundwater seepages in the mines indicate that the coal sequence is a multiaquifer system in which sandstone and coal aquifers alternate with shale aquitards. Based on the hydraulic head data, the groundwater flow is predominantly downwards. Groundwater velocity calculation across the multiaquifer system using the Darcy equation gave a flow velocity of about 1 m day^?1. For groundwater systems, such a calculated velocity is considered high. The high velocity is most probably due to the high fracture porosity as well as the presence of other stratigraphic and structural features such as alluvial fills that provide high hydraulic conductivity pathways across the aquifer system. The pattern of groundwater inflow into the mines is also influenced by these stratigraphie and structural features.     

The origin,typology, spatial distribution and detrimental effects of the sinkholes developed in the alluvial evaporite karst of the Ebro River valley downstream of Zaragoza city (NE Spain)

Three types of sinkhole have been mapped in a 50 km² stretch of the Ebro River valley downstream of Zaragoza: large collapse sinkholes, large shallow subsidence depressions and small cover-collapse sinkholes. The sinkholes relate to the karstification of evaporitic bedrock that wedges out abruptly downstream, giving way to a shale substratum. Twenty-three collapse sinkholes, up to 50 m in diameter by 6 m deep, and commonly hosting saline ponds, have been identified in the floodplain. They have been attributed to the upward stoping of dissolutional cavities formed within the evaporitic bedrock by rising groundwater flows. Twenty-four large shallow subsidence depressions were mapped in the floodplain. These may reach 850 m in length and were formed by structurally controlled interstratal karstification of soluble beds (halite or glauberite? and gypsum) by rising groundwater flow and the progressive settlement of the overlying bedrock and overburden sediments. A total of 447 small cover-collapse, or dropout, sinkholes have been recognized in a perched alluvial level along the southern margin of the valley. These sinkholes result from the upward propagation of voids through the alluvial mantle caused by the downward migration of detrital sediments into dissolutional voids. The majority of these sinkholes, commonly 1·5–2 m in diameter, are induced by human activities. Over the karstic bedrock, there is a significant increase in sinkhole density downstream. This is interpreted as being a result of the evaporitic bedrock wedging out and the convergence of the groundwater flow lines in the karstic aquifer. The collapse sinkholes in this area, locally with a probability of occurrence higher than 140 sinkholes/km²/year, cause substantial damage to the linear infrastructures, buildings and agriculture, and they might eventually cause the loss of human lives. Copyright © 2006 John Wiley & Sons, Ltd.     

Refined conceptualization of TOPMODEL for shallow subsurface flows

The TOPMODEL framework was used to derive expressions that account for saturated and unsaturated flow through shallow soil on a hillslope. The resulting equations were the basis for a shallow‐soil TOPMODEL (STOPMODEL). The common TOPMODEL theory implicitly assumes a water table below the entire watershed and this does not conceptually apply to systems hydrologically controlled by shallow interflow of perched groundwater. STOPMODEL provides an approach for extending TOPMODEL's conceptualization to apply to shallow, interflow‐driven watersheds by using soil moisture deficit instead of water table depth as the state variable. Deriving STOPMODEL by using a hydraulic conductivity function that changes exponentially with soil moisture content results in equations that look very similar to those commonly associated with TOPMODEL. This alternative way of conceptualizing TOPMODEL makes the modelling approach available to researchers, planners, and engineers who work in areas where TOPMODEL was previously believed to be unsuited, such as the New York City Watershed in the Catskills region of New York State. Copyright © 2002 John Wiley & Sons, Ltd.     

Stage‐dependent transient storage of phosphorus in alluvial floodplains

Models for contaminant transport in streams commonly idealize transient storage as a well mixed but immobile system. These transient storage models capture rapid (near‐stream) hyporheic storage and transport, but do not account for large‐scale, stage‐dependent interaction with the alluvial aquifer. The objective of this research was to document transient storage of phosphorus (P) in coarse gravel alluvium potentially influenced by large‐scale, stage‐dependent preferential flow pathways (PFPs). Long‐term monitoring was performed at floodplain sites adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping which was correlated to hydraulic conductivity data, observation wells were installed both in higher hydraulic conductivity and lower hydraulic conductivity subsoils. Water levels in the wells were monitored over time, and water samples were obtained from the observation wells and the stream to document P concentrations at multiple times during high flow events. Contour plots indicating direction of flow were developed using water table elevation data. Contour plots of total P concentrations showed the alluvial aquifer acting as a transient storage zone, with P‐laden stream water heterogeneously entering the aquifer during the passage of a storm pulse, and subsequently re‐entering the stream during baseflow conditions. Some groundwater in the alluvial floodplains had total P concentrations that mirrored the streams' total P concentrations. A detailed analysis of P forms indicated that particulate P (i.e. P attached to particulates greater than 0·45 µm) was a significant portion of the P transport. This research suggests the need for more controlled studies on stage‐dependent transient storage in alluvial systems. Copyright © 2011 John Wiley & Sons, Ltd.