Characterising groundwater–surface water interactions in idealised ephemeral stream systems

Transmission losses from the beds of ephemeral streams are thought to be a widespread mechanism of groundwater recharge in arid and semi-arid regions and support a range of dryland hydro-ecology. Dryland areas cover ~40% of the Earth's land surface and groundwater resources are often the main source of freshwater. It is commonly assumed that where an unsaturated zone exists beneath a stream, the interaction between surface water and groundwater is unidirectional and that groundwater does not exert a significant feedback on transmission losses. To test this assumption, we conducted a series of numerical model experiments using idealised two-dimensional channel-transects to assess the sensitivity and degree of interaction between surface and groundwater for typical dryland ephemeral stream geometries, hydraulic properties and flow regimes. We broaden the use of the term ‘stream–aquifer interactions’ to refer not just to fluxes and water exchange but also to include the ways in which the stream and aquifer have a hydraulic effect on one another. Our results indicate that deep water tables, less frequent streamflow events and/or highly permeable sediments tend to result in limited bi-directional hydraulic interaction between the stream and the underlying groundwater which, in turn, results in high amounts of infiltration. With shallower initial depth to the water table, higher streamflow frequency and/or lower bed permeability, greater ‘negative’ hydraulic feedback from the groundwater occurs which in turn results in lower amounts of infiltration. Streambed losses eventually reach a constant rate as initial water table depths increase, but only at depths of 10s of metres in some of the cases studied. Our results highlight that bi-directional stream–aquifer hydraulic interactions in ephemeral streams may be more widespread than is commonly assumed. We conclude that groundwater and surface water should be considered as connected systems for water resource management unless there is clear evidence to the contrary.     

Heterogeneity influences on stream water–groundwater interactions in a gravel-dominated floodplain

ABSTRACT

Floodplains are composed of complex depositional patterns of ancient and recent stream sediments, and research is needed to address the manner in which coarse floodplain materials affect stream–groundwater exchange patterns. Efforts to understand the heterogeneity of aquifers have utilized numerous techniques typically focused on point-scale measurements; however, in highly heterogeneous settings, the ability to model heterogeneity is dependent on the data density and spatial distribution. The objective of this research was to investigate the correlation between broad-scale methodologies for detecting heterogeneity and the observed spatial variability in stream/groundwater interactions of gravel-dominated alluvial floodplains. More specifically, this study examined the correlation between electrical resistivity (ER) and alluvial groundwater patterns during a flood event at a site on Barren Fork Creek, in the Ozark ecoregion of Oklahoma, USA, where chert gravels were common both as streambed and as floodplain material. Water table elevations from groundwater monitoring wells for a flood event on 1–5 May 2009 were compared to ER maps at various elevations. Areas with high ER matched areas with lower water table slope at the same elevation. This research demonstrated that ER approaches were capable of indicating heterogeneity in surface water–groundwater interactions, and that these heterogeneities were present even in an aquifer matrix characterized as highly conductive. Portions of gravel-dominated floodplain vadose zones characterized by high hydraulic conductivity features can result in heterogeneous flow patterns when the vadose zone of alluvial floodplains activates during storm events.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR X. Chen     

A study of interaction between surface water and groundwater using environmental isotope in Huaisha River basin

The surface water and groundwater are important components of water cycle, and the interaction between surface water and groundwater is the important part in water cycle research. As the effective tracers in water cycle research, environmental isotope and hydrochemistry can reveal the interrelationships between surface water and groundwater effectively. The study area is the Huaisha River basin, which is located in Huairou district, Beijing. The field surveying and sampling for spring, river and well water were finished in 2002 and 2003. The hydrogen and oxygen isotopes and water quality were measured at the laboratory. The spatial characteristics in isotope and evolution of water quality along river lines at the different area were analyzed. The altitude effect of oxygen isotope in springs was revealed, and then using this equation, theory foundation for deducing recharge source of spring was estimated. By applying the mass balance method, the annual mean groundwater recharge rate at the catchment was estimated. Based on the groundwater recharge analysis, combining the hydrogeological condition analysis, and comparing the rainfall-runoff coefficients from the 1960s to 1990s in the Huaisha River basin and those in the Chaobai River basin, part of the runoff in the Huaisha River basin is recharged outside of this basin, in other words, this basin is an un-enclosed basin. On the basis of synthetically analyses, combining the compositions of hydrogen and oxygen isotopes and hydrochemistry, geomorphology, geology, and watershed systems characteristics, the relative contributions between surface water and groundwater flow at the different areas at the catchments were evaluated, and the interaction between surface water and groundwater was re- vealed lastly.     

