RESUMO O objetivo do presente estudo foi avaliar as mudanças no metabolismo osmorregulador de genótipos de algodoeiro submetidos ao déficit hídrico e ao período de recuperação. Em casa de vegetação, seis genótipos de algodoeiro e dois manejos hídricos foram combinados em esquema fatorial 6 x 2, em delineamento inteiramente casualizado, avaliados após 14 e 22 dias em condições de estresse hídrico e recuperação, respectivamente. Foram avaliados o status hídrico, indicadores de crescimento e solutos compatíveis em folhas. Os genótipos BRS Seridó, BRS Aroeira e BRS 7MH suportam melhor a condição de estresse por apresentarem menores variações nos seus potenciais hídricos devido ao ajustamento osmótico e por manterem seu crescimento em condições de deficit hídrico. O BRS 286 é o genótipo mais sensível à condição de deficit hídrico, pois apresentou maiores variações no seu potencial hídrico que repercutiu na paralização do seu crescimento. A reposição hídrica durante o período de recuperação permitiu a retomada do crescimento da maioria dos genótipos, principalmente CNPA 5M e BRS 286. Após o período de recuperação hídrica, as variações nas concentrações dos osmorreguladores indicam plasticidade do algodoeiro em regular os níveis de osmorreguladores em condições favoráveis e desfavoráveis de água.

ABSTRACT The aim of this study was to evaluate the changes in osmoregulatory metabolism of cotton genotypes subjected to water deficit and recovery period. In a greenhouse, six cotton genotypes and two water managements were combined in a 6 x 2 factorial scheme, in a completely randomized design, evaluated after 14 and 22 days under water stress and recovery conditions, respectively. The water status, growth indicators and compatible solutes in the leaves were evaluated. The genotypes BRS Seridó, BRS Aroeira and BRS 7MH withstand better the stress condition because they had lower variations in their water potentials due to the osmotic adjustment and maintained their growth under water deficit conditions. The BRS 286 is the most sensitive genotype to the water deficit condition, since it showed greater variations in its water potential, which ultimately stopped its growth. Water replacement during the recovery period allowed growth resumption in most of the genotypes, mainly CNPA 5M and BRS 286. After the recovery period, variations in the concentrations of osmoregulators indicate the plasticity of cotton in regulating the concentrations of osmoregulators under favorable and unfavorable water conditions.

RESUMO O Índice de Área Foliar (IAF) é um parâmetro essencial na avaliação da resposta espectral da vegetação e estimativa das necessidades nutricionais e hídricas das culturas, porém é difícil obter sua estimativa de forma rápida e precisa em grandes áreas. O objetivo deste estudo foi a estimativa do IAF, da biomassa aérea e da produtividade do algodoeiro irrigado por meio de fotografia hemisférica digital. O experimento consistiu de uma combinação fatorial de quatro doses de nitrogênio (0, 90, 180 e 270 kg ha-1) e quatro de fósforo (0, 120 e 240 e 360 kg ha-1). As fotografias hemisféricas foram obtidas em condições similares de intensidade de brilho no céu aos 60 e 75 dias após a emergência (DAE) utilizando-se o Digital Plant Canopy Imager - CI-110® da CID Inc. O IAF também foi obtido por métodos destrutivos e não-destrutivos com o uso de um integrador de área foliar da LI-COR® modelo LI-3100C e por medições baseadas no comprimento da nervura central de todas as folhas. A análise de imagens hemisféricas obtidas por meio do Digital Plant Canopy Imager CI-110® constitui um método adequado para a estimativa do IAF e da biomassa do algodoeiro e para a estimativa da produtividade são necessários mais estudos para aperfeiçoar o método.

ABSTRACT The Leaf Area Index (LAI) is a key parameter to evaluate the vegetation spectral response, estimating plant nutrition and water requirements. However, in large fields is difficult to obtain accurate data to LAI determination. Therefore, the objective of this study was the estimation of LAI, biomass and yield of irrigated cotton through digital hemispherical photography. The treatments consisted of four nitrogen doses (0, 90, 180 and 270 kg ha-1) and four phosphorus doses (0, 120, 240 and 360 kg ha-1). Digital hemispherical photographs were collected under similar sky brightness conditions at 60 and 75 days after emergence (DAE), performed by the Digital Plant Canopy Imager - CI-110® of CID Inc. Biomass and LAI measurements were made on the same dates. LAI was also determined by destructive and non-destructive methods through a leaf area integrator (LI-COR® -LI-3100C model), and by measurements based on the midrib length of all leaves, respectively. The results indicate that the hemispherical images were appropriate to estimate the LAI and biomass production of irrigated cotton, while for the estimation of yield, more research is needed to improve the method.

