Trace me if you can: the use of intrinsic biogeochemical markers in marine top predators

Human activities have serious impacts on marine apex predators. Inadequate knowledge of the spatial and trophic ecology of these marine animals ultimately compromises the viability of their populations and impedes our ability to use them as environmental biomonitors. Intrinsic biogeochemical markers, such as stable isotopes, fatty acids, trace elements, and chemical pollutants, are increasingly being used to trace the spatial and trophic ecology of marine top predators. Notable advances include the emergence of the first oceanographic"isoscapes" (isotopic geographic gradients), the advent of compound-specific isotopic analyses, improvements in diet reconstruction through Bayesian statistics, and tissue analysis of tracked animals to ground-truth biogeochemical profiles. However, most researchers still focus on only a few tracers. Moreover, insufficient knowledge of the biogeochemical integration in tissues, fractionation and routing processes, and geographic and temporal variability in baseline levels continue to hamper the resolution and potential of these markers in studying the spatial and feeding ecology of top predators.

The Roles of Large Top Predators in Coastal Ecosystems: New Insights from Long Term Ecological Research

During recent human history, human activities such as overhunting and habitat destruction have severely impacted many large top predator populations around the world. Studies from a variety of ecosystems show that loss or diminishment of top predator populations can have serious consequences for population and community dynamics and ecosystem stability. However, there are relatively few studies of the roles of large top predators in coastal ecosystems, so that we do not yet completely understand what could happen to coastal areas if large top predators are extirpated or significantly reduced in number. This lack of knowledge is surprising given that coastal areas around the globe are highly valued and densely populated by humans, and thus coastal large top predator populations frequently come into conflict with coastal human populations. This paper reviews what is known about the ecological roles of large top predators in coastal systems and presents a synthesis of recent work from three coastal eastern US Long Term Ecological Research (LTER) sites where long-term studies reveal what appear to be common themes relating to the roles of large top predators in coastal systems. We discuss three specific themes: (1) large top predators acting as mobile links between disparate habitats, (2) large top predators potentially affecting nutrient and biogeochemical dynamics through localized behaviors, and (3) individual specialization of large top predator behaviors. We also discuss how research within the LTER network has led to enhanced understanding of the ecological roles of coastal large top predators. Highlighting this work is intended to encourage further investigation of the roles of large top predators across diverse coastal aquatic habitats and to better inform researchers and ecosystem managers about the importance of large top predators for coastal ecosystem health and stability.

Small but powerful: Top predator local extinction affects ecosystem structure and function in an intermittent stream.

Top predator loss is a major global problem, with a current trend in biodiversity loss towards high trophic levels that modifies most ecosystems worldwide. Most research in this area is focused on large-bodied predators, despite the high extinction risk of small-bodied freshwater fish that often act as apex consumers. Consequently, it remains unknown if intermittent streams are affected by the consequences of top-predators' extirpations. The aim of our research was to determine how this global problem affects intermittent streams and, in particular, if the loss of a small-bodied top predator (1) leads to a 'mesopredator release', affects primary consumers and changes whole community structures, and (2) triggers a cascade effect modifying the ecosystem function. To address these questions, we studied the topdown effects of a small endangered fish species, Barbus meridionalis (the Mediterranean barbel), conducting an enclosure/exclosure mesocosm experiment in an intermittent stream where B. meridionalis became locally extinct following a wildfire.We found that top predator absence led to 'mesopredator release', and also to 'prey release' despite intraguild predation, which contrasts with traditional food web theory. In addition, B. meridionalis extirpation changed whole macroinvertebrate community composition and increased total macroinvertebrate density. Regarding ecosystem function, periphyton primary production decreased in apex consumer absence. In this study, the apex consumer was functionally irreplaceable; its local extinction led to the loss of an important functional role that resulted in major changes to the ecosystem's structure and function. This study evidences that intermittent streams can be affected by the consequences of apex consumers' extinctions, and that the loss of small-bodied top predators can lead to large ecosystem changes. We recommend the reintroduction of small-bodied apex consumers to systems where they have been extirpated, to restore ecosystem structure and function.

