Interaction of groundwater with Lake Urmia in Iran

Being a large hyper‐saline water body, Lake Urmia in north‐western Iran deals with a gradual decline in its water level. Most of the studies on Lake Urmia have neglected the groundwater issue. In this study, as a direct approach, the interaction between the groundwater level and the lake water level is investigated both in time and space by analysing the groundwater data compiled from observation wells surrounding the lake. Baseflow separation is considered as an indirect approach to understand the groundwater contribution to the river network flowing into the lake. It is determined that about 70% of run‐off in the rivers draining into the lake is born in the form of baseflow. An interaction between the lake and the groundwater storage is clearly seen from the analysis to conclude that groundwater has a potential to recharge the lake. Thus, the shrinkage of water level in Lake Urmia could be expected to accelerate with the drastic use of groundwater, which will be a disaster with no return.     

Determination of groundwater discharge rates and water residence time of groundwater‐fed lakes by stable isotopes of water (18O, 2H) and radon (222Rn) mass balances

Lacustrine groundwater discharge (LGD) and the related water residence time are crucial parameters for quantifying lake matter budgets and assessing its vulnerability to contaminant input. Our approach utilizes the stable isotopes of water (δ¹⁸O, δ²H) and the radioisotope radon (²²²Rn) for determining long‐term average and short‐term snapshots in LGD. We conducted isotope balances for the 0.5‐km² Lake Ammelshainer See (Germany) based on measurements of lake isotope inventories and groundwater composition accompanied by good quality and comprehensive long‐term meteorological and isotopic data (precipitation) from nearby monitoring stations. The results from the steady‐state annual isotope balances that rely on only two sampling campaigns are consistent for both δ¹⁸O and δ²H and suggested an overall long‐term average LGD rate that was used to infer the water residence time of the lake. These findings were supported by the good agreement of the simulated LGD‐driven annual cycles of δ¹⁸O and δ²H lake inventories with the observed lake isotope inventories. However, radon mass balances revealed lower values that might be the result of seasonal LGD variability. For obtaining further insights into possible seasonal variability of groundwater–lake interaction, stable water isotope and radon mass balances could be conducted more frequently (e.g., monthly) in order to use the derived groundwater discharge rates as input for time‐variant isotope balances.     

Multiple‐indicator study of groundwater flow and chemistry and the impacts of river and paddy water on groundwater in the alluvial fan of the Tedori River,Japan

To investigate the source, flow paths, and chemistry of rich resources of high‐quality, shallow groundwater in the alluvial fan between the Tedori and Sai rivers in central Japan, we analysed stable isotope ratios of H, O, and Sr and concentrations of major dissolved ions and trace elements in groundwater, river water, and paddy water. The ⁸⁷Sr/⁸⁶Sr ratios of the groundwater are related to near‐surface geology: groundwater in sediment from the Tedori River has high ⁸⁷Sr/⁸⁶Sr ratios (>0.711), whereas that from the Sai River in the north of the fan has low ⁸⁷Sr/⁸⁶Sr ratios (<0.711). δ²H and δ¹⁸O values and ⁸⁷Sr/⁸⁶Sr ratios indicate that groundwater in the central and southern fans is recharged by the Tedori River, whereas recharge in the north is from the Sai River. Mg²⁺, Ca²⁺, Sr²⁺, HCO₃^?, and SO₄^2? concentrations and δ²H and δ¹⁸O values in the groundwater are high in the central fan and, except for the northern area, tend to increase with distance from the Tedori River. There are linear relationships between ⁸⁷Sr/⁸⁶Sr ratio and the reciprocal concentrations of Sr²⁺, Mg²⁺, and Ca²⁺. These geochemical characteristics suggest that as groundwater recharged from the Tedori River flows towards the central fan, it mixes with waters derived from precipitation and paddy water that have become enriched in these components during downward infiltration. These results are consistent with our hydrological analysis and numerical simulation of groundwater flow, thus verifying the validity of the model we used in our simulation of groundwater flow. Copyright © 2016 John Wiley & Sons, Ltd.     

The importance of incorporating diurnally fluctuating stream discharge in stream temperature energy balance models

