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

A preliminary analysis of the correlation of food-web characteristics with hydrology and nutrient gradients in the southern Everglades

We estimated trophic position and carbon source for three consumers (Florida gar, Lepisosteus platyrhincus; eastern mosquitofish, Gambusia holbrooki; and riverine grass shrimp, Palaemonetes paludosus) from 20 sites representing gradients of productivity and hydrological disturbance in the southern Florida Everglades, U.S.A. We characterized gross primary productivity at each site using light/dark bottle incubation and stem density of emergent vascular plants. We also documented nutrient availability as total phosphorus (TP) in floc and periphyton, and the density of small fishes. Hydrological disturbance was characterized as the time since a site was last dried and the average number of days per year the sites were inundated for the previous 10 years. Food-web attributes were estimated in both the wet and dry seasons by analysis of δ15N (trophic position) and δ13C (food-web carbon source) from 702 samples of aquatic consumers. An index of carbon source was derived from a two-member mixing model with Seminole ramshorn snails (Planorbella duryi) as a basal grazing consumer and scuds (amphipods Hyallela azteca) as a basal detritivore. Snails yielded carbon isotopic values similar to green algae and diatoms, while carbon values of scuds were similar to bulk periphyton and floc; carbon isotopic values of cyanobacteria were enriched in C13compared to all consumers examined. A carbon source similar to scuds dominated at all but one study site, and though the relative contribution of scud-like and snail-like carbon sources was variable, there was no evidence that these contributions were a function of abiotic factors or season. Gar consistently displayed the highest estimated trophic position of the consumers studied, with mosquitofish feeding at a slightly lower level, and grass shrimp feeding at the lowest level. Trophic position was not correlated with any nutrient or productivity parameter, but did increase for grass shrimp and mosquitofish as the time following droughts increased. Trophic position of Florida gar was positively correlated with emergent plant stem density.

Influence of hydrobiogeochemistry on transport of chromium through manganese -containing sediments: Experimental and modeling approaches

Chromium (Cr) is a metal of particular environmental concern, owing to its toxicity and widespread occurrence in groundwater, soil, and soil solution. A combination of hydrological, geochemical, and microbiological processes governs the subsurface migration of Cr. Little effort has been devoted to examining how these biogeochemical reactions combine with hydrologic processes influence Cr migration. This study has focused on the complex problem of predicting the Cr transport in laboratory column experiments. A 1-D reactive transport model was developed and evaluated against data obtained from laboratory column experiments. ^ A series of dynamic laboratory column experiments were conducted under abiotic and biotic conditions. Cr(III) was injected into columns packed with β-MnO 2-coated sand at different initial concentrations, variable flow rates, and at two different pore water pH (3.0 and 4.0). In biotic anaerobic column experiments Cr(VI) along with lactate was injected into columns packed with quartz sand or β-MnO2-coated sand and bacteria, Shewanella alga Simidu (BrY-MT). A mathematical model was developed which included advection-dispersion equations for the movement of Cr(III), Cr(VI), dissolved oxygen, lactate, and biomass. The model included first-order rate laws governing the adsorption of each Cr species and lactate. The equations for transport and adsorption were coupled with nonlinear equations for rate-limited oxidation-reduction reactions along with dual-monod kinetic equations. Kinetic batch experiments were conducted to determine the reduction of Cr(VI) by BrY-MT in three different substrates. Results of the column experiments with Cr(III)-containing influent solutions demonstrate that β-MnO2 effectively catalyzes the oxidation of Cr(III) to Cr(VI). For a given influent concentration and pore water velocity, oxidation rates are higher, and hence effluent concentrations of Cr(VI) are greater, at pH 4 relative to pH 3. Reduction of Cr(VI) by BrY-MT was rapid (within one hour) in columns packed with quartz sand, whereas Cr(VI) reduction by BrY-MT was delayed (57 hours) in presence of β-MnO 2-coated sand. BrY-MT grown in BHIB (brain heart infusion broth) reduced maximum amount of Cr(VI) to Cr(III) followed by TSB (tryptic soy broth) and M9 (minimum media). The comparisons of data and model results from the column experiments show that the depths associated with Cr(III) oxidation and transport within sediments of shallow aquatic systems can strongly influence trends in surface water quality. The results of this study suggests that carefully performed, laboratory column experiments is a useful tool in determining the biotransformation of redox-sensitive metals even in the presence of strong oxidant, like β-MnO2. ^

