Impacts of riparian zone plant water use on temporal scaling of groundwater systems

In this work, we study groundwater system temporal scaling in relation to plant water use and near‐river‐stage fluctuations in riparian zones where phreatophytes exist. Using detrended fluctuation analysis (DFA), we investigate the influence of regular diurnal fluctuations due to phreatophyte water use on temporal scaling properties of groundwater level variations. We found that groundwater use by phreatophytes, at the field site on the Colorado River, USA, results in distinctive crossovers (slope changes when the plots are fitted with straight lines) in the logarithm plots of root‐mean‐square fluctuations of the detrended water level time series versus time scales of groundwater level dynamics. For groundwater levels monitored at wells close to the river, we identified one crossover at ～1 day in the scaling characteristics of groundwater level variations. When time scale exceeds 1 day, the scaling properties decrease from persistent to close to 1/f noise, where f is the frequency. For groundwater levels recorded at wells further away from the river, the slope of the straight line fit (i.e. scaling exponent) is smallest when the time scale is between 1 and 3 days. When the time scale is < 1 day, groundwater variations become persistent. When the time scale is between 1 and 3 days, the variations are close to white noise, but return to persistent when the time scale is > 3 days. Copyright © 2011 John Wiley & Sons, Ltd.     

Vertical surface water–groundwater exchange processes within a headwater floodplain induced by experimental floods

Restoring hydrologic connectivity between channels and floodplains is common practice in stream and river restoration. Floodplain hydrology and hydrogeology impact stream hydraulics, ecology, biogeochemical processing, and pollutant removal, yet rigorous field evaluations of surface water–groundwater exchange within floodplains during overbank floods are rare. We conducted five sets of experimental floods to mimic floodplain reconnection by pumping stream water onto an existing floodplain swale. Floods were conducted throughout the year to capture seasonal variation and each involved two replicate floods on successive days to test the effect of varying antecedent moisture. Water levels and specific conductance were measured in surface water, soil, and groundwater within the floodplain, along with surface flow into and out of the floodplain. Vegetation density varied seasonally and controlled the volume of surface water storage on the floodplain. By contrast, antecedent moisture conditions controlled storage of water in floodplain soils, with drier antecedent moisture conditions leading to increased subsurface storage and slower flood wave propagation across the floodplain surface. The site experienced spatial heterogeneity in vertical connectivity between surface water and groundwater across the floodplain surface, where propagation of hydrostatic pressure, preferential flow, and bulk Darcy flow were all mechanisms that may have occurred during the five floods. Vertical connectivity also increased with time, suggesting higher frequency of floodplain inundation may increase surface water–groundwater exchange across the floodplain surface. Understanding the variability of floodplain impacts on water quality noted in the literature likely requires better accounting for seasonal variations in floodplain vegetation and antecedent moisture as well as heterogeneous exchange flow mechanisms. Copyright © 2016 John Wiley & Sons, Ltd.     

A numerical model for estimating groundwater flux from subsurface temperature profiles

A simple numerical model is presented for estimating vertical groundwater flux from transient subsurface temperature profiles obtained from field measurements. The model developed utilizes the MacCormack scheme, which is based on the Finite Difference Method (FDM), for solving the governing partial differential equation of convection–diffusion heat transport with appropriate initial and boundary conditions within the subsurface. In order to validate the model, numerical solutions obtained for the study area located in the Nagoka plain, Japan are compared with the published measured data and results obtained by others. Results obtained show good agreement and fit the observed data with a correlation coefficient, R², of 0·88. The estimated groundwater flux is 1·85 × 10⁻⁷ m s⁻¹. Sensitivity analyses were also carried out to investigate the effect of variations in groundwater fluxes, thermal properties and the annual thermal variability due to climatic changes on the transient subsurface temperature profiles and to have a better understanding of the subsurface thermal dynamics. A substantial effect of annual climatic variability is observed on the temporal distributions of temperature depth profiles, and a better estimate of thermal parameters is required to estimate vertical groundwater flux. The largest change in subsurface temperature depth profiles due to groundwater flux over a year is within ± 4 °C. The influence of groundwater flux on subsurface temperature distributions in space and time may be more pronounced in areas where the top of the saturated layer fluctuates considerably. Variation in thermal diffusivity results in temperature change up to ± 1·5% and may cause change in groundwater flux estimate by ± 18%. The model presented has merits over analytical solutions (type curve matching techniques) in terms of suitability and applicability to real field problems, and can be a good asset to hydrological models as quantifying groundwater recharge or deducing it from other quantities, such as rainfall, evapotranspiration and runoff, is often complicated. Copyright © 2007 John Wiley & Sons, Ltd.     

