Civil Engineering Civil Engineering major 100% Free Open Access
Chapter 5 • Theory & Derivations

Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing

Comprehensive university textbook unit covering fundamental theory, mathematical derivations, computational hydraulics, and engineering design validation of Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing. Developed for undergraduate civil engineering study.

§5.1 Case study: Wetland flood peak attenuation

Case study: Wetland flood peak attenuation

A freshwater wetland receiving a constant 50 m³/s from an upstream channel initially releases 10 m³/s. A linear storage timescale of 5 hours predicts an increasing outflow but no immediate match with inflow. Downstream bridge operations depend on the peak time rather than just the total event volume.

Investigation task. Calculate outflow at 1, 3 and 6 hours and explain why an actual wetland stage-discharge function may not be linear.

The unit as part of a watershed investigation

This unit asks an engineer to defend storage response and routing as a quantitative argument rather than a number inserted into a generic formula. The natural control volume is the selected hydrologic subsystem and its actual inflow and outflow boundaries. Write the unit of every parameter and specify the temporal support of the observations; an hourly intensity, a daily depth and an annual volume cannot be substituted for each other without integration. A concept sketch should show the site boundaries, measurement stations, surface and subsurface pathways, relevant human infrastructure and the time interval over which the process operates. Interpret the final result relative to the real decision at wetland flood peak attenuation, not merely against an arbitrary calculator threshold.

Field evidence, scale and uncertainty

The most important practical complication in hydrographs, routing and storage is that model inputs are estimates of spatially and temporally varying states. Two instruments reporting the same variable may not have comparable averaging periods or calibration. The impact of an uncertain coefficient depends on the derivative of the output with respect to that coefficient and on where the operating point falls relative to a threshold. If a predicted response appears robust, deliberately test a neighboring process that the simple model omits: channel blockage, antecedent saturation, high tailwater, sensor drift, groundwater exchange or a different storm shape. Record failures of the model openly, rather than editing the inputs until a visually pleasing curve appears.

Professional inference standard

An educational simulation is useful when it generates a falsifiable prediction. Before using any of the four lessons to make a design claim, identify an independently measurable output, a second approximation with different assumptions and a failure consequence. For wetland flood peak attenuation, a field campaign might require time-synchronized precipitation and flow observations, a georeferenced inventory of hydraulic structures, separate dry- and wet-weather site visits and a provenance ledger for every parameter. The exercises below are intentionally synthetic and no parameter set is asserted to describe a real Bangladesh, US or other watershed. Local codes and agency-approved hydrologic methods supersede these teaching screens.

§5.2 Event hydrographs and baseflow separation

Lesson 5.1 — Event hydrographs and baseflow separation

Hydrological mechanism

A stream hydrograph combines storm-responsive discharge and slower sources. Separating a baseflow component is a modeling convention, not direct physical observation of individual water molecules. Peak discharge, time to peak, rising-limb slope, recession and cumulative volume offer complementary information. A hydrograph can have the correct peak but the wrong volume, which matters for downstream storage and floodplain decisions.

Governing expression and dimensions

Teaching relation: Event volume V=sum(Q_k Delta t)-baseflow volume. Every term must be traced to its physical definition before numerical evaluation. A hydrologic depth expressed in millimeters corresponds to one thousand cubic meters for every square kilometer of area. A discharge expressed in cubic meters per second becomes volume only after integrating over seconds. A rate reported over a nominal day or hour is not a spatial average by itself. State whether the equation is a conservation balance, idealization, calibrated empirical law or fitted statistical relation; that classification determines what evidence is needed before transfer to a new catchment.

Units and source audit. For event hydrographs and baseflow separation, list the observed quantities, estimated coefficients and chosen reference conditions separately. Include the dates and methods used to obtain each. Calculate one reference state using all quantities in compatible units and round only after a dimensional check. Then reverse a sign convention or convert to a different set of units and verify that the physical meaning survives. A unit-consistent equation may still be wrong because its assumptions do not match the setting.

Applied decision: Review an upstream and downstream gauge pair during a storm

The central decision is not whether a formula can be evaluated; it is whether the value carries information sufficient for an engineering action. In this situation, first identify the outlet, observation period and decision variable. Establish a low, central and high forcing scenario using values that can be defended by measurements or explicitly labeled hypotheses. Predict how the response changes before running the model. A prediction that contradicts conservation of water, plausible timing or an independently measured constraint requires investigation even when the graph looks smooth.

