Philip Brunner

Hydrogeologic controls on disconnection between surface water and groundwater

[1] Understanding how changes in the groundwater table affect surface water resources is of fundamental importance in quantitative hydrology. If the groundwater table below a stream is sufficiently deep, changes in the groundwater table position effectively do not alter the infiltration rate. This is referred to as a disconnected system. Previous authors noted that a low-conductivity layer below the surface water body is a necessary but not sufficient criterion for disconnection to occur. We develop a precise criterion that allows an assessment of whether surface water-groundwater systems can disconnect or not. We further demonstrate that a disconnected system can be conceptualized by a saturated groundwater mound and the development of a capillary zone above this mound. This conceptualization is used to determine the critical water table position at the point where full disconnection is reached. A comparison of this calculated critical water table position with a measurement of the water table depth in a borehole allows the assessment of the disconnection status. A sensitivity analysis of this critical water table showed that for a given aquifer thickness and river width, the depth to groundwater where the system disconnects is approximately proportional to the stream depth and the hydraulic conductivity of the streambed sediments and inversely proportional to the thickness of these sediments and the hydraulic conductivity of the aquifer. The conceptualization also allows the disconnection problem to be analyzed using both variably saturated and fully saturated groundwater models and provides guidance for numerical and analytical approaches.

Disconnected Surface Water and Groundwater: From Theory to Practice

When describing the hydraulic relationship between rivers and aquifers, the term disconnected is frequently misunderstood or used in an incorrect way. The problem is compounded by the fact that there is no definitive literature on the topic of disconnected surface water and groundwater. We aim at closing this gap and begin the discussion with a short introduction to the historical background of the terminology. Even though a conceptual illustration of a disconnected system was published by Meinzer (1923), it is only within the last few years that the underlying physics of the disconnection process has been described. The importance of disconnected systems, however, is not widely appreciated. Although rarely explicitly stated, many approaches for predicting the impacts of groundwater development on surface water resources assume full connection. Furthermore, management policies often suggest that surface water and groundwater should only be managed jointly if they are connected. However, although lowering the water table beneath a disconnected section of a river will not change the infiltration rate at that point, it can increase the length of stream that is disconnected. Because knowing the state of connection is of fundamental importance for sustainable water management, robust field methods that allow the identification of the state of connection are required. Currently, disconnection is identified by showing that the infiltration rate from a stream to an underlying aquifer is independent of the water table position or by identifying an unsaturated zone under the stream. More field studies are required to develop better methods for the identification of disconnection and to quantify the implications of heterogeneity and clogging processes in the streambed on disconnection.

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.

An Analysis of River Bank Slope and Unsaturated Flow Effects on Bank Storage

Recognizing the underlying mechanisms of bank storage and return flow is important for understanding streamflow hydrographs. Analytical models have been widely used to estimate the impacts of bank storage, but are often based on assumptions of conditions that are rarely found in the field, such as vertical river banks and saturated flow. Numerical simulations of bank storage and return flow in river-aquifer cross sections with vertical and sloping banks were undertaken using a fully-coupled, surface-subsurface flow model. Sloping river banks were found to increase the bank infiltration rates by 98% and storage volume by 40% for a bank slope of 3.4° from horizontal, and for a slope of 8.5°, delay bank return flow by more than four times compared with vertical river banks and saturated flow. The results suggested that conventional analytical approximations cannot adequately be used to quantify bank storage when bank slope is less than 60° from horizontal. Additionally, in the unconfined aquifers modeled, the analytical solutions did not accurately model bank storage and return flow even in rivers with vertical banks due to a violation of the dupuit assumption. Bank storage and return flow were also modeled for more realistic cross sections and river hydrograph from the Fitzroy River, Western Australia, to indicate the importance of accurately modeling sloping river banks at a field scale. Following a single wet season flood event of 12 m, results showed that it may take over 3.5 years for 50% of the bank storage volume to return to the river.

A hydraulic mixing-cell method to quantify the groundwater component of streamflow within spatially distributed fully integrated surface water-groundwater flow models

The complexity of available hydrological models continues to increase, with fully integrated surface water-groundwater flow and transport models now available. Nevertheless, an accurate quantification of streamflow generation mechanisms within these models is not yet possible. For example, such models do not report the groundwater component of streamflow at a particular point along the stream. Instead, the groundwater component of streamflow is approximated either from tracer transport simulations or by the sum of exchange fluxes between the surface and the subsurface along the river. In this study, a hydraulic mixing-cell (HMC) method is developed and tested that allows to accurately determine the groundwater component of streamflow by using only the flow solution from fully integrated surface water-groundwater flow models. By using the HMC method, the groundwater component of streamflow can be extracted accurately at any point along a stream provided the subsurface/surface exchanges along the stream are calculated by the model. A key advantage of the HMC method is that only hydraulic information is used, thus the simulation of tracer transport is not required. Two numerical experiments are presented, the first to test the HMC method and the second to demonstrate that it quantifies the groundwater component of streamflow accurately.

