# Microbially induced potassium enrichment in Paleoproterozoic shales and implications for reverse weathering on early Earth

Jérémie Aubineau, Abderrazak El Abani, Andrey Bekker, Andréa Somogyi, Olabode Bankole, Roberto Macchiarelli, Alain Meunier, Armelle Riboulleau, Jean-Yves Reynaud, Kurt O. Konhauser, et al.

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

## To cite this version:

Jérémie Aubineau, Abderrazak El Abani, Andrey Bekker, Andréa Somogyi, Olabode Bankole, et al. Microbially induced potassium enrichment in Paleoproterozoic shales and implications for reverse weathering on early Earth. Nature Communications, 2019, 10, pp.2370. 10.1038/s41467-019-10620-3. hal-02157763

HAL Id: hal-02157763

https://hal.science/hal-02157763

Submitted on 17 Jun 2019

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_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-019-10620-3

OPEN

# Microbially induced potassium enrichment in Paleoproterozoic shales and implications for reverse weathering on early Earth

Jérémie Aubineau<fnref n="1" />, Abderrazak El Abani<fnref n="1" />, Andrey Bekker², Andrea Somogyi ³, Olabode M. Bankole<fnref n="1" />, Roberto Macchiarelli⁴,⁵, Alain Meunier<fnref n="1" />, Armelle Riboulleau⁶, Jean-Yves Reynaud⁶ & Kurt O. Konhauser⁷

Illitisation requires potassium incorporation into a smectite precursor, a process akin to reverse weathering. However, it remains unclear whether microbes facilitate K⁺ uptake to the sediments and whether illitisation was important in the geological past. The 2.1 billion-year-old Francevillian Series of Gabon has been shown to host mat-related structures (MRS) and, in this regard, these rocks offer a unique opportunity to test whether ancient microbes induced illitisation. Here, we show high K content confined to illite particles that are abundant in the facies bearing MRS, but not in the host sandstone and black shale. This observation suggests that microbial biofilms trapped K⁺ from the seawater and released it into the pore-waters during respiration, resulting in illitisation. The K-rich illite developed exclusively in the fossilized MRS thus provides a new biosignature for metasediments derived from K-feldspar-depleted rocks that were abundant crustal components on ancient Earth.

<footnote n="1">¹UMR 7285 CNRS IC2MP, University of Poitiers, Poitiers 86073, France. ²Department of Earth and Planetary Sciences, University of California, Riverside, CA 92521, USA. ³Nanoscopium Beamline Synchrotron Soleil, BP 48, Saint-Aubin, Gif-sur-Yvette 91192, France. ⁴Department of Geosciences, University of Poitiers, Poitiers 86073, France. ⁵Department of Prehistory, UMR 7194 CNRS, National Museum of Natural History, Paris 75005, France. ⁶UMR 8187 CNRS LOG, University of Lille, ULCO, Villeneuve d'Ascq 59655, France. ⁷Department of Earth and Atmospheric Sciences, University of Alberta, Edmonton, AB T6G 2E3, Canada. Correspondence and requests for materials should be addressed to A.E.A. (email: abder.albani@univ-poitiers.fr)</footnote>
The conversion of smectite to illite–smectite mixed-layer minerals (I–S MLMs) is a multistep process that takes place during diagenesis, and is mainly controlled by temperature, time, and $K^+$ availability$^{1–4}$. The illitisation process progressively converts hydrated smectite-rich layers to illite-rich layers, resulting in a common co-occurrence of smectite-rich, randomly interstratified I–S MLMs and long-range ordered I–S MLMs$^{5,6}$. Two main mechanisms for smectite to illite transformation have been inferred, and these are likely linked to specific environmental conditions that are yet to be defined$^{3,5–8}$. A solid-state transformation characterized by a progressive replacement of smectite with illite via tetrahedral/octahedral Al substitutions and interlayer exchange has been proposed$^{3,7}$. The most likely mechanism involves the dissolution of the smectite lattice and growth of illite crystals due to the Ostwald ripening-like effect$^{5,6,8}$.

Microorganisms are capable of precipitating and transforming clay minerals as a result of their high surface reactivity and metabolic activity$^{9–17}$. Perhaps one of the most striking examples is the transformation of hydrated smectite to illite through a process that involves the microbial release of structural Fe(III) from the smectite lattice via dissimilatory iron reduction (DIR)$^{13–15}$. This results in the negative charge of the octahedral sheet and $K^+$ uptake into the interlayer spaces to balance the structural charge imposed by DIR, thus enhancing the illitisation process. In laboratory experiments, this clay transformation is enhanced by an abundant supply of $K^+$ cations to the medium$^{13,15}$. In natural systems, K-rich fluids and dissolution of K-feldspar may also provide $K^+$ for this reaction to take place$^{17}$. A biologically controlled $K^+$ addition to sediments during diagenesis has not yet been described.

The Francevillian Formation B (FB) marine sedimentary rocks are exquisitely preserved within the Francevillian basin in southeastern Gabon$^{18}$ (Fig. 1a; Supplementary Fig. 1a). The maximum burial depth is estimated to reach about 2 km with minimal burial temperature$^{19–21}$. An assemblage of diagenetic illite, chlorite, I–S MLMs, and K-bentonite characterizes the FB Formation$^{19,22,23}$. These sediments were deposited in an oxygenated environment$^{24}$ ~2.1 billion years (Ga) ago$^{25}$ during a

transgression (FB₁ Member) followed by a sea-level fall (FB₂ Member). The shallowing brought the depositional site to the photic zone and allowed the formation of mat-related structures (MRS) under the dominant influence of cyanobacteria$^{26,27}$. The marine FB₂ Member, which bears the MRS, is essentially a massive sandstone deposit (FB2a unit) overlain by black shales and diagenetic carbonates interbedded with thin siltstone layers (FB2b unit; Supplementary Fig. 1b). Here, we present whole-rock and in situ geochemical analyses, bulk XRD, electron microscopy images, and synchrotron-based elemental distribution maps from MRS and their host sediments. We studied samples from a 15-m-thick, coarse-grained sandstone that conformably underlies a 5-m-thick, thinly-bedded black shale that hosts the oldest large colonial macrofossils with evidence for organism motility previously described from the FB₂ Member of the FB Formation$^{28–30}$ (Supplementary Fig. 1b). In this study, we highlight that K enrichments are localized in illite particles within fossilized MRS, and, most crucially, that this provides insight for the search of biosignatures in early Earth's sediments where stable isotope evidence are controversial. We discuss how these findings have important implications for Earth's climate and ocean chemistry in the Paleoproterozoic.

