170 likes | 322 Views
Kinematic Routing Model and its Parameters Definition. Routing Model. Real HRAP Cell. Hillslope model. Cell-to-cell channel routing. Separate Treatment of Fast and Slow Runoff. HRAP Cell. Hillslope Routing. Kinematic Wave Koren et al. (2004)
E N D
Routing Model Real HRAP Cell Hillslope model Cell-to-cell channel routing
Hillslope Routing • Kinematic Wave • Koren et al. (2004) • Independent routing for each hillslope element • Only routes fast runoff x Grid Pixel Conceptual Hillslope q = discharge per unit area of hillslope h = average overland flow depth Rs = fast runoff from water balance Sh = hillslope slope nh = hillslope roughness D = drainage density Lh = hillslope length Continuity: Momentum:
Channel Routing • Kinematic Wave • Koren et al. (2004) • Routes • fast runoff from hillslope • Slow runoff x Grid Pixel Q = channel discharge A = channel cross-sectional area qLh = overland flow rate at the hillslope outlet Rg = slow runoff component from the water balance Fc = grid cell area Lc = channel length within a cell Continuity: Momentum:
Kinematic Wave Advantages • Require few parameters • Easy to generate fast implicit numerical scheme compared, e.g., to diffusive model • Flexible in selection of simulation time-space increments • Allows selection of larger time increments compared to other models
Kinematic Wave Disadvantages • Lack of attenuation specifically for very flat basins • RDHM defines a simple one-shape channel cross-section • Potential effect on results comparing other models • When channel properties vary in space, attenuation may occur • Any numerical scheme introduces some attenuation. So use, e.g., diffusive model will accelerate physical attenuation by numerical • There are few criteria that allows estimation of potential errors • In headwater basins, wave form change depends mostly on lateral inflow contribution and joining channels, not attenuation • One-shape channel can affect significantly on simulation results if there is a flood plane. It might be difficult to get a reasonable peak timing for high and low floods
HL-RDHM Routing Parameters • There are three basic parameters • Hillslope depth-discharge relationship parameter, qs • Two parameters of channel discharge-cross-section relationship, qo and qm • Parameters have to be defined at each grid cell above selected basin • These parameters are not directly measurable • Combination of local basin properties (topography, soil, vegetation) and an integrated basin response at the outlet (discharge measurement information)
Hillslope Routing Parameter Derivation • Four hillslope property grids have to be defined • Surface slope, Sh • Manning’s roughness coefficient, nh • Channel density as a ratio of total channel length to area, D • Pixel area, f • HL-RDHM calculates the basic parameter qs at each grid cell during run-time from
Channel Routing Parameters Derivation • Two methods are available in RDHM • ‘Rating curve’ method that estimates the parameters q0 and qm directly using hydraulic measurements at an outlet gauging station • ‘Channel shape’ method assumes a simple parabolic channel geometry and uses outlet hydraulic measurements to estimate shape parameters at outlet • then basic routing parameters at outlet • Grids of Sc (channel slope) and nc (channel Manning’s roughness coefficient) should be available above outlet
R Scripts Provided to Assist with Flow Measurement Analysis Outletmeas_manual.R automatically generates several plots and computes reqressions User can specify plotting and regression weight options Directly to Q = q0*Aqm • Generate A =a*Bb • Q = v*A = q0*Aqm
Typical Channel Shape Depending on β b = 1 b = 0 b > 1 b < 1
Assumptions on Derivation Parametric Grids • Two assumptions from channel geometry laws are adopted for interpolation outlet parameters to upstream • The ratio of channel-forming flows at different cells equals ratio of drainage areas above the cells • The ratio of channel cross-sectional areas of different cells is a known function of stream orders • Also, parameters qm (rating curve method) and β (channel shape method) assumed to be constant above selected outlet
Generating Distributed Routing Parameters • Information needed • Parameters estimated at an outlet pixel • Drainage area • Connectivity • Geomorphologic relationship • Channel slope and roughness for the channel shape only method Rating Curve method Channel Shape Method Extrapolate q0 qm Estimate
Distribute Parameters Upstream using Genpar • Features of Genpar • Needs a base grid • Modifies the entire area upstream of an outlet • Able to handle multiple outlets Assign values to entire upstream area Overwrite values for sub-basins
1 Specific discharge grids generated using different number of outlets: 1,2,6 Arkansas river, USA 2 6