Spatial variability of yield and other parameters associated with pear trees

Variabilidade espacial da produtividade e outros parâmetros associados com pés de pera

TOPOGRAPHY, PHOTOGRAMMETRY AND REMOTE SENSING

Spatial variability of yield and other parameters associated with pear trees

Variabilidade espacial da produtividade e outros parâmetros associados com pés de pera

Marcia R. S. Konopatzki^I; Eduardo G. Souza^II; Lúcia H. P. Nóbrega^III; Miguel A. Uribe-Opazo^IV; Grazieli Suszek^V

^IEngª Agrícola, Pós-Graduanda da UNIOESTE/CASCAVEL/CCET/PGEAGRI/SBA, Grupo de Pesquisa GROSAP, Cascavel - PR, [e-mail]. Bolsista da CAPES

^IIEngº Mecânico, Prof. Associado da UNIOESTE/CASCAVEL/CCET/PGEAGRI/SBA, Grupos de Pesquisa GROSAP e GGEA, Cascavel - PR, [e-mail]. Pesquisador de Produtividade do CNPq

^IIIEngª Agronoma, Profa. Associada da UNIOESTE/CASCAVEL/CCET/PGEAGRI/SBA, Grupos de Pesquisa GROSAP e GGEA, Cascavel - PR, [e-mail]. Pesquisador de Produtividade do CNPq

^IVEstatístico, Prof. Associado da UNIOESTE/CASCAVEL/CCET/PGEAGRI/SBA, Grupos de Pesquisa GROSAP e GGEA, PR, [e-mail]. Pesquisador de Produtividade do CNPq

^VEngª Agrícola, Profa. EBTT Instituto Federal do Mato Grosso do Sul (IFMS) - Nova Andradina - MS, [e-mail]

ABSTRACT

In the last few years, precision agriculture has become commonly used with many crops, particularly cereals, and there is also interest in precision horticulture. Pear is a seasonal fruit and well appreciated by Brazilian people, although it is mostly imported. Brazilian farmers are nowadays trying to increase pear production. Thus, this research aimed at mapping the yield of pear trees in order to study the spatial variability of yield as well as its comparison with spatial variability of soil and plant attributes. The experimental field had 146 pear trees, variety 'Pêra d'água', distributed on a 1.24 ha. Four harvests were performed according to the fruit ripening and from each tree; only the ripe fruits were harvested. In each harvest, all the fruits were weighed and the total yield was obtained based on the sum of each harvest. The soil attributes analyzed were P, K, Ca, Mg, pH in CaCl₂, C, Cu, Zn, Fe, Mn and base saturation, and the plant attributes were fruit length, diameter and yield. Yield had low correlation with soil and plant attributes. An index of spatial variability was suggested in this study and helped in classifying levels of spatial dependence of the various soil and plant attributes: very low (fruit length); low (P, fruit diameter), medium (Mg, pH, Cu, Zn, Fe), high (Ca, K, base saturation and yield), and very high (Mn and C).

Keywords: precision agriculture, thematic maps, horticulture.

RESUMO

Nos últimos anos, a agricultura de precisão tornou-se comumente utilizada em muitas culturas, particularmente cereais, mas há também interesse em fruticultura de precisão. Pera é uma fruta sazonal e muito apreciada no Brasil, apesar de ser principalmente importada. Agricultores brasileiros estão atualmente tentando aumentar a produção de peras. Assim, esta pesquisa visou mapear a produtividade de pés de peras com o objetivo de estudar sua variabilidade assim como relacionar a variabilidade espacial dos atributos do solo e da planta. A área experimental tem 146 pés de pera, variedade Pera-d'água, distribuída em uma área de 1,24 ha. Quatro colheitas foram realizadas de acordo com a maturação da fruta e somente as frutas maduras de cada pé foram colhidas. Em cada colheita, todas as frutas foram pesadas, e a produtividade foi obtida baseada na soma de cada colheita. Os atributos do solo analisados foram P, K, Ca, Mg, pH em CaCl₂, C, Cu, Zn, Fe, Mn e saturação de bases, e os atributos da planta foram comprimento e diâmetro do fruto e produtividade. A produtividade teve baixa correlação com os atributos do solo e da planta. Um índice de variabilidade espacial foi proposto neste estudo e auxiliou na classificação dos níveis de dependência espacial de diferentes atributos do solo e da planta: muito baixo (comprimento de fruta); baixo (P, diâmetro do fruto); médio (Mg, pH, Cu, Zn e Fe), alto (Ca, K, saturação de base e produtividade) e muito alto (Mn e C).