For this particular case the dominant caution is that clock alignment and missing recession tails can dominate integrated volume errors. Assess which boundary condition is most likely to be violated. For an event analysis, compare the timing of inputs and outputs; for a seasonal balance, separate fluxes from storage; for statistical estimates, distinguish frequency assumptions from observed counts. Do not claim a coefficient calibrated in one environment is universal. State one missing measurement and the practical difficulty of obtaining it.

Mechanism checks and falsification

Counterfactual A — scale. Suppose the area or response time is changed while keeping some input depth or intensity fixed. Decide whether the output should scale linearly, sublinearly or nonlinearly, and explain which physical assumption controls that conclusion. For event hydrographs and baseflow separation, the relevant limiting cases include a nearly dry condition, a very large forcing and a coefficient approaching the edge of its documented range. A model that produces negative water depth, impossible storage or negative coefficients without a clear warning is not fit for student use.

Counterfactual B — observation. Ask how a stage recorder, pluviometer, soil-moisture probe or properly maintained flow measurement could disprove the chosen interpretation. Compare two alternatives: an error in the input and an error in the process structure. Error in a timestamp may alter the inferred lag without changing storm totals; an unmodeled flow diversion may alter volume without a rainfall error. Use these fingerprints to choose the next data collection action.

Cross-topic connection. This lesson connects to linear reservoir response because a model transfers its assumptions through the wider analysis. It also affects lakes, reservoirs and flood storage, which may rely on the same boundary records but at a different temporal or spatial resolution. A strong report shows these dependencies explicitly, identifies the quantities shared by two calculations and avoids treating correlated errors as independent measurements.

Field, laboratory and spreadsheet workflow

  1. Draw a control-volume diagram for review an upstream and downstream gauge pair during a storm; mark every measured and assumed boundary.
  1. Calculate the result from Event volume V=sum(Q_k Delta t)-baseflow volume where its inputs are available. Show the time-base and area conversions with intermediate numbers rather than only the software answer.
  1. Perform a one-at-a-time sensitivity screen for the most uncertain input. State which physical relationship predicts the direction of the response; explain any threshold or plateau.
  1. Design an independent field check. Compare peak-fit and volume-fit models using the same observations. Explain exactly what would count as disagreement and how you would respond before changing model parameters.
  1. Write a limitation paragraph and a separate decision paragraph. The limitation describes what has not been established; the decision identifies the next action that is justified by the evidence available.

Worked investigative thought experiment

Consider two plausible observations from review an upstream and downstream gauge pair during a storm. In the first, the measured quantity changes noticeably after a storm, whereas in the second, a neighboring instrument records a muted response. One explanation is spatial variability in forcing; another is storage and travel-time filtering. A third is mismatched sampling intervals. Quantify the first explanation with an appropriately scaled synthetic input, then test the second through a changed process coefficient. Keep the third as a data-quality hypothesis until the metadata are checked. This separation protects against false certainty: a model that can fit both instruments with different parameters has not identified the underlying cause.

In-class questions and solutions guidance

Question 1. In the context of review an upstream and downstream gauge pair during a storm, identify all inputs and outputs of the selected hydrologic control volume. Answer guidance: give measured depths, rates, storage changes and known management transfers; justify the time interval and boundary instead of simply listing variables.

Question 2. Under what conditions is Event volume V=sum(Q_k Delta t)-baseflow volume defensible? Answer guidance: specify the governing conservation or calibration assumptions, a physical range of inputs and at least one failure condition associated with clock alignment and missing recession tails can dominate integrated volume errors.

Question 3. How could you distinguish genuine process sensitivity from measurement bias? Answer guidance: formulate competing mechanisms, specify an independent sensor or additional event, and compare timing, volume and directional response.

Question 4. What changes if the same equation is transferred from a small engineered catchment to a large heterogeneous river basin? Answer guidance: evaluate spatial averaging, storage, travel time, changing active area, scale dependence and available calibration evidence.

Extended engineering reflection

A credible result for event hydrographs and baseflow separation needs more than an apparently precise scalar. Provide a unit table, parameter source table, version of the governing equation, uncertainty register and plot of response against at least one input. Compare the output to a basic conservation or range check; no negative runoff depth can be created from positive water inputs without explaining a signed flux convention. Consider whether the decision is reversible, such as adjusting a monitoring schedule, or hard to reverse, such as constructing a hydraulic structure. For high-consequence decisions, missing data should lead to targeted measurement and professional review rather than a higher number of decimal places.

