Spatial and temporal patterns of microphytobenthos communities along the marine-terrestrial boundary in the German Wadden Sea

Yong, Joanne; Moick, Melissa; Dewenter, Jana; Hillebrand, Helmut; Kröncke, Ingrid; Lõhmus, Kertu; Pieck, Daniela; Rohde, Sven; Moorthi, Stefanie

doi:10.3389/fevo.2022.956092

ORIGINAL RESEARCH article

Front. Ecol. Evol., 14 September 2022

Sec. Biogeography and Macroecology

Volume 10 - 2022 | https://doi.org/10.3389/fevo.2022.956092

This article is part of the Research Topic Ecological Dynamics Across the Aquatic-Terrestrial Boundary View all 9 articles

Spatial and temporal patterns of microphytobenthos communities along the marine-terrestrial boundary in the German Wadden Sea

$\nJoanne Yong$ Joanne Yong¹^*

Melissa Moick¹

Jana Dewenter^1,2

Helmut Hillebrand^1,3

Ingrid Kröncke²

Kertu Lõhmus⁴

Daniela Pieck¹

Sven Rohde¹

Stefanie Moorthi¹

¹Institute for Chemistry and Biology of the Marine Environment (ICBM-Terramare), University of Oldenburg, Wilhelmshaven, Germany
²Department of Marine Biology, Senckenberg am Meer, Wilhelmshaven, Germany
³Helmholtz Institute for Functional Marine Biodiversity (HIFMB), Oldenburg, Germany
⁴Institute of Biology and Environmental Sciences (IBU), University of Oldenburg, Oldenburg, Germany

Microphytobenthos (MPBs) are the main primary producers in shallow marine ecosystems, such as the Wadden Sea. We investigated the spatial and temporal dynamics of MPB communities across the marine-terrestrial boundary over three seasons (spring, summer, and fall) on three East Frisian Islands (Norderney, Spiekeroog, and Wangerooge) in the German Wadden Sea. Natural transects were compared with 12 experimental islands (salt marsh vegetated vs. initially bare islands) established on the tidal flats of Spiekeroog for studying dispersal-mediated community assembly. Sediment cores were taken along triplicate transects and on three elevation levels of the experimental islands, corresponding to the pioneer (pio) zone, the lower salt (LS) marsh, and the upper salt (US) marsh. On both the natural transects and the experimental islands, the highest MPB biomass was observed in the pio zone, where vegetation-driven sediment stabilization and high-mud content could have promoted MPB biomass in this marine-terrestrial transition zone. On the experimental islands, MPB biomass and diversity significantly decreased with elevation regardless of the season, indicating that the rarely submerged upper salt marsh level supported minimal MPB growth. The MPB biomass was also higher on initially vegetated than on bare islands, which was the most pronounced on the US level. On the tidal flat transects, the MPB biomass significantly increased with elevation up to the pio zone before decreasing again in the LS marsh. Temperature, sediment water content, and grain size significantly affected transect MPB biomass. MPB diversity, on the other hand, was not related to elevation but was rather determined by temperature, mean grain size, and mud content. Our study suggests that extending MPB studies into the “terrestrial” domain of salt marshes enhances our understanding of the microalgae–plant interaction in this important boundary zone.

Introduction

Tidal flats are soft bottom habitats between the mean high-tide and mean low-tide water lines found worldwide along coastlines (Dyer et al., 2000). They form the transition zone between marine and terrestrial environments, and are highly dynamic landscapes subject to variable environmental factors (e.g., irradiance, temperature, nutrients, salinity, and inundation) driven by tides, as well as by daily and seasonal fluctuations (Miththapala, 2013; Du et al., 2017). The Wadden Sea, a UNESCO Natural World Heritage Site, is the world's largest tidal flat system situated on the backwater side of low-lying barrier islands stretching for 500 km off the coast of the Netherlands to Germany and up to Denmark in the north. In addition to the extensive tidal flats in this region, coastal salt marshes with salt-tolerant plant communities border the tidal flats. European salt marshes follow a clear elevational gradient, with vegetation zonation forming the pioneer (pio), the lower salt (LS) marsh, and the upper salt (US) marsh zones (Balke et al., 2017).

In these soft bottom habitats, benthic microalgae or microphytobentho (MPB) inhabit the upper few millimeters of illuminated sediments. The MPB consists of unicellular photosynthetic eukaryotic algae (such as diatoms, dinoflagellates, euglenids, cryptophytes, chlorophytes, and chrysophytes) and cyanobacteria (MacIntyre et al., 1996; Jesus et al., 2009), often forming visible brown or greenish microbial mats on the surface of tidal flats. The primary production of MPB and biomass can exceed that of phytoplankton in the overlying water column in many shallow aquatic systems (MacIntyre et al., 1996; Underwood, 2001), and contribute up to 50% of the total primary productivity in estuarine environments (Underwood and Kromkamp, 1999). The continuous resuspension of part of the MPB community into the water column leads to transport and dispersal by waves and currents (Baillie and Welsh, 1980; MacIntyre et al., 1996; Mitbavkar and Anil, 2002), and also results in MPB contributing to part of the primary production of the phytoplankton community in the water column (MacIntyre et al., 1996; Lucas et al., 2000). In addition, the MPB play an important role in mediating nutrient fluxes between the sediment and the water column (Cabrita and Brotas, 2000; Dong et al., 2000) and provide an important food source for benthic consumers, commercially important fish and shell fish stocks, as well as migratory bird populations (Middelburg et al., 2000; Hillebrand et al., 2002; Scholz et al., 2014). Furthermore, MPB stabilize sediments through the secretion of extracellular polymeric substances (Paterson, 1989; MacIntyre et al., 1996). Aside from being the dominant primary producers in tidal flats, MPB are also present in salt marsh environments bordering the tidal flats (MacIntyre et al., 1996), being brought in during flood and ebb spring tides (Redzuan and Underwood, 2020).

Overall, a variety of abiotic and biotic factors determine the spatial and temporal dynamics of biomass and composition of MPB, resulting in a patchy distribution on small local scales driven by nutrient distribution (Thornton et al., 2002), light availability (Kühl and Jorgensen, 1994), and grazing by consumers (Kelaher et al., 2003). The North Sea and the back barrier tidal flats of the East Frisian Islands are strongly influenced by a semi-diurnal tidal system with wide tidal ranges greater than 3 m (Dittmann, 2012) inducing highly variable local environmental conditions in the tidal flat systems. However, regionally different hydrodynamic conditions resulting in differential sediment erosion and deposition may lead to large-scale habitat differences among neighboring islands, for instance, in sediment composition, grain size, and slope gradient leading up to the shoreline. These differences, in turn, may also influence MPB distribution and biomass (Chapman et al., 2010).

