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 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.

Investigation of spatial patterns of ground-water exchange with lakes, using a three-dimensional numerical model

Hydrogeologic variables controlling groundwater exchange with inflow and flow-through lakes were simulated using a three-dimensional numerical model (MODFLOW) to investigate and quantify spatial patterns of lake bed seepage and hydraulic head distributions in the porous medium surrounding the lakes. Also, the total annual inflow and outflow were calculated as a percentage of lake volume for flow-through lake simulations. The general exponential decline of seepage rates with distance offshore was best demonstrated at lower anisotropy ratio (i.e., Kh/Kv = 1, 10), with increasing deviation from the exponential pattern as anisotropy was increased to 100 and 1000. 2-D vertical section models constructed for comparison with 3-D models showed that groundwater heads and seepages were higher in 3-D simulations. Addition of low conductivity lake sediments decreased seepage rates nearshore and increased seepage rates offshore in inflow lakes, and increased the area of groundwater inseepage on the beds of flow-through lakes. Introduction of heterogeneity into the medium decreased the water table and seepage ratesnearshore, and increased seepage rates offshore in inflow lakes. A laterally restricted aquifer located at the downgradient side of the flow-through lake increased the area of outseepage. Recharge rate, lake depth and lake bed slope had relatively little effect on the spatial patterns of seepage rates and groundwater exchange with lakes.

Short-term changes in phosphorus storage in an oligotrophic Everglades wetland ecosystem receiving experimental nutrient enrichment

Natural, unenriched Everglades wetlands are known to be limited by phosphorus (P) and responsive to P enrichment. However, whole-ecosystem evaluations of experimental P additions are rare in Everglades or other wetlands. We tested the response of the Everglades wetland ecosystem to continuous, low-level additions of P (0, 5, 15, and 30 μg L−1 above ambient) in replicate, 100 m flow-through flumes located in unenriched Everglades National Park. After the first six months of dosing, the concentration and standing stock of phosphorus increased in the surface water, periphyton, and flocculent detrital layer, but not in the soil or macrophytes. Of the ecosystem components measured, total P concentration increased the most in the floating periphyton mat (30 μg L−1: mean = 1916 μg P g−1, control: mean = 149 μg P g−1), while the flocculent detrital layer stored most of the accumulated P (30 μg L−1: mean = 1.732 g P m−2, control: mean = 0.769 g P m−2). Significant short-term responses of P concentration and standing stock were observed primarily in the high dose (30 μg L−1 above ambient) treatment. In addition, the biomass and estimated P standing stock of aquatic consumers increased in the 30 and 5 μg L−1 treatments. Alterations in P concentration and standing stock occurred only at the upstream ends of the flumes nearest to the point source of added nutrient. The total amount of P stored by the ecosystem within the flume increased with P dosing, although the ecosystem in the flumes retained only a small proportion of the P added over the first six months. These results indicate that oligotrophic Everglades wetlands respond rapidly to short-term, low-level P enrichment, and the initial response is most noticeable in the periphyton and flocculent detrital layer.

New fuel cells system for portable application with low temperature cofired ceramic (LTCC) technology

Miniature direct methanol fuel cells (DMFCs) are promising micro power sources for portable appliction. Low temperature cofired ceramic (LTCC), a competitive technology for current MEMS based fabrication, provides cost-effective mass manufacturing route for miniature DMFCs. Porous silver tape is adapted as electrodes to replace the traditional porous carbon electrodes due to its compatibility to LTCC processing and other electrochemical advantages. Electrochemical evaluation of silver under DMFCs operating conditions demonstrated that silver is a good electrode for DMFCs because of its reasonable corrosion resistance, low passivating current, and enhanced catalytic effect. Two catalyst loading methods (cofiring and postfiring) of the platinum and ruthenium catalysts are evaluated for LTCC based processing. The electrochemical analysis exhibits that the cofired path out-performs the postfiring path both at the anode and cathode. The reason is the formation of high surface area precipitated whiskers. Self-constraint sintering is utilized to overcome the difficulties of the large difference of coefficient of thermal expansion (CTE) between silver and LTCC (Dupont 951) tape during cofiring. The graphite sheet employed as a cavity fugitive insert guarantees cavity dimension conservation. Finally, performance of the membrane electrode assembly (MEA) with the porous silver electrode in the regular graphite electrode based cell and the integrated cofired cell is measured under passive fuel feeding condition. The MEA of the regular cell performs better as the electrode porosity and temperature increased. The power density of 10 mWcm-2 was obtained at ambient conditions with 1M methanol and it increased to 16 mWcm -2 at 50°C from an open circuit voltage of 0.58V. For the integrated prototype cell, the best performance, which depends on the balance methanol crossover and mass transfer at different temperatures and methanol concentrations, reaches 1.13 mWcm-2 at 2M methanol solution at ambient pressure. The porous media pore structure increases the methanol crossover resistance. As temperature increased to 60°C, the device increases to 2.14 mWcm-2.

Cascading ecological effects of low-level phosphorus enrichment in the Florida Everglades

Few studies have examined long-term ecological effects of sustained low-level nutrient enhancement on wetland biota. To determine sustained effects of phosphorus (P) addition on Everglades marshes we added P at low levels (5, 15, and 30 µg L-1 above ambient) for 5 yr to triplicate 100-m flow-through channels in pristine marsh. A cascade of ecological responses occurred in similar sequence among treatments. Although the rate of change increased with dosing level, treatments converged to similar enriched endpoints, characterized most notably by a doubling of plant biomass and elimination of native, calcareous periphyton mats. The full sequence of biological changes occurred without an increase in water total P concentration, which remained near ambient levels until Year 5. This study indicates that Everglades marshes have a near-zero assimilative capacity for P without a state change, that ecosystem responses to enrichment accumulate over time, and that downstream P transport mainly occurs through biota rather than the water column.

Testing the Efficacy of Panax ginseng Extract as a Novel Anti-Quorum Sensing Agent

Quorum sensing (QS) is the phenomenon by which microorganisms regulate gene expression in response to cell-population density. These microorganisms synthesize and secrete small molecules known as autoinducers that increase in concentration as a function of cell density. Once the cell detects the minimal threshold concentration of an autoinducer, the gene expression is altered accordingly. Although the cellular circuitry responding to QS is relatively conserved, the cellular processes regulated by QS, for example are conjugation, motility, sporulation, biofilm formation and production of virulence factors importamt for infection. Since many pathogens have developed resistance to available antibiotics, novel therapies are required to further treat these pathogens. Inhibitors of QS are potential therapeutic candidates since they would inhibit virulence without selecting for antibiotic resistant strains. Previous studies have shown that the Chinese herb Panax ginseng contains compounds that inhibit QS activity. To further characterize this activity, a highly sensitive quantitative liquid assay was developed using a Chromobacterium violaceum AHL mutant, CV026 as a biomonitor strain. After confirmation that P. ginseng aqueous extracts had anti-QS activity, the extract was fractionated on a reverse phase C18 Sep Pak column. Most anti-QS activity was present in the flow through, and compounds eluted with ten percent acetonitrile in water antibacterial activity. We thus conclude that P. ginseng has anti-QS activity and the active compound is water soluble. Compounds from P. ginseng, could be used as a lead structure to design compound with higher anti-QS activity for therapeutic treatments.

Novel methodology for in situ carbon dioxide enrichment of benthic ecosystems

Future climate change will likely represent a major stress to shallow aquatic and coastal marine communities around the world. Most climate change research, particularly in regards to increased pCO2 and ocean acidification, relies on ex situ mesocosm experimentation, isolating target organisms from their environment. Such mesocosms allow for greater experimental control of some variables, but can often cause unrealistic changes in a variety of environmental factors, leading to “bottle effects.” Here we present an in situ technique of altering dissolved pCO2within nearshore benthic communities (e.g., macrophytes, algae, and/or corals) using submerged clear, open-top chambers. Our technique utilizes a flow-through design that replicates natural water flow conditions and minimizes caging effects. The clear, open-top design additionally ensures that adequate light reaches the benthic community. Our results show that CO2 concentrations and pH can be successfully manipulated for long durations within the open-top chambers, continuously replicating forecasts for the year 2100. Enriched chambers displayed an average 0.46 unit reduction in pH as compared with ambient chambers over a 6-month period. Additionally, CO2 and HCO3 – concentrations were all significantly higher within the enriched chambers. We discuss the advantages and disadvantages of this technique in comparison to other ex situ mesocosm designs used for climate change research.

Importance of storm events in controlling ecosystem structure and function in a Florida gulf coast estuary

From 8/95 to 2/01, we investigated the ecological effects of intra- and inter-annual variability in freshwater flow through Taylor Creek in southeastern Everglades National Park. Continuous monitoring and intensive sampling studies overlapped with an array of pulsed weather events that impacted physical, chemical, and biological attributes of this region. We quantified the effects of three events representing a range of characteristics (duration, amount of precipitation, storm intensity, wind direction) on the hydraulic connectivity, nutrient and sediment dynamics, and vegetation structure of the SE Everglades estuarine ecotone. These events included a strong winter storm in November 1996, Tropical Storm Harvey in September 1999, and Hurricane Irene in October 1999. Continuous hydrologic and daily water sample data were used to examine the effects of these events on the physical forcing and quality of water in Taylor Creek. A high resolution, flow-through sampling and mapping approach was used to characterize water quality in the adjacent bay. To understand the effects of these events on vegetation communities, we measured mangrove litter production and estimated seagrass cover in the bay at monthly intervals. We also quantified sediment deposition associated with Hurricane Irene's flood surge along the Buttonwood Ridge. These three events resulted in dramatic changes in surface water movement and chemistry in Taylor Creek and adjacent regions of Florida Bay as well as increased mangrove litterfall and flood surge scouring of seagrass beds. Up to 5 cm of bay-derived mud was deposited along the ridge adjacent to the creek in this single pulsed event. These short-term events can account for a substantial proportion of the annual flux of freshwater and materials between the mangrove zone and Florida Bay. Our findings shed light on the capacity of these storm events, especially when in succession, to have far reaching and long lasting effects on coastal ecosystems such as the estuarine ecotone of the SE Everglades.

Salinity and Chlorophyll a as Performance Measures to Rehabilitate a Mangrove-Dominated Deltaic Coastal Region: the Ciénaga Grande de Santa Marta–Pajarales Lagoon Complex, Colombia

Salinity, water temperature, and chlorophyll a (chl-a) biomass were used as performance measures in the period 1999–2001 to evaluate the effect of a hydrological rehabilitation project in the Ciénaga Grande de Santa Marta (CGSM)–Pajarales lagoon complex, Colombia where freshwater diversions were initiated in 1995 and completed in 1998. The objective of this study was to evaluate how diversions of freshwater into previously hypersaline (>80) environments changed the spatial and temporal distribution of environmental characteristics. Following the diversion, 19 surveys and transects using a flow-through system were surveyed in the CGSM–Pajarales complex to continuously measure selected water quality parameters. Geostatistical analysis indicates that hydrology and salinity regimes and water circulation patterns in the CGSM lagoon are largely controlled by freshwater discharge from the Fundacion, Aracataca, and Sevilla Rivers. Residence times in the CGSM lagoon were similar before (15.5 ± 3.8 days) and after (14.2 ± 2.0 days) the rehabilitation project and indicated that the system is flushed regularly. In contrast, chl-a biomass was highly variable in the CGSM–Pajarales lagoon complex and not related to discharge patterns. Mean annual chl-a biomass (44–250 μg L−1) following the diversion project was similar to values recorded since the 1980s and still remains among the highest reported in coastal systems around the world owing to its unique hydrology regulated by the Magdalena River and Sierra Nevada de Santa Marta watersheds and the high teleconnection to the El Niño Southern Oscillation (ENSO). Our results confirm that the reduction in salinity in the CGSM lagoon and Pajarales complex during 1999–2000 was largely driven by high precipitation (2500 mm) induced by the ENSO–La Niña rather than by the freshwater diversions.

Geochemical determination of the fate and transport of injected fresh wastewater to a deep saline aquifer

Deep well injection into non-potable saline aquifers of treated domestic wastewater has been used in Florida for decades as a safe and effective alternative to ocean outfall disposal. The objectives of this study were to determine the fate and transport of injected wastewater at two deep well injection sites in Miami Dade County, Florida, USA. Detection of ammonium in the Middle Confining units of the Floridan aquifer above the injection zone at both sites has been interpreted as evidence of upward migration of injected wastewater, posing a risk to underground sources of drinking water. Historical water quality data, including ammonia, chloride, temperature, and pH from existing monitoring wells at both sites from 1983 to 2008, major ions collected monthly from 2006 and 2008, and a synoptic sampling event for stable isotopes, tritium, and dissolved gases in 2008, were used to determine the source of ammonium in groundwater and possible migration pathways. Geochemical modeling was used to determine possible effects of injected wastewater on native water and aquifer matrix geochemistry. Injected wastewater was determined to be the source of elevated ammonium concentrations above ambient water levels, based on the results of major ion concentrations, tritium, dissolved noble gases and 15N isotopes analyses. Various possible fluid migration pathways were identified at the sites. Data for the south site suggest buoyancy-driven vertical pathways to overlying aquifers bypassing the confining units, with little mixing of injected wastewater with native water as it migrated upward. Once it is introduced into an aquifer, the injectate appeared to migrate advectively with the regional groundwater flow. Geochemical modeling indicated that CO 2-enriched injected wastewater allowed for carbonate dissolution along the vertical pathways, enhancing permeability along these flowpaths. At the north site, diffusive upward flow through the confining units or offsite vertical pathways were determined to be possible, however no evidence was detected for any on-site confining unit bypass pathway. No evidence was observed at either site of injected wastewater migration to the Upper Floridan aquifer, which is used as a municipal water supply and for aquifer storage and recovery.

Characteristics of surface-water flows in the ridge and slough landscape of Everglades National Park: implications for particulate transport

Over the last one hundred years, compartmentalization and water management activities have reduced water flow to the ridge and slough landscape of the Everglades. As a result, the once corrugated landscape has become topographically and vegetationally uniform. The focus of this study was to quantify variation in surface flow in the ridge and slough landscape and to relate flow conditions to particulate transport and deposition. Over the 2002–2003 and 2003–2004 wet seasons, surface velocities and particulate accumulation were measured in upper Shark River Slough in Everglades National Park. Landscape characteristics such as elevation, plant density and biomass also were examined to determine their impact on flow characteristics and material transport. The results of this study demonstrate that the release of water during the wet season not only increases water levels, but also increased flow speeds and particulate transport and availability. Further, flow speeds were positively and significantly correlated with water level thereby enhancing particulate transport in sloughs relative to ridges especially during peak flow periods. Our results also indicate that the distribution of biomass in the water column, including floating plants and periphyton, affects velocity magnitude and shape of vertical profiles, especially in the sloughs where Utricularia spp. and periphyton mats are more abundant. Plot clearing experiments suggest that the presence of surface periphyton and Utricularia exert greater control over flow characteristics than the identity (i.e., sawgrass or spike rush) or density of emergent macrophytes, two parameters frequently incorporated into models describing flow through vegetated canopies. Based on these results, we suggest that future modeling efforts must take the presence of floating biomass, such as Utricularia, and presence of periphyton into consideration when describing particulate transport.

Short-term changes in phosphorus storage in an oligotrophic Everglades wetland ecosystem receiving experimental nutrient enrichment

Natural, unenriched Evergladeswetlands are known to be limited by phosphorus(P) and responsive to P enrichment. However,whole-ecosystem evaluations of experimental Padditions are rare in Everglades or otherwetlands. We tested the response of theEverglades wetland ecosystem to continuous,low-level additions of P (0, 5, 15, and30 μg L−1 above ambient) in replicate,100 m flow-through flumes located in unenrichedEverglades National Park. After the first sixmonths of dosing, the concentration andstanding stock of phosphorus increased in thesurface water, periphyton, and flocculentdetrital layer, but not in the soil or macrophytes. Of the ecosystem components measured, total P concentration increased the most in the floating periphyton mat (30 μg L−1: mean = 1916 μg P g−1, control: mean =149 μg P g−1), while the flocculentdetrital layer stored most of the accumulated P(30 μg L−1: mean = 1.732 g P m−2,control: mean = 0.769 g P m−2). Significant short-term responsesof P concentration and standing stock wereobserved primarily in the high dose (30 μgL−1 above ambient) treatment. Inaddition, the biomass and estimated P standingstock of aquatic consumers increased in the 30and 5 μg L−1 treatments. Alterationsin P concentration and standing stock occurredonly at the upstream ends of the flumes nearestto the point source of added nutrient. Thetotal amount of P stored by the ecosystemwithin the flume increased with P dosing,although the ecosystem in the flumes retainedonly a small proportion of the P added over thefirst six months. These results indicate thatoligotrophic Everglades wetlands respondrapidly to short-term, low-level P enrichment,and the initial response is most noticeable inthe periphyton and flocculent detrital layer.

Everglades Ridge, Slough, and Tree Island Mosaics: Year 2 Annual Report

Status and history of the Ridge-Slough Mosaic The Florida Everglades is a large subtropical wetland with diverse hydrologic, edaphic, and vegetative characteristics. Historically, a significant portion of this system was a slow moving river originating from the Kissimmee River floodplain, flowing into the vast but shallow Lake Okeechobee, and draining south-southwest over extensive peatlands into Florida Bay (McVoy 2011). Human-induced alterations to the hydrologic regime, including reduction, stabilization, and impoundment of water flow through diversion and compartmentalization of water via canals and levees have degraded pre-drainage vegetation patterns and microtopographic structure (Davis and Ogden 1994, Ogden 2005, McVoy 2011). The Everglades peatland emerged 5,000 years ago with the stabilization of sea level at approximately current elevations (Loveless 1959, Gleason and Stone 1994). This, combined with subtropical rainfalls, allowed a vast mass of water to slowly flow over a limestone bedrock platform 160 km long and 50 km wide at a near uniform descent totaling about 6 m, ultimately reaching Florida Bay (Stephens 1956, Gleason and Stone 1994, McVoy 2011). Vegetation quickly colonized the area, and peat, in the absence of adequate respiration, accumulated on the limestone bedrock to a depth of 3-3.7 m (Gleason and Stone 1994, McVoy et al. 2011). The “River of Grass” referenced by Douglas (1947) alludes to the dually intertwined processes of the historic riverine nature of the Everglades and the vast sawgrass (Cladium jamaicense) communities that have dominated the landscape for about the last 1,000 years (Bernhardt and Willard 2009).

Natural organic matter and colloid-facilitated arsenic transport and transformation in porous soil media

Prediction of arsenic transport and transformation in soil environment requires understanding the transport mechanisms and proper estimation of arsenic partitioning tong all three phases in soil/aquifer systems: mobile colloids, mobile soil solution, and immobile soil solids. The primary purpose of this research is to study natural dissolved organic matter (DOM)/colloid-facilitated transport of arsenic and understand the role of soil derived carriers in the transport and transformation of both inorganic and organoarsenicals in soils. ^ DOM/colloid facilitated arsenic transport and transformation in porous soil media were investigated using a set of experimental approaches including batch experiment, equilibrium membrane dialysis experiment and column experiment. Soil batch experiment was applied to investigate arsenic adsorption on a variety of soils with different characteristics; Equilibrium membrane dialysis was employed to determine the 'free' and 'colloid-bound/complexed' arsenic in water extracts of chosen soils; Column experiments were also set up in the laboratory to simulate arsenic transport and transformation through golf course soils in the presence and absence of soil-derived dissolved substances. ^ The experimental results revealed that organic matter amendments effectively reduced soil arsenic adsorption. The majority of arsenic present in the soil extracts was associated with small substances of molecular weight (MW) between 500 and 3,500 Da, Only a small fraction of arsenic was associated with higher MW substances (MW > 3,500 Da), which was operationally defined as colloidal part in this study. The association of arsenic and DOM in the soil extracts strongly affected arsenic bioavailability, arsenic transport and transformation in soils. The results of column experiments revealed arsenic complicated behavior with various processes occurring in soils studied, including: soil arsenic' adsorption, facilitated arsenic transportation by dissolved substances presented in soil extracts and microorganisms involved arsenic species transformation. ^ Soil organic matter amendments effectively reduce soil arsenic adsorption capability either by scavenging 'soil arsenic adsorption sites or by interactions between arsenic species and dissolved organic chemicals in soil solution. Close attention must be paid for facilitated arsenic transport by dissolved substances presented in soil solution and microorganisms involved arsenic species transformation in arsenic-contaminated soils.^