USGS Science⌕ Search

Geology topics

K.J. Cunningham

Publications and source records attributed to K.J. Cunningham.

12 recordsLinked to original sources

Multiple technologies applied to characterization of the porosity and permeability of the Biscayne aquifer, Florida

Research is needed to determine how seepage-control actions planned by the Comprehensive Everglades Restoration Plan (CERP) will affect recharge, groundwater flow, and discharge within the dual-porosity karstic Biscayne aquifer where it extends eastward from the Everglades to Biscayne Bay. A key issue is whether the plan can be accomplished without causing urban flooding in adjacent populated areas and diminishing coastal freshwater flow needed in the restoration of the ecologic systems. Predictive simulation of groundwater flow is a prudent approach to understanding hydrologic change and potential ecologic impacts. A fundamental problem to simulation of karst groundwater flow is how best to represent aquifer heterogeneity. Currently, U.S. Geological Survey (USGS) researchers and academic partners are applying multiple innovative technologies to characterize the spatial distribution of porosity and permeability within the Biscayne aquifer.

Open-File Report↗

NMR imaging of fluid exchange between macropores and matrix in eogenetic karst

Sequential time-step images acquired using nuclear magnetic resonance (NMR) show the displacement of deuterated water (D 2 O) by fresh water within two limestone samples characterized by a porous and permeable limestone matrix of peloids and ooids. These samples were selected because they have a macropore system representative of some parts of the eogenetic karst limestone of the Biscayne Aquifer in southeastern Florida. The macroporosity, created by the trace fossil Ophiomorpha, is principally well connected and of centimeter scale. These macropores occur in broadly continuous stratiform zones that create preferential flow layers within the hydrogeologic units of the Biscayne. This arrangement of porosity is important because in coastal areas, it could produce a preferential pathway for salt water intrusion. Two experiments were conducted in which samples saturated with D 2 O were placed in acrylic chambers filled with fresh water and examined with NMR. Results reveal a substantial flux of fresh water into the matrix porosity with a simultaneous loss of D 2 O. Specifically, we measured rates upward of 0.001 mL/h/g of sample in static conditions, and perhaps as great as 0.07 mL/h/g of sample when fresh water continuously flows past a sample at velocities less than those found within stressed areas of the Biscayne. These experiments illustrate how fresh water and D 2 O, with different chemical properties, migrate within one type of matrix porosity found in the Biscayne. Furthermore, these experiments are a comparative exercise in the displacement of sea water by fresh water in the matrix of a coastal, karst aquifer since D 2 O has a greater density than fresh water. ?? 2008 National Ground Water Association.

Ground Water↗

Prominence of ichnologically influenced macroporosity in the karst Biscayne aquifer: Stratiform "super-K" zones

A combination of cyclostratigraphic, ichnologic, and borehole geophysical analyses of continuous core holes; tracer-test analyses; and lattice Boltzmann flow simulations was used to quantify biogenic macroporosity and permeability of the Biscayne aquifer, southeastern Florida. Biogenic macroporosity largely manifests as: (1) ichnogenic macroporosity primarily related to postdepositional burrowing activity by callianassid shrimp and fossilization of components of their complex burrow systems (Ophiomorpha); and (2) biomoldic macroporosity originating from dissolution of fossil hard parts, principally mollusk shells. Ophiomorpha-dominated ichno-fabric provides the greatest contribution to hydrologic characteristics in the Biscayne aquifer in a 345 km 2 study area. Stratiform tabular-shaped units of thalassinidean-associated macroporosity are commonly confined to the lower part of upward-shallowing high-frequency cycles, throughout aggradational cycles, and, in one case, they stack vertically within the lower part of a high-frequency cycle set. Broad continuity of many of the macroporous units concentrates groundwater flow in extremely permeable passage-ways, thus making the aquifer vulnerable to long-distance transport of contaminants. Ichnogenic macroporosity represents an alternative pathway for concentrated groundwater flow that differs considerably from standard karst flow-system paradigms, which describe groundwater movement through fractures and cavernous dissolution features. Permeabilities were calculated using lattice Boltzmann methods (LBMs) applied to computer renderings assembled from X-ray computed tomography scans of various biogenic macroporous limestone samples. The highest simulated LBM permeabilities were about five orders of magnitude greater than standard laboratory measurements using air-permeability methods, which are limited in their application to extremely permeable macroporous rock samples. Based on their close conformance to analytical solutions for pipe flow, LBMs offer a new means of obtaining accurate permeability values for such materials. We suggest that the stratiform ichnogenic groundwater flow zones have permeabilities even more extreme (???2-5 orders of magnitude higher) than the Jurassic "super-K" zones of the giant Ghawar oil field. The flow zones of the Pleistocene Biscayne aquifer provide examples of ichnogenic macroporosity for comparative analysis of origin and evolution in other carbonate aquifers, as well as petroleum reservoirs. ?? 2008 Geological Society of America.

