In order to better characterize the population structure of common dolphins (Delphinus delphis) in the Bay of Biscay, a single digest RADseq (SbfI enzyme) protocol was used to obtain paired-end, 150bp NGS sequences on the Illumina NovaSeq 6000 sequencing platform. D. delphis samples from the Western North Atlantic, and samples from three other delphinid species were included as outgroups.

This visualization product displays the density of floating micro-litter per net normalized in grams per km² per year from specific protocols different from research and monitoring protocols. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Before 2021, there was no coordinated effort at the regional or European scale for micro-litter. Given this situation, EMODnet Chemistry proposed to adopt the data gathering and data management approach as generally applied for marine data, i.e., populating metadata and data in the CDI Data Discovery and Access service using dedicated SeaDataNet data transport formats. EMODnet Chemistry is currently the official EU collector of micro-litter data from Marine Strategy Framework Directive (MSFD) National Monitoring activities (descriptor 10). A series of specific standard vocabularies or standard terms related to micro-litter have been added to SeaDataNet NVS (NERC Vocabulary Server) Common Vocabularies to describe the micro-litter. European micro-litter data are collected by the National Oceanographic Data Centres (NODCs). Micro-litter map products are generated from NODCs data after a test of the aggregated collection including data and data format checks and data harmonization. A filter is applied to represent only micro-litter sampled according to a very specific protocol such as the Volvo Ocean Race (VOR) or Oceaneye. Densities were calculated for each net using the following calculation: Density (weight of particles per km²) = Micro-litter weight / (Sampling effort (km) * Net opening (cm) * 0.00001) When the weight of microlitters or the net opening was not filled, it was not possible to calculate the density. Percentiles 50, 75, 95 & 99 have been calculated taking into account data for all years. Warning: the absence of data on the map does not necessarily mean that they do not exist, but that no information has been entered in the National Oceanographic Data Centre (NODC) for this area.

'''DEFINITION''' Estimates of Ocean Heat Content (OHC) are obtained from integrated differences of the measured temperature and a climatology along a vertical profile in the ocean (von Schuckmann et al., 2018). The products used include three global reanalyses: GLORYS, C-GLORS, ORAS5 (GLOBAL_MULTIYEAR_PHY_ENS_001_031) and two in situ based reprocessed products: CORA5.2 (INSITU_GLO_PHY_TS_OA_MY_013_052) , ARMOR-3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012). Additionally, the time series based on the method of von Schuckmann and Le Traon (2011) has been added. The regional OHC values are then averaged from 60°S-60°N aiming i) to obtain the mean OHC as expressed in Joules per meter square (J/m2) to monitor the large-scale variability and change. ii) to monitor the amount of energy in the form of heat stored in the ocean (i.e. the change of OHC in time), expressed in Watt per square meter (W/m2). Ocean heat content is one of the six Global Climate Indicators recommended by the World Meterological Organisation for Sustainable Development Goal 13 implementation (WMO, 2017). '''CONTEXT''' Knowing how much and where heat energy is stored and released in the ocean is essential for understanding the contemporary Earth system state, variability and change, as the ocean shapes our perspectives for the future (von Schuckmann et al., 2020). Variations in OHC can induce changes in ocean stratification, currents, sea ice and ice shelfs (IPCC, 2019; 2021); they set time scales and dominate Earth system adjustments to climate variability and change (Hansen et al., 2011); they are a key player in ocean-atmosphere interactions and sea level change (WCRP, 2018) and they can impact marine ecosystems and human livelihoods (IPCC, 2019). '''CMEMS KEY FINDINGS''' Since the year 2005, the upper (0-700m) near-global (60°S-60°N) ocean warms at a rate of 0.6 ± 0.1 W/m2. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00234

