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

The CDR-derived Wet Tropospheric Correction (WTC) Product V2 is generated from the Level-2+ along-track altimetry products version 2024 (L2P 2024) distributed by AVISO+ (www.aviso.altimetry.fr). It provides a long-term, homogenized estimation of the wet tropospheric correction based on Climate Data Records (CDRs) of atmospheric water vapour combined with high frequencies MWR data. Two independent CDRs datasets are used: - REMSS V7R2 (coverage until 2022) https://www.remss.com/measurements/atmospheric-water-vapor/tpw-1-deg-product/ - HOAPS V5 precursor CDR from EUMETSAT CM SAF (coverage until 2020) HOAPS V4/V5 data available via https://wui.cmsaf.eu Note: the HOAPS V5 precursor is not yet an official CM SAF product; full validation and public release are pending. The MWR/CDR WTC V2 estimates is derived using spatially varying but temporally constant polynomial coefficients (ai). 1. WTC V2 – Along-track L2P Product Data format: The WTC V2 product is delivered in Level-2+ (L2P) format, along the satellite ground track. Each mission is distributed as a compressed archive (.tar.gz) containing one NetCDF4 CF-1.8 file per mission cycle. Archive naming convention: <mission>_WTC_from_WV_CDR_<version>.tar.gz mission: TP (TOPEX/Poseidon), J1, J2, J3 version: product version (currently V2) File naming convention inside archives: <mission>_C<cycle>.nc cycle: 4-digit cycle index (e.g., C0001) Each NetCDF file contains: 1/ Along-track WTC estimate; 2/ Ancillary information; 3/ Space–time coordinates 2. WTC CDR Uncertainties – Gridded Product: A complementary product is provided, delivering regional trend estimates and associated uncertainties from the WTC Climate Data Record. The uncertainty product is distributed as a single NetCDF4 file: wtc_trend_uncertainties.nc . This file contains global gridded fields of WTC CDR trend and uncertainty parameters. Product content: This is the first dedicated version providing both: WTC CDR (HOAPS) linear trends, and Uncertainty estimates on these trends. Uncertainties are expressed as 1-sigma confidence intervals, and propagated using the methodology described in Section 2.3 of the Product User Manual. The product includes: - Total uncertainty on the WTC trend, propagated from all identified uncertainty sources in the WTC–TCWV regression. - Individual contributions of uncertainty sources (Uncertainties on regression coefficients: a0, a1 and their standard deviations; Uncertainties inherited from the HOAPS TCWV CDR) These fields enable users to assess the relative importance of each uncertainty component and recompute uncertainty propagation with alternative methods. Included regression input variables: To ensure transparency and reproducibility, the product provides: 1/ regression coefficients a0, a1; 2/ their associated uncertainties (std of a0, std of a1); 3/additional diagnostic fields required to recompute uncertainties if needed.

The SWISCA Level 2S (L2S) product provides along-track colocation of SWIM wave measuring instrument onto SCAT scatterometer grid, over the global ocean. SWIM and SCAT are both instruments onboard CFOSAT. CFOSAT (Chinese French Ocean SATellite) is a french-chinese mission launched in 2018, whose aim is to provide wind (SCAT instrument) and wave (SWIM instrument) measurements over the sea surface. The SWISCA L2S product is broken down into three different subproducts: - L2A containing the sigma0 of both SWIM and SCAT - L2B containing the wave parameters measured by SWIM and wind vectors measured by SCAT - AUX containing additional ancillary fields such as sea ice concentration (from CERSAT/SSMI), ocean currents (from CMEMS/GlobCurrent), SST and Wind (from ECMWF), rain rate (from IMERG), and WaveWatch3 wave spectra. All SWIM and ancillary observations are resampled onto SCAT scatterometer's geometry (wind vector cells, WVC). The SWISCA level 2S product is generated in delayed mode, a few days after acquisition. It is intended to foster cross analysis of SWIM and SCAT observations, and their combination to improve the retrieval of both wind and wave parameters. The SWISCA L2S product is generated and distributed by Ifremer / CERSAT in the frame of the Ifremer Wind and Wave Operation Center (IWWOC) co-funded by Ifremer and CNES and dedicated to the processing of the delayed mode data of CFOSAT mission.

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

Global wave hindcast (1961-2020) at 1° resolution using CMIP6 wind and sea-ice forcings for ALL (historical), GHG (historical greenhouse-gas-only), AER (historical Anthropogenic-aerosol-only), NAT (historical natural only) scenario.

Raw reads for the assembly of Gambusia holbrooki genome.

