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

'''DEFINITION''' The regional 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 the regional products as distributed by CMEMS, derived by application of the regional chlorophyll algorithms over remote sensing reflectances (Rrs) produced by the Plymouth Marine Laboratory (PML) using the ESA Ocean Colour Climate Change Initiative processor (ESA OC-CCI, Sathyendranath et al., 2018a). '''CONTEXT''' Phytoplankton and chlorophyll concentration as their proxy respond rapidly to changes in their physical environment. In the Mediterranean Sea, these changes are seasonal and are mostly determined by light and nutrient availability (Gregg and Rousseaux, 2014). By comparing annual mean values to the climatology, we effectively remove the seasonal signal at each grid point, while retaining information on peculiar events during the year. In particular, chlorophyll anomalies in the Mediterranean Sea can then be correlated with the North Atlantic Oscillation (NAO) and El Niño Southern Oscillation (ENSO) (Basterretxea et al 2018, Colella et al 2016). '''CMEMS KEY FINDINGS''' The 2019 average chlorophyll anomaly in the Mediterranean Sea is 1.02 mg m-3 (0.005 in log10 [mg m-3]), with a maximum value of 73 mg m-3 (1.86 log10 [mg m-3]) and a minimum value of 0.04 mg m-3 (-1.42 log10 [mg m-3]). The overall east west divided pattern reported in 2016, showing negative anomalies for the Western Mediterranean Sea and positive anomalies for the Levantine Sea (Sathyendranath et al., 2018b) is modified in 2019, with a widespread positive anomaly all over the eastern basin, which reaches the western one, up to the offshore water at the west of Sardinia. Negative anomaly values occur in the coastal areas of the basin and in some sectors of the Alboràn Sea. In the northwestern Mediterranean the values switch to be positive again in contrast to the negative values registered in 2017 anomaly. The North Adriatic Sea shows a negative anomaly offshore the Po river, but with weaker value with respect to the 2017 anomaly map.

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The time series are derived from the regional chlorophyll reprocessed (REP) product as distributed by CMEMS. This dataset, derived from multi-sensor (SeaStar-SeaWiFS, AQUA-MODIS, NOAA20-VIIRS, NPP-VIIRS, Envisat-MERIS and Sentinel3A-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. (2021). 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. '''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). '''CMEMS KEY FINDINGS''' In the Mediterranean Sea, the trend average for the 1997-2020 period is slightly negative (-0.580.62% per year). Due to the change in processing techniques and chlorophyll retrieval, this trend estimate cannot be compared directly to those previously reported. The observations time series (in grey) shows minima values have been quite constant until 2015 and then there is a little decrease up to 2020, when an absolute minimum occurs with values lower than 0.04 mg m-3. Throughout the time series, maxima are variable year by year (with absolute maximum in 2015, >0.14 mg m-3), showing an evident reduction since 2016. In the last years of the series, the decrease of chlorophyll concentrations is also observed in the deseasonalized timeseries (in green) with a marked step in 2020. This attenuation of chlorophyll values in the last years results in an overall negative trend for the Mediterranean Sea. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00259

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' For the North Atlantic and Arctic oceans, the ESA Ocean Colour CCI Remote Sensing Reflectance (merged, bias-corrected Rrs) data are used to compute surface Chlorophyll (mg m-3, 1 km resolution) using the regional OC5CCI chlorophyll algorithm. The Rrs are generated by merging the data from SeaWiFS, MODIS-Aqua, MERIS, VIIRS and OLCI-3A sensors and realigning the spectra to that of the MERIS sensor. The algorithm used is OC5CCI - a variation of OC5 (Gohin et al., 2002) developed by IFREMER in collaboration with PML. As part of this development, an OC5CCI look up table was generated specifically for application over OC- CCI merged daily remote sensing reflectances. The resulting OC5CCI algorithm was tested and selected through an extensive calibration exercise that analysed the quantitative performance against in situ data for several algorithms in these specific regions. L3 products are daily files, while the L4 are monthly composites. ESA-CCI Rrs raw data are provided by PML. These are processed to produce chlorophyll concentration using the same in-house software as in the operational processing. Ocean colour technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so called ocean colour which is affected by the presence of phytoplankton. By comparing reflectances at different wavelengths and calibrating the result against in-situ measurements, an estimate of chlorophyll content can be derived. '''Processing information:''' ESA OC-CCI Rrs raw data are provided by Plymouth Marine Laboratory, currently at 4km resolution globally. These are processed to produce chlorophyll concentration using the same in-house software as in the operational processing. The entire CCI data set is consistent and processing is done in one go. Both OC CCI and the REP product are versioned. Standard masking criteria for detecting clouds or other contamination factors have been applied during the generation of the Rrs, i.e., land, cloud, sun glint, atmospheric correction failure, high total radiance, large solar zenith angle (70deg), large spacecraft zenith angle (56deg), coccolithophores, negative water leaving radiance, and normalized water leaving radiance at 560 nm 0.15 Wm-2 sr-1 (McClain et al., 1995). For the regional products, a variant of the OC-CCI chain is run to produce high resolution data at the 1km resolution necessary. A detailed description of the ESA OC-CCI processing system can be found in OC-CCI (2014e). '''Description of observation methods/instruments:''' Ocean colour technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so called ocean colour which is affected by the presence of phytoplankton. By comparing reflectances at different wavelengths and calibrating the result against in-situ measurements, an estimate of chlorophyll content can be derived. '''Quality / Accuracy / Calibration information:''' Detailed description of cal/val is given in the relevant QUID, associated validation reports and quality documentation. '''Suitability, Expected type of users / uses:''' This product is meant for use for educational purposes and for the managing of the marine safety, marine resources, marine and coastal environment and for climate and seasonal studies. '''DOI (product) :''' https://doi.org/10.48670/moi-00074

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' For the Global Ocean- In-situ observation delivered in delayed mode. This In Situ delayed mode product integrates the best available version of in situ oxygen, chlorophyll / fluorescence and nutrients data '''DOI (product) :''' https://doi.org/10.17882/86207

'''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 '''DEFINITION''' Oligotrophic subtropical gyres are regions of the ocean with low levels of nutrients required for phytoplankton growth and low levels of surface chlorophyll-a whose concentration can be quantified through satellite observations. The gyre boundary has been defined using a threshold value of 0.15 mg m-3 chlorophyll for the Atlantic gyres (Aiken et al. 2016), and 0.07 mg m-3 for the Pacific gyres (Polovina et al. 2008). The area inside the gyres for each month is computed using monthly chlorophyll data from which the monthly climatology is subtracted to compute anomalies. A gap filling algorithm has been utilized to account for missing data. Trends in the area anomaly are then calculated for the entire study period (September 1997 to December 2020). '''CONTEXT''' Oligotrophic gyres of the oceans have been referred to as ocean deserts (Polovina et al. 2008). They are vast, covering approximately 50% of the Earth’s surface (Aiken et al. 2016). Despite low productivity, these regions contribute significantly to global productivity due to their immense size (McClain et al. 2004). Even modest changes in their size can have large impacts on a variety of global biogeochemical cycles and on trends in chlorophyll (Signorini et al. 2015). Based on satellite data, Polovina et al. (2008) showed that the areas of subtropical gyres were expanding. The Ocean State Report (Sathyendranath et al. 2018) showed that the trends had reversed in the Pacific for the time segment from January 2007 to December 2016. '''CMEMS KEY FINDINGS''' The trend in the North Atlantic gyre area for the 1997 Sept – 2020 December period was positive, with a 0.39% year-1 increase in area relative to 2000-01-01 values. This trend has decreased compared with the 1997-2019 trend of 0.45%, and is statistically significant (p<0.05). During the 1997 Sept – 2020 December period, the trend in chlorophyll concentration was positive (0.24% year-1) inside the North Atlantic gyre relative to 2000-01-01 values. This time series extension has resulted in a reversal in the rate of change, compared with the -0.18% trend for the 1997-209 period and is statistically significant (p<0.05). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00226

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The BALTIC_OMI_TEMPSAL_sst_area_averaged_anomalies product includes time series of monthly mean SST anomalies over the period 1993-2021, relative to the 1993-2014 climatology, averaged for the Baltic Sea. The OMI time series runs from Jan 1, 1993 to December 31, 2021 and is constructed by calculating monthly averages from the daily level 4 SST analysis fields of the SST_BAL_SST_L4_REP_OBSERVATIONS_010_016 product from 1993 to 2021. See the Copernicus Marine Service Ocean State Reports (section 1.1 in Von Schuckmann et al., 2016; section 3 in Von Schuckmann et al., 2018) for more information on the OMI product. '''CONTEXT''' Sea Surface Temperature (SST) is an Essential Climate Variable (GCOS), that is an important input for initialis-ing numerical weather prediction models and fundamental for understanding air-sea interactions and moni-toring climate change (GCOS 2010). The Baltic Sea is a region that requires special attention regarding the use of satellite SST records and the assessment of climatic variability (Høyer and She 2007; Høyer and Karagali 2016). The Baltic Sea is a semi-enclosed basin affected bynatural variability, influenced by large-scale atmos-pheric processes and by the vicinity of land. In addition, the Baltic Sea is one of the largest brackish seas in the world. When analysing regional-scale climate variability, all these effects have to be considered, which re-quires dedicated regional and validated SST products. Satellite observations have previously been used to ana-lyse the climatic SST signals in the North Sea and Baltic Sea (BACC II Author Team 2015; Lehmann et al. 2011). Recently, Høyer and Karagali (2016) demonstrated that the Baltic Sea had warmed 1-2 oC from 1982 to 2012 considering all months of the year and 3-5 oC when only July-September months were considered. This was corroborated in the Ocean State Reports (section 1.1 in Von Schuckmann et al., 2016 and section 3 in Von Schuckmann et al., 2018). '''CMEMS KEY FINDINGS''' The basin-average trend of SST anomalies for Baltic Sea region amounts to 0.049±0.006 °C/year over the pe-riod 1993-2021 which corresponds to an average warming of 1.42°C. Adding the North Sea area, the average trend amounts to 0.03±0.003 °C/year over the same period, which corresponds to an average warming of 0.87°C for the entire region since 1993. '''DOI (product):''' https://doi.org/10.48670/moi-00205

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The BALTIC_OMI_TEMPSAL_sst_trend product includes the cumulative/net trend in sea surface temperature anomalies for the Baltic Sea from 1993-2021. The cumulative trend is the rate of change (°C/year) scaled by the number of years (29 years). The SST Level 4 analysis products that provide the input to the trend calculations are taken from the reprocessed product SST_BAL_SST_L4_REP_OBSERVATIONS_010_016 with a recent update to include 2021. The product has a spatial resolution of 0.02 degrees in latitude and longitude. The OMI time series runs from Jan 1, 1993 to December 31, 2021 and is constructed by calculating monthly averages from the daily level 4 SST analysis fields of the SST_BAL_SST_L4_REP_OBSERVATIONS_010_016 from 1993 to 2021. See the Copernicus Marine Service Ocean State Reports for more information on the OMI product (section 1.1 in Von Schuckmann et al., 2016; section 3 in Von Schuckmann et al., 2018). The times series of monthly anomalies have been used to calculate the trend in SST using Sen’s method with confidence intervals from the Mann-Kendall test (section 3 in Von Schuckmann et al., 2018). '''CONTEXT''' SST is an essential climate variable that is an important input for initialising numerical weather prediction models and fundamental for understanding air-sea interactions and monitoring climate change. The Baltic Sea is a region that requires special attention regarding the use of satellite SST records and the assessment of climatic variability (Høyer and She 2007; Høyer and Karagali 2016). The Baltic Sea is a semi-enclosed basin with natural variability and it is influenced by large-scale atmospheric processes and by the vicinity of land. In addition, the Baltic Sea is one of the largest brackish seas in the world. When analysing regional-scale climate variability, all these effects have to be considered, which requires dedicated regional and validated SST products. Satellite observations have previously been used to analyse the climatic SST signals in the North Sea and Baltic Sea (BACC II Author Team 2015; Lehmann et al. 2011). Recently, Høyer and Karagali (2016) demonstrated that the Baltic Sea had warmed 1-2oC from 1982 to 2012 considering all months of the year and 3-5oC when only July- September months were considered. This was corroborated in the Ocean State Reports (section 1.1 in Von Schuckmann et al., 2016; section 3 in Von Schuckmann et al., 2018). '''CMEMS KEY FINDINGS''' SST trends were calculated for the Baltic Sea area and the whole region including the North Sea, over the period January 1993 to December 2021. The average trend for the Baltic Sea domain (east of 9°E longitude) is 0.049 °C/year, which represents an average warming of 1.42 °C for the 1993-2021 period considered here. When the North Sea domain is included, the trend decreases to 0.03°C/year corresponding to an average warming of 0.87°C for the 1993-2021 period. Trends are highest for the Baltic Sea region and North Atlantic, especially offshore from Norway, compared to other regions. '''DOI (product):''' https://doi.org/10.48670/moi-00206

'''Short description: ''' For the Global Ocean - In-situ observation yearly delivery in delayed mode of Ocean surface currents. '''Detailed description: ''' The In Situ delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for Ocean surface currents. The data are collected from the Surface Drifter Data Assembly Centre (SD-DAC at NOAA AOML) completed by European data provided by EUROGOOS regional systems and national systems by the regional INS TAC components. All surface drifters data have been processed to check for drogue loss. Drogued and undrogued drifting buoy surface ocean currents are provided with a drogue presence flag as well as a wind slippage correction for undrogued buoy. '''Processing information: ''' From the near real time INS TAC product validated on a daily and weekly basis for forecasting purposes, and from the SD-DAC quality controlled dataset a scientifically validated product is created . It s a """"reference product"""" updated on a yearly basis. This product has been processed using a method that checks for drogue loss. Altimeter and wind data have been used to extract the direct wind slippage from the total drifting buoy velocities. The obtained wind slippage values have then been analyzed to identify probable undrogued data among the drifting buoy velocities dataset. A simple procedure has then been applied to produce an updated dataset including a drogue presence flag as well as a wind slippage correction. '''Suitability, Expected type of users / uses: ''' The product is designed to be assimilated into or for validation purposes of operational models operated by ocean forecasting centers for reanalysis purposes or for research community. These users need data aggregated and quality controlled in a reliable and documented manner.