Evaluating within‐population variability in behavior and demography for the adaptive potential of a dispersal‐limited species to climate change

Abstract Multiple pathways exist for species to respond to changing climates. However, responses of dispersal‐limited species will be more strongly tied to ability to adapt within existing populations as rates of environmental change will likely exceed movement rates. Here, we assess adaptive capacity in Plethodon cinereus, a dispersal‐limited woodland salamander. We quantify plasticity in behavior and variation in demography to observed variation in environmental variables over a 5‐year period. We found strong evidence that temperature and rainfall influence P. cinereus surface presence, indicating changes in climate are likely to affect seasonal activity patterns. We also found that warmer summer temperatures reduced individual growth rates into the autumn, which is likely to have negative demographic consequences. Reduced growth rates may delay reproductive maturity and lead to reductions in size‐specific fecundity, potentially reducing population‐level persistence. To better understand within‐population variability in responses, we examined differences between two common color morphs. Previous evidence suggests that the color polymorphism may be linked to physiological differences in heat and moisture tolerance. We found only moderate support for morph‐specific differences for the relationship between individual growth and temperature. Measuring environmental sensitivity to climatic variability is the first step in predicting species' responses to climate change. Our results suggest phenological shifts and changes in growth rates are likely responses under scenarios where further warming occurs, and we discuss possible adaptive strategies for resulting selective pressures.

plasticity, or going locally extinct (Lande & Shannon, 1996;Parmesan, 2006;Sinervo et al., 2010). While there is evidence that some species may be able to track their climate niche through time (Tingley, Monahan, Beissinger, & Moritz, 2009), it is unclear how less mobile species will respond to changing conditions. For dispersal-limited species, including many amphibians (Gibbons et al., 2000; but see Smith & Green, 2005), environmental plasticity and evolution are crucial components of adaptive change because potential to respond in the near-term through sufficient movement is limited, especially in fragmented landscapes (Cushman, 2006;Ruiz-Aravena et al., 2014). Without dispersal, persistence will depend on the interplay of local demographic responses to climate and the degree to which negative responses can be minimized by plasticity and evolutionary adaptation.
At a basic level, predicting population persistence under climate change requires an understanding of the degree to which demography-growth, abundance, survival, recruitment, and emigration/immigration (Hanski & Gilpin, 1991;Thomas, 2000)-is influenced by environmental factors. Connecting climate change to shifts in demographic rates can be challenging (McCain, Szewczyk, & Bracy Knight, 2016), but there is growing evidence that climate change can alter vital rates for the worse, resulting in increased risk of population extinction (Barbraud & Weimerskirch, 2001;Rodenhouse, Christenson, Parry, & Green, 2009;Rohr & Palmer, 2012). For instance, Bestion, Teyssier, Richard, Clobert, and Cote (2015) showed experimental warming reduced adult survival in a common European lizard species.
This resulted in lizard life history favoring earlier production of offspring in face of reduced life span. Despite this adaptive life-history shift, a large portion of the species' populations are still predicted to go extinct before mid-century (Bestion et al., 2015).
Phenotypic plasticity can mediate the response wildlife have to anthropogenic stressors (Hendry, Farrugia, & Kinnison, 2008). Behavioral plasticity, particularly in the timing, patterns, or extent of activity, is one of the most rapid phenotypic responses to novel conditions. These responses immediately affect the environmental conditions to which an individual is exposed (Snell-Rood, 2013;Wong & Candolin, 2015). When responses are adaptive, behavioral plasticity can help a species or populations avoid the consequences of climate change. For example, some turtles are able to shift nest-site selection to avoid detrimental warming conditions (Refsnider & Janzen, 2012). Behavioral responses can minimize short-term impacts and may provide a mechanism for future adaptation via generation of novel traits (Gomez-Mestre & Jovani, 2013;Zuk, Bastiaans, Langkilde, & Swanger, 2014).
The red-backed salamander (Plethodon cinereus), a common North American woodland salamander, is an ideal model to examine these two components of adaptive capacity-demography and behavioral plasticity-in a dispersal-limited species. Plethodon cinereus metapopulations exhibit minimal genetic exchange (Cabe et al., 2007;Marsh, Page, & Hanlon, 2008;Marsh et al., 2007). They are also able to modify the environmental conditions to which they are exposed both through horizontal movement for surface microhabitat selection and vertical movement between the surface and underground refugia (Heatwole, 1962;Jaeger, 1980;Spotila, 1972;Taub, 1961).
To evaluate within-population variability, P. cinereus have a genetically inherited (Highton, 1959) color polymorphism that has previously been tied to differences in environmental tolerance. Although the direct mechanism remains elusive, this relationship between color and climate niche may be due to pleiotropy or linkage disequilibrium, but even the exact mode of inheritance is still unknown (Highton, 1959).
Perhaps as a result of these differences, the frequency of both morphs varies geographically, with striped morphs more common in cooler, wetter regions and the lead-backed more common in warmer, drier regions (Fisher-Reid, Engstrom, Kuczynski, Stephens, & Wiens, 2013;Gibbs & Karraker, 2006;Lotter & Scott, 1977). Collectively, these studies suggest striped morphs prefer cool-wet conditions and that leadbacked morphs can better tolerate warm-dry conditions. However, these climate-morph relationships have not been consistent (Petruzzi et al., 2006), and recent work has criticized the use of the color polymorphism for understanding climate relationships (Moore & Ouellet, 2015).
We use P. cinereus as a model organism for investigating withinpopulation variation in demography and behavioral plasticity in response to environmental conditions ( Figure 1). As a proxy for withinpopulation variation in climate tolerance, we use the color polymorphism described above. Characterizing and understanding within-population variation, particularly for traits tied to climate tolerance, should allow us to better understand the adaptive capacity of the species. Our goals were (1) to determine the extent to which P. cinereus surface activity and demography are impacted by environmental variation in temperature and precipitation and (2) to evaluate the validity of the color polymorphism as a mechanism for illuminating within-population variation in climate response. We focus on the interaction between demographic and behavioral responses to environmental conditions. By simultaneously investigating plasticity and demography, we can improve predictions of how a species might adapt through within-population variation, and we can use this relationship to determine how populations might be impacted by predicted climate change.

