REGIONAL RAINFALL INTENSITY - DURATION - FREQUENCY RELATIONSHIPS

Transcription

1

2

3 Module-3 REGIONAL RAINFALL INENSIY - DURAION - FREQUENCY RELAIONSHIPS Dr. N. K. Goel Professor & Head Department of Hydrology Indian institute of echnology, Roorkee & Shibayan Sarkar Research Scholar Department of Hydrology Indian institute of echnology, Roorkee INRODUCION Rainfall intensity of a particular frequency and duration is required for estimation of average annual flood and sediment yield from the catchments. he design engineers do not have simple and reliable method for estimation of rainfall intensity, particularly for short duration. his necessitates going for the regionalisation of IDF relationships. he limitation of such study is the scarcity of data from self-recording raingauges for sufficiently long time. Again, studies for less than 1-hour data are also scarce. However, an effort has been made to briefly put forth the IDF relationships for various regions (mostly in Indian context) in the forthcoming sections. REVIEW OF LIERAURE Sherman in 1931 developed an empirical relationship of the form, a K I = (1) b ( t+ c) Where, t is duration in minutes, is return period, K, a, b, and c are constants depending on geographical location. his is the most common form of IDF relationship, which is still being used widely. Bernard (1932) developed an empirical relationship in the form of: I t I t ao = a 2 t a1 Where, = rainfall intensity having duration 't' and return period '', ao,a 1 & a 2 are constants and depend upon geographical location. Bilham (1935) published his well - known article on the IDF relationship for U.K. and the frequencies were calculated from the formula as (converted in S.I. Units). (2) Organized by: Department of Hydrology, II Roorkee, Roorkee , India 46

4 n = * 10 5 t [ R ] (3) Where, R is depth of rainfall in mm, n is no. of occurrences in 10 years and t is duration of rain in hours. Bilham's work was modified by Holland (1967), who showed that Bilham's equation over estimates the probabilities of high intensity rainfall (> 35 mm/hr). Yarnall (1935) presented such data in the form of maps of a region with isohyetal lines indicating total rainfall depth that may be expected in a time 't', at a frequency once in N-years for United States. In Australia, the Bureau of Meteorology (Hall, 1984) had developed procedures for estimation of IDF values for return period up to 100 years and duration up to 72 hours. General extreme value (GEV) distribution has been widely used for application to rainfall extremes in Australia, East Africa and U.K. his also forms the essence of rainfall studies in volume - II of U.K. Flood studies report (National Environmental Research Council, 1975). Gumbel distribution (a specific case of GEV distribution) has also been used in several reports of U.S. weather Bureau for durations ranging from 2 minutes to 24 hours. Bell (1969) proposed the following the following depth - duration - frequency formula: Where, [ 0.21ln ][ 0.54t ] R R t = 1 (4) For 2 (years) 100 & 5 t (min.) 100 R t 10 R 1 = -year and t-hour rainfall depth in inches = 10-year and 1-hour rainfall depth in inches Baghirathan and Shaw (1978) made rainfall depth-duration-frequency studies for Sri Lanka. Chen (1983) provided a general relationship for rainfall intensity in U.S.A. Raudkivi of New Zealand presented regional relationship on IDF in All these authors used Bell's equation in their studies. Ferreri and Ferro (1990) verified the applicability of Bell's equation for Sicily and Sardinia in the Mediterranean. Neimczynowicz (1982) used Log Pearson ype-iii distribution with method of moments for preparing areal IDF curves for short-term rainfall events in Lund, Sweden. Steel & McGhee (1979) gave the empirical relationship for United States for duration less than 2 hours and for any given frequency as: A I = (unbalanced) (5) t + B Where, I is intensity in inches/hour, t is duration in min and A & B are constants depending on frequency and climatic condition. Rao et al. (1983) obtained relationship between short duration rainfall and 24 hour rainfall as: 2 I (t) = a + b R 24 + C R 24 (6) Where, a,b,& c are constants. Gert et al., (1987) obtained the following relationships for Pennsylvania, U.S.A. I (t) = ( log t 24) R 24 (7) Where, I(t) is rainfall amount for a duration of 't' hours and R 24 is 24 hours rainfall amount. Organized by: Department of Hydrology, II Roorkee, Roorkee , India 47

