Acessibilidade / Reportar erro

Erosividade, coeficiente de chuva, padrões e período de retorno das chuvas de Quaraí, RS

Erosivity, rainfall coefficient and patterns and return period in Quarai, RS, Brazi

Resumos

O conhecimento da potencialidade das chuvas em causar erosão é necessário para planejamento de atividades agrícolas e de engenharia civil. Para a localidade de Quaraí (RS), foram determinados a erosividade da chuva e a relação com a precipitação e o coeficiente de chuva, os padrões hidrológicos e o período de retorno das chuvas. Utilizaram-se dados pluviográficos diários do período 1966-2003. Para cada chuva erosiva, foram separados os segmentos do pluviograma com a mesma intensidade e registrados os dados em planilha. Com o programa Chuveros, foram calculadas as erosividades mensal, anual e média das chuvas pelo índice EI30, no Sistema Internacional de Unidades, e os padrões hidrológicos de chuva, bem como o coeficiente de chuva. Foram realizadas correlações de Pearson e regressões lineares simples entre o índice de erosividade EI30 e os valores médios mensais (p) e anuais (P) de precipitação e do coeficiente de chuva (Rc). Foi calculada a intensidade máxima da chuva pelo método da distribuição extrema tipo 1 para durações de chuva de 1/6, 1/3, 1/2, 1, 2, 4, 8, 12, 24 e 48 h e períodos de retorno da chuva de 2, 5, 10, 20, 50 e 100 anos. Foram ajustadas equações que relacionam a intensidade máxima e a duração da chuva para os períodos de retorno da chuva de 2, 5, 10, 20, 50 e 100 anos, pelo método de regressão linear simples, e construído o gráfico que relaciona essas características da chuva. O valor médio anual de EI30 (fator R da USLE) calculado para Quaraí foi de 9.292 MJ mm ha-1 h-1 ano-1. Obtiveram-se as equações EI30 = -754,37 + 13,50 p (r² = 0,85) e EI30 = -47,35 + 82,72 Rc (R² = 0,84). Em relação ao total das chuvas estudadas, 44 % do número e 90 % do volume foram erosivas. Do número total das chuvas erosivas, 51 % foram do padrão hidrológico avançado, 25 % do intermediário e 24 % do atrasado, ao passo que, do volume total das chuvas erosivas, 57 % foram do padrão avançado, 25 % do intermediário e 18 % do atrasado. Das chuvas erosivas, 57 % da erosividade correspondeu a chuvas do padrão avançado, 25 % a chuvas do padrão intermediário e 18 % a chuvas do padrão atrasado.

energia cinética; intensidade da chuva; duração da chuva; freqüência da chuva; índice EI30


The planning of agricultural and civil engineering activities must be based on knowledge of rainfall erosion potential. For Quarai, RS, Brazil, the rainfall erosivity and its relationship with precipitation and rainfall coefficient, rainfall patterns and rainfall return period were determined. Daily rainfall data from the 1966-2003 period were used. For each erosive rainfall the segments of the rainfall chart with the same intensity were separated and the data recorded in a worksheet. The mean precipitation and rainfall patterns were estimated, as well as the monthly and annual erosivity by the EI30 index (International System of Units), using the software Chuveros. The rainfall coefficient was calculated. Pearson correlations and linear regressions between the erosivity index EI30 and the mean annual precipitation and rainfall coefficient (Rc) were performed. The rainfall maximum intensity was calculated by the method of extreme distribution type I for different rainfall durations (1/6, 1/3, 1/2, 1, 2, 4, 8, 12, 24 and 48 h) and rainfall return periods (2, 5, 10, 20, 50 and 100 years). Equations were adjusted that relate the maximum intensity and rainfall duration to all return periods, by the method of simple linear regression, and the rainfall characteristics related in a graph. The mean annual values of EI30 (R factor of USLE) determined for Quarai, RS, Brazil, was 9,292 MJ mm ha-1 h-1 year-1. The following equations were obtained: EI30 = -754.37 + 13.50 p (r² = 0.85) and EI30 = -47.35 + 82.72 Rc (r² = 0.84). In relation to the total precipitation studied, 44 % of the number of rainfalls and 90 % of the volume were erosive. Of the total rainfalls evaluated, 44 % of the number and 90 % of the volume were erosive. Of the total erosive rainfall events, 51 % had an advanced, 25 % had an intermediate, and 24 % had a delayed hydrological pattern. Of the total volume of erosive rainfalls, 57 % showed and advanced pattern, 25 % intermediate, and 18 % a delayed pattern.

kinetic energy; rainfall; rainfall intensity; rainfall duration; rainfall frequency; EI30 index; rainfall coefficient


NOTA

Erosividade, coeficiente de chuva, padrões e período de retorno das chuvas de Quaraí, RS1 1 Parte da Tese de Mestrado em Ciência do Solo apresentada pelo primeiro autor à Universidade Federal de Santa Maria – UFSM.

Erosivity, rainfall coefficient and patterns and return period in Quarai, RS, Brazil

Marcos Gabriel Peñalva BazzanoI; Flávio Luiz Foletto EltzII; Elemar Antonino CassolIII

IEngenheiro-Agrônomo, MSc. Asociación para la Agricultura Sostenible en base a Siembra Directa. La Paz 302, Col. Alameda, C.P. 38050, Celaya, Guanajuato, México. E-mail: marcos.penalva@gmail.com

IIProfessor Titular do Departamento de Solos, Universidade Federal de Santa Maria – UFSM. CEP 97105-900 Santa Maria – RS. Bolsista do CNPq. feltz@ccr.ufsm.br

IIIProfessor Adjunto do Departamento de Solos, Faculdade de Agronomia, Universidade Federal de Rio Grande do Sul – UFRGS. CEP 90001-970 Porto Alegre (RS). Bolsista do CNPq. E-mail: cassolea@orion.ufrgs.br

