Drivers of habitat availability for terrestrial mammals: Unravelling the role of livestock, land conversion and intrinsic traits in the past 50 years.

The global decline of terrestrial species is largely due to the degradation, loss and fragmentation of their habitats. The conversion of natural ecosystems for cropland, rangeland, forest products and human infrastructure are the primary causes of habitat deterioration. Due to the paucity of data on the past distribution of species and the scarcity of fine-scale habitat conversion maps, however, accurate assessment of the recent effects of habitat degradation, loss and fragmentation on the range of mammals has been near impossible. We aim to assess the proportions of available habitat within the lost and retained parts of mammals' distribution ranges, and to identify the drivers of habitat availability. We produced distribution maps for 475 terrestrial mammals for the range they occupied 50 years ago and compared them to current range maps. We then calculated the differences in the percentage of 'area of habitat' (habitat available to a species within its range) between the lost and retained range areas. Finally, we ran generalized linear mixed models to identify which variables were more influential in determining habitat availability in the lost and retained parts of the distribution ranges. We found that 59% of species had a lower proportion of available habitat in the lost range compared to the retained range, thus hypothesizing that habitat loss could have contributed to range declines. The most important factors negatively affecting habitat availability were the conversion of land to rangeland and high density of livestock. Significant intrinsic traits were those related to reproductive timing and output, habitat breadth and medium body size. Our findings emphasize the importance of implementing conservation strategies to mitigate the impacts caused by human activities on the habitats of mammals, and offer evidence indicating which species have the potential to reoccupy portions of their former range if other threats cease to occur.

Building use-inspired species distribution models: Using multiple data types to examine and improve model performance.

Species distribution models (SDMs) are becoming an important tool for marine conservation and management. Yet while there is an increasing diversity and volume of marine biodiversity data for training SDMs, little practical guidance is available on how to leverage distinct data types to build robust models. We explored the effect of different data types on the fit, performance and predictive ability of SDMs by comparing models trained with four data types for a heavily exploited pelagic fish, the blue shark (Prionace glauca), in the Northwest Atlantic: two fishery dependent (conventional mark-recapture tags, fisheries observer records) and two fishery independent (satellite-linked electronic tags, pop-up archival tags). We found that all four data types can result in robust models, but differences among spatial predictions highlighted the need to consider ecological realism in model selection and interpretation regardless of data type. Differences among models were primarily attributed to biases in how each data type, and the associated representation of absences, sampled the environment and summarized the resulting species distributions. Outputs from model ensembles and a model trained on all pooled data both proved effective for combining inferences across data types and provided more ecologically realistic predictions than individual models. Our results provide valuable guidance for practitioners developing SDMs. With increasing access to diverse data sources, future work should further develop truly integrative modeling approaches that can explicitly leverage the strengths of individual data types while statistically accounting for limitations, such as sampling biases.

A landscape-scale framework to identify refugia from multiple stressors.

From a conservation perspective, quantifying potential refugial capacity has been predominantly focused on climate refugia, which is critical for maintaining the persistence of species and ecosystems. However, protection from other stressors, such as human-induced changes in fire and hydrology, that cause habitat loss, degradation, and fragmentation is also necessary to ensure that conservation efforts focused on climate are not undermined by other threats. Thus, conceptual and methodological advances for quantifying potential refugia from multiple anthropogenic stressors are important to support conservation efforts. We devised a new conceptual approach, the domains of refugia, for assessing refugial capacity that identifies areas where exposure to multiple stressors is low. In our framework, patterns of environmental variability (e.g., increased frequency of warm summers), thresholds of resilience, and extent and intensity of stressors are used to identify areas of potential refugia from a suite of ongoing anthropogenic stressors (e.g., changes in fire regime). To demonstrate its utility, we applied the framework to a Southern California landscape. Sites with high refugial capacity (super-refugia sites) had on average 30% fewer extremely warm summers, 20% fewer fire events, 10% less exposure to altered river channels and riparian areas, and 50% fewer recreational trails than the surrounding landscape. Our results suggest that super-refugia sites (∼8200 km2 ) for some natural communities are underrepresented in the existing protected area network, a finding that can inform efforts to expand protected areas. Our case study highlights how considering exposure to multiple stressors can inform planning and practice to conserve biodiversity in a changing world.

