Home range-and colony-scale habitat models for Pinyon Jays in piñon-juniper woodlands of New Mexico , USA

Despite the Pinyon Jay’s (Gymnorhinus cyanocephalus) importance as a piñon tree seed disperser and its long-term population declines, piñon-juniper habitat use by Pinyon Jays has been little studied. We created home rangeand colony-scale habitat suitability models for Pinyon Jays in piñon-juniper woodlands at Kirtland Air Force Base and White Sands Missile Range, New Mexico, USA. Home range-scale models included mainly Juniper Woodland and Savanna, Piñon-Juniper, and Piñon Woodland vegetation types. Although home range-scale models at the two study sites contained similar vegetation types, we found differences between colonyscale models in habitat use and within-site diversity. Because of topographic and vegetation differences, colony-scale modeling incorporated fewer vegetation types and topographic classes for the White Sands Missile Range Study site than at Kirtland Air Force Base. Predictive colony-scale models overlapped validation colonies by 75.9% and 100% at the two study sites, indicating good model predictability. The large sizes of Pinyon Jay home ranges make year-round, home range-scale habitat management challenging. Habitat conservation focused at more limited spatial, i.e., colony, and temporal scales, revealed by multiscale modeling, could be more feasible and potentially more effective than management at the larger home range scale. Modèles d'habitat élaborés aux échelles du domaine vital et de la colonie chez le Geai des pinèdes dans les boisés de pins et de genévriers au Nouveau Mexique, États-Unis RÉSUMÉ. Malgré l'importance du Geai des pinèdes (Gymnorhinus cyanocephalus) en tant que vecteur de dispersion de graines du pin pignon et la baisse de longue date de ses effectifs, l'utilisation de l'habitat de pins-genévriers par ce geai a été peu étudiée. Nous avons élaboré des modèles de la qualité de l'habitat aux échelles du domaine vital et de la colonie chez le Geai des pinèdes dans les boisés composés de pins et de genévriers à la Kirtland Air Force Base et au White Sands Missile Range au Nouveau-Mexique, ÉtatsUnis. Les modèles à l'échelle du domaine vital incorporaient essentiellement les types de végétation allant des boisés et des savanes dominés par les genévriers, aux boisés mixtes de pins-genévriers et aux pinèdes. Même si les modèles à l'échelle du domaine vital comportaient des types de végétation semblables aux deux sites d'étude, nous avons trouvé des différences entre les modèles à l'échelle de la colonie pour ce qui est de l'utilisation de l'habitat et la diversité à l'intérieur d'un même site. En raison de différences sur le plan de la topographie et de la végétation, la modélisation à l'échelle de la colonie incorporait moins de types de végétation et de classes topographiques pour le site de White Sands comparativement au site de Kirtland. Les prédictions issues des modèles à l'échelle des colonies correspondaient à 75,9 et à 100 % à la localisation des colonies de validation aux deux sites d'étude, ce qui indique une bonne capacité prédictive des modèles. Étant donné les grands domaines vitaux chez ce geai, l'aménagement de l'habitat à cette échelle et à longueur d'année représente un défi de taille. La modélisation multi-échelles révèle que les mesures de conservation de l'habitat fondées sur des échelles temporelle et spatiale plus limitées, c.-à-d. à l'échelle de la colonie dans le cas de cette dernière, seraient plus réalistes et sans doute plus efficaces que l'aménagement à l'échelle plus grande du domaine vital.


