Worldwide, shellfish aquaculture and fisheries in coastal ecosystems represent crucial activities for human feeding. But these biological productions are under the pressure of climate variability and global change. Anticipating the biological processes affected by climate hazards remains a vital objective for species conservation strategies and human activities that rely on. Within marine species, filter feeders like oysters are real key species in coastal ecosystems due to their economic and societal value (fishing and aquaculture) but also due to their ecological importance. Indeed oysters populations in good health play the role of ecosystem engineers that can give many ecosystem services at several scales: building reef habitats that contribute to biodiversity, benthic-pelagic coupling and phytoplankton bloom control through water filtration, living shorelines against coastal erosion… The Pacific oyster, Crassostrea gigas (Thunberg, 1793), which is currently widespread worldwide, was introduced into the Atlantic European coasts at the end of the 19th century for shellfish culture purposes and becomes the main marine species farmed in France (around 100 000 tons) despite severe mortalities crisis. But in the same time and because of warming, natural oysters beds has spread significantly along the French coast and are supposed to have reach approximately 500 000 tons. In that context, Pacific oyster populations (natural and cultivated) in France are the subjects of many scientific projects. Among them, a specific long-term biological monitoring focuses on the reproduction of these populations at a national scale: the VELYGER national program. With more than 8 years of weekly data at many stations in France, this field-monitoring program offers a valuable dataset for studying processes underpinning reproduction cycle of this key-species in relation to environmental parameters, water quality and climate change.   Database content: Larval concentration (number of individuals per 1.5 m3) monitored, since 2008, at several stations in six bays of the French coast (from south to north): Thau Lagoon and bays of Arcachon, Marennes Oléron, Bourgneuf, Vilaine and Brest (see map below).   Methods used to monitor larval concentration: An important volume of seawater (1.5 m3) is pumped twice a week throughout the spawning season (june-september), at one meter below the surface at high tide (+/- 2h) in several sites within each VELYGER ecosystem. Water is filtered trough plankton net fitted with 40 µm mesh. After a proper rinsing of the net, the retained material is transferred into a polyethylene bottle (1 liter) and fixed with alcohol. At laboratory, sample is then gently filtered and rinse again and transferred into eprouvette. Two sub-samples of 1 mL are then taken using a pipette and examined on a graticule slide for microscope. The microscopic examination is made with a conventional binocular optical microscope with micrometer stage at a magnification of 10 X (or above). During the counting, a special care is necessary as larvae of other bivalves are also collected and confusion is possible. Larvae of C. gigas are also classified into four stage of development: - Stage I = D-shaped straight hinge larvae (shell length <105 µm) - Stage II = Early umbo evolved larvae (shell length between 105 and 150 µm) - Stage III = Medium umbo larvae (shell length between 150 and 235 µm) - Stage IV*= Large umbo eyed pediveliger larvae (shell length > 235 µm) * Larvae that are very closed to settle are sometimes identified into a separated 5th stage, but generally this stage is included in stage IV.   Illustrations: Location of the different Velyger sites along the French coast. From south to north: Thau Lagoon and bays of Arcachon, Marennes Oléron, Bourgneuf, Vilaine and Brest.   Legend: Pacific Oyster Larvae (left side) and Natural oyster bed (right side). Photos : © S. Pouvreau/Ifremer

For the 21 years of the study, an examination of trends in chlorophyll concentration revealed a general decline throughout the Gulf over the production period. These trends, extracted from dynamic linear model, also allowed this decline to be quantified. Expressed as a percentage, a large part of the area below the 50 m bathymetric line showed a decrease of at least 10% over the period, corresponding to a value of at least 0.1 µg.l-1. However, the spatial distribution reveals some more local phenomena. In southern Brittany, from Quimper to Vannes, a particular feature appears, with an upward trend over several kilometres along the coast, followed by a pronounced gradient along the coast. This gradient includes a zone where a continuous monotonic increasing trend is observed, then a zone where the trend becomes not significant and finally, about 15 km from the coast, a new zone where a significant continuous monotonic decreasing trend is observed. The increase in chlorophyll a concentration in the very coastal part is greater than 0.1 µg.l-1 over the period. Another peculiarity concerns the central part, located at the edge of the plateau at Cap Ferrat and Pente Aquitaine, where an increase in chlorophyll a was observed, but the variations remained small, being less than 0.1 µg.l-1. About a hundred kilometres south-west of Saint Nazaire, an area of about 40 by 50 km shows a decrease in chlorophyll a of more than 20%, quantified as more than 0.1 µg.l-1 over the period.

