Estudo da erosão na microbacia do Ceveiro (Piracicaba, SP): I - Estimativa das taxas de perda de solo e estudo de sensibilidade dos fatores do modelo EUPS

A erosão acelerada do solo constitui um sério problema global de degradação de terras agrícolas. Por ser a pesquisa de erosão comumente cara e morosa, a utilização de modelos permite estimar a perda de solo em localidades e condições não diretamente representadas nos estudos. O objetivo deste trabalho foi estimar a expectativa de erosão do solo em uma área intensivamente cultivada com cana-de-açúcar no município de Piracicaba (SP), como subsídio ao diagnóstico do uso atual das terras visando ao controle do processo. O trabalho foi realizado em ambiente de sistema de informação geográfica aplicando a equação universal de perda de solo (EUPS). Especificamente pretendeu-se estudar a influência dos fatores K (erodibilidade do solo), LS (topográfico), C (uso e manejo da cultura) e P (práticas de controle da erosão) nas estimativas de perda de solo. Os resultados obtidos mostraram que o processo de degradação dos solos por erosão acelerada é intenso em cerca de dois terços da área avaliada, ocupados principalmente com cana-de-açúcar. Nessas áreas, a perda média estimada de 58 Mg ha-1 ano-1 equivale a cerca de quatro vezes o limite superior da tolerância média indicada. Nas áreas remanescentes, ocupadas com mata, vegetação ciliar, reflorestamento e pastagens, as taxas estimadas de perda de solo são baixas, em torno de 2 Mg ha-1 ano-1 , abaixo dos valores de tolerância. O estudo de sensibilidade dos fatores do modelo revelou que, para uma mesma condição climática (fator R constante), os fatores C e P (cobertura e manejo da cultura e práticas de controle da erosão) definiram a ordem de grandeza das perdas, em unidades, dezenas ou centenas de Mg ha-1 ano-1 . Para uma mesma categoria de uso da terra (fatores C e P constantes), o fator LS foi o que explicou a maior parte da variação observada. Os resultados obtidos a partir de uma simulação mostraram que a influência do fator K (erodibilidade do solo) nas estimativas foi muito menor do que a do fator topográfico. A aplicação da EUPS permitiu estimar a expectativa de erosão do solo na área de estudo. O uso de geotecnologias e de métodos geoestatísticos de análise mostrou-se uma abordagem promissora nos estudos de erosão.

Estimativa da susceptibilidade à compactação e do suporte de carga do solo com base em propriedades físicas de solos do Rio Grande do Sul

O conhecimento das relações entre propriedades físicas e mecânicas do solo pode contribuir no desenvolvimento de funções de pedotransferência, que permitam estimar outras propriedades do solo de difícil mensuração. Neste trabalho, objetivou-se avaliar a relação entre a susceptibilidade à compactação e o suporte de carga com propriedades físicas de solos do sul do Brasil. Foram avaliadas a resistência à penetração, a umidade, a densidade e a compressibilidade de seis solos. A resistência à penetração pode ser estimada pelo modelo que considera a umidade e densidade do solo. Solos com maior densidade inicial apresentaram menor susceptibilidade à compactação e menor deformação, quando submetidos a pressões externas. Quanto maior a resistência do solo à penetração, menor é a deformação e maior é a capacidade de suporte de carga, embora isso não indique solos com qualidade física adequada para as culturas; quanto maior a deformação do solo, maior a susceptibilidade à compactação e menor a capacidade de suporte de carga. A susceptibilidade de um solo à compactação e sua capacidade de suporte de carga podem ser estimadas, respectivamente, pela densidade inicial e pela resistência do solo à penetração.

Estimativa da resistência mecânica à penetração de um latossolo vermelho sob cultura perene no norte do estado do Paraná