Trophic cascades in boreal landscapes: top predator protection on tree seedlings and forest grouse

Changes in the abundance of top predators have brought about notable, cascading effects in ecosystems around the world. In this thesis, I examined several potential trophic cascades in boreal ecosystems, and their separate interspecific interactions. The main aim of the thesis was to investigate whether predators in the boreal forests have direct or indirect cascading effects on the lower trophic levels. First, I compared the browsing effects of different mammalian herbivores by excluding varying combinations of voles, hares and cervids from accessing the seedlings of silver birch (Betula pendula), Scots pine (Pinus sylvestris) and Norway spruce (Picea abies). Additionally, I studied the effect of simulated predation risk on vole browsing by using auditory cues of owls. Moving upwards on the trophic levels, I examined the intraguild interactions between the golden eagle (Aquila chrysaetos), and its mesopredator prey, the red fox (Vulpes vulpes) and the pine marten (Martes martes). To look at an entire potential trophic cascade, I further studied the combined impacts of eagles and mesopredators on the black grouse (Tetrao tetrix) and the hazel grouse (Tetrastes bonasia), predicting that the shared forest grouse prey would benefit from eagle presence. From the tree species studied, birch appears to be the most palatable one for the mammalian herbivores. I observed growth reductions in the presences of cervids and low survival associated with hares and voles, which suggests that they all weaken regeneration in birch stands. Furthermore, the simulated owl predation risk appeared to reduce vole browsing on birches in late summer, although the preferred grass forage is then old and less palatable. Browsing by voles and hares had a negative effect on the condition and survival of Scots pine, but in contrast, the impact of mammalian herbivores on spruce was found to be small, at least when more preferred food is available. I observed that the presence of golden eagles had a negative effect on the abundance of adult black grouse but a positive, protective effect on the proportion of juveniles in both black grouse and hazel grouse. Yet, this positive effect was not dependent on the abundance foxes or martens, nor did eagles seem to effectively decrease the abundance of these mesopredators. Conversely, the protection effect on grouse could arise from fear effects and also be mediated by other mesopredators. The results of this thesis provide important new information about trophic interactions in the boreal food webs. They highlight how different groups of mammalian herbivores vary in their effects on the growth and condition of different tree seedlings. Lowered cervid abundances could improve birch regeneration, which indirectly supports the idea that the key predators of cervids could cause cascading effects also in Fennoscandian forests. Owls seem to reduce vole browsing through an intimidation effect, which is a novel result of the cascading effects of owl vocalisation and could even have applications for protecting birch seedlings. In the third cascade examined in this thesis, I found the golden eagle to have a protective effect on the reproducing forest grouse, but it remains unclear through which smaller predators this effect is mediated. Overall, the results of this thesis further support the idea that there are cascading effects in the forests of Northern Europe, and that they are triggered by both direct and non‐lethal effects of predation.

Mollusca Non-Grata: The influence of top-down control and residence time on the abundance, distribution, and behavior of non-native marine snails in Washington State

Thesis (Ph.D.)--University of Washington, 2016-06

Top down control in a relatively pristine seagrass ecosystem

The loss of large-bodied herbivores and/or top predators has been associated with large-scale changes in terrestrial, freshwater, and marine ecosystems around the world. Understanding the consequences of these declines has been hampered by a lack of studies in relatively pristine systems. To fill this gap, I investigated the dynamics of the relatively pristine seagrass ecosystem of Shark Bay, Australia. I began by examining the seagrass species distributions, stoichiometry, and patterns of nutrient limitation across the whole of Shark Bay. Large areas were N-limited, P-limited, or limited by factors other than nutrients. Phosphorus-limitation was centered in areas of restricted water exchange with the ocean. Nutrient content of seagrasses varied seasonally, but the strength of seasonal responses were species-specific. Using a cafeteria-style experiment, I found that fast-growing seagrass species, which had higher nutrient content experienced higher rates of herbivory than slow-growing species that are dominant in the bay but have low nutrient content. Although removal rates correlated well with nutrient content at a broad scale, within fast-growing species removal rates were not closely tied to N or P content. Using a combination of stable isotope analysis and animal borne video, I found that green turtles (Chelonia mydas)—one of the most abundant large-bodied herbivores in Shark Bay—appear to assimilate little energy from seagrasses at the population level. There was, however, evidence of individual specialization in turtle diets with some individuals foraging largely on seagrasses and others feeding primarily on macroalgae and gelatinous macroplankton. Finally, I used exclusion cages, to examine whether predation-sensitive habitat shifts by megagrazers (green turtles, dugongs) transmitted a behavior-mediated trophic cascade (BMTC) between sharks and seagrasses. In general, data were consistent with predictions of a behavior-mediated trophic cascade. Megaherbivore impacts on seagrasses were large only in the microhabitat where megaherbivores congregate to reduce predation risk. My study highlights the importance of large herbivores in structuring seagrass communities and, more generally, suggests that roving top predators likely are important in structuring communities—and possibly ecosystems—through non-consumptive pathways.

