Scales and magnitude of hyporheic,river–aquifer and bank storage exchange fluxes

Many studies have investigated the exchange processes that occur between rivers and groundwater systems and have successfully quantified the water fluxes involved. Specifically, these exchange processes include hyporheic exchange, river–aquifer exchange (groundwater discharge and river loss) and bank storage exchange. Remarkably, there are relatively few examples of field studies where more than one exchange process is quantified, and as a consequence, the relationships between them are not well understood. To compare the relative magnitudes of these common exchange processes, we have collected data from 54 studies that have quantified one or more of these exchange flux types. Each flux value is plotted against river discharge at the time of measurement to allow the different exchange flux types to be compared. We show that there are positive relationships between the magnitude of each exchange flux type and increasing river discharge across the different studies. For every one order of magnitude increase in river discharge, the hyporheic, river–aquifer and bank storage exchange fluxes increase by factors of 2.7, 2.9 and 2.5, respectively. On average, hyporheic exchange fluxes are almost an order of magnitude greater than river–aquifer exchange fluxes, which are, in turn, approximately four times greater than bank storage exchange fluxes for the same river discharge. Unless measurement approaches that can distinguish between different types of exchange flux are used, there is potential for hyporheic exchange fluxes to be misinterpreted as river–aquifer exchange fluxes, with possible implications for water resource management decisions. Copyright © 2015 John Wiley & Sons, Ltd.     

A smoothed cap model for variably saturated soils and its robust numerical implementation

This paper deals with the numerical implementation of a cap model for unsaturated soils. It provides a brief review of existing cap model approaches, based on which an improved model formulated in terms of generalised effective stress and matric suction is derived and described in detail. Although the proposed model is a multisurface plasticity model, it can efficiently be implemented using only single‐surface projections because of the smoothness of the model, which is obtained by construction. Numerical algorithms are provided for these single‐surface stress projections, using a single‐equation approach whenever possible. The robustness of the utilised single‐equation approaches is enhanced by proposing problem‐fitted start‐up procedures based on investigations of the nonlinear projection equations. A comparison of the model response with extensive material test data is used to validate the model and to demonstrate the robust application of the approach to silty sands and low to medium plasticity clays. Copyright © 2015 John Wiley & Sons, Ltd.     

Improved spatial delineation of streambed properties and water fluxes using distributed temperature sensing

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre‐optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5‐day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one‐dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one‐dimensional analytical heat‐tracing models. The results show that the amplitude ratio‐based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat‐tracing studies by providing a technique for visualizing and comparing results from fibre‐optic distributed temperature sensing installations and testing the robustness of analytical heat‐tracing models. Copyright © 2016 John Wiley & Sons, Ltd.     

Modelling groundwater/surface water interaction in a managed riparian chalk valley wetland

Understanding hydrological processes in wetlands may be complicated by management practices and complex groundwater/surface water interactions. This is especially true for wetlands underlain by permeable geology, such as chalk. In this study, the physically based, distributed model MIKE SHE is used to simulate hydrological processes at the Centre for Ecology and Hydrology River Lambourn Observatory, Boxford, Berkshire, UK. This comprises a 10‐ha lowland, chalk valley bottom, riparian wetland designated for its conservation value and scientific interest. Channel management and a compound geology exert important, but to date not completely understood, influences upon hydrological conditions. Model calibration and validation were based upon comparisons of observed and simulated groundwater heads and channel stages over an equally split 20‐month period. Model results are generally consistent with field observations and include short‐term responses to events as well as longer‐term seasonal trends. An intrinsic difficulty in representing compressible, anisotropic soils limited otherwise excellent performance in some areas. Hydrological processes in the wetland are dominated by the interaction between groundwater and surface water. Channel stage provides head boundaries for broad water levels across the wetland, whilst areas of groundwater upwelling control discrete head elevations. A relic surface drainage network confines flooding extents and routes seepage to the main channels. In‐channel macrophyte growth and its management have an acute effect on water levels and the proportional contribution of groundwater and surface water. The implications of model results for management of conservation species and their associated habitats are discussed. Copyright © 2015 John Wiley & Sons, Ltd.     

Generation of continuous surface soil moisture dataset using combined optical and thermal infrared images

Surface soil moisture (SSM) is a critical variable for understanding water and energy flux between the atmosphere and the Earth's surface. An easy to apply algorithm for deriving SSM time series that primarily uses temporal parameters derived from simulated and in situ datasets has recently been reported. This algorithm must be assessed for different biophysical and atmospheric conditions by using actual geostationary satellite images. In this study, two currently available coarse‐scale SSM datasets (microwave and reanalysis product) and aggregated in situ SSM measurements were implemented to calibrate the time‐invariable coefficients of the SSM retrieval algorithm for conditions in which conventional observations are rare. These coefficients were subsequently used to obtain SSM time series directly from Meteosat Second Generation (MSG) images over the study area of a well‐organized soil moisture network named REMEDHUS in Spain. The results show a high degree of consistency between the estimated and actual SSM time series values when using the three SSM dataset‐calibrated time‐invariable coefficients to retrieve SSM, with coefficients of determination (R²) varying from 0.304 to 0.534 and root mean square errors ranging from 0.020 m³/m³ to 0.029 m³/m³. Further evaluation with different land use types results in acceptable debiased root mean square errors between 0.021 m³/m³ and 0.048 m³/m³ when comparing the estimated MSG pixel‐scale SSM with in situ measurements. These results indicate that the investigated method is practical for deriving time‐invariable coefficients when using publicly accessed coarse‐scale SSM datasets, which is beneficial for generating continuous SSM dataset at the MSG pixel scale.     

Temporal stability of soil moisture in an oasis of northwestern China

