Identifying habitat patches and potential ecological corridors for remnant Asiatic black bear (Ursus thibetanus japonicus) populations in Japan

The Japanese National Biodiversity Strategy 2010 calls for the creation of ecological networks as a biodiversity conservation policy. However, there is an obvious lack of information on the spatial distribution of many species and a lack of scientific methods for examining habitat requirements to establish the need for constructing these networks for target species. This study presents a quantitative method for assessing the need for ecological networks through modeling the potential geographic distributions of species based on a case study of local populations of Asiatic black bear (Ursus thibetanus japonicus) in Fuji and Tanzawa, Japan. A total of 1541 point records of occurrences of Asiatic black bears and 11 potential predictors were analyzed in a GIS environment. After a predictive distributional map was obtained using the Maximum Entropy (MaxEnt) algorithm, a gap analysis was carried out and population size was estimated. Approximately 24% of the bear's predicted habitat area fell within a wildlife protection area, 2% within a nature reserve, and 37% within natural parks. Conservation forest comprised 54% of the total area of predicted habitat; of this, national forest comprised 2%, and private and communal forest comprised 37%. The total estimated Asiatic black bear population in this region was 242, with 179 individuals in the Fuji local population, 26 in the Tanzawa local population, and 37 in the corridor patch between the two local populations. Our study also found a potential corridor connecting the Fuji and Tanzawa local populations, as well as potential habitat corridors in the Fuji region containing subpopulations on Mt. Fuji (119 individuals) and Mt. Kenashi (53 individuals). An additional subpopulation on Mt. Ashitaka (7 individuals) is isolated and not fully protected by a zoning plan. Mt. Furo's subpopulation is considered to be almost extinct, although black bears were observed here until 2002 based on the report by Mochizuki et al. (2005). The total black bear population of the Fuji-Tanzawa region is considered to be “endangered”; thus, an adequate population size might be difficult to maintain even if this region were to be internally connected by means of an ecological network.     

Characterizing Core and Corridor Use by Nubian Ibex in the Negev Desert, Israel

Abstract: Corridors have become a key element of conservation efforts. To identify and characterize corridors used by Nubian ibex ( Capra ibex nubiana ), we analyzed sighting data recorded for the past 20 years in the Israel Nature Reserves Authority data banks. We categorized each recorded sighting as belonging to a corridor or a core zone based on the total number of sightings in its vicinity. We identified three main core populations, a natural corridor connecting two of them, and a corridor that was not evident connecting the third population. Ibex inside and outside core zones were sighted on steeper terrain than expected by chance; this was more pronounced inside the core zones than outside them. We hypothesize that because ibex outside core zones must move rapidly and directionally across unfamiliar habitats, they must use more moderate terrain. Although ibex sightings in core zones were mostly in the vicinity of water, we found no relationship between the location of the sightings and proximity to water sources in the corridor. Hence, water does not appear to be an important factor in movement through corridors. It was more common to observe ibex out of the core zones during the summer. Males, which can be twice the size of females, were found traveling alone in corridors more often than in core zones. There was no difference between males and females in the steepness of terrain in which they were sighted outside core zones. Our data show that protecting ibex habitat in core zones and corridors is important to ibex conservation in Israel's arid zones. In addition, protecting this habitat may benefit other rock-dwelling species in the area.     

Spatial uncertainty analysis of population models

This paper describes an approach for conducting spatial uncertainty analysis of spatial population models, and illustrates the ecological consequences of spatial uncertainty for landscapes with different properties. Spatial population models typically simulate birth, death, and migration on an input map that describes habitat. Typically, only a single “reference” map is available, but we can imagine that a collection of other, slightly different, maps could be drawn to represent a particular species’ habitat. As a first approximation, our approach assumes that spatial uncertainty (i.e., the variation among values assigned to a location by such a collection of maps) is constrained by characteristics of the reference map, regardless of how the map was produced. Our approach produces lower levels of uncertainty than alternative methods used in landscape ecology because we condition our alternative landscapes on local properties of the reference map. Simulated spatial uncertainty was higher near the borders of patches. Consequently, average uncertainty was highest for reference maps with equal proportions of suitable and unsuitable habitat, and no spatial autocorrelation. We used two population viability models to evaluate the ecological consequences of spatial uncertainty for landscapes with different properties. Spatial uncertainty produced larger variation among predictions of a spatially explicit model than those of a spatially implicit model. Spatially explicit model predictions of final female population size varied most among landscapes with enough clustered habitat to allow persistence. In contrast, predictions of population growth rate varied most among landscapes with only enough clustered habitat to support a small population, i.e., near a spatially mediated extinction threshold. We conclude that spatial uncertainty has the greatest effect on persistence when the amount and arrangement of suitable habitat are such that habitat capacity is near the minimum required for persistence.     

