Study of stream temperature dynamics and corresponding heat fluxes within Miramichi River catchments (New Brunswick,Canada)

Water temperature is a key physical habitat determinant in lotic ecosystems as it influences many physical, chemical, and biological properties of rivers. Hence, a good understanding of the thermal regime of rivers and river heat fluxes is essential for effective management of water and fisheries resources. This study dealt with the modelling of river water temperature using a deterministic model. This model calculated the different heat fluxes at the water surface and from the streambed using different hydrometeorological conditions. The water temperature model was applied on two watercourses of different sizes and thermal characteristics, but within a similar meteorological region, namely, the Little Southwest Miramichi River and Catamaran Brook (New Brunswick, Canada). The model was also applied using microclimate data, i.e. meteorological conditions within the river environment (1–2 m above the water surface), for a better estimation of river heat fluxes. Water temperatures at different depths within the riverbed were also used to estimate the streambed heat fluxes. Results showed that microclimate data were essential to get accurate estimates of the surface heat fluxes. Results also showed that for larger river systems, the surface heat fluxes were generally the dominant component of the heat budget with a correspondingly smaller contribution from the streambed. As watercourses became smaller and groundwater contribution more significant, the streambed contribution became important. For instance, approximately 80% of the heat fluxes occurred at the surface for Catamaran Brook (20% from the streambed) whereas the Little Southwest Miramichi River showed values closer to 90% (10% from the streambed). As was reported in previous studies, the solar radiation input dominated the contribution to the heat gain at 63% for Catamaran Brook and 89% for Little Southwest Miramichi River. Copyright © 2011 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.     

Effect of temporal resolution of water level and temperature inputs on numerical simulation of groundwater–surface water flux exchange in a heavily modified urban river

Groundwater interacts with surface water features nearly in all types of landscapes. Understanding these interactions has practical consequences on the quantity and quality of water in either system, because the depletion or contamination of one of the systems will eventually affect the other one. Many studies have shown that the use of heat as natural tracer in conjunction with water level measurements is an effective method for estimating water flow (fluxes) between groundwater and surface water. A number of studies have explored the effects of spatial and temporal variability of groundwater–surface water flux exchanges using temperature and water level measurements; however, the effect of temporal resolution of water level and temperature data on estimating flux remains unexplored. Therefore, this study investigated the effect of temporal resolution of input data on temporal variation of groundwater–surface water flux exchanges. To this end, we calibrated a variably saturated two‐dimensional groundwater flow and heat transport model (VS2DH) at hourly and daily time scales using temperatures measured at multiple depths below the riverbed of the Zenne River, located at a well‐known Belgian brownfield site. Results of the study showed that the computed water flux through the streambed ranged between ?32 mm/day and +25 mm/day using the hourly model and from ?10 mm/day to ?37 mm/day using the daily model. The hourly model resulted in detecting reversal of flow direction inducing short‐term surface water flow into the streambed. However, such events were not captured if daily temperature and water level measurements were used as input. These findings have important implications for understanding contaminant mass flux and their attenuation in the mixing zone of groundwater and surface water. Copyright © 2012 John Wiley & Sons, Ltd.     

Influence of a rock glacier spring on the stream energy budget and cold‐water refuge in an alpine stream

The thermal regimes of alpine streams remain understudied and have important implications for cold‐water fish habitat, which is expected to decline due to climatic warming. Previous research has focused on the effects of distributed energy fluxes and meltwater from snowpacks and glaciers on the temperature of mountain streams. This study presents the effects of the groundwater spring discharge from an inactive rock glacier containing little ground ice on the temperature of an alpine stream. Rock glaciers are coarse blocky landforms that are ubiquitous in alpine environments and typically exhibit low groundwater discharge temperatures and resilience to climatic warming. Water temperature data indicate that the rock glacier spring cools the stream by an average of 3 °C during July and August and reduces maximum daily temperatures by an average of 5 °C during the peak temperature period of the first two weeks in August, producing a cold‐water refuge downstream of the spring. The distributed stream surface and streambed energy fluxes are calculated for the reach along the toe of the rock glacier, and solar radiation dominates the distributed stream energy budget. The lateral advective heat flux generated by the rock glacier spring is compared to the distributed energy fluxes over the study reach, and the spring advective heat flux is the dominant control on stream temperature at the reach scale. This study highlights the potential for coarse blocky landforms to generate climatically resilient cold‐water refuges in alpine streams.     

