Total mercury concentrations in Steller sea lion bone: Variability among locations and elements

Mercury is a global contaminant that bioaccumulates in a tissue-specific manner in long-lived predators such as Steller sea lions (SSL). Bone is a well-preserved material amenable for studying millennial scale trends; however, little is known about the distribution and variability of total mercury concentrations ([THg]) within individual bones and among bone elements in SSL. We assessed SSL bone [THg] variability with respect to physiologic age, bone type, longitudinally within a bone, and among bone elements. Pup bones (mean ± SD; 31.4 ± 13.58 ppb) had greater [THg] than adults (7.9 ± 1.91 ppb). There were greater and more variable [THg] within individual long bones near epiphyses compared to mid-diaphysis. Pup spongy bone in ribs (62.7 ± 44.79 ppb) had greater [THg] than long bones (23.5 ± 8.83 ppb) and phalanges (19.6 ± 10.78 ppb). These differences are likely due to variability in bone composition, growth, and turnover rate. This study informs standardized sampling procedures for [THg] in bone to improve interpretations of mercury variability over time and space.

A multistate Langevin diffusion for inferring behavior-specific habitat selection and utilization distributions

The identification of important habitat and the behavior(s) associated with it is critical to conservation and place-based management decisions. Behavior also links life-history requirements and habitat use, which are key to understanding why animals use certain habitats. Animal population studies often use tracking data to quantify space use and habitat selection, but they typically either ignore movement behavior (e.g., foraging, migrating, nesting) or adopt a two-stage approach that can induce bias and fail to propagate uncertainty. We develop a habitat-driven Langevin diffusion for animals that exhibit distinct movement behavior states, thereby providing a novel single-stage statistical method for inferring behavior-specific habitat selection and utilization distributions in continuous time. Practitioners can customize, fit, assess, and simulate our integrated model using the provided R package. Simulation experiments demonstrated that the model worked well under a range of sampling scenarios as long as observations were of sufficient temporal resolution. Our simulations also demonstrated the importance of accounting for different behaviors and the misleading inferences that can result when these are ignored. We provide case studies using plains zebra (Equus quagga) and Steller sea lion (Eumetopias jubatus) telemetry data. In the zebra example, our model identified distinct "encamped" and "exploratory" states, where the encamped state was characterized by strong selection for grassland and avoidance of other vegetation types, which may represent selection for foraging resources. In the sea lion example, our model identified distinct movement behavior modes typically associated with this marine central-place forager and, unlike previous analyses, found foraging-type movements to be associated with steeper offshore slopes characteristic of the continental shelf, submarine canyons, and seamounts that are believed to enhance prey concentrations. This is the first single-stage approach for inferring behavior-specific habitat selection and utilization distributions from tracking data that can be readily implemented with user-friendly software. As certain behaviors are often more relevant to specific conservation or management objectives, practitioners can use our model to help inform the identification and prioritization of important habitats. Moreover, by linking individual-level movement behaviors to population-level spatial processes, the multistate Langevin diffusion can advance inferences at the intersection of population, movement, and landscape ecology.

Whiskers as a novel tissue for tracking reproductive and stress-related hormones in North Pacific otariid pinnipeds

Keratinized tissues, including whiskers, are ideal for acquiring a record of physiological parameters. Most tissues provide a snapshot of physiological status; however, whiskers may support longitudinal sampling for reproductive and stress-related hormones, if hormones are incorporated as whiskers grow and concentrations change with physiological state. Whiskers from female Steller sea lions (Eumetopias jubatus) and northern fur seals (Callorhinus ursinus) were serially sectioned and pulverized and steroid hormones were extracted. Standard methods were used to validate enzyme immunoassay kits for cortisol, progesterone, 17β-estradiol and testosterone. All hormones were measurable in whisker segments from both species with progesterone concentrations showing cyclical patterns, which appear to signify previous pregnancies or luteal phases. Yearly progesterone concentrations were greater in years a pup was produced compared with years when no pup was observed. Free-ranging female Steller sea lions had reproductive rates between 0 and 1.0 (0.53 ± 0.33, n = 12) using a yearly progesterone concentration of 30 pg/mg or greater to classify a reproductive year as producing a pup and below 30 pg/mg as non-reproductive. Cortisol concentrations were greater near the root and rapidly declined, lacking any obvious patterns, throughout the rest of the whisker. Progesterone and testosterone concentrations were able to help determine sex of unknown individuals. Immunohistochemistry revealed that steroid hormones most likely do not leach out of whiskers based on the deposition patterns of progesterone and cortisol being present throughout the whisker length. Overall, measuring steroid hormones in whiskers can reveal individual reproductive histories over multiple years in sea lions and fur seals. Cyclical patterns of δ15N were useful for identifying periods of up to ~10 years of growth within whiskers, and measuring both stable isotopes and hormones may be useful for differentiating periods of active gestation from diapause and potentially track multi-year reproductive histories of female otariids.