Interactions between dam‐regulated river flow and riparian groundwater: a case study from the Yellow River,China

The Xiaolangdi Dam, completed in 2000, is second in scale in China to the Three Gorges Project. It has generated remarkable economic and social benefits but with profound impacts to the riverine and regional environments. This paper reports field monitoring of riparian groundwater in the Kouma section of the Yellow River to illustrate the interactions between dam‐regulated river flow and riparian groundwater. The results show that the hydrological condition in riparian zones downstream from the dam has changed from a typical wet–dry cycle to a condition of semi‐permanent dryness, resulting in degradation of the typical attributes and functions of the wetland ecosystem. Hydrological processes in the riparian zone have changed from a complex multiple flooding regime to a simple regime of dominant groundwater drainage towards the river, which only reverses temporarily during the water and sediment regulation period of the dam. Data on groundwater level and groundwater quality show that there are two key points, at ca 200 and 400 m from the river bank, which distinguish zones with different sensitivity to changes of river flow and indicate different interactions between river water and groundwater. The shallow groundwater quality also is negatively affected by the intensive agricultural development that has occurred since the dam was completed. Ecological restoration needs to be carried out to construct a protective natural riparian zone within ca 200 m from the river, this being an ecotone, which is key to the protection of both riparian groundwater and the river. The riparian zone from 200 to 400 m also should be treated as a transitional zone. In addition, ecologically sensitive agriculture and ecotourism organized by local communities would be beneficial in the area beyond 400 m. Copyright © 2011 John Wiley & Sons, Ltd.     

Denitrification in a hyporheic riparian zone controlled by river regulation in the Seine river basin (France)

The purpose of this paper is to study denitrification and the conditions for its development in a hyporheic zone. The study site is the riparian zone of a former branch of the Seine River, where the river stage is kept almost constant during the year by hydraulic regulation. Hydrological and geochemical surveys were performed by monitoring four wells, ten shorter piezometers and the river over a 15‐month period. The water fluxes originating from the chalky hillsides and the river converge in a zone parallel to the river that acts as a drainage flow path through the floodplain. The riparian zone between this flow path and the river shows an important depletion of nitrate during the summer and autumn period, which cannot be explained by a simple mixing of waters coming from the river and the chalky hillsides. It can be attributed to denitrification as it occurs when oxygen concentration is below 2 mg l^?1, and goes along with a consumption of dissolved organic carbon and a decrease of redox potential. The river completely controls these hydro‐geochemical conditions. It also keeps the wetness of the riparian zone almost constant, which allowed us to isolate the high temperatures in summer and autumn as an important triggering factor for denitrification through its influence on the reaction rate and oxygen deficits. We also found a small isotopic enrichment of nitrate, suggesting that denitrification occurs after diffusion of nitrate through the sediment and riparian zone matrix, which is consistent with the hyporheic functioning of the study site. Copyright © 2008 John Wiley & Sons, Ltd.     

Evaluating the reliability of time series analysis to estimate variable riparian travel times by numerical groundwater modelling

The transition zones between rivers and adjacent riparian aquifers are locations of high biogeochemical activities that contribute to a removal of potentially hazardous substances in the aquatic system. The potential of the removal processes depends highly on subsurface water travel times, which can be determined by using the propagation of electrical conductivity (EC) signal from the river into the riparian aquifer. Although this method has been applied and verified in many studies, we observe possible limitations for the usage of EC fluctuation analysis. Our findings are based on EC time series analyses during storm events and artificial hydropeaks induced by watermill operations. Travel times derived by cross‐correlation analysis were compared with travel times calculated based on backward particle tracking of a calibrated transient numerical groundwater flow model. The cross‐correlation method produced only reasonable travel times for the artificial hydropeaks. In contrast, cross‐correlation analysis of the EC data during natural storm events resulted in implausibly negative or unrealistically low travel times for the bulk of the data sets. We conclude that the reason for this behaviour is, first, the low EC contrast between river and groundwater in connection with a strong damping of the infiltrating river EC signal into the subsurface during storm events. Second, the existence of old and less‐mineralized riparian water between the river and the monitoring well resulted in bank‐storage‐driven EC breakthrough curves with earlier arrival times and the subsequent estimation of implausible riparian travel times.     

