Sampling grids used to characterise the spatial variability of pH, Ca, Mg and V% in Oxisols

Malhas amostrais utilizadas na caracterização da variabilidade espacial de pH, Ca, Mg e V% em Latossolos

SOIL SCIENCE

ARTIGO CIENTÍFICO

Sampling grids used to characterise the spatial variability of pH, Ca, Mg and V% in Oxisols¹

Malhas amostrais utilizadas na caracterização da variabilidade espacial de pH, Ca, Mg e V% em Latossolos

Maurício Roberto CherubinI, ^*; Antônio Luis Santi^II; Mateus Tonini Eitelwein^III; Clovis Orlando Da Ros^II; Mateus Bortoluzi Bisognin^IV

^IPrograma de Pós-Graduação em Solos e Nutrição de Plantas, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo-ESALQ/USP, Piracicaba-SP, Brasil, [e-mail]

^IIDepartamento de Ciências Agronômicas e Ambientais, Centro de Educação Superior Norte do Rio Grande do Sul, Universidade Federal de Santa Maria, Frederico Westphalen-RS, Brasil, [e-mail], [e-mail]

^IIIPrograma de Pós-Graduação em Engenharia de Sistemas Agrícolas, ESALQ/USP, Piracicaba-SP, Brasil, [e-mail]

^IVGraduando em Agronomia, Centro de Educação Superior Norte do Rio Grande do Sul, Universidade Federal de Santa

ABSTRACT

Knowledge of spatial variability is an important factor to be considered in planning a program of soil sampling and crop management under precision agriculture (PA). In this context, the aim of this work was to evaluate the efficiency of the dimensions of sampling grids used in the state of Rio Grande do Sul (RS), Brazil to characterise the spatial variability of the attributes pH_water, base saturation (V%), calcium (Ca) and magnesium (Mg) levels. The study was carried out on 30 agricultural sites located in the northern region of RS, having soils classified as Oxisols and managed using the tools of PA. The dimensions of the grids under study were: 100 x 100 m (10 areas), 142 x 142 m (10 areas) and 173 x 173 m (10 areas). Soil was collected at a depth of 0.00 to 0.10 m. The data for pH_water, V%, Ca and Mg were subjected to exploratory statistical analysis and to geostatistical analysis by means of semivariograms. The areas showed high Ca (>4.0 cmol_c dm^-3) and Mg (>1.0 cmol_c dm^-3) levels and localised problems of soil acidity (pH_water <5.5 or V<65%), justifying the carrying out of liming at specific sites. For the geostatistical procedures, the sample grids used at the sites of the Oxisols managed under PA in RS are not efficient in capturing the scales of spatial variability of the attributes pH_water, V%, Ca and Mg, which could compromise the accuracy of corrective prescriptions for specific sites.

Key words: Soil acidity. Precision Agriculture. Soil sampling. Geostatistics.

RESUMO

O conhecimento da variabilidade espacial é um importante fator a ser considerado no planejamento de um programa de amostragem de solo e manejo das culturas na agricultura de precisão (AP). Nesse contexto, o objetivo do trabalho foi avaliar a eficiência das dimensões das malhas amostrais utilizadas no Rio Grande do Sul (RS) na caracterização da variabilidade espacial dos atributos pH_água, saturação por bases (V%), teores de cálcio (Ca) e magnésio (Mg). O estudo foi realizado em 30 áreas agrícolas localizadas na região Norte do RS que apresentam solos classificados como Latossolos Vermelhos e que são manejadas com ferramentas de AP. As dimensões das malhas amostrais estudadas foram: 100 x 100 m (10 áreas), 142 x 142 m (10 áreas) e 173 x 173 m (10 áreas). A profundidade de coleta de solo foi de 0,00-0,10 m. Os dados de pH_água, V%, Ca e Mg foram submetidos à análise estatística exploratória e à análise geoestatística por meio de semivariogramas. As áreas apresentaram elevados teores de Ca (> 4,0 cmol_c dm^-3) e Mg (> 1,0 cmol_c dm^-3) e problemas localizados de acidez (pH_água < 5,5 ou V < 65%), justificando a realização de calagem em sítios específicos. Considerando os procedimentos geoestatísticos, as malhas amostrais utilizadas nas áreas de Latossolos Vermelhos manejados com AP no RS não são eficientes para captar as escalas da variabilidade espacial dos atributos de pH_água, V%, teores de Ca e Mg, podendo comprometer a acurácia das prescrições de corretivos em sítios específicos.

