'''Short description:''' This product consists of vertical profiles of the concentration of nutrients (nitrates, phosphates, and silicates) and carbonate system variables (total alkalinity, dissolved inorganic carbon, pH, and partial pressure of carbon dioxide), computed for each Argo float equipped with an oxygen sensor. The method called CANYON is based on a neural network trained using nutrient data (GLODAPv2 database) '''DOI (product) :''' https://doi.org/10.48670/moi-00048

'''DEFINITION''' Estimates of Arctic sea ice extent are obtained from the surface of oceans grid cells that have at least 15% sea ice concentration. These values are cumulated in the entire Northern Hemisphere (excluding ice lakes) and from 1993 up to the year 2019 aiming to: i) obtain the Arctic sea ice extent as expressed in millions of km square (106 km2) to monitor both the large-scale variability and mean state and change. ii) to monitor the change in sea ice extent as expressed in millions of km squared per decade (106 km2/decade), or in sea ice extent loss since the beginning of the time series as expressed in percent per decade (%/decade; reference period being the first date of the key figure b) dot-dashed trend line, Vaughan et al., 2013). These trends are calculated in three ways, i.e. (i) from the annual mean values; (ii) from the March values (winter ice loss); (iii) from September values (summer ice loss). The Arctic sea ice extent used here is based on the “multi-product” (GLOBAL_MULTIYEAR_PHY_ENS_001_031) approach as introduced in the second issue of the Ocean State Report (CMEMS OSR, 2017). Five global products have been used to build the ensemble mean, and its associated ensemble spread. '''CONTEXT''' Sea ice is frozen seawater that floats on the ocean surface. This large blanket of millions of square kilometers insulates the relatively warm ocean waters from the cold polar atmosphere. The seasonal cycle of the sea ice, forming and melting with the polar seasons, impacts both human activities and biological habitat. Knowing how and how much the sea ice cover is changing is essential for monitoring the health of the Earth as sea ice is one of the highest sensitive natural environments. Variations in sea ice cover can induce changes in ocean stratification, in global and regional sea level rates and modify the key rule played by the cold poles in the Earth engine (IPCC, 2019). The sea ice cover is monitored here in terms of sea ice extent quantity. More details and full scientific evaluations can be found in the CMEMS Ocean State Report (Samuelsen et al., 2016; Samuelsen et al., 2018). '''CMEMS KEY FINDINGS''' Since the year 1993 to 2023 the Arctic sea ice extent has decreased significantly at an annual rate of -0.57*106 km2 per decade. This represents an amount of -4.8 % per decade of Arctic sea ice extent loss over the period 1993 to 2023. Over the period 1993 to 2018, summer (September) sea ice extent loss amounts to -1.18*106 km2/decade (September values), which corresponds to -14.85% per decade. Winter (March) sea ice extent loss amounts to -0.57*106 km2/decade, which corresponds to -3.42% per decade. These values slightly exceed the estimates given in the AR5 IPCC assessment report (estimate up to the year 2012) as a consequence of continuing Northern Hemisphere sea ice extent loss. Main change in the mean seasonal cycle is characterized by less and less presence of sea ice during summertime with time. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00190

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