5 Mc Cuen (1989) gave the mathematical representation of IDF curves for computerizing the elements in hydrologic design for Baltimore in Maryland and the equations are: I = A for t 2 hours t+ b I = c t d for t > 2 hours (8) Where I is intensity in inches/hour, t = duration in hours and a, b, c & d are coefficients that vary with frequency. Chen's (1983) represented IDF as: Where, I t I t = 10 2-x x-1 a I 1 log (10 ) c ( t + b) is Rainfall intensity in inches/hour for -year and t-min storm duration, x is depth of frequency ratio ( Rt / Rt ), a, b & c = storm parameters that are dependent on regional ratio (Rt 100 / R 10 t ) Chen showed that the 10-year, 1 hr rainfall R 10 1 used in Bell's equation alone can not measure the geographical variations of rainfall and equation (8) produces more accurate results. A comparison of available IDF relationships for short durations, given by Chow (1964), Raudkivi (1979), Gert et al. (1987), and Chen (1983) reveals that the values of exponents of variables 't' & '' in Bernard's equation do not vary much from place to place for shorter duration rainfalls. he exponent of '' ranges between 0.18 & For 't' the exponent varies from 0.7 to Koutsoyiannis et al. (1998) used several appropriate statistical distributions functions ( Gumbel, Gamma, GEV, Log Pearson III Log Normal Exponential, Pareto etc.) for a reliable parameter estimation of IDF relationships and proposed a generalized empirical IDF relationships Equation (10) facilitating the description of the geographical variability with regionalization of IDF curves. Moreover, this formulae allows incorporating data from non-recording stations. ω i = (10) ( d ν + θ ) η Where,ω,ν,θ and η are non-negative coefficients withνη 1. he inequality is easily derived from the demand that the rainfall depth h = id is an increasing function of d. Based on the partial duration series (PDS) approach Madsen et al. (1998) proposed regional estimation of extreme precipitation from a high resolution rain gauge network in Denmark. For a preliminary assessment of regional homogeneity and identification of a proper regional distribution L-moment analysis is applied. o analyze the regional variability in more detail, a generalized least squares regression analysis is carried out that relates the PDS model parameters to climatic and physiographic characteristics. A regional parent distribution is identified as the Generalized Pareto distribution. Sivapalan and Blöschl (1998) presented an methodology for IDF curve based on the spatial correlation structure of rainfall linking scientific theories of space-time rainfall fields with design methods rather using by the use of empirically-derived areal reduction factors (ARFs). (9) Organized by: Department of Hydrology, II Roorkee, Roorkee , India 48

