Fish Community Responses to the Combined Effects of Decreased Hydroperiod and Nonnative Fish Invasions in a Karst Wetland: Are Everglades Solution Holes Sinks for Native Fishes?

Short-hydroperiod Everglades wetlands have been disproportionately affected by reductions in freshwater inflows, land conversion and biotic invasions. Severe hydroperiod reductions in these habitats, including the Rocky Glades, coupled with proximity to canals that act as sources of invasions, may limit their ability to support high levels of aquatic production. We examined whether karst solution holes function as dry-down refuges for fishes, providing a source of marsh colonists upon reflooding, by tracking fish abundance, nonnative composition, and survival in solution holes throughout the dry season. We paired field surveys with an in situ nonnative predation experiment that tested the effects of predation by the recent invader, African jewelfish (Hemichromis letourneuxi) on native fishes. Over the 3 years surveyed, a large number of the solution holes dried before the onset of the wet season, while those retaining water had low survivorship and were dominated by nonnatives. In the experiment, mortality of eastern mosquitofish (Gambusia holbrooki) in the presence of African jewelfish was greater than that associated with deteriorating water quality. Under current water management, findings suggest that solution holes are largely sinks for native fishes, given the high frequency of drydown, extensive period of fish residence, and predation by nonnative fishes.

Lattice Boltzmann modeling of fluid flow and solute transport in karst aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Spatial, geomorphological, and seasonal variability of CDOM in estuaries of the Florida Coastal Everglades

This paper demonstrates the usefulness of fluorescence techniques for long-term monitoring and assessment of the dynamics (sources, transport and fate) of chromophoric dissolved organic matter (CDOM) in highly compartmentalized estuarine regions with non-point water sources. Water samples were collected monthly from a total of 73 sampling stations in the Florida Coastal Everglades (FCE) estuaries during 2001 and 2002. Spatial and seasonal variability of CDOM characteristics were investigated for geomorphologically distinct sub-regions within Florida Bay (FB), the Ten Thousand Islands (TTI), and Whitewater Bay (WWB). These variations were observed in both quantity and quality of CDOM. TOC concentrations in the FCE estuaries were generally higher during the wet season (June–October), reflecting high freshwater loadings from the Everglades in TTI, and a high primary productivity of marine biomass in FB. Fluorescence parameters suggested that the CDOM in FB is mainly of marine/microbial origin, while for TTI and WWB a terrestrial origin from Everglades marsh plants and mangroves was evident. Variations in CDOM quality seemed mainly controlled by tidal exchange/mixing of Everglades freshwater with Florida Shelf waters, tidally controlled releases of CDOM from fringe mangroves, primary productivity of marine vegetation in FB and diagenetic processes such as photodegradation (particularly for WWB). The source and dynamics of CDOM in these subtropical estuaries is complex and found to be influenced by many factors including hydrology, geomorphology, vegetation cover, landuse and biogeochemical processes. Simple, easy to measure, high sample throughput fluorescence parameters for surface waters can add valuable information on CDOM dynamics to long-term water quality studies which can not be obtained from quantitative determinations alone.

Chemical characterization of dissolved organic nitrogen in an oligotrophic subtropical coastal ecosystem

Water samples were collected from rivers and estuarine environments within the Florida Coastal Everglades (FCE) ecosystem, USA, and ultrafiltered dissolved organic matter (UDOM; 1 kDa) was isolated for characterization of its source, bioavailability and diagenetic state. A combination of techniques, including 15N cross-polarization magic angle spinning nuclear magnetic resonance (15N CPMAS NMR) and X-ray photoelectron spectroscopy (XPS), were used to analyze the N components of UDOM. The concentrations and compositions of total hydrolysable amino acids (HAAs) were analyzed to estimate UDOM bioavailability and diagenetic state. Optical properties (UV–visible and fluorescence) and the stable isotope ratios of C and N were measured to assess the source and dynamics of UDOM. Spectroscopic analyses consistently showed that the major N species of UDOM are in amide form, but significant contributions of aromatic-N were also observed. XPS showed a very high pyridinic-N concentration in the FCE–UDOM (21.7 ± 2.7%) compared with those in other environments. The sources of this aromatic-N are unclear, but could include soot and charred materials from wild fires. Relatively high total HAA concentrations (4 ± 2% UDOC or 27 ± 4% UDON) are indicative of bioavailable components, and HAA compositions suggest FCE–UDOM has not undergone extensive diagenetic processing. These observations can be attributed to the low microbial activity and a continuous supply of fresh UDOM in this oligotrophic ecosystem. Marsh plants appear to be the dominant source of UDOM in freshwater regions of the FCE, whereas seagrasses and algae are the dominant sources of UDOM in Florida Bay. This study demonstrates the utility of a multi-technique and multi-proxy approach to advance our understanding of DON biogeochemistry.

Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

