Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

Timing and tempo of the Great Oxidation Event

View ORCID ProfileAshley P. Gumsley, Kevin R. Chamberlain, Wouter Bleeker, Ulf Söderlund, Michiel O. de Kock, Emilie R. Larsson, and Andrey Bekker
  1. aDepartment of Geology, Lund University, Lund 223 62, Sweden;
  2. bDepartment of Geology and Geophysics, University of Wyoming, Laramie, WY 82071;
  3. cFaculty of Geology and Geography, Tomsk State University, Tomsk 634050, Russia;
  4. dGeological Survey of Canada, Ottawa, ON K1A 0E8, Canada;
  5. eDepartment of Geosciences, Swedish Museum of Natural History, Stockholm 104 05, Sweden;
  6. fDepartment of Geology, University of Johannesburg, Auckland Park 2006, South Africa;
  7. gDepartment of Earth Sciences, University of California, Riverside, CA 92521

See allHide authors and affiliations

PNAS February 21, 2017 114 (8) 1811-1816; first published February 6, 2017; https://doi.org/10.1073/pnas.1608824114
Ashley P. Gumsley
aDepartment of Geology, Lund University, Lund 223 62, Sweden;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ashley P. Gumsley
  • For correspondence: ashley.gumsley@geol.lu.se
Kevin R. Chamberlain
bDepartment of Geology and Geophysics, University of Wyoming, Laramie, WY 82071;
cFaculty of Geology and Geography, Tomsk State University, Tomsk 634050, Russia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wouter Bleeker
dGeological Survey of Canada, Ottawa, ON K1A 0E8, Canada;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ulf Söderlund
aDepartment of Geology, Lund University, Lund 223 62, Sweden;
eDepartment of Geosciences, Swedish Museum of Natural History, Stockholm 104 05, Sweden;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michiel O. de Kock
fDepartment of Geology, University of Johannesburg, Auckland Park 2006, South Africa;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Emilie R. Larsson
aDepartment of Geology, Lund University, Lund 223 62, Sweden;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrey Bekker
gDepartment of Earth Sciences, University of California, Riverside, CA 92521
fDepartment of Geology, University of Johannesburg, Auckland Park 2006, South Africa;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Edited by Mark H. Thiemens, University of California, San Diego, La Jolla, CA, and approved December 27, 2016 (received for review June 11, 2016)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

We present U-Pb ages for the extensive Ongeluk large igneous province, a large-scale magmatic event that took place near the equator in the Paleoproterozoic Transvaal basin of southern Africa at ca. 2,426 Ma. This magmatism also dates the oldest Paleoproterozoic global glaciation and the onset of significant atmospheric oxygenation. This result forces a significant reinterpretation of the iconic Transvaal basin stratigraphy and implies that the oxygenation involved several oscillations in oxygen levels across 10−5 present atmospheric levels before the irreversible oxygenation of the atmosphere. Data also indicate that the Paleoproterozoic glaciations and oxygenation were ushered in by assembly of a large continental mass, extensive magmatism, and continental migration to near-equatorial latitudes, mirroring a similar chain of events in the Neoproterozoic.

Abstract

The first significant buildup in atmospheric oxygen, the Great Oxidation Event (GOE), began in the early Paleoproterozoic in association with global glaciations and continued until the end of the Lomagundi carbon isotope excursion ca. 2,060 Ma. The exact timing of and relationships among these events are debated because of poor age constraints and contradictory stratigraphic correlations. Here, we show that the first Paleoproterozoic global glaciation and the onset of the GOE occurred between ca. 2,460 and 2,426 Ma, ∼100 My earlier than previously estimated, based on an age of 2,426 ± 3 Ma for Ongeluk Formation magmatism from the Kaapvaal Craton of southern Africa. This age helps define a key paleomagnetic pole that positions the Kaapvaal Craton at equatorial latitudes of 11° ± 6° at this time. Furthermore, the rise of atmospheric oxygen was not monotonic, but was instead characterized by oscillations, which together with climatic instabilities may have continued over the next ∼200 My until ≤2,250–2,240 Ma. Ongeluk Formation volcanism at ca. 2,426 Ma was part of a large igneous province (LIP) and represents a waning stage in the emplacement of several temporally discrete LIPs across a large low-latitude continental landmass. These LIPs played critical, albeit complex, roles in the rise of oxygen and in both initiating and terminating global glaciations. This series of events invites comparison with the Neoproterozoic oxygen increase and Sturtian Snowball Earth glaciation, which accompanied emplacement of LIPs across supercontinent Rodinia, also positioned at low latitude.

  • Great Oxidation Event
  • Snowball Earth
  • Paleoproterozoic
  • Kaapvaal Craton
  • Transvaal Supergroup

The early Paleoproterozoic is characterized by dramatic changes in Earth’s atmosphere and oceans, with the transition from anoxic to oxic conditions commonly referred to as the Great Oxidation Event (GOE) (1). It is generally thought that the onset of the GOE was a singular event (2), an assumption rooted in the perceived bistability of atmospheric oxygen (3). However, this inferred bistability in oxygen was challenged through additional modeling (4), allowing for multiple oscillations in atmospheric oxygen during the onset of the GOE. Geological evidence has also established that this transition was broadly coincident with emplacement of numerous large igneous provinces (LIPs) (5) on extensive continental landmasses positioned at low latitudes (6) and glaciations interpreted to reflect global Snowball Earth conditions (7). Models linking these events have been hampered, however, by uncertainties in local and global stratigraphic correlations and age constraints (2). Evidence from the Huronian Supergroup on the Superior Craton in Canada, which hosts three Paleoproterozoic glacial intervals, indicates that the GOE is bracketed in age between ca. 2,460 and 2,308 Ma (8⇓–10). New observations from the critical Transvaal Supergroup in southern Africa indicate that the GOE may have occurred by ca. 2,310 Ma (9) but requires a fourth glaciation at ≤2,250–2,240 Ma (9, 11) before the development of the oldest widely accepted oxygenated paleosols (12) and red-bed sandstones (13). Discrepancies in correlations between these glacially influenced stratigraphic successions and the age of the GOE hinge on a disputed age for the Ongeluk Formation basalts (14⇓–16). By using high-resolution in situ secondary ion mass spectrometry (SIMS) on microbaddeleyite grains coupled with precise isotope dilution thermal ionization mass spectrometry (ID-TIMS) and paleomagnetic studies, we resolve these uncertainties by obtaining accurate and precise ages for the volcanic Ongeluk Formation and related intrusions in South Africa. These ages lead to a more coherent global perspective on the timing and tempo of the GOE and associated global glaciations and LIPs.

Transvaal Supergroup

The Neoarchean to Paleoproterozoic Transvaal Supergroup overlies the Kaapvaal Craton and is preserved in two main subbasins. In the Griqualand West subbasin, the Makganyene Formation consists of a series of glaciomarine diamictites (Fig. 1) (17). Paleomagnetic data for the Ongeluk Formation, which conformably overlies and interfingers with the Makganyene Formation, indicate that these glacial sedimentary rocks were deposited at low latitude, implying a glacial event of global extent (6, 7). The age of the Ongeluk Formation basalts has long been accepted at 2,222 ± 13 Ma based on a whole-rock Pb-Pb isochron date (14) and correlation with the noncontiguous basalts of the ≤2,250–2,240 Ma Hekpoort Formation in the Transvaal subbasin (11, 12, 14, 18), although this correlation has been questioned recently (15, 16, 19). The Ongeluk Formation is overlain by banded iron and manganese deposits of the Hotazel Formation, which in turn, are followed by carbonate rocks of the Mooidraai Formation (Fig. 1). Pb-Pb and U-Pb dating of the Mooidraai Formation carbonates has yielded dates of 2,394 ± 26 and 2,392 ± 23 Ma (15, 20), respectively, in conflict with the 2,222 ± 13 Ma Pb-Pb date on the stratigraphically lower Ongeluk Formation (14). The Hotazel Formation hosts giant Mn deposits with a negative Ce anomaly that are unequivocally interpreted to reflect deposition after the onset of the GOE (7, 21), whereas the Koegas Subgroup underlying the Makganyene Formation contains detrital pyrite and uraninite grains signifying deposition before the GOE (22). In the Transvaal subbasin (Fig. 1), the start of the GOE has been placed in the middle of either the Duitschland Formation or the Rooihoogte Formation by different authors (23, 24), with the age of the upper Duitschland Formation constrained by detrital zircon to ≤2,424 ± 24 Ma (11). A 2,316 ± 7 Ma Re-Os age for diagenetic pyrite (25) and a 2,309 ± 9 Ma U-Pb age for tuff in the lower Timeball Hill Formation (9), which conformably overlies the Rooihoogte Formation, suggest that the GOE began by ca. 2,309 Ma (24). All chronological and redox records for the Transvaal Supergroup are provided in Tables S1 and S2.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Stratigraphic synthesis (SI Methods, Stratigraphic Synthesis) of the Transvaal Supergroup as preserved in its two main subbasins: (A) Griqualand West in the southwest and (B) Transvaal in the northeast. Dated samples and results (bold) are shown in stratigraphic context and collectively, unlock the long-held correlation between the basalts of the Ongeluk and Hekpoort formations (12, 14, 18). A selection of previously published ages is schematically shown (Table S1) along with redox indicators and ranges of carbon isotope values in carbonates (Table S2). The redox records within the two subbasins tracks the rhythm of GOE and reflects at least two O2 oscillations back through 10−5 PAL after the onset of GOE (bold pink arrows and variable intensity pink background shading). Redox indicators requiring more detailed studies are denoted with question marks. All ages are quoted at 2σ uncertainty.

View this table:
  • View inline
  • View popup
Table S1.

Compilation of ages for the Transvaal Supergroup and the Huronian Supergroup

View this table:
  • View inline
  • View popup
Table S2.

Compilation of carbon isotope data and evidence for the redox state of surface environments during deposition of the Transvaal and Huronian Supergroups

Sampling

To test and resolve some of these critical correlations, we have dated by U-Pb isotopic methods a number of dolerite and basalt samples that are linked geologically and paleomagnetically to the Ongeluk Formation basalts (SI Methods, Sampling, Fig. S1, and Table S3). An N-trending dolerite dike from the Griqualand West subbasin (G02-B) (Fig. 1) (26) as well as an intrusive dolerite sheet from the southeastern Kaapvaal Craton (NL-13c) (Fig. 1) (27) were dated using U-Pb ID-TIMS on baddeleyite. Samples TGS-05 and OLL-2, a coarse-grained, thick basalt flow and a dolerite sill, respectively, from near the base of the Ongeluk Formation basalts (Fig. 1) (6) were dated by in situ U-Pb SIMS on microbaddeleyite grains. To couple geochronological and paleomagnetic records for these mafic units, complementary paleomagnetic studies were conducted on specimens from the TGS-05 and MDK-05 sample sites (Fig. 1) using conventional demagnetization techniques.

