Mapping of soil insect infestations sampled by excavation and acoustic methods.

Geostatistical analysis was used to map traditionally and acoustically sampled populations of soil invertebrates at two different times in a hayfield at Grove Hill, AL, and once in an Auburn, AL, hayfield. The distributions of nearly all the soil invertebrates and their sounds were nonrandom in all three mapping studies. The maps constructed by excavation and acoustic sampling methods were compared by correlating the estimated (kriged) soil invertebrate counts with the estimated (kriged) counts of sounds per minute (pulse rate). Acoustic and traditional estimates were positively correlated in the Auburn study. Kriged estimates for green June beetle grub counts overlapped significantly with kriged estimates of sound pulse rate (R2 = 0.47). Overlap with sound pulse rates increased slightly when other soil organisms were counted along with green June beetle grubs: estimates of sound pulse rates were significantly correlated with counts of all white grubs (R2 = 0.50), all white grubs with earthworms (R2 = 0.52), all white grubs with earthworms and earwigs (R2 = 0.59), and total invertebrates (R2 = 0.59). The correlation between acoustic and traditional estimates was not significant at Grove Hill in either year, possibly because of a lack of experience in signal analysis or because the soil invertebrates may not have generated enough sounds to be detected above the background noise levels. These results suggest that acoustic technology is a promising tool for detecting insect pests in soil, but that further study and additional analysis are needed to improve interpretation of acoustic data obtained in the field.     

Spatio-temporal dynamics of corn rootworm, Diabrotica spp., adults and their spatial association with environment

Multiyear and multilocation studies were conducted to investigate the within‐ and the between‐year spatial dynamics of corn rootworms Diabrotica spp. (Coleoptera: Chrysomelidae), adults and their spatial associations with environmental factors in cornfields. Grid‐based spatial sampling was conducted by trapping adults emerging from the soil (i.e., ‘emergence‐trap counts’) and by counting the number of adults in the ear zone of corn plant (i.e., ‘ear‐zone count’). Spatial analysis with distance indices (SADIE) was used to determine spatial distribution patterns and to investigate spatial associations. Ear‐zone counts and emergence‐trap counts were spatially dynamic within a year and more frequently aggregated in the middle of growing seasons and random early and late in the season. However, during the peak population periods, spatial distribution of ear‐zone counts and emergence‐trap counts were significantly consistent between years, indicating predictability of future spatial distributions. Spatial distribution of emergence‐trap counts and ear‐zone counts were positively associated with soil moisture and presence of corn plants with silks, respectively. This study demonstrated that within‐year spatial distribution patterns were dynamic and that there was between‐year spatial consistency of corn rootworm adult distributions. Such information can be used to improve preventative management of corn rootworms.     

Spatial avoidance of patches of polluted chernozem soils by soil invertebrates

Soil invertebrates and heavy metal concentrations are heterogeneously distributed in the soil of steppe plots surrounding an iron mining enterprise in southern Russia. This study assesses whether patches of high soil invertebrate abundance coincide with patches of low concentrations of pollutants. For this aim, spatial analysis by distance indices (SADIE) was applied. Three valleys in Belogorye Nature Reserve were chosen. One valley faced the tailing pond to the north and the other two faced south-east or south-west. Two sampling plots were chosen in each valley, 60 m apart from each other. On every plot 16 soil cores were collected from a grid of 4 × 4 units with a 5-m distance between each sample unit. Each soil core had an area of 76 cm² and was 12–15 cm deep. All macroinvertebrates were hand-sorted and identified to family. Abundance of soil invertebrates was not controlled by patches of metal concentration in the soil. Epigaeic groups, like insects and other invertebrates inhabiting the litter layer, were not directly associated with local parameters of the soil. On the contrary, belowground invertebrate abundance (elaterid larvae and earthworms) showed significant dissociation with some heavy metal (Fe, Pb, Zn) concentrations in the soil. The patchiness of soil pollution may act as a leading factor of belowground soil invertebrate distribution. The spatial structure of animal populations in industrially transformed soils needs further research.     

