SciELO - Scientific Electronic Library Online

vol.86 número209Corporate governance practices in brazilian family construction companiesAssessing the effect of adding maize and rice brans in the development by twin-screw extrusion of a ready-to-eat cereal formulated with flours of quality protein maize and zinc biofortified rice índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados



Links relacionados

  • Em processo de indexaçãoCitado por Google
  • Não possue artigos similaresSimilares em SciELO
  • Em processo de indexaçãoSimilares em Google



versão impressa ISSN 0012-7353

Dyna rev.fac.nac.minas vol.86 no.209 Medellín abr./jun. 2019 


Variabilidad espacial del pH del suelo muestreado por dos metodologías utilizadas en agricultura de precisión en granjas bajo rotación de cultivos

Spatial variability of soil pH sampled by two methodologies used in precision agriculture in farms under crop rotation

Gabriel Araújo e Silva Ferraz a   , Brenon Diennevan Souza Barbosa a   , Étore Francisco Reynaldo b   , Sthéfany Airane dos Santos a   , Jose Roberto Moreira Ribeiro Gonçalves c   , Patrícia Ferreira Ponciano Ferraz a  

a Department of Agricultural Engineering, Federal University of Lavras, Lavras, Brazil.,,,

b Specialist in Equipment at Monsanto, Brazil.

c University Center Augusto Motta, Rio de Janeiro, Brazil.


El objetivo de este estudio fue caracterizar la variabilidad espacial del pH para suelos de dos fincas en el estado de Paraná, Brasil, con base en dos métodos diferentes de muestreo utilizados en agricultura de precisión, mediante análisis geoestaddísticos. El primero método de muestreo del pH consistió en la colecta de muestras de suelo por ele método tradicional (1punto/ha). El segundo método de determinación del pH fue el sensor de suelo dinâmico (200 puntos / ha). El modelo esférico fue más adecuado para la mayoría de los semivariogramas, independientemente del método de muestreo. La mejor distribución espacial fue obtenida por la cualidad del muestreado del sensor de suelo dinâmico. El número de muestras del pH colectadas y el método de muestreo influyeron en la representación visual de la variabilidad del pH.

Palabras-clave: variabilidad espacial; sensor de contacto móvil; muestreo en cuadrícula; geoestadística


This study aimed to characterize the spatial variability of pH in soils of two farms in the state of Paraná, Brazil, based on two different sampling methods used in precision agriculture, by means of geostatistical analyzes. The first method of sampling the pH grid consisted in the collection of soil samples by the traditional method (1 point / ha). The second method of pH determination was by on-the-go soil sensor (200 points / ha). The spherical model was better suited to most semivariograms, regardless of the sampling method. After adjusting the semivariograms for soil pH determination methods, thematic maps were made using normal kriging. The best spatial distribution of pH was obtained where the attribute was sampled by the on-the-go sensor. The number of pH samples collected and the sampling method influenced the visual representation of pH variability.

Key words: spatial variability; mobile contact sensor; grid sampling; geostatistics

1. Introduction

Precision agriculture (PA) can be defined as the use of agricultural techniques based on information technologies for the treatment of spatial variability of soil properties and characteristics from determined area [1].

The rising demand for variability management of soil and plant properties in the field is allied with availability and adoption of PA tools and technology. Together with geospatial tools, using global navigation satellite systems (GNSS) and geographic information systems (GIS), it becomes a highly reliable tool for the spatial characterization of soil properties [17]. However, there is still an obstacle about techniques for data collection on soil properties in a quick, efficient, and accurate way [25].

The method traditionally used for soil sampling may become unviable due to the expensive field work and laboratory as a function of the high number of samples. Moreover, the traditional method is negatively impacted by the higher cost and time to obtain the results, and hence in the subsequent creation of reliable thematic maps, becoming necessary to seek more efficient and advanced methods for collecting and processing these data [15,2].

According to Mulla [18], the future will require massive data collections and analysis on a scale considered not feasible at present, involving stationary or mobile sensors that can accurately measure different plant characteristics in real time. Sensors capable of mapping some soil properties are already available to the farmer. This type of equipment allows performing soil pH sensing in high resolution through antimony ion-selective electrodes [22].

Aiming at this cost saving and searching for reliable results, geostatistics is a technique that works on the characterization of the spatial dependence or not of the elements under analysis, when supplied with a mass of data with aggregated geospatial information. In PA, this technique becomes very useful in planning and mapping the variability of soil physical and chemical properties, interpolating values of these attributes to non-sampled sites, using a mathematical model named kriging [20,2], enabling the creation of thematic maps that characterize the spatial distribution of properties in an area, becoming a basis for the farmer to make decisions on soil management and farming.

