An Ecological Perspective on Dolomite Formation in Great Salt Lake, Utah

Dunham, Eric C.; Fones, Elizabeth M.; Fang, Yihang; Lindsay, Melody R.; Steuer, Christopher; Fox, Nicholas; Willis, Madelyne; Walsh, Alatna; Colman, Daniel R.; Baxter, Bonnie K.; Lageson, David; Mogk, David; Rupke, Andrew; Xu, Huifang; Boyd, Eric S.

doi:10.3389/feart.2020.00024

ORIGINAL RESEARCH article

Front. Earth Sci., 11 February 2020

Sec. Microbiological Chemistry and Geomicrobiology

Volume 8 - 2020 | https://doi.org/10.3389/feart.2020.00024

An Ecological Perspective on Dolomite Formation in Great Salt Lake, Utah

$\r\nEric C. Dunham$ Eric C. Dunham¹

Melody R. Lindsay¹

Nicholas Fox³

Alatna Walsh³

David Lageson³

David Mogk³

Andrew Rupke⁶

Huifang Xu²

Eric S. Boyd^1*

¹Department of Microbiology and Immunology, Montana State University, Bozeman, MT, United States
²Department of Geoscience, University of Wisconsin-Madison, Madison, WI, United States
³Department of Earth Sciences, Montana State University, Bozeman, MT, United States
⁴Department of Land Resources and Environmental Sciences, Montana State University, Bozeman, MT, United States
⁵Great Salt Lake Institute, Department of Biology, Westminster College, Salt Lake City, UT, United States
⁶Utah Geological Survey, Salt Lake City, UT, United States

Protodolomite was detected in sediments from Great Salt Lake (GSL), Utah, United States, that have no history of elevated temperature or pressure, conditions that are thought to promote dolomitization of sedimentary carbonates. Protodolomite was abundant in a non-oolitic sediment core from the South Arm (SA) of GSL at < 17% salinity but was rare in an oolitic sediment core collected from the North Arm (NA) of GSL at > 26% salinity. Protodolomite was also abundant in a non-oolitic NA sediment hand sample yet was absent in a nearby oolitic sediment hand sample in locations that likely receive allochthonous nutrients. Protodolomite was not detected in benthic photosynthetic mats from the SA. However, the mats comprised aragonite with halite and minimal calcite; benthic photosynthetic mats do not form in the NA. To begin to identify potential controls on the formation of protodolomite in the SA and NA of GSL, the composition and abundance of 16S rRNA gene transcripts, carbon cycling activities, and porewater geochemistry of the sediment cores were characterized. Transcripts affiliated with a dominant halophilic, heterotrophic sulfate-reducing bacterium were detected in the uppermost sections of the SA core and their abundance was positively correlated with rates of acetate oxidation/assimilation and concentrations of sulfide. Differences in the quality of organic matter between the SA and NA cores, as indicated by carbon to nitrogen ratios, indicate fresh deposition of photosynthetic biomass at the SA sediment core site but not in the NA sediment core site. Sediment grains from the SA core exhibit micrometer-sized euhedral protodolomite crystals that were not detected in the NA core. Collectively, these observations suggest that deposition of photosynthetic biomass drives the development of a sharp, anoxic lens of heterotrophic sulfate reduction. Sulfide, in turn, may promote dehydration of Mg²⁺-water complexes and primary protodolomite nucleation and growth. The scarcity of dolomite in the NA sediment core may result from constraints imposed by a combination of extreme hypersalinity and depositional environment on phototrophs and sulfate reducers, their activities, and the thermodynamics of protodolomite formation.

Introduction

Dolomite [CaMg(CO₃)₂] is common in ancient sedimentary rocks of marine origin and is generally thought to have formed through diagenetic processes involving high temperature and pressure or through post-depositional replacement of calcite or aragonite (CaCO₃) (Rosenberg and Holland, 1964; Rosenberg et al., 1967; Hardie, 1987). The crystal structure of dolomite differs from that of calcite in that the ordering of calcium (Ca²⁺) and magnesium (Mg²⁺) ions into alternating layers results in reduced crystal symmetry (Gregg et al., 2015). Prior studies have distinguished high- and very-high-magnesium calcite [HMC, ∼4–25% MgCO₃ content; VHMC, 25–36% MgCO₃ content; both without cation ordering (Gregg et al., 2015)] from protodolomite (> 36% MgCO₃, weak cation ordering) and dolomite (> 36% MgCO₃ and strong cation ordering) (Fang and Xu, 2019). It has been demonstrated that VHMC formation always precedes the formation of ordered dolomite in high-temperature laboratory synthesis experiments (Gregg et al., 2015). While common in older rocks, dolomite and protodolomite are rare in modern marine sediments despite waters being supersaturated with respect to these minerals, which should favor their formation (Fairbridge, 1957). Albeit rare, dolomite and/or protodolomite have been detected in some modern day sedimentary environments that have no history of high temperature or pressure, including marginal hypersaline environments (e.g., Alderman, 1958; Vasconcelos et al., 1995; Van Lith et al., 2002; Pace et al., 2016) and evaporitic environments (e.g., Wells, 1962; Bontognali et al., 2010), among others.

Various thermodynamic, kinetic, chemical, biological, and hydrological factors have been proposed to explain the formation of dolomite in modern day environments (Machel and Mountjoy, 1986). For example, it has been suggested that dolomite could form abiotically in saline environments with water chemistries above thermodynamic saturation with respect to dolomite (e.g., many saline to hypersaline environments), that are alkaline due to microbial activity, and that feature elevated temperatures (Machel and Mountjoy, 1986). Others suggest that biology is a key factor in dolomite formation (Vasconcelos et al., 1995) via either metabolic activities (Vasconcelos and Mckenzie, 1997; Van Lith et al., 2002; Moore et al., 2004; Kenward et al., 2009; Lu et al., 2018) or components of cells [e.g., extracellular polymeric substances (EPS) (Bontognali et al., 2010; Zhang et al., 2012b; Shen et al., 2015; Lu et al., 2018)] that create localized geochemical conditions or nucleation points, respectively, that are potentially favorable for dolomite formation. In particular, the activity of sulfate (SO₄^2–) reducing bacteria (SRB) can reduce the availability of SO₄^2–, which decreases ion pairing with Mg²⁺, potentially making it more available for replacement of Ca²⁺ in calcite or aragonite (Baker and Kastner, 1981; Vasconcelos and Mckenzie, 1997; Warthmann et al., 2000; Bontognali et al., 2010). Alternatively, it has been suggested that removal of SO₄^2– has little effect on the potential for dolomite formation, but rather it is the accumulation of sulfide from SRBs that promotes formation of this mineral (Zhang et al., 2012a, 2013; Lu et al., 2018). In this role, sulfide is thought to adsorb to calcite facies, thereby reducing the energy required for dehydration of Mg²⁺-water complexes on the same surfaces. Still others have suggested that salinity [e.g., sodium chloride (NaCl)] can exert strong influence on dolomite formation. For example, salinity may negatively influence the potential for dolomite formation since increased NaCl can promote complexation of Mg²⁺ and Ca²⁺ with SO₄^2– and Cl^– (Machel and Mountjoy, 1986), thereby rendering Mg²⁺ less available for dolomite formation. However, high salinities, in combination with high Ca²⁺ to Mg²⁺ ratios and SO₄^2– reduction (SR) activity, have been associated with the formation of dolomite in several hypersaline systems (e.g., Van Lith et al., 2002).

SO₄^2– reducing bacteria are common to many of the contemporary hypersaline environments where dolomite has been detected, including Great Salt Lake (GSL), Utah, United States (Brandt et al., 2001; Boyd et al., 2017). GSL is a hypersaline, terminal lake in the western United States (Hassibe and Keck, 1991) that exhibits a range of salinities due to localized freshwater input (∼95% of the total) into its southern end (Gwynn, 1996; Belovsky et al., 2011). Between 1900 and 1959, the salinity in GSL ranged, on average, from 20 to 27% (Stephens, 1990). In 1959, a rock and gravel railroad causeway was constructed, dividing the lake into two “arms” termed the North Arm (NA) and the South Arm (SA) (Cannon and Cannon, 2002) (Figure 1). Restricted water flow between these two arms due to the construction of the causeway has resulted in salinities that reach saturation (> 30%) in the NA and are < 20% in the SA, with localized areas of reduced salinity attributable to freshwater input (Baxter et al., 2005).

