Modelling the effects of soil type and root distribution on shallow groundwater resources

Vegetated, shallow groundwater environments typically have high environmental and economic value. A sound understanding of the complex interactions and feedbacks between surface vegetation and groundwater resources is crucial to managing and maintaining healthy ecosystems while responding to human needs. A vegetated shallow groundwater environment was modelled using the software HYDRUS 2D to investigate the effects of several combinations of soil type and root distributions on shallow groundwater resources. Three rainfall regimes coupled to both natural and anthropogenically affected groundwater conditions were used to investigate the effect that combinations of four soil types and five root distributions can have on (a) groundwater level drops, (b) groundwater depletion, (c) groundwater recharge and (d) water stress conditions. Vegetation with roots distributed across the whole unsaturated zone and vegetation with dimorphic root systems (i.e. roots having larger concentrations both near the surface and the capillary fringe) behaved differently from vegetation growing roots mainly near the saturated zone. Specifically, vegetation with roots in the unsaturated zone caused water‐table drops and groundwater depletions that were half the amount due to deep‐rooted vegetation. Vegetation with a large portion of roots near the soil surface benefited from rainfall and was less vulnerable to water‐table lowering; as such, the fraction of the total area of roots affected by water stress conditions could be 40% smaller than in the case with deep‐rooted vegetation. However, roots uniformly distributed in the unsaturated zone could halve groundwater recharge rates observed in bare soils. Our analysis provided insights that can enable the formulation of site‐ and purpose‐specific management plans to respond to both human and ecosystem water requirements. Copyright © 2015 John Wiley & Sons, Ltd.     

Evaluating water table response to rainfall events in a shallow aquifer and canal system

Shallow aquifers typically have greater hydrologic connectivity and response to recharge and changes in surface water management practices than deeper aquifers and are therefore often managed to reduce the risk of flooding. Quantification of the water table elevation response under different management scenarios provides valuable information in shallow aquifer systems to assess indirect influences of such modifications. The episodic master recession method was applied to the 15‐min water table elevation and NEXRAD rainfall data for 6 wells to identify water table response and individual rainfall events. The objectives of this study were to evaluate the effects of rainfall, water table elevation, canal stage, site‐specific characteristics, and canal structure modification/water management practice on the fluctuations in water table elevations using multiple/stepwise multiple linear regression techniques. With the modification of canal structure and operation adjustment, significant difference existed in water table response in the southern wells due to its relative downstream position regarding the general groundwater flow direction and the structural modification locations. On average, water table response height and flood risk were lower after than before the structure modification to canals. The effect of rainfall event size on the height of water table response was greater than the effect of antecedent water table elevation and canal stage on the height of water table response. Other factors including leakance of the canal bed sediment, specific yield, and rainfall on i  ? 1 day had significant effects on the height of water table response as well. Antecedent water table elevation and canal stage had greater and more linear effects on the height of water table response after the management changes to canals. Variation in water table response height/rainfall event size ratio was attributed to difference in S _y , antecedent soil water content, hydraulic gradient, rainfall size, and run‐off ratio. After the structure modification, water table response height/rainfall event size ratio demonstrated more linear and proportional relationship with antecedent water table elevation and canal stage.     

Impacts of model initialization on an integrated surface water–groundwater model

Integrated hydrologic models characterize catchment responses by coupling the subsurface flow with land surface processes. One of the major areas of uncertainty in such models is the specification of the initial condition and its influence on subsequent simulations. A key challenge in model initialization is that it requires spatially distributed information on model states, groundwater levels and soil moisture, even when such data are not routinely available. Here, the impact of uncertainty in initial condition was explored across a 208 km² catchment in Denmark using the ParFlow.CLM model. The initialization impact was assessed under two meteorological conditions (wet vs dry) using five depth to water table and soil moisture distributions obtained from various equilibrium states (thermal, root zone, discharge, saturated and unsaturated zone equilibrium) during the model spin‐up. Each of these equilibrium states correspond to varying computation times to achieve stability in a particular aspect of the system state. Results identified particular sensitivity in modelled recharge and stream flow to the different initializations, but reduced sensitivity in modelled energy fluxes. Analysis also suggests that to simulate a year that is wetter than the spin‐up period, an initialization based on discharge equilibrium is adequate to capture the direction and magnitude of surface water–groundwater exchanges. For a drier or hydrologically similar year to the spin‐up period, an initialization based on groundwater equilibrium is required. Variability of monthly subsurface storage changes and discharge bias at the scale of a hydrological event show that the initialization impacts do not diminish as the simulations progress, highlighting the importance of robust and accurate initialization in capturing surface water–groundwater dynamics. Copyright © 2015 John Wiley & Sons, Ltd.     

