Applying Biodiversity Metrics as Surrogates to A Habitat Conservation Plan

: Unabated urbanization has led to environmental degradation and subsequent biodiver ‐ sity loss across the globe. As an outcome of unmitigated land use, multi ‐ jurisdictional agencies have developed land use plans that attempt to protect threatened or endangered species across selected areas by which some trade ‐ offs between harm to species and additional conservation approaches are allowed among the partnering organizations. Typical conservation plans can be created to focus on single or multiple species, and although they may protect a species or groups of species, they may not account for biodiversity or its protection across the given area. We applied an approach that clustered deductive habitat models for terrestrial vertebrates into metrics that serve as surro ‐ gates for biodiversity and relate to ecosystem services. In order to evaluate this process, we collab ‐ orated with the partnering agencies who are creating a Multi ‐ Species Habitat Conservation Plan in southern California and compared it to the entire Mojave Desert Ecoregion. We focused on total terrestrial vertebrate species richness and taxon groupings representing amphibians, birds, mam ‐ mals, and reptiles, and two special status species using the Normalized Index of Biodiversity (NIB). The conservation planning area had a lower NIB and was less species rich than the Mojave Desert Ecoregion, but the Mojave River riparian corridor had a higher NIB and was more species ‐ rich, and while taxon analysis varied across the geographies, this pattern generally held. Additionally, we analyzed desert tortoise ( Gopherus agassizii ) and desert kit fox ( Vulpes macrotis arsipus ) as umbrella species and determined that both species are associated with increased NIB and large numbers of species for the conservation area. Our process provided the ability to incorporate value ‐ added sur ‐ rogate information into a formal land use planning process and used a metric, NIB, which allowed comparison of the various planning areas and geographic units. Although this process has been applied to Apple Valley, CA, and other geographies within the U.S., the approach has practical application for other global biodiversity initiatives.