RESUMO

O conhecimento da potencialidade das chuvas em causar erosão é necessário para planejamento de atividades agrícolas e de engenharia civil. Para a localidade de Quaraí (RS), foram determinados a erosividade da chuva e a relação com a precipitação e o coeficiente de chuva, os padrões hidrológicos e o período de retorno das chuvas. Utilizaram-se dados pluviográficos diários do período 1966-2003. Para cada chuva erosiva, foram separados os segmentos do pluviograma com a mesma intensidade e registrados os dados em planilha. Com o programa Chuveros, foram calculadas as erosividades mensal, anual e média das chuvas pelo índice EI30, no Sistema Internacional de Unidades, e os padrões hidrológicos de chuva, bem como o coeficiente de chuva. Foram realizadas correlações de Pearson e regressões lineares simples entre o índice de erosividade EI30 e os valores médios mensais (p) e anuais (P) de precipitação e do coeficiente de chuva (Rc). Foi calculada a intensidade máxima da chuva pelo método da distribuição extrema tipo 1 para durações de chuva de 1/6, 1/3, 1/2, 1, 2, 4, 8, 12, 24 e 48 h e períodos de retorno da chuva de 2, 5, 10, 20, 50 e 100 anos. Foram ajustadas equações que relacionam a intensidade máxima e a duração da chuva para os períodos de retorno da chuva de 2, 5, 10, 20, 50 e 100 anos, pelo método de regressão linear simples, e construído o gráfico que relaciona essas características da chuva. O valor médio anual de EI30 (fator R da USLE) calculado para Quaraí foi de 9.292 MJ mm ha-1 h -1 ano-1. Obtiveram-se as equações EI30 = -754,37 + 13,50 p (r2 = 0,85) e EI30 = -47,35 + 82,72 Rc (R2 = 0,84). Em relação ao total das chuvas estudadas, 44 % do número e 90 % do volume foram erosivas. Do número total das chuvas erosivas, 51 % foram do padrão hidrológico avançado, 25 % do intermediário e 24 % do atrasado, ao passo que, do volume total das chuvas erosivas, 57 % foram do padrão avançado, 25 % do intermediário e 18 % do atrasado. Das chuvas erosivas, 57 % da erosividade correspondeu a chuvas do padrão avançado, 25 % a chuvas do padrão intermediário e 18 % a chuvas do padrão atrasado.

Termos de indexação: energia cinética, intensidade da chuva, duração da chuva, freqüência da chuva, índice EI30.

SUMMARY

The planning of agricultural and civil engineering activities must be based on knowledge of rainfall erosion potential. For Quarai, RS, Brazil, the rainfall erosivity and its relationship with precipitation and rainfall coefficient, rainfall patterns and rainfall return period were determined. Daily rainfall data from the 1966-2003 period were used. For each erosive rainfall the segments of the rainfall chart with the same intensity were separated and the data recorded in a worksheet. The mean precipitation and rainfall patterns were estimated, as well as the monthly and annual erosivity by the EI30 index (International System of Units), using the software Chuveros. The rainfall coefficient was calculated. Pearson correlations and linear regressions between the erosivity index EI30 and the mean annual precipitation and rainfall coefficient (Rc) were performed. The rainfall maximum intensity was calculated by the method of extreme distribution type I for different rainfall durations (1/6, 1/3, 1/2, 1, 2, 4, 8, 12, 24 and 48 h) and rainfall return periods (2, 5, 10, 20, 50 and 100 years). Equations were adjusted that relate the maximum intensity and rainfall duration to all return periods, by the method of simple linear regression, and the rainfall characteristics related in a graph. The mean annual values of EI30 (R factor of USLE) determined for Quarai, RS, Brazil, was 9,292 MJ mm ha-1 h -1 year-1. The following equations were obtained: EI30 = -754.37 + 13.50 p (r2 = 0.85) and EI30 = -47.35 + 82.72 Rc (r2 = 0.84). In relation to the total precipitation studied, 44 % of the number of rainfalls and 90 % of the volume were erosive. Of the total rainfalls evaluated, 44 % of the number and 90 % of the volume were erosive. Of the total erosive rainfall events, 51 % had an advanced, 25 % had an intermediate, and 24 % had a delayed hydrological pattern. Of the total volume of erosive rainfalls, 57 % showed and advanced pattern, 25 % intermediate, and 18 % a delayed pattern.

Index terms: kinetic energy, rainfall, rainfall intensity, rainfall duration, rainfall frequency, EI30 index, rainfall coefficient.

INTRODUÇÃO

A erosão hídrica é a conseqüência da interação da chuva com o solo, ou seja, é a resultante do poder da chuva em causar erosão e da capacidade do solo em resistir a esta. Ao potencial da chuva em causar erosão chama-se erosividade.

Na primeira metade do século passado foram iniciados esforços para estimar as perdas de solo sob diferentes condições de clima, solo e manejos, visando a um planejamento conservacionista. Em 1954, surgiu nos Estados Unidos, a partir do aperfeiçoamento de diferentes modelos, a Equação Universal de Perdas de Solo (USLE), também conhecida como equação de Wischmeier & Smith (1978), na qual o fator R, índice de erosividade, expressa o potencial da chuva de provocar a erosão do solo, sendo variável de região para região.

Segundo Wischmeier (1959), o índice de erosividade EI30 expressa a relação entre a intensidade e a energia cinética da chuva e se correlaciona bem com as perdas de solo por erosão. Morais et al. (1988a), estudando três localidades do Rio Grande do Sul, concluíram que o índice de erosividade EI30 (Wischmeier & Smith, 1978) foi adequado para estimar o potencial erosivo da chuva.