Woody buffer effects on water temperature: The role of spatial configuration and daily temperature fluctuations

Water temperature is a key driver for riverine biota and strongly depends on shading by woody riparian vegetation in summer. While the general effects of shading on daily maximum water temperature T_max are well understood, knowledge gaps on the role of the spatial configuration still exist. In this study, the effect of riparian buffer length, width, and canopy cover (percentage of buffer area covered by woody vegetation) on T_max was investigated during summer baseflow using data measured in seven small lowland streams in western Germany (wetted width 0.8–3.7 m). The effect of buffer length on T_max differed between downstream cooling and heating: T_max approached cooler equilibrium conditions after a distance of 0.4 km (~45 min travel-time) downstream of a sharp increase in canopy cover. In contrast, T_max continued to rise downstream of a sharp decrease in canopy cover along the whole 1.6 km stream length investigated. The effect of woody vegetation on T_max depended on buffer width, with changes in canopy cover in a 10 m wide buffer being a better predictor for changes in T_max compared to a 30 m buffer. The effect of woody vegetation on T_max was linearly related to canopy cover but also depended on daily temperature range T_range, which itself was governed by cloudiness, upstream canopy cover, and season. The derived empirical relationship indicated that T_max was reduced by −4.6°C and increased by +2.7°C downstream of a change from unshaded to fully shaded conditions and vice versa. This maximum effect was predicted for a 10 m wide buffer at sunny days in early summer, in streams with large diel fluctuations (large T_range). Therefore, even narrow woody riparian buffers may substantially reduce the increase in T_max due to climate change, especially in small shallow headwater streams with low baseflow discharge and large daily temperature fluctuations.     

Travel time controls the magnitude of nitrate discharge in groundwater bypassing the riparian zone to a stream on Virginia's coastal plain

Groundwater that bypasses the riparian zone by travelling along deep flow paths may deliver high concentrations of fertilizer‐derived NO₃^? to streams, or it may be impacted by the NO₃^? removal process of denitrification in streambed sediments. In a study of a small agricultural catchment on the Atlantic coastal plain of Virginia's eastern shore, we used seepage meters deployed in the streambed to measure specific discharge of groundwater and its solute concentrations for various locations and dates. We used values of Cl^? concentration to discriminate between bypass water recharged distal to the stream and that contained high NO₃^? but low Cl^? concentrations and riparian‐influenced water recharged proximal to the stream that contained low NO₃^? and high Cl^? concentrations. The travel time required for bypass water to transit the 30‐cm‐thick, microbially active denitrifying zone in the streambed determined the extent of NO₃^? removal, and hydraulic conductivity determined travel time through the streambed sediments. At all travel times greater than 2 days, NO₃^? removal was virtually complete. Comparison of the timescales for reaction and transport through the streambed sediments in this system confirmed that the predominant control on nitrate flux was travel time rather than denitrification rate coefficients. We conclude that extensive denitrification can occur in groundwater that bypasses the riparian zone, but a residence time in biologically active streambed sediments sufficient to remove a large fraction of the NO₃^? is only achieved in relatively low‐conductivity porous media. Instead of viewing them as separate, the streambed and riparian zone should be considered an integrated NO₃^? removal unit. Copyright © 2011 John Wiley & Sons, Ltd.     

Identification of groundwater exfiltration,interflow discharge,and hyporheic exchange flows by fibre optic distributed temperature sensing supported by electromagnetic induction geophysics

Water exchange across the sediment–water interface of streams impresses a characteristic thermal pattern at the interface. The use of fibre optic distributed temperature sensing at the sediment–water interface in a small sand‐bed stream identifies such temperature patterns. Groundwater and interflow can be differentiated based on the temporal evolution of temperature patterns. Additionally, sudden temperature changes at the sediment–water interface observed during the transit of floods enable spatial identification of local up and downwelling. Electromagnetic induction geophysics can detect subsurface texture structures that support groundwater–surface water exchange. Our results show that areas of permanent temperature anomalies observed with fibre optic distributed temperature sensing match areas of comparatively homogeneous electrical conductivity. This indicates groundwater discharge and enables differentiating groundwater discharge from interflow and local downwelling.     

Riparian influence on hyporheic‐zone formation downstream of a small dam in the Blackland Prairie region of Texas