Hydrogeophysical characterization of anisotropy in the Biscayne Aquifer using geophysical methods

The anisotropy of the Biscayne Aquifer which serves as the source of potable water for Miami-Dade County was investigated by applying geophysical methods. Electrical resistivity imaging, self potential and ground penetration radar techniques were employed in both regional and site specific studies. In the regional study, electrical anisotropy and resistivity variation with depth were investigated with azimuthal square array measurements at 13 sites. The observed coefficient of electrical anisotropy ranged from 1.01 to 1.36. The general direction of measured anisotropy is uniform for most sites and trends W-E or SE-NW irrespective of depth. Measured electrical properties were used to estimate anisotropic component of the secondary porosity and hydraulic anisotropy which ranged from 1 to 11% and 1.18 to 2.83 respectively. 1-D sounding analysis was used to models the variation of formation resistivity with depth. Resistivities decreased from NW (close to the margins of the everglades) to SE on the shores of Biscayne Bay. Porosity calculated from Archie's law, ranged from 18 to 61% with higher values found along the ridge. Higher anisotropy, porosities and hydraulic conductivities were on the Atlantic Coastal Ridge and lower values at low lying areas west of the ridge. The cause of higher anisotropy and porosity is attributed to higher dissolution rates of the oolitic facies of the Miami Formation composing the ridge. The direction of minimum resistivity from this study is similar to the predevelopment groundwater flow direction indicated in published modeling studies. Detailed investigations were carried out to evaluate higher anisotropy at West Perrine Park located on the ridge and Snapper Creek Municipal well field where the anisotropy trend changes with depth. The higher anisotropy is attributed to the presence of solution cavities oriented in the E-SE direction on the ridge. Similarly, the change in hydraulic anisotropy at the well field might be related to solution cavities, the surface canal and groundwater extraction wells.^

Spatial and Temporal Changes in Groundwater Salinity in South Florida

A combination of statistical and interpolation methods and Geographic Information System (GIS) spatial analysis was used to evaluate the spatial and temporal changes in groundwater Cl− concentrations in Collier and Lee Counties (southwestern Florida), and Miami-Dade and Broward Counties (southeastern Florida), since 1985. In southwestern Florida, the average Cl− concentrations in the shallow wells (0–43 m) in Collier and Lee Counties increased from 132 mg L−1 in 1985 to 230 mg L−1 in 2000. The average Cl− concentrations in the deep wells (>43 m) of southwestern Florida increased from 392 mg L−1 in 1985 to 447 mg L−1 in 2000. Results also indicated a positive correlation between the mean sea level and Cl− concentrations and between the mean sea level and groundwater levels for the shallow wells. Concentrations in the Biscayne Aquifer (southeastern Florida) were significantly higher than those of southwestern Florida. The average Cl− concentrations increased from 159 mg L−1 in 1985 to 470 mg L−1 in 2010 for the shallow wells (<33 m) and from 1360 mg L−1 in 1985 to 2050 mg L−1 in 2010 for the deep wells (>33 m). In the Biscayne Aquifer, wells showed a positive or negative correlation between mean sea level and Cl− concentrations according to their location with respect to the saltwater intrusion line. Wells located inland behind canal control structures and west of the saltwater intrusion line showed negative correlation values, whereas wells located east of the saltwater intrusion line showed positive values. Overall, the results indicated that since 1985, there was a potential decline in the available freshwater resources estimated at about 12–17% of the available drinking-quality groundwater of the southeastern study area located in the Biscayne Aquifer.

Taylor Slough Hydrology

