E ﬀ ects of Consecutive Extreme Weather Events on a Temperate Dystrophic Lake: A Detailed Insight into Physical, Chemical and Biological Responses

: Between May and July 2018, Ireland experienced an exceptional heat wave, which broke long-term temperature and drought records. These calm, stable conditions were abruptly interrupted by a second extreme weather event, Atlantic Storm Hector, in late June. Using high-frequency monitoring data, coupled with fortnightly biological sampling, we show that the storm directly a ﬀ ected the stratiﬁcation pattern of Lough Feeagh, resulting in an intense mixing event. The lake restabilised quickly after the storm as the heatwave continued. During the storm there was a three-fold reduction in Schmidt stability, with a mixed layer deepening of 9.5 m coinciding with a two-fold reduction in chlorophyll a but a three-fold increase in total zooplankton biomass. Epilimnetic respiration increased and net ecosystem productivity decreased. The ratio of total nitrogen:total phosphorus from in-lake versus inﬂow rivers was decoupled, leading to a cascade e ﬀ ect on higher trophic levels. A step change in nitrogen:phosphorus imbalances suggested that the zooplankton community shifted from phosphorus to nitrogen nutrient constraints. Such characterisations of both lake thermal and ecological responses to extreme weather events are relatively rare but are crucial to our understanding of how lakes are changing as the impacts of global climate change accelerate.


Introduction
Extreme-weather-related episodic events are becoming more frequent as a consequence of global climate change [1]. Meteorological extremes can be defined as those occasions when temperature (e.g., heat waves or cold spells), precipitation (e.g., floods or droughts) and/or wind speeds (e.g., hurricanes) break long-term records for a specific site with values well above (or below) a threshold of the studied variable [2,3]. Climate change, in turn, has been shown to influence lake thermal dynamics by altering lake mixing regimes [4]. In addition, weather-related episodic events facilitate the export of dissolved and particulate compounds from the surrounding terrestrial ecosystems [5]. As a consequence, lake ecosystems are excellent sentinels of change as they respond directly to meteorological episodic events, incorporating the catchment effects and integrating responses over time at different trophic levels [6]. However, responses to climate change in lake ecosystems operate on various temporal scales, and disentangling driver-response relationships impacting freshwater aquatic ecosystems is complex [7]. For this reason, understanding the responses of single climatic forcing events in the context of longer-term dynamics (e.g., seasonal patterns) is important and requires ecological variables The Burrishoole catchment is situated on the northwest Atlantic coast of Ireland, in County Mayo (Figure 1). The area of the catchment is approximately 90 km 2 and the geology is primarily composed of granite, gneiss and schist [67]. The system comprises 45 km of low-productive and poorly buffered shallow streams draining into the sea in the northeast corner of Clew Bay, on the Atlantic coast [68]. The streams form two main sub-catchments comprising the two main inflow rivers to Lough Feeagh, the Glenamong and Black rivers (Figure 1). The presence of small amounts of dolomite and wacke in the eastern part of the catchment has a marked effect on the buffering capacity of the Black river compared with the Glenamong [69]. A terminal moraine, marking the boundary between metamorphic and sedimentary rock-types, separates the largest freshwater lake, Lough Feeagh, from the lagoonal estuary, Lough Furnace [70] (Figure 1). The land cover in the catchment comprises 64% blanket bog with a peat depth ranging from 0.5 to 2 m [71]. In addition, 23% of the land is used for commercial afforestation. The remaining area comprises reduced extensions of native oak woodland, natural grasslands and agricultural land [70]. The Burrishoole catchment is situated on the northwest Atlantic coast of Ireland, in County Mayo (Figure 1). The area of the catchment is approximately 90 km 2 and the geology is primarily composed of granite, gneiss and schist [67]. The system comprises 45 km of low-productive and poorly buffered shallow streams draining into the sea in the northeast corner of Clew Bay, on the Atlantic coast [68]. The streams form two main sub-catchments comprising the two main inflow rivers to Lough Feeagh, the Glenamong and Black rivers (Figure 1). The presence of small amounts of dolomite and wacke in the eastern part of the catchment has a marked effect on the buffering capacity of the Black river compared with the Glenamong [69]. A terminal moraine, marking the boundary between metamorphic and sedimentary rock-types, separates the largest freshwater lake, Lough Feeagh, from the lagoonal estuary, Lough Furnace [70] (Figure 1). The land cover in the catchment comprises 64% blanket bog Water 2020, 12, 1411 5 of 25 with a peat depth ranging from 0.5 to 2 m [71]. In addition, 23% of the land is used for commercial afforestation. The remaining area comprises reduced extensions of native oak woodland, natural grasslands and agricultural land [70].
The Marine Institute (MI) maintains an Automatic Water Quality Monitoring Station (AWQMS) on Lough Feeagh (Coordinates; 53.95, −9.58) [72] and Automatic River Monitoring Stations (ARMS) on the Glenamong (Coordinates: 53.96, −9.60) and Black (Coordinates: 53.97, −9.58) rivers (Figure 1). The MI also measures water level on the Glenamong river downstream of the ARMS every 15 minutes using a data logger (Orpheus mini, OTT Hydromet, Germany). An established rating curve has been used to convert water level into stream flow. To measure Black river discharges, the 'drainage area ratio method' was used [73]. Feeagh AWQMS collects high-frequency sensor information every 2 minutes. Vertical temperature profiles were measured using 12 platinum resistance thermometers (PRTs: Lab facility PT100 1/10DIN 4 wire sensor, Labfacility Ltd., UK) at 2. 5,5,8,11,14,16,18,20,22,27,32, and 42 m. The AWQMS has a multiparameter probe (Hydrolab Data Sonde 5x, OTT Hydromet, UK) at 0.9 m from the surface to measure temperature and DO. Sensors on the AWQMS are cleaned fortnightly and DO on the sonde is calibrated once per month. Wind speed was measured in the AWQMS (2.3 m above lake surface level) with an anemometer (Vector Instruments A100L2-WR, Windspeed Ltd., UK). Data from the AWQMS is available from 2004 onwards. More information about the deployed sensors in the buoy can be found in the GitHub Data Repository [74]. In addition, Met Éireann maintains an automatic weather station (Newport Furnace) installed at 22 m above mean sea level, on the shore of Lough Feeagh (Coordinates: 53.92, −9.58), from which minimum and maximum air temperature as well as daily rainfall data have been obtained [75] (Figure 1).

