Skip to content

Channel Lining Design

Open Channel Lining Design

Select and size erosion protection for open channels using the FHWA HEC-15 tractive-force method. This guide covers applied and permissible shear stress, the four major flexible-lining families (rock, vegetative, RECP, gabion), bend-induced shear amplification, and the full design procedure — including a comparison of lining types and integration with Channel Analysis.

Channel lining design is the process of selecting and sizing erosion protection materials for open channels. Whenever water flows in a channel, it exerts shear forces on the channel bed and sides. If these forces exceed the resistance of the channel boundary material, erosion occurs — leading to channel widening, deepening, migration, and potential structural failure of adjacent infrastructure.

Erosion protection is critical for constructed channels (roadside ditches, stormwater conveyances, spillway chutes) and for natural channels downstream of hydraulic structures where flow conditions have been altered. The consequences of inadequate protection range from routine maintenance headaches to catastrophic failures that endanger life and property.

The FHWA HEC-15 tractive-force method is the industry-standard approach for flexible channel-lining design. It provides a rational, physics-based framework for comparing the shear stress that flowing water applies to the channel boundary against the permissible shear stress that a given lining material can withstand.

The tractive-force approach is conceptually simple: the lining is adequate when its resistance exceeds the force applied to it. Formally, this is expressed as a comparison of two shear stresses:

τp    SF×τd\tau_{p} \;\ge\; \mathrm{SF} \times \tau_{d}
Tractive-force design criterion

Where:

  • τd\tau_d = applied (design) shear stress from the flowing water (Pa or N/m²)
  • τp\tau_p = permissible shear stress of the lining material (Pa or N/m²)
  • SF\mathrm{SF} = safety factor (typically 1.0 to 1.5, depending on consequences of failure)

If the permissible shear stress exceeds the factored applied shear stress, the lining is adequate. If not, the designer must choose a more resistant lining, reduce the channel slope, or modify the channel geometry to reduce the applied shear.

  • Constructed channels conveying stormwater, irrigation water, or spillway flows.
  • Roadside ditches and median channels along highways.
  • Channels downstream of culvert outlets, energy dissipators, or dam spillways.
  • Stream-bank stabilisation at bridge abutments and other structures.
  • Any channel where the native material cannot resist the design-flow shear stress.
Consequence of FailureRecommended SF
Low (minor maintenance)1.0
Moderate (property damage)1.2 – 1.3
High (structural damage, public safety)1.3 – 1.5

When water flows in a channel under uniform-flow conditions, the gravitational component of the water weight acting along the channel bed creates a shear force on the channel boundary. This boundary shear stress varies across the cross-section: it is highest on the channel bed and lower on the side slopes.

The maximum shear stress occurs at the channel bed. For uniform flow, the bed shear stress is given by HEC-15 Equation 3.1:

τd  =  γdS0\tau_{d} \;=\; \gamma \cdot d \cdot S_{0}
HEC-15 bed shear stress (Eq. 3.1)

Where:

  • γ\gamma = unit weight of water = 9 810 N/m³ (at standard conditions)
  • dd = maximum flow depth in the channel (m)
  • S0S_0 = longitudinal bed slope of the channel (m/m)

This equation follows directly from the balance of forces on a control volume of water: the weight component along the channel (proportional to γ\gamma, dd, and S0S_0) is balanced by the shear force on the boundary. It is the foundation of all lining design calculations.

The shear stress on the channel side slopes is lower than the bed shear stress. HEC-15 applies a reduction factor K1K_1 that depends on the side-slope geometry:

τside  =  K1γdS0\tau_{\mathrm{side}} \;=\; K_1 \cdot \gamma \cdot d \cdot S_{0}
Side-slope shear stress

The factor K1K_1 is typically in the range 0.7 to 0.8 for trapezoidal channels. For a channel with a bottom-width-to-depth ratio (B/dB/d) of 4, K10.75K_1 \approx 0.75 — meaning the side-slope shear is about 75% of the maximum bed shear. K1K_1 increases toward 1.0 for very narrow channels and decreases for very wide channels.

B/d RatioK1K_1 (Side Shear Factor)
10.77
20.74
40.75
60.76
≥ 80.76 – 0.78

The HEC-15 methodology covers four major categories of flexible channel lining. Each has distinct performance characteristics, cost profiles, and application niches. The following subsections describe each lining type, its permissible-shear calculation, and practical design guidance.

