Influence of upwelling on the shallow water chemistry in a small wetland riparian zone (Basque Country)

Many hydrologic conceptual models in riparian areas assume that the alluvial deposits zone is hydraulically more active than the fractured bedrock below. Therefore, these models undervalue the possible contribution of deeper groundwater from the fractured bedrock system. A hydrochemical study, under various hydrological conditions, has been carried out in a small riparian zone of the Salburua wetland (Basque Country) in order to highlight the conceptual model. This wetland is included in a wide Quaternary aquifer, which has been declared a Nitrate Vulnerable Zone. The results of this study suggest that the fractured bedrock is at least as dynamic as the upper clayey deposits. The presence of more fractured zones, which act as hydraulic “windows”, allow the upwelling of deeper groundwater and, consequently, make the upper alluvial deposits and the fractured bedrock water systems to be cross-connected. Nevertheless, this upwelling is limited to some small areas in the riparian zone. As a result of this local interaction, several chemical reactions have been observed and the hydrochemical characteristics of shallow groundwater undergo seasonal variations. The study shows that a hydrologic conceptual model, which does not consider the hydraulic activity of the fractured bedrock, can be too simplistic. The presence of hydraulic windows could be considered throughout the entire Quaternary aquifer, knowledge of which could help the managers of the Vulnerable Zone and the wetland to take more effective measures for regulation and conservation.     

Using saline tracers to evaluate preferential recharge in fractured rocks,Floyd County,Virginia, USA

Potassium chloride (KCl) and potassium bromide (KBr) tracers were used to explore the role of geologic structure on groundwater recharge and flow at the Fractured Rock Research Site in Floyd County, Virginia, USA. Tracer migration was monitored through soil, saprolite, and fractured crystalline bedrock for a period of 3 months with chemical, physical, and geophysical techniques. The tracers were applied at specific locations on the ground surface to directly test flow pathways in a shallow saprolite and deep fractured-rock aquifer. Tracer monitoring was accomplished with differential electrical resistivity, chemical sampling, and physical monitoring of water levels and spring discharge. KCl, applied at a concentration of 10,000 mg/L, traveled 160 m downgradient through the thrust fault aquifer to a spring outlet in 24 days. KBr, applied at a concentration of 5,000 mg/L, traveled 90 m downgradient through the saprolite aquifer in 19 days. Tracer breakthrough curves indicate diffuse flow through the saprolite aquifer and fracture flow through the crystalline thrust fault aquifer. Monitoring saline tracer migration through soil, saprolite, and fractured rock provided data on groundwater recharge that would not have been available using other traditional hydrologic methods. Travel times and flowpaths observed during this study support preferential groundwater recharge controlled by geologic structure.     

Hydrogeophysical investigation and estimation of groundwater potentials of the Lower Palaeozoic to Precambrian crystalline basement rocks in Keffi area, north-central Nigeria, using resistivity methods

The use of resistivity sounding and two-dimensional (2-D) resistivity imaging was investigated with the aim of delineating and estimating the groundwater potential in Keffi area. Rock types identified are mainly gneisses and granites. Twenty-five resistivity soundings employing the Schlumberger electrode array were conducted across the area. Resistivity sounding data obtained were interpreted using partial curve matching approach and 1-D inversion algorithm, RESIST version 1.0. The 2-D resistivity imaging was also carried out along two traverses using dipole–dipole array, and the data obtained were subjected to finite element method modeling using DIPRO inversion algorithm to produce a two-dimensional subsurface geological model. Interpretation of results showed three to four geoelectrical layers. Layer thickness values were generally less than 2 m for collapsed zone, and ranged from 5 to 30 m for weathered bedrock (saprolite). Two major aquifer units, namely weathered bedrock (saprolite) aquifer and fractured bedrock (saprock) aquifer, have been delineated with the latter usually occurring beneath the former in most areas. Aquifer potentials in the area were estimated using simple schemes that involved the use of three geoelectrical parameters, namely: depth to fresh bedrock, weathered bedrock (saprolite) resistivity and fractured bedrock (saprock) resistivity. The assessment delineated the area into prospective high, medium and low groundwater potential zones.     

