'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' You can find here the new Mercator Ocean (Toulouse, FR) Global Ocean Ensemble Reanalysis: monthly means of Temperature, Salinity, Currents and Ice variables at 1 degree horizontal resolution for 75 vertical levels, starting from 1993 onward. Global ocean reanalyses are homogeneous 3D gridded descriptions of the physical state of the ocean spanning several decades, produced with a numerical ocean model constrained with data assimilation of satellite and in situ observations. 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-GLORS05 from CMCC (It). ; provided as four different time series of global ocean simulations 3D monthly estimates, post-processed to create this Global Reanalysis Ensemble Product (GREP). Available variables are temperature, salinity, velocities and ice variables. These reanalyses are built to be as close as possible to the observations (i.e. realistic) and in agreement with the model physics. The numerical products available for users are monthly mean averages describing the ocean from surface to bottom (5900 m). '''DOI (product) :''' https://doi.org/10.48670/moi-00023

'''DEFINITION''' The OMI_EXTREME_SL_IBI_slev_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 sea level measured by tide gauges along the coast. The use of percentiles instead of annual maximum and minimum values, makes this extremes study less affected by individual data measurement errors. The annual percentiles referred to annual mean sea level are temporally averaged and their spatial evolution is displayed in the dataset omi_extreme_sl_ibi_slev_mean_and_anomaly_obs, 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''' Sea level (SLEV) is one of the Essential Ocean Variables most affected by climate change. Global mean sea level rise has accelerated since the 1990’s (Abram et al., 2019, Legeais et al., 2020), due to the increase of ocean temperature and mass volume caused by land ice melting (WCRP, 2018). Basin scale oceanographic and meteorological features lead to regional variations of this trend that combined with changes in the frequency and intensity of storms could also rise extreme sea levels up to one meter by the end of the century (Vousdoukas et al., 2020, Tebaldi et al., 2021). This will significantly increase coastal vulnerability to storms, with important consequences on the extent of flooding events, coastal erosion and damage to infrastructures caused by waves (Boumis et al., 2023). The increase in extreme sea levels over recent decades is, therefore, primarily due to the rise in mean sea level. Note, however, that the methodology used to compute this OMI removes the annual 50th percentile, thereby discarding the mean sea level trend to isolate changes in storminess. The Iberian Biscay Ireland region shows positive sea level trend modulated by decadal-to-multidecadal variations driven by ocean dynamics and superposed to the long-term trend (Chafik et al., 2019). '''COPERNICUS MARINE SERVICE KEY FINDINGS''' The completeness index criteria is fulfilled by 62 stations in 2023, five more than those available in 2022 (57), recently added to the multi-year product INSITU_GLO_PHY_SSH_DISCRETE_MY_013_053. The mean 99th percentiles reflect the great tide spatial variability around the UK and the north of France. Minimum values are observed in the Irish eastern coast (e.g.: 0.66 m above mean sea level in Arklow Harbour) and the Canary Islands (e.g.: 0.93 and 0.96 m above mean sea level in Gomera and Hierro, respectively). Maximum values are observed in the Bristol Channel (e.g.: 6.25 and 5.78 m above mean sea level in Newport and Hinkley, respectively), and in the English Channel (e.g.: 5.16 m above mean sea level in St. Helier). The annual 99th percentiles standard deviation reflects the south-north increase of storminess, ranging between 1-3 cm in the Canary Islands to 15 cm in Hinkley (Bristol Channel). Negative or close to zero anomalies of 2023 99th percentile prevail throughout the region this year, reaching < -20 cm in several stations of the UK western coast and the English Channel (e.g.: -22 cm in Newport; -21 cm in St.Helier). Significantly positive anomaly of 2023 99th percentile is only found in Arcklow Harbour, in the eastern Irish coast. '''DOI (product):''' https://doi.org/10.48670/moi-00253

'''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 Mediterranean water mass formation rates are evaluated in 4 areas as defined in the Ocean State Report issue 2 (OSR2, von Schuckmann et al., 2018) section 3.4 (Simoncelli and Pinardi, 2018): (1) the Gulf of Lions for the Western Mediterranean Deep Waters (WMDW); (2) the Southern Adriatic Sea Pit for the Eastern Mediterranean Deep Waters (EMDW); (3) the Cretan Sea for Cretan Intermediate Waters (CIW) and Cretan Deep Waters (CDW); (4) the Rhodes Gyre, the area of formation of the so-called Levantine Intermediate Waters (LIW) and Levantine Deep Waters (LDW). Annual water mass formation rates have been computed using daily mixed layer depth estimates (density criteria Δσ = 0.01 kg/m3, 10 m reference level) considering the annual maximum volume of water above mixed layer depth with potential density within or higher the specific thresholds specified in Table 1 then divided by seconds per year. Annual mean values are provided using the Mediterranean 1/24o eddy resolving reanalysis (Escudier et al. 2020, 2021). Time spans from 1987 to the year preceding the current one [-1Y], operationally extended yearly. '''CONTEXT''' The formation of intermediate and deep water masses is one of the most important processes occurring in the Mediterranean Sea, being a component of its general overturning circulation. This circulation varies at interannual and multidecadal time scales and it is composed of an upper zonal cell (Zonal Overturning Circulation) and two main meridional cells in the western and eastern Mediterranean (Pinardi and Masetti 2000). The objective is to monitor the main water mass formation events using the eddy resolving Mediterranean Sea Reanalysis (MEDSEA_MULTIYEAR_PHY_006_004, Escudier et al. 2020, 2021) and considering Pinardi et al. (2015) and Simoncelli and Pinardi (2018) as references for the methodology. The Mediterranean Sea Reanalysis can reproduce both Eastern Mediterranean Transient and Western Mediterranean Transition phenomena and catches the principal water mass formation events reported in the literature. This will permit constant monitoring of the open ocean deep convection process in the Mediterranean Sea and a better understanding of the multiple drivers of the general overturning circulation at interannual and multidecadal time scales. Deep and intermediate water formation events reveal themselves by a deep mixed layer depth distribution in four Mediterranean areas: Gulf of Lions, Southern Adriatic Sea Pit, Cretan Sea and Rhodes Gyre. '''KEY FINDINGS''' The Western Mediterranean Deep Water (WMDW) formation events in the Gulf of Lion appear to be larger after 1999 consistently with Schroeder et al. (2006, 2008) related to the Eastern Mediterranean Transient event. This modification of WMDW after 2005 has been called Western Mediterranean Transition. WMDW formation events are consistent with Somot et al. (2016) and the event in 2009 is also reported in Houpert et al. (2016). The Eastern Mediterranean Deep Water (EMDW) formation in the Southern Adriatic Pit region displays a period of water mass formation between 1988 and 1993, in agreement with Pinardi et al. (2015), in 1996, 1999 and 2000 as documented by Manca et al. (2002). Weak deep water formation in winter 2006 is confirmed by observations in Vilibić and Šantić (2008). An intense deep water formation event is detected in 2012-2013 (Gačić et al., 2014). Last years are characterized by large events starting from 2017 (Mihanovic et al., 2021). Cretan Intermediate Water formation rates present larger peaks between 1989 and 1993 with the ones in 1992 and 1993 composing the Eastern Mediterranean Transient phenomena. The Cretan Deep Water formed in 1992 and 1993 is characterized by the highest densities of the entire period in accordance with Velaoras et al. (2014). The Levantine Deep Water formation rate in the Rhode Gyre region presents the largest values between 1992 and 1993 in agreement with Kontoyiannis et al. (1999). '''DOI (product):''' https://doi.org/10.48670/mds-00318