| Data collection
Between October 2009 and May 2013, we conducted capture-markrecapture surveys for P. cinereus at the Patuxent Wildlife Research Center (Laurel, MD, USA). Three plots were established >20 m apart in lowland-deciduous hardwood forest under similar canopy conditions. Plots contained an array of cover boards (30.5 × 30.5 × 2.54 cm pieces of rough-cut pine) spaced at 1-m intervals, allowing us to effectively monitor P. cinereus populations and movement (Miller Hesed, 2012). Plot I was 20 × 20 m (400 cover boards), and plots II and III were 10 × 10 m (100 cover boards each). Captured salamanders were given individually identifying marks with visual implant elastomer, a technique that provides easily interpretable, long-lasting marks (Gillette & Peterson, 2001;Grant, 2008). Gender, color morph (striped or lead-back), and snout-to-vent length (SVL) were recorded every encounter and were independently determined twice to account for observer error. Environmental conditions (i.e., temperature, rainfall) were gathered from a weather station less than one kilometer away (2009)(2010)(2011)(2012) or from a weather station three kilometers (2012-2013) from the study area. Although this reduced the resolution of our environmental data, the magnitude and direction of changing environmental conditions are highly correlated at such small spatial scales. We opportunistically surveyed plots 3-9 times each spring and autumn, ensuring a minimum of a week between surveys to maximize cover board effectiveness (Marsh & Goicochea, 2003). Plot I was surveyed from autumn 2009 to spring 2011, and plots II and III were surveyed from autumn 2009 to spring 2013.
We measured two types of responses: behavioral plasticity in the timing, duration, and extent of surface use during the spring and autumn and demography, including the rates of individual growth and population survival. Our general approach allowed us to determine (1) the degree to which behavior and demography responded to environmental variability and (2) whether the two color morphs differed in their response in concordance with past research (Table 1). To gain inference on our two responses, we used three quantitative approaches including mark-recapture (behavior and survival), spatial capture-recapture (movement), and nonlinear growth models (individual growth).

