Introduction

Non-migratory butterfly species endemic to montane meadows are particularly susceptible to the effects of habitat loss and climate change because they are: (1) dependent upon specialized plant-insect interactions that are at risk of phenological mismatching (Parmesan 2007; Singer and Parmesan 2010), (2) constrained to potentially shrinking habitats due to forest encroachment (Roland et al. 2000; Matter et al. 2004; Roland and Matter 2007), and (3) vulnerable to genetic isolation as a result of small, isolated populations (Keyghobadi et al. 1999, 2005). Climate scenarios predict that high elevations where montane meadow butterfly populations exist will face larger, more rapid changes and are likely to have the greatest rate of species loss compared to other ecosystems (Kim et al. 2002; Thuiller et al. 2005). Due to the combined and potentially synergistic threats of habitat loss and climate change threatening butterflies worldwide (Warren et al. 2001; McLaughlin et al. 2002; Forister et al. 2010), it is essential to determine occupancy and identify habitat requirements for sensitive montane meadow butterfly species in order to prioritize conservation efforts.

Our study focused on the Clodius Parnassian butterfly (Parnassius clodius; Lepidoptera: Papilionidae) in Grand Teton National Park (GRTE) and the surrounding Bridger-Teton National Forest (BTNF) and John D. Rockefeller, Jr. Memorial Parkway (JODR) in the Rocky Mountains of Wyoming, United States (hereafter Teton region). This species is a range-restricted, high elevation montane meadow butterfly species for which habitat alterations and climate change may cause increasing conservation concerns. Although P. clodius is not currently threatened, it faces similar issues to those of the related European Apollo butterfly (Parnassius apollo), which has been declining throughout the twentieth century due to long-term climatic changes, habitat succession, anthropogenic factors, genetic erosion, and behavioral changes in small demes (Nakonieczny et al. 2007). P. apollo is categorized by the International Union for the Conservation of Nature as Near Threatened (van Swaay et al. 2011) and is considered a high priority for conservation (Todisco et al. 2010). Parnassius species, known as the “invertebrate epitome of the conservation of mountain habitats” in Europe, are considered a valuable study organism and can serve as a flagship species for montane environments (Todisco et al. 2010). By studying the habitat preferences of P. clodius, we may be able to anticipate similar threats P. apollo faces and predict how these threats may impact P. clodius in the future.

Prior to this research, there was limited information on the distribution and abundance of Parnassius populations in the Teton region outside of one well-studied population of P. clodius in the meadow along Pilgrim Creek road (1500 × 300 m in size) within GRTE (Auckland et al. 2004). Mark-recapture studies were conducted on the Pilgrim Creek population, considered to be one of the largest populations of P. clodius in the Teton region, to assess population parameters including sex ratio, population size, percentage of mated females, and emergence dates for males and females (Auckland et al. 2004). Population fluctuations documented by these mark-recapture studies from 1998 to 2000 (Auckland et al. 2004) and 2009–2012 (Jill A. Sherwood and Diane M. Debinski, unpublished data) indicated the need for a better understanding of the ecosystem-wide distribution pattern of P. clodius. Given there are no extensive data on P. clodius outside of Pilgrim Creek, it is unclear whether the current population configuration across the Teton region is a stable one.

To document site occupancy patterns of P. clodius on an ecosystem-wide level and to assess how particular components of habitat influenced patterns of occupancy, we conducted a series of presence-absence butterfly surveys and vegetation surveys in potentially suitable habitat for P. clodius in the Teton region in 2012 and 2013. There were two objectives of this study: (1) use occupancy estimates and detection probabilities of P. clodius to determine where current populations exist in the Teton region and (2) analyze potentially suitable Parnassius butterfly habitat to examine fine-scale habitat preferences of P. clodius. For the second objective, we hypothesized that the occupancy of P. clodius would be positively correlated with abundance of low sagebrush (Artemisia arbuscula) and their preferred nectar source, sulphur flower buckwheat (Eriogonum umbellatum). Field observations indicated that A. arbuscula is generally associated with a higher cover of bare ground, which allows for the growth of their host plant, longhorn steershead (Dicentra uniflora), and important nectaring forbs.

Methods

Species description

