Fish population dynamics in a seasonally varying wetland

Small fishes in seasonally flooded environments such as the Everglades are capable of spreading into newly flooded areas and building up substantial biomass. Passive drift cannot account for the rapidity of observed population expansions. To test the ‘reaction-diffusion’ mechanism for spread of the fish, we estimated their diffusion coefficient and applied a reaction-diffusion model. This mechanism was also too weak to account for the spatial dynamics. Two other hypotheses were tested through modeling. The first—the ‘refuge mechanism’—hypothesizes that small remnant populations of small fishes survive the dry season in small permanent bodies of water (refugia), sites where the water level is otherwise below the surface. The second mechanism, which we call the ‘dynamic ideal free distribution mechanism’ is that consumption by the fish creates a prey density gradient and that fish taxis along this gradient can lead to rapid population expansion in space. We examined the two alternatives and concluded that although refugia may play an important role in recolonization by the fish population during reflooding, only the second, taxis in the direction of the flooding front, seems capable of matching empirical observations. This study has important implications for management of wetlands, as fish biomass is an essential support of higher trophic levels.     

Simple rules for ranking and optimally managing metapopulations

We present two ‘rules of thumb’ for metapopulation management. The first identifies an explicit formula for the persistence time of the population, and thus enables the population manager to form a priority species ranking by identifying those species most at risk of extinction. The second identifies an optimal management strategy that gives direction on how to alter the colonisation rate (creation or improvement of habitat corridors) and local extinction rate (restoring habitat quality or expanding habitat) in order to maximise the persistence time under a budgetary constraint. We employ a simple stochastic version of Levins (1969) metapopulation model, which is first calibrated to a more realistic spatial model. Our rules are tested on computer-generated patch networks and a model for malleefowl (Leipoa ocellata) in the Bakara region of South Australia.     

Predicting local population distributions around a central shelter based on a predation risk-growth trade-off

Animals face trade-offs between predation risk and foraging success depending on their location in the landscape; for example, individuals that remain near a common shelter may be safe from predation but incur stronger competition for resources. Despite a long tradition of theoretical exploration of the relationships among foraging success, conspecific competition, predation risk, and population distribution in a heterogeneous environment, the scenario we describe here has not been explored theoretically. We construct a model of habitat use rules to predict the distribution of a local population (prey sharing a common shelter and foraging across surrounding habitats). Our model describes realized habitat quality as a ratio of density- and location-dependent mortality to density-dependent growth. We explore how the prey distribution around a shelter is expected to change as the parameters governing the strength of density dependence, landscape characteristics, and local abundance vary. Within the range of parameters where prey spend some time away from shelter but remain site-attached, the prey density decreases away from shelter. As the distance at which prey react to predators increases, the population range generally increases. At intermediate reaction distances, however, increases in the reaction distance lead to decreases in the maximum foraging distance because of increased evenness in the population distribution. As total abundance increases, the population range increases, average population density increases, and realized quality decreases. The magnitude of these changes differs in, for example, ‘high-’ and ‘low-visibility’ landscapes where prey can detect predators at different distances.     

Estimating the critical phosphorus loading of shallow lakes with the ecosystem model PCLake: Sensitivity, calibration and uncertainty

