Blooming of a microbial community in an Ediacaran extreme volcanic lake system

Apr 27, 2023

Scientific Reports volume 13, Article number: 9080 (2023) Cite this article

Metrics details

Ancient aquatic sediments are critical archives for studying early microbial life and the types of environments in which they thrived. The recently characterized Amane Tazgart microbialites in the Anti-Atlas, Morocco, are a rare and well-preserved non-marine deposit that evolved in an alkaline volcanic lake setting during the Ediacaran Period. A multiproxy geochemical toolbox reveals evidence pointing to spatio-temporal ecosystem organization and succession related to changing lake water chemistry. This is marked by secular transition from a cold/dry climate, hypersaline alkaline thermophilic and anoxic–oxic community, to a stable state warm/wet climate fully oxygenated fresh to brackish water ecosystem, predominated by oxygenic stromatolites. Extreme dissolved Arsenic concentrations suggest that these polyextremophiles required robust detoxification mechanisms to circumvent arsenic toxicity and phosphate deficiency. We propose that self-sustaining and versatile anoxic to oxic microbial ecosystems thrived in aquatic continental settings during the Ediacaran Period, when complex life co-evolved with a rise in atmospheric oxygen content.

Modern aquatic microbial ecosystems are known to produce distinct spatio-temporal communities along geochemical gradients1,2,3 that tend to precipitate diagnostic biogenic minerals such as carbonates, iron oxides, and pyrite3,4,5,6,7,8,9,10,11,12. Geochemical characterization of these microbial habitats ties shale-normalized rare earth element (REE), yttrium (Y), and redox-sensitive trace element distribution to prevailing environmental redox conditions, chemistry, and mineral deposition11,13. Particularly, REE patterns normalized to Post Archean Australian Shales (PAAS), marked by light REE depletion, negative cerium (Ce) and slightly positive lanthanum (La) anomalies, and high Y to holmium (Y/Ho) ratios, collectively reflect oxygenated modern seawater mass conditions14; and in certain instances, non-marine aquatic ecosystems influenced by elevated alkalinity15. In fact, the firm connection of REE + Y (REY) fractionation to dissolved carbonate ion concentration, serves as a valuable tool for reconstructing alkalinity trends, especially in restricted aquatic ecosystems15.

On the other hand, the relative abundance of europium (Eu) is often associated with hydrothermal fluid enrichment, feldspars, and diagenesis16,17. While negative Eu anomalies may represent the original composition of the aquatic environment from which limestones formed16, and La, gadolinium (Gd), and Y anomalies, mineral complexation phenomena unique to the hydrosphere18. Hence, authentic, uncontaminated marine, estuarine, and lacustrine water chemistries are often deduced from REY systematics.

Various geochemical palaeosalinity proxies utilizing trace elemental ratios in bulk mudstones are among the most popular. They are founded on the observation that certain non-redox-sensitive elements show relatively higher concentrations in seawater such as B, Y, and Sr, or freshwater, like Ga, Ho, and Ba, on a salinity-normalized basis19,20,21,22,23,24. Elemental cross plots allow estimation of expected salinity baselines for each proxy: (1) B/Ga > 6 is typical of marine facies, 3–6 brackishness, and < 3 freshwater20,23; (2) Sr/Ba is > 0.5 characterizes marine deposits, 0.2–0.5 brackishness, and < 0.2 freshwater20; and (3) B/K > 40 μg/g distinguishes marine from brackish/freshwater sediments22,25. Furthermore, several geochemical parameters such as Ni/Co, Cu/Zn, V/Cr, V/(V + Ni), U, U/Th, authigenic trace elements, and Ce/Ce*, are used to provide estimates for paleo-redox depositional conditions26,27,28. However29, in the absence of living biomass and the associated biosignatures that enable comparative characterization of modern aquatic microbial communities, very little evidence exists directly linking microbial community succession patterns to long-term environmental drivers in the Proterozoic biosphere. This has led to considerable difficulties in our ability to unambiguously characterize and tie microbial processes to shifting geochemical gradients, environmental disturbances, and geodynamics in deep-time ecosystems, further weakening attempts at coupling biological evolution to evolving early Earth surface chemistry.

Here, we combine high-resolution multiproxy geochemistry to a sequence of well-preserved microbialites in a continuous sedimentary profile of an Ediacaran alkaline lake ecosystem from Morocco's Anti-Atlas30, to demonstrate a strong link between changing environmental chemistry and the assemblage, organization, and succession of a rare extremophilic microbial community. We identify and characterize the discrete biogeochemical intervals that drove their spatial and temporal differentiation into distinct assemblages.

Major elemental composition of the Amane Tazgart microbialites shows two distinct patterns defined by their carbonate and epiclastic chemistries (Supplementary Table 1). The concentrations of major elements generally change across the succession, with major elements correlating positively with Al2O3 in the carbonates, except for MnO and a negative correlation for CaO (Supplementary Fig. 1). Fe2O3, Na2O, TiO2, and P2O5 equally correlate positively with Al2O3 in the epiclastic samples, while the rest show no relationship (Supplementary Fig. 1). In the carbonates, P2O5 content is consistently below the detection limit of 0.1% but above this threshold in the epiclastic samples. The concentrations of As, Ba, Cr, Cu, Pb, Sb, U, V, and Zr are significantly high, while Mo is below the detection limit of 0.5 ppm in most of the carbonates and epiclastic samples, except for four samples (AT2, AT3-C, AT12-S, and AT27), which had values ranging from 0.63 to 15.1 ppm (Supplementary Table 2). The As content decreased upward in the section with a general variation from 108 to 277 ppm, averaging 189 ppm in the carbonates and 22.7–77.2 ppm with an average of 49.62 ppm in the epiclastic deposits (Supplementary Table 2). The trace element systematics are detailed in Supplementary Note 1.

The epiclastic microbialites are typically more enriched in REY than the carbonates compared to PAAS (Fig. 1). Conversely, the carbonate-dominated samples are more enriched in LREE relative to the HREE compared to the epiclastic-dominated samples (Supplementary Note 2) (Fig. 1). Generally, the sampled sequence exhibits positive Eu anomaly (Eu/Eu*) values, varying from 1.04 to 1.41 with a mean of 1.22, and negative Ce anomalies (Ce/Ce*) range from 0.69 to 0.94 with an average of 0.87 (Supplementary Table 3). Y/Ho ratio averages 27.75 for all samples, with values being mostly > 28 and typically associated with the carbonates, except for four samples (AT3C, AT7, AT15-b, and AT17) showing values similar to those recorded in the epiclastic microbialites (Supplementary Table 3). ΣREY varies from 7.63 to 44.47 in the carbonates, increasing to > 23.21 and up to 152.71 in the epiclastic deposits (Supplementary Table 4).

Rare Earth Element and Yttrium spectra (REE + Y). (a) REY spectra of carbonates. (b) REY spectra of Epiclastic deposits.

The concentration of As in all the samples is considerably high, but more enriched in the carbonate than in the epiclastic rocks (Supplementary Note 1, Supplementary Table 5). To evaluate the distribution, location, and sinks of As within the microbialites, nanoscale Synchrotron radiation X-ray micro-Fluorescence (SR-XRF) imaging was performed on several representative carbonate and epiclastic samples, distinguished by fabric similarity and dissimilarity. The investigated samples included iron- and carbonate-rich areas, corresponding to the main microbial products. XRF imaging shows that As is either colocalized with Ca or can be found within iron-rich microbial laminae (Fig. 2). In the carbonate microbial deposits, thrombolites, and composite microbialites, As is mainly concentrated in the peloidal micritic areas within mesoclots (Fig. 2a,b) but rarely associated with the iron-rich mineral phases that are generally detrital in the case of the thrombolitic and composite microbialites. In contrast, stromatolitic fabrics in both the carbonates and epiclastic rocks tend to concentrate As within the iron-rich microbial laminae but not in sparitic laminae (Fig. 2c–f). Strikingly, As can be found within nano-and micro-crystalline carbonate phases, as shown by the pinkish magenta areas in the RGB images (Fig. 2b), corresponding to Ca (blue)-As (red) colocalization. The gray scale image of the As and Ca distribution in Fig. 3a,b shows that the distribution of As reflects the typical crystal shape of the nano-/micro-carbonates, indicating that As is distributed homogenously within these phases. On the other hand, in the iron-rich microbial laminae the distribution of As is heterogenous at the nanometer and micrometer scale containing micrometer-sized As-rich areas and elongated fiber-like structures (Fig. 2f–d). Raman spectroscopy on different types of iron-rich laminae in the epiclastic stromatolites show that the iron-rich laminae are composed of hematite (Supplementary Fig. 2), as is the bulk iron mineralogy of the studied section32. As speciation in the microbial products by Nano-XANES analysis of two different carbonate samples (Figs. 3c, 4a–c), revealed the presence of As(III-O) and As(V–O) within calcite deposits (Fig. 3d) and predominantly As(V–O) in the iron oxides (Fig. 4d).