RESUMOObjetivando avaliar o efeito da irrigação sobre a qualidade da fibra do algodoeiro, experimento foi conduzido no município de Barbalha, CE, no período de julho a dezembro de 2010. Os tratamentos consistiram da combinação fatorial de duas cultivares de algodoeiro herbáceo (BRS Aroeira e BRS Araripe) e de cinco lâminas de irrigação (260,93; 418,93; 514,21; 711,81 e 894,68 mm). Utilizou-se o delineamento em blocos casualizados dispostos em faixas, com quatro repetições. Foi utilizado um sistema de irrigação por aspersão em linha sendo a lâmina controle (711,81 mm) aplicada em função da evapotranspiração da cultura. As variáveis analisadas foram: percentagem de fibra, comprimento, uniformidade, índice de fibras curtas, resistência, alongamento à ruptura, índice micronaire, maturidade, grau de amarelamento, reflectância e índice de fiabilidade. As lâminas de irrigação influenciaram o comprimento, o índice de fibras curtas, a resistência, o micronaire, a maturidade e a reflectância da fibra. As cultivares influenciaram a porcentagem, o comprimento e a cor da fibra (grau de amarelamento). A cultivar BRS Araripe obteve a melhor qualidade de fibra com a lâmina de irrigação de 514,21 mm e a BRS Aroeira com a de 418,93 mm.

ABSTRACTAiming to evaluate the effect of irrigation depths on fiber quality of upland cotton, an experiment was conducted from July to December 2010 in Barbalha-CE, Brazil. The treatments consisted of a factorial combination of two upland cotton cultivars (BRS Aroeira and BRS Araripe) and five irrigation depths (260.93, 418.93, 514.21, 711.81 and 894.68 mm), arranged in a split-plot design with four replicates. A line-source sprinkler irrigation system was used and irrigation depth in the control treatment was calculated according to the crop evapotranspiration. The analysed fiber quality variables were: fiber percentage, length, uniformity, short-fiber index, resistance, elongation at rupture, micronaire index, maturity, degree of yellowing, reflectance degree and count strength product (CSP) index. The irrigation depths influenced fiber length, short-fiber index, strength, micronaire index, maturity and reflectance degree. The cultivars influenced fiber percentage, length and color (degree of yellowing). The best results of fiber quality were found with irrigation depths of 514.21 and 418.93 mm for the upland cotton cultivars BRS Araripe and BRS Aroeira, respectively.

<title>RESUMO</title><p>O objetivo deste trabalho foi avaliar o efeito de lâminas de irrigação sobre o rendimento e eficiência no uso da água de cultivares de algodoeiro no semiárido brasileiro. Foi conduzido um experimento de campo por dois anos consecutivos, na região de Apodi, RN, com irrigação por aspersão. O experimento consistiu de uma combinação fatorial em parcelas subdivididas, sendo quatro lâminas de irrigação (130; 100; 70 e 40%ETc) e quatro cultivares de algodoeiro (FiberMax 993, BRS 286, BRS 336 e BRS 335), em delineamento de blocos casualizados, com quatro repetições. Os dados foram avaliados através de teste de médias (Tukey) e análise de regressão. A produtividade com a menor lâmina aplicada (40% ETc) foi 48% menor se comparada a maior lâmina (130% ETc). A maior eficiência no uso da água foi encontrada com a lâmina de 70% ETc (0,69 kg m^-3) e a maior produtividade de algodão em caroço foi alcançada com a lâmina de irrigação igual a 130%ETc para todas as cultivares avaliadas. O fator de resposta do rendimento (Ky) foi igual a 0,632; 0,711; 0,784 e 0,858 para as cultivares FiberMarx 993, BRS 286, BRS 335 e BRS 336, respectivamente. As cultivares FiberMax 993 e BRS 286, no geral, apresentaram os melhores resultados, demonstrando serem mais indicadas para o cultivo irrigado na região semiárida.</p>