Algal biomass tends to be higher on the upper tidal flats due to longer emersion duration and shallower water depth during high tide, which results in higher irradiances, warmer temperatures, and higher sediment stability (low resuspension rate) (van der Wal et al., 2010, 2017; Daggers et al., 2020; Jacobs et al., 2021). Algal biomass is further strongly influenced by sediment properties, which affect water content, organic matter content, and light penetration. For instance, muddy sediments were shown to support higher biomass than sandy sediments (Brotas et al., 1995; Herman et al., 2001; Billerbeck et al., 2007). Also, seasonal succession patterns of MPB based on seasonal increases in irradiance and temperature have been demonstrated, with some systems having diatom bloom peaks in spring and/or autumn, while others have summer cyanobacteria bloom peaks (Thornton et al., 2002; Montani et al., 2003; van der Wal et al., 2010; Scholz and Liebezeit, 2012a; Brito et al., 2013). Cyanobacteria dominate MPB communities, especially in summer, when temperatures are high (Watermann et al., 1999), as they usually exhibit lower sensitivities to high temperatures compared to other algal groups like diatoms. Filamentous cyanobacteria especially dominate MPB communities on the upper shore and within salt marshes during the warmer seasons (Kaas, 1987; Tracy and South, 1989), resulting in high algal biomass, but low species richness at higher elevations. High MPB biomass often coincides with low diversity, where usually only one or few species dominated the MPB community (Colijn and Dijkema, 1981; Hillebrand and Sommer, 1997; Forster et al., 2006; Ribeiro et al., 2021).

Complementing the studies described above that investigated the effects of abiotic conditions (temperature, nutrients, and light) and sediment properties (grain size, water content) on MPB communities, numerous studies have been conducted to investigate spatial and temporal dynamics of marine MPB on tidal flats (Facca and Sfriso, 2007; Du et al., 2009; Jesus et al., 2009; Scholz and Liebezeit, 2012a,b; Vogt et al., 2018). However, most of these studies were purely observational and conducted at small local scales, while studies investigating MPB communities over larger regional scales are scarce (Sawai et al., 2016). Additionally, most of these studies only covered the marine gradient from subtidal to intertidal sites. Even in the few studies that have included the supratidal salt marsh (Brotas et al., 1995; Janousek, 2009; Redzuan and Underwood, 2020), MPB communities were compared only between selected tidal flat and salt marsh sites and not along a continuous gradient including the transitional zone across the marine-terrestrial interface.

In the present study, we contribute to the current knowledge mainly in two novel aspects: (1) by extending the analysis into the supratidal zone as we include the vegetated (salt marsh) part of the elevational gradient, and (2) we use the unique opportunity to compare natural transects with data from experimental islands. These twelve experimental islands have been installed in 2014 on the back-barrier tidal flat of Spiekeroog, each consisting of three different elevation levels mirroring the pio, LS, and US marsh, to study spatial processes determining community assembly, species sorting, dispersal, and competition (full details can be found in Balke et al., 2017). Six islands were initially bare and colonized via dispersal from source communities of the salt marsh and tidal flat, and the other six received transplanted LS marsh plant communities. We were, therefore, able to compare MPB communities in initially vegetated vs. unvegetated versions of the pio, LS, and US marsh zones.

We monitored the spatial and temporal distribution of MPB communities on the experimental islands, following a seasonal succession from spring, over summer, to autumn, and compared these patterns to natural transects across the marine-terrestrial interface of three East Frisian Islands (Norderney, Spiekeroog, and Wangerooge) differing in regional hydrodynamics and local environmental factors. We investigated chlorophyll-a (Chl-a) biomass as a proxy for MPB biomass, and also calculated MPB diversity based on pigment diversity, which is to our knowledge, a novel approach for studying MPB diversity in intertidal sediments. Previous studies have mainly used pigment composition data to analyze functional group composition of MPB (Lucas and Holligan, 1999; Méléder et al., 2005; Janousek, 2009).

We formulated our hypotheses as following: (H1) Total algal biomass is negatively correlated to pigment diversity. (H2) Total algal biomass increases with elevation and peaks toward the salt marsh, while pigment diversity decreases with elevation. (H3) Total algal biomass and pigment diversity are constrained by physical factors, such as temperature and sediment properties (water and mud content, mean grain size) that explain patterns in space and time. On the experimental islands, we additionally tested the following hypothesis: (H4) Algal biomass and pigment diversity are higher on the initially vegetated islands compared to the bare islands due to higher sediment water content as a result of vegetation cover.

Materials and methods

Study area

The study area is situated within the Lower Saxony Wadden Sea National Park in Northern Germany and focuses on the back-barrier tidal flats of the East Frisian Islands Norderney, Spiekeroog, and Wangerooge (Figures 1A,B). On each island, three transects (Figures 1A–C) were sampled across the marine-terrestrial interface, starting at the edge of the tidal flats close to the salt marsh and extending from 720 to 1,000 m in seaward direction (Figures 1C–E). Transects were spaced from 300 to 450 m apart, while sampling sites were spaced from 150 to 450 m apart. Site I was the sampling site situated furthest outside on the tidal flat of each island, and the sampling transect extended inwards down to site IV, which is the area right before where the salt marsh boundary zone starts (Figures 1C–E). Transects on Spiekeroog were extended to include the pio zone (site V) and LS marsh (site VI) of the local natural salt marsh (Figure 1C). The US marsh could not be sampled, due to the extremely dry and hard sediment, which is further made impenetrable by a tough layer of plant roots. Natural salt marshes with different vegetation zones were not present on the other two islands, which were thus only sampled up to site IV. The range of elevation (meters above the mean sea level; m AMSL) of the sampling sites differed slightly from one island to the other. Sampling sites on Norderney ranged between −0.4 and 0.4 m AMSL, while the sampled sites on Wangerooge ranged from −0.3 to 0.85 m AMSL. The transects on Spiekeroog covered elevations between 0 m and 0.9 m AMSL for sites I to IV, and up to 1.7 m AMSL for the pio and LS marsh zones. For more precise GPS locations, elevation data, and sampling dates, refer to Tables S1 and S2 in the Supplementary material.

FIGURE 1

Figure 1. (A) Location of the East Frisian Islands on the north-western coast of Germany, in the state of Lower Saxony. (B) Close-up view (marked in red) of the three East Frisian Islands involved in this study (L-R) Norderney, Spiekeroog, and Wangerooge. (C) Close up image of individual sampling sites (I-VI) and transects A to C on Spiekeroog. Included also are the sites of the twelve experimental islands (shown as yellow pins). (D,E) Close up images of sampling sites on Norderney and Wangerooge, respectively. All images are modified Google Earth satellite images. GPS coordinates of individual sampling sites are listed in the Supplementary material.

Experimental islands

Twelve experimental islands made of steel cages measuring 2 m × 6 m were installed in September 2014 on the back-barrier tidal flat of Spiekeroog (for construction design details and experimental aims, refer to Balke et al., 2017). Each island consists of twelve 1 × 1 m metal cages, of which four cages were joined together to form one of the three elevation levels (70 cm, 100 cm, and 130 cm AMSL). These elevation levels correspond to the pio, LS, and US zones of the natural salt marsh on Spiekeroog. Six of these islands were initially non-vegetated (islands with even numbers in Figure 1C) and were filled with sediments from the surrounding tidal flat. The other six islands (islands with odd numbers; Figure 1C) initially received transplantations of LS marsh vegetation sods (to a depth of 30 cm) from the natural lower salt marsh on Spiekeroog (Balke et al., 2017). For the current study, samples were also taken from all the three elevational levels of the 12 experimental islands.

Sampling and sample preservation