Characterization of the shallow groundwater system in an alpine watershed: Handcart Gulch,Colorado, USA

Water-table elevation measurements and aquifer parameter estimates are rare in alpine settings because few wells exist in these environments. Alpine groundwater systems may be a primary source of recharge to regional groundwater flow systems. Handcart Gulch is an alpine watershed in Colorado, USA comprised of highly fractured Proterozoic metamorphic and igneous rocks with wells completed to various depths. Primary study objectives include determining hydrologic properties of shallow bedrock and surficial materials, developing a watershed water budget, and testing the consistency of measured hydrologic properties and water budget by constructing a simple model incorporating groundwater and surface water for water year 2005. Water enters the study area as precipitation and exits as discharge in the trunk stream or potential recharge for the deeper aquifer. Surficial infiltration rates ranged from 0.1–6.2×10^?5 m/s. Discharge was estimated at 1.28×10^?3 km³. Numerical modeling analysis of single-well aquifer tests predicted lower specific storage in crystalline bedrock than in ferricrete and colluvial material (6.7×10^?5–2.0×10^?3 l/m). Hydraulic conductivity in crystalline bedrock was significantly lower than in colluvial and alluvial material (4.3×10^?9–2.0×10^?4 m/s). Water budget results suggest that during normal precipitation and temperatures water is available to recharge the deeper groundwater flow system.     

The relationship between the water table and the surface flow of a losing stream, lower Medano Creek, Great Sand Dunes National Monument, Colorado

 Proposed groundwater withdrawals in the San Luis Valley of Colorado may lower the water table in Great Sand Dunes National Monument. In response, the National Park Service initiated a study that has produced a generalized conceptual model of the hydrologic system in order to assess whether a lowering of the water table might decrease the surface flow of lower Medano Creek. Based upon information obtained during the drilling of several boreholes, there appear to be five important hydrostratigraphic units underlying lower Medano Creek within the upper 30 m of the ground surface: 1. a perched aquifer overlying an aquitard located between about 5 and 6 m below the ground surface; 2. the aquitard itself; 3. an unconfined aquifer located between the upper and lower aquitards; 4. an aquitard located between about 27 and 29 m below the ground surface; and 5. a confined underlying the lower aquitard. Because the areal extent of the aquitards cannot be determined from the borehole data, a detailed conceptual model of the hydrogeologic system underlying lower Medano Creek cannot be developed. However, a generalized conceptual model can be envisioned that consists of a complex system of interlayered aquifers and leaky aquitards, with each aquifer having a unique hydraulic head. Water levels in the perched aquifer rise rapidly to their annual maximum levels in response to the arrival of the flow terminus of Medano Creek during the spring runoff event, and the location of the flow terminus is directly dependent upon the discharge of the creek. Water levels in the deeper, non-perched aquifers do not appear to fluctuate significantly in response to the arrival of the flow terminus, demonstrating that it is unlikely that the proposed groundwater withdrawals will decrease the surface flow of lower Medano Creek. Received: 27 December 1995 · Accepted: 20 February 1996     

Evaluation of the importance of clay confining units on groundwater flow in alluvial basins using solute and isotope tracers: the case of Middle San Pedro Basin in southeastern Arizona (USA)

