'''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:''' For the Global Ocean- Gridded objective analysis fields of temperature and salinity using profiles from the in-situ near real time database are produced monthly. Objective analysis is based on a statistical estimation method that allows presenting a synthesis and a validation of the dataset, providing a support for localized experience (cruises), providing a validation source for operational models, observing seasonal cycle and inter-annual variability. '''DOI (product) :''' https://doi.org/10.48670/moi-00037

'''Short description:''' In wavenumber spectra, the 1hz measurement error is the noise level estimated as the mean value of energy at high wavenumbers (below ~20km in term of wavelength). The 1hz noise level spatial distribution follows the instrumental white-noise linked to the Surface Wave Height but also connections with the backscatter coefficient. The full understanding of this hump of spectral energy (Dibarboure et al., 2013, Investigating short wavelength correlated errors on low-resolution mode altimetry, OSTST 2013 presentation) still remain to be achieved and overcome with new retracking, new editing strategy or new technology. '''DOI (product) :''' https://doi.org/10.48670/moi-00144

"'Short description: ''' Global Ocean - This delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for ocean surface and subsurface currents. Current data from 5 different types of instruments are distributed: * The drifter's near-surface velocities computed from their position measurements. In addition, a wind slippage correction is provided from 1993. Information on the presence of the drogue of the drifters is also provided. * The near-surface zonal and meridional total velocities, and near-surface radial velocities, measured by High Frequency (HF) radars that are part of the European HF radar Network. These data are delivered together with standard deviation of near-surface zonal and meridional raw velocities, Geometrical Dilution of Precision (GDOP), quality flags and metadata. * The zonal and meridional velocities, at parking depth (mostly around 1000m) and at the surface, calculated along the trajectories of the floats which are part of the Argo Program. * The velocity profiles within the water column coming from Acoustic Doppler Current Profiler (vessel mounted ADCP, Moored ADCP, saildrones) platforms * The near-surface and subsurface velocities calculated from gliders (autonomous underwater vehicle) trajectories '''DOI (product) :''' https://doi.org/10.17882/86236

'''Short description:''' The IBI-MFC provides the biogeochemical multi-year (non assimilative) product for the Iberia-Biscay-Ireland region starting in 01/01/1993, extended every year to use available reprocessed upstream data and regularly updated on a monthly basis to cover the period up to month M-4 using an interim processing system. The model system is designed, developed and run by Mercator Ocean International, while the operational product post-processing and interim processing system are run by NOW Systems with the support of CESGA supercomputing centre. The biogeochemical model PISCES is run simultaneously with the ocean physical NEMO model, generating products at 1/36° horizontal resolution. The PISCES model is able to simulate the first levels of the marine food web, from nutrients up to mesozooplankton and it has 24 state variables. The product provides daily, monthly and yearly averages of the main biogeochemical variables. Additionally, climatological parameters (monthly mean and standard deviation) of these variables for the period 1993-2016 are delivered. '''DOI (Product)''': https://doi.org/10.48670/moi-00028

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.

'''Short description:''' The operational global ocean analysis and forecast system of Météo-France with a resolution of 1/12 degree is providing daily analyses and 10 days forecasts for the global ocean sea surface waves. This product includes 3-hourly instantaneous fields of integrated wave parameters from the total spectrum (significant height, period, direction, Stokes drift,...etc), as well as the following partitions: the wind wave, the primary and secondary swell waves. The global wave system of Météo-France is based on the wave model MFWAM which is a third generation wave model. MFWAM uses the computing code ECWAM-IFS-38R2 with a dissipation terms developed by Ardhuin et al. (2010). The model MFWAM was upgraded on november 2014 thanks to improvements obtained from the european research project « my wave » (Janssen et al. 2014). The model mean bathymetry is generated by using 2-minute gridded global topography data ETOPO2/NOAA. Native model grid is irregular with decreasing distance in the latitudinal direction close to the poles. At the equator the distance in the latitudinal direction is more or less fixed with grid size 1/10°. The operational model MFWAM is driven by 6-hourly analysis and 3-hourly forecasted winds from the IFS-ECMWF atmospheric system. The wave spectrum is discretized in 24 directions and 30 frequencies starting from 0.035 Hz to 0.58 Hz. The model MFWAM uses the assimilation of altimeters with a time step of 6 hours. The global wave system provides analysis 4 times a day, and a forecast of 10 days at 0:00 UTC. The wave model MFWAM uses the partitioning to split the swell spectrum in primary and secondary swells. '''DOI (product) :''' https://doi.org/10.48670/moi-00017

