Diatom-based paleolimnological reconstruction of Quaternary environments in a Florida sinkhole lake

Despite lake sensitivity to climate change, few Florida paleolimnological studies have focused on changes in hydrology. Evidence from Florida vegetation histories raise questions about long-term hydrologic history of Florida lakes, and a 25-year limnological dataset revealed recent climate-driven effects on Lake Annie. The objectives of this research are (1) to use modern diatom assemblages to develop methods for reconstruction of climatic and anthropogenic change (2) to reconstruct both long-term and recent histories of Lake Annie using diatom microfossils. Paleoenvironmental reconstruction models were developed from diatom assemblages of various habitat types from modern lakes. Plankton and sediment assemblages were similar, but epiphytes were distinct, suggesting differences in sediment delivery from different parts of the lakes. Relationships between a variety of physical and chemical data and the diatoms from each habitat type were explored. Total phosphorus (TP), pH, and color were found to be the most relevant variables for reconstruction, with sediment and epiphyte assemblages having the strongest relationships to those variables, six calibration models were constructed from the combination of these habitat types and environmental variables. Reconstructions utilizing the weighted averaging models in this study may be used to directly reveal TP, color, and pH changes from a sediment record, which might be suggestive of hydrologic change as well. These variables were reconstructed from the diatom record from both a long-term (11,000 year) and short-term (100 year) record and showed an interaction between climate-driven and local land-use impacts on Lake Annie. The long-term record begins with Lake Annie as a wetland, then the lake filled to a high stand around 4000 years ago. A period of relative stability after that point was interrupted near the turn of the last century by subtle changes in diatom communities that indicate acidification. Abrupt changes in the diatom communities around 1970 AD suggest recovery from acidification, but concurrent hydrologic change intensified anthropogenic effects on the lake. Diatom evidence for alkalization and phosphorus loading correspond to changes seen in the limnological record.

Development of methods for spatial modeling of sediments quality in Lake Okeechobee, Florida

The major objectives of this dissertation were to develop optimal spatial techniques to model the spatial-temporal changes of the lake sediments and their nutrients from 1988 to 2006, and evaluate the impacts of the hurricanes occurred during 1998–2006. Mud zone reduced about 10.5% from 1988 to 1998, and increased about 6.2% from 1998 to 2006. Mud areas, volumes and weight were calculated using validated Kriging models. From 1988 to 1998, mud thicknesses increased up to 26 cm in the central lake area. The mud area and volume decreased about 13.78% and 10.26%, respectively. From 1998 to 2006, mud depths declined by up to 41 cm in the central lake area, mud volume reduced about 27%. Mud weight increased up to 29.32% from 1988 to 1998, but reduced over 20% from 1998 to 2006. The reduction of mud sediments is likely due to re-suspension and redistribution by waves and currents produced by large storm events, particularly Hurricanes Frances and Jeanne in 2004 and Wilma in 2005. Regression, kriging, geographically weighted regression (GWR) and regression-kriging models have been calibrated and validated for the spatial analysis of the sediments TP and TN of the lake. GWR models provide the most accurate predictions for TP and TN based on model performance and error analysis. TP values declined from an average of 651 to 593 mg/kg from 1998 to 2006, especially in the lake’s western and southern regions. From 1988 to 1998, TP declined in the northern and southern areas, and increased in the central-western part of the lake. The TP weights increased about 37.99%–43.68% from 1988 to 1998 and decreased about 29.72%–34.42% from 1998 to 2006. From 1988 to 1998, TN decreased in most areas, especially in the northern and southern lake regions; western littoral zone had the biggest increase, up to 40,000 mg/kg. From 1998 to 2006, TN declined from an average of 9,363 to 8,926 mg/kg, especially in the central and southern regions. The biggest increases occurred in the northern lake and southern edge areas. TN weights increased about 15%–16.2% from 1988 to 1998, and decreased about 7%–11% from 1998 to 2006.

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.

Multidecadal climate oscillations detected in a transparency record from a subtropical Florida lake

Synchronous interannual variability in water transparency observed in neighboring lakes has been linked to regional precipitation and resultant runoff of dissolved organic material, but many climate forcings oscillate over time scales longer than most limnological records can detect. A strong relationship (R2 5 0.86) between transparency and the previous two years’ rainfall and lake stage in a 25-yr record from a Florida lake enabled us to hindcast transparency from a longer 75-yr record of rainfall and lake stage. Predictions revealed a ,30-yr cycle in transparency linked to the Atlantic Multidecadal Oscillation (AMO). Transparency was greatest (4–8 m) in the cool phase of the AMO (,1962–1993) associated with below-average rainfall in south Florida and lowest (0.1– 3.0 m) during two warm phases (,1932–1961, 1994–present) associated with above-average, but more variable, annual rainfall. Models that predict effects of large-scale hydrologic restoration projects on solute export from South Florida’s expansive wetlands need to account for recent entry into a warm AMO phase, where teleconnections between the AMO phases and runoff are opposite of those shown for the U.S. interior.