## Results

**Textures of mat-related structures.** The investigated MRS exhibit a range of surface morphologies including those generated through mat propagation (mat-growth structures) and mat preservation/protection of structures formed independently from mat growth (mat-protected structures)$^{26,27,31}$. The elephant-skin texture (Fig. 1b) is one of the most common mat-growth morphologies showing evidence for biological activity (i.e., undirected gliding motility, phototactic behavior$^{32,33}$), while wrinkle and Kinneyia structures, linear ridges, and nodule-like structures point to the biostabilization of sediment. Mat laminae associated with irregular rounded pits (Fig. 1c) provide evidence for biofilms that were capable of bioweathering the substratum onto which they were attached$^{34}$. Scanning electron

<figure><img src="image_4.png" /><figcaption>**Fig. 1** Geological map and field photographs of the FB₂ Member outcrops, Gabon. **a** Geological map of the Francevillian basin adapted from Bouton et al.$^{18}$. The studied area is the Moulendé Quarry (green star). **b** Elephant-skin texture on the bedding plane of coarse-grained sandstones. Inset box shows reticulate patterns. **c** Micrometer-thick microbial mat laminae (blue arrow) on the bedding plane of sandstones with rounded pits (white arrows)</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Biogenetic fabrics in the mat-related structures. The MRS textures are shown through SEM images. **a** Tufted microbial fabrics developed above the poorly-sorted quartz sandstone. Yellow and blue arrows point to tufts and quartz grains, respectively. **b** A void-filling titanium oxides that may have filled an oxygen bubble produced within the microbial mat. **c** Nearly-circular void filled with titanium oxides at the tip of a cone-like feature (green arrows). Detrital dioctahedral micas (e.g., muscovite) are shown by purple arrows</figcaption></figure>

microscopy-back-scattered electron imaging (SEM-BSE) coupled to energy-dispersive X-ray spectrometry (EDX) displays textures and mineralogical compositions within the mat laminae that are different from host sediments (Fig. 2; Supplementary Figs. 2 and 3)²⁶.

Laminated black shales host pyritized MRS that are primarily composed of 20–30 µm euhedral pyrite grains (Supplementary Fig. 2c, d). The total sulfur content (higher than 6 wt%) within these MRS is associated with low total organic carbon content (Supplementary Table 1). In addition, laminations between pyrite crystals, reflecting successive mat growth, were recently documented²⁶. The non-pyritized MRS (with less than 4 wt% of S) comprise tufted microbial fabrics above coarse-grained sandstones (Fig. 2a) with abundant nearly circular voids within the mat structures filled with titanium oxides, and they are typically enriched in heavy minerals, such as titanium oxides and zircons (Fig. 2b, c; Supplementary Figs. 2e and 3; Supplementary Table 1). Inferred to be gas escape features, the voids show a precompactional formation whereby the clay particles are disorganized around the cone-like structures yet aligned above them.

**Potassium enrichment in the mat-related structures.** Whole-rock geochemical analyses of major elements are presented in Supplementary Table 1. Our data show higher K₂O/SiO₂ ratios for the MRS compared to the surrounding sediments, thus indicating a potassium-rich source specific to this lithology (Fig. 3). The observed difference in the dataset is statistically significant (n = 41; Kruskal–Wallis test: χ² = 33.29, df = 3, p-value < 10⁻⁷). Pairwise comparisons show that the difference between all pairs of groups is also statistically significant (Supplementary Table 2). The subtle, but significant difference in

ratios between pyritized mats and black shales from the FB2b unit might be related to contamination of some black shale sediments during the extraction of pyritized mat structures (Fig. 3). Binary plots of selected major elements were used to compare the MRS and host sediments (Supplementary Note 1; Supplementary Fig. 4), as well as the likely control of mineral composition and the effect of quartz dilution on these sediments.

Synchrotron-based scanning X-ray fluorescence microscopy (XFM) study was performed on millimeter-sized areas of representative samples. This non-destructive, element-specific imaging technique provides the distribution of S, K, Ca, Ba, Mn, Fe, Ni, Cu, Zn, Ga, Ge, and As with micron-scale spatial resolution and high analytical sensitivity. The distributions of S and K were imaged in pyritized MRS (Fig. 4a), black shale (Fig. 4b), and non-pyritized MRS (Fig. 4c). Sulfur is restricted to pyrite microcrystals (Fig. 4a, c), highlighting the difference between pyritized and non-pyritized MRS. Pyrite crystals are not enriched in K. Potassium is homogeneously distributed throughout the MRS and black shales, without significant micro-scale heterogeneities. The K-distribution maps reveal a higher K content in the pyritized and non-pyritized MRS compared to that in the host sediments.

**Potassium-rich clay assemblage patterns.** Eleven representative samples, consisting of five MRS and six host sediments from the FB₂ Member, were further analyzed using the X-ray diffraction (XRD) technique (Fig. 5a, d; Supplementary Figs. 5–7). The stacking mode (R, Reichweite ordering parameter) of I–S MLMs, ranging from R0 (randomly interstratified) to R ≥ 1 (long-range ordered) MLMs, was calculated (see Methods) to characterize the illite content with a higher *R* parameter corresponding to a higher
<figure><img src="image_6.png" /><figcaption>**Fig. 3** K₂O/SiO₂ ratios from mat-related structures and host sediments. Pyritized MRS are hosted by black shales and non-pyritized MRS are observed on bedding surface of sandstones. The data for each lithology are represented as box plots with a red square showing the mean, black diamonds corresponding to individual samples, 50% of the data are shown as a box and whiskers extend to 1.5 times the interquartile range</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 4** High-resolution sulfur and potassium distribution maps. **a** Pyritized MRS. Red arrows point to K enrichment between pyrite crystals and in the host sediment of MRS. **b** Black shales. **c** Non-pyritized MRS. Higher X-ray fluorescence (XRF) intensities correspond to higher S and K contents. For easier comparison, a common intensity scale was chosen for the K distribution maps</figcaption></figure>