§5.3 Linear reservoir response

Lesson 5.2 — Linear reservoir response

Hydrological mechanism

A linear reservoir relates outflow to storage with Q=S/K and dS/dt=I-Q. Under a constant inflow rate, outflow approaches inflow exponentially over characteristic time K. The model is useful for illustrating delay and attenuation but not for backwater-affected floodplains or structures with complicated elevation–discharge curves. Initial conditions can matter as much as the forcing when the modeled interval is short.

Governing expression and dimensions

Teaching relation: Q(t)=I+(Q0-I) exp(-t/K) for constant I. Every term must be traced to its physical definition before numerical evaluation. A hydrologic depth expressed in millimeters corresponds to one thousand cubic meters for every square kilometer of area. A discharge expressed in cubic meters per second becomes volume only after integrating over seconds. A rate reported over a nominal day or hour is not a spatial average by itself. State whether the equation is a conservation balance, idealization, calibrated empirical law or fitted statistical relation; that classification determines what evidence is needed before transfer to a new catchment.

Units and source audit. For linear reservoir response, list the observed quantities, estimated coefficients and chosen reference conditions separately. Include the dates and methods used to obtain each. Calculate one reference state using all quantities in compatible units and round only after a dimensional check. Then reverse a sign convention or convert to a different set of units and verify that the physical meaning survives. A unit-consistent equation may still be wrong because its assumptions do not match the setting.

Applied decision: Explore wetland detention under a steady inflow pulse

The central decision is not whether a formula can be evaluated; it is whether the value carries information sufficient for an engineering action. In this situation, first identify the outlet, observation period and decision variable. Establish a low, central and high forcing scenario using values that can be defended by measurements or explicitly labeled hypotheses. Predict how the response changes before running the model. A prediction that contradicts conservation of water, plausible timing or an independently measured constraint requires investigation even when the graph looks smooth.

For this particular case the dominant caution is that the single storage coefficient cannot capture floodplain thresholds and reverse flow. Assess which boundary condition is most likely to be violated. For an event analysis, compare the timing of inputs and outputs; for a seasonal balance, separate fluxes from storage; for statistical estimates, distinguish frequency assumptions from observed counts. Do not claim a coefficient calibrated in one environment is universal. State one missing measurement and the practical difficulty of obtaining it.

Mechanism checks and falsification

Counterfactual A — scale. Suppose the area or response time is changed while keeping some input depth or intensity fixed. Decide whether the output should scale linearly, sublinearly or nonlinearly, and explain which physical assumption controls that conclusion. For linear reservoir response, the relevant limiting cases include a nearly dry condition, a very large forcing and a coefficient approaching the edge of its documented range. A model that produces negative water depth, impossible storage or negative coefficients without a clear warning is not fit for student use.

Counterfactual B — observation. Ask how a stage recorder, pluviometer, soil-moisture probe or properly maintained flow measurement could disprove the chosen interpretation. Compare two alternatives: an error in the input and an error in the process structure. Error in a timestamp may alter the inferred lag without changing storm totals; an unmodeled flow diversion may alter volume without a rainfall error. Use these fingerprints to choose the next data collection action.

Cross-topic connection. This lesson connects to muskingum routing and numerical stability because a model transfers its assumptions through the wider analysis. It also affects event hydrographs and baseflow separation, which may rely on the same boundary records but at a different temporal or spatial resolution. A strong report shows these dependencies explicitly, identifies the quantities shared by two calculations and avoids treating correlated errors as independent measurements.

Field, laboratory and spreadsheet workflow

  1. Draw a control-volume diagram for explore wetland detention under a steady inflow pulse; mark every measured and assumed boundary.
  1. Calculate the result from Q(t)=I+(Q0-I) exp(-t/K) for constant I where its inputs are available. Show the time-base and area conversions with intermediate numbers rather than only the software answer.
  1. Perform a one-at-a-time sensitivity screen for the most uncertain input. State which physical relationship predicts the direction of the response; explain any threshold or plateau.
  1. Design an independent field check. Vary K and identify the moment when lag becomes operationally important. Explain exactly what would count as disagreement and how you would respond before changing model parameters.
  1. Write a limitation paragraph and a separate decision paragraph. The limitation describes what has not been established; the decision identifies the next action that is justified by the evidence available.

Worked investigative thought experiment

