ENGINEERING HYDROLOGY (ENCE 306) Past Year Question Solution

Humankind interferes with the hydrological cycle in many ways by altering natural water movement and storage:

• Urbanization: Replaces natural land surfaces with impervious materials, reducing infiltration and increasing surface runoff, which often leads to flooding.
• Deforestation & Afforestation: Deforestation decreases interception and transpiration, causing more runoff and soil erosion, while afforestation and irrigation increase evapotranspiration.
• Dams and Reservoirs: Construction of dams modifies natural stream flow patterns, increases evaporation losses, and alters downstream water availability.
• Groundwater Extraction: Excessive groundwater extraction lowers the water table and reduces base flow to rivers, whereas artificial recharge structures can enhance groundwater replenishment.
• Other Activities: Industrial activities, agriculture, and pollution also affect precipitation patterns and degrade water quality.

Overall, human activities disturb the natural balance of the hydrological cycle, resulting in problems such as floods, droughts, and groundwater depletion.

Definition: Hydrology is the scientific study of the movement, distribution, and quality of water on Earth, including the hydrological cycle and water resources. It is the science that deals with the occurrence, circulation, and distribution of water of the earth and earth’s atmosphere. Engineering hydrology bridges the gap between science and civil engineering.

Scope and Applications:

• To measure hydrological data such as rainfall, stream flow, evaporation, and groundwater levels.
• To calculate rainfall, runoff, floods, and drought characteristics.
• To find out the maximum probable flood at dams, reservoirs, bridges, and city drainage systems.
• To forecast floods and provide data for flood routing and early warning systems.
• It is used in the design and operation of hydraulic structures.
• It is used for hydropower generation and planning.
• It determines the water balance for a particular region and finds the relation between catchment surface water and groundwater resources.
• It helps to assess the impact of climate change on rainfall, runoff, and water availability.
• It helps in the control of erosion to minimize the sedimentation of reservoirs.

The study of hydrology provides the scientific foundation needed to understand, manage, and utilize water effectively for hydraulic structures. Its importance includes:

• Estimating Design Floods: Helps to estimate the design flood (peak discharge) required for deciding the capacity, span, and waterway of structures like dams, spillways, and bridges.
• Determining Highest Flood Level (HFL): Necessary for fixing bridge deck levels, freeboards for dams, and safe elevations for other structures.
• Calculating Scour Depth: Helps control the depth of the foundation and dictates the type of foundation used to prevent structural failure.
• Flow Characteristics: Provides understanding of river flow characteristics such as velocity, discharge variation, and flood frequency to ensure safe hydraulic design.
• Durability and Economy: Ensures long-term durability and economy by reducing the risks of flood-induced damage, failure, excessive afflux, or waterlogging.

Importance for Water Resources Projects: The study of hydrology is critical for engineers because it helps understand the behavior of water on the Earth’s surface and below it. This is necessary for planning, designing, and managing water resource projects (especially in topographically diverse regions like Nepal). It helps assess surface water availability, predict hydrologic extremes (floods and droughts) for risk assessment, design hydropower projects, and plan drainage/flood control systems by analyzing rainfall patterns.

Significant Features of Water Balance: The water balance equation is based on the law of conservation of mass. The change in storage per unit time on a control volume equals the sum of inflows minus the sum of outflows. The general water budget equation can be expressed as:

Where: $P$ = Precipitation (Inflow), $R$ = Surface runoff (Outflow), $G$ = Net ground water flow out of the catchment (Outflow), $E$ = Evaporation (Outflow), $T$ = Transpiration (Outflow), $\Delta S$ = Change in storage.

The Hydrologic Cycle: The hydrologic cycle is the continuous, natural process by which water circulates from the Earth’s surface to the atmosphere and back again. It is driven by solar energy and gravity.

(Draw a schematic sketch showing: Ocean → Evaporation → Clouds → Precipitation over land/ocean → Surface Runoff/Infiltration/Groundwater Flow → Return to Ocean.)

Components of the Hydrologic Cycle:

1. Transportational Components:

• Precipitation: Water falling from clouds to the surface (rain, snow).
• Evaporation: Water turning into vapor from oceans, land, and water bodies.
• Transpiration: Release of water vapor into the atmosphere by vegetation.
• Infiltration: Water penetrating into the soil layer from the surface.
• Runoff: Water flowing over the land surface entering stream channels.

2. Storage Components:

• Surface Storage: Depression storage, ponds, lakes, and oceans.
• Soil Moisture Storage: Water held in the unsaturated soil zone.
• Groundwater Storage: Water held in deep aquifers and rock formations.

Water Balance/Budget Equation: It is an equation based on the continuity principle (conservation of mass) applied to a hydrologic system (like a catchment). It equates the difference between water inflow and outflow to the change in water storage over a specific time interval.

Components:

• Inflow: Primarily Precipitation ($P$) — rainfall, snowfall.
• Outflow: Includes Runoff ($R$), Evaporation ($E$), Transpiration ($T$), and Net Groundwater Outflow ($G$).
• Storage ($\Delta S$): Changes in surface water bodies, soil moisture, and groundwater reserves.

Significance in Engineering Hydrology:

• Resource Assessment: Verifies how much water is available, where it comes from, and how it changes over time.
• Estimating Unknowns: If all but one component is known, the equation allows engineers to calculate the missing parameter (e.g., estimating evaporation losses).
• Hydrologic Modeling: Supports accurate modeling and watershed management for water supply, irrigation, and drought prediction.

Causes of Inconsistencies in Rainfall Data:

• Shifting of a Rain Gauge: Moving the station to a new location with different topographical conditions.
• Change in Neighborhood Environment: Growth of tall trees, construction of tall buildings, or other obstructions near the gauge over time.
• Change in Ecosystem/Exposure: A severe forest fire, landslide, or deforestation changing the local microclimate.
• Change in Observational Procedure: Errors or changes in the method of taking readings, or replacement of the observer.
• Change in Instrument: Upgrading or replacing a non-recording gauge with a different type of recording gauge.

