Bat-cave catchment areas: using stable isotopes (deltaD) to determine the probable origins of hibernating bats

The application of stable hydrogen isotope (deltaD) techniques has swiftly advanced our understanding of animal movements, but this progression is dominated by studies of birds and relatively long-distance, north-south migrants. This dominance reflects the challenge of incorporating multiple sources of error into geographic assignments and the nature of spatially explicit deltaD models, which possess greater latitudinal than longitudinal resolution. However, recent progress in likelihood-based assignments that incorporate multiple sources of isotopic error and Bayesian approaches that include additional sources of information may advance finer-scale understanding of animal movements. We develop a stable-isotope method for determining probable origins of bats within hibernacula and show that this method produces spatially explicit, continuous assignments with regional resolution. We outline how these assignments can be used to infer hibernacula connectivity, an application that could inform spatial modeling of white-nose syndrome. Additionally, estimates of seasonal and annual flight distances for many cave-dwelling bat species can be derived from this approach. We also discuss how this application can be used in general to provide insights into variable migratory and foraging strategies within bat populations.     

Efficacy of stable isotope ratios in assigning endangered migrants to breeding and wintering sites.

A primary constraint on effective conservation of migratory animals is our inability to track individuals through their annual cycle. One such animal is the endangered southwestern subspecies of the Willow Flycatcher, which is difficult to distinguish from conspecifics. Identifying wintering regions used by the endangered subspecies would be an important step in formulating an effective conservation strategy. Our objective was to use stable isotope ratios as a means of identifying wintering sites of Southwestern Willow Flycatchers. We analyzed stable isotope ratios of carbon, nitrogen, and hydrogen from feathers of breeding and wintering Willow Flycatchers. Based on winter samples, we document a positive trend in hydrogen isotope ratios across latitude. We also found that Willow Flycatchers use C4 food webs south of 8 degrees N latitude, but we found no evidence of use of C4 food webs farther north. Nitrogen stable isotope ratios of feathers showed no discernable geographic variation. Discriminant function analyses, based on stable isotope ratios of wintering Willow Flycatchers, were only useful (>50% accurate) for assigning individuals to winter regions if the regions were large and the threshold probability for assignment was relatively high. When using these discriminant functions, most breeding samples of Southwestern Willow Flycatchers were assigned to two wintering regions: central Mexico and Ecuador. We think that assignment of Southwestern Willow Flycatchers to Ecuador is unrealistic. Given the large percentages of samples that could not be classified with certainty, we are not confident that these two regions are truly more likely to harbor wintering Southwestern Willow Flycatchers than other winter regions. We think our inconclusive results are due primarily to weak and nonlinear gradients in isotope ratios across the winter range of Willow Flycatchers.     

Range-wide effects of breeding- and nonbreeding-season climate on the abundance of a Neotropical migrant songbird

Geographic variation in the population dynamics of a species can result from regional variability in climate and how it affects reproduction and survival. Identifying such effects for migratory birds requires the integration of population models with knowledge of migratory connectivity between breeding and nonbreeding areas. We used Bayesian hierarchical models with 26 years of Breeding Bird Survey data (1982-2007) to investigate the impacts of breeding- and nonbreeding-season climate on abundance of American Redstarts (Setophaga ruticilla) across the species range. We focused on 15 populations defined by Bird Conservation Regions, and we included variation across routes and observers as well as temporal trends and climate effects. American Redstart populations that breed in eastern North America showed increased abundance following winters with higher plant productivity in the Caribbean where they are expected to overwinter. In contrast, western breeding populations showed little response to conditions in their expected wintering areas in west Mexico, perhaps reflecting lower migratory connectivity or differential effects of winter rainfall on individuals across the species range. Unlike the case with winter climate, we found few effects of temperature prior to arrival in spring (March-April) or during the nesting period (May-June) on abundance the following year. Eight populations showed significant changes in abundance, with the steepest declines in the Atlantic Northern Forest (-3.4%/yr) and the greatest increases in the Prairie Hardwood Transition (4%/yr). This study emphasizes how the effects of climate on populations of migratory birds are context dependent and can vary depending on geographic location and the period of the annual cycle. Such knowledge is essential for predicting regional variation in how populations of a species might vary in their response to climate change.     

Connectivity of wood thrush breeding,wintering, and migration sites based on range‐wide tracking

