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.

This product attempt to follow up on the sea level rise per stretch of coast of the North Atlantic, over past 100 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 absolute sea level trend for 100 years we extract the information from grided sea level reconstruction datasets (using a combination of satellite and tide gauges) and extrapolate it to the nearest strecth of coast. The product is Provided in tabular form and as a map layer.

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The trend map is derived from version 5 of the global climate-quality chlorophyll time series produced by the ESA Ocean Colour Climate Change Initiative (ESA OC-CCI, Sathyendranath et al. 2019; Jackson 2020) and distributed by CMEMS. The trend detection method is based on the Census-I algorithm as described by Vantrepotte et al. (2009), where the time series is decomposed as a fixed seasonal cycle plus a linear trend component plus a residual component. The linear trend is expressed in % year -1, and its level of significance (p) calculated using a t-test. Only significant trends (p < 0.05) are included. '''CONTEXT''' Phytoplankton are key actors in the carbon cycle and, as such, recognised as an Essential Climate Variable (ECV). Chlorophyll concentration is the most widely used measure of the concentration of phytoplankton present in the ocean. Drivers for chlorophyll variability range from small-scale seasonal cycles to long-term climate oscillations and, most importantly, anthropogenic climate change. Due to such diverse factors, the detection of climate signals requires a long-term time series of consistent, well-calibrated, climate-quality data record. Furthermore, chlorophyll analysis also demands the use of robust statistical temporal decomposition techniques, in order to separate the long-term signal from the seasonal component of the time series. '''CMEMS KEY FINDINGS''' The average global trend for the 1997-2020 period was 0.59% per year, with a maximum value of 25% per year and a minimum value of -6.1% per year. Positive trends are pronounced in the high latitudes of both northern and southern hemisphehres. The significant increases in chlorophyll reported in 2016-2017 (Sathyendranath et al., 2018b) for the Atlantic and Pacific oceans at high latitudes continued to be observed after the 2020 extension, as well as the negative trends over the equatorial Pacific and the Indian Ocean Gyre. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00230

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

One product and 3 components were developed in order to fulfill the third objectif ATLANTIC_CH02_Product_5 / Distribution of ocean monitoring systems to assess climate change existing into the MPA network • Physical parameter monitoring • Chemical parameter monitoring • Biological parameter monitoring The aim of the product is the identification of ocean monitoring systems to assess climate change in MPAs.

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The sea level ocean monitoring indicator is derived from the DUACS delayed-time (DT-2021 version) altimeter gridded maps of sea level anomalies based on a stable number of altimeters (two) in the satellite constellation. The product is distributed by the Copernicus Climate Change Service and the Copernicus Marine Service (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). The regional sea level trends are derived from a linear fit of the altimeter sea level maps. The altimeter data have not been corrected for the effect of the Glacial Isostatic Adjustment nor the TOPEX-A instrumental drift during the period 1993-1998. The trend uncertainty is provided in a 90% confidence interval (Prandi et al., 2021). This estimate only considers errors related to the altimeter observation system (i.e., orbit determination errors, geophysical correction errors and inter-mission bias correction errors). The presence of the interannual signal can strongly influence the trend estimation considering to the altimeter period considered (Wang et al., 2021; Cazenave et al., 2014). The uncertainty linked to this effect is not taken into account. '''CONTEXT''' The indicator on sea level trend is a crucial index of climate change, and individual components contribute to sea level rise, including expansion due to ocean warming and melting of glaciers and ice sheets (WCRP Global Sea Level Budget Group, 2018). According to the recent IPCC 6th assessment report, global mean sea level (GMSL) increased by 0.20 (0.15 to 0.25) m over the period 1901 to 2018 with a rate 25 of rise that has accelerated since the 1960s to 3.7 (3.2 to 4.2) mm yr-1 for the period 2006–2018. Human activity was very likely the main driver of observed GMSL rise since 1970 (IPCC WGII, 2021). The weight of the different contributions evolves with time and in the recent decades the mass change has increased, contributing to the on-going acceleration of the GMSL trend (IPCC, 2022a; Legeais et al., 2020; Horwath et al., 2022). At regional scale, sea level does not change homogenously, and regional sea level change is also influenced by various other processes, with different spatial and temporal scales, such as local ocean dynamic, atmospheric forcing, Earth gravity and vertical land motion changes (IPCC WGI, 2021). Rising sea level can strongly affect population and infrastructures in coastal areas, increase their vulnerability and risks for food security, particularly in low lying areas and island states. Adverse impacts from floods, storms and tropical cyclones with related losses and damages have increased due to sea level rise, and increase their vulnerability and increase risks for food security, particularly in low lying areas and island states (IPCC, 2019, 2022b). Adaptation and mitigation measures such as the restoration of mangroves and coastal wetlands, reduce the risks from sea level rise (IPCC, 2022c). '''CMEMS KEY FINDINGS''' The altimeter mean sea level trends over the (1993/01/01, 2021/08/02) period exhibit large-scale variations at rates reaching up to more than +10 mm/yr in regions such as the western tropical Pacific Ocean. In this area, trends are mainly of thermosteric origin (Legeais et al., 2018; Meyssignac et al., 2017) in response to increased easterly winds during the last two decades associated with the decreasing Interdecadal Pacific Oscillation (IPO)/Pacific Decadal Oscillation (e.g., McGregor et al., 2012; Merrifield et al., 2012; Palanisamy et al., 2015; Rietbroek et al., 2016). Prandi et al. (2021) have estimated a regional altimeter sea level error budget from which they determine a regional error variance-covariance matrix and they provide uncertainties of the regional sea level trends. Over 1993-2019, the averaged local sea level trend uncertainty is around 0.83 mm/yr with local values ranging from 0.78 to 1.22 mm/yr. '''DOI (product):''' https://doi.org/10.48670/moi-00238

'''DEFINITION''' Estimates of Arctic sea ice extent are obtained from the surface of oceans grid cells that have at least 15% sea ice concentration. These values are cumulated in the entire Northern Hemisphere (excluding ice lakes) and from 1993 up to the year 2019 aiming to: i) obtain the Arctic sea ice extent as expressed in millions of km square (106 km2) to monitor both the large-scale variability and mean state and change. ii) to monitor the change in sea ice extent as expressed in millions of km squared per decade (106 km2/decade), or in sea ice extent loss since the beginning of the time series as expressed in percent per decade (%/decade; reference period being the first date of the key figure b) dot-dashed trend line, Vaughan et al., 2013). These trends are calculated in three ways, i.e. (i) from the annual mean values; (ii) from the March values (winter ice loss); (iii) from September values (summer ice loss). The Arctic sea ice extent used here is based on the “multi-product” (GLOBAL_MULTIYEAR_PHY_ENS_001_031) approach as introduced in the second issue of the Ocean State Report (CMEMS OSR, 2017). Five global products have been used to build the ensemble mean, and its associated ensemble spread. '''CONTEXT''' Sea ice is frozen seawater that floats on the ocean surface. This large blanket of millions of square kilometers insulates the relatively warm ocean waters from the cold polar atmosphere. The seasonal cycle of the sea ice, forming and melting with the polar seasons, impacts both human activities and biological habitat. Knowing how and how much the sea ice cover is changing is essential for monitoring the health of the Earth as sea ice is one of the highest sensitive natural environments. Variations in sea ice cover can induce changes in ocean stratification, in global and regional sea level rates and modify the key rule played by the cold poles in the Earth engine (IPCC, 2019). The sea ice cover is monitored here in terms of sea ice extent quantity. More details and full scientific evaluations can be found in the CMEMS Ocean State Report (Samuelsen et al., 2016; Samuelsen et al., 2018). '''CMEMS KEY FINDINGS''' Since the year 1993 to 2023 the Arctic sea ice extent has decreased significantly at an annual rate of -0.57*106 km2 per decade. This represents an amount of -4.8 % per decade of Arctic sea ice extent loss over the period 1993 to 2023. Over the period 1993 to 2018, summer (September) sea ice extent loss amounts to -1.18*106 km2/decade (September values), which corresponds to -14.85% per decade. Winter (March) sea ice extent loss amounts to -0.57*106 km2/decade, which corresponds to -3.42% per decade. These values slightly exceed the estimates given in the AR5 IPCC assessment report (estimate up to the year 2012) as a consequence of continuing Northern Hemisphere sea ice extent loss. Main change in the mean seasonal cycle is characterized by less and less presence of sea ice during summertime with time. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00190

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 absolute sea level trend for 10 years we extract the information from grided satellite altimetry data and extrapolate it to the nearest strecth of coast. The product is Provided in tabular form and as a map layer.

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)

'''DEFINITION''' The linear change of zonal mean subsurface temperature over the period 1993-2019 at each grid point (in depth and latitude) is evaluated to obtain a global mean depth-latitude plot of subsurface temperature trend, expressed in °C. The linear change is computed using the slope of the linear regression at each grid point scaled by the number of time steps (27 years, 1993-2019). A multi-product approach is used, meaning that the linear change is first computed for 5 different zonal mean temperature estimates. The average linear change is then computed, as well as the standard deviation between the five linear change computations. The evaluation method relies in the study of the consistency in between the 5 different estimates, which provides a qualitative estimate of the robustness of the indicator. See Mulet et al. (2018) for more details. '''CONTEXT''' Large-scale temperature variations in the upper layers are mainly related to the heat exchange with the atmosphere and surrounding oceanic regions, while the deeper ocean temperature in the main thermocline and below varies due to many dynamical forcing mechanisms (Bindoff et al., 2019). Together with ocean acidification and deoxygenation (IPCC, 2019), ocean warming can lead to dramatic changes in ecosystem assemblages, biodiversity, population extinctions, coral bleaching and infectious disease, change in behavior (including reproduction), as well as redistribution of habitat (e.g. Gattuso et al., 2015, Molinos et al., 2016, Ramirez et al., 2017). Ocean warming also intensifies tropical cyclones (Hoegh-Guldberg et al., 2018; Trenberth et al., 2018; Sun et al., 2017). '''CMEMS KEY FINDINGS''' The results show an overall ocean warming of the upper global ocean over the period 1993-2019, particularly in the upper 300m depth. In some areas, this warming signal reaches down to about 800m depth such as for example in the Southern Ocean south of 40°S. In other areas, the signal-to-noise ratio in the deeper ocean layers is less than two, i.e. the different products used for the ensemble mean show weak agreement. However, interannual-to-decadal fluctuations are superposed on the warming signal, and can interfere with the warming trend. For example, in the subpolar North Atlantic decadal variations such as the so called ‘cold event’ prevail (Dubois et al., 2018; Gourrion et al., 2018), and the cumulative trend over a quarter of a decade does not exceed twice the noise level below about 100m depth. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00244