COURSE NOTES AND STUDY GUIDE FOR GEO 432/532 APPLIED GEOMORPHOLOGY

Size: px
Start display at page:

Download "COURSE NOTES AND STUDY GUIDE FOR GEO 432/532 APPLIED GEOMORPHOLOGY"

Transcription

1 COURSE NOTES AND STUDY GUIDE FOR GEO 432/532 APPLIED GEOMORPHOLOGY STEPHEN T. LANCASTER, ASSOC. PROF., CEOAS 1. General Advice for Exam Preparation I ve included nearly 100 equations in the summary below, not because I expect you to memorize them but, rather, because math is the most straightforward way to represent the concepts, and I expect you to understand those concepts, so I expect you to understand the equations. Many of those equations are variations on a theme, such as the Exner equation. Many others are equations that you have had to use during the term for your assignments or that we have visited repeatedly in class. Regarding what might be considered memorization, things you should know include (but are not necessarily limited to) the basic Exner equation and what it means such that you could, say, make modifications to make it fit some special cases, like some of those in the last few topics covered; the basic form of most bedload sediment transport equations; the most common equations used to calculate bed shear stress; the common implications of the normal flow approximation (i.e., what relationships are valid). Some of the key concepts for the course include (but are not necessarily limited to) sediment and grain size distributions; bankfull channels; normal flow; thresholds for initiation of motion; bed material sediment transport; hydraulic roughness due to skin friction and form drag; morphodynamics and the Exner equation (conservation of mass). Date: Spring

2 2 LANCASTER 2. Course Summary What follows is a summary of the material covered during the term relative to the outline presented in the syllabus. 1. Introduction: Felix Exner and river morphodynamics (Parker (2004), Ch. 1) Morphodynamics: the class of problems for which the flow over a bed interacts strongly with the shape of the bed, both of which evolve in time. Exner equation describes conservation of bed sediment, (1 λ p ) η t = q t x, (1) which is coupled with the St. Venant shallow water equations for conservation of water mass, H t + UH = 0, (2) x and momentum, UH t + U 2 H x = gh H x gh η x C fu 2, (3) in the downstream direction and a relation between sediment transport rate and flow hydraulics, which Exner wrote as q t = q t (U) (4) Equations (1), (2), (3), and (4) are linked via the sediment transport rate, q t, the flow velocity, U, and the bed elevation, η. 2. Sediment (Parker (2004), Ch. 2) a. Grain-size distributions ψ-scale: base-2 logarithmic scale used for grain size distribution (why?): ψ = φ = log 2 (D), (5) where D s units are mm. Important size ranges for this course are sand (D = mm, ψ = 4 1), gravel (D = 2 64 mm, ψ = 1 6), and cobbles (D = mm, ψ =6 8). Characterize distributions in terms of N + 1 boundary sizes, D b,i, such that f f,i is the mass fraction that is finer than D b,i and N grain size ranges with sizes, and fractions, in each size range. ψ i = 1 2 (ψ b,i + ψ b,i+1 ), D i = (D b,i D b,i+1 ) 1/2, (6) D x size at which x% is finer, where, e.g., x = 50, 90, etc. f i = f f,i+1 f f,i, (7)

3 Geometric mean size: Geometric standard deviation: GEO 432/532 COURSE NOTES 3 D g = 2 ψ, ψ = σ g = 2 σ, σ 2 = N ψ i f i. (8) i=1 N ( ψi ψ ) 2 fi. (9) i=1 Submerged specific gravity of sediment, b. Fall velocity R = ρ s ρ 1. (10) Idea is to balance gravitational and drag forces to find terminal velocity: F g = 4 ( ) 3 D 3 πρrg = F D = 1 ( ) 2 D 2 2 πρc D, c D = c D (Re vp), Re vp = v sd 2 ν We can solve for the nondimensional fall velocity, R f = v s / RgD: [ ] 1/2 4 R f = (12) 3c D (Re vp ) We also find that Re vp = R f Re p, where Re p = RgDD/ν, which is the particle Reynolds number. So now the equation gets complicated, partly because the drag coefficient, c D, is a complicated function of Re vp. c. Modes of transport Bed material load: sediment load that exchanges with the bed (and contributes to morphodynamics). Subdivided into bedload and suspended load. Wash load is transported without exchange with the bed (generally silt and clay). Bedload: sliding, rolling, or saltating in ballistic trajectory just above bed, indirect role of turbulence. Suspended load: feels direct dispersive effect of eddies, may be wafted high into water column. 3. Bankfull channel characteristics (Parker (2004), Ch. 3) a. Bankfull concept: Discharge at which flow spills significantly onto floodplain. b. Gravel and sand beds Sand-bed streams have geometric mean (or median) bed sediment size, D sg 2 mm. Gravel-bed streams have D sg > 16 mm, encompasses cobble- and boulder-bed streams. Transitional streams have 2 < D sg 16 mm; few in this category. (11)

4 4 LANCASTER c. Characteristic dimensionless parameters: Dimensionless bankfull discharge: ˆQ = Q bf gds50 D 2 s50 (13) Dimensionless bankfull depth: Good relationship with ˆQ Ĥ = H bf D s50 (14) Dimensionless bankfull width: Good relationship with ˆQ ˆB = B bf D s50 (15) Bankfull Froude number: describes transition between subcritical (Fr < 1) and supercritical (Fr > 1) flow Fr bf = Q bf B bf H bf ghbf = U bf ghbf (16) Bankfull Shields number (estimated): Typically 0.05 for gravel-bed and 1.9 for sand-bed streams τ bf50 = H bfs RD s50 (17) Bankfull Chezy resistance coefficient: note that it s actually an inverse roughness, i.e., larger values mean lower resistance (also note how similar it is to Froude number) Cz bf = Q bf B bf H bf ghbf S = U bf ghbf S Particle Reynolds number (surrogate for grain size): (18) Re p50 = RgDs50 D s50 ν (19) 4. Conservation of bed sediment (more Exner-type equations) (Parker (2004), Ch. 4) Note list of variables and parameters and explanation of coordinate system a. Bedload (uniform grain size): nearly identical to equation (1) (1 λ p ) η t = q b x (20) b. Bedload + suspended load (1 λ p ) η t = q b x + v s (c b E) (21) c. Bedload + suspended load + subsidence or uplift ( ) η (1 λ p ) t + σ = q b x + v s (c b E) (22)

5 d. Sediment in suspension t GEO 432/532 COURSE NOTES 5 ( H 0 ) c dz + q s x = v s (E c b ) (23) e. Total load + subsidence or uplift (1 λ p ) ( ) η t + σ = q t x (24) f. Active layer concept: Active, exchange, or surface layer approximation: g. Mixtures Sediment grains in active layer (thickness L a )have constant, finite probability of entrainment. Sediment grains below active layer (in substrate) have zero probability of entrainment. Size fractions: F i in active layer have no vertical structure f i in substrate do not vary in time interfacial exchange f li at substrate-active layer boundary p bi = q bi /q bt in transport, q bt = N q bi (25) i=1 Conservation of mixture with active layer: [ (1 λ p ) f li t (η L a) + ] t (F il a ) = q bt p bi x Sum over i just gives equation (1), and for evolution of active layer grain size distribution, multiply (1) by f li and subtract from (26): [ F i (1 λ p ) L a t + (F i f li ) L ] a = q bt p bi q bt + f li (27) t x x where the exchange fractions depend on whether there s erosion or deposition: η f i z=η La, f li = t < 0 η (28) αf i + (1 α) p bi, t > 0, 0 α 1 5. River hydraulics (Parker (2004), Ch. 5) (26)

6 6 LANCASTER a. Boundary resistance and steady, uniform flow Hydraulic radius and wide-channel approximation: R h = HB B + 2H = H H H, B H (29) B Shields number: scales as the ratio of drag (impelling) and friction (resisting) forces τ = τ b (30) ρrgd Average velocity: Shear velocity: Dimensionless resistance coefficient: Chezy resistance coefficient: U = u = C f = Q BH τb ρ (31) (32) τ b ρu 2 (33) Cz = U u = C 1/2 f (34) Resistance for hydraulically rough flow, Keulegan: Cz = C 1/2 f = 1 ( κ ln 11 H ) k s (35) Manning-Strickler: Cz = C 1/2 f = α r ( H k s ) 1/6 α r = 8.1, k s n k D s50 for flat bed, n k = 2. b. Normal flow: equilibrium state defined by a perfect balance between downstream gravitational impelling force and resisting bed force; flow is constant in time and in the downstream direction, i.e., uniform and steady. Conservation of mass and momentum reduce to depth-slope product rule for normal flow: τ b = ρghs, u = ghs (37) Combine momentum conservation and resistance forms: ρc f U 2 ghs = ρghs U = (38) C f With mass conservation: H = (36) ( ) Cf qw 2 1/3 (39) gs