FIGURE 1

Figure 1. Map of Great Salt Lake (GSL) with sampling locations indicated by white circles with crosses. The railroad causeway separating the North Arm (NA) from the South Arm (SA), typical NA and SA salinities, the primary rivers draining into GSL, and major faults (hatched areas and dashed lines) are indicated. The distribution of microbialites is depicted based on mapping conducted by Eardley (1938). Structure modified from Colman et al. (2002) and UGS Quaternary Fault and Fold Database (Survey, 2013). Abbreviations: NAN, North Arm north hand sample collection site; SI, Stansbury Island microbial mat sampling site.

SO₄^2– reducing bacteria are abundant components of benthic sediment communities in GSL (Boyd et al., 2017) and can reach densities of 10⁷–10⁸ cells cm^–3 sediment (Brandt et al., 2001). The ecological success of SRB in GSL is likely due, in part, to elevated concentrations of SO₄^2– in the lake that can exceed 30 g L^–1 (Brandt et al., 2001; Gwynn, 2007; Boyd et al., 2017). Additionally, productive phototrophic plankton communities in the overlying water column (Post, 1977; Wurtsbaugh et al., 2011; Lindsay et al., 2017) can deliver fresh biomass to the benthos to support heterotrophic respiration. Finally, thick (∼1 cm) phototrophic, benthic microbial mats are commonplace along the shallow margins of the lake (Bouton et al., 2016; Pace et al., 2016; Lindsay et al., 2017, 2019) and may serve as an additional source of organic carbon capable of supporting heterotrophic SRB (Lindsay et al., 2019).

Previous studies have shown that the rate of SR is three orders of magnitude lower in sediments from GSL at elevated salinity (27% salinity) than in those at lower salinity (12%) (Brandt et al., 2001). Likewise, estimates of the abundance of cells performing SR and the abundance of 16S rRNA transcripts affiliated with SRB in sediments collected from hypersaline GSL locations were several orders of magnitude lower than those in less saline conditions (Brandt et al., 2001; Boyd et al., 2017). Activities of SRB sampled from hypersaline sediments in GSL increased when the salinity was artificially diluted, indicating that these populations are not well adapted to high salt conditions (Brandt et al., 2001). Importantly, benthic phototrophic microbial mats are estimated to contribute up to 30% of the productivity in the SA of GSL (Wurtsbaugh et al., 2011). Benthic phototrophic mats are not known to exist in the NA of GSL (Lindsay et al., 2017), likely due to constraints imposed by salinity on the primary architect of these mats, Cyanobacteria within the genus Euhalothece (Garcia-Pichel et al., 1998; Lindsay et al., 2017, 2019). The absence of benthic phototrophic mats, and therefore low levels of inorganic carbon fixation, in the NA may place constraints on the delivery of organic carbon capable of supporting heterotrophic SR in benthic sediments.

Previous studies have documented precipitation of aragonite in microbial mats (e.g., Reid et al., 2000; Schneider et al., 2013), including those that inhabit the SA of GSL (Pace et al., 2016). In GSL, aragonite precipitation was suggested to be localized to portions of the mat where SRB are actively degrading components of cyanobacterial biomass (Pace et al., 2016). However, other studies based on 16S rRNA gene (Lindsay et al., 2017, 2019) and transcript sequencing (M.R. Lindsay, unpublished data), the latter of which was conducted over a diel cycle, indicate that SRBs comprise only a minor portion (< 2–3%) of the communities associated with mats sampled from Bridger Bay in GSL, very close to where the previous study was conducted. Pace et al. (2016) also suggest that extracellular organic matter within the photosynthetic community may serve as a point of nucleation for aragonite mineral precipitation and growth (Pace et al., 2016). Moreover, it was suggested that partial replacement of aragonite with dolomite may take place within the mat (Pace et al., 2016), leading to dolomite preservation in benthic organo-carbonate structures dated to ∼2000 to 10,000 BP (Bouton et al., 2016) that cover ∼300 and 700 km² of the NA and SA lake bed (Eardley, 1938; Baskin, 2014).

In the present study, we sought to further characterize the role of microbial processes in the formation of dolomite in GSL. Specifically, we hypothesized that an abundant supply of SO₄^2– in sediment porewaters and infusion of nutrient-rich photosynthate and photosynthetic biomass from microbial mats promotes heterotrophic SR activity, which in turn promotes precipitation of VHMC and, potentially, (proto)dolomite. Moreover, it was hypothesized that (proto)dolomite would be more prevalent in the SA when compared to the NA due to constraints imposed by salinity on the productivity of the NA ecosystem (Lindsay et al., 2017, 2019), the activity and abundances of SRB (Brandt et al., 2001; Boyd et al., 2017), the availability of Mg²⁺ ions (Machel and Mountjoy, 1986), or a combination of these factors. To address these guiding hypotheses, a variety of sediment and microbial mat hand samples were collected from the SA and NA of GSL and these were subjected to mineralogical analyses. In addition, half meter sediment cores were collected from the SA and NA of GSL and these were subjected to mineralogical and porewater geochemical analyses. In parallel, rates of microbial inorganic carbon assimilation (primary production), microbial acetate mineralization and assimilation (secondary production), and quantities and compositions of archaeal and bacterial 16S rRNA gene transcripts were determined. Results are discussed in the context of the potential interplay of biological and abiological factors in the formation of carbonate minerals [i.e., aragonite, (proto)dolomite].

Materials and Methods

Sampling Locations

Four sediment cores were collected from GSL for detailed characterization: two from the NA and two from the SA (Figure 1). The NA cores were collected offshore and a few hundred meters northwest of Spiral Jetty (41.43783°, −112.67103°). The NA sampling site is a west facing shoreline. The SA cores were collected offshore of Antelope Island Marina (41.06117°, −112.24531°) in a north facing bay that is protected from the rest of the lake by a small point extending into GSL called Ladyfinger. These locations were selected for several reasons: (i) both sample locations are along the GSL Fault Zone and represent marginal lake facies, (ii) both sample locations have microbialites (Eardley, 1938), although the NA microbialites have previously been inferred to no longer be actively forming (Lindsay et al., 2017), and (iii) both sampling sites are near wave base. However, the SA site is partially sheltered from wave activity by Ladyfinger Point (Figure 1).

Two sediment hand samples (top 5–8 cm) were collected from a playa in the northern most extension of the NA in GSL, termed Spring Bay (41.664801°, −112.836629°), on 9 April 2019 (Figure 1). These sediments were chosen since they, unlike the other NA sampling location, are likely to be influenced by seasonal input of organic matter. However, despite being only approximately 3 m apart, they exhibited different textures, as described below. These samples, termed NA North 1 and 2 (NAN-1 and NAN-2, respectively), were kept cool after transport to the lab where they were stored at −20°C.

Benthic microbial mats were collected from water column depths of no more than 0.5 m in several locations in the SA; benthic mats do not form in the NA. These hand samples include mats collected from Antelope Island 24 m to 8 km from the sediment core sampling site (depicted as SA on map; 41.06132°, −112.24510°) and Stansbury Island (depicted as SI on map; 40.91731°, −112.49111°) on 15 and 16 August 2016, respectively (Figure 1). Mat samples labeled Windward West, Bridger Bay, Buffalo Point, and Ridges were collected from the Antelope Island sampling location. All mat samples were collected aseptically, frozen on site using dry ice, and kept at −80°C prior to analyses.

Sample Collection and Field-Based Analyses

