Postglacial climate changes and rise of three ecotypes of harbour porpoises, Phocoena phocoena, in western Palearctic waters.

Despite no obvious barriers to gene flow in the marine realm, environmental variation and ecological specializations can lead to genetic differentiation in highly mobile predators. Here, we investigated the genetic structure of the harbour porpoise over the entire species distribution range in western Palearctic waters. Combined analyses of 10 microsatellite loci and a 5085 base-pair portion of the mitochondrial genome revealed the existence of three ecotypes, equally divergent at the mitochondrial genome, distributed in the Black Sea (BS), the European continental shelf waters, and a previously overlooked ecotype in the upwelling zones of Iberia and Mauritania. Historical demographic inferences using approximate Bayesian computation (ABC) suggest that these ecotypes diverged during the last glacial maximum (c. 23-19 kilo-years ago, kyrbp). ABC supports the hypothesis that the BS and upwelling ecotypes share a more recent common ancestor (c. 14 kyrbp) than either does with the European continental shelf ecotype (c. 28 kyrbp), suggesting they probably descended from the extinct populations that once inhabited the Mediterranean during the glacial and post-glacial period. We showed that the two Atlantic ecotypes established a narrow admixture zone in the Bay of Biscay during the last millennium, with highly asymmetric gene flow. This study highlights the impacts that climate change may have on the distribution and speciation process in pelagic predators and shows that allopatric divergence can occur in these highly mobile species and be a source of genetic diversity.

Let's go beyond taxonomy in diet description: testing a trait-based approach to prey-predator relationships.

Understanding 'Why a prey is a prey for a given predator?' can be facilitated through trait-based approaches that identify linkages between prey and predator morphological and ecological characteristics and highlight key functions involved in prey selection. Enhanced understanding of the functional relationships between predators and their prey is now essential to go beyond the traditional taxonomic framework of dietary studies and to improve our knowledge of ecosystem functioning for wildlife conservation and management. We test the relevance of a three-matrix approach in foraging ecology among a marine mammal community in the northeast Atlantic to identify the key functional traits shaping prey selection processes regardless of the taxonomy of both the predators and prey. Our study reveals that prey found in the diet of marine mammals possess functional traits which are directly and significantly linked to predator characteristics, allowing the establishment of a functional typology of marine mammal-prey relationships. We found prey selection of marine mammals was primarily shaped by physiological and morphological traits of both predators and prey, confirming that energetic costs of foraging strategies and muscular performance are major drivers of prey selection in marine mammals. We demonstrate that trait-based approaches can provide a new definition of the resource needs of predators. This framework can be used to anticipate bottom-up effects on marine predator population dynamics and to identify predators which are sensitive to the loss of key prey functional traits when prey availability is reduced.

Deep ocean drivers better explain habitat preferences of sperm whales Physeter macrocephalus than beaked whales in the Bay of Biscay.

Species Distribution Models are commonly used with surface dynamic environmental variables as proxies for prey distribution to characterise marine top predator habitats. For oceanic species that spend lot of time at depth, surface variables might not be relevant to predict deep-dwelling prey distributions. We hypothesised that descriptors of deep-water layers would better predict the deep-diving cetacean distributions than surface variables. We combined static variables and dynamic variables integrated over different depth classes of the water column into Generalised Additive Models to predict the distribution of sperm whales Physeter macrocephalus and beaked whales Ziphiidae in the Bay of Biscay, eastern North Atlantic. We identified which variables best predicted their distribution. Although the highest densities of both taxa were predicted near the continental slope and canyons, the most important variables for beaked whales appeared to be static variables and surface to subsurface dynamic variables, while for sperm whales only surface and deep-water variables were selected. This could suggest differences in foraging strategies and in the prey targeted between the two taxa. Increasing the use of variables describing the deep-water layers would provide a better understanding of the oceanic species distribution and better assist in the planning of human activities in these habitats.

Form, function, and divergence of a generic fin shape in small cetaceans.

