'''Short description:''' For The Global Ocean - The GHRSST Multi-Product Ensemble (GMPE) system has been implemented at the Met Office which takes inputs from various analysis production centres on a routine basis and produces ensemble products at 0.25deg.x0.25deg. horizontal resolution. A large number of sea surface temperature (SST) analyses are produced by various institutes around the world, making use of the SST observations provided by the Global High Resolution SST (GHRSST) project. These are used by a number of groups including: numerical weather prediction centres; ocean forecasting groups; climate monitoring and research groups. There is a requirement to develop international collaboration in this field in order to assess and inter-compare the different analyses, and to provide uncertainty estimates on both the analyses and observational products. The GMPE system has been developed for these purposes and is run on a daily basis at the Met Office, producing global ensemble median and standard deviations for SST on a regular 0.25 degree resolution global grid. '''DOI (product) :''' https://doi.org/10.48670/mds-00378

'''Short description:''' For the Baltic Sea- The DMI Sea Surface Temperature L3S aims at providing daily multi-sensor supercollated data at 0.03deg. x 0.03deg. horizontal resolution, using satellite data from infra-red radiometers. Uses SST satellite products from these sensors: NOAA AVHRRs 7, 9, 11, 14, 16, 17, 18 , Envisat ATSR1, ATSR2 and AATSR. '''DOI (product) :''' https://doi.org/10.48670/moi-00154

