Identified areas across the north Atlantic which have been flagged as priority locations for quality bathymetry data, in the context of expanded shipping traffic and port expansions. The reference to determine the priority survey areas in combination with shiping routes and port locations are the bathymetric data sources used for product 2( GEBCO, EMODnet bathymetry, USGS and CHS) and the depth uncertainty derived of Product 2. The adequacy assessment of the input characteristics of Product 3 is limited to the shiping routes and port locations.

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

We took inspiration from a “Matrix of marine activities” (appropriate for each IUCN management category) extracted from IUCN paper, to achieve the first objective by computing 1 product comprising the following 12 components: Product ATLANTIC_CH02_Product_1 / MPA Atlantic network classified in IUCN classification • Traditional fishing area • Sustainable fishing area (industrial) • Leisure fishing area • Leisure activity area (diving, surfing, tourist beaches) • Shipping area (shipping trajectory, aids navigation) • Scientific activity area • Renewable energy generation facility area (ocean energy facilities, wind farms) • Aquaculture area (finfish production, shellfish production) • Shipping infrastructure area (harbours, dredging area...) • Waste discharge area • Mining area (aggregate extraction, hydrocarbon extraction) • Habitation area (urban area) Each geographic information required for the components was compiled into a layer in grid format. These grids were intersected with the MPAs layer to assign each MPA a IUCN category according to the conditional matrix illustrated below : If the MPA area contains : Habitation area (urban area) The IUCN category is :V If the MPA area contains : Mining area (aggregate extraction, hydrocarbon extraction) The IUCN category is V If the MPA area contains : Waste discharge area The IUCN category is : V If the MPA area contains : Shipping infrastructure area (harbours, dredging area...) The IUCN category is IV If the MPA area contains : Aquaculture area (finfish production, shellfish production) The IUCN category is IV If the MPA area contains : Renewable energy generation facility area (ocean energy facilities, wind farms) The IUCN category is IV If the MPA area contains : Leisure fishing area The IUCN category is IV If the MPA area contains : Sustainable fishing area (industrial) The IUCN category is IV If the MPA area contains : Shipping area (shipping trajectory, aids navigation) The IUCN category is II If the MPA area contains : Leisure activity area (diving, surfing, tourist beaches) The IUCN category is Ib If the MPA area contains : Traditional fishing area The IUCN category is Ib If the MPA area contains : Scientific activity area The IUCN category is Ia

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

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

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The indicator of the Kuroshio extension phase variations is based on the standardized high frequency altimeter Eddy Kinetic Energy (EKE) averaged in the area 142-149°E and 32-37°N and computed from the DUACS (https://duacs.cls.fr) delayed-time (reprocessed version DT-2021, CMEMS SEALEVEL_GLO_PHY_L4_MY_008_047) and near real-time (CMEMS SEALEVEL_GLO_PHY_L4_NRT_OBSERVATIONS_008_046) altimeter sea level gridded products. The change in the reprocessed version (previously DT-2018) and the extension of the mean value of the EKE (now 27 years, previously 20 years) induce some slight changes not impacting the general variability of the Kuroshio extension (correlation coefficient of 0.988 for the total period, 0.994 for the delayed time period only). '''CONTEXT''' The long-term mean and trends alone do not give a complete view of the likely changes in position of unstable western boundary current extensions (Kelly et al., 2010). The Kuroshio Extension is an eastward-flowing current in the subtropical western North Pacific after the Kuroshio separates from the coast of Japan at 35°N, 140°E. Being the extension of a wind-driven western boundary current, the Kuroshio Extension is characterized by a strong variability and is rich in large-amplitude meanders and energetic eddies (Niiler et al., 2003; Qiu, 2003, 2002). The Kuroshio Extension region has the largest sea surface height variability on sub-annual and decadal time scales in the extratropical North Pacific Ocean (Jayne et al., 2009; Qiu and Chen, 2010, 2005). Prediction and monitoring of the path of the Kuroshio are of huge importance for local economies as the position of the Kuroshio extension strongly determines the regions where phytoplankton and hence fish are located. '''CMEMS KEY FINDINGS''' The different states of the Kuroshio extension phase have been presented and validated by Bessières et al. (2013) and further reported by Drévillon et al. (2018) in the Copernicus Ocean State Report #2. Two rather different states of the Kuroshio extension are observed: an ‘elongated state’ (also called ‘strong state’) corresponding to a narrow strong steady jet, and a ‘contracted state’ (also called ‘weak state’) in which the jet is weaker and more unsteady, spreading on a wider latitudinal band. When the Kuroshio Extension jet is in a contracted (elongated) state, the upstream Kuroshio Extension path tends to become more (less) variable and regional eddy kinetic energy level tends to be higher (lower). In between these two opposite phases, the Kuroshio extension jet has many intermediate states of transition and presents either progressively weakening or strengthening trends. In 2018, the indicator reveals an elongated state followed by a weakening neutral phase since then. '''DOI (product):''' https://doi.org/10.48670/moi-00222

Data from a number of different sources have been integrated to provide new perspectives on fishing activities. Vessel Monitoring Systems (VMS) record and transmit the position and speed of fishing vessels at intervals of two hours or less. Fishing time can be calculated from the VMS data and combining this parameter with vessel logbook data, maps of fishing effort and intensity at different spatial and temporal scales can be calculated. The statistical software package “R” is used to extract the required information then re-interrogated to produce maps of fishing effort or intensity per month and year. The use of Automatic Identification System (AIS) data was not considered as combining AIS data with fisheries logbook data would pose issues namely; the ability of the AIS system to be switched off, only mandatory on vessels > 15 meters in length, cost involved to purchase data, and confidentiality.

Analysis of tuna stomach contents

Annual time series of eel recruitement, (2005-2014) • Time series of glass and yellow eel for those rivers used in the annual ICES advice to the EU • Location, data availability and long term annual (LTA) eel recruiment per river mouth

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.