Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 28 June 2021
Sec. Aquatic Photosynthetic Organisms
This article is part of the Research Topic Ice and Snow Algae View all 10 articles

Macro-Nutrient Stoichiometry of Glacier Algae From the Southwestern Margin of the Greenland Ice Sheet

  • 1Bristol Glaciology Centre, School of Geographical Sciences, University of Bristol, Bristol, United Kingdom
  • 2School of Biological Sciences, University of Bristol, Bristol, United Kingdom
  • 3Department of Environmental Science, Aarhus University, Aarhus, Denmark

Glacier algae residing within the surface ice of glaciers and ice sheets play globally significant roles in biogeochemical cycling, albedo feedbacks, and melt of the world’s cryosphere. Here, we present an assessment of the macro-nutrient stoichiometry of glacier algal assemblages from the southwestern Greenland Ice Sheet (GrIS) margin, where widespread glacier algal blooms proliferate during summer melt seasons. Samples taken during the mid-2019 ablation season revealed overall lower cellular carbon (C), nitrogen (N), and phosphorus (P) content than predicted by standard microalgal cellular content:biovolume relationships, and elevated C:N and C:P ratios in all cases, with an overall estimated C:N:P of 1,997:73:1. We interpret lower cellular macro-nutrient content and elevated C:N and C:P ratios to reflect adaptation of glacier algal assemblages to their characteristic oligotrophic surface ice environment. Such lower macro-nutrient requirements would aid the proliferation of blooms across the nutrient poor cryosphere in a warming world. Up-scaling of our observations indicated the potential for glacier algal assemblages to accumulate ∼ 29 kg C km2 and ∼ 1.2 kg N km2 within our marginal surface ice location by the mid-ablation period (early August), confirming previous modeling estimates. While the long-term fate of glacier algal autochthonous production within surface ice remains unconstrained, data presented here provide insight into the possible quality of dissolved organic matter that may be released by assemblages into the surface ice environment.

Introduction

Microbial communities that reside on the surfaces of glaciers and ice sheets play globally significant roles in carbon (C) and nutrient cycling and surface ice melt (Hodson et al., 2008; Stibal et al., 2012; Anesio et al., 2017; Williamson et al., 2018, 2020). Of particular importance are Streptophyte “glacier algae” (Williamson et al., 2019), whose presence in surface ice lowers the bare ice albedo, enhances solar energy absorption, and drives surface melt through the process of “biologically driven albedo reduction” (Yallop et al., 2012; Stibal et al., 2017; Williamson et al., 2018, 2020; Cook et al., 2020). On the surface of the Greenland Ice Sheet (GrIS), summer blooms of glacier algae are responsible for widespread albedo decline that has paralleled accelerating surface melt since the early 1990s (Yallop et al., 2012; Tedesco et al., 2016; Stibal et al., 2017; van den Broeke et al., 2017; Cook et al., 2020; Tedstone et al., 2020; Williamson et al., 2020), with recent estimates attributing an additional 5.5–8.0 Gt of runoff to glacier algal growth along the western ice sheet margin; 6–9% of the total runoff (Cook et al., 2020). Given that melt of the GrIS is the single largest cryospheric contributor to global eustatic sea level rise (Bamber et al., 2018), constraining bloom dynamics into the future remains a significant research priority (Williamson et al., 2019, 2020; Cook et al., 2020; Tedstone et al., 2020).

Glacier algal blooms initiate following snow line retreat (Williamson et al., 2018), with population doubling times on the GrIS ranging 3.75–5.5 days, and maximal cell densities approaching 105 cells mL–1 of melt water during major bloom years (Yallop et al., 2012; Stibal et al., 2017; Williamson et al., 2018, 2020). At these cell densities, widespread albedo decline is driven by significant secondary phenolic pigmentation produced by glacier algae to protect their low-light adapted chloroplasts (Williamson et al., 2020). By dissipating the intercepted incident irradiance as heat, this secondary pigmentation also generates liquid water proximal to the cells, driving surface ice ablation while providing access to nutrient and other resources required to promote algal growth (Dial et al., 2018; Williamson et al., 2020). Accordingly, growth proceeds as a function of bare-ice melt duration, such that strong patterning in accumulated biomass is apparent across the ablation zone, with maximal cell densities accumulated within the most marginal regions that experience the longest ablation periods, and a decreasing trend toward the equilibrium line (Williamson et al., 2020). As the ablation period progresses into polar winter, the fate of accumulated biomass remains unknown, although glacier algal species have been shown to overwinter in alpine locations (Remias et al., 2012b), and active glacier algal assemblages have been observed in GrIS shallow surface ice prior to snow line retreat (Nicholes et al., 2019).

One of the biggest questions remaining about glacier algal blooms asks, What factors limit the distribution and magnitude of blooms in space and time? While physical conditions such as snow-pack height, light availability, and temperature produce a first order control on the ability of blooms to form and proliferate in any given year (Tedstone et al., 2017; Williamson et al., 2020), little information currently exists on potential “bottom-up” or “top-down” regulation of blooms in supraglacial environments, restricting abilities to project bloom occurrence into the future (Williamson et al., 2019; McCutcheon et al., 2021). Within the marine environment, for example, top-down pressures such as zooplankton grazing and/or bottom-up availability of nutrients represent fundamental controls on phytoplankton biomass (Huppert et al., 2002). Supraglacial environments are characterized by truncated, microbially dominated trophic structures (Anesio and Laybourn-Parry, 2012; Anesio et al., 2017) and highly oligotrophic conditions (Hawkings et al., 2016; Wadham et al., 2016), indicating that large potential exists for bottom-up limitation of glacier algal blooms.

