Pentadal (5-year average) resolution time-series of bottom temperature for North Atlantic ocean area deeper than 1000m. Calculate the 5 year average bottom temperature at each point on the grid and then calculate the area weighted average.

The data set aims to contribute to a better biological characterization of European marine ecosystems. As such it represents probabilities of EUNIS (EUropean Nature Information System) habitat presence at Level 3 for marine habitats including information on sea ice coverage (this corresponds to EUNIS level 2 for terrestrial habitats). The map combines spatially explicit data on marine bathymetry and sea-bed with non-spatially referenced habitat information of the EUNIS classification. The objective of the data set produced by EEA and its Topic Centre ETC/ULS is to improve the biological description of marine based ecosystem types and their spatial distribution. The work supports Target 2 Action 5 of the implementation of the EU Biodiversity Strategy to 2020, established to achieve the Aichi targets of the Convention of Biological Diversity (CBD). It further addresses the MAES process (Mapping and Assessing of Ecosystems and their Services). The data set represents 2 classes of the MAES classification level 3, namely “Marine inlets and transitional waters” and “Marine”. The dataset comprises the following information: • Sea region (1 – Arctic, 2 – Atlantic, 3 – Baltic, 4 – Mediterranean, 5 – Black Sea) • Sea zone (1 – Littoral, 2 – Infralittoral, 3 – Circalittoral, 4 – Offshore circalittoral, 5 – Upper bathyal, 6 – Lower bathyal, 7 – Abyssal,8 - Coastal Lagoons, 9 - Coastal Lagoons) • Substrate (0 – undetermined substrate, 1 – rock and biogenic, 3 – coarse sediment, 4 – mixed sediment, 5 – sand, 6 – mud) • Sea ice coverage (0 – no sea ice presence, 1 – seasonal sea ice presence, 2 – perennial sea ice presence)

This product attempt to follow up on the sea level rise per stretch of coast of the North Atlantic, over 50 years as follows: • Characterization of absolute sea level trend at annual resolution, along the coasts of EU Member States (including Outermost Regions), Canada, Faroes, Greenland, Iceland, Mexico, Morocco, Norway and USA; The stretchs or coast are defined by the administrative regions of the Atlantic Coast: • from NUTS3** administrative division for EU countries (see Eurostat), and • from GADM*** administrative divisions for non-EU countries. ** Third level of Nomenclature of Territorial Units for Statistics *** Global Administrative Areas For relative sea level trend for 50 years we extract the information from coastal tide gauges data available at each stretch of coast, if there is not a tide gauge there is a data gap. The product is Provided in tabular form and as a map layer.

'''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 attempt to follow up on the sea level rise per stretch of coast of the North Atlantic, over past 10 years as follows: • Characterization of absolute sea level trend at annual resolution, along the coasts of EU Member States (including Outermost Regions), Canada, Faroes, Greenland, Iceland, Mexico, Morocco, Norway and USA; The stretchs or coast are defined by the administrative regions of the Atlantic Coast: • from NUTS3** administrative division for EU countries (see Eurostat), and • from GADM*** administrative divisions for non-EU countries. ** Third level of Nomenclature of Territorial Units for Statistics *** Global Administrative Areas For relative sea level trend for 10 years we extract the information from coastal tide gauge data available for each stretch of coast, if no tide gauge available there is a dat gap. The product is Provided in tabular form and as a map layer.

The impact of fishing on benthic habitats has previously been investigated however; a conclusive classification of potentially sensitive habitats per gear type does not exist. Currently only qualitative estimates of fishery impact using Broad-scale habitat maps are possible. Here a sensitivity matrix using both fishing pressure (fishing Intensity) and habitat sensitivity is employed to define habitat disturbance categories. The predominant fishing activities associated with physical abrasion of the seafloor area are from bottom contacting towed fishing gear. The swept area of the aforementioned gear in contact with the seabed is generally considered a function of gear width, vessel speed and fishing effort (ICES. 2015). The varying characteristics of fishing gear, their interaction with the sea floor and species being targeted; provide scope for differing interactions with subsurface (infaunal) and surface (epifaunal) dwelling communities. An evaluation of the abrasion pressure and habitat sensitivity split into surface and subsurface pressure allows greater insight to the ecological effects. Fishing intensity was calculated annually and based on the area of sea floor being swept (or swept area ratio SAR) by gear type. Calculations are based on SAR’s of gear types per area, per year. Fishing pressure ranks and habitat sensitivity ranks obtained from WGSFD working group (01 WGSFD - Report of the Working Group on Spatial Fisheries Data 2015) can be incorporated within a GIS environment to existing ICES fisheries data to provide habitat disturbance maps (fishing pressure maps+ habitat sensitivity maps) ICES. 2015. Report of the Working Group on Spatial Fisheries Data (WGSFD), 8–12 June 2015, ICES Headquarters, Copenhagen, Denmark. ICES CM 2015/SSGEPI:18. 150 pp.

Map at 1 degree resolution of 50-year linear trend in sea water temperature at 3 levels: surface, 500m, bottom.

Map the occurrence of ice at 1-degree resolution over different periods of the last century (1915-2014, 1965-2014, 2005-2014, 2009-2014). For each entire period (100, 50, 10, 5 years) find and map all cells of the 1 degree grid that experience ice conditions in at least 1 month.

'''DEFINITION''' Meridional Heat Transport is computed by integrating the heat fluxes along the zonal direction and from top to bottom of the ocean. They are given over 3 basins (Global Ocean, Atlantic Ocean and Indian+Pacific Ocean) and for all the grid points in the meridional grid of each basin. 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 PetaWatt (PW). '''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''' After an anusual 2016 year (Bricaud 2016), with a higher global meridional heat transport in the tropical band explained by, the increase of northward heat transport at 5-10 ° N in the Pacific Ocean during the El Niño event, 2017 northward heat transport is lower than the 1993-2014 reference value in the tropical band, for both Atlantic and Indian + Pacific Oceans. At the higher latitudes, 2017 northward heat transport is closed to 1993-2014 values. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00246

Pentadal time-series of the area in the North Atlantic (IHO, 1953) where ice occurred. On a 1 degree grid find all cells that experienced ice in at least 1 month of each 5 year period between 1915 and 2014, and then calculate the total area that these cells covered.