7 GEO 432/532 COURSE NOTES 7 For Shields stress: τ = τ ( ) b ρrgd = Cf qw 2 1/3 S 2/3 g RD Also yields, with Manning-Strickler roughness, g U = α H 2/3 S 1/2 = H2/3 S 1/2 n k 1/6 s (40) (41) c. St. Venant shallow water equations: equations (2) and (3) d. Gradually varied flow and backwater curves Steady, gradually varied flow, St. Venant equations reduce to the backwater equation: dh dx = S S f 1 Fr 2, S = η x, S f = C f Fr 2 (42) Critical flow depth: Normal flow depth: H c = ( ) q 2 1/3 w (43) g H n = ( Cf q 2 w gs ) 1/3 or H n = (k 1/3 s q 2 w α 2 rgs ) 3/10 (44) If H n > H c, (Fr) n < 1: normal flow is subcritical, bed slope is mild. If H n < H c, (Fr) n > 1: normal flow is supercritical, bed slope is steep. For mild slope, three cases: M1: H 1 > H n > H c, flow downstream deepening into standing water, shallows and approaches normal flow upstream M2: H n > H 1 > H c, flow downstream approaches critical depth and waterfall, flow upstream deepens and approaches normal flow M3: H n > H c > H 1, flow downstream goes through hydraulic jump, shallows upstream Also family of steep curves. 6. Initiation of motion (Parker (2004), Ch. 6) a. Angle of repose: downslope gravitational force equals net Coulomb resistive force tan θ r = µ c (45) b. Critical Shields number, τ c : Shields number at threshold of motion from modified Shields diagram In limit of fully rough turbulent flow (large Re p ), τ c 0.03

8 8 LANCASTER c. Law of the wall for turbulent flow: u = 1 ( ) z u κ ln + B k s ( ) u k s ν where u k s /ν is the roughness Reynolds number: k s 8.5, > 8.62 or u k s > 100 (hydraulically rough) δ v ν B = ( ) κ ln u k s k s, < 0.26 or u k s < 3 (hydraulically smooth) ν δ v ν read off the graph, otherwise (transitional flow) Law of the wall for rough turbulent flow and balance of friction and drag forces produces estimate of critical Shields number near d. Significant streamwise slope: τ c = τ c0 cos α ( 1 tan α ) µ c e. Sediment mixtures: coarser grains are harder to move, but not as much harder as their size alone would dictate. Hiding function: τsci τscg = ( Di D sg (46) (47) (48) ) γ, 0.65 γ 0.90 (49) γ = 0: grain indepence, i.e., each size fraction behaves as if part of uniform bed γ = 1: equal threshold (or mobility) condition, i.e., hiding so effective that all grains move at the same critical boundary stress f. Condition for significant suspension: u sus = v s, or τ sus = v2 s RgD = R f 2 (Re p ) (50) g. Really a critical shear stress, or threshold of motion? Use of threshold may underpredict transport at low rates. 7. Bedload transport, mainly gravel on plane beds (Parker (2004), Ch. 7) a. Dimensionless Einstein bedload number: qb q b = (51) RgDD Seek transport as a function of Shields stress, q b = q b (τ ) or q b = q b (τ τ c ) b. Meyer-Peter and Müller bedload transport relation is typical of many (use corrected form of Wong, 2003): q b = 3.97 (τ τ c ) 3/2, τ c = (52) The main problem: data used for MP&M correspond to extreme flood flows in gravelbed rivers. Its threshold for initiation of motion is similar to τbf50 for most gravel-bed

9 GEO 432/532 COURSE NOTES 9 streams, so it predicts little or no transport of larger grain sizes, which do in fact move at bankfull, so MP&M relation, and those like it, are not the best relations for gravel-bed streams. c. Similarity of many relations motivates formulation of K (or W ): ( ) q b u 2 3/2 = K or K = Rgq b RgDD RgD u 3 (53) d. Sheet flow: τ < τsheet : Bedload localized to rolling, sliding, and saltating grains that exchange only with immediate bed surface. τ > τsheet : Bedload layer devolves into sliding layer of grains up to several grains thick. Estimates: τsheet , decreases with increasing Froude number. e. Mechanistic derivations: based on ξ bl : Volume of sediment in transport per unit area; u bl : Mean velocity of bedload particles; E bl : Volume rate of entrainment of bed particles into bedload movement (not suspension) per bed area; L sbl : Average step length of bedload particle, from entrainment to deposition, typically over many saltations. From continuity: q b = ξ bl u bl (cf. Bagnold), q b = E bl L sbl (cf. Einstein) (54) f. Calculations with bedload transport relations: Specify 1) sediment submerged gravity, R; 2) representative, or characteristic, surface grain size, e.g., D sg or D s50, for D in the relation; 3) value for shear velocity, u (and thus τ b ) from flow field with techniques of Ch. 5, e.g., ( ) u 2 k 1/3 3/10 ( ) s qw 2 = g 7/10 S 7/10 τ ks 1/3 3/10 qw 2 S 7/10 = (55) αrg 2 RD Compute α 2 r 1) τ = u2 RgD 2) q b = q b (τ ) 3) q b = RgDDq b

10 10 LANCASTER g. Alternative dimensionless bedload transport: W W is just K from above: W = q b τ = Rgq b 3/2 u 3 Most bedload transport relations can be expressed in terms of W h. Surface-based bedload transport formulation for mixtures: For active layer with thickness, L a = n a D s90, n a 1 2, divide bed material into N grain size ranges with characteristic size D i ; Let F i be fraction of material in surface (active) layer in i th size range; Volume of bedload transport per unit width in i th size range is q bi, total bedload is N q bt = q bi ; (57) i=1 Fraction of bedload in i th size range is p bi, where Define dimensionless terms for i th size range: (56) p bi = q bi q bt (58) τi = τ b = u2 ρrgd i RgD i (59) qbi q bi = RgDi D i F i (60) W i i. Parker bedload relation for gravel mixtures: Must renormalize F i s so that sand is removed. = Rgq bi u 3 F i (61) Based on surface geometric mean size, D sg (8), surface arithmetic standard deviation of ψ, σ s (9), (not geometric standard deviation, 2 σs ) computed from renormalized F i s, calculate transport in terms of a mobility function, G: W i = ( )G(φ i ) (62) ( ) Di φ i = ωφ sg0, φ sg0 = τ sg, τ D sg τssrg sg = u2, τssrg = (63) RgD sg ( ) 4.5, φ > 1.59 φ G(φ) = exp [ 14.2 (φ 1) 9.28 (φ 1) 2] (64), 1 φ 1.59 φ 14.2, φ < 1 σ s ω = 1 + σ 0 (φ sg0 ) [ω 0(φ sg0 ) 1] (65)

11 GEO 432/532 COURSE NOTES 11 ω 0 and σ 0 are set functions of φ sg0, computed numerically. The Parker transport relation is implemented as Acronym1 in the Excel workbook RTe-bookAcronym1.xls. (Note: ACRONYM isn t actually an acronym for anything; it s just a result of Gary s sense of humor.) There are three ways to run the model: Acronym1: User inputs specific gravity of sediment, R + 1, shear velocity, u, and the grain size distribution of the bed (surface) material. The code computes the magnitude and size distribution of the bedload transport. Acronym1 R: User inputs R + 1, flow discharge Q, bed slope S, channel width B, parameter n k relating roughness height k s to surface size D s90, and the grain size distribution of the bed (surface) material. The code computes shear velocity with Manning-Strickler resistance and assumption of normal flow and, then, magnitude and size distribution of the bedload transport. Acronym1 D: Uses the same formulation as Acronym1 R but allows specification of flow-duration curve so that average annual gravel transport rate and grain size distribution are computed. j. Wilcock and Crowe bedload relation for sand-gravel mixtures: Sand is not excluded; method based on D sg and fraction of sand in surface (active) layer, F s : Wi = G(φ i ) (66) 0.002φ, φ < 1.35 ( G(φ) = ) 4.5, φ > 1.35 φ 0.5 (67) φ i = τ sg ( ) b Di, τ τssrg ssrg = exp ( 20F s ), b = D sg exp (1.5 D i /D sg ) Note the strong dependence on sand fraction, F s, in the reference Shields stress, τ ssrg. Also, the exponent, b, on the hiding ratio, D i /D sg, is not constant but, rather, depends on the value of the ratio, and b , for small to large grain sizes. Wilcock & Crowe isn t set up to simply calculate magnitude and size distribution of transport, but the workbook, RTe-bookAgDegNormGravMixPW.xls, can be hacked to provide these outputs, as described in Exercise 2 for this class. k. Relative dimensionless shear stress, φ i : It s a bit tricky to sort out the meaning of the hiding ratio, D i /D sg, in the Parker and Wilcock & Crowe relations. First, note that b γ. So, what does the value of b signify relative to the extremes of equal mobility (γ = 1) and grain independence (γ = 0)? Let s break it down. For Wilcock & Crowe, φ i = τ sg τ ssrg ( Di D sg (68) ) b (69)