Chemical and isotopic assessment of surface water–shallow groundwater interaction in the arid Grande river basin,North-Central Chile

ABSTRACT

This paper analyses the composition of surface water and shallow groundwater in the Grande River basin, North-Central Chile, using this information to characterize water interactions. Chemical and isotopic data for surface water and groundwater (7 and 6 sampling locations, respectively) were obtained from three sampling campaigns performed in March–April (autumn), August–September (late winter) and December (early summer) 2012. Precipitation samples were also collected. Data was processed using spatial distribution charts, Piper and Stiff diagrams, and multivariate analysis. In general, the results for each method converge on a high degree of connectivity between surface water and shallow groundwater in the study area. Furthermore, approximately a 10% of groundwater contribution to the surface flow discharge was estimated for a particular reach. This multi-method approach was useful for the characterization of surface water–groundwater interactions in the Grande River basin, and may become a suitable and replicable scheme for studies in arid and semi-arid basins facing similar water management challenges.

Editor D. Koutsoyiannis; Associate editor B. Dewals     

Integrated surface water–groundwater modelling in the context of increasing water reserves of a regional Sahelian aquifer

Abstract

Despite the Sahelian drought of the 1970s–1990s, the unconfined aquifer in southwest Niger exhibits a multidecadal increase in groundwater reserves. Recent changes in land surface conditions have enhanced runoff and thus indirect groundwater recharge below endorheic ponds. This paper presents a model-based investigation of surface runoff and groundwater recharge at mesoscale (～5000 km²). A new lumped-conceptual runoff model applicable to the large number of ungauged endorheic catchments is specially developed, derived from an existing fine-scale, physically-based hydrologic model. Runoff simulated for sites identified as groundwater recharge sources are used to derive recharge forcing for a Modflow-based model of the aquifer. The rising water table trend and its spatial distribution over the period 1992–2003 are generally well simulated, albeit smoothed year-to-year dynamics. Comparison with alternative methods of recharge estimation suggests, however, that there may presently exist more recharging sites and/or contributing surfaces than those considered so far.

Citation Massuel, S., Cappelaere, B., Favreau, G., Leduc, C., Lebel, T. & Vischel, T. (2011) Integrated surface water–groundwater modelling in the context of increasing water reserves of a regional Sahelian aquifer. Hydrol. Sci. J. 56(7), 1242–1264.     

Simulation and control of the linked systems of water quantity–water quality–socio-economics in the Huaihe River basin

Abstract

This work makes explicit an algebraic expression giving the matrix of transient influence coefficients associated with a one-dimensional semi-confined aquifer model. The domain studied is divided into a series of connected and completely mixed compartments over which the governing equation is discretized. The discrete equations obtained are solved for the compartmental hydraulic head and used to derive the algebraic expression in question. The basic properties of the so-called algebraic influence coefficients are investigated. In particular, their consistency with the exact Green function is highlighted. Finally, the newly derived influence coefficients are applied to a simplified aquifer system in order to formulate and solve the problem of identifying illegal groundwater pumping.     

Groundwater–surface water interactions in a river estuary and the importance of geomorphology: Insights from hydraulic,thermal and geophysical observations

This study presents the groundwater flow and salinity dynamics along a river estuary, the Werribee River in Victoria, Australia, at local and regional scales. Along a single reach, salinity across a transverse section of the channel (~80 m long) with a point bar was monitored using time-lapse electrical resistivity (ER) through a tidal cycle. Groundwater fluxes were concurrently estimated by monitoring groundwater levels and temperature profiles. Regional porewater salinity distribution was mapped using 6-km long longitudinal ER surveys during summer and winter. The time-lapse ER across the channel revealed a static electrically resistive zone on the side of the channel with a pronounced cut bank. Upward groundwater flux and steep vertical temperature gradients with colder temperatures deeper within the sediment suggested a stable zone of fresh groundwater discharge along this cut bank area. Generally, less resistive zones were observed at the shallow portion of the inner meander bank and at the channel center. Subsurface temperatures close to surface water values, vertical head gradients indicating both upward and downward groundwater flux, and higher porewater salinity closer to that of estuary water suggest strong hyporheic circulation in these zones. The longitudinal surveys revealed higher ER values along deep and sinuous segments and low ER values in shallow and straighter reaches in both summer and winter; these patterns are consistent with the local channel-scale observations. This study highlights the interacting effects of channel morphology, broad groundwater–surface water interaction, and hyporheic exchange on porewater salinity dynamics underneath and adjacent to a river estuary.     

