Short Term Effects of Hurricane Irma and Cyanobacterial Blooms on Ammonium Cycling Along a Freshwater–Estuarine Continuum in South Florida

Lacustrine and coastal systems are vulnerable to the increasing number and intensity of tropical storms driven by climate change. Strong winds associated with tropical storms can mobilize nutrients in sediments and alter nitrogen and phosphorus cycling, leading to ampliﬁcation of preexisting conditions, such as eutrophication and cyanobacterial blooms (cyanoHABs). In 2016, Florida declared a State of Emergency within and downstream of Lake Okeechobee (LO) due to toxic cyanobacterial blooms (primarily Microcystis ). The blooms originated in LO, but ﬂood control measures released water from LO to the brackish St. Lucie Estuary (SLE). In September 2017, Hurricane Irma traversed the Florida peninsula with sustained winds exceeding 160 km h − 1 , generating torrential rains over the watershed. We quantiﬁed ammonium (NH + 4 ) regeneration and potential uptake rates, and Microcystis toxin gene ( mcyD ) abundance in LO and SLE during the massive bloom in July 2016, the bloom in August 2017 (2 weeks before Irma), and 10 days after Hurricane Irma landfall. In 2016, cyanoHABs were present in both LO and SLE, and potential NH + 4 uptake rates were high in both systems. In 2017, the bloom was constrained to LO, potential NH + 4 uptake rates in LO exceeded those in SLE, and mcyD gene abundance was greater in LO than SLE. Post Hurricane Irma, potential NH + 4 uptake rates decreased signiﬁcantly in LO and SLE, while mcyD gene abundance decreased in LO and increased slightly in SLE. Average NH + 4 regeneration rates could support 25–40% of water column potential NH + 4 demand in the lake and, when extrapolated to the entire LO water column, exceeded external nitrogen loading. These results emphasize the importance of internal NH + 4 recycling for bloom expansion and toxicity in the lake and downstream estuaries. In 2018, the cyanobacterial bloom in the Okeechobee region was one of the largest recorded and is presumed to be driven by the aftermath of Hurricane Irma. Large-scale blooms have also been observed in SLE, likely due to LO ﬂushing and decreased salinity post-hurricane. Thus, results from this study support predictions that increased frequency and strength of tropical storms will lead to more intense blooms in aquatic systems.


INTRODUCTION
Anthropogenically-driven climate change has major effects on aquatic systems globally. Additionally, anthropogenic alterations to the nitrogen (N) cycle via synthetic fertilizer production (Haber-Bosch process) have quadrupled the amounts of chemically reduced N in the last 60 years (Erisman et al., 2015). Together, increased N runoff and climate change contribute to eutrophication and cyanobacterial harmful algal blooms (cyanoHABs; Paerl et al., 2016;Glibert, 2017) in freshwater and coastal systems. Internal loading of N and phosphorus (P) from sediments and recycling processes within the water column can enhance cyanoHABs (Havens et al., 2001a;McCarthy et al., 2007McCarthy et al., , 2016, especially in shallow lakes (Havens et al., 2001a;Moss et al., 2003;Jeppesen et al., 2007;James et al., 2009). However, external nutrient loading from non-point agricultural runoff and point sources is the major driver of internal nutrient loading/recycling and cyanoHABs in eutrophic lakes (Huisman et al., 2018).
Lake Okeechobee in southern Florida is a large (1800 km 2 ), shallow lake (mean depth = 2.7 m) and has experienced cyanoHABs for decades (Havens et al., 1994). Lake water levels are closely monitored and managed by the US Army Corps of Engineers. These activities allow the lake to fill in the winter dry season, supplying water for agricultural and urban uses, and lake water is released in spring to increase storage capacity for the summer wet season (Julian and Osborne, 2018). The St. Lucie Estuary (SLE) is subjected to high freshwater and nutrient inputs from Lake Okeechobee through the South Fork of the estuary and Indian River Lagoon. Lake Okeechobee and the SLE receive high external N and P loads from agricultural runoff and anthropogenic activities (James et al., 2011;Phlips et al., 2012), and recent toxic cyanoHABs have been welldocumented in both systems (e.g., Phlips et al., 2012;Kramer et al., 2018). In 2016, a large, toxic cyanobacterial bloom occurred in Lake Okeechobee and SLE, leading to a State of Emergency declaration in Florida (Kramer et al., 2018). Increased rainfall in 2016, corresponding with an El Niño event, led to increased N concentrations in the water column and decreased salinity in the estuary, conditions ideal for development of a Microcystis bloom (Kramer et al., 2018). Environmental degradation in SLE (e.g., cyanoHABs, fecal bacteria, degradation of nearshore reefs) has also been attributed to on-site sewage disposal systems (septic tanks; Lapointe et al., 2012Lapointe et al., , 2017, which have led to nutrient enrichment and microbial contamination (Lapointe et al., 2017).
Cyanobacteria community structure in Lake Okeechobee has shifted from N-fixing to non-N fixing taxa (Microcystis) in recent decades (Havens et al., 2003). Eutrophication management efforts in Lake Okeechobee included 40% reductions of P loading and reducing back-pumping from the Everglades Agricultural Area (EAA; James et al., 2011). The Okeechobee Basin Management Action Plan (BMAP) was also implemented in 2014 to reduce daily loads for total P (Zhang et al., 2016), and these approaches will likely also reduce N loads. However, non-point source pollution remains a problem in the Lake Okeechobee watershed, and concerns about water-level and flood control management activities persist (Kramer et al., 2018).
Microcystis is a potential toxin-producing cyanobacterial genus and strong competitor for chemically reduced forms of N (Blomqvist et al., 1994;Yang et al., 2017;Hampel et al., 2018), such as ammonium (NH + 4 ) and urea. Microcystin (MC), synthesized by Microcystis and other cyanobacteria, is a potent hepatotoxin and is N-rich (10 N atoms per molecule; Gobler et al., 2016). In fact, reduced forms of N can yield higher production of MC (Monchamp et al., 2014;Davis et al., 2015;Harke and Gobler, 2015). Ammonium is the preferred source of N for most primary producers, including non-N-fixing cyanobacteria (Blomqvist et al., 1994;Beversdorf et al., 2015). Thus, internal cycling of NH + 4 is of particular importance in lakes affected by non-diazotrophic cyanoHABs, and internal NH + 4 regeneration often sustains late summer blooms (Paerl et al., 2011;McCarthy et al., 2013;Hampel et al., 2019). NH + 4 turnover rates in Lake Okeechobee are rapid, while ambient NH + 4 concentrations are generally low, suggesting that NH + 4 is in high demand (McCarthy et al., 2009;James et al., 2011). Large discrepancies have been observed between inflow and outflow N in previous mass balance studies for Lake Okeechobee (Havens et al., 2001a), suggesting that internal N cycling processes play a critical role in supporting primary productivity and cyanoHABs (James et al., 2011).
Current climate change models forecast increased severity of extreme weather events, including hurricanes and typhoons (Bender et al., 2010;Knutson et al., 2010;Paerl et al., 2018). Hurricane events increase rainfall in affected areas, leading to higher runoff and nutrient concentrations (Havens et al., 2001b;James et al., 2008). Strong winds associated with tropical storms also enhance eutrophication by mobilizing nutrients in sediments via resuspension (Havens et al., 2001b;James et al., 2008;Ding et al., 2012). Previous studies in Lake Okeechobee reported doubled to quadrupled concentrations of total N (TN) and total phosphorus (TP), along with major increases in soluble reactive phosphorus (SRP), NH + 4 , and other forms of N, after hurricanes (James et al., 2008;Ding et al., 2012). Increased nutrient availability after these events leads to amplified productivity, biomass, and chlorophyll levels (Fogel et al., 1999;Paerl et al., 2001).
The objectives of this study were to (1) quantify water column NH + 4 dynamics during a severe cyanoHAB in 2016 that affected both Lake Okeechobee and the SLE and (2) compare those rates to those from 2017, when the cyanoHAB was restricted to Lake Okeechobee. We also aimed to (3) quantify these rates within 2 weeks after a major hurricane in September 2017 to examine post-hurricane NH + 4 cycling rates in both the lake and estuary. Finally, we (4) quantified MC synthetase gene (mcyD) abundance before and after the hurricane to investigate the effects of hurricane passage on MC producing cyanobacteria. We hypothesized that NH + 4 cycling rates in Lake Okeechobee would be higher during the large bloom in 2016 than in 2017 during the smaller bloom. We also hypothesized that NH + 4 uptake and regeneration rates in the estuary in 2016 would decrease along the salinity gradient and be higher than in 2017. Lastly, we anticipated that the hurricane would have a major effect on NH + 4 cycling and Microcystis in both systems, and that regeneration rates would increase post-disturbance, while uptake rates and abundance of MC synthetase would decrease.

Sample Collection
Lake Okeechobee (LO) and SLE water samples were collected on three occasions: July 25-27, 2016;August 22-24, 2017;and September 20-21, 2017. Sampling in July 2016 followed a major cyanobacterial bloom in LO that extended into SLE and reached the Atlantic Ocean (Kramer et al., 2018). The August 2017 sampling occurred shortly before Hurricane Irma passed over the lake as a Category 3 hurricane (sustained winds 178-208 km h −1 at landfall near Marco Island) on September 10. The September 2017 sampling occurred 10 days after Hurricane Irma passed through central Florida and the Okeechobee region (sustained winds ∼40 mph on September 10, DBHydro, SFWMD).
In July 2016, sampling was conducted at two stations in LO and four stations in SLE: L004 at surface and bottom (∼3 m) water depths; LZ40 at surface and bottom (∼4 m) water depths; and SLE80, SLE2, SLE4, and SLE8 (all SLE sampling was surface water; Figure 1). Bottom water samples were collected about 0.5 m above the sediment-water interface. L004 is located in eastern LO, ∼8.5 km from the St. Lucie Canal, and LZ40 is in the center of the lake. Stations in SLE followed a salinity gradient ( Table 2), with SLE80 closest to the canal discharge lock (S308; Figure 1), and SLE8 closest to the Atlantic Ocean. In August 2017, due to poor weather, only L004 (surface and bottom) and a northern LO station (SAV 165, surface water only) were sampled. Similarly, SLE stations (SLE7, SLE5, and SLE80) could only be sampled from land due to poor weather conditions. In September 2017, stations SLE5 and SLE7 were sampled in SLE, and stations L004 (surface), LZ40 (surface), and LOBG in southeastern LO were sampled (Figure 1).
Water for in situ nutrient analyses was filtered in the field using 0.2 µm syringe filters into 15 ml polypropylene tubes and frozen upon return to the laboratory. Geophysical parameters (temperature, dissolved oxygen, pH, salinity, and conductivity) were measured using a YSI multi-parameter sonde in July 2016 and August 2017 and a Manta 2 sonde (Eureka Waterprobes) in September 2017. Water for NH + 4 dynamics experiments was collected into pre-washed and sample-rinsed 6 L cubitainers and returned to the lab within 3 h for incubations.
Samples for chlorophyll a (Chl a), TN, and TP analyses were collected and analyzed by South Florida Water Management District (SFWMD) according to standard USEPA methods (as described in SFWMD, 2002). Dissolved nutrient analyses included NH + 4 , NO − 2 , NO − 3 , orthophosphate (OP), and urea and were analyzed using a Lachat Quikchem 8500 FIA nutrient analyzer according to the manufacturer's instructions.

NH + 4 Regeneration and Potential Uptake
Water column NH + 4 regeneration and potential uptake experiments followed the protocol described in Hampel et al. (2018) with modifications for 15 N analysis. Briefly, 1 L of water collected at each station was amended with 98% 15 NH 4 Cl (Isotec; final concentration added: 16 µM for all stations except SAV165 = 32 µM). Amended water was mixed thoroughly and decanted into six, 125 ml clear polystyrene Nalgene bottles (triplicates for light and dark incubations). Initial samples were filtered (0.2 µm syringe filter) immediately after spiking (T0) into 15 ml clear, polypropylene tubes (for total NH + 4 ) and 12 ml gas-tight exetainers (for 15 N analysis). Dark incubation bottles were wrapped in aluminum foil, and all samples were placed in a standard lab incubator (also used for algal cultures), at nearambient temperature (29 • C) and set to the appropriate diurnal light cycle, for 20 h. After incubation, final samples (Tf) were processed as described for the T0 samples and stored frozen until analysis. Total NH + 4 ( 14+15 N) concentrations were determined using the Lachat Quikchem 8500 FIA nutrient analyzer, and 15 NH + 4 concentrations were determined using the OX-MIMS method, combining membrane inlet mass spectrometry (MIMS; Kana et al., 1994) with NH + 4 reduction to N 2 gas (Yin et al., 2014). Samples for OX-MIMS were treated with 200 µl of hypobromite iodine solution and immediately measured on the MIMS. 15 NH + 4 concentrations were determined using the line equation from the standard curve (0.1-100 µM 15 NH + 4 ) and total 15 N 2 production (Yin et al., 2014). Potential uptake and actual regeneration rates were calculated using the Blackburn/Caperon model (Blackburn, 1979;Caperon et al., 1979). Volumetric NH + 4 uptake rates (light and dark), but not regeneration rates (actual rates; Gardner et al., 2017), reported in this study are qualified as potential rates due to saturating additions of substrate, which can alter steady-state conditions (Glibert, 1988). However, in eutrophic systems, results from saturating-and tracer-level isotope additions tend to converge (Glibert, 1988).

DNA Collection and Extraction
Environmental DNA for gene abundance analysis was collected in August and September 2017 using 0.2 µm Sterivex filters (EMD Millipore, Burlington, MA, United States) and preserved with Ambion RNAlater (Invitrogen, Carlsbad, CA, United States). In August, approximately 120-240 ml of water was pushed through a Sterivex filter. However, due to highly turbid waters after the hurricane in September, only 45-60 ml of water were filtered for stations L004, LZ40, and SLE, and 300 ml for station LOBG. Preserved filters were frozen at −80 • C. DNA was extracted using the Gentra PureGene kit (Qiagen, Inc., United States). Residual RNAlater in the Sterivex filters was removed by pushing 10 ml of Phosphate Buffer Saline 1X Solution (Fisher BioReagents, United States) through the filter. Lysis buffer (0.9 ml) and Proteinase K (10 µl) were added to the filters, followed by 1 h incubation at 55 • C and 1 h incubation at 65 • C (Newell et al., 2011). Concentration and purity of extracted DNA were measured spectrophotometrically (Nanodrop 2000, Thermo Scientific).

qPCR Analysis
The MC synthetase gene was amplified using the mcyD-F2 and mcyD-R2 primers, targeting a 298 bp region of the mcy operon (Kaebernick et al., 2000). qPCR standards were prepared by   Davis et al. (2009) for mcyD (95 • C initial denaturation for 2 min, 95 • C denaturation for 15 s, 50 • C annealing for 1 min, and 60 • C extension for 1 min; 45 cycles) followed by the melting curve. Automatic settings for the thermocycler (Realplex, Eppendorf) were used to determine the threshold cycle (Ct values), efficiency (98%), and a standard curve with R 2 values of 0.99. mcyD gene copies were calculated as (ng * number mol −1 )/(bp * ng g −1 * g mol −1 of bp) and is reported in gene copies ml −1 of sample water.

Statistical Analysis
All statistical analyses were performed using RStudio software (version 1.1.383). Environmental data were checked for normality using the Shapiro-Wilk normality test. After ensuring that the data were not normally distributed, the Kruskal-Wallis test for non-parametric data was used to determine significant differences between sites and time points. The Kendall correlation method for non-parametric data was used between environmental variables and geochemical rates and mcyD gene abundance. Kendall's p-values are usually more accurate for datasets with small sample sizes than Spearman's (Croux and Dehon, 2010). The correlation matrix presented below was constructed showing the Kendall's tau correlation coefficient and p-values (p < 0.05 was considered statistically significant).

St. Lucie Estuary
Water temperatures in SLE also were consistent between sampling events (29.4-31.1 • C; Table 2). Salinity in July 2016 increased from SLE80 (closest to the canal; 0.19) to SLE8 (closest to the Atlantic Ocean; 8.9), and the same general pattern was observed in August 2017, with lower salinity at SLE5 (4.6) than SLE7 (6.75). However, after the hurricane, salinity decreased significantly to 0.19 (SLE5) and 0.22 (SLE7; p < 0.05). Similar to LO, highest conductivity was observed in September 2017 and was significantly different from other sampling events (436 ± 30.5 µs cm −1 ; p < 0.001). Highest mean Chl a values were recorded in August 2017 (18 ± 1.9 µg L −1 ).

Lake Okeechobee
As mentioned previously, light and dark uptake rates are qualified as potential rates because of added substrate. However, regeneration rates are considered actual rates because they represent 14 NH + 4 turnover (Gardner et al., 2017). In July 2016, potential NH + 4 uptake rates in the light ranged from 0.54 to 2.24 µM h −1 , and peak light NH + 4 uptake rates were observed at L004 in surface water (Figure 2A). Light NH + 4 uptake rates in August 2017 (mean = 2.88 ± 0.72 µM h −1 ; Figure 2B) were greater than in July 2016, with peak rates in the central-eastern part of the lake (L004; Figure 2B). After the hurricane (Figure 2C), light uptake rates significantly decreased (mean = 0.66 ± 0.27 µM h −1 ; p < 0.05) and showed spatial variability; highest light uptake was observed at the shore station (LOBG) and lowest in the central lake (LZ40).
Dark NH + 4 uptake rates were significantly lower than light rates on all occasions (p < 0.001). In July 2016, dark NH + 4 uptake rates were comparable between stations and depths (mean = 0.30 ± 0.05 µM h −1 ; Figure 2A). Dark rates in August 2017 were not statistically different from the July 2016 rates. However, dark uptake rates in September 2017 were significantly lower than in July 2016 (mean = 0.12 ± 0.02 µM h −1 ; Figure 2C), with highest dark uptake in the central-eastern lake (L004; Figure 2C) and lowest at the shore station (LOBG).

St. Lucie Estuary
Light NH + 4 uptake rates in the estuary were significantly different between sampling events and years (p < 0.05). Mean light uptake rates in the estuary were highest in July 2016 (mean = 1.76 ± 0.25 µM h −1 ), with highest rates close to the canal (SLE80 and SLE2) and decreasing toward the ocean (Figure 2A). In August 2017, light uptake rates were significantly lower (mean 0.89 ± 0.22 µM h −1 ; p < 0.05) but followed the same pattern, with peak uptake rates at the lock (SLE80) and decreasing at SLE5 and SLE7 (Figure 2B). After the hurricane, light uptake in the estuary decreased even further and was similar between stations SLE5 and SLE7 (mean = 0.22 ± 0.03 µM h −1 ; Figure 2C).
Dark NH + 4 uptake rates in the estuary were not statistically different between sampling events but were significantly lower than the light NH + 4 uptake rates (p < 0.001). In July 2016, dark rates in the estuary ranged from 0.12 to 1.31 µM h −1 , with highest dark uptake near the lock at SLE2 and lowest further down the salinity gradient (SLE8; Figure 2A). In August 2017, dark uptake rates were lower than in 2016 (mean = 0.09 ± 0.03 µM h −1 ) but followed a similar pattern, with peak rates observed at SLE80 (Figure 2B). In contrast to light uptake rates, dark uptake rates in the estuary in September 2017 were slightly greater than before the hurricane, but not statistically different (0.16 ± 0.08 µM h −1 ; Figure 2C).
Light NH + 4 uptake rates in LO and SLE were positively correlated with ambient temperature and negatively correlated with NO − 3 , TN, and TP concentrations ( Table 3). Dark NH + 4 uptake rates were positively correlated with DO concentration and mcyD abundance and negatively correlated with NO − 2 , OP, urea, and TP concentrations ( Table 3).

St. Lucie Estuary
NH + 4 regeneration rates in the estuary were not statistically different between sampling events, and peak regeneration rates were observed in July 2016 (mean = 0.38 ± 0.15 µM h −1 ; Figure 2A). NH + 4 regeneration rates followed the same pattern as uptake rates in July 2016, with highest values observed at sites closest to the canal (SLE2 and SLE80) and decreasing toward the Atlantic Ocean (SLE8). In August 2017, regeneration rates at the lock (SLE80; 0.46 ± 0.08 µM h −1 ) were similar to 2016 regeneration rates and decreased slightly along the salinity gradient (SLE5 and SLE7; Figure 2B). After the hurricane, NH + 4 regeneration rates at SLE5 and SLE7 (mean = 0.16 ± 0.01 µM h −1 ) remained similar to rates in August 2017. Ammonium regeneration rates in LO and SLE were positively correlated with temperature and Chl a and negatively correlated with salinity, NO − 3 , and TP concentrations ( Table 3).

St. Lucie Estuary
Abundance of mcyD in SLE in August 2017 (mean = 9.67 ± 3.82 × 10 3 copies ml −1 ) was significantly lower than in LO (Figure 3A; p < 0.05). Highest abundance of mcyD in August in SLE was observed at SLE80 near the canal, while lowest abundance was observed further along the salinity gradient at SLE5 (Figure 3A). After the hurricane, mcyD abundance in SLE (mean = 2.69 ± 1.14 × 10 4 copies ml −1 ; Figure 3B) was slightly greater than in August, but not statistically different (p > 0.05).

NH + 4 Uptake During 2016 and 2017 Microcystis Blooms
Increased rainfall in 2016, an El Niño year, led to increased water levels in LO, requiring the release of billions of gallons of lake water through the St. Lucie Canal into SLE (Figure 4A; Kramer et al., 2018). Release of high nutrient freshwater into the brackish estuary created ideal conditions for a Microcystis bloom to develop (Kramer et al., 2018). In 2017, the net flow of water from LO and SLE was lower until Hurricane Irma made landfall in September ( Figure 4B). We compared NH + 4 cycling rates in LO and SLE during two summer blooms and 10 days after Hurricane Irma disturbance.
In July 2016 (El Niño year), potential NH + 4 uptake rates were comparable between the lake and the estuary (Figure 2A). Bloom conditions were present in both systems, with high abundances of Microcystis (95% of algal biomass) and high MC and Chl a concentrations (Kramer et al., 2018). In the estuary, potential NH + 4 uptake decreased following a gradient from the lock (SLE80, SLE2) into the saline part of the estuary (SLE8; Figure 3A). This eastward decrease in NH + 4 dynamics corresponded with decreasing Chl a values in this study (Table 1), as well as decreased Chl a and MC concentrations and mcyE abundance reported previously (Kramer et al., 2018). Light NH + 4 uptake rates exceeded dark uptake rates in both systems, suggesting that NH + 4 demand was dominated by photoautotrophs.
In contrast, without water releases from LO in 2017 (Figure 4B), potential NH + 4 uptake rates in the lake were significantly higher than in the estuary ( Figure 2B). However, NH + 4 uptake rates, Chl a concentrations, and mcyD abundances ( Figure 3A) in SLE followed a similar pattern as in 2016, with highest values upstream (SLE80) and decreasing toward the coast. These results reflect the influence of water releases for flood control and support previous work showing that flood prevention practices in the Okeechobee region drive cyanoHABs in the downstream estuary, even when it may be beneficial for the trophic status of the lake (Lapointe et al., 2012;Phlips et al., 2012;Julian and Osborne, 2018).
Interestingly, mcyD abundance showed a strong, positive correlation with dark NH + 4 uptake rates. Microcystis has been shown to synthesize MC at night (Penn et al., 2014;Davenport et al., 2019), suggesting that MC has a role in Microcystis diurnal metabolism. Dark NH + 4 uptake rates observed here may be involved in MC production, but night time MC synthesis has not been extensively studied and merits future research, including diurnal transcriptomics and gene expression experiments.

Hurricane Irma
After Hurricane Irma, potential NH + 4 uptake rates in both systems significantly decreased ( Figure 2C). Uptake rates at pelagic stations (L004 and LZ40) were the lowest measured in LO in this study (<0.5 µM h −1 ) and suggest decreased algal biomass and activity. In contrast, uptake rates close to the shore (LOBG) were significantly higher than in the central lake (1.2 ± 0.07 µM h −1 ), despite high ambient NH + 4 concentrations ( Table 1). Algal biomass may have accumulated along shorelines, driven by winds and seiches. Post-hurricane abundance of mcyD in LO also decreased significantly, and mcyD was most abundant at the shoreline station (LOBG, Figure 3B), where peak NH + 4 uptake rates were also recorded ( Figure 2C).
NH + 4 uptake rates and Chl a concentrations in SLE after the hurricane were significantly lower, but mcyD abundance remained unchanged ( Figure 3B). Following Hurricane Wilma in 2005, the phytoplankton community in LO shifted toward diatoms (James et al., 2008). Under poor light conditions, such as high turbidity and low water column stability, diatoms may have outcompeted cyanobacteria (James et al., 2008). Indeed, turbidity and total suspended solids (TSSs) increased significantly after the hurricane in both LO and SLE (Supplementary Table 1). Cyanobacteria are also susceptible to flushing events , and rapid flushing can also explain the lower Chl a concentrations in SLE ( Table 2; Paerl et al., 2018). Net water flow at S308 increased after Hurricane Irma and remained high for the rest of the year (Figure 4B).
These results represent the short-term (10 days) response of a HAB to a major storm disturbance. Numerous other studies have reported long-term responses of cyanobacteria and phytoplankton to hurricanes/typhoons in other systems (Yannarell et al., 2007;Zhu et al., 2014;Paerl et al., 2018). These studies have reported that increased pulses of nutrient loading, along with wind-induced nutrient resuspension from the sediments, may ultimately stimulate HABs and amplify preexisting eutrophication (James et al., 2008;Zhu et al., 2014). The HAB in the LO and SLE region in 2018 was one of the largest recorded and may have been related to rainfall and high winds from Hurricane Irma and heavy rainfall in spring 2018 (Paerl et al., 2019). Our short-term lake and estuary responses to Hurricane Irma emphasize that studies focusing on both short and long-term effects are necessary for understanding hurricane impacts on shallow, eutrophic ecosystems and cyanoHAB dynamics.

Ammonium Regeneration
Actual NH + 4 regeneration rates in the lake and estuary followed the same general pattern as NH + 4 uptake rates. During the 2016 bloom, internal NH + 4 regeneration rates in the lake were similar to regeneration rates in SLE. In July 2016, actual regeneration could support 43% of the potential light uptake in LO. Similarly, during the 2017 bloom, 38% of light uptake could be supported by regeneration, despite higher uptake rates. Previous NH + 4 dynamics experiments in LO showed that ∼35% of light uptake could be supported by regeneration (James et al., 2011). Following the hurricane, water column regeneration decreased significantly in LO and SLE, and only 24% of NH + 4 uptake in LO could be supported by regeneration, likely due to flushing and decreased biomass post-disturbance.
Total nitrogen (TN) from external loading into LO is estimated at ∼6.36 × 10 3 tons yr −1 . LO can be divided into three distinct regions: littoral zone, nearshore (1-2 m depth), and pelagic zone (offshore; James et al., 2008). Acknowledging the temporal and spatial limitations of the data generated in this study, we extrapolated the measured regeneration rates in the pelagic zone of the lake (stations L004 and LZ40) and nearshore (stations SAV165 and LOBG from this study and L005 from James et al., 2011) to estimate how much NH + 4 may be regenerated in LO annually relative to external N loads. Extrapolating summer bloom values in the pelagic zone of LO (estimated volume = 2.85 km 3 ), 2.86 × 10 5 tons of N yr −1 could be regenerated internally. This value is 45 times greater than the estimated external loading of TN (6.36 × 10 3 tons yr −1 ; James et al., 2009). Nearshore, regeneration could supply 7.94 × 10 3 tons of N yr −1 . These extrapolations of rates measured in summer months are almost certainly an overestimate due to higher productivity and algal biomass. However, the same calculations, using significantly lower regeneration rates from September 2017, showed that 3.37 × 10 4 tons of N yr −1 can be regenerated in the pelagic water column after a major disturbance. While this value is much lower than those reported based on July and August rates, it is still five times greater than annual external TN loading. These results suggest that internal NH + 4 regeneration supplements external TN loading in LO to a large degree, even following a major disturbance, which led to significantly lower regeneration rates. This exercise also illustrates the importance of characterizing internal nutrient dynamics on more extensive temporal and spatial scales to more explicitly constrain annual nutrient budgets.
Results from this study, showing that internal N recycling in the pelagic zone of LO supports a large proportion of primary producer NH + 4 demand, are supported by other studies in hypereutrophic lakes experiencing cyanoHABs, including: Lake Balaton, Hungary (Présing et al., 2001); Lake Biwa, Japan (Haga et al., 1995;Takahashi et al., 1995); Missisquoi Bay, Lake Champlain, United States/Canada (McCarthy et al., 2013); Lake Taihu, China (Paerl et al., 2011;Hampel et al., 2018), and Lake Erie, United States/Canada (Hampel et al., 2019). This study contributes to the growing literature showing that rapid NH + 4 turnover rates can fuel and sustain HABs, despite low NH + 4 concentrations (McCarthy et al., 2007;Hampel et al., 2019). In hypereutrophic Lake Taihu, about 60% of potential light uptake was supported by regeneration during a mid-summer bloom (June), and all potential NH + 4 demand could be supported internally in late summer (August; Hampel et al., 2018). In Sandusky Bay (Lake Erie), summer NH + 4 regeneration in the water column provided bioavailable N equivalent to ∼77% of the annual N load and could support most of the primary producer NH + 4 demand during the HAB season (Hampel et al., 2019).

CONCLUSION
Results from this study provided insight into NH + 4 and HAB dynamics in Lake Okeechobee and St. Lucie Estuary during summer HABs and after Hurricane Irma. We emphasize the importance of NH + 4 cycling on HAB expansion in LO and downstream nutrient pollution in SLE. Results from this study support those from previous studies (Lapointe et al., 2012;Phlips et al., 2012;Kramer et al., 2018), which call for water management strategies that balance flood control and health of downstream estuaries. Previous management strategies (40% P reductions and reduced nutrient inflow from Everglades Agricultural Area; James et al., 2011) have not prevented cyanoHABs in the Okeechobee-St. Lucie system. We show that NH + 4 regeneration plays an important role in Lake Okeechobee in sustaining cyanoHABs, even after a major disturbance. High nutrient loading, combined with internal recycling and climate change, will likely continue to exacerbate cyanoHABs in the Okeechobee-St. Lucie and other freshwaterestuarine systems if further N and P reductions are not implemented. Lastly, short-term lake and estuary responses to Hurricane Irma emphasize that studies focusing on both short and long-term effects are necessary for understanding hurricane impacts on shallow, eutrophic ecosystems and cyanoHAB dynamics.

AUTHOR CONTRIBUTIONS
SN, MM, and JH designed the study and analyzed and interpreted the data. JH and MR collected and analyzed the samples. JH primarily wrote the manuscript. All authors contributed to text and revised the manuscript.

FUNDING
Funding for this work was provided by Florida Sea Grant (Project No. PD-16-10 to MM and SN).