Parnassius clodius is a medium-sized white butterfly that inhabits montane meadows. The species occurs along the northern Rocky Mountains in northwestern United States and western Canada (Scott 1986). Parnassius species are univoltine butterflies that overwinter as eggs, pupate on the soil during spring, and have an adult flight season from mid June to mid July (Scott 1986). The host plant of P. clodius is D. uniflora, a spring ephemeral wildflower. The primary nectar source of P. clodius in the Teton region is E. umbellatum (Auckland et al. 2004). Data collected during 2011–2013 documented that several additional nectar sources are being used by P. clodius in the Teton region, including arrow-leaved balsamroot (Balsamorhiza sagittata), groundsel (Senecio spp.), sticky geranium (Geranium viscosissimum), nettle-leaved giant-hyssop (Agastache urticifolia), spreading dogbane (Apocynum androsaemifolium), several yellow and purple flowering species within the Asteraceae family, snowbrush ceanothus (Ceanothus velutinus), birch-leaved spiraea (Spiraea betulifolia), and common snowberry (Symphoricarpos albus; Kimberly E. Szcodronski and Jill A. Sherwood, unpublished data).

Study area

The Teton region is a relatively pristine, large-scale protected area within the Greater Yellowstone Ecosystem (GYE) that contains a heterogeneous distribution of habitat types along elevation and hydrological gradients that include conifer and aspen forest (Populus spp.), willow shrubland (Salix spp.), sagebrush flats (Artemisia spp.), and montane meadows. Landsat satellite multispectral imagery data were used to classify six montane meadow types in the GYE along a hydrological gradient ranging from hydric willow and sedge (Carex spp.) meadows to mesic meadows with high forb cover to xeric sagebrush meadows (Jakubauskas et al. 1996, 1998; Kindscher et al. 1998; Debinski et al. 1999). From 1997 to 2007, Debinski et al. (2006, 2010, 2013) collected vegetation and butterfly distribution data across this hydrological gradient within 55 montane meadows in the GYE, including 25 meadows in the Teton region. Butterfly abundance counts from these surveys indicated that Parnassius butterflies were most often found in mesic and xeric meadows (Debinski et al. 2013). Furthermore, Auckland et al. (2004) determined P. clodius was located at highest densities in the Teton region in “dry, cobbly, sagebrush meadows where its host plant species, D. uniflora, is abundant.”

Based on the butterfly’s known habitat preferences, study sites for this project were restricted to meadows with montane mesic herbaceous vegetation, montane xeric herbaceous vegetation, and Artemisia spp. cover including big sagebrush (Artemisia tridentata) and A. arbuscula. The 2002–2005 GRTE Vegetation Mapping Project (Cogan et al. 2005) provided geospatial vegetation data that were used in ArcGIS 10.3.1 (ESRI 2015) to identify 61 meadows as potentially suitable Parnassius habitat in the Teton region. Of the 61 sites, 41 sites were surveyed as part of the study in 2013, three sites were not classified as potential habitat when ground-truthed and were removed from the study, nine sites were excluded due to inaccessibility (> 4.5 km from an access road), and the flight season ended before we were able to survey the remaining eight sites. Nineteen of the 41 sites were also surveyed in 2012 during a pilot study. The selected study sites were distributed across ~ 135,072 ha (Supplemental Material, Table A.1; Table A.2). Study sites ranged in size from 1 to 1307 ha and in elevation from 2006–2491 m (Supplemental Material, Table A.1; Table A.2).

Butterfly surveys

Presence–absence butterfly surveys were conducted for P. clodius in 19 study sites in the summer of 2012 (Supplemental Material, Table A.1) and 41 study sites (including all of the 2012 sites) during the summer of 2013 (Supplemental Material, Table A.2). In our analysis “presence” is defined as occupancy and “absence” is defined as both true absence and non-detection. While the 2012 surveys were not included in all components of our final analyses, these data are useful in providing evidence that the results of our 2013 surveys were representative of occupancy patterns across years. Surveys occurred from mid June to late July. Based on delayed vegetation phenology and later butterfly emergence in higher elevation sites, we suspected that P. clodius would emerge at higher elevation sites 1–2 weeks later than lower elevation sites. Therefore, we conducted surveys at lower elevation sites first and moved up along the elevation gradient as the season progressed. The emergence date of P. clodius varies annually; therefore it was important to verify that the species had emerged before surveys began. The Pilgrim Creek population, one of the largest P. clodius populations in the Teton region and among the lower elevation sites, was monitored frequently near the predicted Parnassius adult flight season to document butterfly emergence and determine the start of the survey period.

