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

Identification of a hydrodynamic threshold in karst rocks from the Biscayne Aquifer, south Florida, USA

A hydrodynamic threshold between Darcian and non-Darcian flow conditions was found to occur in cubes of Key Largo Limestone from Florida, USA (one cube measuring 0.2 m on each side, the other 0.3 m) at an effective porosity of 33% and a hydraulic conductivity of 10 m/day. Below these values, flow was laminar and could be described as Darcian. Above these values, hydraulic conductivity increased greatly and flow was non-laminar. Reynolds numbers (Re) for these experiments ranged from

Geochemical determination of the fate and transport of injected fresh wastewater to a deep saline aquifer

Deep well injection into non-potable saline aquifers of treated domestic wastewater has been used in Florida for decades as a safe and effective alternative to ocean outfall disposal. The objectives of this study were to determine the fate and transport of injected wastewater at two deep well injection sites in Miami Dade County, Florida, USA. Detection of ammonium in the Middle Confining units of the Floridan aquifer above the injection zone at both sites has been interpreted as evidence of upward migration of injected wastewater, posing a risk to underground sources of drinking water. Historical water quality data, including ammonia, chloride, temperature, and pH from existing monitoring wells at both sites from 1983 to 2008, major ions collected monthly from 2006 and 2008, and a synoptic sampling event for stable isotopes, tritium, and dissolved gases in 2008, were used to determine the source of ammonium in groundwater and possible migration pathways. Geochemical modeling was used to determine possible effects of injected wastewater on native water and aquifer matrix geochemistry. Injected wastewater was determined to be the source of elevated ammonium concentrations above ambient water levels, based on the results of major ion concentrations, tritium, dissolved noble gases and 15N isotopes analyses. Various possible fluid migration pathways were identified at the sites. Data for the south site suggest buoyancy-driven vertical pathways to overlying aquifers bypassing the confining units, with little mixing of injected wastewater with native water as it migrated upward. Once it is introduced into an aquifer, the injectate appeared to migrate advectively with the regional groundwater flow. Geochemical modeling indicated that CO 2-enriched injected wastewater allowed for carbonate dissolution along the vertical pathways, enhancing permeability along these flowpaths. At the north site, diffusive upward flow through the confining units or offsite vertical pathways were determined to be possible, however no evidence was detected for any on-site confining unit bypass pathway. No evidence was observed at either site of injected wastewater migration to the Upper Floridan aquifer, which is used as a municipal water supply and for aquifer storage and recovery.

Determination of vertical and horizontal pathways of injected fresh wastewater into a deep saline aquifer (Florida, USA) using natural chemical tracers

Two deep-well injection sites in south Florida, USA, inject an average of 430 million liters per day (MLD) of treated domestic fresh wastewater into a deep saline aquifer 900 m below land surface. Elevated levels of NH3 (highest concentration 939 µmol) in the overlying aquifer above ambient concentrations (concentration less than 30 µmol) were evidence of the upward migration of injected fluids. Three pathways were distinguished based on ammonium, chloride and bromide ratios, and temperature. At the South District Wastewater Treatment Plant, the tracer ratios showed that the injectate remained chemically distinct as it migrated upwards through rapid vertical pathways via density-driven buoyancy. The warmer injectate (mean 28°C) retained the temperature signal as it vertically migrated upwards; however, the temperature signal did not persist as the injectate moved horizontally into the overlying aquifers. Once introduced, the injectate moved slowly horizontally through the aquifer and mixed with ambient water. At the North District Wastewater Treatment Plant, data provide strong evidence of a one-time pulse of injectate into the overlying aquifers due to improper well construction. No evidence of rapid vertical pathways was observed at the North District Wastewater Treatment Plant.

Hydrogeophysical characterization of anisotropy in the Biscayne Aquifer using geophysical methods

