The data sets presented here result from the long-term monitoring of individual growth patterns in anchovy and sardine in the Bay of Biscay, from 2000 to 2018. They derived from the PELGAS survey series (Doray et al., 2018), which monitors annually the Bay of Biscay pelagic ecosystem since 2000. The survey is performed in May during the peak spawning of anchovy and main spawning of sardine. Among the many data collected, anchovy and sardine populations are assessed by combining acoustic records with pelagic trawl hauls catches and ICES survey protocoles are used, as detailed in Doray et al. (2021). Briefly, fish acoustic backscatter are recorded along survey transect lines and pelagic trawl hauls undertaken opportunistically to identify echotraces to species and collect fish samples for acquiring biometric data. At each trawl haul and for each species, a random subsample of individuals is taken to establish the species’ length distributions. For anchovy and sardine, this subsample is further subsampled, spanning the whole length range, to take individual fish measurements. These amount to extracting otoliths and measuring individuals’ age, length, weight, sexual maturity and other parameters. Individual measurements are taken on fourty individuals of anchovy and sardine when the species are present in the catch. For each individual fish, the two otolith sagittae are extracted on board, mounted in leukit for age reading on board when permitting and/or on land in the laboratory. Growth patterns in the otoliths were analysed on land with a binocular stereomicroscope under reflected natural light. For anchovy, otoliths’ growth was measured for all individuals in all the hauls. For sardine, trawl hauls were selected and all individual otoliths were measured in each selected haul. The selection was made using the geographical stratification defined in Petitgas et al. (2018) based on the ecosystem spatial structure. An average of two to three hauls in each of the four strata were selected per year. The otoliths mounted in leukit were imaged and growth-at-age in the otoliths was measured with the software TNPC (Traitement numérique des pièces calcifiées: Mahé et al., 2009). Under the binocular microscope and natural light, the otoliths showed hyaline (aragonite-poor) rings corresponding to winter periods of low growth and between the rings, white opaque (aragonite-rich) portions corresponding to annual growth periods. The annual ring determination, the age assignment and the measurement of annual ring diameters followed ICES protocoles and guidelines for anchovy and sardine (ICES, 2010; 2011). The age was estimated as the number of hyaline rings. If the edge was hyaline, it was counted as a ring as a hyaline edge observed within the first half of the year is assumed to represent the last winter. The diameter of each annual ring was measured from middle of the hyaline ring on one side to the middle of the ring on the opposite side along the major elongated axis of the otolith and passing through its centre. The distance was expressed in mm after calibration of the stereomiscroscope and the pixel images. Such diameter corresponded to growth-at-age. Ages 0 to 4 were considered (diameters R1 to R5). The total diameter of the otolith was also measured. The data sets span 19 years, from 2000 to 2018 and comprise 20,186 and 8,624 individual fish analyzed at 535 and 235 trawl hauls for anchovy and sardine, respectively. These data sets were used by Boëns et al. (2021 and 2023) to analyse changes in growth patterns and growth-selective mortality at age in anchovy and sardine under environmental and fishing pressures. References: Doray, M., Boyra, G. and Van Der Kooij, J. (eds) (2021). ICES Survey Protocols – Manual for acoustic surveys coordinated under ICES Working Group on Acoustic and Egg Surveys for Small Pelagic Fish (WGACEGG). 1st Edition. ICES Techniques in Marine Environmental Sciences, 64. https://doi.org/10.17895/ices.pub.7462  Doray, M., Petitgas, P., Romagnan, J.-B., Huret, M., Duhamel, E., Dupuy, Ch., Spitz, J., Authier, M., Sanchez, F., Berger, L., Doremus, G., Bourriau, P., Grellier, P. and Masse, J. (2018). The PELGAS survey: ship-based integrated monitoring of the Bay of Biscay pelagic ecosystem. Progress In Oceanography, 166, 15-29. https://doi.org/10.1016/j.pocean.2017.09.015 ICES (2010). Report of the Workshop on Age reading of European anchovy (WKARA), 9-13 November 2009, Sicily, Italy. ICES CM 2009/ACOM: 43. 122 pp. https://doi.org/10.17895/ices.pub.19280525 ICES (2011). Report of the Workshop on Age Reading of European Atlantic Sardine (WKARAS), 14-18 February 2011, Lisbon, Portugal. ICES CM 2011/ACOM:42. 91 pp. https://doi.org/10.17895/ices.pub.19280855 Petitgas, P., Huret, M., Dupuy, Ch., Spitz, J., Authier, M., Romagnan, J.-B. and Doray, M. (2018). Ecosystem spatial structure revealed by integrated survey data. Progress In Oceanography, 166, 189-198. https://doi.org/10.1016/j.pocean.2017.09.012 Mahe, K., Bellail, R., Dufour, J.-L., Boiron-Leroy, A., Dimeet, J., Duhamel, E., Elleboode, R., Felix, J., Grellier, P., Huet, J., Labastie, J., Le Roy, D., Lizaud, O., Manten, M.-L., Martin, S., Metral, L., Nedelec, D., Verin, Y. and Badts, V. (2009). Synthèse française des procédures d'estimation d'âge / French summary of age estimation procedures. https://archimer.ifremer.fr/doc/00000/7294/ Boëns, A., Grellier, P., Lebigre, Ch. and Petitgas, P. (2021). Determinants of growth and selective mortality in anchovy and sardine in the Bay of Biscay. Fisheries Research, 239, 105947. https://doi.org/10.1016/j.fishres.2021.105947 Boëns, A., Ernande, B., Petitgas, P. and Lebigre, Ch. (2023). Different mechanisms underpin the decline in growth of anchovies and sardines of the Bay of Biscay. Evolutionary Applications, 16: 1393–1411. https://doi.org/10.1111/eva.13564  