6 Yu and Cheng (1998) formulated a generalized regional IDF relationship by pooling annual maximum rainfall series of southern aiwan. Naghettini (2000) investigated the properties of time scale invariance of rainfall applied to intensity-duration-frequency relationships for short-duration rainfall for southeastern Brazil from the statistical characteristics of daily data. he hypothesis of simple scaling implies in direct and empirically-verifiable relations among the moments of several orders and also among the probability distributions of rainfall intensities of different durations. refry et al. (2000) applied method of moments or maximum likelihood procedures to fit a suitable probability distribution to annual maximum or partial duration series data for each gage of Michigan to estimates site-specific IDF curve. Further, isopluvial maps have been developed from these using interpolation procedures and judgment. Davis and Naghettini (2000) aimed to estimate regional intensity-durationfrequency (IDF) curves using partial duration series and for the Brazilian state of Rio de Janeiro, whereas, application of L-Moments have improved parameter and quantile estimates of extreme rainfall intensities for this study. refry and Watkins Jr. (2001) discusse the application of a PDS / Generalized Pareto (GPA) regional index-flood procedure assuming independent peaks and a Poissondistributed number of threshold exceedances for updating rainfall intensity-durationfrequency (IDF) estimates for the State of Michigan. Yu et al. (2004) developed regional Intensity Duration Frequency (IDF) formulas for non-recording sites of northern aiwan based on the scaling theory using annual maximum rainfall series for various durations. he hypothesis of piecewise simple scaling combined with Gumbel distribution was used to develop the IDF scaling formulas. Durrans and Kirby (2004) studied IDF curves and 24-h design storm hyetographs based on US Soil Conservation Service which consist of at-site estimates of generalized extreme value distribution parameters for gauging sites of Alabama. Maurino (2004) compared the generalized rainfall intensity-duration-frequency relationships proposed by Bell in 1969 with results obtained from data registered in different climatic regions of Argentina. Amin and Shaaban (2004) used GEV and EV1 distribution with one step least square method for parameter estimation to propose regionalized IDF relationship for Peninsular Malaysia. Hadadin (2005) constructed IDF curve for Mujib Basin in Jordon which is compared with Gumbel method and Water Authority of Jordon. refry et. Al. (2005) applied a regional frequency analysis approach based on L moments to carry out storm water management plans for the State of Michigan; considering the area as a homogeneous region two regional index flood models were applied: a generalized Pareto distribution fit to partial duration series data (PDS) data, and a generalized extreme value distribution fit to annual maximum series (AMS) data. Ghahraman and Hosseini (2005) carried out a study to test the performance of commonly used IDF models for three synoptic stations in Iran. Guo (2006) extended IDF analysis from previous analysis results for Chicago, considering the effect of changing climate conditions of the first and second halves of the last century. Further this IDF is applied on the design and performance of urban drainage systems. Organized by: Department of Hydrology, II Roorkee, Roorkee , India 49

7 Nhat et al. (2006) constructed IDF curves for the monsoon area of Vietnam and to propose a generalized IDF formula using base rainfall depth, and base return period for Red River Delta (RRD) of Vietnam. Desa et al. (2006) proposed an approach to estimate design rainfall depth at low ARI using the Partial Duration (PD) data series of short duration high quality rainfall data to overcome the drawback of At site frequency analysis for economical and efficient urban drainage design of Sungai. Finally, they showed relatively simple application which is relevant for the at-site and regional IDF relationship. Bougadis and Adamowski (2006) examined scaling properties of extreme rainfall to establish scaling behaviour of statistical non-central moments over different durations, whereas, a scale invariance concept is explored for disaggregation (or downscaling) of rainfall intensity from low to high resolution and is applied to the derivation of scaling IDF curves based on scaling of the generalized extreme value (GEV) and Gumbel probability distributions for the province of Ontario. Regalado and Yuste (2006) proposed an intra-station regionalization, meaning, a regionalization in the same station for Spain. In addition to this, they incorporate GIS (Geographical Information System) application to improve this relationship called MAXIN. For small and/or urban catchments of Scotland, Svensson et al. (2007) compared a number of IDF relationship approaches calculated from each artificially fragmented record which is modified considering the existing data gaps of varying degrees of severity. Raiford et al. (2007) updated the existing intensity-duration-frequency curves for ungauged sites of South Carolina, North Carolina, and Georgia. he L-moment method with X-10 test was used to search for homogeneous regions within the study area. Finally, at-site statistics were calculated to develop frequency relationships. Normal, lognormal, generalized extreme value, Pearson type III, and log Pearson type III probability distribution functions were used to fit the maximum annual precipitation data at each gauging site for each duration. he chi-squared goodness-of-fit test was used to determine the best fit probability distribution. Salas et al. (2007) carried out IDF studies for Spain incorporating regionalization, and more adequate functions, such as the SQR-E max function of distribution. REGIONAL IDF RELAIONSHIPS FOR INDIAN REGIONS Parthsarathy and Singh (1961) prepared rainfall intensity duration frequency curves for India for local drainage design. Based on 15 min tabulations of rainfall for 50 stations and hourly tabulations for 67 self-recording raingauges in India, Ayyar and ripathi (1973,1974) have prepared generalized charts of 2,5,10,25 & 50 year return period values of 15, 30, 45 min, 3, 6, 9, 12 & 15 hours rainfall. Ram Babu et al. (1979), after analysing rainfall characteristics for 42 stations, presented IDF equations and nomographs. With the equations and/ or nomographs, the intensity for any desired duration for a given frequency may be determined. Its general form is: a K I = (11) ( t+ b) n Where, I is intensity in cm/hr, is return period in years, t is storm duration in hours. K, a, b & n are constants developed for various stations and zones of India. he values of K, a, b and n are presented in able 1 for different stations and zones of India. Organized by: Department of Hydrology, II Roorkee, Roorkee , India 50

