Abstract Recent observations suggest that polar bears (Ursus maritimus) are increasingly using land habitats in some parts of their range, where they have minimal access to their preferred prey, likely in response to loss of their sea ice habitat associated with climatic warming. We used location data from female polar bears fit with satellite radio collars to compare land use patterns in the Chukchi Sea between two periods (1986–1995 and 2008–2013) when substantial summer sea-ice loss occurred. In both time periods, polar bears predominantly occupied sea-ice, although land was used during the summer sea-ice retreat and during the winter for maternal denning. However, the proportion of bears on land for > 7 days between August and October increased between the two periods from 20.0% to 38.9%, and the average duration on land increased by 30 days. The majority of bears that used land in the summer and for denning came to Wrangel and Herald Islands (Russia), highlighting the importance of these northernmost land habitats to Chukchi Sea polar bears. Where bears summered and denned, and how long they spent there, was related to the timing and duration of sea ice retreat. Our results are consistent with other studies supporting increased land use as a common response of polar bears to sea-ice loss. Implications of increased land use for Chukchi Sea polar bears are unclear, because a recent study observed no change in body condition or reproductive indices between the two periods considered here. This result suggests that the ecology of this region may provide a degree of resilience to sea ice loss. However, projections of continued sea ice loss suggest that polar bears in the Chukchi Sea and other parts of the Arctic may increasingly use land habitats in the future, which has the potential to increase nutritional stress and human-polar bear interactions.

Citation: Rode KD, Wilson RR, Regehr EV, St. Martin M, Douglas DC, Olson J (2015) Increased Land Use by Chukchi Sea Polar Bears in Relation to Changing Sea Ice Conditions. PLoS ONE 10(11): e0142213. https://doi.org/10.1371/journal.pone.0142213 Editor: Daniel E. Crocker, Sonoma State University, UNITED STATES Received: August 19, 2015; Accepted: October 19, 2015; Published: November 18, 2015 This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication Data Availability: Location data collected from Chukchi Sea female polar bears 1986-1995 are available at doi:10.5066/F7X34VHM and 2008-2013 at doi:10.5061/dryad.n5509. Funding: Funding for this study was provided by the USFWS Marine Mammals Management, USGS Changing Arctic Ecosystem Initiative, the Detroit Zoological Association and the State of Alaska Coastal Impact Assessment Program. Only the USFWS and USGS contributed to study design, data collection and analysis, decision to publish and preparation of the manuscript. In-kind support was received from Teck Alaska, Inc. Competing interests: Teck Alaska Inc provided in-kind support for the study including use of their facilities as a base for field operations. This in-kind support does not alter the authors' adherence to PLOS ONE policies on sharing data and materials.