Water uptake by saltcedar (Tamarix ramosissima) in a desert riparian forest: responses to intra‐annual water table fluctuation

There is considerable interest in naturalizing flow regime on managed rivers to slow the spread of saltcedar (Tamarix ramosissima) invasion in southwestern USA or to preserve riparian forests dominated by saltcedar and other species in northwestern China. However, little is known about the responses of established saltcedar in water sources to frequent intra‐annual fluctuation of water table resulting from this new, more dynamic flow regime. This study investigates how saltcedar at a riparian site in the middle reaches of the Heihe River, northwest China, responds in water sources use to intra‐annual water table fluctuations. Stable oxygen isotope was employed to determine accurate depth at which saltcedar obtains its water supply, and soil moisture monitoring was used to determine sources of plant‐available soil water. We found that the primary zone of water uptake by saltcedar were stable at 25–60 cm depth, but the water sources used by saltcedar switched between groundwater and soil moisture with the water table fluctuations. Saltcedar derived its water from groundwater when water table was at depth less than 60 cm but switched to soil moisture at 25–60 cm depth when water table declined. It is supposed that the well‐developed clay layer at 60–80 cm depth constrained lateral roots of saltcedar to the soil layers above 60 cm, while the fine‐textured soils at this site, which were periodically resaturated by rising groundwater before the stored soil moisture had become depleted, provided an important water reservoir for saltcedar when groundwater dropped below the primary zone of fine roots. The root distribution of saltcedar may also be related to local groundwater history. The quick decline in water table in the early 1980s when the riparian saltcedar had established may strand its roots in the shallow unsaturated zone. We suggested that raising the water table periodically instead of maintaining it invariably above the rooting depth could sustain desired facultative phreatophytes while maximizing water deliveries. Copyright © 2015 John Wiley & Sons, Ltd.     

Water percolation in a thick unsaturated loess layer considering the ground‐atmosphere interaction

The key objective of this paper is to advance our present understanding of how surface water infiltrates in thick unsaturated loess, which is found in arid and semiarid regions of the world, considering the ground‐atmosphere interaction. In situ data for a period of 1 year in thick loess layer at a site in the Loess Plateau of China that has groundwater table at 97.5 m depth were collected for achieving this objective. Climate factors, mainly rainfall and actual evaporation, were measured. In addition, variations of soil temperature and water content at different depths in the unsaturated zone were also measured. The data were used to interpret the water percolation characteristics by dividing the thick unsaturated zone into three zones; namely, (i) surface zone, which constitutes the top 1.0 m, (ii) unsteady zone, which is from 1.0 to 7.0 m, and (iii) steady zone, which is below 7.0 m. In the surface zone, soil temperature and water content are sensitive to climate factors. There is a variation of water content associated with the cumulative influence of infiltration and evaporation in the precipitation and nonprecipitation periods, respectively. In the unsteady zone, the water content is relatively constant; however, temperature varies in different seasons. Water percolation in this zone is both in liquid and vapour phases. In the steady zone, both soil temperature and water content are constant during the entire investigation period. The percolation velocity in this zone is approximately 1.23 × 10^?8 m/s or 0.39 m/year, which suggests that it will take approximately 230.8 years for surface water to pass through the thick unsaturated zone and recharge the groundwater.     

Determining hydrological regimes in an agriculturally used tropical inland valley wetland in Central Uganda using soil moisture,groundwater, and digital elevation data

