Spatial variability in structural and functional aspects of macrofauna communities and their environmental parameters in the Jade Bay (Wadden Sea Lower Saxony, southern North Sea)
© Springer-Verlag and AWI 2012
Received: 20 September 2011
Accepted: 18 April 2012
Published: 5 May 2012
Spatial distribution and functional structure of intertidal benthic macrofauna in relation to environmental variables in the Jade Bay (southern North Sea) were studied and compared with other intertidal areas of the Wadden Sea. A total of 128 stations covering the whole Jade Bay were sampled in summer 2009. A total of 114 taxa were found. Highest species numbers occurred in the subtidal areas, whereas highest mean abundances were found in the upper intertidal areas. Based on species abundance data, six significantly distinct macrofauna communities in the Jade Bay were identified and evaluated with multivariate statistics, univariate correlations and canonical correspondence analysis. Differences in these community patterns were caused by the response of the dominant species (Hydrobia ulvae, Tubificoides benedii, Pygospio elegans, Caulleriella killariensis, Scoloplos armiger, Urothoe poseidonis, Microprotopus maculatus) to prevailing environmental conditions along the gradient from the lower and exposed sandy intertidal areas via intermediate mixed sediments to the upper mudflat areas. Distribution patterns in relation to tidal zonation were best explained by variability in submergence time, Chlorophyll a (chl a) content and sediment composition (mud content), which are proxies for hydrodynamic conditions and food availability. Species inventory and species richness were comparable with other intertidal areas of the Wadden Sea, but the Jade Bay differs from these areas regarding dominant species. Differences in sediment composition and morphological characteristics (macrotidal versus mesotidal Wadden Sea areas) are discussed for comparison of regional differences.
Tidal flat ecosystems are dynamic and complex habitats. They represent highly productive areas for macrobenthic organisms (Heip et al. 1995), thus being used as important feeding and nursery grounds by fishes and shore birds (Kuipers 1977; Widdows et al. 2004). Macrofauna species inhabiting tidal flats are highly adapted to this dynamic environment characterized by tidal fluctuations, changing water currents, wave action, altering periods of submersion and exposure to air from low to high intertidal levels (Newell 1970; Hummel et al. 1994). The hydrodynamic force itself influences sediment characteristics (grain size, organic carbon and chl a contents), topography (intertidal height) and salinity (Snelgrove and Butman 1994; Ysebaert et al. 2003; Widdows et al. 2004; Gray and Elliott 2009). It is well known that the prevailing environmental conditions affect the distribution of tidal flat macrofauna species, resulting in a specific zonation pattern of the living communities (Linke 1939; Michaelis 1987; Beukema 1988; Hertweck 1994; Reise et al. 1994; Ysebaert et al. 2003; van Colen et al. 2009; Puls et al. 2011). In addition, several macrofauna species themselves, such as Arenicola marina or Hydrobia ulvae (e. g. Flach 1992; Orvain et al. 2006), act as bioturbators, which increase sediment erosion, turbidity, nutrient and carbon processing (Herman et al. 1999). In contrast, bio-stabilisers such as seagrass beds, mussel beds and microphytobenthos mats can modify the habitat by reducing turbidity, increasing light penetration and enhancing sedimentation processes (Kröncke 1996; Bergfeld 1999; Widdows and Brinsley 2002; Bos et al. 2007; van Katwijk et al. 2010).
The tidal flat systems in the European Wadden Sea along the coastal region of the North Sea between Den Helder (Netherlands) and Skallingen (Denmark) are classified according to their tidal range to define the geomorphological, hydrological and sedimentological environment (Diekmann et al. 1987). Depending on the tidal range, tidal flats are distinguished as macrotidal (>3.5 m), mesotidal (1–3.50 m) and microtidal (<1 m) (Hayes 1979; Flemming 2011). Mesotidal flats are bordered by barrier islands, whereas macrotidal flats are open to sea where barrier islands are no longer present (Diekmann et al. 1987). With a mean tidal range of 3.8 m (Götschenberg and Kahlfeld 2008), the Jade Bay belongs to a macrotidal flat system, but differs from other macrotidal areas due to its land-locked bay character. In contrast to other sandier mesotidal flats (e.g. Balgzand and Königshafen), the Jade Bay is characterized by a large mudflat belt with a unique, muddier sediment composition.
The spatial, temporal and seasonal variability of macrofauna species of European intertidal flat systems have been intensively studied (e.g. Beukema 1974, 1988; Jensen 1992; Essink et al. 1998; van der Graaf et al. 2009). Especially in the German intertidal Wadden Sea the classification in mud, mixed sediments or sandflats/Arenicola-sandflats, and mussel or seagrass beds for macrofauna studies is still common (e.g. Linke 1939; Michaelis 1987; Hertweck 1994; Reise et al. 1994; Steuwer 2002). In the present study, we have used a quantitative statistical approach to identify macrofauna community patterns, which was also applied for different subtidal areas in the East Frisian Wadden Sea (Reiss and Kröncke 2001; Nehmer and Kröncke 2003) and the Pomeranian Bay (southern Baltic Sea) (Glockzin and Zettler 2008) as well as for intertidal areas of the Belgian and Dutch Wadden Sea (van Colen et al. 2009)
The objectives in the present study were (1) to describe the present spatial distribution and functional structure of the macrofauna communities in the Jade Bay by applying a statistical approach, (2) to investigate their relations to environmental factors by using univariate and multivariate statistical approaches in order to explain the spatial distribution patterns of species and structuring forces for zonation and (3) to compare the Jade Bay with intertidal macrofauna communities of other Wadden Sea areas.
Materials and methods
During high tide, the water surface area is 158 km², whereas at low tide it decreases to 44 km² (H. Westphal unpublished data). Furthermore, the study area contains intertidal Mytilus edulis mussel beds, which covered an area of 2.07 km² in 2009 (G. Millat unpublished data), as well as Zostera noltii seagrass beds, which were mostly found in the eastern and south-western areas and contained 10.9 km² in 2008 (Adolph 2010).
Sampling and laboratory analysis
Macrofauna sampling design
In the intertidal zone, stations were located according to the historical investigations by Linke (1939) in the 1930s and Michaelis (1987) in the 1970s. The study area was divided into different transects (Fig. 1) covering the whole Jade Bay. A total of 128 stations were sampled from April to July 2009 during low tide. Based on maps given by Linke (1939) and Michaelis (1987) the station coordinates were matched as closely as possible. Sampling positions were located by using a hand-held GPS receiver (GPS 72, GARMIN, WGS 84). At each station, 5 replicates were taken with a 10-cm-diameter cylindrical corer (79 cm² total surface) to a depth of 30 cm. A total of 642 samples were taken. Samples were sieved over 0.5 mm mesh size and fixed in 4 % buffered formaldehyde-seawater solution. In the laboratory, samples were stained with Bengal rose, and all macrofauna was sorted and identified to the lowest possible taxonomic level and then counted.
One surface sediment sample (<2 cm depth) was taken in parallel for analysing sediment characteristics (mud content <63 μm, median grain size), total organic carbon and chlorophyll a contents. Grain size distribution was analysed on fresh, ultrasonically suspended samples by means of Fritsch Analysette 22 laser granulometer. Sediment types were classified according to Kramer et al. (1994) based on differences in median grain size. The total organic carbon content (TOC) was determined according to Babu et al. (1999). The difference between total carbon and inorganic carbon was measured with a CS 500 IR analyser (Eltra, Neuss, Germany) and a CO2 coulometer (UIC, Joliet, USA), respectively. Chl a was analysed by high-performance liquid chromatography (HPLC) by the extraction from 5 g of homogenized sediment with 3 ml Dichlormethan-methanol (9 + 1). Chl a contents were calculated in mg/m² according to Brotas et al. (2007). Mean submergence time (given in hours per day) for each station was estimated using the finite-volume model FV COM (Lettmann et al. 2010).
Prior to data analysis, species abundance data were standardized to individuals per m²; species number is given per sample (79 cm2). Species such as Balanidae, which are typical for hard substrate, were excluded from the analysis. Each species was classified into functional feeding groups (surface deposit feeders, subsurface deposit feeders, suspension feeders, omnivores, predators, interface feeders and sand lickers based on available literature (Fauchald and Jumars 1979; Ysebaert et al. 2003; van Colen et al. 2009).
Data were first tested for normal distribution (D´Agostino and Pearson omnibus normality test). Multivariate techniques on fourth root transformed abundance data were used to analyse the community structure. Group average cluster analysis based on Bray-Curtis similarity followed by a similarity profile test (SIMPROF; Clarke and Gorley 2006) and multidimensional scaling (MDS) were calculated. ANOSIM randomization test (Clarke and Green 1988) was used to test significant spatial differences in community structure.
Discriminating species, which significantly (p < 0.05) separated the different communities, were identified using the similarity percentage routine (SIMPER). Grouped clusters were characterized by their mean abundance, species richness, Shannon-Wiener diversity index, eveness index, characteristic species, feeding type composition in terms of abundance and environmental characteristics. Statistical analysis were performed by using GraphPad Prism 5, the multivariate analysis were performed using The Plymouth Routines in multivariate ecological research (PRIMER) package, version 6 (Clarke and Gorley 2006).
Multivariate canonical correspondence analysis (CCA) by using CANOCO 4 software (Microcomputer Power) was performed to detect the variation in the species data set and the relationship between species composition and measured environmental variables. Forward selection was used to determine and rank the importance of environmental variables, which in turn would be used for determining the species composition (Ter Braak and Verdonshot 1995). Monte Carlo permutation test was performed for statistical validity.
Relationships between species abundance and the environmental variables were also examined by using univariate (Spearman rank) correlations.
General characterization of the macrofauna in the Jade Bay
A total of 114 taxa (45 polychaetes, 34 crustaceans, 20 molluscs, 6 oligochaetes and 9 others) were found during the study period. Number of taxa in a single sample (79 cm²) varied between 0 and 32 taxa, with a mean of 12 taxa. Total abundance varied between 0 and 480,380 ind./m², with a mean of 56,203 ± 3,291 ind./m². Gastropods (87 %) were the most abundant taxonomic group, followed by oligochaetes (6 %), polychaetes (6 %) and crustaceans (2 %).
The Jade Bay was generally dominated by Hydrobia ulvae (69 %), Tubificoides benedii (12 %), Caulleriella killariensis (5 %), Pygospio elegans (3 %) and Scoloplos armiger (3 %), contributing 92 % of the total macrofauna abundance. Excluding H. ulvae, species such as T. benedii (29 %), P. elegans (15 %), Retusa obtusa (13 %), C. killariensis (11 %), Cerastoderma edule (6 %), S. armiger (5 %) and Macoma balthica (4 %) were most abundant and made up 83 % of the total macrofauna abundance. Trophic structure in terms of mean abundance was dominated by surface deposit feeders (89 %) followed by subsurface deposit feeders (7 %).
Community structure of the macrofauna
Summary of the six communities, their habitat types and environmental variables
partially covered with seagrass beds
Mixed sediments 1
Mixed sediments 2
(n = 39)
(n = 43)
(n = 15)
(n = 13)
(n = 14)
(n = 4)
Median grain size (μm)
Mud content (%)
Chl a (mg/m²)
Submergence time (h)
S (total/79 cm²)
S (mean/79 cm²)
Mean abundance (ind./m²)
Characteristic species (mean abundance ind./m²)
Surface deposit feeders
Subsurface deposit feeders
The Hydrobia ulvae community (Community 1) included 39 stations (mean similarity of 66 %) located in the eastern and western intertidal flats. The mostly very fine sand environment was characterized by moderate submergence time (8.5 h) and mud, chl a and TOC contents. This community was characterized by a species-rich benthic fauna, highest mean abundance but lowest Shannon-Wiener diversity and eveness. The dominant species was H. ulvae (99 %), which was widespread in high densities in the Jade Bay (Fig. 4f). A clustered pattern with maximum abundance of more than 370,000 ind./m² was found in close proximity to Zostera noltii seagrass beds in the eastern and south-western parts. Due to its dominance, the community structure was characterized by surface deposit feeders. By exclusion of H. ulvae, species such as C. edule (Fig. 3a), M. balthica (Fig. 3b), Retusa obtusa (Fig. 4e), P. elegans (Fig. 2e), Nephtys hombergii (Fig. 2c) and T. benedii dominated the community structure. This was the only community with high abundances of bivalves. Scrobicularia plana (Fig. 3c), commonly known for muddy sediments, occurred only at 29 stations, which were restricted to the southern and south-western upper intertidal areas.
Forty-three stations (mean similarity of 32 %) belonged to the upper intertidal mudflats (Community 2) arranged belt-shaped with lowest median grain size and submergence time (6 h), but with highest mud, TOC and chl a contents. The community revealed the second highest mean abundance and moderate values of Shannon-Wiener diversity and eveness, respectively. The community was numerically dominated by the subsurface deposit feeding T. benedii (Fig. 3f), which reached abundances of more than 40,000 ind./m², and again H. ulvae. Higher abundances were also found for P. elegans, Corophium volutator (Fig. 4b), R. obtusa, Heteromastus filiformis (Fig. 3d) and Juveniles of Nereis spp.. In general, the omnivorous feeding ragworm Hediste diversicolor (Fig. 2b) occurred widespread in this muddier upper intertidal area. The dissimilarity between Community 2 and Communities 1, 3-6 was mainly caused by the two dominant species T. benedii and P. elegans. Classified as subsurface and surface deposit feeders, both species presented also the dominant feeding modes of this community.
Communities 3 and 4 both were referred to as mixed sediments. Community 3 included 15 stations (mean similarity of 33 %), which were widely distributed over the lower intertidal flats of the Jade Bay area predominantly found along the margins of the tidal gullies and tidal channels. Six stations in the eastern Jade Bay were located near Mytilus beds. Community 4 was summarized by 13 stations (mean similarity of 23 %) located in the southern part of the Jappensand (central Jade Bay) and along the tidal inlets in the western part of the Jade Bay. The environmental characteristics of both communities were rather similar according to moderate mud, chl a and TOC contents, but differed in median grain size. Submergence time varied between 10.7 and 14.5 h, respectively. In Community 3, species richness and mean abundance were low, but Shannon-Wiener diversity and eveness high. The community structure was numerically dominated by interface feeding polychaetes, such as C. killariensis (Fig. 3e) and P. elegans, as well as H. ulvae. These three dominant species, which contributed with 85 % to the total macrofauna abundance, were also responsible for the distinction with the five other communities. Surface deposit and interface feeders dominated the community structure in terms of feeding modes. Community 4 was characterized by the lowest species richness and the lowest mean abundance in contrast to the five others. Trophic structure in terms of abundance was dominated by subsurface deposit-, surface deposit- and interface feeders, and predators caused by the characteristic species S. armiger, H. ulvae, C. killariensis and R. obtusa. S. armiger (Fig. 2d) represented the most and widespread polychaete of the lower intertidal area, reaching a maximum abundance of 1,519 ind./m².
The sandflat community (Community 5) included 14 stations (mean similarity of 15 %) characterized by the lowest mud, chl a and TOC contents in the Jade Bay, but highest median grain size. Submergence time was 10.4 h, respectively. Stations were found in the lower intertidal central part of the Jade Bay (Jappensand) and on sandflats along the adjacent channels Stenkentief and Vareler Fahrwasser, which consisted of very fine to medium sands. A total of 26 species were found. Mean abundance, Shannon-Wiener diversity and eveness revealed moderate values. Characteristic species were H. ulvae and S. armiger, several sand-licking amphipods (U. poseidonis, Bathyporeia sarsi; Fig. 4a, c) and Juveniles of the polychaete Nephtys spp.. The two amphipod species occurred exclusively in Community 5, and were attributed to the differences between Community 5 and Communities 1–4 and 6. Additionally, these species contributed to the high percentage of sandlickers in this community. Surface deposit feeders were again the dominant feeding mode.
The subtidal community (Community 6) (mean similarity of 17 %) was located within the adjacent channels (Ahne, Vareler Fahrwasser) characterized by 24 h of submergence time, sediments with high median grain size, low mud, chl a and TOC contents. This community was characterized by the highest species richness, Shannon-Wiener diversity and eveness. The mean abundance was relatively low with 3,991 ± 2,547 ind./m² similar to Community 3. S. armiger, the amphipod M. maculatus (Fig. 4d) and the polychaete Lanice conchilega were the dominant species. Other species numerically dominating in this community were polychaetes, such as P. elegans, Aricidea minuta, Eteone longa (Fig. 2f) and Phyllodoce mucosa, accompanied by R. obtusa and the pantopod Nymphon brevirostre. Dominant feeding modes were subsurface deposit- and interface feeders, followed by predators. Surface deposit feeders were of minor importance in contrast to the five other communities.
Relationship between spatial distribution and environmental variables
Species (e.g. L. conchilega, M. maculatus, S. armiger, A. minuta) related to a longer submergence time were positively correlated to the first axis. Species such as B. sarsi and U. poseidonis, which were positively correlated to the second axis, related to a higher median grain size. In contrast, species such as C. volutator, T. benedii and H. filiformis related to more muddy sediment, were negatively correlated.
Spearman rank correlations between species densities and environmental variables (sediment fractions: >63 μm, <63 μm (mud) expressed as percentage, median grain size (μm), total carbon (%), TOC (%), Chl a (mg/m²) and submergence time (h)
>63 μm (%)
<63 μm% (mud)
Median grain size (μm)
Total carbon (%)
Chl a (mg/m²)
Relationships between spatial distribution of macrofauna communities and environmental variables
The spatial distribution of macrofauna communities in the Jade Bay was best explained by the variability in submergence time, sediment composition (mud content, grain size) and chl a content, all following an intertidal gradient. Our results agree with several studies in other areas of the Wadden Sea that sediment characteristics, food availability (expressed as TOC and/or chl a contents) and intertidal height (submergence time) are the primary causal factors for the distinct distribution patterns of intertidal macrofauna species resulting in characteristic zonation patterns (Beukema 1988; Junoy and Viéitez 1990; Ysebaert et al. 2003; van Colen et al. 2009). Salinity, which is also a major factor affecting the macrofauna community structure in estuaries (Ysebaert et al. 1998, 2003; Ysebaert and Herman 2002), was of rather minor importance in the whole polyhaline Jade Bay with no major freshwater influence.
Although the Spearman rank correlation coefficients were weakly correlated, results agree well with the species positions in the CCA biplot, which confirms the relationships between the macrofauna communities and the predominant gradients. Similarly, weak relationships were found by van Colen et al. (2009), who explained this with the high small-scale patchiness of benthic organisms. According to their sensitivity to withstand desiccation during increasing times of emersion towards the upper muddy intertidal areas and to cope with increasing hydrodynamics towards the more exposed, lower sandy intertidal areas, species are correlated to the two main gradients in the CCA plot. The general decrease from subsurface deposit and deposit feeders towards interface feeders in terms of abundance with increasing submergence time (decreasing tidal height), reflects also the close relationship between food availability and/or food quality for the macrofauna. Our results (Table 1) confirm that higher sedimentation rates result in nutrient and TOC-rich sediments in the upper intertidal areas, in contrast to the sandflats and subtidal areas, where higher hydrodynamics prevent organic enrichment. These relationships were previously described by Pearson and Rosenberg (1978) and Hertweck (1994).
The different statistical methods used revealed similar results and divided the Jade Bay into almost the same zones with characteristic macrofauna community compositions corresponding with different habitat types (subtidal vs. intertidal, sandflat, mixed sediments, mudflat, seagrass beds). These results are in accordance with the still common sediment related classification of intertidal flats in the Wadden Sea (see “Introduction”).
In addition to the above discussed role of environmental parameters on macrofauna community structure, biological interactions (predation, competition), as well as bio-stabilization, bio-irrigation and bioturbation of species can modify the habitat (Schubert and Reise 1986; Rosenberg 1997; van der Meer et al. 2000; Reise 2002; Orvain et al. 2004, 2006) and, thus, have probably an influence on the distribution patterns of species and the structure of the benthic communities in the Jade Bay. Furthermore, no significant differences in community structure caused by seasonal fluctuations primarily triggered by species-specific differences in the timing of recruitment (ANOSIM test, unpublished data) were found between spring and early summer.
Comparison of Jade Bay communities with other Wadden Sea areas
In general, the intertidal Jade Bay was characterized by species such as H. ulvae, T. benedii, C. killariensis, P. elegans, S. armiger, U. poseidonis, R. obtusa, C. edule, M. balthica, C. volutator and H. diversicolor. Although, differences exist in sampling design (area, mesh size), the observed species composition is very similar compared with other areas of the Wadden Sea (e.g. Dörjes 1978; Beukema 1988; Jensen 1992; Hertweck 1994; Reise et al. 1994; van Colen et al. 2009). Reise et al. (1994) found 14 species per sample in Königshafen, while Beukema (1988) for Balgzand and van Colen et al. (2009) for the Zwin nature reserve observed both an average of 11 species per sample. In accordance with these studies, species richness in the Jade Bay revealed an average of 12 species per sample.
Comparison of the dominant species ranked according to highest species abundance found per m² on intertidal mudflats, mixed sediments, sandflats and seagrass beds in different areas of the Wadden Sea
Jade Bay (present study)
Balgzand (Beukema 1988)
Königshafen (Reise et al.1994)
Mesh size/season dominant species (ind./m²)
0.5 mm/late spring/early summer 2009
1 mm/early spring
1 mm/sping and autumn 1990*
Comparison of the coverage of intertidal mudflats, mixed sediments, sandflats and seagrass beds of the entire intertidal (area in %) between different areas of the Wadden Sea
The upper intertidal mudflats (Community 2) were dominated by the pelophilous subsurface-deposit feeding oligochaete T. benedii, common for muddy intertidal flats (Dörjes et al. 1969; Vöge et al. 2008). Supported by the significant correlation between abundance and mud content, this opportunistic species with a short life span, high productivity, small body size and an invasive ability (Pearson and Rosenberg 1978) is adapted to these sediments rich in organic carbon and able to tolerate anoxic conditions (Giere et al. 1999). A similar autecological potential has been recorded for the opportunistic P. elegans, which is well adapted to survive in moderate hypoxic conditions (Kube and Powilleit 1997). Mud offers a favourable substratum for diatoms (de Jonge 1985), expressed by the highest chl a content in this community. This might have explained the high abundance of species such as H. ulvae, P. elegans and Juveniles of Corophium spp., which feed on benthic diatoms (Fenchel et al. 1975; Hertweck 1994).
A comparison between different mudflat areas in the Wadden Sea shows that the numerical dominance of various small, r-selected and opportunistic species is common (Table 3). On Balgzand, the upper muddy coastal areas were characterized by a Corophium volutator-zone dominated by C. volutator and H. ulvae (Beukema 1988, Beukema and Cadée 1997). In contrast, M. balthica, Tubificoides pseudogaster, P. elegans, Tharyx spp. and H. filiformis were dominant species in Königshafen (Reise et al. 1994). However, a Corophium volutator-zone had occurred there in the 1930s (Reise et al. 2008), as also described for the Jade Bay in the 1930s (Linke 1939). Erosion and sedimentation processes, which have been linked to sea level rise, are suggested to be major causes of these observed long-term changes (Reise et al. 2008, Schückel and Kröncke submitted). In contrast to Balgzand and Königshafen, the present high densities of T. benedii in the Jade Bay were related to the higher areal extent of mudflats (43 % of the intertidal area (Meyer and Ragutzki 1999), Table 4), which contained a mean percentage of 76 % sediment mud content. These regional differences of tidal flats in the Wadden Sea are caused by differences of tidal range (macrotidal vs. mesotidal), which affect morphological characteristics (back-barrier island, land-locked bay or open bay) and sedimentary processes (Bungenstock and Weerts 2010). Sedimentation processes in the tidal flats along the Dutch and the German North Sea coasts are mainly controlled (90 %) by the sediment supply from the North Sea (Beets and van der Spek 2000; Hoselmann and Streif 2004). In contrast, Königshafen is predominated by medium sands originating from surrounding dunes (Reise et al. 1994). However, the comparison between the mesotidal Balgzand (Dutch Wadden Sea) and the macrotidal Jade Bay (German Wadden Sea) mudflat areas reveal that size and mud content varied largely. On Balgzand, only 20 % (Table 4) of the intertidal area contained sediments with a mean percentage of only 10–24 % mud content (Beukema 1974, 1988). This is probably caused by a northward transport of sand brought into the back-barrier basin derived from the fronting barrier (Beets and van der Spek 2000), resulting in the dominance of sandy sediments (Beukema 1988). Higher percentages of mud content in the sediments affect the spatial distribution of macrofauna species (Sakamaki and Nishimura 2009). An increase in mud content causes an increase in sediment cohesiveness resulting in an impermeable and high organically concentrated sediment layer, which is harmful for sensitive species suffer from anoxic and sulfidic conditions (Sakamaki and Nishimura 2009), but seems to be advantageous for T. benedii in the Jade Bay.
Two communities typical for mixed sediments were found in the Jade Bay. According to Beukema (1988) and supported by our results, these areas were defined as intermediate tidal flats characterized by mixed sediments and moderate submergence times, where most species of the mudflat and sandflat areas co-occurred. However, the cluster analysis revealed two significant distinct communities different in terms of species richness, mean abundance and species dominance. Differences were mainly attributed to the presence or absence of C. killariensis, P. elegans, H. ulvae and S. armiger. In Community 3 the high abundances of P. elegans’ biogenic structures may have additionally affected the macrofauna. The increase of the psammophilous orbiinid polychaete S. armiger in Community 4 indicated an increase in hydrodynamical stress, expressed by the somewhat higher median grain size and longer submersion times. Stronger tidal currents at stations of this community are also confirmed by the lower TOC content resulting in the lowest macrofauna densities and species richness.
The intermediate tidal flats on Balgzand were, in a similar manner, numerically dominated by H. ulvae (Table 3), but also by suspension-feeding bivalves (C. edule, M. edulis and Mya arenaria) (Beukema 1988), which declined since the 1990s (Beukema and Dekker 2005). As mentioned above for the mudflats, there is evidence that these differences in species composition and abundance between our present and earlier studies (in the Jade Bay and other parts of the Wadden Sea) were related to long-term changes. For instance, several studies discussed the absence of severe winters in the last two decades to be responsible for the lower recruitment success of bivalve species in different areas of the Wadden Sea (e.g. Strasser et al. 2003; Beukema and Dekker 2005; van der Graaf et al. 2009), as well as observed in the Jade Bay (Schückel and Kröncke submitted).
The sandflats (Community 5) located in the central part (Jappensand) and at the margins of the northern entrance of the Jade Bay were dominated by the opportunistic S. armiger and mobile amphipods such as U. poseidonis, B. sarsi and B. pilosa. As commonly known from exposed beaches, S. armiger and the fast swimming and digging amphipods are typical well-adapted inhabitants of unstable, well oxygenated coarser sandy sediments (pers. observation; Degraer et al. 2003), which is supported by the CCA plot. This explains why corresponding community structures were found at the lower and exposed marginal areas of Balgzand (Beukema 1988) and the Lütetsburger Plate (Steuwer 2002), on the sandflats in Spiekeroog (Grotjahn 1990) and Königshafen (Reise et al. 1994), as well as on the emerging flood delta shoals between the islands of Rømø and Sylt (Lackschewitz and Reise 1998). Especially in Königshafen, the large contribution of sandflats (76 %, Table 4) and the predominance of medium sands might be a structuring factor resulting in higher abundances of S. armiger in this part of the Wadden Sea than in the Jade Bay. Higher densities of Bathyporeia were also found in the transition area between intertidal and subtidal, as also reported for the Knechtsand by Michaelis (1969).
The subtidal community (Community 6) was dominated by the cosmopolitan S. armiger in terms of abundance. Adapted to mobile sediments, this species indicated a higher hydrodynamic stress and tidal currents. Population studies on genetic divergence, reproductive strategies and tolerance to hypoxia of S. armiger in subtidal and intertidal areas indicated a sibling species: the intertidal “Type I” and the subtidal “Type S” (Kruse et al. 2003, 2004). Thus, it might be the case that Community 6 was dominated by the subtidal “Type S”, but Communities 4 and 5 by the intertidal “Type I”. Further genetic studies are necessary to clear this assumption. In addition, the tube-building polychaetes L. conchilega and P. elegans, as well the amphipod M. maculatus, which is common for subtidal flats, occurred in high abundances. Comparing this community structure with other shallow subtidal areas, an almost similar species composition was found in the exposed sandy inlets of the Wichter Ee (Nehmer and Kröncke 2003) and the Otzumer Balje (Reiss and Kröncke 2001) in the East Frisian Wadden Sea.
Of course, our communities revealed a large overlap of dominant species in the different communities (Table 1), which is commonly known from other studies (Reise et al. 1994; Ysebaert and Herman 2002; van Colen et al. 2009; Puls et al. 2011), indicating the wide range of ecological tolerance of these species along the intertidal gradient. However, our results revealed significant differences between the macrofauna communities in the Jade Bay compared with other, more homogenous Wadden Sea areas. The main reason for these differences is the large variation in the proportion of mudflats between areas and the higher sediment mud contents in most of the Jade Bay mudflats resulting, for instance, in the dominance of T. benedii. The macrotidal and special bay shaped character of the Jade Bay area seems to cause the regional-specific sedimentation processes that in turn influence the sediment composition, submergence time and food availability for the macrofauna communities, in comparison with other Wadden Sea areas.
We thank Helmo Nicolai and his team of the Institute for Chemistry and Biology of the Marine Environment (ICBM) for their help with fieldwork and providing ship time throughout the year. A special thank goes to Malte Groh and Jonas Arlt for their untiring help in the field. We thank several students and Christian März, Jörn Logemann and Kerstin Wintermann for help with sampling and sorting samples. We are grateful to Gerald Millat for providing aerial images, Karsten Lettmann and Oliver Bleich for providing submergence time data and Martina Schulz for technical assistance. We are also grateful to Karsten Reise and two anonymous reviewers for very helpful comments on the manuscript. Our results are part of the “Jade Bay project”, an interdisciplinary background study addressing to collect historically and new basic data of natural sciences and cultural studies for a coastal database of Lower Saxony. This research was funded by the Ministry for Science and Culture of Lower Saxony (“Niedersächsisches Vorab der Volkswagen Stiftung”).
- Adolph W (2010) Praxistest Monitoring Küste 2008. Seegraskartierung - Gesamtbestandserfassung der eulitoralen Seegrasbestände im Niedersächsichen Wattenmeer und Bewertung nach EG-Wasserrahmenrichtlinie. NLWKN Küstengewässer und Ästuare 2/2010, p 52Google Scholar
- Babu CP, Brumsack H-J, Schnetger B (1999) Distribution of organic carbon in surface sediments along the eastern Arabian sea: a revisit. Mar Geol 162:91–103View ArticleGoogle Scholar
- Beets DJ, van der Spek AJF (2000) The Holocene evolution of the barrier and back-bariier basins of Belgium and the Netherlands as a function of late Weichselian morphology, relative sea-level rise and sediment supply. Geologie en Mijnbouw/Neth J Geosci 79:3–16Google Scholar
- Behre K-H (2004) Coastal development sea-level change and settlement history during the later Holocene in the Clay District of Lower Saxony (Niedersachsen) northern Germany. Quatern Int 112:37–53View ArticleGoogle Scholar
- Bergfeld C (1999) Macrofaunal community pattern in an intertidal sandflat: effects of organic enrichment via biodeposition by mussel beds. First Results. Senck Marit 29:23–27View ArticleGoogle Scholar
- Beukema JJ (1974) Seasonalchanges in the biomass of the macro-benthos of a tidal flat area I the Dutch Wadden Sea. Neth J Sea Res 8(1):94–107View ArticleGoogle Scholar
- Beukema JJ (1988) An evaluation of the ABC-method (abundance/biomass comparison) as applied to macrozoobenthic communities living on tidal flats in the Dutch Wadden Sea. Mar Biol 99:425–433View ArticleGoogle Scholar
- Beukema JJ, Cadée GC (1997) Local differences in macrozoobenthic response to enhanced food supply caused by mild eutrophication in a Wadden Sea area: food is only locally a limiting factor. Limnol Oceanogr 42(6):1424–1435View ArticleGoogle Scholar
- Beukema JJ, Dekker R (2005) Decline of recruitment success in cockles and other bivalves in the Wadden Sea: possible role of climate change, predation on postlarvae and fisheries. Mar Ecol Prog Ser 287:149–167View ArticleGoogle Scholar
- Bos AR, Bouma TJ, de Kort GLJ, van Katwijk MM (2007) Ecosystem engineering by annual intertidal seagrass beds: sediment accretion and modification. Estuar Coast Shelf Sci 74:344–348View ArticleGoogle Scholar
- Brotas V, Mendes CR, Cartaxana P (2007) Microphytobenthic biomass assessment by pigment analysis: comparison of spectrophotometry and High Performance Liquid Chromatography methods. Hydrobiologia 587:19–24View ArticleGoogle Scholar
- Bungenstock F, Weerts HJT (2010) The high-resolution Holocene sea-level curve for Northwest Germany: global signals, local effects or data-artefacts? Int J Earth Sci (Geol Rundsch) 99:1687–1706View ArticleGoogle Scholar
- Cardoso PG, LillebØ AI, Pardal MA, Ferreira SM, Marques JC (2002) The effect of different primary producers on Hydrobia ulvae populations: a case study in a temperate intertidal estuary. J Exp Mar Biol Ecol 277:173–195View ArticleGoogle Scholar
- Cardoso PG, Raffaelli D, Pardal MA (2007) Seagrass beds and intertidal invertebrates: an experimental test of the role of habitat structure. Hydrobiologia 575:221–230View ArticleGoogle Scholar
- Clarke KR, Gorley RN (2006) Primer v6: user manual/tutorial Plymouth UK: Primer-E Plymouth Marine LaboratoryGoogle Scholar
- Clarke KR, Green RH (1988) Statistical design and analysis for a “biological effects” study. Mar Ecol Prog Ser 173:121–127Google Scholar
- De Jonge VN (1985) The occurrence of “epipsammic” diatom populations: a result of interaction between physical sorting of sediment and certain properties of diatom species. Estuar Coast Shelf Sci 21:607–622View ArticleGoogle Scholar
- Degraer S, Volckaert A, Vincx M (2003) Macrobenthic zonation patterns along a morphodynamical continuum of macrotidal low tide bar/rip and ultra-dissipative sandy beaches. Estuar Coast Shelf Sci 56:459–468View ArticleGoogle Scholar
- Diekmann R, Osterthun M, Partenscky H-W (1987) Influence of water-level elevation and tidal range on the sedimentation in a German tidal flat area. Prog Oceanogr 15:151–166View ArticleGoogle Scholar
- Dörjes J (1978) Das Watt als Lebensraum. In: Reineck H-E (ed) Das Watt—Ablagerungs- und Lebensraum, Senckenberg Buch 50, W Kramer Frankfurt am Main, pp 107–143Google Scholar
- Dörjes J, Gadow S, Reineck H-E, Singh IB (1969) Die Rinnen der Jade (südliche Nordsee) Sediment und Makrobenthos. Senck Marit 50:5–62Google Scholar
- Essink K, Beukema JJ, Madsen, PB, Michaelis H, Vedel GR (1998) Long-term development of biomass of intertidal macrozoobenthos in different parts of the Wadden Sea. Governed by nutrient loads? Senckenb Marit 29(1/6):25–35Google Scholar
- Fauchald K, Jumars PA (1979) The diet of worms: a study of polychaete feeding guilds. Oceanogr Mar Biol Ann Rev 17:193–284Google Scholar
- Fenchel T, Kofoed LH, Lappalainen A (1975) Particle size-selection of two deposit feeders: the amphipod Corophium volutator and the Prosobranch Hydrobia ulvae. Mar Biol 30:19–128View ArticleGoogle Scholar
- Flach EC (1992) The influence of four macrozoobenthic species on the abundance of the amphipod Corophium volutator on tidal flats of the Wadden Sea. Neth J Sea Res 29:379–394View ArticleGoogle Scholar
- Flemming BW (2011) Siliciclastic back-barrier tidal flats. In: Davis RA Jr, Dalrymple RW (eds) Tidal depositional systems. Springer, Berlin, pp 1–48Google Scholar
- Giere O, Preusse JH, Dubilier N (1999) Tubificoides benedii (Tubificidae Oligochaeta)—a pioneer in hypoxic and sulfidic environments. An overview of adaptive pathways. Hydrobiologia 406:235–241View ArticleGoogle Scholar
- Glockzin M, Zettler ML (2008) Spatial macrozoobenthic patterns in relation to major environmental factors—a case study from the Pomeranian Bay (southern Baltic Sea)Google Scholar
- Götschenberg A, Kahlfeld A (2008) The Jade. Die Küste 74:263–274Google Scholar
- Gray JS, Elliott M (2009) Ecology of marine sediments—from science to management. 2nd Edition. Oxford University Press, p 225 Google Scholar
- Grotjahn M (1990) Sedimente und Makrofauna der Watten bei der Insel Spiekeroog -Untersuchungen im Rahmen des “Sensitivitätsrasters Deutsche Nordseeküste”. -Jber. 1987, Jber. Forsch.—Stelle Küste 39:97–119Google Scholar
- Haubois A-G, Guarini J-M, Richard P, Blanchard GF, Sauriau P-G (2002) Spatio-temporal differentiation in the population structure of Hydrobia ulvae on an intertidal mudflat (Marennes-Oléron Bay France). J Exp Mar Biol Ecol 371:1–12Google Scholar
- Hayes MO (1979) Barrier island morphology as a function of tidal and wave regime. In: Leatherman SP (ed) Barrier islands, Academic Press, New York, pp 1–27Google Scholar
- Heip CHR, Goosen NK, Herman PMJ, Kromkamp J, Middelburg JJ, Soetaert K (1995) Production and consumption of biological particles in temperate tidal estuaries. Oceanogr Mar Biol Ann Rev 33:1–149Google Scholar
- Herkül K, Kotta J (2009) Effects of eelgrass (Zostera marina) canopy removal and sediment addition on sediment characteristics and benthic communities in the Northern Baltic Sea. Mar Biol 30:74–82Google Scholar
- Herman PMJ, Middelburg JJ, van de Koppel J, Heip CHR (1999) Ecology of estuarine macrobenthos. Adv Ecol Res 29:195–240View ArticleGoogle Scholar
- Herman PMJ, Middelburg JJ, Widdows J, Lucas CH, Heip CHR (2000) Stable isotop labeling experiments confirm the importance of mirophytobenthos as food for macrofauna. Mar Ecol Prog Ser 204:79–92View ArticleGoogle Scholar
- Hertweck G (1994) Zonation of benthos and Lebensspuren in the tidal flats of the Jade Bay, southern North Sea. Senckenb Marit 24(1/6):157–170Google Scholar
- Hoselmann C, Streif H (2004) Holocene sea-level rise and its effect on the mass balance of coastal deposits. Quat Int 112:89–103View ArticleGoogle Scholar
- Hummel H, Fortuin AW, Bogaards RH, Meijboom A, de Wolf L (1994) The effects of prolonged emersion and submersion by tidal manipulation on marine macrobenthos. Hydrobiologia 282(283):219–234View ArticleGoogle Scholar
- Jensen KT (1992) Macrozoobenthos on an intertidal mudflat in the Danish Wadden Sea: comparison of surveys made in the 1930 s, 1940 s and 1980 s. Helgol Meeresunters 46:363–376View ArticleGoogle Scholar
- Junoy J, Viéitez JM (1990) Macrobenthic community structure in the Ría de Foz an intertidal estuary (Galicia North West Spain). Mar Biol 107:329–339View ArticleGoogle Scholar
- Kramer KMJ, Brockmann UH, Warwick RM (1994) Tidal estuaries: manual of sampling and analytical procedures. AABalkema, Rotterdam, p 303Google Scholar
- Kröncke I (1996) Impact of biodeposition on macrofaunal communities in intertidal sandflats. Mar Ecol 17:159–174View ArticleGoogle Scholar
- Kruse I, Reusch TBH, Schneider MV (2003) Sibling species or poecilogony in the polychaete Scolopos armiger? Mar Biol 142:937–947Google Scholar
- Kruse I, Strasser M, Thiermann F (2004) The role of ecological divergence in speciation between intertidal and subtidal Scoloplos armiger (Polychaeta: Orbiniidae) J Sea 51:53–62Google Scholar
- Kube J, Powilleit M (1997) Factors controlling the distribution of Marenzellaria cf viridis Pygospio elegans and Streblospio shrubsoli (Polychaeta: Spionida) in the southern Baltic Sea with special attention for the response to an event of hypoxia. Aquat Ecol 31:187–198View ArticleGoogle Scholar
- Kuipers BR (1977) On the ecology of juvenile plaice on a tidal flat in the Wadden Sea Neth. J Sea Res 11(1):56–91View ArticleGoogle Scholar
- Lackschewitz D, Reise K (1998) Macrofauna on flood delta shoals in the Wadden Sea with an underground association between the lugworm Arenicola marina and the ampipod Urothoe poseidonis. Helgol Meeresunters 52:147–158View ArticleGoogle Scholar
- Lettmann K, Wolff J-O, Liebezeit G, Meier G (2010) Investigation of the spreading and dilution of domestic waste water inputs into a tidal bay using the finite-volume model FVCOM. Abstract for EGU Conference 2010, ViennaGoogle Scholar
- Linke O (1939) Die Biota des Jadebusenwattes. Helgol Meeresunters 1:201–348View ArticleGoogle Scholar
- Meyer C, Ragutzki G (1999) KFKI Forschungsvorhaben Sedimentverteilung als Indikator für morphodynamische Prozesse. Dienstbericht der Forschst Küste Norderney 21:1–43Google Scholar
- Michaelis H (1969) Makrofauna und Vegetation der Knechtsandwatten. -Jber.1967, Forsch.-Stelle f. Insel- u. Küstenschutz 19, Norderney, pp 147–175Google Scholar
- Michaelis H (1987) Bestandsaufnahme des eulitoralen Makrobenthos im Jadebusen in Verbindung mit einer Luftbild-Analyse. Jber 1986 Forschst Küste Norderney 38:13–97Google Scholar
- Nehmer P, Kröncke I (2003) Macrofaunal communities in the Wichter Ee a Channel System in the East Frisian Wadden Sea. Senckenb Marit 32(1/2):1–10View ArticleGoogle Scholar
- Newell RC (1970) Biology of intertidal animals. Logos Press London, p 555Google Scholar
- Orvain F, Sauriau P-G, Sygut A, Joassard L, Hir LP (2004) Interacting effects of Hydrobia ulvae bioturbation and microphytobenhos on the erodibility of mudflat sediments. Mar Ecol Prog Ser 278:205–223View ArticleGoogle Scholar
- Orvain F, Sauriau P-G, Bacher C, Prineau M (2006) The influence of sediment cohesiveness on bioturbation effects due to Hydrobia ulvae on the initial erosion of intertidal sediments: a study combining flume and model approaches. J Sea Res 55:54–73View ArticleGoogle Scholar
- Pearson TH, Rosenberg R (1978) Macrobenthic succession in relation to organic enrichment and pollution of the marine environment. Oceanogr Mar Biol Ann Rev 16:229–311Google Scholar
- Puls W, van Bernem K-H, Eppel D, Kapitza H, Pleskachevsky A, Riethmüller R, Vaessen B (2011) Prediction of benthic community structure from environmental variables in a soft-sediment tidal basin (North Sea). Helgol Mar Res. doi:10.1007/s10152-011-0275-y Google Scholar
- Reise K (2002) Sediment mediated species interactions in coastal waters. J Sea Res 48:127–141View ArticleGoogle Scholar
- Reise K, Herre E, Sturm M (1994) Biomass and abundance inintertidal sediments of Königshafen in the northern Wadden Sea. Helgol Meeresunters 48:201–215View ArticleGoogle Scholar
- Reise K, Herre E, Sturm M (2008) Mudflat biota since the 1930 s: change beyond return? Helgol Mar Res 62:13–22View ArticleGoogle Scholar
- Reiss H, Kröncke I (2001) Spatial and temporal distribution of macrofauna in the Otzumer Balje (East Frisian Wadden Sea). Senckenb Marit 31(2):283–298View ArticleGoogle Scholar
- Rosenberg GG (1997) Bioresuspension and biodeposition: a review. J Mar Syst 11:269–278View ArticleGoogle Scholar
- Sakamaki T, Nishimura O (2009) Is sediment mud content a significant predictor of macrobenthos abundance in low-mud-content tidal flats? Mar Freshw Res 60:160–167View ArticleGoogle Scholar
- Schubert A, Reise K (1986) Predatory effects of Nephtys hombergii on other polychaetes in tidal flats sediments. Mar Ecol Prog Ser 34:117–124View ArticleGoogle Scholar
- Schückel U, Kröncke I (submitted) Long-term changes in the macrofauna communities in the Jade Bay (Wadden Sea of Lower Saxony, southern North Sea). A comparison between the 1930s, 1970s and 2009. Estuar Coast Shelf SciGoogle Scholar
- Snelgrove PVR, Butman CA (1994) Animal-sediment relationships revisited: cause versus effect. Oceanogr Mar Biol Ann Rev 32:111–177Google Scholar
- Steuwer J (2002) Die Makrofauna der Lütetsburger Plate. TMAP-Monitoring 2001. Forschst Küste Norderney 7:1–11Google Scholar
- Strasser M, Dekker R, Essink K, Günther C-P, Jaklin S, Kröncke I, Madsen PB, Michaelis H, Vedel G (2003) How predictable is high bivalve recruitment in the Wadden Sea after a severe winter. J Sea Res 49:47–57View ArticleGoogle Scholar
- Ter Braak CFJ, Verdonshot PFM (1995) Canonical correspondence analysis and related multivariate methods in aquatic ecology. Aquat Sci 57(3):255–289View ArticleGoogle Scholar
- van Colen C, Snoeck F, Struyf K, Vincx M, Degraer S (2009) Macrobenthic community structure and distribution in the Zwin nature reserve (Belgium and the Netherlands). J Mar Biol Ass UK 89(3):431–438View ArticleGoogle Scholar
- Van der Graaf S, de Vlas J, Herlyn M, Voss J, Heyer K, Drent J (2009) Macrozoobenthos. Thematic Report No 10. In: Marencic H, de Vlas J (eds) Quality Status Report 2009. Wadden Sea Ecosystem No 25. Common Wadden Sea Secretariat, Trilateral Monitoring and Assessment Group, Wilhelmshaven, Germany, p 27Google Scholar
- Van der Meer J, Beukema JJ, Dekker R (2000) Population dynamics of two marine polchaetes: the relative role of density dependence predation and winter conditions. ICES J Mar Sci 57:1488–1494View ArticleGoogle Scholar
- Van Katwijk MM, Bos AR, Hermus DCR, Suykerbuyk W (2010) Sediment modification by seagrass beds: muddification and sandification induced by plant cover and environemental conditions. Estuar Coast Shelf Sci 89:175–181View ArticleGoogle Scholar
- Vöge S, Reiss H, Kröncke I (2008) Macrofauna succession in an infilling salt marsh clay pit. Senckenb Marit 38:93–106View ArticleGoogle Scholar
- Widdows J, Brinsley M (2002) Impact of biotic and abiotic processes on sediment dynamics and the consequences to the structure and functioning of the intertidal zone. J Sea Res 48:143–156View ArticleGoogle Scholar
- Widdows J, Blauw A, Heip CHR, Herman PMJ, Lucas CH, Middelburg JJ, Schmidt S, Brinsley MD, Twisk F, Verbeek H (2004) Role of physical and biological processes in sediment dynamics of a tidal flat in Westerschelde Estuary SW Netherlands. Mar Ecol Prog Ser 274:41–56View ArticleGoogle Scholar
- Ysebaert T, Herman PMJ (2002) Spatial and temporal variation in benthic macrofauna and relationships with environmental variables in an estuarine intertidal soft-sediment environment. Mar Ecol Prog Ser 244:105–124View ArticleGoogle Scholar
- Ysebaert T, Meire P, Coosen J, Essink K (1998) Zonation of intertidal macrobenthos in the estuaries of Schelde and Ems. Aquat Ecol 32:53–71View ArticleGoogle Scholar
- Ysebaert T, Herman PMJ, Meire P, Craeymeersch J, Verbeek H, Heip CHR (2003) Large-scale spatial patterns in estuaries: estuarine macrobenthic communities in the Schelde estuary NW Europe. Estuar Coast Shelf Sci 57:335–355View ArticleGoogle Scholar