'''Short description:''' For the Atlantic Ocean - The product contains daily Level-3 sea surface wind with a 1km horizontal pixel spacing using Synthetic Aperture Radar (SAR) observations and their collocated European Centre for Medium-Range Weather Forecasts (ECMWF) model outputs. Products are processed homogeneously starting from the L2OCN products. '''DOI (product) :''' https://doi.org/10.48670/mds-00339

'''This product has been archived''' '''DEFINITION''' The CMEMS IBI_OMI_tempsal_extreme_var_temp_mean_and_anomaly OMI indicator is based on the computation of the annual 99th percentile of Sea Surface Temperature (SST) from model data. Two different CMEMS products are used to compute the indicator: The Iberia-Biscay-Ireland Multi Year Product (IBI_MULTIYEAR_PHY_005_002) and the Analysis product (IBI_ANALYSISFORECAST_PHY_005_001). Two parameters have been considered for this OMI: • Map of the 99th mean percentile: It is obtained from the Multi Year Product, the annual 99th percentile is computed for each year of the product. The percentiles are temporally averaged over the whole period (1993-2021). • Anomaly of the 99th percentile in 2022: The 99th percentile of the year 2022 is computed from the Analysis product. The anomaly is obtained by subtracting the mean percentile from the 2022 percentile. This indicator is aimed at monitoring the extremes of sea surface temperature every year and at checking their variations in space. The use of percentiles instead of annual maxima, makes this extremes study less affected by individual data. This study of extreme variability was first applied to the sea level variable (Pérez Gómez et al 2016) and then extended to other essential variables, such as sea surface temperature and significant wave height (Pérez Gómez et al 2018 and Alvarez Fanjul et al., 2019). More details and a full scientific evaluation can be found in the CMEMS Ocean State report (Alvarez Fanjul et al., 2019). '''CONTEXT''' The Sea Surface Temperature is one of the essential ocean variables, hence the monitoring of this variable is of key importance, since its variations can affect the ocean circulation, marine ecosystems, and ocean-atmosphere exchange processes. As the oceans continuously interact with the atmosphere, trends of sea surface temperature can also have an effect on the global climate. While the global-averaged sea surface temperatures have increased since the beginning of the 20th century (Hartmann et al., 2013) in the North Atlantic, anomalous cold conditions have also been reported since 2014 (Mulet et al., 2018; Dubois et al., 2018). The IBI area is a complex dynamic region with a remarkable variety of ocean physical processes and scales involved. The Sea Surface Temperature field in the region is strongly dependent on latitude, with higher values towards the South (Locarnini et al. 2013). This latitudinal gradient is supported by the presence of the eastern part of the North Atlantic subtropical gyre that transports cool water from the northern latitudes towards the equator. Additionally, the Iberia-Biscay-Ireland region is under the influence of the Sea Level Pressure dipole established between the Icelandic low and the Bermuda high. Therefore, the interannual and interdecadal variability of the surface temperature field may be influenced by the North Atlantic Oscillation pattern (Czaja and Frankignoul, 2002; Flatau et al., 2003). Also relevant in the region are the upwelling processes taking place in the coastal margins. The most referenced one is the eastern boundary coastal upwelling system off the African and western Iberian coast (Sotillo et al., 2016), although other smaller upwelling systems have also been described in the northern coast of the Iberian Peninsula (Alvarez et al., 2011), the south-western Irish coast (Edwars et al., 1996) and the European Continental Slope (Dickson, 1980). '''CMEMS KEY FINDINGS''' In the IBI region, the 99th mean percentile for 1993-2021 shows a north-south pattern driven by the climatological distribution of temperatures in the North Atlantic. In the coastal regions of Africa and the Iberian Peninsula, the mean values are influenced by the upwelling processes (Sotillo et al., 2016). These results are consistent with the ones presented in Álvarez Fanjul (2019) for the period 1993-2016. The analysis of the 99th percentile anomaly in the year 2023 shows that this period has been affected by a severe impact of maximum SST values. Anomalies exceeding the standard deviation affect almost the entire IBI domain, and regions impacted by thermal anomalies surpassing twice the standard deviation are also widespread below the 43ºN parallel. Extreme SST values exceeding twice the standard deviation affect not only the open ocean waters but also the easter boundary upwelling areas such as the northern half of Portugal, the Spanish Atlantic coast up to Cape Ortegal, and the African coast south of Cape Aguer. It is worth noting the impact of anomalies that exceed twice the standard deviation is widespread throughout the entire Mediterranean region included in this analysis. '''DOI (product):''' https://doi.org/10.48670/moi-00254

