How Might Changing Climate Limit Cyanobacteria Growth in Shallow Prairie Lakes? An Empirical Space-For-Time Evaluation of the Potential Role of Increasing Sulfate

Cyanobacteria blooms alter aquatic ecosystems and occur frequently in shallow prairie lakes, which are predicted to increase in salinity as the regional climate becomes hotter and drier. However, flat landscapes that experience depression bottom salinity with high concentrations of sulfate in addition to sodium and chloride, may mitigate nutrient increases or even inhibit cyanobacteria growth. Cyanobacteria can dominate shallow lakes with low N:P ratios because many cyanobacteria species fix dissolved N2, whether due to in-lake denitrification or exchange with the atmosphere, a process that requires molybdenum as an enzyme cofactor. Sulfate can compete with molybdate at cellular uptake sites, potentially limiting the competitive advantage of cyanobacteria. We studied 25 lakes located in a relatively limited geographic region of southern Alberta (Canada) and used a space-for-time analysis to model scenarios of increased sulfate concentrations under changing climate. Monthly, we measured nitrogen, phosphorus, sulfate, molybdenum, and cyanobacterial pigments and used mixed effects models to identify empirical relationships. Sulfate drives conductivity in the region and we found that most saline lakes we sampled are turbid lakes with high nutrients and high cyanobacteria biomass. In addition to phosphorus, molybdenum predicted cyanobacterial pigments in the top two models, showing a positive relationship with cyanobacterial biomass. We also found a negative relationship between sulfate concentration and cyanobacteria pigments, which suggests that as lakes get saltier, even with increased nutrients, there may not be an incremental increase in cyanobacteria biomass. Our results therefore suggest that competition between sulfate and molybdate may limit future cyanobacteria growth in shallow lakes and that with a warmer and drier climate it may not be inevitable that shallow lakes will continue to be dominated by cyanobacterial blooms, a hypothesis that could be tested directly via experimentation.


Introduction
Toxic cyanobacteria are a global public health concern [1,2] and are increasingly found in shallow lakes [3][4][5][6]. Many shallow prairie lakes are endorheic basins that are expected to increase in salinity due to regional climate becoming hotter and drier, which would decrease water inputs, increase water evaporation and decrease water levels, and therefore increase salt concentrations [7][8][9]. Salts, such as sulfate, accumulate from weathering of surrounding soil and bedrock. Sulfate concentrations in shallow lakes may increase [10] as prairie climate becomes hotter (increases of 3-5°C annual mean temperatures) and drier (decreases of 5-10% in precipitation) by the 2050s [9]. In agricultural regions, shallow prairie lakes also have increased nitrogen (N) and phosphorus (P) due to fertilizer used on surrounding croplands over the past 60-70 years [11,12]. The role of macronutrients and their supply ratios in generating cyanobacterial blooms has been well studied [13,14], but the importance of micronutrients, including enzyme co-factors, is much less clear, particularly in shallow lakes. However, increases in salinity may actually inhibit cyanobacteria growth [15].
Phosphorus has been argued to be the major driver of cyanobacteria growth [16,17] with the role of N more contentious [18,19]. Phytoplankton biomass often correlates with total P [20], yet a low N:P ratio can shift the phytoplankton to cyanobacterial dominance [21] because cyanobacterial N-fixation can reduce N-limitation [22]. However, if N-fixation alone provided enough N, primary productivity should never be N-limited [23]. Low productivity biomes, such as open oceans, are in contact with the N 2 -rich atmosphere yet make it clear that N-limitation exists even with abundant N 2 [24].
During times of N-limitation, N-fixation provides a competitive advantage to N-fixing cyanobacteria provided there are sufficient micronutrients, such as molybdenum (Mo) or iron (Fe) that are enzyme cofactors of nitrogenase [25]. Mo is a trace metal predominantly found in oxic water as molybdate (MoO 4 2-) [26] and derives from geologic weathering [27]. Typically, MoO 4 2-concentrations are relatively conserved [28]. Cyanobacteria take up MoO 4 2through specific sites prior to synthesizing nitrogenase [26]. However, in saline water, sulfate (SO 4 2-) can outcompete MoO 4 2at uptake sites because the ions share a similar charge-to-mass ratio and stereochemistry [26]. Although low Mo availability has been shown to limit cyanobacteria growth in saline coastal systems [29], it remains unknown whether this applies to all aquatic systems or whether a threshold exists where SO 4 2restricts the availability of MoO 4 2enough to limit cyanobacteria growth in freshwater systems, including shallow prairie lakes.
The goal of this study was to investigate implications of changing climate on future cyanobacteria growth in shallow prairie lakes. We employed a space-for-time analysis with 25 southern Alberta shallow lakes that span a large range in salinity, yet are within a relatively limited geographical area. We used general linear mixed effects models to test the hypothesis that MoO 4 2availability limits cyanobacteria growth and therefore late summer cyanobacteria biomass. We also looked for empirical evidence of a threshold ratio where SO 4 2outcompetes MoO 4 2and inhibits cyanobacteria growth.

