Assessment of the occurrence and potential effects of pharmaceuticals and personal care products in South Florida waters and sediments

A LLE-GC-MS method was developed to detect PPCPs in surface water samples from Big Cypress National Park, Everglades National Park and Biscayne National Park in South Florida. The most frequently found PPCPs were caffeine, DEET and triclosan with detected maximum concentration of 169 ng/L, 27.9 ng/L and 10.9 ng/L, respectively. The detection frequencies of hormones were less than PPCPs. Detected maximal concentrations of estrone, 17β-estradiol, coprostan-3-ol, coprostane and coprostan-3-one were 5.98 ng/L, 3.34 ng/L, 16.5 ng/L, 13.5 ng/L and 6.79 ng/L, respectively. An ASE-SPE-GC-MS method was developed and applied to the analysis of the sediment and soil area where reclaimed water was used for irrigation. Most analytes were below detection limits, even though some of analytes were detected in the reclaimed water at relatively high concentrations corroborating the fact that PPCPs do not significantly partition to mineral phases. An online SPE-HPLC-APPI-MS/MS method and an online SPE-HPLC-HESI-MS/MS method were developed to analyze reclaimed water and drinking water samples. In the reclaimed water study, reclaimed water samples were collected from the sprinkler for a year-long period at Florida International University Biscayne Bay Campus, where reclaimed water was reused for irrigation. Analysis results showed that several analytes were continuously detected in all reclaimed water samples. Coprostanol, bisphenol A and DEET's maximum concentration exceeded 10 μg/L (ppb). The four most frequently detected compounds were diphenhydramine (100%), DEET (98%), atenolol (98%) and carbamazepine (96%). In the study of drinking water, 54 tap water samples were collected from the Miami-Dade area. The maximum concentrations of salicylic acid, ibuprofen and DEET were 521 ng/L, 301 ng/L and 290 ng/L, respectively. The three most frequently detected compounds were DEET (93%), carbamazepine (43%) and salicylic acid (37%), respectively. Because the source of drinking water in Miami-Dade County is the relatively pristine Biscayne aquifer, these findings suggest the presence of wastewater intrusions into the delivery system or the onset of direct influence of surface waters into the shallow aquifer.

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

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

Sediment Dynamics of a Shallow Hypereutrophic lake: Lake Jesup, Florida, USA

Improved knowledge of sediment dynamics within a lake system is important for understanding lake water quality. This research was focused on an assessment of the vertical sediment flux in Lake Jesup, a shallow (1.3 m average depth) hypereutrophic lake of central Florida. Sediment dynamics were assessed at varying time scales (daily to weekly) to understand the transport of sediments from external forces; wind, waves, precipitation and/or runoff. Four stations were selected within the lake on the basis of water depth and the thicknesses of unconsolidated (floc) and consolidated sediments. At each of these stations, a 10:1 (length to diameter) high aspect ratio trap (STHA) was deployed to collect particulate matter for a one to two week period. The water and sediment samples were collected and analyzed for total carbon (TC), total phosphorus (TP) and total nitrogen (TN). Mass accumulation rates (MAR) collected by the traps varied from 77 to 418 g m-2 d-1 over seven deployments. TN, TP and TC sediment concentrations collected by the traps were consistently higher than the sediments collected by coring the lake bottom and is most likely associated with water column biomass. A yearly nutrient budget was determined from August 2009 to August 2010 with flux calculated as 2,033,882 mt yr-1.

Biogeochemical effects of iron availability on primary producers in a shallow marine carbonate environment

We completed a synoptic survey of iron, phosphorus, and sulfur concentrations in shallow marine carbonate sediments from south Florida. Total extracted iron concentrations typically were 50 μmol g-1 dry weight (DW) and tended to decrease away from the Florida mainland, whereas total extracted phosphorus concentrations mostly were 10 μmol g-1 DW and tended to decrease from west to east across Florida Bay. Concentrations of reduced sulfur compounds, up to 40 μmol g-1 DW, tended to covary with sediment iron concentrations, suggesting that sulfide mineral formation was iron-limited. An index of iron availability derived from sediment data was negatively correlated with chlorophyll a concentrations in surface waters, demonstrating the close coupling of sediment-water column processes. Eight months after applying a surface layer of iron oxide granules to experimental plots, sediment iron, phosphorus, and sulfur were elevated to a depth of 10 cm relative to control plots. Biomass of the seagrass Thalassia testudinum was not different between control and iron addition plots, but individual shoot growth rates were significantly higher in experimental plots after 8 months. Although the iron content of leaf tissues was significantly higher from iron addition plots, no difference in phosphorus content of T. testudinum leaves was observed. Iron addition altered plant exposure to free sulfide, documented by a significantly higher δ34S of leaf tissue from experimental plots relative to controls. Iron as a buffer to toxic sulfides may promote individual shoot growth, but phosphorus availability to plants still appears to limit production in carbonate sediments.

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.

Temporal and Spatial Variability of Mesozooplankton in a Shallow Sub-Tropical Bay: Influence of Top-Down Control

