MATHEMATICAL MODELING OF FLOOD PROPAGATION IN RIVER BASINS

Background

Flooding among most destructive recurrent natural hazards affecting Nigeria Niger-Benue principal drainage network producing severe seasonal flooding repeatedly displaced millions extensive damage farmland infrastructure property most notably 2012 and 2022 flood events. Propagation flood wave along river channel from point peak rainfall-generated runoff through successive downstream reaches to vulnerable settlements fundamentally mathematical problem in theory hyperbolic and parabolic PDEs and effective early warning depends critically on ability to forecast with adequate lead time and accuracy how given upstream hydrograph will transform as travels downstream. Mathematical description unsteady open-channel flow originates Saint-Venant equations coupled nonlinear PDEs expressing conservation mass and momentum for one-dimensional shallow water flow first derived 1871. While full Saint-Venant captures complete dynamics including backwater effects and downstream boundary influences computationally demanding and for many practical flood routing applications can be simplified without significant loss to kinematic wave approximation in which momentum equation reduced to balance between gravitational and frictional forces yielding single first-order hyperbolic PDE relating discharge directly to channel rating curve. Kinematic wave equation while mathematically elegant and computationally efficient has well-known limitation because purely advective predicts flood wave propagates without reduction in peak discharge prediction inconsistent with attenuation and dispersion observed in genuine floods which arises physically from diffusive effects longitudinal pressure gradients and channel storage not captured. Muskingum-Cunge method diffusion-wave routing technique developed by Cunge 1969 reformulating classical Muskingum with physically based channel-derived parameters addresses limitation by incorporating numerical diffusion term calibrated to match physical diffusivity of full dynamic wave equation and remains among most widely used flood routing methods worldwide. Both require numerical solution for all but simplest cases and reliable application depends on solid understanding numerical stability and convergence properties of finite difference schemes governed for hyperbolic kinematic wave by CFL condition first identified Courant Friedrichs Lewy 1928 foundational analysis finite difference schemes hyperbolic PDEs.

Flood modeling and water resources project topics | External: USGS - Flood Routing, NOAA - Muskingum-Cunge Routing, Nigeria Hydrological Services Agency

Statement of Problem

Despite well-established mathematical theory kinematic and diffusion wave flood routing many applied treatments in Nigerian undergraduate mathematics research literature present flood routing methods most commonly Muskingum method in purely empirical calibration-based form without deriving underlying PDE without establishing numerical stability conditions governing chosen solution method and without formally comparing physically important distinction between non-attenuating kinematic wave prediction and attenuating diffusion-wave Muskingum-Cunge prediction. Limits both mathematical rigour and practical decision-relevant completeness for genuine flood forecasting application. Furthermore quantitative flood routing studies explicitly applied to Nigerian river basin conditions using validated numerical methods and physically grounded channel and roughness parameters and explicitly reporting engineering-relevant metrics such as peak attenuation percentage and time-to-peak lag at multiple downstream stations useful for setting flood early warning lead times less commonly presented within undergraduate mathematics tradition which more often addresses flood-related topics through purely statistical frequency analysis rather than mechanistic governing-equation-based modelling. Study addresses both gaps by formulating numerically implementing and rigorously validating kinematic wave and Muskingum-Cunge methods against exact analytical solutions and formal stability theory explicitly quantifying and explaining physical distinction between non-attenuating and attenuating routing predictions and applying validated framework to genuinely applied quantitatively specified flood routing problem along reach representative of Nigerian river system. Problem stated concisely: how can finite difference and diffusion-wave methods for flood routing be formulated numerically validated and applied to forecast flood wave attenuation and travel time along Nigerian river reach in manner directly useful for flood early warning system design?

Aim and Objectives

·         Formulate linear and nonlinear kinematic wave equations for flood routing and derive exact analytical method-of-characteristics solution linear case for validation

