Developing performance measures of mangrove wetlands using simulation models of hydrology, nutrient biogeochemistry, and community dynamics

The goal of mangrove restoration projects should be to improve community structure and ecosystem function of degraded coastal landscapes. This requires the ability to forecast how mangrove structure and function will respond to prescribed changes in site conditions including hydrology, topography, and geophysical energies. There are global, regional, and local factors that can explain gradients of regulators (e.g., salinity, sulfides), resources (nutrients, light, water), and hydroperiod (frequency, duration of flooding) that collectively account for stressors that result in diverse patterns of mangrove properties across a variety of environmental settings. Simulation models of hydrology, nutrient biogeochemistry, and vegetation dynamics have been developed to forecast patterns in mangroves in the Florida Coastal Everglades. These models provide insight to mangrove response to specific restoration alternatives, testing causal mechanisms of system degradation. We propose that these models can also assist in selecting performance measures for monitoring programs that evaluate project effectiveness. This selection process in turn improves model development and calibration for forecasting mangrove response to restoration alternatives. Hydrologic performance measures include soil regulators, particularly soil salinity, surface topography of mangrove landscape, and hydroperiod, including both the frequency and duration of flooding. Estuarine performance measures should include salinity of the bay, tidal amplitude, and conditions of fresh water discharge (included in the salinity value). The most important performance measures from the mangrove biogeochemistry model should include soil resources (bulk density, total nitrogen, and phosphorus) and soil accretion. Mangrove ecology performance measures should include forest dimension analysis (transects and/or plots), sapling recruitment, leaf area index, and faunal relationships. Estuarine ecology performance measures should include the habitat function of mangroves, which can be evaluated with growth rate of key species, habitat suitability analysis, isotope abundance of indicator species, and bird census. The list of performance measures can be modified according to the model output that is used to define the scientific goals during the restoration planning process that reflect specific goals of the project.

Dr. Rene Price and Jeremy Stalker (FIU graduate student) collecting a groundwater sample, Taylor Slough

This material is based upon work supported by the National Science Foundation through the Florida Coastal Everglades Long-Term Ecological Research program under Cooperative Agreements #DBI-0620409 and #DEB-9910514. This image is made available for non-commercial or educational use only.

Development of an integrated surface and subsurface model of Everglades National Park

An integrated surface-subsurface hydrological model of Everglades National Park (ENP) was developed using MIKE SHE and MIKE 11 modeling software. The model has a resolution of 400 meters, covers approximately 1050 square miles of ENP, includes 110 miles of drainage canals with a variety of hydraulic structures, and processes hydrological information, such as evapotranspiration, precipitation, groundwater levels, canal discharges and levels, and operational schedules. Calibration was based on time series and probability of exceedance for water levels and discharges in the years 1987 through 1997. Model verification was then completed for the period of 1998 through 2005. Parameter sensitivity in uncertainty analysis showed that the model was most sensitive to the hydraulic conductivity of the regional Surficial Aquifer System, the Manning's roughness coefficient, and the leakage coefficient, which defines the canal-subsurface interaction. The model offers an enhanced predictive capability, compared to other models currently available, to simulate the flow regime in ENP and to forecast the impact of topography, water flows, and modifying operation schedules.

Crayfish assemblage shifts in a large drought-prone wetland: the roles of hydrology and competition