Double Mass Curve Analysis:

• Principle: A plot of accumulated rainfall of the station in question against the concurrent accumulated average rainfall of a group of surrounding base stations should yield a straight line if the data is consistent.
• Procedure:
  1. Select a group of 5 to 10 reliable neighboring base stations near station “X”.
  2. Arrange the annual rainfall data in reverse chronological order (latest to oldest).
  3. Calculate cumulative rainfall of station X ($\Sigma P_x$) and cumulative mean annual rainfall of base stations ($\Sigma P_{avg}$).
  4. Plot $\Sigma P_{avg}$ on the X-axis and $\Sigma P_x$ on the Y-axis.
• Analysis: A straight line = consistent data. A break in slope at a particular year = inconsistency from that year onwards.
• Correction Factor:
  $$P_{cx} = P_x \times \frac{M_c}{M_a}$$
  Where $M_c$ = Corrected slope, $M_a$ = Original slope.

Different Types of Recording Rain Gauges:

• Tipping Bucket Type: Uses a divided bucket mechanism that tips when a specific volume of water is collected, sending an electrical pulse.
• Weighing Bucket Type: Catches rain in a bucket resting on a spring or weighing scale. The increasing weight translates to rainfall depth recorded on a clock-driven chart. Also suitable for recording snow.
• Natural Syphon (Float) Type: Rain falls into a chamber containing a float. As water level rises, the float rises, and a pen records the depth on a rotating drum. A syphon mechanism empties it when full.

Working Principle of Tipping Bucket Rain Gauge:

• Components: A circular receiver funnel (usually 30 cm diameter) directs rainwater into a pair of small bucket compartments balanced on a horizontal pivot.
• Working Mechanism:
  1. Rain falls through the funnel into one of the two bucket compartments.
  2. When a specific amount (usually 0.25 mm) collects, the bucket becomes top-heavy and tips downward.
  3. This simultaneously empties the water into a reservoir and positions the second, empty bucket under the funnel.
  4. Each tip closes an electrical switch, sending a pulse to an electronic data logger or a pen on a clock-driven drum.
  5. The number of tips per unit time gives rainfall intensity.

Neat Sketch Description: (Draw a funnel at the top. Below the funnel, draw a V-shaped object balanced on a central pivot. Label “Bucket 1” and “Bucket 2”. Connect the pivot to a “Data Logger/Recorder”. Draw a lower “Measuring Tube/Reservoir” beneath.)

The optimum number of rain gauges is determined statistically based on the allowable percentage of error in estimating mean rainfall.

Step 1: Calculate the Mean Rainfall ($\bar{P}$):

Step 2: Calculate the Standard Deviation ($\sigma_{n-1}$):

Step 3: Calculate the Coefficient of Variation ($C_v$):

Step 4: Determine Optimum Number ($N$): If $\epsilon$ is the allowable percentage error (usually 10%):

If $N > n$, additional gauges $(N – n)$ are required. If $N \leq n$, the existing network is adequate.

Missing Precipitation Data: Missing precipitation data refers to gaps in the continuous historical rainfall records at a specific gauge station due to instrument failure, vandalism, absence of the observer, or natural disasters.

1. Simple Arithmetic Mean Method (Index Method):

• Condition: Used if the Normal Annual Precipitation of surrounding index stations is within 10% of the Normal Annual Precipitation of the station with missing data ($N_x$).
• Formula:
  $$P_x = \frac{P_1 + P_2 + … + P_m}{m}$$

2. Normal Ratio Method:

• Condition: Used if the Normal Annual Precipitation of any surrounding station differs from $N_x$ by more than 10%.
• Formula:
  $$P_x = \frac{N_x}{m} \left( \frac{P_1}{N_1} + \frac{P_2}{N_2} + … + \frac{P_m}{N_m} \right)$$
  Where $N_1, N_2… N_m$ are the normal annual precipitations of the $m$ surrounding stations.

1. Arithmetic Mean Method:

• The simplest method. Average rainfall is the sum of depths at all stations divided by the number of stations.
• Formula: $\bar{P} = \frac{\Sigma P_i}{n}$
• Suitable only if the catchment is very flat, gauges are uniformly distributed, and rainfall varies little across the area.

2. Thiessen Polygon Method:

• Assigns a “weight” (area of influence) to each rain gauge by drawing perpendicular bisectors to the lines connecting adjacent gauges.
• Formula: $\bar{P} = \frac{\Sigma P_i A_i}{\Sigma A_i}$
• Better than arithmetic method but does not account for topographic effects.

3. Isohyetal Method:

• The most accurate method. Isohyets (lines of equal rainfall depth) are drawn across the catchment map using interpolation.
• Formula: $\bar{P} = \frac{\Sigma \left( \frac{P_i + P_{i+1}}{2} \right) \times A_i}{\Sigma A_i}$
• Highly accurate because it allows the analyst to consider topographic effects when drawing the isohyets.

Different Forms of Precipitation:

• Rain: Drops of liquid water larger than 0.5 mm in diameter.
• Drizzle: Fine droplets smaller than 0.5 mm in diameter, usually falling from stratus clouds.
• Snow: Ice crystals formed by sublimation of water vapor in the atmosphere.
• Sleet: Frozen raindrops or transparent ice pellets formed when rain falls through a freezing layer of air.
• Hail: Lumps of ice (usually >8 mm diameter) formed in strong thunderstorm updrafts.
• Glaze: Freezing rain that forms a coating of ice when liquid rain strikes a cold surface (below freezing).

Mechanism of Precipitation Formation: Three primary conditions must be met: (1) Cooling of a moist air mass to its dew point, (2) Presence of condensation nuclei (dust, smoke, salt), and (3) Droplet growth through coalescence or the Bergeron ice-crystal process.

Types Based on Lifting Mechanism:

• Convective Precipitation: Caused by natural upward heating of lighter, warmer air near the earth’s surface. Typical of summer thunderstorms — intense, short duration, highly localized.
• Orographic Precipitation: Caused when a moist air mass is forced to rise over a natural topographic barrier (mountain range). The windward side receives heavy rain; the leeward side is a dry “rain shadow”.
• Cyclonic (Frontal) Precipitation: Occurs when air masses with different temperatures clash. Cold front causes heavy, short rain; Warm front causes steady, prolonged rain.

(i) Depth-Area-Duration (DAD) Curves:

• Concept: DAD curves express the relationship between the maximum average depth of rainfall over a specific area for a specific duration of a storm. As the area increases, the average depth decreases. For a given area, depth increases as duration increases.
• Equation: $P_A = P_0 \exp(-K A^n)$, where $P_A$ is the average depth over area $A$, $P_0$ is the maximum point rainfall at the center, and $K, n$ are constants.
• Application: Crucial for determining the maximum design storm over a specific catchment area, required for designing large dams and spillways.

(ii) Intensity-Duration-Frequency (IDF) Relationship:

• Concept: Links the intensity of rainfall (mm/hr), the duration (minutes or hours), and the frequency (Return Period, $T$, in years). Higher intensity storms have shorter durations (for a given return period).
• Equation: $i = \frac{K T^x}{(D + a)^n}$, where $i$ is intensity, $T$ is return period, $D$ is duration, and $K, x, a, n$ are regional constants.
• Application: Fundamental input for the Rational Method, extensively used in the design of urban drainage systems, culverts, and storm sewers.

Probable Maximum Precipitation (PMP): PMP is defined as the greatest depth of precipitation for a given duration that is meteorologically possible over a given size storm area at a particular geographical location. It represents the absolute upper limit of rainfall, used to design the Probable Maximum Flood (PMF) for high-risk structures like spillways of large dams.

Estimation of PMP:

• Meteorological Method (Moisture Maximization): The most severe historical storms recorded in or near the region are analyzed and scaled up (maximized) by a maximization ratio, assuming the atmosphere held the maximum possible precipitable water (based on maximum dew point temperatures).
• Statistical Method (Hershfield Technique): Uses the mean ($\bar{P}$) and standard deviation ($\sigma_n$) of annual maximum rainfall series:
  $$PMP = \bar{P} + K \sigma_n$$
  Where $K$ is an empirical frequency factor (often around 15, based on Hershfield’s envelope curve of worldwide extreme data).

When precipitation falls over a catchment, not all of it contributes to surface runoff. The portion “lost” before reaching the stream includes:

• Interception: Rainfall caught and retained by the foliage of plants, buildings, and other objects before reaching the ground. It eventually evaporates back into the atmosphere.
• Depression Storage: Water trapped in small depressions, puddles, and hollows on the ground surface. This water either evaporates or infiltrates.
• Evaporation: Physical process by which liquid water from free water surfaces (lakes, rivers, reservoirs) is converted into water vapor.
• Transpiration: Process by which water is absorbed by plant roots, transported through the plant, and released as vapor through the stomata of the leaves.
• Infiltration: Water entering the soil surface. The most significant loss component determining the surface runoff.

Evaporation is the physical process by which a liquid (water) changes into a gaseous state (water vapor) at the free surface below its boiling point through the transfer of heat energy.

Physical Mechanism: Water molecules are in constant motion. Evaporation occurs when the kinetic energy of water molecules at the surface is sufficient to overcome intermolecular cohesive forces. According to Dalton’s Law, the rate of evaporation ($E$) is proportional to the vapor pressure deficit:

Where $e_w$ = saturation vapor pressure at the water temperature, $e_a$ = actual vapor pressure in the overlying air, $C$ = a constant.

Meteorological Factors affecting Evaporation:

• Solar Radiation / Temperature: Higher solar radiation increases water temperature and the saturation vapor pressure ($e_w$), exponentially increasing the evaporation rate.
• Vapor Pressure: Evaporation is directly proportional to the vapor pressure deficit ($e_w – e_a$). Dry air promotes high evaporation.
• Wind Speed: Wind removes saturated air above the water surface and replaces it with drier air, maintaining a high vapor pressure gradient.
• Atmospheric Pressure: A decrease in atmospheric pressure slightly increases the evaporation rate.
• Water Quality: Soluble salts reduce the saturated vapor pressure; saline water evaporates about 2-3% less than fresh water.

The energy balance method applies the Law of Conservation of Energy to a water body. The energy balance equation for a lake is:

Where $R_n$ = Net radiation, $H$ = Sensible heat to atmosphere, $\lambda E$ = Latent heat of evaporation, $G$ = Heat conducted to ground, $\Delta S$ = Change in heat storage.

For practical estimation over longer periods, $G$ and $\Delta S$ are often negligible. The sensible heat ($H$) is related to latent heat using the Bowen ratio ($B$), where $B = \frac{H}{\lambda E}$:

Solving for the rate of evaporation:

Where $E$ is in mm/day, $R_n$ in cal/cm²/day, and $\lambda \approx 585$ cal/g.

Evapotranspiration (ET): Also known as Consumptive Use, it is the combined process of water loss from a given land area to the atmosphere. It includes evaporation from the soil surface and transpiration from plants. When sufficient moisture is continuously available to fully meet vegetation needs, the resulting ET is called Potential Evapotranspiration (PET).

Penman’s Equation: A physically-based method that combines the energy balance approach and the aerodynamic (mass transfer) approach to calculate PET:

Explanation of Parameters:

• $PET$ = Potential Evapotranspiration (mm/day).
• $A$ = Slope of the saturation vapor pressure vs. temperature curve at the mean air temperature (mm of Hg / °C).
• $H_n$ = Net incoming solar radiation expressed in mm of evaporable water per day.
• $E_a$ = Aerodynamic parameter incorporating wind speed and vapor pressure deficit: $E_a = f(u) \cdot (e_w – e_a)$.
• $\gamma$ = Psychrometric constant (approx. 0.49 mm of Hg / °C).

