The upper ocean pycnocline (UOP) monthly climatology is based on the ISAS20 ARGO dataset containing Argo and Deep-Argo temperature and salinity profiles on the period 2002-2020. Regardless of the season, the UOP is defined as the shallowest significant stratification peak captured by the method described in Sérazin et al. (2022), whose detection threshold is proportional to the standard deviation of the stratification profile. The three main characteristics of the UOP are provided -- intensity, depth and thickness -- along with hydrographic variables at the upper and lower edges of the pycnocline, the Turner angle and density ratio at the depth of the UOP. A stratification index (SI) that evaluates the amount of buoyancy required to destratify the upper ocean down to a certain depth, is also included. When evaluated at the bottom of the UOP, this gives the upper ocean stratification index (UOSI) as discussed in Sérazin et al. (2022). Three mixed layer depth variables are also included in this dataset, including the one using the classic density threshold of 0.03 kg.m-3, along with the minimum of these MLD variables. Several statistics of the UOP characteristics and the associated quantities are available in 2°×2° bins for each month of the year, whose results were smoothed using a diffusive gaussian filter with a 500 km scale. UOP characteristics are also available for each profile, with all the profiles sorted in one file per month.

The GEBCO_2022 Grid is a global continuous terrain model for ocean and land with a spatial resolution of 15 arc seconds. In regions outside of the Arctic Ocean area, the grid uses as a base Version 2.4 of the SRTM15_plus data set (Tozer, B. et al, 2019). This data set is a fusion of land topography with measured and estimated seafloor topography. Included on top of this base grid are gridded bathymetric data sets developed by the four Regional Centers of The Nippon Foundation-GEBCO Seabed 2030 Project. The GEBCO_2022 Grid represents all data within the 2022 compilation. The compilation of the GEBCO_2022 Grid was carried out at the Seabed 2030 Global Center, hosted at the National Oceanography Centre, UK, with the aim of producing a seamless global terrain model. Outside of Polar regions, the Regional Centers provide their data sets as sparse grids i.e. only grid cells that contain data are populated. These data sets were included on to the base using a remove-restore blending procedure. This is a two-stage process of computing the difference between the new data and the base grid and then gridding the difference and adding the difference back to the existing base grid. The aim is to achieve a smooth transition between the new and base data sets with the minimum of perturbation of the existing base data set. The data sets supplied in the form of complete grids (primarily areas north of 60N and south of 50S) were included using feather blending techniques from GlobalMapper software. The GEBCO_2022 Grid has been developed through the Nippon Foundation-GEBCO Seabed 2030 Project. This is a collaborative project between the Nippon Foundation of Japan and the General Bathymetric Chart of the Oceans (GEBCO). It aims to bring together all available bathymetric data to produce the definitive map of the world ocean floor by 2030 and make it available to all. Funded by the Nippon Foundation, the four Seabed 2030 Regional Centers include the Southern Ocean - hosted at the Alfred Wegener Institute, Germany; South and West Pacific Ocean - hosted at the National Institute of Water and Atmospheric Research, New Zealand; Atlantic and Indian Oceans - hosted at the Lamont-Doherty Earth Observatory, Columbia University, USA; Arctic and North Pacific Oceans - hosted at Stockholm University, Sweden and the Center for Coastal and Ocean Mapping at the University of New Hampshire, USA.