'''DEFINITION''' Ocean heat content (OHC) is defined here as the deviation from a reference period (1993-20210) and is closely proportional to the average temperature change from z1 = 0 m to z2 = 2000 m depth: With a reference density of ρ0 = 1030 kgm-3 and a specific heat capacity of cp = 3980 J/kg°C (e.g. von Schuckmann et al., 2009) Averaged time series for ocean heat content and their error bars are calculated for the Iberia-Biscay-Ireland region (26°N, 56°N; 19°W, 5°E). This OMI is computed using IBI-MYP, GLO-MYP reanalysis and CORA, ARMOR data from observations which provide temperatures. Where the CMEMS product for each acronym is: • IBI-MYP: IBI_MULTIYEAR_PHY_005_002 (Reanalysis) • GLO-MYP: GLOBAL_REANALYSIS_PHY_001_031 (Reanalysis) • CORA: INSITU_GLO_TS_OA_REP_OBSERVATIONS_013_002_b (Observations) • ARMOR: MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012 (Reprocessed observations) The figure comprises ensemble mean (blue line) and the ensemble spread (grey shaded). Details on the product are given in the corresponding PUM for this OMI as well as the CMEMS Ocean State Report: von Schuckmann et al., 2016; von Schuckmann et al., 2018. '''CONTEXT''' Change in OHC is a key player in ocean-atmosphere interactions and sea level change (WCRP, 2018) and can impact marine ecosystems and human livelihoods (IPCC, 2019). Additionally, OHC is one of the six Global Climate Indicators recommended by the World Meterological Organisation (WMO, 2017). In the last decades, the upper North Atlantic Ocean experienced a reversal of climatic trends for temperature and salinity. While the period 1990-2004 is characterized by decadal-scale ocean warming, the period 2005-2014 shows a substantial cooling and freshening. Such variations are discussed to be linked to ocean internal dynamics, and air-sea interactions (Fox-Kemper et al., 2021; Collins et al., 2019; Robson et al 2016). Together with changes linked to the connectivity between the North Atlantic Ocean and the Mediterranean Sea (Masina et al., 2022), these variations affect the temporal evolution of regional ocean heat content in the IBI region. Recent studies (de Pascual-Collar et al., 2023) highlight the key role that subsurface water masses play in the OHC trends in the IBI region. These studies conclude that the vertically integrated trend is the result of different trends (both positive and negative) contributing at different layers. Therefore, the lack of representativeness of the OHC trends in the surface-intermediate waters (from 0 to 1000 m) causes the trends in intermediate and deep waters (from 1000 m to 2000 m) to be masked when they are calculated by integrating the upper layers of the ocean (from surface down to 2000 m). '''CMEMS KEY FINDINGS''' The ensemble mean OHC anomaly time series over the Iberia-Biscay-Ireland region are dominated by strong year-to-year variations, and an ocean warming trend of 0.41±0.4 W/m2 is barely significant. '''DOI (product):''' https://doi.org/10.48670/mds-00316