Spatial patterns of surface soil properties and vegetation in a Mediterranean semi-arid steppe

In arid and semi-arid areas with sparse vegetation cover, the spatial pattern of surface soil properties affects water and nutrient flows, and is a question of considerable interest for understanding degradation processes and establishing adequate management measures. In this study, we investigate the spatial distribution of vegetation and surface soil properties (biological crusts, physical crusts, mosses, rock fragments, earthworm casts, fine root accumulation and below-ground stones) in a semi-arid Stipa tenacissima L. steppe in SE Spain. We applied the combination of spatial analysis by distance indices (SADIE) and geostatistics to assess the spatial pattern of soil properties and vegetation, and correlation analyses to explore how these patterns were related. SADIE analysis detected significant clumped patterns in the spatial distribution of vegetation, mosses, fine root accumulation and below-ground stone content. Contoured SADIE index of clustering maps suggested the presence of patchiness in the distribution of earthworm casts, fine roots, below-ground stone content, mosses and biological crusts. Correlation analyses suggested that spatial pattern of some soil properties such as biological crusts, moss cover, surface rock fragments, physical crusts and fine roots were significantly related with above-ground plant distribution. We discuss the spatial arrangement of surface soil properties and suggest mechanistic explanations for the observed spatial patterns and relationships.     

Spatial analyses of ecological count data: A density map comparison approach

Analysing spatial patterns of population distributions may help to infer the decisive underlying ecological processes. Here we propose a method adapted to the spatial analysis of count data. Named MAPCOMP (MAP COMParison), it is based on the calculation of a formal distance, the Hellinger distance, between the density map of counts and the density map of sampling effort. Statistical tests of spatial homogeneity are based on count permutations across sampling sites and on valuable properties of the Hellinger distance. We assessed the efficiency of MAPCOMP by simulating different types and locations of clusters of individuals and compared its performance to the classical red–blue SADIE method, used as a reference. The two methods were also compared with respect to counts of codling moth larvae in orchards. Thanks to its better theoretical properties than SADIE, MAPCOMP was efficient in detecting spatial inhomogeneity when clusters were located on square or elongated spatial domains and more or less close to the edges, even for small sample sizes. It also appeared not very sensitive to edge effects. Another advantage of MAPCOMP is a bandwidth parameter that allows assessing the spatial extent of heterogeneity, if any.     

The spatial pattern of soil-dwelling termites in primary and logged forest in Sabah, Malaysia

Abstract.  1. Primary and logged lowland dipterocarp forest sites were sampled for subterranean termites using soil pits located on a grid system in order to detect any patchiness in their distribution.
2. A spatial pattern in termite distributions was observed in the primary and logged sites, but the response differed between soil-feeding and non-soil-feeding termites.
3. Spatial analysis showed that soil-feeding termites were homogeneously distributed in the primary forest but significantly aggregated in the logged forest. This pattern was reversed for non-soil-feeding termites and may result from differences in resource provisioning between the two sites.
4. Gaps in termite distribution comprised a greater area than patches for both feeding groups and sites, but gaps dominated the logged site.
5. A significant association between soil-feeding and non-soil-feeding termite distributions occurred at both sites. This arose from an association between patches in the primary forest and between gaps in the logged forest.
6. Termite spatial pattern was optimally observed at a minimum extent of 64 m and lag of 2 m.
7. The spatially explicit SADIE (Spatial Analysis by Distances IndicEs) analyses were more successful than (non-spatially explicit) multivariate analysis (Canonical Correspondence Analysis) at detecting associations between termite spatial distributions and that of other biotic and abiotic variables.     

Analysis of spatial patterns at a geographical scale over north-western Europe from point-referenced aphid count data