INTRODUCTION
Pinyon Jays (Gymnorhinus cyanocephalus) and piñon pines have a mutualism whereby the trees provide mast crops of nutritional seeds that affect the jays' population viability, and the jays cache the seeds in the ground, serving as the tree's main long-distance seed disperser (Marzluff andBalda 1992, Lanner 1996).Pinyon Jays also depend on piñon and juniper trees for nest sites (Johnson et al. 2014(Johnson et al. , 2015)).Despite recent climate impacts to piñon-juniper vegetation and a long-term decline of Pinyon Jay populations (Sauer et al. 2014), there is "no known detailed quantification of habitat anywhere within its range" (Balda 2002:4).To manage and conserve Pinyon Jays and their piñon-juniper habitats in the face of a changing climate and other impacts, it is necessary to understand the habitat requirements of the piñon pine's most important long-distance seed disperser.
Although fire has been found to play a role in structuring some piñon-juniper woodlands in the southwestern USA, birds, insect pests, and drought are also critical in the establishment and demise of these woodlands (Romme et al. 2009).Pinyon Jays serve as short-and long-distance seed dispersers for piñon pines.Adapted for carrying and caching millions of seeds in a few weeks, a Pinyon Jay flock is capable of replanting a woodland decimated by fire, drought, or insect pests (Ligon 1978).The Clark's Nutcracker (Nucifraga columbiana), the only other long-distance seed disperser of piñon seeds, is a less important seed disperser than the Pinyon Jay because the nutcracker forages in smaller flocks, forages in a wide variety of pine forests and woodlands, and often occurs at higher elevations than piñon pine woodlands (Tomback 1998).
Pinyon Jays are highly social, nesting colonially at traditional colony sites, and sometimes breeding cooperatively, with helpers at the nest (Marzluff andBalda 1992, Balda 2002).They are yearround residents of the southwestern USA, their distribution completely overlapping that of piñon trees (Marzluff and Balda 1992).When not nesting, they commonly travel over large landscapes in foraging flocks of up to several hundred individuals (Balda 2002).However, the sizes of nonbreeding flocks in recent anecdotal reports are much smaller than those reported historically (Bailey 1928, Linsdale 1936), as Pinyon Jay populations have been declining for over 40 years (Sauer et al. 2014).
Piñon trees produce mast crops of piñon seeds, which provide critical nutrition that enhances Pinyon Jay reproductive success and survival (Ligon 1978, Marzluff andBalda 1992).The Pinyon Jay's extraordinary physiological and behavioral dependence on piñon seeds is exemplified by reversal of gonadal regression in winter and breeding behavior in 11 months of the year, stimulated by abundant piñon seeds (Ligon 1978).Pinyon Jays also nest in larger than average piñon and juniper trees and in areas of high canopy cover (Johnson et al. 2014(Johnson et al. , 2015)).
Because of Pinyon Jay dependence on piñon pines for food and nest sites (Ligon 1978, Lanner 1996), climate impacts to piñon mast production, tree health, canopy cover, and distribution will likely negatively impact the jay's populations.Decline of the piñon's primary long-distance seed disperser will limit the potential of the trees to re-establish in areas of high mortality, colonize higher elevations, or shift distributions northward in response to climate impacts, thus limiting climate resilience of piñon habitats.Pinyon Jay populations in turn will be further impacted by habitat loss.
Pinyon Jays must make habitat choices on at least three scales (Fig. 1): 1.The home range of a Pinyon Jay flock typically includes thousands of hectares and several land cover or ecosystem types, both characteristics of the landscape scale (Ferrari and Ferrarini 2008).Hence, because of their large sizes and diverse land cover types, Pinyon Jay home range habitat models are at the landscape-scale (hereafter "home range scale" for simplicity).

3.
Each Pinyon Jay pair must then choose a suitable tree and nest site within the nesting colony, on a scale of approximately one hectare.Because home range-and colony-scale models are both commonly based on geospatial data (Ferrari and Ferrarini 2008), we address these two scales here (Fig. 1).Habitat use at the nest scale, beyond the scope of this paper, is addressed separately (K.Johnson, G. Sadoti, and J. Smith, unpublished manuscript).The goal of this study was to create predictive geospatial (geographic information system, GIS) models of Pinyon Jay habitat use at the home range and nesting colony scales in piñonjuniper woodlands, to inform habitat conservation and Pinyon Jay management.Specific objectives were the following: (1) to establish differential habitat use at the scale of the home range (based on radio telemetry and ground occurrence data); (2) to create predictive models at the scale of the nesting colony (delineated based on nest locations within colony sites) to identify a limited subset of the available piñon-juniper habitat, with characteristic vegetation type and topography suitable for nesting; and (3) to understand scale implications in applying the geospatial models to Pinyon Jay conservation and management.

Study areas
The two study areas were chosen because each was home to a Pinyon Jay flock of up to 50 individuals that maintained a home range and nested therein at traditional nesting colonies.Nesting colonies differed in topography and vegetation type and thus presented the opportunity to model a range of nesting colony habitats.

White Sands Missile Range home range-scale habitat use
We visited the study site at WSMR on average 1.1 times per week (range 0 -3) from 24 March through 29 October 2009.On 10 April 2009, we set up a solar-powered, battery-operated automatic feeder (Sweeney Enterprises, Boerne, TX) within the area of Colony WA.
The feeder delivered about two cups of P. edulis seed twice daily.
On five days from June-August 2009, we trapped Pinyon Jays in a modified Australian crow trap set near the feeder and baited with P. edulis seed.We processed captured birds as at KAFB.
We attached radio transmitters to four of the birds captured at Colony WA.We tracked birds with a TRX-1000 receiver, either following the birds until we could see them or triangulating on their location as at KAFB.We attempted to locate the radio-tagged birds weekly from the time we banded them until the last transmitter batteries failed or birds left the area.We mapped locations as at KAFB.We collected bird location data for the modeling study in 2009, but we had additional Pinyon Jay location data from ground surveys in 2005, 2008, 2009, and 2011 at WSMR, which were incorporated into models.

Both study sites
Telemetry points indicated the location of an individual, but because Pinyon Jays are highly social, individuals were typically with other jays; hence, individual locations were indicative of a flock's home range.To delineate home ranges, using Arc GIS 10.2 (Esri 2013) we created minimum bounding geometry polygons around the set of detection points, excluding repetitive points at nests and feeders.
We adjusted the use of habitat type (land cover class) according to availability within the year-round Pinyon Jay home ranges by computing a concentration of use index (Neatherlin and Marzluff 2004), the probability of use for a land cover type in the home range divided by the occurrence of that land cover in the home range.Land cover classes composing < 5% of the breeding, nonbreeding, or year-round home range tended to have inflated concentration of use values relative to the most-used habitat types (one detection in a very small area results in an inflated use value).We therefore report concentration of use only for those land cover types composing > 5% of the home range polygons.