The successional dynamics of lightning-initiated canopy gaps in the mangrove forests of Shark River, Everglades National Park, USA

Gap succession is a significant determinant of structure and development in most forest communities. Lightning strikes are an important source of canopy gaps in the mangrove forest of Everglades National Park. I investigated the successional dynamics of lightning-initiated canopy gaps to determine their influence on forest stand structure of the mixed mangrove forests ( Rhizophora mangle, Laguncularia racemosa, and Avicennia germinans ) of the Shark River. I measured gap size, gap shape, light environment, soil characteristics, woody debris, and fiddler crab abundance. I additionally measured the vegetative composition in a chronosequences of gap successional stages (new, recruiting, and growing gaps). I recorded survivorship, recruitment, growth and soil elevation dynamics within a subset of new and growing gaps. I determined the relationship between intact forest soil elevation and site hydrology in order to interpret the effects of lightning disturbance on soil elevation dynamics. ^ Gap size averaged 289 ± 20 m2 (± 1SE) and light transmittance decreased exponentially as gaps filled with saplings. Fine woody debris was highest in recruiting gaps. Soil strength was lower in the gaps than in the forest. The abundance of large and medium fiddler crab burrows increased linearly with total seedling abundance. Soil surface elevation declined in newly formed lightning gaps; this loss was due to a combination of superficial erosion (8.5 mm) and subsidence (60.9 mm). A distinct two-cohort recruitment pattern was evident in the seedling/sapling surveys, suggesting a partitioning of the succession between individuals present before and after lightning strike. In new gaps, the seedling recruitment rate was twice as high as in forest and the sapling population increased. At the growing gap stage, R. mangle seedling mortality was 10 times greater and sapling mortality was 13 times greater than recruitment. Growing gaps had reduced seedling stem elongation, sapling growth and adult growth. However, a few individuals (R. mangle saplings) were able to recruit into the adult life stage. In conclusion, the high density of R. mangle seedlings and saplings imply that lightning strike disturbances in these mangrove forests favor their recruitment over that of A. germinans and L. racemosa. ^

Quantifying nitrogen cycling rates in freshwater marshes of the southern Everglades using nitrogen-15 tracer techniques

This dissertation research project addressed the question of how hydrologic restoration of the Everglades is impacting the nutrient dynamics of marsh ecosystems in the southern Everglades. These effects were analyzed by quantifying nitrogen (N) cycle dynamics in the region. I utilized stable isotope tracer techniques to investigate nitrogen uptake and cycling between the major ecosystem components of the freshwater marsh system. I recorded the natural isotopic signatures (δ15N and δ 13C) for major ecosystem components from the three major watersheds of the Everglades: Shark River Slough, Taylor Slough, and C-111 basin. Analysis of δ15 N and δ13C natural abundance data were used to demonstrate the spatial extent to which nitrogen from anthropogenic or naturally enriched sources is entering the marshes of the Everglades. In addition, I measured the fluxes on N between various ecosystem components at both near-canal and estuarine ecotone locations. Lastly, I investigated the effect of three phosphorus load treatments (0.00 mg P m-2, 6.66 mg P m-2, and 66.6 mg P m-2) on the rate and magnitude of ecosystem N-uptake and N-cycling. The δ15N and δ13C natural abundance data supported the hypothesis that ecosystem components from near-canal sites have heavier, more enriched δ 15N isotopic signatures than downstream sites. The natural abundance data also showed that the marshes of the southern Everglades are acting as a sink for isotopically heavier, canal-borne dissolved inorganic nitrogen (DIN) and a source for "new" marsh derived dissolved organic nitrogen (DON). In addition, the 15N mesocosm data showed the rapid assimilation of the 15N tracer by the periphyton component and the delayed N uptake by soil and macrophyte components in the southern Everglades.

