Projecting Mammal Distributions in Response to Future Alternative Landscapes in a Rapidly Transitioning Region

: Finding balance between the needs of people and wildlife is an essential component of planning sustainable landscapes. Because mammals make up a diverse and ecologically important taxon with varying responses to human disturbance, we used representative mammal species to examine how alternative land-use policies might a ﬀ ect their habitats and distributions in the near future. We used wildlife detections from camera traps at 1591 locations along a large-scale urban to wild gradient in northern Virginia, to create occupancy models which determined land cover relationships and the drivers of contemporary mammal distributions. From the 15 species detected, we classiﬁed ﬁve representative species into two groups based on their responses to human development; sensitive species (American black bears and bobcats) and synanthropic species (red foxes, domestic cats, and white-tailed deer). We then used the habitat models for the representative species to predict their distributions under four future planning scenarios based on strategic versus reactive planning and high or low human population growth. The distributions of sensitive species did not shrink drastically under any scenario, whereas the distributions of synanthropic species increased in response to anthropogenic development, but the magnitude of the response varied based on the projected rate of human population growth. This is likely because most sensitive species are dependent on large, protected public lands in the region, and the majority of projected habitat losses should occur in non-protected private lands. These ﬁndings illustrate the importance of public protected lands in mitigating range loss due to land use changes, and the potential positive impact of strategic planning in further mitigating mammalian diversity loss in private lands.


Introduction
Effective land planning must balance the needs of people and the natural world [1], particularly because land cover change and habitat loss are among the leading causes of biodiversity loss [2]. Mammals are no exception to the global phenomena of range reduction and habitat loss [3]. Mammals across a range of body sizes and trophic guilds change their behaviors in direct response to human-induced changes in land cover and increased human presence [4,5]. Furthermore, there are apparent dichotomies in mammal responses to changes in the Anthropocene, with species such as omnivorous mesopredators, herbivores, and domesticated species exhibiting clear advantages over obligate carnivores, apex predators, and habitat specialists [6][7][8][9]. For example, some species within carnivore guilds avoid urban areas, whereas others increase their abundances in these anthropogenic environments [10][11][12][13].
Mammals are often difficult to detect even when present, but technological advances in non-invasive sampling (e.g., camera traps, genetic tagging, etc.) have helped accumulate global The eMammal database is a data management platform for camera trap projects with data from both citizen scientists and professionals [24]. We used data from several projects within eMammal in which camera deployments and the data they collected were standardized in camera trap height (~50 cm above ground) and without bait or lures [23,25]. Although these projects differed in their original objectives, all examined wildlife-habitat associations and used passive-infrared motion-activated camera traps (~90% Reconyx RC55, RC800, orc RC900, RECONYX, Inc., Holmen WI, USA or ~10% Bushnell TrophyCam HD, Bushnell, Overland Park, KS, USA) [23,25]. These camera trap data came from sites within the CLI study area, and the greater Washington DC region as well, to capture a fuller range of human development (Figure 1). In total, these data were derived from 1591 camera deployments conducted during the growing seasons from August 2012 through October 2018. Following each deployment, we extracted the raw detection data for each species and created daily detection histories for each site with a minimum of one week of sampling and we truncated the detection histories at 30 days to realize closure for the sites. Although the data were collected over several field seasons, we did not resample any sites in a multi-season framework, and we, therefore, refer to ψ as "site use" [26]; however, we did test for annual and seasonal variation in detection, as described later in the methods.

Covariates and Models
To measure habitat variables, we used ArcGIS 10.7 (Environmental Systems Research Institute [ESRI], Inc., Redlands, CA, USA) to create a land cover map of the camera trapping area based on the reclassification of the 2011 National Land Cover Database (NLCD-Appendix Table A1) with 30 m resolution land cover data modified on the Anderson Level II classification system, from previous CLI research iterations [22,27]. Reclassification was based on generalized landscape characteristics, Following each deployment, we extracted the raw detection data for each species and created daily detection histories for each site with a minimum of one week of sampling and we truncated the detection histories at 30 days to realize closure for the sites. Although the data were collected over several field seasons, we did not resample any sites in a multi-season framework, and we, therefore, refer to ψ as "site use" [26]; however, we did test for annual and seasonal variation in detection, as described later in the methods.

