Next Article in Journal
Pliocene Origin, Ice Ages and Postglacial Population Expansion Have Influenced a Panmictic Phylogeography of the European Bee-Eater Merops apiaster
Next Article in Special Issue
Algal Epibionts as Co-Engineers in Mussel Beds: Effects on Abiotic Conditions and Mobile Interstitial Invertebrates
Previous Article in Journal
Muddy Boots Beget Wisdom: Implications for Rare or Endangered Plant Species Distribution Models
Previous Article in Special Issue
Individual Morphology and Habitat Structure Alter Social Interactions in a Range-Shifting Species
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 11
Issue 1
10.3390/d11010011
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

After the Fall: Legacy Effects of Biogenic Structure on Wind-Generated Ecosystem Processes Following Mussel Bed Collapse

by
John A. Commito
1,2,*,
Brittany R. Jones
1,3,
Mitchell A. Jones
1,
Sondra E. Winders
1 and
Serena Como
4
1
Environmental Studies Department, Gettysburg College, Gettysburg, PA 17325, USA
2
Unità di Biologia Marina ed Ecologia, Dipartimento di Biologia, Università di Pisa, Via Derna 1, 56126 Pisa, Italy
3
College of Fisheries and Ocean Sciences, University of Alaska Fairbanks, Fairbanks, AK 99775, USA
4
Consiglio Nazionale delle Ricerche, Istituto per l’Ambiente Marino Costiero (CNR-IAMC), Località Sa Mardini, Torregrande, 09170 Oristano, Italy
*
Author to whom correspondence should be addressed.
Diversity 2019, 11(1), 11; https://doi.org/10.3390/d11010011
Submission received: 11 December 2018 / Revised: 9 January 2019 / Accepted: 11 January 2019 / Published: 15 January 2019
(This article belongs to the Special Issue Diversity of Ecosystem Engineers in the World Coasts and Oceans)
Browse Figures
Versions Notes

Abstract

:
Blue mussels (Mytilus edulis) are ecosystem engineers with strong effects on species diversity and abundances. Mussel beds appear to be declining in the Gulf of Maine, apparently due to climate change and predation by the invasive green crab, Carcinus maenas. As mussels die, they create a legacy of large expanses of shell biogenic structure. In Maine, USA, we used bottom traps to examine effects of four bottom cover types (i.e., live mussels, whole shells, fragmented shells, bare sediment) and wind condition (i.e., days with high, intermediate, and low values) on flow-related ecosystem processes. Significant differences in transport of sediment, meiofauna, and macrofauna were found among cover types and days, with no significant interaction between the two factors. Wind condition had positive effects on transport. Shell hash, especially fragmented shells, had negative effects, possibly because it acted as bed armor to reduce wind-generated erosion and resuspension. Copepods had the greatest mobility and shortest turnover times (0.15 d), followed by nematodes (1.96 d) and the macrofauna dominant, Tubificoides benedeni (2.35 d). Shell legacy effects may play an important role in soft-bottom system responses to wind-generated ecosystem processes, particularly in collapsed mussel beds, with implications for recolonization, connectivity, and the creation and maintenance of spatial pattern.

Graphical Abstract

1. Introduction

Blue mussels (Mytilus edulis) are ecosystem engineers with strong effects on sediment, species composition, and abundances of meiofauna and macrofauna compared to soft-bottoms without mussels [1,2,3,4,5,6,7,8,9,10]. Mussel beds can experience dramatic short- and long-term changes in abundance [11,12,13,14,15]. They appear to be declining in the Gulf of Maine, apparently due to climate change [16,17,18,19,20,21] and predation by the invasive green crab, Carcinus maenas [22,23,24,25]. When mussels die, they leave behind large expanses of non-living biogenic material in the form of shell hash [3,4,26] (Figure 1A).
Mussel beds are not a simple matrix of live mussels. At our research sites in Maine, USA, mussel beds have a fractal, hierarchical spatial structure [26,27,28,29,30,31]. They comprise a complex patchwork of live mussels, bare sediment, and whole (empty) and fragmented shells [3,4,26] (Figure 1B). The influence of these different bottom cover types on flow-related ecosystem processes depends on a variety of patch attributes [4], but remains poorly known. Live mussels and shell hash are roughness elements that project above the sediment surface, creating irregular surface topography that alters water movement over the bottom [4,28]. They convert laminar flow to turbulent flow, where the boundary layer regime affects the transport and resuspension of particles [32,33,34,35,36].
Live mussels and shell hash have significant effects on sediment flux, including the overall mass moved, grain-size proportions, and organic content, as well on univariate and multivariate measures of meiofauna transport, with a smaller influence on macrofauna [2,4]. Rates of sediment flux and animal transport in bedload tend to be positively correlated with wind velocity and the fetch over which the wind blows when water covers the site [37,38,39,40,41,42,43,44,45,46,47]. How does this wind-forcing interact with live mussel, whole shell, fragmented shell, and bare sediment cover types in mussel beds? Does wind-forcing operate independently of bottom cover type, or does a significant interaction occur between wind and cover type?
The answers are important because this type of local dispersal plays a major role in regulating population dynamics and persistence in spatially heterogeneous systems like intertidal soft-bottoms [39,48,49,50,51,52]. Wind-generated transport carries sediment and organisms to nearby locations on a steadily recurring basis. It provides a mechanism for juveniles and adults to colonize suitable sites on each tidal cycle rather than only after larvae are produced in a reproductive event. Wind-generated transport may be particularly important for benthic species that do not produce free-swimming larvae, including meiofauna and some macrofauna [2,4,38,39,42,43,53,54,55,56,57,58].
For this paper we took advantage of an archived data set from a Maine, USA, soft-bottom, intertidal Mytilus edulis bed that has undergone a recent decline in live mussel abundance [4]. We partitioned bottom trap data from that study into days with different wind condition values and calculated three measures of animal dispersal. We found that both wind-forcing and biogenic structure had significant effects on sediment flux and bottom fauna dispersal. Somewhat surprisingly, we discovered no significant interaction between wind condition and bottom cover type, suggesting that they acted independently during the short time course of this investigation. The results provide evidence that changes in biogenic structure after mussel bed collapse have significant impacts on transport of sediment and animals in soft-bottom systems in the Gulf of Maine and elsewhere.

2. Materials and Methods

2.1. Study Site

The research was conducted at Carrying Place Cove, Harrington, Maine, USA (44.5451° N, −67.7844° W), a relatively protected embayment with a bottom of muddy sand (Figure 2A). As is typical in this region, the intertidal mussel bed extends across the mouth of the cove near the low tide line. Long-term time series of aerial photographs suitable for mussel cover analysis do not exist for Carrying Place Cove. However, our twice-per-year field observations at the site over the last 30 years indicate a precipitous drop in mussel abundance at Carrying Place Cove. At the time of this investigation in 2011, digital analysis of photo-transects at the site revealed a scattered array of small, isolated patches of live mussels accounting for <4% of the bottom cover within the bed [4]. Whole and fragmented shells combined covered approximately 20%, and the remaining portion was bare sediment [4].

2.2. Wind Condition

For each of three successive days when the fieldwork was conducted at the site (15–18 May 2011; see below), we created a qualitative wind condition index based on wind velocity and fetch length over which the wind blew. Wind velocity and direction data were obtained from National Ocean Service Station ATGM1-8413320, Bar Harbor, Maine, ≈50 km from our study site. For 6-minute intervals when water was covering the bed on each day that traps were deployed (15–18 May 2011), wind velocity was calculated as mean ± SE and wind direction was plotted on rose diagrams (Figure 2B). For wind velocity, the ranks (highest to lowest velocity) were day 1 > day 2 > day 3. For fetch length, the ranks (longest to shortest length) were day 3 > day 1 > day 2. Assigning equal weight to both factors, the combined wind condition index ranks (strongest to weakest) were day 1 > day 3 > day 2.

2.3. Field and Laboratory Procedures

