An individual-based model for southern Lake Superior wolves: A tool to explore the effect of human-caused mortality on a landscape of risk

Abstract Gray wolves (Canis lupus) have complex life-histories due, in part, to mating systems that depend on intra-group dominance hierarchies set within an inter-group (pack) social structure linked to philopatric territories. In addition to this spatially oriented social structure, mortality risk associated with interactions with humans varies spatially. We developed an individual-based spatially explicit (IBSE) model for the southern Lake Superior wolf population to better capture the life-history of wolves in a harvest model. Simulated wolves underwent an annual cycle of life-history stage-dependent mate-finding, dispersal, reproduction, and aging on a simulated landscape reflecting spatially explicit state and water boundaries, Indian reservation boundaries and ceded territories, wolf harvest zones, livestock depredation areas, and a spatial mortality risk surface. The latter 3 surfaces were linked to mortality events for simulated wolves. We assessed our IBSE model and conducted a sensitivity analysis of the most uncertain parameters with a categorical calibration of patterns observed at the individual, pack, population, and landscape level. We found that without recreational harvest, the Wisconsin wolf population grew to an average carrying capacity of 1242 wolves after 50 years and breeding pairs persisted for a mean 1.8 years. We simulated 6 recreational harvest scenarios with varying rates and timings of harvest and assessed effects on population size, pack sizes, age ratios, dispersal and immigration rates, and breeding pair tenures of the Wisconsin wolf population. The simulated harvest with rates of 14% which corresponded to the 2012 harvest in Wisconsin reduced the populations 4% in the first year of harvest and equilibrated to the pre-harvest population size after 20 years of harvest, on average. A 30% harvest rate across the simulation on average reduced the populations by 65% after 20 years with some populations going extinct before 100 years. In general, harvest increased the proportion of pups in the simulated populations and decreased breeding pair tenure. Targeted lethal control was more effective than harvest for reducing the number of wolves near known livestock depredation sites. Our model facilitates prediction of important population patterns that is simultaneously dependent on complexities associated with spatially structured life history and mortality.

[1]  M. Gail,et al.  Likelihood calculations for matched case-control studies and survival studies with tied death times , 1981 .

[2]  David J. Mladenoff,et al.  Predicting Human‐Carnivore Conflict: a Spatial Model Derived from 25 Years of Data on Wolf Predation on Livestock , 2004 .

[3]  R. Levins.  The strategy of model building in population biology , 1966 .

[4]  Michael J. North,et al.  Tutorial on agent-based modelling and simulation , 2005, Proceedings of the Winter Simulation Conference, 2005..

[5]  B. White,et al.  Protection from harvesting restores the natural social structure of eastern wolf packs , 2010 .

[6]  J. Gareth Polhill,et al.  The ODD protocol: A review and first update , 2010, Ecological Modelling.

[7]  L. Boitani,et al.  Wolf Social Ecology , 2003 .

[8]  T. Fuller,et al.  Wolf Population Dynamics , 2003 .

[9]  R. Wayne,et al.  The genealogy and genetic viability of reintroduced Yellowstone grey wolves , 2008, Molecular ecology.

[10]  R. Haight,et al.  A Regional Landscape Analysis and Prediction of Favorable Gray Wolf Habitat in the Northern Great Lakes Region , 1995 .

[11]  Frederick F. Knowlton,et al.  Evaluating coyote management strategies using a spatially explicit, individual-based, socially structured population model , 2008 .

[12]  M. Hebblewhite,et al.  A spatially explicit model for an Allee effect: why wolves recolonize so slowly in Greater Yellowstone. , 2006, Theoretical population biology.

[13]  Birgit Müller,et al.  A standard protocol for describing individual-based and agent-based models , 2006 .

[14]  M. Clayton,et al.  Change in Occupied Wolf Habitat in the Northern Great Lakes Region , 2009 .

[15]  O. Liberg,et al.  The Effects of Breeder Loss on Wolves , 2008 .

[16]  Adrian P. Wydeven,et al.  Recovery of Gray Wolves in the Great Lakes Region of the United States: An Endangered Species Success Story , 2009 .

[17]  Winfried Kurth,et al.  Facilitating Parameter Estimation and Sensitivity Analysis of Agent-Based Models: A Cookbook Using NetLogo and 'R' , 2014, J. Artif. Soc. Soc. Simul..

[18]  J. Bart.  Acceptance Criteria for Using Individual‐Based Models to Make Management Decisions , 1995 .

[19]  Luigi Boitani,et al.  Wolves : behavior, ecology, and conservation , 2003 .

[20]  Steven F. Railsback,et al.  Individual-based modeling and ecology , 2005 .

[21]  L. Waits,et al.  Estimating gray wolf pack size and family relationships using noninvasive genetic sampling at rendezvous sites , 2011 .

[22]  Jean-Dominique Lebreton,et al.  DYNAMICAL AND STATISTICAL MODELS FOR EXPLOITED POPULATIONS † , 2005 .

[23]  R. Haight,et al.  Computer simulation of vasectomy for wolf control , 1997 .

[24]  Charles M. Macal,et al.  Tutorial on agent-based modeling and simulation , 2005 .

[25]  H. Pulliam,et al.  Sources, Sinks, and Habitat Selection: A Landscape Perspective on Population Dynamics , 1991, The American Naturalist.

[26]  L. E. Travis,et al.  Computer simulation of wolf-removal strategies for animal damage control , 2002 .

[27]  Adrian P. Wydeven,et al.  Wolf–Human Conflicts and Management in Minnesota, Wisconsin, and Michigan , 2009 .

[28]  P. David.  Ma’iingan and the Ojibwe , 2009 .

[29]  D. Boukal,et al.  Linking the Allee Effect, Sexual Reproduction, and Temperature‐Dependent Sex Determination Via Spatial Dynamics , 2001, The American Naturalist.

[30]  Adrian P. Wydeven,et al.  Dispersal of Gray Wolves in the Great Lakes Region , 2009 .

[31]  Volker Grimm,et al.  Merging validation and evaluation of ecological models to ‘evaludation’: A review of terminology and a practical approach , 2014 .

[32]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[33]  David J. Mladenoff,et al.  Predicting gray wolf landscape recolonization: logistic regression models vs. new field data , 1999 .

[34]  Carl J. Walters,et al.  ECOLOGICAL OPTIMIZATION AND ADAPTIVE MANAGEMENT , 1978 .

[35]  V. Grimm.  Ten years of individual-based modelling in ecology: what have we learned and what could we learn in the future? , 1999 .

[36]  François Bousquet,et al.  Multi-agent simulations and ecosystem management: a review , 2004 .

[37]  Paul W. Box,et al.  An individual-based model of canid populations: modelling territoriality and social structure , 2003 .

[38]  Uta Berger,et al.  Pattern-Oriented Modeling of Agent-Based Complex Systems: Lessons from Ecology , 2005, Science.

[39]  J. Swenson,et al.  Saving large carnivores, but losing the apex predator? , 2013 .

[40]  Adrian P. Wydeven,et al.  History, Population Growth, and Management of Wolves in Wisconsin , 2009 .

[41]  Nicholas J. Gotelli,et al.  A Primer of Ecology , 1995 .

[42]  Steven F Railsback,et al.  Pattern-oriented modelling: a ‘multi-scope’ for predictive systems ecology , 2012, Philosophical Transactions of the Royal Society B: Biological Sciences.