Introduction

The Wadden Sea is a large intertidal ecosystem along the Eastern shore of the North Sea and extends from the north of The Netherlands to Denmark. It is the largest wetland area of Europe and contains the largest unbroken stretch of intertidal flats worldwide. In ecological terms, the Wadden Sea benthos is dominated by filter feeders (Baird et al. 2004), and in historic time oyster beds (Ostrea edulis) and Sabellaria reefs (Sabellaria spinulosa) were regarded as characteristic features of the tidal channels in the Wadden Sea (Hagmeier and Kändler 1927). Until the beginning of the twentieth century the native European oyster (Ostrea edulis) was widely spread in the Wadden Sea, but declined due to overfishing (Reise 1982). The last living bed of European oyster in the international Wadden Sea was found in 1940 (Hagmeier 1941).

The spread of the non-indigenous Pacific oyster (Crassostrea gigas) over The Netherlands started with its introduction for cultivation purposes in the Dutch Oosterschelde in 1964 (Drinkwaard 1999). The Pacific oyster occurs naturally in marine waters of Japan and southeast Asia. Although it was thought that the low water temperatures would prevent natural production in the Oosterschelde, between 1974 and 1976 the first massive spat fall took place. An exponential increase of wild populations occurred in later years, and from there it was spread all over the Dutch coastal waters (Dankers et al. 2006). Since the beginning of the 1980s wild Pacific oysters have been reported from the Dutch Wadden Sea, starting from Texel (Bruins 1983), where they were introduced (Tydeman 2008).

In 1996 a first settlement occurred in the Niedersachsen part of the German Wadden Sea area, which may have dispersed from The Netherlands by natural means (Wehrmann et al. 2000; Nehring 2003). In 1986 Pacific oysters were introduced in the Northern Wadden Sea at Sylt for cultivation (Reise 1998), and in 1991 the first oysters were found outside culture plots in the Northern Wadden Sea.

At the start of the new millennium it became clear that the Pacific oyster occurred throughout the Dutch Wadden Sea, and since 2004 Pacific oysters have occurred to such an extent that yearly increases in standing stock can be roughly quantified on the basis of standard shellfish surveys [cockle (Cerastoderma edule) and mussel (Mytilus edulis)].

We will give an overview of the development of the Pacific oyster in the Dutch Wadden Sea from the early 1980s to the current situation. For this we will combine data of shellfish surveys with data on expansion of individual oyster reefs and yearly growth of individual oysters of the subtidal and intertidal areas of the Dutch Wadden Sea.

Method

Occurrence and distribution of Pacific oysters in the Dutch Wadden Sea

Until now no survey program exists that focuses on the development of the Pacific oyster (Crassostrea gigas) in the Dutch Wadden Sea. The invasion of the Pacific oyster between 1983 and 2003 is reconstructed from reports of private persons in these years and from forms that were distributed by Internet (www.Waddenzee.nl) in 2002 and 2003. Respondents were asked to give information on the location, density, size, substrate and year of first record.

Since 1990, every year in spring (March, April and May), the Wadden Sea has been sampled for the distribution of cockles (Cerastoderma edule) and mussels (Mytilus edulis). Sampling is done either with an adapted mechanical cockle dredge (0.4 m2, 7 cm depth) or a hand-held sampling device (0.1 m2, 7 cm depth). Samples are sieved with a 5 × 5-mm sieve. All data from these samples are calculated to individuals/m2 and total fresh weight/m² (Goudswaard et al. 2007). Between 1,370 and 1,775 samples were taken each year, following a stratified approach based on occurrence of cockle (Cerastoderma edule) or mussel (Mytilus edulis) beds, not on Pacific oyster occurrence. The position of the samples (1,110–4,440 m apart) was based on data of previous surveys and information from fishermen on the occurrence of cockle or mussel beds.

Intertidal mussel beds are also charted with GPS, and occurrence of Pacific oysters has been documented briefly in recent years.

Development of individual Pacific oyster reefs

Pacific oyster (Crassostrea gigas) development on three locations has been followed since 2002 by establishing the contours of the Pacific oyster reefs (congregated oysters in large densities) and by measuring Pacific oyster lengths (in cm) in samples of 0.5 × 0.5 m².