Palavras-chave: Acidez do solo. Agricultura de precisão. Amostragem de solo. Geoestatística.

INTRODUCTION

Oxisols are widely distributed across almost all regions of Brazil, and make up the predominant taxonomic class, occupying 31.61% of the surface of the country (2.69 million km²) (ANJOS et al., 2012). In the state of Rio Grande do Sul (RS), Oxisols occur mainly in the north, a region located in the geomorphic province of the Plateau, formed by a succession of volcanic flows, followed by cooling on the surface of the earth's crust and the formation of predominantly basaltic rocks (STRECK et al., 2008).

These are weathered, deep, well-drained and acidic soils and display high homogeneity in the profile (ANJOS et al., 2012; MONTANARI et al., 2008; STRECK et al., 2008). However, despite being considered homogeneous soils, they may comprise variations in their characteristics over both short and long distances.

Variability of the soil is a result of the interaction of complex pedogenic processes controlled by factors in its formation such as climate, topography, source material and vegetation, and intensified by land-use and management practices (MALLARINO; WITTRY, 2004). Understanding spatial variability in the soil is therefore a factor that should be considered when planning a program of soil sampling and of crop management (KERRY; OLIVER; FROGBROOK, 2010; MALLARINO; WITTRY, 2004; MONTANARI et al., 2008; MONTANARI et al., 2012; NANNI et al., 2011; SOUZA et al., 2004; STEPIEN; GOZDOWSKI; SAMBORSKI, 2013).

With the advent of precision agriculture (PA) in Brazil and especially in the South region, where only in RS state, it is estimated that approximately two million hectares are managed using tools of PA (SANTI et al., 2009), knowledge of the spatial variability of the soil was made possible through systematic sampling by means of regular grids. However, although the use of sampling grids is officially recommended in the South of Brazil (COMISSÃO DE QUÍMICA E FERTILIDADE DO SOLO, 2004), there are still doubts as to the ideal size which would allow capture of the spatial variability of the chemical properties of the soil.

Currently, the sampling grids that are commonly adopted on Brazilian farms have dimensions ranging from 100 to 225 m (NANNI et al., 2011). In the South region, samples are collected at somewhat smaller distances, from 100 to175 m. However, the definition of such dimensions for the grids is primarily based on economic and practical reasons and sometimes neglects the geostatistical principles of spatial dependence (VIEIRA, 2000; WEBSTER; OLIVER, 2007).

The influence that the dimensions of the grids used in collecting soil samples have on the accuracy of the generated information is widely discussed in international literature. Studies involve the use of the principles and tools of geostatistics to direct the sampling of soils (KERRY; OLIVER; FROGBROOK, 2010; WEBSTER; LARK, 2012), the influence of the size of the sampling grid on the quality of data interpolation (KRAVCHENKO 2003; FRANZEN, 2011), the reproducibility of the spatial variability (CHANG et al., 1999; STEPIEN; GOZDOWSKI; SAMBORSKI, 2013) and the management of nutrients at specific sites (MALLARINO; WITTRY, 2004; STEPIEN; GOZDOWSKI; SAMBORSKI, 2013). However, there are few studies on the conditions of Brazilian soils (MONTANARI et al., 2012; NANNI et al., 2011) and none for conditions in the South of Brazil.

In this context, the aim of this work was to evaluate the efficiency of the dimensions of sampling grids used in RS in characterising the spatial variability of the attributes pH_water and V% and Ca and Mg levels in Oxisols under PA.

MATERIAL AND METHODS

The study was carried out in 30 agricultural areas located in 13 municipalities in the north of Rio Grande do Sul, Brazil (Figure 1), covering 2,467.82 ha. The study area is located in the geomorphic province of the Plateau , and extends to the physiographic regions of the Alto Uruguai, Planalto Médio and Missões. The terrain is gently undulating with the predominant soils classified as Oxisols (Latossolos Vermelhos by Brazilian Classification Soil System), having high aluminum and iron levels associated with a low base saturation (alumino ferric, dystroferric or dystrophic) and clayey texture (STRECK et al., 2008).