Small‐order streams have highly variable flows that can result in large temporal and spatial variation of the hyporheic zone. Dam construction along these intermittent headwater streams alters downstream flow and influences the hydrologic balance between stream water and the adjacent riparian zone. A 3‐year site study was conducted along an impounded second‐order stream to determine the water balance between stream, unsaturated zone, groundwater and riparian vegetation. The presence of the upstream impoundment provided near‐perennial water flow in the stream channel. The observed woody plant transpiration accounted for 71% of average annual water loss in the site. The overall contribution of stream water via the hyporheic zone to site water balance was 73 cm, or 44% of total inputs. This exceeded both rainfall and upland subsurface contribution to the site. A highly dynamic hyporheic zone was indicated by high water use from woody plants that fluctuated seasonally with stream water levels. We found leaf area development in the canopy layer to be closely coupled with stream and groundwater fluctuations, indicating its usefulness as a potential indicator of site water balance for small dam systems. The net result of upstream impoundment increased riparian vegetation productivity by influencing movement of stream water to storage in the groundwater system. Copyright © 2006 John Wiley & Sons, Ltd.     

Spatial and temporal variations of nutrient concentration in the groundwater of a floodplain: effect of hydrology,vegetation and substrate

Spatio‐temporal variations in nitrogen and phosphorus concentrations in groundwater were analysed and related to the variations in hydrological conditions, vegetation type and substrate in an alluvial ecosystem. This study was conducted in the Illwald forest in the Rhine Plain (eastern France) to assess the removal of nutrients from groundwater in a regularly flooded area. We compared both forest and meadow ecosystems on clayey‐silty soils with an anoxic horizon (pseudogley) at 1·5–2 m depth (eutric gley soil) and a forest ecosystem on a clayey‐silty fluviosoil rich in organic matter with a gley at 0·5 m depth (calcaric gley soil). Piezometers were used to measure the nutrient concentrations in the groundwater at 2 m depth in the root layer and at 4·5 m depth, below the root layer. Lower concentrations of nitrate and phosphate in groundwater were observed under forest than under meadow, which could be explained by more efficient plant uptake by woody species than herbaceous plants. Thus NO₃‐N inputs by river floods were reduced by 73% in the shallow groundwater of the forested ecosystem, and only by 37% in the meadow. Compared with the superficial groundwater layer, the lowest level of nitrate nitrogen (NO₃‐N) and the highest level of ammonium nitrogen (NH₄‐N) were measured in the deep layer (under the gley horizon at 2·5 m depth), which suggests that the reducing potential of the anoxic horizon in the gley soils contributes to the reduction of nitrate. Nitrate concentrations were higher in the groundwater of the parcel rich in organic matter than in the one poorer in organic matter. Phosphate (PO₄‐P) concentrations in both shallow and deep groundwater are less than 62 to 76% of those found in surface water which can be related to the retention capacity of the clay colloids of these soils. Moreover, the temporal variations in nutrient concentrations in groundwater are directly related to variations in groundwater level during an annual hydrological cycle. Our results suggest that variations in groundwater level regulate spatio‐temporal variations in nutrient concentrations in groundwater as a result of the oxidation–reduction status of soil, which creates favourable or unfavourable conditions for nutrient bioavailability. The hydrological variations are much more important than those concerning substrate and type of vegetation. Copyright © 1999 John Wiley & Sons, Ltd.     

Beaver dams along an agricultural stream in southern Ontario,Canada: their impact on riparian zone hydrology and nitrogen chemistry

The hydrology and nitrogen biogeochemistry of a riparian zone were compared before and after the construction of beaver dams along an agricultural stream in southern Ontario, Canada. The beaver dams increased surface flooding and raised the riparian water table by up to 1·0 m. Increased hydraulic gradients inland from the stream limited the entry of oxic nitrate‐rich subsurface water from adjacent cropland. Permeable riparian sediments overlying dense till remained saturated during the summer and autumn months, whereas before dam construction a large area of the riparian zone was unsaturated in these seasons each year. Beaver dam construction produced significant changes in riparian groundwater chemistry. Median dissolved oxygen concentrations were lower in riparian groundwater after dam construction (0·9–2·1 mg L^?1) than in the pre‐dam period (2·3–3·9 mg L^?1). Median NO₃‐N concentrations in autumn and spring were also lower in the post‐dam (0·03–0·07 mg L^?1) versus the pre‐dam period (0·1–0·3 mg L^?1). In contrast, median NH₄‐N concentrations in autumn and spring months were higher after dam construction (0·3–0·4 mg L^?1) than before construction (0·13–0·14 mg L^?1). Results suggest that beaver dams can increase stream inflow to riparian areas that limit water table declines and increase depths of saturated riparian soils which become more anaerobic. These changes in subsurface hydrology and chemistry have the potential to affect the transport and transformation of nitrate fluxes from adjacent cropland in agricultural landscapes. Copyright © 2009 John Wiley & Sons, Ltd.     

