A geostatistical analysis of the spatial variation of soil mineral nitrogen and potentially available nitrogen within an arable field.

The technology for site-specific applications of nitrogen (N) fertilizer has exposed a gap in our knowledge about the spatial variation of soil mineral N, and that which will become available during the growing season within arable fields. Spring mineral N and potentially available N were measured in an arable field together with gravimetric water content, loss on ignition, crop yield, percentages of sand, silt, and clay, and elevation to describe their spatial variation geostatistically. The areas with a larger clay content had larger values of mineral N, potentially available N, loss on ignition and gravimetric water content, and the converse was true for the areas with more sandy soil. The results suggest that the spatial relations between mineral N and loss on ignition, gravimetric water content, soil texture, elevation and crop yield, and between potentially available N and loss on ignition and silt content could be used to indicate their spatial patterns. Variable-rate nitrogen fertilizer application would be feasible in this field because of the spatial structure and the magnitude of variation of mineral N and potentially available N.

A geostatistical analysis of soil, vegetation, and image data characterizing land surface variation

The elucidation of spatial variation in the landscape can indicate potential wildlife habitats or breeding sites for vectors, such as ticks or mosquitoes, which cause a range of diseases. Information from remotely sensed data could aid the delineation of vegetation distribution on the ground in areas where local knowledge is limited. The data from digital images are often difficult to interpret because of pixel-to-pixel variation, that is, noise, and complex variation at more than one spatial scale. Landsat Thematic Mapper Plus (ETM+) and Satellite Pour l'Observation de La Terre (SPOT) image data were analyzed for an area close to Douna in Mali, West Africa. The variograms of the normalized difference vegetation index (NDVI) from both types of image data were nested. The parameters of the nested variogram function from the Landsat ETM+ data were used to design the sampling for a ground survey of soil and vegetation data. Variograms of the soil and vegetation data showed that their variation was anisotropic and their scales of variation were similar to those of NDVI from the SPOT data. The short- and long-range components of variation in the SPOT data were filtered out separately by factorial kriging. The map of the short-range component appears to represent the patterns of vegetation and associated shallow slopes and drainage channels of the tiger bush system. The map of the long-range component also appeared to relate to broader patterns in the tiger bush and to gentle undulations in the topography. The results suggest that the types of image data analyzed in this study could be used to identify areas with more moisture in semiarid regions that could support wildlife and also be potential vector breeding sites.

Geostatistical analysis of rainfall

Rainfall can be modeled as a spatially correlated random field superimposed on a background mean value; therefore, geostatistical methods are appropriate for the analysis of rain gauge data. Nevertheless, there are certain typical features of these data that must be taken into account to produce useful results, including the generally non-Gaussian mixed distribution, the inhomogeneity and low density of observations, and the temporal and spatial variability of spatial correlation patterns. Many studies show that rigorous geostatistical analysis performs better than other available interpolation techniques for rain gauge data. Important elements are the use of climatological variograms and the appropriate treatment of rainy and nonrainy areas. Benefits of geostatistical analysis for rainfall include ease of estimating areal averages, estimation of uncertainties, and the possibility of using secondary information (e.g., topography). Geostatistical analysis also facilitates the generation of ensembles of rainfall fields that are consistent with a given set of observations, allowing for a more realistic exploration of errors and their propagation in downstream models, such as those used for agricultural or hydrological forecasting. This article provides a review of geostatistical methods used for kriging, exemplified where appropriate by daily rain gauge data from Ethiopia.

Geostatistical radar-raingauge combination with nonparametric correlograms: methodological considerations and application in Switzerland

