Regional assessment of status and trends in PAH concentrations in sediment

  • Overview
  • Summary of individual time series results
  • Regional assessment - methods
  • Regional assessment - availability of time series
  • Regional trends
  • Regional status compared to the Effects Range Low (ERL)
  • Regional status compared to the Background Assessment Concentration (BAC)
  • PAH profiles

Overview

The 2023 OSPAR CEMP assessment describes the status and trends of contaminant concentrations in biota and sediment at monitoring stations in the OSPAR area. Assessments are made for a large number of time series, each of a single contaminant in a single species and tissue (for biota) at a single monitoring station. This document is one of a series that synthesises the results of the individual time series to assess status and trends at the MIME regional level. In particular, it considers PAH concentrations in sediment, where the PAHs are NAP, ACNLE, ACNE, FLE, PA, ANT, DBT, FLU, PYR, BAA, CHR, BAP, BEP, BBKF, DBAHA, PER, BGHIP, ICDP. For simplicity, the term ‘region’ is used throughout to describe MIME regions. OSPAR regions are always referred to as such.

The regional status assessments presented in this document are integrated with regional status assessments for other contaminants and biological effects using the CHASE approach. Details can be found here.

The following sections:

  • summarise the individual time series results
  • describe the methods used to combine the individual time series results to estimate regional status and trends
  • present the estimates of regional trends
  • present the estimates of regional status


Summary of individual time series results

A time series of PAH concentrations is assessed for status if:

  • there is at least one year with data in the period 2016 to 2021
  • there are at least three years of data over the whole time series
  • a parametric model can be fitted to the data and used to estimate the mean concentration in the final monitoring year (or, occasionally, if a non-parametric test of status is applied)

The conditions are more stringent for trends. Specifically, a time series is assessed for trends if:

  • there is at least one year with data in the period 2016 to 2021
  • there are at least five years of data over the whole time series
  • a parametric model can be fitted to the data and used to estimate the trend in mean concentration

Note that all trend assessments for individual time series and most status assessments are based on the fit of a parametric model. This is important because only the parametric results are passed into the regional assessments in the following sections. 9 of the 3132 PAH time series were assessed for status using a non-parametric test.

The first tab below shows all the monitoring stations where there are individual time series assessments of trend or status for PAHs. The purple circles are stations where there is a trend assessment for at least one PAH, and the light blue cirles are stations where there are only status assessments.

The next set of tabs show the trend and status assessments for each PAH in turn. The colours have the following meaning:

  • blue: the mean concentration is significantly (p < 0.05) below the Background Assessment Concentration (BAC)
  • green: the mean concentration is significantly (p < 0.05) below the Effects Range Low (ERL)
  • orange: there is no ERL and the mean concentration is not significantly below the BAC
  • red: the mean concentration is not signficantly below the ERL
  • black: there are no assessment criteria available

And the shapes have the following meaning:

  • downward triangle: the mean concentration is significantly decreasing (p < 0.05)
  • circle: there is no significant change in mean concentration (p > 0.05) or there are too few years to test for trends
  • upward triangle: the mean concentration is significantly increasing (p < 0.05)

The final set of tabs give:

  • the number of time series with upwards, downwards or no trend in each MIME and OSPAR region
  • the same information presented as proportions
  • the number of time series with blue, green, orange or red status in each MIME and OSPAR region
  • the same information presented as proportions

Note that only eleven of the PAHs have assessment criteria: NAP, PA, ANT, DBT, FLU, PYR, BAA, CHR, BAP, BGHIP, ICDP. Of these, DBT only has an ERL and BGHIP and ICDP only have a BAC. (NAP has an ERL and a BAC in all regions apart from the Iberian Sea and Gulf of Cadiz, where only the ERL applies.)


Stations


NAP


ACNLE


ACNE


FLE


PA


ANT


DBT


FLU


PYR


BAA


CHR


BAP


BEP


BBKF


DBAHA


PER


BGHIP


ICDP


Trend summary (numbers)

Number of time series with upwards, downwards or no trend by MIME and OSPAR region

OSPAR region MIME region status NAP ACNLE ACNE FLE PA ANT DBT FLU PYR BAA CHR BAP BEP BBKF DBAHA PER BGHIP ICDP total
2 Northern North Sea upward trend 2 0 0 0 0 1 0 1 1 1 1 1 1 0 0 1 1 1 12
no trend 10 7 12 12 13 12 11 14 15 14 13 13 12 13 15 13 12 14 225
downward trend 4 2 1 4 4 2 3 2 1 2 3 3 4 4 2 3 4 2 50
Southern North Sea upward trend 4 1 1 2 5 8 0 3 2 1 0 0 0 1 1 0 0 2 31
no trend 29 15 20 30 44 43 0 43 42 46 46 46 28 37 46 13 39 36 603
downward trend 3 1 5 4 8 3 0 12 13 9 10 10 6 8 6 2 15 18 133
Channel upward trend 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
no trend 2 3 3 2 2 2 0 10 2 3 3 6 2 6 3 1 5 4 59
downward trend 1 0 0 1 1 1 0 0 1 0 0 0 1 2 0 2 0 0 10
total upward trend 6 1 1 2 5 9 0 4 3 2 1 1 1 1 1 1 1 3 43
no trend 41 25 35 44 59 57 11 67 59 63 62 65 42 56 64 27 56 54 887
downward trend 8 3 6 9 13 6 3 14 15 11 13 13 11 14 8 7 19 20 193
3 Irish and Scottish West Coast upward trend 0 0 0 0 1 0 0 1 2 0 1 1 1 1 1 0 1 1 11
no trend 6 4 5 5 6 7 5 5 5 6 5 6 6 6 6 6 6 6 101
downward trend 1 0 0 1 0 0 0 1 0 1 1 0 0 0 0 1 0 0 6
Irish Sea upward trend 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 3
no trend 10 6 12 11 14 12 6 12 12 13 13 13 11 10 13 7 13 12 200
downward trend 4 5 1 1 1 3 1 3 3 2 2 2 3 4 2 6 2 2 47
Celtic Sea upward trend 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
no trend 1 1 1 1 1 1 0 1 1 1 0 1 1 1 1 1 1 1 16
downward trend 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1
total upward trend 1 0 0 0 1 0 0 1 2 0 1 1 2 1 1 0 1 2 14
no trend 17 11 18 17 21 20 11 18 18 20 18 20 18 17 20 14 20 19 317
downward trend 5 5 1 2 1 3 1 4 3 3 4 2 3 4 2 7 2 2 54
4 Gulf of Cadiz upward trend 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1
no trend 0 0 0 0 12 11 0 9 9 12 11 12 10 13 7 0 12 7 125
downward trend 0 0 0 0 1 1 0 4 4 1 2 1 3 0 0 0 1 6 24
total upward trend 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1
no trend 0 0 0 0 12 11 0 9 9 12 11 12 10 13 7 0 12 7 125
downward trend 0 0 0 0 1 1 0 4 4 1 2 1 3 0 0 0 1 6 24