Derivation of lake mixing and stratification indices from high-resolution lake buoy data

Lake Analyzer is a numerical code coupled with supporting visualization tools for determining indices of mixing and stratification that are critical to the biogeochemical cycles of lakes and reservoirs. Stability indices, including Lake Number, Wedderburn Number, Schmidt Stability, and thermocline depth are calculated according to established literature definitions and returned to the user in a time series format. The program was created for the analysis of high-frequency data collected from instrumented lake buoys, in support of the emerging field of aquatic sensor network science. Available outputs for the Lake Analyzer program are: water temperature (error-checked and/or down-sampled), wind speed (error-checked and/or down-sampled), metalimnion extent (top and bottom), thermocline depth, friction velocity, Lake Number, Wedderburn Number, Schmidt Stability, mode-1 vertical seiche period, and Brunt-Väisälä buoyancy frequency. Secondary outputs for several of these indices delineate the parent thermocline depth (seasonal thermocline) from the shallower secondary or diurnal thermocline. Lake Analyzer provides a program suite and best practices for the comparison of mixing and stratification indices in lakes across gradients of climate, hydro-physiography, and time, and enables a more detailed understanding of the resulting biogeochemical transformations at different spatial and temporal scales.

Diatom-based paleolimnological reconstruction of regional climate and local land-use change from a protected sinkhole lake in southern Florida, USA

Despite their sensitivity to climate variability, few of the abundant sinkhole lakes of Florida have been the subject of paleolimnological studies to discern patterns of change in aquatic communities and link them to climate drivers. However, deep sinkhole lakes can contain highly resolved paleolimnological records that can be used to track long-term climate variability and its interaction with effects of land-use change. In order to understand how limnological changes were regulated by regional climate variability and further modified by local land-use change in south Florida, we explored diatom assemblage variability over centennial and semi-decadal time scales in an ~11,000-yr and a ~150-yr sediment core extracted from a 21-m deep sinkhole lake, Lake Annie, on the protected property of Archbold Biological Station. We linked variance in diatom assemblage structure to changes in water total phosphorus, color, and pH using diatom-based transfer functions. Reconstructions suggest the sinkhole depression contained a small, acidic, oligotrophic pond ~11000–7000 cal yr BP that gradually deepened to form a humic lake by ~4000 cal yr BP, coinciding with the onset of modern precipitation regimes and the stabilization of sea-level indicated by corresponding palynological records. The lake then contained stable, acidophilous planktonic and benthic algal communities for several thousand years. In the early AD 1900s, that community shifted to one diagnostic of an even lower pH (~5.6), likely resulting from acid precipitation. Further transitions over the past 25 yr reflect recovery from acidification and intensified sensitivity to climate variability caused by enhanced watershed runoff from small drainage ditches dug during the mid-twentieth Century on the surrounding property.

Lake metabolism and the diel oxygen technique: State of the science

Significant improvements have been made in estimating gross primary production (GPP), ecosystem respiration (R), and net ecosystem production (NEP) from diel, “free-water” changes in dissolved oxygen (DO). Here we evaluate some of the assumptions and uncertainties that are still embedded in the technique and provide guidelines on how to estimate reliable metabolic rates from high-frequency sonde data. True whole-system estimates are often not obtained because measurements reflect an unknown zone of influence which varies over space and time. A minimum logging frequency of 30 min was sufficient to capture metabolism at the daily time scale. Higher sampling frequencies capture additional pattern in the DO data, primarily related to physical mixing. Causes behind the often large daily variability are discussed and evaluated for an oligotrophic and a eutrophic lake. Despite a 3-fold higher day-to-day variability in absolute GPP rates in the eutrophic lake, both lakes required at least 3 sonde days per week for GPP estimates to be within 20% of the weekly average. A sensitivity analysis evaluated uncertainties associated with DO measurements, piston velocity (k), and the assumption that daytime R equals nighttime R. In low productivity lakes, uncertainty in DO measurements and piston velocity strongly impacts R but has no effect on GPP or NEP. Lack of accounting for higher R during the day underestimates R and GPP but has no effect on NEP. We finally provide suggestions for future research to improve the technique.

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.