Hydrological dynamics between a coastal aquifer and the adjacent estuarine system, Biscayne Bay, South Florida

Geochemical and geophysical approaches have been used to investigate the freshwater and saltwater dynamics in the coastal Biscayne Aquifer and Biscayne Bay. Stable isotopes of oxygen and hydrogen, and concentrations of Sr2+ and Ca2+ were combined in two geochemical mixing models to provide estimates of the various freshwater inputs (precipitation, canal water, and groundwater) to Biscayne Bay and the coastal canal system in South Florida. Shallow geophysical electromagnetic and direct current resistivity surveys were used to image the geometry and stratification of the saltwater mixing zone in the near coastal (less than 1km inland) Biscayne Aquifer. The combined stable isotope and trace metal models suggest a ratio of canal input-precipitation-groundwater of 38%–52%–10% in the wet season and 37%–58%–5% in the dry season with an error of 25%, where most (20%) of the error was attributed to the isotope regression model, while the remaining 5% error was attributed to the Sr2+/Ca2+ mixing model. These models suggest rainfall is the dominate source of freshwater to Biscayne Bay. For a bay-wide water budget that includes saltwater and freshwater mixing, fresh groundwater accounts for less than 2% of the total input. A similar Sr 2+/Ca2+ tracer model indicates precipitation is the dominate source in 9 out of 10 canals that discharge into Biscayne Bay. The two-component mixing model converged for 100% of the freshwater canal samples in this study with 63% of the water contributed to the canals coming from precipitation and 37% from groundwater inputs ±4%. There was a seasonal shift from 63% precipitation input in the dry season to 55% precipitation input in the wet season. The three end-member mixing model converged for only 60% of the saline canal samples possibly due to non-conservative behavior of Sr2+ and Ca2+ in saline groundwater discharging into the canal system. Electromagnetic and Direct Current resistivity surveys were successful at locating and estimating the geometry and depth of the freshwater/saltwater interface in the Biscayne Aquifer at two near coastal sites. A saltwater interface that deepened as the survey moved inland was detected with a maximum interpreted depth to the interface of 15 meters, approximately 0.33 km inland from the shoreline. ^

Estimation of ravine sediment production and sediment dynamics in the lower Le Sueur River watershed, Minnesota

This study focuses on quantifying explicitly the sediment budget of deeply incised ravines in the lower Le Sueur River watershed, in southern Minnesota. High-rate-gully-erosion equations along with the Universal Soil Loss Equation (USLE) were implemented in a numerical modeling approach that is based on a time-integration of the sediment balance equations. The model estimates the rates of ravine width and depth change and the amount of sediment periodically flushing from the ravines. Components of the sediment budget of the ravines were simulated with the model and results suggest that the ravine walls are the major sediment source in the ravines. A sensitivity analysis revealed that the erodibility coefficients of the gully bed and wall, the local slope angle and the Manning’s coefficient are the key parameters controlling the rate of sediment production. Recommendations to guide further monitoring efforts in the watershed and increased detail modeling approaches are highlighted as a result of this modeling effort.

Multi-Scaled Grassland-Woody Plant Dynamics in the Heterogeneous Marl Prairies of the Southern Everglades

The Everglades freshwater marl prairie is a dynamic and spatially heterogeneous landscape, containing thousands of tree islands nested within a marsh matrix. Spatial processes underlie population and community dynamics across the mosaic, especially the balance between woody and graminoid components, and landscape patterns reflect interactions among multiple biotic and abiotic drivers. To better understand these complex, multi-scaled relationships we employed a three-tiered hierarchical design to investigate the effects of seed source, hydrology, and more indirectly fire on the establishment of new woody recruits in the marsh, and to assess current tree island patterning across the landscape. Our analyses were conducted at the ground level at two scales, which we term the micro- and meso-scapes, and results were related to remotely detected tree island distributions assessed in the broader landscape, that is, the macro-scape. Seed source and hydrologic effects on recruitment in the micro- and meso-scapes were analyzed via logistic regression, and spatial aggregation in the macro-scape was evaluated using a grid-based univariate O-ring function. Results varied among regions and scales but several general trends were observed. The patterning of adult populations was the strongest driver of recruitment in the micro- and meso-scape prairies, with recruits frequently aggregating around adults or tree islands. However in the macro-scape biologically associated (second order) aggregation was rare, suggesting that emergent woody patches are heavily controlled by underlying physical and environmental factors such as topography, hydrology, and fire.

Effects of hydrologic and water quality drivers on periphyton dynamics in the southern Everglades

Everglades periphyton mats are tightly-coupled autotrophic (algae and cyanobacteria) and heterotrophic (eubacteria, fungi and microinvertebrates) microbial assemblages. We investigated the effect of water column total phosphorus and nitrogen concentrations, water depth and hydroperiod on periphyton of net production, respiration, nutrient content, and biomass. Our study sites were located along four transects that extended southward with freshwater sheetflow through sawgrass-dominated marsh. The water source for two of the transects were canal-driven and anchored at canal inputs. The two other transects were rain-driven (ombrotrophic) and began in sawgrass-dominated marsh. Periphyton dynamics were examined for upstream and downstream effects within and across the four transects. Although all study sites were characterized as short hydroperiod and phosphorus-limited oligotrophic, they represent gradients of hydrologic regime, water source and water quality of the southern Everglades. Average periphyton net production of 1.08 mg C AFDW−1 h−1 and periphyton whole system respiration of 0.38 mg C AFDW−1 h−1 rates were net autotrophic. Biomass was generally highest at ombrotrophic sites and sites downstream of canal inputs. Mean biomass over all our study sites was high, 1517.30 g AFDW m−2. Periphyton was phosphorus-limited. Average periphyton total phosphorus content was 137.15 μg P g−1 and average periphyton total N:P ratio was 192:1. Periphyton N:P was a sensitive indicator of water source. Even at extremely low mean water total phosphorus concentrations ( ≤ 0.21 μmol l−1), we found canal source effects on periphyton dynamics at sites adjacent to canal inputs, but not downstream of inflows. These canal source effects were most pronounced at the onset of wet season with initial rewetting. Spatial and temporal variability in periphyton dynamics could not solely be ascribed to water quality, but was often associated with both hydrology and water source.

Dynamics of ecosystem metabolism and flocculent detritus transport in estuarine Taylor River

Estuaries and estuarine wetlands are ecologically and societally important systems, exhibiting high rates of primary production that fuel offshore secondary production. Hydrological processes play a central role in shaping estuarine ecosystem structure and function by controlling nutrient loading and the relative contributions of marine and terrestrial influences on the estuary. The Comprehensive Everglades Restoration Plan includes plans to restore freshwater delivery to Taylor Slough, a shallow drainage basin in the southern Everglades, ultimately resulting in increased freshwater flow to the downstream Taylor River estuary. The existing seasonal and inter-annual variability of water flow and source in Taylor River affords the opportunity to investigate relationships between ecosystem function and hydrologic forcing. Estimates of aquatic ecosystem metabolism, derived from free-water, diel changes in dissolved oxygen, were combined with assessments of wetland flocculent detritus quality and transport within the context of seasonal changes in Everglades hydrology. Variation in ecosystem gross primary production and respiration were linked to seasonal changes in estuarine water quality using multiple autoregression models. Furthermore, Taylor River was observed to be net heterotrophic, indicating that an allochthonous source of carbon maintained ecosystem respiration in excess of autochthonous primary production. Wetland-derived detritus appears to be an important vector of energy and nutrients across the Everglades landscape; and in Taylor River, is seasonally flushed into ponded segments of the river where it is then respired. Lastly, seasonal water delivery appears to govern feedbacks regulating water column phosphorus availability in the Taylor River estuary.

Seagrass sediments as a global carbon sink: Isotopic constraints

Seagrass meadows are highly productive habitats found along many of the world's coastline, providing important services that support the overall functioning of the coastal zone. The organic carbon that accumulates in seagrass meadows is derived not only from seagrass production but from the trapping of other particles, as the seagrass canopies facilitate sedimentation and reduce resuspension. Here we provide a comprehensive synthesis of the available data to obtain a better understanding of the relative contribution of seagrass and other possible sources of organic matter that accumulate in the sediments of seagrass meadows. The data set includes 219 paired analyses of the carbon isotopic composition of seagrass leaves and sediments from 207 seagrass sites at 88 locations worldwide. Using a three source mixing model and literature values for putative sources, we calculate that the average proportional contribution of seagrass to the surface sediment organic carbon pool is ∼50%. When using the best available estimates of carbon burial rates in seagrass meadows, our data indicate that between 41 and 66 gC m−2 yr−1 originates from seagrass production. Using our global average for allochthonous carbon trapped in seagrass sediments together with a recent estimate of global average net community production, we estimate that carbon burial in seagrass meadows is between 48 and 112 Tg yr−1, showing that seagrass meadows are natural hot spots for carbon sequestration.

