Hydrological dynamics between a coastal aquifer and the adjacent estuarine system, Biscayne Bay, South Florida

Geochemical and geophysical approaches have been used to investigate the freshwater and saltwater dynamics in the coastal Biscayne Aquifer and Biscayne Bay. Stable isotopes of oxygen and hydrogen, and concentrations of Sr2+ and Ca2+ were combined in two geochemical mixing models to provide estimates of the various freshwater inputs (precipitation, canal water, and groundwater) to Biscayne Bay and the coastal canal system in South Florida. Shallow geophysical electromagnetic and direct current resistivity surveys were used to image the geometry and stratification of the saltwater mixing zone in the near coastal (less than 1km inland) Biscayne Aquifer. The combined stable isotope and trace metal models suggest a ratio of canal input-precipitation-groundwater of 38%–52%–10% in the wet season and 37%–58%–5% in the dry season with an error of 25%, where most (20%) of the error was attributed to the isotope regression model, while the remaining 5% error was attributed to the Sr2+/Ca2+ mixing model. These models suggest rainfall is the dominate source of freshwater to Biscayne Bay. For a bay-wide water budget that includes saltwater and freshwater mixing, fresh groundwater accounts for less than 2% of the total input. A similar Sr 2+/Ca2+ tracer model indicates precipitation is the dominate source in 9 out of 10 canals that discharge into Biscayne Bay. The two-component mixing model converged for 100% of the freshwater canal samples in this study with 63% of the water contributed to the canals coming from precipitation and 37% from groundwater inputs ±4%. There was a seasonal shift from 63% precipitation input in the dry season to 55% precipitation input in the wet season. The three end-member mixing model converged for only 60% of the saline canal samples possibly due to non-conservative behavior of Sr2+ and Ca2+ in saline groundwater discharging into the canal system. Electromagnetic and Direct Current resistivity surveys were successful at locating and estimating the geometry and depth of the freshwater/saltwater interface in the Biscayne Aquifer at two near coastal sites. A saltwater interface that deepened as the survey moved inland was detected with a maximum interpreted depth to the interface of 15 meters, approximately 0.33 km inland from the shoreline. ^

Shifting N and P limitation along a north-south gradient of mangrove estuaries in South Florida

A multivariate statistical analysis was applied to a 10 year, multiparameter data set in an effort to describe the spatial dependence and inherent variation of water quality patterns in the mangrove estuaries of Ten Thousand Islands – Whitewater Bay area. Principal component analysis (PCA) of 16 water quality parameters collected monthly resulted in five groupings, which explained 72.5% of the variance of the original variables. The “Organic” component (PCI) was composed of alkaline phosphatase activity, total organic nitrogen, and total organic carbon; the “Dissolved Inorganic N” component (PCII) contained NO 3 − , NO 2 − , and NH 4 + ; the “Phytoplankton” component (PCIII) was made up of total phosphorus, chlorophyll a, and turbidity; dissolved oxygen and temperature were inversely related (PCIV); and salinity and soluble reactive phosphorus made up PCV. A cluster analysis of the mean and SD of PC scores resulted in the spatial aggregation of the 47 fixed stations into six classes having similar water quality, which we defined as: Mangrove Rivers, Whitewater Bay, Gulf Islands, Coot Bay, Blackwater River, and Inland Waterway. Marked differences in physical, chemical, and biological characteristics among classes were illustrated by this technique. Comparison of medians and variability of parameters among classes allowed large scale generalizations as to underlying differences in water quality in these regions. A strong south to north gradient in estuaries from high N - low P to low N - high P was ascribed to marked differences in landuse, freshwater input, geomorphology, and sedimentary geology along this tract. The ecological significance of this gradient discussed along with potential effects of future restoration plans.

Sediment and Nutrient Deposition Associated with Hurricane Wilma in Mangroves of the Florida Coastal Everglades