The aim of this study was to characterize the magnitude of the spatial variability of soil pH sampled by two different methods used in PA, in order to fit semivariograms from different methods and models, thus obtaining isochore maps interpolated by kriging more reliable, being possible to observe differences among the sampling methodologies.

2. Material and methods

The experiment was performed in the municipalities of Guarapuava and Cantagalo, both in the state of Paraná, located at the geographic coordinates 25° 31’ 22.60” S and 51° 30’ 16.46" W for the Cupim Farm (A1) with an area of 60 ha, and coordinates 25° 16’ 52.25” S 52° 6’ 12.89" W for the Juquiá Farm (A2), with an area of 20 ha. The climate of the region is Cfa, according to Köppen classification. The areas are cultivated under crop rotation of soybeans, oats and maize. The soil of the study areas is characterized as inceptisol, with prominent A horizon, smooth undulating relief and basalt substrate [11], with a textural class varying from clayey to very clayey.

For the determination of soil pH values for both areas, two different methods were used in relation to their determination principles: i) collection in grid sampling, ii) collection through on-the-go soil sensing.

Grid sampling of one point per hectare (100 x 100 m) were developed in both study areas in order to collect the pH, which were correctly georeferenced with the aid of a GNSS signal receiver with high-precision differential correction. Each sample was withdrawn at a depth of 0-20 cm using an soil auger. At each georeferenced point, 11 to 15 subsamples were collected randomly at a radius of 3 m from the center point. These subsamples were homogenized properly generating a single sample representative of the sampling point. The soil samples were sent to the Coodetec soil analysis laboratory, where the necessary chemical analyses were performed.

The second sampling method, also performed in both areas, used an on-the-go sensor for soil pH of the Veris PMC® platform ("P" - pH, "M" - organic matter and "C" - apparent electrical conductivity).

Schirrmann et al. ([22]) describe Soil pH Manager™ which is a module of the Mobile Verisme Platform (SMP) as a sensor that measures soil pH in direct contact with the samples when the equipment is in motion. This sensor is divided into three parts: i) hydraulic system for routing of soil samples, ii) pH electrode measurement system, iii) washing system with distilled water.

During the course of sampling, the hydraulic system is activated so that a metal structure penetrates the soil to a defined foot depth for sample removal. The collected samples are suspended by the hydraulic actuator so that they come in contact with the antimony electrodes. The up and down movement of the soil sampler (hydraulic system) is driven by a proximity sensor. PH is measured in untreated and naturally moist samples. The pH value of the samples is calculated as a function of the average voltage that the electrodes emit to the system that performs the voltage conversion to pH values.

The measurement time of each sample is on average from 7 to 25 s. After pH measurement the electrodes are washed with distilled water stored in the upper part of the equipment by two nozzles located next to each electrode.

According to Schirrmann et al. [22] and in accordance with Veris Technology's calibration recommendations, calibration of the sensor's antimony electrodes is done by a two-point calibration using standard solutions of pH 4 and 7 at the beginning of each experiment. In this calibration the electrodes were inserted into the buffer solution and it was gently shaken for 30 sec to a steady state. The calibration process lasts around 2 min.

The sensor was tractioned by a tractor with an average forward speed of 9 km.h-1, which ranged 30 m spaced transects in the entire area. The density observed for the use of sensor averaged 200 points per hectare. The tractor also had an integrated satellite navigation system, which provided a correct parallelism among transects and the georeferencing of the points sampled by the equipment.

In order to characterize the spatial dependence of the soil pH in both study areas referring to the two soil sampling methodologies, geostatistical analyses were performed by classic and robust semivariograms. The classic semivariogram was estimated according to eq. (1):