Uma das propriedades físicas importantes para o manejo e estudo da qualidade física dos solos é sua resistência à penetração, uma vez que esta propriedade se apresenta relacionada a diversos atributos do solo indicadores do grau de compactação. O objetivo deste trabalho foi estimar a resistência mecânica do solo à penetração de um Latossolo Vermelho eutroférrico sob cultura perene, em três condições de umidade gravimétrica, a partir da densidade, umidade, teor de argila (Arg) e matéria orgânica do solo. A resistência mecânica do solo à penetração foi avaliada mediante o uso do penetrômetro de impacto, e os dados obtidos em impactos dm-1 foram convertidos em resistência dinâmica (MPa). No momento da determinação da resistência mecânica do solo, foram coletadas 30 amostras por profundidade: 0-0,10; 0,1-0,3; 0,3-0,5 m, para determinação de três diferentes teores de umidade gravimétrica (Ug) (Ug1 < 0,20 kg kg-1 , 0,20 < Ug2 < 0,30 kg kg-1 e Ug3 > 0,30 kg kg-1 ), composição granulométrica, teor de matéria orgânica (MO) e densidade do solo (Ds). Foram aplicados testes de distribuição normal para seleção das variáveis participantes dos modelos, por meio do teste W. Em seguida, foram realizadas análises de correlação simples entre as variáveis e análise de regressão múltipla por Stepwise. Para obtenção dos modelos, realizou-se uma análise com quatro variáveis independentes, representadas por: X1 = Ug, X2 = Ds, X3 = MO, X4 = Arg, e a variável dependente (Y), caracterizada pela resistência do solo à penetração (RP). Após obtenção dos modelos finais, testaram-se os efeitos de multicolinearidade por meio do VIF (Fator de Incremento de Variância), sendo os valores de VIF menores que 10 indicativos de que o efeito da multicolineridade não influencia os resultados da regressão. Foi utilizado o nível de significância do valor F da variável de 15 % de probabilidade para a inclusão e exclusão de variáveis no modelo. Os resultados permitiram concluir que as equações que melhor estimaram a resistência do solo à penetração foram: para Ug < 0,20 kg kg-1 : RPE = 1,12 + 2,04Ds + 0,17MO (R² = 0,92); para (0,20 < Ug < 0,30 kg kg-1): RPE = 2,93 - 17,85Ug + 2,98Ds (R² = 0,95); para Ug > 0,30 kg kg-1: RPE = 3,20 - 13,78Ug + 4,43Ds (R² = 0,97). Em todas as situações estudadas, os coeficientes de variação variaram de altos (situação UG1) a muito altos (situação UG2 e Ug3), caracterizando alta variabilidade espacial da RP; logo, a utilização dos modelos selecionados deve ser feita com cautela, pois ocorre importante dispersão de pontos estimados por meio dos modelos, em relação aos valores medidos.

Estimativa do escoamento superficial em uma bacia hidrográfica com base em modelagem dinâmica e distribuída

A estimativa do escoamento superficial em bacias hidrográficas é de suma importância para conservação dos recursos naturais; entretanto, esse é um processo complexo e dinâmico, principalmente no contexto de sua variabilidade espacial. Dessa forma, torna-se adequada a aplicação dos Sistemas de Informações Geográficas (SIG) usando pequenas células de informação, pois assim é possível considerar o comportamento espacial das variáveis associadas à origem do escoamento superficial. Este trabalho teve como objetivo implementar os modelos hidrológicos Curva Número (CN-SCS) e Curva-Número Modificado (CN-MMS), com base na linguagem de programação do SIG PCRaster e em uma base de dados reduzida, de forma distribuída e dinâmica, com o intuito de estimar as lâminas de escoamento superficial geradas numa bacia hidrográfica de Latossolos, localizada no município de Nazareno, região dos Campos das Vertentes, Minas Gerais. Para aplicação do modelo CN-SCS foi preciso desenvolver um mapa com valores de CN no formato do PCRaster, enquanto para o modelo CN-MMS foram necessários os seguintes mapas: umidade volumétrica de saturação do solo, umidade volumétrica inicial do solo e profundidade de solo. Para simulação e avaliação de ambos os modelos, foram aplicados 18 eventos de chuva natural que provocaram escoamento superficial, durante o ano hidrológico 2004-2005, e suas respectivas lâminas de escoamento observadas. A análise do desempenho dos modelos foi feita aplicando-se análise de sensibilidade baseada no erro médio e na Raiz do Erro Quadrático (REQ). Tendo-se como referência essas estatísticas de precisão, pôde-se constatar que o modelo CN-MMS apresentou melhor calibração quando comparado ao modelo CN-SCS, devido à consideração direta da umidade inicial do solo. Contudo, a estruturação dos modelos no SIG PCRaster possibilitou o desenvolvimento de uma ferramenta computacional eficaz e útil para simulação do escoamento superficial, visto que propicia estruturação de rotinas computacionais considerando os problemas associados à variabilidade espacial dos dados de entrada dos modelos.