On a Pacific oyster reef near Texel (location 216) every year three squares at fixed locations are analysed. On a mussel bed (Mytilus edulis) near Ameland (location 502) every year four squares at fixed locations are analysed. On this mussel bed Pacific oysters also occurred in the regular samples taken during mussel surveys. These consisted of five yearly samples of 0.05 m². On the dyke near Oudeschild (Texel) two sample transects (O-10 and O-13) of approximately 6 m length and 0.5 m width were established from high water level to low water level.

The contours of the Pacific oyster reefs are established by satellite positioning (GPS) of the borders of the reef. Only parts that have a cover of more than 5% (Brinkman et al. 2003) are considered part of the reef. Inlets that are bigger than 25 m are recognised, and lumps of Pacific oysters laying 25 m apart from the reef are treated as separate (patch) reefs. The contours of the reefs are analysed with GIS.

Results

Current status of the Pacific oyster in the Dutch Wadden Sea

From data on first occurrences of Pacific oysters (Crassostrea gigas) in the Dutch Wadden Sea, it is clear that a relatively long period elapsed between the first record (1983) and reports on further spread in the area around 1995 (Fig. 1). Only just before the start of the new millennium did it become clear that the Pacific oyster had spread throughout the entire Dutch part of the Wadden Sea.

Fig. 1
figure 1

Reports by private persons on the Interwad website (www.waddenzee.nl) in 2002 and 2003 of Pacific oyster (Crassostrea gigas) occurrences in the Dutch Wadden Sea. Years indicate first reports at specified locations, years between brackets indicate the supposed year of settlement based on the size of the oyster. Coordinate system in Dutch national reference system RD

In the early stages of Pacific oyster development only small areas were covered with Pacific oysters. Information on the substrates where Pacific oysters were found and corresponding abundances came from respondents to the Internet inquiry (Fig. 2). Since 2004, however, the Pacific oyster abundance increased in such a way that it could be estimated through yearly shellfish surveys (Table 1). This survey is based on the occurrence of mussels (Mytilus edulis) and cockles (Cerastoderma edule) and not on the stratification of Pacific oysters. Therefore, the reliability of estimating the standing stock of the Pacific oyster is still far from optimal. Nevertheless, these estimates can be used to assess the increase in Pacific oyster occurrences (Table 1; Fig. 3).

Fig. 2
figure 2

Substrate for settling of Pacific oysters (Crassostrea gigas) in the Dutch Wadden Sea. Based on Internet enquiry in 2002 and 2003

Table 1 Samples containing Pacific oysters (Crassostrea gigas) in yearly shellfish surveys (spring) based on mussel (Mytilus edulis) and cockle (Cerastoderma edule) occurrences
Fig. 3
figure 3

Percentage of samples (grey bars) in the yearly tidal shellfish survey (1,293 samples), based on mussel (Mytilus edulis) and cockle (Cerastoderma edule) occurrences, in which Pacific oyster (Crassostrea gigas) was found and the mean biomass (fresh weight? g/m²) of oysters in these samples (black bars) with standard deviation

Before 2002 Pacific oysters were found in fewer than 4 out of more than 1,500 samples. After 2005 this increased to more than 70. The mean biomass of Pacific oysters within these samples increased from 132 to 2,050 g/m2 fresh weight (Fig. 3). From these results it can be concluded that the Pacific oyster increased explosively between 2003 and 2005, but it seems to have stabilised since then. A comparably fast increase has also been found in the German Wadden Sea (Diederich et al. 2005; Brandt et al. 2008; Schmidt et al. 2008).

Development of individual Pacific oyster reefs

Since 2002 several locations with Pacific oyster development have been followed on a yearly basis. The measurements focus mainly on the length-frequency distribution of Pacific oysters and on the contours of the Pacific oyster reefs.

Pacific oysters at location 216 near Texel (Fig. 4) settled on a dead shell bank/bed (mainly dead Cerastoderma edule) in 1999. In 2002 the Pacific oysters started to take a vertical position, new spat fell onto older individuals (2002 and 2003), and they congregated to large cemented clumps, thus forming a reef. In 2006 this reef was 17.4 ha (Fig. 4) with a mean abundance of 317 oysters/m² on the Pacific oyster-covered parts (Fig. 5a, b). The area of the reef does not seem to have increased much since 2004 (Fig. 4). The reef consists of many different size classes (Fig. 6b).