Marco Conceptual a para Identificar Refugios de Múltiples Amenazas a Escala de Paisaje Resumen Desde la perspectiva de la conservación, la cuantificación de la capacidad potencial de refugio se ha enfocado principalmente en los refugios climáticos, los cuales son críticos para mantener la persistencia de las especies y los ecosistemas. Sin embargo, la protección ante otros factores estresantes, como los cambios inducidos por los humanos en los incendios y la hidrología, que causan la pérdida, degradación y fragmentación del hábitat, también son necesarios para asegurar que los esfuerzos de conservación enfocados en el clima no sean afectados por otras amenazas. Por lo tanto, los avances conceptuales y metodológicos para cuantificar los refugios potenciales ante múltiples factores estresantes causados por el humano son importantes para asegurar que los esfuerzos de conservación logren sus objetivos. Diseñamos una nueva estrategia conceptual, los dominios de los refugios, para evaluar la capacidad de refugio de un paisaje donde la exposición a múltiples factores estresantes es baja. En nuestro marco conceptual usamos los patrones de variabilidad ambiental (p. ej.: incremento en la frecuencia de veranos cálidos), los umbrales de resiliencia y la extensión e intensidad de los factores estresantes para identificar las áreas de refugios potenciales a partir de un conjunto de factores antropogénicos persistentes (p. ej.: cambios en el régimen de incendios). Para demostrar su utilidad, aplicamos el marco conceptual a un paisaje del sur de California. Los sitios con una alta capacidad de refugio (sitios de súper-refugios) tuvieron en promedio un 30% menos veranos extremadamente cálidos, 20% menos eventos de incendios y 50% menos senderos recreativos que el paisaje circundante. Nuestros resultados sugieren que los sitios de súper-refugios (∼ 8,200 km2 ) para algunas comunidades naturales están subrepresentados en la red existente de áreas protegidas, un resultado que puede orientar los esfuerzos por expandir las áreas protegidas. Nuestro estudio de caso resalta que considerar la exposición a múltiples amenazas puede guiar la planificación y la práctica de la conservación de la biodiversidad en un mundo cambiante.

Evaluation of thermal infrared imaging from uninhabited aerial vehicles for arboreal wildlife surveillance.

An important component of wildlife management and conservation is monitoring the health and population size of wildlife species. Monitoring the population size of an animal group can inform researchers of habitat use, potential changes in habitat and resulting behavioral adaptations, individual health, and the effectiveness of conservation efforts. Arboreal monkeys are difficult to monitor as their habitat is often poorly accessible and most monkey species have some degree of camouflage, making them hard to observe in and below the tree canopy. Surveys conducted using uninhabited aerial vehicles (UAVs) equipped with thermal infrared (TIR) cameras can help overcome these limitations by flying above the canopy and using the contrast between the warm body temperature of the monkeys and the cooler background vegetation, reducing issues with impassable terrain and animal camouflage. We evaluated the technical and procedural elements associated with conducting UAV-TIR surveys for arboreal and terrestrial macaque species. Primary imaging missions and analyses were conducted over a monkey park housing approximately 160 semi-free-ranging Japanese macaques (Macaca fuscata). We demonstrate Repeat Station Imaging (RSI) procedures using co-registered TIR image pairs facilitate the use of image differencing to detect targets that were moving during rapid sequence imaging passes. We also show that 3D point clouds may be generated from highly overlapping UAV-TIR image sets in a forested setting using structure from motion (SfM) image processing techniques. A point cloud showing area-wide elevation values was generated from TIR imagery, but it lacked sufficient point density to reliably determine the 3D locations of monkeys.