Home range-scale habitat modeling
We created home range-scale models using Pinyon Jay observation data from KAFB 2009-2011and WSMR 2005, 2008, 2009, and 2011.We first created a land cover map over the KAFB study area.
A vegetation classification, map, and associated plot data (E.Muldavin, G. Harper, P. Neville, Y. Chauvin, and P. Mehlhop, unpublished data) were available for the WSMR study area.

KAFB home range-scale habitat modeling
The Mid-Region Council of Governments provided six-inchresolution, natural color, digital ortho-imagery of Bernalillo County (MRGCOC 2008) flown in March and April 2008 covering all habitats at KAFB.Using this imagery, we delineated 72 polygons that contained Juniper Woodland and Savanna, Piñon-Juniper Woodland, and Piñon Pine Woodland overlapping and surrounding areas where Pinyon Jays had been observed across the study area.We visited 69 of the 72 polygons, where we took photos and made preliminary assignments to vegetation type.Mean polygon area was 8.6 ha (range: 0.12 -67 ha; total polygon area = 597 ha, 6% of final mapped area).We used these data to define map units for a classification of vegetation in the study area.
At scales of 1:3000 -6000, we applied aerial photo-interpretive techniques of visual landscape elements including color, texture, and size, with overlays of elevation contours and field data descriptions, to create the land cover map.Finally, we delineated 325 polygons for potential field checking in areas of known Pinyon Jay activity in areas with access unrestricted by KAFB.
We visited 89 of these (mean area: 3.8 ha, range: 0.03 -21 ha) that were reasonably accessible and concentrated in the habitat types used by the jays.We collected the following data: date, GPS coordinates, aspect of the described slope, percent cover class of the dominant species, canopy cover, and relative cover of piñon and juniper.These data were used to refine the habitat classification.
Dick-Peddie (1993) defined woodlands as having trees with nonoverlapping canopies and distinguished juniper savanna as scattered stands with densities less than 130 trees/acre.Using these definitions as a guideline, we classified those areas where piñon accounts for 25 -50% of the cover as Piñon-Juniper Woodland.
We classified Piñon Woodland as areas where piñon makes up > 50% of the cover (based on our field data, extrapolated to similar areas on the imagery).Juniper-dominated classes contain < 25% piñon.We considered the vegetation types overlapping Pinyon Jay occurrence and home range polygons to be Pinyon Jay habitat but also mapped all nonwoodland land cover types within the study area.

WSMR home range-scale habitat modeling
To model home range-scale habitat use by Pinyon Jays at WSMR, we referenced the existing vegetation map and associated plot data, in combination with 1-m, natural-color, aerial photography acquired in July 2009 (NAIP 2009); existing data on tree species and size; and Pinyon Jay occurrence datasets (Johnson and Smith 2006, 2008, Johnson et al. 2014).We assigned all land cover types within the study area.We used our detection points and home range polygons to assign the vegetation types containing Pinyon Jay occurrences to Pinyon Jay habitat.We used nest locations to define the extent of colonies at KAFB and WSMR, then collected geospatial data described in Appendix 1, Table A1.1, sampled at 100 m² from within the nesting colony boundaries and across the study area.To define biophysical habitat features important at the nesting colony scale for both sites, we developed and applied topographic and satellite-derived data available at colony-and-above scales (Appendix 1).We used 15 topographic predictors derived from digital elevation models, seasonal solar radiation, and vegetation indices derived from Landsat 5 Thematic Mapper satellite data acquired in 2005.