Inadequate knowledge exists on the distribution of soil moisture and shallow groundwater in intensively cultivated inland valley wetlands in tropical environments, which are required for determining the hydrological regime. This study investigated the spatial and temporal variability of soil moisture along 4 hydrological positions segmented as riparian zone, valley bottom, fringe, and valley slope in an agriculturally used inland valley wetland in Central Uganda. The determined hydrological regimes of the defined hydrological positions are based on soil moisture deficit calculated from the depth to the groundwater table. For that, the accuracy and reliability of satellite‐derived surface models, SRTM‐30m and TanDEM‐X‐12m, for mapping microscale topography and hydrological regimes are evaluated against a 5‐m digital elevation model (DEM) derived from field measurements. Soil moisture and depth to groundwater table were measured using frequency domain reflectometry sensors and piezometers installed along the hydrological positions, respectively. Results showed that spatial and temporal variability in soil moisture increased significantly (p < .05) towards the riparian zone; however, no significant difference was observed between the valley bottom and riparian zone. The distribution of soil hydrological regimes, saturated, near‐saturated, and nonsaturated regimes does not correlate with the hydrological positions. This is due to high spatial and temporal variability in depth to groundwater and soil moisture content across the valley. Precipitation strongly controlled the temporal variability, whereas microscale topography, soil properties, distance from the stream, anthropogenic factors, and land use controlled the spatial variability in the inland valley. TanDEM‐X DEM reasonably mapped the microscale topography and thus soil hydrological regimes relative to the Shuttle Radar Topography Mission DEM. The findings of the study contribute to improved understanding of the distribution of hydrological regimes in an inland valley wetland, which is required for a better agricultural water management planning.     

Refined conceptualization of TOPMODEL for shallow subsurface flows

The TOPMODEL framework was used to derive expressions that account for saturated and unsaturated flow through shallow soil on a hillslope. The resulting equations were the basis for a shallow‐soil TOPMODEL (STOPMODEL). The common TOPMODEL theory implicitly assumes a water table below the entire watershed and this does not conceptually apply to systems hydrologically controlled by shallow interflow of perched groundwater. STOPMODEL provides an approach for extending TOPMODEL's conceptualization to apply to shallow, interflow‐driven watersheds by using soil moisture deficit instead of water table depth as the state variable. Deriving STOPMODEL by using a hydraulic conductivity function that changes exponentially with soil moisture content results in equations that look very similar to those commonly associated with TOPMODEL. This alternative way of conceptualizing TOPMODEL makes the modelling approach available to researchers, planners, and engineers who work in areas where TOPMODEL was previously believed to be unsuited, such as the New York City Watershed in the Catskills region of New York State. Copyright © 2002 John Wiley & Sons, Ltd.     

Nitrate attenuation in soil and shallow groundwater under a bottomland hedgerow in a European farming landscape

In many agricultural areas, hedgerows give rise to strong expectations of reducing the inputs of excess nitrate to the groundwater and rivers. This study aims to analyse the spatial and seasonal influences of a hedgerow on nitrate dynamics in the soil and groundwater. Nitrate (NO₃^?) and chloride (Cl^?) concentrations were measured with spatially dense sampling in the unsaturated soil and in the groundwater along a transect intersecting a bottomland oak (Quercus rubor) hedgerow after the growing season and during the dormant season. We explain NO₃^? dynamics by using Cl^? as an index of tree‐root extension and water transfer. At the end of the growing season, NO₃^? is entirely absorbed by the trees over a large and deep volume corresponding to the rooting zone, where, in contrast Cl^? is highly concentrated due to root exclusion. However, these observed patterns in the soil have no influence on the deep groundwater composition at this season. During the dormant season, water transfer processes feeding the shallow groundwater layer are different upslope and downslope from the hedgerow in relation to the thickness of the unsaturated zone. Upslope, the shallow groundwater is fed by rainwater infiltration through the soil which favours Cl^? dilution. Right under the hedge and downslope, the rapid ascent of the groundwater near the ground surface prevents rainwater input and Cl^? dilution. Under the hedgerow the highest concentrations of Cl^? coincide with the absence of NO₃^? in the shallow groundwater layer and with high concentrations of dissolved organic carbon. The absence of NO₃^? during the dormant season seems to be due to denitrification in the hedgerow rooting zone when it is rapidly saturated by groundwater. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling hydrological management for the restoration of acidified floating fens