Tail flukes as well as the dorsal fin are the apomorphic traits of cetaceans which appeared during the evolutionary process of adaptation to the aquatic life. Both appendages present a wing-like shape associated with lift generation and low drag. We hypothesized that the evolution of fins as lifting structures led to a generic wing design, where the dimensionless parameters of the fin cross-sections are invariant with respect to the body length and taxonomy of small cetaceans (Hypothesis I). We also hypothesized that constraints on variability of a generic fin shape are associated with the primary function of the fin as a fixed or flapping hydrofoil (Hypothesis II). To verify these hypotheses, we examined how the variation in the fin's morphological traits is linked to the primary function, species and body length. Hydrodynamic characteristics of the fin cross-sections were examined with the CFD software and compared with similar engineered airfoils. Generic wing design of both fins was found in a wing-like planform and a streamlined cross-sectional geometry optimized for lift generation. Divergence in a generic fin shape both on the planform and cross-sectional level was found to be related with the fin specialization in fixed or flapping hydrofoil function. Cross-sections of the dorsal fin were found to be optimized for the narrow range of small angles of attack. Cross-sections of tail flukes were found to be more stable for higher angles of attack and had gradual stall characteristics. The obtained results provide an insight into the divergent evolutionary pathways of a generic wing-like shape of the fins of cetaceans under specific demands of thrust production, swimming stability and turning control.

Towards a better characterisation of deep-diving whales' distributions by using prey distribution model outputs?

In habitat modelling, environmental variables are assumed to be proxies of lower trophic levels distribution and by extension, of marine top predator distributions. More proximal variables, such as potential prey fields, could refine relationships between top predator distributions and their environment. In situ data on prey distributions are not available over large spatial scales but, a numerical model, the Spatial Ecosystem And POpulation DYnamics Model (SEAPODYM), provides simulations of the biomass and production of zooplankton and six functional groups of micronekton at the global scale. Here, we explored whether generalised additive models fitted to simulated prey distribution data better predicted deep-diver densities (here beaked whales Ziphiidae and sperm whales Physeter macrocephalus) than models fitted to environmental variables. We assessed whether the combination of environmental and prey distribution data would further improve model fit by comparing their explanatory power. For both taxa, results were suggestive of a preference for habitats associated with topographic features and thermal fronts but also for habitats with an extended euphotic zone and with large prey of the lower mesopelagic layer. For beaked whales, no SEAPODYM variable was selected in the best model that combined the two types of variables, possibly because SEAPODYM does not accurately simulate the organisms on which beaked whales feed on. For sperm whales, the increase model performance was only marginal. SEAPODYM outputs were at best weakly correlated with sightings of deep-diving cetaceans, suggesting SEAPODYM may not accurately predict the prey fields of these taxa. This study was a first investigation and mostly highlighted the importance of the physiographic variables to understand mechanisms that influence the distribution of deep-diving cetaceans. A more systematic use of SEAPODYM could allow to better define the limits of its use and a development of the model that would simulate larger prey beyond 1,000 m would probably better characterise the prey of deep-diving cetaceans.

Genetic and historic evidence for climate-driven population fragmentation in a top cetacean predator: the harbour porpoises in European water.

Recent climate change has triggered profound reorganization in northeast Atlantic ecosystems, with substantial impact on the distribution of marine assemblages from plankton to fishes. However, assessing the repercussions on apex marine predators remains a challenging issue, especially for pelagic species. In this study, we use Bayesian coalescent modelling of microsatellite variation to track the population demographic history of one of the smallest temperate cetaceans, the harbour porpoise (Phocoena phocoena) in European waters. Combining genetic inferences with palaeo-oceanographic and historical records provides strong evidence that populations of harbour porpoises have responded markedly to the recent climate-driven reorganization in the eastern North Atlantic food web. This response includes the isolation of porpoises in Iberian waters from those further north only approximately 300 years ago with a predominant northward migration, contemporaneous with the warming trend underway since the 'Little Ice Age' period and with the ongoing retreat of cold-water fishes from the Bay of Biscay. The extinction or exodus of harbour porpoises from the Mediterranean Sea (leaving an isolated relict population in the Black Sea) has lacked a coherent explanation. The present results suggest that the fragmentation of harbour distribution range in the Mediterranean Sea was triggered during the warm 'Mid-Holocene Optimum' period (approx. 5000 years ago), by the end of the post-glacial nutrient-rich 'Sapropel' conditions that prevailed before that time.

Rise of oceanographic barriers in continuous populations of a cetacean: the genetic structure of harbour porpoises in Old World waters.