As groundwater becomes an increasingly important water resource worldwide, it is essential to understand how local geology affects groundwater quality, flowpaths and residence times. This study utilized multiple tracers to improve conceptual and numerical models of groundwater flow in the Middle San Pedro Basin in southeastern Arizona (USA) by determining recharge areas, compartmentalization of water sources, flowpaths and residence times. Ninety-five groundwater and surface-water samples were analyzed for major ion chemistry (water type and Ca/Sr ratios) and stable (¹⁸O, ²H, ¹³C) and radiogenic (³H, ¹⁴C) isotopes, and resulting data were used in conjunction with hydrogeologic information (e.g. hydraulic head and hydrostratigraphy). Results show that recent recharge (<60 years) has occurred within mountain systems along the basin margins and in shallow floodplain aquifers adjacent to the San Pedro River. Groundwater in the lower basin fill aquifer (semi confined) was recharged at high elevation in the fractured bedrock and has been extensively modified by water-rock reactions (increasing F and Sr, decreasing ¹⁴C) over long timescales (up to 35,000 years BP). Distinct solute and isotope geochemistries between the lower and upper basin fill aquifers show the importance of a clay confining unit on groundwater flow in the basin, which minimizes vertical groundwater movement.     

Potential of groundwater contamination by polybrominated diphenyl ethers (PBDEs) in a sensitive bedrock aquifer (Canada)

It is necessary to understand the presence, movement, and persistence of contaminants in aquifers to develop adequate groundwater protection plans. Fractured bedrock aquifers with thin overburden cover are very sensitive to contamination, and little is known about transport processes from the ground surface to depth in this setting. This study was undertaken to investigate the potential of groundwater contamination by polybrominated diphenyl ethers (PBDEs), which are flame retardants, in a natural fractured bedrock aquifer in Canada proven to be sensitive to contamination. PBDEs, which had not been previously measured in groundwater in detail, were detected in the study aquifer at concentrations greater than those observed in surface-water bodies. Potential sources include manure, septic tanks, and the atmosphere. From this scoping study, it is evident that additional surveys of PBDE concentrations in groundwater are warranted, especially in settings with high potential source concentrations coupled with sensitive aquifers.     

Interpreting temporal variations in river response functions: an example from the Arkansas River,Kansas, USA

Groundwater/surface-water interactions can play an important role in management of water quality and quantity, but the temporal and spatial variability of these interactions makes them difficult to incorporate into conceptual models. There are simple methods for identifying the presence of groundwater/surface-water interactions; however, identifying flow mechanisms and pathways can be challenging. More complex methods are available to better identify these mechanisms and pathways but are often too time consuming or costly. In this work, a simple method for interpreting and identifying flow mechanisms and sources using temporal variations of river response functions is presented. This approach is demonstrated using observations from two sites along the Arkansas River in Kansas, USA. A change in flow mechanisms between the rising and falling limbs of river hydrographs was identified, along with a second surface-water source to the aquifer, a finding that was validated with stable isotope analyses.     

Development of a hydrogeological conceptual wetland model in the data-scarce north-eastern region of Kilombero Valley,Tanzania

Understanding groundwater/surface-water interactions in wetlands is crucial because wetlands provide not only a high potential for agricultural production, but also sensitive and valuable ecosystems. This is especially true for the Kilombero floodplain wetland in Tanzania, which represents a data-scarce region in terms of hydrological and hydrogeological data. A comprehensive approach combining hydrogeological with tracer-based assessments was conducted, in order to develop a conceptual hydrogeological wetland model of the area around the city of Ifakara in the north-eastern region of Kilombero catchment. Within the study site, a heterogeneous porous aquifer, with a range of hydraulic conductivities, is underlain by a fractured-rock aquifer. Groundwater chemistry is mainly influenced by silicate weathering and depends on groundwater residence times related to the hydraulic conductivities of the porous aquifer. Groundwater flows from the hillside to the river during most of the year. While floodwater close to the river is mainly derived from overbank flow of the river, floodwater at a greater distance from the river mainly originates from precipitation and groundwater discharge. Evaporation effects in floodwater increase with increasing distance from the river. In general, the contribution of flood and stream water to groundwater recharge is negligible. In terms of an intensification of agricultural activities in the wetland, several conclusions can be drawn from the conceptual model. Results of this study are valuable as a base for further research related to groundwater/surface-water interactions and the conceptual model can be used in the future to set up numerical flow and transport models.     

