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.     

Niche partitioning and the effect of interspecific competition on microhabitat use by two sympatric galaxiid stream fishes

1. Numerous interacting abiotic and biotic factors influence niche use and assemblage structure of freshwater fishes, but the strength of each factor changes with spatial scale. Few studies have examined the role of interspecific competition in structuring stream fish assemblages across spatial scales. We used field and laboratory approaches to examine microhabitat partitioning and the effect of interspecific competition on microhabitat use in two sympatric stream fishes (Galaxias‘southern’ and Galaxias gollumoides) at large (among streams and among sites within streams) and small (within artificial stream channels) spatial scales. 2. Diurnal microhabitat partitioning and interspecific competition at large spatial scales were analysed among three sympatry streams (streams with allotopic and syntopic sites; three separate catchments) and four allopatry streams (streams with only allotopic sites; two separate catchments). Electro‐fishing was used to sample habitat use of fishes at 30 random points within each site by quantifying four variables for each individual: water velocity, depth, distance to nearest cover and substratum size. Habitat availability was then quantified for each site by measuring those variables at each of 50 random points. Diet and stable isotope partitioning was analysed from syntopic sites only. Diel cycles of microhabitat use and interspecific competition at small spatial scales were examined by monitoring water velocity use over 48 h in artificial stream channels for three treatments: (i) allopatric G. ‘southern’ (10 G. ‘southern’); (ii) allopatric G. gollumoides (10 G. gollumoides) and (iii) sympatry (five individuals of each species). 3. One hundred and ninety‐four G. ‘southern’ and 239 G. gollumoides were sampled across all seven streams, and habitat availability between the two species was similar among all sites. Galaxias‘southern’ utilised faster water velocities than G. gollumoides in both the field and in channel experiments. Both species utilised faster water velocities in channels at night than during the day. Diet differences were observed and were supported by isotopic differences (two of three sites). No interspecific differences were observed for the other three microhabitat variables in the field, and multivariate habitat selection did not differ between species. Interspecific competition had no effect on microhabitat use of either species against any variable either in the field (large scale) or in channels (small scale). 4. The results suggest that niche partitioning occurs along a subset of microhabitat variables (water velocity use and diet). Interspecific competition does not appear to be a major biotic factor controlling microhabitat use by these sympatric taxa at any spatial scale. The results further suggest that stream fish assemblages are not primarily structured by biotic factors, reinforcing other studies de‐emphasising interspecific competition.     

Grazing and an invasive grass confound spatial pattern of exotic and native grassland plant species richness

Previous work has shown exotic and native plant species richness are negatively correlated at fine spatial scales and positively correlated at broad spatial scales. Grazing and invasive plant species can influence plant species richness, but the effects of these disturbances across spatial scales remain untested. We collected species richness data for both native and exotic plants from five spatial scales (0.5–3000 m²) in a nested, modified Whittaker plot design from severely grazed and ungrazed North American tallgrass prairie. We also recorded the abundance of an abundant invasive grass, tall fescue (Schedonorus phoenix (Scop.) Holub), at the 0.5-m² scale. We used linear mixed-effect regression to test relationships between plant species richness, tall fescue abundance, and grazing history at five spatial scales. At no scale was exotic and native species richness linearly related, but exotic species richness at all scales was greater in grazed tracts than ungrazed tracts. Native species richness declined with increasing tall fescue abundance at all five spatial scales, but exotic species richness increased with tall fescue abundance at all but the broadest spatial scales. Severe grazing did not reduce native species richness at any spatial scale. We posit that invasion of tall fescue in this working landscape of originally native grassland plants modifies species richness-spatial scale relationships observed in less disturbed systems. Tall fescue invasion constitutes a unique biotic effect on plant species richness at broad spatial scales.     

Spatial variation in the climatic predictors of species compositional turnover and endemism

