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

Moving 6-year analysis of Water body dissolved inorganic nitrogen in the NorthEast Atlantic for each season: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Every year of the time dimension corresponds to the 6-year centred average of each season. 6-year periods span from 1971/1976 until 2016/2021. Observation data span from 1971 to 2021. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0, 400.0, 500.0, 600.0, 700.0, 800.0, 900.0, 1000.0, 1100.0, 1200.0, 1300.0, 1400.0, 1500.0, 1750.0, 2000.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Descrption of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.4, using GEBCO 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 400km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. Signal-to-noise ratio was fixed to 1 for vertical profiles and 0.1 for time series to account for the redundancy in the time series observations. A logarithmic transformation (DIVAnd.Anam.loglin) was applied to the data prior to the analysis to avoid unrealistic negative values. Background field: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

This visualization product displays the total abundance of marine macro-litter (> 2.5cm) per beach per year from non-MSFD monitoring surveys, research & cleaning operations. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Since the beginning of 2018, data of beach litter have been gathered and processed in the EMODnet Chemistry Marine Litter Database (MLDB). The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols and reference lists used on a European scale. Preliminary processing were necessary to harmonize all the data: - Exclusion of OSPAR 1000 protocol: in order to follow the approach of OSPAR that it is not including these data anymore in the monitoring; - Selection of surveys from non-MSFD monitoring, cleaning and research operations; - Exclusion of beaches without coordinates; - Some categories & some litter types like organic litter, small fragments (paraffin and wax; items > 2.5cm) and pollutants have been removed. The list of selected items is attached to this metadata. This list was created using EU Marine Beach Litter Baselines and EU Threshold Value for Macro Litter on Coastlines from JRC (these two documents are attached to this metadata). - Exclusion of surveys without associated length; - Normalization of survey lengths to 100m & 1 survey / year: in some case, the survey length was not 100m, so in order to be able to compare the abundance of litter from different beaches a normalization is applied using this formula: Number of items (normalized by 100 m) = Number of litter per items x (100 / survey length) Then, this normalized number of items is summed to obtain the total normalized number of litter for each survey. Finally, the median abundance for each beach and year is calculated from these normalized abundances per survey. Percentiles 50, 75, 95 & 99 have been calculated taking into account other sources data for all years. More information is available in the attached documents. Warning: the absence of data on the map doesn't necessarily mean that they don't exist, but that no information has been entered in the Marine Litter Database for this area.

10 years of L-Band remote sensing Sea Surface Salinity (SSS) measurements have proven the capability of satellite SSS to resolve large scale to mesoscale SSS features in tropical to subtropical ocean. In mid to high latitude, L-Band measurements still suffer from large scale and time varying biases. Here, a simple method is proposed to mitigate the large scale and time varying biases. First, in order to estimate these biases, an Optimal Interpolation (OI) using a large correlation scale is used to map SMOS and SMAP L3 products and is compared to equivalent mapping of in situ observations. Then, a second mapping is performed on corrected SSS at scale of SMOS/SMAP resolution (~45 km). This procedure allows to correct and merge both products, and to increase signal to noise ratio of the absolute SSS estimates. Using thermodynamic equation of state (TEOS-10), the resulting L4 SSS product is combined with microwave satellite SST products to produce sea surface density and spiciness, useful to fully characterize the surface ocean water masses. The new L4 SSS products is validated against independent in situ measurements from low to high latitudes. The L4 products exhibits a significant improvement in mid-and high latitude in comparison to the existing SMOS and SMAP L3 products. However, in the Arctic Ocean, L-Band SSS retrieval issues such as sea ice contamination and low sensitivity in cold water are still challenging to improve L-Band SSS data.

Hydrodynamic and sediment dynamic hindcast modelling with a resolution of 2.5 km in the Bay of Biscay and Channel, produced by coupling the hydrodynamic model CROCO with the sediment dynamic module MUSTANG.