Species and population specific gene expression in blood transcriptomes of marine turtles.

BACKGROUND: Transcriptomic data has demonstrated utility to advance the study of physiological diversity and organisms' responses to environmental stressors. However, a lack of genomic resources and challenges associated with collecting high-quality RNA can limit its application for many wild populations. Minimally invasive blood sampling combined with de novo transcriptomic approaches has great potential to alleviate these barriers. Here, we advance these goals for marine turtles by generating high quality de novo blood transcriptome assemblies to characterize functional diversity and compare global transcriptional profiles between tissues, species, and foraging aggregations. RESULTS: We generated high quality blood transcriptome assemblies for hawksbill (Eretmochelys imbricata), loggerhead (Caretta caretta), green (Chelonia mydas), and leatherback (Dermochelys coriacea) turtles. The functional diversity in assembled blood transcriptomes was comparable to those from more traditionally sampled tissues. A total of 31.3% of orthogroups identified were present in all four species, representing a core set of conserved genes expressed in blood and shared across marine turtle species. We observed strong species-specific expression of these genes, as well as distinct transcriptomic profiles between green turtle foraging aggregations that inhabit areas of greater or lesser anthropogenic disturbance. CONCLUSIONS: Obtaining global gene expression data through non-lethal, minimally invasive sampling can greatly expand the applications of RNA-sequencing in protected long-lived species such as marine turtles. The distinct differences in gene expression signatures between species and foraging aggregations provide insight into the functional genomics underlying the diversity in this ancient vertebrate lineage. The transcriptomic resources generated here can be used in further studies examining the evolutionary ecology and anthropogenic impacts on marine turtles.

Fisheries bycatch risk to marine megafauna is intensified in Lagrangian coherent structures.

Incidental catch of nontarget species (bycatch) is a major barrier to ecological and economic sustainability in marine capture fisheries. Key to mitigating bycatch is an understanding of the habitat requirements of target and nontarget species and the influence of heterogeneity and variability in the dynamic marine environment. While patterns of overlap among marine capture fisheries and habitats of a taxonomically diverse range of marine vertebrates have been reported, a mechanistic understanding of the real-time physical drivers of bycatch events is lacking. Moving from describing patterns toward understanding processes, we apply a Lagrangian analysis to a high-resolution ocean model output to elucidate the fundamental mechanisms that drive fisheries interactions. We find that the likelihood of marine megafauna bycatch is intensified in attracting Lagrangian coherent structures associated with submesoscale and mesoscale filaments, fronts, and eddies. These results highlight how the real-time tracking of dynamic structures in the oceans can support fisheries sustainability and advance ecosystem-based management.

Rethinking megafauna.

Concern for megafauna is increasing among scientists and non-scientists. Many studies have emphasized that megafauna play prominent ecological roles and provide important ecosystem services to humanity. But, what precisely are 'megafauna'? Here, we critically assess the concept of megafauna and propose a goal-oriented framework for megafaunal research. First, we review definitions of megafauna and analyse associated terminology in the scientific literature. Second, we conduct a survey among ecologists and palaeontologists to assess the species traits used to identify and define megafauna. Our review indicates that definitions are highly dependent on the study ecosystem and research question, and primarily rely on ad hoc size-related criteria. Our survey suggests that body size is crucial, but not necessarily sufficient, for addressing the different applications of the term megafauna. Thus, after discussing the pros and cons of existing definitions, we propose an additional approach by defining two function-oriented megafaunal concepts: 'keystone megafauna' and 'functional megafauna', with its variant 'apex megafauna'. Assessing megafauna from a functional perspective could challenge the perception that there may not be a unifying definition of megafauna that can be applied to all eco-evolutionary narratives. In addition, using functional definitions of megafauna could be especially conducive to cross-disciplinary understanding and cooperation, improvement of conservation policy and practice, and strengthening of public perception. As megafaunal research advances, we encourage scientists to unambiguously define how they use the term 'megafauna' and to present the logic underpinning their definition.