Previous research focusing on broad‐scale or geographically invariant species‐environment dependencies suggest that temperature‐related variables explain more of the variation in reptile distributions than precipitation. However, species–environment relationships may exhibit considerable spatial variation contingent upon the geographic nuances that vary between locations. Broad‐scale, geographically invariant analyses may mask this local variation and their findings may not generalize to different locations at local scales. We assess how reptile–climatic relationships change with varying spatial scale, location, and direction. Since the spatial distributions of diversity and endemism hotspots differ for other species groups, we also assess whether reptile species turnover and endemism hotspots are influenced differently by climatic predictors. Using New Zealand reptiles as an example, the variation in species turnover, endemism and turnover in climatic variables was measured using directional moving window analyses, rotated through 360°. Correlations between the species turnover, endemism and climatic turnover results generated by each rotation of the moving window were analysed using multivariate generalized linear models applied at national, regional, and local scales. At national‐scale, temperature turnover consistently exhibited the greatest influence on species turnover and endemism, but model predictive capacity was low (typically r² = 0.05, P < 0.001). At regional scales the relative influence of temperature and precipitation turnover varied between regions, although model predictive capacity was also generally low. Climatic turnover was considerably more predictive of species turnover and endemism at local scales (e.g., r² = 0.65, P < 0.001). While temperature turnover had the greatest effect in one locale (the northern North Island), there was substantial variation in the relative influence of temperature and precipitation predictors in the remaining four locales. Species turnover and endemism hotspots often occurred in different locations. Climatic predictors had a smaller influence on endemism. Our results caution against assuming that variability in temperature will always be most predictive of reptile biodiversity across different spatial scales, locations and directions. The influence of climatic turnover on the species turnover and endemism of other taxa may exhibit similar patterns of spatial variation. Such intricate variation might be discerned more readily if studies at broad scales are complemented by geographically variant, local‐scale analyses.     

Combined effects of local environment and continental biogeography on the distribution of Ostracoda

Summary 1. The relationship between altitudinal gradients on small spatial scales and latitudinal gradients on broader scales has been repeatedly recognised in the biogeography of animals and plants. However, little is known about this topic in the ecology and biogeography of ostracod communities in Mediterranean flowing waters or the factors underlying these spatial patterns. 2. We analysed the ostracod assemblages of near‐natural headwater streams in the Betic and Pre‐Betic Mountains in the southern Iberian Peninsula to decipher the most important environmental gradients structuring ostracod communities on a local scale. In addition, the European altitudinal and latitudinal distributions of the most commonly found species were analysed with GIS and regression models to compare geographical effects from local to continental scales. 3. Forty sampling sites, distributed among six catchments and ranging in altitude between 150 and 1940 m a.s.l., were sampled seasonally. Limnological and geographical information was also recorded for each sample. Seventeen ostracod species were found, two of which were new findings for the Iberian Peninsula: Potamocypris fulva and Cypria reptans. The most common species were Potamocypris zschokkei, Candona neglecta, Herpetocypris brevicaudata, Cyclocypris ovum, Potamocypris villosa and Pseudocandona albicans. The distribution of these species in 918 European locations was analysed to test the hypothesised change in altitudinal distribution with varying latitude. 4. The best subset of logistic and linear regression models, selected by means of the information‐theoretic approach, found that oxygen content and the variables related with substratum and discharge were the most important variables with a negative influence on ostracod presence, abundance and species richness on a local scale. These findings suggest that the negative effect on benthic invertebrates of physical disturbances relates to high flow velocity and turbulences. 5. Multivariate ordination methods show how altitude and water chemistry are the most important variables to explain the distribution of ostracod assemblages on the small spatial scale. On a larger scale, differences in latitudinal distribution throughout Europe were significant for the six most common species found in Granada. In addition, four of these showed significant negative linear relationships between latitude and altitude in Europe, supporting the important effect of climate on local and continental scale distributions. While ostracod biogeographies are still poorly known, our results indicate the influence of Quaternary climate variability on ostracod dynamic colonisation and extinction in Europe in accordance with species‐specific temperature and water chemistry preferences.     

Habitat scale and biodiversity: influence of catchment, stream reach and bedform scales on local invertebrate diversity