Fig. 4
figure 4

Development of the contours of the Pacific oyster (Crassostrea gigas) reef (216) near Texel from 2002 to 2006

Fig. 5
figure 5

a Mean abundance of Pacific oysters (Crassostrea gigas) per m² on the oyster reef near Texel (Zeeburg, 216). b Mean abundance per length class (in cm) on the Pacific oyster (Crassostrea gigas) reef near Texel (Zeeburg, 216)

Fig. 6
figure 6

Mussel (Mytilus edulis) bed (502) near Ameland, which has been partly transformed into a Pacific oyster (Crassostrea gigas) reef

Several small trenches and pools occur in the reef, which contain, among others, gobies (Gobio gobio) and shrimp (Crangon crangon). Eider ducks (Somateria mollissima), red knots (Calidris canutus), Eurasian spoonbills (Platalea leucorodia) and little egrets (Egretta garzetta) have been observed at this oyster reef and its interspaces (own observations). On the Pacific oyster reef Gutweed (Enteromorpha sp.), sea letuce (Ulva spec.) and the red algae Gracilaria verrucosa occur. On the shells the periwinkle Littorina littorea is very common.

Near Ameland (location 502, Fig. 6) Pacific oysters have settled on a mussel bed (Mytilus edulis). Part of this mussel bed, which originates from 1994, is already transformed into a Pacific oyster reef. In 2000 the first Pacific oyster was found on this mussel bed. Since 2002 the size classes of these Pacific oysters have been measured, and the progress of replacement of mussels (Mytilus edulis) by Pacific oysters has been monitored. Pacific oysters are now mainly concentrated on the eastern side of the mussel bed and cover an area of about 26.5 ha (Fig. 6) with a mean abundance of 336 Pacific oysters/m2 on the patches (Fig. 7a, b). On the mussel bed both Gracilaria and Fucus sp. occur. The abundance and length frequency distribution data (Fig. 7a, b) show a notable increase in both number and size of Pacific oysters on this location. In the five samples of 0.05 m² taken for the mussel inventories, Pacific oysters were observed from 2003 onwards (Table 2). Variation between samples is such that they are unreliable for biomass estimates for the bed as a whole.

Fig. 7
figure 7

a Mean abundance of Pacific oysters (Crassostrea gigas) per m² on a mussel (Mytilus edulis) bed (502) near Ameland. b Mean abundance per length class (in cm) of Pacific oysters (Crassostrea gigas) on a mussel (Mytilus edulis) bed (502) near Ameland

Table 2 Abundance of Pacific oysters (Crassostrea gigas) in samples from mussel (Mytilus edulis) beds

Near Oudeschild (Texel location O-10 and O-13) Pacific oysters have settled on the stones of the dike. First settlement probably took place around 1999 (Figs. 8, 9). It is clear that abundances are higher on the lower parts of the dike (Fig. 10). In 2003 a notable increase in abundance took place, which is a clear result of the spat fall in 2002 (Fig. 8). The spat fall of 2003, which was clearly visible at several locations in the tidal flats in the Wadden Sea, did not result in an increase in oyster densities on the dike in 2004. On one of the transects (O-10) a notable decrease in oyster densities occurred in 2004. The oyster abundance on the two transects of the dike of Oudeschild (O-10 = 43 inds./m² and O-13 = 38 inds./m2) (Figs. 10, 11) seems to be considerably lower than the density on the tidal flats (Zeeburg-216 = 317 inds./m2 and Ameland-502 = 336 inds./m2).

Fig. 8
figure 8

Abundance of Pacific oysters (Crassostrea gigas) per m² on the dike of Texel (location O-10 and O-13)

Fig. 9
figure 9

Size distribution (in cm) of Pacific oysters (Crassostrea gigas) on the dike of Texel (location O-10 and O-13)

Fig. 10
figure 10

Abundance of Pacific oysters (Crassostrea gigas) on the dike of Texel (location O-10 and O-13) with 11 samples from high water level to low water level

Fig. 11
figure 11

