This dataset contains the dynamical outputs of a global ocean simulation coupling dynamics and biogeochemistry at ¼° over the year 2019. The simulation has been performed using the coupled circulation/ecosystem model NEMO/PISCES (https://www.nemo-ocean.eu/), which is here enhanced to perform an ensemble simulation with explicit simulation of modeling uncertainties in the physics and in the biogeochemistry. This dataset is one of the 40 members of the ensemble simulation. This study was part of the Horizon Europe project SEAMLESS (https://seamlessproject.org/Home.html), with the general objective of improving the analysis and forecast of ecosystem indicators.   See Popov et al. (https://os.copernicus.org/articles/20/155/2024/) for more details on the study.

This dataset contains (1) outputs of Species Distribution Models (SDMs) for marine fish and (2) estimation of species richness using those outputs. SDMs use correlative algorithms to link presences of species to the environment recorded in place and time of their observation, calculate their environmental niche, estimate the geographical location suitable for them (habitat suitability) and in fine their geographical distribution. Here, we downloaded presences of marine fish from two open source databases, GBIF and OBIS and 13 environmental predictors known to be relevant in fish ecology (list below). We used the CEPHALOPOD pipeline, a framework allowing the user the compute a lot of species at the same time, with comparable methods and a verification of quality of inputs and outputs at every steps (Schickele et al., 2025). 3,642 fish made it to the final step and have the habitat suitability estimated for 12 month + annual mean, 10 bootstrap to quantify uncertainty and x algorithms. Those estimation are available in the “L2_marine_fish_*.nc” files, organized by water column position (bathydemersal, bathypelagic, benthopelagic, demersal, pelagic-neritic, pelagic-oceanic, reef-associated). We then used their annual mean to estimate their actual geographic distribution by applying (1) a bathymetric filtration and (2) a cutting procedure which removes isolated patches of high suitability (i.e., potential distribution) with no recorded presences (i.e., considered to be outside of the species dispersion range). Those geographic distribution were then stacked to estimate global species richness of every fish and each water column position.

These monthly gridded climatology were produced using MBT, XBT, Profiling floats, Gliders, and ship-based CTD data from different database and carried out in the Med. between 1969 and 2013. The Mixed Layer Depth (MLD) is calculated with a delta T= 0.1 C criterion relative to 10m reference level on individual profiles. The Depth of the Bottom of the Seasonal Thermocline (DBST) is calculated on individual profiles as the maximum value from a vector composed of two elements: 1) the depth of the temperature minimum in the upper 200m; 2) the MLD. This double criterion for the calculation of DBST is necessary in areas where the mixed layer exceed 200m depth. DBST is the integration depth used in the calculation of the upper-ocean Heat Storage Rate. For more details about the data and the methods used, see: Houpert et al. 2015, Seasonal cycle of the mixed layer, the seasonal thermocline and the upper-ocean heat storage rate in the Mediterranean Sea derived from observations, Progress in Oceanography, http://doi.org/10.1016/j.pocean.2014.11.004

This dataset is an aggregation of all availale in situ data from Coriolis and Copernicus in situ data centres, observed in the French DCSMM area. It contains 5167 NetCDF CF files from 1903 to 2017. Each file contains the observations of a specific platform (e.g. vessel, mooring site, sea level station). Observed parameters are temperature, salinity, pressure, oxygen, nitrate, chlorophyll (and other bio-geo-chemicals), current, wave, sea level, river flow.  

This data set provides a monthly time series of the upper limb of the Meridional Overturning Circulation (MOC) intensity at the A25 Greenland-Portugal OVIDE line from 1993 to 2015. The MOC was derived by combining AVISO altimetry with ISAS temperature and salinity data. The reader is referred to Mercier et al. (2015, Progress in Oceanography) for a full description of the method.

This dataset is composed by the climatological seasonal field of the Ocean Salinity Stratification as defined from the Brunt-Vaisala frequency limited to the upper 300 m depth. The details are given in Maes, C., and T. J. O’Kane (2014), Seasonal variations of the upper ocean salinity stratification in the Tropics, J. Geophys. Res. Oceans, 119, 1706–1722, doi:10.1002/2013JC009366.