BACKGROUND: Understanding the role of seascape in shaping genetic and demographic population structure is highly challenging for marine pelagic species such as cetaceans for which there is generally little evidence of what could effectively restrict their dispersal. In the present work, we applied a combination of recent individual-based landscape genetic approaches to investigate the population genetic structure of a highly mobile extensive range cetacean, the harbour porpoise in the eastern North Atlantic, with regards to oceanographic characteristics that could constrain its dispersal. RESULTS: Analyses of 10 microsatellite loci for 752 individuals revealed that most of the sampled range in the eastern North Atlantic behaves as a 'continuous' population that widely extends over thousands of kilometres with significant isolation by distance (IBD). However, strong barriers to gene flow were detected in the south-eastern part of the range. These barriers coincided with profound changes in environmental characteristics and isolated, on a relatively small scale, porpoises from Iberian waters and on a larger scale porpoises from the Black Sea. CONCLUSION: The presence of these barriers to gene flow that coincide with profound changes in oceanographic features, together with the spatial variation in IBD strength, provide for the first time strong evidence that physical processes have a major impact on the demographic and genetic structure of a cetacean. This genetic pattern further suggests habitat-related fragmentation of the porpoise range that is likely to intensify with predicted surface ocean warming.

How many sightings to model rare marine species distributions.

Despite large efforts, datasets with few sightings are often available for rare species of marine megafauna that typically live at low densities. This paucity of data makes modelling the habitat of these taxa particularly challenging. We tested the predictive performance of different types of species distribution models fitted to decreasing numbers of sightings. Generalised additive models (GAMs) with three different residual distributions and the presence only model MaxEnt were tested on two megafauna case studies differing in both the number of sightings and ecological niches. From a dolphin (277 sightings) and an auk (1,455 sightings) datasets, we simulated rarity with a sighting thinning protocol by random sampling (without replacement) of a decreasing fraction of sightings. Better prediction of the distribution of a rarely sighted species occupying a narrow habitat (auk dataset) was expected compared to the distribution of a rarely sighted species occupying a broad habitat (dolphin dataset). We used the original datasets to set up a baseline model and fitted additional models on fewer sightings but keeping effort constant. Model predictive performance was assessed with mean squared error and area under the curve. Predictions provided by the models fitted to the thinned-out datasets were better than a homogeneous spatial distribution down to a threshold of approximately 30 sightings for a GAM with a Tweedie distribution and approximately 130 sightings for the other models. Thinning the sighting data for the taxon with narrower habitats seemed to be less detrimental to model predictive performance than for the broader habitat taxon. To generate reliable habitat modelling predictions for rarely sighted marine predators, our results suggest (1) using GAMs with a Tweedie distribution with presence-absence data and (2) implementing, as a conservative empirical measure, at least 50 sightings in the models.

Assessing cetacean surveys throughout the Mediterranean Sea: a gap analysis in environmental space.

Heterogeneous data collection in the marine environment has led to large gaps in our knowledge of marine species distributions. To fill these gaps, models calibrated on existing data may be used to predict species distributions in unsampled areas, given that available data are sufficiently representative. Our objective was to evaluate the feasibility of mapping cetacean densities across the entire Mediterranean Sea using models calibrated on available survey data and various environmental covariates. We aggregated 302,481 km of line transect survey effort conducted in the Mediterranean Sea within the past 20 years by many organisations. Survey coverage was highly heterogeneous geographically and seasonally: large data gaps were present in the eastern and southern Mediterranean and in non-summer months. We mapped the extent of interpolation versus extrapolation and the proportion of data nearby in environmental space when models calibrated on existing survey data were used for prediction across the entire Mediterranean Sea. Using model predictions to map cetacean densities in the eastern and southern Mediterranean, characterised by warmer, less productive waters, and more intense eddy activity, would lead to potentially unreliable extrapolations. We stress the need for systematic surveys of cetaceans in these environmentally unique Mediterranean waters, particularly in non-summer months.

Ocean sunfish as indicators for the 'rise of slime'.

Overfishing and ocean warming are drastically altering the community composition and size structure of marine ecosystems, eliminating large bodied species [1]. Against a backdrop of such environmental change, the heaviest of all bony fish, the ocean sunfish (Mola mola), seems an improbable survivor. Indeed this indolent giant is killed globally as bycatch, and is listed as 'Vulnerable'[2]. We undertook the most extensive aerial surveys of sunfish ever conducted and found surprisingly high abundances off the Atlantic and Mediterranean coasts of Western Europe. With up to 475 individuals per 100 km2, these figures are one order of magnitude higher than abundance estimates for other areas [3-5]. Using bioenergetic modelling, we estimate that each sunfish requires 71 kg day-1 of jellyfish, a biomass intake more than an order of magnitude greater than predicted for a similarly sized teleost. Scaled up to the population level, this equates to a remarkable 20,774 tonnes day-1 of predated jellyfish across our study area in summer. Sunfish abundance may be facilitated by overfishing and ocean warming, which together cause reduced predation of sunfish by sharks and elevated jellyfish biomass. Our combined survey and bioenergetic data provide the first-ever estimate of spatialized ocean sunfish daily food requirements, and stress the importance of this species as a global indicator for the 'rise of slime'. This hypothesis posits that, in an overfished world ocean exposed to global warming, gelatinous zooplankton should flourish, to the detriment of other mesotrophic species such as small pelagic fish, causing irreversible trophic cascades as well as a series of other environmental and economic issues.