·         Implement explicit upwind and Lax-Wendroff finite difference schemes for kinematic wave equation and validate each against exact analytical solution achieving errors 2.32 and 0.011 m3/s at Cr 0.144

·         Derive CFL stability criterion for upwind scheme Cr = c dt/dx <=1 and demonstrate computationally using localized discharge pulse stable diffusive at 0.80 but unbounded oscillatory growth >20000 times initial within 3h at 1.15

·         Conduct formal convergence study using methodology avoiding single-point tracking artifacts to empirically verify order accuracy schemes 0.79 and 1.26 vs theoretical 1 and 2

·         Apply nonlinear kinematic wave equation formulated using Manning's rating curve to demonstrate quantify characteristic downstream steepening rising limb and estimate kinematic shock formation time

·         Implement Muskingum-Cunge diffusion routing method and directly compare physically realistic peak attenuation against non-attenuating prediction pure kinematic wave 45.4% vs 10.4% over 20km

·         Conduct sensitivity analysis peak outlet discharge with respect to channel roughness and bed slope indices -0.47 and 0.50

·         Apply validated Muskingum-Cunge framework to original case study routing design flood hydrograph along 120km reach representative Niger-Benue reporting peak attenuation and time-to-peak lag

Research Questions

·         How accurately do explicit upwind and Lax-Wendroff schemes reproduce exact analytical solution linear kinematic wave?

·         Does upwind exhibit stable-to-unstable transition predicted by CFL theory as Courant number crosses critical 1 and can demonstrated computationally?

·         What empirically observed orders convergence upwind and Lax-Wendroff and how compare theoretical predictions?

·         How does flood wave rising limb steepen as propagates under nonlinear kinematic wave and at what estimated time kinematic shock form?

·         How does peak attenuation predicted by Muskingum-Cunge compare with non-attenuating prediction pure kinematic wave finite difference?

·         How sensitive is peak outlet discharge to channel roughness and bed slope?

·         What peak attenuation and time-to-peak lag does validated framework predict along 120km reach representative of major Nigerian river system and implications for early warning lead time?

Significance

Significant to several stakeholders. To numerical analysts and applied mathematics students provides rigorously validated treatment finite difference methods for hyperbolic PDE including formal CFL stability demonstration and convergence-order verification using methodology explicitly designed to avoid common pitfalls in numerical convergence studies complementing extending parabolic heat equation methodology developed in related work same research programme. To hydrologists water resources engineers flood management agencies operating within Niger-Benue and other Nigerian river basins demonstrates mathematically transparent numerically validated method for quantifying flood wave attenuation and travel time quantities of direct urgent relevance to design flood early warning systems whose effectiveness depends critically on lead time available between upstream detection and downstream impact. To Nigerian federal and state emergency management agencies including Nigeria Hydrological Services Agency case study results Chapter Four provide concrete quantitatively specified illustration how validated mathematical flood routing can inform design warning lead times along representative long river reach: 120km reach 7.75-hour peak travel time 0.53% attenuation directly relevant to upstream-to-downstream lead times.

Numerical methods project topics | Civil engineering hydraulics topics

Scope and Limitations

Restricted to one-dimensional flood routing along single river channel reach using kinematic wave equation in both linear and nonlinear form and Muskingum-Cunge diffusion routing method. Full dynamic Saint-Venant wave equations which additionally capture backwater effects from downstream boundary conditions tributary inflows and floodplain storage outside scope though diffusion-wave approximation embodied in Muskingum-Cunge captures leading-order attenuation behaviour of full dynamic wave for gradually varied flow conditions typical large lowland river reaches such as those considered applied case study Section 4.8. Channel geometry represented using simplified wide-rectangular cross-section with Manning's equation used derive discharge rating curve; more complex compound or irregular cross-sections not modelled. Applied case study Section 4.8 uses stylized 120km reach parameters chosen broadly representative lower Niger-Benue system rather than specific surveyed river reach with field-measured cross-sectional and roughness data given constraints undergraduate research project. Limitations: kinematic and diffusion wave approximations neglect backwater effects arising from downstream boundary conditions such as tidal influence confluence with another river or downstream reservoir or structure which can be significant in lower reaches large river systems including parts Niger-Benue. Channel geometry idealised wide rectangular; genuine channels have irregular often compound main channel plus floodplain cross-sections modifying discharge rating curve and consequently wave celerity and routing parameters. Manning roughness coefficient and channel geometry parameters drawn from representative literature ranges rather than field survey or gauged data for specific named Nigerian river reach. Applied case study inflow hydrograph synthetic idealised smooth design hydrograph rather than historically observed flood hydrograph from specific Nigerian gauging station. Tributary inflows along routed reach which can significantly alter downstream hydrograph shape and peak in genuine basins not modelled; study considers only propagation transformation single input hydrograph along isolated reach.