Consider two plausible observations from explore wetland detention under a steady inflow pulse. In the first, the measured quantity changes noticeably after a storm, whereas in the second, a neighboring instrument records a muted response. One explanation is spatial variability in forcing; another is storage and travel-time filtering. A third is mismatched sampling intervals. Quantify the first explanation with an appropriately scaled synthetic input, then test the second through a changed process coefficient. Keep the third as a data-quality hypothesis until the metadata are checked. This separation protects against false certainty: a model that can fit both instruments with different parameters has not identified the underlying cause.

In-class questions and solutions guidance

Question 1. In the context of explore wetland detention under a steady inflow pulse, identify all inputs and outputs of the selected hydrologic control volume. Answer guidance: give measured depths, rates, storage changes and known management transfers; justify the time interval and boundary instead of simply listing variables.

Question 2. Under what conditions is Q(t)=I+(Q0-I) exp(-t/K) for constant I defensible? Answer guidance: specify the governing conservation or calibration assumptions, a physical range of inputs and at least one failure condition associated with the single storage coefficient cannot capture floodplain thresholds and reverse flow.

Question 3. How could you distinguish genuine process sensitivity from measurement bias? Answer guidance: formulate competing mechanisms, specify an independent sensor or additional event, and compare timing, volume and directional response.

Question 4. What changes if the same equation is transferred from a small engineered catchment to a large heterogeneous river basin? Answer guidance: evaluate spatial averaging, storage, travel time, changing active area, scale dependence and available calibration evidence.

Extended engineering reflection

A credible result for linear reservoir response needs more than an apparently precise scalar. Provide a unit table, parameter source table, version of the governing equation, uncertainty register and plot of response against at least one input. Compare the output to a basic conservation or range check; no negative runoff depth can be created from positive water inputs without explaining a signed flux convention. Consider whether the decision is reversible, such as adjusting a monitoring schedule, or hard to reverse, such as constructing a hydraulic structure. For high-consequence decisions, missing data should lead to targeted measurement and professional review rather than a higher number of decimal places.

§5.4 Muskingum routing and numerical stability

Lesson 5.3 — Muskingum routing and numerical stability

Hydrological mechanism

Muskingum routing uses storage S=K[X I+(1-X) O] as an empirical representation of channel storage. In a discrete step, O_next=C0 I_next+C1 I_now+C2 O_now, with coefficients depending on K, X and Delta t. For nonnegative coefficients in the conventional scheme, 2 K X <= Delta t <= 2 K(1-X), with 0<=X<=0.5. This stability/positivity range is not proof of fidelity to a specific channel.

Governing expression and dimensions

Teaching relation: C0=(Delta t-2KX)/(2K(1-X)+Delta t). Every term must be traced to its physical definition before numerical evaluation. A hydrologic depth expressed in millimeters corresponds to one thousand cubic meters for every square kilometer of area. A discharge expressed in cubic meters per second becomes volume only after integrating over seconds. A rate reported over a nominal day or hour is not a spatial average by itself. State whether the equation is a conservation balance, idealization, calibrated empirical law or fitted statistical relation; that classification determines what evidence is needed before transfer to a new catchment.

Units and source audit. For muskingum routing and numerical stability, list the observed quantities, estimated coefficients and chosen reference conditions separately. Include the dates and methods used to obtain each. Calculate one reference state using all quantities in compatible units and round only after a dimensional check. Then reverse a sign convention or convert to a different set of units and verify that the physical meaning survives. A unit-consistent equation may still be wrong because its assumptions do not match the setting.

Applied decision: Test the effect of time step on a routed synthetic hydrograph

The central decision is not whether a formula can be evaluated; it is whether the value carries information sufficient for an engineering action. In this situation, first identify the outlet, observation period and decision variable. Establish a low, central and high forcing scenario using values that can be defended by measurements or explicitly labeled hypotheses. Predict how the response changes before running the model. A prediction that contradicts conservation of water, plausible timing or an independently measured constraint requires investigation even when the graph looks smooth.

For this particular case the dominant caution is that poor K and X estimation can hide flood-wave translation and attenuation errors. Assess which boundary condition is most likely to be violated. For an event analysis, compare the timing of inputs and outputs; for a seasonal balance, separate fluxes from storage; for statistical estimates, distinguish frequency assumptions from observed counts. Do not claim a coefficient calibrated in one environment is universal. State one missing measurement and the practical difficulty of obtaining it.

Mechanism checks and falsification