To account for imperfect detection, presence-absence surveys were conducted at study sites in 2012 and 2013 throughout P. clodius’ flight season by two independent observers per visit (MacKenzie et al. 2002, 2006; MacKenzie and Royle 2005). Each survey by an observer was considered an independent detection occasion. In 2012, we surveyed 19 sites with six sites visited twice (four detection occasions/site) and 13 sites visited once (two detection occasions/site). In 2013, we surveyed 41 sites with 39 sites visited twice (four detection occasions/site), one site visited twice with three detection occasions (one visit with one surveyor and one visit with two surveyors), and one site surveyed once by one observer (one detection occasion). The detection pattern obtained across multiple visits within a relatively short time was used to estimate detection probability and to correct for false absences (MacKenzie et al. 2002, 2006). Each observer searched for P. clodius for 30 min per survey. The total number of P. clodius observed at a site was summed across observers and surveys at each location for 2013. A site was considered occupied if at least one P. clodius was verified to be present in the meadow in at least one of the surveys within a year. Detection may have been lower in larger meadows (> 250 ha) due to their size, but additional site visits outside of our surveys support our findings of no occupancy.

The other Parnassius species in the Teton region, Rocky Mountain Parnassian (Parnassius smintheus), has similar wing coloration of P. clodius and is present in low numbers in xeric meadows throughout the region (Debinski et al. 2013). Due to the high similarity in phenotype, each Parnassius butterfly observed was caught with an insect net during surveys to verify species identity and was released unharmed at the end of the survey period. Based on high P. clodius recapture rates (\(\overline{\text{x}}\) = 26%) from Pilgrim Creek mark-recapture studies from 1998 to 2000 (Auckland et al. 2004) and 2009–2012 (Jill A. Sherwood and Diane M. Debinski, unpublished data), careful handling by trained technicians does not have a negative effect on the survival of this species. Butterfly surveys were performed during optimal butterfly flight conditions which occurred mid June–late July between 10:00 and 17:00 h when the temperature was above 21°C, sun was not obscured by clouds, and wind speed was < 16 km/h (Auckland et al. 2004).

Vegetation surveys

We surveyed vegetation at each study site once during the Parnassius butterfly flight season in the summer of 2013, ranging from mid June–late July depending on the elevation of the study site. Study sites varied considerably in size and were therefore classified into three size classes, with the number of vegetation transects proportional to meadow size (< 10 ha = 1 transect, 10–100 ha = 2 transects, > 100 ha = 3 transects). Observations indicate vegetation composition is relatively homogenous within vegetation classifications in our study area, which allowed a small sampled area to be representative of an entire meadow. Vegetation data were collected using frequency of occurrence sampling methods (Bonham 2013). We recorded the presence or absence of vegetation categories along a 100 m transect in 1 × 0.5 m quadrat placed every 5 m (n = 20 quadrats per transect). Frequency of occurrence was calculated as the total number of times a vegetation category occurred in a transect divided by the total number of quadrats surveyed. If multiple transects were collected in a meadow, frequency of occurrence was calculated by averaging across transects. Data were recorded for 47 vegetation categories in the field and were condensed into 14 major vegetation categories based on taxonomic similarities and relevance to P. clodius to simplify model selection and inference (Supplemental Material, Table A.3). Summary statistics for the major vegetation categories were performed in program R version 3.4.0 (R Core Team 2017) using the lsmeans package (Lenth 2013) to compare sites with P. clodius present and absent.

Based on the differences in mean frequency of occurrence between sites with and without P. clodius, we further reduced the 14 vegetation categories to seven site covariates to be used in occupancy modeling. The seven covariates included the two dominant sagebrush species, (1) A. tridentata and (2) A. arbuscula as well as (3) dominant forb, Lupinus spp., (4) dominant forb and nectar source, B. sagittata, (5) P. clodius’ main nectaring source, E. umbellatum, (6) P. smintheus’ host plant, Sedum spp., and (7) secondary nectar sources, which encompassed all nectar sources utilized less frequently by P. clodius (refer to Supplemental Material, Table A.3 for list of nectar sources). It is important to note that a key habitat variable, the host plant of P. clodius, D. uniflora, was not included as a covariate due to sampling logistical constraints. D. uniflora is briefly visible in early spring after snowmelt and senesces before the butterfly flight season. We included the host plant of P. smintheus, Sedum spp., within the selected covariates to check for any possible association between Sedum spp. and P. clodius. Pearson’s product-moment correlation coefficients (r) were calculated in program R for the seven site covariates using the “cor” function and the “cor.test” function was used to test for significance of the correlations (p < 0.05). We did not include any combination of covariates in the same model with |r| > 0.70.

