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

'''Short description:''' The Med MFC physical multiyear product is generated by a numerical system composed of an hydrodynamic model, supplied by the Nucleous for European Modelling of the Ocean (NEMO) and a variational data assimilation scheme (OceanVAR) for temperature and salinity vertical profiles and satellite Sea Level Anomaly along track data. The model horizontal grid resolution is 1/24˚ (ca. 4-5 km) and the unevenly spaced vertical levels are 141. The datasets are extended every year as well as on a monthly basis through one-month extensions in interim mode, reaching one month before present. ''DOI (Product)'': https://doi.org/10.48670/mds-00375 ''DOI (Interim dataset)'': https://doi.org/10.48670/mds-00375

'''This product has been archived''' '''DEFINITION''' Estimates of Ocean Heat Content (OHC) are obtained from integrated differences of the measured temperature and a climatology along a vertical profile in the ocean (von Schuckmann et al., 2018). The regional OHC values are then averaged from 60°S-60°N aiming i) to obtain the mean OHC as expressed in Joules per meter square (J/m2) to monitor the large-scale variability and change. ii) to monitor the amount of energy in the form of heat stored in the ocean (i.e. the change of OHC in time), expressed in Watt per square meter (W/m2). Ocean heat content is one of the six Global Climate Indicators recommended by the World Meterological Organisation for Sustainable Development Goal 13 implementation (WMO, 2017). '''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 ocean shapes our perspectives for the future (von Schuckmann et al., 2020). Variations in OHC can induce changes in ocean stratification, currents, sea ice and ice shelfs (IPCC, 2019; 2021); they set time scales and dominate Earth system adjustments to climate variability and change (Hansen et al., 2011); they are a key player in ocean-atmosphere interactions and sea level change (WCRP, 2018) and they can impact marine ecosystems and human livelihoods (IPCC, 2019). '''CMEMS KEY FINDINGS''' Since the year 2005, the upper (0-2000m) near-global (60°S-60°N) ocean warms at a rate of 1.0 ± 0.1 W/m2. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00235

'''Short description:''' Arctic L3 sea ice product providing concentration, stage-of-development and floe size information retrieved from Sentinel-1 and RCM SAR imagery and GCOM-W AMSR2 microwave radiometer data using a deep learning algorithm and delivered on a 0.5 km grid. '''DOI (product) :''' https://doi.org/10.48670/mds-00343

'''DEFINITION''' Estimates of Ocean Heat Content (OHC) are obtained from integrated differences of the measured temperature and a climatology along a vertical profile in the ocean (von Schuckmann et al., 2018). 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). Additionally, the time series based on the method of von Schuckmann and Le Traon (2011) has been added. The regional OHC values are then averaged from 60°S-60°N aiming i) to obtain the mean OHC as expressed in Joules per meter square (J/m2) to monitor the large-scale variability and change. ii) to monitor the amount of energy in the form of heat stored in the ocean (i.e. the change of OHC in time), expressed in Watt per square meter (W/m2). Ocean heat content is one of the six Global Climate Indicators recommended by the World Meterological Organisation for Sustainable Development Goal 13 implementation (WMO, 2017). '''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 ocean shapes our perspectives for the future (von Schuckmann et al., 2020). Variations in OHC can induce changes in ocean stratification, currents, sea ice and ice shelfs (IPCC, 2019; 2021); they set time scales and dominate Earth system adjustments to climate variability and change (Hansen et al., 2011); they are a key player in ocean-atmosphere interactions and sea level change (WCRP, 2018) and they can impact marine ecosystems and human livelihoods (IPCC, 2019). '''CMEMS KEY FINDINGS''' Since the year 2005, the upper (0-700m) near-global (60°S-60°N) ocean warms at a rate of 0.6 ± 0.1 W/m2. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00234

