he Global ARMOR3D L4 Reprocessed dataset is obtained by combining satellite (Sea Level Anomalies, Geostrophic Surface Currents, Sea Surface Temperature) and in-situ (Temperature and Salinity profiles) observations through statistical methods. References : - ARMOR3D: Guinehut S., A.-L. Dhomps, G. Larnicol and P.-Y. Le Traon, 2012: High resolution 3D temperature and salinity fields derived from in situ and satellite observations. Ocean Sci., 8(5):845–857. - ARMOR3D: Guinehut S., P.-Y. Le Traon, G. Larnicol and S. Philipps, 2004: Combining Argo and remote-sensing data to estimate the ocean three-dimensional temperature fields - A first approach based on simulated observations. J. Mar. Sys., 46 (1-4), 85-98. - ARMOR3D: Mulet, S., M.-H. Rio, A. Mignot, S. Guinehut and R. Morrow, 2012: A new estimate of the global 3D geostrophic ocean circulation based on satellite data and in-situ measurements. Deep Sea Research Part II : Topical Studies in Oceanography, 77–80(0):70–81.

'''DEFINITION''' Volume transport across lines are obtained by integrating the volume fluxes along some selected sections and from top to bottom of the ocean. The values are computed from models’ daily output. The mean value over a reference period (1993-2014) and over the last full year are provided for the ensemble product and the individual reanalysis, as well as the standard deviation for the ensemble product over the reference period (1993-2014). The values are given in Sverdrup (Sv). '''CONTEXT''' The ocean transports heat and mass by vertical overturning and horizontal circulation, and is one of the fundamental dynamic components of the Earth’s energy budget (IPCC, 2013). There are spatial asymmetries in the energy budget resulting from the Earth’s orientation to the sun and the meridional variation in absorbed radiation which support a transfer of energy from the tropics towards the poles. However, there are spatial variations in the loss of heat by the ocean through sensible and latent heat fluxes, as well as differences in ocean basin geometry and current systems. These complexities support a pattern of oceanic heat transport that is not strictly from lower to high latitudes. Moreover, it is not stationary and we are only beginning to unravel its variability. '''CMEMS KEY FINDINGS''' The mean transports estimated by the ensemble global reanalysis are comparable to estimates based on observations; the uncertainties on these integrated quantities are still large in all the available products. At Drake Passage, the multi-product approach (product no. 2.4.1) is larger than the value (130 Sv) of Lumpkin and Speer (2007), but smaller than the new observational based results of Colin de Verdière and Ollitrault, (2016) (175 Sv) and Donohue (2017) (173.3 Sv). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00247

'''This product has been archived''' For operational and online products, please visit https://marine.copernicus.eu '''Short description:''' For the Global Ocean - the OSTIA diurnal skin Sea Surface Temperature product provides daily gap-free maps of: *Hourly mean skin Sea Surface Temperature at 0.25° x 0.25° horizontal resolution, using in-situ and satellite data from infra-red radiometers. The Operational Sea Surface Temperature and Ice Analysis (OSTIA) system is run by the Met Office. A 1/4° (approx. 28 km) hourly analysis of skin Sea Surface temperature (SST) is produced daily for the global ocean. The skin temperature of the ocean is the temperature measured by satellite infra-red radiometers and can experience a large diurnal cycle. The skin SST L4 product is created by combining: 1. the OSTIA foundation SST analysis which uses in-situ and satellite observations; 2. the OSTIA diurnal warm layer analysis which uses satellite observations; and 3. a cool skin model. OSTIA uses satellite data provided by the GHRSST project. '''DOI (product) :''' https://doi.org/10.48670/moi-00167