(a) Thin section photomicrograph of a thrombolite showing a mesoclot characterized by a micritic core (Ms). (b) RGB image of the distribution of As (red), Fe (green), and Ca (blue) in the area indicated by the yellow square in (a), showing the diffused distribution of As within calcite. (c) Thin section photomicrograph of a stromatolite showing an alternation between iron microbial crusts (Ir) and sparitic crusts (Sp). (d) RGB image of the distribution of As, Fe, and Ca in the area indicated by the yellow square in (c), showing the condensed concentration of As within iron microbial crust. (e) Thin section photomicrograph of a stromatolite composed of the alternation of grain-sized microbial laminae (Gs) and sparitic crusts (Sp), (f) RGB images of the distribution of As, Fe, and Ca of grain-sized microbial laminae and a sparitic crust expressing the enrichment of grain-sized laminae in irons and arsenic.

(a) and (b) As and Ca distribution in a carbonate sample, respectively, in (c). Lighter colors indicate higher XRF intensities (higher elemental concentrations). The typical crystal shape of the nano-/micro-carbonates reflected in the carbonate distribution indicates that As is distributed homogenously within these phases. (c) RGB image of the distribution of As, Fe and Ca showing the points where XANES spectra were measured. (d) XANES spectra at the As k-edge were obtained in three representative points from (c).

(a) RGB of As, Fe, and Ca distribution in a carbonate sample showing the association of As within iron, with an inset view of the disseminated distribution of As marked by a white square for (b), corresponding to point 1 of the XANES spectra. (c) RGB image of As, Fe, and Ca distribution measured in a different location of the same sample as (a), showing two points where XANES spectra were performed. (d) XANES spectra at the As K-edge were obtained in three representative points from both (a) and (c).

The projection of the Amane Tazgart samples in the A-CN-K ternary diagram (Fig. 5a) shows CIA values ranging from 35.41 to 51.61 and 38.24 to 58.08 in the carbonates and epiclastic deposits, respectively. Most of the samples deviate from the ideal weathering trend on the A-CN-K diagram, with the tendency to be projected toward the A-K axis suggesting pervasive k-metasomatism may have affected bulk sample composition31,32. The K-corrected CIA values (CIAcorre) range from 34.22 to 58.22 with an average of 47.17 in the carbonates, and from 38.93 to 66 with an average of 55.25 in the epiclastic rocks. They follow the ideal weathering trend (Fig. 5b) and project the samples to a potential andesite protolith. In support of this observation, Th/Sc versus La/Sc discriminant diagrams used to infer parent rock source33,34 (Supplementary Fig. 3), placed the samples in the mafic field, close to an andesite affinity (Supplementary Fig. 3), thus confirming andesite as the likely dominant parent rock source of our samples. This observation is supported by the prevalence of andesitic rocks underlying and overlying the Amane Tazgart succession35 and around the studied lake.

(a) A-CN-K ternary diagram with uncorrected data. (b) A-CN-K ternary diagram with corrected data. CIA–Chemical Index of Alteration, CIAcorr–Corrected CIA, A–Al2O3, CN–CaO* + Na2O, K–K2O (all in molar proportions), CaO*–CaO incorporated in the silicate fraction of the sample. Note: average Upper Continental Crust (UCC; Taylor and McLennan87); Post-Archean average Australian Shale (PAAS: Taylor and McLennan87); Andesite; Granite; TTG and Basalt (Condie88).

The Eu3+/Eu2+ equilibrium in aqueous solutions is viewed as broadly dependent on pressure, pH, speciation, and in particular, temperature, according to theoretical and experimental considerations36,37. Eu3+ dominates in oxidized conditions near the Earth's surface, except in highly reducing and alkaline settings, whereas Eu2+ dominates at temperatures > 250 °C37, where REYN patterns with characteristic positive Eu anomalies are produced38. Therefore, Eu is enriched in Ca-rich and early-magmatic minerals such as plagioclase feldspar39. We produced cross plots of Eu/Eu*, Al, and Y/Ho ratio to discriminate the detrital flux versus hydrothermal origin of the positive Eu anomaly in the studied samples (Supplementary Fig. 4). There is no considerable correlation between Al and Eu/Eu* in our samples (Supplementary Fig. 4a). All the samples plot predominantly in the zone of strong hydrothermal influence with Y/Ho > 25 and Eu/Eu* > 1 (Supplementary Fig. 4b). Together, these observations suggest that the Eu/Eu* anomaly is not overprinted by detrital input and probably derived from hydrothermal fluids. Moreover, the enrichment of light rare earth elements (LREE) over (HREE) in several samples (Supplementary Fig. 5) is consistent with the influence of high temperature hydrothermal fluids characterized by LREE enrichment and positive Eu anomalies40. For comparison, we used a cherty sample from Archean Josefsdal Chert (3.5–3.33 Ga)41, as a reference of microbialites deposited in a restricted shallow marine environment, strongly influenced by hydrothermal activity41,42. The average of carbonate REE + Y patterns displays Josefsdal Chert-like REE + Y trends (Supplementary Fig. 6), supporting the view that hydrothermal fluids were a factor. Although the epiclastic sediments have the same REE + Y pattern, their Eu/Eu* and Ce/Ce* compositions are more attenuated, indicating that hydrothermal fluids have less of an impact on them.

Salinity is an essential chemical feature of water masses that critically regulates homeostatic balance in aquatic biological systems25. It serves as a crucial modulator of microbial biomass volume and biodiversity in alkaline lakes known to possess generally high and variable salinity levels43. Although salinity is challenging to evaluate compared to various palaeoenvironmental parameters in the sedimentary rock record, some elements, such as Y, Ho, B, Ga, Sr, and Ba have been successfully applied as geochemical proxies in the reconstruction of paleoenvironmental salinity condition20,21,22,44,45. In this study, data from the Triassic Zhangjiatan formation46, interpreted as deposited in a fresh lake ecosystem20,46, were employed as a paleosalinity reference.

Y/Ho ratios in chemical sediments such as carbonates are elevated when deposited in saline waters due to the preferential stabilization of Y by anionic salts, resulting in a positive correlation between Y/Ho ratios and salinity45. Our carbonates exhibit an average Y concentration of 5.82 ppm with a tendency to vary from 3.79 to 11.1 ppm, and average Ho content of 0.21 ppm being in the range of 0.127–0.439 ppm. This yielded an average Y/Ho ratio of 28.19 ppm with a spread of 0.07–25.24 ppm (Supplementary Table 3). The epiclastic samples contain an average Y content of 13.3 ppm distributed between 3.79 and 19.5 ppm, being 0.5 ppm for Ho and spanning 0.154–0.705 ppm. This produced a Y/Ho ratio of 26.14 ppm with a spread of 24.57–27.61 (Supplementary Table 3). The difference in the average Y/Ho ratio of up to 2.05 ppm from the carbonates to the epiclastic deposit is consistent with salinity decreasing toward the top of the section. Thus, saltwater conditions are suggested to have been generally prevalent when the carbonates deposited, comparable to seawater conditions. The consistency of the data with the typical high salinity of modern alkaline lakes17,45, points to a temporal transition from saline towards a freshwater ecosystem from the underlying carbonates, up into the stratigraphically younger epiclastic deposit (Fig. 6). There are no Y or Ho data for the Zhangjiatan Formation46.