Counterfactual A — scale. Suppose the area or response time is changed while keeping some input depth or intensity fixed. Decide whether the output should scale linearly, sublinearly or nonlinearly, and explain which physical assumption controls that conclusion. For muskingum routing and numerical stability, the relevant limiting cases include a nearly dry condition, a very large forcing and a coefficient approaching the edge of its documented range. A model that produces negative water depth, impossible storage or negative coefficients without a clear warning is not fit for student use.

Counterfactual B — observation. Ask how a stage recorder, pluviometer, soil-moisture probe or properly maintained flow measurement could disprove the chosen interpretation. Compare two alternatives: an error in the input and an error in the process structure. Error in a timestamp may alter the inferred lag without changing storm totals; an unmodeled flow diversion may alter volume without a rainfall error. Use these fingerprints to choose the next data collection action.

Cross-topic connection. This lesson connects to lakes, reservoirs and flood storage because a model transfers its assumptions through the wider analysis. It also affects linear reservoir response, which may rely on the same boundary records but at a different temporal or spatial resolution. A strong report shows these dependencies explicitly, identifies the quantities shared by two calculations and avoids treating correlated errors as independent measurements.

Field, laboratory and spreadsheet workflow

  1. Draw a control-volume diagram for test the effect of time step on a routed synthetic hydrograph; mark every measured and assumed boundary.
  1. Calculate the result from C0=(Delta t-2KX)/(2K(1-X)+Delta t) where its inputs are available. Show the time-base and area conversions with intermediate numbers rather than only the software answer.
  1. Perform a one-at-a-time sensitivity screen for the most uncertain input. State which physical relationship predicts the direction of the response; explain any threshold or plateau.
  1. Design an independent field check. Check coefficient sum and nonnegativity before trusting a routed point. Explain exactly what would count as disagreement and how you would respond before changing model parameters.
  1. Write a limitation paragraph and a separate decision paragraph. The limitation describes what has not been established; the decision identifies the next action that is justified by the evidence available.

Worked investigative thought experiment

Consider two plausible observations from test the effect of time step on a routed synthetic hydrograph. In the first, the measured quantity changes noticeably after a storm, whereas in the second, a neighboring instrument records a muted response. One explanation is spatial variability in forcing; another is storage and travel-time filtering. A third is mismatched sampling intervals. Quantify the first explanation with an appropriately scaled synthetic input, then test the second through a changed process coefficient. Keep the third as a data-quality hypothesis until the metadata are checked. This separation protects against false certainty: a model that can fit both instruments with different parameters has not identified the underlying cause.

In-class questions and solutions guidance

Question 1. In the context of test the effect of time step on a routed synthetic hydrograph, identify all inputs and outputs of the selected hydrologic control volume. Answer guidance: give measured depths, rates, storage changes and known management transfers; justify the time interval and boundary instead of simply listing variables.

Question 2. Under what conditions is C0=(Delta t-2KX)/(2K(1-X)+Delta t) defensible? Answer guidance: specify the governing conservation or calibration assumptions, a physical range of inputs and at least one failure condition associated with poor k and x estimation can hide flood-wave translation and attenuation errors.

Question 3. How could you distinguish genuine process sensitivity from measurement bias? Answer guidance: formulate competing mechanisms, specify an independent sensor or additional event, and compare timing, volume and directional response.

Question 4. What changes if the same equation is transferred from a small engineered catchment to a large heterogeneous river basin? Answer guidance: evaluate spatial averaging, storage, travel time, changing active area, scale dependence and available calibration evidence.

Extended engineering reflection

A credible result for muskingum routing and numerical stability needs more than an apparently precise scalar. Provide a unit table, parameter source table, version of the governing equation, uncertainty register and plot of response against at least one input. Compare the output to a basic conservation or range check; no negative runoff depth can be created from positive water inputs without explaining a signed flux convention. Consider whether the decision is reversible, such as adjusting a monitoring schedule, or hard to reverse, such as constructing a hydraulic structure. For high-consequence decisions, missing data should lead to targeted measurement and professional review rather than a higher number of decimal places.

§5.5 Lakes, reservoirs and flood storage

Lesson 5.4 — Lakes, reservoirs and flood storage

Hydrological mechanism

Reservoir storage changes according to inflow minus outflow plus other water inputs and losses. A finite-time balance is a useful screening calculation but actual routing needs a stage–storage curve and a stage–discharge relation. Reservoir rules, gates, spillways, safety thresholds and upstream backwater make operational outcomes nonlinear. A single average inlet flow should not be mistaken for a complete flood wave.

Governing expression and dimensions