Operational Definitions

Flood routing: Process computing how flood hydrograph changes shape and timing as propagates along river channel - validated via method of characteristics exact travelling-wave.

Kinematic wave: Approximation unsteady open-channel flow in which momentum equation reduces to balance between gravity and friction yielding purely advective non-attenuating wave equation - error 2.32 upwind 0.011 Lax-Wendroff at Cr0.144.

Wave celerity: Speed at which given discharge value propagates along channel generally increasing with discharge for natural channel rating curves - used in Cr = c dt/dx.

Kinematic shock: Discontinuity forms in kinematic wave solution when faster-moving higher-discharge characteristics overtake slower-moving lower-discharge characteristics emitted earlier - formation time estimated 1.12h rising limb compression 29.9 min to 6.0 min after 20km.

Muskingum-Cunge: Diffusion-wave flood routing method computing physically realistic peak attenuation using channel-derived storage weighting parameters - 45.4% attenuation vs 10.4% numerical diffusion kinematic over 20km.

Peak attenuation: Percentage reduction in flood peak discharge as wave propagates upstream to downstream - 45.4% Muskingum-Cunge vs 10.4% kinematic 20km reach; 0.53% across 120km case.

Time-to-peak lag: Delay between time peak discharge upstream and time peak downstream - 7.75-hour peak travel time across 120km Niger-Benue representative reach directly relevant early warning lead times.

CFL condition: Courant-Friedrichs-Lewy stability criterion requiring Courant number Cr = c dt/dx not exceed scheme-dependent threshold for numerical scheme remain stable - upwind threshold Cr<=1 demonstrated stable diffusive at 0.80 unbounded oscillatory growth >20000x initial within 3h at 1.15 convergence orders 0.79 upwind 1.26 Lax-Wendroff vs theoretical 1 and 2.

Conclusion

Linear kinematic wave solved analytically via method characteristics yielding exact travelling-wave solution against which explicit upwind and Lax-Wendroff validated achieving max absolute errors 2.32 and 0.011 m3/s at Courant 0.144. Von Neumann-type CFL analysis established criterion Cr <=1 for upwind threshold demonstrated computationally using localized pulse stable diffusive at Cr 0.80 but unbounded oscillatory growth >20000 times initial within three hours at Cr 1.15. Convergence study using full spatial-profile comparison at fixed evaluation time avoiding single-point tracking artifacts confirmed empirical orders 0.79 and 1.26 for upwind and Lax-Wendroff respectively consistent ranking though moderately below theoretical first- and second-order accuracy discrepancy attributed finite rather than infinite smoothness realistic flood hydrograph shape functions. Nonlinear kinematic wave formulated using Manning rating curve applied to 20km validation reach reproduced characteristic downstream steepening rising limb predicted theory converging characteristics rising-limb duration compressing 29.9 min upstream to 6.0 min after 20km consistent analytically estimated shock formation time 1.12h. Because kinematic purely advective predicts no reduction peak discharge output compared against Muskingum-Cunge diffusion method incorporating physically realistic peak attenuation: over same 20km Muskingum-Cunge produced 45.4% peak attenuation compared with 10.4% attributable numerical diffusion for kinematic scheme demonstrating physical necessity diffusion rather than pure advection term for realistic forecasting. Sensitivity analysis found peak outlet discharge respond comparable magnitude opposite sign to channel roughness index -0.47 and bed slope index 0.50. Validated Muskingum-Cunge applied to original case routing design flood hydrograph over 120km reach representative Niger-Benue predicting 7.75-hour peak travel time and 0.53% peak attenuation across reach results directly relevant to design upstream-to-downstream flood early warning lead times. Concludes finite difference and diffusion routing methods when rigorously validated against analytical solutions and formal stability theory provide mathematically defensible tools for flood forecasting in Nigerian river basins and recommends incorporation into operational early warning systems.

