A geomorphology‐based integrated stream–aquifer interaction model for semi‐gauged catchments

Many of the existing stream–aquifer interaction models available in the literature are very complex with limited applicability in semi‐gauged and ungauged catchments. In this study, to estimate the influent and effluent subsurface water fluxes under limited geo‐hydrometeorological data availability conditions, a simple stream–aquifer interaction model, namely, the variable parameter McCarthy–Muskingum (VPMM) hillslope‐storage Boussinesq (hsB) model, has been developed. This novel model couples the VPMM streamflow transport with the hsB groundwater flow transport modules in online mode. In this integrated model, the surface water–groundwater flux exchange process is modelled by the Darcian approach with the variable hydraulic heads between the river stage and groundwater table accounting for the rainfall forcing. Considering the exchange fluxes in the hyporheic zone and lateral overland flow contribution, this approach is field tested in a typical 48‐km stretch of the Brahmani River in eastern India to simulate the streamflow and its depth with the minimum Nash–Sutcliffe efficiency of 94% and 88%; the maximum root mean square error of 134 m³/s and 0.35 m; and the minimum index of agreement of 98% and 97%, respectively. This modelling approach could be very well utilized in data‐scarce world‐river basins to estimate the stream–aquifer exchange flux due to rainfall forcings.     

Modeling surface water–groundwater interaction in New Zealand: Model development and application

Most rivers worldwide have a strong interaction with groundwater when they leave the mountains and flow over alluvial plains before flowing into the seas or disappearing in the deserts, and in New Zealand, typically, rivers lose water to the groundwater in the upper plains and generally gain water from the groundwater in the lower plains. Aiming at simulating surface water–groundwater interaction nationally in New Zealand, we developed a conceptual groundwater module for the national hydrologic model TopNet to simulate surface water–groundwater interaction, groundwater flow, and intercatchment groundwater flow. The developed model was applied to the Pareora catchment in South Island of New Zealand, where there are concurrent spot gauged flows. Results show that the model simulations not only fit quite well to flow measurement but also to concurrent spot gauged flows, and compared to the original TopNet, it has a significant improvement in the low flows. Sensitivity analysis shows river flow is sensitive to the river losing/gaining rate instead of groundwater characteristic, while groundwater storage is sensitive to both river losing/gaining rate and groundwater characteristic. This indicates our conceptual approach is promising for nationwide modeling without the large amount of geology and aquifer data typically required by physically‐based modeling approaches.     

Analytical and statistical approach for evaluating the effects of a river barrage on river–aquifer interactions

The construction of a river barrage can increase groundwater levels upstream of the barrage during the rainy season. Analytical and statistical approaches were applied to evaluate the relationship between groundwater and river water at the Changnyeong–Haman river barrage in Korea using time series data of water level and electrical conductivity from June 2011 to September 2014. An artificial neural network based time series model was designed to filter out the effect of rainfall from the groundwater level data in the study area. Aquifer diffusivity and river resistance were estimated from the analytical solution of a one‐dimensional unit step response function by using the filtered groundwater level data. River resistance increased in response to groundwater level fluctuations. Cross‐correlation analyses between the groundwater and the river water showed that the lag time increased during the observation period for both the water level and the electrical conductivity while the cross‐correlation function declined for the same period. The results indicated that a constant river stage maintained at the river barrage can weaken the hydrologic stress and reduce the exchange of material between the river and the adjacent aquifer because of the deposition of fine sediment on the river bottom and walls. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

River–aquifer interactions in a semi‐arid environment stressed by groundwater abstraction

Rivers and aquifers are, in many cases, a connected resource and as such the interactions between them need to be understood and quantified for the resource to be managed appropriately. The objective of this paper is to advance the understanding of river–aquifer interactions processes in semi‐arid environments stressed by groundwater abstraction. This is performed using data from a specific catchment where records of precipitation, evapotranspiration, river flow, groundwater levels and groundwater abstraction are analysed using basic statistics, hydrograph analysis and a simple mathematical model to determine the processes causing the spatial and temporal changes in river–aquifer interactions. This combined approach provides a novel but simple methodology to analyse river–aquifer interactions, which can be applied to catchments worldwide. The analysis revealed that the groundwater levels have declined (~ 3 m) since the onset of groundwater abstraction. The decline is predominantly due to the abstraction rather than climatic changes (r = 0.84 for the relationship between groundwater abstraction and groundwater levels; r = 0.92 for the relationship between decline in groundwater levels and magnitude of seasonal drawdown). It is then demonstrated that, since the onset of abstraction, the river has changed from being gaining to losing during low‐flow periods, defined as periods with flow less than 0.5, 1.0 or 1.5 GL/day (1 GL/day = 1 × 10⁶ m³/day). If defined as < 1.0 GL/day, low‐flow periods constitute approximately 65% of the river flows; the periods where the river is losing at low‐flow conditions are thus significant. Importantly, there was a significant delay (> 10 years) between the onset of groundwater abstraction and the changeover from gaining to losing conditions. Finally, a relationship between the groundwater gradient towards the river and the river flow at low‐flow is demonstrated. The results have important implications for water management as well as water ecology and quality. Copyright © 2012 John Wiley & Sons, Ltd.     

