he Global ARMOR3D L4 Reprocessed dataset is obtained by combining satellite (Sea Level Anomalies, Geostrophic Surface Currents, Sea Surface Temperature) and in-situ (Temperature and Salinity profiles) observations through statistical methods. References : - ARMOR3D: Guinehut S., A.-L. Dhomps, G. Larnicol and P.-Y. Le Traon, 2012: High resolution 3D temperature and salinity fields derived from in situ and satellite observations. Ocean Sci., 8(5):845–857. - ARMOR3D: Guinehut S., P.-Y. Le Traon, G. Larnicol and S. Philipps, 2004: Combining Argo and remote-sensing data to estimate the ocean three-dimensional temperature fields - A first approach based on simulated observations. J. Mar. Sys., 46 (1-4), 85-98. - ARMOR3D: Mulet, S., M.-H. Rio, A. Mignot, S. Guinehut and R. Morrow, 2012: A new estimate of the global 3D geostrophic ocean circulation based on satellite data and in-situ measurements. Deep Sea Research Part II : Topical Studies in Oceanography, 77–80(0):70–81.

'''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. 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_MULTIYEAR_PHY_ENS_001_031): GLORYS, C-GLORS, ORAS5, FOAM • Two observation based products: CORA (INSITU_GLO_PHY_TS_OA_MY_013_052) and ARMOR3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012). 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-2022 at rate of 1.38±0.08 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

'''This product has been archived''' '''DEFINITION''' The temporal evolution of thermosteric sea level in an ocean layer is obtained from an integration of temperature driven ocean density variations, which are subtracted from a reference climatology to obtain the fluctuations from an average field. The regional thermosteric sea level values are then averaged from 60°S-60°N aiming to monitor interannual to long term global sea level variations caused by temperature driven ocean volume changes through thermal expansion as expressed in meters (m). '''CONTEXT''' Most of the interannual variability and trends in regional sea level is caused by changes in steric sea level. At mid and low latitudes, the steric sea level signal is essentially due to temperature changes, i.e. the thermosteric effect (Stammer et al., 2013, Meyssignac et al., 2016). Salinity changes play only a local role. Regional trends of thermosteric sea level can be significantly larger compared to their globally averaged versions (Storto et al., 2018). Except for shallow shelf sea and high latitudes (> 60° latitude), regional thermosteric sea level variations are mostly related to ocean circulation changes, in particular in the tropics where the sea level variations and trends are the most intense over the last two decades. '''CMEMS KEY FINDINGS''' Significant (i.e. when the signal exceeds the noise) regional trends for the period 2005-2019 from the Copernicus Marine Service multi-ensemble approach show a thermosteric sea level rise at rates ranging from the global mean average up to more than 8 mm/year. There are specific regions where a negative trend is observed above noise at rates up to about -8 mm/year 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-00241

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' Oligotrophic subtropical gyres are regions of the ocean with low levels of nutrients required for phytoplankton growth and low levels of surface chlorophyll-a whose concentration can be quantified through satellite observations. The gyre boundary has been defined using a threshold value of 0.15 mg m-3 chlorophyll for the Atlantic gyres (Aiken et al. 2016), and 0.07 mg m-3 for the Pacific gyres (Polovina et al. 2008). The area inside the gyres for each month is computed using monthly chlorophyll data from which the monthly climatology is subtracted to compute anomalies. A gap filling algorithm has been utilized to account for missing data. Trends in the area anomaly are then calculated for the entire study period (September 1997 to December 2020). '''CONTEXT''' Oligotrophic gyres of the oceans have been referred to as ocean deserts (Polovina et al. 2008). They are vast, covering approximately 50% of the Earth’s surface (Aiken et al. 2016). Despite low productivity, these regions contribute significantly to global productivity due to their immense size (McClain et al. 2004). Even modest changes in their size can have large impacts on a variety of global biogeochemical cycles and on trends in chlorophyll (Signorini et al. 2015). Based on satellite data, Polovina et al. (2008) showed that the areas of subtropical gyres were expanding. The Ocean State Report (Sathyendranath et al. 2018) showed that the trends had reversed in the Pacific for the time segment from January 2007 to December 2016. '''CMEMS KEY FINDINGS''' The trend in the North Atlantic gyre area for the 1997 Sept – 2020 December period was positive, with a 0.39% year-1 increase in area relative to 2000-01-01 values. This trend has decreased compared with the 1997-2019 trend of 0.45%, and is statistically significant (p<0.05). During the 1997 Sept – 2020 December period, the trend in chlorophyll concentration was positive (0.24% year-1) inside the North Atlantic gyre relative to 2000-01-01 values. This time series extension has resulted in a reversal in the rate of change, compared with the -0.18% trend for the 1997-209 period and is statistically significant (p<0.05). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00226

'''DEFINITION''' The CMEMS IBI_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 Iberia-Biscay-Ireland Multi Year Product (IBI_MULTIYEAR_PHY_005_002) and the 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 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-2021). • Anomaly of the 99th percentile in 2022: The 99th percentile of the year 2022 is computed from the Analysis product. The anomaly is obtained by subtracting the mean percentile from the 2022 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 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 Sea Surface Temperature 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 Iberia-Biscay-Ireland 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). Also relevant in the region are the upwelling processes taking place in the coastal margins. 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-2021 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 anomaly in the year 2023 shows that this period has been affected by a severe impact of maximum SST values. Anomalies exceeding the standard deviation affect almost the entire IBI domain, and regions impacted by thermal anomalies surpassing twice the standard deviation are also widespread below the 43ºN parallel. Extreme SST values exceeding twice the standard deviation affect not only the open ocean waters but also the easter boundary upwelling areas such as the northern half of Portugal, the Spanish Atlantic coast up to Cape Ortegal, and the African coast south of Cape Aguer. It is worth noting the impact of anomalies that exceed twice the standard deviation is widespread throughout the entire Mediterranean region included in this analysis. '''Figure caption''' Iberia-Biscay-Ireland Surface Temperature extreme variability: Map of the 99th mean percentile computed from the Multi Year Product (left panel) and anomaly of the 99th percentile in 2022 computed from the Analysis product (right panel). Transparent grey areas (if any) represent regions where anomaly exceeds the climatic standard deviation (light grey) and twice the climatic standard deviation (dark grey). '''DOI (product):''' https://doi.org/10.48670/moi-00254

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' Global Ocean- in-situ reprocessed Carbon observations. This product contains observations and gridded files from two up-to-date carbon and biogeochemistry community data products: Surface Ocean Carbon ATlas SOCATv2021 and GLobal Ocean Data Analysis Project GLODAPv2.2021. The SOCATv2022-OBS dataset contains >25 million observations of fugacity of CO2 of the surface global ocean from 1957 to early 2022. The quality control procedures are described in Bakker et al. (2016). These observations form the basis of the gridded products included in SOCATv2020-GRIDDED: monthly, yearly and decadal averages of fCO2 over a 1x1 degree grid over the global ocean, and a 0.25x0.25 degree, monthly average for the coastal ocean. GLODAPv2.2022-OBS contains >1 million observations from individual seawater samples of temperature, salinity, oxygen, nutrients, dissolved inorganic carbon, total alkalinity and pH from 1972 to 2020. These data were subjected to an extensive quality control and bias correction described in Olsen et al. (2020). GLODAPv2-GRIDDED contains global climatologies for temperature, salinity, oxygen, nitrate, phosphate, silicate, dissolved inorganic carbon, total alkalinity and pH over a 1x1 degree horizontal grid and 33 standard depths using the observations from the previous iteration of GLODAP, GLODAPv2. SOCAT and GLODAP are based on community, largely volunteer efforts, and the data providers will appreciate that those who use the data cite the corresponding articles (see References below) in order to support future sustainability of the data products. '''DOI (product) :''' https://doi.org/10.48670/moi-00035

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' We have derived an annual eutrophication and eutrophication indicator map for the North Atlantic Ocean using satellite-derived chlorophyll concentration. Using the satellite-derived chlorophyll products distributed in the regional North Atlantic CMEMS REP Ocean Colour dataset (OC- CCI), we derived P90 and P10 daily climatologies. The time period selected for the climatology was 1998-2017. For a given pixel, P90 and P10 were defined as dynamic thresholds such as 90% of the 1998-2017 chlorophyll values for that pixel were below the P90 value, and 10% of the chlorophyll values were below the P10 value. To minimise the effect of gaps in the data in the computation of these P90 and P10 climatological values, we imposed a threshold of 25% valid data for the daily climatology. For the 20-year 1998-2017 climatology this means that, for a given pixel and day of the year, at least 5 years must contain valid data for the resulting climatological value to be considered significant. Pixels where the minimum data requirements were met were not considered in further calculations. We compared every valid daily observation over 2020 with the corresponding daily climatology on a pixel-by-pixel basis, to determine if values were above the P90 threshold, below the P10 threshold or within the [P10, P90] range. Values above the P90 threshold or below the P10 were flagged as anomalous. The number of anomalous and total valid observations were stored during this process. We then calculated the percentage of valid anomalous observations (above/below the P90/P10 thresholds) for each pixel, to create percentile anomaly maps in terms of % days per year. Finally, we derived an annual indicator map for eutrophication levels: if 25% of the valid observations for a given pixel and year were above the P90 threshold, the pixel was flagged as eutrophic. Similarly, if 25% of the observations for a given pixel were below the P10 threshold, the pixel was flagged as oligotrophic. '''CONTEXT''' Eutrophication is the process by which an excess of nutrients – mainly phosphorus and nitrogen – in a water body leads to increased growth of plant material in an aquatic body. Anthropogenic activities, such as farming, agriculture, aquaculture and industry, are the main source of nutrient input in problem areas (Jickells, 1998; Schindler, 2006; Galloway et al., 2008). Eutrophication is an issue particularly in coastal regions and areas with restricted water flow, such as lakes and rivers (Howarth and Marino, 2006; Smith, 2003). The impact of eutrophication on aquatic ecosystems is well known: nutrient availability boosts plant growth – particularly algal blooms – resulting in a decrease in water quality (Anderson et al., 2002; Howarth et al.; 2000). This can, in turn, cause death by hypoxia of aquatic organisms (Breitburg et al., 2018), ultimately driving changes in community composition (Van Meerssche et al., 2019). Eutrophication has also been linked to changes in the pH (Cai et al., 2011, Wallace et al. 2014) and depletion of inorganic carbon in the aquatic environment (Balmer and Downing, 2011). Oligotrophication is the opposite of eutrophication, where reduction in some limiting resource leads to a decrease in photosynthesis by aquatic plants, reducing the capacity of the ecosystem to sustain the higher organisms in it. Eutrophication is one of the more long-lasting water quality problems in Europe (OSPAR ICG-EUT, 2017), and is on the forefront of most European Directives on water-protection. Efforts to reduce anthropogenically-induced pollution resulted in the implementation of the Water Framework Directive (WFD) in 2000. '''CMEMS KEY FINDINGS''' Some coastal and shelf waters, especially between 30 and 400N showed active oligotrophication flags for 2020, with some scattered offshore locations within the same latitudinal belt also showing oligotrophication. Eutrophication index is positive only for a small number of coastal locations just north of 40oN, and south of 30oN. In general, the indicator map showed very few areas with active eutrophication flags for 2019 and for 2020. The Third Integrated Report on the Eutrophication Status of the OSPAR Maritime Area (OSPAR ICG-EUT, 2017) reported an improvement from 2008 to 2017 in eutrophication status across offshore and outer coastal waters of the Greater North Sea, with a decrease in the size of coastal problem areas in Denmark, France, Germany, Ireland, Norway and the United Kingdom. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00195

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' Global Ocean- in-situ reprocessed Carbon observations. This product contains observations and gridded files from two up-to-date carbon and biogeochemistry community data products: Surface Ocean Carbon ATlas SOCATv2021 and GLobal Ocean Data Analysis Project GLODAPv2.2021. The SOCATv2021-OBS dataset contains >25 million observations of fugacity of CO2 of the surface global ocean from 1957 to early 2021. The quality control procedures are described in Bakker et al. (2016). These observations form the basis of the gridded products included in SOCATv2020-GRIDDED: monthly, yearly and decadal averages of fCO2 over a 1x1 degree grid over the global ocean, and a 0.25x0.25 degree, monthly average for the coastal ocean. GLODAPv2.2021-OBS contains >1 million observations from individual seawater samples of temperature, salinity, oxygen, nutrients, dissolved inorganic carbon, total alkalinity and pH from 1972 to 2019. These data were subjected to an extensive quality control and bias correction described in Olsen et al. (2020). GLODAPv2-GRIDDED contains global climatologies for temperature, salinity, oxygen, nitrate, phosphate, silicate, dissolved inorganic carbon, total alkalinity and pH over a 1x1 degree horizontal grid and 33 standard depths using the observations from the previous iteration of GLODAP, GLODAPv2. SOCAT and GLODAP are based on community, largely volunteer efforts, and the data providers will appreciate that those who use the data cite the corresponding articles (see References below) in order to support future sustainability of the data products. '''DOI (product) :''' https://doi.org/10.48670/moi-00035

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''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) and near real-time (CMEMS SEALEVEL_GLO_PHY_L4_NRT_OBSERVATIONS_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 long-term mean and trends alone do not give a complete view of the likely changes in position of unstable western boundary current extensions (Kelly et al., 2010). 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. '''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