Although many studies have investigated the influence of environmental patterns on local stream invertebrate diversity, there has been little consistency in reported relationships between diversity and particular environmental variables. Here we test the hypothesis that local stream invertebrate diversity is determined by a combination of factors occurring at multiple spatial scales. We developed predictive models relating invertebrate diversity (species richness and equitability) to environmental variables collected at various spatial scales (bedform, reach and catchment, respectively) using data from 97 sampling sites dispersed throughout the Taieri River drainage in New Zealand. Models based on an individual scale of perception (bedform, reach or catchment) were not able to match predictions to observations (r < 0.26, P > 0.01, between observed and predicted equitability and species richness). In contrast, models incorporating all three scales simultaneously were highly significant (P < 0.01; r = 0.55 and 0.64, between observed and predicted equitability and species richness, respectively). The most influential variables for both richness and equitability were median particle size at the bedform scale, adjacent land use at the reach scale, and relief ratio at the catchment scale. Our findings suggest that patterns observed in local assemblages are not determined solely by local mechanisms acting within assemblages, but also result from processes operating at larger spatial scales. The integration of different spatial scales may be the key to increasing model predictability and our understanding of the factors that determine local biodiversity.     

Spatial dependence of phenotype-environment associations for tadpoles in natural ponds

Within natural habitats, phenotypes are shaped by many environmental factors. Consequently, heterogeneity of these factors can promote phenotypic divergence. However, because environments exhibit heterogeneity at different spatial scales, phenotypic divergence should also exhibit such scale-dependence. Using hierarchical linear models, I determined how multiple environmental factors at two spatial scales affected the morphology of wood frog (Rana (Lithobates) sylvatica) tadpoles collected from natural ponds. Among ponds, predation risk intensity and tadpole density were strong predictors of tadpole morphology, while within ponds, other environmental variables such as water depth and leaf litter were more important. Spatial analyses revealed that water depth and leaf litter, but not predation risk intensity or tadpole density, exhibited heterogeneous spatial distributions within ponds, suggesting that spatial properties of environmental variables influenced the scale at which they shaped phenotypes. Furthermore, patterns of phenotypic variation with respect to predation risk intensity and tadpole density among ponds largely matched observations from previous laboratory studies. These results emphasize the importance of considering phenotype-environment associations across multiple spatial scales.     

Fine‐scale spatial distribution of murundus structures in the semi‐arid region of Brazil

Murundus are earth mounds widespread in most landscapes in the semi‐arid region of Brazil. Evidence obtained from predictive modelling has suggested a termite origin for these structures, opening up new opportunities for further research. Distribution of densely packed murundus at larger spatial scales is most related to climatic regime and soil nutrient availability. However, factors and processes underlying their distribution and density at smaller spatial scales are not yet fully understood. In this study, we adopted an approach based on mapping point data using high‐resolution satellite imagery, multi‐scale second‐order analysis and general linear models to examine the fine‐scale spatial distribution and density of murundus. Our results suggest that the distribution of those structures within densely packed areas is regulated by more than one process acting or interacting across multiple spatial scales. All densely packed murundus showed a significant regular distribution at the distance scale of up to 50 m radially and a completely random distribution across all other upper distance scales. We interpret the regular pattern as a result of competition for foraging territories between different termite colonies during the process formation of densely packed murundus. The random pattern at larger distance scales (above 50 m radially) can be attributed to habitat selection preferences by termite species builders of murundus mediated by local environmental resources and conditions (i.e. availability of food resources and nesting and open habitat), which would be randomly distributed in space. Thus, at finer spatial scales murundus distributions are associated with biotic interactions acting on an abiotic template. On the basis of significant linear correlations, we suggest that the density of murundus is strongly related to local temperature regime with soil‐type influencing its effect on the murundus densities. Our findings provide novel evidences that mound‐building termites are involved in the formation of murundus in the semi‐arid region of Brazil.     

Spatial factors contribute to benthic diatom structure in streams across spatial scales: Considerations for biomonitoring

