Chapter 4 Surface Runoff and Flow Measurement 2024.pdf

Stream discharge measurement, current meter, ADCP, slope-area method, salt dilution method, chemical method of discharge measurement, surface velocity method

1. Chapter 4: Surface Runoff and Flow
Measurement
(5 hours)
Prof. Dr. Hari Krishna Shrestha
Contact: [email protected]
Last update: April 18, 2024
Chapter 4: Surface Runoff and Flow Measurement (5 hrs) (17%)
1.Drainage basin and its quantitative characteristics
2.Factors affecting surface runoff
3.Rainfall-runoff correlation (linear)
4.Stream gauging, selection of site, types of gauges and their selection
5.Stream flow measurement
1.Velocity area method, current meters, floats, velocity rods & dilution
techniques
2.Slope area method
6.Development of rating curve and its uses

2. 4.0 Surface Runoff- Process
•Precipitation
•Infiltration
•Percolation
•Sub-surface flow
•Interflow
•Surface flow – runoff
•Sheet flow
•Concentrated flow – rivulet, brook, stream, river

3. Precipitation
Interception
Evaporation
Transpiration
Stem Flow
Runoff
Throughfall
Infiltration
Subsurface Flow
Evaporation
Uptake
Runoff
Streams
Rivers

4. 4.1 Drainage Basins and its Quantitative Characteristics
•A drainage basin (river catchment) is an area of land drained by a river and its
tributaries; when it rains in this area, the water goes towards the main river and
ends up at the river’s mouth.
•Quantitative Characteristics:
1.Catchment area and hypsometric curve
2.Catchment Centroid (x, y) in degrees/meters or northing/easting, z in meters (mamsl)
3.Catchment Slope, maximum and minimum elevations of the catchment divide
4.Land use/Land Cover (LULC) and local depression areas
5.Average infiltration for different return-period rainfall events, and Runoff Coefficient
6.Length and Slope of the Main River and major tributaries
7.Drainage Density and order of the tributaries
8.Surface soil type and the area of each soil type
9.Spatial distribution of rainfall in the catchment
10.Time series data of rainfall in the catchment (average, maximum and minimum)
11.Evapotranspiration from the catchment
12.Time series data of radiation, temperature, humidity, pressure, and wind speed
13.Intensity-Duration-Frequency and Depth-Area-Duration curves of the catchment
14.Groundwater recharge rate and volume in the catchment
15.Long term flow, peak flood flow, minimum flow from the river mouth of the catchment
16.Sediment type and discharge of each sediment type out of the catchment
Additional info: https://www.youtube.com/watch?v=jPF8oXSEx4s&ab_channel=learnsomethingtoday
http://www.ijstm.com/images/short_pdf/1415296455_P39-50.pdf

5. 4.2 Factors Affecting Surface Runoff
Kaligandaki River

6. The braids on the Yarlung Tsangpo (Brahmaputra) near Lhasa before it enters a series of gorges
where the world's largest hydroelectric plant is being built by China.

7. 4.2 Factors Affecting Surface Runoff
•Catchment factors
•Basin size, shape, slope
•Nature of the valley: wide, narrow
•Elevation
•Drainage density
•Infiltration factors
•Land-use and land-cover
•Soil type & geological conditions
•Depression storages
•Channel characteristics
•Cross section
•Roughness:
–river bed, river banks
•Storage capacity
Meteorological Factors:
•Storm Characteristics
•Initial loss
•Evapotranspiration
CatchmentBasin

8. 4.3 Rainfall-Runoff Correlation (linear)
•Various forms of rainfall-runoff relations can be
developed based on specific rainfall and associated
runoff data
•Linear relationship: Q = a P + C
•Non-linear relationship: Exponential/Power/
Polynomial
–Q = discharge (weekly, monthly) from Direct Runoff
Hydrograph
–P = Precipitation

9. 4.3 Rainfall-Runoff Correlation
Linear type relation
•R = a P + b; R = runoff, P = precipitation, a and b are linear regression constants
•a = [N(PR) – (P)(R)]/[N(P
2
)-(P)
2
]
•b = [R – a(P)]/N
•�=
[N(PR) – (P)(R)]
[N(P
2
)−(P)2] [N(R
2
)−(R)2]
;�=??????����??????��??????�� ??????����??????????????????���
Exponential type relation
•R = P
m
or
•ln R = m ln P + ln 
Multiple Regression
Y = b
0 + b
1 x
1 +b
2 x
2 + b
3 x
3 + b
4x
4
Y = runoff, x
1 = baseflow, x
2 = autumn precipitation, x
3= snow water equivalent, x
4 = spring
precipitation, and b
i = regression coefficients.
Unless base flow can be separated and effect of snow melt can be established (if applicable),
the relation between daily data cannot be made. Normally relation between monthly, seasonal
or annual data is made. The ratio of runoff to rainfall is Runoff Coefficient.
Steps:
Convert weekly/monthly/annual flow into depth or volume.
Convert precipitation into depth or volume.
Plot the data and find the relation using appropriate regression method.