Morais et al. (1988a), Lago (2000), Cassol et al. (2004a,b) e Roncato et al. (2004) determinaram o índice EI30 de várias localidades do Rio Grande do Sul a partir de dados pluviográficos. A determinação desse índice é um processo trabalhoso, que implica ter pluviogramas de vários anos para determinar a energia cinética e intensidade máxima das chuvas. Em muitas ocasiões não se contam com registros de pluviógrafo e sim com dados de pluviômetro. A fim de facilitar a determinação do índice de erosividade, Oliveira Júnior & Medina (1990), Morais et al. (1988b, 1991), Silva & Dias (2003), entre outros, relacionaram o índice EI30 com a precipitação média anual (P) e o coeficiente da chuva expresso como Rc (igual ao quadrado da precipitação média mensal, p2, dividido pela precipitação média anual, P), proposto por Fournier (1956) e modificado por Lombardi Neto (1977). Em várias localidades, as correlações entre esses parâmetros não foram significativas (Albuquerque, 1992; Bertol, 1993, 1994).

Perdas de solo e danos em infra-estrutura e obras de engenharia hidráulica ocorrem freqüentemente por causa da chuva. Eficientes dimensionamentos de obras hidráulicas e conservacionistas podem ser realizados quando são considerados aspectos de risco e freqüência da chuva. A análise de distribuição de freqüência permite estimar o período de retorno de um evento máximo anual de chuva, que pode ser definido como o período, em anos, em que o referido evento pode ser igualado ou excedido, pelo menos uma vez, em média (Eltz et al., 1992). Além disso, os padrões hidrológicos de chuva, classificados em padrão avançado, intermediário e atrasado, quando a localização dos picos de maior intensidade está no início, meio ou fim do período de duração das chuvas, respectivamente (Horner & Jens, 1941), podem mudar de uma região para outra, causando diferentes perdas de solo (Mehl et al., 2001), e ter diferentes períodos de retorno (Mehl et al., 2001; Moreti et al., 2003).

O conhecimento das características da chuva permite a planificação mais segura de estruturas de conservação de solo (terraços, curvas de nível) e de práticas agrícolas que visem à conservação do solo por meio de manutenção de sua cobertura, assim como outras obras (barragens, canais escoadouros), e de estruturas hidráulicas de fluxo para águas pluviais, o que justifica sua determinação.

Assim, os objetivos deste trabalho foram determinar: erosividade da chuva e relação com a precipitação e o coeficiente de chuva; padrões hidrológicos da chuva; e período de retorno da chuva para a localidade de Quaraí, no Estado do Rio Grande do Sul.

MATERIAL E MÉTODOS

Foram calculados a erosividade da chuva, o coeficiente de chuva e sua relação com a erosividade, a relação da erosividade com a precipitação, determinados os padrões e os períodos de retorno das chuvas da localidade de Quaraí, no Estado de Rio Grande do Sul, a partir de dados pluviográficos da estação meteorológica de Quaraí, os quais foram cedidos pela Fepagro (Fundação Estadual de Pesquisa Agropecuária), de uma série de 38 anos, correspondente ao período 1966-2003.

Quaraí localiza-se na região fisiográfica denominada Campanha. A estação meteorológica do local está situada a uma altitude de 100 m, latitude sul de 30 º 23 ' e longitude oeste de 56 º 26 '. Segundo a classificação de Köppen, o clima é Cfa, ou seja, temperado úmido sem estiagem, com temperaturas média, máxima e mínima anual de 19,6, 24,5 e 13,7 ºC, respectivamente, e precipitação média anual de 1.513 mm, no período 1971-2000 (Ipagro, 1989).

Conforme sugerido por Cabeda (1976), foram selecionadas as chuvas individuais erosivas, considerando-se chuvas individuais aquelas separadas de outras por intervalo de no mínimo seis horas com precipitação inferior a 1 mm e erosivas quando a precipitação total for igual ou superior a 10 mm ou quando a precipitação for igual ou superior a 6 mm em 15 min ou menos de chuva.

A partir dos pluviogramas diários, cada chuva erosiva individual foi separada em segmentos uniformes, ou seja, segmentos de mesma inclinação (intensidade constante). Os dados das chuvas foram registrados manualmente em planilha, segundo Freitas & Madeira Neto (1980), onde foram anotadas as horas em que houve variação de intensidade e a chuva acumulada até aquele momento. A amplitude de registro dos pluviogramas foi de 10 mm de precipitação, com unidade de 0,1 mm. O tempo de registro foi de 24 h, com unidades de 10 min, realizando-se a leitura da precipitação com unidades de 5 min.

Por meio do programa Chuveros, elaborado por Elemar Antonino Cassol da UFRGS (Universidade Federal do Rio Grande do Sul), e planilhas de cálculo, foram estimadas as erosividades mensal, anual e média das chuvas pelo índice EI30 (Wischmeier & Smith, 1958), no Sistema Internacional de Unidades (Foster et al., 1981). Os cálculos do programa foram realizados da forma explicitada a seguir.

A energia cinética de cada segmento uniforme de uma chuva erosiva individual é determinada por:

em que ECs é a energia cinética no segmento de chuva, em MJ ha-1, EC é a unidade de energia cinética e h é a quantidade de chuva no segmento uniforme, em mm.

A equação para cálculo da energia cinética é a seguinte:

em que EC é a energia cinética por mm de chuva, em MJ ha-1 mm-1, e I é a intensidade da chuva no segmento uniforme (intensidade constante), em mm h-1.

A energia cinética total da chuva erosiva individual foi calculada pela expressão:

em que ECt é a energia cinética total da chuva erosiva individual, em MJ ha-1, e S ECs é o somatório da energia cinética de cada segmento da chuva, em MJ ha-1.