Continuous versus binary representations of landscape heterogeneity in spatially-explicit models of mobile populations

How a landscape is represented is an important structural assumption in spatially-explicit simulation models. Simple models tend to specify just habitat and non-habitat (binary), while more complex models may use multiple levels or a continuum of habitat quality (continuous). How these different representations influence model projections is unclear. To assess the influence of landscape representation on population models, I developed a general, individual-based model with local dispersal and examined population persistence across binary and continuous landscapes varying in the amount and fragmentation of habitat. In binary and continuous landscapes habitat and non-habitat were assigned a unique mean suitability. In continuous landscapes, suitability of each individual site was then drawn from a normal distribution with fixed variance. Populations went extinct less often and abundances were higher in continuous landscapes. Production in habitat and non-habitat was higher in continuous landscapes, because the range of habitat suitability sampled by randomly dispersing individuals was higher than the overall mean habitat suitability. Increasing mortality, dispersal distance, and spatial heterogeneity all increased the discrepancy between continuous and binary landscapes. The effect of spatial structure on the probability of extinction was greater in binary landscapes. These results show that, under certain circumstances, model projections are influenced by how variation in suitability within a landscape is represented. Care should be taken to assess how a given species actually perceives the landscape when conducting population viability analyses or empirical validation of theory.     

A decision-support system in ICZM for protecting the ecosystems: integration with the habitat directive

This paper describes a decision-support system based on landscape ecology and focused on the study of ecosystems’ health. System capabilities are illustrated with three cases of integrated coastal zone management (ICZM), in the Adriatic Sea (Italy): the lagoon of Venice and the Rimini and Ancona coastal areas. Indicators and indices are developed with a focus on sub-regional and local problems in coastal management, with a multi-scale approach based on landscape and seascape ecology. Land-use changes of the coastal areas were detected by analyzing two sets of satellite images. Indices combining satellite imagery, socio-economic and environmental indicators, and landscape and seascape maps were created, showing ecological changes, habitat loss and gaps in conservation policy. The approach used provides means for the identification of conflicts and for the assessment of sustainability. Results show that the lagoon of Venice plays an important role in mitigating and compensating the impacts of human activities, and needs to be protected and restored. The Rimini area shows high ecological footprint and development-intensity and low biocapacity. The Ancona area needs the protection of its natural coastal space from potential sources of anthropogenic impacts to maintain its sustainability. A model of environment changes is critical for formulating effective environmental policies and management strategies. The developed decision-support system provides a suitability map per each area analyzed, which can be used in order to maximize different policy objectives and reduce coastal conflicts.     

Habitat Suitability Models and the Shortfall in Conservation Planning for African Vertebrates

Abstract:  Ongoing loss of biodiversity requires identifying large-scale conservation priorities, but the detailed information on the distribution of species required for this purpose is often missing. We present a systematic reserve selection for 1223 African mammals and amphibians in which habitat suitability models are used as estimates of the area occupied by species. In the framework of the World Conservation Union (IUCN) Global Amphibian Assessment and IUCN Global Mammal Assessment, we collected the geographic range (extent of occurrence) and habitat preferences for each species. We used the latter to build species-specific habitat suitability models inside geographic ranges, and for 181 species we verified the models by comparing suitability levels to presence-absence data collected in the field. We then used the suitable areas as estimators of the area of occupancy and compared the results of systematic reserve selection based on geographic ranges to those based on estimated areas of occupancy. Our results showed that the reserve system would need a 30-100% expansion to achieve minimal conservation targets, concentrated in the tropics, where species richness reaches a maximum. Comparative analyses revealed that using geographic ranges, which overestimate the area occupied by species, underestimates the total amount of area that needs to be conserved. The area selected for conservation doubled when we used the estimated area of occupancy in place of the geographic ranges. This happened because the suitable areas potentially occupied by each species overlapped less than their geographic ranges. As a result, any given protected area contained fewer species than predicted by the analysis of ranges. Because species are more specialized than our estimates of distribution based on extent of occurrence suggest, we propose that this is a general effect in systematic conservation planning.     

Evaluating mortality rates with a novel integrated framework for nonmonogamous species