Seasonal Chemical and Biological Sampling
Lough Feeagh and its main inflows, the Glenamong and Black rivers, were sampled fortnightly from April to October 2018 during daytime hours (10:00−16:00). The lake was sampled at its deepest point next to the buoy, and the rivers were sampled next to the monitoring stations ( Figure 1). Samples for water chemistry in the two rivers were collected in 2 L pre-acid-washed high-density polyethylene bottles (HDPB), pre-rinsed with river water prior to sample collection. In Lough Feeagh, a 6 m integrated water sample was collected using a Plexiglas tube (length: 6 m, diameter: 70 mm). A quantity of 30 L of lake water was collected and placed in three acid-washed 10 L HDPB bottles which were pre-rinsed with lake water. All water samples (filtered and unfiltered) were kept in dark at 4 • C until further analyses, which were carried out within 24-48 h following collection. Unfiltered lake water from the 6-m integrated samples was placed into 120-mL amber glass bottles and preserved with Lugol's iodine solution for further phytoplankton and ciliate counts. Zooplankton samples were collected using a 53-µm mesh net (diameter 0.3 m) with a flow meter (438115, HYDRO-BIOS, Germany) attached. The net was lowered to a 20-m depth and drawn vertically through the water column at an approximate speed of 1 m·s −1 . After collection, zooplankton were transferred to 500-mL storage containers and anaesthetised with half an Alka-Seltzer ® tablet to prevent ejection of eggs. The samples were then preserved in 70% v/v final concentration industrial methylated spirits (IMS). A fourth zooplankton sample was taken without the preservation step and kept alive for later gut evacuation and preparation for nutrient content analyses.
In addition, during the period of peak stratification (from May to August 2018), monthly water samples (n = 3, specifically on the 28th May, the 26th June and the 24th July) for phytoplankton and ciliate composition were collected at two discrete depths at the same time as the integrated samples to examine vertical patterns in plankton community composition. These samples were taken at 12 m (to represent the more stable, stratified zone below the surface mixed layer) and at 25 m (to represent hypolimnetic water below the thermocline) and were collected with a standard peristaltic pump and a pre-rinsed weighted silicone tube.

Chlorophyll a and Nutrient Analyses
Chlorophyll a analysis was carried out on 2 L triplicate samples from the integrated water sample filtered through GF/F filters and extracted with ethanol, with absorbance read in a spectrophotometer (UV-1800, Shimadzu, Germany) at 665 and 750 nm, with a 1 cm cuvette [76]. Total nitrogen samples were analysed in a certified external laboratory (Aquatic Services Unit, University College Cork, Ireland) with a flow injection analysis system (Quickchem ® 8500, Lachat Instruments, Germany) on unfiltered water after a persulfate oxidation digestion step (method detection limit of 0.1 mg N·L −1 ) [77]. Total phosphorus was analysed on unfiltered water samples using the molybdenum blue method, with absorbance read in a spectrophotometer (UV−1800, Shimadzu, Germany) at 882 nm, with a 5 cm cuvette, following a digestion step with potassium persulfate (K 2 S 2 O 8 ) under pressure and temperature (autoclave for 30 min at 15 psi and 121 • C) (method detection limit of 1 µg P·L −1 ) [76,78]. Total phosphorus analysis was carried out in triplicate using standard curves with certified reagents and quality controls using separate reagents to assess the performance of the analytical method.
Seston carbon (C), nitrogen (N) and phosphorus (P) content was analysed from lake water passed through a 53-µm mesh to remove zooplankton and filtered onto pre-combusted and weighted GF/F filters. For zooplankton community nutrient content analyses, the fourth un-preserved zooplankton sample was transferred to 2 L of filtered (0.45 µm) lake water for approximately four hours to evacuate their gut. After gut evacuation, zooplankton fraction was rinsed with Milli-Q water and divided into six equal, well-mixed sub-samples and filtered onto pre-combusted and weighted GF/F filters. All GF/F filters (seston and zooplankton fractions) were immediately dried at 60 • C for 24 h. Triplicate filters for C and N analyses were measured using an elemental analyser (vario EL cube, Elementar ® , Germany). The standard reagent used for this analysis was Acetanilide (C 8 H 9 NO) with a composition of 10.36% of N and 71.09% of C. Triplicate filters for P content were determined as specified for TP water analysis with an additional pre-step of combusting the filters at 550 • C for 2 h to ensure the seston and zooplankton organic matter would break down during the persulfate digestion step [78,79]. A blank filter per replicate was considered for all C, N and P analyses.

