Estimating non-native plant richness with a species-accumulation model along roads

Monitoring non-native plant richness is important for biodiversity conservation and scientific research. The species-area model (SA model) has been used frequently to estimate the total species richness within a region. However, the conventional SA model may not provide robust estimations of non-native plant richness because the ecological processes associated with the accumulation of exotic and native plants may differ. Because roads strongly dictate the distributions of exotic plants, we propose a species-accumulation model along roads (SR model), rather than an SA model, to estimate the non-native plant richness within a region. Using 270 simulated data sets, we compared the differences in performance between the SR and SA models. A decision tree based on prediction accuracy was created to guide model application, which was validated using field data from 3 national nature reserves in 3 different provinces in China. The SR model significantly outperformed the SA model when non-native species were restricted to the roadsides and the proportion of uncommon exotic species was small. More importantly, the SR model accurately estimated the non-native plant richness in all field sites with an error of <1 species per site. We believe our new model meets the practical need to efficiently and robustly estimate non-native plant richness, which may facilitate effective biodiversity conservations and promote research on non-native plant invasion and vegetation dynamics.     

Improving effectiveness of systematic conservation planning with density data

Systematic conservation planning aims to design networks of protected areas that meet conservation goals across large landscapes. The optimal design of these conservation networks is most frequently based on the modeled habitat suitability or probability of occurrence of species, despite evidence that model predictions may not be highly correlated with species density. We hypothesized that conservation networks designed using species density distributions more efficiently conserve populations of all species considered than networks designed using probability of occurrence models. To test this hypothesis, we used the Zonation conservation prioritization algorithm to evaluate conservation network designs based on probability of occurrence versus density models for 26 land bird species in the U.S. Pacific Northwest. We assessed the efficacy of each conservation network based on predicted species densities and predicted species diversity. High‐density model Zonation rankings protected more individuals per species when networks protected the highest priority 10‐40% of the landscape. Compared with density‐based models, the occurrence‐based models protected more individuals in the lowest 50% priority areas of the landscape. The 2 approaches conserved species diversity in similar ways: predicted diversity was higher in higher priority locations in both conservation networks. We conclude that both density and probability of occurrence models can be useful for setting conservation priorities but that density‐based models are best suited for identifying the highest priority areas. Developing methods to aggregate species count data from unrelated monitoring efforts and making these data widely available through ecoinformatics portals such as the Avian Knowledge Network will enable species count data to be more widely incorporated into systematic conservation planning efforts.     

The Role of Roadsides in Plant Invasions: a Demographic Approach

Abstract:  Non-native plant species are common along roadsides, but presence does not necessarily indicate spread along the road axis. Roadsides may serve merely as habitat for a species spreading independently of roads. The potential conduit function of roads depends on the habitat specificity of the spreading species, its dispersal range relative to the spacing of roads in the landscape, and the relative importance of long- and short-range dispersal. We describe a demographic model of the road × species interaction and suggest methods of assessing conduit function in the field based on the model results. A species limited to roadside habitat will be constrained to spread along the road axis unless its long-range dispersal is sufficient to carry it across the intervening unfavorable area to another road. It will propagate along a road corridor at a rate determined by the scale of short-range dispersal. Effective management of an invasion requires distinguishing between the habitat and conduit functions, a distinction difficult to make with only snapshot data. Invasions can be reconstructed by several methods, but none is totally satisfactory. We suggest comparing stem distributions on transects parallel and perpendicular to the road axis, and beside the road, and away from it, with an idealized Gaussian curve. Such comparisons would allow discrimination between pattern determined by habitat suitability and pattern reflecting random and facilitated dispersal.     

Assumption‐versus data‐based approaches to summarizing species’ ranges

