Functional and numerical responses of shrews to competition vary with mouse density

PLOS ONE, Jan 2018

For decades, ecologists have debated the importance of biotic interactions (e.g., competition) and abiotic factors in regulating populations. Competition can influence patterns of distribution, abundance, and resource use in many systems but remains difficult to measure. We quantified competition between two sympatric small mammals, Keen’s mice (Peromyscus keeni) and dusky shrews (Sorex monticolus), in four habitat types on Prince of Wales Island in Southeast Alaska. We related shrew density to that of mice using standardized regression models while accounting for habitat variables in each year from 2010–2012, during which mice populations peaked (2011) and then crashed (2012). Additionally, we measured dietary overlap and segregation using stable isotope analysis and kernel utilization densities and estimated the change in whole community energy consumption among years. We observed an increase in densities of dusky shrews after mice populations crashed in 2012 as expected under competitive release. In addition, competition coefficients revealed that the influence of Keen’s mice was dependent on their density. Also in 2012, shrew diets shifted, indicating that they were able to exploit resources previously used by mice. Nonetheless, increases in shrew numbers only partially compensated for the community energy consumption because, as insectivores, they are unlikely to utilize all food types consumed by their competitors. In pre-commercially thinned stands, which exhibit higher diversity of resources compared to other habitat types, shrew populations were less affected by changes in mice densities. These spatially and temporally variable interactions between unlikely competitors, observed in a relatively simple, high-latitude island ecosystem, highlight the difficulty in assessing the role of biotic factors in structuring communities.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0189471&type=printable

Functional and numerical responses of shrews to competition vary with mouse density