'''DEFINITION''' Ocean acidification is quantified by decreases in pH, which is a measure of acidity: a decrease in pH value means an increase in acidity, that is, acidification. The observed decrease in ocean pH resulting from increasing concentrations of CO2 is an important indicator of global change. The estimate of global mean pH builds on a reconstruction methodology, * Obtain values for alkalinity based on the so called “locally interpolated alkalinity regression (LIAR)” method after Carter et al., 2016; 2018. * Build on surface ocean partial pressure of carbon dioxide (CMEMS product: MULTIOBS_GLO_BIO_CARBON_SURFACE_REP_015_008) obtained from an ensemble of Feed-Forward Neural Networks (Chau et al. 2022) which exploit sampling data gathered in the Surface Ocean CO2 Atlas (SOCAT) (https://www.socat.info/) * Derive a gridded field of ocean surface pH based on the van Heuven et al., (2011) CO2 system calculations using reconstructed pCO2 (MULTIOBS_GLO_BIO_CARBON_SURFACE_REP_015_008) and alkalinity. The global mean average of pH at yearly time steps is then calculated from the gridded ocean surface pH field. It is expressed in pH unit on total hydrogen ion scale. In the figure, the amplitude of the uncertainty (1σ ) of yearly mean surface sea water pH varies at a range of (0.0023, 0.0029) pH unit (see Quality Information Document for more details). The trend and uncertainty estimates amount to -0.0017±0.0004e-1 pH units per year. The indicator is derived from in situ observations of CO2 fugacity (SOCAT data base, www.socat.info, Bakker et al., 2016). These observations are still sparse in space and time. Monitoring pH at higher space and time resolutions, as well as in coastal regions will require a denser network of observations and preferably direct pH measurements. A full discussion regarding this OMI can be found in section 2.10 of the Ocean State Report 4 (Gehlen et al., 2020). '''CONTEXT''' The decrease in surface ocean pH is a direct consequence of the uptake by the ocean of carbon dioxide. It is referred to as ocean acidification. The International Panel on Climate Change (IPCC) Workshop on Impacts of Ocean Acidification on Marine Biology and Ecosystems (2011) defined Ocean Acidification as “a reduction in the pH of the ocean over an extended period, typically decades or longer, which is caused primarily by uptake of carbon dioxide from the atmosphere, but can also be caused by other chemical additions or subtractions from the ocean”. The pH of contemporary surface ocean waters is already 0.1 lower than at pre-industrial times and an additional decrease by 0.33 pH units is projected over the 21st century in response to the high concentration pathway RCP8.5 (Bopp et al., 2013). Ocean acidification will put marine ecosystems at risk (e.g. Orr et al., 2005; Gehlen et al., 2011; Kroeker et al., 2013). The monitoring of surface ocean pH has become a focus of many international scientific initiatives (http://goa-on.org/) and constitutes one target for SDG14 (https://sustainabledevelopment.un.org/sdg14). '''CMEMS KEY FINDINGS''' Since the year 1985, global ocean surface pH is decreasing at a rate of -0.0017±0.0004e-1 per year. '''DOI (product):''' https://doi.org/10.48670/moi-00224

'''Short description: ''' For the '''Global''' Ocean '''Satellite Observations''', ACRI-ST company (Sophia Antipolis, France) is providing '''Bio-Geo-Chemical (BGC)''' products based on the '''Copernicus-GlobColour''' processor. * Upstreams: SeaWiFS, MODIS, MERIS, VIIRS-SNPP & JPSS1, OLCI-S3A & S3B for the '''""multi""''' products, and S3A & S3B only for the '''""olci""''' products. * Variables: Chlorophyll-a ('''CHL'''), Gradient of Chlorophyll-a ('''CHL_gradient'''), Phytoplankton Functional types and sizes ('''PFT'''), Suspended Matter ('''SPM'''), Secchi Transparency Depth ('''ZSD'''), Diffuse Attenuation ('''KD490'''), Particulate Backscattering ('''BBP'''), Absorption Coef. ('''CDM''') and Reflectance ('''RRS'''). * Temporal resolutions: '''daily'''. * Spatial resolutions: '''4 km''' and a finer resolution based on olci '''300 meters''' inputs. * 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 the '''Copernicus-GlobColour''' products in the catalogue, use the search keyword '''""GlobColour""'''. '''DOI (product) :''' https://doi.org/10.48670/moi-00280