REPHYTOX dataset includes long-term time series on phycotoxins in marine bivalve molluscs, since 1987, along the whole French coast. The dataset covers results on lipophilic toxins, PSP toxins, ASP toxins, and palytoxins. REPHYTOX was a full part of the REPHY network until 2015. The whole dataset is available.

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, 2018a). This dataset comprises the biomass (in metric tons) and abundance (in thousands of individuals) of small pelagic fish estimated during the PELGAS survey in the Bay of Biscay in springtime. The dataset also includes the estimation coefficient of variation, derived based on the hydroacoustic methodology described in Doray et al. (2010), and the survey area. Those estimates have been validated by the ICES WGACEGG group and provided to the ICES WGHANSA group for stock assessment purposes. Data have been used in Doray et al., 2018b.

LOCEAN has been in charge of analyzing the isotopic composition of the dissolved inorganic carbon (DIC) in sea water collected during a series of cruises or ships of opportunity mostly in the southern Indian Ocean , the North Atlantic, and the equatorial Atlantic, but also in the Mediterranean Sea and in the equatorial Pacific. The LOCEAN sea-water samples for δ13CDIC were collected in 125/25 ml glass bottles until 2022/since then and poisoned with HgCl2 (1 ml of saturated solution) before storage in a dark room à 4°C until their measurement. The DIC was extracted from the seawater by acidification with phosphoric acid (H3PO4 85%) and CO2 gas that was produced was collected in a vacuum system following the procedure described by Kroopnick (1974). The isotopic composition of CO2 was determined using a dual inlet-isotopic ratio mass spectrometer (SIRA9-VG) by comparing the 13C/12C ratio of the sample to the 13C/12C ratio of a reference material, the Vienna-Pee Dee Belemnite (V-PDB). The isotopic composition is expressed in the δ-unit defined by Craig (1957)(method type 2).  Experience showed that samples older than 3-4 years are likely to have experienced conservation issues and have been dismissed. The mass spectrometer has worked very well until 2014-2015. Afterwards, its aging as well as the aging of the preparation line resulted in more data loss, and often less accurate results. The preparation line was renovated in 2019, and analyses in 2020 were run manually, often repeating the measurement a second time for each sample. Up to 2007-2008, δ13CDIC values have a precision of±0.01 ‰ (Vangriesheim et al.,2009) and a reproducibility of±0.02 ‰. After an interlaboratory comparison exercise led by Claire Normandeau (Dalhousie  University),  results  suggest  that  recent  LOCEAN  samples have a slightly poorer reproducibility (±0.04 ‰ ) as well as an offset of -0.13‰ (details available in Reverdin et al., ESSD 2018) that is confirmed by Becker et al. 2016 work by comparison with other cruises after removing the anthropogenic signal. Recent comparisons in early May 2021 with Orsay GEOPS facility samples suggest that the current offset is much smaller and might be +0.03‰. LOCEAN has installed in 2021 a new measurement device by coupling a Picarro G2131-I cavity ring down spectrometer (CRDS) with a CO2 extractor (Apollo SciTech) that will measure at the same time DIC (method type 3) (Leseurre, 2022). Since then, all water samples have been analyzed on this device. Part of the data set, as well as a scientific context and publications are also presented on the WEB site https://www.locean-ipsl.upmc.fr/oceans13c. Individual files correspond to regional subsets of the whole dataset. The file names are based on two letters for the region followed by (-) the cruise or project name (see below) followed by –DICisotopes, followed by either -s (surface data) or -b (subsurface data), and a version number (-V0, …): example SI-OISO-DICisotopes-s-V0; the highest version number corresponds to the latest update of the cruise/project data set, and can be directly downloaded. Earlier versions can be obtained on request, but are not recommended. The region two letters are the followings:   - SI: station and surface data in the Southern Indian Ocean that include cruises : INDIGO I (1985 – stn) (https://doi.org/10.17600/85000111) CIVA I (1993 – stn & surf) (https://doi.org/10.17600/93000870) (Archambeau et al., JMS 1998) ANTARES (1993 – stn & surf) (https://doi.org/10.17600/93000600) OISO (*) (since 1998 – stn & surf) (https://doi.org/10.18142/228) (Racapé et al., Tellus 2010, Leseurre, 2022)   - EA: station and surface data in the Tropical Atlantic Ocean that include cruises : EQUALANT (1999 & 2000 – surf) (https://doi.org/10.18142/98) EGEE (2005 to 2007 – stn & surf) (https://doi.org/10.18142/95) PIRATA (since 2013 – stn & surf) (https://doi.org/10.18142/14) EUMELI 2 (1991 – stn) (https://doi.org/10.17600/91004011)  (Pierre et al., JMS 1994) BIOZAIRE 3 (2003 – stn & surf ) (https://doi.org/10.17600/3010120) (Vangriesheim et al., DSRII, 2009) TARA-Microbiomes (2021 - stn & surf)   - NA : station and surface data in the North Atlantic Subpolar gyre that include cruises : OVIDE (**) (since 2002 – stn & surf) (https://doi.org/10.17882/46448) (Racapé et al., 2013) RREX (2017 – stn & surf) (https://doi.org/10.17600/17001400) SURATLANT (since 2010 - surf) (https://doi.org/10.17882/54517) (Racapé et al., BG 2014 ; Reverdin et al., ESSD 2018, Leseurre, 2022) NUKATUKUMA (since 2017- surf)   - MS: station data in the Mediterranean sea that include cruises : ALMOFRONT 1 (1991 – stn) (https://doi.org/10.17600/91004211) VICOMED 3 (1990 – stn) (https://doi.org/10.17600/90000711)   - PO: tropical Pacific that include cruises : PANDORA (2012 – stn) (https://doi.org/10.17600/12010050) ALIZE2 (1991 – stn & surf) (https://doi.org/10.17600/91002711) (Laube-Lenfant and Pierre, Oceanologica Acta 1994)   - SO: station and surface data in the Southern Ocean (except OISO) that include cruises: TARA-Microbiomes (2021-2022, stn & surf) AGULHASII-072022 (2022, stn) CONFLUENCE (1993-1994, stn)   - AO: station and surface data in the Arctic Ocean and nearby seas that include cruises: GREENFEEDBACK (2024, stn&surf) TCA (2024, stn) REFUGE ARCTIC (2024, stn) (*) The values for cruises OISO19, 21 and 22 are doubtful (for some, too low) and will require further investigation to find whether adjusted values can be proposed. (**) Some of the OVIDE cruises are also referred to as or GEOVIDE (in 2014), and BOCATS (in 2016). CATARINA, BOCATS1 and BOCATS2 (PID2019-104279GB-C21/AEI/10.13039/501100011033) cruises were funded by the Spanish Research Agency  The values of the OVIDE 2010 stations are doubtful (too low), but no particular error was found, and they have been left in the files.   Data The files are in csv format reported as: - Cruise name, station id, (bottle number), day, month, year, hour, minute, longitude, latitude, pressure (db), depth (m), temperature (°C), temperature qc, salinity (pss-78), salinity qc, d13CDIC, d13CDIC qc, method type - Temperature is an in situ temperature - Salinity is a practical salinity - Method type (1) acid CO2 extraction from helium stripping technique coupled to mass spectrometer, (2) acid CO2 extraction in a vacuum system coupled to mass spectrometer,(3) CO2 extractor (Apollo SciTech) coupled to CRDS measurements. Temperature qc, salinity qc, d13CDIC qc are quality indices equal to: - 0 no quality check (but presumably good data) - 1 probably good data - 2 good data - 3 probably bad data - 4 certainly bad data - 9 missing data (and the missing data are reported with an unlikely missing value)

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  

