Incorporation of uncertainty to improve projections of tidal wetland elevation and carbon accumulation with sea-level rise.

Understanding the rates and patterns of tidal wetland elevation changes relative to sea-level is essential for understanding the extent of potential wetland loss over the coming years. Using an enhanced and more flexible modeling framework of an ecosystem model (WARMER-2), we explored sea-level rise (SLR) impacts on wetland elevations and carbon sequestration rates through 2100 by considering plant community transitions, salinity effects on productivity, and changes in sediment availability. We incorporated local experimental results for plant productivity relative to inundation and salinity into a species transition model, as well as site-level estimates of organic matter decomposition. The revised modeling framework includes an improved calibration scheme that more accurately reconstructs soil profiles and incorporates parameter uncertainty through Monte Carlo simulations. Using WARMER-2, we evaluated elevation change in three tidal wetlands in the San Francisco Bay Estuary, CA, USA along an estuarine tidal and salinity gradient with varying scenarios of SLR, salinization, and changes in sediment availability. We also tested the sensitivity of marsh elevation and carbon accumulation rates to different plant productivity functions. Wetland elevation at all three sites was sensitive to changes in sediment availability, but sites with greater initial elevations or space for upland transgression persisted longer under higher SLR rates than sites at lower elevations. Using a multi-species wetland vegetation transition model for organic matter contribution to accretion, WARMER-2 projected increased elevations relative to sea levels (resilience) and higher rates of carbon accumulation when compared with projections assuming no future change in vegetation with SLR. A threshold analysis revealed that all three wetland sites were likely to eventually transition to an unvegetated state with SLR rates above 7 mm/yr. Our results show the utility in incorporating additional estuary-specific parameters to bolster confidence in model projections. The new WARMER-2 modeling framework is widely applicable to other tidal wetland ecosystems and can assist in teasing apart important drivers of wetland elevation change under SLR.

Applying cumulative effects to strategically advance large-scale ecosystem restoration.

International efforts to restore degraded ecosystems will continue to expand over the coming decades, yet the factors contributing to the effectiveness of long-term restoration across large areas remain largely unexplored. At large scales, outcomes are more complex and synergistic than the additive impacts of individual restoration projects. Here, we propose a cumulative-effects conceptual framework to inform restoration design and implementation and to comprehensively measure ecological outcomes. To evaluate and illustrate this approach, we reviewed long-term restoration in several large coastal and riverine areas across the US: the greater Florida Everglades; Gulf of Mexico coast; lower Columbia River and estuary; Puget Sound; San Francisco Bay and Sacramento-San Joaquin Delta; Missouri River; and northeastern coastal states. Evidence supported eight modes of cumulative effects of interacting restoration projects, which improved outcomes for species and ecosystems at landscape and regional scales. We conclude that cumulative effects, usually measured for ecosystem degradation, are also measurable for ecosystem restoration. The consideration of evidence-based cumulative effects will help managers of large-scale restoration capitalize on positive feedback and reduce countervailing effects.

Can community structure track sea-level rise? Stress and competitive controls in tidal wetlands.

Climate change impacts, such as accelerated sea-level rise, will affect stress gradients, yet impacts on competition/stress tolerance trade-offs and shifts in distributions are unclear. Ecosystems with strong stress gradients, such as estuaries, allow for space-for-time substitutions of stress factors and can give insight into future climate-related shifts in both resource and nonresource stresses. We tested the stress gradient hypothesis and examined the effect of increased inundation stress and biotic interactions on growth and survival of two congeneric wetland sedges, Schoenoplectus acutus and Schoenoplectus americanus. We simulated sea-level rise across existing marsh elevations and those not currently found to reflect potential future sea-level rise conditions in two tidal wetlands differing in salinity. Plants were grown individually and together at five tidal elevations, the lowest simulating an 80-cm increase in sea level, and harvested to assess differences in biomass after one growing season. Inundation time, salinity, sulfides, and redox potential were measured concurrently. As predicted, increasing inundation reduced biomass of the species commonly found at higher marsh elevations, with little effect on the species found along channel margins. The presence of neighbors reduced total biomass of both species, particularly at the highest elevation; facilitation did not occur at any elevation. Contrary to predictions, we documented the competitive superiority of the stress tolerator under increased inundation, which was not predicted by the stress gradient hypothesis. Multifactor manipulation experiments addressing plant response to accelerated climate change are integral to creating a more realistic, valuable, and needed assessment of potential ecosystem response. Our results point to the important and unpredicted synergies between physical stressors, which are predicted to increase in intensity with climate change, and competitive forces on biomass as stresses increase.

Contributions of organic and inorganic matter to sediment volume and accretion in tidal wetlands at steady state.

A mixing model derived from first principles describes the bulk density (BD) of intertidal wetland sediments as a function of loss on ignition (LOI). The model assumes that the bulk volume of sediment equates to the sum of self-packing volumes of organic and mineral components or BD = 1/[LOI/k1 + (1-LOI)/k2], where k1 and k2 are the self-packing densities of the pure organic and inorganic components, respectively. The model explained 78% of the variability in total BD when fitted to 5075 measurements drawn from 33 wetlands distributed around the conterminous United States. The values of k1 and k2 were estimated to be 0.085 ± 0.0007 g cm-3 and 1.99 ± 0.028 g cm-3, respectively. Based on the fitted organic density (k1) and constrained by primary production, the model suggests that the maximum steady state accretion arising from the sequestration of refractory organic matter is ≤ 0.3 cm yr-1. Thus, tidal peatlands are unlikely to indefinitely survive a higher rate of sea-level rise in the absence of a significant source of mineral sediment. Application of k2 to a mineral sediment load typical of East and eastern Gulf Coast estuaries gives a vertical accretion rate from inorganic sediment of 0.2 cm yr-1. Total steady state accretion is the sum of the parts and therefore should not be greater than 0.5 cm yr-1 under the assumptions of the model. Accretion rates could deviate from this value depending on variation in plant productivity, root:shoot ratio, suspended sediment concentration, sediment-capture efficiency, and episodic events.