Covariates and Models
To measure habitat variables, we used ArcGIS 10.7 (Environmental Systems Research Institute [ESRI], Inc., Redlands, CA, USA) to create a land cover map of the camera trapping area based on the reclassification of the 2011 National Land Cover Database (NLCD-Appendix A) with 30 m resolution land cover data modified on the Anderson Level II classification system, from previous CLI research iterations [22,27]. Reclassification was based on generalized landscape characteristics, Remote Sens. 2019, 11, 2482 4 of 16 and thus we do not expect these simplifications of land cover to dramatically affect inferences for mammal responses. We selected habitat covariates based on the known ecology of potential candidate mammals [7,8,10,11,23,28].
We measured the coverage percentages of the four main land covers: (1) forest, (2) grassland, (3) crops, and (4) urban development in 500 and 1000 m radii around each 30 x 30 m pixel of land in the study area. We further measured linear distances (m) to a variety of anthropogenic or natural edges: (5) distance to development, (6) distance to roads, and (7) distance to crops as indices of human presence and/or disturbance; (8) distance to grassland, (9) distance to core forest (forest greater than 60 m from other land class), and (10) distance to protected lands, which included private and public lands protected from future development. We also used a binary covariate (11) to denote sites as within or outside of those protected areas (i.e., lands managed by the National Park Service, National Forest Service, other state, local, and private entities as derived from the PADUS [Protected Areas Database of the United States] layer [29]); however, this covariate does exclude smaller or private protected lands not reported to the PADUS group. Furthermore, protected lands are not necessarily forested, with private protected lands comprised of more agricultural and utilitarian sites than public protected lands [30]. Finally, we used Fragstats to measure applied landscape covariates, such as (12) contagion, (13) total edge, and (14) landscape division indices within the same 500 and 1000 m scales as the land cover percentage measurements [31]. We standardized all continuous covariates to z scores for analyses and excluded any covariates that were highly correlated (r > 0.7).
Because camera traps were deployed across years and seasons, we developed and compared 7 a priori models, including a constant (null) and global (year + season) set, to estimate and account for the effects of season-specific covariates on detection probabilities of each species. We then used the covariates that contained high model support and strong effects on species detection as a constant covariate in the subsequent occurrence models for each species.
For the occurrence models, we developed and compared 30 a priori hypotheses to predict the distributions of mammals in our study area dependent on land cover and our derived covariates (Supplemental Table S1). The approach was hierarchical in that we compared 14 preliminary models with habitat covariates at local and broader scales and then used model selection to determine which scale (e.g., 500 m or 1000 m) was more appropriate to include in further additive and global models for each specific species. We developed and compared occupancy models in the "unmarked" package in R [32,33].

Model Selection
We evaluated the best approximating models based on their Akaike information criterion (AIC) and Akaike weights (ω i ). We considered covariate effects to be substantial if 95% confidence intervals excluded 0. To make habitat distribution models for the contemporary time-period, we evaluated model predictive power using the area under the ROC (Receiver Operating Characteristics) curve (AUC) and considered models to have strong predictive power if the AUC was greater than 0.7. We further used the false positive rate (fpr) as an index of our predictive power and applied conservative (fpr < 0.1) and more liberal thresholds (fpr < 0.3) to generate predicted core (conservative) habitat and suitable (less conservative and inclusive of all core) habitats for each species, respectively. We implemented these analyses within the "ROCR" package [34].