illite content. The mineral assemblage of the MRS (both non-pyritized and pyritized types) includes quartz, mica, illite, chlorite, K-rich R3 I–S MLMs, and pyrite. In contrast, the host sediments contain quartz, mica, illite, chlorite, R0 I–S MLMs, and dolomite (Supplementary Figs. 5–7). Notably, K-feldspars have

been shown to be entirely absent, while plagioclase is present in the FB and FC formation lithologies, potentially reflecting pro-venance composition or depositional rather than climatic (via intense weathering) control¹⁹,³⁵. XRD patterns of the <2 µm fraction from the MRS show peaks at 10.56–10.37 Å in air-dried
<figure><img src="image_8.png" /><figcaption>**Fig. 5** Mineralogical composition of the <2 µm clay-size fractions. X-ray diffraction profiles of oriented preparations after air-drying (black lines) and glycolation (red lines) and their transmission-electron images are given. **a** Microbial mat specimen. **b** Large, well-crystallized lath from microbial mat laminae. The inset shows selected area electron diffraction (SAED) pattern. **hk0** pattern shows a hexagonal structure typical of phyllosilicates and coherently stacked crystals. **c** Hexagonal habit of illite from a mat sample. **d** Associated sandstone sediment. **e** Lathlike and poorly crystallized particles from the host sandstone. **f** Tiny hexagonal-shaped particles from the host sediment. [Green areas correspond to long-range ordered illite-smectite mixed-layer minerals (R3 I/S MLM); Blue areas represent randomly ordered illite-smectite mixed-layer minerals (RO I/S MLM); smectite (S); chlorite (Ch); mixed-layer (ML); illite/mica (I/M); anatase (An); quartz (Q); pyrite (Py); barite (Ba); calcite (Ca)]</figcaption></figure>

(AD) preparations splitting into 11.08–10.82 Å and ~9.9 Å peaks after ethylene glycolation (EG), consistent with characteristics of R3 I–S MLMs (Fig. 5a; Supplementary Fig. 6a, c). Smectite-rich R0 I–S MLMs, mixed with R3 I–S MLMs, in the host sediments are indicated by the shift of the 001 peak from around 14.35–14.1 Å in AD to 17.22–16.92 Å in EG (Fig. 5d; Supplementary Fig. 6b, d). The 002 peak of pure smectite close to 8.5 Å (Fig. 5d) further characterizes the host sediments.

The modeling results by the NEWMOD program³⁶ for the MRS confirm the lack of R0 I–S MLMs. In contrast, the relative abundance of R3 I–S MLMs, together with a relative enrichment in detrital mica and chlorite, are higher than in the host sediments (Supplementary Figs. 8 and 9; Supplementary Table 3), suggesting an intense dissolution of the reactive clay fraction and K input. In addition, smectite layers represent less than 12% of the R3 stacking in both the MRS and sediments, but are more abundant in the R0 I–S MLMs (S = 65%).

Transmission electron microscopy (TEM) was used to image the <2 µm clay-size fraction and, specifically, the illite crystal habit. The illites form large laths, up to 10 µm long and 0.6 µm wide (Fig. 5b; Supplementary Fig. 10a–c), as well as hexagonal particles ~0.2–0.5 µm in size in the MRS (Fig. 5c; Supplementary Fig. 10d–f). The large lath-shaped particles indicate a well-crystallized illite growth through an Ostwald ripening-like process in which the small particles are dissolved, while large particles grow at their expense in order to minimize surface free energy⁸,³⁷,³⁸. The hexagonal particles represent steady-state growth during an advanced stage of illitisation. The particle size of illite is much smaller in the host sediments than in the MRS.

Illite in the host sediments predominantly forms lath-shaped particles ~0.15–0.4 µm long and 0.05 µm wide, with a few hexagonal particles up to 0.2 µm in size (Fig. 5e, f; Supplementary Fig. 10g, h). These tiny lath-shaped particles correspond to poorly crystallized illite³⁸,³⁹.

## Discussion

The textural and geochemical observations of the pyritized MRS suggest that the organic carbon was oxidized via dissimilatory sulfate-reduction (DSR)²⁹. This mechanism is also supported by the common diagenetic mineralization of individual mat layers, arguing for the former presence of organic carbon⁴⁰. Moreover, if the earlier interpretation of the MRS in the FB₂ Member being cyanobacterial in origin is correct²⁶, then the tufted appearance might be related to the escape of oxygen bubbles similar to that described in modern and ancient stromatolites⁴¹⁻⁴³. The differences between non-pyritized MRS and pyritized MRS might reflect a higher sedimentation rate and energy in the depositional setting of the non-pyritized MRS where a rapid burial would not have allowed sulphidic pore-waters to develop. Alternatively, this difference might be due to less reactive iron initially supplied to the depositional setting where the non-pyritized MRS were formed, consistent with the inferred shallow-water and high energy conditions.

The illitisation process results in different abundances of illite and smectite, but points toward the association of the illite end-member with the MRS. The source of K⁺ for diagenetic smectite-to-illite transformation is commonly linked to the dissolution of K-rich tuffs, K-feldspars, and micaceous minerals, as well as to
migrating K-rich fluids³⁹. Petrographic and XRD observations argue against significant contribution of K-rich tuffs to the composition of the FB₂ Member. The survival of R0 I–S MLMs (i.e., smectite-rich) in both the sandstones and black shales lacking the MRS is consistent with incomplete illitisation, attributed to the lack of K-feldspars in the FB Formation sediments¹⁹,³⁵. In our samples, dioctahedral micas (e.g., muscovite) are only mechanically deformed (Fig. 2c; Supplementary Fig. 3), reflecting a limitation to no chemical alteration under diagenetic conditions⁴⁴ and, therefore, they cannot be considered as a significant source of K⁺. Moreover, it is unlikely that an external abiotic supply (e.g., K-rich fluids or K-bearing minerals) would have provided K⁺ only to the µm-thick microbial fabrics, but not to the host sediments. Organic matter could delay illitisation by inhibiting the exchange of interlayer cations⁴⁵; nonetheless non-pyritized MRS show both high organic carbon content (Supplementary Table 1) and an advanced stage of illitisation.

Our observations imply instead that microbes played two roles in the illitisation process. First, the large mineral surface area of smectite may have helped to bind organic compounds within interlayer spaces⁴⁶. The concentration of organic matter is critical because illitisation is promoted via DIR, a process that oxidizes the sorbed organic carbon via reduction of Fe(III) minerals, and ultimately releases Fe(II) to the sediment pore-waters¹⁵. Second, the microbes sequestered K⁺ into the MRS of the FB₂ Member. In the case of the latter, potassium is one of the most essential cations for microorganisms to maintain their cellular machinery⁴⁷, and, as a consequence, it is typically enriched in cells relative to seawater⁴⁸.

