Abstract
Microbialites–layered, organosedimentary deposits–exist in the geologic record and extend back in deep time, including all estimated times of inner core nucleation. Microbialites may preserve magnetic field variations at high-resolution based on their estimated growth rates. Previous studies have shown that microbialites can have a stable magnetization. However, the timing and origin of microbialite magnetization were not well determined, and no study has attempted to evaluate whether actively growing microbialites record the geomagnetic field. Here, we present centimeter-scale magnetization and magnetic property variations within the structure of modern microbialites from Great Salt Lake (GSL), United States, and Laguna Bacalar, Mexico, Pleistocene microbialites from GSL, and a Cambrian microbialite from Mongolia. All samples record field directions close to the expected value. The dominant magnetic carrier has a coercivity of 35–50 mT and unblocking temperatures are consistent with magnetite. A small proportion of additional high coercivity minerals such as hematite are also present, but do not appear to appreciably contribute to the natural remanent magnetization (NRM). Magnetization is broadly consistent along microbialite layers, and directional variations correlate with the internal slope of the layers. These observations suggest that the documented NRM may be primarily detrital in origin and that the timing of magnetization acquisition can be close to that of sediment deposition.
Introduction
Recent work demonstrates the potential of speleothems as paleomagnetic recorders. These laminated cave deposits form by precipitation of carbonate and include detrital or authigenic magnetic particles that appear to document paleosecular variation (PSV) at high temporal resolution (e.g., Lascu and Feinberg, 2011; ; ; Trindade et al., 2018). Stromatolites are also finely laminated carbonate materials, thought to have often formed by a microbially-mediated balance between sedimentation and intermittent lithification (Reid et al., 2000). Compared to speleothems, stromatolites are well represented in the geologic record, and the oldest examples go back at least 3.5 billion years (e.g., ; Riding, 2010; ). This spans all estimated times of inner core nucleation, which range from the late Archean to early Proterozoic (2.4–2.5 Ga) (e.g., Tarduno et al., 2006; Valet et al., 2014), to Mesoproterozoic times (∼1.2 Ga) (e.g., ; Smirnov et al., 2016), to the Ediacaran period (∼570 Ma) (e.g., ). These materials may, therefore, hold promise as paleomagnetic recorders with high temporal resolution, and this type of record in deep time could have potential for capturing geomagnetic field behavior when the core may have been in different states.
Modern stromatolites have widely varying estimated growth rates ranging from 0.05 mm/yr to about 50 mm/yr (Paull et al., 1992; ; ; Petryshyn et al., 2012). Laminations (when present), represent hiatal periods in sediment accretion, are thought to form approximately every 4–10 years in some studied stromatolites (Paull et al., 1992; Petryshyn et al., 2012), with one study finding yearly laminations (Rasmussen et al., 1993). In modern stromatolites (Reid et al., 2000), microbial filaments trap and bind sand and carbonate grains during periods of rapid accretion. During hiatal periods, a continuous exopolymer surface film develops, and bacterial activity promotes formation of a thin micritic crust which may be further modified into a thicker, fused carbonate crust. These laminated crusts therefore represent a sequence of lithified and preserved microbial mat surfaces. Laminations in ancient (>1 billion year old) stromatolites have been interpreted to represent day-night growth cycles (e.g., Vanyo and Awramik, 1982), although these interpretations have been increasingly called into question in light of the work on modern stromatolites mentioned above.
In this work, we look more broadly at microbialites, an umbrella term that includes at least some stromatolites. define microbialites as organosedimentary deposits accreted by a benthic microbial community via trapping and binding detrital sediment and/or precipitation of calcite (abiotic or biologically induced). Microbialites have a relatively wide environmental distribution, found at the sediment-water interface in marine, alkaline freshwater, hypersaline lakes, evaporitic, and hydrothermal environments. They can be classified by different internal structures as stromatolite (laminated structure), thrombolite (clotted structure lacking lamination), dendrolite (dendritic structure), leolite (structureless), or some intergradational combination (Riding, 2000; Riding, 2011). Stromatolite diversity abundance increased from the late Archean (∼2.4 Ga), peaking around ∼1.25 Ga (). Diversity declined up to the Eradican (∼550 Ma) (Walter and Heys, 1985; ), but resurged in the early Phanerozoic and late Cambrian when thrombolites and dendrolites became the predominant type of microbialite (Riding, 2000; Kiessling, 2002). Grazing and metazoan competition have been suggested as a major factor contributing to the decline of stromatolite diversity, but other factors such as seawater chemistry and atmospheric composition could also affect their diversity (Riding, 2006). There are some arguments that not all Archean stromatolites are biogenic in origin (Lowe, 1994; ; McLoughlin et al., 2008) but this interpretation remains controversial (; ; ; Schopf et al., 2007).
A primary magnetization in microbialites may arise from a detrital remanent magnetization (DRM) or a biogenic magnetization, or some combination of the two. Microbialites may trap detrital ferromagnetic (sensu lactu) particles from nearby sources within the microbialite matrix and carry a DRM. Lund et al. (2010) examined coral reef cores and observed detrital titanomagnetite bound within microbialites that had filled the primary cavities of the coral framework. These grains were determined to carry the natural remanent magnetization (NRM) and appeared to accurately record geomagnetic field variability. A DRM may also be biomediated in that the sticky microbial mat may play a role in particle retention and immobilization. Petryshyn et al. (2016) used bulk magnetic susceptibility to differentiate between stromatolites that are biogenic or abiogenic in origin and demonstrated that a cyanobacterial mat traps and binds grains more readily than abiogenic carbonate precipitates.
The other potential source of magnetic particles is biogenic magnetite which can be produced by biologically controlled mineralization (BCM) and biologically induced mineralization (BIM). Magnetotactic bacteria can create single-domain magnetic particles in chains by a matrix mediated BCM process. The presence of biogenic magnetite has been inferred in ancient () and modern (Stolz et al., 1989) microbialite environments by using rock magnetic techniques and transmission electron microscopy (TEM) to identify single domain magnetite with the unique crystal structures typically associated with magnetotactic bacteria. However, it has since been demonstrated that similar crystal morphologies can also be produced abiotically through reductive alteration of nanogoethite (Till et al., 2017) or via thermal decomposition of siderite (). Maher and Taylor (1988) also observed similar fine-grained magnetite in soils and suggested that it may have formed via oxidation of Fe2+ solutions (Taylor et al., 1987). Biogenic magnetite might also be produced extracellularly by iron-reducing bacteria through BIM processes (Lovley, 1991; Konhauser et al., 2011; Konhauser and Riding 2012). Maloof et al. (2007) used TEM and ferromagnetic resonance to infer the presence of superparamagnetic magnetite formed via extracellular processes. Although superparamagnetic material will not contribute to remanence, extracellular magnetite typically has a wide grain size distribution including grains that are large enough to have stable single domain behavior (e.g., ).
Relatively little paleomagnetic work has intentionally focused on microbialites. Vanyo and Awramik (1982) extracted a paleomagnetic direction via alternating field demagnetization from a single specimen of late Proterozoic stromatolite. Thermal demagnetization of Jurassic stromatolites from Poland have shown them to carry a stable magnetization, including a polarity reversal (Muraszko, 2014) also observed in middle-late Jurassic limestones, stromatolites, and radiolarites in Poland (Kądziałko-Hofmokl and Kruczyk, 1987). As noted above, microbialites found in a post Last-Glacial-Maximum and late Pleistocene coral reef framework appear to accurately record the field direction (Lund et al., 2010). Several studies have sampled in or near ancient carbonate beds containing microbialites. The Love’s Creek Member of the Bitter Springs Formation contains stromatolite and microbialite facies and was interpreted to be remagnetized through magnetite authigenesis along with the remainder of the section (Swanson-Hysell et al., 2012). However, the lower Ordovician Oneota dolomite which is locally stromatolitic carries a direction interpreted to be potentially primary (), and Neoproterozoic dolomites containing microbialite laminites are interpreted to carry primary DRM (Trindade et al., 2003; ). Maloof et al. (2007) sampled modern platform carbonate muds in the Bahamas and proposed that early cemented, laminated microbially-bound carbonates may preserve syn-depositional paleomagnetic directions. Few of these studies focus on the timing of magnetization with respect to microbialite formation, and no studies have attempted to evaluate whether actively growing microbialites accurately record the known field.
This study is designed to further assess the feasibility of microbialites as paleomagnetic recorders. We examined two actively growing or recently lithified microbialites, as well as one Cambrian microbialite. Variations in magnetization and magnetic mineralogy were evaluated at the centimeter scale. We also consider the timing of the magnetization with respect to deposition or microbialite growth.
Methods
Sampling and Settings
This study focuses on microbialites from three localities/environments: 1) Great Salt Lake (GSL), Utah, United States (hypersaline lake environment); 2) Laguna Bacalar (LB), Quintana Roo, Mexico (freshwater alkaline environment); 3) the Cambrian Bayan Gol (BG) formation, southwestern Mongolia (marine environment) (Figure 1). Detailed information about locations, sample types, and descriptions are provided in Table 1.
FIGURE 1
TABLE 1
| Location/ID | Latitude | Longitude | Features | Status/Age |
|---|---|---|---|---|
| GSL18/Site 1 (GSL18_01) | 41.06054 | −112.24909 | Ring structure microbialite | Living and recently lithified |
| GSL18/Site 2 (GSL18_02) | 41.05709 | −112.25351 | Ring structure microbialite | Living and recently lithified |
| GSL18/Site 3 (GSL18_03) | 41.06014 | −112.26168 | Dome style thrombolite | Recently lithified |
| GSL18/Site 4 (GSL18_04) | 41.02802 | −112.27246 | Ridge structure and ring structure microbialite | Living and recently lithified |
| GSL19/Site 1 (GSL19_01) | 41.21402 | −112.85211 | Monk-head structure microbialite | Pleistocene |
| GSL19/Site 2 (GSL19_02) | 41.15964 | −112.86058 | Cauliflower structure microbialite | Pleistocene |
| LB19/Site 1 (LB19_01) | 18.64766 | −88.39658 | Thrombolitic stromatolites with cyanobacteria layer | Living and unlithified sediments |
| LB19/Site 2 (LB19_02) | 18.61383 | −88.42595 | Poorly lithified porous microbialite | Holocene or recently lithified? |
| LB19/Site 3 (LB19_03) | 18.55834 | −88.43537 | Poorly lithified porous microbialite with cyanobacteria layer | Living and recently lithified |
| BG14/Site 1 (BG14_01) | 46.4213 | 96.18608 | Dome style laminated stromatolite | Cambrian |
Microbialite sample locations and descriptions.
At the Great Salt Lake, living and recently lithified modern microbialites were sampled from shallow water on the shores of Antelope Island in the southwestern part of the lake (GSL18). Samples were taken from living, low-profile ring-shaped structures and elongated, narrow ridge structures (; ), but we were unable to maintain orientation, and the soft samples easily fell apart. These samples were preserved in lake water during transportation back to the laboratory. One oriented hand sample was collected from the steep side of a lithified, dome-style thrombolite structure that has been exposed by a lake level drop over the past few decades. Pleistocene microbialites are exposed along the western shoreline of the lake at Lakeside (Newell et al., 2017), and oriented hand samples were collected from two locations (GSL19).
Living deposits and lithified giant Holocene microbialites (; ) were collected from shallow water at three sites at Laguna Bacalar. Microbialites from Site 1 are unlithified layered domal structure deposits with a stromatolitic structure at the bottom but transitioning to a more thrombolitic fabric toward the top (Figure 2C). The approximate thickness of the actively growing greenish microbial mat was ∼2 cm from the surface. Site 2 microbialites were lithified but porous. The Site 3 giant microbialites had soft surfaces including living cyanobacteria layers (∼2 cm) and lithified microbialites underneath. Living microbial mats on the surface at Site 1 (LB19_0101) and 3 (LB19_0304) were easily penetrated and collected by pushing standard 2 cm3 plastic cubes) into the mats. One hand sample from Site 1 (LB19_0102) was subsampled with 2 cm3 plastic cubes. Site 1 microbialites were additionally sampled with 5 cm diameter by 14 cm long (LB19_0103) and 7 cm diameter by 30 cm long (LB19_0104) plastic cores by push coring. Lithified samples were collected at Sites 2 and 3 using an electric drill with a 1 inch diamond coring bit. Soft samples collected by push core and directly in cubes were frozen after sampling.
FIGURE 2
A lower Cambrian stromatolite limestone head was taken from the middle of Unit 18 (Khomentovsky and Gibsher, 1996) of the Bayan Gol section within the Bayan Gol Formation. The sample was only oriented with respect to stratigraphic up, and with relatively high uncertainty as the sample was not originally collected for paleomagnetic study. The sample has parallel, but not horizontal, internal laminations (Figure 2F), consistent with convex-upward growth over a domal structure.
In the laboratory, lithified hand samples were sliced perpendicular to laminations or growth directions. Each slice was further subdivided along a grid forming samples of 1–1.5 cm3. These specimens, as well as the unoriented living GSL specimens, were immobilized in 2 cm3 plastic cubes using nonmagnetic paddings. Push cores from LB were split longitudinally while frozen and were then allowed to soften so they could be subsampled vertically every 3.5 cm using 2 cm3 cubes. Four specimens from each vertical horizon were collected. Electric drilled samples were cut into 2 cm long specimens. Types of samples and the number of sub-sampled specimens are listed in Table 2.
TABLE 2
| Sample Id | Type | O | X (n) | X. σ | NRM (n) | NRM. σ | ARM (n) | ARM. σ | IRM (n) | IRM. σ | S-ratio (n) | S. σ |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| GSL18_0101 | HS | U | 2.05E-07 (2) | 3.86E-07 | 4.79E-07 (2) | 3.36E-07 | 9.42E-07 (2) | 6.57E-07 | 6.01E-05 (2) | 3.63E-06 | ||
| GSL18_0102 | HS | U | 1.67E-07 (1) | 2.52E-07 (1) | 4.38E-07 (1) | |||||||
| GSL18_0104 | HS | U | −7.63E-08 (1) | 2.50E-07 (1) | 2.42E-07 (1) | |||||||
| GSL18_0201 | HS | U | 5.97E-07 (1) | 4.84E-07 (1) | 7.42E-07 (1) | |||||||
| GSL18_0202 | HS | U | 9.63E-07 (4) | 6.66E-07 | 3.34E-06 (4) | 4.60E-06 | 8.15E-07 (4) | 1.93E-07 | ||||
| GSL18_0301 | HS | F | 9.90E-07 (55) | 5.07E-07 | 4.42E-07 (55) | 4.42E-07 | 6.58E-07 (55) | 2.92E-07 | 9.05E-05 (19) | 5.87E-05 | 0.96 (19) | 0.03 |
| GSL18_0401 | HS | U | 9.80E-08 (2) | 1.34E-07 | 1.47E-06 (2) | 7.10E-07 | 7.94E-07 (2) | 8.35E-08 | ||||
| GSL18_0402 | HS | U | 9.91E-07 (4) | 3.81E-07 | 7.07E-07 (4) | 5.98E-07 | 9.88E-07 (4) | 3.72E-07 | ||||
| GSL19_0103 | HS | F | 3.12E-07 (45) | 1.03E-06 | 9.76E-07 (45) | 8.96E-07 | 0.92 (35) | 0.01 | ||||
| GSL19_0201 | HS | F | 8.26E-07 (42) | 6.16E-07 | 5.08E-07 (42) | 4.64E-07 | 0.92 (42) | 0.02 | ||||
| GSL19_0202 | HS | F | 9.05E-07 (22) | 4.57E-07 | 4.69E-07 (22) | 3.41E-07 | 1.01E-06 (22) | 3.46E-07 | 7.65E-05 (6) | 2.72E-05 | 0.92 (16) | 0.02 |
| LB19_0101 | BS | F | −5.99E-07 (8) | 2.56E-08 | 2.08E-08 (8) | 9.14E-09 | 8.92E-08 (8) | 8.58E-09 | ||||
| LB19_0102 | BS | F | −5.66E-07 (9) | 2.09E-08 | 4.49E-08 (9) | 3.68E-08 | 1.02E-07 (2) | 2.18E-08 | 5.17E-06 (1) | |||
| LB19_0103 | PS | P | −4.71E-09 (8) | 4.76E-10 | 3.84E-08 (8) | 2.61E-08 | ||||||
| LB19_0104 | PS | P | −5.47E-09 (36) | 7.53E-09 | 2.17E-08 (36) | 2.69E-08 | ||||||
| LB19_0201 | DC | P | 4.87E-08 (3) | 1.62E-08 | ||||||||
| LB19_0301 | DC | P | 5.12E-08 (5) | 7.24E-09 | 1.03E-07 (5) | 3.80E-08 | 3.46E-06 (5) | 3.01E-06 | 1.02 (5) | 0.02 | ||
| LB19_0302 | DC | P | 1.34E-07 (6) | 1.16E-07 | ||||||||
| LB19_0304 | BS | F | −5.90E-07 (6) | 6.48E-08 | 1.87E-07 (6) | 4.13E-07 | 1.08E-07 (2) | 8.07E-09 | 7.14E-06 (2) | 3.34E-06 | ||
| BG14_0101 | HS | P | 3.49E-09 (71) | 1.98E-09 | 3.33E-07 (71) | 5.54E-07 | 9.53E-07 (28) | 4.73E-07 | 2.79E-05 (35) | 1.55E-05 | 0.76 (35) | 0.14 |
Sample types and magnetic properties.
O, orientation, n, number of measured specimens, , mean magnetic susceptibility (m3/kg),. σ, standard deviation of magnetic susceptibility (m3/kg), NRM, mean natural remanent magnetization intensity (Am2/kg), NRM. σ, standard deviation of NRM (Am2/kg), ARM, mean anhysteretic remanent magnetization (Am2/kg), ARM. σ, standard deviation of ARM (Am2/kg), IRM, mean 1000 mT isothermal remanent magnetization (Am2/kg), IRM. σ, standard deviation of IRM (Am2/kg), S-ratio, mean value of -IRM300mT/IRM1000mT, S. σ, standard deviation of S-ratio. HS, hand sample, BS, block sample directly collected from surface, PC, push piston core samples, DC, drilled core samples by portable electric drill, U, unoriented, F, fully oriented, P, partially oriented.
Magnetic Methods
NRM and magnetic susceptibility () were measured on all specimens with the exception of on LB drill core specimens. Most specimens were stepwise alternating field (AF) demagnetized to the point where a specimen lost 85–95% of the initial magnetization intensity. LB push core specimens from samples LB_0103, 0104 were not demagnetized because they dried and orientation was lost following NRM measurement. Frozen living LB samples LB_0101 and 0102 were demagnetized, but specimens thawed during measurement, leading to movement inside the cubes and noisy data that could not be interpreted. Specimen orientation during demagnetization was alternated between + X/+Y/+Z and -X/-Y/-Z directions to reduce the effects of possible unwanted bias field inside the AF demagnetizing device. Additional specimens (GSL19_0103 and BG14_0101) were subjected to stepwise thermal demagnetization to 700°C for comparison. A best-fit characteristic remanent magnetization (ChRM) was calculated using principal component analysis (Kirschvink, 1980) from the higher coercivity or temperature fraction. To calculate sample-average ChRM statistics (), we accepted AF-demagnetized specimens with anchored maximum angular deviation (MAD) < 10°. Data interpretation and sample-level Fisher mean calculations were performed using the PmagPy software (Tauxe et al., 2016).
Following AF demagnetization, a single-step anhysteretic remanent magnetization (ARM) was acquired in a 150 mT AF with 0.05 mT bias field. An isothermal remanent magnetization (IRM) was acquired in steps between 0 mT and 1,000 mT on fully demagnetized specimens. Some specimens had already been exposed to high fields, and in this case a stepwise IRM was acquired starting from a backfield of -1,000 mT IRM. These two methods are symbolized as IRM0mT for the first case and IRM-1000mT for the second case. To calculate the S-ratio, fully AF-demagnetized (to 200 mT) specimens were given an IRM at 300 mT after an exposure to a saturation field at a 1,000 mT in the opposite direction. The S-ratio was calculated based on the formula of Thompson and Oldfieid, 1986 as S-ratio = −IRM300mT/IRM1000mT.
The first derivatives of IRM acquisition curves were unmixed mathmatically to provide information on coercivity distributions and magnetic mineral populations (e.g., Robertson and France, 1994). This was carried out using MaxUnmix (Maxbauer et al., 2016) online software which fits a set of skew normal distributions to the data. For each distribution, the mean coercivity (MC), dispersion parameter (DP), relative proportion, skewness, and standard deviations were obtained.
Finally, to better constrain magnetic mineralogy, selected specimens were subjected to thermal demagnetization of a 3-axis IRM (Lowrie, 1990). IRMs of 1,000, 300, and 100 mT were applied in three orthogonal directions to separate the hard, medium, and soft coercivity fractions, and specimens were then thermally demagnetized to 700°C.
Remanence measurements were made using a 2G Enterprises 755SRMS Superconducting Rock Magnetometer inside a magnetically shielded room. Magnetic susceptibilities were acquired using an AGICO MFK1-FA Multifunction Kappabridge susceptibility bridge. AF demagnetization and ARM acquisition were undertaken with an ASC D-2000 AF demagnetizer. Stepwise thermal demagnetization was undertaken by an ASC Thermal Demagnetizer, and IRM was imparted with an ASC Impulse Magnetizer. These measurements were all conducted at the Paleomagnetism Laboratory at the University of Wisconsin- Milwaukee. Additionally, magnetization vs. applied field (up to 1 T) measurements were conducted on a Princeton Measurements vibrating sample magnetometer at the Institute for Rock Magnetism, University of Minnesota. The samples were all weakly magnetic and had a high diamagnetic:ferrimagnetic or paramagnetic:ferrimagnetic ratio, and hysteresis parameters could not be accurately calculated. We instead report the high-field behavior of the samples.
Results
Magnetization Intensity and Susceptibility Variations
With the exception of LB samples, all samples possess an average NRM intensity 1.00E-07 Am2/kg. Table 2 displays average magnetic susceptibility, NRM, ARM and IRM1000mT with standard deviations of each sample. The NRM intensity measurements for all samples (>1.00E-11 Am2) were not limited by the sensitivity of the cryogenic SQUID magnetometer (∼1.00E-12 Am2 moment sensitivity limit). However, specimens with NRM intensity less than 1.00E-10 Am2 approach the magnetization of sample container (1.00E-13 to 1.00E-11 Am2). measurements are limited by the sensitivity of the MFK1 Kappabridge susceptibility bridge (∼2.00E-08 SI or ∼2.00E-13 m3).
To assess internal spatial variations in magnetization and distribution of magnetic minerals, mass-normalized magnetic susceptibility and NRM intensity of four GSL hand samples, two LB push core samples, and one BG hand sample were contoured by cubic interpolation (Figure 2). Spatial variations in magnetization within these microbialite structures were compared to their internal laminations with the aim of understanding magnetization processes and assessing NRM reliability.
Samples GSL19_0103, GSL19_0201, and GSL19_0202 have a spatial correlation between NRM intensity and magnetic susceptibility. GSL19 microbialite NRM and susceptibility vary along the growth directions and are consistent along internal laminations (Figures 2A,B, Supplementary Figure S1). With the exception of a high NRM intensity spot at the bottom right, LB19_0104 has roughly consistent NRM intensity distribution along internal structures (Figure 2D). The magnetic susceptibility mapping shows high values around 18 cm depth, the transition level where stromatolitic laminations changed to a thrombolytic structure (Figure 2E). BG14_0101 has roughly constant NRM intensity (Figure 2G) and magnetic susceptibilities (Figure 2H) along laminations, with a low value in middle layers compared to outer and bottom layers. Magnetic susceptibility and NRM of GSL18_0301 and LB19_0103 (not shown), which have a width less than 10 cm, have similar magnetic susceptibility and NRM intensity trends. However, no obvious correlations with laminations are observed within these small samples.
Magnetic Mineralogy
Results from all samples are consistent with magnetite as a dominant magnetic carrier. This component has a mean coercivity of ∼29–57 mT (Supplementary Table S1), consistent with single domain (SD) to ‘pseudo-single-domain' (PSD) minerals that have non-uniform spin behavior. Modern samples (e.g., GSL, LB) are dominated by this component and have large negative high-field diamagnetic slopes. The Cambrian BG sample, however, has a mixture of low- and high-coercivity minerals with positive high-field paramagnetic slopes (Supplementary Figure S2).
GSL18: The S-ratio range of 19 lithified GSL18_0301 specimens is 0.92–1.06, with an average of 0.96, indicating that low coercivity ferrimagnetic minerals such as magnetite are dominant. Based on IRM unmixing results, the dominant magnetic components of the GSL18 living and GSL18_0301 lithified specimens have a mean coercivity of ∼56 mT for four IRM0mT unmixing results and ∼38 mT for four IRM-1000mT unmixing results, with similar DPs of 2.7 mT. Coercivity spectra derived from IRM acquisition data for all specimens are shown in Figure 3. There is little variability between living and lithified specimens from the same site. The three-component IRM demagnetization experiment shows a gradual decrease of the soft (71–78%) and medium (17–23%) coercivity fractions to 575°C and flattening after 620°C (Figure 4A). This gradual unblocking is consistent unblocking of domain walls, but may also represent some compositional variability. A high-coercivity fraction contributes <7% to the total IRM, and measurement noise made it difficult to evaluate a maximum unblocking temperature.
FIGURE 3
FIGURE 4