Although stream temperature energy balance models are useful to predict temperature through time and space, a major unresolved question is whether fluctuations in stream discharge reduce model accuracy when not exactly represented. However, high‐frequency (e.g., subdaily) discharge observations are often unavailable for such simulations, and therefore, diurnal streamflow fluctuations are not typically represented in energy balance models. These fluctuations are common due to evapotranspiration, snow pack or glacial melt, tidal influences within estuaries, and regulated river flows. In this work, we show when to account for diurnally fluctuating streamflow. To investigate how diurnal streamflow fluctuations affect predicted stream temperatures, we used a deterministic stream temperature model to simulate stream temperature along a reach in the Quilcayhuanca Valley, Peru, where discharge varies diurnally due to glacial melt. Diurnally fluctuating streamflow was varied alongside groundwater contributions via a series of computational experiments to assess how uncertainty in reach hydrology may impact simulated stream temperature. Results indicated that stream temperatures were more sensitive to the rate of groundwater inflow to the reach compared with the timing and amplitude of diurnal fluctuations in streamflow. Although incorporating observed diurnal fluctuations in discharge resulted in a small improvement in model RMSE, we also assessed other diurnal discharge signals and found that high amplitude signals were more influential on modelled stream temperatures when the discharge peaked at specific times. Results also showed that regardless of the diurnal discharge signal, the estimated groundwater flux to the reach only varied from 1.7% to 11.7% of the upstream discharge. However, diurnal discharge fluctuations likely have a stronger influence over longer reaches and in streams where the daily range in discharge is larger, indicating that diurnal fluctuations in stream discharge should be considered in certain settings.     

How preferential flow delivers pre‐event groundwater rapidly to streams

Groundwater often accounts for a substantial fraction of flood hydrographs, but the processes responsible for this have been unclear. However, many aquifers have preferential flow and this explains how aquifers can be so responsive. In bedrock aquifers, weathering enhances the connectivity and apertures along the most efficient flow paths and hence enhances the permeability. This results in celerities and velocities of the preferential flow in these dual‐porosity aquifers that are two to three orders of magnitude higher than if the aquifers behaved as single‐porosity media. The celerities have been determined from artificial and natural flood pulses, from tidal lags, and from pumping tests. Preferential‐flow velocities have been calculated from tests using applied tracers. Celerities in bedrock aquifers are typically one to two orders of magnitude faster than velocities. The ubiquitous preferential flow in aquifers provides an additional explanation, besides groundwater ridging, for the rapid release of groundwater to streams during storm events.     

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.     

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

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

Coupled groundwater flow and heat transport simulation for estimating transient aquifer–stream exchange at the lowland River Spree (Germany)

Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10^?7 to 10^?6 ms^?1, mean groundwater exfiltration rates of 233 l m^?2 d^?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 John Wiley & Sons, Ltd.     

A tracer‐based simulation approach to quantify seasonal dynamics of surface‐groundwater interactions in the Pantanal wetland

The Pantanal wetland is one of the least explored regions of South America. It is characterized by an outstanding flora and fauna adapted to a seasonal flood pulse controlled by a dry and a wet season within each year. The resulting inundation covers in average an area of approximately 150 000 km² and is seen as the most important driver for ecological integrity. Evaporation from the large floodplain is supposed to influence the climate of the whole continent. The regional groundwater is connected to the surface water and plays an important role for the characteristic flooding regime by regulating the wetland's water table. The water balance assessment of the wetland and the internal water exchange between surface and groundwater is therefore of high relevance for the conservation of the Pantanal biodiversity. Despite of its importance, water balance studies including groundwater–surface water interactions based on field data are rarely undertaken. This is mainly due to the remoteness and difficulty in accessing this area, which results in lack of data. In our study, we developed a new tracer‐based model to simulate the spatio–temporal surface and subsurface fluxes for a range of water bodies. The model was able to simulate these fluxes considering a dynamic simulation of inflow and outflow using a newly collected 2‐year dataset of water levels, stable water isotopes and chloride collected from several water bodies in the northern Pantanal region. Quantitative differences between water bodies according to their location in the floodplain were determined by the flooding regime and connectivity as well as site‐specific characteristics, such as hydraulic conductivity and water depth. Our model simulated water balance fluxes with a Nash–Sutcliffe efficiency of 0.61, whereas it simulated stable water isotopic compositions better than chloride. We present the first study based on field data for the Pantanal, which is able to quantify water balances fluxes. Because their representation in global climate and land cover products is insufficient, our simulation results are valuable for validating large‐scale models. Copyright © 2016 John Wiley & Sons, Ltd.     