Water body phosphate - Monthly Climatology for the European Seas for the period 1960-2020 on the domain from longitude -45.0 to 70.0 degrees East and latitude 24.0 to 83.0 degrees North. Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: The computation was done with the DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.9, using GEBCO 30sec topography for the spatial connectivity of water masses. Horizontal correlation length and vertical correlation length vary spatially depending on the topography and domain. Depth range: 0.0, 5.0, 10.0, 15.0, 20.0, 25.0, 30.0, 35.0, 40.0, 45.0, 50.0, 55.0, 60.0, 65.0, 70.0, 75.0, 80.0, 85.0, 90.0, 95.0, 100.0, 125.0, 150.0, 175.0, 200.0, 225.0, 250.0, 275.0, 300.0, 325.0, 350.0, 375.0, 400.0, 425.0, 450.0, 475.0, 500.0, 550.0, 600.0, 650.0, 700.0, 750.0, 800.0, 850.0, 900.0, 950.0, 1000.0, 1050.0, 1100.0, 1150.0, 1200.0, 1250.0, 1300.0, 1350.0, 1400.0, 1450.0, 1500.0, 1550.0, 1600.0, 1650.0, 1700.0, 1750.0, 1800.0, 1850.0, 1900.0, 1950.0, 2000.0, 2100.0, 2200.0, 2300.0, 2400.0, 2500.0, 2600.0, 2700.0, 2800.0, 2900.0, 3000.0, 3100.0, 3200.0, 3300.0, 3400.0, 3500.0, 3600.0, 3700.0, 3800.0, 3900.0, 4000.0, 4100.0, 4200.0, 4300.0, 4400.0, 4500.0, 4600.0, 4700.0, 4800.0, 4900.0, 5000.0, 5100.0, 5200.0, 5300.0, 5400.0, 5500.0 m. Units: umol/l. The horizontal resolution of the produced DIVAnd analysis is 0.25 degrees.

This product displays for Nickel, positions with values counts that have been measured per matrix for each year and are present in EMODnet regional contaminants aggregated datasets, v2022. The product displays positions for every available year.

This product displays for DDT, DDE, and DDD, median values of the last 6 available years that have been measured per matrix and are present in EMODnet regional contaminants aggregated datasets, v2022. The median values ranges are derived from the following percentiles: 0-25%, 25-75%, 75-90%, >90%. Only "good data" are used, namely data with Quality Flag=1, 2, 6, Q (SeaDataNet Quality Flag schema). For water, only surface values are used (0-15 m), for sediment and biota data at all depths are used.

The main objectives of this dataset is to gather the ocean swells measured by different sensors, including satellite and in-situ sources, that were generated by a given tropical cyclone (TC). This dataset aims at providing characteristics of these swells such as their direction, wavelength (or period) and energy but also the date when they left the influence of the tropical cyclone wind to propagate freely. Wave spectra in tropical cyclones vary strongly per quadrant and provide information about the current and past state of the wave field. However, inside TCs, waves measurements including the wave system direction, energy and wavelength are rare and difficult to obtain with in-situ and remote sensing technics. For this dataset, both moored and drifting buoys are considered as long as they provide wave systems measurements. For the satellite contribution, Synthetic Aperture Radar (SAR) and real aperture radar (RAR) instruments can significantly contribute to the TC-generated waves documentation. Indeed, ocean wave spectra can be derived from modulations of the backscatter in SAR and RAR signal. SAR on board European satellite and in particular the SAR series developed since ERS-1 by ESA and now ESA/Copernicus with Sentinel-1 mission (S-1) are good candidates to provide these ocean waves systems characteristics thanks to the dedicated acquisition mode : the so-called Wave Mode. The wave spectrometer SWIM developed by the French space Agency (CNES) and embedded on the Chinese-French Oceanography SATellite (CFOSAT) has been launched more recently with a new measurement concept relying on a RAR and can certainly complement the S-1 data collection. Although the reasons are different, these two systems are limited for measuring waves generation area within the TC vortex where strong rain rates and wind regimes are observed. Far enough from their source, satellite acquisitions are thus expected to be able to observe these ocean swells during more favorable met-ocean conditions for waves retrieval inversion. As a consequence, our analysis is focused on waves originating from TC but that have been able to propagate far from their source. The analysis of swell measurements far from their area of generation to locate the storm source has been firstly applied to data from one single in-situ wave station (wave energy with frequency and direction) collected 2 miles off shore from San Clemente Island, California and extended to a network of several wave stations in the sixties. More recently, the gathering of swell system observed with SAR far from a storm to characterize the waves properties across the ocean has proven to be efficient in the case of extra-tropical storms. Yet, such analysis is not adapted to Tropical Cyclone whose size is much smaller and currently existing wave datasets do not allow for an accurate monitoring of the tropical cyclones swells. This multi-sensor Level-3 tropical cyclone waves dataset intends to fill this gap and opens for an alternate way of estimating tropical cyclone waves properties over all ocean basins and for all tropical cyclones. This dataset was produced in the frame of the ESA funded Marine Atmosphere eXtreme Satellite Synergy (MAXSS) project. The primary objective of the ESA Marine Atmosphere eXtreme Satellite Synergy (MAXSS) project is to provide guidance and innovative methodologies to maximize the synergetic use of available Earth Observation data (satellite, in situ) to improve understanding about the multi-scale dynamical characteristics of extreme air-sea interaction.