Spatial–temporal variability in water quality and macro-invertebrate assemblages in the Upper Mara River basin,Kenya

Tropical rivers display profound temporal and spatial heterogeneity in terms of environmental conditions. This aspect needs to be considered when designing a monitoring program for water quality in rivers. Therefore, the physico-chemical composition and the nutrient loading of the Upper Mara River and its two main tributaries, the Amala and Nyangores were monitored. Initial daily, and later a weekly monitoring schedule for 4 months spanning through the wet and dry seasons was adopted. Benthic macro-invertebrates were also collected during the initial sampling to be used as indicators of water quality. The aim of the current study was to investigate the physico-chemical status and biological integrity of the Upper Mara River basin. This was achieved by examining trends in nutrient concentrations and analyzing the structure, diversity and abundance of benthic macro-invertebrates in relation to varying land use patterns. Sampling sites were selected based on catchment land use and the level of human disturbance, and using historical records of previous water quality studies. River water pH, dissolved oxygen, electrical conductivity (EC), temperature, and turbidity were determined in situ. All investigated parameters except iron and manganese had concentration values within allowable limits according to Kenyan and international standards for drinking water. The Amala tributary is more mineralized and also shows higher levels of pH and EC than water from the Nyangores tributary. The latter, however, has a higher variability in both the total phosphorus (TP) and total nitrogen (TN) concentrations. The variability in TP and TN concentrations increases downstream for both tributaries and is more pronounced for TN than for TP. Macro-invertebrate assemblages responded to the changes in land use and water quality in terms of community composition and diversity. The study recommends detailed continuous monitoring of the water quality at shorter time intervals and to identify key macro-invertebrate taxa that can be used to monitor changes of the water quality in rivers of the Mara basin as a result of anthropogenic changes.     

Numerical modelling of the impacts of water abstraction for hydraulic fracturing on groundwater–surface water interaction: a case study from northwestern Alberta,Canada

ABSTRACT

The temporal dynamics of groundwater–surface water interaction under the impacts of various water abstraction scenarios are presented for hydraulic fracturing in a shale gas and oil play area (23 984.9 km²), Alberta, Canada, using the MIKE-SHE and MIKE-11 models. Water-use data for hydraulic fracturing were obtained for 433 wells drilled in the study area in 2013 and 2014. Modelling results indicate that water abstraction for hydraulic fracturing has very small (<0.35%) negative impacts on mean monthly and annual river and groundwater levels and stream and groundwater flows in the study area, and small (1–4.17%) negative impacts on environmental flows near the water abstraction location during low-flow periods. The impacts on environmental flow depend on the amount of water abstraction and the daily flow over time at a specific river cross-section. The results also indicate a very small (<0.35%) positive impact on mean monthly and annual groundwater contributions to streamflow because of the large study area. The results provide useful information for planning long-term seasonal and annual water abstractions from the river and groundwater for hydraulic fracturing in a large study area.     

Assessment of variables controlling nitrate dynamics in groundwater: is it a threat to surface aquatic ecosystems?

The impact of fertilised cropping on nitrate-N dynamics in groundwater (GW) was assessed in a catchment from piezometers installed: (i) to different depths, (ii) in different soil types, (iii) on different positions on landscape, and (iv) compared with the Australian and New Zealand Environmental and Conservation Council guideline values provided for different aquatic ecosystems. The GW and NO(3)-N concentration dynamics were monitored in 39 piezometer wells, installed to 5-90 m depth, under fertilized sugarcane (Saccharum officinarum-S) in the Johnstone River Catchment, Australia, from 1999 January through September 2002. The median nitrate-N concentration ranged from 14 to 1511 microg L(-1), and the 80th percentile from 0 to 1341 microg L(-1). In 34 out of the 39 piezometer wells the 80th percentile or 80% of the nitrate-N values were higher than 30 microg L(-1), which is the maximum trigger value provided in the ANZECC table for sustainable health of different aquatic ecosystems. Nitrate-N concentration decreased with increasing well depth, increasing depth of water in wells, and with decreasing relief on landscape. Nitrate-N was higher in alluvial soil profiles than on those formed in-situ. Nitrate-N increased with increasing rainfall at the beginning of the rainy season, fluctuated during the peak rainy period, and then decreased when the rain ceased. The rapid decrease in GW after the rains ceased suggested potential existed for nitrate-N to be discharged as lateral-flow into streams. This may contribute towards the deterioration in the health of down-stream aquatic ecosystems.     