This visualization product displays the cigarette related items abundance of marine macro-litter (> 2.5cm) per beach per year from non-MSFD monitoring surveys, research & cleaning operations without UNEP-MARLIN data. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Since the beginning of 2018, data of beach litter have been gathered and processed in the EMODnet Chemistry Marine Litter Database (MLDB). The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols and reference lists used on a European scale. Preliminary processings were necessary to harmonize all the data: - Exclusion of OSPAR 1000 protocol: in order to follow the approach of OSPAR that it is not including these data anymore in the monitoring; - Selection of surveys from non-MSFD monitoring, cleaning and research operations; - Exclusion of beaches without coordinates; - Selection of plastic bags related items only. The list of selected items is attached to this metadata. This list was created using EU Marine Beach Litter Baselines, the European Threshold Value for Macro Litter on Coastlines and the Joint list of litter categories for marine macro-litter monitoring from JRC (these three documents are attached to this metadata); - Exclusion of surveys without associated length; - Normalization of survey lengths to 100m & 1 survey / year: in some case, the survey length was not 100m, so in order to be able to compare the abundance of litter from different beaches a normalization is applied using this formula: Number of plastic bags related items of the survey (normalized by 100 m) = Number of plastic bags related items of the survey x (100 / survey length) Then, this normalized number of plastic bags related items is summed to obtain the total normalized number of plastic bags related items for each survey. Finally, the median abundance of plastic bags related items for each beach and year is calculated from these normalized abundances of plastic bags related items per survey. Percentiles 50, 75, 95 & 99 have been calculated taking into account plastic bags related items from other sources data for all years. More information is available in the attached documents. Warning: the absence of data on the map does not necessarily mean that they do not exist, but that no information has been entered in the Marine Litter Database for this area.

'''DEFINITION''' The temporal evolution of thermosteric sea level in an ocean layer is obtained from an integration of temperature driven ocean density variations, which are subtracted from a reference climatology to obtain the fluctuations from an average field. The products used include three global reanalyses: GLORYS, C-GLORS, ORAS5 (GLOBAL_MULTIYEAR_PHY_ENS_001_031) and two in situ based reprocessed products: CORA5.2 (INSITU_GLO_PHY_TS_OA_MY_013_052) , ARMOR-3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012). Additionally, the time series based on the method of von Schuckmann and Le Traon (2011) has been added. The regional thermosteric sea level values are then averaged from 60°S-60°N aiming to monitor interannual to long term global sea level variations caused by temperature driven ocean volume changes through thermal expansion as expressed in meters (m). '''CONTEXT''' The global mean sea level is reflecting changes in the Earth’s climate system in response to natural and anthropogenic forcing factors such as ocean warming, land ice mass loss and changes in water storage in continental river basins. Thermosteric sea-level variations result from temperature related density changes in sea water associated with volume expansion and contraction (Storto et al., 2018). Global thermosteric sea level rise caused by ocean warming is known as one of the major drivers of contemporary global mean sea level rise (Cazenave et al., 2018; Oppenheimer et al., 2019). '''CMEMS KEY FINDINGS''' Since the year 2005 the upper (0-2000m) near-global (60°S-60°N) thermosteric sea level rises at a rate of 1.3±0.3 mm/year. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00240

Moving 6-year analysis of Water body dissolved inorganic nitrogen in the NorthEast Atlantic for each season: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Every year of the time dimension corresponds to the 6-year centred average of each season. 6-year periods span from 1971/1976 until 2019/2024. Observation data span from 1971 to 2024. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0, 400.0, 500.0, 600.0, 700.0, 800.0, 900.0, 1000.0, 1100.0, 1200.0, 1300.0, 1400.0, 1500.0, 1750.0, 2000.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 400km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. Signal-to-noise ratio was fixed to 1 for vertical profiles and 0.1 for time series to account for the redundancy in the time series observations. A logarithmic transformation (DIVAnd.loglin) was applied to the data prior to the analysis. Background field: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