Teores de nitrogênio mineral do solo para predição do potencial produtivo de cevada

O teor de N mineral do solo não é usado na recomendação de adubação nitrogenada para cultivos de grãos no Brasil, devido à sua elevada instabilidade por efeito dos fatores ambientais. Seu uso poderia melhorar a eficiência das adubações e reduzir os custos de produção. Foi feita uma análise dos teores de NH4+ e de NO3- do solo, até 15 cm de profundidade, no período da pré-semeadura e do rendimento de grãos em 61 ensaios de cevada, em 20 municípios, em três Regiões Fisiográficas do Rio Grande do Sul, de 2002 a 2005. Os rendimentos foram divididos em classes (muito baixo, baixo, médio, alto e muito alto), que foram relacionadas aos teores de N mineral, usando o modelo de regressão de Gompertz. O ajuste indicou dois marcos de N mineral delimitadores do rendimento relativo: 8 e 18 mg kg-1 de N no solo. Abaixo de 8 mg kg-1 de N houve predomínio de muito baixos rendimentos, evidenciando o efeito da baixa disponibilidade sobre a produtividade. Acima de 18 mg kg-1 de N, este nutriente não foi mais fator decisivo na formação do potencial produtivo, sendo indicativo de potencial produtivo muito alto. Esses valores poderão futuramente ser incorporados como base de um modelo de recomendação de adubação nitrogenada.

Disponibilidade de potássio para as plantas em solos do sul do Brasil estimada por métodos multielementares

A avaliação da disponibilidade de K para as plantas é feita pela estimativa de seus teores na forma trocável. Métodos multielementares aumentam a eficiência nos laboratórios de análises de solo, porém a capacidade de troca de cátions (CTC) do solo pode alterar a predição do K disponível. O objetivo deste trabalho foi estimar o K disponível para o milho e para a soja em solos representativos do Estado do RS, pelos métodos Mehlich-1 (M1), Mehlich-3 (M3), resina de troca iônica (resina) e Acetato de Amônio pH 7,0 (AcNH4), bem como verificar se a CTC influencia a estimativa da disponibilidade de K no solo para as plantas. O experimento foi conduzido em microparcelas a céu aberto e, após 45 dias de cultivo, quantificou-se a matéria seca da parte aérea das plantas e determinou-se o o teor de K no tecido. Foram feitas as correlações entre os teores de K extraídos pelos métodos e a relação com a absorção de K pelas plantas em todos os solos e separando os solos em classes de CTC. As quantidades de K extraídas pela resina foram aproximadamente 50 % menores que as das soluções de M1, M3 e AcNH4. Os coeficientes de correlação entre os teores de K extraído pelos métodos tiveram alto grau de associação. A relação entre os teores de K acumuladas pelas culturas e os teores extraídos pelos métodos não diferiram, independentemente da separação por classes de CTC. A disponibilidade de K para as plantas em solos do RS pode ser feita pelos métodos testados no presente trabalho. Não houve diferença entre os métodos na avaliação da disponibilidade de K nas classes de CTC. O método da resina, embora não diferindo dos demais na predição do teor de K disponível, apresenta desvantagem em relação aos demais pela menor capacidade extrativa, o que resulta em classes de disponibilidade menores e propicia o erro analítico.

Regressões Logísticas Múltiplas: fatores que influenciam sua aplicação na predição de classes de solos

