Geology topics
Research about Florida
Source-linked reports with geographic coverage including Florida.
Flooding-induced failure of an invasive Burmese Python nest in southern Florida
It is important to understand the factors affecting the reproductive success of an invasive species to estimate population size and develop management plans. There remains much we do not understand about the reproductive biology of invasive Burmese Pythons in both their native and invasive range. Oviposition site selection is an important factor in determining reproductive success as nesting is a vulnerable period in the life of reptiles. Flooding can be particularly influential for nesting outcomes and success of developing embryos in habitats that experience periods of heavy rainfall or seasonal flooding. It is not clear how seasonal flooding may impact oviposition site selection by females or hatchling survival of Burmese Pythons in Florida. On 10 May 2023, a radio-transmitted female python oviposited in a hollow log and during incubation all eggs were completely submerged for variable amounts of time. All eggs were collected (n = 78) and dissected to confirm fertilization and assess the stage of embryo development. Two eggs were found to be infertile, while the remaining fertile eggs contained embryos that died at various stages of development. The observed failure of this python nest was a direct result of oviposition site selection by the female. The frequency at which females lay eggs in suboptimal locations in southern Florida is currently unknown. This maladaptive behavior would suggest Burmese Pythons have fitness levels below theoretical optima.
Stony coral tissue loss disease (SCTLD) case definition for wildlife
Diagnostic laboratories receive carcasses and samples for diagnostic evaluation and pathogen/toxin detection. Case definitions bring clarity and consistency to the evaluation process. Their use within and between organizations allows more uniform reporting of diseases and etiologic agents. The intent of a case definition is to provide scientifically based criteria for determining (a) if an individual carcass has a specific disease and degree of confidence in that diagnosis and (b) if there is evidence of a pathogen or toxin in a carcass or sample (for example, swab, tissue sample, skin scraping, blood/serum sample, environmental sample, or other). This case definition is specific to Stony Coral Tissue Loss Disease (SCTLC) and applies to several species of scleractinian corals across seven families: Faviidae, Meandrinidae, Merulinidae, Montastraeidae, Astrocoeniidae, Scleractinia, and Siderastreidae. Other species presumed to be susceptible are in Families Agaridiidae, Faviidae, and Pocilloporidae.
Simulation modeling to assess line transect distance sampling under a range of translocation scenarios
The accuracy of posttranslocation population monitoring methods is critical to assessing long-term success in translocation programs. Translocation can produce unique challenges to monitoring efforts; therefore, it is important to understand the flexibility and robustness of commonly used monitoring methods. In Florida, USA, thousands of gopher tortoises Gopherus polyphemus have been, and continue to be, translocated from development sites to permitted recipient sites. These recipient sites create a broad range of potential monitoring scenarios due to variability in soft-release strategies, habitat conditions, and population demographics. Line transect distance sampling is an effective method for monitoring natural tortoise populations, but it is currently untested for translocated populations. We therefore produced 3,024 individual-based, spatially explicit scenarios of translocated tortoise populations that differed in recipient site and tortoise population properties, based on real-world examples, literature review, and expert opinion. We virtually sampled simulated tortoise populations by using line transect distance sampling methods and built a Bayesian hierarchical model to estimate the population density for each simulation, which incorporated individual-level covariates (i.e., burrow width and burrow occupancy). Line transect distance sampling was largely appropriate for the conditions that typify gopher tortoise recipient sites, particularly when detection probability on the transect lines was greater than or equal to 0.85. Designing the layout of transects relative to the orientation of soft-release pens, to avoid possible sampling biases that lead to extreme outliers in estimates of tortoise densities, resulted in more accurate population estimates. We also suggest that use of individual-level covariates, applied using a Bayesian framework as demonstrated in our study, may improve the applicability of line transect distance sampling surveys in a variety of contexts and that simulation can be a powerful tool for assessing survey design in complex sampling situations.
Linking temperature sensitivity of mangrove communities, populations and individuals across a tropical-temperate transitional zone
Climate change is reshaping coastal wetlands worldwide, driving ecosystem shifts like mangrove poleward expansion into saltmarshes in tropical-temperate transitional zones. Though warming is recognized as the primary driver, a lack of detailed field studies limits our ability to predict mangrove responses to rapid climate warming. Here, we characterized how mangroves vary across a temperature gradient at 18 sites along Florida's Gulf of Mexico coast (USA). We used minimum air temperature ( T min ) derived from daily data from 1989 to 2021 as the independent variable and applied plot-based and synoptic approaches to quantify species-specific mangrove variation at community, population, and individual levels. We then used these results to spatially project future mangrove ecosystem properties under multiple warming scenarios. Across the T min gradient from −10.8 to −1.4°C, mangrove canopy height and coverage ranged from 0.4 to 11.5 m and 15% to 98%, respectively, with both exhibiting sigmoidal increases with T min . Estimated mangrove aboveground biomass ranged from 0 to 496.7 Mg/ha and showed a positive linear relationship with T min due both to the tall tree stratum's increased biomass per tree and higher abundance. While the population abundance and coverage of Rhizophora mangle and Laguncularia racemosa had positive linear relationships with T min , Avicennia germinans exhibited a significant quadratic relationship, reflecting the higher freeze tolerance of this species. Such tolerance may stem from A. germinans ' higher morphological plasticity observed at the individual level, adapting to cold stress by exhibiting a more shrub-like architecture at colder sites. Based on these field-derived quantitative relationships, we projected substantial increases in mangrove coverage and canopy height near current range limits, with tall A. germinans dominating in the north and R. mangle dominating the centre and south of the study region. Synthesis. To better predict the ecological consequences in coastal wetlands under future climate change, it is essential to understand how mangroves respond to winter temperature regimes across a temperature gradient. Collectively, these cross-level and species-specific results advance our understanding of mangrove temperature sensitivity and provide information about the future of coastal wetland structure and function in response to a changing climate.
Modern coral range expansion off southeast Florida falls short of Late Holocene baseline
As thermal stress and disease outbreaks decimate coral reefs throughout the tropics, there is growing evidence that higher latitude marine environments may provide crucial refuges for many at-risk, temperature-sensitive coral species. However, our understanding of how coral populations expand into new areas and sustain themselves over time is constrained by the limited scope of modern observations. Here, we provide geological insights into coral range expansions by reconstructing the composition of a Late Holocene-aged subfossil coral death assemblage on the southeast Florida reef tract and comparing it to modern reefs throughout the region. Our findings show that the Late Holocene coral assemblages were dominated by now critically endangered Acropora species between ~3500 and 1800 years before present, mirroring classic zonation patterns characteristic of healthy pre-1970s Caribbean reefs. In contrast, the modern reefs off southeast Florida are becoming increasingly dominated by stress-tolerant species like Porites astreoides and Siderastrea siderea despite modest expansions of Acropora cervicornis over the past several decades. Our results suggest that ongoing anthropogenic stressors, not present during the Late Holocene, are likely limiting the ability of modern higher latitude reefs in Florida to function as long-term climate refugia.
Insight into sources of benzene, TCE, and PFOA/PFAS in groundwater at Naval Air Station Whiting Field, Florida, through numerical particle-tracking simulations
Past waste-disposal activities at Naval Air Station Whiting Field (NASWF) have led to elevated concentrations of contaminants in the underlying sand and gravel aquifer. Contaminants include two of the most commonly detected chemicals in groundwater in many countries (benzene and trichloroethylene (TCE)) and the “forever chemicals” per- and poly-fluoroalkyl substances (PFAS) such as perfluorooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS). A MODFLOW model (the Whiting Field Groundwater Model (WFGM)) was previously developed for NASWF and the surrounding area to simulate groundwater flow. To obtain insight into groundwater flow pathways for the identification of potential source areas, the MODPATH particle-tracking application was applied to the WFGM for three public supply wells and three monitoring wells at NASWF. The travel time to recharge areas was estimated using concentrations of the groundwater age-dating solutes tritium (as helium ingrowth) and chlorofluorocarbons detected in the monitoring wells. Simulated travel times agree with the groundwater ages and indicate that the calibrated WFGM reasonably represents groundwater flow velocities and pathways. The MODPATH simulations confirm suspected on-base source areas to explain chemical detection in the monitoring wells. In contrast, the particle-tracking simulations indicate that potential source areas to the public supply wells include both on- and off-base sources. This is important because PFAS chemicals can have multiple sources, including land application of sludge-based fertilizers. This approach that combines groundwater age dating with particle-tracking simulations can be applied at similar sites characterized by benzene-, TCE-, and PFAS-contaminated groundwater.
20th century warming in the lower Florida Keys was dominated by increasing winter temperatures
Long-lived Atlantic coral species like Orbicella faveolata are important archives of oceanographic change in shallow, marine environments like the Florida Keys. Not only can coral-based records extend for multiple centuries beyond the limits of the instrumental record, but they can also provide a more accurate representation of in situ conditions than gridded interpolated sea-surface temperature (SST) products for nearshore reef environments. We use the coral Sr/Ca paleothermometer to produce a 150-year (1830–1980 C.E.) monthly SST reconstruction from an O . faveolata colony collected in the Marquesas Keys, FL, USA. An important feature of our record is a significant 20th-century warming trend in winter SSTs. We hypothesize that the winter warming trend was driven partially by a decrease in upwelling associated cyclonic eddies spinning off the Florida Current. A long-term weakening of winter Florida Current transport over the 20th century could be responsible for decreased cyclonic eddy formation in the Florida Straits. Another feature of the record is pronounced multidecadal fluctuations of mean annual warming and cooling in the record, which correspond to Atlantic Multidecadal Variability (AMV), with the AMV lagging behind western Florida Keys temperatures by 5–11 years. Strong coherence between coral-based SST reconstructions in the western Florida Keys with broader scale Atlantic oceanographic trends over the past century suggests a common driver of regional SST variability.
The 3D Elevation Program—Supporting Florida's economy
Introduction Florida has the longest coastline of any State in the contiguous United States, and its coastal resources are one of the main drivers of its economic growth. High-quality elevation data are beneficial for use in emergency management, especially for hurricane response, recovery, and mitigation, as well as for coastal zone management, flood risk management, infrastructure planning, agriculture, forestry, and natural resources management. Having regional or statewide elevation data coverage that was collected at about the same time allows for improved results from in-depth modeling and more meaningful analysis to support the State’s Chief Science Officer, Geographic Information Officer, State agencies, water management districts, and local governments that will use these data for decision making. Critical applications that meet the State’s management needs depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features.
Nutrient and carbonate chemistry patterns associated with Karenia brevis blooms in three West Florida Shelf estuaries 2020-2023
Ocean acidification (OA) driven by eutrophication, riverine discharge, and other threats from local population growth that affect the inorganic carbonate system is already affecting the eastern Gulf of Mexico. Long-term declines in pH of ~ -0.001 pH units yr -1 have been observed in many southwest Florida estuaries over the past few decades. Coastal and estuarine waters of southwest Florida experience high biomass harmful algal blooms (HABs) of the dinoflagellate Karenia brevis nearly every year; and these blooms have the potential to impact and be impacted by seasonal to interannual patterns of carbonate chemistry. Sampling was conducted seasonally along three estuarine transects (Tampa Bay, Charlotte Harbor, Caloosahatchee River) between May 2020 and May 2023 to obtain baseline measurements of carbonate chemistry prior to, during, and following K. brevis blooms. Conductivity, temperature and depth data and discrete water samples for K. brevis cell abundance, nutrients, and carbonate chemistry (total alkalinity, dissolved inorganic carbonate (DIC), p CO 2, and pH T were evaluated to identify seasonal patterns and linkages among carbonate system variables, nutrients, and K. brevis blooms. Karenia brevis blooms were observed during six samplings, and highest p CO 2 and lowest pH T was observed either during or after blooms in all three estuaries. Highest average pH and lowest p CO 2 were observed in Tampa Bay. In all three estuaries, average DIC and pH T were higher and p CO 2 was lower during dry seasons than wet seasons. There was strong influence of net community calcification (NCC) and net community production (NCP) on the carbonate system; and NCC : NCP ratios in Tampa Bay, Charlotte Harbor, and the Caloosahatchee River were 0.83, 0.93, and 1.02, respectively. Linear relationships between salinity and dissolved ammonium, phosphate, and nitrate indicate strong influence of freshwater inflow from river input and discharge events on nutrient concentrations. This study is a first step towards connecting observations of high biomass blooms like those caused by K. brevis and alterations of carbonate chemistry in Southwest Florida. Our study demonstrates the need for integrated monitoring to improve understanding of interactions among the carbonate system, HABs, water quality, and acidification over local to regional spatial scales and event to decadal time scales.
A spatiotemporal optimization engine for prescribed burning in the Southeast US
Many ecosystems in the Southeast US are dependent upon frequent low-intensity surface fires to sustain native biodiversity, ecosystem services, and endangered species populations. Today, landscape-scale prescribed fire is required to manage these systems for conservation objectives and to mitigate wildland fire risk. Successful application of prescribed fire in this region requires careful planning and assessment of the risks and tradeoffs involved when deciding whether or not to conduct a burn. Many of these risks are closely tied to ambient environmental conditions and are reflected in sets of ‘prescription’ parameters that define safe and effective operating conditions to meet objectives or regulatory requirements. To facilitate effective decision making and acknowledge growing uncertainties related to climate change effects on wildland fire operations, we developed a spatiotemporal optimization engine to identify near-term optimal burning opportunities for prescribed fire implementation. By mining historical 3-day numerical weather forecasts and observation-based weather data for 2015–2021, we have developed a Bayesian hierarchical model for forecast verification that provides calibrated daily weather forecasts and joint uncertainty estimates on meteorological variables of interest, with the latter serving as a measure of risk associated with prescribed fire activities. Burn allocation decisions are then optimized by considering this risk jointly with the utility of burning a particular habitat parcel. The initial iteration of the optimization engine is demonstrated through a case study of short-term meteorological conditions for the Eglin Air Force Base, located in Florida, USA. Results indicate agreement between the optimization engine and the observed past decision-making, with the largest divergences likely arising primarily from differences between utility functions presumed important and used to develop the optimization engine versus the true utility functions driving management behavior in practice.
In vitro impacts of glyphosate on manatee lymphocytes
Exposure to contaminants, such as the herbicide glyphosate, can suppress protective immune functions. Glyphosate is the herbicide most used worldwide and has been found in the plasma of more than 50 % of the Florida manatees ( Trichechus manatus latirostris ) and all-year-round in their aquatic environment. Our objectives were to analyze the consequences of glyphosate exposure on their immune responses via T-lymphocyte proliferation assays and transcriptomics. We isolated peripheral blood mononuclear cells (mainly lymphocytes) of free-ranging manatees and performed T-cell proliferation assays. We used transcriptomics to understand the consequences of glyphosate in vitro exposure. The three doses chosen ranged from environmentally relevant concentrations at 10 to 10,000 µg.L −1 that is considered an environmental contamination scenario. Glyphosate caused a dose-dependent reduction in T-lymphocyte proliferation, with a significant mean reduction of 27.3 % at 10,000 µg.L −1 and up to 51.5 % in some individuals. Additionally, T-lymphocyte proliferation was significantly reduced in mid-winter compared to early winter. Transcriptomic analysis of peripheral blood mononuclear cells indicated that all doses of glyphosate (10, 1,000, and 10,000 µg.L −1 ) resulted in up-regulation of genes related to acute phase inflammation and inhibition of the T-lymphocyte proliferation pathway. Exposure to this contaminant along with other environmental stressors, such as extreme winters and red tide, might further affect the adaptive immune response of this threatened species.
Formation of vertical columnar seismic structures and seafloor depressions by groundwater discharge in the drowned Miami Terrace platform and overlying deep-water carbonates, southeastern Florida
The presence of vertical cross-formational fluid migration passageways within sedimentary basins can profoundly impact aquifer and reservoir fluid-flow and their identification is fundamental to informing management of subsurface fluid resources (groundwater, oil, gas). In an onshore and offshore southeastern part of Florida, 2D/3D seismic-reflection and bathymetry data document ∼153 vertical columnar structures composed of reflection disruptions up to 790 m in the height and averaging 360 m in diameter, and ∼219 subcircular to circular seafloor depressions up to 1334 m wide. Our study focuses on these features found within the offshore shallow-marine carbonate Miami Terrace platform, which drowned approximately at the end of the middle Miocene, and within overlying Plio-Quaternary deep-water carbonate slope and drift deposits. Most columnar structures are rooted in stratiform aquifers of the Miami Terrace platform and associated with faults or fault intersections produced by Eocene and circa late Miocene tectonics. The columns commonly terminate within the platform or as subcircular depressions along an amalgamated karstic and drowning unconformity at the platform top. The columns typically stretch upwards from a zone of deep karst cavity collapse through the Miami Terrace platform with upward decreasing sag on internal reflections. Following drowning and Plio-Quaternary partial burial of the Miami Terrace platform by deep-water deposits, the subcircular depressions and faults along the platform top were points of origin for a second phase of column growth upward into the deep-water deposits. The continuation of deep platform cavity collapse and column evolution produced pockmarks along paleo-seafloors within the deep-water deposits and at the present-day sea floor. The Plio-Quaternary pockmarks formed at water depths too deep to suggest an origin related to meteoric karst above or near sea level, but rather their formation is suggested to be related to cyclic sea level falls that drove increased groundwater head and density gradients, and seafloor discharge of offshore freshened groundwater sourced from the underlying platform. Plausibly, mixing of freshened groundwater and seawater at the seafloor discharge sites drove dissolution of the host deep-water deposits, which together with erosion by groundwater venting and current scouring formed the pockmarks. Seaward of the Plio-Quaternary seafloor pockmarks, at the late-middle Miocene upper slope of the Miami Terrace platform and along the regional karst/drowning unconformity is a slope-parallel band of ∼189 densely distributed subcircular seafloor depressions with diameters up to 1334 m at water depths up to ∼660 m. It is plausible that along the upper slope, faults and fractures produced by gravity-driven slope instability and possibly tectonics formed a dense network of fluid passageways that promoted upward artesian freshened groundwater flow to sites of discharge where mixing with seawater generated limestone dissolution and the depressions. But tectonic uplift may have forced emersion and initial meteoric sinkhole formation circa late Miocene with later enhancement by freshened groundwater discharge and bottom current erosion.
Re-evaluating Marine Isotope Stage 5a paleo-sea-level trends from across the Florida Keys reef tract
Unraveling how Global Mean Sea Level (GMSL) fluctuated during past warm periods can improve our understanding of linkages between sea-level fluctuations, orbital forcing, and ice-sheet dynamics. Current estimates of GMSL for Marine Isotope Stages (MIS) 5a and 5c — two warm intervals following the relatively well-documented MIS 5e — contain meters of uncertainty and fewer data due to several challenges. These challenges include concealment of datable in-situ coral facies by MIS 1 deposits and inaccessibility due to submergence by modern sea level. We present a comprehensive dataset based on U–Th dating and stratigraphic correlation of 23 cores totaling over 170 m of recovered coral-reef deposits across the tectonically stable Florida Keys Reef Tract (FKRT). Following detailed facies descriptions, 34 in-situ, minimally altered aragonitic coral samples (≤2.7% calcite) below the Holocene-Pleistocene boundary were targeted for U–Th geochronology. Fourteen closed-system coral U–Th ages from MIS 5a include the commonly used sea-level indicator Acropora palmata , but also the massive coral taxa Pseudodiploria strigosa, Siderastrea siderea, Orbicella spp., and Porites astreoides . Dating yielded ages in the range of 88–81 ka (average 2σ uncertainty of less than 200 years). These ages suggest MIS 5a reef initiation at ∼88 ka BP, a peak near 83 ka with minimum elevations between −6.0 ± 0.5 and −5.6 ± 0.5 m MSL (2σ uncertainty and subsidence-corrected), and reef termination and sea-level fall by ∼81 ka BP. Notably, the range of peak MIS 5a relative sea-level estimates of −6.5 to −5.1 m MSL are more than 2 m shallower (higher) than previous estimates of −11 to −9 m. Our higher resolution regional sea-level reconstruction across four subregions of the Florida Keys reef tract aligns with changes in July insolation at 65° N: a trend that most other records, such as deep-sea sediments, do not have the accuracy and precision to resolve. Three massive coral samples from MIS 5c, consisting of Pseudodiploria clivosa , and Orbicella spp., yielded ages in the range of 104 to 99 ka (average 2σ uncertainty less than 200 years); however, because only one sample met the closed-system criteria, our ability to estimate MIS 5c sea level is relatively limited. More empirical estimates of sea-level from the MIS 5a and MIS 5c intervals based on numerical dating of reliable local sea-level constraints are critical for GMSL calculations and relating changes in sea-level amplitude and timing to global ice volume modeling and glacio-isostatic effects, all of which can improve predictions of future sea-level changes in coastal regions.
Tracking mangrove condition changes using dense Landsat time series
Mangroves in tropical and subtropical coasts are subject to episodic disturbances, notably from severe storms, leading to potential widespread vegetation mortality. The ability of vegetation to recover varies, and with disturbances becoming more frequent and severe, it is vital to track and project vegetation responses to support management and policy decisions. Prior studies have largely focused on binary mangrove mapping (i.e., presence or absence), while tracking conditions and condition change have not received sufficient attention. In this paper, we demonstrate a method based on dense time series Landsat images for continuous monitoring of mangrove conditions, where we track three kinds of post-disturbance mangrove conditions, including disturbed (disturbed, with rebound to the previous state within one growing season), recovering (undergoing natural recovery in longer than one growing season), and declining (showing long-term decline after disturbance). The method starts with disturbance detection using the DEtection and Characterization Of the tiDal wEtland change (DECODE) algorithm, an existing dense time series model designed to detect disturbances in tidal wetlands with adaptation to tidal fluctuations. This algorithm is well suited for the detection of tidal wetland disturbances but does not provide satisfactory post-disturbance monitoring results, due to the substantial variability in post-disturbance Landsat observations. To better monitor post-disturbance conditions, a new time series fitting approach, DECODER (DECODE and Recovery), is proposed for the recovery stage. Additionally, for temporal segments divided by disturbance events, we built a random forest classifier with temporal-spectral variables derived from the time series model to characterize mangrove conditions. Employing this approach in Florida's mangroves, we generated condition maps, such as dieback and recovery, with an overall accuracy of approximately 97.96 ± 0.86- [95 % confidence intervals]. Comparing post-hurricane conditions in Florida revealed that the increased frequency and severity of disturbances are challenging mangrove resilience, potentially diminishing their ability to recover and sustain ecosystem functions.
Assessing microplastics contamination in unviable loggerhead sea turtle eggs
Sea turtles, in comparison with marine mammals, sea birds, and fishes, are the most affected by microplastics in terms of number of individuals impacted and concentration within each organism. The ubiquitous nature and persistence of microplastics in the environment further compromises sea turtles as many species are currently vulnerable, endangered, or critically endangered. The objective of this study was to quantify microplastic contamination in unviable loggerhead sea turtle eggs (Caretta caretta). Eggs were collected from seven locations along the northwest coast of Florida. A total of 70 nests and 350 eggs were examined. Microplastics (n = 510) were found in undeveloped loggerhead sea turtle eggs across all seven sites, suggesting that maternal transference and/or exchange between the internal and external environment were possible. The frequency found was 7.29 ± 1.83 microplastic pieces per nest and 1.46 ± 0.01 per egg. Microplastics were categorized based on color, shape, size, and type of polymer. The predominant color of microplastics were blue/green (n = 236), shape was fibers (n = 369), and length was 10–300 μm (n = 191). Identified fragments, films, beads and one foam (n = 187) had the most common area of 1–10 μm2 (n = 45). Micro-Fourier Transform Infrared (μ-FTIR) spectroscopy analysis demonstrated that polyethylene (11 %) and polystyrene (7 %) were the main polymer types. For the first time microplastics were found in unviable, undeveloped loggerhead sea turtle eggs collected in northwest Florida. This work provides insight into the distribution patterns of microplastic pollutants in loggerhead sea turtle eggs and may extend to other species worldwide.
Detection of periodic peaks in Karenia brevis concentration consistent with the time-delay logistic equation
The logistic equation models single-species population growth with a sigmoid curve that begins as exponential and ends with an asymptotic approach to a final population determined by natural system carrying capacity. But the population of a natural system often does not stabilize as it approaches carrying capacity. Instead, it exhibits periodic change, sometimes with very large amplitudes. The time-delay modification of the logistic equation accounts for this behavior by connecting the present rate of population growth to conditions at an earlier time. The periodic change in population with time can progress from a monotonic approach to the carrying capacity; to oscillation around the carrying capacity; to limit-cycle periodic change; and, finally, to chaotic change.
Microbial diversity, genomics, and phage–host interactions of cyanobacterial harmful algal blooms
The occurrence of cyanobacterial harmful algal blooms (cyanoHABs) is related to their physical and chemical environment. However, less is known about their associated microbial interactions and processes. In this study, cyanoHABs were analyzed as a microbial ecosystem, using 1 year of 16S rRNA sequencing and 70 metagenomes collected during the bloom season from Lake Okeechobee (Florida, USA). Biogeographical patterns observed in microbial community composition and function reflected ecological zones distinct in their physical and chemical parameters that resulted in bloom “hotspots” near major lake inflows. Changes in relative abundances of taxa within multiple phyla followed increasing bloom severity. Functional pathways that correlated with increasing bloom severity encoded organic nitrogen and phosphorus utilization, storage of nutrients, exchange of genetic material, phage defense, and protection against oxidative stress, suggesting that microbial interactions may promote cyanoHAB resilience. Cyanobacterial communities were highly diverse, with picocyanobacteria ubiquitous and oftentimes most abundant, especially in the absence of blooms. The identification of novel bloom-forming cyanobacteria and genomic comparisons indicated a functionally diverse cyanobacterial community with differences in its capability to store nitrogen using cyanophycin and to defend against phage using CRISPR and restriction-modification systems. Considering blooms in the context of a microbial ecosystem and their interactions in nature, physiologies and interactions supporting the proliferation and stability of cyanoHABs are proposed, including a role for phage infection of picocyanobacteria. This study displayed the power of “-omics” to reveal important biological processes that could support the effective management and prediction of cyanoHABs.