The ClimateFish database collates abundance data of 15 fish species proposed as candidate indicators of climate change in the Mediterranean Sea. An initial group of eight Mediterranean indigenous species (Epinephelus marginatus, Thalassoma pavo, Sparisoma cretense, Coris julis, Sarpa salpa, Serranus scriba, Serranus cabrilla and Caranx crysos) with wide distribution, responsiveness to temperature conditions and easy identification were selected by a network of Mediterranean scientists joined under the CIESM programme ‘Tropical Signals’ (https://www.ciesm.org/marine/programs/tropicalization.htm; Azzurro et al. 2010). Soon after, and thanks to the discussion with other expert groups and projects, C. crysos was no longer considered, and Lessepsian fishes (Red Sea species entering the Mediterranean through the Suez Canal) were included, namely: Fistularia commersonii, Siganus luridus, Siganus rivulatus, Pterois miles, Stephanolopis diaspros, Parupeneus forskali, Pempheris rhomboidea and Torquigener flavimaculosus. Considering the trend of increase of these species in the Mediterranean Sea (Golani et al. 2021) and their projected distribution according to climate change scenarios (D’Amen and Azzurro, 2020), more data on these tropical invaders are expected to come in the future implementation of the study. Data were collected according to a simplified visual census methodology (Garrabou et al. 2019) along standard transects of five minutes performed at a constant speed of 10m/min, corresponding approximately to an area of 50x5m. Four different depth layers were surveyed:  0-3m, 5-10 m, 11-20 m, 21-30 m. So far, the ClimateFish database includes fish counts collected along 3142 transects carried out in seven Mediterranean countries between 2009 and 2021, for a total number of 101'771 observed individuals belonging to the 15 fish species. Data were collected by a large team of researchers which joined in a common monitoring strategy supported by different international projects, which are acknowledged below. This database, when associated with climate data, offers new opportunities to investigate spatio-temporal effects of climate change in the Mediterranean Sea and test the effectiveness of each species as a possible climate change indicator.   Contacts: ernesto.azzurro(at)cnr.it   References: Azzurro E., Maynou F., Moschella P. (2010). A simplified visual census methodology to detect variability trends of coastal mediterranean fishes under climate change scenarios. Rapp. Comm. int. Mer Médit., 39. D’Amen, M. and Azzurro, E. (2020). Lessepsian fish invasion in Mediterranean marine protected areas: a risk assessment under climate change scenarios. ICES Journal of Marine Science, 77(1), pp.388-397. Garrabou, J., Bensoussan, N., Azzurro, E. (2019). Monitoring climate-related responses in Mediterranean marine protected areas and beyond: five standard protocols. Golani D.,  Azzurro E.,  Dulčić J.,  Massutí E., Orsi-Relini L.  (2021).  Atlas of Exotic Fishes in the Mediterranean Sea.  2nd edition  [F. Briand, Ed.]  365 pages.  CIESM Publishers, Paris, Monaco. ISBN number  978-92-990003-5-9   

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**

New results acquired in south-Brittany (MD08-3204 CQ core: Bay of Quiberon and VK03-58bis core: south Glénan islands) allow depicting Holocene paleoenvironmental changes from 8.5 ka BP to present through a multi-proxy dataset including sedimentological and palynological data. First, grain-size analyses and AMS-14C dates highlight a common sedimentary history for both study cores. The relative sea level (RSL) slowdown was accompanied by a significant drop of the sedimentation rates between ca. 8.3 and 5.7 ka BP, after being relatively higher at the onset of the Holocene. This interval led to the establishment of a shell-condensed level, identified in core VK03-58bis by the “Turritella layer” and interpreted as a marker for the maximum flooding surface. Palynological data (pollen grains and dinoflagellate cyst assemblages) acquired in core MD08-3204 CQ argue for an amplification of the fluvial influence since 5.7 ka BP; the establishment of the highstand system tract (i.e., mixed marine and fluviatile influences on the platform) then accompanying the slowdown of the RSL rise-rates. On the shelf, the amplification of Anthropogenic Pollen Indicators (API) is then better detected since 4.2 ka BP, not only due to human impact increase but also due to a stronger fluvial influence on the shelf during the Late Holocene. Palynological data, recorded on the 8.5–8.3 ka BP interval along an inshore-offshore gradient, also demonstrate the complexity of the palynological signal such as i) the fluvial influence that promotes some pollinic taxa (i.e., Corylus, Alnus) from proximal areas and ii) the macro-regionalization of palynomorph sources in distal cores. In addition, the comparison of palynological tracers, including API, over the last 7 kyrs, with south-Brittany coastal and mid-shelf sites subjected to northern vs. southern Loire catchment areas, allowed discussing a major hydro-climatic effect on the reconstructed palynological signals. Strengthened subpolar gyre dynamics (SPG), combined with recurrent positive North Atlantic Oscillation (NAO) configurations, appear responsible for increased winter precipitations and fluvial discharges over northern Europe, such as in Brittany. Conversely, weakened SPG intervals, associated with negative NAO-like modes, are characterized by intensified winter fluvial discharges over southern Europe. Interestingly, we record, at an infra-orbital timescale, major peaks of API during periods of strengthened (/weakened) SPG dynamics in sites subjects to Brittany watersheds (/Loire watersheds) inputs.

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.