'''DEFINITION''' The subsurface temperature trends have been derived from regional reanalysis results for the Baltic Sea (product references BALTICSEA_REANALYSIS_PHY_003_011). Horizontal averaging has been conducted over the Baltic Sea domain (13 °E - 31 °E and 53 °N - 66 °N; excluding the Skagerrak strait). The temperature trend has been obtained through a linear fit for each time series of horizontally averaged annual temperature and at each depth level (Mulet et al., 2018). '''CONTEXT''' The Baltic Sea is a semi-enclosed sea in North-Eastern Europe. The temperature of the upper mixed layer of the Baltic Sea is characterized by a strong seasonal cycle driven by the annual course of solar radiation (Leppäranta and Myrberg, 2008). The maximum water temperatures in the upper layer are reached in July and August and the minimum during February, when the Baltic Sea becomes partially frozen (CMEMS OMI Baltic Sea Sea Ice Extent, CMEMS OMI Baltic Sea Sea Ice Volume). Seasonal thermocline, developing in the depth range of 10-30 m in spring, reaches its maximum strength in summer and is eroded in autumn. During autumn and winter the Baltic Sea is thermally mixed down to the permanent halocline in the depth range of 60-80 meters (Matthäus, 1984). The 20–50 m thick cold intermediate layer forms below the upper mixed layer in March and is observed until October within the 15-65 m depth range (Chubarenko and Stepanova, 2018; Liblik and Lips, 2011). The deep layers of the Baltic Sea are disconnected from the ventilated upper ocean layers, and temperature variations are predominantly driven by mixing processes and horizontal advection. A warming trend of the sea surface waters is positively correlated with the increasing trend of diffuse attenuation of light (Kd490) and satellite-detected chlorophyll concentration (Kahru et al., 2016). Temperature increase in the water column could accelerate oxygen consumption during organic matter oxidation (Savchuk, 2018). '''CMEMS KEY FINDINGS''' The subsurface temperature over the 1993-2020 period shows warming trends of about 0.05 °C/year at all depths. The largest warming trend of 0.06 °C/year is recorded at the 20 m depth, which corresponds to seasonal thermocline. Similar positive trend is at the depth of 60-70 meters, which corresponds to the depth of the upper part of the permanent halocline. A positive trend in the sea surface waters has been detected since the 1990s (BACCII Author Team, 2015) as well as a decreasing trend of the start day of the spring phytoplankton bloom (Raudsepp et al., 2019; Kahru et al., 2016). From the measurements Savchuk (2018) has calculated the temperature trend of 0.04◦oC/year since 1979 on average in the deep layers (>60m) of the Baltic Proper. '''DOI (product):''' https://doi.org/10.48670/moi-00208