12 12 LANCASTER Note that, in the first ratio, both numerator and denominator are relative to the geometric mean grain size, D sg. Because we can always multiply by one, we can write φ i = τ sg ( ) γ ( ) 1 γ b D sg Di Di. (70) τssrg D i D sg D sg Then, from (49) and the definition of the Shields number, φ i = τ ( ) si 1 γ b Di (71) τssri D sg Okay, bear with me. If we suppose that a default value of γ is, say, γ = 0.6, then larger values of b effectively get us closer to equal mobility. If smaller grains have smaller γ, then they are more likely to be transported at lower bed shear stresses than are the particles with D sg. If larger grains have larger γ, then they are more likely to be transported at a bed shear stress more similar to that at which particles with D sg are transported. Note that we still haven t specified a value for γ, but we don t actually need to; we can still understand b in these terms: smaller b, less hiding of smaller grains and less protrusion of larger grains; larger b, more hiding of smaller grains and more protrusion of larger grains; if smaller grains have a smaller value of b than larger grains, that means smaller grains are effectively less hidden, and larger grains effectively more exposed. l. Annual bedload yield: Divide annual hydrograph into M bins, Specify characteristic discharge, Q k, for each bin, and Fraction of time, p Qk, the discharge occurs. Sum over discharges and, for mixtures, grain size fractions. 8. Bedforms (Parker (2004), Ch. 8) a. Ripples: b. Dunes: Small-scale bedforms Migrate downstream Asymmetric with gentle stoss and steep lee faces Require reasonably well-defined viscous sublayer, so slow flow, and ripples don t interact with water surface Criterion: Above threshold for motion, but grain diameter less than thickness of viscous sublayer, or D δ v = 11.6ν/u, or τ > τc and ( ) τ < (72) Re p Above criterion says, then, D 0.8 mm; in fact, usually find D 0.6 mm.

13 Common in sand-bed rivers GEO 432/532 COURSE NOTES 13 Wavelengths up to hundreds of meters, heights up to 5 m Asymmetric with gentle stoss and steep lee faces Characteristic of subcritical flow Migrate downstream Interact weakly with water surface, which is out of phase with dunes: flow accelerates over dune crests, where water surface is slightly lower. c. Antidunes: Can occur in sand-bed and gravel-bed rivers High (but not necessarily supercritical) Froude number Most commonly migrate upstream Nearly symmetric Interact strongly with water surface, which is in phase with antidunes: flow typically accelerates down face and into trough, where water surface is lower. Train of symmetrical surface (standing) waves indicative of antidunes on the bed d. Potential flow formulation, Froude number, and the dune-antidune transition Subcritical flow: water surface perturbation is approximately out of phase with bed perturbation Supercritical flow: water surface is approximately in phase with bed perturbation. Define Fr 0 in terms of the average flow depth, H 0, Fr 0 = q 2 w gh0 H 0 (73) and wavenumber of the bedforms, k = 2πH 0 /λ, which is like frequency. According to potential flow theory, the border between subcritical and supercritical flow is function of both Fr 0 and k, i.e., critical flow is defined by Fr 2 0 = tanh k (74) k For long-wavelength bedforms, k 0, and critical flow is at Fr 0 1, which is the shallow-water limit. For typical values of k, the critical flow transition is at Fr 0 < 1, e.g., for λ = 5H 0, k = 1.26, the threshold Fr 0 = Another threshold from potential flow theory defines a regime in which water surface and flow depth are out of phase: Fr 0 2 = 1 k tanh k So, we have three regions defined by Froude number and wavenumber (plus two planebed regimes): (75)

14 14 LANCASTER Fr 0 tanh k/k: Lower-regime plane bed, where subcritical flow is too low to create ripples or dunes; gravel-bed streams usually fall in this regime. Upper end of this regime is τ 2.72τ c. So, for coarse material (Re p 1), dunes will not form unless τ > , which is uncommon. Fr 0 < tanh k/k: Lower-regime, subcritical flow; water surface elevation and flow depth are in phase with one another, and both are out of phase with the bed elevation, and we have downstream-migrating dunes. tanh k/k Fr0 < 1/k tanh k: Upper-regime plane bed, near critical flow, where the competition between the effects of bedload, favoring dunes, and suspended load, favoring antidunes, cancel each other out. tanh k/k < Fr0 < 1/k tanh k: Upper-regime, supercritical flow; water surface elevation, flow depth, and bed elevation are all in phase, and we have upstreammigrating antidunes, although for very low k, cyclic steps may occur in this regime. 1/k tanh k < Fr0 : Upper-regime, supercritical flow; water surface and bed elevation are in phase, and both are out of phase with the flow depth, and we have downstream-migrating antidunes. Rare, but observed. 9. Form drag and sand transport (Parker (2004), Ch. 9) a. Skin friction and form drag: Skin friction: generated by the viscous shear stress acting tangentially to a body. Form drag: generated by the normal stress (mostly pressure) acting on a body. Total shear stress is τ ij = pδ ij + τ v,ij, τ v,ij = ρν ( ui + u ) j x j x i and the drag force on a body is given by integrating the shear stress times the outward normal vector over the surface: D i = τ ji n j ds (77) Drag force is the total of skin friction and form drag, or D i = D si + D fi. Skin friction is from the τ v,ji s, and form drag is from the pressure: ( ui D si = τ v,ji n j ds = ρν + u ) j n j ds (78) x j x i D fi = pn i ds (79) b. Einstein decomposition: subtract out form drag to perform bedload calculations. Given the two basic forms we have for bed shear stress, (76) τ b = ρc f U 2 = ρghs, (80)

15 GEO 432/532 COURSE NOTES 15 For flow with same mean velocity, U, and bed slope, S, with and without bedforms, the friction coefficient and depth due to skin friction alone must be smaller than with bedforms, or C f = C fs + C ff and H = H s + H f, so that τ b = τ bs + τ bf = ρ (C fs + C ff ) U 2 = ρg (H s + H f ) S (81) A decomposition like this is conceptually suspect, but this is generally the way it s done. Note that, while empirical formulations (e.g., shown below) exist for dunes in sand-bed rivers, methods for incorporating other sorts of forms, e.g., logs and boulders, are not well established, and finding the form drag due to a set of several or many such form drag elements is not generally a simple matter of adding together the contributions from each element. Perhaps the best treatment of this subject to date is that of Kean and Smith (2006a,b). c. Skin friction: calculate from grain size distribution and depth, with Keulegan (35) or Manning-Strickler (36). d. Form drag of dunes: Wright and Parker: skin friction as a function of total stress τs = ( ( τ Fr 0.7) ( ) ) HS U = , (82) gh which can be rewritten as [ H = Γ RD s50 S Use with skin friction predictor, C 1/2 fs = RD s50 ( ) ] g /13 ( ) τ 5/4, Γ = s 0.05 (83) U 0.7 U ghs S = 8.32 α strat ( Hs k s ) 1/6, k s = 3D 90, (84) where α strat is a correction for flow stratification, set to 1 without other information. This can be rewritten (without stratification) as U = 8.32 ( ) 1/6 Hs gh s S (85) Computation of depth-discharge relation proceeds as follows: Compute k s from D s90. Assume a value (or series of values) of H s. Compute velocity U from the skin friction predictor (85). Compute Γ and H from (83). Compute the water discharge per unit width, q w = UH. Plot H vs. q w. k s

16 16 LANCASTER e. Prediction of bedload transport with dunes: Substitute for τ with τ s in bedload transport relation, where for normal flow (note typo in Parker (2004)), τs = H ss (86) RD Proceed from depth-discharge relation, e.g., given discharge, slope, and grain size, find H s. There are other predictors, and we can generalize to gradually varied (i.e., not normal) flow. 10. Aggradation and degradation in alluvial streams a. Channel degradation and bank erosion (Parker (2004), Ch. 15) A 1D Exner equation can incorporate the effect of bank erosion on channel degradation. For bed lowering by an amount η, the contribution from the bed to the volume (per distance downstream) of sediment removal is (1 λ p ) B b η. Contribution from each sidewall is (1 λ p ) B s η = (1 λ p ) (η t η) η/s s. Volume rate of supply from both sidewalls per stream length is η I s = 2 (1 λ p ) B s t = 2 (1 λ p) η t η η (87) S s t The Exner equation for sediment concentration for a degrading channel with sidewall erosion and fraction, f bb, contributing to bed material load reduces to (1 λ p ) η t = 1 q t η t η (88) 1 + 2f x bb B b S s where the second term in the denominator of the first ratio on the right-hand side is the extra due to bank erosion. Note that, for a given flux divergence (the partial derivative on the right-hand side), sidewall erosion slows the rate of channel bed degradation. b. Aggradation and channel shift (Parker (2004), Ch. 15) Channel shift: Gradual shift is migration. Sudden shift is avulsion. Sinuosity, Ω: ratio between average down-channel distance x and average down-valley distance x v between the same points. For meandering rivers, typically 1.5 Ω 3. Assuming bed material load is transported within a channel of constant width B and sinuosity Ω, but is deposited across the floodplain width, B f, reach-averaged long-term Exner equation for aggradation is (1 λ p ) η t = Ω B q t B f x As long as Ω < B f /B, which is usually a pretty safe bet, then for a given flux divergence, floodplain deposition slows the rate of channel bed aggradation. Note that we re still only considering bed material load, so floodplain deposition here is really only from lateral accretion of point bars, etc., by a laterally migrating (meandering) channel. (89)