To date, field observations of ambient hydrochemistry during bloom progression have demonstrated a bulk phase shift toward organic over inorganic nutrient resources within surface ice (Holland et al., 2019), which coupled with 28-times lower secondary production relative to primary production (Yallop et al., 2012; Nicholes et al., 2019) has been interpreted to imply inefficient remineralization of inorganic nutrient resources within GrIS surface ice, and a potential mechanism of bottom-up control (Nicholes et al., 2019). Recently, McCutcheon et al. (2021) provided the first evidence for such bottom-up control of glacier algal blooms on the GrIS, highlighting inorganic phosphorus (Pi) limitation and the importance of locally sourced hydroxylapatite in supporting bloom proliferation, and likely important roles for heterotrophic bacterial and fungal communities in accelerating apatite weathering and thus P availability to glacier algal communities. This study did not, however, directly quantify glacier algal abundance within samples, precluding direct calculation of cellular elemental quotas. The importance of this information is illustrated by the findings of Holland et al. (2019) for the same bloom and sampling sites, who concluded that there was sufficient inorganic macro-nutrient availability within surface ice to support the magnitude of glacier algal bloom apparent during the 2016 ablation season when assuming Redfield stoichiometry for glacier algal C, nitrogen (N), and P requirements.

Fundamental to deciphering the importance of bottom-up controls on glacier algal blooms is therefore knowledge of the macronutrient cellular requirements of glacier algal cells. While Redfield stoichiometry provides an important framework against which ambient nutrient concentrations can be contrasted (Redfield, 1958), deviations from this ratio are well documented across several algal lineages (Hecky et al., 1993; Geider and La Roche, 2002; Quigg et al., 2003), with individual cell stoichiometry shown to be dynamic relative to a plethora of drivers (Geider and La Roche, 2002; Quigg et al., 2003; Dickman et al., 2006; Finkel et al., 2016). Here, we provide an estimate of glacier algal macro-nutrient (C, N, and P) cellular stoichiometry in assemblages sampled from the surface of the GrIS. A snapshot of glacier algae stoichiometry was determined for assemblages residing within surface ice sampled from the south western ice sheet margin during the 2019 ablation season in order to provide a first order approximation of the elemental requirements of glacier algal assemblages and to investigate the potential for divergence from Redfield dynamics. This data is important for efforts to project the occurrence of blooms into the future and to constrain bloom impacts to cycles of C and macro-nutrients.

Materials and Methods

Study Area and Sampling Details

Surface ice containing glacier algal assemblages was sampled across August 8–10 (total number of samples = 28) from the marginal ablation zone in the south western GrIS proximal to Point 660. For all samples, the top 2 cm of surface ice was sampled using a pre-cleaned ice saw, transferred into sterile Whirl-Pak® bags (Madison, WI, United States), and melted slowly in the dark over 24 h at 4°C. The melted surface ice was subsequently homogenized and sub-sampled for further analyses including 1 mL sub-sampled into 15 mL centrifuge tubes and fixed immediately with 2% glutaraldehyde final concentration for subsequent algal cell enumeration; triplicate 100–200 mL subsamples filtered onto pre-combusted 25 mm diameter glass microfiber filters (0.7 μm retention; Thermo Fisher Scientific, Pittsburgh, PA, United States), frozen immediately at −20°C for subsequent determination of glacier algal C, N content; and duplicate 100–200 mL subsamples filtered onto pre-combusted 47 mm diameter glass microfiber filters (0.7 μm retention; Thermo Fisher Scientific, Pittsburgh, PA, United States), frozen immediately at −20°C for subsequent determination of glacier algal P content. All samples were transported back to the University of Bristol for subsequent processing.

Algal Cell Enumeration and Cell Biovolume Estimation

Algal abundance (cells mL–1) was determined using methods described by Williamson et al. (2018) on 1 mL aliquots of melted surface ice using a modified Fuchs-Rosenthal hemocytometer (0.2 mm by 1/16 mm2; Hawksley, Lancing, United Kingdom) on a Leica M205 C stereomicroscope (Wetzlar, DE) with attached GXCAM HiChrome-Met HD microscope camera (GT Vision Ltd., Stansfield, United Kingdom). Glacier algal cellular biovolume (μm3) was determined from measurements of cell length and diameter taken using ImageJ software (version 1.52n), for which glacier algal cells were considered to be regular cylinders (Hillebrand et al., 1999).

Algal Carbon and Nitrogen

Using methods outlined by Lorrain et al. (2003), the C (after the removal of carbonates by acid fumigation) and N content of glacier algae retained on 25 mm diameter glass microfiber filters (0.7 μm retention; Thermo Fisher Scientific, Pittsburgh, PA, United States), which had been freeze-dried over 24 h, were determined using an elemental analyzer (Elementar vario PYRO cube®, Langenselbold, Hesse, DE). 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 six replicates of an organic analytical standard (NC Soil Standard 338 40025, cert. 133317, C = 2.29%, N = 0.21%; Elemental Microanalysis Ltd., United Kingdom) were 5.60% and 2.74%, respectively. To determine the quantity of C and N associated with algal biomass (μg mL–1 of surface ice), the bulk glacier algal abundance (cells mL–1) was multiplied by the quantity of cellular C and N (μg cell–1).

Algal Phosphorus

Concentrations of total P (Pt) and Pi associated with glacier algae were determined using methods adapted from Hedley and Stewart (1982) and Stibal et al. (2008). Paired 47 mm diameter glass microfiber filters (0.7 μm retention) were added to separate pre-cleaned 15 mL polypropylene centrifuge tubes (Thermo Fisher Scientific), frozen at −80°C and subsequently freeze-dried for 24 h. After freeze-drying, 1 mL of ethanol-free chloroform (CHCl3; Thermo Fisher Scientific) was added to one from each pair of centrifuge tubes. These tubes were then agitated three times for 10 s using a vortex mixer at 15 min intervals before being placed under a fume hood for 24 h to allow the CHCl3 to evaporate. Both the unfumigated and CHCl3-fumigated filters from each pair were then amended with 6 mL (∼ 1:60 w/v retained sample:extractant) of 1 M magnesium chloride (MgCl2; Thermo Fisher Scientific), agitated at 200 rpm on a reciprocal shaker for 16 h and centrifuged at 1,217 × g at 0°C for 12 min. An aliquot of 1.5 mL from the resulting extracts was transferred to a 9 mL muffled digest tube, mixed with 0.5 mL of an oxidizing agent comprised of 9 g of NaOH and 40 g of potassium persulfate (K2S2O8; Thermo Fisher Scientific) in 1 L of Milli-Q® water, and autoclaved for 60 min at 121°C (digested). The undigested and digested MgCl2 extracts pertaining to both the unfumigated and CHCl3-fumigated filters were then filtered to 0.45 μm using mixed cellulose ester membranes (Whatman®, Maidstone, United Kingdom). Concentrations of Pi and Pt were determined in the undigested and digested MgCl2 extracts, respectively, using a standard orthophosphate colorimetric technique on a photometric meter (GalleryTM Plus Discrete Analyzer, Thermo Fisher Scientific, Waltham, MA, United States). The CV (six replicate mid-range standards) for Pi and Pt were ≤0.20%. In addition, the limits of detection (LoD; three times the standard deviation of six replicate method blanks) for Pi and Pt were 4.73 μg L–1 and 1.93 μg L–1, respectively. Samples were blank corrected when blank concentrations exceeded the detection limits. All reagents used were of analytical grade. Algal Pi and Pt were calculated, respectively, as the differences between P concentrations in the undigested and digested MgCl2 extracts for the CHCl3-fumigated and unfumigated filters within each pair.