'''DEFINITION''' The OMI_EXTREME_WAVE_NORTHWESTSHELF_swh_mean_and_anomaly_obs indicator is based on the computation of the 99th and the 1st percentiles from in situ data (observations). It is computed for the variable significant wave height (swh) measured by in situ buoys. The use of percentiles instead of annual maximum and minimum values, makes this extremes study less affected by individual data measurement errors. The percentiles are temporally averaged, and the spatial evolution is displayed, jointly with the anomaly in the target year. This study of extreme variability was first applied to sea level variable (Pérez Gómez et al 2016) and then extended to other essential variables, sea surface temperature and significant wave height (Pérez Gómez et al 2018). '''CONTEXT''' Projections on Climate Change foresee a future with a greater frequency of extreme sea states (Stott, 2016; Mitchell, 2006). The damages caused by severe wave storms can be considerable not only in infrastructure and buildings but also in the natural habitat, crops and ecosystems affected by erosion and flooding aggravated by the extreme wave heights. In addition, wave storms strongly hamper the maritime activities, especially in harbours. These extreme phenomena drive complex hydrodynamic processes, whose understanding is paramount for proper infrastructure management, design and maintenance (Goda, 2010). In recent years, there have been several studies searching possible trends in wave conditions focusing on both mean and extreme values of significant wave height using a multi-source approach with model reanalysis information with high variability in the time coverage, satellite altimeter records covering the last 30 years and in situ buoy measured data since the 1980s decade but with sparse information and gaps in the time series (e.g. Dodet et al., 2020; Timmermans et al., 2020; Young & Ribal, 2019). These studies highlight a remarkable interannual, seasonal and spatial variability of wave conditions and suggest that the possible observed trends are not clearly associated with anthropogenic forcing (Hochet et al. 2021, 2023). In the North Atlantic, the mean wave height shows some weak trends not very statistically significant. Young & Ribal (2019) found a mostly positive weak trend in the European Coasts while Timmermans et al. (2020) showed a weak negative trend in high latitudes, including the North Sea and even more intense in the Norwegian Sea. For extreme values, some authors have found a clearer positive trend in high percentiles (90th-99th) (Young et al., 2011; Young & Ribal, 2019). '''COPERNICUS MARINE SERVICE KEY FINDINGS''' The mean 99th percentiles showed in the area present a wide range from 2.5 meters in the English Channel with 0.3m of standard deviation (std), 3-5m in the southern and central North Sea with 0.3-0.6m of std, 4 meters in the Skagerrak Strait with 0.6m of std, 6-7m in the northern North Sea with 0.4-0.5m of std to 8 meters in the NorthWest of the British Isles with 0.8-1.0m of std. Results for this year show either low positive or negative anomalies between -0.3m and +0.4m, inside the margin of the standard deviation, in the English Channel, the Skagerrak Strait and the southern and central North Sea except in the station 6200046 with a positive anomaly of 0.8m and a slight negative anomaly (-0.1/-0.5m) inside the margin of the std in the NorthWest of the British Isles and the northern North Sea. '''DOI (product):''' https://doi.org/10.48670/moi-00270

