Demography and space‐use of Eastern Red‐backed Salamanders (Plethodon cinereus) between mature and successional forests

Abstract Space‐use and demographic processes are critical to the persistence of populations across space and time. Despite their importance, estimates of these processes are often derived from a limited number of populations spanning broad habitat or environmental gradients. With increasing appreciation of the role fine‐scale environmental variation in microgeographic adaptation, there is a need and value to assessing within‐site variation in space‐use and demographic patterns. In this study, we analyze 3 years of spatial capture–recapture data on the Eastern Red‐backed Salamander collected from a mixed‐use deciduous forest site in central Ohio, USA. Study plots were situated in both a mature forest stand and successional forest stand separated by <100‐m distance. Our results showed that salamander density was reduced on successional plots, which corresponded with greater distance between nearest neighbors, less overlap in core use areas, greater space‐use, and greater shifts in activity centers when compared to salamanders occupying the mature habitat. By contrast, individual growth rates of salamanders occupying the successional forest were significantly greater than salamanders in the mature forest. These estimates result in successional plot salamanders reaching maturity more than 1 year earlier than salamanders on the mature forest plots and increasing their estimated lifetime fecundity by as much as 43%. The patterns we observed in space‐use and individual growth are likely the result of density‐dependent processes, potentially reflecting differences in resource availability or quality. Our study highlights how fine‐scale, within‐site variation can shape population demographics. As research into the demographic and population consequences of climate change and habitat loss and alteration continue, future research should take care to acknowledge the role that fine‐scale variation may play, especially for abiotically sensitive organisms with limited vagility.

3 years of spatial capture-recapture data on the Eastern Red-backed Salamander collected from a mixed-use deciduous forest site in central Ohio, USA. Study plots were situated in both a mature forest stand and successional forest stand separated by <100-m distance. Our results showed that salamander density was reduced on successional plots, which corresponded with greater distance between nearest neighbors, less overlap in core use areas, greater space-use, and greater shifts in activity centers when compared to salamanders occupying the mature habitat. By contrast, individual growth rates of salamanders occupying the successional forest were significantly greater than salamanders in the mature forest. These estimates result in successional plot salamanders reaching maturity more than 1 year earlier than salamanders on the mature forest plots and increasing their estimated lifetime fecundity by as much as 43%. The patterns we observed in space-use and individual growth are likely the result of density-dependent processes, potentially reflecting differences in resource availability or quality. Our study highlights how fine-scale, within-site variation can shape population demographics. As research into the demographic and population consequences of climate change and habitat loss and alteration continue, future research should take care to acknowledge the role that fine-scale variation may play, especially for abiotically sensitive organisms with limited vagility.