10. https://www.researchgate.net/publication/3190
37166_Water_balance_study_of_Beas_river_Hi
machal_Pradesh_using_ARCGIS_technique_upt
o_Pong_dam/figures?lo=1
Similar Rainfall-
Runoff Relationship
are established on
monthly or seasonal
basis to quickly
estimate the runoff
from rainfall data.

12. 4.4 Stream Gauging (site selection, gauge type)
Runoff gauged indirectly through staff gauging (stage monitoring).
Gauge Type:
A) Discreet (Non-recording)
•Staff Gauge
•Sectional Staff Gauge
•Crest Gauge
B) Continuous (Recording)
•Laser/Radar Gauge (RLS)
•Automatic Staff Gauge
•Pressure Gauge
•Bubble Gauge
Continuous data required for generation and analysis of hydrograph.
The updated number of rivers and rivulets in Nepal, Energy Development
Commission (June 2016), is 11614.
http://www.renewableenergyworld.com/articles/2016/06/nepal-seeks-investors-for-
10-gw-of-electricity-by-2026.html

13. Sample annual gauge
river gauge versus date/year0
100
200
300
400
500
600
7000
200
400
600
800
1000
1200
1400
2017/
Jan
2017/
Apr
2017/
Jul
2017/
Oct
2018/
Jan
2018/
Apr
2018/
Jul
2018/
Oct
2019/
Jan
2019/
Apr
2019/
Jul
dail;y Rainfall [mm/day]
Discharge
[m
3
/s]
rain
obs.Q by H-Q
calc.Q
Measured Q
Snow Depth 0
100
200
300
400
500
600
7000
200
400
600
800
1000
1200
1400
2005/
Jan
2005/
Apr
2005/
Jul
2005/
Oct
2006/
Jan
2006/
Apr
2006/
Jul
2006/
Oct
2007/
Jan
2007/
Apr
2007/
Jul
2007/
Oct
2008/
Jan
2008/
Apr
2008/
Jul
2008/
Oct
dail;y Rainfall [mm/day]
Discharge
[m
3
/s]

14. This technology is gradually getting obsolete, replaced by radar/laser water level sensors.
However, its use for the calibration purpose is still intact.

15. This technology is gradually getting obsolete, replaced by radar/laser water level sensors.

16. This technology is gradually getting obsolete, replaced by radar level sensors (RLS).
Bubble Gauge

17. RLS for continuous river gauging RLS at Rui Khola, Chitwan for EWS

18. 4.4 Stream Gauging (Stage Monitoring)
staff gauge, sectional staff gauge, automatic staff gauge: pressure type, float
type, laser
Staff and Crest Gauge
Automatic Gauge recorder Automatic Gauge recorderLaser type Automatic Gauge recorder
Crest Gauge

20. The U.S. Geological Survey (Rantz et al., 1982) have developed nine criteria for an "ideal" gaging
site. The criteria are:
1.The stream course is straight for about 300 feet upstream and downstream of the gage site.
2.At all stages, the total flow is confined to a single channel. There is also no subsurface or
groundwater flow that bypasses the site.
3.The streambed in the vicinity of the site is not subject to scour and fill. It is also free of
aquatic plants.
4.The banks of the stream channel are permanent. They are free of brush and high enough to
contain floods.
5.The stream channel has unchanging natural controls. These controls are bedrock outcrops or
stable riffle for low flow conditions. During high flows, the controls are channel constrictions or
a cascade or falls that is not submerged at all stages.
6.At extremely low stages, a pool is present upstream from the site. This will ensure the
recording of extremely low flows and avoid the high velocities associated with high stream
flows.
7.The gaging site is far enough removed from the confluence with another stream or from tidal
effects to avoid any possible impacts on the measurement of stream stage.
8.Within the proximity of the gage site, a reach for the measurement of discharge at all
stages is available.
9.The site is accessible for installation and operation and maintenance of the gaging site. The
selection of a gaging site is again a compromise between these criteria.
4.4 Site Selection (for establishing gauge site)

21. Site Selection (for establishing gauge site)
Characteristics of an ideal site for river gauging:
Parameter Reason
Straight reach
No turbulence
No drops
No sharp bends
No backwater effect
Accessible/visible
No direct impact from flow
Stable water surface
Stable cross section
Stable river bed

