Integrated surface-ground water modeling in wetlands with improved methods to simulate vegetative resistance to flow

In topographically flat wetlands, where shallow water table and conductive soil may develop as a result of wet and dry seasons, the connection between surface water and groundwater is not only present, but perhaps the key factor dominating the magnitude and direction of water flux. Due to their complex characteristics, modeling waterflow through wetlands using more realistic process formulations (integrated surface-ground water and vegetative resistance) is an actual necessity. This dissertation focused on developing an integrated surface – subsurface hydrologic simulation numerical model by programming and testing the coupling of the USGS MODFLOW-2005 Groundwater Flow Process (GWF) package (USGS, 2005) with the 2D surface water routing model: FLO-2D (O’Brien et al., 1993). The coupling included the necessary procedures to numerically integrate and verify both models as a single computational software system that will heretofore be referred to as WHIMFLO-2D (Wetlands Hydrology Integrated Model). An improved physical formulation of flow resistance through vegetation in shallow waters based on the concept of drag force was also implemented for the simulations of floodplains, while the use of the classical methods (e.g., Manning, Chezy, Darcy-Weisbach) to calculate flow resistance has been maintained for the canals and deeper waters. A preliminary demonstration exercise WHIMFLO-2D in an existing field site was developed for the Loxahatchee Impoundment Landscape Assessment (LILA), an 80 acre area, located at the Arthur R. Marshall Loxahatchee National Wild Life Refuge in Boynton Beach, Florida. After applying a number of simplifying assumptions, results have illustrated the ability of the model to simulate the hydrology of a wetland. In this illustrative case, a comparison between measured and simulated stages level showed an average error of 0.31% with a maximum error of 2.8%. Comparison of measured and simulated groundwater head levels showed an average error of 0.18% with a maximum of 2.9%. The coupling of FLO-2D model with MODFLOW-2005 model and the incorporation of the dynamic effect of flow resistance due to vegetation performed in the new modeling tool WHIMFLO-2D is an important contribution to the field of numerical modeling of hydrologic flow in wetlands.

Seasonal plant water uptake patterns in the saline southeast Everglades ecotone

The purpose of this study was to determine the seasonal water use patterns of dominant macrophytes coexisting in the coastal Everglades ecotone. We measured the stable isotope signatures in plant xylem water of Rhizophora mangle, Cladium jamaicense, and Sesuvium portulacastrum during the dry (DS) and wet (WS) seasons in the estuarine ecotone along Taylor River in Everglades National Park, FL, USA. Shallow soilwater and deeper groundwater salinity was also measured to extrapolate the salinity encountered by plants at their rooting zone. Average soil water oxygen isotope ratios (δ 18O) was enriched (4.8 ± 0.2‰) in the DS relative to the WS (0.0 ± 0.1‰), but groundwater δ 18O remained constant between seasons (DS: 2.2 ± 0.4‰; WS: 2.1 ± 0.1‰). There was an inversion in interstitial salinity patterns across the soil profile between seasons. In the DS, shallow water was euhaline [i.e., 43 practical salinity units (PSU)] while groundwater was less saline (18 PSU). In the WS, however, shallow water was fresh (i.e., 0 PSU) but groundwater remained brackish (14 PSU). All plants utilized 100% (shallow) freshwater during the WS, but in the DS R. mangle switched to a soil–groundwater mix (δ 55% groundwater) while C. jamaicense and S. portulacastrum continued to use euhaline shallow water. In the DS, based on δ 18O data, the roots of R. mangle roots were exposed to salinities of 25.4 ± 1.4 PSU, less saline than either C. jamaicense(39.1 ± 2.2 PSU) or S. portulacastrum (38.6 ± 2.5 PSU). Although the salinity tolerance of C. jamaicense is not known, it is unlikely that long-term exposure to high salinity is conducive to the persistence of this freshwater marsh sedge. This study increases our ecological understanding of how water uptake patterns of individual plants can contribute to ecosystem levels changes, not only in the southeast saline Everglades, but also in estuaries in general in response to global sea level rise and human-induced changes in freshwater flows.

