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.     

Scale dependence of reproductive failure in fragmented Echinacea populations

I investigated reproduction in a three-year study of Echinacea angustifolia, purple coneflower, growing in a fragmented prairie landscape. I quantified the local abundance of flowering conspecifics at individual-based spatial scales and at a population-based spatial scale. Regression analyses revealed that pollen limitation increased while seed set and fecundity decreased with isolation of individual plants. Isolation, defined as the distance to the k(th) nearest flowering conspecific, was a good predictor of pollen limitation, for all nearest neighbors considered (k = 1-33), but the strength of the relationship, as quantified by R2, peaked at intermediate scales (k = 2-18). The relationship of isolation to seed set and fecundity was similarly strongest at intermediate scales (k = 3-4). The scale dependence of individual density effects on reproduction (density of flowering plants within x meters) resembled that of isolation. Analyses at a population-based scale showed that pollen limitation declined significantly with population size. Seed set and fecundity also declined with population size, but significantly so only in 1998. Whether quantifying local abundance with population- or individual-based measures, reproductive failure due to pollen limitation is a consistent consequence of Echinacea scarcity. However, individual-based measures of local abundance predicted pollen limitation from a wider sample of plants with a simpler model than did population size. Specifically, the largest site, a nature preserve, is composed of plants with intermediate individual isolation and, as predicted, intermediate pollen limitation, but its large population size poorly predicted population mean pollen limitation.     

From population sources to sieves: the matrix alters host-parasitoid source-sink structure

Field experiments that examine the impact of immigration, emigration, or landscape structure (e.g., the composition of the matrix) on the source sink dynamics of fragmented populations are scarce. Here, planthoppers (Prokelisia crocea) and egg parasitoids (Anagrus columbi) were released among host-plant patches that varied in structural (caged, isolated, or in a network of other patches) and functional (mudflat matrix that impedes dispersal vs. brome-grass matrix that facilitates dispersal) connectivity. Planthoppers and parasitoids on caged patches exhibited density-dependent growth rates, achieved high equilibrium densities, and rarely went extinct. Therefore, experimental cordgrass patches were classified as population sources. Because access to immigrants did not result in elevated population densities, source populations were not also pseudosinks, i.e., patches whose densities occur above carrying capacity due to high immigration. Planthoppers and parasitoids in open patches in mudflat had dynamics similar to those in caged patches, but went extinct in 4-5 generations in open patches in brome. Brome-embedded patches leaked emigrants at a rate that exceeded the gains from reproduction and immigration; populations of this sort are known as population sieves. For species whose suitable patches are becoming smaller and more isolated as a result of increased habitat fragmentation, emigration losses are likely to become paramount, a condition favoring the formation of population sieves. An increase in the proportion of patches that are sieves is predicted to destabilize regional population dynamics.     

Connectivity in an Agricultural Landscape as Reflected by Interpond Movements of a Freshwater Turtle

Abstract:  Connectivity is a measure of how landscape features facilitate movement and thus is an important factor in species persistence in a fragmented landscape. The scarcity of empirical studies that directly quantify species movement and determine subsequent effects on population density have, however, limited the utility of connectivity measures in conservation planning. We undertook a 4-year study to calculate connectivity based on observed movement rates and movement probabilities for five age-sex classes of painted turtles ( Chrysemys picta ) inhabiting a pond complex in an agricultural landscape in northern Virginia (U.S.A.). We determined which variables influenced connectivity and the relationship between connectivity and subpopulation density. Interpatch distance and quality of habitat patches influenced connectivity but characteristics of the intervening matrix did not. Adult female turtles were more influenced by the habitat quality of recipient ponds than other age-sex classes. The importance of connectivity on spatial population dynamics was most apparent during a drought. Population density and connectivity were low for one pond in a wet year but dramatically increased as other ponds dried. Connectivity is an important component of species persistence in a heterogeneous landscape and is strongly dependent on the movement behavior of the species. Connectivity may reflect active selection or avoidance of particular habitat patches. The influence of habitat quality on connectivity has often been ignored, but our findings highlight its importance. Conservation planners seeking to incorporate connectivity measures into reserve design should not ignore behavior in favor of purely structural estimates of connectivity.     

Integrating Landscape and Metapopulation Modeling Approaches: Viability of the Sharp-Tailed Grouse in a Dynamic Landscape

