Dart Seq data gathered on Blue Shark in the framework of the PSTBS-IO project supported by funding from FAO, CSIRO Oceans and Atmosphere, AZTI Tecnalia, Institut de recherche pour le développement (IRD), and Research Institute for Tuna Fisheries (RITF) and financial assistance of the European Union (GCP/INT/233/EC – Population structure of IOTC species in the Indian Ocean), and POPSIZE project supported by FEAMP (2014-2020 UE N°508/2014), and Institut français de recherche pour l'Exploitation de la mer (Ifremer).

This visualization product displays the single use plastics (SUP) related items abundance of marine macro-litter (> 2.5cm) per beach per year from Marine Strategy Framework Directive (MSFD) monitoring surveys. 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 processing 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 MSFD surveys only (exclusion of other monitoring, cleaning and research operations); - Exclusion of beaches without coordinates; - Selection of SUP related items only. The list of selected items is attached to this metadata. This list was created using EU Marine Beach Litter Baselines for Macro Litter on Coastlines from JRC (this document is attached to this metadata); - Normalization of survey lengths to 100m & 1 survey / year: in some case, the survey length was not exactly 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 SUP items of the survey (normalized by 100 m) = Number of SUP related items of the survey x (100 / survey length) Then, this normalized number of¨SUP related items is summed to obtain the total normalized number of SUP related items for each survey. Finally, the median abundance of SUP related items for each beach and year is calculated from these normalized abundances of SUP related items per survey. Sometimes the survey length was null or equal to 0. Assuming that the MSFD protocol has been applied, the length has been set at 100m in these cases. Percentiles 50, 75, 95 & 99 have been calculated taking into account SUP related items from MSFD data for all years. More information is available in the attached documents. Warning: the absence of data on the map doesn't necessarily mean that they don't exist, but that no information has been entered in the Marine Litter Database for this area.

'''DEFINITION''' We have derived an annual eutrophication and eutrophication indicator map for the North Atlantic Ocean using satellite-derived chlorophyll concentration. Using the satellite-derived chlorophyll products distributed in the regional North Atlantic CMEMS MY Ocean Colour dataset (OC- CCI), we derived P90 and P10 daily climatologies. The time period selected for the climatology was 1998-2017. For a given pixel, P90 and P10 were defined as dynamic thresholds such as 90% of the 1998-2017 chlorophyll values for that pixel were below the P90 value, and 10% of the chlorophyll values were below the P10 value. To minimise the effect of gaps in the data in the computation of these P90 and P10 climatological values, we imposed a threshold of 25% valid data for the daily climatology. For the 20-year 1998-2017 climatology this means that, for a given pixel and day of the year, at least 5 years must contain valid data for the resulting climatological value to be considered significant. Pixels where the minimum data requirements were met were not considered in further calculations. We compared every valid daily observation over 2021 with the corresponding daily climatology on a pixel-by-pixel basis, to determine if values were above the P90 threshold, below the P10 threshold or within the [P10, P90] range. Values above the P90 threshold or below the P10 were flagged as anomalous. The number of anomalous and total valid observations were stored during this process. We then calculated the percentage of valid anomalous observations (above/below the P90/P10 thresholds) for each pixel, to create percentile anomaly maps in terms of % days per year. Finally, we derived an annual indicator map for eutrophication levels: if 25% of the valid observations for a given pixel and year were above the P90 threshold, the pixel was flagged as eutrophic. Similarly, if 25% of the observations for a given pixel were below the P10 threshold, the pixel was flagged as oligotrophic. '''CONTEXT''' Eutrophication is the process by which an excess of nutrients – mainly phosphorus and nitrogen – in a water body leads to increased growth of plant material in an aquatic body. Anthropogenic activities, such as farming, agriculture, aquaculture and industry, are the main source of nutrient input in problem areas (Jickells, 1998; Schindler, 2006; Galloway et al., 2008). Eutrophication is an issue particularly in coastal regions and areas with restricted water flow, such as lakes and rivers (Howarth and Marino, 2006; Smith, 2003). The impact of eutrophication on aquatic ecosystems is well known: nutrient availability boosts plant growth – particularly algal blooms – resulting in a decrease in water quality (Anderson et al., 2002; Howarth et al.; 2000). This can, in turn, cause death by hypoxia of aquatic organisms (Breitburg et al., 2018), ultimately driving changes in community composition (Van Meerssche et al., 2019). Eutrophication has also been linked to changes in the pH (Cai et al., 2011, Wallace et al. 2014) and depletion of inorganic carbon in the aquatic environment (Balmer and Downing, 2011). Oligotrophication is the opposite of eutrophication, where reduction in some limiting resource leads to a decrease in photosynthesis by aquatic plants, reducing the capacity of the ecosystem to sustain the higher organisms in it. Eutrophication is one of the more long-lasting water quality problems in Europe (OSPAR ICG-EUT, 2017), and is on the forefront of most European Directives on water-protection. Efforts to reduce anthropogenically-induced pollution resulted in the implementation of the Water Framework Directive (WFD) in 2000. '''CMEMS KEY FINDINGS''' The coastal and shelf waters, especially between 30 and 400N that showed active oligotrophication flags for 2020 have reduced in 2021 and a reversal to eutrophic flags can be seen in places. Again, the eutrophication index is positive only for a small number of coastal locations just north of 40oN in 2021, however south of 40oN there has been a significant increase in eutrophic flags, particularly around the Azores. In general, the 2021 indicator map showed an increase in oligotrophic areas in the Northern Atlantic and an increase in eutrophic areas in the Southern Atlantic. The Third Integrated Report on the Eutrophication Status of the OSPAR Maritime Area (OSPAR ICG-EUT, 2017) reported an improvement from 2008 to 2017 in eutrophication status across offshore and outer coastal waters of the Greater North Sea, with a decrease in the size of coastal problem areas in Denmark, France, Germany, Ireland, Norway and the United Kingdom. '''DOI (product):''' https://doi.org/10.48670/moi-00195

This visualization product displays plastic bags density per trawl. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Since the beginning of 2018, data of seafloor litter collected by international fish-trawl surveys 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 (OSPAR and MEDITS protocols) and reference lists used on a European scale. Moreover, within the same protocol, different gear types are deployed during fishing bottom trawl surveys. In cases where the wingspread and/or number of items were unknown, data could not be used because these fields are needed to calculate the density. Data collected before 2011 are affected by this filter. When the distance reported in the data was null, it was calculated from: - the ground speed and the haul duration using this formula: Distance (km) = Haul duration (h) * Ground speed (km/h); - the trawl coordinates if the ground speed and the haul duration were not filled in. The swept area is calculated from the wingspread (which depends on the fishing gear type) and the distance trawled: Swept area (km²) = Distance (km) * Wingspread (km) Densities have been calculated on each trawl and year using the following computation: Density of plastic bags (number of items per km²) = ∑Number of plastic bags related items / Swept area (km²) Percentiles 50, 75, 95 & 99 have been calculated taking into account data for all years. The list of selected items for this product is attached to this metadata. Information on data processing and calculation is detailed in the attached methodology document. Warning: the absence of data on the map doesn't necessarily mean that they don't exist, but that no information has been entered in the Marine Litter Database for this area.

1) Demographic traits These data are published data of age-specific mortality rates, age-specific lengths or weights, length and age at maturity, fecundity-length relationships, and egg size for 84 populations from 49 species of primarily commercial teleost fishes. The populations included are those for which all the life history traits under study have been estimated over a period shorter than 10 years. Traits were estimated from within the ten year window or averaged across it when data were available. Only studies in which reference population, sample size, techniques used for ageing fish and counting eggs, and models used for estimating mortality were reported are included. When only a size or age range was available, the midpoint between the extreme values was used. Raw data were converted into seven demographic traits: - Time-to-5%-survival (T.05): the time elapsed from sexual maturity until 95% of a cohort is dead. T.05 fwas estimated from an exponential mortality model, based on total mortality coefficients estimated by Virtual Population Analysis (age-structured model) in most cases or cohort analysis or catch curves. - Length-at-5%-survival (L.05). In fishes, adult size is difficult to measure because of their indeterminate growth. Adult size reported here is length at time-to-5%-survival. - Age at sexual maturity (Tm): median age at maturity was estimated directly from the data or by fitting a logistic curve to age-specific proportion mature data. When only an age range was available, the midpoint between minimum and maximum is reported. - Length at sexual maturity (Lm): median length at maturity was estimated as age at maturity. - Slope of the fecundity-length relationship (Fb): fish fecundity, defined as the number of eggs present in the ovaries immediately before spawning, is known to increase intraspecifically with the size of females. This increase is usually described by a power-law F = aLb. The exponent of this relationship, b (slope of the log-log fecundity-length regression), accounts for the increase in fecundity with size. - Fecundity at maturity (Fm): fecundity in the year of maturity was estimated from length at maturity, the fecundity-length relationship and the number of spawning bouts per year for batch spawners. - Egg volume (Egg): When information on egg size was unavailable in specific papers, values were borrowed from other studies, using the following criteria in the descending order: from the same period, the same population, the same species. In five species of Perciformes no estimate was available for any population, thus egg volume was estimated from other species of the same family. 2) Fishing pressure Three types of environments with low, moderate and high fishing pressure were defined. - To scale the pressure exerted by fishing to the natural population turn-over, it was expressed as the ratio of fishing mortality to natural mortality rates (F/M). Data were gathered from the literature together with demographic traits. Authors use the following methods to estimate natural mortality rates: intercept of a regression of total mortality on fishing effort, linear relationship known between estimates of natural mortality, growth parameters and the temperature, or multispecies models. Fishing mortality rates were estimated from Virtual Population Analysis or cohort analysis, or as the difference between total and natural mortality. Three levels of fishing pressure were defined: low fishing pressure (fishing mortality lower than natural mortality, F/M < 1), intermediate (1 <= F/M < 2) and high (F/M >= 2).

