An evaluation of surface micro- and mesoplastic pollution in pelagic ecosystems of the Western Mediterranean Sea

This study examines the distribution, abundance and characteristics of surface micro- and mesoplastic debris in the Western Mediterranean Sea. 41 samples were collected in 2011 (summer) and 2012 (summer). Results, firstly, revealed that micro- (<5mm) and mesoplastic debris were widely and uniformly distributed in this area with average concentrations of 130,000 parts/km(2) and 5700 parts/km(2), respectively. Importantly, a strong correlation between micro- and mesoplastic concentrations was identified. Secondly, a classification based on the shape and appearance of microplastics indicated the predominant presence of fragments (73 %) followed by thin films (14 %). Thirdly, the average mass ratio of microplastic to dry organic matter has been measured at 0.5, revealing a significant presence of microplastics in comparison to plankton. Finally, a correction method was applied in order to correct wind mixing effect on microplastics' vertical distribution. This data allows for a comprehensive view, for the first time, of the spatial distribution and nature of plastic debris in the Western Mediterranean Sea.

Radiogenic isotopes in the Mediterranean Sea: water mass exchange and precessional variability during the Messinian

During the late Miocene, exchange between the Mediterranean Sea and Atlantic Ocean changed dramatically, culminating in the Messinian Salinity Crisis (MSC). Understanding Mediterranean-Atlantic exchange at that time could answer the enigmatic question of how so much salt built up within the Mediterranean, while furthering the development of a framework for future studies attempting to understand how changes may have impacted global thermohaline circulation. Due to their association with specific water masses at different scales, radiogenic Sr, Pb, and Nd isotope records were generated from various archives contained within marine deposits to endeavour to understand better late Miocene Mediterranean-Atlantic exchange. The archives used include foraminiferal calcite (Sr), fish teeth and bone (Nd), dispersed authigenic ferromanganese oxyhydroxides (Nd, Pb), and a ferromanganese crust (Pb). The primary focus is on sediments preserved at one end of the Betic corridor, a gateway that once connected the Mediterranean to the Atlantic through southern Spain, although other locations are investigated. The Betic gateway terminated within several marginal sub-basins before entering the Western Mediterranean; one of these is the Sorbas Basin, a well-studied location whose sediments have been astronomically tuned at high temporal resolution, providing the necessary age control for sub-precessional resolution records. Since the climatic history of the Mediterranean is strongly controlled by precessional changes in regional climate, the aim was to produce records at high (sub-precessional) temporal resolution, to be able to observe clearly any precessional cyclicity driven by regional climate which could be superimposed over longer trends. This goal was achieved for all records except the ferromanganese crust record. The 87Sr/86Sr isotope record (Ch. 3) shows precessional frequency excursions away from the global seawater curve. As precessional frequency oscillations are unexpected for this setting, a numerical box model was used to determine the mechanisms causing the excursions. To enable parameterisation of model variables, regional Sr characteristics, data from general circulation model HadCM3L, and new benthic foraminiferal assemblage data are employed. The model results imply that the Sorbas Basin likely had a positive hydrologic budget in the late Miocene, very different to that of today. Moreover, the model indicates that the mechanism controlling the Sr isotope ratio of Sorbas Basin seawater was not restriction, but a lack of density-driven exchange with the Mediterranean. Beyond improving our understanding of how marginal Mediterranean sub-basins may evolve different isotope signatures, these results have implications for astronomical tuning and stratigraphy in the region, findings which are crucial considering the geological and climatic history of the late Miocene Mediterranean is based entirely on marginal deposits. An improved estimate for the Nd isotope signature of late Miocene Mediterranean Outflow (MO) was determined by comparing Nd isotope signatures preserved in the deeper Alborán Sea at ODP Site 978 with literature data as well as the signature preserved in the Sorbas Basin (Ch. 4; -9.34 to -9.92 ± 0.37 εNd(t)). It was also inferred that it is unlikely that Nd isotopes can be used reliably to track changes in circulation within the shallow settings characteristic of the Mediterranean-Atlantic connections; this is significant in light of a recent publication documenting corridor closure using Nd isotopes. Both conclusions will prove useful for future studies attempting to understand changes in Mediterranean-Atlantic exchange. Excursions to high values, with precessional frequency, are also observed in the radiogenic Pb isotope record for the Sorbas Basin (Ch. 5). Widening the scope to include locations further away from the gateways, records were produced for late Miocene sections on Sicily and Northern Italy, and similar precessional frequency cyclicity was observed in the Pb isotope records for these sites as well. Comparing these records to proxies for Saharan dust and available whole rock data indicates that, while further analysis is necessary to draw strong conclusions, enhanced dust production during insolation minima may be driving the observed signal. These records also have implications for astronomical tuning; peaks in Pb isotope records driven by Saharan dust may be easier to connect directly to the insolation cycle, providing improved astronomical tuning points. Finally, a Pb isotope record derived using in-situ laser ablation performed on ferromanganese crust 3514-6 from the Lion Seamount, located west of Gibraltar within the MO plume, has provided evidence that plume depth shifted during the Pliocene. The record also suggests that Pb isotopes may not be a suitable proxy for changes in late Miocene Mediterranean-Atlantic exchange, since the Pb isotope signatures of regional water masses are too similar. To develop this record, the first published instance of laser ablation derived 230Thexcess measurements are combined with 10Be dating.