Contaminant Exposure Linked to Cellular and Endocrine Biomarkers in Southern California Bottlenose Dolphins.

Cetaceans in the Southern California Bight (SCB) are exposed to high levels of halogenated organic contaminants (HOCs), which have previously been linked to impaired reproductive health and immune responses. We used a combination of molecular tools to examine the potential physiological impacts of HOC exposure in two bottlenose dolphin ( Tursiops truncatus) ecotypes in the SCB. We quantified 25 HOCs in the blubber of 22 biopsies collected from males between 2012 and 2016. We then analyzed genome-wide gene expression in skin using RNA-sequencing and measured blubber testosterone to compare HOC exposure with cellular and endocrine biomarkers. We found high levels of HOCs in both ecotypes with significantly higher total polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), tris(4-chlorophenyl)methanol (TCPMOH), and chlordane-related compounds in the coastal ecotype versus the offshore ecotype. We found evidence of PBDE bioaccumulation in both ecotypes, however, the pattern of bioaccumulation or endocrine disruption for other HOCs was different between the ecotypes, suggesting potential endocrine disruption in the coastal ecotype. We also observed correlations between HOCs and gene coexpression networks enriched for xenobiotic metabolism, hormone metabolism, and immune response that could indicate cellular effects from HOC exposure. By integrating measurements of HOC load with both transcriptome profiling and endocrine biomarkers, our approach provides insight into HOC exposure and potential impacts on wild cetacean health in southern California.

Comprehensive Screening Links Halogenated Organic Compounds with Testosterone Levels in Male Delphinus delphis from the Southern California Bight.

While environmental pollutants have been associated with changes in endocrine health in cetaceans, efforts to link contaminant exposure with hormones have largely been limited to a list of known, targeted contaminants, overlooking minimally characterized or unknown compounds of emerging concern. To address this gap, we analyzed a suite of potential endocrine disrupting halogenated organic compounds (HOCs) in blubber from 16 male short-beaked common dolphins ( Delphinus delphis) with known maturity status collected from fishery bycatch in the Southern California Bight. We employed a suspect screening mass spectrometry-based method to investigate a wide range of HOCs that were previously observed in cetaceans from the same region. Potential endocrine effects were assessed through the measurement of blubber testosterone. We detected 167 HOCs, including 81 with known anthropogenic sources, 49 of unknown origin, and 37 with known natural sources. The sum of 11 anthropogenic and 4 unknown HOC classes were negatively correlated with blubber testosterone. Evidence suggests that elevated anthropogenic HOC load contributes to impaired testosterone production in mature male D. delphis. The application of this integrative analytical approach to cetacean contaminant analysis allows for inference of the biological consequences of accumulation of HOCs and prioritization of compounds for future environmental toxicology research.

Fit to predict? Eco-informatics for predicting the catchability of a pelagic fish in near real time.