Biomass (g/m² fresh weight) of Pacific oysters (Crassostrea gigas) on an oyster reef (216) (left) and mussel (Mytilus edulis) bed (502) taken over by oysters (right). Notice different scales

Based on abundances and size classes and an equation developed by Büttger and Nehls (personnel communication), the Pacific oyster biomass was calculated for the samples from beds 216 and 502 (Fig. 11). In the period that the oyster biomass (fresh weight) was increasing in Ameland, it stabilised or even decreased at Zeeburg. Still, the total Pacific oyster biomass (fresh weight) at Zeeburg was considerably higher.

Reproduction and individual growth of Pacific oysters

Although no data are available on the quantity of Pacific oyster recruitment in the Dutch Wadden Sea, from the length-frequency distribution data on these different locations it can be deduced that successful spat fall did occur every year between 2002 and 2006.

Although no study has been performed aiming specifically at Pacific oyster growth under special conditions and at different locations in the Wadden Sea, from several studies on length frequencies at least some information on growth can be deduced. Oysters seem to grow up to 4 cm in the first year (this study; Cadée 2000; Reise 1998; Wehrmann et al. 2000; Schmidt et al. 2008). Cadée (2000) also measured growth of second year oysters in the Mokbaai (Texel). He found an increase from 3 to 8 cm between April and August.

In combination with the data on size distribution presented in this article, it can be deduced that oysters grow about 3–4 cm in the first year and about 2–3 cm in the second year (Figs. 6, 9, and 11).

Takeover of mussel beds?

Mussel bed area (Mytilus edulis) in the tidal area has increased from almost zero in the early 1990s to around 2,000 ha in 2006. Figure 12 indicates the area of Pacific oyster reefs and mussel beds with oysters in the period from 2003 to 2006. Although considerable parts of mussel beds have been invaded, there is still a majority of mussel beds in the Dutch Wadden Sea with no or insignificant amounts of oysters.

Fig. 12
figure 12

Area of tidal mussel (Mytilus edulis) beds and Pacific oyster (Crassostrea gigas) coverage in the Dutch Wadden Sea. The oyster coverage is estimated in percentages by eye, and the survey took place during the intertidal mussel bed survey (spring)

Discussion

Since the new millennium the Pacific oyster seems to have established a solid place in the Wadden Sea ecosystem. Pacific oysters occur on hard substrate (dikes and stones), on other secondary hardsubstrates [e.g., mussel beds (Mytilus edulis)] and on shell banks. They can occur as clusters, patches or reefs.

Although recruitment was irregular at first (Dankers et al. 2004; Reise 1998; Wehrmann et al. 2000), wild populations have increased at a dramatic pace since around the year 2000. Several factors might be of importance for this pattern of slow settlement followed by explosive growth. It always was expected that Pacific oysters would be at the edge of their physiological range in the Wadden Sea and would rely on high late summer temperatures, which may occur at least every 5 years (Diederich et al. 2005). However, there are indications that eco-physiological changes might have occurred, which made it possible for the Pacific oyster to reproduce in an unfavourable environment (Cardoso 2004; Wehrmann et al. 2006). In addition, the rise in temperature during the last decade, improving environmental conditions for natural recruitment and dispersal, might explain the recent explosive growth (Diederich et al. 2005; Nehls et al. 2006; Wehrmann et al. 2006). A positive feedback mechanism might also contribute to the pattern. A study by Diederich (2005) and Wehrmann et al. (2006) has revealed that settlement of Pacific oysters is higher on conspecifics than on any other substrate. This might lead to rapid reef formation after first settlement and thus explosive growth after slow establishment.

The absence of predators might enforce a rapid expansion in a new habitat once it has been established (Colautti et al. 2004). In the Dutch Wadden Sea, Pacific oysters do not suffer from high predation levels. Crabs (Carcinus maenas) are able to open oysters of 1–3.5 cm, but prefer mussels (Mascaro and Seed 2000). Herring gulls (Larus argentatus) and oystercatchers (Haematopus ostralegus) are able to predate on Pacific oysters between 3 and 10 cm (Cadée 2001, 2008a, b, Wehrmann, personnel communication). In the Oosterschelde Pacific oysters suffer from suffocation by the colonial tunicate Didemnum lahillei. This species was found for the first time in the Oosterschelde in 1991, and occurrences peaked between 1996 and 1998 (de Kluijver and Dubbeldam 2003). This species does not (yet) occur in the Wadden Sea.