Sedimentary archives of climate and sea-level changes during the Holocene in the Rhone prodelta (NW Mediterranean Sea)

A 7.38 m-long sediment core was collected from the eastern part of the Rhone prodelta (NW Mediterranean) at 67 m water depth. A multi-proxy study (sedimentary facies, benthic foraminifera and ostracods, clay mineralogy, and major elements from XRF) provides a multi-decadal to century-scale record of climate and sea-level changes during the Holocene. The early Holocene is marked by alternative silt and clay layers interpreted as distal tempestites deposited in a context of rising sea level. This interval contains shallow infra-littoral benthic meiofauna (e.g. Pontocythere elongata, Elphidium spp., Quinqueloculina lata) and formed between ca. 20 and 50 m water depth. The middle Holocene (ca. 8.3 to 4.5 ka cal. BP), is characterized, at the core site, by a period of sediment starvation (accumulation rate of ca. 0.01 cm yr−1) resulting from the maximum landward shift of the shoreline and the Rhone outlet(s). From a sequence stratigraphic point of view, this condensed interval, about 35 cm-thick, is a Maximum Flooding Surface that can be identified on seismic profiles as the transition between delta retrogradation and delta progradation. It is marked by very distinct changes in all proxy records. Following the stabilization of the global sea level, the late Holocene is marked by the establishment of prodeltaic conditions at the core site, as shown by the lithofacies and by the presence of benthic meiofauna typical of the modern Rhone prodelta (e.g. Valvulineria bradyana, Cassidulina carinata, Bulimina marginata). Several periods of increased fluvial discharge are also emphasized by the presence of species commonly found in brackish and shallow water environments (e.g. Leptocythere). Some of these periods correspond to the multi-decadal to centennial late Holocene humid periods recognized in Europe (i.e. the 2.8 ka event and the Little Ice Age). Two other periods of increased runoffs at ca. 1.3 and 1.1 ka cal. BP are recognized, and are likely to reflect periods of regional climate deterioration that are observed in the Rhone watershed.

Characterizing, modelling and understanding the climate variability of the deep water formation in the North-Western Mediterranean Sea