The ocean is a dynamic environment inhabited by a diverse array of highly migratory species, many of which are under direct exploitation in targeted fisheries. The timescales of variability in the marine realm coupled with the extreme mobility of ocean-wandering species such as tuna and billfish complicates fisheries management. Developing eco-informatics solutions that allow for near real-time prediction of the distributions of highly mobile marine species is an important step towards the maturation of dynamic ocean management and ecological forecasting. Using 25 yr (1990-2014) of NOAA fisheries' observer data from the California drift gillnet fishery, we model relative probability of occurrence (presence-absence) and catchability (total catch per gillnet set) of broadbill swordfish Xiphias gladius in the California Current System. Using freely available environmental data sets and open source software, we explore the physical drivers of regional swordfish distribution. Comparing models built upon remotely sensed data sets with those built upon a data-assimilative configuration of the Regional Ocean Modelling System (ROMS), we explore trade-offs in model construction, and address how physical data can affect predictive performance and operational capacity. Swordfish catchability was found to be highest in deeper waters (>1,500 m) with surface temperatures in the 14-20°C range, isothermal layer depth (ILD) of 20-40 m, positive sea surface height (SSH) anomalies, and during the new moon (<20% lunar illumination). We observed a greater influence of mesoscale variability (SSH, wind speed, isothermal layer depth, eddy kinetic energy) in driving swordfish catchability (total catch) than was evident in predicting the relative probability of presence (presence-absence), confirming the utility of generating spatiotemporally dynamic predictions. Data-assimilative ROMS circumvent the limitations of satellite remote sensing in providing physical data fields for species distribution models (e.g., cloud cover, variable resolution, subsurface data), and facilitate broad-scale prediction of dynamic species distributions in near real time.

Global patterns of marine mammal, seabird, and sea turtle bycatch reveal taxa-specific and cumulative megafauna hotspots.

Recent research on ocean health has found large predator abundance to be a key element of ocean condition. Fisheries can impact large predator abundance directly through targeted capture and indirectly through incidental capture of nontarget species or bycatch. However, measures of the global nature of bycatch are lacking for air-breathing megafauna. We fill this knowledge gap and present a synoptic global assessment of the distribution and intensity of bycatch of seabirds, marine mammals, and sea turtles based on empirical data from the three most commonly used types of fishing gears worldwide. We identify taxa-specific hotspots of bycatch intensity and find evidence of cumulative impacts across fishing fleets and gears. This global map of bycatch illustrates where data are particularly scarce--in coastal and small-scale fisheries and ocean regions that support developed industrial fisheries and millions of small-scale fishers--and identifies fishing areas where, given the evidence of cumulative hotspots across gear and taxa, traditional species or gear-specific bycatch management and mitigation efforts may be necessary but not sufficient. Given the global distribution of bycatch and the mitigation success achieved by some fleets, the reduction of air-breathing megafauna bycatch is both an urgent and achievable conservation priority.

Improved estimation of intrinsic growth r(max) for long-lived species: integrating matrix models and allometry.

Intrinsic population growth rate (r(max)) is an important parameter for many ecological applications, such as population risk assessment and harvest management. However, r(max) can be a difficult parameter to estimate, particularly for long-lived species, for which appropriate life table data or abundance time series are typically not obtainable. We describe a method for improving estimates of r(max) for long-lived species by integrating life-history theory (allometric models) and population-specific demographic data (life table models). Broad allometric relationships, such as those between life history traits and body size, have long been recognized by ecologists. These relationships are useful for deriving theoretical expectations for r(max), but r(max) for real populations may vary from simple allometric estimators for "archetypical" species of a given taxa or body mass. Meanwhile, life table approaches can provide population-specific estimates of r(max) from empirical data, but these may have poor precision from imprecise and missing vital rate parameter estimates. Our method borrows strength from both approaches to provide estimates that are consistent with both life-history theory and population-specific empirical data, and are likely to be more robust than estimates provided by either method alone. Our method uses an' allometric constant: the product of r(max) and the associated generation time for a stable-age population growing at this rate. We conducted a meta-analysis to estimate the mean and variance of this allometric constant across well-studied populations from three vertebrate taxa (birds, mammals, and elasmobranchs) and found that the mean was approximately 1.0 for each taxon. We used these as informative Bayesian priors that determine how much to "shrink" imprecise vital rate estimates for a data-limited population toward the allometric expectation. The approach ultimately provides estimates of r(max) (and other vital rates) that reflect a balance of information from the individual studied population, theoretical expectation, and meta-analysis of other populations. We applied the method specifically to an archetypical petrel (representing the genus Procellaria) and to white sharks (Carcharodon carcharias) in the context of estimating sustainable-fishery bycatch limits.