| INTRODUC TI ON
Landscapes are composed of biotic and abiotic features that are heterogeneous in at least one dimension at any scale (Turner & Gardner, 2015). Extant populations of organisms are usually acclimated to their local landscapes (Urban et al., 2016); however, much of our knowledge about a given species is often derived from a limited number of populations and assume demographic processes are spatially invariant. Such generalizations may be particularly misleading for organisms that interact with their environment at fine scales (i.e., tens of meters or less), such as plants, invertebrates, and many amphibians because these organisms typically have limited dispersal and may have physiological and/or behavioral constraints that restrict their ability to actively select habitat (De Bie et al., 2012). It is therefore imperative to evaluate variation in demographic processes at spatial scales relevant to the target organisms to mitigate potential biases of broad-scale generalizations.
An increasing body of evidence indicates that microgeographic adaptation occurs among numerous species and ecosystems (Richardson et al., 2014), influencing distribution, abundance, and individual phenotype (Cicchino et al., 2021). For example, limpets (genus Patella) separated by <2 m experienced significantly different sun exposure and thermal stress depending on the side of the rock that they inhabit (Seabra et al., 2011). Chronic thermal stress can significantly impact growth, reproduction, and overall fitness (Dantzer et al., 2014;Wingfield & Romero, 2011), which likely results in heterogeneous fitness among individuals in close proximity to each other. Similarly, high-elevation hatchlings of the water snake, Natrix maura, are smaller in size and poorer swimmers compared to those at low elevations, likely as a result of hypoxia and lower temperatures (Souchet et al., 2020). Therefore, the evaluation of finescale demographic differences between populations is a key step to understanding spatial nuance that may be important for the development of comprehensive conservation and management goals of organisms with limited habitat selection capability. Space-use, or the amount and extent of a given area used, by individuals within a population shapes demographic patterns and can be influenced by the availability of resources, appropriate microclimate, and intraspecific density (Gaillard et al., 2010;Morales et al., 2010). Population density may be particularly influential because of the many cues it signals. For example, high-density areas may indicate good quality habitats that can support more individuals and may offer cooperation with conspecifics (e.g., antipredator grouping behavior, resource sharing, and cooperative breeding) and consequently promote site philopatry and reduced movement and space-use (Clobert et al., 2009;Le Galliard et al., 2005). Conversely, high-density areas may have higher competition for resources, mates, and more aggressive individuals thus promoting space-use and movement away from the high-density site (Clobert et al., 2009).
Similarly, population demographic parameters such as individual growth rates can vary in density-dependent ways whereby higher densities result in lower individual growth rates, due to fewer resources available to each individual, and allocation of obtained resources toward other processes like aggressive interactions and competition (Getz, 1996). Animal movement and space-use vary across species geographic range (Boyle et al., 2009), but also between populations in close proximity especially when microhabitats differ (Gonzales et al., 2020;Reeve et al., 2009). Yet, measuring life history and space-use at fine scales can be challenging due to cost, labor intensity, time, and methodological constraints. Thus, we lack a comprehensive understanding of fine-scale population dynamics.
Plethodontid salamanders are a particularly well-suited group to evaluate fine-scale variation in key demographic patterns due to their high abundances across a wide geographic range (Petranka, 1998), generally low dispersal rates, and relative ease of repeated captures. Terrestrial lungless salamanders in the genus Plethodon are the most abundant vertebrate animals in many North American forests, accounting for more vertebrate biomass than any other taxa in these ecosystems (Burton & Likens, 1975;Semlitsch et al., 2014).
Plethodon are highly philopatric and rely on cool and moist microhabitat conditions to facilitate cutaneous respiration. The Eastern Red-backed Salamander, Plethodon cinereus, is one of the most abundant and widely distributed species, and while we have some understanding of demographic patterns across their geographic range (Anthony & Pfingsten, 2013;Nagel, 1977;Sayler, 1966), the fine-scale demographic differences between populations separated by small distances are largely unknown. Plethodon cinereus is found across eastern North American forests but has the highest density in forests with greater percent canopy cover, larger trees, and with high densities of well-decayed coarse woody debris (McKenny et al., 2006;Otto et al., 2014;Wilk et al., 2020). Forests with these attributes provide suitable cool and moist microhabitats and higher prey abundance for salamanders. However, habitat heterogeneity can exist across fine scales with variable conditions and resources within a small area. As such, salamander demographic patterns likely differ between local populations (Farallo & Miles, 2016) and targeted research to unveil these differences is necessary.
Here, we use a multiseason spatial capture-mark-recapture study to evaluate fine-scale demographic rates of P. cinereus inhabiting sites with different forest compositions but separated by just 100 m. We predicted that salamanders occupying the two different forest types will exhibit differences in density, space-use, and individual growth K E Y W O R D S demography, fecundity, fine-scale, forest, salamander, succession

T A X O N O M Y C L A S S I F I C A T I O N
Applied ecology, Autecology, Movement ecology, Population ecology, Spatial ecology, Zoology rates. Specifically, we expected that early-successional forest habitat would be suboptimal to mature forest due to reduced stand density promoting warmer and drier microclimate and reduced aboveground biomass that may reduce invertebrate prey (McKenny et al., 2006). As such, successional stands were predicted to support fewer individual salamanders with larger home ranges and lower individual growth rates compared with salamanders occupying the mature forest habitat.

