Bayesian Networks and Adaptive Management of Wildlife Habitat

Abstract: Adaptive management is an iterative process of gathering new knowledge regarding a system's behavior and monitoring the ecological consequences of management actions to improve management decisions. Although the concept originated in the 1970s, it is rarely actively incorporated into ecological restoration. Bayesian networks (BNs) are emerging as efficient ecological decision‐support tools well suited to adaptive management, but examples of their application in this capacity are few. We developed a BN within an adaptive‐management framework that focuses on managing the effects of feral grazing and prescribed burning regimes on avian diversity within woodlands of subtropical eastern Australia. We constructed the BN with baseline data to predict bird abundance as a function of habitat structure, grazing pressure, and prescribed burning. Results of sensitivity analyses suggested that grazing pressure increased the abundance of aggressive honeyeaters, which in turn had a strong negative effect on small passerines. Management interventions to reduce pressure of feral grazing and prescribed burning were then conducted, after which we collected a second set of field data to test the response of small passerines to these measures. We used these data, which incorporated ecological changes that may have resulted from the management interventions, to validate and update the BN. The network predictions of small passerine abundance under the new habitat and management conditions were very accurate. The updated BN concluded the first iteration of adaptive management and will be used in planning the next round of management interventions. The unique belief‐updating feature of BNs provides land managers with the flexibility to predict outcomes and evaluate the effectiveness of management interventions.     

Improving inferences about private land conservation by accounting for incomplete reporting

Private lands provide key habitat for imperiled species and are core components of function protectected area networks; yet, their incorporation into national and regional conservation planning has been challenging. Identifying locations where private landowners are likely to participate in conservation initiatives can help avoid conflict and clarify trade-offs between ecological benefits and sociopolitical costs. Empirical, spatially explicit assessment of the factors associated with conservation on private land is an emerging tool for identifying future conservation opportunities. However, most data on private land conservation are voluntarily reported and incomplete, which complicates these assessments. We used a novel application of occupancy models to analyze the occurrence of conservation easements on private land. We compared multiple formulations of occupancy models with a logistic regression model to predict the locations of conservation easements based on a spatially explicit social–ecological systems framework. We combined a simulation experiment with a case study of easement data in Idaho and Montana (United States) to illustrate the utility of the occupancy framework for modeling conservation on private land. Occupancy models that explicitly accounted for variation in reporting produced estimates of predictors that were substantially less biased than estimates produced by logistic regression under all simulated conditions. Occupancy models produced estimates for the 6 predictors we evaluated in our case study that were larger in magnitude, but less certain than those produced by logistic regression. These results suggest that occupancy models result in qualitatively different inferences regarding the effects of predictors on conservation easement occurrence than logistic regression and highlight the importance of integrating variable and incomplete reporting of participation in empirical analysis of conservation initiatives. Failure to do so can lead to emphasizing the wrong social, institutional, and environmental factors that enable conservation and underestimating conservation opportunities in landscapes where social norms or institutional constraints inhibit reporting.     

Applying network theory to prioritize multispecies habitat networks that are robust to climate and land‐use change

Designing connected landscapes is among the most widespread strategies for achieving biodiversity conservation targets. The challenge lies in simultaneously satisfying the connectivity needs of multiple species at multiple spatial scales under uncertain climate and land‐use change. To evaluate the contribution of remnant habitat fragments to the connectivity of regional habitat networks, we developed a method to integrate uncertainty in climate and land‐use change projections with the latest developments in network‐connectivity research and spatial, multipurpose conservation prioritization. We used land‐use change simulations to explore robustness of species’ habitat networks to alternative development scenarios. We applied our method to 14 vertebrate focal species of periurban Montreal, Canada. Accounting for connectivity in spatial prioritization strongly modified conservation priorities and the modified priorities were robust to uncertain climate change. Setting conservation priorities based on habitat quality and connectivity maintained a large proportion of the region's connectivity, despite anticipated habitat loss due to climate and land‐use change. The application of connectivity criteria alongside habitat‐quality criteria for protected‐area design was efficient with respect to the amount of area that needs protection and did not necessarily amplify trade‐offs among conservation criteria. Our approach and results are being applied in and around Montreal and are well suited to the design of ecological networks and green infrastructure for the conservation of biodiversity and ecosystem services in other regions, in particular regions around large cities, where connectivity is critically low.     

Predicting occurrence of juvenile shark habitat to improve conservation planning

