2008 Joint Meeting of The Geological Society of America, Soil Science Society of America, American Society of Agronomy, Crop Science Society of America, Gulf Coast Association of Geological Societies with the Gulf Coast Section of SEPM

Paper No. 6
Presentation Time: 9:25 AM

Nuclear Magnetic Resonance Imaging of Groundwater Flow within Centimeter-Scale Macroporosity of the Biscayne Aquifer of Southeast Florida


FLOREA, Lee J., Department of Geography and Geology, Western Kentucky University, 1906 College Heights Blvd, Bowling Green, KY 42101-1066, CUNNINGHAM, Kevin J., U.S. Geological Survey, 7500 SW 36th Street, Fort Lauderdale, FL 33314 and ALTOBELLI, Stephen A., New Mexico Resonance, 2301 Yale Blvd. SE, Suite C-1, Albuquerque, NM 87106, Lee.Florea@wku.edu

In this study, we use an innovative, non-invasive technology, nuclear magnetic resonance imaging (NMRI), to visualize the direction and magnitude of ground water flow in field samples of late Pleistocene limestone of the Biscayne aquifer. One set of NMRI experiments investigates the migration of freshwater into the matrix of permeable (gas minipermeameter mean 10-13.5 m2 on four samples) and porous (mean of 44%, on four samples) peloid-oöid grainstone initially saturated with heavy-water (D2O). In a second set of experiments, the focus of this abstract, we use phase-encoded, stimulated-echo imaging to map the advective velocity of water flowing at two rates of specific discharge (0.025 and 0.013 cm/s) through a 10-cm diameter cylindrical, epoxy-resin model. The model interior accurately reproduces a well-connected maze of ichnologically influenced, centimeter-scale, touching-vug macroporosity common within preferred flow zones in parts of the Biscayne aquifer.

Results of the velocity experiments at both values of specific discharge reveal select macropores with concentrated flow. The average axial components of velocity (0.016 cm/s and 0.04 cm/s) are in the range of advective flow expected in the Biscayne aquifer. Maximum Reynolds number (Re) values in the low and high specific discharge experiments are 15 and 30, respectively. Average values of Re are much lower, ranging between 1 and 6. These low Re values indicate laminar flow conditions in the model. Two further observations point to laminar flow conditions: 1) virtually all values of the axial component of velocity are positive, there are no eddies; and 2) the axial component of the velocity field is proportional between the two values of specific discharge. Exponential functions fit the distributions of the axial velocity. Distributions of the transverse component of velocity are symmetric and appear to fit a Cauchy rather than a Gaussian function.