Contrasting patterns of individual specialization and trophic coupling in two marine apex predators

1. Apex predators are often assumed to be dietary generalists and, by feeding on prey from multiple basal nutrient sources, serve to couple discrete food webs. But there is increasing evidence that individual level dietary specialization may be common in many species, and this has not been investigated for many marine apex predators. 2. Because of their position at or near the top of many marine food webs, and the possibility that they can affect populations of their prey and induce trophic cascades, it is important to understand patterns of dietary specialization in shark populations. 3. Stable isotope values from body tissues with different turnover rates were used to quantify patterns of individual specialization in two species of ‘generalist’ sharks (bull sharks, Carcharhinus leucas, and tiger sharks, Galeocerdo cuvier). 4. Despite wide population-level isotopic niche breadths in both species, isotopic values of individual tiger sharks varied across tissues with different turnover rates. The population niche breadth was explained mostly by variation within individuals suggesting tiger sharks are true generalists. In contrast, isotope values of individual bull sharks were stable through time and their wide population level niche breadth was explained by variation among specialist individuals. 5. Relative resource abundance and spatial variation in food-predation risk tradeoffs may explain the differences in patterns of specialization between shark species. 6. The differences in individual dietary specialization between tiger sharks and bull sharks results in different functional roles in coupling or compartmentalizing distinct food webs. 7. Individual specialization may be an important feature of trophic dynamics of highly mobile marine top predators and should be explicitly considered in studies of marine food webs and the ecological role of top predators.

Spatial Relationships among Hydroacoustic, Hydrographic and Top Predator Patterns: Cetacean Distributions in the Mid-Atlantic Bight

Effective conservation and management of top predators requires a comprehensive understanding of their distributions and of the underlying biological and physical processes that affect these distributions. The Mid-Atlantic Bight shelf break system is a dynamic and productive region where at least 32 species of cetaceans have been recorded through various systematic and opportunistic marine mammal surveys from the 1970s through 2012. My dissertation characterizes the spatial distribution and habitat of cetaceans in the Mid-Atlantic Bight shelf break system by utilizing marine mammal line-transect survey data, synoptic multi-frequency active acoustic data, and fine-scale hydrographic data collected during the 2011 summer Atlantic Marine Assessment Program for Protected Species (AMAPPS) survey. Although studies describing cetacean habitat and distributions have been previously conducted in the Mid-Atlantic Bight, my research specifically focuses on the shelf break region to elucidate both the physical and biological processes that influence cetacean distribution patterns within this cetacean hotspot.

In Chapter One I review biologically important areas for cetaceans in the Atlantic waters of the United States. I describe the study area, the shelf break region of the Mid-Atlantic Bight, in terms of the general oceanography, productivity and biodiversity. According to recent habitat-based cetacean density models, the shelf break region is an area of high cetacean abundance and density, yet little research is directed at understanding the mechanisms that establish this region as a cetacean hotspot.

