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.

Effect of Water Management on Interannual Variation in Bulk Soil Properties from the Eastern Coastal Everglades

We examined interannual variation in soil properties from wetlands occurring in adjacent drainage basins from the southeastern Everglades. Triplicate 10-cm soil cores were collected, homogenized, and analyzed during the wet season 2006–2010 from five freshwater sawgrass wetland marshes and three estuarine mangrove forests. Soil bulk density from the Taylor Slough basin ranged from 0.15 gm-cm−3 to 0.5 gm-cm−3, was higher than from the Panhandle basin every year, and generally increased throughout the study period. Organic matter as a percent loss on ignition ranged from 7 % to 12 % from freshwater marshes and from 13 % to 56 % from estuarine mangroves. Extractable iron in soils was similar among drainage basins and wetland types, typically ranging from 0.6 to 2.0 g Fe kg−1. In contrast, inorganic sulfur was on average over four times higher from estuarine soils relative to freshwater, and was positively correlated with soil organic matter. Finally total soil phosphorus (P) was lower in freshwater soils relative to estuarine soils (84 ± 5 versus 326 ± 32 mg P kg−1). Total P from the freshwater marshes in the Panhandle basin rose throughout the study period from 54.7 ± 8.4 to 107 ± 17 mg P kg−1, a possible outcome of differences in water management between drainage basins.

Areal and vertical distribution of total soil mercury and total phosphorusin the southern half of water conservation area 3-A, Everglades, Southern Florida

Total soil-mercury and phosphorus concentrations were determined in 64 sites in the southern half of Water Conservation Area 3A, an area of approximately 500 km2 . Surface soil-Hg concentrations ranged from 117 to 300 ng-g-1;total phosphorus concentrations range from 350 to 850 pg~g-1. No consistent north-south or east-west trends are found in the mercury or phosphorus surface concentrations when they are normalized to soil bulk density. Nine sites were used for the determination of the vertical distribution of soilmercury. Vertical profiles of soil-Hg revealed decreasing concentrations with depth and correlated well with phosphorus in soil profiles. Mercury concentrations in soil profiles may be interpreted as an increase in the rate of deposition of mercury in the region in recent decades and/or as postdepositionalmobilization of mercury to surface layers.