Predicting cetacean habitats from their energetic needs and the distribution of their prey in two contrasted tropical regions.

To date, most habitat models of cetaceans have relied on static and oceanographic covariates, and very few have related cetaceans directly to the distribution of their prey, as a result of the limited availability of prey data. By simulating the distribution of six functional micronekton groups between the surface and ≃1,000 m deep, the SEAPODYM model provides valuable insights into prey distributions. We used SEAPODYM outputs to investigate the habitat of three cetacean guilds with increasing energy requirements: sperm and beaked whales, Globicephalinae and Delphininae. We expected High Energy Requirements cetaceans to preferentially forage in habitats of high prey biomass and/or production, where they might easily meet their high energetic needs, and Low Energy Requirements cetaceans to forage in habitats of either high or low prey biomass and/or production. Cetacean sightings were collected from dedicated aerial surveys in the South West Indian Ocean (SWIO) and French Polynesia (FP). We examined cetacean densities in relation to simulated distributions of their potential prey using Generalised Additive Models and predicted their habitats in both regions. Results supported their known diving abilities, with Delphininae mostly related to prey present in the upper layers of the water column, and Globicephalinae and sperm and beaked whales also related to prey present in deeper layers. Explained deviances ranged from 9% for sperm and beaked whales in the SWIO to 47% for Globicephalinae in FP. Delphininae and Globicephalinae appeared to select areas where high prey biomass and/or production were available at shallow depths. In contrast, sperm and beaked whales showed less clear habitat selection. Using simulated prey distributions as predictors in cetacean habitat models is crucial to understand their strategies of habitat selection in the three dimensions of the ocean.

Combining Methods to Describe Important Marine Habitats for Top Predators: Application to Identify Biological Hotspots in Tropical Waters.

In tropical waters resources are usually scarce and patchy, and predatory species generally show specific adaptations for foraging. Tropical seabirds often forage in association with sub-surface predators that create feeding opportunities by bringing prey close to the surface, and the birds often aggregate in large multispecific flocks. Here we hypothesize that frigatebirds, a tropical seabird adapted to foraging with low energetic costs, could be a good predictor of the distribution of their associated predatory species, including other seabirds (e.g. boobies, terns) and subsurface predators (e.g., dolphins, tunas). To test this hypothesis, we compared distribution patterns of marine predators in the Mozambique Channel based on a long-term dataset of both vessel- and aerial surveys, as well as tracking data of frigatebirds. By developing species distribution models (SDMs), we identified key marine areas for tropical predators in relation to contemporaneous oceanographic features to investigate multi-species spatial overlap areas and identify predator hotspots in the Mozambique Channel. SDMs reasonably matched observed patterns and both static (e.g. bathymetry) and dynamic (e.g. Chlorophyll a concentration and sea surface temperature) factors were important explaining predator distribution patterns. We found that the distribution of frigatebirds included the distributions of the associated species. The central part of the channel appeared to be the best habitat for the four groups of species considered in this study (frigatebirds, brown terns, boobies and sub-surface predators).

Understanding the distribution of marine megafauna in the English channel region: identifying key habitats for conservation within the busiest seaway on earth.

The temperate waters of the North-Eastern Atlantic have a long history of maritime resource richness and, as a result, the European Union is endeavouring to maintain regional productivity and biodiversity. At the intersection of these aims lies potential conflict, signalling the need for integrated, cross-border management approaches. This paper focuses on the marine megafauna of the region. This guild of consumers was formerly abundant, but is now depleted and protected under various national and international legislative structures. We present a meta-analysis of available megafauna datasets using presence-only distribution models to characterise suitable habitat and identify spatially-important regions within the English Channel and southern bight of the North Sea. The integration of studies from dedicated and opportunistic observer programmes in the United Kingdom and France provide a valuable perspective on the spatial and seasonal distribution of various taxonomic groups, including large pelagic fishes and sharks, marine mammals, seabirds and marine turtles. The Western English Channel emerged as a hotspot of biodiversity for megafauna, while species richness was low in the Eastern English Channel. Spatial conservation planning is complicated by the highly mobile nature of marine megafauna, however they are important components of the marine environment and understanding their distribution is a first crucial step toward their inclusion into marine ecosystem management.