Agricultural and forested land use impacts on floodplain shallow groundwater temperature regime

Given the importance of groundwater temperature to the biogeochemical health of aquatic ecosystems, a floodplain study was implemented to improve understanding of rural land use impacts on shallow groundwater (SGW) temperature. Study sites included a historic agricultural field (Ag) and bottomland hardwood forest (BHF), each with nine piezometers in an 80 × 80 m grid. Piezometers were equipped with pressure transducers to monitor SGW temperature and level at 30 min intervals during the 2011, 2012, 2013, and 2014 water years. The study is one of the first to utilize long‐term, continuous, automated, in situ monitoring to investigate rural land use impacts on shallow groundwater temperatures. Average SGW temperature during the study period was 11.1 and 11.2 °C at the Ag and BHF sites, respectively. However, temperature range at the Ag site was 72% greater than at the BHF site. Results indicate a greater responsiveness to seasonal climate fluctuations in Ag site SGW temperature related to absence of forest canopy. Patterns of intra‐site groundwater temperature differences at both study sites illustrate the influence of stream–aquifer thermal conduction and occasional baseflow reversals. Considering similar surface soil temperature amplitudes and low average groundwater flow values at both sites, results suggest that contrasting rates of plant water use, groundwater recharge, and subsurface hydraulic conductivity are likely mechanistic causes for the observed SGW temperature differences. Results highlight the long‐term impact of forest removal on subsurface hydrology and groundwater temperature regime. Copyright © 2015 John Wiley & Sons, Ltd.     

Quantitative analysis of riverbank groundwater flow for the Qinhuai River,China, and its influence factors

Interactions of surface water and groundwater (SW–GW) play an important role in the physical, chemical, and ecological processes of riparian zones. The main objective of this study was to describe the two‐dimensional characteristics of riverbank SW–GW interactions and to quantify their influence factors. The SW–GW exchange fluxes for six sections (S1 to S6) of the Qinhuai River, China, were estimated using a heat tracing method, and field hydrogeological and thermodynamic parameters were obtained via inverse modelling. Global sensitivity analysis was performed to compare the effects of layered heterogeneity of hydraulic conductivity and river stage variation on SW–GW exchange. Under the condition of varied river stage, only the lateral exchange fluxes at S1 apparently decreased during the monitoring period, probably resulting from its relatively higher hydraulic conductivity. Meanwhile, the SW–GW exchanges for the other five sections were quite stable over time. The lateral exchange fluxes were higher than the vertical ones. The riverbank groundwater flow showed different spatial variation characteristics for the six sections, but most of the higher exchange fluxes occurred in the lower area of a section. The section with larger hydraulic conductivity has an apparent dynamic response to surface water and groundwater level differences, whereas lower permeabilities severely reduced the response of groundwater flow. The influence of boundary conditions on SW–GW interactions was restricted to a limited extent, and the impact extent will expand with the increase of peak water level and hydraulic conductivity. The SW–GW head difference was the main influence factors in SW–GW interactions, and the influence of both SW–GW head difference and hydraulic conductivity decreased with an increase of the distance from the surface water boundary. For each layer of riverbank sediment, its hydraulic conductivity had greater influence on its groundwater flow than the other layers, whereas it had negligible effects on its overlying/underlying layers. Consequently, the variations in river stage and hydraulic conductivity were the main factors influencing the spatial and temporal characteristics of riverbank groundwater flow, respectively.     

High frequency stream bed mobility of a low‐gradient agricultural stream with implications on the hyporheic zone