Black band disease of corals: Ecology and physiology of Phormidium corallyticum

Black band disease of corals consists of a microbial community dominated by the cyanobacteriurn Phormidium corallyticum. The disease primarily affects reef-framework coral species, Active black band disease continually opens up new substrate in reef environments by destroying coral tissue as the disease line advances across the surface of infected colonies. A field study was carried out to determine the abundance and distribution of black band disease on the reef building corals in the Florida Keys. During July of 1992 and 1993, up to 0.72% of coral colonies were infected with black band disease. Analysis of the distribution showed that the disease was clumped. Seasonal patters varied, with some coral colonies infected year round, others exhibiting reinfection from summer 1992 to summer 1993, and some colonies infected for one year only. Statistical analysis of black band disease incidence in relation to various environmental parameters revealed that black band disease was associated with relatively shallow water depths, higher temperatures, elevated nitrite levels, and decreased ortho-phosphate levels. Additional field studies determined recovery of scleractinian coral colonies damaged or killed through the activities of black band disease over a five-year period. These studies determined if the newly exposed substrate was recolonized through scleractinian recruitment, if there was overgrowth of the damaged areas by the formerly diseased colony, or if coral tissue destruction continued after the cessation of black band disease activity. Tissue loss continued on all coral colonies with only one colony exhibiting new tissue growth. The majority of recolonization was by non-reef-framework corals and octocorallians, limited recruitment by framework species was observed. Physiological studies of P. corallyticum were carried out to investigate the photosynthetic capacity of this cyanobacterium, and to determine if this species has the ability to fix dinitrogen. The results of this research demonstrated that P. corallyticum reaches maximum photosynthetic rates at very low light intensities (27.9 μE/m/sec), and that P. corallyticum is able to carry out oxygenic photosynthesis in the presence of sulfide, an ability that is uncommon in prokaryotic organisms. ^

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

The White Ibis and Wood Stork as indicators for restoration of the everglades ecosystem

Large numbers of colonially nesting herons, egrets, ibises, storks and spoonbills were one of the defining natural phenomena of the historical Everglades. Reproduction of these species has been tracked over at least a century, and some clear responses to dramatic anthropogenic hydrological alterations have been established. These include a marked decline in nesting populations of several species, and a movement of colonies away from the over-drained estuarine region. Ponding in a large portion of the freshwater marsh has favored species that hunt by sight in deep water (egrets, cf. 25–45 cm), while tactile feeders (ibises and storks) that depend on concentrated prey in shallow water (5–25 cm) have become proportionately much less common. There has been a marked increase in the interval between exceptionally large breeding aggregations of White Ibises (Eudocimus albus). Loss of short hydroperiod wetlands on the margins of the Everglades have delayed nest initiations 1–2 months by Wood Storks (Mycteria americana) resulting in poor nesting success. These responses are consistent with mechanisms that involve foraging, and the availability and production of prey animals, and each of the relationships is highly dependent on hydrology. Here, we define a group of characteristics about wading bird dynamics (= indicators) that collectively track the specific ecological relationships that supported ibises and storks in the past. We suggest four metrics as indicators of restoration success: timing of nesting by storks, the ratio of nesting ibises + storks to Great Egrets, the proportion of all nests located in the estuarine/freshwater ecotone, and the interval between years with exceptionally large ibis nestings. Each of these metrics has historical (e.g., predrainage) data upon which to base expectations for restoration, and the metrics have little measurement error relative to the large annual variation in numbers of nests. In addition to the strong scientific basis for the use of these indicators, wading birds are also a powerful tool for public communication because they have strong aesthetic appeal, and their ecological relationships with water are intuitively understandable. In the interests of communicating with the public and decision-makers, we integrate these metrics into a single-page annual “traffic-light” report card for wading bird responses. Collectively, we believe these metrics offer an excellent chance of detecting restoration of the ecosystem functions that supported historical wading bird nesting patterns.

Lower Aptian Comparative Stratigraphy of the Basco-Cantabrian Region (Spain) and Eastern Cordillera (Colombia): implications for local factors in the depositional record of Oceanic Anoxic Event 1a (OAE-1a)