Three-component IRM technique thermal demagnetization results for microbialite samples (A) GSL18 (B) GSL19 (C) LB19, and (D) BG14. Three-component IRM (Am2/kg) produced with applied field of 100 mT (soft), 300 mT (intermediate), and 1,000 mT (hard). Soft, medium, and hard fractions are shown as blue circle, orange triangle, and green square, respectively. Additional plots are included in Supplementary Figure S4.
GSL19: The average S-ratio value of 93 GSL19 specimens is slightly lower than GSL18, with an average of 0.92 and a range of 0.87–0.96 (Table 2). All three GSL19 samples have a spatially correlated NRM intensity and S-ratio. High NRM correlates with high S-ratio (Supplementary Figures S1, S3). Only a single magnetic component can be confidently identified in six specimens of the Pleistocene GSL19_0202 based on IRM-1000mT unmixing. This dominant component has a mean coercivity of 38 mT with a DP of 3.1 mT, again consistent with SD to ‘PSD’-like magnetite. The soft (68–73%) and medium (22–26%) coercivity fractions flatten out at ∼580°C in the three-component IRM test (Figure 4B). Although a higher-coercivity component was not identified in the relatively noisy IRM unmixing data, hard minerals comprise a small fraction (5–8%) of the three-component IRM and gradually decrease to 675°C, consistent with hematite.
LB19: LB results are similar to those observed in GSL18 specimens. The S-ratio range of five LB19_0301 specimens is 0.99–1.05 with a mean value of 1.02, again consistent with low-coercivity minerals like magnetite. LB19_0102, 0301, and 0304 were subjected to IRM unmixing experiments and all LB microbialites have a consistent dominant component with a mean coercivity of 28–50 mT with a DP of 2.2–2.9. No major differences are observed between living microbial sediments and lithified microbialite samples in terms of component numbers, major MC, or DP. In the three-component IRM experiment, the low-coercivity fraction accounted for the highest proportion (75–81%) of IRM for 4 specimens of LB19_0301. This component gradually decreased and flattened out after 575°C–600°C, indicating magnetite (Figure 4C). With the exception of LB19_0301_C, the medium coercivity (17–22%) components have maximum unblocking temperatures of 550°C, and there is no significant hard coercivity fraction (<5%) (Supplementary Figure S4).
BG14: Unlike the GSL and LB samples, the BG stromatolite has variable magnetic mineralogy with a large S-ratio range from 0.42 to 0.95 (Figure 5). Low coercivity minerals are present in higher concentrations in inner layers and high coercivity minerals occur in outermost layers except two points at the left top. This is also shown in IRM acquisition (Figure 3C) and unmixing curves (Figures 3D, 5) where outer layers have a significant high coercivity component (purple lines in Figure 5). Based on the three-component IRM results (Figure 4D), the soft fraction is the largest (51–69%) and gradually unblocks up to ∼575°C, again consistent with magnetite. The medium fraction (20–25%) decreases to ∼300–350°C for low S-ratio specimens (e.g., BG14_0101_A05 and 12), consistent with titanomagnetite and/or pyrrhotite. The medium fraction in high S-ratio specimens additionally decreases to ∼650°C (e.g., BG14_0101_A12 and 70), consistent with (titano)hematite and/or thermally stable maghemite (Supplementary Figure S4). The hard fraction of low S-ratio and high S-ratio specimens accounts for 20–30% and <12% of the IRM, respectively. This fraction decreases to about 625°C–675°C, consistent with the presence of hematite.
FIGURE 5

