'''DEFINITION''' Heat transport across lines are obtained by integrating the heat fluxes along some selected sections and from top to bottom of the ocean. The values are computed from models’ daily output. The mean value over a reference period (1993-2014) and over the last full year are provided for the ensemble product and the individual reanalysis, as well as the standard deviation for the ensemble product over the reference period (1993-2014). The values are given in PetaWatt (PW). '''CONTEXT''' The ocean transports heat and mass by vertical overturning and horizontal circulation, and is one of the fundamental dynamic components of the Earth’s energy budget (IPCC, 2013). There are spatial asymmetries in the energy budget resulting from the Earth’s orientation to the sun and the meridional variation in absorbed radiation which support a transfer of energy from the tropics towards the poles. However, there are spatial variations in the loss of heat by the ocean through sensible and latent heat fluxes, as well as differences in ocean basin geometry and current systems. These complexities support a pattern of oceanic heat transport that is not strictly from lower to high latitudes. Moreover, it is not stationary and we are only beginning to unravel its variability. '''CMEMS KEY FINDINGS''' The mean transports estimated by the ensemble global reanalysis are comparable to estimates based on observations; the uncertainties on these integrated quantities are still large in all the available products. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00245

'''DEFINITION''' The omi_climate_sst_ibi_trend product includes the Sea Surface Temperature (SST) trend for the Iberia-Biscay-Irish areas over the period 1982-2023, i.e. the rate of change (°C/year). This OMI is derived from the CMEMS REP ATL L4 SST product (SST_ATL_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 SST trend over the Iberia-Biscay-Irish areas. This reprocessed product consists of daily (nighttime) interpolated 0.05° grid resolution SST maps built from the ESA Climate Change Initiative (CCI) (Merchant et al., 2019) and Copernicus Climate Change Service (C3S) initiatives. Trend analysis has been performed by using the X-11 seasonal adjustment procedure (see e.g. Pezzulli et al., 2005), which has the effect of filtering the input SST time series acting as a low bandpass filter for interannual variations. Mann-Kendall test and Sens’s method (Sen 1968) were applied to assess whether there was a monotonic upward or downward trend and to estimate the slope of the trend and its 95% confidence interval. '''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.022 ±0.002 °C/year over the period 1982-2023. '''Figure caption''' Sea surface temperature trend over the period 1982-2023 in the Iberia-Biscay-Irish areas. The trend is the rate of change (°C/year). The trend map in sea surface temperature is derived from the CMEMS SST_ATL_SST_L4_REP_OBSERVATIONS_010_026 product (see e.g. the OMI QUID, http://marine.copernicus.eu/documents/QUID/CMEMS-OMI-QUID-ATL-SST.pdf). The trend is estimated by using the X-11 seasonal adjustment procedure (e.g. Pezzulli et al., 2005) and Sen’s method (Sen 1968). '''DOI (product):''' https://doi.org/10.48670/moi-00257

'''DEFINITION''' The indicator of the Kuroshio extension phase variations is based on the standardized high frequency altimeter Eddy Kinetic Energy (EKE) averaged in the area 142-149°E and 32-37°N and computed from the DUACS (https://duacs.cls.fr) delayed-time (reprocessed version DT-2021, CMEMS SEALEVEL_GLO_PHY_L4_MY_008_047, including “my” (multi-year) & “myint” (multi-year interim) datasets) and near real-time (CMEMS SEALEVEL_GLO_PHY_L4_NRT _008_046) altimeter sea level gridded products. The change in the reprocessed version (previously DT-2018) and the extension of the mean value of the EKE (now 27 years, previously 20 years) induce some slight changes not impacting the general variability of the Kuroshio extension (correlation coefficient of 0.988 for the total period, 0.994 for the delayed time period only). '''CONTEXT''' The Kuroshio Extension is an eastward-flowing current in the subtropical western North Pacific after the Kuroshio separates from the coast of Japan at 35°N, 140°E. Being the extension of a wind-driven western boundary current, the Kuroshio Extension is characterized by a strong variability and is rich in large-amplitude meanders and energetic eddies (Niiler et al., 2003; Qiu, 2003, 2002). The Kuroshio Extension region has the largest sea surface height variability on sub-annual and decadal time scales in the extratropical North Pacific Ocean (Jayne et al., 2009; Qiu and Chen, 2010, 2005). Prediction and monitoring of the path of the Kuroshio are of huge importance for local economies as the position of the Kuroshio extension strongly determines the regions where phytoplankton and hence fish are located. Unstable (contracted) phase of the Kuroshio enhance the production of Chlorophyll (Lin et al., 2014). '''CMEMS KEY FINDINGS''' The different states of the Kuroshio extension phase have been presented and validated by (Bessières et al., 2013) and further reported by Drévillon et al. (2018) in the Copernicus Ocean State Report #2. Two rather different states of the Kuroshio extension are observed: an ‘elongated state’ (also called ‘strong state’) corresponding to a narrow strong steady jet, and a ‘contracted state’ (also called ‘weak state’) in which the jet is weaker and more unsteady, spreading on a wider latitudinal band. When the Kuroshio Extension jet is in a contracted (elongated) state, the upstream Kuroshio Extension path tends to become more (less) variable and regional eddy kinetic energy level tends to be higher (lower). In between these two opposite phases, the Kuroshio extension jet has many intermediate states of transition and presents either progressively weakening or strengthening trends. In 2018, the indicator reveals an elongated state followed by a weakening neutral phase since then. '''Figure caption''' Standardized Eddy Kinetic Energy over the Kuroshio region (following Bessières et al., 2013) Blue shaded areas correspond to well established strong elongated states periods, while orange shaded areas fit weak contracted states periods. The ocean monitoring indicator is derived from the DUACS delayed-time (reprocessed version DT-2021, “my” (multi-year) dataset used when available, “myint” (multi-year interim) used after) completed by DUACS near Real Time (“nrt”) sea level multi-mission gridded products. The vertical red line shows the date of the transition between “myint” and “nrt” products used. '''DOI (product):''' https://doi.org/10.48670/moi-00222

'''DEFINITION''' The OMI_EXTREME_SST_MEDSEA_sst_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 sea surface temperature measured by in situ buoys at depths between 0 and 5 meters. 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''' Sea surface temperature (SST) is one of the essential ocean variables affected by climate change (mean SST trends, SST spatial and interannual variability, and extreme events). In Europe, several studies show warming trends in mean SST for the last years (von Schuckmann et al., 2016; IPCC, 2021, 2022). An exception seems to be the North Atlantic, where, in contrast, anomalous cold conditions have been observed since 2014 (Mulet et al., 2018; Dubois et al. 2018; IPCC 2021, 2022). Extremes may have a stronger direct influence in population dynamics and biodiversity. According to Alexander et al. 2018 the observed warming trend will continue during the 21st Century and this can result in exceptionally large warm extremes. Monitoring the evolution of sea surface temperature extremes is, therefore, crucial. The Mediterranean Sea has showed a constant increase of the SST in the last three decades across the whole basin with more frequent and severe heat waves (Juza et al., 2022). Deep analyses of the variations have displayed a non-uniform rate in space, being the warming trend more evident in the eastern Mediterranean Sea with respect to the western side. This variation rate is also changing in time over the three decades with differences between the seasons (e.g. Pastor et al. 2018; Pisano et al. 2020), being higher in Spring and Summer, which would affect the extreme values. '''COPERNICUS MARINE SERVICE KEY FINDINGS''' The mean 99th percentiles showed in the area present values from 25ºC in Ionian Sea and 26º in the Alboran sea and Gulf of Lion to 27ºC in the East of Iberian Peninsula. The standard deviation ranges from 0.6ºC to 1.2ºC in the Western Mediterranean and is around 2.2ºC in the Ionian Sea. Results for this year show a slight negative anomaly in the Ionian Sea (-1ºC) inside the standard deviation and a clear positive anomaly in the Western Mediterranean Sea reaching +2.2ºC, almost two times the standard deviation in the area. '''DOI (product):''' https://doi.org/10.48670/moi-00267

'''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.03e-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''' Ocean acidification is quantified by decreases in pH, which is a measure of acidity: a decrease in pH value means an increase in acidity, that is, acidification. The observed decrease in ocean pH resulting from increasing concentrations of CO2 is an important indicator of global change. The estimate of global mean pH builds on a reconstruction methodology, * Obtain values for alkalinity based on the so called “locally interpolated alkalinity regression (LIAR)” method after Carter et al., 2016; 2018. * Build on surface ocean partial pressure of carbon dioxide (CMEMS product: MULTIOBS_GLO_BIO_CARBON_SURFACE_REP_015_008) obtained from an ensemble of Feed-Forward Neural Networks (Chau et al. 2022) which exploit sampling data gathered in the Surface Ocean CO2 Atlas (SOCAT) (https://www.socat.info/) * Derive a gridded field of ocean surface pH based on the van Heuven et al., (2011) CO2 system calculations using reconstructed pCO2 (MULTIOBS_GLO_BIO_CARBON_SURFACE_REP_015_008) and alkalinity. The global mean average of pH at yearly time steps is then calculated from the gridded ocean surface pH field. It is expressed in pH unit on total hydrogen ion scale. In the figure, the amplitude of the uncertainty (1σ ) of yearly mean surface sea water pH varies at a range of (0.0023, 0.0029) pH unit (see Quality Information Document for more details). The trend and uncertainty estimates amount to -0.0017±0.0004e-1 pH units per year. The indicator is derived from in situ observations of CO2 fugacity (SOCAT data base, www.socat.info, Bakker et al., 2016). These observations are still sparse in space and time. Monitoring pH at higher space and time resolutions, as well as in coastal regions will require a denser network of observations and preferably direct pH measurements. A full discussion regarding this OMI can be found in section 2.10 of the Ocean State Report 4 (Gehlen et al., 2020). '''CONTEXT''' The decrease in surface ocean pH is a direct consequence of the uptake by the ocean of carbon dioxide. It is referred to as ocean acidification. The International Panel on Climate Change (IPCC) Workshop on Impacts of Ocean Acidification on Marine Biology and Ecosystems (2011) defined Ocean Acidification as “a reduction in the pH of the ocean over an extended period, typically decades or longer, which is caused primarily by uptake of carbon dioxide from the atmosphere, but can also be caused by other chemical additions or subtractions from the ocean”. The pH of contemporary surface ocean waters is already 0.1 lower than at pre-industrial times and an additional decrease by 0.33 pH units is projected over the 21st century in response to the high concentration pathway RCP8.5 (Bopp et al., 2013). Ocean acidification will put marine ecosystems at risk (e.g. Orr et al., 2005; Gehlen et al., 2011; Kroeker et al., 2013). The monitoring of surface ocean pH has become a focus of many international scientific initiatives (http://goa-on.org/) and constitutes one target for SDG14 (https://sustainabledevelopment.un.org/sdg14). '''CMEMS KEY FINDINGS''' Since the year 1985, global ocean surface pH is decreasing at a rate of -0.0017±0.0004e-1 per year. '''DOI (product):''' https://doi.org/10.48670/moi-00224

'''DEFINITION''' The OMI_EXTREME_SL_MEDSEA_slev_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 sea level measured by tide gauges along the coast. The use of percentiles instead of annual maximum and minimum values, makes this extremes study less affected by individual data measurement errors. The annual percentiles referred to annual mean sea level are temporally averaged and their spatial evolution is displayed in the dataset omi_extreme_sl_medsea_slev_mean_and_anomaly_obs, 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''' Sea level (SLEV) is one of the Essential Ocean Variables most affected by climate change. Global mean sea level rise has accelerated since the 1990’s (Abram et al., 2019, Legeais et al., 2020), due to the increase of ocean temperature and mass volume caused by land ice melting (WCRP, 2018). Basin scale oceanographic and meteorological features lead to regional variations of this trend that combined with changes in the frequency and intensity of storms could also rise extreme sea levels up to one meter by the end of the century (Vousdoukas et al., 2020, Tebaldi et al., 2021). This will significantly increase coastal vulnerability to storms, with important consequences on the extent of flooding events, coastal erosion and damage to infrastructures caused by waves (Boumis et al., 2023). The increase in extreme sea levels over recent decades is, therefore, primarily due to the rise in mean sea level. Note, however, that the methodology used to compute this OMI removes the annual 50th percentile, thereby discarding the mean sea level trend to isolate changes in storminess. The Mediterranean Sea shows statistically significant positive sea level trends over the whole basin. However, at sub-basin scale sea level trends show spatial variability arising from local circulation (Calafat et al., 2022; Meli et al., 2023). '''COPERNICUS MARINE SERVICE KEY FINDINGS''' The completeness index criteria is fulfilled in this region by 38 stations, 26 more than in 2021, significantly increasing spatial coverage with new in situ data in the central Mediterranean Sea, primarily from Italian stations. The mean 99th percentiles reflect the spatial variability of the tide, a microtidal regime, along the Spanish, French and Italian coasts, ranging from around 0.20 m above mean sea level in Sicily and the Balearic Islands (e.g.: 0.22 m in Porto Empedocle, 0.23 m in Ibiza)) to around 0.60 m above mean sea level in the Northern Adriatic Sea (e.g.: 0.63 m in Trieste, 0.61 m in Venice). . The annual 99th percentiles standard deviation ranges between 2 cm in Málaga and Motril (South of Spain) to 8 cm in Marseille. . The 2022 99th percentile anomalies present negative values mainly along the Spanish coast (as in 2021) and in the islands of Corsica and Sardinia (Western part of the region), while positive values are observed along the Eastern French Mediterranean coast and at most of the Italian stations (closer to the central part of the region), with values ranging from -4 cm in Málaga and Motril (Spain) to +5 cm in Ancona (Italy). '''DOI (product):''' https://doi.org/10.48670/moi-00265

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The BALTIC_OMI_TEMPSAL_sst_area_averaged_anomalies product includes time series of monthly mean SST anomalies over the period 1993-2021, relative to the 1993-2014 climatology, averaged for the Baltic Sea. The OMI time series runs from Jan 1, 1993 to December 31, 2021 and is constructed by calculating monthly averages from the daily level 4 SST analysis fields of the SST_BAL_SST_L4_REP_OBSERVATIONS_010_016 product from 1993 to 2021. See the Copernicus Marine Service Ocean State Reports (section 1.1 in Von Schuckmann et al., 2016; section 3 in Von Schuckmann et al., 2018) for more information on the OMI product. '''CONTEXT''' Sea Surface Temperature (SST) is an Essential Climate Variable (GCOS), that is an important input for initialis-ing numerical weather prediction models and fundamental for understanding air-sea interactions and moni-toring climate change (GCOS 2010). The Baltic Sea is a region that requires special attention regarding the use of satellite SST records and the assessment of climatic variability (Høyer and She 2007; Høyer and Karagali 2016). The Baltic Sea is a semi-enclosed basin affected bynatural variability, influenced by large-scale atmos-pheric processes and by the vicinity of land. In addition, the Baltic Sea is one of the largest brackish seas in the world. When analysing regional-scale climate variability, all these effects have to be considered, which re-quires dedicated regional and validated SST products. Satellite observations have previously been used to ana-lyse the climatic SST signals in the North Sea and Baltic Sea (BACC II Author Team 2015; Lehmann et al. 2011). Recently, Høyer and Karagali (2016) demonstrated that the Baltic Sea had warmed 1-2 oC from 1982 to 2012 considering all months of the year and 3-5 oC when only July-September months were considered. This was corroborated in the Ocean State Reports (section 1.1 in Von Schuckmann et al., 2016 and section 3 in Von Schuckmann et al., 2018). '''CMEMS KEY FINDINGS''' The basin-average trend of SST anomalies for Baltic Sea region amounts to 0.049±0.006 °C/year over the pe-riod 1993-2021 which corresponds to an average warming of 1.42°C. Adding the North Sea area, the average trend amounts to 0.03±0.003 °C/year over the same period, which corresponds to an average warming of 0.87°C for the entire region since 1993. '''DOI (product):''' https://doi.org/10.48670/moi-00205

'''DEFINITION''' The OMI_EXTREME_SL_NORTHWESTSHELF_slev_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 sea level measured by tide gauges along the coast. The use of percentiles instead of annual maximum and minimum values, makes this extremes study less affected by individual data measurement errors. The annual percentiles referred to annual mean sea level are temporally averaged and their spatial evolution is displayed in the dataset omi_extreme_sl_northwestshelf_slev_mean_and_anomaly_obs, 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''' Sea level (SLEV) is one of the Essential Ocean Variables most affected by climate change. Global mean sea level rise has accelerated since the 1990’s (Abram et al., 2019, Legeais et al., 2020), due to the increase of ocean temperature and mass volume caused by land ice melting (WCRP, 2018). Basin scale oceanographic and meteorological features lead to regional variations of this trend that combined with changes in the frequency and intensity of storms could also rise extreme sea levels up to one metre by the end of the century (Vousdoukas et al., 2020, Tebaldi et al., 2021). This will significantly increase coastal vulnerability to storms, with important consequences on the extent of flooding events, coastal erosion and damage to infrastructures caused by waves (Boumis et al., 2023). The increase in extreme sea levels over recent decades is, therefore, primarily due to the rise in mean sea level. Note, however, that the methodology used to compute this OMI removes the annual 50th percentile, thereby discarding the mean sea level trend to isolate changes in storminess. The North West Shelf area presents positive sea level trends with higher trend estimates in the German Bight and around Denmark, and lower trends around the southern part of Great Britain (Dettmering et al., 2021). '''COPERNICUS MARINE SERVICE KEY FINDINGS''' The completeness index criteria is fulfilled in this region by 34 stations, eight more than in 2021 (26), most of them from Norway. The mean 99th percentiles present a large spatial variability related to the tidal pattern, with largest values found in East England and at the entrance of the English channel, and lowest values along the Danish and Swedish coasts, ranging from the 3.08 m above mean sea level in Immingan (East England) to 0.57 m above mean sea level in Ringhals (Sweden) and Helgeroa (Norway). The standard deviation of annual 99th percentiles ranges between 2-3 cm in the western part of the region (e.g.: 2 cm in Harwich, 3 cm in Dunkerke) and 7-8 cm in the eastern part and the Kattegat (e.g. 8 cm in Stenungsund, Sweden).. The 99th percentile anomalies for 2022 show positive values in Southeast England, with a maximum value of +8 cm in Lowestoft, and negative values in the eastern part of the Kattegat, reaching -8 cm in Oslo. The remaining stations exhibit minor positive or negative values. '''DOI (product):''' https://doi.org/10.48670/moi-00272

'''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. '''Figure caption''' Time series of annual mean area averaged ocean heat content in the Iberia-Biscay-Ireland region (basin wide) and integrated over the 0-2000m depth layer during 1993-2022: ensemble mean (blue line) and ensemble spread (shaded area). The ensemble mean is based on different data products i.e., the IBI Reanalysis, global ocean reanalysis, and the global observational based products CORA, and ARMOR3D. Trend of ensemble mean (dashed line and bottom-right box) with 95% confidence interval computed in the period 1993-2022. Details on the products are given in the corresponding PUM and QUID for this OMI. '''DOI (product):''' https://doi.org/10.48670/mds-00316