Diatoms are widely used in the biological monitoring of streams because they are strong responders to environmental change, but dispersal and spatial factors can play important and potentially confounding roles in the presence, absence, and abundance of species along with characterizing species–environment relationships. To examine how spatial factors affect diatom community structure and biomonitoring, multiple scales were sampled including the Western Allegheny Plateau (n = 58), Leading Creek watershed (n = 18), and the adjacent Shade River watershed (n = 21) in southeast Ohio. Partitioning of spatial, environmental, and spatially-structured environmental variation was conducted on diatom assemblages and on diatom metrics used in biomonitoring. At the regional scale, diatom assemblages and metrics had strong relationships with agricultural (e.g., significant correlations with nutrients, conductivity, and pasture/row crops in the watershed) and alkalinity gradients. Diatom assemblages and metrics in both watersheds were strongly associated with acid mine drainage (AMD) impacts, and when spatial factors were set as covariables in CCAs, relationships with AMD gradients became even stronger, indicating the need to consider how spatial factors could reduce the strength of diatom-environment relationships. Metrics calculated at all scales had very little variation explained exclusively by spatial factors, likely because multiple species are combined into a simplified metric that reduces the effects of species dispersal. Local environmental variables accounted for 57, 42, and 42% of the total variation explained (TVE), and spatial variables accounted for 28, 31, and 37% of the TVE in the regional, Leading Creek, and Shade River datasets, respectively. The amounts of variation in diatom assemblages explained solely by spatial factors at these scales were substantial and similar to what has been reported at continental, national, and large regional (Level I Omernik ecoregions) scales (approximately 1/3 of TVE). Although amounts of variation explained are similar across scales, processes underlying the spatial structure likely differ. In addition to describing ecological patterns, recognizing the potential influence of spatial factors could improve the identification and management of environmental problems at a range of scales, as well as aid in the development of new research questions and hypotheses aimed at exploring factors that could explain portions of the spatially explicit variation.     

Predicting probability of occurrence and factors affecting distribution and abundance of three Ozark endemic crayfish species at multiple spatial scales

相似文献   

Spatial variation in adult sex ratio across multiple scales in the invasive golden apple snail,Pomacea canaliculata

Adult sex ratio (ASR) has critical effects on behavior and life history and has implications for population demography, including the invasiveness of introduced species. ASR exhibits immense variation in nature, yet the scale dependence of this variation is rarely analyzed. In this study, using the generalized multilevel models, we investigated the variation in ASR across multiple nested spatial scales and analyzed the underlying causes for an invasive species, the golden apple snail Pomacea canaliculata. We partitioned the variance in ASR to describe the variations at different scales and then included the explanatory variables at the individual and group levels to analyze the potential causes driving the variation in ASR. We firstly determined there is a significant female‐biased ASR for this species when accounting for the spatial and temporal autocorrelations of sampling. We found that, counter to nearly equal distributed variation at plot, habitat and region levels, ASR showed little variation at the town level. Temperature and precipitation at the region level were significantly positively associated with ASR, whereas the individual weight, the density characteristic, and sampling time were not significant factors influencing ASR. Our study suggests that offspring sex ratio of this species may shape the general pattern of ASR in the population level while the environmental variables at the region level translate the unbiased offspring sex ratio to the female‐biased ASR. Future research should consider the implications of climate warming on the female‐biased ASR of this invasive species and thus on invasion pattern.     

Genetic Analysis of Invasive Plant Populations at Different Spatial Scales

Measuring genetic diversity requires selection of a spatial scale of analysis. Different levels of genetic structuring are revealed at different spatial scales, however, and the relative importance of factors driving genetic structuring varies along the spatial scale continuum. Unequal gene flow is a major factor determining genetic structure in plant populations at the local level, while the effect of selection imposed by environmental heterogeneity increases with the spatial scale of analysis. At a continental and global scale genetic structure of invasive plant populations is significantly affected by founder effect and propagule transport via human vectors. Although genetic analysis at one spatial scale provides only partial information about the invasion process, little published research reports such data for the same species at multiple scales. A multi-faceted approach to investigating the genetic structure of invasive plant species that incorporates sampling at different spatial and temporal scales would provide a more complete picture of the role of genetic forces in invasion.     

Similar biotic factors affect early establishment and abundance of an invasive plant species across spatial scales