'''DEFINITION''' The global yearly ocean CO2 sink represents the ocean uptake of CO2 from the atmosphere computed over the whole ocean. It is expressed in PgC per year. The ocean monitoring index is presented for the period 1985 to year-1. The yearly estimate of the ocean CO2 sink corresponds to the mean of a 100-member ensemble of CO2 flux estimates (Chau et al. 2022). The range of an estimate with the associated uncertainty is then defined by the empirical 68% interval computed from the ensemble. '''CONTEXT''' Since the onset of the industrial era in 1750, the atmospheric CO2 concentration has increased from about 277±3 ppm (Joos and Spahni, 2008) to 412.44±0.1 ppm in 2020 (Dlugokencky and Tans, 2020). By 2011, the ocean had absorbed approximately 28 ± 5% of all anthropogenic CO2 emissions, thus providing negative feedback to global warming and climate change (Ciais et al., 2013). The ocean CO2 sink is evaluated every year as part of the Global Carbon Budget (Friedlingstein et al. 2022). The uptake of CO2 occurs primarily in response to increasing atmospheric levels. The global flux is characterized by a significant variability on interannual to decadal time scales largely in response to natural climate variability (e.g., ENSO) (Friedlingstein et al. 2022, Chau et al. 2022). '''CMEMS KEY FINDINGS''' The rate of change of the integrated yearly surface downward flux has increased by 0.04±0.01e-1 PgC/yr2 over the period 1985 to year-1. The yearly flux time series shows a plateau in the 90s followed by an increase since 2000 with a growth rate of 0.06±0.04e-1 PgC/yr2. In 2021 (resp. 2020), the global ocean CO2 sink was 2.41±0.13 (resp. 2.50±0.12) PgC/yr. The average over the full period is 1.61±0.10 PgC/yr with an interannual variability (temporal standard deviation) of 0.46 PgC/yr. In order to compare these fluxes to Friedlingstein et al. (2022), the estimate of preindustrial outgassing of riverine carbon of 0.61 PgC/yr, which is in between the estimate by Jacobson et al. (2007) (0.45±0.18 PgC/yr) and the one by Resplandy et al. (2018) (0.78±0.41 PgC/yr) needs to be added. A full discussion regarding this OMI can be found in section 2.10 of the Ocean State Report 4 (Gehlen et al., 2020) and in Chau et al. (2022). '''DOI (product):''' https://doi.org/10.48670/moi-00223

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' For the Global ocean, the ESA Ocean Colour CCI surface Chlorophyll (mg m-3, 4 km resolution) using the OC-CCI recommended chlorophyll algorithm is made available in CMEMS format. L3 products are daily files, while the L4 are monthly composites. ESA-CCI data are provided by Plymouth Marine Laboratory at 4km resolution. These are processed using the same in-house software as in the operational processing. Standard masking criteria for detecting clouds or other contamination factors have been applied during the generation of the Rrs, i.e., land, cloud, sun glint, atmospheric correction failure, high total radiance, large solar zenith angle (actually a high air mass cutoff, but approximating to 70deg zenith), coccolithophores, negative water leaving radiance, and normalized water leaving radiance at 555 nm 0.15 Wm-2 sr-1 (McClain et al., 1995). Ocean colour technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so called ocean colour which is affected by the presence of phytoplankton. By comparing reflectances at different wavelengths and calibrating the result against in-situ measurements, an estimate of chlorophyll content can be derived. A detailed description of calibration & validation is given in the relevant QUID, associated validation reports and quality documentation. '''Processing information:''' ESA-CCI data are provided by Plymouth Marine Laboratory at 4km resolution. These are processed using the same in-house software as in the operational processing. The entire CCI data set is consistent and processing is done in one go. Both OC CCI and the REP product are versioned. Standard masking criteria for detecting clouds or other contamination factors have been applied during the generation of the Rrs, i.e., land, cloud, sun glint, atmospheric correction failure, high total radiance, large solar zenith angle (actually a high air mass cutoff, but approximating to 70deg zenith), coccolithophores, negative water leaving radiance, and normalized water leaving radiance at 555 nm 0.15 Wm-2 sr-1 (McClain et al., 1995). '''Description of observation methods/instruments:''' Ocean colour technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so called ocean colour which is affected by the presence of phytoplankton. By comparing reflectances at different wavelengths and calibrating the result against in-situ measurements, an estimate of chlorophyll content can be derived. '''Quality / Accuracy / Calibration information:''' Detailed description of cal/val is given in the relevant QUID, associated validation reports and quality documentation. '''Suitability, Expected type of users / uses:''' This product is meant for use for educational purposes and for the managing of the marine safety, marine resources, marine and coastal environment and for climate and seasonal studies. '''DOI (product) :''' https://doi.org/10.48670/moi-00101

'''Short description:''' For the '''Global''' Ocean '''Satellite Observations''', Brockmann Consult (BC) is providing '''Bio-Geo_Chemical (BGC)''' products based on the ESA-CCI inputs. * Upstreams: SeaWiFS, MODIS, MERIS, VIIRS-SNPP, OLCI-S3A & OLCI-S3B for the '''""multi""''' products. * Variables: Chlorophyll-a ('''CHL'''). * Temporal resolutions: '''monthly'''. * Spatial resolutions: '''4 km''' (multi). * Recent products are organized in datasets called Near Real Time ('''NRT''') and long time-series (from 1997) in datasets called Multi-Years ('''MY'''). To find these products in the catalogue, use the search keyword '''""ESA-CCI""'''. '''DOI (product) :''' https://doi.org/10.48670/moi-00283

''' Short description: ''' For the Mediterranean Sea - the CNR diurnal sub-skin Sea Surface Temperature (SST) product provides daily gap-free (L4) maps of hourly mean sub-skin SST at 1/16° (0.0625°) horizontal resolution over the CMEMS Mediterranean Sea (MED) domain, by combining infrared satellite and model data (Marullo et al., 2014). The implementation of this product takes advantage of the consolidated operational SST processing chains that provide daily mean SST fields over the same basin (Buongiorno Nardelli et al., 2013). The sub-skin temperature is the temperature at the base of the thermal skin layer and it is equivalent to the foundation SST at night, but during daytime it can be significantly different under favorable (clear sky and low wind) diurnal warming conditions. The sub-skin SST L4 product is created by combining geostationary satellite observations aquired from SEVIRI and model data (used as first-guess) aquired from the CMEMS MED Monitoring Forecasting Center (MFC). This approach takes advantage of geostationary satellite observations as the input signal source to produce hourly gap-free SST fields using model analyses as first-guess. The resulting SST anomaly field (satellite-model) is free, or nearly free, of any diurnal cycle, thus allowing to interpolate SST anomalies using satellite data acquired at different times of the day (Marullo et al., 2014). [https://help.marine.copernicus.eu/en/articles/4444611-how-to-cite-or-reference-copernicus-marine-products-and-services How to cite] '''DOI (product) :''' https://doi.org/10.48670/moi-00170

'''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 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 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''' Since the year 2005, the upper (0-700m) near-global (60°S-60°N) ocean warms at a rate of 0.6 ± 0.1 W/m2. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00234

'''DEFINITION''' The sea level ocean monitoring indicator has been presented in the Copernicus Ocean State Report #8. The ocean monitoring indicator on regional mean sea level 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 Irish-Biscay-Iberian (IBI) 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 ensemble 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). In IBI region, the RMSL trend is modulated by decadal variations. As observed over the global ocean, the main actors of the long-term RMSL trend are associated with anthropogenic global/regional warming. Decadal variability is mainly linked to the strengthening or weakening of the Atlantic Meridional Overturning Circulation (AMOC) (e.g. Chafik et al., 2019). The latest is driven by the North Atlantic Oscillation (NAO) for decadal (20-30y) timescales (e.g. Delworth and Zeng, 2016). Along the European coast, the NAO also influences the along-slope winds dynamic which in return significantly contributes to the local sea level variability observed (Chafik et al., 2019). ""KEY FINDINGS "" Over the [1999/02/20 to 2025/10/18] period, the area-averaged sea level in the IBI area rises at a rate of 5.0 ± 0.8 mm/yr with an acceleration of 0.29 ± 0.06 mm/yr². This trend estimation is based on the altimeter measurements corrected from global 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-00252

'''DEFINITION''' The Iberia Biscay Ireland (IBI) Sea Surface Temperature extreme from Reanalysis ocean monitoring indicator (OMI) (OMI_CLIMATE_TEMPSAL_IBI_extreme_var_temp_mean_and_anomaly) is based on the computation of the annual 99th percentile of Sea Surface Temperature (SST) from model data. Two different Copernicus Marine products are used to compute the indicator: The IBI Reanalysis (IBI_MULTIYEAR_PHY_005_002) and the IBI Analysis product (IBI_ANALYSISFORECAST_PHY_005_001). Two parameters have been considered for this OMI: * '''Map of the 99th mean percentile''': It is obtained from the reanalysis product, the annual 99th percentile is computed for each year of the product. The percentiles are temporally averaged over the whole period (1993-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 from the 2024 percentile. This indicator is aimed at monitoring the extremes of sea surface temperature every year and at checking their variations in space. The use of percentiles instead of annual maxima, makes this extremes study less affected by individual data. This study of extreme variability was first applied to the sea level variable (Pérez Gómez et al 2016) and then extended to other essential variables, such as sea surface temperature and significant wave height (Pérez Gómez et al 2018 and Alvarez Fanjul et al., 2019). More details and a full scientific evaluation can be found in the CMEMS Ocean State report (Alvarez Fanjul et al., 2019). '''CONTEXT''' The Sea Surface Temperature (SST) is one of the essential ocean variables, hence the monitoring of this variable is of key importance, since its variations can affect the ocean circulation, marine ecosystems, and ocean-atmosphere exchange processes. As the oceans continuously interact with the atmosphere, trends of sea surface temperature can also have an effect on the global climate. While the global-averaged sea surface temperatures have increased since the beginning of the 20th century (Hartmann et al., 2013) in the North Atlantic, anomalous cold conditions have also been reported since 2014 (Mulet et al., 2018; Dubois et al., 2018). The IBI area is a complex dynamic region with a remarkable variety of ocean physical processes and scales involved. The SST field in the region is strongly dependent on latitude, with higher values towards the South (Locarnini et al. 2013). This latitudinal gradient is supported by the presence of the eastern part of the North Atlantic subtropical gyre that transports cool water from the northern latitudes towards the equator. Additionally, the IBI region is under the influence of the Sea Level Pressure dipole established between the Icelandic low and the Bermuda high. Therefore, the interannual and interdecadal variability of the surface temperature field may be influenced by the North Atlantic Oscillation pattern (Czaja and Frankignoul, 2002; Flatau et al., 2003). Upwelling processes, taking place in the coastal margins, are also relevant in the IBI region. The most referenced one is the eastern boundary coastal upwelling system off the African and western Iberian coast (Sotillo et al., 2016), although other smaller upwelling systems have also been described in the northern coast of the Iberian Peninsula (Alvarez et al., 2011), the south-western Irish coast (Edwars et al., 1996) and the European Continental Slope (Dickson, 1980). '''CMEMS KEY FINDINGS''' In the IBI region, the 99th mean percentile for 1993-2023 shows a north-south pattern driven by the climatological distribution of temperatures in the North Atlantic. In the coastal regions of Africa and the Iberian Peninsula, the mean values are influenced by the upwelling processes (Sotillo et al., 2016). These results are consistent with the ones presented in Álvarez Fanjul (2019) for the period 1993-2016. The analysis of the 99th percentile SST anomaly for the year 2024 reveals that the northeastern Atlantic region, between latitudes 36° N and 48° N, experienced thermal anomalies exceeding twice the standard deviation. Similar anomalies are also observed near the northeastern Iberian Peninsula, suggesting that inshore and coastal areas may have been affected as well. In contrast, the upwelling region west of the Iberian Peninsula shows negative anomalies in maximum SST, indicating an intensification of upwelling processes in this area. '''DOI (product):''' https://doi.org/10.48670/moi-00254

'''This product has been archived''' '''DEFINITION''' The linear change of zonal mean subsurface temperature over the period 1993-2019 at each grid point (in depth and latitude) is evaluated to obtain a global mean depth-latitude plot of subsurface temperature trend, expressed in °C. The linear change is computed using the slope of the linear regression at each grid point scaled by the number of time steps (27 years, 1993-2019). A multi-product approach is used, meaning that the linear change is first computed for 5 different zonal mean temperature estimates. The average linear change is then computed, as well as the standard deviation between the five linear change computations. The evaluation method relies in the study of the consistency in between the 5 different estimates, which provides a qualitative estimate of the robustness of the indicator. See Mulet et al. (2018) for more details. '''CONTEXT''' Large-scale temperature variations in the upper layers are mainly related to the heat exchange with the atmosphere and surrounding oceanic regions, while the deeper ocean temperature in the main thermocline and below varies due to many dynamical forcing mechanisms (Bindoff et al., 2019). Together with ocean acidification and deoxygenation (IPCC, 2019), ocean warming can lead to dramatic changes in ecosystem assemblages, biodiversity, population extinctions, coral bleaching and infectious disease, change in behavior (including reproduction), as well as redistribution of habitat (e.g. Gattuso et al., 2015, Molinos et al., 2016, Ramirez et al., 2017). Ocean warming also intensifies tropical cyclones (Hoegh-Guldberg et al., 2018; Trenberth et al., 2018; Sun et al., 2017). '''CMEMS KEY FINDINGS''' The results show an overall ocean warming of the upper global ocean over the period 1993-2019, particularly in the upper 300m depth. In some areas, this warming signal reaches down to about 800m depth such as for example in the Southern Ocean south of 40°S. In other areas, the signal-to-noise ratio in the deeper ocean layers is less than two, i.e. the different products used for the ensemble mean show weak agreement. However, interannual-to-decadal fluctuations are superposed on the warming signal, and can interfere with the warming trend. For example, in the subpolar North Atlantic decadal variations such as the so called ‘cold event’ prevail (Dubois et al., 2018; Gourrion et al., 2018), and the cumulative trend over a quarter of a decade does not exceed twice the noise level below about 100m depth. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00244