Palavras-chave: agricultura de precisão, mapas temáticos, fruticultura.

INTRODUCTION

Precision agriculture (PA) is the management of agricultural production according to the characteristics and needs of small production units. Sometimes these units may be based on a single plant or animal. PA aims to improve the use of inputs, thus increasing profit for farmers and preserving environmental quality. This practice has several steps, including data collection, management planning and site-specific application of agricultural inputs.

During data collection, the objective is to identify the variability of production factors, as well as of cropping in the field. Usually a yield map is created to visualize spatial variability. However, this map is only one of the steps involved in the whole PA process and represents the combined effect of several sources of spatial and temporal variability. Part of this variability can be attributed to constant factors or to factors which vary slowly. Others are transitory, changing in their importance and/or spatial and temporal distributions from one season to another.

Pear is one of the favorite seasonal fruits in Brazil, but most are imported. In 2005, Brazil produced approximately 16,000 metric tons of pear, in an area of 1,760 ha, while the world production was 20 million tons. In 2007, Brazilian pear imports reached 103,000 tons (FIORAVANÇO, 2007). Therefore, pear production is an attractive agricultural alternative to farmers of the Brazilian Southern Region, which is edaphically and climatically prepared for such a crop. Only apple production is more economically advantageous than pear production in this region.

Although the majority of Brazilian and worldwide researches in PA are applied to grain crops, some researches in horticulture have already been published. In Brazil, FARIAS et al. (2003), MOLIN & MASCARIN (2007), LEÃO et al. (2010), and SIQUEIRA et al. (2010) worked with citrus; SILVA et al. (2007), with coffee; CALDERON et al. (2008), with mango; and MIRANDA et al. (2005), with melon. Internationally, AGGELOPOULOU et al. (2010a, b, c) worked with apples; TISSEYRE & Mcbratney (2008) with grapevine; and PERRY et al. (2009) with pear.

In addition to the lack of knowledge of the root system and of specific fertilizer application techniques, defined criteria for the better recommendation of inputs in perennial plants do not exist yet. The fruit trees explore large soil volumes and they differ from annual plants because their roots, stems, and branches can store nutrients from one year to the next. Generally, attempts to correlate values from soil analysis with the nutrient needs of perennial plants have not been successful. This probably happens because it is almost impossible to collect a sample to represent the soil mass explored by the roots and also because perennial plants have a different capacity for extracting nutrients from the soil than annual plants (FACHINELLO et al., 2008).

Therefore, this experiment mapped the yield of pear trees in order to analyze the spatial variability of yield, as well as to compare it with spatial variability of soil and plant attributes.

MATERIAL AND METHODS

The experimental area (Figure 1) is located in Nova Laranjeiras, State of Paraná, South of Brazil (25º23'22"S; 52º34'15"W; altitude 750 m). The soil is classified as Distroferric Red Latossol (EMBRAPA, 1999) with 70% of clay.

For georeferencing (Figure 2), a Trimble Geo Explorer 3 Geodesic Differential Global Positioning System (GDGPS) with post-processing was used. The experimental area has 146 pear trees, distributed in a 1.24 ha and originally farmed under no-tillage with crops of corn, soybean, oats and sorghum. The orchard was established in 2000, using an 8 m distance between rows and 10 m between plants within the row. The pear variety "Kieffer" was used as rootstock upon which the pear variety "Pêra d'água" was grafted. In this study, yield data from the 2007 harvest was analyzed.

At the orchard establishment, the soil of each planting hole was fertilized with phosphorus, lime and organic manure. Every subsequent year, additional fertilization with N and P was performed. The first harvest was accomplished on February 2006. Once Nova Laranjeiras, State of Paraná, is not considered climatically ideal to pear production, dormancy had to be chemically overcome to enable yield increases.

Four harvests were performed based on fruit maturation, and in each harvest, only the ripe fruits were picked, based on the farmer experience. The weight, diameter and length of fruits from 36 pear trees were computed and the total yield was then obtained by the sum of the individual yield of each plant (total weight). From these 36 trees, 18 were selected from a 30 m regular grid and the others 18 are the closest to these ones.

Soil samples were collected on a grid with a 30 m constant spacing and evaluated for P, K, Cu, Zn, Fe, Mn, Ca, Mg, pH (Calcium Chloride) and C. The base saturation (V%) were calculated. Each sample was composed of four sub-samples, randomly taken within one-meter radius. The sample depth was 20 cm. Before data interpolation, the Pearson linear correlation coefficient (R) and the Spearman rank-order correlation coefficient (R_S) were used to evaluate the extent of linear (R) and non linear (R_S) correlation among variables. In geostatistical analysis, semivariograms were built to verify the spatial dependence among samples. In order to estimate the structures of the experimental semivariance function, the classic Matheron estimator (variables with normal distribution) and Cressie and Hawkins estimator (non-normal distribution variables) were used. The extent of the semivariogram spatial dependence was evaluated using the spatial variability index (SVI, eq.(1)), proposed by the authors.

where, partial sill and nugget effect are parameters of the semivariogram.

The SVI classification proposed by the authors and adopted was: very low for SVI < 20%; low for 20 < SVI < 40%; medium for 40 < SVI < 60%; high for 60 < SVI < 80%; and very high for SVI > 80%. This indicator is similar to the one proposed by CAMBARDELA et al. (1994), with the advantages of having five interpretation levels instead of three and being proportional to the spatial variability (the higher SVI the higher the spatial variability).

The Anderson-Darling and Kolmogorov-Smirnov tests were used to verify data normality at 5% probability. The data were considered normally distributed when fitting at least, one of the tests. Outliers were verified through box-plot graphs. The coefficient of variation (C.V.) was considered low when C.V. < 10%; medium when 10% < CV < 20%; high when 20% < C.V. < 30%; and very high when C.V. > 30% (PIMENTEL-GOMES & GARCIA, 2002).

To help compare the distribution values of each nutrient, it was proposed the nutritional sufficiency index - NSI (%):

where,

OBS - observed nutrient content, and

LIM_INF - lower limit of the 'high' range for the nutrient. This value is frequently defined by soil scientists.

A multiple regression model was developed for yield as a function of chemical attributes of soil and plant characteristics. The best subset variable selection method was used. A 0.05 significance value for the F distribution was used to control entering and removal the effects. The adjusted multiple coefficient of determination (R² adjusted) was used as selection criterion for the best model for the specified effects. After an exploratory analysis, the decision was to fit the data to the following model:

where,

Pr - response variable: yield, kg ha^-1;

P, K, pH, Cu, Fe, Mn and Sat - soil attributes;

Leng - fruit length;

Diam - fruit diameter;

a, b ... r - model parameters to be estimated by least squares method, and

ε - random error.

RESULTS AND DISCUSSION

The computed attributes (table 1) showed a range of variability, including low (pH, length and diameter of fruit), medium (Mg and C), high (Ca and Cu) and very high (P, K, Zn, Fe, Mn and yield) (PIMENTEL-GOMES & GARCIA, 2002). Several attributes were normally distributed (P, K, Ca, Mg, C, Mn and length of fruit), but the others were not (pH, Cu, Zn, Fe, diameter of fruit and yield).

Soil attributes were in the range: very low (K), low (P, K, and base saturation), medium (P, K, Ca, Fe, and base saturation), and high (all attributes), according to the interpretation levels (table 2) suggested by COMISSÃO DE FERTILIDADE DO SOLO (1994). All attributes were at levels of nutrient availability sufficient for plant development, except that 28% of the area was lacking in P and 83% in K.

On average, the spatial variability characterized by the SVI increased with the range (m) (Figure 3). The R² was not good but the linear coefficient was significative at 0.05 significance value for the F test.

The C.V. (Figure 4) also increased with the SVI, and but the R² was low and the linear coefficient was not significative at 0.05 significance value for the F test. This fact should be investigated more carefully during the next years in order to have definitive conclusions.

The proposed nutritional sufficiency index (NSI, %) of soil attributes was displayed in box-plots (Figure 5). In this index, the 100% figure represents the lower value of the "high level" range of each attribute, which was herein considered the ideal level. The attributes were ordered from left to right, according to their increasing nutritional sufficiency. Special attention should be given to Zn, P, V% and Fe, since part of the samples concerning each attribute showed lower responses when compared to the optimum level. Attention should also be given to excess macro and micronutrients due to a possible occurrence of toxic action and interference with other nutrients. A classic example is the negative interaction between P and Zn (PERES, 2002) and the toxic action of some excess micronutrients (VALLADARES et al., 2007).

Pear yield had low linear correlation with P and K (Table 3), i.e., less than 12% (R²=0.12; R=0.35) of the yield variability was explained by these attributes. Furthermore, this correlation was lower for the remaining soil attributes. This low correlation probably occurred because all soil attributes were at least at a medium level of fertility, with the exception of P, which was low in 17% of the area. Fruit diameter explained 27% of yield variability (R²=0.27; R=0.52). Within the set, the most closely correlated attributes were: Ca and Mg (R = 0.86); Ca and pH (R = 0.86); Mg and pH (R = 0.89); base saturation and Ca (R = 0.88); base saturation and pH (R = 0.89).

Table 4

The models for the experimental semivariograms (Table 5, Figure 6) were: spherical (K, pH, C, Zn, length and fruit diameter); exponential (P); and Gaussian (Ca, Mg, Cu, Fe, Mn, base saturation and yield).

The data were all spatially dependent, ranging from very low (fruit length) to low (P, fruit diameter), medium (Mg, pH, Cu, Zn, Fe), high (Ca, K, saturation of bases and yield), and very high (Mn and C), according to the spatial variability index (SVI, %) (Figure 7). According to CAMBARDELLA et al. (1994), the attributes that present a strong spatial dependence are more influenced by intrinsic soil properties, such as texture and mineralogy. On the other hand, those that present a weak dependence are more influenced by external factors, such as fertilizer applications and soil tillage.

The thematic maps (Figure 8) built by ordinary Kriging show that, for each variable, minimum and maximum values were found in different locations, in accordance with the low correlation registered among the variables. The lowest yield was found at the east side of the area, indicating the need of further investigation for determining limiting factors such as soil depth and compaction, management history, drainage problems, and occurrence of pests, diseases and weeds.

Models were developed to describe the influence of soil and plant attributes on pear trees yield (Prod, kg ha^-1) (table 6). Such models allowed an explanation of 54% (R²=0.54) of yield variability, when only soil attributes were considered, and 90% (R²=0.90), when fruit length and diameter were included. The model fit is presented in figure 9.

CONCLUSION

The yield of pear trees had low linear correlation with the measured attributes (P, K, Ca, Mg, Al_3+, pH, CaCl₂, C, Cu, Zn, Fe, Mn, base saturation, length and diameter of fruits).

The variability of the measured attributes was considered very low (K), low (P, K, and base saturation), medium (P, K, Ca, Fe, and base saturation), and high (all attributes), according to the coefficient of variation.

Spatial dependence levels were found to be very low (fruit length), low (P, diameter of fruit), medium (Mg, pH, Cu, Zn, Fe), high (Ca, K, saturation of bases and yield), and very high (Mn and C), according to the spatial variability index that was proposed in this work.

ACKNOWLEDGEMENTS

The authors thank the support provided by the State University of Western Paraná (Unioeste), the Araucária Foundation (Fundação Araucária), the General office of State of Science, Technology and Higher Education - SETI/PR, the Coordination for the Improvement of Higher Level Personnel (CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior), the National Council for Scientific and Technological Development (CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico), and the Government of the State of Paraná.

Recebido pelo Conselho Editorial em: 16-5-2011

Aprovado pelo Conselho Editorial em: 9-11-2011

ERRATUM

In the paper "SPATIAL VARIABILITY OF YIELD AND OTHER PARAMETERS ASSOCIATED WITH PEAR TREES", with DOI number: 10.1590/S0100-69162012000200018, published in the journal Agricultural Engineering 32 (2):381-392, on the page 381:

Where it reads:

MARCIA R. S. KONOPATZKI¹, EDUARDO G. SOUZA², LÚCIA H. P. NÓBREGA³, MIGUEL A. URIBE-OPAZO⁴, GRAZIELE SUSZEK⁵

It should read:

MARCIA R. S. KONOPATZKI¹, EDUARDO G. SOUZA², LÚCIA H. P. NÓBREGA³, MIGUEL A. URIBE-OPAZO⁴, GRAZIELI SUSZEK⁵

• AGGELOPOULOU, A.D., BOCHTIS, D.; FOUNTAS, S.; SWAIN, K.C.; GEMTOS, T.A.; NANOS, G.D. Yield prediction in apple orchards based on image processing. Precision Agriculture Dordrecht, 2010a. Disponível em: <http://www.springerlink.com/content/n67l072m3474vvq1/>. Acesso em: 20 jan. 2011.
• AGGELOPOULOU, K.D.; PATERAS, D.; FOUNTAS, S.; GEMTOS, T.A.; NANOS, G.D. Soil spatial variability and site-specific fertilization maps in an apple orchard. Precision Agriculture Dordrecht, 2010b. Disponível em: <http://www.springerlink.com/content/23786817wj83u041/>. Acesso em: 20 jan. 2011.
• AGGELOPOULOU, K.D.; WULFSOHN, D.; FOUNTAS, S.; GEMTOS, T.A.; NANOS, G.D.; BLACKMORE, S. Spatial variation in yield and quality in a small apple orchard. Precision Agriculture, Dordrecht, v.11, n.5, p.538-556, 2010c.
• CALDERÓN, C.A.V.; MARTÍNEZ, L.J.M.; HENAO, R.G.Variabilidad espacial del suelo y su relacion con el rendimiento de mango (Mangifera indica L.). Revista Brasileira de Fruticultura, Jaboticabal, v.30, n.4, p.1146-1151, dez. 2008.
• CAMBARDELLA, C.A.; MOORMAN, T.B.; NOVAK, J.M.; PARKIN, T.B.; KARLEN, D.L.; TURCO, R.F.; KONOPKA, A.E. Field-scale variability of soil properties in central Iowa soils. Soil Science Society of America Journal, Madison, v.58, n.5, p.1501-1511, 1994.
• COMISSÃO DE FERTILIDADE DO SOLO - RS/SC. Recomendações de adubação e calagem para os estados do Rio Grande do Sul e Santa Catarina 3.ed. Passo Fundo: SBCS - Núcleo Regional Sul/ EMBRAPA/CNPT, 1994. 224 p.
• EMBRAPA - EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA. Serviço Nacional de Levantamento e Conservação de Solos. Sistema Brasileiro de classificação de solos Brasília: Embrapa Produção de Informação, 1999. 412 p.
• FACHINELLO, J.C.; NACHTIGAL, J. C.; KERSTEN, E. Fruticultura: fundamentos e práticas. Pelotas: Embrapa Clima Temperado, 2008. 311 p.
• FARIAS, P.R.S.; NOCITI, L.A.S.; BARBOSA, J.C.; PERECIN, D. Agricultura de precisão: mapeamento de produtividade em pomares cítricos usando geoestatística. Revista Brasileira de Fruticultura, Jaboticabal, v.25, n.2, p.235-241, ago. 2003.
• FIORAVANÇO, J.C. A cultura da pereira no Brasil: situação econômica e entraves para o seu crescimento. Informações Econômicas, São Paulo, v.37, n.3, p.52-60, mar. 2007.
• LEÃO, M.G.A.; MARQUES JÚNIOR, J.; SOUZA, Z.M.; SIQUEIRA, D.S.; PEREIRA, G.T. O relevo na interpretação da variabilidade espacial dos teores de nutrientes em folha de citros. Revista Brasileira de Engenharia Agrícola e Ambiental, Campina Grande, v.14, n.11, p.1152-1159, nov. 2010.
• MIRANDA, N.O.; OLIVEIRA, T.S.; LEVIEN, S.L.A.; SOUZA, E.R. Variabilidade espacial da qualidade de frutos de melão em áreas fertirrigadas. Horticultura Brasileira, Brasília, v.23, n.2, p. 242-249, abr-jun. 2005.
• MOLIN, J.P; MASCARIN, L.S. Colheita de citros e obtenção de dados para mapeamento da produtividade. Revista Engenharia Agrícola, Jaboticabal, v.27, n.1, p.259-266, jan-abr. 2007.
• PERES, L. E. P. Nutrição mineral. 2002. Disponível em: www.biologia.ufc.br/backup/ docentes/JoaoLuiz/Nutricao.Mineral.de.Plantas.pdf NutricaoMineral.pdf> Acesso em: 15 jan. 2011.
• PERRY, E.M.; DEZZANI, R.J.; SEAVERT, C.F.; PIERCE, F.J. Spatial variation in tree characteristics and yield in a pear orchard. Precision Agriculture, Dordrecht, v.11, n.1, p.42-60, 2009.
• PIMENTEL-GOMES, F.; GARCIA, C.H. 2002. Estatística aplicada a experimentos agronômicos e florestais. Piracicaba: FEALQ, 2002. 309 p.
• SILVA, F.M.; SOUZA, Z.M.; FIGUEIREDO, C.A.P.; MARQUES JÚNIOR, J.; MACHADO, R.V. Variabilidade espacial de atributos químicos e de produtividade na cultura do café. Ciência Rural, Santa Maria, v.37, n.2, p.401-407, mar-abr. 2007.
• SIQUEIRA, D.S.; MARQUES JÚNIOR, J.; PEREIRA, G. T. The use of landforms to predict the variability of soil and orange attributes. Geoderma, Amsterdam, v.155, n.1-2, p.55-66, 2010.
• TISSEYRE, B.; MCBRATNEY, A.B. A technical opportunity index based on mathematical morphology for site-specific management: an application to viticulture. Precision Agriculture, Dordrecht, v.9, n.1-2, p.101-113, 2008.
• VALLADARES, G.S.; CAMARGO, O.A.; AZEVEDO, E.C. ; SILVA, A.M.C. Variabilidade espacial e disponibilidade de zinco em solos de vinhedos e adjacências da região vitícola de Jundiaí, SP. In: SIMPÓSIO BRASILEIRO DE SENSORIAMENTO REMOTO, 13., 2007, Florianópolis. Anais… Florianópolis: INPE, 2007. p.449-456.

Publication Dates

• Publication in this collection
  16 July 2012
• Date of issue
  Apr 2012

History

• Received
  16 May 2011
• Accepted
  09 Nov 2011