Moving 6-year analysis of phosphate in the Arctic Ocean, for each season in the period 1965-2024. Every year of the time dimension corresponds to the 6-year centered average for each season. Winter: December-February, Spring: March-May, Summer: June-August, Autumn: September-November. Depth range (IODE standard depths): 0, 5, 10, 20, 30, 40, 50, 75, 100, 125, 150, 200, 250, 300, 400, ..., 1500, 1750, 2000, 2500m. Units: umol/l. Description of DIVA analysis: The computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.13, using GEBCO 30sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps grids is 0.125 degrees. Signal-to-noise ratio was fixed to 3.0, horizontal correlation length varying from 45 km near the coast to 150 km, and vertical correlation length varying between 25 and 1000 m. Logarithmic transformation is applied to the data prior to the analysis. Background field: analysis with signal-to-noise ratio = 20, horizontal correlation length 50-200 km, and vertical correlation length 25-1000 m.

Trends in biomass or diversity of benthic Focal Ecosystem Components across each Arctic Marine Area. STATE OF THE ARCTIC MARINE BIODIVERSITY REPORT - Chapter 4 - Page 179 - Figure 4.3

Status of marine mammal Focal Ecosystem Component stocks by Arctic Marine Area. STATE OF THE ARCTIC MARINE BIODIVERSITY REPORT - <a href="https://arcticbiodiversity.is/findings/marine-mammals" target="_blank">Chapter 3</a> - Page 157 - Figure 3.6.3

Moving 6-year analysis and visualization of Water body phosphate in the North Sea. Four seasons (December-February, March-May, June-August, September-November). Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis) tool, version 2.7.12. results were subjected to the minfield option in DIVAnd to avoid negative/underestimated values in the interpolated results; error threshold masks L1 (0.3) and L2 (0.5) are included as well as the unmasked field. The depth dimension allows visualizing the gridded field at various depths.

Moving 6-year analysis and visualization of Water body silicate in the North Sea. Four seasons (December-February, March-May, June-August, September-November). Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis) tool, version 2.7.12. results were subjected to the minfield option in DIVAnd to avoid negative/underestimated values in the interpolated results; error threshold masks L1 (0.3) and L2 (0.5) are included as well as the unmasked field. The depth dimension allows visualizing the gridded field at various depths.

Moving 6-year analysis of Water body dissolved oxygen concentration in the NorthEast Atlantic for each season: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Every year of the time dimension corresponds to the 6-year centred average of each season. 6-year periods span from 1950/1955 until 2019/2024. Observation data span from 1950 to 2024. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0, 400.0, 500.0, 600.0, 700.0, 800.0, 900.0, 1000.0, 1100.0, 1200.0, 1300.0, 1400.0, 1500.0, 1750.0, 2000.0, 2500.0, 3000.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 400km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. Signal-to-noise ratio was fixed to 1 for vertical profiles and 0.1 for time series to account for the redundancy in the time series observations. A logarithmic transformation (DIVAnd.loglin) was applied to the data prior to the analysis. Background field: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

Trends in kittiwake colonies 2001-2010, based on linear regression with year as the explanatory variable. Slope of the regression is red = negative trend, blue = positive trend; shaded circle = significant trend (at p<0.05), open circle = non-significant trend. Non-significant deviation from zero could imply a stable population, but in some cases was due to low sample size and low power. Provided with permission from Descamps et al. (in prep). STATE OF THE ARCTIC MARINE BIODIVERSITY REPORT - <a href="https://arcticbiodiversity.is/findings/seabirds" target="_blank">Chapter 3</a> - Page 135 - Figure 3.5.3 This figure is compiled from data from researchers working throughout circumpolar regions, primarily members of the Circumpolar Seabird Group, an EN of CAFF/seabirds. Dr. Sebastien Decamps conducted the analysis and produced the original figure; the full results will be available in an article in prep titled: “Descamps et al. in prep. Circumpolar dynamics of black-legged kittiwakes track large-scale environmental shifts and oceans' warming rate.” [expected submission spring 2016]. Colony population trends were analyzed using a linear regression with the year as explanatory variable. Based on slope of the regression (which cannot be exactly 0) colonies are either Declining (Slope of the regression <0) or Increasing (Slope of the regression >0). (Colonies may have had a negative but not significant slope, and could be stable but for some others, the slope is not significant due to small sample size / low power; thus we cannot say that all colonies with a non- significant slope are stable. The threshold was put at 5% to assess the significance of the trend.

Moving 6-year analysis of Water body silicate in the NorthEast Atlantic for each season: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Every year of the time dimension corresponds to the 6-year centred average of each season. 6-year periods span from 1950/1955 until 2019/2024. Observation data span from 1950 to 2024. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0, 400.0, 500.0, 600.0, 700.0, 800.0, 900.0, 1000.0, 1100.0, 1200.0, 1300.0, 1400.0, 1500.0, 1750.0, 2000.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 300km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. Signal-to-noise ratio was fixed to 1 for vertical profiles and 0.1 for time series to account for the redundancy in the time series observations. A logarithmic transformation (DIVAnd.loglin) was applied to the data prior to the analysis. Background field: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

Average September sea ice extent in 1979 (blue) compared with 2016 (white) and the median sea ice extent (yellow line) from 1981 to 2010 (Data: NSDIC 2016). STATE OF THE ARCTIC MARINE BIODIVERSITY REPORT - <a href="https://arcticbiodiversity.is/marine" target="_blank">Chapter 2</a> - Page 27 - Figure 2.4

Moving 6-year analysis of Water body phosphate in the NorthEast Atlantic for each season: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Every year of the time dimension corresponds to the 6-year centred average of each season. 6-year periods span from 1950/1955 until 2019/2024. Observation data span from 1950 to 2024. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0, 400.0, 500.0, 600.0, 700.0, 800.0, 900.0, 1000.0, 1100.0, 1200.0, 1300.0, 1400.0, 1500.0, 1750.0, 2000.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 250km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. Signal-to-noise ratio was fixed to 1 for vertical profiles and 0.1 for time series to account for the redundancy in the time series observations. A logarithmic transformation (DIVAnd.loglin) was applied to the data prior to the analysis. Background field: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.