Heterogeneity in avian richness-environment relationships along the Pacific Crest Trail

Predictions of the responses of montane bird communities to climate change generally presuppose that species and assemblages hold constant relationships to temperature across large study regions. However, comparative studies of avian communities exploring the factors that currently shape species richness patterns rarely analyze relationships across neighboring ecological regions of the same mountain chain. Evaluations of the intrinsic regional differences in species-environment relationships are needed to better inform expectations of how bird communities may be affected by future climate change. In this study, we evaluated the relative importance of three environmental factors (temperature, precipitation, and net primary productivity) in structuring avian richness patterns along a continuous mega-transect. We followed the route of the Pacific Crest Trail (PCT) (32.58° N to 42.00° N, ranging in elevation from 365 to 4020 m) on the California cordillera and completed avian point counts on 3578 systematically established survey plots. We divided this mega-transect into five sections, which corresponded to distinct ecological regions along the mountain chain. Regions differed both for elevation-richness patterns, exhibiting linear and unimodal trends, and for model-supported environmental drivers of patterns, with some richness-environment correlations changing sign across adjacent regions. These results were robust to sampling bias, regional species availability, and spatial autocorrelation. Although seasonal variation in avian movements may have limited influence on our results, we conclude that intrinsic regional environments affect bird species richness differently in each of these sections on the PCT, thus creating region-specific species-environment relationships. Appreciation of regional environmental heterogeneity will only increase in light of forecasted climate change, where regional predictions often diverge greatly from global trends, necessitating a site-specific approach to climate adaptation rather than ‘one size fits all’ strategies. Hétérogénéité des relations richesse aviaire-environnement le long du sentier des Crêtes du Pacifique RÉSUMÉ. Les prédictions des réactions des communautés d’oiseaux de montagne face aux changements climatiques supposent habituellement que les espèces et les assemblages d’oiseaux maintiennent une relation constante avec la température dans de grandes régions d’étude. Pourtant, les études comparatives de communautés aviaires qui explorent les facteurs à l’origine des profils de richesse aviaire analysent rarement les relations entre les régions écologiques voisines d’une même chaîne de montagnes. L’évaluation des différences régionales intrinsèques des relations espèces-environnement est nécessaire si l’on veut mieux prévoir comment les communautés d’oiseaux pourraient être affectées par les futurs changements climatiques. La présente étude évalue l’importance relative de trois facteurs environnementaux (la température, les précipitations et la productivité primaire nette) à structurer les profils de la richesse aviaire le long d’un mégatransect continu. Nous avons parcouru le sentier des Crêtes du Pacifique (en anglais Pacific Crest Trail, PCT) (du 32,58° N. au 42,00° N.; de 365 à 4 020 m d’altitude) sur la cordillère californienne et avons réalisé des points d’écoute dans 3 578 parcelles d’inventaire réparties systématiquement. Nous avons divisé le mégatransect en cinq sections qui correspondent à des régions écologiques distinctes sur la chaîne de montagnes. Les régions ont différé sur deux plans : d’une part, sur celui des profils altitude-richesse, qui montraient des tendances linéaire et unimodale et, d’autre part, sur celui des facteurs environnementaux qui agissaient sur les profils et que les modèles ont validés. Certaines régions adjacentes ont montré des corrélations richesse-environnement opposées. Ces résultats étaient robustes aux biais d’échantillonnage, la disponibilité aviaire régionale et l’autocorrélation spatiale. Même si les variations saisonnières des déplacements aviaires ont pu avoir un effet limité sur nos résultats, nous concluons que l’environnement régional intrinsèque affecte la richesse en espèces d’oiseaux de façon différente dans chacune des sections du PCT et crée ainsi des relations espèces-environnement régionalement spécifiques. L’estimation de l’hétérogénéité environnementale régionale ne s’améliorera qu’en fonction des prévisions relatives aux changements climatiques, sachant que les prédictions régionales divergent souvent grandement des tendances globales, ce qui explique que l’adaptation climatique requiert une approche spécifique au site plutôt que des stratégies « universelles ».


