skip to main content


Title: Metapopulation dynamics of oysters: sources, sinks, and implications for conservation and restoration
Abstract

Metapopulation and source–sink dynamics are increasingly considered within spatially explicit management of wildlife populations, yet the application of these concepts has generally been limited to comparisons of the performance (e.g., demographic rates or dispersal) inside vs. outside protected areas, and at spatial scales that do not encompass an entire metapopulation. In the present study, a spatially explicit, size‐structured matrix model was applied to simulate the dynamics of an Eastern oyster (Crassostrea virginica) metapopulation in the second largest estuary in the United States—the Albemarle‐Pamlico Estuarine System in North Carolina. The model integrated larval dispersal simulations with empirical measures of oyster demographic rates to simulate the dynamics of the entire oyster metapopulation consisting of 646 reefs and five reef types: (1) restored subtidal reefs closed to harvest (i.e., sanctuaries or protected areas;n = 14), (2) restored subtidal reefs open to harvest (n = 53), (3) natural subtidal reefs open to harvest (n = 301), (4) natural intertidal reefs open to harvest (n = 129), and (5) oyster reefs on manmade, hard structures such as seawalls (n = 149). Key findings included (1) an overall stable, yet slightly declining oyster metapopulation, (2) variable reef type‐specific population trajectories, largely dependent on spatiotemporal variation in larval recruitment, (3) a greater relative importance of inter‐reef larval connectivity on metapopulation dynamics than local larval retention processes, and (4) spatiotemporal variation in the source–sink status of reef subpopulations wherein subtidal sanctuaries and reefs located in the northeastern portion of the estuary were frequent sources. From a management perspective, continued protection of oyster sanctuaries is warranted. Sanctuaries represented only 6.2% of the total reef area, however, they harbored 19% (± 2%) of all oysters and produced 25% (± 6%) of all larvae settling within the metapopulation. Additional management priorities should focus on restoration or conservation of subpopulations that serve as frequent source subpopulations (including those with poor demographic rates, but high connectivity potential), and management of harvest from sink subpopulations. The application of a source–sink framework and similar integrated modeling approach could inform management of oysters in other systems, as well as other species that exhibit similar metapopulation characteristics.

 
more » « less
NSF-PAR ID:
10449363
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Ecosphere
Volume:
12
Issue:
7
ISSN:
2150-8925
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Foundation species, such as mangroves, saltmarshes, kelps, seagrasses, and oysters, thrive within suitable environmental envelopes as narrow ribbons along the land–sea margin. Therefore, these habitat‐forming species and resident fauna are sensitive to modified environmental gradients. For oysters, many estuaries impacted by sea‐level rise, channelization, and municipal infrastructure are experiencing saltwater intrusion and water‐quality degradation that may alter reef distributions, functions, and services. To explore decadal‐scale oyster–reef community patterns across a temperate estuary in response to environmental change, we resampled reefs in the Newport River Estuary (NRE) during 2013–2015 that had previously been studied during 1955–1956. We also coalesced historical NRE reef distribution (1880s–2015), salinity (1913–2015), and water‐quality‐driven shellfish closure boundary (1970s–2015) data to document environmental trends that could influence reef ecology and service delivery. Over the last 60–120 years, the entire NRE has shifted toward higher salinities. Consequently, oyster–reef communities have become less distinct across the estuary, manifest by 20%–27% lower species turnover and decreased faunal richness among NRE reefs in the 2010s relative to the 1950s. During the 2010s, NRE oyster–reef communities tended to cluster around a euhaline, intertidal‐reef type more so than during the 1950s. This followed faunal expansions farther up estuary and biological degradation of subtidal reefs as NRE conditions became more marine and favorable for aggressive, reef‐destroying taxa. In addition to these biological shifts, the area of suitable bottom on which subtidal reefs persist (contracting due to up‐estuary intrusion of marine waters) and support human harvest (driven by water quality, eroding from up‐estuary) has decreased by >75% since the natural history of NRE reefs was first explored. This “coastal squeeze” on harvestable subtidal oysters (reduced from a 4.5‐km to a 0.75‐km envelope along the NRE's main axis) will likely have consequences regarding the economic incentives for future oyster conservation, as well as the suite of services delivered by remaining shellfish reefs (e.g., biodiversity maintenance, seafood supply). More broadly, these findings exemplify how “squeeze” may be a pervasive concern for biogenic habitats along terrestrial or marine ecotones during an era of intense global change.

     
    more » « less
  2. Abstract

    Metapopulation models include spatial population dynamics such as dispersion and migration between subpopulations. Integral projection models (IPMs) can include demographic rates as a function of size. Traditionally, metapopulation models do not included detailed populaiton models such as IPMs. In some situations, both local population dynamics (e.g. size‐based survival) and spatial dynamics are important.

    We present a Python package,MetaIPM, which places IPMs into a metapopulation framework, and allow users to readily construct and apply these models that combine local population dynamics within a metapopulation framework.

    MetaIPMincludes an IPM for each subpopulation that is connected to other subpopulations via a metapopulation movement model. These movements can include dispersion, migration or other patterns. The IPM can include for size‐specific demographic rates (e.g. survival, recruitment) as well as management actions, such as length‐based harvest (e.g. gear specific capture sizes, varying slot limits across political boundaries). The model also allows for changes in metapopulation connectivity between locations, such as a fish passage ladders to enhance movement or deterrents to reduce movement. Thus, resource managers can useMetaIPMto compare different management actions such as the harvest gear type (which can be length‐specific) and harvest locations.

    We demonstrate howMetaIPMmay be applied to inform managers seeking to limit the spread of an invasive species in a system with important metapopulation dynamics. Specifically, we compared removal lengths (all length fish versus longer fish only) for an invasive fish population in a fragmented, inland river system.MetaIPMallowed users to compare the importance of harvesting source populations away from the invasion front, as well as species at the invasion front. The model would also allow for future comparisons of different deterrent placement locations in the system.

    Moving beyond our example system, we describe howMetaIPMcan be applied to other species, systems and management approaches. TheMetaIPMpackages includes Jupyter Notebooks documenting the package as well as a second set of JupyterNotebooks showing the application of the package to our example system.

     
    more » « less
  3. Abstract

    Restoration aims to reverse the global declines of foundation species, but it is unclear how project attributes, the physical setting, and antecedent conditions affect restoration success. In coastal seas worldwide, oyster reef restoration is increasing to counter historical habitat destruction and associated declines in fisheries production and biodiversity. Yet, restoration outcomes are highly variable and the factors that enhance oyster production and nekton abundance and diversity on restored reefs are unresolved. To quantify the drivers of oyster restoration success, we used meta‐analysis to synthesize data from 158 restored reefs paired with unstructured habitats along the United States Gulf and Atlantic coasts. The average recovery of oyster production was 65% greater in subtidal (vs. intertidal) zones, 173% greater in polyhaline (vs. mesohaline) environments and increased with tidal range, demonstrating that physical conditions can strongly influence the restoration success of foundation species. Additionally, restoration increased the relative abundance and richness of nektonic fishes and invertebrates over time as reefs aged (at least 8 years post‐construction). Thus, the restoration benefits for provisioning habitat and enhancing biodiversity accrue over time, highlighting that restoration projects need multiple years to maximize ecosystem functions. Furthermore, long‐term monitoring of restored and control sites is needed to assess restoration outcomes and associated drivers. Last, our work reveals data constraints for several potential drivers of restoration outcomes, including reef construction material, reef dimensions, harvest pressure and disease prevalence. More experimental and observational studies are needed to target these factors and measure them with consistent methods across studies. Our findings indicate that the assisted recovery of foundation species yields several enhancements to ecosystem services, but such benefits are mediated by time and environmental conditions.

     
    more » « less
  4. Long‐term monitoring is vital to understanding the efficacy of restoration approaches and how restoration may enhance ecosystem functions. We revisited restored oyster reefs 13 years post‐restoration and quantified the resident and transient fauna that utilize restored reefs in three differing landscape contexts: on mudflats isolated from vegetated habitat, along the edge of salt marsh, and in between seagrass and salt marsh habitat. Differences observed 1–2 years post‐restoration in reef development and associated fauna within reefs restored on mudflats versus adjacent to seagrass/salt marsh and salt marsh‐only habitats persisted more than 10 years post‐restoration. Reefs constructed on open mudflat habitats had the highest densities of oysters and resident invertebrates compared to those in other landscape contexts, although all restored reefs continued to enhance local densities of invertebrate taxa (e.g. bivalves, gastropods, decapods, polychaetes, etc.). Catch rates of juvenile fishes were enhanced on restored reefs relative to controls, but to a lesser extent than directly post‐restoration, potentially because the reefs have grown vertically within the intertidal and out of the preferred inundation regime of small juvenile fishes. Reef presence and landscape setting did not augment the catch rates of piscivorous fishes in passive gill nets, similar to initial findings; however, hook‐and‐line catch rates were greater on restored reefs than non‐reef controls. We conclude that ecosystem functions and associated services provided by restored habitats can vary both spatially and temporally; therefore, a better understanding of how service delivery varies among landscape setting and over time should enhance efforts to model these processes and restoration decision‐making.

     
    more » « less
  5. Restoration of degraded estuarine oyster reefs typically involves deploying recycled oyster shell. In low‐salinity, low‐predation areas of estuaries, high‐volume shell deployments are known to improve flow conditions and thus oyster survival and growth. It is also hypothesized that the physical structure of restored reefs could suppress foraging by oyster predators in high‐salinity, high‐predation zones. That hypothesis is untested. Given limited resources, it is important to determine how much shell is needed for successful restoration and whether there are diminishing returns in shell addition. In Apalachicola Bay, Florida, we manipulated shell volume on an oyster reef to create three 0.4 ha areas of low (no shell addition), moderate (153 m3shell), and high (306 m3shell) habitat structure. We repeated experiments and surveys over 2 years to determine if restoration success increased with habitat structure. Predation on oysters was greater on the non‐shelled area than on the reshelled reefs, but similar between the two reshelled reefs. Oyster larval supply did not differ among the reef areas, but by the end of the experiment, oyster density (per unit area) increased quadratically with habitat structure, plateauing at high levels of structure. Model selection indicated that the most parsimonious explanation for these patterns was that increased habitat structure reduced predation and increased overall recruitment, but that the higher reshelling treatment did not have better outcomes than moderate reshelling. Thus, restoration could be optimized by deploying a moderate amount of shell per unit area.

     
    more » « less