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.

Arctic ecosystem responses to changes in water availability and warming: Short and long-term responses

Arctic soils store close to 14% of the global soil carbon. Most of arctic carbon is stored below ground in the permafrost. With climate warming the decomposition of the soil carbon could represent a significant positive feedback to global greenhouse warming. Recent evidence has shown that the temperature of the Arctic is already increasing, and this change is associated mostly with anthropogenic activities. Warmer soils will contribute to permafrost degradation and accelerate organic matter decay and thus increase the flux of carbon dioxide and methane into the atmosphere. Temperature and water availability are also important drivers of ecosystem performance, but effects can be complex and in opposition. Temperature and moisture changes can affect ecosystem respiration (ER) and gross primary productivity (GPP) independently; an increase in the net ecosystem exchange can be a result of either a decrease in ER or an increase in GPP. Therefore, understanding the effects of changes in ecosystem water and temperature on the carbon flux components becomes key to predicting the responses of the Arctic to climate change. The overall goal of this work was to determine the response of arctic systems to simulated climate change scenarios with simultaneous changes in temperature and moisture. A temperature and hydrological manipulation in a naturally-drained lakebed was used to assess the short-term effect of changes in water and temperature on the carbon cycle. Also, as part of International Tundra Experiment Network (ITEX), I determined the long-term effect of warming on the carbon cycle in a natural hydrological gradient established in the mid 90's. I found that the carbon balance is highly sensitive to short-term changes in water table and warming. However, over longer time periods, hydrological and temperature changed soil biophysical properties, nutrient cycles, and other ecosystem structural and functional components that down regulated GPP and ER, especially in wet areas.

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

Investigation of spatial patterns of ground-water exchange with lakes, using a three-dimensional numerical model

Hydrogeologic variables controlling groundwater exchange with inflow and flow-through lakes were simulated using a three-dimensional numerical model (MODFLOW) to investigate and quantify spatial patterns of lake bed seepage and hydraulic head distributions in the porous medium surrounding the lakes. Also, the total annual inflow and outflow were calculated as a percentage of lake volume for flow-through lake simulations. The general exponential decline of seepage rates with distance offshore was best demonstrated at lower anisotropy ratio (i.e., Kh/Kv = 1, 10), with increasing deviation from the exponential pattern as anisotropy was increased to 100 and 1000. 2-D vertical section models constructed for comparison with 3-D models showed that groundwater heads and seepages were higher in 3-D simulations. Addition of low conductivity lake sediments decreased seepage rates nearshore and increased seepage rates offshore in inflow lakes, and increased the area of groundwater inseepage on the beds of flow-through lakes. Introduction of heterogeneity into the medium decreased the water table and seepage ratesnearshore, and increased seepage rates offshore in inflow lakes. A laterally restricted aquifer located at the downgradient side of the flow-through lake increased the area of outseepage. Recharge rate, lake depth and lake bed slope had relatively little effect on the spatial patterns of seepage rates and groundwater exchange with lakes.

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

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

Sea level rise and South Florida coastal forests

Coastal ecosystems lie at the forefront of sea level rise. We posit that before the onset of actual inundation, sea level rise will influence the species composition of coastal hardwood hammocks and buttonwood (Conocarpus erectus L.) forests of the Everglades National Park based on tolerance to drought and salinity. Precipitation is the major water source in coastal hammocks and is stored in the soil vadose zone, but vadose water will diminish with the rising water table as a consequence of sea level rise, thereby subjecting plants to salt water stress. A model is used to demonstrate that the constraining effect of salinity on transpiration limits the distribution of freshwater-dependent communities. Field data collected in hardwood hammocks and coastal buttonwood forests over 11 years show that halophytes have replaced glycophytes. We establish that sea level rise threatens 21 rare coastal species in Everglades National Park and estimate the relative risk to each species using basic life history and population traits. We review salinity conditions in the estuarine region over 1999–2009 and associate wide variability in the extent of the annual seawater intrusion to variation in freshwater inflows and precipitation. We also examine species composition in coastal and inland hammocks in connection with distance from the coast, depth to water table, and groundwater salinity. Though this study focuses on coastal forests and rare species of South Florida, it has implications for coastal forests threatened by saltwater intrusion across the globe.

Comparison of Belowground Biomass in C3- and C4-Dominated Mixed Communities in a Chesapeake Bay Brackish Marsh

Belowground biomass is a critical factor regulating ecosystem functions of coastal marshes, including soil organic matter (SOM) accumulation and the ability of these systems to keep pace with sea-level rise. Nevertheless, belowground biomass responses to environmental and vegetation changes have been given little emphasis marsh studies. Here we present a method using stable carbon isotopes and color to identify root and rhizomes of Schoenoplectus americanus (Pers.) Volk. ex Schinz and R. Keller (C3) and Spartina patens (Ait.) Muhl. (C4) occurring in C3− and C4-dominated communities in a Chesapeake Bay brackish marsh. The functional significance of the biomass classes we identified is underscored by differences in their chemistry, depth profiles, and variation in biomass and profiles relative to abiotic and biotic factors. C3 rhizomes had the lowest concentrations of cellulose (29.19%) and lignin (14.43%) and the lowest C:N (46.97) and lignin:N (0.16) ratios. We distinguished two types of C3 roots, and of these, the dark red C3 roots had anomalously high C:N (195.35) and lignin:N (1.14) ratios, compared with other root and rhizome classes examined here and with previously published values. The C4-dominated community had significantly greater belowground biomass (4119.1 g m−2) than the C3-dominated community (3256.9 g m−2), due to greater total root biomass and a 3.6-fold higher C3-root:rhizome ratio in the C4-dominated community. C3 rhizomes were distributed significantly shallower in the C4-dominated community, while C3 roots were significantly deeper. Variability in C3 rhizome depth distributions was explained primarily by C4 biomass, and C3 roots were explained primarily by water table height. Our results suggest that belowground biomass in this system is sensitive to slight variations in water table height (across an 8 cm range), and that the reduced overlap between C3 and C4 root profiles in the C4-dominated community may account for the greater total root biomass observed in that community. Given that future elevated atmospheric CO2 and accelerated sea-level rise are likely to increase C3 abundance in Atlantic and Gulf coast marshes, investigations that quantify how patterns of C3 and C4 belowground biomass respond to environmental and biological factors stand to improve our understanding of ecosystem-wide impacts of global changes on coastal wetlands.