In Chapter Two I present the basic physical principles of sound in water and describe the methodology used to categorize opportunistically collected multi-frequency active acoustic data using frequency responses techniques. Frequency response classification methods are usually employed in conjunction with net-tow data, but the logistics of the 2011 AMAPPS survey did not allow for appropriate net-tow data to be collected. Biologically meaningful information can be extracted from acoustic scattering regions by comparing the frequency response curves of acoustic regions to theoretical curves of known scattering models. Using the five frequencies on the EK60 system (18, 38, 70, 120, and 200 kHz), three categories of scatterers were defined: fish-like (with swim bladder), nekton-like (e.g., euphausiids), and plankton-like (e.g., copepods). I also employed a multi-frequency acoustic categorization method using three frequencies (18, 38, and 120 kHz) that has been used in the Gulf of Maine and Georges Bank which is based the presence or absence of volume backscatter above a threshold. This method is more objective than the comparison of frequency response curves because it uses an established backscatter value for the threshold. By removing all data below the threshold, only strong scattering information is retained.

In Chapter Three I analyze the distribution of the categorized acoustic regions of interest during the daytime cross shelf transects. Over all transects, plankton-like acoustic regions of interest were detected most frequently, followed by fish-like acoustic regions and then nekton-like acoustic regions. Plankton-like detections were the only significantly different acoustic detections per kilometer, although nekton-like detections were only slightly not significant. Using the threshold categorization method by Jech and Michaels (2006) provides a more conservative and discrete detection of acoustic scatterers and allows me to retrieve backscatter values along transects in areas that have been categorized. This provides continuous data values that can be integrated at discrete spatial increments for wavelet analysis. Wavelet analysis indicates significant spatial scales of interest for fish-like and nekton-like acoustic backscatter range from one to four kilometers and vary among transects.

In Chapter Four I analyze the fine scale distribution of cetaceans in the shelf break system of the Mid-Atlantic Bight using corrected sightings per trackline region, classification trees, multidimensional scaling, and random forest analysis. I describe habitat for common dolphins, Risso’s dolphins and sperm whales. From the distribution of cetacean sightings, patterns of habitat start to emerge: within the shelf break region of the Mid-Atlantic Bight, common dolphins were sighted more prevalently over the shelf while sperm whales were more frequently found in the deep waters offshore and Risso’s dolphins were most prevalent at the shelf break. Multidimensional scaling presents clear environmental separation among common dolphins and Risso’s dolphins and sperm whales. The sperm whale random forest habitat model had the lowest misclassification error (0.30) and the Risso’s dolphin random forest habitat model had the greatest misclassification error (0.37). Shallow water depth (less than 148 meters) was the primary variable selected in the classification model for common dolphin habitat. Distance to surface density fronts and surface temperature fronts were the primary variables selected in the classification models to describe Risso’s dolphin habitat and sperm whale habitat respectively. When mapped back into geographic space, these three cetacean species occupy different fine-scale habitats within the dynamic Mid-Atlantic Bight shelf break system.

In Chapter Five I present a summary of the previous chapters and present potential analytical steps to address ecological questions pertaining the dynamic shelf break region. Taken together, the results of my dissertation demonstrate the use of opportunistically collected data in ecosystem studies; emphasize the need to incorporate middle trophic level data and oceanographic features into cetacean habitat models; and emphasize the importance of developing more mechanistic understanding of dynamic ecosystems.

Top Down Control in a Relatively Pristine Seagrass Ecosystem

The loss of large-bodied herbivores and/or top predators has been associated with large-scale changes in terrestrial, freshwater, and marine ecosystems around the world. Understanding the consequences of these declines has been hampered by a lack of studies in relatively pristine systems. To fill this gap, I investigated the dynamics of the relatively pristine seagrass ecosystem of Shark Bay, Australia. I began by examining the seagrass species distributions, stoichiometry, and patterns of nutrient limitation across the whole of Shark Bay. Large areas were N-limited, P-limited, or limited by factors other than nutrients. Phosphorus-limitation was centered in areas of restricted water exchange with the ocean. Nutrient content of seagrasses varied seasonally, but the strength of seasonal responses were species-specific. Using a cafeteria-style experiment, I found that fast-growing seagrass species, which had higher nutrient content experienced higher rates of herbivory than slow-growing species that are dominant in the bay but have low nutrient content. Although removal rates correlated well with nutrient content at a broad scale, within fast-growing species removal rates were not closely tied to N or P content. Using a combination of stable isotope analysis and animal borne video, I found that green turtles (Chelonia mydas) – one of the most abundant large-bodied herbivores in Shark Bay – appear to assimilate little energy from seagrasses at the population level. There was, however, evidence of individual specialization in turtle diets with some individuals foraging largely on seagrasses and others feeding primarily on macroalgae and gelatinous macroplankton. Finally, I used exclusion cages, to examine whether predation-sensitive habitat shifts by megagrazers (green turtles, dugongs) transmitted a behavior-mediated trophic cascade (BMTC) between sharks and seagrasses. In general, data were consistent with predictions of a behavior-mediated trophic cascade. Megaherbivore impacts on seagrasses were large only in the microhabitat where megaherbivores congregate to reduce predation risk. My study highlights the importance of large herbivores in structuring seagrass communities and, more generally, suggests that roving top predators likely are important in structuring communities - and possibly ecosystems - through non-consumptive pathways.