Wetlands show a large decline in biodiversity. To protect and restore this biodiversity, many restoration projects are carried out. Hydrology in wetlands controls the chemical and biological processes and may be the most important factor regulating wetland function and development. Hydrological models may be used to simulate these processes and to evaluate management scenarios for restoration. HYDRUS2D, a combined saturated–unsaturated groundwater flow and transport model, is presented. This simulates near‐surface hydrological processes in an acidified floating fen, with the aim to evaluate the effect of hydrological restoration in terms of conditions for biodiversity. In the acidified floating fen in the nature reserve Ilperveld (The Netherlands), a trench system was dug for the purpose of creating a runoff channel for acid rainwater in wet periods and to enable circum‐neutral surface water to enter the fen in dry periods. The model is calibrated against measured conductivity values for a 5 year period. From the model simulations, it was found that lateral flow in the floating raft is limited. Furthermore, the model shows that the best management option is a combination of trenches and inundation, which gave the best soil water quality in the root zone. It is concluded that hydrological models can be used for the calculation of management scenarios in restoration projects. The combined saturated–unsaturated model concept used in this paper is able to incorporate the governing hydrological processes in the wetland root zones. Copyright © 2005 John Wiley & Sons, Ltd.     

A large weighing lysimeter for evapotranspiration and soil‐water–groundwater exchange studies

A large weighing lysimeter was installed at Yucheng Comprehensive Experimental Station, north China, for evapotranspiration and soil‐water–groundwater exchange studies. Features of the lysimeter include the following: (i) mass resolution equivalent to 0·016 mm of water to accurately and simultaneously determine hourly evapotranspiration, surface evaporation and groundwater recharge; (ii) a surface area of 3·14 m² and a soil profile depth of 5·0 m to permit normal plant development, soil‐water extraction, soil‐water–groundwater exchanges, and fluctuations of groundwater level; (iii) a special supply–drainage system to simulate field conditions of groundwater within the lysimeter; (iv) a soil mass of about 30 Mg, including both unsaturated and saturated loam. The soil consists mainly of mealy sand and light loam. Monitoring the vegetated lysimeter during the growing period of winter wheat, from October 1998 through to June 1999, indicated that during the period groundwater evaporation contributed 16·6% of total evapotranspiration for a water‐table depth from 1·6 m to 2·4 m below ground surface. Too much irrigation reduced the amount of upward water flow from the groundwater table, and caused deep percolation to the groundwater. Data from neutron probe and tensiometers suggest that soil‐water‐content profiles and soil‐water‐potential profiles were strongly affected by shallow groundwater. Copyright © 2000 John Wiley & Sons, Ltd.     

Experimental study and mathematical modelling of soluble chemical transfer from unsaturated/saturated soil to surface runoff

A one‐dimensional, two‐layer solute transport model is developed to simulate chemical transport process in an initially unsaturated soil with ponding water on the soil surface before surface runoff starts. The developed mathematical model is tested against a laboratory experiment. The infiltration and diffusion processes are mathematically lumped together and described by incomplete mixing parameters. Based on mass conservation and water balance equations, the model is developed to describe solute transport in a two‐zone layer, a ponding runoff zone and a soil mixing zone. The two‐zone layer is treated as one system to avoid describing the complicated chemical transport processes near the soil surface in the mixing zone. The proposed model was analytically solved, and the solutions agreed well with the experimental data. The developed experimental method and mathematical model were used to study the effect of the soil initial moisture saturation on chemical concentration in surface runoff. The study results indicated that, when the soil was initially saturated, chemical concentration in surface runoff was significantly (two orders of magnitude) higher than that with initially unsaturated soil, while the initial chemical concentrations at the two cases were of the same magnitude. The soil mixing depth for the initially unsaturated soil was much larger than that for the initially saturated soil, and the incomplete runoff mixing parameter was larger for the initially unsaturated soil. The higher the infiltration rate of the soil, the greater the infiltration‐related incomplete mixing parameter. According to the quantitative analysis, the soil mixing depth was found to be sensitive for both initially unsaturated and saturated soils, and the incomplete runoff mixing parameter was sensitive for initially saturated soil but not for the initially unsaturated soil; the incomplete infiltration mixing parameter behaved just the opposite. Some suggestions are made for reducing chemical loss from runoff. Copyright © 2010 John Wiley & Sons, Ltd.     

