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.

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

Wave impact is the primary cause of coastal structure failure. While wave impact is widely studied in controlled environments, in situ measurements of wave impact pressure are rare. The results of a campaign to measure wave impact pressure in situ are summarised here. Data were collected from 2016 to 2019 from anchored pressure gauges on the wall of the Artha breakwater in southwestern France. The acquisition frequency is 10 kHz and 10-minute bursts are recorded every hour. Two databases are published, one by burst and one by impact. The burst database summarises the main parameters describing the 10-minute record, while the impact database contains a list of parameters describing each impact.

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

A prerequisite for a successful development of a multi-mission wind dataset is to ensure good inter-calibration of the different extreme wind datasets to be integrated in the product. Since the operational hurricane community is working with the in-situ dropsondes as wind speed reference, which are in turn used to calibrate the NOAA Hurricane Hunter Stepped Frequency Microwave Radiometer (SFMR) wind data, MAXSS has used the latter to ensure extreme-wind inter-calibration among the following scatterometer and radiometer systems: the Advanced Scatterometers onboard the Metop series (i.e., ASCAT-A, -B, and -C), the scatterometers onboard Oceansat-2 (OSCAT) and ScatSat-1 (OSCAT-2), and onboard the HY-2 series (HSCAT-A, -B); the Advanced Microwave Scanning Radiometer 2 onboard GCOM-W1(AMSR-2), the multi-frequency polarimetric radiometer (Windsat), and the L-band radiometers onboard the Soil Moisture and Ocean Salinity (SMOS) and the Soil Moisture Active Passive (SMAP) missions. In summary, a two-step strategy has been followed to adjust the high and extreme wind speeds derived from the mentioned scatterometer and radiometer systems, available in the period 2009-2020. First, the C-band ASCATs have been adjusted against collocated storm-motion centric SFMR wind data. Then, both SFMR winds and ASCAT adjusted winds have been used to adjust all the other satellite wind systems. In doing so, a good inter-calibration between all the systems is ensured not only under tropical cyclone (TC) conditions, but also elsewhere. This dataset was produced in the frame of the ESA funded Marine Atmosphere eXtreme Satellite Synergy (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.

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.

In order to better characterize the genetic diversity of Cetaceans and especially the common Dolphin from the Bay of Biscay, sequences from the mitochondrial Cytochrome B region were obtained from water samples acquired close to groups of dolphins.

The ESA Sea State Climate Change Initiative (CCI) project has produced global multi-sensor time-series of along-track satellite synthetic aperture radar (SAR) significant wave height (SWH) data (referred to as SAR WV onboard Sentinel-1 Level 2P (L2P) SWH data) with a particular focus for use in climate studies. This dataset contains the Sentinel-1 SAR Remote Sensing Significant Wave Height product (version 1.0), which is part of the ESA Sea State CCI Version 3.0 release. This product provides along-track SWH measurements at 20km resolution every 100km, processed using the Quach et al statistical model , separated per satellite and pass, including all measurements with flags, corrections and extra parameters from other sources. These are expert products with rich content and no data loss. The SAR Wave Mode data used in the Sea State CCI dataset v3 come from Sentinel-1 satellite missions spanning from 2015 to 2021 (Sentinel-1 A, Sentinel-1 B).