'''DEFINITION''' The global annual chlorophyll anomaly is computed by subtracting a reference climatology (1997-2014) from the annual chlorophyll mean, on a pixel-by-pixel basis and in log10 space. Both the annual mean and the climatology are computed employing ESA Ocean Colour Climate Change Initiative (ESA OC-CCI, Sathyendranath et al., 2018a) global products (i.e. using the standard OC-CCI chlorophyll algorithms, OCI) as distributed by CMEMS. '''CONTEXT''' Phytoplankton – and chlorophyll concentration as a proxy for phytoplankton – respond rapidly to changes in their physical environment. Some of those changes are seasonal and are determined by light and nutrient availability (Racault et al., 2012). By comparing annual mean values to a climatology, we effectively remove the seasonal signal, while retaining information on potential events during the year. Chlorophyll anomalies can be correlated to climate indexes in particular regions, such as the ENSO index in the equatorial Pacific (Behrenfeld et al. 2006; Racault et al., 2012) and the IOD index in the Indian Ocean (Brewin et al., 2012). It is important to study chlorophyll anomalies in consonance with sea surface temperature and sea level anomalies, as increases in chlorophyll are generally consistent with decreases in SST and sea level anomalies, suggesting an increase in mixing and vertical nutrient transport (von Schuckmann et al., 2016). '''CMEMS KEY FINDINGS''' The average global chlorophyll anomaly 2019 is -0.02 log10(mg m-3), with a maximum value of 1.7 log10(mg m-3) and a minimum value of -3.2 log10(mg m-3). That is to say that, in average, the annual 2019 mean value is slightly lower (96%) than the 1997-2014 climatological value. The positive signals reported in 2016 and 2017 (Sathyendranath et al., 2018b) in the southern Pacific Ocean could still be observed in the 2019 map, while the significant negative anomalies in the tropical waters of the northern Pacific Ocean were also detected to a lesser extent. Areas showing a change of anomaly sign from 2019 include the southern coast of Japan (no anomaly to positive) and the tropical Atlantic (anomalies close to zero for 2019). A marked increase in chlorophyll concentration was observed during 2019 in the Great Australian Bight, while negative anomalies became stronger in the Guatemala Basin and the region south of the Gulf of Guinea and, with values of chlorophyll reaching as low as 30% of the climatological value (anomaly < -0.5 log10(mg m-3)). The persistent positive anomalies in the higher latitudes of the North Atlantic (> 40°) match the cooling observed in the 2018 and previous years SST anomaly maps.

'''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 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). 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-2023 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 -5 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

'''DEFINITION''' Volume transport across lines are obtained by integrating the volume 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 Sverdrup (Sv). '''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. At Drake Passage, the multi-product approach (product no. 2.4.1) is larger than the value (130 Sv) of Lumpkin and Speer (2007), but smaller than the new observational based results of Colin de Verdière and Ollitrault, (2016) (175 Sv) and Donohue (2017) (173.3 Sv). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00247