Análise da transição entre dias secos e chuvosos por meio da cadeia de Markov de terceira ordem

Analysis of the transition between dry and wet days through third-order Markov chains

a01v41n9

CLIMATOLOGIA

Análise da transição entre dias secos e chuvosos por meio da cadeia de Markov de terceira ordem

Analysis of the transition between dry and wet days through third-order Markov chains

Thadeu Keller Filho^I; Jurandir Zullo Junior^II; Paulo Roberto Schubnell de Rezende Lima^III

^IPontifícia Universidade Católica do Rio de Janeiro, Rua Marquês de São Vicente, nº 225, Sala 210 F, CEP 22453-900 Rio de Janeiro, RJ. E-mail: [e-mail]

^IIUniversidade Estadual de Campinas, Centro de Pesquisas Meteorológicas e Climáticas Aplicadas à Agricultura, Cidade Universitária Zeferino Vaz, CEP 13083-970 Campinas, SP. E-mail: [e-mail]

^IIIEstocástica Consultoria e Serviços Técnicos em Estatística Ltda., SQS 314, Bloco D, apto. 502, CEP 70383-040 Brasília, DF. E-mail: [e-mail]

RESUMO

O objetivo deste trabalho foi verificar se as ocorrências de dias secos e chuvosos são condicionalmente dependentes da seqüência dos três dias secos e chuvosos anteriores, numa zona pluviometricamente homogênea, por meio da cadeia não-homogênea de Markov de terceira ordem. Os resultados mostraram que as probabilidades diárias de transição podem ser adequadamente estimadas, com base em dados agregados bimestralmente, seguidas de interpolação por meio de funções sinusoidais. Além disso, evidenciou-se que, naquela zona, as ocorrências diárias de chuva são condicionalmente dependentes da seqüência de dias secos e chuvosos nos três dias anteriores. A cadeia não-homogênea de Markov de terceira ordem é um importante instrumento para a análise da dependência entre as seqüências de dias secos e chuvosos em determinadas regiões.

Termos para indexação: zonas pluviometricamente homogêneas, ocorrência de chuva, probabilidade de transição.

ABSTRACT

The aim of this work was to verify if the occurrence of dry and wet days are conditionally dependent on the sequences of the dry and wet three preceding days, in a rainfall homogeneous area, using the nonhomogeneous third-order Markov chains. The results showed that daily transition probabilities can be properly estimated from two-month aggregate data, and then adjusted by means of sinusoidal functions. Besides, it was evidenced that everyday rain events in that area are conditionally dependent on the sequences of the dry and wet three days previous to occurrences. The third-order nonhomogeneous Markov chains are an important instrument for the analysis of the dependence between sequences of dry and wet days in certain areas.

Index terms: rainfall homogeneous areas, rainfall occurrence, transition probabilities.

Introdução

A distribuição e a irregularidade da chuva, no tempo e no espaço, são as principais características climáticas de uma região. A precipitação pluvial, assim como outros fatores climáticos, é sujeita a grandes variações, que afetam de modo significativo a produção agrícola.

A produtividade agrícola é fortemente influenciada pela oferta pluviométrica e pela freqüência e duração dos períodos secos. Na maior parte das regiões tropicais, a alternância entre períodos chuvosos e secos apresenta um notável comportamento sazonal. A agricultura, conforme Ayode (1983), é praticada com maior êxito durante a estação chuvosa ou com o uso da irrigação durante a estação seca.

Em algumas regiões, principalmente no Cerrado, a precipitação total do período chuvoso é suficiente para o desenvolvimento da agricultura, porém, é comum a ocorrência de seqüências de dias secos durante a estação chuvosa; são breves períodos conhecidos como veranicos (Souza & Peres, 1998). O estudo da ocorrência de veranicos é de grande importância para a agricultura pois, geralmente, esses breves períodos de estiagem atingem a cultura em sua fase reprodutiva, o que causa redução em sua produtividade (Assad et al., 1993).

A determinação da distribuição de probabilidades da precipitação pluvial e a estimação das probabilidades da transição entre dias secos e chuvosos são de grande importância para a otimização do planejamento das atividades agrícolas, pois contribuem de forma decisiva para a tomada de decisão quanto às épocas de plantio, para a avaliação do risco de deficiência hídrica e para o dimensionamento e manejo dos sistemas irrigados.

O principal instrumento na análise probabilística das transições entre dias secos e chuvosos é a cadeia de Markov (Gordon, 1965), amplamente utilizada no estudo da ocorrência da precipitação pluvial diária, tanto no exterior quanto no Brasil.