8 able 1 Intensity duration return period relationship, India Zone Station K a b n Northern Agra zone Allahabad Amristar Dehradun Jaipur Jodhpur Lucknow New Delhi Srinagar Northern Zone Central Bagra-tawa zone Bhopal Indore Jabalpur Jagdalpur Nagpur Punase Raipur hikrl Central zone Western Aurangabad zone Bhuj Mahabaleshwar Nandurbar Vengurla Veraval Western Zone Eastern Agarthala zone Dumdum Gauhati Gaya Imphal Jamshedpur Jharsuguda North Lakhimpur Sagarisland Shillong Eastern Zone Southern Bangalore zone Hyderabad Kodaikanal Madras Mangalore iruchirapalli rivandrum Visakhapatnam Southern Zone Rambabu et al. (1979) also give monograph (Fig 1) to convert one hour rainfall intensity into rainfall intensities of other durations. For Vasad and Kota, the following relationships have been obtained by Central Soil Water Conservation Research and raining Institute. Organized by: Department of Hydrology, II Roorkee, Roorkee , India 51

9 Figure 1 Locations of raingauge stations and zonal boundaries in India Vasad I = (12) (t + 0.5) Kota I = (13) (t + 0.5) 0.85 Kothyari & Garde (1992) developed a general relationship on IDF after analysing 80 raingauge stations (Fig. 2) in India. hey have made use of the assumption that general properties of convective cells that are associated with short - period rainfalls are similar in different hydrologic regions Raudkivi (1979). he formula is of the form: I t Where, I t = C t (R 2 24 ) 0.33 (14) = rainfall intensity in mm/hr for -year return period & t hours duration Organized by: Department of Hydrology, II Roorkee, Roorkee , India 52

10 C = constant having value 8.31 for the whole of the considered stations (Values for different regions are given in able 2) 2 R 24 = 2-year return period & 24 hr. rainfall in mm. Figure 2 Relation to convert one our rainfall to intensities for other durations able 2 Values of 'C' for different regions of India Geographical region Zones in figure 1 Value of C Northern India Eastern India Central India Western India Southern India In India, number of studies have been carried out to analyse the data of individual stations e.g. Rama and Krishna (1958) for Delhi, Alipore and Madras Stations, Rama and Bandyopadhyay (1969) for Calcutta. CONCLUSIONS Gumbel s extreme value distribution has been most oftenly used for analysis of the short duration rainfall data. he IDF relationships developed for Indian regions are based on limited data of very few stations and there is need to further improve upon the relationships based on larger data base. Organized by: Department of Hydrology, II Roorkee, Roorkee , India 53

