A consistent dataset of bottom trawl survey data spanning 47 years in the Bay of Biscay was assembled. The dataset includes data from the current EVHOE survey from 1987 to 2019 and two previous surveys carried out in 1973 and 1976. The recent EVHOE time-series from 1997 is also available from DATRAS (https://www.ices.dk/data/data-portals/Pages/DATRAS.aspx). The catch in numbers and weight (kg) per haul of all Rajiformes species caught in these surveys is provided. Haul information is provided for all hauls, including those with no catch of Rajiformes. Areas of the sampling strata of the survey and spatial polygones of these strata are provided in separate files.

The anchovy (Engraulis encrasicolus) and sardine (Sardina pilchardus) populations in the Bay of Biscay are jointly surveyed each year in May since 2000 and in September since 2003 by means of acoustic surveys. The integrated survey PELGAS (Doray et al., 2018) is run in May by France and covers the French shelf of the Bay of Biscay. Its objectives are to monitor the Bay of Biscay pelagic ecosystem in springtime and assess the biomass of its small pelagic fish species, including sardine and anchovy. Amongst many information on the ecosystem, the survey PELGAS provides knowledge on the adults of anchovy and sardine during their spawning in spring. The survey JUVENA (Boyra et al., 2013) is run in September by Spain. It has a larger spatial coverage than PELGAS, including part of the Spanish coast and open ocean outside the shelf because it targets juvenile anchovy. It also provides knowledge on sardine as well as other pelagic species. Both surveys are coordinated by the ICES Working Group on Acoustic and Egg Surveys for Small Pelagic Fish (WGACEGG), together with other pelagic surveys in ICES areas 7, 8 and 9. Survey protocoles are detailed in Doray et al. (2021). Briefly, fish backscatter data are recorded along the survey transect lines and pelagic trawl hauls are undertaken opportunistically to identify the echotraces to species and collect fish samples for biometric data. The trawl haul catches have provided the anchovy and sardine length data, from which the maps presented here are derived. At each trawl haul, the catch is sorted by species and weighted. A subsample by species is measured to estimate the species’ length distribution. The four maps presented here correspond to the average maps of anchovy and sardine length distributions in May and September, derived from the PELGAS and JUVENA trawl haul data series. The maps were obtained by kriging, following the procedure explained in Petitgas et al. (2011) for mapping functions instead of variables. For each species, the experimental length distribution at each haul was fitted by a linear combination of Legendre polynomials, the coefficients of which were co-kriged. The number of polynomials varied from 15 to 22 depending on the survey and species, with a higher number for sardine and in autumn. The length histogram at each grid node was then deduced from the mapped coefficients. When the length distribution at a given haul was estimated with less than 40 individual fish, the haul was not taken into account for mapping. This threshold defined presence and absence of the species in the haul data sets. The trawl hauls from 2000 to 2019 were pooled for the PELGAS series (1965 stations) and from 2003 to 2020 for JUVENA (852 stations). The mapping was performed on the same grid for both PELGAS and JUVENA and both species, and with similar moving kriging neighbourhoods. The grid has a mesh size of 0.25 x 0.25 decimal degree square. In addition to mapping the length distribution, presence/absence was mapped for each species by ordinary kriging on the same grid and with the same neighbourhoods as previously. The computations were performed in R (version 4.0.5) with the RGeostats package (version 13.0.1) freely available at http://rgeostats.free.fr. The map data files comprise for each species the following information: the geographical coordinates of the grid points, the probability of presence and the probability of each length class of width 0.5 cm ranging from 2.5 to 27 cm.

The Pélagiques Gascogne (PELGAS, Doray et al., 2000) integrated survey aims at assessing the biomass of small pelagic fish and monitoring and studying the dynamics and diversity of the Bay of Biscay pelagic ecosystem in springtime. PELGAS has been conducted within the EU Common Fisheries Policy Data Collection Framework and Ifremer’s Fisheries Information System. Details on survey protocols and data processing methodologies can be found in Doray et al., (2014, 2018). This dataset comprises the abundance (no. of individuals), biomass (metric tons), mean length (cm), mean weight (g) of marine organisms collected by midwater trawling to identify fish echoes detected during PELGAS surveys (2000-2018). All parameters have been raised to the trawl haul level. Trawl haul metadata and species reference list are also provided.

The bathymetric DEM of the atlantic facade with a resolution of 0.001° (~ 100 m) was prepared as part of the HOMONIM project. It encompasses part of the North Sea, the Channel and the Bay of Biscay. Offshore, this DTM extends beyond the slope to approximately 4800 m depth. The DEM is designed to be used in hydrodynamic models in order to produce high-precision forecasts for coastal water levels and sea conditions and therefore improve the pertinence of the Waves-Submersion monitoring programme. This product is available with the Lowest Astronomic Tide (LAT) or the Mean Sea Level (MSL) as a vertical datum.

This dataset comprises stomach contents of small pelagic fish species on the french shelf of the Bay of Bisacy, in spring, autumn and winter, from 2004 to 2024. The spring data were acquired in May on the pelagic survey series PELGAS from 2004 to 2024, the autumn data in October/Novermber on the demersal survey series EVHOE from 2020 to 2024 and the winter data were acquired on chartered fishing vessels in February 2023 and 2024. The dataset concerns anchovy (Engraulis encrasicolus) and sardine (Sardina pilchardus) in the 3 seasons and also mackerel (Scomber scombrus), sprat (Sprattus, sprattus) and horse mackerel (Trachurus trachurus) in spring for some years. The dataset represents a unique long-term monitoring of stomach contents characterized with a low taxonomic resolution and semi-quantitative abundance quotation.  The pelagic ecosystem survey PELGAS (Doray et al., 2018) is run in each year in May since 2000, to monitor the Bay of Biscay pelagic ecosystem at springtime and assess the biomass of its small pelagic fish species. During the survey, pelagic trawl hauls are undertaken to identify echotraces to species and to measure individual fish traits. All hauls are performed during day time. In 2010, some hauls were undertaken at night to sample stomach contents over the day/night cycle. The fish stomachs are sampled from the haul catch. For a given species, twenty individuals are selected at random from the catch, their stomachs dissected and preserved. This is repeated at three hauls in each of the ten spatial strata defined to cover the entire Biscay shelf. In some years, fish length categories (lower and greater than 14 cm for anchovy and 18 cm for sardine) were also considered when sampling the stomachs. Stomach sampling by species depended on the trawl haul catch and all species were not systematically sampled jointly at the same trawl haul. Also, the number of stations with stomach sampling varied between species and years. The stomachs were preserved in formaline until 2018 and in ethanol since. Anchovy and sardine stomach sampling on the demersal survey EVHOE (Mahe and Poulard, 2005) followed the same protocole as for PELGAS but with fewer stations, depending on the catch of anchovy and sardine in the bottom trawl. In 2020 due to the Covid pandemic, the PELGAS survey was canceled and to compensate, a pair-trawler was chartered in autumn to perform some pelagic trawl hauls during the EVHOE 2020 survey. In winter 2023 and 2024 a pair-trawler was also chartered, for identifying echotraces observed previously on the survey DRIX (Doray et al., 2024) in the area delimited by the Gironde and Loire estuaries, the coast and the 100 m isobath. On the fishing vessels the fish were frozen onboard, the stomachs were dissected on land in the laboratory and preserved in ethanol.  The taxonomic analysis of the stomach contents was performed in the laboratory under a binocular magnifyer by the company LAPHY. A simplified taxonomic resolution was used, which considered five ichtyoplankton groups, two copepod groups, euphausids or mysids, amphipods, two decapod groups, other crustacea, other zooplankton, phytoplankton and pulp. Taxon abundance was defined by a quotation : 0 (absence), 1 (presence : <10 individuals), 2 (abundant : between 10 and 100), 3 (very abundant : > 100). The dataset comprises trawl haul information, information on the quality of the stomach contents and abundance quotes for the list of plankton taxons. A preliminary analysis of the data (Petitgas, 2024) showed a large overlap in stomach contents between species, the importance of small copepods in the diets, and how different drivers such as habitat and length influence the diets. 

The Pélagiques Gascogne (PELGAS, Doray et al., 2000) integrated survey aims at assessing the biomass of small pelagic fish and monitoring and studying the dynamics and diversity of the Bay of Biscay pelagic ecosystem in springtime. PELGAS has been conducted within the EU Common Fisheries Policy Data Collection Framework and Ifremer’s Fisheries Information System. The PELGAS survey model has allowed for the establishment of a long-term time-series of spatially-explicit data of the Bay of Biscay pelagic ecosystem since the year 2000. Main sampled components of the targeted ecosystem are: hydrology, phytoplankton, mesozooplankton, fish and megafauna (cetacean and seabirds). This dataset presents gridded maps of standard pelagic ecosystem parameters collected in the main sampled components during the PELGAS survey. Ecosystem parameters were mapped on a 15km x 15km grid by applying a block averaging procedure (Petitgas et al., 2009, 2014). The dataset also includes the ecologically meaningful survey dates proposed by Huret et al. (2017), mapped on the same grid. Details on survey protocols and data processing methodologies can be found in Doray et al., (2014, 2017a). This dataset was used in Authier et al., 2017; Doray et al., 2017b, 2017c, 2017a; Huret et al., 2017; Petitgas et al., 2017.

These data are outputs of a spatio-temporal model inferring fish distribution. The maps are based on high-resolution catch data (VMS-logbook). They have a montly time resolution and a 0.05° spatial resolution. Four demersal species of the Bay of Biscay are available in the dataset: common sole (Solea solea), megrim (Lepidorhombus whiffiagonis), anglerfish (Lophius spp) and thornback ray (Raja clavata). Maps are provided for year 2008 to 2018 ; they were produced in the context of the MACCO project (https://www.macco.fr/en/accueil-english/), an Ifremer project that aims at proposing alternative management strategies for the mixed demersal fisheries of the Bay of Biscay.

We developed a panel of single nucleotide polymorphism (SNP) markers for thornback ray Raja clavata using a RADSeq protocole. Demultiplexed sequences were aligned to the genome of Leucoraja erinacea which was used as reference genome. From an initial set of 389 483 putative SNPs, 7741 SNPs with the largest minor allele frequency were selected for implementation on an Infinium® XT iSelect-96 SNP-array implemented by LABOGENA DNA. For the array, SNPs [T/C] and [T/G] were replaced by those from the complementary strand [A/G] and [A/C] respectively. For some SNPs, a second SNP was found in the 50 nucleotide bases flanking sequence. In these cases, two SNP probes were developed with each of the two alleles of the second SNP. A SNP probe naming convention was adopted to identify these pairs of probes corresponding to the same SNP locus: “MAJ” or “MIN” followed by the corresponding base was included in the probe name. For some of these pairs, only one of the two markers could be developed, resulting in a total set of 9120 SNP probes, including 6360 single SNP probes, 10 MAJ or MIN probes for which a single probe was successfully developed, and 1375 pairs of probes with MAJ and MIN versions. The 9120 SNP genotypes were then scored using the clustering algorithm implemented in the Illumina® GenomeStudio Genotyping Analysis Module v2.0.3 for 7726 individual samples, including duplicates, mostly from the Bay of Biscay but also from the Mediterranean Sea and West Iberia. Overall, 1643 SNPs failed to be genotyped in all individuals, for 319 markers the minor allele was not found and 7158 markers (including 1974 for 987 MIN-MAJ pairs) produced bi-allelic genotypes. The majority of these SNPs had a minor allele frequency between 0.1 and 0.5. The MIN-MAJ probes can be used for quality checking the genotyping results

The West Gironde Mud Patch (WGMP) is a 420-km2 mud belt in the Bay of Biscay, located 25 km off the mouth of the Gironde estuary. This clay-silt feature of 4 m in thickness extends between 30 and 75m water depth, surrounded by the sands and gravels that cover the North Aquitaine continental shelf. Interface cores were collected during JERICOBent-1 cruise (October 2016; Deflandre (2016) doi.org/10.17600/16010400) along two cross-shelf transects for a total of 9 sites. Each sediment core was carefully extruded every 0.5 cm from the top core to 4 cm and every 1 cm below until the core bottom. The sediment layers were used to determine dry bulk density, grain size and selected radioisotope activities (210Pb, 226Ra, 137Cs, 228Th, K).

The Coriolis Ocean Dataset for Reanalysis for the Ireland-Biscay-Iberia region (hereafter CORA-IBI) product is a regional dataset of in situ temperature and salinity measurements. The latest version of the product covers the period 1950-2014. The CORA-IBI observations comes from many different sources collected by Coriolis data centre in collaboration with the In Situ Thematic Centre of the Copernicus Marine Service (CMEMS INSTAC).  The observations integrated in the CORA-IBI product have been acquired both by autonomous platforms (Argo profilers, fixed moorings, gliders, drifters, sea mammals, fishery observing system from the RECOPESCA program), research or opportunity vessels ( CTDs, XBTs, ferrybox).  This CORA-IBI product has been controlled using an objective analysis (statistical tests) method and a visual quality control (QC). This QC procedure has been developed with the main objective to improve the quality of the dataset to the level required by the climate application and the physical ocean re-analysis activities. It provides T and S individual profiles on their original level with QC flags. The reference level of measurements is immersion (in meters) or pressure (in decibars). It is a subset on the IBI (Iberia-Bay-of-Biscay Ireland) of the CMEMS product referenced hereafter. The main new features of this regional product compared with previous global CORA products are the incorporation of coastal profiles from fishery observing system (RECOPESCA programme) in the Bay of Biscay and the English Channel as well as the use of an historical dataset collected by the Service hydrographique de la Marine (SHOM).