Observing, modelling and understanding the climate-scale variability of the deep water formation (DWF) in the North-Western Mediterranean Sea remains today very challenging. In this study, we first characterize the interannual variability of this phenomenon by a thorough reanalysis of observations in order to establish reference time series. These quantitative indicators include 31 observed years for the yearly maximum mixed layer depth over the period 1980–2013 and a detailed multi-indicator description of the period 2007–2013. Then a 1980–2013 hindcast simulation is performed with a fully-coupled regional climate system model including the high-resolution representation of the regional atmosphere, ocean, land-surface and rivers. The simulation reproduces quantitatively well the mean behaviour and the large interannual variability of the DWF phenomenon. The model shows convection deeper than 1000 m in 2/3 of the modelled winters, a mean DWF rate equal to 0.35 Sv with maximum values of 1.7 (resp. 1.6) Sv in 2013 (resp. 2005). Using the model results, the winter-integrated buoyancy loss over the Gulf of Lions is identified as the primary driving factor of the DWF interannual variability and explains, alone, around 50 % of its variance. It is itself explained by the occurrence of few stormy days during winter. At daily scale, the Atlantic ridge weather regime is identified as favourable to strong buoyancy losses and therefore DWF, whereas the positive phase of the North Atlantic oscillation is unfavourable. The driving role of the vertical stratification in autumn, a measure of the water column inhibition to mixing, has also been analyzed. Combining both driving factors allows to explain more than 70 % of the interannual variance of the phenomenon and in particular the occurrence of the five strongest convective years of the model (1981, 1999, 2005, 2009, 2013). The model simulates qualitatively well the trends in the deep waters (warming, saltening, increase in the dense water volume, increase in the bottom water density) despite an underestimation of the salinity and density trends. These deep trends come from a heat and salt accumulation during the 1980s and the 1990s in the surface and intermediate layers of the Gulf of Lions before being transferred stepwise towards the deep layers when very convective years occur in 1999 and later. The salinity increase in the near Atlantic Ocean surface layers seems to be the external forcing that finally leads to these deep trends. In the future, our results may allow to better understand the behaviour of the DWF phenomenon in Mediterranean Sea simulations in hindcast, forecast, reanalysis or future climate change scenario modes. The robustness of the obtained results must be however confirmed in multi-model studies.

Using SeaDataNet Management System to preserve the XBT Data-Set of the Mediterranean Sea

significant amount of Expendable Bathythermograph (XBT) data has been collected in the Mediterranean Sea since 1999 in the framework of operational oceanography activities. The management and storage of such a volume of data poses significant challenges and opportunities. The SeaDataNet project, a pan-European infrastructure for marine data diffusion, provides a convenient way to avoid dispersion of these temperature vertical profiles and to facilitate access to a wider public. The XBT data flow, along with the recent improvements in the quality check procedures and the consistence of the available historical data set are described. The main features of SeaDataNet services and the advantage of using this system for long-term data archiving are presented. Finally, focus on the Ligurian Sea is included in order to provide an example of the kind of information and final products devoted to different users can be easily derived from the SeaDataNet web portal.

Effects of fishery protection on biometry and genetic structure of two target sea cucumber species from the Mediterranean Sea

Sea cucumber fisheries are now occurring in most of the tropical areas of the world, having expanded from its origin in the central Indo-Pacific. Due to the overexploitation of these resources and the increasing demand from Asian countries, new target species from Mediterranean Sea and northeastern Atlantic Ocean are being caught. The fishery effects on biometry and genetic structure of two target species (Holothuria polii and H. tubulosa) from Turkey, were assessed. The heaviest and largest individuals of H. polii were found into the non-fishery area of Kusadasi, also showing the highest genetic diversity. Similar pattern was detected in H. tubulosa, but only the weight was significantly higher in the protected area. However, the observed differences on the fishery effects between species, could be explained considering the different percentage of catches (80% for H. polii and 20% for H. tubulosa).

Genetic differentiation among Parastichopus regalis populations in the Western Mediterranean Sea: potential effects from its fishery and current connectivity

