Estimating and controlling for spatial structure in the study of ecological communities

Aim Variation partitioning based on canonical analysis is the most commonly used analysis to investigate community patterns according to environmental and spatial predictors. Ecologists use this method in order to understand the pure contribution of the environment independent of space, and vice versa, as well as to control for inflated type I error in assessing the environmental component under spatial autocorrelation. Our goal is to use numerical simulations to compare how different spatial predictors and model selection procedures perform in assessing the importance of the spatial component and in controlling for type I error while testing environmental predictors. Innovation We determine for the first time how the ability of commonly used (polynomial regressors) and novel methods based on eigenvector maps compare in the realm of spatial variation partitioning. We introduce a novel forward selection procedure to select spatial regressors for community analysis. Finally, we point out a number of issues that have not been previously considered about the joint explained variation between environment and space, which should be taken into account when reporting and testing the unique contributions of environment and space in patterning ecological communities. Main conclusions In tests of species‐environment relationships, spatial autocorrelation is known to inflate the level of type I error and make the tests of significance invalid. First, one must determine if the spatial component is significant using all spatial predictors (Moran's eigenvector maps). If it is, consider a model selection for the set of spatial predictors (an individual‐species forward selection procedure is to be preferred) and use the environmental and selected spatial predictors in a partial regression or partial canonical analysis scheme. This is an effective way of controlling for type I error in such tests. Polynomial regressors do not provide tests with a correct level of type I error.     

Spatial autocorrelation and the selection of simultaneous autoregressive models

Aim Spatial autocorrelation is a frequent phenomenon in ecological data and can affect estimates of model coefficients and inference from statistical models. Here, we test the performance of three different simultaneous autoregressive (SAR) model types (spatial error = SAR_err, lagged = SAR_lag and mixed = SAR_mix) and common ordinary least squares (OLS) regression when accounting for spatial autocorrelation in species distribution data using four artificial data sets with known (but different) spatial autocorrelation structures. Methods We evaluate the performance of SAR models by examining spatial patterns in model residuals (with correlograms and residual maps), by comparing model parameter estimates with true values, and by assessing their type I error control with calibration curves. We calculate a total of 3240 SAR models and illustrate how the best models [in terms of minimum residual spatial autocorrelation (minRSA), maximum model fit (R²), or Akaike information criterion (AIC)] can be identified using model selection procedures. Results Our study shows that the performance of SAR models depends on model specification (i.e. model type, neighbourhood distance, coding styles of spatial weights matrices) and on the kind of spatial autocorrelation present. SAR model parameter estimates might not be more precise than those from OLS regressions in all cases. SAR_err models were the most reliable SAR models and performed well in all cases (independent of the kind of spatial autocorrelation induced and whether models were selected by minRSA, R² or AIC), whereas OLS, SAR_lag and SAR_mix models showed weak type I error control and/or unpredictable biases in parameter estimates. Main conclusions SAR_err models are recommended for use when dealing with spatially autocorrelated species distribution data. SAR_lag and SAR_mix might not always give better estimates of model coefficients than OLS, and can thus generate bias. Other spatial modelling techniques should be assessed comprehensively to test their predictive performance and accuracy for biogeographical and macroecological research.     

INTEGRATING LANDSCAPE GENOMICS AND SPATIALLY EXPLICIT APPROACHES TO DETECT LOCI UNDER SELECTION IN CLINAL POPULATIONS

Uncovering the genetic basis of adaptation hinges on the ability to detect loci under selection. However, population genomics outlier approaches to detect selected loci may be inappropriate for clinal populations or those with unclear population structure because they require that individuals be clustered into populations. An alternate approach, landscape genomics, uses individual‐based approaches to detect loci under selection and reveal potential environmental drivers of selection. We tested four landscape genomics methods on a simulated clinal population to determine their effectiveness at identifying a locus under varying selection strengths along an environmental gradient. We found all methods produced very low type I error rates across all selection strengths, but elevated type II error rates under “weak” selection. We then applied these methods to an AFLP genome scan of an alpine plant, Campanula barbata, and identified five highly supported candidate loci associated with precipitation variables. These loci also showed spatial autocorrelation and cline patterns indicative of selection along a precipitation gradient. Our results suggest that landscape genomics in combination with other spatial analyses provides a powerful approach for identifying loci potentially under selection and explaining spatially complex interactions between species and their environment.     

A spatially explicit method for evaluating accuracy of species distribution models

Aim Models predicting the spatial distribution of animals are increasingly used in wildlife management and conservation planning. There is growing recognition that common methods of evaluating species distribution model (SDM) accuracy, as a global overall value of predictive ability, could be enhanced by spatially evaluating the model thereby identifying local areas of relative predictive strength and weakness. Current methods of spatial SDM model assessment focus on applying local measures of spatial autocorrelation to SDM residuals, which require quantitative model outputs. However, SDM outputs are often probabilistic (relative probability of species occurrence) or categorical (species present or absent). The goal of this paper was to develop a new method, using a conditional randomization technique, which can be applied to directly spatially evaluate probabilistic and categorical SDMs. Location Eastern slopes, Rocky Mountains, Alberta, Canada. Methods We used predictions from seasonal grizzly bear (Ursus arctos) resource selection functions (RSF) models to demonstrate our spatial evaluation technique. Local test statistics computed from bear telemetry locations were used to identify areas where bears were located more frequently than predicted. We evaluated the spatial pattern of model inaccuracies using a measure of spatial autocorrelation, local Moran’s I. Results We found the model to have non‐stationary patterns in accuracy, with clusters of inaccuracies located in central habitat areas. Model inaccuracies varied seasonally, with the summer model performing the best and the least error in areas with high RSF values. The landscape characteristics associated with model inaccuracies were examined, and possible factors contributing to RSF error were identified. Main conclusions The presented method complements existing spatial approaches to model error assessment as it can be used with probabilistic and categorical model output, which is typical for SDMs. We recommend that SDM accuracy assessments be done spatially and resulting accuracy maps included in model metadata.     

Model selection and information theory in geographical ecology

Aim   Although parameter estimates are not as affected by spatial autocorrelation as Type I errors, the change from classical null hypothesis significance testing to model selection under an information theoretic approach does not completely avoid problems caused by spatial autocorrelation. Here we briefly review the model selection approach based on the Akaike information criterion (AIC) and present a new routine for Spatial Analysis in Macroecology (SAM) software that helps establishing minimum adequate models in the presence of spatial autocorrelation.
Innovation    We illustrate how a model selection approach based on the AIC can be used in geographical data by modelling patterns of mammal species in South America represented in a grid system ( n  = 383) with 2° of resolution, as a function of five environmental explanatory variables, performing an exhaustive search of minimum adequate models considering three regression methods: non-spatial ordinary least squares (OLS), spatial eigenvector mapping and the autoregressive (lagged-response) model. The models selected by spatial methods included a smaller number of explanatory variables than the one selected by OLS, and minimum adequate models contain different explanatory variables, although model averaging revealed a similar rank of explanatory variables.
Main conclusions    We stress that the AIC is sensitive to the presence of spatial autocorrelation, generating unstable and overfitted minimum adequate models to describe macroecological data based on non-spatial OLS regression. Alternative regression techniques provided different minimum adequate models and have different uncertainty levels. Despite this, the averaged model based on Akaike weights generates consistent and robust results across different methods and may be the best approach for understanding of macroecological patterns.     

Uninformative Parameters and Model Selection Using Akaike's Information Criterion

Abstract: As use of Akaike's Information Criterion (AIC) for model selection has become increasingly common, so has a mistake involving interpretation of models that are within 2 AIC units (ΔAIC ≤ 2) of the top-supported model. Such models are <2 ΔAIC units because the penalty for one additional parameter is +2 AIC units, but model deviance is not reduced by an amount sufficient to overcome the 2-unit penalty and, hence, the additional parameter provides no net reduction in AIC. Simply put, the uninformative parameter does not explain enough variation to justify its inclusion in the model and it should not be interpreted as having any ecological effect. Models with uninformative parameters are frequently presented as being competitive in the Journal of Wildlife Management, including 72% of all AIC-based papers in 2008, and authors and readers need to be more aware of this problem and take appropriate steps to eliminate misinterpretation. I reviewed 5 potential solutions to this problem: 1) report all models but ignore or dismiss those with uninformative parameters, 2) use model averaging to ameliorate the effect of uninformative parameters, 3) use 95% confidence intervals to identify uninformative parameters, 4) perform all-possible subsets regression and use weight-of-evidence approaches to discriminate useful from uninformative parameters, or 5) adopt a methodological approach that allows models containing uninformative parameters to be culled from reported model sets. The first approach is preferable for small sets of a priori models, whereas the last 2 approaches should be used for large model sets or exploratory modeling.     

Accounting for multi‐scale spatial autocorrelation improves performance of invasive species distribution modelling (iSDM)

Aim Analyses of species distributions are complicated by various origins of spatial autocorrelation (SAC) in biogeographical data. SAC may be particularly important for invasive species distribution models (iSDMs) because biological invasions are strongly influenced by dispersal and colonization processes that typically create highly structured distribution patterns. We examined the efficacy of using a multi‐scale framework to account for different origins of SAC, and compared non‐spatial models with models that accounted for SAC at multiple levels. Location We modelled the spatial distribution of an invasive forest pathogen, Phytophthora ramorum, in western USA. Methods We applied one conventional statistical method (generalized linear model, GLM) and one nonparametric technique (maximum entropy, Maxent) to a large dataset on P. ramorum occurrence (n = 3787) to develop four types of model that included environmental variables and that either ignored spatial context or incorporated it at a broad scale using trend surface analysis, a local scale using autocovariates, or multiple scales using spatial eigenvector mapping. We evaluated model accuracies and amounts of explained spatial structure, and examined the changes in predictive power of the environmental and spatial variables. Results Accounting for different scales of SAC significantly enhanced the predictive capability of iSDMs. Dramatic improvements were observed when fine‐scale SAC was included, suggesting that local range‐confining processes are important in P. ramorum spread. The importance of environmental variables was relatively consistent across all models, but the explanatory power decreased in spatial models for factors with strong spatial structure. While accounting for SAC reduced the amount of residual autocorrelation for GLM but not for Maxent, it still improved the performance of both approaches, supporting our hypothesis that dispersal and colonization processes are important factors to consider in distribution models of biological invasions. Main conclusions Spatial autocorrelation has become a paradigm in biogeography and ecological modelling. In addition to avoiding the violation of statistical assumptions, accounting for spatial patterns at multiple scales can enhance our understanding of dynamic processes that explain ecological mechanisms of invasion and improve the predictive performance of static iSDMs.     

Effect of field sampling design on variation partitioning in a dendritic stream network

Variation partitioning is one of the most frequently used method to infer the importance of environmental (niche based) and spatial (dispersal) processes in metacommunity structuring. However, the reliability of the method in predicting the role of the major structuring forces is less known. We studied the effect of field sampling design on the result of variation partitioning of fish assemblages in a stream network. Along with four different sample sizes, a simple random sampling from a total of 115 stream segments (sampling objects) was applied in 400 iterations, and community variation of each random sample was partitioned into four fractions: pure environmentally (landscape variables) explained, pure spatially (MEM eigenvectors) explained, jointly explained by environment and space, and unexplained variance. Results were highly sensitive to sample size. Even at a given sample size, estimated variance fractions had remarkable random fluctuation, which can lead to inconsistent results on the relative importance of environmental and spatial variables on the structuring of metacommunities. Interestingly, all the four variance fractions correlated better with the number of the selected spatial variables than with any design properties. Sampling interval proved to be a fundamentally influential sampling design property because it affected the number of the selected spatial variables. Our findings suggest that the effect of sampling design on variation partitioning is related to the ability of the eigenvectors to model complex spatial patterns. Hence, properties of the sampling design should be more intensively considered in metacommunity studies.     

Developing approaches for linear mixed modeling in landscape genetics through landscape‐directed dispersal simulations

Dispersal can impact population dynamics and geographic variation, and thus, genetic approaches that can establish which landscape factors influence population connectivity have ecological and evolutionary importance. Mixed models that account for the error structure of pairwise datasets are increasingly used to compare models relating genetic differentiation to pairwise measures of landscape resistance. A model selection framework based on information criteria metrics or explained variance may help disentangle the ecological and landscape factors influencing genetic structure, yet there are currently no consensus for the best protocols. Here, we develop landscape‐directed simulations and test a series of replicates that emulate independent empirical datasets of two species with different life history characteristics (greater sage‐grouse; eastern foxsnake). We determined that in our simulated scenarios, AIC and BIC were the best model selection indices and that marginal R² values were biased toward more complex models. The model coefficients for landscape variables generally reflected the underlying dispersal model with confidence intervals that did not overlap with zero across the entire model set. When we controlled for geographic distance, variables not in the underlying dispersal models (i.e., nontrue) typically overlapped zero. Our study helps establish methods for using linear mixed models to identify the features underlying patterns of dispersal across a variety of landscapes.     

基于地理加权回归模型探究环境异质性对秦岭大熊猫空间利用的影响

环境异质性对野生动物分布的影响具有明显的空间不均匀性。传统分析中多采用经典线性回归模型来量化野生动物分布与环境变量之间的关系,难以准确反映物种-环境关系的空间异质特征。地理加权回归(GWR)是近年来提出的一种新的空间分析方法,通过将空间结构嵌入线性回归模型中,以此来探测空间关系的非均匀性。以秦岭大熊猫为例,应用GWR模型分析大熊猫空间分布与环境异质性特征之间的潜在关系,并同经典的全局最小二乘回归法(OLS)进行比较。结果表明,GWR模型的AIC、R~2和校正R~2均显著优于OLS模型,GWR模型的局部回归系数估计能够更加深刻地揭示大熊猫空间分布与环境变量间的复杂空间关系,且GWR模型能够为物种的科学保护提供更加有效的理论支撑。因此,GWR模型可为探究物种-环境关系的空间异质特征提供一种新的方法,在物种栖息地选择与利用研究中具有一定的应用前景。     

Application of a geographically-weighted regression analysis to estimate net primary production of Chinese forest ecosystems

Aim The objective of this paper is to obtain a net primary production (NPP) regression model based on the geographically weighted regression (GWR) method, which includes spatial non‐stationarity in the parameters estimated for forest ecosystems in China. Location We used data across China. Methods We examine the relationships between NPP of Chinese forest ecosystems and environmental variables, specifically altitude, temperature, precipitation and time‐integrated normalized difference vegetation index (TINDVI) based on the ordinary least squares (OLS) regression, the spatial lag model and GWR methods. Results The GWR method made significantly better predictions of NPP in simulations than did OLS, as indicated both by corrected Akaike Information Criterion (AIC_c) and R². GWR provided a value of 4891 for AIC_c and 0.66 for R², compared with 5036 and 0.58, respectively, by OLS. GWR has the potential to reveal local patterns in the spatial distribution of a parameter, which would be ignored by the OLS approach. Furthermore, OLS may provide a false general relationship between spatially non‐stationary variables. Spatial autocorrelation violates a basic assumption of the OLS method. The spatial lag model with the consideration of spatial autocorrelation had improved performance in the NPP simulation as compared with OLS (5001 for AIC_c and 0.60 for R²), but it was still not as good as that via the GWR method. Moreover, statistically significant positive spatial autocorrelation remained in the NPP residuals with the spatial lag model at small spatial scales, while no positive spatial autocorrelation across spatial scales can be found in the GWR residuals. Conclusions We conclude that the regression analysis for Chinese forest NPP with respect to environmental factors and based alternatively on OLS, the spatial lag model, and GWR methods indicated that there was a significant improvement in model performance of GWR over OLS and the spatial lag model.     

On the selection of phylogenetic eigenvectors for ecological analyses

Among the statistical methods available to control for phylogenetic autocorrelation in ecological data, those based on eigenfunction analysis of the phylogenetic distance matrix among the species are becoming increasingly important tools. Here, we evaluate a range of criteria to select eigenvectors extracted from a phylogenetic distance matrix (using phylogenetic eigenvector regression, PVR) that can be used to measure the level of phylogenetic signal in ecological data and to study correlated evolution. We used a principal coordinate analysis to represent the phylogenetic relationships among 209 species of Carnivora by a series of eigenvectors, which were then used to model log‐transformed body size. We first conducted a series of PVRs in which we increased the number of eigenvectors from 1 to 70, following the sequence of their associated eigenvalues. Second, we also investigated three non‐sequential approaches based on the selection of 1) eigenvectors significantly correlated with body size, 2) eigenvectors selected by a standard stepwise algorithm, and 3) the combination of eigenvectors that minimizes the residual phylogenetic autocorrelation. We mapped the mean specific component of body size to evaluate how these selection criteria affect the interpretation of non‐phylogenetic signal in Bergmann's rule. For comparison, the same patterns were analyzed using autoregressive model (ARM) and phylogenetic generalized least‐squares (PGLS). Despite the robustness of PVR to the specific approaches used to select eigenvectors, using a relatively small number of eigenvectors may be insufficient to control phylogenetic autocorrelation, leading to flawed conclusions about patterns and processes. The method that minimizes residual autocorrelation seems to be the best choice according to different criteria. Thus, our analyses show that, when the best criterion is used to control phylogenetic structure, PVR can be a valuable tool for testing hypotheses related to heritability at the species level, phylogenetic niche conservatism and correlated evolution between ecological traits.     

