As part of the marine water quality monitoring of the “Pertuis” and the “baie de l’Aiguillon” (France), commissioned by the OFB and carried out by setec énergie environnement, three monitoring stations were installed. Two of them were set up at the mouths of the Charente and Seudre rivers on February 6 and 27, 2019, respectively, while a third was deployed in the Bay of Aiguillon on March 24, 2021. The dataset presented here concerns the station installed in the Charente estuary. Measurements are organized into .csv files, with one file per year. Data is collected using a SAMBAT multiparameter probe, which records the following parameters: - Temperature (-5 to 35 °C) - Conductivity (0 to 10 mS/cm) - Pressure (0 to 10 m) - Turbidity (0 to 300 NTU) - Dissolved Oxygen (0 to 20 mg/L & 0 to 200 %) - Fluorescence (0 to 50 µg/l) - PH (0/14)

Seasonal climatology of Water body dissolved inorganic nitrogen for Loire river for the period 1971-2024 and for the following seasons: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Observation data span from 1971 to 2024. 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.12, 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): 72000 (winter), 43000 (spring), 215000 (summer), 150000 (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.loglin) was applied to the data prior to the analysis. Background field: the vertically-filtered data mean profile is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

Moving 6-year analysis of Water body phosphate 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 1970-1975 until 2018-2023. Description of DIVA analysis: The computation was done with the DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 30sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps grids is dx=dy=0.125 degrees (around 13.5km and 10.9km accordingly). The vertical resolution is 20 depth levels: [0.,5.,10.,20.,30.,50.,75.,100.,125.,150.,200.,250.,300.,400.,500.,600.,700.,800.,900.,1000.]. The horizontal correlation length is 200km. The vertical correlation length (in meters) was set twices the vertical resolution: [10.,10.,20.,20.,40.,50.,50.,50.,50.,100.,100.,100.,200.,200.,200.,200.,200.,200.,200.,200.]. Duplicates check was performed using the following criteria for space and time: dlon=0.001deg., dlat=0.001deg., ddepth=1m, dtime=1hour, dvalue=0.1. The error variance (epsilon2) was set equal to 1 for profiles and 10 for time series to reduce the influence of close data near the coasts. An anamorphosis transformation was applied to the data (function DIVAnd.Anam.loglin) to avoid unrealistic negative values: threshold value=200. A background analysis field was used for all years (1970-2023) with correlation length equal to 600km and error variance (epsilon2) equal to 20. Quality control of the observations was applied using the interpolated field (QCMETHOD=3). Residuals (differences between the observations and the analysis (interpolated linearly to the location of the observations) were calculated. Observations with residuals outside the minimum and maximum values of the 99% quantile were discarded from the analysis. 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.

The present repository makes available the model, material and outputs of the ISIS-Fish modeling work showcased in the peer-reviewed scientific article by Bastardie et al. 2025. As part of the SEAwise research project (seawiseproject.org), we used an ISIS-Fish database (Mahevas et al 2003, Pelletier et al. 2009, isis-fish.org) previously developed within the MACCO project which describes the mixed demersal fishery in the Bay of Biscay. For this application, the spatial extent of the fishery is the Bay of Biscay, defined here by ICES divisions 8a, 8b and 8d and the resolution chosen is 1/16 ICES statistical rectangle. The biological module (Vajas et al. 2024) includes 7 species of economic interest in the mixed demersal fishery: European hake (Merluccius merluccius), common sole (Solea solea), Norway lobster (Nephrops norvegicus), megrim (Lepidorhombus whiffiagonis), anglerfish (Lophius piscatorius) and two ray species (Raja clavata, Leucoraja naevus). The fishing activities module (Mahevas et al. 2024) is made up of 41 demersal fleets (including all French vessels < 12 meters and > 12 meters fishing in this area, Spanich, UK and Belgium fleets) and 431 métiers (combination of a gear, location and mix of target species) catching these 7 species, as target or bycatch. Monthly effort of a fleet distributes among the possible métiers (those historically practiced). The biological and fishing activity modules are identical to the published version. The original model used here has been calibrated on historical catch data 2015-2018 by tuning accessibility and catchability parameters. In the present application the Bay of Biscay model is used to investigate the spatial- and effort- based fisheries management strategies. Consistently with for a task of the SEAwise project (Bastardie et al. 2024) simulations were conducted from 2021 onwards, projecting the effect of an implementation of 3 different closures from 2022 to 2050, under current fishing effort conditions or in a context of fishing effort reduction. Outcomes of these simulations are averaged over short/medium (10 year horizon) and long-term period (20 year horizon). The data project includes: 1) the database including the biological module and fishing activity module; 2) 8 .properties files, each corresponding to one combination of management measure and closure, to restore the simulations parameters in the ISIS-Fish interface and reproduce the simulation runs; 3) the .java scripts to force effort dynamics and simulate spatio-temporal closures, as well as generate the main output files - they will be called by the ISIS-Fish software once the simulations restored 4) the .rds containing the main outputs of the simulations and the associated .html document displaying the R code to compute the indices of interest at different levels of aggregation and reproduce the figures in Bastardie et al. 2025. All files are provided in the Zip. Associated with this material, a study summary and a readme .docx are provided. The first one provides context on the present work and describes the model and simulations' design. The second provides guidelines to reproduce the simulations and their derived outcomes from the data project material made available in this repository. They are both directly downloadable from this repository and are also copied to the zipped folder containing the data project. All the data are reproducible using isis-fish-4.4.8.1 (isis-fish.org; available at forge.codelutin.com) and R 4.2.0.