Métodos mais eficazes para determinação do padrão de distribuição de classes de solo na paisagem precisam ser avaliados visando suprir a demanda por mapas de solo em escalas regional e global. Neste estudo, Regressões Logísticas Múltiplas foram utilizadas como modelos preditores em uma aplicação de Mapeamento Digital de Solos. Os modelos foram gerados utilizando um mapa de solos existente como variável dependente e atributos de terreno como variáveis independentes, o que possibilitou determinar a probabilidade de encontrar classes de solo na paisagem no primeiro e no segundo nível categórico do SiBCS. A qualidade dos mapas preditos foi verificada por meio da matriz de contingência. A classe dos Argissolos foi predita corretamente, em relação ao mapa original, em aproximadamente 85 %. As classes de solos hidromórficos (Planossolos e Gleissolos) foram preditas corretamente em 75 %. Houve confundimento dos modelos para as classes que ocupam posições muito semelhantes na paisagem. Foi verificado também que classes de solo pouco representativas na paisagem não são adequadamente espacializadas em razão da sensibilidade dos modelos logísticos à proporção relativa das amostras usadas para treinar os modelos.

Influência da densidade do solo na estimativa da umidade em um nitossolo vermelho distroférrico, por meio da técnica de TDR

Entre as técnicas para avaliação da umidade do solo, destaca-se a que utiliza a Reflectometria no Domínio do Tempo (TDR), a qual consiste em um método indireto, não destrutivo, rápido, exato e pouco dependente dos fatores ambientais. Porém, trabalhos desenvolvidos após as suas aplicações sugerem a influência de fatores ligados às características mineralógicas e estruturais dos solos, o que torna necessária a calibração local para o adequado funcionamento da TDR. O objetivo deste trabalho foi estudar a influência da densidade do solo (Ds) na medida da umidade volumétrica do solo (θ), utilizando-se a técnica de TDR, em um Nitossolo Vermelho distroférrico. Um modelo que expressa a relação entre a constante dielétrica do solo (Ka) e θ foi ajustado por meio de regressão polinomial múltipla, com a incorporação da Ds para a estimativa dos seus coeficientes. Os resultados permitiram concluir que, nesse solo, o modelo de calibração deve levar em conta o efeito da Ds, sem o que a qualidade das estimativas de umidade, a partir das medidas de constante dielétrica, fica comprometida. Os resultados evidenciaram também que a sensibilidade da TDR à variação da umidade do solo é elevada, tornando possível que a técnica seja utilizada com exatidão em ações que envolvam o monitoramento da variação da umidade do solo. Os resultados indicaram a necessidade de mais estudos sobre a utilização da TDR em solos tropicais.

Estimativa da capacidade de campo pela curva de retenção e pela densidade de fluxo da água

A capacidade de campo é um parâmetro de inegável relevância para o manejo adequado da irrigação. A partir da determinação confiável do conteúdo de água no solo na capacidade de campo, pode-se otimizar a produtividade das culturas agrícolas, maximizando a eficiência do uso da água pelas plantas e evitando a contaminação do lençol freático por lixiviação de fertilizantes e agroquímicos. Nesse sentido, o trabalho teve como objetivo avaliar a capacidade de campo (a) pelo cálculo da densidade de fluxo da água durante o processo de redistribuição no experimento clássico de determinação desse parâmetro e (b) a partir de uma dada tensão da água na curva de retenção, em um Latossolo Vermelho-Amarelo textura média, localizado em área experimental da Escola Superior de Agricultura Luiz de Queiroz, município de Piracicaba, Estado de São Paulo (Brasil). Para isso, instalaram-se no solo, até a profundidade de 1,9 m, 60 tubos de acesso a uma sonda de nêutrons, numa rede de malha quadrada de 5 m (6 x 10 pontos amostrais), e, adjacente a cada tubo, dois tensiômetros, um a 0,75 m e outro a 0,85 m de profundidade. Nos 60 pontos (locais) amostrais, desenvolveu-se um experimento similar ao experimento clássico de determinação da capacidade de campo, no qual o solo foi devidamente submetido a uma lâmina de infiltração visando à saturação do seu perfil; sua superfície, coberta com lona plástica para evitar fluxo de água através dela; e a redistribuição da água, monitorada até a profundidade de 0,8 m durante 20 dias. Observou-se que o conteúdo de água correspondente a uma densidade de fluxo de 1,0 mm dia-1 é a melhor estimativa da capacidade de campo para esse solo, uma vez que as densidades de fluxo de 0,1 e 0,01 mm dia-1, também recomendadas para estimar a capacidade de campo, são muito baixas, a ponto de não terem sido alcançadas neste estudo. Quanto aos resultados obtidos pelo método baseado na curva de retenção, utilizando-se o conteúdo de água correspondente à tensão de 10 kPa em curvas de retenção elaboradas no campo e no laboratório também para os 60 pontos, na profundidade de 0,8 m, observou-se que os valores obtidos a partir das curvas de retenção elaboradas no campo e no laboratório subestimaram e superestimaram, respectivamente, aqueles baseados na densidade de fluxo da água de 1,0 mm dia-1, e a medida de campo foi mais confiável.

Análise multivariada e variabilidade espacial na estimativa da erodibilidade de um argissolo vermelho-amarelo

A erodibilidade é um fator de extrema importância na caracterização da perda de solo, representando os processos que regulam a infiltração de água e sua resistência à desagregação e o transporte de partículas. Assim, por meio da análise de dependência espacial dos componentes principais da erodibilidade (fator K), objetivou-se estimar a erodibilidade do solo em uma área de nascentes da microbacia do Córrego do Tijuco, Monte Alto-SP, e analisar a variabilidade espacial das variáveis granulométricas do solo ao longo do relevo. A erodibilidade média da área foi considerada alta, e a análise de agrupamento k-means apontou para uma formação de cinco grupos: no primeiro, os altos teores de areia grossa (AG) e média (AM) condicionaram sua distribuição nas áreas planas; o segundo, caracterizado pelo alto teor de areia fina (AF), distribui-se nos declives mais convexos; o terceiro, com altos teores de silte e areia muito fina (AMF), concentrou-se nos maiores declives e concavidades; o quarto, com maior teor de argila, seguiu as zonas de escoamento de água; e o quinto, com alto teor de matéria orgânica (MO) e areia grossa (AG), distribui-se nas proximidades da zona urbana. A análise de componentes principais (ACP) mostrou quatro componentes com 87,4 % das informações, sendo o primeiro componente principal (CP1) discriminado pelo transporte seletivo de partículas principalmente em zonas pontuais de maior declividade e acúmulo de sedimentos; o segundo (CP2), discriminado pela baixa coesão entre as partículas, mostra acúmulo da areia fina nas áreas de menor cota em toda a área de concentração de água; o terceiro (CP3), discriminado pela maior agregação do solo, concentra-se principalmente nas bases de grandes declives; e o quarto (CP4), discriminado pela areia muito fina, distribui-se ao longo das declividades nas maiores altitudes. Os resultados sugerem o comportamento granulométrico do solo, que se mostra suscetível ao processo erosivo devido às condições texturais superficiais e à movimentação do relevo.

Verossimilhança na seleção de modelos para predição espacial

Na modelagem estatística da variabilidade espacial, estimam-se os parâmetros da dependência espacial, que são utilizados na interpolação de valores em locais não amostrados. Para tal, o processo de modelagem deve ser realizado com critérios estatísticos que garantam predições confiáveis e representem a real variabilidade local. Este trabalho avaliou diferentes formulações do modelo geoestatístico gaussiano para reconstituir a superfície que representa o fósforo (P) na área, a partir de medições dos teores de P em 48 parcelas experimentais localizadas em Xanxerê, SC, destacando o método utilizado nas análises. A combinação da presença de covariáveis no modelo e a necessidade de transformação para normalidade dos dados definiram quatro alternativas para modelagem. Utilizou-se a função de correlação de Matèrn, avaliada nos valores 0,5; 1,5; e 2,5 para parâmetro de suavidade. Os modelos foram comparados pelo valor maximizado do logaritmo da função de verossimilhança e também por validação cruzada. O modelo selecionado foi o que incorporou a variável resposta transformada, as coordenadas da área como covariáveis e o valor 0,5 para o parâmetro de suavidade. As medidas de validação cruzada pouco acrescentaram aos resultados de comparação por verossimilhança, que evidenciaram que na modelagem geoestatística, o cuidado com observações globais ou locais atípicas, além da seleção com base em diferentes modelos, deve ser o foco para obter resultados compatíveis com a realidade.

Modelo para estimativa da infiltração de água e perfil de umidade do solo

A adequada modelagem da infiltração de água no solo é fundamental para estimação do movimento de água, erosão hídrica, recarga e contaminação de aquíferos. Este trabalho apresenta um modelo para estimativa da infiltração de água no solo (GAML-c), com base no modelo de Green-Ampt-Mein-Larson, que provê descrição da geometria e do deslocamento da frente de umedecimento no solo. Testes experimentais foram conduzidos em Latossolo Vermelho-Amarelo para avaliar o GAML-c, usando-se quatro diferentes cenários: considerando a condutividade hidráulica do solo igual à taxa de infiltração estável (Tie) e a umidade máxima do solo igual ao teor de água na zona de transmissão (θw) (TW); condutividade hidráulica do solo igual à do solo saturado (K0) e a umidade máxima do solo igual θw (KW); condutividade hidráulica do solo igual à Tie e a umidade máxima do solo igual ao teor de água na saturação (θs) (TS); e condutividade hidráulica do solo igual a K0 e a umidade máxima do solo igual θs (KS). Verificou-se que o GAML-c no cenário TW foi o melhor estimador do perfil de umidade do solo, resultando em aceitáveis estimativas da infiltração de água.

Erodibilidade, fatores cobertura e manejo e práticas conservacionistas em argissolo vermelho-amarelo, sob condições de chuva natural

Os modelos de predição de perda de solo têm se tornado importantes ferramentas no planejamento conservacionista, sendo, para tanto, fundamental a estimativa local dos parâmetros que influenciam o processo erosivo. O objetivo deste trabalho foi determinar a erodibilidade (fator K) e os fatores manejo e cobertura (fator C) e práticas conservacionistas de suporte (fator P) da Equação Universal da Perda de Solo (USLE), em um Argissolo Vermelho-Amarelo. Foram avaliadas as perdas de solo (PS) ocorridas em parcelas-padrão de Wischmeier, no período de 2006 a 2011, com seus respectivos índices de erosividade (EI30). Valores de erodibilidade foram calculados pelo quociente entre PS e EI30 das chuvas anuais (Ka) e total no período (Kt) e pelo coeficiente angular gerado pela análise de regressão linear entre esses mesmos parâmetros (Kci e Kct). Os fatores C e P foram estimados para a cultura do milho, durante três anos de cultivo em nível (MN). Os valores de Kt e Kct obtidos foram de 0,0090 e 0,0106 Mg ha h ha-1 MJ-1 mm-1, respectivamente. Os valores médios calculados para os fatores C e P são de 0,0070 Mg ha Mg-1 ha-1 e de 0,45, respectivamente, indicando redução de 55 % na perda de solo do MN, em relação ao MMA. Os valores dos fatores K, C e P encontrados podem ser usados como primeira aproximação para estudos de manejo e conservação do solo e da água na região.

Estimativa da infiltração de água no solo na escala de bacia hidrográfica