We used archived data from the study site [4]. Ten patches of each of the four cover types were selected within a 10-m wide band that ran from the upper margin to the lower margin of the center of the bed, using Hurlbert’s systematic model [59] to achieve treatment dispersion and reduce the likelihood of segregation, type 1 error, and spurious treatment effects. Large patches of mussels like those used previously [2] could not be found at the site, so patches ≥ 0.5 m × 0.5 m of all four cover types were utilized.
At low tide on 15 May 2011, a 1.3-cm diameter (cross-sectional area = 1.33 cm2) core sample was taken 5-cm deep in each of the 40 patches. Cores of this size have been used successfully to sample sediment characteristics, meiofauna, and macrofauna at similar sites in Maine and elsewhere. They collect approximately the same numbers of macrofaunal species as larger cores and allowed us to utilize a larger sample size with the same sampling and laboratory processing effort [2,4,54,55,57]. A 15-cm tall bottom trap of the same diameter (aspect ratio = 11.5:1) was filled with seawater and inserted flush with the sediment surface into each of the holes where the cores were extracted. Because trap capture efficiency declines as traps fill with sediment, in the original study each trap was replaced with a new trap in the same hole on each successive day, and the cumulative contents constituted a single sample for the three-day deployment. For this paper, we took advantage of that sampling strategy and utilized each new trap on each day as an individual sample. To achieve sample independence, we randomly selected three traps from each of the four cover types on each of the three days so that no trap location was used more than once.
Core and trap contents were stained with rose bengal and fixed in situ with buffered formalin. In the laboratory, samples were wet sieved (0.5 mm mesh for macrofauna and the coarse sediment fraction; 0.063 mm mesh for meiofauna and the fine sediment fraction), and the material that passed through the sieves (silt-clay sediment fraction) was retained. Meiofauna (to family or higher levels, e.g., Copepoda, Nematoda, Foraminifera) and macrofauna (to species level when possible [60]) were identified and counted in gridded Petri dishes. The three sediment size classes were dried at 85 °C for 24 h and weighed. Organic matter was calculated using loss-on-ignition by burning at 500 °C for 4 h [61]. Core contents were analyzed and presented elsewhere [4] and used here only to calculate relative dispersal rate (see below).
Trap contents (Table S1) were used to determine rates of sediment and animal movement [2,4,44,57,62]:
  • Sediment flux rate = g sediment trap−1 d−1
  • Absolute dispersal rate (ADR) = number of individuals trap−1 d−1
  • Relative dispersal rate (RDR) = number of individuals trap−1 d−1 ambient individual−1, calculated by dividing the number of individuals collected in a trap by the number of individuals collected in the core from the same location. Relative dispersal rate normalizes for ambient density. It is equivalent to per capita dispersal. For example, RDR > 1 occurs when the number of individuals in a trap is larger than the number of individuals in its corresponding core. Note that in some cases the number of macrofauna individuals in a core = 0, so RDR was undefined due to the 0 in the denominator. Those samples could not be used, resulting in unequal sample sizes. Our statistical analysis relied on a balanced design, so RDR was not calculated for macrofauna.
  • Bulk dispersal rate (BDR) = number of individuals g sediment−1 trap−1 d−1, calculated by dividing the number of individuals collected in a trap by the sediment mass collected in that trap. It is a measure of dispersal per unit of transported sediment, providing information on how tightly linked animal movement is to passive bedload transport.

2.4. Data Analysis

Analysis of trap sediment included mass for each grain-size, total mass, and percent organic matter. For meiofauna and macrofauna, analysis included taxon richness (for meiofauna) or species richness (for macrofauna), Shannon–Weiner Index (H’) using log base e, and the absolute, relative (meiofauna only), and bulk dispersal rates for each dominant taxon or species and the total number of individuals. Because the oligochaete Tubificoides benedeni often dominates the macrofauna in mussel beds and exhibits abundance patterns different from other macrofaunal species [2,3,4], the non-oligochaetes were also analyzed together as one group. Differences were analyzed using a 2-factor ANOVA design with cover types (4 levels: fixed) and days (3 levels: fixed) as orthogonal factors, with 3 replicates. After application of Cochran’s C-test to check for homogeneity of variances, data were transformed if necessary, as indicated in the ANOVA tables. When significant differences were found (alpha = 0.05), a posteriori comparisons were made using Student-Newman-Keuls (SNK) tests. Cases occurred where SNK could not discriminate among alternative hypotheses despite significant ANOVA differences. Because a posteriori multiple comparisons have less power than the original ANOVA F-test, such cases indicated that the 2 groups with the largest and smallest means within the factor were significantly different, and no further resolution among the groups in that factor could be determined [63]. Analyses were performed with the StatSoft 6.1 computer program [64].
The dispersing assemblages were analyzed with a distance-based permutational multivariable analysis of variance (PERMANOVA) based on Bray–Curtis dissimilarity measures [65,66]. The PERMANOVA used the same design as the ANOVA (see above) and was executed with the “adonis” function in the “vegan” package for R v3.4.2 [67]. The non-metric multidimensional scaling (nMDS) ordination model based on the Bray–Curtis dissimilarity matrix was conducted with the PRIMER v5.2 package [68] and used to plot the differences in meiofauna and macrofauna assemblages among cover types and days.

3. Results

3.1. Sediment

Significant differences in trap sediment variables were found among cover types and days, with no significant interaction between the two factors (Table 1, Figure 3A,B). Among cover types, silt-clay mass was higher in live mussels than in fragmented shells, with no further discrimination possible. Among days, more fine sediment and total sediment was captured on the day with the strongest wind condition rank (day 1) than on the other days, which were not different from each other. The percentage of organic material was highest on the day with the second strongest wind condition rank (day 3), followed in order by day 1 and day 2. Overall, cover types did not matter as much as days, where the highest total sediment flux rate and organic content were on the days with strongest wind condition ranks.

3.2. Meiofauna

The traps collected eight meiofauna taxa and were dominated by copepods (55.6%) and nematodes (39.8%), followed by foraminiferans (3.9%) and much smaller numbers of other taxa (Appendix A). Significant differences in meiofauna variables were found among cover types and days, with no significant interactions between the two factors.
Among cover types (Table 2, Figure 4A,B), H’ was higher in bare sediment than in live mussels, with no further discrimination possible. Absolute dispersal rates (ADRs) for copepods and the total assemblage were lower in fragmented shells than in the other cover types, which were not different from each other. The ADR for nematodes was higher in whole shells than in fragmented shells, with no further discrimination possible. The ADR for foraminiferans was higher in bare sediment than in fragmented shells, with no further discrimination possible. Among days, H’ and ADR for copepods, nematodes, foraminiferans, and the total assemblage were generally highest on the days with the strongest wind condition ranks (day 1 and day 3, especially day 1) and lowest on day 2 (Table 2, Figure 4A,B). Overall, fragmented shells depressed ADR values, and wind increased them.
At the multivariate level, PERMANOVA revealed significant differences in the structure of the dispersing meiofauna assemblage among cover types and days, with no significant interactions between the 2 factors (Table 3). The nMDS plot showed separation among cover types and days (Figure 5A), highlighting the PERMANOVA results. Dissimilarity values were large for fragmented shells compared to the other cover types, as well as for the weakest wind condition day (day 2) compared to day 1 and day 3 with stronger wind condition ranks, consistent with the univariate ANOVA results for ADR.
Relative dispersal rates (RDR) were > 1, indicating high mobility relative to ambient densities, particularly for copepods (Table 4, Figure 4C). RDR among cover types for copepods and the total assemblage were higher in whole shells than in the other cover types, which were not different from each other. RDR was higher for foraminiferans in whole shells than in fragmented shells, with no further discrimination possible, and higher on the day with the strongest wind condition rank (day 1) than on the other 2 days, which were not different from each other. These foraminiferan results should be interpreted with caution because the RDR data for foraminiferans were heterogeneous despite several types of transformation. They are presented because ANOVA is robust to deviations from the assumption of homogeneity of variance when sample sizes are equal [69]. Overall, whole shells and days with high wind condition ranks had the strongest effects on meiofauna RDR. Both increased RDR values.
Bulk dispersal rate (BDR) among cover types for nematodes was higher in bare sediment than in fragmented shells, with no further discrimination possible (Table 5, Figure 4D). Among days, BDR for copepods was lower on the day with the strongest wind condition rank (day 1) than on the other days, which were not different from each other. Overall, cover types and days had weak effects on BDR compared to their effects on ADR and RDR. Each unit of moving sediment was transported with similar numbers of animals, regardless of cover type or day.

3.3. Macrofauna

The traps collected 19 macrofauna species and were dominated by the oligochaete Tubificoides benedeni (42.3%), with smaller numbers of polychaetes, an isopod, and other species (Appendix A). None of the macrofaunal variables showed significant differences among cover types or days except for BDR, where non-oligochaetes and the total assemblage had highest values on the weakest wind condition day (day 2), and the other days were not different from each other (Table 6 and Table 7, Figure 6A–C).
At the multivariate level, PERMANOVA revealed no significant differences in the structure of the dispersing macrofaunal assemblage among cover types or days (Table 3B). The nMDS plot showed no separation among cover types and days (Figure 5B). The multivariate and univariate results were consistent in showing little effect of cover types or days on the dispersing macrofaunal assemblage.

4. Discussion

4.1. Cover Types

Fragmented shells were the cover type with the largest impact on sediment flux and animal dispersal. Of the 10 univariate cases with significant ANOVA differences among cover types, fragmented shells were responsible for the lowest value in seven of them (silt-clay; ADR of total meiofauna, copepods, nematodes, foraminiferans; RDR of foraminiferans; BDR of nematodes). Fragmented shells also had the largest Bray–Curtis dissimilarity values among cover types revealed by the multivariate PERMANOVA.
Why were fragmented shells so important? Some investigations have shown that mollusk shell material can increase abundance and diversity in benthic systems [36,70,71,72,73,74,75,76,77,78]. Other studies found mixed, weak, or no effects [40,79,80,81,82,83]. But the effects on flow-related ecosystem processes have not been well-studied. What is it about fragmented shells in particular that caused its flow-related effects at our site? Shell fragments are often observed with their edges projecting upward from the bed [28]. Rough surface elements like these would be expected to increase erosion and resuspension of sediment and animals due to higher turbulent energy and bed shear stress [35,84,85], with reduced flow velocity leading to higher rates of sediment and animal deposition [36].
But fragmented shells had lower rates, not higher. A possible reason is that shell fragments, like terrestrially-derived pebbles and gravel [86], can form a dense layer that acts as bed armor to protect the sediment beneath it from erosion and transport [4]. Because fragments vary in size and shape, they can pack together tightly, with small pieces fitting into the interstices between larger pieces, as predicted by packing theory [87]. This configuration is not possible with whole shells. Whole shells do not exist in as wide a range of sizes, and they are curved, with gaps between shells even when in contact with each other [28]. Compared to shell fragments, whole shells form a more porous, loosely packed layer with exposed sediment between the biogenic elements. In addition, whole shells tend to lie flat on the bed, with a tendency for the convex side to face upward, thus presenting a relatively smooth surface over which water may flow [28].
This argument is not meant to downplay the effects of whole shells. Of the other three significant ANOVA cases, whole shells had the lowest value in two of them (total meiofauna and copepods). In addition, PERMANOVA showed that whole shells had Bray–Curtis dissimilarity values second only to those for fragmented shells. So, both forms of shell hash had significant depressing effects on flow-related ecosystem processes.
Live mussels have been shown in some other studies to affect ambient sediment and community structure [1,2,3,4,88]. So, it might make sense that they would have effects on flow-related ecosystem processes. Yet live mussels were responsible for only two significant ANOVA differences among cover types, the highest silt-clay mass and lowest value of H’. Mussels produce copious feces and pseudofeces and project up above the bed surface as roughness elements [28], both of which factors may explain the significant silt-clay result. Increased rates of sediment deposition have been observed in beds of Mytilus edulis [2,88,89,90], other mussel species [91,92], and pinnid bivalves [93,94]. Live mussel cover also represents a pool of species different from the other cover types [1,2,3,4,88], which may explain the H’ result. The small overall impact of live mussels may have been due to wind conditions that produced turbulent mixing sufficient to lift sediment and animals away from the bottom and export them elsewhere [95]. That kind of interplay might be expected to produce significant interactions between cover types and days, which we never observed.
Bare sediment was responsible for three significant ANOVA difference among cover types (the highest values of H’, ADR of foraminiferans, and BDR of nematodes). As explained above, these results reflect the depressing effect of other cover types on these three variables compared to bare sediment, the dominant cover type at our study site [4].

