The Application of Neutral Landscape Models in Conservation Biology

Neutral landscape models, derived from percolation theory in the field of landscape ecology, are grid-based maps in which complex habitat distributions are generated by random or fractal algorithms. This grid-based representation of landscape structure is compatible with the raster-based format of geographical information systems (GIS), which facilitates comparisons between theoretical and real landscapes. Neutral landscape models permit the identification of critical thresholds in connectivity, which can be used to predict when landscapes will become fragmented. The coupling of neutral landscape models with generalized population models, such as metapopulation theory, provides a null model for generating predictions about population dynamics in fragmented landscapes. Neutral landscape models can contribute to the following applications in conservation: (1) incorporation of complex spatial patterns in (meta)population models; (2) identification of species' perceptions of landscape structure; (3) determination of landscape connectivity; (4) evaluation of the consequences of habitat fragmentation for population subdivision; (5) identification of the domain of metapopulation dynamics; (6) prediction of the occurrence of extinction thresholds; ( 7) determination of the genetic consequences of habitat fragmentation; and (8) reserve design and ecosystem management. This generalized, spatially explicit framework bridges the gap between spatially implicit, patch-based models and spatially realistic GIS applications which are usually parameterized for a single species in a specific landscape. Development of a generalized, spatially explicit framework is essential in conservation biology because we will not be able to develop individual models for every species of management concern.     

A geostatistical approach for mapping thematic classification accuracy and evaluating the impact of inaccurate spatial data on ecological model predictions

Spatial information in the form of geographical information system coverages and remotely sensed imagery is increasingly used in ecological modeling. Examples include maps of land cover type from which ecologically relevant properties, such as biomass or leaf area index, are derived. Spatial information, however, is not error-free: acquisition and processing errors, as well as the complexity of the physical processes involved, make remotely sensed data imperfect measurements of ecological attributes. It is therefore important to first assess the accuracy of the spatial information being used and then evaluate the impact of such inaccurate information on ecological model predictions. In this paper, the role of geostatistics for mapping thematic classification accuracy through integration of abundant image-derived (soft) and sparse higher accuracy (hard) class labels is presented. Such assessment leads to local indices of map quality, which can be used for guiding additional ground surveys. Stochastic simulation is proposed for generating multiple alternative realizations (maps) of the spatial distribution of the higher accuracy class labels over the study area. All simulated realizations are consistent with the available pieces of information (hard and soft labels) up to their validated level of accuracy. The simulated alternative class label representations can be used for assessing joint spatial accuracy, i.e., classification accuracy regarding entire spatial features read from the thematic map. Such realizations can also serve as input parameters to spatially explicit ecological models; the resulting distribution of ecological responses provides a model of uncertainty regarding the ecological model prediction. A case study illustrates the generation of alternative land cover maps for a Landsat Thematic Mapper (TM) subscene, and the subsequent construction of local map quality indices. Simulated land cover maps are then input into a biogeochemical model for assessing uncertainty regarding net primary production (NPP).     

Dispersal in Spatially Explicit Population Models

Abstract: Ruckelshaus et al. (1997) outlined a simulation model of dispersal between patches in a fragmented landscape. They showed that dispersal success—the proportion of dispersers successfully locating a patch—was particularly sensitive to errors in dispersal mortality and concluded that this limits the utility of spatially explicit population models in conservation biology. I contend that, although they explored error propagation in a simple dispersal model, they did not explore how errors are propagated in spatially explicit population models, as no consideration of population processes was included. I developed a simple simulation model to investigate the effect of varying dispersal success on predictions of patch occupancy and population viability, the conventional outputs of spatially explicit population models. The model simulates births and deaths within habitat patches and dispersal as the transfer of individuals between them. Model predictions were sensitive to changes in dispersal success across a restricted range of within-patch growth rates, which depended on the dispersal initiation mechanism, patch carrying capacities, and number of generations simulated. Predictions of persistence and patch occupancy were generally more sensitive to changes in dispersal success (1) under presaturation rather than saturation dispersal; (2) at lower patch carrying capacities; and (3) over longer time periods. The framework I present provides a means of assessing, quantitatively, the regions of parameter space for which differences in dispersal success are likely to have a large effect on population model outputs. Investigating the effect of the representation of dispersal behavior within the demographic and landscape context provides a more useful assessment of whether our lack of knowledge is likely to cause unacceptable uncertainty in the predictions of spatially explicit population models.     

Graph theory as a proxy for spatially explicit population models in conservation planning.

