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.

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.

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.

Ecohydrology, evapotranspiration and hydrogeochemistry of carbonate mangrove wetlands

Coastal environments can be highly susceptible to environmental changes caused by anthropogenic pressures and natural events. Both anthropogenic and natural perturbations may directly affect the amount and the quality of water flowing through the ecosystem, both in the surface and subsurface and can subsequently, alter ecological communities and functions. The Florida Everglades and the Sian Ka'an Biosphere Reserve (Mexico) are two large ecosystems with an extensive coastal mangrove ecotone that represent a historically altered and pristine environment, respectively. Rising sea levels, climate change, increased water demand, and salt water intrusion are growing concerns in these regions and underlies the need for a better understanding of the present conditions. The goal of my research was to better understand various ecohydrological, environmental, and hydrogeochemical interactions and relationships in carbonate mangrove wetlands. A combination of aqueous geochemical analyses and visible and near-infrared reflectance data were employed to explore relationships between surface and subsurface water chemistry and spectral biophysical stress in mangroves. Optical satellite imagery and field collected meteorological data were used to estimate surface energy and evapotranspiration and measure variability associated with hurricanes and restoration efforts. Furthermore, major ionic and nutrient concentrations, and stable isotopes of hydrogen and oxygen were used to distinguish water sources and infer coastal groundwater discharge by applying the data to a combined principal component analysis-end member mixing model. Spectral reflectance measured at the field and satellite scales were successfully used to estimate surface and subsurface water chemistry and model chloride concentrations along the southern Everglades. Satellite imagery indicated that mangrove sites that have less tidal flushing and hydrogeomorphic heterogeneity tend to have more variable evapotranspiration and soil heat flux in response to storms and restoration. Lastly, water chemistry and multivariate analyses indicated two distinct fresh groundwater sources that discharge to the phosphorus-limited estuaries and bays of the Sian Ka'an Biopshere Reserve; and that coastal groundwater discharge was an important source for phosphorus. The results of the study give us a better understanding of the ecohydrological and hydrogeological processes in carbonate mangrove environments that can be then be extrapolated to similar coastal ecosystems in the Caribbean.

Quantifying the effects of temperature and concentration on variable-density flow in numerical modeling of groundwater systems: Implications for predictive uncertainty and data collection

Groundwater systems of different densities are often mathematically modeled to understand and predict environmental behavior such as seawater intrusion or submarine groundwater discharge. Additional data collection may be justified if it will cost-effectively aid in reducing the uncertainty of a model's prediction. The collection of salinity, as well as, temperature data could aid in reducing predictive uncertainty in a variable-density model. However, before numerical models can be created, rigorous testing of the modeling code needs to be completed. This research documents the benchmark testing of a new modeling code, SEAWAT Version 4. The benchmark problems include various combinations of density-dependent flow resulting from variations in concentration and temperature. The verified code, SEAWAT, was then applied to two different hydrological analyses to explore the capacity of a variable-density model to guide data collection. ^ The first analysis tested a linear method to guide data collection by quantifying the contribution of different data types and locations toward reducing predictive uncertainty in a nonlinear variable-density flow and transport model. The relative contributions of temperature and concentration measurements, at different locations within a simulated carbonate platform, for predicting movement of the saltwater interface were assessed. Results from the method showed that concentration data had greater worth than temperature data in reducing predictive uncertainty in this case. Results also indicated that a linear method could be used to quantify data worth in a nonlinear model. ^ The second hydrological analysis utilized a model to identify the transient response of the salinity, temperature, age, and amount of submarine groundwater discharge to changes in tidal ocean stage, seasonal temperature variations, and different types of geology. The model was compared to multiple kinds of data to (1) calibrate and verify the model, and (2) explore the potential for the model to be used to guide the collection of data using techniques such as electromagnetic resistivity, thermal imagery, and seepage meters. Results indicated that the model can be used to give insight to submarine groundwater discharge and be used to guide data collection. ^

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.

Physical factors influencing the distribution of a top predator in a subtropical oligotrophic estuary

We used longline fishing to determine the effects of distance from the ocean, season, and short-term variation in abiotic conditions on the abundance of juvenile bull sharks (Carcharhinus leucas) in an estuary of the Florida Everglades, U.S.A. Logistic regression revealed that young-of-the-year sharks were concentrated at a protected site 20 km upstream and were present in greater abundance when dissolved oxygen (DO) levels were high. For older juvenile sharks (age 1+), DO levels had the greatest influence on catch probabilities followed by distance from the ocean; they were most likely to be caught at sites with .3.5 mg L21 DO and on the main branch of the river 20 km upstream. Salinity had a relatively small effect on catch rates and there were no seasonal shifts in shark distribution. Our results highlight the importance of considering DO as a possible driver of top predator distributions in estuaries, even in the absence of hypoxia. In Everglades estuaries hydrological drivers that affect DO levels (e.g., groundwater discharge, modification of primary productivity through nutrient fluxes) will be important in determining shark distributions, and the effects of planned ecosystem restoration efforts on bull sharks will not simply be mediated by changing salinity regimes and the location of the oligohaline zone. More generally, variation in DO levels could structure the nature and spatiotemporal pattern of top predator effects in the coastal Everglades, and other tropical and subtropical estuaries, because of interspecific variation in reliance on DO within the top predator guild.

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. ^

Comparative study of dissolved organic matter from groundwater and surface water in the Florida coastal Everglades using multi-dimensional spectrofluorometry combined with multivariate statistics

Dissolved organic matter (DOM) in groundwater and surface water samples from the Florida coastal Everglades were studied using excitation–emission matrix fluorescence modeled through parallel factor analysis (EEM-PARAFAC). DOM in both surface and groundwater from the eastern Everglades S332 basin reflected a terrestrial-derived fingerprint through dominantly higher abundances of humic-like PARAFAC components. In contrast, surface water DOM from northeastern Florida Bay featured a microbial-derived DOM signature based on the higher abundance of microbial humic-like and protein-like components consistent with its marine source. Surprisingly, groundwater DOM from northeastern Florida Bay reflected a terrestrial-derived source except for samples from central Florida Bay well, which mirrored a combination of terrestrial and marine end-member origin. Furthermore, surface water and groundwater displayed effects of different degradation pathways such as photodegradation and biodegradation as exemplified by two PARAFAC components seemingly indicative of such degradation processes. Finally, Principal Component Analysis of the EEM-PARAFAC data was able to distinguish and classify most of the samples according to DOM origins and degradation processes experienced, except for a small overlap of S332 surface water and groundwater, implying rather active surface-to-ground water interaction in some sites particularly during the rainy season. This study highlights that EEM-PARAFAC could be used successfully to trace and differentiate DOM from diverse sources across both horizontal and vertical flow profiles, and as such could be a convenient and useful tool for the better understanding of hydrological interactions and carbon biogeochemical cycling.

Carbon and nitrogen stable isotopic patterns in South Florida coastal ecosystems: Modern and paleoceanographic perspectives

Long term management plans for restoration of natural flow conditions through the Everglades increase the importance of understanding potential nutrient impacts of increased freshwater delivery on coastal biogeochemistry. The present study sought to increase understanding of the coastal marine system of South Florida under modern conditions and through the anthropogenic changes in the last century, on scales ranging from individual nutrient cycle processes to seasonal patterns in organic material (OM) under varying hydrodynamic regime, to century scale analysis of sedimentary records. In all applications, carbon and nitrogen stable isotopic compositions of OM were examined as natural recorders of change and nutrient cycling in the coastal system. ^ High spatial and temporal variability in stable isotopic compositions were observed on all time scales. During a transient phytoplankton bloom, δ 15N values suggested nitrogen fixation as a nutrient source supporting enhanced productivity. Seasonally, particulate organic material (POM) from ten sites along the Florida Reef Tract and in Florida Bay demonstrated variable fluctuations dependent on hydrodynamic setting. Three separate intra-annual patterns were observed, yet statistical differences were observed between groupings of Florida Bay and Atlantic Ocean sites. The POM δ 15N values ranged on a quarterly basis by 7‰, while δ 13C varied by 22‰. From a sediment history perspective, four cores collected from Florida Bay further demonstrated the spatial and temporal variability of the system in isotopic composition of bulk OM over time. Source inputs of OM varied with location, with terrestrial inputs dominating proximal to Everglades freshwater discharge, seagrasses dominating in open estuary cores, and a marine mixture of phytoplankton and seagrass in a core from the boundary zone between Florida Bay and the Gulf of Mexico. Significant shifts in OM geochemistry were observed coincident with anthropogenic events of the 20th century, including railroad and road construction in the Florida Keys and Everglades, and also the extensive drainage changes in Everglades hydrology. The sediment record also preserved evidence of the major hurricanes of the last century, with excursions in geochemical composition coincident with Category 4-5 storms. ^