| Field sampling
We conducted this study in Galena, Ohio, at a 36-ha site consisting of a mix of mature oak-hickory forest (Quercus and Carya spp.), early-to mid-successional mesic hardwoods (Acer spp.), white pine plantations (Pinus strobus), and open field habitats ( Figure 1, Table 1).
The early-successional forests are <40 years in age and are growing in what was previously pastureland used for low-intensity grazing in the 1970s. The mature forests are found on the ravine slopes and bottomlands surrounding a rocky stream that flows through the property. The well-drained, relatively undisturbed upper slopes of this ravine are adjacent to the flat, poorly drained, early-successional upland forests that were historically grazed.
We established four cover board arrays each consisting of wooden boards measuring 30 × 30 × 2.5 cm. We set two arrays in successional habitat and two arrays on the adjacent mature forest habitat. Each pair of arrays was at least 20 m apart and consisted of a 5 × 10-m array of 50 cover boards equally spaced 1 m apart.
Successional and mature forest arrays were 60-100 m apart. This design minimized distance between arrays while capturing perceived differences in abiotic gradients of temperature and moisture. We established all arrays during Fall 2016 and began sampling for salamanders in Spring 2017 after one season of board weathering. We identified all tree species and measured the diameter at breast height (DBH) in a 10-m radius around the center of each array ( Table 1).
Plethodon cinereus in Ohio are most active at the surface during spring and autumn and retreat into deeper soil to avoid desiccation and freezing during summer and winter, respectively (Anthony & Pfingsten, 2013). Thus, we sampled for salamanders three times provides safe and long-lasting marks (Grant, 2008;Northwest Marine Technologies, LLC). Each recaptured individual was identified by reading the visual elastomer implants with a UV flashlight. We also recorded snout-vent length (SVL, mm), tail length (TL), and sex of every individual. Males were determined based on the presence of a mental gland or visual testes, and females were determined by the presence of eggs (Gillette & Peterson, 2001). If neither were visible, we made an expert determination based on the shape of the nose as being either square (males) or round (females) (Dawley, 1992). Individuals were considered juveniles if they were <34 mm SVL (Anthony & Pfingsten, 2013) and were not assigned a sex. We returned salamanders to their board of capture within 4 h of initial capture.
Sampling occurred between 09:00 and 12:00 and during each sampling occasion, we measured weather covariates at each array including soil surface temperature, average leaf litter depth, and temperature. Additionally, we collected soil core samples at the center of each array to determine organic soil depth and obtain soil water content by measuring the difference in mass before and after airdrying soil samples. We used t-tests to determine statistical differences in these survey-level measures between the successional and mature sites. We also retrieved weather covariates for each sampling survey for our study site from the PRISM dataset. For each year, we designated spring and fall as the active seasons, with the remainder of the year as the inactive seasons. We estimated average temperature, precipitation, and days since rain using a 5-day moving window analysis for every day of active and nonactive seasons. We calculated the coefficient of variation for 5-day average temperature and total precipitation by dividing the seasonal mean by the standard deviation.

| SCR model
We investigated survival and space-use parameters using a robustdesign spatial capture-recapture (SCR) adapted from Ergon and Gardner (2014). The robust design describes a sampling structure that divides "primary" sessions and "secondary" sessions within each primary session. The robust design assumes that demographic processes are closed between secondary sessions but open between primary sessions (Pollock, 1982). In our study, fall and spring sampling seasons serve as the primary sessions with 2-3 secondary sessions within each primary session. An SCR differs from a traditional Cormack-Jolly-Seber capture-recapture model by explicitly incorporating spatial capture locations to account for individual movement or dispersal, allowing for a more accurate estimate of true survival (Schaub & Royle, 2014). Dispersal distance is an estimate of the difference between activity centers between seasons and activity centers were assumed to have a uniform distribution and dispersal only occurred between primary sessions. We included the aforementioned PRISM weather covariates in the survival submodel of the SCR; however, null models were better supported and we subsequently only report results from those models.