Stratigraphic profile of enrichment factors (EF) of selected redox trace elements, the chemical index of alteration (CIA), and the Y/Ho ratios as a salinity proxy of Amane Tazgart succession. (EF of Co, Cr, Ni, Th, U, V, and Zn, were calculated relative to the average shale PAAS (Taylor and McLennan87), and As relative to the UCC (McLennan89). The lithostratigraphic column is modified from Chraiki et al.31.

Another proposed bulk mudstone paleosalinity proxy is B versus Ga plots, widely applied in different settings20,21,25, on the assumption that high values of this index are associated with the marine-type origin of sediments, owing to a higher concentration of B in seawater than in freshwater, whereas Ga is unaffected by salinity. Consequently, B/Ga ratios are < 3 in freshwater, 3–6 in brackish, and > 6 in marine-type facies20. The B-Ga biplots for the carbonate and epiclastic samples place the carbonates in the saline and brackish fields and the epiclastic samples in the freshwater field (Fig. 7a), consistent with the Y/Ho results. Comparatively, samples of the Triassic Zhangjiatan Formation plot in the freshwater field, in agreement with the findings of Li et al.46 and Wei, et al.20.

Salinity proxies for Amane Tazgart samples. (a) Biplots of B versus Ga. (b) B versus K (%).

These observations are further supported by the important B/K ratio that has been effectively applied as a salinity proxy on modern sediments22,47 and on deeply buried Paleozoic and Ediacaran fossils22. The B/K ratios of 40 μg/g distinguish marine from non-marine salinities, including riverine and lacustrine deposits, but do not discriminate between brackish and freshwater salinities22,48,49. Our compiled dataset reveals B/K ratios of 17.09–189.96 for carbonates, 2.34–38.8 for epiclastic deposits, and 5.54–23.66 for the lacustrine Zhangjiatan Formation (Supplementary Table 6), suggesting overall that the carbonates were deposited in higher salinity waters, compared to the epiclastic lithologies and the Zhangjiatan Formation that plot in the brackish-freshwater field (Fig. 7b), which agree with the conclusions drawn from the Y/Ho and B/Ga ratios.

Paleosalinity reconstruction based on Sr/Ba ratios is frequently used to differentiate between freshwater, brackish, and full marine facies20,21,46. However, high sediment carbonate content may impact the reliability of this proxy20, as carbonate can raise bulk Sr/Ba ratios in comparison to the clay mineral fraction. The Amane Tazgart samples include up to 48% carbonates, which renders the Sr/Ba ratio unsuitable as a paleosalinity proxy in this case.

Chemical weathering on the continents is assumed to be constrained mainly by variations in atmospheric temperature and humidity. Thus the chemical index of alteration (CIA) was proposed by Nesbitt and Young32 as a proxy for the degree of chemical weathering of soils and sediments, with important inferences for paleoclimate reconstruction50,51. In arid and polar climates, chemical weathering intensity is commonly weak due to cold temperatures and limited precipitation52. Physical weathering or mechanical erosion, for example, due to salt wedging, the formation and thawing of ice, and the grinding of rocks by moving icesheets, prevail over chemical weathering in these settings, resulting in the formation of debris with the same chemical composition as the parent rocks50. On the other hand, chemical reactions between rocks and a solvent, resulting in the selective removal of mobile cations during chemical weathering, are facilitated by warmer and wetter environment conditions. Chemical sediments, therefore, tend to assume chemical compositions that are dissimilar to the parent rocks from which they were sourced and are characterized by higher CIA values.

The K2O content of rocks can be impacted by K-metasomatism, a post-depositional process in which K+ is absorbed from burial fluids by secondary mineral phases in highly altered rocks31,53. The extent of K enrichment can be ascertained from the A-CN-K plot, through deviations from the theoretical alteration trend that must be parallel to the A-CN axis31,50. Any deviation toward the A-K axis represents K addition by K-metasomatism31,50,51.

Chemical weathering significantly affects the mineralogy and major elemental geochemistry of siliciclastic sediments32,54,55. Quantitative parameters, including the CIA32, are therefore theoretically useful in evaluating the degree of chemical weathering. High CIA values generally > 60 reflect the elimination of labile cations (e.g. Ca2+, Na+, K+) compared to stable residual constituents of Al3+, and Ti4+ during weathering31,32,50. Inversely low CIA values being generally < 60, reflect the virtual absence of chemical weathering and thus may indicate cool and/or arid conditions.

The CIA record of the Amane Tazgart succession shows systematic changes closely related to climatic fluctuations (Fig. 6), potentially linked to the end of the preceding Gaskiers glaciations 580–581 million years ago56. While the K-corrected (CIAcorr) for the carbonate deposits presents persistently low values with a mean of 47.17, implying weak chemical weathering typical of a dry and cold climate, the epiclastic deposits show an increasing CIAcorr trend, with values reaching up to 66.42 suggesting warmer and possibly wetter climatic conditions (Figs. 5b, 6, Supplementary Table 7). The predicted arid and wetter climatic conditions in which the carbonates and epiclastic lithologies are deposited, are well correlated with the high salinity lake setting required for the formation of the carbonate microbialites and the brackish to freshwater epiclastic microbialite depositional conditions, respectively (Fig. 6). The preceding cold/dry climatic conditions are consistent with the requirement for intense evaporation to increase water salinity levels during the deposition of carbonate-rich microbialites57. These results are similar to actual modern extreme alkaline lake system conditions where the growth and deposition of carbonate microbialites occur mostly in the dry season linked to higher water evaporation rates57.

The REE composition of carbonates, particularly for microbial and reef carbonates, provides insights into paleo-environmental features58,59. For instance, the Ce/Ce* values of carbonates highlight redox depositional conditions. Nevertheless, Ce anomalies of original aquatic carbonates could be influenced by diagenesis39,58 and detritus58.

Shields and Stille39 state that diagenetic exchange after deposition would influence REE patterns towards becoming progressively more Ce-enriched and Eu-depleted, with decreased DyN/SmN ratios. Such diagenetic exchange would blur or decrease primary features associated with REE distribution patterns and might be recognized by negative correlations between Ce/Ce* versus Eu/Eu* and DyN/SmN39. All samples in this study exhibit small positive Eu anomalies that vary from 1.04 to 1.41 and show a positive correlation with Ce anomalies calculated by using the equation of Lawrence et al.19. The negative correlations between Ce/Ce* versus Eu/Eu* and DyN/SmN, stemming from diagenetic alteration39 do not appear in our samples (Supplementary Fig. 7), suggesting that the calculated Ce anomalies were not significantly altered during diagenesis.

Detrital contamination is one of the most likely causes of deviations in carbonate REE distribution patterns from inherited primary features because of the relative enrichment of silicate derived REE compared to chemically derived carbonate REE content. Such effects could be ascertained by using concentrations and correlations among chemical elements such as bulk rock Al content and Ce/Ce* profiles of leachable carbonates58,59. The Ce/Ce* versus Al biplot demonstrates a moderate positive correlation for the carbonate samples (R2 = 0.6) and no correlation for the epiclastic samples (Fig. 8A). This means that the bulk Ce/Ce* data can be effectively used to infer the redox depositional conditions of the epiclastic samples but less so for the carbonates. Ce/Ce* was calculated and presented according to Bau & Dulski14 (Fig. 8B). As seen, the epiclastic samples cluster tightly in the negative Ce and positive La anomaly fields, with three samples showing deposition in oxic waters.

(A) Biplots of Ce/Ce* eq2 and Al showing a moderate positive correlation for carbonates and no correlation for epiclastic samples. (B) Amane Tazgart carbonates and epiclastic sediments in a graph of Ce/Ce* vs Pr/Pr* (Bau and Dulski15).

Trace element concentrations for U, Th, Cu, Zn, Ni, Co, V, and Cr, as well as Ni/Co, V/V + Ni, V/Cr, Cu/Zn, and U/Th elemental ratios, commonly predict sediment redox depositional condition27,28. Pearson's correlation analysis for the carbonates shows a strong positive correlation between V, Ni, and Al, a moderate positive correlation with Co, Cr, and U, and no correlation with Zn and Cu (Supplementary Table 8). These observations indicate that these elements are not authigenic or of dissolved water column origin, but due to detrital input, except for Zn and Cu. For the epiclastic samples, Pearson´s correlation analysis illustrates a positive correlation between Th, V, Co, and Al, and no correlation with the rest of the redox-sensitive elements (Supplementary Table 9), pointing to the authigenic or water column origin of Cr, Cu, Ni, U, and Zn in the epiclastic samples. Furthermore, and as previously stated, the Amane Tazgart samples are characterized by high EF values for Cu, Zn, and U (Supplementary Table 5), which means that the Cu/Zn ratio can be used as a redox proxy for both the epiclastic and carbonate deposits. The Cu/Zn ratio is proposed as a valid proxy for ancient redox conditions26,60,61,62. In reducing environments, the precipitation of Cu is favored over Zn, while they have about the same solubility in oxidizing environments. Thus, according to Hallberg61, high Cu/Zn ratios indicate reducing depositional conditions, while low ratios suggest highly oxidizing conditions. Except for AT3-C, AT19, and AT35, most of the carbonate samples display low Cu/Zn ratios, ranging from 0 to 0.82 (Fig. 6). These imply likely anoxic environments for AT3-C, AT19, and AT35 and highly oxidizing conditions for the rest of the carbonate samples. The low Cu/Zn ratio in the totality of the epiclastic samples, ranging from 0.05 to 0.64, which is taken to indicate highly oxidizing depositional environments, in full agreement with the Ce anomaly data. In facies formed under anoxic conditions, trace elements frequently show significant enrichment63. Mo EF, Cu EF, Zn EF, Cd EF, and As EF were used as redox proxies in the current study. As EF and Cd EF show the same patterns, expressing high enrichment in trace elements for carbonate samples and moderate enrichment in the epiclastic samples (Fig. 6). Samples with high Cu/Zn ratios are also highly enriched in the cited trace elements compared to the others, consistent with their anoxic depositional environment. The subsurface or pore water conditions were assessed from the early diagenetic phases of samples AT5-V and AT12-S as references (Fig. 9). These cements are highly enriched in redox-sensitive trace elements (Supplementary Table 5) and express high Cu/Zn ratios of 1.6 and 2.51, respectively, which may reflect anoxic pore water interface61,64. Thus, it appears that the carbonate lithologies formed in sub-oxic to oxic conditions and the epiclastic facies in fully oxygenated waters.

Detailed lithostratigraphic column associated with thin section pictures of the different analyzed samples of the Amane Tazgart succession (the lithostratigraphic column is modified from chraiki et al.31). Yellow arrows mark the analyzed parts of the samples (1 corresponds to primary fabrics, 2 correspond to early diagenetic phases).

It is thought that As was a major component of Earth's earliest oceans and that it played a crucial role in biological evolution3,65,66,67. In reduced environments, highly mobile As(III) predominates, conditions believed to have been prevalent in the anoxic Archean biosphere67, a time when As(III) detoxification is thought to have first evolved68. The oxidized As species, As(V), first rose to global prominence during the Paleoproterozoic Great Oxidation Event, (GOE) some 2.45 billion years ago67. Thus, the distribution and prevalence of As(III) and As(V) in the different lithologies have important implications for deducing the depositional environment, redox, and microbial activity.

In this regard, petrographic studies show that the microbialites preserve pristine carbonate phases and primary microbial fabrics69. However, they are slightly affected by early diagenesis, including dissolution and filling vugs with botryoidal and blocky/drusy mosaic calcite and late silicification30. The concentration of As in the studied microbialites presents a diffused pattern observed mainly within the mesoclots affiliated with the thrombolitic and composite microbialites, being concentrated within fine-grained clotted micrite (Fig. 2b). Whereas condensed As patterns are associated mostly with iron-rich stromatolitic microbial mats in the carbonates and epiclastic deposits (Fig. 2c–f). Sparry or botryoidal calcite in the same mapped samples does not concentrate As, which could be explained by the effect of early diagenesis that has likely affected the distribution of As in these diagenetic products, unlike the well preserved primary microbial features.

In line with previous observations above, As speciation in these different mat types is characterized by the co-occurrence of As(III) and As(V) in the calcite patterns (Fig. 3), and the condensed iron-rich As phases mainly in the As(V) oxidation state (Fig. 4). Similarly, microbialites of modern extreme volcanic crater lakes are marked by As(III) enrichment relative to As(V) in diffused calcite patterns, with the predominance of As(V) in globular As-rich phases located within organic matter70. Like the conditions identified in the present study, these crater lakes are characterized by elevated salinity, pH, and elevated As concentrations71. The increased mobility of As(III) in comparison to As(V) under alkaline circumstances was shown to account for this variation in the distribution of As compounds72.

According to laboratory studies, abiotic incorporation of As in calcite structures is mostly associated with As(V)73, while biogenic reactions incorporate As(III) and As(V)74. According to Bardelli et al.,75 bacteria may have been a significant factor in promoting As(III) uptake in the calcite lattice. The presence of As(III) and As(V), linked to the diffused As patterns in the crystalized biogenic calcite particles, strongly indicates that the incorporation of As into the fossilized carbonates was microbially mediated (Figs. 2b, 3). Typically, As(III) enters the cells through aquaporin membrane transporters68,76,77 or can be generated within the cytoplasm when As(V) is trafficked across the cell membrane through phosphate membrane transporters and then reduced to As(III) by As(V) reductases for detoxification67,78. Intracellular As can be further transformed into various methylated arsenicals and associated with biomass to further attenuate toxicity or as biological weapons68,77. In fact, because of its high toxicity, almost all known modern prokaryotes are involved in As detoxification of some form or kind3,71. These include specialized lineages involved in chemoautotrophic and photoautrophic As(III) oxidation, as well as respiratory As(V) reduction to As(III) species71.

Considering this view, the considerably enrichment of As, with an EF of up to 919, relative to UCC in the Ediacaran Amane Tazgart volcanic lake ecosystem can predictably be linked to their hydrothermal depositional environment, redox and biological activity. For example, such extreme As enrichments commonly characterize modern shallow submarine hydrothermal systems where hydrothermal fluids contain > 3000 times more As than seawater79 and alkaline volcanic lakes80. As a consequence, As detoxification is indispensable for microbial communities to succeed in these modern As-rich settings3,71,78,81. Because of the ancient origin of biological As detoxification76 and the abundance of similar As redox species in modern aquatic environments as was then, it is implied that comparable As detoxification mechanisms would have been required for the Amane Tazgart microbialites.

Further still, because As(V) is a strong analog of the phosphate macronutrient, its abundant accumulation in the water column would have posed a serious challenge to extant life, by for example competing with phosphate for uptake into cells through phosphate membrane transporters78. This interaction has been identified as a key cause of phosphate limitation in modern As(V)-rich alkaline volcanic lakes78 and in seawater settings82. Critically, phosphorus is crucial for all life, being the backbone of nucleic acids, and cell membranes, and serves as the cell's main energy trafficking molecule83. In oxic bottom water environments, Fe (III) oxide-oxyhydroxides co-precipitate phosphate from the water column and sediment pore waters at the oxic/anoxic interface84. Below this zone in the anoxic sediments, the Fe (III) oxide-oxyhydroxides are dissolved by reductive processes, and dissolved P is returned to the overlying water column or sediment pore water.

Our study shows that P is absent or below detection in the carbonate samples, which if present should have been accumulated as abundant apatite carbonate minerals, but was instead co-enriched with Fe present principally as hematite30,69 in most of the epiclastic samples (Supplementary Table 1). A well-supported positive linear correlation between Fe2O3 and P2O5 in the epiclastic deposit (Supplementary Fig. 8), suggests P and As sequestration may be linked to Fe (III) oxide-hydroxide precipitation dynamics84, as a function of oxygen availability in the epiclastic stromatolites. The more variable redox conditions in the carbonates imply increased reductive recycling of Fe (III) oxide-oxyhydroxide bound-P and As was likely more common than in the strongly oxygenated epiclastic microbial mats. Consequently, stress resulting from phosphate scarcity would therefore have varied significantly during the deposition of the epiclastic and carbonate microbialites. Particularly, high-affinity phosphate membrane transporters would have been required to scavenge small amounts of available P when dissolved As(V) to phosphate ratios rose dramatically, similar to observations in comparable modern As-rich soda lakes78 and hydrothermal ecosystems3.

Our high-resolution biogeochemical reconstruction demonstrates ecosystem succession by extremophilic Ediacaran microbial communities, driven by climate change, hydrothermal activity, and geochemistry. The underlying thrombolitic microbialites represent extremophiles that subsisted in oxygen-poor As-rich high-temperature hydrothermal saline to brackish waters at a time when the climate was cold/dry. The overlying epiclastic stromatolites evolved in a freshwater body in a warmer/wetter climate. They thrived in an oxygen-rich environment marked by lesser hydrothermal influence. This study provides a rare glimpse into how changing environmental chemistry structured adaptation and shifts in extreme microbial communities towards the end, and possibly, throughout much of Precambrian times.

The Amane Tazgart succession crops out ~ 25 km southward of Ouarzazate city, and in the southern part of the Ediacaran Saghro massif of the Anti-Atlas belt (Supplementary Fig. 9). It is a lens-shaped unit, less than 1 km wide and up to ~ 15 m thick. More details about the geological background of this location can be found in the supplementary file (Supplementary Note 3). To understand the diversity of the Amane Tazgart microbial deposits and the biogeochemical conditions in which they thrived, several primary fabrics from each facies were selected and analyzed from the carbonates (AT2 to AT19, Fig. 9) and epiclastic microbialites (AT37 to AT 28, Fig. 9). The early diagenetic phase was also selected from two samples of carbonate microbialites (AT5 and AT2, Fig. 9) to decipher the physical and chemical influence of pore water during burial diagenesis and also to compare with data from pristine carbonate. Whole rock geochemistry was performed on the epiclastic stromatolites and non-microbialitic samples because they are less affected by diagenesis.

Geochemical analyses were carried out on 13 microbialites and 10 host sediments and analyzed at "Service d’Analyse des Roches et Minéraux" (SARM) of the "Centre de Recherches Pétrographiques et Géochimiques" (CRPG), Nancy, France. Based on a detailed petrographic study, 23 mineralogical phases were cautiously sampled from 21 polished slabs, using a 0.1 to 1 mm diameter microdrill under a stereoscopic microscope to avoid contamination. Each sample was powdered in an agate mortar, and approximately 0.5 g was fused with lithium metaborate (LiBO2) and dissolved in dilute nitric acid (HNO3). Major elements, trace, and rare earth elements data were obtained using inductively coupled plasma atomic emission spectrometry (ICP-AES) and inductively coupled plasma mass spectrometry (ICP-MS), respectively, following the techniques described in Carignan et al.85. Different statistical methods were applied to acquired data, including Pearson's correlation coefficient and principal component analysis (PCA), using the XLSTAT software. Information about the REY data normalization, and the method of calculation of Ce anomalies, Pr anomalies, Eu anomalies, enrichment factor, chemical index of alteration, and statistical analysis can be found in the supplementary file (Supplementary Methods).

Elemental distribution maps were obtained by Synchrotron Radiation scanning X-ray fluorescence (SR-XRF) at the Nanoscopium hard X-ray nanoprobe beamline, Synchrotron Soleil, Paris, France86. Polished sample sections of 0.8–1.2 mm thickness were used for measurements. The monochromatic incident X-ray beam was focused on the sample position by using Kirckpatrick-Baez (KB) nano-focusing mirrors. For XRF imaging, the incident beam energy was tuned to 12.8 keV by the Si(111) double crystal monochromator, above the As K-edge but below the Pb-L3 edge, in order to obtain unambiguous As distribution maps. The correction of the reconstructed elemental distribution maps for XRF spectral overlap, dead-time correction, etc. were performed by Matlab codes developed by the Nanoscopium beamline.

Nano-XANES measurements were performed at the Nanoscopium beamline of Synchrotron Soleil. Two thin Si diodes were used to measure the I0 intensity in front of the KB and the ITr transmitted beam intensity behind the sample, respectively. In addition, a thin Au foil and a dedicated Si diode were inserted behind the Itr monitor. In order to cover the As-K and Au-L3 absorption edges, the energy of the incident X-ray beam was scanned in the 11.825–11.935 keV energy range using 0.5 eV energy steps. This permitted energy calibration of each XANES spectrum compared to the inflection point of the simultaneously measured Au-L3 edge. The As XANES was measured on several characteristic points for samples AT2 and AT4 (Fig. 9).

All data generated or analyzed during this study are included in this published article and supplementary information files.

Bouvier, T. C., Del Giorgio, P. A., Wasniak, T., Sunberg, K. & Alexander, J. Compositional changes in free-living bacterial communities along a salinity gradient in two temperate estuaries. Limnol. Oceanogr. 47(2), 453–470. https://doi.org/10.4319/LO.2002.47.2.0453 (2002).

Article ADS CAS Google Scholar

Herlemann, D. P. R. et al. Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. ISME J. 5(10), 1571–1579. https://doi.org/10.1038/ismej.2011.41 (2011).

Article CAS PubMed PubMed Central Google Scholar

Chi Fru, E. et al. Arsenic and high affinity phosphate uptake gene distribution in shallow submarine hydrothermal sediments. Biogeochemistry 141(1), 41–62. https://doi.org/10.1007/s10533-018-0500-8 (2018).

Article CAS Google Scholar

Allwood, A. C., Kamber, B. S., Walter, M. R., Burch, I. W. & Kanik, I. Trace elements record depositional history of an Early Archean stromatolitic carbonate platform. Chem. Geol. 270(1–4), 148–163. https://doi.org/10.1016/j.chemgeo.2009.11.013 (2010).

Article ADS CAS Google Scholar

Chi, F. E. Iron oxides influence bacterial community structure and the spatial distribution of the aerobic methanotrophs and sulphate reducers in granitic aquifers. Geomicrobiol. J. 26(6), 415–429. https://doi.org/10.1080/01490450902965989 (2009).

Article CAS Google Scholar

Chi, F. E. Copper biogeochemistry: A cornerstone in aerobic methanotrophic bacterial ecology and activity?. Geomicrobiol. J. 28(7), 601–614. https://doi.org/10.1080/01490451.2011.581325 (2011).

Article CAS Google Scholar

Chi Fru, E., Piccinelli, P. & Fortin, D. Insights into the global microbial community structure associated with iron oxyhydroxide minerals deposited in the aerobic biogeosphere. Geomicrobiol. J. 29(7), 587–610. https://doi.org/10.1080/01490451.2011.599474 (2012).

Article CAS Google Scholar

Chi Fru, E. et al. Fossilized iron bacteria reveal a pathway to the biological origin of banded iron formation. Nat. Commun. 4(1), 1–7. https://doi.org/10.1038/ncomms3050 (2013).

Article CAS Google Scholar

Franchi, F. Petrographic and geochemical characterization of the Lower Transvaal Supergroup stromatolitic dolostones (Kanye Basin, Botswana). Precambrian Res. 310, 93–113. https://doi.org/10.1016/j.precamres.2018.02.018 (2018).

Article ADS CAS Google Scholar

Kamber, B. S., Bolhar, R. & Webb, G. E. Geochemistry of late Archaean stromatolites from Zimbabwe: Evidence for microbial life in restricted epicontinental seas. Precambrian Res. 132(4), 379–399. https://doi.org/10.1016/j.precamres.2004.03.006 (2004).

Article ADS CAS Google Scholar

Kamber, B. S. & Webb, G. E. The geochemistry of late Archaean microbial carbonate: Implications for ocean chemistry and continental erosion history. Geochim. Cosmochim. Acta 65(15), 2509–2525. https://doi.org/10.1016/S0016-7037(01)00613-5 (2001).

Article ADS CAS Google Scholar

Deng, S. et al. Microbial dolomite precipitation using sulfate reducing and halophilic bacteria: Results from Qinghai Lake, Tibetan Plateau, NW China. Chem. Geol. 278(3–4), 151–159. https://doi.org/10.1016/J.CHEMGEO.2010.09.008 (2010).

Article ADS CAS Google Scholar

Petrash, D. A., Robbins, L. J., Shapiro, R. S., Mojzsis, S. J. & Konhauser, K. O. Chemical and textural overprinting of ancient stromatolites: Timing, processes, and implications for their use as paleoenvironmental proxies. Precambrian Res. 278, 145–160. https://doi.org/10.1016/j.precamres.2016.03.010 (2016).

Article ADS CAS Google Scholar

Bau, M. & Dulski, P. Distribution of yttrium and rare-earth elements in the Penge and Kuruman iron-formations, Transvaal Supergroup, South Africa. Precambrian Res. 79(1–2), 37–55. https://doi.org/10.1016/0301-9268(95)00087-9 (1996).

Article ADS CAS Google Scholar

Paula-Santos, G. M. et al. Rare earth elements in the terminal Ediacaran Bambuí Group carbonate rocks (Brazil): evidence for high seawater alkalinity during rise of early animals. Precambrian Res. 336, 105506. https://doi.org/10.1016/j.precamres.2019.105506 (2020).

Article ADS CAS Google Scholar

Roy, A., Chakrabarti, G. & Shome, D. Geochemistry of the Neoproterozoic Narji limestone, Cuddapah Basin, Andhra Pradesh, India: implication on palaeoenvironment. Arab. J. Geosci. 11(24), 1–13. https://doi.org/10.1007/s12517-018-4135-9 (2018).

Article CAS Google Scholar

Cao, J. et al. An alkaline lake in the Late Paleozoic Ice Age (LPIA): A review and new insights into paleoenvironment and petroleum geology. Earth-Sci. Rev. 202, 103091. https://doi.org/10.1016/j.earscirev.2020.103091 (2020).

Article CAS Google Scholar

Bau, M. & Dulski, P. Comparing yttrium and rare earths in hydrothermal fluids from the Mid-Atlantic Ridge: Implications for Y and REE behaviour during near-vent mixing and for the Y/Ho ratio of Proterozoic seawater. Chem. Geol. 155(1–2), 77–90. https://doi.org/10.1016/S0009-2541(98)00142-9 (1999).

Article ADS CAS Google Scholar

Lawrence, M. G., Greig, A., Collerson, K. D. & Kamber, B. S. Rare earth element and yttrium variability in South East Queensland waterways. Aquat. Geochem. 12(1), 39–72. https://doi.org/10.1007/s10498-005-4471-8 (2006).

Article CAS Google Scholar

Wei, W. & Algeo, T. J. Elemental proxies for paleosalinity analysis of ancient shales and mudrocks. Geochim. Cosmochim. Acta 287, 341–366. https://doi.org/10.1016/j.gca.2019.06.034 (2020).

Article ADS CAS Google Scholar

Wei, W. et al. Identifying marine incursions into the Paleogene Bohai Bay Basin lake system in northeastern China. Int. J. Coal Geol. 200(May), 1–17. https://doi.org/10.1016/j.coal.2018.10.001 (2018).

Article ADS CAS Google Scholar

Retallack, G., Gates, N., York, N. & Feldshuhorn, M. Boron paleosalinity proxy for deeply buried Paleozoic and Ediacaran fossils. Palaeogeogr. Palaeoclimatol. Palaeoecol. 540, 109536. https://doi.org/10.1016/j.palaeo.2019.109536 (2020).

Article Google Scholar

Wang, A., Wang, Z., Liu, J., Xu, N. & Li, H. The Sr/Ba ratio response to salinity in clastic sediments of the Yangtze River Delta. Chem. Geol. 559, 119923. https://doi.org/10.1016/j.chemgeo.2020.119923 (2021).

Article ADS CAS Google Scholar

Goldschmidt, V. M., Peters, C. L. & Hauptmann, H. Zur Geochemie des Bors. Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Math Klasse 1932, 402–407 (1932).

Google Scholar

Jewuła, K., Środoń, J., Kuligiewicz, A., Mikołajczak, M. & Liivamägi, S. Critical evaluation of geochemical indices of palaeosalinity involving boron. Geochim. Cosmochim. Acta 322, 1–23. https://doi.org/10.1016/j.gca.2022.01.027 (2022).

Article ADS CAS Google Scholar

Abedini, A. & Calagari, A. A. Geochemistry of claystones of the ruteh formation, NW Iran: Implications for provenance, source-Area weathering, and paleo-redox conditions. Neues Jahrb fur Miner. Abhandlungen 194(2), 107–123. https://doi.org/10.1127/njma/2017/0040 (2017).

Article CAS Google Scholar

Jones, B. & Manning, D. A. C. Comparison of geochemical indices used for the interpretation of palaeoredox conditions in ancient mudstones. Chem. Geol. 111(1–4), 111–129. https://doi.org/10.1016/0009-2541(94)90085-X (1994).

Article ADS Google Scholar

Armstrong-Altrin, J. S., Nagarajan, R., Balaram, V. & Natalhy-Pineda, O. Petrography and geochemistry of sands from the Chachalacas and Veracruz beach areas, western Gulf of Mexico, Mexico: Constraints on provenance and tectonic setting. J. S. Am. Earth Sci. 64, 199–216. https://doi.org/10.1016/J.JSAMES.2015.10.012 (2015).

Article CAS Google Scholar

Ionescu, D. et al. Calcium dynamics in microbialite-forming exopolymer-rich mats on the atoll of Kiritimati, Republic of Kiribati, Central Pacific. Geobiology 13(2), 170–180. https://doi.org/10.1111/gbi.12120 (2015).

Article CAS PubMed Google Scholar

Chraiki, I. et al. A 571 million-year-old alkaline volcanic lake photosynthesizing microbial community, the Anti-atlas, Morocco. Geobiology https://doi.org/10.1111/gbi.12425 (2020).

Article PubMed Google Scholar

Fedo, C. M., Nesbitt, H. W. & Young, G. M. Unraveling the effects of potassium metasomatism in sedimentary-rocks and paleosols, with implications for paleoweathering conditions and provenance unraveling the effects of potassium metasomatism in sedimentary rocks and paleosols, with implications. Geology 23(October), 921–924. https://doi.org/10.1130/0091-7613(1995)023%3c0921 (1995).

Article ADS CAS Google Scholar

Nesbitt, H. W. & Young, G. M. Early proterozoic climates and plate motions inferred from major element chemistry of lutites. Nature 299(5885), 715–717. https://doi.org/10.1038/299715a0 (1982).

Article ADS CAS Google Scholar

Bankole, O. M., El Albani, A., Meunier, A., Poujol, M. & Bekker, A. Elemental geochemistry and Nd isotope constraints on the provenance of the basal siliciclastic succession of the middle Paleoproterozoic Francevillian Group, Gabon. Precambrian Res. https://doi.org/10.1016/j.precamres.2020.105874 (2020).

Article Google Scholar

Bhatia, M. R. & Crook, K. A. W. Trace element characteristics of graywackes and tectonic setting discrimination of sedimentary basins. Contrib. Miner. Petrol. 92(2), 181–193. https://doi.org/10.1007/BF00375292 (1986).

Article ADS CAS Google Scholar

Walsh, G. J. et al. Neoproterozoic tectonic evolution of the Jebel Saghro and Bou Azzer—El Graara inliers, eastern and central Anti-Atlas, Morocco. Precambrian Res. 216–219, 23–62. https://doi.org/10.1016/J.PRECAMRES.2012.06.010 (2012).

Article ADS Google Scholar

Möller, P. & Bau, M. Rare-earth patterns with positive cerium anomaly in alkaline waters from Lake Van, Turkey. Earth Planet. Sci. Lett. 117(3–4), 671–676. https://doi.org/10.1016/0012-821X(93)90110-U (1993).

Article ADS Google Scholar

Sverjensky, D. A. Europium redox equilibria in aqueous solution. Earth Planet. Sci. Lett. 67(1), 70–78. https://doi.org/10.1016/0012-821X(84)90039-6 (1984).

Article ADS CAS Google Scholar

Bau, M. Rare-earth element mobility during hydrothermal and metamorphic fluid-rock interaction and the significance of the oxidation state of europium. Chem. Geol. 93(3–4), 219–230. https://doi.org/10.1016/0009-2541(91)90115-8 (1991).

Article ADS CAS Google Scholar

Shields, G. & Stille, P. Diagenetic constraints on the use of cerium anomalies as palaeoseawater redox proxies: An isotopic and REE study of Cambrian phosphorites. Chem. Geol. 175(1–2), 29–48. https://doi.org/10.1016/S0009-2541(00)00362-4 (2001).

Article ADS CAS Google Scholar

Zhao, Y. et al. Rare earth element geochemistry of carbonates as a proxy for deep-time environmental reconstruction. Palaeogeogr. Palaeoclimatol. Palaeoecol. 574, 110443. https://doi.org/10.1016/J.PALAEO.2021.110443 (2021).

Article Google Scholar

Hubert A. Chemical and mineralogical signatures of oxygenic photosynthesis in Archean and Paleoproterozoic sediments (2015).

Westall, F. et al. A hydrothermal-sedimentary context for the origin of life. Astrobiology 18(3), 259–293. https://doi.org/10.1089/ast.2017.1680 (2018).

Article ADS CAS PubMed PubMed Central Google Scholar

Pecoraino, G. & Alessandro, W. D. The other side of the coin: Geochemistry of Alkaline lakes in volcanic areas. Sedimentology https://doi.org/10.1007/978-3-642-36833-2 (2015).

Article Google Scholar

Xia, L. et al. Unsynchronized evolution of salinity and pH of a Permian alkaline lake influenced by hydrothermal fluids: A multi-proxy geochemical study. Chem. Geol. https://doi.org/10.1016/j.chemgeo.2020.119581 (2020).

Article Google Scholar

Lawrence, M. G. & Kamber, B. S. The behaviour of the rare earth elements during estuarine mixing-revisited. Mar. Chem. 100(1–2), 147–161. https://doi.org/10.1016/j.marchem.2005.11.007 (2006).

Article CAS Google Scholar

Li, D. et al. Elemental characteristics and paleoenvironment reconstruction: A case study of the Triassic lacustrine Zhangjiatan oil shale, southern Ordos Basin, China. Acta Geochim. 37(1), 134–150. https://doi.org/10.1007/s11631-017-0193-z (2018).

Article ADS CAS Google Scholar

Walker, C. T. & Price, N. B. Departure curves for computing paleosalinity from boron in Illites and shales. Aapg Bull. 5, 833–841 (1963).

Google Scholar

Chetelat, B. et al. Boron isotopes geochemistry of the Changjiang basin rivers. Geochim. Cosmochim. Acta 73(20), 6084–6097. https://doi.org/10.1016/J.GCA.2009.07.026 (2009).

Article ADS CAS Google Scholar

Ericson, D. B., Ewing, M., Wollin, G. & Heezen, B. C. Atlantic deep-sea sediment cores|GSA Bulletin|GeoScienceWorld. GSA Bull. 72, 193–286. https://doi.org/10.1130/0016-7606(1961)72[193:ADSC]2.0.CO;2 (1961).

Article CAS Google Scholar

Wang, P. et al. The chemical index of alteration (CIA) as a proxy for climate change during glacial-interglacial transitions in Earth history. Earth-Sci. Rev. 201, 103032. https://doi.org/10.1016/j.earscirev.2019.103032 (2020).

Article CAS Google Scholar

Sun, S. et al. Rapid climatic fluctuations during the Guadalupian-Lopingian transition: Implications from weathering indices recorded in acid-insoluble residues of carbonate rocks, South China. J. Asian Earth Sci. 230, 105222. https://doi.org/10.1016/J.JSEAES.2022.105222 (2022).

Article Google Scholar

Nesbitt, S. W. & Zipser, E. J. The Diurnal cycle of rainfall and convective intensity according to three years of TRMM measurements. J. Clim. 16(10), 1456–1475. https://doi.org/10.1175/1520-0442(2003)016 (2003).

Article ADS Google Scholar

Panahi, A., Young, G. M. & Rainbird, R. H. Behavior of major and trace elements (including REE) during Paleoproterozoic pedogenesis and diagenetic alteration of an Archean granite near Ville Marie, Québec, Canada. Geochim. Cosmochim. Acta 64(13), 2199–2220. https://doi.org/10.1016/S0016-7037(99)00420-2 (2000).

Article ADS CAS Google Scholar

Johnsson, M. J., Stallard, R. F. & Meade, R. H. First-cycle quartz arenites in the Orinoco river basin, Venezuela and Colombia. Univ. Chic. Press J. 96(3), 263–277. https://doi.org/10.1086/629219 (1988).

Article Google Scholar

McLennan, S. M. Weathering and global denudation. Univ. Chic. Press J. 101(2), 295–303. https://doi.org/10.1086/648222 (1993).

Article Google Scholar

Pu, J. P. et al. Dodging snowballs: Geochronology of the Gaskiers glaciation and the first appearance of the Ediacaran biota. Geology 44(11), 955–958. https://doi.org/10.1130/G38284.1 (2016).

Article ADS CAS Google Scholar

Sanz-Montero, M. E., Cabestrero, Ó. & Sánchez-Román, M. Microbial Mg-rich carbonates in an extreme alkaline lake (Las Eras, Central Spain). Front. Microbiol. 10, 148. https://doi.org/10.3389/fmicb.2019.00148 (2019).

Article PubMed PubMed Central Google Scholar

Nothdurft, L. D., Webb, G. E. & Kamber, B. S. Rare earth element geochemistry of Late Devonian reefal carbonates, Canning Basin, Western Australia: Confirmation of a seawater REE proxy in ancient limestones. Geochim. Cosmochim. Acta 68(2), 263–283. https://doi.org/10.1016/S0016-7037(03)00422-8 (2004).

Article ADS CAS Google Scholar

Webb, G. E. & Kamber, B. S. Rare earth elements in Holocene reefal microbialites: A new shallow seawater proxy. Geochim. Cosmochim. Acta 64(9), 1557–1565. https://doi.org/10.1016/S0016-7037(99)00400-7 (2000).

Article ADS CAS Google Scholar

Adebayo, O. F., Akinyemi, S. A., Madukwe, H. Y., Aturamu, A. O. & Ojo, A. O. Geochemical characterization and palynological studies of some Agbada Formation deposits of the Niger Delta basin: Implications for paleodepositional environments. Turk. J. Earth Sci. 25(6), 573–591. https://doi.org/10.3906/yer-1512-8 (2016).

Article CAS Google Scholar

Hallberg, R. O. A geochemical method for investigation of paleoredox conditions in sediment. Ambio Spec. Rep. 4(4), 139–147 (1976).

CAS Google Scholar

Nagarajan, R., Madhavaraju, J., Nagendra, R., Selvamony Armstrong-Altrin, J. & Moutte, J. Geochemistry of neoproterozoic shales of the Rabanpalli formation, Bhima Basin, Northern Karnataka, Southern India: Implications for provenance and paleoredox conditions. Rev. Mex. Cienc. Geol. 24, 150–160 (2007).

Google Scholar

Algeo, T. J. & Maynard, J. B. Trace-element behavior and redox facies in core shales of Upper Pennsylvanian Kansas-type cyclothems. Chem. Geol. 206(3–4), 289–318. https://doi.org/10.1016/J.CHEMGEO.2003.12.009 (2004).

Article ADS CAS Google Scholar

Little, S. H., Vance, D., Lyons, T. W. & McManus, J. Controls on trace metal authigenic enrichment in reducing sediments: Insights from modern oxygen-deficient settings. Am. J. Sci. 315(2), 77–119. https://doi.org/10.2475/02.2015.01 (2015).

Article ADS CAS Google Scholar

Chen, S. C. et al. The great oxidation event expanded the genetic repertoire of arsenic metabolism and cycling. Proc. Natl. Acad. Sci. U. S. A. 117(19), 10414–10421. https://doi.org/10.1073/PNAS.2001063117/SUPPL_FILE/PNAS.2001063117.SD03.TXT (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Chi Fru, E. et al. Arsenic stress after the Proterozoic glaciations. Sci. Rep. https://doi.org/10.1038/srep17789 (2015).

Article PubMed Central Google Scholar

Chi Fru, E. et al. The rise of oxygen-driven arsenic cycling at ca. 2.48 Ga. Geology 47(3), 243–246 (2019).

Article ADS Google Scholar

Chen, J. & Rosen, B. P. The Arsenic methylation cycle: How microbial communities adapted methylarsenicals for use as weapons in the continuing war for dominance. Front. Environ. Sci. 8, 43. https://doi.org/10.3389/FENVS.2020.00043/BIBTEX (2020).

Article CAS Google Scholar

Chraiki, I., Bouougri, E. H. & El Albani, A. Microbialites diversity from the Ediacaran of the Anti-Atlas (Morocco): A snapshot of microbial oases thriving in an alkaline volcanic lake. Ann. Paléontol. 108(4), 102584. https://doi.org/10.1016/J.ANNPAL.2022.102584 (2023).

Article Google Scholar

Sancho-Tomás, M. et al. Geochemical evidence for arsenic cycling in living microbialites of a high altitude Andean Lake (Laguna Diamante, Argentina). Chem. Geol. https://doi.org/10.1016/j.chemgeo.2020.119681 (2020).

Article Google Scholar

Rascovan, N., Maldonado, J., Vazquez, M. P. & Eugenia, F. M. Metagenomic study of red biofilms from Diamante Lake reveals ancient arsenic bioenergetics in haloarchaea. ISME J. 10(2), 299–309. https://doi.org/10.1038/ISMEJ.2015.109 (2016).

Article CAS PubMed Google Scholar

Smedley, P. L. & Kinniburgh, D. G. A review of the source, behaviour and distribution of arsenic in natural waters. Appl. Geochem. 17(5), 517–568. https://doi.org/10.1016/S0883-2927(02)00018-5 (2002).

Article ADS CAS Google Scholar

Yokoyama, Y., Tanaka, K. & Takahashi, Y. Differences in the immobilization of arsenite and arsenate by calcite. Geochim. Cosmochim. Acta 91, 202–219. https://doi.org/10.1016/J.GCA.2012.05.022 (2012).

Article ADS CAS Google Scholar

Catelani, T. et al. Arsenic uptake in bacterial calcite. Geochim. Cosmochim. Acta 222, 642–654. https://doi.org/10.1016/J.GCA.2017.11.013 (2018).

Article ADS CAS Google Scholar

Bardelli, F. et al. Arsenic uptake by natural calcite: An XAS study. Geochim. Cosmochim. Acta 75(11), 3011–3023. https://doi.org/10.1016/J.GCA.2011.03.003 (2011).

Article ADS CAS Google Scholar

Chen, S. et al. Simultaneous removal of para-arsanilic acid and the released inorganic arsenic species by CuFe2O4 activated peroxymonosulfate process. Sci. Total Environ. 742, 140587. https://doi.org/10.1016/J.SCITOTENV.2020.140587 (2020).

Article ADS CAS PubMed Google Scholar

Li, Y. et al. Arsenic and antimony co-contamination influences on soil microbial community composition and functions: Relevance to arsenic resistance and carbon, nitrogen, and sulfur cycling. Environ. Int. 153, 106522. https://doi.org/10.1016/J.ENVINT.2021.106522 (2021).

Article CAS PubMed Google Scholar

Elias, M. et al. The molecular basis of phosphate discrimination in arsenate-rich environments. Nature 491(7422), 134–137. https://doi.org/10.1038/nature11517 (2012).

Article ADS CAS PubMed Google Scholar

Breuer, C. & Pichler, T. Arsenic in marine hydrothermal fluids. Chem. Geol. 348, 2–14. https://doi.org/10.1016/J.CHEMGEO.2012.10.044 (2013).

Article ADS CAS Google Scholar

Lloyd, J. R. & Oremland, R. S. Microbial transformations of arsenic in the environment: From Soda lakes to aquifers. Elements 2(2), 85–90. https://doi.org/10.2113/GSELEMENTS.2.2.85 (2006).

Article CAS Google Scholar

Kurth, D. et al. Arsenic metabolism in high altitude modern stromatolites revealed by metagenomic analysis /631/326/47/4112 /631/326/2565/2142 /45/23 /45/22 article. Sci. Rep. 7(1), 1–16. https://doi.org/10.1038/s41598-017-00896-0 (2017).

Article CAS Google Scholar

Hashihama, F. et al. Arsenate and microbial dynamics in different phosphorus regimes of the subtropical Pacific Ocean. Prog. Oceanogr. 176, 102115. https://doi.org/10.1016/J.POCEAN.2019.05.007 (2019).

Article Google Scholar

Smith, E. M. & Prairie, Y. T. Bacterial metabolism and growth efficiency in lakes: The importance of phosphorus availability. Limnol. Oceanogr. 49(1), 137–147. https://doi.org/10.4319/LO.2004.49.1.0137 (2004).

Article ADS CAS Google Scholar

Hemmingsson, C. et al. Evaluation of phosphate-uptake mechanisms by Fe(III) (oxyhydr)oxides in early proterozoic oceanic conditions. Environ. Chem. 15(2), 18–28. https://doi.org/10.1071/EN17124 (2018).

Article CAS Google Scholar

Carignan, J., Hild, P., Mevelle, G., Morel, J. & Yeghicheyan, D. Routine analyses of trace elements in geological samples using flow injection and low pressure on-line liquid chromatography coupled to ICP-MS: A study of geochemical reference materials BR, DR-N, UB-N AN-G and GH. Geostand Newsl. 25(2–3), 187–198. https://doi.org/10.1111/J.1751-908X.2001.TB00595.X (2001).

Article CAS Google Scholar

Somogyi, A. et al. Optical design and multi-length-scale scanning spectro-microscopy possibilities at the Nanoscopium beamline of Synchrotron Soleil. J. Synchrotron Radiat. 22(4), 1118–1129. https://doi.org/10.1107/S1600577515009364 (2015).

Article CAS PubMed Google Scholar

Taylor, S R; & McLennan SM. The Continental Crust: Its Composition and Evolution. United States. united states. https://www.osti.gov/biblio/6582885 (1985).

Condie, K. C. Chemical composition and evolution of the upper continental crust: Contrasting results from surface samples and shales. Chem. Geol. 104, 1–37 (1993).

Article ADS CAS Google Scholar

McLennan, S. M. Relationships between the trace element composition of sedimentary rocks and upper continental crust. Geochem. Geophys. Geosyst. https://doi.org/10.1029/[email protected]/(ISSN)1525-2027.GERM1 (2001).

Article Google Scholar

Download references

This work was supported and carried out under the Scientific project funded by Hassan II Academy of Sciences and Technology- Morocco (Project Leader, Prof. El Albani). Support is acknowledged from "La Région Nouvelle Aquitaine", University of Poitiers, Synchrotron Soleil-CNRS, Paris Saclay, and Cadi Ayyad University, the Faculty of Sciences-Semlalia, Marrakech. We also like to thank Céderic Lebailly and Nadia Guignard from the University of Poitiers for their technical assistance.

Department of Geology, Faculty of Sciences-Semlalia, Cadi Ayyad University, Marrakesh, Morocco

Ibtissam Chraiki & El Hafid Bouougri

Centre for Geobiology and Geochemistry, College of Physical and Engineering Sciences, School of Earth and Environmental Sciences, Cardiff University, Cardiff, UK

Ernest Chi Fru

Nanoscopium Beamline Synchrotron Soleil, Saint-Aubin, 91192, Gif-sur-Yvette, France

Andrea Somogyi

CNRS IC2MP UMR 7285, University of Poitiers, Poitiers, France

Olabode Bankole, Mohamed Ghnahalla & Abderrazak El Albani

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

I.C., E.H.B. and A.E.A. designed the study, I.C. and A.S. carried out and treated Synchrotron analysis, I.C., E.C.F., A.E.A., E.H.B., O.B. and M.G. analyzed the geochemical data, I.C. wrote the manuscript with the contribution of the co-authors.

Correspondence to Abderrazak El Albani.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Chraiki, I., Chi Fru, E., Somogyi, A. et al. Blooming of a microbial community in an Ediacaran extreme volcanic lake system. Sci Rep 13, 9080 (2023). https://doi.org/10.1038/s41598-023-36031-5

Download citation

Received: 05 February 2023

Accepted: 27 May 2023

Published: 05 June 2023

DOI: https://doi.org/10.1038/s41598-023-36031-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.