This dataset consists of metatranscriptomic sequencing reads corresponding to coastal micro-eukaryote communities sampled in Western Europe in 2018 and 2019.

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.

The ESA Sea State Climate Change Initiative (CCI) project has produced global multi-sensor time-series of along-track satellite synthetic aperture radar (SAR) integrated sea state parameters (ISSP) data from Sentinel-1 (referred to as SAR WV onboard Sentinel-1 Level 2P (L2P) ISSP data) with a particular focus for use in climate studies. This dataset contains the Sentinel-1 SAR Remote Sensing Integrated Sea State Parameter product (v1.0), which forms part of the ESA Sea State CCI version 3.0 release. This product provides along-track primary significant wave height measurements and secondary sea state parameters, calibrated with CMEMS model data and reference in situ measurements at 20km resolution every 100km, processed using the Pleskachevsky et. al., 2021 emprical model, separated per satellite and pass, including all measurements with flags and uncertainty estimates. These are expert products with rich content and no data loss. The SAR Wave Mode data used in the Sea State CCI SAR WV onboard Sentinel-1 Level 2P (L2P) ISSP v3 dataset come from the Sentinel-1 satellite missions spanning from 2014 to 2021 (Sentinel-1 A, Sentinel-1 B).

The SARWAVE project is developing a new sea state processor from SAR images to be applied over open ocean, sea ice, and coastal areas, and exploring potential synergy with other microwave and optical EO products.

Understanding the dynamics of species interactions for food (prey-predator, competition for resources) and the functioning of trophic networks (dependence on trophic pathways, food chain flows, etc.) has become a thriving ecological research field in recent decades. This empirical knowledge is then used to develop population and ecosystem modelling approaches to support ecosystem-based management. The TrophicCS data set offers spatialized trophic information on a large spatial scale (the entire Celtic Sea continental shelf and upper slope) for a wide range of species. It combines ingested prey (gut content analysis) and a more integrated indicator of food sources (stable isotope analysis). A total of 1337 samples of large epifaunal invertebrates (bivalve mollusks and decapod crustaceans), zooplankton, fish and cephalopods, corresponding to 114 species, were collected and analyzed for stable isotope analysis of their carbon and nitrogen content. Sample size varied between taxa (from 1 to 52), with an average of 11.72 individuals sampled per species, and water depths ranged from 57 to 516 m. The gut contents of 1026 fish belonging to ten commercially important species: black anglerfish (Lophius budegassa), white anglerfish (Lophius piscatorius), blue whiting (Micromesistius poutassou), cod (Gadus morhua), haddock (Melanogrammus aeglefinus), hake (Merluccius merluccius), megrim (Lepidorhombus whiffiagonis), plaice (Pleuronectes platessa), sole (Solea solea) and whiting (Merlangius merlangus) were analyzed. The stomach content data set contains the occurrence of prey in stomach, identified to the lowest taxonomic level possible. To consider potential ontogenetic diet changes, a large size range was sampled. The TrophicCS data set was used to improve understanding of trophic relationships and ecosystem functioning in the Celtic Sea. When you use the data in your publication, we request that you cite this data paper. If you use the present data set (TrophicCS) for the majority of the data analyzed in your study, you may wish to consider inviting at least one author of the core team of this data paper to become a collaborator /coauthor of your paper.

French benthic invertebrates composition and abundance taxa data are collected during monitoring surveys on the English Channel / Bay of Biscay coasts and Mediterranean coast (Quadrige program code : REBENT_FAU, RSL_FAU). Protocols are implemented in the Water Framework Directive.  Data are transmitted in a Seadatanet format (CDI + ODV) to EMODnet Biology european database. 498 ODV files have been generated from period 01/01/2003 to 31/12/2021.

We genotyped 1680 thornback ray Raja clavata sampled in the Bay of Biscay using a DNA chip described in Le Cam et al. (2019). After quality control 4604 SNPs were retained for identifying potential sex-linked SNPs using three methods: i) identification of excess of heterozygotes in one sex, ii) FST outlier analysis between the two sexes and iii) neuronal net modelling. Genotype coding: 0 homozygous for major allele, 1 heterozygous, 2 homozygous for minor allele. Flanking DNA sequences of SNPs identified with methods i) and ii) are also provided.  

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

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