Lattice Boltzmann modeling of fluid flow and solute transport in karst aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Nonlinear Flow in Karst Formations

The variation of effective hydraulic conductivity as a function of specific discharge in several 0.2-m and 0.3-m cubes of Key Largo Limestone was investigated. The experimental results closely match the Forchheimer equation. Defining the pore-size length scale in terms of Forchheimer parameters, it is demonstrated that significant deviations from Darcian flow will occur when the Reynolds number exceeds 0.11. A particular threshold model previously proposed for use in karstic formations does not show strong agreement with the data near the onset of nonlinear flow.

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

Comparison of snail density, standing stock, and body size between Caribbean karst wetlands and other freshwater ecosystems

Synthesizing data from multiple studies generates hypotheses about factors that affect the distribution and abundance of species among ecosystems. Snails are dominant herbivores in many freshwater ecosystems, but there is no comprehensive review of snail density, standing stock, or body size among freshwater ecosystems. We compile data on snail density and standing stock, estimate body size with their quotient, and discuss the major pattern that emerges. We report data from 215 freshwater ecosystems taken from 88 studies that we placed into nine categories. Sixty-five studies reported density, seven reported standing stock, and 16 reported both. Despite the breadth of studies, spatial and temporal sampling scales were limited. Researchers used 25 different sampling devices ranging in area from 0.0015 to 2.5 m2. Most ecosystem categories had similar snail densities, standing stocks, and body sizes suggesting snails shared a similar function among ecosystems. Caribbean karst wetlands were a striking exception with much lower density and standing stock, but large body size. Disparity in body size results from the presence of ampullariids in Caribbean karst wetlands suggesting that biogeography affects the distribution of taxa, and in this case size, among aquatic ecosystems. We propose that resource quality explains the disparity in density and standing stock between Caribbean karst wetlands and other categories. Periphyton in Caribbean karst wetlands has high carbon-to-phosphorous ratios and defensive characteristics that inhibit grazers. Unlike many freshwater ecosystems where snails are key grazers, we hypothesize that a microbial loop captures much of the primary production in Caribbean karst wetlands.

Fish Community Responses to the Combined Effects of Decreased Hydroperiod and Nonnative Fish Invasions in a Karst Wetland: Are Everglades Solution Holes Sinks for Native Fishes?

Short-hydroperiod Everglades wetlands have been disproportionately affected by reductions in freshwater inflows, land conversion and biotic invasions. Severe hydroperiod reductions in these habitats, including the Rocky Glades, coupled with proximity to canals that act as sources of invasions, may limit their ability to support high levels of aquatic production. We examined whether karst solution holes function as dry-down refuges for fishes, providing a source of marsh colonists upon reflooding, by tracking fish abundance, nonnative composition, and survival in solution holes throughout the dry season. We paired field surveys with an in situ nonnative predation experiment that tested the effects of predation by the recent invader, African jewelfish (Hemichromis letourneuxi) on native fishes. Over the 3 years surveyed, a large number of the solution holes dried before the onset of the wet season, while those retaining water had low survivorship and were dominated by nonnatives. In the experiment, mortality of eastern mosquitofish (Gambusia holbrooki) in the presence of African jewelfish was greater than that associated with deteriorating water quality. Under current water management, findings suggest that solution holes are largely sinks for native fishes, given the high frequency of drydown, extensive period of fish residence, and predation by nonnative fishes.

The Influence of Hydrologic Restoration on Groundwater-Surface Water Interactions in a Karst Wetland, the Everglades (FL, USA)

Efforts to rehydrate and restore surface water flow in karst wetlands can have unintended consequences, as these highly conductive and heterogeneous aquifers create a close connection between groundwater and surface water. Recently, hydrologic restoration efforts in the karstic Taylor Slough portion of the Everglades has changed from point source delivery of canal water (direct restoration), to the use of a series of surface water recharge retention basins (diffuse restoration). To determine the influence of restoration on groundwater-surface water interactions in the Taylor Slough headwaters, a water budget was constructed for 1997–2011 using 70 hydro-meteorological stations. With diffuse restoration, groundwater seepage from the Everglades toward the urban boundary increased, while the downstream delivery of surface water to the main portion of the slough declined. The combined influence of diffuse restoration and climate led to increased intra-annual variability in the volume of groundwater and surface water in storage but supported a more seasonally hydrated wetland compared to the earlier direct tactics. The data further indicated that hydrologic engineering in karst wetland landscapes enhances groundwater-surface water interactions, even those designed for restoration purposes.

