Cumulative effects of human footprint, natural features and predation risk best predict seasonal resource selection by white-tailed deer

Land modified for human use alters matrix shape and composition and is a leading contributor to global biodiversity loss. It can also play a key role in facilitating range expansion and ecosystem invasion by anthrophilic species, as it can alter food abundance and distribution while also influencing predation risk; the relative roles of these processes are key to habitat selection theory. We researched these relative influences by examining human footprint, natural habitat, and predator occurrence on seasonal habitat selection by range-expanding boreal white-tailed deer (Odocoileus virginianus) in the oil sands of western Canada. We hypothesized that polygonal industrial features (e.g. cutblocks, well sites) drive deer distributions as sources of early seral forage, while linear features (e.g. roads, trails, and seismic lines) and habitat associated with predators are avoided by deer. We developed seasonal 2nd -order resource selection models from three years of deer GPS-telemetry data, a camera-trap-based model of predator occurrence, and landscape spatial data to weigh evidence for six competing hypotheses. Deer habitat selection was best explained by the combination of polygonal and linear features, intact deciduous forest, and wolf (Canis lupus) occurrence. Deer strongly selected for linear features such as roads and trails, despite a potential increased risk of wolf encounters. Linear features may attract deer by providing high density forage opportunity in heavily exploited landscapes, facilitating expansion into the boreal north.