Characterization, sources, and transformations of dissolved organic matter (DOM) in Florida coastal Everglades (FCE)

Dissolved organic matter (DOM) is one of the largest carbon reservoirs on this planet and is present in aquatic environments as a highly complex mixture of organic compounds. The Florida coastal Everglades (FCE) is one of the largest wetlands in the world. DOM in this system is an important biogeochemical component as most of the nitrogen (N) and phosphorous (P) are in organic forms. Achieving a better understanding of DOM dynamics in large coastal wetlands is critical, and a particularly important issue in the context of Everglades restoration. In this work, the environmental dynamics of surface water DOM on spatial and temporal scales was investigated. In addition, photo- and bio-reactivity of this DOM was determined, surface-to-groundwater exchange of DOM was investigated, and the size distribution of freshwater DOM in Everglades was assessed. The data show that DOM dynamics in this ecosystem are controlled by both hydrological and ecological drivers and are clearly different on spatial scales and variable seasonally. The DOM reactivity data, modeled with a multi-pool first order degradation kinetics model, found that fluorescent DOM in FCE is generally photo-reactive and bio-refractory. Yet the sequential degradation proved a “priming effect” of sunlight on the bacterial uptake and reworking of this subtropical wetland DOM. Interestingly, specific PARAFAC components were found to have different photo- and bio-degradation rates, suggesting a highly heterogeneous nature of fluorophores associated with the DOM. Surface-to-groundwater exchange of DOM was observed in different regions of the system, and compositional differences were associated with source and photo-reactivity. Lastly, the high degree of heterogeneity of DOM associated fluorophores suggested based on the degradation studies was confirmed through the EEM-PARAFAC analysis of DOM along a molecular size continuum, suggesting that the fluorescence characteristics of DOM are highly controlled by different size fractions and as such can exhibit significant differences in reactivity.

Xavier Zapata (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.

A synthesis of long-term research by the Florida Coastal Everglades LTER Program

This paper synthesizes research conducted during the first 5–6 years of the Florida Coastal Everglades Long-Term Ecological Research Program (FCE LTER). My objectives are to review our research to date, and to present a new central theme and conceptual approach for future research. Our research has focused on understanding how dissolved organic matter (DOM) from upstream oligotrophic marshes interacted with a marine source of the limiting nutrient, phosphorus (P), to control productivity in the oligohaline estuarine ecotone. We have been working along freshwater to marine transects in two drainage basins located in Everglades National Park (ENP). The Shark River Slough transect (SRS) has a direct connection to the Gulf of Mexico, providing this estuarine ecotone with a source of marine P. The oligohaline ecotone along our southern Everglades transect (TS/Ph), however, is separated from this marine P source by the Florida Bay estuary. We originally hypothesized an ecosystem productivity peak in the SRS ecotone, driven by the interaction of marine P and Everglades DOM, but no such productivity peak in the TS/Ph ecotone because of this lack of marine P. Our research to date has tended to show the opposite pattern, however, with many ecosystem components showing enhanced productivity in the TS/Ph ecotone, but not in the SRS ecotone. Water column P concentrations followed a similar pattern, with unexpectedly high P in the TS/Ph ecotone during the dry season. Our organic geochemical research has shown that Everglades DOM is more refractory than originally hypothesized. We have also begun to understand the importance of detrital organic matter production and transport to ecotone dynamics and as the base of aquatic food webs. Our future research will build on this substantial body of knowledge about these oligotrophic estuaries. We will direct our efforts more strongly on biophysical dynamics in the oligohaline ecotone regions. Specifically, we will be focusing on inputs to these regions from four primary water sources: freshwater Everglades runoff, net precipitation, marine inputs, and groundwater. We are hypothesizing that dry season groundwater inputs of P will be particularly important to TS/Ph ecotone dynamics because of longer water residence times in this area. Our organic geochemical, biogeochemical, and ecosystem energetics work will focus more strongly on the importance of detrital organics and will take advantage of a key Everglades Restoration project, scheduled for 2008 or 2009, that will increase freshwater inputs to our SRS transect only. Finally, we will also begin to investigate the human dimensions of restoration, and of a growing population in south Florida that will become increasingly dependent on the Everglades for critical ecosystem services (including fresh water) even as its growth presents challenges to Everglades sustainability.

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.