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.

Trans Terrains: Gendered Embodiments and Religious Landscapes in Yogyakarta, Indonesia

Transgendered Indonesians live in the fourth most populated nation in the world with more Muslims than any other country. This thesis summarizes an ethnography conducted on one religiously oriented male-to-female transgender community known in the city of Yogyakarta as the waria. This study analyzes the waria’s gender and religious identities from an emic and etic perspective, focusing on how individuals comport themselves inside the world’s first transgender mosque-like institution called a pesantren waria. The waria take their name from the Indonesian words wanita (woman) and pria (man). I will chart how this male-to-female population create spaces of spiritual belonging and physical security within a territory that has experienced geo-religio-political insecurity: natural disasters, fundamentalist movements, and toppling dictatorships. This work illuminates how the waria see themselves as biologically male, not men. Anatomy is not what gives the waria their gender, their feminine expression and sexual attraction does. Although the waria self-identity as women/waria, in a religious context they perform as men, not women.

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.

A hermeneutic phenomenology of graduate education students' understandings of instructor power in a higher education classroom

Higher education is a distribution center of knowledge and economic, social, and cultural power (Cervero & Wilson, 2001). A critical approach to understanding a higher education classroom begins with recognizing the instructor's position of power and authority (Tisdell, Hanley, & Taylor, 2000). The power instructors wield exists mostly unquestioned, allowing for teaching practices that reproduce the existing societal patterns of inequity in the classroom (Brookfield, 2000). ^ The purpose of this hermeneutic phenomenological study was to explore students' experiences with the power of their instructors in a higher education classroom. A hermeneutic phenomenological study intertwines the interpretations of both the participants and the researcher about a lived experience to uncover layers of meaning because the meanings of lived experiences are usually not readily apparent (van Manen, 1990). Fifteen participants were selected using criterion, convenience, and snowball sampling. The primary data gathering method were semi-structured interviews guided by an interview protocol (Creswell, 2003). Data were interpreted using thematic reflection (van Manen, 1990). ^ Three themes emerged from data interpretation: (a) structuring of instructor-student relationships, (b) connecting power to instructor personality, and (c) learning to navigate the terrains of higher education. How interpersonal relationships were structured in a higher education classroom shaped how students perceived power in that higher education classroom. Positive relationships were described using the metaphor of family and a perceived ethic of caring and nurturing by the instructor. As participants were consistently exposed to exercises of instructor power in a higher education classroom, they attributed those exercises of power to particular instructor traits rather than systemic exercises of power. As participants progressed from undergraduate to graduate studies, they perceived the benefits of expertise in content or knowledge development as secondary to expertise in successfully navigating the social, cultural, political, and interpersonal terrains of higher education. Ultimately, participants expressed that higher education is not about what you know; it is about learning how to play the game. Implications for teaching in higher education and considerations for future research conclude the study.^

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.

A generalized adaptive mathematical morphological filter for LIDAR data

Airborne Light Detection and Ranging (LIDAR) technology has become the primary method to derive high-resolution Digital Terrain Models (DTMs), which are essential for studying Earth's surface processes, such as flooding and landslides. The critical step in generating a DTM is to separate ground and non-ground measurements in a voluminous point LIDAR dataset, using a filter, because the DTM is created by interpolating ground points. As one of widely used filtering methods, the progressive morphological (PM) filter has the advantages of classifying the LIDAR data at the point level, a linear computational complexity, and preserving the geometric shapes of terrain features. The filter works well in an urban setting with a gentle slope and a mixture of vegetation and buildings. However, the PM filter often removes ground measurements incorrectly at the topographic high area, along with large sizes of non-ground objects, because it uses a constant threshold slope, resulting in "cut-off" errors. A novel cluster analysis method was developed in this study and incorporated into the PM filter to prevent the removal of the ground measurements at topographic highs. Furthermore, to obtain the optimal filtering results for an area with undulating terrain, a trend analysis method was developed to adaptively estimate the slope-related thresholds of the PM filter based on changes of topographic slopes and the characteristics of non-terrain objects. The comparison of the PM and generalized adaptive PM (GAPM) filters for selected study areas indicates that the GAPM filter preserves the most "cut-off" points removed incorrectly by the PM filter. The application of the GAPM filter to seven ISPRS benchmark datasets shows that the GAPM filter reduces the filtering error by 20% on average, compared with the method used by the popular commercial software TerraScan. The combination of the cluster method, adaptive trend analysis, and the PM filter allows users without much experience in processing LIDAR data to effectively and efficiently identify ground measurements for the complex terrains in a large LIDAR data set. The GAPM filter is highly automatic and requires little human input. Therefore, it can significantly reduce the effort of manually processing voluminous LIDAR measurements.

Loose Ends

Loose Ends is a collection of lyric and narrative poems that explores the multiple terrains of identity—individual, cultural, and historical. The poems embrace the essential incoherence of the self, resisting monolithic identity in favor of a multi-faceted, historically complex, imagistic rendering of the inner life. At its heart, the collection seeks to grapple with the gravitas of living: the continual assault of history and nature on human agency, the staggering context of the universe as a backdrop for communal and individual struggle. While single poems may only touch briefly or incompletely on these themes, the collection as a whole presents an admittedly inchoate picture of contemporary American identity.