Spatial and temporal intensification of lateral hyporheic flux in narrowing intra‐meander zones

Meander bends in alluvial rivers morphologically evolve towards meander cut‐off with narrowing intra‐meander necks, and this should steepen hydraulic gradients and intensify intra‐meander hyporheic flux. This research used dye tracking and head loss measurements in a 1 : 500 planimetrically scaled laboratory river table to quantify the spatial and temporal intensification of intra‐meander flux rates at two evolution ages. The younger meander bend, M1, had a sinuosity of 2.3, a river neck width of 0.39 cm, and 0.6% river slope, and the older meander bend, M3, had a sinuosity of 5.2, a river neck width of 0.12 cm, and 0.5% river slope. Flux into and out of the meander bend was estimated along the normalized curvilinear distance s*, with the meander neck at s* = 0.1 and s* = 0.9, the meander centroid at s* = 0.37 and s* = 0.63, and the apex at s* = 0.5. Between the meander centroid and neck, we documented a 60% spatial intensification for M1 and a 90% spatial intensification for M3. Between M1 and M3, we documented a 135% temporal intensification at the neck and a 100% intensification at the centroid. Our empirical spatial and temporal intensification rates involving the M1 and the M3 scenario were one to three times lower than theoretical rates derived from a river evolution model with equivalent M1 and M3 planimetry. Overestimation by the theoretical model was attributed to exaggerated head loss caused by the model neglecting groundwater contributions to river stage. Hyporheic exchange provides critical ecosystem services, and its spatial and temporal variation with meander evolution should be considered in river management. Copyright © 2012 John Wiley & Sons, Ltd.     

Flow‐Through Stream Modeling with MODFLOW and MT3D: Certainties and Limitations

This paper aims to assess MODFLOW and MT3D capabilities for simulating the spread of contaminants from a river exhibiting an unusual relationship with an alluvial aquifer, with the groundwater head higher than the river head on one side and lower on the other (flow‐through stream). A series of simulation tests is conducted using a simple hypothetical model so as to characterize and quantify these limitations. Simulation results show that the expected contaminant spread could be achieved with a specific configuration composed of two sets of parameters: (1) modeled object parameters (hydraulic groundwater gradient, hydraulic conductivity values of aquifer and streambed), and (2) modeling parameters (vertical discretization of aquifer, horizontal refinement of stream modeled with River [RIV] package). The influence of these various parameters on simulation results is investigated, and potential complications and errors are identified. Contaminant spread from stream to aquifer is not always reproduced by MT3D due to the RIV package's inability to simulate lateral exchange fluxes between stream and aquifer. This paper identifies the need for a MODFLOW streamflow package allowing lateral stream‐aquifer interactions and streamflow routine calculations. Such developments could be of particular interest for modeling contaminated flow‐through streams.     

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 Grain to Floodplain: Evaluating heterogeneity of floodplain hydrostatigraphy using sedimentology,geophysics, and remote sensing

Stratigraphy is a fundamental component of floodplain heterogeneity and hydraulic conductivity and connectivity of alluvial aquifers, which affect hydrologic processes such as groundwater flow and hyporheic exchange. Watershed-scale hydrological models commonly simplify the sedimentology and stratigraphy of floodplains, neglecting natural floodplain heterogeneity and anisotropy. This study, conducted in the upper reach of the East River in the East River Basin, Colorado, USA, combines point-, meander-, and floodplain-scale data to determine key features of alluvial aquifers important for estimating hydrologic processes. We compare stratigraphy of two meanders with disparate geometries to explore floodplain heterogeneity and connectivity controls on flow and transport. Meander shape, orientation, and internal stratigraphy affected residence time estimates of laterally exchanged hyporheic water. Although the two meanders share a sediment source, vegetation, and climate, their divergent river migration histories resulted in contrasting meander hydrofacies. In turn, the extent and orientation of these elements controlled the effective hydraulic conductivity and, ultimately, estimates of groundwater transport and hyporheic residence times. Additionally, the meanders’ orientation relative to the valley gradient impacted the hydraulic gradient across the meanders—a key control of groundwater velocity. Lastly, we combine our field data with remotely sensed data and introduce a potential approach to estimate key hydrostratigraphic packages across floodplains. Prospective applications include contaminant transport studies, hyporheic models, and watershed models. © 2019 John Wiley & Sons, Ltd.     