The conservation of wildlife requires management based on quantitative evidence, and especially for large carnivores, unraveling cause‐specific mortalities and understanding their impact on population dynamics is crucial. Acquiring this knowledge is challenging because it is difficult to obtain robust long‐term data sets on endangered populations and, usually, data are collected through diverse sampling strategies. Integrated population models (IPMs) offer a way to integrate data generated through different processes. However, IPMs are female‐based models that cannot account for mate availability, and this feature limits their applicability to monogamous species only. We extended classical IPMs to a two‐sex framework that allows investigation of population dynamics and quantification of cause‐specific mortality rates in nonmonogamous species. We illustrated our approach by simultaneously modeling different types of data from a reintroduced, unhunted brown bear (Ursus arctos) population living in an area with a dense human population. In a population mainly driven by adult survival, we estimated that on average 11% of cubs and 61% of adults died from human‐related causes. Although the population is currently not at risk, adult survival and thus population dynamics are driven by anthropogenic mortality. Given the recent increase of human‐bear conflicts in the area, removal of individuals for management purposes and through poaching may increase, reversing the positive population growth rate. Our approach can be generalized to other species affected by cause‐specific mortality and will be useful to inform conservation decisions for other nonmonogamous species, such as most large carnivores, for which data are scarce and diverse and thus data integration is highly desirable.     

扎龙湿地丹顶鹤巢址空间分布变化及其对环境变化指征

扎龙国家级湿地保护区是我国一级保护鸟类丹顶鹤最重要的集中繁殖栖息地。以扎龙保护区为研究区,利用质心模型和最近邻体模型等方法,模拟1996至2004年间丹顶鹤（Grus japonensis）巢址的空间分布变化,分析这种变化对环境变化的指示作用。结果表明：丹顶鹤种群巢址分布出现向核心区集中的现象,鹤巢平均最近邻体距离大幅减少,巢址的空间分布格局由1996年的均匀分布转变为2004年的成群分布型;丹顶鹤巢址的空间分布变化表明该区环境发生巨大的变化：1996年核心区内栖息资源相对均匀、丰富,而缓冲区和实验区适宜生境呈零星分布格局,此后生境质量经历了由适宜一非适宜,适宜生境空间分布由均匀一斑块化的过程,到2004年各功能区适宜生境数量大幅减少,且核心区内适宜生境斑块化严重。为研究区丹顶鹤种群及其生境的保护和管理提供参考。     

Do resources or natural enemies drive bee population dynamics in fragmented habitats?

The relative importance of bottom-up or top-down forces has been mainly studied for herbivores but rarely for pollinators. Habitat fragmentation might change driving forces of population dynamics by reducing the area of resource-providing habitats, disrupting habitat connectivity, and affecting natural enemies more than their host species. We studied spatial and temporal population dynamics of the solitary bee Osmia rufa (Hymenoptera: Megachilidae) in 30 fragmented orchard meadows ranging in size from 0.08 to 5.8 ha in an agricultural landscape in central Germany. From 1998 to 2003, we monitored local bee population size, rate of parasitism, and rate of larval and pupal mortality in reed trap nests as an accessible and standardized nesting resource. Experimentally enhanced nest site availability resulted in a steady increase of mean local population size from 80 to 2740 brood cells between 1998 and 2002. Population size and species richness of natural enemies increased with habitat area, whereas rate of parasitism and mortality only varied among years. Inverse density-dependent parasitism in three study years with highest population size suggests rather destabilizing instead of regulating effects of top-down forces. Accordingly, an analysis of independent time series showed on average a negative impact of population size on population growth rates but provides no support for top-down regulation by natural enemies. We conclude that population dynamics of O. rufa are mainly driven by bottom-up forces, primarily nest site availability.     

The Effect of Fragment Shape and Species' Sensitivity to Habitat Edges on Animal Population Size