A variabilidade espacial das características do solo torna a estimativa da infiltração de água altamente dependente da escala de análise. Muitos trabalhos necessitam de informações da dinâmica da infiltração para descrever o mecanismo de outros processos associados à água, aos sedimentos ou aos solutos. O objetivo deste trabalho foi utilizar duas estratégias metodológicas para a estimativa da infiltração na escala de bacia. Os parâmetros de dois modelos de infiltração foram ajustados a um conjunto de dados coletados durante dois anos de monitoramento, em 77 eventos de chuva-vazão ocorridos sob diferentes condições de uso e manejo do solo numa bacia hidrográfica rural de 1,2 km². Os resultados dos ajustes para os dois modelos de infiltração na escala de bacia foram comparados com aqueles obtidos para dois métodos pontuais de estimativa da infiltração. Com o conjunto de dados do monitoramento, foi possível definir um modelo de infiltração para a bacia monitorada. Além disso, os resultados indicaram que a comparação realizada entre os valores dos modelos com os pontuais são fortemente dependentes da escala de análise; isto é, as medidas pontuais revelaram valores médios de taxa de infiltração alta para a bacia em estudo. Todavia, observações no campo indicaram baixa infiltração, principalmente nas áreas saturadas onde condicionam o escoamento superficial nessa bacia, ao contrário dos resultados pelos modelos de estimativa, que apresentaram menor ordem de magnitude dos valores estimados de infiltração, em razão, entre outros, da maior incorporação dos fatores controladores do processo em seu método.

ESTIMATIVA DA EROSIVIDADE DE CHUVAS NO ESTADO DO PARANÁ PELO MÉTODO DA PLUVIOMETRIA: ATUALIZAÇÃO COM DADOS DE 1986 A 2008

Três estimativas de erosividade de chuvas no Estado do Paraná foram feitas nas décadas de 1980-90, mas nenhuma outra depois. Os objetivos deste trabalho foram atualizar, numa quarta estimativa, as informações de erosividade anual das chuvas para o Paraná, com dados pluviométricos de 1986 a 2008; e gerar informações de erosividade mensal para identificar a distribuição temporal da erosividade. Foram utilizados os métodos pluviográfico e pluviométrico de determinação de erosividade. O primeiro serviu para revalidar equações prévias de correlação entre erosividade por pluviografia e por pluviometria, isso feito em três localidades; já o segundo, para calcular o coeficiente de chuva em 114 localidades, utilizando-se dados de precipitação média mensal e anual, e então a erosividade, com as equações de correlação revalidadas. Mapas de erosividade anual e mensal foram gerados por interpolações geoestatísticas. A erosividade variou entre 5.449 e 12.581 MJ mm ha-1 h-1 ano-1 e aumentou no sentido oeste, a partir da região metropolitana de Curitiba; e no sentido sul, a partir do noroeste paranaense. Isso resultou em maiores erosividades nas regiões centro-sul (Guarapuava), centro-ocidental (Campo Mourão), oeste (Cascavel) e, principalmente, sudoeste (Pato Branco). Também houve incremento de erosividade a partir da região metropolitana de Curitiba, no sentido leste, indicando o litoral como outra região crítica em termos de erosividade. A maior erosividade mensal ocorreu em janeiro (média estadual de 1394 MJ mm ha-1 h-1 mês-1), seguida por tendência de diminuição até agosto (331 MJ mm ha-1 h-1 mês-1), quando então tornou a aumentar. Porém, ocorrem dois importantes repiques de erosividade: um em maio e outro em outubro (843 e 1.173 MJ mm ha-1 h-1 mês-1, respectivamente), principalmente nas regiões sudoeste e oeste. Esses repiques coincidem com períodos de implantação de culturas de inverno ou de verão nessas regiões. Em relação à última estimativa de erosividade anual, feita em 1993, houve aumento da área delimitada pela linha isoerodente de 11.000 MJ mm ha-1 h-1 ano-1 nas regiões sudoeste e oeste e um deslocamento das linhas isoerodentes de 7.000 e 8.000 MJ mm ha-1 h-1 ano-1 para o sentido leste nas regiões norte pioneiro, centro-oriental e sudeste, sugerindo elevação da erosividade nessas regiões. Tais alterações de erosividade podem estar associadas às eventuais mudanças no clima da região, mas também ao simples fato de este estudo estar considerando maior quantidade de dados em relação aos estudos anteriores. Concluiu-se que as equações de correlação determinadas entre os métodos pluviográfico e pluviométrico, obtidas para o Paraná em trabalhos anteriores, continuam válidas; os meses críticos quanto à erosividade são janeiro, maio e outubro; e as menores erosividades ocorrem nas regiões noroeste e metropolitana do Estado.