Report of the Workshop on Growth-increment Chronologies in Marine Fish: climate-ecosystem interactions in the North Atlantic 2 (WKGIC2). 18-22 April 2016 Esporles, Mallorca, Spain

Over the past several decades, thousands of otoliths, bivalve shells, and scales have been collected for the purposes of age determination and remain archived in European and North American fisheries laboratories. Advances in digital imaging and computer software combined with techniques developed by tree-ring scientists provide a means by which to extract additional levels of information in these calcified structures and generate annually resolved (one value per year), multidecadal time-series of population-level growth anomalies. Chemical and isotopic properties may also be extracted to provide additional information regarding the environmental conditions these organisms experienced.Given that they are exactly placed in time, chronologies can be directly compared to instrumental climate records, chronologies from other regions or species, or time-seriesof other biological phenomena. In this way, chronologies may be used to reconstruct historical ranges of environmental variability, identify climatic drivers of growth, establish linkages within and among species, and generate ecosystem-level indicators. Following the first workshop in Hamburg, Germany, in December 2014, the second workshop on Growth increment Chronologies in Marine Fish: climate-ecosystem interactions in the North Atlantic (WKGIC2) met at the Mediterranean Institute for Advanced Studies headquarters in Esporles, Spain, on 18–22 April 2016, chaired by Bryan Black (USA) and Christoph Stransky (Germany).Thirty-six participants from fifteen different countries attended. Objectives were to i) review the applications of chronologies developed from growth-increment widths in the hard parts (otoliths, shells, scales) of marine fish and bivalve species ii) review the fundamentals of crossdating and chronology development, iii) discuss assumptions and limitations of these approaches, iv) measure otolith growth-increment widths in image analysis software, v) learn software to statistically check increment dating accuracy, vi) generate a growth increment chronology and relate it to climate indices, and vii) initiate cooperative projects or training exercises to commence after the workshop.The workshop began with an overview of tree-ring techniques of chronology development, including a hands-on exercise in cross dating. Next, we discussed the applications of fish and bivalve biochronologies and the range of issues that could be addressed. We then reviewed key assumptions and limitations, especially those associated with short-lived species for which there are numerous and extensive otolith archives in European fisheries labs. Next, participants were provided with images of European plaice otoliths from the North Sea and taught to measure increment widths in image analysis software. Upon completion of measurements, techniques of chronology development were discussed and contrasted to those that have been applied for long-lived species. Plaice growth time-series were then related to environmental variability using the KNMI Climate Explorer. Finally, potential future collaborations and funding opportunities were discussed, and there was a clear desire to meet again to compare various statistical techniques for chronology development using a range existing fish, bivalve, and tree growth-increment datasets. Overall, we hope to increase the use of these techniques, and over the long term, develop networks of biochronologies for integrative analyses of ecosystem functioning and relationships to long-term climate variability and fishing pressure.

Otolith shape variation provides a marker of stock origin for north Atlantic bluefin tuna (Thunnus thynnus)

Two stocks of bluefin tuna (Thunnus thynnus) inhabit the north Atlantic; the western and eastern stocks spawn in the Gulf of Mexico and the Mediterranean Sea respectively. Trans-Atlantic movements occur outside spawning time whereas natal homing maintains stock structure. Commercial fisheries may exploit a mixed assemblage of both stocks. The incorporation of mixing rates into stock assessment is precluded by uncertainties surrounding stock discrimination. Otolith shape descriptors were used to characterise western and eastern stocks of Atlantic bluefin tuna in the present study and to estimate stock composition in catches of unknown origin. Otolith shape varied with length and between locations and years. Within a restricted size range (200-297-cm fork length (FL)) the two stocks were distinguished with an accuracy of 83%. Bayesian stock mixture analysis indicated that samples from the east Atlantic and Mediterranean were predominantly of eastern origin. The proportion assigned to the eastern stock showed slight spatial variation; however, overlapping 95% credible intervals indicated no significant difference (200-297 cm FL: central Atlantic, 73-100%; Straits of Gibraltar, 73-100%; Morocco, 50-99%; Portugal 64-100%). Otolith shape could be used in combination with other population markers to improve the accuracy of mixing rate estimates for Atlantic bluefin tuna.

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.

Response of Core Microbial Consortia to Chronic Hydrocarbon Contaminations in Coastal Sediment Habitats

Traditionally, microbial surveys investigating the effect of chronic anthropogenic pressure such as polyaromatic hydrocarbons (PAHs) contaminations consider just the alpha and beta diversity and ignore the interactions among the different taxa forming the microbial community. Here, we investigated the ecological relationships between the three domains of life (i.e., Bacteria, Archaea, and Eukarya) using 454 pyrosequencing on the 16S rRNA and 18S rRNA genes from chronically impacted and pristine sediments, along the coasts of the Mediterranean Sea (Gulf of Lion, Vermillion coast, Corsica, Bizerte lagoon and Lebanon) and the French Atlantic Ocean (Bay of Biscay and English Channel). Our approach provided a robust ecological framework for the partition of the taxa abundance distribution into 859 core Operational taxonomic units (OTUs) and 6629 satellite OTUs. OTUs forming the core microbial community showed the highest sensitivity to changes in environmental and contaminant variations, with salinity, latitude, temperature, particle size distribution, total organic carbon (TOC) and PAH concentrations as main drivers of community assembly. The core communities were dominated by Gammaproteobacteria and Deltaproteobacteria for Bacteria, by Thaumarchaeota, Bathyarchaeota and Thermoplasmata for Archaea and Metazoa and Dinoflagellata for Eukarya. In order to find associations among microorganisms, we generated a co-occurrence network in which PAHs were found to impact significantly the potential predator – prey relationship in one microbial consortium composed of ciliates and Actinobacteria. Comparison of network topological properties between contaminated and non-contaminated samples showed substantial differences in the network structure and indicated a higher vulnerability to environmental perturbations in the contaminated sediments.