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.

Larger foraminiferal biostratigraphy, systematics and paleoenvironments of the Avon Park Formation and Ocala Limestone, Highlands County, Florida

This study investigates the use of larger foraminifera in determining the biostratigraphy of the Avon Park Formation and the Ocala Limestone in central Florida. Sedimentary rocks of the Avon Park Formation are the oldest exposed deposits in the state of Florida, and together with the Ocala Limestone comprise a part of the confining unit of the Floridan Aquifer, a major source of Florida's water supply. ^ Material from the ROMP 29A core collected by the U.S. Geological Survey was evaluated and compared to previous studies of the biostratigraphy of the formations. The larger foraminifera of the Avon Park Formation were examined in thin section, and those of the Ocala Limestone were free specimens. The larger foraminifera from both units were described and identified, and the biostratigraphy determined. The morphological features of the larger foraminifera of the Ocala Limestone were measured and analyzed at various depths within the ROMP 29A core.^ The Avon Park Formation contains predominantly the shallow-water, conical foraminifera Fallotella cookei, Fallotella floridana, Pseudochrysalidina floridana, Coleiconus christianaensis, Coleiconus sp. A, Coskinolina sp. A, Coskinolina sp. B, Fallotella sp. A, Fallotella sp. B, Fabularia vaughani and larger miliolids. ^ The Ocala Limestone contains a different, deeper water assemblage that included the larger foraminifera Heterostegina ocalana, Lepidocyclina ocalana varieties, Lepidocyclina chaperi, Lepidocyclina pustulosa, Nummulites willcoxi, Nummulites striatoreticulatus, Nummulites floridensis and Pseudophragmina spp. A, B, and C. The age of the Avon Park Formation was corroborated by the occurrence of the biomarker echinoid Neolaganum dalli as Eocene, and the Ocala Limestone also contained Eocene larger foraminifera with Eocene to possibly Oligocene calcareous nannofossils. The distribution of the larger foraminifera of the Avon Park Formation was correlated with the subtidal and peritidal zones of the continental shelf. Analyses of variance showed that the changes in measurements of the morphology in Heterostegina ocalana, Lepidocyclina spp. and Nummulites spp. were correlated with change in the depositional environments.^

Pink shrimp as an indicator for restoration of everglades ecosystems

The pink shrimp, Farfantepenaeus duorarum, familiar to most Floridians as either food or bait shrimp, is ubiquitous in South Florida coastal and offshore waters and is proposed as an indicator for assessing restoration of South Florida's southern estuaries: Florida Bay, Biscayne Bay, and the mangrove estuaries of the lower southwest coast. Relationships between pink shrimp and salinity have been determined in both field and laboratory studies. Salinity is directly relevant to restoration because the salinity regimes of South Florida estuaries, critical nursery habitat for the pink shrimp, will be altered by changes in the quantity, timing, and distribution of freshwater inflow planned as part of the Comprehensive Everglades Restoration Project (CERP). Here we suggest performance measures based on pink shrimp density (number per square meter) in the estuaries and propose a restoration assessment and scoring scheme using these performance measures that can readily be communicated to managers, policy makers, and the interested public. The pink shrimp is an appropriate restoration indicator because of its ecological as well as its economic importance and also because scientific interest in pink shrimp in South Florida has produced a wealth of information about the species and relatively long time series of data on both juveniles in estuarine nursery habitats and adults on the fishing grounds. We suggest research needs for improving the pink shrimp performance measure.

A conceptual model of ecological interactions in the mangrove estuaries of the Florida Everglades