Variability in the vertical hyporheic water exchange affected by hydraulic conductivity and river morphology at a natural confluent meander bend

River confluences and their associated tributaries are key morphodynamic nodes that play important roles in controlling hydraulic geometry and hyporheic water exchange in fluvial networks. However, the existing knowledge regarding hyporheic water exchange associated with river confluence morphology is relatively scarce. On January 14 and 15, 2016, the general hydraulic and morphological characteristics of the confluent meander bend (CMB) between the Juehe River and the Haohe River in the southern region of Xi'an City, Shaanxi Province, China, were investigated. The patterns and magnitudes of vertical hyporheic water exchange (VHWE) were estimated based on a one‐dimensional heat steady‐state model, whereas the sediment vertical hydraulic conductivity (K_v) was calculated via in situ permeameter tests. The results demonstrated that 6 hydrodynamic zones and their extensions were observed at the CMB during the test period. These zones were likely controlled by the obtuse junction angle and low momentum flux ratio, influencing the sediment grain size distribution of the CMB. The VHWE patterns at the test site during the test period mostly showed upwelling flow dominated by regional groundwater discharging into the river. The occurrence of longitudinal downwelling and upwelling patterns along the meander bend at the CMB was likely subjected to the comprehensive influences of the local sinuosity of the meander bend and regional groundwater discharge and finally formed regional and local flow paths. Additionally, in dominated upwelling areas, the change in VHWE magnitudes was nearly consistent with that in K_v values, and higher values of both variables generally occurred in erosional zones near the thalweg paths of the CMB, which were mostly made up of sand and gravel. This was potentially caused by the erosional and depositional processes subjected to confluence morphology. Furthermore, lower K_v values observed in downwelling areas at the CMB were attributed to sediment clogging caused by local downwelling flow. The confluence morphology and sediment K_v are thus likely the driving factors that cause local variations in the VHWE of fluvial systems.     

Investigating Groundwater‐Lake Interactions by Hydraulic Heads and a Water Balance

Discharge of groundwater into lakes (lacustrine groundwater discharge, LGD) can play a major role in water balances of lakes. Unfortunately, studies often neglect this input path because of methodological difficulties in its determination. Direct measurements of LGD are labor‐consuming and prone to error. The present study uses both spatially variable hydraulic‐head data and meteorological data to estimate groundwater input by LGD and lake water output through infiltration. The study sites are two shallow, groundwater‐fed lakes without any surface inflows or outflows. Horizontally interpolated groundwater heads were combined with lake water levels to obtain vertical hydraulic gradients between the aquifer and the lake, which are separated by a thick layer of lake bed sediment which has a much lower hydraulic conductivity than the underlying aquifer. By fitting the hydraulic gradient to the results of a simple mass balance and considering the process of clogging, we were able to estimate the hydraulic conductivity of the lake bed sediments. We calculated groundwater inputs by LGD and lake water outputs by infiltration on an annual basis. Although our method requires several assumptions, the results are reasonable and provide useful information about the exchange between the aquifer and the lake, which can, for example, be used for the calculation of nutrient mass balances.     

Geological controls on the formation of alluvial meanders and floodplain wetlands: the example of the Klip River,eastern Free State,South Africa

Floodplain wetlands are common features of rivers in southern Africa, but they have been little studied from a geological or geomorphological perspective. Study of the upper Klip River, eastern Free State, South Africa, indicates strong geological controls on the formation of alluvial meanders and associated floodplain wetlands. Along this river, pronounced and abrupt changes in valley width are strongly linked to lithological variations. Where weakly cemented sandstone crops out, the Klip has laterally eroded bedrock and carved valleys up to 1500 m wide. In these valleys, the river meanders (sinuosity up to ～1·75) on moderate gradients (<0·001) within extensive floodplains marked by numerous oxbow lakes, backswamps and abandoned channels, many of which host substantial wetlands. In contrast, where highly resistant dolerite crops out, lateral erosion of bedrock is restricted, with the Klip tending instead to erode vertically along joints or fractures. Here, valleys are narrower (<200 m), channel‐bed gradients are steeper (>0·003), the river follows a much straighter course (sinuosity ～1·10–1·34), and floodplains are restricted in width. Long‐term landscape development in the Klip and numerous similar catchments depends on the interaction between fluvial processes in the sandstone and dolerite valleys. In the sandstone valleys, vertical erosion rates are controlled by erosion rates of the more resistant dolerites downstream. Hence, in the short‐ to medium‐term (decades to tens of thousands of years), lateral erosion dominates over vertical erosion, with the river concomitantly planing sandstone in the channel floor and reworking floodplain sediments. The thickness of alluvial fill in the sandstone valleys is limited (<4 m), but the resultant meanders are naturally dynamic, with processes such as point bar deposition, cutoff formation and channel avulsion resulting in an assemblage of fluvial landforms. In the longer term (greater than tens of thousands of years), however, vertical erosion will occur in the sandstone valleys as the downstream dolerites are lowered by erosion, resulting in channel incision, floodplain abandonment, and desiccation of the wetlands. Identification of the geological controls on meander and wetland formation provides information vital for the design of effective management guidelines for these ecologically rich habitats, and also contributes to a better understanding of rivers that are intermediate between fully alluvial and fully bedrock. Copyright © 2002 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.     

Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