Fig. S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S1.

The Kaapvaal Craton in southern Africa with the three subbasins of the Transvaal Supergroup shown along with the Archean basement. Sample localities in the Ongeluk Formation and related intrusions from this study and that of the Westerberg Sill Province (28) are also shown.

View this table:
  • View inline
  • View popup
Table S3.

Sampling locality descriptions of the Ongeluk Formation volcanic rocks and related mafic intrusions

Results

Samples G02-B and NL-13c produce upper intercept baddeleyite dates of 2,421 ± 3 and 2,423 ± 7 Ma, respectively (Fig. 2, SI Methods, Geochronology—ID-TIMS Analysis, Fig. S2, and Table S4), whereas samples TGS-05 and OLL-2 yield upper intercept dates of 2,424 ± 32 and 2,397 ± 22 Ma, respectively (Fig. 2, SI Methods, Geochronology—SIMS Analysis, Fig. S2, and Table S5). The upper intercept dates of 2,424 ± 32, 2,423 ± 7, and 2,421 ± 3 Ma are interpreted as reliable crystallization ages and overlap within 2σ uncertainty, whereas the date obtained from OLL-2 is likely affected by a minor contribution of secondary zircon as indicated by SEM imaging. We, therefore, interpret the 2,397 ± 22 Ma date of OLL-2 as a minimum age. In support of this interpretation, one baddeleyite grain (spot F861b) (Fig. 2) yielded a 207Pb/206Pb date of 2,421 ± 22 Ma.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Weighted mean age of the Ongeluk LIP. Shown is a comparison of upper intercept dates with 2σ uncertainties (red columns) from five samples of the Ongeluk LIP, including the Westerberg Sill Province (samples TGS-01 and M03WA) (28), with a calculated weighted mean age of 2,425.5 ± 2.6 Ma (green bar). The result from a single analysis spot (F861b in OLL-2) is shown for comparison (blue column) (Fig. S2).

Fig. S2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S2.

(A) ID-TIMS concordia diagrams for baddeleyite from samples NL-13c and G02-B, with data ellipses at 2σ uncertainty. (B) In situ SIMS concordia diagrams of microbaddeleyite in samples OLL-2 and TGS-05, with data ellipses at 1σ uncertainty. In addition, (B, Left) polished thin sections (with SIMS analysis points used in the age calculations) were mapped by SEM to locate and identify (B, Right) Zr-bearing accessory minerals (SEM imagery of baddeleyites and zircons analyzed with SIMS and used in the age calculations). Some baddeleyite grains (white in the back-scattered electron images) display minor alteration to zircon (light-gray rims; e.g., spot F755).

View this table:
  • View inline
  • View popup
Table S4.

ID-TIMS U-Pb baddeleyite isotopic data for dolerite samples NL-13c and G02-B

View this table:
  • View inline
  • View popup
Table S5.

Baddeleyite and zircon in situ SIMS U-Pb isotopic data for samples TGS-05 and OLL-2

A structurally corrected virtual geomagnetic pole (VGP) for the dated sample TGS-05 is located at −0.7° N and 288.2° E, whereas the VGP of sample MDK-05 is located at −0.9° N and 283.7° E (SI Methods, Paleomagnetism, Fig. S3, and Table S6). Samples OLL-2, G02-B, and NL-13c already have published shallow magnetization directions similar to those reported for the Ongeluk Formation (6, 26, 27).

Fig. S3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S3.

(A, Inset) Representative Zijderveld diagrams for the demagnetization behavior of (a) TGS-05 and (c) MDK-05, with (b and d) corresponding equal area projections of the characteristic magnetic components. Equal area projections are in situ. In the Zijderveld diagrams, tick marks represent 1 μA⋅m−2. (A) The magnetic components on which the Ongeluk LIP key paleomagnetic pole is based. (B) The Ongeluk LIP key paleomagnetic pole. The Ongeluk LIP is defined with error margins (red circle). In addition, the Paleoproterozoic apparent polar wander path is presented with previously published paleomagnetic poles and error margins (blue circles) (Tables S6 and S7).

View this table:
  • View inline
  • View popup
Table S6.

Paleomagnetic data for TGS-05 and MDK-05, with a summary of published individual VGP data from the ca. 2,426 Ma Ongeluk Formation and associated intrusions, which were then combined to form a grand mean key paleomagnetic pole

Discussion

The Ongeluk LIP.

All samples were collected stratigraphically from within or below the lower Ongeluk Formation (Fig. 1) and represent either the feeder system of dolerite dikes and sills to the Ongeluk Formation basalts or coarse-grained interiors of thicker basalt flows. We interpret the Westerberg Sill Province and the N-trending dolerite dike swarm (Fig. 1), both intruding into the Griqualand West subbasin stratigraphy, as parts of the same short-lived magmatic event based on their temporal, spatial, and stratigraphic proximities and similar paleomagnetic results. The age of the Westerberg Sill itself was defined by an upper intercept date of 2,441 ± 6 Ma composed of five discordant baddeleyite analyses (28). However, excluding the most discordant analysis from this result yields a more probable upper intercept date of 2,428 ± 4 Ma. This reinterpretation is supported by concordant baddeleyite analyses from a second sill dated at 2,426 ± 1 Ma as part of the same study (28). Combining all of these dates, we calculate a weighted mean date of 2,425.5 ± 2.6 Ma (Fig. 2) as the age of a single relatively short-lived magmatic event, which is now clearly distinguished from ≤2,250–2,240 Ma Hekpoort Formation volcanism with which it was previously correlated. Because sample NL-13c is located ∼1,000 km to the east of the other sample sites, our results indicate a defined craton-scale LIP, the Ongeluk LIP.

The Ongeluk Key Paleomagnetic Pole.

The VGPs presented in this study overlap with previously published paleopoles for the Ongeluk Formation and associated intrusions (6, 26⇓–28). Collectively, the combined VGPs from all of these studies for the basalts from the Ongeluk Formation (6) and their intrusive feeders (26⇓–28) define a grand key paleomagnetic pole for the Ongeluk LIP that is near-equatorial (29) at 4.1° N, 282.9° E (Table S7), with an A95 of 5.3° that achieves five of seven on the quality scale by Van der Voo (30).

View this table:
  • View inline
  • View popup
Table S7.

Compilation of paleomagnetic poles for the Neoarchean-Paleoproterozoic Kaapvaal Craton

Revising the Transvaal Supergroup Stratigraphy.

The age for Ongeluk Formation volcanism demands the revision of stratigraphic correlations between the successions of the Griqualand West and Transvaal subbasins of the Transvaal Supergroup (Fig. 1). Previous studies have correlated the Postmasburg and Pretoria Groups using the 2,222 ± 13 Ma Ongeluk Formation and the ≤2,250–2,240 Ma Hekpoort Formation volcanic rocks (11, 12, 14, 18). This correlation is now shown to be incorrect based on the 2,426 ± 3 Ma age constraint provided by the Ongeluk Formation volcanic rocks and associated intrusions. In addition, lithologic and chemostratigraphic data for the Duitschland and Rooihoogte formations (15, 18, 23, 24, 31) coupled with recent age constraints (9, 11, 25) and arguments relying on basin architecture (16) indicate that the Duitschland and Rooihoogte formations are younger than the Postmasburg Group (Fig. 1). A more complete discussion on the Transvaal Supergroup stratigraphy and the proposed revisions is provided in the stratigraphic synthesis (SI Methods, Stratigraphic Synthesis).

Atmospheric Oxygen Oscillations.

This stratigraphic interpretation indicates a dynamic state of atmospheric oxygen levels during the early Paleoproterozoic glacial period. The onset of the GOE occurred in the immediate aftermath of the Makganyene Formation glaciation with deposition of the world’s largest manganese deposit, the Hotazel Formation, with a negative Ce anomaly, both indicative of oxygenation (7, 21). This oxygenation was followed by a return to anoxic atmospheric conditions as indicated by the lack of a Ce anomaly in foreslope carbonates of the Mooidraai Formation (Fig. 1) (21, 32) and a mass-independent fractionation of sulfur (MIF-S) signal recorded by early diagenetic sulfides from the lower Duitschland Formation (Fig. 1) (23). Subsequent oxygenation events occurred during deposition of the upper Duitschland Formation (23) and once again, in the middle of the Rooihoogte Formation (24) as indicated by the reappearance and disappearance of the MIF-S signal (Fig. 1).

Detrital pyrite grains persist above the oldest glacial diamictite in the Mississagi Formation in the Huronian Supergroup of the Superior Craton, Canada (33) and support oscillations in atmospheric oxygen herein inferred from the records of the Transvaal Supergroup (Fig. 3). Our geochronologic and stratigraphic framework argues against a simple, monotonic rise of atmospheric oxygen in the early Paleoproterozoic, a time period further characterized by four glaciations. Instead, the onset of the GOE was followed by oscillations in atmospheric oxygen content across the 10−5 present atmospheric level (PAL) threshold over an ∼200-My interval, adding empirical evidence to atmospheric modeling predictions (4).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

A is a graph illustrating the approximate chronology (Table S1) of the glaciations and atmospheric oxygen oscillations according to the related redox indicators in the Huronian and Transvaal basins (Table S2) using the same symbols as used in Fig. 1. Also shown are dated mafic and felsic magmatic events as well as δ13C ranges for carbonates (Tables S1 and S2) and the extent of stratigraphic records in each basin denoting gaps in records at unconformities and disconformities. (B) The early Paleoproterozoic geography of the Superior, Kola–Karelia, Hearne, and Wyoming cratons as integral parts of the supercraton Superia (5, 36)* with the addition of the Kaapvaal and Pilbara cratons in the supercraton Vaalbara configuration (52). The early Paleoproterozoic basins developed on these cratonic fragments include both ca. 2.51–2.43 Ga volcanic rocks and glacial units, which can be correlated across the cratons. The glacial units in bold denote the glacial deposits likely recording the first glaciation. All of the cratonic fragments also contain dolerite dikes and sills emplaced between ca. 2.51 Ga and 2.43 Ga, showing the extent of the LIPs formed during this time. Available paleomagnetic studies indicate that the majority of the cratonic fragments (as part of supercraton Superia) were positioned near the paleoequator. The arrows denoting present-day true north in the crustal blocks illustrate the rotations necessary to make the reconstruction. (Inset) The hypothesized paleolatitude of these Archean cratons in the early Paleoproterozoic.