An important episode of carbon sequestration, Oceanic Anoxic Event 1a (OAE-1a), characterizes the Lower Aptian worldwide, and is mostly known from deeper-water settings. The present work of two Lower Aptian deposits, Madotz (N Spain) and Curití Quarry (Colombia), is a multiproxy study that includes fossil assemblages, microfacies, X-ray diffraction bulk and clay mineralogy, elemental analyses (major, minor, trace elements), Rock-Eval pyrolysis, biomarkers, inorganic and organic carbon content, and stable carbon isotopes. The results provide baseline evidence of the local and global controlling environmental factors influencing OAE-1a in shallow-water settings. The data also improve our general understanding of the conditions under which organic-carbon-rich deposits accumulate. The sequence at Madotz includes four intervals (Unit 1; Subunits 2a, 2b and 2c) that overlap the times prior to, during and after the occurrence of OAE-1a. The Lower Unit 1(3m thick) is essentially siliciclastic, and Subunit 2a (20m) contains Urgonian carbonate facies that document abruptly changing platform conditions prior to OAE-1a. Subunit 2b (24.4 m) is a mixed carbonate-siliciclastic facies with orbitolinid-rich levels that coincides with OAE-1a δ13C stages C4-C6, and is coeval with the upper part of the Deshayesites forbesi ammonite zone. Levels with pyrite and the highest TOC values (0.4-0.97%), interpreted as accumulating under suboxic conditions, and are restricted to δ13C stages C4 and C5. The best development of the suboxic facies is at the level representing the peak of the transgression. Subunit 2c, within δ13C stage C7, shows a return of the Urgonian facies. The 23.35-m section at Curití includes a 6.3-m interval at the base of the Paja Formation dominated by organic-rich marlstones and shales lacking benthic fossils and bioturbation, with TOC values as high as 8.84%. The interval overlies a level containing reworked and phosphatized assemblages of middle Barremian to lowest Aptian ammonites. The range of values and the overall pattern of the δ13Corg (-22.05‰ to -20.47‰) in the 6.3m-interval is comparable with Lower Aptian δ13C stage C7. Thus, conditions of oxygen depletion at this site also occurred after Oceanic Anoxic Event-1a, which developed between carbon isotope stages C3 and C6. Both sites, Madotz and Curití, attest to the importance of terrigenous and nutrient fluxes in increasing OM productivity that led to episodic oxygen deficiency.

An investigation into equations for estimating water requirements and the development of new equations for predicting total water intake

The primary purpose of this study was to investigate agreement among five equations by which clinicians estimate water requirements (EWR) and to determine how well these equations predict total water intake (TWI). The Institute of Medicine has used TWI as a measure of water requirements. A secondary goal of this study was to develop practical equations to predict TWI. These equations could then be considered accurate predictors of an individual’s water requirement. ^ Regressions were performed to determine agreement between the five equations and between the five equations and TWI using NHANES 1999–2004. The criteria for agreement was (1) strong correlation coefficients between all comparisons and (2) regression line that was not significantly different when compared to the line of equality (x=y) i.e., the 95% CI of the slope and intercept must include one and zero, respectively. Correlations were performed to determine association between fat-free mass (FFM) and TWI. Clinically significant variables were selected to build equations for predicting TWI. All analyses were performed with SAS software and were weighted to account for the complex survey design and for oversampling. ^ Results showed that the five EWR equations were strongly correlated but did not agree with each other. Further, the EWR equations were all weakly associated to TWI and lacked agreement with TWI. The strongest agreement between the NRC equation and TWI explained only 8.1% of the variability of TWI. Fat-free mass was positively correlated to TWI. Two models were created to predict TWI. Both models included the variables, race/ethnicity, kcals, age, and height, but one model also included FFM and gender. The other model included BMI and osmolality. Neither model accounted for more than 28% of the variability of TWI. These results provide evidence that estimates of water requirements would vary depending upon which EWR equation was selected by the clinician. None of the existing EWR equations predicted TWI, nor could a prediction equation be created which explained a satisfactory amount of variance in TWI. A good estimate of water requirements may not be predicted by TWI. Future research should focus on using more valid measures to predict water requirements.^

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

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.