| Individual growth model
We estimated individual growth using Fabens capture-recapture growth model (Fabens, 1965). The growth function for individual i at time t was defined as: where asymptotic size L was allowed to differ by sex and was estimated from a normally distributed prior with a mean of 48 and precision of 0.01, a prior based on estimates from . SVL0 i,t is the size at first capture and follows a uniform distribution with a minimum of 10 and maximum of 60. We removed any observations where SVL at the final time step was less than SVL at the first time step, which would be due to measurement errors.
K represents the individual growth rate, and I is the annual scaling

| Space-use
Using parameters estimated from our fitted SCR model, we assessed space-use and overlap in salamanders occupying the successional and mature habitats. Specifically, we plotted each individual's spatial location in coordinate space and then calculated the probability (p) of each individual (i) using adjacent spatial locations (j) as a function of distance (d) following a negative exponential function: The rate of probability decay in space is governed by σ, which was estimated during the fitting of the SCR model. Probability-ofuse surfaces were created for each individual at 1000 samples (k) of the fitted posterior model distributions. Using the probability surfaces p ijk , we distributed 1000 hypothetical "use" points on the landscape following a random multinomial process. We then calculated kernel density utilization distributions (UD) of these spatially referenced use points for each individual and posterior sample using the R package "adehabitatHR" (Calenge, 2006). Finally, we calculated the probability that the core 50% of each individual i's UD overlapped with all other core 50% UD j calculated as the probability of home range overlap (Fieberg & Kochanny, 2005). We then determined the average number of individuals with overlapping core UDs, as well as the average probability of overlap.

| Population projection
Using parameters from our fitted individual growth model and from the literature, we conducted population projection simulations to understand how differences in individual growth could impact lifetime fecundity. We assume that all individuals are 13.5 mm SVL upon hatching and that sexual maturity is first reached at 34 mm SVL (Anthony & Pfingsten, 2013). However, following Lotter (1978), we assume that individuals between 34 and 43 mm SVL have 56% chance of being gravid, while 94% of females >43 mm SVL are likely to be gravid (Jaworski et al., 2018;Wise & Jaeger, 2021). Regional variation in reproduction has been documented (Lotter, 1978;Sayler, 1966;Werner, 1971), but our data are not sufficient to ascertain frequency of reproduction in our Ohio population and we therefore use the averages reported by Lotter (1978). Similarly, we could not confidently determine the average number of eggs produced by each female, nor whether there was a size by fecundity interaction.
As such, we fit a linear model with a normal distribution to the data reported in Lotter (1978) relating clutch size to SVL using the R package "brms" (Bürkner, 2016). We found that the normally distributed

