共查询到20条相似文献，搜索用时 359 毫秒
1.
Previous experimental studies of capillary barriers have identified highly hysteretic soil moisture retention characteristics in the materials used. In this study, numerical modelling is used to analyse the role of soil moisture hysteresis in capillary barrier functioning. Comparisons between published experimental results and model simulations indicate that soil moisture hysteresis was a necessary inclusion in the modelling approach to adequately reproduce pore water pressure distributions and the timing of breakthrough occurrences. Under hypothetical intermittent infiltration and evaporation conditions, the predicted volumetric water content in the moisture retention layer was significantly different for hysteretic and non‐hysteretic models. The hysteresis effect was found to be dependent on the nature of infiltration–evaporation cycling, although the predicted volume of flow through the hysteretic barrier was lower than that of the non‐hysteretic case, regardless of the nature of the cyclic upper boundary conditions. For practical engineering designs, where the water leakage through the barrier is the primary concern, the inclusion of soil moisture hysteresis in numerical modelling is needed. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
2.
The one-dimensional pressure head profile above a fixed water-table was studied for different steady infiltration rates. As shown in previous studies, when the infiltration rate (qin) is less than the soil saturated hydraulic conductivity (Ks), this profile can be divided into two parts: (1) from the water-table surface (z0) to an elevation zγ, the pressure head varies from 0 to a value hγ such as K(hγ)=−qin; (2) above the elevation zγ the pressure head is constant and equal to hγ. Above the water-table the zone where the pressure head is variable has been called ‘transition zone’. Its height is shown to be highly variable and to depend on soil properties as well as on the infiltration rate. This transition zone is not the ‘capillary fringe’ as defined by Gillham (Gillham R.W., 1984. The capillary fringe and its effect on water-table response. J. Hydrol. 67, 307–324). Numerical experiments performed with HYDRUS-2D® for the case of a water-table drained by parallel ditches have shown that the height of the transition zone is similar in the one-dimensional profile and in the two-dimensional system as long as the local slope of the water-table is small. This result is important since in a two-dimensional system, the transition zone is the place where all the horizontal component of the unsaturated flow occurs. The ratio of the horizontal component of the unsaturated flow vs the total horizontal flow in both the unsaturated and saturated zones has been computed. For a given soil, this ratio decreases as the infiltration rate increases; for a given infiltration rate, the soil with the thinnest transition zone transfers the largest amount of water above the water-table. 相似文献
3.
Non-unique solutions of inverse problems arise from a lack of information that satisfies necessary conditions for the problem to be well defined. This paper investigates these conditions for inverse modeling of water flow through multi-dimensional variably saturated porous media. It shows that in order to obtain a unique estimate of hydraulic parameters, along each streamline of the flow field (1) spatial and temporal head observations must be given; (2) the number of spatial and temporal head observations required should be greater or equal to the number of unknown parameters; (3) the flux boundary condition or the pumping rate of a well must be specified for the homogeneous case and both boundary flux and pumping rate are a must for the heterogeneous case; (4) head observations must encompass both saturated and unsaturated conditions, and the functional relationships for unsaturated hydraulic conductivity/pressure head and for the moisture retention should be given, and (5) the residual water content value also need to be specified a priori or water content measurements are needed for the estimation of the saturated water content.For field problems, these necessary conditions can be collected or estimated but likely involve uncertainty. While the problems become well defined and have unique solutions, the solutions likely will be uncertain. Because of this uncertainty, stochastic approaches are deemed to be appropriate for inverse problems as they are for forward problems to address uncertainty. Nevertheless, knowledge of these necessary conditions is critical to reduce uncertainty in both characterization of the vadose zone and the aquifer, and prediction of water flow and solute migration in the subsurface. 相似文献
4.
This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results. 相似文献
5.
《Advances in water resources》2007,30(6-7):1593-1607
An experimental and numerical investigation was conducted to study the colonization dynamics of a bioluminescent bacterium, Pseudomonas fluorescens HK44, during growth in a porous medium under steady, variably saturated flow conditions. Experiments were conducted in a thin-slab light transmission chamber filled with uniform, translucent quartz sand. Steady, variably saturated flow conditions were established using drip emitters mounted on the top of the chamber, with glucose applied through a central dripper located directly above an inoculated region of the chamber. Periodic pulses of salicylate and a dye tracer were applied to induce bioluminescence of the bacterium to monitor colony expansion and to track changes in the hydraulic and transport properties of the sand. Changes in the apparent water saturation of the sand were quantified by monitoring light transmission through the chamber with a CCD camera. The colonized region expanded laterally by about 15 cm, and upward against the flow by 7–8 cm during the 6-day experiment while apparent saturations in the colonized region decreased by 7–9% and the capillary fringe dropped by ∼5 cm. The observed data were reproduced approximately using a numerical model that accounted for the processes of water flow, solute and bacterial transport, cell growth and accumulation, glucose and oxygen consumption, and gas diffusion and exchange. The results of this study illustrate some of the complexities associated with coupled flow, reactive transport, and biological processes in variably saturated porous media, such as localized desaturation, capillary fringe lowering effects, and upstream movement of bacterial colonization, that may not readily observable using other experimental techniques. 相似文献
6.
This paper describes a stochastic analysis of steady state flow in a bounded, partially saturated heterogeneous porous medium subject to distributed infiltration. The presence of boundary conditions leads to non-uniformity in the mean unsaturated flow, which in turn causes non-stationarity in the statistics of velocity fields. Motivated by this, our aim is to investigate the impact of boundary conditions on the behavior of field-scale unsaturated flow. Within the framework of spectral theory based on Fourier–Stieltjes representations for the perturbed quantities, the general expressions for the pressure head variance, variance of log unsaturated hydraulic conductivity and variance of the specific discharge are presented in the wave number domain. Closed-form expressions are developed for the simplified case of statistical isotropy of the log hydraulic conductivity field with a constant soil pore-size distribution parameter. These expressions allow us to investigate the impact of the boundary conditions, namely the vertical infiltration from the soil surface and a prescribed pressure head at a certain depth below the soil surface. It is found that the boundary conditions are critical in predicting uncertainty in bounded unsaturated flow. Our analytical expression for the pressure head variance in a one-dimensional, heterogeneous flow domain, developed using a nonstationary spectral representation approach [Li S-G, McLaughlin D. A nonstationary spectral method for solving stochastic groundwater problems: unconditional analysis. Water Resour Res 1991;27(7):1589–605; Li S-G, McLaughlin D. Using the nonstationary spectral method to analyze flow through heterogeneous trending media. Water Resour Res 1995; 31(3):541–51], is precisely equivalent to the published result of Lu et al. [Lu Z, Zhang D. Analytical solutions to steady state unsaturated flow in layered, randomly heterogeneous soils via Kirchhoff transformation. Adv Water Resour 2004;27:775–84]. 相似文献
7.
Role of epikarst in near‐surface hydrological processes in a soil mantled subtropical dolomite karst slope: implications of field rainfall simulation experiments 
下载免费PDF全文
下载免费PDF全文 Epikarst exerts a strong control on run‐off generation in karst regions, but it is still unclear in karst regions. Our study aimed to demonstrate the effect of epikarst on near‐surface hydrological processes in a subtropical cockpit karst region of southwest China, using plot‐scale rainfall simulation experiments with different rainfall intensities (low and high) and antecedent moisture conditions (dry and wet). A trench excavated to the epikarst lower boundary allowed identification of flow pathways in the entire soil–epikarst architecture system, thus facilitating the water balance calculations using a conceptual model with the assumption of a two‐stage hydrological evolution. More than 70% of the total rainfall water moved vertically through the shallow soil layer and then was redistributed by the epikarst as subsurface flow occurring on the soil–epikarst interface, depression filling on epikarst surface, water held by epikarst and deep percolation. Epikarst water regulation capacity, defined as the sum of depression filling on epikarst surface, water held by epikarst, epikarst seepage flow and deep percolation, was 58 mm (wet antecedent condition) and 223 mm (dry antecedent condition). Total run‐off from the soil–epikarst system was dominated by saturated subsurface flow showing a threshold process controlled by epikarst storage capacity (storing as much as 181 mm of rainfall water under dry antecedent condition). Our study proved that despite the epikarst being relatively poorly developed and covered by a soil mantle, it still exerted a strong influence on near‐surface hydrological processes and thus should be adequately considered in future modelling of water recharge and depletion dynamics in this integrated soil–epikarst system. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
8.
Stormwater infiltration systems are a popular method for urban stormwater control. They are often designed using an assumption of one‐dimensional saturated outflow, although this is not very accurate for many typical designs where two‐dimensional (2D) flows into unsaturated soils occur. Available 2D variably saturated flow models are not commonly used for design because of their complexity and difficulties with the required boundary conditions. A purpose‐built stormwater infiltration system model was thus developed for the simulation of 2D flow from a porous storage. The model combines a soil moisture–based model for unsaturated soils with a ponded storage model and uses a wetting front‐tracking approach for saturated flows. The model represents the main physical processes while minimizing input data requirements. The model was calibrated and validated using data from laboratory 2D stormwater infiltration trench experiments. Calibrations were undertaken using five different combinations of calibration data to examine calibration data requirements. It was found that storage water levels could be satisfactorily predicted using parameters calibrated with either data from laboratory soils tests or observed water level data, whereas the prediction of soil moistures was improved through the addition of observed soil moisture data to the calibration data set. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
9.
The errors in surface runoff prediction by neglecting the relationship between infiltration rate and overland flow depth 总被引：9，自引：0，他引：9
Many simplifications are used in modeling surface runoff over a uniform slope. A very common simplification is to determine the infiltration rate independent of the overland flow depth and to combine it afterward with the kinematic-wave equation to determine the overland flow depth. Another simplication is to replace the spatially variable infiltration rates along the slope i(x, t) due to the water depth variations h(x,t) with an infiltration rate that is determined at a certain location along the slope. The aim of this study is to evaluate the errors induced by these simplications on predicted infiltration rates, overland flow depths, and total runoff volume. The error analysis is accomplished by comparing a simplified model with a model where the interaction between the overland flow depth and infiltration rate is counted. In this model, the infiltration rate is assumed to vary along the slope with the overland flow depth, even for homogeneous soil profiles. The kinematic-wave equation with interactive infiltration rate, calculated along the slopy by Richard's equation, are then solved by a finite difference scheme for a 100-m-long uniform slope. In the first error analysis, we study the effect of combining an ‘exact’ and ‘approximate’ one-dimensional infiltration rate with the kinematic-wave equation for three different soil surface roughness coefficients. The terms ‘exact’ and ‘approximate’ stand for the solution of Richard's equation with and without using the overland flow depth in the boundary condition, respectively. The simulations showed that higher infiltration rates and lower overland flow depths are obtained during the rising stage of the hydrograph when overland flow depth is used in the upper boundary condition of the one-dimensional Richard's equation. During the recession period, the simplified model predicts lower infiltration rates and higher overland flow depths. The absolute relative errors between the ‘exact’ and ‘approximate’ solutions are positively correlated to the overland flow depths which increase with the soil surface roughness coefficient. For this error analysis, the relative errors in surface runoff volume per unit slope width throughout the storm are much smaller than the relative errors in momentary overland flow depths and discharges due to the alternate signs of the deviations along the rising and falling stages. In the second error analysis, when the spatially variable infiltration rate along the slope i(x, t) is replaced in the kinematic-wave equation by i(t), calculated at the slope outlet, the overland flow depth is underestimated during the rising stage of the hydrograph and overestimated during the falling stage. The deviations during the rising stage are much smaller than the deviations during the falling stage, but they are of a longer duration. This occurs because the solution with i(x, t) recognizes that part of the slope becomes dry after rainfall stops, while overland flow still exists with i(t) determined at the slope outlet. As obtained for the first error analysis, the relative errors in surface runoff volume per unit slope width are also much smaller than the relative errors in momentary overland flow depths and discharges. The relation between the errors in overland flow depth and discharge to different mathematical simplifications enables to evaluate whether certain simplifications are justified or more computational efforts should be used. 相似文献
10.
Macropores are a relatively small proportion of the soil volume, but they play an important role in the movement of water and chemicals owing to occasional rapid fluxes through them. The occurrence of macropore flow does not depend on the water content (or potential) of the bluk matrix unless the soil is close to saturation, but depends instead principally upon surface boundary conditions. Accordingly, three control situations of infiltration are recognized: macropore control, application control, and matrix control. These three situations indicate that the two-domain system may be a proper approach for the simulation of macropore soil. In this conceptualization, macropores are defined as channeling pores of different radii in which the flux density (with unit hydraulic gradient) occurring in the minimum sizes of such pores is greater than or equal to the saturated matrix hydraulic conductivity. Recognizing the two structural domains of the macropore and matrix, and possible water flow situations, three flow regions are suggested: matrix, macropore, and transaction. The matrix and the macropore are the two domains, and the transaction represents the exchange of water between the matrix and the macropore. The classic approach of the Richards equation is applicable to describe the flow in the matrix domain. The Hagen-Poiseuille and the Chezy-Manning equations for tube flow can be applied to represent the relationship between the hydraulic conductivity of the macroporosity and the total macroporosity, where the total macroporosity is defined as the ratio of the summed macropore cross-sectional area and the total soil cross-sectional area. An equation describing water flow in the macropore domain is then obtained. 相似文献
11.
Inflow to a tunnel is a great public concern and is closely related to groundwater hydrology, geotechnical engineering, and mining engineering, among other disciplines. Rapid computation of inflow to a tunnel provides a timely means for quickly assessing the inflow discharge, thus is critical for safe operation of tunnels. Dewatering of tunnels is another engineering practice that should be planned. In this study, an analytical solution of the inflow to a tunnel in a fractured unconfined aquifer is obtained. The solution takes into account either the spherical or slab-shaped matrix block and the unsteady state interporosity flow. The instantaneous drainage water table and anisotropic hydraulic conductivities of the fractures network are also considered. Both uniform flux and uniform head boundary condition are considered to simulate the constant head boundary condition in the tunnel. The effects of the hydraulic parameters of the fractured aquifer on the inflow variation of the tunnel are explored. The application of the presented solution to obtain the optimum location and discharge of the well to minimize the inflow to a tunnel is illustrated. 相似文献
12.
An exact, closed-form analytical solution is developed for calculating ground water transit times within Dupuit-type flow systems. The solution applies to steady-state, saturated flow through an unconfined, horizontal aquifer recharged by surface infiltration and discharging to a downgradient fixed-head boundary. The upgradient boundary can represent, using the same equation, a no-flow boundary or a fixed head. The approach is unique for calculating travel times because it makes no a priori assumptions regarding the limit of the water table rise with respect to the minimum saturated aquifer thickness. The computed travel times are verified against a numerical model, and examples are provided, which show that the predicted travel times can be on the order of nine times longer relative to existing analytical solutions. 相似文献
13.
14.
Subsurface drainage systems have been widely used to deal with soil salinization and waterlogging problems around the world. In this paper, a mathematical model was introduced to quantify the transient behavior of the groundwater table and the seepage from a subsurface drainage system. Based on the assumption of a hydrostatic pressure distribution, the model considered the pore-water flow in both the phreatic and vadose soil zones. An approximate analytical solution for the model was derived to quantify the drainage of soils which were initially water-saturated. The analytical solution was validated against laboratory experiments and a 2-D Richards equation-based model, and found to predict well the transient water seepage from the subsurface drainage system. A saturated flow-based model was also tested and found to over-predict the time required for drainage and the total water seepage by nearly one order of magnitude, in comparison with the experimental results and the present analytical solution. During drainage, a vadose zone with a significant water storage capacity developed above the phreatic surface. A considerable amount of water still remained in the vadose zone at the steady state with the water table situated at the drain bottom. Sensitivity analyses demonstrated that effects of the vadose zone were intensified with an increased thickness of capillary fringe, capillary rise and/or burying depth of drains, in terms of the required drainage time and total water seepage. The analytical solution provides guidance for assessing the capillary effects on the effectiveness and efficiency of subsurface drainage systems for combating soil salinization and waterlogging problems. 相似文献
15.
Effects of rainfall on soil moisture and water movement in a subalpine dark coniferous forest in southwestern China 总被引：3，自引：0，他引：3
Qing Xu Shirong Liu Xianchong Wan Chunqian Jiang Xianfang Song Jingxin Wang 《水文研究》2012,26(25):3800-3809
Water content and movement in soil profile and hydrogen isotope composition (δD) of soil water, rainwater, and groundwater were examined in a subalpine dark coniferous forest in the Wolong National Nature Reserve in Sichuan, China, following rainfall events in 2003–2004. Light rainfall increased water content in the litter and at soil depth of 0–80 cm, but the increased soil water was lost in several days. Heavy rainfall increased soil water content up to 85% at depths of 0–40 cm. Following the light rainfall in early spring, the δD of water from the litter, humus, illuvial, and material layers decreased first and then gradually reached the pre‐rainfall level. In summer, light rainfall reached the litter humus, and illuvial layer, but did not hit the material layer. Heavy rainfall affected δD of water in all layers. The δD of soil interflow slightly fluctuated with rainfall events. The δD of shallow groundwater did not differ significantly among all rainfall events. Light rainfall altered the shape of δD profile curve of water in the upper layer of soil, whereas heavy rainfall greatly affected the shape of δD profile curve of water in all soil layers. Following the heavy rainfall, preferential flow initially occurred through macropores, decayed plant roots, and rocks at different depths of soil profile. With continuing rainfall, the litter and surface soil were nearly saturated or fully saturated, and infiltration became homogeneous and plug‐like. Forest soil water, particularly in deeper soil profile, was slightly affected by rainfall and, thus, can be a source of water supply for regional needs, particularly during dry seasons. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
16.
The dynamic response of a tunnel buried in a two-dimensional poroelastic soil layer subjected to a moving point load was investigated theoretically. The tunnel was simplified as an infinite long Euler–Bernoulli beam, which was placed parallel to the traction-free ground surface. The saturated layer was governed by Biot’s theory. Combined with the specified boundary conditions along the beam and saturated poroelastic layer, the coupled equations of the system were solved analytically in the frequency–wavenumber domain based on Fourier transform. The time domain responses were obtained by the fast inverse Fourier transform. The critical velocity of the considered structure was determined from the dispersion curves. The different dynamic characteristics of the elastic soil medium and the saturated poroelastic medium subjected to the underground moving load were investigated. It is concluded that, for coarse materials or fine materials subjected to the high-velocity loading, models ignoring the coupling effects between the pore fluid and the soil skeleton may cause errors. The shear modulus and the permeability coefficients of the saturated soil as well as the load moving velocity had significant influence on the displacement and pore pressure responses. 相似文献
17.
Characterising deep vadose zone water movement and solute transport under typical irrigated cropland in the North China Plain 下载免费PDF全文
下载免费PDF全文 Understanding the dynamics and mechanisms of soil water movement and solute transport is essential for accurately estimating recharge rates and evaluating the impacts of agricultural activities on groundwater resources. In a thick vadose zone (0–15 m) under irrigated cropland in the piedmont region of the North China Plain, soil water content, matric potential, and solute concentrations were measured. Based on these data, the dynamics of soil water and solutes were analysed to investigate the mechanisms of soil water and solute transport. The study showed that the 0–15‐m vadose zone can be divided into three layers: an infiltration and evaporation layer (0–2 m), an unsteady infiltration layer (2–6 m), and a quasi‐steady infiltration layer (6–15 m). The chloride, nitrate, and sulphate concentrations all showed greater variations in the upper soil layer (0–1 m) compared to values in the deep vadose zone (below 2 m). The average concentrations of these three anions in the deep vadose zone varied insignificantly with depth and approached values of 125, 242, and 116 mg/L. The accumulated chloride, sulphate, and nitrate were 2,179 ± 113, 1,760 ± 383, and 4,074 ± 421 kg/ha, respectively. The soil water potential and solute concentrations indicated that uniform flow and preferential flow both occurred in the deep vadose zone, and uniform flow was the dominant mechanism of soil water movement in this study. The piston‐like flow velocity of solute transport was 1.14 m per year, and the average value of calculated leached nitrate nitrogen was 107 kg/ha?year below the root zone. The results can be used to better understand recharge processes and improve groundwater resources management. 相似文献
18.
The static and seismic sliding limit equilibrium condition of retaining walls is investigated, and analytical solutions for the angle of the active slip surface, the critical acceleration coefficient and the coefficient of active earth pressure are provided for different surcharge conditions. In particular, walls retaining a horizontal backfill without surcharge, walls supporting an extended uniform surcharge applied at different distances from the wall and walls supporting a limited uniform surcharge or linear uniform surcharge parallel to the wall are considered in the analysis.The solutions have been derived in the framework of the limit equilibrium approach, considering the effect of the wall through its weight, and accounting for the shear resistance at the base of the wall and the inertia force arising in the wall under seismic conditions.For the wall without surcharge the effect of the vertical component of the seismic acceleration as well as the effects of the inclination of the wall internal face and of the soil–wall friction were also investigated.The angle of the slip plane, the critical seismic acceleration coefficient and the coefficient of active earth pressure are given as functions of dimensionless parameters and the boundary conditions for the applicability of each solution are specified. The influence of soil weight, surcharge conditions and inertia forces on the active earth pressure coefficient is analysed. 相似文献
19.
COEN J. RITSEMA JANNES STOLTE KLAAS OOSTINDIE ERIK VAN DEN ELSEN PAUL M. VAN DIJK 《水文研究》1996,10(8):1081-1089
Surface runoff may be generated when the rainfall intensity exceeds the infiltration capacity, or when the soil profile is saturated with water. Indications exist that both types of overland flow may occur in hilly agricultural loess regions. Here, for a loessial hillslope under maize in the southern part of The Netherlands, it was shown, with pressure head and runoff measurements, that Hortonian overland flow occurs during typical summer rain events. Surface runoff was initiated after saturation of the top 5–10 cm of the soil. Deeper in the soil, unsaturated conditions prevailed while runoff took place. Peak runoff discharges at the outlet of the subcatchment occurred a few minutes after peak rainfall intensities were measured. It appeared that SWMS_2D, a two-dimensional water flow model, was capable in simulating observed pressure head changes and runoff. Simulated potential runoff for the transect studied was higher by a magnitude of three than the measured areal average. This indicates effects of surface ponding, and the probable location of this particular transect in a region with high runoff production. 相似文献
20.
《Continental Shelf Research》2006,26(5):622-652
This paper reviews research on beach groundwater dynamics and identifies research questions which will need to be answered before swash zone sediment transport and beach profile evolution can be successfully modelled. Beach groundwater hydrodynamics are a result of combined forcing from the tide and waves at a range of frequencies, and a large number of observations exist which describe the shape and elevation of the beach watertable in response to tidal forcing at diurnal, semi-diurnal and spring-neap tidal frequencies. Models of beach watertable response to tidal forcing have been successfully validated; however, models of watertable response to wave forcing are less well developed and require verification. Improved predictions of swash zone sediment transport and beach profile evolution cannot be achieved unless the complex fluid and sediment interactions between the surface flow and the beach groundwater are better understood, particularly the sensitivity of sediment transport processes to flow perpendicular to the permeable bed.The presence of a capillary fringe, particularly when it lies just below the sand surface, has influences on beach groundwater dynamics. The presence of a capillary fringe can have a significant effect on the exchange of water between the ocean and the coastal aquifer, particularly in terms of the storage capacity of the aquifer. Field and laboratory observations have also shown that natural groundwater waves usually propagate faster and decay more slowly in aquifers with a capillary fringe, and observations which suggest that horizontal flows may also occur in the capillary zone have been reported. The effects of infiltration and exfiltration are generally invoked to explain why beaches with a low watertable tend to accrete and beaches with a high watertable tend to erode. However, the relative importance of processes such as infiltration losses in the swash, changes in the effective weight of the sediment, and modified shear stress due to boundary layer thinning, are not yet clear. Experimental work on the influence of seepage flows within sediment beds provides conflicting results concerning the effect on bed stability. Both modelling and experimental work indicates that the hydraulic conductivity of the beach is a critical parameter. However, hydraulic conductivity varies both spatially and temporally on beaches, particularly on gravel and mixed sand and gravel beaches. Another important, but poorly understood, consideration in beach groundwater studies is the role of air encapsulation during the wetting of beach sand. 相似文献