4.2. Days and Wind Condition

Our results revealed rapid sediment flux and animal dispersal responses to changes in wind condition. Wind condition was strongest on day 1, weakest on day 2, and intermediate on day 3. Had it increased (or decreased) monotonically during the 3-day period, any daily wind condition effects could have been smoothed or obscured by the sustained addition to (or removal from) the water column of resuspended sediment and animals. But wind condition changed non-monotonically over the 3-day period. Concomitant with these changing daily wind conditions, we were able to observe significant daily changes in sediment flux and animal transport. Of the 12 univariate cases with significant ANOVA differences among days, day 1 (strongest wind condition) was responsible for the highest value in seven of them (fine sediment; total sediment; meiofauna H’; ADR of nematodes, foraminiferans; RDR of foraminiferans; BDR of copepods). Day 2 (weakest wind condition) was responsible for the lowest value in four of them (percentage of organic material; ADR of total meiofauna, copepods, foraminiferans) and the highest value in two of them, the latter two being the only significant effects observed for macrofauna (BDR of total macrofauna and non-oligochaetes). Day 2 had the largest Bray–Curtis dissimilarity values among cover types revealed by the multivariate PERMANOVA, followed by day 1. Day 3 (intermediate wind condition) had intermediate dissimilarity values, and it never produced any highest or lowest significantly different univariate responses.
Sediment flux and animal transport tend to be positively correlated with wind condition [37,38,40,41,42,43,44,45,46,47,62]. So, our results are consistent with this general pattern. However, bottom hydrodynamics in mussel beds are complex, with laminar flow being converted to turbulent flow over the roughness elements, especially as flow velocity increases [32,35]. Thus, it seemed likely that significant interactions would occur between cover types and days due to varying wind-driven water velocities over the different kinds of roughness elements. Yet neither the ANOVA nor the PERMANOVA results revealed such interactions.

4.3. Modes of Dispersal, Mobility, and Turnover Times

Copepods and nematodes were by far the two most abundant meiofauna taxa in cores and traps at our site, with foraminiferans a distant third. Copepods engage in active emergent behavior and subsequent transport, whereas nematodes are eroded from the bed and transported passively in bedload [53,57]. Relative dispersal rates (RDRs) were >1 individual trap−1 d−1 ambient individual−1 for copepods, indicating that more individuals were caught in a trap in a single day than were located in a core of the same diameter in the same location. Their high per capita mobility contrasts with that of nematodes, which had RDR < 1 individual trap−1 d−1 ambient individual−1. In addition, BDR values were higher for copepods than for nematodes, consistent with the RDR results.
In this study, the overall copepod/nematode RDR ratio = (6.64 individuals trap−1 d−1 ambient individual−1)/(0.51 individuals trap−1 d−1 ambient individual−1) = 13.0. For comparison, we calculated RDR values for copepods and nematodes from data available in other studies. Whether RDR values were larger or smaller than those in this study, their ratios all indicated higher per capita mobility for copepods than for nematodes: intertidal mudflats in Maine [57] (10.09/0.11 = 91.7) and Louisiana, USA [96] (0.97/0.07 = 14.0) and a high energy intertidal sandflat in North Carolina, USA [54] at high tide (0.066/0.013 = 5.14) and flood tide (0.256/0.126 = 2.03).
Macrofauna RDR values were not calculated using individual trap and core values for cover types and days in this study, as explained above. However, rough estimates of RDR based on summed values of traps and cores for combined cover types and days were 0.21 individuals trap−1 d−1 ambient individual−1 for the oligochaete Tubificoides benedeni (by far the dominant macrofauna species) and 0.87 individuals trap−1 d−1 ambient individual−1 for the non-oligochaetes as a group. These estimates are in agreement with RDR values < 1 individual trap−1 d−1 ambient individual−1 for T. benedeni and a variety of macrofauna elsewhere [2,44,45,62]. These values demonstrate on a per capita basis that macrofauna were somewhat less mobile than passively dispersing nematodes and much less mobile than actively emerging copepods. Macrofauna probably have lower rates of erosion and transport than meiofauna because of their larger size and deeper life position within the sediment. Oligochaetes in particular burrow downwards in response to flow-induced sediment erosion [97].
Turnover time is the number of days it would take for the individuals in an ambient community core to be replaced by new individuals arriving by post-larval transport and settling to the bottom. It is mathematically equivalent to the inverse of the RDR. The RDR, hence turnover time, varied among cover types and days at our research site, but turnover times based on overall mean RDR values were 0.15 d for copepods, 1.96 d for nematodes, 0.55 d for foraminiferans, and 2.35 d for T. benedeni. These turnover times point out the highly dynamic nature of our intertidal, soft-bottom research site, despite its location in a low-energy, muddy-sand cove. Meiofauna and macrofauna are on the move with sediment every day, with implications for recolonization, connectivity, and the creation and maintenance of spatial pattern across scales [47,52].

4.4. Space and Time in Collapsed Mussel Beds

Some of the largest, densest mussel beds where we have conducted research since the mid-1970s, including Carrying Place Cove, contained virtually no live mussels at the time of this study in 2011, and we have observed no rebound since then. The Gulf of Maine sea surface temperature has risen faster in the past decade than in 99% of the world ocean [98]. The region has the world’s fastest sea level rise, approximately 3–4 times faster than the world average [99], drowning saltmarshes [100] and posing a serious threat to the very existence of intertidal mudflats [101]. As mussel beds collapse in the Gulf of Maine and elsewhere due to climate change [16,17,18,19,20,21] and predation by the invasive green crab, Carcinus maenas [22,23,24,25], at what density and spatial arrangement of individuals will we decide they have ceased to exist as ecosystem engineers? Is it possible that the formerly robust mussel bed at Carrying Place Cove has been depleted of mussels to such a degree that it can no longer even be called a mussel bed? Will there be a time in the future when the different cover types will be so homogenized that biogenic structure ceases to matter?
These existential questions about this important ecosystem engineer have no firm answers. Maine mussel beds could return to their previous levels in response to long-term oscillating oceanographic conditions in the northwest Atlantic region [102,103,104]. Mussel beds have been observed elsewhere to bounce back relatively quickly after population declines [11,12,13,14,15]. Or mussel beds may shift to a new configuration. In northern Europe where Pacific oysters, Magallana (formerly Crassostrea) gigas, have invaded M. edulis beds, both bivalve species can coexist because the oysters provide mussels with a physically complex spatial refuge from predators [105].
In Maine, where invasive oysters have not been observed, a different path is likely. It appears that mussels are producing high cover values of whole or fragmented shells [4], depending on the mode of death, e.g., one or both valves broken by crushing predators, neither valve broken by physiological senescence or thermal stress [106,107,108]. Over time, shells are exported, buried, or broken down into smaller and smaller pieces, similar in size to terrestrially-derived sand and mud, resulting in bare sediment [106,107,108]. During this taphonomic time course, we predict a commensurate shift in ecosystem processes affected by flow over the bottom, including sediment flux and animal transport. We predict an initial period with lowered rates of these ecosystem processes because of high bottom cover values of shell hash, particularly fragmented shells, and then a rebound to the rates observed over bare sediment. Because mussels can cover extensive areas in soft-bottom systems [14,29], even subtle changes in these ecosystem processes may assume great importance when integrated over the large spatial extent of potentially collapsing mussel beds worldwide.

Supplementary Materials

The following are available online at https://www.mdpi.com/1424-2818/11/1/11/s1, Table S1: Sediment, meiofauna, and macrofauna data for bottom traps from Carrying Place Cove, Maine, USA.

Author Contributions