O índice de erosividade de cada chuva foi calculado por:

em que EI30 é o índice de erosividade da chuva erosiva individual, em MJ mm ha-1 h-1 ; ECt é a energia cinética total da chuva, em MJ ha-1; e I30 é a intensidade máxima da chuva em um período de 30 min, em mm h-1, ou seja, a quantidade máxima de chuva em um período de 30 min multiplicada por 2, determinada a partir dos dados obtidos do pluviograma.

Os valores mensais dos índices de erosividade EI30 são a média dos meses para a localidade. O valor médio anual do índice de erosividade EI30 é a média dos anos do período estudado.

Os valores médios mensais da precipitação e do índice de erosividade foram expressos como percentagens do valor médio anual da precipitação e do índice de erosividade, respectivamente, para obter a curva de distribuição acumulada mensal da precipitação e do índice de erosividade, conforme Wischmeier & Smith (1958).

O coeficiente de chuva (Rc) foi calculado por meio da equação proposta por Lombardi Neto (1977):

em que p representa a precipitação média mensal (mm) e P a precipitação média anual (mm).

Foram realizadas análises de correlação de Pearson e de regressão linear simples entre o índice de erosividade EI30 e cada um dos valores médios anuais de precipitação e coeficiente de chuva, utilizando o programa estatístico SAS (Statistical Analysis System), versão 8.02 para Windows, com o apoio do manual para análises estatísticas de dados experimentais (Rebolledo, 2002).

O padrão da chuva foi determinado pelo programa Chuveros, separando-se as chuvas de acordo com os padrões de chuva avançado, intermediário e atrasado, da classificação proposta por Horner & Jens (1941). Padrão avançado, intermediário e atrasado caracteriza-se quando a localização do pico de maior intensidade está no início (primeiro terço), meio (segundo terço) ou fim (terceiro terço) do período de duração da chuva, respectivamente. Também foi calculada a percentagem do total das chuvas que representa cada padrão.

O período de retorno (T), definido como o período em anos em que um evento máximo anual pode ser igualado ou excedido pelo menos uma vez, em média, a cada T anos, foi determinado segundo o método descrito em Eltz et al. (1992).

De acordo com Bedient & Huber (1988), a função de densidade probabilística pode ser vista como:

em que F é uma função de densidade probabilística; e, a base neperiana (2,712817); a, um parâmetro de escala; e u (modal), um parâmetro de localização. Esses parâmetros se relacionam com a média (µ) e o desvio-padrão (s) da intensidade máxima da chuva por:

e

O valor da intensidade máxima esperada para determinado período de retorno e duração da chuva é dado por:

sendo K o fator de freqüência, que se calcula por:

em que ln é o logaritmo neperiano e T o período de retorno.

O período de retorno, T, desejado é calculado pela equação:

Assim, é possível calcular o período de retorno para uma determinada intensidade e duração da chuva.

Por meio do programa Chuveros, geraram-se os valores das quantidades máximas de chuva ocorridas durante 1/6, 1/3, 1/2, 1, 2 e 4 h, e a partir dos pluviogramas das chuvas de cada ano foram determinadas as quantidades máximas de chuva ocorridas durante 8, 12, 24 e 48 h em qualquer momento do período de tempo analisado. Nos intervalos menores (até 4 h), as durações de chuva consideradas podem representar um período mais curto que a duração total de uma chuva individual, enquanto as durações maiores (oito ou mais horas) podem abranger mais de uma chuva individual. Os dados das intensidades máximas em cada duração de chuva foram testados pela distribuição extrema tipo I (Gumbel, 1958), com o programa Freqrisk, descrito em Kite (1988), que calculou a média (µ), o desvio-padrão (s), o parâmetro a e o período de retorno para 2, 5, 10, 20, 50 e 100 anos. Para cada período de retorno, foram ajustadas equações polinomiais de segunda ordem com a forma:

em que log é o logaritmo base 10; I, a intensidade máxima da chuva, em mm h-1; D, a duração da chuva, em horas; e a, b e c, parâmetros de ajuste. Também foi ajustada a equação utilizada por Denardin et al. (1980), que relaciona as características intensidade da chuva, duração da chuva e período de retorno da chuva:

em que I é a intensidade máxima da chuva, em mm h-1; T, o período de retorno, em anos; t, a duração da chuva, em horas; e a, b, c e d, parâmetros de ajuste. O ajuste da equação (13) foi realizado por regressão múltipla linear após linearização com a forma:

em que log é o logaritmo-base 10. Em virtude da associação com a variável independente t na equação (14), o parâmetro c não pôde ser estimado. Assim, atribuíram-se valores arbitrários ao coeficiente c, estimando-se os coeficientes a, b e d para cada caso, até se encontrar a equação com maior coeficiente de determinação múltipla (R2), considerada a de melhor ajuste.

RESULTADOS E DISCUSSÃO

Erosividade da chuva (índice EI30) e coeficiente de chuva (Rc)

No quadro 1 são apresentados os valores mensais, anuais e médios do índice de erosividade EI30 da chuva, com os respectivos desvios-padrão e coeficientes de variação, bem como a amplitude de variação. O valor do fator R da USLE para Quaraí – valor médio anual da erosividade da chuva – é de 9.292 MJ mm ha-1 h -1 ano-1. No período dezembro-março observou-se a maior erosividade média mensal, com valores do índice EI30 que variaram entre 1.072 e 1.304 MJ mm ha-1 h -1 mês-1. Nos meses de julho e agosto ocorreu, em média, a menor erosividade (286 e 211 MJ mm ha-1 h -1 mês-1, respectivamente). O coeficiente de variação dos valores mensais foi muito alto, variando entre 91 e 169 %, enquanto o coeficiente de variação anual, embora menor, foi também muito alto, alcançando 48 %, indicando que as características da chuva que determinam seu potencial erosivo, principalmente intensidade e quantidade, são muito variáveis. No período em estudo, o máximo valor de erosividade anual, 22.828 MJ mm ha-1 h -1 ano-1, foi observado em 2002, e a menor erosividade, 2.239 MJ mm ha-1 h -1 ano-1, em 1995.