Biologically derived K could have been incorporated into MRS through the K⁺ uptake into the living microbial cells, and/or K⁺ absorption by reactive organic macromolecules—containing carboxyl, phosphate, and hydroxyl groups—as part of the cell envelope or extracellular polymeric substances (EPS) of biofilms. In the first process, the structure of the K⁺ channel (e.g., proteins controlling K⁺ transport across cell membranes) forms a selectivity filter that prevents Na⁺ uptake, but works perfectly for incorporating K cations⁴⁹. Thus, K⁺ contributes to the regulation of osmotic pressure differences, the control of internal pH, and the activation of intracellular enzymes⁴⁷; the thylakoid potassium channel (SynK) even promotes an optimal photosynthesis in cyanobacteria⁵⁰. As a result of K-dependant metabolic reactions, the intracellular concentration of K⁺ in modern cells (100 mM) is one order of magnitude greater than that in the seawater⁴⁸,⁵¹, further supporting the ability of microbes to sequester this macronutrient. The second K⁺ enrichment process involves the electrostatic binding of metal cations from solution to the cell wall and/or EPS, and in turn, the nucleation and mineral precipitation of authigenic mineral phases⁵²,⁵³. Indeed, microbial cell surfaces have been shown to concentrate K⁺ from solution⁵⁴, and importantly, increased K⁺ supply has been shown to stimulate the production of EPS by some bacterial strains⁵⁵. In addition, microbial respiration of EPS during early diagenesis could liberate adsorbed cations into the pore waters, thus inducing K⁺ super-saturation that can additionally promote mineral authigenesis⁵⁶. Regardless of the path for K⁺ uptake, biologically derived K⁺ would have been required to facilitate the smectite-to-illite transformation during both decay of microbial mat organic matter and early diagenesis.

The proposed mechanism of bacterially enhanced illitisation has important implications for the Paleoproterozoic climate and ocean chemistry because the facilitated removal of K⁺ from seawater into the clay fraction of sediments is akin to the modern reverse weathering process⁵⁷ to explain difference in the chemical composition of rivers and seawater. Formation of authigenic clays

in the oceans involves (as an unbalanced reaction): Si(OH)₄ + cations (K⁺, Mg²⁺, Fe²⁺, Al³⁺) + HCO₃⁻ → clay mineral + CO₂ + H₂O. In the Amazon continental shelf, reverse weathering was estimated to account for a sink of ~10% of the annual riv-erine supply of K⁺ to the oceans⁵⁸. Moreover, the release of CO₂ from bicarbonate consumption has been recently inferred to have helped maintain a warm Precambrian Earth at a time of lower solar luminosity and stabilized marine pH at lower values by releasing protons to seawater⁵⁹. While the conversion of smectite to illite is different⁶⁰: smectite + K⁺ → illite + cations (Na⁺, Ca²⁺, Mg²⁺, Fe²⁺) + Si(OH)₄, its importance lies in the enfor-cing ultimate removal of K⁺ from seawater with MRS, which would have strengthened the reverse weathering buffer on global climate during the Precambrian. Notably, in the absence of ter-restrial vegetation and limited clay burial on landmasses in the Precambrian to consume K⁺ released by chemical weathering, terrestrial K⁺ fluxes to the oceans would have been significantly larger than in the Phanerozoic⁶¹. The K content in Precambrian seawater is highly uncertain with estimates ranging from those lower than that in the modern seawater to those that are higher⁶²⁻⁶⁵. Large land masses likely did not emerge until ~2.5 Ga⁶⁶, thus the continental margin repository for marine eva-porites was smaller with the correspondingly larger content of Na and K in seawater⁶²,⁶⁷,⁶⁸. Regardless, the ability of bacteria to facilitate the incorporation of K⁺ into illite at the time when microbial mats were likely more extensive in the oceans in the absence of grazing metazoans⁶⁹ provides a previously unrecog-nized biological feedback on the Earth surface system to control seawater pH and the atmospheric CO₂ content.

The recognition of ancient metabolic activity more specific than biogenic graphite is also of paramount importance in the search for Archean biosignatures. In this regard, we suggest that the process of illitisation in association with the MRS might be evidence of biologically supplied K⁺. Until now, this concept has not been tested because of the common abundance of K-feldspars in Phanerozoic and Precambrian sedimentary basins, which is related to the presence of felsic igneous rocks in the continental crust of all ages⁷⁰. However, the predominance of mafic and ultramafic igneous rocks in the Archean and early Proterozoic greenstone belts⁷¹,⁷² requires a different mechanism to explain the K⁺ enrichment with respect to host sediments that lack potassium-rich detrital minerals and illite-rich authigenic clays. Our finding of K-controlled mineral transformation during diagenesis in clay-rich shale deposits provides a new insight for seeking a potential biosignature in the early Earth's sedimentary environments that lacked granitoids in their provenance and where morphologic and carbon isotope evidence for MRS are often compromised.

## Methods

**Sampling and sample preparation.** We collected more than 100 MRS with their associated host sediments from Moulendé Quarry, Francevillian basin, Gabon. The samples were photographed at the University of Poitiers using a Nikon Europe D610 digital single-lens reflex camera equipped with a Nikon AF-S 24–120 mm f/4G ED VR lens. Textural contrast between mat specimens and host sediments was investigated on polished slab sections with a ZEISS Discovery V8 stereoscope coupled with an Axio Cam ERc 5s microscope camera. Thirty samples covering all rock types were selected for mineralogical analyses. The µm-thick mat laminae (both non-pyritized and pyritized, defined according to petrographic analyses²⁶ and their sulfur content) were carefully removed from their underlying host rocks by using a cutter blade to avoid contamination. Directly below the mat-bearing rocks, sandstone or black shale sediments were gently extracted with a hammer for mineralogical and geochemical analyses. MRS were carefully removed to extract and prepare <2 µm clay fraction and whole-rock powder for mineral and organic carbon analyses.