There is a vast body of knowledge that eutrophication of lakes may cause algal blooms. Among lakes, shallow lakes are peculiar systems in that they typically can be in one of two contrasting (equilibrium) states that are self-stabilizing: a ‘clear’ state with submerged macrophytes or a ‘turbid’ state dominated by phytoplankton. Eutrophication may cause a switch from the clear to the turbid state, if the P loading exceeds a critical value. The ecological processes governing this switch are covered by the ecosystem model PCLake, a dynamic model of nutrient cycling and the biota in shallow lakes. Here we present an extensive analysis of the model, using a three-step procedure. (1) A sensitivity analysis revealed the key parameters for the model output. (2) These parameters were calibrated on the combined data on total phosphorus, chlorophyll-a, macrophytes cover and Secchi depth in over 40 lakes. This was done by a Bayesian procedure, giving a weight to each parameter setting based on its likelihood. (3) These weights were used for an uncertainty analysis, applied to the switchpoints (critical phosphorus loading levels) calculated by the model. The model was most sensitive to changes in water depth, P and N loading, retention time and lake size as external input factors, and to zooplankton growth rate, settling rates and maximum growth rates of phytoplankton and macrophytes as process parameters. The results for the ‘best run’ showed an acceptable agreement between model and data and classified nearly all lakes to which the model was applied correctly as either ‘clear’ (macrophyte-dominated) or ‘turbid’ (phytoplankton-dominated). The critical loading levels for a standard lake showed about a factor two uncertainty due to the variation in the posterior parameter distribution. This study calculates in one coherent analysis uncertainties in critical phosphorus loading, a parameter that is of great importance to water quality managers.     

Is there a double-dividend from anti-sprawl policies?

This paper examines the welfare effects of anti-sprawl policies, such as development tax, in a simple spatial explicit urban model with two market failures - urban decline at the city core and underprovision of open space amenities at the urban fringe - and pre-existing distortionary property tax, used to fund public services and improvements to mitigate urban decline. Consistent with prior double-dividend literature, there is a tax interaction effect that occurs between the development tax and the pre-existing property tax. However, there are two fundamental differences between the tax interaction effect identified here and that of prior literature. Ours one has two components: First, there is a cost-side tax interaction effect that is ‘spatially’ concentrated at the urban fringe, as only agents at the urban fringe alter their behavior in response to the development tax. Second, there is also a benefit-side tax interaction effect, as increases in open space at the urban fringe are capitalized into housing prices throughout the city. In contrast to prior literature, we find that the empirical importance of the combined tax interaction effect is of substantially less importance and, as a consequence, the likelihood of a ‘double-dividend’ is higher than in prior studies. Further, we show that the development tax should be part of the local tax system, even in the absence of open space benefits.     

Appreciating Ecological Complexity: Habitat Contours as a Conceptual Landscape Model

Abstract:  Organisms respond to their surroundings at multiple spatial scales, and different organisms respond differently to the same environment. Existing landscape models, such as the "fragmentation model" (or patch-matrix-corridor model) and the "variegation model," can be limited in their ability to explain complex patterns for different species and across multiple scales. An alternative approach is to conceptualize landscapes as overlaid species-specific habitat contour maps. Key characteristics of this approach are that different species may respond differently to the same environmental conditions and at different spatial scales. Although similar approaches are being used in ecological modeling, there is much room for habitat contours as a useful conceptual tool. By providing an alternative view of landscapes, a contour model may stimulate more field investigations stratified on the basis of ecological variables other than human-defined patches and patch boundaries. A conceptual model of habitat contours may also help to communicate ecological complexity to land managers. Finally, by incorporating additional ecological complexity, a conceptual model based on habitat contours may help to bridge the perceived gap between pattern and process in landscape ecology. Habitat contours do not preclude the use of existing landscape models and should be seen as a complementary approach most suited to heterogeneous human-modified landscapes.     

Modeling the brown bear population in Slovenia: A tool in the conservation management of a threatened species

In this paper, we address three aspects of the brown bear population in Slovenia: its size (and its evolution over time), its spatial expansion out of the core area, and its potential habitat based on natural habitat suitability. Data collected through measurement/observation of the bear population and from the literature are used. A model is developed for each aspect. The results are estimates of population size, a picture of the spatial expansion of the population and maps of its optimal and maximal potential habitat (based on natural suitability). Overall, the brown bear population has been increasing since the establishment of a core protective area and has been expanding outside this area. The habitat suitability maps show that there is room for further expansion. Based on habitat suitability and bear population density, as well as human activity and current damage reports, we recommend that the Alps should be temporarily kept free of the bears, until the necessary mitigation measures regarding human–bear conflicts are carried out. On the other hand it is of crucial importance to adapt human activities and improve bear management in the optimal habitat, with which the goals of successful conservation of the species might be achieved.     

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.     

Downscaling species occupancy from coarse spatial scales

The measurement and prediction of species' populations at different spatial scales is crucial to spatial ecology as well as conservation biology. An efficient yet challenging goal to achieve such population estimates consists of recording empirical species' presence and absence at a specific regional scale and then trying to predict occupancies at finer scales. So far the majority of the methods have been based on particular species' distributional features deemed to be crucial for downscaling occupancy. However, only a minority of them have dealt explicitly with specific spatial features. Here we employ a wide class of spatial point processes, the shot noise Cox processes (SNCP), to model species occupancies at different spatial scales and show that species' spatial aggregation is crucial for predicting population estimates at fine scales starting from coarser ones. These models are formulated in continuous space and locate points regardless of the arbitrary resolution that one employs to study the spatial pattern. We compare the performances of nine models, calibrated at regional scales and demonstrate that a very simple class of SNCP, the Thomas process, is able to outperform other published models in predicting occupancies down to areas four orders of magnitude smaller than the ones employed for the parameterization. We conclude by explaining the ability of the approach to infer spatially explicit information from spatially implicit measures, the potential of the framework to combine niche and spatial models, and the possibility of reversing the method to allow upscaling.     

Hierarchical Bayesian space-time models

Space-time data are ubiquitous in the environmental sciences. Often, as is the case with atmo- spheric and oceanographic processes, these data contain many different scales of spatial and temporal variability. Such data are often non-stationary in space and time and may involve many observation/prediction locations. These factors can limit the effectiveness of traditional space- time statistical models and methods. In this article, we propose the use of hierarchical space-time models to achieve more flexible models and methods for the analysis of environmental data distributed in space and time. The first stage of the hierarchical model specifies a measurement- error process for the observational data in terms of some 'state' process. The second stage allows for site-specific time series models for this state variable. This stage includes large-scale (e.g. seasonal) variability plus a space-time dynamic process for the anomalies'. Much of our interest is with this anomaly proc ess. In the third stage, the parameters of these time series models, which are distributed in space, are themselves given a joint distribution with spatial dependence (Markov random fields). The Bayesian formulation is completed in the last two stages by speci- fying priors on parameters. We implement the model in a Markov chain Monte Carlo framework and apply it to an atmospheric data set of monthly maximum temperature.     

Modeling high-frequency position data of large herbivores with a phase-state model

Understanding the rules and factors that drive the foraging behavior of large herbivores is important to describe their interaction with the landscape at various spatial scales. Some unresolved questions refer to landscape-behavioral interactions that result in oriented or random search in seasonally changing landscapes. Remotely sensed position data indicate that herbivores select local patches of heterogeneous landscapes depending on a complex host of dynamically varying animal and environmental conditions. Since foraging paths consist in successions of relatively short steps, increasing the frequency at which position information is acquired would contribute to entangle the mechanisms resulting in herbivores’ foraging paths. We addressed the question whether herbivores would obtain information at a patch scale that would modify their distribution at a landscape scale based on directed movement or navigation ability. We considered a set of 100,000 high-frequency (1 min intervals) position data of several free-ranging sheep (Ovis aries) at a seasonal-varying range (Patagonian Monte, Argentina) and observed their movements at landscape and at single vegetation patch scales. At a landscape scale, we inspected the spatial co-variation of seasonally varying forage offer and ewes’ movement speeds. At a patch scale, we developed a phase-state (P-S) model of movement cycles based on the occurrence of behavioral phases along foraging paths, and fitted it to the observed daily time series of ewes’ movement speeds. Ewes were preferentially distributed in areas with high forage offer during periods of low forage availability and the reverse occurred during the season of high forage availability. Parameters of the model of activity cycles amenable to control by ewes (duration of speed phases, time elapsed between speed cycles) did not covariate with forage offer, but varied significantly among ewes. The shape (kurtosis) parameter of the model of movement cycles, one which is unlikely under ewes’ control, co-varied significantly with spatial forage offer but did not differ among ewes. We conclude that ewes allocated foraging time along a series of similar movement efforts irrespective of forage availability at small patches. Average forage scarcity at multi-patch level increases the ratio of searching to feeding time. This results in apparent selective time allocation to richer forage areas but does not imply evidence for oriented movement at a landscape scale. We advance a behavioral-based definition of forage patches and discuss its implications in developing foraging theory and models. The P-S model applied to high-frequency position data of large herbivores substantially improves the interpretation of the factors controlling their time allocation in space with respect to previous models of herbivore spatial behavior by discriminating among behavioral-based and environmentally induced components of their movements.     

Bird‐community responses to habitat creation in a long‐term,large‐scale natural experiment

Ecosystem function and resilience are compromised when habitats become fragmented due to land‐use change. This has led to national and international conservation strategies aimed at restoring habitat extent and improving functional connectivity (i.e., maintaining dispersal processes). However, biodiversity responses to landscape‐scale habitat creation and the relative importance of spatial and temporal scales are poorly understood, and there is disagreement over which conservation strategies should be prioritized. We used 160 years of historic post‐agricultural woodland creation as a natural experiment to evaluate biodiversity responses to habitat creation in a landscape context. Birds were surveyed in 101 secondary, broadleaf woodlands aged 10–160 years with ≥80% canopy cover and in landscapes with 0‐17% broadleaf woodland cover within 3000 m. We used piecewise structural equation modeling to examine the direct and indirect relationships between bird abundance and diversity, ecological continuity, patch characteristics, and landscape structure and quantified the relative conservation value of local and landscape scales for bird communities. Ecological continuity indirectly affected overall bird abundance and species richness through its effects on stand structure, but had a weaker influence (effect size near 0) on the abundance and diversity of species most closely associated with woodland habitats. This was probably because woodlands were rapidly colonized by woodland generalists in ≤10 years (minimum patch age) but were on average too young (median 50 years) to be colonized by woodland specialists. Local patch characteristics were relatively more important than landscape characteristics for bird communities. Based on our results, biodiversity responses to habitat creation depended on local‐ and landscape‐scale factors that interacted across time and space. We suggest that there is a need for further studies that focus on habitat creation in a landscape context and that knowledge gained from studies of habitat fragmentation and loss should be used to inform habitat creation with caution because the outcomes are not necessarily reciprocal.     

Transcending scale dependence in identifying habitat with resource selection functions

Multi-scale resource selection modeling is used to identify factors that limit species distributions across scales of space and time. This multi-scale nature of habitat suitability complicates the translation of inferences to single, spatial depictions of habitat required for conservation of species. We estimated resource selection functions (RSFs) across three scales for a threatened ungulate, woodland caribou (Rangifer tarandus caribou), with two objectives: (1) to infer the relative effects of two forms of anthropogenic disturbance (forestry and linear features) on woodland caribou distributions at multiple scales and (2) to estimate scale-integrated resource selection functions (SRSFs) that synthesize results across scales for management-oriented habitat suitability mapping. We found a previously undocumented scale-specific switch in woodland caribou response to two forms of anthropogenic disturbance. Caribou avoided forestry cut-blocks at broad scales according to first- and second-order RSFs and avoided linear features at fine scales according to third-order RSFs, corroborating predictions developed according to predator-mediated effects of each disturbance type. Additionally, a single SRSF validated as well as each of three single-scale RSFs when estimating habitat suitability across three different spatial scales of prediction. We demonstrate that a single SRSF can be applied to predict relative habitat suitability at both local and landscape scales in support of critical habitat identification and species recovery.     

An ecological model of the habitat mosaic in estuarine nursery areas: Part I—Interaction of dispersal theory and habitat variability in describing juvenile fish distributions

Identification of critical habitat in estuarine nursery areas is an important conservation and management objective. Habitat can be viewed as a mosaic of both temporally variable environmental features and spatially variable structural features that combine to define optimal habitat. Effective models of juvenile distributions should account for individual movement, as well as the full suite of habitat variability including both spatial and temporal components. We have extended a terrestrial model of small-scale movement patterns to describe habitat choices of an index juvenile fish in an estuarine nursery system. Movement of small juvenile fishes was found to be influenced by both spatial and temporal patterns in habitat quality, and it was a balanced mix of both that resulted in an optimal distribution. Fishes that perceive habitat on a scale much smaller than the scale of spatial heterogeneity may respond to temporal change as a movement cue allowing for more deterministic outcomes at larger scales despite perceptual limitations. These model outcomes suggest a hierarchical approach is best for describing habitat choice in juvenile fishes and this approach will be used in the future to explore individual and population responses to predictable habitat change.     

Tradeoffs between forestry resource and conservation values under alternate policy regimes: A spatial analysis of the western Canadian boreal plains

An important element of resource management and conservation is an understanding of the tradeoffs between marketed products, such as timber, and measures of environmental quality, such as biodiversity. In this paper, we develop an integrated economic-ecological spatial optimization model that we then apply to evaluate alternate forest policies on a 560,000 km² study region of managed boreal forest in Alberta and British Columbia, Canada. The integrated model incorporates dynamic forest sector harvesting, current levels of oil and gas sector development, coarse-filter or habitat-based old forest indicators, a set of empirical forest bird abundance models, and statistical models of the natural and current fire regimes. Using our integrated model, economic tradeoff curves, or production possibility frontiers, are developed to illustrate the cost of achieving coarse-filter targets by a set time (50 years) within a 100-year time horizon. We found levels of ecological indicators and economic returns from the timber industry could both be increased if spatial constraints imposed by the current policy environment were relaxed; other factors being equal, this implies current policy should be revised. We explore the production possibility frontier's relationship to the range of natural variation of old forest habitat, and show how this range can be used to guide choices of preferred locations along the frontier. We also show that coarse-filter constraints on the abundance of certain habitat elements are sufficient to satisfy some fine-filter objectives, expressed as the predicted abundances of various species of songbirds.     

Influence of grain size on species-habitat models

High resolution remote sensing data facilitate the use of small-scale habitat features such as trees or hedges in the analysis of species-habitat relationships. Such data potentially enable more accurate species-habitat mapping than lower resolution data. Here, for the first time, we systematically investigated this hypothesis by altering the spatial resolution from 1 m up to 1000 m grain size in species-habitat models of 13 bird species. The study area covered the Nidda river catchment in central Germany, a large heterogeneous landscape of 1620 km². A high resolution habitat map of the area was converted to coarser spatial and thematic resolutions in seven steps. We investigated how model performance responded to grain size, and we compared the differential effects of spatial resolution and thematic resolution on model performance. Explained deviance (D²) of the bird models generally decreased with coarser spatial resolution of the data, although it did not decrease monotonically in all species. On average across all species, model D² decreased from 41.5 at 1 m grain size to 15.9 at 1000 m grain size. Ten species were best modelled at 1 m, two species at 3 m and one species at 32 m grain size. Model performance degraded continuously with increasing grain size, both in habitat generalist and habitat specialist bird species, and was systematically lower in habitat generalists. The higher model performance observed at finer grain sizes was most likely caused by the combination of three factors: (1) high spatial accuracy of bird records and (2) a more precise location and delineation of habitat features and, (3) to a lesser degree, by more habitat types differentiated in maps of finer resolution. We conclude that higher spatial and thematic resolution data can be essential for deriving accurate predictions on bird distribution patterns from species-habitat models. Especially for bird species that are sensitive to specific land-use types or to small-scaled habitat features, a grain size of 1-3 m seems most promising.     

