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)

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

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.

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.

The network was initiated by IFREMER from 1993 to 2009 (under the acronym REMORA) to study the rearing performance of the Pacific oyster Crassostrea gigas at a national scale. To do so, the network monitored annually the mortality and growth of standardized batches of 18-month-old oysters. Starting in 1995, the monitoring of the rearing performance of 6-month-old oyster spat was integrated into this network. These sentinel batches were distributed simultaneously each year on 43 sites and were monitored quarterly. These sites were distributed over the main French oyster farming areas and allowed a national coverage of the multiannual evolution of oyster farming performances. Most of the sites were located on the foreshore at comparable levels of immersion. Field studies were carried out by the "Laboratoires Environnement Ressources" (LER) for the sites included in their geographical area of investigation. Following the increase in spat mortality in 2008, the network evolved in 2009 (under the acronym RESCO). From this date, the network selected 13 sites among the 43 sites previously monitored in order to increase the frequency of visits (twice a month) and the number of sentinel batches. More precisely, sentinel batches of oysters corresponding to different origins (wild or hatchery, diploid or triploid) and to two rearing age classes (spat or 18-month-old adults) were selected. The monitoring of environmental variables (temperature, salinity) associated with the 13 sites was also implemented. The actions of the network have thus contributed to disentangle the biotic and abiotic parameters involved in mortality phenomena, taking into account the different compartments (environment / host / infectious agents) likely to interact with the evolution of oyster rearing performance. Finally, since 2015, the network has merged the RESCO and VELYGER networks to adopt the acronym ECOSCOPA. The general objective of this current network is to analyze the causes of spatio-temporal variability of the main life traits (Larval stage - Recruitment - Reproduction - Growth - Survival - Cytogenetic abnormalities) of the cupped oyster in France and to follow their evolution on the long term in the context of climate change. To do this, the network proposes a regular spatio-temporal monitoring of the major proxies of the life cycle of the oyster, organized in three major thematic groups: (1) proxies related to growth, physiological tolerance and survival of experimental sentinel populations over 3 age classes: (2) proxies related to reproduction, larval phase and recruitment of the species throughout its natural range in France, and: (3) proxies related to environmental parameters essential to the species (weather conditions, temperature, salinity, pH, turbidity, chlorophyll a and phytoplankton) at daily or sub-hourly frequencies. Working in a geographical network associating several laboratories, ECOSCOPA provide these monitoring within 8 sites selected among the previous ones to ensure the continuity of the data acquisition. Today, these 8 sites are considered as ecosystems of common interest, contrasted, namely : - The Thau lagoon - The Arcachon basin - The Marennes Oléron basin - The Bourgneuf Bay - The bay of Vilaine - The bay of Brest - The bay of Mont Saint Michel - The bay of Veys The ECOSCOPA network is therefore one of the relevant monitoring tools on a national scale, allowing to objectively measure through different proxies the general state of health of cultivated and wild oyster populations, and this for the different sensitive phases of their life cycle. This network aims at allowing a better evaluation, on the long term, of the biological risks incurred by the sector but also by the ecosystems, in particular under the increasing constraint of climatic and anthropic changes. Figure : Sites monitored by the ECOSCOPA network  

The ICES Working Group on Fisheries Benthic Impact and Trade-offs (WGFBIT) has developed an assessment framework based on the life history trait longevity, to evaluate the benthic impact of fisheries at the regional scale. In order to apply this framework to the Mediterranean sea, several Mediterranean longevity databases were merged together with existing North-East Atlantic ones to develop a common database. Longevity was fuzzy coded into four longevity classes: <1, 1-3, 3-10 and >10 years. Both benthic mega and macrofauna organisms are included in this dataset. Further details about both the purpose and the methodology may be found in ICES (2022) and Cuyvers et al. (2023). The result of the final dataset merging is one dataset containing the fuzzy coded average longevity (and standard deviation) for 2264 taxa and for each, the number of databases used. 

This daily High-Resolution (HR) Level 3 gridded wind product is derived from Copernicus Sentinel-1 SAR (Synthetic Aperture Radar) observations, over the Mediterranean Sea ("MED" area). It is based on the European Space Agency (ESA) Level-2 OCN products at the highest available resolution. Although L2-OCN products already contain wind vectors, those are calculated using the CMOD5.n Geophysical Model Function (GMF) applied to the co-polarized (co-pol) VV channel (emitting in Vertical polarization and receiving in Vertical polarization). This VV GMF was mapped from scatterometer sensors (Hersbach et al., 2007) which are only able to use co-pol measurements. However, these co-pol GMF are known to lose sensitivity for wind above 20 m/s. Therefore, wind based on such GMF alone, are known to under-estimate wind speed (Polverari et al., 2022). For the L3 products winds based on SAR, we take advantage of the available cross-polarized (cross-pol) VH channel (emitting in Vertical polarization and receiving in Horizontal polarization) for which GMF were specifically derived based on C-Band SAR (Mouche et al., 2017, Mouche et al., 2019). Winds estimated from the combination of both the co-pol and cross-pol channels are referred to as dual-polarization (or dual-pol) winds. As shown in Mouche et al. (2019), taking advantage of the dual polarization strongly improves the wind estimation for high wind conditions thanks to the much greater VH channel sensitivity compared to VV. These new wind estimations are then gridded with a 0.012 degree resolution (between 0.5 and 1.2 km in zonal direction depending on the latitude and 1.3 km in meridional direction) using a cylindrical equidistant projection, independently for ascending and descending satellite passes and for each satellite (so 4 wind fields are available per day for two satellites). This dataset is generated over all Sentinel-1 mission time series starting from March 2018 and updated in delayed mode with a 4-months delay. It is also produced for 4 other different European areas. This dataset is produced and disseminated in the frame of Copernicus Marine Service.

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 )