Parastichopus regalis (Cuvier, 1817) is the most expensive seafood product on the Catalonian market (NE Spain), with prices at approximately 130 €/Kg (fresh weight). Despite its ecological and economic importance, biological and genetic information on this sea cucumber species is scarce. Here, we provide both the first insight on the genetic structure of P. regalis using sequences of cytochrome oxidase I (COI) and 16S genes and a morphological description of its population. Individual sea cucumbers were collected in six locations along the Spanish Mediterranean coast, including an area under fishery pressure (Catalonia). We found high haplotype diversity and low nucleotide diversity for both genes, with higher levels of genetic diversity observed in the COI gene. The population pairwise fixation index (FST), AMOVA and correspondence analysis (CA) based on the COI gene revealed significant genetic differentiation among some locations. However, further analysis using nuclear markers (e.g., microsatellites) is necessary to corroborate these results. Moreover, the genetic and morphological data may indicate fishery effects on the Catalonian population with a decrease in the size and weight averages and lower genetic diversity compared with locations that lack fishery pressure. For the appropriate management of this species, we suggest the following: 1) accurately assessing the stock status along the Spanish coasts; 2) studying the reproductive cycle of this target species and the establishment of a closed fishery season according to the reproductive cycle; and 3) establishing protected areas (i.e., not take zones) to conserve healthy populations and favour recruitment in the nearby areas.

West versus East Mediterranean Sea: origin and genetic differentiation of the sea cucumber Holothuria polii

We studied the genetic structure of the sea cucumber Holothuria (Roweothuria) polii (Delle Chiaje 1823) by analysing the mitochondrial DNA variation in two fragments of cytochrome oxidase I (COI) and 16S genes. Individuals were collected in seven locations along the Mediterranean Sea, which cover a wide range of the species distribution. We found high haplotype diversity for COI and moderate diversity for 16S, and low nucleotide diversity for both genes. Our results for the COI gene showed many recent and exclusive haplotypes with few mutational changes, suggesting recent or ongoing population expansion. The Western and Eastern Mediterranean populations exhibited slight but significant genetic differentiation (COI gene) with higher genetic diversity in the East. The most ancient haplotype was not present in the westernmost sampling location (SE Spain). The oldest expansion time was observed in Turkey, corresponding to mid-Pleistocene. Turkey had also the highest genetic diversity (number of total and exclusive haplotypes, polymorphisms, haplotype and nucleotide diversity). This suggests that this region could be the origin of the subsequent colonizations through the Mediterranean Sea, a hypothesis that should be assessed with nuclear markers in future research.

Phylogeography of the Atlanto-Mediterranean sea cucumber Holothuria (Holothuria) mammata: the combined effects of historical processes and current oceanographical pattern

We assessed the genetic structure of populations of the widely distributed sea cucumber Holothuria (Holothuria) mammata Grube, 1840, and investigated the effects of marine barriers to gene flow and historical processes. Several potential genetic breaks were considered, which would separate the Atlantic and Mediterranean basins, the isolated Macaronesian Islands from the other locations analysed, and the Western Mediterranean and Aegean Sea (Eastern Mediterranean). We analysed mitochondrial 16S and COI gene sequences from 177 individuals from four Atlantic locations and four Mediterranean locations. Haplotype diversity was high (H = 0.9307 for 16S and 0.9203 for COI), and the haplotypes were closely related (p = 0.0058 for 16S and 0.0071 for COI). The lowest genetic diversities were found in the Aegean Sea population. Our results showed that the COI gene was more variable and more useful for the detection of population structure than the 16S gene. The distribution of mtDNA haplotypes, the pairwise FST values and the results of exact tests and AMOVA revealed: (i) a significant genetic break between the population in the Aegean Sea and those in the other locations, as supported by both mitochondrial genes, and (ii) weak differentiation of the Canary and Azores Islands from the other populations; however, the populations from the Macaronesian Islands, Algarve and West Mediterranean could be considered to be a panmictic metapopulation. Isolation by distance was not identified in H. (H.) mammata. Historical events behind the observed findings, together with the current oceanographic patterns, were proposed and discussed as the main factors that determine the population structure and genetic signature of H. (H.) mammata

First molecular determination of herpesvirus from two mysticete species stranded in the Mediterranean Sea