The management system that has been adopted in the area is direct planting (no-tillage), rotating soybean and maize crops in the summer and white oats, black oat and forage turnip in the winter. Fertilization is carried out in the sowing furrow for each crop, and liming over the total area every five years or so.

The areas were sampled from the dimensions of the grids that are commonly used at field level in the south of Brazil, these being: 100 x 100 m (one sample for each ha); 142 x 142 m (one sample for every two ha); and 173 x 173 m (one sample for every three ha). Thus, for each grid 10 areas were studied, giving a total of 1,461 soil samples. The areas being analysed were sampled for the first time by georeferencing, so they did not suffer the effects of previous variable-rate applications of lime and fertilizer.

Each soil sample was composed of 10 subsamples collected within 10 m of the georeferenced point using a automatic screw auger. The collection depth was from 0.00 to 0.10 m, as recommended by the Commission for Chemistry and Soil Fertility (COMISSÃO DE QUÍMICA E FERTILIDADE DO SOLO, 2004) for agricultural areas under consolidated no-tillage systems. The samples were later forwarded to the laboratory where basic chemical analyses were conducted following methodologies recommended by the Commission for Chemistry and Soil Fertility (2004). The following attributes were evaluated in the study: active acidity (pH_water, 1:1), base saturation (V%) and the cations related to acidity: calcium (Ca, cmol_c dm^-3) and magnesium (Mg, cmol_c dm^-3).

Data were subjected to descriptive statistical analysis using the Statistical Analysis System - SAS 8.0 software (SAS Inc, Cary, USA), with the aim of verifying position and dispersion of the data. The statistical parameters determined were; minimum; mean; maximum; standard deviation and coefficients of variation (CV%), precision (CP%), skewness (Cs) and kurtosis (Ck). Based on the CV values obtained, the dispersion of the data was classified as low (CV < 15%), moderate (CV 15-35%) and high (CV > 35%) (WILDING; DREES, 1983). In addition, the existence of a central tendency (normality) in the original data was noted using the W-test (p<0.05).

Analysis of the spatial variability of the pH_water data, V%, Ca and Mg, was carried out by semivariogram, adjusted by theoretical mathematical models using the Gamma Design Software - GS+. The models for the semivariogram were adjusted based on the best coefficient of determination (r²) and lowest residual sum of squares (RSS), and were evaluated using cross validation. The parameters of the semivariogram were defined from the adjustment of a mathematical model to the data: nugget effect (C₀), contribution (C₁), sill (C₀ + C1) and range (a). The spatial dependence index (SDI) was calculated using the relationship shown in Equation 1.

C₀/(C₀+C₁)*100 (1)

The degree of spatial dependence (DSD), based on the SDI, was classified as strong for SDI < 25%; moderate for an SDI between 25 and 75%, and weak for SDI> 75% (CAMBARDELLA et al. 1994).

RESULTS AND DISCUSSION

Exploratory statistical analysis of the attributes for pH_água and V% (Table 1) made it possible to verify that 27 areas studied (90.0%) present places with problems of soil acidity, pH <5.5 or V <65.0%, which can restrict normal plant development under no-tillage system (COMISSÃO DE QUÍMICA E FERTILIDADE DO SOLO, 2004). This finding is related to the natural acidity characteristic of soils that have undergone intense weathering, and associated with processes of soil acidification triggered by the management practice adopted (ANJOS et al., 2012; STRECK et al., 2008). In agricultural areas, soil acidity is enhanced by absorption of the base cations by crops and their subsequent export at harvest, by inadequate soil management favouring the erosion and exposure of (more acidic) subsurface horizons, by the use of nitrogen fertilizers and by the oxidation of sulphur and organic matter (SOUZA; MIRANDA; OLIVEIRA, 2007).

The use of smaller sampling grids, allowing for an increase in the number of sampling points, resulted in improvements in the representativeness and accuracy of soil sampling conducted in the study areas (<% CP). However, irrespective of the grid dimensions, it was possible to identify, based on the minimum and maximum values for pH_água and V%, a considerable amplitude relative to the mean, demonstrating that the use of georeferenced soil sampling, as recommended under PA, was characterised as an important strategy in identifying places where soil acidity may become a limiting factor. Thus, PA allows an increase in the technical and economic efficiency of the use of limestone, based on interventions located at a specific site.