Geological Society of America Bulletin↗

First documentation of tidal-channel sponge biostromes (upper Pleistocene, southeastern Florida)

Sponges are not a common principal component of Cenozoic reefs and are more typically dominant in deep-water and/or cold-water localities. Here we report the discovery of extensive upper Pleistocene shallow-marine, tropical sponge biostromes from the Mami Limestone of southeastern Florida built by a new ceractinomorph demosponge. These upright, barrel- to vase-shaped sponges occur in monospecific aggregations constructed within the tidal channels of an oolitic tidal-bar belt similar to modern examples on the Great Bahama Bank. The biostromes appear to have a ribbon-like geometry, with densely spaced sponges populating a paleochannel along a 3.5 km extent in the most lengthy biostrome. These are very large (as high as 2 m and 1.8 m in diameter), particularly well-preserved calcified sponges with walls as hard as concrete. Quartz grains are the most common particles agglutinated in the structure of the sponge walls. Where exposed, sediment fill between the sponges is commonly a highly burrowed or cross-bedded ooid-bearing grainstone and, locally, quartz sand. It is postulated that the dense, localized distribution of these particular sponges was due to a slight edge over competitors for food or energy supply and space in a stressed environment of tidal-influenced salinity and nutrient changes, strong currents, and frequently shifting submarine sand dunes. To our knowledge, this represents the first documentation of sponge biostromes composed of very large upright sponges within high-energy tidal channels between ooid shoals. The remarkably well-preserved accumulations provide an alternative example of sponge reefs for comparative paleoenvironmental studies. ?? 2007 The Geological Society of America.

Geology↗

Assessing the vulnerability of a municipal well field to contamination in a karst aquifer

Proposed expansion of extractive lime-rock mines near the Miami-Dade County Northwest well field and Everglades wetland areas has garnered intense scrutiny by government, public, environmental stakeholders, and the media because of concern that mining will increase the risk of pathogen contamination. Rock mines are excavated to the same depth as the well field's primary producing zone. The underlying karst Biscayne aquifer is a triple-porosity system characterized by (1) a matrix of interparticle porosity and separate vug porosity; (2) touching-vug porosity that forms preferred, stratiform passageways; and, less commonly, (3) conduit porosity formed by thin solution pipes, bedding-plane vugs, and cavernous vugs. Existing ground-water flow and particle tracking models do not provide adequate information regarding the ability the aquifer to limit the advective movement of pathogens and other contaminants. Chemical transport and colloidal mobility properties have been delineated using conservative and microsphere-surrogate tracers for Cryptosporidium parvum. Forced-gradient tests were executed by introducing conservative tracers into injection wells located 100 m (328 ft) from a municipal-supply well. Apparent mean advective velocity between the wells is one to two orders of magnitude greater than previously measured. Touching-vug, stratiform flow zones are efficient pathways for tracer movement at the well field. The effective porosity for a continuum model between the point of injection and tracer recovery ranges from 2 to 4 percent and is an order of magnitude smaller than previously assumed. Existing well-field protection zones were established using porosity estimates based on specific yield. The effective, or kinematic, porosity of a Biscayne aquifer continuum model is lower than the total porosity, because high velocities occur along preferential flow paths that result in faster times of travel than can be represented with the ground-water flow equation. Tracer tests indicate that the relative ease of contaminant movement to municipal supply wells is much greater than previously considered.

Environmental & Engineering Geoscience↗

Application of ground-penetrating radar, digital optical borehole images, and cores for characterization of porosity hydraulic conductivity and paleokarst in the Biscayne aquifer, southeastern Florida, USA