This product contains observations and gridded files from two up-to-date carbon and biogeochemistry community data products: Surface Ocean Carbon ATlas SOCATv2023 and GLobal Ocean Data Analysis Project GLODAPv2.2023. The SOCATv2023-OBS dataset contains >25 million observations of fugacity of CO2 of the surface global ocean from 1957 to early 2023. The quality control procedures are described in Bakker et al. (2016). These observations form the basis of the gridded products included in SOCATv2023-GRIDDED: monthly, yearly and decadal averages of fCO2 over a 1x1 degree grid over the global ocean, and a 0.25x0.25 degree, monthly average for the coastal ocean. GLODAPv2.2023-OBS contains >1 million observations from individual seawater samples of temperature, salinity, oxygen, nutrients, dissolved inorganic carbon, total alkalinity and pH from 1972 to 2021. These data were subjected to an extensive quality control and bias correction described in Olsen et al. (2020). GLODAPv2-GRIDDED contains global climatologies for temperature, salinity, oxygen, nitrate, phosphate, silicate, dissolved inorganic carbon, total alkalinity and pH over a 1x1 degree horizontal grid and 33 standard depths using the observations from the previous major iteration of GLODAP, GLODAPv2. SOCAT and GLODAP are based on community, largely volunteer efforts, and the data providers will appreciate that those who use the data cite the corresponding articles (see References below) in order to support future sustainability of the data products.

Argo is a global array of 3,000 free-drifting profiling floats that measures the temperature and salinity of the upper 2000 m of the ocean. This allows, for the first time, continuous monitoring of the temperature, salinity, and velocity of the upper ocean, with all data being relayed and made publicly available within hours after collection. The array provides 100,000 temperature/salinity profiles and velocity measurements per year distributed over the global oceans at an average of 3-degree spacing. Some floats provide additional bio-geo parameters such as oxygen or chlorophyll. All data collected by Argo floats are publically available in near real-time via the Global Data Assembly Centers (GDACs) in Brest (France) and Monterey (California) after an automated quality control (QC), and in scientifically quality controlled form, delayed mode data, via the GDACs within six months of collection. The BGC-Argo Sprof snapshot is a subset of the global Argo data snapshot.  It is created to ease BGC-Argo data usage.  The content is the same if you are to download the global Argo data snapshot, and then select all the BGC-Argo Sprof files.  Please use the same DOI and citation as the global Argo data snapshot.

Ensemble simulations of the ecosystem model Apecosm (https://apecosm.org) forced by the IPSL-CM6-LR climate model with the climate change scenario SSP5-8.5. The output files contain yearly mean biomass density for 3 communities (epipelagic, mesopelagic migratory and mesopelagic redidents) and 100 size classes (ranging from 0.12cm to 1.96m) The model grid file is also provided. Units are in J/m2 and can be converted in kg/m2 by dividing by 4e6. These outputs are associated with the "Assessing the time of emergence of marine ecosystems from global to local scales using IPSL-CM6A-LR/APECOSM climate-to-fish ensemble simulations" paper from the Earth's Future "Past and Future of Marine Ecosystems" Special Collection.

Ensemble simulations of the ecosystem model Apecosm (https://apecosm.org) forced by the IPSL-CM6-LR climate model with the climate change scenario SSP1-2.6. The output files contain yearly mean biomass density for 3 communities (epipelagic, mesopelagic migratory and mesopelagic redidents) and 100 size classes (ranging from 0.12cm to 1.96m) The model grid file is also provided. Units are in J/m2 and can be converted in kg/m2 by dividing by 4e6. These outputs are associated with the "Assessing the time of emergence of marine ecosystems from global to local scales using IPSL-CM6A-LR/APECOSM climate-to-fish ensemble simulations" paper from the Earth's Future "Past and Future of Marine Ecosystems" Special Collection.