<title>ABSTRACT</title><p>The objective of this study was to evaluate the effect of irrigation depths on seed cotton yield and water-use efficiency of cotton cultivars in the Brazilian semi-arid region. Field experiment was conducted during two consecutive years in the Apodi region, RN, using sprinkler irrigation. The experiment consisted of factorial combination in split-plots, composed of four irrigation depths (130; 100; 70 and 40% ETc and four cotton cultivars - FiberMax 993, BRS 286, BRS 336 and BRS 335), in randomized block design with 4 replicates. Data were evaluated by mean test (Tukey) and regression analysis. Considering the irrigation depths of 40% ETc, cotton yield was 48% lower compared to the higher irrigation depth (130% ETc). The higher water-use efficiency (0.69 kg m^-3) was obtained with 70% ETc irrigation depth. The highest seed cotton yield was achieved with the higher water depth for all evaluated cultivars. Yield response factor (Ky) was equal to 0.632, 0.711, 0.784 and 0.858, for FiberMax 993, BRS 286, BRS 335 and BRS 336 cultivars, respectively. FiberMax 993 and BRS 286 cultivars presented the best performance, showing that they are more suitable for irrigated farming in the semi-arid region.</p>

Os modelos hidrológicos são ferramentas importantes que vêm sendo utilizadas no planejamento e na gestão dos recursos hídricos. Assim, objetivou-se com este trabalho calibrar e validar, numa escala de tempo diária, o modelo SWAT (Soil and Water Assesment Tool) para a bacia hidrográfica do Córrego do Galo, situada no Estado do Espírito Santo. Para a realização do estudo foram utilizados mapas georreferenciados de relevo, tipo e uso do solo, além de séries históricas diárias de clima e vazão da bacia. Na modelagem foram utilizadas séries temporais correspondentes aos períodos de 1º-1- 1995 a 31-12-2000 e 1º-1-2001 a 20-12-2013, para a calibração e validação, respectivamente. A avaliação do desempenho do modelo foi realizada usando os coeficientes de Nash-Sutcliffe (E NS) e percentual de viés (P BIAS). Também foram feitas avaliações do SWAT na simulação das seguintes variáveis hidrológicas: vazões máximas e mínimas diárias anuais; e vazões mínimas de referência Q90 e Q95, com base no erro médio absoluto. Os valores de E NS e P BIAS foram de 0,65 e 7.2% e 0,70 e 14.1%, na calibração e na validação, respectivamente, indicando desempenho satisfatório do modelo. O SWAT simulou adequadamente as vazões mínimas diárias anuais e as de referência para outorga Q90 e Q95 e não foi adequado na simulação das vazões máximas diárias anuais.

Hydrological models are important tools that have been used in water resource planning and management. Thus, the aim of this work was to calibrate and validate in a daily time scale, the SWAT model (Soil and Water Assessment Tool) to the watershed of the Galo creek , located in Espírito Santo State. To conduct the study we used georeferenced maps of relief, soil type and use, in addition to historical daily time series of basin climate and flow. In modeling were used time series corresponding to the periods Jan 1, 1995 to Dec 31, 2000 and Jan 1, 2001 to Dec 20, 2003 for calibration and validation, respectively. Model performance evaluation was done using the Nash-Sutcliffe coefficient (E NS) and the percentage of bias (P BIAS). SWAT evaluation was also done in the simulation of the following hydrological variables: maximum and minimum annual daily flowsand minimum reference flows, Q90 and Q95, based on mean absolute error. E NS and P BIAS were, respectively, 0.65 and 7.2% and 0.70 and 14.1%, for calibration and validation, indicating a satisfactory performance for the model. SWAT adequately simulated minimum annual daily flow and the reference flows, Q90 and Q95; it was not suitable in the simulation of maximum annual daily flows.

