One goal of Habitat Conservation Plans is to protect viable populations of animal species. Management actions included in such plans often focus on vegetation restoration to benefit the target animal species. Yet, such restoration activities are often uninformed by fine-scale animal survey data. Using the California Gnatcatcher (*Polioptila californica californica*), we demonstrate how survey data can guide restoration toward the goal of improving gnatcatcher viability by identifying habitat conditions most favorable for gnatcatcher occupancy. We found that gnatcatcher presence and colonization probabilities tripled as the coverage of coastal sagebrush (*Artemisia californica*) increased from 10% to 40%. Coastal sagebrush was more likely to be present at low elevations, and was most closely associated with soil texture of 5%–20% clay, 40%–70% sand, and 20%–40% silt. Higher gnatcatcher extinction probabilities were associated with closed, dense habitat; optimal conditions were between 30% and 40% openness, creating a slightly closed canopy. Open habitat was associated with southern aspects, shallow slopes, and inland areas. An understanding of the soil types and physical parameters of the environment that affect vegetation, especially the amount of coastal sagebrush needed for high gnatcatcher occupancy and colonization rates, will help to define restoration target conditions. Using multiseason, occupancy-based surveys in conjunction with detailed habitat measurements will allow ornithologists and land managers to design powerful restoration prescriptions, even within narrowly defined ecosystems.

## INTRODUCTION

Implementation of Habitat Conservation Plans is commonly used to aid the recovery and persistence of species of concern (Kareiva et al. 1999). Monitoring of animal population trends is required to assess the effectiveness of these plans (Noss et al. 1997, Case and Fisher 2001). Even so, a disconnect exists between an understanding of the indirect effects of environmental factors that influence vegetation communities and the optimization of conditions for target species. Although habitat suitability models developed using broad-scale data, such as remote sensing, are useful for planning purposes, restoration projects, where the goal is the recovery of populations, are most effective when models are refined using covariates that vary at a smaller scale relative to species occupancy patterns in the field. For example, breeding sites selected for Black Terns (*Chlidonias niger*) using only landscape suitability indices had <20% occupancy (Wyman and Cuthbert 2016). In contrast, Taillie et al. (2015) used both local and landscape level predictors to explain the distribution patterns of Bachman's Sparrows (*Peucaea aestivalis*), and showed that small prescribed burns to control the density of high grass increased the predicted occupancy probability to 75%. This disconnect of scale between landscape-level planning and site-specific restoration often leaves a wide gap between vertebrate conservation and ecological restoration (Cabin 2007, Cabin et al. 2010, DeSimone 2013). Our goal was to lessen this gap by demonstrating an approach that provides guidance for restoring vegetation to benefit animals covered under conservation plans. As a case study, we focused on Habitat Conservation Plans for southern California that aim to preserve or restore habitat for the coastal California Gnatcatcher (*Polioptila californica californica*).

In southern California, the gnatcatcher serves as an umbrella and flagship species for many Habitat Conservation Plans (R.J. Meade Consulting 1996, San Diego County 1998, Dudek Consulting 2006). Past work (Winchell and Doherty 2008, 2014, TAIC 2002) has studied gnatcatcher occupancy as a function of climatic and landscape variables, but this work is of limited utility to restoration ecologists because of the coarse scale of inference. Nonetheless, from these studies we understand that gnatcatcher distribution west of the Peninsular Range Provinces (coastal mountains) is limited by steep slopes, too much precipitation, low temperatures, and high elevations. Thus, colonization is greater on relatively gentle slopes with dry, warm climates that occur mainly at comparatively low elevations. In addition, gnatcatcher occupancy has been found to be strongly dependent on the presence of coastal sagebrush (*Artemisia californica*; Bontrager 1991, Beyers and Wirtz 1997, Atwood et al. 1998).

Preserving and restoring coastal sagebrush has been a key component of actions to improve the viability of the gnatcatcher in southern California, but gnatcatcher monitoring activities have not been directly linked to coastal sagebrush restoration activities (e.g., Winchell and Doherty 2008, 2014). Questions about the optimal spacing and density of coastal sagebrush for gnatcatcher presence, whether nonnative plants affect gnatcatcher presence, and what vegetation characteristics are related to gnatcatcher occupancy, colonization, and extinction all need answers. Identifying suitable locales for restoring coastal sagebrush, and thus gnatcatchers, is the next step. In particular, identifying combinations of aspect, slope, and soil type that are most likely to produce optimal vegetation conditions can focus restoration. Our goal was to guide restoration at a scale relative to habitat usage by nesting pairs of gnatcatchers. We measured biological and environmental parameters to determine what covariates maximized sustained occupancy by gnatcatchers. Our objectives were to analyze gnatcatcher occupancy dynamics (occupancy, local colonization, and extinction) in relation to habitat covariates (e.g., vegetation and soil types) at a finer scale than has been done with previous models. This information can guide land managers toward more effective restoration and maintenance of gnatcatcher habitat.