Introduction
Unabated urbanization has led to environmental degradation and subsequent biodiversity loss across the globe [1]. Conservation biology began as a reaction to these emerging environmental challenges, using ecological principals and applying them to address species, communities, and ecosystems modified by human activities or natural causes [2,3]. The move to a more proactive, rather than reactive, effort has been a goal of the conservation community. There have been several notable conservation successes such as species removed from listing on the United States Endangered Species Act (ESA) and planning to inform habitat conservation versus development [4][5][6].
The ESA was enacted in 1973; the purpose of the ESA is to provide a means whereby the ecosystems upon which endangered and threatened species depend may be conserved. Moreover, it specifically prohibited the "take" of endangered species [7]. Congress added a provision (ESA Section 10(a)1(B)) in 1982 to recognize that take was not always counter to conservation. This provision allowed the agencies responsible for enforcing the ESA, i.e., Department of Interior's (DOI) US Fish and Wildlife Services (USFWS) and National Oceanic and Atmospheric Administration's (NOAA) National Marine Fisheries Service (NMFS), to issue permits for the incidental take of listed species in certain circumstances. Incidental take can occur during a lawful action when a take is attendant to that action [8]. To receive an incidental take permit, applicants must submit and have approved a habitat conservation plan (HCP) based on guidelines provided by USFWS and NMFS [9,10]. As a result, multi-jurisdictional agencies have developed land use plans that attempt to protect threatened or endangered species across a project area by which some trade-offs between harm to species and additional conservation approaches are allowed among the partnering organizations. HCPs have been created to focus on single or multiple species, i.e., Multiple Species Habitat Conservation Plan (MSHCP), and although they may protect a species or groups of species they typically do not account for biodiversity or its protection across the given area [11].
In 2000, the United Nation's (UN) Secretary-General initiated the Millennium Ecosystem Assessment (MEA) [12], engaging 1360 experts worldwide. Their findings provide a state-of-the art scientific appraisal of the condition and trends in the world's ecosystems and the services they provide. Everything within the MEA paradigm is viewed in terms of human experience and values. Humans are regarded as the central element or factor. These ecosystem services represent the fundamental outputs of a biophysical function or process that contribute, directly or indirectly, to human well-being including economy, health, and a sense of security [13]. The central premise of the MEA ecosystem services framework is that all forms of "Life on Earth", i.e., biodiversity, provide the core benefits that humans derive from their environment.
Biodiversity has been defined as "the variability among living organisms from all sources including terrestrial, marine, and other aquatic ecosystems and the ecological complexes of which they are part" [14]. The term includes diversity among multiple levels of organization from genetics to species, between species, and within ecosystems [15]. Because of the complexity of biodiversity, measuring it as a single metric is difficult. The measurement of surrogates such as indicator species and groups of species has been used assuming there are relationships with unmeasured species [16]. Other direct measures focus on subsets of species, but without those assumptions. One measure often used is species richness (the number of different species in a given area) such as all terrestrial vertebrates, amphibians, birds, mammals, and reptiles, which is informative and can provide a baseline to analyze the value these species provide to humans [13,17,18]. Additionally, in some conservation practices, management actions which rely on one umbrella species are assumed to benefit other species that co-occur within the same habitat and range [19,20]. Thus, the umbrella species concept has been used as a surrogate for multi-species management focused on size and type of habitat to conserve [20], whereby the minimum area required by the umbrella species is nearly as comprehensive as the remaining planning area [21].
The citizens of the United States place considerable economic value on wildlife, not only based on actual spending for recreational activities, but also through investments by state, tribal, and federal governments and non-governmental organizations for conservation planning and easements. As early as 2013, Southwick Associates estimated the conservation economy contributed $93 USD billion to economic output, $41 billion to income, and $59.7 billion to the U.S. Gross Domestic Product [22]. Further analysis in 2017 by Hohensee et al. [23] calculated the total estimated annual economic benefit of terrestrial wildlife in the U.S. to be $346 billion. Based on the U.S. Bureau of Economic Analysis' (BEA) macroeconomic and industry statistics, wildlife economics related to food, materials, recreation, and cultural and aesthetic goods and activities are generally comparable to the activity generated by traditional economic sectors such as agriculture, oil and gas, utilities, and computer and electronics. This BEA effort did not include aesthetic values such as wildlife art, festivals, museum visitation, photography, film, magazines, and television [23].
The U.S. Environmental Protection Agency (EPA) along with other federal, state, and non-governmental organizations created the online EnviroAtlas tool by identifying, creating and serving datasets to help characterize ecosystem services [24]. EnviroAtlas gives users the ability to view, analyze, and download geospatial data and other resources. The basic design and purpose of the EnviroAtlas is to inform decision-making, education, and research. Data within the EnviroAtlas are organized into seven ecosystem service benefit categories: clean air; clean and plentiful water; biodiversity conservation; food, fuel, and materials; natural hazard mitigation; climate stabilization; and recreation, culture, and aesthetics. Data include geospatial indicators and indices of ecosystem goods and services. The EnviroAtlas has analytic and interpretive tools to allow users to conduct analyses for one point in time based on based on all or parts of the geospatial data. The Envi-roAtlas also highlights the connections among the natural resources that provide ecosystem goods and services, the benefits or beneficiaries, and the stressors that affect their provision by measuring and mapping ecosystem services for the conterminous U.S. at the 12-digit Hydrologic Unit Code (HUC) [25] resolution (~104 km 2 in size on average).
The U.S. Geological Survey's (USGS) Gap Analysis Project (GAP) provides data for regional and national analyses of terrestrial vertebrate species and plant species conservation status [18]. GAP is a long-term project created under the U.S. DOI designed to identify gaps in representation of biological diversity using a national land cover map, a national protected areas database, and species-specific deductive habitat models [18,26,27]. The project created species-level habitat models for amphibians, birds, mammals, and reptiles across the conterminous U.S. These models identify habitat suitability within the known range of the species and can be combined to identify potential species habitat richness including within national protected areas [28]. Boykin et al. [13] used GAP species habitat data and the EnviroAtlas approach to identify and test a suite of metrics that serve as proxies for biodiversity and represent several wildlife related ecosystem service endpoints in the Southwest United States. The effort was expanded to the national level based on the 1590 individual GAP habitat species models and 26 national biodiversity metrics which were developed for the conterminous U.S. to reflect aspects of biodiversity conservation; food, fiber, and materials; and recreation, culture, and aesthetics, or their combination (Table 1). These data provide the ability to compare sites at multiple scales anywhere within the conterminous U.S. and to be used in decision-making processes.
The Town of Apple Valley and San Bernardino County, California, are preparing a MSHCP and Natural Community Conservation Plan (NCCP) to address impacts of local development to natural communities and federal and state listed species for the next 30 years. The NCCP is the state component of the overall plan and addresses the state requirements for take under the California Endangered Species Act (CESA) [29]. The process covering this planning is based on USFWS guidance [9], however the Apple Valley MSHCP has yet to be finalized [30]. The MSHCP intent is to protect natural resources while allowing economic growth and orderly build-out of the town and county. The MSHCP seeks to establish a landscape-scale conservation approach versus requiring project-by-project mitigation. Implementation of the MSCHP can minimize and mitigate for impacts to sensitive species and natural communities and provide one-stop assessment for project approval and receipt of environmental permits at the local level. MSHCPs are developed in collaboration with agencies such as the USFWS, California Department of Fish and Wildlife (CDFW), and U.S. Bureau of Land Management (BLM). The Apple Valley MSHCP includes a "Linkage Design" that will help connect the management of 6.5 million acres of state, private, and federal lands designated for conservation. The Apple Valley MSHCP proposes to conserve and manage the lands within this aggregate. The MSHCP purpose is to protect quality of life, including biodiversity, and environmental and scenic values for future generations. The delineated MHSCP area has four finer scale planning area zones (PAZs) including the Urban Core, Rural Development Area, the Desert Linkages, and the Mojave River Corridor.
Boykin et al. [13] developed an approach to produce quantitative biodiversity habitat metrics and maps from the USGS GAP deductive habitat models in collaboration with the EPA EnviroAtlas [13,18]. The purpose of the approach was to allow analysis for placebased (e.g., MSHCP), regional, and national assessments to measure and track changes in the environment and biodiversity. These biodiversity metrics provide indicators, such as species richness, that relate to ecosystem service categories of biodiversity, recreation, and food resource endpoints which are also listed as specific benefit categories in the EPA EnviroAtlas.
This case study was designed to test the efficacy and efficiency of incorporating biodiversity metrics with a real-world application of conservation planning. The objectives of this research were to determine the validity of using value-added biodiversity metrics to provide enhanced analysis for a known geography of key biodiversity interest, i.e., the Apple Valley MSHCP. The analysis was performed in two stages. First, by analyzing the Apple Valley MSHCP area and the Apple Valley PAZs, i.e., Urban Core and Rural Development Areas, Desert Linkages, and Mojave River Corridor, in comparison to the entire Mojave Desert Ecoregion (MDE) in terms of biodiversity metrics. The second objective was to compare biodiversity metrics to two potential umbrella species covered under the MSCHP (desert tortoise and desert kit fox) and evaluate the ability of these two species to provide full protection to the land-use planning area. Table 1. National biodiversity metrics that reflect (A) biodiversity conservation; (B) food, fiber, and materials; and (C) recreation, culture, and aesthetics [31]. Also see additional descriptive information about each metric in the associated EnviroAtlas Fact Sheet ((https://www.epa.gov/enviroatlas/enviroatlas-fact-sheets) accessed on 25 March 2021).

Partners in Amphibian and Reptile Conservation (PARC) Species Richness (A)
Species listed as either T&E, G1-G3 or IUCN (n = 167) as measured by modeled habitat present within a pixel in the conterminous U.S.

Study Area
The Apple Valley MSHCP/NCCP (hereafter "Plan") is a multi-jurisdictional area located within the high desert of San Bernardino County, the largest county in the United States. Stakeholders in the study area include federal (35,059 hectares), state (2973 hectares), county (32,068 hectares), private and town (19,084 hectares), and utilities and other special districts (1990 hectares). Additionally, the Plan provides landscape-level linkages that connect over 2.6 million hectares of existing conservation areas within the MDE (Figure 1).

Landscape-Level Linkages
The Town of Apple Valley has proposed a landscape-level linkage conservation strategy [30]. The purpose is to connect large open areas of natural habitat to adjacent open spaces or habitat areas. The effort delineated the overall Plan area comprised of the following PAZs: the urban core and rural development area (Development), a desert linkages zone (Linkages), and a Mojave River corridor zone (Corridor). The desert linkages identified connect the Plan Development Zone to the San Bernardino National Forest to the south and to three desert tortoise critical habitat units to the east and north. The Corridor is an important riparian area within the surrounding arid environment. The river area includes designated critical habitat for southwest willow flycatcher (Empidonax traillii extimus), one of the nine proposed covered species. The proposed covered species are not limited to the Plan area, so the Plan proposes to conserve and manage linkage areas that connect to the Plan's boundaries. The intent of the MSHCP conservation strategy is to protect identified habitat within, and to maintain the functional values of, the Plan's Linkage Design (i.e., its conservation areas). Additionally, the MSHCP/NCCP has identified Linkage Extension Areas outside of its Plan Boundaries for future conservation planning. Once fully conserved, the linkage system identified by the Plan will connect 6.5 million acres of federal lands designated for conservation. We obtained spatial data layers for the Plan boundaries, Corridor, and Linkages from the Town of Apple Valley.

Covered Species
The proposed species covered by the Apple Valley MSHCP include five birds, two mammals, one reptile, and one plant ( Table 2). The federal status of species ranges from BLM-sensitive and status under review to endangered or protected under other federal legislation such as the Migratory Bird Treaty Act of 1918 [32]. The CDFW state status of these species ranged from species of special concern to endangered. One species (desert tortoise) was identified as vulnerable, and the other species were listed as least concern by the International Union of Conservation of Nature (IUCN) red list [33].

Biodiversity Metrics
GAP habitat models (1590 spp.) for the conterminous U.S. are top-down deductive models based on expert opinion, literature, and wildlife habitat relationships [18]. These data use predictive environmental variables (e.g., land cover, elevation, and distance to water) to derive habitat models for each species. These national models are in geotiff format and include 621 birds, 365 mammals, 322 reptiles, and 282 amphibians for the conterminous U.S. They were derived through a partnership between the USGS GAP, New Mexico State University, North Carolina State University, Boise State University, and EPA (EnviroAtlas). Previously, similar deductive habitat models were shown to have ≥60% concordance with actual habitat features on the ground [34]. In addition to habitat models, GAP created Land Ownership (PAD-US) and GAP Land Cover products [18,28,35,36].
Boykin et al. [13] used species richness, or number of species per pixel, to identify and create metrics that represent either biodiversity or some ecosystem service. These datasets spatially combined all 1590 species-level models or subsets of these datasets to create richness maps at a moderate spatial resolution (30 m). Derived datasets included total terrestrial vertebrate species and four major taxonomic classes: mammals, birds, reptiles, and amphibians. Additionally, species lists reflecting other important subsets were derived and metrics created that are associated with taxa richness (winter, summer birds), harvestable species, rare species, and conservation importance (Table 1). These datasets allow analysis across the conterminous U.S. As the GAP habitat models remove areas without suitable habitat (non-habitat) within known range of the species, the patterns of richness are more refined than coarse-resolution species' habitat distribution models discussed in the literature [37]. These models are created from mapped datasets associated with suitable habitat; thus, these datasets provide for analysis of biodiversity, change at ecosystem and landscape scales, and alternative futures.
The 26 national biodiversity metrics (Table 1) are useful indicators of the potential for providing three broad classes of ecosystem services: biodiversity conservation; food, fiber, and materials; and recreation, culture, and aesthetics, or their combination [38]. The metrics were reviewed by peers from state, federal, county agencies, tribes, and non-governmental organizations. Reviews and comments were solicited from scientistic peers via numerous domestic and international meetings in addition to a peer reviewed National Biodiversity Metrics and Mapping Workshop before finalizing the final metrics [38].

Desert Tortoise and Desert Kit Fox
The desert tortoise (Gopherus agassizii) and desert kit fox (Vulpes macrotis arsipus) were identified by the MSHCP team as proposed covered species of conservation interest for further assessment. Both species are wide ranging species within the MDE and have been identified as being umbrella species [39,40]. We used the individual species models for both species developed by GAP [18] to evaluate their efficacy as a conservation tool in comparison to the biodiversity metrics for the designated Plan area.

Analysis
We computed biodiversity metric values for each of the major terrestrial vertebrate taxonomic groups for four conservation zones and two specific species. The four conservation zones were the Urban Core and the Rural Development Area (Development Zone), the Apple Valley Desert Linkages zone (Linkages), Mojave River corridor zone (Corridor), and the entire MDE. The two specific species were the desert tortoise and desert kit fox (habitat vs. non-habitat). We used the Mojave Basin and Range ecoregion dataset available from the EPA [41]. We conducted all spatial analyses using ArcGIS Pro 2.4 [42] and calculated statistics using the ArcGIS Zonal Statistics Tool for each zone including the mean, standard deviation, minimum, and maximum number of species. Zonal statistic calculations are based on all grid cell values within each zone. We also calculated the Normalized Index of Biodiversity (NIB) for the MDE [13]. The NIB provides a value to facilitate comparison of metric values between study areas. NIB normalizes the mean value for each metric for a given study geography relative to the maximum value among all pixels in a reference geography, in this case the MDE. Thus, the NIB metric values ranged from 0 to 1. We then calculated the Zone:MDE ratio by dividing the mean species for each zone by the mean species for the entire MDE. Values greater than 1 indicated the zones were more species rich than the MDE and values below 1 indicate the zone is less species rich.

Results
Habitat for 329 terrestrial vertebrate species was modeled in the Apple Valley study area, including 196 birds, 75 mammals, 51 reptiles, and 7 amphibians. The Plan area (134,161 hectares) represents 1% of the MDE (12,773,240 hectares). The three PAZs represent smaller proportions with the Development zone (45,343 hectares), the Linkages zone covering 44,101 hectares, and the Corridor zone the smallest at 560 hectares.

Apple Valley in Context
Total species richness is higher in the eastern part of San Bernardino County than the western area where the Plan area is located (Figure 2). Total vertebrate species richness had a high of 175 species/pixel within the Plan area and a high of 204 species/pixel within the entire MDE ( Table 3). The NIBTotal provides values for the comparison between areas and identified that the Corridor (0.54) had a higher NIB than the MDE (0.48). The Development (0.47) and Linkages (0.45) were similar and slightly below the MDE. The Corridor zone is typical of species-rich riparian areas in the Southwest U.S. [43] and represents less than 1.0% of the Plan area. The Corridor has higher mean total species richness than the other two zones (mean 109.7; CorridorTotal: MDETotal = 1.12). The Development zone is slightly lower in mean richness than the MDE (DevelopmentTotal: MDETotal = 0.99; Table 3). The Linkages zone was the least rich (mean = 91.2; LinkageTotal: MDETotal = 0.93).
Birds were the most numerous taxa and ranged up to 95 species/pixel. Bird species richness (Figure 3b Amphibian species richness was low in each of the zones and MDE analyzed ranging from 0-4 species/pixel (Table 3). High amphibian species richness (Figure 3a)

Desert Tortoise and Desert Kit Fox
Modeled suitable habitat for the desert tortoise was identified on 66.1% (8,436,915 hectares) of the MDE compared to 33.9% (4,336,331 hectares) without habitat. Similarly, modeled habitat for the desert tortoise occurred on 68.4% (61,569 hectares) of the Plan Area compared to 31.6% (28,433 hectares) without habitat. The desert tortoise habitat covers almost the entire lowland area of the planning area apart from urban and agricultural areas and road conveyances. GAP Modeled habitat for the desert kit fox was identified on 72.2% (9,220,673 hectares) of the MDE compared to 27.8% (3,552,573 hectares) without modeled habitat. Modeled habitat for the desert kit fox was identified on 68.9% (62,047 hectares) of the Plan Area compared to 31.1% (27,956 hectares) without modeled habitat. The desert kit fox has similar overall coverage but can occupy different habitats throughout the study area and Linkages area.
The desert tortoise and desert kit fox have been identified as umbrella species [39,40]. Our analysis found their modeled habitat is pervasive across the Plan area and thus their habitat across the entire MDE nearly covers the co-existence of almost all other species associated with the study area ( Figure 4).
The Normalized Index of Biodiversity allows comparisons between areas and between metrics [13]. The NIB was higher in the areas with habitat for desert tortoise habitats and desert kit fox ( Figure 5; Table 4). This pattern was observed with all the biodiversity metrics. Similar patterns were observed at both the MDE and Plan area ( Table 4) with areas of desert tortoise or desert kit fox habitat having higher NIB than areas of non-habitat.
Comparison of species richness between desert tortoise habitat and non-habitat consistently yielded higher richness within desert tortoise habitat. For example, total species richness within the MDE averaged 107 species per pixel for desert tortoise habitat and 80 species per pixel without habitat (Table 4). Specific to the Apple Valley MSHCP Plan area, this comparison yielded higher richness in areas of desert tortoise habitat (mean = 104.8) than non-habitat (mean = 70.4) (Table 4). Similarly, for desert kit fox, total species richness averaged 108.9 species per pixel on desert kit fox habitat and 69.2 species without habitat. Within the Plan area, desert kit fox had an average of 104.6 species/pixel and non-habitat had 70.3 species/pixel. For desert kit fox, the habitat is similar to that of the desert tortoise, though a few different model variables are identified throughout the Plan area ( Figure 4). Fox habitat included higher elevations and corresponding vegetation types such as woodlands and chaparral.

Discussion
This case study was designed to test the efficacy and efficiency of incorporating biodiversity metrics into conservation planning. We used biodiversity metrics previously created from the USGS Gap Analysis Project and the EPA's EnviroAtlas [13] to test metric application specifically to the Apple Valley Multi-Species Habitat Conservation Plan. Using five biodiversity metrics, i.e., total terrestrial vertebrate, mammal, bird, reptile, and amphibian species richness, we analyzed the Apple Valley PAZs in relation to the entire MDE. We then compared these biodiversity metrics to habitat for two species (desert tortoise and desert kit fox) proposed to be covered within the MSHCP.
Comparison of biodiversity metrics in the Apple Valley PAZs including Development zone (Urban Core and Rural Development Areas), Linkages zone, and Corridor zone varied in relation to the entire MDE. NIB values for the Corridor were consistently higher than the MDE, while Development and Linkage areas were lower than the MDE. This analysis indicates that the Development and Linkage zones were less species-rich and biodiverse than the MDE as a whole and the Corridor was more species rich and biodiverse for total, bird, and amphibian metrics. This was not the case with the mammal and reptile metrics. NIB for each PAZ and MDE were similar for mammals. For reptiles, the PAZs were all higher in NIB than the MDE. Total richness was influenced by bird richness providing a high richness in the Corridor that is typical of southwestern riparian habitats even with the limited spatial representation [44].
Taxa specific outcomes varied. There was high bird species richness along the Corridor and greater richness in a few of the higher mountains (Black Mountain, and Sidewinder Mountain) in the eastern Linkages zone. Riparian areas are often focal areas as they provide water, structure, and cover [44]. Mammal species richness was similar across all zones, including the MDE. There was increased mammal species richness directly north, east, and south of the Town of Apple Valley in the Development zone. For reptiles, species richness was similar within the MSHCP geographies and slightly higher than the MDE as a whole. Increased reptile species richness was identified in the lowland areas directly north, east, and south of the Town of Apple Valley located in the Development and Linkages zones. Reptiles are well suited to the dry desert climate of the MDE. Amphibians are rare in desert regions because of the taxa's life history requirements combined with ambient arid conditions The MDE was slightly more species rich than the Development zone but numbers were low. The ranges of many amphibians within the MDE are peripheral to the Plan area. The few amphibians modeled in the Plan area were largely associated along the perennial reaches of the Mojave River (Corridor).
Only small differences were observed between metrics in the Development and Linkage zones and in some cases, the Development zone values were higher than the Linkage. These similarities are associated with similar habitats and the proximity of both zones. Except for the limited mountains, the habitat is largely associated with desert scrub. The differences, if any, between the zones may be associated with proximity of the Development zone to the more species rich corridor area.
The second objective was to compare biodiversity metrics to two species covered under the Plan (desert tortoise and desert kit fox) and evaluate their ability to provide full protection to the land use planning area. The desert tortoise and desert kit fox are wide ranging species within the MDE and are often identified as being umbrella species [39,40], a species with habitat requirements that are also important for many other species [20]. Desert tortoise use arid areas, particularly creosote bush shrublands with loamy soils while avoiding stony substrates and areas of low plant cover [45,46]. This is the only proposed covered species by the MSHCP listed as vulnerable on the IUCN list ( Table 2). The desert kit fox uses open grassland, scrublands, deserts, and salt flats [40,46]. These habitat features are widespread across the study area and the MDE in general. Our analysis supports the umbrella concept for these species: areas with either tortoise or kit fox habitat were also associated with large numbers of species. This was not just limited to either taxon conspecifics for each species, but for richness across all taxa groups and total richness.
Previous mapping efforts of desert tortoise habitat have identified critical habitat [47] or modeled tortoise species distribution using inductive modeling [48]. Critical habitat was identified adjacent to the project area but was not identified within the area. The species distribution model identified similar areas as the GAP habitat model, though it was done at a much coarser resolution. These distribution models identify suitable habitat, both within and outside the Plan area, that could be protected to maintain species.
We used the Normalized Index of Biodiversity to compare areas and metrics [13]. Our findings are consistent with the intended purpose of the Plan's Linkage Design to support biodiversity including the Corridor. Species richness was higher in the Corridor compared to the entire MDE and the Development and Linkages zones. Other studies [49,50] have used habitats or land cover types to indicate biodiversity, which could provide further context to our analysis. A focus on desert tortoise and/or desert kit fox, species with large area requirements, should benefit not only those species, but afford protection for a high number of other species that share the same habitat, thus allowing the Plan approach of using covered surrogate species for making conservation-related decisions and a simpler way to manage ecological communities [19][20][21]. Overall, the analysis is consistent with the intended purpose of the Plan for both covered species and common species.
A powerful aspect of GAP data is that it allows the ability to conduct analysis on all terrestrial vertebrate species within any geography within the conterminous U.S., and within habitat of other species such as umbrella species. Species can be selected that are identified by decision-makers and comparisons made between individual species or subsets of species. Additionally, the richness of individual land cover types can also be used as an analysis element [27]. Further, this provides the ability to assess species complementarity within study areas or geographic regions [51]. As in all models, errors can exist within the GAP habitat models [34]. Errors in individual models may be lessened because of the summation of multiple species models in creation of metrics. GAP data and mixed combinations of national biodiversity metrics have been used at multiple scales to inform decision-making from counties such as Clark County, Nevada [52] to watersheds including the South Platte, CO and Upper Rio Grande, NM [17,53] to nationwide [18,28]. The analysis can inform decisions, but care and thoughtful application is required at the implementation stage.
The ecosystem services paradigm is a useful organizing framework for characterizing and assessing biodiversity conservation. Along with the deductive modeling, it provides great utility for mapping and quantifying metrics of biodiversity conservation at multiple scales and within a reasonable timeframe [13,17,53]. In this effort, we characterized three specific conservation zones in regard to a habitat conservation plan. We worked with the decision-makers to focus our effort on layers that were deemed important to the specific effort. As the Plan is still in process, we only used our five basic richness layers in this case, but we have the ability to analyze an additional 21 metrics should the need arise ( Table 1). The ability to focus on any specific species was also important in this specific case. The EnviroAtlas can allow users to quickly review similar information [24], though the analysis is better conducted in a Geographic Information Systems platform.
Historically, conservation focused on specific species or specific geographies, i.e., forest management plans, resource management plans, and habitat conservation plans. Data are now available to understand those species and their habitats in context of larger regions and with surrogates for biodiversity. These analyses provide a greater understanding of the overall system for the decision makers. These analyses with the understanding of the value of biodiversity, both in ecological and economic terms, can provide an important step forward in conservation in the context of urbanization [12,13,23,24].

Conclusions
Our analysis provides ecological context for the Apple Valley MSHCP. Decisionmakers were interested in understanding how the five specific metrics (total vertebrate, amphibian, bird, mammal, and reptile species richness) used would reflect the importance of the area within the larger Mojave region. Our analysis provides the context for all three study area zones (Development, Linkages, and Corridor) within the larger Mojave Desert Ecoregion. Our approach incorporates built-in extensibility so that as other metrics are identified or if alternative plans are developed, they can be explored and added to the system. For the Apple Valley MSHCP, our analysis was able to efficiently identify how each of our biodiversity metrics were related to the proposed MSHCP landscapes and provide additional information beyond the few covered species within the Plan. This provides an expanded analysis to identify any further potential conflicts in species conservation within the area.
The analysis described in this paper provides baseline or reference condition for considering and comparing alternative future scenarios (e.g., climate change, urbanization, and trend analysis). It also establishes common sense indicators of ecosystem services for national, regional, tribal, and local end user and decision-maker needs. It is flexible enough to add and test new metrics as they are identified and can now be undertaken anywhere within the conterminous U.S.
Although this process has been applied to Apple Valley and other geographies within the U.S., the approach has practical application for other global biodiversity initiatives, such as the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES). IPBES, for example, is designed to provide consistent indicators that can relate science and policy for the purpose of conducting periodic assessments of biodiversity and ecosystem services at global, regional, and subregional scales [54]. The approach and metrics (Table 1) utilized for the Apple Valley MSHCP could equally be applied to other non-U.S. geographies providing that common technology and spatial datasets (e.g., species distribution models and digital land cover), in addition to local expertise, are available to organize the data into relevant functional groups (biodiversity metrics) using stakeholder involvement.

Institutional Review Board Statement: Not Applicable
Data Availability Statement: Publicly available datasets were analyzed in this study. The individual species models used to generate the richness for this manuscript are available at https://gapanalysis.usgs.gov/apps/species-data-download/ (accessed on 21 July 2021) and through the science base repository (https://doi.org/10.5066/F7V122T2) (accessed on 21 July 2021).