| RE SULTS
Tree species diversity was greater at successional sites (n = 9 species) relative to mature sites (n = 5), and both sites were dominated by maple species. There was no difference in average DBH between successional and mature sites. There were fewer trees at successional sites (n = 57) than at mature sites (n = 76), but the greater number of trees present at mature sites was largely driven by maple saplings ( ( Table 2). Overall, observations of salamanders co-occurring under the same cover board were more than twice as common on mature plots (297) compared with successional plots (135). This equates to 76% of all salamander observations under mature plot cover boards being with another salamander, while only 46% of salamanders were with another salamander under successional plot cover boards. The most frequent co-occurrence observation among mature salamanders at both mature and successional forest plots was female-female; male-male and male-female co-occurrences were 2-3 times less common ( Table 2).

| SCR model
For most parameters estimated in our spatial capture-recapture model, mature and successional forest plots had moderate differences (Table 3). Annual survival was estimated to be >0.99 for both successional and mature forest plots. Activity centers of salamander in successional plots shifted slightly more between primary sample periods when compared to mature forest plots (1.376 m vs. 1.241 m, respectively), and successional plot salamanders exhibited greater space-use (successional = 3.823 m, mature = 3.496 m). The density of salamanders was significantly higher in mature forest plots, which also had significantly higher probability of detection ( Table 3).

| Individual growth estimates and time to maturity
On average, initial salamander mean SVL ± standard deviation was Male-female 32 10 Female-female 63 23 Note: Habitat values are means (±SD) of measurements collected during each survey (n = 17 measurements). Co-occurrence summaries report the number of times, across all surveys and recaptures, that two or more salamanders were found under the same cover board. Co-occurrence observations could be male-male, female-female, or male-female; observations with at least one juvenile salamander are included in "All."

TA B L E 3
Parameter estimates from fitted spatial capture-recapture model for Plethodon cinereus in central Ohio. no observed difference in the overall mean SVL between mature forest (37.61 ± 4.45 mm) and successional forest (36.43 ± 5.16 mm) plots.
However, there were significant differences in asymptotic growth and individual growth rates between males and females and significant differences in individual growth rates between the mature and successional forest locations (Table 4; Figure 2). This results in males in the successional forest reaching sexual maturity in 2.25 years and males in the mature forest maturing in 2.75 years, while females in the successional forest sexually mature in 3.30 years and those in the mature forest sexually reach maturity in 4.30 years (Figure 3).
Despite the greater number of individuals potentially occupying the same space in the mature forest habitat, the average probability of overlap was nearly identical between the two habitats (Table 5).
However, the average distance to the nearest neighboring salamander tended to be less for salamanders occupying mature habitat (0.40 m ± 0.06) compared with salamanders occupying successional habitat (0.50 m ± 0.08; Figure 4).

| Population projection
Given the annual survival rate estimated from , females are estimated to live an average of 5.87 years (± 4.90). Females in the mature habitat are estimated to average 2.2 (±3.7) clutches in their lifetime, equating to a projected mean lifetime fecundity of 13.7 (±27.3). Because maturity is reached earlier in the successional habitat, females are estimated to average 2.8 (±4.1) clutches and produce a mean of 19.6 (±33.5) offspring in their lifetime, which is 43% more than females occupying mature forest (Table 5).

| DISCUSS ION
Demographic processes, population vital rates, and space-use are driven by the abiotic and biotic environment experienced by an organism. While variation in these rates is often expected across broad F I G U R E 2 Density plot showing the posterior distributions for asymptotic size for male and female Eastern Redbacked Salamanders, Plethodon cinereus (a) and posterior distributions for the growth coefficient for males and females occupying mature and successional forest plots (b). With greater than 99% probability, all contrasts indicate that males are smaller than females, males grow faster than females, and males and females on the successional forest plots grow faster than males and females on the mature forest plots. spatial scales (e.g., latitude), our study shows that variation can exist at fine spatial scales between animals occupying different microhabitats and separated by as little as 100 m. In this study, we predicted that there would be differences in salamanders occupying the mature and successional forest plots due to environmental differences between sites. As expected, salamanders in successional forest habitat were present in lower densities and had greater space-use and activity center shifts than salamanders from mature forest habitat.
Salamanders in successional habitat were estimated to be farther from their nearest neighbor and to have less core use overlap with conspecifics than salamanders in mature forest plots. Contrary to our predictions, however, these differences corresponded with salamanders in successional forest growing more rapidly, reaching sexual maturity sooner, and, based on our population projection, having greater projected lifetime fecundity. Ultimately, the processes shaping these patterns remain uncertain as we observed no meaningful differences in soil moisture or air and soil temperature. There were pronounced differences in tree composition and size between sites, but limited differences in leaf litter depth.
We found a higher density of animals with subsequently shorter distances between individuals and greater home range overlap on mature forest plots. The higher proportion of adults on the mature plots, the higher recapture rates on mature plots, and the observed space-use patterns on mature plots relative to the successional plots are suggestive of territorial adults defending their F I G U R E 3 Time to maturity plot, indicating the expected time it would take female (a) and male (b) Eastern Red-backed Salamanders, Plethodon cinereus, to reach sexual maturity (34 mm SVL), given their development in either mature or successional habitats. Starting from a hatching SVL of 13 mm, 50% of juvenile female salamanders occupying successional habitat are expected to reach the minimum size of sexual maturity after 3.30 years compared to 4.30 years for females occupying the mature forest plots (a). By contrast, 50% of males are expected to reach maturity in 2.25 years and 2.75 years in successional and mature forest habitat, respectively (b). Note: Probability of home range overlap (PHR) is reported as both the average for each individual (i) relative to all other individuals (j) within the same survey plot, as well as the maximum probability of core utilization distribution (UD) overlap. Overlap reports the average number of individuals with overlapping core UD and NN summarizes the average distance to the next closest salamander in the plot. All statistics are means (±SD). home range. It is also evident that salamanders encounter each other more frequently at the mature forest plots, but they may not engage in agonistic behaviors due to the energetic costs related to frequent aggressive interactions and allow such overlap to occur.

TA B L E 5
Plethodon cinereus and other terrestrial plethodontids are known to reduce agonistic interactions with familiar conspecifics ("dear enemy hypothesis"), especially in areas with high density (Dalton et al., 2020;Jaeger, 1981;Jaeger & Peterson, 2002). The majority of sexually mature salamander co-occurrence observations in our study were between the same sex and were more prevalent on the higher density mature plots ( Table 2). Our home range analyses, in addition to co-occurrence summaries, suggest differential behaviors of space and cover use at fine scales. In contrast to Hernández-Pacheco et al. (2019) who found that home ranges are not limited by density, we observed that space-use was reduced at the higher density mature forest plots. The forest where our mature plots are situated has remained largely undisturbed relative to other central Ohio sites and has the highest observed density of salamanders per square meter when compared to nine other sites in central Ohio, USA (Wilk et al., 2020). Historically, agricultural land uses were more widespread throughout Ohio and eastern North America, with forest cover often increasing with agricultural abandonment (Drummond & Loveland, 2010;Monsted & Matlack, 2021). Our successional plots were historically used as sheep pasture and likely were not widely occupied by P. cinereus.
As such, mature vs. successional plot differences in our study may reflect habitat quality as well as historical population stability and may be a microcosm of broader trends of forest reversion across Eastern North America and of P. cinereus.
The most prominent effect observed in our study was the significant difference in individual growth rates between our plots ( Figure 3). Salamanders occupying successional plots grew faster and reached maturity >1 year earlier than individuals occupying mature forest plots. There are at least two possible mechanisms for the observed differences. The first may be related to tree community differences between successional and mature sites ( Table 1). Forests across Eastern North America have undergone mesophication and shifted from oak-dominated tree communities to maple-dominated (McEwan et al., 2011;Nowacki & Abrams, 2008). Maple leaf litter is often considered higher quality because of the greater proportion of nitrogen to carbon and greater microbial diversity, which is partially responsible for the faster decomposition rates relative to oak (Laking et al., 2021;Lehmann et al., 2020). Such microbial abundance may support a more diverse quality and quantity of prey items available to salamanders (Rittenhouse et al., 2008;Templer et al., 2003) resulting in the faster individual growth and lifetime fecundity observed on the successional plots. Anurans have also been observed to have higher growth and survival when raised in maple-dominated mesocosms compared with oak (Breslau, 2018). Maple-dominated forests also tend to be wetter and cooler (McEwan et al., 2011), conditions amenable to lungless Plethodon salamanders. While we did not measure notable differences in moisture between our plots, our successional sites had larger maple trees and only one oak, a stark contrast to the 16 large oaks found at the mature sites (Table 1).
A second possible driver, especially in individual growth rate, could be density dependence. Harper and Semlitsch (2007) found that density had a negative effect on survival and growth in metamorphosed American toads (Anaxyrus americanus) and wood frogs (Rana sylvatica), and Berven (2009) reaffirmed these effects in a long-term data set of wood frogs. Numerous other studies have identified density dependence in demographic parameters in larval or aquatic urodeles (e.g., Bendik & Dries, 2018;Ousterhout & Semlitsch, 2016;Semlitsch, 1987;Van Buskirk & Smith, 1991), but there is limited research into how density directly affects population demographic parameters of terrestrial plethodontid salamanders. It is important to note, however, that we do not have any estimates of food availability or quality, which should be a focus of future work to better understand the role of density-dependent processes (Hantak et al., 2016;Kuzmin, 1995).
Regardless of mechanism, the differences in growth rates substantially increase projected lifetime fecundity for females occupying the successional forest plots (Table 5). It is possible there is a greater rate of emigration off successional plots, which may F I G U R E 4 Density plot showing the average number of individual Plethodon cinereus that had overlapping core (50% UD) home ranges with each salamander (a) and the average distance between salamander activity centers (b) in mature and successional forest plots. Salamanders occupying mature habitat tended to have more individuals potentially occurring within their core home range than salamanders occupying successional forest habitat, which coincided with less distance between salamanders occupying mature forest plots.
be reflected in our data as we had ~6% lower recapture rate on successional plots as compared to mature forest plots. Spatial capture-recapture models fit to data collected under a robust design allow for the estimation of true rather than apparent survival (Ergon & Gardner, 2014;Gardner et al., 2010;, but permanent emigration remains an elusive parameter. Emigration in plethodontid salamanders can be particularly challenging as salamanders can temporarily migrate underground or can disperse over land to a new location. During any given survey, only a small fraction of the population is available to be sampled on the surface (Bailey et al., 2004a(Bailey et al., , 2004b. More research is needed to better understand emigration, especially in seasonally active animals such as P. cinereus.

| CON CLUS ION
Demography is the most critical driver affecting population persistence (Hanski & Gilpin, 1991). Low reproductive rates, slow maturation, and longer generation times all increase the susceptibility of a population to stochastic events and the potential for local extinction (McKinney, 1997). However, variation in life history traits can buffer populations when environments change (Anderson et al., 2015).
Plethodon cinereus has proven to be a resilient and adaptable species with a distribution encompassing much of eastern North America and populations frequently persisting in highly altered or urbanized landscapes (Gibbs, 1998;Petranka, 1998;Wilk et al., 2020). The ability to thrive and not just persist in altered or changing habitats may be critical to the species' broad distribution and persistence.
We found P. cinereus occupying successional habitat to have greater growth rates, which are predicted to result in earlier maturation and greater lifetime fecundity. While observed differences in salamander growth rates and the subsequent demographic differences may be driven by habitat variation and density-dependent processes, it is also possible for there to be different selective pressures between habitats, leading to microgeorgraphic adaptation. Our results reinforce the role that fine-scale variation can play in spatial-temporal population processes. Perhaps most notably, these differences occurred between sites <100 m apart, highlighting the importance of accounting for fine-scale, within-site variation when assessing demographic processes. As research into the demographic and population consequences of climate change and habitat loss and alteration continue, future research should take care to acknowledge the role that fine-scale variation in both biotic and abiotic environments may play, especially for organisms with small home ranges or limited mobility.

ACK N OWLED G M ENTS
We would like to thank the many people who aided in the field data collection including Eric Gangloff, Alyssa Baxter, Christy Kunkle, Kasey Foley, Amelia Gray, Renna Wittum, and Aaron D'Amore.
We also thank Harlan Bockbrader, the Overbrook Presbyterian

CO N FLI C T O F I NTE R E S T
The authors do not have any competing interest.