Little Kickapoo Creek (LKC), a low‐gradient stream, mobilizes its streambed–fundamentally altering its near‐surface hyporheic zone–more frequently than do higher‐gradient mountain and karst streams. LKC streambed mobility was assessed through streambed surveys, sediment sampling, and theoretical calculations comparing basal shear stress (τ_b) with critical shear stress (τ_c). Baseflow τ_b is capable of entraining a d₅₀ particle; bankfull flow could entrain a 51·2 mm particle. No particle that large occurs in the top 30 cm of the substrate, suggesting that the top 30 cm of the substrate is mobilized and redistributed during bankfull events. Bankfull events occur on average every 7·6 months; flows capable of entraining d₅₀ and d₈₅ particles occur on average every 0·85 and 2·1 months, respectively. Streambed surveys verify streambed mobility at conditions below bankfull. While higher gradient streams have higher potential energy than LKC, they achieve streambed‐mobilization thresholds less frequently. Heterogeneous sediment redistribution creates an environment where substrate hydraulic conductivity (K) varies over four orders of magnitude. The frequency and magnitude of the substrate entrainment has implications on hyporheic zone function in fluid, solute and thermal transport models, interpretations of hyporheic zone stability, and understanding of LKC's aquatic ecosystem. Copyright © 2008 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.     

Heat as a groundwater tracer in shallow and deep heterogeneous media: Analytical solution,spreadsheet tool,and field applications

Groundwater flow advects heat, and thus, the deviation of subsurface temperatures from an expected conduction‐dominated regime can be analysed to estimate vertical water fluxes. A number of analytical approaches have been proposed for using heat as a groundwater tracer, and these have typically assumed a homogeneous medium. However, heterogeneous thermal properties are ubiquitous in subsurface environments, both at the scale of geologic strata and at finer scales in streambeds. Herein, we apply the analytical solution of Shan and Bodvarsson ( 2004 ), developed for estimating vertical water fluxes in layered systems, in 2 new environments distinct from previous vadose zone applications. The utility of the solution for studying groundwater‐surface water exchange is demonstrated using temperature data collected from an upwelling streambed with sediment layers, and a simple sensitivity analysis using these data indicates the solution is relatively robust. Also, a deeper temperature profile recorded in a borehole in South Australia is analysed to estimate deeper water fluxes. The analytical solution is able to match observed thermal gradients, including the change in slope at sediment interfaces. Results indicate that not accounting for layering can yield errors in the magnitude and even direction of the inferred Darcy fluxes. A simple automated spreadsheet tool (Flux‐LM) is presented to allow users to input temperature and layer data and solve the inverse problem to estimate groundwater flux rates from shallow (e.g., <1 m) or deep (e.g., up to 100 m) profiles. The solution is not transient, and thus, it should be cautiously applied where diel signals propagate or in deeper zones where multi‐decadal surface signals have disturbed subsurface thermal regimes.     

Determination of upward/downward groundwater fluxes using transient variations of soil profile temperature: test of the method with Voyons (Aube,France) experimental data

The temperature variations recorded at several points of a vertical shallow‐depth profile are governed both by conductive and convective heat transfers and can be used to calculate the vertical component of the Darcy velocity and thermal diffusivity in the soil. This paper describes such calculations when transient variations over tens of days are considered and tests them using data collected at Voyons (Aube, France). The temperature was recorded during a year and a half period with a 1 h sampling time step at three different depths: 0·2, 0·4 and 0·75 m. By processing the annual variation of temperature, we obtained a value of the Darcy velocity in good agreement with the value of actual/potential evapotranspiration ratio. By processing transient variations, despite the limitation of the calculations due to the lack of sensitivity of the sensors, results obtained at Voyons were in good correlation with tensiometric data. Copyright © 2005 John Wiley & Sons, Ltd.     

Heat tracing of heterogeneous hyporheic exchange adjacent to in‐stream geomorphic features

There are many field techniques used to quantify rates of hyporheic exchange, which can vary in magnitude and direction spatially over distances of only a few metres, both within and between morphological features. We used in‐stream mini‐piezometers and heat transport modelling of stream and streambed temperatures to quantify the rates and directions of water flux across the streambed interface upstream and downstream of three types of in‐stream geomorphic features: a permanent dam, a beaver dam remnant and a stream meander. We derived hyporheic flux estimates at three different depths at six different sites for a month and then paired those flux rates with measurements of gradient to derive hydraulic conductivity (K) of the streambed sediments. Heat transport modelling provided consistent daily flux estimates that were in agreement directionally with hydraulic gradient measurements and also identified vertical heterogeneities in hydraulic conductivity that led to variable hyporheic exchange. Streambed K varied over an order of magnitude (1·9 × 10^?6 to 5·7 × 10^?5 m/s). Average rates of hyporheic flux ranged from static (q < ±0·02 m/day) to 0·42 m/day. Heat transport modelling results suggest three kinds of flow around the dams and the meander. Copyright © 2010 John Wiley & Sons, Ltd.