Linking Snowball Earth and the GOE.

The Makganyene Formation, deposited in the tropics (6), records the oldest known Snowball Earth event (7). Makganyene Formation diamictites are now constrained to be slightly older than ca. 2,426 Ma and likely correlate with glacial units of the Ramsay Lake Formation in the Huronian Supergroup, Canada and the Campbell Lake Formation from the Snowy Pass Supergroup on the Wyoming Craton in the United States. Other correlative glacial units worldwide include the ca. 2,435 Ma Polisarka Formation on the Kola–Karelia Craton in Fennoscandia (34) as well as possibly, the Meteorite Bore Member on the Pilbara Craton in Australia, defining the wide extent of the oldest Paleoproterozoic glaciation (Fig. 3). This early Paleoproterozoic glaciation is broadly coeval with the onset of the GOE; both events are now tightly bracketed between ca. 2,460 Ma, the age of volcanic rocks near the base of the Huronian Supergroup (8, 10), and ca. 2,426 Ma, the approximate age of the Makganyene Formation glaciation and near-coeval Ongeluk Formation volcanic rocks in the Transvaal Supergroup. In addition, the 2,442 ± 2 Ma Seidorechka Formation and the overlying 2,435 ± 2 Ma Polisarka formations (34, 35) overlying the Kola–Karelia Craton may tightly bracket the oldest Paleoproterozoic glacial event between ca. 2,442 Ma and 2,435 Ma, respectively, implying that it lasted less than 7 My (Fig. 3).

Our results add to the growing evidence for large low-latitude continental landmasses in the early Paleoproterozoic, including the Kaapvaal Craton (6) and the clan of cratons that defines supercraton Superia (5, 36, 37): Superior, Wyoming, Hearne, and Karelia–Kola. These landmasses, at least in part contiguous, record a series of LIP events between 2,510 Ma and 2,440 Ma and include the Mistassini, Kaminak, Baltic, Baggot Rocks, and Matachewan LIPs (Fig. 3) (5) before the ca. 2,426 Ma Ongeluk LIP of the Kaapvaal Craton. Cumulatively, these large juvenile volcanic provinces on extensive low-latitude continental landmasses are likely to have triggered near-equatorial glaciations via enhanced chemical weathering of aerially extensive, nutrient-rich continental flood basalts. This weathering resulted in increased carbon dioxide drawdown (38, 39) and an enhanced flux of phosphorus and other essential nutrients (40) onto extensive continental margins and into intracratonic basins. An enhanced nutrient flux would have greatly increased photosynthetic activity and oxygen production, temporally linked to higher net burial of organic carbon in accumulating sediments, as reflected by δ13C values of carbonates in the upper Tongwane Formation of the Transvaal Supergroup (31). The Tongwane Formation is locally preserved above the ca. 2,480–2,460 Ma iron formations (41) but below the major unconformity documented in the entire Transvaal basin (Fig. 1) (18). Thus, even an incipient rise of free atmospheric oxygen would have led to rapid oxidation of atmospheric methane (42), forcing catastrophic climate change and plunging Earth into a global glaciation (43). Importantly, the dated near-equatorial Ongeluk LIP, conformably overlying and interfingering with the uppermost Makganyene Formation glacial diamictites (Fig. 1), illustrates the dual role of LIPs in these global events, in that they would also have contributed carbon dioxide to rebuilding the greenhouse atmosphere that led to abrupt termination of the first Snowball Earth state of the early Paleoproterozoic.

Comparing the Paleoproterozoic with the Neoproterozoic.

In the Neoproterozoic, a remarkably similar sequence of events occurred, involving successive emplacement of multiple LIPs on the supercontinent Rodinia, a low-latitude position of this supercontinent, and incipient rifting and breakup (44). The massive Franklin LIP at ca. 717 Ma (45) immediately preceded the most dramatic and longest global glaciation of the Neoproterozoic, the Sturtian (46). This overall period is characterized by the second most dramatic change in surface redox conditions linked with Snowball Earth glaciations (47) and accompanied high rates of organic carbon burial (48, 49). Although substantial differences between the Paleo- and Neoproterozoic glacial periods might be expected, for instance, in the triggering mechanisms for the initial global glaciations, because methane was probably a more important greenhouse gas before the Makganyene glaciation (42, 43) than before the Sturtian glaciation (38, 39), there are also uncanny parallels. Examples include supercraton- or supercontinent-size continental landmasses that were capped by continental flood basalts, incipient rifting and/or breakup, and rapid transit to low latitudes. All of these factors enhanced chemical weathering of juvenile basaltic material and greatly increased the flux of bio-limiting nutrients to depositional basins, thus leading to a biotic response of higher organic carbon burial (40, 50) as well as deposition of giant iron and manganese deposits. It seems likely that these similar scenarios are not coincidental but that the critical factors (assembly of large landmasses, LIPs, incipient rifting, and relief enhancement—all resulting in a lithospheric mass anomaly and movement to low latitudes, high rates of organic carbon burial, surface oxygenation, and Snowball Earth glaciations) are mechanistically linked. In this case, a critical link might be true polar wander caused by the lithospheric mass anomaly that nudged the basalt-covered and rifting supercontinental landmasses to the equator (51), where chemical weathering and nutrient fluxes kicked into high gear and triggered the biotic, redox, and climatic responses.

SI Methods

Sampling.

Samples analyzed and interpreted as part of this study were taken from dolerite intrusions (i.e., sills and dikes) that intruded into basement of the western part of the Kaapvaal Craton and the overlying cover succession of the Transvaal Supergroup in the Griqualand West subbasin, South Africa, except for one sample that was taken ∼1,000 km to the east on the southeastern part of the craton (Fig. 1, Fig. S1, and Table S3). However, one sample (TGS-05) is tentatively interpreted as a coarse-grained basalt from the base of the Ongeluk Formation (Table S3). The sampled targets were selected using 1:250,000 map sheets of the Council for Geoscience (TGS-05 and MDK-05) or studied previously (e.g., OLL-2, G02-B, and NL-13c) (6, 26, 27). The previously studied intrusions yielded “Ongeluk-like” paleomagnetic directions. In addition, samples M03WA and TGS-01 (28) were reinterpreted in this study.

Sample TGS-05 (Fig. 1, Fig. S1, and Table S3) was identified in a road cut between Danielskuil and Postmasburg. It is interpreted as a medium-grained green basalt at the very base of the Ongeluk Formation. It is tentatively interpreted as a lava flow based on amygdales and nearby flow-top breccias. The transition into the volcanic breccias was not observed, however, and a shallow dolerite sill cannot be ruled out. As part of our study, it was drilled for complimentary paleomagnetic data. The Ongeluk Formation dips less than 10° in this location. Paleomagnetic samples were also collected from dolerite sill site MDK-05 (Fig. 1, Fig. S1, and Table S3) on the Groenwater Farm between Danielskuil and Postmasburg (near sample site TGS-05). It is a medium-grained green dolerite at the contact of the Asbestos Hill Subgroup and the overlying Makganyene Formation. Contacts and structural attitude were not observed because of poor outcrop, but the geological map indicates a conformable sill that may feed up into the Ongeluk Formation basalts. Nearby units of the Ongeluk Formation dip less than 10°. Sample OLL-2 (Fig. 1, Fig. S1, and Table S3) is from a sill within the lower Ongeluk Formation volcanic rocks exposed along the road from Griekwastad to Campbell. It is a coarse-grained green dolerite occurring within the lower Ongeluk Formation and was paleomagnetically studied previously (6). It is interpreted as a sill, but contacts were not observed. The attitude is unknown. Sample G02-B (Fig. 1, Fig. S1, and Table S3) was taken along the Orange River at the Gewonne Farm, where an exposure of dolomite and limestone of Schmidtsdrif Subgroup mapped within the river bed reveals an ∼80-m-wide subvertical, N-trending dolerite dike with a previously determined Ongeluk-like paleomagnetic direction (26). It is a medium- to coarse-grained green dolerite. The sample NL-13c (Fig. 1, Fig. S1, and Table S3) was taken from a shallow-dipping, medium-grained gray dolerite sheet intruding into Mesoarchean basement below the Pongola Supergroup ∼1,000 km away in the southeastern Kaapvaal Craton. It was described as a sheet intrusion (27), having a structural orientation similar to the regional strike and dip of the nearby basal Pongola Supergroup units. It cross-cuts an east–north–east-trending dyke, which in turn, cross-cuts an south–east-trending dolerite dike in the region. This dolerite sheet had been interpreted previously as a ca. 0.18 Ga intrusion based on a structurally uncorrected ‘”Karoo-like” magnetic remanence (27); however, structural correction reveals the Ongeluk-like magnetic remanence. Sample M03WA (Fig. 1, Fig. S1, and Table S3) is from the Westerberg Sill itself in the Griqualand West subbasin and was sampled for geochronologic and paleomagnetic studies previously along with TGS-01 (28), yielding a 2,441 ± 6 Ma age and Ongeluk-like paleomagnetic direction. M03WA was taken near the old asbestos mine village of Westerberg along the Orange River. It is a medium-grained gray–green dolerite sill intruding the Asbestos Hills Subgroup along the contact between contrasting rock units. TGS-01 (28), a medium-grained gray–green dolerite sample (Fig. 1, Fig. S1, and Table S3), intruded the Asbestos Hill Subgroup much farther northeast into the hinterland of the Kaapvaal Craton nearer to Kuruman. It is probably a sill, but no contacts were seen because of poor exposure. Regional dip of bedding is less than 10°. It was dated at 2,426 ± 1 Ma and also, has an Ongeluk-like paleomagnetic direction (28).

Geochronology—ID-TIMS Analysis.

Recovery of baddeleyite was attempted from all of the mafic rock samples using a water-based separation technique from crushed material at Lund University. Baddeleyite is a common accessory phase in mafic intrusions and coarse-grained basalt flows and a proven U-Pb geochronometer. This technique was successful in separating baddeleyite from the coarser-grained mafic intrusions NL-13c and G02-B. Consequently, these samples were dated by dissolution using U-Pb ID-TIMS (53).