Air/water/sediment temperature contrasts in small streams to identify groundwater seepage locations

The need to identify groundwater seepage locations is of great importance for managing both stream water quality and groundwater sourced ecosystems due to their dependency on groundwater‐borne nutrients and temperatures. Although several reconnaissance methods using temperature as tracer exist, these are subjected to limitations related to mainly the spatial and temporal resolution and/or mixing of groundwater and surface water leading to dilution of the temperature differences. Further, some methods, for example, thermal imagery and fiber optic distributed temperature sensing, although relative efficient in detecting temperature differences over larger distances, these are labor‐intensive and costly. Therefore, there is a need for additional cost‐effective methods identifying substantial groundwater seepage locations. We present a method expanding the linear regression of air and stream temperatures by measuring the temperatures in dual‐depth; in the stream column and at the streambed‐water interface (SWI). By doing so, we apply metrics from linear regression analysis of temperatures between air/stream and air/SWI (linear regression slope, intercept, and coefficient of determination), and the daily water temperature cycle (daily mean temperatures, temperature variance, and the mean diel temperature fluctuation). We show that using metrics from only single‐depth stream temperature measurements are insufficient to identify substantial groundwater seepage locations in a head‐water stream. Conversely, comparing the metrics from dual‐depth temperatures show significant differences; at groundwater seepage locations, temperatures at the SWI merely explain 43–75% of the variation opposed to ? 91% at the corresponding stream column temperatures. In general, at these locations at the SWI, the slopes ( < 0.25) and intercepts ( > 6.5 °C) are substantially lower and higher, respectively, while the mean diel temperature fluctuations ( < 0.98 °C) are decreased compared to remaining locations. The dual‐depth approach was applied in a post‐glacial fluvial setting, where metrics analyses overall corroborated with field measurements of groundwater fluxes and stream flow accretions. Thus, we propose a method reliably identifying groundwater seepage locations along streambeds in such settings.     

Combining isotopic tracers (222Rn and δ13C) for improved modelling of groundwater discharge to small rivers

In regions where aquifers sustain rivers, the location and quantification of groundwater discharge to surface water are important to prevent pollution hazards, to quantify and predict low flows and to manage water supplies. ²²²Rn is commonly used to determine groundwater discharge to rivers. However, using this isotopic tracer is challenging because of the high diffusion capacity of ²²²Rn in open water. This study illustrates how a combination of isotopic tracers can contribute to an enhanced understanding of groundwater discharge patterns in small rivers. The aim of this paper is to combine ²²²Rn and δ¹³C_DIC to better constrain the physical parameters related to the degassing process of these tracers in rivers. The Hallue River (northern France) was targeted for this study because it is sustained almost exclusively by a fractured chalk aquifer. The isotopes ²²²Rn, δ¹³C_DIC, δ²H and δ¹⁸O were analysed along with other natural geochemical tracers. A mass balance model was used to simulate ²²²Rn and δ¹³C_DIC. The results of δ²H and δ¹⁸O analyses prove that evaporation did not occur in the river. The calibration of a numerical model to reproduce ²²²Rn and δ¹³C_DIC provides a best‐fit diffusive layer thickness of 3.21 × 10^?5 m. This approach is particularly useful for small rivers flowing over carbonate aquifers with high groundwater DIC where the evolution of river DIC reflects the competing processes of groundwater inflow and CO₂ degassing. This approach provides a means to evaluate groundwater discharge in small ungauged rivers. Copyright © 2014 John Wiley & Sons, Ltd.     

Field assessment of surface water–groundwater connectivity in a semi‐arid river basin (Murray–Darling,Australia)