| Behavioral plasticity analyses
Plethodon cinereus in our population are largely underground during the summer and winter due to unfavorable environmental conditions (Taub, 1961). Even during peak surface activity, P. cinereus may be unavailable for capture because they retreat to underground refugia (Bailey, Simons, & Pollock, 2004). We used robust design models that estimate within-season detection probabilities and among-season survival rates while accounting for this temporary unavailability (Kendall, Nichols, & Hines, 1997;Pollock, 1982). Detection probabilities reflect both the probability an individual was on the surface and the probability it was captured and identified. We can therefore estimate when, or under what conditions, salamanders are more likely to be on the surface and how this probability varies within each spring and autumn. Detection probability, parameter p, was modeled using a quadratic function to estimate the optimal environmental condition under which surface use peaked. We tested whether optima were different between morphs for three variables (Table 1): calendar day (prediction 1), the 3-day average rainfall (prediction 2), and the 11-day average of air temperature, which roughly characterizes surface soil temperature T A B L E 1 Predictions generated by climate-morph relationships in the literature

Number
Factor

Model parameter
Behavioral plasticity 1

Surface use and timing
Striped: emerge earlier in spring, peak surface use in early spring, retreat earlier into summer, emerge later in autumn, peak surface use later in autumn, and retreat later in winter. Lead-backed: emerge later in spring, peak surface use later in spring, retreat later into summer, emerge earlier in autumn, peak surface use later in autumn, and retreat earlier into winter Predictions one through four relate to morph differences in behavioral plasticity, and predictions five through eight relate to differences in demography.
Predictions are based on evidence that the striped morph is cool-wet-adapted and the lead-backed morph is warm-dry-adapted. For each prediction, a specific model was developed to test the effect of color morph, and the relevant parameter from that model is specified.
(prediction 3; Kang, Kim, Oh, & Lee, 2000). This resulted in three separate models, one for each environmental predictor. These models included a fixed effect of morph, a fixed effect for morph-environment interactions, and a fixed effect for plot to account for site differences (Table 2). Parameters were estimated using closed-population robust design models in program MARK (White & Burnham, 1999).
See Appendix A for further details on model development.
Another aspect of behavioral plasticity is the breadth and extent of horizontal surface use. Salamanders move to forage, find mates, and defend their territories (Petranka, 1998). We investigated whether morphs moved differently depending on season, with autumn being warmer and drier and spring being cooler and wetter.
Spatial capture-recapture models extend traditional mark-recapture models to better estimate population density by accounting for individual movement (Royle, Chandler, Sollmann, & Gardner, 2014).
These models do so using the location of capture events to estimate a spatial parameter, σ. Morphs that exhibited greater breadth in surface use will have a larger estimated σ. Therefore, we would predict σ to be larger for the lead-backed morph in autumn and larger for striped morphs in the spring (prediction 4; Table 1). For location, we used the coordinate of the cover board under which a salamander was found (Muñoz et al., in press;Sutherland, Muñoz, Miller, & Grant, 2016). We ran the spatial capture-recapture model separately for each season, using program R and package "runjags" to call program JAGS (Denwood, 2016;Plummer, 2003; R Core Team, 2014). For details, see Appendix B.