The best-quality baddeleyite grains were transferred to Teflon dissolution capsules and rinsed repeatedly in ultrapure 7 M HNO3 and H2O to decrease the Pb blank. An ultrapure HF-HNO3 mixture was then added to the dissolution capsules together with a 205Pb-233–236U tracer solution. Capsules were then placed in an oven at 190 °C for 3 d until the grains were fully dissolved. After dissolution, the samples were dried down on a hotplate and redissolved in ultrapure 6 M HCl with the addition of ultrapure 0.25 M H3PO4 before being dried down again on a hotplate. The samples were then dissolved in 2 μL silica gel and loaded on outgassed Re filaments. The U and Pb isotopic compositions were measured at the Department of Geosciences at the Swedish Museum of Natural History using a Thermo Finnigan Triton Mass Spectrometer equipped with Faraday detectors and a Secondary Electron Multiplier. Pb isotope analyses were measured in dynamic (peak-jumping) mode using the Secondary Electron Multiplier. After completion of the 204Pb, 205Pb, 206Pb, and 207Pb isotopic measurements (typically in the 1,210 °C to 1,230 °C filament temperature range) for between 50 and 100 cycles, the 233U, 236U, and 238U isotopic measurements were made in dynamic mode at filament temperatures greater than 1,350 °C. Procedural blanks are typically 1.0 pg Pb and 0.1 pg U at the Swedish Museum of Natural History. Mass fractionation is 0.1% per mass unit for Pb determined by replicate analyses of National Bureau of Standards (NBS) reference materials SRM 981 and SRM 983. U fractionation was determined directly from the measured 233U/236U isotopic ratio. Initial Pb compositions were taken from the modeled common Pb evolution curve (54). The decay constants used were 1.55125 × 10−11 (238U) and 9.8485 × 10−10 (235U), with the isotopic composition of U being 238U/235U = 137.88 (55, 56). The listed uncertainties in Pb/U ratios were calculated by propagating the within-run error for measured isotopic ratios with the uncertainties in fractionation (±0.04 for Pb; absolute uncertainties), Pb and U blank concentration (±50%), and Pb blank composition (2.0% for 206Pb/204Pb and 0.2% for 207Pb/204Pb). Analytical results were calculated and plotted using the Microsoft Excel Macro, Isoplot (57).

Four baddeleyite fractions were analyzed from NL-13c (Table S4), each composed of between five and eight moderately brown and minute (<40 μm in length) crystals. One fraction is concordant within error, whereas the other fractions plot between 2.5 and 9.5% discordant (Fig. S2). The upper and lower intercept dates are 2,423 ± 7 and 510 ± 180 Ma, respectively, with a mean square weighted deviation (MSWD) of 0.27. The date of 2,423 ± 7 Ma is interpreted to represent the crystallization age of this sample. For G02-B, five fractions were analyzed, each composed of between one and four grains (Table S4). The grains were dark brown and somewhat frosty because of the thin skins of secondary zircon, and results were between 4 and 12% discordant (Fig. S2). An upper intercept date is calculated at 2,421 ± 3 Ma (MSWD = 2.20), which we interpret as the crystallization age of this dyke. The lower intercept date of 132 ± 89 Ma indicates minor near-0 Ma age Pb loss or isotopic disturbance caused by ca. 0.18 Ga Karoo magmatism.

Geochronology—SIMS Analysis.

Physical separation failed for the two critical finer-grained samples, TGS-05 and OLL-2, in which baddeleyite and zircon grains were generally ≤25 μm in size. For the analysis of such grains, we applied in situ U-Pb isotopic measurements using SIMS (58). The first step in this in situ dating method of “microbaddeleyites” involves identifying baddeleyite grains (and minor zircon if present) in polished thin sections by locating Zr-bearing phases through imaging and mapping using back-scattered electron imaging with brightness triggers and follow-up energy-dispersive spectrometry on a scanning electron microscope (Fig. S2). Mapped accessory Zr phases (baddeleyite, zircon, and zirconolite) are then assessed and ranked based on their size, shape, content of inclusions, and degree of alteration (e.g., secondary zircon rims on baddeleyite). Selected regions of the polished thin section, hosting the best targets, are then cut out, mounted in an epoxy disk together with baddeleyite and zircon reference materials, and gold-coated for analysis. After gold coating, the mounts are reimaged in reflected light microscopy to assist in locating the target grains for SIMS analysis. U and Pb isotopic ratios were determined on in situ baddeleyite grains in thin sections using a 20-µm-diameter primary ion beam on the CAMECA ims1270 Mass Spectrometer at UCLA. The field aperture was adjusted to reduce the effective sampling regions to 4-µm areas to screen out ions from host phases. Analysis of grains as small as 3 μm is possible. Additional details of the operating procedures for the CAMECA 1270 at UCLA are given elsewhere (58). Beams were focused onto the Zr-bearing phases using a 20–25-μm-diameter spot, with the secondary ions extracted at 10 kV. The secondary Pb+ yield was enhanced by a factor of ∼10 through oxygen flooding (∼105 torr) into the sample chamber. Isotopes of 204Pb, 206Pb, 207Pb, and 208Pb were counted in conjunction with 232Th, 238U, 238UO, and 238UO2 using a Secondary Electron Multiplier for between 15 and 20 cycles, taking ∼15 min per spot. The instrumental fractionation of Pb/U was corrected using UO2+/U+ vs. Pb/U for baddeleyite of the Phalaborwa reference material (values from 12.4 to 4.6) and UO+/U+ for zircon using the AS3 reference material (values from 10.4 to 8.2) (59). U concentration was calculated from measured U/94Zr2O of zircon reference material 91,500 with a concentration of 80 ppm U (×1 in zircon and ×2 in baddeleyite). Pb isotopes are not fractionated in SIMS analysis (59). Common Pb was corrected using the 204Pb method (60). Radiogenic 206Pb was typically >95%. Analytical results were also calculated and plotted using the Microsoft Excel Macro, Isoplot (57).

For sample TGS-05, 14 spots on different baddeleyite grains were analyzed (Fig. S2 and Table S5). Six analyses were rejected from the final age calculation because of low radiogenic 206Pb (<98%) and also, unacceptably high Th/U (>0.6), which has been observed to correlate with alteration or contamination. Linear regression of the eight remaining baddeleyite analyses resulted in an upper intercept date of 2,424 ± 32 Ma (MSWD = 0.76), which we interpret as the emplacement age of the volcanic lava flow (Fig. S2). The baddeleyites of OLL-2 suffered more alteration than those of TGS-05 (Fig. S2 and Table S5), which was indicated by secondary zircon rims of variable width around the majority of the grains. Three of nine baddeleyite analyses were rejected because of low radiogenic 206Pb (<98%). Linear regression of the remaining six analyses resulted in an upper intercept date of 2,397 ± 22 Ma with an MSWD = 1.90 (Fig. S2). Analysis of one baddeleyite and one coexisting zircon yielded 207Pb/206Pb dates of 2,421 ± 22 and 2,404 ± 16 Ma, respectively (2σ). Data from two other zircons were extremely discordant and suggest significant Pb loss during more recent alteration and metamorphism. The date of 2,421 ± 22 Ma from a single well-preserved baddeleyite grain (F861a) is interpreted to reflect better the true emplacement age of the sill in light of the evidence for alteration and metamorphism.

Paleomagnetism.

Core specimens from dolerite bodies TGS-05 and MDK-05 were collected using a portable gas-powered drill and oriented using both magnetic and sun compasses. Structural orientation (attitude) was derived from bedding of the Ongeluk Formation volcanic rocks or the Asbestos Hills Subgroup in the vicinity. All measurements of magnetic remanence were made using the superconducting rock magnetometer at the UJ (a vertical 2G Enterprises DC-4K Magnetometer) equipped with an RAPID Consortium Automatic Sample Changer. All specimens were exposed to stepwise alternating-field demagnetization from 2 to 100 mT. The process was stopped when sample intensity dropped below instrument noise level (∼1 pAm2) or samples started to show erratic behavior. Magnetic components were quantified via Zijderveld diagrams and least squares component analysis (61) using the software Paleomag 3.1b2 (62). Linear fits were included in subsequent analyses if they had a mean angular deviation ≤10°. Paleopole calculations are based on the assumption of a geocentric axial dipole field and a stable Earth radius throughout geological time that equals the present day radius. Visualization of pole positions and paleogeographic reconstructions was made with GPlates (63).

A summary of the demagnetization results is provided (Fig. S3 and Table S6). Specimens from sample site TGS-05 display a low-coercivity magnetic component that was removed during the first couple of demagnetization steps from 7.5 to 20 mT. This component is directed north and upward, with six of eight specimens displaying this component [declination (Decl.) = 346.3°, inclination (Incl.) = −46.3°, α95 = 8.55, k = 51.98, n = 6]. At demagnetization steps from 20 to 60 mT, a consistent medium- to high-coercivity component was defined and oriented west and downward, with characteristic remanence directions decaying toward the origin on removal of erratic directions recorded above 60 mT. This measurement defines the characteristic magnetic component at Decl. = 266.7°, Incl. = 6.4°, α95 = 15.0, k = 21.6, and n = 5, which on structural correction, yields Decl. = 266.9° and Incl. = −8.3° (Fig. S3 and Table S6). A structurally corrected VGP of −0.7° N and 288.2° E (Fig. S3 and Table S6) is defined for this characteristic magnetization (dp = 7.6°, dm = 15.1°). The oval area of confidence around the paleopole has principal axes aligned parallel (dp) and perpendicular (dm) to the great circle connecting the paleopole to the sampling site. Samples from site MDK-05 displayed randomly oriented lower-coercivity magnetic components that were removed below 50 mT and either shallow upward and westerly or shallow downward and easterly directed high-coercivity components above 50 mT. The high-coercivity components were identified in all eight samples and yielded a combined VGP mean at Decl. = 263.2°, Incl. = −0.6°, α95 = 11.5, k = 21.3, and n = 8, which on structural correction, becomes Decl. = 264.6° and Incl. = −15.8° (Fig. S3 and Table S6). A structurally corrected VGP of −0.9° N and 283.7° E is revealed (dp = 6.1°, dm = 11.8°) (Fig. S3 and Table S6).

A paleomagnetic pole was calculated from combining 27 tilt-corrected site VGPs of ca. 2,426 Ma intrusive and extrusive units of the Ongeluk LIP, Kaapvaal Craton (Fig. S3 and Table S6), which represent ∼20 independent cooling units. Unit component means with α95 values larger than 20°, and k values smaller than 10 were excluded from the pole calculation. This pole calculation includes the dated site OLL-2 as well as G02-B, TKW-A (M03WA), FYL (TGS-01), and NL-13c (6, 26⇓–28). The resulting paleopole is located at 4.1° N and 282.9° E (with an antipole at −4.1° N and 103.9° E) with A95 = 5.3° and K = 28.1, indicating a paleolatitude of 11° ± 6° (Table S7 and Fig. S3). This pole rates as five of seven possible in terms of quality criteria by Van der Voo (30), only missing criterions 7 (“no resemblance to paleopoles of younger age”) and 6 (“the presence of reversals”) (30). The ca. 2,426 Ma paleomagnetic pole resembles the paleomagnetic pole of a ca. 2.00 Ga paleoweathering surface that is exposed in the Griqualand West subbasin (64). Site OLL-2 did display an opposite polarity compared with all other sites but was not included in the paleopole calculation. The primary nature of the paleomagnetic pole is supported by a positive baked contact test (26) and positive conglomerate test (6). Our ca. 2,426 Ma paleomagnetic pole satisfies the more stringent criteria to be recognized as a key paleomagnetic pole (29).