Initial Loss: The portion of precipitation that occurs before any surface runoff begins — including water required to fulfill interception by vegetation and to fill small surface depressions.

Infiltration: The physical process of water penetrating the soil surface and moving downwards. It directly dictates how much rainfall converts to surface runoff versus groundwater recharge.

Horton’s Infiltration Curve: Horton established an empirical relationship describing how infiltration rate decreases exponentially over time during continuous rainfall:

Parameters:

• $f_t$: Infiltration capacity at any time $t$ (cm/hr or mm/hr).
• $f_o$: Initial Infiltration Capacity at $t = 0$. High because the soil is dry and capillary forces strongly draw water in.
• $f_c$: Ultimate/Equilibrium Infiltration Capacity. The rate at which the capacity levels off once soil pores are saturated. Roughly equal to the saturated hydraulic conductivity of the soil.
• $k$: Decay Constant ($time^{-1}$). Dictates how rapidly the rate drops from $f_o$ to $f_c$. Depends on soil characteristics and vegetation.
• $e$: Base of the natural logarithm.

(A sketch showing ‘t’ on the X-axis and ‘f’ on the Y-axis, starting high at $f_o$ and curving down exponentially to a horizontal asymptote at $f_c$, is expected for full marks.)

Phi-index ($\Phi$-index): The average rate of rainfall intensity above which the volume of rainfall exactly equals the volume of direct surface runoff. It is a constant horizontal line on a hyetograph. All rainfall above this line is considered effective rainfall (i.e., it becomes runoff).

W-index: The average infiltration rate during the time when rainfall intensity exceeds the infiltration capacity. It explicitly separates and removes initial losses (interception and depression storage) from total abstraction:

Where $P$ = Total precipitation, $R$ = Total Runoff, $I_{loss}$ = Initial losses, $t_e$ = duration of rainfall excess.

Which is practically more useful? The $\Phi$-index is more practically useful because:

1. Simplicity: It is very difficult to accurately measure initial losses ($I_{loss}$). The $\Phi$-index bypasses this problem by lumping all losses into one uniform rate.
2. Adequacy for Major Storms: In large flood-producing storms, initial losses are small compared to total infiltration and runoff, so lumping them does not cause significant errors.
3. Ease of Application: It allows engineers to easily separate rainfall excess from a hyetograph to derive Unit Hydrographs and estimate flood peaks.

Runoff: Surface runoff is the portion of precipitation that flows over the land surface toward streams, rivers, or drains. It occurs when rainfall intensity exceeds the infiltration capacity of the soil, or when the soil becomes fully saturated.

Factors Affecting Runoff:

A. Climatic Factors:

• Types of precipitation: Rain produces higher runoff compared to snow.
• Rainfall intensity: Higher intensity results in greater surface runoff.
• Duration of rainfall: As duration increases, soil infiltration capacity decreases, causing more runoff.
• Temperature: Higher temperatures increase evaporation, resulting in less runoff.
• Humidity: Lower humidity leads to more evaporation, reducing runoff.

B. Physiographic Factors:

• Size of basin: Larger watersheds take longer to deliver runoff to the outlet.
• Shape of watershed: Elongated catchments delay peak runoff; circular catchments produce quicker peak flows.
• Soil type: Sandy soils have high infiltration (less runoff); clayey soils have low infiltration (more runoff).
• Land use: Forested areas offer high resistance to flow, leading to less runoff.
• Drainage density: A higher drainage density results in faster runoff and shorter lag time.

1. Linear Regression Method (Rainfall-Runoff Relationship): Fits a linear regression line between Runoff ($R$) and Precipitation ($P$): $R = aP + b$. The reliability is checked using the correlation coefficient ($r$). For larger areas, an exponential relation $R = \alpha P^\beta$ is used.

2. Regional Formulae for Monthly Flows:

• Medium Irrigation Project (MIP) Method: For catchments less than $100 \text{ km}^2$. Requires one low-flow measurement (Nov-Apr) and uses pre-established monthly coefficients to estimate flows for the entire year.
• WECS/DHM Method: Developed by Nepal’s Water and Energy Commission Secretariat. Suitable for catchments $> 100 \text{ km}^2$. Uses catchment area (below 5000 m) and mean monsoon precipitation to calculate monthly flows.
• MHSP Method: Developed under Nepal’s Medium Hydropower Study Project for catchments $> 50\text{-}100 \text{ km}^2$. Divides the catchment by altitude zones to estimate flows for hydropower and irrigation planning.

The Medium Irrigation Project (MIP) method is an empirical technique for estimating the distribution of monthly flows for ungauged locations in Nepal.

Key Characteristics and Methodology:

• Catchment Suitability: Gives better results for smaller catchment areas, specifically those less than $100 \text{ km}^2$.
• Data Requirement: At least one flow measurement during the low-flow period (November to April) is mandatory.
• Use of Coefficients: Pre-established regional coefficients exist for flow measurements taken on the 15th of a particular month.
• Interpolation: If the measurement is taken on any other date, the coefficient must be found using interpolation.

Calculation Steps:

1. Calculate the base “April flow” using the field-measured discharge and the corresponding monthly coefficient:
   $$\text{April flow} = \frac{\text{Measured discharge}}{\text{Coefficient of the measurement month}}$$
2. Estimate the average flow for any other month:
   $$\text{Monthly flow} = \text{April flow} \times \text{Monthly coefficient}$$

Definition: A Flow Duration Curve (FDC) is a cumulative frequency curve that shows the percentage of time a specified discharge is equaled or exceeded during a given period.

Characteristics:

• Data Period Requirement: The period of data used must be an integer multiple of water years to ensure a complete cycle of climatic changes is represented.
• Shape: A steep slope indicates a flashy stream with highly variable flow. A flat slope indicates a steady stream with high groundwater (base flow) contribution.

Practical Applications:

• Hydropower Design: Determines the dependable flow available for electricity generation (e.g., $Q_{90}$ or $Q_{50}$ flows).
• Water Supply Planning: Assesses the reliability of a river for continuous domestic or irrigation water supply.
• Water Quality Management: Determines design flows for waste load allocations and dilution capabilities of a stream.
• Irrigation Projects: Determines the proportion of time adequate water will be available for agricultural commands.

Different Methods of Measuring Streamflow:

1. Direct Methods:

• Area-Velocity Method (using current meters or floats).
• Dilution Techniques (using chemical or radioactive tracers).
• Electromagnetic Method (based on Faraday’s law).
• Ultrasonic Method (using ultrasonic sound waves across the channel).

2. Indirect Methods:

• Hydraulic Structures (weirs, flumes, or notches).
• Slope-Area Method (using Manning’s equation and high-water marks).

The Area-Velocity Method (In Detail): Based on the continuity equation: $Q = A \times V$.

Procedure:

1. Site Selection: A straight, uniform reach with stable banks and bed is selected.
2. Dividing the Cross-Section: The width is divided into 15 to 30 segments so that no single segment carries more than 10% of total discharge.
3. Measuring Width and Depth: Distance from a reference bank and depth at each vertical are measured.
4. Measuring Velocity with a Current Meter:
   • If depth $y < 3m$: Velocity measured at $0.6y$ from the surface ($v_{avg} = v_{0.6}$).
   • If depth $y > 3m$: Velocity measured at $0.2y$ and $0.8y$, and average is taken ($v_{avg} = \frac{v_{0.2} + v_{0.8}}{2}$).
5. Computation: Discharge for each segment is $\Delta Q_i = A_i \times v_i$. Total discharge: $Q = \sum \Delta Q_i$.

The mid-section method is the most common computational technique in the Area-Velocity method.

Principle: The velocity measured at a given vertical represents the mean velocity for a segment of width extending halfway to the adjacent verticals on either side.

Calculation: For the $i^{th}$ vertical at position $x_i$ with depth $y_i$ and mean velocity $v_i$:

Description: The Slope-Area method is an indirect method used post-flood to estimate the peak flood discharge when direct measurement was impossible. It relies on Manning’s equation:

Procedure:

1. Site Selection: A straight reach with uniform cross-section and clearly visible high-water marks (mud lines) is chosen.
2. Field Survey: Cross-sections are surveyed at the upstream and downstream ends.
3. Determining Slope: Water surface elevations at both sections are determined from high-water marks. Energy slope: $S = \frac{Z_1 – Z_2}{L}$.
4. Estimating ‘n’: Manning’s roughness coefficient is estimated based on riverbed material and vegetation.
5. Computation: Average $A$ and $R$ are calculated and plugged into Manning’s equation.

Limitations:

• Accuracy of ‘n’: The most severe limitation. A small error in ‘n’ causes a proportional error in discharge.
• Identification of Marks: High-water marks may be indistinct, washed away, or affected by local wave action.
• Steady-Uniform Flow Assumption: Manning’s equation assumes steady and uniform flow, whereas flood flows are inherently unsteady and non-uniform.
• Low Accuracy: Overall yields rough estimates, often with errors of ±15-25%.

Stage-Discharge Relationship: The relationship between the stage (water surface elevation, $h$) and the discharge ($Q$). Once established, only the daily stage reading is needed to find the corresponding discharge. The mathematical relationship is:

Where $a$ = zero-flow stage, $C_r$ and $\beta$ = rating curve constants.

Rating Curve: The graphical plot of Stage ($h$) on the Y-axis versus Discharge ($Q$) on the X-axis. On arithmetic paper it forms a parabolic curve; on log-log paper it becomes a straight line.

Extrapolation Methods:

1. Logarithmic Method: Because the rating curve is a straight line on log-log paper, the line fitted to the known lower-stage data is simply extended outward to the required high flood stage.
2. Conveyance Method: Conveyance $K$ is defined from Manning’s equation as $K = \frac{1}{n} A R^{2/3}$. Since $K$ is a purely geometric property, it can be calculated at high stages via topographical survey. A curve of Stage vs. $K$ is plotted and extended to high stages. For a given high stage, read $K$ and the ratio $(Q/K)$ to calculate extrapolated $Q$.

Hydrograph: A graphical representation of the discharge (volume of flow per unit time) of a stream at a specific location plotted against time. It reflects the continuous response of a catchment to rainfall or snowfall over time.

Components of a Flood Hydrograph:

• Rising Limb (Concentration Curve): The ascending portion. Represents the increase in discharge as more areas of the catchment contribute to the outlet. Its shape depends on rainfall intensity, duration, and catchment properties.
• Crest Segment (Peak): The highest part containing the peak flow ($Q_p$). Represents the point where the entire catchment is contributing to the flow at the outlet.
• Falling Limb (Recession Curve): The descending portion extending from the peak to the normal base flow level. Represents the withdrawal of stored water. Its shape depends entirely on catchment physical properties.
• Inflection Point: The point on the falling limb where the curve changes concavity. Roughly indicates when direct surface runoff ceases and groundwater becomes primary.
• Base Flow: The delayed flow reaching the stream as groundwater. Sustains river flow during dry periods.

(Sketch: Time on X-axis, Discharge on Y-axis. Draw a bell-shaped curve. Label Rising Limb, Peak/Crest, Falling Limb, Inflection Point. Draw a base line separating Direct Runoff (above) from Base Flow (below).)

Base Flow: The portion of streamflow that comes from groundwater seepage into the channel. It is the delayed flow that sustains the stream during periods without precipitation.

Methods of Base Flow Separation:

1. Straight Line Method (Method I): Draw a straight line connecting the point where the rising limb begins (Point A) to a point on the falling limb (Point B). Point B is located using: $N = 0.83 \times A^{0.2}$ days from peak (where $A$ is drainage area in $km^2$). Suitable for small catchments.
2. Fixed Base Length Method (Method II): Extend the base flow curve from the start of the rising limb to directly under the peak (Point C). From Point C, draw a straight line to Point B on the falling limb (calculated using $N = 0.83 \times A^{0.2}$). More realistic than Method I.
3. Variable Slope Method (Method III): The initial base flow recession curve is extrapolated forward to a point directly under the peak (Point C). The recession curve from the latter part of the falling limb is extrapolated backward to a point under the inflection point (Point D). Points C and D are joined by a smooth curve. The most realistic method.

