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.

Development of a mathematical and computer model to assess water quality impacts of hazardous compounds from minor gasoline spills in the inter-tidal zone of the Miami River, Florida

The objective of this study was to develop a model to predict transport and fate of gasoline components of environmental concern in the Miami River by mathematically simulating the movement of dissolved benzene, toluene, xylene (BTX), and methyl-tertiary-butyl ether (MTBE) occurring from minor gasoline spills in the inter-tidal zone of the river. Computer codes were based on mathematical algorithms that acknowledge the role of advective and dispersive physical phenomena along the river and prevailing phase transformations of BTX and MTBE. Phase transformations included volatilization and settling. ^ The model used a finite-difference scheme of steady-state conditions, with a set of numerical equations that was solved by two numerical methods: Gauss-Seidel and Jacobi iterations. A numerical validation process was conducted by comparing the results from both methods with analytical and numerical reference solutions. Since similar trends were achieved after the numerical validation process, it was concluded that the computer codes algorithmically were correct. The Gauss-Seidel iteration yielded at a faster convergence rate than the Jacobi iteration. Hence, the mathematical code was selected to further develop the computer program and software. The model was then analyzed for its sensitivity. It was found that the model was very sensitive to wind speed but not to sediment settling velocity. ^ A computer software was developed with the model code embedded. The software was provided with two major user-friendly visualized forms, one to interface with the database files and the other to execute and present the graphical and tabulated results. For all predicted concentrations of BTX and MTBE, the maximum concentrations were over an order of magnitude lower than current drinking water standards. It should be pointed out, however, that smaller concentrations than the latter reported standards and values, although not harmful to humans, may be very harmful to organisms of the trophic levels of the Miami River ecosystem and associated waters. This computer model can be used for the rapid assessment and management of the effects of minor gasoline spills on inter-tidal riverine water quality. ^

Design Optimization of Submerged Jet Nozzles for Enhanced Mixing

The purpose of this thesis was to identify the optimal design parameters for a jet nozzle which obtains a local maximum shear stress while maximizing the average shear stress on the floor of a fluid filled system. This research examined how geometric parameters of a jet nozzle, such as the nozzle's angle, height, and orifice, influence the shear stress created on the bottom surface of a tank. Simulations were run using a Computational Fluid Dynamics (CFD) software package to determine shear stress values for a parameterized geometric domain including the jet nozzle. A response surface was created based on the shear stress values obtained from 112 simulated designs. A multi-objective optimization software utilized the response surface to generate designs with the best combination of parameters to achieve maximum shear stress and maximum average shear stress. The optimal configuration of parameters achieved larger shear stress values over a commercially available design.

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.

Derivation of lake mixing and stratification indices from high-resolution lake buoy data

Lake Analyzer is a numerical code coupled with supporting visualization tools for determining indices of mixing and stratification that are critical to the biogeochemical cycles of lakes and reservoirs. Stability indices, including Lake Number, Wedderburn Number, Schmidt Stability, and thermocline depth are calculated according to established literature definitions and returned to the user in a time series format. The program was created for the analysis of high-frequency data collected from instrumented lake buoys, in support of the emerging field of aquatic sensor network science. Available outputs for the Lake Analyzer program are: water temperature (error-checked and/or down-sampled), wind speed (error-checked and/or down-sampled), metalimnion extent (top and bottom), thermocline depth, friction velocity, Lake Number, Wedderburn Number, Schmidt Stability, mode-1 vertical seiche period, and Brunt-Väisälä buoyancy frequency. Secondary outputs for several of these indices delineate the parent thermocline depth (seasonal thermocline) from the shallower secondary or diurnal thermocline. Lake Analyzer provides a program suite and best practices for the comparison of mixing and stratification indices in lakes across gradients of climate, hydro-physiography, and time, and enables a more detailed understanding of the resulting biogeochemical transformations at different spatial and temporal scales.

Distribution of Diatoms and Development of Diatom-Based Models for Inferring Salinity and Nutrient Concentrations in Florida Bay and Adjacent Coastal Wetlands of South Florida (USA)