'''DEFINITION''' The omi_climate_sst_ibi_area_averaged_anomalies product for 2024 includes Sea Surface Temperature (SST) anomalies, given as monthly mean time series starting on 1982 and averaged over the IBI areas. The IBI SST OMI is built from the CMEMS Reprocessed European North West Shelf Iberai-Biscay-Irish areas (SST_MED_SST_L4_REP_OBSERVATIONS_010_026, see e.g. the OMI QUID, http://marine.copernicus.eu/documents/QUID/CMEMS-OMI-QUID-CLIMATE-SST- IBI_v3.pdf), which provided the SSTs used to compute the evolution of SST anomalies over the IBI areas. This reprocessed product consists of daily (nighttime) interpolated 0.05° grid resolution SST maps over the European North West Shelf Iberai-Biscay-Irish areas built from re-processed ESA SST CCI, C3S (Embury et al., 2019). Anomalies are computed against the 1991-2020 reference period. The reference for this OMI can be found in the first and second issue of the Copernicus Marine Service Ocean State Report (OSR), Section 1.1 (Roquet et al., 2016; Mulet et al., 2018). '''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.012 ±0.002 °C/year over the period 1982-2024. '''DOI (product):''' https://doi.org/10.48670/moi-00256

'''Short description:''' Altimeter satellite gridded Sea Level Anomalies (SLA) computed with respect to a twenty-year [1993, 2012] mean. The SLA is estimated by Optimal Interpolation, merging the L3 along-track measurement from the different altimeter missions available. Part of the processing is fitted to the European Sea area. (see QUID document or http://duacs.cls.fr [http://duacs.cls.fr] pages for processing details). The product gives additional variables (i.e. Absolute Dynamic Topography and geostrophic currents (absolute and anomalies)). It serves in delayed-time applications. This product is processed by the DUACS multimission altimeter data processing system. '''DOI (product):''' https://doi.org/10.48670/moi-00141

'''Short description:''' The Mean Dynamic Topography MDT-CMEMS_2020_MED is an estimate of the mean over the 1993-2012 period of the sea surface height above geoid for the Mediterranean Sea. This is consistent with the reference time period also used in the SSALTO DUACS products '''DOI (product) :''' https://doi.org/10.48670/moi-00151

'''DEFINITION''' The indicator of Volume Transport Anomaly in Selected Vertical Sections in the Iberia–Biscay–Ireland (IBI) region (OMI_CIRCULATION_VOLTRANS_IBI_section_integrated_anomalies) is defined as the time series of annual mean volume transport calculated across a set of vertical ocean sections. These sections have been selected to represent the temporal variability of key ocean currents within the IBI domain. The monitored ocean currents include the transport towards the North Sea through the Rockall Trough (RTE) (Holliday et al., 2008; Lozier and Stewart, 2008), the Canary Current (CC) (Knoll et al., 2002; Mason et al., 2011), the Azores Current (AC) (Mason et al., 2011), the Algerian Current (ALG) (Tintoré et al., 1988; Benzohra and Millot, 1995; Font et al., 1998), and the net transport along the 48° N latitude parallel (N48) (see OMI figure). To produce ensemble-based results, six datasets provided by the Copernicus Marine Service have been used: * '''IBI-REA''' & '''IBI-INT''': IBI_MULTIYEAR_PHY_005_002 (reanalysis and interim datasets) * '''GLO-REA''': GLOBAL_MULTIYEAR_PHY_001_030 (reanalysis) * '''ARMOR''': MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012 (reprocessed observations) * '''MED-REA''': MEDSEA_MULTIYEAR_PHY_006_004 (reanalysis) * '''NWS-REA''': NWSHELF_MULTIYEAR_PHY_004_009 (reanalysis) The time series displays the ensemble mean (blue line), the ensemble spread (grey shaded area), and the mean transport with reversed sign (red dashed line), which indicates the threshold of anomaly values corresponding to a reversal in the direction of the current transport. In addition, the trend analysis at the 95% confidence level is shown in the bottom-right corner of each diagram. Further details on the product are provided in the corresponding Product User Manual (de Pascual-Collar et al., 2026a) and Quality Information Document (de Pascual-Collar et al., 2026b), as well as in de Pascual-Collar et al., 2024. '''CONTEXT''' The IBI area is a highly complex region characterized by a remarkable variety of ocean currents. Among them, we can highlight those that originate as a result of the closure of the North Atlantic Drift (Mason et al., 2011; Holliday et al., 2008; Peliz et al., 2007; Bower et al., 2002; Knoll et al., 2002; Pérez et al., 2001; Jia, 2000); the subsurface currents flowing northward along the continental slope (de Pascual-Collar et al., 2019; Pascual et al., 2018; Machín et al., 2010; Fricourt et al., 2007; Knoll et al., 2002; Mazé et al., 1997; White & Bowyer, 1997); and the exchange currents occurring in the Strait of Gibraltar and the Alboran Sea (Sotillo et al., 2016; Font et al., 1998; Benzohra & Millot, 1995; Tintoré et al., 1988). The variability of ocean currents in the IBI domain is relevant to the global thermohaline circulation and other climatic and environmental processes. For example, as discussed by Fasullo and Trenberth (2008), subtropical gyres play a crucial role in the meridional energy balance. The poleward salt transport of Mediterranean water, driven by subsurface slope currents, has significant implications for salinity anomalies in the Rockall Trough and the Nordic Seas, as studied by Holliday (2003), Holliday et al. (2008), and Bozec et al. (2011). The Algerian Current serves as the only pathway for Atlantic Water to reach the Western Mediterranean. '''CMEMS KEY FINDINGS''' The volume transport time series reveal periods during which the monitored currents exhibited notably high or low variability. Specifically, the RTE current shows pronounced variability in 2010 and during 2014–2015; the N48 section between 2012 and 2014; the ALG current in 2006 and 2017; the AC current between 2005–2007 and in 2021; and the CC current between 2005–2007. Furthermore, certain periods display anomalies of sufficient magnitude (in absolute value) to indicate a reversal in the net transport direction of the current. This is the case for the ALG current in 2017 and 2024 (with net transport towards the west), and for the CC current in 2010 (with net transport towards the north). Trend analysis over the period 1993–2023 does not reveal any statistically significant trends for the monitored currents. However, the confidence interval for the trend in the ALG section is close to rejecting the null hypothesis of no trend. '''DOI (product):''' https://doi.org/10.48670/mds-00351

'''DEFINITION''' The Copernicus Marine IBI_OMI_seastate_extreme_var_swh_mean_and_anomaly OMI indicator is based on the computation of the annual 99th percentile of Significant Wave Height (SWH) from model data. Two different CMEMS products are used to compute the indicator: The Iberia-Biscay-Ireland Multi Year Product (IBI_MULTIYEAR_WAV_005_006) and the Analysis product (IBI_ANALYSISFORECAST_WAV_005_005). 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 in the whole period (1980-2023). * Anomaly of the 99th percentile in 2024: The 99th percentile of the year 2024 is computed from the Analysis product. The anomaly is obtained by subtracting the mean percentile to the percentile in 2024. This indicator is aimed at monitoring the extremes of annual significant wave height and evaluate the spatio-temporal variability. The use of percentiles instead of annual maxima, makes this extremes study less affected by individual data. This approach was first successfully applied to 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 Álvarez-Fanjul et al., 2019). Further details and in-depth scientific evaluation can be found in the CMEMS Ocean State report (Álvarez- Fanjul et al., 2019). '''CONTEXT''' The sea state and its related spatio-temporal variability affect dramatically maritime activities and the physical connectivity between offshore waters and coastal ecosystems, impacting therefore on the biodiversity of marine protected areas (González-Marco et al., 2008; Savina et al., 2003; Hewitt, 2003). Over the last decades, significant attention has been devoted to extreme wave height events since their destructive effects in both the shoreline environment and human infrastructures have prompted a wide range of adaptation strategies to deal with natural hazards in coastal areas (Hansom et al., 2015). Complementarily, there is also an emerging question about the role of anthropogenic global climate change on present and future extreme wave conditions (Young and Ribal, 2019). The Iberia-Biscay-Ireland region, which covers the North-East Atlantic Ocean from Canary Islands to Ireland, is characterized by two different sea state wave climate regions: whereas the northern half, impacted by the North Atlantic subpolar front, is of one of the world’s greatest wave generating regions (Mørk et al., 2010; Folley, 2017), the southern half, located at subtropical latitudes, is by contrast influenced by persistent trade winds and thus by constant and moderate wave regimes. The North Atlantic Oscillation (NAO), which refers to changes in the atmospheric sea level pressure difference between the Azores and Iceland, is a significant driver of wave climate variability in the Northern Hemisphere. The influence of North Atlantic Oscillation on waves along the Atlantic coast of Europe is particularly strong in and has a major impact on northern latitudes wintertime (Gleeson et al., 2017; Martínez-Asensio et al. 2016; Wolf et al., 2002; Bauer, 2001; Kushnir et al., 1997; Bouws et al., 1996; Bacon and Carter, 1991). Swings in the North Atlantic Oscillation index produce changes in the storms track and subsequently in the wind speed and direction over the Atlantic that alter the wave regime. When North Atlantic Oscillation index is in its positive phase, storms usually track northeast of Europe and enhanced westerly winds induce higher than average waves in the northernmost Atlantic Ocean. Conversely, in the negative North Atlantic Oscillation phase, the track of the storms is more zonal and south than usual, with trade winds (mid latitude westerlies) being slower and producing higher than average waves in southern latitudes (Marshall et al., 2001; Wolf et al., 2002; Wolf and Woolf, 2006). Additionally, a variety of previous studies have uniquevocally determined the relationship between the sea state variability in the IBI region and other atmospheric climate modes such as the East Atlantic pattern, the Arctic Oscillation, the East Atlantic Western Russian pattern and the Scandinavian pattern (Izaguirre et al., 2011, Martínez-Asensio et al., 2016). In this context, long‐term statistical analysis of reanalyzed model data is mandatory not only to disentangle other driving agents of wave climate but also to attempt inferring any potential trend in the number and/or intensity of extreme wave events in coastal areas with subsequent socio-economic and environmental consequences. '''CMEMS KEY FINDINGS''' The climatic mean of 99th percentile (1980-2023) reveals a north-south gradient of Significant Wave Height with the highest values in northern latitudes (above 8m) and lowest values (2-3 m) detected southeastward of Canary Islands, in the seas between Canary Islands and the African Continental Shelf. This north-south pattern is the result of the two climatic conditions prevailing in the region and previously described. The 99th percentile anomalies in 2024 show that during this period, virtually the entire IBI region was affected by positive anomalies in maximum SWH, which exceeded the standard deviation of the historical record in the waters west of the Iberian Peninsula, the Spanish coast of the Bay of Biscay, and the African coast south of Cape Ghir. Anomalies reaching twice the standard deviation of the time series were also observed in coastal regions of the English Channel. '''DOI (product):''' https://doi.org/10.48670/moi-00249