| Demographic analyses
Closed-population robust design models also estimate apparent survival probabilities (are alive and do not permanently emigrate from study site), Φ, among seasons while accounting for temporary emigration (Kendall et al., 1997;Pollock, 1982). To test predictions relating to demography and environmental conditions (prediction 5 and 6; Table 1), morph-specific survival rates were modeled as a function of season and temperature (Table 2). Prediction 5 predicts that leadback morphs would have higher relative over-summer survival and striped morphs would have higher relative overwinter survival. For prediction 6, we estimate how each morph's seasonal survival relates to mean summer temperature (mean low daily temperature for July and August) and mean winter temperature (mean low temperatures during January and February).
We modified the Faben's (1965) capture-recapture formulation of the von Bertalanffy growth model to estimate individual growth rates of SVL for each color morph (Schofield, Barker, & Taylor, 2013).
Snout-to-vent length is a standard measurement of growth (Leclair, Levasseur, & Leclair, 2006), as salamanders can gain or lose mass rapidly depending on water availability. These models estimate two parameters: a growth coefficient, K, and an asymptotic maximum size, L inf . We allowed growth coefficients to differ by a season by morph interaction (autumn, winter, spring, and summer; Table 2). We would expect the rate of growth in summer and autumn to be higher for lead-backed morphs and rate of growth in winter and spring is to be higher for striped morphs (prediction 7; of the year and to account for the thermal inertia that carries over into the next season. We predicted that lead-backed morphs would grow faster under hotter conditions and that striped morphs would grow faster under cooler conditions. We fit growth models using program R and package "runjags" to call program JAGS (Denwood, 2016;Plummer, 2003;R Core Team, 2014). See Appendix C for model description and JAGS code.

| RESULTS
Parameters are a function of the predictors found within parentheses. All predictors are fixed effects. Parameters not central to predictions found in Table 1 are not included but may be found in the Appendix D. Parameter p is detection probability, Φ is survival, σ is spatial breadth of movement, and K is growth coefficient. Overall, mean striped breadth of movement across all combinations of plot and season was 1.55 m ± 1.24 SD and mean lead-backed was 1.07 m ± 0.514 SD.
We found evidence that environmental conditions may affect survival and growth. Overwinter survival was generally higher than over-summer survival across the three plots ( Figure 4a). There were no significant differences between morphs, but simple comparisons of mean estimates suggest higher overwinter survival by lead-back morph and higher over-summer survival by the striped morph, contradicting prediction 5 (Table 1). We did not find support for a strong effect of temperature on survival probabilities for the population Striped and lead-backed morphs showed the fastest growth during the autumn, followed by less rapid growth in spring. In the winter and summer, growth was severely depressed in both morphs (prediction 7; Figure 5,  Figure 6). Only prediction 8 was supported given that warmer temperatures more negatively impacted striped morphs in the autumn.

| DISCUSSION
Predicting species' responses to climate change require key data on a variety of aspects of an organism's ecology including both demography and behavior (Huey et al., 2012;Urban et al., 2016). Our study shows that the salamander population is clearly influenced by environmental and seasonal conditions both in use of surface habitat and in individual growth. Predicted climate change-warmer temperatures and more variable precipitation (Hayhoe et al., 2007)-will likely be detrimental to P. cinereus populations by shifting the timing and availability of suitable surface conditions. Additionally, we found that warmer temperatures dramatically reduce autumn growth, the most productive season for this species. We also attempted to characterize within-population variation in behavior and demography, given it is a key mechanism for adapting to changing conditions (Barrett & Schluter, 2008). In our study system, we used a color polymorphism as a potential indicator of within-population variation in climate change adaptive capacity.
Multiple lines of evidence suggest that the P. cinereus color polymorphism may be linked to differences in climate niche; however, only one of our eight predictions, temperature-dependent growth, provided moderate support for the climate-morph relationships. While some heterogeneity within the population can be explained by color morph, other ways of characterizing variation in climate tolerance are required.
Behaviorally, warmer temperatures and drier conditions both lead to a reduced presence on the surface for each morph (Figure 2a,b), likely leading to similar patterns in the timing of surface use (Figure 2c).
Although none of the behavioral predictions were supported, contradicting past findings (Anthony et al., 2008;Fisher-Reid et al., 2013;Lotter & Scott, 1977;Moreno, 1989), we did reveal that P. cinereus F I G U R E 1 Plethodon cinereus, the red-backed salamander, is a widely distributed and abundant woodland salamander in eastern North America surface activity is strongly influenced by environmental variables. Two critical aspects of salamander ecology happen on the surface: foraging and courtship (Jaeger, 1980;Petranka, 1998). Our findings suggest that strong seasonal shifts to warmer and drier conditions may limit opportunities for P. cinereus surface activity. As a result, salamanders will need to change the timing of their use to find optimal conditions, increase their reliance on microhabitat refugia, or remain active on the surface under despite likely higher energetic costs (Homyack, Haas, & Hopkins, 2011).
Our seasonal estimates of survival show that mortality was generally greater during the summer than the winter (Figure 4a). In the southern portion of the P. cinereus range where our study takes place, it is likely that desiccation and heat stress in the summer is a greater driver of mortality than cold stress during the winter. Predation, F I G U R E 2 Surface detection as a function of soil temperature (a), rainfall (b), and calendar day (c) for Plethodon cinereus in Laurel, MD, USA. Spring (left) and autumn (right) detection functions are plotted for both morphs. Mean striped morph (solid) and mean lead-backed morph (dashed) estimates are represented by lines. 95% confidence intervals are represented by shaded regions (striped = dark, leadbacked = light). Both temperature and rainfall influence surface detection, leading to bimodal surface activity patterns breeding, and competition also likely contribute to the differences between summer and winter survival. Many salamander predators are in torpor during the winter (e.g., garter snake, Thamnophis sirtalis; Venesky & Anthony, 2007). Plethodon cinereus are also territorial, and antagonistic interactions for desirable microhabitat during the summer may impact demography, as they can result in the loss of a tail (Mathis, 1991;Schieltz, Haywood, & Marsh, 2010). Lastly, breeding may lower summer survival because clutch-laying females may brood their clutch until their energy reserves are depleted (Yurewicz & Wilbur, 2004). Between morphs, we found slight evidence for mean lead-backed summer survival to be lower than mean striped summer survival (Figure 4a). Rather than being climate driven, this is might be because lead-back morphs have poorer quality diets (Anthony et al., 2008), are more submissive to striped morphs (Reiter, Anthony, & Hickerson, 2014), and have poorer quality territories (Paluh, Eddy, Ivanov, Hickerson, & Anthony, 2015). Striped morphs are also more territorial and aggressive, which may prevent lead-backed morphs from finding necessary refugia during the summer (Reiter et al., 2014).
When conditions become stressful during the summer, their survival may be negatively affected the most.
We found large differences in growth rates in relation to season and environmental conditions. Growth is highest in the autumn and is depressed when salamanders are underground in the winter and summer ( Figure 5). Temperature variation during the winter and summer had clear impacts on salamander growth. Most importantly, hot summer temperatures affected growth in the autumn, with warmer summers decreasing growth by 1.5 times compared to our coolest summer observed ( Figure 6). In the autumn, striped morphs responded more negatively to warming temperatures, following prediction 7. During the winter, both morphs responded positively to warmer temperatures. Warmer winter temperatures may increase the chances for opportunistic foraging (Caldwell & Jones, 1973). Contrariwise, hot summer temperatures reduce moisture availability and consequently the leaf-litter invertebrate community (food for salamanders) and may force salamanders underground or to microhabitats where prey and suitable conditions persist (Jaeger, 1972(Jaeger, , 1979. Further, hotter temperatures increase energetic costs, which can slow individual growth (Homyack, Haas, & Hopkins, 2010). Regardless of temperature, striped morphs had higher mean growth rates during surface-active seasons  Second, the environmental conditions we observed may not have been extreme enough to influence demography, or if morph differences did exist, it may not lead to behavioral or demography differences for the factors we measured. Lastly, the P. cinereus color morph may not be a useful indicator in understanding climate tolerances. Our findings lend credence to growing evidence that the polymorphism is not tied to climate (Moore & Ouellet, 2015), but may be maintained by assortative mating (Anthony et al., 2008) or apostatic selection (Fitzpatrick, Shook, & Izally, 2009). Our study suggests that the color morph is an equivocal proxy at best for understanding climate tolerance variability.
One of the most coherent responses to climate change is shifts in species' range distributions (Parmesan, 2006). Dispersal-limited species are less likely to exhibit range shifts and more likely must persist or witness range contractions (Midgley, Hughes, Thuiller, & Rebelo, 2006). Our results on the P. cinereus population may illuminate strategies for how dispersal-limited species at large may persist.
Behaviorally, our results add to the evidence that shifts in phenology, in order to match optimal conditions, are a likely response to climate change (Parmesan, 2006); however, the demographic costs to changes in phenology remain unexplored in many systems (Miller-Rushing, Høye, Inouye, & Post, 2010). Our study suggests that even with possible changes in phenology, increasing summer temperatures will still likely reduce individual growth. Consequently, it could take longer for salamanders to become sexually mature (Nagel, 1977;Sayler, 1966).
Two adaptive strategies arise. First, selection will shift toward behavioral or physiological traits that ensure survival until reproductive size is reached (i.e., demographic buffering hypothesis; Boyce et al., 2006).

F I G U R E 4
Estimates for overwinter ("W") and over-summer ("S") survival probabilities in all three plots ("1", "2", and "3"; a) and survival probability as a function of mean winter temperature (b) for Plethodon cinereus in Laurel, MD, USA. (a) For plots 2 and 3, there were differences between summer and winter survival. Across plots and seasons, there were no clear differences between morphs (striped = black, lead-backed = gray; squares = means, segments = 95% confidence intervals). (b) Striped mean survival (solid, with 95% CI) is not different from lead-backed mean survival (dashed, with 95% CI). Both morphs exhibit relationships not different from zero. There was little summer variation in temperature, so no figure is provided The predicted seasonal growth coefficients from the von Bertalanffy growth analyses for Plethodon cinereus in Laurel, MD, USA. Striped mean growth (solid black) is not different from lead-backed growth (gray) in any of the four seasons ("A" autumn, "W" winter, "Sp" spring, and "Sm" summer). Means are presented as squares with 95% Bayesian credible intervals as segments Our results indicate resource availability will likely be restricted by future suboptimal conditions, so adaptive traits may be those that better secure resources such as suitable microhabitat. For many ectotherms, microhabitat can play an important role in buffering deleterious responses to climate change (Scheffers, Edwards, Diesmos, Williams, & Evans, 2014). However, the minimal variability in behavior in our population suggests physiological traits may become increasingly important. Other systems have shown physiological adaptation to warming conditions (spiders, Krehenwinkel & Tautz, 2013;plankton, Padfield, Yvon-durocher, Buckling, Jennings, & Yvon-durocher, 2015), but given the slower life history of P. cinereus, it is unlikely that adaptation can occur fast enough.
Instead, thermal acclimation may play a central role in how the species mitigate climate-driven restrictions in resource availability (Seebacher, White, & Franklin, 2014; but see Gunderson & Stillman, 2015).
A second adaptive pathway may select life-history strategies that invest in reproduction at smaller sizes and younger ages. Variation in size Growth coefficient K determines the speed at which an individual grows. They were a function of both seasons (autumn, winter, spring, summer), seasonal temperature, and color morph. L inf is the maximum size an individual can reach in Laurel, MD, USA. β represents coefficients from modeling growth as a function of seasonal temperature. Model parameters, parameter description, and the mean estimate (±SE, 95% Bayesian credible interval) are provided.
T A B L E 3 Results from the von Bertalanffy growth models at reproductive maturity already exists in many salamander populations (Peterman, Crawford, & Hocking, 2016;Tilley, 1973); however, "hastening" their life history, while providing more opportunities to reproduce, may reduce overall fecundity as smaller size correlates to fewer eggs per clutch in some salamander populations (Petranka, 1998). Other systems also show similar life-history responses to climate change, including birds (Winkler, Dunn, & McCulloch, 2002), lizards (Bestion et al., 2015), and annual plants (Franks & Weis, 2008). Shifts in life history may be a key response to climate change for dispersal-limited species, but it is unclear whether such a shift is sufficient for populations to persist.
Our goal was to understand how environmental conditions influenced behavior and demography and whether the color polymorphism was useful for understanding within-population variability among those relationships. For organisms like P. cinereus that are dispersallimited, rapid environmental change may overwhelm plastic and adaptive pathways (Chevin, Lande, & Mace, 2010). While changes in populations and distributions are highly idiosyncratic across species (Gibson-Reinemer & Rahel, 2015), our study suggests that projected increases in regional drought and temperature will act as strong negative environmental pressures on P. cinereus population persistence both behaviorally and demographically. Our results also show that the next step is to characterize genetic variability in responses. Genetic variability is a main driver of adaptive capacity (Barrett & Schluter, 2008), and although genomic resources relating genes to phenotypes for many species are undeveloped (Ekblom & Galindo, 2010), P. cinereus is widely studied and will likely have genomic data available in the near future. Our study provides necessary information and insights as to how P. cinereus will be impacted by future climate change (Huey et al., 2012;Urban et al., 2016;Williams, Shoo, Isaac, Hoffmann, & Langham, 2008) and suggests how it, and other dispersal-limited species, may adaptively respond to such impacts.

CONFLICT OF INTEREST
None declared.
F I G U R E 6 The impact of temperature on the growth coefficient for Plethodon cinereus in Laurel, MD, USA. In the autumn, mean growth declines as the previous summer's temperature increases for striped morphs (solid) and lead-backed morphs (dashed). In the winter, increases for both morphs. In the spring, there were no strong relationships. Lastly, summer growth remained constantly low despite warming temperatures. Means are lines with 95% Bayesian credible intervals as shaded regions (striped = dark, lead-backed = light) APPENDIX A

MARK MODELS
Robust design population models use primary seasons (e.g., the autumn and spring seasons in our study) and secondary occasions (e.g., the survey occasions within seasons) to estimate abundance, apparent survival, detection, and temporary emigration (Kendall et al., 1997;Pollock, 1982 show the results from Akaike's information criteria (AIC c ) model selection that determined the final structure of our models. Using the best-supported structure for the nuisance parameters, we then generated five models to test predictions regarding differences in surface use and survival between the color morphs. This was carried out by examining support for an interaction between color morph and our driver of interest. This resulted in five models, each testing a different climate-morph prediction (predictions 1-3, 5, and 6; Table 2).

Predictor in JAGS
The most basic spatial capture-recapture model is referred to as SCR 0 by Royle et al. (2014) . SCR 0 is a single season closed-population model for estimating density, and it can be used to estimate four parameters: abundance, density, detection probability, and breadth of movement.
To estimate abundance, a homogenous binomial point process is used,  The best model was interactive between plots and secondary sampling occasions. This allows us to test our predictions by modeling detection using secondary sampling occasion covariates like rainfall, temperature, and calendar day.
The probability an individual will be captured at a given trap declines with distance from the animal's activity center, often specified using a half-normal encounter probability, where the probability of encountering an individual at location x with activity center s is a function of p 0 the baseline detection probability, σ the breadth of the detection kernel, and ||x j − s i || the Euclidean distance between the location of trap j and the activity center of individual i. The spatial encounter histories (the data), y ijk , are then evaluated. The probability that individual i is caught (y ijk = 1) or not caught (y ijk = 0) in trap j on occasion k follows a Bernoulli trial where y ijk ~ Bernoulli(p ijk ), and p ijk comes from equation A1. In SCR 0 , detection is constant across the study period, so it is possible to "flatten" the encounter history so that captures (y ij = 1, 2, 3, …, n) or noncaptures where α 0 is our baseline detection rate and α 1 represents the coefficient for how fast detection decreases with distance. If we wanted to model detection or space use as a function of biological or environmental data, it is possible to estimate the effects (β 1 , …, β t ) of t covariates (v 1 , …, v t ) on detection.
As with other closed-population abundance models, model SCR 0 assumes demographic closure and some degree of geographic closure.
However, SCR relaxes the assumptions of equal detectability among individuals. For instance, it is not necessary to ensure all individuals have a trap within their home range ("no holes"). Activity centers are assumed to be randomly distributed throughout the state space and are independent of each other. Given that this is rarely true in real animal populations, SCR models are fortunately robust to violations of the uniform distribution assumption. Third, detection is assumed to decline as distance increases from an animal's activity center. Lastly, encounters are assumed to be independent, meaning animals do not exclude one another and that individuals do not become adverse to traps. where i is the interval between marking and recapture, L r i is the estimated size at recapture, L m i is the size when an individual was marked, (A1) p(x,s) = p 0 exp − 1 2σ 2 ||x j − s i || , L inf is the maximum size an individual can reach, K i is the growth coefficient, and Δi is the duration of the interval in days (scaled to year by dividing by 365). Only individuals captured more than once were used in this analysis, and SVL measurements were used for size.
To test the predictions regarding growth, we first modeled the growth coefficient as a function of four seasons to address prediction seven.
where each of the growth coefficients represents the autumn, winter, spring, and summer, respectively. For each growth coefficient, the number of autumn season days, Δf i , winter days, Δw i , spring days, Δsp i , and summer days, Δsm i , were used for each capture interval. Seasons were defined the same across all 4 years, where the spring (March 2-May 16) and autumn (September 6-December 4) always contained all field surveys, and summer and winter were the periods between these surveys.
To test prediction eight, we modeled the growth as a function of mean season temperature. To evaluate the impacts of extreme heat and extreme cold, we modeled the surface-active seasons (autumn and spring) as a function of previous summer and previous winter temperatures. Because the growth coefficient must always be positive to avoid estimation errors, the log of the growth coefficients was modeled as, where β 0 is the mean growth rate in each season, β temp is the coefficient for how growth changes with temperature, and T i is the mean of the season's low temperature in interval i. Temperature effects were ran separately for each of the four seasons. The temperature coefficients will determine whether morphs exhibit differential demography in regard to heat or cool stress.
This model assumes that the growth models start from age zero and that growth rates are conditional on the maximum asymptotic size an animal can reach. We used vague priors for all parameters and kept