Characteristics of preferential flow and groundwater discharge to Shingobee Lake,Minnesota, USA

Small‐scale heterogeneities and large changes in hydraulic gradient over short distances can create preferential groundwater flow paths that discharge to lakes. A 170 m² grid within an area of springs and seeps along the shore of Shingobee Lake, Minnesota, was intensively instrumented to characterize groundwater‐lake interaction within underlying organic‐rich soil and sandy glacial sediments. Seepage meters in the lake and piezometer nests, installed at depths of 0·5 and 1·0 m below the ground surface and lakebed, were used to estimate groundwater flow. Statistical analysis of hydraulic conductivity estimated from slug tests indicated a range from 21 to 4·8 × 10^?3 m day^?1 and small spatial correlation. Although hydraulic gradients are overall upward and toward the lake, surface water that flows onto an area about 2 m onshore results in downward flow and localized recharge. Most flow occurred within 3 m of the shore through more permeable pathways. Seepage meter and Darcy law estimates of groundwater discharge agreed well within error limits. In the small area examined, discharge decreases irregularly with distance into the lake, indicating that sediment heterogeneity plays an important role in the distribution of groundwater discharge. Temperature gradients showed some relationship to discharge, but neither temperature profiles nor specific electrical conductance could provide a more convenient method to map groundwater–lake interaction. These results suggest that site‐specific data may be needed to evaluate local water budget and to protect the water quality and quantity of discharge‐dominated lakes. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of variable‐density groundwater flow on nitrate flux to coastal waters

A two‐dimensional variable‐density groundwater flow and transport model was developed to provide a conceptual understanding of past and future conditions of nitrate (NO₃) transport and estimate groundwater nitrate flux to the Gulf of Mexico. Simulation results show that contaminant discharge to the coast decreases as the extent of saltwater intrusion increases. Other natural and/or artificial surface waters such as navigation channels may serve as major sinks for contaminant loading and act to alter expected transport pathways discharging contaminants to other areas. Concentrations of NO₃ in the saturated zone were estimated to range between 30 and 160 mg?L^?1 as NO₃. Relatively high hydraulic vertical gradients and mixing likely play a significant role in the transport processes, enhancing dilution and contaminant migration to depth. Residence times of NO₃ in the deeper aquifers vary from 100 (locally) to about 300 years through the investigated aquifer system. NO₃ mass fluxes from the shallow aquifers (0 to 5.7 × 10⁴ mg?m^?2?day^?1) were primarily directed towards the navigation channel, which intersects and captures a portion of the shallow groundwater flow/discharge. Direct NO₃ discharge to the sea (i.e. Gulf of Mexico) from the shallow aquifer was very low (0 to 9.0 × 10¹ mg · m^?2?day^?1) compared with discharge from the deeper aquifer system (0 to 8.2 × 10³ mg?m^?2?day^?1). Both model‐calibrated and radiocarbon tracer‐determined contaminant flux estimates reveal similar discharge trends, validating the use of the model for density‐dependent flow conditions. The modelling approach shows promise to evaluate contaminant and nutrient loading for similar coastal regions worldwide. Copyright © 2015 John Wiley & Sons, Ltd.     

Groundwater–surface water interaction in Lake Nasser,Southern Egypt