Em um trabalho pioneiro, Gabriel & Neumann (1962) utilizaram cadeias homogêneas de Markov para descrever a alta persistência de dias secos consecutivos em Tel Aviv, Israel. Outro trabalho pioneiro foi o de Stern & Coe (1982), que definiram um modelo de simulação para a precipitação pluvial diária, no qual a ocorrência de dias chuvosos é gerada por uma cadeia de Markov, e a quantidade de chuva é gerada por uma distribuição gama. Esse modelo foi aperfeiçoado, utilizando-se séries de Fourier para representar as probabilidades de transição (Stern & Coe, 1984).

A partir desses trabalhos, multiplicaram-se as aplicações de cadeias de Markov na análise estatística da ocorrência de dias secos e úmidos. Entre as inúmeras pesquisas recentes, pode-se citar: Muhammad & Nabi (1991), que analisaram a ocorrência de chuva em quatro localidades no Paquistão; Assis (1991), que modelou a ocorrência de dias secos e úmidos em Piracicaba, SP, e Pelotas, RS; Back (1997), que se baseou nos modelos de Stern & Coe (1982, 1984), para ajustar a precipitação pluvial no posto de Urussanga, SC; Punyawardena & Kulasiri (1997), que compararam modelos de cadeias de Markov de primeira e de segunda ordem para a descrição da ocorrência semanal de chuva, na zona seca de Sri Lanka; De Groen & Savenije (2000), que utilizaram cadeias de Markov para simular a ocorrência de chuva, em um modelo de transpiração mensal, com dados de três localidades no Zimbábue; Andrade Júnior et al. (2001), que simularam a precipitação pluvial diária para Parnaíba e Teresina, PI, com cadeias de Markov para gerar seqüências de dias secos e úmidos.

Na maioria dos estudos se utiliza a cadeia de Markov de primeira ou de segunda ordem. Todavia, o uso de cadeias de Markov de ordem superior a dois é de grande importância, pois, em certas regiões, a ocorrência de chuva em um dia pode ser condicionalmente dependente da seqüência de dias secos e chuvosos em vários dias anteriores. Isto pode acontecer, por exemplo, em regiões onde, em certas épocas, é freqüente a ocorrência de vários dias consecutivos de chuva ou longas séries de dias secos.

O objetivo deste trabalho foi verificar se as ocorrências de dias secos e chuvosos são condicionalmente dependentes da seqüência de dias secos e chuvosos nos três dias anteriores, em uma zona pluviometricamente homogênea, por meio de cadeias não-homogêneas de Markov de terceira ordem.

Material e Métodos

O processo gerador de uma seqüência aleatória de dias secos e chuvosos pode ser descrito por uma cadeia de Markov, com as seguintes características:

1) estados da cadeia: S = dia seco; C = dia chuvoso;

2) indicador do estado no dia t:

3) probabilidades de transição de estado:

em que r (r = 1,2,3,...) é a ordem da cadeia, e J_t é o valor do indicador do estado no dia t.

Assim, P(0/001) representa a probabilidade de que um dia seja seco, na hipótese de que os três dias anteriores tenham sido, seqüencialmente, seco, seco e chuvoso.

Como, nesse modelo, as probabilidades de transição não são constantes, pois variam com o dia t, a cadeia de Markov é não-homogênea.

Uma cadeia de Markov é de ordem r (r = 1,2,3,...), se for:

quaisquer que sejam J_t-r-1, J_t-r, J_t-r+1,..., J_t-1, J_t (r = 1,2,3, ...).

Como os indicadores J_t assumem apenas os valores 0 e 1, a igualdade (3) implica:

Confrontando-se (3) e (4), podemos concluir que uma cadeia de Markov é de ordem r (r = 1,2,3,...) se:

quaisquer que sejam J_t-r-1, J_t-r, J_t-r+1,..., J_t-1 (r = 1,2,3,...).

Convém definir, ainda, cadeias de Markov de ordem zero, correspondentes a seqüências não-correlacionadas de dias secos e úmidos. Uma cadeia de Markov é de ordem zero se:

Para se estimarem as probabilidades de transição, em uma cadeia de Markov de ordem r, são necessários dados referentes ao número de ocorrências de todos os 2^r+1 estados possíveis, em r+1 dias consecutivos. Assim, representando-se por N_n (J_t-r, J_t-r+1,..., J_t) o número de ocorrência da seqüência de estados (J_t-r, J_t-r+1,..., J_t), em n dias consecutivos, as estimativas das probabilidades associadas à cadeia de Markov são dadas por:

Uma dificuldade para a estimação dessas probabilidades é a insuficiência de dados. De fato, o número de observações disponíveis da precipitação pluvial, em determinado dia, coincide com o número de anos abrangidos pela série histórica correspondente. Como o número de estados possíveis das seqüências (J_t-r, J_t-r+1,..., J_t) aumenta exponencialmente com a ordem da cadeia de Markov, os dados disponíveis nem sempre são suficientes, para que as estimativas (7) e (8) sejam efetuadas com precisão, se a ordem da cadeia de Markov é maior que três. Por essa razão, as análises deste trabalho limitam-se a cadeias de Markov definidas até a terceira ordem.

A escassez de dados pode ser contornada mediante a agregação espacial dos dados, agrupados em regiões onde o regime pluviométrico possui comportamento assemelhado. A regionalização explora ao máximo as informações existentes, o que permite melhorar as estimativas baseadas em dados insuficientes (Tucci, 2004).

Neste trabalho, a agregação espacial dos dados pluviométricos teve por base os estudos realizados por Keller Filho (1998) e Lima (1998) e os resultados obtidos por Keller Filho et al. (2005), que identificaram 25 zonas pluviometricamente homogêneas, em todo o território brasileiro. Foi escolhida uma dessas zonas homogêneas, que Keller Filho et al. (2005) denominaram Zona A, que abrange o norte do Tocantins, parte do sul de Goiás e grande parte de Mato Grosso. Essa zona, inserida no bioma Cerrado, apresenta uma precipitação pluvial anual relativamente elevada, com chuvas no verão e estiagens prolongadas no inverno.

Foram utilizados dados da precipitação pluvial diária, observada em cada um dos 82 postos pluviométricos da Zona A, cuja distribuição espacial é mostrada na Figura 1. Os dados abrangem um período de 15 a 30 anos e foram obtidos no banco de dados do Sistema de Monitoramento Agroclimatológico (Agritempo), desenvolvido pela Embrapa Informática Agropecuária e pela Unicamp. No Sistema Agritempo, esses dados foram validados e tiveram sua qualidade e consistência analisadas, segundo metodologia descrita por Lima (2003). A definição de dia seco corresponde àquele em que a precipitação pluvial foi inferior a 0,1 mm.

Analisando-se as proporções observadas de dias secos na Zona A, em cada um dos 365 dias de um ano, verifica-se que elas evoluem de forma bastante regular e apresentam um único pico no intervalo anual, sem oscilações secundárias (Figura 2). Isto sugere que as estimativas das probabilidades de transição entre dias secos e chuvosos, ao longo de um ano, possam ser ajustadas por uma função periódica do tipo sinusoidal (Bloomfield, 1976), de periodicidade igual a 365 dias:

em que a, b e g são constantes; t = dia do ano (t = 1,2,3, ........365) e

A regularidade da evolução das proporções de dias secos e chuvosos permite aumentar a precisão das estimativas das probabilidades de transição, mediante agregação temporal dos dados (Kendall, 1976). Verificou-se que diversas seqüências possíveis dos estados de chuva e seca, em dois ou três dias consecutivos, são relativamente raras em períodos inferiores a dois meses. Por isso, serão analisadas apenas agregações em períodos não inferiores a um bimestre.

O critério de escolha foi o de encontrar um período tão longo quanto possível, dentro do qual as probabilidades de transição não apresentem grandes discrepâncias, a fim de se obter adequado grau de precisão nas estimativas dessas probabilidades. Para esse propósito, todas as probabilidades de transição referentes a cadeias de Markov até a terceira ordem foram estimadas com base em agregações sucessivas em períodos bimestrais, trimestrais, quadrimestrais, semestrais e anuais. Utilizou-se, a seguir, a análise de variância para testar a igualdade das probabilidades de transição, dentro de cada período de agregação. Para identificar, aproximadamente, as variações das probabilidades de transição, dentro dos períodos, cada um deles foi subdividido em vários grupos de meses, formados por aqueles que ocupam a mesma ordem dentro de cada período de agregação (Tabela 1). A análise de variância também permitiu testar, em cada período de agregação, a homogeneidade das estimativas das probabilidades de transição, ao longo dos anos de observação.