Conceptualization, J.A.C, B.R.J., M.A.J., S.E.W., S.C.; Methodology, J.A.C, B.R.J., M.A.J., S.E.W., S.C.; Software, S.C.; Validation, J.A.C, B.R.J., M.A.J., S.E.W., S.C.; Formal Analysis, J.A.C, B.R.J., M.A.J., S.E.W., S.C.; Investigation, J.A.C, B.R.J., M.A.J., S.E.W., S.C.; Resources, J.A.C, B.R.J., M.A.J., S.E.W., S.C.; Data Curation, J.A.C, B.R.J., M.A.J., S.E.W., S.C.; Writing-Original Draft Preparation, J.A.C, B.R.J., M.A.J., S.E.W., S.C.; Writing-Review & Editing, JJ.A.C, B.R.J., M.A.J., S.E.W., S.C.; Visualization, J.A.C, B.R.J., M.A.J., S.E.W., S.C.; Supervision, J.A.C.; Project Administration, J.A.C.; Funding Acquisition, J.A.C, B.R.J., M.A.J., S.E.W.

Funding

This research was funded by Howard Hughes Medical Institute grant number 52007540, Gettysburg College Research and Professional Development Program, and Gettysburg College Senior Honors Research Program. The APC was funded by the journal.

Acknowledgments

We thank Angela Commito and Ann Commito for assistance in the field and the University of Pisa and Gettysburg College for administrative and technical support.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Appendix A

Table
Table A1. Macrofauna and meiofauna captured in traps (overall mean + SE individuals per trap) by rank order. For macrofauna: A = Arthropoda, M = Mollusca, N = Nemertea, O = Oligocheata, P = Polychaeta.
Table A1. Macrofauna and meiofauna captured in traps (overall mean + SE individuals per trap) by rank order. For macrofauna: A = Arthropoda, M = Mollusca, N = Nemertea, O = Oligocheata, P = Polychaeta.
Taxon RankMeanSE
Meiofauna
Copepoda 1268.9439.57
Nematoda 2192.0645.47
Foraminifera 319.223.76
Ostracoda 41.530.36
Kinorhyncha 50.420.13
Halacaridae 60.280.18
Cumacea 70.080.06
Acarina 80.030.03
Macrofauna
Tubificoides benedeniO11.080.27
Lepidonotus squamatusP20.420.14
Polychaeta errant juvenile juvenileP30.220.08
Isopoda sp. AA40.170.07
Ampharetidae sp. AP50.110.05
Capitella capitataP60.060.04
Polycheata sp. A juvenile P60.060.04
Mya arenariaM60.060.04
Ampelisca abditaA60.060.04
Carcinus maenasA60.060.04
Gammaridae sp. AA60.060.04
Littorina littoreaM60.060.04
Phyllodocidae sp. AP60.060.04
Lineus viridisN140.030.03
Hetermastus filiformisP140.030.03
Gammarus finmarchicusA140.030.03
Polycheata sp. B juvenile P140.030.03
Gammarus oceanicusA140.030.03
Polycheata sp. C juvenileP140.030.03