'''Short description: ''' For the Global Ocean - In-situ observation yearly delivery in delayed mode of Ocean surface currents. '''Detailed description: ''' The In Situ delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for Ocean surface currents. The data are collected from the Surface Drifter Data Assembly Centre (SD-DAC at NOAA AOML) completed by European data provided by EUROGOOS regional systems and national systems by the regional INS TAC components. All surface drifters data have been processed to check for drogue loss. Drogued and undrogued drifting buoy surface ocean currents are provided with a drogue presence flag as well as a wind slippage correction for undrogued buoy. '''Processing information: ''' From the near real time INS TAC product validated on a daily and weekly basis for forecasting purposes, and from the SD-DAC quality controlled dataset a scientifically validated product is created . It s a """"reference product"""" updated on a yearly basis. This product has been processed using a method that checks for drogue loss. Altimeter and wind data have been used to extract the direct wind slippage from the total drifting buoy velocities. The obtained wind slippage values have then been analyzed to identify probable undrogued data among the drifting buoy velocities dataset. A simple procedure has then been applied to produce an updated dataset including a drogue presence flag as well as a wind slippage correction. '''Suitability, Expected type of users / uses: ''' The product is designed to be assimilated into or for validation purposes of operational models operated by ocean forecasting centers for reanalysis purposes or for research community. These users need data aggregated and quality controlled in a reliable and documented manner.