Modelling spatial covariance is an essential part of all geostatistical methods. Traditionally, parametric semivariogram models are fit from available data. More recently, it has been suggested to use nonparametric correlograms obtained from spatially complete data fields. Here, both estimation techniques are compared. Nonparametric correlograms are shown to have a substantial negative bias. Nonetheless, when combined with the sample variance of the spatial field under consideration, they yield an estimate of the semivariogram that is unbiased for small lag distances. This justifies the use of this estimation technique in geostatistical applications. Various formulations of geostatistical combination (Kriging) methods are used here for the construction of hourly precipitation grids for Switzerland based on data from a sparse realtime network of raingauges and from a spatially complete radar composite. Two variants of Ordinary Kriging (OK) are used to interpolate the sparse gauge observations. In both OK variants, the radar data are only used to determine the semivariogram model. One variant relies on a traditional parametric semivariogram estimate, whereas the other variant uses the nonparametric correlogram. The variants are tested for three cases and the impact of the semivariogram model on the Kriging prediction is illustrated. For the three test cases, the method using nonparametric correlograms performs equally well or better than the traditional method, and at the same time offers great practical advantages. Furthermore, two variants of Kriging with external drift (KED) are tested, both of which use the radar data to estimate nonparametric correlograms, and as the external drift variable. The first KED variant has been used previously for geostatistical radar-raingauge merging in Catalonia (Spain). The second variant is newly proposed here and is an extension of the first. Both variants are evaluated for the three test cases as well as an extended evaluation period. It is found that both methods yield merged fields of better quality than the original radar field or fields obtained by OK of gauge data. The newly suggested KED formulation is shown to be beneficial, in particular in mountainous regions where the quality of the Swiss radar composite is comparatively low. An analysis of the Kriging variances shows that none of the methods tested here provides a satisfactory uncertainty estimate. A suitable variable transformation is expected to improve this.

Multiscale geostatistical analysis of AVHRR, SPOT-VGT, and MODIS global NDVI products

Global NDVI data are routinely derived from the AVHRR, SPOT-VGT, and MODIS/Terra earth observation records for a range of applications from terrestrial vegetation monitoring to climate change modeling. This has led to a substantial interest in the harmonization of multisensor records. Most evaluations of the internal consistency and continuity of global multisensor NDVI products have focused on time-series harmonization in the spectral domain, often neglecting the spatial domain. We fill this void by applying variogram modeling (a) to evaluate the differences in spatial variability between 8-km AVHRR, 1-km SPOT-VGT, and 1-km, 500-m, and 250-m MODIS NDVI products over eight EOS (Earth Observing System) validation sites, and (b) to characterize the decay of spatial variability as a function of pixel size (i.e. data regularization) for spatially aggregated Landsat ETM+ NDVI products and a real multisensor dataset. First, we demonstrate that the conjunctive analysis of two variogram properties – the sill and the mean length scale metric – provides a robust assessment of the differences in spatial variability between multiscale NDVI products that are due to spatial (nominal pixel size, point spread function, and view angle) and non-spatial (sensor calibration, cloud clearing, atmospheric corrections, and length of multi-day compositing period) factors. Next, we show that as the nominal pixel size increases, the decay of spatial information content follows a logarithmic relationship with stronger fit value for the spatially aggregated NDVI products (R2 = 0.9321) than for the native-resolution AVHRR, SPOT-VGT, and MODIS NDVI products (R2 = 0.5064). This relationship serves as a reference for evaluation of the differences in spatial variability and length scales in multiscale datasets at native or aggregated spatial resolutions. The outcomes of this study suggest that multisensor NDVI records cannot be integrated into a long-term data record without proper consideration of all factors affecting their spatial consistency. Hence, we propose an approach for selecting the spatial resolution, at which differences in spatial variability between NDVI products from multiple sensors are minimized. This approach provides practical guidance for the harmonization of long-term multisensor datasets.

Towards more precise sugar beet management based on geostatistical analysis of spatial variability within fields

Within-field variation in sugar beet yield and quality was investigated in three commercial sugar beet fields in the east of England to identify the main associated variables and to examine the possibility of predicting yield early in the season with a view to spatially variable management of sugar beet crops. Irregular grid sampling with some purposively-located nested samples was applied. It revealed the spatial variability in each sugar beet field efficiently. In geostatistical analyses, most variograms were isotropic with moderate to strong spatial dependency indicating a significant spatial variation in sugar beet yield and associated growth and environmental variables in all directions within each field. The Kriged maps showed spatial patterns of yield variability within each field and visual association with the maps of other variables. This was confirmed by redundancy analyses and Pearson correlation coefficients. The main variables associated with yield variability were soil type, organic matter, soil moisture, weed density and canopy temperature. Kriged maps of final yield variability were strongly related to that in crop canopy cover, LAI and intercepted solar radiation early in the growing season, and the yield maps of previous crops. Therefore, yield maps of previous crops together with early assessment of sugar beet growth may make an early prediction of within-field variability in sugar beet yield possible. The Broom’s Barn sugar beet model failed to account for the spatial variability in sugar yield, but the simulation was greatly improved when corrected for early canopy development cover and when the simulated yield was adjusted for weeds and plant population. Further research to optimize inputs to maximise sugar yield should target the irrigation and fertilizing of areas within fields with low canopy cover early in the season.