Abstract:  Habitat fragmentation causes extinction of local animal populations by decreasing the amount of viable "core" habitat area and increasing edge effects. It is widely accepted that larger fragments make better nature reserves because core-dwelling species have a larger amount of suitable habitat. Nevertheless, fragments in real landscapes have complex, irregular shapes. We modeled the population sizes of species that have a representative range of preferences for or aversions to habitat edges at five spatial scales (within 10, 32, 100, 320, and 1000 m of an edge) in a nation-wide analysis of forest remnants in New Zealand. We hypothesized that the irregular shapes of fragments in real landscapes should generate statistically significant correlations between population density and fragment area, purely as a "geometric" effect of varying species responses to the distribution of edge habitat. Irregularly shaped fragments consistently reduced the population size of core-dwelling species by 10–100%, depending on the scale over which species responded to habitat edges. Moreover, core populations within individual fragments were spatially discontinuous, containing multiple, disjunct populations that inhabited small spatial areas and had reduced population size. The geometric effect was highly nonlinear and depended on the range of fragment sizes sampled and the scale at which species responded to habitat edges. Fragment shape played a strong role in determining population size in fragmented landscapes; thus, habitat restoration efforts may be more effective if they focus on connecting disjunct cores rather than isolated fragments.     

Contrasting ecological information content in whaling archives with modern cetacean surveys for conservation planning and identification of historical distribution changes

Many species are restricted to a marginal or suboptimal fraction of their historical range due to anthropogenic impacts, making it hard to interpret their ecological preferences from modern-day data alone. However, inferring past ecological states is limited by the availability of robust data and biases in historical archives, posing a challenge for policy makers . To highlight how historical records can be used to understand the ecological requirements of threatened species and inform conservation, we investigated sperm whale (Physeter macrocephalus) distribution in the Western Indian Ocean. We assessed differences in information content and habitat suitability predictions based on whale occurrence data from Yankee whaling logs (1792–1912) and from modern cetacean surveys (1995–2020). We built maximum entropy habitat suitability models containing static (bathymetry-derived) variables to compare models comprising historical-only and modern-only data. Using both historical and modern habitat suitability predictions  we assessed marine protected area (MPA) placement by contrasting suitability in- and outside MPAs. The historical model predicted high habitat suitability in shelf and coastal regions near continents and islands, whereas the modern model predicted a less coastal distribution with high habitat suitability more restricted to areas of steep topography. The proportion of high habitat suitability inside versus outside MPAs was higher when applying the historical predictions than the modern predictions, suggesting that different marine spatial planning optimums can be reached from either data sources. Moreover, differences in relative habitat suitability predictions between eras were consistent with the historical depletion of sperm whales from coastal regions, which were easily accessed and targeted by whalers, resulting in a modern distribution limited more to steep continental margins and remote oceanic ridges. The use of historical data can provide important new insights and, through cautious interpretation, inform conservation planning and policy, for example, by identifying refugee species and regions of anticipated population recovery.     

A Habitat-Based Metapopulation Model of the California Gnatcatcher

We present an analysis of the metapopulation dynamics of the federally threatened coastal California Gnatcatcher (Polioptila c. californica) for an approximately 850 km² region of Orange County, California. We developed and validated a habitat suitability model for this species using data on topography, vegetation, and locations of gnatcatcher pair observations. Using this habitat model, we calculated the spatial structure of the metapopulation, including size and location of habitat patches and the distances among them. We used data based on field studies to estimate parameters such as survival, fecundity, dispersal, and catastrophes, and combined these parameters with the spatial structure to build a stage-structured, stochastic, spatially-explicit metapopulation model. The model predicted a fast decline and high risk of population extinction with most combinations of parameters. Results were most sensitive to density-dependent effects, the probability of weather-related catastrophes, adult survival, and adult fecundity. Based on data used in the model, the greatest difference in results was given when the simulation's time horizon was only a few decades, suggesting that modeling based on longer or shorter time horizons may underestimate the effects of alternative management actions.     

基于InVEST模型的三峡库区(重庆段)生境质量时空演变分析

为识别生境质量优劣的具体区域,揭示土地利用变化下生境质量时空演变特征,对集大农村、大库区于一体的三峡库区(重庆段)开展生境质量研究.以2000、2010及2020年土地利用数据为基础,运用InVEST模型对三峡库区(重庆段)2000—2020年生境质量进行定量评价,并结合转移矩阵方法挖掘了生境质量退化的具体区域.结果表...     

Identifying conservation priorities for threatened Eastern Himalayan mammals