Assessment of future water resources and water scarcity considering the factors of climate change and social–environmental change in Han River basin,Korea

Water resources are influenced by various factors such as weather, topography, geology, and environment. Therefore, there are many difficulties in evaluating and analyzing water resources for the future under climate change. In this paper, we consider climate, land cover and water demand as the most critical factors affecting change in future water resources. We subsequently introduce the procedures and methods employed to quantitatively evaluate the influence of each factor on the change in future water resources. In order to consider the change in land cover, we apply the Multi-Regression approach from the cellular automata-Markov Chain technique using two independent variables, temperature and rainfall. In order to estimate the variation of the future runoff due to climate change, the data of the SRES A2 climate change scenario were entered in the SLURP model to simulate a total of 70 years, 2021–2090, of future runoff in the Han River basin in Korea. However, since a significant amount of uncertainties are involved in predicting the future runoff due to climate change, 50 sets of daily precipitation data from the climate change scenario were generated and used for the SLURP model to forecast 50 sets of future daily runoff. This process was used to minimize the uncertainty that may occur when the prediction process is performed. For future water balance analysis, the future water demand was divided into low demand, medium demand and high demand categories. The three water demand scenarios and the 50 daily runoff scenarios were combined to form 150 sets of input data. The monthly water balance within the Han River basin was then calculated using this data and the Korean version of Water Evaluation and Planning System model. As a result, the future volume of water scarcity of the Han River basin was predicted to increase in the long term. It is mostly due to the monthly shift in the runoff characteristic, rather than the change in runoff volume resulting from climate change.     

Geochemical tracing and modeling of surface and deep water–rock interactions in elementary granitic watersheds(Strengbach and Ringelbach CZOs,France)

From the study of the Strengbach and Ringelbach watersheds we propose to illustrate the interest of combining the geochemical tracing and geochemical modeling approaches on surface and deep borehole waters,to decipher the diversity of the water flow and the associated water–rock interactions in such elementary mountainous catchments. The results point to a clear geochemical typology of waters depending on the water circulations(deep vs. hypodermic) within the substratum.     

Surface water–groundwater interactions of Xiluodu Reservoir based on the dynamic evolution of seepage,temperature, and hydrochemistry due to impoundment

Reservoir construction greatly affects the regional ecological environment, particularly surface water–groundwater interactions around the reservoir. Xiluodu Reservoir, a representative large-scale reservoir in China, has had substantial impacts on surface water–groundwater interactions at the dam site since impoundment. This study analysed the dynamic characteristics of surface water–groundwater level, temperature, and hydrochemistry to determine the evolution of surface water–groundwater interaction before and after the impoundment. The levels of groundwater and some surface water rose by more than 100 m after impoundment and the water level of saturated limestone gradually stabilized, whereas basalt saturation in the affected area continued to expand. The groundwater temperature did not decrease significantly, whereas the hydrochemical types and ion contents of both surface water and groundwater experienced significant changes. Calculation of the saturation index indicated spatiotemporal changes in the saturation state of minerals. The replenishment source of each type of water and their mutual relationships were determined using cluster analysis and isotope characteristics. The results confirmed continuous, significant, and variable surface water–groundwater interactions at the dam site, which were partially reversed after impoundment. Changes in surface water–groundwater interactions were due to impoundment, the impact of which decreases with distance from the dam, as well as the unique geological conditions and artificial construction.     

FluBiDi – a model to estimate flood based on runoff: validation using extreme and natural basin conditions

FluBiDi is a two-dimensional model created to simulate real events that can take days and months, as well as short events (minutes or hours) and inclusive laboratory tests. To verify the robustness of FluBiDi, it was tested using a previous study with both designed and real digital elevation models. The results highlight good agreement between the models (i.e. Mike Flood, SOBEK, ISIS 2D, and others) tested and FluBiDi (around 90% for a specific instant and 95% for the complete time simulation). In the simulated hydrographs, the discharge peak value, time to peak, and water level results were accurate, reproducing them with an error of less than 5%. The velocity differences observed in a couple of tests in FluBiDi were associated with very short periods of time (seconds). However, FluBiDi is highly accurate for simulating floods under real topographical conditions with differences of around 2 cm when water depth is around 150 cm. The average water depth and velocities are precise, and the model describes with high accuracy the pattern and extent of floods. FluBiDi has the capability to be adjusted to different types of events and only requires limited input data.     