This paper presents examples of ground-penetrating radar (GPR) data from two study sites in southeastern Florida where karstic Pleistocene platform carbonates that comprise the unconfined Biscayne aquifer were imaged. Important features shown on resultant GPR profiles include: (1) upward and lateral qualitative interpretative distribution of porosity and hydraulic conductivity; (2) paleotopographic relief on karstic subaerial exposure surfaces; and (3) vertical stacking of chronostratigraphic high-frequency cycles (HFCs). These characteristics were verified by comparison to rock properties observed and measured in core samples, and identified in digital optical borehole images. Results demonstrate that an empirical relation exists between measured whole-core porosity and hydraulic conductivity, observed porosity on digital optical borehole images, formation conductivity, and GPR reflection amplitudes-as porosity and hydraulic conductivity determined from core and borehole images increases, formation conductivity increases, and GPR reflection amplitude decreases. This relation allows for qualitative interpretation of the vertical and lateral distribution of porosity and hydraulic conductivity within HFCs. Two subtidal HFCs in the uppermost Biscayne aquifer have significantly unique populations of whole-core porosity values and vertical hydraulic conductivity values. Porosity measurements from one cycle has a median value about two to three times greater than the values from the other HFC, and median values of vertical hydraulic-conductivity about three orders of magnitude higher than the other HFC. The HFC with the higher porosity and hydraulic conductivity values is shown as a discrete package of relatively low-amplitude reflections, whereas the HFC characterized by lower porosity and hydraulic-conductivity measurements is expressed by higher amplitude reflections. Porosity and hydraulic-conductivity values measured from whole-core samples, and vuggy porosity identified on digital borehole images from shallowing-upward, peritidal HFCs show that the highest porosity occurs at the base of the cycles, moderate porosity at the middle of the cycles, and lowest porosity occurs at the top of cycles. Hydraulic conductivity is also highest at the base of the peritidal cycles and lowest in the middle to upper parts of cycles. This change in porosity and hydraulic conductivity from bottom to top is visible as an upward variation in reflection amplitude on GPR profiles-lowest amplitudes at the base and highest at the cycle tops. This study demonstrates that GPR can be used to show the qualitative distribution of porosity and hydraulic conductivity within a cycle-stratigraphic framework composed of carbonate HFCs. The distribution of porosity and hydraulic conductivity within HFCs is related to depositional textures. The upward and lateral patterns of the rock facies within the HFCs can be translated to geophysical-log properties and radar facies configurations that could aid in interpretation and prediction of ground-water flow through a carbonate aquifer. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Applied Geophysics↗

New method for quantification of vuggy porosity from digital optical borehole images as applied to the karstic Pleistocene limestone of the Biscayne aquifer, southeastern Florida

Vuggy porosity is gas- or fluid-filled openings in rock matrix that are large enough to be seen with the unaided eye. Well-connected vugs can form major conduits for flow of ground water, especially in carbonate rocks. This paper presents a new method for quantification of vuggy porosity calculated from digital borehole images collected from 47 test coreholes that penetrate the karstic Pleistocene limestone of the Biscayne aquifer, southeastern Florida. Basically, the method interprets vugs and background based on the grayscale color of each in digital borehole images and calculates a percentage of vuggy porosity. Development of the method was complicated because environmental conditions created an uneven grayscale contrast in the borehole images that makes it difficult to distinguish vugs from background. The irregular contrast was produced by unbalanced illumination of the borehole wall, which was a result of eccentering of the borehole-image logging tool. Experimentation showed that a simple, single grayscale threshold would not realistically differentiate between the grayscale contrast of vugs and background. Therefore, an equation was developed for an effective subtraction of the changing grayscale contrast, due to uneven illumination, to produce a grayscale threshold that successfully identifies vugs. In the equation, a moving average calculated around the circumference of the borehole and expressed as the background grayscale intensity is defined as a baseline from which to identify a grayscale threshold for vugs. A constant was derived empirically by calibration with vuggy porosity values derived from digital images of slabbed-core samples and used to make the subtraction from the background baseline to derive the vug grayscale threshold as a function of azimuth. The method should be effective in estimating vuggy porosity in any carbonate aquifer. ?? 2003 Published by Elsevier B.V.

Journal of Applied Geophysics↗

Interplay of late Cenozoic siliciclastic supply and carbonate response on the southeast Florida platform