Implementation of Solute Transport in the Vadose Zone into the “HYDRUS Package for MODFLOW”

The “HYDRUS package for MODFLOW” is an existing MODFLOW package that allows MODFLOW to simultaneously evaluate transient water flow in both unsaturated and saturated zones. The package is based on incorporating parts of the HYDRUS-1D model (to simulate unsaturated water flow in the vadose zone) into MODFLOW (to simulate saturated groundwater flow). The coupled model is effective in addressing spatially variable saturated-unsaturated hydrological processes at the regional scale. However, one of the major limitations of this coupled model is that it does not have the capability to simulate solute transport along with water flow and therefore, the model cannot be employed for evaluating groundwater contamination. In this work, a modified unsaturated flow and transport package (modified HYDRUS package for MODFLOW and MT3DMS) has been developed and linked to the three-dimensional (3D) groundwater flow model MODFLOW and the 3D groundwater solute transport model MT3DMS. The new package can simulate, in addition to water flow in the vadose zone, also solute transport involving many biogeochemical processes and reactions, including first-order degradation, volatilization, linear or nonlinear sorption, one-site kinetic sorption, two-site sorption, and two-kinetic sites sorption. Due to complex interactions at the groundwater table, certain modifications of the pressure head (compared to the original coupling) and solute concentration profiles were incorporated into the modified HYDRUS package. The performance of the newly developed model is evaluated using HYDRUS (2D/3D), and the results indicate that the new model is effective in simulating the movement of water and contaminants in the saturated-unsaturated flow domains.     

Sinuosity evolution along an incising channel: New insights from the Jordan River response to the Dead Sea level fall

The geomorphic evolution of the Jordan River in recent decades indicates that interaction between incision and high-magnitude floods controls sinuosity changes under increasing mouth gradients during base-level fall. The evolution of the river was analyzed based on digital elevation models, remotely sensed imagery, hydrometric data, and a hydraulic model. The response varies along the river. Near the river mouth, where incision rate is high and a deep channel forms, overbank flooding is less likely. There, large floods exert high shear stress within the confined channel, increasing sinuosity. Upstream, near the migrating knickzone channel gradients also increase, incision is more moderate and floods continue to overtop the banks, favoring meander chute cutoffs. The resulting channel has a downstream well-confined meandering segment and an upstream low-sinuosity segment. These new insights regarding spatial differences along an incising channel can improve interpretations of the evolution of ancient planforms and floodplains that responded to base-level decline. © 2018 John Wiley & Sons, Ltd.     

Influence of a thin veneer of low‐hydraulic‐conductivity sediment on modelled exchange between river water and groundwater in response to induced infiltration

A thin layer of fine‐grained sediment commonly is deposited at the sediment–water interface of streams and rivers during low‐flow conditions, and may hinder exchange at the sediment–water interface similar to that observed at many riverbank‐filtration (RBF) sites. Results from a numerical groundwater‐flow model indicate that a low‐permeability veneer reduces the contribution of river water to a pumping well in a riparian aquifer to various degrees, depending on simulated hydraulic gradients, hydrogeological properties, and pumping conditions. Seepage of river water is reduced by 5–10% when a 2‐cm thick, low‐permeability veneer is present on the bed surface. Increasing thickness of the low‐permeability layer to 0·1 m has little effect on distribution of seepage or percentage contribution from the river to the pumping well. A three‐orders‐of‐magnitude reduction in hydraulic conductivity of the veneer is required to reduce seepage from the river to the extent typically associated with clogging at RBF sites. This degree of reduction is much larger than field‐measured values that were on the order of a factor of 20–25. Over 90% of seepage occurs within 12 m of the shoreline closest to the pumping well for most simulations. Virtually no seepage occurs through the thalweg near the shoreline opposite the pumping well, although no low‐permeability sediment was simulated for the thalweg. These results are relevant to natural settings that favour formation of a substantial, low‐permeability sediment veneer, as well as central‐pivot irrigation systems, and municipal water supplies where river seepage is induced via pumping wells. Published in 2011 by John Wiley & Sons, Ltd.     