O valor da erosividade média anual das chuvas (fator "R" da USLE) obtido neste trabalho para Quaraí-RS (9.292 MJ mm ha-1 h -1 ano-1) está dentro da isoerodente estimada por Lago (2000) para a metade sul do Estado do Rio Grande do Sul. Roncato et al. (2004) determinaram para Santa Maria-RS um índice EI30 de 9.026 MJ mm ha-1 h -1 ano-1, com base no período de 1996 a 2004, similar ao do presente estudo, embora a precipitação média anual dessa localidade seja maior (ao redor de 1.750 mm). Por sua vez, Bertol (1993) determinou para Lages-SC um índice EI30 de 5.694 MJ mm ha-1 h-1 ano-1, embora a precipitação média anual desta localidade (1.674 mm) fosse maior que a determinada para Quaraí (1.359 mm). Já Morais et al. (1991), estudando as localidades do sudoeste do Mato Grosso, Barranquinho e Flechas, com precipitações de 1.369 e 1.316 mm respectivamente, similares às do presente trabalho, encontraram valores médios anuais do índice EI30 de 8.493 e 7.830 MJ mm ha-1 h -1 ano-1, menores que os obtidos para Quaraí.

O quadro 2 apresenta a precipitação, o índice de erosividade e o coeficiente de chuva Rc médios mensais e anuais. A tendência de distribuição dos valores percentuais mensais da precipitação e do coeficiente de chuva (Figura 1) é similar pelo fato de o coeficiente de chuva considerar a precipitação mensal e a anual. Dos três parâmetros apresentados na figura 1, há relativa melhor distribuição temporal das precipitações, comparada às distribuições da erosividade e do coeficiente de chuva. Assim, observam-se valores extremos da precipitação de 4,7 a 11,2 %, que são maiores no período janeiro-abril (9,7–11,2 %) e menores em junho-setembro (4,7–6,6 %), enquanto as distribuições da erosividade e do coeficiente de chuva apresentam maior amplitude entre os valores extremos (2,3–14,0 % e 2,6–14,4 %, respectivamente). Em janeiro-abril observam-se os valores maiores da erosividade (9,0–14,0 %) e do coeficiente de chuva (10,8–14,4 %), e os menores, em junho-setembro, com 2,3–5,2 % e 2,6–4,9 % para erosividade e coeficiente de chuva, respectivamente. Os valores começam altos no mês de janeiro (10,7, 12,3 e 13,0 %, para precipitação, erosividade e coeficiente de chuva, respectivamente); logo, atingem um pico no mês de fevereiro, o qual é o maior do ano para precipitação (11,2 %) e coeficiente de chuva (14,4 %). Continuando no decorrer do ano, os valores vão diminuindo até chegar ao mínimo no mês de agosto (4,7 % para a precipitação, 2,3 % para o índice EI30 e 2,6 % para o coeficiente de chuva). A partir desse mês, os valores começam a subir e o índice EI30 alcança seu valor máximo (14,0 %) no mês de dezembro. No período novembro-março a erosividade ultrapassou a precipitação, sugerindo maior potencial erosivo das chuvas nesses meses, o que contrasta com o período abril-outubro, em que a erosividade ficou abaixo da precipitação, com exceção do mês de maio, que teve valores similares de erosividade e precipitação (9,1 e 8,7 %, respectivamente). Isso indica o efeito preponderante da intensidade da chuva em caracterizar o seu potencial erosivo (Wischmeier & Smith, 1958) e quando se deve ter os maiores cuidados na proteção ao solo.



O período de outubro a maio concentrou 85,4 % da erosividade anual (Figura 2), sendo por isso considerado crítico em relação à erosão. Nesse período, práticas de preparo de solo para semeaduras de culturas de verão e de inverno aumentam os riscos de perdas de solo. Assim, o conhecimento da distribuição da erosividade durante o ano contribui para o planejamento adequado das práticas conservacionistas do solo que permitam reduzir os riscos de erosão.


As análises de correlação de Pearson mostraram alta correlação do índice de erosividade EI30 com a precipitação média mensal (r = 0,92; significativo a 0,01 %) e com o coeficiente de chuva (r = 0,92; significativo a 0,01 %) (Quadro 3).


Correlações significativas entre o índice EI30 e o coeficiente de chuva Rc foram observadas também em outras localidades, como os casos de oito regiões do Estado do Paraná (Rufino et al., 1993), Goiânia-GO (Silva et al., 1997), Fortaleza-CE (Silva & Días, 2003), Campinas-SP (Lombardi Neto, 1977), sudoeste de Mato Grosso do Sul (Morais et al., 1991), Selvíria-MS (Carvalho et al., 1997) e Piraju-SP (Roque et al., 2001). Correlações não-significativas foram encontradas para a localidade de Campos Novos-SC (Bertol, 1994).

As regressões lineares obtidas foram altamente significativas, sendo as equações EI30 = -754,37 + 13,50 p; R2 = 0,85 (significativo a 0,01 %), em que p = precipitação média mensal (mm), e EI30 = -47,35 + 82,72 Rc; R2 = 0,84 (significativo a 0,01 %), em que Rc = coeficiente de chuva = p2 P-1, sendo p = precipitação média mensal (mm) e P = precipitação média anual (mm) (Quadro 3). Assim, é possível utilizar as equações para determinar o índice EI30 por meio do uso de dados pluviométricos mensais e anuais. Essa é uma possibilidade de se utilizarem os dados pluviométricos e as expressões ajustadas no quadro 3 para estimar a erosividade das chuvas em locais da mesma região agroclimática, no entorno de Quaraí.