Hydrologic processes on tree islands in the Everglades (Florida, USA): tracking the effects of tree establishment and growth

The hydrodynamics of tree islands during the growth of newly planted trees has been found to be influenced by both vegetation biomass and geologic conditions. From July 2007 through June 2009, groundwater and surface-water levels were monitored on eight recently planted tree islands at the Loxahatchee Impoundment Landscape Assessment (LILA) facility in Boynton Beach, Florida, USA. Over the 2-year study, stand development coincided with the development of a water-table depression in the center of each of the islands that was bounded by a hydraulic divide along the edges. The water-table depression was greater in islands composed of limestone as compared to those composed of peat. The findings of this study suggest that groundwater evapotranspiration by trees on tree islands creates complex hydrologic interactions between the shallow groundwater in tree islands and the surrounding surface water and groundwater bodies.

Tree Island Response to Fire and Flooding in the Short-Hydroperiod Marl Prairie Grasslands of the Florida Everglades

Within the marl prairie grasslands of the Florida Everglades, USA, the combined effects of fire and flooding usually lead to very significant changes in tree island structure and composition. Depending on fire severity and post-fire hydroperiod, these effects vary spatially and temporally throughout the landscape, creating a patchy post-fire mosaic of tree islands with different successional states. Through the use of the Normalized Difference Vegetation Index (NDVI) and three predictor variables (marsh water table elevation at the time of fire, post-fire hydroperiod, and tree island size), along with logistic regression analysis, we examined the probability of tree island burning and recovering following the Mustang Corner Fire (May to June 2008) in Everglades National Park. Our data show that hydrologic conditions during and after fire, which are under varying degrees of management control, can lead to tree island contraction or loss. More specifically, the elevation of the marsh water table at the time of the fire appears to be the most important parameter determining the severity of fire in marl prairie tree islands. Furthermore, in the post-fire recovery phase, both tree island size and hydroperiod during the first year after the fire played important roles in determining the probability of tree island recovery, contraction, or loss.

The Influence of Vegetation on the Hydrodynamics and Geomorphology of a Tree Island in Everglades National Park (Florida, United States)

Transpiration-driven nutrient accumulation has been identified as a potential mechanism governing the creation and maintenance of wetland vegetation patterning. This process may contribute to the formation of nutrient-rich tree islands within the expansive oligotrophic marshes of the Everglades (Florida, United States). This study presents hydrogeochemical data indicating that tree root water uptake is a primary driver of groundwater ion accumulation across one of these islands. Sap flow, soil moisture, water level, water chemistry, and rainfall were measured to identify the relationships between climate, transpiration, and groundwater uptake by phreatophytes and to examine the effect this uptake has on groundwater chemistry and mineral formation in three woody plant communities of differing elevations. During the dry season, trees relied more on groundwater for transpiration, which led to a depressed water table and the advective movement of groundwater and dissolved ions, including phosphorus, from the surrounding marsh towards the centre of the island. Ion exclusion during root water uptake led to elevated concentrations of all major dissolved ions in the tree island groundwater compared with the adjacent marsh. Groundwater was predominately supersaturated with respect to aragonite and calcite in the lower-elevation woody communities, indicating the potential for soil formation. Elevated groundwater phosphorous concentrations detected in the highest-elevation woody community were associated with the leaching of inorganic sediments (i.e. hydroxyapatite) in the vadose zone. Understanding the complex feedback mechanisms regulating plant/groundwater/surface water interactions, nutrient dynamics, and potential soil formation is necessary to manage and restore patterned wetlands such as the Everglades.

Seasonal activity and road mortality of the snakes of the Pa-Hay-Okee wetlands of Everglades National Park

The current study describes the composition and activity of the snake community of the Pa-hay-okee wetlands of Everglades National Park. The study was conducted from January 1987 to January 1989. Sixteen species were observed, with Thamnophis sauritus, Thamnophis sirtalis, Nerodia fasciata pictiventris, and Agkistrodon piscivorus representing 90.2% of the total sample. The seasonal distribution and activity of the snakes were closely related to fluctuations in the water table. Most activity occurred in the winter months as snakes migrated west following the drying water edge of Shark River Slough. Seventy percent of all snakes observed during this study were either injured or dead on the road. Over 50% of annual mortality occurred during migration. The impact that road mortality is having on the local snake community cannot be ignored. Management options are provided to minimize loss. A comparison is made to the snake community of the Long Pine Key Region of Everglades National Park.