Beaver dam analogues drive heterogeneous groundwater–surface water interactions

Beaver dam analogues (BDAs) are a cost-effective stream restoration approach that leverages the recognized environmental benefits of natural beaver dams on channel stability and local hydrology. Although natural beaver dams are known to exert considerable influence on the hydrologic conditions of a stream system by mediating geomorphic processes, nutrient cycling, and groundwater–surface water interactions, the impacts of beaver-derived restoration methods on groundwater–surface water exchange are poorly characterized. To address this deficit, we monitored hyporheic exchange fluxes and streambed porewater biogeochemistry across a sequence of BDAs installed along a central Wyoming stream during the summer of 2019. Streambed fluxes were quantified by heat tracing methods and vertical hydraulic gradients. Biogeochemical activity was evaluated using major ion porewater chemistry and principal component analysis. Vertical fluxes of approximately 1.0 m/day were observed around the BDAs, as was the development of spatially heterogeneous zones of nitrate production, groundwater upwelling, and anaerobic reduction. Strong contrasts in hyporheic zone processes were observed across BDAs of differing sizes. This suggests that structures may function with size-dependent behaviour, only altering groundwater–surface water interactions after a threshold hydraulic step height is exceeded. Patterns of hyporheic exchange and biogeochemical cycling around the studied BDAs resemble those around natural beaver dams, suggesting that BDAs may provide comparable benefits to channel complexity and near-stream function over a 1-year period.     

Influence of water flux and redox conditions on chlorobenzene concentrations in a contaminated streambed

Significant natural attenuation may occur on the passage of groundwater plumes through streambed sediments because of the transition from anaerobic to aerobic conditions and an increased microbial activity. Varying directions and magnitudes of water flow in the streambed may enhance or inhibit the supply of oxygen to the streambed and thus influence the redox zoning. In a field study at a small stream in the industrial area of Bitterfeld‐Wolfen, we observed the variability of hydraulic gradients, streambed temperatures, redox conditions and monochlorobenzene (MCB) concentrations in the streambed over the course of 5 months. During the observation period, the hydrologic conditions changed from losing to gaining. Accordingly, the temperature‐derived water fluxes changed from recharge to discharge. Redox conditions were highly variable between ? 170 and 368 mV in the shallow streambed at a depth of 0·1 m below the streambed surface. Deeper in the streambed, at depths of 0·3 m and 0·5 m, the redox conditions were more stable between ? 198 and ? 81 mV and comparable to those typically found in the aquifer. MCB concentrations in the streambed at 0·3 and 0·5 m depth increased with increasing upward water flux. The MCB concentrations in the shallow streambed at 0·1 m depth appeared to be independent of the hydrologic conditions suggesting that degradation of MCB may have occured. Copyright © 2010 John Wiley & Sons, Ltd.     

Estimating groundwater discharge to surface waters using heat as a tracer in low flux environments: the role of thermal conductivity