References

  1. Commito, J.A. Adult larval interactions: Predictions, mussels, and cocoons. Estuar. Coast. Shelf Sci. 1987, 25, 599–606. [Google Scholar] [CrossRef]
  2. Commito, J.A.; Celano, E.A.; Celico, H.J.; Como, S.; Johnson, C.P. Mussels matter: Postlarval dispersal dynamics altered by a spatially complex ecosystem engineer. J. Exp. Mar. Biol. Ecol. 2005, 316, 133–147. [Google Scholar] [CrossRef]
  3. Commito, J.A.; Como, S.; Grupe, B.M.; Dow, W.E. Species diversity in the soft-bottom intertidal zone: Biogenic structure, sediment, and macrofauna across mussel bed spatial scales. J. Exp. Mar. Biol. Ecol. 2008, 366, 70–81. [Google Scholar] [CrossRef]
  4. Commito, J.A.; Jones, B.R.; Jones, M.A.; Winders, S.E.; Como, S. What happens after mussels die? Biogenic legacy effects on community structure and ecosystem processes. J. Exp. Mar. Biol. Ecol. 2018, 506, 30–41. [Google Scholar] [CrossRef]
  5. Thiel, M.; Ullrich, N. Hard rock versus soft bottom: The fauna associated with intertidal mussel beds on hard bottoms along the coast of Chile, and considerations on the functional role of mussel beds. Helgol. Mar. Res. 2002, 56, 21–30. [Google Scholar] [CrossRef]
  6. Norling, P.; Kautsky, N. Structural and functional effects of Mytilus edulis on diversity of associated species and ecosystem functioning. Mar. Ecol. Prog. Ser. 2007, 351, 163–175. [Google Scholar] [CrossRef]
  7. Norling, P.; Kautsky, N. Patches of the mussel Mytilus sp. are islands of high biodiversity in subtidal sediment habitats in the Baltic Sea. Aquat. Biol. 2008, 4, 75–87. [Google Scholar] [CrossRef]
  8. Bouma, T.J.; Olenin, S.; Reise, K.; Ysebaert, T. Ecosystem engineering and biodiversity in coastal sediments: Posing hypotheses. Helgol. Mar. Res. 2009, 63, 95–106. [Google Scholar] [CrossRef]
  9. Buschbaum, C.; Dittmann, S.; Hong, J.-H.; Hwang, I.-S.; Strasser, M.; Thiel, M.; Valdiva, N.; Yoon, S.-P.; Reise, K. Mytilid mussels: Global habitat engineers in coastal sediments. Helgol. Mar. Res. 2009, 63, 47–58. [Google Scholar] [CrossRef]
  10. Gutiérrez, J.L.; Jones, C.G.; Byers, J.E.; Arkema, K.K.; Berkenbusch, K.; Commito, J.A.; Duarte, C.M.; Hacker, S.D.; Lambrinos, J.G.; Hendriks, I.E.; et al. Physical ecosystem engineers and the functioning of estuaries and coasts. In Treatise on Estuarine and Coastal Science; Wolanski, E., McLusky, D.S., Eds.; Elsevier Academic Press: Waltham, MA, USA, 2011; Volume 7, pp. 53–81. ISBN 9780080878850. [Google Scholar]
  11. Nehls, G.; Thiel, M. Large-scale distribution patterns of the mussel Mytilus edulis in the Wadden Sea of Schleswig–Holstein: Do storms structure the ecosystem? Neth. J. Sea Res. 1993, 31, 181–187. [Google Scholar] [CrossRef]
  12. Reusch, T.B.H.; Chapman, A.R.O. Persistence and space occupancy by subtidal blue mussel patches. Ecol. Monogr. 1997, 67, 65–87. [Google Scholar] [CrossRef]
  13. Hertweck, G.; Liebezeit, G. Historic mussel beds (Mytilus edulis) in the sedimentary record of a back-barrier tidal flat near Spiekeroog Island, southern North Sea. Helgol. Mar. Res. 2002, 56, 51–58. [Google Scholar] [CrossRef]
  14. Folmer, E.O.; Drent, J.; Troost, K.; Büttger, H.; Dankers, N.; Jansen, J.; van Stralen, M.; Millat, G.; Herlyn, M.; Philippart, C.J. Large-scale spatial dynamics of intertidal mussel (Mytilus edulis L.) bed coverage in the German and Dutch Wadden Sea. Ecosystems 2014, 17, 550–566. [Google Scholar] [CrossRef]
  15. Khaitov, V.M.; Lentsman, N.V. The cycle of mussels: Long-term dynamics of mussel beds on intertidal soft bottoms at the White Sea. Hydrobiologia 2016, 781, 161–180. [Google Scholar] [CrossRef]
  16. Helmuth, B.; Mieszkowska, N.; Moore, P.; Hawkins, S.J. Living on the edge of two changing worlds: Forecasting the responses of rocky intertidal ecosystems to climate change. Ann. Rev. Ecol. Evol. Syst. 2006, 37, 373–404. [Google Scholar] [CrossRef]
  17. Jones, S.J.; Lima, F.P.; Wethey, D.S. Rising environmental temperatures and biogeography: Poleward range contraction of the blue mussel, Mytilus edulis L., in the western Atlantic. J Biogeogr. 2010, 37, 2243–2259. [Google Scholar] [CrossRef]
  18. Sorte, C.J.B.; Jones, S.J.; Miller, L.P. Geographic variation in temperature tolerance as an indicator of potential population responses to climate change. J. Exp. Mar. Biol. Ecol. 2011, 400, 209–217. [Google Scholar] [CrossRef]
  19. Sorte, C.J.B.; Davidson, V.E.; Franklin, M.C.; Benes, K.M.; Doellman, R.J.; Etter, R.J.; Hannigan, R.E.; Lubchenco, J.; Menge, B.A. Long-term declines in an intertidal foundation species parallel shifts in community composition. Glob. Chang. Biol. 2016. [Google Scholar] [CrossRef]
  20. Lesser, M.P. Climate change stressors cause metabolic depression in the blue mussel, Mytilus edulis, from the Gulf of Maine. Limnol. Oceanogr. 2016, 61, 1705–1717. [Google Scholar] [CrossRef]
  21. Steeves, L.E.; Filgueira, R.; Guyondet, T.; Chassé, J.; Comeau, L. Past, present, and future: Performance of two bivalve species under changing environmental conditions. Front. Mar. Sci. 2018, 5, 184. [Google Scholar] [CrossRef]
  22. Grosholz, E.D.; Ruiz, G.M. Predicting the impact of introduced marine species: Lessons from the multiple invasions of the European green crab Carcinus maenas. Biol. Conserv. 1996, 78, 59–66. [Google Scholar] [CrossRef]
  23. Whitlow, W.L.; Grabowski, J.H. Examining how landscapes influence benthic community assemblages in seagrass and mudflat habitats in southern Maine. J. Exp. Mar. Biol. Ecol. 2012, 411, 1–6. [Google Scholar] [CrossRef]
  24. Petraitis, P.S.; Dudgeon, S.R. Variation in recruitment and the establishment of alternative community states. Ecology 2015, 96, 3186–3196. [Google Scholar] [CrossRef] [PubMed]
  25. Tan, E.B.P.; Beal, B.F. Interactions between the invasive European green crab, Carcinus maenas (L.), and juveniles of the soft-shell clam, Mya arenaria L., in eastern Maine, USA. J. Exp. Mar. Biol. Ecol. 2015, 462, 62–73. [Google Scholar] [CrossRef]
  26. Commito, J.A.; Dow, W.E.; Grupe, B.M. Hierarchical spatial structure in soft-bottom mussel beds. J. Exp. Mar. Biol. Ecol. 2006, 330, 27–37. [Google Scholar] [CrossRef]
  27. Snover, M.L.; Commito, J.A. The fractal geometry of Mytilus edulis L. spatial distribution in a soft-bottom system. J. Exp. Mar. Biol. Ecol. 1998, 223, 53–64. [Google Scholar] [CrossRef]
  28. Commito, J.A.; Rusignuolo, B.R. Structural complexity in mussel beds: The fractal geometry of surface topography. J. Exp. Mar. Biol. Ecol. 2000, 225, 133–152. [Google Scholar] [CrossRef]
  29. Crawford, T.W.; Commito, J.A.; Borowik, A.M. Fractal characterization of Mytilus edulis L. spatial structure in intertidal landscapes using GIS methods. Landsc. Ecol. 2006, 21, 1033–1044. [Google Scholar] [CrossRef]
  30. Commito, J.A.; Commito, A.E.; Platt, R.V.; Grupe, B.M.; Piniak, W.E.D.; Gownaris, N.J.; Reeves, K.A.; Vissichelli, A.M. Recruitment facilitation and spatial pattern formation in soft-bottom mussel beds. Ecosphere 2014, 5, art160. [Google Scholar] [CrossRef]
  31. Commito, J.A.; Gownaris, N.J.; Haulsee, D.E.; Coleman, S.E.; Beal, B.F. Separation anxiety: Mussels self-organize into similar power-law clusters regardless of predation threat cues. Mar. Ecol. Prog. Ser. 2016, 547, 107–119. [Google Scholar] [CrossRef]
  32. Butman, C.A.; Fréchette, M.; Geyer, W.R.; Starczak, V.R. Flume experiments on food supply to the blue mussel Mytilus edulis L. as a function of boundary-layer flow. Limnol. Oceanogr. 1994, 39, 1755–1768. [Google Scholar] [CrossRef]
  33. Commito, J.A.; Dankers, N. Dynamics of spatial and temporal complexity in European and North American soft-bottom mussel beds. In Ecological Comparisons of Sedimentary Shores; Reise, K., Ed.; Springer: Heidelberg, Germany, 2001; pp. 39–59. ISBN 3-540-41254-9. [Google Scholar]
  34. Widdows, J.; Brinsley, M. Impact of biotic and abiotic processes on sediment dynamics and the consequences to the structure and functioning of the intertidal zone. J. Sea Res. 2002, 48, 143–156. [Google Scholar] [CrossRef]
  35. Widdows, J.; Pope, N.D.; Brinsley, M.D.; Gascoigne, J.; Kaiser, M.J. Influence of self-organised structures on near-bed hydrodynamics and sediment dynamics within a mussel (Mytilus edulis) bed in the Menai Strait. J. Exp. Mar. Biol. Ecol. 2009, 379, 92–100. [Google Scholar] [CrossRef]
  36. Gutiérrez, J.L.; Jones, C.G.; Strayer, D.L.; Iribarne, O.O. Mollusks as ecosystem engineers: The role of shell production in aquatic habitats. Oikos 2003, 101, 79–90. [Google Scholar] [CrossRef]
  37. Emerson, C.W.; Grant, J. The control of soft-shell clam (Mya arenaria) recruitment on intertidal sandflats by bedload sediment transport. Limnol. Oceanogr. 1991, 36, 1288–1300. [Google Scholar] [CrossRef]
  38. Armonies, W. Drifting meio-and microbenthic invertebrates on tidal flats in Königshafen: A review. Helgol. Meeresunters. 1994, 48, 299–320. [Google Scholar] [CrossRef]
  39. Commito, J.A.; Thrush, S.A.; Pridmore, R.D.; Hewitt, J.E.; Cummings, V.J. Dispersal dynamics in a wind-driven benthic system. Limnol. Oceanogr. 1995, 40, 1513–1518. [Google Scholar] [CrossRef] [Green Version]
  40. Turner, S.J.; Grant, J.; Pridmore, R.D.; Hewitt, J.E.; Wilkinson, M.R.; Hume, T.M.; Morrisey, D.J. Bedload and water-column transport and colonization processes by post-settlement benthic macrofauna: Does infaunal density matter? J. Exp. Mar. Biol. Ecol. 1997, 216, 51–75. [Google Scholar] [CrossRef]
  41. Lundquist, C.J.; Thrush, S.F.; Hewitt, J.E.; Halliday, J.; MacDonald, I.; Cummings, V.J. Spatial variability in recolonisation potential: Influence of organism behaviour and hydrodynamics on the distribution of macrofaunal colonists. Mar. Ecol. Prog. Ser. 2006, 324, 67–81. [Google Scholar] [CrossRef]
  42. Hunt, H.L.; Maltais, M.-J.; Fugate, D.C.; Chant, R.J. Spatial and temporal variability in juvenile bivalve dispersal: Effects of sediment transport and flow regime. Mar. Ecol. Prog. Ser. 2007, 352, 145–159. [Google Scholar] [CrossRef]
  43. Jennings, L.B.; Hunt, H.L. Distances of dispersal of juvenile bivalves (Mya arenaria (Linnaeus), Mercenaria mercenaria (Linnaeus), Gemma gemma (Totten)). J. Exp. Mar. Biol. Ecol. 2009, 376, 76–84. [Google Scholar] [CrossRef]
  44. Valanko, S.; Norkko, A.; Norkko, J. Rates of post-larval bedload dispersal in a non-tidal soft-sediment system. Mar. Ecol. Prog. Ser. 2010, 416, 145–163. [Google Scholar] [CrossRef] [Green Version]
  45. Valanko, S.; Norkko, A.; Norkko, J. Strategies of post-larval dispersal in non-tidal soft-sediment communities. J. Exp. Mar. Biol. Ecol. 2010, 384, 51–60. [Google Scholar] [CrossRef]
  46. Pacheco, A.S.; Uribe, R.A.; Thiel, M.; Oliva, M.E.; Riascos, J.M. Dispersal of post-larval macrobenthos in subtidal sedimentary habitats: Roles of vertical diel migration, water column, bedload transport and biological traits’ expression. J. Sea Res. 2013, 77, 79–92. [Google Scholar] [CrossRef]
  47. Green, M.O.; Coco, G. Review of wave-driven sediment resuspension and transport in estuaries. Rev. Geophys. 2014, 52, 77–117. [Google Scholar] [CrossRef] [Green Version]
  48. Levin, S.A.; Pacala, S.W. Theories of simplification and scaling of spatially distributed processes. In Spatial Ecology: The Role of Space in Population Dynamics and Interspecific Interactions; Tilman, D., Kareiva, P., Eds.; Princeton Univ. Press: Princeton, NJ, USA, 1997; pp. 271–295. ISBN 9780691016528. [Google Scholar]
  49. Hanski, I. Metapopulation Ecology; Oxford Univ. Press: Oxford, UK, 1999; ISBN 9780198540656. [Google Scholar]
  50. Hiebeler, D. Populations on fragmented landscapes with spatially structured heterogeneities: Landscape generation and local dispersal. Ecology 2000, 81, 1629–1641. [Google Scholar] [CrossRef]
  51. Thrush, S.F.; Hewitt, J.; Cummings, V.; Green, M.; Funnell, G.A.; Wilkinson, M.R. The generality of field experiments: Interactions between local and broad-scale processes. Ecology 2000, 81, 399–415. [Google Scholar] [CrossRef]
  52. Pilditch, C.A.; Valanko, S.; Norkko, J.; Norkko, A. Post-settlement dispersal: The neglected link in maintenance of soft-sediment biodiversity. Biol. Lett. 2015, 11, 20140795. [Google Scholar] [CrossRef]
  53. Palmer, M.A. Dispersal of marine meiofauna: A review and conceptual model explaining passive transport and active emergence with implications for recruitment. Mar. Ecol. Prog. Ser. 1988, 48, 81–91. [Google Scholar] [CrossRef]
  54. Fegley, S.R. A comparison of meiofaunal settlement onto the sediment surface and recolonization of defaunated sandy sediment. J. Exp. Mar. Biol. Ecol. 1988, 123, 97–113. [Google Scholar] [CrossRef]
  55. DePatra, K.D.; Levin, L.A. Evidence of the passive deposition of meiofauna into fiddler crab burrows. J. Exp. Mar. Biol. Ecol. 1989, 125, 173–192. [Google Scholar] [CrossRef]
  56. Fleeger, J.W.; Yund, P.O.; Sun, B. Active and passive processes associated with initial settlement and postsettlement dispersal of suspended meiobenthic copepods. J. Mar. Res. 1995, 53, 609–645. [Google Scholar] [CrossRef]
  57. Commito, J.A.; Tita, G. Differential dispersal rates in an intertidal meiofauna assemblage. J. Exp. Mar. Biol. Ecol. 2002, 268, 237–256. [Google Scholar] [CrossRef]
  58. Hunt, H.L. Effects of sediment source and flow regime on clam and sediment transport. Mar. Ecol. Prog. Ser. 2005, 296, 143–153. [Google Scholar] [CrossRef] [Green Version]
  59. Hurlbert, S.H. Pseudoreplication and the design of ecological field experiments. Ecol. Monogr. 1984, 54, 187–211. [Google Scholar] [CrossRef]
  60. Pollock, L.W. A Practical Guide to the Marine Animals of Northeastern North America; Rutgers University Press: New Brunswick, NJ, USA, 1998; ISBN 0-8135-2398-2. [Google Scholar]
  61. Dean, W.E. Determination of carbonate and organic matter in calcareous sediments and sedimentary rocks by loss on ignition: A comparison with other methods. J. Sediment. Petrol. 1974, 44, 242–248. [Google Scholar]
  62. Commito, J.A.; Currier, C.A.; Kane, L.R.; Reinsel, K.A.; Ulm, I.M. Dispersal dynamics of the bivalve Gemma gemma in a patchy environment. Ecol. Monogr. 1995, 65, 1–20. [Google Scholar] [CrossRef]
  63. Underwood, A.J. Experiments in Ecology: Their Logical Design and Interpretation Using Analysis of Variance; Cambridge University: Cambridge, UK, 1997; ISBN 0-521-55329-6. [Google Scholar]
  64. StatSoft 6.1. STATISTICA Data Analysis Software System; StatSoft, Inc.: Tulsa, OK, USA, 1994. [Google Scholar]
  65. Anderson, M.J. A new method for non-parametric multivariate analysis of variance. Aust. J. Ecol. 2001, 26, 32–46. [Google Scholar]
  66. McArdle, B.H.; Anderson, M.J. Fitting multivariate models to community data: A comment on distance-based redundancy analysis. Ecology 2001, 82, 290–297. [Google Scholar] [CrossRef]
  67. R Core Development Team. R: A Language and Environment for Statistical Computing 2017; R Foundation for Statistical Computing: Vienna, Austria, 2015; Available online: http://www.R-project.org/ (accessed on 12 June 2015).
  68. Clarke, K.R.; Warwick, R.M. Changes in Marine Communities: An Approach to Statistical Analysis and Interpretation, 2nd ed.; PRIMER-E: Plymouth, UK, 2001. [Google Scholar]
  69. Zar, J.H. Biostatistical Analysis, 4th ed.; Prentice Hall: Upper Saddle River, NJ, USA, 1999; ISBN 0-13-081542-X. [Google Scholar]
  70. Hily, C. Is the activity of benthic suspension feeders a factor controlling water quality in the Bay of Brest? Mar. Ecol. Prog. Ser. 1991, 69, 179–188. [Google Scholar] [CrossRef]
  71. Kraeuter, J.N.; Kennish, M.J.; Dobarro, J.; Fegley, S.R.; Flimlin, G.E. Rehabilitation of the northern quahog (hard clam) (Mercenaria mercenaria) habitats by shelling—11 years in Barnegat Bay, New Jersey. J. Shellfish Res. 2003, 22, 61–67. [Google Scholar]
  72. Guay, M.; Himmelman, J.H. Would adding scallop shells (Chlamys islandica) to the sea bottom enhance recruitment of commercial species? J. Exp. Mar. Biol. Ecol. 2004, 312, 299–317. [Google Scholar] [CrossRef]
  73. Ribeiro, P.D.; Iribarne, O.O.; Daleo, P. The relative importance of substratum characteristics and recruitment in determining the spatial distribution of the fiddler crab Uca uruguayensis Nobili. J. Exp. Mar. Biol. Ecol. 2005, 314, 99–111. [Google Scholar] [CrossRef]
  74. Rodney, W.S.; Paynter, K.T. Comparisons of macrofaunal assemblages on restored and non-restored oyster reefs in mesohaline regions of Chesapeake Bay in Maryland. J. Exp. Mar. Biol. Ecol. 2006, 35, 39–51. [Google Scholar] [CrossRef]
  75. Summerhayes, S.A.; Bishop, M.J.; Leigh, A.; Kelaher, B.P. Effects of oyster death and shell disarticulation on associated communities of epibiota. J. Exp. Mar. Biol. Ecol. 2009, 379, 60–67. [Google Scholar] [CrossRef]
  76. Wilding, T.A.; Nickell, T.D. Changes in benthos associated with mussel (Mytilus edulis L.) farms on the West-Coast of Scotland. PLoS ONE 2013, 8, e68313. [Google Scholar] [CrossRef] [PubMed]
  77. Hubbard, W.A. Benthic studies in upper Buzzards Bay, Massachusetts: 2011–12 as compared to 1955. Mar. Ecol. 2016, 37, 532–542. [Google Scholar] [CrossRef]
  78. Tomatsuri, M.; Kon, K. Effects of dead oyster shells as a habitat for the benthic faunal community along rocky shore regions. Hydrobiologia 2017, 790, 25–232. [Google Scholar] [CrossRef]
  79. Gutiérrez, J.L.; Iribarne, O.O. Role of Holocene beds of the stout razor clam Tagelus plebeius in structuring present benthic communities. Mar. Ecol. Prog. Ser. 1999, 185, 213–228. [Google Scholar] [CrossRef]
  80. Bomkamp, R.E.; Page, H.M.; Dugan, J.E. Role of food subsidies and habitat structure in influencing benthic communities of shell mounds at sites of existing and former offshore oil platforms. Mar. Biol. 2004, 146, 201–211. [Google Scholar] [CrossRef]
  81. Hewitt, J.E.; Thrush, S.F.; Halliday, J.; Duffy, C. The importance of small-scale habitat structure for maintaining beta diversity. Ecology 2005, 86, 1619–1626. [Google Scholar] [CrossRef]
  82. Nicastro, A.; Bishop, M.; Kelaher, B.P.; Benedetti-Cecchi, L. Export of non-native gastropod shells to a coastal lagoon: Alteration of habitat structure has negligible effects on infauna. J. Exp. Mar. Biol. Ecol. 2009, 374, 31–36. [Google Scholar] [CrossRef]
  83. Mann, R.; Southworth, M.; Fisher, R.J.; Wesson, J.A.; Erskine, A.J.; Leggett, T. Oyster planting protocols to deter losses to cownose ray predation. J. Shellfish Res. 2016, 35, 127–136. [Google Scholar] [CrossRef]
  84. Folkard, A.M.; Gascoigne, J. Hydrodynamics of discontinuous mussel beds: Laboratory flume simulations. J. Sea Res. 2009, 62, 250–257. [Google Scholar] [CrossRef]
  85. wa Kangeri, A.K.; Jansen, J.M.; Joppe, D.J.; Dankers, N.M.J.A. In situ investigation of the effects of current velocity on sedimentary mussel bed stability. J. Exp. Mar. Biol. Ecol. 2016, 485, 65–72. [Google Scholar] [CrossRef]
  86. Singer, J.K.; Anderson, J.B. Use of total grain-size distributions to define bed erosion and transport for poorly sorted sediment undergoing simulated bioturbation. Mar. Geol. 1984, 57, 335–359. [Google Scholar] [CrossRef]
  87. Chen, P.; Fu, Z.; Lim, A.; Rodrigues, B. Two-Dimensional Packing for Irregular Shaped Objects. In Proceedings of the 36th Hawaii International Conference on System Sciences, Big Island, HI, USA, USA, 6–9 January 2003; Institute of Electrical and Electronics Engineers: New York, NY, USA, 2003. [Google Scholar] [CrossRef]
  88. Ragnarsson, S.A.; Raffaelli, D. Effects of Mytilus edulis L. on the invertebrate fauna of sediments. J. Exp. Mar. Biol. Ecol. 1999, 241, 31–43. [Google Scholar] [CrossRef]
  89. Albrecht, A.; Reise, K. Effects of Fucus vesiculosus covering intertidal mussel beds in the Wadden Sea. Helgol. Meeresunters. 1994, 48, 243–256. [Google Scholar] [CrossRef]
  90. Albrecht, S.H. Soft bottom versus hard rock: Community ecology of macroalgae on intertidal mussel beds in the Wadden Sea. J. Exp. Mar. Biol. Ecol. 1998, 229, 85–109. [Google Scholar] [CrossRef]
  91. Crooks, J.A. Habitat alteration and community-level effects of an exotic mussel, Musculista senhousia. Mar. Ecol. Prog. Ser. 1998, 162, 137–152. [Google Scholar] [CrossRef]
  92. Crooks, J.A.; Khim, H.S. Architectural vs. biological effects of a habitat-altering, exotic mussel, Musculista senhousia. J. Exp. Mar. Biol. Ecol. 1999, 240, 53–75. [Google Scholar] [CrossRef]
  93. Norkko, A.; Hewitt, J.E.; Thrush, S.F.; Funnell, G.A. Benthic–pelagic coupling and suspension-feeding bivalves: Linking site-specific sediment flux and biodeposition to benthic community structure. Limnol. Oceanogr. 2001, 46, 2067–2072. [Google Scholar] [CrossRef]
  94. Hewitt, J.E.; Thrush, S.F.; Legendre, P.; Cummings, V.J.; Norkko, A. Integrating heterogeneity across spatial scales: Interactions between Atrina zelandica and benthic macrofauna. Mar. Ecol. Prog. Ser. 2002, 239, 115–128. [Google Scholar] [CrossRef]
  95. Donadi, S.; van der Heide, T.; van der Zee, E.M.; Eklöf, J.S.; van de Koppel, J.; Weerman, E.J.; Piersma, T.; Olff, H.; Eriksson, B.K. Cross-habitat interactions among bivalve species control community structure on intertidal flats. Ecology 2013, 94, 489–498. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Sun, B.; Fleeger, J.W. Field experiments on the colonization of meiofauna into sediment depressions. Mar. Ecol. Prog. Ser. 1994, 110, 167–175. [Google Scholar] [CrossRef]
  97. Zühlke, R.; Reise, K. Response of macrofauna to drifting tidal sediments. Helgol. Meeresunters. 1994, 48, 277–289. [Google Scholar] [CrossRef] [Green Version]
  98. Pershing, A.J.; Alexander, M.A.; Hernandez, C.M.; Kerr, L.A.; Le Bris, A.; Mills, K.E.; Nye, J.A.; Record, N.R.; Scannell, H.A.; Scott, J.D.; et al. Slow adaptation in the face of rapid warming leads to collapse of the Gulf of Maine cod fishery. Science 2015, 350, 809–812. [Google Scholar] [CrossRef]
  99. Sallenger, A.H.; Doran, K.S.; Howd, P.A. Hotspot of accelerated sea-level rise on the Atlantic coast of North America. Nat. Clim. Chang. 2012, 2, 884–888. [Google Scholar] [CrossRef] [Green Version]
  100. Donnelly, J.P.; Bertness, M.D. Rapid shoreward encroachment of salt marsh cordgrass in response to accelerated sea-level rise. Proc. Natl. Acad. Sci. USA 2001, 98, 14218–14223. [Google Scholar] [CrossRef] [Green Version]
  101. van der Wegen, M.; Jaffe, B.; Foxgrover, A.; Roelvink, D. Mudflat morphodynamics and the impact of sea level rise in south San Francisco Bay. Estuar. Coasts 2017, 40, 37–49. [Google Scholar] [CrossRef]
  102. Beaugrand, G.; Edwards, M.; Brander, K.; Luczak, C.; Ibanez, F. Causes and projections of abrupt climate–driven ecosystem shifts in the North Atlantic. Ecol. Lett. 2008, 11, 157–1168. [Google Scholar] [CrossRef]
  103. Edwards, M.; Beaugrand, G.; Helaouët, P.; Alheit, J.; Coombs, S. Marine ecosystem response to the Atlantic multidecadal oscillation. PLoS ONE 2013, 8, e57212. [Google Scholar] [CrossRef]
  104. Greene, C.H.; Meyer-Gutbrod, E.; Monger, B.C.; McGarry, L.P.; Pershing, A.J.; Belkin, I.M.; Fratantoni, P.S.; Mountain, D.G.; Pickart, R.S.; Proshutinsky, A.; et al. Remote climate forcing of decadal-scale regime shifts in Northwest Atlantic shelf ecosystems. Limnol. Oceanogr. 2013, 58, 803–816. [Google Scholar] [CrossRef] [Green Version]
  105. Reise, K.; Buschbaum, C.; Büttger, H.; Wegner, K.M. Invading oysters and native mussels: From hostile takeover to compatible bedfellows. Ecosphere 2017, 8, e01949. [Google Scholar] [CrossRef]
  106. Kidwell, S.M.; Jablonski, D. Taphonomic feedback: Ecological consequences of shell accumulation. In Biotic Interactions in Recent and Fossil Benthic Communities; Tevesz, M.J.S., McCall, P.L., Eds.; Plenum Press: New York, New York, USA, 1983; pp. 195–248. ISBN 978-1-4757-0742-7. [Google Scholar]
  107. Reise, K. Sediment mediated species interactions in coastal waters. J. Sea Res. 2002, 48, 127–141. [Google Scholar] [CrossRef]
  108. Powell, E.N.; Staff, G.M.; Callender, W.R.; Ashton-Alcox, K.A.; Brett, C.E.; Parsons-Hubbard, K.M.; Walker, S.E.; Raymond, A. Taphonomic degradation of molluscan remains during thirteen years on the continental shelf and slope of the northwestern Gulf of Mexico. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2011, 312, 209–232. [Google Scholar] [CrossRef]