The distribution of mangrove biomass and forest structure along Shark River estuary in the Florida Coastal Everglades (FCE) has been correlated with elevated total phosphorus concentration in soils thought to be associated with storm events. The passage of Hurricane Wilma across Shark River estuary in 2005 allowed us to quantify sediment deposition and nutrient inputs in FCE mangrove forests associated with this storm event and to evaluate whether these pulsing events are sufficient to regulate nutrient biogeochemistry in mangrove forests of south Florida. We sampled the spatial pattern of sediment deposits and their chemical properties in mangrove forests along FCE sites in December 2005 and October 2006. The thickness (0.5 to 4.5 cm) of hurricane sediment deposits decreased with distance inland at each site. Bulk density, organic matter content, total nitrogen (N) and phosphorus (P) concentrations, and inorganic and organic P pools of hurricane sediment deposits differed from surface (0–10 cm) mangrove soils at each site. Vertical accretion resulting from this hurricane event was eight to 17 times greater than the annual accretion rate (0.30± 0.03 cm year−1) averaged over the last 50 years. Total P inputs from storm-derived sediments were equivalent to twice the average surface soil nutrient P density (0.19 mg cm−3). In contrast, total N inputs contributed 0.8 times the average soil nutrient N density (2.8 mg cm−3). Allochthonous mineral inputs from Hurricane Wilma represent a significant source of sediment to soil vertical accretion rates and nutrient resources in mangroves of southwestern Everglades. The gradient in total P deposition to mangrove soils from west to east direction across the FCE associated with this storm event is particularly significant to forest development due to the P-limited condition of this carbonate ecosystem. This source of P may be an important adaptation of mangrove forests in the Caribbean region to projected impacts of sea-level rise.

A conceptual model of ecological interactions in the mangrove estuaries of the Florida Everglades

A brackish water ecotone of coastal bays and lakes, mangrove forests, salt marshes, tidal creeks, and upland hammocks separates Florida Bay, Biscayne Bay, and the Gulf of Mexico from the freshwater Everglades. The Everglades mangrove estuaries are characterized by salinity gradients that vary spatially with topography and vary seasonally and inter-annually with rainfall, tide, and freshwater flow from the Everglades. Because of their location at the lower end of the Everglades drainage basin, Everglades mangrove estuaries have been affected by upstream water management practices that have altered the freshwater heads and flows and that affect salinity gradients. Additionally, interannual variation in precipitation patterns, particularly those caused to El Nin˜o events, control freshwater inputs and salinity dynamics in these estuaries. Two major external drivers on this system are water management activities and global climate change. These drivers lead to two major ecosystem stressors: reduced freshwater flow volume and duration, and sea-level rise. Major ecological attributes include mangrove forest production, soil accretion, and resilience; coastal lake submerged aquatic vegetation; resident mangrove fish populations; wood stork (Mycteria americana) and roseate spoonbill (Platelea ajaja) nesting colonies; and estuarine crocodilian populations. Causal linkages between stressors and attributes include coastal transgression, hydroperiods, salinity gradients, and the ‘‘white zone’’ freshwater/estuarine interface. The functional estuary and its ecological attributes, as influenced by sea level and freshwater flow, must be viewed as spatially dynamic, with a possible near-term balancing of transgression but ultimately a long-term continuation of inland movement. Regardless of the spatio-temporal timing of this transgression, a salinity gradient supportive of ecologically functional Everglades mangrove estuaries will be required to maintain the integrity of the South Florida ecosystem.

Chilling damage in a changing climate in coastal landscapes of the subtropical zone: a case study from south Florida

Freeze events significantly influence landscape structure and community composition along subtropical coastlines. This is particularly true in south Florida, where such disturbances have historically contributed to patch diversity within the mangrove forest, and have played a part in limiting its inland transgression. With projected increases in mean global temperatures, such instances are likely to become much less frequent in the region, contributing to a reduction in heterogeneity within the mangrove forest itself. To understand the process more clearly, we explored the dynamics of a Dwarf mangrove forest following two chilling events that produced freeze-like symptoms, i.e., leaf browning, desiccation, and mortality, and interpreted the resulting changes within the context of current winter temperatures and projected future scenarios. Structural effects from a 1996 chilling event were dramatic, with mortality and tissue damage concentrated among individuals comprising the Dwarf forest's low canopy. This disturbance promoted understory plant development and provided an opportunity for Laguncularia racemosa to share dominance with Rhizophora mangle. Mortality due to the less severe 2001 event was greatest in the understory, probably because recovery of the protective canopy following the earlier freeze was still incomplete. Stand dynamics were static over the same period in nearby unimpacted sites. The probability of reaching temperatures as low as those recorded at a nearby meteorological station (≤3 °C) under several warming scenarios was simulated by applying 1° incremental temperature increases to a model developed from a 42-year temperature record. According to the model, the frequency of similar chilling events decreased from once every 1.9 years at present to once every 3.4 and 32.5 years with 1 and 4 °C warming, respectively. The large decrease in the frequency of these events would eliminate an important mechanism that maintains Dwarf forest structure, and promotes compositional diversity.

Coastal groundwater discharge – an additional source of phosphorus for the oligotrophic wetlands of the Everglades

