Evaluating the effects of Natural Resources Conservation Service project implementation on the disturbance-dependent avian community with implications for Blue-winged Warblers

. The Blue-winged Warbler ( Vermivora cyanoptera ) is a songbird that breeds in eastern deciduous forests of North America. The species is declining, partially due to declines in forest disturbances. According to the umbrella species concept, management actions implemented to benefit other critically declining disturbance-dependent species like the Cerulean ( Setophaga cerulea ) and Golden-winged ( Vermivora chrysoptera ) warblers should positively affect Blue-winged Warbler site occupancy and species richness of shrubland and grassland birds. Similarly, determining if the umbrella concept is supported by relating species richness of disturbance-dependent avian guilds would support continued funding for species-specific conservation and management. Our goal was to evaluate if Natural Resources Conservation Service (NRCS) projects in West Virginia implemented for Cerulean and Golden-winged warblers also positively affected Blue-winged Warbler site occupancy and the disturbance-dependent avian community. We hypothesized that Blue-winged Warbler single-season occupancy and species richness for shrubland and grassland bird species would be greater on treated sites than on untreated sites. We also included other vegetation variables (i.e., percent cover of grasses, forbs, etc.) and spatial variables (i.e., elevation (m), ecoregion, etc.) that could affect Blue-winged Warbler site occupancy. We conducted point count surveys at 341 total locations distributed among 20 private properties managed for Golden-winged Warblers (n = 147); 19 private properties managed for Cerulean Warblers (n = 197); and two properties managed for both species during 2019–2020. Treatments included a variety of management practices (i.e., brush management) following specific guidelines to improve Cerulean and Golden-winged warbler habitat. We identified and defined untreated sites as either pre-treatment sites with planned management that had not yet occurred, or as reference sites, which were outside of treatment areas and representative of pre-treatment vegetation structure conditions. Contrary to our hypotheses, treated points had lower Blue-winged Warbler site occupancy than untreated points by 34–44% depending on ecoregion (Central Appalachians, Ridge and Valley, Western Allegheny Plateau), and shrubland and grassland avian guild richness were not different at untreated and treated locations. Thus, NRCS conservation project implementation for Cerulean and Golden-winged warblers did not meaningfully affect Blue-winged Warbler site occupancy or associated shrubland and grassland bird avian richness. We detected Blue-winged Warblers across the range of elevations surveyed (244–917 m), suggesting that their breeding distribution is continuing to expand into higher elevations in the Central Appalachians. Additionally, Blue-winged Warbler site occupancy was positively correlated with shrubland conditions within 100 m of survey points and decreased with increasing basal area within 100 m of survey points. Thus, management that increases the amount of shrubland in the Central Appalachians has potential as a conservation action to benefit Blue-winged Warbler site occupancy.