Taylor Slough, in Everglades National Park, has experienced an evolution of water management infrastructure since drainage activities arrived in South Florida. This has included the excavation of canals, installation of large capacity pump stations, and a variety of operational strategies focused on resolving the conflict between managing the water level for developed areas while providing water supply for Everglades National Park. This study provides a review of water management practices and the concurrent hydrologic conditions in the Taylor Slough basin and adjacent canal system from 1961 through 2010. Analyses of flow, water level and rainfall data were divided into time periods that correspond to significant changes in structural features and operational plans. In the early 1960s, Taylor Slough was disconnected from the greater Everglades system by the construction of levees upstream. As water supply for Taylor Slough became more urgent, the Slough was connected to the regional water supply system via a network of canals and pump stations to relieve over-drained conditions. The increased water supply and pump capacity succeeded in raising water level and increasing flow and hydroperiod in the marsh.

Diatom-based Models for Inferring Hydrology and Periphyton Abundance in a Subtropical Karstic Wetland: Implications for Ecosystem-Scale Bioassessment

We developed diatom-based prediction models of hydrology and periphyton abundance to inform assessment tools for a hydrologically managed wetland. Because hydrology is an important driver of ecosystem change, hydrologic alterations by restoration efforts could modify biological responses, such as periphyton characteristics. In karstic wetlands, diatoms are particularly important components of mat-forming calcareous periphyton assemblages that both respond and contribute to the structural organization and function of the periphyton matrix. We examined the distribution of diatoms across the Florida Everglades landscape and found hydroperiod and periphyton biovolume were strongly correlated with assemblage composition. We present species optima and tolerances for hydroperiod and periphyton biovolume, for use in interpreting the directionality of change in these important variables. Predictions of these variables were mapped to visualize landscape-scale spatial patterns in a dominant driver of change in this ecosystem (hydroperiod) and an ecosystem-level response metric of hydrologic change (periphyton biovolume). Specific diatom assemblages inhabiting periphyton mats of differing abundance can be used to infer past conditions and inform management decisions based on how assemblages are changing. This study captures diatom responses to wide gradients of hydrology and periphyton characteristics to inform ecosystem-scale bioassessment efforts in a large wetland.

A Review of the Effects of Altered Hydrology and Salinity on Vertebrate Fauna and Their Habitats in Northeastern Florida Bay

Estuarine productivity is highly dependent on the freshwater sources of the estuary. In Florida Bay, Taylor Slough was historically the main source of fresh water. Beginning in about 1960, and culminating with the completion of the South Dade Conveyance System in 1984, water management practice began to change the quantity and distribution of flow from Taylor Slough into Northeastern Florida Bay. These practices altered salinity and hydrologic parameters that had measurable negative impacts on vertebrate fauna and their habitats. Here, I review those impacts from published and unpublished literature and anecdotal observations. Almost all vertebrates covered in this review have shown some form of population decline since 1984; most of the studies implicate declines in food resources as the main stressor on their populations. My conclusion is that the diversion of fresh water resulted in an ecological cascade starting with hydrologic stresses on primary then secondary producers culminating in population declines at the top of the food web.

The Influence of Hydrologic Restoration on Groundwater-Surface Water Interactions in a Karst Wetland, the Everglades (FL, USA)

Efforts to rehydrate and restore surface water flow in karst wetlands can have unintended consequences, as these highly conductive and heterogeneous aquifers create a close connection between groundwater and surface water. Recently, hydrologic restoration efforts in the karstic Taylor Slough portion of the Everglades has changed from point source delivery of canal water (direct restoration), to the use of a series of surface water recharge retention basins (diffuse restoration). To determine the influence of restoration on groundwater-surface water interactions in the Taylor Slough headwaters, a water budget was constructed for 1997–2011 using 70 hydro-meteorological stations. With diffuse restoration, groundwater seepage from the Everglades toward the urban boundary increased, while the downstream delivery of surface water to the main portion of the slough declined. The combined influence of diffuse restoration and climate led to increased intra-annual variability in the volume of groundwater and surface water in storage but supported a more seasonally hydrated wetland compared to the earlier direct tactics. The data further indicated that hydrologic engineering in karst wetland landscapes enhances groundwater-surface water interactions, even those designed for restoration purposes.