17 c. Longitudinal river profiles (Parker (2004), Ch. 25) Profile shapes: GEO 432/532 COURSE NOTES 17 A river profile is upward concave if slope decreases in the downstream direction, i.e., S x = η 2 x < 0 (90) 2 Otherwise, i.e., slope increases downstream, the profile is upward convex. Transient long profiles: If sediment inflow is greater than sediment outflow, i.e., if ( Q t / x) < 0, profile is upward-concave during aggradation toward a new equilibrium. If sediment inflow is less than sediment outflow, i.e., if ( Q t / x) > 0, profile is upward-convex during degradation toward a new equilibrium. Many rivers have upward-concave quasi-equilibrium long profiles ( quasi in the sense that sediment inflow and outflow are not equal). Reasons: Subsidence: filling hole left by subsidence makes ( Q t / x) < 0 over long term, so profile is upward-concave over long term. Sea level rise: Sea level rise of 120 m in last 12,000 yrs has same effect as subsidence. Delta progradation: forces long-term aggradation. Downstream sorting of sediment: Finer sediment is easier to transport, so same volume can be transported at a lower slope. Sediment is typically sorted such that median grain size decreases downstream. Finer sediment downstream lower slope downstream means upward-concave profile. Abrasion of sediment: Gravel is broken into finer particles, including sand, silt, and smaller gravel, all of which can be transported in the same volume at lower slopes, so effect is similar to downstream sorting. Effect of tributaries: Tributaries typically add proportionally more water than sediment supply (this is partly a way of saying that tributaries tend not to all provide high sediment concentrations at the same time), so sediment concentration decreases downstream and requires lower slopes for transport, and the tendency is for upward-concave profiles. For river with constant sea level rise: Reach length, L; floodplain width B f B; sand-bed grain size D; Transport Q tbf only during fraction I f of year, so annual yield is Q btf I f. Deposition over width of floodplain; Λ units of wash load (or floodplain material load) deposited per unit of bed material load deposited; no mass balance for wash load; assumed adequate supply.

