Estimativa de dimensões de bulbo molhado em irrigação por gotejamento superficial aplicando modelo de superfície de resposta

Estimate of wetted bulb dimensions in surface drip irrigation using response surface model

ARTIGOS CIENTÍFICOS

ENGENHARIA RURAL

Estimativa de dimensões de bulbo molhado em irrigação por gotejamento superficial aplicando modelo de superfície de resposta

Estimate of wetted bulb dimensions in surface drip irrigation using response surface model

Celsemy Eleutério Maia^I,1; Sérgio Luiz Aguilar Levien^I

^IDepartamento de Ciências Ambientais e Tecnológicas, Universidade Federal Rural do Semiárido (UFERSA), 59625-900, Mossoró, RN, Brasil. E-mail: [e-mail]

RESUMO

As informações de profundidade e largura do bulbo molhado na irrigação localizada são importantes tanto para o dimensionamento, quanto para o manejo da irrigação, auxiliando na determinação da quantidade de água aplicada para a planta. Com o objetivo de avaliar o modelo de superfície de resposta, para estimar as dimensões do bulbo molhado em função da vazão do emissor e do tempo de aplicação de água, foram instalados experimentos em seis solos. Os resultados obtidos mostraram que o modelo proposto pode ser utilizado para estimar as dimensões do bulbo molhado.

Palavras-chave: modelagem de água e solo, infiltração de água no solo, irrigação localizada.

ABSTRACT

The information of depth and width in wetted zone soil under drip irrigation plays great significance in design and management in the irrigation system for the amount of water required to the plant. For that, six experiments were installed in different soils with the objective of evaluating the dimensions of the wet bulb in the drip irrigation depending on the emitter discharged and of the water application time, using the power model. The obtained results showed that the proposed model can be used to estimate the dimensions of the wetted bulb for the all soils.

Key words: soil and water modeling, soil and water infiltration, localized irrigation.

INTRODUÇÃO

A irrigação localizada é um sistema de irrigação que tem como principais vantagens a aplicação de água perto das raízes das plantas e a economia de água pela maior eficiência de aplicação. Entre os sistemas de irrigação localizada, o gotejamento superficial é o mais comum e baseia-se na aplicação de água na superfície do solo junto à planta. Esse sistema permite molhar o solo, formando um volume de solo molhado, denominado de bulbo molhado ou bulbo úmido, onde se concentram as raízes das plantas. Molhando-se apenas uma porcentagem do solo, as perdas de água por evaporação e percolação são reduzidas, como também a lixiviação dos nutrientes (YOUNGS et al., 1999). Assim, o manejo da irrigação por gotejamento requer conhecimento do volume de solo molhado para evitar a sobreposição excessiva que pode diminuir a eficiência na aplicação de água e fertilizantes, principalmente em solos de textura arenosa que apresentam maior infiltração (KELLER & BLIESNER, 1990).

Nas últimas décadas, verificou-se um avanço no desenvolvimento de equipamentos, no dimensionamento e no manejo da irrigação localizada. Isso se deve principalmente ao melhor entendimento do movimento de água no solo quando se utiliza um ponto de emissão de água na superfície. Uma das condições para o melhor dimensionamento da irrigação localizada com ponto de emissão superficial é o conhecimento da distribuição da umidade no solo para diferentes vazões dos emissores e tempo de irrigação, sendo que a vazão e o tempo juntamente com o tipo de solo influenciam o movimento de água (LUBANA & NARDA, 1998). Resultados empíricos e simulados da distribuição de água no solo sob irrigação localizada enfatizam duas características práticas: a) quando o solo seco é irrigado, ocorre a formação do volume de solo molhado, em que a frente de molhamento separa o solo seco do solo molhado; b) a maior fração do volume de solo úmido tem um relativo teor uniforme de água. Esses dois fatores são utilizados na fertirrigação, em que a frente de molhamento define o volume de solo onde a planta emite raiz, e o teor de água e a concentração de nutriente na solução do solo podem ser mantidos próximos do ótimo, contribuindo para a maior eficiência na aplicação de fertilizante (BAR-YOSEF, 1999). No entanto, COELHO et al. (1999) destacam que os altos investimentos necessários para a implantação desse sistema de irrigação poderão não ser compensados se não forem utilizadas técnicas adequadas de manejo de irrigação que visem à racionalização do uso da água e ao aumento da produtividade.