The Role of Recharge and Evapotranspiration as Hydraulic Drivers of Ion Concentrations in Shallow Groundwater on Everglades Tree Islands, Florida (USA)

Recently, evapotranspiration has been hypothesized to promote the secondary formation of calcium carbonate year-round on tree islands in the Everglades by influencing groundwater ions concentrations. However, the role of recharge and evapotranspiration as drivers of shallow groundwater ion accumulation has not been investigated. The goal of this study is to develop a hydrologic model that predicts the chloride concentrations of shallow tree island groundwater and to determine the influence of overlying biomass and underlying geologic material on these concentrations. Groundwater and surface water levels and chloride concentrations were monitored on eight constructed tree islands at the Loxahatchee Impoundment Landscape Assessment (LILA) from 2007 to 2010. The tree islands at LILA were constructed predominately of peat, or of peat and limestone, and were planted with saplings of native tree species in 2006 and 2007. The model predicted low shallow groundwater chloride concentrations when inputs of regional groundwater and evapotranspiration-to-recharge rates were elevated, while low evapotranspiration-to-recharge rates resulted in a substantial increase of the chloride concentrations of the shallow groundwater. Modeling results indicated that evapotranspiration typically exceeded recharge on the older tree islands and those with a limestone lithology, which resulted in greater inputs of regional groundwater. A sensitivity analysis indicated the shallow groundwater chloride concentrations were most sensitive to alterations in specific yield during the wet season and hydraulic conductivity in the dry season. In conclusion, the inputs of rainfall, underlying hydrologic properties of tree islands sediments and forest structure may explain the variation in ion concentration seen across Everglades tree islands.

Monitoring of Tree Island Conditions in the Southern Everglades: The Effects of Hurricanes and Hydrology on the Status and Population Dynamics of Sixteen Tropical Hardwood Hammock Tree Islands

In 2005 we began a multi-year intensive monitoring and assessment study of tropical hardwood hammocks within two distinct hydrologic regions in Everglades National Park, under funding from the CERP Monitoring and Assessment Program. In serving as an Annual Report for 2010, this document, reports in detail on the population dynamics and status of tropical hardwood hammocks in Shark Slough and adjacent marl prairies during a 4-year period between 2005 and 2009. 2005-09 was a period that saw a marked drawdown in marsh water levels (July 2006 - July 2008), and an active hurricane season in 2005 with two hurricanes, Hurricane Katrina and Wilma, making landfall over south Florida. Thus much of our focus here is on the responses of these forests to annual variation in marsh water level, and on recovery from disturbance. Most of the data are from 16 rectangular permanent plots of 225-625 m2 , with all trees mapped and tagged, and bi-annual sampling of the tree, sapling, shrub, and herb layer in a nested design. At each visit, canopy photos were taken and later analyzed for determination of interannual variation in leaf area index and canopy openness. Three of the plots were sampled at 2-month intervals, in order to gain a better idea of seasonal dynamics in litterfall and litter turnover. Changes in canopy structure were monitored through a vertical line intercept method.

Application of biomarkers and compound specific stable isotopes for the assessment of hydrology as a driver of organic matter dynamics in the Everglades ecosystem

The Everglades is a sub-tropical coastal wetland characterized among others by its hydrological features and deposits of peat. Formation and preservation of organic matter in soils and sediments in this wetland ecosystem is critical for its sustainability and hydrological processes are important divers in the origin, transport and fate of organic matter. With this in mind, organic matter dynamics in the greater Florida Everglades was studied though various organic geochemistry techniques, especially biomarkers, bulk and compound specific δ13C and δD isotope analysis. The main objectives were focused on how different hydrological regimes in this ecosystem control organic matter dynamics, such as the mobilization of particulate organic matter (POM) in freshwater marshes and estuaries, and how organic geochemistry techniques can be applied to reconstruct Everglades paleo-hydrology. For this purpose organic matter in typical vegetation, floc, surface soils, soil cores, and estuarine suspended particulates were characterized in samples selected along hydrological gradients in the Water Conservation Area 3, Shark River Slough and Taylor Slough. ^ This research focused on three general themes: (1) Assessment of the environmental dynamics and source-specific particulate organic carbon export in a mangrove-dominated estuary. (2) Assessment of the origin, transport and fate of organic matter in freshwater marsh. (3) Assessment of historical changes in hydrological conditions in the Everglades (paleo-hydrology) though biomarkes and compound specific isotope analyses. This study reports the first estimate of particulate organic carbon loss from mangrove ecosystems in the Everglades, provides evidence for particulate organic matter transport with regards to the formation of ridge and slough landscapes in the Everglades, and demonstrates the applicability of the combined biomarker and compound-specific stable isotope approach as a means to generate paleohydrological data in wetlands. The data suggests that: (1) Carbon loss from mangrove estuaries is roughly split 50/50 between dissolved and particulate carbon; (2) hydrological remobilization of particulate organic matter from slough to ridge environments may play an important role in the maintenance of the Everglades freshwater landscape; and (3) Historical changes in hydrology have resulted in significant vegetation shifts from historical slough type vegetation to present ridge type vegetation. ^