The spatial analysis by distance indices (SADIE) technique was developed to evaluate the spatial pattern of point-referenced count data as well as the spatial association between two sets of data sharing the same point locations. This paper presents an analysis of spatial patterns in aphid count data and the association of these data with climate across north-west Europe. The paper tests the applicability of the technique to large geographical areas. Aggregation and cluster indices were calculated for the total annual abundance of the peach-potato aphid Myzus persicae (Sulzer) and for the annual mean rainfall and temperature at aphid monitoring sites. Association indices demonstrated the stability in time of aphid spatial structures and the correlation between aphid density and climate patterns. Groups of relatively large numbers of aphids, termed patches, and groups of relatively small numbers of aphids, termed gaps, were located and their mean size estimated. The aphid patterns were quite stable in time and the spatial patterns of temperature and rainfall were weakly associated with M. persicae annual abundance. Similarities were observed between the results of SADIE and those from the more widely used technique of spatial autocorrelation (SAC). However, the SADIE association index has the advantage of quantifying the possible associations between aphid data and the factors that determine population distribution. Thus, high temperature and low rainfall were identified as environmental factors that were positively associated with aphid abundance across north-west Europe.     

Spatially explicit analyses of anopheline mosquitoes indoor resting density: implications for malaria control

Background

The question of sampling and spatial aggregation of malaria vectors is central to vector control efforts and estimates of transmission. Spatial patterns of anopheline populations are complex because mosquitoes'' habitats and behaviors are strongly heterogeneous. Analyses of spatially referenced counts provide a powerful approach to delineate complex distribution patterns, and contributions of these methods in the study and control of malaria vectors must be carefully evaluated.

Methodology/Principal Findings

We used correlograms, directional variograms, Local Indicators of Spatial Association (LISA) and the Spatial Analysis by Distance IndicEs (SADIE) to examine spatial patterns of Indoor Resting Densities (IRD) in two dominant malaria vectors sampled with a 5×5 km grid over a 2500 km² area in the forest domain of Cameroon. SADIE analyses revealed that the distribution of Anopheles gambiae was different from regular or random, whereas there was no evidence of spatial pattern in Anopheles funestus (Ia = 1.644, Pa<0.05 and Ia = 1.464, Pa>0.05, respectively). Correlograms and variograms showed significant spatial autocorrelations at small distance lags, and indicated the presence of large clusters of similar values of abundance in An. gambiae while An. funestus was characterized by smaller clusters. The examination of spatial patterns at a finer spatial scale with SADIE and LISA identified several patches of higher than average IRD (hot spots) and clusters of lower than average IRD (cold spots) for the two species. Significant changes occurred in the overall spatial pattern, spatial trends and clusters when IRDs were aggregated at the house level rather than the locality level. All spatial analyses unveiled scale-dependent patterns that could not be identified by traditional aggregation indices.

Conclusions/Significance

Our study illustrates the importance of spatial analyses in unraveling the complex spatial patterns of malaria vectors, and highlights the potential contributions of these methods in malaria control.     

Spatial relationships between two Agriotes click-beetle species and wireworms in agricultural fields

1 The production of new insect pheromones for pest monitoring proceeds at a greater rate than their evaluation, with the consequential possibility of premature introduction. Fundamental to their successful deployment is the determination of a consistent relationship between adult male pheromone trap catches and pest damage. In the present study, adult pheromone traps and larval bait traps were used to examine spatial relationships between two species of Agriotes beetle and wireworms at the field scale. 2 The spatial distributions of adult male Agriotes lineatus and Agriotes obscurus in two fields were determined and compared with the distribution of their larvae. Data were assembled as spatially referenced trap counts, and analysed for evidence of aggregation and clustering using Spatial Analysis by Distance IndicEs (SADIE) methodology. Spatial stabilities of adult populations between sampling dates were tested using association tests. Spatial and quantitative linkages between adult and larval trap catches were also tested. Moreover, a new way of adapting SADIE methodologies is presented for situations where two datasets within an area do not share the same sampling points. 3 There was no significant difference in variance : mean relationships for the two species but there were differences in their spatial distributions, and this is a definitive example of the general argument stating that it is important to consider spatial as well as count data in ecological studies. The spatial distribution of A. lineatus varied between sampling occasions at both sites whereas A. obscurus had consistently significant SADIE indices over time at one site, and adult catches could also be linked to larval distributions and counts. It is proposed that observed differences between the two species can be explained by interference between traps and dissimilar movement rates. There was some evidence of an edge effect at the field boundaries. 4 The distance between pheromone traps is related to the time that elapses before adjacent traps interfere with trap captures and this limits the detection of statistically significant spatial patterns. It is shown that the current practice of adding trap counts for different Agriotes species and treating them as numerically equivalent is insufficiently robust to be recommended at this stage. 5 The implications for the use of sex pheromone traps in wireworm pest management are considered. It is concluded that pheromone traps, as currently used, will not reliably indicate where wireworms occur in a field, and that the complexity of interpreting adult male trap counts limits quantitative predictions of population size.     

