Wildlife population dynamics are modulated by abiotic and biotic factors, typically climate, resource availability, density-dependent effects, and predator-prey interactions. Understanding if human-caused disturbances shape these processes is needed for the conservation and management of ecological communities within increasingly human-dominated landscapes. Garnering this understanding is difficult due to the lack of long-term longitudinal data on wildlife populations. Concurrently, the interplay between long-term human-mediated disturbances, climate, and predator density on ungulate population dynamics has been under-studied. Using a 50-year time series (1962-2012) on mule deer (Odocoileus hemionus) demographics, seasonal weather, predator density, and oil and gas development patterns from the North Dakota Badlands to investigate the long-term effects of landscape-level disturbance, we aimed to evaluate if harsh weather conditions in-combination with energy development and predators affected fall mule deer recruitment. We found that density-dependent effects and harsh seasonal weather primarily drove recruitment in the North Dakota Badlands. Recruitment was further shaped by interacting effects of harsh seasonal weather and predator presence in the form of high coyote density. Additionally, we found that fall recruitment was subtly modulated by interactions between seasonal weather and energy development (i.e., lower recruitment when harsher weather was combined with higher density of active oil and gas wells), and that the combined effect of predator density and energy development was not interactive but rather additive. Our analysis demonstrates the effect of energy development by modulating mule deer recruitment fluctuations concurrent with main recruitment drivers being biotic (density-dependency, habitat, predation) and abiotic (harsh seasonal weather, woody vegetation encroachment). A pattern emerges of density dependence, presumably due to limited quality habitat, being the primary factor influencing fall fawn recruitment in mule deer. Secondarily, stochastic weather events periodically cause dramatic declines in recruitment. Finally, the interactions between human disturbance and predation can be additive to the aforementioned drivers of recruitment and subsequently cause further declines.

Understanding what variables affect ungulate neonate survival is imperative to successful conservation and management of the species. Predation is commonly cited as a cause-specific source of mortality and ecological covariates often influence neonate survival. However, variation in survival estimates related to capture methodology has been documented with opportunistically captured neonates generally displaying greater survival than those captured via aid of vaginal implant transmitters (VITs), likely because of increased left truncation observed in the opportunistically captured datasets. Our goal was to assess if 3- and 6-month survival estimates varied by capture method while simultaneously assessing if capture method affected model selection and interpretation of ecological covariates for white-tailed deer neonates captured from three study sites in North Dakota and South Dakota, USA. We found survival varied by capture method for 3-month neonate survival with opportunistically captured neonates displaying up to 26% greater survival than their counterparts captured via VITs; however, this relationship was not present for 6-month survival. We also found model selection and subsequent interpretation of ecological covariates varied when analyzing datasets comprised of neonates captured via VITs, neonates captured opportunistically, and all neonates combined regardless of capture method. When interpreting results from our VIT only analysis for 3-month survival, we found survival varied by three time intervals and was lowest in the first two weeks of life. Capture method did not affect 6-month survival which was most influenced by total precipitation occurring during 3 – 8 weeks of a neonate’s life and percent canopy cover found at a neonate’s capture site. Our results support previous research that capture method must be accounted for when deriving survival estimates for ungulate neonates as it can impact derived estimates and subsequent interpretation of results.