Geostatistical analysis of the spatial distribution of Rotylenchulus reniformis on cotton cultivated under crop rotation

The spatial distribution of Rotylenchulus reniformis on cotton cultivated in crop rotation with sorghum-peanut-velvetbean was studied using geostatistics. The experimental field, which had been continuously cropped with cotton for 20 years, comprised two 32 x 48 m-grids, each divided in sixty-four 4 x 6 in sampling plots. For all crops, 300 cm(3) soil samples were taken at the center of each plot at crop germination (Pi) and again at harvest (Pf), from which the numbers of nematodes were determined. The results revealed that the spatial distribution of R. reniformis was highly aggregated and with the aid of geostatistical techniques the nematode intensities were mapped and the risk areas accurately identified. Consequently, geostatistics is here considered a useful tool for planning nematode control strategies, particularly in precision agriculture.

Using geostatistical analysis to evaluate the presence of Rotylenchulus reniformis in cotton crops in Brazil: Economic implications

In recent years, the productivity of cotton in Brazil has been progressively decreasing, often the result of the reniform nematode Rotylenchulus reniformis. This species call reduce crop productivity by up to 40%. Nematodes can be controlled by nematicides but, because of expense and toxicity, application of nematicides to large crop areas may be undesirable. In this Work. a methodology using geostatistics for quantifying the risk of nematicide application to small crop areas is proposed. This risk, in economic terms, can be compared to nematicide cost to develop an optimal strategy for Precision Farming, Soil (300 cm(3)) was sampled in a regular network from a R. reniformis-infested area that was a cotton monoculture for 20 years. The number of nematodes in each sample was counted. The nematode number per volume of soil was characterized using geostatistics, and 100 conditional simulations were conducted. Based on the simulations, risk maps were plotted showing the areas where nematicide should be applied in a Precision Farming context. The methodology developed can be applied to farming in countries that ale highly dependent on agriculture, with useful economic implications.

Geostatistical characterization of the spatial distribution of Xylella fastidiosa sharpshooter vectors on citrus

A distribuição espacial das espécies de cigarrinhas (Dilobopterus costalimai Young, Acrogonia sp. e Oncometopia facialis Signoret), vetoras da Xylella fastidiosa, agente causal da Clorose Variegada dos Citros, foi estudada com o uso da geoestatística. As avaliações foram feitas em um pomar comercial de laranja 'Pêra' (Citrus sinensis [L.] Osb.), objetivando estabelecer meios para melhor controle dos vetores e da doença. O monitoramento da ocorrência das cigarrinhas no pomar foi feito através de amostragens mensais, utilizando-se armadilhas adesivas amarelas de 3 x 5, distribuídas uniformemente em 50 pontos na área, dispostas em laranjeiras à altura de 1,5 m do solo e substituídas mensalmente. Acrogonia sp. foi a espécie prevalente nas amostragens. Os resultados possibilitaram ajustar modelos aos semivariogramas da distribuição espacial das três espécies no pomar estudado. Durante os três anos consecutivos de amostragem, as populações de Acrogonia sp., D. costalimai e O. facialis apresentaram modelos de distribuição agregada somente nos meses de verão, inverno e primavera, respectivamente, mostrando a necessidade de monitoramento constante desses vetores para reduzir a sua população em épocas favoráveis ao seu desenvolvimento. Através de parâmetros geoestatísticos foi possível calcular a área de agregação das cigarrinhas no pomar. A espécie Acrogonia sp. apresentou área média de agregação de 15.760 m², enquanto para O. facialis e D. costalimai foi possível constatar áreas médias de agregação de 11.555 m² e 10.980 m², respectivamente. Esses resultados indicaram que para um levantamento seguro de cigarrinhas é necessário pelo menos dispor de uma armadilha por hectare.