INTRODUCTION
Mountainous regions are hotspots of avian diversity (Ruggiero and Hawkins 2008) that are expected to be highly sensitive to climate change (La Sorte and Jetz 2010).In California, studies of ongoing climate change indicate that montane birds are already responding by shifting their elevation ranges (Tingley et al. 2012) and that long-term effects include richness declines and high community turnover (Tingley and Beissinger 2013).Despite evidence for location-specific sensitivity of species to different climatic factors (Tingley et al. 2009(Tingley et al. , 2012)), range predictions of montane birds, given climate change, generally assume species and communities respond universally to climate across their range (Sekercioglu et al. 2008, Stralberg et al. 2009, Wiens et al. 2009, Anderson et al. 2013).Some regional studies of avian speciesenvironment relationships appear to provide partial support of this assumption, particularly for large regions with relatively low topographic relief and homogenous habitat (DesGranges and LeBlanc 2012).On mountains, environmental variables change with elevation over relatively short distances.These steep gradients provide opportunities to conduct observational and comparative studies of the environmental drivers of avian diversity (Körner 2007).These types of studies can inform avian conservation efforts because they can identify both the relative importance of different environmental factors in driving species diversity as well as the universality, or regional heterogeneity, of environmental drivers of diversity among neighboring regions (Rowe 2009).Despite the potential contributions of studies in mountain systems to increase our understanding of avian ecology and conservation, systematic efforts that gather field data at a macroecological scale have been rare in these remote regions.Here, we survey avian richness continuously along the Pacific Crest Trail (PCT) through five distinct mountain regions that together form the California cordillera.We ask how speciesenvironment relationships change along the PCT through each of these sections, providing inference on how climate may differentially structure communities across adjacent ecological regions.
The topography of California is dominated by a massive northsouth oriented cordillera, which supports a rich avifauna because of the presence of strong, overlapping gradients in temperature (south to north with latitude and low to high with elevation) and precipitation (south to north and east to west with a rain shadow effect).Significant changes in temperature and precipitation are projected for California by the end of the 21st century.Temperature is projected to increase by 2 to 5°C across the state (Snyder et al. 2002) and although there is less agreement among climate models on the direction of precipitation change (IPCC 2007), there is more agreement that the proportion of precipitation that falls as rain will increase, while snowfall will decrease.These climatic changes may have a profound effect on birds and their habitats (Siegel et al. 2014), by altering plant productivity, vegetation structure, and plant species composition, including a projected loss of 50 to 90% of subalpine and alpine habitats by the year 2099 (Hayhoe et al. 2004).Some climaticinduced habitat changes are already occurring on California's mountains (Thorne et al. 2008, Dolanc et al. 2013).The magnitude of these climate projections makes more urgent the need for studies that clarify species-environment relationships, particularly in California's mountain regions.
Our ability to test for region-specific species-environment relationships in California is informed by the strong role that mountains have played globally in testing environmental hypotheses for elevation species richness patterns (Sanders and Rahbek 2012).Described patterns of elevation richness include increasing, decreasing, and most commonly, unimodal trends (Rahbek 2005), where maximum species richness occurs at some intermediate elevation.Ecologists and biogeographers have sought to elucidate species-environment relationships on mountains and found that many factors may correlate with avian richness, including past and present climate, geologic and evolutionary history, topography, biological interactions, habitat distribution, and human disturbance patterns (Lee et al. 2004, Ruggiero and Hawkins 2008, Sanders and Rahbek 2012).Perhaps the most widespread support can be found for temperature, net primary productivity, and precipitation as correlates of avian richness on mountains (McCain 2009).Temperature can limit bird distributions by exceeding physiological tolerances and both temperature and precipitation strongly limit plant growth and influence the amount of biomass available to consumers in an ecosystem.Although the relationship between temperature, precipitation, and productivity with avian richness has been examined at the global scale (Hawkins et al. 2003, McCain 2009) and several energy-related underlying mechanisms for this relationship have been proposed and tested (Evans et al. 2005, Honkanen et al. 2010), it remains unknown how these three factors differentially relate to avian richness along latitudinal and elevation gradients among adjacent mountain regions.
Given that region-specific environmental characteristics can differentially correlate with richness (Rowe 2009), our goal was to examine whether and how species-environment relationships change across adjacent regions.We surveyed local avian richness on the PCT through five adjoining mountain regions.Our objectives were (1) to quantify avian richness patterns along elevation gradients on the PCT in each region and (2) to compare the relative strength of these three factors as correlates with richness among the five regions, while simultaneously accounting for several inferential challenges, including sampling effects, spatial autocorrelation, and the regional species-pool.We used our observations to test the null hypothesis, typical of climate change predictions, that species richness-environment relationships show generality across neighboring ecoregions.