It was found that for dispersion of the data, based on the coefficient of variation (CV%), the values for pH_water showed low dispersion (CV: <15%), with an amplitude of from 2.06 to 7.76%. Whereas the values for %V ranged from 3.92 to 25.18%, classifying the dispersion as from low to moderate (CV: <35%). These results agree with several studies made into Oxisols (BOTTEGA et al., 2013; CHERUBIN et al., 2011; DALCHIAVON et al., 2012; LIMA; SILVA; SILVA, 2013; MONTANARI et al., 2008; SANTI et al., 2012; SOUZA et al., 2004). However, the low CVs seen for the values of pH_water are due to its logarithmic scale of expression, requiring caution in comparison with other variables (SOUZA; MIRANDA; OLIVEIRA, 2007).

Although the data for pH_water and V% displayed low and moderate dispersion, 56.66 and 26.66% of the areas studied did not follow normal frequency distributions, being confirmed by coefficients of skewness, offset to the left (Cs <0) or to the right (Cs> 0), and by coefficients of kurtosis with platykurtic (Ck <0) and leptokurtic (Ck> 0) distributions. Results like these, that indicate a lack of normality in the data, were also seen in areas of Oxisols by Cavalcante et al. (2007), Cherubin et al. (2011), Dalchiavon et al. (2012) and Montanari et al. (2008). However, it is important to note that normality of the data is not a requirement of geostatistics; but the presence of a skewed distribution having many anomalous values should be considered, since kriging is a linear estimator (WEBSTER; OLIVER, 2007). In geostatistics, more important than the normality of the data is whether or not the so-called proportional effect occurs, in which the mean and variance of the data may not be constant, nor be a function of spatial positioning in the study area (CAVALCANTE et al., 2007).

For the Ca and Mg cations (Table 2), it was found that the Oxisols studied display on average levels considered as high for both of these macronutrients (Ca > 4.0 cmol_c dm^-3 and Mg > 1.0 cmol_c dm^-3) (COMISSÃO DE QUÍMICA E FERTILIDADE DO SOLO, 2004).

In studies conducted by Amado et al. (2009), Cherubin et al. (2011) and Santi et al. (2012) into Oxisol in RS, high levels of Ca and Mg were also found, agreeing with the results of the present work. These high levels are due to the weathering of minerals rich in these elements, such as plagioclase, pyroxene and olivine, in the makeup of the basalt rocks which are the predominant source material of the Oxisol occurring in the study region (STRECK et al., 2008). Furthermore, the addition of dolomitic limestone over the years, associated with the low mobility of the Ca and Mg in the soil, may also have significantly contributed to the maintenance and rise in levels of these macronutrients.

Dispersion of the Ca and Mg data was classified as low to moderate, with CV amplitudes that ranged from 11.94 to 34.08% and from 14.86 to 35.00% respectively, agreeing with Amado et al. (2009), Bottega et al. (2013), Cavalcanti et al. (2007), Cherubin et al. (2011), Lima, Silva and Silva (2013), Santi et al. (2012) and Souza et al. (2004). Normal frequency distributions for Ca and Mg were observed in 66.66 and 50.00% of the areas respectively. However, in areas with non-normal distribution, the data showed no discrepant values that could limit the use of geostatistics.

The spatial variability of the attributes pH_water and V%, and of the bases Ca and Mg, was analysed from the results of the geostatistical analysis (Tables 3 and 4). In the four attributes studied, a large difference in the structure of spatial variability was seen for the areas under study, even when the samples were collected with the same dimensions of sampling grid. This variation can be evidenced by the amplitude observed for the radius of spatial dependence between samples (ranges). For all the grids under study, areas were noted where the values for pH_water, V%, Ca and Mg showed no spatial dependence, featuring random distributions (pure nugget effect), and where geostatistics therefore could not be applied (VIEIRA, 2000). Note that there is a trend towards an increase in the occurrence of random distributions as the size of the sampling grid increases, confirming assumptions of Cambardella et al. (1994), Vieira (2000) and Webster and Oliver (2007), in which the sampling grid used does not have enough points to detect the dependence, if any, which will be expressed at distances which are smaller than the smallest spacing between samples.