Where: N (h) is the number of experimental pairs of observations Z(xi) and Z(xi + h) separated by a distance h. The semivariogram is represented by the graph ("γ ") ̂"(h)" versus h. From the fit of a mathematical model to the calculated semivariation values (("γ “) ̂"(h)"), the coefficients of the theoretical model were estimated for the semivariogram called: nugget effect (C0), sill (C0 + C1), and range (a), as described by Vieira et al. ([24]).

According to Cressie and Hawkins [7], the robust estimator of semivariogram values is less susceptible to the influence of mass data values than the classic estimator. Thus, the robust estimator is described by eq. (2).

The robust estimator assumes that the differences Z(s+h) - Z(s) are distributed normally for all pairs (s+h, s). The transformation of the square root of differences is shown as having moments similar to those from the normal distribution and the denominator of the equation is the bias correction [9].

The fitting semivariogram models were selected according to the ordinary least squares (OLS) and weighted least squares (WLS) estimated by classic and robust mode and maximum likelihood (ML) and restricted maximum likelihood (REML) using the classic estimator. For all methods, the spherical, exponential and Gaussian models will be tested, totaling 18 semivariograms for each variable. For the choice of semivariogram fitting methods and models, data validation will be considered (Fig. 1) [8,9], assuming the stationarity of the intrinsic hypothesis of omnidirectional behavior. [2]

Validation is an error estimation technique that compare predicted values with sampled ones [12]. The sample value, at a certain location Z(si), is temporarily discarded from the data set, and then a kriging prediction is performed on the location, using the remaining samples.

Some values will be very useful for the best method choice, such as: Mean error (ME), the standard deviation of mean error (SDME), reduced mean error (RE) and the standard deviation of reduced mean errors (SDRE). Eq. (3) represents the mean error by validation (ME):

Where n is the data number, Z(si) is the value observed at point si, and is the value predicted by ordinary kriging at point si, excluding the observation Z(si) (Faraco et al., 2008).

The reduced mean error (RE) is defined by eq.(4) :

Where is the standard deviation of kriging at point si, excluding the observation Z(si). The standard deviation of reduced errors (SDRE) is obtained from Eq. (5) :

According to Cressie [6] and McBratney and Webster [16], the models can be evaluated by reduced mean error (ER), standard deviation of mean errors (SDME) and standard deviation of reduced errors (SDRE).

As closer to zero the average difference among the values, the best the estimator. The selection criteria based on validation should find the value of ME and RE closest to zero, the value of SDME should be the lowest, and the value of SDRE should be closest to one.

The choice of the best theoretical fitting model of empirical semivariogram of pH for each area and for each sampling method was performed, considering some criteria based on validation, such as the mean error (ME) and the reduced mean error (RE) closest to zero, the value of standard deviation of mean error (SDME) should be the lowest among the methods, and the value of standard deviation of reduced mean error (SDRE) should be the closest to one. Thus, the fitting model and method of semivariogram that met the highest number of validation requirements would be used. However, if this is not met, the model and the method where the sill (the sum of the nugget effect plus the contribution) is closest to the variance will be used as the decision criterion ([10]). According to Webster and Oliver [27], the sill value must be close to the value of data variance if the semivariogram cloud has a sill.

After fitting the mathematical models of semivariogram and the conjugation that promoted the best validation statistics, the data interpolation observed by ordinary kriging with an interpolation radius equivalent to the range of the fitted semivariogram was performed in order to allow prediction of pH values at not sampled locations [13]. For modeling purposes of the experimental semivariograms, an omnidirectional pattern was assumed.

The georeferencing of the area was performed based on the geographic coordinates obtained by a differential GNSS signal receiver in the demarcation phase of sampling points. The data processing for the geostatistical analysis was performed in the R statistical software, through the geoR library [19], and isochore maps were generated in ArcGIS 10.0 software.

3. Results and discussion

The results referring to the descriptive statistics applied to the pH data collected through the two different methodologies are described in Table 1. A difference/amplitude of the maximum and minimum pH values is observed through on-the-go soil sensor in the order of 1.6 pH unit (5.2-6.8) referring to sampling for area A1. In contrast, this difference was 2.25 pH units (5.15-7.40) for area A2.

Based on the analysis of descriptive statistics presented in Table 1, it can be observed that the averages of pH values sampled by grids and sampled by sensor in both areas were different.

Table 1. Descriptive statistics of pH sampled by grid (pH grid) sampling and by on-the-go soil sensor (pH sensor 

Min - minimum value of the variable, Max - Maximum value of the variable, Av - Averange, Me - Median, SD - Standard deviation, Var - Variance, C.V. - Coefficient of variation, k- Coefficient of kurtosis, Asy - Asymmetry.

Source: The Authors.

This can be explained mainly because the on-the- go soil sensor obtains a larger number of samples than the traditional method. The pH values observed by the sensor were higher than the values obtained by the traditional method in area A2.

Chirmman et al. (2011) reported an estimate of pH values by on-the-go soil sensor higher than the data from the laboratory for soils of Germany. This sensor behavior was not observed for area A1.

The results from the CV analysis of both pH sampling methods were lower than 10%, so that the scatter of values around the mean is considered homogeneous, according to Warrick and Nielsen [26]. However, it should be emphasized that the occurrence of highest and lowest values of pH in the areas was not characterized spatially, being necessary to use a geoestatistical analysis of these values.

Based on the methodology for determination of geostatistical analyses, it was possible determine the parameters (C0, C0+C1, a) of semivariograms (Tables 2,3) fitted by different methods and models for characterization of spatial dependence of pH obtained by different sampling methods in both study areas.

Table 2 Methods, models and estimated parameters of experimental semivariograms for the pH sampled by grid and by sensor in the Cupim Farm (A1) 

C0 - Nugget effect, C1 - Contribution, C0+C1 - Sill, a - range, a'- practical range, ME - Mean error, SDME - Standard deviation of mean error, RE - Reduced mean error, SDRE - Standard deviation of reduced mean error, WLS - C - Weighted least squares classic estimator, ML - Maximum likelihood, OLS - R - Ordinary least squares robust estimator, REML - Restricted maximum likelihood, Sph - Spherical, Exp - Exponential.

Source: The Authors.

Table 3 Methods, models and estimated parameters of experimental semivariograms for the pH sampled by grid and by sensor in the Juquiá Farm (A2) 

C0 - Nugget effect, C1 - Contribution, C0+C1 - Sill, a - range, a'- practical range, ME - Mean error, SDME - Standard deviation of mean error, RE - Reduced mean error, SDRE - Standard deviation of reduced mean error, WLS - C - Weighted least squares classic estimator, ML - Maximum likelihood, OLS - R - Ordinary least squares robust estimator, REML - Restricted maximum likelihood, Sph - Spherical, Exp - Exponential.

Source: The Authors.

Based on the validation criteria, for the pH data referring grid sampling method (Tables 3,, semivariograms (Fig. 1 a , 2 a) were obtained for both areas fitted using the robust WLS method, indicating the spherical model (Fig. 1 a ) and exponential model (Fig. 2b). After the theoretical fitting, it was possible to observe that the range (a) determined for both areas ranged from 92 to 100 m. The range (a) is important to determine the spatial dependency threshold of the attribute with other points, which can also be comprehended as the threshold distance at which a sampled point correlates with the points around it, indicating that points sampled in a radius equal to the range tend to be more homogeneous [6,20].

Source: The Authors

Figure 1 Semivariograms of pH sampled by grid (a) and sampled by on-the-go soil sensor (b) on the Cupim Farm area (A1). 

For the sampling method by on-the-go soil sensor, the experimental semivariograms allowed fitting through the REML (Table 2) and OLS method (Table 3), and it was evidenced that the spherical model was satisfactory according to validation criteria for both studied areas (Figs. 2 b , 3b). However, the ranges (a) found for semivariograms referring to pH sampling by sensor (Figs. 2b, 3 b) were 81 m and 364 m, respectively. According to Souza et al. [23], changes in range values (a) may be related to management techniques applied in soil, which may be distinct for different areas.

Source: The Authors.

Figure 2 Semivariograms of pH sampled by grid (a) and sampled by on-the-go soil sensor (b) on the Juquiá Farm area (A2). 

The semivariograms referring to grid sampling method (Figs. 2 a-3 a) show C0 values equal to 0 (Table 2). However, for semivariograms representing the semivariance of pH data collected by sensor (Tables 2, 3), the values are equal to 0.065 and 0.183 for areas A1 and A2, respectively. According to McBratney and Webster (1986), the nugget effect (C0) is a semivariogram parameter that can indicate how much the data variance cannot be explained for a distance smaller than the radius among samples and for measuring errors. However, it is not possible to detect which of these factors has greater weight in the discontinuity of variance values.

Sana et al. ([21]) reported that determining the range of semivariogram fitted according to a theoretical model is essential, since it allows interpreting how far the prediction of a value is influenced by a sampled point in another nearby place. These authors found in their study a low CV for the traditional sampling method of pH, for a grid sample of 0.64 points per hectare, and later generated pH maps through kriging, observing that the best fitting model of semivariogram for pH data was the spherical model, which corroborates with the results obtained in this study for the traditional method of pH determination.

Cherubin et al. ([3]) reported that a reliable fitting model cannot be generalized to evaluate the variability of attributes for a latosol in Rio Grande do Sul, Brazil, because the variability of each attribute is conditioned to intrinsic and extrinsic factors from each analyzed area, in which the semivariogram referring to each analyzed attribute and area can be fitted by different theoretical models.

In this study, it was observed that the range of semivariogram was changed, with higher values due to the determination method of pH values (Figs. 1, 2, Tables 2, 3). According to Corá et al. ([4]), the prediction quality of values in non-sampled points is strongly influenced by range values, and the ordinary kriging results are more reliable when higher range values are adopted, approaching the field reality on thematic maps referring to the spatial distribution.

After the best method and fitting model of each semivariogram was chosen, the ordinary kriging of pH values was performed, making it possible to generate thematic maps of the spatial distribution of pH sampled by grid and sampled by on-the-go soil sensor. To provide a visual comparison of maps, seven division classes of pH values were used for both study areas (Figs. 3-4).

In Figs. 3 a,4 a, the maps of the spatial variability of pH on area A1 are observed, where it is noted a difference in pH values between the soil sampling methods performed in the area. Observing Fig. 3a, it can be observed that the pH values of the area ranged from 5.5 to 6.0, and the map is mostly darkest gray (whose pH value is 5.5), totaling an area of 18.42 ha. By analyzing Fig. 3b, it is noted that the map represents spatial variability greater than the map of Fig. 3a, and that the value of 6 pH units corresponds to a significant area of the map (11.49 ha).

Source: The Authors.

Figure 3: Grid sampling (a) and on-the-go soil sensor (b) of the area belonging to the Cupim Farm (A1). 

Figs. 4 a, 4 b show the spatial variability maps of the pH from A2 referring sampling methods by grid and by sensor. From the analysis of Fig. 4 a , it can be realized that a large part of the area had pH values equal to 5.5 (darkest gray), representing an area equal to 39.51 ha. In Fig. 4 b , it was observed that a large part of the terrain showed light grey color (pH of 6.5), totaling an area of 29.28 ha. By Fig. 4 b , it is possible to observe a greater variation of pH in relation to Fig. 4 a , in which a greater number of pH classes cannot be observed visually.

Source: The Authors.

Figure 4 Grid sampling (a) and on-the-go soil sensor (b) of the area belonging to the Juquiá Farm (A2). 

In Figs. 3b, 4b are shown points of higher pH values that were not observed in grid sampling (Figs. 3 a, 4 a), but it was not possible to observe points of lower pH values through maps of sampling by sensor, which were detected in the sampling maps by soil auger (Fig. 4 a ). A possible explanation for this difference in the readings using the sensor is that the sample was collected possibly more accurately in the pre-established layer, which in our case was 15 cm depth. Thus, the long-time under no-till system in the area (about 15 years) contributed to the existence of a layer with higher pH.

Outliers were observed in the spatial distribution of pH among the results obtained by on-the-go soil sensor (Fig4 a ) and by grid sampling (Fig4 a ), performed at common points for both methods. In Fig4 a , it is verified that pH values are higher than the value informed by laboratory analyses (Fig4 a ), indicating that the obtained results were influenced possibly by some external factors.

The variations in the pH values obtained by sensor and values determined in the laboratory can be explained, according to Chirmman et al. (2011), by the interference of external factors, increasing the reading error of the on-the-go soil sensor in relation to the traditional method of laboratory analysis. These authors emphasize that a best calibration of the equipment is needed, but also make reservations to laboratory errors, in which the mixtures of soils from different layers and locations to compose a single sample can mask the result, besides the number of samplings performed for each determination method.

Corassa et al. [5] analyzed the correlation of soil electrical conductivity by on-the-go soil sensor, and observed that this process of determining real-time soil properties requires specific care, such as calibration for each type of soil and crop. However, when there is a high spatial resolution, i.e., a high number of sampling points, associated with oriented or traditional sampling, the use of sensor can lead to satisfactory results and savings in the application of inputs.

In Figs. 3b, 4 b, it was observed that the maps generated through interpolated data from on-the-go soil sensor are more sensitive to spatial variation of pH than maps generated through pH data determined by grid sampling method and from laboratory analyses (Figs. 3 a, 4 a). Cherubin et al. [3] reported that the reduced grid sample improves the accuracy and the variability characterization of soil chemical properties through thematic maps. Souza et al. [23] studied the intensity of sampling in red latosols in relation to the accuracy of values estimated by kriging and reported that the number of sampling points influenced the prediction of values for points not sampled by kriging, where the error was increased whether the number of samples was lower than 100. These authors also emphasize that a minimum of 100 samples is required for geostatistical analysis, in order to prepare thematic maps to guide agronomic management.

Kerry and Oliver [14] also recommend that practitioners of PA should not use the kriging process with variograms calculated based on a low sampling number and with large distances among points, since they may lead to misinterpretations in the application of agricultural inputs and pesticides at variable rate, masking the benefits from the PA. Kamimura et al. [13] conclude that the kriging method for point interpolation in non-sampled regions is a valid process.

Schirrmann et al. [22] report that prediction errors of 0.1 pH units may result in an under or over application of 400 kg.ha-1 of CaO, according to local liming recommendations, which may result in losses for the producer. Although these authors have been successful in correlating pH soil values in northeastern Germany through on-the-go soil sensor and values determined by laboratory, these authors also report that calibration process of the equipment is extremely important, reducing errors in the application of inputs in the area and production costs.

The use of kriging, after fitting semivariograms, proved to be a very valuable process, allowing a best prediction of pH values by on-the-go soil sensor. These information become very useful for the producer, assisting him in the decision making, delimitating management zones in order to perform processes of input application at variable rate, which would not happen whether the map generated from pH values was obtained by traditional analysis.

4. Conclusions

The semivariograms allowed characterizing the magnitude of the spatial variability of pH sampled by grid and by on-the-go soil sensor.

The test of different methods and models allowed identifying which of them best fitted the variable in each sampling method studied.

The interpolation by kriging allowed generating isocolors maps, providing the observation of spatial variability.

The geostatistical analysis of both different sampling methods of soil pH allowed suggesting that the on-the-go soil sensor provided a best representation of soil pH variability in both studied areas and that the sampling number has limited the visual representation of soil pH variability when it is determined by grid sampling method.


[1] Bernardi, A.C.C., Bettiol, G.M., Grego, C.R., Andrade, R.G., Rabello, L.M. and Inamasu, R.Y., Ferramentas de agricultura de precisão como auxílio ao manejo da fertilidade do solo (online). Cadernos de Ciência & Tecnologia, Brasília, [online]. 32(1/2), pp. 205-221. 2015. [consulted: January 12th of 2018]. Available at: Available at: ]

[2] Bottega, E.L., de Queiroz, D.M., de Carvalho-Pinto, F.D.A. and de Souza, C.M.A., Variabilidade espacial de atributos do solo em sistema de semeadura direta com rotação de culturas no cerrado brasileiro. Revista Ciência Agronômica, [online]. 44(1), pp. 1-9. 2012. [consulted: November 10th of 2018]. Available at: Available at: ]

[3] Cherubin, M.R., Santi, A.L., Eitelwein, M.T., Amado, T.J.C., Simon, D.H. and Damian, J.M., Dimensão da malha amostral para caracterização da variabilidade espacial de fósforo e potássio em Latossolo Vermelho. Pesquisa Agropecuária Brasileira, 50(2), pp.168-177, 2015. DOI: 10.1590/S0100-204X2015000200009. [ Links ]

[4] Corá, J.E., Araújo, A.V., Pereira, G.T. and Beraldo, J.M.G., Variabilidade espacial de atributos do solo para adoção do sistema de agricultura de precisão na cultura de cana-de-açúcar. Revista Brasileira de Ciência do Solo, [online]. 28(6), pp. 1013-1021. 2004. [date of reference: January 16th of 2018]. Available at: Available at: ]

[5] Corassa, G.M., Amado, T.J.C., Tabaldi, F.M., Schwalbert, R.A., Reimche, G.B., Dalla Nora, D., & Horbe, T.A.N., Espacialização em alta resolução de atributos da acidez de Latossolo por meio de sensoriamento em tempo real. Pesquisa Agropecuária Brasileira, 51(9), pp. 1306-1316, 2016. DOI: 10.1590/s0100-204x2016000900030 [ Links ]

[6] Cressie, N.A., Statistics for spatial data. John Willey & Sons. Inc., New York. [online]. [Google Scholar].1993. [ Links ]

[7] Cressie, N. and Hawkins, D.M., Robust estimation of the variogram: I. Journal of the International Association for Mathematical Geology, 12(2), pp.115-125, 1980. DOI: /10.1007/BF01035243 [ Links ]

[8] Faraco, M.A., Uribe-Opazo, M.A., Alves da Silva, E.A., Johann, J.A. and Borssoi, J.A., Seleção de modelos de variabilidade espacial para elaboração de mapas temáticos de atributos físicos do solo e produtividade da soja. Revista Brasileira de Ciência do Solo, 32(2), pp. 436-476, 2008 DOI: 10.1590/S0100-06832008000200001 [ Links ]

[9] Ferraz, G.A. e S. Silva , F.M.da , Alves , M. de C , Bueno., R. de L , and da Costa P.A.N ,., Geostatistical analysis of fruit yield and detachment force in coffee. Precision Agriculture, 13, pp.76-89, 2012. DOI: 10.1007/s11119-011-9223-8 [ Links ]

[10] Ferraz, G.A. e S .Silva , F.M.da , .Oliveira , M.S. de , Silva , F.C. da , and Bueno. R. de L ,. Variabilidade espacial da força de desprendimento de frutos do cafeeiro. Engenharia Agrícola, 34(6), pp. 1210-1223, 2014 DOI: 10.1590/S0100-69162014000600016 [ Links ]

[11] Fontoura, S.M.V., Vieira, R.C.B., Viero, F., Bayer, C., Anghinoni, I. and Moraes, R.P., Fertilidade do solo e seu manejo em plantio direto do Centro-Sul do Paraná. 1nd Guarapuava: Fundação Agrária de Pesquisa Agropecuária, Press, 2015. [ Links ]

[12] Isaaks, E.H. and Srivastava, R.M., Applied Geostatistics. Oxford University Press, New York, 1989. [ Links ]

[13] Kamimura, K.M., Rodrigues dos Santos, G., Silva de Oliveira, M., de Souza-Dias Jr, M. and Gontijo-Guimarães, P.T., Variabilidade espacial de atributos físicos de um latossolo vermelho-amarelo, sob lavoura cafeeira. Revista Brasileira de Ciência do Solo, 37(4), pp 877-888, 2013 DOI: 10.1590/S0100-06832013000400006 [ Links ]

[14] Kerry, R. and Oliver, M.A., Determining nugget: sill ratios of standardized variograms from aerial photographs to krige sparse soil data. Precision Agriculture, 9(1-2), pp.33-56. 2008. DOI: 10.1007/s11119-008-9058-0 [ Links ]

[15] King, J.A., Dampney, P.M.R., Lark, R.M., Wheeler, H.C., Bradley, R.I. and Mayr, T.R., Mapping potential crop management zones within fields: use of yield-map series and patterns of soil physical properties identified by electromagnetic induction sensing. Precision Agriculture, 6(2), pp.167-181, 2005. DOI: 10.1007/s11119-005-1033-4 [ Links ]

[16] McBratney, A.B. and Webster, R., Choosing functions for semi‐variograms of soil properties and fitting them to sampling estimates. European Journal of Soil Science, 37(4), pp. 617-639, 1986. DOI: 10.1111/j.1365-2389.1986.tb00392.x [ Links ]

[17] Moral, F.J., Terrón, J.M. and Da Silva, J.M., Delineation of management zones using mobile measurements of soil apparent electrical conductivity and multivariate geostatistical techniques. Soil and Tillage Research, 106(2), pp. 335-343, 2010. DOI: 10.1016/j.still.2009.12.002 [ Links ]

[18] Mulla, D.J., Twenty five years of remote sensing in precision agriculture: key advances and remaining knowledge gaps. Biosystems Engineering, 114(4), pp. 358-371, 2013. DOI: 10.1016/j.biosystemseng.2012.08.009 [ Links ]

[19] Ribeiro Jr, P.J. and Diggle, P.J. geoR: a package for geostatistical analysis. R news, [Online]. 1(2), pp.14-18. 2001. [consulted: January 12th of 2018] Available at: Available at: ]