Ambient Chemistry

Concentrations of dissolved organic C (DOC) and dissolved total N (DTN) in filtered ice samples were determined on a total organic C analyzer (Shimadzu TOC-LCPH, Kyoto, Japan) coupled with a total N measuring unit (Shimadzu TNM-L). The ammonium-N (NH4+-N), nitrate-N (NO3-N), and phosphate-P (PO43–-P) contents of the same samples were quantified on a photometric meter (GalleryTM Plus Discrete Analyzer). According to six replicate standards, the CV for the DOC, DTN, NH4+-N, NO3-N, and PO43–-P analyses were 3.71%, 7.38%, 0.20%, 0.26%, and 0.35%, respectively. The corresponding detection limits for the same analyses were 51.40 μg L–1, 39.12 μg L–1, 1.19 μg L–1, 3.13 μg L–1, and 4.95 μg L–1. Dissolved organic N (DON) was calculated for each sample as the difference between concentrations of DTN and dissolved inorganic N (DIN; sum of NH4+-N and NO3-N).

Data Analyses

All plots and statistical analyses were performed using R version 3.4.1 in accordance with methods outlined by Crawley (2005).

Results

Ambient Chemistry of Surface Ice

Major C and macro-nutrient phases were determined for GrIS surface ice from which glacier algal assemblages were sampled to provide contextual hydrochemistry information to contrast with glacier algal cellular stoichiometries. Dissolved nutrient phases measured from melted surface ice were dominated by organic phases (Table 1), with the majority comprised of DOC, which exhibited concentrations ∼ 18 times greater than DON. In contrast, the concentrations of dissolved inorganic nutrients fell at or below the LoD. Dissolved inorganic N was only measurable in the form of NH4+ and accounted for ∼ 4% of the DTN, while PO43–-P concentrations were below the LoD in all cases.

TABLE 1
www.frontiersin.org

Table 1. Concentrations of dissolved organic carbon (DOC), total dissolved nitrogen (DTN), organic nitrogen (DON), ammonium-nitrogen (NH4+-N), nitrate-nitrogen (NO3-N), and phosphate-phosphorus (PO43–-P) in surface ice from the south western margin of the Greenland Ice Sheet (GrIS) during the 2019 ablation season (mean ± standard error, n = 12).

Assemblage Characteristics

Glacier algal assemblages ranged in abundance from 0 to 8.97 × 105 cells mL–1 and were comprised predominantly of Ancylonema cf. nordenskiöldii (67.9%) and cf. Mesotaenium berggrenii (31.7%; Figure 1A). Mean cellular biovolumes for A. cf. nordenskiöldii and cf. M. berggrenii were 2,788 ± 121 μm3 cell–1 and 2,197 ± 131 μm3 cell–1, respectively (Figure 1B). In contrast, fewer than 1% of cells were identified as Cylindrocystis brebissonii, though mean cellular biovolumes of these cells were an order of magnitude larger than A. cf. nordenskiöldii and cf. M. berggrenii (Figure 1B).

FIGURE 1
www.frontiersin.org

Figure 1. The abundance (A) and cellular biovolume (B) of Ancylonema cf. nordenskiöldii (Ancy., n = 19 abundance, n = 215 biovolume), cf. Mesotaenium berggrenii (Meso., n = 20 abundance and n = 140 biovolume), Cylindrocystis brebissonii (Cylin., n = 2 abundance, n = 5 biovolume), and the total glacier algal assemblage (All, n = 21 abundance n = 360 biovolume) sampled in surface ice from the south western margin of the Greenland Ice Sheet (GrIS) during the mid-2019 ablation season. Plots show median ± interquartile range.

Cellular C and N contents of glacier algal assemblages averaged 106 ± 35 pg C cell1 and 4.5 ± 1.3 pg N cell1, respectively (Figures 2A,B), with cellular atomic C:N ratios ranging from 15.6 to 40.8. Using measured values of cellular C and N, we estimate an average of 27.4 ± 6.0 μg C mL1 and 1.2 ± 0.25 μg N mL1 to be contained within glacier algal biomass within the surface ice environment during our sampling period. In contrast to cellular C and N, quantifying algal P presented challenges as concentrations fell below the limit of detection (LoD), with Pt and Pi only detectable in seven and four of 28 samples, respectively. Across these, mean Pt is 0.14 ± 0.09 pg cell1 (Figure 2C). If organic P is defined as the difference between mean concentrations of Pt and Pi, we determined that up to 86% of algal P may be contained within the organic phase (∼ 0.12 pg cell1). Using the mean Pt, we estimate the atomic C:N:P stoichiometry of glacier algal assemblages to be 1,997:73:1; however, given the difficulty in detecting cellular P, this ratio should be considered with caution.

FIGURE 2
www.frontiersin.org

Figure 2. Cellular concentrations of (A) carbon (C, n = 28), (B) nitrogen (N, n = 28), (C) total phosphorus (Pt, n = 4), and (C) inorganic phosphorus (Pi, n = 4) of glacier algal assemblages comprised of Ancylonema cf. nordenskiöldii, cf. Mesotaenium berggrenii, and Cylindrocystis brebissonii in surface ice from the south western margin of the Greenland Ice Sheet (GrIS) during the mid-2019 ablation season. Plots show median ± interquartile range.