基于不同空间模型的马尾松林生态系统碳密度研究

以江西省马尾松林生态系统为研究对象,基于样地调查及样品碳含量测定结果计算其碳密度,并选取立地、植被及气象等方面的15个因子,采用多元线性逐步回归方法筛选出对生态系统碳密度影响显著的因子,然后分别利用最小二乘模型(OLS)、空间误差模型(SEM)、空间滞后模型(SLM)和地理加权回归模型(GWR)构建生态系统碳密度与其影响因子之间的关系模型,筛选出最优的拟合模型。结果表明:对马尾松林生态系统碳密度影响显著的因子分别为海拔、坡度、土层厚度、胸径、年均温度和年均降水量。4种模型拟合结果均显示碳密度与坡度呈负相关,与海拔、土层厚度、胸径呈正相关。模型的决定系数(R~2)由大到小分别为GWR(0.8043)SEM(0.6371)SLM(0.6364)OLS(0.6321),模型均方误差(MSE)与赤池信息准则(AIC)最大的均为OLS模型,最小的均为GWR模型;残差检验表明GWR模型能有效降低模型残差的空间自相关性。综合分析得出GWR模型的拟合效果最优,更适用于江西省马尾松林生态系统碳密度的估测。     

森林害虫发生的空间模式研究进展