In this manuscript we define a new term we call coastal groundwater discharge (CGD), which is related to submarine groundwater discharge (SGD), but occurs when seawater intrudes inland to force brackish groundwater to discharge to the coastal wetlands. A hydrologic and geochemical investigation of both the groundwater and surface water in the southern Everglades was conducted to investigate the occurrence of CGD associated with seawater intrusion. During the wet season, the surface water chemistry remained fresh. Enhanced chloride, sodium, and calcium concentrations, indicative of brackish groundwater discharge, were observed in the surface water during the dry season. Brackish groundwaters of the southern Everglades contain 1–2.3μM concentrations of total phosphorus (TP). These concentrations exceed the expected values predicted by conservative mixing of local fresh groundwater and intruding seawater, which both have TPμM. The additional source of TP may be from seawater sediments or from the aquifer matrix as a result of water–rock interactions (such as carbonate mineral dissolution and ion exchange reactions) induced by mixing fresh groundwater with intruding seawater. We hypothesize that CGD maybe an additional source of phosphorus (a limiting nutrient) to the coastal wetlands of the southern Everglades.

Seasonal and spatial variation in the stable isotopic composition (δ18O and δD) of precipitation in south Florida

Precipitation data collected from five sites in south Florida indicate a strong seasonal and spatial variation in δ18O and δD, despite the relatively limited geographic coverage and low-lying elevation of each of the collection sites. Based upon the weighted-mean stable isotope values, the sites were classified as coastal Atlantic, inland, and lower Florida Keys. The coastal Atlantic sites had weighted-mean values of δ18O and δD of −2.86‰ and −12.8‰, respectively, and exhibited a seasonal variation with lower δ18O and δD values in the summer wet-season precipitation (δ18O = −3.38‰, δD = −16.5‰) as compared to the winter-time precipitation (δ18O = −1.66‰, δD = −3.2‰). The inland site was characterized as having the highest d-excess value (+13.3‰), signifying a contribution of evaporated Everglades surface water to the local atmospheric moisture. In spite of its lower latitude, the lower Keys site located at Long Key had the lowest weighted-mean stable isotope values (δ18O = −3.64‰, δD = −20.2‰) as well as the lowest d-excess value of (+8.8‰). The lower δD and δ18O values observed at the Long Key site reflect the combined effects of oceanic vapor source, fractionation due to local precipitation, and slower equilibration of the larger raindrops nucleated by a maritime aerosol. Very low δ18O and δD values (δ18O < −6‰, δD < −40‰) were observed just prior to the passage of hurricanes from the Gulf of Mexico as well as during cold fronts from the north-west. These results suggest that an oceanic vapor source region to the west, may be responsible for the extremely low δD and δ18O values observed during some tropical storms and cold fronts.

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.

Paleoenvironmental change in wetlands of the Florida Everglades, southeast USA

In this special issue, we report on efforts to reconstruct paleoclimate/paleolimnology of the Florida Everglades, applying a wide range of techniques including sedimentological, micropaleontological and biogeochemical approaches. The papers included here describe results obtained by studies conducted in Everglades National Park and the greater South Florida Everglades by Florida Coastal Everglades Long Term Ecological Research Program (FCE LTER) collaborators. This multi-investigator project contrasts nutrient dynamics in two inland-to-marine transects aligned along separate drainages in southern Florida that differ in their susceptibility to coastal pressures and in volume of freshwater delivery. This effort focuses on the paleoecological aspects of FCE LTER research that address scales of ecosystem transformations driven by climate variability and change and human activities. The central question addressed by this body of work is “How is the shape of the freshwater-to-marine gradient in the Florida coastal Everglades controlled by changes in climate, freshwater inflow (i.e. through human activities), and disturbance (i.e. sea level rise, hurricanes, fire)?”

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.

Environmental Dynamics of Dissolved Black Carbon in Wetlands

Wetlands are ecosystems commonly characterized by elevated levels of dissolved organic carbon (DOC), and although they cover a surface area less than 2 % worldwide, they are an important carbon source representing an estimated 15 % of global annual DOC flux to the oceans. Because of their unique hydrological characteristics, fire can be an important ecological driver in pulsed wetland systems. Consequently, wetlands may be important sources not only of DOC but also of products derived from biomass burning, such as dissolved black carbon (DBC). However, the biogeochemistry of DBC in wetlands has not been studied in detail. The objective of this study is to determine the environmental dynamics of DBC in different fire-impacted wetlands. An intensive, 2-year spatial and temporal dynamics study of DBC in a coastal wetland, the Everglades (Florida) system, as well as one-time sampling surveys for the other two inland wetlands, Okavango Delta (Botswana) and the Pantanal (Brazil), were reported. Our data reveal that DBC dynamics are strongly coupled with the DOC dynamics regardless of location, season or recent fire history. The statistically significant linear regression between DOC and DBC was applied to estimate DBC fluxes to the coastal zone through two main riverine DOC export routes in the Everglades ecosystem. The presence of significant amounts of DBC in these three fire-impacted ecosystems suggests that sub-tropical wetlands could represent an important continental-ocean carrier of combustion products from biomass burning. The discrimination of DBC molecular structure (i.e. aromaticity) between coastal and terrestrial samples, and between samples collected in wet and dry season, suggests that spatially-significant variation in DBC source strength and/or degree of degradation may also influence DBC dynamics.

