Manejo da irrigação em pastagem irrigada por pivô-central

Irrigated pasture: water management under center pivot irrigation

MANEJO DE ÁGUA E SOLO

Manejo da irrigação em pastagem irrigada por pivô-central

Irrigated pasture: water management under center pivot irrigation

Alexandre C. Xavier^I; Rubens D. Coelho^II; Leandro F. Lourenço^III; Ronalton E. Machado^IV

^IUFES/CCA/DER. CP 16, CEP 29500-000, Alegre, ES. Fone (28) 3552-1400. E-mail: [e-mail]

^IIUSP/ESALQ, Depto. de Engenharia Rural, CP 9 - CEP:13418-900 - Piracicaba, SP. Fone: (19) 3429-4217 ramal 223 E-mail: [e-mail]

^IIIBolsista CAPES. Fone: (19) 3429-4217 ramal 223. E-mail: [e-mail]

^IVUSP/ESALQ. Fone (19) 3429-4143 ramal 208. E-mail: [e-mail]

RESUMO

A aplicação de lâminas de irrigação em pastagem irrigada sob pivô-central é, de maneira geral, realizada sem um critério técnico pertinente ao sistema, pois se deve considerar que para um mesmo período a pastagem se encontra em diferentes estádios de desenvolvimento em cada parcela, apresentando taxas de evapotranspiração diferenciadas dentro da área irrigada; todavia, usualmente se aplica uma única lâmina para toda a área. Neste trabalho foi desenvolvido um modelo para aplicação de lâminas de irrigação distinta para cada parcela do pivô o qual, de modo geral, considera: i) a capacidade do pivô-central de aplicar lâminas distintas na área; ii) o nível de desenvolvimento da cultura em cada parcela; iii) o período de retorno do gado a determinada parcela (ciclo de pastejo); e iv) o potencial de desenvolvimento da pastagem de certa região. Para modelar o coeficiente de cultura (Kc) foram utilizadas duas metodologias, a primeira com taxa de variação do Kc constante com o número de dias em que a parcela está em descanso (k), e a segunda, com taxa de variação do Kc na forma senoidal com k. O modelo foi aplicado para pastagens hipotéticas nas regiões de Piracicaba e Pereira Barreto, para avaliação e, como resultado, observou-se que o modelo se mostrou sensível ao nível de desenvolvimento de cada parcela e às condições de variação do clima de cada região.

Palavras-chave: capim Tanzânia, manejo da água, coeficiente de cultura

ABSTRACT

The application of irrigation depths in irrigated pasture under center pivot machines, in a general way, is accomplished without a pertinent technical criterion, because it should be considered that for any time period, the pasture plots are at different development stages (rotary pasture), presenting different evapotranspiration rates inside the irrigated area. Furthermore, farmers usually apply a single irrigation depth for the whole area. In this study a model was developed for the application of different irrigation depths in each portion of the pivot (pizza shape plots). The model proposed in this study considers: i) the capacity of the pivot-center to apply different water depths in the area; ii) development stage of crop in each plot iii) cattle return period for grazing in-determined plot (pasture cycle fixed); and iv) potential of pasture development during the agronomic year. To compute the crop coefficient (Kc) variable in the model, two methodologies were used: a) constant variation rate of Kc, according to the resting days the portion is submitted and b) sinusoidal variation rate of Kc, according to the resting period. The model was applied for hypothetical pastures conditions under rural areas of Piracicaba and Pereira Barreto municipalities, in the State of São Paulo, Brazil. The model's sensitivity for the irrigation depths applied was observed for the different stages of Tanzania grass development under the center-pivot, according to the variable climatic conditions of each region.

Key words: Tanzânia grass, irrigation management, crop coefficient

INTRODUÇÃO

Pecuaristas de várias regiões do Brasil vêm adotando o sistema de irrigação por aspersão tipo pivô-central, em pastagens formadas por gramíneas do gênero Panicum (Capim Tanzânia e Monbaça) e Brachiaria (Capim Braquiária), aliado ao uso intensivo de adubações de reposição (N, P e K). Estas áreas são manejadas de forma rotacionada, com ciclo variando entre 28 a 36 dias, sendo o consumo do capim diretamente por pastejo ou através do armazenamento em silagem, utilizada para o confinamento. Todavia, só a irrigação e nutrientes não são condições suficientes para o aumento significativo da produção mas, sim, a associação destas com condições climatológicas adequadas, por exemplo, temperatura adequada para o desenvolvimento (Müller et al., 2002).