Introduction The Arctic is experiencing rapid and extensive sea ice loss that is projected to continue as a result of increasing levels of greenhouse gas emissions [1]. Among Arctic marine mammals, the response of polar bears (Ursus maritimus) to ecological change has received the most attention relative to global warming impacts [2]. Much of the information about polar bears and sea ice loss, however, is based on only a few of the 19 recognized subpopulations [3], particularly the Western Hudson Bay, Southern Beaufort Sea, and Barents Sea subpopulations. Western Hudson Bay is seasonally ice free and the entire subpopulation comes onshore during the summer where they survive primarily by metabolizing accumulated stores of fat and lean mass [4]. In the southern Beaufort and Barents seas, the ice retreats northward in the summer. Some bears in both subpopulations come onshore but the majority appear to move north with the ice [5, 6]. With recent sea ice loss, bears in these three subpopulations have increasingly spent time on shore in the summer [5, 7, 8, 9]. However, both annual sea ice dynamics and polar bear life history strategies vary across the Arctic, and in many regions, baseline data on land use and its relationship with changing sea ice conditions are unknown. The degree to which increased land use is a common response of polar bears to sea ice loss is relevant to the state of Arctic ecosystems, humans living and working in the Arctic, and polar bear conservation and management. Several studies in areas where bears are coming onshore earlier and spending more time on shore have documented localized but intense nest predation of sea birds and geese [9, 10, 11, 12, 13]. On Wrangel Island, Russia, polar bear predation resulted in the death of at least 226 walruses (Odobenus rosmarus) between 1989 and 2006, of which 81% died as a result of stampedes caused by polar bears [14]. Grizzly bears (Ursus arctos) occupy terrestrial habitats in many parts of the polar bear’s range, such that increased land use could increase sympatry and potential competition [15]. Hybridization between grizzly bears and polar bears has been documented in northwestern Canada [16]. Further, increases in terrestrial habitat use can lead to increased interactions with humans that may result in safety concerns for humans or in disturbance or death for polar bears [17]. In addition to the effects listed above, increased land use is likely to have nutritional consequences for polar bears. While on land, polar bears have little or no access to ice seals, their primary prey. Although numerous studies have documented that polar bears consume terrestrially-derived foods [18], scavenge beached marine mammal carcasses, eat remains of Native subsistence harvests [5], and actively hunt hauled-out walruses [14], most terrestrial feeding is opportunistic and of lower energetic content compared to the marine prey that polar bears consume on the sea ice [19, 20]. Thus, terrestrial foods are likely capable of supporting only a small proportion of current polar bear numbers. In seasonally ice-free habitats of eastern Canada, such as western Hudson Bay, analyses of polar bear blood chemistry [21], feeding habits via stable isotopes [22, 23] and body mass and fat changes indicate that polar bears are primarily fasting while on land [4]. Longer periods fasting on shore has been suggested as the primary mechanism by which body condition, reproduction, and ultimately polar bear survival in Western Hudson Bay have declined in conjunction with sea ice loss [7, 24, 25]. A large portion of the polar bear’s circumpolar range occurs in areas where sea ice retreats towards the Arctic basin and away from land masses during the summer [26], requiring bears to make decisions as to whether to stay on land or remain with the retreating ice both during the summer and for denning. Little is known about the feeding patterns of polar bears that remain on the sea ice over the Arctic basin throughout the summer. Prey densities in deep waters beyond the continental shelf are thought to be low [27], and polar bears that summer on the sea ice in the southern Beaufort Sea appear to have energy expenditures typical of fasting, non-hibernating mammals [28]. Although a majority of bears in the southern Beaufort Sea remain with the sea ice during the summer, pregnant females primarily den on land [29]. In areas such as this where land-based denning has increased [29] or changed in distribution [30], the ability of bears to maintain access to terrestrial denning areas without compromising foraging opportunities before or after denning may also influence cub survival [31]. There is some evidence that important spring foraging habitats near terrestrial denning areas in Svalbard have declined [32]. Thus, the decisions a bear makes as to whether to use land prior to and during maternal denning may have nutritional implications for pregnant females [33]. Because increased land use by polar bears has implications for Arctic ecosystem ecology, bear-human interactions, and polar bear health, here, we examined land use by adult female polar bears in the Chukchi Sea; members of a subpopulation which range between Alaska and Russia. We compared land use during two time periods (1986–1995 and 2008–2013), when substantial reductions in sea ice habitat occurred [34]. A portion of bears in this subpopulation have typically summered [14, 35] and denned on Wrangel and Herald Islands [36] and the northeast coast of Chukotka [37]. Sea-ice dynamics are similar to the adjacent southern Beaufort Sea, as well as regions throughout Russia and Norway, where the sea ice retreats towards the pole during the annual sea ice minimum [26]. Thus, we sought to determine if the Chukchi Sea subpopulation has similarly increased duration on land during the summer and increased use of terrestrial habitats for denning, which would further support that increased land use may be a common response of polar bears to sea ice loss. Our objectives were to 1) determine which land masses were used by female polar bears; 2) compare the proportion of female polar bears in the Chukchi Sea subpopulation using land, and the amount of time spent on land, between the two time periods; and 3) evaluate environmental factors that influence the location, timing, and duration of land use.