森林昆虫种群表现出多样的时空模式,空间同步性是其中最常见的。回顾了森林昆虫空间同步性的特点、形成机制及研究方法方面的进展。森林害虫发生的同步性是广泛存在的,但不同昆虫种类的同步性大小不同。空间同步性常随距离的增大而下降,还与时间尺度有关。Moran效应和扩散是解释空间同步性的两种主要机制,通常Moran效应的影响要比扩散大。从虫害发生数据的获取、同步性的度量及成因3个方面介绍了空间同步性的研究方法方面的进展。利用树轮生态学原理重建森林虫害发生历史的方法可在事后获取可靠的数据,很值得国内研究者借鉴和应用。在空间自相关度量上,空间统计学方法和地统计学方法都是非常有力的手段,但由于不能处理多时间点数据而限制了其在同步性研究中的应用。在同步性成因研究中,利用变异分解将基于距离的Moran特征向量图(dbMEM)作为空间变量研究害虫发生的驱动力是比较新颖的研究方法。     

Using adaptive traits to consider potential consequences of temporal variation in selection: male guppy colour through time and space

Temporal variation in selection is typically evaluated by estimating and comparing selection coefficients in natural populations. Meta‐analyses of these coefficients have yielded important insights, but selection coefficients are limited in several respects, including low statistical power, imperfect fitness surrogates, and uncertainty regarding consequences for trait change. A complementary approach without these limitations is to examine temporal variation in adaptive traits themselves, which is mechanistically easier and more directly relevant to evolutionary consequences. We illustrate this approach by analyzing the colour patterns of male guppies, Poecilia reticulata, from each of six sites in Trinidad in each of 6 years. This system is particularly appropriate for our study because key aspects of colour variation are genetically‐based and responsive to selection. However, although spatial patterns of colour variation have been extensively considered in this system, no study has yet formally assessed annual temporal variation in non‐manipulated populations. Matching previous conclusions for the guppy system, we find that guppies from different sites manifest different colour patterns in association with different predation regimes. We here add the new finding that, although some temporal variation is present, spatial patterns of colour variation are generally consistent across years. These results suggest that, when considering adaptive traits, spatial variation is more important than temporal variation, although our study system might be exceptional in this regard. Additional studies examining spatiotemporal variation in adaptive traits could help to improve our understanding of the role that spatiotemporal variation in selection plays in the evolutionary process. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 112 , 108–122.     

