Geostatistical analysis of adult Rhyzopertha dominica (F.) (Coleoptera: Bostrichidae) in wheat stored at constant temperatures

Insect monitoring and sampling programmes are used in the stored grains industry for the detection and estimation of insect pests. At the low pest densities dictated by economic and commercial requirements, the accuracy of both detection and abundance estimates can be influenced by variations in the spatial structure of pest populations over short distances. Geostatistical analysis of Rhyzopertha dominica populations in 2 dimensions showed that, in both the horizontal and vertical directions and at all temperatures examined, insect numbers were positively correlated over short (0-5cm) distances, and negatively correlated over longer (≥10cm) distances. Analysis in 3 dimensions showed a similar pattern, with positive correlations over short distances and negative correlations at longer distances. At 35°C, insects were located significantly further from the grain surface than at 25 and 30°C. Dispersion metrics showed statistically significant aggregation in all cases. This is the first research using small sample units, high sampling intensities, and a range of temperatures, to show spatial structuring of R. dominica populations over short distances. This research will have significant implications for sampling in the stored grains industry.

Nassau grouper (Epinephelus striatus) spawning aggregations: hydroacoustic surveys and geostatistical analysis

With the near extinction of many spawning aggregations of large grouper and snapper throughout the Caribbean, Gulf of Mexico, and tropical Atlantic, we need to provide baselines for their conservation. Thus, there is a critical need to develop techniques for rapidly assessing the remaining known (and unknown) aggregations. To this end we used mobile hydroacoustic surveys to estimate the density, spatial extent, and total abundance of a Nassau grouper spawning aggregation at Little Cayman Island, Cayman Islands, BWI. Hydroacoustic estimates of abundance, density, and spatial extent were similar on two sampling occasions. The location and approximate spatial extent of the Nassau grouper spawning aggregation near the shelf-break was corroborated by diver visual observations. Hydroacoustic density estimates were, overall, three-times higher than the average density observed by divers; however, we note that in some instances diver-estimated densities in localized areas were similar to hydroacoustic density estimates. The resolution of the hydroacoustic transects and geostatistical interpolation may have resulted in over-estimates in fish abundance, but still provided reasonable estimates of total spatial extent of the aggregation. Limitations in bottom time for scuba and visibility resulted in poor coverage of the entire Nassau grouper aggregation and low estimates of abundance when compared to hydroacoustic estimates. Although the majority of fish in the aggregation were well off bottom, fish that were sometimes in close proximity to the seafloor were not detected by the hydroacoustic survey. We conclude that diver observations of fish spawning aggregations are critical to interpretations of hydroacoustic surveys, and that hydroacoustic surveys provide a more accurate estimate of overall fish abundance and spatial extent than diver observations. Thus, hydroacoustics is an emerging technology that, when coupled with diver observations, provides a comprehensive survey method for monitoring spawning aggregations of fish.

Geostatistical analysis of chlorophyll a in freshwater ecosystems

Horizontal spatial patterns of chlorophyll a in Meiziya Reservoir, Hubei Province, China were analyzed once each month during May, June and July 1997. Two geostatistical techniques, semivariance and fractal analysis, were used to determine variation in chlorophyll a over the whole study area (isotropic) and in different directions (anisotropic). Both techniques provided useful information for detecting and assessing spatial pattern changes of chlorophyll a in freshwater environments. Based on our case study, the distribution of chlorophyll a shifted from aggregated to random distribution in the case of small rainfall event, and then returned to the aggregated distribution after a large rainfall event. On the other hand, the distribution of chlorophyll a became more heterogeneous or random in the direction of water flow (S-N direction) when rainfall events occurred, which was enhanced by rainfall intensity. In contrast, the influence of water flow on the spatial patterns was weak in the E-W direction, and thus the distribution of chlorophyll a remained aggregate with a moderate spatial heterogeneity.

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.

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.