Investigating life-history traits of Steller sea lions with multistate hidden Markov mark-recapture models: Age at weaning and body size effects

The duration of offspring care is critical to female fitness and population resilience by allowing flexibility in life-history strategies in a variable environment. Yet, for many mammals capable of extended periods of maternal care, estimates of the duration of offspring dependency are not available and the relative importance of flexibility of this trait on fitness and population viability has rarely been examined. We used data from 4,447 Steller sea lions Eumetopias jubatus from the Gulf of Alaska and multistate hidden Markov mark-recapture models to estimate age-specific weaning probabilities. Maternal care beyond age 1 was common: Weaning was later for animals from Southeast Alaska (SEAK) and Prince William Sound (PWS, weaning probabilities: 0.536-0.648/0.784-0.873 by age 1/2) compared with animals born to the west (0.714-0.855/0.798-0.938). SEAK/PWS animals were also smaller than those born farther west, suggesting a possible link. Females weaned slightly earlier (+0.080 at age 1 and 2) compared with males in SEAK only. Poor survival for weaned versus unweaned yearlings occurred in southern SEAK (female survival probabilities: 0.609 vs. 0.792) and the central Gulf (0.667 vs. 0.901), suggesting poor conditions for juveniles in these areas. First-year survival increased with neonatal body mass (NBM) linearly in the Gulf and nonlinearly in SEAK. The probability of weaning at age 1 increased linearly with NBM for SEAK animals only. Rookeries where juveniles weaned at earlier ages had lower adult female survival, but age at weaning was unrelated to population trends. Our results suggest the time to weaning may be optimized for different habitats based on long-term average conditions (e.g., prey dynamics), that may also shape body size, with limited short-term plasticity. An apparent trade-off of adult survival in favor of juvenile survival and large offspring size in the endangered Gulf of Alaska population requires further study.

Regional variations and drivers of mercury and selenium concentrations in Steller sea lions

Mercury (Hg) can be neurotoxic to mammals and impact reproduction, whereas selenium (Se) is an important antioxidant known to ameliorate some adverse effects of Hg. Total Hg concentrations ([THg]) were measured in lanugo (pelage grown in utero) of 812 Steller sea lion (Eumetopias jubatus) pups across Alaska and Russia to assess fetal exposure during late gestation. The molar ratio of total Se to THg (TSe:THg) was determined in whole blood collected from 291 pups. Stable isotope ratios of carbon and nitrogen were measured in sections of vibrissae (whiskers, n = 498) and in lanugo (n = 480) of pups grown during late gestation to track diet variations among adult females that can drive Hg and Se exposure during this critical fetal development period. Lanugo [THg] ranged from 1.4 to 73.7 μg/g dry weight with the lowest median [THg] in Southeast Alaska. Pups from the Western Aleutian Islands had higher median lanugo [THg] than pups from other metapopulations in Alaska. Over 25% of pups in the Western Aleutian Islands had [THg] above published risk thresholds (20 μg/g) for other mammals. Whole blood molar TSe:THg was significantly lower in the Western Aleutian Islands and in some parts of the Central Aleutian Islands with higher molar ratios found in the Eastern Aleutian Islands and Central Gulf of Alaska. This suggests a limitation on potential protective functions of Se in the western regions with the highest relative [THg]. The Central Aleutian Island pups with [THg] over 20 μg/g had higher δ¹⁵N ratios than pups with lower [THg] suggesting dams consuming higher trophic level prey is a key driver for Hg exposure. However, regional differences likely reflect variability in diet of the dam during gestation and in Hg food web dynamics between oceanic regimes east and west of key passes in the Aleutian Islands.