Study Area
Lakes sampled were 60 to 160 kilometers east of Calgary, Alberta, in the semi-arid prairie pothole region characterized by mixed fescue grasses, black and brown chernozemic and solonetzic soils where there are natural saline conditions [30]. The dominant salts are sodium and magnesium sulfates, and salts primarily derive from bedrock and glacial till [31]. Groundwater seepages that concentrate salts are primarily local, rather than regional, and primarily depression-bottom salinity [31]. A variety of crops, such as spring and winter (non-durum) wheat and canola dominate the land cover [32] in addition to livestock operations. In 2011, 88% of total cropland in the local area had commercial fertilizer applied to it [33]. Livestock manure is typically redistributed on fields [33].

Lake Selection
Sampled lakes ( Figure 1) were purposely chosen to include a salinity gradient, which was estimated from conductivity measured in the field before sampling began. Eleven lakes had been previously sampled [34,35], which provided a known range of salinity and history of cyanobacteria blooms for those lakes. Table 1 summarizes selected physical and chemical attributes of all study lakes.  Table 1 Selected characteristics of 25 shallow prairie lakes sampled to determine relationships between cyanobacteria pigments, sulfate and molybdate across a broad salinity gradient. Values are the mean +/-SEM of three samples taken in August 2016. Lake location and selected biological and chemical characteristics of 25 shallow prairie lakes sampled in Alberta. *Z max was taken from previous studies [34,35]. Chlorophyll-a Cyanopigment (mg · L -1 ) (mg· L -1 ) TP (mg · L -1 ) (mg · L -1 ) (mg · L -1 ) (mg · L -1 ) Two Bar 7,898 (+/-83) 1,920 (+/-527)

Lake Sampling
Twenty-five lakes were sampled four times (once each in June, July, August and September) between June 13 and September 6, 2016. At each lake, three samples were taken in acid-bathed (minimum 3 hours in 25% HCl and triple DDW rinsed) 1L Nalgene bottles by wading 1 to 4 meters from the shore and by using a bottle holder to avoid debris suspended from wading. At each site, the bottle was triple-rinsed with lake water and then the sample was taken from immediately under the water surface, then stored in a dark cooler on ice for transport back to the lab. The depth at the sampling site never exceeded 1.5 m and occasionally was less than 30 cm. The samples were purposely taken to include a "worst-case" scenario by visually inspecting the lake's surface for suspected cyanobacteria abundance, and sampling in the bloom, if seen. The three sampling sites were at least 30 m apart, if possible. Turbidity was measured with a HACH Portable turbidity meter (model 2100Q). Conductivity and temperature were measured with a ThermoScientific Orion Star (model A325) pH/conductivity portable multi-parameter meter. The pH of each sample was measured in the lab using a Mettler Toledo, FiveEasy Plus (model FP20) pH meter.

Ion and Nutrient Analysis of Lake Water
Soluble reactive phosphorus, ammonia, nitrate + nitrite, sulfate and molybdenum samples were filtered using Pall Corporation Life Sciences GN Metricel Grid 0.45 m filters within 7 hours of collection and stored in sterile 50 mL polypropylene conical tubes overnight at 4C. All nutrient samples were analyzed within 24 hours of collection, except where noted.
Total nitrogen was measured on a Shimadzu TOC-L Combustion Analyzer with TNM-L module with a Shimadzu ASI-L auto sampler with dilution, if necessary. Nitrate and nitrite were measured using ion chromatography on a Metrohm 940 Professional IC Vario equipped with a Metrohm 858 Professional sample processor.
Ammonia, sulfate, total phosphorus and soluble reactive phosphorus were measured on a WestCo Scientific Instruments Inc. Smart Chem (Model 170) discrete analyzer. Ammonia was measured with the Berthelot reaction following Method AMM-001-A [36]. Total phosphorus and soluble reactive phosphorus were measured with the molybdate blue method [37] using Method PHO-004-A [38]. Sulfate was measured using Method SUL-002-A [39]. Samples were auto-diluted by the Smart Chem or manually diluted beforehand and auto-diluted when sulfate levels were above the calibration curve.
Molybdenum samples were filtered through a Pall Corporation Life Sciences GN Metricel Grid 0.45 m and refrigerated at 4°C for a maximum of 5 months. At analysis, samples were acidified to 1% HNO 3to prevent precipitation and measured using an Agilent Technologies 8800 ICP-MS Triple Quad with an ASX-500 Series ICP-MS auto sampler. Molybdenum 95 and 98 were measured and the molybdenum 98 isotope measurements were used. Typically, triplicate measurements were made from a single vial for a given sample. In five lakes (Brush, Black, Whey, Horse and Dawson), the salt concentrations created matrix effects and in these instances, the method of standard addition was performed to obtain an accurate measurement. Indium was used as the internal standard to correct for any signal drift during analysis. Total molybdenum (Mo) was measured because molybdenum is predominantly found in oxic waters as the thermodynamically stable oxyanion molybdate [26], which is the Mo form taken up by cyanobacteria.

Pigment Analyses
Chlorophyll-a and accessory pigment measurements were performed on phytoplankton collected on 47 mm VWR glass microfiber filters (model 696; 1.2 m pore size) under low light and then frozen in aluminum packets at -20 o C for up to 5 months for later analyses. Filters were first freeze-dried in the dark over 48 hours in a Labconco, FreeZone 6 freeze drier to remove remaining water, then the pigments were extracted under low light by placing the filter in 10 mL of 98% ethanol, vortexing and letting the samples extract for 24 hours. Samples were then centrifuged and measured in small batches using a Molecular Devices SpectraMax M2 plate reader at 300 wavelengths from 400 -700 nm. The data were then analyzed using pigment-based chemotaxonomy and a Gaussian peak function, which predicts the range of pigments present in the sample [40,41]. The software estimates pigments from cyanobacteria, diatoms, dinoflagellates and green algae. We used chlorophyll-a, plus myxoxanthin, canthaxanthin and echinenone, the latter three which are specific to cyanobacteria, to create a parameter called "cyanopigment" by selecting the highest concentration out of the three pigments from each sample. Adding the three cyanobacteria pigments together would over-estimate the concentrations because all three of the pigments can occur in the same cyanobacterial species [42]. However, to use only one pigment would underestimate the cyanobacteria biomass because not every species produces all three pigments. Use of the highest pigment still underestimates the total concentration if two species in a sample contain different accessory pigments and we only use the highest accessory pigment to quantify the sample.

Statistical Analyses
All analyses were performed with R version 3.3.2 [43], and RStudio version 1.0.136 [44], using the lme4 [45], lmerTest [46] and effects [47] packages. General linear mixed effects models were analyzed with the lmer function in the lme4 R package [45] and generalized linear mixed effects models were analyzed with the lme4 function in the glmer package [48]. Stepwise regression was performed with the lme4::drop1 function [45]. P-values were calculated with the lme4::drop1 function using the Chi-squared test and with the lmerTest package using the Satterthwaite method. Figures were made with MS Excel, ggplot2 [49] and ggmap [50].
Cyanopigment was the dependent variable of three models: one with three extreme lakes removed (lme), one with all lakes included (lme) and one with the binomial family (glmm). For the two linear mixed effects models, lake and month were included as random effects. The random effect structure reflected the study design by allowing a random slope of month, within lake, with correlated intercept [51] and also calculated covariance. The assumptions of normality and homoscedasticity in the residuals were met by visually inspecting the residuals vs. fitted plot [52]. Only observations of cyanopigment above the detectable limits were included in the linear mixed effects models. To minimize the spread in the fixed effects' ranges and to reduce heterogeneity, sulfate was log 10 transformed and the other fixed effects and the dependent variable were natural log (ln) transformed. To test the effect of extreme values on the models, we ran the model with and without the values and report both models. In the first model, two lakes with extreme sulfate concentrations were removed, which lowered the range of sulfate from 13-17,324 mg · L -1 to 13-7,856 mg · L -1 (mean: 1460, median: 854) and one lake with extreme molybdenum concentrations was also removed, which lowered the range from 0.4-132 g · L -1 to 0.4-38 g · L -1 . Truncated data were included in the third model where cyanopigment presence/absence was coded with a categorical variable (0s, 1s) to determine significant effects for detection and below detection. The binomial model was a generalized linear mixed effects model with the family as binomial and the link function as logit. This model included all lakes and tested the same fixed effects as the first model; however the interaction term was not included. In this model, lake was included as a random effect and month was included as a fixed effect. Month was not included as a random effect on its own because there were insufficient levels (five to six levels are recommended at minimum; [51]). Variables were scaled to account for covariance.

General Patterns of the Study Lakes
The 25 shallow lakes sampled included broad ranges of conductivity, sulfate, total nitrogen and total phosphorus while in contrast, molybdenum varied less (Table 1). Conductivity had a relatively even and broad range (330.9-20,760 µS/cm - Figure 2) as lakes were purposely selected to include a conductivity gradient. Sulfate encompassed three orders of magnitude (13-17,324 mg · L -1 - Figure  2) with higher sulfate levels further east (Figure 1), except for the two large lakes farthest east (Seiu Lake and Crawling Valley). Brush Lake (Figure 2) had the highest conductivity (18.17-20.7 mS/cm) yet almost the lowest sulfate (13-93 mg . L -1 ). Molybdenum concentrations varied less across the 25 lakes (median 3.0 g · L -1 ; range 0.4-132 g · L -1 ). Goat Lake had a maximum molybdenum concentration 2.5 times greater than the next highest lake (132 g · L -1 compared to 37 g · L -1 ), which accounts for the high range and mean. Every lake had detectable molybdenum concentrations.

Figure 2
Relationship between conductivity (mS . cm -1 ), sulfate (mg . L -1 ) and cyanopigment (µg . L -1 ) for 25 shallow prairie lakes sampled in July 2016. Points represent the mean of three samples. The horizontal black dashed line represents 2 mS . cm -1 , the approximate cut off between fresh and brackish water. The dotted line represents the regression line. The full range (a) and (b) details for conductivity < 14 mS . cm -1 and sulfate < 7,000 mg . L -1 are shown.
Three lakes had low conductivity (<0.800 mS . cm -1 ), low sulfate and low cyanobacteria biomass ( Figure 2). These lakes also had low N and P, and macrophytes were visually observed during sampling. Another group of lakes between ~1.0-2.0 mS . cm -1 (170 -730 mg·L -1 sulfate) consistently had higher cyanobacteria pigment (>1 µg·L -1 ), higher nutrients and higher sulfate. Previous studies have found that sulfate concentrations of 768 mg·L -1 (8 mM) and greater start to inhibit molybdate uptake [54]. Some lakes above 2.0 mS . cm -1 consistently had very low or no cyanobacterial pigments; however a few lakes had mean cyanobacterial pigments between 1-5 µg · L -1 and two saline lakes had mean cyanobacterial pigments above 14 µg · L -1 .
All lakes were relatively small and shallow, except for Crawling Valley, which is a large reservoir with riverine characteristics. The lakes were slightly alkaline (pH 8.09-10.68) and the water temperature ranged from 11.5-26.0 °C. Turbidity ranged from 1.03-800 NTU. Although not confirmed in this study, some of the shallow lakes (e.g. Mushroom) have historically been fishless ( [34]; Jackson unpublished data, August 2018). Barnett Lake has Brook stickleback (Culaea inconstans), West Lake has Prussian Carp (Carassius gibelio) and presumably East Lake and Long Lake do too, based on their proximity and connecting culverts (Jackson, unpublished  Joshi & Jackson, Figure 3 pike (Esox lucius) [55]. Crawling Valley has been stocked with sport fish. Three species of sucker and two species of minnows have also been caught in the reservoir [55].

Patterns of Cyanobacteria Pigments with Sulfate and Molybdenum
In general, we found higher cyanobacteria pigments (our proxy for biomass) in lakes with TN:TP < 16:1 (Figure 4). Lakes with TN:TP below 16:1 typically also had low sulfate:molybdenum. The first model found that ln cyanopigment was significantly correlated to ln(TN), ln(TP) in July and September, and ln(Mo). Log 10 (SO 4 2-), although not significant (p = 0.055), was included in the model. The linear mixed effects model used 22 lakes; the two lakes with the highest sulfate and one lake with the highest molybdenum (and greatest variability in the sulfate and molybdedum measurements) were removed from this model and only samples with cyanobacteria pigment >0 were used (n = 165). The model found that ln(cyanobacteria pigment) was not significantly predicted by log 10 (SO 4 2-) (p = 0.055), with a negative relationship that had an estimate of -0.89 +/-0.45 (intercept was 1.14 [52]). Although this relationship was not significant at α = 0.05 using the Satterthwaite method, it was closer, yet still not significant using the chi-square test (p = 0.051). There was a significant positive relationship between ln(cyanobacteria pigments) and ln(Mo) (p = 0.005), that had an estimate of 0.81 +/-0.29, and a strong, significant relationship between ln(cyanobacteria pigments) and ln(TN) (p = 0.005), with an estimate of 1.18 +/-0.41 (Supplementary Figure S1). There was also a significant interaction by month between ln(cyanobacteria pigments) and ln(TP), with a significant relationship in July  Model 3 found that cyanopigment presence was significantly predicted by ln(TP). To account for the dependent variable being truncated in Models 1 and 2, a generalized linear mixed effects model using the binomial distribution to analyze detection and non-detection in cyanobacteria pigment was analyzed to complement the previous two models.
All lakes and samples were included in the model (n = 300), which found ln total phosphorus could significantly predict the presence (detection) of cyanobacterial pigments (p = 0.02) with an estimate of 0.41 +/-0.18 (intercept was 1.44 [52]). The other fixed effects in Model 1 were tested in this model but were not significant: scaled log 10 of sulfate (p = 0.79), scaled ln molybdenum (p = 54) and scaled ln total nitrogen (p = 31). These terms were removed from the model.

Discussion
We found a significant positive relationship between cyanobacteria pigment and molybdenum concentration and negative coefficients between cyanobacteria pigment (our proxy for cyanobacteria biomass; [56]) and sulfate concentration. These relationships suggest that increasing sulfate and decreasing molybdenum may limit cyanobacterial growth by out-competing molybdate uptake needed for the Mo-nitrogenase enzyme system responsible for biological N-fixation [57]. Because sulfate and conductivity were strongly correlated in the shallow prairie lakes sampled and regional climate is predicted to become hotter and drier, sulfate will increase in these shallow endorheic lakes and could limit cyanobacterial growth that is projected to increase over the next 30-60 years from anthropogenic eutrophication.
N-fixing cyanobacteria have been shown to dominate low N:P lakes through their ability to fix atmospherically derived N 2 because other phytoplankton groups cannot [16,21,22]. Some, but not all, lakes had higher cyanobacteria growth at low N:P ratios (Figure 4), and the low N:P lakes with the highest cyanobacterial pigments also had low sulfate:molybdenum ratios, which supports the notion that when molybdate availability does not limit Mo-nitrogenase production, cyanobacteria can dominate in low N:P lakes. However, some studies suggest that the total amounts of N and P better predict cyanobacterial dominance [58,59]. Total nitrogen can be a proxy for phytoplankton biomass, so a strong relationship between cyanobacterial pigment and total nitrogen is not surprising, yet shows that the amount of total nitrogen had the strongest and most consistent correlation to the biomass of the cyanobacterial community. The effect of total phosphorus depended on the month, indicating either a seasonal response to nutrients, or an unmeasured process, such as zooplankton grazing, that affected cyanobacterial biomass. Our results provide correlations between cyanopigment, SO 4 2and Mo (our estimate of MoO 4 2-); experimental manipulations would be required to confirm causality.
High concentrations of N and P come primarily from fertilizers applied to the surrounding cropland [11] although this region also has high geologic P inputs from glacial till. Most synthetic nitrogen applied to crops is used by the plants, yet isotopic tracing reveals that nitrogen can have long residency times in the soil and that 8-12% of nitrogen can enter aquatic ecosystems and groundwater [11]. Nitrogen would leave an endorheic shallow lake water column mainly through denitrification and sedimentation; however, inputs tend to exceed losses, so nitrogen and phosphorus concentrations tend to accumulate in shallow lakes. The concentrations of total nitrogen and total phosphorus in the lakes we sampled (Table 1) even exceeded measurements in other shallow lakes in Europe or Asia [3,60], yet the concentrations we measured are similar to previous measurements in these lakes [35]. Most of the lakes range between eutrophic and highly eutrophic and therefore have favourable nutrient conditions for cyanobacterial blooms. However as total nitrogen strongly correlated with conductivity ( Figure 3) these highly eutrophic lakes are also more likely to have high sulfate concentrations and provided an opportunity to evaluate the interacting effects of N and sulfate.
Cyanobacteria require molybdenum as an essential cofactor for the predominant form of nitrogenase [61]. Microscopy on phytoplankton communities in many of our study lakes over the last 20 years reveals that the same cyanobacterial species appear in summer blooms in these lakes. Microcystis sp. typically appears in late June or July. As summer progresses, the Cyanobacterial community is typically dominated (numerically) by N-fixing species such as Aphanizomenon flosaquae, Oscillatoria sp., Anabaena sp., Nostoc sp., Lyngbya sp., Gleotrichia sp., and Gleocapsa sp.. Blooms of N-fixing capable cyanobacteria suggests that the lakes have transitioned to N-limitation, a hypothesis that could be tested by N addition to mesocosms. Molybdenum predicted cyanobacterial pigments in both models and through all months, which supported the hypothesis that molybdate is correlated to cyanobacteria biomass. Sulfate competes with molybdate for uptake sites and showed a negative (though non-significant; p = 0.055) relationship with cyanobacterial growth. Although sulfate was not statistically significant in the model, the borderline non-significant statistical result may reveal a biologically significant relationship (additional samples from other lakes would help clarify the generality of this finding). Sulfate would derive from localized areas in the region with high salinity where salts accumulate in groundwater and then pool at the surface in patchy distributions, making some lakes highly saline and other lakes fresh [31]. Nitrogen is also assumed to come from leaching through soil and groundwater, and the combined results of sulfate, molybdenum, N and P together show that currently shallow lakes do produce large cyanobacterial blooms, some with associated toxins [52]. Although sulfate can limit cyanobacterial growth, not every lake above 2.0 mS/cm (~750 mg·L -1 sulfate) had low cyanobacterial biomass, presumably because these lakes are eutrophic and have N from sources other than the atmosphere (e.g. NH 3 , urea, dissolved organic nitrogen). Our results are consistent with previous studies that have shown that sulfate and molybdate affect nitrogen fixation in estuaries [25,62,63] where rates of N-fixation are low because molybdenum is scarce due to high sulfate [25].
Seawater is about 29 mM SO 4 2- [64] and about 107 nM Mo [65], which is approximately 270,000:1 SO 4 2-:Mo. In the study lakes we calculate an average ratio about one order of magnitude higher at 1.308 x 10 6 even if we assume Mo has the molecular mass of MoO 4 2because molybdate is the predominant form of Mo. The ratio of SO 4 2-:MoO 4 2therefore is higher in these shallow saline lakes than in seawater. Cole [64] suggested that, in marine systems, sulfate would inhibit MoO 4 2uptake at about 5% seawater concentrations or about 1.45 mM sulfate. Our lakes average (+/-1 sd) 17.65 (20.8) mM sulfate. If cyanobacteria are similarly affected by sulfate in marine and freshwater systems, sulfate would limit Mo availability in all but four lakes and would be right at limitation in one lake. This could be tested experimentally by adding eg, Na 2 MoO 4 to lower the SO 4 2-:MoO 4 2and measuring the corresponding response of cyanobacteria.
Mo has historically been thought to be essential for N 2 fixation [66]; however, recent evidence indicates that some terrestrial N-fixers have V-and Fe-only nitrogenases [67]. Together with the Mo-nitrogenase, these forms are collectively encoded by a complex system of over 80 [61] Nif genes that produce the necessary gene products for fully functional nitrogenases. In some laboratory growth studies, N-fixation has occurred despite no free Mo [68,69] and many genes responsible for V-and Fe-nitrogenases have now been identified [70,71]. Whether or not the same V-and Fenitrogenase forms identified in soil bacteria exist in freshwater cyanobacteria is not known; however, if they are present and produced they could potentially lead to weaker correlations between Mo and N-fixing cyanobacteria in aquatic systems and could provide a mechanism for a compensatory response if increasing sulfate effectively makes molybdate less available. Furthermore, Fe form and therefore its availability has also been shown to create Fe-limitation in some systems, including saline prairie lakes [72]. In shallow, well-mixed polymictic lakes, such as those we've sampled, Fe 2+ released from anoxic sediments would be rapidly oxidized but cyanobacteria do produce siderophones to aid in Fe acquisition [73].
Some of the high TN, TP, chlorophyll-a and cyanopigments we measured would be the result of the worst-case scenario sampling we employed. There would be different values of TN and TP from an integrated sample taken at various depths because of patchiness in the algal biomass in some lakes. The measurements would also be affected by variation in annual phytoplankton phenologyspring (2016) was about 3 weeks earlier than 'normal' with about three times higher precipitation in July 2016 (142 mm) than the monthly average for July (49 mm from 2004-2015) [74]. Even though the September sampling was at the very beginning of the month, most of the cyanobacterial blooms seen in July and August were not visible by September. The lakes were located within 100 km of each other and should share regional climate and geology. The models we ran indicate that even when interactions between independent variables are included, the bivariate relationships remain significant. Furthermore, the removal of outliers from the data in Model 1 did not lead to different results when compared to Model 2.

Thresholds and Alternate Stable States in the Lakes
Shallow lakes and aquatic ecosystems are among the most altered ecosystems on the planet [75] and continually respond to anthropogenic nutrient loading. Nutrient inputs can shift lakes from a clear, low-nutrient, high-macrophyte state to a turbid, high-nutrient, high-phytoplankton state, the latter often containing cyanobacteria [34,[76][77][78]. We did not identify statistical thresholds that would suggest bimodality in the relationships between sulfate or molybdenum and cyanobacteria and therefore the relationships between sulfate, molybdenum, N, P and cyanobacteria are continuous in the lakes and year we sampled.
Despite lack of evidence of nutrient-cyanobacteria thresholds, the lakes can still be grouped into three categories. In the two lakes with low nutrients, low salinity and low cyanobacterial biomass ( Figure 2) we also noted high macrophyte abundance and low turbidity. A second group of lakes had higher nutrients, were more turbid and had more cyanobacteria biomass, and while they had increased conductivity, it was below 2.0 mS/cm, which is near the threshold of molybdate uptake inhibition of 768 mg·L -1 sulfate found by Marino et al. [54] and also a rough distinction between fresh and brackish water. The third group of lakes had conductivity above 2.0 mS/cm and also had increasing sulfate, increasing turbidity and less cyanobacteria biomass. However, this last group of lakes also had the highest N ( Figure 3). There were two lakes above this threshold where cyanobacteria bloomed above 14 µg·L -1 , showing that this threshold did not entirely limit cyanobacterial growth. One of these lakes (Dog, which had extreme blooms) may have an unidentified nutrient source while the other lake (Bow) has abundant Microcystis sp., which does not fix N. Other brackish lakes had cyanobacteria pigment between 1-5 µg · L -1 periodically through the season. Cyanobacterial growth in brackish lakes shows that other conditions can override sulfate inhibition; for example, high N concentrations can support cyanobacteria growth that does not require N-fixation.
Most of the lakes we studied appear to be in a turbid state. A level of ~4 NTU has been identified in these shallow lakes as a threshold between clear and turbid states [34]. Three lakes were consistently below this level, while 22 lakes had values higher than this threshold. This strongly suggests that the majority of saline lakes in this region are turbid lakes with high nutrients and higher phytoplankton growth. An increase in turbidity leads to a decline in macrophytes as lakes become more eutrophied [79]; a survey in Europe showed that macrophytes did not grow above 2.5 mg/L of TN (Jeppesen et al. 2007). This turbidity has implications for management because clear lakes are typically more desirable, and are associated with higher water quality. Turbid lakes can also be challenging for management because shallow lakes tend to resist returning to a clear state even when nutrient levels are reduced [80], or nutrient levels need to be reduced lower than the level that tipped the lake into the turbid state originally. Furthermore, turbid lakes might exist as either intermediate turbidity, mixed-assemblage phytoplankton communities or high turbidity, cyanobacteria-dominated communities [80]. Cyanobacterial blooms can also produce microcystins, which pose additional challenges for drinking water, recreational use, watering livestock and management.
Zooplankton help to stabilize clear lakes by grazing phytoplankton; however, zooplankton have historically been thought to be unable to handle large colonial cyanobacteria and lack critical nutrients from grazing cyanobacteria [81][82][83]. More recent research suggests that zooplankton can graze large filamentous and colonial cyanobacteria [84] and that successive generations of zooplankton living in a cyanobacteria-dominated lake could co-exist with, but not control, increased cyanobacterial blooms [85]. With increasing salinity, the species richness in the zooplankton community may decrease, however, if salinity increases to a level where fish cannot survive, zooplankton may positively respond to a lack of predators [86]. We did not measure zooplankton abundance, yet their grazing could affect our results because it would reduce the phytoplankton biomass and specifically cyanobacteria biomass we measured, which of course would affect relationships with other variables. This unmeasured factor could also affect the significance of interaction terms by month because the zooplankton populations' predation pressure may rise and fall over the summer months, as evidenced by zooplankton-related phenomena such as the spring clear-water phase [87].

Conclusions
The negative relationship between sulfate and cyanobacteria biomass in our space-for-time analysis suggests that as prairie lakes become more saline, even with increased nutrients, cyanobacteria biomass will decrease. Climate change is predicted to increase annual mean temperatures across the prairies by 3-5°C and reduce precipitation by 5-10% by the 2050s [9]. While warm temperatures can favour cyanobacterial growth [88], warm temperatures can decrease water levels in evaporation basins, increase sulfate concentrations and potentially limit cyanobacteria growth. Knowledge gained from this space-for-time analysis suggests that increasing sulfate affects cyanobacteria growth, although not all lakes with high sulfate will have low cyanobacteria.