IRM unmixing curve and S-ratio distribution for sample BG14_0101. IRM results were plotted through the MaxUnmix online software (Maxbauer et al., 2016) on a log scale. Gray dots represent actual data, and yellow lines represent total IRM best-fit based on real data. Colored lines (blue and purple) represent different coercivity components. The shadings represent 95% confidence intervals calculated using a resampling algorithm with N = 300. An asterisk next to a specimen name indicates that the specimen went through the -1000mT IRM process.
Directional Analysis
GSL: Fifteen measured (unoriented) actively growing cyanobacteria layers are fully demagnetized by around 40–100 mT, and a stable paleomagnetic direction is isolated up to 60 mT (Figure 6A). A small overprint is typically removed by ∼7.5 mT. Of the lithified (oriented) samples GSL18_0301, GSL19_0103, and GSL19_0202, 47% (26/55), 100% (10/10) and 35% (8/23) respectively have stable AF demagnetization behavior (Figure 6C), and 80% (8/10) of GSL19_0103 specimens have stable thermal demagnetization behavior (Figure 6D). Sample-level Fisher mean declinations (Figure 7A) of all samples are close to that predicted for a geocentric axial dipole (GAD) or to models of paleosecular variation for the past 10 kyr (CALS10k.2;
FIGURE 6

Vector endpoint diagrams for representative microbialite specimens. (A) AF demagnetization pattern for living Antelope Island GSL microbialite specimen GSL18_0202_01 (B) AF demagnetization pattern for living lithified LB Site 3 microbialite specimen LB19_0301_A (C) AF demagnetization pattern for specimen Lakeside lithified GSL microbialite GSL19_0103_15 (D) thermal demagnetization of specimen GSL19_0103_04 (E) AF demagnetization for BG stromatolite specimen BG14_0101_48; and (F) thermal demagnetization for specimen BG14_0101_06.
FIGURE 7

Equal area stereographic projection of magnetization components observed from (A) GSL (B) LB, and (C) BG samples. The mean declination and inclination with α95 ellipses of NRM or ChRM. Where declination was not available (LB19_0104) mean inclination is shown as a circle around the center with α95. Open and closed symbols represent directions in the upper and lower hemisphere, respectively. A gray star in (A) GSL (B) LB panel represents GAD value of each site. The black line on GSL and LB stereographic projections represent the CALSk10.2 model based on the latitude and longitude of GSL sampling site from 10 ka to present (
LB: Living sediments (e.g., LB19_0101, 0102, 0304) and lithified microbialites (e.g., LB19_0301) were fully demagnetized by 25–60 mT and 60–70 mT, respectively. Five specimens of lithified microbialite LB19_0301 have stable AF-demagnetization, and magnetic directions are isolated between 12.5 mT and 60 mT (Figure 6B). Most of the LB sites were located in moving water, which may complicate interpretation of DRM directions. LB19 Site 1 was in a stagnant environment, so only specimens from this site were used for directional analysis. These samples were not demagnetized (see Magnetic Methods), we report only NRM directions for this site. The
TABLE 3
| D | I | N* | N | R | k | α95 | |
|---|---|---|---|---|---|---|---|
| GSL18_0301 (ChRM) | 356.8 | 43.9 | 26 | 55 | 24.5 | 16.9 | 7.1 |
| GSL19_0103 (ChRM) | 3.76 | 36.7 | 18 | 20 | 17.6 | 44.8 | 5.2 |
| GSL19_0202 (ChRM) | 350.2 | 62.9 | 8 | 23 | 7.9 | 67.0 | 6.8 |
| LB19_0101 (NRM) | 346.6 | 43.6 | 8 | 6 | 3.5 | 34.9 | |
| LB19_0102 (NRM) | 40.4 | 36.6 | 8 | 6.5 | 4.5 | 29.3 | |
| LB19_0104 (NRM) | 36.8 | 36 | 31.5 | 7.8 | 7.15 | ||
| BG14_0101 (ChRM) | 18.5 | −66.3 | 22 | 37 | 21.7 | 70.2 | 3.7 |
Sample-average Fisher statistics.
D, mean declination; I, mean inclination; N*, number of MAD <10 specimens; N, number of total specimens; R, the length of the resultant vector; k, the precision parameter; α95; the circle of 95% confidence angle about the mean.
BG: BG14_0,101 was only oriented with respect to up, but specimens in the structure were mutually oriented. AF demagnetization was more effective in extracting stable magnetizations (59.4%), and an overprint is typically removed by ∼30 mT. None of the thermally processed specimens (N = 10) had a coherent demagnetization pattern (Figures 6E,F). Pilot thermomagnetic data, c(T), show the production of magnetite at temperatures > ∼400°C (Supplementary Figure S5), consistent with the thermal decomposition of siderite (e.g., Pan et al., 2000). ChRM Fisher mean inclinations for 22 out of 37 AF-demagnetized specimens (−66.3° ± 3.7°) overlap with the inclination (−62 ± 4°) obtained from siliclastic sediments with calcareous interbeds in the Bayan Gol formation (
FIGURE 8

(A) Inclination mapping of the BG14_0101 sample slice with contouring by cubic interpolation. The black arrow and dashed lines represent up directions and laminations. (B) Equal area stereographic projection of reoriented magnetization components of the sample BG14_0101 by rotating total mean declination to the magnetic north. Specimen declinations and inclinations from the left and right sides with respect to the stratigraphic up direction (black arrow) are symbolized as red and blue circles, respectively. The mean direction of each group is marked as a square symbol with α95 ellipses. The gray shading annulus on the BG panel (inclination: −62 ± 4°) represents the expected directions derived from sedimentary beds in the same formation (
Discussion
Magnetic Minerology
While S-ratio values for all modern microbialites (i.e., GSL18, LB19) are high (i.e., magnetite), ancient samples (i.e., GSL19, BG14) have lower values that indicate the presence of additional high coercivity materials. IRM unmixing results and three-component IRM tests show that all modern samples are dominated by a component with a mean coercivity of 37 mT for IRM-1000mT and 48 mT for IRM0mT, with the DP of 2.7, consistent with magnetite. The BG microbialite has a high relative proportion of high coercivity and high unblocking-temperature (up to 675°C) minerals consistent with hematite. A signature from high coercivity fractions is observed in IRM unmixing results for weakly magnetic living GSL18 and LB19 cyanobacteria layer samples. This might represent a small proportion of high coercivity material, but experimental noise means this component is poorly constrained. While three-component IRM results contain a hint of hard fractions (<7%), we infer that these signatures are experimental noise based on S-ratio results.
This difference in magnetic mineralogy between young and old microbialite samples may be attributed to different depositional material in the original environment. It is also possible that the high-coercivity minerals formed during diagenesis, leading to concerns of a partial or complete overprint in the ancient samples. However, these high-coercivity minerals do not appear to contribute to the NRM, which is fully AF-demagnetized at < 200 mT and is held entirely by low-coercivity minerals.
The coercivity distribution for nearly all GSL and LB specimens is remarkably similar, although GSL specimens have slightly higher coercivity (Figure 3). These distributions are similar to IRM acquisition for some remagnetized Devonian carbonates (
In Figure 3B, we also compare our coercivity spectra derived from the IRM0mT with spectra derived from AF demagnetization of a 300 mT (
The BG samples have a similar component consistent with detrital and/or extracellular magnetite, in addition to the high coercivity component. Due to the age of the sample, we cannot exclude the additional possibility of authigenic magnetite (e.g.,
Origin and Timing of Magnetization
Directional information from microbialites with moments <1.00E-10 Am2 is apparently less reliable, likely due to a small but significant contribution from the chosen sample containers (1.00E-13 to 1.00E-11 Am2). As these samples demagnetize, they also approach the stated sensitivity of the magnetometer (∼1.00E-12 Am2). These effects could be reduced in future studies in part by careful choice of sample holders.
Our results for microbialite ChRMs with sufficiently high magnetization appear to record directions close to the expected field direction, with some caveats. The Laguna Bacalar NRM directions are consistent with the expected field direction but have large uncertainties. We interpret these uncertainties to arise from the relatively small number of specimens measured and the lack of AF cleaning which would remove any overprints. While the inclination recorded by the Bayan Gol formation stromatolite is consistent with that reported by
Despite these limitations, an attraction of using ancient microbialites is to potentially recover high temporal resolution variations of geomagnetic field, and it is therefore important to assess the timing and origin of the acquired magnetization. Is the magnetization detrital, locked in soon after deposition and incorporated into the microbialite matrix? Or is it a chemical remanence (biomediated or otherwise) acquired soon after deposition? Or is it related to later diagenesis?
The observed relatively homogenous magnetization along microbialite laminations and variations across laminations suggests a magnetic mineral concentration and magnetization distribution associated with deposition. All studied microbialites seem to preserve this magnetization pattern.
The observed deviation of paleomagnetic directions away from expected values (declination ±10°, inclination ±20°) may be caused by deposition on sloping beds. King (1955) and Kodama (2012) proposed that when a bed dips (shallowly) in the same direction as the magnetic field, the magnetization vector would be steeper by the degree of dip gradient. If the bedding slope is opposite to the Earth’s magnetic field, the inclination can be shallower. Many microbialites have domal internal structures, hence, magnetic directions recorded by different parts of the structure might be disturbed by the effect of deposition on a sloping bed, although there are some obvious differences between normal sedimentary deposition and incorporation into a filamentous microbial matrix. Sample GSL18_0301 was from the (very) steep side of the columnar mound structure, with a slope opposite to the magnetic field. The sample Fisher mean inclination of 43.9° is 16.1° shallower than the geocentric axial dipole (GAD) inclination, and the magnetic vectors may have experienced a bedding-related shallowing. Data from partially oriented sample BG14_0101 were rotated (Figure 8B) so that the mean declination corresponds to geographic south (expected declination for reverse polarity). In this orientation, the left side of the stromatolite (Figure 8A) dips down toward geographic south, opposite to the field direction which is up (negative inclination) and has an inclination shallower than total mean of −66.3° (Figure 8B). By contrast, the right side on Figure 8A dips down to geographic north, roughly parallel to the field, and has a steeper inclination. Inclinations from right side of the structure are −75.0° to −65.0° (Fisher mean = −69.0° ± 3.3°), while the specimens from left side are −65.0° to −55.0° (Fisher mean = −63.1° ± 8.3°), which might support the bedding slope effect. However, Watson’s F-test (Watson, 1956) and Vw-test (Watson, 1983) both show that these means are not statistically distinct at the 95% confidence level. In contrast to these two samples which show inclination deviations, all GSL 19 samples are small structures (<40 cm diameter), and specimens from throughout the structure were used and averaged out. The push core sample LB19_0104 was from the middle of the structure with nearly horizontal layering, and here the Fisher mean inclination of 36.8° ± 7.2° is close to the expected GAD inclination of 33.94°. While not conclusive, these inclination variations associated with the slope of the microbialite layering are intriguing and point to preservation of a primary DRM.
DRM is a relatively inefficient magnetization process, so the DRM/ARM ratio should be relatively low. While not definitive, a ratio of NRM/ARM <1 is consistent with a DRM. Clastic sediments have a wide NRM/ARM range from ∼0.1 to 0.8 (e.g., Levi and Banerjee, 1976; King et al., 1983;
Figure 9 shows NRM vs ARM intensity data for this study, along with NRM/ARM bounds for clastic sediments (0.1–0.8) and marine biogenic magnetite (<0.05). The geometric mean of the NRM/ARM of GSL18 living microbialites, GSL18_0301, and GSL19_0202, are 0.98, 0.57, and 0.37, respectively. LB19 living microbialites and LB19_0301 have a geometric mean of 0.41 and 0.52, respectively. BG14_0101 has a geometric mean NRM/ARM ratio of 0.19, the lowest value among all microbialites. Most specimens are within the clastic detrital range with only one BG14_0101 specimen less than 0.05. These data are fully consistent with a predominantly clastic detrital origin, but some uncertainty in ratio for a biogenic contribution(s) and/or a CRM leaves room for other interpretations. While more work needs to be done to examine the magnetization process and the possible effects of diagenesis, the internal distribution of magnetization with respect to microbialite layering, NRM/ARM ratios, and the inclination variations correlated with bedding slopes suggest that the timing of sediment deposition and acquisition of magnetization samples are close and may be consistent with a primary DRM.
FIGURE 9

NRM intensity (x-axis) vs. ARM intensity (y-axis) of microbialites on log scales. The sample group of GSL18_living (GSL18 living microbialites with cyanobacteria layer collected from Antelope Island Site 1, 2 and 4), GSL18_0301 (lithified microbialite from Antelope Island Site 3), GSL19_0202 (lithified modern microbialite from Lakeside Site 2), LB19_living (surface living microbialite samples collected from LB19 Site 1 and 3), LB19_0301 (lithified porous microbialite sample from LB19 Site 3), BG14_0101 (Cambrian Bayan Gol formation stromatolite) are shown with red, orange, yellow, green, blue, and purple circles, and their geometric means are symbolized with stars. The black solid line represents NRM to ARM ratio of 1 (x = y), and measurements plotted on the left side of this line (ARM >> NRM) is consistent with a DRM origin. The gray shading shows the NRM/ARM range of clastic sediments (0.1–0.8) (e.g., Levi and Banerjee, 1976;
Conclusion
We investigated the potential of modern and ancient microbialites as paleomagnetic recorders. Microbialites with NRM >1.00E-10 Am2 record recoverable, stable magnetic directions. All microbialites appear to record average directions close to the expected field direction. The dominant NRM carrier in all samples is 35–50 mT SD to ‘PSD’ magnetite. The Cambrian BG sample also has a contribution from high coercivity minerals, which do not contribute to the NRM.
We interpret the NRM to be at least partially detrital in origin. This interpretation is based on the magnetization distribution in the structure, the NRM/ARM ratios consistent with detrital sediments, and inclination variations that correlate with changes in bedding slopes. In some cases, the NRM to ARM ratio is higher, suggesting another magnetization process. This may be due to biomineralization of magnetic minerals and/or diagenesis soon after deposition.
To reduce inclination variations resulting from sloping beds, samples from large domal structures should be collected avoiding steep sides or by sampling all sides. For smaller scale hand samples (<50 cm), the sample slice should be cut parallel to the expected direction of the magnetic declination. Then, specimens from different lamination angles should be collected to average out the deviations generated by sloping beds.
Statements
Data availability statement
Raw data supporting the conclusions of this article are available in the Magnetics Information Consortium (MagIC) database at earthref.org/MagIC/17087.
Author contributions
J-IJ participated in GSL and LB field work, and collected and analyzed the data. JB designed the study and participated in GSL field work. Both authors contributed to the writing and editing of the manuscript.
Funding
This work was supported by a Research and Creative Activities Support (RACAS) award from the University of Wisconsin–Milwaukee (project AAD3431).
Acknowledgments
We thank Bonnie Baxter and Jaimi Butler of Westminster College who provided assistance and education in sampling microbialites on Antelope Island. These samples were collected under permit provided to Westminster College by the Division of Parks and Recreation/Antelope Island State Park. We thank Tim Grundl for help with sampling at Laguna Bacalar, Steve Dornbos who provided the BG sample, and Joshua Marquardt for preliminary measurements on the BG sample. Many thanks to Andrew Roberts and Nick Swanson‐Hysell whose comments helped improve this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2021.603805/full#supplementary-material.
References
1
AllwoodA. C.WalterM. R.KamberB. S.MarshallC. P.BurchI. W. (2006). Stromatolite reef from the early archaean era of Australia. Nature441, 714–718. 10.1038/nature04764
2
AwramikS. M. (2006). Palaeontology: respect for stromatolites. Nature441, 700–701. 10.1038/44170010.1038/441700a
3
AwramikS. M.GreyK. (2005). Stromatolites: biogenicity, biosignatures, and bioconfusion. Proc. SPIE3906, 5906P1–5906P9. 10.1117/12.625556
4
AwramikS. M.SprinkleJ. (1999). Proterozoic stromatolites: the first marine evolutionary biota. Hist. Biol.13, 241–253. 10.1080/08912969909386584
5
AwramikS. M. (1992). The oldest records of photosynthesis. Photosyn Res.33, 75–89. 10.1007/bf00039172
6
BaumgartnerR. J.KranendonkM. J.WaceyD.FiorentiniM. L.SaundersM.CarusoS.et al (2019). Nano−porous pyrite and organic matter in 3.5-billion-year-old stromatolites record primordial life. Geology47, 1039–1043. 10.1130/g46365.1
7
BerelsonW. M.CorsettiF. A.Pepe-RanneyC.HammondD. E.BeaumontW.SpearJ. R. (2011). Hot spring siliceous stromatolites from Yellowstone National Park: assessing growth rate and laminae formation. Geobiology9, 411–424. 10.1111/j.1472-4669.2011.00288.x
8
BergM. D. (2019). Domes, rings, ridges, and polygons: characteristics of microbialites from Utah's Great Salt Lake. Sediment. Rec.17, 4–10. 10.2110/sedred.2019.1.4
9
BigginA. J.PiispaE. J.PesonenL. J.HolmeR.PatersonG. A.VeikkolainenT.et al (2015). Paleomagnetic field intensity variations suggest Mesoproterozoic inner-core nucleation. Nature526, 245–248. 10.1038/nature15523
10
BonoR. K.TardunoJ. A.NimmoF.CottrellR. D. (2019). Young inner core inferred from Ediacaran ultra-low geomagnetic field intensity. Nat. Geosci.12, 143–147. 10.1038/s41561-018-0288-0
11
BourneM. D.FeinbergJ. M.StraussB. E.HardtB.ChengH.RoweH. D.et al (2015). Long-term changes in precipitation recorded by magnetic minerals in speleothems. Geology43, 595–598. 10.1130/g36695.1
12
BradyA. L.SlaterG.LavalB.LimD. S. (2009). Constraining carbon sources and growth rates of freshwater microbialites in Pavilion Lake using (14)C analysis. Geobiology7, 544–555. 10.1111/j.1472-4669.2009.0021510.1111/j.1472-4669.2009.00215.x
13
BrasierM.GreenO.LindsayJ.SteeleA. (2004). “Earth’s oldest (approximately 3.5 Ga) fossils and the “Early Eden hypothesis” questioning the evidence” in Origins of Life and Evolution of the Biosphere. Editor LuisiP. L. (Dordrecht: Kluwer), 257–269. 10.1023/B:ORIG.0000009845.62244.d3
14
BurneR. V.MooreL. S. (1987). Microbialites: organosedimentary deposits of benthic microbial communities. Palaios2, 241. 10.2307/3514674
15
Castro-ContrerasS. I.GingrasM. K.PecoitsE.AubetN. R.PetrashD.Castro-ContrerasS. M.et al (2014). Textural and geochemical features of freshwater microbialites from Laguna bacalar, Quintana Roo, Mexico. Palaios29, 192–209. 10.2110/palo.2013.063
16
ChangS. R.StolzJ. F.KirschvinkJ. L.AwramikS. M. (1989). Biogenic magnetite in stromatolites. II. Occurrence in ancient sedimentary environments. Precambrian Res.43, 305–315. 10.1016/0301-9268(89)90062-4
17
ChidseyT. C.BergM. D.EbyD. E. (2015). Petrography and characterization of microbial carbonates and associated facies from modern Great Salt Lake and Uinta Basin's Eocene Green River Formation in Utah, USA. London: Geological Society, London, Special Publications,418, 261-286. 10.1144/sp418.6
18
ChouY. M.JiangX.LiuQ.HuH. M.WuC. C.LiuJ.et al (2018). Multidecadally resolved polarity oscillations during a geomagnetic excursion. Proc. Natl. Acad. Sci. U S A.115, 8913–8918. 10.1073/pnas.1720404115
19
ConstableC. G. (1985). Eastern Australian geomagnetic field intensity over the past 14000 yr. Geophys. J. Int.81, 121–130. 10.1111/j.1365-246x.1985.tb01354.x
20
ConstableC.KorteM.PanovskaS. (2016). Persistent high paleosecular variation activity in southern hemisphere for at least 10 000 years. Earth Planet. Sci. Lett.453, 78–86. 10.1016/j.epsl.2016.08.015
21
EgliR. (2004). Characterization of individual rock magnetic components by analysis of remanence curves, 1. Unmixing Natural Sediments. Studia Geophysica et Geodaetica48, 391–446. 10.1023/B:SGEG.0000020839.45304.6d
22
EgliR.LowrieW. (2002). Anhysteretic remanent magnetization of fine magnetic particles. J. Geophys. Res.107, 2209. 10.1029/2001JB000671
23
EgliR.ZhaoX. (2015). Natural remanent magnetization acquisition in bioturbated sediment: general theory and implications for relative paleointensity reconstructions. Geochem. Geophys. Geosystems16, 995–1016. 10.1002/2014GC005672
24
EvansD. A.ZhuravlevA. Y.BudneyC. J.KirschvinkJ. L. (1996). Paleomagnetism of the bayan Gol formation, western Mongolia. Geol. Mag.133, 487–496. 10.1017/s0016756800007615
25
FisherR. (1953). Dispersion on a sphere. Proc. R. Soc. A: Math. Phys. Eng. Sci.217, 295–305. 10.1098/rspa.1953.0064
26
FontE.TrindadeR. I.NédélecA. (2005). Detrital remanent magnetization in haematite-bearing Neoproterozoic Puga cap dolostone, Amazon craton: a rock magnetic and SEM study. Geophys. J. Int.163, 491–500. 10.1111/j.1365-246x.2005.02776.x
27
GischlerE.GibsonM. A.OschmannW. (2008). Giant Holocene freshwater microbialites, Laguna bacalar, Quintana Roo, Mexico. Sedimentology55, 1293–1309. 10.1111/j.1365-3091.2007.00946.x
28
GoldenD.MingD.MorrisR.BrearleyA.LauerH.TreimanA.et al (2004). Evidence for exclusively inorganic formation of magnetite in Martian meteorite ALH84001. Am. Mineral.89, 681–695. 10.2138/am-2004-5-602
29
JacksonM. (1990). Diagenetic sources of stable remanence in remagnetized paleozoic cratonic carbonates: a rock magnetic study. J. Geophys. Res.95, 2753. 10.1029/jb095ib03p02753
30
JacksonM.Swanson-HysellN. L. (2012). Rock magnetism of remagnetized carbonate rocks: another look. Geol. Soc. Lond. Spec. Publ371, 229–251. 10.1144/sp371.3
31
JacksonM.Van der VooR. (1985). A lower ordovician paleomagnetic pole from the Oneida dolomite, upper Mississippi river valley. J. Geophys. Res.90, 10449. 10.1029/jb090ib12p10449
32
KatariK.TauxeL. (2000). Effects of pH and salinity on the intensity of magnetization in redeposited sediments. Earth Planet. Sci. Lett.181, 489–496. 10.1016/S0012-821X(00)00226-0
33
KhomentovskyV. V.GibsherA. S. (1996). The Neoproterozoic–lower Cambrian in northern Govi-Altay, western Mongolia: Regional setting, lithostratigraphy and biostratigraphy. Geolog. Mag,133, 371-390. 10.1017/s001675680000755xSociety for Sedimentary Geology.
34
KiesslingW. (2002). Secular variations in the phanerozoic reef ecosystem. Soc. Sediment. Geology.72, 625–690. 10.2110/pec.02.72.0625
35
KingJ. W.BanerjeeS. K.MarvinJ. (1983). A new rock-magnetic approach to selecting sediments for geomagnetic paleointensity studies: application to paleointensity for the last 4000 years. J. Geophys. Res. Solid Earth88, 5911–5921. 10.1029/jb088ib07p05911
36
KingR. F. (1955). The remanent magnetism of artificially deposited sediments. Geophys. J. Int.7, 115–134. 10.1111/j.1365-246x.1955.tb06558.x
37
KirschvinkJ. L. (1980). The least-squares line and plane and the analysis of paleomagnetic data. Geophys. J. Int.62, 699–718. 10.1111/j.1365-246x.1980.tb02601.x
38
KodamaK. P. (2012). Paleomagnetism of sedimentary rocks: process and interpretation. Chichester: Wiley-Blackwell.
39
KonhauserK. O.KapplerA.RodenE. E. (2011). Iron in microbial metabolisms. Elements7, 89–93. 10.2113/gselements.7.2.89
40
KonhauserK.RidingR. (2012). “Bacterial biomineralization,” in Fundamentals of geobiology. Editors KnollA. H.CanfieldD. E.KonhauserK. O. (Chichester: Wiley-Blackwell)), 105–130. 10.1002/9781118280874.ch8
41
Kądziałko-HofmoklM.KruczykJ. (1987). Paleomagnetism of Middle-Late Jurassic sediments from Poland and implications for the polarity of the geomagnetic field. Tectonophysics139, 53–66. 10.1016/0040-1951(87)90197-1
42
LascuI.FeinbergJ. M. (2011). Speleothem magnetism. Quat. Sci. Rev.30, 3306–3320. 10.1016/j.quascirev.2011.08.004
43
LeviS.BanerjeeS. K. (1976). On the possibility of obtaining relative paleointensities from lake sediments. Earth Planet. Sci. Lett.29, 219–226. 10.1016/0012-821x(76)90042-x
44
LovleyD. R. (1991). “Magnetite formation during microbial dissimilatory iron reduction,” in Iron biominerals. Editors FrankelR. B.BlakemoreR. P. (Boston, MA: Springer), 151–166. 10.1007/978-1-4615-3810-3_11
45
LoweD. R. (1994). Abiological origin of described stromatolites older than 3.2 Ga. Geology22, 387. 10.1130/0091-7613(1994)0222.3.co;2
46
LowrieW. (1990). Identification of ferromagnetic minerals in a rock by coercivity and unblocking temperature properties. Geophys. Res. Lett.17, 159–162. 10.1029/gl017i002p00159
47
LundS.PlatzmanE.ThouvenyN.CamoinG.CorsettiF.BerelsonW. (2010). Biological control of paleomagnetic remanence acquisition in carbonate framework rocks of the Tahiti coral reef. Earth Planet. Sci. Lett.298, 14–22. 10.1016/j.epsl.2010.07.010
48
MaherB.TaylorR. (1988). Formation of ultrafine-grained magnetite in soils. Nature336, 368–370. 10.1038/336368a0
49
MaloofA. C.KoppR. E.GrotzingerJ. P.FikeD. A.BosakT.ValiH.et al (2007). Sedimentary iron cycling and the origin and preservation of magnetization in platform carbonate muds, Andros Island, Bahamas. Earth Planet. Sci. Lett.259, 581–598. 10.1016/j.epsl.2007.05.021
50
MaxbauerD. P.FeinbergJ. M.FoxD. L. (2016). MAX UnMix: a web application for unmixing magnetic coercivity distributions. Comput. Geosciences95, 140–145. 10.1016/j.cageo.2016.07.009
51
McCabeC.ChannellJ. E. (1994). Late Paleozoic remagnetization in limestones of the Craven Basin (northern England) and the rock magnetic fingerprint of remagnetized sedimentary carbonates. J. Geophys. Res. Solid Earth99, 4603–4612. 10.1029/93jb02802
52
McFaddenP. L.ReidA. B. (1982). Analysis of paleomagnetic inclination data. Geophys. J. Int.69, 307–319. 10.1111/j.1365-246X.1982.tb04950.x\
53
McLoughlinN.WilsonL. A.BrasierM. D. (2008). Growth of synthetic stromatolites and wrinkle structures in the absence of microbes - implications for the early fossil record. Geobiology6, 95–105. 10.1111/j.1472-4669.2007.00141.x
54
MitraR.TauxeL. (2009). Full vector model for magnetization in sediments. Earth Planet. Sci. Lett.286, 535–545. 10.1016/j.epsl.2009.07.019
55
MuraszkoJ. R. (2014). Magnetic properties of Jurassic stromatolites from selected localities in Poland. [master’s thesis]. [Poland]: University of Warsaw.
56
NewellD. L.JensenJ. L.FrantzC. M.BergM. D. (2017). Great Salt Lake (Utah) Microbialite δ13C, δ18O, and δ15N record fluctuations in lake biogeochemistry since the late Pleistocene. Geochem. Geophys. Geosystems18, 3631–3645. 10.1002/2017gc007078
57
PanY.ZhuR.BanerjeeS. K.GillJ.WilliamsQ. (2000). Rock magnetic properties related to thermal treatment of siderite: behavior and interpretation. J. Geophys. Res. Solid Earth105, 783–794. 10.1029/1999jb900358
58
PaullC.NeumannA.BeboutB.ZabielskiV.ShowersW. (1992). Growth rate and stable isotopic character of modern stromatolites from San Salvador, Bahamas. Palaeogeogr. Palaeoclimatol. Palaeoecol.95, 335–344. 10.1016/0031-0182(92)90149-y
59
PetryshynV. A.CorsettiF. A.BerelsonW. M.BeaumontW.LundS. P. (2012). Stromatolite lamination frequency, Walker Lake, Nevada: implications for stromatolites as biosignatures. Geology40, 499–502. 10.1130/g32675.1
60
PetryshynV. A.CorsettiF. A.FrantzC. M.LundS. P.BerelsonW. M. (2016). Magnetic susceptibility as a biosignature in stromatolites. Earth Planet. Sci. Lett.437, 66–75. 10.1016/j.epsl.2015.12.016
61
PickT.TauxeL. (1991). Chemical remanent magnetization in synthetic magnetite. J. Geophys. Res.96, 9925–9936. 10.1029/91JB00706
62
RasmussenK. A.MacintyreI. G.PrufertL. (1993). Modern stromatolite reefs fringing a brackish coastline, Chetumal Bay, Belize. Geology21, 199. 10.1130/0091-7613(1993)0212.3.co;2
63
ReidR. P.VisscherP. T.DechoA. W.StolzJ. F.BeboutB. M.DuprazC.et al (2000). The role of microbes in accretion, lamination and early lithification of modern marine stromatolites. Nature406, 989–992. 10.1038/35023158
64
RidingR. (2006). Microbial carbonate abundance compared with fluctuations in metazoan diversity over geological time. Sediment. Geology.185, 229–238. 10.1016/j.sedgeo.2005.12.01
65
RidingR. (2000). Microbial carbonates: the geological record of calcified bacterial-algal mats and biofilms. Sedimentology47, 179–214. 10.1046/j.1365-3091.2000.00003.x
66
RidingR. (2011). “Microbialites, stromatolites, and thrombolites,” in Encyclopedia of geobiology encyclopedia of Earth sciences series. Editors ReitnerJ.ThielV. (Dordrecht: Springer), 635–654. 10.1007/978-1-4020-9212-1_196
67
RidingR. (2010). The nature of stromatolites: 3,500 million years of history and a century of research. Adv. Stromatolite Geobiology Lecture Notes Earth Sci.131, 29–74. 10.1007/978-3-642-10415-2_3
68
RobertsA. P.ChangL.HeslopD.FlorindoF.LarrasoañaJ. C. (2012). Searching for single domain magnetite in the “pseudo-single-domain” sedimentary haystack: implications of biogenic magnetite preservation for sediment magnetism and relative paleointensity determinations. J. Geophys. Res. Solid Earth117. 10.1029/2012jb009412
69
RobertsA. P. (2015). Magnetic mineral diagenesis. Earth-Science Rev.151, 1–47. 10.1016/j.earscirev.2015.09.010
70
RobertsonD.FranceD. (1994). Discrimination of remanence-carrying minerals in mixtures, using isothermal remanent magnetisation acquisition curves. Phys. Earth Planet. Interiors82, 223–234. 10.1016/0031-9201(94)90074-4
71
RodelliD.JovaneL.GiorgioniM.RegoE. S.CornaggiaF.BenitesM.et al (2019). Diagenetic fate of biogenic soft and hard magnetite in chemically stratified sedimentary environments of mamanguá ría, Brazil. J. Geophys. Res. Solid Earth124, 2313–2330. 10.1029/2018jb016576
72
RyanW. B.CarbotteS. M.CoplanJ. O.O'haraS.MelkonianA.ArkoR.et al (2009). Global multi-resolution topography (GRMT) synthesis. Geochem. Geophys. Geosystems10. 10.1029/2008gc002332
73
SchopfJ. W.KudryavtsevA. B.CzajaA. D.TripathiA. B. (2007). Evidence of archean life: stromatolites and microfossils. Precambrian Res.158, 141–155. 10.1016/j.precamres.2007.04.009
74
SmirnovA. V.TardunoJ. A.KulakovE. V.McenroeS. A.BonoR. K. (2016). Palaeointensity, core thermal conductivity and the unknown age of the inner core. Geophys. J. Int.205, 1190–1195. 10.1093/gji/ggw080
75
SnowballI.SandgrenaP. (2002). Geomagnetic field variations in northern Sweden during the Holocene quantified from varved lake sediments and their implications for cosmogenic nuclide production rates. The Holocene12, 517–530. 10.1191/0959683602hl562rp
76
StolzJ. F.ChangS. R.KirschvinkJ. L. (1989). Biogenic magnetite in stromatolites. I. Occurrence in modern sedimentary environments. Precambrian Res.43, 295–304. 10.1016/0301-9268(89)90061-2
77
SukD.HalgedahlS. L. (1996). Hysteresis properties of magnetic spherules and whole rock specimens from some Paleozoic platform carbonate rocks. J. Geophys. Res. Solid Earth101, 25053–25075. 10.1029/96jb02271
78
SukD.PeacorD. R.VooR. V. (1990a). Replacement of pyrite framboids by magnetite in limestone and implications for paleomagnetism. Nature345, 611–613. 10.1038/345611a0
79
SukD.VooR. V.PeacorD. R. (1990b). Scanning and transmission electron microscope observations of magnetite and other iron phases in Ordovician carbonates from east Tennessee. J. Geophys. Res. Solid Earth95, 12327–12336. 10.1029/JB095iB08p1232
80
Swanson-HysellN. L.MaloofA. C.KirschvinkJ. L.EvansD. A.HalversonG. P.HurtgenM. T. (2012). Constraints on neoproterozoic paleogeography and paleozoic orogenesis from paleomagnetic records of the bitter Springs Formation, amadeus basin, central Australia. Am. J. Sci.312, 817–884. 10.2475/08.2012.01
81
TardunoJ. A.CottrellR. D.SmirnovA. V. (2006). The paleomagnetism of single silicate crystals: recording geomagnetic field strength during mixed polarity intervals, superchrons, and inner core growth. Rev. Geophys.44. 10.1029/2005rg000189
82
TauxeL.ShaarR.JonestraskL.Swanson‐HysellN. L.MinnettR.KoppersA. A.et al (2016). PmagPy: software package for paleomagnetic data analysis and a bridge to the magnetics information Consortium (MagIC) database. Geochem. Geophys. Geosystems17, 2450–2463. 10.1002/2016gc006307
83
TaylorR. M.MaherB. A.SelfP. G. (1987). Magnetite in soils: 1. The synthesis of single-domain and superparamagnetic magnetite. Clay Minerals22, 411–422. 10.1180/claymin.1987.022.4.05
84
ThompsonR.OldfieidF. (1986). Environmental magnetism. Dordrecht: Springer. 10.1007/978-94-011-8036-8
85
TillJ. L.GuyodoY.LagroixF.MorinG.MenguyN.Ona-NguemaG. (2017). Presumed magnetic biosignatures observed in magnetite derived from abiotic reductive alteration of nanogoethite. Comptes Rendus Geosci.349, 63–70. 10.1016/j.crte.2017.02.001
86
TrindadeR. I.FontE.D'agrella-FilhoM. S.NogueiraA. C.RiccominiC. (2003). Low-latitude and multiple geomagnetic reversals in the Neoproterozoic Puga cap carbonate, Amazon craton. Terra Nova15, 441–446. 10.1046/j.1365-3121.2003.00510.x
87
TrindadeR. I.JaquetoP.Terra-NovaF.BrandtD.HartmannG. A.FeinbergJ. M.et al (2018). Speleothem record of geomagnetic south atlantic anomaly recurrence. Proc. Natl. Acad. Sci.115, 13198–13203. 10.1073/pnas.1809197115
88
ValetJ.BesseJ.KumarA.Vadakke-ChanatS.PhilippeE. (2014). The intensity of the geomagnetic field from 2.4 Ga old Indian dykes. Geochem. Geophys. Geosystems15, 2426–2437. 10.1002/2014gc005296
89
VanyoJ. P.AwramikS. M. (1982). Length of day and obliquity of the ecliptic 850 MA ago: preliminary results of a stromatolite growth model. Geophys. Res. Lett.9, 1125–1128. 10.1029/gl009i010p01125
90
WalterM. R.HeysG. R. (1985). Links between the rise of the metazoan and the decline of stromatolites. Precambrian Res.29, 149–174. 10.1016/0301-9268(85)90066-X
91
WatsonG. S. (1956). Analysis of dispersion on a sphere. Geophys. J. Int.7, 153–159. 10.1111/j.1365-246x.1956.tb05560.x
92
WatsonG. S. (1983). Large sample theory of the Langevin distribution. J. Stat. Plann. Inference8, 245–256. 10.1016/0378-3758(83)90043-5
93
WeilA. B.Van der VooR. (2002). Insights into the mechanism for orogen‐related carbonate remagnetization from growth of authigenic Fe‐oxide: a scanning electron microscopy and rock magnetic study of Devonian carbonates from northern Spain. J. Geophys. Res. Solid Earth107, 1–14. 10.1029/2001JB000200
94
ZhaoX.EgliR.GilderS.MüllerS. (2016). Microbial assisted recording of the Earth's magnetic field in sediment. Nat. Commun.7, 10673. 10.1038/ncomms10673
Summary
Keywords
microbialites, stromatolites, paleomagnetism, rock magnetism, paleosecular variation
Citation
Jung J-I and Bowles JA (2021) A Feasibility Study of Microbialites as Paleomagnetic Recorders. Front. Earth Sci. 9:603805. doi: 10.3389/feart.2021.603805
Received
08 September 2020
Accepted
25 January 2021
Published
05 March 2021
Volume
9 - 2021
Edited by
Eric C. Ferre, University of Louisiana at Lafayette, United States
Reviewed by
Andrew Philip Roberts, Australian National University, Australia
Nicholas L. Swanson-Hysell, University of California, Berkeley, United States
Updates

Check for updates
Copyright
© 2021 Jung and Bowles.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ji-In Jung, jiinjung@stanford.edu
This article was submitted to Geomagnetism and Paleomagnetism, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.