Sampling for total algal biomass, pigment composition, and sediment properties (water and mud content, mean grain size) was performed in spring, summer, and autumn of 2019 (for exact dates, see Supplementary Table S1) during low tide. Cut-off syringes (ϕ 26 mm) were used to extract sediment cores from the tidal flats and the different elevation levels of the experimental islands. Sediment core samples were taken in triplicates per site, and the upper 2 cm of the three extracted cores, containing the majority of phototrophic cells (De Jonge and Colijn, 1994), were pooled and transferred into labeled plastic zipper bags, resulting in a total sediment volume of 31.85 cm³ per site. When sampling within the salt marsh or on experimental islands, extra care was taken to exclude any visible pieces of plant material or detritus from the sampled sediment cores. Sediment samples in the field were kept in the dark and cooled with ice packs in a cooler box. Upon return to the laboratory, sediment samples for pigment extraction were immediately transferred to −80°C, while the samples for other parameters were stored at −20°C.

Temperature

With the exception of spring temperature data from Spiekeroog, all in-situ sediment temperatures at all sampling sites were measured with a hand-held sensor (Fibox 4; PreSens Precision Sensing GmbH, Regensburg, Germany). Due to a technical malfunction with our hand-held sensor during the spring sampling campaign on Spiekeroog, in-situ sediment temperatures were recorded using four HOBO^® temperature/light loggers (HOBO^® UA-002-64 pendant, Onset^® Corporation, Bourne, MA, USA) instead. The loggers were deployed at 2 cm sediment depth, anchored securely into the tidal flat sediment at sites I to IV of transect B, and left in the field to capture light and temperature fluctuations of at least one semidiurnal tidal cycle before retrieval (24 h). Temperature data from the pio zone (site V) and LS marsh (site VI) on Spiekeroog in spring were derived from a long-term monitoring program with permanently installed temperature loggers (Meier et al., 2020). The data were recorded with DEFI-T loggers (JFE Advantech Co., Ltd., Tokyo, Japan), which were installed at 5 cm depth in the salt marsh soil. On the experimental islands, no temperature data could be recorded in spring due to the technical malfunction.

Sediment characteristics

Sediment water content was analyzed in triplicates and determined by drying 5 g of wet sediment at 60°C for 48 h. Water content was defined as the weight ratio before and after drying. Average grain size and mud content was determined using a laser diffraction particle size analyzer (Horiba LA-950; Horiba Europe GmbH, Oberursel, Germany). Prior to analysis, the samples were suspended in a sodium metaphosphate solution (4% w/v) and coarse particles (>2 mm) were removed by sieving. Samples from the salt marsh and experimental islands were additionally treated with sodium peroxide (30% w/v) and heated to remove all organic plant material. The same treatment was not required for samples from the unvegetated tidal flat. On the experimental islands, sediment water content was the only parameter analyzed for all three elevation zones of the 12 islands.

Sediment pigment analysis