22. 4.5 Stream Flow Measurement by V-A methods
4.5.1 Current Meter: vertical axis, horizontal axis
4.5.2 Surface Float
4.5.3 Velocity Rods
4.5.4 Dilution (salt dilution)
•Hydraulic Structures
•Electromagnetic
•Acoustic Doppler Current Profiler
(ADCP)

23. ← Vertical Axis Current Meter
Horizontal Axis Current Meter ↓
4.5.1 Current Meters

24. 4.5.1 Stream Flow Measurements: Cup type current meter

25. Stream Discharge Measurements using cup-type current meter
Head phone to count number of revolutions Lead weight to keep hand-line vertical and
current meter stationary
Cable car for Q measurement in bigger rivers Stilling well for stage hydrograph

26. 4.5.2:
4.5.3: Velocity

27. Stream Discharge Measurements
Winch to control position of cable car
A-frame to elevate cable of cable car
Surface velocity measurement

28. Acoustic Doppler
Current profiler
(ADCP)

30. Discharge (m
3
/s, cumec)
•Discharge = cross sectional area of flow × average river flow velocity
•Q = V × A [unit: m
2
* m/s = m
3
/s] cusec?
V = ƒ (number of rotation, time) (traditional cup-type or propeller type
current meter); higher the velocity higher the rotation speed
Stream Flow Measurement methods (Direct and Indirect):
•Current Meter: vertical axis, horizontal axis (based on the axis along
which the cups/blades of current meter rotates)
•Surface Float: Average flow velocity = surface velocity * K
•Chemical – dilution, salt dilution: mostly for small rivers with high
turbulence. Why?
•Electromagnetic
•Acoustic/Acoustic Doppler Current Profilers (ADCP)
•Hydraulic (or gated) Structures: Q = ƒ (h) = C h
n
(weir/ notch)
•Slope Area Method (Section 4.6)

31. Some newer types of current meters display the flow velocity directly, without having to
count the number of revolutions. However, regular calibration of the current meters is
needed to ensure accuracy of the discharge measurement.

35. Total Discharge = 7.065 m
3
/s
Total area = 20.6 m
2

Average Velocity = 0.343 m/s

36. Distance
from bank
(m)
Width
W(m)
Effective
Width (m)
Depth
(m)
Revolutio
n
Time
(sec)
Area
(m
2
)
Velocity
(m/sec)
Discharge
(m
3
/sec)
0.0 0
0.6 0.60.675 1 15 500.6750.290.196
1.2 0.60.70 4 39 552.80.6741.887
2.0 0.80.90 5.5 50 504.950.854.208
3.0 1 0.90 6.5 56 505.850.9465.534
3.8 0.80.75 4.5 39 503.3750.6742.275
4.5 0.70.60 2.5 35 501.50.610.915
5.0 0.50.602 1 20 500.6020.370.223
5.6 0.6 0 0
19.750.7715.237
Effective width (1
st
) = (W
1 + (W
2/2))
2
/(2W
1)
Effective width (last) = (W
n+ (W
n-1/2))
2
/(2W
n)
V = aN+b, a = 0.016, b = 0.05
Average (stream flow) velocity = Q / A
Numerical example of a stream discharge measurement calculation.

37. V= a N + b
a = 0.0113; b = 0.0059
N= no. of revolution in 60 sec
(in this example)

38. Sample Discharge measurement and
calculation sheet.