Modeling tidal marsh distribution with sea-level rise: evaluating the role of vegetation, sediment, and upland habitat in marsh resiliency.

Tidal marshes maintain elevation relative to sea level through accumulation of mineral and organic matter, yet this dynamic accumulation feedback mechanism has not been modeled widely in the context of accelerated sea-level rise. Uncertainties exist about tidal marsh resiliency to accelerated sea-level rise, reduced sediment supply, reduced plant productivity under increased inundation, and limited upland habitat for marsh migration. We examined marsh resiliency under these uncertainties using the Marsh Equilibrium Model, a mechanistic, elevation-based soil cohort model, using a rich data set of plant productivity and physical properties from sites across the estuarine salinity gradient. Four tidal marshes were chosen along this gradient: two islands and two with adjacent uplands. Varying century sea-level rise (52, 100, 165, 180 cm) and suspended sediment concentrations (100%, 50%, and 25% of current concentrations), we simulated marsh accretion across vegetated elevations for 100 years, applying the results to high spatial resolution digital elevation models to quantify potential changes in marsh distributions. At low rates of sea-level rise and mid-high sediment concentrations, all marshes maintained vegetated elevations indicative of mid/high marsh habitat. With century sea-level rise at 100 and 165 cm, marshes shifted to low marsh elevations; mid/high marsh elevations were found only in former uplands. At the highest century sea-level rise and lowest sediment concentrations, the island marshes became dominated by mudflat elevations. Under the same sediment concentrations, low salinity brackish marshes containing highly productive vegetation had slower elevation loss compared to more saline sites with lower productivity. A similar trend was documented when comparing against a marsh accretion model that did not model vegetation feedbacks. Elevation predictions using the Marsh Equilibrium Model highlight the importance of including vegetation responses to sea-level rise. These results also emphasize the importance of adjacent uplands for long-term marsh survival and incorporating such areas in conservation planning efforts.

Evaluating tidal marsh sustainability in the face of sea-level rise: a hybrid modeling approach applied to San Francisco Bay.

BACKGROUND: Tidal marshes will be threatened by increasing rates of sea-level rise (SLR) over the next century. Managers seek guidance on whether existing and restored marshes will be resilient under a range of potential future conditions, and on prioritizing marsh restoration and conservation activities. METHODOLOGY: Building upon established models, we developed a hybrid approach that involves a mechanistic treatment of marsh accretion dynamics and incorporates spatial variation at a scale relevant for conservation and restoration decision-making. We applied this model to San Francisco Bay, using best-available elevation data and estimates of sediment supply and organic matter accumulation developed for 15 Bay subregions. Accretion models were run over 100 years for 70 combinations of starting elevation, mineral sediment, organic matter, and SLR assumptions. Results were applied spatially to evaluate eight Bay-wide climate change scenarios. PRINCIPAL FINDINGS: Model results indicated that under a high rate of SLR (1.65 m/century), short-term restoration of diked subtidal baylands to mid marsh elevations (-0.2 m MHHW) could be achieved over the next century with sediment concentrations greater than 200 mg/L. However, suspended sediment concentrations greater than 300 mg/L would be required for 100-year mid marsh sustainability (i.e., no elevation loss). Organic matter accumulation had minimal impacts on this threshold. Bay-wide projections of marsh habitat area varied substantially, depending primarily on SLR and sediment assumptions. Across all scenarios, however, the model projected a shift in the mix of intertidal habitats, with a loss of high marsh and gains in low marsh and mudflats. CONCLUSIONS/SIGNIFICANCE: Results suggest a bleak prognosis for long-term natural tidal marsh sustainability under a high-SLR scenario. To minimize marsh loss, we recommend conserving adjacent uplands for marsh migration, redistributing dredged sediment to raise elevations, and concentrating restoration efforts in sediment-rich areas. To assist land managers, we developed a web-based decision support tool (www.prbo.org/sfbayslr).

Diversity-function relationships changed in a long-term restoration experiment.

The central tenet of biodiversity-ecosystem function (BEF) theory, that species richness increases function, could motivate restoration practitioners to incorporate a greater number of species into their projects. But it is not yet clear how well BEF theory predicts outcomes of restoration, because it has been developed through tests involving short-run and tightly controlled (e.g., weeded) experiments. Thus, we resampled our 1997 BEF experiment in a restored salt marsh to test for long-term effects of species richness (plantings with 1, 3, and 6 species per 2 x 2 m plot), with multiple ecosystem functions as response variables. Over 11 years, 1- and 6-species assemblages converged on intermediate richness (mean = 3.9 species/ 0.25-m2 plot), and composition changed nonrandomly throughout the site. While three species became rare, the two most productive species became co-dominant. The two dominants controlled and increased shoot biomass, which appeared to decrease species richness. Diversity-function relationships became less positive over 11 years and differed significantly with (a) the species-richness metric (planted vs. measured), and (b) the indicator of function (shoot biomass, height, and canopy layering). The loss of positive relationships between species richness and function in our restored site began soon after we stopped weeding and continued with increasing dominance by productive species. Where species-rich plantings are unlikely to ensure long-term restoration of functions, as in our salt marsh, we recommend dual efforts to establish (1) dominant species that provide high levels of target functions, and (2) subordinate species, which might provide additional functions under current or future conditions.