Trend summary (proportions)

Proportion of time series with upwards, downwards or no trend by MIME and OSPAR region

OSPAR region MIME region status NAP ACNLE ACNE FLE PA ANT DBT FLU PYR BAA CHR BAP BEP BBKF DBAHA PER BGHIP ICDP total
2 Northern North Sea upward trend 13 0 0 0 0 7 0 6 6 6 6 6 6 0 0 6 6 6 4
no trend 62 78 92 75 76 80 79 82 88 82 76 76 71 76 88 76 71 82 78
downward trend 25 22 8 25 24 13 21 12 6 12 18 18 23 24 12 18 23 12 18
Southern North Sea upward trend 11 6 4 6 9 15 5 3 2 0 0 0 2 2 0 0 4 4
no trend 81 88 77 83 77 80 74 74 82 82 82 82 81 87 87 72 64 79
downward trend 8 6 19 11 14 5 21 23 16 18 18 18 17 11 13 28 32 17
Channel upward trend 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
no trend 67 100 100 67 67 67 100 67 100 100 100 67 75 100 33 100 100 86
downward trend 33 0 0 33 33 33 0 33 0 0 0 33 25 0 67 0 0 14
total upward trend 11 4 3 4 6 13 0 5 4 3 1 1 2 1 1 3 1 4 4
no trend 74 86 83 80 77 79 79 79 77 83 82 82 78 79 88 77 74 70 79
downward trend 15 10 14 16 17 8 21 16 19 14 17 17 20 20 11 20 25 26 17
3 Irish and Scottish West Coast upward trend 0 0 0 0 14 0 0 14 29 0 14 14 14 14 14 0 14 14 9
no trend 86 100 100 83 86 100 100 72 71 86 72 86 86 86 86 86 86 86 86
downward trend 14 0 0 17 0 0 0 14 0 14 14 0 0 0 0 14 0 0 5
Irish Sea upward trend 7 0 0 0 0 0 0 0 0 0 0 0 7 0 0 0 0 7 1
no trend 66 55 92 92 93 80 86 80 80 87 87 87 73 71 87 54 87 80 80
downward trend 27 45 8 8 7 20 14 20 20 13 13 13 20 29 13 46 13 13 19
Celtic Sea upward trend 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
no trend 100 100 100 100 100 100 100 100 100 0 100 100 100 100 100 100 100 94
downward trend 0 0 0 0 0 0 0 0 0 100 0 0 0 0 0 0 0 6
total upward trend 4 0 0 0 5 0 0 4 9 0 4 4 9 5 4 0 4 9 4
no trend 74 69 95 89 91 87 92 78 78 87 78 87 78 77 87 67 87 82 82
downward trend 22 31 5 11 4 13 8 18 13 13 18 9 13 18 9 33 9 9 14
4 Gulf of Cadiz upward trend 0 8 0 0 0 0 0 0 0 0 0 0 1
no trend 92 84 69 69 92 85 92 77 100 100 92 54 83
downward trend 8 8 31 31 8 15 8 23 0 0 8 46 16
total upward trend 0 8 0 0 0 0 0 0 0 0 0 0 1
no trend 92 84 69 69 92 85 92 77 100 100 92 54 83
downward trend 8 8 31 31 8 15 8 23 0 0 8 46 16


Status summary (numbers)

Number of time series with each status by MIME and OSPAR region