Key physico-chemical parameters (salinity, temperature, turbidity and dissolved oxygen) were measured in surface water during longitudinal transects in the Loire and Gironde estuaries in summers 2017 and 2018. This objective of this work was to determine the distribution of the dissolved oxygen and to detect potential severe desoxygenation. The transects were scheduled in order to begin the measurements at high tide from a site located upstream of an area where severe deoxygenation have been already been reported. Then, the transect was realised by sailing at low speed downstream with a multiparameter probe SAMBAT, maintained at 0.5 m below the surface, that collected a measurement every 2 minutes.

Moving 6-year analysis of Oxygen at Atlantic Sea 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 centered average of each season. 6-year periods span from 1960-1965 until 2015-2020. Observational data span from 1960 to 2020. Depth range (IODE standard depths): -3000.0, -2500.0, -2000.0, -1750, -1500.0, -1400.0, -1300.0, -1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -20.0, -10.0, -5.0, -0.0 Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5 DIVA settings. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Signal to noise ratio was fixed to 1. Background field: the data mean value is subtracted from the data. Detrending of data: no, Advection constraint applied: no. Units: umol/l

Moving 6-year analysis of Phosphate at Atlantic Sea 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 1967-1972 until 2015-2020 (winter), - from 1960-1965 until 2015-2020 (spring), - from 1968-1973 until 2015-2020 (summer), - from 1961-1966 until 2015-2020 (autumn). Observational data span from 1960 to 2020. Depth range (IODE standard depths): -2000.0, -1750, -1500.0, -1400.0, -1300.0, -1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -20.0, -10.0, -5.0, -0.0 Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5 DIVA settings. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Signal to noise ratio was fixed to 1. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no, Advection constraint applied: no. Units: umol/l

Moving 6-year analysis of Water_body_silicate in the Mediterranean Sea 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 centered average of the season. 6-years periods span from 1965-1970 until 2014-2019. Observational data span from 1960 to 2019. Depth range (IODE standard depths): 1500.0, 1400.0, 1300.0, 1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0, -30.0, -20.0, -10.0, -5.0, -0.0. Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. Profiles were interpolated at standard depths using weighted parabolic interpolation algorithm (Reiniger and Ross, 1968). GEBCO 1min topography is used for the contouring preparation. Analysed filed masked using relative error threshold 0.3 and 0.5. DIVA settings: A constant value for signal-to-noise ratio was used equal to 1. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. 'log(data)-exp(analysis' transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no. Advection constraint applied: no. Units: umol/l. The entire set of related maps can be found in the viewing service: http://ec.oceanbrowser.net/emodnet/ . Originators of Italian data sets-List of contributors: - Brunetti Fabio (OGS) - Cardin Vanessa, Bensi Manuel doi:10.6092/36728450-4296-4e6a-967d-d5b6da55f306 - Cardin Vanessa, Bensi Manuel, Ursella Laura, Siena Giuseppe doi:10.6092/f8e6d18e-f877-4aa5-a983-a03b06ccb987 - Cataletto Bruno (OGS) - Cinzia Comici Cinzia (OGS) - Civitarese Giuseppe (OGS) - DeVittor Cinzia (OGS) - Giani Michele (OGS) - Kovacevic Vedrana (OGS) - Mosetti Renzo (OGS) - Solidoro C.,Beran A.,Cataletto B.,Celussi M.,Cibic T.,Comici C.,Del Negro P.,De Vittor C.,Minocci M.,Monti M.,Fabbro C.,Falconi C.,Franzo A.,Libralato S.,Lipizer M.,Negussanti J.S.,Russel H.,Valli G., doi:10.6092/e5518899-b914-43b0-8139-023718aa63f5 - Celio Massimo (ARPA FVG) - Malaguti Antonella (ENEA) - Fonda Umani Serena (UNITS) - Bignami Francesco (ISAC/CNR) - Boldrini Alfredo (ISMAR/CNR) - Marini Mauro (ISMAR/CNR) - Miserocchi Stefano (ISMAR/CNR) - Zaccone Renata (IAMC/CNR) - Lavezza, R., Dubroca, L. F. C., Ludicone, D., Kress, N., Herut, B., Civitarese, G., Cruzado, A., Lefèvre, D., Souvermezoglou, E., Yilmaz, A., Tugrul, S., and Ribera d'Alcala, M.: Compilation of quality controlled nutrient profiles from the Mediterranean Sea, doi:10.1594/PANGAEA.771907, 2011.

This dataset contains bin-averaged optical particle measurements from Biogeochemical Argo floats. Full description of data and methodology is contained in the manuscript submitted to Science entitled "Particle fragmentation exerts strong control on biological sequestration of CO2 by the oceans". Optical measurements are particulate backscattering and chlorophyll fluorescence, and each have been partitioned into large (>100 µm) and small (<100 µm) size classes for use in estimating the rate of fragmentation of large, sinking partiles. The data cover 34 high-latitude open-ocean mesopelagic sinking particle plumes in the supolar North Atlantic and the Southern Ocean observed at 1 m vertical resolution and 2-5 day temporal resolution by profiling floats over a 30-day period. Each 30-day period is divided into five temporal bins of six days each. Vertical binning is at 50-meter intervals from the 250-950 m.