Pigment extraction and analysis were performed according to the method by Thrane et al. (2015) with slight modifications as described here: sediment samples were freeze-dried for at least 48 h. Dried samples were homogenized with a spatula and subsamples (~5 g) were directly weighed into pre-tared 15 ml centrifuge tubes. 10 ml of analytical grade 96% ethanol (Art No.: 147194.1212, PanReac AppliChem, Darmstadt, Germany) was added to each sample, followed by mixing samples thoroughly with a vortex mixer, and lastly, ultrasonication on ice for 30 min. Samples were stored at 4°C in the dark for 20–24 h, before centrifugation at 4,000 rpm, 4°C for 30 min. The supernatants (pigment extracts) were measured in triplicates of 330 μl each. Each sample replicate was pipetted into individual wells of a 96-well flat bottom clear polypropylene plate (Greiner bio one #655201). Pigment absorbances were measured using a Synergy MX plate reader (BioTek instruments, Vermont, USA). Spectral scans were performed between 300 and 800 nm with a resolution of 1 nm to determine Chl-a concentration as a proxy for total algal biomass and pigment composition to calculate pigment diversity as a proxy for MPB taxa diversity (see Section Data analysis).

Data analysis

Spectral scans from the plate reader were analyzed in R (R Core Team, 2021) using the script provided by Thrane et al. (2015), based on a modified Gauss-peak spectra (GPS) method that represents individual pigment spectra as weighted sums of Gaussian functions. The method detects biomass concentration of the following pigments (in microgram per gram dry sediment): alloxanthin, beta-carotene, 9′-cis-neoxanthin, trans-canthaxanthin, chlorophyll-a, -b, -c1, and -c2, trans-diadinoxanthin, diatoxanthin, dinoxanthin, trans-echinenone, fucoxanthin, lutein, myxoxanthophyll, peridinin, pheophytin a and b, and violaxanthin. Several of these are commonly used as diagnostic pigments for certain algal taxonomic groups (Supplementary Table S4). Pigment concentrations were then used to calculate the values of effective number of species (ENS; based on pigment and not species diversity), using the R package vegan (Oksanen et al., 2013), where ENS is derived as an inverse Simpson diversity index (1/D). In our case, D = ∑(n/N)², where “n” represents the concentration of a particular pigment, and “N” represents the total concentration of all pigments measured in a sample. ENS is generally defined as the number of species that would be observed if all of the species in a sample were equally abundant (Jost, 2006). In other words, it is the minimum number of species that is required to achieve the same entropy (diversity) as the observed one. We provide a simple illustrative example here: if pigment A comprises 99% of the concentration and 1% is contributed by pigment B, then pigment richness would be 2, but pigment ENS would be barely larger than 1 as effectively only one pigment dominates. If pigments A and B each comprises 50% of the concentration, pigment ENS would be 2 and the richness is equal. The ENS is a robust diversity measure with regards to sampling and dominance issues, as it weighs elements by their frequency without disproportionately favoring either rare or common elements (Chase and Knight, 2013).

We chose to use the modified GPS method by Thrane et al. (2015) over other methods, because it is a fast, inexpensive, and relatively simple to operate high throughput method for analyzing algal pigment mixtures. Additionally, all aspects from model fitting to data output works well with open-source software (e.g., R). Since the method has already been successfully tested on lake phytoplankton and lake sediments in the original study, we adopted it for analyzing MPB communities in intertidal sediments.

Prior to statistical analyses, Chl-a data was logarithmically transformed and tested with a Shapiro–Wilks test to ensure normality. The hypothesis, H1 was tested using a non-parametric Spearman rank order correlation for algal biomass and pigment diversity. Hypotheses H2–H4 were tested with linear mixed-effects models (LMMs) on predictor variables, using the lmer() function from the R package lme4 (Bates et al., 2014). For H2, the effect of elevation was tested against Chl-a biomass or ENS. The model for the natural transects data set included the transect unit (A–C) nested into island location (Norderney, Spiekeroog, and Wangerooge), and season as random effects. For the experimental islands, H2 and H4 were tested with a single model that included island type (bare or vegetated) and elevation as interacting fixed effects, whereas island number and season were included as random effects. H3 was tested in a separate LMM for Chl-a biomass and ENS, which included temperature and three sediment parameters (mean grain size, sediment water, and mud content) as fixed effects. The random effects were determined as described above. Since the model output for ENS generated a singular fit caused by the lack of spatial variables-associated variance, we simplified the random effects structure down to season alone. Due to the fact that we were unable to obtain all environmental parameters from the experimental islands, the LMM for testing H3 on the experimental islands only included temperature and sediment water content as fixed effects. Additionally, since spring temperature data for the experimental islands was unavailable, only data from two seasons (summer, autumn) was used in the statistical analysis for H3. All statistical analyses were run separately for the natural transects and experimental islands data sets. Marginal and conditional R² based on Johnson (2014) were calculated to estimate the respective model fit.

Results

In-situ environmental parameters

In-situ sediment temperatures on the tidal flats ranged between 8.8°C (in spring) and 27°C (in summer) and were within a similar range in spring and autumn. Spring temperatures in the salt marsh (sites V and VI) on Spiekeroog were considerably lower than those on the tidal flats, resulting in a high temperature variability across the marine-terrestrial gradient (Supplementary Figure S1). On the experimental islands, sediment temperatures ranged between 11.6°C (in autumn) and 22.7°C (in summer) and were on average higher and more variable on bare islands (Supplementary Figure S2).

The tidal flat on Norderney is characterized by muddy sediments, while the one on Wangerooge is predominantly sandy. On Spiekeroog, the tidal flat is characterized by a mixture of muddy and sandy zones, with sites in the middle of the transect being predominantly sandy, while sites on both the ends of the transect are predominantly muddy (Supplementary Table S3).

Overall, the mean sediment grain size was negatively correlated to increasing elevation (r = −0.26, p = 0.0061, Supplementary Figure S3A). This pattern was mainly driven by the Spiekeroog transect, which had the longest elevational gradient where the mean grain size decreased toward the salt marsh (Supplementary Figure S3B, r = −0.68, p < 0.001), and was much finer within the salt marsh than in the marine tidal flat. On the other two islands, the correlation coefficient between the mean grain size and elevation was weaker on Wangerooge (r = −0.35, p = 0.038) and positively correlated on Norderney (Supplementary Figure S3B, r = 0.37, p = 0.028).

The mean grain size and sediment water content of both the sediments exhibited the highest variation on Spiekeroog and lower variation on the other two islands (Supplementary Figures S4, S5). Results showed that the two parameters were inversely correlated to one another. Across transects, the mud content was not significantly correlated to increasing elevation (r = 0.16, p = 0.1, Supplementary Figure S6A). Individual correlation tests performed for each island indicated that the mud content was positively correlated with elevation only on Spiekeroog (r = 0.71, p < 0.001), mainly driven by the sites sampled in the salt marsh. The other two islands revealed negative and weaker correlations between elevation and mud content (Supplementary Figure S6B, r = −0.35, p = 0.047 for Norderney, r = −0.45, p = 0.0061 for Wangerooge).

Relationship between total algal biomass and pigment diversity (H1)

Algal biomass and ENS showed opposing correlations on the transects and experimental islands, which were additionally marginally non-significant. On the tidal flats, the ENS non-significantly decreased with increasing biomass (Figure 2A, r = −0.17, p = 0.082). The lowest ENS values were observed in summer especially on Wangerooge, and the highest in autumn. On the experimental islands, algal biomass and ENS were non-significantly positively correlated (Figure 2B, r = 0.17, p = 0.084).

FIGURE 2

Figure 2. Correlation test results between logarithmic-transformed Chl-a biomass in micrograms per gram dry weight sediment (μg/g dw) and pigment diversity [effective number of species (ENS)] on natural tidal flat transects (A) and on experimental islands (B). Shapes represent seasons and colored symbols represent island location or island type (bare or vegetated experimental islands).

Effect of elevation on total algal biomass and pigment diversity on the natural transects (H2)

Overall, the MPB biomass significantly increased with the increasing elevation (Figure 3A, Table 1), but elevation alone explained only 9.2% of the variance in the biomass, whereas a substantial variation was captured by seasonal and spatial random effects. The elevational effect was mainly driven by high Chl-a concentrations up to 1.25 m on Spiekeroog, but we also observed a decline at the highest elevation (Figure 3A) which matches the decline at the corresponding elevation on the experimental islands (Figure 3B, see below for details).

FIGURE 3

Figure 3. Relationship between logarithmic-transformed Chl-a biomass in micrograms per gram dry weight sediment (μg/g dw) and pigment diversity (ENS) against elevation (meters above mean sea level; m AMSL) on natural tidal flats [left panel; (A,C)] and on experimental islands (bare and vegetated) [right panel; (B,D)].

TABLE 1

Table 1. Linear mixed-effects model (LMM) output from the analysis of the response variables.

We found no significant effect of elevation on pigment diversity, but strong variation between seasons (Table 1, Figure 3C). Seasons stood as the reason for the majority of the random effects on ENS, which was the lowest in summer and the highest in autumn.

Relationship of algal biomass and pigment diversity to temperature and sediment characteristics on the natural transects (H3)

Sediment Chl-a significantly increased with temperature (Figure 4A, p = 0.011) and sediment water content (Figure 4E, p = 0.017). The mean grain size had a negative, but only marginally significant effect (Figure 4C, p = 0.07), while the mud content had a positive but non-significant effect (Figure 4G, p = 0.808) on sediment Chl-a biomass. These fixed effects explained 48.8% of the variance in biomass (Table 2). ENS, on the other hand, significantly decreased with temperature (Figure 4B, p < 0.001) and the mean grain size (Figure 4D, p = 0.019), but increased with the sediment mud content (Figure 4H, p = 0.003). The increase in ENS with increasing water content was marginally non-significant (Figure 4F, p = 0.057). These fixed effects accounted for 27.5% of the variance in ENS (Table 2).

FIGURE 4

Figure 4. Relationship between Chl-a biomass in micrograms per gram dry weight sediment (μg /g dw) (left panel) and pigment diversity (ENS, right panel) with physical abiotic parameters: temperature (A,B), sediment mean grain size (C,D), sediment water content (E,F) and sediment mud content (G,H) on natural tidal flat transects.

TABLE 2

Table 2. Linear mixed effects model output from the analysis of the response variables chlorophyll-a and pigment diversity (ENS) against the predictor variables: temperature, mean grain size, water content, and mud content on natural transects.

Effect of elevation and vegetation on the total algal biomass and pigment diversity on experimental islands (H2, H4)

On the experimental islands, sediment Chl-a biomass was more variable than on the tidal flats and was significantly affected by elevation, island type (vegetated vs. bare), and an interaction of both (p < 0.001, Table 1, marginal R² = 57.5 %). Biomass significantly decreased with the increasing elevation, being the highest in the pioneer zone and the lowest in the US level (Figure 3B). Comparison of the values between transects (with a maximum elevation of 1.69 m AMSL) and the islands (with a minimum elevation of 1.38 m AMSL) revealed similar total biomass for this section of the elevation gradient present in both data sets (Figures 3A,B). Additionally, the experimental islands continued the negative trend of biomass with elevation that appeared at the upward end of the transects as well. The significant interaction between island type and elevation reflected that the decline in the elevation was much stronger for initially bare islands than the vegetated, which led to overall lower biomass in the US zone of bare islands than the vegetated ones (Figure 3B).

Similarly, ENS on the islands was also in the same range as on the transects, and significantly decreased with the elevation on the experimental islands (Figure 3D, p = 0.002, Table 1). However, this pattern was not affected by interaction with the island type, which was not significantly different (Table 1, marginal R² = 16.8%).

Relationship of algal biomass and pigment diversity to temperature and sediment water content on experimental islands (H3)

On the experimental islands, Chl-a biomass significantly increased with sediment water content (Figure 5C, p = <0.001) while temperature did not have a significant effect (Figure 5A, p = 0.366). The fixed effects accounted for 39.4% of the variance in biomass and together with the random effect season accounted for 60.2% of the variance (Table 3). Though non-significant, ENS tended to decrease with the increasing temperature (Figure 5B) but increased with the increasing sediment water content (Figure 5D). However, the random effect season contributed substantially to the conditional R² covering 32.5% of the variance in ENS (Table 3).

FIGURE 5

Figure 5. Relationship between Chl-a biomass in micrograms per gram dry weight sediment (μg /g dw) (left panel) and pigment diversity (ENS, right panel) with physical abiotic parameters: temperature (A,B) and sediment water content (C,D) on experimental islands.

TABLE 3

Table 3. Linear mixed effects model output from the analysis of the response variables chlorophyll-a and pigment diversity (ENS) against the predictor variables: temperature and sediment water content on experimental islands.

Pigment composition on natural transects and experimental islands

On the natural transects, the most commonly detected pigments were alloxanthin, beta-carotene, Chl-b, canthaxanthin, diatoxanthin, echinenone, fucoxanthin, peridinin, and the chlorophyll degradation products, pheophytin a and b (Figure 6). Pheophytin a and b were ubiquitous and abundant during all three seasons and at all sites, as were Chl-b, diatoxanthin, and echinenone. The algal biomarker pigments, such as alloxanthin, fucoxanthin, and peridinin displayed seasonal distribution patterns. Alloxanthin was detected in summer and autumn on Spiekeroog and Wangerooge. Fucoxanthin was present in varying percentages on all islands and sites in summer and was only present at certain sites during other seasons. Peridinin was present on Wangerooge at sites II to IV in summer, and at site IV in autumn. A small percentage was only present at site IV in summer on Spiekeroog. Canthaxanthin exhibited a spatial distribution pattern and was found only in the salt marsh (sites V and VI) on Spiekeroog, as well as at site IV in summer. Additionally, low concentrations of canthaxanthin were also detected at sites III and IV on Wangerooge in autumn.

FIGURE 6

Figure 6. Seasonal pigment:Chl-a biomass ratios (expressed as percentage) calculated for 18 pigments (based on Thrane et al., 2015) from samples sampled from Norderney, Spiekeroog, and Wangerooge (Pigment abbreviations: Allo, Alloxanthin; bb.Car, beta-carotene; c.Neo, 9'-cis-neoxanthin; Cantha, trans-canthaxanthin; Chl.b, chlorophyll-b; Chl.c1, chlorophyll-c1; Chl.c2, chlorophyll-c2; Diadino, trans-diadinoanthin; Diato, Diatoxanthin; Dino, Dinoxanthin; Echin, trans-echinenone; Fuco, Fucoxanthin; Lut, lutein; Myxo, myxoxanthophyll; Peri, Peridinin; Phe.a, pheophytin a; Phe.b, pheophytin b; Viola, violaxanthin).

On the experimental islands, pigment richness was lower than on the natural transects. The chlorophyll degradation products, pheophytin a and b were also ubiquitous across all seasons and elevation levels, accounting for up to 60% of the total pigment biomass (Figure 7). Similar to the transects, Chl-b, diatoxanthin, and echinenone were also ubiquitous. Alloxanthin and peridinin were not detected on the experimental islands. Fucoxanthin represented only a very small percentage of the detected algal pigments, mainly in the LS and US zones in spring and autumn. Canthaxanthin was present in all the three seasons mainly in the US zone of vegetated islands, and in small percentages in the pio and LS zones of vegetated islands in autumn.

FIGURE 7

Figure 7. Seasonal pigment:Chl-a ratios (expressed as percentage) on all three elevational levels: pioneer (pio), lower salt marsh (LS), and upper salt marsh (US) of the experimental islands of Spiekeroog (Pigment abbreviations: Allo, Alloxanthin; bb.Car, beta-carotene; c.Neo, 9'-cis-neoxanthin; Cantha, trans-canthaxanthin; Chl.b, chlorophyll-b; Chl.c1, chlorophyll-c1; Chl.c2, chlorophyll-c2; Diadino, trans-diadinoanthin; Diato, Diatoxanthin; Dino, Dinoxanthin; Echin, trans-echinenone; Fuco, Fucoxanthin; Lut, lutein; Myxo, myxoxanthophyll; Peri, Peridinin; Phe.a, pheophytin a; Phe.b, pheophytin b; Viola, violaxanthin).

Discussion

The current study found a weak inverse relationship between the algal biomass and ENS on the natural transects, and an opposing weak trend on the experimental islands. The opposing signs of marginally non-significant trends falsifies hypothesis H1. We found significant effects of elevation on MPB biomass on both the natural transects and experimental islands. The correlations were different (positive on transects, negative on islands, and partly refuting hypothesis H2), but since both covered different parts of the elevational gradient, they jointly converged with an increasing biomass trend up to the pio zone, and then declined toward the upper salt marsh zones. ENS, on the other hand, did not vary with the elevation on the natural transects, but decreased significantly with elevation on the experimental islands (partly refuting hypothesis H2). Environmental factors affected Chl-a biomass and ENS on the natural transects (accepting hypothesis H3). Biomass was higher at higher temperature and in water content and lower grain size (marginally non-significant).ENS significantly decreased with increasing temperature and mean grain size but increased with increasing mud and water content (marginally non-significant). On the experimental islands, vegetation history played a role in biomass (partly supporting hypothesis H4) but not in the ENS (partly refuting hypothesis H4). Chl-a biomass declined with elevation more on bare islands than on initially vegetated ones, leading to lower overall biomass on the former in the US zone (partly supporting hypothesis H4). Pigment composition results showed some spatial and temporal distribution of certain algal pigments on the natural transects, while both spatial and seasonal variations were minimal on the experimental islands.

MPB biomass and pigment diversity patterns

While the inverse relationship of MPB Chl-a biomass and pigment (community) diversity was weak in our study, the general pattern has been observed before in other microphytobenthos communities (Hillebrand and Sommer, 1997; Forster et al., 2006; Ribeiro et al., 2021). It should be noted, though, that the diversity estimates in the cited studies were based mainly on microscopic analyses of diatom assemblages in contrast to the current study using the overall pigment diversity and encompassing all MPB taxa, but at a much lower resolution. Both Colijn and Dijkema (1981) and Forster et al. (2006) demonstrated that high MPB biomass at certain sites was usually associated with single or few dominant (>50% of population) species. However, in our study this relationship was not significant and inconsistent, especially in the pio and LS zones on Spiekeroog, where high Chl-a biomass and high ENS were observed. MPB communities in this transitional zone across the marine-terrestrial interface are more protected from both marine stressors on the tidal flats (inundation, high irradiance, temperature stress, resuspension, turbidity, etc.) and from terrestrial stressors in the salt marsh (mainly from desiccation and high salinity stress) (Pomeroy, 1959; Janousek, 2009). Hence, such transition zones are considered hotspots for MPB biomass and diversity as they combine high nutrient/water retention properties and high organic matter content (Traut, 2005; Janousek, 2009; Ferro and Morrone, 2014). On the experimental islands, terrestrial stressors increase at higher elevational levels, causing a deterioration in conditions for MPB growth and allowing for only highly adapted species to thrive in this environment. We acknowledge that the use of Chl-a as a proxy for algal biomass in sediments has its pitfalls (Baulch et al., 2009; Ramaraj et al., 2013). Sediment Chl-a may not be completely of in-situ MPB origin but it can be derived from allochthonous sources, such as macroalgae, detritus from macrophytes, sedimenting tycho-/phytoplankton, and terrestrial inputs (Lucas and Holligan, 1999; Steinman et al., 2017; Daggers et al., 2020). This can lead to an overestimation of MPB biomass, as might have been the case in the salt marsh samples due to plant debris. Furthermore, degradation products of Chl-a, such as pheophytin a and pheophorbide a can also absorb light under very similar wavelengths as Chl-a (Steinman et al., 2017), and contribute up to 60% of the measured Chl-a content in freshwater samples (Marker, 1980). Additionally, the amount of Chl-a within cells may vary from species to species, and this is further affected by cell physiology, temperature, light, and nutrient supply, as well as other changes in environmental conditions (Schlüter et al., 2006; Barlow et al., 2008; Steinman et al., 2017). The same applies for other accessory pigments, whose ratios within cells may change not only with changing environmental conditions, but also with senescence and grazing (Wright and Jeffrey, 2006; Chao et al., 2016; Lima et al., 2019), leading to changes in pigment composition and diversity without actual changes in cell abundances within the community. Despite these intricacies, Chl-a remains a simple biomass proxy and pigment composition has successfully been used as a proxy for MPB community structure (Lucas and Holligan, 1999; Méléder et al., 2005; Steinman et al., 2017) and diversity (Janousek, 2009).

Currently, three methods are commonly used for determining algal pigment concentrations: spectrophotometry, fluorometry, and high performance liquid chromatography (HPLC) (Steinman et al., 2017). The latter is considered the “golden standard” for determining pigment concentrations at a high resolution in plants and algae, resolving most chlorophylls, carotenoids, and their degradation products. However, this method is costly in terms of time and the equipment required (Mantoura and Llewellyn, 1983; Wright et al., 1991). The original GPS method by Küpper et al. (2000) is conceptually sound, but relies on the use of single-cuvette spectrophotometry and a separate compilation of data and pigment concentration determination, making the method cumbersome and time-consuming (Thrane et al., 2015). The modified GPS method of Thrane et al. (2015) is a faster and relatively simple to operate high throughput method that has been shown to accurately estimate concentrations of total chlorophylls and total carotenoids in both plankton and sediment lake samples. However, individual pigment concentrations (especially carotenoids) can be difficult to resolve due to similarity between their absorbance spectra, also known as aliasing. To alleviate this problem, the authors assessed aliasing by Monte-Carlo simulation of the identification error rate of each pigment and used the information to come up with a “core set” of correctly identifiable pigments that are common in natural samples and less prone to aliasing. However, the potential for aliasing between certain pigments still exists when using any GPS method, so that the interpretation of the estimated pigment concentrations should be regarded with care, also in our study.

HPLC-based methods, such as CHEMTAX, which enable the calculation of the contribution from different taxa to total Chl-a based on pigment:Chl-a ratios, are more comprehensive and sensitive compared to the modified GPS method used in our study. However, the latter is still able to accurately estimate the total Chi-a concentration in samples and the concentrations of eighteen other pigments (with lower accuracy). Pigment:Chl-a mass ratios could still be generated, which we used to gain an insight into the MPB community structure. We understand that the absorption spectra of some phytopigments may overlap and many pigments are not unique to just a single taxon (Wright and Jeffrey, 2006). Therefore, MPB community composition was not a main focus in our study and pigment-based MPB diversity was likely underestimated in our samples. Another aspect to consider is the fact that a large portion (sometimes >50%) of the pigments measured consisted of chlorophyll degradation products (pheophytin a and b) that were also included into the calculations of ENS. However, these two pigments were ubiquitous in all our samples, increasing the pigment diversity uniformly in the entire data set. Despite its pitfalls, we believe that the modified GPS method was suitable to detect the spatial and seasonal diversity changes in major taxonomic groups within MPB communities. However, we would like to stress that a deeper insight into MPB community structure at a higher taxonomic resolution requires other methods, such as microscopic examinations and molecular barcoding to complement any chemotaxonomic method (Paerl et al., 2003; Wright and Jeffrey, 2006).

The use of pigment concentrations for estimating pigment diversity and relating it to the diversity of taxa within a community is an idea adopted from the field of ecology where ENS is used as a diversity measure. ENS is regarded a robust measure for comparing within-community and among-community diversity, as it overcomes several limitations in terms of both diversity indices and species richness (S) (Jost, 2006; Cao and Hawkins, 2019). The use of ENS as a diversity measure allows us to substitute species abundances with concentrations of detected pigments, and the sum of all species abundances for total pigment concentrations. This method seems relevant in all cases where a microscopic analysis of the MPB communities is not included into the scope of the study. To our knowledge, only one study was based on diversity metrics on pigments and the inverse Simpson index (Janousek, 2009), but with lower number of pigments.

Elevation, MPB biomass, and pigment diversity patterns on natural transects

The significant increase in MPB biomass with elevation in the tidal flats matched the findings of previous studies of other European tidal flats (Brotas et al., 1995; van der Wal et al., 2010; Daggers et al., 2020). Upper tidal flats fall dry for a longer period of time during low tide, resulting in higher stability (low resuspension rate) of the sediment, and longer exposure to sunlight, which therefore increases the photosynthetic period (van der Wal et al., 2017; Jacobs et al., 2021). We did not find a significant relationship between ENS and elevation, which is in contrast to previous studies that observed a decrease in MPB diversity with increasing elevation, driven by local in-situ factors, such as sediment grain size and mud content (Scholz and Liebezeit, 2012b; Ribeiro et al., 2021). Other studies related the opposite pattern (increasing diversity with elevation) to stronger wave action and hydrodynamic conditions at the seaward side of the elevation gradient, potentially causing resuspension of MPB, and the mixing of sediments and thus burying MPB cells into deeper sediment layers (Underwood and Kromkamp, 1999; van der Wal et al., 2010; Jacobs et al., 2021). Janousek (2009) also reported higher pigment diversity in the salt marsh compared to the adjacent bare tidal flats, due to the various conditions provided by salt marsh vegetation already described earlier. These studies and our results indicate that elevation patterns across marine-terrestrial transects cannot straightforwardly be generalized for different systems, as these gradients are strongly associated with in-situ environmental and hydrodynamic factors. Therefore, elevation itself is still likely to be a strong determinant in many systems (including our study system) because hydrodynamics and the gradients of many in-situ environmental factors, such as sediment grain size (and in turn water content), organic matter content (not assessed in this study), and temperature, are closely coupled to the elevational gradient.

Environmental factors, MPB biomass, and pigment diversity patterns on natural transects

Positive relationships between MPB biomass with temperature, and water content have frequently been observed beyond our study system (Colijn and Dijkema, 1981; Lucas and Holligan, 1999; Scholz and Liebezeit, 2012b; Ribeiro et al., 2021). Seasonal variability in MPB biomass and productivity has been related to changes in irradiance and temperature, nutrient availability, and sediment composition in other systems (Colijn and Van Buurt, 1975; De Jonge and Colijn, 1994; Montani et al., 2003; van der Wal et al., 2010), indicating that the determinants for MPB biomass are highly system-specific. Since some of the sediment characteristics vary in a multicollinear way, we suggest some caution in interpreting the LMM results (Gareth et al., 2013). Comparing the LMM results (Table 2) to the bivariate patterns (Figure 4) reveals that the apparent positive association between biomass and mud content did not appear in the LMM, potentially reflecting a similar relationship for water content that is correlated with mud content.

A combination of several forces, such as tidal currents, waves and associated winds, and drainage influences the hydrodynamics on intertidal flats (Le Hir et al., 2000). Hydrodynamics can influence sediment erosion and deposition processes, the distribution of sediment grain sizes and types, and organic matter accumulation. Additionally, factors, such as sediment input and bathymetry can also impact sedimentation rates (Morris et al., 2002; Fagherazzi et al., 2014). In general, sediment grain size decreases with increasing elevation due to reduced wave action, which allows for fine particles like silt and clay to settle close to the salt marshes (Brotas et al., 1995; van der Wal et al., 2017; Jacobs et al., 2021). Finer muddier sediments are associated with higher organic matter, water, and nutrient content, which support higher MPB biomass (Pomeroy, 1959; Brotas et al., 1995; Barranguet et al., 1998; Thornton et al., 2002). Coarser sediments, in contrast, tend to be deposited in areas characterized by higher tidal and wave action. The stronger hydrodynamic forces result in higher instability, export and resuspension of sediment and MPB material, and lower organic matter content (Barranguet et al., 1998). Consequently, Chl-a biomass is often lower at high energy sites.

Tidal currents resuspend and transport sediments to the salt marsh during high slack water, where fine cohesive particles get deposited due to the strong friction of the salt marsh vegetation canopy (Zhang et al., 2019). Salt marsh plants maintain the salt marsh elevation not only by trapping sediment during periods of inundation, but also by building biomass and accumulating organic matter, constantly modifying the elevation in response to the mean sea level (Morris et al., 2002). Climate change induced increases in flooding and storm events will likely lead to longer inundation duration in salt marshes, resulting in either positive or negative responses in vegetation biomass and organic matter accretion, depending on marsh height and associated sediment loads. Hence, future salt marshes will either be able to keep up with the sea level rise, or ultimately drown and be lost to the sea (Kirwan and Guntenspergen, 2012). Future changes in hydrodynamics, potentially resulting in altered elevation gradients and salt marsh heights will also likely impact MPB community dynamics in salt marsh-mudflat systems.

In addition to the conditions associated with finer muddier sediments, the presence of plants in the salt marshes also creates a favorable environment for MPB growth. Sediment is stabilized by the root systems of these plants and is therefore not as erosion-prone as on the open tidal flats (Chen et al., 2012). This positive relationship between MPB biomass and vegetated habitats in contrast to unvegetated tidal flats has also been found before (Pomeroy, 1959; Facca and Sfriso, 2007; Janousek, 2009). Whitcraft and Levin (2007) showed that experimental plots with salt marsh vegetation exhibited lower sediment temperature, as well as higher water content due to shading. Since the US marsh is rarely inundated, and if it occurs, only briefly during the highest spring tides or storm events; the US marsh sediments can become very dry and subject to high salinity stress. Thus, the US marsh becomes unsuitable for algal growth, especially during the hotter and drier summer months (Zedler, 1980).

Effective number of species was also significantly influenced by temperature and sediment characteristics, i.e., mean grain size and mud content. The negative relationship with temperature supported previous findings that show that warmer temperatures often favor dominance of certain MPB groups (Sin et al., 2009; Cartaxana et al., 2015; Ribeiro et al., 2021), whereas diatoms thrive only at lower temperatures (Watermann et al., 1999; Underwood, 2001; Scholz and Liebezeit, 2012b). Janousek (2009) and Semcheski et al. (2016) reported higher MPB diversity in muddy sediments than in sandy sediments, a result of the depository nature of more protected muddy sites, and higher resistance of MPB biofilms on muddy sediments against erosion (Cahoon et al., 1999; Herman et al., 2001). Likewise, Scholz and Liebezeit (2012b) demonstrated that diversity decreases in coarser sediments, while other studies found higher diversity in sandbars than in mudflats from mixed assemblages (Paterson and Hagerthey, 2001; Underwood and Barnett, 2006; Jesus et al., 2009; Rivera-Garcia et al., 2018). These studies indicate that diversity patterns in relation to sediment types may differ, depending on how diversity was defined in the respective study (species or pigments), what communities were targeted (diatoms vs. mixed communities), and how samples were analyzed (chemotaxonomy, microscopy, and counting methods). Janousek (2009) dealt with pigment composition, whilst Semcheski et al. (2016) used the Utermöhl method on sediment samples and Scholz and Liebezeit (2012b) focused on the non-diatom MPB fractions. Sediment type and grain size can also greatly affect light availability for MPB cells within sediments, as well as porewater nutrient content and water retention capabilities during low tide, which in turn determine the MPB composition and diversity (Barranguet et al., 1998; Watermann et al., 1999; Billerbeck et al., 2007).

Aside from abiotic conditions, benthic herbivore grazing may largely control MPB biomass and diversity in tidal flat habitats, especially in summer (Admiraal and Peletier, 1980; Colijn and Dijkema, 1981; MacIntyre et al., 1996; Daggers et al., 2020). Grazing pressure resulted in higher MPB diversity in some studies (Sommer, 1999; Hagerthey et al., 2002; Balvanera et al., 2006) and was shown to be higher in sandy sediments than in muddier sediments (Underwood and Kromkamp, 1999; Herman et al., 2001), indicating the interactive effects of sediment structure and grazing pressure on MPB.

Environmental factors, MPB biomass, and pigment diversity patterns on the experimental islands

Both MPB biomass and ENS significantly decreased with increasing elevation on the experimental islands. This decline in MPB biomass from the pio zone to the US marsh has been well-documented in other studies (Pomeroy, 1959; Zedler, 1980; Thornton and Visser, 2009). The experimental islands support these results and the decrease in biomass toward the US zone is in line with observations from the upper end of the natural transect. When comparing MPB biomass patterns on natural transects and experimental islands, the pio zone promotes the highest MPB biomass. In the seaward direction, marine stressors increase (leading to the positive relationship between MPB biomass and elevation), and toward the salt marsh, the terrestrial stressors increase (leading to the negative relationship between elevation and MPB biomass on the experimental islands). Additionally, the experimental islands show a clear interaction between vegetation cover history and the biomass-elevation correlation, reflected by a significant interaction term in the LMM. We postulate that the initial vegetation buffered some of the negative impacts of high elevation described above (desiccation and salinity stress), leading to higher biomass in vegetated than bare islands on the highest elevational US zone. The lack of a significant effect of temperature on Chl-a biomass suggests that in addition to the other benefits already discussed, vegetation cover also has a buffering effect on temperature conditions on the islands. In addition, the positive and strong correlation between sediment water content and MPB biomass, supported the patterns on the natural transects. Though non-significant, the relationship trends of temperature and sediment water content with ENS were also similar to those on the natural transects. Hence, these results indicate that conditions on the experimental islands were comparable to those on the natural transects. Furthermore, we saw that there was no significant difference in the ENS between vegetated and bare islands, which also hints toward the buffering effect of vegetation on abiotic conditions on the experimental islands. Nevertheless, we acknowledge that the environmental parameters data set for the experimental islands was incomplete and consider the conclusions to be preliminary as the LMMs can change with the increasing number of observations and variables.

Pigment composition patterns on natural transects and experimental islands

Pigment composition revealed some distinct spatial and temporal patterns, which were more pronounced on the natural transects than on the experimental islands. Our results showed that canthaxanthin, which has been used as a marker for heterocystous cyanobacteria (Charette and Prepas, 2003; Krajewska et al., 2019), was only detected in the salt marsh (sites V and VI) and at the marine-terrestrial boundary zone (site IV). This indicates that heterocystous cyanobacteria were more distributed in salt marsh MPB communities, an observation also commonly reported for other temperate salt marsh systems (Tracy and South, 1989; Calvo et al., 1999; Janousek, 2009). However, higher plants can also be a source of canthaxanthin (Naz et al., 2021); therefore, our results should be regarded with care despite our efforts to remove all visible plant material from sampled cores. The presence of canthaxanthin only on Spiekeroog was likely due to the fact that only the Spiekeroog transects incorporated saltmarsh sites in contrast to the other two islands comprising only tidal flat samples. Alloxanthin and peridinin were present only in summer and autumn and suggest the succession of cryptophytes and dinoflagellates in the MPB communities. The succession of these algal groups in the Wadden Sea MPB communities during the warmer months was also reported by Scholz and Liebezeit (2012b), using microscopic analyses. The prevalence of fucoxanthin in our samples indicating the presence of diatoms mostly in summer did not match the results of other studies focusing on MPB pigment composition, where fucoxanthin dominated the MPB community all year round (Méléder et al., 2005; Ribeiro et al., 2021). Other microscopy-based studies found diatoms to dominate MPB communities in spring and autumn (Barranguet et al., 1997; Scholz and Liebezeit, 2012b), or in spring and summer (Semcheski et al., 2016) (together with other algae). However, the carotenoid diatoxanthin formed a substantial portion of the pigment mixture across all seasons and locations in our study, indicating the presence of diatoms, dinoflagellates, and/or haptophytes. Since diatoms are the dominant group in the Wadden Sea MPB communities, it is feasible that the diatoxanthin measured in our samples was mostly of diatom origin. In that case, diatoms indeed formed a large percentage of our MPB communities during all three seasons and at most sites, similar to the studies mentioned above. A few studies (Lucas and Holligan, 1999; Méléder et al., 2005; Jesus et al., 2009) found an increase in cyanobacteria pigments in the tidal flat sediments during summer; however, this was not observed in our study, as echinenone was present in all the three seasons. Overall, our pigment composition results showed that MPB taxonomic groups varied not just seasonally and among islands, but also along transects, which was also observed in studies by Colijn and Nienhuis (1977) and Scholz and Liebezeit (2012b).

On the experimental islands, seasonal variations in pigment composition for alloxanthin, fucoxanthin, and peridinin were not observed as on the natural transects. This suggests that MPB communities on the experimental islands were likely not of the same composition as on the transects and also not much influenced by seasons, although the LMM results showed that seasons had a substantial (accounted for >30% of the variance) effect on ENS. Pigment composition results here were non-conclusive and offered little information about the seasonal and elevation differences in MPB communities.

Overall, we would like to emphasize that the pigment composition results presented here provide only a rough insight into the community structure and are not intended to serve as a definitive species inventory. Most MPB pigment composition studies use the HPLC method and thus hampers their comparison with our data. As we can draw only limited conclusions regarding algal group composition from pigment composition using the modified GPS method, we explicitly did not include MPB community structure in our analysis, but only ENS.

Outlook

The current study enabled us to compare spatial and temporal MPB patterns in natural tidal flat transects and the salt marsh with artificial islands, and their specific manipulation of initial plant cover. Our major findings comprised of the decline in algal biomass with more stressful marine and terrestrial conditions, leading to a convergence of peak biomass in the pioneer zone of the salt marsh, right at the marine-terrestrial transition zone. We further showed that environmental stress in the terrestrial environment can be mediated by vegetation cover. Nevertheless, our study can only be a first survey generating new hypotheses that need to be tested with more detailed analyses of species composition and their traits in relation to the environmental variables. Despite some pitfalls in the modified GPS method, we believe that this method was sufficient and suitable for detecting spatial and seasonal changes in MPB biomass and the diversity of major taxonomic groups within MPB communities.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

SM, HH, SR, and IK contributed to conception and design of the study. JY and HH performed the statistical analysis. JD performed the sediment characteristics analysis. JY, MM, JD, KL, and DP collected data for the manuscript. JY wrote the first draft of the manuscript. SM, HH, and MM wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.

Funding

The DynaCom Research Unit was funded by German Research Foundation (DFG FOR 2715), the microphytobenthos subproject under the funding number MO 1931/5-1, and Bibliotheks- und Informationssystem (BIS) of the University of Oldenburg partially funded the publishing costs for this publication.

Acknowledgments

We would like to express our gratitude to the workshop staff of the Institute for Chemistry and Biology of the Marine Environment (ICBM), especially Helmo Nicolai, Gerrit Behrens, and Waldemar Siewert for their assistance in the logistics, provision of specially modified equipment and materials during our field sampling campaigns. In addition, special thanks to the staff from the Lower Saxony Wadden Sea National Park, Nationalpark-Haus Wittbülten am Internat Hermann Lietz-Schule on Spiekeroog, and the NLWKN Betriebstelle Norderney for granting us permission and access to the field, and for their logistic support and advice.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2022.956092/full#supplementary-material

Supplementary Figure S1. In-situ sediment temperatures measured along sampling transects (zone I to IV on Norderney and Wangerooge, and zone I to VI on Spiekeroog). Values were averaged values from three transects. ^*Due to a technical fault with the temperature probe, temperatures for spring on Spiekeroog could not be measured in triplicates. Data shown were derived from backup HOBO temperature loggers.

Supplementary Figure S2. In-situ sediment temperatures measured on the experimental islands of Spiekeroog, on all three elevational levels [pioneer (pio), lower salt marsh (LS) marsh, upper salt marsh (US) marsh], during the summer and autumn of 2019. Temperature data are not available for spring due to a technical fault in our measurement device.

Supplementary Figure S3. (A) General relationship between sediment mean grain size (in micrometers) and elevation (meters above mean sea level) on Spiekeroog, Norderney, and Wangerooge. (B) Correlations between elevation (meters above mean sea level) and sediment mean grain size (in micrometers) correlated individually for each of the sampled East Frisian Islands.

Supplementary Figure S4. Overview of sediment mean grain size (averaged between three transects) along the sampling sites (zones I–IV) on Norderney and Wangerooge, and along sampling sites (zones I–VI) on Spiekeroog.

Supplementary Figure S5. Overview of sediment water content (averaged between three transects) along the sampling sites (zones I–IV) on Norderney and Wangerooge, and along sampling sites (zones I-VI) on Spiekeroog.

Supplementary Figure S6. (A) General relationship between mud content (%) and elevation (meters above mean sea level) on Spiekeroog, Norderney. and Wangerooge. (B) Correlations between elevation (meters above mean sea level; m AMSL) and mud content (%) correlated individually for each of the sampled East Frisian Islands.

Supplementary Table S1. Sampling dates of the sampling campaigns in 2019.

Supplementary Table S2. List of all sampling sites on the three East Frisian Islands and corresponding GPS coordinates (in decimal degrees). Exp. Isl, experimental island; pio, pioneer zone; LS, lower saltmarsh; US, upper saltmarsh.

Supplementary Table S3. Mean and standard deviation values for abiotic environmental parameters measured and analyzed on all three islands and seasons. Standard deviation values for spring temperature data from Spiekeroog, and autumn temperature data from Norderney are unavailable, due to a technical fault in our measurement device in the field. pio, pioneer zone; LS, lower saltmarsh.

References

Admiraal, W., and Peletier, H. (1980). Influence of seasonal variations of temperature and light on the growth rate of cultures and natural populations of intertidal diatoms. Mar. Ecol. Prog. Ser. 2, 35–43. doi: 10.3354/meps002035