Smallmammals Yosemiti Sci2008
Smallmammals Yosemiti Sci2008
Smallmammals Yosemiti Sci2008
If you wish to distribute this article to others, you can order high-quality copies for your
colleagues, clients, or customers by clicking here.
Permission to republish or repurpose articles or portions of articles can be obtained by
following the guidelines here.
Updated information and services, including high-resolution figures, can be found in the online
version of this article at:
http://www.sciencemag.org/cgi/content/full/322/5899/261
Supporting Online Material can be found at:
http://www.sciencemag.org/cgi/content/full/322/5899/261/DC1
A list of selected additional articles on the Science Web sites related to this article can be
found at:
http://www.sciencemag.org/cgi/content/full/322/5899/261#related-content
This article cites 24 articles, 2 of which can be accessed for free:
http://www.sciencemag.org/cgi/content/full/322/5899/261#otherarticles
This article has been cited by 7 article(s) on the ISI Web of Science.
This article has been cited by 8 articles hosted by HighWire Press; see:
http://www.sciencemag.org/cgi/content/full/322/5899/261#otherarticles
This article appears in the following subject collections:
Ecology
http://www.sciencemag.org/cgi/collection/ecology
Science (print ISSN 0036-8075; online ISSN 1095-9203) is published weekly, except the last week in December, by the
American Association for the Advancement of Science, 1200 New York Avenue NW, Washington, DC 20005. Copyright
2008 by the American Association for the Advancement of Science; all rights reserved. The title Science is a
registered trademark of AAAS.
The following resources related to this article are available online at www.sciencemag.org
(this information is current as of May 5, 2010 ):
REPORTS
A
1
documented (710). Models of future climatechange scenarios predict large range shifts, high
global extinction rates, and reorganized communities (11, 12), but model outcomes are also highly uncertain (13, 14). Most studies of species
responses span only a few decadestypically
from the 1960 or 1970s, which was a relatively
cool period, to the present. Such results can be
confounded by decadal-scale climate oscillations
(15) and landscape modification (8, 16). Furthermore, range shifts are uncertain when confounded
by false absences due to limited historic sampling
and inability to control for changes in detectability
between sampling periods (17, 18).
www.sciencemag.org
SCIENCE
VOL 322
10 OCTOBER 2008
261
Future warming is predicted to cause substantial turnover of species within North American
National Parks, including Yosemite (25). Given
marked regional warming over the past century,
we predicted that species ranges should have
shifted upward (5, 10). This should manifest as
upward contraction of the lower range limit for
mid- to high-elevation species, upward shift of
the entire range or expansion of the upper limit
for low- to mid-elevation species, and altered community composition within elevational bands (9).
Elevational ranges of species and their habitats differed markedly between the gradual western and steep eastern slopes of the transect (19)
(Fig. 1). On the west slope, we trapped small
mammals at 121 sites compared to 56 in Grinnells
time (table S1), but overall effort and elevational
range (~50 to 3300 m) were comparable (22).
There were fewer sites on the east side in both
time periods (9 for Grinnell, 12 for resurveys)
because of limited extent (Fig. 1). Our analyses
of richness and turnover focused on species detectable by standardized trapping (37 species) or
by observation (6 species; table S2). To test for
elevational shifts, we applied occupancy modeling (22, 23) to the 23 west slope taxa with sufficient trapping records to estimate detectability
in both periods (tables S1 and S2 and Fig. 2). The
best detection model in a set of 36 (table S3) was
used to calculate the probability of a false absence (Pfa) across trapping sites, where a species
was not observed in one sampling period but was
in the other (Table 1). Range shifts were significant if Pfa 0.05. For each species we evaluated
eight hypothesized relationships of occupancy,
era, and elevation (fig. S2) using the 14 best detection models (table S3) to create model-averaged
occupancy-elevation profiles (Fig. 2 and fig. S3).
Conservatively, we excluded shifts that were statistically significant but biologically trivial (Fig.
3). In most cases where the Pfa test indicated an
elevation shift, occupancy models agreed (Table
1 and fig. S3). Exceptions occurred when occupancy models were weak (i.e., insufficient data)
or detected changes in occupancy at elevations
other than range limits, or when nonstandard data
(i.e., records from ad hoc collecting) were included in Pfa tests but not in occupancy models.
Elevation limits shifted mostly upward (Table
1 and Fig. 3A), and this occurred more frequently
for lower than upper limits (c2 = 4.26, df = 1, P =
0.039). Twelve of 28 (43%) west slope species
showed significant shifts in lower limits, of which
10 increased (mean = +475 m) and two, both
shrews, decreased (mean = 744 m). In contrast,
upper limits changed significantly in only seven
instances, with similar numbers of upward (n = 4,
mean = +501 m) and downward shifts (n = 3,
mean = 309 m).
High-elevation species typically experienced
range contractions, whereas low-elevation species
expanded their ranges upward (c2 = 8.8, df = 2,
P = 0.012), a pattern expected with increased
temperature. Lower range limits contracted in
50% of the high-elevation species but in only
262
Fig. 1. Map of surveyed sites in Grinnell (Historic) and Current surveys relative to the Yosemite
National Park boundary and life zones (upper panel), and to an averaged elevational profile (lower
panel).
10 OCTOBER 2008
VOL 322
SCIENCE
www.sciencemag.org
REPORTS
REPORTS
survey, we provide an unbiased comparison of
changes in species ranges at the centennial scale.
Because much of the transect spans a longprotected National Park, confounding effects of
land-use change are minimized. Even so, vegetation has changed within YNP over this period,
in part due to fire suppression (22). The park was
hardly pristine in the early 20th century, with
ranching of introduced herbivores in Yosemite
Valley and the high country recovering from historical overgrazing. As examples, expansion by
C. californicus and west slope P. truei are associated with fire-related conversion of conifer to
shrub habitats, whereas the downward shift of
S. monticola could reflect recovery of their
preferred wet meadow habitats. Increased prevalence of mesic small mammals following cessation of grazing has also been reported for an
Table 1. Analyses of elevation change for 28 west slope species. Given are
average detectability per site for Grinnell [P(G)] and current [P(C)] periods, original
elevation range, changes in upper (U) and lower (L) range limit that are significant
by the Pfa tests, the best supported form of the occupancy model (Elev, elevation;
NA, not analyzed), the cumulative Akaikes Information Criterion weight for all
No.
1
2
Species
P(G)
P(C)
Original
elevation
range (m)
0.81
0.58
571160
models with those terms (w), and original Lifezone classification (18), where L and H
refer, respectively, to species with mostly low- to mid-elevation ranges (<2000 m)
and mid- to high-elevation ranges (>2000 m) in Grinnells time; P. maniculatus
covered the entire transect. Values in bold are further supported by occupancy
models. See fig. S4 for elevation plots and models of individual species.
Best
occupancy
model
Range limit
change (m)
+112 U
+589 U, +468 L
0.50
0.99
0.32
0.74
0.37
0.48
0.74
0.53
0.48
0.78
0.39
0.83
0.32
0.56
NA
0.72
NA
0.89
0.62
0.60
0.40
0.36
0.87
0.93
571160
1831220
5
6
0.28
0.32
0.99
0.19
0.93
0.97
193914
549914
22123287
7
8
9
10
11
12
13
14
15
16
Dipodomys heermanni
Microtus longicaudus
Zapus princeps
Tamias senex
Spermophilus lateralis
Sorex palustris
Neotoma cinerea*
Spermophilus beldingi*
Tamias alpinus
Ochotona princeps
0.16
0.99
0.98
0.95
0.70
0.39
0.90
0.98
0.92
NA
0.98
0.98
0.90
0.71
0.89
0.23
0.71
0.98
0.95
NA
571025
6233287
12913185
14022743
16463200
16583155
17983287
22863287
23073353
23773871
17
Peromyscus
maniculatus*
Thomomys bottae
Spermophilus beecheyi
Neotoma macrotis
Peromyscus boylii
Sorex trowbridgii
Microtus montanus*
Tamiasciurus
douglasi*
Tamias
quadrimaculatus
Tamias speciosus*
Thomomys monticola
Marmota flaviventris
0.99
NA
0.50
0.90
0.98
0.71
0.81
0.99
NA
0.82
0.91
0.97
0.88
0.98
573287
571676
612734
1831646
1832469
11602286
12173155
No change
No change
250 U
+67 U
122 L
No change
No change
Era*(Elev + Elev2)
NA
Era*(Elev + Elev2)
Elev + Elev2
Elev + Elev2
Elev + Elev2
Elev + Elev2
NA
NA
12293185
No change
NA
0.95
1.00
NA
NA
0.85
1.00
NA
NA
14942210
17683155
19053155
24693353
+50 U
+128 L, +65 U
No change
No change
Era*(Elev + Elev2)
Era*(Elev + Elev2)
NA
NA
25
26
27
28
* Similar trends are observed for east-side populations (see fig. S4).
0.36
0.99
0.99
18
19
20
21
22
23
24
Range expansions
+505 U
Elev
Microtus californicus
Reithrodontomys
megalotis
Peromyscus truei*
Chaetodippus
californicus
Sorex ornatus
Sorex monticolus
3
4
Elev
Era*(Elev + Elev2)
+800 U
Era*(Elev + Elev2)
485 L
Era *(Elev + Elev2)
1003 L
Era + Elev + Elev2
Range contractions
+63 L, 293 U
Era*Elev
+614 L
Era + Elev + Elev2
+159 L, 64 U
Era + Elev + Elev2
+1007 L, 334 U
Elev +Elev2
+244 L
Era*(Elev + Elev2)
+512 L
Era + Elev + Elev2
+609 L, 719 U
Era*(Elev + Elev2)
+355 L
Elev
+629 L
Era + Elev
+153 L
NA
No change
NA
0.78
1.00
NA
NA
TransitionHudsonian (H)
TransitionCanadian (H)
CanadianHudsonian (H)
CanadianHudsonian (H)
CanadianArctic-Alpine (H)
These species were encountered by observation and/or specialized trapping and were not subject to occupancy analyses.
www.sciencemag.org
SCIENCE
VOL 322
10 OCTOBER 2008
263
REPORTS
Fig. 2. Example elevation plots
from the west slope transect of upward range expansion (T. alpinus
and P. truei) (A and C), and
range collapse (N. cinerea) (B).
Shown are occupied (black) and
unoccupied (gray) sites, probability of false absence (Pfa),
and model-averaged occupancyelevation profiles (table S3 and
fig. S2). P. truei colonized high
elevations west of the Sierra crest
from the eastern slope. Red
marks for historical elevation
profile of T. alpinus refer to ad
hoc records.
264
10 OCTOBER 2008
VOL 322
SCIENCE
www.sciencemag.org