Stratigraphic Synthesis.

Our samples and U-Pb dating results, together with all other relevant details, are summarized in Fig. 1, which represents our synthesis of the complex stratigraphy of the Transvaal Supergroup in its two main structural subbasins, the Griqualand West subbasin in the southwest and the Transvaal subbasin in the northeast (Fig. S1). These basin remnants occur on either side of a positive basement feature, the Vryburg Arch (16, 17). A third erosional remnant across the border in Botswana, the Kanye subbasin, can be considered as the western continuation of the Transvaal subbasin and is not further discussed.

The stratigraphy as shown in Fig. 1 generally follows the South African Committee on Stratigraphy and other reviews (9, 12) but with some important modifications. Our U-Pb age for the Ongeluk magmatism, at 2,426 ± 3 Ma, unlocks the long-held Ongeluk Formation basalts (Griqualand West subbasin) to Hekpoort Formation basalts (Transvaal subbasin) correlation (11, 12, 14, 18) and eliminates the ∼200-My hiatus at the base of the Postmasburg Group. With no major time hiatus at the base of the Makganyene Formation to contend with, we recognize that the transition from the Koegas Subgroup, characterized by alternating Mn-rich iron formations and sandstones (16, 17, 65), to the diamictites of the Makganyene Formation is essentially conformable (17) over distances of 50–200 km. This observation precludes significant deformation and basin inversion having occurred before deposition of the lower Postmasburg Group. Erosional effects and very low-angle disconformable relationships where Makganyene Formation diamictites directly overlie the upper Asbestos Hills Subgroup iron formations on the Ghaap Plateau (18, 66) can be largely attributed to a major sea-level fall associated with the low latitude and therefore, global (6, 7) Makganyene Formation glaciation. Basin upward shallowing started during final deposition of the iron formations. Less than 1° warping of the platform edge, perhaps associated with minor activity on syn-sedimentary basin-margin faults (e.g., the Griquatown fault system) (17, 65, 66), could also have contributed to the gradual disappearance of the Koegas Subgroup to the northeast. We thus consider the base of the Makganyene Formation to be a disconformity largely attributed to the glaciation (Fig. 1). Interestingly, the Tongwane Formation, which has a similar stratigraphic position with respect to the underlying iron formations as the Koegas Subgroup, also contains moderate Mn enrichment (41).

With no major time hiatus and considering the new ages that have emerged in recent years (15, 19, 20), a next logical step is to consider the Makganyene Formation diamictites and conformably overlying basalts of the Ongeluk Formation as a significant but temporary interruption of the carbonate- and iron formation-dominated deposition of the entire Griqualand West subbasin stratigraphy (17). After these climatic and magmatic perturbations decayed, Mn-rich iron formation deposition resumed with the Hotazel Formation (7, 17), continuing a theme first established in the underlying Koegas Subgroup. Deposition then transitioned once again to a carbonate platform, the foreslope Mooidraai Formation (15, 21, 32). Lack of compelling litho- and chronostratigraphic links of any of the Postmasburg Group formations with the upper Chuniespoort Group of the Transvaal subbasin strongly suggests that the entire Postmasburg Group is younger than the upper Chuniespoort Group (i.e., Penge and Tongwane formations) and simply has no preserved time-equivalent record in the Transvaal subbasin.

Sometime after deposition of the Mooidraai Formation carbonates, a significant tectonic disturbance led to the demise of the carbonate and iron formation platform in the Griqualand West subbasin, resulting in folding, uplift, and erosion. When deposition resumed, to the extent that it is preserved, the depocenter had shifted to the northeastern Transvaal subbasin, where the Duitschland Formation lies unconformably on folded iron formations of the Penge Formation and locally preserved iron-rich carbonates of the Tongwane Formation (31, 41), in a fundamentally different (successor?) basin with a different depocenter (18). Therefore, we consider the Duitschland Formation and the lower Pretoria Group (starting with the Rooihoogte Formation) of the Transvaal subbasin to be entirely younger than the depositional record of the Griqualand West subbasin as displayed in Fig.1.

This view is further supported by the different depositional and lithological character of the Duitschland Formation and lower Pretoria Group compared with strata of the Ghaap and Postmasburg groups and contrasting carbon isotope values for carbonates in these units (31, 67). Thick submarine volcanic rocks of the Ongeluk Formation overlying the Makganyene Formation diamictites have no equivalent in the Duitschland Formation above its basal diamictite, although a dolerite sheet of Ongeluk age is now recognized in the basement of the southeasternmost Kaapvaal Craton in this study; also, there is no lithological equivalent in the Duitschland Formation of the Hotazel Formation Mn- and Fe-rich deposits. Cap carbonates of the lower Duitschland Formation, with negative carbon isotope values, and carbonates of the upper Duitschland Formation, with positive carbon isotope values (31), also have no equivalents in the Postmasburg Group. Vice versa, carbonates of the Mooidraai Formation with near-zero carbon isotope values recording seawater composition (15) lack an equivalent in the Duitschland Formation. Our overall interpretation is fully compatible with all of the existing isotopic age constraints (Fig. 1 and Tables S1 and S2).

The Rooihoogte Formation has been considered by some authors (2, 11, 24) to be a time equivalent of the Duitschland Formation, with the Duitschland Formation restricted to the northeastern part of the Transvaal subbasin and the Rooihoogte Formation developed to the southwest. However, this interpretation is not universally adopted (18). There are important sedimentologic and provenance differences between the Duitschland and Rooihoogte formations that preclude their correlation. Although the Rooihoogte Formation is thin (up to ∼220-m thick and increases in thickness to the southwest) and has a gradational contact with the Timeball Hill Formation in the southwestern part of the Transvaal subbasin, the Duitschland Formation is up to ∼1,000 m in thickness and has a thick chert breccia on the contact with the Timeball Hill Formation, indicating exposure and an unconformity (18). In contrast, thickness of the Timeball Hill Formation does not change significantly along strike (18). There is also a dramatic change in provenance from the Duitschland Formation to the Timeball Hill Formation, with some Archean detrital zircon distribution modes disappearing, although early Paleoproterozoic sources appear for the first time in the provenance (11). The Rooihoogte Formation is, therefore, genetically linked with the Timeball Hill Formation, whereas the Duitschland Formation is unconformably overlain by the Timeball Hill Formation. The Rooihoogte and Timeball Hill formations were deposited in a different basin after a period of erosion that removed most of the Duitschland Formation and sourced from a distinctly different provenance. It is, therefore, reasonable to infer that the Rooihoogte and Duitschland formations are not correlative, with the Duitschland Formation being older.

We display the stratigraphy of the Transvaal subbasin up to the Dwaalheuvel Formation, which has the first widely accepted terrestrial red-bed sandstones in southern Africa, deposited in a fluvial to fluviodeltaic setting (12, 13). These red beds overlie oxic paleosols developed on the upper contact of the subaerial Hekpoort basalts (13) that are younger than ca. 2,250–2,240 Ma based on the age of a small number of detrital zircons in intercalated volcanoclastic sandstone (11) and the 2,256 ± 6 Ma U-Pb zircon age of a tuff in the underlying upper Timeball Hill Formation (9). Well-developed hematitic ooliths and pisoliths in fluviodeltaic deposits of the middle Timeball Hill Formation (68) constrain a significant rise in atmospheric oxygen level to have occurred well before extrusion of the Hekpoort Formation basalts and development of its overlying paleosol (12).

Fig. 1 also summarizes the record of mass-independent sulfur isotope fractionations being either present or absent (MIF-S vs. no MIF-S, respectively) as currently available in the literature (23, 24, 69⇓–71). Combining this record with that of other redox indicators (Table S2), we conclude that the onset of the GOE (1) was approximately coeval with the oldest Paleoproterozoic glaciation and the submarine extrusion of the up to ∼800-m thick Ongeluk Formation basalts (i.e., ca. 2,426 Ma) and postdates the last robust record of well-rounded detrital pyrite and rare uraninite grains in the Koegas Subgroup (22). Interestingly and puzzlingly, the upper Koegas Subgroup might also contain the potentially oldest red beds with detrital grains coated by hematite (65), although this interpretation remains to be further tested. Given the notable recurrence of an MIF-S signature in the lower Duitschland and Rooihoogte formations (23, 24), we conclude that this record resolves at least two oscillations back through the threshold of 10−5 PAL of oxygen (42). Lack of a Ce anomaly in carbonates of the upper Mooidraai Formation, following a negative anomaly in Mn-rich deposits of the underlying Hotazel Formation (21), is also inconsistent with a simple monotonic rise in atmospheric oxygen levels after the initial onset of GOE.

If our stratigraphic interpretations are correct, the onset of GOE is now geochronologically well-constrained to have occurred between ca. 2,460 Ma (72) and our precise U-Pb age of 2,426 ± 3 Ma for the Ongeluk Formation in the Transvaal Supergroup. Adding constraints from the Huronian Supergroup in Canada, where detrital pyrite grains persist well above the now well-dated ca. 2,460–2,455 Ma lower Huronian volcanic rocks (8, 10), tighten the age constraint closer to our Ongeluk Formation age. The tempo and rhythm of the GOE were likely complex, perhaps driven by various nonlinear processes, with current data resolving at least two major oscillations of atmospheric oxygen levels through the MIF-S threshold. Global interbasinal correlations based largely on the assumption of a single MIF-S to no MIF-S transition (2) might be erroneous in the absence of independent stratigraphic and geochronologic constraints.

Methods

Sampling.

Samples in this study were taken from dolerite intrusions that intruded into the basement of the Kaapvaal Craton and the overlying cover succession of the Transvaal Supergroup in the Griqualand West subbasin, South Africa, except for one sample taken ∼1,000 km to the east (Fig. 1, SI Methods, Fig. S1, and Table S3). One sample is tentatively interpreted as a coarse-grained basalt from the base of the Ongeluk Formation.

Geochronology—ID-TIMS Analysis.