Spatial autocorrelation of allozyme and quantitative markers within a natural population of Centaurea jacea (Asteraceae)

This paper compares the fine‐scale genetic structure of quantitative traits and allozyme markers within a natural population of Centaurea jacea s.l. To that end, a spatial autocorrelation approach is developed based on pairwise correlation coefficients between individuals and using sib families. Statistical properties of the proposed statistics are investigated with numerical simulations. Our results show that most quantitative traits have a significant spatial structure for their genetic component. On average, allozyme markers and the genetic component of quantitative traits have similar patterns of spatial autocorrelation that are consistent with a neutral model of isolation by distance. We also show evidence that environmental heterogeneity generates a spatial structure for the environmental component of quantitative traits. Results are discussed in terms of mechanisms generating spatial structure and are compared with those obtained on a large geographical scale.     

Modelling geographical patterns in species richness using eigenvector-based spatial filters

Aim To test the mechanisms driving bird species richness at broad spatial scales using eigenvector‐based spatial filtering. Location South America. Methods An eigenvector‐based spatial filtering was applied to evaluate spatial patterns in South American bird species richness, taking into account spatial autocorrelation in the data. The method consists of using the geographical coordinates of a region, based on eigenanalyses of geographical distances, to establish a set of spatial filters (eigenvectors) expressing the spatial structure of the region at different spatial scales. These filters can then be used as predictors in multiple and partial regression analyses, taking into account spatial autocorrelation. Autocorrelation in filters and in the regression residuals can be used as stopping rules to define which filters will be used in the analyses. Results Environmental component alone explained 8% of variation in richness, whereas 77% of the variation could be attributed to an interaction between environment and geography expressed by the filters (which include mainly broad‐scale climatic factors). Regression coefficients of environmental component were highest for AET. These results were unbiased by short‐scale spatial autocorrelation. Also, there was a significant interaction between topographic heterogeneity and minimum temperature. Conclusion Eigenvector‐based spatial filtering is a simple and suitable statistical protocol that can be used to analyse patterns in species richness taking into account spatial autocorrelation at different spatial scales. The results for South American birds are consistent with the climatic hypothesis, in general, and energy hypothesis, in particular. Habitat heterogeneity also has a significant effect on variation in species richness in warm tropical regions.     