A cross‐sectional model, based on the two dimensional groundwater flow equation of Edelman, was applied at seven transects distributed over four geological cross sections to estimate groundwater heads and recharge from/or groundwater discharge to Lake Nasser. The lake with a length of 500 km and an average width of 12 km was created over the period 1964–1970, the time for constructing the Aswan High Dam (AHD). The model, constrained by regional‐scale groundwater flow and groundwater head data in the vicinity of the lake, was successfully calibrated to timeseries of piezometeric heads collected at the cross sections in the period 1965–2004. Inverse modeling yielded high values for the horizontal hydraulic conductivity in the range of 6.0 to 31.1 m day^?1 and storage coefficient between 0.01 and 0.40. The results showed the existence of a strong vertical anisotropy of the aquifer. The calibrated horizontal permeability is systematically higher than the vertical permeability (≈1000:1). The calibrated model was used to explore the recharge from/or groundwater discharge to Lake Nasser at the seven transects for a 40‐year period, i.e. from 1965 to 2004. The analysis for the last 20‐year period, 1985–2004, revealed that recharge from Lake Nasser reduced by 37% compared to the estimates for the first 20‐year period, 1965–1984. In the period 1965–2004, seepage of Lake Nasser to the surrounding was estimated at 1.15 × 10⁹ m³ year^?1. This led to a significant rise of the groundwater table. Variance‐based sensitivity and uncertainty analysis on the Edelman results were conducted applying quasi‐Monte Carlo sequences (Latin Hypercube sampling). The maximum standard deviation of the total uncertainty on the groundwater table was 0.88 m at Toshka (west of the lake). The distance from the lake, followed by the storage coefficient and hydraulic conductivity, were identified as the most sensitive parameters. Copyright © 2012 John Wiley & Sons, Ltd.     

Innovative Contaminant Mass Flux Monitoring in an Aquifer Subject to Tidal Effects

Exposure from groundwater contamination to aquatic receptors residing in receiving surface water is dependent upon the rate of contaminated groundwater discharge. Characterization of groundwater fluxes is challenging, especially in coastal environments where tidal fluctuations result in transient groundwater flows towards these receptors. This can also be further complicated by the high spatial heterogeneity of subsurface deposits enhanced by anthropogenic influences such as the mixing of natural sediments and backfill materials, the presence of subsurface built structures such as sheet pile walls or even occurrence of other sources of contaminant discharge. In this study, the finite volume point dilution method (FVPDM) was successfully used to characterize highly transient groundwater flows and contaminant mass fluxes within a coastal groundwater flow system influenced by marked tides. FVPDM tests were undertaken continuously for more than 48 h at six groundwater monitoring wells, in order to evaluate groundwater flow dynamics during several tide cycles. Contaminant concentrations were measured simultaneously which allowed calculating contaminant mass fluxes. The study highlighted the importance of the aquifer heterogeneity, with groundwater fluxes ranging from 10⁻⁷ to 10⁻³ m/s. Groundwater flux monitoring enabled a significant refinement of the conceptual site model, including the fact that inversion of groundwater fluxes was not observed at high tide. Results indicated that contaminant mass fluxes were particularly higher at a specific monitoring well, by more than three orders of magnitude, than at other wells of the investigated aquifer. This study provided crucial information for optimizing further field investigations and risk mitigation measures.     

From streambed temperature measurements to spatial‐temporal flux quantification: using the LPML method to study groundwater–surface water interaction

Knowledge on groundwater–surface water interaction and especially on exchange fluxes between streams and aquifers is an important prerequisite for the study of transport and fate of contaminants and nutrients in the hyporheic zone. One possibility to quantify groundwater–surface water exchange fluxes is by using heat as an environmlental tracer. Modern field equipment including multilevel temperature sticks and the novel open‐source analysis tool LPML make this technique ever more attractive. The recently developed LPML method solves the one‐dimensional fluid flow and heat transport equation by combining a local polynomial method with a maximum likelihood estimator. In this study, we apply the LPML method on field data to quantify the spatial and temporal variability of vertical fluxes and their uncertainties from temperature–time series measured in a Belgian lowland stream. Over several months, temperature data were collected with multilevel temperature sticks at the streambed top and at six depths for a small stream section. Long‐term estimates show a range from gaining fluxes of ?291 mm day^?1 to loosing fluxes of 12 mm day^?1; average seasonal fluxes ranged from ?138 mm day^?1 in winter to ?16 mm day^?1 in summer. With our analyses, we could determine a high spatial and temporal variability of vertical exchange fluxes for the investigated stream section. Such spatial and temporal variability should be taken into account in biogeochemical cycling of carbon, nutrients and metals and in fate analysis of contaminant plumes. In general, the stream section was gaining during most of the observation period. Two short‐term high stream stage events, seemingly caused by blockage of the stream outlet, led to a change in flow direction from gaining to losing conditions. We also found more discharge occurring at the outer stream bank than at the inner one indicating a local flow‐through system. With the conducted analyses, we were able to advance our understanding of the regional groundwater flow system. Copyright © 2015 John Wiley & Sons, Ltd.     