In areas where the attributes pH_water, V%, Ca and Mg present a defined structure of spatial variability, it was generally found that the degree of spatial dependency was moderate to strong, showing that under such conditions, these attributes were more affected by intrinsic properties of the soil (CAMBARDELLA et al., 1994). Furthermore, management practices, such as the uniform surface liming carried out over the period of agricultural land use, contribute to the low spatial variability of the attributes of acidity over shorter distances.

These results however warrant caution, since the semivariograms were generated using a limited number of points (n). Webster and Lark (2012) explain that the use of 30 to 50 pairs of points over each distance, as indicated by some authors (LANDIM, 2006), means having less than 50 sampling points in two-dimensional grids and almost inevitably leads to estimates which are only slightly accurate and to semivariograms with high levels of error. According to these authors, for reliable estimates to be generated when the variation is isotropic, at least 100 sampling points are required and ideally 150 to 200.

Moreover, although the number of sampling points is considered to be important in the literature, it is not subject to generalisation for all situations, and may be too small in areas where the variables display their variability on smaller scales (short distances), and too large for areas where their variability is presented on larger scales (long distances) (KERRY; OLIVER; FROGBROOK, 2010).

When larger sampling grids are used therefore, with more widely-spaced points (usually ≥ 100 m), together with a reduced size for the areas (a situation commonly seen in the areas under PA in southern Brazil) the comprehension and representation of the different scales of variability for the attributes of acidity and related cations (Ca and Mg) are still limited, even with the use of geostatistics. It can therefore be inferred that the sampling grids used in the Oxisol areas in RS managed under PA are not efficient in capturing the different scales of spatial variability for pH_water, V%, Ca and Mg, especially when expressed over short distances. However, it is important to note that the insertion of PA and consequently of soil sampling with grids, even if displaying limited levels of reliability in some situations, has enabled progress to be made in diagnosing soil fertility which is unprecedented in the agriculture of Rio Grande do Sul and in the rest of Brazil, redeeming the role of sampling in the management of soil fertility and agricultural production.

The inefficiency of sampling grids used commercially is also pointed out by Franzen (2011), Kerry and Oliver (2008) and Kerry, Oliver and Frogbrook (2010) in the United States, where generally grids of 100 m (one sample per hectare) are employed in areas under PA. For Brazilian conditions, the results of the present study agree with those obtained by Corá and Beraldo (2006) and Nanni et al. (2011) in which the collection of one sample per hectare was not sufficient to capture the spatial variability of V%, with the use of denser grids being necessary in order to make more accurate recommendations for soil correctors possible.

Thus, in the absence of spatial dependence of the data, or with dependencies of limited reliability, it is suggested that users of PA exercise caution when employing kriging interpolation, otherwise the thematic maps of soil attributes, used to determine the variable rates of correctors and fertilizers, will not express the main patterns of spatial variability present in the area (KERRY; OLIVER, 2008). Under such conditions, the use of simple interpolation (linear or polynomial) such as inverse distance, can be just as effective as the kriging method (KRAVCHENKO, 2003; FRANZEN, 2011). According to Corá and Beraldo (2006), map accuracy is dependent on the interpolation method used to estimate the values for non-sampled locations, and in turn, the interpolation method is dependent on the density of sampled points per area.

These results endorse what is currently being practised in the field, where users of PA, during the preparation process for thematic maps, employ computer software available in the market, but which do not take into account the spatial dependence of the attributes being analysed in estimating values at non-sampled locations (CORÁ; BERALDO, 2006). However, it is important to state again that the ability to generate thematic maps for the attributes pH_water, V%, Ca and Mg in the absence of any spatial dependence of the data verified by geostatistics, does not exempt the planning of sampling from recommending denser sampling, which would be capable of actually detecting the different scales of spatial variability of the analysed attributes. As good as the interpolation method is, it will never be able to predict values with accuracy, compared to the values obtained by sampling in the field.