KAFB colony-scale predictive modeling
We performed a clustering analysis within Colonies KA and KB, using the ISODATA algorithm in ERDAS Imagine software (Hexagon Geospatial, ERDAS 2011), which creates an unsupervised classification based on pixel similarities.The unsupervised classification yielded 20 classes.
Some of the 15 geospatial layers were likely correlated with one another.To determine if a subset of the geospatial layers would be sufficient, we used the Jeffries-Matusita distance formula, (Swain and Davis 1978), which eliminated five layers while still maintaining high separability (mean = 1400.82,0 -1414 possible) between the 20 classes.Eliminating additional layers resulted in a significant drop in separability to 1200 or less.The five eliminated layers were: July NDVI (Normalized Difference Vegetation Index), January Moisture, October Moisture, Planiform Curvature (rate of elevation change perpendicular to slope), and Profile Curvature (rate of elevation change in direction of slope).
We then classified the study area, using a minimum-distance classifier.The supervised classification was based on statistics derived from the unsupervised classification of the nesting colonies.The result was a 20-class supervised classification based on the 10 geospatial measures.The 20 classes fell naturally into three groups with similar landscape position (low, medium, and high slopes) and spatial contiguity.To include the influence of both geospatial features and vegetation types, we merged the three geospatial classes with the three primary vegetation types, to make nine new classes that included vegetation and topographic information: high to moderate to low slopes in Piñon Pine Woodland, Piñon-Juniper Woodland, and Juniper Woodland and Savanna.
Next, we ran a neighborhood analysis for diversity over the entire study area, using a 32 X 32 matrix (of the 10 x 10 m cells), approximately equal to the area of our smallest known nesting colony.This analysis identified polygons with class diversity and area similar to that of known colonies (diversity values from 7 -10, out of 10).We converted the raster file containing diversity values ≥ 7 to vector format and calculated the percent cover value of Juniper Woodland and Savanna, Piñon-Juniper Woodland, Piñon Woodland, and no value for each polygon, for a total of 981 polygons.No-value polygons represent nontarget map classes such as shrubland.The percent cover values were input into the ArcGIS 10.2 (Esri 2013) Grouping Analysis tool to combine adjacent, categorically similar polygons.A spatial constraint allowed only contiguous polygon features of the same class that shared edges or vertices to be included in a group.The grouping analysis resulted in 15 polygons.We retained four additional, isolated but otherwise-qualified polygons eliminated by the Grouping Analysis tool.
We then created rules for retaining only polygons that conformed to the vegetation profiles of the two known colonies.These rules were the following: (1) retained polygons must contain ≥ 10%

WSMR colony-scale predictive modeling
We used the same approach to colony-scale modeling at WSMR.Because WSMR differed in topography and vegetation from KAFB, some steps used for KAFB did not work well for WSMR and were eliminated, as noted below.For WSMR, we started with 16 geospatial layers (Appendix 1, Table A1.1) and assembled the geospatial data within the study area.We attempted an unsupervised classification of Colonies WA and WB using ISODATA clustering with the same number of classes as at KAFB, but the resulting classification was highly pixelated and did not cluster well.We repeated the clustering analysis, specifying fewer classes until it produced a five-class classification that was not highly pixelated.The resulting matrix from the unsupervised classification performed within the two colonies (WA and WB) was noninvertible (invertibility is a requirement for the Jeffries-Matusita distance formula to reduce the number of geospatial measures).This left all 16 geospatial measures for the unsupervised classification.We then used the minimum-distance classifier to create a five-class, supervised classification of the study area outside the colonies, based on the statistics generated by the unsupervised classification of Colonies WA and WB.
A diversity focal analysis like that performed at KAFB was unsuccessful at WSMR.First, many WSMR nests were near the escarpment edge, which resulted in nesting colony-sized polygons that contained more nonhabitat than suitable habitat.Second, the WSMR colonies used for modeling were situated on a gentle, east-facing slope; hence, the geospatial measures (slope, aspect, elevation, etc.) showed very low variation relative to those at Of the five geospatial classes produced by the supervised classification, we retained the four having Pinyon Jay nests.We then merged the extent of the remaining geospatial classes with the land cover map of the study area and eliminated cells with vegetation types that did not contain Pinyon Jay nests, which were only found in Piñon Woodland.We eliminated small, isolated polygons and retained contiguous polygons that included all four nest-related classes.The resulting predictive colony-scale model for WSMR identifies areas having similar vegetation and geospatial characteristics to those of Colonies WA and WB.

Validation of nesting colony models
In April 2013, we field surveyed all except four small KAFB colony site polygons predicted by the model.We walked through each polygon, listening and watching for Pinyon Jays.We assumed we were able to hear jays vocalizing within 200 m and therefore covered each model polygon with walking transects approximately 400 m apart.For validation of the KAFB colonyscale model, we computed the percent overlap of the minimum bounding polygon of Colony KC (12 nests) with the KAFB colony model.For validation of the WSMR colony-scale model, we did the same for Colony WC (10 nests).The breeding-season home range polygon at KAFB covered 4418.9 ha, and the nonbreeding season home range was 4599.4 ha.The year-round range covered 5978.0 ha.

Home range-scale habitat use Kirtland Air Force Base home range-scale habitat use
Colony KA was situated in steep Piñon Woodland and Piñon-Juniper Woodland.Colony KB was centered on slopes of Piñon-Juniper Woodland extending north to an area of Juniper Woodland and Savanna (Fig. 3).Our banding and telemetry data indicate that individual jays moved between the north and south ends of the KAFB home range, and the two nesting colonies were occupied by members of the same nonbreeding season flock.We observed jays caching seeds in two areas.In early spring 2009, they retrieved caches in open habitat at the top of a previously burned hill (Fig. 3, north caching area).The south caching area (Fig. 3), used in both breeding and nonbreeding seasons, is a mostly south-facing, 18-ha, sparsely vegetated, gentle hillslope dominated by grasses and shrubs with scattered juniper trees.