Geochemical indicators of interbasin groundwater flow within the southern Rio Grande Valley, southwestern USA

Increased groundwater withdrawals for the growing population in the Rio Grande Valley and likely alteration of recharge to local aquifers with climate change necessitates an understanding of the groundwater connection between the Jornada del Muerto Basin and the adjoining and more heavily used aquifer in the Mesilla Basin. Separating the Jornada and Mesilla aquifers is a buried bedrock high from Tertiary intrusions. This bedrock high or divide restricts and/or retards interbasin flow from the Jornada aquifer into the Mesilla aquifer. The potentiometric surface of the southern Jornada aquifer near part of the bedrock high indicates a flow direction away from the divide because of a previously identified damming effect, but a groundwater outlet from the southern Jornada aquifer is necessary to balance inputs from the overall Jornada aquifer. Differences in geochemical constituents (major ions, δD, δ¹⁸O, δ³⁴S, and ⁸⁷Sr/⁸⁶Sr) indicate a deeper connection between the two aquifers through the Tertiary intrusions where Jornada water is geochemically altered because of a geothermal influence. Jornada groundwater likely is migrating through the bedrock high in deeper pathways formed by faults of the Jornada Fault Zone, in addition to Jornada water that overtops the bedrock high as previously identified as the only connection between the two aquifers. Increased groundwater withdrawals and lowering of the potentiometric surface of the Jornada aquifer may alter this contribution ratio with less overtopping of the bedrock high and a continued deeper flowpath contribution that could potentially increase salinity values in the Mesilla Basin near the divide.     

Natural spatial and temporal variations in groundwater chemistry in fractured,sedimentary rocks: scale and implications for solute transport

Natural tracers (major ions, δ¹⁸O, and O₂) were monitored to evaluate groundwater flow and transport to a depth of 20 m below the surface in fractured sedimentary (primarily shale and limestone) rocks. Large temporal variations in these tracers were noted in the soil zone and the saprolite, and are driven primarily by individual storm events. During nonstorm periods, an upward flow brings water with high TDS, constant δ¹⁸O, and low dissolved O₂ to the water table. During storm events, low TDS, variable δ¹⁸O, and high dissolved O₂ water recharges through the unsaturated zone. These oscillating signals are rapidly transmitted along fracture pathways in the saprolite, with changes occurring on spatial scales of several meters and on a time scale of hours. The variations decreased markedly below the boundary between the saprolite and less weathered bedrock. Variations in the bedrock units occurred on time scales of days and spatial scales of at least 20 m. The oscillations of chemical conditions in the shallow groundwater are hypothesized to have significant implications for solute transport. Solutes and colloids that adsorb onto aquifer solids can be released into solution by decreases in ionic strength and pH. The decreases in ionic strength also cause thermodynamic undersaturation of the groundwater with respect to some mineral species and may result in mineral dissolution. Redox conditions are also changing and may result in mineral dissolution/precipitation. The net result of these chemical variations is episodic transport of a wide range of dissolved solutes or suspended particles, a phenomenon rarely considered in contaminant transport studies.     

Hydrogeology of a montane headwater groundwater system downgradient of a coal-mine waste rock dump: Elk Valley,British Columbia,Canada

Steelmaking-coal waste rock placed in mountain catchments in the Elk Valley, British Columbia, Canada, drain constituents of interest (CIs) to surface water downgradient of the waste rock dumps. The role of groundwater in transporting CIs in the headwaters of mountain catchments is not well understood. This study characterizes the physical hydrogeology of a portion of a 10-km² headwater catchment (West Line Creek) downgradient of a 2.7-km² waste rock dump placed over a natural headwater valley-bottom groundwater system. The study site was instrumented with 13 monitoring wells. Drill core samples were collected to determine subsurface lithology and geotechnical properties. The groundwater system was characterized using field testing and water-level monitoring. The valley-bottom sediments were composed of unconsolidated glacial and meltwater successions (<64 m thick) deposited as a series of cut and fill structures overlying shale bedrock. An unconfined basal alluvial aquifer located above fractured bedrock was identified as the primary conduit for groundwater flow toward Line Creek (650 m from the toe of the dump). Discharge through the basal alluvial aquifer was estimated using the geometric mean hydraulic conductivity (±1 standard deviation). These calculations suggest groundwater discharge could account for approximately 15% (ranging from 2 to 60%) of the total water discharged from the watershed. The residence time from the base of the waste rock dump to Line Creek was estimated at <3 years. The groundwater system was defined as a snowmelt (i.e., nival) regime dominated by direct recharge (percolation of precipitation) across the catchment.     

