Mapping soil properties, using geostatistical methods in support of precision agriculture and related activities, requires a large number of samples. To reduce soil sampling and measurement time and cost, a combination of field proximal soil sensors was used to predict and map laboratory-measured soil properties in a 3.4-ha pasture field in southeastern Brazil. Sensor soil properties were measured in situ on a 10 × 10-m dense grid (377 samples) using apparent electrical conductivity meters, apparent magnetic susceptibility meter, gamma-ray spectrometer, water content reflectometer, cone penetrometer, and portable X-ray fluorescence spectrometer (pXRF). Soil samples were collected on a 20 × 20-m thin grid (105 samples) and analyzed in the laboratory for organic C, sum of bases, cation exchange capacity, clay content, soil volumetric moisture, and bulk density. Another 25 samples collected throughout the area were also analyzed for the same soil properties and used for independent validation of models and maps. To test whether the combination of sensors enhances soil property predictions, stepwise multiple linear regression (MLR) models of the laboratory soil properties were derived using individual sensor covariate data versus combined sensor data—except for the pXRF data, which were evaluated separately. Then, to test whether a denser grid sample boosted by sensor-based soil property predictions enhances soil property maps, ordinary kriging of the laboratory-measured soil properties from the thin grid was compared to ordinary kriging of the sensor-based predictions from the dense grid, and ordinary cokriging of the laboratory properties aided by sensor covariate data. The combination of multiple soil sensors improved the MLR predictions for all soil properties relative to single sensors. The pXRF data produced the best MLR predictions for organic C content, clay content, and bulk density, standing out as the best single sensor for soil property prediction, whereas the other sensors combined outperformed the pXRF sensor for the sum of bases, cation exchange capacity, and soil volumetric moisture, based on independent validation. Ordinary kriging of sensor-based predictions outperformed the other interpolation approaches for all soil properties, except organic C content, based on validation results. Thus, combining soil sensors, and using sensor-based soil property predictions to increase the sample size and spatial coverage, leads to more detailed and accurate soil property maps.

The scientific field of precision agriculture employs increasingly innovative techniques to optimize inputs, maximize profitability, and reduce environmental impact. However, obtaining a high number of soil samples is challenging in order to make precision agriculture viable. There is a trade-off between the amount of data needed and the time and resources spent to obtain these data compared to the accuracy of the maps produced with more or fewer points. In the present study, the research was based on an exhaustive dataset of apparent electrical conductivity (aEC) containing 3906 points distributed along 26 transects with spacing between each of up to 40 m, measured by the proximal soil sensor EM38-MK2, for a grain-producing area of 72 ha in São Paulo, Brazil. A second sparse dataset was simulated, showing only four transects with a 400 m distance and, in the end, only 162 aEC points. The aEC map via ordinary kriging (OK) from the grid with 26 transects was considered the reference, and two other mapping approaches were used to map aEC via sparse grid: kriging with external drift (KED) and geographically weighted regression (GWR). These last two methods allow the increment of auxiliary variables, such as those obtained by remote sensors that present spatial resolution compatible with the pivot scale, such as data from the Landsat-8, Aster, and Sentinel-2 satellites, as well as ten terrain covariates derived from the Alos Palsar digital elevation model. The KED method, when used with the sparse dataset, showed a relatively good fit to the aEC data (R² = 0.78), with moderate prediction accuracy (MAE = 1.26, RMSE = 1.62) and reasonable predictability (RPD = 1.76), outperforming the GWR method, which had the weakest performance (R² = 0.57, MAE = 1.78, RMSE = 2.30, RPD = 0.81). The reference aEC map using the exhaustive dataset and OK showed the highest accuracy with an R² of 0.97, no systematic bias (ME = 0), and excellent precision (RMSE = 0.56, RPD = 5.86). Management zones (MZs) derived from these maps were validated using soil texture data from clay samples measured at 0–10 cm depth in a grid of 72 points. The KED method demonstrated the highest potential for accurately defining MZs for irrigation, producing a map that closely resembled the reference MZ map, thereby providing reliable guidance for irrigation management.