Multi‐scale assimilation of root zone soil water predictions

When hydrology model parameters are determined, a traditional data assimilation method (such as Kalman filter) and a hydrology model can estimate the root zone soil water with uncertain state variables (such as initial soil water content). The simulated result can be quite good. However, when a key soil hydraulic property, such as the saturated hydraulic conductivity, is overestimated or underestimated, the traditional soil water assimilation process will produce a persistent bias in its predictions. In this paper, we present and demonstrate a new multi‐scale assimilation method by combining the direct insertion assimilation method, particle swarm optimisation (PSO) algorithm and Richards equation. We study the possibility of estimating root zone soil water with a multi‐scale assimilation method by using observed in situ data from the Wudaogou experiment station, Huaihe River Basin, China. The results indicate there is a persistent bias between simulated and observed values when the direct insertion assimilation surface soil water content is used to estimate root zone soil water contents. Using a multi‐scale assimilation method (PSO algorithm and direct insertion assimilation) and an assumed bottom boundary condition, the results show some obvious improvement, but the root mean square error is still relatively large. When the bottom boundary condition is similar to the actual situation, the multi‐scale assimilation method can well represent the root zone soil water content. The results indicate that the method is useful in estimating root zone soil water when available soil water data are limited to the surface layer and the initial soil water content even when the soil hydraulic conductivities are uncertain. Copyright © 2011 John Wiley & Sons, Ltd.     

Evidence for climatic and hillslope‐aspect controls on vadose zone hydrology and implications for saprolite weathering

Through the delivery of water in snowmelt, climate should govern the rate and extent of saprolite formation in snow‐dominated mountain watersheds, yet the mechanisms by which water flows deeply into regolith are largely unexplored. In this study we link rainfall, snow depth, and water content data from both soil and shallow saprolite to document vadose zone dynamics in two montane catchments over 2 years. Measurements of snow pack thickness and soil moisture reveal strong contrasts between north‐ and south‐facing slopes in both the timing of meltwater delivery and the duration of significant soil wetting in the shallow vadose zone. Despite similar magnitudes of snowmelt recharge, north‐facing slopes have higher sustained soil moisture compared to south‐facing slopes. To help interpret these observations, we use a 2D numerical model of vadose zone dynamics to calculate the expected space–time moisture patterns on an idealized hillslope under two wetting scenarios: a single sustained recharge pulse versus a set of short pulses. The model predicts that the duration of the recharge event exerts a stronger control on the depth and residence time of water in the upper unsaturated zone than the magnitude of the recharge event. Model calculations also imply that water should move more slowly through the subsurface and downward water flux should be substantially reduced when water is applied in several pulses rather than in one sustained event. The results suggest that thicker soil and more deeply weathered rock on north‐facing slopes may reflect greater water supply to the deep subsurface. Copyright © 2015 John Wiley & Sons, Ltd.     

How hydrological factors initiate instability in a model sandy slope

