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

Ibtissam Chraiki, Ernest Chi Fru, Andrea Somogyi, El Hafid Bouougri, Olabode Bankole, Mohamed Ghnahalla, Abderrazak El Albani

► To cite this version:

Ibtissam Chraiki, Ernest Chi Fru, Andrea Somogyi, El Hafid Bouougri, Olabode Bankole, et al. Blooming of a microbial community in an Ediacaran extreme volcanic lake system. Scientific Reports, 2023, 13 (1), pp.9080. 10.1038/s41598-023-36031-5. hal-04157601

HAL Id: hal-04157601

https://hal.science/hal-04157601

Submitted on 12 Jul 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

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

Ibtissam Chraiki<fnref n="1" />, Ernest Chi Fru², Andrea Somogyi³, El Hafid Bouougri<fnref n="1" />, Olabode Bankole⁴, Mohamed Ghnahalla⁴ & Abderrazak El Albani⁴✉

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 gradients<fnref n="1" />⁻³ that tend to precipitate diagnostic biogenic minerals such as carbonates, iron oxides, and pyrite³⁻¹². 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 deposition¹¹,¹³. 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 conditions¹⁴; and in certain instances, non-marine aquatic ecosystems influenced by elevated alkalinity¹⁵. 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 ecosystems¹⁵.

On the other hand, the relative abundance of europium (Eu) is often associated with hydrothermal fluid enrichment, feldspars, and diagenesis¹⁶,¹⁷. While negative Eu anomalies may represent the original composition of the aquatic environment from which limestones formed¹⁶, and La, gadolinium (Gd), and Y anomalies, mineral complexation phenomena unique to the hydrosphere¹⁸. 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 basis¹⁹⁻²⁴. 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 freshwater²⁰,²³; (2) Sr/Ba is > 0.5 characterizes marine deposits, 0.2–0.5 brackishness, and <0.2 freshwater²⁰; and (3) B/K > 40 µg/g distinguishes marine from brackish/freshwater sediments²²,²⁵. 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 conditions²⁶⁻²⁸. However²⁹, 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

<footnote n="1">¹Department of Geology, Faculty of Sciences-Semlalia, Cadi Ayyad University, Marrakesh, Morocco. ²Centre for Geobiology and Geochemistry, College of Physical and Engineering Sciences, School of Earth and Environmental Sciences, Cardiff University, Cardiff, UK. ³Nanoscopium Beamline Synchrotron Soleil, Saint-Aubin, 91192 Gif-sur-Yvette, France. ⁴CNRS IC2MP UMR 7285, University of Poitiers, Poitiers, France. ✉email: abder.albani@univ-poitiers.fr</footnote>
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-Atlas³⁰, 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.

## Results

**Whole rock geochemistry.** 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 Al₂O₃ in the carbonates, except for MnO and a negative correlation for CaO (Supplementary Fig. 1). Fe₂O₃, Na₂O, TiO₂, and P₂O₅ equally correlate positively with Al₂O₃ in the epiclastic samples, while the rest show no relationship (Supplementary Fig. 1). In the carbonates, P₂O₅ 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).

**Arsenic intercorrelation with other elements.** 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

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Rare Earth Element and Yttrium spectra (REE + Y). (a) REY spectra of carbonates. (b) REY spectra of Epiclastic deposits.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** (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.</figcaption></figure>

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).
<figure><img src="image_5.png" /><figcaption>**Figure 3.** (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).</figcaption></figure>

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 section³². 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).