Discussion

Understanding the macronutrient requirements of glacier algal cells is an important first step on the way to constraining potential bottom-up controls, which may ultimately restrict bloom magnitude and extent (McCutcheon et al., 2021). The present study examined the C, N, and P content of glacier algal assemblages sampled during the mid-2019 ablation season (early August) from the southwestern margin of the GrIS, providing a direct measurement of glacier algal cellular macro-nutrient content and stoichiometry. Through this we identify deviation from standard cell-constituent:volume ratios and classic Redfield stoichiometry, and we are able to estimate glacier-algal-associated macro-nutrient reservoirs within this highly oligotrophic environment.

The absolute cellular macronutrient content of glacier algal assemblages fell within the range of values reported across a diversity of algal lineages (Montagnes et al., 1994), though they deviated from cell-constituent:volume relationships established across microalgal taxa. For example, glacier algal cellular C and N content averaged 106 ± 35 pg C cell–1 and 4.5 ± 1.3 pg N cell–1 for cell volumes that ranged ∼ 2,000–3,000 μm3. While absolute cellular contents were within the range reported across a host of microalgae (Montagnes et al., 1994; Finkel et al., 2016), values were conspicuously lower than would be predicted using established cell-constituent:volume ratios for cells of their size, whereby 2,000–3,000 μm3 cells would be expected to contain 203–304 pg C cell–1 and 40–62 pg N cell–1, respectively (Montagnes et al., 1994). Glacier algal cellular C was thus approximately half that expected based on cell volume, with cellular N an order of magnitude lower.

Laboratory and field studies have consistently identified species-level differences in microalgal elemental requirements that reflect evolutionary histories and acclimation to environmental conditions (Geider and La Roche, 2002; Finkel et al., 2016; Garcia et al., 2018). Such differences are a product of the adaptations of cellular architecture and biochemistry that correspond to changes in macromolecular composition, with protein the primary reservoir of cellular N, phospholipids, polyphosphates, and nucleic acids the reservoirs of cellular P, and cellular C largely determined by the combination of protein, lipid, and carbohydrate (Geider and La Roche, 2002; Finkel et al., 2016; Barcyte et al., 2020). In this respect, lower cellular requirements for C and N may reflect adaptation of Streptophyte glacier algae to their oligotrophic surface ice environment; documented here and in numerous previous works to be highly deficient in inorganic nutrient supplies (Table 1; Hawkings et al., 2016; Wadham et al., 2016; Holland et al., 2019). It may also reflect the lower growth rates of glacier algal communities from icy environments (∼ 5 days doubling time; Stibal et al., 2017; Williamson et al., 2018) compared to more temperate microalgal taxa (e.g., Garcia et al., 2018), whereby species with lower growth rates have associated lower protein, and thus N, content (Finkel et al., 2016). Finally, divergence from typical cell volume:constituent relationships may also reflect the highly vacuolized nature of glacier algal cells, which results in an overall higher water content compared to many other algae (Remias et al., 2009, 2012a). Separating the phylogenetic signature of microalgal elemental requirements from e.g., dynamic responses to in situ nutrient regimes (see below) is ideally achieved through the assessment of actively growing microalgal cells cultured under nutrient replete conditions (Geider and La Roche, 2002; Finkel et al., 2016). To date, however, an inability to culture glacier algal taxa ex situ has hindered advances in our understanding of their physiology and interactions with key environmental stressors (Williamson et al., 2019), necessitating field-based approaches as here.

In contrast to C and N contents, measurement of cellular P from glacier algal assemblages proved challenging, with few samples yielding detectable concentrations. This was consistent with our aqueous hydrochemistry data whereby all dissolved PO43–-P concentrations from melted surface ice samples fell below the LoD (LoD = 4.95 μg L–1). For glacier algae, just four of 28 samples yielded quantifiable P, with a mean Pt content of 0.14 ± 0.09 pg P cell–1 across samples. To compare to published values for major marine phyla (Ho et al., 2003), a glacier algal cell volume of 2,500 μm3 was applied to convert to units of mmol P per liter cell volume (mmol L–1 cell volume), yielding ∼ 1.50 mmol P L–1 cell volume for our glacier algal assemblages. This is substantially lower than values reported for 15 marine eukaryote species (∼ 9.1–250 mmol P L–1 cell volume; see Table 2 in Ho et al., 2003), though likely reflects the nutrient deplete conditions under which glacier algal assemblages grow. For example, the cyanobacterium Microcystis aeruginosa demonstrated much lower cellular P content ranging from 0.2 to 1.2 pg P cell–1 when cultured under ambient P conditions that ranged 0–256 μg P L–1, with lower intracellular content associated with lower ambient P availability (Ghaffar et al., 2017). Similarly, for more closely related planktonic desmid species (Zygnematophyceae), Spijkerman and Coesel (1996) demonstrated how Cosmarium abbreviatum var. planctonicum originating from an oligo-meso-trophic lake was consistently capable of greater biomass production for a given amount of P (across the range 0.5–10 μmol P l–1) than Staurastrum pingue, which originated from a eutrophic lake. Cellular P content ranged 0.2–2.17 pg P cell–1, comparable to values reported here for glacier algae, with C. planctonicum shown to be an affinity species, possessing a competitive advantage in an environment where growth is permanently P limited (Spijkerman and Coesel, 1996). While limited in sample size, our data demonstrate how glacier algal macro-nutrient content broadly reflects ambient inorganic nutrient availability within surface ice habitats of the GrIS, and suggests an overall lower cellular macro-nutrient requirement likely advantageous in the ultra-oligotrophic supraglacial. These findings are consistent with recent advances in understanding potential bottom-up controls on glacier algal bloom proliferation (McCutcheon et al., 2021).

