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

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' For the Global Ocean- In-situ observation delivered in delayed mode. This In Situ delayed mode product integrates the best available version of in situ oxygen, chlorophyll / fluorescence and nutrients data '''DOI (product) :''' https://doi.org/10.17882/86207

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

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The ocean monitoring indicator on regional mean sea level is derived from the DUACS delayed-time (DT-2021 version) altimeter gridded maps of sea level anomalies based on a stable number of altimeters (two) in the satellite constellation. These products are distributed by the Copernicus Climate Change Service and the Copernicus Marine Service (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). The mean sea level evolution estimated in the Irish-Biscay-Iberian (IBI) region is derived from the average of the gridded sea level maps weighted by the cosine of the latitude. The annual and semi-annual periodic signals are removed (least square fit of sinusoidal function) and the time series is low-pass filtered (175 days cut-off). The curve is corrected for the regional mean effect of the Glacial Isostatic Adjustment (GIA) using the ICE5G-VM2 GIA model (Peltier, 2004). During 1993-1998, the Global men sea level (hereafter GMSL) has been known to be affected by a TOPEX-A instrumental drift (WCRP Global Sea Level Budget Group, 2018; Legeais et al., 2020). This drift led to overestimate the trend of the GMSL during the first 6 years of the altimetry record (about 0.04 mm/y at global scale over the whole altimeter period). A correction of the drift is proposed for the Global mean sea level (Legeais et al., 2020). Whereas this TOPEX-A instrumental drift should also affect the regional mean sea level (hereafter RMSL) trend estimation, currently this empirical correction is currently not applied to the altimeter sea level dataset and resulting estimated for RMSL. Indeed, the pertinence of the global correction applied at regional scale has not been demonstrated yet and there is no clear consensus achieved on the way to proceed at regional scale. Additionally, the estimation of such a correction at regional scale is not obvious, especially in areas where few accurate independent measurements (e.g. in situ)- necessary for this estimation - are available. The trend uncertainty is provided in a 90% confidence interval (Prandi et al., 2021). This estimate only considers errors related to the altimeter observation system (i.e., orbit determination errors, geophysical correction errors and inter-mission bias correction errors). The presence of the interannual signal can strongly influence the trend estimation considering to the altimeter period considered (Wang et al., 2021; Cazenave et al., 2014). The uncertainty linked to this effect is not taken into account. '''CONTEXT''' The indicator on area averaged sea level is a crucial index of climate change, and individual components contribute to sea level rise, including expansion due to ocean warming and melting of glaciers and ice sheets (WCRP Global Sea Level Budget Group, 2018). According to the recent IPCC 6th assessment report, global mean sea level (GMSL) increased by 0.20 (0.15 to 0.25) m over the period 1901 to 2018 with a rate 25 of rise that has accelerated since the 1960s to 3.7 (3.2 to 4.2) mm yr-1 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 RMSL rise can also be 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). Rising sea level can strongly affect population and infrastructures in coastal areas, increase their vulnerability and risks for food security, particularly in low lying areas and island states. Adverse impacts from floods, storms and tropical cyclones with related losses and damages have increased due to sea level rise, and increase their vulnerability and increase risks for food security, particularly in low lying areas and island states (IPCC, 2022b). Adaptation and mitigation measures such as the restoration of mangroves and coastal wetlands, reduce the risks from sea level rise (IPCC, 2022c). In IBI region, the RMSL trend is modulated by decadal variations. As observed over the global ocean, the main actors of the long-term RMSL trend are associated with anthropogenic global/regional warming. Decadal variability is mainly linked to the strengthening or weakening of the Atlantic Meridional Overturning Circulation (AMOC) (e.g. Chafik et al., 2019). The latest is driven by the North Atlantic Oscillation (NAO) (e.g. Delworth and Zeng, 2016). Along the European coast, the NAO also influences the along-slope winds dynamic which in return significantly contributes to the local sea level variability observed (Chafik et al., 2019). '''CMEMS KEY FINDINGS''' Over the [1993/01/01, 2021/08/02] period, the basin-wide RMSL in the IBI area rises at a rate of 3.8  0.82 mm/year. '''DOI (product):''' https://doi.org/10.48670/moi-00252

'''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'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The ibi_omi_tempsal_sst_area_averaged_anomalies product for 2021 includes Sea Surface Temperature (SST) anomalies, given as monthly mean time series starting on 1993 and averaged over the Iberia-Biscay-Irish Seas. The IBI SST OMI is built from the CMEMS Reprocessed European North West Shelf Iberai-Biscay-Irish Seas (SST_MED_SST_L4_REP_OBSERVATIONS_010_026, see e.g. the OMI QUID, http://marine.copernicus.eu/documents/QUID/CMEMS-OMI-QUID-ATL-SST.pdf), which provided the SSTs used to compute the evolution of SST anomalies over the European North West Shelf Seas. This reprocessed product consists of daily (nighttime) interpolated 0.05° grid resolution SST maps over the European North West Shelf Iberai-Biscay-Irish Seas built from the ESA Climate Change Initiative (CCI) (Merchant et al., 2019) and Copernicus Climate Change Service (C3S) initiatives. Anomalies are computed against the 1993-2014 reference period. '''CONTEXT''' Sea surface temperature (SST) is a key climate variable since it deeply contributes in regulating climate and its variability (Deser et al., 2010). SST is then essential to monitor and characterise the state of the global climate system (GCOS 2010). Long-term SST variability, from interannual to (multi-)decadal timescales, provides insight into the slow variations/changes in SST, i.e. the temperature trend (e.g., Pezzulli et al., 2005). In addition, on shorter timescales, SST anomalies become an essential indicator for extreme events, as e.g. marine heatwaves (Hobday et al., 2018). '''CMEMS KEY FINDINGS''' The overall trend in the SST anomalies in this region is 0.011 ±0.001 °C/year over the period 1993-2021. '''DOI (product):''' https://doi.org/10.48670/moi-00256

'''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 SOCATv2022-OBS dataset contains >25 million observations of fugacity of CO2 of the surface global ocean from 1957 to early 2022. 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.2022-OBS contains >1 million observations from individual seawater samples of temperature, salinity, oxygen, nutrients, dissolved inorganic carbon, total alkalinity and pH from 1972 to 2020. 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

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

'''Short description:''' Arctic L3 sea ice product providing concentration, stage-of-development and floe size information retrieved from Sentinel-1 and RCM SAR imagery and GCOM-W AMSR2 microwave radiometer data using a deep learning algorithm and delivered on a 0.5 km grid. '''DOI (product) :''' https://doi.org/10.48670/mds-00343