Water-based separation of baddeleyite was attempted for all rock samples at Lund University. Samples NL-13c and G02-B yielded baddeleyite grains using this method. Grains were selected and analyzed using U-Pb ID-TIMS on a Thermo Finnigan Triton Mass Spectrometer at the Department of Geosciences at the Swedish Museum of Natural History. Additional details of the operating procedures for the ID-TIMS analyses at the Swedish Museum of Natural History and results are given in SI Methods, Fig. S2, and Table S4.

Geochronology—SIMS Analysis.

Zr-bearing phases were imaged, mapped using SEM in conjunction with energy-dispersive spectrometry in samples TGS-05 and OLL-2, and then, ranked (SI Methods and Fig. S2). U-Pb isotopic data using SIMS were determined in situ on the microbaddeleyite grains in the mapped thin sections using the CAMECA ims1270 Mass Spectrometer at the University of California, Los Angeles (UCLA). Additional details of the operating procedures for the SIMS at UCLA and results are given in SI Methods, Fig. S2, and Table S5.

Paleomagnetism.

Measurements of magnetic remanence from samples TGS-05 and MDK-05 were made using the superconducting vertical 2G Enterprises DC-4K Rock Magnetometer at the University of Johannesburg (UJ). All specimens were exposed to stepwise alternating-field demagnetization. Additional details of the sampling and operating procedures for the magnetometer at the UJ and results are given in SI Methods, Fig. S3, and Table S6.

Acknowledgments

We thank the staff of the SIMS and the TIMS laboratories at UCLA and the Department of Geosciences, Swedish Museum of Natural History, respectively, for all of the assistance during analyses. We thank N. Beukes for providing sample TGS-05. Funding was provided by grants from the Royal Physiographic Society in Lund (A.P.G.) and the Swedish Research Council (U.S.); M.O.d.K. and A.B. acknowledge support from the South African Department of Science and Technology and the National Research Foundation (DST-NRF)–funded Centre of Excellence for Integrated Mineral and Energy Resource Analysis (CIMERA). This article is a contribution to International Geoscience Programme (IGCP) 648: Supercontinents and Global Geodynamics.

Footnotes

  • ↵1To whom correspondence should be addressed. Email: ashley.gumsley{at}geol.lu.se.
  • Author contributions: A.P.G., U.S., and M.O.d.K. designed research; A.P.G., K.R.C., W.B., U.S., M.O.d.K., and E.R.L. performed research; K.R.C., W.B., U.S., and M.O.d.K. contributed new reagents/analytic tools; A.P.G., K.R.C., W.B., U.S., M.O.d.K., E.R.L., and A.B. analyzed data; and A.P.G., K.R.C., W.B., U.S., M.O.d.K., E.R.L., and A.B. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • ↵*Bleeker W, Chamberlain K, Kamo S, Kilian T, Buchan K (2016) Kaapvaal, Superior and Wyoming: Nearest neighbours in supercraton Superia. Proceedings of the 35th International Geological Congress.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1608824114/-/DCSupplemental.

Freely available online through the PNAS open access option.

