'''DEFINITION''' The global yearly ocean CO2 sink represents the ocean uptake of CO2 from the atmosphere computed over the whole ocean. It is expressed in PgC per year. The ocean monitoring index is presented for the period 1985 to year-1. The yearly estimate of the ocean CO2 sink corresponds to the mean of a 100-member ensemble of CO2 flux estimates (Chau et al. 2022). The range of an estimate with the associated uncertainty is then defined by the empirical 68% interval computed from the ensemble. '''CONTEXT''' Since the onset of the industrial era in 1750, the atmospheric CO2 concentration has increased from about 277±3 ppm (Joos and Spahni, 2008) to 412.44±0.1 ppm in 2020 (Dlugokencky and Tans, 2020). By 2011, the ocean had absorbed approximately 28 ± 5% of all anthropogenic CO2 emissions, thus providing negative feedback to global warming and climate change (Ciais et al., 2013). The ocean CO2 sink is evaluated every year as part of the Global Carbon Budget (Friedlingstein et al. 2022). The uptake of CO2 occurs primarily in response to increasing atmospheric levels. The global flux is characterized by a significant variability on interannual to decadal time scales largely in response to natural climate variability (e.g., ENSO) (Friedlingstein et al. 2022, Chau et al. 2022). '''CMEMS KEY FINDINGS''' The rate of change of the integrated yearly surface downward flux has increased by 0.04±0.01e-1 PgC/yr2 over the period 1985 to year-1. The yearly flux time series shows a plateau in the 90s followed by an increase since 2000 with a growth rate of 0.06±0.04e-1 PgC/yr2. In 2021 (resp. 2020), the global ocean CO2 sink was 2.41±0.13 (resp. 2.50±0.12) PgC/yr. The average over the full period is 1.61±0.10 PgC/yr with an interannual variability (temporal standard deviation) of 0.46 PgC/yr. In order to compare these fluxes to Friedlingstein et al. (2022), the estimate of preindustrial outgassing of riverine carbon of 0.61 PgC/yr, which is in between the estimate by Jacobson et al. (2007) (0.45±0.18 PgC/yr) and the one by Resplandy et al. (2018) (0.78±0.41 PgC/yr) needs to be added. A full discussion regarding this OMI can be found in section 2.10 of the Ocean State Report 4 (Gehlen et al., 2020) and in Chau et al. (2022). '''DOI (product):''' https://doi.org/10.48670/moi-00223

"'Short description: ''' Global Ocean - This delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for ocean surface and subsurface currents. Current data from 5 different types of instruments are distributed: * The drifter's near-surface velocities computed from their position measurements. In addition, a wind slippage correction is provided from 1993. Information on the presence of the drogue of the drifters is also provided. * The near-surface zonal and meridional total velocities, and near-surface radial velocities, measured by High Frequency (HF) radars that are part of the European HF radar Network. These data are delivered together with standard deviation of near-surface zonal and meridional raw velocities, Geometrical Dilution of Precision (GDOP), quality flags and metadata. * The zonal and meridional velocities, at parking depth (mostly around 1000m) and at the surface, calculated along the trajectories of the floats which are part of the Argo Program. * The velocity profiles within the water column coming from Acoustic Doppler Current Profiler (vessel mounted ADCP, Moored ADCP, saildrones) platforms * The near-surface and subsurface velocities calculated from gliders (autonomous underwater vehicle) trajectories '''DOI (product) :''' https://doi.org/10.17882/86236

'''Short description:''' Global Ocean - near real-time (NRT) in situ quality controlled observations, hourly updated and distributed by INSTAC within 24-48 hours from acquisition in average. Data are collected mainly through global networks (Argo, OceanSites, GOSUD, EGO) and through the GTS '''DOI (product) :''' https://doi.org/10.48670/moi-00036