The effect of transient storage on nitrate uptake lengths in streams: an inter‐site comparison

Surface water–groundwater interaction in the hyporheic zone may enhance biogeochemical cycling in streams, and it has been hypothesized that streams exchanging more water with the hyporheic zone should have more rapid nitrate utilization. We used simultaneous conservative solute and nitrate addition tracer tests to measure transient storage (which includes hyporheic exchange and in‐stream storage) and the rate of nitrate uptake along three reaches within the Red Canyon Creek watershed, Wyoming. We calibrated a one‐dimensional transport model, incorporating transient storage (OTIS‐P), to the conservative solute breakthrough curves and used the results to determine the degree of transient storage in each reach. The nitrate uptake length was quantified from the exponential decrease in nitrate concentration with distance during the tracer tests. Nitrate uptake along the most downstream reach of Red Canyon Creek was rapid (turnover time K^?1_c = 32 min), compared with nitrate uptake reported in other studies (K^?1_c = 12 to 551 min), but other sites within the watershed showed little nitrate retention or loss. The uptake length S_w‐NO^?₃ for the most downstream reach was 500 m and the mass transfer coefficient V_f‐NO^?₃ was 6·3 m min^?1. Results from 15 other nitrate‐addition tracer tests were used to create a regression model relating transient storage and measures of stream flow to nitrate uptake length. The model, which includes specific discharge and transient storage area, explains almost half the variability in nitrate uptake length (adjusted R² = 0·44) and is most effective for comparing sites with very different stream characteristics. Although large differences in specific discharge and storage zone area explain inter‐site differences in nitrate uptake, other unmeasured variables, such as available organic carbon and microbial community composition, are likely important for predicting differences in nitrate uptake between sites with similar specific discharge rates and storage zone areas, such as when making intra‐site comparisons. Copyright © 2007 John Wiley & Sons, Ltd.     

Impact of redox and transport processes in a riparian wetland on stream water quality in the Fichtelgebirge region,southern Germany

Biologically mediated redox processes in the riparian zone, like denitrification, can have substantially beneficial impacts on stream water quality. The extent of these effects, however, depends greatly on the hydrological boundary conditions. The impact of hydrological processes on a wetland's nitrogen sink capacity was investigated in a forested riparian fen which is drained by a first‐order perennial stream. Here, we analysed the frequency distributions and time‐series of pH and nitrogen, silica, organic carbon and oxygen concentrations in throughfall, soil solution, groundwater and stream water, and the groundwater levels and stream discharges from a 3‐year period. During baseflow conditions, the stream was fed by discharging shallow, anoxic groundwater and by deep, oxic groundwater. Whereas the latter delivered considerable amounts of nitrogen (～0·37 mg l^?1) to the stream, the former was almost entirely depleted of nitrogen. During stormflow, near‐surface runoff in the upper 30 cm soil layer bypassed the denitrifying zone and added significant amounts to the nitrogen load of the stream. Nitrate‐nitrogen was close to 100% of deep groundwater and stream‐water nitrogen concentration. Stream‐water baseflow concentrations of nitrate, dissolved carbon and silica were about 1·6 mg l^?1, 4 mg l^?1 and 7·5 mg l^?1 respectively, and >3 mg l^?1, >10 mg l^?1 and <4 mg l^?1 respectively during discharge peaks. In addition to that macroscale bypassing effect, there was evidence for a corresponding microscale effect: Shallow groundwater sampled by soil suction cups indicated complete denitrification and lacked any seasonal signal of solute concentration, which was in contrast to piezometer samples from the same depth. Moreover, mean solute concentration in the piezometer samples resembled more that of suction‐cup samples from shallower depth than that of the same depth. We conclude that the soil solution cups sampled to a large extent the immobile soil‐water fraction. In contrast, the mobile fraction that was sampled by the piezometers exhibited substantially shorter residence time, thus being less exposed to denitrification, but predominating discharge of that layer to the stream. Consequently, assessing the nitrogen budget based on suction‐cup data tended to overestimate the nitrogen consumption in the riparian wetland. These effects are likely to become more important with the increased frequency and intensity of rainstorms that are expected due to climate change. Copyright © 2006 John Wiley & Sons, Ltd.     