References

  1. ↵
    1. Holland HD
    (2002) Volcanic gases, black smokers, and the Great Oxidation Event. Geochim Cosmochim Acta 66(21):3811–3826.
    .
    OpenUrlCrossRef
  2. ↵
    1. Hoffman PF
    (2013) The Great Oxidation and a Siderian snowball Earth: MIF-S based correlation of Paleoproterozoic glacial epochs. Chem Geol 362:143–156.
    .
    OpenUrlCrossRef
  3. ↵
    1. Goldblatt C,
    2. Lenton TM,
    3. Watson AJ
    (2006) Bistability of atmospheric oxygen and the Great Oxidation. Nature 443(7112):683–686.
    .
    OpenUrlCrossRefPubMed
  4. ↵
    1. Daines SJ,
    2. Lenton TM
    (2016) The effect of widespread early aerobic marine ecosystems on methane cycling and the Great Oxidation. Earth Planet Sci Lett 434:42–51.
    .
    OpenUrl
  5. ↵
    1. Ernst R,
    2. Bleeker W
    (2010) Large igneous provinces (LIPs), giant dyke swarms, and mantle plumes: Significance for breakup events within Canada and adjacent regions from 2.5 Ga to the Present. Can J Earth Sci 47(5):695–739.
    .
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Evans DA,
    2. Beukes NJ,
    3. Kirschvink JL
    (1997) Low-latitude glaciation in the Palaeoproterozoic era. Nature 386(6622):262–266.
    .
    OpenUrlCrossRef
  7. ↵
    1. Kirschvink JL, et al.
    (2000) Paleoproterozoic snowball earth: Extreme climatic and geochemical global change and its biological consequences. Proc Natl Acad Sci USA 97(4):1400–1405.
    .
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Ketchum KY,
    2. Heaman LM,
    3. Bennett G,
    4. Hughes DJ
    (2013) Age, petrogenesis and tectonic setting of the Thessalon volcanic rocks, Huronian Supergroup, Canada. Precambrian Res 233:144–172.
    .
    OpenUrlCrossRef
  9. ↵
    1. Rasmussen B,
    2. Bekker A,
    3. Fletcher IR
    (2013) Correlation of Paleoproterozoic glaciations based on U-Pb zircon ages for tuff beds in the Transvaal and Huronian Supergroups. Earth Planet Sci Lett 382:173–180.
    .
    OpenUrl
  10. ↵
    1. Bleeker W,
    2. Kamo SL,
    3. Ames DE,
    4. Davis D
    (2015) New field observations and U-Pb ages in the Sudbury area: Toward a detailed cross-section through the deformed Sudbury Structure. Geol Surv Canada Open File 7856:151–166.
    .
    OpenUrl
  11. ↵
    1. Schröder S,
    2. Beukes NJ,
    3. Armstrong RA
    (2016) Detrital zircon constraints on the tectonostratigraphy of the Paleoproterozoic Pretoria Group, South Africa. Precambrian Res 278:362–393.
    .
    OpenUrl
  12. ↵
    1. Beukes NJ,
    2. Dorland H,
    3. Gutzmer J,
    4. Nedachi M,
    5. Ohmoto H
    (2002) Tropical laterites, life on land, and the history of atmospheric oxygen in the Paleoproterozoic. Geology 30(6):491–494.
    .
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Eriksson PG,
    2. Cheney ES
    (1992) Evidence for the transition to an oxygen-rich atmosphere during the evolution of red beds in the lower Proterozoic sequences of southern Africa. Precambrian Res 54(2-4):257–269.
    .
    OpenUrlCrossRef
  14. ↵
    1. Cornell DH,
    2. Schütte SS,
    3. Eglington BL
    (1996) The Ongeluk basaltic andesite formation in Griqualand West, South Africa: Submarine alteration in a 2222 Ma Proterozoic sea. Precambrian Res 79(1-2):101–123.
    .
    OpenUrlCrossRef
  15. ↵
    1. Bau M,
    2. Romer RL,
    3. Lüders V,
    4. Beukes NJ
    (1999) Pb, O, and C isotopes in silicified Mooidraai dolomite (Transvaal Supergroup, South Africa): Implications for the composition of Paleoproterozoic seawater and “dating” the increase of oxygen in the Precambrian atmosphere. Earth Planet Sci Lett 174(1-2):43–57.
    .
    OpenUrl
  16. ↵
    1. Moore JM,
    2. Tsikos H,
    3. Polteau S
    (2001) Deconstructing the Transvaal Supergroup, South Africa: Implications for Palaeoproterozoic palaeoclimate models. J Afr Earth Sci 33(3-4):437–444.
    .
    OpenUrl
  17. ↵
    1. Polteau S,
    2. Moore JM,
    3. Tsikos H
    (2006) The geology and geochemistry of the Palaeoproterozoic Makganyene diamictite. Precambrian Res 148(3-4):257–274.
    .
    OpenUrl
  18. ↵
    1. Johnson MJ,
    2. Anhaeusser CR,
    3. Thomas RJ
    1. Eriksson PG,
    2. Altermann W,
    3. Hartzer FJ
    (2006) The Transvaal Supergroup. The Geology of South Africa, eds Johnson MJ, Anhaeusser CR, Thomas RJ (Geological Society of South Africa/Council for Geoscience, Johannesburg), pp 237–260.
    .
  19. ↵
    1. Moore JM,
    2. Polteau S,
    3. Armstrong RA,
    4. Corfu F,
    5. Tsikos H
    (2012) The age and correlation of the Postmasburg Group, southern Africa: Constraints from detrital zircon grains. J Afr Earth Sci 64:9–19.
    .
    OpenUrl
  20. ↵
    1. Fairey B,
    2. Tsikos H,
    3. Corfu F,
    4. Polteau S
    (2013) U-Pb systematics in carbonates of the Postmasburg Group, Transvaal Supergroup, South Africa: Primary versus metasomatic controls. Precambrian Res 231:194–205.
    .
    OpenUrlCrossRef
  21. ↵
    1. Bau M,
    2. Alexander B
    (2006) Preservation of primary REE patterns without Ce anomaly during dolomitization of Mid-Paleoproterozoic limestone and the potential re-establishment of marine anoxia immediately after the “Great Oxidation Event.” South Afr J Geol 109(1-2):81–86.
    .
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Johnson JE,
    2. Gerpheide A,
    3. Lamb MP,
    4. Fischer WW
    (2014) O2 constraints from Paleoproterozoic detrital pyrite and uraninite. Geol Soc Am Bull 126(5-6):813–830.
    .
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Guo Q, et al.
    (2009) Reconstructing Earth’s surface oxidation across the Archean-Proterozoic transition. Geology 37(5):399–402.
    .
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Luo G, et al.
    (2016) Rapid oxygenation of Earth’s atmosphere 2.33 billion years ago. Sci Adv 2(5):e1600134.
    .
    OpenUrlFREE Full Text
  25. ↵
    1. Hannah JL,
    2. Bekker A,
    3. Stein HJ,
    4. Markey RJ,
    5. Holland HD
    (2004) Primitive Os and 2316 Ma age for marine shale: Implications for Paleoproterozoic glacial events and the rise of atmospheric oxygen. Earth Planet Sci Lett 225:43–52.
    .
    OpenUrl
  26. ↵
    1. de Kock MO, et al.
    (2009) Paleomagnetism of a Neoarchean-Paleoproterozoic carbonate ramp and carbonate platform succession (Transvaal Supergroup) from surface outcrop and drill core, Griqualand West region, South Africa. Precambrian Res 169(1-2):80–99.
    .
    OpenUrlCrossRef
  27. ↵
    1. Lubnina N,
    2. Ernst R,
    3. Klausen M,
    4. Söderlund U
    (2010) Paleomagnetic study of NeoArchean-Paleoproterozoic dykes in the Kaapvaal Craton. Precambrian Res 183(3):523–552.
    .
    OpenUrlCrossRef
  28. ↵
    1. Kampmann TC,
    2. Gumsley AP,
    3. de Kock MO,
    4. Söderlund U
    (2015) U-Pb geochronology and paleomagnetism of the Westerberg Sill Suite, Kaapvaal Craton - Support for a coherent Kaapvaal-Pilbara Block (Vaalbara) into the Paleoproterozoic? Precambrian Res 269:58–72.
    .
    OpenUrl
  29. ↵
    1. Buchan KL
    (2013) Key paleomagnetic poles and their use in Proterozoic continent and supercontinent reconstructions: A review. Precambrian Res 238:93–110.
    .
    OpenUrlCrossRef
  30. ↵
    1. Van der Voo R
    (1990) The reliability of paleomagnetic data. Tectonophysics 184(1):1–9.
    .
    OpenUrlCrossRef
  31. ↵
    1. Bekker A, et al.
    (2001) Chemostratigraphy of the Paleoproterozoic Duitschland Formation, South Africa: Implications for coupled climate change and carbon cycling. Am J Sci 301(3):261–285.
    .
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Kunzmann M,
    2. Gutzmer J,
    3. Beukes NJ,
    4. Halverson GP
    (2014) Depositional environment and lithostratigraphy of the Paleoproterozoic Mooidraai Formation, Kalahari Manganese Field, South Africa. South Afr J Geol 117(2):173–192.
    .
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Ulrich T,
    2. Long DGF,
    3. Kamber BS,
    4. Whitehouse MJ
    (2011) In situ trace element and sulfur isotope analysis of pyrite in a Paleoproterozoic Gold Placer Deposit, Pardo and Clement Townships, Ontario, Canada. Econ Geol 106(4):667–686.
    .
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Brasier AT, et al.
    (2013) Earth’s earliest global glaciation? Carbonate geochemistry and geochronology of the Polisarka Sedimentary Formation, Kola Peninsula, Russia. Precambrian Res 235:278–294.
    .
    OpenUrl
  35. ↵
    1. Amelin YV,
    2. Heaman LM,
    3. Semenov VS
    (1995) U-Pb geochronology of layered mafic intrusions in the eastern Baltic Shield: Implications for the timing and duration of Paleoproterozoic continental rifting. Precambrian Res 75(1-2):31–46.
    .
    OpenUrlCrossRef
  36. ↵
    1. Hanski E,
    2. Mertanen S,
    3. Rämö T,
    4. Vuollo J
    1. Bleeker W,
    2. Ernst R
    (2006) Short-lived mantle generated magmatic events and their dyke swarms: The key unlocking Earth’s paleogeographic record back to 2.6 Ga. Dyke Swarms—Time Markers of Crustal Evolution, eds Hanski E, Mertanen S, Rämö T, Vuollo J (Taylor and Francis, London), pp 3–26.
    .
  37. ↵
    1. Bleeker W
    (2003) The late Archean record: A puzzle in ca. 35 pieces. Lithos 71(2-4):99–134.
    .
    OpenUrlCrossRef
  38. ↵
    1. Marshall HG,
    2. Walker JCG,
    3. Kuhn WR
    (1988) Long-term climate change and the geochemical cycle of carbon. J Geophys Res Atmos 93(D1):791–801.
    .
    OpenUrl
  39. ↵
    1. Donnadieu Y,
    2. Goddéris Y,
    3. Ramstein G,
    4. Nédélec A,
    5. Meert J
    (2004) A ‘snowball Earth’ climate triggered by continental break-up through changes in runoff. Nature 428(6980):303–306.
    .
    OpenUrlCrossRefPubMed
  40. ↵
    1. Horton F
    (2015) Did phosphorus derived from the weathering of large igneous provinces fertilize the Neoproterozoic ocean? Geochem Geophys Geosyst 16(6):1723–1738.
    .
    OpenUrl
  41. ↵
    1. Schröder S,
    2. Warke MR
    (2016) Termination of BIF deposition in the Paleoproterozoic: The Tongwane Formation, South Africa. South Afr J Geol 119(2):329–346.
    .
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Pavlov AA,
    2. Kasting JF
    (2002) Mass-independent fractionation of sulfur isotopes in Archean sediments: Strong evidence for an anoxic Archean atmosphere. Astrobiology 2(1):27–41.
    .
    OpenUrlCrossRefPubMed
  43. ↵
    1. Bekker A,
    2. Kaufman AJ
    (2007) Oxidative forcing of global climate change: A biogeochemical record across the oldest Paleoproterozoic ice age in North America. Earth Planet Sci Lett 258(3-4):486–499.
    .
    OpenUrl
  44. ↵
    1. Li ZX, et al.
    (2008) Assembly, configuration, and break-up history of Rodinia: A synthesis. Precambrian Res 160(1-2):179–210.
    .
    OpenUrlCrossRef
  45. ↵
    1. Macdonald FA, et al.
    (2010) Calibrating the Cryogenian. Science 327(5970):1241–1243.
    .
    OpenUrlAbstract/FREE Full Text
  46. ↵
    1. Rooney AD,
    2. Strauss JV,
    3. Brandon AD,
    4. Macdonald FA
    (2015) A Cryogenian chronology: Two long-lasting synchronous Neoproterozoic glaciations. Geology 43(5):459–462.
    .
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Pogge von Strandmann PAE, et al.
    (2015) Selenium isotope evidence for progressive oxidation of the Neoproterozoic biosphere. Nat Commun 6:10157.
    .
    OpenUrlCrossRefPubMed
  48. ↵
    1. Thomson D,
    2. Rainbird RH,
    3. Planavsky NJ,
    4. Lyons TW,
    5. Bekker A
    (2015) Chemostratigraphy of the Shaler Supergroup, Victoria Island, NW Canada: A record of ocean composition prior to the Cryogenian glaciations. Precambrian Res 263:232–245.
    .
    OpenUrlCrossRef
  49. ↵
    1. Turner EC,
    2. Bekker A
    (2016) Thick sulfate evaporite accumulations marking a mid-Neoproterozoic oxygenation event (Ten Stone Formation, Northwest Territories, Canada). Geol Soc Am Bull 128(1-2):203–222.
    .
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Cox GM, et al.
    (2016) Continental flood basalt weathering as a trigger for Neoproterozoic Snowball Earth. Earth Planet Sci Lett 446:89–99.
    .
    OpenUrl
  51. ↵
    1. Li Z-X,
    2. Zhong S
    (2009) Supercontinent-superplume coupling, true polar wander and plume mobility: Plate dominance in whole-mantle tectonics. Phys Earth Planet Inter 176(3-4):143–156.
    .
    OpenUrl
  52. ↵
    1. de Kock MO,
    2. Evans DAD,
    3. Beukes NJ
    (2009) Validating the existence of Vaalbara in the Neoarchean. Precambrian Res 174(1-2):145–154.
    .
    OpenUrlCrossRef
  53. ↵
    1. Gumsley A, et al.
    (2015) Precise U-Pb baddeleyite age dating of the Usushwana Complex, southern Africa: Implications for the Mesoarchaean magmatic and sedimentological evolution of the Pongola Supergroup, Kaapvaal Craton. Precambrian Res 267:174–185.
    .
    OpenUrl
  54. ↵
    1. Stacey JS,
    2. Kramers JD
    (1975) Approximation of terrestrial lead isotope evolution by a two-stage model. Earth Planet Sci Lett 26(2):207–221.
    .
    OpenUrl
  55. ↵
    1. Jaffey AH,
    2. Flynn KF,
    3. Glendenin LE,
    4. Bentley WC,
    5. Essling AM
    (1971) Precision measurement of half-lives and specific activities of U235 and U238. Phys Rev C Nucl Phys 4:1889–1906.
    .
    OpenUrlCrossRef
  56. ↵
    1. Steiger RH,
    2. Jäger E
    (1977) Subcommission on geochronology: Convention on the use of decay constants in geo- and cosmochronology. Earth Planet Sci Lett 36(3):359–362.
    .
    OpenUrl
  57. ↵
    1. Ludwig KR
    (2003) User’s Manual for Isoplot 3.00—A Geochronological Toolkit for Microsoft Excel (Berkeley Geochronology Center, Berkeley, CA), Special Publ 4.
    .
  58. ↵
    1. Chamberlain KR, et al.
    (2010) In situ U-Pb SIMS (IN-SIMS) micro-baddeleyite dating of mafic rocks: Method with examples. Precambrian Res 183(3):379–387.
    .
    OpenUrlCrossRef
  59. ↵
    1. Compston W,
    2. Williams IS,
    3. Meyer C
    (1984) U-Pb geochronology of zircons form lunar Breccia 73217 using a sensitive high mass-resolution ion microprobe. J Geophys Res Solid Earth 89(S02):B525–B534.
    .
    OpenUrl
  60. ↵
    1. Ireland TR,
    2. Williams IS
    (2003) Considerations in zircon geochronology by SIMS. Rev Miner Geochem 53(1):215–241.
    .
    OpenUrlFREE Full Text
  61. ↵
    1. Kirschvink JL
    (1980) The least-squares line and plane and the analysis of paleomagnetic data. Geophys J Int 62(3):699–718.
    .
    OpenUrlAbstract/FREE Full Text
  62. ↵
    1. Jones CH
    (2002) User-driven integrated software lives: “Paleomag” paleomagnetics analysis on the Macintosh. Comput Geosci 28(10):1145–1151.
    .
    OpenUrlCrossRef
  63. ↵
    1. Williams SE,
    2. Müller RD,
    3. Landgrebe TCW,
    4. Whittaker JM
    (2012) An open-source software environment for visualizing and refining plate tectonic reconstructions using high-resolution geological and geophysical data sets. GSA Today 22(4):4–9.
    .
    OpenUrlCrossRef
  64. ↵
    1. Evans DAD,
    2. Beukes NJ,
    3. Kirschvink JL
    (2002) Paleomagnetism of a lateritic paleoweathering horizon and overlying Paleoproterozoic red beds from South Africa: Implications for the Kaapvaal apparent polar wander path and a confirmation of atmospheric oxygen enrichment. J Geophys Res Solid Earth 107(B12):EPM 2-1–EPM 2-22.
    .
    OpenUrl
  65. ↵
    1. Schröder S,
    2. Bedorf D,
    3. Beukes NJ,
    4. Gutzmer J
    (2011) From BIF to red beds: Sedimentology and sequence stratigraphy of the Paleoproterozoic Koegas Subgroup (South Africa). Sediment Geol 236(1-2):25–44.
    .
    OpenUrlCrossRef
  66. ↵
    1. Beukes NJ
    (1984) Sedimentology of the Kuruman and Griquatown iron-formations, Transvaal Supergroup, Griqualand West, South Africa. Precambrian Res 24(1):47–84.
    .
    OpenUrlCrossRef
  67. ↵
    1. Bekker A, et al.
    (2008) Fractionation between inorganic and organic carbon during the Lomagundi (2.22–2.1 Ga) carbon isotope excursion. Earth Planet Sci Lett 271(1-4):278–291.
    .
    OpenUrl
  68. ↵
    1. Coetzee LL,
    2. Beukes NJ,
    3. Gutzmer J,
    4. Kakegawa T
    (2006) Links of organic carbon cycling and burial to depositional depth gradients and establishment of a snowball Earth at 2.3 Ga. Evidence from the Timeball Hill Formation,Transvaal Supergroup, South Africa. South Afr J Geol 109(1-2):109–122.
    .
    OpenUrlAbstract/FREE Full Text
  69. ↵
    1. Farquhar J,
    2. Bao H,
    3. Thiemens M
    (2000) Atmospheric influence of Earth’s earliest sulfur cycle. Science 289(5480):756–759.
    .
    OpenUrlAbstract/FREE Full Text
  70. ↵
    1. Bekker A, et al.
    (2004) Dating the rise of atmospheric oxygen. Nature 427(6970):117–120.
    .
    OpenUrlCrossRefPubMed
  71. ↵
    1. Johnson JE, et al.
    (2013) Manganese-oxidizing photosynthesis before the rise of cyanobacteria. Proc Natl Acad Sci USA 110(28):11238–11243.
    .
    OpenUrlAbstract/FREE Full Text
  72. ↵
    1. Pickard AL
    (2003) SHRIMP U-Pb zircon ages for the Palaeoproterozoic Kuruman Iron Formation, northern Cape Province, South Africa: Evidence for simultaneous BIF deposition on Kaapvaal and Pilbara cratons. Precambrian Res 125(3-4):275–315.
    .
    OpenUrlCrossRef
    1. Sumner DY,
    2. Bowring SA
    (1996) U-Pb geochronologic constraints on deposition of the Campbellrand Subgroup, Transvaal Supergroup, South Africa. Precambrian Res 79(1-2):25–35.
    .
    OpenUrlCrossRef
    1. Barton ES,
    2. Altermann W,
    3. Williams IS,
    4. Smith CB
    (1994) U-Pb zircon age for a tuff in the Campbell Group, Griqualand West sequence, South Africa: Implications for early Proterozoic rock accumulation rates. Geology 22(4):343–346.
    .
    OpenUrlAbstract/FREE Full Text
    1. Altermann W,
    2. Nelson DR
    (1998) Sedimentation rates, basin analysis and regional correlations of three Neoarchaean and Palaeoproterozoic sub-basins of the Kaapvaal craton as inferred from precise U-Pb zircon ages from volcaniclastic sediments. Sediment Geol 120(1-4):225–256.
    .
    OpenUrlCrossRef
    1. Nelson DR,
    2. Trendall AF,
    3. Altermann W
    (1999) Chronological correlations between the Pilbara and Kaapvaal cratons. Precambrian Res 97(3-4):165–189.
    .
    OpenUrlCrossRef
    1. Swart QD
    (1999) Carbonate Rocks of the Paleoproterozoic Pretoria and Postmasburg Groups, Transvaal Supergroup (Rand Afrikaans Univ, Johannesburg).
    .
    1. Bekker A,
    2. Karhu JA,
    3. Kaufman AJ
    (2006) Carbon isotope record for the onset of the Lomagundi carbon isotope excursion in the Great Lakes area, North America. Precambrian Res 148(1-2):145–180.
    .
    OpenUrlCrossRef
    1. Sekine Y, et al.
    (2011) Manganese enrichment in the Gowganda Formation of the Huronian Supergroup: A highly oxidizing shallow-marine environment after the last Huronian glaciation. Earth Planet Sci Lett 307(1-2):201–210.
    .
    OpenUrl
    1. Chandler FW
    (1988) Diagenesis of sabkha-related, sulphate nodules in the early Proterozoic Gordon Lake formation, Ontario, Canada. Carbonates Evaporites 3(1):75–94.
    .
    OpenUrl
    1. Fedo CM,
    2. Nesbitt HW,
    3. Young GM
    (1995) Unraveling the effects of potassium metasomatism in sedimentary rocks and paleosols, with implications for paleoweathering conditions and provenance. Geology 23(10):921–924.
    .
    OpenUrlAbstract/FREE Full Text
    1. Papineau D,
    2. Mojzsis SJ,
    3. Schmitt AK
    (2007) Multiple sulfur isotopes from Paleoproterozoic Huronian interglacial sediments and the rise of atmospheric oxygen. Earth Planet Sci Lett 255(1-2):188–212.
    .
    OpenUrl
    1. Bekker A,
    2. Kaufman AJ,
    3. Karhu JA,
    4. Eriksson KA
    (2005) Evidence for Paleoproterozoic cap carbonates in North America. Precambrian Res 137(3-4):167–206.
    .
    OpenUrlCrossRef
    1. Murakami T,
    2. Utsunomiya S,
    3. Imazu Y,
    4. Prasad N
    (2001) Direct evidence of late Archean to early Proterozoic anoxic atmosphere from a product of 2.5 Ga old weathering. Earth Planet Sci Lett 184(2):523–528.
    .
    OpenUrl
    1. Roscoe SM
    (1969) Huronian rocks and uraniferous conglomerates in the Canadian Shield. Geol Surv Canada Open File 68:1–205.
    .
    OpenUrl
    1. Gumsley A,
    2. Rådman J,
    3. Söderlund U,
    4. Klausen M
    (2016) U–Pb baddeleyite geochronology and geochemistry of the White Mfolozi Dyke Swarm: Unravelling the complexities of 2.70–2.66 Ga dyke swarms across the eastern Kaapvaal Craton, South Africa. GFF 138(1):115–132.
    .
    OpenUrl
    1. Reischmann T
    (1995) Precise U/Pb age determination with baddeleyite (ZrO2), a case study from the Phalaborwa igneous complex, South Africa. South Afr J Geol 98(1):1–14.
    .
    OpenUrlAbstract
    1. Letts S,
    2. Torsvik TH,
    3. Webb SJ,
    4. Ashwal LD
    (2011) New Palaeoproterozoic palaeomagnetic data from the Kaapvaal Craton, South Africa. Geol Soc Lond Spec Publ 357:9–26.
    .
    OpenUrl
    1. Zeh A,
    2. Ovtcharova M,
    3. Wilson AH,
    4. Schaltegger U
    (2015) The Bushveld Complex was emplaced and cooled in less than one million years—results of zirconology, and geotectonic implications. Earth Planet Sci Lett 418:103–114.
    .
    OpenUrl
    1. Letts S,
    2. Torsvik TH,
    3. Webb SJ,
    4. Ashwal LD
    (2009) Palaeomagnetism of the 2054 Ma Bushveld Complex (South Africa): Implications for emplacement and cooling. Geophys J Int 179(2):850–872.
    .
    OpenUrlAbstract/FREE Full Text
    1. Dorland HC,
    2. Beukes NJ,
    3. Gutzmer J,
    4. Evans DAD,
    5. Armstrong RA
    (2006) Precise SHRIMP U-Pb zircon age constraints on the lower Waterberg and Soutpansberg Groups, South Africa. South Afr J Geol 109(1-2):139–156.
    .
    OpenUrlAbstract/FREE Full Text
    1. de Kock MO,
    2. Evans DAD,
    3. Dorland HC,
    4. Beukes NJ,
    5. Gutzmer J
    (2006) Paleomagnetism of the lower two unconformity-bounded sequences of the Waterberg Group, South Africa: Towards a better-defined apparent polar wander path for the Paleoproterozoic Kaapvaal Craton. South Afr J Geol 109(1-2):157–182.
    .
    OpenUrlAbstract/FREE Full Text
    1. Kamo SL,
    2. Reimold WU,
    3. Krogh TE,
    4. Colliston WP
    (1996) A 2.023 Ga age for the Vredefort impact event and a first report of shock metamorphosed zircons in pseudotachylitic breccias and Granophyre. Earth Planet Sci Lett 144(3-4):369–387.
    .
    OpenUrl
    1. Carporzen L,
    2. Gilder SA,
    3. Hart RJ
    (2005) Palaeomagnetism of the Vredefort meteorite crater and implications for craters on Mars. Nature 435(7039):198–201.
    .
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Timing and tempo of the Great Oxidation Event
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Timing and tempo of the Great Oxidation Event
Ashley P. Gumsley, Kevin R. Chamberlain, Wouter Bleeker, Ulf Söderlund, Michiel O. de Kock, Emilie R. Larsson, Andrey Bekker
Proceedings of the National Academy of Sciences Feb 2017, 114 (8) 1811-1816; DOI: 10.1073/pnas.1608824114

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Timing and tempo of the Great Oxidation Event
Ashley P. Gumsley, Kevin R. Chamberlain, Wouter Bleeker, Ulf Söderlund, Michiel O. de Kock, Emilie R. Larsson, Andrey Bekker
Proceedings of the National Academy of Sciences Feb 2017, 114 (8) 1811-1816; DOI: 10.1073/pnas.1608824114
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Physical Sciences
  • Earth, Atmospheric, and Planetary Sciences
Proceedings of the National Academy of Sciences: 114 (8)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Transvaal Supergroup
    • Sampling
    • Results
    • Discussion
    • SI Methods
    • Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Model of the Amazon forest
News Feature: A sea in the Amazon
Did the Caribbean sweep into the western Amazon millions of years ago, shaping the region’s rich biodiversity?
Image credit: Tacio Cordeiro Bicudo (University of São Paulo, São Paulo, Brazil), Victor Sacek (University of São Paulo, São Paulo, Brazil), and Lucy Reading-Ikkanda (artist).
Syrian archaeological site
Journal Club: In Mesopotamia, early cities may have faltered before climate-driven collapse
Settlements 4,200 years ago may have suffered from overpopulation before drought and lower temperatures ultimately made them unsustainable.
Image credit: Andrea Ricci.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Birds nestling on tree branches
Parent–offspring conflict in songbird fledging
Some songbird parents might improve their own fitness by manipulating their offspring into leaving the nest early, at the cost of fledgling survival, a study finds.
Image credit: Gil Eckrich (photographer).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490