The composition and distribution of diatom algae inhabiting estuaries and coasts of the subtropical Americas are poorly documented, especially relative to the central role diatoms play in coastal food webs and to their potential utility as sentinels of environmental change in these threatened ecosystems. Here, we document the distribution of diatoms among the diverse habitat types and long environmental gradients represented by the shallow topographic relief of the South Florida, USA, coastline. A total of 592 species were encountered from 38 freshwater, mangrove, and marine locations in the Everglades wetland and Florida Bay during two seasonal collections, with the highest diversity occurring at sites of high salinity and low water column organic carbon concentration (WTOC). Freshwater, mangrove, and estuarine assemblages were compositionally distinct, but seasonal differences were only detected in mangrove and estuarine sites where solute concentration differed greatly between wet and dry seasons. Epiphytic, planktonic, and sediment assemblages were compositionally similar, implying a high degree of mixing along the shallow, tidal, and storm-prone coast. The relationships between diatom taxa and salinity, water total phosphorus (WTP), water total nitrogen (WTN), and WTOC concentrations were determined and incorporated into weighted averaging partial least squares regression models. Salinity was the most influential variable, resulting in a highly predictive model (r apparent 2  = 0.97, r jackknife 2  = 0.95) that can be used in the future to infer changes in coastal freshwater delivery or sea-level rise in South Florida and compositionally similar environments. Models predicting WTN (r apparent 2  = 0.75, r jackknife 2  = 0.46), WTP (r apparent 2  = 0.75, r jackknife 2  = 0.49), and WTOC (r apparent 2  = 0.79, r jackknife 2  = 0.57) were also strong, suggesting that diatoms can provide reliable inferences of changes in solute delivery to the coastal ecosystem.

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.

Using Optical Properties to Quantify Fringe Mangrove Inputs to the Dissolved Organic Matter (DOM) Pool in a Subtropical Estuary

Dissolved organic carbon (DOC) concentration and dissolved organic matter (DOM) optical properties were analyzed along two estuarine river transects during the wet and dry seasons to better understand DOM dynamics and quantify mangrove inputs. A tidal study was performed to assess the impacts of tidal pumping on DOM transport. DOM in the estuaries showed non-conservative mixing indicative of mangrove-derived inputs. Similarly, fluorescence data suggest that some terrestrial humic-like components showed non-conservative behavior. An Everglades freshwater-derived fluorescent component, which is associated with soil inputs from the Northern Everglades, behaved conservatively. During the dry season, a protein-like component behaved conservatively until the mid-salinity range when non-conservative behavior due to degradation and/or loss was observed. The tidal study data suggests mangrove porewater inputs to the rivers following low tide. The differences in quantity of DOM exported by the Shark and Harney Rivers imply that geomorphology and tidal hydrology may be a dominant factor controlling the amount of DOM exported from the mangrove ecotone, where up to 21 % of the DOC is mangrove-derived. Additionally, nutrient concentrations and other temporal factors may control DOM export from the mangroves, particularly for the microbially derived fluorescent components, contributing to the seasonal differences. The wet and dry season fluxes of mangrove DOM from the Shark River is estimated as 0.27 × 109 mg C d−1 and 0.075 × 109 mg C d−1, respectively, and the Harney River is estimated as 1.9 × 109 mg C d−1 and 0.20 × 109 mg C d−1.

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.

Jay Sah crossing tidal creek to measure tree damage by Hurricane Wilma, Harney River

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.

Analyzing Tidal Fluctuations in the Big Pine Key Freshwater Lens with Time-Lapse Resistivity

The tidal influence on the Big Pine Key saltwater/freshwater interface was analyzed using time-lapse electrical resistivity imaging and shallow well measurements. The transition zone at the saltwater/freshwater interface was measured over part of a tidal cycle along three profiles. The resistivity was converted to salinity by deriving a formation factor for the Miami Oolite. A SEAWAT model was created to attempt to recreate the field measurements and test previously established hydrogeologic parameters. The results imply that the tide only affects the groundwater within 20 to 30 m of the coast. The effect is small and caused by flooding from the high tide. The low relief of the island means this effect is very sensitive to small changes in the magnitude. The SEAWAT model proved to be insufficient in modeling this effect. The study suggests that the extent of flooding is the largest influence on the salinity of the groundwater.