[20] Ribeiro, L.S., Oliveira, I.R., Dantas, J.S., Silva, C.V., Silva, G.B. and Azevedo J.R ., Variabilidade espacial de atributos físicos de solo coeso sob sistemas de manejo convencional e de plantio direto. Pesq Agropec Bras. 51(6), pp. 1699-1702, 2016. DOI:10.1590/S0100-204X2016000900071 [ Links ]

[21] Sana, R.S., Anghinoni, I., Brandão, Z.N. and Holzschuh, M.J., Variabilidade espacial de atributos físico-químicos do solo e seus efeitos na produtividade do algodoeiro. Revista Brasileira de Engenharia Agrícola e Ambiental, 18(10), pp. 994-1002, 2014. DOI:10.1590/1807-1929/agriambi.v18n10p994-1002 [ Links ]

[22] Schirrmann, M., Gebbers, R., Kramer, E. and Seidel, J., Soil pH mapping with an on-the-go sensor. Sensors, 11(1), pp.573-598, 2011 [ Links ]

[8] Souza, Z.M. de., Souza, G.S. de., Marques Jr, J. and Pereira, G.T., Número de amostras na análise geoestatística e na krigagem de mapas de atributos do solo. Ciência Rural, 44(2), pp.261-268, 2014. DOI: 10.1590/S010384782014000200011. [ Links ]