'''DEFINITION''' The OMI_EXTREME_WAVE_MEDSEA_swh_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 significant wave height (swh) measured by in situ buoys. The use of percentiles instead of annual maximum and minimum values, makes this extremes study less affected by individual data measurement errors. The percentiles are temporally averaged, and the spatial evolution is displayed, 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''' Projections on Climate Change foresee a future with a greater frequency of extreme sea states (Stott, 2016; Mitchell, 2006). The damages caused by severe wave storms can be considerable not only in infrastructure and buildings but also in the natural habitat, crops and ecosystems affected by erosion and flooding aggravated by the extreme wave heights. In addition, wave storms strongly hamper the maritime activities, especially in harbours. These extreme phenomena drive complex hydrodynamic processes, whose understanding is paramount for proper infrastructure management, design and maintenance (Goda, 2010). In recent years, there have been several studies searching possible trends in wave conditions focusing on both mean and extreme values of significant wave height using a multi-source approach with model reanalysis information with high variability in the time coverage, satellite altimeter records covering the last 30 years and in situ buoy measured data since the 1980s decade but with sparse information and gaps in the time series (e.g. Dodet et al., 2020; Timmermans et al., 2020; Young & Ribal, 2019). These studies highlight a remarkable interannual, seasonal and spatial variability of wave conditions and suggest that the possible observed trends are not clearly associated with anthropogenic forcing (Hochet et al. 2021, 2023). For the Mediterranean Sea an interesting publication (De Leo et al., 2024) analyses recent studies in this basin showing the variability in the different results and the difficulties to reach a consensus, especially in the mean wave conditions. The only significant conclusion is the positive trend in extreme values for the western Mediterranean Sea and in particular in the Gulf of Lion and in the Tyrrhenian Sea. '''COPERNICUS MARINE SERVICE KEY FINDINGS''' The mean 99th percentiles showed in the area present a range from 1.5-3.5 in the Gibraltar Strait and Alboran Sea with 0.25-0.6 of standard deviation (std), 2-5m in the East coast of the Iberian Peninsula and Balearic Islands with 0.2-0.4m of std, 3-4m in the Aegean Sea with 0.4-0.6m of std to 2-5m in the Gulf of Lyon with 0.3-0.5m of std. Results for this year show a slight negative anomaly in the Gibraltar Strait reaching -0.95m and the Gulf of Lyon (-0.3/-0.7m) slightly over the std in the respective areas, close to zero anomaly in the Aegean Sea (-0.1m) and slight positive or negative anomalies in the East coast of the Iberian Peninsula and Balearic Islands (-0.4/+0.3m) inside the margin of the std. '''DOI (product):''' https://doi.org/10.48670/moi-00263