Occupancy modeling analyses

We used program PRESENCE version 11.8 (Hines 2006) to analyze occupancy patterns for P. clodius by estimating detection probabilities (p) and the probability of a site being occupied (ψ). Single-species, single-season models (MacKenzie et al. 2002, 2006) were used for the analyses. To examine sampling and site covariate effects on detection probability and occupancy, a priori models were developed and evaluated using an information-theoretic approach (Burnham and Anderson 2002).

In order to quantify the consistency in P. clodius occupancy patterns across years, we tested occupancy models for a year effect with the combined 2012–2013 data. The model that allowed occupancy to vary by year [ψ(year), p(.)] had a ∆AIC c value of 2.01 compared to the top model that kept occupancy constant across years [ψ(.), p(.)]. Furthermore, the beta coefficient for year in the logistic model overlapped one (β = 1.13, 95% CI 0.71–1.54), indicating a non-informative parameter (Arnold 2010). While the 2012 data were essential for addressing between-year variation in occupancy and detection probabilities, the data from 2012 were insufficient for further occupancy analyses in relation to habitat variables and therefore the remainder of this paper focuses on the results of the 2013 surveys.

A two-stage approach was used to analyze occupancy patterns. First, models were compared by evaluating the effects of sampling and site covariates on detection probability while holding occupancy constant (MacKenzie et al. 2002, 2006). Six models were chosen to evaluate effects on detection probability: meadow size [ψ(.), p(meadow size)], elevation [ψ(.), p(elevation)], observer [ψ(.), p(observer)], survey [ψ(.), p(survey)], fully time-varying model where there are different detection probabilities for each detection occasion [ψ(.), p(t)], and null constant parameter model where p and ψ are held constant as a comparison [ψ(.), p(.)]. We were careful to control for weather conditions during butterfly surveys; therefore factors of wind speed, temperature, and time of day were excluded from model comparisons. Akaike’s Information Criterion with small sample bias correction (AIC c ) and Akaike’s weight (wi) were used to rank the detection probability models and the best model was selected based on these criteria (Burnham and Anderson 2002). The probability of occupancy given that the species is not detected at a site was calculated by: [Pr(species present and not detected)]/[Pr(species not detected)] (MacKenzie et al. 2006).

Next, 27 additive covariate models were created to evaluate the influence of the seven site covariates on occupancy probability. We followed the recommendations of MacKenzie et al. (20022006) by using the estimated detection probability value to correct for false absences and estimate occupancy as a function of detection probability. Thus, for each occupancy model, the best-supported detection model was used to account for detection probability in all of the occupancy models. We interpreted the top-supported occupancy models (ΔAIC c  ≤ 2; Burnham and Anderson 2002) with vegetation variables as covariates as the most important factors influencing the habitat suitability of P. clodius. Model-averaged parameter estimates were used to plot the relationship of P. clodius with A. tridentata using the RMark package (Laake 2013) in program R.

Results

P. clodius was present at 12 (63%) out of 19 sites surveyed in 2012 (Supplemental Material, Table A.1; Table A.4) and 32 (78%) of the 41 sites surveyed in 2013 (Fig. 1; Supplemental Material, Table A.2). Of the sites surveyed, P. clodius was found across the entire elevation gradient from 2006 to 2491 m (Supplemental Material, Table A.1; Table A.2). In 2013, 21 populations had low total observations of P. clodius during the occupancy surveys, with less than nine individuals summed across the four presence-absence surveys (\(\overline{\text{x}}\) = 3 butterflies observed/site; Supplemental Material, Table A.5). The remaining 11 populations had total observations ranging from 14 to 35 individuals across the presence-absence surveys (\(\overline{\text{x}}\) = 25 butterflies observed/site; Supplemental Material, Table A.5).

Fig. 1
figure 1

Vegetation map of Grand Teton National Park, WY and surrounding territories displaying occupancy of Parnassius clodius in meadow sites (n = 41) surveyed from mid June–late July 2013. Sites were considered occupied if at least one butterfly was verified to be present in at least one of the presence-absence butterfly surveys; sites were considered to be unoccupied if no butterflies were observed across all presence-absence surveys