The stranding anomaly as population indicator: the case of harbour porpoise Phocoena phocoena in North-Western Europe.

Ecological indicators for monitoring strategies are expected to combine three major characteristics: ecological significance, statistical credibility, and cost-effectiveness. Strategies based on stranding networks rank highly in cost-effectiveness, but their ecological significance and statistical credibility are disputed. Our present goal is to improve the value of stranding data as population indicator as part of monitoring strategies by constructing the spatial and temporal null hypothesis for strandings. The null hypothesis is defined as: small cetacean distribution and mortality are uniform in space and constant in time. We used a drift model to map stranding probabilities and predict stranding patterns of cetacean carcasses under H0 across the North Sea, the Channel and the Bay of Biscay, for the period 1990-2009. As the most common cetacean occurring in this area, we chose the harbour porpoise Phocoena phocoena for our modelling. The difference between these strandings expected under H0 and observed strandings is defined as the stranding anomaly. It constituted the stranding data series corrected for drift conditions. Seasonal decomposition of stranding anomaly suggested that drift conditions did not explain observed seasonal variations of porpoise strandings. Long-term stranding anomalies increased first in the southern North Sea, the Channel and Bay of Biscay coasts, and finally the eastern North Sea. The hypothesis of changes in porpoise distribution was consistent with local visual surveys, mostly SCANS surveys (1994 and 2005). This new indicator could be applied to cetacean populations across the world and more widely to marine megafauna.

Assessing the impact of bycatch on dolphin populations: the case of the common dolphin in the eastern North Atlantic.

Fisheries interactions have been implicated in the decline of many marine vertebrates worldwide. In the eastern North Atlantic, at least 1000 common dolphins (Delphinus delphis) are bycaught each year, particularly in pelagic pair-trawls. We have assessed the resulting impact of bycatch on this population using a demographic modeling approach. We relied on a sample of females stranded along the French Atlantic and western Channel coasts. Strandings represent an extensive source of demographic information to monitor our study population. Necropsy analysis provided an estimate of individual age and reproductive state. Then we estimated effective survivorship (including natural and human-induced mortality), age at first reproduction and pregnancy rates. Reproductive parameters were consistent with literature, but effective survivorship was unexpectedly low. Demographic parameters were then used as inputs in two models. A constant parameter matrix proposed an effective growth rate of -5.5±0.5%, corresponding to the current situation (including bycatch mortality). Subsequently, deterministic projections suggested that the population would be reduced to 20% of its current size in 30 years and would be extinct in 100 years. The demographic invariant model suggested a maximum growth rate of +4.5±0.09%, corresponding to the optimal demographic situation. Then, a risk analysis incorporating Potential Biological Removal (PBR), based on two plausible scenarii for stock structure suggested that bycatch level was unsustainable for the neritic population of the Bay of Biscay under a two-stock scenario. In depth assessment of stock structure and improved observer programs to provide scientifically robust bycatch estimates are needed. Effective conservation measures would be reducing bycatch to less than 50% of the current level in the neritic stock to reach PBR. Our approach provided indicators of the status and trajectory of the common dolphin population in the eastern North Atlantic and therefore proved to be a valuable tool for management, applicable to other dolphin populations.

Cost of living dictates what whales, dolphins and porpoises eat: the importance of prey quality on predator foraging strategies.

Understanding the mechanisms that drive prey selection is a major challenge in foraging ecology. Most studies of foraging strategies have focused on behavioural costs, and have generally failed to recognize that differences in the quality of prey may be as important to predators as the costs of acquisition. Here, we tested whether there is a relationship between the quality of diets (kJ · g(-1)) consumed by cetaceans in the North Atlantic and their metabolic costs of living as estimated by indicators of muscle performance (mitochondrial density, n = 60, and lipid content, n = 37). We found that the cost of living of 11 cetacean species is tightly coupled with the quality of prey they consume. This relationship between diet quality and cost of living appears to be independent of phylogeny and body size, and runs counter to predictions that stem from the well-known scaling relationships between mass and metabolic rates. Our finding suggests that the quality of prey rather than the sheer quantity of food is a major determinant of foraging strategies employed by predators to meet their specific energy requirements. This predator-specific dependence on food quality appears to reflect the evolution of ecological strategies at a species level, and has implications for risk assessment associated with the consequences of changing the quality and quantities of prey available to top predators in marine ecosystems.