Distance measures and optimization spaces in quantitative fatty acid signature analysis

Quantitative fatty acid signature analysis has become an important method of diet estimation in ecology, especially marine ecology. Controlled feeding trials to validate the method and estimate the calibration coefficients necessary to account for differential metabolism of individual fatty acids have been conducted with several species from diverse taxa. However, research into potential refinements of the estimation method has been limited. We compared the performance of the original method of estimating diet composition with that of five variants based on different combinations of distance measures and calibration‐coefficient transformations between prey and predator fatty acid signature spaces. Fatty acid signatures of pseudopredators were constructed using known diet mixtures of two prey data sets previously used to estimate the diets of polar bears Ursus maritimus and gray seals Halichoerus grypus, and their diets were then estimated using all six variants. In addition, previously published diets of Chukchi Sea polar bears were re‐estimated using all six methods. Our findings reveal that the selection of an estimation method can meaningfully influence estimates of diet composition. Among the pseudopredator results, which allowed evaluation of bias and precision, differences in estimator performance were rarely large, and no one estimator was universally preferred, although estimators based on the Aitchison distance measure tended to have modestly superior properties compared to estimators based on the Kullback–Leibler distance measure. However, greater differences were observed among estimated polar bear diets, most likely due to differential estimator sensitivity to assumption violations. Our results, particularly the polar bear example, suggest that additional research into estimator performance and model diagnostics is warranted.     

Survival and abundance of polar bears in Alaska’s Beaufort Sea, 2001–2016

The Arctic Ocean is undergoing rapid transformation toward a seasonally ice‐free ecosystem. As ice‐adapted apex predators, polar bears (Ursus maritimus) are challenged to cope with ongoing habitat degradation and changes in their prey base driven by food‐web response to climate warming. Knowledge of polar bear response to environmental change is necessary to understand ecosystem dynamics and inform conservation decisions. In the southern Beaufort Sea (SBS) of Alaska and western Canada, sea ice extent has declined since satellite observations began in 1979 and available evidence suggests that the carrying capacity of the SBS for polar bears has trended lower for nearly two decades. In this study, we investigated the population dynamics of polar bears in Alaska''s SBS from 2001 to 2016 using a multistate Cormack–Jolly–Seber mark–recapture model. States were defined as geographic regions, and we used location data from mark–recapture observations and satellite‐telemetered bears to model transitions between states and thereby explain heterogeneity in recapture probabilities. Our results corroborate prior findings that the SBS subpopulation experienced low survival from 2003 to 2006. Survival improved modestly from 2006 to 2008 and afterward rebounded to comparatively high levels for the remainder of the study, except in 2012. Abundance moved in concert with survival throughout the study period, declining substantially from 2003 and 2006 and afterward fluctuating with lower variation around an average of 565 bears (95% Bayesian credible interval [340, 920]) through 2015. Even though abundance was comparatively stable and without sustained trend from 2006 to 2015, polar bears in the Alaska SBS were less abundant over that period than at any time since passage of the U.S. Marine Mammal Protection Act. The potential for recovery is likely limited by the degree of habitat degradation the subpopulation has experienced, and future reductions in carrying capacity are expected given current projections for continued climate warming.     

Habitat degradation affects the summer activity of polar bears

Understanding behavioral responses of species to environmental change is critical to forecasting population-level effects. Although climate change is significantly impacting species’ distributions, few studies have examined associated changes in behavior. Polar bear (Ursus maritimus) subpopulations have varied in their near-term responses to sea ice decline. We examined behavioral responses of two adjacent subpopulations to changes in habitat availability during the annual sea ice minimum using activity data. Location and activity sensor data collected from 1989 to 2014 for 202 adult female polar bears in the Southern Beaufort Sea (SB) and Chukchi Sea (CS) subpopulations were used to compare activity in three habitat types varying in prey availability: (1) land; (2) ice over shallow, biologically productive waters; and (3) ice over deeper, less productive waters. Bears varied activity across and within habitats with the highest activity at 50–75% sea ice concentration over shallow waters. On land, SB bears exhibited variable but relatively high activity associated with the use of subsistence-harvested bowhead whale carcasses, whereas CS bears exhibited low activity consistent with minimal feeding. Both subpopulations had fewer observations in their preferred shallow-water sea ice habitats in recent years, corresponding with declines in availability of this substrate. The substantially higher use of marginal habitats by SB bears is an additional mechanism potentially explaining why this subpopulation has experienced negative effects of sea ice loss compared to the still-productive CS subpopulation. Variability in activity among, and within, habitats suggests that bears alter their behavior in response to habitat conditions, presumably in an attempt to balance prey availability with energy costs.     

Simultaneous estimation of diet composition and calibration coefficients with fatty acid signature data