To augment mammal conservation in the Eastern Himalayan region, we assessed the resident 255 terrestrial mammal species and identified the 50 most threatened species based on conservation status, endemism, range size, and evolutionary distinctiveness. By using the spatial analysis package letsR and the complementarity core‐area method in the conservation planning software Zonation, we assessed the current efficacy of their protection and identified priority conservation areas by comparing protected areas (PAs), land cover, and global ecoregion 2017 maps at a 100 × 100 m spatial scale. The 50 species that were most threatened, geographically restricted, and evolutionarily distinct faced a greater extinction risk than globally nonthreatened and wide‐ranging species and species with several close relatives. Small, medium‐sized, and data‐deficient species faced extinction from inadequate protection in PAs relative to wide‐ranging charismatic species. There was a mismatch between current PA distribution and priority areas for conservation of the 50 most endangered species. To protect these species, the skewed regional PA distribution would require expansion. Where possible, new PAs and transboundary reserves in the 35 priority areas we identified should be established. There are adequate remaining natural areas in which to expand current Eastern Himalayan PAs. Consolidation and expansion of PAs in the EH requires strengthening national and regional transboundary collaboration, formulating comprehensive regional land‐use plans, diversifying conservation funding, and enhancing information sharing through a consolidated regional database.     

A Voronoi diagram based population model for social species of wildlife

A population model is presented that accounts for spatial structure within habitat patches. It is designed for social species of wildlife that form social group home ranges that are much smaller than patch size. The model represents social group home ranges by Voronoi regions that tessellate a patch to form a Voronoi diagram. Neighbouring social groups are linked with habitat-confined shortest paths and form a dispersal network. The model simulates population dynamics and makes use of Voronoi diagrams and dispersal networks as a spatial component. It then produces density maps as outputs. These are maps that show predicted animal densities across the patches of a landscape. A construction procedure for the particular Voronoi diagram type used by the model is described. As a test case, the model is run for the squirrel glider (Petaurus norfolcensis), a small arboreal marsupial native to Australia. A time series of density maps are produced that show squirrel glider density changing across a landscape through time.     

A dispersal-constrained habitat suitability model for predicting invasion of alpine vegetation

Developing tools to predict the location of new biological invasions is essential if exotic species are to be controlled before they become widespread. Currently, alpine areas in Australia are largely free of exotic plant species but face increasing pressure from invasive species due to global warming and intensified human use. To predict the potential spread of highly invasive orange hawkweed (Hieracium aurantiacum) from existing founder populations on the Bogong High Plains in southern Australia, we developed an expert-based, spatially explicit, dispersal-constrained, habitat suitability model. The model combines a habitat suitability index, developed from disturbance, site wetness, and vegetation community parameters, with a phenomenological dispersal kernel that uses wind direction and observed dispersal distances. After generating risk maps that defined the relative suitability of H. aurantiacum establishment across the study area, we intensively searched several locations to evaluate the model. The highest relative suitability for H. aurantiacum establishment was southeast from the initial infestations. Native tussock grasslands and disturbed areas had high suitability for H. aurantiacum establishment. Extensive field searches failed to detect new populations. Time-step evaluation using the location of populations known in 1998-2000, accurately assigned high relative suitability for locations where H. aurantiacum had established post-2003 (AUC [area under curve] = 0.855 +/- 0.035). This suggests our model has good predictive power and will improve the ability to detect populations and prioritize areas for ongoing monitoring.     

The importance of incorporating functional habitats into conservation planning for highly mobile species in dynamic systems

The distribution of mobile species in dynamic systems can vary greatly over time and space. Estimating their population size and geographic range can be problematic and affect the accuracy of conservation assessments. Scarce data on mobile species and the resources they need can also limit the type of analytical approaches available to derive such estimates. We quantified change in availability and use of key ecological resources required for breeding for a critically endangered nomadic habitat specialist, the Swift Parrot (Lathamus discolor). We compared estimates of occupied habitat derived from dynamic presence‐background (i.e., presence‐only data) climatic models with estimates derived from dynamic occupancy models that included a direct measure of food availability. We then compared estimates that incorporate fine‐resolution spatial data on the availability of key ecological resources (i.e., functional habitats) with more common approaches that focus on broader climatic suitability or vegetation cover (due to the absence of fine‐resolution data). The occupancy models produced significantly (P < 0.001) smaller (up to an order of magnitude) and more spatially discrete estimates of the total occupied area than climate‐based models. The spatial location and extent of the total area occupied with the occupancy models was highly variable between years (131 and 1498 km²). Estimates accounting for the area of functional habitats were significantly smaller (2–58% [SD 16]) than estimates based only on the total area occupied. An increase or decrease in the area of one functional habitat (foraging or nesting) did not necessarily correspond to an increase or decrease in the other. Thus, an increase in the extent of occupied area may not equate to improved habitat quality or function. We argue these patterns are typical for mobile resource specialists but often go unnoticed because of limited data over relevant spatial and temporal scales and lack of spatial data on the availability of key resources. Understanding changes in the relative availability of functional habitats is crucial to informing conservation planning and accurately assessing extinction risk for mobile resource specialists.     

