Geostatistical modelling of within-field soil and yield variability for management zones delineation: a case study in a durum wheat field

The paper proposes a geostatistical approach for delineating management zones (MZs) based on multivariate geostatistics, showing the use of polygon kriging to compare durum wheat yield among the different MZs (polygons). The study site was a durum wheat field in southern Italy and yield was measured over three crop seasons. The first regionalized factor, calculated with factorial cokriging, was used to partition the field into three iso-frequency classes (MZs). For each MZ, the expected value and standard deviation of yield were estimated with polygon kriging over the three crop seasons. The yield variation was only in part related to soil properties but most of it might be ascribable to different patterns of meteorological conditions. Both components of variation (plant and soil) in a cropping system should then be taken into account for an effective management of rainfed durum wheat in precision agriculture. The proposed approach proved multivariate Geostatistics to be effective for MZ delineation even if further testing is required under different cropping systems and management.     

A simple index to determine if within-field spatial production variation exhibits potential management effects: application in vineyards using yield monitor data

Precision Agriculture - As agricultural data and information becomes more abundant, diagnostics are needed to quickly and efficiently interrogate these data. Indices exist to identify sensor data...     

Leaf area index estimation in vineyards using a ground-based LiDAR scanner

Estimation of grapevine vigour using mobile proximal sensors can provide an indirect method for determining grape yield and quality. Of the various indexes related to the characteristics of grapevine foliage, the leaf area index (LAI) is probably the most widely used in viticulture. To assess the feasibility of using light detection and ranging (LiDAR) sensors for predicting the LAI, several field trials were performed using a tractor-mounted LiDAR system. This system measured the crop in a transverse direction along the rows of vines and geometric and structural parameters were computed. The parameters evaluated were the height of the vines (H), the cross-sectional area (A), the canopy volume (V) and the tree area index (TAI). This last parameter was formulated as the ratio of the crop estimated area per unit ground area, using a local Poisson distribution to approximate the laser beam transmission probability within vines. In order to compare the calculated indexes with the actual values of LAI, the scanned vines were defoliated to obtain LAI values for different row sections. Linear regression analysis showed a good correlation (R ² = 0.81) between canopy volume and the measured values of LAI for 1 m long sections. Nevertheless, the best estimation of the LAI was given by the TAI (R ² = 0.92) for the same length, confirming LiDAR sensors as an interesting option for foliage characterization of grapevines. However, current limitations exist related to the complexity of data process and to the need to accumulate a sufficient number of scans to adequately estimate the LAI.     

Evaluation of the use of LIDAR laser scanner to map pruning wood in vineyards and its potential for management zones delineation

Vine vigour assessment has been a major concern of precision viticulture studies in order to identify areas of uniform vine performance within vineyards. Moreover, the counting and weighing of winter dormant canes is considered as the most informative measurement to indicate vine balance and is commonly performed manually by grape growers for management purposes. The main concern of this measurement is that it is time consuming and laborious and it cannot accommodate detailed sampling density. In the present study, the potential of using laser scanner technology as an automated, easy and rapid way to perform mapping of the winter pruning wood across the vineyard was investigated. The study was conducted during 2010 and 2011, in a one hectare commercial vineyard in central Greece, planted with cv. Agiorgitiko, a traditional Greek variety for the production of red wine. Parameters of topography, soil depth, soil texture, canopy properties (NDVI), yield, and grape quality were mapped and analysed in conjunction to winter canes weighing at pruning time. The mapping of the dormant canes was carried out using a 2D laser scanner sensor prior to pruning and manually measuring the pruning weight on a 10 × 20 m grid. Laser scanner measurements showed significant relationship in both 2010 and 2011 with pruning weight (r = 0.809 and r = 0.829 respectively, p < 0.001), yield and early season NDVI, showing the potential of using laser scanner measurements to assess variability in vine vigour within vineyards. These results suggest that laser scanners offer great promise to characterize within field variability in vine performance.     

Vineyard zone delineation by cluster classification based on annual grape and vine characteristics

This study describes a method for vineyard zone delineation based on spatial interpolation of data on annual monitoring of grape and vine growth from 2007 to 2012 for four commercial vines (Cabernet Sauvignon, Mencía, Merlot and Tempranillo) located in the Bierzo Denomination of Origen (NW Spain). A sampled grid of 20 × 29 m (14 vines/ha) was defined for each vineyard and data were collected for ten soil, six grape composition, three grape production and five vine vigour variables. Continuous maps of each variable were created by spatial interpolation from the sampled points. Several zone delineations were obtained by clustering—using the iterative self-organizing data analysis (ISODATA) algorithm—according to different combinations of the studied variables. The resulting zone delineations were analysed (ANOVA) in order to determine whether the variables in the two cluster classifications for two or three zones were statistically different from each other. The selected delineation was the cluster that included total soluble solids, titratable acidity, total phenolic content, pH, mean cluster weight and length of the internode in two zones. The results point to the feasibility of this approach to vineyard zone delineation. Further research is necessary to confirm the effectiveness of this approach for other locations and evaluate the usefulness of introducing new grape and vine variables.     

Mapping within-field leaf chlorophyll content in agricultural crops for nitrogen management using Landsat-8 imagery

Precision Agriculture - Spatial information on crop nutrient status is central for monitoring vegetation health, plant productivity and managing nutrient optimization programs in agricultural...     

Investigating geostatistical methods to model within-field yield variability of cranberries for potential management zones

Cranberry harvesting methods give only one yield value per field making characterization of within-field variation, the usual first step in precision farming, difficult. Time-consuming berry count yield and fruit rot estimations are the best “ground truth” indication of yield variation within fields. Correlations and coincidence of binary classifications based on less expensive methods such as enhanced vegetation index (EVI) from imagery, and area to point (AtoP) kriging of useable, poor quality and trash yields were compared with this “ground truth”. In general AtoP kriged values gave higher correlations and kappa statistic values with berry counts and fruit rot than EVI. Geostatistical disaggregation of per field yield totals using AtoP kriging with EVI as an external drift (AtoP^KED) was also investigated. Factorial kriging was used to separate the several scales of variation in “ground truth” and EVI data and determine which ones were most spatially coherent/manageable and which related best to the AtoP kriged data. The spatial trend component of pre-harvest berry counts and AtoP kriging of yields both gave a good initial definition of spatially coherent, relatively permanent management zones. They were related to topography and depth of water table in the soil which are key factors governing cranberry yield. AtoP kriging or AtoP^KED are recommended for defining management zones as they are less expensive than berry counts. The value of AtoP kriging to precision farmers for other crops to map soils at the farm scale with some imagery and just one bulked soil sample per field or use nutrient levels associated with each polygon of traditional soil survey maps is discussed in the conclusions.     

Mapping within-field variability in wheat yield and biomass using remote sensing vegetation indices

Precision Agriculture - This paper explored the ability of remote sensing (RS) and meteorological data to map the variability of yield/biomass in cultivated wheat (Triticum aestivum). The...     

Near infrared spectroscopy as a tool for intensive mapping of vineyards soil

The wine industry has always been particularly interested in the influence of the terroir characteristics on the features of a wine. Over the past few years, a growing interest has spurred on the mechanisms by which a particular soil influences the vine’s growth, grape variety characteristics and ultimately wine quality. Near-infrared spectroscopy (NIRS) is a rapid, non-destructive, low-cost and robust analytical method for chemical and physical property determination. Its use for soil characterization, discrimination and compound determination is rapidly increasing. In this work, NIRS data were collected in two vineyards, one in the Dão Wine Region (centre of Portugal) and one in the Vinhos Verdes Wine Region (North of Portugal) previously characterized in terms of soils. Wet, dried and dried-ground soil samples collected from specific vineyard locations were scanned on a Fourier-transform near infrared spectrometer (FTLA 2000, ABB, Quebec, Canada) in diffuse reflectance mode. Spectra were analysed with chemometric tools, namely principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA). Results revealed that dried, ground soil samples gave better results, but not substantially so when compared with wet or dried samples. Discriminant models showed that the NIRS method is able to discriminate the different vineyard soil types, reproducing very accurately the mapping generated by pedology methods. Variations within the same soil type (present at different locations in the vineyard) were also detected by NIRS. The NIRS technology was shown to be suitable for correlating, complementing and perhaps eventually replacing costly, time-consuming vineyard soil mapping methods.     

Bayesian analysis of within-field variability of corn yield using a spatial hierarchical model

Understanding relationships of soil and field topography to crop yield within a field is critical in site-specific management systems. Challenges for efficiently assessing these relationships include spatially correlated yield data and interrelated soil and topographic properties. The objective of this analysis was to apply a spatial Bayesian hierarchical model to examine the effects of soil, topographic and climate variables on corn yield. The model included a mean structure of spatial and temporal co-variates and an explicit random spatial effect. The spatial co-variates included elevation, slope and apparent soil electrical conductivity, temporal co-variates included mean maximum daily temperature, mean daily temperature range and cumulative precipitation in July and August. A conditional auto-regressive (CAR) model was used to model the spatial association in yield. Mapped corn yield data from 1997, 1999, 2001 and 2003 for a 36-ha Missouri claypan soil field were used in the analysis. The model building and computation were performed using a free Bayesian modeling software package, WinBUGS. The relationships of co-variates to corn yield generally agreed with the literature. The CAR model successfully captured the spatial association in yield. Model standard deviation decreased about 50% with spatial effect accounted for. Further, the approach was able to assess the effects of temporal climate co-variates on corn yield with a small number of site-years. The spatial Bayesian model appeared to be a useful tool to gain insights into yield spatial and temporal variability related to soil, topography and growing season weather conditions.

Using RapidEye imagery to identify within-field variability of crop growth and yield in Ontario,Canada

Precision Agriculture - Remote sensing has been recognized as a cost-effective way to detect the spatial and temporal variability of crop growth and productivity. In this study, multispectral...     

Double seismic zone for deep earthquakes in the izu-bonin subduction zone

A double seismic zone for deep earthquakes was found in the Izu-Bonin region. An analysis of SP-converted phases confirms that the deep seismic zone consists of two layers separated by approximately 20 kilometers. Numerical modeling of the thermal structure implies that the hypocenters are located along isotherms of 500 degrees to 550 degrees C, which is consistent with the hypothesis that deep earthquakes result from the phase transition of metastable olivine to a high-pressure phase in the subducting slab.     

A comparison of different algorithms for the delineation of management zones

One approach to the application of site-specific techniques and technologies in precision agriculture is to subdivide a field into a few contiguous homogenous zones, often referred to as management zones (MZs). Delineating MZs can be based on some sort of clustering, however there is no widely accepted method. The application of fuzzy set theory to clustering has enabled researchers to account better for the continuous variation in natural phenomena. Moreover, the methods based on non-parametric density estimation can detect clusters of unequal size and dispersion. The objectives of this paper were to: (1) compare different procedures for creating management zones and (2) determine the relation of the MZs delineated with potential yield. One hundred georeferenced point measurements of soil and crop properties were obtained from a 12 ha field cropped with durum wheat for two seasons. The trial was carried out at the experimental farm of CRA-CER in Foggia (Italy). All variables were interpolated on a 1 × 1 m grid using the geostatistical techniques of kriging and cokriging. The techniques compared to identify MZs were: (1) the ISODATA method, (2) the fuzzy c-means algorithm and (3) a non-parametric density algorithm. The ISODATA method, which was the simplest, subdivided the field into three distinct classes of suitable size for uniform management, whereas the other two methods created two classes. The non-parametric density algorithm characterized the edge properties between adjacent clusters more efficiently than the fuzzy method. The clusters from the non-parametric density algorithm and yield maps for three seasons (2005–2006, 2006–2007 and 2007–2008) were compared and agreement measures were computed. The kappa coefficients for the three seasons were negative or small positive values which indicate only slight agreement. These results illustrate the importance of temporal variation in spatial variation of yield in rainfed conditions, which limits the use of the MZ approach.     

Data fusion techniques for delineation of site-specific management zones in a field in UK

Fusion of different data layers, such as data from soil analysis and proximal soil sensing, is essential to improve assessment of spatial variation in soil and yield. On-line visible and near infrared (Vis–NIR) spectroscopy have been proved to provide high resolution information about spatial variability of key soil properties. Multivariate geostatistics tools were successfully implemented for the delineation of management zones (MZs) for precision application of crop inputs. This research was conducted in a 18 ha field to delineate MZs, using a multi-source data set, which consisted of eight laboratory measured soil variables (pH, available phosphorus (P), cation exchange capacity, total nitrogen (TN), total carbon (TC), exchangeable potassium (K), sand, silt) and four on-line collected Vis–NIR spectra-based predicted soil variables (pH, P, K and moisture content). The latter set of data was predicted using the partial least squares regression (PLSR) technique. The quality of the calibration models was evaluated by cross-validation. Multi-collocated cokriging was applied to the soil and spectral data set to produce thematic spatial maps, whereas multi-collocated factor cokriging was applied to delineate MZ. The Vis–NIR predicted K was chosen as the exhaustive variable, because it was the most correlated with the soil variables. A yield map of barley was interpolated by means of the inverse distance weighting method and was then classified into 3 iso-frequency classes (low, medium and high). To assess the productivity potential of the different zones of the field, spatial association between MZs and yield classes was calculated. Results showed that the prediction performance of PLSR calibration models for pH, P, MC and K were of excellent to moderate quality. The geostatistical model revealed good performance. The estimates of the first regionalised factor produced three MZs of equal size in the studied field. The loading coefficients for TC, pH and TN of the first factor were highest and positive. This means that the first factor can be assumed as a synthetic indicator of soil fertility. The overall spatial association between the yield classes and MZs was about 40 %, which reveals that more than 50 % of the yield variation can be attributed to more dynamic factors than soil parameters, such as agro-meteorological conditions, plant diseases and nutrition stresses. Nevertheless, multivariate geostatistics proved to be an effective approach for site-specific management of agricultural fields.     

A Web-based system for vineyards management, relating inventory data, vectors and images

A Web-based system is presented, integrating spatial information from remote sensing images, GPS measurements and inventory data. Monitoring, research and management of the grape production at Rio Grande do Sul State, Brazil (624,000 metric tons in 2006) can be done through a system entirely based in open-source codes. Information from three different sources are integrated: Information on production comes from the State’s Viticultural Inventory, detailing the regional yearly grape production; positional information comes from field GPS measurements of vine parcels and the system’s capabilities of making maps; imagery information comes from aerial or satellite images. Further capabilities on image classification leads to the identification of vine areas; this allows crossing this data with information provided by owners. This system allows both the monitoring of grape production for administrative purposes and investigations either on regional land and soil cover, or in other applications derived from image classification. The system concept can be readily extended to other applications in land use monitoring.     

A flexible approach to managing variability in grain yield and nitrate leaching at within-field to farm scales

We up-scaled the APSIM simulation model of crop growth, water and nitrogen dynamics to interpret and respond to spatial and temporal variations in soil, season and crop performance and improve yield and decrease nitrate leaching. Grain yields, drainage below the maximum root depth and nitrate leaching are strongly governed by interaction of plant available soil water storage capacity (PAWC), seasonal rainfall and nitrogen supply in the water-limited Mediterranean-type environment of Western Australia (WA). APSIM simulates the interaction of these key system parameters and the robustness of its simulations has been rigorously tested with the results of several field experiments covering a range of soil types and seasonal conditions in WA. We used yield maps, soil and weather data for farms at two locations in WA to determine spatial and temporal patterns of grain yield, drainage below the maximum root depth and nitrate leaching under a range of weather, soil and nitrogen management scenarios. On one farm, we up-scaled APSIM simulations across the whole farm using local weather and fertiliser use data and the average PAWC values of soil type polygons. On a 70 ha field on another farm, we used a linear regression of apparent soil electrical conductivity (EC_a) measured by EM38 against PAWC to transform an EC_a map of the field into a high resolution (5 m grid) PAWC map. We then used regressions of simulated yields, drainage below the maximum root depth and nitrate leaching on PAWC to upscale the APSIM simulations for a range of weather and fertiliser management scenarios. This continuous mapping approach overcame the weakness of the soil polygons approach, which assumed uniformity in soil properties and processes within soil type polygons. It identified areas at greatest financial and environmental risks across the field, which required focused management and simulated their response to management interventions. Splitting nitrogen applications increased simulated wheat yields at all sites across the field and decreased nitrate leaching particularly where the water storage capacity of the soil was small. Low water storage capacity resulted in both low wheat yields and large leaching loss. Another management option to decrease leaching may be to grow perennial vegetation that uses more water and loses less by drainage.Paper from the 5th European Conference on Precision Agriculture (5ECPA), Uppsala, Sweden, 2005