1. Faster growing, larger and/or more aggressive crayfish species are predicted to dominate permanent waterbodies. We tested this prediction using a 9 year dataset for two species of crayfish (Procambarus alleni and Procambarus fallax) co-existing in a sub-tropical flowing slough in southern Florida. Using a series of laboratory and mesocosm experiments we also compared life history traits and performance of the respective species to test mechanisms that could explain dominance shifts in the local crayfish assemblages. 2. Over the 9-year period, P. alleni densities were the greatest in shallower, shorterhydroperiod areas bordering the slough, while P. fallax densities were higher in deeper, longer-hydroperiod central areas. These areas were separated by 0.8–2 km of continuous wetland with no apparent barriers to movement between them. 3. Density of P. fallax was not strongly affected by any measures of hydrological variation, while P. alleni density increased with more severe drought conditions. Following the strongest droughts, P. alleni colonized areas in the centre of the slough where they had been absent or scarce in wetter years. 4. We conducted experiments to compare growth rates, drought tolerance, and competitive dominance of these species. P. alleni survived drought conditions better, had higher growth rates, and was the dominant competitor for space and food. While drought probably limits P. fallax in the drier slough habitats, neither drought sensitivity nor interspecific competition with P. fallax can explain decreases of P. alleni with wetter conditions. 5. Our results indicate that a competition-colonization tradeoff cannot explain the crayfish compositional dynamics in this wetland because P. alleni is both the best competitor and the best at surviving in and colonizing areas with the strongest droughts. Future attention should focus on the potential for selective effects of predators that co-vary with hydrology. 6. The traits (large size, fast growth, competitive dominance) exhibited by P. alleni, which is absent in long-hydroperiod wetlands, are those exhibited by dominant crayfish in permanent lakes and streams containing fish. Although these traits make crayfish less vulnerable to fish in some lakes and streams, life-history models of community structure across permanence gradients suggest the opposite traits should be favoured for co-existence with fish.

Use of tritium and helium to define groundwater flow conditions in Everglades National Park

The concentrations of tritium (3H) and helium isotopes (3He and4He) were used as tracers of groundwater flow in the surficial aquifer system (SAS) beneath Everglades National Park (ENP), south Florida. From ages determined by 3H/3He dating techniques, groundwater within the upper 28 m originated within the last 30 years. Below 28 m, waters originated prior to 30 years before present with evidence of mixing at the interface. Interannual variation of the 3H/3He ages within the upper 28 m was significant throughout the 3 year investigation, corresponding with varying hydrologic conditions. In the region of Taylor Slough Bridge, younger groundwater was consistently detected below older groundwater in the Biscayne Aquifer, suggesting preferential flow to the lower part of the aquifer. An increase in 4He with depth in the SAS indicated that radiogenic 4He produced in the underlying Hawthorn Group migrates into the SAS by diffusion. Higher Δ4He values in brackish groundwaters compared to fresh waters from similar depths suggested a possible enhanced vertical transport of4He in the seawater mixing zone. Groundwater salinity measurements indicated the presence of a wide (6–28 km) seawater mixing zone. Comparison of groundwater levels with surface water levels in this zone indicated the potential for brackish groundwater discharge to the overlying Everglades surface water.

Coastal groundwater discharge – an additional source of phosphorus for the oligotrophic wetlands of the Everglades

In this manuscript we define a new term we call coastal groundwater discharge (CGD), which is related to submarine groundwater discharge (SGD), but occurs when seawater intrudes inland to force brackish groundwater to discharge to the coastal wetlands. A hydrologic and geochemical investigation of both the groundwater and surface water in the southern Everglades was conducted to investigate the occurrence of CGD associated with seawater intrusion. During the wet season, the surface water chemistry remained fresh. Enhanced chloride, sodium, and calcium concentrations, indicative of brackish groundwater discharge, were observed in the surface water during the dry season. Brackish groundwaters of the southern Everglades contain 1–2.3μM concentrations of total phosphorus (TP). These concentrations exceed the expected values predicted by conservative mixing of local fresh groundwater and intruding seawater, which both have TPμM. The additional source of TP may be from seawater sediments or from the aquifer matrix as a result of water–rock interactions (such as carbonate mineral dissolution and ion exchange reactions) induced by mixing fresh groundwater with intruding seawater. We hypothesize that CGD maybe an additional source of phosphorus (a limiting nutrient) to the coastal wetlands of the southern Everglades.

Hydrological Conditions Control P Loading and Aquatic Metabolism in an Oligotrophic, Subtropical Estuary

