Multiscale habitat relationships of snowshoe hares (Lepus americanus) in the mixed conifer landscape of the Northern Rockies, USA: Cross-scale effects of horizontal cover with implications for forest management

Abstract

Snowshoe hares (Lepus americanus) are an ecologically important herbivore because they modify vegetation through browsing and serve as a prey resource for multiple predators. We implemented a multiscale approach to characterize habitat relationships for snowshoe hares across the mixed conifer landscape of the northern Rocky Mountains, USA. Our objectives were to (1) assess the relationship between horizontal cover and snowshoe hares, (2) estimate how forest metrics vary across the gradient of snowshoe hare use and horizontal cover, and (3) model and map snowshoe hare occupancy and intensity of use. Results indicated that both occupancy and intensity of use by snowshoe hares increased with horizontal cover and that the effect became stronger as intensity of use increased. This underscores the importance of dense horizontal cover to achieve high use, and likely density, of snowshoe hares. Forest structure in areas with high snowshoe hare use and horizontal cover was characterized as multistoried with dense canopy cover and medium-sized trees (e.g., 12.7–24.4 cm). The abundance of lodgepole pine (Pinus contorta) was associated with snowshoe hare use within a mixed conifer context, and the only species to increase in abundance with horizontal cover was Engelmann spruce (Picea engelmannii) and subalpine fir (Abies lasiocarpa). Our landscape-level modeling produced similar patterns in that we observed a positive effect of lodgepole pine and horizontal cover on both occupancy and use by snowshoe hares, but we also observed a positive yet parabolic effect of snow depth on snowshoe hare occupancy. This work is among the first to characterize the multiscale habitat relationships of snowshoe hares across a mixed conifer landscape as well as to map their occupancy and intensity of use. Moreover, our results provide stand- and landscape-level insights that directly relate to management agencies, which aids in conservation efforts of snowshoe hares and their associated predators.

