Length scales and statistical characteristics of outer bank roughness for large elongate meander bends: The influence of bank material properties,floodplain vegetation and flow inundation

This article explores the length scales and statistical characteristics of form roughness along the outer banks of two elongate bends on a large meandering river through investigation of topographic variability of the bank face. The analysis also examines how roughness varies over the vertical height of the banks and when the banks are exposed subaerially and inundated during flood stage. Detailed data on the topography of the outer banks were obtained subaerially using terrestrial LiDAR during low flow conditions and subaqueously using multibeam echo sounding (MBES) during near‐bankfull conditions. The contributions of various length scales of topographic irregularity to roughness for subaerial conditions were evaluated for different elevation contours on the bank faces using Hilbert–Huang Transform (HHT) spectral analysis. Statistical characteristics for discrete areas on the bank faces were determined by calculating the root‐mean‐square of normal distances from a triangulated irregular network (TIN) surface. Results of the HHT analysis show that the characteristics of roughness along bank faces composed primarily of non‐cohesive sediment, and eroding into cropland, vary with bank elevation and exhibit a dominant range of roughness length scales (~15–50 m). However, bank faces composed predominantly of cohesive material and carved into a forested floodplain have relatively uniform topographic roughness characteristics over the vertical extent of the bank face and do not exhibit a dominant roughness length scale or range of length scales. Additionally, comparison between local surface roughness for subaerial versus subaqueous conditions shows that roughness decreases considerably when the banks are submerged, most likely because of the removal of vegetation and eradication of small‐scale erosional features in non‐cohesive bank materials by flow along the bank face. Thus, roughness appears to be linked to the hydraulic conditions affecting the bank, at least relative to conditions that develop when banks are exposed subaerially. Copyright © 2017 John Wiley & Sons, Ltd.     

Assessment of geomorphological bank evolution of the alluvial threshold rivers based on entropy concept parameters

The complex stream bank profiles in alluvial channels and rivers that are formed after reaching equilibrium has been a popular topic of research for many geomorphologists and river engineers. The entropy theory has recently been successfully applied to this problem. However, the existing methods restrict the further application of the entropy parameter to determine the cross-section slope of the river banks. To solve this limitation, we introduce a novel approach in the extraction of the equation based on the calculation of the entropy parameter (λ) and the transverse slope of the bank profile at threshold channel conditions. The effects of different hydraulic and geometric parameters are evaluated on a variation of the entropy parameter. Sensitivity analysis on the parameters affecting the entropy parameter shows that the most effective parameter on the λ-slope multiplier is the maximum slope of the bank profile and the dimensionless lateral distance of the river banks.     

Seepage erosion in layered stream bank material

Current stream restoration practices often require anthropogenic manipulation of natural field soils to reconstruct stream banks in the absence of stabilizing vegetation. For this study, researchers conducted laboratory experiments on reconstructed, non‐vegetated stream banks with layered soils experiencing seepage. The objective of the study was to determine the effect of seepage, pore water pressure, and bank geometry on erosion and bank stability of layered streambanks. The experimental design consisted of an intermediate‐size soil lysimeter packed with a sandy clay loam top soil and an underlying fine sand layer at three bank slopes (90°, 45° and 26°). Shallow groundwater flow and seepage resulted in bank failure of geometrically stable banks. Pop out failures, liquid deformation, and piping were all observed failure mechanisms in the underlying sand material, dependent on the bank angle. Groundwater seepage processes created small‐scale failures of the underlying sand leading to larger‐scale failures of the overlying sandy clay loam. The underlying sand layer eroded according to the initial bank angle and change in overburden loading. The overlying loam layer failed along linear failure planes. The gradually sloped bank (i.e. 26° slope) failed faster, hypothesized to be due to less confining pressure and greater vertical seepage forces. Researchers analyzed the laboratory experiments using the Bank Stability and Toe Erosion Model, version 4·1. The model calculated an accurate shear surface angle similar to the failure angle observed in the lysimeter tests. The model predicted failure only for the undercut 90° bank slope, and indicated stable conditions for the other geometries. Steeper initial bank slopes and undercut banks decreased the bank factor of safety. The observed failure mechanisms and measured saturation data indicated an interaction between overburden pressure, seepage forces, and bank slope on bank stability. Future bank stability modeling would benefit by incorporating lateral seepage erosion and soil liquefaction prediction calculations. Copyright © 2009 John Wiley & Sons, Ltd.     

