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.

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

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' For the Global ocean, the ESA Ocean Colour CCI surface Chlorophyll (mg m-3, 4 km resolution) using the OC-CCI recommended chlorophyll algorithm is made available in CMEMS format. Phytoplankton functional types (PFT) products provide daily chlorophyll concentrations of three size-classes, consisting of nano, pico and micro-phytoplankton. L3 products are daily files, while the L4 are monthly composites. ESA-CCI data are provided by Plymouth Marine Laboratory at 4km resolution. These are processed using the same in-house software as in the operational processing. Standard masking criteria for detecting clouds or other contamination factors have been applied during the generation of the Rrs, i.e., land, cloud, sun glint, atmospheric correction failure, high total radiance, large solar zenith angle (actually a high air mass cutoff, but approximating to 70deg zenith), coccolithophores, negative water leaving radiance, and normalized water leaving radiance at 555 nm 0.15 Wm-2 sr-1 (McClain et al., 1995). 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. By comparing reflectances at different wavelengths and calibrating the result against in-situ measurements, an estimate of chlorophyll content can be derived. A detailed description of calibration & validation is given in the relevant QUID, associated validation reports and quality documentation. '''Processing information:''' ESA-CCI data are provided by Plymouth Marine Laboratory at 4km resolution. These are processed using the same in-house software as in the operational processing. The entire CCI data set is consistent and processing is done in one go. Both OC CCI and the REP product are versioned. Standard masking criteria for detecting clouds or other contamination factors have been applied during the generation of the Rrs, i.e., land, cloud, sun glint, atmospheric correction failure, high total radiance, large solar zenith angle (actually a high air mass cutoff, but approximating to 70deg zenith), coccolithophores, negative water leaving radiance, and normalized water leaving radiance at 555 nm 0.15 Wm-2 sr-1 (McClain et al., 1995). '''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. By comparing reflectances at different wavelengths and calibrating the result against in-situ measurements, an estimate of chlorophyll content can be derived. '''Quality / Accuracy / Calibration information:''' Detailed description of cal/val is given in the relevant QUID, associated validation reports and quality documentation.''' '''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. '''DOI (product) :''' https://doi.org/10.48670/moi-00097

'''This product has been archived''' For operational and online products, please visit https://marine.copernicus.eu '''Short description:''' For The Global Ocean - The GHRSST Multi-Product Ensemble (GMPE) system has been implemented at the Met Office which takes inputs from various analysis production centres on a routine basis and produces ensemble products at 0.25deg.x0.25deg. horizontal resolution. A large number of sea surface temperature (SST) analyses are produced by various institutes around the world, making use of the SST observations provided by the Global High Resolution SST (GHRSST) project. These are used by a number of groups including: numerical weather prediction centres; ocean forecasting groups; climate monitoring and research groups. There is a requirement to develop international collaboration in this field in order to assess and inter-compare the different analyses, and to provide uncertainty estimates on both the analyses and observational products. The GMPE system has been developed for these purposes and is run on a daily basis at the Met Office, producing global ensemble median and standard deviations for SST on a regular 1/4 degree resolution global grid. '''DOI (product) :''' https://doi.org/10.48670/moi-00166

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' Marine primary production corresponds to the amount of inorganic carbon which is converted into organic matter during the photosynthesis, and which feeds upper trophic layers. The daily primary production is estimated from satellite observations with the Antoine and Morel algorithm (1996). This algorithm modelized the potential growth in function of the light and temperature conditions, and with the chlorophyll concentration as a biomass index. The monthly area average is computed from monthly primary production weighted by the pixels size. The trend is computed from the deseasonalised time series (1998-2019), following the Vantrepotte and Mélin method. More details are provided in the Ocean State Reports 4 (Cossarini et al. ,2020). '''CONTEXT''' Marine primary production is at the basis of the marine food web and produce about 50% of the oxygen we breath every year (Behrenfeld et al., 2001). Study primary production is of paramount importance as ocean health and fisheries are directly linked to the primary production (Pauly and Christensen, 1995, Fee et al., 2019). Changes in primary production can have consequences on biogeochemical cycles, and specially on the carbon cycle, and impact the biological carbon pump intensity, and therefore climate (Chavez et al., 2011). Despite its importance for climate and socio-economics resources, primary production measurements are scarce and do not allow a deep investigation of the primary production evolution over decades. Satellites observations and modelling can fill this gap. However, depending of their parametrisation, models can predict an increase or a decrease in primary production by the end of the century (Laufkötter et al., 2015). Primary production from satellite observations present therefore the advantage to dispose an archive of more than two decades of global data. This archive can be assimilated in models, in addition to direct environmental analysis, to minimise models uncertainties (Gregg and Rousseaux, 2019). In the Ocean State Reports 4, primary production estimate from satellite and from modelling are compared at the scale of the Mediterranean Sea. This demonstrate the ability of such a comparison to deeply investigate physical and biogeochemical processes associated to the primary production evolution (Cossarini et al., 2020) '''CMEMS KEY FINDINGS''' The trend for the global ocean is negative over the period 1998-2019 with a decline in primary production of about 0.67 mgC.m-2.yr-1 or equivalently 0.2 %.yr-1. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00225

'''DEFINITION''' The global annual chlorophyll anomaly is computed by subtracting a reference climatology (1997-2014) from the annual chlorophyll mean, on a pixel-by-pixel basis and in log10 space. Both the annual mean and the climatology are computed employing ESA Ocean Colour Climate Change Initiative (ESA OC-CCI, Sathyendranath et al., 2018a) global products (i.e. using the standard OC-CCI chlorophyll algorithms, OCI) as distributed by CMEMS. '''CONTEXT''' Phytoplankton – and chlorophyll concentration as a proxy for phytoplankton – respond rapidly to changes in their physical environment. Some of those changes are seasonal and are determined by light and nutrient availability (Racault et al., 2012). By comparing annual mean values to a climatology, we effectively remove the seasonal signal, while retaining information on potential events during the year. Chlorophyll anomalies can be correlated to climate indexes in particular regions, such as the ENSO index in the equatorial Pacific (Behrenfeld et al. 2006; Racault et al., 2012) and the IOD index in the Indian Ocean (Brewin et al., 2012). It is important to study chlorophyll anomalies in consonance with sea surface temperature and sea level anomalies, as increases in chlorophyll are generally consistent with decreases in SST and sea level anomalies, suggesting an increase in mixing and vertical nutrient transport (von Schuckmann et al., 2016). '''CMEMS KEY FINDINGS''' The average global chlorophyll anomaly 2019 is -0.02 log10(mg m-3), with a maximum value of 1.7 log10(mg m-3) and a minimum value of -3.2 log10(mg m-3). That is to say that, in average, the annual 2019 mean value is slightly lower (96%) than the 1997-2014 climatological value. The positive signals reported in 2016 and 2017 (Sathyendranath et al., 2018b) in the southern Pacific Ocean could still be observed in the 2019 map, while the significant negative anomalies in the tropical waters of the northern Pacific Ocean were also detected to a lesser extent. Areas showing a change of anomaly sign from 2019 include the southern coast of Japan (no anomaly to positive) and the tropical Atlantic (anomalies close to zero for 2019). A marked increase in chlorophyll concentration was observed during 2019 in the Great Australian Bight, while negative anomalies became stronger in the Guatemala Basin and the region south of the Gulf of Guinea and, with values of chlorophyll reaching as low as 30% of the climatological value (anomaly < -0.5 log10(mg m-3)). The persistent positive anomalies in the higher latitudes of the North Atlantic (> 40°) match the cooling observed in the 2018 and previous years SST anomaly maps.

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

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

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