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

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).

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).

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. **This version is now superseded by the version 4 with higher spatial and temporal resolution**

Particularly suited to the purpose of measuring the sensitivity of benthic communities to trawling, a trawl disturbance indicator (de Juan and Demestre, 2012, de Juan et al. 2009) was proposed based on benthic species biological traits to evaluate the sensibility of mega- and epifaunal community to fishing pressure known to have a physical impact on the seafloor (such as dredging and bottom trawling). The selected biological traits were chosen as they determine vulnerability to trawling: mobility, fragility, position on substrata, average size and feeding mode that can easily be related to the fragility, recoverability and vulnerability ecological concepts. The five categories retained are functional traits that were selected based on the knowledge of the response of benthic taxa to trawling disturbance (de Juan et al., 2009). They reflect respectively the possibility to avoid direct gear impact, to benefit from trawling for feeding, to escape gear, to get caught by the net and to resist trawling/dredging action, each of these characteristics being either advantageous or sensitive to trawling. To expand this approach to that proposed by Certain et al. (2015), the protection status of certain species was also indicated. To enable quantitative analysis, a score was assigned to each category: from low sensitivity (0) to high sensitivity (3). Biological traits of species have been defined, from the BIOTIC database (MARLIN, 2014) and from information given by Garcia (2010), Le Pape et al. (2007) and Brind’Amour et al. (2009). For missing traits, additional information from literature has been considered. The protection status of each taxa was also scored: Atlantic species listed in OSPAR List of Threatened and/or Declining Species and Habitats (https://www.ospar.org/work-areas/bdc/species-habitats/list-of-threatened-declining-species-habitats) and Mediterranean species listed in Vulnerable Marine Ecosystems (FAO, 2018 and Oceana, 2017) were scored 3 and other species were scored 1. The scores of 1085 taxa commonly found in bottom trawl by-catch in the southern North Sea, English Channel and north-western Mediterranean were described.

French intertidal and subtidal Macroalgae taxa data are collected during monitoring surveys on the English Channel / Bay of Biscay coasts.  Protocols are implemented in the Water Framework Directive. Data are transmitted in a Seadatanet format (CDI + ODV) to EMODnet Biology european database. 131 ODV files have been generated from period 01/01/2006 to 31/12/2021.

This data set is related to the article "Improving the robustness of dissipation rate estimates from microstructure shear data processing in ocean turbulence" (submitted to JTech, AMS). It provides the raw data files (with extension .p) from a vertical microstructure profiler VMP-6000 used for the published study. Raw data files are provided since the study precisely report on the data processing of raw microstructure horizontal velocity shear data to get dissipation rates of turbulent kinetic energy using the manufacturer software. The software (ODAS) can be downloaded from the manufacturer website at: https://rocklandscientific.com. The profiles were collected in the Mediterranean Sea in 2013 (French VAD cruise) and 2014 (Italian MEDOCC cruise), and in the North Atlantic Ocean in 2021 (MoMAR cruise). More details on the profiles are given in the related publication.

The Programme Ocean Multidisciplinaire Meso Echelle (POMME) was designed to describe and quantify the role of mesoscale processes in the subduction of mode waters in the Northeast Atlantic. Intensive situ measurements were maintained during 1 year (September 2000 - October 2001), over a 8 degrees square area centered on 18 degrees W, 42 degrees N. In order to synthesized the in-situ physical observations, and merge them with satellite altimetry and surface fluxes datasets, a simplified Kalman filter has been designed. Daily fields of temperature, salinity, and stream function were produced on a regular grid over a full seasonal cycle. We propose here the gridded fields (KA_ files) and the in-situ datasets used by the analysis (Data_ files).

This dataset gathers isotopic ratios (carbon and nitrogen) and concentrations of both priority (mercury species and polychlorinated biphenyls congeners) and emerging (musks and sunscreens) micropollutants measured in a host-parasite couple (hake Merluccius merluccius muscle and in its parasite Anisakis sp) from the south of Bay of Biscay in 2018. In addition, the hake infection degree measured as the number of Anisakis sp. larvae was added for each hake collected.

The Level 4 merged microwave wind product is a complete set of hourly global 10-m wind maps on a 0.25x0.25 degree latitude-longitude grid, spanning 1 Jan 2010 through the end of 2020. The product combines background neutral equivalent wind fields from ERA5, daily surface current fields from CMEMS, and stress equivalent winds obtained from several microwave passive and active sensors to produce hourly surface current relative stress equivalent wind analyses. The satellite winds include those from recently launched L-band passive sensors capable of measuring extreme winds in tropical cyclones, with little or no degradation from precipitation. All satellite winds used in the analyses have been recalibrated using a large set of collocated satellite-SFMR wind data in storm-centric coordinates. To maximize the use of the satellite microwave data, winds within a 24-hour window centered on the analysis time have been incorporated into each analysis. To accomodate the large time window, satellite wind speeds are transformed into deviations from ERA5 background wind speeds interpolated to the measurement times, and then an optical flow-based morphing technique is applied to these wind speed increments to propagate them from measurement to analysis time. These morphed wind speed increments are then added to the background wind speed at the analysis time to yield a set of total wind speeds fields for each sensor at the analysis time. These individual sensor wind speed fields are then combined with the background 10-m wind direction to yield vorticity and divergence fields for the individual sensor winds. From these, merged vorticity and divergence fields are computed as a weighted average of the individual vorticity and divergence fields. The final vector wind field is then obtained directly from these merged vorticity and divergence fields. Note that one consequence of producing the analyses in terms of vorticity and divergence is that there are no discontinuities in the wind speed fields at the (morphed) swath edges. There are two important points to be noted: the background ERA5 wind speed fields have been rescaled to be globally consistent with the recalibrated AMSR2 wind speeds. This rescaling involves a large increase in the ERA5 background winds beyond about 17 m/s. For example, an ERA5 10 m wind speed of 30 m/s is transformed into a wind speed of 41 m/s, and a wind speed of 34 m/s is transformed into a wind speed of about 48 m/s. Besides the current version of the product is calibrated for use within tropical cyclones and is not appropriate for use elsewhere. This dataset was produced in the frame of ESA MAXSS project. The primary objective of the ESA Marine Atmosphere eXtreme Satellite Synergy (MAXSS) project is to provide guidance and innovative methodologies to maximize the synergetic use of available Earth Observation data (satellite, in situ) to improve understanding about the multi-scale dynamical characteristics of extreme air-sea interaction.