Abstract
Glacier algal blooms dominate the surfaces of glaciers and ice sheets during summer melt seasons, with larger blooms anticipated in years that experience the greatest melt. Here, we characterize the glacier algal bloom proliferating on Morteratsch glacier, Switzerland, during the record 2022 melt season, when the Swiss Alps lost three times more ice than the decadal average. Glacier algal cellular abundance (cells ml−1), biovolume (μm3 cell−1), photophysiology (Fv/Fm, rETRmax), and stoichiometry (C:N ratios) were constrained across three elevations on Morteratsch glacier during late August 2022 and compared with measurements of aqueous geochemistry and outputs of nutrient spiking experiments. While a substantial glacier algal bloom was apparent during summer 2022, abundances ranged from 1.78 × 104 to 8.95 × 105 cells ml−1 of meltwater and did not scale linearly with the magnitude of the 2022 melt season. Instead, spatiotemporal heterogeneity in algal distribution across Morteratsch glacier leads us to propose melt-water-redistribution of (larger) glacier algal cells down-glacier and presumptive export of cells from the system as an important mechanism to set overall bloom carrying capacity on steep valley glaciers during high melt years. Despite the paradox of abundant glacier algae within seemingly oligotrophic surface ice, we found no evidence for inorganic nutrient limitation as an important bottom-up control within our study site, supporting our hypothesis above. Fundamental physical constraints may thus cap bloom carrying-capacities on valley glaciers as 21st century melting continues.
Introduction
Streptophyte glacier algae (Procházková et al., ) represent the dominant primary producers in glacier surface ice environments, forming large-scale blooms during summer melt seasons that can approach densities of 103-105 cells per ml of meltwater under optimal growth conditions (Yallop et al., ; Williamson et al., ; Di Mauro et al., ; Irvine-Fynn et al., ). At these densities, glacier algal blooms lend a distinctive brownish/purple colouration to the ice surface because of the cells' heavy investment in secondary phenolic pigmentation (Remias et al., ; Williamson et al., ); shown to be a key adaptation to life in surface ice that protects against high irradiance (Williamson et al., ). Given the substantial ice surface darkening and resulting albedo decline associated with the growth of these heavily pigmented microalgal cells (Yallop et al., ; Di Mauro et al., ; Halbach et al., ), glacial algal blooms hold significant potential to exacerbate the melting of glacial ice across the cryosphere: in Greenland, the presence of algae was found to increase melting of high biomass regions by as much as 26% (Cook et al., ). As liquid water is considered a pre-requisite for glacier algal growth within ice (Williamson et al., ), the potential exists for establishment of a positive feedback loop between glacier algal growth, ice surface darkening, liquid water generation, and continued bloom proliferation. Accordingly, larger blooms are anticipated in years that experience the greatest overall melt (Tedstone et al., ; Williamson et al., ; Cook et al., ; Di Mauro et al., ).
Record melt years have recently been witnessed across the cryosphere, with the European Alps experiencing a record amount of glacier melt during the summer of 2022, when 5.9% of the total ice volume was lost (World Glacier Monitoring Service, ). In the Ötztal Alps of Austria, for example, the winter mass balance of Hintereisferner glacier was 43% below the decadal average of 2011–2020, paving the way for an extreme glacial mass loss during the summer season (Voordendag et al., ). Melting started earlier than normal, with the first positive degree days occurring in early May, initiating a nearly uninterrupted melt season that resulted in loss of 3,319 kg of ice per m−2, 3.4–5 times greater than the two proceeding years and 3.3 times greater than the 1991–2020 average (Voordendag et al., ). In the Swiss Alps, the average storage change for all Swiss glaciers over the 2022 melt season was estimated at 3.63 ± 0.26 km3 water equivalent, roughly 60% greater (three times higher) than the average annual loss over the past decade (Cremona et al., ). This was similarly driven by exceptionally low winter accumulation and high summer melting, the latter caused by a record number of 23 extreme heat wave events (defined as at least three consecutive days with a daily average maximum temperature >30°C) that occurred during summer 2022 (GLAMOS, ; Cremona et al., ). The consequences of these rapid melt events have a multitude of impacts on the surrounding environment and present many forms of threat to the safety of local communities (Beniston et al., ; Farinotti et al., ; Trautmann et al., ).
Within this context, 2022 was an important year to study glacier algal blooms, allowing to examine whether expected outcomes in bloom magnitude were apparent given the extreme melt conditions experienced. To this end, we report on a glacier algal bloom on Morteratsch glacier, Switzerland, during the 2022 summer melt season. Morteratsch glacier is the largest in the Bernina Range of Switzerland and is known to host extensive glacier algal blooms during most summer seasons (Di Mauro et al., , ). Previous studies have identified a long-term darkening trend of surface ice at Morteratsch glacier (Oerlemans et al., ; Wientjes et al., ), with subsequent studies highlighting the role of glacier algae in albedo decline at this location (Di Mauro et al., , ). The presence of both dominant glacier algal species, Ancylonema nordenskiöldii and A. alaskana, is known at Morteratsch glacier, with abundances reported at ~104 cells ml−1 during the 2016 melt season (Di Mauro et al., ). Here, we provide a characterization of the glacier algal bloom at Morteratsch glacier during August 2022 to examine whether record melt conditions resulted in record bloom proliferation within this location.
Methods
Field site and sampling locations
Characterization of the glacier algal bloom proliferating on surface ice of Morteratsch glacier, Switzerland, was undertaken over a 3-day period in late August 2022 (21st August 2022–23rd August 2022). Surface ice samples were collected across three main sites (one site sampled per day) distributed from the terminus of the glacier at 2,100 m up to 2,300 m elevation above sea level (Figure 1), hereafter referred to as lower (46.418N, 9.933E; Figure 2D), middle (46.415N, 9.934E; Figure 2C), and upper (46.414N, 9.934E; Figure 2B) sampling locations. At each location, n = 10 randomly selected areas of supraglacial surface ice were sampled at the start (10:00–12:00), middle (14:00–16:00), and end (17:00–19:00) of the day to capture potential diurnal dynamism in the biogeochemistry associated with blooms, with n = 30 samples per lower/middle/upper sampling location. All sampling was undertaken during clear, sunny-sky conditions (Figure 2). For each individual ice sample a 20 × 20 × 2 cm depth patch of surface ice was sampled using an ice saw and metal scoop, pre-sterilised with 10% ethanol, directly into a sterile Whirl-Pak bag using standard sampling techniques (Williamson et al., , ). Before collection of each ice sample, an image with a scale was taken, and the slope, aspect, elevation with respect to sea level, and GPS coordinates were recorded (directly adjacent to sampling patch to avoid contamination) using the Phyphox app on an Apple iPhone 14 (Staacks et al., ).
Figure 1
Figure 2
Sample processing
Ice samples were melted in the dark over 1–2 days at 4°C and homogenized prior to further processing. An initial 1.5 ml of each homogenized sample was placed into an individual 1.5 ml Eppendorf tube and fixed using Lugol's solution for subsequent cell counts and biovolume estimates at the University of Bristol, UK. A known volume of sample was then filtered through a pre-combusted (450°C for 5 h) 13-mm diameter GF/A filter (1.6 μm retention; Cytiva Whatman, Maidstone, UK), which was placed into an individual 1.5 ml Eppendorf tube and maintained at −20°C for subsequent cellular stoichiometry determination. Filtrates were collected into pre-acid-washed high density polyethylene (HDPE) bottles and maintained at −20°C for subsequent inorganic (ammonium, nitrate, nitrite, and phosphate) and organic (total organic carbon and nitrogen) aqueous geochemistry quantification. Samples were then driven to the University of Bristol, UK, and maintained within portable fridges at 4°C or freezers at −20°C.
Nutrient addition experiments
To examine potential bottom-up control of glacier algal blooms at Morteratsch glacier, a nutrient spiking experiment was conducted ex situ with glacier algal communities sampled from surface ice. Fifteen samples were collected as above from lower Morteratsch glacier (46.418N, 9.933E) on 26th August 2022 and transported in the dark at 1°C within a portable fridge to Photon Systems Instruments, Czech Republic, where they were melted over 2 days as above. Approximately 7 days after sample collection, an incubation experiment was established with melted surface ice samples and their constituent glacier algal communities. Melted samples were initially pooled and then 50 ml distributed into n = 5 replicates across each of five treatments: control (no nutrient addition), +ammonium (+, 10 μmol L−1), +nitrate (+, 10 μmol L−1), +phosphate (+, 1 μmol L−1), and +ALL (+ and + at 10 μmol L−1 and + at 2 μmol L−1). , , and spikes were added as -N, NO3-N, and -P standards, respectively (Thermo Fisher Scientific, UK). Treatments were chosen to provide ~10 times the reported ambient inorganic nitrogen and phosphorus concentrations in supraglacial ice (Wadham et al., ; Smith et al., ) and maintain a 10:1 nitrogen:phosphorus ratio within our “+ALL” treatment (McCutcheon et al., ). Incubations were conducted within 50 ml Corning culture flasks with vented caps (Corning, New York, USA) and maintained at 4°C under 500 μmol photons m−2 s−1 of light on a 16:8 Light:Dark cycle provided by warm white fluorescent tubes with no UV provision. Destructive sampling of n = 5 replicate incubations per treatment proceeded at 0, 24, 72, and 120 h for the determination of cellular abundance/biovolume, algal photophysiology, and cellular stoichiometry as described below.
Abundance and biovolume determination
The cellular abundance (cells ml−1) and average biovolume per species (μm3 cell−1) were assessed for all samples collected in the field and used for incubation experiments. Cellular abundance was measured by counting cells on a modified Fuchs Rosenthal Hemocytometer (0.2 mm by 1/16 mm2; Hawksley, Lancing, United Kingdom) using a bright field Olympus BX41 microscope (Germany). Images of each sample were taken at 10 × and 40 × magnification with a MicroPublisher 6 CCD camera attachment (Teledyne Photometrics, USA), from which cells were counted and the width and radius measured for each sample using ImageJ software (Schneider et al., ). Final cellular biovolumes were calculated assuming species to be an average cylinder (Hillebrand et al., ). Power analysis was applied to determine the number of cells to measure to achieve a margin of error < 250 μm3 per species/sample (Jones et al., ). The total algal biovolume per sample (μm3 ml−1) was calculated by multiplying species' cellular abundances per sample (cells ml−1) by their corresponding mean biovolumes (μm3 cell−1).
Photophysiology determination
Rapid light response curves (RLCs; Perkins et al., ) were performed using Pulse Amplitude Modulation (PAM) fluorometry to characterize the photophysiological states of glacier algal communities as sampled from Morteratsch glacier and throughout nutrient spiking experiments. Measurements were conducted on 3 ml subsamples using a Walz Water-PAM fluorometer with attached red-light emitter/detector cuvette system and stirrer (Walz GmBH, Germany). Each sample was dark adapted for a minimum of 5 min prior to RLC measurement. RLCs consisted of nine sequential light steps of 20 s duration ranging in irradiance from 0 to 3,000 μmol photons m−2 s−1 of Photosynthetically Active Radiation (PAR). The maximum quantum efficiency (Fv/Fm) was calculated from minimum (F0) and maximum (Fm) fluorescence yields measured in the dark-adapted state during the initial RLC step of 20 s darkness as Fv/Fm = Fm – Fo/Fm (Consalvey et al., ). Electron transport through photosystem II (PSII) was calculated across subsequent light steps in relative units (rETR) assuming an equal division of light between PSI and PSII as rETR = Y (PSII) × PAR × 0.5 (Consalvey et al., ). Analysis of all RLC data (rETR vs. PAR) followed Eilers and Peeters () for calculation of the relative maximum electron transport rate (rETRmax).
Aqueous geochemistry determination
Aqueous concentrations of , , , and were derived for all melted ice samples spectrophotometrically using a Gallery Plus Discrete Photometric Analyser (Thermo Fisher Scientific, UK). The limit of detection (LoD) for all nutrients was determined by the mean concentration plus three times the standard deviation of calibration blanks (n = 3). LoDs were 0.04 μmol L−1 (), 0.02 μmol L−1 (), 0.15 μmol L−1 −(-), and 0.15 μmol L−1 (). Precisions were ±2.0% (), ±1.0% (NO2), ±1.3% (), and ±2.0% () as determined by comparison with diluted 71.43 mmol L−1-N, -N, and -N and diluted 32.29 mmol L−1-P certified stock standards to a concentration of 3.6 μmol L−1 ( and ), 2.9 μmol L−1 (), and 6.1 μmol L−1 () (Sigma TraceCERT®).
Filtrate was also analyzed for total organic carbon (TOC) and total nitrogen (TN) concentrations via a TOC/TN Organic Carbon Analyser (Shimadzu, UK). Non-purgeable organic carbon (NPOC) was measured after the acidification of samples with 9N sulfuric acid and catalytic combustion at 720°C as CO2. TN was measured as NO by chemiluminescence. The LoD was 18.3 μmol L−1 (TOC) and 30 μmol L−1 (TN), precision was ±2.0% as determined by comparison with diluted 41.7 mmol L−1 TOC certified stock standards to a concentration of 8.3 mmol L−1 and ±1.7% by comparison with diluted 14.3 mmol L−1 TN certified stock standards to a concentration of 3.6 mmol L−1 (Sigma TraceCERT®).
Cellular stoichiometry
Glacier algal cellular carbon (C) and nitrogen (N) contents were determined for n = 5 samples from Morteratsch glacier for each location and time point. For this, GF/A filters were freeze-dried for 24 h to remove water, re-weighed and wrapped in individual 16 mm tin disks prior to elemental analysis using a Vario PYRO cube® (Elementar, Stockport, UK). The detection limits of elemental concentrations were 0.001% for both elements measured, and the coefficient of variation (CV) for C and N according to 12 replicates of an organic analytical standard (NC Soil Standard 338 40025, cert. 341506, C = 2.31%, N = 0.23%; ThermoFisher Scientific, Bremen, Germany) were 1.76% and 1.13%, respectively. The molar content of carbon and nitrogen per sample was derived from the total recorded carbon and nitrogen area as follows:
where Arx is the relative atomic mass of nutrient x (e.g., carbon), x [%] is the derived percentage of nutrient x present in a processed sample, w is the total weight of the processed sample. Data are presented as molar C:N ratios.
Data analysis
Analysis and plotting of data were completed using R v.4.2.1 (R Core Team, ). Each data series collected (cell abundance, biovolume, photophysiological, aqueous geochemistry, and stoichiometry) was tested for homogeneity of variance and normality of distribution. Statistical comparisons of photophysiology, aqueous geochemistry, and stoichiometry between sample locations and time points were investigated using one-way or two-way analysis of variance (ANOVA) with post-hoc Tukey Honest Significant Difference (HSD) analysis applied to all significant ANOVA results. ANOVA and Tukey HSD were also applied to comparisons between nutrient spiking treatments. To compare cell abundances and biovolumes between sampling locations and sampling times, Kruskal–Wallis rank sum tests (for non-parametric datasets) were employed. Linear regression analysis was employed to investigate potential relationships between cell abundance or biovolume and sample site slope. Aspect of the slope was treated as circular data (Von Mises distribution) and tested for correlation with cell abundance or biovolume using the circular correlation coefficient equivalent for Pearson's test for correlation (Fisher and Lee, ; Rao Jammalamadaka and Ramakrishna Sarma, ). Data were plotted using the ggplot2 package v.3.4.4.
Results and discussion
A substantial glacier algal bloom was apparent during August 2022 at Morteratsch glacier, Switzerland (Figure 1), producing a conspicuous discoloration of the glacier's surface ice, particularly toward the terminus (Figure 2A). Sampling throughout the diurnal cycle across three elevations on the glacier allowed a broad characterization of the bloom, highlighting compositional and spatiotemporal heterogeneities in glacier algal biomass across the ice surface. These demonstrated how bloom magnitude does not scale linearly with melt season magnitude but rather upper limits on bloom carrying capacity may be determined by meltwater-driven redistribution and eventual export of glacier algal cells from the glacier.
Cell abundance does not scale linearly with meltwater availability between seasons
Glacier algal communities were composed of both cosmopolitan glacier algal species at Morteratsch during 2022, with Ancylonema alaskana showing the greatest overall cellular abundance in all sampling locations as compared with A. nordenskiöldii (85% A. alaskana relative abundance vs. 15% A. nordenskiöldii across all data; Figure 3), though a more even split between species was apparent in our lowest sampling location based on total biovolume datasets (52% A. alaskana vs. 48% A. nordenskiöldii; Figure 4). These data are consistent with the study by Di Mauro et al. () who reported an overall greater cellular abundance of A. alaskana relative to A. nordenskiöldii for the 2016 bloom in the same site; though contrasting findings based on DNA recovery. At present, we hold little-to-no knowledge on the factors that determine glacier algal community composition within or between sites and seasons; a key facet of bloom ecology yet to be understood. Indeed, it is likely that we do not yet have a full understanding of the true diversity of glacier algae across the cryosphere (Remias et al., ).
Figure 3
Figure 4
Glacier algal cellular abundance ranged from 1.78 × 104 to 8.95 × 105 cells ml−1 during August 2022 on Morteratsch glacier; at the higher end of abundances reported across the cryosphere during large bloom years (Stibal et al., ; Williamson et al., ; Di Mauro et al., ) but not substantially greater than previous records. For example, Di Mauro et al. () reported an average of 2.03 × 104 cells ml−1 from the same site in September 2016, during which melt may have been as little as a third of that in 2022 (World Glacier Monitoring Service, , ). On the southwestern Greenland Ice Sheet (GrIS), known to harbor extensive glacier algal blooms, maximal abundances have been reported as 8.5 × 104 cells ml−1 (Stibal et al., ), 1.6 × 104 cells ml−1 (Williamson et al., ), and 3.0 × 105 cells ml−1 at the margin (Yallop et al., ); comparable to our observations in Morteratsch. This suggests a more complex relationship between melt and glacier algal abundance and/or the presence of other controls on overall bloom carrying capacity.
Spatiotemporal distribution of glacier algal cells
Spatiotemporal heterogeneity in glacier algal abundance recorded here indicated the potential melt-water-driven redistribution of (larger) glacier algal cells to down-glacier locations and presumptive export of biomass from the glacial system at Morteratsch; the latter a putative mechanism to set upper limits on bloom carrying capacity. At the diurnal scale, cellular abundance datasets recovered consistently higher abundances earlier in the day, when the volume of meltwater on the glacier surface was at its lowest (Figure 3). As melt increased throughout the day, corresponding estimates of abundance declined, reflecting dilution of assemblages by the higher melt fraction. The timing of sampling relative to the daily hydrograph is thus important for final estimates of abundance and should be considered in future sampling efforts.
Across our study locations, the abundance of A. nordenskiöldii was significantly higher in the lower region of the glacier (Figure 3), with a lower and more heterogenous abundance at mid and upper sites evidenced by a range of 0 cells ml−1 to 8.0 × 104 cells ml−1 across samples. While the more abundant A. alaskana (cells ml−1) was more evenly distributed across the glacier (Figure 3), when cellular biovolume (μm3 cell−1) was incorporated into estimates, a clear spike in total biovolume (μm3 ml−1) was apparent at our lower site, highlighting the presence of larger A. alaskana and A. nordenskiöldii cells in this region of the glacier (Figure 4). Interestingly, this heterogeneity in glacier algal distribution was independent of surface ice slope or aspect. Data thus reflected our previous findings (Williamson et al., ), whereby true spatial trends in GrIS glacier algal blooms at the km scale were reflected only by total biovolume datasets rather than estimates of cellular abundance alone.
The presence of more and larger cells at our lowest sampling location irrespective of slope or aspect could reflect preferential growth within this region of the glacier and/or preferential transport of larger cells down glacier by meltwater, serving to concentrate assemblages at the terminus, with subsequent export to downstream environments. At larger spatial scales, the duration since snow-line retreat is a known factor controlling the progression of blooms, such that areas of glacial ice uncovered the earliest (i.e., at lower altitudes) provide more time for glacier algal blooms to develop (Tedstone et al., ; Williamson et al., ). This dynamic has been noted across the ablation zone of the GrIS (Tedstone et al., ; Williamson et al., , ) as well as for valley glaciers (Takeuchi, ). However, sampling sites targeted here covered a relatively small altitudinal range (~100 meters of elevation gain) during a record melt year with rapid snow line retreat (GLAMOS, ; Cremona et al., ), and it is unlikely that substantial differences in exposure duration were experienced across our transect during 2022. PAM analysis did not indicate differences in community photophysiological state across our sampling locations (Figure 5), e.g., Fv/Fm remained consistently high across all locations at 0.66 ± 0.12 (Figure 5B). Additionally, aqueous geochemistry data and nutrient spiking incubations (see below) gave no indication of preferential conditions across our sampling locations.
Figure 5
Alternatively, larger cells (i.e., the filaments of A. nordenskiöldii and larger cells of A. alaskana) may have been preferentially transported over the weathering crust, concentrating in the lowest glacier sampling site. At the macro-scale, a clear concentration of glacier algal biomass was apparent at the glacier terminus throughout our sampling campaign (Figure 2). We also observed obvious down-glacier cell transport within meltwater flowing over the weathering crust via microscopic observations, indicating potential down-glacier concentration of cells (data not shown). Preferential movement of larger cells contrasts a previous report by Irvine-Fynn et al. (), which suggested that smaller cells are more likely to be transported through the weathering crust of the GrIS. However, it is likely that hydrological controls of cell distribution and weathering crust properties on steep valley glaciers differ from those of the flatter “dark zone” of the ice sheet; hence, more information is required on the meltwater-driven redistribution of glacier algal cells across the diversity of surface ice environments and its role in regulating bloom expansion vs. driving loss of cells from the system.
Geochemistry of the heterogenous glacier algal bloom
Characterization of the bulk-phase aqueous geochemistry (Figure 6) once again demonstrated the paradox of highly abundant glacier algal communities thriving within oligotrophic surface ice (Williamson et al., , ). Inorganic macro-nutrient concentrations ranged 0.22 – 4.36 μmol L−1 of , 0.04 – 1.38 μmol L−1 NO3−, 0.09 – 0.29 μmol L−1, and 0.22 – 0.85 μmol L−1 (Figure 6), with a shift to dominance of organic speciation apparent by our sampling dates, i.e., 142.70 – 811.36 μmol L−1 dissolved organic carbon (DOC) and 0.23 – 15.34 μmol L−1 dissolved organic nitrogen (DON), reflecting known phase-shifts in the dominant geochemistry associated with blooms (Holland et al., ). Previously, Williamson et al. () argued that glacier algae themselves may be adapted to low ambient nutrient concentrations through overall lower cellular N and P requirements. Indeed, we recorded similarly high C:N ratios here across sampling locations (average 12.28 ± 0.70 across all locations; Figure 7), with elevated C:N at our lowest sampling location corresponding to the presence of larger cells of both species; likely containing more carbon-rich phenolic pigmentation, rather than reflecting nutrient limitation at the terminus. To confirm this, algal communities sampled from our lower location were spiked with a series of nutrient additions (~10 times ambient , , and or a combination of all) and incubated for up to 120 h. Even after this 5-day duration, no clear indication of inorganic nutrient limitation was apparent for our sampled assemblages (Figure 8), with no difference across treatments in the maximum quantum yield in the dark-adapted state (Fv/Fm; Figure 8B) and minimal variance in rETRmax (Figure 8C).
Figure 6
Figure 7
Figure 8
These findings, taken together with previous studies (McCutcheon et al., ; Williamson et al., ), indicate that certainly at the point of sampling, glacier algae do not appear to be inorganic nutrient limited within surface ice environments of both mountain glaciers (here) or larger ice sheet systems (McCutcheon et al., ). This may reflect a true adaptation to low nutrient conditions (Williamson et al., ), efficient uptake and storage of those nutrients available within the system (Barcyté et al., ), and/or poor characterization of the available nutrients within the thin meltwater film that glacier algae reside. To date, researchers have typically employed the “bulk” method used here of ice sampling with subsequent melt and aqueous geochemistry characterization. This potentially dilutes the true nutrient available to glacier algae within their habitat, which may be enhanced due to ice solute exclusion and/or preferential elution of solute-rich fractions when glacier ice melts (Ginot et al., ; Holland et al., ). Thus, more work on true glacier algal nutrient requirements and nutrient availability within their habitat is needed to tackle the oligotrophic-bloom-paradox that is ubiquitous across the cryosphere.
Conclusion
Characterization of a glacier algal bloom during a record melt year demonstrated how bloom proliferation does not scale linearly with meltwater availability on a steep valley glacier. Instead, a threshold appears to have been crossed, whereby meltwater shifts from being a promotor of glacier algal growth to a regulator of overall carrying capacity within surface ice via lateral transport of cells in surface meltwater. Abundances apparent during the record 2022 melt season at Morteratsch were high but not notably above previous estimates across the cryosphere. In contrast, spatiotemporal patterning in glacier algal biomass taken together with both macro- and micro-observations of bloom dynamics highlighted the potential for significant melt-water-driven, down-glacier transport of (larger) cells and presumptive export from the glacier, providing a first-order control on bloom carrying capacity. The findings were corroborated by aqueous geochemistry, algal eco-physiology, and nutrient spiking datasets that failed to reveal any biogeochemical rationale for the glacier algal distribution observed, once again highlighting the oligotrophic-bloom-paradox often reported across the cryosphere. Bottom-up controls on bloom proliferation will likely remain minimal as more of the glacial cryosphere becomes unlocked by climate warming in the future. However, our findings indicate that exponential proliferation of blooms as meltwater increases through time is likewise unlikely in valley glacier settings, where unknown thresholds of meltwater availability are surpassed and physical transport and export of cells from the system dominate. Many key facets of bloom ecology remain unknown as we head into the peak melt century.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
JM: Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing—original draft, Writing—review & editing. EB: Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing—original draft, Writing—review & editing. ML: Investigation, Validation, Writing—review & editing. AB: Writing—review & editing. AT: Writing—review & editing. CW: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Writing—original draft, Writing—review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. Funding for this project was provided by the Leverhulme Trust (iDAPT RPG-2020-199 to CW). Funding was provided to CW from the Natural Environmental Research Council (NE/Y002636/1). Also AT was supported by the European Research Council award 818994 – Cassandra.
Acknowledgments
The authors would like to acknowledge doctoral training funding of EB and ML from the University of Bristol.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
BarcytéD.PilátováJ.MojzešP.NedbalováL. (2019). The arctic Cylindrocystis (Zygnematophyceae, Streptophyta) green algae are genetically and morphologically diverse and exhibit effective accumulation of polyphosphate. J. Phycol.56, 217–232. 10.1111/jpy.12931
2
BenistonM.StoffelM.HillM. (2011). Impacts of climatic change on water and natural hazards in the Alps: can current water governance cope with future challenges? Examples from the European “ACQWA” project. Environ. Sci. Policy14, 734–743. 10.1016/j.envsci.2010.12.009
3
ConsalveyM.PerkinsR. G.PatersonD. M.UnderwoodG. J. C. (2005). PAM fluorescence: a beginners guide for benthic diatomists. Diatom Res.20, 1–22. 10.1080/0269249X.2005.9705619
4
CookJ. M.TedstoneA. J.WilliamsonC.McCutcheonJ.HodsonA. J.DayalA.et al. (2020). Glacier algae accelerate melt rates on the south-western Greenland Ice Sheet. Cryosphere14, 309–330. 10.5194/tc-14-309-2020
5
CremonaA.HussM.LandmannJ. M.BornerJ.FarinottiD. (2023). European heat waves 2022: contribution to extreme glacier melt in Switzerland inferred from automated ablation readings. Cryosphere17, 1895–1912. 10.5194/tc-17-1895-2023
6
Di MauroB.BaccoloG.GarzonioR.GiardinoC.Massab,òD.PiazzalungaA.et al. (2017). Impact of impurities and cryoconite on the optical properties of the Morteratsch Glacier (Swiss Alps). Cryosphere11, 2393–2409. 10.5194/tc-11-2393-2017
7
Di MauroB.GarzonioR.BaccoloG.FranzettiA.PittinoF.LeoniB.et al. (2020). Glacier algae foster ice-albedo feedback in the European Alps. Sci. Rep.10, 4739. 10.1038/s41598-020-61762-0
8
EilersP. H. C.PeetersJ. C. H. (1988). A model for the relationship between light intensity and the rate of photosynthesis in phytoplankton. Ecol. Modell.42, 199–215. 10.1016/0304-3800(88)90057-9
9
FarinottiD.CremonaA.van TielM.HussM. (2023). Glaciers' Role as Water Resource in the Swiss Alps.Vienna: EGU General Assembly.
10
FisherN. I.LeeA. J. (1983). A correlation coefficient for circular data. Biometrika70, 327–332. 10.1093/biomet/70.2.327
11
GinotF.LenavetierT.DedovetsD.DevilleS. (2020). Solute strongly impacts freezing under confinement. Appl. Phys. Lett. 116. 10.1063/5.0008925
12
GLAMOS (2022). Available online at: https://doi.glamos.ch/data/massbalance/massbalance_2022_r2022.html (accessed November 2, 2023).
13
HalbachL.ChevrollierL.-A.DotingE. L.CookJ. M.JensenM. B.BenningL. G.et al. (2022). Pigment signatures of algal communities and their implications for glacier surface darkening. Sci. Rep.12, 17643. 10.1038/s41598-022-22271-4
14
HillebrandH.DürselenC.-D.KirschtelD.et al. (1999). Biovolume calculation for pelagic and benthic microalgae. J. Phycol.35403–424. 10.1046/j.1529-8817.1999.3520403.x
15
HollandA. T.WilliamsonC. J.SgouridisF.TedstoneA. J.McCutcheonJ.CookJ. M.et al. (2019). Dissolved organic nutrients dominate melting surface ice of the Dark Zone (Greenland Ice Sheet). Biogeosciences16, 3283–3296. 10.5194/bg-16-3283-2019
16
HollandA. T.WilliamsonC. J.TedstoneA. J.AnesioA. M.TranterM. (2022). Dissolved nitrogen speciation and concentration during spring thaw in the greenland ice sheet dark zone: evidence for microbial activity. Front. Earth Sci. Chin.10, 711560. 10.3389/feart.2022.711560
17
Irvine-FynnT. D. L.EdwardsA.StevensI. T.MitchellA. C.BuntingP.BoxJ. E.et al. (2021). Storage and export of microbial biomass across the western Greenland Ice Sheet. Nat. Commun.12, 3960. 10.1038/s41467-021-24040-9
18
JonesS. R.CarleyS.HarrisonM. (2003). An introduction to power and sample size estimation. Emerg. Med. J.20, 453–458. 10.1136/emj.20.5.453
19
McCutcheonJ.LutzS.WilliamsonC.CookJ. M.TedstoneA. J.VanderstraetenA.et al. (2021). Mineral phosphorus drives glacier algal blooms on the Greenland Ice Sheet. Nat. Commun.12, 570. 10.1038/s41467-020-20627-w
20
OerlemansJ.GiesenR. H.Van Den BroekeM. R. (2009). Retreating alpine glaciers: increased melt rates due to accumulation of dust (Vadret da Morteratsch, Switzerland). J. Glaciol.55, 729–736. 10.3189/002214309789470969
21
PerkinsR. G.MougetJ.-L.LefebvreS.LavaudJ. (2006). Light response curve methodology and possible implications in the application of chlorophyll fluorescence to benthic diatoms. Mar. Biol.149, 703–712. 10.1007/s00227-005-0222-z
22
ProcházkováL.RezankaT.NedbalováL.RemiasD. (2021). Unicellular versus Filamentous: The Glacial Alga Ancylonema alaskana comb. et stat. nov. and Its Ecophysiological Relatedness to Ancylonema nordenskioeldii (Zygnematophyceae, Streptophyta). Microorganisms 9, 1103. 10.3390/microorganisms9051103
23
R Core Team (2022). A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
24
Rao JammalamadakaS.Ramakrishna SarmaY. (1988). A Correlation Coefficient for Angular Variables. Available online at: dx.doi.org/ (accessed November 15, 2023).
25
RemiasD.HolzingerA.AignerS.LützC. (2012). Ecophysiology and ultrastructure of Ancylonema nordenskiöldii (Zygnematales, Streptophyta), causing brown ice on glaciers in Svalbard (high arctic). Polar Biol.35, 899–908. 10.1007/s00300-011-1135-6
26
RemiasD.ProcházkováL.NedbalováL.BenningL. G.LutzS. (2023). Novel insights in cryptic diversity of snow and glacier ice algae communities combining 18S rRNA gene and ITS2 amplicon sequencing. FEMS Microbiol. Ecol. 99, fiad134. 10.1093/femsec/fiad134
27
SchneiderC. A.RasbandW. S.EliceiriK. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat. Methods9, 671–675. 10.1038/nmeth.2089
28
SmithH.FosterR.McKnightD.LisleJ.LittmannS.KuypersM.et al. (2017). Microbial formation of labile organic carbon in Antarctic glacial environments. Nat. Geosci.10, 356–359. 10.1038/ngeo2925
29
StaacksS.HützS.HeinkeH.StampferC. (2018). Advanced tools for smartphone-based experiments: phyphox. Phys. Educ.53, 045009. 10.1088/1361-6552/aac05e
30
StibalM.BoxJ. E.CameronK. A.LangenP. L.YallopM. L.MottramR. H.et al. (2017). Algae drive enhanced darkening of bare ice on the greenland ice sheet. Geophys. Res. Lett.44, 463–471. 10.1002/2017GL075958
31
TakeuchiN. (2013). Seasonal and altitudinal variations in snow algal communities on an Alaskan glacier (Gulkana glacier in the Alaska range). Environ. Res. Lett. 8, 035002. 10.1088/1748-9326/8/3/035002
32
TedstoneA. J.BamberJ. L.CookJ. M.WilliamsonC. J.FettweisX.HodsonA. J.et al. (2017). Dark ice dynamics of the south-west Greenland Ice Sheet. Cryosphere11, 2491–2506. 10.5194/tc-11-2491-2017
33
TrautmannS.KnoflachB.StötterJ.ElsnerB.IllmerP.GeitnerC. (2023). Potential impacts of a changing cryosphere on soils of the European Alps: a review. Catena232, 107439. 10.1016/j.catena.2023.107439
34
VoordendagA.PrinzR.SchusterL.KaserG. (2023). Brief communication: the Glacier Loss Day as indicator for extreme glacier melt in 2022. Cryosphere2023, 1–7. 10.5194/tc-2023-49
35
WadhamJ. L.HawkingsJ.TellingJ.ChandlerD.AlcockJ.O'DonnellE.et al. (2016). Sources, cycling and export of nitrogen on the Greenland Ice Sheet. Biogeosciences13, 6339–6352. 10.5194/bg-13-6339-2016
36
WientjesI. G. M.WalR. S. W. V.de ReichartG. J.SluijsA.OerlemansJ. (2011). Dust from the dark region in the western ablation zone of the Greenland ice sheet. Cryosphere5, 589–601. 10.5194/tc-5-589-2011
37
WilliamsonC. J.AnesioA. M.CookJ.TedstoneA.PonieckaE.HollandA.et al. (2018). Ice algal bloom development on the surface of the Greenland Ice Sheet. FEMS Microbiol. Ecol. 94, fiy025. 10.1093/femsec/fiy025
38
WilliamsonC. J.CameronK. A.CookJ. M.ZarskyJ. D.StibalM.EdwardsA. (2019). Glacier algae: a dark past and a darker future. Front. Microbiol.10, 524. 10.3389/fmicb.2019.00524
39
WilliamsonC. J.CookJ.TedstoneA.YallopM.McCutcheonJ.PonieckaE.et al. (2020). Algal photophysiology drives darkening and melt of the Greenland Ice Sheet. Proc. Natl. Acad. Sci. U. S. A.117, 5694–5705. 10.1073/pnas.1918412117
40
WilliamsonC. J.Turpin-JelfsT.NicholesM. J.YallopM. L.AnesioA. M.TranterM. (2021). Macro-nutrient stoichiometry of glacier algae from the southwestern margin of the Greenland Ice Sheet. Front. Plant Sci.12:673614. 10.3389/fpls.2021.673614
41
World Glacier Monitoring Service (2021). WGMS 2021: Fluctuations of Glaciers Database. Zurich: World Glacier Monitoring Service.
42
World Glacier Monitoring Service (2022). WGMS 2022: Fluctuations of Glaciers Database. Zurich: World Glacier Monitoring Service.
43
YallopM. L.AnesioA. M.PerkinsR. G.CookJ.TellingJ.FaganD.et al. (2012). Photophysiology and albedo-changing potential of the ice algal community on the surface of the Greenland ice sheet. ISME J.6, 2302–2313. 10.1038/ismej.2012.107
Summary
Keywords
glacier algae, streptophyte, Morteratsch, glacier melt, Alps, climate change
Citation
Millar JL, Broadwell ELM, Lewis M, Bowles AMC, Tedstone AJ and Williamson CJ (2024) Alpine glacier algal bloom during a record melt year. Front. Microbiol. 15:1356376. doi: 10.3389/fmicb.2024.1356376
Received
15 December 2023
Accepted
08 January 2024
Published
20 February 2024
Volume
15 - 2024
Edited by
Mark Alexander Lever, The University of Texas at Austin, United States
Reviewed by
Marek Stibal, Charles University, Czechia
Jarishma K. Gokul, University of Pretoria, South Africa
Updates
Copyright
© 2024 Millar, Broadwell, Lewis, Bowles, Tedstone and Williamson.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Christopher J. Williamson c.williamson@bristol.ac.uk
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.