BACKGROUND Herpesvirus can infect a wide range of animal species: mammals, birds, reptiles, fish, amphibians and bivalves. In marine mammals, several alpha- and gammaherpesvirus have been identified in some cetaceans and pinnipeds species. To date, however, this virus has not been detected in any member of the Balaenoptera genus. CASE PRESENTATION Herpesvirus was determined by molecular methods in tissue samples from a male fin whale juvenile (Balaenoptera physalus) and a female common minke whale calf (Balaenoptera acutorostrata) stranded on the Mediterranean coast of the Region of Valencia (Spain). Samples of skin and penile mucosa from the fin whale and samples of skin, muscle and central nervous system tissue from the common minke whale tested positive for herpesvirus based on sequences of the DNA polymerase gene. Sequences from fin whale were identical and belonged to the Alphaherpesvirinae subfamily. Only members of the Gammaherpesvirinae subfamily were amplified from the common minke whale, and sequences from the muscle and central nervous system were identical. Sequences in GenBank most closely related to these novel sequences were viruses isolated from other cetacean species, consistent with previous observations that herpesviruses show similar phylogenetic branching as their hosts. CONCLUSIONS To our knowledge, this is the first molecular determination of herpesvirus in the Balaenoptera genus. It shows that herpesvirus should be included in virological evaluation of these animals.

Unusual striped dolphin mass mortality episode related to cetacean morbillivirus in the Spanish Mediterranean sea

BACKGROUND In the last 20 years, Cetacean Morbillivirus (CeMV) has been responsible for many die-offs in marine mammals worldwide, as clearly exemplified by the two dolphin morbillivirus (DMV) epizootics of 1990-1992 and 2006-2008, which affected Mediterranean striped dolphins (Stenella coeruleoalba). Between March and April 2011, the number of strandings on the Valencian Community coast (E Spain) increased. CASE PRESENTATION Necropsy and sample collection were performed in all stranded animals, with good state of conservation. Subsequently, histopathology, immunohistochemistry, conventional reverse transcription polymerase chain reaction (RT-PCR) and Universal Probe Library (UPL) RT-PCR assays were performed to identify Morbillivirus. Gross and microscopic findings compatible with CeMV were found in the majority of analyzed animals. Immunopositivity in the brain and UPL RT-PCR positivity in seven of the nine analyzed animals in at least two tissues confirmed CeMV systemic infection. Phylogenetic analysis, based on sequencing part of the phosphoprotein gene, showed that this isolate is a closely related dolphin morbillivirus (DMV) to that responsible for the 2006-2008 epizootics. CONCLUSION The combination of gross and histopathologic findings compatible with DMV with immunopositivity and molecular detection of DMV suggests that this DMV strain could cause this die-off event.

The mesophotic zone of the Mediterranean Sea: spatial definition, biodiversity, and predictive models

The mesophotic zone is frequently defined as ranging between 30-40 and 150 m depth. However, these borders are necessarily imprecise due to variations in the penetration of light along the water column related to local factors. Moreover, density of data on mesophotic ecosystems vary along geographical distance, with temperate latitudes largely less explored than tropical situations. This is the case of the Mediterranean Sea, where information on mesophotic ecosystems is largely lower with respect to tropical situations. The lack of a clear definition of the borders of the mesophotic zone may represent a problem when information must be transferred to the policy that requires a coherent spatial definition to plan proper management and conservation measures. The present thesis aims at providing information on the spatial definition of the mesophotic zone in the Mediterranean Sea, its biodiversity and distribution of its ecosystems. The first chapter analyzes information on mesophotic ecosystems in the Mediterranean Sea to identify gaps in the literature and map the mesophotic zone in the Mediterranean Sea using light penetration estimated from satellite data. In the second chapter, different visual techniques to study mesophotic ecosystems are compared to identify the best analytical method to estimate diversity and habitat extension. In the third chapter, a set of Remotely Operated vehicles (ROV) surveys performed on mesophotic assemblages in the Mediterranean Sea are analyzed to describe their taxonomic and functional diversity and environmental factors influencing their structure. A Habitat Suitability Model is run in the fourth chapter to map the distribution of areas suitable for the presence of deep-water oyster reefs in the Adriatic-Ionian area. The fifth chapter explores the mesophotic zone in the northern Gulf of Mexico providing its spatial and vertical extension of the mesophotic zone and information on the diversity associated with mesophotic ecosystems.