Sediment Dynamics of a Shallow Hypereutrophic lake: Lake Jesup, Florida, USA

Improved knowledge of sediment dynamics within a lake system is important for understanding lake water quality. This research was focused on an assessment of the vertical sediment flux in Lake Jesup, a shallow (1.3 m average depth) hypereutrophic lake of central Florida. Sediment dynamics were assessed at varying time scales (daily to weekly) to understand the transport of sediments from external forces; wind, waves, precipitation and/or runoff. Four stations were selected within the lake on the basis of water depth and the thicknesses of unconsolidated (floc) and consolidated sediments. At each of these stations, a 10:1 (length to diameter) high aspect ratio trap (STHA) was deployed to collect particulate matter for a one to two week period. The water and sediment samples were collected and analyzed for total carbon (TC), total phosphorus (TP) and total nitrogen (TN). Mass accumulation rates (MAR) collected by the traps varied from 77 to 418 g m-2 d-1 over seven deployments. TN, TP and TC sediment concentrations collected by the traps were consistently higher than the sediments collected by coring the lake bottom and is most likely associated with water column biomass. A yearly nutrient budget was determined from August 2009 to August 2010 with flux calculated as 2,033,882 mt yr-1.

Tracking rates of ecotone migration due to salt-water encroachment using fossil mollusks in coastal South Florida

We determined the rate of migration of coastal vegetation zones in response to salt-water encroachment through paleoecological analysis of mollusks in 36 sediment cores taken along transects perpendicular to the coast in a 5.5 km2 band of coastal wetlands in southeast Florida. Five vegetation zones, separated by distinct ecotones, included freshwater swamp forest, freshwater marsh, and dwarf, transitional and fringing mangrove forest. Vegetation composition, soil depth and organic matter content, porewater salinity and the contemporary mollusk community were determined at 226 sites to establish the salinity preferences of the mollusk fauna. Calibration models allowed accurate inference of salinity and vegetation type from fossil mollusk assemblages in chronologically calibrated sediments. Most sediments were shallow (20–130 cm) permitting coarse-scale temporal inferences for three zones: an upper peat layer (zone 1) representing the last 30–70 years, a mixed peat-marl layer (zone 2) representing the previous ca. 150–250 years and a basal section (zone 3) of ranging from 310 to 2990 YBP. Modern peat accretion rates averaged 3.1 mm yr)1 while subsurface marl accreted more slowly at 0.8 mm yr)1. Salinity and vegetation type for zone 1 show a steep gradient with freshwater communities being confined west of a north–south drainage canal constructed in 1960. Inferences for zone 2 (pre-drainage) suggest that freshwater marshes and associated forest units covered 90% of the area, with mangrove forests only present along the peripheral coastline. During the entire pre-drainage history, salinity in the entire area was maintained below a mean of 2 ppt and only small pockets of mangroves were present; currently, salinity averages 13.2 ppt and mangroves occupy 95% of the wetland. Over 3 km2 of freshwater wetland vegetation type have been lost from this basin due to salt-water encroachment, estimated from the mollusk-inferred migration rate of freshwater vegetation of 3.1 m yr)1 for the last 70 years (compared to 0.14 m yr)1 for the pre-drainage period). This rapid rate of encroachment is driven by sea-level rise and freshwater diversion. Plans for rehydrating these basins with freshwater will require high-magnitude re-diversion to counteract locally high rates of sea-level rise.

Factors Controlling Surface Water Flow in a Low-gradient Subtropical Wetland

