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.     

Assessing the ‘two water worlds’ hypothesis and water sources for native and exotic evergreen species in south‐central Chile

Recent studies using water‐stable isotopes (δ¹⁸O and δ²H) have suggested an ecohydrological separation of water flowing to streams or recharging groundwater and water used by trees, known as the ‘two water worlds’ (TWW) hypothesis. In this study, we measured water isotopic composition in precipitation [open field and throughfall, i.e. local meteoric water line (LMWL)] and the mobile water compartment (i.e. stream and soil solution), bulk soil water and xylem water over a period of 1.5 years in two headwater catchments: NF, covered with old growth native evergreen forest (Aetoxicon punctatum, Laureliopsis philippiana and Eucriphya cordifolia), and EP, covered with 4 and 16‐year‐old Eucalyptus nitens stands. Our results show that precipitation, stream and soil solution plot approximately along the LMWL, while xylem waters from all studied tree species plot below the LMWL, supporting the TWW hypothesis. However, we also found evidence of ecohydrological connectivity during the wet season, likely controlled by the amount of antecedent precipitation. These observations hold for all investigated tree species. On both sites, a different precipitation source for stream and xylem water was observed. However, in EP, bulk soil showed a similar precipitation source as xylem water from both E. nitens stands. This suggests that E. nitens may use water that is recharging the bulk soil compartment. We conclude that under a rainy temperate climate, the TWW hypothesis is temporal and does not apply during wet seasons. Copyright © 2016 John Wiley & Sons, Ltd.     

Energy partitioning over a semi‐arid shrubland in northern China

During the last decade, the widely distributed shrublands in northern China have shown significant signs of recovery from desertification, the result of widespread conservation practices. However, to support the current efforts in conservation, more knowledge is needed on surface energy partitioning and its biophysical controls. Using eddy‐covariance measurements made over a semi‐arid shrubland in northwest China in 2012, we examined how surface energy‐balance components vary on diurnal and seasonal scales, and how biophysical factors control bulk surface parameters and energy exchange. Sensible heat flux (H) exceeded latent heat flux (λE) during most of the year, resulting in an annual Bowen ratio (β, i.e. H/λE) of 2.0. λE exceeded H only in mid‐summer when frequent rainfall co‐occurred with the seasonal peak in leaf area index (LAI). Evapotranspiration reached a daily maximum of 3.3 mm day^?1, and summed to 283 mm yr^?1. The evaporative fraction (EF, i.e. λE/R_n), Priestley–Taylor coefficient (α), surface conductance (g_s) and decoupling coefficient (Ω) were all positively correlated with soil water content (SWC) and LAI. The direct enhancement of λE by high vapour pressure deficit (VPD) was buffered by a concurrent suppression of g_s. The g_s played a direct role in controlling EF and α by mediating the effects of LAI, SWC and VPD. Our results highlight the importance of adaptive plant responses to water scarcity in regulating ecosystem energy partitioning, and suggest an important role for revegetation in the reversal of desertification in semi‐arid areas. Copyright © 2015 John Wiley & Sons, Ltd.     

A new approach to estimate canopy evaporation and canopy interception capacity from evapotranspiration and sap flow measurements during and following wetting

Canopy interception and its evaporation into the atmosphere during irrigation or a rainfall event are important in irrigation scheduling, but are challenging to estimate using conventional methods. This study introduces a new approach to estimate the canopy interception from measurements of actual total evapotranspiration (ET) using eddy covariance and estimation of the transpiration from measurements of sap flow. The measurements were conducted over a small‐scale sprinkler‐irrigated cotton field before, during and after sprinkler irrigation. Evaporation and sap flow dynamics during irrigation show that the total ET during irrigation increased significantly because of the evaporation of free intercepted water while transpiration was suppressed almost completely. The difference between actual ET and transpiration (sap flow) during and immediately following irrigation (post irrigation) represents the total canopy evaporation while the canopy interception capacity was calculated as the difference between actual ET and transpiration (sap flow) during drying (post irrigation) following cessation of the irrigation. The canopy evaporation of cotton canopy was calculated as 0.8 mm, and the interception capacity was estimated to be 0.31 mm of water. The measurement uncertainty in both the non‐dimensional ET and non‐dimensional sap flow was shown to be very low. Copyright © 2015 John Wiley & Sons, Ltd.     

Dissolved inorganic carbon evolution in neutral discharge from mine tailings piles

We measured the concentrations of dissolved inorganic carbon (DIC) and major ions and the stable carbon isotope ratios of DIC (δ¹³C_DIC) in two creeks discharging from carbonate‐rich sulphide‐containing mine tailings piles. Our aim was to assess downstream carbon evolution of the tailings discharge as it interacted with the atmosphere. The discharge had pH of 6.5–8.1 and was saturated with respect to carbonates. Over the reach of one creek, the DIC concentrations decreased by 1.1 mmol C/l and δ¹³C_DIC increased by ~4.0‰ 200 m from the seep source. The decrease in the DIC concentrations was concomitant with decreases in the partial pressure of CO_2(aq) because of the loss of excess CO_2(aq) from the discharge. The corresponding enrichment in the δ¹³C_DIC is because of kinetic isotope fractionation accompanying the loss of CO_2(g). Over the reach of the other creek, there was no significant decrease in the DIC concentrations or notable changes in the δ¹³C_DIC. The insignificant change in the DIC concentrations and the δ¹³C_DIC is because the first water sample was collected 160 m away from the discharge seep, not accessible during this research. In this case, most of the excess CO_2(aq) was lost before our first sampling station. Our results indicate that neutral discharges from tailings piles quickly lose excess CO_2(aq) to the atmosphere and the DIC becomes enrich in ¹³C. We suggest that a significant amount of carbon cycling in neutral discharges from tailings piles occur close to the locations where the discharge seeps to the surface. Copyright © 2015 John Wiley & Sons, Ltd.     

Modelling atmospheric and hydrologic processes for assessment of meadow restoration impact on flow and sediment in a sparsely gauged California watershed

The restoration of meadowland using the pond and plug technique of gully elimination was performed in a 9‐mile segment along Last Chance Creek, Feather River Basin, California, in order to rehabilitate floodplain functions such as mitigating floods, retaining groundwater, and reducing sediment yield associated with bank erosion and to significantly alter the hydrologic regime. However, because the atmospheric and hydrological conditions have evolved over the restoration period, it was difficult to obtain a comprehensible evaluation of the impact of restoration activities by means of field measurements. In this paper, a new use of physically based models for environmental assessment is described. The atmospheric conditions over the sparsely gauged Last Chance Creek watershed (which does not have any precipitation or weather stations) during the combined historical critical dry and wet period (1982–1993) were reconstructed over the whole watershed using the atmospheric fifth‐generation mesoscale model driven with the US National Center for Atmospheric Research and US National Center for Environmental Prediction reanalysis data. Using the downscaled atmospheric data as its input, the watershed environmental hydrology (WEHY) model was applied to this watershed. All physical parameters of the WEHY model were derived from the existing geographic information system and satellite‐driven data sets. By comparing the prerestoration and postrestoration simulation results under the identical atmospheric conditions, a more complete environmental assessment of the restoration project was made. Model results indicate that the flood peak may be reduced by 10–20% during the wet year and the baseflow may be enhanced by 10–20% during the following dry seasons (summer to fall) in the postrestoration condition. The model results also showed that the hydrologic impact of the land management associated with the restoration mitigates bank erosion and sediment discharge during winter storm events. Copyright © 2013 John Wiley & Sons, Ltd.     

Integrating isolated and riparian wetland modules in the PHYSITEL/HYDROTEL modelling platform: model performance and diagnosis

Mathematical modelling is a well‐accepted framework to evaluate the effects of wetlands on stream flow and watershed hydrology in general. Although the integration of wetland modules into a distributed hydrological model represents a cost‐effective way to make this assessment, the added value brought by landscape‐specific modules to a model's ability to replicate basic hydrograph characteristics remains unclear. The objectives of this paper were the following: (i) to present the adaptation of PHYSITEL (a geographic information system) to parameterize isolated and riparian wetlands; (ii) to describe the integration of specific isolated wetland and riparian wetland modules into HYDROTEL, a distributed hydrological model; and (iii) to evaluate the performance of the updated modelling platform with respect to the capacity of replicating various hydrograph characteristics. To achieve this, two sets of simulations were performed (with and without wetland modules), and the added value was assessed at three river segments of the Becancour River watershed, Quebec, Canada, using six general goodness‐of‐fit indicators and 14 water flow criteria. A sensitivity analysis of the wetland module parameters was performed to characterize their impact on stream flows of the modelled watershed. Results of this study indicate the following: (i) integration of specific wetland modules can slightly increase the capacity of HYDROTEL to replicate basic hydrograph characteristics; and (ii) the updated modelling platform allows for the explicit assessment of the impact of wetlands (e.g. typology and location) on watershed hydrology. Copyright © 2015 John Wiley & Sons, Ltd.