Propôs-se, neste trabalho, avaliar a variabilidade espacial dos atributos de fertilidade do solo em Latossolo Vermelho Amarelo em plantio direto para o algodoeiro em Cristalina, GO, visando à adoção da agricultura de precisão no manejo da área. A grade amostral foi constituída de 90 pontos (80 x 80 m) distribuídos numa área de 57,6 ha. Foram determinados os teores de P, K, Ca, Mg, S, Na, Cu, Mn, B, Zn, Fe, matéria orgânica (MO), pH, soma de bases (SB), saturação por bases (V), saturação por alumínio (m) e condutividade elétrica aparente do solo (CE) a 30 e 90 cm de profundidade. Efetuou-se a análise descritiva clássica e em seguida foram modelados semivariogramas para todos os atributos e mapas de krigagens. A dependência espacial das propriedades do solo analisadas indica que as variações espaciais devem ser consideradas no planejamento de coleta de amostras de solo e nas práticas de manejo do solo. Para a adoção da agricultura de precisão e aplicação de insumos em taxa variável, deve-se realizar um estudo criterioso da variabilidade da fertilidade do solo, gerando mapas por classe de fertilidade de cada atributo, o que permite identificar para quais atributos esta prática de manejo deve ser adotada.

The objective of this work was to analyze the spatial variability of fertility attributes in an Oxisol under tillage to cotton cultivation in Cristalina, GO, aiming the adoption of precision agriculture in the area management. The soil was collected using a sample grid, which consisted of 90 points (80 x 80 m), distributed into 57.6 ha. Were determined contents of P, K, Ca, Mg, S, Na, Cu, Mn, B, Zn, Fe, organic matter (OM), pH, sum of bases (SB), base saturation (V), aluminum saturation (m), as well as the apparent electrical conductivity of soil (EC) measured at 30 and 90 cm depth. Classic descriptive analysis was conducted and then, semivariograms were modeled for all attributes and krigging maps. The spatial dependence of analyzed soil properties indicates that the spatial variability should be considered during the planning of soil sampling and also to soil management practices. Concerning to precision agriculture adoption and the application of inputs in variable rates, it is required to do a careful study of the soil fertility variability, creating maps based on fertility classes for each attribute, allowing to identify for which attributes this management practice should be adopted.

A escassez de herbicidas registrados para a cultura da mamoneira ainda é fator limitante para sua expansão. Objetivou-se então, com este trabalho, avaliar a tolerância da mamoneira ao herbicida trifluralin em solos com diferentes texturas. O experimento foi conduzido em casa de vegetação, em delineamento inteiramente casualizado com cinco repetições, em esquema fatorial com 8 doses do herbicida trifluralin aplicado em pré-emergência (0; 187,5; 375; 750; 1.500; 3.000; 6.000 e 12.000 g ha-1 i.a.) e quatro tipos de substrato (areia lavada; franco-arenoso; franco argiloso-arenoso e franco-argiloso). Aos 28 dias após a semeadura foram determinadas as variáveis: altura de plantas, área foliar, massa seca da parte aérea, volume e massa seca do sistema radicular. Concluiu-se que a dose do herbicida trifluralin tolerada pela mamoneira é influenciada pelas características químicas e físicas do solo. Em solos de textura arenosa e franco-arenosa doses do herbicida trifluralin abaixo daquelas comumente utilizadas para controle de plantas daninhas, são suficientes para ocasionar redução no crescimento das plantas de mamoneira.

The scarcity of registered herbicides for cultivation of castor beans is still a limiting factor for its expansion. Thus, the aim of this study was to evaluate the tolerance of castor bean to trifluralin herbicide in soils of different textures. The experiment was conducted in a greenhouse in a completely randomized design with five replications in a factorial with 8 doses of the trifluralin herbicide (0; 187.5; 375; 750; 1.500; 3.000; 6.000 and 12.000 g ha-1 i.a.) applied at pre-emergence and four types of substrates (washed sand; sandy loam; sandy clay loam; and clay loam soil). At 28 days after sowing, the variables plant height, leaf area, shoot dry weight, volume and dry mass of roots were measured. It was concluded that the tolerance of castor bean to trifluralin herbicide is influenced by chemical and physical characteristics of soil. In sandy and sandy loam soils, doses of trifluralin below those commonly used for weed control are sufficient to cause the decrease in the growth of castor bean plants.