Demographic Analyses of a Hunted Black Bear Population with Access to a Refuge

Studying one of two bear species not experiencing widespread population decline, provides insight into the population responses of the six bear species that are in decline and into responses of other long-lived species for which data are difficult to collect. Black bear ( Ursus americanus ) sanctuaries were established in North Carolina (U.S.) in 1971 to protect core populations of bears and to provide dispersing bears for hunting. Population index values, derived from counts of bears visiting bait stations, were significantly greater inside the Pisgah Bear Sanctuary than outside and were greater along trails than along roads. Survivorship of bears outfitted with transmitter collars was greater for sanctuary bears alone than for sanctuary plus non-sanctuary bears. Monte Carlo analyses of Leslie matrices showed that the bear population in the sanctuary would be stable if cub survivorship, p ₀, was about O.7, and the population in the sanctuary plus the surrounding area would be stable if p ₀ was about 0.83. Estimates of litter survivorship in North Carolina indicate, however, that p ₀ can not exceed O.71. Overall, the matrix analyses indicated an ultimate population decline in the total bear population (sanctuary plus surrounding area). The population index of the bait station did not show a discernible decline. The Pisgah Bear Sanctuary provides dispersing bears for hunters and provides some protection for the resident bears. The sanctuary may not, however, provide resident bears with enough protection to maintain a viable breeding population within its boundaries. Reducing human access to bears and their habitat appears crucial, either by making large sanctuaries or by eliminating roads.     

Application of least-cost path model to identify a giant panda dispersal corridor network after the Wenchuan earthquake—Case study of Wolong Nature Reserve in China

With growing levels of human-activity and frequent natural disturbances throughout the world, it is increasingly important that both research and management efforts take into account the widespread landscape fragmentation and its consequences for biodiversity conservation. The magnitude 5.12 Wenchuan earthquake in China caused dramatic impacts on giant panda (Ailuropoda melanoleuca) habitat in the nature reserves within Minshan and Qionglai mountains. With the combined stresses of the natural disaster and the extensive human activities during postquake reconstruction, giant panda habitat in this region may become more fragmented in the future. In order to preserve the giant panda population after the earthquake and protect the species against habitat fragmentation, this article explores a method of identifying giant panda migration corridors involving habitat suitability assessments and a least-cost path model. Focusing on postquake Wolong Nature Reserve, our results demonstrate that it contains 430.3 km² of suitable habitat (21.1% of total area), 463.8 km² of marginally suitable habitat (22.8%) and 1141.9 km² of unsuitable habitat (50.1%). We further show that several giant panda dispersal corridors exist in the reserve, including four corridor groups that cross the provincial highway and five corridors that do not intersect areas of human activity. This study will contribute to management and conservation efforts in Wolong Nature Reserve and beyond after the Wenchuan earthquake.     

Modeling range dynamics in heterogeneous landscapes: invasion of the hemlock woolly adelgid in eastern North America

Range expansion by native and exotic species will continue to be a major component of global change. Anticipating the potential effects of changes in species distributions requires models capable of forecasting population spread across realistic, heterogeneous landscapes and subject to spatiotemporal variability in habitat suitability. Several decades of theory and model development, as well as increased computing power and availability of fine-resolution GIS data, now make such models possible. Still unanswered, however, is the question of how well this new generation of dynamic models will anticipate range expansion. Here we develop a spatially explicit stochastic model that combines dynamic dispersal and population processes with fine-resolution maps characterizing spatiotemporal heterogeneity in climate and habitat to model range expansion of the hemlock woolly adelgid (HWA; Adelges tsugae). We parameterize this model using multiyear data sets describing population and dispersal dynamics of HWA and apply it to eastern North America over a 57-year period (1951-2008). To evaluate the model, the observed pattern of spread of HWA during this same period was compared to model predictions. Our model predicts considerable heterogeneity in the risk of HWA invasion across space and through time, and it suggests that spatiotemporal variation in winter temperature, rather than hemlock abundance, exerts a primary control on the spread of HWA. Although the simulations generally matched the observed current extent of the invasion of HWA and patterns of anisotropic spread, it did not correctly predict when HWA was observed to arrive in different geographic regions. We attribute differences between the modeled and observed dynamics to an inability to capture the timing and direction of long-distance dispersal events that substantially affected the ensuing pattern of spread.