## METHODS

We worked within San Diego County, California, USA, on lands designated as public or quasi-public, and on some military lands (Figure 1). In general, the area included preserve lands regulated under Habitat Conservation Plans, other public lands designated as open space and regulated by local jurisdictions, and lands separately managed by Marine Corps Air Station Miramar or Naval Weapons Station Seal Beach, Detachment Fallbrook. All lands had some component that was distinctly managed for gnatcatcher conservation. Marine Corps Base Camp Pendleton, private lands, and Tribal lands were excluded from the sampling frame because of access limitations, logistic difficulties, and liability concerns.

We collected data on gnatcatchers and their habitat in 2004, 2007, and 2009 as part of the San Diego Multiple Species Conservation Program (San Diego County 1998). Each study year recorded below average rainfall for San Diego County. During our entire survey time frame, only 2005 had above average rainfall. Because we sought to use fine-scale data for our occupancy analyses, we focused on areas that were within strata labeled as high or very high quality by TAIC (2002), were generally within the preferred climatic and landscape conditions identified by previous analyses (Winchell and Doherty 2014), and were located where gnatcatchers were more likely to occur. These strata represented habitat covering 22,230 ha. In 2007, our study area experienced a massive wildfire, reducing our area to 14,150 ha and our sample size to 76 survey point locations that were unburned and thus still likely to have gnatcatchers. Bird survey details have been previously described (Winchell and Doherty 2008, 2014); in brief, point locations were surveyed for gnatcatchers for 15 min during 5–6 visits each year during the breeding season. We included gnatcatcher observations within 50 m of survey points, with distance determined using a rangefinder. From these data, multiple season robust design occupancy encounter histories were constructed (MacKenzie et al. 2003), from which gnatcatcher occupancy, local colonization, local extinction, and detection probabilities were estimated.

In 2007, we collected habitat data along two 50-m transects centered on each of the 76 survey points. The first transect was along a random azimuth, establishing a 50-m radius, whereas the second transect was along a random azimuth at least ±45° from the first transect. We conducted habitat measurements every 2 m along each transect for a total of 50 sampling locations associated with each gnatcatcher survey point. At each of the 50 locations, we recorded the presence of any vegetation, as well as the presence of individual grass, herb, shrub, and tree species. We summarized these data as the proportion of the 50 locations at which no plant species occurred, the proportion at which a particular plant species or growth form (grass, herb, shrub, tree) occurred, as well as whether a plant species was detected at any of the locations. We also collected a soil sample at each survey point, from which we derived the percentage of clay, sand, and silt (Bouyoucos 1962, Gee and Bauder 1979). Additionally, we collected soil samples at the end of each transect to determine whether soil texture was homogeneous within our 50-m radius area.

### Predictions and Analysis

We first tested predictions concerning gnatcatcher occupancy, colonization, extinction, and detection probabilities. Specifically, we thought that gnatcatcher occupancy and colonization probabilities would be positively associated with, and extinction probability would be negatively associated with, the proportion of coastal sagebrush coverage (ARCA; Figure 2). We also anticipated an optimal intermediate level of coastal sagebrush coverage, so we included a quadratic term for the proportion of coastal sagebrush (ARCA^{2}; Figure 2). Similarly, we predicted that higher proportional coverage of invasive nonnative grasses or herbs, such as black mustard (*Brassica nigra*) and stork's bill (*Erodium* spp.), would be negatively associated with gnatcatcher occupancy and colonization probabilities, and positively associated with extinction probability, because invasives may compete for water with germinating native plants and may provide cover for predators. We predicted that higher proportional coverage of black sage (*Salvia mellifera*; SAME) would be negatively associated with gnatcatcher occupancy and colonization probabilities but positively associated with extinction probability were gnatcatchers to be averse to using black sage due to its dense structure and limited forage (Robinson and Holmes 1982, Roach 1989). We predicted that the presence of laurel sumac (*Malosma laurina*; MALA) would be positively associated with gnatcatcher occupancy and colonization, but negatively with gnatcatcher extinction, because laurel sumac is frost intolerant (Riggan et al. 1994) and gnatcatchers are sensitive to temperatures below 4.5°C (Root 1988, Mock 1998); thus, we viewed MALA as a fine-scale indicator of temperatures. We predicted that gnatcatchers might prefer more openness (proportion of locations with no trees or shrubs; Atwood and Bontrager 2001) in their territories (Figure 2). Thus, we thought that openness would be positively associated with occupancy and colonization probabilities and negatively associated with extinction probability, as gnatcatchers frequently forage for tiny insects by flying between shrubs rather than hopping from shrub to shrub, and can probably not forage as efficiently in dense vegetation. None of these predictor variables were correlated (*r* > |0.47|), and, since all were proportions or binomial, we did not rescale the variables prior to analysis. Based on past analyses (Winchell and Doherty 2014), we predicted that detection probability would be constant among years (*t*); we tested this prediction by including constant and annually variable detection structures in our modeling efforts.