Spatial and temporal distribution of Empoasca vitis within a vineyard

1 The intra‐plot spatial distribution of the green leafhopper Empoasca vitis (Goethe) (Cicadellidae, Empoascini) was assessed over three successive growing seasons in a Bordeaux vineyard. Weekly measurements of adult trappings and nymphal counts were performed in a 1.7‐ha plot on 130 sampling points. 2 Statistical and geostatistical analyses revealed inconsistent spatial distributions of adults and nymphs in spring, but consistent spatial distributions in summer, which were identical every year. 3 Similarities between spatial distributions of the mobile adults and the sedentary and aggregated nymphs strongly suggest that adults disperse inside the plot to areas preferred for egg deposition. 4 The similarity of summer populations among years suggests that this insects distribution is based on (perennial) differences in plot characteristics.     

The spatial distribution of larvae of Culicoides impunctatus biting midges

The spatial distribution of Culicoides impunctatus Goetghebuer (Diptera: Ceratopogonidae) larvae was sampled at a site in western Scotland, and geostatistical analyses were used to quantify spatial dependencies. Nested sampling and analysis indicated that most of the spatial variance occurred within distances of 160-640 cm, levelling off at distances >640 cm. Semivariograms for transformed data from three 100 m x 100 m grids showed similar, isotropic patterns for larval counts, soil organic content, soil water content and the distribution of Juncus spp. rushes, with the variance increasing with separation distance. High levels of significance were associated with power models used to describe the semivariograms, which was indicative of the absence of a plateau (or 'sill') in the respective data. Correlation analysis of transformed data revealed significantly positive relationships between larval counts and soil pH, soil percentage organic content, soil percentage water content and also the distribution of Sphagnum spp., Juncus spp. and Myrica gale. There were also significantly negative relationships between larval counts and the distribution of Pteridium aquilinum and all mosses other than Sphagnum spp. The results suggest a far more structured and predictable pattern of C. impunctatus larval sites than previous studies and are discussed in relation to their application in localized Culicoides control and to studies of the mechanisms determining the spatial distribution of C. impunctatus larvae.     

Considerations for the use of SADIE statistics to quangify spatial patterns

The SADIE methodology has recently been developed to quangify spatial patterns, particularly in determining the randomness of observed patterns, estimating patches and gaps, and quangifying correlation between species. We used argificially generated data sets to investigate the relationship of the SADIE pattern and association statistics with patch number and spatial positions of the patches. Results showed that SADIE statistics are critically dependent on the absolute positions of quadrat counts, and not just on their relative positions. In particular, degree of patchiness increased as the patch was displaced away from the centre. Therefore, the dependence of SADIE statistics on absolute sampling positions should be taken into consideration when using and interpreting SADIE-based results.     

Spatio-temporal distribution of Ricania shantungensis (Hemiptera: Ricaniidae) in chestnut fields: Implications for site-specific management