Dentre os fatores que influenciam a formação do bulbo molhado, pode-se destacar o tipo de solo e a sua estratificação, a vazão do emissor e o tempo de irrigação (PIZARRO, 1996). Outros fatores que influenciam as dimensões do bulbo molhado são: estrutura do solo, sistema de irrigação, altura do emissor em relação ao solo, umidade inicial do solo e presença de concreções e cascalho (BELL et al., 1990; FOLCH & FÁBREGA, 1999). Nos solos mais argilosos, em razão da menor taxa de infiltração, verifica-se a formação de bulbo com maior dimensão horizontal, explicado pela maior influência da capilaridade sobre a gravidade. Em solos que apresentam camadas estratificadas, verificam-se camadas com diferentes porosidades, que afetam o fluxo e a retenção de água e, consequentemente, o bulbo molhado. Isso implica que, quando a frente de molhamento atinge camadas de diferentes texturas, esta se comporta como uma barreira, principalmente em solos que apresentam gradiente textural no perfil do solo. Para diferentes vazões de emissor e para mesmo tempo de irrigação, menores vazões tenderão a formar bulbos molhados mais profundos e com menor raio superficial; já para mesma vazão de emissores quanto maior o tempo de irrigação, maior o raio superficial.

O objetivo deste trabalho foi avaliar as dimensões do bulbo molhado na irrigação por gotejamento superficial em função da vazão do emissor e do tempo de aplicação de água, usando o modelo de superfície de resposta.

MATERIAL E MÉTODOS

Para uma avaliação do comportamento do bulbo molhado, foi feito um levantamento dos principais solos utilizados na fruticultura irrigada do Rio Grande do Norte, mais precisamente no Agropólo Assu-Mossoró. A região é caracterizada, segundo a classificação de Gaussen, com bioclima da região de 4ath tropical quente de seca a acentuada, com índice xerotérmico entre 200 e 150 e sete a oito meses secos. Segundo a classificação de Köppen, o clima da região é do tipo BSwh, isto é, seco, muito quente e com estação chuvosa no verão, atrasando-se para o outono e apresentando temperatura média anual de 27,4°C, precipitação pluviométrica anual bastante irregular, com média de 673,9mm e umidade relativa de 68,9% (CARMO FILHO et al., 1991).

Foram selecionados seis solos com diferentes características físicas e mineralógicas. De acordo com EMBRAPA (1999), os solos selecionados para o referido estudo foram classificados como Luvissolo Crômico, Argissolo Vermelho-Amarelo, Latossolo Vermelho, Cambissolo Háplico, Neossolo Quartizarênico e Neossolo Flúvico. Antes da montagem do sistema de irrigação, os solos foram preparados de forma a simular as reais condições de plantio, com aração e passagem da grade de disco de forma cruzada.

Para determinar com maior exatidão o tamanho do bulbo molhado (diâmetro da área molhada e profundidade), foi desenvolvido um sistema de irrigação portátil, em que foram dispostos oito emissores por linha, distribuídos aleatoriamente na parcela, ao longo de uma tubulação de polietileno de 16 mm de diâmetro, abastecida por dois depósitos de água (1000L e 500L). Para regular a vazão e a pressão, foram utilizados reservatórios reguladores, conectados aos primeiros, que mantinham um nível constante por meio de um sistema de bóia. As diferentes vazões foram obtidas pelos diferentes comprimentos dos microtubos inseridos na tubulação de polietileno. Estes, juntamente com a carga hidráulica constante, geraram diferentes vazões (1; 2; 4 e 8L h^-1) em suas saídas e eram fechados quando cada tempo de irrigação era concluído. O espaçamento entre os emissores foi estabelecido para que não ocorresse sobreposição na formação dos bulbos, ou seja, interferência dos emissores laterais.