'''DEFINITION''' Significant wave height (SWH), expressed in metres, is the average height of the highest third of waves. This OMI provides global maps of the seasonal mean and trend of significant wave height (SWH), as well as time series in three oceanic regions of the same variables and their trends from 2002 to 2020, calculated from the reprocessed global L4 SWH product (WAVE_GLO_PHY_SWH_L4_MY_014_007). The extreme SWH is defined as the 95th percentile of the daily maximum SWH for the selected period and region. The 95th percentile is the value below which 95% of the data points fall, indicating higher than normal wave heights. The mean and 95th percentile of SWH (in m) are calculated for two seasons of the year to take into account the seasonal variability of waves (January, February and March, and July, August and September). Trends have been obtained using linear regression and are expressed in cm/yr. For the time series, the uncertainty around the trend was obtained from the linear regression, while the uncertainty around the mean and 95th percentile was bootstrapped. For the maps, if the p-value obtained from the linear regression is less than 0.05, the trend is considered significant. '''CONTEXT''' Grasping the nature of global ocean surface waves, their variability, and their long-term interannual shifts is essential for climate research and diverse oceanic and coastal applications. The sixth IPCC Assessment Report underscores the significant role waves play in extreme sea level events (Mentaschi et al., 2017), flooding (Storlazzi et al., 2018), and coastal erosion (Barnard et al., 2017). Additionally, waves impact ocean circulation and mediate interactions between air and sea (Donelan et al., 1997) as well as sea-ice interactions (Thomas et al., 2019). Studying these long-term and interannual changes demands precise time series data spanning several decades. Until now, such records have been available only from global model reanalyses or localised in situ observations. While buoy data are valuable, they offer limited local insights and are especially scarce in the southern hemisphere. In contrast, altimeters deliver global, high-quality measurements of significant wave heights (SWH) (Gommenginger et al., 2002). The growing satellite record of SWH now facilitates more extensive global and long-term analyses. By using SWH data from a multi-mission altimetric product from 2002 to 2020, we can calculate global mean SWH and extreme SWH and evaluate their trends, regionally and globally. '''KEY FINDINGS''' From 2002 to 2020, positive trends in both Significant Wave Height (SWH) and extreme SWH are mostly found in the southern hemisphere (a, b). The 95th percentile of wave heights (q95), increases faster than the average values, indicating that extreme waves are growing more rapidly than average wave height (a, b). Extreme SWH’s global maps highlight heavily storms affected regions, including the western North Pacific, the North Atlantic and the eastern tropical Pacific (a). In the North Atlantic, SWH has increased in summertime (July August September) but decreased in winter. Specifically, the 95th percentile SWH trend is decreasing by 2.1 ± 3.3 cm/year, while the mean SWH shows a decrease of 2.2 ± 1.76 cm/year. In the south of Australia, during boreal winter, the 95th percentile SWH is increasing at 2.6 ± 1.5 cm/year (c), with the mean SWH increasing by 0.5 ± 0.66 cm/year (d). Finally, in the Antarctic Circumpolar Current, also in boreal winter, the 95th percentile SWH trend is 3.2 ± 2.14 cm/year (c) and the mean SWH trend is 1.7 ± 0.84 cm/year (d). These patterns highlight the complex and region-specific nature of wave height trends. Further discussion is available in A. Laloue et al. (2024). '''DOI (product):''' https://doi.org/10.48670/mds-00352

EMODnet Chemistry aims to provide access to marine chemistry datasets and derived data products concerning eutrophication, acidity and contaminants. The importance of the selected substances and other parameters relates to the Marine Strategy Framework Directive (MSFD). This aggregated dataset contains all unrestricted EMODnet Chemistry data on potential hazardous substances, despite the fact that some data might not be related to pollution (e.g. collected by deep corer). Temperature, salinity and additional parameters are included when available. It covers the Northeast Atlantic Ocean (40W). Data were harmonised and validated by '‘IFREMER / IDM / SISMER - Scientific Information Systems for the SEA’ in France. The dataset contains water (profiles), sediment (profiles and timeseries) and biota (timeseries). The temporal coverage is 1974–2018 for water measurements, 1966–2022 for sediment measurements and 1979–2023 for biota measurements. Regional datasets concerning contaminants are automatically harvested and the resulting collections are harmonised and validated using ODV Software and following a common methodology for all sea regions ( https://doi.org/10.6092/8b52e8d7-dc92-4305-9337-7634a5cae3f4 ). Parameter names are based on P01 vocabulary, which relates to BODC Parameter Usage Vocabulary and is available at: https://vocab.nerc.ac.uk/search_nvs/P01/ . The harmonised dataset can be downloaded as as an ODV spreadsheet, which is composed of a metadata header followed by tab separated values. This spreadsheet can be imported into ODV Software for visualisation (more information can be found at: https://www.seadatanet.org/Software/ODV ). In addition, the same dataset is offered also as a txt file in a long/vertical format, in which each P01 measurement is a record line. Additionally, there are a series of columns that split P01 terms into subcomponents (substance, CAS number, matrix...).This transposed format is more adapted to worksheet applications (e.g. LibreOffice Calc).