Phytoplankton Community Identification and Enumeration
Identification and counting of preserved phytoplankton samples was carried out following the technique developed by Utermöhl [80]. Preserved phytoplankton samples were well mixed before transferring them to a 25-mL sedimentation chamber. Chambers were settled for at least 18 h. Identification and enumeration were performed under an inverted microscope (Leica ® Microsystems, Ireland). A camera linked with the microscope enabled accurate phytoplankton measurements from images by using a calibrated stage micrometre available in the Leica Application Suite imaging software (LAS V4.12). Two counts were carried out on each sample, a count of large species at x40 magnification, and a high magnification count at x400 for small species. At x40 magnification, the whole chamber was counted to include large species or colonies which may not be numerous enough to be captured by the subsequent count at high magnification. Counts at x400 magnification were carried out over two transects of the chamber base. Samples were identified to genus and species level when possible [81,82]. Biovolumes were then calculated in order to compare the relative contributions of different algae samples [83]. This was done by multiplying the number of cells for each species present in a sample by its average cell volume. Average cell biovolume was calculated by measuring 10-20 individuals of each species from each sample. Fewer individuals were measured for rare species.

Zooplankton Community Identification and Enumeration
Preserved zooplankton samples were sieved into a 53-µm mesh and washed with deionized water to remove IMS. The zooplankton fraction was then transferred into a graduated volumetric cylinder and diluted to 40-100 ml, depending on the density of the sample. A minimum of 10% of the final volume was then identified and counted by transferring a 5-mL sub-sample, using a wide-bore 5-mL pipette, to a Ward rotary chamber [84]. Sub-samples were counted under a stereoscope microscope (SZX16, Olympus ® , Germany). Zooplankton individuals were identified to groups (cladocerans, copepods and rotifers) and to genus level or species level where possible [85,86]. For biomass calculations, the lengths of, at least, the first 50 individuals encountered in a sub-sample were measured using a calibrated ocular micrometre. Copepod total body length was measured from the top of the head to the base of the caudal rami. Cladoceran total body length was measured from the top of the head to the tip of the abdomen, not including spines or projections. Taxa abundance was converted to biomass using regression formulae relating length to weight, calculated for individuals from the most abundant taxa in Lough Feeagh [87], using formulae from previous studies in Irish lakes [88], or from the literature [84,89].

Data Analyses
For the high-frequency data collected on Feeagh buoy (AWQMS), a QA/QC was applied to account for data gaps, calibration periods or any sensor problems [74]. Data was then transformed from sub-daily (recorded every few minutes) to daily data (daily averages) using the subdaily2daily function from the hydroTSM package in R [90]. Schmidt stability was calculated using lake water temperature profile data, Feeagh bathymetry and wind speed as input files to the rLakeAnalyzer package in R [91]. Schmidt stability was used as a physical index indicative of the resistance of Lough Feeagh to mixing (e.g., a value of 0 implies that a water column is fully mixed with units of amount of energy required to fully mix a water column per lake area) [92]. Mixed layer depth (MLD) was defined as the depth where the water density was 0.1 km·m −3 greater than the shallowest sensor depth, which, based on the findings from Wilson et al. (2020), is a robust method for temporal analysis in Lough Feeagh [93]. Quantifying the changes of the estimated MLD of Lough Feeagh over the study period was used to better characterise the physical impacts of the extreme-weather-related events in the lake which, in turn, had implications for chemical and planktonic community dynamics [94]. Based on the summer 2018 lake stability pattern, the period under study was divided into three phases (pre-, intermediate and post-Storm Hector). Pre-storm conditions were delineated by the initial heat wave (28th May-13th June); this was followed by an intermediate period dominated by the arrival and aftereffects of Storm Hector (14th June-3rd July) and, finally, the post-storm period was represented by a return to stable conditions, with high Schmidt stability values (4th-15th of July).
Metabolic estimates were calculated using the maximum likelihood estimation method (mle) from the LakeMetabolizer package in R [95]. Two-minute time step changes of Lough Feeagh epilimnetic DO (0.9 m), water temperature, photosynthetic active radiation (PAR) and wind speed for the study period (April-October 2018) were used for producing daily GPP (positive values), R (negative values) and NEP (GPP + R) metabolism estimates. Daily values of air-water oxygen exchange velocity (k600) were calculated following Cole and Caraco (1998) [96]. Regarding daily metabolic estimates calculated from DO values at 0.9 m of Lough Feeagh, only the sub-period from the 28th May to the 15th July was used, comprising similar days for the pre-, intermediate and post-storm event phases. Days with negative values of GPP or positive values of R were removed from the dataset. The non-parametric Kruskal-Wallis test and a following pairwise multiple comparison Dunn's test for determining which groups differ from others [97,98] were performed to test for statistical differences between periods (pre-, intermediate and post-storm event) on metabolic estimates using the FSA package in R [99].
Point loadings (L) from the main inflow rivers were calculated by multiplying the concentration from the biweekly grab sampling nutrient concentrations by the daily average flow on the day the sample was collected [100] (Equation (1)). In turn, biweekly nutrient ratios of total inflow loadings were calculated as the sum of the Glenamong and Black rivers' nutrient loadings, converted to molar concentrations to make it possible to calculate nutrient ratios. where L is the nutrient loading in tonnes (t) or kilograms (kg) per day; K is the factor to convert milligram or microgram concentrations to t or kg; V d is the daily cumulative streamflow in litres (flow in L·s −1 multiplied by 86,400 s·day −1 ); and C d is the grab sampling nutrient concentration measured on the streams each specific day. Spearman's rank correlation coefficient and significance of the correlation between inflow and in-lake TN:TP ratios were carried out with the Chart.correlation function of the PerformanceAnalytics package in R [101]. Nutrient N:P imbalances were calculated as the N:P of the seston fraction minus the N:P of the zooplankton community fraction. Elemental nutrient imbalances can provide an indication of the magnitude of potential food quality constraints for the zooplankton community. Therefore, elemental imbalances close to 0 would indicate no food quality constraints, whereas values at variance from 0 (positive or negative values) indicate possible food quality constrains for the zooplankton community [61]. Shannon diversity index for phytoplankton and ciliate community was calculated with the diversity function of the vegan package in R [102] using abundance data based on the lowest taxonomic group possible [103]. All R analyses were performed in version 3.6.1 [104].

Physical and Nutrient Dynamics
During summer 2018, high pressure (>1013 mbar) dominated over Ireland and, as a consequence, the country experienced unusually dry and warm conditions from the middle of May to the 12th June compared to the previous 15-year average during this time period (Figure 2A,B). Between the 13th and the 19th June, Storm Hector brought heavy rain and strong winds across Ireland as well as a significant drop in mean daily air temperatures ( Figure 2). There was an accumulation of 64 mm of precipitation during these 7 days, which occurred after 11 days and was followed by 20 days with no rain (0 mm·day −1 ) ( Figure 2B). On the 14th June 2018, coinciding with the arrival of Storm Hector, wind speed reached a daily mean value of 9.6 m·s −1 ( Figure 2C), with maximum gusts of 26.91 m·s −1 [75]. After Storm Hector, from the 21st June to the 12th July, high pressure stretching from the Azores to Scandinavia brought very warm and dry conditions, during which a maximum air temperature of 29.8 • C on the 27th June was recorded at the study site (Newport Furnace, Met Éireann weather station) ( Figure 2A).
Thermal stratification was recorded in Feeagh from mid-May as a result of the increasing air temperatures and calmer weather conditions (low precipitation and wind speeds) recorded at that time. This resulted in a sudden increase in Schmidt stability (St), reaching values above 400 J·m −2 from the 1st June and a maximum of 548.25 J·m −2 recorded on the 7th June ( Figure 3A). With the arrival of Storm Hector on the 14th June, strong winds and a decrease in surface heat input from the atmosphere led to a decline in St values to 192.52 J·m −2 (three-fold reduction). After two weeks, stable stratification was again observed in the Feeagh water column, with the maximum St value of 669.41 J·m −2 reached on the 30th June. Schmidt stability remained well above the 15-year average throughout July 2018 (>510 J·m −2 ) ( Figure 3A). The mixed layer deepened from 6.5 m before the storm to 16 m on the 21st of June, highlighting a deepening of the mixed surface layer to a maximum of 9.5 m during the intermediate period ( Figure 3B). From July onwards, the mixed layer gradually deepened from 2 m on the 3rd July to 42 m by the end of September, when the lake became fully mixed again and there was a complete breakdown of the water column stratification ( Figure 3B). J·m −2 reached on the 30th June. Schmidt stability remained well above the 15-year average throughout July 2018 (>510 J·m −2 ) ( Figure 3A). The mixed layer deepened from 6.5 m before the storm to 16 m on the 21st of June, highlighting a deepening of the mixed surface layer to a maximum of 9.5 m during the intermediate period ( Figure 3B). From July onwards, the mixed layer gradually deepened from 2 m on the 3rd July to 42 m by the end of September, when the lake became fully mixed again and there was a complete breakdown of the water column stratification ( Figure 3B).   (Figure 4). Respiration increased while NEP decreased, corresponding to a decrease in GPP during the mixing event linked with extreme-weather-related conditions. The resulting NEP estimates were nearly always negative, but clearly higher during the periods with calm and hot sunny days (early June and early July) (Figure 4). days (n) for each specific period. Pre-and post-storm values of R and NEP were significantly different from the intermediate period (Kruskal-Wallis test and Dunn's post hoc test, p < 0.05) (Figure 4). Respiration increased while NEP decreased, corresponding to a decrease in GPP during the mixing event linked with extreme-weather-related conditions. The resulting NEP estimates were nearly always negative, but clearly higher during the periods with calm and hot sunny days (early June and early July) (Figure 4). The nutrient ratios (TN:TP) of Lough Feeagh and the inflowing streams ranged from 44 to 243 from April to October 2018 ( Figure 5A), and were strongly correlated (ρ = 0.64, p < 0.05, n = 14). However, a decoupling of nutrient ratios between the inflowing streams and Lough Feeagh was observed to begin to occur during Storm Hector in June, with TN:TP concentrations in Lough Feeagh much higher relative to TN:TP loading from the inflows ( Figure 5A). At the same time, nutrient imbalances between seston and zooplankton were positively different from 0, indicating potential P limitation for the zooplankton community during the pre-storm period ( Figure 5B). Later on, following the Storm Hector mixing event, there was a step change from P to N nutrient constraints for primary consumers, when nutrient imbalances were negative on the 10th July and remained closer to 0 from the end of July until October ( Figure 5B).
The nutrient ratios (TN:TP) of Lough Feeagh and the inflowing streams ranged from 44 to 243 from April to October 2018 ( Figure 5A), and were strongly correlated (ρ = 0.64, p < 0.05, n = 14). However, a decoupling of nutrient ratios between the inflowing streams and Lough Feeagh was observed to begin to occur during Storm Hector in June, with TN:TP concentrations in Lough Feeagh much higher relative to TN:TP loading from the inflows (Figure 5A). At the same time, nutrient imbalances between seston and zooplankton were positively different from 0, indicating potential P limitation for the zooplankton community during the pre-storm period ( Figure 5B). Later on, following the Storm Hector mixing event, there was a step change from P to N nutrient constraints for primary consumers, when nutrient imbalances were negative on the 10th July and remained closer to 0 from the end of July until October ( Figure 5B).

Phytoplankton and Ciliate Dynamics
Chlorophyll a concentrations in Lough Feeagh ranged from 0.76 to 2.61 µg Chl a·L −1 , Shannon diversity ranged from 1.02 to 2.66, and phytoplankton community biovolume ranged from 97.81 to 539.74 mm 3 ·m −3 during the study period (April to October 2018). Increases in phytoplankton biovolume occurred as soon as St increased above 400 J·m −2 , which was a consequence of the warm and calm meteorological conditions. In addition, from the end of May to mid-June (pre-storm period), Chl a, phytoplankton and ciliates biovolume quickly increased to maximum values of 1.93 µg Chl a·L −1 and 325.97 mm 3 ·m −3 , respectively. However, phytoplankton diversity decreased 2.5-fold during the first heatwave at the beginning of June ( Figure 6). From April to the end of May (mean ± SD, n = 5), the system was clearly dominated by Bacillariophyta (diatom bloom), with a biovolume of 134.48 ± 73.48 mm 3 ·m −3 (71.96 ± 13.86% of relative biovolume) in the top 6 m water column ( Figure 6A), but also at deeper layers (12 m and 25 m), with an average biovolume of 259.11 mm 3 ·m −3 (86.95%) for that specific day (Figure 7, 28th May pre-storm). On the 12th June, just before Storm Hector, and at the end of the first heat wave, Cryptophyta, Euglenophyta and ciliates increased their biovolume to 92.28 mm 3 ·m −3 (28.31%), 43.95 mm 3 ·m −3 (13.48%) and 69.08 mm 3 ·m −3 (21.19%), respectively. In turn, Bacillariophyta biovolume decreased to 105.09 mm 3 ·m −3 (32.24%). Specifically, there was a relatively large and rapid increase in Plagioselmis nannoplanctica (Crysophyta), Trachelomonas volvocina (Euglenophyta) and Strobilidium sp. (large ciliates), decreasing Shannon diversity index from 2.51 to 1.02 ( Figure 6A).
After the Storm Hector event (12 days post-storm), although there was a recovery in phytoplankton Shannon diversity to maximum values of 2.35, Chl a, phytoplankton and ciliate biovolume decreased, dropping to the minimum values for the study period (0.98 µg Chl a·L −1 and 171.1 mm 3 ·m −3 , respectively; Figure 6). By comparing phytoplankton community changes between the 12th (pre-storm period) and the 26th (intermediate period) June, it was observed that Cryptophyta and ciliates decreased to biovolumes of 24.72 mm 3 ·m −3 (14.45%) and 19.50 mm 3 ·m −3 (11.40%), respectively, while Chlorophyta biovolume started to increase to 44.00 mm 3 ·m −3 (25.71%) in the top 6 m layers ( Figure 6A).
From July onwards, the system restabilised and remained stratified until the autumn turnover period at the end of September (Figure 3). These repeated warm and calm conditions at a later stage of the phytoplankton seasonal dynamics resulted in another increase in Chl a concentrations, reaching a maximum of 2.61 µg Chl a·L −1 in August, but coinciding with another Shannon diversity minimum (1.22) on the 10th July (Figure 6), when the system was dominated by Plagioselmis nannoplantica in the Cryptophyta group and Botrycoccus brunni in the Chlorophyta group. Immediately following the restabilisation period (10th July), there was an increase in Chlorophyta biovolume to 113.03 mm 3 ·m −3 (33.39%) and, again, Cryptophyta biovolume increased to 110.30 mm 3 ·m −3 (32.58%). On the 10th July, total biovolumes reached similar values to pre-storm conditions (338.51 mm 3 ·m −3 ). In turn, on the 24th July, total biovolume reached maximum values of 535.49 mm 3 ·m −3 coinciding with Dinophyta (Gymnodinium sp. and Ceratium sp.) and cyanobacteria (Woronichinia naegeliana and Anabaena flos-aquae) biovolume increases to 113.38 mm 3 ·m −3 (21.17%) and 104.29 mm 3 ·m −3 (19.47%), respectively ( Figure 6A).
In contrast to what was observed on the 26th June, with similar total biovolumes at different depths (Figure 7, 26th June intermediate), on the 24th July (post-storm period), phytoplankton cells were more concentrated within the first top 6 m from the surface, representing a decrease from 535.49 to 254.63 mm 3 ·m −3 of total biovolume at 12 m and to 67.97 mm 3 ·m −3 at 25 m, respective to the 6 m integrated water sample taken on the same sampling day. Despite the lower biovolumes recorded at the bottom layers relative to the top 6 m layer, the bottom layers were highly dominated by ciliates and Cryptophyta taxa (Figure 7, 24th July post-storm). On the 7th August, cyanobacteria reached their maximum biovolume, corresponding to 210.46 mm 3 ·m −3 , an unprecedentedly high value for this lake [10], and this was mainly a consequence of large aggregates of Woronichinia naegeliana. In turn, from late July onwards, the Shannon diversity index increased again to a maximum value of 2.65 for the study period on the 3rd September ( Figure 6A).

Discussion
The work described here gives an insight into how a deep monomictic dystrophic lake situated on the Atlantic coast of Europe responded to sequential extreme weather events. To our knowledge, it is the first study reporting on planktonic community responses relating to the 2018 European heatwave, when Western Europe, in particular, experienced notably increased summer surface lake water temperature [52,64]. In turn, for the study site, Storm Hector was of particular interest, as it abruptly interrupted the heat wave conditions, with high temperatures and drought well outside the ranges normally recorded in Ireland. The physical structure of the Feeagh water column responded to the arrival of the first heat wave in the early summer season (late May-early June), with a rapid stratification of the water column (high Schmidt stability). Consequently, once Storm Hector arrived, Feeagh experienced a characteristic mixing event with a surface mixed layer deepening of 9.5 m at the end of June, never previously reported at this site [63]. Following Storm Hector, the anomalous pre-storm heat wave conditions returned, allowing an unprecedented isolated summer mixing episodic event in the middle of extremely high lake stability, for the site, to be examined.
The unusual lake stratification pattern presented in this study differed from another episodic flood event which occurred on the 2nd of July 2009 at the same site. Compared to the 2018 event, in 2009 there was a loss of Schmidt stability for the rest of the year and, despite the extreme nature of

Zooplankton Dynamics
From April to October 2018, total zooplankton biomass in Lough Feeagh ranged from 7.33 to 151.92 µg DW·L −1 ( Figure 6B). From mid-May to just before Storm Hector, total zooplankton biomass started to increase, coinciding with the initial peak in Chl a at this time. Nevertheless, total crustacean biomass remained below 20 µg DW·L −1 , with rotifers dominating zooplankton biomass ( Figure 6B).
In contrast, after Storm Hector at the end of June (intermediate period, 26th June), all zooplankton groups increased their biomass, achieving maximum values of 151.92 µg DW·L −1 . In particular, calanoids (mainly Eudiaptumus gracilis) and Daphnia increased their biomass to 46.52 µg DW·L −1 (30.62%) and 18.89 µg DW·L −1 (12.44%), respectively ( Figure 6B). Rotifers' total biomass also increased, with a direct effect on Keratella and Tricocherca groups (specific rotifers genus not shown). This peak in total zooplankton biomass coincided with a decrease in both phytoplankton and ciliates biovolumes, and Chl a concentration. (Figure 6).

Discussion
The work described here gives an insight into how a deep monomictic dystrophic lake situated on the Atlantic coast of Europe responded to sequential extreme weather events. To our knowledge, it is the first study reporting on planktonic community responses relating to the 2018 European heatwave, when Western Europe, in particular, experienced notably increased summer surface lake water temperature [52,64]. In turn, for the study site, Storm Hector was of particular interest, as it abruptly interrupted the heat wave conditions, with high temperatures and drought well outside the ranges normally recorded in Ireland. The physical structure of the Feeagh water column responded to the arrival of the first heat wave in the early summer season (late May-early June), with a rapid stratification of the water column (high Schmidt stability). Consequently, once Storm Hector arrived, Feeagh experienced a characteristic mixing event with a surface mixed layer deepening of 9.5 m at the end of June, never previously reported at this site [63]. Following Storm Hector, the anomalous pre-storm heat wave conditions returned, allowing an unprecedented isolated summer mixing episodic event in the middle of extremely high lake stability, for the site, to be examined.
The unusual lake stratification pattern presented in this study differed from another episodic flood event which occurred on the 2nd of July 2009 at the same site. Compared to the 2018 event, in 2009 there was a loss of Schmidt stability for the rest of the year and, despite the extreme nature of the event, the overall lake impact was considered relatively small [10]. In contrast, results from this study show that the mixing event induced by extreme continuous weather events drove a shift towards an increase in R and a decrease in NEP, despite the already natural heterotrophic conditions of the system [66]. In turn, a two-fold reduction in Chl a coincided with a three-fold reduction in lake Schmidt stability. In contrast, pre-and post-storm heat waves clearly encouraged fast growth of phytoplankton cells in the top water layers. As a consequence, an unexpected and atypical increase in cyanobacteria was observed at the end of July and beginning of August. In addition, observations of N:P imbalances for primary consumers in Lough Feeagh indicated a step change from P to N nutrient constrains following the mixing event in June. Furthermore, an increase in total zooplankton biomass during the intermediate period indicated a quick response to potential availability of more favourable food for primary consumers.
Metabolism estimates can provide a concise summation of the many effects of extreme meteorological events. In this study, we observed a significant increase in ecosystem R and a decrease in NEP during the intermediate mixing period, compared to pre-and post-storm warm and stable periods, which resulted in a slight decrease in GPP following Storm Hector. Similar to our study, some authors have also reported a decrease in GPP and an increase in R after flooding events, coinciding with an increase in bacterial biomass and a decrease in phytoplankton biomass and Chl a [10,44,46]. In contrast, results from thermocline-deepening experimental mesocosms, simulating the effects of storm events, showed that surface water GPP responded more strongly than R as a consequence of increased autotrophy driven by the proliferation of filamentous cyanobacteria released from light and nutrient limitation, resulting in an increase in NEP [43].
Phytoplankton primary production and bacterioplankton secondary production are influenced by bioavailable nutrients and, in turn, both processes support fluxes of energy and matter to higher trophic compartments [105]. In freshwater ecosystems, a range of organic and inorganic species of N and P co-occur. Both the amounts and stoichiometry of the total bioavailable resources determine nutrient limitation patterns in planktonic communities [23]. Taking into account that bulk TN:TP resources often vastly overestimates the corresponding bioavailable resources [105], in this study, we found similar in-lake and inflow river nutrient ratios, suggesting that lake water chemistry was highly influenced by nutrient loadings from the surrounding catchment, as observed in other studies [100]. In contrast, following Storm Hector, there was a decoupling period of TN:TP patterns of in-lake and inflowing streams lasting for 2 months (from mid-June to mid-August), with higher TN:TP values recorded in the lake compared with the inflowing rivers. Previous studies have concluded that, in general, TN:TP ratios are relatively higher in lakes than in their inflow streams, in part likely owing to enhanced water column primary production and microbial activity [106]. In addition, Sadro and Malack (2012) reported a disproportionally larger increase in nitrate loading during an extreme rain event, causing a substantial increase of in-lake dissolved inorganic nitrogen:total phosphorus ratio (DIN:TP) [46]. Another study suggested that storm events disproportionately mobilize P compared to N, contributing to reduction in N:P loading to the receiving water bodies [19]. Nevertheless, as Storm Hector in Feeagh mainly influenced the stratification pattern, instead of being characterised by an event that increased catchment nutrient runoff, the enhanced difference between inflow and in-lake total nutrient ratios seemed to be driven by internal biological dynamics accentuated, in turn, by the loss of lake stability, as suggested by Maranger et al. (2018) [106].
Phytoplankton community composition appeared to follow the natural seasonal succession expected for a dystrophic low-productivity lake such as Feeagh [65]. However, while cyanobacteria blooms have previously been linked to increases in water temperature and high stability [107][108][109][110], directly associated with extreme weather events [31,111], and associated with the depletion of available nitrogen [111][112][113][114], the higher cyanobacteria biomass that we found has rarely been seen before in Feeagh [10]. Woronichinia naegeliana appeared to take advantage of post-storm warm and calm conditions in Lough Feeagh. This cyanobacteria species is one of the most commonly associated with toxic cyanobacteria blooms [115]. It is not capable of fixing atmospheric nitrogen, but it is very efficient at taking up other dissolved fractions, with cells usually taking advantage when N concentrations are limiting to other algae [116]. A decrease in turbulent diffusion was shown to be the dominant environmental parameter determining the timing of phytoplankton development, together with an increase in water temperatures when light is limited (e.g., brown lakes such as Feeagh) [25]. Therefore, despite storm events disrupting the stability of the systems, if the persistent calm conditions of this study occur more frequently [4], a quick return to thermal stratification will potentially increase cyanobacteria biomass, even in lakes with low inorganic nutrient concentrations and high water colour [117].
In addition to cyanobacterial increases, a wide range of effects on phytoplankton community composition after physico-chemical lake changes driven by episodic meteorological events have been reported [118]. In accordance with our study, Havens et al. (2011) observed that, unlike previous documented effects of hurricanes, where increased N inputs stimulated primary productivity, a substantial reduction in phytoplankton biomass occurred, coinciding with a large decline in mean irradiance in the mixed layer and increases in phytoplankton grazers [55]. In this study, an increase in the biomass of grazers also simultaneously occurred with a decrease in phytoplankton biomass at different depths, particularly the Cryptophyta group, which is considered a high-quality food for primary consumers [119,120], and is highly selected by copepods [121]. In addition, ciliates, which can be a significant part of the copepod diet [122], also declined, which may have potentially led to a release of predation on bacteria resulting in an increase in ecosystem respiration [123]. Therefore, the increase in total zooplankton biomass during the intermediate period would support the idea that phytoplankton and ciliates experienced higher grazing pressure during the mixing event, in addition to a dilution effect at the top layers (above 6 m) and less favourable conditions for phytoplankton cells (decrease of euphotic depth and water temperature above 6 m).
Observations that calanoid copepods, in particular Eudiaptomus gracilis, experienced the highest increase in biomass compared with pre-storm conditions would agree with other published results [53,56,124]. This increase may be a result of higher-quality food availability due to higher turbulence enhancing predator-prey encounters [54,125] and the fact that mixed layer deepening brought food to lower depths allowing zooplankton to feed while still avoiding visual fish predators [126,127]. Moreover, despite a decrease in rotifers' relative biomass compared to the pre-storm period, rotifers total biomass also increased, as also reported in other studies [27,57]. In addition, an increase in temperature at deeper layers (below 6 m from the surface compared to pre-storm conditions) due to the mixing event after Storm Hector could have also played an important role in the increase of zooplankton biomass, owing to an increase in metabolic rates driven, at the same time, by an increase in water temperature [50].
After Storm Hector, two peaks in copepod abundance potentially accounted for the shift from P to N zooplankton nutrient constraints during July and September. N:P imbalances were negative during both occasions, suggesting that zooplankton N:P body content was higher than seston N:P, which can be accounted for by the dominance of copepods in the system, leading to higher seston N uptake [61,128]. This is in accordance with ecological stoichiometric theory, which states that taxa such as copepods require more N in their diet than high P demanding taxa such as Daphnia [23]. The N:P imbalance remained close to 0 (no nutrient imbalances) for the remainder of the study period (August-October), suggesting a potential nutrient constraint release for primary consumers, driven by the weather induced mixing event, and maintained for the rest of the summer period.
In contrast, zooplankton biomass declined once the lake restabilised, despite this coinciding with a peak in Chl a concentration. These observations can be explained by the fact that phytoplankton community composition at that time was dominated by less nutritious groups (e.g., large Chlorophyta and cyanobacteria), known to be low in essential fatty acids [129,130], which were also less available (e.g., concentrated at top 6 m water layers) [131]. Therefore, a decrease in zooplankton biomass when Chl a values were maximum would suggest that extremely stable conditions are not linked to higher food availability for primary consumers, as well as a potential higher risk of fish predation at the top layers where food was in higher quantities. Minimum zooplankton biomass values observed during April could be linked with potentially higher fish predation as a consequence of the movement of salmon (Salmo salar) smolts downstream to the sea at that time of the year [132].

Conclusions
More frequent summer heat waves, leading to highly stable lake conditions, could enhance atypical cyanobacteria increases, even in nutrient-and light-limited aquatic ecosystems (e.g., dystrophic lakes). In addition, this study suggests that lake mixing disturbances driven by storms, especially in lakes experiencing more intense summer stratification due to warming, can enhance ecosystem respiration and primary consumer biomass and, in turn, reduce zooplankton nutrient constraints and phytoplankton biomass. Future scenarios predict more summers such as the one presented in this paper. Characterisations of both lake thermal, nutrient and ecological responses to consecutive extreme weather events are relatively rare, but are crucial to understand how lakes are changing, inform forecasting tools and therefore protect such ecosystems, as the impacts of global climate change accelerate.