[24] Vieira, S.R. and Novais, R.D. Geoestatística em estudos de variabilidade espacial do solo. Tópicos em ciência do solo. Viçosa: Sociedade Brasileira de Ciência do Solo, Brasil, 2000, pp. 1-54. [ Links ]

[25] Wang, Y., Huang, T., Liu, J., Lin, Z., Li, S., Wang, R. and Ge, Y., Soil pH value, organic matter and macronutrients contents prediction using optical diffuse reflectance spectroscopy. Computers and Electronics in Agriculture, 111, pp. 69-77, 2015. DOI: 10.1016/j.compag.2014.11.019 [ Links ]

Warrick, A.W. and Nielsen, D.R., Spatial variability of soil physical properties in the field. Applications of Soil Physics, pp. 319-344, 1980. DOI: 10.1016/b978-0-12-348580-9.50018-3. [ Links ]

Webster, R. and Oliver, M.A., Geostatistics for Environmental Scientists, 2nd, John Wiley and Sons, United Kingdom, 2007. [ Links ]

How to cite: Ferraz, G.A.S., Barbosa, B.D.S., Reynaldo, É.F., Santos, S.A., Gonçalves, J.R.M.R. and Ferraz, P.F.P., Spatial variability of soil pH sampled by two methodologies used in precision agriculture in farms under crop rotation. DYNA, 86(209), pp. 289-297, April - June, 2019.