Fuzzy Logic and Geostatistical Techniques for Spatialization of Soil Texture in Region with Rough Terrains

From the geotechnical standpoint, it is interesting to analyse the soil texture in regions with rough terrain due to its relation with the infiltration and runoff processes and, consequently, the effect on erosion processes. The purpose of this paper is to present a methodology that provides the soil texture spatialization by using Fuzzy logic and Geostatistic. The results were correlated with maps drawn specifically for the study area. The knowledge of the spatialization of soil properties, such as the texture, can be an important tool for land use planning in order to reduce the potential soil losses during rain seasons. (c) 2011 Published by Elsevier Ltd. Selection and peer-review under responsibility of Spatial Statistics 2011

Comparison between univariate and bivariate geostatistical models for estimating catch per unit of effort (cpue): A simulation study

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Diagnosis of a compacted area for mining through the mechanical resistance to the penetration using geostatistical methods in the Amazon Forest

The mining process promotes land modification and complete landscape alteration. Those alterations in the surface are shown more obviously in the aesthetical aspect as the visual elements of form, texture, climbs, complexity and color which composes the landscape. As a consequence, mining has impacts on the topography, in the soil, in the vegetation and in the area's drainage, with a direct influence on the enterprise. A quite common problem in the recovery of degraded areas in mineral exploration is the compaction of the soil due to the intense traffic of machines and earth movement. The most common problem of the compaction of a degraded surface is an increase of the mechanical resistance to the penetration of plant roots, a reduction of the aeration, an alteration of the flow of water and heat, also in the availability of water and nutrients. Thus, the present work had the basic objective of diagnosing the compaction of an area degraded by mining in a spacial way, through the mechanical resistance and the penetration, to guide the future subsoiling in the area requiring recovery. Through the studies, it was concluded that the krigagem method in agreement with the space variation allows the division of the area under study into sub areas facilitating a future work to reduce costs and unnecessary interference to the atmosphere. The method was shown to be quite appropriate and it can be used in the diagnosis of compaction in a degraded area by mining, foreseeing the subsoiling requirement.

Correlating soil porosity and respective geological unit in Paraíba do sul Valley, Brazil: a geostatistical methodology proposal

This manuscript aims proposing a methodology for correlating soil porosity to the respective geological units using geostatistical analysis techniques, including interpolation data by kriging. The site studied was in Lorena municipality, Paraíba do Sul Valley, southeastern Brazil. Specifically all studies were carried out within an area of 12 km2 located at Santa Edwirges farm. The database comprehended 41 soil samples taken at different geological and geomorphologic units at three different depths: surface, 50 cm and 100 cm depth. The geostatistical analyses results were correlated to a geological mapping specifically elaborated for the site. This mapping accounts for two different geological formations and a geological contact characterized by a shearing zone. The results indicate the existence of a significant relationship between the soil porosity and the respective geological units. The studies revealed that the residual soils from weathered granitic rocks tend to have higher porosities than the residual soils from weathered biotite gneiss rocks, while the soil porosity within the shearing zone is relatively un-sensitive to the respective geological formation. The spatial patterns observed were efficient to evaluate the relationship between the soil porosity, geology unit and the and geomorphology showing a good potential for correlating with others soil properties such as hydraulic conductivity, soil water retention curves and erosion potentials.

COMPARISON MEASURES OF MAPS GENERATED BY GEOSTATISTICAL METHODS

This study uses several measures derived from the error matrix for comparing two thematic maps generated with the same sample set. The reference map was generated with all the sample elements and the map set as the model was generated without the two points detected as influential by the analysis of local influence diagnostics. The data analyzed refer to the wheat productivity in an agricultural area of 13.55 ha considering a sampling grid of 50 x 50 m comprising 50 georeferenced sample elements. The comparison measures derived from the error matrix indicated that despite some similarity on the maps, they are different. The difference between the estimated production by the reference map and the actual production was of 350 kilograms. The same difference calculated with the mode map was of 50 kilograms, indicating that the study of influential points is of fundamental importance to obtain a more reliable estimative and use of measures obtained from the error matrix is a good option to make comparisons between thematic maps.