Spatially explicit population models (SEPMs) are often considered the best way to predict and manage species distributions in spatially heterogeneous landscapes. However, they are computationally intensive and require extensive knowledge of species' biology and behavior, limiting their application in many cases. An alternative to SEPMs is graph theory, which has minimal data requirements and efficient algorithms. Although only recently introduced to landscape ecology, graph theory is well suited to ecological applications concerned with connectivity or movement. This paper compares the performance of graph theory to a SEPM in selecting important habitat patches for Wood Thrush (Hylocichla mustelina) conservation. We use both models to identify habitat patches that act as population sources and persistent patches and also use graph theory to identify patches that act as stepping stones for dispersal. Correlations of patch rankings were very high between the two models. In addition, graph theory offers the ability to identify patches that are very important to habitat connectivity and thus long-term population persistence across the landscape. We show that graph theory makes very similar predictions in most cases and in other cases offers insight not available from the SEPM, and we conclude that graph theory is a suitable and possibly preferable alternative to SEPMs for species conservation in heterogeneous landscapes.     

Local population size in a flightless insect: importance of patch structure-dependent mortality

In spatially heterogeneous systems, utilizing population models to integrate the effects of multiple population rates can yield powerful insights into the relative importance of the component rates. The relative importance of demographic rates and dispersal in shaping the distribution of the western tussock moth (Orgyia vetusta) among patches of its host plant was explored using stage-structured population models. Tussock moth dispersal occurs passively in first-instar larvae and is poor or absent in all other life stages. Spatial surveys suggested, however, that moth distribution is not well explained by passive dispersal; moth populations were greater on small patches and on isolated ones. Further analysis showed that several local demographic rates varied significantly with patch characteristics. Two mortality factors in particular may explain the observed patterns. First, crawler mortality both increased with patch size and was density-dependent. A single-patch difference equation model showed mortality related to patch size is strong enough to overcome the homogenizing effect of density dependence; greater equilibrium densities were predicted for smaller patches. Second, although three rates were found to vary with local patch density, only pupal parasitism by a chalcid wasp could potentially account for higher moth abundances on isolated patches. A spatially explicit simulation model of the multiple-patch system showed that spatial variation in pupal parasitism is indeed strong enough to generate such a pattern. These results demonstrate that habitat spatial structure can affect multiple population processes simultaneously, and even relatively low attack rates imposed on a reproductively valuable life stage of the host can have a dominant effect on population distribution among habitat patches.     

Modeling Species' Distributions to Improve Conservation in Semiurban Landscapes: Koala Case Study

Abstract:  Models of species' distributions are commonly used to inform landscape and conservation planning. In urban and semiurban landscapes, the distributions of species are determined by a combination of natural habitat and anthropogenic impacts. Understanding the spatial influence of these two processes is crucial for making spatially explicit decisions about conservation actions. We present a logistic regression model for the distribution of koalas (  Phascolarctos cinereus ) in a semiurban landscape in eastern Australia that explicitly separates the effect of natural habitat quality and anthropogenic impacts on koala distributions. We achieved this by comparing the predicted distributions from the model with what the predicted distributions would have been if anthropogenic variables were at their mean values. Similar approaches have relied on making predictions assuming anthropogenic variables are zero, which will be unreliable if the training data set does not include anthropogenic variables close to zero. Our approach is novel because it can be applied to landscapes where anthropogenic variables are never close to zero. Our model showed that, averaged across the study area, natural habitat was the main determinant of koala presence. At a local scale, however, anthropogenic impacts could be more important, with consequent implications for conservation planning. We demonstrated that this modeling approach, combined with the visual presentation of predictions as a map, provides important information for making decisions on how different conservation actions should be spatially allocated. This method is particularly useful for areas where wildlife and human populations exist in close proximity.     

Aphid population response to agricultural landscape change: A spatially explicit,individual-based model

A spatially explicit individual-based simulation model has been developed to represent aphid population dynamics in agricultural landscapes. The application of the model to Rhopalosiphum padi (L.) population dynamics is detailed, including an outline of the construction of the model, its parameterisation and validation. Over time, the aphids interact with the landscape and with one another. The landscape is modified by varying a simple pesticide regime, and the multi-scale spatial and temporal implications for a population of aphids is analysed. The results show that a spatial modelling approach that considers the effects on the individual of landscape properties and factors such as wind speed and wind direction provides novel insight into aphid population dynamics both spatially and temporally. This forms the basis for the development of further simulation models that can be used to analyse how changes in landscape structure impact upon important species distributions and population dynamics.     

Mapping legal authority for terrestrial conservation corridors along streams

