Print version ISSN 0100-0683
Rev. Bras. Ciênc. Solo vol.36 no.3 Viçosa May/June 2012
Duration-frequency relationships of heavy rainfall in Santa Catarina, Brazil1
Relações duração-frequência de chuvas intensas no estado de Santa Catarina
Álvaro José BackI; José Luiz Rocha OliveiraII; Alan HennIII
IEpagri Researcher, Experimental Station of Urussanga, Rod.SC 446, km 16, Zip Code 88840-000, Urussanga (SC), Brazil. E-mail: firstname.lastname@example.org
IIEpagri Researcher, Information Center of Hydrometeorology and Environmental Resources of Santa Catarina/Ciram, Rod Admar Gonzaga 1347, Itacorubi, PO Box 502, Zip Code 88034-901, Florianópolis (SC), Brazil, E-mail: email@example.com
IIIEpagri Researcher, Information Center of Hydrometeorology and Environmental Resources of Santa Catarina/Ciram, Rod Admar Gonzaga 1347, Itacorubi, PO Box 502, Zip Code 88034-901, Florianópolis (SC), Brazil, E-mail: firstname.lastname@example.org
The purpose of this study was to adjust equations that establish relationships between rainfall events with different duration and data from weather stations in the state of Santa Catarina, Brazil. In this study, the relationships between different duration heavy rainfalls from 13 weather stations of Santa Catarina were analyzed. From series of maximum annual rainfalls, and using the Gumbel-Chow distribution, the maximum rainfall for durations between 5 min and 24 h were estimated considering return periods from 2 to 100 years. The data fit to the Gumbel-Chow model was verified by the Kolmogorov-Smirnov test at 5 % significance. The coefficients of Bell's equation were adjusted to estimate the relationship between rainfall duration t (min) and the return period T (y) in relation to the maximum rainfall with a duration of 1 hour and a 10 year return period. Likewise, the coefficients of Bell's equation were adjusted based on the maximum rainfall with a duration of 1 day and a 10 year return period. The results showed that these relationships are viable to estimate short-duration rainfall events at locations where there are no rainfall records.
Index terms: design rainfall, hydrology, drainage.
O objetivo deste trabalho foi ajustar equações que relacionam as precipitações com diferentes durações aos dados medidos em estações pluviográficas do Estado de Santa Catarina. Foram estudadas as relações entre precipitações intensas de diferentes durações de 13 estações pluviográficas do Estado de Santa Catarina. A partir das séries de chuvas máximas anuais, foram estimadas, por meio da distribuição de Gumbel-Chow, as chuvas máximas com durações variando de 5 min a 24 h, considerando-se os períodos de retorno de 2 a 100 anos. A aderência dos dados ao modelo de Gumbel-Chow foi comprovada por meio do teste de Kolmogorov-Smirnov, com nível de significância de 5 %. Ajustaram-se os coeficientes da equação proposta por Bell a fim de estimar a relação entre a chuva com duração t (min) e o período de retorno T (anos), em função da chuva máxima com período de retorno de 10 anos e duração de 1 h. Do mesmo modo, foram ajustados os coeficientes da equação de Bell, cuja base foi a chuva máxima com período de retorno de 10 anos e duração de um dia. Os valores obtidos mostram que o uso dessas relações constitui-se em alternativa viável para estimar as chuvas de curta duração em locais sem registros pluviográficos.
Termos de indexação: chuva de projeto, hidrologia, drenagem.
Knowledge on rainfall characteristics allows a safe planning and determination of the dimensions of soil conservation structures such as terraces and contour farming. Likewise, the dimensions of hydraulic structures to retain or channel surface runoff after heavy rain, e.g., drainage channels, gullies, retention basins and others are determined based on the maximum rainfall intensity with a certain probability of being equaled or exceeded and defined by the return period.
The dimensions of surface drainage works are determined by adopting values of rainfall events that occur with return periods between 5 and 100 years to design rainfall. To obtain these rains, the annual maximum series observed in more than 10 years is generally used and theoretical probability distributions are adjusted to estimate the maximum rainfall with a given return period or probability of being equaled or exceeded. For heavy rainfall events, several studies showed that the Gumbel distribution adjusts well and has therefore been widely used, although a series of annual maximum values must be available for its application (Kite, 1978; Back, 2001).
In drainage projects, it is often necessary to obtain the intensity of short-duration rainfall, usually shorter than one day. Typically, the rainfall intensity is estimated based on Intensity-Duration-Frequency (IDF) relationships of rainfall events. In this sense, for locations with rainfall data, it is tried to adjust heavy rainfall equations, also called IDF equations, as in the studies of Fendrich (1998), Costa & Brito (1999), Oliveira et al. (2000), Silva et al. (2002), Back (2002), and Bazzano et al. (2007). IDF relationships are obtained by statistical analyses of long data series from rain gauges. In Brazil, it is relatively easy to obtain 1 day rainfall data. Shorter rainfalls are however rarely recorded, due to the lack of measuring equipment and when available, the measurement series are relatively short and contain numerous data errors. Lima et al. (2005), Oliveira et al. (2008), Cassol et al. (2008), and Back (2009) report the lack of long series of precipitation data, mainly of rainfall.
Bell (1969) established empirical relationships between rainfall events with different durations based on partial data series of rain in the USA, Australia, USSR, Puerto Rico, Alaska, South Africa, and Hawaii. The theoretical basis of this study is the existence of similarity between the mechanisms of these rainfall events. The maximum value of these rainfall events is associated to convective cells that have similar characteristics in many parts of the world and for this reason Bell´s equation is used to estimate maximum precipitation within the specified limits. Bell´s work is a classic of the field, but somewhat restricted, e.g., by the fact that the equations were obtained from rainfall data from several parts of the world, so that their results are based on averages rather than specific values for a location. A second limitation of Bell's method is that the value of maximum rainfall obtained is valid only for durations from 5 to 120 min, and a third is the need to know the maximum 1 h rainfall and 10 year return period, which requires data from a rain-gauge station. Uehara et al. (1980) adjusted the coefficients of Bell's equation for Brazil and for São Paulo State. These equations are limited by the fact that they are based on old data. According to Genovez & Zuffo (2000), the methods based on the relationships between heavy rainfall events of different duration have regional validity, although the average values of these relationships are very similar for several parts of the world. For local estimates, new coefficients should be established, related to the local climate characteristics. Thus, the purpose of this study was to calculate the coefficients of Bell's (1969) equations considering the annual maximum rainfall for rain events of 1 hour and of 1 day for a 10 year return period for the weather stations in the state of Santa Catarina, Brazil.
MATERIAL AND METHODS
Daily graphs of meteorological stations of the Brazilian Institute of Meteorology (Inmet) and of an agricultural institute of research and rural extension of Santa Catarina (Epagri) from locations listed in table 1 were used (Figure 1).
Rainfall charts were digitized and annual maximum rainfall series were determined (durations of 5, 10, 15, 20, 25, 30, 45, 50, 55, 60, 75, 90, 105, 120, 180, 240, 300, 360, 420, 480, 600, 720, 840, 960, 1080, 1200, 1320, and 1440 min) as well as the maximum 1 day rainfall, defined as the total rainfall of one day, accumulated until 9 a.m. For each duration, rainfall values for return periods of 2, 5, 10, 20, 25, 50, and 100 years were also estimated, applying the Gumbel-Chow distribution as:
where: PT - maximum precipitation with return period T (years, mm); - average of annual maximum series (mm); S - standard deviation of annual maximum series (mm); Yn - average of reduced variable; Sn - standard deviation of reduced variable; and Y - reduced variable estimated by:
where: T is the return period (years).
Coefficients of the equation used by Bell (1969) were adjusted, with the height of the estimated rainfall events, in relation to the maximum rainfall for a return period of 10 years and duration of 1 h (60 min) (P1060), as:
where: a, b, c, d, e - empirical coefficients; - estimated precipitation for return period T (years) and duration t (min); - observed precipitation for a return period of 10 years and duration of 1 h (60 min).
Results were compared to the estimates obtained with the coefficients used by Bell (1969) given by:
and the equation proposed by Uehara et al. (1980) for Brazil (based on data of 26 weather stations) and São Paulo with the respective mathematical expressions:
The great majority of rainfall data is measured in pluviometers, which provide only the daily total precipitation. To enable the estimation of design rainfalls with surface runof from maximum rainfall events measured in pluviometers, the coefficients of Bell's equation were also adjusted in relation to the maximum rainfall with a return period of 10 years and duration of 1 day, given by:
where: - estimated precipitation for a return period of 10 years and duration of 1 day.
The coefficients were adjusted according to the methodology described by Back (2009), minimizing the function S, given by:
where: S - sum of squared deviations between precipitation estimated by the equation 3 or 7 and the observed precipitation; t - rainfall duration, min; T - return period, years; G - estimated precipitation using the Gumbel-Chow distribution for duration t and return period T.
The coefficients of equation were adjusted for each one of the weather stations and the coefficients of the general equation were also obtained covering all weather stations of Santa Catarina. Since the results indicated different relationships, the stations of Itajaí, Florianópolis and Laguna were grouped in "Coastal stations", for lying at a distance of less than 10 km from the Atlantic Ocean. The others were all grouped as "Interior stations". Thus, the general equations were adjusted and evaluated for the groups of the stations of the coastal area and the interior of the State.
RESULTS AND DISCUSSION
The relationship between 1 h precipitation and 24 h precipitation (Table 2) varied from 0.22 in Florianópolis and Laguna to 0.45 at the station of Nova Itá, with an average value of 0.35. The Environmental Sanitation Technology Company (CETESB, 1986) indicates an average relationship of 0.42 for several weather stations of Brazil, which is higher than that observed for most stations of Santa Catarina. These differences can be explained by the fact at that some stations extreme values of heavy rainfall with 24 h duration were measured, resulting in a lower relationship between 1 h and 24 h rainfall. The relationship between maximum rainfall of 24 h and rainfall of 1 day varied from 1.30 in Porto União to 1.07 in Chapecó, with a mean of 1.20. Occhipinti & Santos (1966), using a data series from 1928 to 1965, of São Paulo, observed a mean relationship of 1.14. Genovez & Pegoraro (1998), analyzing data of 23 rain gauges of São Paulo State, obtained a minimum value of 1.03 and maximum value of 1.17, with an average value of 1.13. The environmental sanitation technology company - Cetesb (1986) indicated the value of 1.14 based on data of several rainfalls of Brazil.
The coefficients of the equations for each station were adjusted for the estimate of rainfall relationship for return period T and duration t, in terms of the estimated rainfall for 10 year return period and 60 min duration (P1060) (Table 3). The coefficients of the general equation for Santa Catarina were also adjusted considering data from all stations. The relationships between rainfall events with a 10 year return period and duration t min (P10t) obtained in this study for Santa Catarina were compared with the global values obtained by Bell (1969) and the values proposed for Brazil and São Paulo by Uehara et al. (1980) (Figure 2). The standard error values of rainfall estimates were adjusted with each equation and also standard error of estimate in relation to the values estimated by the equations with global validity for Brazil and the general equation for Santa Catarina as well as the equations valid for the coastal area (Itajaí, Florianópolis and Laguna) and for the interior (other stations) of Santa Catarina (Table 4). Using the equation adjusted for Santa Catarina and for 10 min rainfalls, the expression estimates rainfall around 16 % to 19 % lower than Bell's estimates, for return periods from 5 to 100 years, respectively.
When compared to the values adjusted for Brazil, the differences of precipitation are lower than 3 %. Genovez et al. (1994) compared five generalized IDF equations, including the ones proposed by Bell (1969) and by Uehara et al. (1980) for Brazil and São Paulo, with data from nine weather stations of São Paulo State, and verified that the results of generalized IDF equations were similar and in general, the absolute deviations from the estimated values by the equations with the observed values were lower than 15 %, although for some stations this deviation reached 45 %. When adjusting the equations for each weather station, a reduction of the standard error of estimate higher than 50 % was observed at most stations. The data (Table 4) confirm that the establishment of local equations results in more accurate rainfall estimates than the regional or global equations. On the other hand, despite the differences between the values of rainfall relationships, the behavior pattern is similar (Figure 2).
In figure 3 there are the relationships between rainfall of duration t min and 10 year return period with rainfall P1060, where a differentiated behavior between the stations is observed. For durations up to 60 min the obtained coefficients for the stations of Itajaí, Florianópolis and Laguna estimate lower relationships than the ones obtained for the other stations and for durations higher than 60 min. Those stations measured values higher than the other ones (Figure 3a). For 10 year return period and 10 min rainfalls, the relationships vary from 0.31 in Florianópolis to 0.49 in Porto União. For 120 min rainfalls the relationships vary from 1.58 in Florianópolis to 1.20 in Urussanga. For rainfalls longer than 120 min the highest relationships were obtained for the stations of Florianópolis and Itajaí with differences above 30 % of the average value obtained for Santa Catarina (Figure 3b). These variations can be caused by the fact that some stations measure extreme values of heavy rainfall for 24 h duration, resulting in a lower relationship between 1 h rainfall and 24 h rainfall, as observed for the stations of Itajaí and Florianópolis (Table 2).
The coefficients adjusted for rainfall with a 10 year return period and 1 day duration (P101day) (Table 5). As the values of standard error of estimate and coefficients R2 obtained were equivalent to those obtained in Table 3, it was observed that the coefficients of the equation presented in Table 5 can be used in locations where only information of rainfall for 1 day duration is available. These adjusted equations have practical application for estimates of heavy rainfall of short duration in locations where there are only data of daily rainfall observed with rain gauges.
1. The relationships between heavy rainfalls lasting between 10 and 120 min for Santa Catarina show differences below 20 % of the relationships obtained by Bell (1969) and differences below 3 % of the relationships indicated by Uehara et al. (1980) for Brazil.
2. Equations of Itajaí, Florianópolis and Laguna located in the coastal area of Santa Catarina present differentiated relationships from the other stations of the interior of the State.
3. The relationship of maximum rainfall for a 10 year return period and 24 h duration and of the maximum rainfall in 1 day of the different stations varied between 1.07 and 1.30, with an average of 1.20.
4. The relationships based on the maximum rainfall for 1 day duration and 10 year return period present a standard error of estimate equivalent to that obtained based on rainfall for 1 h duration and 10 year return period and represents an alternative to estimate heavy rainfalls at locations where no rainfall data are available.
BACK, A.J. Relações entre precipitações intensas de diferentes durações ocorridas no município de Urussanga, SC. R. Bras. Eng. Agríc. Amb., 13:170-175, 2009. [ Links ]
BACK, A.J. Seleção de distribuição de probabilidade para chuvas diárias extremas do Estado de Santa Catarina. R. Bras. Meteorol., 16:211-222, 2001. [ Links ]
BACK, Á.J. Chuvas intensas e chuva de projeto de drenagem superficial para o Estado de Santa Catarina. Florianópolis, Epagri, 2002. 56p. [ Links ]
BAZZANO, M.G.P.; ELTZ, F.L.F. & CASSOL, E.A. Erosividade, coeficiente de chuvas, padrões e período de retorno das chuvas de Quaraí, RS. R. Bras. Ci. Solo, 31:1205-1217, 2007. [ Links ]
BELL, F.G. Generalized rainfall-duration-frequency relationships. J. Hydr. Division - ASCE, 95:311-327, 1969. [ Links ]
CASSOL, A.E.; ELTZ, F.L.F.; MARTINS, D.; LEMOS, A.M.; LIMA, V.S. & BUENO, A.C. Erosividade, padrões hidrológicos, período de retorno e probabilidade de ocorrências das chuvas em São Borja, RS. R. Bras. Ci. Solo, 32:1239-1251, 2008. [ Links ]
COMPANHIA DE TECNOLOGIA DE SANEAMENTO AMBIENTAL - CETESB. Drenagem urbana - manual de projeto. São Paulo, 1986. 464p. [ Links ]
COSTA, A.R. & BRITO, V.F. Equações de chuva intensa para Goiás e Sul de Tocantins. In: SIMPÓSIO BRASILEIRO DE RECURSOS HÍDRICOS, 13. Belo Horizonte, 1999. Anais... Belo Horizonte: ABRH, 1999. CD ROM. [ Links ]
FENDRICH, R. Chuvas intensas no Estado do Paraná (cidades de Palmas e São Miguel do Iguaçu). Curitiba, PUCPR, 1998. 57p. [ Links ]
GENOVEZ, A.M. & PEGORARO, R.C.F. Análise e avaliação de equações de chuvas intensas generalizadas. In: CONGRESSO LATINOAMERICANO DE HIDRÁULICA, 18., Oaxaca, 1998. Anais... Oaxaca, IAHR, 1998. v.1. p.369-379. [ Links ]
GENOVEZ, A.M.; ZUFFO, A.C. & BORRI GENOVEZ, A.I. Relação entre chuvas intensas de diferente duração e avaliação das equações de chuvas generalizadas. In: CONGRESSO LATINOAMERICANO DE HIDRÁULICA, 16., Santiago, 1994. Anais... Santiago, IAHR, 1994. v.3. p.279-290. [ Links ]
GENOVEZ, A.M. & ZUFFO, A.C. Chuvas intensas no Estado de São Paulo: Estudos existentes e análise comparativa. R. Bras. Rec. Hídr., 5:45-58, 2000. [ Links ]
KITE, G.H. Frequency and risk analyses in hydrology. Fort Collins, Water Resources Publications, 1978. 224p. [ Links ]
LIMA, H.M.F.; MATA, I.P. & LIMA, A.V.F. Aplicação e validação de um simulador estocástico de variáveis climáticas: o caso da precipitação. Ing. Água, 12:27-37, 2005. [ Links ]
OCCHIPINTI, A.G. & SANTOS, P.M. Relações entre as precipitações máximas de um dia e de 24 horas na cidade de São Paulo. São Paulo, IAU/USP, 1966. 10p. [ Links ]
OLIVEIRA, L.F.C.; ANTONINI, J.C.; GRIEBELER, N.; FIOREZE, A.P. & SILVA, M.A.S. Métodos de estimativa de precipitação máxima para o Estado de Goiás. R. Bras. Eng. Agríc. Amb., 12:620-625, 2008. [ Links ]
OLIVEIRA. L.F.; CORTES, F.C.; BARBOSA, F.O.A.; ROMÃO, P.A. & CARVALHO, D.F. Estimativa das equações de chuvas intensas para algumas localidades de Goiás pelo método da desagregação de chuvas. Pesq. Agropec. Trop., 30:23-27, 2000. [ Links ]
SILVA, D.D.; GOMES FILHO, R.R.; PRUSKI, F.F.; PEREIRA, S.B. & NOVAES, L.F. Chuvas intensas no Estado da Bahia. R. Bras. Eng. Agric. Amb., 6:362-367, 2002. [ Links ]
UEHARA, K.; ZAHED FILHO, K.; SILVEIRA, L.N.L. & PALERMO, M.A. Pequenas bacias hidrográficas do Estado de São Paulo: Estudo de vazões médias e máximas. São Paulo, FDTE/Escola Politécnica da USP, 1980. v.2. 780p. [ Links ]
1 Received for publication in October 4, 2011 and approved in February 14, 2012.