Assessment of glacier algal cellular macro-nutrient quotients allowed for a first order approximation of organic C and N stocks within assemblages inhabiting our sampling location on the marginal southwestern GrIS. Based on cellular C and N contents and abundances recorded across our samples, we calculated an average of 27.4 ± 6.0 μg C mL–1 and 1.2 ± 0.25 μg N mL–1 to be contained within glacier algal cells within our marginal sampling location. Upscaling from these point observations to the km2 scale (i.e., μg C or N mL–1 × 1.061 = kg C or N km2, after Williamson et al., 2018) provides an estimate of ∼ 29.0 kg C km2 and ∼ 1.2 kg N km2 stored within glacier algal assemblages inhabiting surface ice at the mid-ablation season (early August). Previously, Williamson et al. (2018) modeled glacier algal net production in southwestern Greenland forced by the number of days since snowline retreat. They estimated an overall average of 15.82 ± 8.14 kg C km2 produced by glacier algal assemblages across their 8.24 × 104 km2 model region during the 2016 ablation period, with spatial variability ranging from <10 kg C km2 toward the equilibrium line up to ∼ 40–50 kg C km2 at their most marginal model regions by the end of the 2016 ablation season (Williamson et al., 2018). Our measurements of glacier algal cellular macro-nutrient content corroborate modeling estimates by Williamson et al. (2018) and extend these trends into the most marginal zone of the ice sheet, omitted by this previous study given uncertainties regarding microbial activity within such locations (Hodson et al., 2010; Stibal et al., 2012). We confirm here that glacier algal cells represent a significant producer of autochthonous organic C and N within the most marginal zone of the GrIS, likely driving the high rates of net production previously recorded for surface ice in this region (e.g., Musilova et al., 2017).

Additional to the cellular quotient of macro-nutrients, our data further demonstrate strong divergence in the relative abundance, i.e., cellular stoichiometry, of glacier algal C, N, and P relative to Redfield stoichiometry (106:16:1 C:N:P; Redfield, 1958), which we interpret here to represent the true stoichiometry of glacier algal communities under normal surface ice conditions as opposed to a dynamic response to e.g., nutrient limitation. Atomic C:N ratios of glacier algal assemblages ranged 15.6–40.8 across all samples, approximately two to six times greater than the Redfield C:N ratio of 6.6, and while P measurements proved challenging (see above), we estimate here an overall glacier algal C:N:P stoichiometry in the range 1,997:73:1. Under high light, low nutrient conditions, algal photosynthesis is often uncoupled from growth, with surplus photosynthate reallocated to storage carbohydrates and lipids that require minimal or no N and/or P inputs (Geider and La Roche, 2002; James et al., 2011; Pichrtová et al., 2014; Talmy et al., 2014). This results in a C-rich but N- and P-poor assemblage (Sterner et al., 1997), higher C:N and C:P ratios than predicted by Redfield stoichiometry, and typically conclusion that the community under study is nutrient limited (Berman-Frank, 1999; Dickman et al., 2006). For glacier algae residing within the high-light surface ice environment, it is likely that their abundant secondary phenolic pigmentation contributes to the elevated C:N and C:P ratios recorded here. This C-rich pigmentation is required to protect the cells against excess irradiance and provide a mechanism to convert abundant light energy into heat in order to liberate meltwater adjacent to the cell (Remias et al., 2012b; Dial et al., 2018; Williamson et al., 2018, 2020), likely also acting as an effective sink for excess photosynthate. However, given the characteristic oligotrophy of surface ice environments (Holland et al., 2019), and the contrasting ability of glacier algae to form widespread blooms under such conditions cryosphere wide (Williamson et al., 2019), we argue that our low absolute macro-nutrient content and elevated C:N and C:P ratios also likely signify overall reduced requirements for N and P resources by glacier algal communities, a potentially key adaptation to life in oligotrophic surface ice. This assertion is consistent with the findings of McCutcheon et al. (2021), who though concluding mineral phosphorus to be a first order control on glacier algal bloom presence, demonstrated that responses to nutrient addition (specifically phosphorous) were only evident in incubated field assemblages after 5 days of incubation, suggesting sufficient nutrient resources for glacier algal growth under ambient conditions. For context, this study reported highly comparable C:N:P stoichiometry to the present study, with POC:N:Porg ranging 690:48:1 to 2,615:196:1 across eight samples dominated by glacier algae at their main study location (McCutcheon et al., 2021). This putative lower macro-nutrient requirement would allow glacier algal assemblages to continue to progress into new oligotrophic areas of bare ice as the climate warms and the GrIS ablation zone expands into the future (Ryan et al., 2019), exacerbating local albedo decline and surface ablation in a positive feedback mechanism (Cook et al., 2020). It may also form a key consideration for ongoing efforts to culture glacier algal taxa ex situ under laboratory conditions.

To aid the balancing of internal stoichiometry, microalgae can actively release excess C into the ambient environment via the production and release of exopolymeric substances (EPS) (Berman-Frank, 1999; Hessen and Anderson, 2008; Palmucci et al., 2011) and/or divert excess photosynthate into the production of abundant secondary pigmentation as discussed above (Berman-Frank, 1999). For glacier algae blooming within ablating surface ice environments, EPS production has been postulated as a potentially important mechanism to aid retention of cells within the melting surface ice matrix (Williamson et al., 2019). While glacier algal EPS release remains unquantified, higher cellular C:N and C:P ratios likely impact the quality of dissolved organic matter available within surface ice environments, with consequences for other functional groups. For example, the dissolved organic nutrient phases that dominated the ambient hydrochemistry data presented here showed a DOC:DON ratio of ∼18, and a DOC:DOP ratio in the range 1,000–5,000 (Table 1). These data are highly consistent with Holland et al. (2019) who documented molar dissolved organic nutrient ratios in surface ice dominated by a glacier algal bloom further into the southwestern GrIS, with DOC:DON ranging 16–17 and DOC:DOP ranging 800–2,000. For that same 2016 glacier algal bloom, Nicholes et al. (2019) further demonstrated how bacterial secondary production was on average 28 times lower than rates of glacier algal primary production, indicating an inefficient microbial loop potentially limited by N and P resources (Holland et al., 2019; Nicholes et al., 2019). It is thus likely that glacier algal cellular stoichiometry is reflected in the quality of dissolved organic matter available within surface ice environments during blooms and may potentially limit the activity of associated communities. Constraining the short-term (within season) and long-term (inter-annual) fate of glacier algal autochthonous organic carbon and macro-nutrient production within surface ice will be the next step toward understanding their role as the dominant primary producer within supraglacial systems, and their wider impacts to biogeochemical cycling.