Lattice boltzmann modeling of fluid flow to determine the permeability of a karst specimen

A limestone sample was scanned using computed tomography (CT) and the hydraulic conductivity of the 3D reconstructed sample was determined using Lattice- Boltzmann methods (LBM) at varying scales. Due to the shape and size of the original sample, it was challenging to obtain a consistent rectilinear test sample. Through visual inspection however, 91 mm and 76 mm samples were digitally cut from the original. The samples had porosities of 58% and 64% and produced hydraulic conductivity values of K= 13.5 m/s and K=34.5 m/s, respectively. Both of these samples were re-sampled to 1/8 and 1/64 of their original size to produce new virtual samples at lower resolutions of 0.542 mm/lu and 1.084 mm/lu, while still representing the same physical dimensions. The hydraulic conductivity tended to increase slightly as the resolution became coarser. In order to determine an REV, the 91 mm sample was also sub-sampled into blocks that were 1/8 and 1/64 the size of the original. The results were consistent with analytical expectations such as those produced by the Kozeny-Carman equation. A definitive REV size was not reached, however, indicating the need for a larger sample. The methods described here demonstrate the ability of LBM to test rock structures and sizes not normally attainable.

Geophysical Analysis of a Central Florida Karst Terrain using Light Detection and Ranging (LIDAR) and Ground Penetrating Radar (GPR) Derived Surfaces

Airborne LIDAR (Light Detecting and Ranging) is a relatively new technique that rapidly and accurately measures micro-topographic features. This study compares topography derived from LIDAR with subsurface karst structures mapped in 3-dimensions with ground penetrating radar (GPR). Over 500 km of LIDAR data were collected in 1995 by the NASA ATM instrument. The LIDAR data was processed and analyzed to identify closed depressions. A GPR survey was then conducted at a 200 by 600 m site to determine if the target features are associated with buried karst structures. The GPR survey resolved two major depressions in the top of a clay rich layer at ~10m depth. These features are interpreted as buried dolines and are associated spatially with subtle (< 1m) trough-like depressions in the topography resolved from the LIDAR data. This suggests that airborne LIDAR may be a useful tool for indirectly detecting subsurface features associated with sinkhole hazard.

Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Lattice Boltzmann Modeling and Specialized Laboratory Techniques to Determine the Permeability of Megaporous Karst Rock

The Pleistocene carbonate rock Biscayne Aquifer of south Florida contains laterally-extensive bioturbated ooltic zones characterized by interconnected touching-vug megapores that channelize most flow and make the aquifer extremely permeable. Standard petrophysical laboratory techniques may not be capable of accurately measuring such high permeabilities. Instead, innovative procedures that can measure high permeabilities were applied. These fragile rocks cannot easily be cored or cut to shapes convenient for conducting permeability measurements. For the laboratory measurement, a 3D epoxy-resin printed rock core was produced from computed tomography data obtained from an outcrop sample. Permeability measurements were conducted using a viscous fluid to permit easily observable head gradients (~2 cm over 1 m) simultaneously with low Reynolds number flow. For a second permeability measurement, Lattice Boltzmann Method flow simulations were computed on the 3D core renderings. Agreement between the two estimates indicates an accurate permeability was obtained that can be applied to future studies.

Hydrologic Dynamics of the Ground-Water-Dependent Sian Ka’an Wetlands, Mexico, Derived from InSAR and SAR Data

The 5,280 km2 Sian Ka’an Biosphere Reserve includes pristine wetlands fed by ground water from the karst aquifer of the Yucatan Peninsula, Mexico. The inflow through underground karst structures is hard to observe making it difficult to understand, quantify, and predict the wetland dynamics. Remotely sensed Synthetic Aperture Radar (SAR) amplitude and phase observations offer new opportunities to obtain information on hydrologic dynamics useful for wetland management. Backscatter amplitude of SAR data can be used to map flooding extent. Interferometric processing of the backscattered SAR phase data (InSAR) produces temporal phase-changes that can be related to relative water level changes in vegetated wetlands. We used 56 RADARSAT-1 SAR acquisitions to calculate 38 interferograms and 13 flooding maps with 24 day and 48 day time intervals covering July 2006 to March 2008. Flooding extent varied between 1,067 km2 and 2,588 km2 during the study period, and main water input was seen to take place in sloughs during October–December. We propose that main water input areas are associated with water-filled faults that transport ground water from the catchment to the wetlands. InSAR and Landsat data revealed local-scale water divides and surface water flow directions within the wetlands.

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.