Simulating wildlife habitat dynamics over the next century to help inform best management strategies for biodiversity in the Lake Tahoe Basin, California

. Many forests of the western United States have undergone over one hundred years of anthropogenic impacts that have led to increased tree density, homogenization in forest structure, and accumulation of woody material, which combined with a changing climate pose threats to valued social and ecological features. In California, recent waves of tree mortality and unprecedented large and destructive fires have led to rising concerns about the impact of these disturbances on biodiversity and how forest management actions can mitigate negative impacts. To better understand the degree to which different management scenarios could mitigate the negative impacts of these disturbances on biodiversity, we used a spatially explicit modeling platform to model forest management impacts on habitat for terrestrial vertebrate species in the Lake Tahoe Basin of California and Nevada. Specifically, we modeled how 5 different management scenarios that differed in the type of fuel reduction treatment (e.g., fire and mechanical removal of vegetation) and extent of area treated influenced the amount, value, and distribution of reproductive habitat for the 159 species present in the study area. Our model results suggested that within the study area forest growth was predicted to out-pace disturbance leading to a higher percentage of late seral conditions; however, choice of management strategy impacted the composition and structure of the forested landscape leading to different trajectories for wildlife. In general, scenarios that allowed for more extensive use of fire led to a more equitable distribution of habitat types, whereas extensive thinning by hand and mechanical methods resulted in future forest structure that provided better outcomes in terms of reproductive habitat for wildlife. Our modeling results also suggested that low to moderate management strategies were not likely to change the current trajectory to more dense forests dominated by fewer species.