'''DEFINITION''' The OMI_EXTREME_SST_NORTHWESTSHELF_sst_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 surface temperature measured by in situ buoys at depths between 0 and 5 meters. 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''' Sea surface temperature (SST) is one of the essential ocean variables affected by climate change (mean SST trends, SST spatial and interannual variability, and extreme events). In Europe, several studies show warming trends in mean SST for the last years (von Schuckmann, 2016; IPCC, 2021, 2022). An exception seems to be the North Atlantic, where, in contrast, anomalous cold conditions have been observed since 2014 (Mulet et al., 2018; Dubois et al. 2018; IPCC 2021, 2022). Extremes may have a stronger direct influence in population dynamics and biodiversity. According to Alexander et al. 2018 the observed warming trend will continue during the 21st Century and this can result in exceptionally large warm extremes. Monitoring the evolution of sea surface temperature extremes is, therefore, crucial. The North-West Self area comprises part of the North Atlantic, where this refreshing trend has been observed, and the North Sea, where a warming trend has been taking place in the last three decades (e.g. Høyer and Karagali, 2016). '''COPERNICUS MARINE SERVICE KEY FINDINGS''' The mean 99th percentiles showed in the area present a range from 14-15ºC in the North of the British Isles, 16-19ºC in the West of the North Sea to 19-20ºC in the Helgoland Bight. The standard deviation ranges from 0.7-0.8ºC in the North of the British Isles, 0.6-2ºC in the West of the North Sea to 0.8-3ºC in in the Helgoland Bight. Results for this year show positive moderate anomalies (+0.3/+1.0ºC) in all the positions except in one station in the West of the Noth Sea where the anomaly is negative (-0.3ºC), all of them inside the standard deviation margin. '''DOI (product):''' https://doi.org/10.48670/moi-00274

'''Short description:''' This product consists of 3D fields of Particulate Organic Carbon (POC), Particulate Backscattering coefficient (bbp), Chlorophyll-a concentration (Chla), Downwelling Photosynthetic Available Radiation (PAR) and downwelling irradiance, at 0.25°x0.25° resolution from the surface to 1000 m. A neural network estimates the vertical distribution of Chla and bbp from surface ocean color measurements with hydrological properties and additional drivers. The SOCA-light models is used to integrate light. '''DOI (product):''' https://doi.org/10.48670/moi-00046

'''DEFINITION''' The OMI_EXTREME_WAVE_IBI_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). In the North Atlantic, the mean wave height shows some weak trends not very statistically significant. Young & Ribal (2019) found a mostly positive weak trend in the European Coasts while Timmermans et al. (2020) showed a weak negative trend in high latitudes, including the North Sea and even more intense in the Norwegian Sea. For extreme values, some authors have found a clearer positive trend in high percentiles (90th-99th) (Young, 2011; Young & Ribal, 2019). '''COPERNICUS MARINE SERVICE KEY FINDINGS''' The mean 99th percentiles showed in the area present a wide range from 2-3.5m in the Canary Island with 0.1-0.3 m of standard deviation (std), 3.5m in the Gulf of Cadiz with 0.5m of std, 3-6m in the English Channel and the Irish Sea with 0.5-0.6m of std, 4-7m in the Bay of Biscay with 0.4-0.9m of std to 8-10m in the West of the British Isles with 0.7-1.4m of std. Results for this year show slight negative anomalies in the Canary Island (-0.4/0.0m) and in the Gulf of Cadiz (-0.8m) barely out of the standard deviation range in both areas, slight positive or negative anomalies in the West of the British Isles (-0.6/+0.4m) and in the English Channel and the Irish Sea (-0.6/+0.3m) but inside the range of the standard deviation and a general positive anomaly in the Bay of Biscay reaching +1.0m but close to the limit of the standard deviation. '''DOI (product):''' https://doi.org/10.48670/moi-00250

'''DEFINITION''' The time series are derived from the regional chlorophyll reprocessed (MY) product as distributed by CMEMS (OCEANCOLOUR_MED_BGC_L3_NRT_009_141). This dataset, derived from multi-sensor (SeaStar-SeaWiFS, AQUA-MODIS, NOAA20-VIIRS, NPP-VIIRS, Envisat-MERIS and Sentinel3-OLCI) Rrs spectra produced by CNR using an in-house processing chain, is obtained by means of the Mediterranean Ocean Colour regional algorithms: an updated version of the MedOC4 (Case 1 (off-shore) waters, Volpe et al., 2019, with new coefficients) and AD4 (Case 2 (coastal) waters, Berthon and Zibordi, 2004). The processing chain and the techniques used for algorithms merging are detailed in Colella et al. (2023). Monthly regional mean values are calculated by performing the average of 2D monthly mean (weighted by pixel area) over the region of interest. The deseasonalized time series is obtained by applying the X-11 seasonal adjustment methodology on the original time series as described in Colella et al. (2016), and then the Mann-Kendall test (Mann, 1945; Kendall, 1975) and Sens’s method (Sen, 1968) are subsequently applied to obtain the magnitude of trend. This OMI has been introduced since the 2nd issue of Ocean State Report in 2017. '''CONTEXT''' Phytoplankton and chlorophyll concentration as a proxy for phytoplankton respond rapidly to changes in environmental conditions, such as light, temperature, nutrients and mixing (Colella et al. 2016). The character of the response depends on the nature of the change drivers, and ranges from seasonal cycles to decadal oscillations (Basterretxea et al. 2018). Therefore, it is of critical importance to monitor chlorophyll concentration at multiple temporal and spatial scales, in order to be able to separate potential long-term climate signals from natural variability in the short term. In particular, phytoplankton in the Mediterranean Sea is known to respond to climate variability associated with the North Atlantic Oscillation (NAO) and El Niño Southern Oscillation (ENSO) (Basterretxea et al. 2018, Colella et al. 2016). '''KEY FINDINGS''' In the Mediterranean Sea, the average chlorophyll trend for the 1997–2024 period is slightly negative, at -0.77 ± 0.59% per year, reinforcing the findings of the previous releases. This result contrasts with the analysis by Sathyendranath et al. (2018), which reported increasing chlorophyll concentrations across all European seas. From around 2010–2011 onward, excluding the 2018–2019 period, a noticeable decline in chlorophyll levels is evident in the deseasonalized time series (green line) and in the observed maxima (grey line), particularly from 2015. This sustained decline over the past decade contributes to the overall negative trend observed in the Mediterranean Sea. '''DOI (product):''' https://doi.org/10.48670/moi-00259

'''Short description:''' For the Mediterranean Sea - The product contains daily Level-3 sea surface wind with a 1km horizontal pixel spacing using Near Real-Time Synthetic Aperture Radar (SAR) observations and their collocated European Centre for Medium-Range Weather Forecasts (ECMWF) model outputs. Products are updated several times daily to provide the best product timeliness. '''DOI (product) :''' https://doi.org/10.48670/mds-00334