Rock lining is the most widely used flexible erosion-protection method for channels with moderate to high shear stresses. It consists of a layer of graded stone (riprap) placed on the channel bed and sides, often over a geotextile filter fabric to prevent fine material from piping through the rock voids.

When to use:

  • Channels with applied shear stresses exceeding 50 Pa.
  • Steep channels (slopes > 2%) where vegetative linings fail.
  • Permanent protection with a design life of 50+ years.
  • Channels requiring immediate protection (no establishment period).
  • Spillway chutes, dam toes, and culvert outlet protection.

The permissible shear stress for rock lining is based on the Shields stability criterion, relating the critical shear to the stone size and submerged weight:

τp  =  F(γsγw)D50\tau_{p} \;=\; F^{*} \cdot (\gamma_{s} - \gamma_{w}) \cdot D_{50}
Shields permissible shear for rock lining

Where:

  • FF^{*} = Shields parameter (dimensionless stability coefficient)
  • γs\gamma_s = unit weight of rock (typically 25 500 – 26 500 N/m³)
  • γw\gamma_w = unit weight of water (9 810 N/m³)
  • D50D_{50} = median stone diameter (m) — the size for which 50% of the rock by weight is smaller
Rock TypeF* RangeTypical F*
Crushed angular0.047 – 0.100.082
Angular to sub-angular0.06 – 0.120.094
Rounded (river cobble)0.03 – 0.060.047

Rock shape affects the angle of repose (ϕ\phi), which is the steepest angle at which loose rock remains stable. This angle determines the side-slope shear-reduction factor.

Rock ShapeAngle of Repose (ϕ\phi)
Very angular, crushed40° – 42°
Angular38° – 40°
Sub-angular to sub-rounded34° – 38°
Rounded30° – 34°

The side-slope stability factor (K2K_2) reduces the permissible shear on side slopes to account for the gravitational component acting down the slope:

K2  =  1sin2θsin2ϕK_{2} \;=\; \sqrt{ 1 - \frac{\sin^{2}\theta}{\sin^{2}\phi} }
Side-slope stability factor

where θ\theta is the channel side-slope angle (from horizontal) and ϕ\phi is the angle of repose. The side-slope permissible shear becomes τp,side=K2τp\tau_{p,\mathrm{side}} = K_{2} \cdot \tau_{p}.

The roughness of rock lining varies with the ratio of flow depth to stone size. HEC-15 provides the following relationship:

n  =  αD501/6n \;=\; \alpha \cdot D_{50}^{1/6}
Manning's n for rock lining

where α0.0395\alpha \approx 0.0395 (SI units with D50D_{50} in metres). The roughness coefficient typically ranges from 0.025 for small gravel (D50=25D_{50} = 25 mm) to 0.060 for large riprap (D50=500D_{50} = 500 mm). Larger rock creates more flow resistance, so the same channel carries less water per unit area.

In South African practice, rock for erosion protection is classified by nominal stone diameter following SANRAL specifications. Common classes include:

ClassD50D_{50} (mm)Typical Application
Gravel25 – 75Low-shear channels, filter layers
Light riprap100 – 200Roadside ditches, minor channels
Medium riprap200 – 350Stormwater channels, culvert outlets
Heavy riprap350 – 500Spillway chutes, high-velocity channels
Very heavy riprap500 – 750Dam toes, extreme shear applications

Vegetative linings use established grass or other ground cover to protect the channel boundary. The root system binds the soil, and plant stems increase roughness to slow the flow and reduce shear on the underlying soil. Vegetation is the most common lining for low-gradient channels in residential and agricultural settings.

When to use:

  • Channels with applied shear stresses below approximately 100 Pa (well-established grass).
  • Low to moderate channel slopes (typically < 5%).
  • Applications where aesthetics and environmental integration are important.
  • Situations where long-term maintenance (mowing) is feasible.
  • Temporary erosion protection during construction phases.

The permissible shear for vegetative linings depends on both the soil type and the vegetation coverage factor. HEC-15 uses a combined approach:

τp  =  τp,soilCf\tau_{p} \;=\; \tau_{p,\mathrm{soil}} \cdot C_{f}
Vegetative lining permissible shear

Where:

  • τp,soil\tau_{p,\mathrm{soil}} = base permissible shear of the underlying soil
  • CfC_f = vegetation coverage factor (enhancement multiplier, typically 1.5 to 5.0)

HEC-15 classifies vegetation into five retardance classes (A through E) based on stem height and density. These classes determine the Manning’s nn value, which varies with flow depth:

ClassDescriptionStem Heightn Range
AVery tall, dense grass (e.g. weeping lovegrass)> 750 mm0.10 – 0.40
BTall grass (e.g. bermuda, tall fescue)300 – 600 mm0.06 – 0.20
CMedium grass (e.g. mowed bermuda)150 – 300 mm0.04 – 0.12
DShort grass (e.g. closely mowed lawn)50 – 150 mm0.03 – 0.08
EVery short, sparse vegetation< 50 mm0.02 – 0.05

The underlying soil type is critical because vegetation does not fully protect the soil during extreme events when stems are flattened. Base permissible shear stresses for bare soil:

Soil Typeτp,soil\tau_{p,\mathrm{soil}} (Pa)
Non-cohesive fine sand (D75<1.3D_{75} < 1.3 mm)1 – 3
Non-cohesive coarse sand / gravel3 – 10
Cohesive lean clay (PI = 10 – 20)5 – 10
Cohesive stiff clay (PI > 20)10 – 20

Vegetative linings require ongoing maintenance to remain effective:

  • Regular mowing to maintain the design retardance class.
  • Re-seeding of bare or damaged areas before the wet season.
  • Weed and invasive-species control.
  • Fertilisation and irrigation during droughts (especially during establishment).
  • Sediment removal where significant deposition occurs.

RECP Lining — Rolled Erosion Control Products

Section titled “RECP Lining — Rolled Erosion Control Products”

RECPs are manufactured products that are rolled out on the channel surface to provide immediate erosion protection. They range from temporary biodegradable blankets to permanent turf reinforcement mats (TRMs) that work in combination with established vegetation.

When to use:

  • Temporary protection during vegetation establishment (construction phase).
  • Channels with moderate shear stresses (20 – 300 Pa, depending on product type).
  • Steep slopes where grass seed alone would wash away before germination.
  • Environmentally sensitive areas where vegetation is the long-term goal.
  • Applications requiring immediate protection without heavy equipment access.
TypeDescriptionτp\tau_p (Pa)Lifespan
Open-weave textile (OWT)Jute, coir, or straw netting20 – 556 – 24 months
Erosion control blanket (ECB)Straw, coconut, or excelsior fibre in a degradable net45 – 1001 – 3 years
Turf reinforcement mat (TRM)Permanent UV-stable polymer matrix100 – 300+Permanent
TRM-W (vegetated TRM)TRM with fully established grass300 – 500+Permanent

RECP permissible shear depends on the product type and underlying soil. For unvegetated RECPs, the permissible shear is the product’s rated threshold. For vegetated RECPs (TRMs with established grass), the combined permissible shear is:

τp  =  τp,product+τp,veg\tau_{p} \;=\; \tau_{p,\mathrm{product}} + \tau_{p,\mathrm{veg}}
Combined permissible shear for vegetated RECPs

Where τp,product\tau_{p,\mathrm{product}} is the manufacturer’s tested permissible shear for the bare product, and τp,veg\tau_{p,\mathrm{veg}} is the additional shear resistance provided by the established vegetation. TRMs with fully established grass can achieve combined permissible shear values exceeding 300 Pa.

Gabion mattresses are flat, rectangular wire-mesh baskets filled with rock, placed as a continuous lining on the channel bed and sides. The wire mesh provides structural confinement that significantly increases the erosion resistance compared to loose riprap of the same stone size. Gabions also allow vegetation to grow through them, improving long-term aesthetics and ecological value.

When to use:

  • Steep channels with high velocity where loose riprap may be displaced.
  • Applications requiring thinner lining than equivalent loose riprap.
  • Sites where smaller, locally available stone is available (gabion confinement allows use of smaller rock).
  • Channels where ecological integration and habitat creation are priorities.
  • Curved channels where gabion mattresses can conform to the bank shape.

The permissible shear for gabion-mattress lining is based on both the mattress thickness and the stone-fill size. The HEC-15 formula is:

τp  =  0.0091(γsγw)(t+D50)\tau_{p} \;=\; 0.0091 \cdot (\gamma_{s} - \gamma_{w}) \cdot (t + D_{50})
Gabion mattress permissible shear

Where:

  • tt = gabion mattress thickness (m)
  • D50D_{50} = median stone-fill diameter (m)
  • γs\gamma_s = unit weight of rock fill (N/m³)
  • γw\gamma_w = unit weight of water (9 810 N/m³)
  • 0.0091 = empirical coefficient from gabion stability testing