Knowledge of the mechanisms of rain‐induced shallow landslides can improve the prediction of their occurrence and mitigate subsequent sediment disasters. Here, we examine an artificial slope's subsurface hydrology and propose a new slope stability analysis that includes seepage force and the down‐slope transfer of excess shear forces. We measured pore water pressure and volumetric water content immediately prior to a shallow landslide on an artificial sandy slope of 32°: The direction of the subsurface flow shifted from downward to parallel to the slope in the deepest part of the landslide mass, and this shift coincided with the start of soil displacement. A slope stability analysis that was restricted to individual segments of the landslide mass could not explain the initiation of the landslide; however, inclusion of the transfer of excess shear forces from up‐slope to down‐slope segments improved drastically the predictability. The improved stability analysis revealed that an unstable zone expanded down‐slope with an increase in soil water content, showing that the down‐slope soil initially supported the unstable up‐slope soil; destabilization of this down‐slope soil was the eventual trigger of total slope collapse. Initially, the effect of apparent soil cohesion was the most important factor promoting slope stability, but seepage force became the most important factor promoting slope instability closer to the landslide occurrence. These findings indicate that seepage forces, controlled by changes in direction and magnitude of saturated and unsaturated subsurface flows, may be the main cause of shallow landslides in sandy slopes. Copyright © 2013 John Wiley & Sons, Ltd.     

Changes in water table level influence solute transport in uniform porous media

Changes in the water table level result in variable water saturation and variable hydrological fluxes at the interface between the unsaturated and saturated zone. This may influence the transport and fate of contaminants in the subsurface. The objective of this study was to examine the impact of a decreasing and an increasing water table on solute transport. We conducted tracer experiments at downward flow conditions in laboratory columns filled with two different uniform porous media under static and transient flow conditions either increasing or decreasing the water table. Tracer breakthrough curves were simulated using a mobile–immobile transport model. The resulting transport parameters were compared to identify dominant transport processes. Changes in the water table level affected dispersivities and mobile water fractions depending on the direction of water table movement and the grain size of the porous media. In fine glass beads, the water flow velocity was similar to the decline rate of the water table, and the mobile water fraction was decreased compared with steady‐state saturated conditions. However, immobile water was negligible. In coarse glass beads, water flow was faster because of fingered flow in the unsaturated part, and the mobile water fraction was smaller than in the fine material. Here, a rising water table led to an even smaller mobile water fraction and increased solute spreading because of diffusive interaction with immobile water. We conclude that changes of the water table need to be considered to correctly simulate transport in the subsurface at the transition of the unsaturated–saturated zone. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Simulating soil‐water movement under a hedgerow surrounding a bottomland reveals the importance of transpiration in water balance

The objective of this study was to quantify components of the water balance related to root‐water uptake in the soil below a hedgerow. At this local scale, a two‐dimensional (2D) flow domain in the x–z plane 6 m long and 1·55 m deep was considered. An attempt was made to estimate transpiration using a simulation model. The SWMS‐2D model was modified and used to simulate temporally and spatially heterogeneous boundary conditions. A function with a variable spatial distribution of root‐water uptake was considered, and model calibration was performed by adjusting this root‐water uptake distribution. Observed data from a previous field study were compared against model predictions. During the validation step, satisfactory agreement was obtained, as the difference between observed and modelled pressure head values was less than 50 cm for 80% of the study data. Hedge transpiration capacity is a significant component of soil‐water balance in the summer, when predicted transpiration reaches about 5·6 mm day^?1. One of the most important findings is that hedge transpiration is nearly twice that of a forest canopy. In addition, soil‐water content is significantly different whether downslope or upslope depending on the root‐water uptake. The high transpiration rate was mainly due to the presence of a shallow water table below the hedgerow trees. Soil‐water content was not a limiting factor for transpiration in this context, as it could be in one with a much deeper water table. Hedgerow tree transpiration exerts a strong impact not only on water content within the vadose zone but also on the water‐table profile along the transect. Results obtained at the local scale reveal that the global impact of hedges at the catchment scale has been underestimated in the past. Transpiration rate exerts a major influence on water balance at both the seasonal and annual scales for watersheds with a dense network of hedgerows. Copyright © 2007 John Wiley & Sons, Ltd.     

Water content and porosity estimated from ground-penetrating radar and resistivity