**Chemical weathering and sediment provenance.** 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 composition³¹,³². The K-corrected CIA values (CIA<sub>correct</sub>) 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 source³³,³⁴ (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

<figure><img src="image_6.png" /><figcaption>**Figure 4.** (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).</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** (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–Al₂O₃, CN–CaO* + Na₂O, K–K₂O (all in molar proportions), CaO*–CaO incorporated in the silicate fraction of the sample. Note: average Upper Continental Crust (UCC; Taylor and McLennan⁸⁷); Post-Archean average Australian Shale (PAAS: Taylor and McLennan⁸⁷); Andesite; Granite; TTG and Basalt (Condie⁸⁸).</figcaption></figure>

by the prevalence of andesitic rocks underlying and overlying the Amane Tazgart succession³⁵ and around the studied lake.

Discussion

The Eu³⁺/Eu²⁺ equilibrium in aqueous solutions is viewed as broadly dependent on pressure, pH, speciation, and in particular, temperature, according to theoretical and experimental considerations³⁶,³⁷. Eu³⁺ dominates in oxidized conditions near the Earth's surface, except in highly reducing and alkaline settings, whereas Eu²⁺ dominates at temperatures > 250 °C³⁷, where REYₙ patterns with characteristic positive Eu anomalies are produced³⁸. Therefore, Eu is enriched in Ca-rich and early-magmatic minerals such as plagioclase feldspar³⁹. 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 anomalies⁴⁰. For comparison, we used a cherty sample from Archean Josefsdal Chert (3.5–3.33 Ga)⁴¹, as a reference of microbialites deposited in a restricted shallow marine environment, strongly influenced by hydrothermal activity⁴¹,⁴². 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 systems²⁵. It serves as a crucial modulator of microbial biomass volume and biodiversity in alkaline lakes known to possess generally high and variable salinity levels⁴³. 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 condition²⁰⁻²²,⁴⁴,⁴⁵. In this study, data from the Triassic Zhangjiatan formation⁴⁶, interpreted as deposited in a fresh lake ecosystem²⁰,⁴⁶, 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 salinity⁴⁵. 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 lakes¹⁷,⁴⁵, points to a temporal transition from
<figure><img src="image_8.png" /><figcaption>**Figure 6.** 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 McLennan<sup>87</sup>), and As relative to the UCC (McLennan<sup>89</sup>). The lithostratigraphic column is modified from Chraiki et al.<sup>31</sup></figcaption></figure>

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 Formation<sup>46</sup>.

Another proposed bulk mudstone paleosalinity proxy is B versus Ga plots, widely applied in different
settings<sup>20,21,25</sup>, on the assumption that high values of this index are associated with the marine-type origin of sedi-
ments, 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 facies<sup>20</sup>. 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.<sup>46</sup> and Wei, et al.<sup>20</sup>.

These observations are further supported by the important B/K ratio that has been effectively applied as a salinity proxy on modern sediments<sup>22,47</sup> and on deeply buried Paleozoic and Ediacaran fossils<sup>22</sup>. 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 salinities<sup>22,48,49</sup>. 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 facies<sup>20,21,46</sup>. However, high sediment carbonate content may impact the reliability of

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Salinity proxies for Amane Tazgart samples. (a) Biplots of B versus Ga. (b) B versus K (%).</figcaption></figure>
this proxy²⁰, 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 Young³² as a proxy for the degree of chemical weathering of soils and sediments, with important inferences for paleoclimate reconstruction⁵⁰,⁵¹. In arid and polar climates, chemical weathering intensity is commonly weak due to cold temperatures and limited precipitation⁵². 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 rocks⁵⁰. 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 K₂O 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 rocks³¹,⁵³. 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 axis³¹,⁵⁰. Any deviation toward the A-K axis represents K addition by K-metasomatism³¹,⁵⁰,⁵¹.

Chemical weathering significantly affects the mineralogy and major elemental geochemistry of siliciclastic sediments³²,⁵⁴,⁵⁵. Quantitative parameters, including the CIA³², are therefore theoretically useful in evaluating the degree of chemical weathering. High CIA values generally > 60 reflect the elimination of labile cations (e.g. Ca²⁺, Na⁺, K⁺) compared to stable residual constituents of Al³⁺, and Ti⁴⁺ during weathering³¹,³²,⁵⁰. 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 ago⁵⁶. While the K-corrected (CIA<sub>corr</sub>) 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 CIA<sub>corr</sub> 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 microbialites⁵⁷. 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 rates⁵⁷.

The REE composition of carbonates, particularly for microbial and reef carbonates, provides insights into paleo-environmental features⁵⁸,⁵⁹. For instance, the Ce/Ce* values of carbonates highlight redox depositional conditions. Nevertheless, Ce anomalies of original aquatic carbonates could be influenced by diagenesis³⁹,⁵⁸ and detritus⁵⁸.

Shields and Stille³⁹ state that diagenetic exchange after deposition would influence REE patterns towards becoming progressively more Ce-enriched and Eu-depleted, with decreased Dy<sub>N</sub>/Sm<sub>N</sub> 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 Dy<sub>N</sub>/Sm<sub>N</sub>³⁹. 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.¹⁹. The negative correlations between Ce/Ce* versus Eu/Eu* and Dy<sub>N</sub>/Sm<sub>N</sub>, stemming from diagenetic alteration³⁹ 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 carbonates⁵⁸,⁵⁹. The Ce/Ce* versus Al biplot demonstrates a moderate positive correlation for the carbonate samples (R² = 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 & Dulski¹⁴ (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.

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 condition²⁷,²⁸. 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 conditions²⁶,⁶⁰⁻⁶². In reducing environments, the precipitation of Cu is favored over Zn, while they have about the same solubility
<figure><img src="image_10.png" /><figcaption>**Figure 8.** (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 Dulski<sup>15</sup>).</figcaption></figure>

in oxidizing environments. Thus, according to Hallberg<sup>61</sup>, 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 enrichment<sup>63</sup>. 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 interface<sup>61,64</sup>. Thus, it appears that the carbonate lithologies formed in sub-oxic to oxic conditions and the epiclastic facies in fully oxygenated waters.

It is thought that As was a major component of Earth's earliest oceans and that it played a crucial role in biological evolution<sup>3,65-67</sup>. In reduced environments, highly mobile As(III) predominates, conditions believed to have been prevalent in the anoxic Archean biosphere<sup>67</sup>, a time when As(III) detoxification is thought to have first evolved<sup>68</sup>. The oxidized As species, As(V), first rose to global prominence during the Paleoproterozoic Great Oxidation Event, (GOE) some 2.45 billion years ago<sup>67</sup>. 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 fabrics<sup>69</sup>. However, they are slightly affected by early diagenesis, including dissolution and filling vugs with botryoidal and blocky/drusy mosaic calcite and late silicification<sup>30</sup>. 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 matter<sup>70</sup>. Like the conditions identified in the present study, these crater lakes are characterized by elevated salinity, pH, and elevated As concentrations<sup>71</sup>. 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 compounds<sup>72</sup>.

According to laboratory studies, abiotic incorporation of As in calcite structures is mostly associated with As(V)<sup>73</sup>, while biogenic reactions incorporate As(III) and As(V)<sup>74</sup>. According to Bardelli et al.,<sup>75</sup> 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 transporters<sup>68,76,77</sup> or can be generated within the cytoplasm when As(V)
<figure><img src="image_11.png" /><figcaption>**Legende**

Carbonate

Calcarenite

Sandstones

Conglomerates

Sandy limestones

Silty sands

Clastic stromatolites

Peperites</figcaption></figure>

<caption>**Figure 9.** 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.³¹). Yellow arrows mark the analyzed parts of the samples (1 corresponds to primary fabrics, 2 correspond to early diagenetic phases).</caption>
is trafficked across the cell membrane through phosphate membrane transporters and then reduced to As(III) by As(V) reductases for detoxification⁶⁷,⁷⁸. Intracellular As can be further transformed into various methylated arsenicals and associated with biomass to further attenuate toxicity or as biological weapons⁶⁸,⁷⁷. In fact, because of its high toxicity, almost all known modern prokaryotes are involved in As detoxification of some form or kind³,⁷¹. These include specialized lineages involved in chemoautotrophic and photoautrophic As(III) oxidation, as well as respiratory As(V) reduction to As(III) species⁷¹.

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 seawater⁷⁹ and alkaline volcanic lakes⁸⁰. As a consequence, As detoxification is indispensable for microbial communities to succeed in these modern As-rich settings³,⁷¹,⁷⁸,⁸¹. Because of the ancient origin of biological As detoxification⁷⁶ 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 transporters⁷⁸. This interaction has been identified as a key cause of phosphate limitation in modern As(V)-rich alkaline volcanic lakes⁷⁸ and in seawater settings⁸². 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 molecule⁸³. In oxic bottom water environments, Fe (III) oxide-oxyhydroxides co-precipitate phosphate from the water column and sediment pore waters at the oxic/anoxic interface⁸⁴. 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 hematite³⁰,⁶⁹ in most of the epiclastic samples (Supplementary Table 1). A well-supported positive linear correlation between Fe₂O₃ and P₂O₅ in the epiclastic deposit (Supplementary Fig. 8), suggests P and As sequestration may be linked to Fe (III) oxide-hydroxide precipitation dynamics⁸⁴, 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 lakes⁷⁸ and hydrothermal ecosystems³.

## Conclusions

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.

## Methods

**Geological setting and sample preparation.** The Amane Tazgart succession crops out ~25 km southward of Ouarzazate city, and in the southern part of the Ediacaran Saghrro 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.

**Whole-rock geochemistry.** 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 (LiBO₂) and dissolved in dilute nitric acid (HNO₃). 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 tech-
niques described in Carignan et al.<sup>85</sup>. 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).

**Synchrotron radiation nano-XRF imaging.** 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, France<sup>86</sup>. 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.

**X-ray absorption near edge structure spectroscopy (nano-XANES).** Nano-XANES measurements were performed at the Nanoscopium beamline of Synchrotron Soleil. Two thin Si diodes were used to measure the I₀ 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-L₃ 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).