'''DEFINITION''' The northwestshelf_omi_tempsal_sst_trend product includes the Sea Surface Temperature (SST) trend for the European North West Shelf Seas over the period 1993-2020, 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-ATL-SST.pdf, Mulet et al., 2018), which provided the SSTs used to compute the SST trend over the European North West Shelf Seas. 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 due to its role for climate variability and change (Deser et al., 2010; IPCC, 2021). 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''' Over the period 1993-2020, the European North West Shelf Seas area experiences overall surface warming, particularly in the eastern part of the area such as the Kattegat region. '''DOI (product):''' https://doi.org/10.48670/moi-00276

'''DEFINITION''' The sea level ocean monitoring indicator has been presented in the Copernicus Ocean State Report #8. The sea level ocean monitoring indicator is derived from the DUACS delayed-time (DT-2024 version, “my” (multi-year) dataset used when available) sea level anomaly maps from satellite altimetry 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 time series of area averaged anomalies correspond to the area average of the maps in the Baltic Sea weighted by the cosine of the latitude (to consider the changing area in each grid with latitude) and by the proportion of ocean in each grid (to consider the coastal areas). The time series are corrected from regional mean GIA correction (weighted GIA mean of a 27 ensembles model following Spada et Melini, 2019). The time series are adjusted for seasonal annual and semi-annual signals and low-pass filtered at 6 months. Then, the trends/accelerations are estimated on the time series using ordinary least square fit. The trend uncertainty is provided in a 90% confidence interval. It is calculated as the weighted mean uncertainties in the region from 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 considered. '''CONTEXT''' Change in mean sea level is an essential indicator of our evolving climate, as it reflects both the thermal expansion of the ocean in response to its warming and the increase in ocean mass due to the melting of ice sheets and glaciers (WCRP Global Sea Level Budget Group, 2018). At regional scale, sea level does not change homogenously. It is 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). The adverse effects of floods, storms and tropical cyclones, and the resulting losses and damage, have increased as a result of rising sea levels, increasing people and infrastructure vulnerability and food security risks, particularly in low-lying areas and island states (IPCC, 2022a). Adaptation and mitigation measures such as the restoration of mangroves and coastal wetlands, reduce the risks from sea level rise (IPCC, 2022b). The Baltic Sea is a relatively small semi-enclosed basin with shallow bathymetry. Different forcings have been discussed to trigger sea level variations in the Baltic Sea at different time scales. In addition to steric effects, decadal and longer sea level variability in the basin can be induced by sea water exchange with the North Sea, and in response to atmospheric forcing and climate variability (e.g., the North Atlantic Oscillation; Gräwe et al., 2019). '''KEY FINDINGS''' Over the [1999/02/20 to 2023/12/31] period, the area-averaged sea level in the Baltic Sea rises at a rate of 4.8 ± 0.8 mm/year with an acceleration of 0.12 ± 0.07 mm/year2. This trend estimation is based on the altimeter measurements corrected from regional GIA correction (Spada et Melini, 2019) to consider the ongoing movement of land. The TOPEX-A is no longer included in the computation of regional mean sea level parameters (trend and acceleration) with version 2024 products due to potential drifts, and ongoing work aims to develop a new empirical correction. Calculation begins in February 1999 (the start of the TOPEX-B period). '''DOI (product):''' https://doi.org/10.48670/moi-00202

'''DEFINITION''' The sea level ocean monitoring indicator 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 are also available in the CMEMS catalogue (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). The mean sea level evolution estimated in the Black Sea 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 adjusted and the time series is low-pass filtered. The curve is corrected for the effect of the Glacial Isostatic Adjustment using the ICE5G-VM2 GIA model (Peltier, 2004). '''CONTEXT''' Mean sea level evolution has a direct impact on coastal areas and is a crucial index of climate change since it reflects both the amount of heat added in the ocean and the mass loss due to land ice melt (e.g. IPCC, 2013; Dieng et al., 2017). Long-term and inter-annual variations of the sea level are observed at global and regional scales. They are strongly related to the internal variability observed at basin scale and these variations can strongly affect population living in coastal areas. '''CMEMS KEY FINDINGS''' Using the latest reprocessed altimeter sea level products, it is possible to estimate the sea level rise in the Black Sea since 01/1993 (see the proposed figure of the indicator for the updated trend value). The associated uncertainty is provided in a 90% confidence interval and only errors related to the altimeter observation system have been considered in the sea level trend uncertainty (Prandi et al., 2021). The uncertainty due to the sea level internal variability of the observed ocean is not included and should be considered additionally. The Black Sea is a relatively small semi-enclosed basin with shallow bathymetry, which explains the high level of inter annual variability observed in the sea level record compared to large, deeper and open ocean areas (North West Shelf and Iberian-Biscay-Ireland regions). '''DOI (product):''' https://doi.org/10.48670/moi-00215

'''DEFINITION''' The sea level ocean monitoring indicator 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 are also available in the CMEMS catalogue (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). The regional sea level trends are derived from a linear fit of the altimeter sea level maps. The altimeter data have not been corrected for the effect of the Glacial Isostatic Adjustment nor the TOPEX-A instrumental drift during the period 1993-1998. '''CONTEXT''' Mean sea level evolution has a direct impact on coastal areas and is a crucial index of climate change since it reflects both the amount of heat added in the ocean and the mass loss due to land ice melt (e.g. IPCC, 2013; Dieng et al., 2017). Long-term and inter-annual variations of the sea level are observed at global and regional scales. They are strongly related to the internal variability observed at basin scale and these variations can strongly affect population living in coastal areas. '''CMEMS KEY FINDINGS''' The altimeter mean sea level trends since 1993 exhibit large-scale variations with amplitudes reaching up to +8 mm/yr in regions such as the western tropical Pacific Ocean. In this area, trends are mainly of thermosteric origin (Legeais et al., 2016; Meyssignac et al., 2017) in response to increased easterly winds during the last two decades associated with the decreasing Interdecadal Pacific Oscillation (IPO)/Pacific Decadal Oscillation (e.g. McGregor et al. 2012; Merrifield et al. 2012; Palanisamy et al. 2014; Han et al. 2010; Rietbroek et al. 2016). Prandi et al. (2021) have estimated a regional altimeter sea level error budget from which they determine a regional error variance-covariance matrix and they provide uncertainties of the regional sea level trends. Only the contribution of the measurement system is considered (the contribution from the ocean natural variability is not considered). Over 1993-2019, the averaged local sea level trend uncertainty is around 0.83 mm/yr with local values ranging from 0.78 to 1.22 mm/yr. '''DOI (product):''' https://doi.org/10.48670/moi-00238

'''DEFINITION''' The ibi_omi_tempsal_sst_trend product includes the Sea Surface Temperature (SST) trend for the Iberia-Biscay-Irish Seas over the period 1993-2020, 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-ATL-SST.pdf), which provided the SSTs used to compute the SST trend over the Iberia-Biscay-Irish Seas. 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 identified as Essential Climate Variable (GCOS 2010), as it is used to analyze and monitor climate variability and change (e.g. Deser et al., 2010). In addition, SST anomalies are used for the analysis of extreme events (marine heatwaves, Hobday et al., 2018). '''CMEMS KEY FINDINGS''' Over the period 1993-2020, most of the Iberia-Biscay-Irish Seas area shows overall surface warming, particularly in the northeastern part of the region. '''DOI (product):''' https://doi.org/10.48670/moi-00257

'''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. '''DOI (product):''' https://doi.org/10.48670/moi-00245