The duplicate NA sediment cores were collected on 9 May 2016 and the SA cores on 12 May 2016. Cores were collected ∼3 m apart at each location, and each core was taken from an area of sediment submerged ∼10 cm beneath the surface of the water and ∼5 m from the shoreline. Cores were taken from areas interstitial to existing microbialite structures. PVC pipes (6.3 cm diameter) were driven into the sediment to a depth of 50 cm from the sediment–water interface. Cores were removed while minimizing sediment loss. One PVC pipe used at each sample site had pre-drilled holes stoppered with butyl rubber stoppers at 5 cm intervals. These cores were stored upright on ice within a cooler during transport to the laboratory where they were stored at 4°C and used to assess microbial activities within 1 week of collection (described below). The second core from each site was frozen upright with dry ice on site and remained frozen during transport to the lab, where it was transferred to a −80°C freezer. These cores were sectioned into 5 cm segments (i.e., pucks) in a −20°C facility using a sterilized band saw. Core pucks were then halved, with one half reserved for stratigraphic, mineralogical, and image analyses and the other half reserved for molecular microbiological and pore water chemical analyses. Half-pucks were stored in sterile Whirl-Pak sampling bags at −80°C.

Field measurements on water overlaying the core sampling sites were conducted prior to core collection. Temperature and pH were measured using a YSI pH100CC-01 temperature compensated pH meter (YSI Inc.). Conductivity and total dissolved solids (TDS) were measured using a temperature-compensated YSI EC300 conductivity meter. Field measurements were not made for waters overlaying NAN-1 and NAN-2 sediments or SA microbial mats.

Core Descriptions

Half-cores reserved for mineralogy were subsampled for lithologic description. Individual samples were cut in centimeter increments with an X-Acto knife. Cuttings were then examined using a Nikon SMZ645 Model C-DS optical microscope.

Mineralogical Analysis

The NA and SA sediment columns were subsampled with 1.1 mm diameter Durafilm Kapton Tubing at 1–3 cm intervals. Multiple samples were taken at each depth without drying the columns to minimize sampling bias and avoid further halite precipitation. Individual sample-containing Kapton tubes were sealed at both ends with playdough to prevent dehydration during XRD analysis. Samples were then analyzed with the Rigaku Rapid II X-ray diffraction system (Rigaku Corporation, Tokyo, Japan) with Mo Kα radiation. This XRD instrument uses a 2-D image-plate detector for signal collection and integration using Rigaku’s 2DP software. X-ray diffraction (XRD) experiments were run at 50 kV with a 100-μm diameter collimator. Mineral phases were identified using Jade 9.0 software. Mineral assemblages were calculated based on Rietveld refinement using TOPAS software with crystal structures from the American Mineralogist Crystal Database (AMCSD). Pearson VII peak were used for all refinements. Scanning electron microscopy (SEM) and energy-dispersive X-ray spectroscopy (EDS) were performed using a Hitachi S3400 (Hitachi, Ltd., Tokyo, Japan) with the Oxford EDS system (Oxford Instruments, Abingdon, United Kingdom). Samples for SEM were prepared by putting powder on carbon tape and then coating with 1 nm iridium.

For the purposes of this study, the CaMg(CO₃)₂ terminology laid out in Fang and Xu (2019) is used. To this end, the ordering state parameter s^∗ was calculated for each analyzed sample:

s^{*} = 2 x_{M g} - x_{M g - m a x}

Where x_Mg is the occupancy of Mg at Mg sites in the mineral lattice and x_Mg–max is the maximum occupancy allowed by the mol% MgCO₃ of the sample. Here, VHMC will refer to Ca-Mg carbonate with 25–36% MgCO₃ and no cation ordering (s^∗ = 0), and “protodolomite” will refer to Ca-Mg carbonate with > 36% MgCO₃ and weak-to-moderate cation ordering (0 < s^∗ < x_Mg–max). “Dolomite” will be reserved for reference to Ca-Mg carbonate with >36% MgCO₃ and strong cation ordering (s^∗ = x_Mg–max for Ca-rich dolomite, where x_Mg–max ≤ 1; s^∗ = x_Ca–max for Mg-rich dolomite, where x_Ca–max ≤ 1).

Total Organic Carbon (TOC) and Total Nitrogen (TN)

Subsamples of sediment hand samples and select core subsections were thawed at 4°C overnight and then subsamples (∼0.2 g of wet weight sediment) were incubated at 90°C for 12 h to remove excess water. Dried samples were homogenized with a mortar and pestle and fumed with 12 N HCl for 24 h in a hermetically sealed glass container. Samples were dried once more at 90°C for 12 h prior to transfer to individual tin cups, weighing, and analysis on a Thermo Finnigan FlashEA 1112 nitrogen and carbon analyzer (Macdonald et al., 2018). TOC and TON were normalized to grams dry weight sediment.

Sediment Porewater Characterization

Subsamples of ∼5 g were removed from each frozen core section using a scroll saw in a −20°C facility. The wet mass of one such sample was determined, and then the tube was uncapped and placed in a 90°C drying oven. After 48 h in the oven, the sample was weighed again to obtain its dry mass. Water content was calculated as percent of mass lost during drying.

pH and NaCl concentrations were measured in a second subsample of each core. Core sediments formed a thick mud when they were thawed making it difficult to separate a supernatant. Therefore, subsamples were added to 15 mL conical tubes containing 3 mL diH₂O. Thawed sediment slurries were mixed by vortexing, followed by centrifugation at 4000 × g for 15 min. The refractive index of one aliquot of each clarified supernatant was measured with a Reichert AR200 handheld digital refractometer (Reichert Technologies). Since GSL is primarily a NaCl lake (Oren, 2013), the refractive indices were converted to NaCl concentrations using a NaCl standard curve. Another aliquot from the supernatant described above was used to measure pH using an Accumet AB150 pH probe. Sulfate concentrations were determined using the barium sulfate turbidity method (Burlage et al., 1998). Sulfate standards were prepared using sodium sulfate decahydrate (Na₂SO₄ ⋅ 10H₂O) in solutions of NaCl at concentrations specific to each sediment core.

Sulfide concentrations were determined in a third subsample of each frozen core section added to a 50 mL conical tube. To minimize sulfide oxidation, tubes were flushed with 60 mL N₂ gas, capped, and stored in a −80°C freezer prior to analysis. To measure sulfide in sediment porewaters while minimizing turbidity that interferes with the measurements, tubes were brought into an anaerobic chamber where 5 mL anoxic water was added to each to thaw and suspend sediments. As soon as samples had completely thawed, they were centrifuged at 4600 × g for 5 min. Once the spin was complete, the concentration of total dissolved sulfide in the supernatant was immediately determined by the methylene blue reduction method (Fogo and Popowsky, 1949). Standard curves were generated with sodium sulfide (Na₂S) solutions made at NaCl concentrations specific to each sediment core. Sulfide concentration was calculated based on the water content of sediments and dilution by the added diH₂O. The amount of sulfide reported is likely a minimum estimate considering the potential for degassing and/or oxidation during the preparation of the sample and does not account for mineral sulfides that might have been present in the particulate fraction of the sediment cores.

Dissolved gases other than sulfide were measured using a fourth subsample added to a 60 mL syringe sealed with a stopcock, using previously described methods (Fendinger and Adams, 1986). Briefly, the headspace of the syringe was flushed with high purity nitrogen (N₂) gas five times before 10 mL anoxic zinc acetate solution was added (final concentration of 0.1 M) to precipitate any sulfide that might be present and to provide sufficient supernatant for equilibration of analytes between liquid and gas phases. A 10 mL volume of high purity N₂ gas was then added to the syringe. The subsample was thawed and was then shaken vigorously for 1 min by hand. The 10 mL gas volume was then removed and transferred to a pre-evacuated Cali Bond gas-sampling bag (Calibrated Instruments, Inc.) that had been flushed three times with high purity N₂ gas.