The anisotropy of the Biscayne Aquifer which serves as the source of potable water for Miami-Dade County was investigated by applying geophysical methods. Electrical resistivity imaging, self potential and ground penetration radar techniques were employed in both regional and site specific studies. In the regional study, electrical anisotropy and resistivity variation with depth were investigated with azimuthal square array measurements at 13 sites. The observed coefficient of electrical anisotropy ranged from 1.01 to 1.36. The general direction of measured anisotropy is uniform for most sites and trends W-E or SE-NW irrespective of depth. Measured electrical properties were used to estimate anisotropic component of the secondary porosity and hydraulic anisotropy which ranged from 1 to 11% and 1.18 to 2.83 respectively. 1-D sounding analysis was used to models the variation of formation resistivity with depth. Resistivities decreased from NW (close to the margins of the everglades) to SE on the shores of Biscayne Bay. Porosity calculated from Archie's law, ranged from 18 to 61% with higher values found along the ridge. Higher anisotropy, porosities and hydraulic conductivities were on the Atlantic Coastal Ridge and lower values at low lying areas west of the ridge. The cause of higher anisotropy and porosity is attributed to higher dissolution rates of the oolitic facies of the Miami Formation composing the ridge. The direction of minimum resistivity from this study is similar to the predevelopment groundwater flow direction indicated in published modeling studies. Detailed investigations were carried out to evaluate higher anisotropy at West Perrine Park located on the ridge and Snapper Creek Municipal well field where the anisotropy trend changes with depth. The higher anisotropy is attributed to the presence of solution cavities oriented in the E-SE direction on the ridge. Similarly, the change in hydraulic anisotropy at the well field might be related to solution cavities, the surface canal and groundwater extraction wells.^

Trends in Water Quality within the Broward County Portion of the Biscayne Aquifer

Continuous and reliable monitoring of contaminants in drinking water, which adversely affect human health, is the main goal of the Broward County Well Field Protection Program. In this study the individual monitoring station locations were used in a yearly and quarterly spatiotemporal Ordinary Kriging interpolation to create a raster network of contaminant detections. In the final analysis, the raster spatiotemporal nitrate concentration trends were overlaid with a pollution vulnerability index to determine if the concentrations are influenced by a set of independent variables. The pollution vulnerability factors are depth to water, recharge, aquifer media, soil, impact to vadose zone, and conductivity. The creation of the nitrate raster dataset had an average RMS Standardized error close to 1 at 0.98. The greatest frequency of detections and the highest concentrations are found in the months of April, May, June, July, August, and September. An average of 76.4% of the nitrate intersected with cells of the pollution vulnerability index over 100.

Hydrogeophysical Characterization of Anisotropy in the Biscayne Aquifer Using Geophysical Methods

The anisotropy of the Biscayne Aquifer which serves as the source of potable water for Miami-Dade County was investigated by applying geophysical methods. Electrical resistivity imaging, self potential and ground penetration radar techniques were employed in both regional and site specific studies. In the regional study, electrical anisotropy and resistivity variation with depth were investigated with azimuthal square array measurements at 13 sites. The observed coefficient of electrical anisotropy ranged from 1.01 to 1.36. The general direction of measured anisotropy is uniform for most sites and trends W-E or SE-NW irrespective of depth. Measured electrical properties were used to estimate anisotropic component of the secondary porosity and hydraulic anisotropy which ranged from 1 to 11% and 1.18 to 2.83 respectively. 1-D sounding analysis was used to models the variation of formation resistivity with depth. Resistivities decreased from NW (close to the margins of the everglades) to SE on the shores of Biscayne Bay. Porosity calculated from Archie's law, ranged from 18 to 61% with higher values found along the ridge. Higher anisotropy, porosities and hydraulic conductivities were on the Atlantic Coastal Ridge and lower values at low lying areas west of the ridge. The cause of higher anisotropy and porosity is attributed to higher dissolution rates of the oolitic facies of the Miami Formation composing the ridge. The direction of minimum resistivity from this study is similar to the predevelopment groundwater flow direction indicated in published modeling studies. Detailed investigations were carried out to evaluate higher anisotropy at West Perrine Park located on the ridge and Snapper Creek Municipal well field where the anisotropy trend changes with depth. The higher anisotropy is attributed to the presence of solution cavities oriented in the E-SE direction on the ridge. Similarly, the change in hydraulic anisotropy at the well field might be related to solution cavities, the surface canal and groundwater extraction wells.

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

