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.

This dataset contains the pictures used for morphometric measurements, as well as the elemental compositon and production rates data, of planktonic Rhizaria. Specimens were collected in the bay of Villefranche-sur-Mer in May 2019 and during the P2107 cruise in the California Current in July-August 2021. Analyses of the data can be found at https://github.com/MnnLgt/Elemental_composition_Rhizaria.

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

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.

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

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.  

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.

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