Changes in mass and nutrient content of wood during decomposition in a south Florida mangrove forest

1. Large pools of dead wood in mangrove forests following disturbances such as hurricanes may influence nutrient fluxes. We hypothesized that decomposition of wood of mangroves from Florida, USA (Avicennia germinans, Laguncularia racemosa and Rhizophora mangle), and the consequent nutrient dynamics, would depend on species, location in the forest relative to freshwater and marine influences and whether the wood was standing, lying on the sediment surface or buried. 2. Wood disks (8–10 cm diameter, 1 cm thick) from each species were set to decompose at sites along the Shark River, either buried in the sediment, on the soil surface or in the air (above both the soil surface and high tide elevation). 3. A simple exponential model described the decay of wood in the air, and neither species nor site had any effect on the decay coefficient during the first 13 months of decomposition. 4. Over 28 months of decomposition, buried and surface disks decomposed following a two-component model, with labile and refractory components. Avicennia germinans had the largest labile component (18 ± 2% of dry weight), while Laguncularia racemosa had the lowest (10 ± 2%). Labile components decayed at rates of 0.37–23.71% month−1, while refractory components decayed at rates of 0.001–0.033% month−1. Disks decomposing on the soil surface had higher decay rates than buried disks, but both were higher than disks in the air. All species had similar decay rates of the labile and refractory components, but A. germinans exhibited faster overall decay because of a higher proportion of labile components. 5. Nitrogen content generally increased in buried and surface disks, but there was little change in N content of disks in the air over the 2-year study. Between 17% and 68% of total phosphorus in wood leached out during the first 2 months of decomposition, with buried disks having the greater losses, P remaining constant or increasing slightly thereafter. 6. Newly deposited wood from living trees was a short-term source of N for the ecosystem but, by the end of 2 years, had become a net sink. Wood, however, remained a source of P for the ecosystem. 7. As in other forested ecosystems, coarse woody debris can have a significant impact on carbon and nutrient dynamics in mangrove forests. The prevalence of disturbances, such as hurricanes, that can deposit large amounts of wood on the forest floor accentuates the importance of downed wood in these forests.

Summertime Influences of Tidal Energy Advection on the Surface Energy Balance in a Mangrove Forest

Mangrove forests are ecosystems susceptible to changing water levels and temperatures due to climate change as well as perturbations resulting from tropical storms. Numerical models can be used to project mangrove forest responses to regional and global environmental changes, and the reliability of these models depends on surface energy balance closure. However, for tidal ecosystems, the surface energy balance is complex because the energy transport associated with tidal activity remains poorly understood. This study aimed to quantify impacts of tidal flows on energy dynamics within a mangrove ecosystem. To address the research objective, an intensive 10-day study was conducted in a mangrove forest located along the Shark River in the Everglades National Park, FL, USA. Forest–atmosphere turbulent exchanges of energy were quantified with an eddy covariance system installed on a 30-m-tall flux tower. Energy transport associated with tidal activity was calculated based on a coupled mass and energy balance approach. The mass balance included tidal flows and accumulation of water on the forest floor. The energy balance included temporal changes in enthalpy, resulting from tidal flows and temperature changes in the water column. By serving as a net sink or a source of available energy, flood waters reduced the impact of high radiational loads on the mangrove forest. Also, the regression slope of available energy versus sink terms increased from 0.730 to 0.754 and from 0.798 to 0.857, including total enthalpy change in the water column in the surface energy balance for 30-min periods and daily daytime sums, respectively. Results indicated that tidal inundation provides an important mechanism for heat removal and that tidal exchange should be considered in surface energy budgets of coastal ecosystems. Results also demonstrated the importance of including tidal energy advection in mangrove biophysical models that are used for predicting ecosystem response to changing climate and regional freshwater management practices.