Five mL of gas from each Cali Bond gas sampling bag was injected into an SRI 8610C gas chromatograph (SRI instruments) as previously described (Lindsay et al., 2018). Hydrogen (H₂), carbon dioxide (CO₂), and methane (CH₄) were separated on a 4.5 m × 0.125” OD Hayesep DB 100/120 packed column (Valco Instruments Company, Inc.). Helium was used as the carrier gas. Gas sample concentrations were calculated using a standard curve (1–10,000 ppm calibration range) generated daily using certified (± 5%) standards (American Gas Group). Gas phase concentrations were converted to dissolved aqueous phase concentrations according to Henry’s law, as previously described (Spear et al., 2005). Concentrations were normalized to mean sediment pore water content for each core (35.4 and 13.4% water by mass for SA and NA, respectively).

Carbonate Saturation Indices

Saturation indices for aragonite and dolomite were calculated for both the NA and SA over a range of historical geochemical conditions in GSL using PHREEQc (Saini-Eidukat and Yahin, 1999; Parkhurst and Appelo, 2013). Historical GSL water column alkalinity measurements were collected from prior publications (Stephens and Gillespie, 1976; Post, 1977; Gwynn, 1980; Post, 1981; Jones et al., 2009; Pace et al., 2016) and the minimum, maximum, and average values were determined for use in modeling (Supplementary Table S1). Historical GSL water column chemical data were obtained from the USGS GSL Brine Chemistry Database (Rupke and Mcdonald, 2012) using the 30 May 2017 revision, and the minimum, maximum, and average values for select sites were determined for use in modeling (Supplementary Table S2).

The maximum, minimum, and average historical saturation indices in the NA and SA of GSL (Supplementary Table S3) were calculated with the maximum, minimum, and average alkalinity and bulk water chemistry values (Supplementary Tables S1, S2). Data from Utah Geological Survey sites LVG4 (−112.7616, 41.3240) in the NA and AS2 (−112.2558, 40.8349) in the SA were selected because they possessed the longest and most complete datasets over the period of 1966–2017. Saturation indices for aragonite and dolomite in the bulk water column were calculated at 32°C (historic maximum), 15°C (historic median), and 0°C (near freezing). The pH of the water column was fixed at 7.7 for LVG5 (NA) and 8.2 for AS2 (SA), as recorded previously (Post, 1977).

The median alkalinity for the NA or the SA was used along with the chemistry of overlying waters at Utah Geological Survey sampling sites SJ-1 (−112.6684, 41.4380) in the NA and AIS (−112.2544, 41.0146) in the SA on 1 July 2016 (Supplementary Table S4) to calculate aragonite and dolomite saturation indices as a function of bulk water column pH (Supplementary Table S5). These sites were selected for calculations of aragonite and dolomite saturation as a function of pH because they were located closest to our sampling sites and data were collected near to the time that samples analyzed in this study were collected. A surface water temperature of 26°C was invoked, as in a previous study (Pace et al., 2016), for calculations of aragonite and dolomite saturation indices.

Determination of Microbial Primary and Secondary Production/Activity in Core Sediments

Rates of dissolved inorganic carbon (DIC) assimilation and acetate assimilation/mineralization were determined as proxies for microbial primary production and secondary production/activity, respectively. Sediments were subsampled from the unfrozen SA core at the sediment-water interface and at settled depths scaling to 11.3, 22.5, 33.8, and 45.0 cm. Likewise, sediments were subsampled from the unfrozen NA core at the sediment–water interface and at depths scaling to 11.7, 23.3, 29.3, and 35.0 cm. Roughly 100 mg aliquots of sediment from each depth interval were added to sterile 24 mL serum bottles. After sediment addition, the bottles were capped with butyl rubber stoppers and purged with N₂ for ∼5 min. Sediments were then overlaid with 10 mL filter-sterilized, N₂ purged GSL water collected from the water column overlying each sampling location. Serum bottles and their contents were then degassed with N₂ for an additional 30 min. For each depth interval at each site, 15 replicate assays were prepared for measurement of DIC assimilation and 15 were prepared to measure acetate mineralization (see below).

Dissolved inorganic carbon assimilation rate potential assays were initiated by injecting 1.0 μCi of [¹⁴C] sodium bicarbonate (NaH¹⁴CO₃) to each serum bottle to achieve a final concentration of 2 μM NaH¹⁴CO₃. Acetate mineralization and assimilation rate potential assays were initiated by injecting 0.5 μCi of [¹⁴C] sodium acetate (¹⁴CH₃COONa) to each serum bottle to achieve a final concentration of 1 μM. Fifteen assays for each ¹⁴C source were set up for each depth interval and were incubated at temperatures corresponding to the overlying water at the time of sampling. Triplicate assays were removed following 1, 2, 4, 8, and 12 h of incubation at the specified temperature. The assays were then terminated by freezing at −20°C prior to analysis. Assays were processed for assimilation of inorganic carbon to biomass, assimilation of acetate to biomass, and oxidation of acetate to CO₂ using scintillation counting as previously described (Urschel et al., 2015). Radioactivity was then measured in counts per minute (CPM) on a Beckman LS 6500 liquid scintillation counter (Beckman Coulter, Inc., Indianapolis, IN, United States) and converted to disintegrations per minute (DPM) using a quench curve. The average and standard error of the mean of three replicate assimilation assays are presented for the pair of incubation time points that showed the largest difference in substrate transformation. Separate samples of the inoculum were dried and weighed to determine the dry mass of the sediment in the sediment slurry that was added to each assay.

RNA Extraction and cDNA Synthesis

Duplicate samples of sediment were subsampled from frozen core sections with a sterile spatula after first scraping off the top layer (∼3 mm). Samples were allowed to thaw before they were subjected to RNA extraction, DNAse treatment, quantification, and conversion to cDNA, as described previously (Lindsay et al., 2018). cDNA was subjected to amplification of archaeal, bacterial, and eukaryal small subunit rRNA gene transcripts using the primer sets 344F/915R, 1100F/1492R, and A7F/570R, respectively, as previously described (Hamilton et al., 2013). Quantitative PCR (qPCR) of 16S/18S rRNA gene transcripts was performed as previously described (Lindsay et al., 2018) using plasmid standards previously constructed using amplicons from GSL (Lindsay et al., 2017). Negative control reactions were performed in the absence of added cDNA.

Eukaryal amplicons were not detected in cDNA via PCR or qPCR. Archaeal and bacterial 16S cDNA amplicons were barcoded and sequenced with the 2 × 300 MiSeq Reagent Kit v3 following the manufacturer’s instructions, as described previously (Lindsay et al., 2018). Post sequence processing was performed with Mothur (ver. 1.36.1) (Schloss et al., 2009) as previously described (Hamilton et al., 2013) after merging the paired reads. Merging was only conducted for bacterial reads; archaeal reads were not merged due to poor overlap. Further details of how sequences were processed can be found in our previous publication (Lindsay et al., 2018). Following processing, a normalized size of 533 and 1809 16S cDNA sequences was obtained for each of the archaeal and bacterial libraries, respectively. Sequences specific for each operational taxonomic unit (OTU, classified at 97% sequence similarity) were classified using the Bayesian classifier of Mothur (Wang et al., 2007) and the RDP database (training set 14), with manual verification using BLASTn. The compositions of the forward and reverse reads for archaeal 16S rRNA gene transcript libraries were similar, and thus only the results of the forward reads are presented. Raw reads, quality scores, and mapping files for the archaeal and bacterial 16S rRNA gene transcript libraries have been deposited in the NCBI short reads archive under BioProject accession number PRJNA481529.

Results

Sediment Sample Characteristics

Two sediment hand samples were collected from Spring Bay in the NA of GSL, NAN-1 and NAN-2. NAN-1 consisted primarily of rounded oolitic sand and playa mud while NAN-2 was primarily playa mud. In addition, sediment cores were collected for use in characterizing the chemical, mineralogical, and physical attributes of NA and SA sediments.