Generating soil electrical conductivity maps at regional level by integrating measurements on the ground and remote sensing data

In arid and semi‐arid areas, salinization of soil and water resources is one of the major threats to irrigated agriculture. For management purposes, quantifying both the extent and distribution of salinization is important, but accurate data with sufficient spatial resolution are often not available. Commonly used techniques such as soil sampling and geophysical methods are time‐consuming and yield only point data. A method is described in which multispectral remote sensing images can be used to regionalize point data measured on the field. Field data consist of measurements of electrical conductivity and are obtained by the combination of geophysical methods and the analysis of field soil samples. Uncalibrated salinity maps were calculated with spectral correlation mapping using image‐based reference spectra of saline areas. As an alternative indicator for soil salinity, the NDVI was used. The method was verified in the Yanqi Basin, northwestern China. Correlations between field data and the uncalibrated salinity maps were found over non‐irrigated sites for all images. Good correlations (R 2 up to 0.85) resulted for images collected during the winter months. The high correlation coefficients allow the uncalibrated salinity maps to be scaled to electrical conductivity maps.

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.

Rapid identification of transience in streambed conductance by inversion of floodwave responses

Streambed conductance controls the interaction between surface and groundwater. However, the streambed conductance is often subject to transience. Directly measuring hydraulic properties in a river yields only point values, is time-consuming and therefore not suited to detect transience of physical properties. Here, we present a method to continuously monitor transience in streambed conductance. Input data are time series of stream stage and near stream hydraulic head. The method is based on the inversion of floodwave responses. The analytical model consists of three parameters: x, the distance between streambank and an observation well, α, the aquifer diffusivity, and a the retardation coefficient that is inversely proportional to the streambed conductance. Estimation of a is carried out over successive time steps in order to identify transience in streambed conductance. The method is tested using synthetic data and is applied to field data from the Rhone River and its alluvial aquifer (Switzerland). The synthetic method demonstrated the robustness of the proposed methodology. Application of the method to the field data allowed identifying transience in streambed properties, following flood events in the Rhone. This method requires transience in the surface water, and the river should not change its width significantly with a rising water level. If these conditions are fulfilled, this method allows for a rapid and effective identification of transience in streambed conductance.

Integrating hydrological modelling, data assimilation and cloud computing for real-time management of water resources

Online data acquisition, data assimilation and integrated hydrological modelling have become more and more important in hydrological science. In this study, we explore cloud computing for integrating field data acquisition and stochastic, physically-based hydrological modelling in a data assimilation and optimisation framework as a service to water resources management. For this purpose, we developed an ensemble Kalman filter-based data assimilation system for the fully-coupled, physically-based hydrological model HydroGeoSphere, which is able to run in a cloud computing environment. A synthetic data assimilation experiment based on the widely used tilted V-catchment problem showed that the computational overhead for the application of the data assimilation platform in a cloud computing environment is minimal, which makes it well-suited for practical water management problems. Advantages of the cloud-based implementation comprise the independence from computational infrastructure and the straightforward integration of cloud-based observation databases with the modelling and data assimilation platform. A cloud-based real-time modelling and data assimilation framework is established.Can be connected to a cloud-based data acquisition and monitoring module.HydroGeoSphere is used as the hydrological forward model.The tiltedV-catchment problem is used as a benchmark for the system.

When can inverted water tables occur beneath streams

Decline in regional water tables (RWT) can cause losing streams to disconnect from underlying aquifers. When this occurs, an inverted water table (IWT) will develop beneath the stream, and an unsaturated zone will be present between the IWT and the RWT. The IWT marks the base of the saturated zone beneath the stream. Although a few prior studies have suggested the likelihood of an IWT without a clogging layer, most of them have assumed that a low-permeability streambed is required to reduce infiltration from surface water to groundwater, and that the IWT only occurs at the bottom of the low-permeability layer. In this study, we use numerical simulations to show that the development of an IWT beneath an unclogged stream is theoretically possible under steady-state conditions. For a stream width of 1 m above a homogeneous and isotropic sand aquifer with a 47 m deep RWT (measured in an observation point 20 m away from the center of the stream), an IWT will occur provided that the stream depth is less than a critical value of 4.1 m. This critical stream depth is the maximum water depth in the stream to maintain the occurrence of an IWT. The critical stream depth decreases with stream width. For a stream width of 6 m, the critical stream depth is only 1 mm. Thus while theoretically possible, an IWT is unlikely to occur at steady state without a clogging layer, unless a stream is very narrow or shallow and the RWT is very deep.