'''Short description: ''' For the '''Atlantic''' Ocean '''Satellite Observations''', ACRI-ST company (Sophia Antipolis, France) is providing '''Bio-Geo-Chemical (BGC)''' products based on the '''Copernicus-GlobColour''' processor. * Upstreams: SeaWiFS, MODIS, MERIS, VIIRS-SNPP & JPSS1, OLCI-S3A & S3B for the '''""multi""''' products, and S3A & S3B only for the '''""olci""''' products. * Variables: Chlorophyll-a ('''CHL'''), Gradient of Chlorophyll-a ('''CHL_gradient'''), Phytoplankton Functional types and sizes ('''PFT'''), Suspended Matter ('''SPM'''), Secchi Transparency Depth ('''ZSD'''), Diffuse Attenuation ('''KD490'''), Particulate Backscattering ('''BBP'''), Absorption Coef. ('''CDM''') and Reflectance ('''RRS'''). * Temporal resolutions: '''daily'''. * Spatial resolutions: '''1 km''' and a finer resolution based on olci '''300 meters''' inputs. * Recent products are organized in datasets called Near Real Time ('''NRT''') and long time-series (from 1997) in datasets called Multi-Years ('''MY'''). To find the '''Copernicus-GlobColour''' products in the catalogue, use the search keyword '''""GlobColour""'''. '''DOI (product) :''' https://doi.org/10.48670/moi-00286

'''DEFINITION''' Based on daily, global climate sea surface temperature (SST) analyses generated by the Copernicus Climate Change Service (C3S) (product SST-GLO-SST-L4-REP-OBSERVATIONS-010-024). Analysis of the data was based on the approach described in Mulet et al. (2018) and is described and discussed in Good et al. (2020). The processing steps applied were: 1. The daily analyses were averaged to create monthly means. 2. A climatology was calculated by averaging the monthly means over the period 1991 - 2020. 3. Monthly anomalies were calculated by differencing the monthly means and the climatology. 4. The time series for each grid cell was passed through the X11 seasonal adjustment procedure, which decomposes a time series into a residual seasonal component, a trend component and errors (e.g., Pezzulli et al., 2005). The trend component is a filtered version of the monthly time series. 5. The slope of the trend component was calculated using a robust method (Sen 1968). The method also calculates the 95% confidence range in the slope. '''CONTEXT''' Sea surface temperature (SST) is one of the Essential Climate Variables (ECVs) defined by the Global Climate Observing System (GCOS) as being needed for monitoring and characterising the state of the global climate system (GCOS 2010). It provides insight into the flow of heat into and out of the ocean, into modes of variability in the ocean and atmosphere, can be used to identify features in the ocean such as fronts and upwelling, and knowledge of SST is also required for applications such as ocean and weather prediction (Roquet et al., 2016). '''CMEMS KEY FINDINGS''' Warming trends occurred over most of the globe between 1982 and 2024, with the strongest warming in the Northern Pacific and Atlantic Oceans. However, there were cooling trends in parts of the Southern Ocean and the South-East Pacific Ocean. '''DOI (product):''' https://doi.org/10.48670/moi-00243