Unit Hydrograph (UH): The Direct Runoff Hydrograph (DRH) resulting from 1 cm (or 1 inch) of effective rainfall generated uniformly over the entire catchment area at a uniform rate during a specified duration ($D$ hours).

Basic Assumptions:

1. Time Invariance: The shape of the UH is independent of the time of the storm.
2. Linear Response (Proportionality): If 1 cm of ER produces a certain DRH, then ‘n’ cm in the same duration will produce a DRH with ordinates ‘n’ times the UH ordinates.
3. Superposition: The total DRH is the sum of the individual DRHs produced by each burst of rain, appropriately shifted in time.
4. Uniform Intensity: Effective rainfall intensity is uniform throughout the specified duration.
5. Uniform Distribution: Effective rainfall is distributed uniformly over the entire catchment area.

Limitations:

• Not applicable for very large catchments ($> 5000 \text{ km}^2$) or very small catchments ($< 2 \text{ km}^2$).
• The theory fails if rainfall is highly concentrated in one part of the basin.
• Cannot be applied to runoff caused by snowmelt.
• In reality, hydrologic systems are somewhat non-linear; larger floods may travel faster than predicted.

To derive a D-hour Unit Hydrograph from an isolated storm:

1. Select a suitable isolated storm: Choose a storm with relatively constant intensity for a duration $D$ that occurred somewhat uniformly over the catchment. Obtain the streamflow hydrograph.
2. Base Flow Separation: Subtract base flow ordinates from the total hydrograph to get the Direct Runoff Hydrograph (DRH).
3. Calculate Volume of Direct Runoff:
   $$V = \sum O_i \times \Delta t$$
4. Calculate Effective Rainfall Depth ($R_e$):
   $$R_e = \frac{V}{A} \text{ (in cm)}$$
5. Calculate Unit Hydrograph Ordinates: Divide each DRH ordinate by the effective rainfall depth:
   $$U_i = \frac{O_i}{R_e}$$
6. Determine Duration (D): From the rainfall hyetograph, identify the duration of the effective rainfall period.

S-Curve: The hydrograph of direct surface runoff resulting from a continuous, uniform effective rainfall at a constant rate of 1 cm/hr for an infinite duration. It eventually reaches a maximum equilibrium discharge ($Q_e$) where runoff rate equals rainfall rate.

Construction:

1. Start with a known $D$-hour Unit Hydrograph.
2. Shift the $D$-hour UH successively by $D$ hours (at $t=0, t=D, t=2D$, etc.).
3. Sum the ordinates of all these shifted UHs at any given time $t$: $S_t = u_t + S_{t-D}$.

Application in Deriving UHs of Different Durations:

1. Derive the S-curve ($S_A$) from the given $D$-hour UH.
2. Shift the entire S-curve by the desired duration $T$ to get a second S-curve ($S_B$).
3. Subtract: $\Delta y = S_A – S_B$. This represents a DRH from ER of depth $T/D$ cm.
4. Multiply by $D/T$ to obtain the $T$-hour Unit Hydrograph:
   $$T\text{-hour UH ordinate} = \Delta y \times \frac{D}{T}$$

Synthetic Unit Hydrograph: For catchments where historical rainfall and streamflow records are insufficient (ungauged catchments), empirical equations are used to generate a unit hydrograph based on the basin’s physical and topographical characteristics.

Snyder’s Method (Key Parameters):

• Basin Lag ($t_p$):
  $$t_p = C_t \left(L \times L_{ca}\right)^{0.3}$$
  Where $L$ = length of main stream to divide (km), $L_{ca}$ = distance to catchment centroid (km), $C_t$ = regional constant.
• Standard Duration ($t_r$): $t_r = t_p / 5.5$
• Peak Discharge ($Q_p$):
  $$Q_p = \frac{2.78 \times C_p \times A}{t_p}$$
  Where $A$ = catchment area ($km^2$), $C_p$ = regional constant.
• Base Time ($T_b$): $T_b = 3 + t_p / 8$ (in days).
• Adjusted Lag for Non-Standard Duration $t_R$: $t’_{p} = t_p + (t_R – t_r)/4$. Adjusted peak: $Q’_{p} = 2.78 C_p A / t’_{p}$.
• Widths ($W_{50}$ and $W_{75}$):
  $$W_{50} = \frac{5.87}{q^{1.08}}, \quad W_{75} = \frac{W_{50}}{1.75}$$
  Where $q = Q_p/A$.

Time of Concentration ($t_c$): The time required for a drop of water to travel from the most hydraulically remote point of the catchment to the outlet. It represents the time taken for the entire catchment area to contribute to the flow at the outlet.

Effect on Hydrograph Shape and Peak:

• When Rainfall Duration $\ge t_c$ ($t_r \ge t_c$): The storm lasts long enough for runoff from the farthest point to reach the outlet while rain is still falling. The entire catchment contributes simultaneously. Result: Maximum possible peak runoff. Steep rising limb with a sharp, high peak.
• When Rainfall Duration $< t_c$ ($t_r < t_c$): Rain stops before runoff from the most distant parts reaches the outlet. Only a partial area contributes at any given time. Result: Attenuated (lower) peak discharge. Flatter rising limb, lower peak, and wider time base.

1. Design Flood: The flood discharge adopted for the design of a specific hydraulic structure to ensure its safety and economic viability. Selected after balancing the cost of construction against the potential risk and damage if exceeded. Different structures use different return periods based on their importance.
2. Standard Project Flood (SPF): The flood resulting from the most severe combination of meteorological and hydrological conditions considered reasonably characteristic of the geographic region. Excludes extremely rare combinations. Typically used for major hydraulic structures and is generally about 40–60% of the PMF.
3. Probable Maximum Flood (PMF): The extreme maximum flood that can ever physically occur in a river basin. Results from the most severe possible combination of meteorological and hydrological conditions. Used strictly for spillway design of very large dams where structural failure is unacceptable.