A multi‐tracer approach to quantifying groundwater inflows to an upland river; assessing the influence of variable groundwater chemistry

Understanding the behaviour and variability of environmental tracers is important for their use in estimating groundwater discharge to rivers. This study utilizes a multi‐tracer approach to quantify groundwater discharge into a 27 km upland reach of the Gellibrand River in southwest Victoria, Australia. Ten sampling campaigns were conducted between March 2011 and June 2012, and the distribution of ²²²Rn activities, Cl and ³H concentrations imply the river receives substantial groundwater inflows. Mass balances based on ²²²Rn, Cl and ³H yield estimates of groundwater inflows that agree to within ± 12%, with cumulative inflows in individual campaigns ranging from 24 346 to 88 467 m³/day along the studied river section. Groundwater discharge accounts for between 10 and 50% of river flow dependent on the time of year, with a high proportion (>40 %) of groundwater sustaining summer flows. Groundwater inflow is largely governed by regional groundwater flowpaths; between 50 and 90% of total groundwater inflows occur along a narrow 5–10 km section where the river intersects the Eastern View Formation, a major regional aquifer. Groundwater ²²²Rn activities over the 16 month period were spatially heterogeneous across the catchment, ranging between 2000 Bq/m³ and 16 175 Bq/m³. Although groundwater ²²²Rn activities display temporal variation, spatial variation in groundwater ²²²Rn is a key control on ²²²Rn mass balances in river catchments where groundwater and river ²²²Rn activities are within an order of magnitude of each other. Calculated groundwater discharges vary from 8.4 to 15 m³/m/day when groundwater ²²²Rn activities are varied by ± 1 σ. Copyright © 2013 John Wiley & Sons, Ltd.     

Screening for Heat Transport by Groundwater in Closed Geothermal Systems

Heat transfer due to groundwater flow can significantly affect closed geothermal systems. Here, a screening method is developed, based on Peclet numbers for these systems and Darcy's law. Conduction‐only conditions should not be expected where specific discharges exceed 10^?8 m/s. Constraints on hydraulic gradients allow for preliminary screening for advection based on rock or soil types. Identification of materials with very low hydraulic conductivity, such as shale and intact igneous and metamorphic rock, allow for analysis with considering conduction only. Variability in known hydraulic conductivity allows for the possibility of advection in most other rocks and soil types. Further screening relies on refinement of estimates of hydraulic gradients and hydraulic conductivity through site investigations and modeling until the presence or absence of conduction can be confirmed.     

A comparison of fibre‐optic distributed temperature sensing to traditional methods of evaluating groundwater inflow to streams

There are several methods for determining the spatial distribution and magnitude of groundwater inputs to streams. We compared the results of conventional methods [dye dilution gauging, acoustic Doppler velocimeter (ADV) differential gauging, and geochemical end‐member mixing] to distributed temperature sensing (DTS) using a fibre‐optic cable installed along 900 m of Ninemile Creek in Syracuse, New York, USA, during low‐flow conditions (discharge of 1·4 m³ s^?1). With the exception of differential gauging, all methods identified a focused, contaminated groundwater inflow and produced similar groundwater discharge estimates for that point, with a mean of 66·8 l s^?1 between all methods although the precision of these estimates varied. ADV discharge measurement accuracy was reduced by non‐ideal conditions and failed to identify, much less quantify, the modest groundwater input, which was only 5% of total stream flow. These results indicate ambient tracers, such as heat and geochemical mixing, can yield spatially and quantitatively refined estimates of relatively modest groundwater inflow even in large rivers. DTS heat tracing, in particular, provided the finest spatial characterization of groundwater inflow, and may be more universally applicable than geochemical methods, for which a distinct and consistent groundwater end member may be more difficult to identify. Copyright © 2011 John Wiley & Sons, Ltd.     