First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach

Distributed hydrologic models capable of simulating fully‐coupled surface water and groundwater flow are increasingly used to examine problems in the hydrologic sciences. Several techniques are currently available to couple the surface and subsurface; the two most frequently employed approaches are first‐order exchange coefficients (a.k.a., the surface conductance method) and enforced continuity of pressure and flux at the surface‐subsurface boundary condition. The effort reported here examines the parameter sensitivity of simulated hydrologic response for the first‐order exchange coefficients at a well‐characterized field site using the fully coupled Integrated Hydrology Model (InHM). This investigation demonstrates that the first‐order exchange coefficients can be selected such that the simulated hydrologic response is insensitive to the parameter choice, while simulation time is considerably reduced. Alternatively, the ability to choose a first‐order exchange coefficient that intentionally decouples the surface and subsurface facilitates concept‐development simulations to examine real‐world situations where the surface‐subsurface exchange is impaired. While the parameters comprising the first‐order exchange coefficient cannot be directly estimated or measured, the insensitivity of the simulated flow system to these parameters (when chosen appropriately) combined with the ability to mimic actual physical processes suggests that the first‐order exchange coefficient approach can be consistent with a physics‐based framework. Copyright © 2009 John Wiley & Sons, Ltd.     

Hydrochemical variations during flood pulses in the south‐west China peak cluster karst: impacts of CaCO3–H2O–CO2 interactions

High‐resolution measurements of rainfall, water level, pH, conductivity, temperature and carbonate chemistry parameters of groundwater at two adjacent locations within the peak cluster karst of the Guilin Karst Experimental Site in Guangxi Province, China, were made with different types of multiparameter sonde. The data were stored using data loggers recording with 2 min or 15 min resolution. Waters from a large, perennial spring represent the exit for the aquifer's conduit flow, and a nearby well measures water in the conduit‐adjacent, fractured media. During flood pulses, the pH of the conduit flow water rises as the conductivity falls. In contrast, and at the same time, the pH of groundwater in the fractures drops, as conductivity rises. As Ca²⁺ and HCO₃^? were the dominant (>90%) ions, we developed linear relationships (both r² > 0·91) between conductivity and those ions, respectively, and in turn calculated variations in the calcite saturation index (SI_C) and CO₂ partial pressure (P) of water during flood pulses. Results indicate that the P of fracture water during flood periods is higher than that at lower flows, and its SI_C is lower. Simultaneously, P of conduit water during the flood period is lower than that at lower flows, and its SI_C also is lower. From these results we conclude that at least two key processes are controlling hydrochemical variations during flood periods: (i) dilution by precipitation and (ii) water–rock–gas interactions. To explain hydrochemical variations in the fracture water, the water–rock–gas interactions may be more important. For example, during flood periods, soil gas with high CO₂ concentrations dissolves in water and enters the fracture system, the water, which in turn has become more highly undersaturated, dissolves more limestone, and the conductivity increases. Dilution of rainfall is more important in controlling hydrochemical variations of conduit water, because rainfall with higher pH (in this area apparently owing to interaction with limestone dust in the lower atmosphere) and low conductivity travels through the conduit system rapidly. These results illustrate that to understand the hydrochemical variations in karst systems, considering only water–rock interactions is not sufficient, and the variable effects of CO₂ on the system should be evaluated. Consideration of water–rock–gas interactions is thus a must in understanding variations in karst hydrochemistry. Copyright © 2004 John Wiley & Sons, Ltd.     

Identification of groundwater recharge sources and processes in a heterogeneous alluvial aquifer: results from multi‐level monitoring of hydrochemistry and environmental isotopes in a riverside agricultural area in Korea

We evaluated sources and pathways of groundwater recharge for a heterogeneous alluvial aquifer beneath an agricultural field, based on multi‐level monitoring of hydrochemistry and environmental isotopes of a riverside groundwater system at Buyeo, Korea. Two distinct groundwater zones were identified with depth: (1) a shallow oxic groundwater zone, characterized by elevated concentrations of NO₃^? and (2) a deeper (>10–14 m from the ground surface) sub‐oxic groundwater zone with high concentrations of dissolved Fe, silica, and HCO₃^?, but little nitrate. The change of redox zones occurred at a depth where the aquifer sediments change from an upper sandy stratum to a silty stratum with mud caps. The δ¹⁸O and δ²H values of groundwater were also different between the two zones. Hydrochemical and δ¹⁸O? δ²H data of oxic groundwater are similar to those of soil water. This illustrates that recharge of oxic groundwater mainly occurs through direct infiltration of rain and irrigation water in the sandy soil area where vegetable cropping with abundant fertilizer use is predominant. Oxic groundwater is therefore severely contaminated by agrochemical pollutants such as nitrate. In contrast, deeper sub‐oxic groundwater contains only small amounts of dissolved oxygen (DO) and NO₃^?. The ³H contents and elevated silica concentrations in sub‐oxic groundwater indicate a somewhat longer mean residence time of groundwater within this part of the aquifer. Sub‐oxic groundwater was also characterized by higher δ¹⁸O and δ²H values and lower d‐excess values, indicating significant evaporation during recharge. We suggest that recharge of sub‐oxic groundwater occurs in the areas of paddy rice fields where standing irrigation and rain water are affected by strong evaporation, and that reducing conditions develop during subsequent sub‐surface infiltration. This study illustrates the existence of two groundwater bodies with different recharge processes within an alluvial aquifer. Copyright © 2009 John Wiley & Sons, Ltd.     