A mamoneira é uma espécie de crescimento inicial lento. A ausência de herbicidas registrados para esta cultura, ainda é fator limitante para sua expansão; assim, objetivou-se, com este trabalho, avaliar a tolerância da mamoneira ao herbicida pendimethalin em solos com diferentes capacidades de adsorção. O experimento foi conduzido em casa de vegetação, em delineamento inteiramente casualizado com cinco repetições e consistiu em uma combinação fatorial de oito doses do herbicida pendimethalin aplicado em pré-emergência (0; 187,5; 375; 750; 1.500; 3.000; 6.000 e 12.000 g ha-1 i.a.) e quatro tipos de substrato, sendo três solos de textura: franco-argilo-arenosa, franco-argilosa e franco-arenosa, além do substrato arenoso. Aos 22 dias após a semeadura determinaram-se as variáveis altura de plantas, área foliar, massa fresca e seca da parte aérea, volume e massa seca do sistema radicular. Concluiu-se que a dose do herbicida pendimethalin tolerado pela mamoneira é influenciada pela capacidade de adsorção do solo. O herbicida mostrou-se seletivo à cultura não ocasionando redução do crescimento nas doses utilizadas para controle de plantas daninhas.

The castor bean is a species of slow initial growth. The absence of registered herbicides for this crop is still a limiting factor for its expansion. The aim of this study was to evaluate the castor bean tolerance to the herbicide pendimethalin in soils with different adsorption capacities. The experiment was conducted under greenhouse conditions in a completely randomized design with five replications and consisted of a factorial combination of eight doses of the herbicide pendimethalin applied as pre-emergence (0; 187,5; 375; 750; 1500; 3000; 6000 and 12000 g ha-1 i.a.), and four types of substrate (clay loam, sandy clay loam, sandy loam and sandy soil). At 22 days after sowing the variables: plant height, leaf area, fresh and dry mass of shoots, root system volume and dry weight were determined. It was concluded that the dose of the herbicide pendimethalin tolerated by the castor bean is influenced by the adsorption capacity of the soil. The herbicide was found to be very selective to this crop, not causing significant reduction of plant growth at the doses indicated to be used for weed control.

RESUMO Este trabalho objetivou avaliar o desempenho do modelo de Green-Ampt-Mein-Larson (GAML) na simulação do processo de infiltração em solos com encrostamento superficial. Os ensaios de infiltração foram realizados num Cambissolo com a utilização de simulador de chuvas, em solo nu. O valor da condutividade hidráulica na zona de transmissão (Kt) foi considerado igual a Kt* que, por sua vez, é igual à taxa de infiltração estável (Tie) multiplicada pelo fator f, que representa o decaimento da Tie em função da energia cinética acumulada da chuva (Ec), ou seja, o valor de Kt variando ao longo do tempo. O modelo GAML com o valor de Kt constante ao longo do tempo não obteve bom desempenho superestimando, na maioria dos casos, os valores da taxa de infiltração (Ti) enquanto com o uso de Kt* o modelo apresentou bom desempenho, em que os melhores resultados foram obtidos com a combinação de Kt* com o potencial matricial (jf) calculado com uso da equação de Cecílio et al. (2007). O modelo de GAML com valor da condutividade hidráulica na zona de transmissão (Kt) variando no tempo apresentou bons resultados na simulação do processo de infiltração em solos sujeitos ao encrostamento superficial.

ABSTRACT The soil crust affects significantly the soil water infiltration rate. Thus, the infiltration simulation models must somehow consider the effect of crust layer to obtain good results.The objective of this work was to evaluate the performance of the Green-Ampt-Mein-Larson (GAML) model to simulate the soil water infiltration for crusting soils. The trials were carried out in a Podzol Tb distrophic Haplic Cambisol, using a rain simulator, on a bare soil. The GAML model parameters were determined , being proposed for the hydraulic conductivity of transmission zone (Kt) value the use of Kt*, which is equal to Tie x f, where Tie is stable infiltration rate and f is a decrease factor of the Tie as a function of cumulative kinetic energy of rainfall (Ec), i.e., Kt value varying over time. The GAML model with constant Kt value over time did not provide good performance, overestimating the values of infiltration rate (Ti) in most cases, where as the use of Kt* the GAML model showed good performance, being the best results obtained by combination of Kt* with matric potential (jf) values calculated by Cecílio et al. (2007) equation. The GAML model with the hydraulic conductivity in the transmission zone (Kt) value varying over time showed good results in simulation of the infiltration process in soils subjected to crusting.