**X-ray diffraction.** Approximately the same quantity of the whole-rock powder and clay mineral fraction (<2 µm) from 11 samples were analyzed with a Bruker D8 ADVANCE diffractometer at the University of Poitiers using CuKα radiation
operating at 40 kV and 40 mA. The <2 µm clay fraction was separated by disper-
sion of gently hand-crushed bulk samples in deionized water with an Elma S60
ultrasonic agitation device without any chemical pre-treatment<sup>73</sup>. The dispersed
particles were let to settle under gravity at a controlled room temperature of 20 °C,
and centrifuged to separate <2 µm clay fraction. Oriented slides were prepared by
drying ~1 mL of suspension on glass slides at room temperature. The <2 µm
fraction of some smectite-rich and illite-rich samples has been Ca-saturated after
three dispersions in 1 M CaCl₂ solution followed by several washings with deio-
nized water. This cation exchange homogenizes the sheets of swelling clays to keep
stable interlayer properties. Analysis of whole-rock powder samples was performed
over an angular range of 2–65° 2θ with a step size of 0.025° 2θ per 3 s. Oriented
slides of the <2 µm clay fraction were analyzed at a step size of 0.02° 2θ per 3 s
counting time and over 2–30° 2θ angular range after successive air drying (AD) and
ethylene glycol (EG) saturation. Oriented Ca-saturated samples were then prepared
and analyzed in AD and EG states. Background stripping, indexing of peaks, and
mineral identification were done using Bruker Eva software by comparing with
International Centre for Diffraction Data (ICDD) files. Illite polytypes were
observed on randomly oriented powders (<2 µm clay fraction) and analyzed at a
step size of 0.01° 2θ per 2 s counting time and over 19–33° 2θ angular range. The
results were compared with reference data<sup>74–76</sup>. Illite with 2M₁ and 1Mᵢ polytypes
represent detrital and diagenetic illite, respectively<sup>77</sup>.

**X-ray diffraction profile modeling.** The stacking mode (R, Reichweite ordering parameter) of I–S MLMs ranging from randomly interstratified (R = 0) MLMs to long-range ordered MLMs (R = 1, 2, and ≥3) as well as the proportion of smectite and illite layers were determined by fitting experimental Ca-saturated samples in AD and EG states (over 2–30° 2θ CuKα) using the NEWMOD 2.0 program simulation³⁶. When the same solution is obtained for the AD and EG states, we assumed the fit is acceptable. The instrumental and experimental factors such as horizontal and vertical beam divergences, goniometer radius, length, and thickness of the oriented slides were introduced without further adjustment⁷⁸. Sigmaster (degree of orientation of the clay particles on the slide) was set between 12° and 18° and the mass absorption coefficient (μ*) to 45 cm²/g, as suggested by Moore and Reynolds⁷³. Angular domains containing non-clay mineral reflections were not included in the profile-fitting process.

**Electron microscopy.** Polished slab sections of selected samples were carbon coated and imaged with a FEI Quanta 200 scanning electron microscope (SEM) equipped with an energy dispersive X-ray spectrometer (EDX) at the University of Lille. Mineral identification and documentation of textural relationships were acquired in back-scattered electron mode (BSE) operated at accelerating voltage of 15 kV, 1 nA beam current, and a working distance of 10.5 mm. TEM analysis of the <2 µm clay fraction was performed with a JEOL 2100 UHR microscope working at 200 kV accelerating voltage with a LaB6 emitter and equipped with a JEOL EDX detector at the University of Poitiers. Samples were prepared for TEM analysis by suspending ~1 mg of each sample in distilled water. A droplet of diluted suspension after ultrasonic treatment was deposited on a C-coated Au grid. Lath- and hexagon-shaped particles were selected from the grids because they only correspond to illitic minerals³⁷⁻³⁹. Selected-area electron diffraction (SAED) patterns were collected to image the hexagonal net pattern of the analyzed particles in (hk0) reflections⁷⁹⁻⁸¹.

**Whole rock analysis.** Whole-rock geochemical analyses of major elements were carried out on 18 biofilms and 23 host sediments at Service d'Analyse des Roches et Minéraux (SARM) of the Centre de Recherches Pétrographiques et Géochimiques (CRPG), Nancy, France. Each sample was powdered in agate mortar and approximately 1 g was fused with lithium metaborate (LiBO₂) and dissolved in nitric acid. Major element concentrations were obtained by inductively-coupled plasma atomic emission spectrometry (ICP-AES).

**Statistics.** For descriptive treatment (boxplot) and statistical analyses, we used the R version 3.5.1<sup>82</sup>. The values for each group do not fulfill the assumption of normality and the Bartlett test of homogeneity of variances was computed prior statistical analyses. Our data returned a significant result (Bartlett test, $\chi^2 = 26.92$, df = 3, *p*-value < 10<sup>-6</sup>), indicating overall inter-group heteroscedasticity of the dataset. According to these results, a multiple factor non-parametric test (i.e., Kruskal-Wallis test) was used to calculate whether the difference in K₂O/SiO₂ ratios is statistically significant or not among lithologies. Finally, post hoc pairwise multiple comparison tests according to Conover were applied to figure out which pairs of groups (i.e., lithologies) are different.

**Carbon analysis.** MRS and their host sediments were powdered for the mea-
surement of organic matter content. We used a FlashEA 1112 (ThermoFisher
Scientific) CHNS analyzer for flash dynamic combustion at 970 °C under a
constant helium flow. An Eager 300 software was used for data acquisition. A
calibration curve was performed with aspartic acid and nicotinamide before each
analysis. Carbon content includes both inorganic and organic carbon, but whole-
rock geochemical data show that CaO content is low.