Natural organic matter and colloid-facilitated arsenic transport and transformation in porous soil media

Prediction of arsenic transport and transformation in soil environment requires understanding the transport mechanisms and proper estimation of arsenic partitioning tong all three phases in soil/aquifer systems: mobile colloids, mobile soil solution, and immobile soil solids. The primary purpose of this research is to study natural dissolved organic matter (DOM)/colloid-facilitated transport of arsenic and understand the role of soil derived carriers in the transport and transformation of both inorganic and organoarsenicals in soils. ^ DOM/colloid facilitated arsenic transport and transformation in porous soil media were investigated using a set of experimental approaches including batch experiment, equilibrium membrane dialysis experiment and column experiment. Soil batch experiment was applied to investigate arsenic adsorption on a variety of soils with different characteristics; Equilibrium membrane dialysis was employed to determine the 'free' and 'colloid-bound/complexed' arsenic in water extracts of chosen soils; Column experiments were also set up in the laboratory to simulate arsenic transport and transformation through golf course soils in the presence and absence of soil-derived dissolved substances. ^ The experimental results revealed that organic matter amendments effectively reduced soil arsenic adsorption. The majority of arsenic present in the soil extracts was associated with small substances of molecular weight (MW) between 500 and 3,500 Da, Only a small fraction of arsenic was associated with higher MW substances (MW > 3,500 Da), which was operationally defined as colloidal part in this study. The association of arsenic and DOM in the soil extracts strongly affected arsenic bioavailability, arsenic transport and transformation in soils. The results of column experiments revealed arsenic complicated behavior with various processes occurring in soils studied, including: soil arsenic' adsorption, facilitated arsenic transportation by dissolved substances presented in soil extracts and microorganisms involved arsenic species transformation. ^ Soil organic matter amendments effectively reduce soil arsenic adsorption capability either by scavenging 'soil arsenic adsorption sites or by interactions between arsenic species and dissolved organic chemicals in soil solution. Close attention must be paid for facilitated arsenic transport by dissolved substances presented in soil solution and microorganisms involved arsenic species transformation in arsenic-contaminated soils.^

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.

Larger foraminiferal biostratigraphy, systematics and paleoenvironments of the Avon Park Formation and Ocala Limestone, Highlands County, Florida

This study investigates the use of larger foraminifera in determining the biostratigraphy of the Avon Park Formation and the Ocala Limestone in central Florida. Sedimentary rocks of the Avon Park Formation are the oldest exposed deposits in the state of Florida, and together with the Ocala Limestone comprise a part of the confining unit of the Floridan Aquifer, a major source of Florida's water supply. ^ Material from the ROMP 29A core collected by the U.S. Geological Survey was evaluated and compared to previous studies of the biostratigraphy of the formations. The larger foraminifera of the Avon Park Formation were examined in thin section, and those of the Ocala Limestone were free specimens. The larger foraminifera from both units were described and identified, and the biostratigraphy determined. The morphological features of the larger foraminifera of the Ocala Limestone were measured and analyzed at various depths within the ROMP 29A core.^ The Avon Park Formation contains predominantly the shallow-water, conical foraminifera Fallotella cookei, Fallotella floridana, Pseudochrysalidina floridana, Coleiconus christianaensis, Coleiconus sp. A, Coskinolina sp. A, Coskinolina sp. B, Fallotella sp. A, Fallotella sp. B, Fabularia vaughani and larger miliolids. ^ The Ocala Limestone contains a different, deeper water assemblage that included the larger foraminifera Heterostegina ocalana, Lepidocyclina ocalana varieties, Lepidocyclina chaperi, Lepidocyclina pustulosa, Nummulites willcoxi, Nummulites striatoreticulatus, Nummulites floridensis and Pseudophragmina spp. A, B, and C. The age of the Avon Park Formation was corroborated by the occurrence of the biomarker echinoid Neolaganum dalli as Eocene, and the Ocala Limestone also contained Eocene larger foraminifera with Eocene to possibly Oligocene calcareous nannofossils. The distribution of the larger foraminifera of the Avon Park Formation was correlated with the subtidal and peritidal zones of the continental shelf. Analyses of variance showed that the changes in measurements of the morphology in Heterostegina ocalana, Lepidocyclina spp. and Nummulites spp. were correlated with change in the depositional environments.^

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

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