Escolhido o período de agregação, as estimativas das probabilidades de transição foram atribuídas às épocas correspondentes ao meio dos respectivos períodos. Embora esses períodos tenham sido escolhidos de tal forma que dentro deles as probabilidades de transição não apresentassem grandes discrepâncias, estabeleceu-se que as probabilidades diárias de transição deveriam apresentar pequenas diferenças ao longo de dias sucessivos. A acumulação inevitável dessas pequenas diferenças, em um período do tempo suficientemente longo, faz com que entre os diversos períodos de agregação as diferenças temporais entre as probabilidades de transição sejam consideráveis. Assim, as estimativas das probabilidades diárias de transição foram obtidas por interpolação de uma função, ajustada às estimativas das probabilidades de transição alocadas ao meio de cada período. Para isso, utilizou-se a função sinusoidal, definida pela equação (9). Como essa equação é uma função linear de seus parâmetros, utilizou-se o critério dos mínimos quadrados (Bloomfield, 1976) para estimar esses parâmetros.

A escassez de dados somente permite utilizar cadeias de Markov até a terceira ordem. A escolha da ordem baseou-se no teste das igualdades (5) e (6), com r = 0,1,2 (Tabela 2). Estruturou-se esses testes de modo que permitissem verificar se as estimativas das probabilidades de transição são significativamente diferentes, entre os períodos anuais abrangidos pelos dados observados. Para efetuar esses testes, utilizou-se a análise de variância, por sua versatilidade na decomposição da variabilidade dos dados em mais de uma fonte de variação.

Resultados e Discussão

Os resultados iniciais do trabalho referem-se à escolha de um período adequado para a agregação dos dados. Mediante o uso dos estimadores (7) e (8), foram estimadas todas as probabilidades de transição referentes a cadeias de Markov até a terceira ordem, correspondentes a agregações em períodos bimestrais, trimestrais, quadrimestrais, semestrais e anuais. Para essas estimativas, na Tabela 3 estão descritos os P-Valores referentes à aplicação da análise de variância, para testar a igualdade das probabilidades de transição dentro dos grupos de meses definidos na Tabela 1, o que corresponde a verificar se as probabilidades de transição não são muito discrepantes, dentro de cada período de agregação.

Examinando-se a Tabela 3, verifica-se que a hipótese de igualdade das probabilidades de transição, dentro de cada período considerado, somente pode ser sustentada, quando a agregação é feita em períodos de dois meses.

A análise de variância permitiu, ainda, que fosse testada a homogeneidade das estimativas de cada probabilidade de transição, ao longo dos anos de observação. Os P-Valores desse teste, aplicado às estimativas obtidas nos diversos períodos de agregação, estão apresentados na Tabela 4. Verifica-se que as estimativas das diversas probabilidades de transição apresentam maior homogeneidade ao longo dos anos de observação, quando os dados são agregados em períodos bimestrais. Pode-se observar que nas agregações bimestrais todos os P-Valores são superiores a 5%, exceto o correspondente à probabilidade de transição de dia chuvoso para seco, que é igual a 4,51%. Por essas razões, decidiu-se agregar os dados com base em períodos de dois meses.

Na Tabela 5, são descritas as estimativas das probabilidades de ocorrência de dia seco e as probabilidades de transição entre dias secos e chuvosos, para cadeias de Markov de ordens um, dois e três, obtidas com os estimadores (7) e (8), com base em dados agregados em períodos de dois meses.

Os valores dessas estimativas foram centrados em cada período bimestral e depois ajustados pela equação (9), utilizando-se o método dos mínimos quadrados. Na Tabela 6 apresentam-se, para cada probabilidade de transição, as estimativas dos parâmetros a, b e g da equação (9), bem como os valores dos respectivos coeficientes de determinação e erros quadráticos médios.

Estimados seus parâmetros, as funções ajustantes permitiram obter aproximações para todas as probabilidades diárias de transição, nas cadeias de Markov de ordem zero, um, dois e três, por meio de interpolação. A representação de probabilidades de transição por funções periódicas, definidas por séries de Fourier, é muito utilizada nos estudos climáticos (Stern & Coe, 1984).

Na Figura 3, está descrita a evolução, ao longo de um período anual, das funções ajustantes das probabilidades de transição para dia seco, nas cadeias de Markov de ordens zero e dois. Analogamente, a Figura 4 descreve a evolução dessas funções ajustantes, quando as cadeias de Markov são de ordem três. Nas Figuras 3 e 4 também estão destacadas as probabilidades de transição que foram atribuídas a cada período bimestral.