Modelização de distribuição de cetáceos: comparando conjuntos de dados e criando mapas realistas de distribuição

Jornadas "Ciência nos Açores - que futuro?", Biblioteca Pública e Arquivo Regional de Ponta Delgada, Largo do Colégio, Ponta Delgada, 7-8 de junho.

Rede ecológica para a deteção e estudo de aves marinhas arrojadas REDEMAR-AVES (costa norte portuguesa)

Dissertação de mestrado em Ecologia

Relación entre anchoveta y otros componentes del ecosistema

Se hipotetiza que la mortalidad por pesca ocasiona efectos ecosistémicos, no sólo sobre la anchoveta, sino también sobre otros componentes del ecosistema, tales como los depredadores tope. El objetivo es realizar simulaciones con un modelo ecotrófico multiespecífico cubriendo el decenio de los años 2000, cambiando la mortalidad por pesca y analizando las variaciones en la biomasa de anchoveta, aves guaneras y lobos marinos. Se usó el Índice de Oscilación Peruano (IOP) para incluir una mediación que afecte la vulnerabilidad de las presas de la anchoveta. Se comparó el ajuste de los datos observados, usando dos tipos de dieta para anchoveta (fitoplanctófaga y zooplanctófaga). Se realizaron proyecciones de la biomasa, cambiando la mortalidad por pesca de 0,0 a 1,4 año-1. Las simulaciones con la dieta zooplanctófaga, que eleva el nivel trófico de la anchoveta de 2,35 a 3,36, produjo un menor ajuste entre los datos observados y simulados. La relación inversa entre la mortalidad por pesca y la biomasa desovante de la anchoveta, mostró que mortalidades por pesca entre 0,8 y 1,4 año-1 estarían asociadas a una biomasa desovante mínima de anchoveta, tomando en consideración sus relaciones multiespecíficas. También se encontró una relación inversa entre la mortalidad por pesca y las poblaciones de aves guaneras y lobos marinos.

Trophic dynamics in the northern Humboldt Current system: insights from stable isotopes and stomach content analyses