Assessing aquifer susceptibility to and severity of atrazine contamination at a field site in south-central Wisconsin, USA

 A field study from October 1989 through July 1992, conducted on a 4.1-km² area in south-central Wisconsin, USA, examined the distributions of atrazine and its chlorinated metabolites in groundwater and related those distributions to the groundwater flow system. MODFLOW and PATH3D were used to assess bedrock-aquifer susceptibility to contamination. Estimated travel time from water table to bedrock surface ranges from <0.25 to >512 yr. Spatial distribution of the estimates demonstrates that increased travel time to bedrock can result from the presence of shallow surface-water bodies, greater depths to bedrock, and smaller hydraulic conductivities. Estimated travel times to local domestic wells are inversely related to atrazine and desethylated atrazine concentrations observed in water from those wells. The potential impact of long-term atrazine use on aquifer water quality was investigated using MT3D in two best-case scenarios. Uncertainties associated with predicted atrazine concentrations at various depths and times were estimated. For shallow groundwater, widespread violations of Wisconsin's current preventive action limit were predicted, but with large uncertainty stemming from uncertain estimates of input parameter values. The simulations indicate, however, that moderate inputs at the water table are very unlikely to produce violations of Wisconsin's standards deeper in the aquifer. Received, October 1997 Revised, July 1998 Accepted, July 1998     

The role of buried bedrock valleys on the development of karstic aquifers in flat-lying carbonate bedrock: insights from Guelph, Ontario, Canada

A spatial relationship between high capacity municipal production wells (>5,000 m³/day), completed in a deep bedrock aquifer, and a buried bedrock valley was recognized in the city of Guelph, southwestern Ontario, Canada. Most production wells are completed in a discrete zone, ～60 m below ground surface, within flat-lying dolostones of the Silurian Amabel Formation. Thick overburden and limited subsurface data make characterization of the karstic aquifer difficult. This study integrates hydrogeologic data with models of karst formation, deriving a conceptual model of porosity development as it relates to valley incision. Bedrock valley incision likely occurred prior to the early Wisconsinan age (>60–75 ka). Incision created steep hydraulic gradients within the flat-lying bedrock, and provided the driving force required to integrate regional groundwater flow into karst conduits that drained at the base of the valley. Dissolution in production zone dolostones was favoured over dissolution in shallower bedrock due to abundant bedding plane partings and fossiliferous facies with high intercrystalline porosity. Burial of the valley during subsequent ice advances reduced the valley’s hydraulic influence and the efficacy of the flow system to cause dissolution. The high capacity municipal wells near the buried bedrock valley tap into the now dormant karst aquifer system.     

Stream-aquifer interactions in the Maules Creek catchment, Namoi Valley, New South Wales, Australia

The interaction between surface-water streams and groundwater in the Maules Creek catchment of northern New South Wales, Australia has been investigated using a wide range of techniques. Zones of groundwater discharge were mapped by measuring the temperature and fluid electrical-conductivity distribution in bores and surface water. Zones where surface water appears to be recharging the aquifer were investigated by measuring the vertical head gradient between the stream and adjacent bores and by estimates of the decreasing surface flow. Geological heterogeneity appears to be the most significant factor in controlling exchange. Lithological information was assembled using geophysical logging of existing bores, supplemented by the results of electrical resistivity imaging. A preliminary water balance was assembled from the available State records of groundwater abstraction for irrigation, rainfall, evapotranspiration and flow gauging in Maules Creek and the adjacent Namoi River. The analysis has demonstrated the complexity of these coupled systems and gives an indication of the most efficient techniques to be deployed in the field to investigate these complex but important systems.     