Thickness (mm)Typical D50D_{50} (mm)Approx. τp\tau_p (Pa)Application
15075 – 10035 – 40Low shear, channel banks
230100 – 15055 – 65Moderate channels
300100 – 20070 – 85Major stormwater channels
500150 – 250100 – 120High-velocity channels, spillways
  • Stone must be hard, durable, and weather-resistant (not susceptible to slaking or freeze-thaw degradation).
  • D50D_{50} should be between 1.0 and 2.0 times the mesh opening size.
  • D85/D15D_{85}/D_{15} ratio should be between 1.5 and 3.5 for dense packing.
  • Minimum stone specific gravity of 2.5 (preferably > 2.6).
  • Wire mesh should be galvanised and PVC-coated for corrosive or marine environments.

Channel bends create secondary currents that redirect flow toward the outer bank, significantly increasing the local shear stress. This is one of the most common causes of channel-lining failure in practice. HEC-15 addresses this through a bend shear multiplier.

The bend multiplier increases the applied shear stress on the outer bank of a bend:

τd,bend  =  Kbτd\tau_{d,\mathrm{bend}} \;=\; K_{b} \cdot \tau_{d}
Bend-adjusted applied shear

The multiplier KbK_b depends on the ratio of the channel centreline bend radius (RcR_c) to the top width of flow (TT):

Rc/TR_c / TKbK_bBend Classification
< 22.0Very tight bend
2 – 31.8 – 2.0Tight bend
3 – 51.5 – 1.8Moderate bend
5 – 101.2 – 1.5Gentle bend
> 101.0Effectively straight

The increased shear from a bend persists for a significant distance downstream. HEC-15 recommends extending the bend protection for a length LpL_p downstream of the bend exit:

Lp  =  KbRcTL_{p} \;=\; K_{b} \cdot \frac{R_{c}}{T}
Downstream bend-protection length

For tight bends (Rc/T<3R_c/T < 3), the elevated shear may persist for 5 to 10 channel widths downstream. It is common practice to extend the heavier lining for at least this distance before transitioning to the standard lining.

Flow around a bend creates a water surface that is higher on the outside and lower on the inside. The super-elevation (Δy\Delta y) must be considered for freeboard calculations:

Δy  =  V2TgRc\Delta y \;=\; \frac{V^{2} \cdot T}{g \cdot R_{c}}
Super-elevation at bends

Where:

  • VV = average flow velocity (m/s)
  • TT = top width of flow (m)
  • gg = gravitational acceleration (9.81 m/s²)
  • RcR_c = centreline bend radius (m)

The following step-by-step procedure outlines the complete channel-lining design process using the HEC-15 tractive-force method. The Channel Lining Design calculator automates these steps, but understanding the procedure is essential for interpreting and validating results.

  1. Define channel geometry and flow conditions. Establish the channel cross-section (shape, dimensions, side slopes), longitudinal bed slope, and design flow rate. The design flow is typically derived from a hydrological analysis for a specified return period (e.g. 1:50 year for major stormwater channels). Use Channel Analysis to solve for flow depth and velocity.

  2. Calculate applied shear stress. Compute the maximum bed shear (τd=γdS0\tau_d = \gamma \cdot d \cdot S_0) and the side-slope shear (τside=K1τd\tau_{\mathrm{side}} = K_1 \cdot \tau_d). If the channel has bends, apply the bend multiplier KbK_b to the outer-bank shear.

  3. Select lining type. Based on the magnitude of applied shear stress, project constraints (cost, aesthetics, maintenance), and material availability, select a candidate lining type: rock, vegetation, RECP, or gabion. Use the comparison table below as a starting point.

  4. Calculate permissible shear stress. Using the appropriate formula for the selected lining type, compute the permissible shear stress. For rock: use the Shields equation with the chosen D50D_{50}. For vegetation: combine the soil base shear with the coverage factor. For RECP: use the manufacturer’s tested value. For gabion: use the mattress-thickness formula.

  5. Apply safety factor. Multiply the applied shear stress by the safety factor appropriate for the consequences of failure (see the Safety Factor Selection table above). The factored applied shear is the demand that the lining must resist.

  6. Check adequacy. Verify that τpSFτd\tau_p \ge \mathrm{SF} \cdot \tau_d for both the bed and side slopes. Check both the bed lining and side-slope lining independently, as the controlling condition may differ. The design ratio (τp/(SFτd)\tau_p / (\mathrm{SF} \cdot \tau_d)) should be at least 1.0; values between 1.0 and 1.3 are considered efficient.

  7. Iterate if needed. If the lining is inadequate, increase the stone size (rock), use a thicker mattress (gabion), select a higher-rated product (RECP), or change the lining type entirely. If the lining is excessively over-designed (design ratio > 2.0), consider downsizing for cost savings. Repeat the check.

  8. Consider bend effects. If the channel has bends, check the outer-bank shear with the bend multiplier KbK_b. A heavier lining may be required for bend sections. Determine the protection length downstream of each bend and check freeboard for super-elevation.

  9. Determine freeboard. Add freeboard above the design water surface to account for wave action, flow surges, debris, and super-elevation at bends. Typical freeboard is 0.15 to 0.60 m depending on channel size and the consequences of overtopping. Extend the lining to above the freeboard level.