INTRODUCTION
In the seasonally dry forests of the western United States, fire is a natural disturbance that shapes habitat composition, structure, and distribution (Larson and Churchill 2012, Knapp et al. 2013, Lutz et al. 2018). In the mid-elevation forests of the Sierra Nevada Mountains in California, fire historically maintained a speciesrich mixture of tree densities and size classes (Barbour et al. 2002, Beaty and Taylor 2007, North et al. 2009), but due to past practices such as logging, grazing, and highly effective fire suppression, mixed-conifer forests have become denser and simplified in structure and composition (Knapp et al. 2013, Stephens et al. 2013. Coupled with rapid changes in climate, these changes in forest structure have contributed to large-scale susceptibility of trees to bark beetle attacks during drought (Young et al. 2017) and massive, high-severity fires (Lydersen et al. 2014, Young et al. 2017. In California, theses large-scale disturbances threaten the State's economy because these forests not only provide woody material, clean water, and abundant recreational opportunities, but also serve as habitat for a diverse community of plants and animals (Mayer and Laudenslayer 1988). The rapid change in these forests underscores the need to increase the pace and scale of forest treatments designed to reduce woody material and increase forest resiliency to future, large-scale disturbances , Stephens et al. 2013.
Forest treatments designed to reduce the risk of high-severity fire commonly focus on reducing the amount and connectivity of woody material on the forest floor (e.g., branches, logs, and shrubs), between the forest floor and the canopy (e.g., shrubs and small trees), and within the canopy (e.g., smaller, subdominant canopy trees; Agee and Skinner 2005). To implement large-scale forest restoration projects will require not only mechanical thinning of trees and shrubs, but due to the relatively small percentage of forest that can be mechanically treated because of logistical and fiscal constraints (North et al. 2015), will likely require the use of prescribed burning and managed wildfire . Although data collected from forested areas receiving both mechanical/hand removal of vegetation and prescribed burning treatments have been shown to have neutral to positive effects on many wildlife species (Fontaine andKennedy 2012, Stephens et al. 2014), many species thrive in forests burning with moderate to high-severity fire effects (Smucker et al. 2005, Tingley et al. 2016, White et al. 2016, higher cover of shrubs (White et al. 2013), or complex, early seral habitat conditions (Swanson et al. 2011). Collectively, these studies confirm that wildlife response to fuel reduction treatments will depend on the type of disturbance or restoration action and on the resulting structure.
Forest structure, and its composition, provide the environmental conditions needed to support essential food, resting, and breeding requirements for wildlife. Changes to forest structure inherently impact the degree to which environmental conditions provide suitable habitat, i.e., improving habitat conditions for some species, while reducing it for others (Fontaine and Kennedy 2012, White et al. 2013, Stephens et al. 2014. The California Wildlife Habitat Relationships (CWHR) database (https://www.wildlife. ca.gov/Data/CWHR) classifies the suitability of habitat for vertebrate species based on predominate vegetation (tree, shrub, herbaceous, etc.), species composition, and stage of vegetation development (Fig. 1). For each combination of these conditions, CWHR provides a habitat suitability classification (high, We build on previous efforts to better understand how the type and intensity of management treatments are likely to influence wildlife habitat (e.g., Fontaine and Kennedy 2012, White et al. 2013, Stephens et al. 2014) over large spatial and temporal scales. Although providing for habitat connectivity and a mosaic of forest conditions provides general principles for landscape scale wildlife management, specific targets are rarely quantified and will differ with species and existing forest conditions (North and Manley 2012). Using a modeling framework, we assess how different types of forest restoration treatments (i.e., fire suppression, mechanical and hand thinning, prescribed burning, and managed wildfire) and placement, i.e., wildland-urban interface (WUI) or general forest, impact the quality and distribution of reproductive habitat for terrestrial vertebrate species in the Lake Tahoe Basin of California over the next century. We predicted that management scenarios incorporating managed wildfire would lead to greater structure complexity in forest types and subsequently provide greater benefits to wildlife than management scenarios focusing on mechanical treatments in the wildland-urban interface. This effort is part of a larger initiative led by the Lake Tahoe West Restoration Partnership to restore the social-ecological resilience of forests, watersheds, and communities on 59,000 acres of Lake Tahoe's west shore (Abelson et al 2022).

METHODS
The Lake Tahoe Basin, located on the eastern crest of the Sierra Nevada Mountains straddling the states of California and Nevada, supports high levels of biodiversity (Fig. 1). The location of the basin between two major biogeographic provinces, the Sierra Nevada and the Great Basin, creates a juxtaposition of habitats, each province of which supports a unique array of species. In addition, the basin's steep elevational gradients (ranging from 1900 m at Lake Tahoe's surface to 3400 m at the highest mountain peak) span lower montane, upper montane, subalpine, and alpine habitats. Within these habitats, differences in forest structure and composition further influence the provisioning of resources used by wildlife for reproduction, foraging, and cover. Our assessment focused on the habitat of 159 vertebrate species that have suitable habitat for breeding in the Lake Tahoe Basin.

Dynamic modeling of habitat
To model habitat change over time, we used the landscape disturbance and succession model LANDIS-II (https://www. landis-ii.org/) to simulate forest change as a function of growth and succession (Mladenoff 2004, Scheller et al. 2007). Within LANDIS-II, trees and shrubs are modeled as species-age cohorts; each species has its own unique life history attributes and as such, each responds uniquely to disturbance. Every cell (1-ha resolution) across the landscape has a unique combination of cohorts, live biomass, dead material, and soils as determined by its simulated history of succession and disturbance type and intensity. In addition to fire and management, several species of bark beetles, including the Jeffrey pine beetle (Dendroctonus jeffreyi), Mountain pine beetle (D. ponderosae), and fir engraver (Scolytus ventralis) are a frequent disturbance in the study area that was modeled. Due to previous efforts, LANDIS-II has been calibrated to conditions in the Lake Tahoe Basin (Loudermilk et al. 2013, Scheller et al. 2019).
Using the CWHR classification scheme, we developed a rule set for LANDIS-II biomass outputs that would approximate CWHR habitats (Appendix 1). For this study, we modeled change in treedominated and shrub-dominated CWHR habitats (Fig. 1A). Areas that would be classified as aquatic, developed, or barren following CWHR classification rules were excluded from the landscape analysis because these areas were not projected to change over time. In addition, LANDIS-II does not model the herbaceous community, and we classified herbaceous-dominated habitat structure according to the woody vegetation component (i.e., tree or shrub-dominated) or as NULL following recent disturbances that had removed all woody vegetation. The shrub component of different habitat types can vary, and we classified an area as shrub dominated when ≥75% of the total biomass was attributed to shrubs (Appendix 1). This biomass threshold was chosen based on LANDIS-II biomass estimates of the current landscape compared to current classification of CWHR shrubdominated habitats. Species composition in tree-dominated habitats followed the CWHR classification scheme but relied on species-specific biomass estimates instead of estimates of canopy cover. Species-specific biomass estimates produced by LANDIS-II to classify composition included: white fir (Abies concolor), red fir (A. magnifica), lodgepole pine (Pinus contorta), Jeffrey pine (P. jeffreyi), incense-cedar (Calocedrus decurrens), sugar pine (Pinus lambertiana), mountain hemlock (Tsuga mertensiana), western white pine (P. monticola), and whitebark pine (P. albicaulis). Aspen (Populus tremuloides) is the only hardwood species that occurs regularly within the study area and was the only treedominated habitat modeled that was not classified as conifer. Based on this rule set, we identified eight CWHR habitat types in our study area: White Fir (WFR), Red Fir (RFR), Jeffrey Pine (JPN), Lodgepole Pine (LPN), Sierra Mixed Conifer (SMC), Subalpine Conifer (SCN), Aspen (ASP), and Montane Chaparral (MCP).
For conifer dominated habitat types, we assigned each habitat type a seral stage based on CWHR size class descriptions (Fig.  1B). There are five CHWR size classes that influence habitat suitability for wildlife based on the mean diameter at breast height (DBH) of predominant trees. Seedling trees are the first size class and are defined by trees with a mean diameter < 2.5 cm, followed by sapling trees (size class 2: 2.5-15 cm), pole trees (size class 3: 15-28 cm), small trees (size class 4: 28-61 cm), and med/large trees (size class 5: > 61 cm). We estimated these size classes by combining estimates from LANDIS-II for the biomass of different species age cohorts with data from Forest Inventory and Analysis (FIA; Burrill et al. 2021) on the distribution of stem sizes (for equivalent species and age) in trees that have been measured in the northern Sierra Nevada ecoregion where Lake Tahoe is situated. We fit the Weibull distribution to FIA tree data for species-age groups that were equivalent to each LANDIS-II cohort. We assumed that the LANDIS-II cohorts would follow a similar size distribution to what was observed in the FIA data and distributed the LANDIS-II biomass estimates among twoinch DBH height bins (i.e., into a 0-5 cm bin, a 5-10 cm bin, etc.). To assign a pixel to a CHWR size class, we used the DBH of the 95 th percentile tree in our estimated distribution and compared that to the CHWR size classes as described above (classes 1 through 5 being < 2.5 cm; 2.5-15 cm; 15-28 cm; 28-61 cm; and > 61 cm). We assessed different approaches, specifically, we also completed the classification using the size of the 90 th percentile tree and the 99 th percentile tree, and we found that using the 95 th percentile yielded a CHWR in year 1 (2010) of the LANDIS-II model run that most closely reflected our baseline habitat distribution in the Lake Tahoe Basin. We considered this method reliable in producing three distinct tree-dominated seral stages: early (tree size approximating size classes 1-3), mid (tree size approximating size class 4), and late (tree size approximating size class 5).
Canopy cover in early seral habitat did not have a predictable relationship with canopy cover estimates, but for mid and late conifer-dominated habitats, we additionally assigned each pixel a canopy cover estimate (open, moderate, closed) based on CWHR canopy classes (Fig. 1C). Using the biomass bins described above, we again used the FIA data for average stem to biomass relationships at different size classes to translate biomass values into stem counts for each size bin (Evans et al. 2022). This process gave us an estimated size distribution of trees on each pixel as well as the estimated number of trees per pixel, which we translated into an estimate of trees per acre (TPA) to estimate canopy cover. We based our canopy thresholds on estimates of TPA whereby mid and late seral habitats with < 50 TPA were categorized as open, 51-500 TPA as moderate, and > 500 TPA as closed. Using these methods, we assigned each 1 ha pixel a CWHR habitat (composition, size class, canopy class). Although the CWHR classification scheme separates open from sparsely covered habitats, we combined these two into one category creating three canopy classes that are ecologically important and are used to establish guidelines for management. Descriptions of CWHR habitats modeled and the number of hectares of each habitat combination produced are presented in Appendix 2.

Simulated management scenarios
Management scenarios simulated in this study were selected by the Lake Tahoe Basin West Restoration Partnership to better understand the societal benefits of increasing the pace and scale of forest restoration (Table 1). Under the suppression-only (Suppress) scenario no forest management activities were modeled other than active fire suppression (i.e., if a fire started, it was suppressed). Under the wildland-urban interface (WUI) scenario, in addition to fire suppression, forest thinning (both mechanical and hand treatment) was modeled in the vegetated area adjoining structures and evacuation routes (which accounts for ~58% of the landscape). Similar to the WUI scenario, the thin scenario (Thin) modeled high levels of forest thinning in the WUI, as well as some treatments in the general forest and surrounding wilderness areas. The WUI and thin scenarios thus differed in the Understand the societal benefits of fire suppression.
Wildland-urban interface (WUI) In addition to fire suppression, approximately 1.8% of the vegetated area was treated annually in the WUI and defense zone.
Understand the societal benefits of current management actions focused on reducing fuels around human-built structures and evacuation zones. Thin In addition to fire suppression and WUI-focused fuel reduction treatments, this scenario treated approximately 6.7% of the vegetated area annually including the forest beyond the WUI.
Understand the societal benefits of expanding the use of mechanical and hand-thinning into the general forest.

Fire
Management actions were focused on reducing fuels via prescribed fire and managed wildfire with approximately 4% of the vegetated area treated each year.
Understand the societal benefits of using fire-driven fuel reduction treatments in addition to mechanical and handthinning mechanisms. Fire + Management actions were focused on reducing fuels via prescribed fire and managed wildfire with approximately 7.2% of the vegetated area treated each year.
Increase the amount of vegetated area treated each year via fire-driven mechanisms (75% of all fuel reductions) to better differentiate between the societal benefits of fuel reduction mechanisms and the area treated. amount and extent of treatment. The two remaining scenarios were modeled with the use of fire as a management tool; both modeled moderate levels of thinning in the WUI while allowing for some prescribed burns and managed wildfire ignitions outside of the WUI. The fire scenario (Fire) however, simulated moderate levels of prescribed and managed fire, whereas the Fire + scenario simulated high use of prescribed and managed fire. Information on the parameterization of the model is available here: https:// github.com/LANDIS-II-Foundation/Project-Lake- Tahoe-2017 For each of the five scenarios, annual change in forest conditions were simulated from 2010-2110 using climate data from the CanESM2 global change model based on the relative concentration pathway 4.5. Each scenario run was simulated 10 times (i.e., 10 replicates) to quantify the mean and variation in change in habitat. At each 10-year timestep, we assigned CWHR habitat types to each spatial output (5 scenarios X 10 replicates X 10 decadal time steps) resulting in 500 spatially explicit habitat maps at 1-ha resolution.

Simulating wildlife-habitat outcomes
We extracted data from CWHR (https://www.wildlife.ca.gov/ Data/CWHR) on wildlife-habitat relationships for each vertebrate species known to occur in the Lake Tahoe Basin (Schlesinger and Romsos 2000). We excluded any species from the analysis that relied primarily on aquatic habitats (e.g., amphibians and waterfowl). This resulted in a list of 135 birds, 51 mammals, and 5 reptiles. Of these 191 species, we were unable to simulate change in reproductive habitat for 23 birds and 2 mammals because they are not known to breed in the study area, or only breed in urban environments. Additionally, we excluded species in which their reproductive habitat was rare (< 5% of the landscape) including the Loggerhead Shrike (Lanius ludovicianus), Bullock's Oriole (Icterus bullockii), White-crowned Sparrow (Zonotrichia leucophrys), Tree Swallow (Tachycineta bicolor), Lark Sparrow (Chondestes grammacus), and Gray-crowned Rosy Finch (Leucosticte tephrocotis). This resulted in 159 species in which a change in the value of habitat for reproductive purposes was modeled (Appendix 3).
For each modeled species and every habitat map (N = 500), we assigned a habitat-suitability ranking of high (1.00), moderate (0.66), low (0.33), or unsuitable (0.00) for reproduction based on the CWHR database for each 1-ha pixel ( Fig. 2A). Because our translation from LANDIS II outputs (biomass by species by age cohort) to CWHR habitat resulted in the decision to combine certain CWHR habitats (i.e., sparse and open canopy classes were combined), we assigned the maximum value of reproductive habitat (hereafter habitat) for the combined categories (Appendix 3). From each map, we generated three metrics: (1) the total sum of habitat suitability values on the landscape for each species, (2) the proportion of the landscape classified as high-value habitat (i.e., suitability = 1; Fig. 2B), and (3) the average size of habitat patches of high reproductive value using ClassStat function in package SDMTools (McGarigal et al. 2002) in Program R version 3.6.1 (R Core Team 2018). For each species*scenario*replicate* timestep, we compared the computed metric to baseline starting conditions in 2010. For each metric (overall habitat quality, proportion of landscape considered high-quality habitat, average patch size of high-quality habitat) at each scenario*replicate* timestep, we summed the number of species in which each metric was predicted to increase and decrease. Because these metrics always fluctuate due to the dynamic nature of the model, we chose an arbitrary threshold of > 30% change to quantify more significant impacts on species and quantified these species separately. Similar to the predicted outcomes for habitat type composition, management scenarios had little influence on the proportion of the landscape comprised of early, mid, and late seral stages over time (Fig. 3). Of the conifer dominated habitat types, 63% of the cells were classified initially as mid-seral, 26% as late seral, and 11% as early seral (the latter of which primarily occurred in White Fir and Red Fir habitat types). Under all management scenarios, the amount of early and mid-seral forest decreased whereas late seral forest increased. In the Fire + scenario, conversion of more areas to Montane Chapparal resulted in a lower proportions of conifer forest, with late seral habitat conditions plateauing at 22,500 ha and the smallest percentage of early seral conditions.
Compared to habitat type composition and seral stage, forest cover was more strongly influenced by management actions that increased the rate and extent of fuel reduction activities (Fig. 3). Compared to the management scenario simulating fire suppression (Suppress), fuel reduction treatments concentrated in the WUI (WUI) and with only moderate fire use (Fire) had a stabilizing effect on the proportions of the landscape in which forest cover was classified as moderate and closed. Increasing the pace and scale of fuel reduction treatments (Thin) was projected to reduce the amount of closed habitat conditions relative to the Suppress scenario but resulted in the largest proportion of the landscape being comprised of moderate forest cover. With the simulated increase in use of fire in the Fire + scenario, forest cover across the landscape was projected to become sparser, with the most limited area of closed forest conditions. For the remaining species the value of the landscape for reproductive purposes, both in the amount and configuration, differed by scenario. In terms of the overall value of reproductive habitat, all scenarios except the Fire + scenario, led to outcomes that resulted in slightly more species experiencing increases in this metric than decreases (Fig. 4A, B). Under the Fire + scenario, more species were predicted to experience declines in the overall amount of habitat for reproduction, although under this scenario the number of species experiencing significant change (i.e., > 30%) was minimized relative to the other scenarios ( Fig. 4B).

Changes in suitability of reproductive habitat for wildlife
In terms of high-quality reproductive habitat, all scenarios tended to have relatively neutral impacts on the proportion of their habitat categorized as high value. In the first half of the century, the Fire + scenario resulted in increases in this metric for the largest number of species, whereas the Suppress scenario resulted in the least number of species (Fig. 4C). In the second half of the century, differences in this metric between scenarios were minimal and largely neutral. Similar to the results for the impact of scenarios on the overall habitat quantity of individual species, the Fire + scenario tended to reduce the number of species experiencing more significant fluctuation in this metric during the first half of the century (Fig. 4D). For those species with > 30% fluctuation in the proportion of their high-quality reproductive habitat, more species were predicted to see increases in this metric relative to decreases, for all scenarios (Fig. 4D). Although the overall proportion of high-quality habitat was predicted to increase for as many species as it decreased, the average patch size of individual patches of this habitat was predicted to decrease for most of the species modeled, most notably under the Fire + scenario (Fig. 4E, F). Differences between all other scenarios was minor and the variance in replicates overlapped, but the Thin scenario was predicted to lead to the most favorable biodiversity outcomes in the size of high-quality patches (Fig. 4C).
As patches of high-quality reproductive habitat were predicted to decrease over time regardless of scenario, we grouped modeled species based on their size and mobility and compared how each scenario was projected to impact the median size of high-quality reproductive patches for each grouping in 2040 (early century) and in 2070 the mid-century (Table 2, Appendix 3). For species with low mobility that tended to occupy small to moderate sized home ranges (≤ 40 ha; Manley and Schlesinger 2000), the median size of habitat patches decreased over the century regardless of scenario. These decreases were most dramatic under the Fire + scenario both early and mid-century. Relative to the other scenarios, the Thin scenario was predicted to result in larger median patch sizes in the near term. In comparison, for largerbodied species, the median patch size was predicted to be the smallest relative to other species groups and more similarly across time periods, with the Thin and Fire + scenario leading to higher median patch sizes. For highly mobile species that typically occupy home ranges > 40 ha, patch size decreased over time regardless of scenario with more favorable outcomes again under the Thin and Fire + scenario, both scenarios of which treated larger proportions of the landscape.

DISCUSSION
Recent extreme fire years in forested ecosystems across the globe have raised concerns regarding the potential impact of these climate-driven fires on biodiversity (Barlow et al. 2020, Kelly et al. 2020, Ward et al. 2020) and actions that can be taken to increase resilience of forests to future fires (Wintle et al. 2020, Stephens et al. 2021). In the seasonally dry forests of the western United States, forest restoration strategies that reduce tree densities and fuel loading have been shown to alter fire behavior and reduce fire severity Moghaddas 2005, Tubbesing et al. 2019). Our model results suggest that management regimes that incorporated more extensive use of fuel reduction treatments and fire can alter forest trajectories and habitat for biodiversity.
With the exception of a management regime incorporating more prescribed and managed wildfire (Fire + ), our models indicated that forests surrounding Lake Tahoe will experience an increase in White Fir and Sierra Mixed Conifer and a decrease in stands dominated by Jeffrey pine and montane chapparal. This prediction is supported by previous research indicating that Jeffrey pine is adapted to a frequent-fire disturbance regime whereby frequent removal of small diameter trees and fireintolerant tree species by low-intensity fire is needed to maintain this species on the landscape (Fulé et al. 2012). Our modeling results suggest that higher use of fire is needed to maintain Jeffrey Pine habitat in our study area. Previous research has also shown that red fir seedling establishment and growth are strongly related to periodic disturbance (Chappell and Agee 1996) and this may explain why the Fire + scenario was also the only scenario predicted to maintain Red Fir habitat. Climate change is predicted to have a negative impact on red fir (Meng et al. 2015) as establishment of red fir seedlings in influenced by soil moisture and post-fire temperatures and maintaining Red Fir habitat will be a future management challenge. Although maintaining a more even mix of habitat types, increased use of fire under the Fire + scenario led to the highest proportion of montane chapparal. Many nonconifer vegetation communities have been negatively impacted by fire suppression including montane chapparal and riparian communities, including aspen (White and Long 2019), and our modeling results indicate that higher use of fire may be one way to ensure persistence of these habitat types in the Lake Tahoe Basin.
Whereas frequent use of fire appeared to be important in maintaining a diversity of vegetative communities, the area being treated led to larger differences among scenarios in the density of forest canopy. Increasing levels of disturbance resulted in landscapes with a higher percentage of the landscape comprised of moderate tree densities and associated canopy (Thin scenario) or low tree densities and open canopy (Fire + scenario). Over the century, present day areas with open canopy were predicted to transition to areas with moderate canopy cover consistent with forest growth and succession with limited disturbance and competition for resources. This trajectory is similar to that observed over that last century with increasing stand densities of smaller-sized trees and less structural variability (Taylor 2004). Landscape restoration strategies incorporating more extensive fuel reduction resulted in forests with more moderate tree densities that are less prone to stand-replacing fires (Agee and Skinner 2005) and more resilient to a changing climate (North et al. 2019), but that also may be less conducive to species preferring closed canopy conditions.
One of our main objectives was to understand the impact of different landscape management scenarios on wildlife in the Lake Tahoe Basin based on wildlife-habitat associations. Our modeling results suggest that over the next century changes in wildlife Ecology and Society 27(2): 31 https://www.ecologyandsociety.org/vol27/iss2/art31/ habitat will occur in quantifiable ways, but predictions for most species indicate that they will experience habitat fluctuations that are less than 30% of their current habitat regardless of management approach, suggesting that over the long-term the trajectory of forest change will have an overall neutral effect on biodiversity. Although the Fire + scenario was predicted to promote greater representation of different habitat types, it was predicted to result in a reduction in the overall value of the habitat for many species despite maintaining more consistent amounts of high-quality reproductive habitat than the other scenarios (Fig.  4). This suggests that in our study area, structural variables of forest habitat were more influential in maintaining habitat for biodiversity than forest composition.
The reduced frequency and extent of fire have caused an increase in forest density, a shift to more fire-sensitive species and largescale homogenization of forest conditions (Taylor 2001, Safford and. Over the next century, our simulations project that disturbances altering forest structure will break up the landscape regardless of scenario, reducing the current size of large contiguous forest patches. The distribution of habitat across landscape will impact species differently depending on the distribution of needed resources and their ability to transverse landscapes (Kindlmann and Burel 2008). Based on the size of high-valued reproductive habitat patches, we attempted to predict how different management approaches were likely to impact different groupings of wildlife based on their typical home-range size and mobility (Manley and Schlesinger 2000). We acknowledge that the habitat matrix (i.e., inter-patch habitat) can be as important as the size of habitat patches (Tischendorf and Fahrig 2000), but this is beyond the scope of this study. In general, we found that patches of contiguous high-quality reproductive habitat would decrease for the majority of species, most notably under the Fire + scenario (Fig. 4D, E). Management scenarios that were simulated at greater extents (i.e., the Thin and Fire + scenarios) tended to result in larger patch sizes for the more mobile species and those with larger home ranges. This result in combination with the metric of overall habitat quality for wildlife indicates that management strategies are likely to have idiosyncratic effects of wildlife species, but biodiversity in this disturbance-prone system is fairly resilient to change (see also Sollmann et al. 2015, Kelt et al. 2017. However, in our model we assume that habitat suitability is strongly linked to forest composition and structure, and we do not account for other abiotic and biotic factors that determine the capacity of a given habitat to support the needs of a species.

CONCLUSIONS
Change is predicted, and management can have an impact on the trajectory and speed of change for habitat types and associated species. Choice of best management approach is a societal decision (Abelson et al. 2022) and will depend on prioritization of resources and outcomes, of which biodiversity is just one. Our results are consistent with other studies (Fontaine and Kennedy 2012, Stephens et al. 2014, Stevens et al. 2016 suggesting that different management actions will alter habitat suitability for biodiversity, but that the effects will be species-dependent with compensatory effects on species richness. However, our simulations are based on a representative concentration pathway (rcp 4.5) that does not match our current global emissions (Schwalm et al. 2020). It is likely that use of a higher greenhouse gas emissions scenario would increase the amount of fire on the landscape, leading to more frequent disturbance and landscape conditions most similar to the Fire + scenario. We have only considered one metric for measuring biodiversity without considering the functional role of individual species (Caro 2010, Slauson et al. 2022), which may be as important in sustaining ecosystem function as maintaining current species richness. Ensuring the sustainability of forested ecosystems into the future will continue to be a challenge for managers and society, and integrated frameworks are needed to assess the value of future landscapes (Stevens et al. 2016, Abelson et al. 2022

Acknowledgments:
We acknowledge the indigenous Native American people, the Washoe, who originally lived around the Lake Tahoe Basin and pay our respect to them, their culture, and their Elders, past, present, and future. Our research was informed by continual discussions with many dedicated people who actively participated the in the Lake Tahoe West Restoration Partnership, and it is through these conversations that progress can be best made. In particular, we acknowledge the contributions of P. Manley, J. Long, S. Di Vittori to project management, vision, and oversight  The biomass composition of each pixel in the LANDIS-II output is evaluated to approximated predefined habitat types including montane chaparral (MCP), aspen (ASP), white fir (WFR), red fir (RFR), lodgepole pine (LPN), Jeffrey pine (JPN), Sierra mixed conifer (SMC), and subalpine conifer (SCN). After habitat type is determined, conifer-dominated pixels are classified as early, mid or late seral conditions according to average tree size estimates. Finally, canopy cover class is estimated in mid to late seral forests using estimates of tree density. https://www.ecologyandsociety.org/vol27/iss2/art31/ APPENDIX       LPN2M  LPN2D  LPN3S  LPN3M  LPN3D  RFR1  RFR2S  RFR2M  RFR2D  1  1  1  1  1  1  1  1  1  1