**Data availability**

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

Received: 5 February 2023; Accepted: 27 May 2023

Published online: 05 June 2023

**References**

1. 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).

2. 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).

3. 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).

4. 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).

5. 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).

6. 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).

7. 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).

8. 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).

9. 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).

10. 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).

11. 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).

12. 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).

13. 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).

14. 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).

15. 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).

16. 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).

17. 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).

18. 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).
19. 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).

20. 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).

21. 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).

22. 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).

23. 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).

24. 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).

25. Jewula, K., Środoń, J., Kuligiewicz, A., Mikolajczak, 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).

26. 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).

27. 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).

28. 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).

29. 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).

30. 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).

31. 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).

32. 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).

33. 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).

34. 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).

35. 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).

36. 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).

37. 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).

38. 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).

39. Shields, G. & Stille, P. Diagenetic constraints on the use of cerium anomalies as palaeosewater 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).

40. 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).

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

42. 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).

43. 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).

44. 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).

45. 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).

46. 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).

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

48. 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).

49. 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).

50. 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).

51. 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).

52. 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).

53. 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).

54. 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).

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

56. Pu, J. *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).
57. 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).

58. 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).

59. 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).

60. Adebayo, O. F., Akyinyemi, 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).

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

62. 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).

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

64. 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).

65. 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).

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

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

68. 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).

69. Chraiki, I., Bououagri, 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).

70. 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).

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

72. 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).

73. 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).

74. 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).

75. 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).

76. Chen, S. et al. Simultaneous removal of para-arsanilic acid and the released inorganic arsenic species by CuFe₂O₄ activated per-oxymonosulfate process. *Sci. Total Environ.* **742**, 140587. https://doi.org/10.1016/J.SCITOTENV.2020.140587 (2020).

77. 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).

78. 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).

79. 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).

80. 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).

81. 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).

82. 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).

83. 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).

84. 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).

85. Carignan, J., Hild, P., Mevelle, G., Morel, J. & Yeghichayan, 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).

86. 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).

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

88. 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).

89. McLennan, S. M. Relationships between the trace element composition of sedimentary rocks and upper continental crust. *Geochem. Geophys. Geosyst.* https://doi.org/10.1029/2000GC000109@10.1002/(ISSN)1525-2027.GERM1 (2001).

## Acknowledgements

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.
**Author contributions**

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.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-023-36031-5.

**Correspondence** and requests for materials should be addressed to A.E.A.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_12.png" /></figure>

**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/.

© The Author(s) 2023