Wildlife corridors aim to promote species’ persistence by connecting habitat patches across fragmented landscapes. Their implementation is limited by patterns of land ownership and complicated by differences in the jurisdictional and regulatory authorities under which lands are managed. Terrestrial corridor conservation requires coordination across jurisdictions and sectors subject to site-specific overlapping sources of legal authority. Mapping spatial patterns of legal authority concurrent with habitat condition can illustrate opportunities to build or leverage capacity for connectivity conservation. Streamside areas provide pragmatic opportunities to leverage existing policy mechanisms for riverine and terrestrial habitat connectivity across boundaries. Conservation planners and practitioners can make use of these opportunities by harmonizing actions for multiple conservation outcomes. We formulated an integrative, data-driven method for mapping multiple sources of legal authority weighted by capacity for coordinating terrestrial habitat conservation along streams. We generated a map of capacity to coordinate streamside corridor protections across a wildlife habitat gap to demonstrate this approach. We combined values representing coordination capacity and naturalness to generate an integrated legal-ecological resistance map for connectivity modeling. We then computed least-cost corridors across the integrated map, masking the terrestrial landscape to focus on streamside areas. Streamside least-cost corridors in the integrated, local-scale model diverged (∼25 km) from national-scale least-cost corridors based on naturalness. Spatial categories comparing legal- and naturalness-based resistance values by stream reach highlighted potential locations for building or leveraging existing capacity through spatial coordination of policy mechanisms or restoration actions. Agencies or nongovernmental organizations intending to restore or maintain habitat connectivity across fragmented landscapes can use this approach to inform spatial prioritization and build coordination capacity. Article impact statement: Combined mapping of legal authority and habitat condition reveals capacity to coordinate actions along streams for clean water and wildlife.     

Integrating functional connectivity and fire management for better conservation outcomes

Globally, the mean abundance of terrestrial animals has fallen by 50% since 1970, and populations face ongoing threats associated with habitat loss, fragmentation, climate change, and disturbance. Climate change can influence the quality of remaining habitat directly and indirectly by precipitating increases in the extent, frequency, and severity of natural disturbances, such as fire. Species face the combined threats of habitat clearance, changing climates, and altered disturbance regimes, each of which may interact and have cascading impacts on animal populations. Typically, conservation agencies are limited in their capacity to mitigate rates of habitat clearance, habitat fragmentation, or climate change, yet fire management is increasingly used worldwide to reduce wildfire risk and achieve conservation outcomes. A popular approach to ecological fire management involves the creation of fire mosaics to promote animal diversity. However, this strategy has 2 fundamental limitations: the effect of fire on animal movement within or among habitat patches is not considered and the implications of the current fire regime for long-term population persistence are overlooked. Spatial and temporal patterns in fire history can influence animal movement, which is essential to the survival of individual animals, maintenance of genetic diversity, and persistence of populations, species, and ecosystems. We argue that there is rich potential for fire managers to manipulate animal movement patterns; enhance functional connectivity, gene flow, and genetic diversity; and increase the capacity of populations to persist under shifting environmental conditions. Recent methodological advances, such as spatiotemporal connectivity modeling, spatially explicit individual-based simulation, and fire-regime modeling can be integrated to achieve better outcomes for biodiversity in human-modified, fire-prone landscapes. Article impact statement: Land managers may conserve populations by using fire to sustain or enhance functional connectivity.     

Metapopulation Extinction Risk under Spatially Autocorrelated Disturbance

Abstract:  Recent extinction models generally show that spatial aggregation of habitat reduces overall extinction risk because sites emptied by local extinction are more rapidly recolonized. We extended such an investigation to include spatial structure in the disturbance regime. A spatially explicit metapopulation model was developed with a wide range of dispersal distances. The degree of aggregation of both habitat and disturbance pattern could be varied from a random distribution, through the intermediate case of a fractal distribution, all the way to complete aggregation (single block). Increasing spatial aggregation of disturbance generally increased extinction risk. The relative risk faced by populations in different landscapes varied greatly, depending on the disturbance regime. With random disturbance, the spatial aggregation of habitat reduced extinction risk, as in earlier studies. Where disturbance was spatially autocorrelated, however, this advantage was eliminated or reversed because populations in aggregated habitats are at risk of mass extinction from coarse-scale disturbance events. The effects of spatial patterns on extinction risk tended to be reduced by long-distance dispersal. Given the high levels of spatial correlation in natural and anthropogenic disturbance processes, population vulnerability may be greatly underestimated both by classical (nonspatial) models and by those that consider spatial structure in habitat alone.     

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.     