Using community-level metrics to monitor the effects of marine protected areas on biodiversity.

Marine protected areas (MPAs) are used to protect species, communities, and their associated habitats, among other goals. Measuring MPA efficacy can be challenging, however, particularly when considering responses at the community level. We gathered 36 abundance and 14 biomass data sets on fish assemblages and used meta-analysis to evaluate the ability of 22 distinct community diversity metrics to detect differences in community structure between MPAs and nearby control sites. We also considered the effects of 6 covariates-MPA size and age, MPA size and age interaction, latitude, total species richness, and level of protection-on each metric. Some common metrics, such as species richness and Shannon diversity, did not differ consistently between MPA and control sites, whereas other metrics, such as total abundance and biomass, were consistently different across studies. Metric responses derived from the biomass data sets were more consistent than those based on the abundance data sets, suggesting that community-level biomass differs more predictably than abundance between MPA and control sites. Covariate analyses indicated that level of protection, latitude, MPA size, and the interaction between MPA size and age affect metric performance. These results highlight a handful of metrics, several of which are little known, that could be used to meet the increasing demand for community-level indicators of MPA effectiveness.

Global hidden spillover effects among concurrent green initiatives.

Concurrently implemented green initiatives to combat global environmental crises may be curtailed or even sacrificed given the ongoing global economic contraction. We collected empirical data and information about green initiatives from 15 sites or countries worldwide. We systematically explored how specific policy, intended behaviors, and gains of given green initiative may interact with those of other green initiatives concurrently implemented in the same geographic area or involving the same recipients. Surprisingly, we found that spillover effects were very divergent: one initiative could reduce the gain of another by 22 % âˆ¼ 100 %, representing alarming losses, while in other instances, substantial co-benefits could arise as one initiative can increase the gain of another by 9 % âˆ¼ 310 %. Leveraging these effects will help countries keep green initiatives with significant co-benefits but stop initiatives with substantial spillover losses in the face of widespread budget cuts, better meeting the United Nations' sustainable development goals.

Widespread habitat loss and redistribution of marine top predators in a changing ocean.

The Northwest Atlantic Ocean and Gulf of Mexico are among the fastest warming ocean regions, a trend that is expected to continue through this century with far-reaching implications for marine ecosystems. We examine the distribution of 12 highly migratory top predator species using predictive models and project expected habitat changes using downscaled climate models. Our models predict widespread losses of suitable habitat for most species, concurrent with substantial northward displacement of core habitats >500 km. These changes include up to >70% loss of suitable habitat area for some commercially and ecologically important species. We also identify predicted hot spots of multi-species habitat loss focused offshore of the U.S. Southeast and Mid-Atlantic coasts. For several species, the predicted changes are already underway, which are likely to have substantial impacts on the efficacy of static regulatory frameworks used to manage highly migratory species. The ongoing and projected effects of climate change highlight the urgent need to adaptively and proactively manage dynamic marine ecosystems.

Shifting the life-history paradigm: discovery of novel habitat use by hawksbill turtles.

Adult hawksbill turtles (Eretmochelys imbricata) are typically described as open-coast, coral reef and hard substrate dwellers. Here, we report new satellite tracking data on female hawksbills from several countries in the eastern Pacific that revealed previously undocumented behaviour for adults of the species. In contrast to patterns of habitat use exhibited by their Caribbean and Indo-Pacific counterparts, eastern Pacific hawksbills generally occupied inshore estuaries, wherein they had strong associations with mangrove saltwater forests. The use of inshore habitats and affinities with mangrove saltwater forests presents a previously unknown life-history paradigm for adult hawksbill turtles and suggests a potentially unique evolutionary trajectory for the species. Our findings highlight the variability in life-history strategies that marine turtles and other wide-ranging marine wildlife may exhibit among ocean regions, and the importance of understanding such disparities from an ecological and management perspective.

Dynamic habitat models: using telemetry data to project fisheries bycatch.

