Flow and water budgets in the C-111 and L-31W canals near the Everglades National Park, Dade County, Florida

Acoustic velocity meter (AVM) sites, located both distant and adjacent to canal water control structures, were constructed and calibrated in L-31W borrow canal and Canal 111 (C-111) to measure canal water velocity. Data were used to compute monthly discharge volumes and overall water budgets for several canal reaches from August 1994 to May 1996. The water budgets indicated extensive aquifer inflows in L-31W associated, in part, with S-332 pump station return flows. Canal and groundwater piezometer data showed 5 distinct hydrologic scenarios (distinguished by the direction and magnitude of hydraulic gradients) in the important Frog Pond area on the eastern boundary of the Everglades National Park. Most of the water lost from C-111 was via surface water losses near the outlet of the system, close to Florida Bay. The distribution of flows during the study suggest an alteration of the present South Dade Conveyance System modification plan to improve water deliveries to Taylor Slough and the Eastern Panhandle of the Everglades National Park. ^

Lattice Boltzmann method for flow and heat transfer in microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

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 variability in Florida Bay particulate organic matter composition: combining flow cytometry with stable isotope analyses

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 Florida Bay biogeochemistry. Planktonic communities respond quickly to changes in water quality, thus spatial variability in community composition and relationships to nutrient parameters must be understood in order to evaluate future downstream impacts of modifications to Everglades hydrology. Here we present initial results combining flow cytometry analyses of phytoplankton and bacterial populations (0.1–50 μm size fraction) with measurements of δ13C and δ15N composition and dissolved inorganic nutrient concentrations to explore proxies for planktonic species assemblage compositions and nutrient cycling. Particulate organic material in the 0.1–50 μm size fraction was collected from five stations in Northeastern and Western Florida Bay to characterize spatial variability in species assemblage and stable isotopic composition. A dense bloom of the picocyanobacterium, Synechococcus elongatus, was observed at Western Florida Bay sites. Smaller Synechococcus sp. were present at Northeast sites in much lower abundance. Bacteria and detrital particles were also more abundant at Western Florida Bay stations than in the northeast region. The highest abundance of detritus occurred at Trout Creek, which receives freshwater discharge from the Everglades through Taylor Slough. In terms of nutrient availability and stable isotopic values, the S. elongatus population in the Western bay corresponded to low DIN (0.5 μM NH 4 + ; 0.2 μM NO 3 − ) concentrations and depleted δ15N signatures ranging from +0.3 to +0.8‰, suggesting that the bloom supported high productivity levels through N2-fixation. δ15N values from the Northeast bay were more enriched (+2.0 to +3.0‰), characteristic of N-recycling. δ13C values were similar for all marine Florida Bay stations, ranging from −17.6 to −14.4‰, however were more depleted at the mangrove ecotone station (−25.5 to −22.3‰). The difference in the isotopic values reflects differences in carbon sources. These findings imply that variations in resource availability and nutrient sources exert significant control over planktonic community composition, which is reflected by stable isotopic signatures.

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.

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.

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.

Relating Freshwater Flow with Estuarine Water Quality in the Southern Everglades Mangrove Ecotone

Florida Bay is more saline than it was historically, and reduced freshwater flows may lead to more phosphorus inputs to the mangrove ecotone from the marine end-member. This is important given plans to restore freshwater flow into eastern Florida Bay. We investigated the relationships between salinity, nutrients, and hydrologic variables in the mangrove ecotone of Taylor Slough. We expected that total phosphorus (TP) would increase with salinity, reflecting a downstream marine source, while total nitrogen (TN) would increase with flow in the mangrove ecotone. Despite expectations of increased flows improving the ecological health of lower Taylor Slough and Florida Bay, total nitrogen (TN) and total phosphorus (TP) dynamics may shift in response to new conditions of flow and salinity as well as organic carbon, N, and P availability. Our results showed that TP concentrations are more discharge-driven while TN is more variable and potentially derived from different sources along the flow path from the freshwater Everglades marshes to Florida Bay. Increased flow of freshwater through Taylor Slough will likely decrease TP concentrations in this historically oligotrophic and P-limited ecosystem. However, more studies along the mangrove ecotone is needed to understand how increased flows will affect nitrogen dynamics relative to phosphorus.

Taylor Slough Hydrology

Taylor Slough, in Everglades National Park, has experienced an evolution of water management infrastructure since drainage activities arrived in South Florida. This has included the excavation of canals, installation of large capacity pump stations, and a variety of operational strategies focused on resolving the conflict between managing the water level for developed areas while providing water supply for Everglades National Park. This study provides a review of water management practices and the concurrent hydrologic conditions in the Taylor Slough basin and adjacent canal system from 1961 through 2010. Analyses of flow, water level and rainfall data were divided into time periods that correspond to significant changes in structural features and operational plans. In the early 1960s, Taylor Slough was disconnected from the greater Everglades system by the construction of levees upstream. As water supply for Taylor Slough became more urgent, the Slough was connected to the regional water supply system via a network of canals and pump stations to relieve over-drained conditions. The increased water supply and pump capacity succeeded in raising water level and increasing flow and hydroperiod in the marsh.

Benthic Exchange of C, N, and P Along the Estuarine Ecotone of Lower Taylor Slough, Florida (USA): Effect of Seasonal Flows and Phosphorus Availability

The southern Everglades and Florida Bay have experienced a nearly 50 % reduction in freshwater flow resulting in increased salinity and landward expansion of mangrove forest. Given the marine end-member is a natural source of P to this region, it is necessary to understand the interactions between inflows and P availability in controlling the exchange of materials across the mangrove ecotone. From 2007 to 2008, we used sediment core incubations to quantify fluxes of dissolved inorganic N and P and dissolved organic carbon (DOC) in three ecotone areas (dwarf mangrove, pond, and bay). Experiments were repeated seasonally over 2 years involving P-enriched surface water as a factor. We saw consistent uptake of soluble reactive P (SRP), DOC, and nitrate + nitrite (N+N) by the soils/sediments and release of ammonium (NH4 +) from soils/sediments to the water column across all sites and seasons. P enrichment had no discernible effect on DIN or DOC flux, suggesting that rapid P uptake may have been more geochemically mediated. However, uptake of added P occurred across all sites and seasons, reflecting high uptake capacity in this carbonate system and the potential of the mangrove ecotone to sequester P as it becomes more available.

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.

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).