Predicting under Land-Planning Scenarios
Using planning scenarios developed for the region via a stakeholder-driven scenario planning process [19], we projected differences in land cover use across the CLI study area. During two regional meetings, local stakeholders defined four scenarios by crossing the poles of two drivers of change: human population growth and the degree of strategic planning. The results are as follows: (HS) high human population growth and strategic planning, (HR) high population growth and reactive planning, (LR) low population growth and reactive planning, and (LS) low population growth and Remote Sens. 2019, 11, 2482 5 of 16 strategic planning ( Figure 2). The stakeholders also provided estimates of the amount of difference they expected to see between scenarios that we used in the quantification process described below.
Remote Sens. 2019, 11, x FOR PEER REVIEW 5 of 17 difference they expected to see between scenarios that we used in the quantification process described below. Figure 2. Schematic of the four quadrants defining our scenario plans along population growth and strategic planning gradients for the CLI study area in northern Virginia.
We simulated these scenarios with the spatially-explicit modeling platform Dinamica EGO (Environment for Geoprocessing Objects) version 3.41 [35]. Dinamica uses a Markov chain/cellular automata framework to project change at one location on a landscape, based on information about the surrounding area. We used that framework to project landscape configurations 50 years into the future, using primary inputs derived from two patterns of land use change observed between our calibration landscapes; NLCD 2001 and NLCD 2011. We parameterized each model by making modifications to the patterns observed.
The first input pattern was the amount of area that changed for each class-to-class transition for each of the four classes (i.e., both forest that becomes grass and grass that becomes forest). This value controls the number of changes Dinamica makes to the initial landscape. We did not track transitions from development to other classes because this was a rare event. We increased the rate of change for development in the HR and HS scenarios to represent the increased need for housing and services created by higher human population density [36]. The amount of change was quantified by first matching observed rates of development growth to observed population growth and then projecting that relationship to decadal population values 24% higher than those projected by the Weldon Cooper Center [37]. For the LR and LS scenarios, we kept the observed rate of development transitions constant to represent low population growth. We then applied a stakeholder-driven, percent-based reduction to both the modified HS rates and unedited LS rates to represent the smaller footprint caused by centralized high density housing expected under strategic planning. We also made modifications for transitions between forest, grass, and crops based on stakeholder input regarding agricultural change [19].
The second input was the "weights of evidence" values calculated for each class-to-class transition in relation to a suite of environmental and socioeconomic variables. These values control where Dinamica makes changes to the initial landscape through maps of the accumulated weights of all variables for each transition type. To distinguish between the centralized growth plans expected in strategic scenarios and the more sprawled growth expected in reactive scenarios, we modified the transition to development weight maps such that the accumulated weights near urban centers were higher for HS and LS scenarios and lower for HR and LR scenarios. We simulated these scenarios with the spatially-explicit modeling platform Dinamica EGO (Environment for Geoprocessing Objects) version 3.41 [35]. Dinamica uses a Markov chain/cellular automata framework to project change at one location on a landscape, based on information about the surrounding area. We used that framework to project landscape configurations 50 years into the future, using primary inputs derived from two patterns of land use change observed between our calibration landscapes; NLCD 2001 and NLCD 2011. We parameterized each model by making modifications to the patterns observed.
The first input pattern was the amount of area that changed for each class-to-class transition for each of the four classes (i.e., both forest that becomes grass and grass that becomes forest). This value controls the number of changes Dinamica makes to the initial landscape. We did not track transitions from development to other classes because this was a rare event. We increased the rate of change for development in the HR and HS scenarios to represent the increased need for housing and services created by higher human population density [36]. The amount of change was quantified by first matching observed rates of development growth to observed population growth and then projecting that relationship to decadal population values 24% higher than those projected by the Weldon Cooper Center [37]. For the LR and LS scenarios, we kept the observed rate of development transitions constant to represent low population growth. We then applied a stakeholder-driven, percent-based reduction to both the modified HS rates and unedited LS rates to represent the smaller footprint caused by centralized high density housing expected under strategic planning. We also made modifications for transitions between forest, grass, and crops based on stakeholder input regarding agricultural change [19].
The second input was the "weights of evidence" values calculated for each class-to-class transition in relation to a suite of environmental and socioeconomic variables. These values control where Dinamica makes changes to the initial landscape through maps of the accumulated weights of all variables for each transition type. To distinguish between the centralized growth plans expected in strategic scenarios and the more sprawled growth expected in reactive scenarios, we modified the transition to development weight maps such that the accumulated weights near urban centers were higher for HS and LS scenarios and lower for HR and LR scenarios.

Species' Habitat Projections
We used the future landscapes for each scenario to generate scenario-specific versions of each of the model covariates. We then ran the top supported occupancy model on each suite of scenario covariates to generate predicted distributions (suitable and core habitat) for each species. We calculated the predicted change in area for both core and suitable habitats between the present and 2060 for each of the planning scenarios. We further divided those habitats into protected and non-protected areas (e.g., protected from future development), to more broadly evaluate the importance of the protected areas network in the study area.

