Using integrated step selection to determine effects of predation risk on bison habitat selection and movement

Animal movement is a fundamental mechanism that shapes communities and ecosystems. Ungulates alter the ecosystems they inhabit and understanding their movements and distribution is critical for linking habitat with population dynamics. Predation risk has been shown to strongly influence ungulate movement patterns, such that ungulates may select habitat where predation risk is lower (refugia), adjust movement rates, temporal patterns, or selection of cover variables in areas with greater predation risk. We evaluated potential predation avoidance behavior in a population of plains bison inhabiting the north rim of Grand Canyon National Park (GRCA) and adjacent Kaibab National Forest (KNF). The KNF has year-round hunting managed by Arizona Game and Fish Department, whereas hunting is not allowed in GRCA. Human-maintained water sources on the KNF are particularly important resources for bison wherein they may be exposed to higher predation risk to access these resources. We used 2-h GPS locations for three years from 31 bison (n = 9 males; n = 22 females), and integrative step selection analysis to test four hypotheses about the potential for bison to reduce their risk from human predation by avoiding areas of high predation risk; moving faster in areas with high predation risk; entering high-risk areas at night when risk is reduced; and entering high-risk areas in habitats that provide cover (coniferous forest). The highest performing model indicated bison movement was 1.3 times faster per 2-h step interval than in areas with no hunting across all vegetation classes (coniferous forest, shrub, quaking aspen, grass-forb meadow) and across all topography classes (valley, slope, ridge). Bison moved more slowly in grass-forb meadows than all other vegetation types, and in valleys relative to slopes and ridges. Several radio-collared individuals had no GPS locations in KNF for the duration of the study. Bison avoided predation risk using two strategies: moving faster while in the KNF, and fully avoiding high-risk areas by remaining within GRCA. Management that manipulates or reduces timing of hunting seasons may reduce perceived predation risk and encourage bison to distribute into the KNF and across a broader range of available habitat.