For conservation decision making, species’ geographic distributions are mapped using various approaches. Some such efforts have downscaled versions of coarse‐resolution extent‐of‐occurrence maps to fine resolutions for conservation planning. We examined the quality of the extent‐of‐occurrence maps as range summaries and the utility of refining those maps into fine‐resolution distributional hypotheses. Extent‐of‐occurrence maps tend to be overly simple, omit many known and well‐documented populations, and likely frequently include many areas not holding populations. Refinement steps involve typological assumptions about habitat preferences and elevational ranges of species, which can introduce substantial error in estimates of species’ true areas of distribution. However, no model‐evaluation steps are taken to assess the predictive ability of these models, so model inaccuracies are not noticed. Whereas range summaries derived by these methods may be useful in coarse‐grained, global‐extent studies, their continued use in on‐the‐ground conservation applications at fine spatial resolutions is not advisable in light of reliance on assumptions, lack of real spatial resolution, and lack of testing. In contrast, data‐driven techniques that integrate primary data on biodiversity occurrence with remotely sensed data that summarize environmental dimensions (i.e., ecological niche modeling or species distribution modeling) offer data‐driven solutions based on a minimum of assumptions that can be evaluated and validated quantitatively to offer a well‐founded, widely accepted method for summarizing species’ distributional patterns for conservation applications.     

An integrated approach for predicting fates of reintroductions with demographic data from multiple populations

We devised a novel approach to model reintroduced populations whereby demographic data collected from multiple sites are integrated into a Bayesian hierarchical model. Integrating data from multiple reintroductions allows more precise population-growth projections to be made, especially for populations for which data are sparse, and allows projections that account for random site-to-site variation to be made before new reintroductions are attempted. We used data from reintroductions of the North Island Robin (Petroica longipes), an endemic New Zealand passerine, to 10 sites where non-native mammalian predators are controlled. A comparison of candidate models that we based on deviance information criterion showed that rat-tracking rate (an index of rat density) was a useful predictor of robin fecundity and adult female survival, that landscape connectivity and a binary measure of whether sites were on a peninsula were useful predictors of apparent juvenile survival (probably due to differential dispersal away from reintroduction sites), and that there was unexplained random variation among sites in all demographic rates. We used the two best supported models to estimate the finite rate of increase (λ) for populations at each of the 10 sites, and for a proposed reintroduction site, under different levels of rat control. Only three of the reintroduction sites had λ distributions completely >1 for either model. At two sites, λ was expected to be >1 if rat-tracking rates were <5%. At the other five reintroduction sites, λ was predicted to be close to 1, and it was unclear whether growth was expected. Predictions of λ for the proposed reintroduction site were less precise than for other sites because distributions incorporated the full range of site-to-site random variation in vital rates. Our methods can be applied to any species for which postrelease data on demographic rates are available and potentially can be extended to model multiple species simultaneously.     

Use of Community‐Composition Data to Predict the Fecundity and Abundance of Species

Abstract: Species distribution models are critical tools for the prediction of invasive species spread and conservation of biodiversity. The majority of species distribution models have been built with environmental data. Community ecology theory suggests that species co‐occurrence data could also be used to predict current and potential distributions of species. Species assemblages are the products of biotic and environmental constraints on the distribution of individual species and as a result may contain valuable information for niche modeling. We compared the predictive ability of distribution models of annual grassland plants derived from either environmental or community‐composition data. Composition‐based models were built with the presence or absence of species at a site as predictors of site quality, whereas environment‐based models were built with soil chemistry, moisture content, above‐ground biomass, and solar radiation as predictors. The reproductive output of experimentally seeded individuals of 4 species and the abundance of 100 species were used to evaluate the resulting models. Community‐composition data were the best predictors of both the site‐specific reproductive output of sown individuals and the site‐specific abundance of existing populations. Successful community‐based models were robust to omission of data on the occurrence of rare species, which suggests that even very basic survey data on the occurrence of common species may be adequate for generating such models. Our results highlight the need for increased public availability of ecological survey data to facilitate community‐based modeling at scales relevant to conservation.     

Use of Coarse‐Resolution Models of Species’ Distributions to Guide Local Conservation Inferences

Abstract: Distribution models are used increasingly for species conservation assessments over extensive areas, but the spatial resolution of the modeled data and, consequently, of the predictions generated directly from these models are usually too coarse for local conservation applications. Comprehensive distribution data at finer spatial resolution, however, require a level of sampling that is impractical for most species and regions. Models can be downscaled to predict distribution at finer resolutions, but this increases uncertainty because the predictive ability of models is not necessarily consistent beyond their original scale. We analyzed the performance of downscaled, previously published models of environmental favorability (a generalized linear modeling technique) for a restricted endemic insectivore, the Iberian desman (Galemys pyrenaicus), and a more widespread carnivore, the Eurasian otter (Lutra lutra), in the Iberian Peninsula. The models, built from presence–absence data at 10 × 10 km resolution, were extrapolated to a resolution 100 times finer (1 × 1 km). We compared downscaled predictions of environmental quality for the two species with published data on local observations and on important conservation sites proposed by experts. Predictions were significantly related to observed presence or absence of species and to expert selection of sampling sites and important conservation sites. Our results suggest the potential usefulness of downscaled projections of environmental quality as a proxy for expensive and time‐consuming field studies when the field studies are not feasible. This method may be valid for other similar species if coarse‐resolution distribution data are available to define high‐quality areas at a scale that is practical for the application of concrete conservation measures.     

How decisions about fitting species distribution models affect conservation outcomes

Species distribution models (SDMs) are increasingly used in conservation and land-use planning as inputs to describe biodiversity patterns. These models can be built in different ways, and decisions about data preparation, selection of predictor variables, model fitting, and evaluation all alter the resulting predictions. Commonly, the true distribution of species is unknown and independent data to verify which SDM variant to choose are lacking. Such model uncertainty is of concern to planners. We analyzed how 11 routine decisions about model complexity, predictors, bias treatment, and setting thresholds for predicted values altered conservation priority patterns across 25 species. Models were created with MaxEnt and run through Zonation to determine the priority rank of sites. Although all SDM variants performed well (area under the curve >0.7), they produced spatially different predictions for species and different conservation priority solutions. Priorities were most strongly altered by decisions to not address bias or to apply binary thresholds to predicted values; on average 40% and 35%, respectively, of all grid cells received an opposite priority ranking. Forcing high model complexity altered conservation solutions less than forcing simplicity (14% and 24% of cells with opposite rank values, respectively). Use of fewer species records to build models or choosing alternative bias treatments had intermediate effects (25% and 23%, respectively). Depending on modeling choices, priority areas overlapped as little as 10–20% with the baseline solution, affecting top and bottom priorities differently. Our results demonstrate the extent of model-based uncertainty and quantify the relative impacts of SDM building decisions. When it is uncertain what the best SDM approach and conservation plan is, solving uncertainty or considering alterative options is most important for those decisions that change plans the most.     

Extrapolating cetacean densities to quantitatively assess human impacts on populations in the high seas

As human activities expand beyond national jurisdictions to the high seas, there is an increasing need to consider anthropogenic impacts to species inhabiting these waters. The current scarcity of scientific observations of cetaceans in the high seas impedes the assessment of population‐level impacts of these activities. We developed plausible density estimates to facilitate a quantitative assessment of anthropogenic impacts on cetacean populations in these waters. Our study region extended from a well‐surveyed region within the U.S. Exclusive Economic Zone into a large region of the western North Atlantic sparsely surveyed for cetaceans. We modeled densities of 15 cetacean taxa with available line transect survey data and habitat covariates and extrapolated predictions to sparsely surveyed regions. We formulated models to reduce the extent of extrapolation beyond covariate ranges, and constrained them to model simple and generalizable relationships. To evaluate confidence in the predictions, we mapped where predictions were made outside sampled covariate ranges, examined alternate models, and compared predicted densities with maps of sightings from sources that could not be integrated into our models. Confidence levels in model results depended on the taxon and geographic area and highlighted the need for additional surveying in environmentally distinct areas. With application of necessary caution, our density estimates can inform management needs in the high seas, such as the quantification of potential cetacean interactions with military training exercises, shipping, fisheries, and deep‐sea mining and be used to delineate areas of special biological significance in international waters. Our approach is generally applicable to other marine taxa and geographic regions for which management will be implemented but data are sparse.     