'''Short description:''' Altimeter satellite gridded Sea Level Anomalies (SLA) computed with respect to a twenty-year [1993, 2012] mean. The SLA is estimated by Optimal Interpolation, merging the L3 along-track measurement from the different altimeter missions available. Part of the processing is fitted to the European Sea area. (see QUID document or http://duacs.cls.fr [http://duacs.cls.fr] pages for processing details). The product gives additional variables (i.e. Absolute Dynamic Topography and geostrophic currents (absolute and anomalies)). It serves in delayed-time applications. This product is processed by the DUACS multimission altimeter data processing system. '''DOI (product):''' https://doi.org/10.48670/moi-00141

'''Short description:''' The Mean Dynamic Topography MDT-CMEMS_2024_EUR is an estimate of the mean over the 1993-2012 period of the sea surface height above geoid for the European Seas. This is consistent with the reference time period also used in the SSALTO DUACS products '''DOI (product) :''' https://doi.org/10.48670/mds-00337

'''DEFINITION''' Ocean salt content (OSC) is defined and represented here as the volume average of the integral of salinity in the Mediterranean Sea from z1 = 0 m to z2 = 300 m depth: ¯S=1/V ∫V S dV Time series of annual mean values area averaged ocean salt content are provided for the Mediterranean Sea (30°N, 46°N; 6°W, 36°E) and are evaluated in the upper 300m excluding the shelf areas close to the coast with a depth less than 300 m. The total estimated volume is approximately 5.7e+5 km3. '''CONTEXT''' The freshwater input from the land (river runoff) and atmosphere (precipitation) and inflow from the Black Sea and the Atlantic Ocean are balanced by the evaporation in the Mediterranean Sea. Evolution of the salt content may have an impact in the ocean circulation and dynamics which possibly will have implication on the entire Earth climate system. Thus monitoring changes in the salinity content is essential considering its link to changes in: the hydrological cycle, the water masses formation, the regional halosteric sea level and salt/freshwater transport, as well as for their impact on marine biodiversity. The OMI_CLIMATE_OSC_MEDSEA_volume_mean is based on the “multi-product” approach introduced in the seventh issue of the Ocean State Report (contribution by Aydogdu et al., 2023). Note that the estimates in Aydogdu et al. (2023) are provided monthly while here we evaluate the results per year. Six global products and a regional (Mediterranean Sea) product have been used to build an ensemble mean, and its associated ensemble spread. The reference products are: • The Mediterranean Sea Reanalysis at 1/24°horizontal resolution (MEDSEA_MULTIYEAR_PHY_006_004, DOI: https://doi.org/10.25423/CMCC/MEDSEA_MULTIYEAR_PHY_006_004_E3R1, Escudier et al., 2020) • Four global reanalyses at 1/4°horizontal resolution (GLOBAL_REANALYSIS_PHY_001_031, GLORYS, C-GLORS, ORAS5, FOAM, DOI: https://doi.org/10.48670/moi-00024, Desportes et al., 2022) • Two observation-based products: CORA (INSITU_GLO_TS_REP_OBSERVATIONS_013_001_b, DOI: https://doi.org/10.17882/46219, Szekely et al., 2022) and ARMOR3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012, DOI: https://doi.org/10.48670/moi-00052, Grenier et al., 2021). Details on the products are delivered in the PUM and QUID of this OMI. '''CMEMS KEY FINDINGS''' The Mediterranean Sea salt content shows a positive trend in the upper 300 m with a continuous increase over the period 1993-2021 at rate of 7.0*10-3 ±3.1*10-4 psu yr-1. The overall ensemble mean of different products is 38.57 psu. During the early 1990s in the entire Mediterranean Sea there is a large spread in salinity with the observational based datasets showing a higher salinity, while the reanalysis products present relatively lower salinity. The maximum spread between the period 1993–2021 occurs in the 1990s with a value of 0.12 psu, and it decreases to as low as 0.02 psu by the end of the 2010s. '''DOI (product):''' https://doi.org/10.48670/mds-00325

'''Short description:''' This product provides daily (nighttime), gap-free (Level-4, L4) maps of foundation Sea Surface Temperature (SST) - that is, the SST free from diurnal warming - over the Mediterranean Sea, at high (HR, 1/16°) and ultra-high (UHR, 1/100°) spatial resolutions, covering the period from 2008 to present. Each map represents nighttime SST values (centered at 00:00 UTC) and is produced by the Italian National Research Council – Institute of Marine Sciences (CNR-ISMAR). L4 maps are generated by selecting only the highest-quality SST observations from upstream Level-2 (L2) data acquired within a short local nighttime window, in order to minimize cloud contamination and avoid the effects of the diurnal cycle. The main L2 sources currently ingested include SLSTR from Sentinel-3A and -3B, VIIRS from NOAA-21, NOAA-20, and Suomi-NPP, AVHRR from Metop-B and -C, and SEVIRI. A two-step algorithm allows to interpolate SST data at high and ultra-high spatial resolution, applying statistical techniques (Buongiorno Nardelli et al., 2013; Buongiorno Nardelli et al., 2015). Additionally, from 2024 onwards, an improved first-guess field has been used in the generation of the MED UHR L4 data, enhancing the product's spatial resolution of SST features and the accuracy of SST gradients via machine learning techniques (Fanelli et al., 2024). '''DOI (product) :''' https://doi.org/10.48670/moi-00172

'''Short description:''' The global ocean objective analysis (gridded) fields of temperature and salinity are produced from the in situ profiles available in the multiparameter near–real-time in situ product INSITU_GLO_PHYBGCWAV_DISCRETE_MYNRT_013_030, with the exception of moorings, thermosalinographs and surface drifters. The objective analysis is based on the ISAS method (Gaillard et al. 2009), a statistical estimation approach that enables the mapping of ocean in situ profiles onto three-dimensional gridded fields. The resulting gridded product has a spatial resolution of 0.5° in latitude and 0.5° in longitude at the equator and includes 187 vertical levels. It provides monthly temperature and salinity fields centered the 15th of the month, with the analysis of the previous month delivered on the 8th of each month. A monthly file (data file) gathering the observed profiles used to calculate the analysis, which are interpolated on the vertical grid, is also provided. The product covers a two-year sliding time window. Older data are available in the corresponding delayed-mode product, INSITU_GLO_PHY_TS_OA_MY_013_052. '''DOI (product) :''' https://doi.org/10.48670/moi-00037

'''This product has been archived''' '''Short description:''' Near-Real-Time multi-mission global satellite-based spectral integral parameters. Only valid data are used, based on the L3 corresponding product. Included wave parameters are partition significant wave height, partition peak period and partition peak or principal direction. Those parameters are propagated in space and time at a 3-hour timestep and on a regular space grid, providing information of the swell propagation characteristics, from source to land. One file gathers one swell system, gathering observations originating from the same storm source. This product is processed by the WAVE-TAC multi-mission SAR data processing system to serve in near-real time the main operational oceanography and climate forecasting centers in Europe and worldwide. It processes data from the following SAR missions: Sentinel-1A and Sentinel-1B. All the spectral parameter measurements are optimally interpolated using swell observations belonging to the same swell field. The SAR data processing system produces wave integral parameters by partition (partition significant wave height, partition peak period and partition peak or principal direction) and the associated standard deviation and density of propagated observations. '''DOI (product) :''' https://doi.org/10.48670/moi-00175

'''Short description:''' In wavenumber spectra, the 1hz measurement error is the noise level estimated as the mean value of energy at high wavenumbers (below ~20km in term of wavelength). The 1hz noise level spatial distribution follows the instrumental white-noise linked to the Surface Wave Height but also connections with the backscatter coefficient. The full understanding of this hump of spectral energy (Dibarboure et al., 2013, Investigating short wavelength correlated errors on low-resolution mode altimetry, OSTST 2013 presentation) still remain to be achieved and overcome with new retracking, new editing strategy or new technology. '''DOI (product) :''' https://doi.org/10.48670/moi-00144