Two‐dimensional flow response to tidal fluctuation in a heterogeneous aquifer‐aquitard system

In alluvial coastal aquifers, finer sediments are preferentially deposited along the downstream direction, so the hydraulic conductivity is generally heterogeneous and changes with distance from the coastline. To investigate the influence of aquifer heterogeneity on seawater‐groundwater interaction, a new two‐dimensional model characterising groundwater flow in an aquifer‐aquitard system was developed assuming that the hydraulic conductivity of the aquifer linearly increases with the distance from the coastline along the inland direction. A closed‐form analytical solution was derived using the separation‐of‐variables method. Comparing the new solution with the numerical solution by comsol Multiphysics (Sweden) based on the finite‐element method, one can see that the new solution agreed with the numerical solution very well except at the early time. We found that both aquitard leakance and the heterogeneity factor (b) could result in the propagation bias. The propagation bias represents the inconsistency between the theoretical calculation and the observed strong attenuation and small time lag between the head and tide fluctuations. The attenuation decreased with perpendicular distance from the coastline (x‐axis), whereas the time lag increased with distance along the x‐axis. The relationship between the time lag and the distance along the x‐axis seemed to be linear when b was 0.001 m^?1, whereas it obeyed a power function when b was greater than 0.01 m^?1. Copyright © 2014 John Wiley & Sons, Ltd.     

Rainfall infiltration‐derived submarine groundwater discharge from multi‐layered aquifer system terminating at the coastline

This article investigates the quantity of submarine groundwater discharge (SGD) from a coastal multi‐layered aquifer system in response to constant rainfall infiltration. The system comprises an unconfined aquifer, a leaky confined aquifer and an aquitard between them and terminates at the coastline. An approximate analytical solution is derived based on the following assumptions: (i) flow is horizontal in the aquifers and vertical in the aquitard, and (ii) flow in the unconfined aquifer is described by nonlinear Boussinesq equation. The analytical solution is compared with numerical solutions of the strictly two‐dimensional nonlinear model to validate the model assumptions used for the analytical solution. The SGD from the leaky confined aquifer increases with the inland rainfall infiltration recharge and the specific leakage of aquitard. The maximum SGD ranges from 1·87 to 10·37 m³ per day per meter of shoreline when rainfall infiltration ranges from 18·2 to 182 mm/year and the specific leakage of aquitard varies from 10^?9 to 10^?1 l/day. Copyright © 2011 John Wiley & Sons, Ltd.     

Groundwater – the disregarded component in lake water and nutrient budgets. Part 2: effects of groundwater on nutrients