Accuracy of Short‐Term Demographic Data in Projecting Long‐Term Fate of Populations

Short‐term surveys are useful in conservation of species if they can be used to reliably predict the long‐term fate of populations. However, statistical evaluations of reliability are rare. We studied how well short‐term demographic data (1999–2002) of tartar catchfly (Silene tatarica), a perennial riparian plant, projected the fate and growth of 23 populations of this species up to the year 2010. Surveyed populations occurred along a river with natural flood dynamics and along a regulated river. Riparian plant populations are affected by flooding, which maintains unvegetated shores, while forest succession proceeds in areas with little flooding. Flooding is less severe along the regulated river, and vegetation overgrowth reduces abundance of tartar catchfly on unvegetated shores. We built matrix models to calculate population growth rates and estimated times to population extinction in natural and in regulated rivers, 13 and 10 populations, respectively. Models predicted population survival well (model predictions matched observed survival in 91% of populations) and accurately predicted abundance increases and decreases in 65% of populations. The observed and projected population growth rates differed significantly in all but 3 populations. In most cases, the model overestimated population growth. Model predictions did not improve when data from more years were used (1999–2006). In the regulated river, the poorest model predictions occurred in areas where cover of other plant species changed the fastest. Although vegetation cover increased in most populations, it decreased in 4 populations along the natural river. Our results highlight the need to combine disturbance and succession dynamics in demographic models and the importance of habitat management for species survival along regulated rivers. Precisión de Datos Demográficos de Corto Plazo en la Proyección del Destino de Poblaciones a Largo Plazo     

Validating graph-based connectivity models with independent presence–absence and genetic data sets

Habitat connectivity is a key objective of current conservation policies and is commonly modeled by landscape graphs (i.e., sets of habitat patches [nodes] connected by potential dispersal paths [links]). These graphs are often built based on expert opinion or species distribution models (SDMs) and therefore lack empirical validation from data more closely reflecting functional connectivity. Accordingly, we tested whether landscape graphs reflect how habitat connectivity influences gene flow, which is one of the main ecoevolutionary processes. To that purpose, we modeled the habitat network of a forest bird (plumbeous warbler [Setophaga plumbea]) on Guadeloupe with graphs based on expert opinion, Jacobs’ specialization indices, and an SDM. We used genetic data (712 birds from 27 populations) to compute local genetic indices and pairwise genetic distances. Finally, we assessed the relationships between genetic distances or indices and cost distances or connectivity metrics with maximum-likelihood population-effects distance models and Spearman correlations between metrics. Overall, the landscape graphs reliably reflected the influence of connectivity on population genetic structure; validation R² was up to 0.30 and correlation coefficients were up to 0.71. Yet, the relationship among graph ecological relevance, data requirements, and construction and analysis methods was not straightforward because the graph based on the most complex construction method (species distribution modeling) sometimes had less ecological relevance than the others. Cross-validation methods and sensitivity analyzes allowed us to make the advantages and limitations of each construction method spatially explicit. We confirmed the relevance of landscape graphs for conservation modeling but recommend a case-specific consideration of the cost-effectiveness of their construction methods. We hope the replication of independent validation approaches across species and landscapes will strengthen the ecological relevance of connectivity models.     

A Behaviorally Explicit Demographic Model Integrating Habitat Selection and Population Dynamics in California Sea Lions

Abstract: Although there has been a call for the integration of behavioral ecology and conservation biology, there are few tools currently available to achieve this integration. Explicitly including information about behavioral strategies in population viability analyses may enhance the ability of conservation biologists to understand and estimate patterns of extinction risk. Nevertheless, most behavioral‐based PVA approaches require detailed individual‐based data that are rarely available for imperiled species. We present a mechanistic approach that incorporates spatial and demographic consequences of behavioral strategies into population models used for conservation. We developed a stage‐structured matrix model that includes the costs and benefits of movement associated with 2 habitat‐selection strategies (philopatry and direct assessment). Using a life table for California sea lions (Zalophus californianus), we explored the sensitivity of model predictions to the inclusion of these behavioral parameters. Including behavioral information dramatically changed predicted population sizes, model dynamics, and the expected distribution of individuals among sites. Estimated population sizes projected in 100 years diverged up to 1 order of magnitude among scenarios that assumed different movement behavior. Scenarios also exhibited different model dynamics that ranged from stable equilibria to cycles or extinction. These results suggest that inclusion of behavioral data in viability models may improve estimates of extinction risk for imperiled species. Our approach provides a simple method for incorporating spatial and demographic consequences of behavioral strategies into population models and may be easily extended to other species and behaviors to understand the mechanisms of population dynamics for imperiled populations.     

Predicting Bird Species Distributions in Reconstructed Landscapes

Abstract:  Landscape optimization for biodiversity requires prediction of species distributions under alternative revegetation scenarios. We used Bayesian model averaging with logistic regression to predict probabilities of occurrence for 61 species of birds within highly fragmented box–ironbark forests of central Victoria, Australia. We used topographic, edaphic, and climatic variables as predictors so that the models could be applied to areas where vegetation has been cleared but may be replanted. Models were evaluated with newly acquired, independent data collected in large blocks of remnant native vegetation. Successful predictions were obtained for 18 of 45 woodland species (40%). Model averaging produced more accurate predictions than "single best" models. Models were most successful for smaller-bodied species that probably depend on particular vegetation types. Predictions for larger, generalist species, and seasonal migrants were less successful, partly because of changes in species distributions between model building (1995–1997) and validation (2004–2005) surveys. We used validated models to project occurrence probabilities for individual species across a 12,000-km² region, assuming native vegetation was present. These predictions are intended to be used as inputs, along with landscape context and temporal dynamics, into optimization algorithms to prioritize revegetation. Longer-term data sets to accommodate temporal dynamics are needed to improve the predictive accuracy of models.     

Assessing the effects of pseudo-absences on predictive distribution model performance

Modelling species distributions with presence data from atlases, museum collections and databases is challenging. In this paper, we compare seven procedures to generate pseudo-absence data, which in turn are used to generate GLM-logistic regressed models when reliable absence data are not available. We use pseudo-absences selected randomly or by means of presence-only methods (ENFA and MDE) to model the distribution of a threatened endemic Iberian moth species (Graellsia isabelae). The results show that the pseudo-absence selection method greatly influences the percentage of explained variability, the scores of the accuracy measures and, most importantly, the degree of constraint in the distribution estimated. As we extract pseudo-absences from environmental regions further from the optimum established by presence data, the models generated obtain better accuracy scores, and over-prediction increases. When variables other than environmental ones influence the distribution of the species (i.e., non-equilibrium state) and precise information on absences is non-existent, the random selection of pseudo-absences or their selection from environmental localities similar to those of species presence data generates the most constrained predictive distribution maps, because pseudo-absences can be located within environmentally suitable areas. This study shows that if we do not have reliable absence data, the method of pseudo-absence selection strongly conditions the obtained model, generating different model predictions in the gradient between potential and realized distributions.     

Forest Roads as Partial Barriers to Terrestrial Salamander Movement

Abstract:  Roads can fragment animal populations by disrupting movement among formerly continuous habitats. Although models have demonstrated that disrupted movement can contribute to long-term extinction, there are few empirical data on the effects of roads on animal movement. We used displacement and homing experiments to determine whether forest roads are barriers to the movement of terrestrial salamanders. We displaced 1471 red-backed salamanders ( Plethodon cinereus ) across five forest roads and compared return rates to those of salamanders displaced equal distances toward the forest interior. Roads significantly reduced the return rate of salamanders, with a mean reduction of 51%. Steep roadside verges further reduced return rates, particularly for salamanders moving downhill across verges. The permeability of roads to salamander movement did not appear to be related to road surface type. Gravel roads had both the highest and lowest observed permeability with the two paved roads intermediate between these. We conclude that narrow forest roads are partial barriers to salamander movement and that steep roadside verges may exacerbate these effects.     

Ecological‐Niche Modeling and Prioritization of Conservation‐Area Networks for Mexican Herpetofauna

Abstract: One of the most important tools in conservation biology is information on the geographic distribution of species and the variables determining those patterns. We used maximum‐entropy niche modeling to run distribution models for 222 amphibian and 371 reptile species (49% endemics and 27% threatened) for which we had 34,619 single geographic records. The planning region is in southeastern Mexico, is 20% of the country's area, includes 80% of the country's herpetofauna, and lacks an adequate protected‐area system. We used probabilistic data to build distribution models of herpetofauna for use in prioritizing conservation areas for three target groups (all species and threatened and endemic species). The accuracy of species‐distribution models was better for endemic and threatened species than it was for all species. Forty‐seven percent of the region has been deforested and additional conservation areas with 13.7% to 88.6% more native vegetation (76% to 96% of the areas are outside the current protected‐area system) are needed. There was overlap in 26 of the main selected areas in the conservation‐area network prioritized to preserve the target groups, and for all three target groups the proportion of vegetation types needed for their conservation was constant: 30% pine and oak forests, 22% tropical evergreen forest, 17% low deciduous forest, and 8% montane cloud forests. The fact that different groups of species require the same proportion of habitat types suggests that the pine and oak forests support the highest proportion of endemic and threatened species and should therefore be given priority over other types of vegetation for inclusion in the protected areas of southeastern Mexico.     

Distribution of a Community of Mammals in Relation to Roads and Other Human Disturbances in Gabon,Central Africa

We present the first community‐level study of the associations of both roads and other human disturbances with the distribution of mammals in Gabon (central Africa). Our study site was in an oil concession within a littoral mosaic landscape. We conducted surveys along 199 line transects and installed camera traps on 99 of these transects to document mammal presence and abundance. We used generalized linear mixed‐effect models to document associations between variables related to the ecosystem (land cover, topography, and hydrology), roads (coating, width of rights of way, condition, type of vehicle used on the road, traffic level, affiliation of users, and general type of road), and other human disturbances (urbanization, agriculture, hunting, logging, gathering, and industrial activities) and the abundance or presence of 17 species or groups of mammals including elephant (Loxodonta cyclotis), buffalo (Syncerus caffer), sitatunga (Tragelaphus spekei), red river hog (Potamochoerus porcus), smaller ungulates, gorilla (Gorilla gorilla), chimpanzee (Pan troglodytes), side‐striped jackal (Canis adustus), carnivores, monkeys, and large rodents. Some types of roads and other human disturbances were negatively associated with the abundance or presence of elephants, buffalos, gorillas, sitatungas, some monkeys, and duikers. The pattern of associations of mammals with roads and other human disturbances was diverse and included positive associations with road presence (red river hog, some monkeys, and duikers), agriculture (sitatunga, small carnivores, and large rodents) and industrial activities (sitatunga, red river hog, red duikers, and side‐striped jackal). Our results suggest that the community of mammals we studied was mostly affected by hunting, agriculture, and urbanization, which are facilitated by road presence. We recommend increased regulation of agriculture, hunting, and road building in the area. Distribución de una Comunidad de Mamíferos en Relación a Carreteras y Otras Perturbaciones Humanas en Gabón, Africa Central     

