'''DEFINITION''' The time series are derived from the regional chlorophyll reprocessed (REP) products as distributed by CMEMS which, in turn, result from the application of the regional chlorophyll algorithms over remote sensing reflectances (Rrs) provided by the ESA Ocean Colour Climate Change Initiative (ESA OC-CCI, Sathyendranath et al. 2019; Jackson 2020). Daily regional mean values are calculated by performing the average (weighted by pixel area) over the region of interest. A fixed annual cycle is extracted from the original signal, using the Census-I method as described in Vantrepotte et al. (2009). The deseasonalised time series is derived by subtracting the mean seasonal cycle from the original time series, and then fitted to a linear regression to, finally, obtain the linear trend. '''CONTEXT''' Phytoplankton – and chlorophyll concentration as a proxy for phytoplankton – respond rapidly to changes in environmental conditions, such as temperature, light and nutrients availability, and mixing. The response in the North Atlantic ranges from cyclical to decadal oscillations (Henson et al., 2009); it is therefore of critical importance to monitor chlorophyll concentration at multiple temporal and spatial scales, in order to be able to separate potential long-term climate signals from natural variability in the short term. In particular, phytoplankton in the North Atlantic are known to respond to climate variability associated with the North Atlantic Oscillation (NAO), with the initiation of the spring bloom showing a nominal correlation with sea surface temperature and the NAO index (Zhai et al., 2013). '''CMEMS KEY FINDINGS''' While the overall trend average for the 1997-2021 period in the North Atlantic Ocean is slightly positive (0.16 ± 0.12 % per year), an underlying low frequency harmonic signal can be seen in the deseasonalised data. The annual average for the region in 2021 is 0.25 mg m-3. Though no appreciable changes in the timing of the spring and autumn blooms have been observed during 2021, a lower peak chlorophyll concentration is observed in the timeseries extension. This decrease in peak concentration with respect to the previous year is contributing to the reduction trend. '''DOI (product):''' https://doi.org/10.48670/moi-00194