Unveiling human-assisted dispersal mechanisms in invasive alien insects: Integration of spatial stochastic simulation and phenology models

Capturing the spread of biological invasions in heterogeneous landscapes is a complex modelling task where information on both dispersal and population dynamics needs to be integrated. Spatial stochastic simulation and phenology models have rarely been combined to assist in the study of human-assisted long-distance dispersal events.Here we develop a process-based spatially explicit landscape-extent simulation model that considers the spread and detection of invasive insects. Natural and human-assisted dispersal mechanisms are modelled with an individual-based approach using negative exponential and negative power law dispersal kernels and gravity models. The model incorporates a phenology sub-model that uses daily temperature grids for the prediction and timing of the population dynamics in each habitat patch. The model was applied to the study of the invasion by the important maize pest western corn rootworm (WCR) Diabrotica virgifera ssp. virgifera in Europe. We parameterized and validated the model using maximum likelihood and simulation methods from the historical invasion of WCR in Austria.WCR was found to follow stratified dispersal where international transport networks in the Danube basin played a key role in the occurrence of long-distance dispersal events. Detection measures were found to be effective and altitude had a significant effect on limiting the spread of WCR. Spatial stochastic simulation combined with phenology models, maximum likelihood methods and predicted versus observed regression showed a high degree of flexibility that captured the salient features of WCR spread in Austria. This modelling approach is useful because it allows to fully exploit and the often limited and heterogeneous information available regarding the population dynamics and dispersal of alien invasive insects.     

Crustacean habitat potential mapping in a tidal flat using remote sensing and GIS

To make a macrofaunal (crustacean) habitat potential map, the spatial distribution of ecological variables in the Hwangdo tidal flat, Korea, was explored. Spatial variables were mapped using remote sensing and a geographic information system (GIS) combined with field observations. A frequency ratio (FR) and logistic regression (LR) model were employed to map the macrofauna potential area for the Ilyoplax dentimerosa, a crustacean species. Spatial variables affecting the tidal macrofauna distribution were selected based on abundance and biomass and used within a spatial database derived from remotely sensed data of various types of sensors. The spatial variables included the intertidal digital elevation model (DEM), slope, distance from a tidal channel, tidal channel density, surface sediment facies, spectral reflectance of the near infrared (NIR) bands and the tidal exposure duration. The relation between the I. dentimerosa and each spatial variable was calculated using the FR and LR. The species was randomly divided into a training set (70%) to analyse habitat potential using FR and LR and a test set (30%) to validate the predicted habitat potential map. The relations were overlaid to produce a habitat potential map with the species potential index (SPI) value for each pixel. The potential habitat maps were compared with the surveyed habitat locations such as validation data set. The comparison results showed that the LR model (accuracy is 85.28%) is better in prediction than the FR (accuracy is 78.96%) model. The performance of models gave satisfactory accuracies. The LR provides the quantitative influence of variables on a potential habitat of species; otherwise, the FR shows the quantitative influence of a class in each variable. The combination of a GIS-based frequency ratio and logistic regression models and remote sensing with field observations is an effective method to determine locations favorable for macrofaunal species occurrences in a tidal flat.     

Spatial selection and inheritance: applying evolutionary concepts to population dynamics in heterogeneous space

Organisms in highly suitable sites generally produce more offspring, and offspring can inherit this suitability by not dispersing far. This combination of spatial selection and spatial inheritance acts to bias the distribution of organisms toward suitable sites and thereby increase mean fitness (i.e., per capita population increase). Thus, population growth rates in heterogeneous space change over time by a process conceptually analogous to evolution by natural selection, opening avenues for theoretical cross-pollination between evolutionary biology and ecology. We operationally define spatial inheritance and spatial selective differential and then combine these two factors in a modification of the breeder's equation, derived from simple models of population growth in heterogeneous space. The modified breeder's equation yields a conservative criterion for persistence in hostile environments estimable from field measurements. We apply this framework for understanding gypsy moth population persistence amidst abundant predators and find that the predictions of the modified breeder's equation match initial changes in population growth rate in independent simulation output. The analogy between spatial dynamics and natural selection conceptually links ecology and evolution, provides a spatially implicit framework for modeling spatial population dynamics, and represents an important null model for studying habitat selection.     

Comparing spatial conservation prioritization methods with site- versus spatial dependency-based connectivity

Larval dispersal is an important component of marine reserve networks. Two conceptually different approaches to incorporate dispersal connectivity into spatial planning of these networks exist, and it is an open question as to when either is most appropriate. Candidate reserve sites can be selected individually based on local properties of connectivity or on a spatial dependency-based approach of selecting clusters of strongly connected habitat patches. The first acts on individual sites, whereas the second acts on linked pairs of sites. We used a combination of larval dispersal simulations representing different seascapes and case studies of biophysical larval dispersal models in the Coral Triangle region and the province of Southeast Sulawesi, Indonesia, to compare the performance of these 2 methods in the spatial planning software Marxan. We explored the reserve design performance implications of different dispersal distances and patterns based on the equilibrium settlement of larvae in protected and unprotected areas. We further assessed different assumptions about metapopulation contributions from unprotected areas, including the case of 100% depletion and more moderate scenarios. The spatial dependency method was suitable when dispersal was limited, a high proportion of the area of interest was substantially degraded, or the target amount of habitat protected was low. Conversely, when subpopulations were well connected, the 100% depletion was relaxed, or more habitat was protected, protecting individual sites with high scores in metrics of connectivity was a better strategy. Spatial dependency methods generally produced more spatially clustered solutions with more benefits inside than outside reserves compared with site-based methods. Therefore, spatial dependency methods potentially provide better results for ecological persistence objectives over enhancing fisheries objectives, and vice versa. Different spatial prioritization methods of using connectivity are appropriate for different contexts, depending on dispersal characteristics, unprotected area contributions, habitat protection targets, and specific management objectives. Comparación entre los métodos de priorización de la conservación espacial con sitio y la conectividad espacial basada en la dependencia     

The spatial links tool: Automated mapping of habitat linkages in variegated landscapes

The removal, alteration and fragmentation of habitat in many parts of the world has led to a loss of biodiversity. Within the prevailing societal limitations the process is not easily reversed. Attempts are being made to minimise the fragmentation of remaining habitat by strategically reversing or managing habitat loss. Although their relative usefulness is a topic of debate among ecologists, habitat corridors are seen as one way of maintaining spatially dependent ecological processes within landscapes where habitat has been seriously depleted. Corridors can only be effective if they significantly contribute to the species sustaining processes of gene flow, resource access or the colonisation of vacant patches. We present a spatial habitat modelling methodology for evaluating the contribution and potential contribution of connecting paths to landscape connectivity. We have developed the spatial links tool (SLT), which maps link value across a region. The SLT combines connectivity measures from metapopulation ecology with the least cost path algorithm from graph theory, and can be applied to continuously variable landscape data. Combined with expert judgement, link value maps can be used to delineate habitat corridors. The approach capitalises on some synergies between ecological relevance and computational efficiency to produce an easily applied heuristic tool that has been successfully applied in NSW Australia.     

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.     

Landscape evaluation in conservation: molecular sampling and habitat modeling for the Iberian lynx.

Conservation of endangered species requires comprehensive understanding of their distribution and habitat requirements, in order to implement better management strategies. Unfortunately, this understanding is often difficult to gather at the short term required by rapidly declining populations of many rare vertebrates. We present a spatial habitat modeling approach that integrates a molecular technique for species detection with landscape information to assess habitat requirements of a critically endangered mammalian carnivore, the Iberian lynx (Lynx pardinus), in a poorly known population in Spain. We formulated a set of model hypotheses for habitat selection at the spatial scale of home ranges, based on previous information on lynx requirements of space, vegetation, and prey. To obtain the required data for model selection, we designed a sampling protocol based on surveys of feces and their molecular analysis for species identification. After comparing candidate models, we selected a parsimonious one that allowed (1) reliable assessment of lynx habitat requirements at the scale of home ranges, (2) prediction of lynx distribution and potential population size, and (3) identification of landscape management priorities for habitat conservation. This model predicted that the species was more likely to occur in landscapes with a higher percentage of rocky areas and higher cover of bushes typical of mature mediterranean shrubland mosaics. Its accuracy for discriminating lynx presence was approximately 85%, indicating high predictive performance. Mapping model predictions showed that only 16% of the studied areas constitute potential habitat for lynx, even though the region is dominated by large extents of well-preserved native vegetation with low human interference. Habitat was mostly clumped in two nearby patches connected by vegetation adequate for lynx dispersal and had a capacity for 28-62 potential breeding territories. The lynx population in Sierra Morena is probably the largest persisting today, but it is still critically small for optimism about its long-term persistence. Model results suggest habitat conservation and restoration actions needed for preserving the species, including reconciliation of hunting management with preservation of mature shrubland over large areas (particularly in rocky landscapes). The approach presented here can be applied to many other species for which the ecological information needed to develop sound habitat conservation strategies is lacking.