Os experimentos foram instalados em delineamento em blocos inteiramente casualizados, no esquema de parcela subdividida, com três repetições, cujos tratamentos foram compostos pelos fatores tempo de aplicação de água e vazão do emissor, sendo a parcela formada pelo tempo e a subparcela formada pela vazão. As repetições constaram de duas linhas de irrigação, com quatro diferentes tempos (1; 2; 4 e 7h) e, dentro de cada tempo, as vazões de 1; 2; 4 e 8L h^-1 obtidas por microtubos de tamanhos diferenciados que proporcionam as vazões desejadas sob as mesmas condições de carga hidráulica.

Para cada vazão, após os tempos pré-determinados de 1; 2; 4 e 7h, foram abertas trincheiras no centro do bulbo molhado, abaixo do emissor, onde foram realizadas as medições das dimensões do bulbo com fita métrica, com precisão de 1,0mm. Para avaliar o efeito da vazão do emissor e do tempo de aplicação de água na formação do bulbo molhado, avaliaram-se as seguintes características: diâmetro superficial (D_s), diâmetro máximo (D_max), profundidade máxima (Z_max) e profundidade onde ocorre D_max (Z_Dmax). Além dessas medidas, foram medidos os diâmetros a cada 5 ou 10cm de profundidade para avaliação do formato do bulbo molhado. Para avaliar as dimensões do bulbo molhado em função de q e t, foram realizadas análises de regressão a fim de avaliar o comportamento do bulbo pelo modelo de superfície de resposta. Para tanto, observou-se, em outros trabalhos, que as dimensões do bulbo em função de t ou q isoladamente podem ser estimadas com o modelo potencial (THORBURN et al., 2003; LI et al., 2004); assim, a partir de y = α₁.tβ¹ e y = α₁.qβ¹ e da mesma medida do bulbo, tem-se y = α₁.tβ^1.α₂.qβ². Sendo α= α₁. α₂, tem-se y = α.tβ¹.qβ² ou y = a.t^b.q^c, em que t é o tempo de aplicação de água (h), q é a vazão do emissor (L h^-1) e a, b e c são parâmetros do modelo.

Para avaliar o uso de modelo de superfície de resposta para estimativa das dimensões do bulbo em função de q e t, foi utilizada a metodologia proposta por GAUCH et al. (2003), baseada no Quadrado Médio dos Desvios (MSD), que é decomposto em Quadrado da Bias (SB), Quadrado Médio para Declividade não Unitária (NU) e Quadrado Médio para Falta de Correlação (LC), sendo MSD = SB+NU+LC. O SB está relacionado com o intercepto da reta ajustada entre os valores estimados e observados, o NU está relacionado com a declividade da reta e o LC, com a correlação entre os dados observados e estimados pelos modelos e que quanto menor o valor do MSD, melhor será o modelo.

RESULTADOS E DISCUSSÃO

Pelo modelo de superfície de resposta ajustado aos dados (Tabela 1), observa-se que, com exceção de Z_max, Z_Dmax e D_max, para o Neossolo Flúvico, os coeficientes de determinação foram superiores a 90% para todas as dimensões dos solos avaliados. Levando-se em consideração que as dimensões do bulbo podem ser estimadas isoladamente em função de q ou t e utilizando-se o modelo potencial, uma das vantagens da superfície de resposta é que apenas 24 equações são suficientes para abranger todas as combinações de t e q do experimento, ao contrário do modelo potencial, em que esse número seria de 192 equações.

Os valores estimados das dimensões do bulbo molhado em função dos valores observados indicam que o modelo de superfície de resposta estimou bem essas dimensões e, além dos altos coeficientes de determinação, observa-se que os coeficientes angulares e os interceptos da reta ajustada foram estatisticamente iguais a 1 e 0, respectivamente (Figura 1).