Analytical modelling of heat transport was used to address effects of uncertainty in thermal conductivity on groundwater–surface water exchange. In situ thermal conductivities and temperature profiles were measured in a coastal lagoon bed where groundwater is known to discharge. The field site could be divided into three sediment zones where significant spatial changes in thermal conductivity on metre to centimetre scale show that spatial variability connected to the sediment properties must be considered. The application of a literature‐based bulk thermal conductivity of 1.84 Wm^?1 °C^?1, instead of field data that ranged from 0.62 to 2.19 W m^?1 °C^?1, produced a mean overestimation of 2.33 cm d^?1 that, considering the low fluxes of the study area, represents an 89% increase and up to a factor of 3 in the most extreme cases. Incorporating the uncertainty due to sediment heterogeneities leads to an irregular trend of the flux distribution from the shore towards the lagoon. The natural variability of the thermal conductivity associated with changes in the sediment composition resulted in a mean variation of ±0.66 cm d^?1 in fluxes corresponding to a change of ±25.4%. The presence of organic matter in the sediments, a common situation in the near‐shore areas of surface water bodies, is responsible for the decrease of thermal conductivity. The results show that the natural variability of sediment thermal conductivity is a parameter to be considered for low flux environments, and it contributes to a better understanding of groundwater–surface water interactions in natural environments. Copyright © 2015 John Wiley & Sons, Ltd.     

Convergence of Groundwater Discharge through the Hyporheic Zone of Streams

Significant attention has been given to hyporheic water fluxes induced by hydromorphologic processes in streambeds and the effects they have on stream ecology. However, the impact of hyporheic fluxes on regional groundwater flow discharge zones as well as the interaction of these flows are much less investigated. The groundwater-hyporheic interactive flow not only governs solute mass and heat transport in streams but also controls the retention of solute and contamination following the discharge of deep groundwater, such as naturally occurring solutes and leakage from geological waste disposal facilities. Here, we applied a physically based modeling approach combined with extensive hydrologic, geologic and geographical data to investigate the effect of hyporheic flow on groundwater discharge in the Krycklan catchment, located in a boreal landscape in Sweden. Regional groundwater modeling was conducted using COMSOL Multiphysics by considering geologic heterogeneity and infiltration constraint of the groundwater circulation intensity. Moreover, the hyporheic flow was analyzed using an exact spectral solution accounting for the fluctuating streambed topography and superimposed with the regional groundwater flow. By comparing the discharge flow fields with and without consideration of hyporheic flows, we found that the divergence of the discharge was substantially enhanced and the distribution of the travel times of groundwater was significantly shifted toward shorter times due to the presence of hyporheic flow. Particularly important is that the groundwater flow paths contract near the streambed interface due to the hyporheic flow, which leads to a phenomenon that we name “fragmentation” of coherent areas of groundwater upwelling in pinhole-shaped stream tubes.     

Spatial and temporal water column and streambed temperature dynamics within an alpine catchment: implications for benthic communities

Despite the significant influence of temperature upon alpine stream benthic communities, thermal regimes of the water column and hyporheic zone of these mountain streams have received limited attention. This paper reports upon a detailed spatio‐temporal study of water column and streambed temperatures undertaken within the Taillon–Gabiétous catchment, French Pyrénées, that aims: (1) to characterize the nature and dynamics of alpine stream water column and streambed thermal patterns; (2) to investigate stream thermal variability under a range of hydroclimatological conditions; and (3) to consider the implications of (1) and (2) for alpine stream benthic communities. The catchment contains four highly dynamic hydrological sources and pathways: (1) two cirque glaciers (Taillon and Gabiétous); (2) seasonal snowpacks; (3) a karst groundwater system; and (4) hillslope aquifers. Water column temperatures were monitored continuously at four sites located along the Taillon glacial stream and at three groundwater springs (two karstic and one hillslope) over the 2002 summer melt season. An eighth site (Tourettes) was established on a predominantly groundwater‐fed stream with limited meltwater input. Bed temperatures (0·05, 0·20 and 0·40 m depth) and river discharge were measured at three sites: (1) the Taillon stream; (2) the Tourettes stream; and (3) below the confluence of (1) and (2). Air temperatures, incoming short‐wave radiation and precipitation were recorded to characterize atmospheric conditions. Glacial stream water column temperatures increased downstream, although groundwater tributaries punctuated longitudinal patterns. Karstic groundwater streams were cooler and more thermally stable than the glacial stream (except at the glacier snout). Hillslope groundwater stream temperatures were most variable and, on average, the warmest of all sites. Streambed temperatures in the glacial stream were coldest and most variable whilst the warmest and least variable streambed temperatures were recorded in an adjacent groundwater tributary. Temperature variability was strongly related to: (1) dynamic water source and pathway contributions; (2) proximity to source; and (3) prevailing hydroclimatological conditions. The high thermal heterogeneity within this catchment may sustain relatively diverse benthic communities, including some endemic Pyrénéan macroinvertebrate taxa. Copyright © 2004 John Wiley & Sons, Ltd.     

Improved spatial delineation of streambed properties and water fluxes using distributed temperature sensing

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre‐optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5‐day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one‐dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one‐dimensional analytical heat‐tracing models. The results show that the amplitude ratio‐based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat‐tracing studies by providing a technique for visualizing and comparing results from fibre‐optic distributed temperature sensing installations and testing the robustness of analytical heat‐tracing models. Copyright © 2016 John Wiley & Sons, Ltd.     

Reducing monitoring gaps at the aquifer–river interface by modelling groundwater–surface water exchange flow patterns

This study investigates spatial patterns and temporal dynamics of aquifer–river exchange flow at a reach of the River Leith, UK. Observations of sub‐channel vertical hydraulic gradients at the field site indicate the dominance of groundwater up‐welling into the river and the absence of groundwater recharge from surface water. However, observed hydraulic heads do not provide information on potential surface water infiltration into the top 0–15 cm of the streambed as these depths are not covered by the existing experimental infrastructure. In order to evaluate whether surface water infiltration is likely to occur outside the ‘window of detection’, i.e. the shallow streambed, a numerical groundwater model is used to simulate hydrological exchanges between the aquifer and the river. Transient simulations of the successfully validated model (Nash and Sutcliff efficiency of 0·91) suggest that surface water infiltration is marginal and that the possibility of significant volumes of surface water infiltrating into non‐monitored shallow streambed sediments can be excluded for the simulation period. Furthermore, the simulation results show that with increasing head differences between river and aquifer towards the end of the simulation period, the impact of streambed topography and hydraulic conductivity on spatial patterns of exchange flow rates decreases. A set of peak flow scenarios with altered groundwater‐surface water head gradients is simulated in order to quantify the potential for surface water infiltration during characteristic winter flow conditions following the observation period. The results indicate that, particularly at the beginning of peak flow conditions, head gradients are likely to cause substantial increase in surface water infiltration into the streambed. The study highlights the potential for the improvement of process understanding of hyporheic exchange flow patterns at the stream reach scale by simulating aquifer‐river exchange fluxes with a standard numerical groundwater model and a simple but robust model structure and parameterization. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

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.     

Assimilation of historical head data to estimate spatial distributions of stream bed and aquifer hydraulic conductivity fields

Management of water resources in alluvial aquifers relies mainly on understanding interactions between hydraulically connected streams and aquifers. Numerical models that simulate this interaction often are used as decision support tools for water resource management. However, the accuracy of numerical predictions relies heavily on unknown system parameters (e.g., streambed conductivity and aquifer hydraulic conductivity), which are spatially heterogeneous and difficult to measure directly. This paper employs an ensemble smoother to invert groundwater level measurements to jointly estimate spatially varying streambed and alluvial aquifer hydraulic conductivity along a 35.6‐km segment of the South Platte River in Northeastern Colorado. The accuracy of the inversion procedure is evaluated using a synthetic experiment and historical groundwater level measurements, with the latter constituting the novelty of this study in the inversion and validation of high‐resolution fields of streambed and aquifer conductivities. Results show that the estimated streambed conductivity field and aquifer conductivity field produce an acceptable agreement between observed and simulated groundwater levels and stream flow rates. The estimated parameter fields are also used to simulate the spatially varying flow exchange between the alluvial aquifer and the stream, which exhibits high spatial variability along the river reach with a maximum average monthly aquifer gain of about 2.3 m³/day and a maximum average monthly aquifer loss of 2.8 m³/day, per unit area of streambed (m²). These results demonstrate that data assimilation inversion provides a reliable and computationally affordable tool to estimate the spatial variability of streambed and aquifer conductivities at high resolution in real‐world systems.     

Quantitative guidance for efficient vertical flow measurements at the sediment–water interface using temperature–depth profiles

Upward discharge to surface water bodies can be quantified using analytical models based on temperature–depth (T-z) profiles. The use of sediment T-z profiles is attractive as discharge estimates can be obtained using point-in-time data that are collected inexpensively and rapidly. Previous studies have identified that T-z methods can only be applied at times of the year when there is significant difference between the streambed–water interface and deeper sediment temperatures (e.g., winter and summer). However, surface water temperatures also vary diurnally, and the influence of these variations on discharge estimates from T-z methods is poorly understood. For this study, synthetic T-z profiles were generated numerically using measured streambed interface temperature data to assess the influence of diurnal temperature variations on discharge estimation and provide insight into the suitable application of T-z methods. Results show that the time of day of data collection can have a substantial influence on vertical flux estimates using T-z methods. For low groundwater discharge fluxes (e.g., 0.1 m d⁻¹), daily transience in streambed temperatures led to relatively large errors in estimated flow magnitude and direction. For higher discharge fluxes (1.5 m d⁻¹), the influence of transient streambed temperatures on discharge estimates was strongly reduced. Discharge estimates from point-in-time T-z profiles were most accurate when the uppermost point in the T-z profile was near the bed interface daily mean (two time periods daily). Where temperature time series data are available, daily averaged T-z profiles can produce accurate discharge estimates across a wide range of discharge rates. Seasonality in shallow groundwater temperature generally had a negligible influence on vertical flow estimates. These findings can be used to plan field campaigns and provide guidance on the optimal application of T-z methods to quantify vertical groundwater discharge to surface water bodies.     

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.     

Delineating and quantifying ground water discharge zones using streambed temperatures

Streambed temperature mapping, hydraulic testing using minipiezometers, and geochemical analyses of interstitial water of the streambed were used to delineate the pattern of ground water discharge in a sandy streambed and to develop a flux-based conceptual model for ground water/surface water interactions. A new and simple empirical method was used to relate fluxes obtained from minipiezometer data to streambed temperatures. The relationship allowed flux to be calculated at locations where only streambed temperature measurements were made. Slug testing and potentiomanometer measurements at 34 piezometers indicated ground water discharge ranged from 0.03 to 446 L/m2/day (and possibly as high as 7060 L/m2/day) along a 60 m long by 11 to 14 m wide reach of river. Complex but similar plan-view patterns of flux were calculated for both summer and winter using hundreds of streambed temperatures measured on a 1 by 2 m grid. The reach was dominated by ground water discharge and 5% to 7% of the area accounted for approximately 20% to 24% of the total discharge. < 12% of the total area consisted of recharge zones or no-discharge zones. A conceptual model for ground water/surface water interactions consisting of five different behaviors was developed based on the magnitude and direction of flux across the surface of the streambed. The behaviors include short-circuit discharge (e.g., high-flow springs), high discharge (e.g., preferential flowpaths), low to moderate discharge, no discharge (e.g., horizontal hyporheic or ground water flow), and recharge. Geological variations at depth played a key role in determining which type of flow behavior occurred in the streambed.     

Estimating streambed parameters for a disconnected river

Evaluation of stream–aquifer interaction and water balance for a catchment often requires specific information on streambed parameters, such as streambed hydraulic conductivity, seepage flux across the streambed and so on. This paper describes a simple, inexpensive instrument that is used to measure these streambed parameters under the condition of a stream disconnected from groundwater. Our method includes a seepage cylinder for simulation of river water depth. The proposed method was applied to estimate the vertical hydraulic conductivity of a streambed and the changes in vertical seepage rate from stream to groundwater with varied stream water depth in the Manasi River of Xinjiang Uygur Autonomous Region, China. The vertical hydraulic conductivities of the streambed determined from 12 sites along the Manasi River vary from 1.01 to 29.m/day where the stream disconnects from the groundwater. The experimental results suggest that there are two kinds of relations between the vertical seepage rate and the simulated stream water depth. One is a linear relation between the two variables with low Reynolds numbers (less than 10); the other is a nonlinear relation (exponential relation) between the two variables with larger Reynolds numbers (greater than 10). This second relationship is quite different from the traditional model that usually calculates the vertical seepage rate from stream to groundwater under the condition of disconnection using a linear relation (Darcy's Law). Our results suggest that a linear relation can only be used for a limited range of river water depth. This method gives a convenient tool for rapidly estimating the streambed hydraulic conductivity and the changes in the vertical seepage rate across streambed with varied stream water depths for the case of a stream disconnected from groundwater. Copyright © 2013 John Wiley & Sons, Ltd.     

A multi-technique approach to determine temporal and spatial variability of groundwater–stream water exchange

Characterizing the spatio-temporal distribution of groundwater–surface water (GW–SW) exchange fluxes is of paramount importance in understanding catchment behavior. A wide range of field-based techniques are available for such characterization. The objective of this study is to quantify the spatio-temporal distribution of the exchange fluxes along the Çakıt stream (Niğde, Turkey) through coupling a set of geophysical techniques and in-stream measurements in a hierarchical manner. First, geological and water quality information were combined at the catchment scale to determine key areas for reach-scale focus. Second, electromagnetic induction (EMI) surveys were conducted along the reach to pinpoint potential groundwater upwelling locations. EMI anomalies guided our focus to a 665 m-long reach of the stream. Along this selected reach, a fibre-optic distributed temperature sensing (FO-DTS) system was utilized to investigate streambed temperature profiles at fine spatial and temporal scales. Furthermore, vertical hydraulic gradients and exchange fluxes were investigated using nested piezometers and vertical temperature profiles, respectively, at two potential upwelling locations and a potential downwelling location identified by previous surveys. The results of the study reveal heterogeneity of vertical water-flow components with seasonal variability. The EMI survey was successful in identifying a localized groundwater upwelling location. FO-DTS measurements revealed a warm temperature anomaly during cold air temperature and low streamflow conditions at the same upwelling site. Our point-based methods, namely vertical temperature profiles and vertical hydraulic gradient estimates, however, did not always provide consistent results with each other and with EMI and FO-DTS measurements. This study, therefore, highlights the opportunities and challenges in incorporating multi-scale observations in a hierarchical manner in characterization of the GW–SW exchange processes that are known to be highly heterogeneous in time and space. Overall, a combination of different methods helps to overcome the limitations of each single method and increases confidence in the obtained results.     

Hydropeaking induces losses from a river reach: observations at multiple spatial scales

In humid regions, where gaining river conditions generally prevail, daily hydroelectric dam releases alter downstream surface water–groundwater interactions by reversing the head gradient between river and adjacent groundwater. Previously, it has been noted that artificial stage changes due to dam releases enhance hyporheic exchange. Here we investigate the regulated Deerfield River in northwestern Massachusetts at multiple scales to evaluate how changing downstream geologic conditions along the river mediate this artificial hyporheic pumping. Water budget analysis indicates that roughly 10% of bank‐stored water is permanently lost from the 19.5‐km river reach, likely as a result of transpiration by bank vegetation. An adjacent reference stream with similar dimensions and geomorphology, but without hydropeaking, shows predictable gaining conditions. Field observations from streambed piezometers and thermistors show that water losses are not uniform throughout the study reach. Riparian aquifer transmissivity in river sub‐reaches largely determines the magnitude of surface water–groundwater exchange as well as net water loss from the river. These newly documented losses from hydropeaking river systems should inform decisions by river managers and hydroelectric operators of additional tradeoffs of oscillatory dam‐release river management. Copyright © 2015 John Wiley & Sons, Ltd.