Quantifying the relationship between mesozooplankton and water quality parameters identifies the factors that structure the mesozooplankton community and can be used to generate hypotheses regarding the mechanisms that control the mesozooplankton population and potentially the trophic network. To investigate this relationship, mesozooplankton and water quality data were collected in Florida Bay from 1994 to 2004. Three key characteristics were found in the mesozooplankton community structure: (1) there are significant differences between the four sub-regions of Florida Bay; (2) there is a break in May of 1997 with significant differences before and after this date; and (3) there is a positive correlation between mesozooplankton abundance and salinity. The latter two characteristics are closely correlated with predator abundance, indicating the importance of top-down control. Hypersaline periods appear to provide a refuge from predators, allowing mesozooplankton to increase in abundance despite the increased physiological stress.

The role of recharge and evapotranspiration as hydraulic drivers of ion concentrations in shallow groundwater on Everglades tree islands, Florida (USA)

Recently, evapotranspiration has been hypothesized to promote the secondary formation of calcium carbonate year-round on tree islands in the Everglades by influencing groundwater ions concentrations. However, the role of recharge and evapotranspiration as drivers of shallow groundwater ion accumulation has not been investigated. The goal of this study is to develop a hydrologic model that predicts the chloride concentrations of shallow tree island groundwater and to determine the influence of overlying biomass and underlying geologic material on these concentrations. Groundwater and surface water levels and chloride concentrations were monitored on eight constructed tree islands at the Loxahatchee Impoundment Landscape Assessment (LILA) from 2007 to 2010. The tree islands at LILA were constructed predominately of peat, or of peat and limestone, and were planted with saplings of native tree species in 2006 and 2007. The model predicted low shallow groundwater chloride concentrations when inputs of regional groundwater and evapotranspiration-to-recharge rates were elevated, while low evapotranspiration-to-recharge rates resulted in a substantial increase of the chloride concentrations of the shallow groundwater. Modeling results indicated that evapotranspiration typically exceeded recharge on the older tree islands and those with a limestone lithology, which resulted in greater inputs of regional groundwater. A sensitivity analysis indicated the shallow groundwater chloride concentrations were most sensitive to alterations in specific yield during the wet season and hydraulic conductivity in the dry season. In conclusion, the inputs of rainfall, underlying hydrologic properties of tree islands sediments and forest structure may explain the variation in ion concentration seen across Everglades tree islands.

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

Spatial and Temporal Changes in Groundwater Salinity in South Florida

A combination of statistical and interpolation methods and Geographic Information System (GIS) spatial analysis was used to evaluate the spatial and temporal changes in groundwater Cl− concentrations in Collier and Lee Counties (southwestern Florida), and Miami-Dade and Broward Counties (southeastern Florida), since 1985. In southwestern Florida, the average Cl− concentrations in the shallow wells (0–43 m) in Collier and Lee Counties increased from 132 mg L−1 in 1985 to 230 mg L−1 in 2000. The average Cl− concentrations in the deep wells (>43 m) of southwestern Florida increased from 392 mg L−1 in 1985 to 447 mg L−1 in 2000. Results also indicated a positive correlation between the mean sea level and Cl− concentrations and between the mean sea level and groundwater levels for the shallow wells. Concentrations in the Biscayne Aquifer (southeastern Florida) were significantly higher than those of southwestern Florida. The average Cl− concentrations increased from 159 mg L−1 in 1985 to 470 mg L−1 in 2010 for the shallow wells (<33 m) and from 1360 mg L−1 in 1985 to 2050 mg L−1 in 2010 for the deep wells (>33 m). In the Biscayne Aquifer, wells showed a positive or negative correlation between mean sea level and Cl− concentrations according to their location with respect to the saltwater intrusion line. Wells located inland behind canal control structures and west of the saltwater intrusion line showed negative correlation values, whereas wells located east of the saltwater intrusion line showed positive values. Overall, the results indicated that since 1985, there was a potential decline in the available freshwater resources estimated at about 12–17% of the available drinking-quality groundwater of the southeastern study area located in the Biscayne Aquifer.

The Role of Recharge and Evapotranspiration as Hydraulic Drivers of Ion Concentrations in Shallow Groundwater on Everglades Tree Islands, Florida (USA)

Recently, evapotranspiration has been hypothesized to promote the secondary formation of calcium carbonate year-round on tree islands in the Everglades by influencing groundwater ions concentrations. However, the role of recharge and evapotranspiration as drivers of shallow groundwater ion accumulation has not been investigated. The goal of this study is to develop a hydrologic model that predicts the chloride concentrations of shallow tree island groundwater and to determine the influence of overlying biomass and underlying geologic material on these concentrations. Groundwater and surface water levels and chloride concentrations were monitored on eight constructed tree islands at the Loxahatchee Impoundment Landscape Assessment (LILA) from 2007 to 2010. The tree islands at LILA were constructed predominately of peat, or of peat and limestone, and were planted with saplings of native tree species in 2006 and 2007. The model predicted low shallow groundwater chloride concentrations when inputs of regional groundwater and evapotranspiration-to-recharge rates were elevated, while low evapotranspiration-to-recharge rates resulted in a substantial increase of the chloride concentrations of the shallow groundwater. Modeling results indicated that evapotranspiration typically exceeded recharge on the older tree islands and those with a limestone lithology, which resulted in greater inputs of regional groundwater. A sensitivity analysis indicated the shallow groundwater chloride concentrations were most sensitive to alterations in specific yield during the wet season and hydraulic conductivity in the dry season. In conclusion, the inputs of rainfall, underlying hydrologic properties of tree islands sediments and forest structure may explain the variation in ion concentration seen across Everglades tree islands.

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.