'''DEFINITION''' Ocean heat content (OHC) is defined here as the deviation from a reference period (1993-2014) and is closely proportional to the average temperature change from z1 = 0 m to z2 = 700 m depth: OHC=∫_(z_1)^(z_2)ρ_0 c_p (T_yr-T_clim )dz [1] with a reference density of = 1030 kgm-3 and a specific heat capacity of cp = 3980 J kg-1 °C-1 (e.g. von Schuckmann et al., 2009). Time series of annual mean values area averaged ocean heat content is provided for the Mediterranean Sea (30°N, 46°N; 6°W, 36°E) and is evaluated for topography deeper than 300m. '''CONTEXT''' Knowing how much and where heat energy is stored and released in the ocean is essential for understanding the contemporary Earth system state, variability and change, as the oceans shape our perspectives for the future. The quality evaluation of MEDSEA_OMI_OHC_area_averaged_anomalies is based on the “multi-product” approach as introduced in the second issue of the Ocean State Report (von Schuckmann et al., 2018), and following the MyOcean’s experience (Masina et al., 2017). Six global products and a regional (Mediterranean Sea) product have been used to build an ensemble mean, and its associated ensemble spread. The reference products are: • The Mediterranean Sea Reanalysis at 1/24 degree horizontal resolution (MEDSEA_MULTIYEAR_PHY_006_004, DOI: https://doi.org/10.25423/CMCC/MEDSEA_MULTIYEAR_PHY_006_004_E3R1, Escudier et al., 2020) • Four global reanalyses at 1/4 degree horizontal resolution (GLOBAL_MULTIYEAR_PHY_ENS_001_031): GLORYS, C-GLORS, ORAS5, FOAM • Two observation based products: CORA (INSITU_GLO_PHY_TS_OA_MY_013_052) and ARMOR3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012). Details on the products are delivered in the PUM and QUID of this OMI. '''CMEMS KEY FINDINGS''' The ensemble mean ocean heat content anomaly time series over the Mediterranean Sea shows a continuous increase in the period 1993-2022 at rate of 1.38±0.08 W/m2 in the upper 700m. After 2005 the rate has clearly increased with respect the previous decade, in agreement with Iona et al. (2018). '''DOI (product):''' https://doi.org/10.48670/moi-00261

'''This product has been archived''' '''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''' The global yearly ocean CO2 sink represents the ocean uptake of CO2 from the atmosphere computed over the whole ocean. It is expressed in PgC per year. The ocean monitoring index is presented for the period 1985 to year-1. The yearly estimate of the ocean CO2 sink corresponds to the mean of a 100-member ensemble of CO2 flux estimates (Chau et al. 2022). The range of an estimate with the associated uncertainty is then defined by the empirical 68% interval computed from the ensemble. '''CONTEXT''' Since the onset of the industrial era in 1750, the atmospheric CO2 concentration has increased from about 277±3 ppm (Joos and Spahni, 2008) to 412.44±0.1 ppm in 2020 (Dlugokencky and Tans, 2020). By 2011, the ocean had absorbed approximately 28 ± 5% of all anthropogenic CO2 emissions, thus providing negative feedback to global warming and climate change (Ciais et al., 2013). The ocean CO2 sink is evaluated every year as part of the Global Carbon Budget (Friedlingstein et al. 2022). The uptake of CO2 occurs primarily in response to increasing atmospheric levels. The global flux is characterized by a significant variability on interannual to decadal time scales largely in response to natural climate variability (e.g., ENSO) (Friedlingstein et al. 2022, Chau et al. 2022). '''CMEMS KEY FINDINGS''' The rate of change of the integrated yearly surface downward flux has increased by 0.04±0.01e-1 PgC/yr2 over the period 1985 to year-1. The yearly flux time series shows a plateau in the 90s followed by an increase since 2000 with a growth rate of 0.06±0.04e-1 PgC/yr2. In 2021 (resp. 2020), the global ocean CO2 sink was 2.41±0.13 (resp. 2.50±0.12) PgC/yr. The average over the full period is 1.61±0.10 PgC/yr with an interannual variability (temporal standard deviation) of 0.46 PgC/yr. In order to compare these fluxes to Friedlingstein et al. (2022), the estimate of preindustrial outgassing of riverine carbon of 0.61 PgC/yr, which is in between the estimate by Jacobson et al. (2007) (0.45±0.18 PgC/yr) and the one by Resplandy et al. (2018) (0.78±0.41 PgC/yr) needs to be added. A full discussion regarding this OMI can be found in section 2.10 of the Ocean State Report 4 (Gehlen et al., 2020) and in Chau et al. (2022). '''DOI (product):''' https://doi.org/10.48670/moi-00223

'''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) and 5Hz (~1km) sampling. It serves in near-real 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, HY-2B). The system exploits the most recent datasets available based on the enhanced OGDR/NRT+IGDR/STC production. All the missions are homogenized with respect to a reference mission. Part of the processing is fitted to the European Seas. (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 http://marine.copernicus.eu/services-portfolio/access-to-products/?option=com_csw&view=details&product_id=SEALEVEL_GLO_PHY_NOISE_L4_STATIC_008_033 [http://marine.copernicus.eu/services-portfolio/access-to-products/?option=com_csw&view=details&product_id=SEALEVEL_GLO_PHY_NOISE_L4_STATIC_008_033] 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-00140

'''Short description:''' Near-Real-Time multi-mission global satellite-based spectral integral parameters. Only valid data are used, based on the L3 corresponding products. Included wave parameters are partition significant wave height, partition peak period and partition peak or principal direction. Those parameters are propagated in space and time at a 3-hour timestep and on a regular space grid, providing information of the swell propagation characteristics, from source to land. The ouput products corresponds to one file per month gathering all the swell systems at a global scale. This product is processed by the WAVE-TAC multi-mission SAR and CFOSAT/SWIM data processing system to serve in near-real time the main operational oceanography and climate forecasting centers in Europe and worldwide. It processes data from the following missions: SAR (Sentinel-1A and Sentinel-1B) and CFOSAT/SWIM. All the spectral parameter measurements are optimally interpolated using swell observations belonging to the same swell field. The spectral data processing system produces wave integral parameters by partition (partition significant wave height, partition peak period and partition peak or principal direction) and the associated standard deviation and density of propagated observations. '''DOI (product) :''' https://doi.org/10.48670/moi-00175

'''This product has been archived''' '''Short description:''' Near-Real-Time multi-mission global satellite-based spectral integral parameters. Only valid data are used, based on the L3 corresponding product. Included wave parameters are partition significant wave height, partition peak period and partition peak or principal direction. Those parameters are propagated in space and time at a 3-hour timestep and on a regular space grid, providing information of the swell propagation characteristics, from source to land. One file gathers one swell system, gathering observations originating from the same storm source. This product is processed by the WAVE-TAC multi-mission SAR data processing system to serve in near-real time the main operational oceanography and climate forecasting centers in Europe and worldwide. It processes data from the following SAR missions: Sentinel-1A and Sentinel-1B. All the spectral parameter measurements are optimally interpolated using swell observations belonging to the same swell field. The SAR data processing system produces wave integral parameters by partition (partition significant wave height, partition peak period and partition peak or principal direction) and the associated standard deviation and density of propagated observations. '''DOI (product) :''' https://doi.org/10.48670/moi-00175