OSPAR region MIME region status NAP PA ANT DBT FLU PYR BAA CHR BAP BGHIP ICDP total
2 Northern North Sea blue 13 12 12 0 11 10 10 10 11 12 9 110
green 4 5 7 21 10 11 11 11 10 0 0 90
orange 0 0 0 0 0 0 0 0 0 9 12 21
red 4 4 2 0 0 0 0 0 0 0 0 10
Southern North Sea blue 0 4 1 0 3 2 4 6 13 43 46 122
green 36 56 62 5 61 63 60 59 51 0 0 453
orange 0 0 0 0 0 0 0 0 0 21 19 40
red 9 5 2 0 1 0 1 0 1 0 0 19
Channel blue 2 4 0 0 4 2 2 2 5 12 21 54
green 20 30 25 4 37 40 34 40 38 0 0 268
orange 0 0 0 0 0 0 0 0 0 36 27 63
red 10 15 18 0 8 7 13 7 6 0 0 84
total blue 15 20 13 0 18 14 16 18 29 67 76 286
green 60 91 94 30 108 114 105 110 99 0 0 811
orange 0 0 0 0 0 0 0 0 0 66 58 124
red 23 24 22 0 9 7 14 7 7 0 0 113
3 Irish and Scottish West Coast blue 1 4 5 0 6 5 4 3 6 8 8 50
green 8 5 4 6 3 4 4 6 3 0 0 43
orange 0 0 0 0 0 0 0 0 0 1 1 2
red 0 0 0 0 0 0 1 0 0 0 0 1
Irish Sea blue 0 0 0 0 0 0 0 0 0 2 4 6
green 16 13 15 10 16 15 14 14 15 0 0 128
orange 0 0 0 0 0 0 0 0 0 15 13 28
red 1 4 2 0 1 2 3 3 2 0 0 18
Celtic Sea blue 0 1 0 0 1 1 1 1 2 2 3 12
green 3 4 3 1 4 4 4 4 3 0 0 30
orange 0 0 0 0 0 0 0 0 0 3 2 5
red 1 0 1 0 0 0 0 0 0 0 0 2
total blue 1 5 5 0 7 6 5 4 8 12 15 68
green 27 22 22 17 23 23 22 24 21 0 0 201
orange 0 0 0 0 0 0 0 0 0 19 16 35
red 2 4 3 0 1 2 4 3 2 0 0 21
4 Iberian Sea blue 0 4 3 0 6 7 4 9 11 9 8 61
green 0 21 19 0 22 23 20 19 19 0 0 143
orange 0 0 0 0 0 0 0 0 0 25 25 50
red 0 9 10 0 6 4 9 6 4 0 0 48
Gulf of Cadiz blue 0 18 2 0 18 16 13 18 16 14 18 133
green 0 6 21 0 6 8 11 6 7 0 0 65
orange 0 0 0 0 0 0 0 0 0 12 7 19
red 0 1 1 0 1 1 1 1 1 0 0 7
total blue 0 22 5 0 24 23 17 27 27 23 26 194
green 0 27 40 0 28 31 31 25 26 0 0 208
orange 0 0 0 0 0 0 0 0 0 37 32 69
red 0 10 11 0 7 5 10 7 5 0 0 55


Status summary (proportions)

Proportion of time series with each status by MIME and OSPAR region

OSPAR region MIME region status NAP PA ANT DBT FLU PYR BAA CHR BAP BGHIP ICDP total
2 Northern North Sea blue 62 57 57 0 52 48 48 48 52 57 43 48
green 19 24 33 100 48 52 52 52 48 0 0 39
orange 0 0 0 0 0 0 0 0 0 43 57 9
red 19 19 10 0 0 0 0 0 0 0 0 4
Southern North Sea blue 0 6 2 0 5 3 6 9 20 67 71 19
green 80 86 95 100 94 97 92 91 78 0 0 72
orange 0 0 0 0 0 0 0 0 0 33 29 6
red 20 8 3 0 1 0 2 0 2 0 0 3
Channel blue 6 8 0 0 8 4 4 4 10 25 44 12
green 63 61 58 100 76 82 69 82 78 0 0 57
orange 0 0 0 0 0 0 0 0 0 75 56 13
red 31 31 42 0 16 14 27 14 12 0 0 18
total blue 15 15 10 0 13 10 12 13 22 50 57 21
green 61 67 73 100 80 85 78 82 73 0 0 61
orange 0 0 0 0 0 0 0 0 0 50 43 9
red 24 18 17 0 7 5 10 5 5 0 0 9
3 Irish and Scottish West Coast blue 11 44 56 0 67 56 45 33 67 89 89 52
green 89 56 44 100 33 44 44 67 33 0 0 45
orange 0 0 0 0 0 0 0 0 0 11 11 2
red 0 0 0 0 0 0 11 0 0 0 0 1
Irish Sea blue 0 0 0 0 0 0 0 0 0 12 24 3
green 94 76 88 100 94 88 82 82 88 0 0 71
orange 0 0 0 0 0 0 0 0 0 88 76 16
red 6 24 12 0 6 12 18 18 12 0 0 10
Celtic Sea blue 0 20 0 0 20 20 20 20 40 40 60 25
green 75 80 75 100 80 80 80 80 60 0 0 61
orange 0 0 0 0 0 0 0 0 0 60 40 10
red 25 0 25 0 0 0 0 0 0 0 0 4
total blue 3 16 17 0 23 19 16 13 26 39 48 21
green 90 71 73 100 74 74 71 77 68 0 0 62
orange 0 0 0 0 0 0 0 0 0 61 52 11
red 7 13 10 0 3 7 13 10 6 0 0 6
4 Iberian Sea blue 12 9 17 20 12 26 32 26 24 20
green 62 60 65 68 61 56 56 0 0 47
orange 0 0 0 0 0 0 0 74 76 17
red 26 31 18 12 27 18 12 0 0 16
Gulf of Cadiz blue 72 8 72 64 52 72 67 54 72 59
green 24 88 24 32 44 24 29 0 0 29
orange 0 0 0 0 0 0 0 46 28 9
red 4 4 4 4 4 4 4 0 0 3
total blue 37 9 41 39 29 46 46 38 45 37
green 46 71 47 53 54 42 45 0 0 40
orange 0 0 0 0 0 0 0 62 55 13
red 17 20 12 8 17 12 9 0 0 10


Regional assessment - methods

Tabulating the number of time series with each status category by region provides a quick summary of the individual time series results. However, it does not provide an objective regional assessment of status. Similarly, tabulating the number of time series with an upward or downward trend does not provide an objective regional assessment of trend. This section describes how the individual time series results can by synthesised in a meta-analysis to assess both status and trend at the regional level.

Trends

For a regional trend assessment, the trend in each time series is summarised by the estimated change in log concentration over the last twenty years (or shorter if the time series doesn’t extend that far back). Regional trends are then estimated by fitting the following linear mixed model by restricted maximum likelihood:

  • response: trend (yearly change in log concentration)
  • fixed model: region
  • random model: PAH + region.PAH + station + trend estimation variation + residual variation

The fixed model means that a trend is estimated for each region. Trends in the different PAHs are assumed to be broadly similar, because inputs and control measures are likely to affect all PAHs in a similar way. Each regional trend can thus be interpreted as the common trend in PAH concentrations in that region (i.e. the trend averaged across all PAHs).

The random model has five terms:

  • PAH allows for variation in trend between PAHs common to all regions
  • region.PAH allows for additional variation in trend between PAHs within regions
  • station allows for variation in trend between stations common to all PAHs
  • trend estimation variation is the variance of the trend estimates from the individual time series analysis, assumed known and fixed
  • residual variation is the variation that cannot be explained by any of the fixed effects or the other random effects

The meta-analysis is restricted to time series from monitoring stations that are classified as ‘representative’ or ‘baseline’ (near pristine conditions or only subject to very diffuse inputs). ‘Impacted’ stations (those close to a point source) are omitted because trends there will likely reflect changes due to the point source and will be ‘unrepresentative’ at the regional level.

The meta-analysis is further restricted to:

  • regions with at least three trend stations with good geographic spread. Three stations is considered the minimum required to provide some sort of evidence base at the regional level
  • PAHs with trend information for at least three stations (over all regions): PAHs are treated as random effects in the meta-anlaysis, so information about trends in one PAH can be ‘borrowed’ from the trend information in the other PAHs

Status

For a regional status assessment, the status of each time series is summarised by the difference between the estimated mean log concentration in the final monitoring year and the log assessment concentration. This ensures that status is always measured on the same scale, even though the assessment concentration might vary between PAHs and regions. Essentially the same linear mixed model as for trends is then fitted:

  • response: status (mean log concentration - log assessment concentration)
  • fixed model: region
  • random model: PAH + region.PAH + station + status estimation variation + residual variation

where status estimation variation is the variation in the status estimates from the individual time series analysis, assumed known and fixed.

There are no restrictions on the time series used in the status meta-analysis based on the classification of the monitoring station; time series from baseline, representative and impacted stations are all included. However, the few time series with a non-parametric assessment of status must be excluded, because there is no summary measure of status to use in the mixed model.

Again, the meta-analysis is restricted to regions with at least three status stations with good geographic spread and to PAHs with status information for at least three stations (over all regions).

Presentation

The results of both the trend and status regional assessments are back-transformed for presentation. The estimated regional trend is then interpreted as the percentage yearly change in concentration and the estimated regional status as the ratio of the mean regional concentration to the assessment concentration. To illustrate the latter, a value of 1 indicates that the mean regional concentration is equal to the assessment concentration; a value of 0.5 indicates that it is half the assessment concentration, and a value of 2 indicates that it is twice the assessment concentration.


Regional assessment - availability of time series

This section provides more detail on the number and geographic spread of the time series available for the meta-analyses. The tabs show the stations with parametric trend and status assessments for each OSPAR region in turn, having excluded trend assessments at impacted stations. The purple circles are stations where there is a trend assessment for at least one PAH, and the light blue cirles are stations where there are only status assessments.

The following regions have an insufficient number of time series or spread of stations and are excluded from the regional assessment:

  • trend: Celtic Sea
  • status:


OSPAR region 1


OSPAR region 2


OSPAR region 3


OSPAR region 4


Regional trends

The first 4 tabs show:

  • Regional trends: estimates of the regional trend (averaged over PAHs) with pointwise 95% confidence intervals
  • Regional trends by PAH: regional trend estimates for each PAH with pointwise 95% confidence intervals plotted by PAH
  • Regional trends by region: regional trend estimates for each PAH with pointwise 95% confidence intervals plotted by region
  • Individual time series: the trend estimates from the individual time series on which the regional assessment is based. The estimates for regions with insufficient geographic spread (and excluded from the meta-analysis) are also shown.

The symbols in all these plots have the following interpretation:

  • downward triangle: the mean concentration is significantly decreasing (p < 0.05)
  • circle: there is no significant change in mean concentration (p > 0.05)
  • upward triangle: the mean concentration is significantly increasing (p < 0.05)

The final 2 tabs show the estimates of the regional trend (averaged over PAHs) and the regional trend by PAH


Regional trends


Regional trends by PAH


Regional trends by region


Individual time series


Estimates by region

This table shows the estimates of the regional trend:

  • trend: the estimated mean yearly change in log concentration across the region (multiplied by 100 for presentation)
  • se: the corresponding standard error
  • lower, upper: the corresponding pointwise 95% confidence limits
  • % yearly change: the estimated mean % yearly change in concentration across the region
  • %yc lower, %yc upper: the corresponding pointwise 95% confidence limits


MIME region trend se lower upper % yearly change %yc lower %yc upper
Northern North Sea -0.96 0.86 -2.64 0.73 -0.95 -2.61 0.73
Southern North Sea -0.67 0.51 -1.67 0.33 -0.67 -1.65 0.33
Channel -3.11 1.95 -6.93 0.71 -3.06 -6.69 0.71
Irish and Scottish West Coast 2.27 1.31 -0.30 4.84 2.29 -0.30 4.96
Irish Sea -0.72 0.97 -2.62 1.19 -0.71 -2.59 1.20
Gulf of Cadiz -1.98 0.98 -3.91 -0.05 -1.96 -3.83 -0.05


Estimates by region and PAH

This table shows the estimates of the regional trend by determinand:

  • trend: the estimated mean yearly change in log concentration across the region (multiplied by 100 for presentation)
  • se: the corresponding standard error
  • lower, upper: the corresponding pointwise 95% confidence limits
  • % yearly change: the estimated mean % yearly change in concentration across the region
  • %yc lower, %yc upper: the corresponding pointwise 95% confidence limits


MIME region PAH trend se lower upper % yearly change %yc lower %yc upper
Northern North Sea NAP -1.88 0.98 -3.80 0.03 -1.86 -3.73 0.03
ACNLE -1.51 1.21 -3.87 0.85 -1.50 -3.80 0.85
ACNE -0.68 1.00 -2.64 1.27 -0.68 -2.60 1.27
FLE -1.38 0.93 -3.21 0.45 -1.37 -3.16 0.45
PA -0.95 0.91 -2.74 0.84 -0.95 -2.71 0.84
ANT -1.05 0.97 -2.95 0.84 -1.05 -2.91 0.85
DBT -1.05 1.00 -3.01 0.92 -1.04 -2.97 0.92
FLU -0.72 0.92 -2.51 1.08 -0.71 -2.48 1.09
PYR -0.66 0.92 -2.46 1.14 -0.66 -2.43 1.15
BAA -0.79 0.93 -2.61 1.04 -0.78 -2.58 1.04
CHR -1.56 0.94 -3.40 0.28 -1.55 -3.34 0.28
BAP -0.80 0.92 -2.61 1.01 -0.80 -2.58 1.01
BEP -0.63 0.92 -2.44 1.17 -0.63 -2.41 1.18
BBKF -1.68 0.94 -3.52 0.15 -1.67 -3.46 0.15
DBAHA -0.96 0.98 -2.88 0.96 -0.95 -2.84 0.96
PER -1.00 0.98 -2.92 0.91 -1.00 -2.88 0.91
BGHIP -0.29 0.93 -2.11 1.53 -0.29 -2.09 1.54
ICDP 0.36 0.94 -1.49 2.20 0.36 -1.48 2.22
Southern North Sea NAP 0.51 0.53 -0.52 1.55 0.52 -0.52 1.56
ACNLE 0.86 0.80 -0.71 2.43 0.86 -0.71 2.46
ACNE -1.42 0.60 -2.59 -0.25 -1.41 -2.56 -0.25
FLE -0.08 0.54 -1.14 0.99 -0.08 -1.14 0.99
PA -0.14 0.47 -1.08 0.79 -0.14 -1.07 0.79
ANT 0.64 0.48 -0.30 1.59 0.65 -0.30 1.60
FLU -0.94 0.47 -1.86 -0.02 -0.94 -1.84 -0.02
PYR -1.06 0.47 -1.97 -0.14 -1.05 -1.95 -0.14
BAA -0.92 0.47 -1.84 0.00 -0.91 -1.83 0.00
CHR -0.81 0.47 -1.74 0.11 -0.81 -1.72 0.11
BAP -1.02 0.48 -1.96 -0.09 -1.02 -1.94 -0.09
BEP -0.66 0.48 -1.61 0.29 -0.66 -1.60 0.29
BBKF -1.05 0.48 -1.98 -0.12 -1.05 -1.97 -0.12
DBAHA -0.60 0.50 -1.57 0.37 -0.60 -1.56 0.38
PER -1.39 0.71 -2.78 -0.00 -1.38 -2.74 -0.00
BGHIP -1.50 0.47 -2.42 -0.59 -1.49 -2.39 -0.59
ICDP -1.81 0.47 -2.74 -0.89 -1.80 -2.70 -0.89
Channel NAP -3.28 2.11 -7.41 0.86 -3.22 -7.14 0.86
ACNLE -2.25 2.16 -6.49 1.99 -2.22 -6.28 2.01
ACNE -2.52 2.16 -6.76 1.71 -2.49 -6.53 1.72
FLE -3.01 2.09 -7.12 1.09 -2.97 -6.87 1.10
PA -3.14 2.07 -7.20 0.93 -3.09 -6.95 0.93
ANT -3.06 2.10 -7.17 1.06 -3.01 -6.92 1.07
FLU -3.03 2.12 -7.18 1.13 -2.98 -6.93 1.13
PYR -3.58 2.09 -7.68 0.53 -3.51 -7.39 0.53
BAA -2.97 2.16 -7.20 1.26 -2.93 -6.95 1.27
CHR -3.27 2.16 -7.50 0.97 -3.21 -7.23 0.98
BAP -3.36 2.10 -7.47 0.76 -3.30 -7.20 0.76
BEP -3.53 2.08 -7.61 0.55 -3.47 -7.33 0.55
BBKF -4.19 2.08 -8.26 -0.11 -4.10 -7.93 -0.11
DBAHA -2.77 2.11 -6.91 1.37 -2.73 -6.67 1.37
PER -4.06 2.06 -8.10 -0.02 -3.98 -7.78 -0.02
BGHIP -2.52 2.08 -6.60 1.55 -2.49 -6.39 1.56
ICDP -2.32 2.10 -6.43 1.79 -2.29 -6.23 1.81
Irish and Scottish West Coast NAP 1.72 1.46 -1.14 4.58 1.74 -1.13 4.69
ACNLE 2.30 1.62 -0.88 5.49 2.33 -0.88 5.64
ACNE 2.37 1.50 -0.57 5.30 2.39 -0.57 5.45
FLE 2.38 1.44 -0.44 5.20 2.40 -0.44 5.33
PA 2.72 1.43 -0.09 5.52 2.75 -0.09 5.67
ANT 2.46 1.47 -0.41 5.34 2.49 -0.41 5.48
DBT 2.88 1.43 0.08 5.68 2.92 0.08 5.84
FLU 1.69 1.42 -1.09 4.47 1.70 -1.09 4.57
PYR 2.44 1.43 -0.35 5.24 2.47 -0.35 5.38
BAA 1.77 1.45 -1.08 4.61 1.78 -1.07 4.72
CHR 2.36 1.43 -0.43 5.16 2.39 -0.43 5.29
BAP 2.28 1.43 -0.53 5.08 2.30 -0.52 5.21
BEP 3.17 1.43 0.37 5.97 3.22 0.37 6.15
BBKF 1.72 1.41 -1.04 4.49 1.74 -1.03 4.59
DBAHA 2.00 1.41 -0.75 4.76 2.02 -0.75 4.87
PER 1.58 1.43 -1.22 4.37 1.59 -1.21 4.47
BGHIP 2.45 1.42 -0.33 5.23 2.48 -0.33 5.37
ICDP 2.52 1.42 -0.25 5.30 2.56 -0.25 5.44
Irish Sea NAP -0.23 1.04 -2.28 1.82 -0.23 -2.25 1.83
ACNLE -0.68 1.13 -2.90 1.54 -0.68 -2.86 1.55
ACNE -0.24 1.11 -2.42 1.94 -0.24 -2.39 1.96
FLE -0.68 1.07 -2.77 1.42 -0.67 -2.73 1.43
PA 0.23 1.05 -1.83 2.29 0.23 -1.82 2.31
ANT -1.33 1.07 -3.42 0.77 -1.32 -3.36 0.77
DBT -0.67 1.11 -2.84 1.51 -0.66 -2.80 1.52
FLU -0.34 1.03 -2.36 1.69 -0.34 -2.33 1.70
PYR -0.23 1.03 -2.25 1.79 -0.23 -2.23 1.81
BAA -0.99 1.04 -3.03 1.06 -0.98 -2.99 1.06
CHR -0.39 1.05 -2.46 1.67 -0.39 -2.43 1.69
BAP -0.23 1.05 -2.29 1.83 -0.23 -2.27 1.85
BEP -0.69 1.05 -2.75 1.36 -0.69 -2.71 1.37
BBKF -1.71 1.05 -3.78 0.36 -1.70 -3.71 0.36
DBAHA -1.34 1.05 -3.40 0.72 -1.33 -3.35 0.72
PER -3.18 1.05 -5.24 -1.12 -3.13 -5.10 -1.11
BGHIP -0.15 1.04 -2.19 1.89 -0.15 -2.17 1.91
ICDP -0.04 1.04 -2.07 1.99 -0.04 -2.05 2.01
Gulf of Cadiz PA -1.37 1.04 -3.41 0.68 -1.36 -3.35 0.68
ANT 1.69 1.12 -0.52 3.89 1.70 -0.52 3.97
FLU -4.22 1.06 -6.30 -2.14 -4.13 -6.10 -2.12
PYR -4.28 1.08 -6.40 -2.17 -4.19 -6.20 -2.15
BAA -1.23 1.11 -3.40 0.95 -1.22 -3.35 0.95
CHR -2.00 1.09 -4.13 0.13 -1.98 -4.05 0.13
BAP -2.27 1.08 -4.39 -0.16 -2.25 -4.29 -0.16
BEP -2.21 1.12 -4.40 -0.02 -2.18 -4.31 -0.02
BBKF -1.04 1.09 -3.18 1.11 -1.03 -3.13 1.12
DBAHA -0.59 1.20 -2.95 1.78 -0.58 -2.90 1.79
BGHIP -2.01 1.10 -4.18 0.15 -1.99 -4.09 0.15
ICDP -4.23 1.06 -6.29 -2.16 -4.14 -6.10 -2.13


Regional status compared to the Effects Range Low (ERL)

The first 4 tabs show:

  • Regional status: estimates of the regional status (averaged over PAHs) with pointwise upper one-sided 95% confidence intervals
  • Regional status by PAH: regional status estimates for each PAH with pointwise upper one-sided 95% confidence intervals plotted by PAH
  • Regional status by region: regional status estimates for each PAH with pointwise upper one-sided 95% confidence intervals plotted by region
  • Individual time series: the status estimates from the individual time series on which the regional assessment is based. The estimates for regions with insufficient geographic spread (and excluded from the meta-analysis) are also shown.

The symbols in all the plots have the following interpretation:

  • blue: the mean concentration is significantly (p < 0.05) below the Background Assessment Concentration (BAC)
  • green: the mean concentration is significantly (p < 0.05) below the ERL
  • red: the mean concentration is not significantly below the ERL

Note that to see why some points might be coloured blue, it is necessary to look ahead to the next section, where status is compared to the BAC. Unfortunately, it isn’t easy to present the assessment of status relative to the ERL with the assessment of status relative to the BAC because the ratio of the BAC to the ERL can vary between PAHs and regions.

The final 2 tabs show the estimates of the regional status (averaged over PAHs) and the regional status trend by PAH


Regional status


Regional status by PAH


Regional status by region


Individual time series


Estimates by region

This table shows the estimates of the regional status:

  • status: the estimated mean log concentration minus the log ERL
  • se: the corresponding standard error
  • upper: the corresponding pointwise upper one-sided 95% confidence limit
  • concentration: the estimated mean concentration relative to the ERL
  • conc upper: the corresponding pointwise upper one-sided 95% confidence limit


MIME region status se upper concentration conc upper
Northern North Sea -3.28 0.31 -2.77 0.04 0.06
Southern North Sea -2.25 0.24 -1.86 0.11 0.16
Channel -2.43 0.25 -2.02 0.09 0.13
Irish and Scottish West Coast -3.10 0.42 -2.42 0.04 0.09
Irish Sea -1.75 0.33 -1.21 0.17 0.30
Celtic Sea -2.89 0.53 -2.02 0.06 0.13
Iberian Sea -2.71 0.27 -2.26 0.07 0.10
Gulf of Cadiz -4.85 0.29 -4.37 0.01 0.01


Estimates by region and PAH

This table shows the estimates of the regional trend by determinand:

  • status: the estimated mean log concentration minus the log ERL
  • se: the corresponding standard error
  • upper: the corresponding pointwise upper one-sided 95% confidence limit
  • concentration: the estimated mean concentration relative to the ERL
  • conc upper: the corresponding pointwise upper one-sided 95% confidence limit