Quantifying bank storage of variably saturated aquifers

Numerical simulations were conducted to quantify bank storage in a variably saturated, homogenous, and anisotropic aquifer abutting a stream during rising stream stage. Seepage faces and bank slopes ranging from 1/3 to 100/3 were simulated. The initial conditions were assumed steady-state flow with water draining toward the stream. Then, the stream level rose at a constant rate to the specified elevation of the water table given by the landward boundary condition and stayed there until the system reached a new steady state. This represents a highly simplified version of a real world hydrograph. For the specific examples considered, the following conclusions can be made. The volume of surface water entering the bank increased with the rate of stream level rise, became negligible when the rate of rise was slow, and approached a positive constant when the rate was large. Also, the volume decreased with the dimensionless parameter M (the product of the anisotropy ratio and the square of the domain's aspect ratio). When M was large (>10), bank storage was small because most pore space was initially saturated with ground water due to the presence of a significant seepage face. When M was small, the seepage face became insignificant and capillarity began to play a role. The weaker the capillary effect, the easier for surface water to enter the bank. The effect of the capillary forces on the volume of surface water entering the bank was significant and could not be neglected.     

Stream water infiltration, bank storage, and storage zone changes due to stream-stage fluctuations

During a flood period, stream-stage increases induce infiltration of stream water into an aquifer; subsequent declines in stream stage cause a reverse motion of the infiltrated water. This paper presents the results of the water exchange rate between a stream and aquifer, the storage volume of the infiltrated stream water in the surrounding aquifer (bank storage), and the storage zone. The storage zone is the part of aquifer where groundwater is replaced by stream water during the flood. MODFLOW was used to simulate stream–aquifer interactions and to quantify rates of stream infiltration and return flow. MODPATH was used to trace the pathlines of the infiltrated stream water and to determine the size of the storage zone. Simulations were focused on the analyses of the effects of the stream-stage fluctuation, aquifer properties, the hydraulic conductivity of streambed sediments, regional hydraulic gradients, and recharge and evapotranspiration (ET) rates on stream–aquifer interactions. Generally, for a given stream–aquifer system, larger flow rates result from larger stream-stage fluctuations; larger storage volumes and storage zones are produced by larger and longer-lasting fluctuations. For a given stream-stage hydrograph, a lower-permeable streambed, an aquitard, or an anisotropic aquifer of low vertical hydraulic conductivity can significantly reduce the rate of infiltration and limit the size of the storage zone. The bank storage solely caused by the stage fluctuation differs slightly between gaining and losing streams. Short-term rainfall recharge and ET loss in the shallow groundwater slightly influence on the flow rate, but their effects on bank storage in a larger area for a longer period can be considerable.     

Study on the stability of non-cohesive river bank

Frequent bank collapse in nature highlights the need to study the mechanism of bank stability. This paper presents a theoretical analysis and a flume experimental study on the interaction of hydrodynamic conditions and non-cohesive banks of meandering and straight rivers. No bank collapse occurs if the bank angle is smaller than a critical value. The critical angle is a function of a dimensionless parameter K_UD, which is directly proportional to the square of flow velocity near river bed and inversely proportional to the median diameter of bank material. Furthermore, the critical angle reduces with flow velocity and is higher in meandering rivers than in straight rivers. Formulas for estimating the threshold of bank stability are obtained by curve fitting method with experimental data. The results agree with the data measured in the middle and lower reaches of the Yangtze River.     

A model of equilibrium bed topography for meander bends with erodible banks

Channel curvature produces secondary currents and a transverse sloping channel bed, along which the depth increases towards the outer bank. As a result deep pools tend to form adjacent to the outer bank, promoting bank collapse. The interaction of sediment grains with the primary and secondary flow and the transverse sloping bed also causes meanders to move different grain sizes in different proportions and directions, resulting in a consistent sorting pattern. Several models have been developed to describe this process, but they all have the potential to over‐predict pool depth because they cannot account for the influence of erodible banks. In reality, bank collapse might lead to the development of a wider, shallower cross‐section and any resulting flow depth discrepancy can bias associated predictions of flow, sediment transport, and grain‐size sorting. While bed topography, sediment transport and grain sorting in bends will partly be controlled by the sedimentary characteristics of the bank materials, the magnitude of this effect has not previously been explored. This paper reports the development of a model of flow, sediment transport, grain‐size sorting, and bed topography for river bends with erodible banks. The model is tested via intercomparison of predicted and observed bed topography in one low‐energy (5·3 W m^?2 specific stream power) and one high‐energy (43·4 W m^?2) study reach, namely the River South Esk in Scotland and Goodwin Creek in Mississippi, respectively. Model predictions of bed topography are found to be satisfactory, at least close to the apices of bends. Finally, the model is used in sensitivity analyses that provide insight into the influence of bank erodibility on equilibrium meander morphology and associated patterns of grain‐size sorting. The sensitivity of meander response to bank cohesion is found to increase as a function of the available stream power within the two study bends. Copyright © 2002 John Wiley & Sons, Ltd.     

Stochastic erosion of composite banks in alluvial river bends

The erosion of composite river banks is a complex process involving a number of factors including fluvial erosion, seepage erosion, and cantilever mass failure. To predict the rate of bank erosion with these complexities, a stochastic bank erosion model is suitable to define the probability distribution of the controlling variables. In this study, a bank erosion model in a river bend is developed by coupling several bank erosion processes with an existing hydrodynamic and morphological model. The soil erodibility of cohesive bank layers was measured using a submerged jet test apparatus. Seasonal bank erosion rates for four consecutive years at a bend in the Brahmaputra River, India, were measured by repeated bankline surveys. The ability of the model to predict erosion was evaluated in the river bend that displayed active bank erosion. In this study, different monsoon conditions and the distribution functions of two variables were considered in estimating the stochastic bank erosion rate: the probability of the soil erodibility and stochastic stage hydrographs for the nth return period river stage. Additionally, the influences of the deflection angle of the streamflow, longitudinal slope of river channel, and bed material size on bank erosion rate were also investigated. The obtained stochastic erosion predictions were compared with the observed distribution of the annual‐average bank erosion rate of 45 river bends in the Brahmaputra River. The developed model appropriately predicted the short‐term morphological dynamics of sand‐bed river bends with composite banks. Copyright © 2014 John Wiley & Sons, Ltd.     

Channel Representation in Physically Based Models Coupling Groundwater and Surface Water: Pitfalls and How to Avoid Them

Recent models that couple three‐dimensional subsurface flow with two‐dimensional overland flow are valuable tools for quantifying complex groundwater/stream interactions and for evaluating their influence on watershed processes. For the modeler who is used to defining streams as a boundary condition, the representation of channels in integrated models raises a number of conceptual and technical issues. These models are far more sensitive to channel topography than conventional groundwater models. On all spatial scales, both the topography of a channel and its connection with the floodplain are important. For example, the geometry of river banks influences bank storage and overbank flooding; the slope of the river is a primary control on the behavior of a catchment; and at the finer scale bedform characteristics affect hyporheic exchange. Accurate data on streambed topography, however, are seldom available, and the spatial resolution of digital elevation models is typically too coarse in river environments, resulting in unrealistic or undulating streambeds. Modelers therefore perform some kind of manual yet often cumbersome correction to the available topography. In this context, the paper identifies some common pitfalls, and provides guidance to overcome these. Both aspects of topographic representation and mesh discretization are addressed. Additionally, two tutorials are provided to illustrate: (1) the interpolation of channel cross‐sectional data and (2) the refinement of a mesh along a stream in areas of high topographic variability.     

Flow separation at the inner (convex) and outer (concave) banks of constant‐width and widening open‐channel bends

There is a paucity of data and insight in the mechanisms of, and controls on flow separation and recirculation at natural sharply‐curved river bends. Herein we report on successful laboratory experiments that elucidate flow structure in one constant‐width bend and a second bend with an outer‐bank widening. The experiments were performed with both a flat immobile gravel bed and mobile sand bed with dominant bedload sediment transport. In the constant‐width bend with immobile bed, a zone of mainly horizontal flow separation (vertical rotational axis) formed at the inner bank that did not contain detectable flow recirculation, and an outer‐bank cell of secondary flow with streamwise oriented rotational axis. Surprisingly, the bend with widening at the outer bank and immobile bed did not lead to a transverse expansion of the flow. Rather, flow in the outer‐bank widening weakly recirculated around a vertical axis and hardly interacted with the inner part of the bend, which behaved as a constant‐width bend. In the mobile bed experiment, downstream of the bend apex a pronounced depositional bar developed at the inside of the bend and pronounced scour occurred at the outside. Moreover the deformed bed promoted flow separation over the bar, including return currents. In the constant‐width bend, the topographic steering impeded the generation of an outer‐bank cell of secondary flow. In the bend with outer‐bank widening, the topographic steering induced an outward expansion of the flow, whereby the major part of the discharge was conveyed in the central part of the widening section. Flow in the outer‐bank widening was highly three dimensional and included return currents near the bottom. In conclusion, the experiments elucidated three distinct processes of flow separation common in sharp bends: flow separation at the inner bank, an outer‐bank cell of secondary flow, and flow separation and recirculation in an outer‐bank widening. Copyright © 2012 John Wiley & Sons, Ltd.     

Lateral circulation driven by boundary mixing and the associated transport of sediments in idealized partially mixed estuaries

A three-dimensional, hydrostatic, primitive equation numerical model with modern turbulence closures is used to explore lateral circulation and the associated transport of sediments in idealized, moderately to highly stratified estuaries. The model results suggest that boundary mixing on a sloping bottom can drive a significant amount of lateral circulation. This mechanism has received little attention to date in the estuarine literature. Good agreement with an analytical solution and similar vertical structures of lateral flows to observations from the Hudson River estuary support the importance of the boundary mixing mechanism. Boundary mixing is at least as important as differential advection for the modeled scenarios, when the two mechanisms are evaluated using the salt balance equation for model runs without rotation. Linearly superposing analytical solutions for lagged boundary mixing lateral flow and Ekman-forced lateral flow yields a good representation of the near-bottom lateral flow from the model with rotation. The 2 h lag required for the boundary mixing solution is roughly equal to the vertical diffusion time scale, indicating that lateral flow adjustment depends on development of a bottom mixed layer. Sediment dynamics at cross sections seaward and landward of the salt intrusion are very different. Seaward of the salt intrusion, sediments are eroded in the channel and preferentially deposited on the right slope (looking seaward), mainly due to the combination of high sediment concentration in the channel during flood with strong up-slope transport on that side (tidal pumping). Lateral sediment re-distribution landward of the salt intrusion is negligible due to weak residual lateral circulation.     

Intra‐event scale bar–bank interactions and their role in channel widening

Channel bars and banks strongly affect the morphology of both braided and meandering rivers. Accordingly, bar formation and bank erosion processes have been greatly explored. There is, however, a lack of investigations addressing the interactions between bed and bank morphodynamics, especially over short timescales. One major implication of this gap is that the processes leading to the repeated accretion of mid‐channel bars and associated widenings remain unsolved. In a restored section of the Drau River, a gravel‐bed river in Austria, mid‐channel bars have developed in a widening channel. During mean flow conditions, the bars divert the flow towards the banks. One channel section exhibited both an actively retreating bank and an expanding mid‐channel bar, and was selected to investigate the morphodynamic processes involved in bar accretion and channel widening at the intra‐event timescale. We repeatedly surveyed riverbed and riverbank topography, monitored riverbank hydrology and mounted a time‐lapse camera for continuous observation of riverbank erosion processes during four flow events. The mid‐channel bar was shown to accrete when it was submerged during flood events, which at the subsequent flow diversion during lower discharges narrowed the branch along the bank and increased the water surface elevation upstream from the riffle, which constituted the inlet into the branch. These changes of bed topography accelerated the flow along the bank and triggered bank failures up to 20 days after the flood events. Four analysed flow events exhibited a total bar expansion from initially 126 m² to 295 m², while bank retreat was 6 m at the apex of the branch. The results revealed the forcing role of bar accretion in channel widening and highlighted the importance of intra‐event scale bed morphodynamics for bank erosion, which were summarized in a conceptual model of the observed bar–bank interactions. Copyright © 2015 John Wiley & Sons, Ltd.     

Bank erosion and channel width change in a tropical catchment

Catchment sediment budget models are used to predict the location and rates of bank erosion in tropical catchments draining to the Great Barrier Reef lagoon, yet the reliability of these predictions has not been tested due to a lack of measured bank erosion data. This paper presents the results of a 3 year field study examining bank erosion and channel change on the Daintree River, Australia. Three different methods were employed: (1) erosion pins were used to assess the influence of riparian vegetation on bank erosion, (2) bench‐marked cross‐sections were used to evaluate annual changes in channel width and (3) historical aerial photos were used to place the short term data into a longer temporal perspective of channel change (1972–2000). The erosion pin data suggest that the mean erosion rate of banks with riparian vegetation is 6·5 times (or 85%) lower than that of banks without riparian vegetation. The changes measured from cross‐section surveys suggest that channel width has increased by an average of 0·74 (±0·47) m a^?1 over the study period (or ～0·8% yr^?1). The aerial photo results suggest that over the last 30 years the Daintree River has undergone channel contraction of the order of 0·25 m a^?1. The cross‐section data were compared against modelled SedNet bank erosion rates, and it was found that the model underestimated bank erosion and was unable to represent the variable erosion and accretion processes that were observed in the field data. The reach averaged bank erosion rates were improved by the inclusion of locally derived bed slope and discharge estimates; however, the results suggest that it will be difficult for catchment scale sediment budget models to ever accurately predict the location and rate of bank erosion due to the variation in bank erosion rates in both space and time. Copyright © 2008 John Wiley & Sons, Ltd.     

Influence of Groundwater Hydraulic Gradient on Bank Storage Metrics

The hydraulic gradient between aquifers and rivers is one of the most variable properties in a river/aquifer system. Detailed process understanding of bank storage under hydraulic gradients is obtained from a two‐dimensional numerical model of a variably saturated aquifer slice perpendicular to a river. Exchange between the river and the aquifer occurs first at the interface with the unsaturated zone. The proportion of total water exchanged through the river bank compared to the river bed is a function of aquifer hydraulic conductivity, partial penetration, and hydraulic gradient. Total exchange may be estimated to within 50% using existing analytical solutions provided that unsaturated zone processes do not strongly influence exchange. Model‐calculated bank storage is at a maximum when no hydraulic gradient is present and increases as the hydraulic conductivity increases. However, in the presence of a hydraulic gradient, the largest exchange flux or distance of penetration does not necessarily correspond to the highest hydraulic conductivity, as high hydraulic conductivity increases the components of exchange both into and out of an aquifer. Flood wave characteristics do not influence ambient groundwater discharge, and so in large floods, hydraulic gradients must be high to reduce the volume of bank storage. Practical measurement of bank storage metrics is problematic due to the limitations of available measurement technologies and the nested processes of exchange that occur at the river‐aquifer interface. Proxies, such as time series concentration data in rivers and groundwater, require further development to be representative and quantitative.     

Propagation of seasonal temperature signals into an aquifer upon bank infiltration

Infiltrating river water carries the temperature signal of the river into the adjacent aquifer. While the diurnal temperature fluctuations are strongly dampened, the seasonal fluctuations are much less attenuated and can be followed into the aquifer over longer distances. In one-dimensional model with uniform properties, this signal is propagated with a retarded velocity, and its amplitude decreases exponentially with distance. Therefore, time shifts in seasonal temperature signals between rivers and groundwater observation points may be used to estimate infiltration rates and near-river groundwater velocities. As demonstrated in this study, however, the interpretation is nonunique under realistic conditions. We analyze a synthetic test case of a two-dimensional cross section perpendicular to a losing stream, accounting for multi-dimensional flow due to a partially penetrating channel, convective-conductive heat transport within the aquifer, and heat exchange with the underlying aquitard and the land surface. We compare different conceptual simplifications of the domain in order to elaborate on the importance of different system elements. We find that temperature propagation within the shallow aquifer can be highly influenced by conduction through the unsaturated zone and into the underlying aquitard. In contrast, regional groundwater recharge has no major effect on the simulated results. In our setup, multi-dimensionality of the flow field is important only close to the river. We conclude that over-simplistic analytical models can introduce substantial errors if vertical heat exchange at the aquifer boundaries is not accounted for. This has to be considered when using seasonal temperature fluctuations as a natural tracer for bank infiltration.     

Influence of shallowness,bank inclination and bank roughness on the variability of flow patterns and boundary shear stress due to secondary currents in straight open-channels

Boundary shear stress and flow variability due to its interaction with main flow and secondary currents were investigated under conditions that extend previous research on trapezoidal channels. Secondary currents that scale with the flow depth were found over the entire width in all experiments. These findings contradict the widespread perception that secondary currents die out at a distance of 2.5 times the flow depth from the bank, a perception which is largely based on experiments with smooth boundaries. The reported results indicate that a stable pattern of secondary currents over the entire channel width can only be sustained over a fixed horizontal bed if the bed's roughness is sufficient to provide the required transverse oscillations in the turbulent shear stresses. Contrary to laboratory flumes, alluvial river bed always provide sufficient roughness. The required external forcing of this hydrodynamic instability mechanism is provided by the turbulence-generated near-bank secondary currents. The pattern of near-bank secondary currents depends on the inclination and the roughness of the bank. In all configurations, secondary currents result in a reduction of the bed shear stress in the vicinity of the bank and a heterogeneous bank shear stress that reaches a maximum close to the toe of the bank. Moreover, these currents cause transverse variability of 10–15% for the streamwise velocities and 0.2u_*²–0.3u_*?² for the bed shear stress. These variations are insufficient to provide the flow variability required in river restoration projects, but nevertheless must be accounted for in the design of stable channels.     

The effects of ice cover on flow characteristics in a subarctic meandering river

The effects of ice cover on flow characteristics in meandering rivers are still not completely understood. Here, we quantify the effects of ice cover on flow velocity, the vertical and spatial flow distribution, and helical flow structure. Comparison with open‐channel low flow conditions is performed. An acoustic doppler current profiler (ADCP) is used to measure flow from up to three meander bends, depending on the year, in a small sandy meandering subarctic river (Pulmanki River) during two consecutive ice‐covered winters (2014 and 2015). Under ice, flow velocities and discharges were predominantly slower than during the preceding autumn open‐channel conditions. Velocity distribution was almost opposite to theoretical expectations. Under ice, velocities reduced when entering deeper water downstream of the apex in each meander bend. When entering the next bend, velocities increased again together with the shallower depths. The surface velocities were predominantly greater than bottom/riverbed velocities during open‐channel flow. The situation was the opposite in ice‐covered conditions, and the maximum velocities occurred in the middle layers of the water columns. High‐velocity core (HVC) locations varied under ice between consecutive cross‐sections. Whereas in ice‐free conditions the HVC was located next to the inner bank at the upstream cross‐sections, the HVC moved towards the outer bank around the apex and again followed the thalweg in the downstream cross‐sections. Two stacked counter‐rotating helical flow cells occurred under ice around the apex of symmetric and asymmetric bends: next to the outer bank, top‐ and bottom‐layer flows were towards the opposite direction to the middle layer flow. In the following winter, no clear counter‐rotating helical flow cells occurred due to the shallower depths and frictional disturbance by the ice cover. Most probably the flow depth was a limiting factor for the ice‐covered helical flow circulation, similarly, the shallow depths hinder secondary flow in open‐channel conditions. Copyright © 2016 John Wiley & Sons, Ltd.     

Impacts of river bed gas on the hydraulic and thermal dynamics of the hyporheic zone

Despite the presence of gas in river beds being a well known phenomenon, its potential feedbacks on the hydraulic and thermal dynamics of the hyporheic zone has not been widely studied. This paper explores hypotheses that the presence of accumulated gas impacts the hydraulic and thermal dynamics of a river bed due to changes in specific storage, hydraulic conductivity, effective porosity, and thermal diffusivity. The hypotheses are tested using data analysis and modelling for a study site on the urban River Tame, Birmingham, UK. Gas, predominantly attributed to microbial denitrification, was observed in the river bed up to around 14% by volume, and to at least 0.8 m depth below river bed. Numerical modelling indicates that, by altering the relative hydraulic conductivity distribution, the gas in the river bed leads to an increase of groundwater discharge from the river banks (relative to river bed) by a factor of approximately 2 during river low flow periods. The increased compressible storage of the gas phase in the river bed leads to an increase in the simulated volume of river water invading the river bed within the centre of the channel during storm events. The exchange volume can be more than 30% greater in comparison to that for water saturated conditions. Furthermore, the presence of gas also reduces the water-filled porosity, and so the possible depth of such invading flows may also increase markedly, by more than a factor of 2 in the observed case. Observed diurnal temperature variations within the gaseous river bed at 0.1 and 0.5 m depth are, respectively, around 1.5 and 6 times larger than those predicted for saturated sediments. Annual temperature fluctuations are seen to be enhanced by around 4 to 20% compared to literature values for saturated sediments. The presence of gas may thus alter the bulk thermal properties to such a degree that the use of heat tracer techniques becomes subject to a much greater degree of uncertainty. Although the likely magnitude of thermal and hydraulic changes due to the presence of gas for this site have been demonstrated, further research is needed into the origins of the gas and its spatial and temporal variability to enable quantification of the significance of these changes for chemical attenuation and hyporheic zone biology.     

A coupled vertically integrated model to describe lateral exchanges between surface and subsurface in large alluvial floodplains with a fully penetrating river

This paper presents a vertically averaged model for studying water and solute exchanges between a large river and its adjacent alluvial aquifer. The hydraulic model couples horizontal 2D Saint Venant equations for river flow and a 2D Dupuit equation for aquifer flow. The dynamic coupling between river and aquifer is provided by continuity of fluxes and water level elevation between the two domains. Equations are solved simultaneously by linking the two hydrological system matrices in a single global matrix in order to ensure the continuity conditions between river and aquifer and to accurately model two‐way coupling between these two domains. The model is applied to a large reach (about 36 km²) of the Garonne River (south‐western France) and its floodplain, including an instrumented site in a meander. Simulated hydraulic heads are compared with experimental measurements on the Garonne River and aquifer in the floodplain. Model verification includes comparisons for one point sampling date (27 piezometers, 30 March 2000) and for hydraulic heads variations measured continuously over 5 months (5 piezometers, 1 January to 1 June 2000). The model accurately reproduces the strong hydraulic connections between the Garonne River and its aquifer, which are confirmed by the simultaneous variation of the water level in the river and in piezometers located near the river bank. The simulations also confirmed that the model is able to reproduce groundwater flow dynamics during flood events. Given these results, the hydraulic model was coupled with a solute‐transport component, based on advection‐dispersion equations, to investigate the theoretical dynamics of a conservative tracer over 5 years throughout the 36 km² reach studied. Meanders were shown to favour exchanges between river and aquifer, and although the tracer was diluted in the river, the contamination moved downstream from the injection plots and affected both river banks. Copyright © 2008 John Wiley & Sons, Ltd.     

Aquifer response to stream-stage and recharge variations. II. Convolution method and applications

In this second of two papers, analytical step-response functions, developed in the companion paper for several cases of transient hydraulic interaction between a fully penetrating stream and a confined, leaky, or water-table aquifer, are used in the convolution integral to calculate aquifer heads, streambank seepage rates, and bank storage that occur in response to stream-stage fluctuations and basinwide recharge or evapotranspiration. Two computer programs developed on the basis of these step-response functions and the convolution integral are applied to the analysis of hydraulic interaction of two alluvial stream–aquifer systems in the northeastern and central United States. These applications demonstrate the utility of the analytical functions and computer programs for estimating aquifer and streambank hydraulic properties, recharge rates, streambank seepage rates, and bank storage. Analysis of the water-table aquifer adjacent to the Blackstone River in Massachusetts suggests that the very shallow depth of water table and associated thin unsaturated zone at the site cause the aquifer to behave like a confined aquifer (negligible specific yield). This finding is consistent with previous studies that have shown that the effective specific yield of an unconfined aquifer approaches zero when the capillary fringe, where sediment pores are saturated by tension, extends to land surface. Under this condition, the aquifer's response is determined by elastic storage only. Estimates of horizontal and vertical hydraulic conductivity, specific yield, specific storage, and recharge for a water-table aquifer adjacent to the Cedar River in eastern Iowa, determined by the use of analytical methods, are in close agreement with those estimated by use of a more complex, multilayer numerical model of the aquifer. Streambank leakance of the semipervious streambank materials also was estimated for the site. The streambank-leakance parameter may be considered to be a general (or lumped) parameter that accounts not only for the resistance of flow at the river–aquifer boundary, but also for the effects of partial penetration of the river and other near-stream flow phenomena not included in the theoretical development of the step-response functions.