Surface water flow patterns in wetlands play a role in shaping substrates, biogeochemical cycling, and ecosystem characteristics. This paper focuses on the factors controlling flow across a large, shallow gradient subtropical wetland (Shark River Slough in Everglades National Park, USA), which displays vegetative patterning indicative of overland flow. Between July 2003 and December 2007, flow speeds at five sites were very low (s−1), and exhibited seasonal fluctuations that were correlated with seasonal changes in water depth but also showed distinctive deviations. Stepwise linear regression showed that upstream gate discharges, local stage gradients, and stage together explained 50 to 90% of the variance in flow speed at four of the five sites and only 10% at one site located close to a levee-canal combination. Two non-linear, semi-empirical expressions relating flow speeds to the local hydraulic gradient, water depths, and vegetative resistance accounted for 70% of the variance in our measured speed. The data suggest local-scale factors such as channel morphology, vegetation density, and groundwater exchanges must be considered along with landscape position and basin-scale geomorphology when examining the interactions between flow and community characteristics in low-gradient wetlands such as the Everglades.

Biogeochemical effects of iron availability on primary producers in a shallow marine carbonate environment

We completed a synoptic survey of iron, phosphorus, and sulfur concentrations in shallow marine carbonate sediments from south Florida. Total extracted iron concentrations typically were 50 μmol g-1 dry weight (DW) and tended to decrease away from the Florida mainland, whereas total extracted phosphorus concentrations mostly were 10 μmol g-1 DW and tended to decrease from west to east across Florida Bay. Concentrations of reduced sulfur compounds, up to 40 μmol g-1 DW, tended to covary with sediment iron concentrations, suggesting that sulfide mineral formation was iron-limited. An index of iron availability derived from sediment data was negatively correlated with chlorophyll a concentrations in surface waters, demonstrating the close coupling of sediment-water column processes. Eight months after applying a surface layer of iron oxide granules to experimental plots, sediment iron, phosphorus, and sulfur were elevated to a depth of 10 cm relative to control plots. Biomass of the seagrass Thalassia testudinum was not different between control and iron addition plots, but individual shoot growth rates were significantly higher in experimental plots after 8 months. Although the iron content of leaf tissues was significantly higher from iron addition plots, no difference in phosphorus content of T. testudinum leaves was observed. Iron addition altered plant exposure to free sulfide, documented by a significantly higher δ34S of leaf tissue from experimental plots relative to controls. Iron as a buffer to toxic sulfides may promote individual shoot growth, but phosphorus availability to plants still appears to limit production in carbonate sediments.

Groundwater-surface water interactions on tree islands in the Everglades, south Florida

The marked decline in tree island cover across the Everglades over the last century, has been attributed to landscape-scale hydrologic degradation. To preserve and restore Everglades tree islands, a clear understanding of tree island groundwater-surface water interactions is needed, as these interactions strongly influence the chemistry of shallow groundwater and the location and patterns of vegetation in many wetlands. The goal of this work was to define the relationship between groundwater-surface water interactions, plant-water uptake, and the groundwater geochemical condition of tree islands. Groundwater and surface water levels, temperature, and chemistry were monitored on eight constructed and one natural tree island in the Everglades from 2007–2010. Sap flow, diurnal water table fluctuations and stable oxygen isotopes of stem, ground and soil water were used to determine the effect of plant-water uptake on groundwater-surface water interactions. Hydrologic and geochemical modeling was used to further explore the effect of plant-groundwater-surface water interactions on ion concentrations and potential mineral formation.^

Temporal and Spatial Variability of Mesozooplankton in a Shallow Sub-Tropical Bay: Influence of Top-Down Control

Quantifying the relationship between mesozooplankton and water quality parameters identifies the factors that structure the mesozooplankton community and can be used to generate hypotheses regarding the mechanisms that control the mesozooplankton population and potentially the trophic network. To investigate this relationship, mesozooplankton and water quality data were collected in Florida Bay from 1994 to 2004. Three key characteristics were found in the mesozooplankton community structure: (1) there are significant differences between the four sub-regions of Florida Bay; (2) there is a break in May of 1997 with significant differences before and after this date; and (3) there is a positive correlation between mesozooplankton abundance and salinity. The latter two characteristics are closely correlated with predator abundance, indicating the importance of top-down control. Hypersaline periods appear to provide a refuge from predators, allowing mesozooplankton to increase in abundance despite the increased physiological stress.