Interaction between groundwater, the hyporheic zone and a Chalk stream: a case study from the River Lambourn, UK

Understanding the processes controlling groundwater/surface-water interaction is essential for effective resource management and for protecting sensitive ecosystems. Through intensive monitoring of Chalk groundwater, shallow gravel groundwater and surface water in the River Lambourn, UK, using a combination of hydrochemical and hydrophysical techniques, a complex pattern of interactions has been elucidated. The river is broadly in hydraulic contact with the streambed sediments and adjacent gravels and sands, but these deposits are mainly hydraulically separate from the underlying Chalk at the site. The hydraulic relationship between the river and underlying alluvium is variable, involving components of groundwater flow both parallel and transverse to the river and with both effluent and influent behaviour seen. While the gravel aquifer is significant in controlling groundwater/surface-water interaction, its importance as a route for flow down the catchment is likely to be modest compared with river discharge. The hydrological complexity revealed in a geological setting typical of lowland UK Chalk streams has implications both for investigation methods and for management such as in the setting of environmental objectives in the European Water Framework Directive.     

Potential contaminant transport in the regional Carbonate Aquifer beneath Yucca Mountain, Nevada, USA

Yucca Mountain, Nevada is the site of the proposed US geologic repository for spent nuclear fuel and high-level radioactive waste. The repository is to be a mine, sited approximately 300 m below the crest of the mountain, in a sequence of variably welded and fractured mid-Miocene rhylolite tuffs, in the unsaturated zone, approximately 300 m above the water table. Beneath the proposed repository, at a depth of 2 km, is a thick sequence of Paleozoic carbonate rocks that contain the highly transmissive Lower Carbonate Aquifer. In the area of Yucca Mountain the Carbonate Aquifer integrates groundwater flow from north of the mountain, through the Amargosa Valley, through the Funeral Mountains to Furnace Creek in Death Valley, California where the groundwater discharges in a set of large springs. Data that describe the Carbonate Aquifer suggest a concept for flow through the aquifer, and based upon the conceptual model, a one-layer numerical model was constructed to simulate groundwater flow in the Carbonate Aquifer. Advective transport analyses suggest that the predicted travel time of a particle from Yucca Mountain to Death Valley through the Carbonate Aquifer might be as short as 100 years to as long 2,000 years, depending upon the porosity.     

Arsenic variability and groundwater age in three water supply wells in southeast New Hampshire

Three wells in New Hampshire were sampled bimonthly over three years to evaluate the temporal variability of arsenic concentrations and groundwater age.All samples had measurable concentrations of arsenic throughout the entire sampling period and concentrations in individual wells had a mean variation of more than 7 μg/L.The time series data from this sampling effort showed that arsenic concentrations ranged from a median of 4 μg/L in a glacial aquifer well(SGW-65)to medians of 19μg/L and37 μg/L in wells(SGW-93 and KFW-87)screened in the bedrock aquifer,respectively.These high arsenic concentrations were associated with the consistently high pH(median≥8)and low dissolved oxygen(median0.1 mg/L)in the bedrock aquifer wells,which is typical of fractured crystalline bedrock aquifers in New Hampshire.Groundwater from the glacial aquifer often has high dissolved oxygen,but in this case was consistently low.The pH also is generally acidic in the glacial aquifer but in this case was slightly alkaline(median = 7.5).Also,sorption sites may be more abundant in glacial aquifer deposits than in fractured bedrock which may contribute to lower arsenic concentrations.Mean groundwater ages were less than 50 years old in all three wells and correlated with conservative tracer concentrations,such as chloride;however,mean age was not directly correlated with arsenic concentrations.Arsenic concentrations at KFW-87 did correlate with water levels,in addition,there was a seasonal pattern,which suggests that either the timing of or multiple sampling efforts may be important to define the full range of arsenic concentrations in domestic bedrock wells.Since geochemically reduced conditions and alkaline pHs are common to both bedrock and glacial aquifer wells in this study,groundwater age correlates less strongly with arsenic concentrations than geochemical conditions.There also is evidence of direct hydraulic connection between the glacial and bedrock aquifers,which can influence arsenic concentrations.Correlations between arsenic concentrations and the age of the old fraction of water in SGW-65 and the age of the young fraction of water in SGW-93 suggest that water in the two aquifers may be mixing or at least some of the deeper,older water captured by the glacial aquifer well may be from a similar source as the shallow young groundwater from the bedrock aquifer.The contrast in arsenic concentrations in the two aquifers may be because of increased adsorption capacity of glacio-fluvial sediments,which can limit contaminants more than fractured rock.In addition,this study illustrates that long residence times are not necessary to achieve more geochemically evolved conditions such as high pH and reduced conditions as is typically found with older water in other regions.     