The NA sediment core measured approximately 40 cm in length (Supplementary Figure S1). Most of this core comprised rounded ooids, 0.5 mm or less in diameter, and enveloped in calcium carbonate and halite. Ooids had a dark, spotted texture beneath a coating of halite and calcium carbonate. These sediments were grain supported, loosely cemented, and lacked any matrix. The pH, temperature, specific conductivity, and TDS content of the water overlying the NA sediment column at the time of sampling were 7.50, 22.4°C, 207.2 mS, and 134.9 g L^–1, respectively.

The SA sediment core measured approximately 45 cm in length (Supplementary Figure S1). The core comprised planar mm–cm scale beds of alternating dolomitic to calcareous carbonate mud, clay, and silt particles with abundant salts. Occasional zones of stylolites occurred at bedding boundaries and black, dendritic precipitates existed throughout. Sediments were matrix supported, well-cemented, and contained orange-brown ooids typically less than 0.1 mm in diameter. The pH, temperature, specific conductivity, and TDS content of water overlying the SA sediment column at the time of sampling were 7.17, 30.5°C, 176.3 mS, and 103.9 g L^–1, respectively.

Sediment Column and Microbial Mat Mineralogy

X-ray diffraction revealed very different mineral assemblages among the sediment cores, sediment hand samples, and microbial mats (Figures 2,3 and Supplementary Table S6). The NA sediment core was dominated by aragonite with varied amounts of VHMC, protodolomite, calcite, halite, and smectite, as was the NAN-1 hand sample (Figures 2B, 3A and Supplementary Table S6). In contrast, sediment from NAN-2, collected only a few meters away from NAN-1, was dominated by protodolomite (Supplementary Table S6). The SA sediment core contained predominantly protodolomite with moderate amounts of quartz and aragonite and trace amounts of calcite, halite, and smectite (Figures 2B, 3B). Microbial mat hand samples collected from the SA comprised mostly aragonite with halite and Mg-calcite (Figure 2A and Supplementary Table S6). None of the mat samples contained dolomite or even HMC. The protodolomite-dominated sediment samples (SA core and NAN-2) were each composed of very fine grain muds, while the aragonite-dominated samples (NA core and NAN-1) comprised micron-scale oolites.

FIGURE 2

Figure 2. (A) X-ray diffraction (XRD) patterns of microbial mats from the South Arm (SA) of Great Salt Lake (GSL). Mats are dominated by aragonite (A) and halite (H) with trace amount of calcite (C). (B) Representative XRD patterns of sediments from the North Arm (NA) and SA cores at indicated depths below the sediment–water interface. Thin red lines are refined results overlaid on black original data. Thin black lines below each pattern are residues from Rietveld refinement. NA sediments are dominated by aragonite, halite, and calcite, with trace amounts of high magnesian calcite/protodolomite (D). SA sediments are dominated by high magnesian calcite/protodolomite with small amounts of halite, calcite, and quartz (Q).

FIGURE 3

Figure 3. Mineralogical composition of sediment columns from the North Arm (NA, A) and South Arm (SA, B), as determined via X-ray diffraction (XRD) data analysis with Jade 9 software. Representative spectra are depicted in Figure 2.

Values of mineral lattice spacings (i.e., d₁₀₄ values) of the protodolomite also varied with sediment depth, ranging from 2.883 to 2.931 Å in the NA and from 2.887 to 2.905 Å in the SA. Lower boundaries of mol% MgCO₃ were calculated based on an empirical d₁₀₄ versus mol% MgCO₃ curve for complete ordered dolomite (Graf and Goldsmith, 1956; Reeder and Sheppard, 1984; Zhang et al., 2010; Fang and Xu, 2019). Based on this relationship, the NA core comprised at least 33 mol% MgCO₃ and the SA comprised at least 44% MgCO₃.

SEM

North Arm sediments were homogeneous under SEM throughout the 40 cm core and a substantial portion of grains had cubic halite crystals precipitated on their surfaces. These sediments contained well-rounded grains of ooids and pellets with crater-like textures (Figures 4a,b). Grain size ranged from 2–6 μm (methods and results described in the Supplementary Material) and EDS revealed grain composition to be aragonite with VHMC or protodolomite regions cemented within (Supplementary Figure S2).

FIGURE 4

Figure 4. Scanning electron microscopy (SEM) images of sediments from North Arm (NA) and South Arm (SA) cores. (a) Aragonite ooite from 1 cm depth in the NA core. (b) Increased zoom portion of a showing disordered dolomitic region in the aragonite ooite. (c,d) Euhedral protodolomite crystals from SA core, 0 and 24 cm in depth, respectively. Plus signs and labels indicate locations of energy dispersive X-ray spectroscopy (EDS) analyses, the results of which are presented in Supplementary Figures S2, S3. Abbreviations: Dol, dolomite; Ara, aragonite.

In contrast, SA core sediments exhibited euhedral crystals under SEM (Figures 4c,d), which were found to be consistent with aragonite and protodolomite by EDS analysis (Supplementary Figure S3). Rather than being cemented within oolites, protodolomite in the SA core formed micrites. These rhombohedral protodolomite grains were relatively free of accessory mineral growth on their crystal faces.

Dispersive X-ray spectroscopy results showed the majority of the protodolomite in each core to be Mg-rich. The NAN-1, NAN-2, SA-0 cm, and SA-24 cm samples had averages of 54.4, 47.9, 53.3, and 57.4 mol% MgCO₃, respectively. Together with the d₁₀₄ values acquired through XRD analysis, these compositions indicate a consistent ordering state of s^∗ = ∼0.25 across all samples, as determined by previously published methods (Fang and Xu, 2019).

Sediment Porewater Geochemistry

The water content of the NA sediment column was always lower than that of the SA column, regardless of depth (Supplementary Table S7). While water content varied little with depth in the NA column (11–16% by mass), it was highly variable in the SA column (∼20–50% by mass). The pH of sediment porewater was higher in the NA core than the SA core at all depth intervals sampled (Supplementary Table S7). Porewater pH ranged from 8.04 to 8.54 in the NA core and from 7.64 to 8.01 in the SA core. The pH gradually decreased in the NA between 0 and 30 cm in depth and then increased markedly at 35 and 40 cm. In contrast, SA porewater pH was largely consistent throughout the sediment core.

North Arm porewaters were saturated with respect to NaCl while salinity in the SA core ranged from 13.1 to 16.7% g mL^–1 (Supplementary Table S7). Listed values for the NA greater than the maximum solubility of NaCl resulted from the addition of deionized water to generate enough supernatant to measure with a digital refractometer. This addition allowed for dissolution of halite crystals (visible in the sediment core before manipulation).

Porewater dissolved gas profiles varied markedly as a function of depth both within and between sediment cores (Supplementary Table S7 and Figure 5). In the NA, total dissolved sulfide increased with depth from 6 μM at the sediment–water interface to 110 μM at 35 cm. Porewater CH₄ concentrations in the NA core were uniformly low along the depth transect and never exceeded 2 nM.

FIGURE 5

Figure 5. Select geochemical measurements juxtaposed with rate potentials of biological CO₂ assimilation and acetate transformation and the abundances of sulfate reducing bacteria (SRB) plotted by depth for the North Arm (NA, top) and South Arm (SA, bottom). Not shown are measured pH values, which were 8.3 ± 0.2 throughout the NA core and 7.9 ± 0.1 throughout the SA core (see Supplementary Table S7). Standard error for gene copy number measurements can be found in Supplementary Table S9. Abbreviations: GDWS: gram dry-weight sediment. DIC: dissolved inorganic carbon. DPM: disintegrations per minute. C:N: total organic carbon: total nitrogen ratio.

Total dissolved sulfide in SA sediment porewater was 500 μM at the sediment–water interface and decreased precipitously to below 4 μM at 10 cm and lower, with a spike of 65 μM at 35 cm. CH₄ was detected at concentrations ranging from 15 to 20 nM between the sediment water interface and 10 cm depth, there after falling quickly to sub-nM concentrations by 25 cm and remaining at low levels.