A utilização da irrigação em gramíneas, conduzidas sob pastejo rotacionado, é uma técnica mais complexa do que se imaginava inicialmente. Segundo Lourenço et al. (2001) a evapotranspiração do capim Tanzânia posterior ao corte foi de 4 a 5 vezes inferior àquela apresentada ao final do ciclo. Tal resultado foi associado ao índice de área foliar (IAF) de 1,0 a 1,5 logo após o corte e de 4,0 a 5,5 ao final do ciclo. Os mesmos autores verificaram, também, que o coeficiente de cultura (Kc) variou de 0,30 a 0,40 na fase inicial de desenvolvimento e de 1,10 a 1,40 na fase final, para pastagem com ciclo de 36 dias. Tal variação do Kc influenciou diretamente no manejo de irrigação. Para exemplificar o fato, seja o sistema de pastagem rotacionada em pivô-central apresentada na Figura 1 com a rotação dos animais ocorrendo no sentido anti-horário. A parcela de número três é suposta como sendo a parcela em que os animais se encontram em pastejo; logo, seu IAF será mínimo ao final do pastejo nesta área, implicando em um coeficiente de cultivo mínimo (Kc_min), apresentando necessidade de irrigação também mínima; já a parcela quatro se encontra em pré-pastejo com o número máximo de dias de descanso (k_pmax) de acordo com o ciclo local e, daí, apresentará um desenvolvimento (IAF_max) maior, implicando em maior coeficiente de cultivo (Kc_max) e na necessidade de maior lâmina de irrigação. Verifica-se, seguindo-se o raciocínio acima, que as parcelas posteriores a quatro (p = 5, 6, ..., n, 1, 2) apresentarão um número menor de dias de descanso quando comparadas com a parcela 4, ou seja, menores IAF, Kc e taxas de evapotranspiração, quando comparada com a parcela quatro.

Todavia, na prática se verifica que é aplicada uma lâmina de irrigação uniforme em toda a área do pivô, sem que haja discriminação do grau de desenvolvimento de cada parcela, podendo acarretar perdas/excesso de água e energia, caso o solo não consiga atuar como reservatório de alta eficiência. Desta forma, o objetivo do presente trabalho foi desenvolver um modelo para estimativa da lâmina de irrigação necessária para suprir as necessidades de água do capim Tanzânia (Panicum maximun Jacq) de acordo com a determinação do Kc de uma parcela de controle e se considerando as condições climatológicas locais.

DESENVOLVIMENTO DO MODELO

Assume-se que o pivô-central tem capacidade finita de aplicar diferentes lâminas na área a ser irrigada podendo-se, então, dividir a área do pivô em p parcelas (Figura 2). Os animais pastejam no sentido anti-horário, retornando a determinada parcela depois de certo número de dias de descanso (k_pmax) conforme ciclo estabelecido para a região. Segundo Lourenço et al. (2001) a necessidade de irrigação de cada parcela varia de acordo com o estado de desenvolvimento da pastagem que, por outro lado, é dependente do número de dias que se encontra em descanso após o pastejo (k). No modelo se considera p = 1 como sendo a parcela de controle e que a lâmina a ser irrigada nesta parcela será obtida com auxílio de dados de tensão da água no solo, coletados através de baterias de tensiômetros, a diferentes profundidades para um monitoramento da umidade no perfil do solo.

O modelo que será apresentado tem, como suposição para fins práticos, uma relação direta entre lâmina de irrigação a ser aplicada e o seu respectivo coeficiente da cultura (Kc) e, desta forma, com o número de dias em que cada parcela está em descanso (k).

Lâmina de irrigação na parcela de controle monitorada por tensiometria

A lâmina bruta total a ser aplicada na parcela 1 (L₁) é calculada pela equação:

(1)