Examination of model predictions at different horizontal grid resolutions

While fluctuations in meteorological and air quality variables occur on a continuum of spatial scales, the horizontal grid spacing of coupled meteorological and photochemical models sets a lower limit on the spatial scales that they can resolve. However, both computational costs and data requirements increase significantly with increasing grid resolution. Therefore, it is important to examine, for any given application, whether the expected benefit of increased grid resolution justifies the extra costs. In this study, we examine temperature and ozone observations and model predictions for three high ozone episodes that occurred over the northeastern United States during the summer of 1995. In the first set of simulations, the meteorological model RAMS4a was run with three two-way nested grids of 108/36/12 km grid spacing covering the United States and the photochemical model UAM-V was run with two grids of 36/12 km grid spacing covering the eastern United States. In the second set of simulations, RAMS4a was run with four two-way nested grids of 108/36/12/4 km grid spacing and UAM-V was run with three grids of 36/12/4 km grid spacing with the finest resolution covering the northeastern United States. Our analysis focuses on the comparison of model predictions for the finest grid domain of the simulations, namely, the region overlapping the 12 km and 4 km domains. A comparison of 12 km versus 4 km fields shows that the increased grid resolution leads to finer texture in the model predictions; however, comparisons of model predictions with observations do not reveal the expected improvement in the predictions. While high-resolution modeling has scientific merit and potential uses, the currently available monitoring networks, in conjunction with the scarceness of highly resolved spatial input data and the limitations of model formulation, do not allow confirmation of the expected superiority of the high-resolution model predictions.The U.S. Governments right to retain a non-exclusive royalty-free licence in and to any copyright is acknowledged.     

Fragmentation profiles for real and simulated landscapes

When a natural landscape is represented by a series of categorical raster maps of varying resolution, a multiresolution characterization of spatial pattern can be obtained in which entropy is computed at each resolution conditional on the next coarser resolution. The series of entropy values is plotted as a function of resolution, resulting in a multiresolution profile of fragmentation pattern in the landscape. If a categorical raster map is available at a single resolution only, a series of degraded maps at increasingly coarser resolutions is generated and the fragmentation profile is computed for this series. An algorithm has been developed for obtaining the profile directly from the single resolution map without having to generate and store the coarser resolution maps. A hierarchical stochastic model is described for simulating categorical raster maps and the fragmentation profile of the generating process is obtained in terms of the model parameters. These process profiles provide benchmarks for assessing empirical profiles obtained from raster maps of actual landscapes. Methods of the paper are applied to several watersheds of Pennsylvania using landcover maps derived from satellite imagery. These examples indicate that characteristic landscape types induce characteristic features in their fragmentation profiles.     

Linking chronic wasting disease to mule deer movement scales: a hierarchical Bayesian approach.

Observed spatial patterns in natural systems may result from processes acting across multiple spatial and temporal scales. Although spatially explicit data on processes that generate ecological patterns, such as the distribution of disease over a landscape, are frequently unavailable, information about the scales over which processes operate can be used to understand the link between pattern and process. Our goal was to identify scales of mule deer (Odocoileus hemionus) movement and mixing that exerted the greatest influence on the spatial pattern of chronic wasting disease (CWD) in northcentral Colorado, USA. We hypothesized that three scales of mixing (individual, winter subpopulation, or summer subpopulation) might control spatial variation in disease prevalence. We developed a fully Bayesian hierarchical model to compare the strength of evidence for each mixing scale. We found strong evidence that the finest mixing scale corresponded best to the spatial distribution of CWD infection. There was also evidence that land ownership and habitat use play a role in exacerbating the disease, along with the known effects of sex and age. Our analysis demonstrates how information on the scales of spatial processes that generate observed patterns can be used to gain insight when process data are sparse or unavailable.