Fishing and habitat degradation have increased the extinction risk of sharks, and conservation strategies recognize that survival of juveniles is critical for the effective management of shark populations. Despite the rapid expansion of marine protected areas (MPAs) globally, the paucity of shark‐monitoring data on large scales (100s–1000s km) means that the effectiveness of MPAs in halting shark declines remains unclear. Using data collected by baited remote underwater video systems (BRUVS) in northwestern Australia, we developed generalized linear models to elucidate the ecological drivers of habitat suitability for juvenile sharks. We assessed occurrence patterns at the order and species levels. We included all juvenile sharks sampled and the 3 most abundant species sampled separately (grey reef [Carcharhinus amblyrhynchos], sandbar [Carcharhinus plumbeus], and whitetip reef sharks [Triaenodon obesus]). We predicted the occurrence of juvenile sharks across 490,515 km² of coastal waters and quantified the representation of highly suitable habitats within MPAs. Our species‐level models had higher accuracy (? ≥ 0.69) and deviance explained (≥48%) than our order‐level model (? = 0.36 and deviance explained of 10%). Maps of predicted occurrence revealed different species‐specific patterns of highly suitable habitat. These differences likely reflect different physiological or resource requirements between individual species and validate concerns over the utility of conservation targets based on aggregate species groups as opposed to a species‐focused approach. Highly suitable habitats were poorly represented in MPAs with the most restrictions on extractive activities. This spatial mismatch possibly indicates a lack of explicit conservation targets and information on species distribution during the planning process. Non‐extractive BRUVS provided a useful platform for building the suitability models across large scales to assist conservation planning across multiple maritime jurisdictions, and our approach provides a simple for method for testing the effectiveness of MPAs.     

Evaluating the Effects of Anthropogenic Stressors on Source‐Sink Dynamics in Pond‐Breeding Amphibians

Although interwetland dispersal is thought to play an important role in regional persistence of pond‐breeding amphibians, few researchers have modeled amphibian metapopulation or source‐sink dynamics. Results of recent modeling studies suggest anthropogenic stressors, such as pollution, can negatively affect density and population viability of amphibians breeding in isolated wetlands. Presumably population declines also result in reduced dispersal to surrounding (often uncontaminated) habitats, potentially affecting dynamics of nearby populations. We used our data on the effects of mercury (Hg) on the American toad ( Bufo americanus) as a case study in modeling the effects of anthropogenic stressors on landscape‐scale amphibian dynamics. We created a structured metapopulation model to investigate regional dynamics of American toads and to evaluate the degree to which detrimental effects of Hg contamination on individual populations can disrupt interpopulation dynamics. Dispersal from typical American toad populations supported nearby populations that would otherwise have been extirpated over long time scales. Through support of such sink populations, dispersal between wetland‐associated subpopulations substantially increased overall productivity of wetland networks, but this effect declined with increasing interwetland distance and decreasing wetland size. Contamination with Hg substantially reduced productivity of wetland‐associated subpopulations and impaired the ability of populations to support nearby sinks within relevant spatial scales. Our results add to the understanding of regional dynamics of pond‐breeding amphibians, the wide‐reaching negative effects of environmental contaminants, and the potential for restoration or remediation of degraded habitats. Evaluación de los Efectos de Estresantes Antropogénicos sobre la Dinámica Fuente‐Vertedero en Anfibios que se Reproducen en Charcas     

A spatial approach to combatting wildlife crime

Poaching can have devastating impacts on animal and plant numbers, and in many countries has reached crisis levels, with illegal hunters employing increasingly sophisticated techniques. We used data from an 8‐year study in Savé Valley Conservancy, Zimbabwe, to show how geographic profiling—a mathematical technique originally developed in criminology and recently applied to animal foraging and epidemiology—can be adapted for use in investigations of wildlife crime. The data set contained information on over 10,000 incidents of illegal hunting and the deaths of 6,454 wild animals. We used a subset of data for which the illegal hunters’ identities were known. Our model identified the illegal hunters’ home villages based on the spatial locations of the hunting incidences (e.g., snares). Identification of the villages was improved by manipulating the probability surface inside the conservancy to reflect the fact that although the illegal hunters mostly live outside the conservancy, the majority of hunting occurs inside the conservancy (in criminology terms, commuter crime). These results combined with rigorous simulations showed for the first time how geographic profiling can be combined with GIS data and applied to situations with more complex spatial patterns, for example, where landscape heterogeneity means some parts of the study area are less likely to be used (e.g., aquatic areas for terrestrial animals) or where landscape permeability differs (e.g., forest bats tend not to fly over open areas). More broadly, these results show how geographic profiling can be used to target antipoaching interventions more effectively and more efficiently and to develop management strategies and conservation plans in a range of conservation scenarios.     

Living at the Edge: Local versus Positional Factors in the Long‐Term Population Dynamics of an Endangered Orchid