em que h é o número de camadas monitoradas e h é a eficiência do sistema de irrigação; q_CCé a umidade da capacidade de campo, cm³ cm^-3; q_obsé a umidade observada no solo, cm³ cm^-3; e Z a altura da camada a ser considerada, mm. q_obs será calculada por meio da seguinte equação (van Genuchten, 1980):

(2)

em que: q_sé a umidade de saturação do solo, cm³ cm^-3; q_ré a umidade residual do solo, cm³ cm^-3; a, n e m são os parâmetros; sendo

(3)

em que Y_j é a tensão média lida na superfície do solo (cm.c.a) observada das l-ésimas medidas dos tensiômetros (Y_i) na profundidade j.

Assumindo-se que a taxa de variação da lâmina de irrigação para uma parcela (L_p) com o seu respectivo Kc (Kc_p) é constante (e), tem-se, após integração, que

(4)

em que C₁é uma constante.

Na parcela de controle tanto a lâmina a ser aplicada como o coeficiente de cultura são conhecidos, logo L(Kc₁) = L₁. Segue, também, que quando Kc for zero, a lâmina a ser irrigada será zero; então, L(0) = 0; desta forma a Eq. 4 pode ser reescrita como:

(5)

Estimativa do coeficiente de cultura para cada parcela

O coeficiente de cultura para cada parcela (Kc_p) será modelado de duas maneiras devido a duas situações hipotéticas para o comportamento do Kc verificado por Lourenço et al. (2001):

a. Modelo linear para variação do Kc_p: Este modelo pressupõe que a taxa de variação dos Kc_ps é constante (f) com a variação k_p ou seja,

(6)

integrando esta equação,

(7)

Para resolver o problema de valor inicial, considera-se que determinada parcela (p) terá, no dia zero (início do pastejo) um consumo mínimo de água, a que se denomina coeficiente de cultura mínimo (Kc_min); desta forma, Kc(0) = Kc_min. Também, de acordo com as condições climatológicas do respectivo período analisado para o máximo número de dias de descanso (pré-pastejo) da parcela p, o Kc será máximo; isto é, Kc(k_pmax) = Kc_max. A Eq. 7 pode, então, ser reescrita da seguinte forma:

(8)

b. Modelo senoidal para variação do Kc_p: Para o modelo senoidal diz-se que a taxa de variação do Kc pode não ser constante para os números de dias de descanso de determinada parcela, mas ela será da forma senoidal, com uma taxa de variação inicial (dKc_p(0)/dy=d ) crescente até um máximo, que ocorreria na metade dos dias de descanso; na segunda metade dos dias de descanso ocorre uma queda na taxa de variação do Kc; assim sendo, y=f(k_p)=pk_p/ k_pmax, pode-se escrever:

(9)

integrando esta equação obtém-se

(10)

Resolvendo o problema de valor inicial com as condições conforme mostradas no modelo linear (Kc_p(0)=Kc_min e Kc_p(k_pmax)=Kc_max), tem-se que:

(11)

Cálculo do Kc_max ao longo do ano

A produção de matéria seca em pastagem é influenciada pela temperatura do ar (Müller et al., 2002; Lourenço et al., 2001). Villa Nova et al. (1983) propuseram utilizar a Unidade Fototérmica (UF) como maneira de mensurar a energia em determinada época. Para isto, sugerem a equação:

(12)

em que: são os graus-dia médios do período de n dias; Nf, o fotoperíodo (horas e décimos) no final do período de crescimento, e Ni o fotoperíodo, no início do período de crescimento.

Para modelar o Kc_max, com base na UF, utilizar-se-á um modelo logístico, com capacidade de tolerância para o coeficiente de cultura. Neste modelo, a capacidade de aumento de Kc diminui com o aumento da UF. Desta forma, pode-se escrever:

(13)

em que: T é a capacidade de tolerância para o Kc e t é uma constante. Nota-se, pela Eq. 13, que quando Kc é pequeno comparado com T, então Kc/T será próximo de zero; logo, dKc/dUF » tKc, ou seja, Kc tende a crescer de forma exponencial na direção de T; no entanto, se Kc®T (coeficiente de cultura tende à sua capacidade de tolerância), então Kc/T®1; logo, dKc/dUF®0, representando o equilíbrio com a capacidade do sistema, tendendo a permanecer estável. Resolvendo a Eq.13, tem-se:

(14)

sendo

(15)

A capacidade máxima de tolerância de Kc para o sistema será de 1,55 e o Kc₀ igual a 0,35 (Lourenço et al., 2001). Também, segundo Mello et al. (2001), Kc(6000) = 1,00. Substituindo-se esses valores na Eq. 15, tem-se

(16)

O modelo logístico pode ser utilizado para se determinar o Kc máximo diário em qualquer região para capim Tanzânia, desde que se tenha as séries de temperaturas mensais (utilizadas para calcular graus-dias) e a quantidade de dias em que a parcela permanece em descanso (k_pmax).

Dias de descanso de cada parcela

Tem-se, agora, o problema da determinação do número de dias de descanso para cada parcela (k_p) do pivô-central. A parcela de controle é a única em que o número de dias de descanso é acompanhado (k₁); a partir da segunda parcela, o número de dias de descanso será determinado por:

(17)

(18)

sendo p o número da parcela em análise.

A porção "k_pmax /p" é um fator de proporção para cada mudança de parcela, sendo que, quanto maior a quantidade de lâminas distintas que o pivô pode aplicar, considerando-se k_pmax constante, menor será esta porção e mais fracionada a distribuição de aplicação da água na área. Em equipamentos pivô-central com painel digital é possível se estabelecer de 8 a 12 setores molhados por operação automática.

Análise do modelo

Será apresentada a aplicação do modelo proposto para as condições de pastagem irrigada por pivô-central nas regiões de Piracicaba (21º 42" de latitude sul e 47º 38" de longitude oeste) e Pereira Barreto (20º 41" de latitude sul e 51º 7"de longitude oeste), no Estado de São Paulo; essas duas regiões foram escolhidas por apresentarem condições climatológicas bem distintas permitindo, assim, que se avalie o comportamento do modelo em situações diferentes. O clima da região de Piracicaba é Cwa, subtropical úmido, conforme classificação de Köppen, com verão chuvoso e inverno seco, sendo de 1230 mm a precipitação média anual; as temperaturas médias mensais variam de 25,4 ºC no verão a 17,9 ºC no inverno, com média anual de 21,6 ºC (LES-ESALQ/USP, 2002); já a região de Pereira Barreto apresenta precipitação média anual de 1128 mm, com temperaturas máximas e mínimas de 26,8 ºC (fevereiro) e 21,2 ºC (junho) com média anual de 24,7 ºC. A temperatura média mensal em Pereira Barreto é sempre superior à de Piracicaba, enquanto a precipitação em Piracicaba, na maioria dos meses é superior à de Pereira Barreto.

O manejo da irrigação por pivô-central em pastagem nas regiões de Piracicaba e Pereira Barreto, terá as seguintes características: i) pastagem com ciclo total de 36 dias; ii) irrigações realizadas em duas épocas, no dia primeiro de fevereiro (J=32, J é o dia juliano) e no dia primeiro de agosto (J=213) representando períodos com altas e baixas temperaturas, respectivamente; iii) irrigações baseadas em pivô-central com capacidade de se aplicar 12 lâminas distintas (p) na área, e iv) lâmina calculada para irrigação na parcela controle será conhecida, sendo L₁ = 10 mm, encontrando-se no 16º dia de descanso.

RESULTADOS E DISCUSSÃO

A Unidade Fitotérmica (UF) foi calculada para as duas regiões em estudo, com base em pastagem com o ciclo de 36 dias (UF₃₆), utilizando-se as séries históricas de temperatura média de cada região e localização geográfica. O comportamento geral da UF₃₆ de cada região é apresentado na Figura 3. Para a região de Piracicaba, a máxima UF₃₆é observada no final de fevereiro, enquanto a mínima ocorre em julho; já para Pereira Barreto, a UF₃₆ apresenta máxima em outubro e mínima no final de julho. Pode-se observar também que, em geral a UF₃₆ de Pereira Barreto é maior que a de Piracicaba durante quase todo o ano, resultado principalmente de suas maiores temperaturas médias mensais, visto que as latitudes são próximas, o que influencia pouco nos respectivos fotoperíodos das regiões.

