A VISUAL MODEL FOR PREDICTING STREAM RESPONSE OF ALLUVIAL GRAVEL-BED RIVERS by Nerissa Yarmila Moscote B.Sc.E., Queen's University, Kingston, 1997 A thesis submitted in partial fulfillment of the requirement for the degree of Master of Applied Science in The Faculty of Graduate Studies Department of Civil Engineering We accept this thesis as conforming to the required standard T H E U N I V E R S I T Y OF BRITISH C O L U M B I A December, 2002 © Nerissa Y l Moscote, 2002 In presenting this thesis in partial fulfilment of the requirements degree at the University for an advanced of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that publication of this thesis for financial gain shall not be allowed without permission. Department of CAV) I £i/ng IhejPJKl^^ The University of British Columbia Vancouver, Canada Date JCASIA DE-6 (2/88) extensive .Xll ; 3003 copying or my written ABSTRACT In this thesis, a visual optimization model, RiverMod, has been developed to predict the channel geometry of alluvial gravel-bed rivers. The model is an extension of a previous model developed by Millar and Quick (1993b). The model is based on equilibrium theory, and also includes a bank stability analysis. The adjustment of the dependent variables of a river reach is quantified according to changes incurred by the independent variables. The primary dependent variables are channel width, depth and slope, and the primary independent variables are the discharge, flow resistance parameters, sediment size distribution and bank stability parameters. RiverMod, written in Visual Basic, provides a user interface that promotes clear and simple model usage for multiple runs with different sediment transport and flow resistance equations. The theory behind the Millar and Quick (1993b) model is discussed, as well as Millar's (2000) meandering-braiding transition and its incorporation into RiverMod. Four sediment transport equations and four flow resistance equations are included as part of the model. The model consists of fixed-channel-slope and variable-channel slope versions. The fixed slope version is equivalent to an experiment where the slope is fixed, and the channel width, depth and sediment transport rate adjust to the discharge. The variable slope version more closely approximates natural stream conditions. Each of these models can be applied to streams with cohesive or non-cohesive bank sediment. The model is applied to two gravel-bed rivers, located in British Columbia. The bank stability of both rivers has been decreased due to logging along the banks. RiverMod is used to quantify the impact of such a disturbance by analyzing past and present channel conditions. Restoration methods are suggested based on model output. This thesis provides a full description of theory underlying the model development and a description of model usage; therefore, this thesis is a complete user manual for the RiverMod optimization model. iii Table of Contents Abstract ii Table of Contents iv List of Tables viii List of Figures ix List of Symbols xii Acknowledgments xiii 1 INTRODUCTION 1.1 INTRODUCTION 1 1.2 USES OF THE MODEL 2 1.3 EQUILIBRIUM 3 1.4 TEMPORAL AND SPATIAL SCALES 3 1.5 INDEPENDENT AND DEPENDENT VARIABLES 4 1.5.1 6 1.6 Adjustment of the Dependent Variables STREAM IMPACTS AND CHANNEL RESPONSES 7 1.6.1 Urbanization 8 1.6.2 Logging and Agriculture 10 1.6.3 Dams and Reservoirs 11 1.6.4 Mineral Extraction 12 iv 2 MODEL FORMULATION 2.1 2.2 2.3 2.4 2.5 THEORETICAL BASIS 21 2.1.1 Extremal Hypothesis 22 2.1.2 Constraints 22 SEDIMENT TRANSPORT 23 2.2.1 Einstein (1942), Einstein-Brown (1950) 25 2.2.2 Meyer-Peter and Miiller (1948) 27 2.2.3 Ackers 2.2.4 Parker, Klingeman and McLean (1982) 32 2.2.5 Discussion of Sediment Transport Equations 34 and White (1973) 29 FLOW RESISTANCE 35 2.3.1 Manning (1891) 36 2.3.2 Composite Manning (Horton, 1933) 38 2.3.3 Keulegan (1938) 38 2.3.4 Jarrett (1984) 40 BANK STABILITY 41 2.4.1 Non-Cohesive Bank Sediment 42 2.4.2 Cohesive Bank Sediment 44 2.4.3 Effect of Vegetation on Bank Stability 46 PLANFORM GEOMETRY 47 v 3 TYPE OF ANALYSIS 3.1 INTRODUCTION 58 3.2 FIXED-CHANNEL-SLOPE MODEL 58 3.3 VARIABLE-CHANNEL-SLOPE MODEL 59 4 USING RIVERMOD 4.1 INSTALLATION 64 4.2 FILE SETUP SCREEN 64 4.2.1 Open an Existing File 64 4.2.2 Creating a New File 65 4.3 MODEL SCREEN 65 4.4 E Q U A T I O N SCREEN 66 4.5 C O M P U T A T I O N SCREEN 66 4.5.1 Model Input 66 4.5.2 Model Output 68 4.5.3 Drop-Down Menu 70 4.5.3.1 File 70 4.5.3.2 Edit 71 4.5.3.3 Output 72 4.5.3.4 M i c r o s o f t Excel 74 vi 4.6 PLANFORM SCREEN 76 5 SAMPLE APPLICATIONS 5.1 INTRODUCTION 88 5.2 SLESSE CREEK 88 5.2.1 Investigation 89 5.2.2 Analysis 91 5.2.3 Results 92 NARROWLAKE CREEK 94 5.3.1 Investigation 95 5.3.2 Analysis 97 5.3.3 Results 98 5.3 References 111 vii List of Tables Table 2.1 Bedload transport formulae available in RIVERMOD (adapted from Gomez and Church, 1989) 34 Table 2.2 Values of Manning's roughness coefficient n for natural streams (Chow, 1959) 36 Table 2.3 Summary of "best-fit" k values derived from reach-averaged hydraulic geometry (Millar, 1999) 40 Table 3.1 Results of Millar (1994) active bank stability constraint investigation 61 Table 4.1 Minimum and maximum values for input variables 69 Table 5.1 Slesse Creek channel geometry (MacVicar, 1999) 89 Table 5.2 Input variables for Slesse Creek (adapted from MacVicar, 1999) 91 Table 5.3 Narrowlake Creek channel geometry 95 Table 5.4 Input variables for Narrowlake Creek (adapted from Wilson, 2001) 96 s viii List of Figures Figure 1.1 Schematic representation of the graded condition in two time scales 13 Figure 1.2 Statistical model of a river system illustrating the influence of the independent variables on the dependent variables (Hey, 1976 as cited by Hey, 1982) 14 Figure 1.3 A schematic representation of two solution curves showing the optimal 15 geometry (Millar, 1994) Figure 1.4 Illustration of a dependent variable vs. time. 16 Figure 1.5 Idealized watershed showing three zones of erosion, transport and deposition (Kondolf, 1997) 17 Figure 1.6 Drainage net of Rock Creek in 1913, before modern urbanization and in 18 1964 after modern urbanization (Dunne and Leopold, 1978) Figure 1.7 Channel cross-sectional area at bankf ull vs. drainage area of both rural 19 and urban streams in Pennsylvania (Dunne and Leopold, 1978) Figure 1.8 Reduction in sediment supply from the catchment of the San Luis Rey River due to the construction of the Henshaw Dam (Kondolf, 1997) 20 Figure 2.1 Definition Sketch of a simplified trapezoidal channel for Millar and Quick (1993) model (MacVicar 1999) 50 Figure 2.2 A schematic representation of the armour layer in a gravel-bed channel 51 Figure 2.3 Comparison between observed and calculated data for Meyer-Peter and Muller, Ackers and White, Einstein, and Parker, Klingeman and McLean formulae and Elbow River data (Gomez and Church, 1989) 52 Figure 2.4 Relation between grain-size diameter and Manning's roughness coefficient for channels in which the bottom and sides have different roughness characteristics (Horton, 1933) 53 Figure 2.5 Relation of Manning's Roughness Coefficient to Hydraulic Radius (Jarrett, 1984) 54 IX Figure 2.6 Rivers of increasing sinuosity from straight to meandering (Thorne, 1997) 55 Figure 2.7 Flowchart for the Planform Model 56 Figure 2.8 Definition sketch of one meander arc length (Millar, 1997; Sediment Transport (CIVL 546) lecture notes - University of British Columbia) 57 Figure 3.1 Flowchart for Millar and Quick (1993b) Fixed-Slope Model (MacVicar, 1999) 62 Figure 3.2 Flowchart for Millar and Quick (1993b) Variable-Slope Model (MacVicar, 1999) 63 Figure 4.1 File Setup Screen 78 Figure 4.2 File Setup Screen with selected file name 79 Figure 4.3 Model Screen 80 Figure 4.4 Equation Screen 81 Figure 4.5 Computation Screen for Non-Cohesive Sediment with Einstein-Brown and Keulegan equations in SI Units, using the Fixed Slope Model 82 Figure 4.6 Computation Screen for Cohesive Sediment with Ackers and White and Composite Manning equations in SI Units using the Fixed Slope Model 83 Figure 4.7 Supplementary Output Screen for (a) Non-Cohesive Sediment Model; (b) Cohesive Sediment Model 84 Figure 4.8 Search Screen showing Cohesive Sediment Model Options 85 Figure 4.9 Planform Screen 86 Figure 4.10 Microsoft Excel Spreadsheet with values appended from RIVERMOD 87 Figure 5.1 Location of Slesse Creek Watershed (MacVicar, 1999) 100 Figure 5.2 Airphoto of Slesse Creek (MacVicar, 1999) 101 x Figure 5.3 Slesse Creek lower watershed, showing study reach D in 1936 (MacVicar, 1999) 102 Figure 5.4 Slesse Creek calibration, width (W) vs. modified friction angle (<()') 103 Figure 5.5 Restoration of Slesse Creek, altering <t>' (Dsobank = 0.133 m) 104 Figure 5.6 Restoration of Slesse Creek, altering Dsobank 0t>' = 4 0 ° ) 105 Figure 5.7 Location of Narrowlake Creek Watershed 106 Figure 5.8 Airphoto of Narrowlake Creek (Wilson, 2001) 107 Figure 5.9 Narrowlake Creek calibration, width (W) vs. modified friction angle (<(>') 108 Figure 5.10 Restoration of Narrowlake Creek, altering (j)' (Dsobank 0.047 m) 109 Figure 5.11 Restoration of Narrowlake Creek, altering Dsobank ( f = 4 0 ° ) 110 = xi List of Symbols A - A = function of grain diameter (Ackers and White) a = empirical coefficient C - c = soil cohesion C cross-sectional area function of grain diameter (Ackers and White) = constant (Keulegan) x d = D = sediment grain size D35 = thirty-fifth percentile bed surface grain diameter O50 = median bed surface grain diameter dso = median sub-surface grain diameter C>50bank depth median bank grain diameter of non-cohesive sediment = D90 = ninetieth percentile bed surface grain diameter D a = average bed surface grain diameter D gr = grain diameter (Ackers and White) D x = characteristic grain diameter, where x is % of sediment finer F = fall parameter (Einstein-Brown) f = F gr Darcy-Weisbach friction factor = sediment mobility number (Ackers and White) xii F> = FSH = factor of safety for mass failure FS T Froude number = factor of safety for fluvial erosion g = gravitational acceleration (assumed = 9.81 m/s ) <5b = sediment transport capacity gb* = dimensionless transport rate &bca\c - sediment transport capacity Gbf = bankfull sediment transport capacity G = sediment transport (Ackers and White) 2 gr H = vertical bank height H = Hcrit = maximum stable value of H, I = i'b = specific bed load transport rate (dry weight) (Gomez and Church, 1989) KB = component of the reducing multiplier (Meyer-Peter and Muller) KG = component of the reducing multiplier (Meyer-Peter and Muller) k = equivalent roughness ' = distance along valley axis L = m = function of grain diameter (Ackers and White) n = function of grain diameter (Ackers and White) s Hammer number input of sediment at upstream end of channel distance along meander arc length xiii n = Manning's roughness coefficient tibank = Manning's roughness coefficient for the bank sediment nbed = Manning's roughness coefficient for the bed sediment N = dimensionless stability number s O = output of sediment at downstream end of channel Pbank = bank perimeter Pbed = wetted perimeter of bed Q = discharge Qbf = bankfull discharge q s = sediment discharge c = radius of curvature r Rh = hydraulic radius S = sediment storage S = channel slope S* = transition slope SF = shear force s = specific gravity of sediment (assumed = 2.65) s S v = valley slope U* = shear velocity U,u = mean velocity W = width of channel xiv X = sediment f l u x (Ackers and W h i t e ) Y = depth of channel A = change O = maximum deviation angle along channel G>* = transport parameter <j>' = modified f r i c t i o n angle <(>5o = shear s t r e s s ratio (Parker, Klingeman and McLean) y = unit weight of water = unit weight of sediment Yt = saturated unit weight of soil A, = meander wavelength u = dynamic viscosity v = kinematic viscosity 6 = bank angle Pf = fluid density p = sediment density x = shear stress x* = dimensionless bed shear s t r e s s tbank = mean bank shear s t r e s s Tbed = mean bed shear s t r e s s tent = critical shear s t r e s s f o r fluvial erosion of bank sediment y s s (Einstein-Brown) xv Tdso = co* = dimensionless bedload transport rate (Parker, Klingeman and McLean) £ = sinuosity ^ = sinuosity calculated using the sine-generated curve c v|/ = bed shear stress for the median sub-surface grain size (Parker, Klingeman and McLean) ratio of forces acting on a particle xvi ACKNOWLEDGEMENTS I would like to thank my supervisor Dr. Robert Millar for his contribution to this thesis. Many other people at the University of British Columbia provided valuable assistance. Edmond Yu answered all my questions about Visual Basic, Andrew Wilson provided all sorts of data concerning Narrowlake Creek and Aleteia Greenwood helped me find all the references used in this thesis. I am grateful to Rafael Ayax Moscote, Norma Moscote and Diana Lavery for their constant support. I am grateful to Mary Ellen Mulligan for her support and excellent editing. Thank you Waverley Mei for working on your thesis while I worked on mine. Thank you Reimy for being born at just the right time. Thank you Ryan for the juice, for believing and so much more. —xiii-- Chapter 1 INTRODUCTION 1.1 INTRODUCTION This project is an extension of a hydraulic geometry model developed by Millar and Quick in 1991, 1993 and 1998. Millar developed "a static analytical model based on previous optimization formulations (Millar (1991), Chang (1980), White et al. (1982))". The model is an optimization model which incorporates the effects of bank-stability on the hydraulic geometry of alluvial gravel-bed rivers. The original model, written using Basic, is not very flexible, is difficult to use and is therefore used primarily as a research tool. The computer program described in this thesis, named RiverMod, was written using Microsoft Visual Basic. I t combines four separate, related models and includes several new features. I t is visual, user-friendly and includes: (i) Four sediment transport relations (ii) Four flow resistance equations (iii) A planform model (iv) A user-interface for input and output of data (v) Compatibility with Microsoft Excel (vi) The ability to save, rewrite and group files 1 RiverMod applies to Single-thread alluvial gravel-bed rivers with mobile beds that have the capacity to modify their channel dimensions (Millar, 1994). The assumption that the model treats the river as an equilibrium system forms the basis for the optimization model. 1.2 USES OF THE MODEL The model has two primary applications: I t has proven to be both an effective teaching and research tool. The following is a list of possible uses: (a) River adjustments: RiverMod can help interpret river adjustments by helping researchers and students gain a further understanding into the behaviour of rivers and also predict future changes. (b) Activities that impact rivers: The model can aid in monitoring the effects of disturbances such as urbanization, agricultural development, the construction of dams and reservoirs and logging on river channel geometry. The effects of natural processes such as wild fire or climate change can also be studied. (c) Sensitivity Analysis: The model can be used to run a sensitivity analysis in order to determine how sensitive a variable is to a particular activity. (d) Calibration: The program can also be calibrated to the observed geometry. This can assist in the determination of a variable which is hard to measure with field evaluations. 2 1.3 EQUILIBRIUM RiverMod analysis can be applied to a reach of a river considered to be in equilibrium. For t h e purposes o f this study, a channel will be considered to be in equilibrium when t h e following conditions are met (Millar, 1994): (i) " T h e mean hydraulic geometry of t h e channel reach remains unchanged over an appropriate time scale f o r which a steady-state equilibrium can be assumed. (ii) T h e r e is no net erosion or deposition along t h e reach. (iii) Any perturbations from equilibrium geometry will be o f f s e t and t h e equilibrium geometry restored." 1.4 TEMPORAL A N D SPATIAL SCALES A reach of a river can be said to be in equilibrium when continuity is conserved during specific temporal and spatial scales. T h r e e principal time scales should be considered: Engineering time, which usually spans between 100 and 200 years; geomorphic, or graded time, which is usually less than 10,000 years and; geologic time, which is usually greaier than 10,000 years. For example, over a geologic time frame, t h e slope may appear to be constantly changing i.e. adjusting toward equilibrium. But in t h e short-term, t h e slope will appear to have reached a state of approximate equilibrium (see Figure 1.1). 3 Adjustments needed to achieve grade usually occur over several decades or over hundreds o f years. RiverMod analysis applies to engineering time scales over which a reach or a river can be considered to be in an approximate equilibrium. Equilibrium o f a specific system depends on t h e spatial scale. Generally, an entire river will not be in equilibrium, but certain reaches of t h e river may be. T h e t e r m "reach of a river" is usually used because d i f f e r e n t reaches o f t h e same river may exhibit d i f f e r e n t channel patterns and t h e r e f o r e d i f f e r e n t ultimate graded conditions at t h e same point in time. 1.5 INDEPENDENT AND DEPENDENT VARIABLES I n order to study the way in which rivers adjust their f o r m and dimensions, the variables defining t h e hydraulic geometry o f the channel must be identified and t h e independent and dependent variables must be specified (Hey, 1982). An independent variable is imposed externally on the system; a dependent variable is determined by the value of t h e independent variables (Vanoni, 1975). I^ac\/\car (1999) defines t h e independent variables as those imposed on t h e reach, namely: climate, geology, runoff, vegetation type, vegetation density, and relief. Climate, geology and runoff determine t h e discharge, sediment size and sediment yield. Bank stability is controlled by t h e vegetation type and density and relief determines t h e valley slope. A t the decade time scale, the independent variables become: discharge o f water, 4 sediment discharge, sediment size, bank stability (bank vegetation and soil properties) and valley slope. The dependent variables are hydraulic geometry variables which respond to changes in the imposed independent variables. The primary dependent variables can be considered to be width, depth and bed slope. Other secondary dependent variables, which can be determined from the primary ones are velocity, sinuosity, meander arc length, maximum depth and height of bedforms (MacVicar, 1999). Figure 1.2 is a statistical model of a river system which illustrates the relationship between the dependent and independent variables (Hey, 1976 as cited by Hey, 1982). Any significant change in the independent variables will produce a new regime channel geometry which corresponds to the modified values of the controlling variables. When the channel is again stable, it will be defined by the new independent variables (Hey and Thorne, 1986). When using RiverMod to analyze a river reach, the water and the sediment discharges constitute primary independent variables. The channel slope, width and depth constitute primary dependent variables (Vanoni, 1975). RiverMod attempts to quantify channel changes. The independent variables are input into the model and the corresponding dependent variables are estimated. In this way changes in both the dependent and independent variables of the river system can be simulated. 5 1.5.1 Adjustment of the Dependent Variables The dependent variables are mutually dependent and a river may adjust its width, depth and slope simultaneously (Henderson, 1966). There are three main governing equations for a river system associated with the independent variables. They are the flow resistance equation, the sediment transport equations and a bank stability relation (Henderson, 1966). Although there are an infinite number of combinations of width, depth and slope values for the given independent variables, Gilbert (1914) has shown through flume experiments that these may still be tending toward an optimum determined by the independent variables. The optimum is assumed here to represent an equilibrium condition. Figure 1.3 is a schematic representation of two solution curves illustrating the theory behind Gilbert's (1914) flume experiments. In Figure 1.3(a) discharge and slope are held constant and an optimum width develops which corresponds to the maximum sediment transporting capacity, similar to Gilbert's flume experiments, described above (Gilbert, 1914). In Figure 1.3(b) discharge and sediment transport are held constant and the optimum is now at the minimum slope. This second case mimics the behaviour of most natural rivers over engineering time scales (Millar, 1994). When a channel responds to a disturbance or a change in the independent variables, generally, the fastest changing dependent variables are the depth and the width. The slope is the slowest of these three variables to change (Booth, 1990). Figure 1.4 illustrates the dependent variables of the river initially fluctuating about the 6 equilibrium condition, until a large change is imposed. During adjustment to a new, postdisturbance equilibrium, the channel passes through a transient adjustment phase. The non-equilibrium transient phase could persist for years, or decades. Figure 1.4 also shows that at any point during the pre-disturbance stage, the dependent variables may not be in constant equilibrium, but may be in equilibrium in a statistical or average sense. 1.6 STREAM IMPACTS AND CHANNEL RESPONSES The river system, viewed as an ideal watershed, can be divided into three zones of erosion, transport and deposition (see Figure 1.5). The transport of sediment is continuous and any changes , or disturbances, which affect any reach within the catchment may affect the river downstream (Kondolf, 1997). Occasionally, cessation of the disturbance will return the stream back to its natural conditions. However, more often, "although riparian pressures have been removed, continued disturbance on a watershed scale may limit the ecological conditions and override the influence of channel adjustments" (Magilligan and MacDowell, 1997). The following sections discuss the independent variables, stream impacts and subsequent channel responses . (a) Peak flow discharge: Channels adjust to the dominant discharge, which is at or near bankfull discharge (about the 2-year flood), when most sediment transport occurs (Werritty, 1997). 7 (b) T h e sediment load and calibre: The sediment supply and t h e size o f t h e bed material will likely be a f f e c t e d by any change t h a t produces an instability in t h e flow regime (Hey and Thorne, 1986). (c) Bank stability: T h e bank type (cohesive or non-cohesive) determines whether t h e channel is more susceptible t o fluvial erosion or mass failure. Both fluvial erosion and mass failure may be functions o f parameters such as vegetation (Millar, 1994; M a c V i c a r , 1999). (d) T h e valley slope: T h e valley slope is considered relatively constant in engineering time scales; t h e r e f o r e , channel changes and impacts will only be discussed as they pertain to bank-full discharge, t h e sediment load and calibre, and t h e bank stability. 1.6.1 Urbanization Urbanization near a stream o f t e n signifies more impervious areas in t h e vicinity o f t h e stream. Typical examples o f urban development include t h e construction o f paved s t r e e t s and sidewalks, sewers, houses and other buildings, and parking lots. S u r f a c e cover might initially increase 1978) and later decrease t h e sediment load during construction (Dunne and Leopold, t h e availability o f sediment and increase t h e amount o f water reaching t h e s t r e a m , causing increased channel erosion or stream incision (Macklin and Lewin, 1997). 8 Although there may be an mcrease in the volume of runoff reaching the stream, many tributaries intended to carry runoff toward a channel, may get buried or eliminated. This might result in the disappearance of many smaller channels that under natural conditions assisted in keeping both sediment and runoff distributed. Rock Creek located in the Washington, D.C. Area, provides an example of such an event. Figure 1.6 shows the drainage net for this creek in 1913 and 1964, both before and after modern urbanization (Dunne, and Leopold, 1978). Figure 1.7 shows a graph of channel cross-section vs. drainage area for both rural and urban streams in Pennsylvania. I t shows that rural streams have, on average, smaller cross-sections than their urban counterparts of similar size (Dunne and Leopold, 1978). Another recent study comparing rural and urban catchments in southeastern Pennsylvania found that, overall, urban streams were wider and straighter than their rural counterparts (Pizzuto, et al., 2000). The bankfull sediment load may Increase or decrease depending on the type of development. During the development stages, effects such as land clearing may temporarily increase the sediment load, but once development is completed, if the area is well paved, the sediment load may decrease. The median grain-size (D50) is usually utilized as an indicator of the distribution. I t is a very difficult parameter to predict because the nature in which it is affected depends very highly on the type of development and the D50 of any new materials introduced to the system, or any materials removed from the system. Therefore in most cases it is assumed to remain constant (Personal 9 Communication, Millar, 1998; Human Impacts on River Systems (CIVL 598) lecture notes University of British Columbia). The effect of urbanization on bank stability is also difficult to predict. Postdevelopment stabilization of the banks would increase the bank stability. Without it, with the removal of riparian vegetation it would decrease. 1.6.2 Logging and Agripulture In forested drainage basins, trees are normally cut down in order to use the wood or for farming or grazing to take place (Macklin and Lewin, 1997). (a) Logging: The most obvious effect of logging activity is the increase in sediment load conveyed into the channel. Large amounts of disturbed soil from increased activity, coupled with the upheaval of tree-root systems, cause channel instability by increasing the occurrence of mass soil movements and landslides, which can then increase the sediment load (Dunne and Leopold, 1978). The sediment size would probably remain unchanged, since any incoming sediment would consist of a larger quantity of the same source. In most areas affected by logging, bank stability is decreased. Buffer strips of adequate size, left adjacent to the streambanks, might reduce the impact of deforestation. Logging usually reduces the amount of large woody debris (LWD) in a 10 channel, further reducing the bank stability, since LWD can partially armour banks, add to channel roughness and act as a sediment trap. (b) Agriculture: Land is often cleared of trees to facilitate cultivation or grazing; therefore, the initial stages of agricultural development may have similar effects on a channel as logging. If land clearance spreads to steeper slopes, gullying and mass wasting can increase the proportion of coarse material entering the channel (Clark and Wilcock, 2000), and therefore change the grain-size distribution of the stream. This can adversely affect those species of wildlife dependent on the stream. Magilligan and McDowell (1997) conducted a study which investigated the geomorphic adjustment following the removal of cattle grazing. They observed that in those catchments where cattle grazing had been removed for at least 10 years the bankfull and low flow widths had narrowed and the bed had been remobilized into more pool area. 1.6.3 Dams and Reservoirs The primary effect the construction of a dam would have on a channel is the interruption of the sediment transport of the river. The bedload and all or some of the sediment load is trapped upstream of the dam in the reservoir. The water released by the dam may then erode the bed and banks downstream of the dam, decreasing the bank stability and reducing the total sediment load of the system. Figure 1.8 shows the 11 predicted natural yield of a river and the actual yield after dam construction. The actual yield is markedly less than the predicted natural value (Kondolf, 1997). Increased erosion could possibly change the particle size of the bed. In a gravelbed river this could signify larger grains downstream of the dam. This can threaten salmonid spawning sites if the bed becomes so coarse that the fish can no longer move the gravel (Kondolf, 1997). 1.6.4 Mineral Extraction (a) Instream Gravel Mining: Instream gravel mining - the extraction of gravel directly from the stream bed - may cause channel incision upstream and downstream of the mining site, due to sediment shortage, coarsening of the bed, and lateral channel instability. (b) Metal Mining: Before mining legislation, waste from mine sites was discharged directly into the nearest channel, altering the entire fluvial system. Changes include: hindering vegetation growth along the banks, negatively affecting the bank stability; changes in channel pattern over short periods of time (Macklin and Lewin, 1997); and the presence of toxic materials in the fluvial system due to the extraction and processing of metal ores for up to thousands of years (Macklin, 1996). 12 equilibrium 10,000 time (years) assumed equilibrium 200 time (years) Figure 1.1. Schematic representation of the graded condition in two time scales, (a) The dotted line represents the graded condition the river is trying to adjust to through changes in its slope, (b) This is a close-up of the box shown in (a). The scale is much and over 200 years the river may be said to be in equilibrium (represented by the dotted line). 13 VALLEY BED AND SLOPE CHANNEL SINUOSITY VELOCITY + + HYDRAULIC + + T DISCHARGE SLOPE RADIUS BANK SEDIMENT SEDIMENT OUTPUT OUTPUT/INPUT + It WETTEn PFRIMETER SEDIMENT INPUT Figure 1.2. Statistical model of a river system illustrating how the dependent variables are influenced by the independent variables. Valley Slope, Bed and Bank Sediment, Discharge and Sediment Input are the independent variables (located in the four corner boxes). A "+" indicates the two variables are directly related. A "-" indicates the two variables are inversely related (Hey, 1976 as cited by Hey, 1982). 14 (a) (b) Width Figure 1.3. A schematic representation of two solution curves showing optimal geometry, (a) Discharge and slope are held constant and an optimum width develops which corresponds to the maximum sediment transporting capacity, (b) Discharge and sediment load are held constant and the optimum occurs at the minimum slope (Millar, 1994). 15 dependent variable pre-disturbance transient post-disturbance phase time Figure 1.4. Illustration of a dependent variable vs. time. The equilibrium pre-disturbance stage, transient phase and post-disturbance stage are visible. The fluctuations about the equilibrium condition can be seen in both the pre- and post-disturbance stages. It is difficult to predict the precise path of the transient phase. 16 Figure 1.5. Idealized watershed showing three zones of erosion, transport and deposition (Kondolf, 1997). 17 Figure 1.6. Drainage net of Rock Creek in 1913, before modern urbanization and in 1964 after modern urbanization (Dunne and Leopold, 1978). 18 Figure 1.7. Channel cross-sectional area at bankfull vs. drainage area of both rural and urban streams in Pennsylvania (Dunne and Leopold, 1978). 19 •J} to N A T U R A L YIELD e \ o 2 S> 6h c E A C T U A L Y I E L D (WITH D A M <D EFFECT) i CO 01 .3 o 1940 1950 I960 1970 Figure 1.8. Reduction in sediment supply from the catchment of the San Luis Rey River due to the contruction of the Henshaw Dam (Kondolf, 1997). 20 Chapter 2 MODEL FORMULATION In this chapter the analytical basis for RiverMod is discussed, and the algorithms used to determine sediment transport, flow resistance, bank stability and planform geometry are formulated. 2.1 THEORETICAL BASIS The model presented in this thesis is an extension of a previous model developed by Millar and Quick (1993, 1998) and Millar (2000). I t models the hydraulic geometry of an alluvial gravel-bed channel. The main contribution of this model, as compared with previous ones, is the fact that it assesses the bank stability of the channel. The model operates under the assumption that a river will tend toward an equilibrium geometry. The dependent variables adjust according to the independent variables in order to reach an optimum geometry which Satisfy the discharge, the bank stability and the bedload constraint. The channel is assumed to have a trapezoidal crosssection (see Figure 2.1). 21 Extremal Hypothesis 2.1.1 T h e r e are seven channel geometry equations mentioned in M i l l a r and Q u i c k (1993b) and eight principal dependent variables. The equations are f o r flow resistance, continuity, velocity, mean bank shear stress, mean bed shear s t r e s s , bank stability and sediment transport. The variables are bed width (W) or perimeter (P), channel depth (Y), channel slope (5), f r i c t i o n f a c t o r (f), mean velocity (u), mean bed shear stress (ibed), mean bank shear s t r e s s (thank), and bank angle (6) As t h e r e are more variables than t h e r e are equations to solve f o r t h e variables, an additional relation is necessary. An "extremal hypothesis" is introduced which s t a t e s t h a t "the channel geometry will adjust until t h e sediment transport capacity of t h e channel is equal to t h e value supplied f r o m upstream" (Millar and Quick, 1993b). The e x t r e m a l hypothesis allows an optimum solution to be obtained. The main function of t h e model is to solve t h e seven equations with t h e extremal hypothesis in order to obtain an optimum channel geometry, p e r f o r m e d by a computer routine. 2.1.2 Constraints T h e r e are t h r e e constraints that must be s a t i s f i e d in o r d e r to f i n d an optimum hydraulic geometry f o r a given channel. These are bedload, discharge and bank stability constraints. The bedload constraint ensures t h a t t h e channel is in equilibrium, and t h e 22 sediment load is transported through t h e reach without net deposition or erosion (Millar, 1994). T h e discharge constraint ensures t h e discharge capacity of t h e channel is equivalent to t h e bankfull discharge (Qbf), given by: UA = Q (2.1) bf where U is average velocity, A is cross-sectional area and Qbf is bankfull discharge. T h e bank stability constraint ensures that the banks are stable with respect to both mass failure and fluvial erosion. 2.2 SEDIMENT TRANSPORT Hans A l b e r t Einstein defined a bedload formula as "an equation linking t h e rate of bedload transportation with t h e properties of t h e grain and of t h e flow causing t h e movement" (Einstein, 1942). Bedload is fluvial sediment transported along t h e bed of a channel by rolling, sliding or saltation. I t is very difficult to accurately measure t h e bedload transport rate of gravel-bed rivers, and we are generally f o r c e d to rely on predictive relationships. A t present, no universal bedload transport equation exists (Reid eta/., Einstein (1950) sums up one main reason: 1997). "sediment movement and river behaviour are inherently complex natural phenomena involving a great many variables". So far, t h e development o f bedload transport equations has relied on empirical and experimental work, mostly carried out in flume studies employing uniform bed materials (Reid et al., 23 1997) even though most natural gravel-bed rivers tend to have non-uniform bed material size distributions. Another obstacle in the quest for a universal bedload transport formula is the fact that the existing formulae can only be tested with field data (Reid et al, 1997). RiverMod contains four established, well-known sediment transport equations. Einstein-Brown (1950), Meyer-Peter and Muller (1948), Ackers and White (1973) and Parker, Klingeman and McLean (1982). These particular equations were included in RiverMod because: 1. The parameters required as input variables are usually readily available, or easy to measure; 2. The parameters they require are similar to each other, and therefore it is easy for the user to use more than one equation with available data; 3. The equations are all well-established and well known and have each been included in several comprehensive reviews of bedload transport equations; 4. The user can compare the different answers obtained by using the different sediment transport relationships; 5. The conditions under which a particular equation was derived might be similar to the conditions of the river reach in question. In the fixed slope model within RiverMod, the bedload transport equation is used to determine the sediment transport capacity (in kg/s or Ibs./s) of the stream. In the variable slope model, the sediment transport capacity becomes an input variable and it is 24 used in conjunction with the selected transport equation to determine the slope of the stream. Gomez and Church (1989) expertly summarize the general assumptions, which are usually followed when applying bedload transport formulae; (i) that the flow, sediment properties and bedload transport rate are constant for the period in question; (ii) that the bedload transport rate is a unique function flow and sediment parameters; (iii) that the maximum amount of bedload is being transported. 2.2.1 Einstein (1942). Einstein-Brown (1950) Einstein originally developed an empirical bedload equation in 1942, which he replaced in 1950 with an analytical version (Gomez and Church, 1989). Also in 1950, Rouse, Boyer and Laursen modified Einstein's 1950 function and named it after Brown, author of Chapter XII of the book Engineering Hydraulics(Rouse (ed), 1950, as cited in ASCE Task Committee, 1971). Most bedload transport equations are based on the theory that transport is a function of the excess of a flow quantity above the threshold value for initiation of motion of the sediment particles in a channel bed. The Einstein formula is an important exception to this trend. Based on numerous experiments on bedload, Einstein concluded that "a distinct condition for the beginning of transportation does not seem to exist" 25 (Einstein, 1942). There is a continuous relationship between bedload transport intensity and flow intensity (Reid et al, 1997). The original equation was developed as an empirical relation based on the probability of particle movement (Gomez and Church, 1989). F= - + * 3 6 1 / - 2 1 q F(p,~ P )9t*\P.- 6 1 / 2 Pf s y = 3 (2.3) Pf f 9& A P,-l,_g_ P 5 f R (2 2) ( 2 4 ) H Where F is a parameter for settling velocity, gb* is the dimensionless transport rate, u/ is the ratio of the forces acting on a particle, v is kinematic viscosity, g is gravitational acceleration, pf is fluid density, p is sediment density, D is the sediment s grain size, q is sediment discharge, S is channel slope and Rh is hydraulic radius. s Note that the dimensionless bed shear stress (T*) and u/ are inversely proportional to each other. The dimensionless transport rate (gb*) is calculated using: 26 g* = 0 / / r* < 0 (2.5) / / 0 < r* < 0.093 g* = 2.15e~°' ^ (2.6) ifr*> #,* = 40(r*) (2.7) 0.093 39 3 The sediment transport capacity (in kg/s) is determined using: *L = 2.2.2 * - ^Ao (2-8) M e y e r - P e t e r and Muller (1948) The Meyer-Peter and Muller formula is an extension of the Meyer-Peter formula derived in 1934 (Gomez and Church, 1989). I t is an empirical law of bedload transport based on experimental data. This formula is exclusively for the movement of bedload and suspended load is not considered (Meyer-Peter and Muller, 1948). YS - 0.047 (y - r) s /Q . r -r s 27 r J (2.9) 1 26 (2.10) n 1 u y*A S X (2.11) n where ib is specific bedload transport rate, y is the unit weight of water, y is the unit s weight of the sediment, Q B / Q is a reducing multiplier which accounts for the amount of discharge acting on the bedload, KB / K S is a reducing multiplier which accounts for form resistance, V is the depth of the channel, D90 is the ninetieth percentile bed surface grain diameter, u is mean velocity, n is Manning's roughness coefficient and n' is Manning's grain roughness coeffcient. Meyer-Peter and Muller concluded that the initiation of motion is dependent on a certain magnitude of the shear stress. They also utilize two different grain diameters in their calculations, D90 and D . This is to account for the difference in composition a between the grain composition of the stationary bed and that of the mobile bed. The description of the procedure used by RiverMod to calculate the sediment transport capacity follows. First, the grain fAonnmq's n (n'), is calculated using: (2.12) 26 28 Next the dimensionless bed shear is determined, applying: f/0 r T \ (2.13) {rO (s -l)j a s Then, the value of t* is evaluated, in order to determine the dimensionless bedload (gb*) using: //*-*'< 0.047 g* //*-*'> 0.047 _> * = 8(r*'-0.047) ^ b (2.14) = 0 (2.15) 3 6 Equation 2.9 can be reduced to Equation 2.15 by assuming KB/KG = 1 for planar beds and Q B / Q = 1 for wide channels and by substituting values of y - 26000 N/m , y = 9800 N/m 3 s and g = 9.8 m/s (Gomez and Church, 1989). 2 And finally, the sediment transport capacity (in kg/s) can be determined using: *L 2.2.3 (2-16) * V^V^W = Ackers and White (1973) Ackers and White developed a new framework for the analysis of transport data (Ackers and White, 1973). I t is based on both dimensional analysis and physical arguments. They derived three key dimensionless parameters F number), D gr (grain diameter) and G gr gr (sediment mobility (sediment transport) (Ackers and White, 1973). Their formulae can be seen below. 29 3 u 4 = GgrYs gr / gr 0 v / A gr (2.17) -1 (2.18) - J (2.19) =D 9<- ,1-n u: u 32log(10% (2.20) (Gomez and Church, 1989). Where u* is shear velocity, and A, C, m and n are all dimensionless functions of the grain diameter. For coarse sediments (D > 60), where gr only bedload is being transported, the four parameters are all constant: n = 0; A = 0.170; m = 1.50; C = 0.025 For mixed transport (1 < D < 60), parameters are calculated using: gr (2.21) n = 1 - 0.56logO,gr 0.23/ A= V / m = 9.66/ £ _ 0.14 (2.22) + 1.34 (2.23) 9r J jQ(2.86logO^.-(log^) -3.53) 2 30 (2.24) During field analysis, Ackers and White concluded that the D35 bedload value was the "most appropriate when making predictions of total sediment load" (Ackers and White, 1973). The Ackers and White method for predicting bedload transport is useful for grain sizes larger than 0.04 mm. I t seems to be fairly accurate when predicting initiation of motion, although the authors caution against its use in unsteady flow conditions. The description of the procedure used by RiverMod to calculate the sediment transport capacity follows. The dimensionless grain diameter (D ) is first calculated, gr using Equation 2.19. Then the sediment mobility number (F ) A, C, m and n are calculated gr using Equations 2.20-2.24. Next the sediment flux (X) is calculated using: X = 5 A ^ r (2.25) v U , G gr is determined using Equation 2.18 and U*, the shear velocity and calculated using: (2.26) U* = V_^5 And finally, the sediment transport capacity (in kg/s) can be determined using: e , b =XQ 31 (2.27) 2.2.4 Parker, Klingeman and McLean (1982) The Parker, Klingeman and McLean formula takes the form: 4 (2.28) = 14001 Y s // </> < 0.95 a* = 0 Q 5 // 0.95 < ^ if K > I- 5 0 < 1.65 (2.29) co* = 0.0025exp[l4.2(<z) 50 co' 6 5 11.2 -1) - 9.28(^ 50 0.822^ f V ^50 -1) ] 2 (2.30) (2.31) J (2.32) "50 (2.33) 0.0876 (Gomez and Church, 1989). Where <b is a shear stress ratio, co* is the dimensionless 5 0 bedload transport rate, T 5o d is the bed shear stress for the median sub-surface grain size, i is shear stress, dso is the median sub-Surface grain diameter and S is the Specific s gravity of the sediment. This formula was the first to account for an armour layer and its effect on bedload transport rates, which is present in almost all gravel rivers (see Figure 2.2). 32 Bakke et al. (1999) define the layer termed "armour" or "pavement" as the layer extending to "the bottom of the largest surface particle". The "subarmour" or "subpavement" was identified as the layer just below the pavement of equal thickness. Based solely on field work, the empirical relations for the dimensionless bedload transport (co*) were derived from bedload measurements from Oak Creek (Parker et al., 1982), a small, steep gravel-bed stream in Oregon. Parker et al. (1982) theorized that the armour layer varies with shear stress and sediment load and that for conditions exceeding the critical stress of the pavement, the bedload grain size distribution of the pavement is approximately the same as that of the subpavement, so that: D50 ~ dso They developed the "equal mobility hypothesis", which states that surface coarsening develops to render all size fractions equally mobile (Parker et al, 1982) once the critical condition for movement was exceeded (Gomez and Church, 1989). I f an equilibrium condition develops, a balance has been achieved between the bed material supply and the transporting capacity of the reach (Bakke et al, 1999). The description of the procedure used by RiverMod to calculate the sediment transport capacity follows. First, the shield's stress for median diameter of the subpavement (T 5o) d is calculated using Equation 2.32. Next the shear stress ratio (q> ) 50 is determined, using Equation 2.33. Then, Equations 2.29 - 2.31 are used to compute the dimensionless bedload transport (co*). And finally, the sediment transport capacity (in kg/s) can be determined using: 33 * L = * ^Aed^g{ss-\)d^ M (2.34) 2.2.5 Discussion of Sediment Transport Equations Table 2.1 summarizes the sediment transport equations discussed in the preceding sections. Table 2.1. Bedload transport formulae used by RiverMod (adapted from Gomez and Church, 1989) Formula Particle Size Range (mm) employed in derivation Comments Meyer-Peter and Muller (1948) 0.40 - 28.65: uniform sediments, mixtures and light-weight materials "should be used for pure bedload transport. Solid material rolling or jumping along the bed of a river" (Meyer-Peter and Muller, 1948) Parker, 0.60 - 102.0: "restricted to small-to-medium-sized Klingeman and natural mixture (field paved gravel bed streams with steep or McLean data) moderate slopes, which are not (1982) dominated by sand load" (Parker et al., 1982) Ackers and 0.04 - 4.94: uniform White (1973) sediments and light- arguments were used to develop this weight materials. general function (Ackers and White, dimensional analysis and physical 1973) Einstein 0.785-28.65: uniform restricted to steady, uniform flow, (1942,1950) sediments and light plane bed conditions (Gomez and materials Church, 1989) In 1989 Gomez and Church conducted an extensive review of 12 well-known sediment transport formulae for gravel-bed rivers. The review tested each equation with four sets of river data and three sets of flume data with maximum consistency in all tests conducted. The survey illustrated the fact that very different results can be obtained, depending on which sediment transport equation is used. 34 Figure 2.3 shows a comparison between observed and calculated data using the data from the Elbow River, a gravel-bed channel. I t also demonstrates that any sediment transport equation used will perform better if the conditions of the reach in question are similar to the conditions under which the particular formula was derived. These predictive relationships should be used in conjunction with another form of analysis, including, but not limited to, laboratory experiments, the use of air photos and field work. 2.3 FLOW RESISTANCE Many equations have been proposed to relate average velocity in open channels to flow resistance (Hey, 1979). RiverMod uses four flow resistance relationships to establish the average velocity of the flow through and iterative process, ultimately determining the width and slope of the channel. The four relationships are Manning (1891), Composite Manning (Horton, 1933), Keulegan (1938), and Jarrett (1984). Because each channel has its own characteristic bedforms, profiles, flow resistance and sediment transport, it is important to use a resistance equation which has been designed for use with a similar channel to the study channel (Bathurst, 1978). In order to obtain the velocity and ultimately the channel dimensions, the bankfull discharge, channel slope, grain size distribution and roughness coefficient (except for Jarrett (1984) Equations) are needed (Bathurst, 1997). 35 2.3.1 Manning (1891) The Manning formula, which was presented in 1891, is the most widely used of all uniform-flow formulae for open-channel flow computations (Chow, 1959). I t is known in many countries under several combinations of the names Sauckler, Hagen, Manning and Strickler (Williams, 1970). I t can be written as: U = -R?S* n (2.35) where n is Manning's roughness coefficient. In order to use the Manning formula, one must determine an adequate value for Manning's roughness coefficient, n, because its value cannot be directly measured. This value is normally estimated either from experience, by using an empirical relationship or from documented cases (Bathurst, 1982). Table 2.2 supplies average values of Manning's n for numerous types of natural streams. Table 2.2. Values of Manning's roughness coefficient nior natural streams (adapted from Chow, 1959) Type of Channel and Description n I Top (a) Streams on a (i) clean, straight, full stage, no riffles or pools 0.020 width < plain (ii) same as above, but more stones and weeds 0.024 (iii) clean, winding, some pools and shoals 0.027 (iv) same as above, but some weeds and stones 0.030 (v) same as above, lower stages, more 0.032 30 m ineffective slopes and sections (vi) same as (iv) but more stones 0.034 (vii) sluggish reaches, weedy, deep pools 0.047 36 (viii) very weedy reaches, deep pools, or 0.067 f loodways with heavy stand of timber and underbrush (b) Mountain (i) bottom: gravel, cobbles, and a few boulders 0.027 streams, no (ii) bottom: cobbles with large boulders 0.034 (a) Pasture, no brush (i) short grass 0.020 (ii) high grass 0.024 (b) Cultivated areas (i) no crop 0.020 (ii) mature row crop 0.024 (iii) mature field crop 0.027 (i) scattered brush, heavy weeds 0.034 (ii) light brush and trees, in winter 0.034 (iii) light brush and trees, in summer 0.040 (iv) medium to dense brush, in winter 0.047 (v) medium to dense brush, in summer 0.067 (i) dense willows, summer, straight 0.101 (ii) cleared land with tree stumps, no sprouts 0.027 (iii) same as above, but with heavy growth of 0.040 vegetation on banks, usually steep I I Flood plains (c) Brush (d) Trees sprouts (iv) heavy stand of timber, a few down trees, 0.067 little undergrowth, flood stage below branches (v) same as above, but with flood stage 0.081 reaching branches I I I Top (a) Regular section with no boulders or brush width > (b) Irregular rough section 30 m In order to use documented cases as reliable sources of values for Manning's n, one must first consider the factors affecting it, such as channel surface roughness, vegetation, irregularities, alignment, silting and scouring, and obstructions (Chow, 1959). The user inputs a value for Manning's roughness coefficient (n) and the model calculates the velocity using Equation 2.35 37 2.3.2 Composite Manning (Horton, 1933) Originally developed by Horton (1933), this equation is useful when determining the dimensions or velocity of a channel whose bed and banks display differing degrees of roughness (Horton, 1933). An equivalent roughness is derived using: (2.36) where nbank and nbed are Manning's roughness coefficients for the bank and bed respectively, Pbank is the wetted perimeter of the bank and Pbed is the bed perimeter. Horton (1933) assumed the velocity was constant over the area and found that the values of both nbed and nbank vary with depth. Figure 2.4 shows the variation in equivalent n with depth for rectangular and trapezoidal channels. The user inputs a value of n for the banks and one for the bed and the model calculates a composite value for n using Equation 2.36 and velocity using Equation 2.35. 2.3.3 Keulegan (1938) Keulegan (1938) applied principles already established for use in circular pipes to turbulent flow in open channels and derived a formula for velocity distribution and hydraulic resistance. 38 The Keulegan flow resistance equation, together with the Darcy-Weisbach equation, is one of the most well-known and widely-used equations for open channels. The two equations are given below: f = 2.03log 12^ (2.37) (2.38) where f is the friction factor and k is the equivalent roughness height. Although k is s s actually a measure of the effect of any roughness element on the flow (Bathurst, 1982). Originally, the value of k was determined experimentally for each channel type. s Others have since set it equal to a certain grain-size diameter, commonly D50, D65 and D90 (Bray, 1982). Most studies of k find that it is dependent on the grain-size distribution of s the sediment. For non-uniform conditions, k is usually expressed as: s (2.39) where C is a constant, D is the characteristic grain diameter of which x percentage of x x sediment is finer (Millar, 1999). As shown in Table 2.3, below, values of k can vary s greatly. 39 Table 2.3. Summary o f "best-fit" k values derived f r o m reach-averaged hydraulic s geometry (Millar, 1999) Source k Bray (1980,1982) s 6.8D 50 3.5D 4 8 Charlton e t al.. (1978) 3.5D90 G r i f f i t h s (1981) 5.0D50 Hey (1979) 3.5D Leopold e t al. (1964) 3.9D 4 Limerinos (1970) 3.2D84 Millar (1999) 2.9D 84 8 84 I n spite of t h e great uncertainty associated with these values, according to Millar (1999) they "still represent t h e best approach to estimating velocity in ungauged gravelbed rivers." The user inputs t h e value of k , which is t h e equivalent roughness height f o r t h e s reach. T h e model, then calculates f using Equation 2.37 and uses this value to determine the velocity using Equation 2.38. 2.3.4 J a r r e t t (1984) The J a r r e t t equation was developed to predict Manning's roughness coefficient, n, f o r high-gradient streams, with a slope greater than 0.002 ( J a r r e t t , 1984). I t is as follows: n = 0.39/?, 5 016 40 038 (2.40) Jorrett (1984) found that for uniform flow, Equation 2.40 could be substituted directly into Manning's velocity equation (Equations 2.35). Jarrett (1984) found that most guidelines available for high-gradient streams do not take into account how the roughness elements change, with a change in gradient and therefore depth. As part of Jarrett's study, 21 natural streams were used, which represented a wide range of channel type, width, depth, slope, roughness and bed-material size. Figure 2.5 shows the relationship between Manning's roughness coefficient and hydraulic radius, from which Jarrett deduced that the roughness decreases with flow depth (Jarrett, 1984). Equation 2.40 should be used for natural streams with: 1. stable bed and banks; 2. slopes from 0.002 - 0.04; 3. hydraulic radii which do not include wetted perimeter of bed particles, and; 4. little or no suspended sediment (Jarrett, 1984). The model calculates a value of n using Equation 2.40 and then input into Equation 2.35. 2.4 BANK STABILITY There are two types of erosion affecting the banks of gravel-bed rivers: mass failure and fluvial erosion. The precise way in which the banks and consequently the river geometry are affected depends on the type of erosion as well as the sediment 41 composition of the banks. This thesis deals with either strictly non-cohesive or cohesive bank sediment. A soil is considered cohesive if its properties are influenced primarily by clay and Silt Size particles, where cohesive forces exist between particles. Sediment whose properties are influenced by sand and gravel size particles are referred to as cohesionless, or non-cohesive (Craig, 1992). 2.4.1 Non-Cohesive Bank Sediment Mass failure of non-cohesive bank sediment occurs when the bank angle (0) exceeds the friction angle (<|)), which is not subject to fluid shear stresses. The factor of safety is given by: F =^ * tan6? (2.41) M where Tb k is the mean bank shear stress, an y is the unit weight of water, S is the specific s gravity of the sediment, and Dsobank is median bank grain diameter. Loosening of individual clasts through weathering can reduce the friction angle (Lawler eta/., 1997). Oversteepening of the banks due to erosion of the bank toe can increase the bank angle. Both of these processes may cause mass failure in banks composed of non-cohesive sediment. Assessing the stability of undrained and drained 42 banks is similar with the addition that failure may result from an increase in pore water pressure under submerged conditions (Thorne, 1982). The flow of water and sediment in the main river channel produces shear on the bed and banks. I f this stress, stresses which is proportional to the velocity gradient, is increased sufficiently, fluvial entrainment of the sediment will occur (Thorne, 1982). In order for a particular grain to remain stable, its frictional forces must resist gravity and the forces exerted by the fluid. The limiting bank stability is: (2.42) where Tbankc and tbedc are the critical shear stresses acting on a grain located on the bank and bed respectively (Chow, 1959; Henderson, 1966; as cited by Millar, 1994). In the original relation, can only achieve a maximum of 40°, but if the effects of bank vegetation and consolidation are included, § can be replaced with an in situ value which can reach a maximum of 90° (Millar, 1994). The final relation is: (2.43) This equation is used in order to determine whether the banks of a stream composed of non-cohesive bank sediment are stable with respect to both mass failure and fluvial erosion which accounts for the effects of bank vegetation and consolidation of the bank sediment (Millar, 1994). 43 2.4.2 Cohesive Bank Sediment In order to assess the stability of a channel bank with respect to mass failure, the motivating and resisting forces acting on the bank must be balanced. The motivating force is due to the weight of the soil and the resisting force is composed of the cohesion (c) and internal friction angle (<(>) of the soil (Millar, 1994). ^Aass failure occurs when the critical bank height is reached (Lawler eta/., 1997). The factor of Safety with respect to mass failure (FSH) must be greater than or equal to one. It is calculated using: = ^ FS ^ crif " = * N H (2 44) C Hy 1 ] t where H is the vertical bank height; H ht is the maximum height at which the banks C are stable; y is the saturated unit weight of soil; and N is a dimensionless stability s t number, given by: N s = 3.83 + 0.052 (90 - 0) - 0.0001 (90 - 6f where 9 is the bank angle (Millar and Quick, 1998). 44 (2.45) There are several ways mass wasting can occur. The mode of bank failure depends on the bank angle, material properties such as grain size and soil cohesion and the location of the groundwater table or channel water level. Fluvial erosion of cohesive banks usually signifies parcels of soil, not individual grains are being eroded. This is because of the strong cohesive forces that exist between the grains (Lawler et al., 1997). Erodibility is a function of many factors, including mineralogy, particle size, temperature of the water and moisture content of the soil. For example, hard, dry banks are more resistant than wet banks, which are more easily eroded (Thorne, 1982). In order for the channel banks to be considered stable, the factor of safety with respect to fluvial erosion (FS ) must be greater than or equal to one. I t is calculated T using: f=5 x = (2.46) bank 1 where T b Q n k is the mean bank shear stress and Quick, 1998). T ank b and T it can be calculated using: 45 c r is the critical bank shear stress (Millar 1 bank ys SFbank (2.47) 4K r where fp SF nk ba = 1-766 ^ +1.5 \ P \ bank -1.4026 (2.48) r 2.4.3 Effect of Vegetation on Bank Stability T h e r e have been many studies conducted to determine the influence o f riparian vegetation on t h e bank stability o f gravel-bed rivers (Hey and Thorne, 1986; Millar and Quick, 1993; Beeson and Doyle, 1995; Millar, 2000). Beeson and Doyle's (1995) study o f four streams in southern British Columbia compared bank erosion in vegetated and non-vegetated channel bends and found that riparian vegetation decreases bank erosion. They also concluded that t h e denser t h e vegetation, t h e more e f f e c t i v e it was at reducing bank erosion and t h e r e f o r e increasing bank stability. T h e theory presented by Millar and Quick (1993b) f o r channels with non-cohesive banks, was t e s t e d using published data values as input, and the output was compared to the observed channel geometry. T h e modified friction angle ((j)') was used t o quantify the influence o f bank vegetation and 40° is a reasonable value f o r non-vegetated banks (Millar, 2000). Millar and Quick (1998) t e s t e d t h e theory f o r channels with cohesive 46 banks and concluded that bank vegetation has an effect on the critical shear stress (tent), which has a large influence on the stable channel width. RiverMod incorporates a measure of bank vegetation by utilizing the internal bank friction angle ((()') as an indicator of bank vegetation strength (Millar and Quick, 1993). It increases with increased density of bank vegetation (Millar, 2000). 2.5 PLANFORM G E O M E T R Y Rivers are often classified based on their planform geometry. Braided channels are those with relatively stable alluvial islands and therefore two or more separate channels (Leopold and Wolman, 1957). A meandering river has a single winding and sinuous channel. The planform of a river is dependant on the bed slope and discharge. For a given discharge, meanders will occur on smaller slopes than braids (Leopold and Wolman, 1957). River slopes can adjust through the process of meandering. As a river creates meanders and the sinuosity of those meanders increases, the channel slope decreases. Figure 2.6 shows three streams with increasing sinuosity from straight to meandering. The valley slope is built-up over time through aggradation. The channel slope adjusts to a steady-state condition. Figure 2.7 is a flowchart showing the computations in the planform model. RiverMod determines the planform geometry of the river by comparing the bedslope to the valley slope and the transitional slope. I f the bedslope is greater than the valley slope, the channel is unstable and is said to be aggrading. I f the bedslope is less than the 47 valley slope, the channel is stable and is either braided or meandering. The transition slope is determined using: (2.49) where W is the river width and F is the Froude number (Parker, 1976, as cited by Millar, r 2000). The transition slope is compared to the bed slope to determine if the river is braided or meandering: 5 > S* = braided (2.50a) S < S* = meandering (2.50b) I f the river is meandering, the sinuosity (£) and meander wavelength (A.) are determined using empirical relations developed by Leopold and Wolman (1960). These empirical relations ordinarily use a meander wavelength multiplier between 10 and 12, the model uses a default value of 11. The model also utilizes Langbein and Leopold's (1966) sine-generated curve to calculate values for the maximum deviation angle along the channel ( O ) and the radius of curvature (r ). Langbein and Leopold (1966, as cited by Thorne, 1997), found that the c sine-generated curve resembled an idealized river, by approximating the path of least resistance in flowing around a bend. The sine-generated curve is defined by: (2.51) 0(0 = $ sin 48 where 6 is the channel deviation angle, cD is the maximum value of 0,1 is the distance along the valley axis and L is the distance along the meander arc length (see Figure 2.8). An iterative technique is used to solve the sine-generated curve for O , such that the river sinuosity (£j ) determined using 0 and O is equal to the sinuosity (cj) determined c from the bedslope and valley slope. Ultimately the maximum deviation angle is calculated and presented as output (See Figure 2.7). Equation 2.51 can be used to graph the idealized planform geometry of the river. The idealized meander parameters are useful for river restoration, or attempting to return a channel to its natural form. The planform geometry analysis may assist the redesign of sections of channel that have been straightened or channelized (Millar and MacVicar, 1997). 49 w Actual Channel Cross Section Simplified Trapezoidal Cross Section Figure 2.1. Definition Sketch of a simplified trapezoidal channel for Millar and Quick (1993) model (MacVicar 1999). 50 flow coarse upper layer fine lower layer • . • ' - >•••.'• . .v - • , : ; ? / S « J - subarmour Figure 2.2. A Schematic representation of the armour layer in a gravel-bed channel. 51 OBSERVED UNIT BEDLOAD TRANSPOHT RATE (kfl/m §-') Figure 2.3. Comparison between observed and calculated data for Meyer-Peter and Muller, Ackers and White, Einstein, and Parker, Klingeman and McLean formulae and Elbow River data (Gomez and Church, 1989). 52 Figure 2.4. Relation between grain-size diameter (D) in feet and Manning's roughness coefficient (n) in imperial values for rectangular and trapezoidal channels in which the bottom and sides have different roughness characteristics (Horton, 1933). 53 Figure 2.5. Relation of Manning's Roughness Coefficient to Hydraulic Radius for four of the streams Jarrett used in his study (Jarrett, 1984). 54 Figure 2.6. Rivers of increasing sinuosity from a straight channel (A) to a meandering channel (C) (Thorne, 1997). 55 Bedslope S, from Fixed or Variable Slope M o d e l Valley Slope S , Input by User v Calculate Transitional Slope S* Aggrading and Unstable Meandering Planform Braided Planform Calculate Sinuosity \ = S /S v % < 1.1 Straight Channel 1.1 > 3.1 Meandering £>3.1 Calculate Wavelength A. and Radius of Curvature Figure 2.7. Flowchart for the Planform Model. 56 X 4* =0 1 1=0 1 -Q=L/4 -R=L/2 I I 1 0=0 1 Figure 2.8. Definition sketch of one meander arc. 57 — ^=3174 1 i ^ - *i e= i Chapter 3 TYPE O F ANALYSIS 3.1 INTRODUCTION There are two versions of the original (Millar and Quick, 1993b) model and therefore two versions of RiverMod: a fixed-channel-slope model and a variable-channelslope model. Each has been formulated for channels with either non-cohesive or cohesive banks. 3.2 FIXED-CHANNEL-SLOPE MODEL The fixed-channel-slope version of RiverMod is equivalent to an experiment where the slope is fixed, and the channel width, depth and sediment transport rate adjust to the discharge. I t is useful for examining the effect of bank stability on the channel geometry (Millar, 1994). Discharge and bank stability are the only two constraints that need to be satisfied in the fixed-channel-slope version of the model. Since the value of the slope is input by the user, it is treated as an independent variable. The bedload constraint is only used in the variable-slope version (Millar, 1994). Figure 3.1 is a flowchart for the Millar and Quick (1993b) fixed-slope model. The model starts by selecting a trial value of the bed perimeter and the bank angle (6). (Pbed), bank perimeter (Pbank) These are the primary dependent variables. The other dependent 58 variables can be calculated from these three. Only one of the three primary dependent variables is varied at a time. Once the discharge constraint is satisfied, the channel dimensions are calculated and then the bank stability constraint is assessed. Finally the optimal solution is found at the point of maximum sediment transport (Millar, 1994). Details regarding the dependent and independent variables are discussed in Chapter 4. The fixed-channel-slope model has been formulated for channels with non-cohesive bank sediment and channels with cohesive bank sediment. The discharge constraint is the same for both of these channel types; however, each channel type has its own bank stability constraint (see Section 2.4). Millar (1994) noted that in channels with non-cohesive bank sediment, the bank vegetation tended to increase <(>' and therefore bank stability. In channels with cohesive bank sediment, bank vegetation tended to increase 3.3 T c r i t. VARIABLE-CHANNEL-SLOPE MODEL The variable-channel-slope model more correctly approximates the behaviour of d natural gravel-bed river, since the slope is variable. In the variable-channel-slope version of the model, channel slope (S) is treated as a dependent variable. The bankfull sediment discharge capacity (&bf) is considered an independent variable. This version is similar to the fixed-slope version, with an additional bedload constraint, which requires the bedload transporting capacity of the channel to be equal to the sediment load (Gb) (Millar and Quick, 1993b). 59 Figure 3.2 is a flowchart showing the steps that make up the variable-channelslope model. I t is similar to Figure 3.1, except that the channel slope is varied for trial values of PbedA strong relationship between the dependent variables is noticed when using the variable-slope model. The final output of the model are the dependent variables: width, depth and slope. They would change in response to changes the independent variables caused by alterations in the stream environment. As the independent variables change, the dependent variables will also adjust in accordance to the new conditions. The dependent variables adjust their values at different rates. According to Booth (1990) the depth and width of the channel adjust at a much quicker rate than the slope. Even though they adjust at different rates, the width, depth and slope of a river are all inter-related and it is unlikely one will change without the others also modifying their values. For channels with non-cohesive channel banks, Millar (1994) notes that if the width of the channel changes, the depth changes in order to satisfy continuity and the slope adjusts to satisfy the bedload constraint. For channels with cohesive channel banks, the model can be used to determine whether a stream is bank-height or bank-shear constrained. In general, only one of these bank stability constraints is active at a time. Millar (1994) documented the change in the active bank stability constraint by varying Qbf, 60 keeping t9 f b constant and then varying Gbf, and keeping Qbf constant. In both cases, the remaining independent variables were kept constant for both investigations. The results are shown in Table 3.1 Table 3.1. Results of Millar (1994) active bank stability constraint investigation. Bank-full Discharge (Q f) Bankfull Sediment Transport Active Bank Stability in m /s Capacity (Gbf) in kg/s Constraint b 3 Qbf < 100 100 < Qbf < 250 bank-shear constrained Gf > 5 b both bank-height and bank- 3 < Gbf < 5 shear constrained Qbf > 250 bank-height constrained Gbf < 3 Whichever bank stability constraint is not active, does not seem to affect the channel geometry (Millar, 1994). In summary, the f ixed-charinel-slope and variable-channel-slope models show that the bank stability greatly affects the geometry of the channel. <j>' is useful as a single parameter for assessing the effects of vegetation on the bank stability in streams with non-cohesive banks. In streams with cohesive banks, purpose. 61 T c r i t can be used for the same Select Trial/ , 3 bed Adjust Select Trial 0 bod Select Trial P» Adjust bank Figure 3.1. Flowchart for Millar and Quick (1993b) Fixed-Slope Model (MacVicar, 1999). 62 Input Independent X Variables J Figure 3.2. Flowchart for Millar and Quick (1993b) Variable-Slope Model (MacVicar, 1999). 63 Chapter 4 USING RIVERMOD This chapter functions as a user manual, by guiding the user while running RiverMod. The following sections illustrate a step-by-step process for the successful operation of the model. 4.1 INSTALLATION RiverMod is installed by inserting the CD and selecting "setup.exe". The installation program will guide the user through the setup process. 4.2 FILE SETUP SCREEN RiverMod prompts a choice between two options: the user can either open an existing file or create a new record, (see Figure 4.1). 4.2.1 Open an Existing File The user can choose to open a file saved earlier, or to return to the opening screen by clicking "Cancel". Once the user chooses which file to open, the name of the file will appear on the screen. Clicking "Open" (see Figure 4.2) will open the file and load the Computation Screen. 64 4.2.2 Creating a New File The user will be guided through the RiverMod screens, which enable a new file to be created. This sequence is demonstrated in the following sections. 4.3 MODEL SCREEN Once the user has chosen to create a new file, the Model Screen is loaded (see Figure 4.3). The Model Screen allows the user to select the Model Type, Bank Type and Unit Type. There are two different model types: the fixed slope model and the variable slope model. There are two different bank types to choose from. The model has been formulated for gravel-bed rivers with either non-cohesive or cohesive bank sediment. The model operates using the SI (Systeme Internationale d'Unites) system of units. However, the user may select between two systems of units: the British Gravitational (BG) system, when working with variables in imperial units, and the S I system when dealing with metric units. Once the user has selected the model, bank and unit type, the Equation Screen is loaded. 65 4.4 EQUATION SCREEN The Equation Screen presents the user with four sediment transport models and four flow resistance Einstein-Brown equations (see Figure 4.4). The sediment transport models are (1950), Ackers and White (1973), Meyer-Peter and Muller ( 1 9 4 8 ) and Parker, Klingeman and McLean (1982). The flow resistance equations are Keulegan (1938), Manning (1891), Composite Manning's (Horton, 1 9 3 3 ) and Jarrett (1984). The combination of equations will affect the model's ability to accurately predict channel changes. 4.5 COMPUTATION SCREEN This screen provides an interface to the nucleus of the model. I t is divided into an input section and an output section (see Figures 4.5 and 4.6). The options selected on the Model and Equation Screens determine which variables are needed as input and which will appear as output. 4.5.1 Model Input (a) Discharge Bankfull discharge ( Q f ) is used in either m /s or ft /s. Often Q f is taken to be 3 b 3 b the discharge with a two-year return period (Q2). (b) Roughness Coefficient The roughness coefficient depends on which flow resistance equation has been selected. 66 1. Keulegan - equivalent roughness (k ) in metres or feet. s 2. Manning's n - Manning's roughness coefficient (n). I t is unitless. 3. Composite Manning's n - Manning's roughness coefficient for the bed of the channel (nbed) and for the channel banks (nbank)4. Jarrett - this relation calculates a value for Manning's roughness coefficient; therefore no specific roughness input is needed. (c) Sediment Transport Each sediment transport relation requires a particular grain diameter to represent the sediment grain-size distribution. The units should be in metres or feet. 1. The Einstein-Brown model requires the median grain diameter of the bedload sediment (Dso). 2. The Ackers and White model requires the thirty-fifth percentile grain diameter of the bedload sediment (D35). 3. The Meyer-Peter and Muller model requires both the ninetieth percentile grain diameter (D90) and the average grain diameter (D ) of the bedload sediment. a 4. The Parker, Klingeman and McLean model requires the median sub-surface grain diameter (dso). (d) Bank Stability The bank stability constraint has been formulated for both cohesive and noncohesive bank sediments (Millar, 1994). For channels with non-cohesive bank sediment, the median grain diameter of the bank material (Dsobank) and the internal friction angle 67 (<|)') are necessary For channels as input data. Dsobank is in metres or feet and <$>' is measured in degrees. with cohesive bank sediment, soil cohesion (c), critical shear stress (t t) and cri the unit weight of bank sediment (y ) are needed as input data. Soil cohesion and critical t shear stress are measured in Pascals or lb/in (psi), and y is measured in N/m or lb/ft . 2 3 3 t (e) Channel Slope I f the fixed-slope model is used, the reach-averaged bedslope (S) appears on the input screen as an independent variable. I t is unitless and is the ratio of the vertical elevation change of the channel to the horizontal distance along the bed. In the variableslope version of the model, the slope appears only as an output and is therefore a dependent variable. (f) Sediment Transport Capacity I f the variable-slope model is used, the sediment transport capacity (Gb) appears as input data, measured in kg/s or Ib/s. In the fixed-slope version, it appears as the output data. 4.5.2 Model Output The surface width (W), mean depth (Y) and bed slope (S) are considered the three primary dependent variables. Surface width and mean depth are calculated in either metres or feet, while bed slope is unitless. The three primary dependent variables appear on the main Computation Screen. Bank height (in metres or feet), bank angle, (in degrees), and sediment transport capacity (in kg/s or Ib/s) further describe the 68 geometry of the channel and are considered secondary dependent variables. They may be viewed on the Supplementary Output Screen by clicking "More Details" (see Figure 4.7) For channels with cohesive bank sediment, two additional output variables appear on the Supplementary Output Screen, are H / H a - i t and They are both unitless TbankAcnt- and further indicators of bank stability. Table 4.1 shows the minimum and maximum values for each of the variables in order for RiverMod to run properly for both S I and Imperial Units. Some of the limits may extend beyond the normal range for a certain variable in a small gravel-bed river. Table 4.1. Minimum and maximum values for input variables SI Units Input Variable 0 - 25000 m /s 0.002 - 1 m 0.01 - 0.15 0.01 - 0.15 0.01-0.15 0.0001 - 0.5 m 0.0001 - 0.5 m 0.0001 - 0.5 m 0.0001 - 0.5 m 0.0001 - 0.5 m 0.0001 - 0.5 m 20 - 90 0 - 150 Pa 0 - 100,000 Pa 15, 000 - 40,000 N/m 0.0001 - 0.05 0 - 5000 kg/s Bankfull Discharge, Qbf Equivalent Roughness, k s Manning's Roughness Coefficient, n Bank Manning's Roughness Coefficient, nbank Bed Manning's Roughness Coefficient, nbed Thirty-fifth Percentile Grain Diameter, D35 Median Gram Diameter, D50 Ninetieth Percentile Grain Diameter, D 9 0 Average Grain Diameter, D Median Sub-surface Grain Diameter, dso a Median Bank Grain Diameter, Dsobank Modified Friction Angle, <)>' Critical Shear Stress, Imperial Units 0 - 883,000 f t / s 0.0066 - 3.28 f t 0.01 - 0.15 0.01 - 0.15 0.01-0.15 0.00033 - 1.64 ft 0.00033 -1.64 ft 0.00033 -1.64 ft 0.00033 -1.64 ft 0.00033 -1.64 ft 0.00033 -1.64 ft 20 - 90 0 - 0.02175 psi 0 - 14.5 psi 3 0 0 tcriticai Soil Cohesion, c Unit Weight of Bank Sediment, y 3 t Reach-Averaged Bed Slope, S Sediment Transport Capacity, Gb 69 3 95.5 - 254.6 lb/ft 0.0001 - 0.05 0 -11000 Ib/s 3 4.5.3 Drop-down Menu The drop-down menu at the top of the Computation Screen provides options for data manipulation. Each of these options is detailed below. 4.5.3.1 File (a) Entering Data Data is entered into the Computation Screen and is grouped into records. This enables the user to save several series together as one file. Each record contains the same sediment transport and flow resistance equations as well as the same model type, bank type and unit type. The "New", "Next" and "Previous" buttons enable navigation from one record to another. (b) Saving Data As data is entered into the Computation Screen, it is automatically saved into a file called "RIVERFILE.DAT". Once the user chooses to save the file, RIVERFTLE.DAT is renamed with the user-selected name. Input and output data are saved together as a "DAT' file. The save function operates similarly to other windows-based programs. (c) Open Existing File There are two ways to open an existing file. The first is to select this option on the File Setup Screen when the model is loaded. The second is to choose this option from the File menu, on the Computation Screen. 70 RiverMod's capacity to check for compatibility between saved options and new user-determined options is a powerful and innovative feature that extends previous models. For example, perhaps the file was originally saved using the Einstein-Brown sediment transport relation, and the user wishes to open it using the Ackers and White sediment transport relation. The program would open the files and all of the data common to both options would appear on the screen. The data not held in common by the chosen equations - in this case, the grain diameter - would not appear on the screen. 4.5.3.2 Edit (a) Create a New Record Choosing this option from the Edit Menu enables the user to add a new blank record to the file. The data in the new record is related to the data in the previous record because they share the same model options, but is input to a completely separate model run. The data contained in the text boxes of the new record can be entirely different from the data in previous or subsequent records. (b) Delete and Clear Current Record Deleting a record allows the user to delete the current record from the series. The series will then contain one less record. Clearing a record will leave all the text boxes in the current record blank, but the number of records in the series will remain the same. 71 (c) Search for a Record The "Search for a Record" option under the Edit Menu, permits the user to search for a particular record. This can be done by searching for a particular value of any of the variables contained within the series of records. Using a series of user-friendly screen prompts (see Figure 4.8), the model will locate the first record containing that value. This feature is particularly useful for a user working with a series containing a large number of records. 4.5.3.3 Output (a) Run Choosing "Run" from the Output menu executes the optimization model. Before running, the model automatically scans for three potential errors that can occur during data entry: • empty text boxes; • non-numeric • values outside the permissible range. characters: I f an error is located, the program pauses, alerts the user, and allows the user to correct the mistake. The user can then choose the "Run" option again, and the program will proceed with the execution of the optimization model. 72 There are two slightly different paths the routine can follow, one for the fixedslope model and one for the variable-slope model. Each is explained in more detail in Chapter 3 . Figure 3.1 shows the fixed-slope version of the optimization model and Figure 3.2 shows the variable-slope version of the model. The two methods are very similar and follow this general procedure (Millar and Quick, 1998): (i) Input the independent variables. (ii) The minimum and maximum bounds for (iii) The midpoints of (iv) Flow resistance and velocity are calculated, using the equation chose by the user, Pbed, 9 and 9 and Pbank are set. are determined. Pb k Qn at the midpoint. The bounds of Pbed, Pbank are adjusted until the discharge constraint is satisfied (UA = Qbf) and the flow resistance and velocity values are also adjusted for each new midpoint value of P (v) ba nk- Bank stability: a. For non-cohesive banks, Equation 2.43 must be satisfied for the banks to be stable, b. For cohesive banks, the bank angle is adjusted until F S H = 1 (Equation 2.44) and FS = 1 (Equation 2.46). I f the banks are not stable, the channel is too narrow, the T bounds for Pbed are reset and the procedure returns to step (iv). I f the banks are stable, the routine proceeds to step (vi). 73 (vi) The bedload transport is calculated using the equation chosen by the user. I f the Pbed value corresponds to the maximum bed load transport, then the optimum solution has been reached and the dependent variables width (W), depth (Y) and bank angle (0) are calculated as output. I f the Pbed value does not represent the optimum solution, the procedure returns to step (iii). The final solution represents the steady-state equilibrium geometry of the reach (Millar, 1994). (b) Planform Geometry This option takes the user to the Planform Screen (see Figure 4.9 and Section 2.5). RiverMod will determine the planform geometry for the record which is open when this option is chosen. This feature of the model is useful for river restoration or design purposes (see Section 4.6). 4.5.3.4 Microsoft Excel (a) Open Microsoft Excel Choosing this option from the Output menu opens the spreadsheet program Microsoft Excel (provided it is available on the user's computer). In order for it to function properly, Microsoft Excel must be closed before this feature is executed. Integration of RiverMod with Microsoft Excel adds power and flexibility, allowing the user to: • create graphs and charts; 74 • visually compare the relationship of one variable to another; • visually compare the results obtained using different equations or model options. The model will automatically identify the options chosen by the user and transfer them to the Microsoft Excel worksheet. The worksheet will have a column for each variable with headings, which include the name and units of each variable. Figure 4.10 shows an example where the discharge was increased by 25 m /s for 5 runs, while all 3 other variables were held constant. Once the user has opened Microsoft Excel, all of its functions are independent of those of RiverMod. I f the user wishes to create graphs or move and delete data, he/she must do so from within the Microsoft Excel program. In order to save any new worksheets, graphs or files created in Microsoft Excel, the user must also do this from within Microsoft Excel and not by using RiverMod. As soon as RiverMod is closed, Microsoft Excel, opened by RiverMod will also be closed, so any file handling within Microsoft Excel should be done before RiverMod is terminated. (b) Append to Microsoft Excel Spreadsheet There are two separate append features built into the model. In order to use these features, Microsoft Excel must already have been opened through RiverMod. Both "append" commands add to the spreadsheet the input and output data, found on both the Computation Screen and the Supplementary Output Screen. 75 • Append Current Record to Spreadsheet - this feature is useful for importing the data in the current record into the Microsoft Excel Spreadsheet. • Append All Records to Spreadsheet - this feature allows transfer of the data in all the records into the spreadsheet simultaneously. The "append all" feature can only be executed from Record 1, so that the model can run through the records in order. 4.6 PLANFORM S C R E E N The user can access the Planform Screen (see Figure 4.9) from the Computation Screen. RiverMod will use the valley slope, entered by the user, in conjunction with a formulation for the transitional slope (S*) to determine the geometry of the river. I f the channel slope is greater than the valley slope, the river is aggrading and unstable. Otherwise, a meandering-braiding transition formulation is used to determine the river planform. Equation 2.49, developed by Parker (1976, as cited in Millar 2000), is used to determine the transitional slope. The transition slope is compared to the bed slope to determine if the river is braided or meandering, using Equations 2.50a and 2.50b. I f the river is meandering, the sinuosity (£) and meander wavelength (X) are determined using empirical relations developed by Leopold and Wolman (1960). These empirical relations ordinarily use a meander wavelength multiplier between 10 and 12. The model uses 11 as the default value, but the user may enter an alternate value. 76 The model also utilizes Langbein and Leopold's (1966) sine-generated curve (see Section 2.5) to calculate values for the maximum deviation angle along the channel (<D) and the radius of curvature (r ). c This information is useful for river restoration, which attempts to return the channel to its "natural" form. This may require the redesign of sections of channel that have been straightened or channelized. 77 78 79 model Type and Units UlOlxl Model Type OK (* Fixed Slope Model Cancel C Variable Slope Model C Calibration Model Bank Type- Units <* Non-cohesive Bank Sediment C SI Units Cohesive Bank Sediment C Figure 4.3. Model Screen 80 Imperial Units Row Resistance and Sediment Transport Options Choose a Sediment Transport Model Einstein-Brown (1950 i Ackers and White (1973) Meyer-Peter and Muller (1948) Parker, Klingeman and McLean (1982) Choose a Flow Resistance Equation Keulegan [1938 Manning's n (1891) Composite Manning's n (Horton, 1933) Jarrett (1984) Back OK Figure 4.4 Equation Screen 81 Record 1/ 3 - C:\s-rreaml.dat Output Microsoft Excel N on-Cohesive M odel I nput: I ndependent Variables SI Units Discharge (Keulegan Equation] Bankfull Discharge, Q. Equivalent Roughness, k s Sediment Transport (Einstein-Brov-in Median Grain Diameter, D 5& Bank Stability Median Bank Grain Diameter, D. M odified Friction Angle, ^ 1 Channel Slope Reach-Averaged Bed Slope, S Jon-Cohesive Model Output: Reach-Averaged Geometry Surface Width, W More Details > Mean Depth, Y Bed Slope, S Change Options Figure 4.5. Computation Screen for Non-Cohesive Sediment with Einstein-Brown and Keulegan equations in SI Units, using the Fixed Slope Model. 82 Record Edit 1/1 Output Microsoft Excel ohesive Model Input: Independent Variables SI Units Discharge (Composite Manning's Equation Bankfull Discharge, & M Bank Manning's n, Bed Manning's n, Sediment Transport (Ackers and White] Thirty-fifth Percentile Grain Diameter, Pss Bank Stability Critical Shear Stress, Soil Cohesion, c Unit Weight of Bank Sediment, Yt Channel Slope 0.001 Reach-Averaged Bed Slope, S Cohesive Model Output: Reach-Averaged Geometry Surface Width, W Mean Depth, V Bed Slope, S Change Options Figure 4.6. Computation Screen for Cohesive Sediment with Ackers and White and Composite Manning equations in SI Units using the Fixed Slope Model. 83 (a) Figure 4.7. Supplementary Output Screen for: (a) Non-Cohesive Sediment Model, with Einstein-Brown and Keulegan equations in SI Units, using the Fixed Slope Model; and (b) Cohesive Sediment Model with Ackers and White and Composite Manning equations in SI units, using the Fixed Slope Model. 84 J-^TjSelect a Variable: p. jBankfull j jDisdh^gd E ^ Grain Diameter r- Critical ~ Shear Stress ( Bank Manning's n Bed Manning's n r r Soil Cohesion C Unit Weight of Bank Sediment 0 , S l o p e Figure 4.8. Search Screen showing Cohesive Sediment Model Options. 85 Planform Model '"Determine Planform Characteristics Bed Slope I Valley Slope Corinpute'Rivef] Geometru 0 0 0 1 Meander Wavelength Multiplier The river is meandering in planform Sinuosity, t; r Maximum deviation angle along the channel $ mm F Wavelength, X 620 m Radius of Curvature, rc 127.6 m Figure 4.9. Planform Screen. This example shows a meandering stream. 86 Figure 4.10. Microsoft Excel Spreadsheet with values appended from RIVERMOD for ai example using Non-Cohesive Sediment with Einstein-Brown and Keulegan equations in SI Units, using the Fixed Slope Model. In this example, the discharge was increased by 25 m /s for 5 runs, while all other variables remained constant. 87 Chapter 5 SAMPLE APPLICATIONS 5.1 INTRODUCTION The following sections present two examples for using RiverMod. The model is applied to two rivers in British Columbia: Slesse Creek and Narrowlake Creek. Both creeks have been affected by logging, resulting in changes to the bank stability, river width and planform geometry. 5.2 S L E S S E CREEK Slesse Creek is a mountain stream situated in southwestern British Columbia and northwestern Washington State (see Figure 5.1). The upper catchment of the river, located on the U.S. side, is in a protected, pristine forested area. The lower reaches of the river, on the Canadian side, have been subjected to clear-cut logging along the banks and floodplain. As a result of this disturbance on the banks, the stream's geometry has changed with active widths increasing by over 3 times from pre-logging to post-logging conditions (MacVicar, 1999). The morphology has changed from a single channel meandering planform to a wide unstable braided planform (see Figure 5.2). This river provides an excellent case Study for RiverMod. The nearby forestry activity has not affected water and sediment Supply, but riparian logging has caused a decrease in the bank stability, due to the loss of bank vegetation. Therefore, the only 88 independent variable to suffer significant change is the internal friction angle, <)>', which decreased from 73° to 40° from 1936-1993 (MacVicar, 1999). 5.2.1 Investigation In order to apply RiverMod to Slesse Creek, the input parameters must be defined for past and present values. f^acV\car (1999) concentrated his analysis on Reach D of Slesse Creek (see Figure 5.3) because the reach is alluvial, the morphology is relatively homogeneous, fish habitat has been greatly decreased and recent restoration efforts have been focused on this reach (Babikaiff and Associates, 1997 as cited by MacVicar, 1999). Site surveys and historic aerial photograph analysis have shown that from 19361993 the river width has increased from 28 m to 145 m (see Figure 5.2) while the slope has remained constant at around 0.02 (MacVicar, 1999). The constant slope indicates the fixed-channel-slope version of the model should be used. Table 5.1 shows the reach-averaged channel geometry acquired from air photos (MacVicar, 1999). Table 5.1. Slesse Creek channel geometry (MacVicar, 1999) Year 1936 1973 1993 Width (m) 28 Sinuosity 1.15 21 1.12 1.06 Slope 0.019 0.020 0.021 Planform Single thread, wandering Single thread, wandering Braided 89 145 The following list describes the method for determining input data, as discussed by MacVicar (1999). The input parameters used are shown in Table 5.2. • Bankfull discharge (Qbf) - Gauge records for Slesse Creek were used to determine the bankfull discharge. For the years 1993 and 1973, the short-term mean of flows was used and for 1936, the long-term mean was used. • Flow resistance (k , n) - No roughness information was available for Slesse Creek. s Since the Jarrett (1984) Equation calculates a value of Manning's n (see Equation 2.40) based on hydraulic radius and slope, this option was selected while using RiverMod. • Sediment size ( D ) - A pebble count was taken, where x D50 = D50 bonk and the sediment size was assumed to remain constant over the years because no historical particle size distribution information was available. • Modified internal friction angle (<)>')- For vegetated banks, <t)i936and ^1973 were determined by calibrating RiverMod. This was done by using the known channel geometry for each year and entering varying values of ty' until the modelled channel geometry matched the observed channel geometry. For unvegetated banks, ^1993 = 40° was used (Millar and Quick, 1993b). • Channel slope (S) - S 1 9 9 3 was determined through field surveys. S1936 and S1973 were calculated using sinuosities determined from air photos and by comparison with S1993. 90 Table 5.2. Input variables for Slesse Creek (adapted from MacVicar, 1999) Variable 1936 1973 1993 Qbf (m /s) 92 67 117 0.11 0.11 0.11 (m) D (m) 0.133 0.133 0.133 0.166 0.166 0.166 D o (m) 0.33 0.33 0.33 DsObank (m) <t>' (°) calibrated 0.133 0.133 0.133 70 73 - - - 40 0.019 0.020 0.021 3 D (m) 3 5 D50 Q 9 using RiverMod f (°) estimated S 5.2.2 Analysis The fixed-channel-slope model of RiverMod for rivers with non-cohesive bank sediment was used to analyze Slesse Creek. The selected sediment transport relation was the Einstein-Brown Equation (1950) and the selected flow resistance formula was the Jarrett Equation (1984). The input values listed above were used to run several tests. All input values were held constant and was varied between 40 and 90° in order to establish any influence on river width. This was done for 1936, 1973 and 1993 values. The results can be seen in Figure 5.4. For each of the three years, as <)>', and therefore bank stability, decreases, the width of the river increases. For 1993, the value was extrapolated to less than 40° because =40° resulted in a width of 119 m. The observed width for 1993 was 145 m, which according to RiverMod, 91 corresponds to a value of 35°. Figure 5.4 shows the observed geometry for each of the years. For the rest of the analysis, however, <)>' = 40° was used, the generally accepted minimum value that applies to non-vegetated banks In order to make restoration recommendations, <(>' was plotted against width and depth using 1993 values (see Figure 5.5). The Millar (2000) meandering-braiding transition slope equation, was used to determine the value of necessary to change the river morphology back to a meandering pattern from its present-day braided pattern. I t is given by: 5* = 0 . 0 0 0 2 £ ° y 6 1 7 5 0 Q° Z5 b (5.1) where S* is the transitional slope, D50 is the median grain diameter, ()>' is the modified internal friction angle and Qbf is the bankfull discharge. Millar (2000) improved on the Parker (1976, as cited by Millar 2000) criterion (see Equation 2.49) by including the effects of bank Stability. Dsobank was also plotted against width and depth for the 1993 value (see Figure 5.6). Using the Millar (2000) meandering-braiding transition slope equation, the Dsobank necessary 5.2.3 to change the stream back to meandering was determined. Results The RiverMod results show that the bank stability parameter (((>') has effectively decreased from 70° in 1936 to 40° in 1993. This decrease, caused by logging on the 92 f loodplain and banks, has led to channel widening. Figure 5.4 shows the variation of <|>' with W for the bankfull discharge of each test year. River restoration could consist of increasing <)>' and/or Dsobank to achieve a narrower, deeper, more stable, single-thread, channel. Holding all other variables constant, including Dsobank at 0.133, it was established that increasing <(>' to 57° would decrease the river width to 62.2 m. A similar analysis, which held <)>' constant at 40°, determined that increasing Dsobank to 0.36 m would decrease the width to 47.8. These two width values, calculated using Equation 5.1, were established as the transitional width between a meandering and braided stream according to the above conditions. Either of these two strategies should change the river back to meandering, over time. Riparian planting and re-vegetation of the f loodplain is a way of increasing <$>' and therefore stabilizing the banks (MacVicar, 1999). This analysis shows the influence the bank stability parameter, <)>' and Dsobank have on W, assuming all other variables are held constant. Several errors are associated with this, since several variables including the sediment size distribution and roughness data were not available for past years. These values were assumed to be the same as the 1993 measured values. The model predicts a change in planform from meandering in 1936 to braided in 1993. Since this has been determined from air photo analysis (MacVicar, 1999), RiverMod is a good predictor of planform geometry in the case of Slesse Creek and can successfully be used for restoration recommendations. 93 5.3 N A R R O W L A K E CREEK Narrowlake Creek is located in the central interior of British Columbia about 80 km southeast of Prince George (see Figure 5.7). I t is a tributary to the Willow River and its f loodplain was extensively logged in the 1960's and 1970's (Wilson, 2001). 35% of the watershed was harvested, including 80% of the riparian bank vegetation (Berry, 1996 as cited by Wilson eta/., 2001). Air photos from 1946 and 1997 (see Figure 5.8) clearly show that the morphology of Narrowlake Creek has changed from a meandering channel before logging to an almost braided channel after logging. The change in Narrowlake Creek's planform geometry over time has been forced by logging on the banks. Accompanying the change to a "semi-braided" morphology is a marked change in channel width. Channel width has doubled from 1946 to 1997, increasing from 29 m to 58 m (Wilson, 2001). Narrowlake Creek is similar to Slesse Creek in that both creeks have undergone a change in planform geometry following logging and disruption of the banks. However, according to Wilson (2001), Narrowlake Creek has experienced a much higher level of disturbance, including the cumulative impacts of road development, stream crossings, vegetation removal and channel avulsions. As a result, sediment supply and upstream conditions of Narrowlake Creek have changed. This is unlike Slesse Creek, where the upper watershed conditions are still pristine. 94 5.3.1 Investigation In order to apply RiverMod to Narrowlake Creek, the input parameters must be defined for past and present values. The following analysis was conducted on treatment Reach 3 as described by Wilson (2001). Historic aerial photographs have shown that from 1946-1997 the river width has increased from 29 m to 58 m while the slope has remained constant at around 0.007 (Wilson, 2001). The constant slope indicates the fixed-channel-slope version of the model should be used. Table 5.3 shows the reach-averaged channel geometry acquired from air photos. Table 5.3. Narrowlake Creek channel geometry. Year 1946 1997 Width (m) 29 58 Slope 0.007 0.007 Planform Single thread, wandering Braided The following list describes the method for determining input data, derived from Wilson (2001). The input parameters used are shown in Table 5.4. • Bankfull discharge (Qbf) - Gauge records do not exist for Narrowlake Creek. Using gauge information from the Willow River, which Narrowlake Creek flows into, an interpolation was performed in order to determine the bankfull discharge (Personal Communication, Andrew Wilson, 2002). 95 • Flow resistance (k , n) - No roughness information was available for Narrowlake s Creek. Since the Jarrett (1984) Equation calculates a value of Manning's n (see Equation 2.40) based on hydraulic radius and slope, this option was selected while using RiverMod. • Sediment size ( D ) - Sediment size information was only available for 1997. With x no historical sediment size data, it was assumed that D50 • and Dsobank, Values for both D1946 = D1997. were obtained. Modified internal friction angle (<)>') - For vegetated banks, c > j9 i4 6 was determined by calibrating RiverMod. This was done by using the model inputs along with the known channel geometry for 1946 and entering varying values of ((>' until the modelled channel geometry matched the observed channel geometry. For unvegetated banks, 4)1997 = 40° was used (Millar and Quick, 1993b). • Channel slope (S) - S 1 9 9 7 was measured through field surveys. determined using air photos and by comparison with S1946 was 51997. Table 5.4. Input variables for Narrowlake Creek (adapted from Wilson, 2001) Variable 1946 1997 Q b f (m /s) 31.1 31.1 (m) 0.055 0.055 (m) 0.058 0.058 (m) 0.093 0.093 (m) 0.088 0.088 0.047 0.047 3 D 3 5 D o 5 D a D90 Dsobank (m) <|>' (°) calibrated 55 using RiverMod 96 (°) estimated f S 40 0.007 0.007 5.3.2 Analysis The fixed-channel-slope model for rivers with non-cohesive bank sediment of RiverMod was used to analyze Narrowlake Creek. The selected sediment transport equation was the Einstein-Brown (1950) and the selected flow resistance equation was Jarrett (1984). The values listed in Table 5.4 were input to RiverMod. ()>' was varied between 4 0 ° and 90° to determine its influence on river width, for 1946 and 1997 data. The results can be seen in Figure 5.8. As , and therefore bank stability decreases, the width of the river increases. For 1997, the <)>' value was extrapolated beyond 4 0 ° because 4 0 ° resulted in a width of 50.8 m. This is a reasonable approximation of the observed W1997 of 58 m. Figure 5.9 shows markers for the observed geometry for each of the curves. In order to make restoration recommendations, <)>' was plotted against width and depth for the 1997 values (see Figure 5.10). The Millar (2000) meandering-braiding transition formulation was used to determine the value of ^' necessary to change the river morphology back to a meandering pattern from its present-day braided pattern. Dsobank w a s plotted against width and depth for the 1993 value (see Figure 5.11). Using the Millar (2000) meandering-braiding transition formulation, the Dsobank necessary to restore the creek's meandering planform was determined. 97 5.3.3 Results RiverMod predicts a meandering channel f o r 1997. Although the present channel is no longer meandering, it is not quite braided either. Since 1946, the transition slope has been approaching the limit between meandering and braided rivers. T h e transition slope may approximate conditions necessary f o r t h e gradual change from one planform to another. In Wilson (2001) it is assumed that <)>' = 70° f o r 1946 and f = 40° f o r 1997, but RiverMod calibration finds <\>' = 55° f o r 1946. These refined values give the same results as those obtained by Wilson (2001) which show the creek to be meandering in both years, when actually it is close to braided in 1997. Once calibrated, the model still shows a t r e n d towards a braided river. According to what can be seen from air photos and field surveys, t h e river has begun to change from meandering to braided over t h e 30 year period between 1946 and 1997, primarily due to logging activity. Although RiverMod does not show this same result, the t r e n d towards a braided river is evident. Since the modified friction angle (<)>') has decreased due to logging and has led to channel widening, any river restoration program should consist of increasing <|>' and/or Dsobank in order to stabilize the banks. These two methods would aid in t h e development of a narrower, deeper, more stable, single-thread, channel. ()>' could be increased by 98 riparian planting and re-vegetation of the f loodplain. Dsobank could be increased by the addition of larger, stabilizing gravel along the banks of the river. Since RiverMod did not accurately predict the change in planform for Narrowlake Creek, credible values of <t/ and Dsobank were not obtained. According to results obtained by using RiverMod, the bank stability values at the meandering-braiding transition are <)>' = 34° and Dsobank = 0.035 m. These values are lower than the 1997 values, and would most likely decrease the bank stability further, if used in restoration planning. Errors contributing to this could be the lack of historical sediment Size data and roughness values. Since no gauge exists on Narrowlake Creek, the bankfull discharge was extrapolated, and therefore may not be the true bankfull discharge. The level of disturbance on Narrowlake Creek is much higher than that of Slesse Creek (Wilson, 2001). This could account for the higher degree of success in the application of RiverMod to Slesse Creek over Narrowlake Creek. The higher level of disturbance on Narrowlake Creek could also signify a larger level of uncertainty associated with input data. I t could also mean that there are more factors influencing the decrease in bank stability than RiverMod can account for. The model should predict a change in planform from meandering in 1946 to braided in 1997, since this has been determined from air photo analysis (Wilson, 2001), RiverMod is a good predictor of planform geometry in the case of Slesse Creek and has provided useful information in the case of Narrowlake Creek. 99 Figure 5.1. (a) Location of Slesse Creek Watershed; (b) Slesse Creek Watershed (MacVicar, 1999) 100 (a) Figure 5.2. Airphoto* of Slesse Creek: (a) 1936,1:22000 scale; (b) 1973,1:19050 scale; (c) 1993,1:17650 scale (MacVicar, 1999). 101 Chilliwack River Scale 1:22,200 Figure 5.3. Slesse Creek lower watershed, showing study reach D in 1936 (MacVicar, 1999). 102 a >_n o a o ON c_ o 4 >£. +W £ o cn o oo -t- o £_ "a <j "a -tO o a s o V) _v> o £ in O Ol •o cn c a C O o H— "O W o I D O £ > vO CO ON *—i II O 0) £_ a cr c o '•»a t_ ro .a ON "a <J CD £_ <o o cn c a CD V) J— V) ro . (ui) M+p.iM in" ^n—t cn 103 ON ON o L. CO I CM 1 00 (ui) qjdaa CM ^1 O T-ii 1 1 104 1 O u P (ui) u,4.daq o vO o -w—t o o Cd O y—t f~i o o 00 O (UJ) H+pjM 105 Eo „ f * o 00 o w l_ ^ LL Figure 5.7. Location of Narrowlake Creek W a t e r s h e d 106 Figure 5.8. Airphotos of Narrowlake Creek: (a) 1946; (b) 1997 (Wilson, 2001). 107 II o o >- 4- £ o cn T3 o +a 00 £_ £ "a o "a 3 -t<J a 5 o o o W o cn •o £ >lO _o cn c a o c M— "O o in O £ 2 's c o o a i_ •4- .D "a o i_ VJ o CO ("0 H+P.'M 108 o s s_ _ a Z ON T—I II I_ a 3 in CT 3 ^ ON (UJ) u,j,daq REFERENCES Ackers, P., and White, W.R., 1973: "Sediment transport: New approach and analysis." Journal of the Hydraulics Division, ASCE, 99 (HY11), 2041-2060. American Society of Civil Engineers (ASCE) Task Force of the Committee of Hydromechanics, 1963: Friction factors in open channels, Journal of the Hydraulics Division, ASCE, 89 (HY2), 97-143. American Society of Civil Engineers (ASCE) Task Committee, Sediment transportation mechanics, 1971: Sediment discharge formulas, Journal of the Hydraulics Divisiom, ASCE, 97 (HY4), 523-567. Bakke, P.D., Basdekas, P.O., Dawdy, D.R., Klingeman, P.C., 1999: "Calibrated ParkerKlingeman model for gravel transport." Journal of Hydraulic Engineering, 125 (6), 657660. Bathurst, J.C., 1978: "Flow resistance of large-scale roughness." Journal of Hydraulics Division, ASCE, 104 (HY12), 1587-1603. the Bathurst, J.C., 1982: "Theoretical aspects of flow resistance." In Gravel Bed Rivers, R.D. Hey, J.C. Bathurst, and CR. Thorne (Eds.), John Wiley & Sons, Chichester, 83-108. Bathurst, J.C, 1997: "Environmental river flow hydraulics." In Applied Fluvial Geomorphology for River Engineering and Management, Thorne, C.R., Hey, R.D., and Newson, M.D., (Eds.). John Wiley and Sons, 69-93. Beeson, C.E., and Doyle, P.F., 1995: "Comparison of bank erosion at vegetated and nonvegetated channel bends." Water Resources Bulletin, 31 (6), 983-990. Blench, T., 1957: Regime behavious or rivers. Butterworths, London, England, 138p. Booth, D.B., 1990: "Stream-channel incision following drainage-basin urbanization." Water Resources Bulletin, 26 (3), 407-417. Bray, D.I., 1979: "Estimating average velocity in gravel-bed rivers." Journal of the Hydraulics Division, ASCE, 105 (HY9), 1103-1122. Bray, D.I., 1982: "Flow resistance in gravel-bed rivers." In Gravel-Bed Rivers, R.D. Hey, J.C. Bathurst, and CR. Thorne (Eds.), John Wiley & Sons, Chichester, 109-137. Ill Chang, H.H., 1980: "Geometry of gravel streams." Journal ASCE, of the Hydraulics Division, 106 (9), 1443-1456. Chow, V.T., 1959: Open Channel Hydraulics. McGraw-Hill Book Company Inc., New York, 580p. Clark, J.J., and Wilcock, P.R., 2000: "Effects of land-use change on channel morphology in northeastern Puerto Rico." Geological Craig, R.F., 1992: Soil Mechanics. Society of America 112 (12), 1763-1777. Bulletin, Chapman & Hall, London, 427p. Dunne, T., and Leopold, L.B., 1978: "Channel changes." In Water in Environmental Planning, W.H. Freeman and Company, San Francisco, 687-710. Einstein, H.A., 1942: "Formulas fo the transportation of bed load." ASCE, Transactions, 107, 561-597. Einstein, H.A., 1950: "The bed-load funciton for sediment transport in open channel flows." Technical Bulletin, 1026, U.S. Department of Agriculture, Washington, D.C, 1-71. Gilbert, G.K., 1914: "The transportation of debris by running water." Professional Paper 86, U.S. Geological Survey, 263p. Gomez, B., and Church, M., 1989: "An assessment of bed load sediment transport formulae for gravel bed rivers." Water Graf, W.H., 1971: Hydraulics Resources of Sediment Research. Transport. 25 (6), 1161-1186. McGraw-Hill, New York, 139 p. Griffiths, G.A., 1981: "Flow resistance in coarse gravel bed rivers." Journal Hydraulics Division, ASCE, 107 (HY7), 899-918. Hammer, T.R., 1972: "Stream channel enlargement due to urbanization." Resources Research, Water 8 (6), 1530-1540. Henderson, F.M., 1966: Open Channel Flow. Prentice-Hall Canada, Hey, R.D., 1979: "Flow resistance in gravel-bed rivers." Journal Division, of the Inc., of the Toronto, 521p. Hydraulics ASCE, 105 (HY4), 365-379. Hey, R.D., 1982: "Gravel-bed rivers: form and processes." In Gravel-Bed Rivers, Hey, J.C. Bathurst, and CR. Thorne (Eds.), John Wiley & Sons, Chichester, 5-13. 112 R.D. Hey, R.D., and Thorne, C.R., 1986: "Stable channels with mobile gravel beds." Journal Hydraulic Engineering, of 112 (8), 671-689. Horton, R.E., 1933: "Separate roughness coefficients for channel bottom and sides." Engineering News Record, 111 (22), 652-653. Huang, H.Q., and Nanson, G.C., 2000: "Hydraulic geometry and maximum flow efficiency as products of the principle of least action." Earth Surface Processes 25, andLandforms, 1-16. Jarrett, R.D., 1984: "Hydraulics of high-gradient streams." Journal Engineering, of Hydraulic 110 (11), 1519-1539. Keulegan, G.H., 1938: "Laws of turbulent flow in open channels." Journal of Research the National Bureau of Standards, 21, Research Paper, RP1151, 707-741. of Knapp, D.J., I960: "Analysis of new bed-load transport data in terms of Meyer-PeterMuller, regime "theory" and H.A. Einstein functions." National Canada: Division of Mechaniccal Research Council of 1-19. Engineering, Kondolf, G.M., 1997: "Hungry water: Effects of dams and gravel mining on river channels." Environmental Management, 21 (4), 533-551. Lawler, D.M., Thorne, C.R., and Hooke, J.M., 1997: "Bank erosion and instability." In Applied Fluvial Geomorphology for River Engineering and Management, Thorne, C.R., Hey, R.D., and Newson, M.D., (Eds.). John Wiley and Sons, 137-172. Leopold, L.B., and Wolman, M.G., 1957: "River channel patterns: braided, meandering and straight." ProfesionalPaper 282-B, U.S. Geological Survey, 85p. Lindley, E.S., 1919: "Regime channels." Proceedings Punjab Mackin, J.H., 1948: "Concept of the graded river." Bulletin, Engineering Conference. U.S. Geological Survey, 59, 463-512. Macklin, M.G., 1996: "Fluxes and storage of sediment-associated heavy metals in f loodplain systems: assessment and river basin management issues at a time of rapid environmental change." In Floodplain Processes, Anderson, M.G., Walling, D.E., and Bates, P.D., (Eds.). John Wiley and Sons, Chichester, 441-460. 113 Macklin, M.G., and Lewin, J., 1997: "Channel, floodplain and drainage basin response to environmental change." In Applied Fluvial Geomorphology for River Engineering and Thorne, C.R., Hey, R.D., and Newson, M.D., (Eds.). John Wiley and Sons, 15- Management, 45. MacVicar, B. J . , 1999: predicting Application for a rational model of stream equilibrium for Unpublished M.A.Sc. thesis, University of British channel adjustments. Columbia, Vancouver, British Columbia, 203p. Magilligan, F.J., and McDowell, P.F., 1997: "Stream channel adjustments following elimination of cattle grazing." Journal of the American Water Resources Association, 33 (4), 867-877. Manning, R., 1891: "On the flow of water in open channels and pipes." Transactions, Institution of Civil Engineers of Ireland, 20, 161-207. McLean, D.6., 1985: "Sensitivity analysis of bedload equations." In: Canadian Civil Engineers Annual Conference, Society of SK., 1-15. Saskatoon, Meyer-Peter, E., and Muller, R., 1948: "Formulas for bed-load transport." Report on the IAHR, Stokholm, Sweden, 39-64. Second Meeting, Millar, R.G., 1991: Development of an analytical model of river response. Unpublished M.A.Sc. thesis, University of British Columbia, Vancouver, British Columbia, 136p. Millar, R.G., and Quick, M.C., 1993a: "An optimization approach to river adjustments." In Hydraulic Engineering '93 : Proceedings of the 1993 Conference, ASCE, 1635-1640. Millar, R.G., and Quick, M.C., 1993b: "Effect of bank stability on geometry of gravel rivers." Journal of Hydraulic Millar, R.G., 1994: river channels. Engineering, An optimization model ASCE, 119 (12), 1343-1363. for the development and response of alluvial Unpublished Ph.D. thesis, The University of British Columbia, Vancouver, British Columbia, 320p. Millar, R.G., and Quick, M.C., 1998: "Stable width and depth of gravel-bed rivers with cohesive banks." Journal of Hydraulic Engineering, 124 (10), 1005-1013. Millar, R.G., 1999: "Grain and form resistance in gravel-bed rivers." Journal Research, 37 (3), 303-312. 114 of Hydraulic Millar, R.G., 2000: "Influence of bank vegetation on alluvial channel patterns." Resources Research, Water 36 (4), 1109-1118. Parker, G., Klingeman, P.C., and McLean, D.G., 1982: "Bedload and size distribution in paved gravel-bed streams." Journal of the Hydraulics Division, ASCE, 108 (HY4), 544-571. Pizzuto, J.E., Hession, W.C, McBride, M., 2000: "Comparing gravel-bed rivers in paired urban and rural catchments of southeastern Pennsylvania." Geology, 28 (1), 79-82. Reid, I., Bathurst, J.C, Carting, P. A., Walling, D.E., and Webb, B.W., 1997: "Sediment erosion, transport and deposition." In Applied Engineering and Management, Fluvial Geomorphology for River Thorne, C.R., Hey, R.D., and Newson, M.D., (Eds.). John Wiley and Sons, 115-135. Rinaldi, M., and Johnson, P.A., 1997: "Stream meander restoration." Journal American Water Resources Association, 33 (4), 855-866. Rouse, H., 1965: "critical analysis of open channel Division, ASCE, of the res\s-far\ce." Journal of the Hydraulics 9 (HY4), 1-25. Streeter, V.L., and Wylie, E.B., 1981: Fluid Mechanics. 562p. McGraw-Hill Ryerson Limited, Thorne, C.R., 1982: "Processes and menchanisms of bank erosion." In Gravel-Bed Rivers, R.D. Hey, J.C. Bathurst, and CR. Thorne (Eds.), John Wiley A Sons, Chichester, 227-271. Thorne, C.R., 1997: "Channel types and morphological classification." In Applied Geomorphology for River Engineering and Management, Fluvial Thorne, C.R., Hey, R.D., and Newson, M.D., (Eds.). John Wiley and Sons, 175-222. Vanoni, V.A., 1975: Sedimentation engineering. ASCE Manuals and Reports on Engineering Practice, 745p. Werritty, A., 1997: "Short-term changes in channel stability." In Applied Geomorphology for River Engineering and Management, Fluvial Thorne, C.R., Hey, R.D., and Newson, M.D., (Eds.). John Wiley and Sons, 47-65. Wilcock, P.R., 2001: "Toward a practical method for estimating sediment-transport rates in gravel-bed rivers." Earth Surface Processes and Land forms, 26, 1395-1408. White, W.R., Bettress, R., and Paris, E., 1982: "An analytical approach to river regime." Journal of Hydraulic Engineering, 108 (10), 1179-1193. 115 Williams, &.?., 1970: "Manning formula - a misnomer?" Journal of the Hydraulics Division, ASCE, 96 (HY1), 193-200. Wilson, A., 2001: "Predicting channel change in Narrowlake Creek in the Central Interior: A tool for watershed protection and restoration." Watershed Bulletin: Streamline, Restoration Technical 5 (4), 8-12. Wilson, A., Nesbit, B. and Pillipow, R., 2001: "Riparian restoration at the Narrowlake Creek Demonstration Watershed." Streamline, Watershed 6 (3), 24-29. 116 Restoration Technical Bulletin:
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- A visual model for predicting stream response of alluvial...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
A visual model for predicting stream response of alluvial gravel-bed rivers Moscote, Nerissa Yarmilla 2002
pdf
Page Metadata
Item Metadata
Title | A visual model for predicting stream response of alluvial gravel-bed rivers |
Creator |
Moscote, Nerissa Yarmilla |
Date Issued | 2002 |
Description | In this thesis, a visual optimization model, RiverMod, has been developed to predict the channel geometry of alluvial gravel-bed rivers. The model is an extension of a previous model developed by Millar and Quick (1993b). The model is based on equilibrium theory, and also includes a bank stability analysis. The adjustment of the dependent variables of a river reach is quantified according to changes incurred by the independent variables. The primary dependent variables are channel width, depth and slope, and the primary independent variables are the discharge, flow resistance parameters, sediment size distribution and bank stability parameters. RiverMod, written in Visual Basic, provides a user interface that promotes clear and simple model usage for multiple runs with different sediment transport and flow resistance equations. The theory behind the Millar and Quick (1993b) model is discussed, as well as Millar's (2000) meandering-braiding transition and its incorporation into RiverMod. Four sediment transport equations and four flow resistance equations are included as part of the model. The model consists of fixed-channel-slope and variable-channel slope versions. The fixed slope version is equivalent to an experiment where the slope is fixed, and the channel width, depth and sediment transport rate adjust to the discharge. The variable slope version more closely approximates natural stream conditions. Each of these models can be applied to streams with cohesive or non-cohesive bank sediment. The model is applied to two gravel-bed rivers, located in British Columbia. The bank stability of both rivers has been decreased due to logging along the banks. RiverMod is used to quantify the impact of such a disturbance by analyzing past and present channel conditions. Restoration methods are suggested based on model output. This thesis provides a full description of theory underlying the model development and a description of model usage; therefore, this thesis is a complete user manual for the RiverMod optimization model. |
Extent | 4464268 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Date Available | 2009-10-17 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0063942 |
URI | http://hdl.handle.net/2429/13956 |
Degree |
Master of Applied Science - MASc |
Program |
Civil Engineering |
Affiliation |
Applied Science, Faculty of Civil Engineering, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 2002-11 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-ubc_2003-0101.pdf [ 4.26MB ]
- Metadata
- JSON: 831-1.0063942.json
- JSON-LD: 831-1.0063942-ld.json
- RDF/XML (Pretty): 831-1.0063942-rdf.xml
- RDF/JSON: 831-1.0063942-rdf.json
- Turtle: 831-1.0063942-turtle.txt
- N-Triples: 831-1.0063942-rdf-ntriples.txt
- Original Record: 831-1.0063942-source.json
- Full Text
- 831-1.0063942-fulltext.txt
- Citation
- 831-1.0063942.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0063942/manifest