Finding an ideal sampling design is a crucial stage in detailed soil mapping to assure reasonable accuracy of resulting soil property maps. This study aimed to evaluate the influence of sampling designs and sample sizes on the quality of soil apparent electrical conductivity (ECa) maps from an electromagnetic sensor survey. Twenty-six (26) parallel transects were gathered in a 72-ha plot in Southeastern Brazil. Soil ECa measurements using an on-the-go electromagnetic induction sensor were taken every second using sensor vertical orientation. Two approaches were used to reduce the sample size and simulate kriging interpolations of soil ECa. Firstly, the number of transect lines was reduced by increasing the distance between them; thus, 26 transects with 40 m spacing; 13 with 80 m; 7 with 150 m; and 4 with 300 m. Secondly, random point selection and Douglas-Peucker algorithms were used to derive four reduced datasets by removing 25, 50, 75, and 95% of the points from the ECa survey dataset. Soil ECa was interpolated at 5 m output spatial resolution using ordinary kriging and the four datasets from each simulation (a total of twelve datasets). Map uncertainty was assessed by root mean square error and mean error metrics from 400 random samples previously selected for external map validation. Maps were evaluated on their uncertainty and spatial structure of variation. The transect elimination approach showed that maps produced with transect spacing up to 150 m could preserve the spatial structure of ECa variations. Douglas-Peucker results showed lower nugget values than random point simulations for all selected sample densities, except for a 95% point reduction. The soil ECa maps derived from the 75% reduced dataset (by random sampling or Douglas-Peucker) or from 13 transect lines (80 m spacing) showed reasonable accuracy (RMSE of validation circa 0.7) relative to the map interpolated from all survey points (RMSE of 0.5), suggesting that transect spacing of 80 m and reading intervals greater than one second can be used for improving the efficiency of on-the-go soil ECa surveys.

The precision agriculture scientific field employs increasingly innovative techniques to optimize inputs, maximize profitability, and reduce environmental impacts. Therefore, obtaining a high number of soil samples to make precision agriculture feasible is challenging. This data bottleneck has been overcome by identifying sub-regions based on data obtained through proximal soil sensing equipment. These data can be combined with freely available remote sensing data to create more accurate maps of soil properties. Furthermore, these maps can be optimally aggregated and interpreted for soil heterogeneity through management zones. Thus, this work aimed to create and combine soil management zones from proximal soil sensing and remote sensing data. To this end, data on electrical conductivity and magnetic susceptibility, both apparent, were measured using the EM38-MK2 proximal soil sensor and the contents of the thorium and uranium elements, both equivalent, via the Medusa MS1200 proximal soil sensor for a 72-ha grain-producing area in São Paulo, Brazil. The proximal soil sensing attributes were mapped using ordinary kriging (OK). Maps were also made using kriging with external drift (KED), and the proximal soil sensor attributes data, combined with remote sensing data, such as Landsat-8, Aster, and Sentinel-2 images, in addition to 10 terrain covariables derived from the digital elevation model Alos Palsar. As a result, three management zone maps were produced via the k-means clustering algorithm: using data from proximal sensors (OK), proximal sensors combined with remote sensors (KED), and remote sensors. Seventy-two samples (0–10 cm in depth) were collected and analyzed in a laboratory (1 sample per hectare) for concentrations of clay, calcium, organic carbon, and magnesium to assess the capacity of the management zone maps created using analysis of variance. All zones created using the three data groups could distinguish the different treatment areas. The three data sources used to map management zones produced similar map zones, but the zone map using a combination of proximal and remote data did not show an improvement in defining the management zones, and using only remote sensing data lowered the significance levels of differentiating each zone compared to the OK and KED maps. In summary, this study not only underscores the global applicability of proximal and remote sensing techniques in precision agriculture but also sheds light on the nuances of their integration. The study’s findings affirm the efficacy of these advanced technologies in addressing the challenges posed by soil heterogeneity, paving the way for more nuanced and site-specific agricultural practices worldwide.