Using high-resolution measures of aquatic ecosystem metabolism and water quality, we investigated the importance of hydrological inputs of phosphorus (P) on ecosystem dynamics in the oligotrophic, P-limited coastal Everglades. Due to low nutrient status and relatively large inputs of terrestrial organic matter, we hypothesized that the ponds in this region would be strongly net heterotrophic and that pond gross primary production (GPP) and respiration (R) would be the greatest during the “dry,” euhaline estuarine season that coincides with increased P availability. Results indicated that metabolism rates were consistently associated with elevated upstream total phosphorus and salinity concentrations. Pulses in aquatic metabolism rates were coupled to the timing of P supply from groundwater upwelling as well as a potential suite of hydrobiogeochemical interactions. We provide evidence that freshwater discharge has observable impacts on aquatic ecosystem function in the oligotrophic estuaries of the Florida Everglades by controlling the availability of P to the ecosystem. Future water management decisions in South Florida must include the impact of changes in water delivery on downstream estuaries.

A tidal creek water budget: Estimation of groundwater discharge and overland flow using hydrologic modeling in the Southern Everglades

Taylor Slough is one of the natural freshwater contributors to Florida Bay through a network of microtidal creeks crossing the Everglades Mangrove Ecotone Region (EMER). The EMER ecological function is critical since it mediates freshwater and nutrient inputs and controls the water quality in Eastern Florida Bay. Furthermore, this region is vulnerable to changing hydrodynamics and nutrient loadings as a result of upstream freshwater management practices proposed by the Comprehensive Everglades Restoration Program (CERP), currently the largest wetland restoration project in the USA. Despite the hydrological importance of Taylor Slough in the water budget of Florida Bay, there are no fine scale (∼1 km2) hydrodynamic models of this system that can be utilized as a tool to evaluate potential changes in water flow, salinity, and water quality. Taylor River is one of the major creeks draining Taylor Slough freshwater into Florida Bay. We performed a water budget analysis for the Taylor River area, based on long-term hydrologic data (1999–2007) and supplemented by hydrodynamic modeling using a MIKE FLOOD (DHI,http://dhigroup.com/) model to evaluate groundwater and overland water discharges. The seasonal hydrologic characteristics are very distinctive (average Taylor River wet vs. dry season outflow was 6 to 1 during 1999–2006) with a pronounced interannual variability of flow. The water budget shows a net dominance of through flow in the tidal mixing zone, while local precipitation and evapotranspiration play only a secondary role, at least in the wet season. During the dry season, the tidal flood reaches the upstream boundary of the study area during approximately 80 days per year on average. The groundwater field measurements indicate a mostly upwards-oriented leakage, which possibly equals the evapotranspiration term. The model results suggest a high importance of groundwater contribution to the water salinity in the EMER. The model performance is satisfactory during the dry season where surface flow in the area is confined to the Taylor River channel. The model also provided guidance on the importance of capturing the overland flow component, which enters the area as sheet flow during the rainy season. Overall, the modeling approach is suitable to reach better understanding of the water budget in the mangrove region. However, more detailed field data is needed to ascertain model predictions by further calibrating overland flow parameters.

Phosphorus budgets in Everglades wetland ecosystems: the effects of hydrology and nutrient enrichment

The Florida Everglades is a naturally oligotrophic hydroscape that has experienced large changes in ecosystem structure and function as the result of increased anthropogenic phosphorus (P) loading and hydrologic changes. We present whole-ecosystem models of P cycling for Everglades wetlands with differing hydrology and P enrichment with the goal of synthesizing existing information into ecosystem P budgets. Budgets were developed for deeper water oligotrophic wet prairie/slough (‘Slough’), shallower water oligotrophic Cladium jamaicense (‘Cladium’), partially enriched C. jamaicense/Typha spp. mixture (‘Cladium/Typha’), and enriched Typha spp. (‘Typha’) marshes. The majority of ecosystem P was stored in the soil in all four ecosystem types, with the flocculent detrital organic matter (floc) layer at the bottom of the water column storing the next largest proportion of ecosystem P pools. However, most P cycling involved ecosystem components in the water column (periphyton, floc, and consumers) in deeper water, oligotrophic Slough marsh. Fluxes of P associated with macrophytes were more important in the shallower water, oligotrophic Cladium marsh. The two oligotrophic ecosystem types had similar total ecosystem P stocks and cycling rates, and low rates of P cycling associated with soils. Phosphorus flux rates cannot be estimated for ecosystem components residing in the water column in Cladium/Typha or Typha marshes due to insufficient data. Enrichment caused a large increase in the importance of macrophytes to P cycling in Everglades wetlands. The flux of P from soil to the water column, via roots to live aboveground tissues to macrophyte detritus, increased from 0.03 and 0.2 g P m−2 yr−1 in oligotrophic Slough and Cladium marsh, respectively, to 1.1 g P m−2 yr−1 in partially enriched Cladium/Typha, and 1.6 g P m−2 yr−1 in enriched Typha marsh. This macrophyte translocation P flux represents a large source of internal eutrophication to surface waters in P-enriched areas of the Everglades.

Water source utilization and foliar nutrient status differs between upland and flooded plant communities in wetland tree islands

Tree islands in the Everglades wetlands are centers of biodiversity and targets of restoration, yet little is known about the pattern of water source utilization by the constituent woody plant communities: upland hammocks and flooded swamp forests. Two potential water sources exist: (1) entrapped rainwater in the vadose zone of the organic soil (referred to as upland soil water), that becomes enriched in phosphorus, and (2) phosphorus-poor groundwater/surface water (referred to as regional water). Using natural stable isotope abundance as a tracer, we observed that hammock plants used upland soil water in the wet season and shifted to regional water uptake in the dry season, while swamp forest plants used regional water throughout the year. Consistent with the previously observed phosphorus concentrations of the two water sources, hammock plants had a greater annual mean foliar phosphorus concentration over swamp forest plants, thereby supporting the idea that tree island hammocks are islands of high phosphorus concentrations in the oligotrophic Everglades. Foliar nitrogen levels in swamp forest plants were higher than those of hammock plants. Linking water sources with foliar nutrient concentrations can indicate nutrient sources and periods of nutrient uptake, thereby linking hydrology with the nutrient regimes of different plant communities in wetland ecosystems. Our results are consistent with the hypotheses that (1) over long periods, upland tree island communities incrementally increase their nutrient concentration by incorporating marsh nutrients through transpiration seasonally, and (2) small differences in micro-topography in a wetland ecosystem can lead to large differences in water and nutrient cycles.

A Hydrological Budget (2002–2008) for a Large Subtropical Wetland Ecosystem Indicates Marine Groundwater Discharge Accompanies Diminished Freshwater Flow

Water budget parameters are estimated for Shark River Slough (SRS), the main drainage within Everglades National Park (ENP) from 2002 to 2008. Inputs to the water budget include surface water inflows and precipitation while outputs consist of evapotranspiration, discharge to the Gulf of Mexico and seepage losses due to municipal wellfield extraction. The daily change in volume of SRS is equated to the difference between input and outputs yielding a residual term consisting of component errors and net groundwater exchange. Results predict significant net groundwater discharge to the SRS peaking in June and positively correlated with surface water salinity at the mangrove ecotone, lagging by 1 month. Precipitation, the largest input to the SRS, is offset by ET (the largest output); thereby highlighting the importance of increasing fresh water inflows into ENP for maintaining conditions in terrestrial, estuarine, and marine ecosystems of South Florida.

Determining Spatial and Temporal Inputs of Freshwater, Including Submarine Groundwater Discharge, to a Subtropical Estuary Using Geochemical Tracers, Biscayne Bay, South Florida

Geochemical mixing models were used to decipher the dominant source of freshwater (rainfall, canal discharge, or groundwater discharge) to Biscayne Bay, an estuary in south Florida. Discrete samples of precipitation, canal water, groundwater, and bay surface water were collected monthly for 2 years and analyzed for salinity, stable isotopes of oxygen and hydrogen, and Sr2+/Ca2+ concentrations. These geochemical tracers were used in three separate mixing models and then combined to trace the magnitude and timing of the freshwater inputs to the estuary. Fresh groundwater had an isotopic signature (δ 18O = −2.66‰, δD −7.60‰) similar to rainfall (δ 18O = −2.86‰, δD = −4.78‰). Canal water had a heavy isotopic signature (δ 18O = −0.46‰, δD  = −2.48‰) due to evaporation. This made it possible to use stable isotopes of oxygen and hydrogen to separate canal water from precipitation and groundwater as a source of freshwater into the bay. A second model using Sr2+/Ca2+ ratios was developed to discern fresh groundwater inputs from precipitation inputs. Groundwater had a Sr2+/Ca2+ ratio of 0.07, while precipitation had a dissimilar ratio of 0.89. When combined, these models showed a freshwater input ratio of canal/precipitation/groundwater of 37%:53%:10% in the wet season and 40%:55%:5% in the dry season with an error of ±25%. For a bay-wide water budget that includes saltwater and freshwater mixing, fresh groundwater accounts for 1–2% of the total fresh and saline water input.

Relationships Between Hydrology and Soils Describe Vegetation Patterns in Seasonally Flooded Tree Islands of the Southern Everglades, Florida

Tree island ecosystems are important and distinct features of Florida Everglades wetlands. We described the inter-relationships among abiotic factors describing seasonally flooded tree islands and characterized plant–soil relationships in tree islands occurring in a relatively unimpacted area of the Everglades. We used Principal Components Analysis (PCA) to reduce our multi-factor dataset, quantified forest structure and vegetation nutrient dynamics, and related these vegetation parameters to PCA summary variables using linear regression analyses. We found that, of the 21 abiotic parameters used to characterize the ecosystem structure of seasonally flooded tree islands, 13 parameters were significantly correlated with four principal components, and they described 78% of the variance among the study islands. Most variation was described by factors related to soil oxidation and hydrology, exemplifying the sensitivity of tree island structure to hydrologic conditions. PCA summary variables describing tree island structure were related to variability in Chrysobalanus icaco (L.) canopy cover, Ilex cassine (L.) and Salix caroliniana (Michx.) canopy cover, Myrica cerifera (L.) plot frequency, litter turnover, % phosphorus resorption of co-dominant species, and nitrogen nutrient-use efficiency. This study supported findings that vegetation characteristics can be sensitive indicators of variability in tree island ecosystem structure. This study produced valuable, information which was used to recommend ecological targets (i.e. restoration performance measures) for seasonally flooded tree islands in more impacted regions of the Everglades landscape.

Maintaining tree islands in the Florida Everglades: nutrient redistribution is the key

The Florida Everglades is an oligotrophic wetland system with tree islands as one of its most prominent landscape features. Total soil phosphorus concentrations on tree islands can be 6 to 100 times greater than phosphorus levels in the surrounding marshes and sloughs, making tree islands nutrient hotspots. Several mechanisms are believed to redistribute phosphorus to tree islands: subsurface water flows generated by evapotranspiration of trees, higher deposition rates of dry fallout, deposition of guano by birds and other animals, groundwater upwelling, and bedrock mineralization by tree exudates. A conceptual model is proposed, in which the focused redistribution of limiting nutrients, especially phosphorus, onto tree islands controls their maintenance and expansion. Because of increased primary production and peat accretion rates, the redistribution of phosphorus can result in an increase in both tree island elevation and size. Human changes to hydrology have greatly decreased the number and size of tree islands in parts of the Everglades. The proposed model suggests that the preservation of existing tree islands, and ultimately of the Everglades landscape, requires the maintenance of these phosphorus redistribution mechanisms.