Chestnut is one of most important forest products in Korea. However, Ricania shantungensis, an invasive species, has been increasingly causing economic damage to chestnut. To increase management efficiency of R. shantungensis in chestnut fields, its spatio-temporal distributions were analyzed using spatial analysis by distance indices (SADIE) and semivariogram. Experiments were conducted by observing selected, marked, and geo-referenced trees in two commercial chestnut fields for each developmental stage (i.e., spring egg, nymph, adult, and fall egg) of R. shantungensis from 2017 to 2019. Spatial distributions of R. shantungensis were statistically (P < 0.05) aggregated except for its nymphal stage. Spatial associations of its distributions also showed the statistically (P < 0.05) positive associations regardless of years or developmental stages except for the spatial relationship between egg and nymphal stages. These results indicated that site-specific management for R. shantungensis would be applicable. In application of site-specific management for R. shantungensis, the distance of spatial dependence for nymphs, 60 m, should be considered to minimize reoccurrence possibility and additional samplings. Incorporating site-specific management into pest control program of R. shantungensis could increase its control efficiency in chestnut fields.     

Fish Species and Community Distributions as Proxies for Seafloor Habitat Distributions: The Stellwagen Bank National Marine Sanctuary Example (Northwest Atlantic, Gulf Of Maine)

Defining the habitats of fishes and associated fauna on outer continental shelves is problematic given the paucity of data on the actual types and distributions of seafloor habitats. However many regions have good data on the distributions of fishes from resource surveys or catch statistics because of the economic importance of the fisheries. Fish distribution data (species or communities) have been used as a proxy for the distribution of habitats to develop precautionary conservation strategies for habitat protection (e.g., marine protected areas, fishing gear restrictions). In this study we assessed the relationships between the distributions of fish communities and species derived from trawl survey data with the spatial distribution of sediment types determined by sampling and acoustic reflectance derived from multibeam sonar surveys in Stellwagen Bank National Marine Sanctuary. Fish communities were correlated with reflectance values but all communities did not occur in unique sediment types. This suggests that use of community distributions as proxies for habitats should include the caveat that a greater number of communities within an area could indicate a greater range of habitat types. Single species distributions showed relationships between abundance and reflectance values. Trawl catches with low abundances had wide variations in reflectance values while those with high abundances had narrower ranges indicating habitat affinities. Significant non-random frequency-dependent relationships were observed for 17 of 20 species although only 12 of 20 species had significant relationships based on rank correlation. These results suggest that species distributions based on trawl survey data can be used as proxies for the distribution of seafloor habitats. Species with known habitat associations can be used to infer habitat requirements of co-occurring species and can be used to identify a range of habitat types.     

Spatiotemporal patterns of the fish assemblages downstream of the Gezhouba Dam on the Yangtze River

An explicit demonstration of the changes in fish assemblages is required to reveal the influence of damming on fish species. However, information from which to draw general conclusions regarding changes in fish assemblages is insufficient because of the limitations of available approaches. We used a combination of acoustic surveys, gillnet sampling, and geostatistical simulations to document the spatiotemporal variations in the fish assemblages downstream of the Gezhouba Dam, before and after the third impoundment of Three Gorges Reservoir (TGR). To conduct a hydroacoustic identification of individual species, we matched the size distributions of the fishes captured by gillnet with those of the acoustic surveys. An optimum threshold of target strength of ?50 dB re 1 m² was defined, and acoustic surveys were purposefully extended to the selected fish assemblages (i.e., endemic Coreius species) that was acquired by the size and species selectivity of the gillnet sampling. The relative proportion of fish species in acoustic surveys was allocated based on the composition (%) of the harvest in the gillnet surveys. Geostatistical simulations were likewise used to generate spatial patterns of fish distribution, and to determine the absolute abundance of the selected fish assemblages. We observed both the species composition and the spatial distribution of the selected fish assemblages changed significantly after implementation of new flow regulation in the TGR, wherein an immediate sharp population decline in the Coreius occurred. Our results strongly suggested that the new flow regulation in the TGR impoundment adversely affected downstream fish species, particularly the endemic Coreius species. To determine the factors responsible for the decline, we associated the variation in the fish assemblage patterns with changes in the environment and determined that substrate erosion resulting from trapping practices in the TGR likely played a key role.     

Spatial patterns and interspecific relations analysis help to better understand species distribution patterns in a Mediterranean high mountain grassland