In semi‐arid and arid river basins, understanding the connectivity between rivers and alluvial aquifers is one of the key challenges for the management of groundwater resources. The type of connection present (gaining, losing‐connected, transitional and losing‐disconnected) was assessed at 12 sites along six Murray–Darling Basin river reaches. The assessments were made by measuring the hydraulic head in the riparian zone near the rivers to evaluate if the water tables intersected the riverbeds and by measuring fluid pressure (ψ) in the riverbeds. The rationale for the latter was that ψ will always be greater than or equal to zero under connected conditions (either losing or gaining) and always lesser than or equal to zero under losing‐disconnected conditions. A mixture of losing‐disconnected, losing‐connected and gaining conditions was found among the 12 sites. The losing‐disconnected sites all had a riverbed with a lower hydraulic conductivity than the underlying aquifer, usually in the form of a silty clay or clay unit 0.5–2 m in thickness. The riparian water tables were 6 to 25 m below riverbed level at the losing‐disconnected sites but never lower than 1 m below riverbed level at the losing‐connected ones. The contrast in water table depth between connected and disconnected sites was attributed to the conditions at the time of the study, when a severe regional drought had generated a widespread decline in regional water tables. This decline was apparently compensated near losing‐connected rivers by increased infiltration rates, while the decline could not be compensated at the losing‐disconnected rivers because the infiltration rates were already maximal there. Copyright © 2012 John Wiley & Sons, Ltd.     

Directions of preferential flow in a hillslope soil, 1. Quasi‐steady flow

Preferred infiltration is mainly perceived as vertically down whereas subsurface storm flow is thought to occur parallel to slopes. The transition from vertical to lateral flow in a layered hillslope soil is the focus of the contribution. Transient flow is assumed to move as a wetting front. Three time‐domain reflectometry (TDR) wave‐guides, each 0·15 m long, were mounted in the shape of a truncated tetrahedron with its peak pointing down. Each wave‐guide focuses the front velocity along its axis. The three front‐velocity vectors are decomposed into their x, y and z components, which are then assembled to the resultant velocity vector. The volume density flux of preferred flow is the product of the front velocity and the mobile water content. The latter is the amplitude of transient soil moisture measured with each wave‐guide. The resultant vector of the volume flux density is computed similarly to the velocity vector. The experimental approach allows for the rapid assessment of transient flows without relying on the variation of water potentials. The experiments indicate that the directions of the resultant vectors of velocity and volume flux density can be estimated if the moisture variations of the three TDR wave‐guides are strongly correlated during the passing of the wetting front. Copyright © 2004 John Wiley & Sons, Ltd.     

Net groundwater inflow in an enclosed lake: from synoptic variations to climatic projections

Using lake Stechlin in northeastern Germany as an example of a small groundwater‐feed lake without surface inflows and outflows, we estimated the temporal scales and the variability ranges of the net groundwater contribution to the lake water budget. High‐resolution water level measurements by a bottom‐mounted pressure logger provided the background for the estimation of the total lake water budget. This method has demonstrated reliability for estimation of lake level variations during periods ranging from subdiurnal to perennial. The typical amplitudes of the synoptic‐to‐perennial variability characterizing the groundwater climate of lake Stechlin are estimated by comparing the two subsequent years 2006 and 2007; one of these years shows an extremely high, and the other an extremely low, annual precipitation–evaporation balance. The net groundwater flow, estimated as the difference between the total water budget and the precipitation–evaporation balance at the surface, revealed synoptic effects of lake water exfiltration into the groundwater aquifer following strong precipitation events. Perennial variations between wet and dry years superimposed seasonal oscillations. The probable origin of the latter is seasonality in the groundwater level on the watershed, although the exact amplitudes are subject to further quantification on account of seasonality in the evaporation estimation error. The results emphasize the non‐stationary behaviour of groundwater flow on timescales shorter than climatic ones. The analysis yielded a net quantitative relationship between groundwater flow and water balance at the lake surface: The water level changes in the lake due to evaporation and precipitation are damped to 60% because of the lake–groundwater exchange by means of intermittent infiltration and exfiltration events. Assuming the remaining 40% of the surface water budget may potentially result in perennial water level variability, we estimated an effect of the precipitation decrease on the lake water budget as predicted by the regional climate scenarios for the next century. Copyright © 2012 John Wiley & Sons, Ltd.     

