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.

Dynamics of ecosystem metabolism and flocculent detritus transport in estuarine Taylor River

Estuaries and estuarine wetlands are ecologically and societally important systems, exhibiting high rates of primary production that fuel offshore secondary production. Hydrological processes play a central role in shaping estuarine ecosystem structure and function by controlling nutrient loading and the relative contributions of marine and terrestrial influences on the estuary. The Comprehensive Everglades Restoration Plan includes plans to restore freshwater delivery to Taylor Slough, a shallow drainage basin in the southern Everglades, ultimately resulting in increased freshwater flow to the downstream Taylor River estuary. The existing seasonal and inter-annual variability of water flow and source in Taylor River affords the opportunity to investigate relationships between ecosystem function and hydrologic forcing. Estimates of aquatic ecosystem metabolism, derived from free-water, diel changes in dissolved oxygen, were combined with assessments of wetland flocculent detritus quality and transport within the context of seasonal changes in Everglades hydrology. Variation in ecosystem gross primary production and respiration were linked to seasonal changes in estuarine water quality using multiple autoregression models. Furthermore, Taylor River was observed to be net heterotrophic, indicating that an allochthonous source of carbon maintained ecosystem respiration in excess of autochthonous primary production. Wetland-derived detritus appears to be an important vector of energy and nutrients across the Everglades landscape; and in Taylor River, is seasonally flushed into ponded segments of the river where it is then respired. Lastly, seasonal water delivery appears to govern feedbacks regulating water column phosphorus availability in the Taylor River estuary.

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.

Historic primary producer communities linked to water quality and hydrologic changes in the northern Everglades

The northern Everglades Water Conservation Areas have experienced recent ecological shifts in primary producer community structure involving marl periphyton mats and dense Typha-dominated macrophyte stands. Multiple investigations have identified phosphorus (P) as a driver of primary producer community structure, but effects of water impoundment beginning in the 1950s and changes in water hardness [e.g., (CaCO3)] have also been identified as a concern. In an effort to understand pre-1950, primary producer community structure and identify community shifts since 1950, we measured pigment proxies on three sediment cores collected in Water Conservation Area-2A (WCA-2A) along a phosphorus enrichment gradient. Photosynthetic pigments, sediment total phosphorus content (TP), organic matter, total organic carbon and nitrogen were used to infer historic primary producer communities and changes in water quality and hydrology regulating those communities. Excess 210Pb was used to establish historic dates for the sediment cores. Results indicate the northern area of WCA-2A increased marl deposition and increased algal abundance ca. 1920. This increase in (presumably) calcareous periphyton before intensive agriculture and impoundment suggest canal-derived calcium inputs and to some extent early drainage effects played a role in initiating this community shift. The northern area community then shifted to Typha dominance around 1965. The areas to the south in WCA-2A experienced increased marl deposition and algal abundance around or just prior to 1950s impoundment, the precise timing limited by core age resolution. Continued increases in algal abundance were evident after 1950, coinciding with impoundment and deepening of canals draining into WCA-2A, both likely increasing water mineral and nutrient concentrations. The intermediate site developed a Typha-dominated community ca. 1995 while the southern-most core site WCA-2A has yet to develop Typha dominance. Numerous studies link sediment TP >650 mg P/kg to marsh habitat degradation into Typha-dominance. The northern and intermediate cores where Typha is currently support this previous research by showing a distinct shift in the sediment record to Typha dominance corresponding to sediment TP between 600 and 700 mg P/kg. These temporal and spatial differences are consistent with modern evidence showing water-column gradients in mineral inputs (including Ca, carbonates, and phosphorus) altering primary producer community structure in WCA-2A, but also suggest hydroperiod has an effect on the mechanisms regulating periphyton development and Typha dominance.

Impact of Water Supply and Sanitation on Diarrheal Morbidity among Young Children in the Socioeconomic and Cultural Context of Rwanda (Africa)

This project studied the frequency and of water contamination at the source, during transportation, and at home to determine the causes of contamination and its impact on the health of children aged 0 to 5 years. The methods used were construction of the infrastructure for three sources of potable water, administration of a questionnaire about socioeconomic status and sanitation behavior, anthropometric measurement of children, and analysis of water and feces. The contamination, first thought to be only a function of rainfall, turned out to be a very complex phenomenon. Water in homes was contaminated (43.4%) with more than 1100 total coliforms/100 ml due to the use of unclean utensils to transport and store water. This socio-economic and cultural problem should be ad- dressed with health education about sanitation, The latrines (found in 43.8% of families) presented a double-edged problem. The extremely high population density reduced the surface area of land per family, which resulted in a severe nutritional deficit (15% of the children) affecting mainly young children, rendering them more susceptible to diarrhea (three episodes/child/year).