White Sands Missile Range home range-scale habitat use
In The WSMR breeding season home range polygon covers 2041.7 ha (Fig. 3).The nonbreeding home range includes 2059.8 ha, and the combined breeding and nonbreeding range covers 3580.1 ha (Fig. 3).The WSMR flock left the nesting area (Colony WA and WB, Fig. 3) after the breeding season every year from 2005 to 2010.We collected several observations in September and October, but we do not know where they winter; hence, the breeding and nonbreeding home ranges combined are not yearround ranges for WSMR, as they are for KAFB.
At WSMR, we observed Pinyon Jays caching piñon seeds in August and September in only two years: 2006, when a moderately sized cone crop occurred, and 2011, when a small crop was produced.Birds cached in Piñon Woodland and Piñon-Juniper Woodland in both years but not in large, open areas as at KAFB (Fig. 3).

Home range-scale habitat models KAFB home range-scale habitat model
The land cover types within the KAFB study area, covering 9424.6 ha, include Ponderosa Pine Woodland, Piñon Pine Woodland, Piñon-Juniper Woodland, Juniper Woodland and Savanna, Shrubby Juniper Savanna, Upland Shrubland, Lowland Shrubland, Grassland, Burn, Built-up Land, and Road (Fig. 3).
A subset of the map units in the study area is primary habitat for Pinyon Jays at KAFB.This 7733.5-halandscape ranges in elevation from approximately 1742 m to 2439 m and includes Piñon Pine Woodland, Piñon-Juniper Woodland, and Juniper Woodland and Savanna.We have observed Pinyon Jays in the Ponderosa Pine Woodland east of KAFB, and it is likely that KAFB jays occasionally ranged into that habitat.However, we do not include Ponderosa Pine in this model, because we have no location data indicating that Pinyon Jays from the KAFB flock ranged into the Ponderosa Pine Woodland.

WSMR home range-scale habitat model
The land cover map for Pinyon Jays in the WSMR study area covers 22,465.9ha of potential habitat ranging in elevation from approximately 1800 to 2430 m.The land cover types in this map are Piñon Pine Woodland, Piñon-Juniper Woodland, Juniper Woodland and Savanna, Shrubland, Grassland, Built-up Land, and Road.Piñon Pine Woodland, Piñon-Juniper Woodland, and Juniper Woodland and Savanna cover 19,598.6 ha (Fig. 3).

Colony-scale predictive models KAFB colony-scale predictive model
The
The modeling process at WSMR unfolded quite differently from that at KAFB.Because the unsupervised classification resulted in a noninvertible matrix, evaluation of divergence among the classes prevented further reduction of the dimensionality of the data, and we were left with the original 16 layers for the classification.Because colonies on WSMR were much less diverse than KAFB colonies in terms of floristics and topography, attempts to apply focal analysis to aggregate classes based on diversity at the 10-ha scale were unsatisfactory.The resulting model includes large areas in which the four classes of interest were relatively equally distributed.Because we were not able to distinguish nesting colony sites from the surrounding Piñon Woodland based on different biophysical characteristics, we concluded that the entire area similar to the colony sites should be classified as two predicted nesting colony polygons, our study site and another more southerly site within the Oscura Mountains (Fig. 6).

Validation of colony-scale models
We found validation colonies during ground surveys in 2012.We found no previously unknown, active colonies during the surveys of the nesting colony model area.The validation colony at KAFB fell entirely within the colony-scale predictive model (Fig. 5).The validation colony at WSMR overlapped the WSMR colony-scale habitat model by 75.9% (Fig. 6).Balda (2002) noted that habitat use by Pinyon Jays had not been reported, and this has remained true until now.The lack of information on Pinyon Jay habitat use is surprising, given the long-term decline of this species (Sauer et al. 2014) and climate impacts to its habitat.Home range sizes, habitat types, differential habitat use within home ranges, and the approach to predictive colony-scale models presented here can provide a basis for conservation and management for Pinyon Jays in piñon-juniper habitats.