Uncertainty estimation of a complex water quality model: The influence of Box–Cox transformation on Bayesian approaches and comparison with a non-Bayesian method

In urban drainage modelling, uncertainty analysis is of undoubted necessity. However, uncertainty analysis in urban water-quality modelling is still in its infancy and only few studies have been carried out. Therefore, several methodological aspects still need to be experienced and clarified especially regarding water quality modelling. The use of the Bayesian approach for uncertainty analysis has been stimulated by its rigorous theoretical framework and by the possibility of evaluating the impact of new knowledge on the modelling predictions. Nevertheless, the Bayesian approach relies on some restrictive hypotheses that are not present in less formal methods like the Generalised Likelihood Uncertainty Estimation (GLUE). One crucial point in the application of Bayesian method is the formulation of a likelihood function that is conditioned by the hypotheses made regarding model residuals. Statistical transformations, such as the use of Box–Cox equation, are generally used to ensure the homoscedasticity of residuals. However, this practice may affect the reliability of the analysis leading to a wrong uncertainty estimation. The present paper aims to explore the influence of the Box–Cox equation for environmental water quality models. To this end, five cases were considered one of which was the “real” residuals distributions (i.e. drawn from available data). The analysis was applied to the Nocella experimental catchment (Italy) which is an agricultural and semi-urbanised basin where two sewer systems, two wastewater treatment plants and a river reach were monitored during both dry and wet weather periods. The results show that the uncertainty estimation is greatly affected by residual transformation and a wrong assumption may also affect the evaluation of model uncertainty. The use of less formal methods always provide an overestimation of modelling uncertainty with respect to Bayesian method but such effect is reduced if a wrong assumption is made regarding the residuals distribution. If residuals are not normally distributed, the uncertainty is over-estimated if Box–Cox transformation is not applied or non-calibrated parameter is used.     

Complexity analysis of precipitation using the Lempel–Ziv algorithm and a multi-scaling approach: a case study in Jilin province,China

Precipitation is an important part of the hydrologic cycle, and its complexity is closely related to surface runoff and changing groundwater dynamics, which in turn influences the accuracy of precipitation forecasts. In this study, we used the Lempel–Ziv algorithm (LZA) and a multi-scaling approach to assess precipitation complexity for 1958–2011 by analyzing time series data from 28 gauging stations located throughout Jilin province, China. The spatial distribution of normalized precipitation complexity was measured by LZA, a symbolic dynamics algorithm, and by a multi-scaling approach, which is described by fractals. In addition, the advantages and limitations of these two methods were investigated. The results indicate that both methods are applicable and consistent for calculating precipitation complexity, and that the degree of relief is a primary factor controlling precipitation complexity in the mountainous area; in the plain terrain, however, the prominent influencing factor is climate.     

Impact of different sea surface roughness on surface gravity waves using a coupled atmosphere–wave model: a case of Hurricane Isaac (2012)

Ocean Dynamics - The complicated pattern of the chaotic ocean surface depends strongly on the interaction between wind and waves. An accurate representation of momentum and energy exchange at...     

Assessment of calanchi and rill–interrill erosion susceptibilities using terrain analysis and geostochastics: A case study in the Oltrepo Pavese,Northern Apennines,Italy

Soil erosion is one of the most important environmental problems distributed worldwide. In the last decades, numerous studies have been published on the assessment of soil erosion and the related processes and forms using empirical, conceptual and physically based models. For the prediction of the spatial distribution, more and more sophisticated stochastic modelling approaches have been proposed – especially on smaller spatial scales such as river basins. In this work, we apply a maximum entropy model (MaxEnt) to evaluate badlands (calanchi) and rill–interrill (sheet erosion) areas in the Oltrepo Pavese (Northern Apennines, Italy). The aim of the work is to assess the important environmental predictors that influence calanchi and rill–interrill erosion at the regional scale. We used 13 topographic parameters derived from a 12 m digital elevation model (TanDEM-X) and data on the lithology and land use. Additional information about the vegetation is introduced through the normalized difference vegetation index based on remotely sensed data (ASTER images). The results are presented in the form of susceptibility maps showing the spatial distribution of the occurrence probability for calanchi and rill–interrill erosion. For the validation of the MaxEnt model results, a support vector machine approach was applied. The models show reliable results and highlight several locations of the study area that are potentially prone to future soil erosion. Thus, coping and mitigation strategies may be developed to prevent or fight the soil erosion phenomenon under consideration. © 2020 John Wiley & Sons, Ltd.