'''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) sampling. It serves in delayed-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, Jason-1, Jason-2, Topex/Poseidon, ERS-1, ERS-2, Envisat, Geosat Follow-On, HY-2A, HY-2B, etc.). The system exploits the most recent datasets available based on the enhanced GDR/NTC production. All the missions are homogenized with respect to a reference mission. 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 (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 https://resources.marine.copernicus.eu/product-detail/SEALEVEL_GLO_PHY_NOISE_L4_STATIC_008_033/INFORMATION 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-00146

'''This product has been archived''' '''DEFINITION''' Estimates of Ocean Heat Content (OHC) are obtained from integrated differences of the measured temperature and a climatology along a vertical profile in the ocean (von Schuckmann et al., 2018). The regional OHC values are then averaged from 60°S-60°N aiming i) to obtain the mean OHC as expressed in Joules per meter square (J/m2) to monitor the large-scale variability and change. ii) to monitor the amount of energy in the form of heat stored in the ocean (i.e. the change of OHC in time), expressed in Watt per square meter (W/m2). Ocean heat content is one of the six Global Climate Indicators recommended by the World Meterological Organisation for Sustainable Development Goal 13 implementation (WMO, 2017). '''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 ocean shapes our perspectives for the future (von Schuckmann et al., 2020). Variations in OHC can induce changes in ocean stratification, currents, sea ice and ice shelfs (IPCC, 2019; 2021); they set time scales and dominate Earth system adjustments to climate variability and change (Hansen et al., 2011); they are a key player in ocean-atmosphere interactions and sea level change (WCRP, 2018) and they can impact marine ecosystems and human livelihoods (IPCC, 2019). '''CMEMS KEY FINDINGS''' Regional trends for the period 2005-2019 from the Copernicus Marine Service multi-ensemble approach show warming at rates ranging from the global mean average up to more than 8 W/m2 in some specific regions (e.g. northern hemisphere western boundary current regimes). There are specific regions where a negative trend is observed above noise at rates up to about -5 W/m2 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-00236

'''This product has been archived''' '''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''' '''DEFINITION''' Estimates of Ocean Heat Content (OHC) are obtained from integrated differences of the measured temperature and a climatology along a vertical profile in the ocean (von Schuckmann et al., 2018). The regional OHC values are then averaged from 60°S-60°N aiming i) to obtain the mean OHC as expressed in Joules per meter square (J/m2) to monitor the large-scale variability and change. ii) to monitor the amount of energy in the form of heat stored in the ocean (i.e. the change of OHC in time), expressed in Watt per square meter (W/m2). Ocean heat content is one of the six Global Climate Indicators recommended by the World Meterological Organisation for Sustainable Development Goal 13 implementation (WMO, 2017). '''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 ocean shapes our perspectives for the future (von Schuckmann et al., 2020). Variations in OHC can induce changes in ocean stratification, currents, sea ice and ice shelfs (IPCC, 2019; 2021); they set time scales and dominate Earth system adjustments to climate variability and change (Hansen et al., 2011); they are a key player in ocean-atmosphere interactions and sea level change (WCRP, 2018) and they can impact marine ecosystems and human livelihoods (IPCC, 2019). '''CMEMS KEY FINDINGS''' Since the year 2005, the near-surface (0-300m) near-global (60°S-60°N) ocean warms at a rate of 0.4 ± 0.1 W/m2. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00233

'''DEFINITION''' Significant wave height (SWH), expressed in metres, is the average height of the highest third of waves. This OMI provides global maps of the seasonal mean and trend of significant wave height (SWH), as well as time series in three oceanic regions of the same variables and their trends from 2002 to 2020, calculated from the reprocessed global L4 SWH product (WAVE_GLO_PHY_SWH_L4_MY_014_007). The extreme SWH is defined as the 95th percentile of the daily maximum SWH for the selected period and region. The 95th percentile is the value below which 95% of the data points fall, indicating higher than normal wave heights. The mean and 95th percentile of SWH (in m) are calculated for two seasons of the year to take into account the seasonal variability of waves (January, February and March, and July, August and September). Trends have been obtained using linear regression and are expressed in cm/yr. For the time series, the uncertainty around the trend was obtained from the linear regression, while the uncertainty around the mean and 95th percentile was bootstrapped. For the maps, if the p-value obtained from the linear regression is less than 0.05, the trend is considered significant. '''CONTEXT''' Grasping the nature of global ocean surface waves, their variability, and their long-term interannual shifts is essential for climate research and diverse oceanic and coastal applications. The sixth IPCC Assessment Report underscores the significant role waves play in extreme sea level events (Mentaschi et al., 2017), flooding (Storlazzi et al., 2018), and coastal erosion (Barnard et al., 2017). Additionally, waves impact ocean circulation and mediate interactions between air and sea (Donelan et al., 1997) as well as sea-ice interactions (Thomas et al., 2019). Studying these long-term and interannual changes demands precise time series data spanning several decades. Until now, such records have been available only from global model reanalyses or localised in situ observations. While buoy data are valuable, they offer limited local insights and are especially scarce in the southern hemisphere. In contrast, altimeters deliver global, high-quality measurements of significant wave heights (SWH) (Gommenginger et al., 2002). The growing satellite record of SWH now facilitates more extensive global and long-term analyses. By using SWH data from a multi-mission altimetric product from 2002 to 2020, we can calculate global mean SWH and extreme SWH and evaluate their trends, regionally and globally. '''KEY FINDINGS''' From 2002 to 2020, positive trends in both Significant Wave Height (SWH) and extreme SWH are mostly found in the southern hemisphere (a, b). The 95th percentile of wave heights (q95), increases faster than the average values, indicating that extreme waves are growing more rapidly than average wave height (a, b). Extreme SWH’s global maps highlight heavily storms affected regions, including the western North Pacific, the North Atlantic and the eastern tropical Pacific (a). In the North Atlantic, SWH has increased in summertime (July August September) but decreased in winter. Specifically, the 95th percentile SWH trend is decreasing by 2.1 ± 3.3 cm/year, while the mean SWH shows a decrease of 2.2 ± 1.76 cm/year. In the south of Australia, during boreal winter, the 95th percentile SWH is increasing at 2.6 ± 1.5 cm/year (c), with the mean SWH increasing by 0.5 ± 0.66 cm/year (d). Finally, in the Antarctic Circumpolar Current, also in boreal winter, the 95th percentile SWH trend is 3.2 ± 2.14 cm/year (c) and the mean SWH trend is 1.7 ± 0.84 cm/year (d). These patterns highlight the complex and region-specific nature of wave height trends. Further discussion is available in A. Laloue et al. (2024). '''DOI (product):''' https://doi.org/10.48670/mds-00352

'''DEFINITION''' The sea level ocean monitoring indicator has been presented in the Copernicus Ocean State Report #8. The sea level ocean monitoring indicator is derived from the DUACS delayed-time (DT-2024 version, “my” (multi-year) dataset used when available) sea level anomaly maps from satellite altimetry based on a stable number of altimeters (two) in the satellite constellation. The product is distributed by the Copernicus Climate Change Service and the Copernicus Marine Service (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). At each grid point, the trends/accelerations are estimated on the time series corrected from regional GIA correction (GIA map of a 27 ensemble model following Spada et Melini, 2019) and adjusted from annual and semi-annual signals. Regional uncertainties on the trends estimates can be found in Prandi et al., 2021. '''CONTEXT''' Change in mean sea level is an essential indicator of our evolving climate, as it reflects both the thermal expansion of the ocean in response to its warming and the increase in ocean mass due to the melting of ice sheets and glaciers(WCRP Global Sea Level Budget Group, 2018). According to the IPCC 6th assessment report (IPCC WGI, 2021), global mean sea level (GMSL) increased by 0.20 [0.15 to 0.25] m over the period 1901 to 2018 with a rate of rise that has accelerated since the 1960s to 3.7 [3.2 to 4.2] mm/yr for the period 2006–2018. Human activity was very likely the main driver of observed GMSL rise since 1970 (IPCC WGII, 2021). The weight of the different contributions evolves with time and in the recent decades the mass change has increased, contributing to the on-going acceleration of the GMSL trend (IPCC, 2022a; Legeais et al., 2020; Horwath et al., 2022). At regional scale, sea level does not change homogenously, and regional sea level change is also influenced by various other processes, with different spatial and temporal scales, such as local ocean dynamic, atmospheric forcing, Earth gravity and vertical land motion changes (IPCC WGI, 2021). The adverse effects of floods, storms and tropical cyclones, and the resulting losses and damage, have increased as a result of rising sea levels, increasing people and infrastructure vulnerability and food security risks, particularly in low-lying areas and island states (IPCC, 2019, 2022b). Adaptation and mitigation measures such as the restoration of mangroves and coastal wetlands, reduce the risks from sea level rise (IPCC, 2022c). '''KEY FINDINGS''' The altimeter sea level trends over the [1999/02/20 to 2025/10/18] period exhibit large-scale variations with trends up to +10 mm/yr in regions such as the western tropical Pacific Ocean. In this area, trends are mainly of thermosteric origin (Legeais et al., 2018; Meyssignac et al., 2017) in response to increased easterly winds during the last two decades associated with the decreasing Interdecadal Pacific Oscillation (IPO)/Pacific Decadal Oscillation (e.g., McGregor et al., 2012; Merrifield et al., 2012; Palanisamy et al., 2015; Rietbroek et al., 2016). Prandi et al. (2021) have estimated a regional altimeter sea level error budget from which they determine a regional error variance-covariance matrix and they provide uncertainties of the regional sea level trends. Over 1993-2019, the averaged local sea level trend uncertainty is around 0.83 mm/yr with local values ranging from 0.78 to 1.22 mm/yr. '''DOI (product):''' https://doi.org/10.48670/moi-00238

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