Integrated models to support multiobjective ecological restoration decisions

Many objectives motivate ecological restoration, including improving vegetation condition, increasing the range and abundance of threatened species, and improving species richness and diversity. Although models have been used to examine the outcomes of ecological restoration, few researchers have attempted to develop models to account for multiple, potentially competing objectives. We developed a combined state‐and‐transition, species‐distribution model to predict the effects of restoration actions on vegetation condition and extent, bird diversity, and the distribution of several bird species in southeastern Australian woodlands. The actions reflected several management objectives. We then validated the models against an independent data set and investigated how the best management decision might change when objectives were valued differently. We also used model results to identify effective restoration options for vegetation and bird species under a constrained budget. In the examples we evaluated, no one action (improving vegetation condition and extent, increasing bird diversity, or increasing the probability of occurrence for threatened species) provided the best outcome across all objectives. In agricultural lands, the optimal management actions for promoting the occurrence of the Brown Treecreeper (Climacteris picumnus), an iconic threatened species, resulted in little improvement in the extent of the vegetation and a high probability of decreased vegetation condition. This result highlights that the best management action in any situation depends on how much the different objectives are valued. In our example scenario, no management or weed control were most likely to be the best management options to satisfy multiple restoration objectives. Our approach to exploring trade‐offs in management outcomes through integrated modeling and structured decision‐support approaches has wide application for situations in which trade‐offs exist between competing conservation objectives.     

Northward shifts of the distributions of Spanish reptiles in association with climate change

It is predicted that climate change will drive extinctions of some reptiles and that the number of these extinctions will depend on whether reptiles are able to change their distribution. Whether the latitudinal distribution of reptiles may change in response to increases in temperature is unknown. We used data on reptile distributions collected during the 20th century to analyze whether changes in the distributions of reptiles in Spain are associated with increases in temperature. We controlled for biases in sampling effort and found a mean, statistically significant, northward shift of the northern extent of reptile distributions of about 15.2 km from 1940-1975 to 1991-2005. The southern extent of the distributions did not change significantly. Thus, our results suggest that the latitudinal distributions of reptiles may be changing in response to climate change.     

Effects of vegetation cover, presence of a native ant species, and human disturbance on colonization by Argentine ants

The spread of non-native invasive species is affected by human activity, vegetation cover, weather, and interaction with native species. We analyzed data from a 17-year study of the distribution of the non-native Argentine ant (Linepithema humile) and the native winter ant (Prenolepis imparis) in a preserve in northern California (U.S.A.). We conducted logistic regressions and used model selection to determine whether the following variables were associated with changes in the distribution of each species: presence of conspecifics at neighboring sites, distance to development (e.g., roads, buildings, and landscaped areas), proportion of vegetation cover taller than 0.75 m, elevation, distance to water, presence of both species at a site, temperature, and rainfall. Argentine ants colonized unoccupied sites from neighboring sites, but the probability of appearance and persistence decreased as distance to development, vegetation cover, and elevation increased. Winter ants appeared and persisted in sites with relatively high vegetation cover (i.e., highly shaded sites). Presence of the 2 species was negatively associated in sites with high vegetation cover (more winter ants) and sites near development (more Argentine ants). Probability of colonization of Argentine ants decreased where winter ants were most persistent. At sites near development within the preserve, abundant Argentine ant populations may be excluding winter ants. The high abundance of Argentine ants at these sites may be due to immigration from suburban areas outside the preserve, which are high-quality habitat for Argentine ants. In the interior of the preserve, distance from development, low-quality habitat, and interaction with winter ants may in combination exclude Argentine ants. Interactions among the variables we examined were associated with low probabilities of Argentine ant colonization in the preserve.