Teaching relation: S_next=S_now+(I-O) Delta t for consistent units. Every term must be traced to its physical definition before numerical evaluation. A hydrologic depth expressed in millimeters corresponds to one thousand cubic meters for every square kilometer of area. A discharge expressed in cubic meters per second becomes volume only after integrating over seconds. A rate reported over a nominal day or hour is not a spatial average by itself. State whether the equation is a conservation balance, idealization, calibrated empirical law or fitted statistical relation; that classification determines what evidence is needed before transfer to a new catchment.

Units and source audit. For lakes, reservoirs and flood storage, list the observed quantities, estimated coefficients and chosen reference conditions separately. Include the dates and methods used to obtain each. Calculate one reference state using all quantities in compatible units and round only after a dimensional check. Then reverse a sign convention or convert to a different set of units and verify that the physical meaning survives. A unit-consistent equation may still be wrong because its assumptions do not match the setting.

Applied decision: Assess flood-storage reserve before a forecast storm

The central decision is not whether a formula can be evaluated; it is whether the value carries information sufficient for an engineering action. In this situation, first identify the outlet, observation period and decision variable. Establish a low, central and high forcing scenario using values that can be defended by measurements or explicitly labeled hypotheses. Predict how the response changes before running the model. A prediction that contradicts conservation of water, plausible timing or an independently measured constraint requires investigation even when the graph looks smooth.

For this particular case the dominant caution is that outlet release constraints and stage-storage nonlinearity limit linear forecasts. Assess which boundary condition is most likely to be violated. For an event analysis, compare the timing of inputs and outputs; for a seasonal balance, separate fluxes from storage; for statistical estimates, distinguish frequency assumptions from observed counts. Do not claim a coefficient calibrated in one environment is universal. State one missing measurement and the practical difficulty of obtaining it.

Mechanism checks and falsification

Counterfactual A — scale. Suppose the area or response time is changed while keeping some input depth or intensity fixed. Decide whether the output should scale linearly, sublinearly or nonlinearly, and explain which physical assumption controls that conclusion. For lakes, reservoirs and flood storage, the relevant limiting cases include a nearly dry condition, a very large forcing and a coefficient approaching the edge of its documented range. A model that produces negative water depth, impossible storage or negative coefficients without a clear warning is not fit for student use.

Counterfactual B — observation. Ask how a stage recorder, pluviometer, soil-moisture probe or properly maintained flow measurement could disprove the chosen interpretation. Compare two alternatives: an error in the input and an error in the process structure. Error in a timestamp may alter the inferred lag without changing storm totals; an unmodeled flow diversion may alter volume without a rainfall error. Use these fingerprints to choose the next data collection action.

Cross-topic connection. This lesson connects to event hydrographs and baseflow separation because a model transfers its assumptions through the wider analysis. It also affects muskingum routing and numerical stability, which may rely on the same boundary records but at a different temporal or spatial resolution. A strong report shows these dependencies explicitly, identifies the quantities shared by two calculations and avoids treating correlated errors as independent measurements.

Field, laboratory and spreadsheet workflow

  1. Draw a control-volume diagram for assess flood-storage reserve before a forecast storm; mark every measured and assumed boundary.
  1. Calculate the result from S_next=S_now+(I-O) Delta t for consistent units where its inputs are available. Show the time-base and area conversions with intermediate numbers rather than only the software answer.
  1. Perform a one-at-a-time sensitivity screen for the most uncertain input. State which physical relationship predicts the direction of the response; explain any threshold or plateau.
  1. Design an independent field check. Close a water balance for three sequential time steps. Explain exactly what would count as disagreement and how you would respond before changing model parameters.
  1. Write a limitation paragraph and a separate decision paragraph. The limitation describes what has not been established; the decision identifies the next action that is justified by the evidence available.

Worked investigative thought experiment

Consider two plausible observations from assess flood-storage reserve before a forecast storm. In the first, the measured quantity changes noticeably after a storm, whereas in the second, a neighboring instrument records a muted response. One explanation is spatial variability in forcing; another is storage and travel-time filtering. A third is mismatched sampling intervals. Quantify the first explanation with an appropriately scaled synthetic input, then test the second through a changed process coefficient. Keep the third as a data-quality hypothesis until the metadata are checked. This separation protects against false certainty: a model that can fit both instruments with different parameters has not identified the underlying cause.

In-class questions and solutions guidance