Inter‐annual variations and factors controlling evapotranspiration in a temperate Japanese cypress forest

This study was undertaken to evaluate the effects of climatic variability on inter‐annual variations in each component of evapotranspiration (ET) and the total ET in a temperate coniferous forest in Japan. We conducted eddy covariance flux and meteorological measurements for 7 years and parameterized a one‐dimensional multi‐layer biosphere‐atmosphere model (Kosugi et al., 2006 ) that partitions ET to transpiration (Tr), wet‐canopy evaporation (E_wet), and soil evaporation (E_soil). The model was validated with the observed flux data. Using the model, the components of ET were estimated for the 7 years. Annual precipitation, ET, Tr, E_wet, and E_soil over the 7 years were 1536 ± 334 mm, 752 ± 29 mm, 425 ± 37 mm, 219 ± 34 mm, and 108 ± 10 mm, respectively. The maximum inter‐annual fluctuation of observed ET was 64 mm with a coefficient of variance (CV) of 2.7%, in contrast to relatively large year‐to‐year variations in annual rainfall (CV = 20.1%). Tr was related to the vapour pressure deficit, incoming radiation, and air temperature with relatively small inter‐annual variations (CV = 8.2%). E_soil (CV = 8.6%) was related mainly to the vapour pressure deficit. E_wet was related to precipitation with large inter‐annual variations (CV = 14.3%) because of the variability in precipitation. The variations in E_wet were counterbalanced by the variations in Tr and E_soil, producing the small inter‐annual variations in total ET. Copyright © 2016 John Wiley & Sons, Ltd.     

Implementation of a second‐order paraxial boundary condition for a water‐saturated layered half‐space in plane strain

A half‐space finite element and a transmitting boundary are developed for a water‐saturated layered half‐space using a paraxial boundary condition. The exact dynamic stiffness of a half‐space in plane strain is derived and a second‐order paraxial approximation of the stiffness is obtained. A half‐space finite element and a transmitting boundary are then formulated. The development is verified by comparison of the dynamic stiffness of impermeable and permeable rigid strip foundations with other published results. The advantage of using the paraxial boundary condition in comparison with the rigid boundary condition is examined. It is shown that the paraxial boundary condition offers significant gain and the resulting half‐space finite element and transmitting boundary can represent the effects of a water‐saturated layered half‐space with good accuracy and efficiency. In addition, the numerical method described herein maintains the strengths and advantages of the finite element method and can be easily applied to demanding problems of soil–structure interaction in a water‐saturated layered half‐space. Copyright © 2010 John Wiley & Sons, Ltd.     

Groundwater–surface water interactions in a large semi‐arid floodplain: implications for salinity management

Flow regulation and water diversion for irrigation have considerably impacted the exchange of surface water between the Murray River and its floodplains. However, the way in which river regulation has impacted groundwater–surface water interactions is not completely understood, especially in regards to the salinization and accompanying vegetation dieback currently occurring in many of the floodplains. Groundwater–surface water interactions were studied over a 2 year period in the riparian area of a large floodplain (Hattah–Kulkyne, Victoria) using a combination of piezometric surface monitoring and environmental tracers (Cl⁻, δ²H, and δ¹⁸O). Despite being located in a local and regional groundwater discharge zone, the Murray River is a losing stream under low flow conditions at Hattah–Kulkyne. The discharge zone for local groundwater, regional groundwater and bank recharge is in the floodplain within ∼1 km of the river and is probably driven by high rates of transpiration by the riparian Eucalyptus camaldulensis woodland. Environmental tracers data suggest that the origin of groundwater is principally bank recharge in the riparian zone and a combination of diffuse rainfall recharge and localized floodwater recharge elsewhere in the floodplain. Although the Murray River was losing under low flows, bank discharge occurred during some flood recession periods. The way in which the water table responded to changes in river level was a function of the type of stream bank present, with point bars providing a better connection to the alluvial aquifer than the more common clay‐lined banks. Understanding the spatial variability in the hydraulic connection with the river channel and in vertical recharge following inundations will be critical to design effective salinity remediation strategies for large semi‐arid floodplains. Copyright © 2005 John Wiley & Sons, Ltd.