Introduction
Authigenic phyllosilicates commonly form in continental aqueous sedimentary environments characterized by alkaline-saline conditions developed under arid climates either by neoformation or transformation (Meunier Reference Meunier2005; Bristow and Milliken Reference Bristow and Milliken2011; Deocampo Reference Deocampo, Larson, Egenhoff and Fishman2015). Concentration of these lake waters by evaporation leads to a complex series of precipitation reactions as the brines evolve and control the mineral assemblages that are present. Phyllosilicates forming from direct precipitation in alkaline lakes tend to be Mg-rich minerals such as sepiolite, stevensite, kerolite, saponite, and hectorite (Jones Reference Jones and Mumpton1986; Meunier Reference Meunier2005). The major geochemical factors controlling the phyllosilicate minerals that precipitate may include: pH, alkalinity, element ratios, the presence or absence of detrital material, dissolved CO2 from biological activity, salinity, and the activities of Mg2+, Si4+, Al3+, Ca2+, and Fe3+ in solution (Jones and Galán Reference Jones, Galán and Bailey1988; Weaver Reference Weaver1989; Deocampo Reference Deocampo2005; Bristow and Milliken Reference Bristow and Milliken2011). The evolutionary pathway followed by brines evaporating in closed-basin lakes depends on the concentrations of HCO3 – and Ca2+ during calcite precipitation (Hardie and Eugster Reference Hardie and Eugster1970). During this critical first step, solutions evolve during calcite precipitation towards bicarbonate-rich brines if [HCO3 – ] > [Ca2+] and towards sulfate- and/or chloride-rich brines if [HCO3 – ] < [Ca2+] (Hardie and Eugster Reference Hardie and Eugster1970).
The authigenic phyllosilicates formed in waters that evolve toward bicarbonate-rich brines have been studied extensively in Lake Chad, Africa (Gac et al. Reference Gac, Droubi, Fritz and Tardy1977); Lake Yao, Africa (Darragi and Tardy Reference Darragi and Tardy1987); Abert Lake, Oregon, USA (Jones and Weir Reference Jones and Weir1983; Banfield et al. Reference Banfield, Jones and Veblen1991); Searles Lake, California, USA (Hay et al. Reference Hay, Guldman, Matthews, L&er, Duffin and Kyser1991); Jurassic Lake T’oo’dichi, Colorado, USA (Turner and Fishman Reference Turner and Fishman1991); the Pleistocene Amboseli basin, Kenya (Stoessell and Hay Reference Stoessell and Hay1978); the Miocene Vallecas-Vicálvaro and Cabañas-Yunclillos deposits, Spain (Galán and Castillo Reference Galán, Castillo, Singer and Galán1984); and the Miocene Eskişehir basin, Turkey (Ece and Çoban Reference Ece and Çoban1994; Kadir et al. Reference Kadir, Erkoyun, Erin, Hugettt and Önalgil2016). Jones and Galán (Reference Jones, Galán and Bailey1988) summarized these bicarbonate-rich systems with regard to silicate minerals outlining the conditions that favor the formation of certain Mg-rich clays over others during evaporative concentration. Brackish waters with a neutral to slightly alkaline pH and high Si activity favor the precipitation of palygorskite in detritus-rich environments. Examples include alluvial-lacustrine margins and sepiolite where reactive detritus is limited and Al is immobilized, such as in calcareous soils and near springs (Jones and Galán Reference Jones, Galán and Bailey1988; Birsoy Reference Birsoy2002; Bristow and Milliken Reference Bristow and Milliken2011). With increasing pH, decreasing activity of silica, and/or high concentrations of Mg2+ in solution, trioctahedral smectites form where the mineralogy is determined by the amount of detrital clay and geochemistry of the solutions (Badaut and Risacher Reference Badaut and Risacher1983; Calvo et al. Reference Calvo, Blanc-Valleron, Rodrigue-Ar&ia, Rouchy, Sanz, Thiry and Simon-Coincon1999; Larsen Reference Larsen2008). Lastly, in the most saline conditions near basin centers and areas with repeated wetting and drying cycles, high K+/Mg2+ activity ratios and high pH result in K uptake into smectite interlayers and the formation of illite-like domains as well as the neoformation of 10 Å clays (Turner and Fishman Reference Turner and Fishman1991; Hay and Kyser Reference Hay and Kyser2001; Hover and Ashley Reference Hover and Ashley2003; Furquim et al. Reference Furquim, Barbiéro, Graham, Neto, Ferreira and Furian2010).
The other half of the geochemical divide produces sulfate- and chloride-rich brines. For these systems the pH of the evaporating brines can be neutral to slightly alkaline or can be acidic. Moderately alkaline, sulfate-rich salars have been found on the Bolivian Altiplano where Mg-rich smectites occur if the pH is >8.5 (Badaut and Risacher Reference Badaut and Risacher1983). For the sulfate-alkaline systems, many of the factors controlling authigenic clay formation are already described by the high pH, bicarbonate-rich systems. Conversely, acid-saline, sulfate-rich lakes have been analyzed in south Western Australia where authigenic kaolinite and Fe oxides are favored to precipitate (Bowen and Benison Reference Bowen and Benison2009; Story et al. Reference Story, Bowen, Benison and Schulze2010). These sulfate/chloride systems are high-priority analogs for understanding past environments on Mars as many layered deposits on the planet have been found with both phyllosilicates and sulfates such as Gale Crater, Mawth Vallis, and Columbus crater, all of which could have hosted paleocrater lakes (Farrand et al. Reference Farr, Glotch, Rice, Hurowitz and Swayze2009; Milliken et al. Reference Milliken, Grotzinger and Thomson2010; Wray et al. Reference Wray, Squyres, Roach, Bishop, Mustard and Dobrea2010, Reference Wray, Milliken, Dundas, Swayze, rews-Hanna, Baldridge, Chojnacki, Bishop, Ehlmann, Murchie, Clark, Seelos, Tornabene and Squyres2011). Unfortunately, modern terrestrial environments with sulfate- and chloride-rich waters are relatively rare and have not been studied extensively (Bristow and Milliken Reference Bristow and Milliken2011). The deposits that have been found should, therefore, be analyzed in greater detail.
While the aforementioned generalizations can be helpful in establishing geochemical controls for clay-mineral formation, the complex environments created by evaporating brines across basins can make predicting clay mineral assemblages difficult. Many of the environmental conditions that favor magnesium clay formation overlap for sepiolite and Mg-smectite, especially for sulfate/chloride-rich systems in closed-basin lakes. For example, both sepiolite and Mg-smectite are favored at intermediate pH (8–9.5), intermediate alkali salinity, and groundwater-dominated input (Jones and Galán Reference Jones, Galán and Bailey1988; Galán and Pozo Reference Galán, Pozo, Galán and Singer2011). One of the main differences in the factors of formation between these two clays is the activity of silica and Mg in solution; where sepiolite is favored at high silica activities (Mg/Si ≤1 or log[aSiO2(aq)] > –4.5) and Mg-smectite at low silica activities (Mg/Si ≥1 or log[aSiO2(aq)] < –4.5) (Jones and Galán Reference Jones, Galán and Bailey1988; Galán and Pozo Reference Galán, Pozo, Galán and Singer2011; Tosca and Masterson Reference Tosca and Masterson2014; Tosca Reference Tosca, Pozo and Galan2015). For more complex systems, equilibrium activity diagrams for the MgO-CaO-Al2O3-SiO2-H2O-CO2-HCl system calculated by Birsoy (Reference Birsoy2002) found sepiolite to be favored when the pH is between 8.5–9.5 and log[aH2SiO4] > –4.75. Therefore, understanding the controls of dissolved silica in solution can aid in explaining Mg-clay formation in these intermediate environments such as the salars on the Bolivian Altiplano.
Previous analysis of the surface waters suggests the sulfate-rich brines for several of the salars are brackish with pH values ranging from neutral to alkaline (7–10) and saturated with respect to dissolved Si — where the log [aH4SiO4] is > –3.5 and well within the zone of sepiolite formation (Badaut and Risacher Reference Badaut and Risacher1983; Servant-Vildary and Roux Reference Servant-Vildary and Roux1990; Risacher and Fritz Reference Risacher and Fritz1991; Dejoux Reference Dejoux1993). Characterization of the clay minerals present in the surface waters by Badaut and Risacher (Reference Badaut and Risacher1983), however, found a highly soluble authigenic Mg-smectite resembling stevensite replacing the diatom frustules. This authigenic Mg-smectite formed only in salars where the pH of the surface water was >8.2, regardless of salinity (Badaut and Risacher Reference Badaut and Risacher1983). Diatom frustules have been cited as a source of silica for authigenic clay minerals and as a buffer for maintaining dissolved silica concentrations (Hay Reference Hay1970; Stoessell and Hay Reference Stoessell and Hay1978; Badaut and Risacher Reference Badaut and Risacher1983; Banfield et al. Reference Banfield, Jones and Veblen1991; Chahi et al. Reference Chahi, Fritz, Duplay, Weber and Lucas1997; Deocampo Reference Deocampo2005; Bristow and Milliken Reference Bristow and Milliken2011; Huggett et al. Reference Huggett, Cuadros, Gale, Wray and Adetunji2016). The overall effect of the algae on water chemistry with respect to clay mineral assemblages has not been examined thoroughly, however.
The goal of the present study was to examine further the salars of the Bolivian Altiplano extending out into the mudflats to aid in understanding the processes of phyllosilicate formation in sulfate-rich environments related to dissolved silica. The study included the analysis of the mudflat sediments from six salars located in the Sur Lipez region of the Bolivian Altiplano. These salars provide a spectrum of varying water chemistry, mineralogy, pH, and amount of detritus.
Geological Setting
The Bolivian Altiplano is a large intermontane basin (200,000 km2) with an average elevation of 3,650 m bounded by the Cordillera Occidental to the west and Cordillera Oriental to the east (Risacher and Fritz Reference Risacher and Fritz1991; Allmendinger et al. Reference Allmendinger, Jordan, Kay and Isacks1997). The climate is classified as arid, with a mean annual precipitation of <0.2 m y–1 and a mean evaporation rate of 1.5 m y–1 (Badaut and Risacher Reference Badaut and Risacher1983; Argollo and Mourguiart Reference Argollo and Mourguiart2000). The annual air temperature for the plateau is estimated to be between 5 and 10°C; the daily range may be up to 40°C, with extreme seasonal temperatures reaching –40°C, however in the winter and +25°C in the summer (Risacher and Fritz Reference Risacher and Fritz1991; Cabrol et al. Reference Cabrol, Grin, Chong, Minkley, Hock, Yu, Bebout, Fleming, Häder, Demergasso, Gibson, Escudero, Dorador, Lim, Woosley, Morris, Tambly, Gaete, Galvez, Smith, Uskin-Peate, Salazar, Dawidowicz and Majerowicz2009). The salars are remnants of large paleolakes that occupied the intravolcanic basin including the Pleistocene Lake Ballivián, Paleolake Minchín (35,000 years BP), and Paleolake Tauca (15,000 years BP) (Argollo and Mourguiart Reference Argollo and Mourguiart2000). The salars included in the current study occupy a number of small basins in the Sur Lipez region surrounded by stratovolcanoes, andesite lava flows, and ignimbrites (Fig. 1). The six salars analyzed in detail for the current study were selected based on overlap with previous research on the Altiplano. The salars chosen were: Laguna Canapa (CAN), Laguna Hedionda (HED), Laguna Ramaditas (RAM), Laguna Honda (HON), Laguna Chiar Khota (CHKT), and Laguna Chulluncani (CHU) (Fig. 1). Detailed descriptions of the salars can be found in Servant-Vildary and Roux (Reference Servant-Vildary and Roux1990) and Risacher and Fritz (Reference Risacher and Fritz1991).
Materials and Methods
Sampling
The salars were sampled in June of 2014 along transects extending from the outer boundary of the mudflat to the sediment–water interface under the surface water. For each sampling site, the sediments were excavated to the water table or ice layer (max 80 cm), with each sedimentary layer collected separately. Depending on the size, dryness, stability, and heterogeneity of sediments across the mudflat, the number of pits at each salar varied. Several dry depressions and slopes surrounding the basins were also sampled for comparison between pedogenic weathering and lacustrine processes for clay-mineral formation. Surface water samples were collected for each salar. The location, pH, and salinity of the surface waters from each lake are included in Table 1. These data also include salinity and pH measurements from previous studies in 1983, 1990, 1991, and 1993. Note that the surface waters of the salars were covered with ice each morning.
X-ray Diffraction
Subsamples were air-dried and ground with a mortar and pestle to pass through a 0.125 mm sieve for bulk X-ray diffraction (XRD) analysis. Due to the soluble and fine-grained nature of the samples, the micronizing mill was not used. Corundum (Al2O3) was mixed with each sample (10 wt.%) as an internal standard to enable estimation of 'X-ray amorphous' phases in randomly oriented powders. For clay mineral analysis, the <2 μm fraction was separated from each sample (wet-sieved to 45 μm) by centrifugation (Moore and Reynolds Reference Moore and Reynolds1989). Due to the low density of diatom frustules (<0.12–0.25 g cm–3), separating the siliceous skeletons out of the <2 μm fraction was difficult; therefore, contamination of the <2 μm fraction with frustules was assumed. Sodium hexametaphosphate was used during particle separation as a deflocculant. During the clay separation, each size fraction (>45 μm, >2 μm, and <2 μm) was weighed to determine the particle-size distribution, while the >45 μm and 45–2 μm fractions were analyzed as randomly oriented powders when required to clarify the quantitative results from the bulk analyses. The <2 μm fractions were prepared as oriented mounts via the pipette method and analyzed in their air-dried, ethylene glycol-solvated (overnight at 60°C), and heat-treated (400°C and 550°C) states. All XRD data were collected using a PANalytical X’Pert Pro diffractometer (Almelo, Netherlands) in Bragg-Brentano geometry operated at 40 kV and 45 mA using Fe-filtered CoKα radiation over the range 4–90°2θ with a 1/4° divergence slit, and 1/2° anti-scatter slit. Randomly oriented powders were prepared as back-loaded mounts and analyzed with 0.2 mm Soller slits with an effective counting time of 90 s/step; the oriented clay mounts were analyzed using 0.4 mm Soller slits with an effective counting time of 15 s/step, repeated three times to maximize intensity. For comparative analysis, 13 references clays were prepared and analyzed with the same procedures. The references were from the Source Clays Repository of The Clay Minerals Society, Iowa State University, and Pacific Northwest National Laboratory and included: K-saturated hectorite, SHCa-1 hectorite, SWa-1 ferruginous smectite, PFl-1 palygorskite, SepSp-1 sepiolite, nontronite (Washington), nontronite (Pennsylvania), saponite, bentonite (South Panther Creek), bentonite volclay, API no. 25 Upton montmorillonite, montmorillonite (Arkansas), and montmorillonite (Mississippi). All XRD data were analyzed using the Highscore Suite (Degen et al. Reference Degen, Sadki, Bron, Köng and Nénert2014). Semi-quantitative modal mineralogy was determined via Rietveld refinement of the diffraction data collected from randomly oriented powders.
Visible and Near-Infrared Spectrometry
Near-infrared reflectance spectra were measured on the <2 μm separates relative to a USP white reference standard using an ASD TerraSpec® 4 Hi-Res Mineral Analyzer (Boulder, Colorado, USA) over the range 350–2500 nm and RS 3 Spectral Acquisition Software. Contact reflectance spectra were collected on air-dried, ground powders at an emission angle of 0°. Reflectance data were analyzed using The Spectral Geologist (TSG®) developed by CSIRO Earth Science and Resource Engineering (Canberra, Australian Capital Territory, Australia).
Sediment and Water Characteristics
Sediment plugs of selected samples were mounted in epoxy and prepared as thin sections for petrographic analysis. Subsamples of each bulk fraction were treated to dissolve evaporite minerals and aid in the identification and weight of the insoluble minerals. First, gypsum was dissolved in 10% (NH4)2CO3(aq), then filtered with P4 filter paper to remove the resulting (NH4)2SO4(aq) (Eq. 1).
The filtrate was discarded and the remaining sediment and calcite sludge was then treated with 0.5 M HCl, filtered, rinsed with Milli-Q water, and dried at 60°C (Eq. 2).
Major and minor elements of surface-water samples and selected <2 μm fraction subsamples (digested in aqua regia) were analyzed for a 31 element suite using inductively coupled plasma optical emission spectroscopy, ICP-OES) (Agilent Pro-axial CCD spectrometer, Santa Clara, California, USA) by the Analytical Services Unit at Queen’s University, Kingston, Ontario. The analytical accuracy was verified using the standard reference material MESS-3 for sediment samples, and the EU-H-4 standard and Cranberry-05 CRM (for Si) for the water samples. Surface water pH was measured in situ during sampling while sediment pH was measured in duplicate via the saturated paste method described by Kalro (Reference Kalro1995) using a Vernier Tris-Compatable Flat pH Sensor and LabQuest Interface (Vernier Software and Technology, Beaverton, Oregon, USA). The morphology and qualitative chemical composition of selected <2 μm fractions were analyzed using an FEI Quanta 650 field-emission scanning electron microscope (FEG-SEM, Thermo Fisher Scientific, Hillsboro, Oregon, United States). Samples were mounted on carbon tape, carbon coated, and observed at 15 or 25 kV.
Results
Three types of sediment were identified from the basins based on location, field notes, texture, and sedimentary environment: (1) regolith sediment, (2) detritus-rich mudflat sediment, and (3) mudflat sediments composed primarily of authigenic minerals. The mineralogy, texture, and abundance of surface crusts and efflorescence varied between each salar and were analyzed via XRD. The weights from the particle-size separation (data not shown) and XRD data were used to determine the mineralogy of the sediments.
X-ray Diffraction
For all the bulk sediment samples, the XRD patterns contained a broad background centered on 25°2θ (4.13 Å) representative of 'X-ray amorphous' phases such as biogenic opal (diatom frustules), volcanic glass, allophane, or imogolite (Elzea et al. Reference Elzea, Odom and Miles1994; Bishop et al. Reference Bishop, Rampe, Bish, Abidin, Baker, Matsue and Henmi2013). Major (>10%) and minor (<10%) minerals present in each sediment type for each salar are presented in Table 2. The '1-Regolith' sediments were sampled from dry depressions and margins surrounding the salars and are composed of 'X-ray amorphous' material, clay minerals, and igneous minerals consistent with the andesitic to rhyolitic surroundings (quartz, cristobalite, biotite, hornblende, albite-anorthite, clinopyroxene, and orthopyroxene) (Table 2). Sediments ranged from fine-grained (<45 μm) silt and clay in dry depressions to coarse (0.5–2 mm), poorly sorted, sub-angular grains surrounding the salar margins. The mineral suite comprising the '2-Detrital' sediments is similar to the '1-Regolith' sediment with the addition of evaporites including calcite, gypsum, halite, and sylvite (Table 2). The sediments are well sorted, sub-angular grains deposited in alluvial to deltaic environments. The '3-Authigenic' mudflat sediments are composed of 'X-ray amorphous' material, clay minerals, and a complex suite of evaporite minerals including calcite, celestine, dolomite, eugsterite, gypsum, mirabilite, authigenic quartz, authigenic cristobalite, sylvite, and ulexite (Table 2). They consist of an extremely fine-grained matrix with evaporite minerals ranging from <2 μm to 1 cm crystals of mirabilite. Surface crusts and efflorescence across the salars vary considerably; however, they are predominantly composed of Na-sulfate minerals (thenardite/mirabilite) with small quantities of analcite, aphthitalite, calcite, cristobalite, gypsum, eugsterite, halite, howlite, sylvite, and syngenite (Table 2).
* Moisture content and field descriptions identified that mirabilite is present with thenardite; italicized minerals are major constituents in specific layers. Numbers correspond to sediment type outlined in text.
† X-ray amorphous material dominated by siliceous diatom frustules identified in SEM images.
‡ For clay mineral descriptions, see Table 3.
§ Petrographic examinations of thin sections from sediment plugs for selected samples confirm the presence of biotite; feldspars range from anorthite to albite; pyroxenes include clinopyroxene and orthopyroxene.
X-ray Diffraction of Clay Minerals
The XRD patterns for nearly all of the oriented, <2 μm fractions displayed evidence of swelling clays and X-ray amorphous material except for sample 45, which contained X-ray amorphous phases only (Table 3).
1-Regolith.
XRD patterns for the glycolated regolith clay fractions all contain a low-angle reflection (15.1 to 17.3 Å) indicative of smectite or interstratified illite-smectite (I-S), discrete illite (10 Å and 5.0 Å), and a ~7.2 Å peak, probably from chlorite with weak (001) and (003) reflections given the arid environment (Fig. 2a). Interstratified clays were identified based on irrational basal spacings and a glycolated d 001 spacing <16.9 Å that collapsed to 10 Å with heat treatment (400°C) (Table 3) (Moore and Reynolds Reference Moore and Reynolds1989). When the (002/003) I-S peak could be resolved, most samples were estimated to have between 20 and 50% randomly interstratified illite layers. These regolith sediments represent basin filling Quaternary alluvium.
2-Detrital.
The majority of the XRD patterns from glycolated samples of the detritus-rich mudflat sediments were similar to the regolith clay patterns which showed detritus-rich mudflat sediment containing low-angle reflections with glycolated d 001 spacings ranging from 14.4 to 17.2 Å mixed with discrete illite and chlorite (Table 3) (Fig. 2b). The percentage of illite interstratification, estimated from the (002/003) I-S reflections, ranged from 30 to 80%. This increase in illite interstratification indicated a general increase in layer charge compared to the regolith clays common with diagenetic alteration in saline or alkaline solutions (Deocampo et al. Reference Deocampo, Behrensmeyer and Potts2010; Hugettt et al. 2016). Discrete illite was not present in the XRD patterns from samples 31, 37, and 38 (sampled from below the sediment–water interface); instead they contained a broad, asymmetric shoulder at 10 Å resembling kerolite (Table 3) (Brindley et al. Reference Brindley, Bish and Wan1977; Léveillé et al. Reference Léveillé, Longstaffe and Fyfe2002).
3-Authigenic.
Authigenic mudflat sediments produced oriented XRD patterns significantly different from the regolith and detrital-rich mudflat sediments. All of the XRD patterns from glycolated samples displayed a very broad, low-angle reflection that swelled with ethylene glycol and either collapsed with heat (400°C) treatment or was destroyed (Fig. 2c). The d spacings for the glycolated, low-angle reflections ranged from 14.8 to 17.6 Å (Table 3). These clays were identified as poorly formed, low-swelling smectite indicative of high layer charge materials with some illite-like layers based on the broad reflections and lack of a 10 Å peak with heat treatment (Deocampo et al. Reference Deocampo, Behrensmeyer and Potts2010). The only other features present (aside from micron-sized calcite) were weak peaks at 10 Å and 4.5 Å that did not swell with ethylene glycol and resembled kerolite (Brindley et al. Reference Brindley, Bish and Wan1977; Léveillé et al. Reference Léveillé, Longstaffe and Fyfe2002). Aside from the reduced d-spacings of the (001) peak for some samples, higher-angle reflections indicative of I-S or kerolite-stevensite interstratification were not observed in the XRD patterns. The samples extracted from the sediment–water interface under the surface water (Nos. 44, 45, and 79) exhibited XRD patterns for glycolated samples dominated by broad reflections consistent with biogenic silica.
* Peak parameters for the (001) reflection were measured on glycolated smectite peaks ~6°2θ. Sme = smectite; Sme* = low-swelling smectites; Ilt = illite; I-S=interstratified illite-smectite; Chl = chlorite; Ker = kerolite; and Mnt = montmorillonite.
Due to interferences in the XRD patterns from diatom frustules and the broad 10 Å shoulder from kerolite in many samples, the coherent scattering domain for the glycolated (001) smectite reflection could not be calculated accurately. Instead, the integral breadth (β i ) of the (001) smectite reflection was used as a proxy to compare the nature of the smectite clays between sediment classes and the Source Clays and clay references analyzed (Table 3; Fig. 3). Broadening of the (001) reflection, resulting in a higher β i value, occurs when materials have limited coherent domains from small crystallite size and randomly interstratified layers (Martin de Vidales et al. Reference Martin De Vidales, Pozo, Alia, Garcia-Navarro and Rull1991). Analysis of box and whisker plots of the β i from (001) smectite reflections by sediment class revealed a systematic increase in both the overall broadness and range of the β i measured from the reference Source Clays (used as a benchmark for well crystallized clay) to the Bolivian clays in the 3-Authigenic sediments (Fig. 3). The Source Clays had the smallest values for β i , ranging from 0.6 to 1.5 with the lowest average ( = 1.0), while the authigenic clays displayed the greatest range in values for β i , 2.4–10.5, and the highest average ( = 5.0). The β i values of the regolith and detrital-rich mudflat clays were fairly similar, ranging from 2.1 to 3.9 ( = 3.0) for the regolith samples and from 2.4 to 3.9 ( = 3.3) for the detrital samples.
Visible and Near-Infrared Spectroscopy
The major absorption bands identified in the samples were located around 1410 nm (7092 cm–1) from OH and H2O overtones (OH2v and H2O2v), 1910 nm (5236 cm–1) from H2O combinations (H2Ov+δ), 2210 nm (4525 cm–1) from Al–OH combinations, and/or 2310 nm (4329 cm–1) from Mg–OH combinations (Bishop et al. Reference Bishop, Pieters and Edwards1994, Reference Bishop, Rampe, Bish, Abidin, Baker, Matsue and Henmi2013). A flat, broad reflectance band between 2210 and 2310 nm was also observed in most of the mudflat samples characteristic of opal or hydrated silica (Fig. 4) (Milliken et al. Reference Milliken, Swayze, Arvidson, Bishop, Clark, Ehlmann, Green, Grotzinger, Morris, Murchie, Mustard and Weitz2008). This Si–OH feature was excluded from Figs. 5, 6, and 7 for clarity unless it was the most prominent feature in the reflectance spectra. Subtle shifts in the 1410 nm band (OH2v), the position of which depends on the corresponding cation (Al, Mg, or Fe), could not be distinguished due to the presence of hydrated silica in all of the samples, which contributed a very broad absorption over the region (Bishop et al. Reference Bishop, Rampe, Bish, Abidin, Baker, Matsue and Henmi2013).
1-Regolith.
Reflectance spectra of the regolith <2 μm fractions all contained diagnostic absorption bands for dioctahedral smectite located at 1414 nm (7072 cm–1), 1908 nm (5241 cm–1), and 2208 nm (4529 cm–1) (Fig. 5). The absence of any Mg–OH absorptions suggested that most of the Mg2+ released from the weathering of pyroxene and hornblende in the regolith was not incorporated into the pedogenic smectite but, instead, contributed to Mg enrichment of the surface water and sediments in the salars. These Al-rich clays are available for transport into the mudflats of the salars.
2-Detrital.
Reflectance spectra from the <2 μm fractions from detrital mudflat sediments were characterized by metal–OH absorption bands for both dioctahedral clay (Al–OH) at 2208 nm (4529 cm–1) and trioctahedral clay (Mg–OH) at 2313 nm (4323 cm–1) (Fig. 6). Spectral analysis revealed that these NIR data resembled montmorillonite-beidellite mixed with magnesite; bulk XRD patterns provided no evidence of magnesite, dolomite, or Mg-calcite, however. Samples 31, 37, and 38 were collected from beneath standing water and exhibited the broad absorption bands from hydrated silica, 2200–2300 nm and, in the case of 37 and 38, Mg–OH absorptions were absent. The clay minerals in the detrital-rich mudflat are a mix of detrital dioctahedral smectite and authigenic trioctahedral smectite. Substitution of Mg into the octahedral sheet (AlMg–OH) should shift the Al–OH band to higher wavelengths near 2212 nm (4520 cm–1) but this was not observed, suggesting that the Al–OH and Mg–OH phases are distinct (Bishop et al. Reference Bishop, Madejová, Komadel and Fröschl2002a).
3-Authigenic.
The NIR spectra of authigenic mudflat sediment displayed absorption features consistent with trioctahedral clay with the main absorption band located at 2312 nm (4325 cm–1) from Mg–OH (Fig. 7). A broad feature was identified between 2200 and 2300 nm (4545–4348 cm–1) in many of the samples and resembles hydrated silica. For samples 58, 60, and 74 the H2O combination band at 1910 nm (5236 cm–1) was shifted to longer wavelengths (1938–1942 nm) due to the presence of vegetation in the samples noted as green and pink algae in the field (Clark Reference Clark and Rencz1999). No evidence for Al–OH absorptions was observed in the NIR spectra, indicating that the dioctahedral clay from the regolith and detrital mudflat sediment is absent. Samples that were dominated by hydrated silica (42, 44, 45, 79) were fresh sediments sampled from under the surface water similar to the detrital sediment samples 31, 37, and 38.
Scanning Electron Microscopy
1-Regolith.
Clay minerals separated from the regolith sediments displayed a mix of well defined, hexagonal, spherical, and platy minerals with small, flaky particles consistent with XRD results from dioctahedral smectite (montmorillonite or beidellite), illite, and chlorite (Fig. 8a). The EDS spectra of the flaky particles with classic smectite textures confirmed that the clays are Al-rich (data not shown).
2-Detrital.
The detrital-rich sediments displayed platy and spherical minerals similar to those present in the regolith, as well as poorly-preserved remnants of diatom fragments and flaky masses (Fig. 8b). Within the flaky masses, individual particles were still distinguished by their well defined edges.
3-Authigenic.
The authigenic mudflat sediments are composed predominantly of diatom frustules with small platy particles mixed with fluffy aggregates (Fig. 8c–f). Within the authigenic sediments the diatom frustules display evidence of dissolution and coating/replacement with fine-grained material as the delicate morphology of the microfossils is obscured for many of the samples (Fig. 8c,d). The NIR reflectance spectra and EDS analysis (data not shown) confirmed that the fine-grained aggregates are rich in Mg and Si. The <2 μm wide scalenohedral crystals were of calcite, confirmed in the XRD spectra and EDS analysis of the grains (data not shown).
Samples taken from below the sediment–water interface (Nos. 42, 44, 45, 79) are composed of well preserved diatom frustules still displaying the delicate morphology of the stria and raphe on the porous structure (Fig. 8e,f). Sparse aggregates of fluffy masses are still present in these samples, representing poorly formed minerals. The dominance of diatom frustules within these sediments is consistent with the NIR spectra and XRD profiles, which exhibited a strong absorption for hydrated silica and a broad hump caused by amorphous material at 25°2θ. These samples represent the most recent sedimentary deposits from the salar, confirming that endogenic calcite and diatom frustules are the main materials settling out of the surface water, and that neoformed phyllosilicate minerals are not precipitating directly from solution in appreciable amounts. The fine-grained, fluffy material in the rest of the samples is a product of frustule dissolution and alteration rather than neoformed phyllosilicates.
Chemical Analysis
The Al:Fe:Mg ternary diagrams for selected <2 μm fractions (Fig. 9a) revealed a pattern of increasing Mg content from the regolith clays to the mudflat clays where: 1-Regolith < 2-Detrital < 3-Authigenic. These data are consistent with the NIR spectra for each sediment class where the regolith clays are dominated by Al-dioctahedral smectite and the authigenic clays are dominated by Mg-trioctahedral smectite. The Fe:K:Li ternary diagrams (Fig. 9b) showed that both the 1-Regolith and 2-Detrital clays contain greater amounts of Fe (>60%) than of K and Li; some of the detrital clays displayed a slight increase in K, however, consistent with the increase in illite layers, or the decrease in the glycolated d 001, observed in the XRD patterns (Fig. 2b). The Mg-rich clays sampled from the 3-Authigenic sediments, conversely, displayed increasing concentrations of K and Li relative to Fe. Given the small, glycolated d 001 spacings for many of the 3-Authigenic clays, K must be present as an interlayer cation contributing to illite-like layers deriving from the low-temperature illitization of smectite. The large Li concentrations in the 3-Authigenic clays could be from hectorite or Li-rich saponite/stevensite in the evaporative basin (Fig. 9b). Minor/trace elements and heavy metals were either below detection limits or enriched in specific salars with no apparent pattern to the sediment class and were not included, therefore.
Water Chemistry
Surface-water chemistry data from studies published over the period 1983 to 2014 (and including the current study) were compiled to understand the range of geochemical conditions experienced within each salar (Badaut and Risacher Reference Badaut and Risacher1983; Servant-Vildary and Roux Reference Servant-Vildary and Roux1990; Risacher and Fritz Reference Risacher and Fritz1991; Dejoux Reference Dejoux1993; Sylvestre et al. Reference Sylvestre, Servant-Vildary and Roux2001) (Fig. 10). Even though the data span ~30 y and are highly variable, the general trend in the Mg:Si ratio of the surface water compared to the inflow water is consistent across studies. The variability in the data is expected for closed basins where ionic concentrations fluctuate with evaporation and depend on the time and season of sampling. Compared to inflow waters that are rich in dissolved Si (Mg:Si < 1), the surface waters of the salars are significantly enriched in Mg2+ with respect to dissolved Si (Mg:Si > 1) (Fig. 10). This Mg2+ enrichment is especially pronounced in HED, CHKT, and CHU where the mean Mg:Si ratios were 130, 98, and 73, respectively. Surface-water pH was also highly variable for all the salars over the 30 y period, fluctuating between neutral and alkaline for most salars, although the pH would depend on the season when sampling was carried out (i.e. wet vs dry) (Table 1; Fig. 11). For HED, RAM, and CHKT the pH values measured/reported were consistently < 8.5, which is the threshold consistently identified for Mg-silicate precipitation from solution (Badaut and Risacher Reference Badaut and Risacher1983; Tosca et al. Reference Tosca, Macdonald, Strauss, Johnston and Knoll2011; Tosca Reference Tosca, Pozo and Galan2015).
Discussion
1-Regolith.
Clay minerals identified in the <2 μm fraction of regolith samples are consistent with illite, chlorite, and randomly interstratified I-S rich in Al-dioctahedral smectite. These samples were all characterized by Al–OH absorptions near 2208 nm (4528 cm–1) in the NIR spectra, chemical composition dominated by Al, and low-angle reflections in the XRD data that swell with ethylene glycol treatment. This suite of dioctahedral phyllosilicates is expected from pedogenic weathering in the soils (regolith) surrounding the salars. Weathering of biotite, feldspar, pyroxene, and hornblende present in the surrounding andesites and rhyodacites has been shown to produce Al-rich phases while releasing Fe and Mg to solution (Hover and Ashley Reference Hover and Ashley2003; Deocampo et al. Reference Deocampo, Cuadros, Wing-Dudek, Olives and Amouric2009). Clay minerals formed via weathering of detrital minerals also tend to have larger crystallite sizes and produce sharper reflections in XRD, both of which were observed for the regolith clays. These dioctahedral phyllosilicates represent potential detrital sources for clay minerals entering the salars.
2-Detrital.
Clay minerals in the detritus-rich mudflats are a mix of detrital (illite, chlorite, and I-S-rich in Al-dioctahedral smectite) and authigenic (Mg-smectite) phyllosilicates. All detritus-rich mudflat sediments displayed absorptions from both Al–OH (2208 nm; 4528 cm–1) and Mg–OH (2313 nm; 4323 cm–1) combinations in the NIR spectra with an overall shift toward Mg-rich chemistries compared to the 1-Regolith clays. The XRD patterns also exhibited evidence of broadening of the glycolated (001) reflections and increased layer charge (smaller glycolated d spacings) of the smectite minerals compared to the regolith smectite clays. Without any shifts in the Al–OH band to higher wavelengths near 2212 nm (4520 cm–1) signifying Mg–OH contributions, these observations suggest that authigenic Mg-smectite is forming on, or replacing, dioctahedral smectite within the mudflats. The authigenic clay minerals are characterized by small crystallite sizes, structural disorder, and high layer charge consistent with low-temperature illitization due to increasing octahedral Mg in the saline environment evident in the broad, asymmetric reflections and Mg content of the minerals (Figs. 2, 3, and 9). Similar broadening of the (001) reflection was measured by Cuevas et al. (Reference Cuevas, Vigil de la Villa, Ramirez, Petit, Meunier and Leguey2003) who used peak deconvolution analysis to conclude that combining a broad 18 Å reflection from authigenic stevensite with a sharper 17 Å reflection from detrital dioctahedral smectite would produce the broad reflections observed. The transformation of detrital Al-dioctahedral smectite to authigenic Mg-trioctahedral smectite has been demonstrated in Abert Lake, Oregon (Jones and Weir Reference Jones and Weir1983; Banfield et al. Reference Banfield, Jones and Veblen1991); in Miocene lacustrine units in the Madrid basin, Spain (Bellanca et al. Reference Bellanca, Calvo, Censi, Neri and Pozo1992); in the Salina Ometepec, Baja, California (Hover et al. Reference Hover, Walter, Peacor and Martini1999); in Neogene Vicalvaro deposits, Madrid, Spain (Cuevas et al. Reference Cuevas, Vigil de la Villa, Ramirez, Petit, Meunier and Leguey2003); in Lake Eyasi and Ngorongoro Crater, Tanzania (Hover and Ashley Reference Hover and Ashley2003; Deocampo Reference Deocampo2005); and in the Pliocene Olduvai Gorge, Tanzania (Deocampo et al. Reference Deocampo, Cuadros, Wing-Dudek, Olives and Amouric2009). The process of the transformation has been theorized by several studies and includes: topotactic growth of Mg-rich layers onto existing detrital clays producing di- and tri-octahedral crystals, complete dissolution of Al-dioctahedral smectite and re-precipitation of Mg-trioctahedral smectite, and solid state transformation of dioctahedral clays to trioctahedral clays producing coexisting domains within the same sheet or an intermediate composition of the octahedral sheet (Banfield et al. Reference Banfield, Jones and Veblen1991; Cuevas et al. Reference Cuevas, Vigil de la Villa, Ramirez, Petit, Meunier and Leguey2003; Hover and Ashley Reference Hover and Ashley2003; Deocampo et al. Reference Deocampo, Cuadros, Wing-Dudek, Olives and Amouric2009). Due to the lack of evidence for octahedral substitution of Mg for Al in the NIR reflectance spectra (i.e. distinct bands for Al–OH and Mg–OH), the most likely process for authigenic clay formation in Bolivia appears to be topotactic growth of the Mg-rich phases onto existing Al-rich phases.
3-Authigenic.
Clay minerals in the authigenic sediments are consistent with Mg-smectite (probably stevensite or saponite), and kerolite forming on/replacing diatom frustules. The enrichment of Li in the authigenic clays in CHKT could indicate that hectorite is forming; however, analytical electron microscopic analysis coupled with transmission electron microscopic analysis (AEM/TEM) of individual clay particles would be required to establish if enough Li is present in any individual particles to identify them as hectorite. If hectorite were present, however, a shift in the trioctahedral combination band to 2306 nm (4335 cm–1) should be observed but was not (Bishop et al. Reference Bishop, Murad and Dyar2002b). The NIR spectra for these sediments are dominated by hydrated silica and/or Mg–OH absorptions (2312 nm; 4325 cm–1). X-ray diffraction analyses produced very broad basal reflections the positions of which changed with ethylene glycol treatment and weak reflections from kerolite (10 and 4.5 Å). Broad basal reflections are indicators of authigenic formation for clay minerals (Deocampo et al. Reference Deocampo, Cuadros, Wing-Dudek, Olives and Amouric2009, Reference Deocampo, Behrensmeyer and Potts2010). Within the mudflats, kerolite is precipitating from the Mg2+ -enriched solutions, characterized by an excessively high Mg2+/SiO2 activity ratio, which leach through the diatomaceous sediment (Galán and Pozo Reference Galán, Pozo, Galán and Singer2011; Tosca and Masterson Reference Tosca and Masterson2014). Subsequent removal of Mg2+ via kerolite precipitation shifts the solutions into the stability field of stevensite, allowing the smectite to precipitate (Galán and Pozo Reference Galán, Pozo, Galán and Singer2011; Tosca and Masterson Reference Tosca and Masterson2014).
The samples that exhibited glycolated d 001 spacings of <17 Å are consistent with a high octahedral layer charge coupled with moderate K+ fixation (Deocampo et al. Reference Deocampo, Cuadros, Wing-Dudek, Olives and Amouric2009). A similar precursor smectite was described in the Tagus basin (Spain) that had variable interlayer spacings (10–18 Å) and a high layer charge from deficiencies in the octahedral cations (Santiago de Buey et al. Reference Santiago de Buey, Barrios, Romero and Montoya2000). The absence of I-S interstratification in the XRD data could be due to incomplete K+ fixation in the Na-rich solutions. Deocampo et al. (Reference Deocampo, Cuadros, Wing-Dudek, Olives and Amouric2009) suggested layer-charge increases and K+ fixation are decoupled, and that K+ fixation occurs under freshwater conditions when the K/Na ratio is higher. The rest of the Mg-rich samples that swelled to ~17 Å resemble authigenic kerolite-stevensite similar to the clays from palustrine environments from the Neogene Madrid basin (Spain) (Pozo and Casas Reference Pozo and Casas1999). The Neogene kerolite is thought to have originated via precipitation from solution or a gel-like medium (neoformation) based on the colloform texture, chemical purity, and XRD reflections; while early diagenetic processes (increasing salinity via evaporation and drying) converted the kerolite to kerolite-stevensite prior to pure stevensite formation (Pozo and Casas Reference Pozo and Casas1999). Although diatoms were not included in the observations by Pozo and Casas (Reference Pozo and Casas1999) diatomaceous marls are present in the Miocene sediments of the Madrid basin, indicating the green algae could be affecting the mineral assemblages present (Bellanca et al. Reference Bellanca, Calvo, Censi, Neri and Pozo1992).
The replacement of biogenic silica with Mg-smectite is consistent with the observations of Badaut and Risacher (Reference Badaut and Risacher1983) who reported stevensite replacing diatom frustules in water samples when the pH was above >8.5 in the same salars as the present study; the pervasiveness of the diatom frustules throughout the mudflats (~1 m of diatomaceous sediment) has not been reported previously, however. Authigenic smectite formation within the mudflats, however, was not restricted by the pH of the surface water as found by Badaut and Risacher, (Reference Badaut and Risacher1983) because all of the sediments exhibited evidence of phyllosilicate formation due to the concentration of solutes within the pore fluids. The mudflats adjacent to saline lakes are dynamic environments subject to periodic wet/dry cycles and strong capillary forces that concentrate solutes, especially Mg2+ after calcite and gypsum precipitation.
Widespread evidence of sepiolite was absent in the XRD patterns and NIR reflectance spectra despite the silicate-rich, brackish environments expected to favor the neoformation of sepiolite from evaporating waters. Direct precipitation of sepiolite from solution has been assumed previously in the dolomitic sepiolite beds of the Amargosa Flat playa, Nye County, Nevada (Papke Reference Papke1972), in the Allou Kagne deposit, Senegal (García-Romero et al. Reference García-Romero, Suárez, Santaré and Alvarez2007), and in the Serinhisar-Acipayam basin, Turkey (Akbulut and Kadir Reference Akbulut and Kadir2003). Based on the Mg:Si ratios measured during the past 30 y, waters within the salars clearly become significantly enriched in Mg, but not silica, during evaporative concentration. From SEM analysis of the authigenic mudflat sediments not obscured by detritus, evidence is clear that these sediments are dominated by diatom frustules. Diatoms are capable of surviving under wide ranges of salinities and pH, and effectively remove silica from the water column via biologically mediated uptake to form their siliceous skeletons (DeMaster Reference DeMaster2003). The ultimate fate of the silica depends on silica solubility and the geochemical conditions present following the death of the biota: i.e. preservation vs dissolution. Factors that affect silica solubility and diatom dissolution have been summarized by Hecky and Kilham (Reference Hecky and Kilham1973); DeMaster (Reference DeMaster2003), Ryves et al. (Reference Ryves, Battarbee, Juggins, Fritz and erson2006), Flower and Ryves (Reference Flower and Ryves2009), Loucaides et al. (Reference Loucaides, Van Cappellen, Roubeix, Moriceau and Ragueneau2012), and Tutolo and Tosca (Reference Tutulo and Tosca2018) and include:
• pH: dissolution increases with increasing pH, due to deprotonation of silanol groups that helps promote hydrolysis of the siloxane bonds, especially above pH >9 where dissolution increases exponentially; the pKa for silicic acid at 25°C is 9.8.
• Temperature: The rate of frustule dissolution increases linearly with increasing temperature from 4 to 30°C.
• Salinity: alkali salts increase frustule dissolution until solutions become hypersaline; further increases in salinity decrease silica solubility, especially when divalent cations are present (Ca2+, Mg2+, and Fe2+).
• Depth of the water column: a deeper water column allows more time for the frustules to dissolve before sedimentation. In deep lakes and ocean environments, many diatoms dissolve within the water column before settling.
• Sinking velocities: the sinking velocity is affected by the size (dependent on the nutrient status of the water) and aggregation of diatom cells where larger cells/groups sink faster, allowing less time for dissolution.
• Al incorporation: soluble Al present during diatom growth can be incorporated into the siliceous skeleton and ultimately decreases the solubility of the frustules.
In the Bolivian salars, temporal variations in the neutral-alkaline solutions and brackish-hypersaline conditions create a range of biogeochemical conditions within each salar and between salars. For most of the salars, however, the slightly alkaline waters (pH <9), cold temperatures (0–10°C), and shallow water depths (<30 cm) that promote rapid rates of sedimentation act to slow the dissolution of biogenic silica. These conditions promote the preservation of diatom frustules, especially in HED, RAM, and CHKT where silica is removed from solution and cannot contribute to authigenic sepiolite formation. Under more optimal conditions for dissolution, Michalopoulos et al. (Reference Michalopouos, Aller and Reeder2000) measured the complete conversion of diatom frustules to aluminosilicate clays in 23 months. In the Bolivian salars, conversely, the sedimentation of the frustules has built-up diatomaceous muds >1 m thick in some of the salars. Instead of dissolving in the water column, the deposited frustules become a slow-release source of Si as Mg-rich fluids leach through sediment pores allowing kerolite and Mg-smectite to form during evaporative concentration in the mudflats. Samples from the current study included transects across the mudflats of the salars; however, samples were not taken where inflowing waters mix with surface waters. Deposits of sepiolite might be present where Si-rich inflow waters react with saline solutions in the salars before silica uptake by the diatoms occurs.
Instead of neoformation (precipitation from solution), in many other basins sepiolite seems to form from the dissolution of trioctahedral Mg-smectite, within soils and lake margins, and from the dissolution of dolomite (Chahi et al. Reference Chahi, Fritz, Duplay, Weber and Lucas1997; Galán and Pozo Reference Galán, Pozo, Galán and Singer2011). In sediments from the Jbel Rhassoul, Morocco, sepiolite was found replacing neoformed stevensite during climatic periods characterized by low salinity and high concentrations of dissolved silica from fresh water entering the basin (Chahi et al. Reference Chahi, Fritz, Duplay, Weber and Lucas1997). In the Ahiler area of the Eskişehir sepiolites, Turkey, late-stage diagenetic events caused the dissolution of dolomite, and diatomites favored the formation of sepiolite (Ece and Çoban Reference Ece and Çoban1994). In the Madrid basin (Spain), weathering of granites and the dissolution of siliciclastic grains during calcretization in alluvial fans produced Si- and Mg-enriched groundwater that led to the precipitation of sepiolite followed by opal within the detrital lake sediments (Bustillo and Bustillo Reference Bustillo and Bustillo2000). Additionally, in the Vicalvaro sepiolite deposit, Madrid Neogene Basin, Spain, sepiolite initially interpreted to have formed via neoformation by precipitation is now interpreted as forming from the diagenetic alteration of Mg-rich smectite (Cuevas, et al. Reference Cuevas, Vigil de la Villa, Ramirez, Petit, Meunier and Leguey2003; Pozo and Calvo Reference Pozo and Calvo2018). In the Sariyer Formation (Turkey), the Intermediate Unit in the Madrid basin (Spain), and the Tuzgölü basin (Turkey), sepiolite-palygorskite were found replacing dolomite during dry periods coupled with sulfide oxidation in the presence of disordered silica (Kadir et al. Reference Kadir, Eren and Atabey2010; Gürel and Özcan Reference Gürel and Özcan2016). In all of these cases, sepiolite-palygorskite apparently may be restricted to subsurface environments and post-depositional alteration of Mg-smectite and dolomite. In subsurface environments (pore fluids of detrital mudflats, soils, and deltas) the sunlight required for algal growth is absent, which inhibits the growth of diatoms and allows dissolved Si concentrations to increase and sepiolite-palygorskite to precipitate. For the Bolivian salars, post-depositional alteration could convert the smectites in the mudflats to sepiolite if freshwater conditions returned or dissolved silica increased significantly.
Conclusions
Clay minerals forming on the Bolivian Altiplano follow a pattern of increasing Mg-enrichment from detrital, Al-dioctahedral phyllosilicates to authigenic, Mg-trioctahedral phyllosilicates transforming detrital clays and diatom frustules. Despite the range of geochemical environments across the salars, from brackish to hypersaline and neutral to alkaline, all of the authigenic clays forming are kerolite or Mg-smectite, probably stevensite. The pH of evaporating solutions defines an important threshold for the neoformation of clays from solution at pH 8.5 (both sepiolite and Mg-smectite); however, the Mg:Si ratio of the pore solutions has a greater impact on the mineralogy of the clays that are forming. Very high Mg:Si ratios favor kerolite, especially in detritus-poor sediments, while decreasing Mg:Si ratios favor Mg-smectite. Authigenic clays in the mudflats are all characterized by high layer charge, small crystallite sizes, and/or stacking disorders producing very broad, low-angle X-ray reflections. Aside from evaporites and phyllosilicates, the mudflats are dominated by diatom frustules, explaining the depletion of dissolved Si in the evaporating waters compared to the inflowing waters. Diatomaceous sediments >1 m thick have accumulated due to slow rates of diatom dissolution as a result of cold temperatures, shallow water depths, salinity, and neutral–slightly alkaline pH solutions. This removal of Si from solution by diatoms prevents sepiolite formation, which requires high concentrations of dissolved Si. In modern sediments, the activity of diatom frustules may even restrict sepiolite-palygorskite to subsurface environments where sunlight is absent and diatom growth is inhibited, thereby allowing dissolved Si to accumulate in the pore fluids. Burial of the Bolivian mudflat sediments could favor the conversion of smectite to sepiolite as the diatom frustules will act as a source of silica during potential, future, epigenetic and diagenetic reactions. Overall, the control of Si concentrations in solution by diatoms has a large impact on the types of clay minerals that will form, and may have been underestimated when analyzing ancient sediments.
ACKNOWLEDGMENTS
Many thanks to the Analytical Services Unit at Queen’s University for the chemical analysis by ICP-OES. The authors are very grateful to James Amonette from Pacific Northwest National Laboratory and Michael Thompson from Iowa State University for the clay mineral references. They also thank Jenny Huggett, Benjamin Tutolo, an anonymous reviewer, and the Editor-In-Chief, Joseph Stucki, for helpful and constructive comments. Field guides, Nelson and Marie from Betto Tours in Uyuni, Bolivia, are thanked for helping to navigate the southern Altiplano. This research was supported by the Mineralogical Association of Canada Travel Grant, the Marie Mottashed Graduate Scholarship, and an NSERC Discovery Grant to R. Peterson.