'''DEFINITION''' The CMEMS NORTHWESTSHELF_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 North-West Shelf Multi Year Product (NWSHELF_MULTIYEAR_PHY_004_009) and the Analysis product (NORTHWESTSHELF_ANALYSIS_FORECAST_PHY_004_013). Two parameters are included on 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-2019). * Anomaly of the 99th percentile in 2020: The 99th percentile of the year 2020 is computed from the Analysis product. The anomaly is obtained by subtracting the mean percentile from the 2020 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''' This domain comprises the North West European continental shelf where depths do not exceed 200m and deeper Atlantic waters to the North and West. For these deeper waters, the North-South temperature gradient dominates (Liu and Tanhua, 2021). Temperature over the continental shelf is affected also by the various local currents in this region and by the shallow depth of the water (Elliott et al., 1990). Atmospheric heat waves can warm the whole water column, especially in the southern North Sea, much of which is no more than 30m deep (Holt et al., 2012). Warm summertime water observed in the Norwegian trench is outflow heading North from the Baltic Sea and from the North Sea itself. '''CMEMS KEY FINDINGS''' The 99th percentile SST product can be considered to represent approximately the warmest 4 days for the sea surface in Summer. Maximum anomalies for 2020 are up to 4oC warmer than the 1993-2019 average in the western approaches, Celtic and Irish Seas, English Channel and the southern North Sea. For the atmosphere, Summer 2020 was exceptionally warm and sunny in southern UK (Kendon et al., 2021), with heatwaves in June and August. Further north in the UK, the atmosphere was closer to long-term average temperatures. Overall, the 99th percentile SST anomalies show a similar pattern, with the exceptional warm anomalies in the south of the domain. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product)''' https://doi.org/10.48670/moi-00273

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

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

'''DEFINITION''' The ocean monitoring indicator on mean sea level has been presented in the Copernicus Ocean State Report #8. The ocean monitoring indicator on mean sea level is derived from the DUACS delayed-time (DT-2024 version, “my” (multi-year) dataset used when available, “myint” (multi-year interim) used after) 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 by 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 Global Ocean 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 global GIA correction of -0.3mm/yr (common global GIA correction, see Spada, 2017). 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 of 0.3 mm/yr is provided at 90% confidence level using altimeter error budget (Quet et al 2024 [in prep.]). 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 depending on the 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). According to the recent IPCC 6th assessment report (IPCC WGI, 2021), global mean sea level (GMSL) increased by 0.20 [0.15 to 0.25] m over the period 1901 to 2018 with a rate of rise that has accelerated since the 1960s to 3.7 [3.2 to 4.2] mm/yr 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). 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, 2022b). Adaptation and mitigation measures such as the restoration of mangroves and coastal wetlands, reduce the risks from sea level rise (IPCC, 2022c). ""KEY FINDINGS "" Over the [1999/02/21, 2023/12/31] period, global mean sea level rises at an average rate of 3.6  0.3 mm/year. This trend estimation is based on the altimeter measurements corrected from the global GIA correction (Spada, 2017) 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). The observed global trend agrees with other recent estimates (Oppenheimer et al., 2019; IPCC WGI, 2021). About 30% of this rise can be attributed to ocean thermal expansion (WCRP Global Sea Level Budget Group, 2018; von Schuckmann et al., 2018), 60% is due to land ice melt from glaciers and from the Antarctic and Greenland ice sheets. The remaining 10% is attributed to changes in land water storage, such as soil moisture, surface water and groundwater. From year to year, the global mean sea level record shows significant variations related mainly to the El Niño Southern Oscillation (Cazenave and Cozannet, 2014). '''DOI (product):''' https://doi.org/10.48670/moi-00237