References Abele, S. L., Wirsing, A. J., & Murray, D. L. (2013). Precommercial forest thinning alters abundance but not survival of snowshoe hares. Journal of Wildlife Management, 77, 84–92. Allard-Duchêne, A., Pothier, D., Dupuch, A., & Fortin, D. (2014). Temporal changes in habitat use by snowshoe hares and red squirrels during post-fire and post-logging forest succession. Forest Ecology and Management, 313, 17–25. Ausband, D. E., & Baty, G. R. (2005). Effects of precommercial thinning on snowshoe hare habitat use during winter in low-elevation montane forests. Canadian Journal of Forestry Research, 35, 206–210. Barnett, T. P., Adam, J. C., & Lettenmaier, D. P. (2005). Potential impacts of a warming climate on water availability in snow-dominated regions. Nature, 438, 303–309. Berg, N. D., & Gese, E. M. (2010). Relationship between fecal pellet counts and snowshoe hare density in western Wyoming. Journal of Wildlife Management, 74, 1745–1751. Berg, N. D., Gese, E. M., Squires, J. R., & Aubry, L. M. (2012). Influence of forest structure on the abundance of snowshoe hares in western Wyoming. Journal of Wildlife Management, 76, 1480–1488. Boulangeat, I., Gravel, D., & Thuiller, W. (2012). Accounting for dispersal and biotic interactions to disentangle the drivers of species distributions and their abundances. Ecology Letters, 15, 584–593. Breiman, L. (2001). Random Forests. Machine Learning, 45, 5–32. Brennan, A., Cross, P. C., & Creel, S. (2015). Managing more than the mean: Using quantile regression to identify factors related to large elk groups. Journal of Applied Ecology, 52, 1656–1664. Bunn, A. G. (2009). The rock and ice problem in national parks: An opportunity for monitoring climate change impacts. Park Science, 26, 17–21. Cade, B. S., & Noon, B. R. (2003). A gentle introduction to quantile regression for ecologists. Frontiers in Ecology and the Environment, 1, 412–420. Cheng, E., Hodges, K. E., & Mills, L. S. (2015). Impacts of fire on snowshoe hares in Glacier National Park Montana, USA. Fire Ecology, 11, 119–136. Clow, D. W., Nanus, L., Verdin, K. L., & Schmidt, J. (2012). Evaluation of SNODAS snow depth and snow water equivalent estimates for the Colorado Rocky Mountains, USA. Hydrological Processes, 26, 2583–2591. Crookston, N. L., & Dixon, G. E. (2005). The forest vegetation simulator: A review of its structure, content, and applications. Computers and Electronics in Agriculture, 49, 60–80. Cushman, S. A., Gutzweiler, K., Evans, J. S., & McGarigal, K. (2010). The gradient paradigm: A conceptual and analytical framework for landscape ecology. In S. Cushman, & F. Huettmann (Eds.), Spatial complexity, informatics, and wildlife conservation (pp. 83–108). New York, NY, USA: Springer. Cutler, D. R., Edwards Jr, T. C., Beard, K. H., Cutler, A., Hess, K. T., Gibson, J., & Lawler, J. J. (2007). Random forests for classification in ecology. Ecology, 88, 2783–2792. Ellsworth, E., Wirsing, A. J., Shipley, L. A., & Murray, D. L. (2013). Do measures of plant intake and digestibility from captive feeding trials align with foraging patterns of free-ranging snowshoe hares. Wildlife Research, 40, 349–357. Elton, C., & Nicholson, M. (1942). The ten-year cycle in numbers of the lynx in Canada. Journal of Animal Ecology, 11, 215–244. ESRI (Environmental Systems Research Institute). (2011). ArcGIS desktop: Release 10. Redlands, CA, USA: ESRI. Evans, J. S., & Cushman, S. A. (2009). Gradient modeling of conifer species using random forests. Landscape Ecology, 24, 673–683. Evans, J. S., Oakleaf, J., Cushman, S. A., & Theobald, D. (2014). An ArcGIS toolbox for surface gradient and geomorphometric modeling, version2.0-0. Retrieved from http://evansmurphy.wix.com/evansspatial. Feierabend, D., & Kielland, K. (2015). Seasonal effects of habitat on sources and rates of snowshoe hare predation in Alaskan boreal forests. PLoS One, 10, e0143543. Fiala, A. C. S., Garman, S. L., & Gray, A. N. (2006). Comparison of five canopy cover estimation techniques in the western Oregon cascades. Forest Ecology and Management, 232, 188–197. Fuller, A. K., & Harrison, D. J. (2013). Modeling the influence of forest structure on microsite habitat use by snowshoe hares. International Journal of Forestry Research, 2013, Article ID 892327. Gessler, P. E., Moore, I. D., McKenzie, N. J., & Ryan, P. J. (1995). Soil-landscape modelling and spatial prediction of soil attributes. International Journal of Geographical Information Systems, 9, 421–432. Griffin, P. C., & Mills, L. S. (2007). Precommercial thinning reduces snowshoe hare abundance in the short term. Journal of Wildlife Management, 71, 559–564. Griffin, P. C., & Mills, L. S. (2009). Sinks without borders: Snowshoe hare dynamics in a complex landscape. Oikos, 118, 1487–1498. Grinnell, J. (1917). Field tests of theories concerning distributional control. The American Naturalist, 51, 115–128. Guisan, A., Weiss, S. B., & Weiss, A. D. (1999). GLM versus CCA spatial modeling of plant species distribution. Plant Ecology, 143, 107–122. Hijmans, R. (2015). Raster: Geographic data analysis and modeling. R package version 2.5-2. Retrieved from https://CRAN.R-project.org/package=raster. Hirzel, A. H., & Le Lay, G. (2008). Habitat suitability modelling and niche theory. Journal of Applied Ecology, 45, 1372–1381. Hodges, K. E. (2000a). The ecology of snowshoe hares in northern boreal forests. In L. F. Ruggiero, K. B. Aubry, S. W. Buskirk, G. M. Koehler, C. J. Krebs, K. S. McKelvey, & J. R. Squires (Eds.), Ecology and conservation of lynx in the United States (pp. 117–161). Boulder, CO, USA: University Press of Colorado. Hodges, K. E. (2000b). Ecology of snowshoe hares in southern boreal and montane forests. In L. F. Ruggiero, K. B. Aubry, S. W. Buskirk, G. M. Koehler, C. J. Krebs, K. S. McKelvey, & J. R. Squires (Eds.), Ecology and conservation of lynx in the United States (pp. 163–206). Boulder, CO, USA: University Press of Colorado. Hodges, K. E., Mills, L. S., & Murphy, K. M. (2009). Distribution and abundance of snowshoe hares in Yellowstone National Park. Journal of Mammalogy, 90, 870–878. Hodges, K. E., & Sinclair, A. R. E. (2003). Does predation risk cause snowshoe hares to modify their diets? Canadian Journal of Zoology, 81, 1973–1985. Hodges, K. E., & Sinclair, A. R. E. (2005). Brows site selection by snowshoe hares: Effects of food supply and predation risk. Canadian Journal of Zoology, 83, 280–292. Hodson, J., Fortin, D., & Bélanger, L. (2010). Fine-scale disturbances shape space-use patterns of a boreal forest herbivore. Journal of Mammalogy, 91, 607–619. Hodson, J., Fortin, D., & Bélanger, L. (2011). Changes in relative abundance of snowshoe hares (Lepus americanus) across a 265-year gradient of boreal forest succession. Canadian Journal of Zoology, 89, 908–920. Holbrook, J. D., Arkle, R. S., Rachlow, J. L., Vierling, K. T., Pilliod, D. S., & Wiest, M. W. (2016). Occupancy and abundance of predator and prey: Implications of the fire-cheatgrass cycle in sagebrush ecosystems. Ecosphere, 7, e01307. Homyack, J. A., Harrison, D. J., & Krohn, W. B. (2007). Effects of precommercial thinning on snowshoe hares in Maine. Journal of Wildlife Management, 71, 4–13. Hope, R. M. 2013. Rmisc: Ryan miscellaneous. R package version 1.5. Retrieved from https://CRAN.R-project.org/package=Rmisc. Huang, C., Wylie, B., Yang, L., Homer, C., & Zylstra, G. (2002). Derivation of a tasseled cap transformation based on Landsat 7 at-satellite reflectance. International Journal of Remote Sensing, 23, 1741–1748. Ivan, J. S., & Shenk, T. M. (2016). Winter diet and hunting success of Canada lynx in Colorado. Journal of Wildlife Management, 80, 1049–1058. Ivan, J. S., White, G. C., & Shenk, T. M. (2014). Density and demography of snowshoe hares in central Colorado. Journal of Wildlife Management, 78, 580–594. Jenness, J. S. (2004). Calculating landscape surface area from digital elevation models. Wildlife Society Bulletin, 32, 829–839. Jenness, J. (2013). DEM surface tools for ArcGIS (surface_area.exe). Jenness Enterprises. Retrieved from http://www.jennessent.com/arcgis/surface_area.htm. Klos, P. Z., Link, T. E., & Abatzoglou, J. T. (2014). Extent of the rain-snow transition zone in the western U.S. under historic and projected climate. Geophysical Research Letters, 41, 4560–4568. Koehler, G. M. (1990). Population and habitat characteristics of lynx and snowshoe hares in north central Washington. Canadian Journal of Zoology, 68, 845–851. Koenker, R. (2016). quantreg: Quantile regression. R package version 5.21. Retrieved from https://CRAN.R-project.org/package=quantreg. Koenker, R., & Bassett, G. (1978). Regression quantiles. Econometrica, 46, 33–50. Krebs, C. J. (2011). Of lemmings and snowshoe hares: The ecology of northern Canada. Proceedings of the Royal Society of London. Series B, Biological Sciences, 278, 481–489. Krebs, C. J., Boonstra, R., Nams, V., O'Donoghue, M., Hodges, K. E., & Boutin, S. (2001). Estimating snowshoe hare population density from pellet counts: A further evaluation. Canadian Journal of Zoology, 79, 14. Krebs, C. J., Boutin, S., & Boonstra, R. (2001). Ecosystem dynamics of the boreal forest: The Kluane project. New York, NY, USA: Oxford University Press. Krebs, C. J., Gilbert, B. S., Boutin, S., & Boonstra, R. (1987). Estimation of snowshoe hare population density from turn transects. Canadian Journal of Zoology, 65, 565–567. Kuhn, M., Wing, J., Weston, S., Williams, A., Keefer, C., Engelhardt, A., … Canda, C. (2016). Caret: Classification and regression training. R package version 6.0-64. Retrieved from https://CRAN.R-project.org/packages=caret. LANDFIRE. (2012). Existing vegetation cover layer. U.S. Department of Agriculture and U.S. Department of the Interior. Retrieved from http://landfire.cr.usgs.gov/viewer/. Lewis, C. W., Hodges, K. E., Koehler, G. M., & Mills, L. S. (2011). Influence of stand and landscape features on snowshoe hare abundance in fragmented forests. Journal of Mammalogy, 92, 561–567. Liaw, A., & Wiener, M. (2002). Classification and regression by randomForest. R News, 2, 18–22. Lindsey, A. A., Barton Jr, J. D., & Miles, S. R. (1958). Field efficiencies of forest sampling methods. Ecology, 39, 428–444. Litvaitis, J. A., Sherburne, J. A., & Bissonette, J. A. (1985). Influence of understory characteristics on snowshoe hare habitat use and density. Journal of Wildlife Management, 49, 866–873. Liu, Y., Goodrick, S. L., & Stanturf, J. A. (2013). Future U.S. wildfire potential trends projected using a dynamically downscaled climate change scenario. Forest Ecology and Management, 294, 120–135. MacKenzie, D., Nichols, J. D., Lachman, G. B., Droege, S., Royle, J. A., & Lantimm, C. A. (2002). Estimating site occupancy rates when detection probabilities are less than one. Ecology, 83, 2248–2255. McCann, K., Hastings, A., & Huxel, G. R. (1998). Weak trophic interactions and the balance of nature. Nature, 395, 794–798. McCune, B., & Keon, D. (2002). Equations for potential annual direct incident radiation and heat load. Journal of Vegetation Science, 13, 603–606. McGarigal, K., Tagil, S., & Cushman, S. A. (2009). Surface metrics: An alternative to patch metrics for the quantification of landscape structure. Landscape Ecology, 24, 433–450. McGill, B. J. (2010). Matters of scale. Science, 328, 575–576. McKelvey, K. S., Copeland, J. P., Schwartz, M. K., Littell, J. S., Aubry, K. B., Squires, J. R., … Mauger, G. S. (2011). Climate change predicted to shift wolverine distributions, connectivity, and dispersal corridors. Ecological Applications, 21, 2882–2897. Milburn, A., Bollenbacker, B., Manning, M., & Bush, R. (2015). Region 1 existing and potential vegetation groupings used for broad-level analysis and monitoring. Report 15-4 v1.0. Mills, L. S., Griffin, P. C., Hodges, K. E., McKelvey, K., Ruggiero, L., & Ulizio, T. (2005). Pellet count indices compared to mark-recapture estimates for evaluating snowshoe hare density. Journal of Wildlife Management, 69, 1053–1062. Mills, L. S., Zimova, M., Oyler, J., Running, S., Abatzoglou, J. T., & Lukacs, P. M. (2013). Camouflage mismatch in seasonal coat color due to decreased snow duration. Proceedings of the National Academy of Sciences of the United States of America, 110, 7099–7528. Murphy, M. A., Evans, J. S., & Storfer, A. (2010). Quantifying Bufo boreas connectivity in Yellowstone National Park with landscape genetics. Ecology, 91, 252–261. Murray, D., Ellsworth, E., & Zack, A. (2005). Assessment of potential bias with snowshoe hare fecal pellet-plot counts. Journal of Wildlife Management, 69, 385–395. NOHRSC (National Operational Hydrologic Remote Sensing Center). (2004). Snow data assimilation system (SNODAS) data products at NSIDC, 1 April 2012–1 April 2013. Boulder, CO, USA: National Snow and Ice Data Center. doi:10.7265/N5TB14TC Pettorelli, N., Olav Vik, J., Mysterud, A., Gillard, J. M., Tucker, C. J., & Stenseth, N. C. (2005). Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends in Ecology and Evolution, 20, 503–510. Picotte, J. J., Peterson, B., Meier, G., & Howard, S. M. (2016). 1984–2010 trends in fire burn severity and area for the conterminous US. International Journal of Wildland Fire, 25, 413–420. Pietz, P. J., & Tester, J. R. (1983). Habitat selection by snowshoe hares in north central Minnesota. Journal of Wildlife Management, 47, 686–696. Poole, K. G. (2003). A review of Canada Lynx, Lynx canadensis, in Canada. The Canadian Field-Naturalist, 117, 360–376. Prugh, L. R., & Krebs, C. J. (2004). Snowshoe hare pellet-decay rates and aging in different habitats. Wildlife Society Bulletin, 32, 386–393. R Core Team (2015). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://R-project.org/. Roth, J. D., Marshall, J. D., Murray, D. L., Nickerson, D. M., & Steury, T. D. (2007). Geographical gradients in diet affect population dynamics of Canada lynx. Ecology, 88, 2736–2743. Savage, S. L., Lawrence, R. L., & Squires, J. R. (2015). Predicting relative species composition within mixed conifer forests pixels using zero-inflated models and Landsat imagery. Remote Sensing of Environment, 171, 326–336. Squires, J. R., DeCesare, N. J., Kolbe, J. A., & Ruggiero, L. F. (2010). Seasonal resource selection of Canada lynx in managed forests of the northern Rocky Mountains. Journal of Wildlife Management, 74, 1648–1660. Squires, J. R., DeCesare, N. J., Olson, L. E., Kolbe, J. A., Hebblewhite, M., & Parks, S. A. (2013). Combining resource selection and movement behavior to predict corridors for Canada lynx at their southern range periphery. Biological Conservation, 157, 187–195. Squires, J. R., & Ruggiero, L. F. (2007). Winter prey selection of Canada lynx in northwestern Montana. Journal of Wildlife Management, 71, 310–315. Stavros, E. N., Abatzoglou, J. T., McKenzie, D., & Larkin, N. K. (2014). Regional projections of the likelihood of very large wildland fires under a changing climate in the contiguous western United States. Climatic Change, 126, 455–468. Sultaire, S. M., Pauli, J. N., Martin, K. J., Meyer, M. W., Notaro, M., & Zuckerberg, B. (2016). Climate change surpasses land-use change in the contracting range boundary of a winter-adapted mammal. Proceedings of the Royal Society of London. Series B, Biological Sciences, 283, 20153104. Sultaire, S. M., Pauli, J. N., Martin, K. J., Meyer, M. W., & Zuckerberg, B. (2016). Extensive forests and persistent snow cover promote snowshoe hare occupancy in Wisconsin. Journal of Wildlife Management, 80, 894–905. Thornton, D. H., Wirsing, A. J., Roth, J. D., & Murray, D. L. (2013). Habitat quality and population density drive occupancy dynamics of snowshoe hare in variegated landscapes. Ecography, 36, 610–621. USFWS (United States Fish and Wildlife Service) (2000). Determination of threatened status for the contiguous US distinct population segment of the Canada lynx and related rule: final rule. US Federal Register, 65, 16052–16086. USFWS (United States Fish and Wildlife Service) (2013). Revised designation of critical habitat for the contiguous US distinct population segment of the Canada lynx and revised distinct population segment boundary. US Federal Register, 78, 59429–59474. Wickham, H. (2009). ggplot2: Elegant graphics for data analysis. New York, NY, USA: Springer-Verlag. Wirsing, A. J., Steury, T. D., & Murray, D. L. (2002). A demographic analysis of a southern snowshoe hare population in a fragmented habitat: evaluating the refugium model. Canadian Journal of Zoology, 80, 169–177. Zimova, M., Mills, L. S., & Nowak, J. J. (2016). High fitness costs of climate change-induced camouflage mismatch. Ecology Letters, 19, 299–307.