39. 4.5.2 Stream Flow Computation by Slope Area Method
https://pubs.usgs.gov/twri/twri3-a2/pdf/twri_3-A2_a.pdf
https://www.yourarticlelibrary.com/water/river-training/slope-area-method-concept-and-selection-of-reach/60961
Indirect method of flow estimation, assuming uniform flow, using Manning’s Formula
Energy equation:
Z + y = h (water surface elevation above the datum; h
L = h
f + h
e
ℎ
1+
??????
1
2
2�
=ℎ
2+
??????
2
2
2�
+ℎ
�+ℎ
� → ℎ
�=ℎ
1−ℎ
2+
??????
1
2
2�
−
??????
2
2
2�
− ℎ
�
�
1+??????
1+
??????
1
2
2�
=�
2+??????
2+
??????
2
2
2�
+ℎ
??????
ℎ
�
�
=�
�=�����?????? �??????���=
�
2
�
2
; �= ??????ℎ����?????? ??????��??????�??????��??????�=
1
�
??????�
ൗ
2
3
Average K for the reach = (k
1*k
2)
0.5
h
f = fall + (V
1
2
/2g - v
2
2
/2g) when h
e ≈ 0
ℎ
�= �
�
??????
1
2
2 �
−
??????
2
2
2 �
K
e = 0.3 for gradual expansion and
0.1 for gradual contraction

40. During a flood flow the depth of water in a 10m wide rectangular channel was found to
be 3.0 m and 2.9 m at two sections 200 m apart. The drop in the water-surface elevation
was found to be 0.12 m. Taking Manning's coefficient to be 0.025, estimate the flood
discharge through the channel.
L = 200m W = 10m
h1 = 3m h2 = 2.9m
head loss=0.12m Manning's n =0.025Flood discharge = ?
A
1
= 30 A
2
=29 A = W h
P
1
= 16 P
2
=15.8 P = W + 2 h
R
1
= 1.875 R
2
=1.835 R = A/P
K
1
= 1824.7 K
2
=1738.9K = (1/n)A R
2/3

Average K for the reach = (k
1
*k
2
)
0.5
= 1781.3
h
f
= fall + (V
1
2
/2g - v
2
2
/2g) = 0.12 + (V
1
2
/2g - V
2
2
/2g)
Trial h
f
S
f
QV
1
2
/2gV
2
2
/2g h
f
1 0.120.0006043.630.10780.11540.1124
2 0.11240.0005642.230.10100.10810.1129
3 0.11290.0005642.320.10140.10860.1129
Note: h
f value for Trial 1 is assumed to be same as head loss value of 0.12m. For
subsequent trials, it is copied from the last column. The final Q value is achieved
when h
f values does not change.

41. 4.6 Rating Curves and
its Uses
• Rating Curve: relationship
between stage and discharge
• Developed by conducting a
series of stream discharge
measurement
• The curve can change if the
cross section of the river change
due to scouring, river bed
aggradation, or bank cutting
Hysteresis in rating curve: different discharge for same gauge
height and vice versa, When a flood wave propagates through a river
corresponding to same stage higher discharges are observed during
rising stage than in falling stages resulting in looped rating curves. This
affect is popularly known as hysteresis in stage–discharge relationship

42. Data to plot rating curve of a river
SNDischargeGH
140.6041.94
2 54.552.11
3 59.842.14
462.0052.18
573.6882.28
6 76.062.27
780.1122.29
886.4852.40
989.8082.44
10100.6212.53
11116.852.64
12129.4432.73
13155.852.92
14221.5013.28
15559.323.80
16599.4784.00
17654.7614.40
18919.5345.28
GH= 0.6453Q
0.2946
R² = 0.9831
0
0.5
1
1.5
2
2.5
3
3.5
4
4.5
5
5.5
6
0 100 200 300 400 500 600 700 800 9001000
G.H (m)
Discharge (m
3
/sec)
GH= 0.6453Q
0.2946
R² = 0.9831
1
10
1 100
G.H (m)
Discharge (m
3
/sec)

44. 4.6 Uses of Rating Curves
•Convert stage into discharge
•Flood Disaster Risk Management: Early Warning
System
•Flood inundation area demarcation
•Gate operation in water resources development
(WRD) projects (hydropower/Irrigation, …)
•Sediment diversion during flood events
•Monitoring minimum flow release from WRD
projects for navigation and socio-cultural and
environmental uses.

45. Discuss the practical uses of rating curve of a river section. How is a rating curve
developed? Why a same river section can have multiple rating curves?
Taking the rating of current meter as V = 0.03 + 0.8 N, where V is in m/sec and N is the
number of revolutions/sec, compute the stream flow (Q) from the given data. If the rating
curve of the river section can be approximated by log Q = 1.1 log S + 0.2, calculate the
discrepancy % in measurement, given the river stage (S) during flow measurement is 1.6 m.
Distance from bank (m)0.00.6 1.5 2.5 5.0 7.07.5
Depth (m) 0 0.3 0.75 1.2 1.2 0.30
Revolution 0 9 25 1530 163016 50
Time (sec) 0 45 90 801009010080400

46. (h) What are the reasons for a hysteresis loop in a rating curve of a river section?
(i) Discuss the role of infiltration capacity index in determining the excess runoff depth from
precipitation data.
(k) Given the stage (GH) versus discharge (Q) data, develop a rating curve equation and
calculate the discharge when GH = 3.5 m.
GH
(m)
1.942.112.142.182.282.272.292.402.442.532.642.732.923.283.804.004.405.28
Q
(m
3
/s)
40.6054.5559.8462.0073.6976.0680.1186.4889.81100.6116.8129.4155.8221.5559.3599.5654.8919.5

47. Expected skills from this chapter:
1.Develop rainfall-runoff relation (equation) from a set of rainfall and runoff data.
2.Stream discharge measurement using a current meter.
3.Calculate stream discharge from discharge measurement data (distance from edge,
depth, number of revolutions and time).
4.Develop rating curve from given set of stage and discharge data, using MS Excel’s
Solver facility.