We incorporated the above predictions into models using the dynamic occupancy model of Mackenzie et al. (2003) in program MARK (White and Burnham 1999). This model permits estimates of occupancy, colonization, extinction, and detection probabilities. If occupancy at time *t*, colonization, and extinction are estimated, then occupancy at time *t* + 1 can be derived. Thus, we estimated the initial occupancy probability (in 2004), as well as time-specific colonization and extinction probabilities, and derived future occupancy estimates. Our surveys did not occur at equal intervals (i.e. 2004, 2007, and 2009), so we standardized the modeling and estimation of colonization and extinction between surveys on an annual basis. We had to assume that annual colonization and extinction probabilities were identical within survey intervals because of the lack of intervening survey data. For example, we had to assume that the colonization probabilities from 2007 to 2008 and from 2008 to 2009 were the same. Considering the above parameters, hypotheses, and limitations, we constructed 92 submodels for colonization, 92 submodels for extinction, 46 submodels for occupancy, and 2 submodels for detection.

Model building and selection is challenging when using models with numerous submodels (e.g., occupancy, colonization, extinction, and detection) and many covariates due to the large number of possible models (Doherty et al. 2012, Bromaghin et al. 2013). To reduce the number of models, we used a strategy similar to that of Bromaghin et al. (2013) by focusing on a set of plausible submodels and variables and combining these submodels in all possible combinations to examine the effects of the remaining covariates. We first focused on each parameter (i.e. occupancy, colonization, extinction, and detection) by holding the other 3 parameters constant. This procedure forced the maximum amount of explainable variation into the focal parameter. For each focal parameter, we constructed models representing all possible additive combinations of predictor variables for a balanced model set and calculated cumulative Akaike model weights for each predictor variable (Burnham and Anderson 2002). Barbieri and Berger (2004) implied that variables with >0.5 cumulative Akaike weights have meaningful relationships with parameters of interest; at this stage in our analysis, we used a more liberal cutoff of 0.4 cumulative Akaike weight to avoid eliminating variables. Using the reduced set of predictor variables for each parameter, we then constructed all possible additive combinations of models across all parameters for a final, balanced model set and used Akaike's information criterion with small sample size correction (AIC_{c}) for model selection and for calculating cumulative Akaike weights for each predictor variable. Following Barbieri and Berger (2004), and maintaining our liberal 0.4 cumulative Akaike weight cutoff for consistency, we defined the model that contained all independent variables with ≥0.4 cumulative Akaike weight as our predictor model from which to produce estimates.