Lacustrine groundwater discharge (LGD) transports nutrients from a catchment to a lake, which may fuel eutrophication, one of the major threats to our fresh waters. Unfortunately, LGD has often been disregarded in lake nutrient studies. Most measurement techniques are based on separate determinations of volume and nutrient concentration of LGD: Loads are calculated by multiplying seepage volumes by concentrations of exfiltrating water. Typically low phosphorus (P) concentrations of pristine groundwater often are increased due to anthropogenic sources such as fertilizer, manure or sewage. Mineralization of naturally present organic matter might also increase groundwater P. Reducing redox conditions favour P transport through the aquifer to the reactive aquifer‐lake interface. In some cases, large decreases of P concentrations may occur at the interface, for example, due to increased oxygen availability, while in other cases, there is nearly no decrease in P. The high reactivity of the interface complicates quantification of groundwater‐borne P loads to the lake, making difficult clear differentiation of internal and external P loads to surface water. Anthropogenic sources of nitrogen (N) in groundwater are similar to those of phosphate. However, the environmental fate of N differs fundamentally from P because N occurs in several different redox states, each with different mobility. While nitrate behaves essentially conservatively in most oxic aquifers, ammonium's mobility is similar to that of phosphate. Nitrate may be transformed to gaseous N₂ in reducing conditions and permanently removed from the system. Biogeochemical turnover of N is common at the reactive aquifer‐lake interface. Nutrient loads from LGD were compiled from the literature. Groundwater‐borne P loads vary from 0.74 to 2900 mg PO₄‐P m^?2 year^?1; for N, these loads vary from 0.001 to 640 g m^?2 year^?1. Even small amounts of seepage can carry large nutrient loads due to often high nutrient concentrations in groundwater. Large spatial heterogeneity, uncertain areal extent of the interface and difficult accessibility make every determination of LGD a challenge. However, determinations of LGD are essential to effective lake management. Copyright © 2014 John Wiley & Sons, Ltd.     

Effect of bed topography on soil aggregates transport by rill flow

After its formation, a rill may remain in the field for months, often receiving lower flow rates than the formative discharge. The objective of this work was to evaluate the rill flow transport capacity of soil aggregates at discharges unable to erode the rill, and to analyse the influence of the rill macro‐roughness on this transport process. A non‐erodible rill was built in which roughness was reproduced in detail. In order to assess only the rill macro‐roughness, a flat channel with a similar micro‐roughness to that in the rill replica was built. Rill and channel experiments were carried out at a slope of 8 and at six discharges (8·3 × 10^?5 to 5·2 × 10^?4 m³ s^?1) in the rill, and eight discharges (1·6 × 10^?5 to 5·2 × 10^?4 m³ s^?1) in the channel. Non‐erodible aggregates of three sizes (1–2, 3–5 and 5–10 mm) were released at the inlet of the rill/channel. The number of aggregates received at the outlet was registered. The number and position of the remaining aggregates along the rill/channel were also determined. The rill flow was a major sediment transport mechanism only during the formation of the rill, as during that period the power of the flow was great enough to overcome the influence of the macro‐roughness of the rill bed. At lower discharges the transport capacity in the previously formed rill was significantly less than that in the flat channel under similar slope and discharge. This was determined to be due to local slowing of flow velocities at the exit of rill pools. Copyright © 2006 John Wiley & Sons, Ltd.     

Estimation of vertical water and heat fluxes in the semi-confined aquifers in tokyo metropolitan area,japan

Groundwater circulation is known to be one of the agents responsible for the redistribution of geothermal energy by acting as a source or sink in the course of its movement through porous media. Heat transport in groundwater systems is considered to be a coupled process and the theory based on this was used to analyse temperature profiles of 30 thermally stable observation wells in a deep, semi-confined aquifer system in the Tokyo Metropolitan area. Vertical water fluxes in the semi-confined aquifers and the associated upward heat fluxes were estimated from a heat flux equation that describes convection and conduction processes of heat transport in one dimension. The vertical downward water fluxes in Shitamachi lowland, Musashino and Tachikawa terraces were 0.69.26.91 × 10^?9, 1.46-70.92 × 10^?9 and 2.61.2204 × 10^?9 m/s, respectively. A vertical upward water flux of 1.80-33.60 × 10^?9 m/s was estimated in Shitamachi lowland. The water flux generally decreased with increasing depth for observation wells which intercepted more than one semi-confining layer. The estimated upward heat fluxes for Shitamachi lowland, Musashino and Tachikawa terraces were 0.32-1.12, 0.49-1.21 and 1.00-11.62 W/m², respectively. The heat flux was highest in Tachikawa terrace where a major fault, the Tachikawa fault, is located. Generally, the estimated heat flux was higher in the semi-confining layers than in the aquifers. Areas with heat sources and sinks as well as groundwater flow patterns in the semi-confined aquifers were revealed by heat flux and temperature distributions in the study area.     

