Trawling activity decreases sharply at both stations studied from 1991 to 2015 (fig. 2). Whereas at the beginning of the sampling period up to almost 150,000 hours per year were fished in the corresponding ICES rectangles (fig. 1), in recent years only about 20,000 hours per year are fished in the same areas with beam and otter trawls.
For the correlation of temperature with the number of benthic species, a stratified analysis of water temperatures in the summer months (June 01 – August 31) in all years at a section of the study area was performed with Ocean Data View 5.5.2 (fig. 3). In summer, the surface temperature is different from the one in deeper layers. At about 25 m, a thermocline was detected and for the following analyses, the temperature data from 25 m and deeper were used.
An evaluation of the mean water temperatures from 1991-2021 for the summer months of June-August for stations 40 (ICES transects 40F3, 40F4, 39F3, and 39F4) and station 36 (ICES transects 40F2, 40F3, 39F2 and 39F3), below 25 m depth, shows that the mean temperature was increasing about 4 °C with time (fig. 4).
3.1 Species composition
A total of 115 species were recorded during the years 1991-2021 from station 40, consisting of the phyla Chordata (31 species), Arthropoda (26 species), Mollusca (21 species), Cnidaria (11 species), Echinodermata (10 species), Annelida (8 species), Bryozoa (6 species) and Porifera (2 species). From station 36, a total of 72 species were identified during the sampling period 1991-2020. These belonged to the phyla Chordata (24 species), Mollusca (14 species), Arthropoda (12 species), Echinodermata (9 species), Cnidaria (4 species), Bryozoa (3 species), Annelida (2 species), Nemertea (1 species) and Porifera (1 species). For both station 36 and 40, the most common (in every year present) species were Arnoglossus laterna (Walbaum, 1792), Asterias rubens Linnaeus, 1758, Buglossidium luteum Chabanaud, 1930, Limanda limanda Linnaeus, 1758 and Pagurus bernhardus (Linnaeus, 1758) (fig. 5).
The species discovery (or so-called power-) curves of the two investigated stations showed rather similar results (fig. 6). While no new species have been recorded at station 40 since 2017, the species discovery curve at station 36 is not yet fully asymptotic, as additional species have also been identified in recent years. This difference comes mainly due to the high number of repetitive samples at station 40, as explained in the methods section. In spite of their difference in the number of (detected + new) species, both curves are showing a rather similar shape.
The number of individuals divided by the number of samples is higher at station 36 than at station 40 (fig. 7). However, the decreasing trend at station 36 is also steeper than at station 40.
3.2 Biodiversity
An analysis of the Shannon-Wiener biodiversity indices over the period 1991-2021 showed that the biodiversity index for both stations generally shows a slightly negative trend (fig. 8). However, in 2021 a strong increase of the biodiversity index can be seen again which was even slightly above the starting value from the year 1991 at station 36. The means of the two stations are not significantly different (p= 0.3837), indicating no difference in mean Shannon-Wiener biodiversity.
Correlation of Shannon-Wiener indices with trawling activity in the studied areas showed no statistically significant relationship. At both stations the determination of Pearson's correlation coefficient was clearly above the 0.05 mark (station 36: p=0.4068, station 40: p=0.1679).
A Bray-Curtis biodiversity analysis of the two stations over the whole sampling period shows that biodiversity changes over time at both stations, while both stations also differ from each other (fig. 9). The stress value with 0.08 is relatively low.
When correlating certain species with trawling intensity at station 40 it is noticeable that some species show a positive correlation with increasing fishing pressure, others relate negatively. The dominant, small flatfish species Buglossidium luteum and Arnoglossus laterna are among the species that show a negative correlation with increasing trawling activity. The opposite is observed for Asterias rubens and Pagurus bernhardus.
A comparison of the mean values with an ANOVA shows with a p-value of < 0.05 that the mean values of the temperatures of the two stations differ. A determination of the Spearman correlation coefficient clearly shows that the temperature increases within the years 1991-2021 (station 36: p=0.0001138, station 40: p=0.00379).
Dogger Bank sampling 2021
In 2021, the benthic fauna of Dogger Bank was surveyed at 25 stations (fig. 10). The stations were chosen within and close to the focus area of the MGF-North Sea project. Furthermore, the stations are located around the historically long-term sampled station 40. This was to ensure that a comparison of a station over time can be counted representatively for the Dogger Bank as it is assumed that the benthic fauna is comparatively uniform. Because of organisatoric difficulties, two different mesh sizes were used in the beam trawl in the stert. 14 stations were sampled with a mesh size of 0.8 cm, 11 stations with a mesh size of 1 cm as in the last decades (fig. 10).
A Bray-Curtis analysis showed that the benthic biodiversity of the sampled stations is close together except for two stations (fig. 11). The conspicuous stations 1 and 13 are both located at the outermost corners of the investigated area on Dogger Bank (fig. 10). The two clusters for the remaining stations come because two different mesh sizes at the cod end (0,8 mm at stations marked in red and 1,0 cm at stations marked in blue) were used in order to compare the beam trawl catching efficiency. Thus, the abundances varied because of the used methodology. The species observed overlapped 47 % at all stations, regardless of the different mesh sizes. However, it is noticeable that a mesh size of 0.8 cm ensures a higher species diversity in the catch, 38 % of all recorded species are only found in the catches with the mesh size of 0.8 cm, while only 15 % of the species were found in the usual 1 cm mesh sizes.
A comparison of the Shannon-Wiener indices (Hs) of the stations also expresses that station 1 as well as station 13 shows a slight difference in biodiversity compared to the stations with a lower Hs. In general, there are only very slight differences in the Hs of the stations studied (fig. 12).