The above analyses identified which vegetation parameters were correlated with gnatcatcher occupancy, colonization, and extinction probabilities. Understanding how vegetation is correlated with soil type and other variables can help to identify ideal locations for restoration, so we focused on explaining variation in the presence of coastal sagebrush, laurel sumac, and openness (i.e. lack of trees or shrubs). We focused on the 5 vegetation sampling locations along each transect nearest to our center-point soil sample, instead of all sampling locations, because we found that soil texture differed at the end of at least one transect at 17% of our survey points. We viewed the 10 locations as repeated vegetation measures associated with each center point and used a repeated measures logistic regression in which the dependent variable was the presence or absence of the plant species of interest. As independent variables, we considered soil texture, aspect, slope, elevation, and distance to the coast of the center point. Because soil texture has 3 components (percent clay, percent sand, percent silt) that sum to 100%, the 3 values are not independent. We followed Aitchison's (1986) advice for analyzing such compositional data: We used a log-ratio transformation of percent clay and percent sand relative to percent silt because percent silt showed the least variation across our sampled points. Aspect (degrees) can influence vegetation through the amount of solar radiation (Kirkby et al. 1990), so we transformed our aspect measures by taking the cosine of the radian to address the circular nature of the degree measurement, thus equalizing the gradient in aspect in both directions of north (Beers et al. 1966, McCune and Keon 2002). We thought that elevation (m), slope (degrees), and distance to the coast (m; because of the fog that can develop) could predict where particular plant species grow. We constructed all possible additive combinations of predictor variables and relied upon AIC_{c} for model selection and for calculating cumulative Akaike weights for each predictor variable. We report estimates and standard errors (SE) unless otherwise noted.

## RESULTS

Using gnatcatcher survey data (summarized in Appendix Table 4), our first set of analyses focused on each parameter (occupancy, colonization, extinction, and detection) separately and helped to reduce the number of predictor variables for further modeling. We identified ARCA, ARCA^{2}, openness, and MALA as relevant for further analysis of occupancy; ARCA and year for analysis of colonization; ARCA, openness, MALA, and year for analysis of extinction; and year for analysis of detection (Table 1). We then constructed all possible combinations of these variables for our final model set and calculated cumulative Akaike model weights for each variable (Table 2). Influential predictors were ARCA, ARCA^{2}, openness, and MALA for occupancy; ARCA and year for colonization; and openness and MALA for extinction (Table 2).

## TABLE 1.