11 REFERENCES Amin, M. Z. M., Shaaban, A. J. (2004). he rainfall intensity-duration-frequency (IDF) relationship for Ungauged sites in peninsular Malaysia using a Mathematical Formulation. 1st international conference on managing rivers in the 21st century: issue and challenges, Rivers 04 Ayyar, P. S. H & ripathi, N. (1973) - One day rainfall of India for different return periods, No.5, IMD, N.Delhi. Ayyar, P. S. H. & ripathi, N. (1974) - Rainfall frequency maps of India, No.-6, IMD, N. Delhi. Baghirathan, V. R. & Shaw, E. M. (1978)- Rainfall depth - duration- frequency studies for SriLanka, Journal of Hydrology, 37(3), Bell, F. C. (1969) - Generalized rainfall-duration-frequency relationships, Journal of Hydraulics Engineering, ASAE, 95(1), Bernard, M. M. (1932) - Formulas for rainfall intensities of long duration, rans. ASCE, Vol. 96, pp Bilham, E. G. (1935) - he classification of heavy falls of rain in short periods, H.M.S.O., London, 1962 (republished). Bougadis, J., Adamowski, K. (2006). Scaling model of a rainfall intensity-durationfrequency relationship, Hydrol. Process. 20, Chen, C. L. (1983) - Rainfall intensity - duration - frequency formula, J. of Hydraulics Engineering, ASCE, 109(12), Chow, V.. (Ed.) (1964) - Handbook of applied hydrology., Mc Graw Hill, New York. Davis, E. G., Naghettini, M. (2000). Regional Analysis of Intensity-Duration-Frequency of Heavy Storms over the Brazilian State of Rio de Janeiro, Journal of Hydrologic Engineering,138 Desa, M. N. M., Amin, M. Z. M., Asnor, M. I., Mohamed, F. (2006). Developing a Low Return Period Regional IDF Relationship using Generalized Pareto (GPA) Distribution, National Conference Water for Sustainable Development owards a Developed Nation by 2020, July, Guoman Resort Port Dickson Durrans, S. R., Kirby, J.. (2004). Regionalization of extreme precipitation estimates for the Alabama rainfall atlas, Journal of Hydrology, 295, Ferreri, G. B. & Ferro, V. (1990) - Short duration rainfalls in Sicily, Journal of Hydraulics Engineering, ASCE, 116(3), Gert, A., Wall, D. J., White, E.L. & Dunn, C.N. (1987) - Regional rainfall intensity - duration-frequency curves for Pennsylvania, Water Res. Bulletin, 23(3), Ghahraman, B., Hosseini, S. M. (2005). A new investigation on the performance of rainfall IDF models, Iranian Journal of Science & echnology, ransaction B, Engineering, 29 (B3), Guo, Y., (2006). Updating Rainfall IDF Relationships to Maintain Urban Drainage Design Standards, Journal of Hydrologic Engineering, 11(5), Organized by: Department of Hydrology, II Roorkee, Roorkee , India 54

12 Hadadin, A. N. (2005). Rainfall Intensity Duration-Frequency Relationship in the Mujib Basin in Jordan. Journal of Applied sciences, 8(10), Hall, A. J. (1984) - Hydrology in tropical Australia & Papua New Guinea, Hydrological Sciences Bulletin, Vol.-29(4). Holland, D. J. (1967) - Rainfall intensity - frequency relationships in Britain, British Rainfall 1961, H.M.S.O. Kothyari, U.C., and Garde, R. J. (1992), - Rainfall intensity - duration-frequency formula for India, Journal of Hydraulics Engineering, ASCE, 118(2). Koutsoyiannis, D., Kozonis, D., Manetas, A. (1998). A mathematical framework for studying rainfall intensity-duration-frequency relationships.journal of Hydrology, 206, Madsen, H., Mikkelsen, P.S., Rosbjerg, D. and Harremoës, P. (1998). Estimation of regional intensity-duration-frequency curves for extreme precipitation. Water Science and echnology, 37(11), Maurino, M. F. (2004). Generalized Rainfall-Duration-Frequency Relationships: Applicability in Different Climatic Regions of Argentina.Journal of Hydrologic Engineering, 9(4), Mc. Cuen, R. H. (1989) - Hydrological analysis and design, Englewood Cliffs, New Jersey,: Prentice Hall Naghettini, M. (2000). A Study of the Properties of Scale Invariance as Applied to Intensity-Duration-Frequency Relationships of Heavy Storms. Journal of Hydrologic Engineering. 139 National Environmental Research Council (1975)- Flood studies report, Vol.-1, Hydrologic studies, NERC, London. Neimczynowicz, J. (1982) - Areal intensity-duration-frequency curves for short term rainfall events, Nordic Hydrology, Vol.13(4), Nhat, Le M., achikawa,y., akara, K. (2006). Establishment of Intensity-Duration- Frequency Curves for Precipitation in the Monsoon Area of Vietnam. Annuals of Disas. Prev. Res. Inst., Kyoto Univ., No. 49 B, Parthasarathy, K. & Singh G., (1961) - Rainfall intensity-duration-frequency for India for local drainage design - Indian Journal of Meteorology & Geophysics, Vol.-12(2), Raiford, J. P., Aziz, N. M., Khanand, A. A., Powell, D. N. (2007). Rainfall Depth-Duration- Frequency Relationships for South Carolina, North Carolina, and Georgia, American Journal of Environmental Sciences, 3 (2), Ram Babu, ejwani, K. K., Agrawal, M. C. & Bhusan, L. S. (1979) - Rainfall intensityduration-return period equations & nomographs of India, CSWCRI, ICAR, Dehradun, India. Raman, P. K. & Krishnan, A. (1958) - Intensity-duration-frequecy relationships of rainfall at selected stations in India, Proceedings Symp. on Meteorological and Hydrological aspects of floods & droughts in India, Raman, V. & Bandyopadhya, M. (1969) - Frequency analysis of rainfall intensities for Calcutta, Journal of Sanitary Engineering Dn., ASCE, 95(6), Organized by: Department of Hydrology, II Roorkee, Roorkee , India 55