High-resolution seismic-reflection data collected along the length of the Caloosahatchee River in southwestern Florida have been correlated to nannofossil biostratigraphy and strontium-isotope chemostratigraphy at six continuously cored boreholes. These data are interpreted to show a major Late Miocene(?) to Early Pliocene fluvial-deltaic depositional system that prograded southward across the carbonate Florida Platform , interrupting nearly continuous carbonate deposition since early in the Cretaceous. Connection of the platform top to a continental source of siliciclastics and significant paleotopography combined to focus accumulation of an immense supply of siliciclastics on the southeastern part of the Florida Platform . The remarkably thick (> 100 m), sand-rich depositional system, which is characterized by clinoformal progradation, filled in deep accommodation, while antecedent paleotopography directed deltaic progradation southward within the middle of the present-day Florida Peninsula. The deltaic depositional system may have prograded about 200 km southward to the middle and upper Florida Keys, where Late Miocene to Pliocene siliciclastic form the foundation of the Quaternary carbonate shelf and shelf margin of the Florida Keys. These far-traveled siliciclastic deposits filled accommodation on the southeastern part of the Florida Platform so that paleobathymetry was sufficiently shallow to allow Quaternary recovery of carbonate sedimentation in the area of southern peninsular Florida and the Florida Keys.

Florida↗

Copepod communities from surface and ground waters in the everglades, south Florida

We studied species composition and individual abundance of copepods in the surficial aquifer northeast of Everglades National Park. We identified the spatial distribution of subsurface habitats by assessing the depth of the high porosity layers in the limestone along a canal system, and we used copepods to assess the exchange between surface water and ground water along canal banks, at levels in the wells where high porosity connections to the canals exist. Surface- and ground-water taxa were defined, and species composition was related to areal position, sampling depth, and time. Subsurface copepod communities were dominated by surface copepods that disperse into the aquifer following the groundwater seepage along canal L-31N. The similarities in species composition between wells along canal reaches, suggest that copepods mainly enter ground water horizontally along canals via active and passive dispersal. Thus, the copepod populations indicate continuous connections between surface- and ground waters. The most abundant species were Orthocyclops modestus, Arctodiaptomus floridanus, Mesocyclops edax, and Thermocyclops parvus, all known in literature from surface habitats; however, these species have been collected in ground water in ENP. Only two stygophiles were collected: Diacylcops nearcticus and Diacyclops crassicaudis brachycercus. Restoration of the Everglades ecosystem requires a mosaic of data to reveal a complete picture of this complex system. The use of copepods as indicators of seepage could be a tool in helping to assess the direction and the duration of surface and ground water exchange.

Florida↗

Sequence-stratigraphic analysis of the Regional Observation Monitoring Program (ROMP) 29A test corehole and its relation to carbonate porosity and regional transmissivity in the Floridan aquifer system, Highlands County, Florida