A conceptual–numerical model to simulate hydraulic head in aquifers that are hydraulically connected to surface water bodies

In this paper, we present a conceptual‐numerical model that can be deduced from a calibrated finite difference groundwater‐flow model, which provides a parsimonious approach to simulate and analyze hydraulic heads and surface water body–aquifer interaction for linear aquifers (linear response of head to stresses). The solution of linear groundwater‐flow problems using eigenvalue techniques can be formulated with a simple explicit state equation whose structure shows that the surface water body–aquifer interaction phenomenon can be approached as the drainage of a number of independent linear reservoirs. The hydraulic head field could be also approached by the summation of the head fields, estimated for those reservoirs, defined over the same domain set by the aquifer limits, where the hydraulic head field in each reservoir is proportional to a specific surface (an eigenfunction of an eigenproblem, or an eigenvector in discrete cases). All the parameters and initial conditions of each linear reservoir can be mathematically defined in a univocal way from the calibrated finite difference model, preserving its characteristics (geometry, boundary conditions, hydrodynamic parameters (heterogeneity), and spatial distribution of the stresses). We also demonstrated that, in practical cases, an accurate solution can be obtained with a reduced number of linear reservoirs. The reduced computational cost of these solutions can help to integrate the groundwater component within conjunctive use management models. Conceptual approximation also facilitates understanding of the physical phenomenon and analysis of the factors that influence it. A simple synthetic aquifer has been employed to show how the conceptual model can be built for different spatial discretizations, the parameters required, and their influence on the simulation of hydraulic head fields and stream–aquifer flow exchange variables. A real‐world case was also solved to test the accuracy of the proposed approaches, by comparing its solution with that obtained using finite‐difference MODFLOW code. Copyright © 2011 John Wiley & Sons, Ltd.     

Modeling Surface Water-Groundwater Interaction with MODFLOW: Some Considerations

The accuracy with which MODFLOW simulates surface water-groundwater interaction is examined for connected and disconnected losing streams. We compare the effect of different vertical and horizontal discretization within MODFLOW and also compare MODFLOW simulations with those produced by HydroGeoSphere. HydroGeoSphere is able to simulate both saturated and unsaturated flow, as well as surface water, groundwater and the full coupling between them in a physical way, and so is used as a reference code to quantify the influence of some of the simplifying assumptions of MODFLOW. In particular, we show that (1) the inability to simulate negative pressures beneath disconnected streams in MODFLOW results in an underestimation of the infiltration flux; (2) a river in MODFLOW is either fully connected or fully disconnected, while in reality transitional stages between the two flow regimes exist; (3) limitations in the horizontal discretization of the river can cause a mismatch between river width and cell width, resulting in an error in the water table position under the river; and (4) because coarse vertical discretization of the aquifer is often used to avoid the drying out of cells, this may result in an error in simulating the height of the groundwater mound. Conditions under which these errors are significant are investigated.     

Improving Groundwater Model Calibration with Repeat Microgravity Measurements

Groundwater-flow models depend on hydraulic head and flux observations for evaluation and calibration. A different type of observation—change in storage measured using repeat microgravity—can also be used for parameter estimation by simulating the expected change in gravity from a groundwater model and including the observation misfit in the objective function. The method is demonstrated using new software linked to MODFLOW input and output files and field data from the vicinity of the All American Canal in southeast California, USA. Over a 10-year period following lining of the previously highly permeable canal with concrete, gravity decreased by over 100 μGal (equivalent to about 2.5 m of free-standing water) at some locations as seepage decreased and the remnant groundwater mound dissipated into the aquifer or was removed by groundwater pumping. Simulated gravity from a MODFLOW model closely matched observations, and repeat microgravity data proved useful for constraining both hydraulic conductivity and specific yield estimates. Specific yield estimated using the infinite-horizontal slab approximation agreed well with model-derived values, and the departure from the linear, flat-water-table approximation was small, less than 2%, despite relatively large and dynamic water-table slope. First-order second-moment parameter uncertainty analysis shows reduction in uncertainty for all hydraulic conductivity and specific yield parameter estimates with the addition of repeat microgravity data, as compared to drawdown data alone.     