The table below provides a side-by-side comparison of the four flexible-lining types across key design criteria. Use this as a starting point for lining selection, then refine based on project-specific constraints.

CriterionRock RiprapVegetativeRECPGabion
τp\tau_p range (Pa)50 – 500+5 – 10020 – 30035 – 120
Manning’s n range0.025 – 0.0600.020 – 0.4000.020 – 0.0450.025 – 0.040
Relative costMedium – HighLowLow – MediumMedium – High
Durability50+ yearsIndefinite (with maintenance)6 months – permanent25 – 50 years
AestheticsNatural stone appearanceExcellent (green channel)Good (once vegetated)Good (vegetation through mesh)
Ecological valueModerate (habitat in voids)HighModerate – HighHigh (diverse habitat)
InstallationHeavy equipment neededSeeding / soddingManual roll-out + anchoringManual assembly + fill
MaintenanceLow (periodic inspection)High (mowing, reseeding)Low – MediumLow (wire inspection)
Immediate protectionYesNo (6 – 12 month establishment)YesYes

The Channel Lining Design calculator integrates with the Channel Analysis tool. This allows you to define your channel hydraulics first, then import the results directly into the lining-design workflow.

  1. Channel Analysis: Define your channel geometry, Manning’s nn, and slope. Solve for normal depth at your design flow rate. Save the calculation.
  2. Import: In the Channel Lining Design calculator, use the import feature to pull in the channel geometry, flow depth, velocity, and slope from your saved Channel Analysis calculation.
  3. Lining design: With the hydraulic parameters pre-filled, select your lining type, specify lining parameters (stone size, vegetation class, etc.), and run the shear-stress comparison.
  4. Iterate: If the lining is inadequate, adjust parameters and re-check without returning to the Channel Analysis tool. The imported hydraulic data remains available for reference.
  • Channel Analysis — compute the uniform-flow depth, velocity, and bed shear that feed directly into the lining design.
  • GVF Profiles — for non-uniform reaches where the shear stress varies along the channel.
  • Weir Analysis — for backwater-producing structures that increase upstream lining demand.
  • Culvert Designer — for inlet and outlet protection where roadway channels transition into culverts.
  • FHWA HEC-15: Kilgore, R.T. and Cotton, G.K. (2005). Design of Roadside Channels with Flexible Linings. Hydraulic Engineering Circular No. 15, Third Edition. Federal Highway Administration, Publication No. FHWA-NHI-05-114.
  • FHWA HDS-4: Schall, J.D., Richardson, E.V., and Morris, J.L. (2008). Introduction to Highway Hydraulics. Hydraulic Design Series No. 4, Fourth Edition. Federal Highway Administration.
  • USACE. (1994). Hydraulic Design of Flood Control Channels. Engineer Manual EM 1110-2-1601. US Army Corps of Engineers, Washington DC. Complementary guidance on rigid and flexible channel design.
  • Chow, V.T. (1959). Open-Channel Hydraulics. McGraw-Hill, New York. Classical derivation of the tractive-force concept (Chapter 7).
  • SANRAL Drainage Manual: South African National Roads Agency Limited (2013). Drainage Manual, Sixth Edition. Chapter 8 — Erosion Protection of Open Channels.
  • Shields, A. (1936). Application of Similarity Principles and Turbulence Research to Bed-Load Movement. California Institute of Technology, Pasadena (translated from German). Original derivation of the incipient-motion stability parameter used for rock lining.
  • CIRIA C683: CIRIA, CUR, CETMEF (2007). The Rock Manual: The Use of Rock in Hydraulic Engineering. Second edition. CIRIA Publication C683. London.
  • ASTM D6460: Standard Test Method for Determination of Erosion Control Blanket (ECB) Performance in Protecting Earthen Channels from Stormwater-Induced Erosion. ASTM International, West Conshohocken, PA.

Open Channel Lining Design