Conclusion

We provide here a snapshot of glacier algal cellular macro-nutrient content and stoichiometry from assemblages sampled during the mid-ablation season at the ice margin of the southwestern GrIS. Our findings highlight comparatively lower cellular C, N, and P absolute content in glacier algal cells that deviates from standard cellular content:biovolume relationships, and high C:N and C:P ratios that likely reflect adaptation to this highly oligotrophic surface ice environment as well as the dominance of abundant secondary phenolic pigmentation within glacier algal cells. Based on these observations, we confirm the role of glacier algal blooms in the most marginal region of the ice sheet in driving autochthonous macro-nutrient accumulation, estimating that approximately 29 kg C km2 and 1.2 kg N km2 was amassed by glacier algal assemblages by the middle of the 2019 ablation season. While the long-term fate of these nutrient reservoirs remains unconstrained, our findings indicate the potential for comparatively low quality (i.e., low N and P content) dissolved organic matter release into the surface ice by glacier algal assemblages, with consequences for associated heterotrophic activity.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

CW led the design of the study and completed all fieldwork. TT-J carried out all laboratory work. TT-J and MN performed the data analysis. All authors contributed to the manuscript.

Funding

This work was funded as part of the UK Natural Environment Research Council Consortium Grant “Black and Bloom” (NE/M021025/1). CW acknowledges support from The Leverhulme Trust Research Project Grant “iDAPT” (RPG-2020-199).

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.

Acknowledgments

We would like to recognize Dr. Fotis Sgouridis for providing guidance relating to elemental concentrations.

References

Anesio, A. M., and Laybourn-Parry, J. (2012). Glaciers and ice sheets as a biome. Trends Ecol. Evol. 27, 219–225.

Google Scholar

Anesio, A. M., Lutz, S., Chrismas, N. A. M., and Benning, L. G. (2017). The microbiome of glaciers and ice sheets. Npj Biofilms Microbiomes 3:10.

Google Scholar

Bamber, J. L., Westaway, R. M., Marzeion, B., and Wouters, B. (2018). The land ice contribution to sea level during the satellite Era. Environ. Res. Lett. 13:063008. doi: 10.1088/1748-9326/aac2f0

CrossRef Full Text | Google Scholar

Barcyte, D., Pilátová, J., Mojzeš, P., and Nedbalová, L. (2020). The arctic Cylindrocystis (Zygnematophyceae, Streptophyta) green algae are genetically and morphologically diverse and exhibit effective accumulation of polyphosphate. J. Phycol. 56, 217–232. doi: 10.1111/jpy.12931

PubMed Abstract | CrossRef Full Text | Google Scholar

Berman-Frank, I. (1999). Balanced growth in aquatic plants: myth or reality? phytoplankton use the imbalance between carbon assimilation and biomass production to their strategic advantage. BioScience 49, 29–37. doi: 10.2307/1313491

CrossRef Full Text | Google Scholar

Cook, J. M., Tedstone, A. J., Williamson, C. J., McCutcheon, J., Hodson, A., Dayal, A., et al. (2020). Glacier algae accelerate melt rates on the south-western greenland ice sheet. Cryosphere 14, 309–330. doi: 10.5194/tc-14-309-2020

CrossRef Full Text | Google Scholar

Crawley, M. J. (2005). Statistics: an Introduction Using R. Chichester, UK: John Wiley & Sons, Ltd.

Google Scholar

Dial, R. J., Ganey, G. Q., and Skiles, M. (2018). What color should glacier algae be? an ecological role for red carbon in the cryosphere. FEMS Microbiol. Ecol. 94:fiy007.

Google Scholar

Dickman, E. M., Vanni, M. J., and Horgan, M. J. (2006). Interactive effects of light and nutrients on phytoplankton stoichiometry. Oecologia 149, 676–689. doi: 10.1007/s00442-006-0473-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Finkel, Z. V., Follows, M. J., Liefer, J. D., Brown, C. M., Benner, I., and Irwin, A. J. (2016). Phylogenetic diversity in the macromolecular composition of microalgae. PLoS One 11:e0155977. doi: 10.1371/journal.pone.0155977

PubMed Abstract | CrossRef Full Text | Google Scholar

Garcia, N. S., Sexton, J., Riggins, T., Brown, J., Lomas, M. W., and Martiny, A. C. (2018). High Variability In Cellular Stoichiometry Of Carbon, Nitrogen, And Phosphorous Within Classes Of Marine Eukaryotic Phytoplankton Under Sufficient Nutrient Conditions. Front. Microbiol. 9:543.

Google Scholar

Geider, R. J., and La Roche, J. (2002). Redfield revisited: variability of C:N:P in marine microalgae and its biochemical basis. Eur. J. Phycol. 37, 1–17. doi: 10.1017/s0967026201003456

CrossRef Full Text | Google Scholar

Ghaffar, S., Stevenson, R. J., and Khan, Z. (2017). Effect of phosphorous stress on Microcystis aeruginosa growth and phosphorous uptake. PLoS One 12:e0174349. doi: 10.1371/journal.pone.0174349

PubMed Abstract | CrossRef Full Text | Google Scholar

Hawkings, J., Wadham, J., Tranter, M., Telling, J., and Bagshaw, E. (2016). The greenland ice sheet as a hot spot of phosphorus weathering and export in the arctic. Global Biogeochem. Cycles 30, 191–210. doi: 10.1002/2015gb005237

CrossRef Full Text | Google Scholar

Hecky, R. E., Campbell, P., and Hendzel, L. L. (1993). The stoichiometry of carbon, nitrogen, and phosphorus in particulate matter of lakes and oceans. Limnol. Oceanogr. 38, 709–724. doi: 10.4319/lo.1993.38.4.0709

CrossRef Full Text | Google Scholar

Hedley, M. J., and Stewart, J. W. B. (1982). Method to measure microbial phosphate in soils. Soil Biol. Biochem. 14, 377–385. doi: 10.1016/0038-0717(82)90009-8