Moving toward stock units identification based on spatial population structure of marine species in the Mediterranean Sea and adjacent waters through multidisciplinary and holistic approaches for the sustainability of fisheries resources

Investigating stock identity of marine species in a multidisciplinary holistic approach can reveal patterns of complex spatial population structure and signatures of potential local adaptation. The population structure of common sole (Solea solea) in the Mediterranean Sea was delineated using genomic and otolith data, including single nucleotide polymorphisms (SNPs) markers and otolith data. SNPs were correlated with environmental and spatial variables to evaluate the impact of these features on the actual genetic population structure. Integrated holistic approach was applied to combine the tracers with different spatio-temporal scales. SNPs data was also used to illustrate the population structure of European hake (Merluccius merluccius) within the Alboran Sea, extending into the neighboring Mediterranean Sea and Atlantic Ocean. The aim was to identify patterns of neutral and potential adaptive genetic variation by applying seascape genomic framework. Results from both genetic and otolith data suggested significant divergence among putative populations of common sole, confirming a clear separation between Western, Adriatic Sea and Eastern Mediterranean Sea. Evidence of fine-scale population structure in the Western Mediterranean Sea was observed at outlier loci level and in the Adriatic. Our study not only indicates that separation among Mediterranean sole population is led primarily by neutral processes, but it also suggests the presence of local adaptation influenced by environmental and spatial factors. The holistic approach by considering the spatio-temporal scales of variation confirmed that the same pattern of separation between these geographical sites is currently occurring and has occurred for many generations. Results showed the occurrence of population structure in Merluccius merluccius by detecting westward–eastward differentiation among populations and distinct subgroups at a fine geographical scale using outlier SNPs. These results enhance the knowledge of the population structure of commercially relevant species to support the application of spatial stock assessment models, including a redefinition of fishery management units.

Toward an understanding of the atmospheric forcing uncertainties in the Mediterranean Sea

The thesis has extensively investigated for the first time the statistical distributions of atmospheric surface variables and heat fluxes for the Mediterranean Sea. After retrieving a 30-year atmospheric analysis dataset, we have captured the spatial patterns of the probability distribution of the relevant atmospheric variables for ocean atmospheric forcing: wind components (U,V), wind amplitude, air temperature (T2M), dewpoint temperature (D2M) and mean sea-level pressure (MSL-P). The study reveals that a two-parameter PDF is not a good fit for T2M, D2M, MSL-P and wind components (U,V) and a three parameter skew-normal PDF is better suited. Such distribution captures properly the data asymmetric tails (skewness). After removing the large seasonal cycle, we show the quality of the fit and the geographic structure of the PDF parameters. It is found that the PDF parameters vary between different regions, in particular the shape (connected to the asymmetric tails) and the scale (connected to the spread of the distribution) parameters cluster around two or more values, probably connected to the different dynamics that produces the surface atmospheric fields in the Mediterranean basin. Moreover, using the atmospheric variables, we have computed the air-sea heat fluxes for a 20-years period and estimated the net heat budget over the Mediterranean Sea. Interestingly, the higher resolution analysis dataset provides a negative heat budget of –3 W/m2 which is within the acceptable range for the Mediterranean Sea heat budget closure. The lower resolution atmospheric reanalysis dataset(ERA5) does not satisfy the heat budget closure problem pointing out that a minimal resolution of the atmospheric forcing is crucial for the Mediterranean Sea dynamics. The PDF framework developed in this thesis will be the basis for a future ensemble forecasting system that will use the statistical distributions to create perturbations of the atmospheric ocean forcing.