Conforme se pode observar nas Figuras 3 e 4, a evolução ao longo de um ano da função ajustante das probabilidades de ocorrência de dia seco e de transição para dia seco, nas cadeias de Markov até a terceira ordem, apresentam um mínimo no mês de janeiro e um máximo no mês de julho. Esse resultado é consistente com o fato de que, na região do Cerrado, existem duas estações bem distintas, em que o verão é essencialmente quente e chuvoso e o inverno é menos quente e seco (Vianello & Alves, 1991). No Planalto Central brasileiro, o clima é caracterizado por uma época de seca de maio até setembro (Wolf, 1977).

A escolha da ordem da cadeia de Markov, mais adequada para a Zona A, foi efetuada utilizando-se a análise de variância, para testar a igualdade das diversas probabilidades de transição para dia seco, segundo os estados possíveis dos dias anteriores, de acordo com o roteiro descrito na Tabela 2. Os resultados dessa análise, para a qual se adotou um nível de significância de 5%, estão apresentados na Tabela 7. Verifica-se que uma cadeia de Markov de ordem zero não é adequada, uma vez que as probabilidades de transição para dia seco diferem significativamente nos casos em que o dia anterior é seco ou é chuvoso.

Verifica-se, ainda, que as cadeias de Markov de primeira ordem também não são adequadas, pois as probabilidades de transição para dia seco diferem significativamente quando se consideram os diversos estados possíveis nos dias anteriores.

A mesma conclusão se impõe para a cadeia de Markov de segunda ordem, tendo em vista que as probabilidades de transição para dia seco diferem significativamente quando se consideram todos os estados possíveis nos três dias anteriores.

A Tabela 7 permite constatar, também, que as estimativas de todas as probabilidades de transição diferem significativamente entre os seis períodos bimestrais definidos para a agregação dos dados.

Esses resultados evidenciam que na Zona A, do Planalto Central brasileiro, a ocorrência de um dia chuvoso é condicionalmente dependente da seqüência de dias secos e chuvosos nos três dias anteriores.

Portanto, as cadeias não-homogêneas de Markov de terceira ordem podem ser utilizadas, com vantagem, na análise probabilística da seqüência de dias secos e chuvosos, em regiões de clima assemelhado ao da Zona A do Planalto Central brasileiro.

Além de sua finalidade descritiva, as cadeias de Markov podem servir de base para a definição de um processo gerador, para simular seqüências de dias secos e chuvosos. O uso da simulação permite estimar distribuições de probabilidades de eventos, para os quais existem poucas observações (Koutsoyiannis & Zarris, 1999). Como exemplo, a simulação de seqüências de dias secos e chuvosos, por meio de cadeias não-homogêneas de Markov, pode ser utilizada para a estimação e a análise das seguintes distribuições de probabilidades: número de estiagens de determinada duração, em certo período do tempo; número e duração de veranicos em determinada estação chuvosa; duração da maior estiagem em determinado período do tempo; duração de dias chuvosos consecutivos em determinada região; determinação do tempo de retorno de estiagens prolongadas.

Conclusões

1. Em uma grande zona pluviometricamente homogênea do Planalto Central brasileiro, a ocorrência de chuva, em cada dia, é condicionalmente dependente da seqüência de dias secos e chuvosos nos três dias anteriores.

2. A escassez de dados para se estimarem as probabilidades de transição, entre dias secos e chuvosos, pode ser contornada por meio da adequada agregação espacial e temporal dos dados.

3. As probabilidades diárias de transição entre dias secos e chuvosos podem ser obtidas por interpolação das probabilidades de transição, estimadas com dados agregados temporalmente, mediante o uso de uma função sinusoidal.

Recebido em 18 de maio de 2005 e aprovado em 9 de junho de 2006