Ao se fazer uma simulação para a cultura do melão com os resultados obtidos, para uma profundidade do bulbo molhado de 40cm e utilizando-se a equação Z_max = ƒ(t,q), tem-se o modelo t = ƒ(q), que combina os valores de t e q para obter a profundidade máxima de 40cm. Utilizando-se o critério de para melhor combinação de t e q, obtém-se o valor de q, com o qual, sendo substituído em t = ƒ(q), determina-se o tempo t para atingir Z_max igual a 40cm. Os valores de D_max, D_s e Z_Dmax foram obtidos pela substituição dos valores de t e q em D_max = ƒ(t,q), D_s = ƒ(t,q) e Z_max = ƒ(t,q), respectivamente. A partir do critério proposto por KELLER & BLIESNER (1990), multiplicando-se o valor de D_max por 0,8, tem-se o espaçamento entre os emissores.

Assim, sem considerar o Neossolo Flúvico, os valores de q para um Z_max de 40cm variaram de 3,019 a 3,446L h^-1, para o Argissolo Vermelho-Amarelo e o Cambissolo Háplico, respectivamente, e o tempo de irrigação variou de 2,864 a 3,975h, para o Argissolo Vermelho-Amarelo e o Luvissolo Crômico, respectivamente (Tabela 2). Os valores obtidos para o Neossolo Flúvico não foram comparados em razão da profundidade máxima observada no campo na maior vazão e do tempo de aplicação de água ter sido de aproximadamente 30cm. Desse modo, os dados desse solo apresentados na tabela 2 são extrapolações, ou seja, os resultados estão fora do espaço experimental, e a principal explicação para os valores elevados de t e q para se alcançar os 40cm é a baixa condutividade hidráulica desse solo.

Na avaliação do volume de água aplicado para a formação do bulbo, observa-se que, para atingir os 40cm de profundidade, o volume de água foi diferente para os solos estudados, sendo maior no Luvissolo Crômico e menor no Argissolo Vermelho-Amarelo, sem considerar o Neossolo Flúvico. Apesar do volume de água aplicado para atingir um Z_max de 40cm ter variado entre os solos, quando se verifica a lâmina de água aplicada pelo sistema (mm h^-1), os valores se aproximam, com média de 3,788mm h^-1. Isso se deve a diferença no espaçamento entre os emissores, considerado como sendo igual a 0,8 D_max, e este difere entre os solos.

Na avaliação do modelo da superfície de resposta pelo MSD, apenas para o Neossolo Flúvico os valores foram elevados, principalmente para D_max, D_s e Z_Dmax. Para os demais solos, observa-se que houve bom ajuste do modelo aos dados observados e que o principal componente afetado é a falta de correlação, porém com valores todos superiores a 93%, com exceção do Neossolo Flúvico (Tabela 3).

CONCLUSÃO

O modelo superfície de resposta pode ser utilizado para estimar os valores de D_max, Z_max, D_s e Z_Dmax em função da vazão do emissor e tempo de aplicação de água na irrigação por gotejamento superficial.

Recebido para publicação 16.11.09

Aprovado em 05.05.10

Devolvido pelo autor 18.06.10

CR- 2807