Modelling mid‐channel bars in meandering channels

An analytical modelling framework is proposed to reproduce the frequently observed but poorly studied occurrence of mid‐channel bars in meandering channels. Mid‐channel bars occur in meanders and may characterize transitional morphologies between pure meandering and braided rivers. Based on existing field and experimental observations, we propose that two different mechanisms can generate central topographical patterns in meanders. A former mechanism (‘width‐forced’) is related to spatial width oscillations which determine a laterally symmetrical bed shear stress pattern that promotes mid‐channel bars. A second mechanism (‘curvature‐forced’) can take place also in curvilinear equiwidth streams since also longitudinal variations of channel curvature can produce laterally symmetrical alterations of the sediment transport capacity. A perturbation approach is employed to model both mechanisms within a common framework, allowing reproduction, at least qualitatively, of several observed features. While width‐forced mid‐channel bars are a symmetric linear altimetric response, to reproduce curvature‐forced mid‐channel bars requires modelling nonlinear flow‐bed topography interactions at the second order of the perturbation expansion. Hypotheses on how these mechanisms operate are further discussed through an application to field cases. The amplitude of the nonlinear response can be relevant compared to that of the point bar in equiwidth meanders and the location of mid‐channel bars seldom coincides with bend apexes, mainly depending upon the intrinsic meander wavelength. Central bars tend to symmetrically divert the flow against the two banks, a process which is proposed as a possible cause of cross‐sectional overwidening, along with the asymmetry between the rates of bank erosion and of the opposite bank accretion. The outcomes of this first modelling step on the subject allow discussion of the mutual feedback processes that characterize interactions between mid‐channel bars and width variations in river meanders. Copyright © 2010 John Wiley & Sons, Ltd.     

Impact of groundwater flow on meandering; example from the Geul River,The Netherlands

The Geul River, located in the south‐eastern part of The Netherlands, is a meandering river with a planform shape characterized by large loops consisting of multiple bends. We evaluate the effect(s) of groundwater flow on the shapes of meanders as a possible explanation for the multi‐bend loops, using a combined meandering–groundwater computer model. In the model seeping groundwater enhances bank erodibility. Based on the simulation results, we present a conceptual, generalized model for groundwater–meandering interaction, based on wavelength selection and fixation effects. Wavelength selection occurs because of the positive feedback between growing meander bends and groundwater flow patterns and velocities. The promoted wavelengths have the same spatial scale as the groundwater flow system in the aquifer underlying the floodplain. In the case of the Geul River these wavelengths are of the order of 100 m. Since groundwater flow velocities are largest close to the recharging hill‐slopes, the seepage‐enhanced bank erodibilities are at a maximum near the floodplain limits. At these locations the difference in erodibility between banks facing the floodplain and those facing the hill slopes is large, so it is difficult for the river to migrate away from the floodplain limits. This causes long stretches of the river to be aligned along the floodplain limits, which we term a fixation effect. This mechanism best explains the multi‐bend loops of the Geul River. The general interaction between groundwater flow and meandering is site specific since it depends on climatic, fluvial and hydrogeological parameters. The Geul is characterized by a wide floodplain and steep hill‐slopes, and it is underlain by coarse‐grained deposits with good aquifer properties, favoring an important groundwater system. Since this kind of river frequently occurs, our results could apply to many other river systems. Copyright © 2008 John Wiley & Sons, Ltd.     

A New Package in MODFLOW to Simulate Unconfined Groundwater Flow in Sloping Aquifers