The northern Humboldt Current system (NHCS) off Peru is one of the most productive world marine regions. It represents less than 0.1% of the world ocean surface but presently sustains about 10% of the world fish catch, with the Peruvian anchovy or anchoveta Engraulis ringens as emblematic fish resource. Compared with other eastern boundary upwelling systems, the higher fish productivity of the NHCS cannot be explained by a corresponding higher primary productivity. On another hand, the NHCS is the region where El Niño, and climate variability in general, is most notable. Also, surface oxygenated waters overlie an intense and extremely shallow Oxygen Minimum Zone (OMZ). In this context, the main objective of this study is to better understand the trophic flows in the NHCS using both stomach content and stable isotope analyses. The study focuses on a variety of organisms from low trophic levels such as zooplankton to top predators (seabirds and fur seals). The approach combines both long-term and specific studies on emblematic species such as anchoveta, and sardine Sardinops sagax and a more inclusive analysis considering the 'global' food web in the recent years (2008 – 2012) using stable isotope analysis. Revisiting anchovy and sardine we show that whereas phytoplankton largely dominated anchoveta and sardine diets in terms of numerical abundance, the carbon content of prey items indicated that zooplankton was by far the most important dietary component. Indeed for anchovy euphausiids contributed 67.5% of dietary carbon, followed by copepods (26.3%). Selecting the largest prey, the euphausiids, provide an energetic advantage for anchoveta in its ecosystem where oxygen depletion imposes strong metabolic constrain to pelagic fish. Sardine feed on smaller zooplankton than do anchoveta, with sardine diet consisting of smaller copepods and fewer euphausiids than anchoveta diet. Hence, trophic competition between sardine and anchovy in the northern Humboldt Current system is minimized by their partitioning of the zooplankton food resource based on prey size, as has been reported in other systems. These results suggest an ecological role for pelagic fish that challenges previous understanding of their position in the foodweb (zooplanktophagous instead of phytophagous), the functioning and the trophic models of the NHCS. Finally to obtain a more comprehensive vision of the relative trophic position of NHCS main components we used stable isotope analyses. For that purpose we analyzed the δ13C and δ15N stable isotope values of thirteen taxonomic categories collected off Peru from 2008 - 2011, i.e., zooplankton, fish, squids and air-breathing top predators. The δ15N isotope signature was strongly impacted by the species, the body length and the latitude. Along the Peruvian coast, the OMZ get more intense and shallow south of ~7.5ºS impacting the baseline nitrogen stable isotopes. Employing a linear mixed-effects modelling approach taking into account the latitudinal and body length effects, we provide a new vision of the relative trophic position of key ecosystem components. Also we confirm stomach content-based results on anchoveta Engraulis ringens and highlight the potential remarkable importance of an often neglected ecosystem component, the squat lobster Pleuroncodes monodon. Indeed, our results support the hypothesis according to which this species forage to some extent on fish eggs and larvae and can thus predate on the first life stages of exploited species. However, the δ13C values of these two species suggest that anchoveta and squat lobster do not exactly share the same habitat. This would potentially reduce some direct competition and/or predation.

Adéquation et choix de modele de marche aléatoire des trajectoires des navires de peches

L'étude du mouvement des organismes est essentiel pour la compréhension du fonctionnement des écosystèmes. Dans le cas des écosystèmes marins exploités, cela amène à s'intéresser aux stratégies spatiales des pêcheurs. L'une des approches les plus utilisées pour la modélisation du mouvement des prédateurs supé- rieurs est la marche aléatoire de Lévy. Une marche aléatoire est un modèle mathématique composé par des déplacements aléatoires. Dans le cas de Lévy, les longueurs des déplacements suivent une loi stable de Lévy. Dans ce cas également, les longueurs, lorsqu'elles tendent vers l'in ni (in praxy lorsqu'elles sont grandes, grandes par rapport à la médiane ou au troisième quartile par exemple), suivent une loi puissance caractéristique du type de marche aléatoire de Lévy (Cauchy, Brownien ou strictement Lévy). Dans la pratique, outre que cette propriété est utilisée de façon réciproque sans fondement théorique, les queues de distribution, notion par ailleurs imprécise, sont modélisée par des lois puissances sans que soient discutées la sensibilité des résultats à la dé nition de la queue de distribution, et la pertinence des tests d'ajustement et des critères de choix de modèle. Dans ce travail portant sur les déplacements observés de trois bateaux de pêche à l'anchois du Pérou, plusieurs modèles de queues de distribution (log-normal, exponentiel, exponentiel tronqué, puissance et puissance tronqué) ont été comparés ainsi que deux dé nitions possible de queues de distribution (de la médiane à l'in ni ou du troisième quartile à l'in ni). Au plan des critères et tests statistiques utilisés, les lois tronquées (exponentielle et puissance) sont apparues les meilleures. Elles intègrent en outre le fait que, dans la pratique, les bateaux ne dépassent pas une certaine limite de longueur de déplacement. Le choix de modèle est apparu sensible au choix du début de la queue de distribution : pour un même bateau, le choix d'un modèle tronqué ou l'autre dépend de l'intervalle des valeurs de la variable sur lequel le modèle est ajusté. Pour nir, nous discutons les implications en écologie des résultats de ce travail.