**Synchrotron-based scanning X-ray microscopy.** Polished slab sections of representative samples were chosen for scanning XRM analysis. The XRM measurements were performed on several areas within the studied samples at the Nanoscopium hard X-ray nanoprobe beamline<sup>83</sup> of Synchrotron Soleil (L'Orme des Merisiers Saint-Aubin, France). The monochromatic X-ray beam of 12 keV energy was focused on the sample area by a Kirkpatrick-Baez nano-focusing mirror. In order to obtain micrometer resolution for elemental maps of mm²-sized sample areas, we used the fast continuous scanning (FLYSCAN) technique<sup>84</sup>. Full XRF spectra were collected for each pixel of the scans by two Si-drift detectors (VITUS H50, KETEK GmbH) in order to increase the solid angle of detection. The XRF spectra of the two detectors were added and this sum was used for calculating the distributions of S, K, Ca, Ba, Mn, Fe, Ni, Cu, Zn, Ga, Ge, and As in the measured sample regions. Each elemental map was normalized to 10 ms/pixel exposure time.

**Data availability**

The authors declare that the data supporting this study are available within the paper and
its Supplementary Information files.

Received: 27 September 2018 Accepted: 20 May 2019

Published online: 17 June 2019

**References**

1. Velde, B., Suzuki, T. & Nicot, E. Pressure-temperature-composition of illite/ smectite mixed-layer minerals: Niger delta mudstones and other examples. *Clays Clay Miner.* **34**, 435–441 (1986).

2. Perry, E. A. Jr. & Hower, J. Burial diagenesis in Gulf Coast pelitic sediments. *Clays Clay Miner.* **18**, 165–177 (1970).

3. Hower, J., Eslinger, E. V., Hower, M. E. & Perry, E. A. Jr. Mechanism of burial metamorphism of argillaceous sediment: 1. Mineralogical and chemical evidence. *Geol. Soc. Am. Bull.* **87**, 725 (1976).

4. Srodon, J. Mixed-layer illite-smectite in low-temperature diagenesis: data from the Miocene of the Carpathian Foredeep. *Clay Miner.* **19**, 205–215 (1984).

5. Ferrage, E., Vidal, O., Mosser-Ruck, R., Cathelineau, M. & Cuadros, J. A reinvestigation of smectite illitization in experimental hydrothermal conditions: results from X-ray diffraction and transmission electron microscopy. *Am. Mineral.* **96**, 207–223 (2011).

6. Larson, B., Sakharov, B. A., Claret, F. & Drits, V. A. Diagenetic smectite-to-illite transition in clay-rich sediments: a reappraisal of X-ray diffraction results using the multi-specimen method. *Am. J. Sci.* **309**, 476–516 (2009).

7. Drits, V. A., Sakharov, B. A., Lindgreen, H. & Salyn, A. Sequential structure transformation of illite-smectite-vermiculite during diagenesis of Upper Jurassic shales from the North Sea and Denmark. *Clay Miner.* **32**, 351–371 (1997).

8. Nadeau, P. H., Wilson, M. J., McHardy, W. J. & Tait, J. M. The conversion of smectite to illite during diagenesis: evidence from some illitic clays from bentonites and sandstones. *Mineral. Mag.* **49**, 393–400 (1985).

9. Fortin, D., Ferris, F. G. & Scott, S. D. Formation of Fe-silicates and Fe-oxides on bacterial surfaces in samples collected near hydrothermal vents on the Southern Explorer Ridge in the northeast Pacific Ocean. *Am. Mineral.* **83**, 1399–1408 (1998).

10. Wacey, D. et al. Enhanced cellular preservation by clay minerals in 1 billion-year-old lakes. *Sci. Rep.* **4**, 5841 (2014).

11. Pace, A. et al. Microbial and diagenetic steps leading to the mineralisation of Great Salt Lake microbialites. *Sci. Rep.* **6**, 31495 (2016).

12. Konhauser, K. O. & Urrutia, M. M. Bacterial clay authigenesis: a common biogeochemical process. *Chem. Geol.* **161**, 399–413 (1999).

13. Kim, J., Dong, H., Seabaugh, J., Newell, S. W. & Eberl, D. D. Role of microbes in the smectite-to-illite reaction. *Science* **303**, 830–832 (2004).

14. Vorhies, J. S. & Gaines, R. R. Microbial dissolution of clay minerals as a source of iron and silica in marine sediments. *Nat. Geosci.* **2**, 221–225 (2009).

15. Zhang, G., Kim, J., Dong, H. & Sommer, A. J. Microbial effects in promoting the smectite to illite reaction: role of organic matter intercalated in the interlayer. *Am. Mineral.* **92**, 1401–1410 (2007).

16. O'Reilly, S. E., Watkins, J. & Furukawa, Y. Secondary mineral formation associated with respiration of nontronite, NAu-1 by iron reducing bacteria. *Geochem. Trans.* **6**, 67–78 (2005).

17. Cuadros, J. Clay minerals interaction with microorganisms: a review. *Clay Miner.* **52**, 235–261 (2017).

18. Bouton, P. et al. *Carte géologique de la République du Gabon à 1/200 000, Feuille Franceville - Boumango* (DGMG—Ministère des Mines, du Pétrole, des Hydrocarbures, 2009).

19. Ngombi-Pemba, L., El Albani, A., Meunier, A., Grauby, O. & Gauthier-Lafaye, F. From detrital heritage to diagenetic transformations, the message of clay minerals contained within shales of the Palaeoproterozoic Francevillian basin (Gabon). *Precambrian Res.* **255**, 63–76 (2014).
20. Gauthier-Lafaye, F. & Weber, F. The Francevillian (Lower Proterozoic) uranium ore deposits of Gabon. *Econ. Geol.* **84**, 2267–2285 (1989).

21. Gauthier-Lafaye, F. & Weber, F. Natural nuclear fission reactors: time constraints for occurrence, and their relation to uranium and manganese deposits and to the evolution of the atmosphere. *Precambrian Res.* **120**, 81–100 (2003).

22. Ossa Ossa, F. et al. Exceptional preservation of expandable clay minerals in the ca. 2.1 Ga black shales of the Francevillian basin, Gabon and its implication for atmospheric oxygen accumulation. *Chem. Geol.* **362**, 181–192 (2013).

23. Bankole, O. M. et al. Earth's oldest preserved K-bentonites in the ca. 2.1 Ga Francevillian Basin, Gabon. *Am. J. Sci.* **318**, 409–434 (2018).

24. Canfield, D. E. et al. Oxygen dynamics in the aftermath of the Great Oxidation of Earth's atmosphere. *Proc. Natl Acad. Sci. USA* **110**, 16736–16741 (2013).

25. Horie, K., Hidaka, H. & Gauthier-Lafaye, F. U-Pb geochronology and geochemistry of zircon from the Franceville series at Bidoudouma, Gabon. in *Proceedings of the 15th Annual Goldschmidt Conference* (2005).

26. Aubineau, J. et al. Unusual microbial mat-related structural diversity 2.1 billion years ago and implications for the Francevillian biota. *Geobiology* **16**, 476–497 (2018).

27. Reynaud, J.-Y. et al. Depositional setting of the 2.1 Ga Francevillian macrobiota (Gabon): rapid mud settling in a shallow basin swept by high-density sand flows. *Sedimentology* **65**, 670–701 (2017).

28. El Albani, A. et al. Large colonial organisms with coordinated growth in oxygenated environments 2.1 Gyr ago. *Nature* **466**, 100–104 (2010).

29. El Albani, A. et al. The 2.1 Ga old Francevillian biota: biogenicity, taphonomy and biodiversity. *PLoS ONE* **9**, e99438 (2014).

30. El Albani, A. et al. Organism motility in an oxygenated shallow-marine environment 2.1 billion years ago. *Proc. Natl Acad. Sci. USA* **116**, 3431–3436 (2019).

31. Sarkar, S., Bose, P., Samanta, P., Sengupta, P. & Eriksson, P. Microbial mat mediated structures in the Ediacaran Sonia Sandstone, Rajasthan, India, and their implications for Proterozoic sedimentation. *Precambrian Res.* **162**, 248–263 (2008).

32. Shepard, R. N. & Sumner, D. Y. Undirected motility of filamentous cyanobacteria produces reticulate mats: motility produces reticulate mats. *Geobiology* **8**, 179–190 (2010).

33. Reyes, K. et al. Surface orientation affects the direction of cone growth by *Leptolyngbya* sp. strain C1, a likely architect of coniform structures Octopus Spring (Yellowstone National Park). *Appl. Environ. Microbiol.* **79**, 1302–1308 (2013).

34. Flemming, H.-C. et al. Biofilms: an emergent form of bacterial life. *Nat. Rev. Microbiol.* **14**, 563–575 (2016).

35. Bankole, O. M., El Albani, A., Meunier, A. & Gauthier-Lafaye, F. Textural and paleo-fluid flow control on diagenesis in the Paleoproterozoic Franceville Basin, South Eastern, Gabon. *Precambrian Res.* **268**, 115–134 (2015).

36. Reynolds, R. C. Jr. & Reynolds III, R. C. NEWMOD for Windows. *The Calculation of One-dimensional X-ray Diffraction Patterns of Mixed-layered Clay Minerals* (1996).

37. Inoue, A. & Kitagawa, R. Morphological characteristics of illitic clay minerals from a hydrothermal system. *Am. Mineral.* **79**, 700–711 (1994).

38. Inoue, A., Velde, B., Meunier, A. & Touchard, G. Mechanism of illite formation during smectite-to-illite conversion in a hydrothermal system. *Am. Mineral.* **73**, 1325–1334 (1988).

39. Meunier, A. & Velde, B. *Illite. Origins, Evolution and Metamorphism* (Springer-Verlag, Berlin, 2004).

40. Schieber, J. Microbial mats in terrigenous clastics: the challenge of identification in the rock record. *Palaios* **14**, 3–12 (1999).

41. Bosak, T., Liang, B., Sim, M. S. & Petroff, A. P. Morphological record of oxygenic photosynthesis in conical stromatolites. *Proc. Natl Acad. Sci. USA* **106**, 10939–10943 (2009).

42. Bosak, T. et al. Formation and stability of oxygen-rich bubbles that shape photosynthetic mats. *Geobiology* **8**, 45–55 (2010).

43. Sallstedt, T., Bengtson, S., Broman, C., Crill, P. M. & Canfield, D. E. Evidence of oxygenic phototrophy in ancient phosphatic stromatolites from the Paleoproterozoic Vindhyan and Aravalli Supergroups, India. *Geobiology* **16**, 139–159 (2018).

44. Bassett, W. A. Role of hydroxyl orientation in mica alteration. *Geol. Soc. Am. Bull.* **71**, 449–456 (1960).

45. Li, Y., Cai, J., Song, M., Ji, J. & Bao, Y. Influence of organic matter on smectite illitization: a comparison between red and dark mudstones from the Dongying Depression, China. *Am. Mineral.* **101**, 134–145 (2016).

46. Kennedy, M. J., Pevear, D. R. & Hill, R. J. Mineral surface control of organic carbon in black shale. *Science* **295**, 657–660 (2002).

47. Epstein, W. The roles and regulation of potassium in bacteria. *Prog. Nucleic Acid Res. Mol. Biol.* **75**, 293–320 (2003).

48. Mulkidjanian, A. Y., Bychkov, A. Y., Dibrova, D. V., Galperin, M. Y. & Koonin, E. V. Origin of first cells at terrestrial, anoxic geothermal fields. *Proc. Natl Acad. Sci. USA* **109**, E821–E830 (2012).

49. Doyle, D. A. The structure of the potassium channel: molecular basis of K⁺ conduction and selectivity. *Science* **280**, 69–77 (1998).

50. Checchetto, V. et al. Thylakoid potassium channel is required for efficient photosynthesis in cyanobacteria. *Proc. Natl Acad. Sci. USA* **109**, 11043–11048 (2012).

51. Williams, R. J. P. & Da Silva, J. J. R. F. *The Chemistry of Evolution: The Development of Our Ecosystem* (Elsevier, Amsterdam, 2006).

52. Dupraz, C. & Visscher, P. T. Microbial lithification in marine stromatolites and hypersaline mats. *Trends Microbiol.* **13**, 429–438 (2005).

53. Konhauser, K. O. Diversity of bacterial iron mineralization. *Earth-Sci. Rev.* **43**, 91–121 (1998).

54. Beveridge, T. J. & Fyfe, W. S. Metal fixation by bacterial cell walls. *Can. J. Earth Sci.* **22**, 1893–1898 (1985).

55. Wilkinson, J. F., Duguid, J. P. & Edmunds, P. N. The distribution of polysaccharide production in *Aerobacter* and *Escherichia* strains and its relation to antigenic character. *J. Gen. Microbiol.* **11**, 59–72 (1954).

56. Dupraz, C. et al. Processes of carbonate precipitation in modern microbial mats. *Earth-Sci. Rev.* **96**, 141–162 (2009).

57. Mackenzie, F. T. & Garrels, R. M. Chemical mass balance between rivers and oceans. *Am. J. Sci.* **264**, 507–525 (1966).

58. Michalopoulos, P. & Aller, R. C. Rapid clay mineral formation in Amazon delta sediments: reverse weathering and oceanic elemental cycles. *Science* **270**, 614–617 (1995).

59. Isson, T. T. & Planavsky, N. J. Reverse weathering as a long-term stabilizer of marine pH and planetary climate. *Nature* **560**, 471–475 (2018).

60. Boles, J. R. & Franks, S. G. Clay diagenesis in Wilcox sandstones of southwest Texas: implications of smectite diagenesis on sandstone cementation. *J. Sediment. Petrol.* **49**, 55–70 (1979).

61. McMahon, W. J. & Davies, N. S. Evolution of alluvial mudrock forced by early land plants. *Science* **359**, 1022–1024 (2018).

62. Holland, H. D. *The Chemical Evolution of the Atmosphere and Oceans* (Princeton University Press, Princeton, 1984).

63. Weiershäuser, L. & Spooner, E. T. C. Seafloor hydrothermal fluids, Ben Nevis area, Abitibi Greenstone Belt: Implications for Archean (~2.7 Ga) seawater properties. *Precambrian Res.* **138**, 89–123 (2005).

64. Spear, N. et al. Analyses of fluid inclusions in Neoproterozoic marine halite provide oldest measurement of seawater chemistry. *Geology* **42**, 103–106 (2014).

65. Marty, B., Avice, G., Bekaert, D. V. & Broadley, M. W. Salinity of the Archaean oceans from analysis of fluid inclusions in quartz. *C. R. Geosci.* **350**, 154–163 (2018).

66. Bindeman, I. N. et al. Rapid emergence of subaerial landmasses and onset of a modern hydrologic cycle 2.5 billion years ago. *Nature* **557**, 545–548 (2018).

67. Knauth, L. P. Salinity history of the Earth's early ocean. *Nature* **395**, 554–555 (1998).

68. Knauth, L. P. Temperature and salinity history of the Precambrian ocean: implications for the course of microbial evolution. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **219**, 53–69 (2005).

69. Seilacher, A. & Pflüger, F. From biomats to benthic agriculture: a biohistoric revolution. In *Biostabilization of Sediments* (eds Krumbein, W. E. et al.) 97–105 (Bibliotheks und Informationssystem der Carl von Ossietzky Universität Oldenberg: (BIS)-Verlag, 1994).

70. Greber, N. D. et al. Titanium isotopic evidence for felsic crust and plate tectonics 3.5 billion years ago. *Science* **357**, 1271–1274 (2017).

71. Arndt, N., Barnes, S. J. & Lesher, C. M. *Komatiite* (Cambridge University Press, Cambridge, 2008).

72. Moyen, J.-F. & Martin, H. Forty years of TTG research. *Lithos* **148**, 312–336 (2012).

73. Moore, D. M. & Reynolds, R. C. Jr. *X-ray Diffraction and the Identification and Analysis of Clay Minerals* (Oxford University Press, New York, 1997).

74. Bailey, S. W. Structure of layer silicates. in *Crystal Structure of Clay Minerals and Their X-ray Identification* (eds. Brindley, G. W. & Brown, G.) 1–123 (Mineralogical Society, London, 1980).

75. Drits, V. A. X-ray identification of one-layer illite varieties: application to the study of illites around uranium deposits of Canada. *Clays Clay Miner.* **41**, 389–398 (1993).

76. Haines, S. H. & van der Pluijm, B. A. Clay quantification and Ar–Ar dating of synthetic and natural gouge: application to the Miocene Sierra Mazatán detachment fault, Sonora, Mexico. *J. Struct. Geol.* **30**, 525–538 (2008).

77. Grathoff, G. H., Moore, D. M., Hay, R. L. & Wemmer, K. Origin of illite in the lower Paleozoic of the Illinois basin: evidence for brine migrations. *Geol. Soc. Am. Bull.* **113**, 1092–1104 (2001).

78. Drits, V. A. & Tchoubar, C. X-ray Diffraction by Disordered Lamellar Structures: *Theory and Applications to Microdivided Silicates and Carbons* (Springer-Verlag, Berlin, 1990).

79. Whitney, G. & Velde, B. Changes in particle morphology during illitization: an experimental study. *Clays Clay Miner.* **41**, 209–218 (1993).

80. Kumai, M. Identification of nuclei and concentrations of chemical species in snow crystals sampled at the South Pole. *J. Atmos. Sci.* **33**, 833–841 (1976).
81. Veblen, D. R., Gutherie, G. D. Jr., Livi, K. J. T. & Reynolds, R. C. Jr. High-resolution transmission electron microscopy and electron diffraction of mixed-layer illite/smectite: experimental results. *Clays Clay Miner.* **38**, 1–13 (1990).

82. R. Core Team. *A Language and Environment for Statistical Computing* (2018).

83. Somogyi, A. et al. Optical design and multi-length-scale scanning spectroscopy possibilities at the Nanoscopium beamline of Synchrotron Soleil. *J. Synchrotron Radiat.* **22**, 1118–1129 (2015).

84. Medjoubi, K. et al. Development of fast, simultaneous and multi-technique scanning hard X-ray microscopy at Synchrotron Soleil. *J. Synchrotron Radiat.* **20**, 293–299 (2013).

## Acknowledgements

The authors are grateful to the Gabonese Government, CENAREST, General Direction of Mines and Geology, and Agence Nationale des Parcs Nationaux de Gabon for logistic support. This work was supported by CNRS, the University of Poitiers, Nouvelle Aquitaine Region, and the French Embassy Libreville, Gabon. The authors would like to thank Prof. P. Mouguiama Daouda for his support. The authors also acknowledge E. Chi Fru, B. Dazas, F. Hubert, K. Medjoubi, and L. Pallas for scientific discussions and C. Fontaine, C. Laforest, S. Pronier, and P. Recourt for laboratory support at the Universities of Poitiers and Lille.

## Author contributions

A.E.A. conceived the project. J.A., A.E.A., A.B., and K.O.K. designed research. J.A., A.E.A., O.M.B., and J.-Y.R. performed field research. J.A. performed sample preparation and XRD analyses. A.S. carried out the synchrotron experiments and processed the data. A.M., A.R., R.M., J.-Y.R., A.B., and K.O.K. provided critical input to the manuscript. J.A., A.E.A., A.B., and K.O.K. wrote the manuscript with input from all co-authors.

## Additional information

**Supplementary Information** accompanies this paper at https://doi.org/10.1038/s41467-019-10620-3.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Peer review information:** *Nature Communications* thanks Warren Huff and other anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

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

<figure><img src="image_9.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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019