An analysis was made to describe and interpret the lithology of a part of the Upper Floridan aquifer penetrated by the Regional Observation Monitoring Program (ROMP) 29A test corehole in Highlands County, Florida. This information was integrated into a one-dimensional hydrostratigraphic model that delineates candidate flow zones and confining units in the context of sequence stratigraphy. Results from this test corehole will serve as a starting point to build a robust three-dimensional sequence-stratigraphic framework of the Floridan aquifer system. The ROMP 29A test corehole penetrated the Avon Park Formation, Ocala Limestone, Suwannee Limestone, and Hawthorn Group of middle Eocene to Pliocene age. The part of the Avon Park Formation penetrated in the ROMP 29A test corehole contains two composite depositional sequences. A transgressive systems tract and a highstand systems tract were interpreted for the upper composite sequence; however, only a highstand systems tract was interpreted for the lower composite sequence of the deeper Avon Park stratigraphic section. The composite depositional sequences are composed of at least five high-frequency depositional sequences. These sequences contain high-frequency cycle sets that are an amalgamation of vertically stacked high-frequency cycles. Three types of high-frequency cycles have been identified in the Avon Park Formation: peritidal, shallow subtidal, and deeper subtidal high-frequency cycles. The vertical distribution of carbonate-rock diffuse flow zones within the Avon Park Formation is heterogeneous. Porous vuggy intervals are less than 10 feet, and most are much thinner. The volumetric arrangement of the diffuse flow zones shows that most occur in the highstand systems tract of the lower composite sequence of the Avon Park Formation as compared to the upper composite sequence, which contains both a backstepping transgressive systems tract and a prograding highstand systems tract. Although the porous and permeable layers are not thick, some intervals may exhibit lateral continuity because of their deposition on a broad low-relief ramp. A thick interval of thin vuggy zones and open faults forms thin conduit flow zones mixed with relatively thicker carbonate-rock diffuse flow zones between a depth of 1,070 and 1,244 feet below land surface (bottom of the test corehole). This interval is the most transmissive part of the Avon Park Formation penetrated in the ROMP 29A test corehole and is included in the highstand systems tract of the lower composite sequence. The Ocala Limestone is considered to be a semiconfining unit and contains three depositional sequences penetrated by the ROMP 29A test corehole. Deposited within deeper subtidal depositional cycles, no zones of enhanced porosity and permeability are expected in the Ocala Limestone. A thin erosional remnant of the shallow marine Suwannee Limestone overlies the Ocala Limestone, and permeability seems to be comparatively low because moldic porosity is poorly connected. Rocks that comprise the lower Hawthorn Group, Suwannee Limestone, and Ocala Limestone form a permeable upper zone of the Upper Floridan aquifer, and rocks of the lower Ocala Limestone and Avon Park Formation form a permeable lower zone of the Upper Floridan aquifer. On the basis of a preliminary analysis of transmissivity estimates for wells located north of Lake Okeechobee, spatial relations among groups of relatively high and low transmissivity values within the upper zone are evident. Upper zone transmissivity is generally less than 10,000 feet squared per day in areas located south of a line that extends through Charlotte, Sarasota, DeSoto, Highlands, Polk, Osceola, Okeechobee, and St. Lucie Counties. Transmissivity patterns within the lower zone of the Avon Park Formation cannot be regionally assessed because insufficient data over a wide areal extent have not been compiled.

Florida↗

New Tertiary stratigraphy for the Florida Keys and southern peninsula of Florida

Seven lithologic formations, ranging in age from Oligocene to Pleistocene, were recently penetrated by core holes in southernmost Florida. From bottom to top, they are the early Oligocene Suwannee Limestone; late-early Oligocene-to-Miocene Arcadia Formation, basal Hawthorn Group; late Miocene Peace River Formation, upper Hawthorn Group; newly proposed late Miocene-to-Pliocene Long Key and Stock Island Formations; and Pleistocene Key Largo and Miami Limestones. The rocks of the Suwannee Limestone form a third-order sequence. Although the entire thickness was not penetrated, 96 m of Suwannee core from one well contains at least 50 vertically stacked, exposure-capped limestone cycles, presumably related to rapid eustatic fluctuations while experiencing tropical to subtropical conditions. The Arcadia Formation is a composite sequence containing four high-frequency sequences composed of multiple vertically stacked carbonate cycles. Most cycles do not show evidence of subaerial exposure and were deposited under more temperate conditions, relative to the Suwannee Limestone. The Arcadia Formation in southernmost Florida is bounded by regional unconformities representing third-order sequence boundaries. Post-Arcadia transgression produced a major backstepping of sediment accumulation above the upper sequence boundary of the Arcadia Formation. The Peace River Formation, composed of diatomaceous mudstones, has been identified only beneath the Florida peninsula and is not present beneath the Florida Keys. Deposition occurred during marine transgressive to high-stand conditions and a local phosphatization event (recorded in northeast Florida). The transgression is possibly related to a global rise in sea level, which resulted in upwelling of relatively cooler, relatively nutrient-rich water masses onto the Florida Platform. It is proposed that the absence of Peace River sediments beneath the Keys is due to sediment bypass of the upper surface of the Arcadia, a result of sediment sweeping by an ancestral Florida current. During late Miocene to Pliocene time in the Florida Keys, siliciclastics of the Long Key Formation and fine-grained carbonates of the Stock Island Formation prograded toward the southern edge of the Florida Platform and downlapped onto the regional unconformity at the top of the Arcadia. Shallow-marine Pleistocene limestones (Key Largo and Miami Limestones), deposited during tropical to subtropical conditions, drape over accretionary successions of the Long Key and Stock Island Formations.

Florida↗