Camera Trap Summary Statistics
Between 2012 and 2018, we obtained 24,866 independent detections of 15 mammals from 36,443 camera trap days (Supplemental Table S2). We evaluated occupancy models for all 15 species, but excluded 10 species from further analyses beyond the occupancy model comparisons. We excluded three species (woodchucks (Marmota monax), striped skunks (Mephitis mephitis), and gray foxes (Urocyon cinereoargenteus)) due to sparse data and five species (eastern gray squirrels (Sciurus carolinensis), eastern fox squirrels (Sciurus niger), Virginia opossums (Didelphis virginiana), northern raccoons (Procyon lotor), and coyotes (Canis latrans)) due to low predictive power of the model (AUC < 0.70). We excluded two additional species (eastern chipmunks (Tamias semistriatus) and eastern cottontails (Sylvilagus floridanus)) due to incongruence in distribution predictions from likely, subtle site-specific detection biases associated with their small size and fine-grained habitat preferences. This process resulted in five representative species for all further analyses: American black bears, bobcats, red foxes, domestic cats, and white-tailed deer (Supplemental Table S3).

Occupancy Model Results and Habitat Associations
The global models were most supported for black bear and bobcat occurrence probabilities, each receiving 100% of the Akaike weight. The AUC values for each of these predictive models were also strong with AUC = 0.81 and 0.80 for bears and bobcats, respectively. As predicted, bears and bobcats were positively associated with natural land covers (forest and grasslands) and were positively associated with increasing distances from anthropogenic land covers of development, roads, and crops. However, contrary to our a priori predictions, the binary protected area covariate was not a strong predictor, with a negative, but non-significant, association (Table 1). The anthropogenic models were most supported for domestic cats and red foxes with 70% and 88% of the Akaike weights, respectively. The AUC values for each of these predictive models were also reliable, with AUC = 0.76 and 0.74 for cats and red foxes, respectively. Habitat associations of these two mesopredators followed our a priori predictions, with both species being positively associated with increasing human development and negatively associated with increasing distances from roads and development and increasing core forest ( Table 2). White-tailed deer's occurrence was most supported by the global model (99% of the Akaike weight) with AUC = 0.73, suggesting the species was negatively associated with increasing forest cover, distances from human development, roads, and crops, and positively associated with grassland habitats and protected area status ( Table 2). Table 2. Estimated beta coefficients (β) with standard errors (SE), 95% lower and upper confidence intervals (LCI, UCI), p-values, and a priori predictions for covariate effects from the top-ranking occupancy model, explaining variation in occupancy and detection effects of land cover and anthropogenic covariates on the probability of site use by synanthropic mammals from 1591 camera trap locations across an urban to rural gradient, northern Virginia, 2012-2018. Habitat associations from our occupancy models of the other species detected generally followed our a priori predictions, with species such as eastern chipmunks, eastern fox squirrels, gray foxes, striped skunks, and coyotes responding positively to increasing forest cover, while eastern cottontails and woodchucks responded positively to increasing grassland cover (Supplemental Table S4). The top-ranking occupancy model to predict raccoon occurrence included the derived landscape division index as the only covariate. Additionally, most species (i.e., coyotes, striped skunks, Virginia opossums, eastern gray squirrels, and eastern cottontails) responded positively to increasing distances from development, whereas only eastern fox squirrel models suggested a negative association. Coyotes, opossums, and gray squirrels were negatively associated with increasing distances from crops, while striped skunks, fox squirrels, and cottontails responded positively to increasing distance from crops (Supplemental Table S4).

Predicted Contemporary Distributions
Contemporary black bear habitat estimates ranged from 4067.7 km 2 for core habitat to 10,547.5 km 2 for suitable habitat in the study area, with 56.9% and 38.9% of those areas in protected areas, respectively (Figure 3a). The contemporary bobcat habitat was more restricted, with estimates ranging from 2268.9 km 2 for core habitat and 5613.3 km 2 for suitable habitat, with the majority of habitat in protected areas, to 75.5% core and 58.7% suitable habitat (Figure 3b).
Contemporary domestic cat and red fox distributions were largely limited in the study area, with > 90% of habitats occurring in unprotected private lands. The distribution for cats ranged from 5660.1 km 2 (core) to 14,583.4 km 2 (suitable), and the distribution of red foxes ranged from 1927.8 km 2 (core) to 11,999.6 km 2 (suitable) (Figure 3c,d). White-tailed deer were the most ubiquitous species in the study area with a 6867.9 km 2 core and an 23,213.5 km 2 suitable contemporary habitat (Figure 3e). distances from crops, while striped skunks, fox squirrels, and cottontails responded positively to increasing distance from crops (Supplemental Table S4).