• ANDRADE JÚNIOR, A.S.; FRIZZONE, J.A.; SENTELHAS, P.C. Simulação da precipitação diária para Parnaíba e Teresina, PI, em planilha eletrônica. Revista Brasileira de Engenharia Agrícola e Ambiental, v.5, p.271-278, 2001.
• ASSAD, E.D.; SANO, E.E.J.; MASUTOMO, R.; CASTRO, L.H.; SILVA, F.A.M. Veranicos na região dos cerrados brasileiros: freqüências e probabilidades de ocorrência. Pesquisa Agropecuária Brasileira, v.28, p.993-1003, 1993.
• ASSIS, F.N. Modelagem da ocorrência e da quantidade de chuva e de dias secos em Piracicaba-SP e Pelotas-RS 1991. 134p. Tese (Doutorado) - Escola Superior de Agricultura Luiz de Queiroz, Piracicaba.
• AYODE, J.O. Introdução à climatologia para os trópicos Rio de Janeiro: Bertrand Brasil, 1983. 332p.
• BACK, A.J. Aplicação de cadeias markovianas à modelagem da precipitação diária. Revista de Tecnologia e Ambiente, v.4, p.17-28, 1997.
• BLOOMFIELD, P. Fourier analysis of time series: an introduction. New York: Wiley, 1976. 258p.
• DE GROEN, M.M.; SAVENIJE, H.H.G. Markov chains at daily timesteps, key to transpiration at monthly timesteps. In: WARFSA WATERNET SYMPOSIUM, 1., 2000, Maputo. Sustainable use of water resources: proceedings. Maputo: WARFSA, 2000. p.1-12.
• GABRIEL, K.R.; NEUMANN, J. A Markov chain model for daily rainfall occurrence at Tel Aviv. Quarterly Journal of the Royal Meteorological Society, v.88, p.90-95, 1962.
• GORDON, P. Théorie des chaînes de Markov finies et ses applications Paris: Dunod, 1965. 130p.
• KELLER FILHO, T. Redução do risco climático na agricultura - uma abordagem probabilística. 1. Metodologia. In: SIMPÓSIO INTERNACIONAL DE SECURIDADE E ZONEAMENTO AGRÍCOLA DO MERCOSUL, 1., 1998, Brasília. Anais Brasília: Ministério da Agricultura e Abastecimento, 1998. p.40-47.
• KELLER FILHO, T.; ASSAD, E.D.; LIMA, P.R.S. de R. Regiões pluviometricamente homogêneas no Brasil. Pesquisa Agropecuária Brasileira, v.40, p.311-322, 2005.
• KENDALL, M. Time series 2^nd ed. London: Griffin, 1976. 197p.
• KOUTSOYIANNIS, D.; ZARRIS, D. Simulation of rainfall events for design purposes with inadequate data. In: GENERAL ASSEMBLY OF THE EUROPEAN GEOPHYSICAL SOCIETY, 24., 1999, The Hague. Proceedings The Hague: European Geophysical Society, 1999. p.35-44.
• LIMA, J.G.S. Gerenciamento de dados climatológicos heterogêneos em agricultura 2003. 103p. Dissertação (Mestrado) - Universidade Estadual de Campinas, Campinas.
• LIMA, P.R.S. de R. Redução do risco climático na agricultura - uma abordagem probabilística - Resultados obtidos. In: SIMPÓSIO INTERNACIONAL DE SECURIDADE E ZONEAMENTO AGRÍCOLA DO MERCOSUL, 1., 1998, Brasília. Anais Brasília: Ministério da Agricultura e Abastecimento, 1998. p.49-57.
• MUHAMMAD, F.; NABI, G. Markov chain model for rainfall occurrence in Pakistan. Journal of Islamic Academy of Sciences, v.4, p.173-178, 1991.
• PUNYAWARDENA, B.V.R.; KULASIRI, D. On development and comparative study of two Markov models of rainfall in the dry zone of Sri Lanka. In: JOINT INTERNATIONAL CONFERENCE ON AGRICULTURAL ENGINEERING AND TECHNOLOGY, 1997, Dhaka. Proceedings Dhaka: ASAE, 1997. p.231-238.
• SOUZA, S.A.V. de; PERES, F.C. Programa computacional para simulação da ocorrência de veranicos e queda de rendimento. Pesquisa Agropecuária Brasileira, v.33, p.1951-1956, 1998.
• STERN, R.D.; COE, R. A model fitting analysis of daily rainfall data. Journal of the Royal Statistical Society, v.147, p.1-34, 1984.
• STERN, R.D.; COE, R. The use of rainfall models in agricultural planning. Agricultural Meteorology, v.26, p.35-50, 1982.
• TUCCI, C.E.M. Hidrologia: ciência e aplicação. Porto Alegre: Ed. da UFRGS, 2004. 943p.
• VIANELLO, R.L.; ALVES, A.R. Meteorologia básica e aplicações Viçosa: Imprensa Universitária, Universidade Federal de Viçosa, 1991. 449p.
• WOLF, J.M. Probabilidade de ocorrência de períodos secos na estação chuvosa de Brasília. Pesquisa Agropecuária Brasileira, v.12, p.141-150, 1977.

Datas de Publicação

• Publicação nesta coleção
  25 Out 2006
• Data do Fascículo
  Set 2006

Histórico

• Recebido
  18 Maio 2005
• Aceito
  09 Jun 2006