Given the results of this exploratory study, the assumptions already outlined by Chang et al. (1999) and Stepien, Gozdowski and Samborski (2013) can be confirmed, in which the spatial variability of the attributes pH_water, V%, Ca and Mg is unique for each area, being conditioned at different scales by intrinsic and extrinsic soil factors. Therefore, from sampling grids with dimensions ≥ 100 m, such as those that have been used in the south of Brazil, it is not possible to generalise a reliable model of the spatial variability of these attributes for Oxisols. However, the study does give an overview of the scale of variation for the soil attributes under study, and becomes an important reference that can assist in planning future soil sampling strategies to be adopted in areas under PA in RS.

A suggestion for the future would be carrying out detailed studies in a microregion, farm or field scale to test different dimensions for sampling grids, with a view to better understanding the variability of the different soil chemical attributes, and the possible definition of more efficient strategies for fertility management.

CONCLUSIONS

1. The sampling grids used in Oxisols under PA in RS, generally employing geostatistical procedures, are not efficient in capturing the scales of spatial variability of the attributes pH_water, V%, Ca and Mg, and may lead less accurate liming;

2. Spatial variability over short distances should be taken into consideration in future plans for soil sampling, aiming for the characterization and the site specific management of soil acidity in Oxisols areas managed under PA.

Recebido para publicação em 26/03/2013

Aprovado em 24/06/2014