FAQs

How accurate are upwind and Lax-Wendroff schemes for linear kinematic wave?

Validated against exact method-of-characteristics travelling-wave solution max absolute errors 2.32 m3/s upwind and 0.011 m3/s Lax-Wendroff at Courant 0.144; second-order Lax-Wendroff substantially more accurate.

What is CFL condition and was it demonstrated?

Criterion Cr = c dt/dx <=1 for upwind scheme; computational demonstration using localized pulse remained stable diffusive at Cr0.80 but produced unbounded oscillatory growth >20000 times initial disturbance within 3h at Cr1.15 demonstrating stable-to-unstable transition crossing threshold 1.

What are empirical convergence orders?

Convergence study using full spatial-profile comparison at fixed evaluation time avoiding single-point tracking artifacts confirmed 0.79 upwind and 1.26 Lax-Wendroff consistent ranking though moderately below theoretical 1 and 2 due finite not infinite smoothness realistic hydrograph shape functions.

How does nonlinear kinematic wave steepen rising limb?

Manning rating curve formulation applied to 20km validation reach reproduces characteristic downstream steepening rising limb: 29.9 min upstream boundary compressing to 6.0 min after 20km propagation consistent analytically estimated kinematic shock formation time 1.12h when faster higher-discharge characteristics overtake slower earlier ones.

How does Muskingum-Cunge compare to pure kinematic wave attenuation?

Kinematic purely advective predicts no peak reduction; over 20km Muskingum-Cunge produced 45.4% physically realistic peak attenuation vs 10.4% attributable to numerical diffusion for kinematic scheme demonstrating necessity diffusion term for realistic peak forecasting.

What is sensitivity to roughness and slope?

Peak outlet discharge responds comparable magnitude opposite sign to channel roughness sensitivity index -0.47 and bed slope 0.50 - increasing roughness reduces peak, increasing slope increases peak.

What does 120km Niger-Benue case predict?

Validated Muskingum-Cunge routing design flood hydrograph over 120km stylized reach representative Niger-Benue predicts 7.75-hour peak travel time and 0.53% peak attenuation across reach directly relevant to designing upstream-to-downstream early warning lead times.

What is difference between kinematic and diffusion wave?

Kinematic approximation reduces momentum equation to gravity-friction balance yielding purely advective non-attenuating wave; Muskingum-Cunge diffusion-wave reformulates Muskingum with physically based channel-derived parameters incorporating numerical diffusion calibrated to match physical diffusivity of full dynamic Saint-Venant wave capturing leading-order attenuation.

What are limitations of study?

Neglects backwater effects tidal confluence reservoir, idealised wide rectangular cross-section not compound main plus floodplain, Manning and geometry from literature ranges not field survey, synthetic design hydrograph not historical gauged, no tributary inflows, full dynamic Saint-Venant outside scope though diffusion approximation captures leading-order attenuation for gradually varied large lowland reaches.

Why is this useful for flood early warning?

Quantifies flood wave attenuation and travel time mathematically transparent numerically validated against analytical solutions and formal stability theory providing concrete lead times e.g., 7.75h across 120km useful for setting warning thresholds for NIHSA and emergency management agencies.