'''DEFINITION''' The indicator of the Kuroshio extension phase variations is based on the standardized high frequency altimeter Eddy Kinetic Energy (EKE) averaged in the area 142-149°E and 32-37°N and computed from the DUACS delayed-time (CMEMS SEALEVEL_GLO_PHY_L4_MY_008_047) and near real-time (CMEMS SEALEVEL_GLO_PHY_L4_NRT _008_046) altimeter sea level gridded products. ""CONTEXT"" The Kuroshio Extension is an eastward-flowing current in the subtropical western North Pacific after the Kuroshio separates from the coast of Japan at 35°N, 140°E. Being the extension of a wind-driven western boundary current, the Kuroshio Extension is characterized by a strong variability and is rich in large-amplitude meanders and energetic eddies (Niiler et al., 2003; Qiu, 2003, 2002). The Kuroshio Extension region has the largest sea surface height variability on sub-annual and decadal time scales in the extratropical North Pacific Ocean (Jayne et al., 2009; Qiu and Chen, 2010, 2005). Prediction and monitoring of the path of the Kuroshio are of huge importance for local economies as the position of the Kuroshio extension strongly determines the regions where phytoplankton and hence fish are located. Unstable (contracted) phase of the Kuroshio enhance the production of Chlorophyll (Lin et al., 2014). ""CMEMS KEY FINDINGS"" The different states of the Kuroshio extension phase have been presented and validated by (Bessières et al., 2013) and further reported by Drévillon et al. (2018) in the Copernicus Ocean State Report #2. Two rather different states of the Kuroshio extension are observed: an ‘elongated state’ (also called ‘strong state’) corresponding to a narrow strong steady jet, and a ‘contracted state’ (also called ‘weak state’) in which the jet is weaker and more unsteady, spreading on a wider latitudinal band. When the Kuroshio Extension jet is in a contracted (elongated) state, the upstream Kuroshio Extension path tends to become more (less) variable and regional eddy kinetic energy level tends to be higher (lower). In between these two opposite phases, the Kuroshio extension jet has many intermediate states of transition and presents either progressively weakening or strengthening trends. In 2018, the indicator reveals an elongated state followed by a weakening neutral phase since then. '''DOI (product):''' https://doi.org/10.48670/moi-00222

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

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' This product includes the Mediterranean Sea satellite chlorophyll trend map from 1997 to 2020 based on regional chlorophyll reprocessed (REP) product as distributed by CMEMS OC-TAC. This dataset, derived from multi-sensor (SeaStar-SeaWiFS, AQUA-MODIS, NOAA20-VIIRS, NPP-VIIRS, Envisat-MERIS and Sentinel3A-OLCI) (at 1 km resolution) Rrs spectra produced by CNR using an in-house processing chain, is obtained by means of the Mediterranean Ocean Colour regional algorithms: an updated version of the MedOC4 (Case 1 (off-shore) waters, Volpe et al., 2019, with new coefficients) and AD4 (Case 2 (coastal) waters, Berthon and Zibordi, 2004). The processing chain and the techniques used for algorithms merging are detailed in Colella et al. (2021). The trend map is obtained by applying Colella et al. (2016) methodology, where the Mann-Kendall test (Mann, 1945; Kendall, 1975) and Sens’s method (Sen, 1968) are applied on deseasonalized monthly time series, as obtained from the X-11 technique (see e. g. Pezzulli et al. 2005), to estimate, trend magnitude and its significance. The trend is expressed in % per year that represents the relative changes (i.e., percentage) corresponding to the dimensional trend [mg m-3 y-1] with respect to the reference climatology (1997-2014). Only significant trends (p < 0.05) are included. '''CONTEXT''' Phytoplankton are key actors in the carbon cycle and, as such, recognised as an Essential Climate Variable (ECV). Chlorophyll concentration - as a proxy for phytoplankton - respond rapidly to changes in environmental conditions, such as light, temperature, nutrients and mixing (Colella et al. 2016). The character of the response depends on the nature of the change drivers, and ranges from seasonal cycles to decadal oscillations (Basterretxea et al. 2018). The Mediterranean Sea is an oligotrophic basin, where chlorophyll concentration decreases following a specific gradient from West to East (Colella et al. 2016). The highest concentrations are observed in coastal areas and at the river mouths, where the anthropogenic pressure and nutrient loads impact on the eutrophication regimes (Colella et al. 2016). The the use of long-term time series of consistent, well-calibrated, climate-quality data record is crucial for detecting eutrophication. Furthermore, chlorophyll analysis also demands the use of robust statistical temporal decomposition techniques, in order to separate the long-term signal from the seasonal component of the time series. '''CMEMS KEY FINDINGS''' Chlorophyll trend in the Mediterranean Sea, for the period 1997-2020, is negative over most of the basin. Positive trend areas are visible only in the southern part of the western Mediterranean basin, in the Gulf of Lion, Rhode Gyre and partially along the Croatian coast of the Adriatic Sea. On average the trend in the Mediterranean Sea is about -0.5% per year. Nevertheless, as shown by Salgado-Hernanz et al. (2019) in their analysis (related to 1998-2014 satellite observations), there is not a clear difference between western and eastern basins of the Mediterranean Sea. In the Ligurian Sea, the trend switch to negative values, differing from the positive regime observed in the trend maps of both Colella et al. (2016) and Salgado-Hernanz et al. (2019), referred, respectively, to 1998-2009 and 1998-2014 time period, respectively. The waters offshore the Po River mouth show weak negative trend values, partially differing from the markable negative regime observed in the 1998-2009 period (Colella et al., 2016), and definitely moving from the positive trend observed by Salgado-Hernanz et al. (2019). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00260

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''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:''' The product MULTIOBS_GLO_PHY_SSS_L4_MY_015_015 is a reformatting and a simplified version of the CATDS L4 product called “SMOS-OI”. This product is obtained using optimal interpolation (OI) algorithm, that combine, ISAS in situ SSS OI analyses to reduce large scale and temporal variable bias, SMOS satellite image, SMAP satellite image, and satellite SST information. Kolodziejczyk Nicolas, Hamon Michel, Boutin Jacqueline, Vergely Jean-Luc, Reverdin Gilles, Supply Alexandre, Reul Nicolas (2021). Objective analysis of SMOS and SMAP Sea Surface Salinity to reduce large scale and time dependent biases from low to high latitudes. Journal Of Atmospheric And Oceanic Technology, 38(3), 405-421. Publisher's official version : https://doi.org/10.1175/JTECH-D-20-0093.1, Open Access version : https://archimer.ifremer.fr/doc/00665/77702/ '''DOI (product) :''' https://doi.org/10.48670/mds-00369

'''DEFINITION''' The temporal evolution of thermosteric sea level in an ocean layer is obtained from an integration of temperature driven ocean density variations, which are subtracted from a reference climatology to obtain the fluctuations from an average field. The products used include three global reanalyses: GLORYS, C-GLORS, ORAS5 (GLOBAL_MULTIYEAR_PHY_ENS_001_031) and two in situ based reprocessed products: CORA5.2 (INSITU_GLO_PHY_TS_OA_MY_013_052) , ARMOR-3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012). The regional thermosteric sea level values are then averaged from 60°S-60°N aiming to monitor interannual to long term global sea level variations caused by temperature driven ocean volume changes through thermal expansion as expressed in meters (m). '''CONTEXT''' Most of the interannual variability and trends in regional sea level is caused by changes in steric sea level. At mid and low latitudes, the steric sea level signal is essentially due to temperature changes, i.e. the thermosteric effect (Stammer et al., 2013, Meyssignac et al., 2016). Salinity changes play only a local role. Regional trends of thermosteric sea level can be significantly larger compared to their globally averaged versions (Storto et al., 2018). Except for shallow shelf sea and high latitudes (> 60° latitude), regional thermosteric sea level variations are mostly related to ocean circulation changes, in particular in the tropics where the sea level variations and trends are the most intense over the last two decades. '''CMEMS KEY FINDINGS''' Significant (i.e. when the signal exceeds the noise) regional trends for the period 2005-2023 from the Copernicus Marine Service multi-ensemble approach show a thermosteric sea level rise at rates ranging from the global mean average up to more than 8 mm/year. There are specific regions where a negative trend is observed above noise at rates up to about -5 mm/year such as in the subpolar North Atlantic, or the western tropical Pacific. These areas are characterized by strong year-to-year variability (Dubois et al., 2018; Capotondi et al., 2020). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00241

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' The Global Ocean Satellite monitoring and marine ecosystem study group (GOS) of the Italian National Research Council (CNR), in Rome operationally produces Level-4 product includes monthly averaged datasets of the diffuse attenuation coefficient of light at 490 nm (Kd490) for multi-sensor (MODIS-AQUA, NOAA20-VIIRS, NPP-VIIRS, Sentinel3A-OLCI at 300m of resolution) (at 1 km resolution) and Sentinel3A-OLCI observations (at 300m resolution). Kd490 is the diffuse attenuation coefficient of light at 490 nm, and is a measure of the turbidity of the water column, i.e., how visible light in the blue-green region of the spectrum penetrates the water column. It is directly related to the presence of absorbing and scattering matter in the water column and is estimated through the ratio between Rrs at 490 and 555 nm. For the multi-sensor dataset, single sensor Rrs fields are band-shifted, over the SeaWiFS native bands (using the QAAv6 model, Lee et al., 2002) and merged with a technique aimed at smoothing the differences among different sensors. This technique is developed by the GOS. The QAA allows the inversion of the radiative transfer equations to compute the Inherent Optical Properties. Level-4 product includes monthly averages along with the standard deviation and the number of observations in the period of integration. '''Processing information:''' Multi-sensor products are constituted by MODIS-AQUA, NOAA20-VIIRS, NPP-VIIRS and Sentinel3A-OLCI. For consistency with NASA L2 dataset, BRDF correction was applied to Sentinel3A-OLCI prior to band shifting and multi sensor merging. Hence, the single sensor OLCI data set is also distributed after BRDF correction. Single sensor NASA Level-2 data are destriped and then all Level-2 data are remapped at 1 km spatial resolution (300m for Sentinel3A-OLCI) using cylindrical equirectangular projection. Afterwards, single sensor Rrs fields are band-shifted, over the SeaWiFS native bands (using the QAAv6 model, Lee et al., 2002) and merged with a technique aimed at smoothing the differences among different sensors. This technique is developed by The Global Ocean Satellite monitoring and marine ecosystem study group (GOS) of the Italian National Research Council (CNR, Rome). Then geophysical fields (i.e. chlorophyll, kd490, bbp, aph and adg) are estimated via state-of-the-art algorithms for better product quality. Time averages are computed on the delayed-time data. '''Description of observation methods/instruments:''' Ocean colour technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so-called ocean colour which is affected by the presence of phytoplankton. '''Quality / Accuracy / Calibration information:''' A detailed description of the calibration and validation activities performed over this product can be found on the CMEMS web portal. '''Suitability, Expected type of users / uses:''' This product is meant for use for educational purposes and for the managing of the marine safety, marine resources, marine and coastal environment and for climate and seasonal studies. '''Dataset names :''' *dataset-oc-med-opt-multi-l4-kd490_1km_monthly-rt-v02 *dataset-oc-med-opt-olci-l4-kd490_300m_monthly-rt '''Files format:''' *CF-1.4 *INSPIRE compliant '''DOI (product) :''' https://doi.org/10.48670/moi-00117

'''DEFINITION''' The linear change of zonal mean subsurface temperature over the period 1993-2019 at each grid point (in depth and latitude) is evaluated to obtain a global mean depth-latitude plot of subsurface temperature trend, expressed in °C. The linear change is computed using the slope of the linear regression at each grid point scaled by the number of time steps (27 years, 1993-2019). A multi-product approach is used, meaning that the linear change is first computed for 5 different zonal mean temperature estimates. The average linear change is then computed, as well as the standard deviation between the five linear change computations. The evaluation method relies in the study of the consistency in between the 5 different estimates, which provides a qualitative estimate of the robustness of the indicator. See Mulet et al. (2018) for more details. '''CONTEXT''' Large-scale temperature variations in the upper layers are mainly related to the heat exchange with the atmosphere and surrounding oceanic regions, while the deeper ocean temperature in the main thermocline and below varies due to many dynamical forcing mechanisms (Bindoff et al., 2019). Together with ocean acidification and deoxygenation (IPCC, 2019), ocean warming can lead to dramatic changes in ecosystem assemblages, biodiversity, population extinctions, coral bleaching and infectious disease, change in behavior (including reproduction), as well as redistribution of habitat (e.g. Gattuso et al., 2015, Molinos et al., 2016, Ramirez et al., 2017). Ocean warming also intensifies tropical cyclones (Hoegh-Guldberg et al., 2018; Trenberth et al., 2018; Sun et al., 2017). '''CMEMS KEY FINDINGS''' The results show an overall ocean warming of the upper global ocean over the period 1993-2019, particularly in the upper 300m depth. In some areas, this warming signal reaches down to about 800m depth such as for example in the Southern Ocean south of 40°S. In other areas, the signal-to-noise ratio in the deeper ocean layers is less than two, i.e. the different products used for the ensemble mean show weak agreement. However, interannual-to-decadal fluctuations are superposed on the warming signal, and can interfere with the warming trend. For example, in the subpolar North Atlantic decadal variations such as the so called ‘cold event’ prevail (Dubois et al., 2018; Gourrion et al., 2018), and the cumulative trend over a quarter of a decade does not exceed twice the noise level below about 100m depth. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00244

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' Experimental altimeter satellite along-track sea surface heights anomalies (SLA) computed with respect to a twenty-year [1993, 2012] mean with a 5Hz (~1.3km) sampling. All the missions are homogenized with respect to a reference mission (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 Atmosphic Correction, Ocean Tides, Long Wavelength Errors, Internal tide, …) that can be used to change the physical content for specific needs This product was generated as experimental products in a CNES R&D context. It was processed by the DUACS multimission altimeter data processing system. '''DOI (product) :''' https://doi.org/10.48670/moi-00137

'''DEFINITION''' The product OMI_IBI_CURRENTS_VOLTRANS_section_integrated_anomalies is defined as the time series of annual mean volume transport calculated across a set of vertical ocean sections. These sections have been chosen to be representative of the temporal variability of various ocean currents within the IBI domain. The currents that are monitored include: transport towards the North Sea through Rockall Trough (RTE) (Holliday et al., 2008; Lozier and Stewart, 2008), Canary Current (CC) (Knoll et al. 2002, Mason et al. 2011), Azores Current (AC) (Mason et al., 2011), Algerian Current (ALG) (Tintoré et al, 1988; Benzohra and Millot, 1995; Font et al., 1998), and net transport along the 48ºN latitude parallel (N48) (see OMI Figure). To provide ensemble-based results, four Copernicus products have been used. Among these products are three reanalysis products (GLO-REA, IBI-REA and MED-REA) and one product obtained from reprocessed observations (GLO-ARM). • GLO-REA: GLOBAL_MULTIYEAR_PHY_001_030 (Reanalysis) • IBI-REA: IBI_MULTIYEAR_PHY_005_002 (Reanalysis) • MED-REA: MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012 (Reprocessed observations) • MED-REA: MEDSEA_MULTIYEAR_PHY_006_004MEDSEA_MULTIYEAR_PHY_006_004 (Reanalysis) The time series comprises the ensemble mean (blue line), the ensemble spread (grey shaded area), and the mean transport with the sign reversed (red dashed line) to indicate the threshold of anomaly values that would entail a reversal of the current transport. Additionally, the analysis of trends in the time series at the 95% confidence interval is included in the bottom right corner of each diagram. Details on the product are given in the corresponding Product User Manual (de Pascual-Collar et al., 2024a) and QUality Information Document (de Pascual-Collar et al., 2024b) as well as the CMEMS Ocean State Report: de Pascual-Collar et al., 2024c. '''CONTEXT''' The IBI area is a very complex region characterized by a remarkable variety of ocean currents. Among them, Podemos destacar las que se originan como resultado del closure of the North Atlantic Drift (Mason et al., 2011; Holliday et al., 2008; Peliz et al., 2007; Bower et al., 2002; Knoll et al., 2002; Pérez et al., 2001; Jia, 2000), las corrientes subsuperficiales que fluyen hacia el norte a lo largo del talud continental (de Pascual-Collar et al., 2019; Pascual et al., 2018; Machin et al., 2010; Fricourt et al., 2007; Knoll et al., 2002; Mazé et al., 1997; White & Bowyer, 1997). Y las corrientes de intercambio que se producen en el Estrecho de Gibraltar y el Mar de Alboran (Sotillo et al., 2016; Font et al., 1998; Benzohra and Millot, 1995; Tintoré et al., 1988). The variability of ocean currents in the IBI domain is relevant to the global thermohaline circulation and other climatic and environmental issues. For example, as discussed by Fasullo and Trenberth (2008), subtropical gyres play a crucial role in the meridional energy balance. The poleward salt transport of Mediterranean water, driven by subsurface slope currents, has significant implications for salinity anomalies in the Rockall Trough and the Nordic Seas, as studied by Holliday (2003), Holliday et al. (2008), and Bozec et al. (2011). The Algerian current serves as the sole pathway for Atlantic Water to reach the Western Mediterranean. '''CMEMS KEY FINDINGS''' The volume transport time series show periods in which the different monitored currents exhibited significantly high or low variability. In this regard, we can mention the periods 1997-1998 and 2014-2015 for the RTE current, the period 2012-2014 in the N48 section, the years 2006 and 2017 for the ALG current, the year 2021 for the AC current, and the period 2009-2012 for the CC current. Additionally, periods are detected where the anomalies are large enough (in absolute value) to indicate a reversal of the net transport of the current. This is the case for the years 1999, 2003, and 2012-2014 in the N48 section (with a net transport towards the north), the year 2017 in the ALC current (with net transport towards the west), and the year 2010 in the CC current (with net transport towards the north). The trend analysis of the monitored currents does not detect any significant trends over the analyzed period (1993-2022). However, the confidence interval for the trend in the RTE section is on the verge of rejecting the hypothesis of no trend. '''DOI (product):''' https://doi.org/10.48670/mds-00351