Fisheries bycatch is a recognized threat to marine megafauna. Addressing bycatch of pelagic species however is challenging owing to the dynamic nature of marine environments and vagility of these organisms. In order to assess the potential for species to overlap with fisheries, we propose applying dynamic habitat models to determine relative probabilities of species occurrence for specific oceanographic conditions. We demonstrate this approach by modelling habitats for Laysan (Phoebastria immutabilis) and black-footed albatrosses (Phoebastria nigripes) using telemetry data and relating their occurrence probabilities to observations of Hawaii-based longline fisheries in 1997-2000. We found that modelled habitat preference probabilities of black-footed albatrosses were high within some areas of the fishing range of the Hawaiian fleet and such preferences were important in explaining bycatch occurrence. Conversely, modelled habitats of Laysan albatrosses overlapped little with Hawaii-based longline fisheries and did little to explain the bycatch of this species. Estimated patterns of albatross habitat overlap with the Hawaiian fleet corresponded to bycatch observations: black-footed albatrosses were more frequently caught in this fishery despite being 10 times less abundant than Laysan albatrosses. This case study demonstrates that dynamic habitat models based on telemetry data may help to project interactions with pelagic animals relative to environmental features and that such an approach can serve as a tool to guide conservation and management decisions.

Landscape Seroprevalence of Three Hemorrhagic Disease-Causing Viruses in a Wild Cervid.

Disease plays a major role in shaping wildlife populations worldwide, and changes in landscape conditions can significantly influence risk of pathogen exposure, a threat to vulnerable wild species. Three viruses that cause hemorrhagic disease affect cervid populations in the USA (Odocoileus hemionus adenovirus, bluetongue virus, and epizootic hemorrhagic disease virus), but little is known of their distribution and prevalence in wild populations. We explored the distribution and co-occurrence of seroprevalence of these three pathogens in southern mule deer (Odocoileus hemionus fuliginatus), a subspecies of conservation concern and a harvested species native to southern California, to evaluate the distribution of exposure to these pathogens relative to landscape attributes. We found that habitat type, level of development, and proximity to livestock may affect hemorrhagic disease seroprevalence in southern mule deer. Continued monitoring of hemorrhagic disease-causing viruses in areas where deer are in proximity to cattle and human development is needed to better understand the implications of future outbreaks in wild populations and to identify opportunities to mitigate disease impacts in southern mule deer and other cervid species.

Foraging in marine habitats increases mercury concentrations in a generalist seabird.

Methylmercury concentrations vary widely across geographic space and among habitat types, with marine and aquatic-feeding organisms typically exhibiting higher mercury concentrations than terrestrial-feeding organisms. However, there are few model organisms to directly compare mercury concentrations as a result of foraging in marine, estuarine, or terrestrial food webs. The ecological impacts of differential foraging may be especially important for generalist species that exhibit high plasticity in foraging habitats, locations, or diet. Here, we investigate whether foraging habitat, sex, or fidelity to a foraging area impact blood mercury concentrations in western gulls (Larus occidentalis) from three colonies on the US west coast. Cluster analyses showed that nearly 70% of western gulls foraged primarily in ocean or coastal habitats, whereas the remaining gulls foraged in terrestrial and freshwater habitats. Gulls that foraged in ocean or coastal habitats for half or more of their foraging locations had 55% higher mercury concentrations than gulls that forage in freshwater and terrestrial habitats. Ocean-foraging gulls also had lower fidelity to a specific foraging area than freshwater and terrestrial-foraging gulls, but fidelity and sex were unrelated to gull blood mercury concentrations in all models. These findings support existing research that has described elevated mercury levels in species using aquatic habitats. Our analyses also demonstrate that gulls can be used to detect differences in contaminant exposure over broad geographic scales and across coarse habitat types, a factor that may influence gull health and persistence of other populations that forage across the land-sea gradient.