A brackish water ecotone of coastal bays and lakes, mangrove forests, salt marshes, tidal creeks, and upland hammocks separates Florida Bay, Biscayne Bay, and the Gulf of Mexico from the freshwater Everglades. The Everglades mangrove estuaries are characterized by salinity gradients that vary spatially with topography and vary seasonally and inter-annually with rainfall, tide, and freshwater flow from the Everglades. Because of their location at the lower end of the Everglades drainage basin, Everglades mangrove estuaries have been affected by upstream water management practices that have altered the freshwater heads and flows and that affect salinity gradients. Additionally, interannual variation in precipitation patterns, particularly those caused to El Nin˜o events, control freshwater inputs and salinity dynamics in these estuaries. Two major external drivers on this system are water management activities and global climate change. These drivers lead to two major ecosystem stressors: reduced freshwater flow volume and duration, and sea-level rise. Major ecological attributes include mangrove forest production, soil accretion, and resilience; coastal lake submerged aquatic vegetation; resident mangrove fish populations; wood stork (Mycteria americana) and roseate spoonbill (Platelea ajaja) nesting colonies; and estuarine crocodilian populations. Causal linkages between stressors and attributes include coastal transgression, hydroperiods, salinity gradients, and the ‘‘white zone’’ freshwater/estuarine interface. The functional estuary and its ecological attributes, as influenced by sea level and freshwater flow, must be viewed as spatially dynamic, with a possible near-term balancing of transgression but ultimately a long-term continuation of inland movement. Regardless of the spatio-temporal timing of this transgression, a salinity gradient supportive of ecologically functional Everglades mangrove estuaries will be required to maintain the integrity of the South Florida ecosystem.

Development of an integrated surface and subsurface model of Everglades National Park

An integrated surface-subsurface hydrological model of Everglades National Park (ENP) was developed using MIKE SHE and MIKE 11 modeling software. The model has a resolution of 400 meters, covers approximately 1050 square miles of ENP, includes 110 miles of drainage canals with a variety of hydraulic structures, and processes hydrological information, such as evapotranspiration, precipitation, groundwater levels, canal discharges and levels, and operational schedules. Calibration was based on time series and probability of exceedance for water levels and discharges in the years 1987 through 1997. Model verification was then completed for the period of 1998 through 2005. Parameter sensitivity in uncertainty analysis showed that the model was most sensitive to the hydraulic conductivity of the regional Surficial Aquifer System, the Manning's roughness coefficient, and the leakage coefficient, which defines the canal-subsurface interaction. The model offers an enhanced predictive capability, compared to other models currently available, to simulate the flow regime in ENP and to forecast the impact of topography, water flows, and modifying operation schedules.

Chilling damage in a changing climate in coastal landscapes of the subtropical zone: a case study from south Florida

Extensive portions of the southern Everglades are characterized by series of elongated, raised peat ridges and tree islands oriented parallel to the predominant flow direction, separated by intervening sloughs. Tall herbs or woody species are associated with higher elevations and shorter emergent or floating species are associated with lower elevations. The organic soils in this “Ridge-and-Slough” landscape have been stable over millennia in many locations, but degrade over decades under altered hydrologic conditions. We examined soil, pore water, and leaf phosphorus (P) and nitrogen (N) distributions in six Ridge and Slough communities in Shark Slough, Everglades National Park. We found P enrichment to increase and N to decrease monotonically along a gradient from the most persistently flooded sloughs to rarely flooded ridge environments, with the most dramatic change associated with the transition from marsh to forest. Leaf N:P ratios indicated that the marsh communities were strongly P-limited, while data from several forest types suggested either N-limitation or co-limitation by N and P. Ground water stage in forests exhibited a daytime decrease and partial nighttime recovery during periods of surface exposure. The recovery phase suggested re-supply from adjacent flooded marshes or the underlying aquifer, and a strong hydrologic connection between ridge and slough. We therefore developed a simple steady-state model to explore a mechanism by which a phosphorus conveyor belt driven by both evapotranspiration and the regional flow gradient can contribute to the characteristic Ridge and Slough pattern. The model demonstrated that evapotranspiration sinks at higher elevations can draw in low concentration marsh waters, raising local soil and water P concentrations. Focusing of flow and nutrients at the evapotranspiration zone is not strong enough to overcome the regional gradient entirely, allowing the nutrient to spread downstream and creating an elongated concentration plume in the direction of flow. Our analyses suggest that autogenic processes involving the effects of initially small differences in topography, via their interactions with hydrology and nutrient availability, can produce persistent physiographic patterns in the organic sediments of the Everglades.