Padrões de chuva

Aproximadamente metade das chuvas erosivas (51 %) foi classificada dentro do padrão avançado (Quadro 4). O restante das chuvas esteve, quase por partes iguais, classificado no padrão intermediário (25 %) e no padrão atrasado (24 %). Em relação ao número total de chuvas anuais, as chuvas erosivas representam 44 %, enquanto o volume da precipitação erosiva representa 90 % do volume da precipitação total. Isso quer dizer que 66 % do número de chuvas contribui com apenas 10 % do volume de precipitação. Considerando o volume das chuvas erosivas anuais, a maior proporção correspondeu a chuvas do padrão avançado (57 %), seguido do padrão intermediário (25 %) e do padrão atrasado (18 %). Em relação ao padrão hidrológico das chuvas erosivas, 57 % da erosividade correspondeu a chuvas de padrão avançado, 25 % a chuvas de padrão intermediário e 18 % a chuvas de padrão atrasado.


A distribuição das chuvas concentrada no padrão avançado permite esperar menores perdas de solo devido ao fato de que no momento do pico da chuva este estaria menos úmido que no caso dos outros padrões; dessa forma, a desagregação, o selamento e o transporte de solo seriam menores (Mehl et al., 2001). Eltz et al. (2001), trabalhando com simulador de chuva em um Argissolo Vermelho-Amarelo, observaram maior perda de solo em chuvas de padrão atrasado em relação às perdas de solo das chuvas de padrão intermediário e avançado, entre as quais não houve diferenças significativas, o que foi explicado pelas condições de alteração da superfície e pela umidade do solo no decorrer da chuva. No padrão atrasado, quando ocorre o pico de máxima intensidade, o solo está com maior umidade (em relação aos outros padrões), favorecendo a desagregação, o selamento superficial e o transporte de suas partículas, ocorrendo assim maior perda de solo.

Período de retorno

Os valores de precipitação máxima observados para durações da chuva de 1/6 e 1/3 de hora foram de 40,0 e 50,0 mm, respectivamente, no ano de 1973, ao passo que no ano de 2002 ocorreram as máximas precipitações para durações da chuva de 1/2, 1, 2, 4, 8, 12, 24 e 48 h, alcançando, respectivamente, 63,5; 99,5; 151,4; 216,8; 262,8; 265,5; 266,4; e 287,1 mm (Quadro 5).


Em razão do método empregado para obter os dados de precipitação máxima para as diferentes durações da chuva, os valores aqui apresentados poderiam diferir dos obtidos por Ipagro (1989), a qual considera a quantidade de chuva entre um período fixo de tempo, ou seja, 24 h, para valores tanto de pluviômetro quanto de pluviógrafo. No presente trabalho, é considerada a precipitação máxima para a duração estabelecida em qualquer período de tempo, ou seja, é considerada a data do início de uma chuva erosiva individual, independentemente de quantos dias possa ter durado.

No quadro 6 são apresentados os valores de intensidade média, desvio-padrão e parâmetro a para diferentes durações da chuva da localidade de Quaraí, calculados mediante o programa Freqrisk (Kite, 1988).


A intensidade da chuva para durações e períodos de retorno utilizados neste trabalho, analisados pela distribuição extrema tipo I (Gumbel, 1958), é apresentada no quadro 7. Um teste de Kolmogorov-Smirnov realizado mostrou que essa distribuição se ajusta aos dados observados. Os parâmetros das regressões (a, b e c) entre o log da intensidade e o log da duração da chuva podem ser observados no quadro 8. Obteve-se alta significância estatística (probabilidade menor de 0,01 %) e alto coeficiente de determinação (R2 = 0,9998) nas seis regressões dos períodos de retorno considerados. Isso indica que a variação da intensidade como variável da duração é explicada em mais de 99 % pelo modelo quadrático considerado (log I = a + b (log D) + c (log D)2), significando alta confiabilidade nos resultados quando são utilizadas essas equações. Eltz et al. (1992) obtiveram resultados também altamente significativos para a localidade de Santa Maria-RS.



As curvas que geraram os coeficientes das equações anteriores estão representadas na figura 3. Com o uso dessa figura, é possível estimar diretamente a intensidade, duração e período de retorno da chuva. Considerando duração e período de retorno específico, pode-se determinar a intensidade da chuva para essas condições. De forma similar, partindo de dados de duração e intensidade, chega-se ao período de retorno, e, com os valores de intensidade e período de retorno pode-se encontrar a duração da chuva. Dentro das amplitudes de duração, intensidade e períodos de retorno apresentados no presente estudo, podem ser feitas interpolações. Para extrapolações, a confiabilidade é menor. Em razão de a origem de dados corresponder a uma quantidade de anos limitada (38 anos), os eventos pluviais máximos reais podem ser maiores do que os observados.


Foi ajustada a seguinte equação, que relaciona as características intensidade, duração e período de retorno da chuva:

em que I é a intensidade máxima da chuva, em mm h-1; T, o período de retorno, em anos; e t, o período de duração da chuva, em horas.

Essa equação teve alta significância estatística (probabilidade menor que 0,01 %), alto coeficiente de determinação (R2 = 0,9958) e baixo coeficiente de variação (CV = 2,30 %). Os valores da intensidade da chuva para durações e períodos de retorno selecionados calculados pela equação 15 são apresentados no quadro 9. Os valores são similares aos obtidos no quadro 7, o qual indicaria que se pode utilizar qualquer uma das equações apresentadas.