January Functional and numerical responses of shrews to competition vary with mouse density Carolyn A. Eckrich 1 2 3 4 Elizabeth A. Flaherty 1 2 4 Merav Ben-David 0 1 2 4 0 Department of Zoology and Physiology and Program in Ecology, University of Wyoming , Laramie, WY , United States of America 1 Funding: Funding for the project was provided by the USDA Forest Service Tongass National Forest under cooperative agreement 10JV-11100500-021 ( 2 Data Availability Statement: Data included in this manuscript is available at the University of Wyoming Data Archive at the following DOI: https://doi.org/10.15786/M24T00 , USA 3 Oregon Department of Fish and Wildlife, La Grande, OR, United States of America, 2 Department of Forestry and Natural Resources, Purdue University , West Lafayette, IN , United States of America 4 Editor: Aaron W. Reed, University of Missouri Kansas City , UNITED STATES For decades, ecologists have debated the importance of biotic interactions (e.g., competition) and abiotic factors in regulating populations. Competition can influence patterns of distribution, abundance, and resource use in many systems but remains difficult to measure. We quantified competition between two sympatric small mammals, Keen's mice (Peromyscus keeni) and dusky shrews (Sorex monticolus), in four habitat types on Prince of Wales Island in Southeast Alaska. We related shrew density to that of mice using standardized regression models while accounting for habitat variables in each year from 2010±2012, during which mice populations peaked (2011) and then crashed (2012). Additionally, we measured dietary overlap and segregation using stable isotope analysis and kernel utilization densities and estimated the change in whole community energy consumption among years. We observed an increase in densities of dusky shrews after mice populations crashed in 2012 as expected under competitive release. In addition, competition coefficients revealed that the influence of Keen's mice was dependent on their density. Also in 2012, shrew diets shifted, indicating that they were able to exploit resources previously used by mice. Nonetheless, increases in shrew numbers only partially compensated for the community energy consumption because, as insectivores, they are unlikely to utilize all food types consumed by their competitors. In pre-commercially thinned stands, which exhibit higher diversity of resources compared to other habitat types, shrew populations were less affected by changes in mice densities. These spatially and temporally variable interactions between unlikely competitors, observed in a relatively simple, high-latitude island ecosystem, highlight the difficulty in assessing the role of biotic factors in structuring communities. Introduction The relative importance of biotic interactions and abiotic factors in regulating population dynamics and community structure has been at the core of ecological investigations for over a design, data collection and analysis, decision to publish, or preparation of the manuscript. century [ 1,2 ]. Ideas such as divergence of sympatric species [ 3,4 ], fundamental and realized niches [ 5 ], the neutral theory of biodiversity and biogeography [ 6 ], and compensatory ecological (or zero-sum) dynamics [ 7,8 ], all derive from a vast body of literature documenting the effects of biotic interactions (largely competition) on community composition, diversity, and species abundance. Nonetheless, recent re-analyses of published works suggest that in a majority of cases evaluated, abiotic factors may override the effects of biotic interactions [9]. Theory suggests that in cases with substantial overlap of fundamental niches (i.e., the set of resources potentially used by a species) and asymmetric competition, the superior competitor may drive the other species to local extinction [ 10,11 ]. However, when niche overlap is partial, the weaker competitor may be restricted to a portion of the space, which will represent its realized niche [5]. To allow coexistence of potential competitors, the realized niches of sympatric species may separate in three dimensions: time, space, and food [ 12,13 ]. For example, in Arizona, Brown [14] found that the coexistence of desert granivorous rodents was best explained by seasonal rotation in foraging strategies, where each species maximized its efficiency during a different part of the year. In such instances, competitive release could result in expansion of the realized niche of the weaker competitor [ 15 ]. Numerous empirical studies have shown that competitive interactions among species influence spatial and temporal patterns of distribution, abundance, and resource use [ 16,17 ]. Yet interspecific competition is difficult to measure in the field where multiple processes can interact to mask its effects. Most prominent is the impact of predation, which may produce apparent competition through indirect effects [18] or may conceal competitive interactions of prey populations [ 19 ]. In a meta-analysis, Gurevitch et al. [ 20 ] concluded that the presence of predators lessened the effect of competition. Further, competition may be obscured by unrelated evolutionary pressures. For example, Ben-David et al. [ 21 ] demonstrated that coastal river otters (Lontra canadensis) and American mink (Neovison vison) spatially segregate and coexist in the marine environment largely because of limitations on diving abilities of mink, potentially concealing past or ongoing competition [ 22 ]. Further, documenting biotic interactions on islands can be complicated by island size [ 23 ]. Thus, the role of interspecific competition in observed spatial and temporal partitioning is not always clear [ 24,25 ]. Small mammals (<1kg) represent a widespread and diverse assemblage of animals that are particularly suited for experimental study of competition [ 26,27 ]. For instance, removal and exclusion of the larger D. heermanni in California increased the numbers of the federally endangered D. nitratoides nitratoides by 500% within 12 months [28]. However, even experimental studies can fail to demonstrate the effects of competition when abiotic factors overshadow subtle biotic interactions [ 29 ]. As an alternative, negative correlations between the abundance of sympatric species (i.e., negative covariances) have been used to infer interspecific competition [ 13,16 ]. In such analyses, two sympatric species with similar ecological niches will likely exhibit higher competition coefficients (α) than a comparable pair exhibiting niche separation [30]. Luo et al. [ 31 ] found that when abundance was low and the velvet-furred rat (Rattus lutreolus velutinus) and the long-tailed mouse (Pseudomys higginsi) could segregate spatially, their α coefficients were small. During breeding and dispersal, when abundant populations precluded spatial segregation, their α coefficients were higher. Populations of Keen's mice (Peromyscus keeni) and dusky shrews (Sorex monticolus) are widely distributed across many of the islands in the Tongass National Forest (TNF) in Southeast Alaska [ 32 ]. Keen's mice are a medium-sized Peromyscus sp. (18±24 g) considered habitat and dietary generalists [ 33,34 ]. Dusky shrews (6±8 g) are found in a variety of habitats and are primarily insectivores but are also known to consume mushrooms and lichens [ 35,36 ]. Both species are mostly nocturnal and populations fluctuate annually with mice usually reaching higher densities than shrews [ 37 ]. Although not commonly considered competitors, these 2 / 21 species comprise the base of the food web in the coastal rainforests of Southeast Alaska and have the potential to influence each other's vital rates through competitive interactions. Our objective was to quantify competition between mice and shrews in the spatially heterogeneous landscape of the TNF. Because previous studies of Keen's mice have shown that these small mammals exhibit variable dynamics in different seral stages of regenerating forests [ 33 ], we compared competitive interactions among different habitat types in our study area. We hypothesized that due to their larger body size and generalist foraging strategies, Keen's mice would be a stronger competitor, exerting negative effects on the abundance of dusky shrews. We evaluated annual interspecific competition using the standardized regression technique [ 38 ] in 2010±2012. We hypothesized that coexistence of mice and shrews would be largely facilitated through partitioning of food resources and quantified dietary overlap and segregation between the two species using stable isotope analysis [ 39 ]. Finally, because theory predicts that whole community energy consumption will remain relatively static when interspecific competition alters species abundances (zero-sum dynamics; [ 40 ]), we predicted that shrew numbers would increase to energetically compensate for the competitive release from mice. Materials and methods Study area The TNF encompasses most of the Alexander Archipelago, a naturally fragmented system of large and small islands (Fig 1). Since the 1950's, large-scale logging created over 174,000 ha of young stands in the TNF [ 41 ]. In 2001, the Tongass-Wide Young-Growth Studies (TWYGS) were initiated to evaluate the response of plants and trees to pre-commercial thinning treatments [ 41 ]. The TWYGS comprise a randomized complete block design with replicates of four experimental treatments distributed across the TNF and established from 2002 to 2004. Manipulated sites on POW consisted of an un-thinned control (young-growth) stand paired with two thinned stands [ 41 ]. In 2010, we selected six study sites across northern Prince of Wales Island (POW), Alaska (55.9Ê N, 132.9Ê W; Fig 1); the largest island in the Alexander Archipelago at 6,674 km2. Our sites included seven TWYGS treatments, six un-thinned young-growth stands, three old-growth stands and three clearcuts (for a total of 19 stands). In 2011 and 2012 we added one old-growth and one clearcut grids (21 grids total). This sampling design encompassed a wide range of habitat conditions and successional stages across POW. Pre-commercially thinned stands had a robust understory with remaining trees spaced according to thinning treatment (e.g. 14 × 14 or 18 × 18 ft). Young-growth stands had little spatial variation and were composed of densely regenerated, even-aged, small-diameter trees with sparse understory. Old-growth stands ranged from low-elevation, high-productivity to highelevation, mesic sites. A variety of tree size classes and ages contributed to a more heterogeneous canopy in this stand type [ 42 ]. Clearcuts ranged in age from 1±15 years post-harvest, contained few remaining overstory trees, and were characterized by a dense layer of residual timber slash mixed with regenerating shrubs and tree seedlings. Annual precipitation on POW ranges from 254±508 cm and average temperatures range from 10±17ÊC in summer to 0±6ÊC in winter. Elevation at our study sites varied from 0 to 305 m. Areas lower than ~600 m are characterized by temperate, coniferous rainforest [ 43 ], dominated by Sitka spruce (Picea sitchensis) and western hemlock (Tsuga heterophylla). The understory consists primarily of Vaccinium spp., false azalea (Menziesia ferruginea), and salmonberry (Rubus spectabilis). The herbaceous layer includes skunk cabbage (Lysichiton americanum) and bunchberry (Cornus canadensis), among various other forbs, mosses, and ferns. Muskegs, or peatland bogs, compose the majority of natural non-forested areas on POW. 3 / 21 Fig 1. Study area and small mammal trapping grid locations on Prince of Wales Island, Alaska, USA, 2010± 2012. Trapping grids were located in seven thinned, six un-thinned young-growth, four old-growth and four clearcut stands. Similar to other high latitude islands, POW is characterized by a depauperate fauna with mink (Neovison vison) and ermine (Mustela erminea celenda) the only native mesocarnivores and few small mammals including Keen's mice, long-tailed voles (Microtus longicaudas), dusky shrews, and northern flying squirrels (Glaucomys sabrinus griseifrons) [ 32 ]. Long-tailed vole densities, which were high in the 1970s [ 33 ], are currently so low that in multiple years of attempts by several different research teams, only a handful have been captured ([ 44,45 ], S. MacDonald, pers. comm, this study). American marten (Martes americana) were introduced to POW in the 1930s to provide trapping opportunities for rural residents. The avian fauna of the island is also limited, with dark-eyed junco (Junco heymalis) the only ground-nesting songbird. Northern goshawk (Accipiter gentilis), a raptor that rarely preys on small mammals, is the main avian predator. Small mammal trapping We live-trapped small mammals from May to August 2010±2012 using 40 to 52 Sherman livetraps (H.B. Sherman Traps, Inc, Tallahassee, FL, USA) per grid set at 25-m intervals and configured to fit within each forest stand. We baited traps with a mix of oats, molasses and peanut butter. Polyester bedding was provided to aid in thermoregulation. Each year, trapping occurred for five consecutive nights three times in late spring, early summer and late summer/ fall following a robust-design capture-recapture protocol [ 46 ]. We checked traps once in the 4 / 21 morning. Each trapped individual was weighed, sexed, aged, and assessed for reproductive status. Mice were marked with a passive integrated transponder (PIT) tag (Biomark, Boise, ID, USA) for permanent identification after a brief exposure to Isoflurane (Piramal Healthcare Limited, Andhra Pradesh, India). After immobilization, we collected blood via tail clipping [ 47 ] from all mice at first capture. Live-captured shrews were released without processing. Carcasses of all incidental mortalities (255 mice, 1330 shrews) were deposited in the Vertebrate Museum at the University of Wyoming (UW) Berry Center for Biodiversity Conservation. All procedures were approved by the Institutional Animal Care and Use Committee at UW (see S1 Text) and adhere to the guidelines of the American Society of Mammalogists [ 48 ]. Our trapping permit was obtained from the Alaska Department of Fish and Game and study methods were approved by the US Forest Service, Tongass National Forest, Thorne Bay District. This study did not involve endangered or protected species. Population modeling We generated estimates of survival and abundance for Keen's mice from mark-recapture data using the robust-design population model in Program MARK [ 49 ]. We constructed several a priori competing models (Table A in S1 Text) and selected among them based on Akaike's information criterion values corrected for sample size (AICc), and the significance of parameter estimates. Models were developed separately for each grid for the three years and 45 occasions (three primary and five secondary occasions in each sampling year). We used unequal time intervals (1 month between the summer primary sessions and 10 months for the overwinter periods) to generate survival estimates. For dusky shrews we used a dead-and-alive framework to generate abundance estimates. After removing the trap mortalities from the dataset for each grid, we calculated the abundance of the live animals with the Horvitz-Thompson estimator [ 46 ]. Because we had no marked individuals, we used estimates of capture probability (and 95% confidence intervals) for Sorex shrews from a study by Otto and Roloff [ 50 ]. We treated the number of live captures as minimum number known alive (MNKA). We then added the number of trap mortalities to the live estimates to generate the total abundance estimates. To assess the potential effect of trap mortality on the shrew population, we calculated growth rate (lambda) from abundance estimates by dividing values in sessions 2 and 3 by sessions 1 and 2 in each year. We then regressed the resulting lambda values against trap mortality in the preceding session. Data for additional population indices (MNKA and captures per 100 trap nights [100TN]), as well as the relationships among them for mice (S1 Fig) and shrews (S2 Fig) are provided in the supplementary materials. We used a traditional boundary-strip method to estimate mouse densities. We calculated the mean maximum distance moved (MMDM; average 52 m) by measuring the distances traveled between subsequent captures [ 51 ]. We calculated the effective trapping area (ETA) by adding a 50 m buffer to the area covered by live-traps. We then divided the abundance estimate by the ETA for each trapping session per grid. Shrew densities were calculated by dividing the estimated abundance by the effective trapping area calculated for each grid based on mice MMDM. On a subset of grids we used spatially explicit capture-recapture functions [ 52 ] to estimate density of mice with the package secr in the R statistical environment ([ 53,54 ]; S1 Text). Similar to observations by Gerber and Parmenter [55], our models yielded higher estimates compared to MMDM (S3 Fig). Given these results and our interest in vital rates we did not further pursue estimates using the Efford and Fewster [ 52 ] methods. 5 / 21 Habitat measurements We conducted vegetation sampling three times each summer at eight randomly selected trapping stations in each grid in 2010 and twice each summer at nine trapping stations in 2011 and 2012. For understory components, two 20-m transects were placed in the four cardinal directions, crossing at the center of the trap station. We used line-intercept methods to record percent cover of shrub species, coarse woody debris, and herbaceous plants along transects in all years. In 2010, at the end of each transect, we sampled a 1 × 1 m plot to obtain occurrence for shrubs and herbaceous species. In all years, we collected and weighed all epigeous fungi, truffles, and earthworms, and counted the number of conifer cones in the same plots. In 2010, pitfall traps were used to sample soil macro-invertebrates (mostly beetles) every 5 m along one randomly selected transect. All invertebrates were later identified to family. To characterize the overstory, we measured leaf-area index (LAI) at each sampled trap station in 2010 [ 42 ]. Averages of understory and overstory variables represent stand-level values for each trapping grid. We obtained the elevation at each grid from a digital elevation model using ArcGIS software [ 56 ]. Each grid was characterized by 11 variables (Table B in S1 Text). We reduced the dimensionality of the habitat variables by non-metric multidimensional scaling (NMDS; [ 57 ]) in Program R. We used variable scores (r > 0.40) to interpret the dimensions. Loading of these variables on NMDS1 and NMDS2 are provided in Table C in S1 Text. Competition analysis We used densities of mice and shrews and habitat variables to estimate competition coefficients (α). We first standardized density estimates for each species to a mean of zero and a standard deviation of one [ 38 ]. This eliminated any dependence of α on the density variances [ 58 ]. We included NMDS scores to represent habitat variables. We also included elevation as a separate predictor variable in candidate models because of high variation. The effect of Keen's mouse density and habitat variables on dusky shrew density can be estimated by: Ys ˆ a ‡ bpXp ‡ b1X1 ‡ ‡ bnXn …1† where Ys is the density of shrews, b1. . .bn and X1. . .Xn are the coefficients and associated predictor variables, respectively, Xp is the density of mice, and bp is the α coefficient, estimating the effect of mice on shrews. In these analyses we used density estimates for both species of small mammals from each session and overstory variables from 2010 (assuming overstory cover changed little over the course of the study) and understory variables derived for each grid each summer. Because animal abundance is a function of the net of reproduction and mortality, we also constructed models with the density of shrews in the previous session as one of the potential explanatory variables [i.e., auto-regressive models; 59]: Ys ˆ a ‡ bs…t 1†Xs…t 1† ‡ bpXp ‡ b1X1 ‡ ‡ bnXn …2† where bs(t-1) is the regression coefficient and Xs(t-1) is the density of shrews in the previous session. We only included density estimates from sessions 2 and 3 of each year as the dependent variables and session 1 and 2 as the independent ones. Because mouse density varied by year, we performed separate regressions for 2010, 2011, and 2012. We used generalized linear models (glm) in R to construct models and compare among them using AICc, R2, and parsimony. We created 15 regressive (Table D in S1 Text) and 16 auto-regressive models containing all configurations of NMDS dimensions, elevation, and densities of each species. 6 / 21 Stable isotope analysis To assess dietary overlap of shrews and mice, we used stable isotopes of carbon (δ13C) and nitrogen (δ15N). Other methods were not applicable because stomachs and feces of both species largely contained bait. We analyzed whole blood collected from mice and muscle tissue from mice and shrew mortalities. We sampled diet items including macro-invertebrates (beetles and earthworms), conifer seeds, berries, lichen, and fungi from all grids. We dried whole blood (mice), muscle (mice and shrews), and potential food samples at 60ÊC for 48 h and then homogenized them in a ball mill (Mixer Mill MM200, Retsch Inc., Newtown, PA, USA). Samples were weighed in duplicate and sent to the University of Wyoming Stable Isotope Facility. Analyses of δ13C and δ15N were conducted with a Costech 410 elemental analyzer (Costech Analytical, Valencia, CA, USA) connected to a Thermo Finnigan DeltaPLUS XP Continuous Flow Isotope Ratio Mass Spectrometer (Thermo Fisher Scientific, Inc., Waltham, MA, USA). Results, reported as parts per thousand (½) in relation to internal laboratory standards, were accepted (and averaged) when the variance between duplicate samples was lower than that of the standards within each run. We compared the isotope values of whole blood and muscle tissue from mice using a Knearest neighbor randomization test [ 60 ] in R to ensure they were equivalent. We also tested for differences among the diet items and compared the δ13C and δ15N of mice and shrews by habitat type (young growth, thinned, old growth and clearcut) each year using the KNN test. The isotope data were converted to estimates of relative contributions of dietary items separately for mice and shrews in each habitat type using the multiple source, dual-isotope linear mixing-model package SISUS in R [ 61 ]. Because the variation in shrew signatures was low we used population means per grid for shrews. In contrast, we calculated dietary contribution for individual mice because of the large variation associated with those data. We corrected source data for consumer-diet discrimination by 1½ for δ13C and 3½ for δ15N for mice [ 62,63 ] and 1.3½ for δ13C and 2.5½ for δ15N for shrews [64]. We implemented concentration-dependent models to account for differences in C:N of sources. We then calculated the overlap in diet between each mouse and the population of shrews in each grid and considered it significant when both species included a contribution > 20% of a diet item. To describe the isotopic niches of the two species we used kernel utilization density methods (KUD; [ 65 ]). We extracted the 50%, 75% and 95% contours with the rKIN package in R [ 66 ] from the probability surface (S4 Fig). We then calculated niche size and percent overlap for mice and shrews with rKIN in 2010±2011 and 2012. Because the number of samples available in 2012 was low (21 mice and 64 shrews; as a result of fewer caught mice and lower trap mortality in shrews) we randomly subsampled the 2010±2011 data (280 mice and 231 shrews) with the number of points in 2012. We calculated niche size and percent overlap for this subsample and compared it to the original, full dataset. Whole community energy consumption Theoretically, when released from competition, the weaker competitor should increase in numbers such that the change in abundance is proportional to the inverse of average metabolic rate [ 40 ]. Because the metabolic rate of shrews is higher than expected based on allometric expectations (i.e., of body mass; [ 67 ]), it is likely that their numerical increase will be lower than predicted by body mass alone. To assess the effects of declines in Keen's mice, the dominant competitor, on dusky shrews we estimated the whole community energy consumption by multiplying the abundance of each species by the average metabolic rate of individuals in each year of study. We calculated this value for the last trapping session in each year (i.e., August) when the abundance of both species reached a maximum. For shrews we used values for field 7 / 21 metabolic rate (FMR; Kj per day) published by Ochocińska and Taylor [ 68 ] for common shrews (Sorex araneus). Because common shrews are slightly larger (8.2g) than dusky shrews we converted their estimates to mass specific FMR and used the average body mass of shrews in each year to generate values for individuals. For mice we used FMR reported by Degen [ 69 ] for Peromyscus maniculatus, which is closely related to P. keeni [ 70 ]. The average body mass of Keen's mice in 2010, 2011, and 2012 was 21.2g, 20.3g, and 20.9g, respectively. We then calculated the total energy consumption of the community by summing the values of mice and shrews in each year. Results Abundance and density We captured a total of 2,543 individual Keen's mice and 1,744 dusky shrews in all habitats and study sites in all years. In contrast, we captured only 18 long-tailed voles and 36 northern flying squirrels. There was strong correlation between abundance estimates and indices of population size for both species (supplementary materials S1 and S2 Figs). Abundance estimates show that mouse populations in our study area peaked in 2011 and crashed in 2012 (Fig 2). Concomitantly, high overwinter mortality of dusky shrews was followed by steady increases during the breeding season in all years, reaching the highest abundance in 2012 (Fig 2). This increase was pronounced in all habitat types except thinned stands. Shrew body mass was significantly lower and trap mortality higher in 2010 than in 2012 (Fig 3). Nonetheless, trap mortality had little effect on population growth rate of shrews (Fig 3). In contrast, mortality of mice in traps was low (2010: 13.5%, 2011: 6.5%, 2012: 13.0%) and unrelated to body mass or population growth. There was a significant inverse relationship of Keen's mice and dusky shrew densities in 2010 and 2011 (p = 0.03 and p < 0.001, respectively; Fig 4), with low shrew densities in areas characterized by high mice densities (Fig 4). There was no relation between mice and shrew densities in 2012 (p = 0.16). Competition models Habitat variables were reduced to two dimensions (Table C in S1 Text) representing overstory and understory cover, and species diversity, respectively. These two dimensions separated the four habitat types with most of the divergence occurring on NMDS1 and some on NMDS2, especially between old growth and clearcuts (Fig 5). There was a positive relationship between shrew densities in time t (sessions 2 or 3) and those in time t-1 (sessions 1 or 2; p < 0.001; Fig 6). Additionally, the auto-regressive models exhibited higher R2 values than the regressive ones (Table 1 and Table D in S1 Text), suggesting that models excluding information on shrews did not sufficiently explain the variation in their densities. In 2010 and 2011, the top auto-regressive models only included the density of mice as an added predictor (Table 1). The associated α coefficients indicate a large increase in the influence of mice on shrews in 2011 (α = -0.14 and -0.90, respectively). In 2012, the best fitting model only included shrew density in t-1 and elevation, with densities increasing with altitude (Table 1). Diet composition and niche overlap We found no significant differences in the δ13C and δ15N of whole blood and muscle tissue from mice (p > 0.05) and combined these samples. All diet items in each habitat type exhibited 8 / 21 Fig 2. Abundance estimates (± 95% confidence intervals) for Keen's mice and dusky shrews. Abundance was estimated in each of four habitat types on Prince of Wales Island, Alaska. Every grid was trapped during three sessions (S1-S3) each summer from 2010±2012. significant differences in δ13C and δ15N (p < 0.05). Macro-invertebrates were pooled into a single group with the exception of earthworms, which were isotopically distinct (Fig 7). The variance in diets of mice was double that observed in shrews (Fig 7). Earthworms and other macro-invertebrates had the largest contributions to the diets of these small mammals (29±58% and 60±71%, respectively; Table 2). Additionally, in young-growth stands, mushrooms substantially contributed (34%) to the assimilated diet of dusky shrews. The proportion of mice that included > 20% of invertebrates, mushrooms, or earthworms in their diet varied, with the highest overlap with shrew diets occurring in young-growth and old-growth stands (Table 2). Mice niche (KUD) was nearly 5× larger within the 95% contour than within the 50% contour, whereas for shrews the 95% niche was only twice that of the 50% (Table 3), indicating that their diets were more uniform. Subsampling of the 2010±2011 dataset with the number of data point available in 2012 demonstrated that the difference between the two periods was largely a result of contraction of the dietary niche in mice and a shift in the niche of shrews rather than an artifact of sample size (Table 3). At the 50% level, the niche size of mice overlapped largely with that of shrews expect in old-growth stands; at 95% shrew niches were nearly completely encompassed by mice in all habitats (Fig 8). Between 2010±2011 and 2012, niche size of mice decreased by 32±33% in all contours (Table 3). Concurrently the percent overlap between the two species increased from 13% of the dusky shrews 95% contour to a considerable 53% in the core 50% (Table 3). Dusky shrew core niche (50%) in 2012 only overlapped with that of 2010±2011 by 53.4%, suggesting a large shift in resource use between these 9 / 21 Fig 3. Body mass, mortality of trapped shrews, and relationship to population growth rate. Average body mass (bars; ± 95% confidence intervals), trap mortality rate (line; top), and the relation between trap mortality and population growth rate (bottom) of dusky shrews captured in 2010±2012 on Prince of Wales Island, Alaska. Letters indicate significant differences. 10 / 21 Fig 4. Relationship of shrew and mouse densities. Density of dusky shrews on Prince of Wales Island, Alaska from 2010±2012 in relation to the density of Keen's mice. two periods. This was also noticeable in the 75% (64.5% overlap) and 95% contour (69.1% overlap). Whole community energy consumption The total energy consumption of Keen's mice and dusky shrews increased between 2010 and 2011 (Fig 9) mirroring the increase in mice abundance in all habitats (Fig 2). It then declined in 2012 and reached approximately 69 and 68% of the 2010 and 47 and 45% of the 2011 values in thinned and young-growth stands, respectively. In old growth, total energy consumption reached 88% of 2010 and 63% of 2011 levels, whereas in clearcuts these values were 95% and 67%, respectively. This suggests that in these two habitats, shrews compensated for the declines in mice abundance more significantly than in thinned and young-growth stands. Discussion Our results provide several lines of evidence that interspecific competition influences the densities of mice and shrews in Southeast Alaska, but that the effects of competition are spatially and temporally dynamic. Similar to removal experiments, we observed the expected increase in densities of dusky shrews after mouse populations crashed in 2012. Although shrew abundance in 2012 increased above 2010±2011 levels in all habitat types, this was most pronounced in clearcut and old-growth stands. Indeed, in these two habitats the increase in shrews partially compensated energetically for the declines in mice. In addition, α coefficients revealed that the influence of Keen's mice, the dominant competitor, varied congruent with their densities. The numerical responses of shrews occurred in conjunction with habitat-specific segregation in 11 / 21 Fig 5. Values of non-metric multidimensional scaling (NMDS). NMDS1 and NMDS2 values for habitat and food variables in each of four habitat types on Prince of Wales Island, Alaska from 2010±2012. diet. When mouse populations were low, niche overlap between mice and shrews increased as the diet of shrews shifted, indicating that in the absence of their competitor, shrews were able to exploit previously unavailable resources. Our results clearly illustrate that competition can occur even among species that differ in size and diet and that the evolutionary adaptations of one competitor may limit the ability of the community to compensate energetically. These complex interactions, observed in a relatively simple ecological system, highlight the difficulty in assessing the role of biotic factors in structuring communities and may explain the recent finding of strong abiotic masking of competition and negative covariances [ 9 ]. In 2010 estimated α coefficients suggested weak competitive interactions concurrent with high dietary overlap between mice and shrews. Intraspecific competition for food regulates small mammal populations in other systems [ 71 ]. Although we were unable to measure food limitation per se, our habitat variables included measures of food availability. Biomass of earthworms and numbers of beetles exhibited low scores on the two NMDS axes, both of which explained little of the variation in shrew densities. In addition, in all years, there was a positive relationship between shrew density in the current session and that in the following one, negating expectations of intraspecific competition. Nonetheless, in 2010, shrew body condition was low. It is possible that competition with mice for specific foods that enhance body condition in 12 / 21 Fig 6. Relationship between shrew density in the current session vs the previous one. Density of dusky shrews in the current trapping session (2 or 3) on Prince of Wales Island, Alaska, in relation to their density in the previous session (1 or 2) in each year from 2010±2012. The relationship is shown for all years combined. shrews and that were less available in 2010 resulted in competitive effects at the individual, rather than the population level. Alternatively, other independent processes such as parasite load or disease, [ 72,73 ] may have influenced shrew populations in 2010. Although habitat variables did not exert a strong influence on competition, the four habitat types differed greatly in canopy and horizontal cover (NMDS1) and food availability (NMDS2). In addition, niche overlap and the numerical responses of shrews to competition For all years, the dependent variable was shrew density (SOMO) at time t. Independent variables were shrew density at time t-1 (trapping sessions 1 or 2), Keen's mice density (PEKE) at time t (trapping sessions 2 or 3), elevation at the trapping grid, NMDS1 and NMDS2. Only models with considerable support (ΔAIC 2) are shown. In bold are top models based on AICc, R2 and parsimony (also see Fig 4). 13 / 21 Fig 7. Stable carbon and nitrogen values for Keen's mice, dusky shrews and the potential foods. Values of δ13C and δ15N (average ± 95% confidence intervals) for blood and muscle tissues of Keen's mice and dusky shrews captured on Prince of Wales Island, Alaska 2010±2012. Potential diet items were corrected for consumer-diet discrimination by using a trophic enrichment factor of 1.3½ for δ13C and 2.5½ for δ15N in relation to shrew diets. with mice differed among habitats. In old-growth stands in 2010±2011, there was little overlap in diets of mice and shrews. Concurrently, in these grids shrew density was low until mice numbers crashed, after which they substantially increased. In contrast, in clearcuts, the isotopic niche of shrews was nearly fully encompassed by that of mice while their numerical responses were similar to old growth. This type of habitat-dependent competition has also been shown in Australian rodents where Rattus lutreolus and Pseudomys gracilicaudatus Proportion in the diet and percent of individuals was estimated from δ13C and δ15N of shrews and mice and their potential prey using a dual-isotope linear mixing model. Diet items were corrected using a discrimination factor of 1.3½ for δ13C and 2.5½ for δ15N for shrews and 1½ for δ13C and 3½ for δ15N for mice. Sample sizes are provided in parentheses. Bolded food items constituted > 20% of the assimilated diet of shrews. 14 / 21 Data from 2010±2011 were randomly subsampled 5 times for the sample size of 2012. Confidence intervals (95%) are reported in parentheses. coexist in both wetter and drier habitats [ 26 ]. In drier sites, they displayed symmetric competition but within wetter sites, Rattus reduced the density of Pseudomys with no mutual effect. Despite the vast differences in habitat structure (NMDS1) and food availability (NMDS2) between thinned and young-growth stands, the numerical response of shrews to the crash of Fig 8. Isotopic niche size and percent overlap. Niche size and percent overlap of Keen's mice and dusky shrews on Prince of Wales Island, Alaska 2010±2011 (bottom) from the 50% and 95% kernel utilization density contours (see also supplementary materials S4 Fig) in each of four habitat types. 15 / 21 Fig 9. Total energy consumption of mice and shrews (Kj/day). Energy consumption was calculated based on abundance and species-specific field metabolic rates in four habitat types on Prince of Wales Island, Alaska 2010±2012. The increase in shrews partially compensated for the decline in mice in 2012. mice was dampened and the mice isotopic niche overlapped less with that of shrews, especially in thinned stands. In these habitats, the increase in shrew abundance was insufficient to compensate energetically for the decline in mice, especially compared with 2011. It is possible that in young-growth stands shrews had low food availability even in the absence of mice, which limited their potential increase. In contrast, thinned stands encompass much of the diversity, and thus food availability, observed in the other three habitat types. Nonetheless, although we and other studies [ 36 ] have shown that dusky shrews consume mushrooms and lichens, as insectivores they are unlikely to take advantage of all food types consumed by mice (e.g., berries). This may also explain our observation that shrews only partially compensated energetically for the decline in mice even in old-growth stands and clearcuts. Others have suggested that limitations on resource use could lead to differences in energy consumption in ecosystems despite similarities in species richness [ 74 ]. Alternatively, winter declines and small populations in spring may have prevented shrews from reaching numbers that compensate for the absence of mice. Previous studies have also demonstrated that shrews usually reach lower densities than mice [ 37 ]. Although our study was largely observational, the severe decline in mice densities in 2012 acted as a pseudo-removal experiment. The observed changes in the degree of overlap of the dietary niches of mice and shrews and within shrews between 2010±2011 and 2012 indicate a large shift in resource use. Our subsampling of the 2010±2011 dataset (280 mice and 231 shrews) with the number of mice and shrews represented in 2012 (21 mice and 64 shrews), provides support for this conclusion. Therefore, the differences in sample size could not explain the dietary shift observed in dusky shrews. Indeed, a shift rather than an increase in the size of the dietary niche is consistent with our observation that the numerical increase of 16 / 21 dusky shrews only partially compensated for the changes in community energy consumption following the crash of Keen's mice. The crash of Keen's mouse populations on POW followed a winter with relatively high snow accumulation. During the 2011±2012 winter, accumulation was 4× higher than in 2009± 2010 and 2.5× higher than 2010±2011 (http://www.wrcc.dri.edu/cgi-bin/cliMAIN.pl?ak2227). Severe winters can reduce small mammal survival by limiting access to food and subnivean spaces, which may increase the costs of thermoregulation and predation risk [ 75,76 ]. However, snow accumulation alone cannot explain the lack of recovery of mice during the following breeding season. Also, if snow caused the crash of mice, we would expect to see the same trend in shrews. We suspect that predation, possibly facilitated by high snow accumulation, was the main cause of the demise of Keen's mice. Our companion study of mesocarnivores revealed that abundance of both predators increased after 2010, with ermines peaking in 2011 and martens in 2012 [77]. Because martens rarely prey on shrews [ 78 ], they largely targeted mice in 2012 [ 77 ] and likely inhibited their recovery over the summer. Regular predator-prey cycles of martens, mice and long-tailed voles have been observed in Southeast Alaska [ 79 ]. Given the size of POW and recent findings that competition on large islands is low compared to medium size islands [ 23 ], our observation of the effect of Keen's mice abundance on dusky shrews numbers and diet is significant. Schoener et al. [ 23 ] suggest that low effect size of competitive interactions on large islands is caused by spatial heterogeneity of habitat features and food availability. Nonetheless, despite its size, as a high latitude island POW is characterized by depauperate flora and fauna where heterogeneity largely stems from anthropogenic disturbance. It is likely that biotic interactions, even between unlikely competitors, may be more pronounced in depauperate systems and thus more readily detected. The combination of temporal and spatial variation in competitive interactions we documented illustrates the importance of management practices in the TNF on community dynamics. We observed that in old-growth stands, shrews were better able to energetically compensate for the decline in mice. Although similar increases occurred in clearcuts, forest succession quickly converts them to young-growth stands. Thus, continued old-growth logging in the TNF will result in reduction in the capacity of shrews to recover from competitive pressures imposed by mice. In contrast, in pre-commercially thinned stands, shrew populations were less affected by changes in mice densities. Because shrews constitute the main prey for the endemic POW ermine during scarcity of alternative prey [ 77 ], their relative stability in thinned, open-canopy stands ensures food supply for these predators. This suggests that thinning may reverse the homogeneity resulting from clearcut logging, increase habitat quality, and help to conserve species diversity, especially of ground-dwelling small mammals. Our results support earlier studies by Carey and Wilson [ 80 ], Sullivan and Sullivan [ 37 ], and Suzuki and Hayes [ 81 ], which concluded that pre-commercial thinning or selective harvesting increased habitat quality and abundance for several small mammal species. At this juncture in the development of the new Tongass Land Management Plan, our findings may provide an incentive for accelerating the transition from old-growth clearcut logging to commercial harvest of young-growth stands. Supporting information S1 Text. Supplementary text and tables for the manuscript. (DOCX) S1 Fig. Comparison of Keen's mice abundance estimates. Abundance estimates of Keen's mice generated from the robust-design population model in Program MARK (top) and captures per 100 trap nights (100TN; bottom) in relation to the minimum number known alive 17 / 21 (MNKA) on Prince of Wales Island, Alaska from 2010±2012. (PNG) S2 Fig. Comparison of dusky shrew abundance estimates. Abundance estimates of dusky shrews generated from a dead-and-alive framework (top) and captures per 100 trap nights (100TN; bottom) in relation to the minimum number known alive (MNKA) on Prince of Wales Island, Alaska from 2010±2012. (PNG) S3 Fig. Comparison of Keen's mice density estimates. Density (number per ha) estimates of Keen's mice generated from the robust-design population model in Program MARK in relation to estimates generated from spatially explicit capture-recapture (SECR) on Prince of Wales Island, Alaska from 2010±2012. (PNG) S4 Fig. Niche size and overlap with kernel density estimates. 50%, 75%, and 95% contours of the isotopic niches of Keen's mice and dusky shrews on Prince of Wales Island, Alaska 2010± 2011 estimated with kernel density estimators. Mice are depicted in circles and shades of orange and shrews in triangles and purple. (PNG) Acknowledgments We thank K. Anderson, C. Bickford, V. Clarkston, S. Curtis, J. Dillon, J. Furlaud, M. Hobbs, A. La Porte, P. Mitsdarfer, S. Opitz, E. Reuber, T. Salzman, and M. Youngblood for help in the field. M. Cady and R. Slayton provided logistical support. J. Evans, J. Goheen, D. Legg, R. Liu, and M. Murphy provided comments on early versions of the manuscript. Funding was provided by the United States Forest Service Tongass National Forest, the National Fish and Wildlife Foundation, the UW Biodiversity Institute and Department of Zoology and Physiology. Author Contributions Data curation: Carolyn A. Eckrich. Formal analysis: Carolyn A. Eckrich, Merav Ben-David. Funding acquisition: Elizabeth A. Flaherty, Merav Ben-David. Methodology: Elizabeth A. Flaherty. Writing ± original draft: Carolyn A. Eckrich. Writing ± review & editing: Elizabeth A. Flaherty, Merav Ben-David. 18 / 21 19 / 21 20 / 21 1. MacArthur R , Levins R. The limiting similarity, convergence, and divergence of coexisting species . Am Nat . 1967 ; 377 ± 385 . 2. Adler PB , HilleRisLambers J , Levine JM . A niche for neutrality . Ecol Lett . 2007 ; 10 : 95 ± 104 . https://doi. org/10.1111/j.1461- 0248 . 2006 . 00996 . x PMID : 17257097 3. Brown WL , Wilson EO. Character Displacement . Syst Zool . 1956 ; 5 : 49 ± 64 . https://doi.org/10.2307/ 2411924 4. Cavender-Bares J , Kozak KH , Fine PVA , Kembel SW . The merging of community ecology and phylogenetic biology . Ecol Lett . 2009 ; 12 : 693 ± 715 . https://doi.org/10.1111/j.1461- 0248 . 2009 . 01314 . x PMID : 19473217 5. Hutchinson GE. Concluding Remarks . Cold Spring Harb Symp Quant Biol . 1957 ; 22 : 415 ± 427 . 6. Volkov I , Banavar JR , Hubbell SP , Maritan A . Neutral theory and relative species abundance in ecology . Nature . 2003 ; 424 : 1035 ± 1037 . https://doi.org/10.1038/nature01883 PMID: 12944964 7. Ernest SKM , Brown JH . Homeostasis and compensation: the role of species and resources in ecosystem stability . Ecology . 2001 ; 82 : 2118 ± 2132 . https://doi.org/10. 1890 /0012- 9658 ( 2001 ) 082 [2118: HACTRO] 2.0 .CO; 2 8. Goheen JR , White EP , Ernest SKM , Brown JH . Intra-guild compensation regulates species richness in desert rodents . Ecology . 2005 ; 86 : 567 ± 573 . https://doi.org/10. 1890 /04- 1475 9. Houlahan JE , Currie DJ , Cottenie K , Cumming GS , Ernest SKM , Findlay CS , et al. Compensatory dynamics are rare in natural ecological communities . Proc Natl Acad Sci . 2007 ; 104 : 3273 ± 3277 . https://doi.org/10.1073/pnas.0603798104 PMID: 17360637 10. Gause GF . The struggle for existence . New York, NY: Hafner Press; 1934 . 11. Valone TJ , Brown JH. Effects of competition, colonization, and extinction on rodent species diversity . Science . 1995 ; 267 : 880 ± 883 . https://doi.org/10.1126/science.7846530 PMID: 7846530 12. Schoener TW . Resource partitioning in ecological communities . Science . 1974 ; 185 : 27 ± 39 . https://doi. org/10.1126/science.185.4145.27 PMID: 17779277 13. Eccard JA , YloÈnen H . Interspecific competition in small rodents: from populations to individuals . Evol Ecol . 2003 ; 17 : 423 ± 440 . https://doi.org/10.1023/A:1027305410005 14. Brown JS . Desert rodent community structure: A test of four mechanisms of coexistence . Ecol Monogr . 1989 ; 59 : 1± 20 . https://doi.org/10.2307/2937289 15. Abramsky Z , Rosenzweig ML , Brand S. Habitat selection of Israel desert rodents: Comparison of a traditional and a new method of analysis . Oikos . 1985 ; 45 : 79 ± 88 . https://doi.org/10.2307/3565225 16. Songer MA , Lomolino MV , Perault DR . Niche dynamics of deer mice in a fragmented, old-growth-forest landscape . J Mammal . 1997 ; 1027 ± 1039 . 17. LemaÃõtre J , Fortin D , Morris DW , Darveau M. Deer mice mediate red-backed vole behaviour and abundance along a gradient of habitat alteration . Evol Ecol Res . 2010 ; 12 : 203 . 18. Holt RD . Predation, apparent competition, and the structure of prey communities . Theor Popul Biol . 1977 ; 12 : 197 ± 229 . https://doi.org/10.1016/ 0040 - 5809 ( 77 ) 90042 - 9 PMID: 929457 19. Hanski I , Henttonen H. Predation on competing rodent species: A simple explanation of complex patterns . J Anim Ecol . 1996 ; 65 : 220 ± 232 . https://doi.org/10.2307/5725 20. Gurevitch J , Morrison Janet A. , Hedges LV . The interaction between competition and predation: A meta-analysis of field experiments . Am Nat . 2000 ; 155 : 435 ± 453 . https://doi.org/10.1086/303337 PMID: 10753073 21. Ben-David M , Bowyer RT , Faro JB . Niche separation by mink and river otters: coexistence in a marine environment . Oikos . 1996 ; 75 : 41 ± 48 . https://doi.org/10.2307/3546319 22. Morris DW . A haunting legacy from isoclines: Mammal coexistence and the ghost of competition . J Mammal . 1999 ; 80 : 375 ± 384 . https://doi.org/10.2307/1383286 23. Schoener TW , Spiller DA , Piovia-Scott J . Variation in ecological interaction strength with island area: theory and data from the Bahamian archipelago . Glob Ecol Biogeogr . 2016 ; 25 : 891 ± 899 . https://doi. org/10.1111/geb.12297 24. Fretwell SD , Lucas HL Jr. On territorial behavior and other factors influencing habitat distribution in birds . Acta Biotheor . 1970 ; 19 : 16 ± 36 . 25. Ziv Y , Abramsky Z , Kotler BP , Subach A . Interference competition and temporal and habitat partitioning in two gerbil species . Oikos . 1993 ; 66 : 237 ± 246 . https://doi.org/10.2307/3544810 26. Morris DW , Fox BJ , Luo J , Monamy V . Habitat-dependent competition and the coexistence of Australian heathland rodents . Oikos . 2000 ; 91 : 294 ± 306 . 27. Eccard JA , Fey K , Caspers BA , YloÈnen H . Breeding state and season affect interspecific interaction types: indirect resource competition and direct interference . Oecologia . 2011 ; 167 : 623 ± 633 . https://doi. org/10.1007/s00442-011-2008-y PMID: 21597944 28. Tennant EN , Germano DJ . Competitive interactions between Tipton and Heermann's kangaroo rats in the San Joaquin Valley , California. Southwest Nat. 2013 ; 58 : 258 ± 264 . https://doi.org/10. 1894 /0038- 4909- 58 .2. 258 29. Meserve PL , Gutierrez JR , Yunger JA , Contreras LC , Jaksic FM . Role of biotic interactions in a small mammal assemblage in semiarid Chile . Ecology . 1996 ; 77 : 133 . 30. Rosenzweig ML . A theory of habitat selection . Ecology . 1981 ; 62 : 327 ± 335 . https://doi.org/10.2307/ 1936707 31. Luo J , Monamy V , Fox BJ . Competition between two Australian rodent species: A regression analysis . J Mammal . 1998 ; 79 : 962 ± 971 . 32. MacDonald SO , Cook JA . Mammals and amphibians of Southeast Alaska . Museum of Southwestern Biology , University of New Mexico; 2007 . 33. VanHorne B. Niches of adult and juvenile deer mice (Peromyscus maniculatus) in seral stages of coniferous forest . Ecology . 1982 ; 63 : 992 ± 1003 . https://doi.org/10.2307/1937239 34. Smith WP , Nichols JV . Demography of two endemic forest-floor mammals of Southeastern Alaskan temperate rain forest . J Mammal . 2004 ; 85 : 540 ± 551 . https://doi.org/10.1644/BEH-003 35. Ingles LG . A quantitative study on the activity of the dusky shrew (Sorex vagrans obscurus) . Ecology . 1960 ; 41 : 656 . https://doi.org/10.2307/1931797 36. Gunther PM , Horn BS , Babb GD . Small mammal populations and food selection in relation to timber harvest practices in the western Cascade Mountains . Northwest Sci . 1983 ; 57 : 32 ± 44 . 37. Sullivan TP , Sullivan DS . Influence of variable retention harvests on forest ecosystems. II. Diversity and population dynamics of small mammals . J Appl Ecol . 2001 ; 38 : 1234 ± 1252 . 38. Fox BJ , Luo J . Estimating competition coefficients from census data: a re-examination of the regression technique . Oikos . 1996 ; 77 : 291 ± 300 . https://doi.org/10.2307/3546068 39. Ben-David M , Flaherty EA . Stable isotopes in mammalian research: a beginner's guide . J Mammal . 2012 ; 93 : 312 ± 328 . https://doi.org/10.1644/11-MAMM-S- 166 . 1 40. Ernest SKM , Brown JH , Thibault KM , White EP , Goheen JR . Zero sum, the niche, and metacommunities: Long-term dynamics of community assembly . Am Nat . 2008 ; 172 : E257±E269 . https://doi.org/10. 1086/592402 PMID: 18947326 41. McClellan MH . Adaptive management of young stands on the Tongass National Forest . U S For Serv Gen Tech Rep Pnw . 2008 ; 733 : 225 . 42. Eckrich CA , Flaherty EA , Ben-David M . Estimating leaf area index in Southeast Alaska: A comparison of two techniques . PLoS ONE . 2013 ; 8: e77642 . https://doi.org/10.1371/journal.pone. 0077642 PMID: 24223718 43. Alaback PB . Dynamics of understory biomass in Sitka spruce-western hemlock forests of Southeast Alaska . Ecology . 1982 ; 63 : 1932 ± 1948 . https://doi.org/10.2307/1940131 44. Flynn RW , Schumacher TV , Ben-David M . Abundance, prey availability and diets of American martens: Implications for the design of old-growth reserves in Southeast Alaska . Juneau, Alaska: Alaska Department of Fish and Game; 2004 p. 48 . 45. Pauli JN , Whiteman JP , Riley MD , Middleton AD . Defining noninvasive approaches for sampling of vertebrates . Conserv Biol . 2010 ; 24 : 349 ± 352 . https://doi.org/10.1111/j.1523- 1739 . 2009 . 01298 . x PMID : 19624526 46. Amstrup SC , McDonald TL , Manly BF . Handbook of capture-recapture analysis . Princeton, N.J: Princeton University Press; 2010 . 47. Abatan OI , Welch KB , Nemzek JA . Evaluation of saphenous venipuncture and modified tail-clip blood collection in mice . J Am Assoc Lab Anim Sci . 2008 ; 47 : 8± 15 . 48. Sikes RS , Bryan JA , Byman D , Danielson BJ , Eggleston J , Gannon MR , et al. 2016 Guidelines of the American Society of Mammalogists for the use of wild mammals in research and education . J Mammal . 2016 ; 97 : 663 ± 688 . 49. Cooch EG , White GC . Program MARK : A gentle introduction . 16th edition. [Internet] . 2015 . Available: http://www.phidot.org/software/mark/docs/book/ 50. Otto CRV , Roloff GJ . Using multiple methods to assess detection probabilities of forest-floor wildlife . J Wildl Manag . 2011 ; 75 : 423 ± 431 . https://doi.org/10.1002/jwmg.63 51. Wilson KR , Anderson DR . Evaluation of two density estimators of small mammal population size . J Mammal . 1985 ; 66 : 13 . https://doi.org/10.2307/1380951 52. Efford MG , Fewster RM . Estimating population size by spatially explicit capture±recapture . Oikos . 2013 ; 122 : 918 ± 928 . 53. Efford MG . secr: Spatially explicit capture-recapture models [Internet]. 2013 . Available: http://CRAN.Rproject.org/package=secr 54. R Core Team. R: A language and environment for statistical computing . [Internet]. Vienna, Austria: R Foundation for Statistical Computing; 2013 . Available: http://www.R-project.org/ 55. Gerber BD , Parmenter RR . Spatial capture-recapture model performance with known small-mammal densities . Ecol Appl . 2015 ; 25 : 695 ± 705 . PMID: 26214915 56. ESRI. ArcGIS Desktop: 10 . Redlands, CA: Environmental Systems Research Institute; 2011 . 57. Kruskal JB . Nonmetric multidimensional scaling: a numerical method . Psychometrika . 1964 ; 29 : 115 ± 129 . 58. Rosenzweig ML , Abramsky Z , Kotler B . Can interaction coefficients be determined from census data? Oecologia . 1985 ; 66 : 194 ± 198 . https://doi.org/10.1007/BF00379854 PMID: 28311589 59. Shumway RH . Time series analysis and its applications: with R examples [Internet] . 3rd ed. New York: Springer; 2011 . Available: http://libproxy.uwyo.edu/login/?url=http://link.springer.com/10.1007/978-1- 4419 -7865-3 60. Rosing MN , Ben-David M , Barry RP . Analysis of Stable Isotope Data: A K Nearest-Neighbors Randomization Test. J Wildl Manag . 1998 ; 62 : 380 ± 388 . https://doi.org/10.2307/3802302 61. Erhardt EB , Wolf BO , Ben-David M , Bedrick EJ . Stable isotope sourcing using sampling . Open J Ecol . 2014 ; 04 : 289 ± 298 . https://doi.org/10.4236/oje. 2014 .46027 62. DeNiro MJ , Epstein S. Influence of diet on the distribution of nitrogen isotopes in animals . Geochim Cosmochim Acta . 1981 ; 45 : 341 ± 351 . https://doi.org/10.1016/ 0016 - 7037 ( 81 ) 90244 - 1 63. Newsome SD , Wolf N , Peters J , Fogel ML . Amino acid δ13C analysis shows flexibility in the routing of dietary protein and lipids to the tissue of an omnivore . Integr Comp Biol . 2014 ; 54 : 890 ± 902 . https://doi. org/10.1093/icb/icu106 PMID: 25104856 64. Ben-David M , Hanley TA , Klein DR , Schell DM . Seasonal changes in diets of coastal and riverine mink: the role of spawning Pacific salmon . Can J Zool . 1997 ; 75 : 803 ± 811 . https://doi.org/10.1139/z97- 102 65. Eckrich CA , Albeke SE , Flaherty EA , Bowyer RT , Ben-David M. rKIN: Kernel-based method for estimating isotopic niche size and overlap . J Anim Ecol. In review; 66. Albeke , S. E. rKIN: (Kernel) Isotope Niche Estimator FunctionsÐrKIN . [Internet]. 2016 . Available: http:// github.com/salbeke/rKIN 67. Nagy KA . Field metabolic rate and food requirement scaling in mammals and birds . Ecol Monogr . 1987 ; 112 ± 128 . 68. Ochocińska D , Taylor JRE . Living at the Physiological Limits: Field and Maximum Metabolic Rates of the Common Shrew (Sorex araneus) . Physiol Biochem Zool Ecol Evol Approaches . 2005 ; 78 : 808 ± 818 . https://doi.org/10.1086/431190 PMID: 16096983 69. Degen AA . Ecophysiology of Small Desert Mammals [Internet]. Berlin, Heidelberg: Springer; 1997 . Available: http://link.springer.com/10.1007/978-3- 642 -60351-8 70. Bradley RD , Durish ND , Rogers DS , Miller JR , Engstrom MD , Kilpatrick CW . Toward a molecular phylogeny for Peromyscus: Evidence from mitochondrial Cytochrome-b sequences . J Mammal . 2007 ; 88 : 1146 ± 1159 . https://doi.org/10.1644/06-MAMM-A -342R.1 PMID: 19924266 71. Krebs CJ , Cowcill K , Boonstra R , Kenney AJ . Do changes in berry crops drive population fluctuations in small rodents in the southwestern Yukon? J Mammal. 2010 ; 91 : 500 ± 509 . 72. Arai S , Bennett SN , Sumibcay L , Cook JA , Song J-W , Hope A , et al. Phylogenetically distinct hantaviruses in the masked shrew (Sorex cinereus) and dusky shrew (Sorex monticolus) in the United States . Am J Trop Med Hyg . 2008 ; 78 : 348 ± 351 . PMID: 18256444 73. Greiman SE , Tkach VV , Cook JA . Description and molecular differentiation of a new Staphylocystoides (Cyclophyllidea: Hymenolepididae) from the dusky shrew (Sorex monticolus) in Southeast Alaska . J Parasitol . 2013 ; 99 : 1045 ± 1049 . https://doi.org/10.1645/ 13 - 302 .1 PMID: 23919726 74. Thibault KM , White EP , Ernest SKM . Temporal dynamics in the structure and composition of a desert rodent community . Ecology . 2004 ; 85 : 2649 ± 2655 . https://doi.org/10. 1890 /04- 0321 75. Aars J , Ims RA . Intrinsic and climatic determinants of population demography: the winter dynamics of tundra voles . Ecology . 2002 ; 83 : 3449 ± 3456 . https://doi.org/10. 1890 /0012- 9658 ( 2002 ) 083 [3449: IACDOP] 2.0 .CO; 2 76. Duchesne D , Gauthier G , Berteaux D. Habitat selection, reproduction and predation of wintering lemmings in the Arctic . Oecologia . 2011 ; 167 : 967 ± 980 . https://doi.org/10.1007/s00442-011 -2045-6 PMID: 21701915 77. Mayes KM , Eckrich CA , Flaherty EA , Ben-David M . Spatial segregation of martens and ermines in Southeast Alaska depends on small mammal availability . J Mammal. In prep; 78. Ben-David M , Flynn RW , Schell DM . Annual and seasonal changes in diets of martens: evidence from stable isotope analysis . Oecologia . 1997 ; 111 : 280 ± 291 . https://doi.org/10.1007/s004420050236 PMID: 28308005 79. Flynn RW , Schumacher TV . Temporal changes in population dynamics of American martens . J Wildl Manag . 2009 ; 73 : 1269 ± 1281 . https://doi.org/10.2193/2008- 169 80. Carey AB , Wilson SM . Induced spatial heterogeneity in forest canopies: responses of small mammals . J Wildl Manag . 2001 ; 1014 ± 1027 . 81. Suzuki N , Hayes JP . Effects of thinning on small mammals in Oregon coastal forests . J Wildl Manag . 2003 ; 352 ± 371 .


This is a preview of a remote PDF: http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0189471&type=printable

Carolyn A. Eckrich, Elizabeth A. Flaherty, Merav Ben-David. Functional and numerical responses of shrews to competition vary with mouse density, PLOS ONE, 2018, DOI: 10.1371/journal.pone.0189471