'''DEFINITION''' The OMI_EXTREME_WAVE_IBI_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, 2011; Young & Ribal, 2019). '''COPERNICUS MARINE SERVICE KEY FINDINGS''' The mean 99th percentiles showed in the area present a wide range from 2-3.5m in the Canary Island with 0.1-0.3 m of standard deviation (std), 3.5m in the Gulf of Cadiz with 0.5m of std, 3-6m in the English Channel and the Irish Sea with 0.5-0.6m of std, 4-7m in the Bay of Biscay with 0.4-0.9m of std to 8-10m in the West of the British Isles with 0.7-1.4m of std. Results for this year show close to zero anomalies in the Canary Island (-0.2/+0.1m), the Gulf of Cadiz (-0.2m) and the English Channel and the Irish Sea (-0.1/+0.1), a general slight negative anomaly in the Bay of Biscay reaching -0.7m but inside the range of the standard deviation, and a positive anomaly (+1.0/+1.55m) in the West of the British Isles, barely out of the standard deviation range in the area. '''DOI (product):''' https://doi.org/10.48670/moi-00250

'''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''' 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 CMEMS MEDSEA_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 (MEDSEA_MULTIYEAR_WAV_006_012) and the Analysis product (MEDSEA_ANALYSIS_FORECAST_WAV_006_017). Two parameters have been considered for this OMI: * Map of the 99th mean percentile: It is obtained from the Multy Year Product, the annual 99th percentile is computed for each year of the product. The percentiles are temporally averaged in the whole period (1993-2019). * Anomaly of the 99th percentile in 2020: The 99th percentile of the year 2020 is computed from the Analysis product. The anomaly is obtained by subtracting the mean percentile to the percentile in 2020. 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 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., 2014). Complementarily, there is also an emerging question about the role of anthropogenic global climate change on present and future extreme wave conditions. The Mediterranean Sea is an almost enclosed basin where the complexity of its orographic characteristics deeply influences the atmospheric circulation at local scale, giving rise to strong regional wind regimes (Drobinski et al. 2018). Therefore, since waves are primarily driven by winds, high waves are present over most of the Mediterranean Sea and tend to reach the highest values where strong wind and long fetch (i.e. the horizontal distance over which wave-generating winds blow) are simultaneously present (Lionello et al. 2006). Specifically, as seen in figure and in agreement with other studies (e.g. Sartini et al. 2017), the highest values (5 – 6 m in figure, top) extend from the Gulf of Lion to the southwestern Sardinia through the Balearic Sea and are sustained southwards approaching the Algerian coast. They result from northerly winds dominant in the western Mediterranean Sea (Mistral or Tramontana), that become stronger due to orographic effects (Menendez et al. 2014), and act over a large area. In the Ionian Sea, the northerly Mistral wind is still the main cause of high waves (4-5 m in figure, top). In the Aegean and Levantine Seas, high waves (4-5 m in figure, top) are caused by the northerly Bora winds, prevalent in winter, and the northerly Etesian winds, prevalent in summer (Lionello et al. 2006; Chronis et al. 2011; Menendez et al. 2014). In general, northerly winds are responsible for most high waves in the Mediterranean (e.g. Chronis et al. 2011; Menendez et al. 2014). In agreement with figure (top), studies on the eastern Mediterranean and the Hellenic Seas have found that the typical wave height range in the Aegean Sea is similar to the one observed in the Ionian Sea despite the shorter fetches characterizing the former basin (Zacharioudaki et al. 2015). This is because of the numerous islands in the Aegean Sea which cause wind funneling and enhance the occurrence of extreme winds and thus of extreme waves (Kotroni et al. 2001). Special mention should be made of the high waves, sustained throughout the year, observed east and west of the island of Crete, i.e. around the exiting points of the northerly airflow in the Aegean Sea (Zacharioudaki et al. 2015). This airflow is characterized by consistently high magnitudes that are sustained during all seasons in contrast to other airflows in the Mediterranean Sea that exhibit a more pronounced seasonality (Chronis et al. 2011). '''CMEMS KEY FINDINGS''' In 2020 (bottom panel), higher-than-average values of the 99th percentile of Significant Wave Height are seen over most of the northern Mediterranean Sea, in the eastern Alboran Sea, and along stretches of the African coast (Tunisia, Libya and Egypt). In many cases they exceed the climatic standard deviation. Regions where the climatic standard deviation is exceeded twice are the European and African coast of the eastern Alboran Sea, a considerable part of the eastern Spanish coast, the Ligurian Sea and part of the east coast of France as well as areas of the southern Adriatic. These anomalies correspond to the maximum positive anomalies computed in the Mediterranean Sea for year 2020 with values that reach up to 1.1 m. Spatially constrained maxima are also found at other coastal stretches (e.g. Algeri, southeast Sardinia). Part of the positive anomalies found along the French and Spanish coast, including the coast of the Balearic Islands, can be associated with the wind storm “Gloria” (19/1 – 24/1) during which exceptional eastern winds originated in the Ligurian Sea and propagated westwards. The storm, which was of a particularly high intensity and long duration, caused record breaking wave heights in the region, and, in return, great damage to the coast (Amores et al., 2020; de Alfonso et al., 2021). Other storms that could have contributed to the positive anomalies observed in the western Mediterranean Sea include: storm Karine (25/2 – 5/4), which caused high waves from the eastern coast of Spain to the Balearic Islands (Copernicus, Climate Change Service, 2020); storm Bernardo (7/11 – 18/11) which also affected the Balearic islands and the Algerian coast and; storm Hervé (2/2 – 8/2) during which the highest wind gust was recorded at north Corsica (Wikiwand, 2021). In the eastern Mediterranean Sea, the medicane Ianos (14/9 – 21/9) may have contributed to the positive anomalies shown in the central Ionian Sea since this area coincides with the area of peak wave height values during the medicane (Copernicus, 2020a and Copernicus, 2020b). Otherwise, higher-than-average values in the figure are the result of severe, yet not unusual, wind events, which occurred during the year. Negative anomalies occur over most of the southern Mediterranean Sea, east of the Alboran Sea. The maximum negative anomalies reach about -1 m and are located in the southeastern Ionian Sea and west of the south part of mainland Greece as well as in coastal locations of the north and east Aegean They appear to be quite unusual since they are greater than two times the climatic standard deviation in the region. They could imply less severe southerly wind activity during 2020 (Drobinski et al., 2018). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00262