Diversity 11 00011 g001 550
Figure 1. (A) Collapsed mussel beds in eastern Maine, USA. Shell hash is the primary component of the surface of the beds pictured here. (B) The four bottom cover types commonly observed at robust and collapsed mussel beds.
Figure 1. (A) Collapsed mussel beds in eastern Maine, USA. Shell hash is the primary component of the surface of the beds pictured here. (B) The four bottom cover types commonly observed at robust and collapsed mussel beds.
Diversity 11 00011 g001
Diversity 11 00011 g002 550
Figure 2. (A) Study site at Carrying Place Cove, Harrington, Maine, USA. Longest fetch is from the east, indicated by red arrow. (B) Wind condition index (Day 1 > Day 3 > Day 2) created by combining wind velocity and fetch values when the bottom was covered with water during the three days of the investigation. Map data copyright 2015 Google.
Figure 2. (A) Study site at Carrying Place Cove, Harrington, Maine, USA. Longest fetch is from the east, indicated by red arrow. (B) Wind condition index (Day 1 > Day 3 > Day 2) created by combining wind velocity and fetch values when the bottom was covered with water during the three days of the investigation. Map data copyright 2015 Google.
Diversity 11 00011 g002
Diversity 11 00011 g003 550
Figure 3. (A) Mass values of silt-clay, fine sand, and coarse material in traps. For clarity, error bars are presented only for total sediment mass. (B) Percent organic matter. Values are mean ± 1 SE. Live = live mussels, Whole = whole shells, Frag = fragmented shells, Bare = bare sediment. 1 = Day 1, 2 = Day 2, 3 = Day 3.
Figure 3. (A) Mass values of silt-clay, fine sand, and coarse material in traps. For clarity, error bars are presented only for total sediment mass. (B) Percent organic matter. Values are mean ± 1 SE. Live = live mussels, Whole = whole shells, Frag = fragmented shells, Bare = bare sediment. 1 = Day 1, 2 = Day 2, 3 = Day 3.
Diversity 11 00011 g003
Diversity 11 00011 g004 550
Figure 4. (A) Meiofauna values for number of taxa and H’ in traps. (B) Meiofauna absolute dispersal rates (ADR). (C) Meiofauna relative dispersal rates (RDR). (D) Meiofauna bulk dispersal rates (BDR). Values are mean ± 1 SE. Live = Live mussels, Whole = Whole shells, Frag = Fragmented shells, Bare = Bare sediment. 1 = Day 1, 2 = Day 2, 3 = Day 3.
Figure 4. (A) Meiofauna values for number of taxa and H’ in traps. (B) Meiofauna absolute dispersal rates (ADR). (C) Meiofauna relative dispersal rates (RDR). (D) Meiofauna bulk dispersal rates (BDR). Values are mean ± 1 SE. Live = Live mussels, Whole = Whole shells, Frag = Fragmented shells, Bare = Bare sediment. 1 = Day 1, 2 = Day 2, 3 = Day 3.
Diversity 11 00011 g004
Diversity 11 00011 g005 550
Figure 5. Nonmetric multidimensional scaling (nMDS) plots for traps: (A) meiofauna, (B) macrofauna.
Figure 5. Nonmetric multidimensional scaling (nMDS) plots for traps: (A) meiofauna, (B) macrofauna.
Diversity 11 00011 g005
Diversity 11 00011 g006 550
Figure 6. (A) Macrofauna values for number of taxa and H’ in traps. (B) Macrofauna absolute dispersal rates (ADRs). (C) Meiofauna bulk dispersal rates (BDRs). Values are mean ± 1 SE. Live = live mussels, Whole = whole shells, Frag = fragmented shells, Bare = bare sediment. 1 = Day 1, 2 = Day 2, 3 = Day 3.
Figure 6. (A) Macrofauna values for number of taxa and H’ in traps. (B) Macrofauna absolute dispersal rates (ADRs). (C) Meiofauna bulk dispersal rates (BDRs). Values are mean ± 1 SE. Live = live mussels, Whole = whole shells, Frag = fragmented shells, Bare = bare sediment. 1 = Day 1, 2 = Day 2, 3 = Day 3.
Diversity 11 00011 g006
Table
Table 1. ANOVA results for sediment in traps: mass values and percent total organic material. Transformations indicated where applied. Significant differences (P < 0.05) are in bold and Student-Newman-Keuls (SNK) results are presented. NAH = no alternative hypothesis; the two groups with the largest and smallest means within the factor were significantly different, and no further resolution among the groups in that factor could be determined. L = live mussels, F = fragmented shells. D1 = Day 1, D2 = Day 2, D3 = Day 3.
Table 1. ANOVA results for sediment in traps: mass values and percent total organic material. Transformations indicated where applied. Significant differences (P < 0.05) are in bold and Student-Newman-Keuls (SNK) results are presented. NAH = no alternative hypothesis; the two groups with the largest and smallest means within the factor were significantly different, and no further resolution among the groups in that factor could be determined. L = live mussels, F = fragmented shells. D1 = Day 1, D2 = Day 2, D3 = Day 3.
Sediment
Coarse (g)Fine (g)Silt-Clay (g)
SourceDFFPFPFP
Day22.180.1947.130.0262.470.165
Cover32.050.1340.700.5605.260.006
Day × Cover60.700.6520.920.5001.350.275
Residual24
Total35
Transformation Ln(x + 1) Ln(x + 1)
SNK Day D1 > D2 = D3
SNK Cover NAH (L > F)
Total (g)Organic Matter (%)
SourceDFFPFP
Day27.010.02711.340.0003
Cover31.210.3280.510.680
Day × Cover60.680.6641.020.438
Residual24
Total35
Transformation Ln(x + 1)Ln(x + 1)
SNK Day D1 > D3 = D2D3 > D1 > D2
SNK Cover
Table
Table 2. ANOVA results for meiofauna in traps: number of species, H’, and absolute dispersal rates (ADR). Transformations indicated where applied. Significant differences (P < 0.05) are in bold and SNK results are presented. NAH = no alternative hypothesis; the two groups with the largest and smallest means within the factor were significantly different, and no further resolution among the groups in that factor could be determined. L = live mussels, W = whole shells, F = fragmented shells, B = bare sediment. D1 = Day 1, D2 = Day 2, D3 = Day 3.
Table 2. ANOVA results for meiofauna in traps: number of species, H’, and absolute dispersal rates (ADR). Transformations indicated where applied. Significant differences (P < 0.05) are in bold and SNK results are presented. NAH = no alternative hypothesis; the two groups with the largest and smallest means within the factor were significantly different, and no further resolution among the groups in that factor could be determined. L = live mussels, W = whole shells, F = fragmented shells, B = bare sediment. D1 = Day 1, D2 = Day 2, D3 = Day 3.
Meiofauna
No. of TaxaH’Total Assemblage
SourceDFFPFPFP
Day20.720.4954.170.0286.360.006
Cover31.200.3335.700.0044.340.014
Day × Cover60.540.7722.320.0660.290.936
Residual24
Total35
Transformation Ln(x + 1)
SNK Day D1 > D2 = D3D1 = D3 > D2
SNK Cover NAH (B > L)W = L = B > F
CopepodsNematodesForaminiferans
SourceDFFPFPFP
Day23.720.03910.990.00043.570.044
Cover33.670.0264.040.0194.390.013
Day × Cover60.520.7870.380.8861.470.229
Residual24
Total35
Transformation Ln(x + 1)Ln(x + 1)Ln(x + 1)
SNK Day D1 = D3 > D2D1 > D3 = D2NAH (D1 > D3)
SNK Cover W = L = B > FNAH (W > F)NAH (B > F)
Table
Table 3. PERMANOVA results for meiofauna and macrofauna in traps. Significant differences (P < 0.05) are in bold. For parameters with significant differences among Days or Cover Types, the mean Bray-Curtis dissimilarity values (%) within and among Days and Cover Types are presented. L = Live mussels, W = Whole shells, F = Fragmented shells, B = Bare sediment. D1 = Day 1, D2 = Day 2, D3 = Day 3.
Table 3. PERMANOVA results for meiofauna and macrofauna in traps. Significant differences (P < 0.05) are in bold. For parameters with significant differences among Days or Cover Types, the mean Bray-Curtis dissimilarity values (%) within and among Days and Cover Types are presented. L = Live mussels, W = Whole shells, F = Fragmented shells, B = Bare sediment. D1 = Day 1, D2 = Day 2, D3 = Day 3.
MeiofaunaMacrofauna
SourceDFFPFP
Day23.690.0061.340.179
Cover33.030.0111.350.153
Day × Cover60.560.9341.410.076
Residual24
Total35
Within Days:
D140.41
D248.92
D341.65
Among Days:
D1 vs. D253.39
D1 vs. D343.87
D2 vs. D347.84
Within Cover Types:
L45.32
W47.15
F49.51
B30.04
Among Cover Types:
B vs. L38.69
B vs. W43.06
B vs. F52.83
L vs. W48.8
L vs. F56.21
W vs. F48.66
Table
Table 4. ANOVA results for meiofauna in traps: relative dispersal rates (RDR). Transformations indicated where applied. Significant differences (P < 0.05) are in bold and SNK results are presented. NAH = no alternative hypothesis; the 2 groups with the largest and smallest means within the factor were significantly different, and no further resolution among the groups in that factor could be determined. L = Live mussels, W = Whole shells, F = Fragmented shells, B = Bare sediment. D1 = Day 1, D2 = Day 2, D3 = Day 3.
Table 4. ANOVA results for meiofauna in traps: relative dispersal rates (RDR). Transformations indicated where applied. Significant differences (P < 0.05) are in bold and SNK results are presented. NAH = no alternative hypothesis; the 2 groups with the largest and smallest means within the factor were significantly different, and no further resolution among the groups in that factor could be determined. L = Live mussels, W = Whole shells, F = Fragmented shells, B = Bare sediment. D1 = Day 1, D2 = Day 2, D3 = Day 3.
Meiofauna
Total AssemblageCopepods
SourceDFFPFP
Day22.320.1201.400.267
Cover35.210.0073.290.038
Day × Cover61.020.4370.960.474
Residual24
Total35
Transformation Ln(x + 1)Ln(x + 1)
SNK Day
SNK Cover W > L = F = BW > F = L = B
NematodesForaminiferans
SourceDFFPFP
Day21.750.1964.060.030
Cover33.000.0513.870.022
Day × Cover60.570.7520.890.521
Residual24
Total35
Transformation
SNK Day D1 > D3 = D2
SNK Cover NAH (W > F)
Table
Table 5. ANOVA results for meiofauna in traps: bulk dispersal rates (BDRs). Transformations indicated where applied. Significant differences (P < 0.05) are in bold. For parameters with significant differences, the SNK results are presented. NAH = no alternative hypothesis; the two groups with the largest and smallest means within the factor were significantly different, and no further resolution among the groups in that factor could be determined. F = fragmented shells, B = bare sediment. D1 = Day 1, D2 = Day 2, D3 = Day 3.
Table 5. ANOVA results for meiofauna in traps: bulk dispersal rates (BDRs). Transformations indicated where applied. Significant differences (P < 0.05) are in bold. For parameters with significant differences, the SNK results are presented. NAH = no alternative hypothesis; the two groups with the largest and smallest means within the factor were significantly different, and no further resolution among the groups in that factor could be determined. F = fragmented shells, B = bare sediment. D1 = Day 1, D2 = Day 2, D3 = Day 3.
Meiofauna
Total AssemblageCopepods
SourceDFFPFP
Day20.980.3893.770.038
Cover32.220.1122.200.114
Day × Cover60.590.7340.810.571
Residual24
Total35
Transformation
SNK Day D2 = D3 > D1
SNK Cover
NematodesForaminiferans
SourceDFFPFP
Day20.620.5440.920.412
Cover33.610.0281.400.268
Day × Cover60.290.9381.020.438
Residual24
Total35
Transformation Ln(x + 1)Ln(x + 1)
SNK Day
SNK Cover NAH (B > F)
Table
Table 6. ANOVA results for macrofauna in traps: number of species, H’, and absolute dispersal rates (ADRs). Transformations indicated where applied.
Table 6. ANOVA results for macrofauna in traps: number of species, H’, and absolute dispersal rates (ADRs). Transformations indicated where applied.
Macrofauna
No. of SpeciesH’Total Assemblage
SourceDFFPFPFP
Day21.610.2210.980.3911.570.229
Cover30.810.5020.300.8251.340.284
Day × Cover60.340.9100.240.9580.490.809
Residual24
Total35
Transformation Ln(x + 1)
OligochaetesNon-Oligochaetes
SourceDFFPFP
Day21.960.1631.880.175
Cover31.290.3020.480.698
Day × Cover60.760.6050.240.959
Residual24
Total35
Transformation Ln(x + 1)
Table
Table 7. ANOVA results for macrofauna in traps: bulk dispersal rates (BDRs). Transformation indicated where applied. Significant differences (P < 0.05) are in bold and SNK results are presented. D1 = Day 1, D2 = Day 2, D3 = Day 3.
Table 7. ANOVA results for macrofauna in traps: bulk dispersal rates (BDRs). Transformation indicated where applied. Significant differences (P < 0.05) are in bold and SNK results are presented. D1 = Day 1, D2 = Day 2, D3 = Day 3.
Macrofauna Parameter
NOligochaetesNon-Oligochaetes
SourceDFFPFPFP
Day24.650.0200.940.4054.910.016
Cover32.090.1281.680.1970.620.612
Day × Cover61.080.4030.910.5020.640.699
Residual24
Total35
Transformation Ln(x + 1)
SNK Day D2 > D3 = D1 D2 > D3 = D1

Share and Cite

MDPI and ACS Style

Commito, J.A.; Jones, B.R.; Jones, M.A.; Winders, S.E.; Como, S. After the Fall: Legacy Effects of Biogenic Structure on Wind-Generated Ecosystem Processes Following Mussel Bed Collapse. Diversity 2019, 11, 11. https://doi.org/10.3390/d11010011

AMA Style

Commito JA, Jones BR, Jones MA, Winders SE, Como S. After the Fall: Legacy Effects of Biogenic Structure on Wind-Generated Ecosystem Processes Following Mussel Bed Collapse. Diversity. 2019; 11(1):11. https://doi.org/10.3390/d11010011

Chicago/Turabian Style

Commito, John A., Brittany R. Jones, Mitchell A. Jones, Sondra E. Winders, and Serena Como. 2019. "After the Fall: Legacy Effects of Biogenic Structure on Wind-Generated Ecosystem Processes Following Mussel Bed Collapse" Diversity 11, no. 1: 11. https://doi.org/10.3390/d11010011

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop