Geology topics

Kevin J. Cunningham

Publications and source records attributed to Kevin J. Cunningham.

At least 19 recordsLinked to original sources

USGS2024

Formation of vertical columnar seismic structures and seafloor depressions by groundwater discharge in the drowned Miami Terrace platform and overlying deep-water carbonates, southeastern Florida

The presence of vertical cross-formational fluid migration passageways within sedimentary basins can profoundly impact aquifer and reservoir fluid-flow and their identification is fundamental to informing management of subsurface fluid resources (groundwater, oil, gas). In an onshore and offshore southeastern part of Florida, 2D/3D seismic-reflection and bathymetry data document ∼153 vertical columnar structures composed of reflection disruptions up to 790 m in the height and averaging 360 m in diameter, and ∼219 subcircular to circular seafloor depressions up to 1334 m wide. Our study focuses on these features found within the offshore shallow-marine carbonate Miami Terrace platform, which drowned approximately at the end of the middle Miocene, and within overlying Plio-Quaternary deep-water carbonate slope and drift deposits. Most columnar structures are rooted in stratiform aquifers of the Miami Terrace platform and associated with faults or fault intersections produced by Eocene and circa late Miocene tectonics. The columns commonly terminate within the platform or as subcircular depressions along an amalgamated karstic and drowning unconformity at the platform top. The columns typically stretch upwards from a zone of deep karst cavity collapse through the Miami Terrace platform with upward decreasing sag on internal reflections. Following drowning and Plio-Quaternary partial burial of the Miami Terrace platform by deep-water deposits, the subcircular depressions and faults along the platform top were points of origin for a second phase of column growth upward into the deep-water deposits. The continuation of deep platform cavity collapse and column evolution produced pockmarks along paleo-seafloors within the deep-water deposits and at the present-day sea floor. The Plio-Quaternary pockmarks formed at water depths too deep to suggest an origin related to meteoric karst above or near sea level, but rather their formation is suggested to be related to cyclic sea level falls that drove increased groundwater head and density gradients, and seafloor discharge of offshore freshened groundwater sourced from the underlying platform. Plausibly, mixing of freshened groundwater and seawater at the seafloor discharge sites drove dissolution of the host deep-water deposits, which together with erosion by groundwater venting and current scouring formed the pockmarks. Seaward of the Plio-Quaternary seafloor pockmarks, at the late-middle Miocene upper slope of the Miami Terrace platform and along the regional karst/drowning unconformity is a slope-parallel band of ∼189 densely distributed subcircular seafloor depressions with diameters up to 1334 m at water depths up to ∼660 m. It is plausible that along the upper slope, faults and fractures produced by gravity-driven slope instability and possibly tectonics formed a dense network of fluid passageways that promoted upward artesian freshened groundwater flow to sites of discharge where mixing with seawater generated limestone dissolution and the depressions. But tectonic uplift may have forced emersion and initial meteoric sinkhole formation circa late Miocene with later enhancement by freshened groundwater discharge and bottom current erosion.

Florida↗

USGS2022

New larger benthic foraminifera from the subsurface Lower to Middle Eocene Oldsmar Formation of southeastern Florida (USA)

We describe two larger benthic foraminiferal taxa collected from wells drilled in the subsurface Eocene rocks of southeastern Florida that are new to peninsular Florida and the Caribbean region. Saudia floridana n.sp. is characteristic of a foraminiferal assemblage, along with Helicostegina gyralis , wide forms of the Cushmania americana group, and Gunteria floridana , in an upper part of the Oldsmar Formation. Globogypsinoides browardensis n.gen. n.sp. occurs in a second foraminiferal assemblage, along with Borelis cf. floridanus , Coskinolina cf. yucatanensis , and as-yet undescribed large rotaliids, in a middle part of the Oldsmar Formation. The foraminiferal assemblage of the middle Oldsmar unit is ascribed an Ypresian age and the assemblage of the upper Oldsmar unit a Lutetian age. These two assemblages indicate inner shelf water depths of 40 m or less on the Florida Platform during the Early to Middle Eocene deposition of the middle to upper part of the Oldsmar Formation.

Florida↗

USGS2019

Hydrogeologic characterization of part of the Lower Floridan aquifer at the South District Wastewater Treatment Plant, Miami-Dade County, Florida

The South District Wastewater Treatment Plant in southeastern Miami-Dade County, Florida, includes a Class I treated wastewater injection well system. The detection of ammonia in monitoring zones above the injection zone in the Lower Floridan aquifer has elicited a need to understand the nature of confinement within the Lower Floridan aquifer as it pertains to the vertical migration of injectate out of the injection zone upward into the Underground Source of Drinking Water in the upper part of the Floridan aquifer system. Geologic and geophysical data, borehole video imagery, and aquifer performance data were used to refine and clarify the geologic and hydrogeologic frameworks of part of the Lower Floridan aquifer at the treatment plant. The data provide evidence for zones of enhanced dissolution permeability, extensive secondary porosity, fractures, karst collapse structures, and faults that could provide vertical cross-formational fluid pathways that transect the Lower Floridan aquifer.

Florida↗

USGS2018

Three-dimensional seismic characterization of karst in the Floridan aquifer system, southeastern Miami-Dade County, Florida

Two three-dimensional seismic surveys totaling 3.4 square miles were acquired in southeastern Miami-Dade County during 2015 as part of an ongoing broad regional investigation by the U.S. Geological Survey, in cooperation with the Miami-Dade Water and Sewer Department, that includes mapping and karst characterization of the Floridan aquifer system in southeastern Florida. Twenty columniform seismic-sag structures were identified in the three-dimensional seismic data acquired at the South Miami Heights study area in Miami-Dade County, Florida. The seismic-sag structures are commonly composed of a lower and upper seismic facies, where the lower facies is interpreted to be a high permeability column-shaped volume of intense karstification composed of faults, fractures, and collapsed caves, and the upper facies is interpreted to be a lower permeability column-shaped volume composed of sagging suprastratal reflections having little or no faulting, fracturing, or caves. The seismic-sag structures are columniform karst-collapse structures that formed by cave collapse or stoping (upward void migration caused primarily by the mechanical process of ceiling collapse) or both in carbonate rocks of the early Eocene Oldsmar Formation and in some cases in the lowermost part of the middle Eocene Avon Park Formation. Columniform subsidence and sagging of overburden succeeded karst collapse in the lower part of the structures. At the study area, there may be relatively higher potential for the columniform karst-collapse structures to provide cross-formational fluid migration upward from the Boulder Zone into the lower part of middle confining unit 2 as compared to a lower potential for cross-formational fluid migration from the upper part of middle confining unit 2 upward to the Upper Floridan aquifer.

Florida↗

USGS2018

Statistical approach to neural network imaging of karst systems in 3D seismic reflection data

The current lack of a robust, standardized technique for geophysical mapping of karst systems can be attributed to both the complexity of the environment and prior technological limitations. Abrupt lateral variations in physical properties that are inherent to karst systems generate significant geophysical noise, challenging conventional seismic signal processing and interpretation. Modern application of neural networks to multi-attribute seismic interpretation now provide a semiautomated method for identifying and leveraging the nonlinear relationships exhibited among seismic attributes. The ambiguity generally associated with designing neural networks for seismic object detection can be reduced via statistical analysis of the extracted attribute data. A data-driven approach to selecting the appropriate set of input seismic attributes, as well as the locations and minimum number of training examples, provides a more objective and computationally efficient method for identifying karst systems using reflection seismology. This statistically optimized neural network technique is thoroughly demonstrated using three-dimensional seismic reflection data collected from the southeastern portion of the Florida carbonate platform. Several dimensionality reduction methods are applied and the resulting karst probability models are evaluated relative to one another based on both quantitative and qualitative criteria. Comparing the preferred model, using quadratic discriminant analysis, to previously available seismic object detection workflows demonstrates the karst-specific nature of the tool. Results suggest that the karst multi-attribute workflow presented is capable of approximating the structural boundaries of karst systems with more accuracy and efficiency than a human counterpart or previously presented seismic interpretation schemes. This objective technique, using solely three-dimensional seismic reflection data, likely represents the most practical approach to mapping karst systems for subsequent hydrogeological modeling.

Interpretation↗

USGS2017

Sequence stratigraphy, seismic stratigraphy, and seismic structures of the lower intermediate confining unit and most of the Floridan aquifer system, Broward County, Florida

Deep well injection and disposal of treated wastewater into the highly transmissive saline Boulder Zone in the lower part of the Floridan aquifer system began in 1971. The zone of injection is a highly transmissive hydrogeologic unit, the Boulder Zone, in the lower part of the Floridan aquifer system. Since the 1990s, however, treated wastewater injection into the Boulder Zone in southeastern Florida has been detected at three treated wastewater injection utilities in the brackish upper part of the Floridan aquifer system designated for potential use as drinking water. At a time when usage of the Boulder Zone for treated wastewater disposal is increasing and the utilization of the upper part of the Floridan aquifer system for drinking water is intensifying, there is an urgency to understand the nature of cross-formational fluid flow and identify possible fluid pathways from the lower to upper zones of the Floridan aquifer system. To better understand the hydrogeologic controls on groundwater movement through the Floridan aquifer system in southeastern Florida, the U.S. Geological Survey and the Broward County Environmental Planning and Community Resilience Division conducted a 3.5-year cooperative study from July 2012 to December 2015. The study characterizes the sequence stratigraphy, seismic stratigraphy, and seismic structures of the lower part of the intermediate confining unit aquifer and most of the Floridan aquifer system. Data obtained to meet the study objective include 80 miles of high-resolution, two-dimensional (2D), seismic-reflection profiles acquired from canals in eastern Broward County. These profiles have been used to characterize the sequence stratigraphy, seismic stratigraphy, and seismic structures in a 425-square-mile study area. Horizon mapping of the seismic-reflection profiles and additional data collection from well logs and cores or cuttings from 44 wells were focused on construction of three-dimensional (3D) visualizations of eight sequence stratigraphic cycles that compose the Eocene to Miocene Oldsmar, Avon Park, and Arcadia Formations. The mapping of these seismic-reflection and well data has produced a refined Cenozoic sequence stratigraphic, seismic stratigraphic, and hydrogeologic framework of southeastern Florida. The upward transition from the Oldsmar Formation to the Avon Park Formation and the Arcadia Formation embodies the evolution from (1) a tropical to subtropical, shallow-marine, carbonate platform, represented by the Oldsmar and Avon Park Formations, to (2) a broad, temperate, mixed carbonate-siliciclastic shallow marine shelf, represented by the lower part of the Arcadia Formation, and to (3) a temperate, distally steepened carbonate ramp represented by the upper part of the Arcadia Formation. In the study area, the depositional sequences and seismic sequences have a direct correlation with hydrogeologic units. The approximate upper boundary of four principal permeable units of the Floridan aquifer system (Upper Floridan aquifer, Avon Park permeable zone, uppermost major permeable zone of the Lower Floridan aquifer, and Boulder Zone) have sequence stratigraphic and seismic-reflection signatures that were identified on cross sections, mapped, or both, and therefore the sequence stratigraphy and seismic stratigraphy were used to guide the development of a refined spatial representation of these hydrogeologic units. In all cases, the permeability of the four permeable units is related to stratiform megaporosity generated by ancient dissolution of carbonate rock associated with subaerial exposure and unconformities at the upper surfaces of carbonate depositional cycles of several hierarchical scales ranging from high-frequency cycles to depositional sequences. Additionally, interparticle porosity also contributes substantially to the stratiform permeability in much of the Upper Floridan aquifer. Information from seismic stratigraphy allowed 3D geomodeling of hydrogeologic units—an approach never before applied to this area. Notably, the 3D geomodeling provided 3D visualizations and geocellular models of the depositional sequences, hydrostratigraphy, and structural features. The geocellular data could be used to update the hydrogeologic structure inherent to groundwater flow simulations that are designed to address the sustainability of the water resources of the Floridan aquifer system. Two kinds of pathways that could enable upward cross-formational flow of injected treated wastewater from the Boulder Zone have been identified in the 80 miles of high-resolution seismic data collected for this study: a near-vertical reverse fault and karst collapse structures. The single reverse fault, inferred to be of tectonic origin, is in extreme northeastern Broward County and has an offset of about 19 feet at the level of the Arcadia Formation. Most of the 17 karst collapse structures identified manifest as columniform, vertically stacked sagging seismic reflections that span early Eocene to Miocene age rocks equivalent to much of the Floridan aquifer system and the lower part of the overlying intermediate confining unit. In some cases, the seismic-sag structures extend upward into strata of Pliocene age. The seismic-sag structures are interpreted to have a semicircular shape in plan view on the basis of comparison to (1) other seismic-sag structures in southeastern Florida mapped with two 2D seismic cross lines or 3D data, (2) comparison to these structures located in other carbonate provinces, and (3) plausible extensional ring faults detected with multi-attribute analysis. The seismic-sag structures in the study area have heights as great as 2,500 vertical feet, though importantly, one spans about 7,800 feet. Both multi-attribute analysis and visual detection of offset of seismic reflections within the seismic-sag structures indicate faults and fractures are associated with many of the structures. Multi-attribute analysis highlighting chimney fluid pathways also indicates that the seismic-sag structures have a high probability for potential vertical cross-formational fluid flow along the faulted and fractured structures. A collapse of the seismic-sag structures within a deep burial setting evokes an origin related to hypogenic karst processes by ascending flow of subsurface fluids. In addition, paleo-epigenic karst related to major regional subaerial unconformities within the Florida Platform generated collapse structures (paleo-sinkholes) that are much smaller in scale than the cross-formational seismic-sag structures.

Florida↗

USGS2017

Lithofacies and sequence stratigraphic description of the upper part of the Avon Park Formation and the Arcadia Formation in U.S. Geological Survey G–2984 test corehole, Broward County, Florida

Rock core and sediment from U.S. Geological Survey test corehole G–2984 completed in 2011 in Broward County, Florida, provide an opportunity to improve the understanding of the lithostratigraphic, sequence stratigraphic, and hydrogeologic framework of the intermediate confining unit and Floridan aquifer system in southeastern Florida. A multidisciplinary approach including characterization of sequence stratigraphy, lithofacies, ichnology, foraminiferal paleontology, depositional environments, porosity, and permeability was used to describe the geologic samples from this test corehole. This information has produced a detailed characterization of the lithofacies and sequence stratigraphy of the upper part of the middle Eocene Avon Park Formation and Oligocene to middle Miocene Arcadia Formation. This enhancement of the knowledge of the sequence stratigraphic framework is especially important, because subaerial karst unconformities at the upper boundary of depositional cycles at various hierarchical scales are commonly associated with secondary porosity and enhanced permeability in the Floridan aquifer system.

Florida↗

USGS2015

Geostatistical borehole image-based mapping of karst-carbonate aquifer pores

Quantification of the character and spatial distribution of porosity in carbonate aquifers is important as input into computer models used in the calculation of intrinsic permeability and for next-generation, high-resolution groundwater flow simulations. Digital, optical, borehole-wall image data from three closely spaced boreholes in the karst-carbonate Biscayne aquifer in southeastern Florida are used in geostatistical experiments to assess the capabilities of various methods to create realistic two-dimensional models of vuggy megaporosity and matrix-porosity distribution in the limestone that composes the aquifer. When the borehole image data alone were used as the model training image, multiple-point geostatistics failed to detect the known spatial autocorrelation of vuggy megaporosity and matrix porosity among the three boreholes, which were only 10 m apart. Variogram analysis and subsequent Gaussian simulation produced results that showed a realistic conceptualization of horizontal continuity of strata dominated by vuggy megaporosity and matrix porosity among the three boreholes.

Groundwater↗

USGS2015

Seismic-sequence stratigraphy and geologic structure of the Floridan aquifer system near "Boulder Zone" deep wells in Miami-Dade County, Florida

The U.S. Geological Survey, in cooperation with the Miami-Dade Water and Sewer Department, acquired, processed, and interpreted seismic-reflection data near the North and South District “Boulder Zone” Well Fields to determine if geologic factors may contribute to the upward migration of injected effluent into that upper part of the Floridan aquifer system designated by the U.S. Environmental Protection Agency as an underground source of drinking water. The depth of the Boulder Zone at the North and South District “Boulder Zone” Well Fields ranges from about 2,750 to 3,300 feet below land surface (ft bls), whereas overlying permeable zones used as alternative drinking water supply range in depth from about 825 to 1,580 ft bls at the North and South District “Boulder Zone” Well Fields. Seismic-sequence stratigraphy and geologic structures imaged on seismic-reflection profiles created for the study describe the part of the Floridan aquifer system overlying and within the Boulder Zone. Features of the Floridan aquifer system underlying the Boulder Zone were not studied because seismic-reflection profiles acquired near the North and South District “Boulder Zone” Well Fields lacked adequate resolution at such depths. Stratigraphic analysis of seismic-reflection data collected from the study area was mainly applied to the Floridan aquifer system and used to identify four provisional seismic sequences, which extend vertically from near the base of the Floridan aquifer system upward to the lower part of the intermediate confining unit. These four seismic sequences compose a framework in which each sequence includes a major permeable unit of the Floridan aquifer system; from shallowest to deepest, these units are the Upper Floridan aquifer, Avon Park permeable zone, uppermost major permeable zone of the Lower Floridan aquifer, and Boulder Zone. The relations between seismic-sequence stratigraphy and hydrostratigraphy allow for detailed mapping of permeable zones and semiconfining units of the Floridan aquifer system at a level of resolution never before accomplished using well data alone. In addition to the preceding seismic-reflection analysis, interpretation of geophysical well log data from four effluent injection wells at the North District “Boulder Zone” Well Field delineated a narrow karst collapse structure beneath the injection facility that extends upward about 900 ft from the top of the Boulder Zone to about 125 ft above the top of the uppermost major permeable zone of the Lower Floridan aquifer. No karst collapse structures were identified in the seismic-reflection profiles acquired near the North District “Boulder Zone” Well Field. However, karst collapse structures at the level of the lowermost major permeable zone of the Lower Floridan aquifer at the South District “Boulder Zone” Well Field are present at three locations, as indicated by seismic-reflection data acquired in the C–1 Canal bordering the south side of the injection facility. Results from the North District “Boulder Zone” Well Field well data indicate that a plausible hydraulic connection between faults and stratiform permeability zones may contribute to the upward transport of effluent, terminating above the base of the deepest U.S. Environmental Protection Agency designated underground source of drinking water at the North District “Boulder Zone” Well Field.

Florida↗

USGS2014

Lattice Boltzmann methods applied to large-scale three-dimensional virtual cores constructed from digital optical borehole images of the karst carbonate Biscayne aquifer in southeastern Florida

Digital optical borehole images at approximately 2 mm vertical resolution and borehole caliper data were used to create three-dimensional renderings of the distribution of (1) matrix porosity and (2) vuggy megaporosity for the karst carbonate Biscayne aquifer in southeastern Florida. The renderings based on the borehole data were used as input into Lattice Boltzmann methods to obtain intrinsic permeability estimates for this extremely transmissive aquifer, where traditional aquifer test methods may fail due to very small drawdowns and non-Darcian flow that can reduce apparent hydraulic conductivity. Variogram analysis of the borehole data suggests a nearly isotropic rock structure at lag lengths up to the nominal borehole diameter. A strong correlation between the diameter of the borehole and the presence of vuggy megaporosity in the data set led to a bias in the variogram where the computed horizontal spatial autocorrelation is strong at lag distances greater than the nominal borehole size. Lattice Boltzmann simulation of flow across a 0.4 × 0.4 × 17 m (2.72 m 3 volume) parallel-walled column of rendered matrix and vuggy megaporosity indicates a high hydraulic conductivity of 53 m s −1 . This value is similar to previous Lattice Boltzmann calculations of hydraulic conductivity in smaller limestone samples of the Biscayne aquifer. The development of simulation methods that reproduce dual-porosity systems with higher resolution and fidelity and that consider flow through horizontally longer renderings could provide improved estimates of the hydraulic conductivity and help to address questions about the importance of scale.

Water Resources Research↗

USGS2014

Characterization of the porosity distribution in the upper part of the karst Biscayne aquifer using common offset ground penetrating radar, Everglades National Park, Florida

The karst Biscayne aquifer is characterized by a heterogeneous spatial arrangement of porosity and hydraulic conductivity, making conceptualization difficult. The Biscayne aquifer is the primary source of drinking water for millions of people in south Florida; thus, information concerning the distribution of karst features that concentrate the groundwater flow and affect contaminant transport is critical. The principal purpose of the study was to investigate the ability of two-dimensional ground penetrating radar (GPR) to rapidly characterize porosity variability in the karst Biscayne aquifer in south Florida. An 800-m-long GPR transect of a previously investigated area at the Long Pine Key Nature Trail in Everglades National Park, collected in fast acquisition common offset mode, shows hundreds of diffraction hyperbolae. The distribution of diffraction hyperbolae was used to estimate electromagnetic (EM) wave velocity at each diffraction location and to assess both horizontal and vertical changes in velocity within the transect. A petrophysical model (complex refractive index model or CRIM) was used to estimate total bulk porosity. A set of common midpoint surveys at selected locations distributed along the common-offset transect also were collected for comparison with the common offsets and were used to constrain one-dimensional (1-D) distributions of porosity with depth. Porosity values for the saturated Miami Limestone ranged between 25% and 41% for common offset GPR surveys, and between 23% and 39% for common midpoint GPR surveys. Laboratory measurements of porosity in five whole-core samples from the saturated part of the aquifer in the study area ranged between 7.1% and 41.8%. GPR estimates of porosity were found to be valid only under saturated conditions; other limitations are related to the vertical resolution of the GPR signal and the volume of the material considered by the measurement methodology. Overall, good correspondence between GPR estimates and the direct porosity values from the whole-core samples confirms the ability of GPR common offset surveys to provide rapid characterization of porosity variability in the Biscayne aquifer. The common offset survey method has several advantages: (1) improved time efficiency in comparison to other GPR acquisition modes such as common midpoints; and (2) enhanced lateral continuity of porosity estimates, particularly when compared to porosity measurements on 1-D samples such as rock cores. The results also support the presence of areas of low EM wave velocity or high porosity under saturated conditions, causing velocity pull-down areas and apparent sag features in the reflection record. This study shows that GPR can be a useful tool for improving understanding of the petrophysical properties of highly heterogeneous systems such as karst aquifers, and thus may assist with the development of more accurate groundwater flow models, such as those used for restoration efforts in the Everglades.

Florida↗

USGS2014

Integration of seismic-reflection and well data to assess the potential impact of stratigraphic and structural features on sustainable water supply from the Floridan aquifer system, Broward County, Florida

The U.S. Geological Survey and Broward County water managers commenced a 3.5-year cooperative study in July 2012 to refine the geologic and hydrogeologic framework of the Floridan aquifer system (FAS) in Broward County. A lack of advanced stratigraphic knowledge of the physical system and structural geologic anomalies (faults and fractures originating from tectonics and karst-collapse structures) within the FAS pose a risk to the sustainable management of the resource. The principal objective of the study is to better define the regional stratigraphic and structural setting of the FAS in Broward County. The objective will be achieved through the acquisition, processing, and interpretation of new seismic-reflection data along several canals in Broward County. The interpretation includes integration of the new seismic-reflection data with existing seismic-reflection profiles along Hillsboro Canal in Broward County and within northeast Miami-Dade County, as well as with data from nearby FAS wellbores. The scope of the study includes mapping the geologic, hydrogeologic, and seismic-reflection framework of the FAS, and identifying stratigraphic and structural characteristics that could either facilitate or preclude the sustainable use of the FAS as an alternate water supply or a treated effluent repository. In addition, the investigation offers an opportunity to: (1) improve existing groundwater flow models, (2) enhance the understanding of the sensitivity of the groundwater system to well-field development and upconing of saline fluids, and (3) support site selection for future FAS projects, such as Class I wells that would inject treated effluent into the deep Boulder Zone.

Florida↗

USGS2014

Origins and delineation of saltwater intrusion in the Biscayne aquifer and changes in the distribution of saltwater in Miami-Dade County, Florida

Intrusion of saltwater into parts of the shallow karst Biscayne aquifer is a major concern for the 2.5 million residents of Miami-Dade County that rely on this aquifer as their primary drinking water supply. Saltwater intrusion of this aquifer began when the Everglades were drained to provide dry land for urban development and agriculture. The reduction in water levels caused by this drainage, combined with periodic droughts, allowed saltwater to flow inland along the base of the aquifer and to seep directly into the aquifer from the canals. The approximate inland extent of saltwater was last mapped in 1995. An examination of the inland extent of saltwater and the sources of saltwater in the aquifer was completed during 2008–2011 by using (1) all available salinity information, (2) time-series electromagnetic induction log datasets from 35 wells, (3) time-domain electromagnetic soundings collected at 79 locations, (4) a helicopter electromagnetic survey done during 2001 that was processed, calibrated, and published during the study, (5) cores and geophysical logs collected from 8 sites for stratigraphic analysis, (6) 8 new water-quality monitoring wells, and (7) analyses of 69 geochemical samples. The results of the study indicate that as of 2011 approximately 1,200 square kilometers (km 2 ) of the mainland part of the Biscayne aquifer were intruded by saltwater. The saltwater front was mapped farther inland than it was in 1995 in eight areas totaling about 24.1 km 2 . In many of these areas, analyses indicated that saltwater had encroached along the base of the aquifer. The saltwater front was mapped closer to the coast than it was in 1995 in four areas totaling approximately 6.2 km 2 . The changes in the mapped extent of saltwater resulted from improved spatial information, actual movement of the saltwater front, or a combination of both. Salinity monitoring in some of the canals in Miami-Dade County between 1988 and 2010 indicated influxes of saltwater, with maximum salinities ranging from 1.4 to 32 practical salinity units (PSU) upstream of the salinity control structures. Time-series electromagnetic induction log data from monitoring wells G–3601, G–3608, and G–3701, located adjacent to the Biscayne, Snapper Creek, and Black Creek Canals, respectively, and upstream of the salinity control structures, indicated shallow influxes of conductive water in the aquifer that likely resulted from leakage of brackish water or saltwater from these canals. The determination that saltwater influxes were recent is supported by the similarity in the oxygen and hydrogen stable isotope composition in samples from the Snapper Creek Canal, 1.6 kilometers (km) inland of a salinity control structure, and in samples from well G–3608, which is adjacent to the canal, as well as by the relative ages of the water sampled from well G–3608 and other wells open to the aquifer below the saltwater interface. Historical and recent salinity information from the Card Sound Road Canal, monitoring well FKS8 located adjacent to the canal, and the 2001 helicopter electromagnetic survey indicated that saltwater may occasionally leak from this canal as far inland as 15 km. This leakage may be prevented or reduced by a salinity control structure that was installed in May 2010. Saltwater also may have leaked from the Princeton Canal. Results of geochemical sampling and analysis indicate a close correspondence between droughts and saltwater intrusion. Tritium/helium-3 apparent (piston-flow) ages determined from samples of saltwater with chloride concentrations of about 1,000 milligrams per liter (mg/L) or greater generally corresponded to a period during which droughts were frequent. Comparison of average daily air temperatures in Miami, Florida, with estimates of recharge temperatures determined from the dissolved gas composition in water samples indicated that saltwater likely entered the aquifer in April or early May when water levels are typically at their lowest during the year. Conversely, most of the samples of freshwater with chloride concentrations less than about 1,000 mg/L indicate recharge temperatures corresponding to air temperatures in mid to late May when rainfall and water levels in the aquifer increase, and the piston-flow ages of these samples correspond to wet years. The piston-flow ages of freshwater samples generally were younger than ages of samples of saltwater. Saltwater samples that were depleted in boron, magnesium, potassium, sodium, and sulfate, and enriched in calcium relative to the concentrations theoretically produced by freshwater/seawater mixing, generally were found to be associated with areas where saltwater had recently intruded. The calcium to (bicarbonate + sulfate) molar ratios (Ca/(HCO 3 +SO 4 )) of these samples generally were greater than 1. Saltwater samples from some of the monitoring wells, however, indicated little or no enrichment or depletion of these ions relative to the theoretical freshwater/seawater mixing line, and the Ca/(HCO 3 +SO 4 ) molar ratios of these samples generally were less than 1. Results indicated that aquifer materials are approaching equilibrium with seawater at these well locations.

Florida↗

USGS2014

Hydrogeologic framework and salinity distribution of the Floridan aquifer system of Broward County, Florida

Concerns about water-level decline and seawater intrusion in the surficial Biscayne aquifer, currently the principal source of water supply to Broward County, prompted a study to refine the hydrogeologic framework of the underlying Floridan aquifer system to evaluate its potential as an alternative source of supply. This report presents cross sections that illustrate the stratigraphy and hydrogeology in eastern Broward County; maps of the upper surfaces and thicknesses of several geologic formations or units within the Floridan aquifer system; and maps of two of the potentially productive water-bearing zones within the system, the Upper Floridan aquifer and the Avon Park permeable zone. An analysis of data on rock depositional textures, associated pore networks, and flow zones in the Floridan aquifer system shows that groundwater moves through the system in two ways. These data support a conceptual, dual-porosity model of the system wherein groundwater moves either as concentrated flow in discrete, thin bedding-plane vugs or zones of vuggy megaporosity, or as diffuse flow through rocks with primarily interparticle and moldic-particle porosity. Because considerable exchange of groundwater may occur between the zones of vuggy and matrix-dominated porosity, understanding the distribution of that porosity and flow zone types is important to evaluating the suitability of the several units within the Floridan aquifer system for managing the water through practices such as aquifer storage and recovery (ASR). The salinity of the water in the Floridan aquifer system is highest in the central part of the study area, and lower toward the north and south. Although salinity generally increases with depth, in the western part of the study area a zone of relatively high saline water is perched above water of lower salinity in the underlying Avon Park permeable zone. Overall, the areas of highest salinity in the aquifer system coincide with those with the lowest estimated transmissivity, so that the occurrence of perched saline water in the system may be the consequence of incompletely flushed connate water or intruded seawater. A seismic reflection profile along the Hillsboro Canal, at the northern edge of the study area, shows seven seismic-sag structures that are interpreted as downward deformation of overlying strata into collapsed deep cave systems. These structures may compromise the integrity of the confinement created by the underlying strata by allowing upconing of saline water from depth, which has implications for successful application of ASR and use of the Floridan aquifer system as an alternative water supply.

Florida↗

USGS2014

Geologic and hydrogeologic frameworks of the Biscayne aquifer in central Miami-Dade County, Florida

Evaluations of the lithostratigraphy, lithofacies, paleontology, ichnology, depositional environments, and cyclostratigraphy from 11 test coreholes were linked to geophysical interpretations, and to results of hydraulic slug tests of six test coreholes at the Snapper Creek Well Field (SCWF), to construct geologic and hydrogeologic frameworks for the study area in central Miami-Dade County, Florida. The resulting geologic and hydrogeologic frameworks are consistent with those recently described for the Biscayne aquifer in the nearby Lake Belt area in Miami-Dade County and link the Lake Belt area frameworks with those developed for the SCWF study area. The hydrogeologic framework is characterized by a triple-porosity pore system of (1) matrix porosity (mainly mesoporous interparticle porosity, moldic porosity, and mesoporous to megaporous separate vugs), which under dynamic conditions, produces limited flow; (2) megaporous, touching-vug porosity that commonly forms stratiform groundwater passageways; and (3) conduit porosity, including bedding-plane vugs, decimeter-scale diameter vertical solution pipes, and meter-scale cavernous vugs. The various pore types and associated permeabilities generally have a predictable vertical spatial distribution related to the cyclostratigraphy. The Biscayne aquifer within the study area can be described as two major flow units separated by a single middle semiconfining unit. The upper Biscayne aquifer flow unit is present mainly within the Miami Limestone at the top of the aquifer and has the greatest hydraulic conductivity values, with a mean of 8,200 feet per day. The middle semiconfining unit, mainly within the upper Fort Thompson Formation, comprises continuous to discontinuous zones with (1) matrix porosity; (2) leaky, low permeability layers that may have up to centimeter-scale vuggy porosity with higher vertical permeability than horizontal permeability; and (3) stratiform flow zones composed of fossil moldic porosity, burrow related vugs, or irregular vugs. Flow zones with a mean hydraulic conductivity of 2,600 feet per day are present within the middle semiconfining unit, but none of the flow zones are continuous across the study area. The lower Biscayne aquifer flow unit comprises a group of flow zones in the lower part of the aquifer. These flow zones are present in the lower part of the Fort Thompson Formation and in some cases within the limestone or sandstone or both in the uppermost part of the Pinecrest Sand Member of the Tamiami Formation. The mean hydraulic conductivity of major flow zones within the lower Biscayne aquifer flow unit is 5,900 feet per day, and the mean value for minor flow zones is 2,900 feet per day. A semiconfining unit is present beneath the Biscayne aquifer. The boundary between the two hydrologic units is at the top or near the top of the Pinecrest Sand Member of the Tamiami Formation. The lower semiconfining unit has a hydraulic conductivity of less than 350 feet per day. The most productive zones of groundwater flow within the two Biscayne aquifer flow units have a characteristic pore system dominated by stratiform megaporosity related to selective dissolution of an Ophiomorpha-dominated ichnofabric. In the upper flow unit, decimeter-scale vertical solution pipes that are common in some areas of the SCWF study area contribute to high vertical permeability compared to that in areas without the pipes. Cross-hole flowmeter data collected from the SCWF test coreholes show that the distribution of vuggy porosity, matrix porosity, and permeability within the Biscayne aquifer of the SCWF is highly heterogeneous and anisotropic. Groundwater withdrawals from production well fields in southeastern Florida may be inducing recharge of the Biscayne aquifer from canals near the well fields that are used for water-management functions, such as flood control and well-field pumping. The SCWF was chosen as a location within Miami-Dade County to study the potential for such recharge to the Biscayne aquifer from the C–2 (Snapper Creek) canal that roughly divides the well field in half. Geologic, hydrogeologic, and hydraulic information on the aquifer collected during construction of monitoring wells within the SCWF could be used to evaluate the groundwater flow budget at the well-field scale.

Florida↗

USGS2013

Evaluation of permeability and non-Darcy flow in vuggy macroporous limestone aquifer samples with lattice Boltzmann methods

Lattice Boltzmann flow simulations provide a physics-based means of estimating intrinsic permeability from pore structure and accounting for inertial flow that leads to departures from Darcy's law. Simulations were used to compute intrinsic permeability where standard measurement methods may fail and to provide better understanding of departures from Darcy's law under field conditions. Simulations also investigated resolution issues. Computed tomography (CT) images were acquired at 0.8 mm interscan spacing for seven samples characterized by centimeter-scale biogenic vuggy macroporosity from the extremely transmissive sole-source carbonate karst Biscayne aquifer in southeastern Florida. Samples were as large as 0.3 m in length; 7–9 cm-scale-length subsamples were used for lattice Boltzmann computations. Macroporosity of the subsamples was as high as 81%. Matrix porosity was ignored in the simulations. Non-Darcy behavior led to a twofold reduction in apparent hydraulic conductivity as an applied hydraulic gradient increased to levels observed at regional scale within the Biscayne aquifer; larger reductions are expected under higher gradients near wells and canals. Thus, inertial flows and departures from Darcy's law may occur under field conditions. Changes in apparent hydraulic conductivity with changes in head gradient computed with the lattice Boltzmann model closely fit the Darcy-Forchheimer equation allowing estimation of the Forchheimer parameter. CT-scan resolution appeared adequate to capture intrinsic permeability; however, departures from Darcy behavior were less detectable as resolution coarsened.

Water Resources Research↗

USGS2013

Preliminary stratigraphic and hydrogeologic cross sections and seismic profile of the Floridan aquifer system of Broward County, Florida

To help water-resource managers evaluate the Floridan aquifer system (FAS) as an alternative water supply, the U.S. Geological Survey initiated a study, in cooperation with the Broward County Environmental Protection and Growth Management Department, to refine the hydrogeologic framework of the FAS in the eastern part of Broward County. This report presents three preliminary cross sections illustrating stratigraphy and hydrogeology in eastern Broward County as well as an interpreted seismic profile along one of the cross sections. Marker horizons were identified using borehole geophysical data and were initially used to perform well-to-well correlation. Core sample data were integrated with the borehole geophysical data to support stratigraphic and hydrogeologic interpretations of marker horizons. Stratigraphic and hydrogeologic units were correlated across the county using borehole geophysical data from multiple wells. Seismic-reflection data were collected along the Hillsboro Canal. Borehole geophysical data were used to identify and correlate hydrogeologic units in the seismic-reflection profile. Faults and collapse structures that intersect hydrogeologic units were also identified in the seismic profile. The information provided in the cross sections and the seismic profile is preliminary and subject to revision.

Florida↗

USGS2013

Integrating seismic-reflection and sequence-stratigraphic methods to characterize the hydrogeology of the Floridan aquifer system in southeast Florida

The Floridan aquifer system (FAS) is receiving increased attention as a result of regulatory restrictions on water-supply withdrawals and treated wastewater management practices. The South Florida Water Management District’s Regional Water Availability Rule, adopted in 2007, restricts urban withdrawals from the shallower Biscayne aquifer to pre-April 2006 levels throughout southeast Florida. Legislation adopted by the State of Florida requires elimination of ocean outfalls of treated wastewater by 2025. These restrictions have necessitated the use of the more deeply buried FAS as an alternate water resource to meet projected water-supply shortfalls, and as a repository for the disposal of wastewater via Class I deep injection wells and injection of reclaimed water. Some resource managers in Broward County have expressed concern regarding the viability of the FAS as an alternative water supply due to a lack of technical data and information regarding its long-term sustainability. Sustainable development and management of the FAS for water supply is uncertain because of the potential risk posed by structural geologic anomalies (faults, fractures, and karst collapse structures) and knowledge gaps in the stratigraphy of the system. The integration of seismic-reflection and borehole data into an improved geologic and hydrogeologic framework will provide a better understanding of the structural and stratigraphic features that influence groundwater flow and contaminant transport.

Florida↗