Adsorption and desorption of phosphate on limestone in experiments simulating seawater intrusion

This study investigates the potential release of from carbonate aquifers exposed to seawater intrusion. Adsorption and desorption of in the presence of deionized water (DIW) and seawater were conducted on a large block of Pleistocene age limestone to simulate the effects of seawater intrusion into a coastal carbonate aquifer at the laboratory scale. The limestone showed strong adsorption of in DIW, while adsorption was significantly less in the presence of seawater. Dissolution of CaCO3 was found to prevent adsorption at salinities less than 30 psu. Adsorption of was limited at higher salinities (30–33 psu), due to competition with ions for adsorption sites. At a salinity3 precipitated. Concentrations of between 2 and 5 μmol/L were released by desorption when the limestone was exposed to seawater. The results of this study suggest that as seawater intrudes into an originally freshwater coastal aquifer, adsorbed may be released into the groundwater. Consequently, adsorbed is expected to be released from coastal carbonate aquifers world-wide as sea level continues to rise exposing more of the freshwater aquifer to seawater.

Coastal groundwater discharge – an additional source of phosphorus for the oligotrophic wetlands of the Everglades

In this manuscript we define a new term we call coastal groundwater discharge (CGD), which is related to submarine groundwater discharge (SGD), but occurs when seawater intrudes inland to force brackish groundwater to discharge to the coastal wetlands. A hydrologic and geochemical investigation of both the groundwater and surface water in the southern Everglades was conducted to investigate the occurrence of CGD associated with seawater intrusion. During the wet season, the surface water chemistry remained fresh. Enhanced chloride, sodium, and calcium concentrations, indicative of brackish groundwater discharge, were observed in the surface water during the dry season. Brackish groundwaters of the southern Everglades contain 1–2.3μM concentrations of total phosphorus (TP). These concentrations exceed the expected values predicted by conservative mixing of local fresh groundwater and intruding seawater, which both have TPμM. The additional source of TP may be from seawater sediments or from the aquifer matrix as a result of water–rock interactions (such as carbonate mineral dissolution and ion exchange reactions) induced by mixing fresh groundwater with intruding seawater. We hypothesize that CGD maybe an additional source of phosphorus (a limiting nutrient) to the coastal wetlands of the southern Everglades.

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.

Student Leadership Summit, 10th Anniversary Time Capsule Burial Ceremony, February 2, 2013

Time capsule table at the Student Leadership Summit 10th Anniversary. Time capsule table at the Student Leadership Summit 10th Anniversary. The Annual FIU Student Leadership Summit is held each February on the Biscayne Bay Campus. The Summit is a one-day conference for current student leaders. The Summit offers our students the opportunity to learn from the vast expertise of our faculty and administrators, to share their leadership experiences with each other and to establish a network of support and cooperation within the university. On Feb. 2, 2013, we celebrated the 10th anniversary of holding the Student Leadership Summit. In honor of this occasion, we buried a time capsule containing materials from the day as well as messages from participants to the participants of 2023 when the time capsule is to be opened.

Student Leadership Summit, 10th Anniversary Time Capsule Burial Ceremony, February 2, 2013

Time capsule and burial plaque. The Annual FIU Student Leadership Summit is held each February on the Biscayne Bay Campus. The Summit is a one-day conference for current student leaders. The Summit offers our students the opportunity to learn from the vast expertise of our faculty and administrators, to share their leadership experiences with each other and to establish a network of support and cooperation within the university. On Feb. 2, 2013, we celebrated the 10th anniversary of holding the Student Leadership Summit. In honor of this occasion, we buried a time capsule containing materials from the day as well as messages from participants to the participants of 2023 when the time capsule is to be opened.