Abstract: Populations at the margin of geographic ranges of distribution have been considered more vulnerable than central ones, but recent reviews have caste doubt on this generalization. We examined the reproductive and demographic performance of a rare Euroasiatic orchid (Cypripedium calceolus) at its southwesterly range limit and compared our findings with those of previous studies of nine central populations at the center of the orchid's range. We sought to test the central‐marginal model and to evaluate factors involved in long‐term performance of forest Eurosiberian species with peripheral populations in southern European mountains. We characterized (structure, temporal fluctuations, herbivory, reproductive success, and recruitment at different habitats) four Pyrenean populations of C. calceolus of different sizes (5–3500 ramets) and monitored three of them for up to 13 years. Two quantitative stochastic models (count data and matrix models) were used to assess population trends and viability and the effect of herbivory. Contrary to expectations, and despite the negative effect of sporadic events of herbivory, the peripheral populations we studied (except the smallest one) performed similarly or better than populations occurring in central part of the species’ range in terms of reproductive success and population growth rates. Landscape changes over the last 50 years suggest that natural reforestation could be involved in the success of this plant at its southern limit. Forest expansion in the mountain regions of southern Europe may provide new opportunities for plants with geographic distributions centered mainly at higher latitudes and give some hope for their recovery in future scenarios dominated by biodiversity loss.     

Integrating intraseasonal grassland dynamics in cross-scale distribution modeling to support waterbird recovery plans

Despite much discussion about the utility of remote sensing for effective conservation, the inclusion of these technologies in species recovery plans remains largely anecdotal. We developed a modeling approach for the integration of local, spatially measured ecosystem functional dynamics into a species distribution modeling (SDM) framework in which other ecologically relevant factors are modeled separately at broad scales. To illustrate the approach, we incorporated intraseasonal water-vegetation dynamics into a cross-scale SDM for the Common Snipe (Gallinago gallinago), which is highly dependent on water and vegetation dynamics. The Common Snipe is an Iberian grassland waterbird characteristic of European agricultural meadows and a member of one of the most threatened bird guilds. The intraseasonal dynamics of water content of vegetation were measured using the standard deviation of the normalized difference water index time series computed from bimonthly images of the Sentinel-2 satellite. The recovery plan for the Common Snipe in Galicia (northwestern Iberian Peninsula) provided an opportunity to apply our modeling framework. Model accuracy in predicting the species’ distribution at a regional scale (resulting from integration of downscaled climate projections with regional habitat–topographic suitability models) was very high (area under the curve [AUC] of 0.981 and Boyce's index of 0.971). Local water-vegetation dynamic models, based exclusively on Sentinel-2 imagery, were good predictors (AUC of 0.849 and Boyce's index of 0.976). The predictive power improved (AUC of 0.92 and Boyce's index of 0.98) when local model predictions were restricted to areas identified by the continental and regional models as priorities for conservation. Our models also performed well (AUC of 0.90 and Boyce's index of 0.93) when projected to updated water-vegetation conditions. Our modeling framework enabled incorporation of key ecosystem processes closely related to water and carbon cycles while accounting for other factors ecologically relevant to endangered grassland waterbirds across different scales, allowed identification of priority areas for conservation, and provided an opportunity for cost-effective recovery planning by monitoring management effectiveness from space.     

Effects of spatial autocorrelation and sampling design on estimates of protected area effectiveness

Estimating the effectiveness of protected areas (PAs) in reducing deforestation is useful to support decisions on whether to invest in better management of areas already protected or to create new ones. Statistical matching is commonly used to assess this effectiveness, but spatial autocorrelation and regional differences in protection effectiveness are frequently overlooked. Using Colombia as a case study, we employed statistical matching to account for confounding factors in park location and accounted for for spatial autocorrelation to determine statistical significance. We compared the performance of different matching procedures—ways of generating matching pairs at different scales—in estimating PA effectiveness. Differences in matching procedures affected covariate similarity between matched pairs (balance) and estimates of PA effectiveness in reducing deforestation. Independent matching yielded the greatest balance. On average 95% of variables in each region were balanced with independent matching, whereas 33% of variables were balanced when using the method that performed worst. The best estimates suggested that average deforestation inside protected areas in Colombia was 40% lower than in matched sites. Protection significantly reduced deforestation, but PA effectiveness differed among regions. Protected areas in Caribe were the most effective, whereas those in Orinoco and Pacific were least effective. Our results demonstrate that accounting for spatial autocorrelation and using independent matching for each subset of data is needed to infer the effectiveness of protection in reducing deforestation. Not accounting for spatial autocorrelation can distort the assessment of protection effectiveness, increasing type I and II errors and inflating effect size. Our method allowed improved estimates of protection effectiveness across scales and under different conditions and can be applied to other regions to effectively assess PA performance.     

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

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