'''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. 2021) 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.0021, 0.0024] pH unit (see Quality Information Document for more details). The variation on the trend estimation amounts to 0.0006 pH unit 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 scientific publication is in preparation for this indicator. '''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 contributes to SDG14 (https://sustainabledevelopment.un.org/sdg14). '''CMEMS KEY FINDINGS ''' Since the year 1985, global ocean surface pH is decreasing at a rate of -0.00160.0006 per year. '''DOI (product):''' https://doi.org/10.48670/moi-00224

'''DEFINITION''' The subsurface salinity trends have been derived from regional reanalysis and forecast modelling results of the CMEMS BAL MFC group for the Baltic Sea (product references BALTICSEA_REANALYSIS_PHY_003_011). The salinity trend has been obtained through a linear fit for each time series of horizontally averaged (13 °E - 31 °E and 53 °N - 66 °N; excluding the Skagerrak strait) annual salinity and at each depth level (Mulet et al., 2018). '''CONTEXT''' The Baltic Sea is a brackish semi-enclosed sea in North-Eastern Europe. The surface salinity varies horizontally from ~10 near the Danish Straits down to ~2 at the northernmost and easternmost sub-basins of the Baltic Sea. The halocline, a vertical layer with rapid changes of salinity with depth that separates the well-mixed surface layer from the weakly stratified layer below, is located at the depth range of 60-80 meters (Matthäus, 1984). The bottom layer salinity below the halocline depth varies from 15 in the south down to 3 in the northern Baltic Sea (Väli et al., 2013). The long-term salinity is determined by net precipitation and river discharge as well as saline water inflows from the North Sea (BACCII Author Team, 2015). Long-term salinity decrease may reduce the occurrence and biomass of the Fucus vesiculosus - Idotea balthica association/symbiotic aggregations (Kotta et al., 2019). Changes in salinity and oxygen content affect the survival of the Baltic cod eggs (Raudsepp et al, 2019; von Dewitz et al., 2018). '''CMEMS KEY FINDINGS''' The subsurface salinity over the 1993-2020 period shows no trend in the surface layer of 40-m. Then, the trend increases monotonically until the depth of 80-m, where it reaches the value of about 0.04 per year. Below that depth the salinity trend is uniform within the uncertainty limits. '''DOI (product):''' https://doi.org/10.48670/moi-00207

'''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 (Faizal and Rafiuddin, 2011). 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_REANALYSIS_PHY_001_031): GLORYS, C-GLORS, ORAS5, FOAM Two observation based products: CORA (INSITU_GLO_TS_REP_OBSERVATIONS_013_001_b) and ARMOR3D (GLOBAL_REP_PHY_001_021). 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-2018 at rate of 1.5±0.2 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