Home range-scale models
Pinyon Jay flock home ranges reported in the literature from various habitats vary widely: 1600 ha (Balda 2002), 2300 ha (Marzluff and Balda 1992), 2890 ha (Ligon 1971), and 6400 ha (Balda 2002).This variation among studies occurs in part because earlier estimates were not based on radio telemetry, and it is nearly impossible to closely follow a Pinyon Jay flock that is ranging widely.When satellite transmitters are developed that are small enough for a Pinyon Jay to safely carry, more accurate home range size estimates will be possible.However, even with better technology, home range sizes will likely vary because of differences in season, flock size, flock behavior, and habitat.For example, the KAFB breeding season home range is more than twice the area of the WSMR breeding season range, despite the smaller area of available piñon-juniper habitats at KAFB.This difference likely occurs because the KAFB flock nested at two widely-separated colony sites each year, compared to only one nesting colony site each year at WSMR.Balda (2002) presents a similar difference between home range sizes of two flocks.
For comparison, home range areas of breeding Clark's Nutcrackers, which are approximately 24% heavier than Pinyon Jays (Tomback 1998, Balda 2002), averaged only 200 ha, and those of nonbreeding individuals 400 ha (Schaming 2016).Like Pinyon Jays, Clark's Nutcrackers have a mutualism with a pine species, whitebark pine (P.albicaulis, Lanner 1996), but they are territorial breeders, much less social than Pinyon Jays, and eat a larger variety of pine seeds (Tomback 1998, Schaming 2016).The large sizes and various vegetation types characteristic of Pinyon Jay home ranges are probably a result of their mutualism with piñon trees, unpredictability of piñon mast crops, and perhaps sociality, although for sociality the causal relationship may be reversed.Other North American bird species that depend on mast also make large winter movements: Blue Jay (Cyanocitta cristata), Red-headed Woodpecker (Melanerpes erythrocephalus), and Black-capped Chickadee (Poecile atricapillus; Smith 1986).
Within their large home ranges, Pinyon Jays showed seasonal differences in habitat use, with lower-elevation Juniper Woodland and Savanna and Grassland (KAFB only) habitats used in winter and higher elevation habitat types with denser trees used for nesting.This pattern is more apparent at KAFB, where we have more nonbreeding-season occurrence data.Nonetheless, WSMR Pinyon Jays used high-elevation Piñon Woodland heavily in the early nonbreeding season, compared to KAFB jays.This result could be due to higher warm-season temperatures at WSMR at similar elevations.This example indicates that Pinyon Jay use of vegetation types within home ranges can be highly variable.
Radio telemetry enabled us to collect a substantial number of point locations over several months, but on several occasions, we were unable to locate some transmitters known to be functional.Our inability to find jays at these times could be due to the birds' mobility, their large home ranges, or the limitations of telemetry equipment in rough terrain.The home range for both study sites may therefore be an underestimate, particularly in winter when the species may wander widely (Marzluff and Balda 1992).Our home ranges may also under-represent the higher elevation vegetation types, because even with telemetry, it is difficult to cover the higher elevation habitats where roads are sparse or nonexistent.The home range polygons were based on occurrence data from years with no piñon mast crop.In the fall months of a mast year, Pinyon Jays would be expected to be harvesting seeds from large trees in piñon-dominated woodland.This would result in the inclusion of more Piñon Woodland in the nonbreeding home ranges.

Colony-scale models
Habitat use analyses for birds typically involve creation of linear statistical models incorporating ground measurements at nests (Drever et al. 2015), geospatial data at larger scales such as forest patches (Girvetz and Greco 2009), or a combination of data from different scales (Lawler and Edwards 2006, Doherty et al. 2010, Drever et al. 2015, Jedlikowski et al. 2016).Although existing multiscale models often incorporate land cover data, most focus on nesting habitat only (not year-round home ranges), and even data considered landscape scale are at smaller scales than those required for Pinyon Jays (Girvetz and Greco 2009, Doherty et al. 2010, Jedlikowski et al. 2016).
Because Pinyon Jays nest colonially, the number of available colonies for statistical analysis is limited, and a single nesting colony covers tens of hectares, encompassing wide variation in geospatial measures.A nesting colony therefore cannot serve as a single data point for statistical modeling, as could a nest or territory.These attributes of nesting colonies constrained our ability to conduct analyses within a linear modeling framework, and a different approach to modeling at the colony scale was necessary.We used geospatial modeling techniques and land cover models to classify the study area, then identified sites having suites of topographic and vegetation features similar to those of existing colonies.The resulting geospatial models predict where Pinyon Jays are expected to establish colonies.The usefulness of this approach is indicated by the validation colonies, which fell within the predictive model at both study sites.
The polygons in the KAFB colony-scale model include a mix of Piñon-Juniper Woodland and Piñon Woodland habitat (Fig. 5), whereas the nesting colony model at WSMR contains only Piñon Woodland (Fig. 6).As at the home range scale, differences between the nesting colony models at the two sites and the topography and vegetation variation within KAFB colonies indicate that suitable Pinyon Jay nesting colony habitat includes a fairly large range of vegetation types and topographies, even though colony models cover only 10-15% of the available habitat.
Colony-scale models at both study sites were based on nesting data from March-July, but these models likely apply to other months in which Pinyon Jays have been reported nesting.When Pinyon Jays nest in colder months, incidence of solar radiation at the nest tree is influential, because it can reduce energetic costs to nesting females (Marzluff 1988).Given that habitat choice for nest temperature occurs at the tree, rather than the colony scale, colony-scale attributes are probably applicable to winter nesting, but this question requires further investigation.

Scale considerations
This study highlights the importance of scale in habitat modeling.Although the two study sites are separated by a little over 1° of latitude (~135 km, Fig. 1), the elevation at these sites within Pinyon Jay habitat is similar, differing by less than 80 m at the lowest and 20 m at the highest elevations.A comparison of WSMR and KAFB colony-scale models suggests two useful conclusions.First, Pinyon Jay nesting colony habitat characteristics can vary broadly across large landscapes, as can within-site diversity of biophysical and vegetation features.Second, variation in topography and vegetation between and within study sites is only evident at the nesting colony scale.
What happens when a third, finer-scale model is considered?Nestscale models suggest that general rules for choosing nest sites may be similar between study sites, while habitat specifics differ.For example, Pinyon Jays at several colonies nested in larger-than average trees with denser cover over the nest, but the sizes and species of nest trees varied according to the available trees at each site (Johnson et al. 2014(Johnson et al. , 2015)).The only other studies on choice of nest site were conducted in urban ponderosa pines (Gabaldon 1979, Marzluff 1988), but results have some similarities to our nest-scale habitat study.Gabaldon (1979) found that Pinyon Jays also nested in taller ponderosa pine trees with denser foliage over the nest.