The aim of this study is to investigate distribution patterns of species occurring at definite spatial scales and to address the main ecological factors that structure a Mediterranean high mountain grassland community. Following the protocols of the GLORIA long-term study, four summits were sampled during the year 2006 and 2007 in Sierra de Guadarrama (Spain). The data recorded on two of those summits were analysed for spatial pattern of species, interspecific associations and environmental relationships were determined from the recorded data. This was done by multivariate analyses (ordinations), SADIE analyses (distance indices) and Pearson χ² test, respectively. Results showed that chamaephyte species were more abundant in less disturbed situations whilst caespitose hemicryptophyte species and mosses were more abundant in disturbed ones. High mountain species were more abundant in less disturbed situations. Higher spatial heterogeneity was detected on northern and eastern slopes, which could be related to the increased environmental severity of northern exposures and to the greater intensity of the disturbances on eastern exposures. The frequent aggregated distribution of a few species may be related to facilitation processes due to their frequent participation in positive interactions. Finally we observed that positive pairwise associations of species were more frequent than negative associations. The relevance of the studied relationships lied in a better understanding of those threats on high mountain biodiversity induced by climate warming.     

Influence of land-use intensity on the spatial distribution of N-cycling microorganisms in grassland soils

A geostatistical approach using replicated grassland sites (10 m × 10 m) was applied to investigate the influence of grassland management, i.e. unfertilized pastures and fertilized mown meadows representing low and high land-use intensity (LUI), on soil biogeochemical properties and spatial distributions of ammonia-oxidizing and denitrifying microorganisms in soil. Spatial autocorrelations of the different N-cycling communities ranged between 1.4 and 7.6 m for ammonia oxidizers and from 0.3 m for nosZ-type denitrifiers to scales >14 m for nirK-type denitrifiers. The spatial heterogeneity of ammonia oxidizers and nirS-type denitrifiers increased in high LUI, but decreased for biogeochemical properties, suggesting that biotic and/or abiotic factors other than those measured are driving the spatial distribution of these microorganisms at the plot scale. Furthermore, ammonia oxidizers (amoA ammonia-oxidizing archaea and amoA ammonia-oxidizing bacteria) and nitrate reducers (napA and narG) showed spatial coexistence, whereas niche partitioning was found between nirK- and nirS-type denitrifiers. Together, our results indicate that spatial analysis is a useful tool to characterize the distribution of different functional microbial guilds with respect to soil biogeochemical properties and land-use management. In addition, spatial analyses allowed us to identify distinct distribution ranges indicating the coexistence or niche partitioning of N-cycling communities in grassland soil.     

Characterising spatial patterns of neglected tropical disease transmission using integrated sero-surveillance in Northern Ghana

BackgroundAs prevalence decreases in pre-elimination settings, identifying the spatial distribution of remaining infections to target control measures becomes increasingly challenging. By measuring multiple antibody responses indicative of past exposure to different pathogens, integrated serological surveys enable simultaneous characterisation of residual transmission of multiple pathogens.Methodology/Principal findingsHere, we combine integrated serological surveys with geostatistical modelling and remote sensing-derived environmental data to estimate the spatial distribution of exposure to multiple diseases in children in Northern Ghana. The study utilised the trachoma surveillance survey platform (cross-sectional two-stage cluster-sampled surveys) to collect information on additional identified diseases at different stages of elimination with minimal additional cost. Geostatistical modelling of serological data allowed identification of areas with high probabilities of recent exposure to diseases of interest, including areas previously unknown to control programmes. We additionally demonstrate how serological surveys can be used to identify areas with exposure to multiple diseases and to prioritise areas with high uncertainty for future surveys. Modelled estimates of cluster-level prevalence were strongly correlated with more operationally feasible metrics of antibody responses.Conclusions/SignificanceThis study demonstrates the potential of integrated serological surveillance to characterise spatial distributions of exposure to multiple pathogens in low transmission and elimination settings when the probability of detecting infections is low.