Spatial autocorrelation analysis allows disentangling the balance between neutral and niche processes in metacommunities

One of the most popular approaches for investigating the roles of niche and neutral processes driving metacommunity patterns consists of partitioning variation in species data into environmental and spatial components. The logic is that the distance decay of similarity in communities is expected under neutral models. However, because environmental variation is often spatially structured, the decay could also be attributed to environmental factors that are missing from the analysis. Here, we use a spatial autocorrelation analysis protocol, previously developed to detect isolation‐by‐distance in allele frequencies, to evaluate patterns of species abundances under neutral dynamics. We show that this protocol can be linked with variation partitioning analyses. Moreover, in an attempt to test the neutral model, we derive three predictions to be applied both to original species abundances and to abundances predicted by a pure spatial model species abundances will be uncorrelated; Moran's I correlograms will reveal similar short‐distance autocorrelation patterns; an increasing degree of non‐neutrality will tend to generate patterns of correlation among abundances within groups of species with similar correlograms (i.e. within species with neutral and non‐neutral dynamics). We illustrate our protocol by analyzing spatial patterns in abundance of 28 terrestrially breeding anuran species from Central Amazonia. We recommend that researchers should investigate spatial autocorrelation patterns of abundances predicted by pure spatial models to identify similar patterns of spatial autocorrelation at short distances and lack of correlation between species abundances. Therefore, the hypothesis that spatial patterns in abundances are primarily due to pure neutral dynamics (rather than to missing spatiallystructured environmental factors) can be confirmed after taking environmental variables into account.     