Many migratory animals are experiencing rapid population declines, but migration data with the geographic scope and resolution to quantify the complex network of movements between breeding and nonbreeding regions are often lacking. Determining the most frequently used migration routes and nonbreeding regions for a species is critical for understanding population dynamics and making effective conservation decisions. We tracked the migration of individual Wood Thrushes (Hylocichla mustelina) (n = 102) from across their range with light‐level geolocators and, for the first time, quantified migration routes and wintering regions for distinct breeding populations. We identified regional and species‐level migratory connectivity networks for this declining songbird by combining our tracking results with range‐wide breeding abundance estimates and forest cover data. More than 50% of the species occupied the eastern wintering range (Honduras to Costa Rica), a region that includes only one‐third of all wintering habitat and that is undergoing intensive deforestation. We estimated that half of all Wood Thrushes in North America migrate south through Florida in fall, whereas in spring approximately 73% funnel northward through a narrow span along the central U.S. Gulf Coast (88–93°W). Identifying migratory networks is a critical step for conservation of songbirds and we demonstrated with Wood Thrushes how it can highlight conservation hotspots for regional populations and species as a whole. Conectividad de Sitios de Reproducción, Invierno y Migración del Zorzal con Base en Rastreo de Cobertura Amplia     

Escape migration decisions in Eurasian Woodcocks: insights from survival analyses using large-scale recovery data

During unpredictable adverse conditions, endotherms can engage in emergency behaviors (movement, torpor, hyperphagia) to maintain energy balance and reduce mortality hazards. Bird “escape migration” is one of the most visible of these behaviors. In this study, we focus on a Eurasian Woodcock Scolopax rusticola population. Seasonal migrations bring this population from its breeding grounds in Eastern and Northern Europe to its wintering grounds in France. A varying number of these birds are also regularly reported from Spain, supposedly during additional escape movements that occur in winter. Using models that account for the imperfect detection rate of individuals and a large (>44,000 individuals) dataset combining information from the wintering and breeding ranges, we show that severe winters significantly reduced survival probability, but that migration to Spain increased only during the most intense cold spell that occurred over the 20-year study period. This suggests that the decision to resume migration during the winter is submitted to a threshold mechanism, which we discuss in the light of current models of migratory behavior.     

Are Declines in North American Insectivorous Songbirds Due to Causes on the Breeding Range?

Previous studies have reported a recent decline in breeding populations of migratory songbirds in eastern and central North America. Several explanations have been suggested: deforestation on the wintering grounds in the tropics and habitat loss, increased predation pressure, and increased cowbirds parasitism on the breeding range. We used these factors to assign 47 species of insectivorous passerines to groups with contrasting vulnerability, and then used the North America Breeding Bird Survey to analyze population trends in these groups on a large continental scale. Variables indexing susceptibility to predation on the breeding ground were most strongly correlated with population trends form 1968 to 1987. During the period from 1978 to 1987, migratory status was also significantly associated with population trends long-distance migrants to the neotropics exhibited a small, nonsignificant decreasing trend, whereas residents and short-distance migrants increased strongly. During the same time period, the group of species with low nest location, open nest, and high cowbird parasitism declined significantly. Although it is difficult to separate the effects of multiple factors, our analyses suggest that predation on the breeding ground in North America has played a larger role in the decline of migratory songbirds than deforestation on the wintering grounds in the tropics.     

Effective Population Size in Red-Cockaded Woodpeckers: Population and Model Differences

Loss of genetic variability in isolated populations is an important issue for conservation biology. Most studies involve only a single population of a given species and a single method of estimating rate of loss. Here we present analyses for three different Red-cockaded Woodpecker ( Picoides borealis ) populations from different geographic regions. We compare two different models for estimating the expected rate of loss of genetic variability, and test their sensitivity to model parameters. We found that the simpler model (Reed et al. 1988) consistently estimated a greater rate of loss of genetic variability from a population than did the Emigh and Pollak (1979) model. The ratio of effective population size (which describes the expected rate of loss of genetic variability) to breeder population size varied widely among Red-cockaded Woodpecker populations due to geographic variation in demography. For this species, estimates of effective size were extremely sensitive to survival parameters, but not to the probability of breeding or reproductive success. Sensitivity was sufficient that error in estimating survival rates in the field could easily mask true population differences in effective size. Our results indicate that accurate and precise demographic data are prerequisites to determining effective population size for this species using genetic models, and that a single estimate of rate of loss of genetic variability is not valid across populations.     