Cumulative Akaike model weights for each variable used to model the occupancy, colonization, extinction, and detection rates of California Gnatcatchers in San Diego County, California, USA, in relation to habitat parameters. We constructed all possible additive combinations of variables for each parameter while holding the other 3 parameters constant. Variables with cumulative weights >0.4 were used in further modeling. ARCA is the proportion of 50 transect intersections (‘locations') associated with each gnatcatcher survey point (‘point') with coastal sagebrush (*Artemisia californica*) as a linear term, and ARCA^{2} is a quadratic model; Openness is the proportion of locations at a point with no vegetation; Invasive is the proportion of locations at a point with nonnative grasses (black mustard [*Brassica nigra*] or stork's bill [*Erodium* spp.]); SAME is the proportion of locations at a point with black sage (*Salvia mellifera*); MALA is the proportion of locations at points with laurel sumac (*Malosma laurina*) present (0 = absent, 1 = present). The range of values appears in parentheses after the variable name. Asterisks (*) denote that variables were not considered for a particular parameter. Appendix Table 5 shows the top 10 models for each parameter.

## TABLE 2.

Cumulative Akaike model weights for each variable used to model the occupancy, colonization, extinction, and detection rates of California Gnatcatchers San Diego County, California, USA, in relation to habitat conditions, from models with all possible additive combinations of the variables. See Table 1 for variable codes. Asterisks (*) denote that variables were not considered for a particular parameter. The model containing all variables with ≥0.4 cumulative Akaike weight was the predictor model used in Figures 2–4. Appendix Table 6 shows the models with ΔAIC_{c} < 2.

We found that gnatcatcher occupancy probability increased with proportional coastal sagebrush coverage (Figure 3), with >30% coastal sagebrush coverage resulting in high occupancy probabilities. Occupancy was slightly negatively associated with the proportion of openness (i.e. lack of trees and shrubs; β = −1.95 ± 1.46) and with the presence of laurel sumac (β = −0.92 ± 0.68). Colonization probability increased with proportional coastal sagebrush coverage (Figure 4). Gnatcatcher colonization probability was greater in 2007–2009 than in 2004–2007. Extinction probability was negatively correlated with openness (Figure 5) and increased with the presence of laurel sumac (β = 1.16 ± 0.76). The detection probability of gnatcatchers was 0.42 ± 0.03 per survey visit.

Following the above analysis, we investigated how the presence of coastal sagebrush, openness, and laurel sumac depended on soil texture, aspect, slope, elevation, and distance to the coast. We found that soil texture and elevation were useful predictors of the amount of coastal sagebrush coverage. Aspect, slope, and distance to the coast all had cumulative Akaike weights greater than 0.40 for predicting openness. None of our variables were useful predictors of laurel sumac presence (Table 3).

## TABLE 3.

Cumulative Akaike model weights for independent variables used to investigate the presence of coastal sagebrush (*Artemisia californica*; ARCA), laurel sumac (*Malosma laurina*; MALA), and the absence of trees and shrubs (Openness) at California Gnatcatcher survey points in San Diego County, California, USA. We constructed all possible additive combinations of independent variables for each dependent variable for a balanced model set. The range of values appears in parentheses after each variable name. Appendix Table 7 shows the top 10 models for each dependent variable.

We primarily found coastal sagebrush in soil textures ranging from 5% to 20% clay, 40% to 70% sand, and 20% to 40% silt. These soils were more likely to be present at low elevations, although the relationship with elevation was not strong (β = −0.003 ± 0.003). Openness (absence of trees and shrubs) declined on steep slopes (Figure 6) and had a weak positive relationship with aspect (β = 0.305 ± 0.248). Upon further examination, distance to the coast was found to be an uninformative parameter (Arnold 2010; β = 0.269E-4 ± 0.201E-4).

## DISCUSSION

Our results support the observation that coastal sagebrush strongly influences gnatcatcher presence (Figure 3; Atwood and Bontrager 2001). In turn, coastal sagebrush presence is influenced by soil texture and elevation. We also found that gnatcatcher colonization probabilities were most closely associated with coastal sagebrush (Figure 4). These results give some insight into the dynamic processes underlying gnatcatcher population dynamics. We found that gnatcatcher occupancy (Figure 3) and colonization (Figure 4) increased with coastal sagebrush coverage and were highest with coastal sagebrush coverage of 30%–40% of the plant community, the maximum values recorded at our sites. Restorationists might try to artificially increase coastal sagebrush coverage above 40% through plant pallet composition and irrigation, but we caution against such an approach as this may not be a situation found in natural communities (Sproul et al. 2011). Instead, coastal sagebrush coverage of at least 30% might be a good target for restoration activities. In addition, we found that coastal sagebrush presence was most closely associated with a soil texture of low clay, moderate sand, and moderate silt, a type of soil most likely to be found at low elevations.

Identifying highly suitable locations for sagebrush and testing soil texture prior to restoration will likely lead to increased probabilities of coastal sagebrush restoration success. Gnatcatcher extinction probability was highest in dense vegetation (Figure 5). Restoration activities close to the coast and on slopes less than 20° (Figure 6) are likely to lead to vegetation characteristics preferred by gnatcatchers. Sites with these characteristics have extinction probabilities <50%. We recommend that land managers evaluate the physical parameters of the ecosystem, namely soil, slope, and elevation, and use these factors to predict where coastal sagebrush will best thrive and thereby provide habitat for gnatcatchers. Although we focused on occupancy and related parameters to define important habitat used by birds, other demographic parameters such as abundance, survival, and reproduction should also be considered. However, demographic parameters are often difficult to estimate. In the absence of such parameter estimates, especially when management funds are limited and disjointed, we believe that relying on occupancy-related parameters can lead to successful restoration.

One of our original predictions was that laurel sumac would be positively correlated with occupancy because gnatcatchers are frost intolerant. We found the opposite: Gnatcatcher extinction probability increased with the presence of laurel sumac. This result may have been a factor of scale and location. Gnatcatcher occupancy increased closer to the coast. Near the coast, laurel sumac is located on moist, cool, northern-aspect slopes (Holland 1986, Bartolome et al. 2007) and is a primary component of chaparral-type communities, where gnatcatchers have low occupancy probabilities (Winchell and Doherty 2008). Inland, laurel sumac sparsely populates midelevation valleys where coastal sage scrub occurs and frost is a rare event. Laurel sumac may be a good indicator of gnatcatcher habitat in areas away from the coast; however, such areas are limited as interior land was converted to avocado farms early in the 20^{th} century (Coit 1963, Shepherd and Bender 2001).

The composition of other shrubs may not be as significant as the presence of coastal sagebrush for gnatcatcher habitat. We found that our predictions pertaining to invasive plants and black sage were not supported, perhaps because other plants in the communities dominated by coastal sagebrush can be quite variable For example, Westman (1981) found a high degree of variability in species diversity in coastal sage scrub plant communities. Eastern Mojave buckwheat (*Eriogonum fasciculatum*), the second most dominant plant in coastal sage scrub (Roach 1989), is not necessarily present in sites used by gnatcatchers (Atwood 1990). However, all of Atwood's (1990) sites contained coastal sagebrush. Decreased diversity of perennials may increase nest success, at least within plots with a certain amount of perennial heterogeneity (Braden et al. 1997).

Although we found that gnatcatcher occupancy probability slightly increased as the proportion of openness decreased and cover became denser, a stronger relationship existed with extinction probability: Gnatcatcher extinction probability was highest in the least open sites, and exceeded 50% with openness <30% (Figure 5). Bontrager (1991) and Braden et al. (1997) reported that areas of dense vegetation had fewer gnatcatchers, and that fledging success increased with increasing cover of grasses and forbs. For restoration, it may be important to balance extinction probability and occupancy probability with a shrub density between 30% and 40% openness, maintaining a slightly closed canopy. Managing openness requires constant maintenance, such as cutting or prescribed fire. Once the desired shrub density is established, we recommend maintaining 40% coastal sagebrush composition within the shrub community.

Coastal sage scrub is associated with steeper topography, with grasslands in flatter areas (Westman 1982, Bartolome et al. 2007). This relationship may explain the association that we found between slope and openness, wherein steeper slopes were less open. Areas with slopes >20° had <50% probability of being open (Figure 6). None of our sites had slopes greater than 35%; steep slopes in Mediterranean climates tend to be highly erodible (Morgan 2005) and unable to support stable plant communities such as coastal sage scrub. We did not sample these steep slopes, as they did not include scrub vegetation.

Restoration projects in coastal sage scrub habitats should start with targeting loam or sandy loam soils. Loamy sands should be given a lower priority and clay soils (clay loam, sandy clay loam, and sandy clay) should generally be avoided. Next, low-elevation sites with southern exposures should be given priority. Sites with elevations <100 m have a >50% chance of coastal sagebrush being present. If low-elevation sites, which are typically nearer to the coast, are not available, then sites with sparse laurel sumac could be restored, as gnatcatchers are limited by winter temperatures below 1.1°C (Root 1988). Finally, areas with greater topographic relief should be selected over flatter areas because steeper areas are associated with coastal sage scrub. Communities dominated by coastal sagebrush should be a goal of restoration for gnatcatchers. Areas should be managed to maintain 30%–40% coastal sagebrush coverage, thus revegetative plantings should be dominated by coastal sagebrush, but include a diversity of other shrubs found in occupied sites, such as Eastern Mojave buckwheat, torhleaf goldeneye (*Viguiera laciniata*), common deerweed (*Lotus scoparius*), coyotebrush (*Baccharis pilularis*), laural sumac, or black sage.

Habitat Conservation Plans (HCP) are predicated on maintaining viable populations of animal species of concern, with management actions focused on supporting this goal. In southern California, HCPs have been established to maintain flagship or umbrella species that represent the rich biodiversity of the California Floristic Province (Myers et al. 2000). As in many other systems, vegetation restoration is a common management action, but is seldom informed by animal survey data. Our results demonstrate how occupancy surveys combined with vegetation and abiotic measurements can inform both site selection and management actions that may be beneficial to California Gnatcatchers and potentially the diverse suite of species under the umbrella of their HCP.

## ACKNOWLEDGMENTS

We thank the >50 people who were involved with this project, unselfishly giving their time and talents.

**Funding statement:** The majority of this work was funded under San Diego Association of Governments Agreement No. 5000793, with funds provided via local TRANSNET tax. The findings and conclusions in this article are those of the author(s) and do not necessarily represent the views of the U.S. Fish and Wildlife Service (USFWS).

**Ethics statement:** The authors complied with USFWS 212 FW 1 Standards of Conduct while conducting all aspects of this work.

**Author contributions:** Both authors contributed to all aspects of the study.

**Data deposits:** Data are deposited with the California Department of Fish and Wildlife BIOS database, USGS regional database, San Diego Association of Governments website, and the USFWS Carlsbad Fish and Wildlife Office. Office data files are available for public access upon request.

## LITERATURE CITED

*Polioptila californica*). Manomet Bird Observatory, Manomet, MA, USA. Google Scholar