• BAR-YOSEF, B. Advances in fertigation. Advances in Agronomy, v.65, p.1-77, 1999.
• BELL, J.P. et al. Soil water status: a concept for characterizing soil water conditions beneath a drip irrigated row crop. Agricultural Water Management, v.17, p.171-187, 1990. Disponível em: <http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6T3X-4888H6h-1V&_user=684796&_coverDate=01%2F3 1%2F1990&_rdoc=16&_fmt=high&_orig=browse&_srch=doc-info(%23toc%234958%231990%23999829998%23411130%23FL P%23display%23Volume)&_cdi=4958&_sort=d&_docanch or=&_ct=35&_acct=C000036759&_version=1&_urlVe rsion=0&_userid=684796&md5=b89a4eed25d8637586c15fc600c25db4>. Acesso em: 01 abr. 2010. doi: 10.1016/0378-3774(90)90064-6.
• CARMO FILHO, F. et al. Dados meteorológicos de Mossoró (Jan. de 1988 à Dez. de 1990). Mossoró: ESAM/FGD, 1991. 121p. (Coleção Mossoroense).
• COELHO, E.F. et al. Avaliação de parâmetros hidráulicos para modelos de distribuição de água no solo sob gotejamento. Pesquisa Agropecuária Brasileira, v.34, n.4, p.651-657, 1999.
• EMBRAPA. Centro Nacional de Pesquisas de Solos. Sistema brasileiro de classificação de solos Rio de Janeiro, 1999. 412p.
• FOLCH, J.R.G.; FÁBREGA, J.A.G. El volumen húmedo del suelo en el riego localizado: importancia y evaluación. In: CARPENA, R.M. et al. (Eds). Estudios de la zona no saturada del suelo Tenerife: ICIA, 1999. p.11-17.
• KELLER, J.; BLIESNER, R.D. Sprinkle and trickle irrigation New York: AVI Book, 1990. 640p.
• LI, J. et al. Wetting patterns and nitrogen distributions as affected by fertigation strategies from a surface point source. Agricultural Water Management, v.67, p.89-104, 2004. Disponivel em: <http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6T3X-4C53PYC-2&_user=684796&_coverDate=06%2F15%2F 2004&_rdoc=3&_fmt=high&_orig=browse&_srch=doc-info(%23toc%234958%232004%23999329997%2350 1774%23FLA%23display%23Volume)&_cdi=4958&_sort=d&_docanchor=&_ct=7&_acct=C000036759&_version=1&_urlV ersion=0&_userid=684796&md5=16cdd650450ea83124060ca769e95d44>. Acesso em: 01 abr. 2010. doi: 10.1016/j.agwat.2004.02.002.
• LUBANA, P.P.S.; NARDA, N.K. Soil water dynamics model for trickle irrigated tomatoes. Agricultural Water Management, v.37, p.145-161, 1998. Disponivel em: <http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6T3X-3VGM9SB-5&_user=684796&_coverDate=07%2F01%2F19 98&_rdoc=5&_fmt=high&_orig=browse&_srch=doc-info(%23toc%234958%231998%23999629997%2346069%2 3FLA%23display%23Volume)&_cdi=4958&_sort=d&_docanchor=&_ct =8&_acct=C000036759&_version=1&_urlVersion=0&_userid=68 4796&md5=e35f0b7679e8e90c59ff08212a4dc19d>. Acesso em: 01 abr. 2010. doi: 10.1016/S0378-3774(98)00035-3.
• PIZARRO, F. Riegos localizados de alta frecuencia (RLAF): goteo, microaspersión, exudación. Madrid: Mundi-Prensa, 1996. 511p.
• SCHWARTZMAN, M.; ZUR, B. Emitter spacing and geometry of wetted soil volume. Journal of Irrigation and Drainage Engineering, v.112, n.3, p.242-253, 1986.
• THORBURN, P.J. et al. Soil-dependent wetting from trickle emitters: implications for system design and management. Irrigation Science, v.22, p.121-127, 2003. Disponivel em <http://www.springerlink.com.w10013.dotlib.com.br/content/14g42wquu4w0896r/fulltext.html>. Acesso em: 01 abr. 2010. doi: 10.1007/s00271-003-0077-3.
• VAN GENUCHTEN, M.T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of American Journal, v.44, p.892-898, 1980.
• YOUNGS, E.G. et al. Surface ponding of coarse-textured soil under irrigation with a line of surface emitters. Journal of Agricultural Engineering Research, v.73, p.95-100, 1999.
• ZUR, B. Wetted soil volume as a design objective in trickle irrigation. Irrigation Science, v.16, p.101-105, 1996.

1

Autor para correspondência.

Datas de Publicação

• Publicação nesta coleção
  09 Jul 2010
• Data do Fascículo
  Jun 2010

Histórico

• Recebido
  16 Nov 2009
• Aceito
  05 Maio 2010