Geographic assignment of seabirds to their origin: combining morphologic, genetic, and biogeochemical analyses.

Longline fisheries, oil spills, and offshore wind farms are some of the major threats increasing seabird mortality at sea, but the impact of these threats on specific populations has been difficult to determine so far. We tested the use of molecular markers, morphometric measures, and stable isotope (delta15N and delta13C) and trace element concentrations in the first primary feather (grown at the end of the breeding period) to assign the geographic origin of Calonectris shearwaters. Overall, we sampled birds from three taxa: 13 Mediterranean Cory's Shearwater (Calonectris diomedea diomedea) breeding sites, 10 Atlantic Cory's Shearwater (Calonectris diomedea borealis) breeding sites, and one Cape Verde Shearwater (C. edwardsii) breeding site. Assignment rates were investigated at three spatial scales: breeding colony, breeding archipelago, and taxa levels. Genetic analyses based on the mitochondrial control region (198 birds from 21 breeding colonies) correctly assigned 100% of birds to the three main taxa but failed in detecting geographic structuring at lower scales. Discriminant analyses based on trace elements composition achieved the best rate of correct assignment to colony (77.5%). Body measurements or stable isotopes mainly succeeded in assigning individuals among taxa (87.9% and 89.9%, respectively) but failed at the colony level (27.1% and 38.0%, respectively). Combining all three approaches (morphometrics, isotopes, and trace elements on 186 birds from 15 breeding colonies) substantially improved correct classifications (86.0%, 90.7%, and 100% among colonies, archipelagos, and taxa, respectively). Validations using two independent data sets and jackknife cross-validation confirmed the robustness of the combined approach in the colony assignment (62.5%, 58.8%, and 69.8% for each validation test, respectively). A preliminary application of the discriminant model based on stable isotope delta15N and delta13C values and trace elements (219 birds from 17 breeding sites) showed that 41 Cory's Shearwaters caught by western Mediterranean long-liners came mainly from breeding colonies in Menorca (48.8%), Ibiza (14.6%), and Crete (31.7%). Our findings show that combining analyses of trace elements and stable isotopes on feathers can achieve high rates of correct geographic assignment of birds in the marine environment, opening new prospects for the study of seabird mortality at sea.     

A genoscape-network model for conservation prioritization in a migratory bird

Migratory animals are declining worldwide and coordinated conservation efforts are needed to reverse current trends. We devised a novel genoscape-network model that combines genetic analyses with species distribution modeling and demographic data to overcome challenges with conceptualizing alternative risk factors in migratory species across their full annual cycle. We applied our method to the long distance, Neotropical migratory bird, Wilson's Warbler (Cardellina pusilla). Despite a lack of data from some wintering locations, we demonstrated how the results can be used to help prioritize conservation of breeding and wintering areas. For example, we showed that when genetic, demographic, and network modeling results were considered together it became clear that conservation recommendations will differ depending on whether the goal is to preserve unique genetic lineages or the largest number of birds per unit area. More specifically, if preservation of genetic lineages is the goal, then limited resources should be focused on preserving habitat in the California Sierra, Basin Rockies, or Coastal California, where the 3 most vulnerable genetic lineages breed, or in western Mexico, where 2 of the 3 most vulnerable lineages overwinter. Alternatively, if preservation of the largest number of individuals per unit area is the goal, then limited conservation dollars should be placed in the Pacific Northwest or Central America, where densities are estimated to be the highest. Overall, our results demonstrated the utility of adopting a genetically based network model for integrating multiple types of data across vast geographic scales and better inform conservation decision-making for migratory animals.     

Tracing the geographic origin of traded leopard body parts in the indian subcontinent with DNA‐based assignment tests

Illicit trade in wildlife products is rapidly decimating many species across the globe. Such trade is often underestimated for wide‐ranging species until it is too late for the survival of their remaining populations. Policing this trade could be vastly improved if one could reliably determine geographic origins of illegal wildlife products and identify areas where greater enforcement is needed. Using DNA‐based assignment tests (i.e., samples are assigned to geographic locations), we addressed these factors for leopards (Panthera pardus) on the Indian subcontinent. We created geography‐specific allele frequencies from a genetic reference database of 173 leopards across India to infer geographic origins of DNA samples from 40 seized leopard skins. Sensitivity analyses of samples of known geographic origins and assignments of seized skins demonstrated robust assignments for Indian leopards. We found that confiscated pelts seized in small numbers were not necessarily from local leopards. The geographic footprint of large seizures appeared to be bigger than the cumulative footprint of several smaller seizures, indicating widespread leopard poaching across the subcontinent. Our seized samples had male‐biased sex ratios, especially the large seizures. From multiple seized sample assignments, we identified central India as a poaching hotspot for leopards. The techniques we applied can be used to identify origins of seized illegal wildlife products and trade routes at the subcontinent scale and beyond.     

