Climate change is having widespread ecological effects, including loss of Arctic sea ice. This has led to listing of the polar bear (Ursus maritimus) and other ice-dependent marine mammals under the U.S. Endangered Species Act (ESA). Methods are needed to evaluate the effects of climate change on population persistence to inform recovery planning for listed species. For polar bears, this includes understanding interactions between climate and secondary factors, such as subsistence harvest, which provide economic, nutritional, or cultural value to humans.

We developed a matrix-based demographic model for polar bears that can be used for population viability analysis and to evaluate the effects of human-caused removals. This model includes density-dependence (the potential for a declining environmental carrying capacity), density-independent limitation, and sex- and age-specific harvest vulnerabilities. We estimated values of adult female survival (0.93–0.96), recruitment (number of yearling cubs per adult female; 0.1–0.3), and carrying capacity (>250 animals) that must be maintained for a hypothetical population to achieve a 90-percent probability of persistence over 100 years.

We also developed a state-dependent management framework, based on harvest theory and the potential biological removal method, by linking the demographic model to simulated population assessments. This framework can be used to estimate the maximum sustainable rate of human-caused removals, including subsistence harvest, which maintains a population at its maximum net productivity level. The framework also can be used to calculate a recommended sustainable harvest rate, which generally is lower than the maximum sustainable rate and depends on management objectives, the precision and frequency of population data, and risk tolerance. The historical standard 4.5-percent harvest rate for polar bears, at a 2:1 male-to-female ratio, is reasonable under many biological and management conditions, although lower or higher rates may be appropriate in some cases.

Our modeling results suggest that harvest of polar bears is unlikely to accelerate population declines that result from declining carrying capacity caused by sea-ice loss, provided that several conditions are met: (1) the sustainable harvest rate reflects the population’s intrinsic growth rate, and the corresponding harvest level is obtained by applying this rate to an estimate of population size; (2) the sustainable harvest rate reflects the quality of population data (e.g., lower harvest when data are poor); and (3) the level of human-caused removals can be adjusted. Finally, our results suggest that stopgap measures (e.g., further reduction or cessation of harvest when the population size is less than a critical threshold) may be necessary to minimize the incremental risk associated with harvest, if environmental conditions are deteriorating rapidly. We suggest that the demographic model and approaches presented here can serve as a template for conservation planning for polar bears and other species facing similar challenges.