Chilling damage in a changing climate in coastal landscapes of the subtropical zone: a case study from south Florida

Extensive portions of the southern Everglades are characterized by series of elongated, raised peat ridges and tree islands oriented parallel to the predominant flow direction, separated by intervening sloughs. Tall herbs or woody species are associated with higher elevations and shorter emergent or floating species are associated with lower elevations. The organic soils in this “Ridge-and-Slough” landscape have been stable over millennia in many locations, but degrade over decades under altered hydrologic conditions. We examined soil, pore water, and leaf phosphorus (P) and nitrogen (N) distributions in six Ridge and Slough communities in Shark Slough, Everglades National Park. We found P enrichment to increase and N to decrease monotonically along a gradient from the most persistently flooded sloughs to rarely flooded ridge environments, with the most dramatic change associated with the transition from marsh to forest. Leaf N:P ratios indicated that the marsh communities were strongly P-limited, while data from several forest types suggested either N-limitation or co-limitation by N and P. Ground water stage in forests exhibited a daytime decrease and partial nighttime recovery during periods of surface exposure. The recovery phase suggested re-supply from adjacent flooded marshes or the underlying aquifer, and a strong hydrologic connection between ridge and slough. We therefore developed a simple steady-state model to explore a mechanism by which a phosphorus conveyor belt driven by both evapotranspiration and the regional flow gradient can contribute to the characteristic Ridge and Slough pattern. The model demonstrated that evapotranspiration sinks at higher elevations can draw in low concentration marsh waters, raising local soil and water P concentrations. Focusing of flow and nutrients at the evapotranspiration zone is not strong enough to overcome the regional gradient entirely, allowing the nutrient to spread downstream and creating an elongated concentration plume in the direction of flow. Our analyses suggest that autogenic processes involving the effects of initially small differences in topography, via their interactions with hydrology and nutrient availability, can produce persistent physiographic patterns in the organic sediments of the Everglades.

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.

Adsorption and desorption of phosphate on limestone in experiments simulating seawater intrusion

This study investigates the potential release of from carbonate aquifers exposed to seawater intrusion. Adsorption and desorption of in the presence of deionized water (DIW) and seawater were conducted on a large block of Pleistocene age limestone to simulate the effects of seawater intrusion into a coastal carbonate aquifer at the laboratory scale. The limestone showed strong adsorption of in DIW, while adsorption was significantly less in the presence of seawater. Dissolution of CaCO3 was found to prevent adsorption at salinities less than 30 psu. Adsorption of was limited at higher salinities (30–33 psu), due to competition with ions for adsorption sites. At a salinity3 precipitated. Concentrations of between 2 and 5 μmol/L were released by desorption when the limestone was exposed to seawater. The results of this study suggest that as seawater intrudes into an originally freshwater coastal aquifer, adsorbed may be released into the groundwater. Consequently, adsorbed is expected to be released from coastal carbonate aquifers world-wide as sea level continues to rise exposing more of the freshwater aquifer to seawater.