G.A.S. Ferraz, received the BSc. in Agricultural Engineering in 2008, the MSc. in Agricultural Engineering in 2010, and the PhD in Agricultural Engineering in 2012, all of them from the Universidade Federal de Lavras, Lavras, Brazil. From 2011 to 2014, he worked for Universidade Federal Rural do Rio de Janeiro, Brazil and since 2014 for the Universidade Federal de Lavras where he is a full professor in the Agricultural Engineering Department. His research interests include: precision agriculture, geostatistics, agricultural power machinery and agricultural machines. ORCID: 0000-0001-6403-2210

B.D.S. Barbosa, received the BSc. in Agricultural and Enviroment Engineering in 2014 from the Universidade Federal of Minas Gerais, Brazil, MSc. in Water Resources in Agriculte Systems in 2016 from the University Federal of Lavras, Brazil, and the PhD student in Agricultural Engineering since 2016, in University Federal of Lavras, Brazil. His research interests include: precision agriculture, irrigation, engineering irrigation, water resources, energy in irrigated agriculture systems, unmanned aerial vehicles in agriculture. ORCID: 0000-0001-6791-2504

É.F. Reynaldo, received the BsSc in Agronomic Engineer in 2006 from the Universidade Federal do Tocantins, Brazil, the MSc. in Agricultural Machinery at the Universidade de São Paulo - ESALQ - USP (2009) and PhD in Energy in Agriculture at Universidade Estadual Paulista, Brazil. He have work in the areas of agricultural mechanization, application technology for phytosanitary products, precision agriculture and planning and management of mechanized systems focused on corn, soybean, wheat, barley, cotton and sugarcane crops. Currently he is a specialist in equipment at Monsanto in Brazil, acting directly in the elaboration of an Agricultural Maintenance Plan, seeking to maximize the operations of all machines and implements, ensuring the best operational efficiency of the available structure. It manages the information coming from the field in parallel with the coming of the workshops, for the composition of the maintenance planning in the predictive, preventive and corrective plans, proposing the best maintenance alternatives to be implemented in the network. ORCID: 0000-0002-2184-7177