A hydrologic assessment of a saline‐spring fen in the Athabasca oil sands region,Alberta, Canada – a potential analogue for oil sands reclamation

Canada's post‐mined oil sands will have a higher concentration of salts compared with freshwater peatlands that dominate the landscape. While rare, naturally occurring saline wetlands do exist in Alberta's Boreal Plains and may function as analogues for reclamation, however, little is known about their hydrology. This paper investigates the geochemical and hydrologic characteristics of a natural saline‐spring peatland in Alberta's oil sands region. The fen is located within a saline groundwater discharge area connected to the erosional edge of the Grand Rapids Formation. Na⁺ (195–25,680 mgl^?1) and Cl^? (1785–56,249 mg l^?1) were the dominant salts, and the fen transitioned sharply to freshwater along its margins because in part of subsurface mineral ridges that restricted shallow groundwater exchange. Salinity decreased from hypersaline to brackish along the local groundwater flow path but no active spring outlets were observed over the two‐year study. Vertical groundwater discharge was minimal because of the very low permeability of the underlying sediments. Subsurface storage was exceeded during periods of high flow, resulting in flooding and surface runoff that was enhanced by the ephemerally connected pond network. These findings have implications for reclamation, as mechanisms such as subsurface mineral ridges may function as effective saline groundwater‐control structures in the post‐mined environment. Incorporating saline wetlands into regional monitoring networks will help to better quantify natural discharge, which has implications for belowground wastewater storage related to in situ bitumen extraction. Copyright © 2015 John Wiley & Sons, Ltd.     

Formation waters discharge to rivers near oil sands projects

Groundwater discharges in the western Canadian oil sands region impact river water quality. Mapping groundwater discharges to rivers in the oil sands region is important to target water quality monitoring efforts and to ensure injected wastewater and steam remain sequestered rather than eventually resurfacing. Saline springs composed of Pleistocene‐aged glacial meltwater exist in the region, but their spatial distribution has not been mapped comprehensively. Here we show that formation waters discharge into 3 major rivers as they flow through the Athabasca Oil Sands Region adjacent to many active oil sands projects. These discharges increase river chloride concentrations from river headwaters to downstream reaches by factors of ~23 in the Christina River, ~4 in the Clearwater River, and ~5 in the Athabasca River. Our survey provides further evidence for the substantial impact of formation water discharges on river water quality, even though they comprise less than ~2% of total streamflow. Geochemical evidence supporting formation water discharges as the leading control on river salinity include increases in river chloride concentrations, Na/(Na + Ca) ratios, Cl/(Cl + SO₄) ratios and decreases in ⁸⁷Sr/⁸⁶Sr ratios; each mixing trend is consistent with saline groundwater discharges sourced from Cretaceous or Devonian aquifers. These regional subsurface‐to‐surface connections signify that injected wastewater or steam may potentially resurface in the future, emphasizing the critical importance of mapping groundwater flow paths to understand present‐day streamflow quality and to predict the potential for injected fluids to resurface.     

Evaluation of Hydraulic Properties of Aquitards Using Earthquake‐Triggered Groundwater Variation

