This document describes a regional assessment of the status and trends of PAH concentrations in sediment in Scottish marine waters. It builds on the 2015 assessment of data from the UK’s Clean Seas Environmental Monitoring Programme conducted by the Clean Safe Seas Evidence Group. In this, assessments of status and trends were made for a large number of PAH time series, each of a single PAH at a single monitoring station. The regional assessment synthesises the results of the individual time series to assess status and trends at the biogeographic regional level.
The regional assessment considers the following PAHs: PA, ANT, FLU, PYR, BAA, CHR, BAP, BGHIP, ICDP.
The regional assessment synthesises data from monitoring stations that are in estuarine, coastal and offshore waters. The results might therefore differ from regional assessments for Marine Strategy Framework Directive purposes, which exclude estuarine stations.
The following sections:
A time series of PAH concentrations is assessed for status if:
The conditions are more stringent for trends. Specifically, a time series is assessed for trends if:
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. 3 of the 360 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:
And the shapes have the following meaning:
The final set of tabs give:
Note that there are no ERLs for BGHIP and ICDP. Time series for these PAHs are only assessed against the BAC, and are coloured blue if significantly below the BAC and orange otherwise.
Number of time series with upwards, downwards or no trend by biogeographic region
region | status | PA | ANT | FLU | PYR | BAA | CHR | BAP | BGHIP | ICDP | total |
---|---|---|---|---|---|---|---|---|---|---|---|
Northern North Sea | upward trend | 0 | 0 | 3 | 2 | 3 | 2 | 2 | 2 | 2 | 16 |
no trend | 20 | 14 | 15 | 16 | 14 | 16 | 14 | 14 | 16 | 139 | |
downward trend | 0 | 1 | 2 | 2 | 1 | 0 | 2 | 2 | 0 | 10 | |
Scottish Continental Shelf | upward trend | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
no trend | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | |
downward trend | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 9 | |
Minches & W Scotland | upward trend | 0 | 0 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 3 |
no trend | 6 | 4 | 6 | 5 | 5 | 5 | 4 | 6 | 6 | 47 | |
downward trend | 1 | 1 | 1 | 1 | 0 | 1 | 1 | 1 | 1 | 8 | |
Irish Sea | upward trend | 1 | 0 | 1 | 3 | 0 | 2 | 0 | 0 | 0 | 7 |
no trend | 10 | 10 | 8 | 7 | 11 | 6 | 9 | 7 | 8 | 76 | |
downward trend | 1 | 1 | 2 | 2 | 1 | 3 | 3 | 5 | 4 | 22 | |
total | upward trend | 1 | 0 | 4 | 6 | 4 | 4 | 3 | 2 | 2 | 26 |
no trend | 36 | 28 | 29 | 28 | 30 | 27 | 27 | 27 | 30 | 262 | |
downward trend | 3 | 4 | 6 | 6 | 3 | 5 | 7 | 9 | 6 | 49 |
Proportion of time series with upwards, downwards or no trend by biogeographic region
region | status | PA | ANT | FLU | PYR | BAA | CHR | BAP | BGHIP | ICDP | total |
---|---|---|---|---|---|---|---|---|---|---|---|
Northern North Sea | upward trend | 0 | 0 | 15 | 10 | 17 | 11 | 11 | 11 | 11 | 10 |
no trend | 100 | 93 | 75 | 80 | 78 | 89 | 78 | 78 | 89 | 84 | |
downward trend | 0 | 7 | 10 | 10 | 5 | 0 | 11 | 11 | 0 | 6 | |
Scottish Continental Shelf | upward trend | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
no trend | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | |
downward trend | 100 | 100 | 100 | 100 | 100 | 100 | 100 | 100 | 100 | 100 | |
Minches & W Scotland | upward trend | 0 | 0 | 0 | 14 | 17 | 0 | 16 | 0 | 0 | 5 |
no trend | 86 | 80 | 86 | 72 | 83 | 83 | 67 | 86 | 86 | 81 | |
downward trend | 14 | 20 | 14 | 14 | 0 | 17 | 17 | 14 | 14 | 14 | |
Irish Sea | upward trend | 9 | 0 | 9 | 25 | 0 | 18 | 0 | 0 | 0 | 7 |
no trend | 83 | 91 | 73 | 58 | 92 | 55 | 75 | 58 | 67 | 72 | |
downward trend | 8 | 9 | 18 | 17 | 8 | 27 | 25 | 42 | 33 | 21 | |
total | upward trend | 2 | 0 | 10 | 15 | 11 | 11 | 8 | 5 | 5 | 8 |
no trend | 90 | 88 | 74 | 70 | 81 | 75 | 73 | 71 | 79 | 78 | |
downward trend | 8 | 12 | 16 | 15 | 8 | 14 | 19 | 24 | 16 | 14 |
Number of time series with each status by biogeographic region
region | status | PA | ANT | FLU | PYR | BAA | CHR | BAP | BGHIP | ICDP | total |
---|---|---|---|---|---|---|---|---|---|---|---|
Northern North Sea | blue | 10 | 9 | 9 | 8 | 8 | 8 | 9 | 10 | 10 | 81 |
green | 8 | 5 | 11 | 11 | 11 | 12 | 10 | 1 | 7 | 76 | |
red | 2 | 6 | 0 | 1 | 1 | 0 | 1 | 9 | 3 | 23 | |
Scottish Continental Shelf | blue | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 9 |
green | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | |
red | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | |
Minches & W Scotland | blue | 4 | 5 | 5 | 4 | 3 | 3 | 6 | 6 | 6 | 42 |
green | 3 | 2 | 2 | 3 | 4 | 4 | 1 | 0 | 1 | 20 | |
red | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | |
Irish Sea | blue | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
green | 6 | 9 | 11 | 10 | 7 | 10 | 10 | 0 | 6 | 69 | |
red | 6 | 3 | 1 | 2 | 5 | 2 | 2 | 12 | 6 | 39 | |
total | blue | 15 | 15 | 15 | 13 | 12 | 12 | 16 | 17 | 17 | 132 |
green | 17 | 16 | 24 | 24 | 22 | 26 | 21 | 1 | 14 | 165 | |
red | 8 | 9 | 1 | 3 | 6 | 2 | 3 | 22 | 9 | 63 |
Proportion of time series with each status by biogeographic region
region | status | PA | ANT | FLU | PYR | BAA | CHR | BAP | BGHIP | ICDP | total |
---|---|---|---|---|---|---|---|---|---|---|---|
Northern North Sea | blue | 50 | 45 | 45 | 40 | 40 | 40 | 45 | 50 | 50 | 45 |
green | 40 | 25 | 55 | 55 | 55 | 60 | 50 | 5 | 35 | 42 | |
red | 10 | 30 | 0 | 5 | 5 | 0 | 5 | 45 | 15 | 13 | |
Scottish Continental Shelf | blue | 100 | 100 | 100 | 100 | 100 | 100 | 100 | 100 | 100 | 100 |
green | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | |
red | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | |
Minches & W Scotland | blue | 57 | 71 | 71 | 57 | 43 | 43 | 86 | 86 | 86 | 67 |
green | 43 | 29 | 29 | 43 | 57 | 57 | 14 | 0 | 14 | 32 | |
red | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 14 | 0 | 1 | |
Irish Sea | blue | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
green | 50 | 75 | 92 | 83 | 58 | 83 | 83 | 0 | 50 | 64 | |
red | 50 | 25 | 8 | 17 | 42 | 17 | 17 | 100 | 50 | 36 | |
total | blue | 38 | 38 | 38 | 32 | 30 | 30 | 40 | 42 | 43 | 37 |
green | 42 | 40 | 60 | 60 | 55 | 65 | 52 | 3 | 35 | 46 | |
red | 20 | 22 | 2 | 8 | 15 | 5 | 8 | 55 | 22 | 17 |
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.
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:
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:
The meta-analysis is restricted to time series from monitoring stations that are classified as ‘representative’. ‘Baseline’ stations (those with near pristine conditions or only subject to very diffuse inputs) are omitted because any trends there will likely be caused by different processes to those at representative stations. Similarly ‘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:
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:
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).
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.
This section provides more detail on the number and geographic spread of the time series available for the meta-analyses. The the stations with parametric trend and status assessments having excluded trend assessments at baseline and 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 or time series or spread of stations and are excluded from the regional assessment:
The first 4 tabs show:
The symbols in all these plots have the following interpretation:
The final 2 tabs show the estimates of the regional trend (averaged over PAHs) and the regional trend by PAH
This table shows the estimates of the regional trend:
region | trend | se | lower | upper | % yearly change | %yc lower | %yc upper |
---|---|---|---|---|---|---|---|
Northern North Sea | -2.22 | 4.60 | -11.24 | 6.81 | -2.19 | -10.63 | 7.04 |
Minches & W Scotland | 2.72 | 7.71 | -12.39 | 17.84 | 2.76 | -11.65 | 19.53 |
Irish Sea | 0.38 | 6.47 | -12.29 | 13.05 | 0.38 | -11.57 | 13.94 |
This table shows the estimates of the regional trend by determinand:
region | PAH | trend | se | lower | upper | % yearly change | %yc lower | %yc upper |
---|---|---|---|---|---|---|---|---|
Northern North Sea | PA | -2.41 | 4.57 | -11.36 | 6.54 | -2.38 | -10.74 | 6.76 |
ANT | -3.37 | 4.58 | -12.34 | 5.60 | -3.31 | -11.61 | 5.76 | |
FLU | -1.45 | 4.57 | -10.40 | 7.50 | -1.44 | -9.88 | 7.79 | |
PYR | -2.06 | 4.57 | -11.01 | 6.89 | -2.04 | -10.42 | 7.14 | |
BAA | -0.95 | 4.57 | -9.91 | 8.01 | -0.95 | -9.43 | 8.34 | |
CHR | -2.15 | 4.57 | -11.11 | 6.81 | -2.13 | -10.51 | 7.05 | |
BAP | -2.06 | 4.57 | -11.02 | 6.89 | -2.04 | -10.44 | 7.14 | |
BGHIP | -2.75 | 4.57 | -11.71 | 6.20 | -2.72 | -11.05 | 6.40 | |
ICDP | -2.76 | 4.57 | -11.72 | 6.19 | -2.72 | -11.06 | 6.39 | |
Minches & W Scotland | PA | 2.53 | 7.69 | -12.54 | 17.60 | 2.56 | -11.78 | 19.25 |
ANT | 1.57 | 7.70 | -13.51 | 16.65 | 1.58 | -12.64 | 18.12 | |
FLU | 3.49 | 7.69 | -11.58 | 18.56 | 3.55 | -10.94 | 20.40 | |
PYR | 2.89 | 7.69 | -12.19 | 17.96 | 2.93 | -11.47 | 19.67 | |
BAA | 3.99 | 7.69 | -11.09 | 19.07 | 4.07 | -10.49 | 21.01 | |
CHR | 2.79 | 7.69 | -12.28 | 17.87 | 2.83 | -11.56 | 19.57 | |
BAP | 2.88 | 7.69 | -12.20 | 17.95 | 2.92 | -11.48 | 19.67 | |
BGHIP | 2.19 | 7.69 | -12.89 | 17.26 | 2.21 | -12.09 | 18.84 | |
ICDP | 2.18 | 7.69 | -12.89 | 17.25 | 2.20 | -12.10 | 18.83 | |
Irish Sea | PA | 0.19 | 6.44 | -12.43 | 12.81 | 0.19 | -11.69 | 13.66 |
ANT | -0.77 | 6.44 | -13.41 | 11.86 | -0.77 | -12.55 | 12.59 | |
FLU | 1.14 | 6.44 | -11.48 | 13.77 | 1.15 | -10.84 | 14.76 | |
PYR | 0.54 | 6.44 | -12.08 | 13.16 | 0.54 | -11.38 | 14.07 | |
BAA | 1.64 | 6.44 | -10.98 | 14.27 | 1.66 | -10.40 | 15.34 | |
CHR | 0.45 | 6.44 | -12.18 | 13.07 | 0.45 | -11.47 | 13.97 | |
BAP | 0.53 | 6.44 | -12.09 | 13.16 | 0.53 | -11.39 | 14.06 | |
BGHIP | -0.16 | 6.44 | -12.78 | 12.46 | -0.16 | -12.00 | 13.27 | |
ICDP | -0.17 | 6.44 | -12.79 | 12.46 | -0.17 | -12.00 | 13.26 |
The first 4 tabs show:
The symbols in all the plots have the following interpretation:
Note that to see why some points are 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
This table shows the estimates of the regional status:
region | status | se | upper | concentration | conc upper |
---|---|---|---|---|---|
Northern North Sea | -2.28 | 0.19 | -1.97 | 0.10 | 0.14 |
Minches & W Scotland | -3.11 | 0.26 | -2.68 | 0.04 | 0.07 |
Irish Sea | -0.80 | 0.22 | -0.44 | 0.45 | 0.64 |
This table shows the estimates of the regional trend by determinand:
region | PAH | status | se | upper | concentration | conc upper |
---|---|---|---|---|---|---|
Northern North Sea | PA | -2.22 | 0.16 | -1.95 | 0.11 | 0.14 |
ANT | -2.72 | 0.16 | -2.45 | 0.07 | 0.09 | |
FLU | -2.77 | 0.16 | -2.51 | 0.06 | 0.08 | |
PYR | -2.94 | 0.16 | -2.68 | 0.05 | 0.07 | |
BAA | -2.49 | 0.16 | -2.23 | 0.08 | 0.11 | |
CHR | -2.66 | 0.16 | -2.40 | 0.07 | 0.09 | |
BAP | -2.71 | 0.16 | -2.45 | 0.07 | 0.09 | |
BGHIP | -0.55 | 0.16 | -0.28 | 0.58 | 0.76 | |
ICDP | -1.49 | 0.16 | -1.22 | 0.23 | 0.29 | |
Minches & W Scotland | PA | -2.89 | 0.25 | -2.47 | 0.06 | 0.08 |
ANT | -3.64 | 0.26 | -3.22 | 0.03 | 0.04 | |
FLU | -3.57 | 0.25 | -3.15 | 0.03 | 0.04 | |
PYR | -3.79 | 0.25 | -3.37 | 0.02 | 0.03 | |
BAA | -3.37 | 0.25 | -2.95 | 0.03 | 0.05 | |
CHR | -3.44 | 0.25 | -3.03 | 0.03 | 0.05 | |
BAP | -3.60 | 0.25 | -3.18 | 0.03 | 0.04 | |
BGHIP | -1.39 | 0.25 | -0.97 | 0.25 | 0.38 | |
ICDP | -2.29 | 0.25 | -1.87 | 0.10 | 0.15 | |
Irish Sea | PA | -0.59 | 0.20 | -0.26 | 0.55 | 0.77 |
ANT | -1.20 | 0.20 | -0.87 | 0.30 | 0.42 | |
FLU | -1.23 | 0.20 | -0.90 | 0.29 | 0.41 | |
PYR | -1.36 | 0.20 | -1.02 | 0.26 | 0.36 | |
BAA | -0.97 | 0.20 | -0.64 | 0.38 | 0.53 | |
CHR | -1.17 | 0.20 | -0.84 | 0.31 | 0.43 | |
BAP | -1.20 | 0.20 | -0.87 | 0.30 | 0.42 | |
BGHIP | 0.75 | 0.20 | 1.08 | 2.11 | 2.94 | |
ICDP | -0.21 | 0.20 | 0.12 | 0.81 | 1.13 |
The first 4 tabs show:
The symbols in all the plots have the following interpretation:
The final 2 tabs show the estimates of the regional status (averaged over PAHs) and the regional status trend by PAH
This table shows the estimates of the regional status:
region | status | se | upper | concentration | conc upper |
---|---|---|---|---|---|
Northern North Sea | -0.04 | 0.14 | 0.20 | 0.96 | 1.22 |
Minches & W Scotland | -0.86 | 0.23 | -0.49 | 0.42 | 0.62 |
Irish Sea | 1.45 | 0.18 | 1.74 | 4.26 | 5.71 |
This table shows the estimates of the regional trend by determinand:
region | PAH | status | se | upper | concentration | conc upper |
---|---|---|---|---|---|---|
Northern North Sea | PA | -0.19 | 0.16 | 0.07 | 0.83 | 1.07 |
ANT | 0.08 | 0.16 | 0.35 | 1.08 | 1.41 | |
FLU | -0.06 | 0.16 | 0.20 | 0.94 | 1.23 | |
PYR | 0.29 | 0.16 | 0.56 | 1.34 | 1.74 | |
BAA | 0.24 | 0.16 | 0.51 | 1.28 | 1.66 | |
CHR | 0.23 | 0.16 | 0.49 | 1.26 | 1.64 | |
BAP | -0.06 | 0.16 | 0.20 | 0.94 | 1.22 | |
BGHIP | -0.36 | 0.16 | -0.10 | 0.70 | 0.91 | |
ICDP | -0.52 | 0.16 | -0.26 | 0.60 | 0.77 | |
Minches & W Scotland | PA | -0.85 | 0.25 | -0.44 | 0.43 | 0.64 |
ANT | -0.86 | 0.25 | -0.45 | 0.42 | 0.64 | |
FLU | -0.86 | 0.25 | -0.45 | 0.42 | 0.64 | |
PYR | -0.60 | 0.25 | -0.19 | 0.55 | 0.83 | |
BAA | -0.66 | 0.25 | -0.25 | 0.52 | 0.78 | |
CHR | -0.58 | 0.25 | -0.17 | 0.56 | 0.84 | |
BAP | -0.95 | 0.25 | -0.54 | 0.39 | 0.58 | |
BGHIP | -1.14 | 0.25 | -0.73 | 0.32 | 0.48 | |
ICDP | -1.26 | 0.25 | -0.85 | 0.28 | 0.43 | |
Irish Sea | PA | 1.44 | 0.20 | 1.77 | 4.24 | 5.87 |
ANT | 1.59 | 0.20 | 1.92 | 4.92 | 6.81 | |
FLU | 1.48 | 0.20 | 1.81 | 4.40 | 6.10 | |
PYR | 1.86 | 0.20 | 2.19 | 6.44 | 8.92 | |
BAA | 1.75 | 0.20 | 2.08 | 5.77 | 7.99 | |
CHR | 1.71 | 0.20 | 2.04 | 5.53 | 7.66 | |
BAP | 1.45 | 0.20 | 1.78 | 4.27 | 5.91 | |
BGHIP | 0.96 | 0.20 | 1.29 | 2.62 | 3.63 | |
ICDP | 0.79 | 0.20 | 1.11 | 2.19 | 3.04 |
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:
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.