As informações que proporcionam a análise de freqüência apresentada no gráfico de intensidade-duração-freqüência (Figura 3) e a equação 15 poderão ser levadas em consideração para dimensionamento de obras de conservação de solo, como canais de escoamento e espaçamento entre terraços, assim como diferentes obras civis. O dimensionamento mais ajustado ou superestimado de uma determinada obra implica, respectivamente, menor ou maior segurança da obra para resistir aos eventos pluviais, sendo as últimas mais onerosas. Assim, um planejamento eficaz deverá levar em consideração esses aspectos.

CONCLUSÕES

1. O valor do fator R da equação universal de perdas de solo para a localidade de Quaraí é de 9.292 MJ mm ha-1 h-1 ano-1 .

2. Obteve-se alta correlação do índice de erosividade EI30 com a precipitação e com o coeficiente de chuva, podendo-se utilizar as equações obtidas para estimar o fator R da USLE.

3. Em Quaraí, RS, o número, o volume das chuvas e a erosividade (EI30) corresponderam, em sua maioria, a chuvas classificadas no padrão avançado (51, 57 e 57 %), seguido do intermediário (25, 25 e 25 %) e atrasado (24, 18 e 18 %), respectivamente.

4. Em relação ao total das chuvas estudadas, 44% do número e 90 % do volume foram erosivas.

5. Os períodos de retorno da chuva em Quaraí podem ser calculados por meio das equações afins, utilizando os valores dos parâmetros encontrados, ou pelo gráfico das curvas de intensidade-duração-freqüência.

AGRADECIMENTOS

À FEPAGRO, pela cessão dos dados, e aos observadores meteorológicos da FEPAGRO, que, no trabalho anônimo, cotidiano, fundamental e imprescindível, realizaram a leitura e anotação das precipitações, possibilitando assim o uso dos dados e a realização deste trabalho.

LITERATURA CITADA