Question 1. In the context of assess flood-storage reserve before a forecast storm, identify all inputs and outputs of the selected hydrologic control volume. Answer guidance: give measured depths, rates, storage changes and known management transfers; justify the time interval and boundary instead of simply listing variables.

Question 2. Under what conditions is S_next=S_now+(I-O) Delta t for consistent units defensible? Answer guidance: specify the governing conservation or calibration assumptions, a physical range of inputs and at least one failure condition associated with outlet release constraints and stage-storage nonlinearity limit linear forecasts.

Question 3. How could you distinguish genuine process sensitivity from measurement bias? Answer guidance: formulate competing mechanisms, specify an independent sensor or additional event, and compare timing, volume and directional response.

Question 4. What changes if the same equation is transferred from a small engineered catchment to a large heterogeneous river basin? Answer guidance: evaluate spatial averaging, storage, travel time, changing active area, scale dependence and available calibration evidence.

Extended engineering reflection

A credible result for lakes, reservoirs and flood storage needs more than an apparently precise scalar. Provide a unit table, parameter source table, version of the governing equation, uncertainty register and plot of response against at least one input. Compare the output to a basic conservation or range check; no negative runoff depth can be created from positive water inputs without explaining a signed flux convention. Consider whether the decision is reversible, such as adjusting a monitoring schedule, or hard to reverse, such as constructing a hydraulic structure. For high-consequence decisions, missing data should lead to targeted measurement and professional review rather than a higher number of decimal places.

§5.6 Applied project 5 — Wetland flood peak attenuation

Applied project 5 — Wetland flood peak attenuation

A freshwater wetland receiving a constant 50 m³/s from an upstream channel initially releases 10 m³/s. A linear storage timescale of 5 hours predicts an increasing outflow but no immediate match with inflow. Downstream bridge operations depend on the peak time rather than just the total event volume.

Project deliverable. Calculate outflow at 1, 3 and 6 hours and explain why an actual wetland stage-discharge function may not be linear. Include an explicit data dictionary, local measurement geometry, at least two competing model structures, a water-balance closure or proxy and a sensitivity analysis. Use the two paired laboratories for this unit as exploratory tools only. A peer should be able to reproduce the results using the included synthetic inputs and code without requesting undocumented assumptions.

Unit review and source-of-error audit

Before claiming success, ask whether the result is constrained by water conservation, the measurement period, the sampling interval and the domain of the empirical relation. Identify how a missing observation could change a decision, not merely how it changes the last decimal digit. The strongest student reports explain why a model is unsuitable for a particular site, include a second defensible interpretation and propose a feasible experiment to choose between them.

THREE WORKED QUESTIONS

Step-by-Step Hydrology Solutions

Three original questions connect calculations, definitions and interpretation. Open each solution to follow the reasoning.

Solved Problem Example 5.1: Worked Problem 5.1: Energy and Mass Conservation Balance

A civil hydraulic/water resources system is analyzed under Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing conditions. Evaluate head loss, pressure differentials, and flow discharge along the hydraulic path. Formulate the governing equations, carry out dimensional verification, and verify limiting behavior.

Step 1: Problem Geometry and Input Parameters Identify cross-sectional dimensions, fluid density $\rho$, gravity $g$, boundary elevations, and roughness coefficients for Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing.

Step 2: Governing Flow and Energy Equations Apply the governing equations of open-channel hydraulics or pipe fluid mechanics:

$$Q = \frac{1}{n} A R^{2/3} S^{1/2}, \quad h_f = f \frac{L}{D} \frac{v^2}{2g}$$

Step 3: Numerical Resolution in Base SI Units Compute intermediate values and resolve the target hydraulic parameters with rigorous dimensional tracking:

$$\text{Computed Hydraulic Parameter for Unit } 5 = \text{Consistently Verified against Governing Equations}$$

Step 4: Boundary Consistency & Engineering Validation Confirm that velocities remain within non-scouring and non-silting limits, Reynolds number verifies the flow regime, and energy conservation is fully satisfied.

Solved Problem Example 5.2: Worked Problem 5.2: Hydraulic Gradient & Water Surface Profiles

A civil hydraulic/water resources system is analyzed under Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing conditions. Compute specific energy, critical depth, and classify water surface profile curvature. Formulate the governing equations, carry out dimensional verification, and verify limiting behavior.

Step 1: Problem Geometry and Input Parameters Identify cross-sectional dimensions, fluid density $\rho$, gravity $g$, boundary elevations, and roughness coefficients for Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing.

Step 2: Governing Flow and Energy Equations Apply the governing equations of open-channel hydraulics or pipe fluid mechanics:

$$Q = \frac{1}{n} A R^{2/3} S^{1/2}, \quad h_f = f \frac{L}{D} \frac{v^2}{2g}$$

Step 3: Numerical Resolution in Base SI Units Compute intermediate values and resolve the target hydraulic parameters with rigorous dimensional tracking:

$$\text{Computed Hydraulic Parameter for Unit } 5 = \text{Consistently Verified against Governing Equations}$$

Step 4: Boundary Consistency & Engineering Validation Confirm that velocities remain within non-scouring and non-silting limits, Reynolds number verifies the flow regime, and energy conservation is fully satisfied.

Solved Problem Example 5.3: Worked Problem 5.3: Hydrodynamic Forces & Momentum Conservation

A civil hydraulic/water resources system is analyzed under Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing conditions. Determine thrust forces, hydrodynamic momentum flux, and anchor block resistance. Formulate the governing equations, carry out dimensional verification, and verify limiting behavior.

Step 1: Problem Geometry and Input Parameters Identify cross-sectional dimensions, fluid density $\rho$, gravity $g$, boundary elevations, and roughness coefficients for Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing.

Step 2: Governing Flow and Energy Equations Apply the governing equations of open-channel hydraulics or pipe fluid mechanics:

$$Q = \frac{1}{n} A R^{2/3} S^{1/2}, \quad h_f = f \frac{L}{D} \frac{v^2}{2g}$$

Step 3: Numerical Resolution in Base SI Units Compute intermediate values and resolve the target hydraulic parameters with rigorous dimensional tracking:

$$\text{Computed Hydraulic Parameter for Unit } 5 = \text{Consistently Verified against Governing Equations}$$

Step 4: Boundary Consistency & Engineering Validation Confirm that velocities remain within non-scouring and non-silting limits, Reynolds number verifies the flow regime, and energy conservation is fully satisfied.

Solved Problem Example 5.4: Worked Problem 5.4: Hydraulic Efficiency, Loss Coefficients & Safety Margins

A civil hydraulic/water resources system is analyzed under Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing conditions. Assess minor losses, cavitation risk, and design safety margins against overflow. Formulate the governing equations, carry out dimensional verification, and verify limiting behavior.

Step 1: Problem Geometry and Input Parameters Identify cross-sectional dimensions, fluid density $\rho$, gravity $g$, boundary elevations, and roughness coefficients for Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing.

Step 2: Governing Flow and Energy Equations Apply the governing equations of open-channel hydraulics or pipe fluid mechanics:

$$Q = \frac{1}{n} A R^{2/3} S^{1/2}, \quad h_f = f \frac{L}{D} \frac{v^2}{2g}$$

Step 3: Numerical Resolution in Base SI Units Compute intermediate values and resolve the target hydraulic parameters with rigorous dimensional tracking:

$$\text{Computed Hydraulic Parameter for Unit } 5 = \text{Consistently Verified against Governing Equations}$$

Step 4: Boundary Consistency & Engineering Validation Confirm that velocities remain within non-scouring and non-silting limits, Reynolds number verifies the flow regime, and energy conservation is fully satisfied.

Solved Problem Example 5.5: Worked Problem 5.5: Extreme Flood Events, Dynamic Routing & System Resilience

A civil hydraulic/water resources system is analyzed under Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing conditions. Evaluate routing performance under peak hydrograph conditions and emergency capacities. Formulate the governing equations, carry out dimensional verification, and verify limiting behavior.

Step 1: Problem Geometry and Input Parameters Identify cross-sectional dimensions, fluid density $\rho$, gravity $g$, boundary elevations, and roughness coefficients for Unit Hydrographs, Storage Reservoirs and Muskingum Flood Routing.

Step 2: Governing Flow and Energy Equations Apply the governing equations of open-channel hydraulics or pipe fluid mechanics:

$$Q = \frac{1}{n} A R^{2/3} S^{1/2}, \quad h_f = f \frac{L}{D} \frac{v^2}{2g}$$

Step 3: Numerical Resolution in Base SI Units Compute intermediate values and resolve the target hydraulic parameters with rigorous dimensional tracking:

$$\text{Computed Hydraulic Parameter for Unit } 5 = \text{Consistently Verified against Governing Equations}$$

Step 4: Boundary Consistency & Engineering Validation Confirm that velocities remain within non-scouring and non-silting limits, Reynolds number verifies the flow regime, and energy conservation is fully satisfied.