Uma simplificação das curvas mostradas na Figura 3 por meio da utilização regressão múltipla de modelos polinomiais foi realizada, em que a UF₃₆ e o dia Juliano são as variáveis dependente e independente, respectivamente; desta forma, tem-se uma estimativa da UF₃₆ de maneira contínua no ano e não discreta. Verificou-se que o dia Juliano (J) utilizado nos modelos polinomiais do 4º grau explicou mais de 95% da variação da UF₃₆ para as duas regiões em estudo (Tabela 1).

Pode-se calcular, agora, o Kc_max de maneira contínua, para o pastejo com 36 dias de ciclo, com base no modelo logístico para Piracicaba (Eq. 19) e para Pereira Barreto (Eq. 20).

(19)

(20)

sendo T = - 0,000304(16957,00-104,65 J + 2,48 J² - 0,029 J³ - 1,23.10^-4J^-4

O modelo linear ou senoidal pode, então, ser utilizado para calcular o Kc de cada parcela em determinada época do ano. Considerando os Kc_max obtido pela Eq. 19 e 20, Kc_min igual a 0,35 e d=0,01 (Lourenço et al., 2001), obtiveram-se as Eqs. 21 e 22 para Piracicaba, e as Eqs. 23 e 24, para Pereira Barreto, referentes aos modelos linear e senoidal, respectivamente.

(21)

(22)

(23)

(24)

Nas Figuras 4 e 5 pode-se observar o comportamento geral do Kc em qualquer dia Juliano para o ciclo da pastagem nas duas regiões em estudo, a partir dos modelos linear e senoidal, respectivamente constatando-se, sobretudo, o contraste entre a variação do Kc para as duas regiões em ambos os modelos. Na região de Piracicaba há queda considerável do Kc entre aproximadamente os 150º e os 250º dias Juliano ao longo de todo o ciclo da cultura, com valores próximos ao Kc mínimo (0,35), enquanto para Pereira Barreto a queda do Kc é bem menor, cujo resultado é reflexo direto da baixa UF disponível para a cultura neste período, em Piracicaba, o mesmo não ocorrendo para a região de Pereira Barreto. De forma geral verifica-se, para a região de Piracicaba durante o inverno, uma necessidade menor de irrigação e se espera que a produtividade seja bastante afetada neste período, uma vez que a quantidade de energia disponível para o desenvolvimento da pastagem é limitada. Quanto à comparação entre os modelos (linear e senoidal), o modelo senoidal parece ser mais adequado ao período de altas temperaturas, apresentando taxas de variação distintas para o Kc, com taxa de crescimento menor no início e no final do desenvolvimento, enquanto o modelo linear responderia melhor em período de baixas temperaturas, conforme observações de campo (Lourenço et al, 2001); nota-se, também, que os valores de Kc no início e no final de um ciclo serão os mesmos para os dois modelos em cada região, cuja diferença estará apenas na taxa de variação do Kc.

Lâminas de irrigação (Estudo de caso)

Na simulação dos dias de descanso de cada parcela tem-se que os animais se encontram em pastejo na parcela 6, enquanto a parcela com maior quantidade de dias de descanso é a seguinte (parcela 7) (Tabela 2). Os Kc_maxs calculados por meio do modelo logístico para Piracicaba e Pereira Barreto, foram de 1,52 e 1,54 para fevereiro (J = 32), e de 0,42 e 1,30 para agosto (J = 213), respectivamente. Observa-se, mais uma vez, que o Kc_max de agosto é baixo para Piracicaba (0,42), resultado da baixa quantidade de energia disponível no período.

A Figura 6 apresenta a lâmina de irrigação pela Eq. 5, com Kc_p obtido pelos modelos linear e senoidal para a região de Piracicaba, em fevereiro. Como o cálculo da lâmina se dá de maneira discreta, não se observa variação significativa entre as estimativas a partir dos dois modelos, mas que o modelo linear apresenta o mesmo padrão de variação do modelo senoidal.

Na simulação para as épocas de fevereiro e agosto e se considerando a região de Piracicaba, o modelo mostrou-se sensível à capacidade de desenvolvimento de cada período, aqui relacionada com a UF₃₆ (Figura 7A). A temperatura média observada no mês em que a cultura esteve em desenvolvimento (julho) é de 17,9 ºC, enquanto em janeiro é de 24,4 ºC, implicando em que o Kc em agosto variou apenas de 0,35 a 0,42 (modelo senoidal) enquanto em fevereiro o Kc variou de 0,35 a 1,54 durante o ciclo total da pastagem; também, é bastante clara a sensibilidade do modelo quando foram consideradas as duas regiões para o mês de agosto (Figura 7B). A região de Pereira Barreto, apesar de também se encontrar em um mês de inverno, apresenta temperaturas médias mais elevadas que a região de Piracicaba para este período (21,2 contra 17,9 ºC), refletindo na maior variabilidade das lâminas; ainda para Piracicaba, como há pouco desenvolvimento para parcelas com maior ou menor número de dias de descanso, a tendência é se aplicar uma lâmina bastante próxima com aquela calculada na parcela de controle (10 mm).

A Figura 8 representa a diferença entre uma lâmina de 10 mm que seria aplicada em toda a área do pivô, desconsiderando-se a variação de desenvolvimento da pastagem entre as parcelas, com a estimativa de lâmina a irrigar realizada neste trabalho, para as regiões de Piracicaba e Pereira Barreto. Observa-se que nas primeiras parcelas haveria excesso de água, enquanto a partir da sétima parcela (inclusive) ocorreria um déficit de água em ambas as regiões prejudicando, assim, o desenvolvimento da pastagem, além de um desperdício de água e energia.

CONCLUSÕES

Foi apresentado, neste trabalho, um modelo para cálculo de lâminas de irrigação em pastagem irrigada sob pivô-central que, teoricamente, respondeu de maneira coerente, tanto no que diz respeito à aplicação de lâminas distintas de irrigação de acordo com o estádio de desenvolvimento da parcela como, também, à capacidade de desenvolvimento da pastagem de determinada região e período, baseada na unidade fototérmica.

LITERATURA CITADA

Protocolo 39 - 25/2/2003 - Aprovado em 6/2/2004

• LES-ESALQ/USP. Base de dados agroclimáticos Piracicaba, SP. http://ce.esalq.usp.br/bhbrasil/Saopaulo/ 15 de Nov. de 2002.
• Lourenço, L.F.; Coelho, R.D.; Soria, L.G.T.; Pinheiro, V.D.; Corsi, M. Coeficiente de cultura (Kc) do capim Tanzânia (Panicum maximum Jacq) irrigado por pivô central. In: Reunião Anual da Sociedade Brasileira de Zootecnia, 38, 2001, Piracicaba. Anais. Piracicaba: SBZ, 2001. p.1013-1017.
• Mello, A.C.L.; Lourenço, L.F.; Medeiros, H.R.; Pedreira, C.G.S.; Corsi, M. Dinâmica do índice de área foliar de capim Tanzânia em função do fotoperíodo e da temperatura do ar. In: Fisiologia de plantas no novo milênio: desafios e perspectivas, 2001, Ilhéus. Anais. Congresso Brasileiro de Fisiologia Vegetal, 8. Ilhéus, 2001. p. 130.
• Müller, M.S.; Fancelli, A.L.; Dourado Neto, D.; Garcia, A.G.; Ovejero, R.F.L. Produtividade do Panicum maximum cv. Mombaça irrigado, sob pastejo rotacionado. Scientia Agrícola, Piracicaba, v.59, p.427-433, 2002.
• van Genuchten, M. T. A closed form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society America Journal, Madison, v.44, p.892-898, 1980.
• Villa Nova, N.A.; Carreteiro, M.V.; Scardua, R. Um modelo de avaliação do crescimento de cana-de-açúcar (Saccharum spp.) em termos da ação combinada do fotoperíodo e da temperatura média do ar. In: Congresso Brasileiro de Agrometeorologia, 2, 1983, Campinas. Anais. Campinas: Sociedade Brasileira de Agrometeorologia, 1983. p.31-48.

Datas de Publicação

• Publicação nesta coleção
  26 Jan 2005
• Data do Fascículo
  Dez 2004

Histórico

• Recebido
  25 Fev 2003
• Aceito
  06 Fev 2004