SOMLIT (Service d'Observation en Milieur Littoral) : a French Coastal Monitoring Network Coastal zones are where land, ocean and atmosphere interact. They are important for the exchange of matter and energy, and play a key role in (biogeo)chemical cycles at global scale. These environments are characterised by significant spatial and temporal variability of their physico-chemical and biological parameters due to local and seasonal meteorological drivers which are exacerbated by large-scale climate drivers (e.g. global warming, modification of the wind regime) and local-scale anthropogenic drivers (e.g. nutrient cycle changes linked to the use of fertilisers or the construction of large installations such as dams). These driving mechanisms are often interconnected. In the context of global warming (due to­­ climate and human-induced changes), the identification and understanding of their impact on coastal marine and littoral ecosystems is essential. The scientific objective of SOMLIT is to 1) characterise the multi-decadal evolution of coastal marine and littoral ecosystems, and 2) determine the climatic and anthropogenic drivers. In order to meet this objective, a nationally coordinated multi-site monitoring system was set up in the mid-1990s. The observation strategy is the same for each of the 12 monitored ecosystems with fortnightly sampling and/or measurements, at high tide (for sites subject to tides): 1) in surface-water for a range of 15 parameters (temperature, salinity, dissolved oxygen, pH, nitrate, nitrite, ammonium, phosphate, silicate, suspended particulate matter, chlorophyll a, particulate organic carbon and nitrogen and stable isotopes of particulate organic carbon and nitrogen), 2) in surface-water for a range of 26 parameters of numbering and optical characteristics of pico- and nanoplankton), and 3) along the water column for temperature, salinity, fluorescence and PAR (vertical profiles of multi-parameter probes). SOMLIT’s activities are carried out under a quality assurance / quality control process based on the ISO 17025 standard. SOMLIT’s service provision objectives are to provide data and logistical support for research and other observation activities. SOMLIT has been officially accredited since 1996 as one of the CNRS (French National Centre for Scientific Research) National Observation Services (SNO). SOMLIT’s coordination is hosted by the Observatoire Aquitain des Sciences de l'Univers (University of Bordeaux / CNRS) and the service relies on strong partnerships with nine other institutions (University of Lille, University of the Littoral Opal Coast, University of Caen Normandy, Sorbonne University, University of Western Brittany, La Rochelle University, University of Montpellier, Aix Marseille University, National Museum of Natural History). SOMLIT is one of the nine networks that compose France’s Coastal Research Infrastructure (ILICO).  SOMLIT has strong ties with ILICO’s other networks such as the SNOs MOOSE (Mediterranean Ocean Observing System on Environment), PHYTOBS (microphytoplankton monitoring) and COAST-HF (Coastal Ocean Observing System - High Frequency).

Understanding the spatial and temporal preferences of toxic phytoplankton species is of paramount importance in managing and predicting harmful events in aquatic ecosystems. In this study we address the realised niche of the species Alexandrium minutum, Pseudo-nitzschia fraudulenta and P. australis. We used them to highlight distribution patterns at different scales and determine possible drivers. To achieve this, we have developed original procedures coupling niche theory and habitat suitability modelling using abundance data in four consecutive steps: 1) Estimate the realised niche applying kernel functions. 2) Assess differences between the species’ niche as a whole and at the local level. 3) Develop habitat and temporal suitability models using niche overlap procedures. 4) Explore species temporal and spatial distributions to highlight possible drivers. Data used are species abundance and environmental variables collected over 27 years (1988-2014) and include 139 coastal water sampling sites along the French Atlantic coast. Results show that A. minutum and P. australis niches are very different, although both species have preference for warmer months. They both respond to decadal summer NAO but in the opposite way. P. fraudulenta realised niche lies in between the two other species niches. It also prefers warmer months but does not respond to decadal summer NAO. The Brittany peninsula is now classified as an area of prevalence for the three species. The methodology used here will allow to anticipate species distribution in the event of future environmental challenges resulting from climate change scenarios.

