Research Article Volume 1 Issue 6
Department of Water Resources Development & Management, Indian Institute of Technology Roorkee, India
Correspondence: Kansal ML, Department of Water Resources Development & Management, Indian Institute of Technology Roorkee, India
Received: August 19, 2017 | Published: November 27, 2017
Citation: Kansal ML, Thakur A. Flood estimation for a cloudburst event in an ungauged western himalayan catchment. Int J Hydro. 2017;1(6):163-168. DOI: 10.15406/ijh.2017.01.00030
Cloud burst is one of the reasons of flash flood in hilly regions. Such events have caused heavy loss of life and property due to its devastating nature in terms of landslides and boulder movement along with shooting velocity of water which results in heavy change of momentum and hence the force. Often, these flash floods are associated with very small areal extent, less duration of storm and remoteness of location. It is difficult to predict the location and impact of cloud burst particularly when most of these hilly areas are ungauged. Therefore, such events sometime go unnoticed and even the events which are reported, have scanty data due to the poor ground monitoring mechanism. In India, there are 3 major states in North and 7 states in north-east which are termed as ‘Himalayan states’. During the last few years, these states have experienced a number of cloud burst events which may be attributed to the increased anthropogenic activities and the climate change phenomena. Keeping this in mind, the present study aims at estimating the outflow flood hydrograph at the outlet of an unguaged Himalayan catchment using Synthetic Unit Hydrograph. Some of the historical cloudburst events in Northern Himalaya are studied and accordingly the runoff and basin lag corrections are suggested. The results are compared with the previous study reported in literature. Such study is expected to help in disaster preparedness. Further, it highlights the need for comprehensive data base of such events over the entire Himalayan range for effective disaster management in the region.
Keywords: cloudburst, flash flood, snyder’s method
Cloud burst is one of the reasons of flash flood in hilly regions. Such events are devastating as these induce landslides and boulder movement along with shooting velocity of water in hilly areas. This results in huge change of momentum and hence the large force. This force along with flood movement causes heavy loss of life and property. Cloudburst is a sudden aggressive rainstorm of small duration (few minutes to few hours) with rainfall intensity of more than 100mm/hr.1 Some of the extreme rainfall events recorded in India range from 900mm/day to 1,040mm/day.2 In India, there are 3 major states (Jammu and Kashmir, Himachal Pradesh, and Uttarakhand) in North and 7 states Assam, Arunachal Pradesh, Manipur, Meghalaya, Mizoram, Nagaland and Tripura) in north-east which are termed as ‘Himalayan states’. During the last few years, these states have experienced a number of cloudburst events. Cloudbursts in Himalayan region represent the rapid convective lifting of moist air mass under the conditions of steep orography and thermodynamic instability. The warm and humid air moves uphill due to orography of the area. As this air mass continues to rise up, it forms large clouds. Lack of upper air at such great heights prevents the dissipation and water concentration in the cloud keeps on increasing which finally result in sudden localized downpour of water (Figure 1). Generally, cloud-burst is a localized weather phenomena concentrated over a small area (not exceeding 20-30 km2). Researchers have attributed the cloud burst phenomena to the increased anthropogenic activities and the climate change. Generally, it occurs in the rainy (monsoon) season (June to September). Often, these events are associated with very small areal extent, less duration of storm and remoteness of location. It is difficult to predict the location and impact of cloud burst particularly when most of these hilly areas are ungauged. Therefore, such events sometime go unnoticed and even the events which are reported, have scanty data due to the poor ground monitoring mechanism or sometimes even these also get washed away. Only data which remain available for such events is the flood water scars on the banks of streams, uprooted trees, and visual observations by the local people. These large number of eye witnesses provide crucial and reasonably accurate information of flood timing.3 Some of the recent cloud burst events occurred in the state of Uttarakhand in India and their estimated impact are as follows:
Keeping the importance of cloud burst in mind and its role in flash flood, it is desired to assess the quantum of flood caused by such an event in an unguaged catchment. This study focuses on Synthetic Unit Hydrograph approach for the north Himalayan regions of India. The suggested approach is validated through one of the cloud burst event. In future the proposed methodology can be extended to other regions and its reliability can be improved with availability of more data.
India has been classified in seven hydro-meteorological homogenous zones with different sub-zones. The Western Himalayan belt falls under zone 7 (Figure 2). This covers the states of Uttarakhand, Himachal Pradesh and Jammu & Kashmir. There is very less of no vegetation in such areas to offer resistance to flow. Also the land surface in these catchments generally has less infiltration index. So the flood hydrographs for these catchments are generally high peaked and narrow with very less time of peak. In this study, Leh cloudburst (2010) has been discussed and obtained results for the same have been compared with results of Renoj et al.4 Leh is a small town in Ladhakh region of the state of Jammu and Kashmir. Ladhakh covers 52.6% area of state of Jammu and Kashmir 4 and shares India’s northern boundary with Pakistan and China. The area has a long history of glaciations, and every subsequent glaciations in the region was smaller than the previous one.5 Due to adverse climatic conditions in the region, the settlements are sparsely situated. District Leh is situated roughly between 32 to 36 degree north latitude and 75 to 80 degree East longitude and altitude ranging from 2300 m to 5000 m above sea level. District Leh with an area of 45,100 km2 makes it 2nd largest district in the country after Kutch (Gujrat) with an area 45,652 km2 in terms of area. The district is bounded by Pakistan occupied Kashmir in the west and china in the north and eastern part and Lahul Spiti of Himachal Pradesh in the south east. Ladhakh lies on the rain shadow side of the Himalayan, where dry monsoon winds reaches Leh after being robbed of its moisture in plains and the Himalayas mountain the district combines the condition of both arctic and desert climate. Therefore Ladhakh is often called “COLD DESERT”. There is a wide diurnal and seasonal fluctuation in temperature with -40°C in winter and +35°C in summer. Precipitation in this region is very low with annual precipitation of 10cm mainly in form of snow. Air is very dry and relative humidity range from 6-24%.6
Quantum of flood in a catchment can be assessed in terms of peak, peak time and the time distribution of flow in terms of flood hydrograph. Flood estimation can be studied under following categories of methods:
Empirical method uses regional formulae based on correlation between discharge and catchment characteristics. Almost all the empirical formulae represent discharge as a function of catchment area. These formulae are regionally valid and give approximate results when applied to other catchments. For example, Dicken’s formula7 is widely used in India as the first approximate value of peak flood in a catchment, which is given as:
(1)
Where
Qp: Peak discharge in cumec; Cd: Varies from 6 to 30 (11.5 for northern India); A: Catchment area (km2)
Rational method another semi-empirical formula used for peak flood estimation in small catchments. It considers a rainfall of uniform intensity occurring over small catchments such that duration of rainfall is greater than the time of concentration for catchment. Runoff is assumed to increase gradually to a peak value at the time of concentration or time of peak and stay at the peak as long as the rainfall continues. Mathematically, flood peak and peak time (Kirpich equation) are represented as:
(2)
(3)
WhereIn flood frequency method, one uses the observed maximum annual flow values and estimates the statistical parameters of the probability distribution. This information is for estimation of likelihood of flood as a function of recurrence interval. But as this method is based on historical data, it cannot be used in case of cloudburst event as there is very scanty or almost no historical data available for such events. A unit hydrograph represents the lumped response of catchment to a unit rainfall excess of D-h duration to produce direct runoff hydrograph. The factors that affect the shape of unit hydrograph have been studied under two groups i.e. physiographic factors and storm characteristics. Physiographic factors include basin characteristics (shape, size, slope, nature of the valley, drainage density, elevation etc.), infiltration characteristics (land use and land cover, lakes and other storages, soil type and geological conditions) and channel characteristics (cross- sections, roughness, storage capacity etc.). Storm characteristics include precipitation, intensity, duration, magnitude and movement of storm. The effect of various factors on hydrograph is as follow: The shape of basin influence the time of concentration i.e. the time taken for water from the remotest part of catchment to reach the outlet. Fan shaped catchments gives high peak and narrow hydrographs while elongated catchments give narrow peaked and broad hydrographs. Size of a basin affects its behaviour. In small catchments channel flows are predominant by overland flows. Hence intensity of rainfall and land use has greater impact on peak flood. In larger catchments channel flow is more predominant. Time base for larger catchments will be larger as compare to corresponding hydrographs for smaller catchments. Slope of the stream controls the velocity of flow in channel. Slope has significant effect on recession limb of hydrograph as it represents the depletion of storage. The basin slope for smaller catchments is more important where overland flow is predominant. Steeper slopes result in larger peak discharges. Ratio of total channel length to total drainage area is known as drainage density. Large drainage density results in quick conduction of water down the channel that reflects as pronounced peak discharge. In basins with small drainage densities, the overland flow is more predominant and hence the resulting hydrograph is squat with a slowly rising limb. Vegetal cover offers resistance to the flow of water and increases the infiltration and storage capacities of soil. This results in reduced peak of hydrograph. This effect is very pronounced for catchments with smaller areas. Intensity, duration and direction of storm movement are the most important climatic factors that affect the shape of a flood hydrograph. Peak and volume of surface runoff are directly proportional to intensity of rainfall for a given duration. Duration of a storm of given intensity also has direct proportional effect on volume of runoff. There will be quicker concentration of flow if a storm moves from upstream end of catchment to downstream end. This results in high peak flood. This effect is further attenuated by shape of catchments. Hydrographs for long and narrow catchments are most sensitive to direction of storm movements. Detailed information about rainfall and resulting flood hydrographs is not available for all catchments, especially in remote and small catchments. To construct unit hydrographs for such catchments empirical relationships of regional validity has been given by a number of researchers.8 A number of methods are reported in literature for developing synthetic hydrographs. But these methods being based on empirical relationships have regional limitations and they should not be considered as general relationships for all regions. Central Water Commission (CWC) has classified India in seven hydro- meteorologically homogenous zones. Each zone is described by a different set of relationships for the synthetic unit hydrograph. For example, western Himalayan region lies in the zone 7. However, the recommended synthetic unit hydrograph is only useful for catchment area greater than 25km².9 Therefore, it cannot be used for derivation of unit hydrograph for cloudburst events as it is a localized and restricted to area less than 25km2. Snyder (1938) developed a set of empirical equations for synthetic unit hydrograph based on the study for a large number of catchments in USA. These equations are as follow:
(4)
(5)
(6)
(7)
(8)
(9)
(10)
Where(11)
(12)
WhereValue of Ct in Snyder’s study ranged from 1.35 to 1.65 but many other studies have shown that its value depends on the region under study. Wide variation in the value of Ct (0.3- 6.0) has been reported. Value of Cp is considered to be a representation of retention and storage capacity of catchment and ranges from 0.56 to 0.69 for Snyder’s study. But depending on the region under study these coefficient may have values outside these ranges also. Linsley et al.10 have suggested modified formula for basin lag as:
(13)
Where,
S: Basin slope; Ctl, Basin constant; n: Regional constant (0.38 for study of catchments in USA).
US Soil Conservation services (SCS) has suggested equation to estimate the value of peak discharge Qp and time of peak Tp which is very popular method used in watershed development activities especially in small catchments. The equations that are developed based on the study of large no of catchments are as follow:
(14)
(15)
(16)
(17)
Where
: Time of concentration; : Duration of effective rainfall; : Lag time; : Time of peak; : Peak discharge in cumecs; A: Area of catchment in km2
SCS dimensionless unit hydrograph facilitates the construction of unit hydrograph. In this ordinate is Q/Qp that is the discharge Q expressed as the ratio to peak discharge Qp and the abscissa is t/tp which is the time t expressed as the ratio of time to peak tp (Table 1). This unit hydrograph has same percentage of volume on rising side as the SCS triangular unit hydrograph. Since most of the Himalayan catchments facing cloudbursts are very remote and of smaller areal extent, they are generally unguaged. There is no single recommended method for flood estimation of such events. Therefore, in this study, applicability of these methods is examined for the assessment of flood hydrograph and its characteristics by using empirical, Rational, Snyder’s Synthetic Unit Hydrograph and SCS method. Various coefficients that represent the catchment characteristics in these methods are based on normal conditions before the occurrence of a storm. However, in case of cloudburst, generally the soil is fully saturated and depressions are already full. This causes a significant rise in the peak of outflow hydrograph and decreases the time of peak. Also, a major part of precipitation appears as surface runoff because the initial abstractions are very less. So the coefficients for runoff are higher as compare to those in normal cases. Further, considering steep slopes, small size of catchment, lack of vegetation, very small initial abstraction and saturation of soil make it imperative to select different coefficients. As a case study, cloud burst event which occurred in Leh (J & K State in north Himalaya) during night of 5th August 2010. Gridded rainfall data of India Meteorological Department (IMD) suggested the storm intensity of 120mm/day on 5th August, 2010, whereas Kumar et al.11 suggest occurrence of 70 mm precipitation in 3h on 5th August. Area of catchment is about 0.842 km2 and length of main channel till the outlet of the catchment area is about 1.25km. This event is analysed by Renoj et al.4 This event is analysed by various methods in this study and is discussed as follows: Value of Dicken’s coefficient for northern India is taken as 11.5. Considering this value, the peak flood for Leh catchment for unit rainfall is estimated to be 10.11m3/s. Attempt has been made to use the different values of runoff coefficient, slope and time of peak and the flood is estimated using Rational method. Further, in India, it is a recommended common practice to increase the peak flood by 25% for cloudburst event. The results are shown in Table 2. Characteristics of unit hydrograph using Snyder’s SUH approach have been derived for different sets of with ranging from 0.15 to 0.30 and ranging from 0.75 to 0.95 as shown in Table 3. SCS triangular Unit hydrograph method is used to estimate peak flood and peak time for the catchments and then the hydrograph (Figure 5) is constructed by SCS dimensionless unit hydrograph method using the peak food and base time. Same equation has been used for time of peak (tp) as given by Snyder’s method. Peak flow is estimated to be 7.429m3s/and time of peak as 0.236 hrs. Base time for the hydrograph is estimated to be 1.28 hrs. The results obtained from above three methods are checked for the compliance with requirements of unit hydrograph and then compared with the results reported by Renoj et al.4 estimated by using atmospheric and hydrological modelling approach.
t/Tp |
q/Qp |
t/Tp |
q/Qp |
0 |
0 |
1.4 |
0.75 |
0.1 |
0.02 |
1.5 |
0.66 |
0.2 |
0.08 |
1.6 |
0.56 |
0.3 |
0.16 |
1.8 |
0.42 |
0.4 |
0.28 |
2 |
0.32 |
0.5 |
0.43 |
2.2 |
0.24 |
0.6 |
0.6 |
2.4 |
0.18 |
0.7 |
0.77 |
2.6 |
0.13 |
0.8 |
0.89 |
2.8 |
0.1 |
0.9 |
0.97 |
3 |
0.07 |
1 |
1 |
3.5 |
0.04 |
1.1 |
0.98 |
4 |
0.02 |
1.2 |
0.92 |
4.5 |
0.01 |
1.3 |
0.84 |
5 |
0 |
Table 1 Coordinates of SCS dimensionless unit hydrograph
Slope |
Runoff |
tp |
Qp |
|
|
hr |
m³/s |
0.1 |
0.85 |
0.25 |
4.98 |
0.1 |
0.9 |
0.25 |
5.26 |
0.1 |
0.95 |
0.25 |
5.55 |
0.07 |
0.85 |
0.3 |
4.96 |
0.07 |
0.9 |
0.3 |
5.26 |
0.07 |
0.95 |
0.3 |
5.55 |
0.05 |
0.85 |
0.33 |
4.98 |
0.05 |
0.9 |
0.33 |
5.26 |
0.05 |
0.95 |
0.33 |
5.55 |
Table 2 Peak flood and time of peak for various values of runoff coefficient and slope estimated using rational method
Ct |
Cp |
Qp(m3/s) |
tp hr |
Tb hr |
0.3 |
0.75 |
5.49 |
0.32 |
2.09 |
0.3 |
0.8 |
5.86 |
0.32 |
2.09 |
0.3 |
0.85 |
6.23 |
0.32 |
2.09 |
0.3 |
0.9 |
6.59 |
0.32 |
2.09 |
0.3 |
0.95 |
6.96 |
0.32 |
2.09 |
0.25 |
0.75 |
6.4 |
0.28 |
1.87 |
0.25 |
0.8 |
6.82 |
0.28 |
1.87 |
0.25 |
0.85 |
7.23 |
0.28 |
1.87 |
0.25 |
0.9 |
7.67 |
0.28 |
1.87 |
0.25 |
0.95 |
8.1 |
0.28 |
1.87 |
0.2 |
0.75 |
7.65 |
0.23 |
1.64 |
0.2 |
0.8 |
8.16 |
0.23 |
1.64 |
0.2 |
0.85 |
8.67 |
0.23 |
1.64 |
0.2 |
0.9 |
9.18 |
0.23 |
1.64 |
0.2 |
0.95 |
9.68 |
0.23 |
1.64 |
0.15 |
0.75 |
9.51 |
0.19 |
1.42 |
0.15 |
0.8 |
10.14 |
0.19 |
1.42 |
0.15 |
0.85 |
10.78 |
0.19 |
1.42 |
0.15 |
0.9 |
11.41 |
0.19 |
1.42 |
0.15 |
0.95 |
12.04 |
0.19 |
1.42 |
Table 3 Unit hydrograph characteristics for various values of Ct and Cp
Unit Hydrographs characteristics for different values of coefficients have been checked for compliance with basic requirements of unit hydrograph theory (Figure 3). Unit hydrograph obtained using SCS method and Snyder’s method (Ct =0.25, Cp=.85) have almost equal peak and also have same volume of runoff. Rational formula gives small peak as compare to other methods as well as estimated by Renoj et al.4 Therefore, it is recommended that flood hydrograph obtained by Snyder’s method with Ct =0.25; and Cp=0.85(Figure 4). However, one can also use the SCS method (Figure 5) as it is less data intensive with reasonable level of accuracy. The results of Snyder’s (Ct =0.25, Cp=.85) and SCS methods are compared with Renoj et al.4 (Table 4) and found to hold good considering the uncertainties involved in estimation12,13 Based on the identified parameters of the Snyder’s and SCS approach, unit hydrograph characteristics for different size catchments have been assessed for 15 and 30 minutes storm durations as shown in Table 5.
Method |
Dicken’s |
Snyder |
SCS |
SCS* |
tp(hr) |
- |
0.258 |
0.236 |
- |
tb(hr) |
- |
1.87 |
1.28 |
0.92 |
Qp(m3/s) |
168.8 |
120.9 |
124.06 |
107 |
Method |
Duration |
Snyder |
Scs |
|
|
area (km²) |
tr(hr) |
tp' (hr) |
Qp(m3/s ) |
tp |
Qp(m3/s ) |
1 |
0.5 |
0.34 |
6.82 |
0.23 |
8.96 |
1 |
0.25 |
0.28 |
8.32 |
0.23 |
8.96 |
5 |
0.5 |
0.49 |
23.67 |
0.39 |
26.52 |
5 |
0.25 |
0.43 |
27.05 |
0.39 |
26.52 |
10 |
0.5 |
0.65 |
36.25 |
0.55 |
37.69 |
10 |
0.25 |
0.58 |
40.1 |
0.55 |
37.69 |
Table 5 SUH characteristics by Snyder and SCS methods
It is concluded that cloud burst contributes significantly in flood disasters particularly in hilly areas. Rational formula gives the underestimated values of peak flood whereas Dicken’s method overestimates the same and both these methods provide peak flood value only. On the other hand flood estimation using techniques of unit-hydrograph for an ungauged catchment can help in estimating the temporal variation of cloud burst induced flood in hilly isolated areas. Further, results show that both SCS Triangular Unit Hydrograph method and Snyder’s SUH method (for Ct= 0.25 and Cp=0.85) can be used for peak flood estimation for cloudburst event in unguaged Himalayan catchments. However since Snyder’s method takes into account more number of parameters and hence give more reliable results. So Snyder’s method is suggested for estimation of flood hydrograph characteristics using corrected values of coefficients (Ct= 0.25 and Cp=0.85). SCS method may be used for rough estimations as it is easy and needs no adjustments. Further, there is need for comprehensive data base of such events over the entire Himalayan range for research work and hence effective disaster management in the region.14,15
None.
Authors declare there is no conflict of interest in publishing the article.
©2017 Kansal, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.