sd_channel_sediment3 Subroutine

subroutine sd_channel_sediment3()

Uses

  • proc~~sd_channel_sediment3~~UsesGraph proc~sd_channel_sediment3 sd_channel_sediment3 module~basin_module basin_module proc~sd_channel_sediment3->module~basin_module module~channel_module channel_module proc~sd_channel_sediment3->module~channel_module module~channel_velocity_module channel_velocity_module proc~sd_channel_sediment3->module~channel_velocity_module module~climate_module climate_module proc~sd_channel_sediment3->module~climate_module module~gwflow_module gwflow_module proc~sd_channel_sediment3->module~gwflow_module module~hru_module hru_module proc~sd_channel_sediment3->module~hru_module module~hydrograph_module hydrograph_module proc~sd_channel_sediment3->module~hydrograph_module module~reservoir_module reservoir_module proc~sd_channel_sediment3->module~reservoir_module module~sd_channel_module sd_channel_module proc~sd_channel_sediment3->module~sd_channel_module module~time_module time_module proc~sd_channel_sediment3->module~time_module module~utils utils proc~sd_channel_sediment3->module~utils module~water_body_module water_body_module proc~sd_channel_sediment3->module~water_body_module module~hydrograph_module->module~basin_module module~hydrograph_module->module~time_module

calculate channel sed and nutrient processes if inflow > 0 peak daily flow rate - m3/s using the Fuller and Peck (1974) equation for peak to mean flow ratio interpolate rating curve using peak rate channel velocity from the interpolated rating curve (rating-curve vel avoids the unrealistically high velocities at low flow from continuity Q/A) add precip to inflow - km * m * 1000 m/km * ha/10000 m2 = ha m3 = 10. * mm * ha compute overbank flow and flood plain deposition overbank volume is assumed = flood plain volume at peak rate (from the rating curve) trap efficiency from Dynamic SedNet Component Model Reference Guide: Update 2017 deposit Particulate P and N in the floodplain trap nitrate and sol P in flood plain - when not simulating flood plain interactions? if flood plain link - fill wetlands to emergency fill wetland storage to emergency volume calculate overbank volume left to fill wetland calc bank erosion calculate radius of curvature bank erosion m/yr

write for Peter

mass of sediment eroded -> t = 1000 * bankcut (mm) * depth (m) * lengthcut (m) * bd (t/m3) arc length = 0.33 * meander wavelength * sinuosity

calculate associated nutrients calculate channel deposition based on fall velocity - SWRRB book assume particle size = 0.03 mm -- median silt size

assume bankfull flow depth

calculate channel deposition as the bedload fraction of bank erosion calculate bed erosion no downcutting below equilibrium slope calc critical shear and shear on bottom of channel Pa = N/m^2 * m * m/m critical shear function of d50 bed erosion m/yr calc mass of sediment eroded -> t = m * width (m) * length (km) * 1000 m/km * bd (t/m3) calculate associated nutrients set outputs for sediment budget width and depth at end of the day - m sediment flood plain deposition - tons and mm mm = t / (5.bdwl) -> assume fp width = 5chw; len(m)=1000.km; bd=1.0 t/m3; mm=1000.m daily bank erosion in t and m - calculate w/yr later daily bed erosion in t and m - calculate w/yr later add to ordered sediment budget basin sediment budget outputs - don't sum during skip years

Arguments

None

Calls

proc~~sd_channel_sediment3~~CallsGraph proc~sd_channel_sediment3 sd_channel_sediment3 amax1 amax1 proc~sd_channel_sediment3->amax1 gwflow_floodplain gwflow_floodplain proc~sd_channel_sediment3->gwflow_floodplain proc~exp_w exp_w proc~sd_channel_sediment3->proc~exp_w rcurv_interp_flo rcurv_interp_flo proc~sd_channel_sediment3->rcurv_interp_flo tracebackqq tracebackqq proc~exp_w->tracebackqq