Sampling variability and estimates of density dependence: a composite-likelihood approach

It is well known that sampling variability, if not properly taken into account, affects various ecologically important analyses. Statistical inference for stochastic population dynamics models is difficult when, in addition to the process error, there is also sampling error. The standard maximum-likelihood approach suffers from large computational burden. In this paper, I discuss an application of the composite-likelihood method for estimation of the parameters of the Gompertz model in the presence of sampling variability. The main advantage of the method of composite likelihood is that it reduces the computational burden substantially with little loss of statistical efficiency. Missing observations are a common problem with many ecological time series. The method of composite likelihood can accommodate missing observations in a straightforward fashion. Environmental conditions also affect the parameters of stochastic population dynamics models. This method is shown to handle such nonstationary population dynamics processes as well. Many ecological time series are short, and statistical inferences based on such short time series tend to be less precise. However, spatial replications of short time series provide an opportunity to increase the effective sample size. Application of likelihood-based methods for spatial time-series data for population dynamics models is computationally prohibitive. The method of composite likelihood is shown to have significantly less computational burden, making it possible to analyze large spatial time-series data. After discussing the methodology in general terms, I illustrate its use by analyzing a time series of counts of American Redstart (Setophaga ruticilla) from the Breeding Bird Survey data, San Joaquin kit fox (Vulpes macrotis mutica) population abundance data, and spatial time series of Bull trout (Salvelinus confluentus) redds count data.     

Use of trace element analysis of feathers as a tool to track fine-scale dispersal in birds

Tracking dispersal and migratory movements of animals over small and large spatial scales is a challenge. In birds, a promising yet underutilized tool is the trace elemental composition of feathers. The elemental profile of a feather may reveal information about the geographic origin of a bird provided that molting occurs on the breeding grounds and that elemental differences exist between breeding areas. Here, we explore the use of trace elemental composition in body feathers of the Puget Sound white-crowned sparrow Zonotrichia leucophrys pugetensis as a tool (1) to discriminate among birds collected in four different song dialect populations along a 400-km stretch of the Pacific Northwest coast and (2) to assign males singing nonlocal dialects in one population to potential natal populations. Inductively coupled plasma mass spectrometry detected 34 trace elements in sampled feathers and in a discriminant function analysis seven of these elements differed among the four source populations. Half of the six nonlocal dialect singers, who were likely to have immigrated into the focal population, were assigned to a population that matched their song dialect. Our study suggests that feather microchemistry is a promising tool for identifying geographic origins of dispersing birds over small geographic scales and in combination with other markers, such as song, may give insight into ecological and evolutionary processes.     

Efficacy of Land-Cover Models in Predicting Isolation of Marbled Salamander Populations in a Fragmented Landscape

Abstract:  Amphibians worldwide are facing rapid declines due to habitat loss and fragmentation, disease, and other causes. Where habitat alteration is implicated, there is a need for spatially explicit conservation plans. Models built with geographic information systems (GIS) are frequently used to inform such planning. We explored the potential for using GIS models of functional landscape connectivity as a reliable proxy for genetically derived measures of population isolation. We used genetic assignment tests to characterize isolation of marbled salamander populations and evaluated whether the relative amount of modified habitat around breeding ponds was a reliable indicator of population isolation. Using a resampling analysis, we determined whether certain land-cover variables consistently described population isolation. We randomly drew half the data for model building and tested the performance of the best models on the other half 100 times. Deciduous forest was consistently associated with lower levels of population isolation, whereas salamander populations in regions of agriculture and anthropogenic development were more isolated. Models that included these variables and pond size explained 65–70% of variation in genetically inferred isolation across sites. The resampling analysis confirmed that these habitat variables were consistently good predictors of isolation. Used judiciously, simple GIS models with key land-cover variables can be used to estimate population isolation if field sampling and genetic analysis are not possible.     

The Impact of Variation in Stopover Habitat Quality on Migrant Rufous Hummingbirds