RESUMO Neste trabalho a meta foi avaliar o modelo de Green-Ampt-Mein-Larson (GAML) e GAML modificado por Moore (1981) na simulação da infiltração de água em solos com encrostamento superficial. Os ensaios de infiltração foram realizados em Cambissolo sem cobertura, com quatro perfis de precipitação aos quais foram aplicadas três precipitações sucessivas para cada perfil, espaçadas 24 h. Para a inclusão da camada de encrostamento superficial no modelo de GAML, o valor da condutividade hidráulica na zona de transmissão (Kt) foi igual à taxa de infiltração estável (Tie) obtida em solo nu; no modelo de GAML modificado por Moore (1981), foram propostas uma camada de crosta com espessura de 5 mm e condutividade hidráulica da crosta saturada (Kc) determinada a partir dos dados experimentais e igual a 10 e 25% da condutividade hidráulica do solo saturado (K0). Nas segunda e terceira precipitações as combinações que não consideraram o encrostamento superficial não apresentaram bons resultados, superestimando a taxa de infiltração ao longo do tempo. Para o modelo de GAML os melhores resultados foram obtidos com uso de Kt = Tie determinada em solo nu, enquanto no modelo GAML se obtive, para solos estratificados,valor de Kc = 10 e 25% de K0.

ABSTRACT The aim of this study was to evaluate the Green-Ampt-Mein-Larson (GAML) model and GAML model modified by Moore (1981) for layered soils, in the simulation of water infiltration process in crusting soil. The infiltration tests were performed on a Cambisol without coverage and three successive precipitations were applied at intervals of 24 h. For the inclusion of the soil crust layer in the input parameters of the GAML model, the value of hydraulic conductivity in the transmission zone (Kt) was equal to stable infiltration rate (Tie) obtained in bare soil, and in the GAML model for layered soils, a surface layer (crust) were proposed with a thickness of 5 mm and soil crust saturated hydraulic conductivity (Kc) determined from experimental data, and equal to 10 and 25% of the saturated hydraulic conductivity (K0). In the second and third rainfall, the combinations without consideration of crust did not achieve good results and overestimated infiltration rate over time. For the model of GAML, the best results were obtained using Kt = Tie determined in bare soil, while in Model GAML for layered soils, were obtained with a value of Kc = 10 and 25% of K0.

O modelo de Green-Ampt-Mein-Larson (GAML) é um dos mais empregados na modelagem da infiltração; entretanto, alguns de seus parâmetros não condizem com a situação real do processo de infiltração da água no solo. Assim, o presente trabalho teve como objetivo avaliar o desempenho do modelo de GAML, testando diferentes combinações de metodologias de obtenção dos parâmetros de entrada. Os ensaios foram realizados em um Cambissolo Háplico Tb Distrófico Podzólico, com utilização de simulador de chuvas, em solo com e sem cobertura. Realizaram-se simulações baseadas na combinação de duas metodologias para a determinação da umidade na zona de transmissão (θt); duas para condutividade hidráulica na zona de transmissão (Kt) e três para o potencial matricial na frente de umedecimento (ψf). O modelo GAML com seus parâmetros originais não obteve bom desempenho, superestimando os valores de taxa de infiltração (Ti) e infiltração acumulada (I). As combinações em que se utilizaram valores de ψf calculados pela equação de Risse et al. (1995) subestimaram, durante todo o tempo, os valores de Ti e I, em ambas as condições de superfície. A combinação K0, θs e equação de Cecílio (2005) apresentou os melhores resultados.

Among the several equations that are proposed for simulation of infiltration, the Green-Ampt-Mein-Larson (GAML) model is one of the most used in hydrological studies as it is based on physical processes that occur in the soil during infiltration. However, some of its parameters do not match the real situation of the process of water infiltration in the soil. Thus, the objective of this work was to evaluate the performance of the GAML model, testing different methodologies for obtaining its parameters. The trials were carried out in a Podzol Tb Distrophic Haplic Cambisol, using a rain simulator, under field conditions, with and without mulch. Simulations were performed based on the combination of two methodologies for determination of soil water content in the transmission zone (θt), two for hydraulic conductivity in the transmission zone (θt) and three for the matric suction at the wetting front (ψf). The GAML model with its original parameters did not achieve good performance, overestimating the values of infiltration rate (Ti) and accumulated infiltration (I). Combinations using ψf values calculated by the equation of Risse et al. (1995) underestimated the whole time the values of Ti and I, under all surface conditions. The combination K0, θs and Cecílio (2005) equation showed the best results.