Conservation
Extensive Pinyon Jay home ranges mean that only land managers with jurisdiction over very large piñon-juniper landscapes can effectively manage for the year-round habitat needs of even one Pinyon Jay flock.Home range-and larger scale conservation must include various piñon-juniper subtypes for variable use among seasons: lower-elevation Juniper Woodland and Savanna (and sometimes Grassland) in winter and higher-elevation Piñon-Juniper and Piñon Woodlands for nesting and seed harvesting.
Even given recent piñon die-offs, suitable piñon-juniper nesting habitats currently appear to be abundant across the western USA.However, we found suitable nesting colony habitat in only about 10% and 15% of the study areas containing potential nesting habitat.As the quality of that nesting habitat declines with changing climate (K.Johnson, G. Sadoti, and J. Smith, unpublished manuscript), patches of suitable nesting habitat large enough to accommodate Pinyon Jay colonies are likely to become rare.If one model (McDowell et al. 2016) proves correct, most of the Pinyon Jay's nesting and foraging habitat in the southwestern USA will be lost by the year 2100.
In addition to affecting availability of nesting habitat, climate change may be impacting the size and frequency of mast crops (Zlotin and Parmenter 2008, Redmond et al. 2012) important to Pinyon Jay population viability (Marzluff and Balda 1992).Cone production in piñon trees is quite spatially and temporally variable under the best of circumstances (Lanner 1981).If climate change is affecting mast crops, habitat conservation on a large landscape scale, to account for spatial and temporal variation and reduction in mast crops, would be crucial for Pinyon Jay conservation.
Conservation of large landscapes for Pinyon Jays may be effective but is probably not feasible, especially under the multiple-use policies of most land management agencies in the USA.Understanding of a species' habitat needs at smaller scales offers potential for conservation action at more limited and feasible spatial and temporal scales.For example, thinning for fuels reduction within nesting colonies may negatively impact the quality of Pinyon Jay nest habitat (Johnson et al. 2014).Avoiding thinning or development in predicted nesting colony sites (just 15.35% of the study area mapped at KAFB and 10.3% at WSMR) would protect all used and potential nesting habitat at those study areas.Restricting other human impacts such as recreation or military activities for just three to four months of the breeding season within the colony-scale model would provide significant support for nesting jays.Piñon-producing areas and nearby caching sites could be set aside for Pinyon Jays and other wildlife during the fall season in mast years only.
Another climate-mediated cause of decline could be limited availability of free water (Gutzler 2013).Our colony-scale data reveal that Pinyon Jays used natural or artificial water sources near 10 nesting colonies we have studied.The jays established a nesting colony in atypical Juniper Woodland and Savanna habitat near our feeder, suggesting that Pinyon Jay nesting colony locations can be influenced by the presence of feeders.Marzluff and Balda (1992) also showed Pinyon Jays nesting near feeders in northern Arizona.Providing supplemental food and water near traditional nesting colony sites and elsewhere within flock home ranges is a simple, inexpensive, colony-scale management action that could ameliorate climate-related fluctuations in piñon seed and water availability and thereby increase Pinyon Jay climate resilience.Strategic placement of supplemental food and water could potentially assist managers to influence nesting colony location or movement.These colony-scale, spatially and temporally focused management actions would not be apparent from home range-scale models alone.
Responses to this article can be read online at: http://www.ace-eco.org/issues/responses.php/890

Fig. 1 .
Fig. 1.Scales of Pinyon Jay (Gymnorhinus cyanocephalus) habitat use.Home ranges (landscape scale) cover thousands of hectares.Nesting colonies (ecosystem scale) within home ranges cover tens of hectares; multiple nesting colonies can occur within one flock home range.Within nesting colonies, individual nest sites (not addressed in this paper) are smaller than one hectare; multiple nests occur within a nesting colony.Land cover data (background) over each study site were used to create both home range-and colony-scale models.

Fig. 3 .
Fig. 3. Pinyon Jay (Gymnorhinus cyanocephalus) breeding season, nonbreeding season, and year-round observations and home ranges at Kirtland Air Force Base (KAFB) and White Sands Missile Range (WSMR).Home ranges created by minimum bounding geometry (minimum convex) polygons.Background shows land cover classes.

At
KAFB from 2009 to 2011, we captured 66 Pinyon Jays and attached transmitters to 14. From the transmitters, we obtained 23 useful intersecting bearings resulting in 11 point locations (from 3 transmitters, 2 -6 locations each).Ground surveys yielded 126 additional point locations.Of the 137 bird locations, 103 (non-nests inside the mapped area) were used for assessing differential habitat use.We used 44 nest locations to delineate colonies (Colony KA = 22, KB = 10) and validate the nesting colony model (KC = 12 nests).