A multidisciplinary-based conceptual model of a fractured sedimentary bedrock aquitard: improved prediction of aquitard integrity

A hydrogeologic conceptual model that improves understanding of variability in aquitard integrity is presented for a fractured sedimentary bedrock unit in the Cambrian-Ordovician aquifer system of midcontinent North America. The model is derived from multiple studies on the siliciclastic St. Lawrence Formation and adjacent strata across a range of scales and geologic conditions. These studies employed multidisciplinary techniques including borehole flowmeter logging, high-resolution depth-discrete multilevel well monitoring, fracture stratigraphy, fluorescent dye tracing, and three-dimensional (3D) distribution of anthropogenic tracers regionally. The paper documents a bulk aquitard that is highly anisotropic because of poor connectivity of vertical fractures across matrix with low permeability, but with ubiquitous bed parallel partings. The partings provide high bulk horizontal hydraulic conductivity, analogous to aquifers in the system, while multiple preferential termination horizons of vertical fractures serve as discrete low vertical hydraulic conductivity intervals inhibiting vertical flow. The aquitard has substantial variability in its ability to protect underlying groundwater from contamination. Across widespread areas where the aquitard is deeply buried by younger bedrock, preferential termination horizons provide for high aquitard integrity (i.e. protection). Protection is diminished close to incised valleys where stress release and weathering has enhanced secondary pore development, including better connection of fractures across these horizons. These conditions, along with higher hydraulic head gradients in the same areas and more complex 3D flow where the aquitard is variably incised, allow for more substantial transport to deeper aquifers. The conceptual model likely applies to other fractured sedimentary bedrock aquitards within and outside of this region.     

Understanding groundwater processes by representing aquifer heterogeneity in the Maules Creek Catchment, Namoi Valley (New South Wales, Australia)

A FEFLOW three-dimensional (3D) groundwater model is developed to enhance the understanding of groundwater processes in the complex alluvial stratigraphy of Maules Creek Catchment (New South Wales, Australia). The aquifer vertical heterogeneity is replicated by indexing 204 lithological logs into units of high or low hydraulic conductivity, and by developing a 3D geological conceptual model with a vertical resolution based on the average lithological unit thickness for the region. The model mesh is populated with the indexed geology using nearest neighbour gridding. The calibrated model is successful in simulating the observed flow dynamics and in quantifying the important water-budget components. This indicates that the lateral groundwater flow from the mountainous region is the main inflow component of the system. Under natural conditions, the Namoi River acts as a sink of water, but groundwater abstraction increasingly removes a large amount of water each year causing dewatering of the system. The pumping condition affects the river–aquifer interaction by reversing the flow, from gaining to losing river conditions during the simulation period. The procedure is relevant for the development of groundwater models of heterogeneous systems in order to improve the understanding of the interplay between aquifer architecture and groundwater processes.