Water level observations from unmanned aerial vehicles for improving estimates of surface water–groundwater interaction

Integrated hydrological models are usually calibrated against observations of river discharge and piezometric head in groundwater aquifers. Calibration of such models against spatially distributed observations of river water level can potentially improve their reliability and predictive skill. However, traditional river gauging stations are normally spaced too far apart to capture spatial patterns in the water surface, whereas spaceborne observations have limited spatial and temporal resolution. Unmanned aerial vehicles can retrieve river water level measurements, providing (a) high spatial resolution; (b) spatially continuous profiles along or across the water body, and (c) flexible timing of sampling. A semisynthetic study was conducted to analyse the value of the new unmanned aerial vehicle‐borne datatype for improving hydrological models, in particular estimates of groundwater–surface water (GW–SW) interaction. Mølleåen River (Denmark) and its catchment were simulated using an integrated hydrological model (MIKE 11–MIKE SHE). Calibration against distributed surface water levels using the Differential Evolution Adaptive Metropolis algorithm demonstrated a significant improvement in estimating spatial patterns and time series of GW–SW interaction. After water level calibration, the sharpness of the estimates of GW–SW time series improves by ~50% and root mean square error decreases by ~75% compared with those of a model calibrated against discharge only.     

Tracer and hydrometric study of preferential flow in large undisturbed soil cores from the Georgia Piedmont,USA

We studied the temporal patterns of tracer throughput in the outflow of large (30 cm diameter by 38 cm long) undisturbed cores from the Panola Mountain Research Watershed, Georgia. Tracer breakthrough was affected by soil structure and rainfall intensity. Two rainfall intensities (20 and 40 mm hr⁻¹) for separate Cl⁻ and Br⁻ amended solutions were applied to two cores (one extracted from a hillslope soil and one extracted from a residual clay soil on the ridge). For both low and high rainfall intensity experiments, preferential flow occurred in the clay core, but not in the hillslope core. The preferential flow is attributed to well‐developed interpedal macrochannels that are commonly found in structured clay soils, characteristic of the ridge site. However, each rainfall intensity exceeded the matrix infiltration capacity at the top of the hillslope core, but did not exceed the matrix infiltration capacity at the middle and bottom of the hillslope core and at all levels in the clay core. Localized zones of saturation created when rainfall intensity exceeds the matrix infiltration capacity may cause water and tracer to overflow from the matrix into macrochannels, where preferential flow occurs to depth in otherwise unsaturated soil. Copyright © 1999 John Wiley & Sons, Ltd.     

Groundwater level responses in temperate mountainous terrain: regime classification,and linkages to climate and streamflow

Groundwater responses in temperate mountainous terrain are assessed using groundwater, hydrometric and climatic data from southern British Columbia, Canada. Well and stream hydrographs are analysed using a series of diagnostic tools including time series plots, hysteresis plots, and cross‐correlation plots. Characterizing the seasonal timing of the response requires consideration of the hydroclimatology of the region: rainfall‐dominated (pluvial), snowmelt‐dominated (nival) or hybrid (mixture of rain and snow). The magnitude and timing of the recharge and discharge response of the groundwater system was shown to depend on the storage and permeability characteristics of the aquifer and whether the system is stream‐driven or recharge‐driven. These two dominant stream‐aquifer system types were defined based on classifying different aquifer types found in the southwest portion of the province. The classification scheme and diagnostic tools have the potential to provide a framework for evaluating the responses of wells in other mountainous regions. Using this framework, the potential consequences of future climate change may then be better understood based on the interactions between the hydrogeological and hydroclimatic settings of these aquifers. Copyright © 2010 Her Majesty the Queen in right of Canada. Published by John Wiley & Sons. Ltd     

Hydraulic and geochemical interactions between surface water and sediment pore water in seasonal hypersaline Maharlu Lake,Iran