This product displays for DDT, DDE, and DDD, median values of the last 6 available years that have been measured per matrix and are present in EMODnet regional contaminants aggregated datasets, v2022. The median values ranges are derived from the following percentiles: 0-25%, 25-75%, 75-90%, >90%. Only "good data" are used, namely data with Quality Flag=1, 2, 6, Q (SeaDataNet Quality Flag schema). For water, only surface values are used (0-15 m), for sediment and biota data at all depths are used.

This visualization product displays the location of all the surveys present in the EMODnet marine litter database (MLDB). The different fishing gears used are represented by different colors. 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. Unlike other EMODnet seafloor litter products, all trawls surveyed since 2006 are included in this map even if the wingspread and/or the distance are unknown. Only surveys with an unknown number of items were excluded from this product. More information on data processing and calculation are detailed in the document attached. 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.

Moving 6-year analysis of Water body phosphate 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 1950/1955 until 2016/2021. Observation data span from 1950 to 2021. 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. Descrption of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.4, 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 250km 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.Anam.loglin) was applied to the data prior to the analysis to avoid unrealistic negative values. 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.

This visualization product displays the spatial distribution of seafloor litter 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 (number of items per km²) = ∑Number of items / Swept area (km²) Then a grid with 30km x 30km cells is used to calculate the weighted mean of densities in each cell from the formula : Weighted mean (number of items per km²) = ∑ (Distance (km) * Density (number of items per km²)) / ∑ Distance (km) Percentiles 50, 75, 95 & 99 have been calculated taking into account data for all years. More information on data processing and calculation are detailed in the document attached. 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. This work is based on the work presented in the following scientific article: O. Gerigny, M. Brun, M.C. Fabri, C. Tomasino, M. Le Moigne, A. Jadaud, F. Galgani, Seafloor litter from the continental shelf and canyons in French Mediterranean Water: Distribution, typologies and trends, Marine Pollution Bulletin, Volume 146, 2019, Pages 653-666, ISSN 0025-326X, https://doi.org/10.1016/j.marpolbul.2019.07.030.

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 eutrophication and acidity, and covers the Mediterranean Sea. Data were aggregated and quality controlled by the 'Hellenic Centre for Marine Research, Hellenic National Oceanographic Data Centre (HCMR/HNODC)' in Greece. ITS-90 water temperature and water body salinity variables have also been included ('as are') to complete the eutrophication and acidity data. If you use these variables for calculations, please refer to SeaDataNet for the quality flags: https://www.seadatanet.org/Products/Aggregated-datasets. Regional datasets concerning eutrophication and acidity are automatically harvested, and the resulting collections are aggregated and quality controlled using ODV Software and following a common methodology for all sea regions (https://doi.org/10.13120/8xm0-5m67). Parameter names are based on P35 vocabulary, which relates to EMODnet Chemistry aggregated parameter names and is available at: https://vocab.nerc.ac.uk/search_nvs/P35/. When not present in original data, water body nitrate plus nitrite was calculated by summing all nitrate and nitrite parameters. The same procedure was applied for water body dissolved inorganic nitrogen (DIN), which was calculated by summing all nitrate, nitrite, and ammonium parameters. Concentrations per unit mass were converted to a unit volume using a constant density of 1.025 kg/L. The aggregated dataset can also be downloaded as an ODV collection and spreadsheet, which is composed of a metadata header followed by tab separated values. This spreadsheet can be imported to ODV Software for visualisation (more information can be found at: https://www.seadatanet.org/Software/ODV).