Research in community invasibiliy has focused on biotic and abiotic factors that influence the establishment of invasive species and whether such factors vary with spatial scale. Here, we investigate the role of both biotic and abiotic factors associated with the initial establishment of Lespedeza cuneata (L. cuneata) and its abundance at three spatial scales: neighborhoods (9-m² plots), communities (50-m² transect) and old fields (5,000–70,000 m²). We asked: (1) Do resource availability and community structure affect the establishment of L. cuneata?, and (2) Are resource availability and community structure associated with patterns of L. cuneata abundance from neighborhood scales to old-field scales? To investigate the first question, we manipulated soil nitrogen (N) availability at three levels in an existing old-field community and tracked emergence and persistence of L. cuneata seedlings, as well as total plant biomass of the community, availability of light, and soil moisture content. To address the second question, we performed surveys in which we estimated L. cuneata foliar cover at community scales (50-m² belt transects) and old-field scales (total area of 28 ha), and assessed the same biotic and abiotic variables as in the field experiment. The experiment revealed that establishment and persistence by L. cuneata seedlings were 15× and 5× lower in N-added plots than in N-reduced plots. Total plant community biomass was 30% greater in N-added plots than in N-reduced plots. Conversely, light and soil moisture were 60 and 20% lower in N-added plots than in N-reduced plots. Surveys of old fields indicated that community biomass was positively associated with L. cuneata cover at old-field scales likely resulting from greater soil N input from nitrogen fixation in fields with greater L. cuneata cover. In sum, these results indicate that biotic factors associated with establishment of a Rank 1 invasive plant species at the community scale are also related to its distribution at the old-field scale, but the direction of such associations changed across scales.     

Hierarchical scales of spatial variation in the smaller surface and near-surface macrobenthos of a subtropical intertidal seagrass system in Moreton Bay,Queensland

Core samples were taken along a 4 km stretch of intertidal seagrass on North Stradbroke Island, eastern Australia, at nested scales of 1 m (stations), 150 m (sites), and 2 km (localities) to investigate the extent to which abundance, diversity, and assemblage composition of the dominant smaller members (<10 mm) of the intertidal seagrass macrobenthos vary spatially and over what scales. Gastropods and polychaetes dominated both the 91 species present and, together with decapods, also the numbers of individuals. Abundance was low (mean < 2000 individuals m⁻²) but species diversity was high (overall Simpson’s index of diversity 0.91), with 44% of species occurring only as one or two individuals, and with only two species contributing >10% to the total numbers (the microgastropod Calopia imitata and crab Enigmaplax littoralis, both little known, rarely recorded endemics). On average, a species only occurred at 6% of stations and only four occurred at >25%. Assemblages at the three localities did not vary significantly in gross ecological features (levels of species richness, faunal abundance and species diversity per component site) (ANOVA P ≫ 0.05), but did vary markedly in their composition at all spatial scales (PERMANOVA P < 0.05). Variance partitioning showed that components of total variance were least at the largest spatial scale (locality 15.9%) and greatest at the smallest scale (station 59.3%). The commoner individual species all showed random distributions at small spatial scales but clumped distributions at large spatial scales.     

Multivariate geostatistical methods for analysis of relationships between ecological indicators and environmental factors at multiple spatial scales

As all biodiversity-related variables, ecological indicators are influenced by environmental factors working at different spatial scales. However, assessing the relationship between environmental factors and ecological indicators is limited to a set of spatial scales determined a priori. This a priori assumption can hide important relationships, especially for ecological indicators with a complex spatial structure that can be driven, for example, by the influence of multiple pollutants with different dispersion ranges or by the influence of local and regional factors such as land-cover and climate. To relate ecological indicators and environmental factors without assuming a priori spatial scales of analysis, we used a Linear Model of Coregionalization. This method has been used in literature to analyze the joint distribution of biodiversity variables. Here we show that it can be used to gain insight into spatial patterns of relationships between ecological indicators and underlying environmental factors. We applied this method to a region of south-west Europe, relating data from land-cover, altitude and climate with an ecological indicator, the abundance of fruticose lichen species, known to be very sensitive to multiple environmental factors. Based on variogram analysis we identified distinct spatial scales of relationships between the ecological indicator and environmental factors. For each spatial scale we described relationships using Principal Component Analysis applied to the coregionalization matrices. This way we could assess how strong the relationship between each environmental factor and ecological indicator at each spatial scale was: at medium scales (c. 15 km) open spaces areas (a proxy for particle emissions) were more important; at larger scales (c. 45 km) open spaces, artificial areas (a proxy for gaseous pollutants) and also climate were preponderant. Thus, multivariate geostatistics provided a tool to improve knowledge on relationships between ecological indicators and environmental factors at multiple spatial scales without setting a priori spatial scales of analysis.     