Abstract:  The lack of management experience at the landscape scale and the limited feasibility of experiments at this scale have increased the use of scenario modeling to analyze the effects of different management actions on focal species. However, current modeling approaches are poorly suited for the analysis of viability in dynamic landscapes. Demographic (e.g., metapopulation) models of species living in these landscapes do not incorporate the variability in spatial patterns of early successional habitats, and landscape models have not been linked to population viability models. We link a landscape model to a metapopulation model and demonstrate the use of this model by analyzing the effect of forest management options on the viability of the Sharp-tailed Grouse (  Tympanuchus phasianellus ) in the Pine Barrens region of northwestern Wisconsin (U.S.A.). This approach allows viability analysis based on landscape dynamics brought about by processes such as succession, disturbances, and silviculture. The landscape component of the model (LANDIS) predicts forest landscape dynamics in the form of a time series of raster maps. We combined these maps into a time series of patch structures, which formed the dynamic spatial structure of the metapopulation component (RAMAS). Our results showed that the viability of Sharp-tailed Grouse was sensitive to landscape dynamics and demographic variables such as fecundity and mortality. Ignoring the landscape dynamics gave overly optimistic results, and results based only on landscape dynamics (ignoring demography) lead to a different ranking of the management options than the ranking based on the more realistic model incorporating both landscape and demographic dynamics. Thus, models of species in dynamic landscapes must consider habitat and population dynamics simultaneously.     

Computer simulation of barren-ground caribou dynamics

This paper describes a simulation model of a Canadian caribou population. The model was constructed by an interdisciplinary team of field biologists, managers, and systems ecologists. Population dynamics are represented in terms of age structure, with age-dependent survival and fecundity. Biomass dynamics of the major food species are also simulated, and the food dynamics interact with the caribou population through a foraging submodel that explicitly considers snow depth, seasonal migrations, and total area of useable winter habitat. The model was used to examine two hypotheses regarding the abundance of barren-ground caribou. It was shown that there is no reason to suspect that food supply currently limits the population size; instead, hunting pressure appears to be the critical variable. The implications of this finding for population management are discussed.     

Modeling the spatio-temporal dynamics and interactions of households,landscapes, and giant panda habitat

Human modification of land-cover has been a leading cause of floral and faunal species extirpation and loss of local and global biodiversity. As natural areas are impacted, habitat and populations can become fragmented and isolated. This is particularly evident in the mountainous areas of southwestern China that support the remaining populations of giant pandas (Ailuropoda melanoleuca). Giant panda populations have been restricted to remnants of habitat from extensive past land use and land-cover change. Households are a basic socio-economic unit that continues to impact the remaining habitat through activities such as fuelwood consumption and new household creation. Therefore, we developed a spatio-temporal model of human activities and their impacts by directly integrating households into the landscape. The integrated model allows us to examine the landscape factors influencing the spatial distribution of household activities and household impacts on habitat. As an example application, we modeled household activities in a giant panda reserve in China and examined the spatio-temporal dynamics of households, the landscape, and their mutual interactions. Human impacts are projected to result in the loss of up to 16% of all existing habitat within the reserve over the next 30 years. In addition, we found that accessibility largely controls the spatial distribution of household activities and considerable changes in management and household activities will be required to maintain the current level of habitat within the reserve.     

Landscape Patterns of Florida Scrub Jay Habitat Use and Demographic Success

Remote sensing and geographical information systems are used to analyze landscapes important to species of conservation concern. The accuracy of the methods depends on how closely habitat mapping classes are linked to population demography. Habitat use by Florida Scrub Jays ( Aphelocoma c. coemlescens ) was quantified using circular plots. Habitat variation was mapped using high-resolution aerial photography on a site where all Florida Scrub Jays were color-banded. Nest site selection, nest success, yearling production, and breeder survival were measured within Florida Scrub Jay territories. Habitat use was lowest in areas without scrub oaks or areas within 136 m from forests. Open oak, dominated by scrub oaks and open sandy areas, had the highest use and nest success among habitats. Open oak occurred as narrow patches ( <20 m wide) in landscapes dominated by matrix habitat (palmetto-lyonia and swale marshes). Most wide patches (>50 m) of open oak were potential population sources, where reproduction exceeded mortality. Areas with patches of open oak of less than 1 ha were usually population sinks. Open oak occurred as less flammable patches in a landscape subject to frequent fires. Population sources varied temporally and spatially with fires and site potential to support scrub oaks (soils). Analyses of landscape patterns and dynamics indicated that habitat mapping should not only include patches of currently optimal habitat but should also include landscapes associated with open oak. The influences of landscape patterns on habitat use, reproductive success, survival, and territory size can be quantified at different scales starting with attributes associated with habitat patches, nest sites, and territories. Potential mapping errors occur, however, when habitats are used to quantify the areal extent of sources and sinks and similar population attributes important for species persistence.     