Simulating everglades national park hydrology and phosphorus transport under existing and future scenarios using numerical modeling

The Florida Everglades has a long history of anthropogenic changes which have impacted the quantity and quality of water entering the system. Since the construction of Tamiami Trail in the 1920's, overland flow to the Florida Everglades has decreased significantly, impacting ecosystems from the wetlands to the estuary. The MIKE Marsh Model of Everglades National Park (M3ENP) is a numerical model, which simulates Everglades National Park (ENP) hydrology using MIKE SHE/MIKE 11software. This model has been developed to determine the parameters that effect Everglades hydrology and understand the impact of specific flow changes on the hydrology of the system. ^ As part of the effort to return flows to the historical levels, several changes to the existing water management infrastructure have been implemented or are in the design phase. Bridge construction scenarios were programed into the M3ENP model to review the effect of these structural changes and evaluate the potential impacts on water levels and hydroperiods in the receiving Northeast Shark Slough ecosystem. These scenarios have shown critical water level increases in an area which has been in decline due to low water levels. Results from this work may help guide future decisions for restoration designs. ^ Excess phosphorus entering Everglades National Park in South Florida may promote the growth of more phosphorus-opportunistic species and alter the food chain from the bottom up. Two phosphorus transport methods were developed into the M3ENP hydrodynamic model to determine the factors affecting phosphorus transport and the impact of bridge construction on water quality. Results showed that while phosphorus concentrations in surface waters decreased overall, some areas within ENP interior may experience an increase in phosphorus loading which the addition of bridges to Tamiami Trail. Finally, phosphorus data and modeled water level data was used to evaluate the spectral response of Everglades vegetation to increasing phosphorus availability using Landsat imagery.^

Integrated Surface-Ground Water Modeling in Wetlands With Improved Methods to Simulate Vegetative Resistance to Flow

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

Quantitative Comparison of Plant Community Hydrology Using Large-Extent, Long-Term Data

Large-extent vegetation datasets that co-occur with long-term hydrology data provide new ways to develop biologically meaningful hydrologic variables and to determine plant community responses to hydrology. We analyzed the suitability of different hydrological variables to predict vegetation in two water conservation areas (WCAs) in the Florida Everglades, USA, and developed metrics to define realized hydrologic optima and tolerances. Using vegetation data spatially co-located with long-term hydrological records, we evaluated seven variables describing water depth, hydroperiod length, and number of wet/dry events; each variable was tested for 2-, 4- and 10-year intervals for Julian annual averages and environmentally-defined hydrologic intervals. Maximum length and maximum water depth during the wet period calculated for environmentally-defined hydrologic intervals over a 4-year period were the best predictors of vegetation type. Proportional abundance of vegetation types along hydrological gradients indicated that communities had different realized optima and tolerances across WCAs. Although in both WCAs, the trees/shrubs class was on the drier/shallower end of hydrological gradients, while slough communities occupied the wetter/deeper end, the distribution ofCladium, Typha, wet prairie and Salix communities, which were intermediate for most hydrological variables, varied in proportional abundance along hydrologic gradients between WCAs, indicating that realized optima and tolerances are context-dependent.