Based on the null constant parameter model, P. clodius had high occupancy in 2013 of ψ = 0.78 ± 0.07 SE (95% CI 0.63–0.88) with a detection probability of p = 0.75 ± 0.04 SE (95% CI 0.67–0.82). The top-ranked detection probability model for P. clodius was [ψ(.), p(meadow size)] (Table 1); so this model was used for modeling occupancy. Meadow size was a non-informative parameter for the detection of P. clodius (β = 0.01 ± 0.01 SE; Arnold 2010). The probability of occupancy given that the species was not detected at a site for P. clodius was 0.01, indicating a low probability of false negatives in our survey data.

Table 1 Summary of model selection statistics examining variables associated with probability of detection (p) of Parnassius clodius in Grand Teton National Park and Bridger-Teton National Forest, WY from mid June–late July 2013

Of the 27 occupancy models assessed for P. clodius, one model had strong support (ΔAIC c  ≤ 2, Table 2). The model with the highest level of support showed that P. clodius was negatively correlated with the frequency of occurrence of A. tridentata (β = − 21.39 ± 21.10 SE) and Lupinus spp. (β = − 20.03 ± 21.24 SE). The model with strong support (ΔAIC c  ≤ 2) accounted for 84% of the AIC c model weight out of all models considered (Table 2). The nine models containing A. tridentata as a covariate were the top nine ranked models out of the 27 models assessed (Table 2). A. tridentata and Lupinus spp. are the only variables with evidence of difference between sites where P. clodius was present and sites where P. clodius was absent, both with strong evidence of a difference (p < 0.0001; Fig. 2; Supplemental Material, Table A.6). The occurrence of P. clodius decreased as the frequency of occurrence of A. tridentata increased (Fig. 3).

Table 2 Summary of model selection statistics examining variables associated with the probability of occupancy (ψ) of Parnassius clodius in Grand Teton National Park and Bridger-Teton National Forest, WY from mid June – late July 2013
Fig. 2
figure 2

Average frequency of occurrence for vegetation variables in sites where Parnassius clodius was present (n = 32) and sites where Parnassius clodius was absent (n = 9) in Grand Teton National Park and Bridger-Teton National Forest, WY during the summer of 2013. Error bars indicate standard error for each vegetation category and a star above a bar indicates a significant p-value (p < 0.05)

Fig. 3
figure 3

Influence of the frequency of occurrence of Artemisia tridentata on the probability of Parnassius clodius occupancy in Grand Teton National Park and Bridger-Teton National Forest, WY, during the summer of 2013. The shaded area depicts the 95% confidence interval

Correlation coefficients for the seven site covariates are provided in the Supplemental Material, Table A.7. Interestingly, there was no significant negative correlation between the two sagebrush species, A. tridentata and A. arbuscula (r = − 0.28, p = 0.08). A. arbuscula showed a significant positive correlation with E. umbellatum (r = 0.41, p = 0.01) and Sedum spp. (r = 0.48, p = 0.001). A. tridentata was significantly positively correlated with Lupinus spp. (r = 0.38, p = 0.01).

Discussion

P. clodius was more broadly distributed across the Teton region than expected based on previous studies in the area (Debinski et al. 2006, 2013), occupying 78% of the meadows we predicted would be suitable in 2013. These results are encouraging in the context of long-term population viability. However, it is also important to consider the total habitat available to P. clodius in the Teton region and the population size of P. clodius at these meadows. Using the GRTE geospatial vegetation data (Cogan et al. 2005), we calculated the proportion of area within the GRTE and JODR boundary that could serve as suitable habitat for P. clodius based on previous research in the study area (montane mesic herbaceous vegetation, montane xeric herbaceous vegetation, A. tridentata shrubland, and A. arbuscula shrubland). We found 21.7% of the study area was made up of these four habitat types. However, the results of this study indicate meadows containing A. tridentata were negatively associated with the presence of P. clodius so we removed this habitat type from our calculations. Based on this final classification, 8469 ha (6.3% of the study area) offers suitable habitat for P. clodius in the Teton region.

Furthermore, 21 out of the 32 locations where P. clodius was found in 2013 had low numbers of individuals observed across the four surveys (Supplemental Material, Table A.5). Mark-recapture studies would need to be conducted to obtain estimates of actual population sizes for P. clodius, but our data suggest that few populations within this Teton region are large (i.e., maximum number observed > 25 butterflies summed across all four surveys).