'''DEFINITION''' Estimates of Ocean Heat Content (OHC) are obtained from integrated differences of the measured temperature and a climatology along a vertical profile in the ocean (von Schuckmann et al., 2018). The regional OHC values are then averaged from 60°S-60°N aiming i) to obtain the mean OHC as expressed in Joules per meter square (J/m2) to monitor the large-scale variability and change. ii) to monitor the amount of energy in the form of heat stored in the ocean (i.e. the change of OHC in time), expressed in Watt per square meter (W/m2). Ocean heat content is one of the six Global Climate Indicators recommended by the World Meterological Organisation for Sustainable Development Goal 13 implementation (WMO, 2017). '''CONTEXT''' Knowing how much and where heat energy is stored and released in the ocean is essential for understanding the contemporary Earth system state, variability and change, as the ocean shapes our perspectives for the future (von Schuckmann et al., 2020). Variations in OHC can induce changes in ocean stratification, currents, sea ice and ice shelfs (IPCC, 2019; 2021); they set time scales and dominate Earth system adjustments to climate variability and change (Hansen et al., 2011); they are a key player in ocean-atmosphere interactions and sea level change (WCRP, 2018) and they can impact marine ecosystems and human livelihoods (IPCC, 2019). '''CMEMS KEY FINDINGS''' Regional trends for the period 2005-2019 from the Copernicus Marine Service multi-ensemble approach show warming at rates ranging from the global mean average up to more than 8 W/m2 in some specific regions (e.g. northern hemisphere western boundary current regimes). There are specific regions where a negative trend is observed above noise at rates up to about -5 W/m2 such as in the subpolar North Atlantic, or the western tropical Pacific. These areas are characterized by strong year-to-year variability (Dubois et al., 2018; Capotondi et al., 2020). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00236

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