'''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 Global Ocean. (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 near-real time applications. This product is processed by the DUACS multimission altimeter data processing system. '''DOI (product) :''' https://doi.org/10.48670/moi-00149

'''Short description:''' The IBI-MFC provides the ocean physical reanalysis multi year product for the Iberia-Biscay-Ireland (IBI) region starting in 01/01/1993, extended on yearly basis by using available reprocessed upstream data and regularly updated on monthly basis to cover the period up to month M-4 from present time using an interim processing system. The model system is designed, implemented and run by Mercator Ocean International, while the operational product post-processing and interim system are run by NOW Systems with the support of CESGA supercomputing centre. The IBI numerical core is based on the NEMO v3.6 ocean general circulation model, run at 1/36° horizontal resolution. Altimeter data, in situ temperature and salinity vertical profiles and satellite sea surface temperature are assimilated. The product offers 3D and 2D daily, monthly and yearly physical ocean fields, as well as hourly mean fields for surface variables. Additionally, climatological parameters (monthly mean and standard deviation) of these variables for the period 1993-2016 are delivered. '''DOI (Product)''': https://doi.org/10.48670/moi-00029

'''Short description:''' Altimeter satellite along-track sea surface heights anomalies (SLA) computed with respect to a twenty-year [1993, 2012] mean with a 1Hz (~7km) and 5Hz (~1km) sampling. It serves in near-real time applications. This product is processed by the DUACS multimission altimeter data processing system. It processes data from all altimeter missions available (e.g. Sentinel-6A, Jason-3, Sentinel-3A, Sentinel-3B, Saral/AltiKa, Cryosat-2, HY-2B). The system exploits the most recent datasets available based on the enhanced OGDR/NRT+IGDR/STC production. All the missions are homogenized with respect to a reference mission. Part of the processing is fitted to the European Seas. (see QUID document or http://duacs.cls.fr [http://duacs.cls.fr] pages for processing details). The product gives additional variables (e.g. Mean Dynamic Topography, Dynamic Atmospheric Correction, Ocean Tides, Long Wavelength Errors) that can be used to change the physical content for specific needs (see PUM document for details) '''Associated products''' A time invariant product http://marine.copernicus.eu/services-portfolio/access-to-products/?option=com_csw&view=details&product_id=SEALEVEL_GLO_PHY_NOISE_L4_STATIC_008_033 [http://marine.copernicus.eu/services-portfolio/access-to-products/?option=com_csw&view=details&product_id=SEALEVEL_GLO_PHY_NOISE_L4_STATIC_008_033] describing the noise level of along-track measurements is available. It is associated to the sla_filtered variable. It is a gridded product. One file is provided for the global ocean and those values must be applied for Arctic and Europe products. For Mediterranean and Black seas, one value is given in the QUID document. '''DOI (product) :''' https://doi.org/10.48670/moi-00140

'''Short description:''' MEDSEA_MULTIYEAR_WAV_006_012 is the multi-year wave product of the Mediterranean Sea Waves forecasting system (Med-WAV). It contains a Reanalysis dataset and a monthly climatological dataset (reference period 1993-2016). The Reanalysis dataset is a multi-year wave reanalysis starting from January 1985, composed by hourly wave parameters at 1/24° horizontal resolution, covering the Mediterranean Sea and extending up to 18.125W into the Atlantic Ocean. The dataset is extended every year as well as on a monthly basis through one-month extensions in interim mode, reaching one month before present. The Med-WAV modelling system is based on wave model WAM 4.6.2 and has been developed as a nested sequence of two computational grids (coarse and fine) to ensure that swell propagating from the North Atlantic (NA) towards the strait of Gibraltar is correctly entering the Mediterranean Sea. The coarse grid covers the North Atlantic Ocean from 75°W to 10°E and from 70° N to 10° S in 1/6° resolution while the nested fine grid covers the Mediterranean Sea from 18.125° W to 36.2917° E and from 30.1875° N to 45.9792° N with a 1/24° resolution. The modelling system resolves the prognostic part of the wave spectrum with 24 directional and 32 logarithmically distributed frequency bins. The wave system also includes an optimal interpolation assimilation scheme assimilating significant wave height along track satellite observations available through CMEMS and it is forced with daily averaged currents from Med-Physics and with 1-h, 0.25° horizontal-resolution ERA5 reanalysis 10m-above-sea-surface winds from ECMWF. ''DOI (Product)'': https://doi.org/10.48670/mds-00376 ''DOI (Interim dataset)'': https://doi.org/10.25423/ CMCC/MEDSEA_MULTIYEAR_WAV_006_012_MEDWAM3I ''DOI (climatological dataset)'': https://doi.org/10.25423/ CMCC/MEDSEA_MULTIYEAR_WAV_006_012_CLIM

'''DEFINITION''' Ocean heat content (OHC) is defined here as the deviation from a reference period (1993-2014) and is closely proportional to the average temperature change from z1 = 0 m to z2 = 700 m depth: OHC=∫_(z_1)^(z_2)ρ_0 c_p (T_yr-T_clim )dz [1] with a reference density of = 1030 kgm-3 and a specific heat capacity of cp = 3980 J kg-1 °C-1 (e.g. von Schuckmann et al., 2009). Time series of annual mean values area averaged ocean heat content is provided for the Mediterranean Sea (30°N, 46°N; 6°W, 36°E) and is evaluated for topography deeper than 300m. '''CONTEXT''' Knowing how much and where heat energy is stored and released in the ocean is essential for understanding the contemporary Earth system state, variability and change, as the oceans shape our perspectives for the future. The quality evaluation of MEDSEA_OMI_OHC_area_averaged_anomalies is based on the “multi-product” approach as introduced in the second issue of the Ocean State Report (von Schuckmann et al., 2018), and following the MyOcean’s experience (Masina et al., 2017). 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 degree 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 degree horizontal resolution (GLOBAL_MULTIYEAR_PHY_ENS_001_031): GLORYS, C-GLORS, ORAS5, FOAM • Two observation based products: CORA (INSITU_GLO_PHY_TS_OA_MY_013_052) and ARMOR3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012). Details on the products are delivered in the PUM and QUID of this OMI. '''CMEMS KEY FINDINGS''' The ensemble mean ocean heat content anomaly time series over the Mediterranean Sea shows a continuous increase in the period 1993-2022 at rate of 1.38±0.08 W/m2 in the upper 700m. After 2005 the rate has clearly increased with respect the previous decade, in agreement with Iona et al. (2018). '''DOI (product):''' https://doi.org/10.48670/moi-00261

'''Short description''': The data are provided weekly over a regular grid at 1/4° horizontal resolution, from the surface to 1500 m depth (representative of each Wednesday). The velocities are obtained by solving a diabatic formulation of the Omega equation, starting from ARMOR3D data (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012 ) and ERA5 surface fluxes. '''DOI (product) :''' https://doi.org/10.48670/moi-00053

'''Short description:''' For the Mediterranean Sea - The product contains daily Level-3 sea surface wind with a 1km horizontal pixel spacing using Near Real-Time Synthetic Aperture Radar (SAR) observations and their collocated European Centre for Medium-Range Weather Forecasts (ECMWF) model outputs. Products are updated several times daily to provide the best product timeliness. '''DOI (product) :''' https://doi.org/10.48670/mds-00334

'''Short description:''' For the Baltic Sea- The DMI Sea Surface Temperature L3S aims at providing daily multi-sensor supercollated data at 0.03deg. x 0.03deg. horizontal resolution, using satellite data from infra-red radiometers. Uses SST satellite products from these sensors: NOAA AVHRRs 7, 9, 11, 14, 16, 17, 18 , Envisat ATSR1, ATSR2 and AATSR. '''DOI (product) :''' https://doi.org/10.48670/moi-00154

'''Short description:''' For the Mediterranean Sea - The product contains daily Level-3 sea surface wind with a 1km horizontal pixel spacing using Synthetic Aperture Radar (SAR) observations and their collocated European Centre for Medium-Range Weather Forecasts (ECMWF) model outputs. Products are processed homogeneously starting from the L2OCN products. '''DOI (product) :''' https://doi.org/10.48670/mds-00342