The nonhorizontal‐model‐layer (NHML) grid system is more accurate than the horizontal‐model‐layer grid system to describe groundwater flow in an unconfined sloping aquifer on the basis of MODFLOW‐2000. However, the finite‐difference scheme of NHML was based on the Dupuit‐Forchheimer assumption that the streamlines were horizontal, which was acceptable for slope less than 0.10. In this study, we presented a new finite‐difference scheme of NHML based on the Boussinesq assumption and developed a new package SLOPE which was incorporated into MODFLOW‐2000 to become the MODFLOW‐SP model. The accuracy of MODFLOW‐SP was tested against solution of Mac Cormack (1969). The differences between the solutions of MODFLOW‐2000 and MODFLOW‐SP were nearly negligible when the slope was less than 0.27, and they were noticeable during the transient flow stage and vanished in steady state when the slope increased above 0.27. We established a model considering the vertical flow using COMSOL Multiphysics to test the robustness of constrains used in MODFLOW‐SP. The results showed that streamlines quickly became parallel with the aquifer base except in the narrow regions near the boundaries when the initial flow was not parallel to the aquifer base. MODFLOW‐SP can be used to predict the hydraulic head of an unconfined aquifer along the profile perpendicular to the aquifer base when the slope was smaller than 0.50. The errors associated with constrains used in MODFLOW‐SP were small but noticeable when the slope increased to 0.75, and became significant for the slope of 1.0.     

In situ quantification of spatial and temporal variability of hyporheic exchange in static and mobile gravel‐bed rivers

Seepage meters modified for use in flowing water were used to directly measure rates of exchange between surface and subsurface water in a gravel‐ and cobble bed river in western Pennsylvania, USA (Allegheny River, Q_mean = 190 m³/s) and a sand‐ and gravel‐bed river in Colorado, USA (South Platte River, Q_mean = 9·7 m³/s). Study reaches at the Allegheny River were located downstream from a dam. The bed was stable with moss, algae, and river grass present in many locations. Median seepage was + 0·28 m/d and seepage was highly variable among measurement locations. Upward and downward seepage greatly exceeded the median seepage rate, ranging from + 2·26 (upward) to ? 3·76 (downward) m/d. At the South Platte River site, substantial local‐scale bed topography as well as mobile bedforms resulted in spatial and temporal variability in seepage greatly in exceedence of the median groundwater discharge rate of 0·24 m/d. Both upward and downward seepage were recorded along every transect across the river with rates ranging from + 2·37 to ? 3·40 m/d. Despite a stable bed, which commonly facilitates clogging by fine‐grained or organic sediments, seepage rates at the Allegheny River were not reduced relative to those at the South Platte River. Seepage rate and direction depended primarily on measurement position relative to local‐ and meso‐scale bed topography at both rivers. Hydraulic gradients were small at nearly all seepage‐measurement locations and commonly were not a good indicator of seepage rate or direction. Therefore, measuring hydraulic gradient and hydraulic conductivity at in‐stream piezometers may be misleading if used to determine seepage flux across the sediment‐water interface. Such a method assumes that flow between the well screen and sediment‐water interface is vertical, which appears to be a poor assumption in coarse‐grained hyporheic settings. Copyright © 2011 John Wiley & Sons, Ltd.     

Application of an integrated surface water‐groundwater model to multi‐aquifers modeling in Choushui River alluvial fan,Taiwan

This research proposes a combination of SWAT and MODFLOW, MD‐SWAT‐MODFLOW, to address the multi‐aquifers condition in Choushui River alluvial fan, Taiwan. The natural recharge and unidentified pumping/recharge are separately estimated. The model identifies the monthly pumping/recharge rates in multi‐aquifers so that the daily streamflow can be simulated correctly. A multi‐aquifers condition means a subsurface formation composed of at least the unconfined aquifer, the confined aquifer, and an in‐between aquitard. In such a case, the variation of groundwater level is related to pumping/recharge activities in vertically adjacent aquifer and the river‐aquifer interaction. Both factors in turn affect the streamflow performance. Results show that MD‐SWAT‐MODFLOW performs better than SWAT alone in terms of simulated streamflow, especially during low flow period, when pumping/recharge rates are properly estimated. A sensitivity analysis of individual parameter suggests that the vertical leakance may be the most sensitive among all investigated MODFLOW parameters in terms of the estimated pumping/recharge among aquifers, and the Latin‐Hypercube‐One‐factor‐At‐a‐Time sensitivity analysis indicates that the hydraulic conductivity of channel is the most sensitive to the model performance. It also points out the necessity to simultaneously estimate pumping/recharge rates in multi‐aquifers. The estimated net pumping rate can be treated as a lower bound of the actual local pumping rate. As a whole, the model provides the spatio‐temporal groundwater use, which gives the authorities insights to manage groundwater resources. Copyright © 2012 John Wiley & Sons, Ltd.