13 Rao, D. V. L. N., Goyal, S. C. & Kathuria, S. N. (1983) - Interpolation of short duratin rainfall from 24 hour rainfall in Lower Godavari basin, Mausam Journal, Meteorological Soc., N. Delhi, 34(3): Raudkivi, A. J. (1979) - An advance introduction to Hydrological Processes and Modelling, Oxford Pergamon Press, N. York, N.Y. Regalado, L. de S., Yuste, J. A. F. (2006). Maximum rainfall intensity analysis using l- moments in Spain, Proceedings of the 7 International Conference on HydroScience and Engineering Philadelphia, USA September 10-13, 2006, Aug, 2007 Salas, L. de, Carrero, L., Fernández, J.A. (2007). MAXIN:Gis application to estimaterainfall Intensity-Duration-Frequency laws in the Spanish peninsular area. MS7/ECAM8 Abstracts, Vol. 4, 7th EMSAnnual Meeting / 8th ECAM,EMS2007-A Sivapalan, M., Blöschl, G. (1998). ransformation of point rainfall to areal rainfall: Intensity-duration-frequency curves. Journal of Hydrology, 204, Steel, E. W. & McGhee,. J. (1979) - Water Supply and Sewerage, 5th ed., McGraw-Hill Book Com., New York. Svensson, C., Clarke, R.., Jones, D. A. (2007). An experimental comparison of methods for estimating rainfall intensity-duration-frequency relations from fragmentary records. Journal of Hydrology, 341, refry C. M., Watkins D. W. Jr. (2001). Application of a Partial Duration Series Model for Regional Rainfall Frequency Analysis in Michigan, Journal of Hydrologic Engineering. 41 refry, C. M., Watkins, D. Jr., Johnson, D. (2005). Regional Rainfall Frequency Analysis for the State of Michigan, Journal of Hydrologic Engineering, 10(6), refry, C. M., Watkins, D. W., Johnson, D. L. (2000). Development of Regional Rainfall Intensity-Duration-Frequency Estimates for the State of Michigan, Journal of Hydrologic Engineering, 140 Yarnall, D. L. (1935) - Rainfall intensity frequency data, USDA, Misc. Publ. No. 204, Washington, D. C. Yu, P. S., Lin, C. S., Yang,. C. (2004). Regional rainfall intensity formulas based on scaling property of rainfall. Journal of Hydrology, 295, Organized by: Department of Hydrology, II Roorkee, Roorkee , India 56