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.

In southeastern Florida, the majority of drinking water comes from the Biscayne aquifer. This aquifer is comprised of heterogeneous limestones, sandstones, sand, shell and clayey sand with zones of very high permeability. Visualizing the spatial variations in lithology, porosity and permeability of heterogeneous aquifers, like the Biscayne, can be difficult using traditional methods of investigation. Using the Roxar IRAP RMS software multi-layered 3D conceptual geomodels of the lithology, cyclostratigraphy and porosity were created in a portion of the Biscayne aquifer. The models were built using published data from borehole geophysical measurements, core samples, and thin sections. Spatial relations between lithology, cyclostratigraphy, porosity, and preferential flow zones were compared and contrasted to better understand how these geologic features were inter-related. The models show local areas of differing porosity within and cross-cutting different cycles and lithologies. Porosity in the Biscayne aquifer study area follows a hierarchy attributed to lithofacies with a pattern of increasing porosity for the high frequency cycles. This modeling improves understanding of the distribution and interconnectedness of preferential flow zones, and is thus an invaluable tool for future studies of groundwater flow and groundwater contamination in the Biscayne aquifer. / Includes bibliography. / Thesis (M.S.)--Florida Atlantic University, 2014. / FAU Electronic Theses and Dissertations Collection

The karst Biscayne aquifer is characterized by a heterogeneous spatial arrangement of porosity, making hydrogeological characterization difficult. In this dissertation, I investigate the use of ground penetrating radar (GPR), for understanding the spatial distribution of porosity variability in the Miami Limestone presented as a compilation of studies where scale of measurement is progressively increased to account for varying dimensions of dissolution features. In Chapter 2, GPR in zero offset acquisition mode is used to investigate the 2-D distribution of porosity and dielectric permittivity in a block of Miami Limestone at the laboratory scale (< 1.0 m). Petrophysical models based on fully saturated and unsaturated. water conditions are used to estimate porosity and solid dielectric permittivity of the limestone. Results show a good correspondence between analytical and GPR-based porosity estimates and show variability between 22.0-66.0 %. In Chapter 3, GPR in common offset and common midpoint acquisition mode are used to estimate bulk porosity of the unsaturated Miami Limestone at the field scale (10.0-100.0 m). Estimates of porosity are based on the assumption that the directly measured water table reflector is flat and that any deviation is attributed to changes in velocity due to porosity variability. Results show sharp changes in porosity ranging between 33.2-60.9 % attributed to dissolution areas. In Chapter 4, GPR in common offset mode is used to characterize porosity variability in the saturated Biscayne aquifer at 100-1000 m field scales. The presence of numerous diffraction hyperbolae are used to estimate electromagnetic wave velocity and asses both horizontal and vertical changes in porosity after application of a petrophysical model. Results show porosity variability between 23.0-41.0 % and confirm the presence of isolated areas that could serve as enhanced infiltration or recharge. This research allows for the identification and delineation areas of macroporosity areas at 0.01 m lateral resolution and shows variability of porosity at different scales, reaching 37.0 % within 1.3 m, associated with areas of enhanced dissolution. Such improved resolution of porosity estimates can benefit water management efforts and transport modelling and help to better understand small scale relationships between ground water and surface water interactions. / Includes bibliography. / Dissertation (Ph.D.)--Florida Atlantic University, 2014. / FAU Electronic Theses and Dissertations Collection

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.

The Biscayne Aquifer is a primary source of water supply in Southeast Florida. As a coastal aquifer, it is threatened by saltwater intrusion (SWI) when the natural groundwater flow is altered by over-pumping of groundwater. SWI is detrimental to the quality of fresh groundwater sources, making the water unfit for drinking due to mixing and reactions with aquifer minerals. Increasing water demand and complex environmental issues thus force water utilities in South Florida to sustainably manage saltwater intrusion and develop alternative water supplies (e.g., aquifer storage and recovery, ASR). The objectives of this study were to develop and use calibrated geochemical models to estimate water quality changes during saline intrusion and during ASR in south Florida. A batch-reaction model of saltwater intrusion was developed and important geochemical reactions were inferred. Additionally, a reactive transport model was developed to assess fate and transport of major ions and trace metals (Fe, As) at the Kissimmee River ASR. Finally, a cost-effective management of saltwater intrusion that involves using abstraction and recharge wells was implemented and optimized for the case of the Biscayne Aquifer. Major processes in the SWI areas were found to be mixing and dissolution-precipitation reactions with calcite and dolomite. Most of the major ions (Cl, Na, K, Mg, SO4) behaved conservatively during ASR while Ca and alkalinity were affected by carbonate reactions and cation exchange. A complex set of reactions involving thermodynamic equilibrium, kinetics and surface complexation reactions was required in the ASR model to simulate observed concentrations of Fe and As. The saltwater management model aimed at finding optimal locations and flow rates for abstraction and recharge wells. Optimal solutions (i.e., minimum total salt and total cost Pareto front) were produced for the Biscayne Aquifer for scenarios of surface recharge induced by climate change-affected precipitation. In general, abstraction at the maximum rate near the coast and artificial recharge at locations much further inland were found to be optimal. Knowledge developed herein directly supports the understanding of SWI caused by anthropogenic stressors, such as over-pumping and sea level rise, on coastal aquifers.