Predicted Contemporary Distributions
Contemporary black bear habitat estimates ranged from 4067.7 km 2 for core habitat to 10,547.5 km 2 for suitable habitat in the study area, with 56.9% and 38.9% of those areas in protected areas, respectively (Figure 3a). The contemporary bobcat habitat was more restricted, with estimates ranging from 2268.9 km 2 for core habitat and 5613.3 km 2 for suitable habitat, with the majority of habitat in protected areas, to 75.5% core and 58.7% suitable habitat (Figure 3b).
Contemporary domestic cat and red fox distributions were largely limited in the study area, with > 90% of habitats occurring in unprotected private lands. The distribution for cats ranged from 5660.1 km 2 (core) to 14,583.4 km 2 (suitable), and the distribution of red foxes ranged from 1927.8 km 2 (core) to 11,999.6 km 2 (suitable) (Figures 3c,d). White-tailed deer were the most ubiquitous species in the study area with a 6867.9 km 2 core and an 23,213.5 km 2 suitable contemporary habitat (Figure 3e).

Projected Land Cover and Suitable Habitat Availability Changes under Different Scenarios
Overall, changes under the four planning scenarios are driven by human population growth, with high human population density models (HS and HR) exhibiting the most urban development expansion (+381.8-445.2 km 2 ) and highest loss of forest (−277.4-325.6 km 2 ), grasslands (−32.5-36.0 km 2 ), and cropland (−71.8−83.6 km 2 ) ( Table 3). The LS scenario with low population expansion and strategic planning represents the least expansion of urban development (+251.0 km 2 ) and the lowest reduction in forest cover loss (−227.9 km 2 ), compared to LR with low population and reactive planning competing (288.0 km 2 of development gain and 271.4 km 2 of forest loss).

Projected Land Cover and Suitable Habitat Availability Changes under Different Scenarios
Overall, changes under the four planning scenarios are driven by human population growth, with high human population density models (HS and HR) exhibiting the most urban development expansion (+381.8-445.2 km 2 ) and highest loss of forest (−277.4-325.6 km 2 ), grasslands (−32.5-36.0 km 2 ), and cropland (−71.8−83.6 km 2 ) ( Table 3). The LS scenario with low population expansion and strategic planning represents the least expansion of urban development (+251.0 km 2 ) and the lowest reduction in forest cover loss (−227.9 km 2 ), compared to LR with low population and reactive planning competing (288.0 km 2 of development gain and 271.4 km 2 of forest loss).
Predicted core habitat losses for bears and bobcats ranged from low under the LS planning scenario with 2.2% and 3.6% losses, respectively, to relatively high under the HR (4.7% loss) and HS (8.6% loss) scenarios, respectively. However, when the study area is stratified by protected status, it becomes apparent that habitat loss occurs mostly in unprotected habitat, with as high as 9.2% and 23.7% losses for bears and bobcats in scenarios HR and HS, respectively.
Predicted core and suitable habitats increased for domestic cats and red foxes under all planning scenarios. Substantial expansion, such as 39.3% for domestic cats' core habitat and 88.8% core habitat growth for red foxes, was predicted under the HR planning scenario, whereas under the LS planning scenario, habitat expansion was lowest, ranging from 3.0% to 17.4% and 3.2% to 50.6% for suitable versus core habitat expansion for domestic cats and red foxes, respectively. Habitat for white-tailed deer does not vary strongly under any planning scenarios, ranging from a loss of 1.6% under HS to 2.2% gain under the LR planning scenario. Table 3. Estimated core and suitable habitat area and percent changes over time for five representative mammals in protected and non-protected areas during the contemporary period (2011) and under four land planning scenarios into 2060. HS represents high human population growth and strategic planning, HR represents high human population growth and reactive planning, LR represents low population growth and reactive planning, and LS represents low population growth and strategic planning. Habitat areas were predicted from the top-ranking occupancy models with effects of land cover and anthropogenic covariates from 1591 camera trap locations across an urban to rural gradient, northern Virginia, 2012-2018.