The most supported occupancy model (AIC c  = 147.42, wi = 0.84; Burnham and Anderson 2002) in our analysis allowed us to understand factors influencing habitat suitability of P. clodius. We expected that the occupancy of P. clodius would be positively correlated with the abundance of A. arbuscula and their preferred nectar source, E. umbellatum. Our hypothesis was not supported as occupancy modeling results showed no significant relationships with A. arbuscula or E. umbellatum (Table 2). Instead, our analyses indicated that two dominant constituents of montane meadow plant communities, A. tridentata and Lupinus spp., were negatively associated with the presence of P. clodius (Table 2). A. tridentata meadows cover a much larger percentage of the habitat (15.4%) in GRTE and JODR than A. arbuscula (0.8%). Even though the negative correlation between the two sagebrush species, A. tridentata and A. arbuscula, was not significant (Supplemental Material, Table A.7), our field observations at a larger spatial scale indicate that meadows with high percent cover of A. tridentata often have limited cover of A. arbuscula. We speculate that higher percentage of bare ground associated with A. arbuscula is likely to be correlated with P. clodius’ host plant D. uniflora, but this relationship may not have been detected due to the frequency of occurrence sampling methods used to survey vegetation in our study. Notably, we found a significant positive correlation between A. arbuscula and E. umbellatum (Supplemental Material, Table A.7), the two species that we predicted would influence the occupancy of P. clodius.

It may be possible that the presence of A. tridentata could affect our ability to see P. clodius during occupancy surveys, however we do not believe this was an issue in our study area. Although A. tridentata is considerably taller (0.6–1.8 m; Tilley et al. 2012) than A. arbuscula (0.2–0.4 m; Tilley and St. John 2012), P. clodius is a highly visible butterfly with a flight pattern that is above A. tridentata. Additionally, there is ample space between A. tridentata shrubs such that observation is not limited in the Teton region. The negative association between P. clodius and A. tridentata warrants further investigation to determine whether it is a direct relationship or whether A. tridentata is negatively associated with other important environmental factors that affect habitat conditions for P. clodius. Greater sampling intensity in meadows with high cover of A. tridentata could assist in determining whether A. tridentata is an indicator of low habitat quality for P. clodius or if it reduces detectability.

Our ecosystem-wide surveys have elucidated aspects of population structure that could serve as the foundation for future P. clodius metapopulation research in this study area. The Teton region is a complex landscape in the context of P. clodius habitat because the region is characterized by potentially suitable meadows surrounded by unsuitable habitats, including meadows that lack the proper resources for occupancy and conifer forests, which are considered a barrier to movement for Parnassius butterflies (Roland et al. 2000; Matter et al. 2004; Roland and Matter 2007). A metapopulation study would be important to help investigate whether the current metapopulation dynamics are stable and if P. clodius is at risk due to potential annual weather variability and climate fluctuations. Colonization and extinction rates would be needed to classify the Teton region in the context of a source-sink or core-satellite metapopulation model. Understanding forest edge responses of P. clodius and the connectivity among suitable meadows, as has been done for P. smintheus in Canada (Roland et al. 2000; Matter et al. 2004; Ross et al. 2005), would also aid in assessing the long-term viability of this population.

Our research provides a greater understanding of the distribution of P. clodius across the Teton region and valuable insight into habitat characteristics that may be influencing the occupancy of P. clodius in montane meadows. P. clodius was found to inhabit a small total number of hectares in the Teton region. This information is useful for managers and scientists in the GYE for monitoring species of concern and predicting the effects of habitat changes on regional plant-insect associations and biodiversity. Furthermore, our research is relevant to conservation efforts for other Parnassius species that are threatened or endangered worldwide, such as P. apollo, because Parnassius species have similar life history patterns and habitat requirements globally. Suitable habitat for butterflies is decreasing in different parts of the world (Warren et al. 2001; Roland and Matter 2007) and butterflies that exist in small, isolated populations scattered across a small portion of the landscape, such as P. clodius, could be at a greater risk of decline or extinction than other species. Given the global conservation concerns associated with Parnassius species and the relatively small number of P. clodius individuals observed during many of our occupancy surveys, it would be advisable to assess the population trends and explore the metapopulation dynamics of P. clodius in this ecosystem to better evaluate conservation concerns for this butterfly in the context of potential future habitat changes.