Study area
The PCT extends across a series of mountain ranges from the international borders with Mexico to Canada, across the states of California, Oregon, and Washington.The California section (~2736 km, 32.58° N to 42.00° N) served as a mega-transect along which we surveyed birds at regular intervals (Fig. 1a).The PCT is a recreational hiking and equestrian trail through remote habitats, spanning elevations from 365 m in southern California to 4020 m in the southern Sierra Nevada.Along the California cordillera, total annual precipitation generally increases from south to north (Fig. 1c).It follows a contouring route along the relatively wet and windward, western-facing slopes of the cordillera, and to a lesser extent, it also occurs on drier and leeward, eastern-facing slopes.In southern California, it crosses portions of the Mojave and Colorado Deserts.Thus, the PCT route encounters a wide range of habitat types, including desert scrub, grassland, chaparral, woodland, forested, and alpine habitats (California Department of Fish and Wildlife 2008).We divided the mega-transect into sections that corresponded to five, climatically and geologically distinct, ecological regions: southern California (SoCa), southern Sierra Nevada (SoSN), northern Sierra Nevada (NoSN), southern Cascades (Casc), and Klamath Mountains (Klam).With some modification, we defined our sections based on the boundaries for California ecological regions defined by Miles and Goudey (1997).The sampled length of the elevation gradients varied among these sections and was determined by the elevation profile of the PCT route itself through each region as well as the overall topographic relief of each mountain range.

Bird surveys
Birds were surveyed according to a spatially and temporally standardized protocol at survey plots on a single visit.One observer (MCM) systematically stopped at 10-minute walking intervals (approximately, every 500 -700 m) along the PCT and conducted 5-minute point counts of birds on fixed-radius, 50-m circular plots.All bird species and individuals seen or heard within each plot were counted.A total of 3578 point counts were completed across the 2736 km length of the PCT in California in a single field season (2006), resulting in local richness estimates for each of the plots.A second observer (AMM) assisted with data recording and conducted a rapid habitat assessment.
We followed the methodology for fixed-radius point counts (Ralph et al. 1995), except in two key respects.First, instead of restricting counts to the morning hours, surveys were conducted from dawn to dusk.Second, we surveyed from 2 April to 8 September 2006, when most California migrant birds were present in summer habitats, but before the return of migrants to wintering grounds (McGrann and Thorne 2014).Surveys along each of the five regional sections of the PCT were completed within a narrower time frame: SoCa (2 April -27 May), SoSN (29 May -20 June), NoSN (8 August -8 September), Casc (26 June -20 July), and Klam (20 July -4 August).Although surveying throughout the day and throughout the summer served to maximize data collection in remote landscapes, several montane birds in California have complex diurnal and seasonal behavior and movements (McGrann and Thorne 2014).We explicitly accounted for time-of-day and time-of-year when building our models because these factors may influence the number of species detected.