Compositional Effects of Sea-Level Rise in a Patchy Landscape: The Dynamics of Tree Islands in the Southeastern Coastal Everglades

The landscape structure of emergent wetlands in undeveloped portions of the southeastern coastal Everglades is comprised of two distinct components: scattered forest fragments, or tree islands, surrounded by a low matrix of marsh or shrub-dominated vegetation. Changes in the matrix, including the inland transgression of salt-tolerant mangroves and the recession of sawgrass marshes, have been attributed to the combination of sea level rise and reductions in fresh water supply. In this study we examined concurrent changes in the composition of the region’s tree islands over a period of almost three decades. No trend in species composition toward more salt-tolerant trees was observed anywhere, but species characteristic of freshwater swamps increased in forests in which fresh water supply was augmented. Tree islands in the coastal Everglades appear to be buffered from some of the short term effects of salt water intrusion, due to their ability to build soils above the surface of the surrounding wetlands, thus maintaining mesophytic conditions. However, the apparent resistance of tree islands to changes associated with sea level rise is likely to be a temporary stage, as continued salt water intrusion will eventually overwhelm the forests’ capacity to maintain fresh water in the rooting zone.

Remote-Sensing Monitoring of Tide Propagation Through Coastal Wetlands

Tide propagation through coastal wetlands is a complex phenomenon affected by vegetation, channels, and tidal conditions. Generally, tidal flow is studied using stage (water level) observations, which provide good temporal resolution, but they are acquired in limited locations. Here, a remote-sensing technique, wetland InSAR (interferometric synthetic aperture radar), is used to detect tidal flow in vegetated coastal environments over broad spatial scales. The technique is applied to data sets acquired by three radar satellites over the western Everglades in south Florida. Interferometric analysis of the data shows that the greatest water-level changes occur along tidal channels, reflecting a high velocity gradient between fast horizontal flow in the channel and the slow flow propagation through the vegetation. The high-resolution observations indicate that the tidal flushing zone extends 2–3 km on both sides of tidal channels and can extend 3–4 km inland from the end of the channel. The InSAR observations can also serve as quantitative constraints for detailed coastal wetland flow models.

Development and evaluation of a geographic information system based method to estimate flooding susceptibility in an area of Broward County, Florida

The objective of this study was to develop a GIS-based multi-class index overlay model to determine areas susceptible to inland flooding during extreme precipitation events in Broward County, Florida. Data layers used in the method include Airborne Laser Terrain Mapper (ALTM) elevation data, excess precipitation depth determined through performing a Soil Conservation Service (SCS) Curve Number (CN) analysis, and the slope of the terrain. The method includes a calibration procedure that uses "weights and scores" criteria obtained from Hurricane Irene (1999) records, a reported 100-year precipitation event, Doppler radar data and documented flooding locations. Results are displayed in maps of Eastern Broward County depicting types of flooding scenarios for a 100-year, 24-hour storm based on the soil saturation conditions. As expected the results of the multi-class index overlay analysis showed that an increase for the potential of inland flooding could be expected when a higher antecedent moisture condition is experienced. The proposed method proves to have some potential as a predictive tool for flooding susceptibility based on a relatively simple approach.

Barrier Island Response to Sea Level Rise in North Carolina

The state of North Carolina is home to some of the most spectacular barrier islands in the world. These features are constantly shifting, impacted by waves, tides, and wind. Studies of the Outer Banks, North Carolina have resulted in varied results, but a detailed analysis of the barrier system as a whole is lacking. Using historic topographic surveys (T-sheets) from the 19th, the positions of various barrier segments were analyzed in relation to modern imagery. Changes in area, width, and center line locations were evaluated over the past 150 years. In total, 74 percent of modern transects have decreased in area. Total reductions in size were 130 km2 for the study period. Mean centerlines as a function of migration showed that 53 percent of segments were demonstrating directional movement away from the ocean. The average movement towards the bay between modern and historic centerlines was 8 meters. Thusly, barrier islands in North Carolina are demonstrating both decreases in total area and directional movement inland in response to sea level rise.