We examined the impact of variation in habitat quality on migrating Rufous Hummingbirds ( Selasphorus rufus ) in the California Sierra Nevada. As do other migratory species, these birds depend on "stopover" habitats en route for feeding and replenishing depleted energy stores. During seven years of study, the quality of the stopover habitat (assessed in terms of the density of nectar food resources) varied widely due to natural variation in flowering. In years when stopover habitat quality was poor, incoming body masses were low and stopover durations were long. Population densities of migrant hummingbirds at the study site were coupled to habitat quality both within and among years. These observations demonstrate important effects of stopover habitat variation on the physiological, behavioral, and population ecology of migrating hummingbirds. High-quality stopover habitats are critical links between breeding and wintering areas for many species, and their preservation should be considered an essential component of strategies aiming to conserve migratory bird populations.     

Sulfur isotopes link overwinter habitat use and breeding condition in Double-crested Cormorants.

North American Double-crested Cormorant (Phalacrocorax auritus) populations have increased greatly. Both breeding and overwintering ground factors have likely contributed to these increases. However, demonstrating how overwintering conditions may affect breeding birds has not been possible because of the difficulty in linking breeding birds to their wintering grounds. Here, we demonstrate the utility of stable sulfur isotopes to elucidate overwintering habitat use by cormorants breeding on Lake Erie. Sulfur isotopes in feathers grown on overwintering grounds provided insights into the degree to which birds used freshwater vs. marine environments. The proportion of birds utilizing freshwater habitats increased through time. This change may have reflected increases in freshwater aquaculture (i.e., catfish) in the U.S. south. Examination of body condition in birds returning to breed on Lake Erie indicated that those individuals that solely used marine habitats for at least a portion of the overwintering period were in poorer condition than birds using freshwater. Enhanced foraging opportunities at aquaculture facilities may improve the fitness of individuals that have returned to breed after overwintering at such locations. This study is the first to demonstrate a linkage between overwinter habitat use and breeding ground parameters in Double-crested Cormorants. These results underscore that factors throughout the Mississippi flyway are likely acting together to regulate cormorant populations.     

Estimating species-specific survival and movement when species identification is uncertain

Incorporating uncertainty in the investigation of ecological studies has been the topic of an increasing body of research. In particular, mark-recapture methodology has shown that incorporating uncertainty in the probability of detecting individuals in populations enables accurate estimation of population-level processes such as survival, reproduction, and dispersal. Recent advances in mark-recapture methodology have included estimating population-level processes for biologically important groups despite the misassignment of individuals to those groups. Examples include estimating rates of apparent survival despite less than perfect accuracy when identifying individuals to gender or breeding state. Here we introduce a method for estimating apparent survival and dispersal in species that co-occur but that are difficult to distinguish. We use data from co-occurring populations of meadow voles (Microtus pennsylvanicus) and montane voles (M. montanus) in addition to simulated data to show that ignoring species uncertainty can lead to biased estimates of population processes. The incorporation of species uncertainty in mark-recapture studies should aid future research investigating ecological concepts such as interspecific competition, niche differentiation, and spatial population dynamics in sibling species.     

Combining multistate capture-recapture data with tag recoveries to estimate demographic parameters

Matrix population models that allow an animal to occupy more than one state over time are important tools for population and evolutionary ecologists. Definition of state can vary, including location for metapopulation models and breeding state for life history models. For populations whose members can be marked and subsequently reencountered, multistate mark-recapture models are available to estimate the survival and transition probabilities needed to construct population models. Multistate models have proved extremely useful in this context, but they often require a substantial amount of data and restrict estimation of transition probabilities to those areas or states subjected to formal sampling effort. At the same time, for many species, there are considerable tag recovery data provided by the public that could be modeled in order to increase precision and to extend inference to a greater number of areas or states. Here we present a statistical model for combining multistate capture-recapture data (e.g., from a breeding ground study) with multistate tag recovery data (e.g., from wintering grounds). We use this method to analyze data from a study of Canada Geese (Branta canadensis) in the Atlantic Flyway of North America. Our analysis produced marginal improvement in precision, due to relatively few recoveries, but we demonstrate how precision could be further improved with increases in the probability that a retrieved tag is reported.     

Measuring terrestrial subsidies to aquatic food webs using stable isotopes of hydrogen