Materials and Methods Ethics Statement This research was approved under the Marine Mammal Protection Act and Endangered Species Act with U.S. Fish and Wildlife Service (USFWS) permit number MA046081. Capture protocols were approved by the U.S. Geological Survey, (USGS) Alaska Science Center and USFWS Region 7 Institutional Animal Care and Use Committees. Polar bear capture and collaring Polar bears were captured and released on sea ice or land between mid-March and early May in the Chukchi and northern Bering Seas 1986–1995 and 2008–2013 (Fig 1). The study area ranged from 58 to 83°N and -160 to 166°W. In 1986–1995, 74 capture events occurred on or near Wrangel Island or Herald Island, 34 on sea ice off the Alaskan mainland coast, 13 on sea ice off of St. Lawrence Island, Alaska, and 4 on sea ice off the northeast Chukotkan coast. In 2008–2013, all polar bear captures occurred on sea ice off the Alaskan mainland coast (Fig 1). Polar bears were located from a helicopter and immobilized with a dart containing either M-99 in 1986 or zolazepam-tiletamine (Telazol® or Zoletil®) [38] in 1986–2013. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Fig 1. Locations where adult female polar bears were captured and fitted with satellite radio collars in 1986–1995 and 2008–2013. The area used for ice analysis is outlined in black. https://doi.org/10.1371/journal.pone.0142213.g001 Only adult female bears (≥ 4 years old) can safely be fitted with collars because they are no longer growing. Adult males were not collared because their necks are larger than their heads, preventing collar retention. Between 1986 and 1995, 103 adult females were fitted with satellite radio collars (Telonics Inc., Mesa, AZ), some of which were captured multiple times and/or wore collars for multiple years resulting in 127 bear-years of data. Collars deployed between 1986 and 1995 transmitted to the Argos System (Argos, www.argos-system.org) every 3 days. Between 2008 and 2013, 47 adult females were fitted with satellite radio collars providing 62 bear-years of data. Argos-only collar deployments in 2008–2010 provided Argos locations every one to three days; Global Positioning System (GPS) collar deployments 2010–2013 provided GPS locations four times per day. A bear-year was defined from the time of capture in the spring through at least the following January, so that summer habitat use and denning behavior could be identified. Estimating land use from location data Prior to 2010, all radio collar locations were determined by the Argos System with accuracy from < 250 m to > 1500 m (see http://www.argos-system.org/web/en/78-faq.php#faq-theme-55). From 2010 to 2013, GPS locations were transmitted via Argos. Location data were filtered to remove implausible locations using the Douglas Argos-Filter algorithm [39], which retained all standard quality class locations (classes 3, 2, and 1), rejected all class Z locations, and retained auxiliary class locations (0, A, and B) if they were corroborated by a consecutive location within 10 km, or if movement rates were < 10 km hr−1 and turning angles were not extremely acute. Instances of dropped collars were identified based on activity and temperature sensors, and data collected post-drop were removed. To integrate tracking data obtained from GPS and Argos, we employed the continuous time correlated random walk (CRAWL) model [40] to predict locations on a regularized, daily time interval (daily). The CRAWL model accounts for variable location quality and sampling intervals. For Argos locations, we used accuracy estimates provided by Telonics, Inc. (i.e., 3: 150 m, 2: 350 m, 1: 1000 m, 0: 1500 m; http://www.telonics.com/technotes/argosintro.php). Because location accuracies were not provided by Telonics for auxiliary location classes A or B, we prescribed conservative location accuracies of 5,000 and 10,000 m, respectively. We assigned locations obtained from GPS collars an accuracy of 30 m [41]. Using the observed locations and associated accuracies, we used the R statistical computing [42] package ‘CRAWL’ [43] to predict daily polar bear locations from 1 Jul to at least 31 Jan the following year, a period that encompasses the annual sea ice minimum and the initiation of denning. Because the uncertainty of CRAWL location estimates increases as the time between observed locations increases [44], we excluded predicted locations that occurred between observed locations separated by >14 days. The majority of data gaps filled with the CRAWL model were less than 5 days apart (1986–1995: 7.2% of locations were >4 days apart; 2008–2013: 1.3% of locations were > 4 days part). To validate that the CRAWL model was able to predict land use when there was a 14 day gap between observations, we found all observed land use locations (i.e., midpoints in a 14 day series of locations) where there was another observed location 7 days before and 7 days after that location. We then removed all points in the 14 day period, and used CRAWL to predict whether the midpoint occurred on land. We found that the majority of predicted locations occurred on land (mean = 0.74, SE = 0.03), which validated our inclusion of location data with up to 14-day gaps. We also performed cross-validation by resampling observed data, to evaluate what proportion of predicted locations would be misclassified as occurring on land or ice. For each individual’s path, we withheld five observed locations and their associated land status (on land or on ice). We then predicted the location of the animal at the same time as the withheld location and determined whether the predicted point occurred on land or ice. For data collected 1986–1995, 4.4% of withheld locations were misclassified, whereas the percent of misclassified locations for data collected 2008–2013 was 2.1%. Thus, the CRAWL model performed well at correctly predicting the land-use status of withheld locations. We considered a bear to be on land if its predicted location was within 5 km of land as identified by the Global Self-consistent, Hierarchical, High-resolution, Geographic Database (GSHHG version 2.3.4; http://www.soest.hawaii.edu/pwessel/gshhg/). The 5 km buffer was used to encompass small barrier islands that may receive heavy use by polar bears in the summer [5] but were not depicted as land in the GSHHG, and to account for low accuracy of some locations. Our 5 km buffer might have resulted in some offshore bears being classified as on land, and vice versa, but this was less likely to occur during the focal time period of our analysis (August through October) because landfast ice is largely absent during this period [45] and the pack ice (the nearly continuous mass of ice that retreats northward during the summer) has generally receded >>5km north of the coast. Bears within 5 km of the coast during this time were likely to have been on land. Because ARGOS does not transmit when a bear is in the water, it is unlikely that locations received within 5 km on shore occurred in water. In the case of denning, we used a control chart-based algorithm [46] that classified denning events based exclusively on collar temperature. We defined denning as a consecutive bout of higher than expected temperatures (determined using temperatures from non-denning bears) lasting 35 days or more during winter months (October-May). This control chart-based algorithm for identifying denning behavior agreed with 94.5% of denning events confirmed via observations during VHF radio tracking (n = 73)[46] and 96.3% of den identifications based on temperature, activity, and location data in Fischbach et al. [30]. Because the physical features of dens in snow banks often attenuate transmission, location data were typically unavailable or of low quality during the denning period. If at least one observed location at the start of, during, or at the end of the denning period occurred on land, it was assumed that the den occurred on land. Den events identified could have been successful maternity dens, unsuccessful maternity dens that were abandoned early, or shelter dens used during periods of poor weather or reduced food availability [47]. We compared the proportion of bears that spent more than 7 days onshore and more than 21 days onshore (which included the bears that spent > 7 days onshore) between August and October using a binary logistic regression (i.e., 0 = did not use land > 7 or > 21 days, 1 = used land > 7 or > 21 days). We used the threshold of 21 days to indicate longer-term summering behavior that has been documented elsewhere [5] and 7 days as an indicator of shorter-term land use. Shorter-term land use was included because little is known about the duration of land use by Chukchi Sea polar bears. We also used binary logistic regression to compare substrate use (i.e., land versus ice habitats) during denning between the two time periods of our study. The number of days bears spent onshore (for bears spending > 7 or > 21 days onshore) between August and October was compared using an ANOVA because no individual bear was repeated in the data set (i.e., although some bears provided location data for more than one year, no individual provided data indicating land use > 7 or > 21 days for more than one year). Due to the potential for differences in capture location between time periods to affect the number of bears coming onshore, where they came to shore for denning and summering, and how long they spent there, we also compared these metrics using only bears that were captured in the same portion of the Alaska coast in 1986–1995 as where captured in 2008–2013. To understand the degree to which bears may use land during times of the year other than summer, we also used an ANOVA to compare the percent of locations among entire bear-years that were spent on land for bears with a minimum of 300 days of location data post-deployment. For linear models and logistic regressions, β-values are reported with standard errors (SE). For ANOVA comparisons, means are reported with standard deviations (stdev). Environmental factors: Ice metrics We were interested in the influence of summer sea ice conditions on where bears come to shore, how long they spend on shore, and whether they denned on land or sea ice. The sea ice minimum in the Arctic and in the Chukchi Sea occurs in September. Therefore, we focused on ice conditions during the months on either side of September (i.e., Aug-Oct) using daily sea ice concentration data from the National Snow and Ice Data Center (accessed 16 Oct 2014) 1986–2013 [48]. Several studies suggest that bears will remain on sea ice until the concentration drops below 15–30% [8, 49] and that they prefer ice habitats over shallow, continental shelf waters <300 m deep [8, 49, 50]. We hypothesized that the availability of ice between August and October would determine how long bears spent onshore. We determined the average daily proportion of the continental shelf in the ice analysis area (outlined in black in Fig 1) that was covered by ice of ≥15% concentration between August and October [34], and used this in a regression analysis with annual average days spent onshore as the response variable. We also hypothesized that sea ice conditions may affect where bears summer and den on land among the multiple land masses available within the Chukchi Sea, including the northeast Chukotkan coast, northwest Alaskan coast, and Wrangel and Herald Islands. To determine whether ice conditions were related to where bears came ashore in summer and for land-based denning (i.e., the response variable), we estimated the dates that sea ice (of ≥15% concentration) retreated from (hereafter “ice retreat date”) and returned to (hereafter “ice return date”) within <50 km from 25 x 25 km pixels of rasterized coastal segments along the Chukotkan coast, Alaskan coast, and Wrangel and Herald Islands. From the daily grids of sea ice concentration, we created grids of the average weekly sea ice concentration. We then used the R [42] package ‘SDMTools’ [51] to identify the largest, contiguous, area of weekly average ice with ≥ 15% concentration likely to be indicative of the mass of pack ice that retreats northward in the summer. If pack ice never retreated >50 km from a coastal segment, we assigned retreat and return date of 365 and 1, respectively. We chose 50 km from shore to determine ice retreat and return dates because it corresponds to a “long distance swim” for polar bears [52], which we hypothesized would correspond to an ecologically significant distance for individual animals. To determine if either ice retreat or return dates were related to where polar bears used land, either for denning or during the summer for >7 days, we used conditional logistic regression [53] to relate used areas to areas where use was not observed. For this analysis, we included all bears that were observed > 7 days onshore between August and October regardless of whether they had data gaps during that period. We used the ‘coxph’ function in the survival package [54] for R [42] to determine the potential relationship between ice formation and retreat dates and areas of summering and denning in a given year (i.e., in logistic regression 0 = unused and 1 = used). We built separate models for denning locations and summer land use locations. The time series of annual ice retreat dates was highly correlated with ice return dates (r = -0.80). Therefore, these variables were examined independently in models to determine which variable was most explanatory based on Akaike Information Criterion (AIC) scores (S1 Table)[55]. Each model included ice retreat or ice return date and an interaction term between time period (1986–1995 and 2008–2013) and the ice variable. The model with sea ice retreat was more explanatory than the model with ice return date (S1 Table), thus we used this model in our analysis. Because of the simplicity of the models, we did not implement a model selection process but rather assessed whether the model coefficients were significant at P ≤ 0.05. To ensure that there were no biases in denning and summer land use locations (northwest Alaskan coast, Wrangel/Herald Islands, northeast Chukotkan coast) based on where a bear was captured (i.e., northwestern Alaskan coast, St. Lawrence Island, Wrangel Island, northeast Chukotkan coast), we performed contingency table analyses [56] comparing capture location to ultimate denning or summering location. For examinations of the relationship between ice retreat and return dates and locations of dens on land, we included 6 land-based maternal den locations that were identified by direct observation rather than satellite tracking.