The hydraulic properties of aquitards are not easily obtained because monitoring wells are usually installed in aquifers for groundwater resources management. Earthquake‐induced crust stress (strain) triggers groundwater level variations over a short period of time in a large area. These groundwater anomalies can be used to investigate aquifer systems. This study uses a poroelastic model to fit the postseismic variations of groundwater level triggered by the Chi‐Chi earthquake to evaluate the hydraulic properties of aquitards in the Jhoushuei River alluvial fan (JRAF), Taiwan. Six of the adopted eight wells with depths of 70 to 130 m showed good agreement with the recovery theory. The mean hydraulic conductivities (K) of the aquifers for the eight wells are 1.62 × 10^?4 to 9.06 × 10^?4 m/s, and the thicknesses are 18.8 to 46.1 m. The thicknesses of the aquitards are 11.3 to 42.0 m. Under the isotropic assumption for K, the estimated values of K for the aquitards are 3.0 × 10^?8 to 2.1 × 10^?6 m/s, corresponding to a silty medium. The results match the values obtained for the geological material of the drilling core and those reported in previous studies. The estimated values were combined with those given in previous studies to determine the distribution of K in the first two aquitards in the JRAF. The distribution patterns of the aquitards reflect the sedimentary environments and fit the geological material. The proposed technique can be used to evaluate the K value of aquitards using inverse methods. The inversion results can be used in hydrogeological analyses, contaminant modeling, and subsidence evaluation.     

The role of bromide and iodide ions in the salinization mapping of the aquifer of Glafkos River basin (northwest Achaia,Greece)

Seawater intrusion causes many problems for groundwater quality, whereas natural remediation is time consuming. However, in cases where groundwater replenishment is feasible, groundwater quality remediation is possible and rapid. The alluvial aquifer in the lowland of the Glafkos River basin, which extends south of Patras city, was for over 30 years the major water source supplying the broader area. Groundwater quality has been degraded due to seawater intrusion, caused by overpumping and generally by inappropriate groundwater management. During the last decade, groundwater quality has been remedied due to diminished groundwater abstractions. The remediation rate was further higher because of rapid discharge of the brackish groundwater, through wells with freely flowing water in the coastal area, where, however, groundwater quality remains low. This paper deals with the hydrogeochemical processes that take place in the area. It is ascertained that ion exchange and mineral dilution processes are dominant. The ion relations between chloride, bromide and iodide, as well as the distribution maps of their concentrations, were used to determine the spatial distribution of the seawater intrusion front. In the lower part of the area in a distance from 1000 and 1500 m from the coast, the rBr^?/rCl^? ratio showed low values (<2·5 × 10^?3) similar to those of seawater. The rI^?/rCl^? ratio also presented low values (<7 × 10^?5), with the lowest one (2·7 × 10^?5) detected along the coastline. In the upper part of the area, a gradual change of those ratios was observed upstream, until they receive values similar to those of the surface waters of Glafkos River. Copyright © 2007 John Wiley & Sons, Ltd.     

Indirect estimation of ungauged peak discharges in a bedrock channel with reference to design discharge selection

Geomorphological evidence and recent trash lines were used as stage indicators in a step-backwater computer model of high discharges through an ungauged bedrock channel. The simulation indicated that the peak discharge from the 26.7 m² catchment was close to 150m³s^?1 during the passage of Hurricane Charlie in August 1986. This estimate can be compared with an estimate of 130–160 m³s^?1 obtained using the Flood Studies Report (FSR) unit hydrograph methodology. Other palaeostage marks indicate that higher stages have occurred at an earlier time associated with a discharge of 200 m³s^?1. However, consideration of both the geometry of a plunge pool and transport criteria for bedrock blocks in the channel indicates that floods since 1986 have not exceeded 150 m³s^?1. Given that the estimated probable maximum flood (PMF) calculated from revised FSR procedure is at least 240 m³s^?1, it is concluded that compelling evidence for floods equal to the PMF is lacking. Taking into consideration the uncertainty of the discharge estimation, the 1986 flood computed using field evidence has a minimum return period of 100 years using the FSR procedure. This may be compared with a return period for the same event in the neighbouring gauged River Greta of > 100 years and a rainfall return period of 190 years. In as much as discharges of similar order to FSR estimates are indicated, it is concluded (a) that regional geomorphological evidence and flood simulation within ungauged catchments may be useful as a verification for hydrological estimates of recent widespread flood magnitude and (b) that palaeohydraulic computation can be useful in determining the magnitude of the local maximum [historic] flood when determining design discharges for hydraulic structures within specific catchments.