15 years of L-Band remote sensing Sea Surface Salinity (SSS) measurements have proven the capability of satellite SSS to resolve large scale to mesoscale SSS features in tropical to subtropical ocean. In mid to high latitude, L-Band measurements still suffer from large scale and time varying biases. Here, a simple method is proposed to mitigate the large scale and time varying biases. First, in order to estimate these biases, an Optimal Interpolation (OI) using a large correlation scale is used to map SMOS and SMAP L3 products and is compared to equivalent mapping of in situ observations. Then, a second mapping is performed on corrected SSS at scale of SMOS/SMAP resolution (~45 km). This procedure allows to correct and merge both products, and to increase signal to noise ratio of the absolute SSS estimates. Using thermodynamic equation of state (TEOS-10), the resulting L4 SSS product is combined with microwave satellite SST products to produce sea surface density and spiciness, useful to fully characterize the surface ocean water masses. The new L4 SSS products is validated against independent in situ measurements from low to high latitudes. The L4 products exhibits a significant improvement in mid-and high latitude in comparison to the existing SMOS and SMAP L3 products.

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.

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

This data set contains the gridded hydrographic and transport data for the biennial Go-Ship A25 Greenland–Portugal OVIDE section from 2002 to 2012. The properties and transports are mapped on a 7km x 1m grid. Using a common grid facilitates the comparison between the different occupations of the line and the averaging. This data set was used in Daniault et al. (2016, Progress in Oceanography) to which the reader is referred for a description of the gridding method.

The In Situ delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for ocean surface currents. The data are collected from the Surface Drifter Data Assembly Centre (SD-DAC at NOAA AOML). All surface drifters data have been processed to check for drogue loss. Drogued and undrogued drifting buoy surface ocean currents are provided with a drogue presence flag as well as a wind slippage correction for undrogued buoys. Altimeter and wind data have been used to extract the direct wind slippage from the total drifting buoy velocities. This product is designed to be assimilated into or for validation purposes of operational models operated by ocean forecasting centers for reanalysis purposes or for research community. These users need data aggregated and quality controlled in a reliable and documented manner.

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

Observations of Sea surface temperature and salinity are now obtained from voluntary sailing ships using medium or small size sensors. They complement the networks installed on research vessels or commercial ships. The delayed mode dataset proposed here is upgraded annually as a contribution to GOSUD (http://www.gosud.org )

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.

This dataset comprises two netcdf files. The first file contains the six global two-dimensional maps necessary to implement the tidal mixing parameterization presented in de Lavergne et al. (2020). Four power fields (E_wwi, E_sho, E_cri and E_hil) represent depth-integrated internal tide energy dissipation, with units of Watts per square meter. Each power field corresponds to a specific dissipative process and associated vertical structure of turbulence production. The two remaining fields, H_cri and H_bot, are decay heights (with units of meters) that enter the vertical structures of the E_cri and E_hil components, respectively. The second file contains three-dimensional fields of turbulence production (with units of Watts per kilogram) obtained by application of the parameterization to the WOCE global hydrographic climatology. The file includes the total turbulence production (epsilon_tid), its four components (epsilon_wwi, epsilon_sho, epsilon_cri, epsilon_hil), and the underlying hydrographic fields, as a function of longitude, latitude and depth. All maps have a horizontal resolution of 0.5º. Detailed documentation of the parameterization can be found in the following publication: de Lavergne, C., Vic, C., Madec, G., Roquet, F., Waterhouse, A.F., Whalen, C.B., Cuypers, Y., Bouruet-Aubertot, P., Ferron, B., Hibiya, T. A parameterization of local and remote tidal mixing. Journal of Advances in Modeling Earth Systems, 12, e2020MS002065 (2020). https://doi.org/10.1029/2020MS002065

The continuously updated version of Copernicus Argo floats realtime currents product is distributed from Copernicus Marine catalogue: - https://resources.marine.copernicus.eu/?option=com_csw&view=details&product_id=INSITU_GLO_UV_NRT_OBSERVATIONS_013_048 The Argo current product generated by Copernicus in situ TAC is derived from the original trajectory data from Argo GDAC (Global Data Assembly Center) available at: - Argo float data and metadata from Global Data Assembly Centre (Argo GDAC). SEANOE. https://doi.org/10.17882/42182 In 2021, the GDAC distributes data from more than 15,000 Argo floats. Deep ocean current is calculated from floats drift at parking depth, surface current is calculated from float surface drift. An Argo float drifts freely in the global ocean, performing regular observation cycles. An observation cycle usually spreads over 10 days :  - a surface descent to a parking depth (generally 1500 meters deep) - a 10-day drift at this parking depth - an ascent to the surface (vertical profile) - A short surface drift for data transmission The data transmitted at each cycle contain temperature, salinity observations (and additional biogeochemical parameters if applicable), positions (gps or argos), technical data. The ocean current product contains a NetCDF file for each Argo float. It is updated daily in real time by automated processes. For each cycle it contains the surface and deep current variables: - Date (time, time_qc) - Position  (latitude, longitude, position_qc) - Pressure (pres, pres_qc, representative_park_pressure for parking drift, 0 decibar for surface drift) - Current (ewct, ewct_qc, nsct, nsct_qc; the current vector is positioned and dated at the last position of the N-1 cycle) - Duration (days) of the current variable sampling (time_interval) - Grounded indicator - Positions and dates have a QC 1 (good data). Positions and dates that do not have a QC 1 are ignored. The positions are measured during the surface drift (Argos or GPS positioning). For the deep current of cycle N, we take the last good position of cycle N-1 and the first good position of cycle N. For the surface current of cycle N, we take the first and last good position of the N cycle.