Discussion
We were able to model the predicted response of five representative mammal species to four planning scenarios for the study area. Although there was some variation in species-specific responses among the scenarios, we still presented unified patterns for use by land managers and development planners. Future landscapes would be beneficial and useful to both land managers and animal communities if they retained substantial habitat for sensitive species, and minimized expansion by the synanthropic species. Based on those criteria, it is not surprising that our results suggest that LS would be the most useful planning policy, with strategic planning for development and agriculture (and low population growth) retaining the most forest cover and the associated mammal communities.
Although bears and bobcats are sensitive to habitat loss and fragmentation [12,38], our results suggest that future habitat loss in this region is relatively limited, regardless of the planning policies adopted. This result highlights the observation that bears and bobcats are already largely restricted to the two largest forested areas (Shenandoah NP and George Washington NF) and that these areas will not lose significant habitat-appropriate areas in future projections due to their protected status. Counter to our predictions, the β coefficients from the occupancy models suggested a negative but weak association between these species' distribution and protected areas. However, the relationship reveals that, aside from these two large forested areas, most other protected areas remain small or isolated, particularly among private protected lands [30], and their protected status alone is not enough to support these carnivores. In the case of bobcats, habitat losses are projected to be five times greater in unprotected areas under the high human growth scenarios. Bears exhibit less pronounced habitat losses in unprotected areas, likely due to their adaptability to living in exurban areas [39]. These results highlight the importance of protected areas for biodiversity preservation in developing landscapes [40], and their importance will only increase in the near future, as many small private forest patches are lost to development under each of the scenario models.
Both free-roaming domestic cats and red foxes are predicted to expand their distributions in the study area with the expansion of development and anthropogenic land uses under all scenarios. Those results correspond with other studies that show these mesopredators are highly associated with humans [7,8,10]. The impacts of range expansion by both species could be detrimental to biodiversity (e.g., small mammals, songbirds, and herpetofauna), particularly if high densities of cats and foxes are sustained due to subsidies from humans [41][42][43]. However, our predictive mesopredator models do not include interspecific interactions, which might further influence these predictions [43]. Specifically, coyotes were fairly common in our study area and they are known to influence the site use, abundance and behavior of both domestic cats and foxes [10,44,45]. Yet, our occupancy models did not have strong predictive power to accurately assess coyote distributions, and hence we were unable to directly incorporate these trophic interactions that might continue to limit the future expansion of mesopredators.
While our species are representative of relatively sensitive or synanthropic species, it is important to note that these five species and the additional 10 species detected are only a subset of historical mammal diversity from the study area. Apex predators (i.e., wolves (Canis lupus) and mountain lions (Puma concolor)) were extirpated, and other obligate carnivores, such as mustelids, are rare and were likely extirpated from much of the study area. Mesopredators have ascended to be the de facto top predators in many anthropogenic landscapes, yet we are only now learning the roles of many of these medium-sized predators in the context of trait-mediated cascades and their effectiveness in regulating herbivores and lower trophic taxa [28,46]. As such, our results are applicable for a typical eastern US landscape but might not represent the community dynamics of areas with intact apex predator guilds, such as the upper Midwest or Western US.
Our study was one of the first to pair contemporary camera trap data with occupancy models to predict current species distributions and then explore how those distributions might change under alternative planning scenarios in a rapidly developing area. However, we were unable to accurately predict the distributions of several mammalian species in contemporary times, and hence unable to predict potential changes in the future. Many of those species are common and highly associated with humans; e.g., northern raccoons and eastern gray squirrels. These synanthropes are likely responding to more fine-scale covariates, as exemplified by raccoons, most positively associated with housing densities and the associated human refuse, outdoor pet food, bird feeders, and other resource subsidies in anthropogenic and heterogeneous landscapes [47][48][49]. These unmodeled species will likely interact with other taxa, and therefore, warrant further examination with higher resolution habitat and anthropogenic covariates. Furthermore, our land cover predictions and planning scenarios operated over a short time span and did not incorporate potential climate change. Indeed, contemporary studies have shown that small mammals and plant communities are stratified and shift along elevational gradients, with prey responding variably to these synchronous plant community shifts in response to climate change [50,51]. However, since we cannot account for these changes in medium and large mammals or our land cover predictions, we do not expect these changes to affect medium and large mammals strongly beyond the scope of our predicted models, although we might consider our models conservative with respect to climate uncertainty.

Conclusions
Our results, based on contemporary mammal distributions and land cover associations, provide useful generalized patterns of planning scenarios. Those highlight the value of large and contiguous public lands and their continued protection from development and expansion as a buffer to rapidly changing environments. Strategic planning will also benefit mammal communities and multiple stakeholders by retaining habitats for sensitive species and minimizing the expansion of synanthropic species. Moreover, we find that planning and policy decisions will be particularly influential in private landscapes, where there is ample opportunity for increasing protected lands area and targeting landscape features that benefit mammal communities.