In recent years, large datasets of in situ marine carbonate system parameters (partial pressure of CO2 (pCO2), total alkalinity, dissolved inorganic carbon and pH) have been collated. These carbonate system datasets have highly variable data density in both space and time, especially in the case of pCO2, which is routinely measured at high frequency using underway measuring systems. This variation in data density can create biases when the data are used, for example for algorithm assessment, favouring datasets or regions with high data density. A common way to overcome data density issues is to bin the data into cells of equal latitude and longitude extent. This leads to bins with spatial areas that are latitude and projection dependent (eg become smaller and more elongated as the poles are approached). Additionally, as bin boundaries are defined without reference to the spatial distribution of the data or to geographical features, data clusters may be divided sub-optimally (eg a bin covering a region with a strong gradient). To overcome these problems and to provide a tool for matching in situ data with satellite, model and climatological data, which often have very different spatiotemporal scales both from the in situ data and from each other, a methodology has been created to group in situ data into ‘regions of interest’, spatiotemporal cylinders consisting of circles on the Earth’s surface extending over a period of time. These regions of interest are optimally adjusted to contain as many in situ measurements as possible. All in situ measurements of the same parameter contained in a region of interest are collated, including estimated uncertainties and regional summary statistics. The same grouping is done for each of the other datasets, producing a dataset of matchups. About 35 million in situ datapoints were then matched with data from five satellite sources and five model and re-analysis datasets to produce a global matchup dataset of carbonate system data, consisting of 287,000 regions of interest spanning 54 years from 1957 to 2020. Each region of interest is 100 km in diameter and 10 days in duration. An example application, the reparameterisation of a global total alkalinity algorithm, is shown. This matchup dataset can be updated as and when in situ and other datasets are updated, and similar datasets at finer spatiotemporal scale can be constructed, for example to enable regional studies. This dataset was funded by ESA Satellite Oceanographic Datasets for Acidification (OceanSODA) project which aims at developing the use of satellite Earth Observation for studying and monitoring marine carbonate chemistry.

Reef-building species are recognized as having an important ecological role and as generally enhancing the diversity of benthic organisms in marine habitats.  However, although these ecosystem engineers have a facilitating role for some species, they may exclude or compete with others. The honeycomb worm Sabellaria alveolata (Linnaeus, 1767) is an important foundation species, commonly found from northwest Ireland to northern Mauritania (Curd et al., 2020), whose reef structures increase the physical complexity of the marine benthos, supporting high levels of biodiversity. Local patterns and regional differences in taxonomic and functional diversity were examined in honeycomb worm reefs from ten sites along the northeastern Atlantic to explore variation in diversity across biogeographic regions and the potential effects of environmental drivers. To characterize the functional diversity at each site, a biological trait analysis (BTA) was conducted (Statzner et al., 1994). Here we present the functional trait database used for the benthic macrofauna found to live in association with honeycomb worm reefs. Eight biological traits (divided into 32 modalities) were selected (Table 1), providing information linked to the ecological functions performed by the associated macrofauna. The selected traits provide information on: (i) resource use and availability (by the trophic group of species, e.g. Thrush et al. 2006); (ii) secondary production and the amount of energy and organic matter (OM) produced based on the life cycle of the organisms (including longevity, maximum size and mode of reproduction, e.g. (Cusson and Bourget, 2005; Thrush et al., 2006) and; (iii) the behavior of the species in general [i.e. how these species occupy the environment and contribute to biogeochemical fluxes through habitat, movement, and bioturbation activity at different bathymetric levels, e.g. (Solan et al., 2004; Thrush et al., 2006; Queirós et al., 2013). Species were scored for each trait modality based on their affinity using a fuzzy coding approach (Chevenet et al., 1994), where multiple modalities can be attributed to a species if appropriate, and allowed for the incorporation of intraspecific variability in trait expression. The information concerning polychaetes was derived primarily from Fauchald et al (1979) and Jumars et al (2015). Information on other taxonomic groups was obtained either from databases of biological traits (www.marlin.ac.uk/biotic) or publications (Naylor, 1972; King, 1974; Caine, 1977; Lincoln, 1979; Holdich and Jones, 1983; Smaldon et al., 1993; Ingle, 1996; San Martín, 2003; Southward, 2008; Gil, 2011; Leblanc et al., 2011; Rumbold et al., 2012; San Martín and Worsfold, 2015; Jones et al., 2018). Map indicating the locations of the 10 study sites in the UK, France and Portugal within the four biogeographic provinces defined by Dinter (2001). (All sites were sampled in 8 different stations, except for UK4 where 5 stations were sampled).

WGS for Iatlantic projet ( ) for assessing past and present connectivity

The diet and stable isotopic (i.e. δ15N and δ13C values) compositions of eels have been studied during each season of 2019 with a fyke net in six estuaries located along the French coast of the eastern English Channel (Slack, Wimereux, Liane, Canche, Authie and Somme estuaries) (10.1371/journal.pone.0270348).

Raw reads for the assembly of Gambusia holbrooki genome.

The Mediterranean Sea is generally described as an oligotrophic area where primary productivity is limited to a few coastal environments with nutrient-enriched fluvial input. However, several studies have revealed that the hydrology of the western Mediterranean has major seasonal productive patterns linked either to significant riverine input or to seasonal upwelling cells. This study aims to: i) discuss organic microfossils (i.e. pollen and dinoflagellate cyst assemblages, as well as other non-pollen palynomorphs) from two different productive areas of the western Mediterranean Sea, and ii) examine the importance of the interconnections between marine and continental influences responsible for modern palynomorph distributions. Based on 25 samples from the Gulf of Lion (GoL) and Algerian Margin, this study key findings are: i) that GoL marine productivity is driven by the combination of discharges from the Rhône River and seasonal upwelling mechanisms, ii) that the strong productive pattern of the northern African coast is driven by water density front mixings and related upwellings. These two patterns are discussed in the light of major links that provide a better understanding of the signatures of marine and continental bio-indicators. The dinocyst Lingulodinium machaerophorum can be considered as a tracer of Rhône River plume influence in the GoL. Brigantedinium taxa are shown to be upwelling-sensitive in both studied areas. Typical differences in vegetation across the north–south climate gradient in the western Mediterranean Basin are highlighted by the larger ratio of Euro-Siberian to Mediterranean pollen taxa in the northern sector. Synoptic maps also illustrate the complex interactions of environmental drivers determining the distributions of continental and marine palynomorphs in the western Mediterranean Sea.

