Impacts of hurricanes on surface water flow within a wetland

Between 2001 and 2005, seven category 3 or higher major hurricanes made landfall within the US. The hydrologic impacts of these distinct climatic phenomena frequently occurring in wetland watersheds, however, are not well understood. The focus of this study was to evaluate the impacts of hurricane wind and rainfall conditions on water velocity and water elevations within the study wetland, the Florida Everglades. Specifically water velocity data was measured near two tree islands (Gumbo Limbo (GL) and Satin Leaf (SL)) and wind speed, water elevation, and rainfall were obtained from nearby wind observation stations. During the direct impacts of the hurricanes (Hurricanes Katrina and Wilma), water speed, flow direction, and hydraulic gradients were altered, and the extent of variation was positively related to wind characteristics, with significant alterations in flow direction at depth during Hurricane Wilma due to higher wind speeds. After the direct impacts, the longer lasting effect of hurricanes (time scale of a few days) resulted in altered flow speeds that changed by 50% or less. These longer lasting changes in flow speeds may be due to the redistribution of emergent vegetation.

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.

Smoke on the water: the interplay of fire and water flow on Everglades restoration

Recent research makes clear that much of the Everglade’s flora and fauna have evolved to tolerate or require frequent fires. Nevertheless, restoration of the Everglades has thus far been conceptualized as primarily a water reallocation project. These two forces are directly linked by the influence of water flows on fire fuel moisture content, and are indirectly linked through a series of complex feedback loops. This interaction is made more complex by the alteration and compartmentalization of current water flows and fire regimes, the lack of communication between water and fire management agencies, and the already imperiled state of many local species. It is unlikely, therefore, that restoring water flows will automatically restore the appropriate fire regimes, leaving the prospect of successful restoration in some doubt. The decline of the Cape Sable seaside sparrow, and its potential for recovery, illustrates the complexity of the situation.

Dr. Price measuring surface water flow at SRS-1

http://digitalcommons.fiu.edu/fce_lter_photos/1338/thumbnail.jpg

Graduate student David Lagomasino downloading surface water flow data at SRS-1

http://digitalcommons.fiu.edu/fce_lter_photos/1339/thumbnail.jpg

A tidal creek water budget: Estimation of groundwater discharge and overland flow using hydrologic modeling in the Southern Everglades