Recebido para publicação em janeiro de 2006 e aprovado em maio de 2007.

  • ALBUQUERQUE, A.W. Determinação da erosividade de Caruaru-PE: 1. Correlação com perdas de solo, 2. Distribuição de probabilidade de ocorrência. Piracicaba, Escola Superior de Agricultura Luiz de Queiroz, 1992. 100p. (Tese de Mestrado)
  • BEDIENT, P.B. & HUBER, W.C. Hydrology and food plain analysis. New York, Addison Wesley, 1988. 650p.
  • BERTOL, I. Avaliação da erosividade da chuva na localidade de Campos Novos (SC) no período de 1981-1990. Pesq. Agropec. Bras., 29:1453-1458, 1994.
  • BERTOL, I. Índice de erosividade (EI30) para Lages (SC) 1Ş aproximação. Pesq. Agropec. Bras., 28:515-521, 1993.
  • CABEDA, M.S.V. Computation of storm EI values. West Lafayette, Purdue University, 1976. 6p. Não Publicado.
  • CARVALHO, M.P.; CATANEO, A. & LOMBARDI NETO, F. Parâmetros de erosividade da chuva e da enxurrada correlacionados com as perdas de solo e determinação da erodibilidade de um Podzólico Vermelho-Amarelo de Pindorama. R. Bras. Ci. Solo, 21:279-286, 1997.
  • CASSOL, E.A.; MARTINS, D.; ELTZ, F.L.F. & FALLEIRO, R.M. Erosividade das chuvas em Taquari, RS, determinada pelo índice EI30, no período de 1963 a 1999. In: REUNIÃO BRASILEIRA DE MANEJO E CONSERVAÇÃO DO SOLO E DA ÁGUA; MANEJO: INTEGRANDO A CIÊNCIA DO SOLO NA PRODUÇÃO DE ALIMENTOS, 15., Santa Maria, 2004. Resumos expandidos. Santa Maria, Universidade Federal de Santa Maria, 2004b. CD-ROM
  • CASSOL, E.A.; MAZURANA, J.; ELTZ, F.L.F. & DULLIUS, M.V. Erosividade das chuvas em Porto Alegre, RS, determinada pelo índice EI30, com base no período de 1974 a 2003. In: REUNIÃO SUL BRASILEIRA DE CIÊNCIA DO SOLO, 5., Florianópolis, 2004. Anais. Florianópolis, 2004a. v.1. p.1-4.
  • DENARDIN, J.E.; FREITAS, P.L.; WÜNSCHE, W.A. & WENDT, W. Características fundamentais da chuva no Brasil. I. Rio Grande do Sul. Pesq. Agropec. Bras., 15:419-421, 1980.
  • ELTZ, F.L.F.; MEHL, H.U. & REICHERT, J.M. Perdas de solo e água em entressulcos em um Argissolo Vermelho-Amarelo submetido a quatro padrões de chuva. R. Bras. Ci. Solo, 25:485-493, 2001.
  • ELTZ, F.L.F.; REICHERT, J.M. & CASSOL, E.A. Período de retorno de chuvas em Santa Maria, RS. R. Bras. Ci. Solo, 16:265-269, 1992.
  • FOSTER, G.R.; McCOOL, D.K.; RENARD, K.G. & MOLDENHAUER, W.C. Conversion of the universal soil loss equation to SI units. J. Soil Water Conserv., 36:355-359, 1981.
  • FOURNIER, F. The effect of climatic factors on soil erosion: Estimates of solids transported in suspension in runoff. [S.L.] Association Hydrologic. Int. Public., 1956. v.38. 6p.
  • FREITAS, P.L. & MADEIRA NETO, R.S. Um programa computacional para cálculos dos índices de erosividade das chuvas. In: ENCONTRO NACIONAL DE PESQUISA SOBRE CONSERVAÇÃO DE SOLO, 3., Recife, 1980. Resumos. Recife, SUDENE, 1980.
  • GUMBEL, E.J. Statistics of extremes. New York, Columbia University Press, 1958. 375p.
  • HORNER, W.W. & JENS, S.W. Surface runoff determination from rainfall without using coefficients. Trans.ASAE, 107:1039-1117, 1941.
  • INSTITUTO DE PESQUISAS AGRONÔMICAS - IPAGRO. Atlas agroclimático do Estado do Rio Grande do Sul. Porto Alegre, Secretaria da Agricultura e Abastecimento, 1989. v.1. 102p.
  • KITE, G.W. Frequency and risk analysis in hydrology. Littleton, Water Resources Publications, 1988. 257p.
  • LAGO, J.C. Erosividade das chuvas na metade sul do Rio Grande do Sul. Pelotas, Universidade Federal de Pelotas, 2000. 123p. (Tese de Doutorado)
  • LOMBARDI NETO, F. Rainfall erosivity Its distribution and relationship with soil loss at Campinas, Brasil. West Lafayette, Purdue University, 1977. 53p. (Tese de Mestrado)
  • MEHL, H.U.; ELTZ, F.L.F.; REICHERT, J.M. & DIDONÈ, I.A. Caracterização de padrões de chuvas ocorrentes em Santa Maria (RS). R. Bras. Ci. Solo, 25:475-483, 2001.
  • MORAIS, L.F.B.; MUTTI, L.S.M. & ELTZ, F.L.F. Índices de erosividade correlacionados com perdas de solo no Rio Grande do Sul. R. Bras. Ci. Solo, 12:281-284, 1988a.
  • MORAIS, L.F.B.; MUTTI, L.S.M. & ELTZ, F.L.F. Relações entre características físicas da chuva e perdas de solo no Rio Grande do Sul. R. Bras. Ci. Solo, 12:285-288, 1988b.
  • MORAIS, L.F.B.; SILVA, V.; NASCHENVENG, T.M.C.; HARDOIN, P.C.; ALMEIDA, J.E.L.; WEBER, O.L.S.; BOEL, E. & DURIGON, E. Índice EI30 e sua relação com o coeficiente de chuva do sudoeste do Mato Grosso. R. Bras. Ci. Solo, 15:339-344, 1991.
  • MORETI, D.; CARVALHO, M.P.; MANNIGEL, A.R. & MEDEIROS, L.R. Importantes características de chuva para a conservação do solo e da água no município de São Manuel (SP). R. Bras. Ci. Solo, 27:713-725, 2003.
  • OLIVEIRA JÚNIOR, R.C. & MEDINA, B.F. A erosividade das chuvas em Manaus (AM). R. Bras. Ci. Solo, 14:235-239, 1990.
  • REBOLLEDO R., H.H. Manual SAS por computadora: Análisis estadístico de datos experimentales. México, Trillas, 2002. 208p.
  • RONCATO, M.L.; ELTZ, F.L.F.; GRAMINHO, D.H.; STEFANELO, C.; FIGUEIREDO, J.V. & PEDROSO, R.F. Erosividade mensal das chuvas de Santa Maria de abril 1996 a março de 2004. In: REUNIÃO BRASILEIRA DE MANEJO E CONSERVAÇÃO DO SOLO E DA ÁGUA; MANEJO: INTEGRANDO A CIÊNCIA DO SOLO NA PRODUÇÃO DE ALIMENTOS, 15., Santa Maria, 2004. Resumos expandidos. Santa Maria, Universidade Federal de Santa Maria, 2004. CD-ROM
  • ROQUE, C.G.; CARVALHO, M.P. & PRADO, R.M. Fator erosividade da chuva de Piraju (SP): Distribuição, probabilidade de ocorrência, período de retorno e correlação com o coeficiente de chuva. R. Bras. Ci. Solo, 25:147-156, 2001.
  • RUFINO, R.L.; BISCAIA, R.C.M. & MERTEN, G.H. Determinação do potencial erosivo da chuva do Estado do Paraná através de pluviometria: Terceira aproximação. R. Bras. Ci. Solo, 17:439-444, 1993.
  • SILVA, J.R.C. & DIAS, A.S. A erosividade das chuvas em Fortaleza (CE). II - Correlação com o coeficiente de chuva e atualização do fator R no período de 1962 a 2000. R. Bras. Ci. Solo, 27:347-354, 2003.
  • SILVA, M.L.N.; FREITAS, P.L.; BLANCANEAUX, P. & CURI, N. Índices de erosividade das chuvas da região de Goiânia, GO. Pesq. Agropec. Bras., 32:515-521, 1997.
  • WISCHMEIER, W.H. A rainfall erosion index for a universal soil-loss equation. Soil Sci. Soc. Am. Proc., 23:246-249, 1959.
  • WISCHMEIER, W.H. & SMITH, D.D. Predicting rainfall erosion losses: A guide to conservation planning. Washington, USDA-Agricultural Research, 1978. 58p. (Agricultural Handbook, 537)
  • WISCHMEIER, W.H. & SMITH, D.D. Rainfall energy and its relationship to soil loss. Trans. Am. Geophys. Union, 39:285-291, 1958.
  • 1
    Parte da Tese de Mestrado em Ciência do Solo apresentada pelo primeiro autor à Universidade Federal de Santa Maria – UFSM.
  • Datas de Publicação

    • Publicação nesta coleção
      13 Dez 2007
    • Data do Fascículo
      Out 2007

    Histórico

    • Recebido
      Jan 2006
    • Aceito
      Maio 2007
    Sociedade Brasileira de Ciência do Solo Secretaria Executiva , Caixa Postal 231, 36570-000 Viçosa MG Brasil, Tel.: (55 31) 3899 2471 - Viçosa - MG - Brazil
    E-mail: sbcs@ufv.br