Environmental data
We used the spatial analyst extension in ArcGIS (version 10, Environmental Systems Research Institute, Redlands, California, USA) to extract elevation, temperature, precipitation, and productivity variables from gridded data sources at survey points.Elevations came from the National Elevation Dataset (10 m resolution; U. S. Geological Survey, http://ned.usgs.gov).Spatial grids (1 × 1 km, or 30 arc-second resolution) of temperature and precipitation for 2006 were obtained for California using the Parameter-elevation Relationships on Independent Slopes Model (PRISM).PRISM data are particularly good at predicting the diverse climates across mountainous regions (Daly et al. 2008).From the spatial grids, we extracted average annual maximum, minimum, and mean temperature (°C) and total precipitation (mm) in 2006 for each point count location.We used the satellite-derived Enhanced Vegetation Index (EVI) and Normalized Difference Vegetation Index (NDVI) as proxies of net primary productivity.Spatial grids (250-m resolution) of 16-day composites of EVI and NDVI for 2006 were derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) and obtained from the Land Processes Distributed Active Archive Center (http://lpdaac.usgs.gov).NDVI is calculated from the reflectance of red and near-infrared spectral bands, which is associated with healthy green vegetation.EVI includes adjustments for soil and atmospheric aerosol scattering and both EVI and NDVI serve as good proxies of plant biomass (Jensen 2007).For each bird point count location, we extracted two EVI and NDVI variables: the specific value from the 16-day composite grid that corresponded with the date of the point count survey, and the maximum value over the entire field season in which the surveys were conducted.

Models of richness
We tested the relative importance of temperature, precipitation, and net primary productivity variables as predictors of observed bird richness using Generalized Linear Models (GLM) with a negative binomial distribution to account for overdispersed count data.All analyses were conducted in R version 3.0.3(R Development Core Team 2014) and all models were run using the package "MASS" (Venables and Ripley 2002).Variable selection was conducted during an initial exploratory phase, followed by multimodel comparison.In all cases, model or variable support was measured by the Akaike Information Criterion (AIC, Burnham and Anderson 2002).Although attempts were made to reduce multicollinearity of variables in models, model coefficients of all top models were also inspected for variance inflation and stability in the presence of other covariates.
In the exploratory phase, richness models were constructed to compare support for different predictor variables within the same environmental hypothesis class, i.e., temperature, precipitation, and net primary productivity.Within each hypothesis class, data from the entire study area were combined to test the relative importance of each variable across all five regions.Temperature models compared average annual minimum temperature (T min ), maximum temperature (T max ), and mean temperature (T mean ), as well as squared effects.Precipitation models compared total annual precipitation (P) and squared effects.Productivity models compared maximum EVI (EVI max ) and NDVI (NDVI max ) for the entire field season and EVI (EVI t ) and NDVI (NDVI t ) on the survey date as well as squared effects.Each exploratory model contained a single environmental predictor variable (or two variables, for squared effects).In addition to the temperature, precipitation, or net primary productivity variables, all exploratory models contained time-of-year, time-of-day (including squared effects), and regional species pool as independent variables.By accounting for time-of-day and timeof-year explicitly within models, we acknowledge that these observation data attributes may affect the number of species observed at a single survey point.Regional species pool was included in models because the available pool of species is generally accepted to structure local richness estimates (Ricklefs 2000).Following White and Hurlbert (2010), the regional species pool at each location was estimated by intersecting survey points with range maps from the California Wildlife Habitat Relationships System (California Department of Fish and Wildlife 2008) for all 155 detected bird species (Appendix 1).Independent variables were normalized to a mean of zero and a standard deviation of one.Within each environmental class of variables (temperature, precipitation, or net primary productivity), the individual variable parameterization that was best supported by the data (lowest AIC) was used in the following multimodel comparison to understand region-specific differences in richness relationships.
In the second analysis phase, richness was modeled as a function of the three environmental hypothesis classes for each of the five regional sections on the PCT.Eight GLM models were constructed, which independently and additively tested the three best-supported parameterizations from the exploratory analysis, as well as a null model, i.e., no environmental covariates.All models contained the same sampling-based independent variables (time-of-day, time-of-year, and regional species pool) as the exploratory analysis and used a negative binomial error distribution with a log-link.
To investigate how the three environmental hypothesis variables may be differentially supported among each of the five regions, region-specific models were compared using AIC and AIC weights (w i ).Impact factors (IF), which are cumulative AIC weights in support of an individual variable, were calculated for each of the three classes of variables for each region, to test for shifting support of temperature, precipitation, or net primary productivity predictors across all five regions of the PCT (Burnham and Anderson 2002).IF scores of 1 indicated universal support, scores above 0.8 indicated strong support, and scores below 0.5 indicated weak support.
Because differing support for environmental relationships within regions may arise from incomplete overlapping of sampled environmental gradients across regions, a third analysis phase focused on region-specific environmental relationships within a model that included data from the entire Californian megatransect.An initial, baseline model included temperature, NDVI, and precipitation, matching the "full" model for the second phase of the analysis but combining data from all regions.Eight subsequent models tested for region-specific differences of environmental relationships within this baseline model: (1) a model with region added as a categorical variable; (2-4) models that separately tested for interactions between each region and temperature, precipitation, and NDVI; (5-7) models that tested for interactions between region and two of the three environmental variables; and (8) a model that included interactions between each region and all three environmental variables.Similar to previous analysis stages, all models were compared with AIC and the top model was examined for significant covariate effects.
To explore relationships of elevation with richness in our dataset, we estimated observed elevation richness maxima via GLMs.For each region, we modeled observed richness as a Poissondistributed function of elevation, including quadratic effect.Elevations of maximum richness were identified from these modeled relationships.

Tests of spatial autocorrelation and inferential challenges
We tested the degree to which spatial autocorrelation was present in both our raw richness values and our model residuals (White and Hurlbert 2010).For each of the five ecological regions, we created spatial correlograms by calculating Moran's I (Zuur et al. 2007) with distance lags of 50 km.Significant correlations (p < 0.01) were determined through resampling 1000 permutations.Analyses were conducted in R using the package "ncf " version 1. 1-4 (Bjornstad 2012).
In addition to spatial autocorrelation, a number of other factors can influence counts of local richness, potentially biasing parameter estimates.In many cases, patterns derived from observed richness may differ from true patterns because of sampling bias (Gotelli and Colwell 2001).Although true richness could not be estimated because multiple samples were not collected for most sites, we accounted for potential sampling bias directly by including sampling variables in our models, e.g., timeof-day, time-of-year.Additionally, in a subsequent year (2007), a subsample of high-elevation sites were revisited and surveyed in three consecutive increments for 120 minutes total.This supplemental data allowed a partial evaluation of the completeness of the 5-minute surveys (Appendix 2).

RESULTS
Bird surveys detected 155 bird species, comprising 21,599 individuals and totaling 13 orders, of which 68% of the species belonged to Passeriformes.The exploratory phase in the analysis revealed strongly supported richness parameterizations (Table 1) within each of the three classes of environmental variables (temperature, precipitation, and productivity).T max with squared effect was the best supported temperature model (w i = 1.0), while P (linear-only) was the best supported precipitation model (w i = 0.73), and NDVI max with squared effect was the best supported productivity model (w i = 0.70).These three parameterizations were used in subsequent model comparisons of richness relationships.222.13 0.00 † Environmental hypothesis parameters were tested using both linear and quadratic forms.Within quadratic models (signified by ²), linear terms were also included.‡ All models fit using Generalized Linear Models (negative binomial error distribution, log-link).Other model covariates included time-ofday (including squared), day-of-year (including squared), and size of regional species pool.
Species richness models that tested multiple environmental hypotheses simultaneously showed support for multiple factors as correlates with richness that changed inter-regionally (Tables 2, 3).In the SoCa, Casc, and Klam regions along the PCT, temperature was strongly supported as a factor structuring richness patterns.The relationship of richness with temperature for these regions was generally negative: as temperatures decreased with increasing elevation, richness increased.In the two Sierran regions as well as SoCa, productivity (NDVI max ) received full model support with maximum richness associated with values  negative in the adjoining regions of NoSN (β P = -0.22 ± 0.07) and Casc (β P = -0.16± 0.03), yet was positive in Klam (β P = 0.12 ± 0.04), despite the three regions having similar variation in precipitation across the sampled gradient (Figs. 1c,2c).
Within a single model using data from all regions, multimodel analysis confirmed that there were significant region-specific relationships of richness with each environmental variable.Among the nine models examining region-specific interactions, the top model contained interaction effects for all three variables: T max , NDVI max , and P (Table 4).Close examination of this particular model showed significant interactions, relative to a baseline contrast of Casc, of SoCa and SoSN with T max , of Klam with P, and of SoCa with NDVI max (Table 5).The modeled relationships indicated that the two southernmost regions showed increases at high T max while the two northernmost regions had decreases with T max (Fig. 2A), that SoCa had richness that increased at high NDVI max unlike other regions (Fig. 2B), and that richness increased with precipitation in Klam, whereas it decreased for Casc (Fig. 2C).Although not all regions were sampled on the exact same extent of environmental gradients, significant region-specific interactions contrasted with at least one other region that did overlap in environmental space (Fig. 2).212.9 0.00 † Environmental parameters signified by ² also include linear terms.‡ All models fit using Generalized Linear Models (negative binomial error distribution, log-link).Other model covariates included time-ofday (including squared), day-of-year (including squared), and size of regional species pool.
Observed avian richness exhibited region-specific variance in form (Fig. 3).Unlike modeled relationships with environmental gradients, which were generally quadratic in form (Fig. 2), relationships with elevation within each region were monotonic with the exception of SoSN (Fig. 3).For three regions, the trend of observed richness with elevation peaked at the highest sampled sites (SoCa = 3194 m, Casc = 2312 m, and Klam = 2331 m).The two Sierra Nevadan regions showed richness maxima in the middle and lower-middle of their elevation sampling profiles (SoSN = 2425 m; NoSN = 2117 m).
Our top richness model, including region-specific interactions, showed little effect of spatial autocorrelation on residuals (Fig. 4).In comparison, raw richness values showed significant patterns of spatial autocorrelation.Although the richness model reduced spatial autocorrelation at all distance classes, significant correlation in residuals remained at 50 km, and between 600 -800 km, the scale at which biogeographical differences can influence richness.Overall, magnitudes of correlation coefficients were low, however, with no model of residuals having a significant correlation exceeding ± 0.1.0.058 0.175 † For regional effects, Casc was used as the baseline versus which all other regional contrasts were compared.Lines are bolded where parameter-specific P-values were less than 0.05.‡ Klam: Klamath Mountains, NoSN: northern Sierra Nevada, SoCa: southern California mountains, SoSN: southern Sierra Nevada, Casc: southern Cascade Range.

Environmental relationships
We found that avian elevation richness patterns along the PCT have different correlative structures with environmental factors for each of the five montane ecoregions in our comparison (Fig. 2, Table 3).These environmental relationships were independent of generally similar elevation-richness patterns that show midelevation peaks in 4 out of 5 regions between 2000 and 2500 m (Fig. 3).The mega-transect route captures a wide range of environmental gradients within each region allowing us to elucidate how certain environmental factors may be strongly related to richness in some regions but not others.Temperature, for example, had universal model support in the three regions with the lowest elevations (SoCa, Casc, and Klam).Physiological adaptations to specific temperature ranges may limit the number of species able to tolerate extremes at either ends of this temperature gradient (Hawkins et al. 2003, McCain 2009).Interestingly, this relationship may differ across neighboring regions depending on community composition and communitylevel thermal affinities.The comparison of Casc and Klam to SoSN and SoCa are important cases.All four regions extensively sampled hot, low-elevation areas (Fig. 1b), yet the modeled association of richness with temperature showed a significantly positive relationship for SoCa (and a marginal one for SoSN) but a negative relationship with the northern regions of Casc and Klam (Fig. 2a).Thus, although temperature is an important correlational driver of richness in three regions, the nature of that relationship differs.Net primary productivity can also limit species by influencing the vegetation composition and structure of habitats as well as available food resources for birds (Evans et al. 2005).In our results, NDVI held a strong relationship with richness in the three regions that surveyed the highest elevations (SoCa, SoSN, NoSN).This relationship held a primarily quadratic, i.e., unimodal, form, although SoCa, and to a lesser extent SoSN, showed divergence from this trend with a failure to level off in richness at high NDVI values (Fig. 2b).The southern California mountains, with their shallow precipitation gradient and high aridity, appear to not show a saturation effect of richness with increasing NDVI.
Of the environmental gradients explored, precipitation showed the most interesting region-specific relationships.Although precipitation received strong model support in three regions (NoSN, Casc, Klam), the direction of the relationship reversed across the Casc-Klam boundary (Fig. 2c), from positive to negative.It is worth noting that the gradient of precipitation sampled was nearly identical in both Casc and Klam.Precipitation is often used as a surrogate for net primary productivity (Mittelbach et al. 2001), but in certain contexts, precipitation may not adequately reflect productivity, particularly within climates where water is not a limiting factor, or across a wide range of climate types (Cusens et al. 2012).The flipped relationship between Casc and Klam could reflect such changes, or could represent the changing nature of precipitation as one moves out of the rain shadow to the east, and toward the coast, across the ecoregional boundary (Fig. 1a,c).It is worth noting that we did not account for the form of precipitation, i.e., rain versus snow, or environmental metrics, such as climatic water deficit, which combine precipitation with other aspects of the environment and which have been shown to be important in structuring changes to plant distributions over time (Crimmins et al. 2011).
Models forecasting the effects of climate warming on species occurrences and richness generally assume that species will respond similarly to climate across their range (Sekercioglu et al. 2008, Stralberg et al. 2009, Wiens et al. 2009, Siegel et al. 2014).Although not all regions in our study area are sampled along the exact same extent in environmental space, regionspecific interactions contrasted with at least one other region where environmental gradients did overlap (Figs.2a-c).Taken together, these results indicate subtly changing correlates of local avian species richness as one heads north through the mountain ranges traversed by the PCT.This lack of commonality in pattern further reinforces Rowe's (2009) findings on small mammal richness, where it was also observed that no single environmental factor best explained patterns among neighboring mountain regions.Such region-specific relationships with species richness suggest that models that predict diversity changes over time, as with climate change, may not capture local-or region-scale environmental interactions if conducted at the continental scale.Although we found interactions at the scale of approximately 1.5-2° latitude, it is unknown whether such interactions are more generally expected at larger or smaller scales.Understanding the scale at which diversity interacts with environmental gradients appears to be an important venue for future research.

Addressing inferential challenges and alternative hypotheses
In any study of drivers of richness, it is critical to account for the many factors that can influence either the formation or the detection of patterns.We explicitly examined sampling effects because of the timing of the surveys, namely time-of-day and time-of-year, and found strong relationships (Fig 2d-f).Observed richness was highest in the early morning hours when bird singing activity was highest.Observed richness also increased toward the end of the field season, perhaps because of increased observer familiarity or the detection of juvenile birds.We show in a separate assessment of this dataset that of 74 common species, just 12 species may have been detected in migration in certain regions (McGrann and Thorne 2014).Despite these relationships with the timing of the surveys, the influence of the environmental variables was still discernable after accounting for these sampling effects in our models.
The species-area relationship may also play a role in shaping species richness patterns.As area generally decreases in a mountain range with increasing elevation, richness is expected to decrease as well.Gamma diversity, or region-wide richness, may be highly influenced by area, but alpha diversity, particularly at small spatial units of this study, is less likely to be influenced by area (McCain 2007).Nevertheless, we explicitly accounted for the available regional richness (Appendix 1) at each site in our models and did not find the relationship to be significant for three of the five regions.Thus, we did not find a strong or consistent effect of the regional richness pool on our observed values of local richness.
We also acknowledge that temperature, precipitation, and net primary productivity are not the only forces potentially shaping diversity patterns along the PCT.Habitat alteration (Forister et al. 2010) has played a significant role in restructuring butterfly communities along the California cordillera.In contrast, the PCT traverses many of the most remote and intact habitats left in California.Although the trail crosses busy highways, our survey methodology ceased counts along the very short sections of trail that passed through anthropogenically developed areas.

CONCLUSIONS
Across our five ecological regions along the PCT on the California cordillera, we found that three environmental factors (temperature, precipitation, and net primary productivity) varied in importance as correlates of avian richness.Avian richness responded differently along elevation gradients in the semiarid climate of Southern California, in the more mesic climates of the Klamath and Cascade ranges, and in the two high-elevation Sierra Nevada regions.We found no single pattern that explained correlative structures of environment and avian richness across latitude, and some correlations between richness and environmental factors had opposite signs in adjacent regions (e. g., Casc and Klam).
Montane birds are particularly sensitive to projected shifts in temperature and precipitation, and the area-elevation relationship of mountains is expected to leave less room for birds to track the climate trends of the future if increasing temperatures shift them upslope (Sekercioglu et al. 2008, La Sorte andJetz 2010).Evidence has demonstrated, however, that elevation shifts of the recent past have been both upslope and downslope as bird species locally track temperature and precipitation changes, often in directions that show regional-scale variability in associations with temperature versus precipitation (Tingley et al. 2009(Tingley et al. , 2012)).Our results of current richness-environment relationships support this observation that local-or regional-scale climatic associations can strongly structure bird communities.
Our study area engages a wide range of potential stakeholders and actors for regional conservation efforts.The ownership of lands along the PCT is diverse, including private lands, state and federal parks, national forests, and wilderness areas.Our results suggest that a "one size fits all" approach to climate change adaptation for birds on mountains is inappropriate.If climatic associations with communities can switch over short geographic distances, then species distribution modeling, conservation planning, species prioritization, and adaptation strategies are best done at the local to regional scale with coordination among agencies, conservation organizations, and private citizens (Heller and Zavaleta 2009), and not applied uniformly across broad areas, even contiguous mountain chains.We suggest that it is critical to consider the intrinsic local-and regional-scale environmental factors, such as the specific climate of a mountain region, in conservation efforts anticipating the effects of future climate change on avian communities.

Fig. 1 .
Fig. 1.Relationships of environmental variables with latitude and elevation across the Pacific-Crest Trail (PCT) within California.(a) Map showing the extent of the mega-transect in California and locations of the five regional sections of the PCT.The black line indicates the route of the PCT.Graphs depict underlying environmental gradients in (b) average maximum annual temperature (Tmax), (c) total annual precipitation (Precip), and (d) maximum seasonal Normalized Difference Vegetation Index (NDVI max ) across all point count sites.Dotted lines in graphs indicate latitudinal cutoffs for each region.SoCa: southern California, SoSN: southern Sierra Nevada, NoSN: northern Sierra Nevada, Casc: southern Cascades, and Klam: Klamath Mountains.

Fig. 2 .
Fig. 2. Modeled relationship and effect size of environmental (A-C) and sampling (D-F) factors associated with observed avian richness at point counts across each sampled region along the Pacific-Crest Trail.Relationships derived from each region's best supported Generalized Linear Models (GLM) and show parameter means while controlling for mean conditions of all other covariates.SoCa: southern California mountains, SoSN: southern Sierra Nevada, NoSN: northern Sierra Nevada, Casc: southern Cascade Range, and Klam: Klamath Mountains.

Fig. 3 .
Fig. 3. Relationships of observed avian richness with elevation, by Californian region of the Pacific Crest Trail.Elevation trends of richness (solid black lines) were modeled for each region separately with the elevation at which richness peaks (dotted black lines) were extracted from these models.Klam: Klamath Mountains, Casc: southern Cascade Range, NoSN: northern Sierra Nevada, SoSN: southern Sierra Nevada, SoCa: southern California mountains.

Fig. 4 .
Fig. 4. Patterns of spatial autocorrelation (Moran's I) for raw richness data (red) and residuals of top model including environmental and observational factors (blue).Solid circles represent significant correlations (p < 0.01, Bonferroni corrected for number of distance classes).Dashed line represents zero autocorrelation.

Table 1 .
Parameter selection for each of three environmental hypothesis categories, ranked by AIC weight (w i ).

Table 2 .
Individual model AIC weights, by region, for additive and single-variable models testing the ability of different environmental hypotheses to explain observed avian richness.

Table 3 .
Impact factors † (cumulative AIC weights) indicating weight of evidence in support for different richness-structuring hypotheses for each region.

Table 4 .
Multimodel comparison of different combinations of models that test for interactions between environmental variables and Sierra Nevadan regions.

Table 5 .
Parameter estimates, standard errors, and statistical significance for the top model explaining observed bird richness as a function of environmental factors and regional interactions.