Moving 6-year analysis and visualization of Water body phosphate in the North Sea. Four seasons (December-February, March-May, June-August, September-November). Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis) tool, version 2.7.9. results were subjected to the minfield option in DIVAnd to avoid negative/underestimated values in the interpolated results; error threshold masks L1 (0.3) and L2 (0.5) are included as well as the unmasked field. The depth dimension allows visualizing the gridded field at various depths.

This product displays for Benzo(a)pyrene, median values of the last 6 available years that have been measured per matrix and are present in EMODnet regional contaminants aggregated datasets, v2022. The median values ranges are derived from the following percentiles: 0-25%, 25-75%, 75-90%, >90%. Only "good data" are used, namely data with Quality Flag=1, 2, 6, Q (SeaDataNet Quality Flag schema). For water, only surface values are used (0-15 m), for sediment and biota data at all depths are used.

This visualization product displays the spatial distribution of fishing related items 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 fishing related items (number of items per km²) = ∑Number of fishing related items / Swept area (km²) Then a grid with 30km x 30km cells is used to calculate the weighted mean of densities in each cell from the formula : Weighted mean (number of items per km²) = ∑ (Distance (km) * Density (number of items per km²)) / ∑ Distance (km) Percentiles 50, 75, 95 & 99 have been calculated taking into account data for all years. More information on data processing and calculation are detailed in the document attached. 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. This work is based on the work presented in the following scientific article: O. Gerigny, M. Brun, M.C. Fabri, C. Tomasino, M. Le Moigne, A. Jadaud, F. Galgani, Seafloor litter from the continental shelf and canyons in French Mediterranean Water: Distribution, typologies and trends, Marine Pollution Bulletin, Volume 146, 2019, Pages 653-666, ISSN 0025-326X, https://doi.org/10.1016/j.marpolbul.2019.07.030.

Seasonal climatology of Water body chlorophyll-a for Loire river for the period 1971-2021 and for the following seasons: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Observation data span from 1971 to 2021. Depth levels (m): [0.0, 2.0, 4.0, 6.0, 8.0, 10.0, 15.0, 20.0, 25.0, 30.0, 35.0, 40.0, 45.0, 50.0, 60.0, 70.0, 80.0, 90.0, 100.0, 110.0, 120.0, 130.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.4, using GEBCO 15 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.01 degrees. Horizontal correlation length is defined seasonally (in meters): 14000 (winter), 52000 (spring), 42000 (summer), 125000 (autumn). Vertical correlation length was optimized and vertically filtered and a seasonally-averaged profile was used (DIVAnd.fitvertlen). 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.Anam.loglin) was applied to the data prior to the analysis to avoid unrealistic negative values. Background field: the vertically-filtered data mean profile is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: mg/m3.