*Polioptila californica*), version 2.0. InThe Birds of North America ( A. F. Pooleand F. B. Gill, Editors). Cornell Lab of Ornithology, Ithaca, NY, USA. https://doi.org/10.2173/bna.574 Google Scholar

*Polioptila california*) in south Orange County, California. Unpublished technical report prepared for the Santa Margarita Company, Rancho Santa Margarita, CA, USA. Google Scholar

*Polioptila melanura californica*. M.S. thesis, San Diego State University, San Diego, CA, USA. Google Scholar

*Polioptila californica californica*). The Condor: Ornithological Applications 116:538–545. Google Scholar

*Chlidonias niger*) in the U.S. Great Lakes region. The Condor: Ornithological Applications 118:613–623. Google Scholar

## Appendices

## APPENDIX TABLE 4.

Raw detection data (0 = not detected, 1 = detected) for coastal California Gnatcatchers (*Polioptila californica californica*; CAGN) in 2004, 2007, and 2009 at survey points in San Diego County, California, USA, along with the associated Proportion Openness (the proportion of vegetation transect intercepts [‘locations'] at a gnatcatcher survey point [‘point'] with no vegetation); Proportion invasive (the proportion of locations at a point with nonnative grasses, i.e. black mustard [*Brassica nigra*] or stork's bill [*Erodium* spp.]); Proportion ARCA (proportion of locations at a point with coastal sagebrush [*Artemisia californica*]); Proportion SAME (proportion of locations at a point with black sage [*Salvia mellifera*]); and Presence of MALA (indicates whether laurel sumac [*Malosma laurina*] was present at any locations associated with points).

## APPENDIX TABLE 5.

Table of results for the first stage of model construction and selection to determine the influence of habitat variables on California Gnatcatcher occupancy, colonization, extinction, and detection in coastal California, USA. Each parameter (occupancy, colonization, extinction, and detection) was first considered separately while keeping the other 3 parameters constant. Variables included in models were: ARCA (the proportion of 50 transect intersections (‘locations') associated with each gnatcatcher survey point (‘point') with coastal sagebrush [*Artemisia californica*] as a linear term), and ARCA^{2} as a quadratic model; Openness (the proportion of locations at a point with no vegetation); Invasive (the proportion of locations at a point with nonnative grasses, black mustard [*Brassica nigra*] or stork's bill [*Erodium* spp.]); SAME (the proportion of locations at a point with black sage [*Salvia mellifera*]); MALA (the proportion of locations at a point with laurel sumac [*Malosma laurina*]); and Year. For each parameter, we show the top 10 models as well as the intercept only model and the global model. For each model, we show the difference in Akaike's information criterion with small sample size correction from the focal model and the top model (ΔAIC_{c}), the Akaike weight (*w _{i}*), the model likelihood, the number of parameters in the model (

*K*), and the deviance associated with each model.

## APPENDIX TABLE 6.

Models with ΔAIC_{c} < 2 from the second stage of model construction and selection in which the variables and parameters from the first step were combined in all possible additive combinations to determine the influence of habitat variables on California Gnatcatcher occupancy, colonization, extinction, and detection in coastal California, USA. Parameters included occupancy (ψ), extinction (ϵ), colonization (γ), and detection (*p*). Variables included ARCA (the proportion of 50 transect intersections (‘locations') associated with each gnatcatcher survey point (‘point') with coastal sagebrush [*Artemisia californica*] as a linear term), and ARCA^{2} as a quadratic model; Openness (the proportion of locations at a point with no vegetation); Invasive (the proportion of locations at a point with nonnative grasses, black mustard [*Brassica nigra*] or stork's bill [*Erodium* spp.]); SAME (the proportion of locations at a point with black sage [*Salvia mellifera*]); MALA (the proportion of locations at a point with laurel sumac [*Malosma laurina*]); and Year. For each model, we show the difference in Akaike's information criterion with small sample size correction from the focal model and the top model (ΔAIC_{c}), the Akaike weight (*w _{i}*), the model likelihood, the number of parameters in the model (

*K*), and the Deviance associated with each model.

## APPENDIX TABLE 7.

Top 10 models from a repeated measures logistic regression used to investigate the presence of coastal sagebrush (*Artemisia californica*; ARCA), laurel sumac (*Malosma laurina*; MALA), and the absence of trees and shrubs (Openness) as a function of soil texture, aspect, slope, elevation, and distance to the coast. For each model, we show the difference in Akaike's information criterion with small sample size correction from the focal model and the top model (ΔAIC_{c}), the Akaike weight (*w _{i}*), the model likelihood, the number of parameters in the model (

*K*), and the Deviance associated with each model.

*Polioptila californica californica*)," The Condor 120(3), 581-595, (11 July 2018). https://doi.org/10.1650/CONDOR-17-221.1