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.

'''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 operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' This RRS product is defined as the ratio of upwelling radiance and downwelling irradiance at 412, 443, 490, 510, 560 and 665 nm wavebands (corresponding to MERIS), and can also be expressed as the ratio of normalized water leaving Radiance (nLw) and the extra-terrestrial solar irradiance (F0). The ESA Climate Change Initiative is a 2-part programme aiming to produce “climate quality” merged data records from multiple sensors. The Ocean Colour project within this programme has a primary focus on chlorophyll in open oceans, using the highest quality Rrs merging process to date. This uses a combination of bandshifting to a reference sensor and temporally-weighted bias correction to align independent sensors into a coherent and minimally-biased set of reflectances. These are derived from level 2 data produced by SeaDAS l2gen (SeaWiFS) and Polymer (MODIS, VIIRS, MERIS and OLCI-3A) , and the resulting Rrs bias corrected. '''Processing information:''' ESA-CCI Rrs raw data are provided by Plymouth Marine Laboratory, currently at 4km resolution. These are processed to produce CMEMS representations 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 (70deg), large spacecraft zenith angle (56deg), coccolithophores, negative water leaving radiance, and normalized water leaving radiance at 560 nm 0.15 Wm-2 sr-1 (McClain et al., 1995). For the regional products, a variant of the OC-CCI chain is run to produce high resolution data at the 1km resolution necessary. '''DOI (product) :''' https://doi.org/10.48670/moi-00077

'''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). Additionally, the time series based on the method of von Schuckmann and Le Traon (2011) has been added. 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''' The global mean sea level is reflecting changes in the Earth’s climate system in response to natural and anthropogenic forcing factors such as ocean warming, land ice mass loss and changes in water storage in continental river basins. Thermosteric sea-level variations result from temperature related density changes in sea water associated with volume expansion and contraction (Storto et al., 2018). Global thermosteric sea level rise caused by ocean warming is known as one of the major drivers of contemporary global mean sea level rise (Cazenave et al., 2018; Oppenheimer et al., 2019). '''CMEMS KEY FINDINGS''' Since the year 2005 the upper (0-2000m) near-global (60°S-60°N) thermosteric sea level rises at a rate of 1.3±0.3 mm/year. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00240

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

'''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 distributes Remote Sensing Reflectances (Rrs) and diffuse attenuation coefficient of light at 490 nm (kd490) data. These datasets derived from Rrs multi-sensor (MODIS-AQUA, NOAA20-VIIRS, NPP-VIIRS, Sentinel3A-OLCI) spectra at the state-of-the-art algorithms for multi-sensor merging. 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. Reprocessed (multi-year) products are consistent and homogeneous in terms of format, algorithms and processing software. Rrs is defined as the ratio of upwelling radiance and downwelling irradiance at any wavelength (412, 443, 490, 555, and 670 nm). Kd490 is defined as 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 within the water column. It is directly related to the presence of scattering particles in the water column and is estimated through the ratio between Rrs at 490 and 555 nm. Kd490 is achieved via Mediterranean regional algorithm developed by GOS on the basis of MedBiOp in situ dataset (Volpe et al., 2019). The current day data temporal consistency is evaluated as Quality Index (QI): QI=(CurrentDataPixel-ClimatologyDataPixel)/STDDataPixel where QI is the difference between current data and the relevant climatological field as a signed multiple of climatological standard deviations (STDDataPixel). Inherent Optical Properties (aph443, adg443 and bbp443 at 443nm) are derived via QAAv6 model. '''Processing information:''' Multi-sensor product is 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. Single sensor NASA Level-2 data are destriped and then all Level-2 data are remapped at 1 km spatial resolution 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 and kd490, bbp, aph and adg) are estimated via state-of-the-art algorithms for better product quality. The entire data set is consistent and processed in one-shot mode (with an unique software version and identical configurations). '''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-l3-rrs412_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-rrs443_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-rrs490_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-rrs510_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-rrs555_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-rrs670_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-kd490_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-bbp443_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-adg443_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-aph443_1km_daily-rep-v02 '''Files format:''' *CF-1.4 *INSPIRE compliant '''DOI (product) :''' https://doi.org/10.48670/moi-00116

'''Short description: ''' For the Global Ocean - In-situ observation yearly delivery in delayed mode of Ocean surface currents. '''Detailed description: ''' The In Situ delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for Ocean surface currents. The data are collected from the Surface Drifter Data Assembly Centre (SD-DAC at NOAA AOML) completed by European data provided by EUROGOOS regional systems and national systems by the regional INS TAC components. All surface drifters data have been processed to check for drogue loss. Drogued and undrogued drifting buoy surface ocean currents are provided with a drogue presence flag as well as a wind slippage correction for undrogued buoy. '''Processing information: ''' From the near real time INS TAC product validated on a daily and weekly basis for forecasting purposes, and from the SD-DAC quality controlled dataset a scientifically validated product is created . It s a """"reference product"""" updated on a yearly basis. This product has been processed using a method that checks for drogue loss. Altimeter and wind data have been used to extract the direct wind slippage from the total drifting buoy velocities. The obtained wind slippage values have then been analyzed to identify probable undrogued data among the drifting buoy velocities dataset. A simple procedure has then been applied to produce an updated dataset including a drogue presence flag as well as a wind slippage correction. '''Suitability, Expected type of users / uses: ''' The product is designed to be assimilated into or for validation purposes of operational models operated by ocean forecasting centers for reanalysis purposes or for research community. These users need data aggregated and quality controlled in a reliable and documented manner.

'''This product has been archived''' '''Short description''': You can find here the OMEGA3D observation-based quasi-geostrophic vertical and horizontal ocean currents developed by the Consiglio Nazionale delle RIcerche. 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_REP_015_002 which corresponds to former version of MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012) and ERA-Interim surface fluxes. '''DOI (product) :''' https://doi.org/10.25423/cmcc/multiobs_glo_phy_w_rep_015_007

'''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 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). Additionally, the time series based on the method of von Schuckmann and Le Traon (2011) has been added. 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''' Ocean heat content (OHC) is defined here as the deviation from a reference period (1993-20210) and is closely proportional to the average temperature change from z1 = 0 m to z2 = 2000 m depth: With a reference density of ρ0 = 1030 kgm-3 and a specific heat capacity of cp = 3980 J/kg°C (e.g. von Schuckmann et al., 2009) Averaged time series for ocean heat content and their error bars are calculated for the Iberia-Biscay-Ireland region (26°N, 56°N; 19°W, 5°E). This OMI is computed using IBI-MYP, GLO-MYP reanalysis and CORA, ARMOR data from observations which provide temperatures. Where the CMEMS product for each acronym is: • IBI-MYP: IBI_MULTIYEAR_PHY_005_002 (Reanalysis) • GLO-MYP: GLOBAL_REANALYSIS_PHY_001_031 (Reanalysis) • CORA: INSITU_GLO_TS_OA_REP_OBSERVATIONS_013_002_b (Observations) • ARMOR: MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012 (Reprocessed observations) The figure comprises ensemble mean (blue line) and the ensemble spread (grey shaded). Details on the product are given in the corresponding PUM for this OMI as well as the CMEMS Ocean State Report: von Schuckmann et al., 2016; von Schuckmann et al., 2018. '''CONTEXT''' Change in OHC is a key player in ocean-atmosphere interactions and sea level change (WCRP, 2018) and can impact marine ecosystems and human livelihoods (IPCC, 2019). Additionally, OHC is one of the six Global Climate Indicators recommended by the World Meterological Organisation (WMO, 2017). In the last decades, the upper North Atlantic Ocean experienced a reversal of climatic trends for temperature and salinity. While the period 1990-2004 is characterized by decadal-scale ocean warming, the period 2005-2014 shows a substantial cooling and freshening. Such variations are discussed to be linked to ocean internal dynamics, and air-sea interactions (Fox-Kemper et al., 2021; Collins et al., 2019; Robson et al 2016). Together with changes linked to the connectivity between the North Atlantic Ocean and the Mediterranean Sea (Masina et al., 2022), these variations affect the temporal evolution of regional ocean heat content in the IBI region. Recent studies (de Pascual-Collar et al., 2023) highlight the key role that subsurface water masses play in the OHC trends in the IBI region. These studies conclude that the vertically integrated trend is the result of different trends (both positive and negative) contributing at different layers. Therefore, the lack of representativeness of the OHC trends in the surface-intermediate waters (from 0 to 1000 m) causes the trends in intermediate and deep waters (from 1000 m to 2000 m) to be masked when they are calculated by integrating the upper layers of the ocean (from surface down to 2000 m). '''CMEMS KEY FINDINGS''' The ensemble mean OHC anomaly time series over the Iberia-Biscay-Ireland region are dominated by strong year-to-year variations, and an ocean warming trend of 0.41±0.4 W/m2 is barely significant. '''DOI (product):''' https://doi.org/10.48670/mds-00316