Extinction Dynamics in the American Marten (Martes americana)

We constructed a model of marten population dynamics and used it to investigate extinction processes across a wide range of parameter values. The model was based on rules governing the behavior and physiology of individual martens and focused on energy balance. Spatial dynamics and demographic and environmental stochasticity were incorporated. The outcome was the probability of extinction and quasiextinction (20 females remaining) over 500 years. Three qualitative forms of extinction were delineated. The first was deterministic extinction, associated with those parameter combinations leading to a negative population growth rate. The second was probabilistic extinction in systems with a strong positive growth rate but restricted population size due to habitat constraint. The transition from 100% persistence to 100% quasiextinction, as the input habitat size was decreased, was abrupt. The final form of extinction was in systems with a growth rate of approximately zero. Prey availability maintained an upper limit on these populations, but otherwise fluctuations in population size were essentially random, leading to nontrivial probabilities of extinction in even relatively large populations. A number of issues requiring further empirical research were identified. These included the relationship between habitat quality and marten reproduction, dispersal patterns and dispersal mortality, the effect of habitat edge on marten reproduction and mortality, and the characterization of the severity and frequency of catastrophic mortality as experienced by marten populations.     

Connected fragmented habitats facilitate stable coexistence dynamics

In this paper we endeavor to test the controversial ideas that exist about the role of fragmentation in a conservation context. In line with earlier understanding, we find that habitat fragmentation alone results in a strong detrimental effect (especially for the predator population). Connecting the fragmented habitats facilitates predator survival and hence prey survival as compared to the unconnected fragmented case. Our main result is counterintuitive: in the presence of a high quality predator, connected fragmented habitats ensure a better chance for coexistence than does even the unfragmented case. We explain why a connected fragmented habitat might thus be beneficial for the stabilization of the system, and how connections between sub-habitats are able to protect prey population from over-exploitation. In the model, habitat fragmentation is separated from the effects of habitat destruction, in order to better understand how populations react to habitat transformation.     

Population viability of the Siberian Tiger in a changing landscape: Going, going and gone?

The Amur tiger (Panthera tigris altaica) is a flagship species of the boreal forest ecosystem in northeastern China and Russia Far East. During the past century, the tiger population has declined sharply from more than 3000 to fewer than 600 individuals, and its habitat has become much smaller and greatly fragmented. Poaching, habitat degradation, habitat loss, and habitat fragmentation have been widely recognized as the primary causes for the observed population decline. Using a population viability analysis tool (RAMAS/GIS), we simulated the effects of poaching, habitat degradation, habitat loss, and habitat fragmentation on the population dynamics and extinction risk of the Amur tiger, and then explored the relative effectiveness of three conservation strategies involving improving habitat quality and establishing movement corridors in China and Russia. A series of controlled simulation experiments were performed based on the current spatial distribution of habitat and field-observed vital rates. Our results showed that the Amur tiger population could be viable for the next 100 years if the current habitat area and quality were well-maintained, with poaching strictly prohibited of the tigers and their main prey species. Poaching and habitat degradation (mainly prey scarcity) had the largest negative impacts on the tiger population persistence. While the effect of habitat loss was also substantial, habitat fragmentation per se had less influence on the long-term fate of the tiger population. However, to sustain the subpopulations in both Russia and China would take much greater conservation efforts. The viability of the Chinese population of tigers would rely heavily on its connectivity with the largest patch on the other side of the border. Improving the habitat quality of small patches only or increasing habitat connectivity through movement corridors alone would not be enough to guarantee the long-term population persistence of the Amur tiger in both Russia and China. The only conservation strategy that allowed for long-term persistence of tigers in both countries required both the improvement of habitat quality and the establishment of a transnational reserve network. Our study provides new insights into the metapopulation dynamics and persistence of the Amur tiger, which should be useful in landscape and conservation planning for protecting the biggest cat species in the world.     

Using experimental reintroductions to resolve the roles of habitat quality and metapopulation dynamics on patch occupancy in fragmented landscapes

Declines of species in fragmented landscapes can potentially be reversed either by restoring connectivity or restoring local habitat quality. Models fitted to snapshot occupancy data can be used to predict the effectiveness of these actions. However, such inferences can be misleading if the reliability of the habitat and landscape metrics used is unknown. The only way to unambiguously resolve the roles of habitat quality and metapopulation dynamics is to conduct experimental reintroductions to unoccupied patches so that habitat quality can be measured directly from data on vital rates. We, therefore, conducted a 15-year study that involved reintroducing a threatened New Zealand bird to unoccupied forest fragments to obtain reliable data on their habitat quality and reassess initial inferences made by modeling occupancy against habitat and landscape metrics. Although reproductive rates were similar among fragments, subtle differences in adult survival rates resulted in λ (finite rate of increase) estimations of <0.9 for 9 of the 12 fragments that were previously unoccupied. This was the case for only 1 of 14 naturally occupied fragments. This variation in λ largely explained the original occupancy pattern, reversing our original conclusion from occupancy modeling that this occupancy pattern was isolation driven and suggesting that it would be detrimental to increase connectivity without improving local habitat quality. These results illustrate that inferences from snapshot occupancy should be treated with caution and subjected to testing through experimental reintroductions in selected model systems.     

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.     

Fragmention by agriculture influences reproductive success of birds in a shrubsteppe landscape.

Shrubsteppe communities are among the most imperiled ecosystems in North America as a result of conversion to agriculture and other anthropogenic changes. In the Intermountain West of the United States, these communities support a unique avifauna, including several species that are declining and numerous others that are of conservation concern. Extensive research in the eastern and central United States and in Scandinavia suggests that fragmentation of formerly continuous forests and grasslands adversely affects reproductive success of birds, yet little is known of the potential effects on avian communities in Western shrublands. I used multi-model inference to evaluate the potential effects of local and landscape variables on nest predation and brood parasitism, and behavioral observations of color-banded birds to evaluate the potential effects of habitat fragmentation on seasonal reproductive success of passerines in the shrubsteppe of eastern Washington State, USA. Reproductive success of shrubsteppe-obligate passerines was lower in landscapes fragmented by agriculture than in continuous shrubsteppe landscapes. Daily survival rates for nests of Brewer's Sparrows (Spizella breweri; n=496) and Sage Thrashers (Oreoscoptes montanus; n=128) were lower in fragmented landscapes, and seasonal reproductive success (percentage of pairs fledging young) of Sage Sparrows (Amphispiza belli; n=146) and Brewer's Sparrows (n=59) was lower in fragmented landscapes. Rates of parasitism by Brown-headed Cowbirds (Molothrus ater) overall were low (4%) but were significantly greater in fragmented landscapes for Brewer's Sparrows, and parasitism resulted in fewer young fledged from successful nests. Simple models of population growth using landscape-specific fecundity and estimates of adult survival derived from return rates of banded male Sage Sparrows and Brewer's Sparrows suggest that fragmented shrubsteppe in Washington may be acting as a population sink for some species. Immediate conservation needs include halting further fragmentation of shrubsteppe, restoring low-productivity agricultural lands and annual grasslands to shrubsteppe where possible, and convincing the public of the intrinsic value of these imperiled ecosystems.     

Conservation of Fragmented Populations

In this paper we argue that landscape spatial structure is of central importance in understanding the effects of fragmentation on population survival. Landscape spatial structure is the spatial relationships among habitat patches and the matrix in which they are embedded. Many general models of subdivided populations make the assumptions that (1) all habitat patches are equivalent in size and quality and (2) all local populations (in the patches) are equally accessible by dispersers. Models that gloss over spatial details of landscape structure can be useful for theoretical developments but will almost always be misleading when applied to real-world conservation problems. We show that local extinctions of fragmented populations are common. From this it follows that recolonization of local extinctions is critical for regional survival of fragmented populations. The probability of recolonization depends on (1) spatial relationships among landscape elements used by the population, including habitat patches for breeding and elements of the inter-patch matrix through which dispersers move, (2) dispersal characteristics of the organism of interest, and (3) temporal changes in the landscape structure. For endangered species, which are typically restricted in their dispersal range and in the kinds of habitat through which they can disperse, these factors are of primary importance and must be explicitly considered in management decisions.     

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.     

Resource distributions among habitats determine solitary bee offspring production in a mosaic landscape.

Within mosaic landscapes, many organisms depend on attributes of the environment that operate over scales ranging from a single habitat patch to the entire landscape. One such attribute is resource distribution. Organisms' reliance on resources from within a local patch vs. those found among habitats throughout the landscape will depend on local habitat quality, patch quality, and landscape composition. The ability of individuals to move among complementary habitat types to obtain various resources may be a critical mechanism underlying the dynamics of animal populations and ultimately the level of biodiversity at different spatial scales. We examined the effects that local habitat type and landscape composition had on offspring production and survival of the solitary bee Osmia lignaria in an agri-natural landscape in California (U.S.A.). Female bees were placed on farms that did not use pesticides (organic farms), on farms that did use pesticides (conventional farms), or in seminatural riparian habitats. We identified pollens collected by bees nesting in different habitat types and matched these to pollens of flowering plants from throughout the landscape. These data enabled us to determine the importance of different plant species and habitat types in providing food for offspring, and how this importance changed with landscape and local nesting-site characteristics. We found that increasing isolation from natural habitat significantly decreased offspring production and survival for bees nesting at conventional farms, had weaker effects on bees in patches of seminatural habitat, and had little impact on those at organic farm sites. Pollen sampled from nests showed that females nesting in both farm and seminatural habitats relied on pollen from principally native plant species growing in seminatural habitat. Thus connectivity among habitats was critical for offspring production. Females nesting on organic farms were buffered to isolation effects by switching to floral resources growing at the farm site when seminatural areas were too distant. Overall local habitat conditions (farm management practices) can help bolster pollinators, but maintaining functional connectivity among habitats will likely be critical for persistence of pollinator populations as natural habitats are increasingly fragmented by human activities.     

Simulating population variation and movement within fragmented landscapes: An application to the gopher tortoise (Gopherus polyphemus)

As the human activity footprint grows, land-use decisions play an increasing role in determining the future of plant and animal species. Studies have shown that urban and agricultural development cannot only harm species populations directly through habitat destruction, but also by destroying the corridors that connect habitat patches and populations within a metapopulation. Without these pathways, populations can encounter inbreeding depression and degeneration, which can increase death rates and lower rates of reproduction. This article describes the development and application of the FRAGGLE model, a spatial system dynamics model designed to calculate connectivity indices among populations. FRAGGLE can help planners and managers identify the relative contribution of populations associated with habitat patches to future populations in those patches, taking into account the importance of interstitial land to migration success. The model is applied to the gopher tortoise (Gopherus polyphemus), a threatened species whose southeastern U.S. distribution has diminished significantly within its native range due to agricultural and urban development over the last several decades. This model is parameterized with life history and movement traits of the gopher tortoise in order to simulate population demographics and spatial distribution within an area in west-central Georgia that supports a significant tortoise population. The implications of this simulation modeling effort are demonstrated using simple landscape representations and a hypothetical on land-use management scenario. Our findings show that development resulting in even limited habitat losses (10%) may lead to significant increases in fragmentation as measured by a loss in the rate of dispersions (31%) among area subpopulations.     

A Framework for Conceptualizing Human Effects on Landscapes and Its Relevance to Management and Research Models

Abstract: The concept of habitat fragmentation is limited in its ability to describe the range of possible landscape configurations created by a variety of disturbances. This limitation is especially problematic in landscapes where human use of the habitat matrix occurs at multiple levels and where habitat modification may be a more important consideration than a simple binary classification of habitat versus nonhabitat. We propose a synthesizing scheme that places intact, variegated, fragmented, and relictual landscape states on a continuum, depending on the degree of habitat destruction. At a second level, the scheme considers the patterns of habitat modification that are imposed on remaining habitats. Management for conservation involves halting and sometimes reversing the trends of habitat destruction and modification. Conservation strategies will differ according to the state of alteration of the landscape, but all strategies include some consideration of the degree of modification of the matrix in determining habitat viability. It is convenient for biologists to assess landscape alteration state in terms of the persistence of large structural elements such as trees. Because animal species use habitats differently, however, they also experience the landscape differently. A landscape considered structurally fragmented by humans may be functionally variegated to other species. Therefore, it is necessary to consider the extent to which the entire landscape, including the matrix, is accessible and utilized by organisms with different spatial scales of resource use.