Development of a mathematical and computer model to assess water quality impacts of hazardous compounds from minor gasoline spills in the inter-tidal zone of the Miami River, Florida

The objective of this study was to develop a model to predict transport and fate of gasoline components of environmental concern in the Miami River by mathematically simulating the movement of dissolved benzene, toluene, xylene (BTX), and methyl-tertiary-butyl ether (MTBE) occurring from minor gasoline spills in the inter-tidal zone of the river. Computer codes were based on mathematical algorithms that acknowledge the role of advective and dispersive physical phenomena along the river and prevailing phase transformations of BTX and MTBE. Phase transformations included volatilization and settling. ^ The model used a finite-difference scheme of steady-state conditions, with a set of numerical equations that was solved by two numerical methods: Gauss-Seidel and Jacobi iterations. A numerical validation process was conducted by comparing the results from both methods with analytical and numerical reference solutions. Since similar trends were achieved after the numerical validation process, it was concluded that the computer codes algorithmically were correct. The Gauss-Seidel iteration yielded at a faster convergence rate than the Jacobi iteration. Hence, the mathematical code was selected to further develop the computer program and software. The model was then analyzed for its sensitivity. It was found that the model was very sensitive to wind speed but not to sediment settling velocity. ^ A computer software was developed with the model code embedded. The software was provided with two major user-friendly visualized forms, one to interface with the database files and the other to execute and present the graphical and tabulated results. For all predicted concentrations of BTX and MTBE, the maximum concentrations were over an order of magnitude lower than current drinking water standards. It should be pointed out, however, that smaller concentrations than the latter reported standards and values, although not harmful to humans, may be very harmful to organisms of the trophic levels of the Miami River ecosystem and associated waters. This computer model can be used for the rapid assessment and management of the effects of minor gasoline spills on inter-tidal riverine water quality. ^

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

High-resolution cyclostratigraphy of the Cenomanian-Turonian interval: Paleoceanographic implications, a comparison between deep sea drilling project site 386 and 144