MIME region PAH status se upper concentration conc upper
Northern North Sea NAP -3.23 0.25 -2.82 0.04 0.06
PA -2.50 0.25 -2.09 0.08 0.12
ANT -3.30 0.25 -2.89 0.04 0.06
DBT -4.61 0.25 -4.21 0.01 0.01
FLU -3.19 0.25 -2.79 0.04 0.06
PYR -3.48 0.25 -3.07 0.03 0.05
BAA -2.98 0.25 -2.57 0.05 0.08
CHR -3.13 0.25 -2.72 0.04 0.07
BAP -3.07 0.25 -2.66 0.05 0.07
BGHIP
ICDP
Southern North Sea NAP -1.74 0.14 -1.51 0.18 0.22
PA -1.51 0.14 -1.28 0.22 0.28
ANT -1.87 0.14 -1.64 0.15 0.19
DBT -3.42 0.17 -3.14 0.03 0.04
FLU -2.16 0.14 -1.93 0.12 0.14
PYR -2.54 0.14 -2.31 0.08 0.10
BAA -2.14 0.14 -1.91 0.12 0.15
CHR -2.36 0.14 -2.13 0.09 0.12
BAP -2.51 0.14 -2.28 0.08 0.10
BGHIP
ICDP
Channel NAP -2.93 0.17 -2.64 0.05 0.07
PA -1.79 0.16 -1.52 0.17 0.22
ANT -1.98 0.17 -1.71 0.14 0.18
DBT -3.89 0.20 -3.56 0.02 0.03
FLU -2.23 0.17 -1.96 0.11 0.14
PYR -2.57 0.16 -2.30 0.08 0.10
BAA -1.94 0.16 -1.67 0.14 0.19
CHR -2.28 0.16 -2.00 0.10 0.13
BAP -2.26 0.16 -1.99 0.10 0.14
BGHIP
ICDP
Irish and Scottish West Coast NAP -2.82 0.38 -2.19 0.06 0.11
PA -2.18 0.38 -1.56 0.11 0.21
ANT -3.31 0.38 -2.69 0.04 0.07
DBT -4.08 0.38 -3.46 0.02 0.03
FLU -3.11 0.38 -2.49 0.04 0.08
PYR -3.45 0.38 -2.83 0.03 0.06
BAA -2.99 0.38 -2.37 0.05 0.09
CHR -2.92 0.38 -2.30 0.05 0.10
BAP -3.08 0.38 -2.46 0.05 0.09
BGHIP
ICDP
Irish Sea NAP -1.55 0.28 -1.09 0.21 0.34
PA -0.92 0.27 -0.47 0.40 0.62
ANT -1.66 0.27 -1.20 0.19 0.30
DBT -3.11 0.28 -2.66 0.04 0.07
FLU -1.71 0.27 -1.26 0.18 0.28
PYR -1.89 0.27 -1.44 0.15 0.24
BAA -1.48 0.27 -1.03 0.23 0.36
CHR -1.69 0.27 -1.24 0.18 0.29
BAP -1.74 0.27 -1.29 0.18 0.28
BGHIP
ICDP
Celtic Sea NAP -2.58 0.51 -1.74 0.08 0.18
PA -2.12 0.50 -1.29 0.12 0.27
ANT -2.66 0.51 -1.81 0.07 0.16
DBT -4.16 0.52 -3.31 0.02 0.04
FLU -2.78 0.50 -1.95 0.06 0.14
PYR -3.15 0.51 -2.32 0.04 0.10
BAA -2.65 0.51 -1.82 0.07 0.16
CHR -2.93 0.51 -2.10 0.05 0.12
BAP -3.01 0.51 -2.17 0.05 0.11
BGHIP
ICDP
Iberian Sea PA -1.92 0.19 -1.60 0.15 0.20
ANT -2.20 0.19 -1.88 0.11 0.15
FLU -2.54 0.19 -2.22 0.08 0.11
PYR -2.93 0.19 -2.61 0.05 0.07
BAA -2.41 0.19 -2.09 0.09 0.12
CHR -2.97 0.20 -2.65 0.05 0.07
BAP -2.94 0.19 -2.62 0.05 0.07
BGHIP
ICDP
Gulf of Cadiz PA -4.32 0.22 -3.95 0.01 0.02
ANT -4.32 0.23 -3.95 0.01 0.02
FLU -4.86 0.22 -4.49 0.01 0.01
PYR -5.06 0.22 -4.69 0.01 0.01
BAA -4.40 0.22 -4.03 0.01 0.02
CHR -4.98 0.22 -4.61 0.01 0.01
BAP -4.91 0.22 -4.54 0.01 0.01
BGHIP
ICDP


Regional status compared to the Background Assessment Concentration (BAC)

The first 4 tabs show:

  • Regional status: estimates of the regional status (averaged over PAHs) with pointwise upper one-sided 95% confidence intervals
  • Regional status by PAH: regional status estimates for each PAH with pointwise upper one-sided 95% confidence intervals plotted by PAH
  • Regional status by region: regional status estimates for each PAH with pointwise upper one-sided 95% confidence intervals plotted by region
  • Individual time series: the status estimates from the individual time series on which the regional assessment is based. The estimates for regions with insufficient geographic spread (and excluded from the meta-analysis) are also shown.

The symbols in all the plots have the following interpretation:

  • blue: the mean concentration is significantly (p < 0.05) below the BAC
  • green: the mean concentration is significantly (p < 0.05) below the Effects Range Low (ERL)
  • orange: there is no ERL and the mean concentration is not significantly below the BAC
  • red: the mean concentration is not significantly below the ERL

The final 2 tabs show the estimates of the regional status (averaged over PAHs) and the regional status trend by PAH


Regional status


Regional status by PAH


Regional status by region


Individual time series


Estimates by region

This table shows the estimates of the regional status:

  • status: the estimated mean log concentration minus the log BAC
  • se: the corresponding standard error
  • upper: the corresponding pointwise upper one-sided 95% confidence limit
  • concentration: the estimated mean concentration relative to the BAC
  • conc upper: the corresponding pointwise upper one-sided 95% confidence limit


MIME region status se upper concentration conc upper
Northern North Sea -0.34 0.27 0.11 0.71 1.11
Southern North Sea 0.40 0.19 0.72 1.50 2.06
Channel 0.29 0.21 0.64 1.34 1.89
Irish and Scottish West Coast -0.32 0.38 0.31 0.72 1.36
Irish Sea 0.96 0.30 1.44 2.60 4.22
Celtic Sea -0.19 0.50 0.63 0.83 1.88
Iberian Sea 1.23 0.23 1.62 3.43 5.05
Gulf of Cadiz -0.87 0.26 -0.45 0.42 0.64


Estimates by region and PAH

This table shows the estimates of the regional trend by determinand:

  • status: the estimated mean log concentration minus the log BAC
  • se: the corresponding standard error
  • upper: the corresponding pointwise upper one-sided 95% confidence limit
  • concentration: the estimated mean concentration relative to the BAC
  • conc upper: the corresponding pointwise upper one-sided 95% confidence limit


MIME region PAH status se upper concentration conc upper
Northern North Sea NAP -0.22 0.24 0.18 0.80 1.20
PA -0.45 0.24 -0.05 0.64 0.95
ANT -0.46 0.24 -0.06 0.63 0.94
DBT
FLU -0.44 0.24 -0.04 0.64 0.96
PYR -0.14 0.24 0.26 0.87 1.30
BAA -0.15 0.24 0.25 0.86 1.28
CHR -0.13 0.24 0.27 0.88 1.30
BAP -0.39 0.24 0.01 0.68 1.01
BGHIP -0.48 0.24 -0.08 0.62 0.92
ICDP -0.55 0.24 -0.15 0.58 0.86
Southern North Sea NAP 1.29 0.14 1.52 3.62 4.55
PA 0.51 0.14 0.73 1.66 2.08
ANT 0.96 0.14 1.19 2.62 3.28
DBT
FLU 0.57 0.14 0.79 1.76 2.21
PYR 0.77 0.14 1.00 2.17 2.71
BAA 0.65 0.14 0.87 1.91 2.39
CHR 0.60 0.14 0.82 1.82 2.27
BAP 0.15 0.14 0.38 1.16 1.46
BGHIP -0.61 0.14 -0.39 0.54 0.68
ICDP -0.84 0.14 -0.62 0.43 0.54
Channel NAP -0.01 0.17 0.27 0.99 1.31
PA 0.22 0.16 0.48 1.24 1.62
ANT 0.85 0.16 1.12 2.34 3.06
DBT
FLU 0.50 0.16 0.76 1.64 2.14
PYR 0.75 0.16 1.01 2.11 2.75
BAA 0.85 0.16 1.11 2.33 3.04
CHR 0.67 0.16 0.94 1.96 2.56
BAP 0.41 0.16 0.68 1.51 1.96
BGHIP -0.59 0.16 -0.32 0.55 0.72
ICDP -0.69 0.16 -0.43 0.50 0.65
Irish and Scottish West Coast NAP 0.12 0.37 0.73 1.13 2.08
PA -0.19 0.37 0.42 0.83 1.51
ANT -0.58 0.37 0.03 0.56 1.03
DBT
FLU -0.42 0.37 0.19 0.66 1.21
PYR -0.18 0.37 0.42 0.83 1.52
BAA -0.25 0.37 0.35 0.78 1.42
CHR 0.02 0.37 0.63 1.02 1.87
BAP -0.45 0.37 0.16 0.64 1.17
BGHIP -0.61 0.37 -0.01 0.54 0.99
ICDP -0.71 0.37 -0.10 0.49 0.90
Irish Sea NAP 1.43 0.27 1.87 4.19 6.52
PA 1.07 0.27 1.51 2.93 4.55
ANT 1.15 0.27 1.59 3.17 4.93
DBT
FLU 1.00 0.27 1.44 2.72 4.22
PYR 1.41 0.27 1.85 4.11 6.38
BAA 1.27 0.27 1.71 3.56 5.53
CHR 1.23 0.27 1.67 3.41 5.29
BAP 0.91 0.27 1.36 2.50 3.88
BGHIP 0.12 0.27 0.56 1.13 1.75
ICDP -0.05 0.27 0.39 0.95 1.48
Celtic Sea NAP 0.45 0.51 1.28 1.57 3.61
PA -0.09 0.49 0.72 0.92 2.06
ANT 0.15 0.51 0.98 1.16 2.66
DBT
FLU -0.05 0.49 0.76 0.95 2.13
PYR 0.14 0.49 0.95 1.15 2.59
BAA 0.11 0.49 0.92 1.12 2.52
CHR -0.00 0.49 0.81 1.00 2.25
BAP -0.36 0.50 0.46 0.70 1.59
BGHIP -1.11 0.50 -0.30 0.33 0.74
ICDP -1.13 0.49 -0.31 0.32 0.73
Iberian Sea PA 1.58 0.19 1.89 4.85 6.61
ANT 1.66 0.19 1.98 5.28 7.22
FLU 1.20 0.19 1.51 3.32 4.53
PYR 1.15 0.19 1.46 3.16 4.30
BAA 1.21 0.19 1.52 3.36 4.59
CHR 0.86 0.19 1.18 2.37 3.27
BAP 1.01 0.19 1.32 2.75 3.75
BGHIP 1.06 0.19 1.37 2.88 3.93
ICDP 1.03 0.19 1.34 2.80 3.82
Gulf of Cadiz PA -0.84 0.22 -0.48 0.43 0.62
ANT -0.45 0.22 -0.09 0.64 0.92
FLU -1.13 0.22 -0.77 0.32 0.46
PYR -0.98 0.22 -0.62 0.37 0.54
BAA -0.79 0.22 -0.43 0.45 0.65
CHR -1.12 0.22 -0.76 0.33 0.47
BAP -0.96 0.22 -0.60 0.38 0.55
BGHIP -0.75 0.22 -0.39 0.47 0.68
ICDP -1.11 0.22 -0.75 0.33 0.47


PAH profiles

The profile of concentrations across PAHs can sometimes be useful for examining sources, and for comparing concentrations across regions where assessment criteria are not available. Profiles were examined by fitting the following linear mixed model by restricted maximum likelihood:

  • response: mean log concentration
  • fixed model: region
  • random model: PAH + region.PAH + station + estimation variation + residual variation

where, for each time series, the mean log concentration is the estimated mean log concentration in the final monitoring year and the estimation variance is the variance associated with this estimate.

The estimated PAH concentration profiles by region are shown below with pointwise 95% confidence limits. The units are ug/kg dw normalised to 2.5% organic carbon (apart from in the Iberian Sea and Gulf of Cadiz where they are non-normalised).