CrossRef Full Text | Google Scholar

Hessen, D. O., and Anderson, T. R. (2008). Excess carbon in aquatic organisms and ecosystems: physiological, ecological, and evolutionary implications. Limnol. Oceanogr. 53, 1685–1696. doi: 10.4319/lo.2008.53.4.1685

CrossRef Full Text | Google Scholar

Hillebrand, H., Dürselen, C. D., Kirschtel, D., Pollingher, U., and Zohary, T. (1999). Biovolume calculation for pelagic and benthic microalgae. J. Phycol. 35, 403–424. doi: 10.1046/j.1529-8817.1999.3520403.x

CrossRef Full Text | Google Scholar

Ho, T. Y., Quigg, A., Finkel, Z. V., Milligan, A. J., Wyman, K., Falkowski, P. G., et al. (2003). The elemental composition of some marine phytoplankton. J. Phycol. 39, 1145–1159.

Google Scholar

Hodson, A., Anesio, A. M., Tranter, M., Fountain, A., Osborn, M., Priscu, J., et al. (2008). Glacial ecosystems. Ecol. Monogr. 78, 41–67.

Google Scholar

Hodson, A., Boggild, C., Hanna, E., Hyubrechts, P., Langfort, H., Cameron, K., et al. (2010). The cryoconite ecosystem on the greenland ice sheet. Ann. Glaciol. 51, 123–129. doi: 10.3189/172756411795931985

CrossRef Full Text | Google Scholar

Holland, A. T., Williamson, C. J., Sgouridis, F., Tedstone, A. J., McCutcheon, J., Cook, J. M., et al. (2019). Dissolved organic nutrients dominate melting surface ice of the dark zone (Greenland ice sheet). Biogeosciences 16, 3283–3296. doi: 10.5194/bg-16-3283-2019

CrossRef Full Text | Google Scholar

Huppert, A., Blasius, B., and Stone, L. (2002). A model of phytoplankton blooms. Am. Nat. 159, 156–171.

Google Scholar

James, G. O., Hocart, C. H., Hillier, W., Chen, H., Kordbacheh, F., Price, G. D., et al. (2011). Fatty acid profiling of Chlamydomonas reinhardtii under nitrogen deprivation. Bioresour. Technol. 102, 3343–3351. doi: 10.1016/j.biortech.2010.11.051

PubMed Abstract | CrossRef Full Text | Google Scholar

Lorrain, A., Savoye, N., Chauvaud, L., Paulet, Y. M., and Naulet, N. (2003). Decarbonation and preservation method for the analysis of organic C and N contents and stable isotope ratios of low-carbonated suspended particulate material. Anal. Chim. Acta 491, 125–133. doi: 10.1016/s0003-2670(03)00815-8

CrossRef Full Text | Google Scholar

McCutcheon, J., Lutz, S., Williamson, C., Cook, J. M., Tedstone, A. J., Vanderstraeten, A., et al. (2021). Mineral phosphorous drives glacier algal blooms on the greenland ice sheet. Nat. Commun. 12:570.

Google Scholar

Montagnes, D. J. S., Berges, J. A., Harrison, P. J., and Taylor, F. J. R. (1994). Estimating Carbon, Nitrogen, Protein, and Chlorophyll a from volume in marine phytoplankton. Limnol. Oceanogr. 39, 1044–1060. doi: 10.4319/lo.1994.39.5.1044

CrossRef Full Text | Google Scholar

Musilova, M., Tranter, M., Wadham, J., Telling, J., Tedstone, A., and Anesio, A. M. (2017). Microbially driven export of labile organic carbon from the greenland ice sheet. Nat. Geosci. 10, 360–365. doi: 10.1038/ngeo2920

CrossRef Full Text | Google Scholar

Nicholes, M. J., Williamson, C. J., Tranter, M., Holland, A., Poniecka, E., Yallop, M. L., et al. (2019). Bacterial dynamics in supraglacial habitats of the greenland ice sheet. Front. Microbiol. 10:1366.

Google Scholar

Palmucci, M., Ratti, S., and Giordano, M. (2011). Ecological and evolutionary implications of carbon allocation in marine phytoplankton as a function of nitrogen availability: a fourier transform infrared spectroscopy approach. J. Phycol. 47, 313–323. doi: 10.1111/j.1529-8817.2011.00963.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pichrtová, M., Kulichová, J., and Holzinger, A. (2014). Nitrogen limitation and slow drying induce desiccation tolerance in conjugating green algae (Zygnematophyceae, Streptophyta) from polar habitats. PLoS One 9:e113137. doi: 10.1371/journal.pone.0113137

PubMed Abstract | CrossRef Full Text | Google Scholar

Quigg, A., Finkel, Z. V., Irwin, A. J., Rosenthal, Y., Ho, T.-Y., Reinfelder, J. R., et al. (2003). The evolutionary inheritance of elemental stoichiometry in marine phytoplankton. Nature 425, 291–294. doi: 10.1038/nature01953

PubMed Abstract | CrossRef Full Text | Google Scholar

Redfield, A. C. (1958). The biological control of chemical factors in the environment. Am. Sci. 46, 205–221.

Google Scholar

Remias, D., Holzinger, A., Aigner, S., and Lutz, C. (2012a). Ecophysiology and ultrastructure of Ancylonema nordenskiöldii (Zygnematales, Streptophyta), causing brown ice on glaciers in Svalbard (high arctic). Polar Biol. 35, 899–908. doi: 10.1007/s00300-011-1135-6

CrossRef Full Text | Google Scholar

Remias, D., Holzinger, A., and Lütz, C. (2009). Physiology, ultrastructure and habitat of the ice alga Mesotaenium berggrenii (Zygnematophyceae, Chlorophyta) from glaciers in the European Alps. Phycologia 48, 302–312. doi: 10.2216/08-13.1

CrossRef Full Text | Google Scholar