The purpose of this study was to determine the extent to which oceanic anoxic events (OAE's) are recorded in deep-water deposits of the former western Tethyan Sea, by investigating the Cenomanian-Turonian time interval characterized by the worldwide OAE 2 event. The study improved our knowledge of the possible controlling mechanisms that triggered this event at these sites, and furthered our understanding of this global phenomenon. This was examined by high-resolution, multi-proxy analyses of sediments at DSDP Sites 386 and 144, including sedimentology, scanning electron microscopy, stable isotopes, bulk and clay mineralogy, major and trace element geochemistry, biomarkers, and paleontological data. ^ The results provide a better stratigraphic resolution for the Cenomanian-Turonian, which allowed for more precise determination of chronologic boundaries, sedimentation rates at DSDP Site 386, and a more accurate calculation of the frequency of the cycles recorded in the sequence, which fall predominantly within the precession (∼23 kyr) and short eccentricity (∼100 kyr) ranges. The combined proxies allow assessment of the correlation of δ13Corg, and major and trace elements with the predominance of cyanobacteria. These organisms were the main producers of the organic matter during the dysoxic and euxinic conditions of OAE 2 at DSDP Site 386. A huge amount of microcrystalline quartz of eolian origin is also associated with OAE 2. The geochemical proxies further provide evidence that OAE 2 was linked to increased volcanism outside the deep water of the proto-Atlantic Ocean. The clays in the Turonian sediments are terrigenous and derived predominantly from eolian transport. Comparing DSDP Site 386 and 144 with stratotype sections, the δ13C org and TOC data indicate that OAE 2 seems diachronous throughout the proto-Atlantic Ocean. ^ This study concludes that the development of anoxic conditions in the deep water of the Atlantic during the latest Cenomanian-Turonian resulted from a combination of factors related to local oceanic setting and mitigated by global tectonism and climate. The data provide a more comprehensive view of the interacting factors that led to sustained high productivity of the cyanobacteria and photosynthetic protists that produced organic-carbon-rich deposits in the world's oceans. ^

An integrated flow and transport model to study the impact of mercury remediation strategies for East Fork Poplar Creek watershed, Oak Ridge, Tennessee

An integrated flow and transport model using MIKE SHE/MIKE 11 software was developed to predict the flow and transport of mercury, Hg(II), under varying environmental conditions. The model analyzed the impact of remediation scenarios within the East Fork Poplar Creek watershed of the Oak Ridge Reservation with respect to downstream concentration of mercury. The numerical simulations included the entire hydrological cycle: flow in rivers, overland flow, groundwater flow in the saturated and unsaturated zones, and evapotranspiration and precipitation time series. Stochastic parameters and hydrologic conditions over a five year period of historical hydrological data were used to analyze the hydrological cycle and to determine the prevailing mercury transport mechanism within the watershed. Simulations of remediation scenarios revealed that reduction of the highly contaminated point sources, rather than general remediation of the contaminant plume, has a more direct impact on downstream mercury concentrations.

Water resources assessment and geographic information system (GIS)-based stormwater runoff estimates for artificial recharge of freshwater aquifers in New Providence, Bahamas

The Bahamas is a small island nation that is dealing with the problem of freshwater shortage. All of the country’s freshwater is contained in shallow lens aquifers that are recharged solely by rainfall. The country has been struggling to meet the water demands by employing a combination of over-pumping of aquifers, transport of water by barge between islands, and desalination of sea water. In recent decades, new development on New Providence, where the capital city of Nassau is located, has created a large area of impervious surfaces and thereby a substantial amount of runoff with the result that several of the aquifers are not being recharged. A geodatabase was assembled to assess and estimate the quantity of runoff from these impervious surfaces and potential recharge locations were identified using a combination of Geographic Information Systems (GIS) and remote sensing. This study showed that runoff from impervious surfaces in New Providence represents a large freshwater resource that could potentially be used to recharge the lens aquifers on New Providence.

Effects of climate variability on transparency and thermal structure in subtropical, monomictic Lake Annie, Florida

Lake Annie is a small (37 ha), relatively deep (21 m) sinkhole lake on the Lake Wales Ridge (LWR) of central Florida with a long history of study, including monthly limnological monitoring since June, 1983. The record shows high variability in Secchi disc transparency, which ranged from < 1 to 15 m with a trend toward decreasing values over the latter decade of record. We examined available regional meteorological, groundwater and limnological data to determine the drivers and thermal consequences of variability in water transparency. While total nutrient concentrations and chlorophyll-a were highest during years of low transparency, stepwise regression showed that none of these had a signifi cant effect on transparency after water color was taken into account. Repeated years of high precipitation between 1993–2005 caused an increase in water table height, increasing the transport of dissolved substances from the vegetated watershed into the lake. Groundwater stage explained 73 % of the interannual variability in water transparency. Transparency, in turn, explained 85 % of the interannual variability in the heat budget for the lake, which ranged from 1.8 × 108 to 4.1 × 108 Joules m–2 yr–1, encompassing the range reported across Florida lakes. While surface water temperature was not affected by transparency, depths below 5 m warmed faster during the stratifi ed period during years having a lower rate of light extinction. We show that an increase in precipitation of 20 cm per year reduces the depth of the summer euphotic zone and thermocline by 1.9 and 1.6 m, respectively, and causes a 1-month reduction in the duration of winter mixing in this monomictic lake. Because biota have been shown to respond to shifts in light and heat distribution of much smaller magnitude than exhibited here, our work suggests that subtle changes in precipitation linked to climate fl uctuations may have signifi cant physical as well as biotic consequences.

Multidecadal climate oscillations detected in a transparency record from a subtropical Florida lake

Synchronous interannual variability in water transparency observed in neighboring lakes has been linked to regional precipitation and resultant runoff of dissolved organic material, but many climate forcings oscillate over time scales longer than most limnological records can detect. A strong relationship (R2 5 0.86) between transparency and the previous two years’ rainfall and lake stage in a 25-yr record from a Florida lake enabled us to hindcast transparency from a longer 75-yr record of rainfall and lake stage. Predictions revealed a ,30-yr cycle in transparency linked to the Atlantic Multidecadal Oscillation (AMO). Transparency was greatest (4–8 m) in the cool phase of the AMO (,1962–1993) associated with below-average rainfall in south Florida and lowest (0.1– 3.0 m) during two warm phases (,1932–1961, 1994–present) associated with above-average, but more variable, annual rainfall. Models that predict effects of large-scale hydrologic restoration projects on solute export from South Florida’s expansive wetlands need to account for recent entry into a warm AMO phase, where teleconnections between the AMO phases and runoff are opposite of those shown for the U.S. interior.

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.

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.