Discussion Polar bears in the Chukchi Sea relied heavily on sea ice habitats with greater than 45% of adult females in the subpopulation spending 95% or more of their time on the sea ice (n = 110). Land, however, was also necessary for critical life history functions, with >84% of denning of this subpopulation occurring on land (n = 70). Our observation that there was no change in the proportion of bears denning on land between the two time periods suggests that, despite differences in capture location (e.g., that a higher proportion of bears were captured on Wrangel Island in the earlier period), land-based denning was and continues to be important to the Chukchi Sea subpopulation. This result was consistent when comparing denning location of bears captured in similar areas for both time periods. Furthermore, land appears to be becoming an increasingly important summer habitat for bears as evidenced by a doubling in the proportion of bears coming to shore for an extended period, and an increase of 30 days in the amount of time spent onshore. The earlier arrival of some bears onshore is consistent with observational studies conducted on Wrangel Island [14, 57]. The relationship between land tenure and sea ice conditions suggests that increased land use is likely a result of a longer period of reduced ice conditions. Ice retreat dates in the northern part of the study area were found to occur 20 to 40 days earlier during 2008–2013 compared to 1986–1995. Other observational studies have similarly noted the relationship between the timing of bears arrival and departure on land and sea ice conditions [57, 58]. Pregnant females commonly came onshore during the summer prior to denning in both the earlier (37.5%; n = 16) and latter period (57.1%; n = 14). This differs from the southern Beaufort Sea where bears that came onshore to feed on whale remains from subsistence harvests during the summer did not subsequently den on land [30]. Thus, summering and denning on land may be associated in the Chukchi subpopulation but is less so in the adjacent southern Beaufort Sea subpopulation. These subpopulations also differ in that a much greater proportion of the Chukchi Sea subpopulation has historically used land for denning, whereas in the southern Beaufort Sea, land-based denning increased from 38 to 63% between 1985–1995 and 1998–2004 in apparent response to reduced sea ice availability or quality [30]. Instead of increasingly using land for denning with reduced sea ice, pregnant females in the Chukchi Sea appear to be selecting for different land masses. We found that denning on the Alaskan coast declined by 15% while denning on Wrangel and Herald Islands and the Chukotkan coast increased by 15%. This result is consistent with traditional and local ecological knowledge collected among Alaska Native hunters who reported that in the past 10 years dens have only been observed at the village of Point Lay and to the north, whereas historically some denning was observed south of Point Lay ([59]; this study). Our observation that ice was absent from southerly coastlines 126 days longer than more northerly coastlines suggests that denning in more southerly areas results in longer periods away from foraging habitats; and that perhaps a threshold in the southerly areas has been crossed that has decreased its use as a denning habitat. In land masses currently unused or minimally used for denning (e.g., the southern Alaska coast and far eastern Chukotkan coast) ice retreat begins approximately June 24 compared to more northern regions where ice retreat begins approximately Aug 6. Because sea ice loss in the Chukchi Sea is projected to continue [60], more northern land habitats, particularly Wrangel and Herald islands, may increasingly be used for summering and denning. Wrangel and Herald Islands currently support 68% of denning Chukchi Sea polar bears (this study; total proportion of all dens), and although Wrangel Island is designated as a Nature Reserve, it is currently being developed with both tourism and military infrastructure on some parts of the island. The viability of these habitats for denning may depend on human activities and continuation of the ice retreat date occurring after June 24. Our study is not the first to document changes in the distribution of den sites on land due to changes in the timing of sea ice melt and freeze [31]. Derocher et al. [31] found that in years with late ice advancement in autumn, polar bears in Svalbard, Norway, were much less likely to make use of a specific island to den than when ice was present earlier in autumn. Polar bears have been regularly observed onshore in Chukotka in recent years despite none of the bears in our study spending more than 7 days onshore there between August and October. This could indicate that bears are using the Chukotka coast in the summer less than they have in the past because ice retreat occurs earlier there (15% of bears captured 1986–1995 summered in Chukotka). Nineteen percent of our sample of bears that denned on land denned on Chukotka. We did observe in our study that some bears summered on one land mass and moved to a different land mass for denning, including one bear that came onshore in Alaska, and one that came onshore on Wrangel Island, and both then later denned on the Chukotkan coast. The nutritional implications of increased land use by polar bears in the Chukchi Sea depend on the energetic balance they achieve throughout the year. Given that polar bears primarily consume energetically-rich marine prey on the sea ice, long-term declines in sea ice availability have been associated with negative energetic and eventually demographic implications [24, 61, 62], although such relationships can be expected to exhibit spatial, temporal, and ecological complexity [34]. Increased durations on land can only be accommodated if bears come onshore in sufficient body condition to withstand longer periods of food deprivation, or obtain increased access to food while on shore. Although bears can reduce energetic requirements by reducing activity on land, they cannot reduce metabolic energy costs [28, 29]. On Wrangel Island, polar bears are reported to spend the majority of their time resting or walking [57]. But polar bears also congregate at walrus haulouts that occur on Wrangel Island and on the Chukotkan coast in some years, where they have been documented to hunt and kill live walruses as well as trigger stampedes and scavenge the resulting carcasses. These walrus haulouts occur every year in some locations on the Chukotkan coast, but in other areas, such as Wrangel Island, do not occur every year [58]. Polar bears in this region have also been observed opportunistically feeding on marine and terrestrial mammal carcasses, such as grey whales (Eschrichtius robustus), muskox (Ovibos moschatus) and reindeer (Rangifer tarandus)[58]. During peak lemming (Synaptomys spp.) years, polar bears actively dig them out of burrows as do grizzly bears [58]. Although land-based foraging on marine mammals and terrestrial plants and animals has been documented for Chukchi Sea polar bears, the contribution of terrestrial foods to energetic requirements during the onshore period is unknown. The lack of a change in the body condition and reproduction of Chukchi Sea polar bears during the time period of this study [34, 59] suggest that Chukchi Sea polar bears either come onshore with sufficient body fat or they are finding sufficient food resources on land (marine or terrestrial) to offset increased durations on land. The duration in which Chukchi Sea polar bears can remain on land before a threshold is reached and nutritional consequences occur is unknown. Timing of departure of polar bears from land post-denning in this study (April 7 in 1986–1995 and March 27 in 2008–2013) was consistent with previous studies which reported den emergence occurring between March 20th to April 5th with bears spending an additional 0.5 to 7 days at the den site [63]. In recent years, 3 bears departed from dens prior to March 1 which could suggest that some bears may be compensating for coming to land earlier by departing early; however, whether these maternity denning attempts were successful is unknown. Bears could continue to successfully produce cubs by shifting the denning period earlier, or conversely they could be abandoning dens early if energetic reserves run out. Bears that emerge early from dens could face low food availability and colder temperatures. The increased frequency and duration of land use we observed in this study and its association with reduced summer sea ice suggest that Chukchi Sea polar bears may continue to increase the amount of time they spend on land in the future. This pattern of increased land use in response to sea ice loss is consistent with other studies [5, 7, 9] and predictions for polar bears [17, 32] suggesting that one of the primary ways in which polar bears will respond to continued sea ice loss will be via increased use of Arctic coastal terrestrial habitats. This has implications for both human safety and polar bear conservation and, in light of increased human activity in the Arctic, will likely require proactive management to maintain terrestrial habitats for polar bears and mitigate human-polar bear interactions [17].

Supporting Information S1 Table. Akaike Information Criteria (AIC) scores from conditional logistic regression examining the effects of time period (1986–1995 versus 2008–2013), and ice retreat or return date on the location of summering and denning polar bears in the Chukchi Sea. https://doi.org/10.1371/journal.pone.0142213.s001 (DOCX)

Acknowledgments G. Garner (deceased) led capture efforts in the Chukchi Sea in the 1980s and 1990s. G. York, C. Perham, S. Miller, J. Wilder, T. DeBruyn, K. Titus, M. Lockhart, N. Ovsyanikov and many others assisted with data collection. The communities of Kotzebue and Point Hope, the Selawik National Wildlife Refuge, Teck Alaska, Inc, and NANA Regional Corporation provided significant field support. We thank our excellent pilots for ensuring safe capture operations. The findings and conclusions in this paper are those of the authors and do not necessarily represent the views of the USFWS. This paper was reviewed and approved by USGS under their Fundamental Science Practices policy (http://www.usgs.gov/fsp). Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

Author Contributions Conceived and designed the experiments: KDR EVR RRW MSM. Performed the experiments: KDR RRW MSM EVR. Analyzed the data: KDR RRW MSM JO DCD. Wrote the paper: KDR RRW DCD EVR.