Remias, D., Schwaiger, S., Aigner, S., Leya, T., Stuppner, H., and Lütz, C. (2012b). Characterization of an UV- and VIS-Absorbing, purpurogallin-derived secondary pigment new to algae and highly abundant in Mesotaenium berggrenii (Zygnematophyceae, Chlorophyta), an extremophyte living on glaciers. FEMS Microbiol. Ecol. 79, 638–648. doi: 10.1111/j.1574-6941.2011.01245.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryan, J. C., Smith, L. C., van As, D., Cooley, S. W., Cooper, M. G., Pitcher, L. H., et al. (2019). Greenland ice sheet surface melt amplified by snowline migration and bare ice exposure. Sci. Adv. 5:eaav3738. doi: 10.1126/sciadv.aav3738

PubMed Abstract | CrossRef Full Text | Google Scholar

Spijkerman, E., and Coesel, P. F. M. (1996). Phosphorous uptake and growth kinetics of two planktonic desmid species. Eur. J. Phycol. 31, 53–60. doi: 10.1080/09670269600651191

CrossRef Full Text | Google Scholar

Sterner, R. W., Elser, J. J., Fee, E. J., Guildford, S. J., and Chrzanowski, T. H. (1997). The light:nutrient ratio in lakes: the balance of energy and materials affects ecosystem structure and process. Am. Nat. 150, 663–684. doi: 10.1086/286088

PubMed Abstract | CrossRef Full Text | Google Scholar

Stibal, M., Box, J. E., Cameron, K. A., Langen, P. L., Yallop, M. L., Mottram, H., et al. (2017). Algae drive enhanced darkening of bare ice on the greenland ice sheet. Geophys. Res. Lett. 44, 11463–11471.

Google Scholar

Stibal, M., Telling, J., Cook, J., Mak, K. M., Hodson, A., Anesio, A. M., et al. (2012). Environmental controls on microbial abundance and activity on the greenland ice sheet: a multivariate analysis approach. Microb. Ecol. 63, 74–84. doi: 10.1007/s00248-011-9935-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Stibal, M., Tranter, M., Telling, J., and Benning, L. G. (2008). Speciation, phase association and potential bioavailability of phosphorus on a svalbard glacier. Biogeochemistry 90, 1–13. doi: 10.1007/s10533-008-9226-3

CrossRef Full Text | Google Scholar

Talmy, D., Blackford, J., Hardman-Mountford, N. J., Polimene, L., Follows, M. J., and Geider, R. J. (2014). Flexible C: N ratio enhances metabolism of large phytoplankton when resource supply is intermittent. Biogeosciences 11, 4881–4895. doi: 10.5194/bg-11-4881-2014

CrossRef Full Text | Google Scholar

Tedesco, M., Doherty, S., Fettweis, X., Alexander, P., Jeyaratnam, J., and Stroeve, J. (2016). The darkening of the greenland ice sheet: trends, drivers, and projections (1981-2100). Cryosphere 10, 477–496. doi: 10.5194/tc-10-477-2016

CrossRef Full Text | Google Scholar

Tedstone, A. J., Bamber, J. L., Cook, J. M., Williamson, C. J., Fettweis, X., Hodson, A. J., et al. (2017). Dark ice dynamics of the south-west Greenland Ice Sheet. Cryosphere 11, 2491–2506. doi: 10.5194/tc-11-2491-2017

CrossRef Full Text | Google Scholar

Tedstone, A. J., Cook, J. M., Williamson, C. J., Hofer, S., McCutcheon, J., Irvine-Fynn, T., et al. (2020). Algal growth and weathering crust state drive variability in Western Greenland Ice sheet ice albedo. Cryosphere 14, 521–538. doi: 10.5194/tc-14-521-2020

CrossRef Full Text | Google Scholar

van den Broeke, M., Box, J., Fettweis, X., Hanna, E., Noel, B., Tedesco, M., et al. (2017). Greenland ice sheet surface mass loss: recent developments in observation and modelling. Curr. Climate Chang. Rep. 3, 345–356. doi: 10.1007/s40641-017-0084-8

CrossRef Full Text | Google Scholar

Wadham, J. L., Hawkings, J., Telling, J., Chandler, D., Alcock, J., and Lawson, E. (2016). Sources, cycling and export of nitrogen on the Greenland ice sheet. Biogeoscience 13, 6339–6352. doi: 10.5194/bg-13-6339-2016

CrossRef Full Text | Google Scholar

Williamson, C. J., Anesio, A. M., Cook, J., Tedstone, A., Poniecka, E., Holland, A., et al. (2018). Ice algal bloom development on the surface of the Greenland ice sheet. FEMS Microbiol. Ecol. 94:fiy025.

Google Scholar

Williamson, C. J., Cameron, K. A., Cook, J. M., Zarsky, J. D., Stibal, M., and Edwards, A. (2019). Glacier algae: a dark past and a darker future. Front. Microbiol. 10:524.

Google Scholar

Williamson, C. J., Cook, J. M., Tedstone, A., Yallop, M. L., McCutcheon, J., Poniecka, E., et al. (2020). Algal photophysiology drives darkening and melt of the Greenland ice sheet. Proc. Natl. Acad. Sci. USA 117, 5694–5705. doi: 10.1073/pnas.1918412117

PubMed Abstract | CrossRef Full Text | Google Scholar

Yallop, M. L., Anesio, A. M., Perkins, R. G., Cook, J. M., Telling, J., Fagan, D., 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. doi: 10.1038/ismej.2012.107

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: glacier algae, C:N:P, Greenland Ice Sheet, stoichiometry, supraglacial

Citation: Williamson CJ, Turpin-Jelfs T, Nicholes MJ, Yallop ML, Anesio AM and Tranter M (2021) Macro-Nutrient Stoichiometry of Glacier Algae From the Southwestern Margin of the Greenland Ice Sheet. Front. Plant Sci. 12:673614. doi: 10.3389/fpls.2021.673614

Received: 27 February 2021; Accepted: 23 April 2021;
Published: 28 June 2021.

Edited by:

Linda Nedbalová, Charles University, Czechia

Reviewed by:

Elly Spijkerman, University of Potsdam, Germany
Daniel Remias, University of Applied Sciences Upper Austria, Austria

Copyright © 2021 Williamson, Turpin-Jelfs, Nicholes, Yallop, Anesio and Tranter. 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.