Fig. 4 .
Fig. 4. Concentration of use (CU) by Pinyon Jays (Gymnorhinus cyanocephalus) in breeding, nonbreeding, and combined season home ranges at Kirtland Air Force Base and White Sands Missile Range.CU = probability of use for a land cover type in the home range divided by the occurrence of that land cover in the home range.CUs for land cover types composing > 5% of home range polygons are presented.
At WSMR, combined breeding and nonbreeding Pinyon Jay observations are from Piñon Pine Woodland (87%, n = 62), Juniper Woodland and Savanna (5.6%, n = 4), and Piñon-Juniper Woodland (1.4%, n = 1).Breeding season use was heavily concentrated in Pinyon Pine Woodland: CU = 1.47 vs. 0.14 in Juniper Woodland and Savanna (Fig.4).Nonbreeding season CUs were: Pinyon Pine Woodland = 1.78,Pinyon-Juniper Woodland = 0.1, and Juniper Woodland and Savanna = 0.28 (Fig.4), indicating slightly increased use of lower-elevation Pinyon-Juniper Woodland and Juniper Woodland and Savanna in the nonbreeding season.CUs for both seasons combined were Pinyon Pine Woodland = 2.44, Juniper Woodland and Savanna = 0.14, and Pinyon-Juniper Woodland = 0.06.Numbers of detections were higher in the breeding season, because birds left the study area in the fall.Nonbreeding season CUs likely underestimate the use of Juniper Woodland and Savanna in the nonbreeding season at WSMR, because of the incompleteness of our spatial data for wintering Pinyon Jays at WSMR.

Fig. 5 .
Fig. 5. Kirtland Air Force Base Pinyon Jay (Gymnorhinus cyanocephalus) colony-scale predictive habitat model and nesting colonies.Colony model is shown in bright colors; other land cover is shown in muted colors.Geospatial and vegetation data from within the KA and KB Colonies were used to create the colonyscale model.Colony KC area was used for model validation.100% of KC colony area fell within predictive model.

Fig. 6 .
Fig. 6. White Sands Missile Range Pinyon Jay (Gymnorhinus cyanocephalus) colony-scale predictive habitat model.Inset shows nesting colonies.Colony model is shown in bright green; other land cover is shown in muted colors.Geospatial and vegetation data from within Colonies WA and WB were used to create the colony-scale model.Colony WC area was used for validation.Of the WC Colony area, 75.9% fell within the predictive model.

Home range-scale habitat use Kirtland Air Force Base home range-scale habitat use During
the spring and summer months from March 2009-May 2011, we surveyed for Pinyon Jays in piñon-juniper habitats (Juniper Woodland and Savanna, Piñon-Juniper Woodland, and Piñon Woodland) at KAFB.In 2009, we also surveyed in all other months to establish the nonbreeding range of the Pinyon Jays.We visited areas frequented by the birds and collected Global Positioning System (GPS) coordinates where we found jays (accuracy ±3-4 m).When detections were in restricted areas or we could not reach locations on hillsides, we recorded approximate locations on a map within the GPS unit.
Piñon-juniper habitats at KAFB occur primarily on the western slopes and bajadas of the Manzanita Mountains, a north-south chain that connects the higher Sandia Mountains to the north (elevation 3255 m) and Manzano Mountains (elevation 2802 m) to the south.KAFB has approximately 7733 ha of juniper and piñon-juniper habitats ranging in elevation from 1742 to 2439 m.The study area extends from lower-elevation juniper woodlands through mixed piñonjuniper woodland; piñon-dominated woodland; and, at highest elevations, ponderosa pine (P.ponderosa) woodland.Precipitation comes primarily during the summer months, in the form of heavy, short-duration thunderstorms.Annual precipitation varies from 20.3 cm in arid valleys and mesas to 76.2 cm in the Sandia Mountains (Kirtland Air Force Base 2007).At the Albuquerque Airport weather station, the average monthly temperature ranges from -4.7°C in January to 33.2°C in July (Western Regional Climate Center 2010).Fig. 2. Map showing study areas in New Mexico, USA: Kirtland Air Force Base (Manzanita Mountains) and White Sands Missile Range (Oscura Mountains).
Balda http://www.ace-eco.org/vol11/iss2/art6/2002).From a distance, we watched birds constructing nests and returned to the area after nesting was complete to collect GPS coordinates of known nests and search the colony area for additional nests, some of which remained from previous years (distinguishable from those of other species that nest at the study site).
Ligon 1977,e models were created in 2013 but employed nest location data from 2010-2012 and the land cover map created in 2010 at KAFB and nest location data from 2007, 2009, and 2010-2012 at WSMR.We located nesting colonies at both study sites by surveying for Pinyon Jays during April and May from 2007 -2012.We noted locations of jays showing breeding behaviors such as display flights, rattle calls, stick manipulation, silent sitting, or nest construction (Berger andLigon 1977,