Study of interactions between surface-water and pore-water in lakes is complicated due to spatio-temporal heterogeneities in flow condition across the sediment–water interface. In this study, seasonal hypersaline Maharlu Lake was investigated by collecting surface-water and pore-water samples from four nests of multilevel piezometers installed at different distances from the inflow of rivers to the lake. The hydraulic heads in the piezometers as well as vertical profiles of Mg⁺², Na/Cl, and Br/Cl were used to investigate both hydraulic and geochemical interactions between surface-water and pore-water in the lake. Depletion of lake surface water and pore water with respect to B, Br, Li⁺, K⁺, Mg²⁺ and the absence of Mg-K chlorides and sulphates in the lake bed sediments is probably due to leakage of highly evaporated residual brine from the lake. Hydraulic gradients in the multilevel piezometric nests indicate that a general downward flow from surface-water to pore-water occurs across sediment–water interface. Vertical profiles of Br/Cl, Mg²⁺, and Na/Cl showed that the maximum flow rate was more than 1 m/yr close to the mouth of the inflowing rivers. The downward vertical flow was limited in the area far from the inflowing rivers due to the presence of an impermeable confining halite layer which interrupts the hydraulic connection between shallow pore water (less than 50 cm deep) and deeper zones. The hydraulic and geochemical interactions between surface-water and pore-water across sediment–water interface in the Maharlu Lake are of interest to find out the fate of pollutants and their distribution in the lake.     

Geochemical evidence of hydrothermal recharge in Lake Baringo,central Kenya Rift Valley

Lake Baringo, a freshwater lake in the central Kenya Rift Valley, is fed by perennial and ephemeral rivers, direct rainfall, and hot springs on Ol Kokwe Island near the centre of the lake. The lake has no surface outlet, but despite high evaporation rates it maintains dilute waters by subsurface seepage through permeable sediments and faulted lavas. New geochemical analyses (major ions, trace elements) of the river, lake, and hot spring waters and the suspended sediments have been made to determine the main controls of lake water quality. The results show that evaporative concentration and the binary mixing between two end members (rivers and thermal waters) can explain the hydrochemistry of the lake waters. Two zones are recognized from water composition. The southern part of the lake near sites of perennial river inflow is weakly influenced by evaporation, has low total dissolved species (TDS), and has a seasonally variable load of mainly detrital suspended sediments. In contrast, waters of the northern part of the lake show evidence for strong evaporation (TDS of up to eight times inflow). Authigenic clay minerals and calcite may be precipitating from those more concentrated fluids. The subaerial hot‐spring waters have a distinctive chemistry and are enriched in some elements that are also present in the lake water. Comparison of the chemical composition of the inflowing surface waters and lake water shows (1) an enrichment of some species (HCO₃^?, Cl, SO₄^2?, F, Na, B, V, Cr, As, Mo, Ba and U) in the lake, (2) a depletion in SiO₂ in the lake, and (3) a possible hydrothermal origin for most F. The rare earth element distribution and the F/Cl and Na/Cl ratios give valuable information on the rate of mixing of the river and hydrothermal fluids in the lake water. Calculations imply that thermal fluids may be seeping upward locally into the lake through grid‐faulted lavas, particularly south of Ol Kokwe Island. Copyright © 2006 John Wiley & Sons, Ltd.     

Assessing the importance of seepage and springs to nitrate flux in a stream network in the Wisconsin sand plains

Evaluating the flow paths that contribute to solute flux in stream networks can lead to greater understanding of the linkages between biogeochemistry and hydrology. We compared the contributions of groundwater in spring brooks and in seepage through the streambed to nitrate flux in the Emmons Creek network in the Wisconsin sand plains. We predicted that spring brooks would contribute disproportionately to nitrate flux due to the presumed higher advection rates in springs and less opportunity for nitrate removal relative to seeps. Nitrate flux was measured in 15 spring brooks that entered Emmons Creek. Nitrate flux from seepage was measured at the locations of 30 piezometers, based on Darcy's Law, and by a reach‐scale injection of Rhodamine water tracing (RWT). When seepage discharge was estimated from the RWT release, groundwater inputs from seepage and springs accounted for the discharge gain in the Emmons Creek channel. Springs brooks and seepage (based on the RWT release) contributed 37% and 63%, respectively, to nitrate flux inputs in the study reach. Contrary to our prediction, seeps contributed disproportionately to nitrate flux relative to their discharge. Relatively high rates of seepage discharge and higher than anticipated nitrate concentrations in the shallow pore water at seepage locations contributed to the unanticipated result.