Student Leadership Summit, 10th Anniversary Time Capsule Burial Ceremony, February 2, 2013

Time capsule table at the Student Leadership Summit 10th Anniversary. The Annual FIU Student Leadership Summit is held each February on the Biscayne Bay Campus. The Summit is a one-day conference for current student leaders. The Summit offers our students the opportunity to learn from the vast expertise of our faculty and administrators, to share their leadership experiences with each other and to establish a network of support and cooperation within the university. On Feb. 2, 2013, we celebrated the 10th anniversary of holding the Student Leadership Summit. In honor of this occasion, we buried a time capsule containing materials from the day as well as messages from participants to the participants of 2023 when the time capsule is to be opened.

Comparison of community structure of mangrove and riprap rock habitats along urban coastlines of south Florida

Globally, mangrove ecosystems have substantially declined, largely a result of human impacts. Mangroves provide a number of ecosystem services such as shoreline stabilization and nursery habitat for fish species. As declines continue, many of these ecosystem services are lost or altered. The need for shoreline stabilization has become increasingly apparent when chronic erosion wear away coastlines once mangroves are removed. Limestone boulders called riprap have been employed to offset continued erosion associated with mangrove clearing. In urban coastal areas adjacent to Biscayne Bay, Florida, as much as 80 percent of mangroves have been lost. More recently, riprap has been used in conjunction with mangroves to restore wetlands throughout the Bay. This riprap-mangrove habitat provides structure for marine organisms to colonize. However, fish assemblages and benthic composition could vary between this hybridized habitat and natural mangrove systems. Comparisons of fish and benthic community structure were made, to determine if abundance, species richness, and overall diversity differed between the two habitat types. Visual census and benthic quadrat surveys were conducted in vi mangrove and mangrove-riprap sites within two regions of Biscayne Bay. Total fish abundance was greater in mangroves, but the effect of habitat type on species richness varied between regions. The community structure of fishes and benthic composition differed significantly between mangroves and riprap habitats. Because species composition is so distinct, it is likely that the two communities do no function in the same manner. In areas with cleared shorelines, it may be important to consider the function of added anthropogenic structure for ecological communities.

Hydro-physical characterization of media used in agricultural systems to develop the best management practices for operation of an environmentally sustainable agricultural enterprise

Florida is the second leading horticulture state in the United States with a total annual industry sale of over $12 Billion. Due to its competitive nature, agricultural plant production represents an extremely intensive practice with large amounts of water and fertilizer usage. Agrochemical and water management are vital for efficient functioning of any agricultural enterprise, and the subsequent nutrient loading from such agricultural practices has been a concern for environmentalists. A thorough understanding of the agrochemical and the soil amendments used in these agricultural systems is of special interest as contamination of soils can cause surface and groundwater pollution leading to ecosystem toxicity. The presence of fragile ecosystems such as the Everglades, Biscayne Bay and Big Cypress near enterprises that use such agricultural systems makes the whole issue even more imminent. Although significant research has been conducted with soils and soil mix, there is no acceptable method for determining the hydraulic properties of mixtures that have been subjected to organic and inorganic soil amendments. Hydro-physical characterization of such mixtures can facilitate the understanding of water retention and permeation characteristics of the commonly used mix which can further allow modeling of soil water interactions. The objective of this study was to characterize some of the locally and commercially available plant growth mixtures for their hydro-physical properties and develop mathematical models to correlate these acquired basic properties to the hydraulic conductivity of the mixture. The objective was also to model the response patterns of soil amendments present in those mixtures to different water and fertilizer use scenarios using the characterized hydro-physical properties with the help of Everglades-Agro-Hydrology Model. The presence of organic amendments helps the mixtures retain more water while the inorganic amendments tend to adsorb more nutrients due to their high surface area. The results of these types of characterization can provide a scientific basis for understanding the non-point source water pollution from horticulture production systems and assist in the development of the best management practices for the operation of environmentally sustainable agricultural enterprise