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

'''Short description:''' You can find here the biogeochemistry non assimilative hindcast simulation GLOBAL_REANALYSIS_BIO_001_018 at 1/4° over period 1998 - 2016. Outputs are delivered as monthly mean files with Netcdf format (CF/COARDS 1.5 convention) on the native tripolar grid (ORCA025) at ¼° resolution with 75 vertical levels. This simulation is based on the PISCES biogeochemical model. It is forced offline at a daily frequency by the equivalent of the GLOBAL-REANALYSIS-PHYS-001-009 physics product but without data assimilation. '''Detailed description: ''' There are 8 different datasets: * dataset-global-nahindcast-bio-001-018-no3 containing : nitrate concentration * dataset-global-nahindcast-bio-001-018-po4 containing : phosphate concentration * dataset-global-nahindcast-bio-001-018-si containing : silicate concentration * dataset-global-nahindcast-bio-001-018-o2 containing : dissolved oxygen concentration * dataset-global-nahindcast-bio-001-018-fe containing : iron concentration * dataset-global-nahindcast-bio-001-018-chl containing : chlorophyll concentration * dataset-global-nahindcast-bio-001-018-phyc containing : carbon phytoplankton biomass * dataset-global-nahindcast-bio-001-018-pp containing : primary production The horizontal grid is the standard ORCA025 tri-polar grid (1440 x 1021 grid points). The three poles are located over Antarctic, Central Asia and North Canada. The ¼ degree resolution corresponds to the equator. The vertical grid has 75 levels, with a resolution of 1 meter near the surface and 200 meters in the deep ocean.Biogeochemical and physical simulations start at rest (cold start) in December 1991. The spin-up period consists of 5 years of interannual simulation between 1992 and 1997. The simulation period covers the ocean color era (1998 – 2016).The biogeochemical model used is PISCES (Aumont, in prep). It is a model of intermediate complexity designed for global ocean applications (Aumont and Bopp, 2006) and is part of NEMO modeling platform. It has 24 prognostic variables and simulates biogeochemical cycles of oxygen, carbon and the main nutrients controlling phytoplankton growth (nitrate, ammonium, phosphate, silicic acid and iron). The model distinguishes four plankton functional types based on size: two phytoplankton groups <nowiki>(small = nanophytoplankton and large = diatoms)</nowiki> and two zooplankton groups <nowiki>(small = microzooplankton and large = mesozooplankton).</nowiki>Prognostic variables of phytoplankton are total biomass in C, Fe, Si (for diatoms) and chlorophyll and hence the Fe/C, Si/C, Chl/C ratios are variable. For zooplankton, all these ratios are constant and total biomass in C is the only prognostic variable. The bacterial pool is not modeled explicitly. PISCES distinguishes three non-living pools for organic carbon: small particulate organic carbon, big particulate organic carbon and semi-labile dissolved organic carbon. While the C/N/P composition of dissolved and particulate matter is tied to Redfield stoichiometry, the iron, silicon and carbonate contents of the particles are computed prognostically. Next to the three organic detrital pools, carbonate and biogenic siliceous particles are modeled. Besides, the model simulates dissolved inorganic carbon and total alkalinity. In PISCES, phosphate and nitrate + ammonium are linked by constant Redfield ratio <nowiki>(C/N/P = 122/16/1)</nowiki>, but cycles of phosphorus and nitrogen are decoupled by nitrogen fixation and denitrification. Biogeochemical model PISCES (NEMO3.5) is forced offline by daily fields of the physical model NEMO (OPA module in the NEMO platform) without any assimilation of physical data. The main features of this dynamical ocean are: * NEMO 3.1 * Atmospheric forcings from 3-hourly ERA-Interim reanalysis products, CORE bulk formulation * Vertical diffusivity coefficient is computed by solving the TKE equation * Tidal mixing is parameterized according to the works of Bessières et al. (2008) and Koch-Larrouy et al, (2006). * Sea-Ice model: LIM2 with the Elastic-Viscous-Plastic rheology * Initial conditions: Levitus 98 climatology for temperature and salinity, patched with PHC2.1 for the Arctic regions, and Medatlas for the Mediterranean Sea. A special treatment is done on vertical diffusivity coefficient (Kz): the daily mean is done on Log10(Kz) after a filtering of enhanced convection (Kz increased artificially to 10 m2.s-1 when the water column is unstable). The purpose of this Log10 is to average the orders of magnitudes and to give more weight to small values of vertical diffusivity. The atmospheric forcing fields are daily averages from ERA-Interim reanalysis product (CORE bulk formulation). Boundary fluxes account for nutrient supply from three different sources: Atmospheric deposition (Aumont et al., 2008), rivers for nutrients, dissolved inorganic carbon and alkalinity (Ludwig et al., 1996; Mayorga et al., 2010) and inputs of Fe from marine sediments. Nutrient and freshwater inflows by rivers are colocalized. River and dust inputs are balanced with sediment trapping of NO3, Si and Carbon. An annual and global value of atmospheric carbon dioxide is imposed at sea surface.

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

'''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'''                For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' You can find here the CMEMS Global Ocean Ensemble Reanalysis product at ¼ degree resolution : monthly means of Temperature, Salinity, Currents and Ice variables for 75 vertical levels, starting from 1993 onward. Global ocean reanalyses are homogeneous 3D gridded descriptions of the physical state of the ocean covering several decades, produced with a numerical ocean model constrained with data assimilation of satellite and in situ observations. These reanalyses are built to be as close as possible to the observations (i.e. realistic) and in agreement with the model physics The multi-model ensemble approach allows uncertainties or error bars in the ocean state to be estimated. The ensemble mean may even provide for certain regions and/or periods a more reliable estimate than any individual reanalysis product. The four reanalyses, used to create the ensemble, covering “altimetric era” period (starting from 1st of January 1993) during which altimeter altimetry data observations are available: * GLORYS2V4 from Mercator Ocean (Fr); * ORAS5 from ECMWF; * GloSea5 from Met Office (UK); * and C-GLORSv7 from CMCC (It); These four products provided four different time series of global ocean simulations 3D monthly estimates. All numerical products available for users are monthly or daily mean averages describing the ocean. '''DOI (product) :''' https://doi.org/10.48670/moi-00024

Sea Surface Salinity (SSS) is an essential climate variable that requires long term in situ observation. The French SSS Observation Service (SSS-OS) manages a network of Voluntary Observing Ships equipped with thermosalinographs (TSG) since 2002. The network is global though more concentrated in the tropical Pacific and North Atlantic oceanic basins. The acquisition system is autonomous with real time transmission and is regularly serviced at harbor calls. The high resolution data retrieved from the acquisition system during ship calls is processed through a dedicated software (freely available) for attribution of data quality flags by visual inspection, and correction of TSG time series by comparison with climatology, onboard daily water samples and collocated Argo data. Details can be found in the reference below. The validated delayed time data collected from TSG, together with some bucket samples mostly collected before 2002, are made available for educational and research purposes through an interactive web interface.

The binned Sea Surface Salinity, Temperature and Density data set covers regularly sampled ship-of-opportunity lines. It is based on data collected from 1993 to 2018 from Voluntary Observing Ships subsequently validated. This monthly product is binned with monthly total values deviations and then deviations from climatology in each bin, with additional 1-2-1 time filter applied on the averages. Along B-AX01, some gaps were filled with additional data located one degree north or south of the grid boxes.

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

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