INTRODUCTION
Through the umbrella species concept, habitat management for focal species has been proposed as a mechanism to indirectly benefit co-occurring species and make better use of limited conservation funding (Roberge and Angelstam 2004). Habitat management efforts are underway in priority areas of the Appalachian Mountains region to benefit Golden-winged (Vermivora chrysoptera) and Cerulean (Setophaga cerulea) warblers through Natural Resources Conservation Service (NRCS) projects (Appendix 1). Management actions implemented in these projects follow science-based Goldenwinged and Cerulean warbler habitat management guidelines (Bakermans et al. 2011, Roth et al. 2012, Golden-winged Warbler Working Group 2013, Wood et al. 2013, which may benefit cooccurring species. Many disturbance-dependent wildlife species in eastern North America have experienced population declines resulting from land use change and suppression of disturbance regimes (Hunter et al. 2001, King andSchlossberg 2014). Prior to European settlement, unique vegetation conditions in eastern North America were created via disturbances such as natural fire, intentional fire by Native Americans, windthrow, flooding, and activities by other wildlife species (e.g., North American beaver [Castor canadensis], elk [Cervus canadensis], and American bison [Bison bison]; Lorimer 2001). While the utility of the umbrella species concept as a conservation tool has been criticized, an opportunity exists to determine if Golden-winged and Cerulean warbler-specific habitat management on private lands is positively affecting other declining disturbance-dependent species, such as the Blue-winged Warbler (Vermivora cyanoptera).
The Blue-winged Warbler, a disturbance-dependent generalist (Confer and Knapp 1981), is a declining and understudied North American songbird (Gill et al. 2020). Besides studies on hybridization and competition with the imperiled Golden-winged Warbler (Buehler et al. 2007, Confer et al. 2003, Shapiro et al. 2004, little published research exists on Blue-winged Warbler breeding habitat characteristics or response to habitat management. Their breeding habitat has been described as a dense matrix of herbaceous, shrub, and sapling growth near forest edges where woody vegetation dominates (Confer and Knapp 1981, Askins et al. 2007, Confer et al. 2003, Gill et al. 2020. Landscape-scale forest clearing in the United States during the late 18 th and early 19 th century for timber and agriculture, likely benefited Blue-winged Warblers (Johnson and Govatski 2013). Now they breed in a broad range of early to midsuccessional vegetation communities predominantly resulting from anthropogenic disturbances such as utility corridors, reclaimed surface mines, abandoned agricultural fields, and clearcut harvests in early stages of succession (Confer and Knapp 1981, Confer et al. 2003, Gill et al. 2020. While Blue-winged Warbler breeding distribution has continued to expand northward and into higher elevations in the Appalachians (Canterbury et al. 1993, Patton et al. 2010, Bailey and Rucker 2021, populations declined at a rate of -0.76%/year in West Virginia during 1966-2019 due to habitat loss on breeding, migration, and wintering grounds (Rosenberg et al. 2016, Gill et al. 2020, Sauer et al. 2020. They undoubtedly benefit from habitat management techniques used to create or maintain shrubland and young forest communities (e.g., heavy timber harvest, grazing, prescribed fire, and old field management), however, their response to these techniques has not been well quantified. While Blue-winged Warbler data have been collected as a component of larger studies, the species has rarely been the primary focus because other avian species, such as the Golden-winged Warbler, have experienced more precipitous population declines. Morris et al. (2013) evaluated avian response to silviculture in the Ozark Mountains of Missouri and found that Blue-winged Warbler density significantly increased following even-and uneven-aged harvests. At sites managed for Cerulean Warblers using even-aged harvest with varying levels of residual basal area in the Central Appalachians, Blue-winged Warbler density peaked at low residual basal area (3-5 m²/ha; Sheehan et al. 2014). With Bluewinged Warbler populations continuing to decline, evaluating the effects of various sources of targeted disturbance on the species needs to be considered.
Our objective was to evaluate the response of Blue-winged Warblers and other early successional avian species to NRCS conservation projects in West Virginia. Additionally, we aimed to identify key habitat characteristics associated with Blue-winged Warbler site occupancy. We addressed these objectives by surveying the avian community, vegetation, and habitat characteristics on pre-treatment and post-treatment sites. Because Blue-winged and Golden-winged warblers are often documented nesting in the same areas (Confer et al. 2020), we reviewed previous studies of Golden-winged Warblers to select spatial and habitat characteristics that we hypothesized would influence Bluewinged Warbler site occupancy. We predicted that Blue-winged Warbler site occupancy would be higher at sites managed for Golden-winged and Cerulean warblers. Additionally, we predicted that species richness for shrubland and grassland guilds would be greater at post-treatment sites enrolled in the Goldenwinged Warbler Project, as these species are similarly disturbancedependent. Results from our study may inform whether habitat management for Golden-winged Warblers benefits Blue-winged Warblers and the early successional species guilds or if targeted efforts may be needed in the future. Additionally, research into Blue-winged Warbler site occupancy may help inform their expanding Appalachian distribution, which could have useful implications for Golden-winged Warbler conservation.

Study Area
We sampled 20 sites enrolled in the Golden-winged Warbler Project, 19 sites enrolled in the Cerulean Warbler project, and 2 sites enrolled in both projects on private lands in West Virginia during 2019-2020 ( Fig. 1). Among the sites (n = 41) surveyed in these two projects, we surveyed post-treatment sites (n = 29) that had been certified by NRCS as completed and pre-treatment sites (n = 12) where habitat management had been planned by NRCS but not yet initiated. We considered treatment the species-specific conservation practices administered by NRCS (Appendix 1). Treatments followed recommendations in the habitat management guidelines for Golden-winged and Cerulean warblers (Bakermans et al. 2011, Roth et al. 2012, Golden-winged Warbler Working Group 2013, Wood et al. 2013). There were 12 conservation treatments which included a range of manipulations designed to improve habitat for the target species. Treatments were not mutually exclusive (sites could have received more than one treatment) but treatments were designed specifically for the Appalachians. Although the treatments were different and variable, we hypothesized that they would positively affect site occupancy of Blue-winged Warblers and other disturbancedependent species because they all disturb the forest vegetation structure in a way that should elicit a positive site occupancy response. We did not collect pre-and post-treatment data on the same sites because contract implementation occurred over a 1-to 6-year period. We surveyed sites that were 1-8 years posttreatment to account for multi-year lags in avian response (Sheehan et al. 2014, Bakermans et al. 2015, Aldinger 2018. We placed sample points in three defined site categories: 1) pretreatment, which refers to sites with planned NRCS management but have not been initiated, 2) post-treatment sites that were certified by NRCS and treatments that had been completed 1-8 years prior to this research, 3) reference sites, which are areas outside of treated sites but within the property boundary and should be representative of pre-treatment conditions. For data analyses, we created a categorical variable by combining points on pre-treatment and reference sites into a single "reference" category and compared those to treatment site points. We also evaluated years post-treatment as a continuous variable, with values ranging from 0 (untreated) to 8 years. We verified the number of years since treatment through communication with landowners and NRCS personnel. Total area of treated areas was 1.2-36.3 ha (mean ± standard error [SE] = 9.3 ± 1.5 ha). Treatment areas on 11 of 29 post-treatment sites included multiple noncontiguous patches within the property boundary. These noncontiguous patches were 1.0-36.3 ha (mean ± SE = 4.9 ± 0.96 ha) and varied in shape, elevation, slope position, and aspect.

Fig. 1. Site locations enrolled in Natural Resources
Conservation Service practices for Cerulean (Setophaga cerulea) or Golden-winged (Vermivora chrysoptera) Warblers in West Virginia sampled during 2019-2020. Site locations were randomly offset within a 5-km radius of the actual location to protect personally identifiable information. CERW = Cerulean Warbler and GWWA = Golden-winged Warbler.
All sites were located within the Central Appalachians portion of the Appalachian Mountains Bird Conservation Region (NABCI Committee 2000) and occurred within the Central Appalachians, Western Allegheny Plateau, and Ridge and Valley level III ecoregions of West Virginia (Omernik and Griffith 2014; Fig. 1). The Central Appalachians ecoregion is a high elevation and rugged plateau with extensive forest cover dominated by northern hardwood and mixed-mesophytic forests (Woods et al. 1999). The Western Allegheny Plateau is characterized by sharp ridge tops and narrow valleys and is dominated by Appalachian oak (Quercus spp.) forests. The Ridge and Valley is characterized by elongated alternating forested ridges and agricultural valleys and is dominated by dry oak-pine (Pinus spp.) forests. These ecoregions contain variable soil conditions, geology, climate, elevation, and precipitation (Woods et al. 1999 , and yellow poplar (Liriodendron tulipifera). Shrublands in the study area primarily result from localized anthropogenic disturbances such as livestock grazing, timber harvesting, mechanical cutting for oil and gas wells (e.g., mowing or brushhogging), and surface mine reclamation. Shrublands on our study sites were primarily the result of active/abandoned livestock grazing or timber harvesting, which commonly occur throughout the study area. Common shrub species in the study area include autumn olive (Elaeagnus umbellate), blackberry and raspberry (Rubus spp.), hawthorn (Crataegus spp.), and multiflora rose (Rosa multiflora) (Aldinger 2018). Scattered areas of red spruce forest (Picea rubens) occur at the highest elevations.

Sample point placement
Using ArcGIS, we systematically placed sampling points 200-943 m apart within treatment and reference areas to avoid double counting individual birds (Ralph et al. 1995). For small and/or irregular shaped treatment areas, we placed points in the approximate center of the treatment area for maximum coverage. We did not use random point placement; rather post-treatment points were placed >4 m from the edge of treated boundaries in treatment areas. Though this is an unconventional way to place points, we wanted our sample inference frame (100 m) to include treatment effects which could extend beyond the physical treatments (Paton 1994), and we were limited by the size of treatments. All point centers were >200 m apart, thus reference points were >100 m from treatment areas. This was to ensure no opportunity for overlap of treatment effects on reference points. For sites enrolled in the Golden-winged Warbler Project, we placed reference points in existing early to mid-successional vegetation communities on the same property that contained elements of Vermivora habitat (i.e., shrubland, young forest, and forest edges). For sites enrolled in the Cerulean Warbler project, we placed reference points in adjacent forested areas on the same property with similar slope, aspect, elevation, and forest community type as treated points, whenever possible. Consequently, pre-treatment and reference points had similar habitat conditions and we assumed conditions at reference points were representative of the pre-treatment condition at treated points.

Point count survey protocol
We surveyed the avian community from 7 May to 1 July in 2019 and 2020 to coincide with peak songbird breeding season in the Central Appalachians (Canterbury et al. 1993, Bakermans et al. 2015, Aldinger 2018. We surveyed 202 points in 2019, 105 points in 2020, and 34 points in both 2019 and 2020 (Table 1). We surveyed all points twice per year, with at least one week between surveys. We began surveys 15 minutes before official sunrise and ended 4 hours after official sunrise on days without consistent rain or sustained winds >19 km/hour (i.e., >3 on the Beaufort scale). Cloud cover, wind, and start time were recorded for each survey to incorporate into detection models. A two-person field crew conducted the avian sampling each year (three total observers). All observers were experienced in bird identification, and all were trained in local bird identification, distance estimation with range finders and sampling protocols before surveys began. We recorded all first detections seen or heard within a 10-minute period. All individuals detected were recorded once and monitored throughout the remainder of the survey to avoid double counting. We recorded detections by distance band of the bird from the observer (0-25 m, 25-50 m, 50-75 m, 75-100 m); we did not record or include any detections from beyond 100 m. We recorded detection type (song, call, visual, or fly-over) and sex for each detection, if possible.

GIS-derived covariates
We first used 2018 National Agriculture Imagery Program (NAIP) imagery (1-m cell size; data collected during the growing season) to manually digitize and identify land cover types deemed important for Vermivora breeding habitat in previous studies (Oliver 2021; Appendix 2). Additionally, we did ground truthing during field work to verify cover type accuracy from GIS imagery.  Gill et al. 2020). When delineating shrubland and young forest, the amount of woody cover was determined visually using 2018 NAIP imagery (Bakermans et al. 2015, Aldinger et al. 2018. We also included the sum of young forest (of trees ≤8 cm dbh) and shrubland (≥30% shrub cover) within 100-m of each point in analysis. We used a 1:24,000 digital elevation model (US Geological Survey 1999) in ArcGIS version 10.7 (ESRI 2011) to calculate elevation within a 100-m radius centered on each sampling point by obtaining the majority from the raster layer in ArcGIS, using the "zonal statistic as table" tool.
We used 5-m resolution land cover data generated from 2016 NAIP imagery (Maxwell et al. 2019) to calculate percent forest cover (deciduous, mixed, and conifer combined) and edge density (km/km²; Roth et al. 2012) within 1 km of each sampling point location (Oliver 2021). We included percent forest cover and edge density because Vermivora breeding habitat is associated with forested landscapes, and Vermivora populations will not persist in highly fragmented landscapes (Roth et al. 2012). We created 1km buffers around each sampling point and evaluated percent forest by dividing the total forested area (ha) within each buffer by the area (ha) of each buffer and multiplying by 100. We chose a 1-km buffer following the guidelines recommended for landscape management for Vermivora, and because others have documented significant landscape-level effects (Bakermans et al. 2015, Aldinger 2018. Within each 1-km buffer, we also estimated edge density (km/km²) by dividing the length of forest edge within each buffer by the area of buffer.

Vegetation covariates
We measured vegetation at the point level to characterize vegetation communities at each site for structural metrics deemed important for Vermivora breeding (Confer and Knapp 1981, Askins et al. 2007, Gill et al. 2020. We conducted vegetation sampling once each year after the second point count survey to ensure vegetation growth was maximized within the survey period. We measured percent cover of grass, forbs, Solidago, Rubus, shrubs, saplings (1-8 cm dbh and ≥1 m tall), and canopy trees (>8 cm dbh) within each plot based on ocular tube "hits" at 10 points spaced 10 m apart along each of three 100-m transects, radiating from plot center at 0°, 120°, and 240° (adapted from James andShugart 1970 andMcNeil 2019;Oliver 2021). We recorded the occurrence of each vegetation class when looking straight up and down through the vertical ocular tube. Percent cover of each vegetation class was calculated as the number of "hits" for a class divided by the total number of ocular tube sampling points (n = 30) and multiplied by 100. We recorded maximum vegetation height to the nearest meter by taking a reading with a handheld range finder, looking straight up to the approximate tallest tree or shrub. We calculated the coefficient of variation as the standard deviation divided by the mean for the 30 maximum vegetation height values. We estimated coefficient of variation for maximum vegetation height for each sampling point because Vermivora nest sites have more variable vegetation height than random locations within the same territory (Aldinger 2018).
We measured canopy tree basal area (Golden-winged Warbler Working Group 2013, Wood et al. 2013) at four prism plots per sampling point location placed at the point center and 70 m from point center along the 0°, 120°, and 240° transects. Variable radius prism plots were completed using a wedge prism (equivalent to a 2.2296-BAF metric prism) to tally live trees and snags ≥8 cm dbh at each prism plot (Nareff et al. 2019). The observer tallied trees that were "in" the plot by sighting through the prism held over plot center and rotating around the point. For each "in" live tree, we recorded tree species or group (e.g., hickory group, red oak group) and dbh measured to the nearest centimeter (cm) using a dbh tape. For each "in" snag ≥2 m in height, we recorded dbh. All borderline (i.e., trees that were unclear whether they were "in" or "out" of the prism plot) live trees and snags were counted, but every other borderline live tree was removed for basal area calculations (Bell and Alexander 1957). We calculated basal area at each plot by multiplying the number of "in" live trees by a basal area factor of 2.296. Mean basal area across all four prism plots (100-m basal area) reflects basal area across the 100-m radius plot and, therefore, may include tree tallies from within and outside the treatment area.

Species richness and Habitat guilds
We evaluated species richness of shrubland and grassland bird species on sites managed for Golden-winged and Cerulean warblers. We classified songbird species a priori into these habitat guilds based on breeding biology and previous guild studies from the region (McDermott and Wood 2009; Appendix 3). The shrubland guild consisted of 26 species that breed in early successional vegetation communities (e.g., regenerating forest following recent disturbance, old fields). The grassland guild consisted of 4 species that breed in grasslands (e.g., hayfields, pastures). We included all auditory and visual detections, except flyovers of adult males within a 100-m radius of each point count in species richness analysis. We defined species richness at each point as the total number of bird species within each habitat guild detected across the two visits per year. We chose these guilds because they include disturbance-dependent species, which continue to experience population declines (Rosenberg et al. 2016) and are groups of birds most likely to be positively affected by the management on sites.

Statistical Analysis
We estimated Blue-winged Warbler site occupancy using singleseason occupancy models. Single-season occupancy models use spatially and temporally repeated surveys to estimate the proportion of sites occupied by a species of interest as a product of ecological processes and imperfect detection by linking two sub-models (MacKenzie et al. 2002). For the 34 points on sites surveyed in 2019 and 2020, we included only the 2020 data in analyses because including data from both years was not compatible with our site occupancy analyses and models would not converge using only 2019 data due to low detections. We included all detection types except flyovers (e.g., singing, visual, and call) of adult males within 100 m of sampling point in analysis to relate site occupancy more closely to vegetation and topographic characteristics measured within 100-m radius of each sampling point. Point count data were formatted as detection (1) and non-detection (0) per visit for estimation of detection probability across each survey. We used the occu function in package unmarked version 0.11-0 (Fiske and Chandler 2011) in program R version 3.5.2 (R Development Core Team 2019) for all hierarchical modeling. We considered models with the lowest Akaike's Information Criterion score adjusted for small sample sizes (ΔAIC c ; Burnham and Anderson 2002) to be the best supported model given the data and any models with ΔAIC c <2.0 were considered plausible (Burnham and Anderson 2002).
We first included each detection covariate (observer, ordinal date, time-since-sunrise, sky, wind, and canopy; Table 2 and Appendix 2) in its own detection model to identify which detection covariates were best supported using AIC c (Burnham and Anderson 2002). We then checked confidence intervals of beta estimates, and if confidence intervals did not overlap zero, we created an additive model with those covariates (Arnold 2010). We tested three occupancy model sets. We chose to use a stepwise modeling approach because we did not want treatment effects to be confounded by landscape or vegetation variables, but we wanted to account for those other potential relationships. Table 3. Single-season occupancy model selection procedure results for Blue-winged Warbler (Vermivora cyanoptera) occupancy (ψ) and detection probability (p) using 341 pointcount locations on 41 private properties enrolled in Natural Resources Conservation Service practices for Cerulean (Setophaga cerulea) or Golden-winged Warblers (Vermivora chrysoptera) in West Virginia sampled during 2019-2020. K is the number of parameters in each model, AICc is the Akaike's Information Criterion value for small sample sizes, ΔAICc is the difference between each model's AICc value and the lowest AICc value in the candidate set, and wi is the Akaike weight of each model in relation to each model set. Codes for covariates are defined in Appendix 2.  Proportions of cover types or vegetation structure can affect early successional bird site occupancy differently and be scale dependent (Thogmartin 2010, West et al. 2016. In the first occupancy model set, we evaluated Blue-winged Warbler site occupancy at reference and treated locations to evaluate the conservation project's effectiveness (Table 3 and Appendix 2). Our first model set had seven models using contract treatment status (treated vs. reference) and years-post-treatment as the predictor variables and included ecoregion and elevation to account for inherent differences in site occupancy among regions. According to data present in the inclusive thesis (Oliver 2021), we did not include interactive effects between contract treatment status and vegetation covariates because very few vegetation metrics in our Table 4. Beta (β) estimates, standard errors (SE), and 95% lower and upper confidence intervals (CI) on the logit-scale from models (Table 3)  study differed among contract treatment categories, indicating that management implementation did not meaningfully influence vegetation structure and composition on treated sites. In the second model set, we developed six models to examine the influence of GIS-derived variables (Table 3 and Appendix 2). In the third model set, we developed 10 models to examine the influence of vegetation covariates (Table 3 and Appendix 2). We evaluated Goodness of Fit for top models using the MacKenzie and Bailey (2004) fit statistics for site occupancy. If models did not fit well (p < 0.05), we considered alternative distributions and analyses. Site occupancy predictions were derived using the "predict" function in program R and differences were assessed based on 95% confidence intervals of beta estimates.
We used a one-way analysis of variance (ANOVA) including site name as a random effect to compare observed species richness of shrubland and grassland avian guilds reference and treated locations categories using the "lme" and "summarySE" functions in program R. For the 34 points on sites surveyed in 2019 and 2020, we included only the 2020 data in analyses. For significant ANOVAs, we then identified which contract treatment status categories differed using a post hoc Tukey-Kramer test. We set statistical significance at p <0.05.

RESULTS
We detected 49 Blue-winged Warblers (on 18 points at four Golden-winged Warbler sites and 16 points at nine Cerulean Warbler sites) across 341 points surveyed during the study. Naïve occupancy, which is unadjusted using detection probability, was highest on reference points at sites enrolled in the Golden-winged Warbler Project (Table 1).
Observer (Obs1 β = 3.59 ± 2.57; Obs2 β = 4.10 ± 1.00; Obs3 β = 2.95 ± 0.87) and date (β = -0.044 ± 0.17) were the only supported detection covariates and were used as the detection covariates in models for all subsequent analyses (Table 2). All models passed goodness-of-fit tests (p > 0.05; c-hat ≤ 1.40). For model set one, model selection results indicated that a model for ecoregion and contract treatment status and a model with only ecoregion were supported (Table 3). Treated points had lower site occupancy than reference points by 34-44% depending on ecoregion (Table 4 and Fig. 2). Site occupancy in the Central Appalachians ecoregion was 300% greater (0.25 +/-0.07) than in the Ridge and Valley (0.06 +/-0.03; Fig. 3). From model set 2, models examining GIS variables containing 100-m shrubland and 100-m young forest/ shrubland in the immediate landscape were best supported (Table  3). Site occupancy increased by 880% as 100-m shrubland increased to 100% and by 922% as 100-m young forest/shrubland increased to 100% (Table 4 and Figs. 4 and 5). For model set 3, a model containing 100-m basal area was best supported (Table 3).
Site occupancy decreased by 93% as 100-m basal area increased from 0 to 30 m²/ha (Table 4 and Fig. 6). When comparing site means of covariates between reference and treated points, 100-m shrubland were the same and 100-m young forest was 6% greater on treated points (Appendix 4).

Fig. 2. Predicted occupancy and 95% confidence intervals of Blue-winged Warblers at 341 point-count surveys on 41 private properties enrolled in the Natural Resources Conservation
Service practices for Cerulean (Setophaga cerulea) or Goldenwinged Warblers (Vermivora cyanoptera) in West Virginia sampled during 2019-2020. Presented is occupancy related to treatment status and observer + date for detection probability from top models in Model Set 1. Differences are indicated by alphabetic notation at α = 0.05.
We detected four species in the grassland guild and 26 in the shrubland guild (Appendix 3), across 147 point-count surveys on sites managed or planned to be managed for Golden-winged and Cerulean warblers. There were no differences among the status categories for shrubland (F = 1.56, p = 0.19) or grassland (F = 1.39, p = 0.24) richness (Figs. 7 and 8). Presented is occupancy related to ecoregion and observer + date for detection probability from top models in Model Set 1. Differences are indicated by alphabetic notation at α = 0.05.  Presented is occupancy related to percentage of 100-m young forest/shrubland and observer + date for detection probability from top models in Model Set 2.

DISCUSSION
Private lands conservation projects are being implemented to benefit many imperiled species of conservation concern, but the effects of implementation on these co-occurring species are seldom evaluated. Evaluating the response of co-occurring species is critical to assess the potential of these conservation projects to provide "umbrella" benefits to other wildlife species.
We did not observe positive responses for Blue-winged Warbler site occupancy to Golden-winged and Cerulean warbler management, nor species richness in the grassland or shrubland avian guilds to Golden-winged Warbler management, suggesting the disturbance-dependent avian community did not indirectly benefit from habitat management. We suspect that the lack of avian response was influenced by private lands variability, site selection, small treatment areas, and challenges associated with conservation project implementation resulting in insufficient vegetation structure and composition. For disturbancedependent species specifically, there was low shrubland cover on treated sites enrolled in the Golden-winged Warbler Project, likely resulting from high basal area retained during implementation.
While Blue-winged Warblers were not positively correlated with implementation of habitat management for Golden-winged and Cerulean warblers, we did identify important habitat characteristics that have implications for habitat management. Blue-winged Warbler site occupancy was positively associated with shrubland. We defined shrubland cover as a digitized area with ≥30% shrub cover generally dominated by shrub species with scattered canopy trees and herbaceous understory. The importance of shrubland cover for the closely related Goldenwinged Warbler has been repeatedly demonstrated (Hanowski 2002, Bulluck and Buehler 2008, Roth et al. 2012, Aldinger and Wood 2014, Aldinger 2018. Shrubs are important to Vermivora because nest sites are often located at ecotones between forest edges, shrubs, and herbaceous vegetation, and foraging usually occurs on arthropods from the peripheral branches of trees and shrubs (Gill et al. 2020). While shrubs are a preferred Blue-winged Warbler habitat component (Confer and Knapp 1981, Confer et al. 2003, Gill et al. 2020, the importance and amount of shrubland cover had not been quantified for Blue-winged Warblers. Although habitat management on our sites was not targeted for Blue-winged Warblers and our study was limited by sample size, spatial extent, and treatment implementation, we found that 50-100% shrubland cover at the 100-m radius scale produced highest site occupancy for Blue-winged Warblers. Shrubland conditions can be achieved by setting back succession through the implementation of disturbance (e.g., timber management, prescribed burning, grazing, and mechanical clearing) or advancing succession by reducing disturbance to allow for shrub establishment in existing early successional communities.
We predicted that Blue-winged Warblers would benefit from the guidelines for Golden-winged and Cerulean warblers, however, we believe that the Golden-winged Warbler guidelines will be more effective than Cerulean Warbler guidelines when managing for Blue-winged Warblers. Blue-winged Warbler occupancy on our sites was negatively associated with basal area. We detected only 4 Blue-winged Warblers at 4 sites treated for Cerulean Warblers, and all were detected on points with basal area from 2.3-18.4 m²/ha and contained shrubland and/or young forest  cover. Thus, similar to Sheehan et al. (2014), management for Cerulean Warblers (recommended residual basal area of 9.2-20.7 m²/ha) is unlikely to benefit Blue-winged Warblers. The minimum 100-m basal area for points where Blue-winged Warblers were detected was 2.0 m²/ha (Oliver 2021), which is near the minimum basal area recommendation for Golden-winged Warblers of 1.9 m²/ha. Site occupancy decreased by 25% above 4.6 m²/ha (Fig.  2E). Therefore, we recommend residual basal area of 1.9-4.6 m²/ ha when managing for Blue-winged Warblers in West Virginia. This recommendation is consistent with previous research indicating that Blue-winged Warblers use early and midsuccessional communities (Confer and Knapp 1981, Confer et al. 2003, Gill et al. 2020). Our basal area finding suggests that residual canopy trees are not necessary to incorporate in habitat management for Blue-winged Warblers, unlike recommendations in the Golden-winged Warbler habitat management guidelines (Bakermans et al. 2011, Roth et al. 2012, Golden-winged Warbler Working Group 2013. However, managers should use caution when drastically reducing basal area in large areas for Bluewinged Warblers, as territories normally include a forest edge (Patton et al. 2010). The target basal area range can be achieved by implementing silviculture treatments (e.g., clearcutting, seed tree harvests, and overstory removal with residuals) in existing forest or shrubland management practices (e.g., mechanical clearing) to create or enhance existing early successional vegetation communities.
Ecoregion strongly influenced occupancy on our sites, as occupancy predictions were highest in the Central Appalachians and lowest in the Ridge and Valley. Although 36% of the 341point count locations were in the Ridge and Valley ecoregion, Blue-winged Warblers (n = 4) were only detected on one site. Bluewinged Warbler density is low in West Virginia's Eastern Panhandle (Bailey and Rucker 2021), as this area is near the eastern edge of the species' breeding distribution (Gill et al. 2020). The highest elevations in West Virginia occur in the Central Appalachians ecoregion. Historically, Blue-winged Warblers were restricted to lower elevations (<600 m, Buckelew and Hall 1994), but the species' breeding distribution has expanded into higher elevations (maximum elevation detected = 856 m; Bailey and Rucker 2021). We detected Blue-winged Warblers across most of the elevation range surveyed (244-917 m; Appendix 4). Our findings add to the growing body of research indicating that their breeding distribution is continuing to expand to higher elevations. Populations will require future monitoring to document continued upslope expansion to possibly identify priority areas for habitat management should population declines accelerate.
Our study is one of only a handful focusing on Blue-winged Warbler breeding habitat characteristics especially investigating response to habitat management for Golden-winged and Cerulean warblers. Habitat characteristics important to Bluewinged Warblers are consistent with those previously identified for Golden-winged Warblers, which is unsurprising as both species regularly nest in the same areas (Confer et al. 2020). However, our basal area recommendation of 1.9-4.6 m²/ha encompasses a wider range than suggested for Golden-winged Warblers (~1.9-3.7 m²/ha; Bakermans et al. 2011, Roth et al. 2012, Golden-winged Warbler Working Group 2013, suggesting that Blue-winged Warblers are found in a broader range of vegetation communities Knapp 1981, Confer et al. 2003). Based on important breeding habitat characteristics that we identified, we believe that habitat management has potential as a conservation action to benefit populations in the Central Appalachians, Western Allegheny Plateau, and Ridge and Valley of West Virginia, given the species' reliance on disturbance. Based on our results, we recommend 50-100% shrubland cover at the 100-m radius scale and residual basal area of 1.9-4.6 m²/ha when managing habitat for Blue-winged Warblers in West Virginia. Though these results are limited by the scope of our study area and sample design, we are confident in their applicability as a guide and baseline for Blue-winged Warbler management to improve site occupancy in the region. Future studies with improved management implementation and larger sample sizes across the species' breeding range will be required to further inform the potential of habitat management as a range-wide conservation action for Blue-winged Warblers.