This dataset contains all satellite altimeter wave heights above 9 m, from the following satellite missions: ERS-1, ERS-2, Topex-Poseidon (Topex only), Envisat, SARAL, Jason-1, Jason-2, Jason-3, Sentinel-3A, Sentinel-3B, Sentinel-6A, Cryosat-2, CFOSAT, SWOT. Storm event identification used the DetectHsStorm package developed by M. De Carlo and F. Ardhuin (  https://github.com/ardhuin/) . This data can be combined with modeled storm tracks (see F. Ardhuin, M. De Carlo, Storm tracks based on wave heights from LOPS WAVEWATCH III hindcast and ERA5 reanalysis, years 1991-2024, SEANOE (2025). doi: 10.17882/105148 )

This dataset comprises energy density and proximal composition (water, ash, lipid and protein contents) for anchovy (Engraulis encrasicolus) and sardine (Sardina pilchardus) from the Bay of Biscay, the English Channel and the southern North Sea between 2014 and 2017. Fish were sampled throughout various seasons thanks to the PELGAS (May, Bay of Biscay), EVHOE (October-November, Bay of Biscay), CGFS and CAMANOC (September-October, English Channel) surveys conducted by IFREMER on the RV “Thalassa”, to the JUVENA (September-October, Bay of Biscay) survey conducted by AZTI, and from commercial landings within the European Data Collection Framework (DCF) during the CAPTAIN project (France Filière Pêche). During the surveys, pelagic (PELGAS and JUVENA) or demersal (EVHOE and CGFS) trawl hauls are undertaken to identify species and measure individual fish traits. Professional sampling was performed from pelagic trawl or purse-seine catches. From the various surveys, a sub-sampling of the trawls was performed to cover as much as possible the spatial extent of the surveys along the french coast. From the various selected trawls, a sub-sampling of 5 fish per size class (when possible) was performed to cover the size range of each species, based on the following size classes : sardine (1 : <15 cm ; 2 : 15-20 cm ; 3 : >20cm), anchovy (1 : <10cm ; 2 : 10-14 cm ; 3 : >14cm). Each fish was individually measured to the nearest tenth of a centimeter and weighted to the nearest tenth of a gram. These measurements were taken either at sea or later in the laboratory. The collected fish were frozen individually at -20°C before laboratory processing. In the laboratory, maturity stages were determined following ICES guidelines (ICES, 2008) based on macroscopic gonads observation and using a six-stage key as follows: stages 1 & 2 indicate immature and developing individuals, stages 3–5 indicate three steps of increasing gonad development and the spawning period (stage 3: partial pre-spawning; stage 4: spawning (hydrated); stage 5: partial post-spawning), and stage 6 features the final post-spawning period. Fish characterised by maturity stages 3, 4 or 5 were considered as being in an active reproductive period as opposed to fish in stages 1, 2 or 6. Fish were then ground and freeze-dried during at least 48 hours. Water content of the entire fish was determined from dry mass and wet mass ratio. Then, fish were ground again to obtain fine homogeneous dry powder for subsequent analysis. Energy density measurements were performed following the protocols of Gatti et al. (2018). Two subsamples of fish powder were placed in an adiabatic bomb calorimeter (IKA C-4000 adiabatic bomb calorimeter, IKA-WerkeGmbh & co. KG) for energy measurements. The energy density (ED, kJ.g-1 dry mass) was determined by measuring the heat released through the combustion of a small subsample, approximately 200 mg. If the coefficient of variation between the two measurements exceeded 3%, a third measurement was made. Finally, ED subsamples measurements were averaged and assigned to each individual fish. Energy density analyses were conducted on 503 individuals for anchovy and 976 individuals for sardine. Ash content was determined gravimetrically by combusting dried tissue in a muffle furnace at 550°C for six hours. Lipids and proteins were analysed by a certified laboratory (Labocea, Plouzané, France). Protein content was estimated using the Kjeldahl method. It consists in first determining the quantity of nitrogen contained in the sample, and to convert it into protein content using a conversion factor (6.25 here). Lipid content was determined by hydrolysis, using petroleum ether as an organic solvent. Carbohydrates represent less than 1% of fish mass and were thus neglected. Protein, lipid and ash content did not exactly sum to 1 in DW (anchovy: mean = 0.91, sd = 0.04; sardine: mean = 0.90, sd = 0.04). This discrepancy may arise from residual water, measurement uncertainties, or to a lesser extent the exclusion of carbohydrates. Body component contents have been normalised by dividing each component by the sum of lipids, proteins and ash content, to sum to one, enabling comparisons between fishes, assuming proportional errors across the components. A total of 116 and 104 proximate composition analyses were performed for anchovy and sardine, respectively.