18 18 LANCASTER Constant rate of sea level rise, ξ d ; constant sinuosity, Ω; normal flow. Exner (mass balance) equation: (1 λ p ) η t = ΩI f(1 + Λ) which differs from (89) only by the intermittency factor, I f, and the addition of the wash load ratio, Λ, in the numerator on the right-hand side. Relation for sediment transport in sand-bed stream as a function of channel-forming Shields number, τ form : B f Q t x (91) Q tbf = α EHτform Q RC 1/2 bf S (92) f A linear function of slope! Which means the Exner equation turns into a linear diffusion equation: η t = κ 2 η d x, κ 2 d = I fω(1 + Λ) α EH τform Q bf (1 λ p ) B f RC 1/2 f Decomposing (91) into constant sea level rise and deviations from the steady-state solution, we eliminate the deviations part by supposing that there is a steady-state solution to be found, and we get (the now rather simple ordinary differential equation): (93) which integrates to dq tbf dx = (1 λ p) B f I f Ω(1 + Λ) ξ d (94) Q tbf = Q tbf,feed (1 λ p) B f I f Ω(1 + Λ) ξ dx (95) From this solution, some interesting things can be shown: Slope decreases downstream, S = (1 βˆx)s u. Channel width decreases downstream, B bf Q tbf. Channel depth increases downstream, H bf Q 1 tbf. 11. Meandering and alluvial channel planform morphology (Lancaster and Bras, 2002) (see a. Meandering, or lateral bend migration, is a fundamental process, whatever the planform (e.g.,, combined with mid-channel bar formation and chute cut-off formation, makes braided river planform). b. Appears to require sufficient bedload to create point bars (Shepherd and Schumm, 1974). c. Point bars produce topographic steering, which transfers streamwise momentum laterally across the channel to the opposite bank.

19 GEO 432/532 COURSE NOTES 19 d. Transverse slope from point bar to pool, S T, can be explained in terms of helical flow, induced by centrifugal force, and is proportional to channel curvature, C and average flow depth, H: where, for gravel-bed channels (Ikeda, 1989), ( ) τ K = Cf S T = KHC (96) τ c (97) or, for sand-bed channels (Lancaster and Bras, 2002), K = τ sf τ [ ( sf τ ) τ τc κ ln sf H τ D s50 ] (98) e. Lateral transfer of downstream momentum means the high-velocity flow moves across the channel and impinges on the outer bank and creates a large lateral gradient of downstream velocity and, thus, a large shear stress on the bank. With distance downstream, the gradient and shear stress are dissipated by development of the rough turbulent boundary layer. f. General approach to modeling lateral migration: Changes in near-bank velocity proportional to channel curvature Effect on near-bank velocity at a point along the channel centerline, s, is proportional to a weighted average of local and upstream curvature: u b(s) = F 1 (U, b, C f ) C(s) + F 2 (U, H, b, C f ) s 0 C(s )e 2C f H (s s ) ds (99) Note that the model that simulates the 2D planform is essentially a 1D morphodynamic model. (Model of Lancaster and Bras (2002) is substantially different in its formulation of the shear stress on the bank but similar in integrating upstream effects of curvature.) Bank erosion modeled as proportional to u b (s); deposition on point bar assumed to keep pace with outer bank erosion. Different models generally similar in that Small bends migrate downstream Medium bends migrate outward Large bends migrate upstream Process model selects a characteristic scale for the meander wavelength (twice the distance from crossover to crossover)

20 20 LANCASTER References Ikeda, S. (1989), Sediment transport and sorting at bends, in River Meandering, vol. Water Resources Monograph 12, edited by S. Ikeda and G. Parker, pp , American Geophysical Union, Washington. Kean, J. W., and J. D. Smith (2006a), Form drag in rivers due to small-scale natural topographic features: 1. regular sequences, Journal of Geophysical Research, 111 (F4), doi: /2006JF Kean, J. W., and J. D. Smith (2006b), Form drag in rivers due to small-scale natural topographic features: 2. irregular sequences, Journal of Geophysical Research, 111 (F4), doi: /2006jf Lancaster, S. T., and R. L. Bras (2002), A simple model of river meandering and its comparison to natural channels, Hydrological Processes, 16 (1), Parker, G. (2004), 1D Sediment Transport Morphodynamics with Applications to Rivers and Turbidity Currents, University of Illinois, parkerg//morphodynamics_e-book.htm, visited April, Shepherd, R. G., and S. A. Schumm (1974), Experimental study of river incision, Geological Society of America Bulletin, 85 (2), Note that the list of references does not include many of those that are referenced in Parker (2004).

National Center for Earth-surface Dynamics: Renesse 2003: Non-cohesive Sediment Transport

National Center for Earth-surface Dynamics: Renesse 2003: Non-cohesive Sediment Transport National Center or Earth-surace Dynamics: Summary o Lectures on Transport o Non-Cohesive Sediment What is Morphodynamics? Sediment Properties Modes o Transport o Sediment Equations or Conservation o Bed

More information

SYLLABUS, GEO 432/532 APPLIED GEOMORPHOLOGY

SYLLABUS, GEO 432/532 APPLIED GEOMORPHOLOGY SYLLABUS, GEO 432/532 APPLIED GEOMORPHOLOGY Spring 2013 College of Earth, Ocean, and Atmospheric Sciences Oregon State University 3 credits T Th 8:00 9:20 am, Wlkn 210 Stephen Lancaster Wlkn 142, 7-9258,

More information

Sediment transport and river bed evolution

Sediment transport and river bed evolution 1 Chapter 1 Sediment transport and river bed evolution 1.1 What is the sediment transport? What is the river bed evolution? System of the interaction between flow and river beds Rivers transport a variety

More information

What? River response to base level rise. The morphodynamic system. Why? Channel-forming discharge. Flow. u = What s in a name. Flow Sediment transport

What? River response to base level rise. The morphodynamic system. Why? Channel-forming discharge. Flow. u = What s in a name. Flow Sediment transport River response to base level rise and other boundary conditions Dr. Maarten Kleinhans Summer course climate change and fluvial systems Course materials of Prof. Gary Parker Flow Sediment transport What?

More information

BED LOAD SEDIMENT TRANSPORT

BED LOAD SEDIMENT TRANSPORT BED LOAD SEDIMENT TRANSPORT Kamal EL KADI ABDERREZZAK EDF-R&D, Laboratoire National d Hydraulique et Environnement (LNHE) 1 17-19 September 2009 UNL, Santa Fe, Argentina OUTLINE I. Bed load II. Settling

More information

Aqueous and Aeolian Bedforms

Aqueous and Aeolian Bedforms Aqueous and Aeolian Bedforms 1 Further reading & review articles R.A. Bagnold, 1941, The physics of blown sand and desert dunes Charru et al., 2013, Sand ripples and dunes, Ann. Review of Fluid Mech. 2

More information

(3) Sediment Movement Classes of sediment transported

(3) Sediment Movement Classes of sediment transported (3) Sediment Movement Classes of sediment transported Dissolved load Suspended (and wash load ) Important for scouring algae Bedload (5-10% total load Moves along bed during floods Source of crushing for

More information

1.3.1.1 Incipient Motion Particle movement will occur when the instantaneous fluid force on a particle is just larger than the instantaneous resisting force related to the submerged particle weight and

More information

Sediment continuity: how to model sedimentary processes?

Sediment continuity: how to model sedimentary processes? Sediment continuity: how to model sedimentary processes? N.M. Vriend 1 Sediment transport The total sediment transport rate per unit width is a combination of bed load q b, suspended load q s and wash-load

More information

(3) Sediment Movement Classes of sediment transported

(3) Sediment Movement Classes of sediment transported 9/17/15 (3) Sediment Movement Classes of sediment transported Dissolved load Suspended load Important for scouring algae Bedload (5-10% total load) Moves along bed during floods Source of crushing for

More information

Università degli Studi di Napoli Federico II Facoltà di Ingegneria

Università degli Studi di Napoli Federico II Facoltà di Ingegneria Università degli Studi di Napoli Federico II Facoltà di Ingegneria Dottorato di Ricerca XX ciclo in Ingegneria dei Sistemi Idraulici, di Trasporto e Territoriali A procedure to store and access the stratigraphy

More information

Rivers T. Perron

Rivers T. Perron 1 Rivers T. Perron 12.001 After our discussions of large-scale topography, how we represent topography in maps, and how topography interacts with geologic structures, you should be frothing at the mouth

More information

PART 2:! FLUVIAL HYDRAULICS" HYDROEUROPE

PART 2:! FLUVIAL HYDRAULICS HYDROEUROPE PART 2:! FLUVIAL HYDRAULICS" HYDROEUROPE 2009 1 HYDROEUROPE 2009 2 About shear stress!! Extremely complex concept, can not be measured directly!! Computation is based on very primitive hypotheses that

More information

Geomorphology Geology 450/750 Spring Fluvial Processes Project Analysis of Redwood Creek Field Data Due Wednesday, May 26

Geomorphology Geology 450/750 Spring Fluvial Processes Project Analysis of Redwood Creek Field Data Due Wednesday, May 26 Geomorphology Geology 450/750 Spring 2004 Fluvial Processes Project Analysis of Redwood Creek Field Data Due Wednesday, May 26 This exercise is intended to give you experience using field data you collected

More information

Open Channel Flow Part 2. Ch 10 Young, notes, handouts

Open Channel Flow Part 2. Ch 10 Young, notes, handouts Open Channel Flow Part 2 Ch 10 Young, notes, handouts Uniform Channel Flow Many situations have a good approximation d(v,y,q)/dx=0 Uniform flow Look at extended Bernoulli equation Friction slope exactly

More information

Linear Analysis of Coupled Equations for Sediment Transport

Linear Analysis of Coupled Equations for Sediment Transport Theme B of the XXVII IAHR Congress, San Francisco, 1-15 August, 1997, 156-161. Linear Analysis of Coupled Equations for Sediment Transport YANTAO CUI and GARY PARKER St. Anthony Falls Laboratory, University

More information

15. Physics of Sediment Transport William Wilcock

15. Physics of Sediment Transport William Wilcock 15. Physics of Sediment Transport William Wilcock (based in part on lectures by Jeff Parsons) OCEAN/ESS 410 Lecture/Lab Learning Goals Know how sediments are characteried (sie and shape) Know the definitions

More information

EXAMPLES (SEDIMENT TRANSPORT) AUTUMN 2018

EXAMPLES (SEDIMENT TRANSPORT) AUTUMN 2018 EXAMPLES (SEDIMENT TRANSPORT) AUTUMN 2018 Q1. Using Cheng s formula estimate the settling velocity of a sand particle of diameter 1 mm in: (a) air; (b) water. Q2. Find the critical Shields parameter diameter

More information

NATURAL RIVER. Karima Attia Nile Research Institute

NATURAL RIVER. Karima Attia Nile Research Institute NATURAL RIVER CHARACTERISTICS Karima Attia Nile Research Institute NATURAL RIVER DEFINITION NATURAL RIVER DEFINITION Is natural stream of water that flows in channels with ih more or less defined banks.

More information

Calculation of Stream Discharge Required to Move Bed Material

Calculation of Stream Discharge Required to Move Bed Material Calculation of Stream Discharge Required to Move Bed Material Objective: Students will map two sections of a stream and calculate the depth, velocity, and discharge of flows required to move the stream

More information

Do you think sediment transport is a concern?

Do you think sediment transport is a concern? STREAM RESTORATION FRAMEWORK AND SEDIMENT TRANSPORT BASICS Pete Klingeman 1 What is Your Restoration Project Like? k? Do you think sediment transport is a concern? East Fork Lewis River, WA Tidal creek,

More information

Aquifer an underground zone or layer of sand, gravel, or porous rock that is saturated with water.

Aquifer an underground zone or layer of sand, gravel, or porous rock that is saturated with water. Aggradation raising of the streambed by deposition that occurs when the energy of the water flowing through a stream reach is insufficient to transport sediment conveyed from upstream. Alluvium a general

More information

Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.

Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative

More information

Module 2. The Science of Surface and Ground Water. Version 2 CE IIT, Kharagpur

Module 2. The Science of Surface and Ground Water. Version 2 CE IIT, Kharagpur Module The Science of Surface and Ground Water Lesson Sediment Dynamics in Alluvial Rivers and Channels Instructional Objectives On completion of this lesson, the student shall be able to learn the following:.

More information

INTRODUCTION TO SEDIMENT TRANSPORT AUTUMN 2018

INTRODUCTION TO SEDIMENT TRANSPORT AUTUMN 2018 INTRODUCTION TO SEDIMENT TRANSPORT AUTUMN 2018 1. OVERVIEW 1.1 Introduction 1.2 Particle properties 1.2.1 Diameter, d 1.2.2 Specific gravity, s 1.2.3 Settling velocity, w s 1.2.4 Porosity, P 1.2.5 Angle

More information

* Chapter 9 Sediment Transport Mechanics

* Chapter 9 Sediment Transport Mechanics Chapter 9 Sediment Transport Mechanics Section I Introduction 9-1. Definition Sedimentation embodies the processes of erosion, entrainment, transportation, deposition, and compaction of sediment. These

More information

Overview of fluvial and geotechnical processes for TMDL assessment

Overview of fluvial and geotechnical processes for TMDL assessment Overview of fluvial and geotechnical processes for TMDL assessment Christian F Lenhart, Assistant Prof, MSU Research Assoc., U of M Biosystems Engineering Fluvial processes in a glaciated landscape Martin

More information

MATHEMATICAL MODELING OF FLUVIAL SEDIMENT DELIVERY, NEKA RIVER, IRAN. S.E. Kermani H. Golmaee M.Z. Ahmadi

MATHEMATICAL MODELING OF FLUVIAL SEDIMENT DELIVERY, NEKA RIVER, IRAN. S.E. Kermani H. Golmaee M.Z. Ahmadi JOURNAL OF ENVIRONMENTAL HYDROLOGY The Electronic Journal of the International Association for Environmental Hydrology On the World Wide Web at http://www.hydroweb.com VOLUME 16 2008 MATHEMATICAL MODELING

More information

Stream Entrainment, Erosion, Transportation & Deposition

Stream Entrainment, Erosion, Transportation & Deposition Lecture 12 Zone 2 of the Fluvial System, Continued Stream Entrainment, Erosion, Transportation & Deposition Erosion in a Fluvial Landscape Corrosion Chemical Erosion Corrasion Mechanical Weathering Cavitation

More information

Water quality needs: Flow, velocity. Fish biologists need: Critical depth or velocity. Hydrology gives flows m 3 /s or day

Water quality needs: Flow, velocity. Fish biologists need: Critical depth or velocity. Hydrology gives flows m 3 /s or day Environmental Water Allocation Hydraulics Dr L Beevers Heriot Watt University, it UK l.beevers@hw.ac.uk Overview Why hydraulics in EWA? Different types of flows Theory Case studies Review Why do we need

More information

G433. Review of sedimentary structures. September 1 and 8, 2010

G433. Review of sedimentary structures. September 1 and 8, 2010 G433 Review of sedimentary structures September 1 and 8, 2010 Fluid Parameters The three main parameters that determine the stable bedform in unidirectional flow conditions are: grain size flow velocity

More information

Lecture Notes - Sediment Transport The flow problem Overview

Lecture Notes - Sediment Transport The flow problem Overview Lecture Notes - Sediment Transport The flow problem Overview In the last lecture, we emerged with a transport model in which the rate of transport q s depends on the shear stress τ. Recall that, for the

More information

Streams. Water. Hydrologic Cycle. Geol 104: Streams

Streams. Water. Hydrologic Cycle. Geol 104: Streams Streams Why study streams? Running water is the most important geologic agent in erosion, transportation and deposition of sediments. Water The unique physical and chemical properties of water make it

More information

7. Basics of Turbulent Flow Figure 1.

7. Basics of Turbulent Flow Figure 1. 1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds

More information

Sedimentation Scour Model Gengsheng Wei, James Brethour, Markus Grünzner and Jeff Burnham August 2014; Revised October 2014

Sedimentation Scour Model Gengsheng Wei, James Brethour, Markus Grünzner and Jeff Burnham August 2014; Revised October 2014 Flow Science Report 03-14 Sedimentation Scour Model Gengsheng Wei, James Brethour, Markus Grünzner and Jeff Burnham August 2014; Revised October 2014 1. Introduction The three-dimensional sediment scour

More information

Dedicated to Major Contributors to the Concepts of Flow of Water and Sediment in Alluvial Channels:

Dedicated to Major Contributors to the Concepts of Flow of Water and Sediment in Alluvial Channels: Geomorphic, Hydrologic, Hydraulic and Sediment Concepts Applied To Alluvial Rivers By Daryl B. Simons, Ph.D., P.E., D.B. Simons & Associates, Inc.; Everett V. Richardson, Ph.D., P.E., Ayres Associates,

More information

Technical Memorandum. To: From: Copies: Date: 10/19/2017. Subject: Project No.: Greg Laird, Courtney Moore. Kevin Pilgrim and Travis Stroth

Technical Memorandum. To: From: Copies: Date: 10/19/2017. Subject: Project No.: Greg Laird, Courtney Moore. Kevin Pilgrim and Travis Stroth Technical Memorandum To: From: Greg Laird, Courtney Moore Kevin Pilgrim and Travis Stroth 5777 Central Avenue Suite 228 Boulder, CO 80301 www.otak.com Copies: [Electronic submittal] Date: 10/19/2017 Subject:

More information

Channel Pattern. Channel Pattern, Meanders, and Confluences. Description of Channel Pattern. Bridge (2003)

Channel Pattern. Channel Pattern, Meanders, and Confluences. Description of Channel Pattern. Bridge (2003) Channel Pattern Channel Pattern, Meanders, and Confluences Outline Description of channel pattern Alternate bars Channel pattern continua and evolution Controls of channel pattern Description of Channel

More information

Final Report for TWDB Contract No

Final Report for TWDB Contract No Final Report for TWDB Contract No. 1004831127 Sediment Transport Modeling of Channel Scale Geomorphic Processes J.K. Haschenburger University of Texas at San Antonio July 31, 2012 1 Introduction This study

More information

ES 105 Surface Processes I. Hydrologic cycle A. Distribution % in oceans 2. >3% surface water a. +99% surface water in glaciers b.

ES 105 Surface Processes I. Hydrologic cycle A. Distribution % in oceans 2. >3% surface water a. +99% surface water in glaciers b. ES 105 Surface Processes I. Hydrologic cycle A. Distribution 1. +97% in oceans 2. >3% surface water a. +99% surface water in glaciers b. >1/3% liquid, fresh water in streams and lakes~1/10,000 of water

More information

LECTURE 2: FLUID MECHANICS

LECTURE 2: FLUID MECHANICS LECTURE 2: FLUID MECHANICS Introduction Conservation of mass and momentum General types of flow Laminar vs. turbulent flow Shear Stress Reach-average shear stress Bed roughness and reach average flow velocity

More information

3 Theoretical Basis for SAM.sed Calculations

3 Theoretical Basis for SAM.sed Calculations 3 Theoretical Basis for SAM.sed Calculations Purpose Sediment transport functions can be used to calculate the bed material portion of the sediment discharge rating curve. This rating curve can then be

More information

Fluvial Dynamics. M. I. Bursik ublearns.buffalo.edu October 26, Home Page. Title Page. Contents. Page 1 of 18. Go Back. Full Screen. Close.

Fluvial Dynamics. M. I. Bursik ublearns.buffalo.edu October 26, Home Page. Title Page. Contents. Page 1 of 18. Go Back. Full Screen. Close. Page 1 of 18 Fluvial Dynamics M. I. Bursik ublearns.buffalo.edu October 26, 2008 1. Fluvial Dynamics We want to understand a little of the basic physics of water flow and particle transport, as so much

More information

National Center for Earth-surface Dynamics: Renesse 2003: Non-cohesive Sediment Transport

National Center for Earth-surface Dynamics: Renesse 2003: Non-cohesive Sediment Transport Introduction to Morphodynamics For the original references on the work of Exner see Graf, W., 1971, Hydraulics of Sediment Transport, McGraw Hill, New York, 513 p. Sediment Properties Dietrich, E. W.,

More information

GLG362/GLG598 Geomorphology K. Whipple October, 2009 I. Characteristics of Alluvial Channels

GLG362/GLG598 Geomorphology K. Whipple October, 2009 I. Characteristics of Alluvial Channels I. Characteristics of Alluvial Channels Self-formed morphology set by entrainment, transport, and deposition They move unconsolidated sedimentary materials present in the valley fill flood plain/bank flow

More information

Erosion Surface Water. moving, transporting, and depositing sediment.

Erosion Surface Water. moving, transporting, and depositing sediment. + Erosion Surface Water moving, transporting, and depositing sediment. + Surface Water 2 Water from rainfall can hit Earth s surface and do a number of things: Slowly soak into the ground: Infiltration

More information

1.060 Engineering Mechanics II Spring Problem Set 8

1.060 Engineering Mechanics II Spring Problem Set 8 1.060 Engineering Mechanics II Spring 2006 Due on Monday, May 1st Problem Set 8 Important note: Please start a new sheet of paper for each problem in the problem set. Write the names of the group members

More information

Hydraulics Prof. Dr. Arup Kumar Sarma Department of Civil Engineering Indian Institute of Technology, Guwahati

Hydraulics Prof. Dr. Arup Kumar Sarma Department of Civil Engineering Indian Institute of Technology, Guwahati Hydraulics Prof. Dr. Arup Kumar Sarma Department of Civil Engineering Indian Institute of Technology, Guwahati Module No. # 04 Gradually Varied Flow Lecture No. # 07 Rapidly Varied Flow: Hydraulic Jump

More information

Prediction of bed form height in straight and meandering compound channels

Prediction of bed form height in straight and meandering compound channels Water Resources Management III 311 Prediction of bed form height in straight and meandering compound channels R. D. Karamisheva, J. F. Lyness, W. R. C. Myers, J. O Sullivan & J. B. C. Cassells School of

More information

Lecture 10: River Channels

Lecture 10: River Channels GEOG415 Lecture 10: River Channels 10-1 Importance of channel characteristics Prediction of flow was the sole purpose of hydrology, and still is a very important aspect of hydrology. - Water balance gives

More information

Fluvial Processes in River Engineering

Fluvial Processes in River Engineering Fluvial Processes in River Engineering Howard H. Chang San Diego State University... A WILEY-INTERSCIENCE PUBLTCATION John Wiley & Sons New York Chicbester Brisbane Toronto Singapore CONTENTS PARTI FLUVIAL

More information

Sand Ripple Dynamics on the Inner Shelf

Sand Ripple Dynamics on the Inner Shelf Sand Ripple Dynamics on the Inner Shelf Donald N. Slinn Department of Civil and Coastal Engineering, University of Florida Gainesville, FL 32611-6590, Phone: (352) 392-9537 x 1431 Fax: (352) 392-3466 E-mail:

More information

ESTIMATION OF MORPHOLOGICAL IMPACT OF GROYNE LENGTHENING I. RÁTKY, ÉVA RÁTKY

ESTIMATION OF MORPHOLOGICAL IMPACT OF GROYNE LENGTHENING I. RÁTKY, ÉVA RÁTKY ESTIMATION OF MORPHOLOGICAL IMPACT OF GROYNE LENGTHENING I. RÁTKY, ÉVA RÁTKY Abstract. Hydraulic-morphological calculations in open channel flows still cause problems for modellers, partially because of

More information

The Effect of Bedform-induced Spatial Acceleration on Turbulence and Sediment Transport

The Effect of Bedform-induced Spatial Acceleration on Turbulence and Sediment Transport The Effect of Bedform-induced Spatial Acceleration on Turbulence and Sediment Transport S. McLean (1) (1) Mechanical and Environmental Engineering Dept., University of California, Santa Barbara, CA 93106,

More information

Growing and decaying processes and resistance of sand waves in the vicinity of the Tone River mouth

Growing and decaying processes and resistance of sand waves in the vicinity of the Tone River mouth Advances in River Sediment Research Fukuoka et al. (eds) 2013 Taylor & Francis Group, London, ISBN 978-1-138-00062-9 Growing and decaying processes and resistance of sand waves in the vicinity of the Tone

More information

Advanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati

Advanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati Advanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati Module - 2 Uniform Flow Lecture - 1 Introduction to Uniform Flow Good morning everyone,

More information

MODELING THE TRANSPORT OF SAND AND MUD IN THE MINNESOTA RIVER CHUAN LI THESIS

MODELING THE TRANSPORT OF SAND AND MUD IN THE MINNESOTA RIVER CHUAN LI THESIS MODELING THE TRANSPORT OF SAND AND MUD IN THE MINNESOTA RIVER BY CHUAN LI THESIS Submitted in partial fulfillment of the requirements for the degree of Master of Science in Civil Engineering in the Graduate

More information

compare to Mannings equation

compare to Mannings equation 330 Fluid dynamics Density and viscosity help to control velocity and shear in fluids Density ρ (rho) of water is about 700 times greater than air (20 degrees C) Viscosity of water about 55 times greater

More information

Riparian Assessment. Steps in the right direction... Drainage Basin/Watershed: Start by Thinking Big. Riparian Assessment vs.

Riparian Assessment. Steps in the right direction... Drainage Basin/Watershed: Start by Thinking Big. Riparian Assessment vs. Riparian Assessment vs. Monitoring Riparian Assessment What is a healthy stream? Determine stream/riparian health Determine change or trend, especially in response to mgmt Classification = designation

More information

B-1. Attachment B-1. Evaluation of AdH Model Simplifications in Conowingo Reservoir Sediment Transport Modeling

B-1. Attachment B-1. Evaluation of AdH Model Simplifications in Conowingo Reservoir Sediment Transport Modeling Attachment B-1 Evaluation of AdH Model Simplifications in Conowingo Reservoir Sediment Transport Modeling 1 October 2012 Lower Susquehanna River Watershed Assessment Evaluation of AdH Model Simplifications

More information

MEANDER MIGRATION MODEL ASSESSMENT FOR THE JANUARY 2005 STORM, WHITMAN PROPERTY, SAN ANTONIO CREEK, VENTURA COUNTY, CALIFORNIA

MEANDER MIGRATION MODEL ASSESSMENT FOR THE JANUARY 2005 STORM, WHITMAN PROPERTY, SAN ANTONIO CREEK, VENTURA COUNTY, CALIFORNIA MEANDER MIGRATION MODEL ASSESSMENT FOR THE JANUARY 2005 STORM, WHITMAN PROPERTY, SAN ANTONIO CREEK, VENTURA COUNTY, CALIFORNIA Prepared by Eric Larsen, Ph.D. Mark Rains, Ph.D. October 2006 INTRODUCTION

More information

Diego Burgos. Geology 394. Advisors: Dr. Prestegaard. Phillip Goodling

Diego Burgos. Geology 394. Advisors: Dr. Prestegaard. Phillip Goodling Sediment Transport into an Urban Tributary Junction Diego Burgos Geology 394 Advisors: Dr. Prestegaard Phillip Goodling 1 Abstract Tributary junctions are an important component of stream morphology and

More information

Analysis of antidune migration direction

Analysis of antidune migration direction JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.109/010jf001761, 011 Analysis of antidune migration direction Francisco Núñez González 1 and Juan P. Martín Vide 1 Received 6 May 010; revised 7 January

More information

NATURE OF RIVERS B-1. Channel Function... ALLUVIAL FEATURES. ... to successfully carry sediment and water from the watershed. ...dissipate energy.

NATURE OF RIVERS B-1. Channel Function... ALLUVIAL FEATURES. ... to successfully carry sediment and water from the watershed. ...dissipate energy. 1 2 Function... Sevier River... to successfully carry sediment and water from the watershed....dissipate energy. 3 ALLUVIAL FEATURES 4 CHANNEL DIMENSION The purpose of a stream is to carry water and sediment

More information

Modelling of flow and sediment transport in rivers and freshwater deltas Peggy Zinke

Modelling of flow and sediment transport in rivers and freshwater deltas Peggy Zinke 1 Modelling of flow and sediment transport in rivers and freshwater deltas Peggy Zinke with contributions from Norwegian and international project partners 2 Outline 1. Introduction 2. Basic ideas of flow

More information

FUTURE MEANDER BEND MIGRATION AND FLOODPLAIN DEVELOPMENT PATTERNS NEAR RIVER MILES 241 TO 235, SACRAMENTO RIVER

FUTURE MEANDER BEND MIGRATION AND FLOODPLAIN DEVELOPMENT PATTERNS NEAR RIVER MILES 241 TO 235, SACRAMENTO RIVER FUTURE MEANDER BEND MIGRATION AND FLOODPLAIN DEVELOPMENT PATTERNS NEAR RIVER MILES 241 TO 235, SACRAMENTO RIVER Eric W. Larsen University of California, Davis With the assistance of Evan Girvetz REPORT

More information

EFFECT OF SEEPAGE-INDUCED NON-HYDROSTATIC PRESSURE DISTRIBUTION ON BEDLOAD TRANSPORT AND BED MORPHODYNAMICS

EFFECT OF SEEPAGE-INDUCED NON-HYDROSTATIC PRESSURE DISTRIBUTION ON BEDLOAD TRANSPORT AND BED MORPHODYNAMICS EFFECT OF SEEPAGE-INDUCED NON-HYDROSTATIC PRESSURE DISTRIBUTION ON BEDLOAD TRANSPORT AND BED MORPHODYNAMICS Simona Francalanci 1, Gary Parker 2 and Luca Solari 3 Abstract: The Shields number, or Shields

More information

Uniform Channel Flow Basic Concepts. Definition of Uniform Flow

Uniform Channel Flow Basic Concepts. Definition of Uniform Flow Uniform Channel Flow Basic Concepts Hydromechanics VVR090 Uniform occurs when: Definition of Uniform Flow 1. The depth, flow area, and velocity at every cross section is constant 2. The energy grade line,

More information

Towards the prediction of free-forming meander formation using 3D computational fluid dynamics

Towards the prediction of free-forming meander formation using 3D computational fluid dynamics Wasserbaukolloquium 2006: Strömungssimulation im Wasserbau 31 Dresdner Wasserbauliche Mitteilungen Heft 32 Towards the prediction of free-forming meander formation using 3D computational fluid dynamics

More information

21/08/2014. Introduction MSc course Fluvial Systems GEO Fluvial systems. Four scales. Four scales (1) Dr. Maarten G.

21/08/2014. Introduction MSc course Fluvial Systems GEO Fluvial systems. Four scales. Four scales (1) Dr. Maarten G. 21/08/2014 A B Faculty of Geosciences River and delta morphodynamics Introduction MSc course Fluvial Systems GEO4-4436 Dr. Maarten G. Kleinhans C D 2 Fluvial systems geomorphology sedimentology / geology

More information

Sediment Transport in Open Channels

Sediment Transport in Open Channels 35 Sediment Transport in Open Channels D. A. Lyn Purdue University 35.1 Introduction 35.2 The Characteristics of Sediment Density, Size, and Shape Size Distribution Fall (or Settling) Velocity Angle of

More information

Hydromechanics: Course Summary

Hydromechanics: Course Summary Hydromechanics: Course Summary Hydromechanics VVR090 Material Included; French: Chapters to 9 and 4 + Sample problems Vennard & Street: Chapters 8 + 3, and (part of it) Roberson & Crowe: Chapter Collection

More information

Flow and Bed Topography in a 180 Curved Channel

Flow and Bed Topography in a 180 Curved Channel Flow and Bed Topography in a 180 Curved Channel Jae Wook Jung 1, Sei Eui Yoon 2 Abstract The characteristics of flow and bed topography has been analyzed by changing the bed materials in a 180-degree,

More information

DOWNSTREAM SORTING OF SEDIMENT (additional control on channel width, depth and slope)

DOWNSTREAM SORTING OF SEDIMENT (additional control on channel width, depth and slope) DOWNSTREAM SORTING OF SEDIMENT (additional control on channel width, depth and slope) Image removed due to copyright concerns As the gravel gets finer, it is transported at lower slopes. The result is

More information

The importance of sediment supply data to modelling river morphodynamics

The importance of sediment supply data to modelling river morphodynamics The importance of sediment supply data to modelling river morphodynamics Wilbers & Ten Brinke (2003) bed elevation (m NAP) Pannerden channel Astrid Blom Delft University of Technology Netherlands Waal

More information

Geomorphic Importance of Winter Peak Flows and Annual Snowmelt Hydrographs in a Sierra Nevada Boulder-Bedrock River

Geomorphic Importance of Winter Peak Flows and Annual Snowmelt Hydrographs in a Sierra Nevada Boulder-Bedrock River Geomorphic Importance of Winter Peak Flows and Annual Snowmelt Hydrographs in a Sierra Nevada Boulder-Bedrock River Scott McBain and Bill Trush McBain & Trush, Inc. Clavey River and Cherry Creek vicinity

More information

Stream Classification

Stream Classification Stream Classification Why Classify Streams? Communication Tool Describe Existing Conditions & Trends Describe Restoration Goals Research Tool Morphologic Stream Classification Systems Schumm (1977) Alluvial

More information

In preparation for submission to Journal of Geophysical Research Earth Surface February, 2006

In preparation for submission to Journal of Geophysical Research Earth Surface February, 2006 QUASI-UNIVERSAL RELATIONS FOR BANKFULL HYDRAULIC GEOMETRY OF SINGLE-THREAD GRAVEL-BED RIVERS Gary Parker 1, Peter Wilcock 2, Chris Paola 3 and William E. Dietrich 4 1 Department of Civil & Environmental

More information

Landscape Development

Landscape Development Landscape Development Slopes Dominate Natural Landscapes Created by the interplay of tectonic and igneous activity and gradation Deformation and uplift Volcanic activity Agents of gradation Mass wasting

More information

Numerical Hydraulics

Numerical Hydraulics ETH Zurich, Fall 2017 Numerical Hydraulics Assignment 2 Numerical solution of shallow water wave propagation (www.surfertoday.com) 1 Introduction 1.1 Equations Understanding the propagation of shallow

More information

Department of Hydro Sciences, Institute for Urban Water Management. Urban Water

Department of Hydro Sciences, Institute for Urban Water Management. Urban Water Department of Hydro Sciences, Institute for Urban Water Management Urban Water 1 Global water aspects Introduction to urban water management 3 Basics for systems description 4 Water transport 5 Matter

More information

SEDIMENTATION AND ITS COUNTERMEASURE AT THE OFF-TAKE AREA OF NEW DHALESWARI RIVER

SEDIMENTATION AND ITS COUNTERMEASURE AT THE OFF-TAKE AREA OF NEW DHALESWARI RIVER SEDIMENTATION AND ITS COUNTERMEASURE AT THE OFF-TAKE AREA OF NEW DHALESWARI RIVER Tanjir Saif AHMED* MEE15634 Supervisors: Prof. EGASHIRA Shinji** Assoc. Prof. YOROZUYA Atsuhiro*** ABSTRACT Present study

More information

UNIFORM FLOW CRITICAL FLOW GRADUALLY VARIED FLOW

UNIFORM FLOW CRITICAL FLOW GRADUALLY VARIED FLOW UNIFORM FLOW CRITICAL FLOW GRADUALLY VARIED FLOW Derivation of uniform flow equation Dimensional analysis Computation of normal depth UNIFORM FLOW 1. Uniform flow is the flow condition obtained from a

More information

Submitted to Journal of Hydraulic Research, April 18, 2003 MODELING DOWNSTREAM FINING IN SAND-BED RIVERS II: APPLICATION

Submitted to Journal of Hydraulic Research, April 18, 2003 MODELING DOWNSTREAM FINING IN SAND-BED RIVERS II: APPLICATION Submitted to Journal of Hydraulic Research, April 18, 23 MODELING DOWNSTREAM FINING IN SAND-BED RIVERS II: APPLICATION Scott Wright Graduate Student, St. Anthony Falls Laboratory, University of Minnesota

More information

2.3 The Turbulent Flat Plate Boundary Layer

2.3 The Turbulent Flat Plate Boundary Layer Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The

More information

COMPARISON OF TRANSPORT AND FRICTION OF MONO- SIZED AND TWO-SPECIES SEDIMENT IN UPPER PLANE BED REGIME

COMPARISON OF TRANSPORT AND FRICTION OF MONO- SIZED AND TWO-SPECIES SEDIMENT IN UPPER PLANE BED REGIME ISBN 978-83-927084-8-3 ISSN 0867-7964 COMPARISON OF TRANSPORT AND FRICTION OF MONO- SIZED AND TWO-SPECIES SEDIMENT IN UPPER PLANE BED REGIME Štěpán Zrostlík, Vojtěch Bareš, Jan Krupička, Tomáš Picek, Václav

More information

ADDRESSING GEOMORPHIC AND HYDRAULIC CONTROLS IN OFF-CHANNEL HABITAT DESIGN

ADDRESSING GEOMORPHIC AND HYDRAULIC CONTROLS IN OFF-CHANNEL HABITAT DESIGN ADDRESSING GEOMORPHIC AND HYDRAULIC CONTROLS IN OFF-CHANNEL HABITAT DESIGN Conor Shea - Hydrologist U.S. Fish and Wildlife Service Conservation Partnerships Program Arcata, CA Learning Objectives Examine

More information

Streams. Stream Water Flow

Streams. Stream Water Flow CHAPTER 14 OUTLINE Streams: Transport to the Oceans Does not contain complete lecture notes. To be used to help organize lecture notes and home/test studies. Streams Streams are the major geological agents

More information

Factors affecting confluence scour

Factors affecting confluence scour & Wang (eds) River Sedimentation 1999., Balkema, Rotterdam. ISBN 9 9 3. 17 19 Factors affecting confluence scour R. B. Rezaur & A. W. Jayawardena. Department of Civil Engineering, The University of Hong

More information

EXAMPLE SHEET FOR TOPIC 3 AUTUMN 2013

EXAMPLE SHEET FOR TOPIC 3 AUTUMN 2013 EXAMPLE SHEET FOR TOPIC ATMN 01 Q1. se dimensional analysis to investigate how the capillary rise h of a liquid in a tube varies with tube diameter d, gravity g, fluid density ρ, surface tension σ and

More information

U.S. Army Corps of Engineers Detroit District. Sediment Trap Assessment Saginaw River, Michigan

U.S. Army Corps of Engineers Detroit District. Sediment Trap Assessment Saginaw River, Michigan U.S. Army Corps of Engineers Detroit District December 2001 December 2001 This report has been prepared for USACE, Detroit District by: W.F. BAIRD & ASSOCIATES LTD. 2981 YARMOUTH GREENWAY MADISON, WISCONSIN

More information

Erosion Rate is a Function of Erodibility and Excess Shear Stress = k ( o - c ) From Relation between Shear Stress and Erosion We Calculate c and

Erosion Rate is a Function of Erodibility and Excess Shear Stress = k ( o - c ) From Relation between Shear Stress and Erosion We Calculate c and Equilibrium, Shear Stress, Stream Power and Trends of Vertical Adjustment Andrew Simon USDA-ARS, Oxford, MS asimon@msa-oxford.ars.usda.gov Non-Cohesive versus Cohesive Materials Non-cohesive: sands and

More information

The Equilibrium Channel & Channel Change. Peter Wilcock 3 August 2016

The Equilibrium Channel & Channel Change. Peter Wilcock 3 August 2016 The Equilibrium Channel & Channel Change Peter Wilcock 3 August 2016 1 The search for common empirical attributes of streams Luna B. Leopold M. Gordon ( Reds ) Wolman Watts Branch, MD A meandering stream

More information

Morphology and Sedimentology of Panther Creek, Montgomery County Preserve

Morphology and Sedimentology of Panther Creek, Montgomery County Preserve Morphology and Sedimentology of Panther Creek, Montgomery County Preserve, Ph.D. Houston, Texas February, 2014 This guide is for educational use. The author gives permission to reproduce this guide for

More information

Figure 34: Coordinate system for the flow in open channels.

Figure 34: Coordinate system for the flow in open channels. OE466 redging Processes 5. SCOUR 5.. Steady uniform flow in open channels This chapter is written with a view to bottom scour. The main outcome is the scour velocity as a function of the particle diameter.

More information

BZ471, Steam Biology & Ecology Exam 1

BZ471, Steam Biology & Ecology Exam 1 BZ471, Exam1, p.1 BZ471, Steam Biology & Ecology Exam 1 Name Matching. Select the single best (most precise) answer for each of the following. Use each answer no more than once. No partial credit is given.

More information

Uniform Channel Flow Basic Concepts Hydromechanics VVR090

Uniform Channel Flow Basic Concepts Hydromechanics VVR090 Uniform Channel Flow Basic Concepts Hydromechanics VVR090 ppt by Magnus Larson; revised by Rolf L Feb 2014 SYNOPSIS 1. Definition of Uniform Flow 2. Momentum Equation for Uniform Flow 3. Resistance equations

More information

Sediment Transport, Numerical Modeling and Reservoir Management some Concepts and Applications

Sediment Transport, Numerical Modeling and Reservoir Management some Concepts and Applications Sediment Transport, Numerical Modeling and Reservoir Management some Concepts and Applications CEMRACS 2013 August 6 th Magali Jodeau EDF R&D LNHE magali.jodeau@edf.fr Overview of the presentation What

More information

Lectures Hydrology & Fluvial Geomorphology. Gauley River Images. Ancients' (= Biblical) Model of Water (Hydrologic) Cycle

Lectures Hydrology & Fluvial Geomorphology. Gauley River Images. Ancients' (= Biblical) Model of Water (Hydrologic) Cycle Lectures 11-13 13 Hydrology & Fluvial Geomorphology Gauley River Images http://www.youtube.com/watch?v=eulmuyegtz4&feature=related Ancients' (= Biblical) Model of Water (Hydrologic) Cycle Stream Water

More information

LAB-SCALE INVESTIGATION ONBAR FORMATION COORDINATES IN RIVER BASED ON FLOW AND SEDIMENT

LAB-SCALE INVESTIGATION ONBAR FORMATION COORDINATES IN RIVER BASED ON FLOW AND SEDIMENT LAB-SCALE INVESTIGATION ONBAR FORMATION COORDINATES IN RIVER BASED ON FLOW AND SEDIMENT Mat Salleh M. Z., Ariffin J., Mohd-Noor M. F. and Yusof N. A. U. Faculty of Civil Engineering, University Technology

More information