This database contains hauls collated from 1965 to 2019, from fisheries dependent and independent data, from across eastern Atlantic waters and French Mediterranean waters. From this data diadromous fish spatio-temporal data was cleaned and standardised.

This dataset contains bio-optical measurements from BioGeoChemical-Argo (BGC-Argo) profiling floats complemented with ocean-colour satellite matchups of variables related to the detection of coccolithophore blooms dominated by Emiliania huxleyi. BGC-Argo float data cover the global ocean from November 2012 to December 2018 and include measurements of the particulate backscattering coefficient (BBP_float in m-1), the concentration of Chlorophyll-a (CHLA_float in mg m-3), and the particulate beam attenuation coefficient (CP_float in m-1) with data processing and quality control described in the manuscript entitled “Detection of coccolithophore blooms with BioGeoChemical-Argo floats” submitted to Geophysical Research Letters. The data represent near-surface ocean conditions, calculated as the average value in the top 15m of the water column. Daily ocean-colour satellite data were downloaded from the GlobColour project (ftp://ftp.hermes.acri.fr) with a spatial resolution of 4km and matched with every BGC-Argo float observation by using a 5x5 pixel box and a 9-day temporal window. For each float observation, we extracted concurrent satellite data of the concentrations of Particulate Inorganic Carbon (PIC_sat in mmol m-3) and Particulate Organic Carbon (POC_sat in mmol m-3), from which we derived the proportion of PIC_sat to the total particulate carbon concentration (PIC_POC_sat in % and defined as PIC_sat / [PIC_sat+POC_sat]). Coccolithophore bloom periods were identified using annual times series of PIC_sat and PIC_POC_sat at each profile location as described in the submitted manuscript, and the column “inside_coccolithophore_bloom” reports the float observations occurring inside such blooms.