Taylor Slough is one of the natural freshwater contributors to Florida Bay through a network of microtidal creeks crossing the Everglades Mangrove Ecotone Region (EMER). The EMER ecological function is critical since it mediates freshwater and nutrient inputs and controls the water quality in Eastern Florida Bay. Furthermore, this region is vulnerable to changing hydrodynamics and nutrient loadings as a result of upstream freshwater management practices proposed by the Comprehensive Everglades Restoration Program (CERP), currently the largest wetland restoration project in the USA. Despite the hydrological importance of Taylor Slough in the water budget of Florida Bay, there are no fine scale (∼1 km2) hydrodynamic models of this system that can be utilized as a tool to evaluate potential changes in water flow, salinity, and water quality. Taylor River is one of the major creeks draining Taylor Slough freshwater into Florida Bay. We performed a water budget analysis for the Taylor River area, based on long-term hydrologic data (1999–2007) and supplemented by hydrodynamic modeling using a MIKE FLOOD (DHI,http://dhigroup.com/) model to evaluate groundwater and overland water discharges. The seasonal hydrologic characteristics are very distinctive (average Taylor River wet vs. dry season outflow was 6 to 1 during 1999–2006) with a pronounced interannual variability of flow. The water budget shows a net dominance of through flow in the tidal mixing zone, while local precipitation and evapotranspiration play only a secondary role, at least in the wet season. During the dry season, the tidal flood reaches the upstream boundary of the study area during approximately 80 days per year on average. The groundwater field measurements indicate a mostly upwards-oriented leakage, which possibly equals the evapotranspiration term. The model results suggest a high importance of groundwater contribution to the water salinity in the EMER. The model performance is satisfactory during the dry season where surface flow in the area is confined to the Taylor River channel. The model also provided guidance on the importance of capturing the overland flow component, which enters the area as sheet flow during the rainy season. Overall, the modeling approach is suitable to reach better understanding of the water budget in the mangrove region. However, more detailed field data is needed to ascertain model predictions by further calibrating overland flow parameters.

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.

Vegetation:environment relationships and water management in Shark Slough, Everglades National Park

The hydrologic regime of Shark Slough, the most extensive long hydroperiod marsh in Everglades National Park, is largely controlled by the location, volume, and timing of water delivered to it through several control structures from Water Conservation Areas north of the Park. Where natural or anthropogenic barriers to water flow are present, water management practices in this highly regulated system may result in an uneven distribution of water in the marsh, which may impact regional vegetation patterns. In this paper, we use data from 569 sampling locations along five cross-Slough transects to examine regional vegetation distribution, and to test and describe the association of marsh vegetation with several hydrologic and edaphic parameters. Analysis of vegetation:environment relationships yielded estimates of both mean and variance in soil depth, as well as annual hydroperiod, mean water depth, and 30-day maximum water depth within each cover type during the 1990’s. We found that rank abundances of the three major marsh cover types (Tall Sawgrass, Sparse Sawgrass, and Spikerush Marsh) were identical in all portions of Shark Slough, but regional trends in the relative abundance of individual communities were present. Analysis also indicated clear and consistent differences in the hydrologic regime of three marsh cover types, with hydroperiod and water depths increasing in the order Tall Sawgrass , Sparse Sawgrass , Spikerush Marsh. In contrast, soil depth decreased in the same order. Locally, these differences were quite subtle; within a management unit of Shark Slough, mean annual values for the two water depth parameters varied less than 15 cm among types, and hydroperiods varied by 65 days or less. More significantly, regional variation in hydrology equaled or exceeded the variation attributable to cover type within a small area. For instance, estimated hydroperiods for Tall Sawgrass in Northern Shark Slough were longer than for Spikerush Marsh in any of the other regions. Although some of this regional variation may reflect a natural gradient within the Slough, a large proportion is the result of compartmentalization due to current water management practices within the marsh.We conclude that hydroperiod or water depth are the most important influences on vegetation within management units, and attribute larger scale differences in vegetation pattern to the interactions among soil development, hydrology and fire regime in this pivotal portion of 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.

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.

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.

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.

Phosphorus enrichment effects on interactions among the ecosystem components in a long-hydroperiod oligotrophic marsh in Everglades National Park

Hydrology and a history of oligotrophy unite the massive landscape comprising freshwater marsh in Everglades National Park. With restoration of water flow to the Everglades, phosphorus (P) enrichment, both from agricultural and domestic sources, may increase nutrient load to the marsh ecosystem. Previous research of P enrichment of Everglades soil, periphyton, and macrophytes revealed each of these ecosystem components responds to increased P loads with increased production and nutrient content. Interactions among these ecosystem components and how P affects the magnitude and direction of interaction are poorly understood and are the focus of my research. Here I present results of a two-year, two-factor experiment of P enrichment and manipulation in Everglades National Park. I quantified biomass, nutrient content, and production for periphyton and macrophyes and found macrophyte removal drives change in nutrient content, biomass, and production of periphyton. Periphyton removal did not appear to control macrophyte dynamics. Soil chemical and physical characteristics were explained primarily by site differences but there was an enrichment effect of soil porewater nitrite + nitrate, nitrite, and soluble reactive phosphorus. Flocculent materials production and depth were significantly affected by macrophyte removal where depth and production were significantly greater with the no macrophyte treatment. The dominant macrophyte of the marsh, Eleocharis cellulosa, increased more in the unenriched marsh than in the enriched marsh. The combination of these findings suggests that dynamics in floc and periphyton are controlled primarily by the presence of periphyton and that this relationship is significantly affected by low-level P enrichment. These results may be valuable in their application to both managers and policy makers who are involved in the Everglades restoration process. ^

Organic matter dynamics in the Florida coastal Everglades: A molecular marker and isotopic approach

Everglades National Park (ENP) is about to undergo the world's largest wetland restoration with the aim of improving the quality, timing and distribution of water flow. The changes in water flow are hypothesized to alter the nutrient fluxes and organic matter (OM) dynamics within ENP, especially in the estuarine areas. This study used a multi-proxy approach of molecular markers and stable δ 13C isotope measurements, to determine the present day OM dynamics in ENP. ^ OM dynamics in wetland soils/sediments have proved to be difficult to understand using traditional geochemical approaches. These are often inadequate to describe the multitude of OM sources (e.g. higher land plant, emergent vegetation, submerged vegetation) to the soils/sediments and the complex diagenetic processes that can alter the OM characteristics. A multi-proxy approach, however, that incorporates both molecular level and bulk parameter information is ideal to comprehend complex OM dynamics in aquatic environments. Therefore, biomass-specific molecular markers or proxies can be useful in tracing the sources and processing of OM. This approach was used to examine the OM dynamics in the two major drainage basins, Shark River Slough and Taylor River Slough, of ENP. Freshwater to marine transects were sampled in both systems for soils/sediments and suspended particulate organic matter (SPOM) to be characterized through bulk OM analyses, lipid biomarker determinations (e.g. sterols, fatty acids, hydrocarbons and triterpenoids) and compound-specific stable carbon isotope (δ 13C) determinations. ^ One key accomplishment of the research was the assessment of a molecular marker proxy (Paq) to distinguish between emergent/higher plant vegetation from submerged vegetation within ENP. This proxy proved to be quite useful at tracing OM inputs to the soils/sediments of ENP. A second key accomplishment was the development of a 3-way model using vegetation specific molecular markers. This novel, descriptive model was successfully applied to the estuarine areas of Taylor and Shark River sloughs, providing clear evidence of mixing of freshwater, estuarine and marine derived OM in these areas. In addition, diagenetic transformations of OM in these estuaries were found to be quite different between Taylor and Shark Rivers, and are likely a result of OM quality and hydrological differences. ^

Hydrologic measurements and implications for tree island formation within Everglades National Park

Tree islands in the Shark River Slough of the Everglades National Park (ENP), in the southern state of Florida in the United States, are part of a wetland system of densely vegetated ridges interspersed within relatively open sloughs. Human alteration of this system has had dramatic negative effects on the landscape of the region and restoration efforts will require adjusting the hydrology of the region to assure the preservation of these important ecologic features. The primary objectives of this study were to document the hydrology in the vicinity of tree islands in ENP by measuring velocities in time and space and by characterizing suspended sediments. The results of such measurements were interpreted with respect to factors that may limit tree island growth. The measurements were conducted in the vicinity of three tree islands known as Black Hammock (BH), Gumbo Limbo (GL), and an unnamed island that was named for this study as Satin Leaf (SL). Acoustical Doppler Velocity (ADV) meters were used for measuring the low velocities of the Everglades water flow. Properties of suspended sediments were characterized through measurements of particle size distribution, turbidity, concentration and particle density. Mean velocities observed at each of the tree islands varied from 0.9 to 1.4 cm/s. Slightly higher mean velocities were observed during the wet season (1.2–1.6 cm/s) versus the dry season (0.8–1.3 cm/s). Maximum velocities of more than 4 cm/s were measured in areas of Cladium jamaicense die-off and at the hardwood hammock (head) of the islands. At the island’s head, water is channelized around obstructions such as tree trunks in relatively rapid flow, which may limit the lateral extent of tree island growth. Channelization is facilitated by shade from the tree canopy, which limits the growth of underwater vegetation thereby minimizing the resistance to flow and limiting sediment deposition. Suspended sediment concentrations were low (0.5–1.5 mg/L) at all study sites and were primarily of organic origin. The mean particle size of the suspended sediments was 3 μm with a distribution that was exponential. Critical velocities needed to cause re-suspension of these particles were estimated to be above the actual velocities observed. Sediment transport within the water column appears to be at a near steady state during the conditions evaluated with low rates of sediment loss balanced by presumably the release of equivalent quantities of particles of organic origin. Existing hydrologic conditions do not appear to transport sufficient suspended sediments to result in the formation of tree islands. Of interest would be to collect hydrologic and sediment transport data during extreme hydrologic events to determine if enough sediment is transported under these conditions to promote sufficient sediment accumulations.