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 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.

This delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for ocean surface currents and current vertical profiles. It concerns three delayed time datasets dedicated to near-surface currents measurements coming from three platforms (Lagrangian surface drifters, High Frequency radars and Argo floats) and velocity profiles within the water column coming from the Acoustic Doppler Current Profiler (ADCP, vessel mounted only). The latest version of Copernicus surface and sub-surface water velocity product is also distributed from Copernicus Marine catalogue.

This dataset provides a global Look-Up Table (LUT) of physiological ratios for the real-time adjustment of chlorophyll-a fluorescence measured by biogeochemical Argo (BGC-Argo) profiling floats. The physiological ratios aim to account for the global variability in the relationship between fluorescence and chlorophyll-a concentration, as influenced by phytoplankton physiology. The LUT was developed using two different gap-filled observational Argo-based products (SOCA machine learning-based methodology ; Sauzède et al., 2016; Sauzède et al., 2024). The first product provides gap-filled chlorophyll-a data derived from fluorescence corrected for dark signal and non-photochemical quenching (NPQ) following Schmechtig et al. (2023), while the second product provides chlorophyll-a concentrations derived from light attenuation. The latter is based on the downward irradiance at 490 nm (ED490) derived from the SOCA-light method (Renosh et al., 2023). From this, the diffuse attenuation coefficient (KD490) is computed, which is subsequently used to estimate the chlorophyll-a concentration through the bio-optical relationships described by Morel et al. (2007). These two products, based on fluorescence and radiometry, enable the derivation of spatially varying correction factors, or physiological ratios. These ratios provide a validated grounded framework for adjusting real-time fluorescence observations from OneArgo floats into chlorophyll-a concentrations. The LUT is distributed in NetCDF format and is provided on a regular 1°×1° latitude–longitude grid covering the global ocean. Each grid cell contains the temporal mean, averaged over the water column (from the surface to 1.5 times the euphotic depth), of the physiological ratio. The file also includes metadata describing variable definitions, units, and other relevant information. Variables included: - physiological_ratio — fluorescence-to-radiometry-based chlorophyll correction factor (dimensionless) - physiological_ratio_sd — temporal standard deviation (over the twelve months) of the fluorescence-to-radiometry-based chlorophyll correction factor (dimensionless) - lat, lon — spatial coordinates (degrees north/east) - Global attributes — dataset description, reference citation, and contact information

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.

Mesoscale eddy detection from 2000 to 2021 are computed using the AMEDA algorithm applied on AVISO L4 absolute dynamic topography at 1/8th degree. Eddy numbers correspond to tracks referenced in the DYNED atlas (https://doi.org/10.14768/2019130201.2). Detection is based on AVISO delyed-time product from 2000 to 2019 and on day+6 near-real-time altimetry from 2020 to 2021. Colocalisation is then made with available in situ profiles from Coriolis Ocean Dataset for Reanalysis (CORA) delayed-time data (113486 profiles) and Copernicus near-real-time profiles (43567).

The upper ocean pycnocline (UOP) monthly climatology is based on the ISAS20 ARGO dataset containing Argo and Deep-Argo temperature and salinity profiles on the period 2002-2020. Regardless of the season, the UOP is defined as the shallowest significant stratification peak captured by the method described in Sérazin et al. (2022), whose detection threshold is proportional to the standard deviation of the stratification profile. The three main characteristics of the UOP are provided -- intensity, depth and thickness -- along with hydrographic variables at the upper and lower edges of the pycnocline, the Turner angle and density ratio at the depth of the UOP. A stratification index (SI) that evaluates the amount of buoyancy required to destratify the upper ocean down to a certain depth, is also included. When evaluated at the bottom of the UOP, this gives the upper ocean stratification index (UOSI) as discussed in Sérazin et al. (2022). Three mixed layer depth variables are also included in this dataset, including the one using the classic density threshold of 0.03 kg.m-3, along with the minimum of these MLD variables. Several statistics of the UOP characteristics and the associated quantities are available in 2°×2° bins for each month of the year, whose results were smoothed using a diffusive gaussian filter with a 500 km scale. UOP characteristics are also available for each profile, with all the profiles sorted in one file per month.

This product integrates observations aggregated and validated from the Regional EuroGOOS consortium (Arctic-ROOS, BOOS, NOOS, IBI-ROOS, MONGOOS and Black Sea GOOS) as well as from National Data Centers (NODCs) and JCOMM global systems (Argo, GOSUD, OceanSITES, GTSPP, DBCP) and the Global telecommunication system (GTS) used by the Met Offices. Data are available in a dedicated directory to waves (INSITU_GLO_WAV_REP_OBSERVATIONS_013_045) of GLOBAL Distribution Unit in one file per platform. This directory is updated twice a year. Data are distributed in two datasets, one with original time sampling and the other with hourly data and rounded timestamps. The information distributed includes wave parameters and wave spectral information. The latest version of Copernicus delayed-mode wave product is distributed from Copernicus Marine catalogue. Additional credits: The American wave data are collected from US NDBC (National Data Buoy Center). The Australian wave data are collected from Integrated Marine Observing System (IMOS); IMOS is enabled by the National Collaborative Research Infrastructure Strategy (NCRIS); It is operated by a consortium of institutions as an unincorporated joint venture, with the University of Tasmania as Lead Agent. The Canadian data are collected from Fisheries and Oceans Canada.

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.

Satellite altimetry missions provide a quasi-global synoptic view of sea level over more than 25 years. The satellite altimetry constellation is used to build sea level maps and regional sea level indicators such as trends and accelerations. Estimating realistic uncertainties on these quantities is crucial to address some current climate science questions such as climate change detection and attribution or regional sea level budget closure for example. Previous studies have estimated the uncertainty for the global mean sea level (GMSL), but no uncertainty information is available at regional scales. In this study we estimate a regional satellite altimetry error budget and use it to derive maps of confidence intervals for local sea rise rates and accelerations. We analyze 27 years of satellite altimetry maps and derive the satellite altimetry error variance-covariance matrix at each grid point, prior to the estimation of confidence intervals on local trends and accelerations at the 90% confidence level using extended least squares estimators. Over 1993–2019, we find that the average local sea level trend uncertainty is 0.83 mm.yr-1 with local values ranging from 0.78 to 1.22 mm.yr-1. For accelerations, uncertainties range from 0.057 to 0.12 mm.yr-2, with a mean value of 0.063 mm.yr-2.   Change history: - 2020/07/08: initial dataset submission over 1993-2018 - 2020/10/21: 1993-2019 update and addition of error levels