Both ground-penetrating radar and the resistivity method have proven to be useful tools for exploring water content variations, since related parameters such as dielectric constant and the resistivity of rocks and sediments are highly dependent on moisture content. These methods were used independently to estimate volumetric water content in the unsaturated zone and porosity in the saturated zone in a 100-m sandy section. Two sample sites along the profile were also chosen for a shallow geophysical investigation and soil sampling, to enable the calibration and verification of the indirect geophysical methods. The grain distribution at these sites is dominated by medium-sized sand (0.25–0.5 mm). The water content was 6.9 vol.% and calculated porosities are 37% and 40% respectively. At each of these sites the mean water content values calculated from resistivity are within one percentage unit of measured water content while those calculated from ground-penetrating radar give higher values by as much as 2.9 percentage units. The water contents in the unsaturated zone in the section, estimated from resistivity and ground-penetrating radar, show very similar trends, although that deduced from ground-penetrating radar is generally somewhat larger, consistent with the results from the sample sites. The mean porosity values obtained from the two methods in the saturated zone are in good agreement.     

Influence of seepage undercutting on the stability of root‐reinforced streambanks

Several mechanisms contribute to streambank failure including fluvial toe undercutting, reduced soil shear strength by increased soil pore‐water pressure, and seepage erosion. Recent research has suggested that seepage erosion of noncohesive soil layers undercutting the banks may play an equivalent role in streambank failure to increased soil pore‐water pressure. However, this past research has primarily been limited to laboratory studies of non‐vegetated banks. The objective of this research was to utilize the Bank Stability and Toe Erosion Model (BSTEM) in order to determine the importance of seepage undercutting relative to bank shear strength, bank angle, soil pore‐water pressure, and root reinforcement. The BSTEM simulated two streambanks: Little Topashaw Creek and Goodwin Creek in northern Mississippi. Simulations included three bank angles (70° to 90°), four pore‐water pressure distributions (unsaturated, two partially saturated cases, and fully saturated), six distances of undercutting (0 to 40 cm), and 13 different vegetation conditions (root cohesions from 0·0 to 15·0 kPa). A relative sensitivity analysis suggested that BSTEM was approximately three to four times more sensitive to water table position than root cohesion or depth of seepage undercutting. Seepage undercutting becomes a prominent bank failure mechanism on unsaturated to partially saturated streambanks with root reinforcement, even with undercutting distances as small as 20 cm. Consideration of seepage undercutting is less important under conditions of partially to fully saturated soil pore‐water conditions. The distance at which instability by undercutting became equivalent to instability by increased soil pore‐water pressure decreased as root reinforcement increased, with values typically ranging between 20 and 40 cm at Little Topashaw Creek and between 20 and 55 cm at Goodwin Creek. This research depicts the baseline conditions at which seepage undercutting of vegetated streambanks needs to be considered for bank stability analyses. Copyright © 2008 John Wiley & Sons, Ltd.     

A MODFLOW Infiltration Device Package for Simulating Storm Water Infiltration

This article describes a MODFLOW Infiltration Device (INFD) Package that can simulate infiltration devices and their two‐way interaction with groundwater. The INFD Package relies on a water balance including inflow of storm water, leakage‐like seepage through the device faces, overflow, and change in storage. The water balance for the device can be simulated in multiple INFD time steps within a single MODFLOW time step, and infiltration from the device can be routed through the unsaturated zone to the groundwater table. A benchmark test shows that the INFD Package's analytical solution for stage computes exact results for transient behavior. To achieve similar accuracy by the numerical solution of the MODFLOW Surface‐Water Routing (SWR1) Process requires many small time steps. Furthermore, the INFD Package includes an improved representation of flow through the INFD sides that results in lower infiltration rates than simulated by SWR1. The INFD Package is also demonstrated in a transient simulation of a hypothetical catchment where two devices interact differently with groundwater. This simulation demonstrates that device and groundwater interaction depends on the thickness of the unsaturated zone because a shallow groundwater table (a likely result from storm water infiltration itself) may occupy retention volume, whereas a thick unsaturated zone may cause a phase shift and a change of amplitude in groundwater table response to a change of infiltration. We thus find that the INFD Package accommodates the simulation of infiltration devices and groundwater in an integrated manner on small as well as large spatial and temporal scales.