You are an expert at summarizing long articles. Proceed to summarize the following text: 
unlike water , a layer of sand will not flow unless its surface is inclined beyond a characteristic angle , known as the maximum angle of stability @xcite . this simple fact translates into a host of threshold phenomena wherever granular material is found . many such phenomena play a crucial role in the erosion of earth s surface , and very likely manifest themselves in the richness of the patterns exhibited by drainage networks .    depending on geological , hydrological , and climatological properties , erosion by water is mainly driven either by overland flow or subsurface flow . the former case occurs when the shear stress imposed by a sheet flow exceeds a threshold @xcite . erosion in the latter case  known as seepage erosion , or sapping  occurs when a subsurface flow emerges on the surface . here the eroding stresses derive not only from the resulting sheet flow but also the process of seepage itself @xcite . the onset of erosion for both overland flow and seepage is threshold - dependent , but the additional source of stress in the case of seepage has the potential to create significantly different erosive dynamics .    here we study the seepage case . whereas the case of horton overland flow has been extensively studied @xcite , seepage erosion has received less attention . @xcite suggests that erosive stresses due to seepage are more widespread in typical environments than commonly assumed . he also provides a detailed description of seepage erosion in the field , together with a discussion of the various factors that influence its occurrence . another focus of attention has been the controversial possibility that many erosive features on mars appear to have resulted from subsurface flows @xcite . although the importance of seepage stresses in erosion have been realized by @xcite and @xcite , comprehensive quantitative understanding is difficult to obtain . the complexity arises from the interdependent motion of the sediment and fluid  the `` two - phase phenomenon '' @xcite which , of course , is common to _ all _ problems of erosion .    to further understand seepage erosion , we proceed from experiments @xcite . questions concerning the origin of ancient martian channels have motivated considerable experimental work in the past @xcite . the process of seepage erosion has also been studied as an example of drainage network development @xcite . our experiments , following those of @xcite and others , are designed to enable us to construct a predictive , quantitative theory . consequently , they stress simplicity and completeness of information . although our setup greatly simplifies much of nature s complexity , we expect that at least some of our conclusions will improve general understanding , and therefore be relevant to real , field - scale problems . a previous paper by @xcite provided a qualitative overview of the phenomenology in our experiment . it described the main modes of sediment mobilization : channelization , slumping , and fluidization . here we provide quantitative understanding of the onset and transitions between these modes . our emphasis is on the threshold phenomena associated with the onset of erosion , which we will ultimately characterize in the same way that others @xcite have characterized the onset of dry granular flow beyond the maximum angle of stability . this involves a construction of a generalized shields criterion @xcite valid in the presence of seepage through an inclined surface . a major conclusion is that the onset of erosion driven by seepage is significantly different from the onset of erosion driven by overland flow . we find that there is a critical slope @xmath0 , significantly smaller than the maximum angle of stability , above which the threshold disappears . therefore any slope greater than @xmath0 is unstable to erosion if there is seepage through it . this result is similar to well - known conclusions for the stability to frictional failure of slopes with uniform seepage @xcite . an important distinction in our work , however , concerns the mode of sediment mobilization and its local nature . the existence of the critical slope for seepage erosion may provide a useful quantitative complement to the qualitative distinctions between seepage and overland flow that have already been identified @xcite .    the remaining modes of sediment mobilization , fluidization and slumping , are modeled using well established ideas @xcite . the result of applying these ideas together with the generalized shields criterion provides a theoretical prediction of the outcomes of the experiment , i.e. , a phase diagram . agreement between theory and experiment is qualitative rather than quantitative . we nevertheless believe that our theoretical approach is fundamentally sound and that better agreement would follow from improved experimental procedures . in our experimental setup , first introduced by @xcite , a pile of identical cohesionless glass beads @xmath1 mm in diameter is saturated with water and compacted to create the densest possible packing . it is then shaped into a trapezoidal wedge inclined at an angle @xmath2 with slope @xmath3 as shown in fig.[fig : expt ] . the downslope length of the wedge is @xmath4 cm , its width across the slope is @xmath5 cm , and its height in the middle is approximately @xmath6 cm . water enters the sandpile underneath through a fine metal mesh and exits at the lower end of the pile through the same kind of mesh . a constant head at the inlet is maintained by keeping a constant water level @xmath7 in the reservoir behind the sandbox with the help of an outflow pipe . the slope @xmath8 of the pile and the water level @xmath7 are the control parameters of the experiment . the degree of packing of the granular pile is the variable most difficult to control . our particular method of feeding water into the sandpile , similar to that of @xcite , can be motivated in three ways . the most important justification is the fact that the amount of water flowing on the surface can be finely controlled in our geometry . this feature is essential in probing the onset of erosion . second , our setup allows us to access heads @xmath7 larger than the height of the pile , which therefore allows us to explore dynamic regimes unavailable if water enters the pile through a mesh in the back . third , a similar seepage water flow geometry can exist in the field wherever water travels beneath an impermeable layer that terminates . we have performed two types of experiments : steady and non - steady . for a fixed water level and in absence of sediment motion , water flow reaches steady state . by monitoring the total water flux through the system we estimate the time to reach steady state to be approximately ten minutes . to explore the onset of sediment motion , we raised the water level @xmath7 in small increments , waiting each time for steady state to be established . due to the particular shape of the bulk flow in our experiment , surface flow exists over a finite region of the surface . the width of this seepage face and therefore the depth of the surface flow can be tuned by changing @xmath7 . because of the finite extent of surface flow , its depth and therefore the viscous shear stress reaches a maximum at a certain location . thus , by increasing @xmath7 we can continuously tune the maximum shear stress experienced by the surface grains . the maximum shear stress reaches a critical value for the onset of sediment motion in a certain location on the slope . as we show below , we can compute where the maximum shear stress occurs and thus can reliably detect the onset of sediment motion visually because we focus our attention on this location . once sediment begins to move , channels form almost immediately . these channels grow in length , width , and depth . an example of the evolving channel network is shown in fig.[fig : channels ] . depending on the slope , as the channels deepen , the pile becomes unstable to fluidization or slumping . for slopes lower than approximately 0.05 , the fluidization threshold is reached before sediment is mobilized on the surface . m. the slope of the pile is @xmath9 and the water level @xmath10 cm.,width=316 ]    we also explored the non - steady evolution of the bulk and surface water flow and resulting sediment motion by raising the water level @xmath7 to some higher value from zero . in this case one of three things can happen . the pile can be fluidized within a few seconds or fail by slumping as shown in fig.[fig : slump ] . if this does not occur , the water emerges on the surface just above the inlet . a sheet of water then washes down the slope of the pile . during this initial wash , sediment is mobilized and incipient channels form . these channels grow during subsequent relaxation of the bulk water flow towards steady state . because of the initial wash s erosive power , channels are able to form and grow for lower water pressures than in steady experiments .     angle to the slope after the water flow has been stopped . the width of the imaged region is approximately @xmath11 m . slumping happens along a convex upward arc which looks darker because it is deeper and therefore wetter.,width=316 ]    outcomes of a large number of non - steady experiments and several steady experiments for varying slope @xmath8 and the water level @xmath7 are summarized in the phase diagram in fig.[fig : phase ] . each symbol in the plot represents one experiment . the sediment is either immobile ( stable seepage ) , or it is mobilized on the surface where channels form ( channelization ) or in the bulk ( slumping or fluidization ) . in several experiments , slumping or fluidization happened after channels formed and grew . in the following sections we describe the computations that allow us to construct the theoretical boundaries between the three different modes of sediment mobilization in our experiment .    ; those that produced channels are indicated by @xmath12 ; and those that produced fluidization and/or slumping within one hour of the beginning of the experiment are represented by @xmath13 . the straight line and gray - shaded curves are theoretical predictions for the boundaries separating the four regions indicated by their labels . the thickness of the lines indicates uncertainty in the theory . the boundary between the uneroded and channelized states is reasonably well approximated by our theory . the theoretical boundaries for fluidization and slumping , however , appear to overestimate the critical water level , possibly as a result of inhomogeneities , dynamic changes in the sandpile s shape , or from the assumption of a steady state.,width=316 ] whereas steady - state flow can be readily characterized quantitatively , non - steady flow characterization requires knowledge of the water - table dynamics . however , the theory of the water - table dynamics is less well established than that of the flow through the bulk of a porous medium . also , our steady - state experiments probe all aspects of sediment dynamics . we can therefore focus on the quantitative characterization of the steady - state flow .    to study the onset of erosion quantitatively we need to be able to establish a correspondence between the experimentally measurable quantities such as the slope @xmath8 , the water level @xmath7 , the size of the seepage face , and the water fluxes . the seepage and surface fluxes are the most difficult to measure . in this section we set up their computation . the computation is designed to enable us to infer water fluxes indirectly by measuring the size of the seepage face . in the following sections we will use this computation to quantify the onset of erosion and to compute the slumping and liquefaction boundaries of the channelization phase diagram shown in fig.[fig : phase ] . fig.[fig : expt ] specifies the key quantities and coordinate systems we use in computing the fluxes . the flow profile is independent of the @xmath14-coordinate across the slope of the sandpile except near the side walls of the box . we therefore treat the box as if it were infinitely wide . flow is then two - dimensional and the specific discharge vector @xmath15 is in the @xmath16-@xmath17 plane . we will use two coordinate systems . as shown in fig.[fig : expt ] , the @xmath18 coordinate is measured vertically from the bottom of the box while the @xmath17 coordinate is the normal distance away from the surface of the pile . the flow is governed by darcy s law , @xmath19 where @xmath20 is the scalar hydraulic conductivity , and @xmath21 is the total hydraulic head of the pore water . both @xmath15 and @xmath20 have units of velocity while the scaled pore pressure @xmath22 has units of length . here @xmath23 is the density of water and @xmath24 is the magnitude of the acceleration of gravity . we have measured @xmath20 via a -tube relaxation experiment . to do so we created a water level difference @xmath25 between the two arms of a transparent -shaped tube of width @xmath26 partially filled with glass beads . the rate of change of @xmath25 is given by @xmath27 . by measuring the rate of change of @xmath25 we deduced the value of the hydraulic conductivity @xmath28 mm / s ( @xcite ) . hydraulic conductivity is sensitive to the packing of the grains and is the variable most difficult to control in our experiment . water incompressibility implies @xmath29 , therefore yielding laplace s equation , @xmath30    to compute the pore pressure @xmath31 , boundary conditions must be specified . the walls of the box are impenetrable . therefore the discharge vector is parallel to the walls . in other words , the flux @xmath32 in the direction @xmath33 normal to the walls vanishes . thus @xmath34 . because the glass beads in our experiment are small , capillarity is important . when a tube filled with glass beads is lowered into a reservoir of water , the porous bead - pack fully saturates in a region that extends above the surface of the water by a capillary rise @xmath35 . we measured @xmath36 mm for our material . the capillary rise is a measure of the average radius of the water menisci at the edge of the fully saturated zone . the pore pressure at the edge of the fully saturated zone is @xmath37 ( without loss of generality we set the atmospheric pressure to zero ) . water can rise above the fully saturated zone through the smaller pores and narrower throats . thus a partially saturated capillary fringe exists above the fully saturated zone . however , in this fringe the water is effectively immobile since it is confined to the smaller pores and narrower throats . since water flows only in the fully saturated zone , we define the water table to be at its edge . thus , the pore pressure at the water table is equal to the negative capillary rise @xmath38 . in steady state the discharge vector is parallel to the water table . this extra condition allows us to determine the location of the water table in steady state . we neglect the pressure drop across the inlet mesh . therefore , the pore pressure at the inlet mesh is @xmath39 . the boundary conditions at the surface of the sandpile and at the outlet mesh are more subtle . when no water seeps out , i.e. , when the discharge vector is parallel to the surface , the curvature of the water menisci between grains can freely adjust so that the pressure @xmath31 can vary between zero and @xmath37 . therefore when @xmath40 , no seepage occurs . otherwise , the pore pressure equals the atmospheric pressure @xmath41 ( we neglect the pressure exerted by the thin layer of water on the surface ) , and the discharge vector has a component normal to the surface , i.e. , there is either exfiltration or infiltration .    , or the pressure at the water table when it is below the surface of the pile . note that seepage occurs only where the pore pressure reaches atmospheric pressure . slope is @xmath9 , water level @xmath42 cm.,width=316 ]    to obtain the steady state location of the water table , we guess its position and solve laplace s equation   with the @xmath43 boundary condition on the water table . we then move the water table in the direction of the local discharge vector by an amount proportional to its length . iteration of this procedure converges to the steady - state position of the water table . an example is shown in figure [ fig : table_pressure_overland ] . once the steady flow pattern is known , we can calculate the overland water flux @xmath44 by integrating the one - dimensional continuity condition which states that the downslope derivative of the overland flux is equal to the seepage flux : @xmath45 in this section we assume , based on direct observation , that the onset of channel incision coincides with the onset of erosion ( i.e. , we never observed a homogeneously eroding state ) . in other words , when the overland water flux becomes strong enough to carry grains , the flow of sediment becomes immediately unstable to perturbations transverse to the downslope direction and incipient channels form @xcite . using this assumption and the calculation of the overland water flux we can deduce the threshold condition for the onset of erosion . it is universally assumed after @xcite that the hydrodynamic stresses exerted on the sandpile by the fluid flowing on its surface determine whether cohesionless granular material is entrained . in the limit of laminar flow , the dominant hydrodynamic stress is the viscous shear stress @xmath46 . appropriately scaled this shear stress is termed the shields number @xcite , defined by @xmath47 where @xmath48 is the density of the granular material , @xmath49 is the grain diameter ( @xmath50 mm in our experiment ) , and the surface is not inclined . the conventional shields number   is the ratio between the horizontal force exerted by the flow and vertical force due to grain s weight . to generalize the notion of the shields number to the situation with seepage through an inclined surface , we make two changes in eq .  . we first add the tangential component of the seepage force density @xmath51 acting over a length @xmath49 to the numerator of  . the numerator thus becomes @xmath52 . note that we did not include the tangential component of the grain s weight to the numerator . defined in this way , the generalized shields number measures the effect of the fluid : both the bulk as well as the surface flows . -axis of the resultant of the grain s weight and the seepage force both scaled by @xmath53.,width=220 ]    second , we replace the denominator of eq .   by the resultant ( vectorial sum ) of the seepage force on a grain and its submerged weight , as shown in fig.[fig : mod_shields ] , both scaled by @xmath53 ( @xmath54 to obtain stress as in the numerator and @xmath55 for agreement with the conventional shields number ) , projected onto the @xmath17-axis . according to @xcite the grains on the surface of the bed experience a seepage force roughly half as large as the grains several layers deep . consequently , we assume that the seepage force is reduced by a factor of @xmath56 ; therefore @xmath57 where @xmath2 is the inclination angle of the surface . the importance of the seepage stresses for the criterion for the onset of erosion was previously realized by @xcite . it can be shown that their equation ( 10 ) expressing marginal stability of a surface grain is equivalent to writing @xmath58 , the tangent of the angle of internal friction . the generalized shields number eq .   is a measure of the relative importance of the tangential and normal forces acting on a grain at the surface of the sandpile . therefore , we expect @xmath59 to be a control parameter for erosion . in other words , there exists a critical shields number @xmath60 , such that when @xmath61 , surface grains are immobile , and when @xmath62 , sediment is mobilized . note that ( [ eq : shields ] ) reduces to the classical definition of the shields number for a flat surface without seepage . also note that since we did not include the tangential component of grain s weight , the critical shields number at the onset of sediment motion vanishes when the inclination angle @xmath2 reaches the maximum angle of stability . although we obtain the seepage force density @xmath63 as a result of computing the pore - water pressure , to calculate the boundary shear stress @xmath46 we must estimate the thickness of the surface water layer @xmath64 . since this thickness changes slowly in the downslope direction , we can approximate the surface flow by the steady flow of a uniform layer of viscous fluid . also , the surface water flux is small enough for turbulence to be of no importance . the thickness @xmath65 of laminar surface flow for a given flux @xmath66 is @xcite @xmath67 where @xmath68 is the viscosity , while the viscous shear stress exerted on the sandpile is @xmath69 the particle reynolds number can then be calculated using the bottom shear stress   and shear velocity @xmath70 as @xmath71 where @xmath72 is the kinematic viscosity of water . we estimate that in our experiments , this particle reynolds number varies between 5 and 20 depending on the slope @xmath8 of the pile and the water level @xmath7 . we verify this estimate of the reynolds number by a direct measurement of the thickness of the surface flow . we find that this thickness is several grain diameters . this justifies the laminar flow assumption used in obtaining eq .  . using  , the shields number   can now be conveniently rewritten as @xmath73 this expression can be further simplified by noting that @xmath41 along the seepage face . therefore @xmath74 at the surface wherever there is overland flow . we arrive at the final expression for the modified shields number which depends on the surface flow thickness @xmath65 , the normal component @xmath75 of the seepage force density at the surface , and the seepage force reduction factor @xmath76 @xmath77    in our geometry , both the surface and the seepage water fluxes reach a maximum somewhere along the slope . therefore the shields number has a maximum value as well . below we calculate this maximum shields number in steady state for a given slope @xmath8 and water level @xmath7 . for a pile of slope @xmath9 . at @xmath78 water first seeps through the surface and the shields number jumps to a nonzero value . afterwards it increases rapidly as @xmath79 ( solid line ) . inset : corresponding size of the computed seepage face.,width=316 ]    we now explore the consequences of seepage for the phenomenology of the onset of erosion . because of the additional force on the surface grains , seepage flow is more erosive than overland flow . this notion is reflected quantitatively in the generalized shields number . let us examine how the maximum shields number varies with the water level @xmath7 in our experiment . a representative graph of the maximum shields number versus the water level is shown in fig.[fig : shields_h ] . below a water level @xmath80 that is a function of the slope , no water seeps out to the surface of the pile . even though the water table may be at the surface , the pressure at the water table is below atmospheric pressure and capillarity prevents seepage . when @xmath81 , i.e. , exactly at the onset of seepage , the pressure reaches @xmath41 at some point on the surface . since the seepage flux is still zero , @xmath82 along the wet part of the surface . therefore , just above the seepage onset , when the water layer thickness and the seepage flux are both infinitesimally small , the maximum shields number is @xmath83 in contrast to overland flow , the consequence of seepage is that as soon as the water emerges on the surface , the maximum shields number is some non - zero value which depends on the slope . this also implies that there exists a critical slope @xmath0 such that @xmath84 is equal to the critical shields number @xmath60 , i.e. , @xmath85 for slopes greater than @xmath0 seepage is always erosive . note that for low - density particles this critical slope can be arbitrarily small . the expression for the critical slope for seepage erosion in eq .   is analogous to well - known formulas for stability of slopes to coulomb failure due to uniform seepage @xcite . our result applies locally to the point where non - uniform seepage first emerges on the surface . in this situation , the pile is generally stable to coulomb failure and the sediment is eroded only locally on the surface .     as a function of the downslope coordinate @xmath16 at the onset of seepage @xmath86 and just above.,width=220 ]    we now show that above @xmath87 , the maximum shields number increases rapidly as a @xmath88 power of the water level excess @xmath89 . at the onset of seepage , i.e. , when @xmath81 , the pressure at the water table reaches atmospheric pressure @xmath41 at some point @xmath90 located at @xmath91 and @xmath92 on the surface . even though the water table is at the surface , there is no seepage anywhere , i.e. , @xmath93 . because the pressure is smooth , it can be approximated by a quadratic function near this point so that @xmath94 , where @xmath95 and @xmath96 are constants with appropriate dimensions . when the water level is raised by a small increment @xmath97 , the lowest order change in the head at the water table is an increase of @xmath98 with the exception of the region where this increase would lead to a positive pressure . as illustrated in fig.[fig : increment ] , in this region the pore pressure @xmath31 is set to @xmath99 , and thus this region becomes the seepage face . the width of the seepage face @xmath100 scales like the square root of @xmath98 , i.e. , @xmath101 as seen in the inset of fig.[fig : shields_h ] . the seepage flux can be estimated by noting that the hydraulic head is modified by an amount @xmath98 over a vertical region of order @xmath100 . therefore we obtain @xmath102 . the total surface flux therefore scales like the product of the seepage flux and the width @xmath100 of the seepage face , i.e. , @xmath103 . the lowest order change in the maximum shields number is due to the change of the surface flow depth @xmath104 . thus as we claimed above , just above the water level @xmath87 for the onset of seepage , @xmath105 where the constant @xmath106 is a function of the slope . variation with water level of the computed maximum shields number shown in fig.[fig : shields_h ] is consistent with eq .  . in the previous sections we have detailed the way of calculating the bulk and surface water fluxes in our experiment and the resulting maximum generalized shields number . in this section we use this calculation to examine the onset of the sediment flow and channelization . our first goal is to measure the threshold or critical shields number required for the mobilization of sediment . we then use this measured value of the critical shields number to predict the outcome of steady - state experiments for various values of the slope and the water level and thus compute the channelization boundary in the phase diagram in fig.[fig : phase ] . the actual maximum shields number in the experiment differs from the quantity calculated in eq . ( [ eq : shields - prefinal1 ] ) . in addition to random errors in the measurements of the pile dimensions and water level , there are several sources of systematic error . for example , the pressure drop across the inlet mesh results in a lower effective hydraulic head . also , our measurement of the capillary rise @xmath35 is dependent on a visual estimate of the fully saturated zone and thus can be a source of systematic error . we indeed find that the size of the seepage face calculated for a particular water level @xmath7 is greater than measured in the experiment . however , the size of the seepage face translates directly into the surface water flux and therefore the maximum shields number . the inset of figure [ fig : shields_h ] shows the typical dependence of the size of the seepage face on the water level . the variation of the maximum shields number with the size of the seepage face is shown in fig.[fig : shields_seep ] for three different slopes . we use this computed correspondence between the size of the seepage face and the maximum shields number to infer the maximum shields number in the experiment by measuring the size of the seepage face .    to measure the critical shields number we raise the water level @xmath7 in increments of a few millimeters at a time . each time the water level is increased , the seepage flow is allowed to reach a steady state . in each of these steady states we measure the seepage face size and infer the corresponding maximum shields number . eventually , sediment is mobilized and we record the size of the seepage face and compute the corresponding maximum shields number . this number is an upper bound on the critical shields number for our granular material at that particular slope . the lower bound on the critical shields number is obtained from the largest seepage face at which no sediment is moving or sediment motion is only transient . averaging over several experiments with slope @xmath9 we estimate the critical shields number to be @xmath107 it is not obvious that the generalized critical shields number for the onset of seepage driven erosion should coincide with the critical shields number for overland flow . however our measured value of the critical generalized shields number is within the scatter of the existing data for overland flow summarized in @xcite . our measurement the critical generalized shields number is equivalent to measuring the angle of internal friction due to the correspondence of our definition of @xmath59 and howard and mclane s equation ( 10 ) . deviations from flatness of the pile s surface result in the fluctuations of the thickness of the surface water film . as a result , the maximum bottom shear stress in the experiment is systematically greater than that calculated at a given size of the seepage face . thus the shields number calculated for a particular size of the seepage face is the lower bound on the actual shields number in the experiment .    in principle , the critical shields number should vary with the slope of the pile . evidence for this is the fact that at the maximum angle of stability any additional forcing from the water flowing over the bed mobilizes sediment . since it is reasonable to assume that the critical shields number is continuous and monotonic , we arrive at the notion that it decreases monotonically with slope and vanishes at the maximum angle of stability . for small slopes the critical shields number is expected to decrease as the cosine of the inclination angle since this is the lowest order change in the stabilizing effect of gravity . for most slopes in our experiments , @xmath108 is within a few percent of unity and thus we can ignore the variation of the critical shields number with slope . this assumption allows us to predict the water level at which erosion and therefore channelization should commence in our experiment . in fig.[fig : phase ] a boundary is drawn between regions where sediment is expected to mobilize and remain immobile . to obtain this line we computed for each slope the water level @xmath109 at which the shields number is equal to the critical shields number . below this water level , i.e. , when @xmath110 , the maximum shields number is below critical and thus sediment is immobile . conversely , for @xmath111 , the maximum shields number is above critical and thus sediment is mobilized and channels form . the channelization boundary is widened because of the uncertainty in the critical shields number . qualitative agreement of the channelization boundary with experiments is perhaps due to the opposite action of two effects . first , channelization occurs for lower water levels in non - steady experiments . this happens because in non - steady experiments the maximum shields number overshoots its steady - state value . the overshoot is greatest for small slopes . second , a pressure drop across the inlet mesh and the compacted region of sand close to it has an opposite effect which increases the water level needed for channelization . these two effects , though small , could together affect the accuracy of our predictions . since these effects act in opposite ways , our the predictions of the calculated channelization water level @xmath109 agree qualitatively with the experiments . having computed the channelization boundary in the phase diagram , we now pursue a quantitative description of the other two modes of sediment mobilization exhibited by our sandpile . higher water pressures can cause the sandpile to fail in one of two ways . first , an upward seepage force can lift sand and result in a fluidization or quicksand instability . second , the pile can become unstable to slipping , slumping , or sliding . both failure mechanisms have been discussed by a number of studies , e.g. , those of @xcite or @xcite . and height @xmath26 and compare it with its weight.,width=144 ]    fluidization occurs when at some point @xmath90 in the sandpile the pore pressure is larger than the total hydrostatic pressure due to the weight of the sand and the water above @xmath90 . to see this we compute the total seepage force acting on a slice of sand of width @xmath112 between point @xmath90 and point @xmath113 on the surface of the pile directly above @xmath90 . the vertical component of this force is ( see fig.[fig : fluidize ] ) @xmath114 where @xmath26 is the height of the slice . when this force exceeds the submerged weight @xmath115 of the slice , the slice is lifted and the bed is fluidized . here @xmath116 is the total density of the saturated sand , which for our sand is approximately 2 g/@xmath117 . thus fluidization occurs when there exists points @xmath90 and @xmath113 on the surface directly above @xmath90 such that @xmath118 for uniform seepage this condition is equivalent to those in @xcite and @xcite .    to construct the fluidization boundary in the phase diagram ( fig.[fig : phase ] ) , we find the water level @xmath119 above which there exists at least one point in the pile for which condition   is satisfied . below this fluidization water level this condition is not satisfied for any point in the pile . in addition to fluidization the sandpile can fail by slumping . this can happen in one of two ways . frictional failure can occur in the bulk of the pile due to the seepage stresses . alternatively , surface avalanching can occur . to establish an upper bound on the water level at which the sandpile slumps via either mechanism we use the criterion developed by @xcite for determining when a slope is destabilized by uniform groundwater seepage . essentially it requires calculating the vectorial sum of the seepage and gravity forces acting on a small element of soil near the surface . when the angle between this total force and the downward normal to the surface , which we will call the effective inclination angle , exceeds the maximum angle of stability , the surface grains are destabilized . we measured the maximum angle of stability to be @xmath120 for dry glass beads . the slumping boundary in the phase diagram ( fig.[fig : phase ] ) is constructed by computing the effective inclination angle along the surface of the pile and noting the water level @xmath121 , at which the effective inclination angle reaches the maximum angle of stability @xmath122 at some point of the surface . figure [ fig : phase ] shows the critical water level at which fluidization and slumping should occur according to the criteria above . failure occurs at systematically lower water levels in the experiment . there are several effects which can account for this difference . first , any irregularities in the construction of the pile such as voids or surface height fluctuations make the pile more unstable to fluidization and slumping . second , we compute the instability of an uneroded pile , whereas in most experiments , the pile failed after erosion had changed the shape of the pile . the decrease of pile s height due to erosion increases the head gradient in the bulk and thus makes the pile more prone to slumping and/or fluidization .     and three different water levels . for the highest water level , a region on the surface has an effective angle above the maximum angle of stability and thus the slope is unstable to slumping . the inset shows the plateau value of the effective inclination angle as a function of slope . when the plateau value reaches the maximum angle of stability , even a small amount of seepage destabilizes the pile to slumping.,width=316 ]    at @xmath123 , a jump in the slumping water level is observed in both the experiment and the model . this jump is a purely geometric effect . slumping occurs when , somewhere along the slope , the effective inclination angle , which includes the effect of the seepage force , exceeds the maximum angle of stability . as shown in fig.[fig : slump_jump ] , for slopes smaller than @xmath124 , the effective inclination angle is flat and develops a peak under the water inlet as the water pressure is increased . when the top of this peak crosses the value of the maximum angle of stability , the pile slumps . when the slope exceeds @xmath124 , however , the value of the plateau in the effective inclination angle is above the maximum angle of stability . therefore , for these slopes , the pile will be unstable to slumping as soon as the water emerges on the surface . this article reports on our progress in understanding seepage erosion of a simple non - cohesive granular material in a laboratory - scale experiment introduced in @xcite . our ultimate goal is to construct a quantitative predictive theory of the onset and growth of the channel network observed in this experiment . this goal requires a complete sediment transport model as well as the calculation of the relevant water fluxes . here we obtain the latter and focus on the onset of erosion . prediction of the onset of erosion based on the generalized shields conjecture explains qualitatively the channelization boundary in the experimental phase diagram . by invoking well established simple ideas we also roughly explain the fluidization and slumping boundaries in the phase diagram . greater discrepancy with the experiment for these boundaries indicates that better understanding of the slumping / fluidization mechanisms particular to our experiment is needed . the central result of our exploration is the introduction of the generalized shields criterion for seepage erosion . as a consequence of seepage forces on the surface grains , the threshold for the onset of erosion driven by seepage is slope dependent . the threshold disappears at a critical slope @xmath0 determined by the critical shields number for overland flow and the density contrast between the granular material and water . in most cases this critical slope is significantly smaller than the maximum angle of stability . we find , therefore , that slopes above this critical slope are unstable to any amount of seepage . as a consequence , slopes that sustain seepage must be inclined at an angle smaller than the critical angle for seepage erosion . this behavior contrasts strongly with the threshold phenomena in erosion by overland flow , and therefore provides a mechanistic foundation for distinguishing the two types of erosion . this work was supported by a doe grants de - fg02 - 99er15004 and de - fg02 - 02er15367 . aharonson , o. , m.  t. zuber , d.  h. rothman , n.  schorghofer , and k.  x. whipple , drainage basins and channel incision on mars , _ proceedings of the national academy of sciences usa _ , _ 99 _ , 17801783 , 2002 . baker , v.  r. , spring sapping and valley network development , in _ groundwater geomorphology : the role of subsurface water in earth - surface processes and landforms _ , edited by c.  higgins and d.  coates , chap .  11 , geological society of america , boulder , colorado , 1990 . buffington , j.  m. , and d.  r. montgomery , a systematic analysis of eight decades of incipient motion studies , with special reference to gravel - bedded rivers , _ water resources research _ , _ 33_(8 ) , 19932029 , 1997 . dunne , t. , k.  x. whipple , and b.  f. aubry , microtopography and hillslopes and initiation of channels by horton overland flow , in _ natural and anthropogenic influences in fluvial geomorphology _ , pp . 2744 , american geophysical union , 1995 . howard , a.  d. , groundwater sapping experiments and modelling , in _ sapping features of the colorado plateau , a comparative planetary geology field guide _ , edited by a.  d. howard , r.  c. kochel , and h.  e. holt , pp . 7183 , nasa scientific and technical information division , washington , d.c . , 1988 . kochel , r.  c. , a.  d. howard , and c.  mclane , channel networks developed by groundwater sapping in fine - grained sediments : analogs to some martian valleys , in _ models in geomorphology _ , edited by m.  j. woldenberg , pp . 313341 , allen and unwin , boston , 1985 . kochel , r.  c. , d.  w. simmons , and j.  f. piper , groundwater sapping experiments in weakly consolidated layered sediments : a qualitative summary , in _ sapping features of the colorado plateau , a comparative planetary geology field guide _ , edited by a.  d. howard , r.  c. kochel , and h.  e. holt , pp . 8493 , nasa scientific and technical information division , washington , d.c . , 1988 . shields , a. , _ anwendung der hnlichkeitsmechanik und der turbulenzforschung auf die geschiebebewegung _ , heft 26 , mitteilung der preussischen versuchsanstalt fr wasserbau und schiffbau , berlin , germany , ( in german ) , 1936 .
we study channelization and slope destabilization driven by subsurface ( groundwater ) flow in a laboratory experiment . 
 the pressure of the water entering the sandpile from below as well as the slope of the sandpile are varied . 
 we present quantitative understanding of the three modes of sediment mobilization in this experiment : surface erosion , fluidization , and slumping . 
 the onset of erosion is controlled not only by shear stresses caused by surfical flows , but also hydrodynamic stresses deriving from subsurface flows . 
 these additional forces require modification of the critical shields criterion . whereas surface flows alone 
 can mobilize surface grains only when the water flux exceeds a threshold , subsurface flows cause this threshold to vanish at slopes steeper than a critical angle substantially smaller than the maximum angle of stability . 
 slopes above this critical angle are unstable to channelization by any amount of fluid reaching the surface .