S.A. Santos, received the BSc. In Agricultural Engineering from Universidade Federal de Lavras, Brazil. Currently she is a Master degree student in Agricultural Engineering at Universidade Federal de Lavras, Lavras, Brazil. ORCID: 0000-0003-2262-2663

J.R.M.R. Gonçalves, received the BSc. in Environmental and Sanitary Egineering from Centro Universitário Celso Lisboa in 2011, and recived BSc. in Technology in Environmental Management by the Faculty of Technological Education of the state of Rio de Janeiro, Brazil in 2006. He received a MSc. degree in Agricultural and Environmental Engineering in Universidade Federal Rural do Rio de Janeiro in 2017. Currently he is professor at Centro Universitário Augusto Motta. ORCID: 0000-0003-3321-8818

P.F.P. Ferraz, received the BSc. in Animal Science in 2009, the MSc. in Agricultural Engineering in 2011, and the PhD in Agricultural Engineering in 2013, all of them from the Universidade Federal de Lavras, Lavras, Brazil. Currently she is professor at the Universidade Federal de Lavras. Her research interests include: spatial statistics, mathematical modeling, building environment, bioclimatology, animal buildings and thermal confort. ORCID: 0000-0002-9708-0259

Recebido: 08 de Março de 2018; Revisado: 28 de Fevereiro de 2019; Aceito: 05 de Março de 2019

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons