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. 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-00101

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

'''This product has been archived''' "''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-2022), following the Vantrepotte and Mélin (2009) method. The trend estimate is not shown because the length of the time series does not allow to completely differentiate the climate trend to the natural variability of the primary production. 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 presents 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 demonstrates 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''' Global primary production does not show specific trend and remain relatively constant over the archive 1998-2022. The temporal variability of the primary production appears to be mainly driven by the seasonal variation. However, some specific inter-annual event may induce noticeable increase or decrease in primary production, as for example in the second part of 2011. '''DOI (product):''' https://doi.org/10.48670/moi-00225

'''DEFINITION''' The time series are derived from the regional chlorophyll reprocessed (MY) product as distributed by CMEMS (OCEANCOLOUR_MED_BGC_L3_NRT_009_141). This dataset, derived from multi-sensor (SeaStar-SeaWiFS, AQUA-MODIS, NOAA20-VIIRS, NPP-VIIRS, Envisat-MERIS and Sentinel3-OLCI) 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. (2023). Monthly regional mean values are calculated by performing the average of 2D monthly mean (weighted by pixel area) over the region of interest. The deseasonalized time series is obtained by applying the X-11 seasonal adjustment methodology on the original time series as described in Colella et al. (2016), and then the Mann-Kendall test (Mann, 1945; Kendall, 1975) and Sens’s method (Sen, 1968) are subsequently applied to obtain the magnitude of trend. This OMI has been introduced since the 2nd issue of Ocean State Report in 2017. '''CONTEXT''' Phytoplankton and 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). Therefore, it is 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 Mediterranean Sea is known to respond to climate variability associated with the North Atlantic Oscillation (NAO) and El Niño Southern Oscillation (ENSO) (Basterretxea et al. 2018, Colella et al. 2016). '''KEY FINDINGS''' In the Mediterranean Sea, the average chlorophyll trend for the 1997–2024 period is slightly negative, at -0.77 ± 0.59% per year, reinforcing the findings of the previous releases. This result contrasts with the analysis by Sathyendranath et al. (2018), which reported increasing chlorophyll concentrations across all European seas. From around 2010–2011 onward, excluding the 2018–2019 period, a noticeable decline in chlorophyll levels is evident in the deseasonalized time series (green line) and in the observed maxima (grey line), particularly from 2015. This sustained decline over the past decade contributes to the overall negative trend observed in the Mediterranean Sea. '''DOI (product):''' https://doi.org/10.48670/moi-00259

'''DEFINITION''' The Iberia-Biscay-Ireland (IBI) Ocean Heat Content (OHC) indicator, OMI_CLIMATE_OHC_IBI_area_averaged_anomalies, provides estimates of OHC anomalies computed over the reference period 1993–2024. The values are integrated from the surface down to 2000 m depth, using a reference density of ρ₀ = 1030 kg·m⁻³ and a specific heat capacity of cₚ = 3980 J·kg⁻¹·°C⁻¹ (e.g., von Schuckmann et al., 2009). This variable is directly proportional to the average temperature change in the ocean. Averaged time series of OHC anomalies and their associated uncertainties are computed for the IBI region (26° N–56° N; 19° W–5° E) using the following Copernicus Marine products: * '''IBI-MYP''' & '''IBI-INT''': IBI_MULTIYEAR_PHY_005_002 (reanalysis and interim datasets) * '''GLO-MYP''': GLOBAL_REANALYSIS_PHY_001_031 (reanalysis) * '''CORA''': INSITU_GLO_PHY_TS_OA_MY_013_052 (in situ observations) * '''ARMOR''': MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012 (reprocessed observations) The figure displays the ensemble mean (blue line) and the ensemble spread (grey shading). Further details on the indicator and data processing are provided in the corresponding Product User Manual (de Pascual-Collar et al., 2026) and in de Pascual-Collar et al. (2023), von Schuckmann et al. (2016), and 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 Meteorological Organization (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; Potter and Lozier, 2004), 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). Among the different periods of interannual variability identified by the indicator, a sustained increase in OHC from 2023 onwards is particularly noteworthy. This short-term trend results in 2024 exhibiting the highest OHC value recorded in the time series. '''CMEMS KEY FINDINGS''' The ensemble mean OHC anomaly time series over the Iberia–Biscay–Ireland region is characterized by marked interannual variability and a a statistically significant ocean warming trend of 0.55 ± 0.3 W m⁻² (99% confidence interval). In addition, the final year of the time series (2024) exhibits the highest OHC value recorded, following a period of sustained warming that began in 2023. '''DOI (product):''' https://doi.org/10.48670/mds-00316

'''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''' '''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''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' Global Ocean- in-situ reprocessed Carbon observations. This product contains observations and gridded files from two up-to-date carbon and biogeochemistry community data products: Surface Ocean Carbon ATlas SOCATv2021 and GLobal Ocean Data Analysis Project GLODAPv2.2021. The SOCATv2021-OBS dataset contains >25 million observations of fugacity of CO2 of the surface global ocean from 1957 to early 2021. The quality control procedures are described in Bakker et al. (2016). These observations form the basis of the gridded products included in SOCATv2020-GRIDDED: monthly, yearly and decadal averages of fCO2 over a 1x1 degree grid over the global ocean, and a 0.25x0.25 degree, monthly average for the coastal ocean. GLODAPv2.2021-OBS contains >1 million observations from individual seawater samples of temperature, salinity, oxygen, nutrients, dissolved inorganic carbon, total alkalinity and pH from 1972 to 2019. These data were subjected to an extensive quality control and bias correction described in Olsen et al. (2020). GLODAPv2-GRIDDED contains global climatologies for temperature, salinity, oxygen, nitrate, phosphate, silicate, dissolved inorganic carbon, total alkalinity and pH over a 1x1 degree horizontal grid and 33 standard depths using the observations from the previous iteration of GLODAP, GLODAPv2. SOCAT and GLODAP are based on community, largely volunteer efforts, and the data providers will appreciate that those who use the data cite the corresponding articles (see References below) in order to support future sustainability of the data products. '''DOI (product) :''' https://doi.org/10.48670/moi-00035

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