Concentrations of dissolved H₂ in sediment porewaters from both the NA and the SA were highly variable, with higher concentrations typically observed at depth when compared to surface transects (Supplementary Table S7 and Figure 5). In the NA, dissolved H₂ concentrations were similar along the depth profile and always below 75 nM. SA sediment porewater H₂ was present at ∼30 nM at the sediment–water interface and rose to ∼140 nM at 15 cm in depth. Between 15 and 20 cm porewater H₂ dropped to just 40 nM before rising gradually to over 250 nM at 40 cm.

The TOC to TN (C:N) ratio in sediments from the NA core did not fluctuate and remained ∼28 throughout (Figure 5). In contrast, C:N was low (4–5) in surface sediments of the SA core but then increased to ∼30 in deeper sediments. The increase in the C:N ratio with depth in the SA core was not due to large changes in TOC, but rather to decreases in TN (Supplementary Table S7). Overall, TOC was similar in the two cores, while TN was higher in the SA, especially at shallow depths. TOC and TN were 1.64 and 0.14%, respectively, in the NAN-1 hand sample and were 0.15 and 0.02% in the NAN-2 sample. The C:N ratio was 7 in NAN-1 and 11 in NAN-2 sediments.

Abundance and Affiliations of Bacterial 16S rRNA Gene Transcripts

Bacterial 16S rRNA gene transcript abundances were similar in the NA and SA and relatively consistent throughout the sediment columns (Supplementary Table S8). In the NA column, bacterial transcripts were highest in the 0 cm section, declining gradually with depth to about 30-fold lower in the deepest (35 cm) section. SA transcripts were also highest in the 0 cm section but varied by less than an order of magnitude throughout.

A total of 1809 bacterial 16S rRNA gene transcript sequences were subsampled from libraries from each core section (Supplementary Figure S4 and Supplementary Table S8). At the genus level, the bacterial communities in both cores at 10 cm depth and below were dominated by OTUs affiliated with Sphingomonas (Supplementary Figure S4). However, in the topmost sediments (0 and 5 cm core sections), Sphingomonas sequences were only minor components. The 0 and 5 cm core sections from the NA were dominated (83 and 61% of total reads, respectively) by OTUs that exhibited 93–96% sequence identity to sequences from uncultivated organisms affiliated with the phylum Acetothermia. A minor component (2.6%) of the OTUs from the NA 0 cm core section was affiliated (93–94% sequence identity) with the SRB Desulfonauticus autotrophicus. In the SA, OTUs accounting for 61 and 47% of the reads in the 0 and 5 cm core sections, respectively, exhibited 94–100% sequence identity to the SRB Desulfohalobium utahense.

We compiled OTUs that exhibited close affiliation with known SRBs and multiplied their relative abundances by the total abundance of bacterial transcripts at each depth with which that OTU was affiliated to determine the absolute abundance of templates affiliated with dominant OTUs (Supplementary Table S8 and Figure 5). In the NA core, the predominant SRBs were affiliated with Desulfonauticus (10⁶ copies gdws^–1) and, to a lesser extent, Desulfohalobium (10³ copies gdws^–1). Templates affiliated with Desulfonauticus were only identified in the top section (0 cm) of the sediment column. In the SA, the predominant SRBs were affiliated with Desulfohalobium (10⁶ copies gdws^–1) and to a lesser extent Desulfosalsimonas (10⁴ copies gdws^–1), and these were most abundant in the 0–5 cm core sections. Data for the abundance and composition of archaeal transcripts, which will not be discussed in detail here, can be found in Supplementary Figure S5 and Supplementary Table S8.

Rate Potentials of Microbial Primary and Secondary Production

Assimilation of DIC (primary production) was detected in communities associated with sediments from each NA and SA core section (Supplementary Table S9 and Figure 5). Rates of DIC assimilation were low and did not change significantly with depth in the NA sediment core. In the SA, potential rates of DIC assimilation were up to four times higher than in the NA. These were highest at the sediment–water interface and at 10.0 cm depth, nearly double the rates measured in the bottom three depth intervals (22.5, 33.8, and 45.0 cm). As such, rates of DIC assimilation were inversely correlated with depth (Pearson R = −0.86).

Assimilation and mineralization of acetate (secondary production) were also detected in communities associated with sediments from each assayed core section (Supplementary Table S9 and Figure 5). In the NA core, acetate mineralization was lowest at the sediment–water interface and increased with depth by several orders of magnitude, with the highest rate of acetate mineralization detected at a depth of 29.2 cm and the highest rate of acetate assimilation measured at a depth of 35.0 cm. Thus NA sediments exhibited positive relationships between acetate mineralization/assimilation and depth (Pearson R = 0.85 and 0.64, respectively). Contrastingly, in the SA sediment core, rates of acetate mineralization and assimilation were highest in the top two depth intervals sampled (the sediment–water interface and at 10 cm depth). Rates then decreased with depth, with the lowest rate of mineralization measured at 40 cm and the lowest rate of assimilation at 30 cm. SA sediments thus exhibited inverse relationships between acetate mineralization/assimilation and depth (Pearson R = −0.95 and −0.94, respectively).

Aragonite and Dolomite Saturation Indices

The saturation indices of aragonite and dolomite were calculated based on lake conditions near sediment core collection sites. Figure 6A shows indices calculated using geochemical data collected July 1, 2016. Saturation indices were calculated as functions of pH across a range meant to capture conditions that might develop on a microscale in sediment porewaters.

FIGURE 6

Figure 6. (A) Calculated saturation states of aragonite and dolomite as functions of pH in the North Arm (NA) and South Arm (SA) of Great Salt Lake (GSL) as modeled by PHREEQc using chemical data from the 1966–2016 GSL Brine Database (Rupke and Mcdonald, 2012). (B) Historical maximum, minimum, and average saturation states of aragonite and dolomite calculated at maximum (32°C), minimum (0°C), and average (15°C) temperatures for the NA and SA of GSL modeled by PHREEQc using bulk water column solute data from the 1966–2016 GSL Brine Database (Rupke and Mcdonald, 2012). In A, note that aragonite is undersaturated (i.e., aragonite dissolution is favored) with respect to aqueous solute composition in solutions with pH < 7.0 in both the NA and SA, while dolomite remains saturated (i.e., dolomite precipitation is favored) in solutions with pH > 5.5 in the NA and pH > 6.0 in the SA. In B, note that aragonite only becomes undersaturated under extreme conditions (minimum alkalinity and low temperatures), while dolomite remains saturated under all observed conditions. Values plotted in A are presented in Supplementary Table S5 and values plotted in B are presented in Supplementary Table S3.

These calculations indicate that GSL was saturated with respect to both aragonite and dolomite (favoring mineral formation) through the whole pH range (7.64–8.54) observed in bulk NA and SA porewaters. However, saturation indices for both minerals decline with decreasing pH. The saturation index of aragonite falls below zero (undersaturation, favoring mineral dissolution) at pH < 6.7 in the NA and < 7.1 in the SA. The dolomite saturation index falls below zero at pH < 5.5 in the NA and < 6.0 in the SA. Thus, in the pH range 5.5–6.7 (NA) or 6.0–7.1 (SA), dolomite is thermodynamically stable while aragonite is predicted to dissolve.

To encompass historical variation in the geochemical conditions of GSL and to assess its effect on aragonite and dolomite saturation indices, data were compiled from locations in the NA and SA over the time frame 1966–2017 (Figure 6B). Historical ranges of water solute concentrations were used to calculate a minimum, average, and maximum saturation index for each mineral at the minimum, average, and maximum recorded temperatures. These calculations indicate that geochemical conditions have favored dolomite formation throughout this time frame in both the NA and SA. Conditions have also likely favored aragonite formation throughout this period, with its saturation index falling below zero only at the lowest recorded alkalinities and during the coldest intervals.

Discussion

Proposed Model for Protodolomite Formation in GSL

A mechanism has been previously proposed to explain (proto)dolomite formation in phototrophic microbial mats in Bridger Bay, GSL (Pace et al., 2016). In this model, the activity of SRB drives aragonite precipitation within the mat structure. This biomass-associated aragonite is then deposited to sediments where dolomitization of aragonite takes place due to localized increases in acidity induced by fermentative processes. However, 16S rRNA gene (Lindsay et al., 2017, 2019) and metagenomic sequence (E.S. Boyd, unpublished data) do not reveal abundant populations inferred to be SRBs (< 2% of total community) in benthic phototrophic microbial mats associated with microbialites from Bridger Bay. We therefore propose a modification to the model put forth by Pace et al. (2016) which may also help to explain the distribution of aragonite and protodolomite in GSL (see Figure 7). Like the previous model, we suggest aragonite precipitation takes place in photosynthetic microbial mats due to a localized increase in alkalinity and this aragonite is then deposited into benthic sediments. However, the increase in alkalinity is proposed to be induced by photosynthetic CO₂ consumption rather than the activity of low abundance SRBs. This interpretation is consistent with previous studies showing that oxygenic photosynthetic organisms can precipitate aragonite through localized increases in alkalinity (Merz, 1992; Al-Horani et al., 2003; Kranz et al., 2010).

FIGURE 7

Figure 7. Proposed model for the formation of aragonite and dolomite in mats and sediments in Great Salt Lake (GSL). Photosynthesis in benthic cyanobacterial mats generates a local increase in alkalinity, promoting the precipitation of aragonite on extracellular polymeric substances, as well as generating large quantities of biomass and O₂. Aerobic heterotrophs oxidize sufficient biomass and photosynthate to deplete local O₂, generating a niche for anaerobic fermentative bacteria and sulfate-reducing bacteria (SRB). Organic acids produced by fermenters both lower porewater pH, promoting aragonite dissolution, and provide substrates supporting heterotrophic SRB in underlying sediments. SRB activity promotes dolomite precipitation in mildly acidic microenvironments where aragonite precipitation is thermodynamically disfavored. Abbreviations: Ar, aragonite; Do, dolomite; Cyano, cyanobacteria; Het, aerobic heterotrophic bacteria; Ferm, fermentative bacteria; SRB, sulfate-reducing bacteria.

Also like the previous model, we suggest that partial dissolution of aragonite occurs as a result of the activities of fermentative organisms inhabiting benthic mats and sediments. These organisms produce organic acids as waste, lowering local pH to a level that favors dissolution of aragonite, but not dolomite (Figures 6A, 7). Based on available data, it is unclear whether fermentation takes place in the mats or the sediment column. However, numerous species of Clostridia, a genus commonly known for fermentation (Rainey et al., 2015), have been detected in previous molecular characterizations of benthic cyanobacterial mats from the SA of GSL (Lindsay et al., 2017, 2019). It has been suggested that Mg²⁺ ions also need to be released to allow for dolomite to form (Pace et al., 2016). However, the Mg²⁺:Ca²⁺ ratios of waters in both the SA and NA do vary seasonally and have been > 20 since 1966 (Rupke and Mcdonald, 2012), which should be favorable for (proto)dolomite formation (Folk and Land, 1975).

Finally, in the slightly acidic microenvironments created by the activities of fermentative bacteria, where aragonite precipitation is thermodynamically unfavorable, the activities of SRBs promote the formation of VHMC and protodolomite (Pace et al., 2016). We suggest that the dolomite in the SA core is likely authigenic and primary, since it has never experienced deep burial and the transformation of aragonite to dolomite is sluggish at temperatures common to GSL (0–32°C; Figure 6B) due to kinetic constraints (Machel and Mountjoy, 1986). This is supported by SEM images from the SA sediment core (Figure 4) that show that the protodolomite is euhedral and displays no aragonite pseudomorph, features that suggest it is not a dolomitization product of aragonite.

Potential Controls on the Formation of Protodolomite in GSL

While elucidating the precise mechanisms of protodolomite precipitation is beyond the scope of this study, some inferences can be made from the data presented. A commonly proposed model for dolomite formation in brines rich in MgSO₄ involves depletion of SO₄^2– by SRB, which is thought to reduce ion pairing with Mg²⁺ and make the latter more readily available for incorporation into carbonate lattices (Baker and Kastner, 1981; Vasconcelos and Mckenzie, 1997; Warthmann et al., 2000; Bontognali et al., 2010). However, our SA sediment core, which is rich in protodolomite, exhibits SO₄^2– concentrations exceeding 40 mM throughout (Supplementary Table S7 and Figure 5), indicating that substantial depletion of SO₄^2– via SR activity is unlikely to be responsible for protodolomite formation in this system.

Alternatively, most current hypotheses for sedimentary dolomite formation require the presence of a catalyst such as dioxane (Oomori and Kitano, 1987), methane (Xu, 2010; Zhang et al., 2010), polysaccharides (Zhang et al., 2012b; Shen et al., 2015), EPS (Shen et al., 2015; Huang et al., 2019), dissolved sulfide (Zhang et al., 2012a, 2013; Shen et al., 2014; Lu et al., 2018), or clays (Liu et al., 2019). The presence of a catalyst may help lower the dehydration energy of Mg²⁺-water complexes and allow crystal nucleation and growth (Lippmann, 1973; De Leeuw and Parker, 2001; Gregg et al., 2015). Importantly, the lack of (proto)dolomite in all of our microbial mat samples suggests polysaccharides and EPS are unlikely to control micritic dolomite precipitation in GSL, as these molecules are abundant in such mats (Bolhuis et al., 2014). The clay smectite was detected in both the NA and SA sediment column while protodolomite was largely restricted to the SA core. This set of observations discounts the probability of clays as a potential catalyst. Another potential catalyst is dissolved sulfide, present in high concentrations (∼500 μM) in the top most sections of the SA core. Previous studies demonstrate that dissolved sulfide can promote Mg²⁺ incorporation and dolomite nucleation (Zhang et al., 2012a, 2013; Shen et al., 2014).

Sulfide concentrations of 500 μM, as found in the top section of the SA sediment core (Figure 5), are sufficient to lower the dehydration energy of Mg²⁺-water complexes and promote fast nucleation of HMC in Mg-rich solutions (Zhang et al., 2013). The zone of elevated sulfide in the SA core corresponds with a high abundance of 16S rRNA gene transcripts affiliated with the SRB D. utahense (98% nucleotide identity), previously isolated from the NA of GSL (Jakobsen et al., 2006) (Figure 5). D. utahense is a heterotroph that can use a variety of organic substrates to support SR, including the fermentation product acetate when H₂ is also provided (Jakobsen et al., 2006). Consistent with an abundance of D. utahense transcripts in the top section of the SA core, rates of acetate assimilation/mineralization were highest in the top section of the SA core and dissolved H₂ concentrations increased with depth in the SA core, suggesting a net sink for H₂ in the near surface sediments (Figure 5). Finally, the concentration of sulfide in the SA core was positively correlated with the abundance of transcripts attributable to D. utahense (Pearson R = 0.99, p < 0.01). This is consistent with other hypersaline intertidal areas, where active zones of SR are localized near the sediment–water interface (Van Lith et al., 2002). Together, these observations strongly suggest the presence of an active SR zone near the SA sediment–water interface, where protodolomite was detected in appreciable quantities (≥20% of sediment) despite being absent from overlying mats (Figure 2).

In contrast to the SA, evidence for an active SR zone in the NA core is lacking and protodolomite is scarce (Figure 2B and Supplementary Table S6). The primary SRB in the NA was related to D. autotrophicus, which requires H₂ or formate as an electron donor, but can grow either heterotrophically or autotrophically (Mayilraj et al., 2009). Despite the presence of this putatively autotrophic and hydrogenotrophic SRB, variance in sulfide concentration and CO₂ assimilation activity with depth was not observed in NA sediment porewaters (Figure 5). Moreover, the abundance of transcripts attributable to D. autotrophicus was not correlated with the concentration of sulfide and these sediments did not appear to be a sink for H₂ (Figure 5). It is possible that these cells are minimally active and are in a maintenance mode of metabolism due to salinity stress, a conclusion that would be consistent with previously compiled data on the activity and abundance of SRB in the NA versus the SA (Brandt et al., 2001; Boyd et al., 2017). In any case, the lack of an active SR zone in the NA corresponds with the relative scarcity of protodolomite in the NA sediment core.

Thus, even though the mechanisms of protodolomite precipitation remain unclear, we infer that the activity of SRB is at least one important control on this process in GSL sediments. The biotic and abiotic controls on the activities of SRB may then be expected to affect the formation of protodolomite as well. One such control that may be important is the abundance of nutrient-rich organic carbon, and by extension the activity of benthic photosynthetic communities. Benthic phototrophic microbial mats are ubiquitous in the SA of GSL, either associated with microbialites or in shallow marginal areas of the lake without microbialites. Overall, these mats contribute up to 30% of the total productivity in the SA (Wurtsbaugh et al., 2011). In contrast to the SA, phototrophic mats are conspicuously absent in the NA of GSL (Lindsay et al., 2017). Moreover, the total biomass associated with microbialites sampled from the SA (as estimated by ribosomal RNA gene abundances) was three orders of magnitude greater than in the NA. These observations were attributed to constraints imposed by high salinity on the primary architect and the primary producer of these mats, Cyanobacteria (Wurtsbaugh et al., 2011) within the genus Euhalothece (Lindsay et al., 2017, 2019) which grow sub-optimally or not at all at salinities that exceed ∼20% (Garcia-Pichel et al., 1998).

Not only the abundance of organic matter but also its quality may be important in promoting or constraining the growth of SRB. Organic matter in the top-most sediments in the SA sediment core had a C:N ratio of ∼4 (Figure 5), which is consistent with fresh photosynthetic biomass. In contrast, the C:N ratio of organic matter in the NA core was ∼25 throughout. High C:N ratios have been suggested to indicate input of terrestrial soil or vascular plant material (Prahl et al., 1994). However, since these ratios were not observed in the SA sediment core, where most riverine sediment input takes place, we suggest that in the NA they more likely reflect substantial processing by microorganisms and slow-to-negligible fresh biomass deposition. Thus, salinity may place constraints on the abundance of primary producers and their activity, which in turn influences production of biomass that can support heterotrophic activity, such as fermentation and SR. High salinities also inhibit the growth of D. utahense, which grows optimally at a salinity of 8-10% NaCl, near to that of SA sediment porewaters (13.1–16.5%). Differences in salinity between the NA and SA of GSL may thus contribute to differential growth of both SRB and the Cyanobacteria on which they depend for organic carbon.

Unlike the NA sediment core, both the NAN-1 and NAN-2 hand samples were high in nitrogen-rich organic matter (C:N ratios of 7 and 11). Given the absence of benthic phototrophic mats at this sampling location and in the NA in general, this organic matter is interpreted to represent allochthonous input from down slope surface runoff to this playa location. This suggests that sediment characteristics imposed by depositional environment may represent an additional control on the formation of protodolomite. One such characteristic of both the NA core and NAN-1 sediments that provides information about their depositional environment is the abundance of ooilitic sand. Oolitic sands tend to accumulate in wave agitated waters (Pedone and Norgauer, 2002). Together with hypersalinity, high-energy (e.g., wave agitated) depositional environments may physically constrain the development of robust phototrophic mats and SR zones capable of promoting protodolomite formation. Moreover, oolitic sediments such as those identified in NAN-1 and the NA sediment core also tend to have high porosity (Choquette and Pray, 1970). This could more easily allow infusion of oxic waters to greater sediment column depths, which would inhibit the development of an SR zone due to the O₂ sensitivity of SRB (Postgate, 1979). The contrasting sediment compositions of NAN-1 and NAN-2 indicate that slight differences in location can lead to substantial differences in wave energy and sediment characteristics. Moreover, the presence of protodolomite in the NAN-2 sample suggests these differences may be as important as salinity in constraining the formation of protodolomite, given that salinity is expected to be relatively uniform on the scale (meters) that separated the NAN-1 and NAN-2 samples. Further investigation of the distribution of protodolomite and the microbiology and porewater geochemistry of sediments from a range of depositional settings in GSL may help in further evaluating this hypothesis.

Conclusion

The detection of aragonite as the primary carbonate mineral in photosynthetic mats from the SA of GSL, without the presence of abundant SRB, suggests that oxygenic photosynthesis and corresponding alkalinity increases are responsible for aragonite formation. The detection of protodolomite rather than aragonite as the primary carbonate mineral in sediment cores from the SA suggests primary protodolomite precipitation in benthic sediments, where the lowest organic C:N ratios (indicative of deposition of fresh photosynthetic biomass), the highest rates of acetate transformation, the highest abundance of SRB-associated gene transcripts, and the highest concentrations of sulfide occur. This suggests that deposition of photosynthetic biomass through a zone of active heterotrophic SR in underlying sediments is likely to be a primary variable favoring the formation of protodolomite. Following the model of Pace et al. (2016) and consistent with the general function of microbial mat ecosystems, the deposited organic matter is first oxidized and/or fermented, generating anoxia, organic acids, and localized zones of acidity. This leads to two outcomes, partial dissolution of detrital aragonite and stimulation of SRBs, that together promote protodolomite formation. The comparative lack of protodolomite in the NA sediment core is suggested to result from constraints imposed by extreme hypersalinity and depositional environment that together limit the development of productive phototrophic mats and, in turn, suppress the activities of fermenters and SRBs. Importantly, the detection of protodolomite in the NAN-2 hand sample, which had a high abundance of fresh organic carbon, suggests that the photosynthate that drives the development of an SR zone and that is putatively responsible for the formation of this mineral may not need to be autochthonous.

The protodolomite identified in the SA sediment core exhibits euhedral microtextures nucleating on aragonite surfaces. While SRB were present in the uppermost segments of the NA core, high sulfide concentrations and similar euhedral dolomite microtextures were not observed, suggesting that the presence of SRB alone is insufficient to promote formation of these mineral structures, at least at low abundances such as those observed in the NA core. We suggest that additional key factors promoting protodolomite formation arise due to the productivity of photosynthetic benthic mats, including heterotrophy-driven anoxia, fermentation-derived acidity decreasing the saturation index of aragonite, and heterotrophic SR with accompanying high concentrations of dissolved sulfide. Importantly, such conditions favorable for protodolomite formation may not develop in high energy environments due to the potential to physically remove photosynthetic mats or to create conditions (e.g., high porosity oolitic sands) that otherwise do not favor the development of an SR zone.

Data Availability Statement

Raw reads, quality scores, and mapping files for the archaeal and bacterial 16S rRNA gene transcript libraries have been deposited in the NCBI Short Reads Archive under BioProject accession number PRJNA481529.

Author Contributions

ED, EF, YF, ML, CS, NF, MW, AW, and DC collected data. DM and HX oversaw XRD analyses. BB and AR assisted with sample collection. DL and EB oversaw the project. All authors contributed to data analysis and interpretation. ED, EF, and EB wrote the manuscript.

Funding

The authors thank the Head of the Department of Microbiology and Immunology at Montana State University, Mark Jutila, for supporting this Precambrian Biosphere graduate class project. The authors acknowledge partial support for this project from the NASA Astrobiology Institute grants NNA15BB02A (EB) and NNA13AA94A (EB, HX, and YF). ED acknowledges support from a National Science Foundation Graduate Research Fellowship. ML thanks the Doyle Stephens Scholarship Fund and the Friends of Great Salt Lake for partial support of this work. ML also acknowledges support from the NASA Earth and Space Sciences Fellowship Program.

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2020.00024/full#supplementary-material

References

Alderman, A. R. (1958). Aspects of carbonate sedimentation. J. Geol. Soc. Aust. 6, 1–10. doi: 10.1080/00167615808728491