Temporal stability of soil moisture has been widely used in hydrological monitoring since it emerged. However, the spatial analysis of temporal stability at the landscape scale is often limited because of insufficient sampling numbers. This work made an effort to investigate the spatial variations of temporal stability of soil moisture in an oasis landscape. The specific objectives of the study were to explore the spatial patterns of temporal stability and to determine the controlling factors of temporal stability in the desert oasis. A time series of soil moisture measurements were gathered on 23 occasions at 118 locations over 3 years in a rectangular transect of approximately 100 km². The nonparametric Spearman's rank correlation coefficient, standard deviation of relative difference (SDRD), and mean absolute bias error (MABE) were used to quantify the temporal stability of soil moisture. Results showed that the temporal stability of soil moisture was depth dependent and season dependent. The spatial pattern of soil moisture in a deep soil layer and between two same seasons generally had a high temporal stability. SDRD and MABE were spatially autocorrelated and exhibited strong spatial structures in the geographic space. The concept of temporal stability can be extended to describe the time‐stable areas of soil moisture with geostatistics. There were great differences between SDRD and MABE in describing the temporal stability of soil moisture and in identifying the controlling factors of temporal stability. In this case, MABE was a better alternative to estimate the areal mean soil moisture using representative locations than SDRD. Land use type, soil moisture condition, and soil particle composition were the dominant controls of temporal stability in the oasis. These insights could help to better understand the essence of temporal stability of soil moisture in arid regions.     

Groundwater recharge in an arid grassland as indicated by soil chloride profile and multiple tracers

Previous studies have shown that shallow groundwater in arid regions is often not in equilibrium with near‐surface boundary conditions due to human activities and climate change. This is especially the case where the unsaturated zone is thick and recharge rate is limited. Under this nonequilibrium condition, the unsaturated zone solute profile plays an important role in estimating recent diffuse recharge in arid environments. This paper combines evaluation of the thick unsaturated zone with the saturated zone to investigate the groundwater recharge of a grassland in the arid western Ordos Basin, NW China, using the soil chloride profiles and multiple tracers (²H, ¹⁸O, ¹³C, ¹⁴C, and water chemistry) of groundwater. Whereas conventional water balance and Darcy flux measurements usually involve large errors in recharge estimations for arid areas, chloride mass balance has been widely and generally successfully used. The results show that the present diffuse recharge beneath the grassland is 0.11–0.32 mm/year, based on the chloride mass balance of seven soil profiles. The chloride accumulation age is approximately 2,500 years at a depth of 13 m in the unsaturated zone. The average Cl content in soil moisture in the upper 13 m of the unsaturated zone ranges from 2,842 to 7,856 mg/L, whereas the shallow groundwater Cl content ranges from 95 to 351 mg/L. The corrected ¹⁴C age of shallow groundwater ranges from 4,327 to 29,708 years. Stable isotopes show that the shallow groundwater is unrelated to modern precipitation. The shallow groundwater was recharged during the cold and wet phases of the Late Pleistocene and Holocene humid phase based on palaeoclimate, and consequently, the groundwater resources are nonrenewable. Due to the limited recharge rate and thick unsaturated zone, the present shallow groundwater has not been in hydraulic equilibrium with near‐surface boundary conditions in the past 2,500 years.     

Taxonomic analysis of the Quaternary archaeofauna found at The Lapa do Santo site,Lagoa Santa region,Brazil

The Lapa do Santo archaeological site, located in the Lagoa Santa region, Brazil, represents an important hunter-gatherer occupation dated from the Early and Middle Holocene. Prior studies of archaeofauna dating to this period are few and most of them only provide basic faunal identification with limited information on taphonomic processes. The main goal of this study is to identify the archaeofauna, record the taphonomic processes, and make inferences about its natural (interpreted as the death of the animal in the rockshelter area due to natural causes or due to predation) or anthropic origins. No extinct species (including megamastofauna) were identified. Taxonomic analysis indicated that most faunal remains consist of taxa that may be of mixed origin (natural or anthropic), such as microvertebrates and carnivorous mammals. The Cervidae are the main vertebrate family found in the studied material and these are likely the result of hunting. The Cervidae Ozotoceros support the presence of a fauna that is typical of savannah environments, although the assemblage contains mostly generalist fauna that can thrive in a variety of environments. Very few remains showed clear evidence for taphonomic processes related to human interaction, such as burnt bones.     

Influence of heterogeneity on unsaturated hydraulic properties: (1) local heterogeneity and scale effect

Spatial heterogeneity is ubiquitous in nature, which may significantly affect the soil hydraulic property curves. The models of a closed‐form functional relationship of soil hydraulic property curves (e.g. VG model or exponential model) are valid at point or local scale based on a point‐scale hydrological process, but how do scale effects of heterogeneity have an influence on the parameters of these models when the models are used in a larger scale process? This paper uses a two‐dimensional variably saturated flow and solute transport finite element model (VSAFT2) to simulate variations of pressure and moisture content in the soil flume under a constant head boundary condition. By changing different numerical simulation block sizes, a quantitative evaluation of parameter variations in the VG model, resulting from the scale effects, is presented. Results show that the parameters of soil hydraulic properties are independent of scale in homogeneous media. Parameters of α and n in homogeneous media, which are estimated by using the unsaturated hydraulic conductivity curve (UHC) or the soil water retention curve (WRC), are identical. Variations of local heterogeneities strongly affect the soil hydraulic properties, and the scale affects the results of the parameter estimations when numerical experiments are conducted. Furthermore, the discrepancy of each curve becomes considerable when moisture content becomes closer to a dry situation. Parameters estimated by UHC are totally different from the ones estimated by WRC. Copyright © 2012 John Wiley & Sons, Ltd.