The highest “predation pressure” comes from humans, collecting Pacific oysters for consumption. Pacific oysters are only collected by hand in the Dutch Wadden Sea, and only smaller specimens lying individually are collected. Pacific oysters from reefs are not suitable for collection as they form inseparable clusters of individuals in different sizes and shapes. Until now no pests or diseases have been observed on Pacific oysters in the Wadden Sea (Cadée 2004).

There are, however, several reports of (mass) mortality of Pacific oysters (Cadée 2004; Soletchnik et al. 1999; Dankers et al. 2006; Wehrmann et al. 2006). Suffocation underneath large bundles of green algae (Ulva sp.) (Cadée 2004) or silt does occur. In France summer mortality occurs regularly (Soletchnik et al. 1999). In the summer of 2004 (August), several private persons reported mass mortality (>50%) of oysters in the Dutch Wadden Sea. Especially near Ameland, but also near Texel, Terscheling and the dike of Lauwersoog high percentages of oysters died. There were also locations, however, where mass mortality did not occur (e.g., Balgzand and Eemshaven). In the summer of 2005 mass mortality occurred in several locations of the German Wadden Sea (Marc Herlyn, personnel communication).

Although mortality can reach up to 90% in some locations, so far the cause of mass mortality is not clear. There are no indications, however, that this mortality would bring the invasion of the Pacific oyster to a stop. The shells of oysters remain after mass mortality and form an excellent substrate for new settlement of Pacific oysters.

Possible effects of the increase of Pacific oysters in the Wadden Sea consist of a decrease in space for other organisms, a decrease in carrying capacity and habitat change.

Pacific oysters settle regularly in mussel beds (Mytilus edulis) and compete for space with this species (May 2006). Although Pacific oysters have taken over several mussel beds in the Dutch Wadden Sea, this is still considered only relevant at a small scale (Nehls et al. 2006). There are also several reports of mussel spat settling on Pacific oyster reefs. Thus, although Pacific oysters can overgrow mussel beds, they themselves again form substrates for mussels. This might in the end result in combined reefs.

In the German Wadden Sea almost all mussel beds are now considered oyster reefs (Nehls et al. 2006; Wehrmann et al. 2007). This may be due to the fact that in The Netherlands rejuvenation of mussel beds occurs, while in Germany most mussel beds are senescent.

As Pacific oysters are filter feeders, they may compete for food with other filter feeders, like mussels (Mytilus edulis) and cockles (Cerastoderma edule). In the Oosterschelde there seem to be indications that this is the case (Geurts van Kessel et al. 2003). In the 1990s the Pacific oyster biomass increased significantly in this area, while the mussel and cockle biomass decreased. However, effects of climate change on this change in bivalve species compensation cannot be ruled out. Nehls et al. (2006) consider competition for food and space with the Pacific oyster to be far less important for species displacement than climate change.

Pacific oysters create a hard three-dimensional substrate by forming reefs. Because of this, habitats change from soft to hard substrate. Fauna that prefer a soft substrate, such as cockles (Cerastoderma edule), might suffer from this habitat change. Fauna that prefer a hard substrate can benefit from the occurrence of Pacific oyster reefs. As a result, oyster reefs on tidal flats seem to harbour more species compared to the bare flats (van Broekhoven 2005; Markert 2006). Pacific oysters that attach to hard substrates like dikes and stones also induce an increase in species richness; however, when the covering percentage exceeds 50%, the species richness declines again.

Although part of the effects of Pacific oyster growth might be seen as negative, it should be kept in mind that native oysters (Ostrea edulis) have always occurred in the Wadden Sea, and, although this concerns another species, the introduction of the Pacific oyster might in the end result in a (re-)establishment of reefs as habitat in the Wadden Sea.

As can be deduced from this overview of data on the development of the Pacific oyster in the Dutch Wadden Sea, there are limitations in the accuracy of quantitative data on actual biomass. As the Pacific oyster might cause ecological and economical problems in the Dutch Wadden Sea, it is important for Dutch policy making that quantitative monitoring of the biomass development and distribution of the Pacific oysters takes place.