• AMADO, T. J. C. et al Atributos químicos e físicos de Latossolos e sua relação com os rendimentos de milho e feijão irrigados. Revista Brasileira de Ciência do Solo, v. 33, n. 4, p. 831-843, 2009.
• ANJOS, L. H. C. dos. et al Sistema Brasileiro de Classificação de Solo. In: KER, J. C. et al (Ed.). Pedologia: Fundamentos. Viçosa-MG: Sociedade Brasileira de Ciência do Solo, 2012. p. 303-343.
• BOTTEGA, E. L. et al. Variabilidade espacial de atributos do solo em sistema de semeadura direta com rotação de culturas no cerrado brasileiro. Revista Ciência Agronômica, v. 44, n. 1, p. 1-9, 2013.
• CAMBARDELLA, C. A. et al Field-scale variability of soil properties in central Iowa soils. Soil Science Society American Journal, v. 58, n. 5, p. 1501-1511, 1994.
• CAVALCANTE, E. G. S. et al Variabilidade espacial de atributos químicos do solo sob diferentes usos e manejos. Revista Brasileira de Ciência do Solo, v. 31, n. 6, p. 1329-1339, 2007.
• CHANG, J. et al Precision Farming Protocols: Part 1. Grid distance and soil nutrient impact on the reproducibility of spatial variability measurements. Precision Agriculture, v. 1, n. 3, p. 277-289, 1999.
• CHERUBIN, M. R. et al Caracterização e estratégias de manejo da variabilidade espacial dos atributos químicos do solo utilizando a análise dos componentes principais. Enciclopédia Biosfera, v. 7, n. 13, p. 196-210, 2011.
• COMISSÃO DE QUÍMICA E FERTILIDADE DO SOLO. Manual de adubação e calagem para os estados do Rio Grande do Sul e Santa Catarina Porto Alegre: Sociedade Brasileira de Ciência do Solo - Núcleo Regional Sul, 2004. 400 p.
• CORÁ, J. E.; BERALDO, J. M. G. Variabilidade espacial de atributos do solo antes e após calagem e fosfatagem em doses variadas na cultura de cana-de-açúcar. Engenharia Agrícola, v. 26, n. 2, p. 374-387, 2006.
• DALCHIAVON, F. C. et al Variabilidade espacial de atributos da fertilidade de um Latossolo Vermelho Distroférrico sob Sistema Plantio Direto. Revista Ciência Agronômica, v. 43, n. 3, p. 453-461, 2012.
• FRANZEN, D. W. Collecting and analyzing soil spatial information using kriging and inverse distance. In: CLAY, D. E.; SHANAHAN, J. F. GIS Applications in Agriculture: Nutrients Management for Energy Efficiency, USA: CRC Press, 2011. p. 61-80.
• KERRY, R.; OLIVER, M. A. Determining nugget:sill ratios of standardized variograms from aerial photographs to krige sparse soil data. Precision Agriculture, v. 9, n. 1/2, p. 33-56, 2008.
• KERRY, R.; OLIVER, M. A.; FROGBROOK, Z. L. Sampling in Precision Agriculture. In: OLIVER, M. A. (Org.). Geostatistical Applications for Precision Agriculture. Heidelberg: Springer-Verlag, 2010. p. 35-63.
• KRAVCHENKO, A. N. Influence of spatial structure on accuracy of interpolation methods. Soil Science Society of America Journal, v. 67, n. 5, p. 1564-1571, 2003.
• LANDIM, P. M. B. Sobre Geoestatística e mapas. TERRÆ DIDATICA, v. 2, n. 1, p. 19-33, 2006.
• LIMA, J. S. S.; SILVA, S. A.; SILVA, J. M. Variabilidade espacial de atributos químicos de um Latossolo Vermelho-Amarelo cultivado em plantio direto. Revista Ciência Agronômica, v. 44, n. 1, p. 16-23, 2013.
• MALLARINO, A. P.; WITTRY, D. Efficacy of grid and zone soil sampling approaches for site-specific assessment of phosphorus, potassium, pH, and organic matter. Precision Agriculture, v. 5, n. 2, p. 131-144, 2004.
• MONTANARI, R. et al Variabilidade espacial de atributos químicos em Latossolo e Argissolos. Ciência Rural, v. 38, n. 5, p. 1266-1272, 2008.
• MONTANARI, R. et al The use of scaled semivariograms to plan soil sampling in sugarcane fields. Precision Agriculture, v. 13, n. 5, p. 542-552, 2012.
• NANNI, M. R. et al Optimum size in grid soil sampling for variable rate application in site-specific management. Scientia Agrícola, v. 68, n. 3, p. 386-392, 2011.
• SANTI, A. L. et al Análise de componentes principais de atributos químicos e físicos do solo limitante à produtividade de grãos. Pesquisa Agropecuária Brasileira, v. 47, n. 9, p. 1346-1357, 2012.
• SANTI, A. L. et al É chegada a hora da integração do conhecimento. Revista Plantio Direto, v. 129, n. 1, p. 24-30, 2009.
• SOUZA, D. M. G.; MIRANDA, L. N.; OLIVEIRA, S. A. Acidez do solo e sua correção. In: NOVAIS, R. F. et al. (Ed.). Fertilidade do Solo Viçosa-MG: Sociedade Brasileira de Ciência do Solo, 2007. p. 206-274.
• SOUZA, Z. M. de. et al Variabilidade espacial do pH, Ca, Mg e V% do solo em diferentes formas do relevo sob cultivo de cana-de-açúcar. Ciência Rural, v. 34, n. 6, p. 1763-1771, 2004.
• STEPIEN, M.; GOZDOWSKI, D.; SAMBORSKI, S. A case study on the estimation accuracy of soil properties and fertilizer rates for different soil-sampling grids. Journal of Plant Nutrition and Soil Science, v. 176, n. 1, p. 57-68, 2013.
• STRECK, E. V. et al Solos do Rio Grande do Sul 2. ed. Porto Alegre: EMATER/RS-ASCAR, 2008. 222 p.
• VIEIRA, S. R. Geoestatística em estudos de variabilidade espacial do solo. In: NOVAIS, R. F.; ALVARES V., V. H.; SCHAEFFER, C. E. G. R. (Ed.). Tópicos em Ciência do Solo. Sociedade Brasileira de Ciência do Solo, p. 1-54, 2000. v. 1.
• WEBSTER, R.; LARK, M. Field Sampling for Environmental Science and Management London: Routledge, 2012. 200 p.
• WEBSTER, R.; OLIVER, M. A. Geostatistics for Environmental Scientists 2. ed. Chichester: John Wiley & Sons Ltd. 2007. 330 p.
• WILDING, L. P.; DREES, L. R. Spatial variability and pedology. In: WILDING, L. P.; SMECK, N. E.; HALL, G. F. (Ed). Pedogenesis and Soil Taxonomy: Concepts and Interactions. Amsterdam: Elsevier, 1983. p. 83-116.

*

Autor para correspondência

1

Parte da Dissertação de Mestrado do primeiro autor, apresentada ao Programa de Pós-Graduação em Agronomia: Agricultura e Ambiente da Universidade de Santa Maria, campus Frederico Westphalen, RS

Publication Dates

• Publication in this collection
  26 Aug 2014
• Date of issue
  Dec 2014

History

• Received
  26 Mar 2013
• Accepted
  24 June 2014