Knowledge of animal diets provides essential insights into their life history and ecology, although diet estimation is challenging and remains an active area of research. Quantitative fatty acid signature analysis (QFASA) has become a popular method of estimating diet composition, especially for marine species. A primary assumption of QFASA is that constants called calibration coefficients, which account for the differential metabolism of individual fatty acids, are known. In practice, however, calibration coefficients are not known, but rather have been estimated in feeding trials with captive animals of a limited number of model species. The impossibility of verifying the accuracy of feeding trial derived calibration coefficients to estimate the diets of wild animals is a foundational problem with QFASA that has generated considerable criticism. We present a new model that allows simultaneous estimation of diet composition and calibration coefficients based only on fatty acid signature samples from wild predators and potential prey. Our model performed almost flawlessly in four tests with constructed examples, estimating both diet proportions and calibration coefficients with essentially no error. We also applied the model to data from Chukchi Sea polar bears, obtaining diet estimates that were more diverse than estimates conditioned on feeding trial calibration coefficients. Our model avoids bias in diet estimates caused by conditioning on inaccurate calibration coefficients, invalidates the primary criticism of QFASA, eliminates the need to conduct feeding trials solely for diet estimation, and consequently expands the utility of fatty acid data to investigate aspects of ecology linked to animal diets.     

bels: backward elimination locus selection for studies of mixture composition or individual assignment

Methods of evaluating loci in studies of mixture composition or individual assignment are largely based on performance characteristics of individual loci. Synergisms between loci are not exploited. Loci are often evaluated based on their ability to resolve individual populations, even though multipopulation aggregations are more commonly of interest. In addition, measures of locus performance may indirectly relate to investigative objectives. A new computer program, bels, offers an alternative that addresses these limitations and may be preferable to existing methods in some applications. The algorithm is illustrated using Yukon River chum salmon (Oncorhynchus keta) data.     

TRENDS IN AGE STRUCTURE AND PRODUCTIVITY OF PACIFIC WALRUSES HARVESTED IN THE BERING STRAIT REGION OF ALASKA, 1952–2002

Pacific walruses ( Odobenus rosmarus divergens ) are harvested by subsistence hunters in Alaska as they migrate north through the Bering Strait in the spring. Harvest records and biological specimens have been collected from the Bering Strait communities of Little Diomede, Gambell, and Savoonga since the 1950s. Harvest levels in the Bering Strait region peaked in the late 1980s and declined thereafter; however, there was considerable variation in the size and composition of the harvests among communities and over time. The relationships among characteristics of the community harvests and the presence of temporal trends were investigated using generalized linear models. The proportion of females in the catch increased over time in all three communities, while the proportion pregnant among harvested females declined over the range of sample years. The ages of harvested walruses increased over time in all three communities through the 1980s, after which trends in age stabilized or began to decline. The age of first reproduction was significantly older for females sampled in 1975–1985 than for females sampled between 1952 and 1962 or 1992 and 1998. Factors thought to have influenced the size and composition of the catch over the past 50 yr include hunter preferences, harvest management regimes, environmental conditions, and changes in the population itself.     

Should fatty acid signature proportions sum to 1 for diet estimation?

Knowledge of predator diets, including how diets might change through time or differ among predators, provides essential insights into their ecology. Diet estimation therefore remains an active area of research within quantitative ecology. Quantitative fatty acid signature analysis (QFASA) is an increasingly common method of diet estimation. QFASA is based on a data library of prey signatures, which are vectors of proportions summarizing the fatty acid composition of lipids, and diet is estimated as the mixture of prey signatures that most closely approximates a predator’s signature. Diets are typically estimated using proportions from a subset of all fatty acids that are known to be solely or largely influenced by diet. Given the subset of fatty acids selected, the current practice is to scale their proportions to sum to 1.0. However, scaling signature proportions has the potential to distort the structural relationships within a prey library and between predators and prey. To investigate that possibility, we compared the practice of scaling proportions with two alternatives and found that the traditional scaling can meaningfully bias diet estimators under some conditions. Two aspects of the prey types that contributed to a predator’s diet influenced the magnitude of the bias: the degree to which the sums of unscaled proportions differed among prey types and the identifiability of prey types within the prey library. We caution investigators against the routine scaling of signature proportions in QFASA.     

Body size and condition of coyotes in southern Texas

Body size and condition of coyotes (Canis latrans) from a high-density population in Webb County, Texas (USA) were analyzed for age, sex and seasonal differences during 1980 to 1986. Mean body mass was progressively greater for juvenile, yearling and adult coyotes. Males were heavier and longer than females in each age class. Indices of intraperitoneal fat deposits were similar between sexes. Juveniles continued growth from fall to spring. Adults and yearlings both lost intraperitoneal fat overwinter. Mean body mass of adults decreased overwinter but mass of yearlings did not differ significantly between fall and spring. Territorial and transient female coyotes did not differ in mean body mass, length or indices of subcutaneous fat deposits.