1. Return Period ($T$): Also known as the recurrence interval, it is the average time interval (usually in years) between the occurrence of a flood of a specified magnitude and an equal or larger flood. For example, a “100-year flood” is expected to be equaled or exceeded on average once every 100 years.
2. Probability of Exceedance ($P$): The probability that a flood of a specified magnitude will be equaled or exceeded in any single, given year. It is the reciprocal of the return period:
   $$P = \frac{1}{T}$$
   For a 100-year flood, $P = 0.01$ or 1% chance in any given year.

1. Risk ($R$): The probability that a flood of a given magnitude will be equaled or exceeded at least once during the entire design life of the structure:
   $$R = 1 – (1 – P)^n$$
   Where $P$ = probability of exceedance in a single year, $n$ = design life in years.
2. Reliability ($Re$): The probability that the structure will perform without failure over its design life:
   $$Re = 1 – R = (1 – P)^n$$
3. Safety Factor: An additional margin of capacity (such as adding “freeboard” to a dam’s height) incorporated into the final design to account for unforeseen extreme conditions, data inaccuracies, or changing climatic behaviors beyond the statistical calculation.

Definition: A statistical method used to estimate the magnitude of floods corresponding to specific return periods by fitting observed historical peak discharge data to a probability distribution (such as Gumbel or Log-Pearson).

Significance:

• Predictive Design: The primary tool to estimate design floods required for engineering structures.
• Objective Approach: Relies on actual observed historical streamflow data of the specific catchment.
• Economic Assessment: Allows evaluation of cost-benefit ratios of different structural sizes based on probabilistic risks.

Limitations:

• Data Dependency: Requires a long, reliable, and continuous record (usually >30 years). Short records lead to erroneous estimates.
• Stationarity Assumption: Assumes underlying hydrological conditions remain constant over time. Fails to account for future climate changes or land-use changes.
• Extrapolation Risks: Extrapolating a 1000-year flood from only 40 years of data is highly uncertain.

1. Data Collection: Extract the maximum peak discharge recorded in each year to form the “Annual Maximum Series” ($x_1, x_2, …, x_N$).

2. Calculation of Statistical Parameters:

3. Application of Chow’s General Frequency Equation:

Where $K$ = Frequency Factor (depends on return period $T$ and chosen probability distribution).

4. Selection of Probability Distribution: Choose Gumbel’s Extreme Value Distribution or Log-Pearson Type III Distribution to determine $K$.

5. Calculation of Design Flood ($x_T$): Substitute $K$ into Chow’s equation to yield the design flood magnitude.

Gumbel’s method uses Chow’s general frequency equation: $X_T = \bar{X} + K \cdot \sigma_{n-1}$.

Steps and Components:

1. Calculate mean $\bar{X}$ and standard deviation $\sigma_{n-1}$ from the annual maximum series.
2. Calculate Frequency Factor ($K$):
   $$K = \frac{y_T – \bar{y}_n}{S_n}$$
   Where:
   • $y_T$ = Reduced Variate: $y_T = -\ln \left[ -\ln \left( 1 – \frac{1}{T} \right) \right]$
   • $\bar{y}_n$ = Reduced Mean (from Gumbel tables, depends on $N$; as $N \to \infty$, $\bar{y}_n \to 0.577$)
   • $S_n$ = Reduced Standard Deviation (from Gumbel tables; as $N \to \infty$, $S_n \to 1.2825$)
3. Substitute $K$ into the main equation to find the design flood magnitude $X_T$.

The Log-Pearson Type III distribution accounts for the “skewness” of flood data, which Gumbel’s method ignores.

Procedure:

1. Logarithmic Transformation: Convert all observed annual maximum flood values to base-10 logarithms: $z = \log_{10}(x)$.
2. Calculate Statistical Parameters: Calculate mean ($\bar{z}$), standard deviation ($\sigma_z$), and coefficient of skewness ($C_s$) of the transformed series:
   $$C_s = \frac{N \sum(z – \bar{z})^3}{(N-1)(N-2)(\sigma_z)^3}$$
3. Apply the General Equation:
   $$z_T = \bar{z} + K_z \cdot \sigma_z$$
   Where $K_z$ is the frequency factor obtained from standard Log-Pearson Type III tables based on $T$ and $C_s$.
4. Final Computation: Take the antilog to find the actual flood magnitude:
   $$X_T = 10^{z_T}$$

The Rational Method: A widespread empirical equation to estimate peak surface runoff rate from urban areas and small rural catchments:

Where $Q$ = Peak discharge ($m^3/s$), $C$ = Runoff coefficient, $I$ = Average rainfall intensity (mm/hr) for a duration equal to $t_c$, $A$ = Catchment area in hectares.

Fundamental Assumptions:

1. Rainfall intensity is uniform over the entire catchment area.
2. Rainfall intensity remains constant for a duration equal to or greater than $t_c$.
3. Peak discharge occurs when the entire catchment is contributing to the outlet (at $t = t_c$).
4. The runoff coefficient ($C$) remains constant throughout the storm.

Limitations:

• Only applicable to small catchments (typically under $50 \text{ km}^2$) because uniform rainfall is physically unrealistic over large areas.
• Gives only the peak discharge, not the entire flood hydrograph.
• Selecting the correct Runoff Coefficient ($C$) is highly subjective.

Flow Routing (Flood Routing): A mathematical or empirical procedure used to predict the changes in the magnitude, speed, and shape of a flood wave (hydrograph) as it propagates down a river channel or through a reservoir. Given an upstream inflow hydrograph, routing determines the outflow hydrograph at a point downstream.

Difference between Routing and Flood Tracking:

• Flood Routing: A rigorous hydrological or hydraulic process that computes the exact shape, peak magnitude, attenuation, and time base of the hydrograph at a downstream location using mathematical models. Gives a complete picture of flow variations over time.
• Flood Tracking (Flood Tracing): Generally refers to observing or tracking the movement of a flood wave (specifically the flood crest) geographically over time. Often observational or relies on simple travel-time correlations for early warning systems, rather than mathematically transforming the entire hydrograph.

Basic Equations:

Hydrologic Routing — Continuity Equation:

Hydraulic Routing — Saint-Venant Equations:

1. Continuity: $\frac{\partial A}{\partial t} + \frac{\partial Q}{\partial x} = q$

2. Momentum: $\frac{\partial y}{\partial x} + \frac{V}{g}\frac{\partial V}{\partial x} + \frac{1}{g}\frac{\partial V}{\partial t} = S_0 – S_f$

• Attenuation: As a flood wave travels downstream, a portion of incoming water is temporarily stored in the channel or reservoir, causing the peak of the outflow hydrograph to be lower than the peak of the inflow hydrograph. Attenuation = $I_p – O_p$.
• Lag of Peak: Because it takes time for the reservoir/channel to fill up and release stored water, the occurrence of peak outflow is delayed compared to peak inflow. This time difference is the “lag of peak.”

(Sketch: Time on X-axis, Discharge on Y-axis. Draw a steep “Inflow” hydrograph curve and a second flatter “Outflow” curve shifted to the right. Mark the highest point of Inflow as $I_p$ and Outflow as $O_p$. The vertical difference between $I_p$ and $O_p$ is “Attenuation”. The horizontal time difference between the peaks is the “Lag of Peak.”)

Starting from the continuity equation: $I – O = \frac{dS}{dt}$. Integrating over a time interval $\Delta t = t_2 – t_1$:

Rearranging to group known quantities (left) and unknown quantities (right):

This is the standard reservoir routing equation used in the Modified Puls Method. The right-hand side contains all unknowns at time $t_2$; the left-hand side is entirely known.

Hydrologic Reservoir Routing: Tracks a flood wave as it passes through a reservoir using the continuity equation. Since the water surface is assumed horizontal (level pool routing), storage $S$ is purely a function of the water surface elevation ($h$), and outflow $O$ is also a function of $h$. Therefore, $S = f(O)$.

Modified Pul’s Method Procedure:

1. Data Preparation: Obtain Elevation-Storage and Elevation-Outflow relationships. Combine them to get a Storage vs. Outflow ($S$ vs $O$) relationship.
2. Routing Curve: Select a routing time interval $\Delta t$. For various values of $O$, compute the corresponding $S$. Construct a table/plot of $O$ versus $\left( S + \frac{O\Delta t}{2} \right)$.
3. Initial Conditions: At $t=0$, $I_1, O_1,$ and $S_1$ are known. Compute: $\left( S_1 – \frac{O_1 \Delta t}{2} \right)$.
4. Routing Step: Calculate: $\left( S_2 + \frac{O_2 \Delta t}{2} \right) = \left( \frac{I_1 + I_2}{2} \right) \Delta t + \left( S_1 – \frac{O_1 \Delta t}{2} \right)$.
5. Finding Outflow: Read the corresponding outflow $O_2$ from the routing curve.
6. Update and Repeat: Compute $\left( S_2 – \frac{O_2 \Delta t}{2} \right) = \left( S_2 + \frac{O_2 \Delta t}{2} \right) – O_2 \Delta t$ and continue for the next time step.

Muskingum Method: Models the “wedge” storage in a channel as both inflow and outflow dependent:

Where $K$ = Storage time constant, $x$ = Weighting factor ($0 \le x \le 0.5$). Substituting into the discrete continuity equation yields the routing equation:

Where:

(Note: $C_0 + C_1 + C_2 = 1$)

Determination of ‘K’ and ‘x’: Observed inflow and outflow hydrographs for a historical flood in the given reach are required.

1. Compute cumulative storage $S$ at each time step using continuity.
2. Assume trial values of $x$ (from 0.1 to 0.5).
3. For each trial $x$, compute the weighting term $[xI + (1-x)O]$ at every time step.
4. Plot $[xI + (1-x)O]$ on the X-axis versus Storage ($S$) on the Y-axis.
5. The correct value of $x$ causes the loop to collapse as closely as possible to a single straight line.
6. The slope of that resulting straight line is the value of $K$.

Linear Reservoir: A conceptual storage element where the storage is directly and linearly proportional to outflow: $S = kO$, where $k$ is the storage constant. Using the continuity equation, the routing equation for a linear reservoir is:

Clark’s Unit Hydrograph Concept: Clark’s method separates two main physical processes:

1. Translation (Lag): Movement of water from origin to outlet, simulated using the Time-Area Method.
2. Attenuation (Storage): Temporary retention of water, simulated by routing the translated flow through a single conceptual linear reservoir at the catchment outlet.

Procedure:

1. Construct the Time-Area Histogram: Draw isochrones (lines of equal travel time) on the catchment map at intervals of $\Delta t$. Measure the inter-isochrone areas ($A_1, A_2, …$) to form the Time-Area Histogram.
2. Generate the Inflow Hydrograph ($I$): Apply 1 cm of instantaneous effective rainfall uniformly over the catchment. Convert the time-area histogram into a hypothetical inflow hydrograph. Inflow: $I = \frac{2.78 \cdot A}{\Delta t}$ $m^3/s$ (for $A$ in $km^2$, $\Delta t$ in hours).
3. Route through a Linear Reservoir: Determine the storage constant $K$ of the catchment from the recession limb of an observed hydrograph. Route using:
   $$O_2 = C \cdot I_{avg} + (1 – C) \cdot O_1, \quad \text{where } C = \frac{2 \Delta t}{2K + \Delta t}$$
   The resulting outflow hydrograph $O$ is the Instantaneous Unit Hydrograph (IUH).
4. Convert IUH to Unit Hydrograph: Apply the S-curve technique to get a specific $D$-hour Unit Hydrograph.