Geographical ecology of the palms (Arecaceae): determinants of diversity and distributions across spatial scales

Background

The palm family occurs in all tropical and sub-tropical regions of the world. Palms are of high ecological and economical importance, and display complex spatial patterns of species distributions and diversity.

Scope

This review summarizes empirical evidence for factors that determine palm species distributions, community composition and species richness such as the abiotic environment (climate, soil chemistry, hydrology and topography), the biotic environment (vegetation structure and species interactions) and dispersal. The importance of contemporary vs. historical impacts of these factors and the scale at which they function is discussed. Finally a hierarchical scale framework is developed to guide predictor selection for future studies.

Conclusions

Determinants of palm distributions, composition and richness vary with spatial scale. For species distributions, climate appears to be important at landscape and broader scales, soil, topography and vegetation at landscape and local scales, hydrology at local scales, and dispersal at all scales. For community composition, soil appears important at regional and finer scales, hydrology, topography and vegetation at landscape and local scales, and dispersal again at all scales. For species richness, climate and dispersal appear to be important at continental to global scales, soil at landscape and broader scales, and topography at landscape and finer scales. Some scale–predictor combinations have not been studied or deserve further attention, e.g. climate on regional to finer scales, and hydrology and topography on landscape and broader scales. The importance of biotic interactions – apart from general vegetation structure effects – for the geographic ecology of palms is generally underexplored. Future studies should target scale–predictor combinations and geographic domains not studied yet. To avoid biased inference, one should ideally include at least all predictors previously found important at the spatial scale of investigation.     

Species turnover in Amazonian frogs: low predictability and large differences among forests

The factors explaining species turnover at different spatial scales have been intensively studied, but most work in Amazonia has mainly focused on plants. For animals, it is not as obvious which environmental variables most affect differences in species composition among sites. We sought to identify what causes anuran turnover in Amazonian terra firme forests, and how the effectiveness of these factors varies among regions and across spatial scales. We sampled frogs in 56 plots along ponds and streams distributed over three terra firme forest areas in Eastern Amazonia. Using multiple regressions on distance matrices, we partitioned the variation in species turnover into components explained by variation in environmental and spatial distances. This was done in parallel for each area separately, and for all areas together, to assess the consistency of results between scales and across areas at the same scale. Each community seemed to respond to a set of factors specific to that area, and the identity of the variables that emerged as significant differed among areas and scales. Both geographical distances and environmental differences had larger explanatory power at the regional scale than at the local scale. The large differences among results from different areas caution against making broad generalizations about species turnover patterns from a single community, as real differences may exist among areas.     

Invasion genetics of the Ciona intestinalis species complex: from regional endemism to global homogeneity

Determining the degree of population connectivity and investigating factors driving genetic exchange at various geographical scales are essential to understanding population dynamics and spread potential of invasive species. Here, we explore these issues in the highly invasive vase tunicate, Ciona intestinalis, a species whose invasion history has been obscured by its poorly understood taxonomy and population genetics. Recent phylogenetic and comparative genomic studies suggest that C. intestinalis is a cryptic species complex consisting of at least three species. We reconstructed phylogenies based on both mitochondrial (cytochrome c oxidase subunit 3—NADH dehydrogenase subunit 1 region and NADH dehydrogenase subunit 4 gene) and nuclear (internal transcribed spacer 1) sequences, results of which support four major phylogroups corresponding to the previously reported spA, spB and Ciona spp. (spC) as well as an undescribed cryptic species (spD). While spC and spD remain restricted to their native ranges in the Mediterranean Sea and Black Sea, respectively, the highly invasive species (spA and spB) have disjunct global distributions. Despite extensive interspecific divergences, we identified low phylogeographical structure within these two invasive species. Haplotype network analyses revealed comparatively limited mutation steps among haplotypes within each species. Population genetic analyses based on two mtDNA fragments and eight unlinked microsatellites illustrated relatively low population differentiation and high population connectivity at both regional and continental scales in the two invasive species. Human‐mediated dispersal coupled with a high potential for natural dispersal is probably responsible for the observed genetic homogeneity.