210Pb, 226Ra and 137Cs were measured by non-destructive gamma spectrometry on marine sediment cores, collected during RIKEAU 2002 cruise on board r/v Thalia, on the shelf of the Bay of Biscay

Moving 6-year analysis of Water body dissolved inorganic nitrogen (DIN) 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 1990-1995 until 2018-2023. Description of DIVA analysis: The computation was done with the DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 30sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps grids is dx=dy=0.125 degrees (around 13.5km and 10.9km accordingly). The vertical resolution is 21 depth levels: [0.,5.,10.,20.,30.,50.,75.,100.,125.,150.,200.,250.,300.,400.,500.,600.,700.,800.,900.,1000.,1100.]. The horizontal correlation length is 200km. The vertical correlation length (in meters) was set twices the vertical resolution: [10.,10.,20.,20.,40.,50.,50.,50.,50.,100.,100.,100.,200.,200.,200.,200.,200.,200.,200.,200.,200.]. Duplicates check was performed using the following criteria for space and time: dlon=0.001deg., dlat=0.001deg., ddepth=1m, dtime=1hour, dvalue=0.1. The error variance (epsilon2) was set equal to 1 for profiles and 10 for time series to reduce the influence of close data near the coasts. An anamorphosis transformation was applied to the data (function DIVAnd.Anam.loglin) to avoid unrealistic negative values: threshold value=200. A background analysis field was used for all years (1990-2023) with correlation length equal to 600km and error variance (epsilon2) equal to 20. Quality control of the observations was applied using the interpolated field (QCMETHOD=3). Residuals (differences between the observations and the analysis (interpolated linearly to the location of the observations) were calculated. Observations with residuals outside the minimum and maximum values of the 99% quantile were discarded from the analysis. 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 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 bottom trawl surveys. In cases where the wingspread and/or number of items were/was unknown, it was not possible to use the data because these fields are needed to calculate the density. Data collected before 2011 are concerned by this filter. When the distance reported in the data was null, it was calculated from: - the ground speed and the haul duration using the following 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 was calculated from the wingspread (which depends on the fishing gear type) and the distance trawled: Swept area (km²) = Distance (km) * Wingspread (km) Densities were 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 was 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 were calculated taking into account data for all years. More information on data processing and calculation are detailed in the attached methodology document. 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. 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.

This visualization product displays the size of litter in percent per net per year 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 research and monitoring protocols as MSFD monitoring. To calculate percentages for each size, formula applied is: Size (%) = (∑number of particles of each size)*100 / (∑number of particles of all size) When the number of microlitters was not filled or was equal to zero, it was not possible to calculate the percentage. Standard vocabularies for microliter sizes are taken from Seadatanet's H03 library (https://vocab.seadatanet.org/v_bodc_vocab_v2/search.asp?lib=H03 ). Different protocols with different degrees of precision were used to classify the sampled micro-litters. Consequently, on the map, the distribution of micro-litter in the size classes depends on the protocol applied during the survey. 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.

In October 2019 we chose 15 sites from the 2019 EVHOE survey for environmental DNA (eDNA) sampling. The French international EVHOE bottom trawl survey is carried out annually during autumn in the BoB to monitor demersal fish resources. At each site, we sampled seawater using Niskin bottles deployed with a circular rosette. There were nine bottles on the rosette, each of them able to hold ∼5 l of water. At each site, we first cleaned the circular rosette and bottles with freshwater, then lowered the rosette (with bottles open) to 5 m above the sea bottom, and finally closed the bottles remotely from the boat. The 45 l of sampled water was transferred to four disposable and sterilized plastic bags of 11.25 l each to perform the filtration on-board in a laboratory dedicated to the processing of eDNA samples. To speed up the filtration process, we used two identical filtration devices, each composed of an Athena® peristaltic pump (Proactive Environmental Products LLC, Bradenton, Florida, USA; nominal flow of 1.0 l min–1 ), a VigiDNA 0.20 μm filtration capsule (SPYGEN, le Bourget du Lac, France), and disposable sterile tubing. Each filtration device filtered the water contained in two plastic bags (22.5 l), which represent two replicates per sampling site. We followed a rigorous protocol to avoid contamination during fieldwork, using disposable gloves and single-use filtration equipment and plastic bags to process each water sample. At the end of each filtration, we emptied the water inside the capsule that we replaced by 80 ml of CL1 conservation buffer and stored the samples at room temperature following the specifications of the manufacturer (SPYGEN, Le Bourget du Lac, France). We processed the eDNA capsules at SPYGEN, following the protocol proposed by Polanco-Fernández et al., (2020). Half of the extracted DNA was processed by Sinsoma using newly developped ddPCR assays for European seabass (Dicentrachus labrax), European hake (Merluccius merluccius) and blackspot seabream (Pagellus bogaraveo).  The other half of the extracted DNA was analysed using metabarcoding with teleo primer. The raw metabarcoding data set is available at https://www.doi.org/10.16904/envidat.442 Bottom trawling using a GOV trawl was carried out before or after water sampling. The catch was sorted by species and catches in numbers and weight were recorded. No blackspot seabream individuals were caught.   Data content: * ddPCR/: contains the ddPCR counts and DNA concentrations for each sample and species. * SampleInfo/: contains the filter volume for each eDNA sample. * StationInfo/: contains metadata related to the data collected in the field for each filter. * Metabarcoding/: contains metabarcoding results for teleoprimer. * Trawldata/: contains catch data in numbers and weight (kg).      

Moving 6-year analysis of Water body dissolved oxygen concentration 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 1970-1975 until 2019-2024. Description of DIVA analysis: The computation was done with the DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 30sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps grids is dx=dy=0.125 degrees (around 13.5km and 10.9km accordingly). The vertical resolution is 27 depth levels: [0.,5.,10.,20.,30.,50.,75.,100.,125.,150.,200.,250.,300.,400.,500.,600.,700.,800.,900.,1000.,1100.,1200.,1300.,1400.,1500.,1750.,2000.]. The horizontal correlation length is 200km. The vertical correlation length (in meters) was set twices the vertical resolution: [10.,10.,20.,20.,40.,50.,50.,50.,50.,100.,100.,100.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,500.,500.,500.]. Duplicates check was performed using the following criteria for space and time: dlon=0.001deg., dlat=0.001deg., ddepth=1m, dtime=1hour, dvalue=0.1. The error variance (epsilon2) was set equal to 1 for profiles and 10 for time series to reduce the influence of close data near the coasts. An anamorphosis transformation was applied to the data (function DIVAnd.Anam.loglin) to avoid unrealistic negative values: threshold value=200. A background analysis field was used for all years (1970-2024) with correlation length equal to 600km and error variance (epsilon2) equal to 20. Quality control of the observations was applied using the interpolated field (QCMETHOD=3). Residuals (differences between the observations and the analysis (interpolated linearly to the location of the observations) were calculated. Observations with residuals outside the minimum and maximum values of the 99% quantile were discarded from the analysis. 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.