Understanding river food webs requires distinguishing energy derived from primary production in the river itself (autochthonous) from that produced externally (allochthonous), yet there are no universally applicable and reliable techniques for doing so. We compared the natural abundance stable isotope ratios of hydrogen (deltaD) of allochthonous and autochthonous energy sources in four different aquatic ecosystems. We found that autochthonous organic matter is uniformly far more depleted in deuterium (lower deltaD values) than allochthonous: an average difference of approximately 100% per hundred. We also found that organisms at higher trophic levels, including both aquatic invertebrates and fish, have deltaD values intermediate between aquatic algae and terrestrial plants. The consistent differences between leaves and algae in deltaD among these four watersheds, along with the intermediate values in higher trophic levels, indicate that natural abundance hydrogen isotope signatures are a powerful tool for partitioning energy flow in aquatic ecosystems.     

Estimating age from recapture data: integrating incremental growth measures with ancillary data to infer age-at-length

Estimating the age of individuals in wild populations can be of fundamental importance for answering ecological questions, modeling population demographics, and managing exploited or threatened species. Significant effort has been devoted to determining age through the use of growth annuli, secondary physical characteristics related to age, and growth models. Many species, however, either do not exhibit physical characteristics useful for independent age validation or are too rare to justify sacrificing a large number of individuals to establish the relationship between size and age. Length-at-age models are well represented in the fisheries and other wildlife management literature. Many of these models overlook variation in growth rates of individuals and consider growth parameters as population parameters. More recent models have taken advantage of hierarchical structuring of parameters and Bayesian inference methods to allow for variation among individuals as functions of environmental covariates or individual-specific random effects. Here, we describe hierarchical models in which growth curves vary as individual-specific stochastic processes, and we show how these models can be fit using capture-recapture data for animals of unknown age along with data for animals of known age. We combine these independent data sources in a Bayesian analysis, distinguishing natural variation (among and within individuals) from measurement error. We illustrate using data for African dwarf crocodiles, comparing von Bertalanffy and logistic growth models. The analysis provides the means of predicting crocodile age, given a single measurement of head length. The von Bertalanffy was much better supported than the logistic growth model and predicted that dwarf crocodiles grow from 19.4 cm total length at birth to 32.9 cm in the first year and 45.3 cm by the end of their second year. Based on the minimum size of females observed with hatchlings, reproductive maturity was estimated to be at nine years. These size benchmarks are believed to represent thresholds for important demographic parameters; improved estimates of age, therefore, will increase the precision of population projection models. The modeling approach that we present can be applied to other species and offers significant advantages when multiple sources of data are available and traditional aging techniques are not practical.     

Climate change threatens polar bear populations: a stochastic demographic analysis

The polar bear (Ursus maritimus) depends on sea ice for feeding, breeding, and movement. Significant reductions in Arctic sea ice are forecast to continue because of climate warming. We evaluated the impacts of climate change on polar bears in the southern Beaufort Sea by means of a demographic analysis, combining deterministic, stochastic, environment-dependent matrix population models with forecasts of future sea ice conditions from IPCC general circulation models (GCMs). The matrix population models classified individuals by age and breeding status; mothers and dependent cubs were treated as units. Parameter estimates were obtained from a capture-recapture study conducted from 2001 to 2006. Candidate statistical models allowed vital rates to vary with time and as functions of a sea ice covariate. Model averaging was used to produce the vital rate estimates, and a parametric bootstrap procedure was used to quantify model selection and parameter estimation uncertainty. Deterministic models projected population growth in years with more extensive ice coverage (2001-2003) and population decline in years with less ice coverage (2004-2005). LTRE (life table response experiment) analysis showed that the reduction in lambda in years with low sea ice was due primarily to reduced adult female survival, and secondarily to reduced breeding. A stochastic model with two environmental states, good and poor sea ice conditions, projected a declining stochastic growth rate, log lambdas, as the frequency of poor ice years increased. The observed frequency of poor ice years since 1979 would imply log lambdas approximately - 0.01, which agrees with available (albeit crude) observations of population size. The stochastic model was linked to a set of 10 GCMs compiled by the IPCC; the models were chosen for their ability to reproduce historical observations of sea ice and were forced with "business as usual" (A1B) greenhouse gas emissions. The resulting stochastic population projections showed drastic declines in the polar bear population by the end of the 21st century. These projections were instrumental in the decision to list the polar bear as a threatened species under the U.S. Endangered Species Act.