Spatial autocorrelation and the analysis of invasion processes from distribution data: a study with the crayfish <Emphasis Type="Italic">Procambarus clarkii</Emphasis>

Complex spatial dynamics are frequent in invasive species; analyzing distribution patterns can help to understand the mechanisms driving invasions. We used different spatial regression techniques to evaluate processes determining the invasion of the red swamp crayfish Procambarus clarkii. We evaluated four a priori hypotheses on processes that may determine crayfish invasion: landscape alteration, connectivity, wetland suitability for abiotic and biotic features. We assessed the distribution of P. clarkii in 119 waterbodies in a recently invaded area. We used spatially explicit statistical techniques (spatial eigenvector mapping, generalized additive models, Bayesian intrinsic conditional autoregressive models) within an information-theoretic framework to assess the support of hypotheses; we also analyzed the pattern of spatial autocorrelation of data, model residuals, and eigenvectors. We found strong agreement between the results of spatial eigenvector mapping and Bayesian autoregressive models. Procambarus clarkii was significantly associated with the largest, permanent wetlands. Additive models suggested also association with human-dominated landscapes, but tended to overfit data. The results indicate that abiotic wetlands features and landscape alteration are major drivers of the species’ distribution. Species distribution data, residuals of ordinary least squares regression, and spatial eigenvectors all showed positive and significant spatial autocorrelation at distances up to 2,500 m; this may be caused by the dispersal ability of the species. Our analyses help to understand the processes determining the invasion and to identify the areas most at risk where screening and early management efforts can be focused. The comparison of multiple spatial techniques allows a robust assessment of factors determining complex distribution patterns.