We developed a panel of single nucleotide polymorphism (SNP) markers for thornback ray Raja clavata using a RADSeq protocole. Demultiplexed sequences were aligned to the genome of Leucoraja erinacea which was used as reference genome. From an initial set of 389 483 putative SNPs, 7741 SNPs with the largest minor allele frequency were selected for implementation on an Infinium® XT iSelect-96 SNP-array implemented by LABOGENA DNA. For the array, SNPs [T/C] and [T/G] were replaced by those from the complementary strand [A/G] and [A/C] respectively. For some SNPs, a second SNP was found in the 50 nucleotide bases flanking sequence. In these cases, two SNP probes were developed with each of the two alleles of the second SNP. A SNP probe naming convention was adopted to identify these pairs of probes corresponding to the same SNP locus: “MAJ” or “MIN” followed by the corresponding base was included in the probe name. For some of these pairs, only one of the two markers could be developed, resulting in a total set of 9120 SNP probes, including 6360 single SNP probes, 10 MAJ or MIN probes for which a single probe was successfully developed, and 1375 pairs of probes with MAJ and MIN versions. The 9120 SNP genotypes were then scored using the clustering algorithm implemented in the Illumina® GenomeStudio Genotyping Analysis Module v2.0.3 for 7726 individual samples, including duplicates, mostly from the Bay of Biscay but also from the Mediterranean Sea and West Iberia. Overall, 1643 SNPs failed to be genotyped in all individuals, for 319 markers the minor allele was not found and 7158 markers (including 1974 for 987 MIN-MAJ pairs) produced bi-allelic genotypes. The majority of these SNPs had a minor allele frequency between 0.1 and 0.5. The MIN-MAJ probes can be used for quality checking the genotyping results

French Zostera Marina et Zostera Noltei abundance data are collected during monitoring surveys on the English Channel / Bay of Biscay coasts. Protocols are impletmented in the Water Framework Directive. Data are transmitted in a Seadatanet format (CDI + ODV) to EMODnet Biology european database. 35 ODV files have been generated from period 01/01/2004 to 31/12/2021 for Z. Marina and from 01/01/2011 to 31/12/2021 for Z. Noltei.