Precambrian Research 137 (2005) 167–206 Evidence for Paleoproterozoic cap carbonates in North America A. Bekker a,b,∗ , A.J. Kaufman c , J.A. Karhu d , K.A. Eriksson a a Department of Geological Sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA b Department of Earth and Planetary Sciences, Harvard University, Cambridge, MA 02138, USA c Department of Geology, University of Maryland at College Park, College Park, MD 20742, USA d Department of Geology, University of Helsinki, FIN-00014, Finland Accepted 1 March 2005 Abstract The early Paleoproterozoic Snowy Pass Supergroup of the Medicine Bow Mountains and Sierra Madre, Wyoming, USA and Huronian Supergroup, Ontario, Canada were deposited along the present-day southern flank of the Wyoming and Superior cratons. Whereas three discrete levels of glacial diamictite are developed in both successions, carbonate strata are known only directly above the middle diamictite (Vagner and Bottle Creek formations in Medicine Bow Mountains and Sierra Madre, respectively, and Espanola Formation in southern Ontario) in these thick correlative siliciclastics-dominated strata. The carbonates from each succession record negative ␦13 C values (−4.0 to −0.8‰, V-PDB) and attenuated carbon isotopic difference between organic and inorganic phases. Oxygen in carbonates is strongly depleted in 18 O suggesting exchange with hot fluids, which is consistent with pervasive recrystallization of carbonates and remobilization of elements. However, the stratigraphic coherence of carbon isotopic compositions and the general lack of correlation between ␦13 C and either ␦18 O values or trace element concentrations supports a primary origin for 13 C-depleted carbonates, which are interpreted here to reflect anomalous oceanic compositions. The intimate association of thick carbonate units containing abundant carbonate debris flows with immediately underlying glacial strata indicates that chemical precipitation resulted from a rapid flux of carbonate alkalinity onto ocean margins during post-glacial transgression. Although these early Paleoproterozoic carbonates are similar to Neoproterozoic ‘cap dolomites’ in stratigraphic position and carbon isotopic compositions, the older post-glacial accumulations begin with limestone and lack many of the sedimentary structures typical of Neoproterozoic deposits. Furthermore, it is not understood why carbonates only occur above the middle of the three glacial horizons whereas these deposits are ubiquitous above Neoproterozoic diamictites. The differences might reflect lower overall carbonate saturation in early Paleoproterozoic oceans which contrasts sharply with Archean and later Paleoproterozoic intervals and higher siliciclastic inputs in rift environments, which shut down carbonate deposition. Geological and geochemical indicators suggest a stepwise increase in atmospheric oxygen across the Paleoproterozoic glacial epoch. The tempo and mode of atmospheric oxygen rise has significant consequences for the abundance of the important greenhouse gases CH4 and CO2 and hence for oceanic acidity. If we accept that atmospheric oxidation of methane to carbon ∗ Corresponding author at: Geophysical Laboratory, Carnegie Institution of Washington, 5251 Broad Branch Rd., N.W., Washington, DC 20015, USA. Fax: +1 202 478 8901. E-mail address: [email protected] (A. Bekker). 0301-9268/$ – see front matter © 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.precamres.2005.03.009 168 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 dioxide resulted in each of the three discrete glaciations, it implies that atmospheric CH4 remained high throughout the interval and that pulsed oxidation events, plausibly linked to higher primary productivity and lower hydrothermal activity, led to surface refrigeration. If correct, the unique presence of cap carbonate above the middle Paleoproterozoic diamictite may reflect an appropriate balance of CO2 and CH4 sufficient to provide enough alkalinity to seawater through silicate weathering, but not so high that carbonate preservation would be inhibited by enhanced acidity. © 2005 Elsevier B.V. All rights reserved. Keywords: Chemostratigraphy; Paleoproterozoic; Cap carbonates; Huronian Supergroup; Snowy Pass Supergroup; Rise of atmospheric oxygen 1. Introduction At the beginning and end of the Proterozoic Eon (2500–543 Ma), the Earth was covered on several separate occasions by large continental ice sheets that may have extended to sea-level at low latitudes (e.g. Evans et al., 1997; Williams and Schmidt, 1997; Schmidt and Williams, 1999; Sohl et al., 1999; Kempf et al., 2000). A comparative analysis of the two glaciogenic intervals reveals striking similarities, in particular both are coincident with supercontinent rifting (Young, 1988), and differences, most notably the unusual presence of carbonates and iron formations in the younger, but their general absence in the older. Neoproterozoic glacial diamictites deposited during the rifting and breakup of the Neoproterozoic supercontinent known as Rodinia are commonly enriched in iron oxides with local accumulations of thick iron formations in basinal facies, and are commonly overlain by enigmatic “cap” carbonates (e.g. Kaufman et al., 1991; Kennedy, 1996; Hoffman and Schrag, 2002; Young, 2002). Geochronologic and chemostratigraphic constraints argue for at least three discrete ice ages during the Neoproterozoic (Kaufman et al., 1997; Bowring et al., 2003) which has been considered as “snowball Earth” events (Kirschvink, 1992; Hoffman et al., 1998). Initiation of these glacial episodes has been attributed to various processes, including: (1) CO2 sequestration due to high rates of primary productivity and subsequent organic carbon burial (Kaufman et al., 1997), (2) intense weathering of the high-standing supercontinent during rifting and breakup (Hoffman and Schrag, 2000; Goddéris et al., 2003; Donnadieu et al., 2004), (3) oxidation of a pre-glacial atmosphere dominated by methane (Schrag et al., 2002; Halverson et al., 2002), or (4) changes in ocean circulation associated with the low latitudinal position of Rodinia (Smith and Pickering, 2003). Early Paleoproterozoic ice ages (≤2.45 to >2.22 Ga) are believed to be associated with assembly and rifting of a Late Archean supercontinent (Kenorland of Williams et al., 1991), which was likely positioned in low latitudes (Evans et al., 1997). Although up to three discrete episodes of glaciation are recorded, there are no true iron formations in this interval and very few carbonates. Based on paleomagnetic, chemostratigraphic, and lithostratigraphic data, at least one of these ice ages has been interpreted as an older snowball Earth event potentially related to the transition from an anoxic methane-rich atmosphere to an oxygenated atmosphere with high CO2 levels (Evans et al., 1997; Kirschvink et al., 2000; Pavlov et al., 2000; Bekker et al., 2001). Lacking thick, well-preserved carbonates, little is known of temporal variations in carbon and strontium isotope composition of seawater during the early Paleoproterozoic glacial epoch. However, chemostratigraphic studies of later Paleoproterozoic carbonates from Fennoscandia and elsewhere defined a remarkable positive carbon isotope excursion with an age between ca. 2.22 and 2.10 Ga (Karhu, 1993; Karhu and Holland, 1996; Bekker et al., 2003). The early Paleoproterozoic Huronian and Snowy Pass supergroups in southern Ontario and Wyoming preserve lithologic and geochemical (chemical index of alteration of siliciclastic units) evidence of three separate glaciations with intervening periods of warm climate (Fig. 1; Young, 1991). In both successions, carbonates occur only above the middle glacial diamictite. Veizer et al. (1992) studied post-glacial carbonates of the Espanola Formation, Huronian Supergroup, Ontario and interpreted their negative ␦13 C values as reflecting a lacustrine environment of deposition. In contrast, carbonates above the glacial diamictite of the lower Duitschland Formation, which is sandwiched between two glacial deposits of the early Paleoproterozoic Pretoria Group, South Africa have negative carbon isotope values interpreted to reflect seawater A. Bekker et al. / Precambrian Research 137 (2005) 167–206 169 Fig. 1. Correlation of the Snowy Pass Supergroup, Medicine Bow Mountains, WY, USA, the Huronian Supergroup, ON, Canada, and the Transvaal Supergroup, South Africa based on tectonomagmatic event at ca. 2.48–2.45 Ga, glacial diamictites, carbonates with negative carbon isotope values (based on this paper) and mature Al- and hematite-rich quartzites overlying the middle and upper glacial diamictites respectively, and chemostratigraphy of post-glacial carbonates (Bekker et al., 2001, 2003). References to age constraints: 1: U–Pb TIMS baddeleyite age (Andrews et al., 1986); 2: U–Pb TIMS zircon age for volcanic rocks (Krogh et al., 1984); 3: U–Pb TIMS zircon and baddeleyite ages (Premo and Van Schmus, 1989; Snyder et al., 1995; Cox et al., 2000); 4: U–Pb TIMS detrital zircon age (Premo and Van Schmus, 1989); 5: Pb–Pb whole rock isochron age (Cornell et al., 1996); 6: Re–Os isochron age of early diagenetic pyrite (Hannah et al., 2004); 7: U–Pb SHRIMP zircon age of ash beds (Pickard, 2003); 8: U–Pb SHRIMP zircon age of ash beds (Martin et al., 1998) (Bekker et al., 1996; Buick et al., 1998; Master et al., 1993; Schidlowski et al., 1975; Swart, 1999). composition (Bekker et al., 2001). The carbonates of the lower Duitschland Formation are possibly correlative with carbonates immediately above the middle glacial diamictite in the North American successions based on lithostratigraphic and chemostratigraphic data (Bekker et al., 2001, 2003), but geochronologic constraints are lacking. Notably, the upper Duitschland Formation provides the only example of a positive carbon isotope excursion in Paleoproterozoic pre-glacial carbonates. Worldwide, the end of the Paleoproterozoic glacial epoch is marked by mature, Al-rich quartzites suggesting climate change to greenhouse conditions favoring strong chemical weathering (Young, 1973; Karlstrom et al., 1983; Marmo, 1992; Ojakangas, 1997; Bekker et al., 2001), quite likely under a high CO2 atmosphere. The ca. 2.22–2.10 Ga carbon isotope excursion (Karhu and Holland, 1996) also follows the ice ages, but the temporal relationship between the beginning of the carbon isotope excursion and 170 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 the onset of extreme weathering is poorly constrained. Whereas glacial deposits are known throughout Earth’s long history, only Proterozoic glacial diamictites appear to be immediately overlain by unique and isotopically anomalous carbonates deposited during post-glacial transgression. Neoproterozoic cap carbonates have been additionally defined by their unusual sedimentary structures (Corsetti and Kaufman, 2003), which may have resulted from anomalously high carbonate alkalinity in the ocean driving rapid carbonate precipitation. With the exception of the unusual sedimentary structures, the known post-glacial Paleoproterozoic carbonates may also be accepted as “cap carbonates”, and are considered so in this paper, which describes sedimentology and chemostratigraphy of the carbonates in the early Paleoproterozoic Snowy Pass Supergroup of Wyoming and Huronian Supergroup of Ontario. Results of this study are presented to (1) evaluate the environment of deposition of the carbonate rocks; (2) provide a record of temporal carbon isotope variations during the Paleoproterozoic glacial epoch and its relationship to climatic changes, which may provide constraints on global correlations, and (3) suggest a model for deposition of the Paleoproterozoic postglacial carbonates, which are linked to changes in the redox state of the atmosphere. 2. Regional setting and stratigraphy 2.1. Huronian Supergroup The best exposed and studied early Paleoproterozoic successions in North America are the thick, predominantly siliciclastic packages preserved along the southern margin of the Superior Craton. The Huronian Supergroup outcrops along the north shore of Lake Huron, Ontario, Canada and in a series of large outliers east to Cobalt and north to Timmins (Young, 1991; Fig. 2). The lower part of the Huronian Supergroup was deposited in a rift setting possibly connected with an ocean to the east, whereas the upper part starting with the Gowganda Formation is interpreted to represent a passive margin succession dominated by siliciclastic sediments (Zolnai et al., 1984; Young et al., 2001). Local disconformities and unconformities at the base of the Serpent and Gowganda formations (Young and Church, 1966) as well as abundant clastic dikes, breccias, and syn-depositional faults contemporaneous with deposition of the Espanola Formation provide evidence for tectonic reactivation (Eisbacher, 1970; Young, 1972), likely related to extension accompanying rift-drift transition. Further evidence of extension and related episodic large-scale gravitational displacement includes ‘zig-zag’ or diffuse boundaries between units, chaotic folding and incorporation of large fragments of sediments into adjacent units (Long et al., 1999). The Huronian Supergroup is subdivided by unconformities into four groups; the upper three are climatically controlled cycles with basal glacial diamictites overlain by deltaic shale or carbonate and thick fluvial sandstones (Young et al., 2001). The only thick and extensive carbonate unit of the Huronian Supergroup, the Espanola Formation, occurs directly above the middle glacial diamictite (Bruce Formation) and is overlain by fluvial sandstones of the Serpent Formation. The basal Huronian Supergroup rests unconformably, with locally preserved paleosol, on basement rocks of the Archean Superior Province (Prasad and Roscoe, 1996). Basal sediments above the reduced paleosol contain conglomerates with detrital pyrite and uraninite, suggesting low atmospheric pO2 (Livingstone Creek and Matinenda formations; Fig. 1). Interlayered mafic volcanic rocks and intrusive contacts with the Murray and Creighton granites constrain the age of the basal Huronian to 2.49–2.45 Ga (Krogh et al., 1984, 1996; Smith and Heaman, 1999), whereas the whole Huronian Supergroup is cut by the 2217.5 ± 1.6 Ma Nipissing sills and dikes (Andrews et al., 1986). Red sandstone and shale in the upper part of the Gowganda Formation and the oxidized Ville Marie paleosol below the red beds of the Lorrain Formation provide the first evidence for an oxidizing atmosphere (Wood, 1979; Rainbird et al., 1990; Panahi et al., 2000; Young et al., 2001). Pseudomorphs after anhydrite in the Gordon Lake Formation (Chandler, 1988) are consistent with an intensification of the oxidative part of the S cycle. Al-rich quartzites of the Lorrain Formation above the glacial diamictite of the Gowganda Formation imply rigorous chemical weathering and therefore rapid climate amelioration after the glaciation (Young, 1973; Long et al., 1999). Facies change dramatically across the latitudinal Murray Fault Zone (M.T.Z.; see Fig. 2), which separates shallow water sequences with A. Bekker et al. / Precambrian Research 137 (2005) 167–206 171 Fig. 2. Location of Huronian rocks and studied sections of the Espanola Formation (modified from Fedo et al., 1997). Insets show position of the Huronian Basin in North America and along the north shore of the Lake Huron on the southern margin of the Archean Superior Province. Arrows within the Huronian outcrop represent generalized paleocurrent directions for the Huronian Supergroup, M.F.Z. is Murray Fault Zone. subgreenschist to lower greenschist facies of metamorphism (Card, 1978) and open folds to the north from deep-water, thicker sequences with higher metamorphic grade and tight upright folds to the south. Paleocurrent indicators in the lower Huronian sandstones along with facies analysis and thickness changes suggest a paleoslope towards the southeast (Fralick and Miall, 1989; Rousell and Long, 1998) locally complicated by basement highs (Robertson, 1976, 1977). 2.2. Snowy Pass Supergroup, Wyoming Paleoproterozoic sedimentary and subordinate volcanic rocks are well exposed along the southeastern margin of the Wyoming Craton in the Sierra Madre and Medicine Bow Mountains, Wyoming (see inset in Fig. 3). These sequences were deformed during the 1.78–1.74 Ga Medicine Bow Orogeny when the passive margin of the Wyoming Craton collided with island arcs that accreted to the southern margin of Laurentia (Karlstrom and Houston, 1984; Chamberlain, 1998). The early Paleoproterozoic sequences of the Medicine Bow Mountains experienced greenschist facies grade of metamorphism reaching biotite isograde and are preserved in broad open folds (Houston et al., 1981). The maximum age of the supergroup is constrained by U–Pb TIMS dating of detrital zircons from the Magnolia Formation in the Sierra Madre to <2460 Ma (Premo and Van Schmus, 1989). The Snowy Pass Supergroup either unconformably overlies or is in structural contact 172 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 with Late Archean successions. The supergroup is subdivided by rotated thrusts into the Deep Lake Group, and the lower and upper parts of the Libby Creek Group (Figs. 1 and 3; Houston et al., 1992). Both the Deep Lake Group and the lower Libby Creek Group contain glaciogenic sediments. The Deep Lake Group and the lower Libby Creek Group were deposited in rift and passive margin settings, whereas the upper Libby Creek Group is considered to reflect a transition from passive margin to foredeep (Karlstrom et al., 1983) or dissection of a mature passive margin (Bekker and Eriksson, 2003). The supergroup was intruded by mafic sills and dikes that range in composition from tholeiitic to calc-alkaline. Tholeiitic gabbros that intruded prefolded Deep Lake Group are likely related to rifting (Karlstrom et al., 1981) and were dated in the Sierra Madre and Laramie Mountains at 2.09–2.01 Ga (U–Pb TIMS zircon and baddeleyite ages, Premo and Van Schmus, 1989; Snyder et al., 1995; Cox et al., Fig. 3. (A and B) Detailed maps of two areas in the Medicine Bow Mountains (after Houston and Karlstrom, 1992) with sampled localities of Vagner carbonates shown. Smaller insets on Fig. 3A show the location of Wyoming in North America and Sierra Madre and Medicine Bow Mountains with exposed early Paleoproterozoic sedimentary successions north of the Cheyenne Belt (C.B.) in Wyoming. Larger inset on Fig. 3A (after Houston and Graff, 1995) shows location of Fig. 3A and B in the northern part of the Medicine Bow Mountains. A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Fig. 3. (A and B) (Continued ). 173 174 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 2000). Pegmatitic metagabbro intruding the Cascade Quartzite in the Sierra Madre provides a minimum age for the Deep Lake Group in the Medicine Bow Mountains of ca. 2.09 Ga (U–Pb TIMS zircon age, Premo and Van Schmus, 1989). The Snowy Pass Supergroup commences with fluvial rift deposits, including quartz-pebble conglomerates with detrital pyrite and uraninite (Magnolia Formation) and contains three glacial units sandwiched between quartzites, metapelites, and a thin carbonate unit that overlies the middle glacial horizon (Vagner Formation). A transition from intracratonic rift to open marine settings coincides with deposition of the Vagner Formation (Karlstrom et al., 1983). The youngest glacial diamictite is overlain by thick, mature, Alrich and hematite-bearing quartzites (Medicine Peak Quartzite). The Paleoproterozoic succession in the Sierra Madre is strongly deformed and dismembered by a late cataclastic event (Houston, 1993), but, based primarily on lithostratigraphic similarities, it is considered to be correlative with the Snowy Pass Supergroup in the Medicine Bow Mountains (Houston et al., 1992). Possible correlative horizons in these successions and in the Huronian Supergroup include a ca. 2.45 Ga siliciclastic rift sequence with basal conglomerates containing detrital uraninite and pyrite, three glacial diamictites, carbonates above the middle glacial diamictite, Alrich hematite-bearing mature quartzites, and carbonate units above the upper glacial diamictite with similar 13 C-enrichment (Fig. 1; Bekker et al., 2003). 3. Analytical methods The carbonate-dominated sections were measured and fresh carbonate samples lacking shear, foliation, or fracture were collected from carbonate intervals. Samples collected for whole rock analysis were cleaned, cut, and powdered; these are marked by asterisk in Appendix A. Other samples were petrographically characterized and least altered (lacking veins, discoloration, weathering rinds, and silicification) and finest-grained portions of polished thick sections were microdrilled with 1 mm diamond drills. Thin sections contain xenotopic mosaics of anhedral, coarsely crystalline to very coarsely crystalline crystals of either dolomite or calcite with pseudospar and veins and frac- tures filled with quartz, very coarsely crystalline carbonate and muscovite. Whereas micrite and spar did not survive through neomorphism and recrystallization during diagenesis and metamorphism, gross primary textures (e.g. laminations) are still present. Mn, Sr, and Fe concentrations and all isotopic measurements were completed and quantified using procedures outlined in detail by Bekker et al. (2003). Samples with carbonate content below 50% were not analyzed for major and trace elements. 4. Huronian Supergroup 4.1. Espanola Formation The post-glacial Espanola Formation outcrops across a wide region in both the western and southern parts of the Huronian Basin. However, due to poor exposure and erosional contacts with overlying units, a complete section was only measured in a drill core from Quirke Lake (proximal section; Fig. 4) and at Moose Point in the Whitefish area (distal section; Fig. 5). In addition, partial sections of the Espanola Formation were measured and sampled in proximal sections (Geneva Lake, Echo Lake, and outcrops in the Quirke Lake area; Figs. 6–8, respectively). In studied sections, the Espanola Formation either rests unconformably on Archean basement or comformably on the glacial diamictites of the Bruce Formation. When exposed, the lower contact with the Bruce diamictite is sharply overlain by an up to 3 m-thick transgressive shale interval. The thickness of the formation is highly variable over the outcrop area ranging from as little as 15 m in proximal sections to over 565 m in distal sections. The upper contact is either comformable and gradational with the Serpent Formation or disconformable with the Gowganda Formation. The Espanola Formation in the Quirke Lake area (Fig. 2; see Appendix B for detailed description of all studied sections) is informally subdivided – in ascending order – into limestone, siltstone, and dolostone members (Young, 1973a). The formation thickens to the south of the Murray Fault Zone with a notable increase in thickness of the siltstone member. In the southern area, the upper dolostone member is discontinuous; it is included in the siltstone member which is overlain by a heterolithic member (Bernstein and A. Bekker et al. / Precambrian Research 137 (2005) 167–206 175 Fig. 4. Measured section of the Espanola Formation from the Kerr-McGee Corp. drill core 150/1 collared in the Bouck Township 150. ␦13 C values of carbonate and TOC are shown. Young, 1990). The lower two members were deposited in a subtidal settings whereas the heterolithic member was tide- and storm-influenced (Bernstein and Young, 1990). Stromatolites and other evidence for biological activity are conspicuously lacking in this formation except that crustose and columnar stromatolites were found in loose blocks in the Denison Mine area, Quirke Lake (Hofmann et al., 1980) and stratiform stromatolites in the Geneva Lake area (Fig. 9A). The lower limestone member in proximal sections includes, in addition to limestone, thin graded beds of sandstone and siltstone with erosional bases and crossbedding, shales with convolute bedding, rare stratiform stromatolites, intraformational flat-pebble conglomerates (Fig. 9B) that are sometimes confined to channels, and rare meter-size slumped blocks. Synsedimentary folding and faulting, soft sediment deformation, injection and flame structures are abundant in this member and indicate episodically high sedimentation rates and slope instability. In addition, distal sections contain laminations and climbing ripples, suggesting rapid accumulation. The gradationally overlying siltstone member consists of cross-bedded and laminated argillaceous siltstone and shale with graded beds, ball-and-pillow structures (Fig. 9C) and Bouma beds, intraformational conglomerates with limestone and shale pebbles, and thin carbonate beds and lenses. Shales at the top of sandy beds in distal section have flaser and wavy bedding and both symmetrical and asymmetrical wave ripples with internal cross-laminations. Both this member and the overlying dolomite member 176 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Fig. 5. Measured section of the Espanola Formation on the Moose Point (see Appendix B for location) with ␦13 C variations shown. Note that only the basal part of the heterolithic member is shown as the rest of the member is covered by water. A. Bekker et al. / Precambrian Research 137 (2005) 167–206 177 Fig. 6. Measured sections of the undivided Espanola Formation at Localities 2 and 3 in the Geneva Lake area (see Appendix B for location) with ␦13 C variations shown. Sawtooth pattern at the base of the Espanola Formation indicates unconformity. Due to poor outcrop only partial section of the Espanola Formation was measured at Locality 2. contain soft sediment deformation and injection structures. Synsedimentary faults overlain by disrupted beds and 35 m thick intraformational breccia-conglomerate with boulders up to 1 m in size are present at Moose Point. The uppermost dolomite member in the proximal sections consists of graded beds of carbonate, shale, and intraclasts. Scour surfaces, “beach rosettes”, abundant wave ripples with internal cross-laminations (Fig. 9D), mudcracks, cross-bedding, rare fenestral textures, and peloids indicate a shallow-water environment above wave base. The heterolithic member at Moose Point consists of upward-fining cycles of sandstone, siltstone, shale, and thin impure carbonates. Small conglomerate channels, herringbone cross-stratification, wave ripples, and trough and large low-angle cross-bedding indicate a tidally influenced shallow-water environment. 4.2. Serpent and Gowganda formations Carbonates and calcareous siliciclastic beds have been previously described from units that overlie the Espanola Formation (e.g. Card et al., 1977; Long, 1976). For this study, carbonates of the Serpent Formation were sampled on highway 6 south of Loon Lake, 178 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Fig. 7. Measured section of the limestone member, Espanola Formation at Marble Point of Echo Lake (see Appendix B for location) with ␦13 C variations shown. 100–200 m south of La Chance Drive (Appendix A). Yellow carbonates form two massive and resedimented layers with soft sediment deformation structures. The lower layer is 1 m thick and contains in its upper part shale and sandstone clasts and carbonate rip-ups in sandy carbonate and sandstone matrix. Carbonate layers and lenses of the Serpent Formation were also sampled on the eastern side of the East Simpson Island in the Bay of Islands. Carbonate lenses in laminated shales of the Gowganda Formation were sampled on the portage between Ishmael and Low Lakes. On close inspection, most carbonates of both the Serpent and Gowganda formations appear to be carbonate cements in sandstone or post-sedimentary features likely related to hydrothermal activity. The only known exception is sedimentary carbonates of the Serpent Formation near Loon Lake. 4.3. Geochemical data All geochemical data for the Espanola Formation and the overlying Serpent and Gowganda formations are tabulated in Appendix A. Carbon and oxygen isotope data for Espanola Formation limestone and dolomite samples show a large range with consistently negative carbon isotope values in all sampled sections A. Bekker et al. / Precambrian Research 137 (2005) 167–206 179 Fig. 8. Measured section of the Espanola Formation in the Denison Mine area, Quirke Lake (see Appendix B for location). Note that due to poor outcrop, sample heights and thickness of the members are approximate. However, relative stratigraphic position of samples is well established. The siltstone member was not logged in this area due to limited outcrop. 180 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Fig. 9. Sedimentary structures in carbonate rocks of the Espanola Formation. A—stratiform stromatolites in the Geneva Lake area, coin is 3.5 cm in diameter; B—intraformational conglomerate, Marble Point of Echo Lake, scale is 10 cm in length; C—ball-and-pillow structure in the siltstone member of the Espanola Formation, Moose Point area; D—ripple marks with unidirectional cross-laminations and bundled upbuilding indicating wave activity, coin is 2 cm in diameter. A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Fig. 10. Sedimentary structures in the Vagner Formation. A—angular lonestones in coarse-grained subarkosic matrix, pen is 14 cm long; B—finely laminated carbonates at the base of the section at the BE-C-95-8 locality, hammer is 70 cm long; C—flat-pebble conglomerate at the top of the section at the BE-C-95-8 locality, pencil for scale. 181 182 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 varying from −6.7 to −0.8‰ V-PDB. Oxygen isotope values have a much larger range (∼10‰) than those of carbon isotopes and are significantly depleted in 18 O; all samples have ␦18 O values of less than −10‰. Manganese, iron, and Sr contents in these samples also range widely. Total organic carbon (TOC) contents are low in most samples from the Espanola Formation, ranging from 0.03 to 1.23 mg C/g sample for the whole dataset with most samples having TOC content below 0.2 mg C/g sample. Carbon isotope values of organic matter are highly variable. Measured differences between organic and carbonate carbon yield δ values of less than 26.6‰; the average for the best-preserved drill core samples is around 20‰. Carbonates of the Serpent Formation in the Loon Lake area have highly variable ␦13 C values ranging from −6.3 to +0.9‰ whereas strongly negative ␦18 O values show much smaller range (−21.0 to −19.8‰; Appendix A). These samples have also high Mn (4170–4406 ppm) and Sr (254–463 ppm) contents and high Mn/Sr ratios (9–17.3). 5. Snowy Pass Supergroup, Wyoming 5.1. Vagner Formation The Vagner Formation overlies the Cascade Quartzite and is up to 800 m thick (Karlstrom et al., 1983; Fig. 1). A three-fold subdivision is recognized, including (1) a basal diamictite up to 300 m thick; (2) a middle marble unit with variable thickness from 5 to 60 m; and (3) an upper quartz-rich phyllite and fine-grained quartzite unit (Karlstrom, 1977; Karlstrom and Houston, 1979). An unconformity at the base of the Vagner Formation is inferred from the variable thickness of the underlying Cascade Quartzite, the presence of channels along the contact, and pebbles of Cascade Quartzite in the basal Vagner Formation (Karlstrom, 1977; Bekker, 1998). The diamictite includes angular to subangular clasts in a greenish to reddish coarse-grained subarkosic matrix (Fig. 10A). The basal diamictite has been interpreted as a glaciomarine deposit, based on the presence of dropstones in laminated phyllite, faint stratification in some conglomerate layers, and major element analyses of finegrained matrix from the diamictite (Sylvester, 1973; Houston et al., 1981). The contact with the overlying carbonate unit was not observed in the field but is considered to be conformable (Karlstrom et al., 1983). The marble unit consists of 1–5 cm-thick rhythmically interlayered massive, fine-grained, grey–blue to brown carbonate and brown coarse-grained siliciclastics-rich carbonate. Contacts between these carbonates are transitional. These carbonates are recrystallized into a coarse to very coarse calcite mosaic. Convolute bedding, thin laminations, normal grading, small-scale hummocky cross-stratification, lensoidal bedding, and flat-pebble conglomerate are restricted to the upper part of this unit. The marble unit locally includes basal rhythmites of sandstone–siltstone– carbonate (e.g. BE-C-95-24 locality; Fig. 3A). Phyllite with lonestones of granitic composition ranging from 1 to 5 cm in diameter is interlayered with carbonate in some localities. The only relatively complete section of carbonates (BE-C-95-8 locality; Fig. 3B) is in close proximity to a mafic sill. This section has a shallowing-upward trend; laminated carbonates (Fig. 10B) with stratiform stromatolites and small oölites occur at the base and grade upsection into massive and siliciclastics-rich carbonates with flat-pebble conglomerate (Fig. 10C) at the top. Quartzites in the overlying phyllite-quartzite unit contain climbing ripples, plane and wavy bedding, and rare cross-beds and ripples. Carbonate clasts in this unit are rarely present; but neither warping nor piercing of layers beneath clasts was observed. Due to poor exposure, only partial sections of the marble unit were measured and sampled (Fig. 11). Additional samples were collected across the whole outcrop area of the marble unit (Fig. 3A and B) from siliciclastics-poor carbonates to test variability of stable isotope compositions. Samples from the marble unit of the Vagner Formation consist of low-Mg calcites (Bekker, 1998; Appendix A) and have variable contents of Mn (295–1461 ppm) and Sr (21–217 ppm), and a wide range of Mn/Sr ratios (1.8–16). Oxygen isotope values are strongly depleted in 18 O ranging from −20.9 to −14.5‰. Carbon isotopes are uniformly depleted in 13 C with ␦13 C values ranging from −3.8 to −1.5‰. Analyzed samples have very low TOC content (0.05–0.3 mg C/g sample) which is highly 13 Cenriched (−19.6 to −11.0‰) resulting in small measured difference in ␦13 C values between carbonate and A. Bekker et al. / Precambrian Research 137 (2005) 167–206 183 Fig. 11. Two closely spaced partial sections (BW-series samples in Appendix B) of the marble unit, Vagner Formation sampled near the small lake to the north of Dipper Lake (41◦ 22 23 N, 106◦ 20 39 W; see Fig. 3A for location) with ␦13 C variations shown. organic phases. Two samples with relatively high Sr contents and low Mn/Sr ratios were analyzed for Sr isotopes; the lowest 87 Sr/86 Sr value of ∼0.7081 is significantly lower than those measured from the correlative Espanola Formation (0.71128, Veizer et al., 1992). 5.2. Bottle Creek Formation, Snowy Pass Group, Sierra Madre The Bottle Creek Formation (sensu Houston et al., 1992; previously named as Vagner Formation in Graff, 1978) is considered correlative with the Vagner to Heart formations of the Medicine Bow Mountains based on lithostratigraphic similarity (Fig. 1). It includes diamictite that is overlain successively by marble, quartzite, a second diamictite, and quartzite (Houston et al., 1992). Carbonate is confined to the western part of the Sierra Madre (see Graff, 1978). The lower diamictite and greenish, fine-grained, finely-laminated marble, which is up to 30 m thick and interlayered with sandstone and siltstone, are considered as glaciomarine offshore facies equivalents of the Vagner Formation in the Medicine Bow Mountains (Houston et al., 1992). Samples were collected about 1.6 and 1.9 km to the east along the Big Sandstone Creek from the intersection with the Rawlings Road (see Appendix A). 184 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Like their Vagner Formation counterparts, oxygen isotope values of marble samples from the Bottle Creek Formation are also highly 18 O depleted, ranging from −16.1 to −12.5‰ V-PDB. Most carbon isotope values vary from −3.4 to −2.2‰ V-PDB. 6. Discussion 6.1. Diagenesis and metamorphism It is clear from the highly recrystallized textures of the Espanola, Vagner, and Bottle Creek formation carbonates as well as the significantly depleted 18 O values and elevated 87 Sr/86 Sr compositions (see Appendix A; Veizer et al., 1992) that these lonesome early Paleoproterozoic carbonates have been altered by secondary processes, thereby hampering interpretation of carbon isotope values. The low ␦18 O values, in particular, suggest equilibration with hot metamorphic or diagenetic fluids (Veizer, 1983), which must have also affected trace element concentrations. Similarly, the low abundance of TOC and the decreased fractionation between inorganic and organic carbon isotope compositions is likely related to metamorphic alteration of carbon, which can greatly alter 13 C abundances in TOC (Hayes et al., 1983). If correct, some degree of 13 C depletion in carbonate may reflect either the addition of newly-formed phases subsequent to organic matter degradation or metamorphic devolatization. However, only metamorphism of siliceous carbonates can result in 13 C depletion of carbonate through Rayleigh distillation associated with decarbonation reactions (Valley, 1986; Baumgartner and Valley, 2001) and the formation of certain metamorphic minerals (e.g. wollastonite). The lack of these minerals in most of our samples, which are predominantly pure carbonate, argues against a metamorphic explanation for the low ␦13 C values of the Espanola, Vagner, and Bottle Creek carbonates. Furthermore, mass balance considerations suggest that “organic” additions to these organic-lean carbonates should not affect the ␦13 C values by more than a per mil (cf. Pelechaty et al., 1996) and argue for a primary origin of these negative ␦13 C values. Lastly, the general lack of correlation between 13 C abundances in these carbonates and those of 18 O or elemental compositions (Figs. 12 and 13; Appendix A) support this notion. Carbon isotope values of the least altered samples of the Espanola carbonates range from −4.0 to −0.8‰, but stratigraphic trends are lacking except in the Quirke Lake drill core. In this section we note a poorly defined trend of 13 C-enrichment in carbonate and 13 C-depletion in organic matter upsection. Carbon isotope values of Vagner carbonates also do not show a stratigraphic trend (Fig. 11). In combination with regional coverage provided by the rest of samples and excluding the most altered samples, carbon isotope values of the Vagner carbonates vary from −2.6 to −1.5‰. 6.2. Paleoproterozoic cap carbonates Carbon isotope data for Espanola carbonates from both distal and proximal parts of the Huronian Basin and for the whole outcrop area of Vagner and Bottle Creek carbonates suggest that dissolved seawater bicarbonate was moderately depleted in 13 C after the second early Paleoproterozoic glacial event. A marine rather than lacustrine (cf. Veizer et al., 1992) origin is supported by: (1) the transgressive nature of the deposits; (2) the long range lithostratigraphic correlation of these unique North American carbonate accumulations; (3) the presence of herringbone cross-bedding recording tidal influence (Bernstein and Young, 1990); and (4) the size of the basins as well as the inferred tectonic setting for both units being a failed rift arm of a triple junction extending to the ocean (Young, 1983; Karlstrom et al., 1983). Insofar as all Late Archean and early Paleoproterozoic carbonates record ␦13 C values near 0‰ or higher (e.g. Bekker et al., 2003a), we relate this Paleoproterozoic negative carbon isotope anomaly to environmental changes in response to global glaciation. Support for this view comes from our proposed correlation of the North American post-glacial carbonates with carbonates from the lower Duitschland Formation, South Africa (Bekker et al., 2001). The lower Duitschland carbonate was deposited above glacial diamictite on an open-marine margin during oceanic transgression and records very similar degrees of 13 Cdepletion in carbonates with ␦13 C values ranging from −3.7 to 0.1‰, as well as 13 C-enrichment in organic matter (Bekker et al., 2001). Other Paleoproterozoic successions worldwide apparently do not contain the record of the second glacial event and overlying cap carbonate. A. Bekker et al. / Precambrian Research 137 (2005) 167–206 185 Fig. 12. ␦13 C vs. ␦18 O values of carbonates of the Espanola Formation analyzed in this study; data are grouped by locality. Analyses for outcrop and drill core samples from the Quirke Lake area are shown separately. Note that all datasets show no correlation between these two parameters. Fig. 13. ␦13 C vs. ␦18 O scatter diagram for carbonates of the Vagner Formation. Note lack of correlation between these parameters. 186 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Our interpretation of Paleoproterozoic carbonates in North America and South Africa as post-glacial “cap carbonates” relies, in part, on comparison with known Neoproterozoic occurrences. For example, both Neoproterozoic and Paleoproterozoic glaciogenic successions were deposited in similar tectonic settings. Both glaciogenic intervals coincide with supercontinent rifting (Kenorland and Rodinia for the Paleoproterozoic and Neoproterozoic eras, respectively) resulting in thick packages of siliciclastic sediments including glacial diamictites (Young, 1988; Hoffman, 1991), and variable thicknesses of carbonate deposited during post-glacial transgression. The post-glacial carbonates of both intervals are typically subtidal deposits, recognized by the abundance of fine laminations, hummocky cross-bedding, flat-pebble conglomerates, and graded event beds; many of them represent carbonate formed on shelves that was reworked and deposited in thin graded units further offshore (Kennedy, 1996). The cap carbonates of the Neoproterozoic and Paleoproterozoic are also similar in their depletion of 13 C in carbonates and 13 C-enrichment in co-existing organic matter (e.g. Kaufman and Knoll, 1995; Kaufman et al., 1997; Bekker et al., 2001). In contrast, carbon isotope values of the Vagner, Bottle Creek, and Espanola carbonates do not show convincing stratigraphic or facies-related trends in ␦13 C values, whereas heavily-sampled thick carbonatedominated Neoproterozoic successions record rapid changes in ␦13 C values with time (Kaufman et al., 1991, 1997; Kennedy, 1996; Hoffman et al., 1998; Halverson et al., 2002; Xiao et al., 2004); furthermore, on Neoproterozoic carbonate platforms ␦13 C compositions are more negative on the slope and increase across the shelf (Kennedy, 1996; James et al., 2001). Locally, where the uppermost strata of Neoproterozoic carbonate platform were not eroded below the diamictites, carbon isotope values show a trend upsection from highly positive to highly negative values (Kaufman and Knoll, 1995; Kaufman et al., 1997; Halverson et al., 2002; McKirdy et al., 2001). These isotopic trends are less obvious in Neoproterozoic siliciclastics-dominated sequences where the isotopic records are truncated. Most (but not all) of the Neoproterozoic cap carbonates are notably dolomitic at their base and are overlain by limestones that formed during maximum transgression, or by siliciclastic rocks, depending on sediment supply and other factors; the dolomite-limestone boundary is marked by the presence of a thin barite interval in Mackenzie Mountains, northwest Canada (Hoffman and Schrag, 2002). Above this limestone interval, Neoproterozoic cap carbonates become increasingly dolomitic as accommodation space filled and the basins shallowed. On the other hand, the Paleoproterozoic cap carbonates of North America and South Africa begin with deep-water limestone above thin transgressive deep-water shale. As relative sea-level fell, the limestone accumulations were overlain by shallow-water dolomite or siliciclastic sediments. Furthermore, the Paleoproterozoic post-glacial carbonates lack some of the notable sedimentary structures recognized in many Neoproterozoic examples, including: (1) giant wave ripples (Allen and Hoffman, 2005), (2) pseudomorphs after aragonite fans, (3) finelylaminated inversely-graded micro- and macropeletal dolomites with roll-up structures, (4) sheet-crack cements, and (5) tube-like stromatolites (e.g. Hoffman and Schrag, 2002; Corsetti and Kaufman, 2003). Some of these structures were related to precipitation from seawater oversaturated with carbonate alkalinity and Mg (James et al., 2001; Corsetti and Kaufman, 2003) and seismicity during isostatic rebound following the glaciation (Nogueira et al., 2003). Whereas some of these structures may have been present in the Paleoproterozoic cap carbonates and are now obscured by recrystallization, others were likely not developed at all. Some of the sedimentary structures common to Neoproterozoic cap carbonates but lacking in Paleoproterozoic cap carbonates are likely related to high rates of inorganic carbonate precipitation in shallow-marine conditions from the post-glacial ocean supersaturated with carbonate alkalinity and Mg, but with low oxygen and sulfate contents. Giant wave ripples might reflect the higher storm and cyclone intensity in response to higher surface temperature gradients in the aftermath of the snowball Earth (Emanuel, 1987; Hoffman and Schrag, 2002). Extreme alkalinity in the Neoproterozoic postglacial ocean has been related to ocean overturn at the end of glaciation (e.g. Kaufman et al., 1991; Grotzinger and Knoll, 1995; Hoffman et al., 1998) and high rates of carbonate weathering on continents (Hoffman and Schrag, 2002). Notably, some of the unusual sedimentary structures in the Neoproterozoic cap carbonates (i.e. pseudomorphs after aragonite fans, cement-filled A. Bekker et al. / Precambrian Research 137 (2005) 167–206 sheet cracks, marine cement laminae, and roll-up structures) are also common on Late Archean carbonate platforms and have been related to high carbonate alkalinity in the Archean ocean (Simonson et al., 1993; Sumner and Grotzinger, 2000). Their apparent absence in the Paleoproterozoic cap carbonates, suggesting lower overall oceanic carbonate saturation in the glaciogenic interval, might reflect that either atmospheric pCO2 levels were low, thereby affecting continental weathering, or high, causing oceanic pH to decline and carbonate to dissolve. Studies of the chemical index of alteration of mature siliciclastic sediments between diamictites in the Huronian and Snowy Pass supergroups (Sylvester, 1973; Houston et al., 1981; Nesbitt and Young, 1982; Bekker, 1998) suggest extreme weathering of continental protolith. This process ultimately produces carbonate alkalinity in the oceans through reaction with atmospheric carbon dioxide. Due to the general absence of carbonates in underlying Paleoproterozoic successions, continent derived alkalinity could have only been provided by silicate weathering, which apparently led to generally undersaturated conditions in the Paleoproterozoic post-glacial oceans and may explain the absence of the unusual sedimentary structures in the Espanola and equivalent cap carbonates. On the other hand, if atmospheric CO2 was very high, silicate weathering might have delivered ample alkalinity to seawater but the induced drop in oceanic pH and the rise in the carbonate compensation depth would delay carbonate accumulation. 6.3. A depositional model for the Paleoproterozoic cap carbonates of North America In our depositional model, post-glacial flooding of alkalinity-laden seawater temporally shut down siliciclastic influx leading to accumulation of carbonate on a relatively flat siliciclastic shelf. Based on the general lack of traction-induced sedimentary structures except near their tops, Paleoproterozoic cap carbonates in North America were deposited below wave base as reworked carbonates originally inorganically precipitated in shallower waters; this interpretation is supported by the scarcity of stromatolites and oncolites in both distal and proximal settings. The abundance of carbonate breccia in both distal and proximal sections of the Espanola Formation suggests that the carbonate 187 platform rapidly developed an oversteepened rim, thus shedding carbonate debris into deeper environments due to high rates of carbonate accumulation. The buildup of the platform caused a relative sea-level fall allowing for siliciclastic deposition to resume in marginal marine and, later, in deltaic and fluvial settings (Karlstrom et al., 1983). At this point, some margins returned to ramp-like morphology, which favored synsedimentary to early diagenetic dolomitization of the overlying delta plain deposits (dolomite member of the Espanola Formation), which were overlain by prograding fluvial sandstones of the Serpent Formation, the Huronian Supergroup (Long, 1976). 6.4. Coupled early Paleoproterozoic change of the atmospheric redox state and glaciations Recent geochemical models of the early Paleoproterozoic glaciation have focused on the oxidative conversion of atmospheric methane, a remarkably strong greenhouse gas, to carbon dioxide, which is more transparent to infrared radiation and hence less efficient as a greenhouse gas, during the rise of atmospheric oxygen (Pavlov et al., 2000). This atmospheric transition may have resulted in surface refrigeration at a time of lower solar luminosity (∼87% of modern), but stratigraphic and sedimentologic evidence argues for three discrete glacial events separated by long time intervals with greenhouse conditions. A fully integrated explanation for these ice ages then must account for three glacial events in the early Paleoproterozoic, as well as the unusual appearance of “cap carbonates” above only the middle diamictite. We hypothesize that the glacial events were driven by a misbalance between oxygen production (through photosynthesis) and consumption (through oxidation of reduced compounds in seawater), and its effects on atmospheric CH4 and CO2 concentrations; with their changes directly affecting surface temperatures. Since methane likely acted as the dominant greenhouse gas before the rise of atmospheric oxygen, the glacial epoch might represent climatic adjustments during the transition to an atmosphere with CO2 acting as a major greenhouse gas. Empirical constraints on the atmospheric CO2 and CH4 contents across the early Paleoproterozoic glacial epoch, however, are virtually nonexistent. The single geochemical estimate for 188 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 atmospheric CO2 content around this time is based on the absence of siderite in paleosols older than 2.2 Ga, which suggests an upper limit of 10−1.4 atm (∼100 PAL; Rye et al., 1995). Using the same constraint in a one-dimensional, radiative-convective climate model, Pavlov et al. (2000) suggested that Archean atmosphere should have contained up to 1000 ppmv methane to overcome significantly lower solar luminosity at that time. There are no constraints on atmospheric CO2 content in the Paleoproterozoic after the rise of atmospheric oxygen but since the methane contents should have decreased to less than 100 ppmv in an oxidizing atmosphere (Pavlov et al., 2003), atmospheric pCO2 must have increased dramatically to keep the Earth from freezing over again. Accepting that oxidation of atmospheric methane to carbon dioxide resulted in each of the three discrete glaciations requires that atmospheric pCH4 remained high throughout the interval. Constraints on the atmospheric oxygen level across the Paleoproterozoic glacial epoch may ultimately provide a test for the Pavlov et al. (2000) model. The correlation between Paleoproterozoic successions of North America and South Africa (Fig. 1) is based on their stratigraphic and temporal relationship with respect to the 2.48–2.45 Ga tectonomagmatic event, three glacial events, climatic amelioration after the last glacial event, and >2.2–2.1 Ga carbon isotope excursion. If this correlation is accepted then we can combine constraints for atmospheric oxygen level from these successions (Fig. 14). Geological and geochemical indicators of atmospheric redox state indicate that the atmosphere remained anoxic across the 2.48–2.45 Ga tectonomagmatic event (Fig. 14). These include (1) reduced paleosols at the base of the Huronian succession; (2) well-documented detrital grains of pyrite and uraninite in the Matinenda and Magnolia formations of the Huronian and Snowy Pass supergroups, respectively; (3) diagenetic sulfides with small range of ␦34 S values and non-mass dependent fractionation of sulfur isotopes, and (4) the absence of a Ce anomaly in iron formations of the Ghaap Group, South Africa. Similarly, units overlying the oldest glacial diamictites in both North American successions contain radioactive and pyritiferous conglomerates and sandstones (Roscoe, 1969; Meyn, 1970; Karlstrom and Houston, 1979a; Meyn and Matthews, 1980, 1991; Mossman and Harron, 1983; Long, 1987) suggesting that the atmospheric oxygen content remained low between the first and second glacial events. It was recently confirmed based on three S isotope study of pyrite grains from the pyritiferous conglomerates and sandstones of the Missassagi Formation, Huronian Supergroup that at least some of these grains are detrital in origin (Bekker et al., 2005). This observation constrains the maximum level of atmospheric oxygen after the first glacial event, whereas sulfur isotope data for sulfides extracted from shales and siltstones that sandwich the oldest Huronian glacial diamictite reveal a small range of 33 S values from −0.06 to +0.27‰ V-CDT (Farquhar et al., 2000; Wing et al., 2002, 2004), suggesting that immediately before the oldest Paleoproterozoic glaciation the atmosphere became more oxidizing than in the latest Neoarchean (Farquhar et al., 2000; Ono et al., 2003). There are several lines of evidence suggesting significant amounts of oxygen in the atmosphere between the second and third Paleoproterozoic glacial events including hematitic ironstone, highly negative ␦34 S values, and the lack of non-mass dependent fractionation in S isotopes in South Africa (Bekker et al., 2004a, 2004b). Puzzlingly, the Noomut Formation, below the glaciogenic Padlei Formation of the Hurwitz Group, Canada that is considered correlative with the upper glacial diamictite of the Huronian Supergroup contains variegated beds and detrital pyrite (Aspler and Chiarenzelli, 1997). Carbonates of the upper Duitschland Formation, South Africa that were also deposited between the second and third glacial events have highly positive ␦13 C compositions, suggesting an excess in organic carbon burial and a release of photosyntheticallyproduced oxygen to the atmosphere (Bekker et al., 2001). Units overlying the third and youngest Paleoproterozoic glacial diamictite contain multiple evidence for the oxygenated atmosphere (Fig. 14). Put together, these observations indicate that pCO2 rose dramatically during the early Paleoproterozoic glacial epoch, but likely fluctuated across the three glaciations. We speculate here for a stepwise increase in atmospheric oxygen across these climatic events that (1) commenced the glaciations through oxidation of atmospheric methane, and (2) explains why cap carbonate is present only above the middle of three glacial diamictites. The lack of carbonate platforms after the ca. 2.5 Ga and before the beginning of the A. Bekker et al. / Precambrian Research 137 (2005) 167–206 189 Fig. 14. Indicators of atmospheric and ocean redox state (green—reduced; red—oxidized) arranged with respect to 2.48–2.45 Ga tectonothermal events, three glaciogenic events (triangles, question marks indicate uncertainty on their age), and carbon isotope variations. Carbon isotope composition of the early Paleoproterozoic carbonates (grey curve) is from Karhu and Holland (1996) with additional data (dashed bold line) from Bekker et al. (2001) and this paper. Indicators of atmospheric and ocean redox state are organized in somewhat subjective way from the most sensitive to the oxygen level at the bottom to the least sensitive at the top. The Paleoproterozoic glacial epoch is bracketed between 2.42 and 2.22 Ga; 2.32 Ga Re–Os isochrone age of pyrites from black shales deposited before the last glacial event (Hannah et al., 2004) provide another age constraint for these glacial events. References: detrital pyrite, siderite, and uraninite (Roscoe, 1969, 1981, 1996; Roscoe and Minter, 1993; Meyn, 1970; Meyn and Matthews, 1980, 1991; Karlstrom and Houston, 1979a; Mossman and Harron, 1983; Long, 1987; Coetzee, 2001; Bekker et al., 2005); red/variegated beds, copper stratiform deposits (Wood, 1979; Rainbird and Donaldson, 1988; Chandler, 1989; Houston et al., 1992; Aspler and Chiarenzelli, 1997; Dorland, 1999; Beukes et al., 2002); ␦34 S record (Cameron, 1982, 1983; Hattori et al., 1983); 33 S record (Farquhar et al., 2000; Bekker et al., 2004a, 2004b; Wing et al., 2002, 2004); paleosols (Rainbird et al., 1990; Marmo, 1992; Rye and Holland, 1998; Panahi et al., 2000; Beukes et al., 2002a; Yang and Holland, 2003); sulfate evaporites (Chandler, 1988); Mn deposits (Beukes, 1983; Tsikos and Moore, 1998); Ce anomaly (Alibert and McCulloch, 1993; Bau et al., 1998; Murakami et al., 2001; Maynard, 2004). (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of the article) Paleoproterozoic glacial epoch might be related to the supercontinent assembly. In this case, enhanced chemical weathering would have progressively drawn down atmospheric CO2 (cf. Nance et al., 1988) eventually attenuating alkalinity fluxes, whereas terrestrial silici- clastic fluxes would increase precluding carbonate accumulation. The resulting methane-rich and extremely CO2 -poor atmosphere would have been highly unstable as photosynthetically-produced oxygen started to build up in the atmosphere. It is possible that decrease 190 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 in hydrothermal and terrestrial flux of reduced compounds during the supercontinent assembly allowed oxidation of the surface environments. As the cumulative greenhouse power of the atmospheric methane and carbon dioxide fell below threshold limits, global temperatures dropped and the first glaciation ensued. No carbonate was deposited after the first Paleoproterozoic glaciation because of persistent carbonate undersaturation. Assuming that volcanic CO2 and biogenicallyderived methane continued to be released to the atmosphere during and after the first glaciation, surface temperatures would have warmed and enhanced silicate weathering and continental fluxes of alkalinity. A following increase in the photosynthetically-produced oxygen flux resulting in drawdown of atmospheric methane and decrease in the atmospheric CO2 levels due to chemical weathering led to the second glaciation. In this case, carbonate alkalinity fluxes were sufficient to result in the formation of carbonate in the aftermath of glaciation. A further reduction of the atmospheric methane conceivably associated with the positive carbon isotope excursion in seawater composition recorded in the upper Duitschland Formation carbonates (Bekker et al., 2001) required even more elevated atmospheric CO2 levels to overcome low albedo effect and to warm surface temperatures in the aftermath of the third glaciation. Carbonate deposition was apparently delayed as evidenced by the absence or scarcity of carbonates associated with deeply weathered sediments above the ultimate Paleoproterozoic tillite until carbonate alkalinity levels in the ocean were sufficient to counteract dissolution induced by extremely high atmospheric CO2 levels. We therefore infer that the atmospheric CO2 levels in the aftermath of the second glaciation reached a level sufficient to provide enough carbonate alkalinity to seawater through chemical weathering but not so high that carbonate preservation would be inhibited by higher overall acidity. Carbonate again becomes abundant in the marine sedimentary record during the 2.22–2.10 Ga positive carbon isotope excursion. We speculate that warming of the oceans and enhanced biological uptake at that time (Karhu and Holland, 1996) drew down oceanic CO2 so that alkalinity levels in the ocean were sufficient to counteract dissolution. The concurrent rise of pO2 associated with the ca. 2.22–2.10 Ga carbon isotope excursion likely swept most of the remaining methane from the atmosphere, leaving CO2 as the dominant greenhouse gas thereafter (cf. Kaufman and Xiao, 2003). We conclude that although Neoproterozoic and Paleoproterozoic glacial epochs have many similarities, the early Paleoproterozoic glacial epoch has certain unique characteristics that likely reflect the transition from an anoxic methane-rich atmosphere to the oxygenated atmosphere with carbon dioxide being a dominant greenhouse gas. This speculative model of stepwise or pulsed atmospheric oxidation provides a plausible explanation for multiple early Paleoproterozoic glaciations, the general absence of carbonate throughout the glaciogenic interval, and the singular presence of a cap carbonate in the middle of the glacial epoch. Acknowledgements This project started when A.B. was a MS student at the University of Minnesota, Duluth under the direction of R.W. Ojakangas. D. Long, G. Bennett, G.M. Young, and P.J. Graff helped in field and with logistics. D. Long pointed to outcrops in the Geneva Lake area with stratiform stromatolites. K. Chamberlain provided helpful discussions on the tectonic setting and geochronology of the Snowy Pass Supergroup. D. Rimstidt provided access to and helped with AAS work. Support came from Department of Geology, University of Minnesota, Duluth; GSA and southeastern section of GSA; Colorado Scientific Society; and the graduate school and Department of Geological Sciences, VPI and SU. Funding for participation by AJK and isotopic analyses at the University of Maryland was provided by NSF grants EAR-98-17348 and EAR-0126378 and NASA Exobiology grant NAG 512337. JAK thanks the Geological Survey of Finland for permission to publish data produced in their laboratories and A. Henttinen for assistance. AB gratefully acknowledges financial support from PRF/ACS 270953 (37194-AC2) grant to H. D. Holland, Harvard University. Constructive reviews by R. H. Rainbird and an anonymous reviewer are gratefully acknowledged. AB dedicates this paper to fond memories of Olga M. Haring. A. Bekker et al. / Precambrian Research 137 (2005) 167–206 191 Appendix A. Supplementary materials Carbon, oxygen, and strontium isotope values and trace and major element contents of studied carbonates Sample number Height above base (m) Espanola Formation Geneva Lake Locality 1 GE-4 4.2 GE-5 8.7 Mineral/ rock ␦13 C, ‰ V-PDB ␦18 O, ‰ V-PDB Calc. si-ne Calc. si-ne −5.8 −18.6 Calcite Calcite Calcite Cal. sa-ne Calcite Calcite Calcite Calcite −1.3 −0.8 −1.5 −1.3 −2.1 −1.0 −0.9 −14.8 −14.6 −14.8 −12.8 Locality 2 GE-8 GE-9 GE-11 GE-12 GE-13 GE-14 GE-16 GE-17 0.0 1.6 3.0 3.3 5.3 7.6 9.2 10.2 Locality 3 GE-22 GE-23 GE-24 GE-25 2.7 3.2 5.4 5.7 Calcite Calcite Calcite Calcite −1.0 −1.2 −1.4 −2.3 −16.6 −11.7 −14.1 −11.7 Echo Lake EL-5 EL-6 EL-7 EL-8 EL-9 EL-10 EL-11 EL-12 EL-13 EL-14 EL-15 EL-16 EL-17 EL-18 EL-19 EL-20 EL-21 EL-22 EL-23 EL-24 EL-25 EL-26 El-27 EL-28 EL-29 EL-30 EL-31 EL-32 EL-33 EL-34 EL-34A 1.0 2.2 2.7 3.7 5.2 6.7 8.2 9.7 11.0 12.0 13.0 14.5 16.0 17.2 18.0 19.5 21.0 22.5 24.0 25.5 27.0 28.5 30.0 31.0 32.0 33.5 35.6 37.7 39.8 40.0 40.0 Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite −2.3 −2.0 −2.1 −2.2 −2.6 −1.9 −2.0 −2.1 −1.7 −1.6 −2.1 −2.2 −1.9 −2.1 −2.0 −1.9 −2.1 −1.8 −1.8 −2.0 −2.1 −2.2 −2.3 −2.2 −2.0 −1.9 −2.0 −1.2 −2.1 −2.0 −1.1 −14.4 −15.2 −15.3 −15.7 −16.5 −15.4 −16.0 −15.0 −15.7 −16.4 −13.9 −16.8 −16.1 −15.9 −14.6 −13.1 −14.6 −14.1 −14.2 −14.6 −14.0 −11.8 −14.8 −12.7 −12.2 −15.1 −15.7 −14.5 −14.5 −14.8 −14.1 −11.1 −12.9 −12.7 TOC, mg C/g sample ␦13 CTOC , ‰ V-PDB 0.08 −11.7 δ Mn (ppm) Sr (ppm) 2053 tr. 1301 1378 1987 50 tr. 4 527.2 19.3 363 226 182 437 240 118 113 20 1.5 1.9 1.6 22.2 10396 9963 6047 8092 tr. 40 9 tr. 248.8 687.4 261 166 188 165 154 223 281 173 337 99 262 248 361 277 130 174 251 65 2.6 0.6 0.8 0.5 0.6 1.7 1.6 0.7 5.1 294 221 205 235 221 332 212 171 163 176 83 319 234 93 458 153 168 141 172 208 151 143 156 224 329 154 98 165 198 191 185 114 68 168 57 46 209 280 316 327 200 372 2.1 1.4 0.9 0.7 1.4 3.4 1.3 0.9 0.9 1.0 0.7 4.7 1.4 1.6 10.0 0.7 0.6 0.4 0.5 1.0 0.4 0.17 0.12 0.56 0.05 −21.9 19.6 0.06 −15.7 13.4 0.05 0.08 0.21 0.11 0.07 0.07 0.30 0.06 0.38 0.11 0.11 −15.5 13.4 −23.0 20.4 −20.2 −14.5 −17.8 −23.6 −19.2 −11.1 −15.8 −22.3 −17.9 18.0 12.9 16.2 21.6 17.0 9.2 13.6 20.3 16.0 −16.1 14.3 −12.8 −20.6 10.7 18.4 −12.4 10.2 −15.7 −10.3 −18.2 13.8 8.4 17.1 −12.2 −21.2 10.2 20.1 0.08 0.07 0.11 0.12 0.07 0.10 0.07 0.06 0.11 Mn/Sr 25.8 −22.2 −20.2 0.20 0.12 0.04 0.07 0.09 Fe (ppm) 87 Sr/86 Sr 192 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Appendix A (Continued ) Sample number EL-35 EL-35A EL-36 El-37 El-38 EL-39 EL-40 EL-41 EL-42 EL-42A EL-43 EL-44 EL-45 EL-46 EL-47 EL-48 EL-49 EL-50M EL-50C El-51 EL-52 EL-53 EL-54 El-55 EL-56 EL-57 EL-58 EL-59 EL-60 EL-61 El-62 EL-63 EL-64 EL-65 Height above base (m) 42.0 42.0 43.5 44.3 45.0 47.0 48.1 49.6 51.1 51.2 51.6 53.1 54.6 55.6 56.5 58.4 59.9 61.7 61.7 63.2 65.5 66.3 66.9 68.1 69.6 71.1 72.6 74.1 75.6 77.4 78.9 80.4 81.9 83.4 Mineral/ rock ␦13 C, ‰ V-PDB ␦18 O, ‰ V-PDB TOC, mg C/g sample ␦13 CTOC , ‰ V-PDB δ Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite −2.1 −2.0 −1.6 −2.1 −1.6 −2.1 −1.9 −2.1 −2.0 −1.6 −2.1 −1.8 −1.9 −2.1 −2.2 −1.7 −2.0 −2.0 −1.7 −1.7 −1.7 −1.5 −2.1 −2.1 −1.9 −1.8 −1.5 −2.2 −2.1 −2.1 −1.7 −2.0 −2.1 −2.3 −14.4 −15.4 −14.0 −14.7 −13.5 −13.7 −14.0 −13.6 −14.9 −14.1 −14.4 −13.7 −14.3 −14.3 −14.1 −14.3 −15.8 −14.2 −14.4 −14.9 −13.5 −14.4 −13.9 −15.8 −15.7 −15.9 −14.1 −15.6 −15.8 −15.5 −15.5 −15.0 −13.9 −15.9 0.05 0.10 0.05 0.20 0.04 0.06 0.06 0.06 0.08 0.07 0.05 0.06 0.08 −14.4 12.3 −12.0 10.3 −12.5 10.9 −11.4 9.3 −8.3 6.5 −10.0 −8.7 8.3 6.7 −23.1 21.4 −8.6 6.7 −8.4 6.8 −22.4 −15.7 20.3 13.5 −12.9 10.9 0.05 0.11 0.09 0.11 0.06 0.16 0.06 0.05 0.06 0.08 0.06 0.06 0.14 0.17 0.20 0.15 0.07 0.05 Quirke Lake (Kerr-McGee Corp. Drill core 150/1, Bouck Township 150, 46◦ 27 41 , 82◦ 36 04 ) Dolostone Member DM-1722 524.9 m. calcite −2.1 −12.9 0.11 DM-1735 528.8 m. dolomite −1.9 −14.0 0.24 −23.9 DM-1748 532.8 m. calcite −2.5 −16.4 0.49 −23.3 DM-1749 533.1 m. calcite −2.6 −16.7 0.1 −19.4 DM-1755 A 534.9 m. calcite −2.7 −16.5 0.11 DM-1765 538.0 m. dolomite −2.1 −12.3 0.23 −26.1 DM1772 540.1 m. calcite −2.2 −15.1 0.12 −14.8 DM-1782 6 543.2 m. calcite −2.9 −15.5 0.14 −19.1 DM-1791 545.9 m. calcite −2.2 −13.0 0.1 DM-1796 547.4 m. calcite −2.5 −12.7 0.14 −17.0 DM-1806 550.5 m. calcite −2.3 −12.0 0.29 −24.7 DM-1816 553.5 m. dolomite −2.2 −13.4 −13.7 DM-1829 557.5 calc. mud-ne 0.15 −23.6 DM-1842 561.4 m. calcite −3.0 −15.9 0.45 −25.4 DM-1852 564.5 m. dolomite −2.5 −14.0 0.36 −24.8 DM-1862 567.5 m. calcite −3.1 −15.5 0.3 −22.2 DM-1872 570.6 m. dolomite −2.2 −12.9 0.08 −27.5 DM-1882 573.6 m. dolomite −2.7 −15.9 0.33 −21.3 DM-1892 576.7 m. calcite −3.3 −13.9 0.26 −14.3 DM-1902 579.7 m. dolomite −2.3 −12.1 0.23 −24.3 DM-1912 4 582.8 Calcite −3.1 −16.5 0.25 −23.2 Fe (ppm) Mn (ppm) Sr (ppm) Mn/Sr 196 199 170 164 181 164 197 182 241 335 335 263 333 242 326 223 144 138 0.6 0.6 0.6 0.5 0.7 0.5 0.9 1.3 1.7 155 249 161 267 161 188 169 213 258 237 188 265 196 237 211 98 222 302 0.7 1 0.6 1.4 0.7 0.9 1.7 1.0 0.9 257 264 287 234 159 210 251 268 234 227 194 143 126 182 357 290 313 299 257 257 1.3 1.8 2.3 1.3 0.4 0.7 0.8 0.9 0.9 0.9 289 226 399 226 92 140 1.3 2.5 2.8 46 67 46 54 115 tr. tr. 75 52 41 27 19 22.8 45.4 18.3 15.2 9.2 14.5 22.4 11.5 1036 3030 843 817 1060 895 696 760 2304 2026 1139 1584 22.4 22.2 19.2 25.3 18.6 11.0 22.0 20.1 2090 2792 934 1882 1023 1575 3370 2075 112 tr. 122 tr. tr. tr. 48 367 18.7 22.0 20.7 16.8 23.9 12.6 16.2 10.1 44.5 49.6 42.9 82.3 7.7 70.5 5.6 87 Sr/86 Sr A. Bekker et al. / Precambrian Research 137 (2005) 167–206 193 Appendix A (Continued ) Sample number DM1917 DM-1919 DM-1929 DM-1934 6 DM-1944 6 DM-1959 DM-1968 Height above base (m) Mineral/ rock ␦13 C, ‰ V-PDB ␦18 O, ‰ V-PDB 584.3 584.9 588.0 589.5 592.5 597.1 599.8 m. dolomite m. dolomite m. dolomite m. calcite m. calcite m. dolomite m. dolomite −3.9 −2.5 −3.3 −3.4 −3.6 −4.0 −3.5 −17.7 −12.9 −15.5 −16.2 −15.7 −16.2 −16.1 Siltstone Member DM1977 602.6 DM1990 606.6 DM2000 609.6 DM2010 612.6 DM2020 612.6 DM2030 618.7 DM2040 621.8 DM2050 624.8 DM2060 627.9 DM2080 634.0 DM2111 643.4 DM2127 648.3 DM2138 651.7 DM2148 654.7 DM2151 10 655.6 DM2183 665.4 DM2193 668.4 DM2203 671.5 DM2210 673.6 DM2223 677.6 DM2227 678.8 DM2240 682.8 ∗ DM2250 685.8 DM2283 695.9 DM2302 701.6 DM2308 7 703.5 Calc. si-ne Graywacke Calc. si-ne Calc. si-ne Calc. si-ne Mud-ne Calc. si-ne Mud-ne Calc. si-ne Calc. si-ne Calc. si-ne Calc. si-ne Calc. si-ne Calc. si-ne Calc. si-ne Mud-ne Calc. si-ne Calc. si-ne Calc. si-ne Mud-ne Mud-ne Calc. si-ne Calc. si-ne Calc. si-ne −3.0 Calc. si-ne Calc. si-ne Limestone Member DM2317 4 706.2 DM2318 8 706.5 DM2323 3 708.1 DM2328 709.6 DM2332 9 710.8 DM2339 712.9 DM2341 713.5 DM2343 9 714.1 Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite −2.9 −3.5 −1.5 −2.2 −1.8 −2.9 −2.4 −2.7 −16.5 −18.3 −14.3 −18.8 −17.1 −19.3 −18.8 −19.7 Denison Mine, Quirke Lake Limestone Member QE-1A 0.0 QE-2 6.5 QE-3 7.5 QE-4 9.0 QE-5 10.5 QE-6 17.0 QE-7 22.0 QE-8 27.0 QE-9 30.0 QE-10 36.0 QE-11 42.0 Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite −1.4 −3.9 −2.1 −1.2 −1.0 −1.4 −0.9 −1.4 −1.7 −1.1 −1.7 −19.3 −20.1 −20.5 −18.7 −19.2 −19.2 −18.9 −18.0 −17.2 −15.5 −16.4 −18.0 TOC, mg C/g sample ␦13 CTOC , ‰ V-PDB δ 0.12 0.1 0.59 0.59 0.27 0.49 −21.0 18.6 −23.9 −24.9 −22.8 −22.6 20.5 21.3 18.8 19.2 0.20 0.21 0.25 0.19 0.17 0.41 1.23 0.11 0.63 0.21 0.10 0.15 0.13 0.06 0.14 0.12 0.14 0.12 0.18 0.07 0.09 0.07 0.11 0.35 0.23 0.15 −23.4 −26.8 −15.0 −17.6 −24.7 −12.7 −18.2 −16.1 −13.8 −17.7 −17.8 −17.0 −11.5 −15.1 −19.1 −14.1 −16.7 −17.3 −24.4 −11.7 −15.8 −17.2 −12.8 0.11 0.14 0.06 0.16 0.19 0.07 0.09 0.16 Fe (ppm) Mn (ppm) Sr (ppm) 3086 589 399 1265 744 792 83 tr. tr. 115 tr. 63 475 tr. Mn/Sr 37.0 11.0 12.6 −11.9 −15.0 −11.4 −15.9 7.9 14.4 1342 106 12.6 −21.9 −29.5 20.2 26.6 2162 210 10.3 −17.3 14.6 1767.0 3168.8 1759.1 1830.6 1796.0 1346.5 1169.2 625.4 516.4 332.7 307.2 308.8 413.6 117.6 241.5 250.4 85.3 82.3 117.7 146.4 62.2 174.5 5.7 7.7 15.0 7.6 7.2 15.8 14.2 5.3 3.5 5.3 1.8 0.60 0.11 0.26 −31.6 −21.4 −21.2 30.6 20.0 20.3 0.16 0.21 −19.6 −20.2 17.9 19.1 87 Sr/86 Sr 194 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Appendix A (Continued ) Mineral/ rock ␦13 C, ‰ V-PDB ␦18 O, ‰ V-PDB Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite −1.5 −1.4 −1.5 −2.1 −1.4 −1.5 −1.5 −1.5 −1.3 −2.9 −2.0 −1.0 −15.2 −15.4 −15.1 −15.2 −14.5 −14.5 −14.9 −15.2 −14.8 −17.1 −15.6 −11.9 Dolostone Member QE-24 0.0 QE-25 3.0 QE-26 5.0 QE-27 8.0 QE-29 16.5 QE-30 23.5 QE-31 24.0 QE-32 24.5 QE-33 27.1 QE-34 28.4 QE-35 30.1 QE-37 39.0 QE-38 50.4 QE-39 52.7 QE-40 57.1 QE-41 59.0 QE-42 60.0 QE-44 65.0 Esp-1-3∗ Esp-1-9∗ Esp-2-1∗ Esp-3-1∗ Esp-5-1∗ Calcite Dolomite Calcite Calcite Calcite Calcite Dolomite Dolomite Calcite Dolomite Calcite Calcite Calcite Dolomite Dolomite Dolomite Dolomite Dolomite Dolomite Dolomite Dolomite Dolomite Dolomite −3.0 −2.5 −2.7 −2.7 −2.0 −2.8 −2.4 −2.9 −2.3 −2.5 −2.8 −2.7 −2.3 −2.2 −2.0 −2.3 −2.5 −1.8 −2.5 −2.2 −2.1 −2.1 −2.5 −18.5 −14.3 −13.0 −13.9 −14.8 −16.4 −15.9 −18.2 −14.4 −17.8 −16.0 −15.6 −17.7 −13.5 −13.4 −16.5 −15.1 −14.9 −16.5 −16.7 −16.0 −15.7 −13.6 Moose Point MP-4 MP-7 MP-9 MP-10 MP-11 MP-12 MP-27 MP-28 MP-29 MP-30 MP-31 MP-32 MP-33 MP-34 MP-35 MP-36 MP-38 MP-39 MP-40 Cal. si-ne Dolomite Calcite Dolomite Dolomite Dolomite Calcite Calcite Calc. si-ne Calc. si-ne Calcite Calcite Calcite Calcite Calcite Calcite Calc. si-ne Calcite Calcite −3.2 −3.1 −3.1 −2.8 −3.7 −3.6 −4.2 −12.9 −13.5 −14.5 −14.1 −15.1 −16.4 −18.2 Sample number QE-12 QE-12A QE-13 QE-13A QE-14 QE-15 QE-16 QE-17 QE-18 QE-20 QE-21 QE-22 Height above base (m) 46.0 50.8 52.0 52.0 53.1 54.3 55.3 56.3 57.0 80.0 79.0 77.8 2.1 5.1 7.1 8.1 9.1 10.1 140.0 148.5 155.7 159.0 170.6 178.8 182.2 191.0 194.0 198.7 204.9 207.4 212.3 −5.1 −6.7 −3.7 −3.5 −3.5 −3.9 −5.2 −3.9 −3.7 −18.1 −18.0 −17.8 −15.9 −16.8 −16.1 −17.6 −15.8 −15.4 TOC, mg C/g sample 0.18 0.04 0.20 0.03 0.19 ␦13 CTOC , ‰ V-PDB δ −18.4 16.3 −14.6 13.1 Fe (ppm) 0.17 0.36 0.09 0.11 43100 38000 37800 34900 0.07 0.06 0.13 −16.0 −16.8 −21.4 13.6 18.3 0.46 −17.4 13.7 0.08 0.13 −17.3 −18.4 0.07 0.09 Mn (ppm) Sr (ppm) Mn/Sr 282.4 191.7 220.8 199.3 179.0 210.4 249.6 281.4 183.4 641.7 119.0 1882.7 161.0 184.6 194.3 209.6 174.7 172.1 156.8 180.8 235.8 32.4 257.1 tr. 1.8 1.0 1.1 1.0 1.0 1.2 1.6 1.6 0.8 19.8 0.5 677.3 1882.7 3566.1 2642.0 3762.1 2099.3 3819.8 1238.0 2963.2 3360.7 5171.9 661.7 1534.7 4077.2 4707.0 2294.5 2649.9 606.0 118.0 tr. 65.8 38.1 113.1 307.9 55.6 25.7 55.5 22.8 101.0 tr. tr. 56.1 67.2 tr. 48.2 tr. 5.7 5825.8 5289.0 4685.4 5044.0 186.2 182.1 153.9 131.9 5935 5005 5882 4931 5049 2609 1852 1 10 11 19 27 233 125 8551.3 520.5 515.6 260.1 187.4 11.2 14.8 1756 2771 2927 2970 3321 3358 2756 4261 5477 27 162 214 159 189 173 81 272 237 65.3 17.1 13.7 18.6 17.5 19.5 34.1 15.6 23.1 54.2 69.3 33.3 6.8 68.7 48.2 53.4 147.1 51.2 72.7 70.1 55.0 31.3 29.0 30.4 38.2 87 Sr/86 Sr A. Bekker et al. / Precambrian Research 137 (2005) 167–206 195 Appendix A (Continued ) Height above base (m) Mineral/ rock ␦13 C, ‰ V-PDB ␦18 O, ‰ V-PDB 216.4 227.3 238.6 241.0 265.0 267.8 270.0 274.1 278.1 297.0 Calcite Calc. si-ne Calcite Calcite Calc. si-ne Calcite Calcite Calcite Calcite Calc. si-ne −3.7 −15.9 −3.8 −3.8 −3.3 −3.0 −3.0 −3.0 −3.2 −3.1 −15.9 −16.0 −15.3 −15.2 −12.6 −15.9 −15.4 −16.4 Serpent Formation SE-1-1∗ SE-2.3∗ SE-2.4a∗ SE-2.4b∗ SE-2.5a∗ SE-2.6∗ ES-1 ES-5 SE.2.1 SE.2.2 SA-98-1-A Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite −2.9 0.7 −6.3 −1.0 −0.7 0.7 −8.0 −0.1 −2.5 0.9 −0.5 −20.8 −21.0 −21.0 −20.9 −20.7 −20.3 −16.2 −7.3 −20.2 −19.8 −20.3 Gowganda Formation IS-GO Calcite −8.8 −9.7 Vagner Formation BE-C-95-1∗ BE-C-95-2∗ BE-C-95-3∗ BE-C-95-4∗ BE-C-95-7∗ BE-C-95-8∗ BE-C-95-9∗ BE-C-95-10a∗ BE-C-95-10b∗ BE-C-95-11∗ BE-C-95-12∗ BE-C-95-20∗ BE-C-95-21∗ BE-C-95-24∗ BW99.1 BW99.2 BW99.3 BW99.4 BW99.5 BW99.6 BW99.7 BW99.8 BW99.9 BW99.10 BW99.11 5 6 8 Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite Calcite −2.3 −2.5 -2 −2.5 −1.9 −1.5 −2.6 −2.3 −2.2 −1.9 −1.7 −1.7 −2.2 −2.2 −2.3 −2.1 −2.4 −2.0 −2.6 −2.1 −2.0 −2.0 −1.9 −2.1 −1.9 −2.9 −3.6 −3.8 −19.8 −19.4 −20.9 −18.2 −18.8 −19.9 −16.4 −17.4 −19.6 −19.4 −18.9 −19.0 −17.7 −16.7 −15.3 −15.7 −16.1 −15.4 −16.7 −15.9 −14.5 −16.8 −15.8 −17.9 −18.3 −18.6 −18.1 −17.2 Sample number MP-41 MP-43 MP-46 MP-47 MP-52 MP-52A MP-53 MP-54 MP-54A MP-55B TOC, mg C/g sample ␦13 CTOC , ‰ V-PDB 0.25 −16.0 δ Fe (ppm) 0.25 504 Mn (ppm) Sr (ppm) Mn/Sr 4655 262 17.8 3580 4980 915 2122 2534 1577 1854 1677 231 188 4 113 28 106 113 8 15.5 26.4 212.4 18.7 91.5 14.9 16.4 201.9 998 2374 4170 4406 4227 35 36 463 254 305 28.7 66.9 9.0 17.3 13.8 1055 128 8.3 87 Sr/86 Sr 0.05 −17.7 15.3 0.14 0.23 0.09 −20.4 −32.8 −19.6 92 17.9 94 30.9 6306 18.1 496 443 1461 117 187 181 4.2 2.4 0.708075 8.1 0.1 0.08 −12.6 −14.8 393 520 831 630 217 187 217 187 1.8 0.715257 2.8 3.8 3.4 0.08 0.07 0.3 −18.8 −14.1 −19.1 10.3 288 12.6 390 2228 1829 17.1 11.9 16.9 748 529 295 352 453 523 527 513 571 568 468 408 480 343 548 566 145 59 61 59 48 59 95 125 161 73 52 62 21 42 36 3.7 5.0 5.8 7.7 10.9 8.9 5.4 4.6 3.5 6.4 7.8 7.8 16.3 13.0 15.7 196 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Appendix A (Continued ) Sample number Height above base (m) Mineral/ rock ␦13 C, ‰ V-PDB ␦18 O, ‰ V-PDB Bottle Creek Formation, Sierra Madre SE1/4, SE1/4, SE1/4, Sec. 17, T14N, R87W BS-1 6.5 Calcite −2.2 BS-2 5 Calcite −2.4 BS-3 0 Calcite −3.4 BS-4A 13.5 Calcite −3.1 BS-4B 13.5 Calcite −2.8 BS-4C 13.5 Calcite −3.6 BS-5 19 Calcite −4.8 −14.2 −13.9 −14.5 −14.1 −12.5 −15.1 −16.1 SE1/4, SE1/4, SE1/4, Sec. 16, T14N, R87W BS-6 Calcite −2.5 −13.4 TOC, mg C/g sample ␦13 CTOC , ‰ V-PDB δ Fe (ppm) Mn (ppm) Sr (ppm) Mn/Sr 87 Sr/86 Sr Samples marked with the asterisk were analyzed at the Geological Survey of Finland using whole rock powders, for the rest of samples C, O, and Sr isotope analyses of carbonates were done at the University of Maryland and Mn and Sr contents were measured at VPI & SU. TOC content was measured at the University of Maryland; carbon isotope values of organic carbon were measured at the Mountain Mass Spectrometry. M: cement, C: clast. Appendix B. Detailed description of studied sections B.1. Geneva Lake In the Geneva Lake area, the Espanola Formation rests unconformably on Archean basement or, locally, conformably on conglomerates of the Bruce Formation. Its thickness ranges from less than 15 m to over 90 m with an average around 45 m (Card and Innes, 1981). The upper contact with the Serpent Formation is conformable and gradational (Card and Innes, 1981). The lower part of the Espanola Formation includes a 15–60 m thick limestone-dolomite member with thin graded beds of sandstone and siltstone with carbonate cement and mudchips, shales, and carbonates (Locations 2 and 3 in Fig. 15). Intraformational flat-pebble conglomerates with imbricated rounded pebbles of laminated carbonate up to 2–3 cm in diameter are locally present (Loc. 2; Fig. 6). Graded sandstones have erosional base and cross-bedding, whereas overlying shales have convolute bedding and other features of soft-sediment deformation. The lower part of the section in Locality 3 is marked by a basal conglomerate overlain by white to pink carbonates with shale seams, stratiform stromatolites (Fig. 9A) and convolute bedding (Location 3, Fig. 6). Shaly carbonate with sandstone caps the lower part of the section. The upper part of the Espanola Formation in this area is 15–30 m thick and consists of laminated siltstones and sandstones with carbonate cement and thin carbonate beds (Location 1 in Fig. 15). Parallel, cross-, and lensoidal bedding and graded beds are common in the upper part. B.2. Echo Lake The lower contact of the Espanola Formation is not exposed in the sampled section (Fig. 16). Here, the Gowganda Formation disconformably overlies the limestone member of the Espanola Formation with the higher members and the Serpent Formation missing. Three facies are interlayered in this section (Fig. 7): (1) intraformational conglomerate with carbonate clasts up to 8 cm in size (Fig. 9B); (2) finely interlayered carbonate, marl, and shale; and (3) massive carbonate with thin shale beds. Some intraformational conglomerates include thin beds of shaly and massive carbonates. The shale-dominated facies have convolute bedding and include channelized beds of intraformational conglomerate and lenses and beds of massive carbonates ranging in thickness from 0.5 to 60 cm. There are also meter-size slumped blocks in these facies. Massive carbonate beds increase in thickness and abundance upsection. This section displays evidence for synsedimentary faulting and folding and was later intricately folded forming sheath folds. A. Bekker et al. / Precambrian Research 137 (2005) 167–206 197 Fig. 15. Geological map of the Huronian Supergroup in the Geneva Lake area (modified from Card and Innes, 1981) with sampled localities shown. See Fig. 2 for location of the Geneva Lake area. B.3. Quirke Lake Two sections of the Espanola Formation were logged and sampled in the Quirke Lake area; one is the Kerr-McGee drill core 150/1 from the Bouck (150) Township and the other is along the Quirke Lake at Denison Mine (Fig. 17). The contact with the diamictite of the Bruce Formation in the drill core is sharply overlain by thin shale interval (not shown on Fig. 4) succeeded by white limestones with thin shaly lenses and beds (Fig. 4). The limestone interval contains thin beds of limestone pebbles in a siliciclastics-rich matrix with soft sediment deformation, injection, and flame structures indicating high sedimentation rates and slope instability. A thin bed with spherules up to 2 mm in diameter filled with megaquartz with undulose radial extinction occurs in the upper part of this member. These spherules are similar to those formed after anhydrite in Phanerozoic and Neoproterozoic successions (Milliken, 1979). The overlying siltstone member consists of dark argillaceous siltstone with graded beds and intraformational conglomerates with limestone and shale pebbles, and thin carbonate beds and lenses. Graded beds have sharp erosional lower contact. This member also contains soft sediment deformation and injection structures. The upper member is dolomitic in the field but in the drill core is identified by higher carbonate content in silty beds within thick mudstone 198 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Fig. 16. Geological map of the Echo Lake area (modified from Bennett, 1982). The Marble Point section was sampled in this study. See Fig. 2 for location of the Echo Lake area. A. Bekker et al. / Precambrian Research 137 (2005) 167–206 199 Fig. 17. Geological map of the Quirke Lake area (modified from Robertson, 1968). See Fig. 2 for location of the Quirke Lake area. Sampled section is at the shore of Quirke Lake between Quarry Bay and No. 2 Shaft and sampled drill hole 150/1 is along the southern border of the map. 200 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 intervals containing soft sediment deformation and injection structures. This member has graded beds of sandy carbonate, carbonate, and shale and contains shale and carbonate intraclasts. Sandy carbonates have basal scour surfaces and cross-bedding whereas overlying carbonates contain rare fenestral textures and peloids. The limestone member in the outcrop section starts with a 1.5–3.0 m thick shale overlain by thin-bedded grey siltstone and limestone (Fig. 9). Limestone beds are thinly laminated and contain soft sediment deformation structures. The overlying siltstone member contains ball-and-pillow structures and Bouma beds. The contact between the limestone and siltstone members ranges from gradational to erosional. Conglomerate beds with closely packed pebbles of carbonate, rare granite, and chert occur above this contact. Intraformational breccias are common at the top of the member. The overlying dolostone member consists of interlayered ferruginous dolomite, calcareous siltstone, and limestone. The member contains flat-pebble conglomerate beds that resemble “beach rosettes” suggesting deposition in oscillatory flow. Other sedimentary structures of the dolostone member are small-scale and sigmoidal cross-bedding, abundant wave ripples with internal laminations (Fig. 9D), convolute bedding, thin laminations, rare ball-and-pillow structures, and graded beds. The contact with the underlying member is gradational. Several meter thick upward-deepening cycles contain wave-rippled dolomite at the base grading to massive dolomite and overlain by calcareous siltstone. The uppermost part of the member contains mudcracks, symmetrical and ladder ripples, and wave ripples indicating a shallow-water depositional setting. The contact with the overlying Serpent Formation is conformable in this area. Fig. 18. Geological map of Moose Point (modified from Bernstein and Young, 1990) with sampled profile shown. See Fig. 2 for location of the Moose Point area. A. Bekker et al. / Precambrian Research 137 (2005) 167–206 B.4. Moose Point near Whitefish Falls The Espanola Formation in this section (Fig. 18) shows middle to upper greenschist facies grade metamorphism (Card, 1978). The limestone member is 6–17 m thick (Bernstein and Young, 1990) and sharply overlies diamictite of the Bruce Formation (Fig. 5). The basal part contains laminated sandy shales overlain by interlayered grey limestones and dolostones, silty limestones and dolostones, and calcareous siltstones. The following sedimentary structures indicate a lowenergy, deep-water environment with episodic highenergy fluxes due to storm events, turbidity currents, or slope failure resulting from fluidization: laminations, rare climbing ripples, graded beds, convolute and crossbedded layers, and synsedimentary folds. A 0.5-m thick conglomerate present in the upper part consists of imbricated clasts up to 30 cm in size of non-calcareous to slightly calcareous shales, calcareous siltstones, and silty limestones. The gradationally overlying 230–310 m thick siltstone member consists of cross-bedded and laminated calcareous siltstone and sandy shale with ball-andpillow structures (Fig. 9C) and lenses and beds of impure carbonates. Flaser and wavy bedding and symmetrical and asymmetrical wave ripples with internal cross-laminations are present in shales at the top of sandy beds. The member also contains syndepositional faults overlain by disrupted beds. The siltstone member contains at the base a 35 m thick intraformational breccia-conglomerate with pebbles and boulders up to 1 m in size of calcareous and non-calcareous siltstone and sandstone. The 212–238 m thick heterolithic member consists of upward-fining cycles of coarse-grained sandstone, siltstone, shale, and thin impure dolostone and limestone. Sandstones have trough and large low-angle cross-bedding, ball-and-pillow structures, herringbone cross-stratification and climbing, symmetrical, and asymmetrical ripples. Small channels with granite and quartzite pebbles up to 5 cm in diameter occur in this member. Fine-grained rocks have convolute bedding, thin laminations, starved ripples, and slump and load structures. The Espanola Formation at Moose Point records a regressive cycle with the limestone member deposited below the wave base, the siltstone member above the wave base, and the heterolithic member 201 influenced by tidal processes (Bernstein and Young, 1990). References Alibert, C., McCulloch, M.T., 1993. Rare element and neodymium isotopic compositions of the banded iron-formations and associated shales from Hamersley, western Australia. Geochim. Cosmochim. Acta 57, 187–204. Allen, P.A., Hoffman, P.F., 2005. Extreme winds and waves in the aftermath of a Neoproterozoic glaciation. Nature 433, 123– 127. Andrews, A.J., Masliwec, A., Morris, W.A., Owsiacki, L., York, D., 1986. The silver deposits at Cobalt and Gowganda, Ontario. II: An experiment in age determinations employing radiometric and paleomagnetic measurements. Can. J. Earth Sci. 23, 1507– 1518. Aspler, L.B., Chiarenzelli, J.R., 1997. Initiation of ∼2.45–2.1 Ga intracratonic basin sedimentation of the Hurwitz Group, Keewatin Hinterland, Northwest Territories, Canada. Precambrian Res. 81, 265–297. Bau, M., Beukes, N.J., Romer, R.L., 1998. Increase of oxygen in the Earth’s atmosphere and hydrosphere between −2.5 and −2.4 Ga B.P. Min. Mag. 62A, 127–128. Baumgartner, L.P., Valley, J.W., 2001. Stable Isotope Transport and Contact Metamorphic Fluid Flow. Rev. Mineral. Geochem. 43, 415–467. Bekker, A., 1998. Chemostratigraphy and climatostratigraphy of the Paleoproterozoic Snowy Pass Supergroup, Wyoming and its application for correlation with other sequences in North America. M.S. thesis, University of Minnesota, Duluth, 104 pp. Bekker, A., Eriksson, K.A., 2003. Paleoproterozoic drowned carbonate platform on the southeastern margin of the Wyoming Craton: a record of the Kenorland breakup. Precambrian Res. 120, 327–364. Bekker, A., Holland, H.D., Wang, P.-L., Rumble, D., Stein, H.J., Hannah, J.L., Coetzee, L.L., Beukes, N.J., 2004a. Dating the rise of the atmospheric oxygen. Nature 427, 117–120. Bekker, A., Holmden, C., Patterson, W., Eglington, B.L., Coetzee, L.L., Beukes, N.J., 2004. Chemostratigraphy of early Paleoproterozoic carbonates of South Africa. GSA Abstracts with Programs 36(6), A-341. Bekker, A., Karhu, J.A., Bennett, G., 1996. Lomagundi event in North America: application for stratigraphic correlation. Geological Society of America, Annual Meeting, Abstracts with Programs 28(7), A-229. Bekker, A., Karhu, J.A., Eriksson, K.A., Kaufman, A.J., 2003. Chemostratigraphy of Paleoproterozoic carbonate successions of the Wyoming Craton: tectonic forcing of biogeochemical change? Precambrian Res. 120, 279–325. Bekker, A., Kaufman, A.J., Karhu, J.A., Beukes, N.J., Swart, Q.D., Coetzee, L.L., Eriksson, K.A., 2001. Chemostratigraphy of the Paleoproterozoic Duitschland Formation, South Africa: Implications for coupled climate change and carbon cycling. Am. J. Sci. 301, 261–285. 202 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Bekker, A., Ono, S., Rumble, D., 2005. Low atmospheric pO2 in the aftermath of the oldest Paleoproterozoic glaciation. Astrobiology 5(2), 244. Bekker, A., Sial, A.N., Karhu, J.A., Ferrerira, V.P., Noce, C.M., Kaufman, A.J., Romano, A.W., Pimentel, M.M., 2003a. Chemostratigraphy of carbonates from the Minas Supergroup, Quadrilátero Ferrı́fero (Iron Quadrangle), Brazil: a stratigraphic record of early Proterozoic atmospheric, biogeochemical and climatic change. Am. J. Sci. 303, 865–904. Bennett, G., 1982. Geology of the Two Horse Lake Area, District of Algoma. Ontario Geological Survey, Report, vol. 210, 63 pp. Bernstein, L., Young, G.M., 1990. Depositional environments of the Early Proterozoic Espanola Formation, Ontario, Canada. Can. J. Earth Sci. 27, 539–551. Beukes, N.J., Dorland, H.C., Gutzmer, J., 2002. Pisolitic ironstone and ferricrete in the 2.22–2.4 Ga Timeball Hill Formation, Transvaal Supergroup: implications for the history of atmospheric oxygen. GSA Abstracts with Programs 34(6), A-283. Beukes, N.J., Dorland, H., Gutzmer, J., Nedachi, M., Ohmoto, H., 2002a. Tropical laterites, life on land, and the history of atmospheric oxygen in the Paleoproterozoic. Geology 30, 491–494. Beukes, N.J., 1983. Paleoenvironmental setting of iron-formations in the depositional basin of the Transvaal Supergroup, South Africa. In: Trendall, A.F., Morris, R.C. (Eds.), Iron-Formations: Facts and Problems. Elsevier, Amsterdam, Netherlands, pp. 131–209. Bowring, S.A., Myrow, P.M., Landing, E., Ramezani, J., Condon, D., Hofmann, K.H., 2003. Geochronologic constraints on Neoproterozoic glaciations and the rise of metazoans. GSA Abstracts with Programs 35(6), 516. Buick, I.S., Uken, R., Gibson, R.L., Wallmach, T., 1998. High␦13 C Paleoproterozoic carbonates from the Transvaal Supergroup, South Africa. Geology 26, 875–878. Cameron, E.M., 1983. Evidence from early Proterozoic anhydrite for sulfur isotopic partitioning in Precambrian oceans. Nature 304, 54–56. Cameron, E.M., 1982. Sulphate and sulphate reduction in early Precambrian oceans. Nature 296, 145–148. Card, K.D., 1978. Metamorphism of the Middle Precambrian (Aphebian) rocks of the eastern southern province, Metamorphism in the Canadian Shield. Geol. Surv. Can., Paper 78-10, pp. 269–282. Card, K.D., Innes, D.G., 1981. Geology of the Benny Area, District of Sudbury, Ontario Geological Survey, Report 206, 117 pp. Card, K.D., Innes, D.G., Debicki, R.L., 1977. Stratigraphy, sedimentology, and petrology of the Huronian Supergroup in the SudburyEspanola Area, Ontario Div, Mines, GS16, 99 pp. Chamberlain, K.R., 1998. Medicine Bow orogeny: Timing of deformation and model of crustal structure produced during continentarc collision, ca. 1.78 Ga, southeastern Wyoming. Rocky Mount. Geol. 33 (2), 259–277. Chandler, F.W., 1989. Lower Proterozoic sabkha-related copper mineralization, paleoenvironment and diagenesis, Cobre Lake, Ontario. In: Boyle, R.W., Brown, A.C., Jefferson, C.W., Jowett, E.C., Kirkham, R.V. (Eds.), Sediment-Hosted Stratiform Copper Deposits. Geol. Ass. Canada Special Paper 36, pp. 225–244. Chandler, F.W., 1988. Diagenesis of sabkha-related, sulphate nodules in the Early Proterozoic Gordon Lake Formation, Ontario, Canada. Carbon. Evap. 3 (1), 75–94. Coetzee, L.L., 2001. Genetic stratigraphy of the Paleoproterozoic Pretoria Group in the Western Transvaal. M.Sc. thesis, Rand Afrikaans University, 212 pp. Cornell, D.H., Schütte, S.S., Eglington, B.L., 1996. The Ongeluk basaltic andesite formation in Griqualand West, South Africa: submarine alteration in a 2222 Ma Proterozoic sea. Precambrian Res. 79, 101–123. Corsetti, F.A., Kaufman, A.J., 2003. Stratigraphic investigations of carbon isotope anomalies and Neoproterozoic ice ages in Death Valley, California. GSA Bull. 115, 916–932. Cox, D.M., Frost, C.D., Chamberlain, K.R., 2000. 2.01-Ga Kennedy dike swarm, southeastern Wyoming: record of a rifted margin along the southern Wyoming province. Rocky Mount. Geol. 35 (1), 7–30. Donnadieu, Y., Godderis, Y., Ramstein, G., Nédélec, A., Meert, J., 2004. A ‘snowball Earth’ climate triggered by continental breakup through changes in runoff. Nature 428, 303–306. Dorland, H.C., 1999. Paleoproterozoic laterites, red beds and ironstones of the Pretoria Group with reference to the history of atmospheric oxygen. Unpubl. M.Sc. thesis. Rand Afrikaans University, Johannesburg, South Africa, 147 pp. Eisbacher, G.H., 1970. Contemporaneous faulting and clastic intrusions in the Quirke Lake Group, Elliot Lake, Ontario. Can. J. Earth Sci. 7, 215–225. Emanuel, K.A., 1987. The dependence of hurricane intensity on climate. Nature 326, 483–485. Evans, D.A., Beukes, N.J., Kirschvink, J.L., 1997. Low-latitude glaciation in the Paleoproterozoic era. Nature 386, 262–266. Farquhar, J., Bao, H., Thiemens, M., 2000. Atmospheric influence of Earth’s earliest sulfur cycle. Science 289, 756–758. Fedo, C.M., Young, G.M., Nesbitt, H.W., 1997. Paleoclimatic control on the composition of the Paleoproterozoic Serpent Formation, Huronian Supergroup, Canada: a greenhouse to icehouse transition. Precambrian Res. 86, 201–223. Fralick, P.W., Miall, A.D., 1989. Sedimentology of the Lower Huronian Supergroup (Early Proterozoic), Elliot Lake area, Ontario, Canada. Sed. Geol. 63, 127–153. Goddéris, Y., Donnadieu, Y., Nédélec, A., Dupré, B., Dessert, C., Grard, A., Ramstein, G., François, L.M., 2003. The Sturtian ’snowball’ glaciation: fire and ice. Earth Plan. Sci. Lett. 211, 1–12. Graff, P.J., 1978. Geology of the Early Proterozoic Snowy Pass Supergroup, Sierra Madre, Wyoming-with chapters on Proterozoic regional tectonics and uraniferous quartz-pebble conglomerates. Ph.D. thesis, The University of Wyoming, Laramie, Wyoming, 85 pp. Grotzinger, J.P., Knoll, A.H., 1995. Anomalous Carbonate Precipitates: Is the Precambrian the Key to the Permian? Palaios 10, 578–596. Halverson, G.P., Hoffman, P.F., Schrag, D.P., Kaufman, A.J., 2002. A major perturbation of the carbon cycle before the Ghaub glaciation (Neoproterozoic) in Namibia: Prelude to snowball Earth? Geochem. Geophys. Geosyst. 3 (6), June. Hannah, J.L., Bekker, A., Stein, H.J., Markey, R.J., Holland, H.D., 2004. Primitive Os and 2316 Ma age for marine shale: implications for Paleoproterozoic glacial events and the rise of atmospheric oxygen. Earth Planetary Sci. Lett. 225, 43–52. A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Hattori, K., Krouse, H.R., Campbell, F.A., 1983. The start of sulfur oxidation in continental environments: about 2.2 × 109 years ago. Science 221, 549–551. Hayes, J.M., Kaplan, I.R., Wedeking, K.W., 1983. Precambrian organic geochemistry, preservation of the record. In: Schopf, J.W. (Ed.), The Earth’s Earliest Biosphere: Its Origin and Evolution. Princeton University Press, Princeton, N.J, pp. 93–134. Hoffman, P.F., 1991. Did the breakup of Laurentia turn Gondwanaland inside-out? Science 252, 1409–1412. Hoffman, P.F., Schrag, D.P., 2002. The snowball Earth hypothesis: testing the limits of global change. Terra Nova 14, 129–155. Hoffman, P.F., Schrag, D.P., 2000. Snowball Earth. Scientific Am. 282 (1), 68–75. Hoffman, P.F., Kaufman, A.J., Halverson, G.P., Schrag, D.P., 1998. A Neoproterozoic snowball Earth. Science 281, 1342–1346. Hofmann, H.J., Pearson, D.A.B., Wilson, B.H., 1980. Stromatolites and fenestral fabric in Early Proterozoic Huronian Supergroup, Ontario. Can. J. Earth Sci. 17, 1351–1357. Houston, R.S., 1993. Late Archean and Paleoproterozoic geology of southeastern Wyoming. In: Snoke, A.W., Steidtmann, J.R., Roberts, S.M. (Eds.), Geology of Wyoming, vol. 5. Geological Survey of Wyoming Memoir, pp. 78–116. Houston, R.S., Graff, P.J., 1995. Geologic Map of Precambrian rocks of the Sierra Madre, Carbon County, Wyoming, and Jackson and Routt Counties, Colorado. U.S. Department of the Interior U.S. Geol. Surv., Misc. Invest. Ser. I-2452. Houston, R.S., Karlstrom, K.E., 1992. Geological map of Precambrian metasedimentary rocks of the Medicine Bow Mountains, Albany and Carbon Counties, Wyoming. Unites States Geological Survey, Misc. Inv. Ser. Map I-2280, 1:50,000. Houston, R.S, Karlstrom, K.E., Graff, P.J., Flurkey, A.J., 1992. New stratigraphic subdivisions and redefinition of subdivisions of Late Archean and Paleoproterozoic metasedimentary and metavolcanic rocks of the Sierra Madre and Medicine Bow Mountains, southern Wyoming. U.S. Geological Survey Professional Paper 1520, 50 pp. Houston, R.S., Lanthier, L.R., Karlstrom, K.E., Sylvester, G., 1981. Paleoproterozoic diamictite of southern Wyoming. In: Hambrey, M.J., Harland, W.B. (Eds.), Earth’s Pre-Pleistocene Glacial Record. Cambridge University Press, New York, pp. 795–799. James, N.P., Narbonne, G.M., Kyser, T.K., 2001. Late Neoproterozoic cap carbonates: Mackenzie Mountains, northwestern Canada: precipitation and global glacial meltdown. Can. J. Earth Sci. 38, 1229–1262. Karhu, J.A., 1993. Paleoproterozoic evolution of the carbon isotope ratios of sedimentary carbonates in the Fennoscandian Shield. Geol. Surv. Finland Bull. 371. Espoo, 87 pp. Karhu, J.A., Holland, H.D., 1996. Carbon isotopes and the rise of atmospheric oxygen. Geology 24 (10), 867–870. Karlstrom, K.E., 1977. Geology of the Proterozoic Deep Lake Group, central Medicine Bow Mountains, Wyoming. M.S. thesis, University of Wyoming, Laramie, 116 pp. Karlstrom, K.E., Houston, R.S., 1984. The Cheyenne belt: analysis of a Proterozoic suture in southern Wyoming. Precambrian Res. 25, 415–446. Karlstrom, K.E., Houston, R.S., 1979. Stratigraphy of the Phantom Lake Metamorphic Suite and Deep Lake Group and a review of 203 the Precambrian tectonic history of the Medicine Bow Mountains. Contr. Geol. 17 (2), 111–133. Karlstrom, K.E., Houston, R.S., 1979a. Stratigraphy and uranium potential of Early Proterozoic metasedimentary rocks in the Medicine Bow Mountains, Wyoming. The Geological Survey of Wyoming, Report of Investigations No. 13, Laramie, Wyoming, 45 pp. Karlstrom, K.E., Flurkey, A.J., Houston, R.S., 1983. Stratigraphy and depositional setting of the Proterozoic Snowy Pass Supergroup, southeastern Wyoming: Record of a Paleoproterozoic Atlantictype cratonic margin. Geol. Soc. Am. Bull. 94, 1257–1274. Karlstrom, K.E., Houston, R.S., Flurkey, A.J., Coolidge, C.M., Kratochvil, A.L., Sever, C.K., 1981. A summary of the geology and uranium potential of Precambrian conglomerates in southeastern Wyoming, U.S. Dept. Energy GJBX-139-81. Grand Junction, Colorado, 537 pp. Kaufman, A.J., Knoll, A.H., 1995. Neoproterozoic variations in the C-isotopic composition of seawater: stratigraphic and biogeochemical implications. Precambrian Res. 73, 27–49. Kaufman, A.J., Xiao, S., 2003. High CO2 levels in the Proterozoic atmosphere estimated from analyses of individual microfossils. Nature 425, 279–282. Kaufman, A.J., Hayes, J.M., Knoll, A.H., Germs, G.J.B., 1991. Isotope compositions of carbonates and organic carbon from upper Proterozoic successions in Namibia: Stratigraphic variation and the effects of diagenesis and metamorphism. Precambrian Res. 49, 301–327. Kaufman, A.J., Knoll, A.H., Narbonne, G.M., 1997. Isotopes, ice ages, and terminal Proterozoic Earth history. Natl. Acad. Sci. Proc. 94, 6600–6605. Kempf, O., Kellerhals, P., Lowrie, W., Matter, A., 2000. Paleomagnetic directions in Late Precambrian glaciomarine sediments of the Mirbat Sandstone Formation, Oman. Earth Plan. Sci. Lett. 175, 181–190. Kennedy, M.J., 1996. Stratigraphy, sedimentology, and isotopic geochemistry of Australian Neoproterozoic postglacial cap dolostones: deglaciation ␦13 C excursions, and carbonate precipitation. J. Sed. Res. 66 (6), 1050–1064. Kirschvink, J.L., 1992. Late Proterozoic low-latitude global glaciation: the snowball earth. In: Schopf, J.W., Klein, C. (Eds.), The Proterozoic Biosphere. Cambridge University Press, Cambridge, pp. 51–52. Kirschvink, J.L., Gaidos, E.J., Bertani, L.E., Beukes, N.J., Gutzmer, J., Maepa, L.N., Steinberger, R.E., 2000. The Paleoproterozoic snowball Earth: extreme climatic and geochemical global change and its biological consequences. Proc. Natl. Acad. Sci. 97, 1400–1405. Krogh, T.E., Davis, D.W., Corfu, F., 1984. Precise U-Pb zircon and badelleyite ages for the Sudbury area. In: Pye, E.G., Naldrett, A.J., Giblin, P.E. (Eds.), The Geology and Ore Deposits of the Sudbury Structure. Ontario Geol. Surv. Spec. 1, pp. 431– 446. Krogh, T.E., Kamo, S.L., Bohor, B.F., 1996. Shock metamorphosed zircons with correlated U–Pb discordance and melt rocks with concordant protolith ages indicate an impact origin for the Sudbury structure. Earth Processes: Reading the isotopic code. Geophys. Monograph 95. Am. Geophys. Union, pp. 343–353. 204 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Long, D.G.F., 1987. Sedimentary framework of uranium deposits in the southern Cobalt embayment, Ontario, Canada, Proterozoic quartz-pebble conglomerates. International Atomic Agency, IAEA-TECDOC-427 Vienna, pp. 155–188. Long, D.G.F., 1976. The stratigraphy and sedimentology of the Huronian (lower Aphebian) Mississagi and Serpent Formations. Ph.D. thesis, University of Western Ontario, London, 291 pp. Long, D.G.F., Young, G.M., Rainbird, R.H., Fedo, C.M., 1999. Actualistic and non-actualistic Precambrian sedimentary styles: examples from the Proterozoic, north shore of Lake Huron. Field Trip B5 Guidebook, Geol. Ass. Can. -Min. Ass. Can., Joint Annual Meeting, 50 pp. Marmo, J.S., 1992. The lower Proterozoic Hokkalampi paleosol in North Karelia, eastern Finland. In: Schidlowski, M., Golubic, S., Kimberley, M.M., Trudinger, P.A. (Eds.), Early Organic Evolution: Implications for Mineral and Energy Resources. Springer, Berlin Heidelberg, pp. 41–66. Martin, D.M., Clendenin, C.W., Krapez, B., McNaughton, N.J., 1998. Tectonic and geochronological constraints on late Archaean and Palaeoproterozoic stratigraphic correlation within and between the Kaapvaal and Pilbara Cratons. J. Geol. Soc. Lond. 155, 311–322. Master, S., Verhagen, B.T., Bassot, J.P., Beukes, N.J., Lemoine, S., 1993. Stable isotopic signatures of Paleoproterozoic carbonate rocks from Guinea, Sinegal, South Africa and Zimbabwe: constrains on the timing of the ca. 2 Ga “Lomagundi” ␦13 C excursion. Symposium: Early Proterozoic geochemical and structural constraints—metallogeny, vol. 23. Publication Occasionnelle, Dacar, Sénégal, pp. 38–41. Maynard, J.B., 2004. Manganiferous sediments, rocks, and ores. In: McKenzie, F.T. (Ed.), Treatise of Geochemistry. ElsevierPergamon, Oxford, pp. 289–308. McKirdy, D.M., Burgess, J.M., Lemon, N.M., Yu, X., Cooper, A.M., Gostin, V.A., Jenkins, J.F., Both, R.A., 2001. A chemostratigraphic overview of the late Cryogenian interglacial sequence in the Adelaide Fold-Thrust Belt, South Australia. Precambrian Res. 106, 149–186. Meyn, H.D., 1970. Geology of Hutton and Parkin Townships, Ontario Department of Mines, Geological Report 80. Toronto, 78 pp. Meyn, H.D., Matthews, M.K., 1991. A study of some uranium occurrences at the base of the Huronian Supergroup in an area north of Lake Wanapitei, District of Sudbury, Ontario. Ontario Geological Survey, Open File Report 5793, 58 pp. Meyn, H.D., Matthews, M.K., 1980. No. S26 Uranium deposits of the Cobalt Embayment. In: Milne, V.G., White, O.L., Barlow, R.B., Robertson, J.A., Colvine, A.C. (Eds.), Summary of Fieldwork, 1980 by the Ontario Geological Survey. Ontario Geological Survey, Ministry of Natural Resources Miscellaneous Paper 96, pp. 195–199. Milliken, K.L., 1979. The silicified evaporite syndrome-two aspects of silicification history of former evaporite nodules from Southern Kentucky and Northern Tennessee. J. Sed. Petr. 49, 245–256. Mossman, D.J., Harron, G.A., 1983. Origin and distribution of gold in the Huronian Supergroup, Canada—the case for Witwatersrandtype paleoplacers. Precambrian Res. 20, 543–583. Murakami, T., Utsunomiya, S., Imazu, Y., Prasad, N., 2001. Direct evidence of late Archean to early Proterozoic anoxic atmosphere from a product of 2.5 Ga old weathering. Earth Planetary Sci. Lett. 184, 523–528. Nance, R.D., Worsley, T.R., Moody, J.B., 1988. The Supercontinent cycle. Scientific Am. 256 (7), 72–79. Nesbitt, H.W., Young, G.M., 1982. Early Proterozoic climates and plate motions inferred from major element chemistry of lutites. Nature 299, 715–717. Nogueira, A.C.R., Riccomini, C., Sial, A.N., Moura, C.A.V., Fairchild, T.R., 2003. Soft-sediment deformation at the base of the Neoproterozoic Puga cap carbonate (southwestern Amazon craton, Brazil): confirmation of rapid icehouse to greenhouse transition in snowball Earth. Geology 31, 613–616. Ojakangas, R.W., 1997. Correlative sequence within the Marquette Range Supergroup (Michigan) and the Huronian Supergroup (Ontario); Glaciogenics, paleosols, and orthoquartzites (Abs), 43rd Inst. Lake Superior Geol., pp. 47–48. Ono, S., Eigenbrode, J.L., Pavlov, A.A., Kharecha, P., Rumble III, D., Kasting, J.F., Freeman, K.H., 2003. New insights into Archean sulfur cycle from mass-independent sulfur isotope records from the Hamersley Basin, Australia. Earth Planetary Sci. Lett. 213, 15–30. Panahi, A., Young, G.M., Rainbird, R.H., 2000. Behaviour of major and trace elements (including REE) during Paleoproterozoic pedogenesis and diagenetic alteration of an Archean granite near Ville Marie, Quebec, Canada. Geochim. Cosmochim. Acta 64, 2199–2220. Pavlov, A.A., Hurtgen, M.T., Kasting, J.F., Arthur, M.A., 2003. Metahne-rich Proterozoic atmosphere? Geology 31, 87–90. Pavlov, A.A., Kasting, J.F., Brown, L.L., 2000. Greenhouse warming by CH4 in the atmosphere of early Earth. J. Geophys. Res. 105 (E5), 11981–11990. Pelechaty, S.M., Kaufman, A.J., Grotzinger, J.P., 1996. Evaluation of ␦13 C chemostratigraphy for intrabasinal correlation: Vendian strata of northeast Siberia. GSA Bull. 108, 992–1003. Pickard, A.L., 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, 275– 315. Prasad, N., Roscoe, S.M., 1996. Evidence of anoxic to oxic atmospheric change during 2.45–2.22 Ga from lower and upper subHuronian paleosols, Canada. Catena 27, 105–121. Premo, W.R., Van Schmus, W.R., 1989. Zircon geocronology of Precambrian rocks in southeastern Wyoming and nothern Colorado. In: Grambling, J.A., Tewksbury, B.J. (Eds.), Proterozoic Geology of the Southern Rocky Mountains. Geol. Soc. Am. Spec. Pap. 235, pp. 1–12. Rainbird, R.H., Donaldson, J.A., 1988. Nonglaciogenic deltaic deposits in the early Proterozoic Gowganda Formation, Cobalt Basin, Ontario. Can. J. Earth Sci. 25, 710–724. Rainbird, R.H., Nesbitt, H.W., Donaldson, J.A., 1990. Formation and diagenesis of a sub-Huronian saprolith: comparison with a modern weathering profile. J. Geol. 98, 801–822. Robertson, J.A., 1968. Geology of Township 149 and Township 150, Ontario Dep. Mines GR 57. Toronto, 162 pp. Robertson, J.A., 1977. Geology of the Cutler area, District of Algoma. Ontario Div. Mines GR 147, 73 pp. A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Robertson, J.A., 1976. Geology of the Massey Area, Districts of Algoma, Manitoulin and Sudbury. Ontario, Ontario Div. Mines GR 136, 129 pp. Roscoe, S.M., 1996. Paleoplacer uranium, gold. In: Eckstrand, O.R., Sinclair, W.D., Thorpe, R.I. (Eds.), Geology of Canadian Mineral Deposit Types. Geological Society of America, The Geology of North America P-1, pp. 10–23. Roscoe, S.M., 1981. Temporal and other factors affecting deposition of uraniferous conglomerates. USGS Prof. Paper 1161-W, pp. W1-W17. Roscoe, S.M., 1969. Huronian rocks and uraniferous conglomerates in the Canadian Shield. Geol. Surv. Canada Paper 68-40, 205 pp. Roscoe, S.M., Minter, W.E.L., 1993. Pyritic paleoplacer gold and uranium deposits. In: Kirkham, R.V., Sinclair, W.D., Thorpe, R.I., Duke, J.M. (Eds.), Mineral Deposit Modeling. Geol. Ass. Canada Special Paper vol. 40, pp. 103–124. Rousell, D.H., Long, D.G.F., 1998. Are outliers of the Huronian Supergroup preserved in structures associated with the collapse of the Sudbury impact crater? J. Geol. 106, 407–419. Rye, R., Holland, H.D., 1998. Paleosols and the evolution of atmospheric oxygen: a critical review. Am. J. Sci. 298, 621–672. Rye, R., Kuo, P.H., Holland, H.D., 1995. Atmospheric carbon dioxide concentrations before 2.2 billion years ago. Nature 378, 603–605. Schidlowski, M., Eichmann, R., Junge, C.E., 1975. Precambrian Sedimentary Carbonates: carbon and oxygen isotope geochemistry and implications for the terrestrial oxygen budget. Precambrian Res. 2, 1–69. Schmidt, P.W., Williams, G.E., 1999. Paleomagnetism of the Paleoproterozoic hematitic breccia and paleosol at Ville-Marie, Québec: further evidence for the low paleolatitude of Huronian glaciation. Earth Planetary Sci. Lett. 172, 273–285. Schrag, D.P., Berner, R.A., Hoffman, P.F., Halverson, G.P., 2002. On the initiation of a snowball Earth. Geochem. Geophys. Geosyst., 3(6): 10.1029/2001GC000219 (available at http://www.g-cubed.org). Simonson, B.M., Schubel, K.A., Hassler, S.W., 1993. Carbonate sedimentology of the early Precambrian Hamersley Group of Western Australia. Precambrian Res. 60, 287–335. Smith, A.G., Pickering, K.T., 2003. Oceanic gateways as a critical factor to initiate icehouse Earth. J. Geol. Soc. Lond. 160, 337–340. Smith, M.D., Heaman, L.M., 1999. Constraints on the timing of felsic magmatism associated with the Matachewan igneous events: preliminary results for the Creighton granite, Ontario, Geol. Ass. Canada, Min. Ass. Canada Joint Ann. Meeting Abstract 24, p. 119. Snyder, G.L., Siems, D.F., Grossman, J.N., Ludwig, K.R., Nealey, L.D., 1995. Geology, petrochemistry, and geochronology of the Precambrian rocks of the Fletcher Park-Johnson Mountain area, Albany and Platte counties, Wyoming with a section on the Fletcher Park shear zone by Frost, B.R., Chamberlain, K.R., and Snyder, G.L., US Geol. Surv. Misc. Inv. Ser. Ma I-2233, scale 1:24,000, 3 sheets. Sohl, L.E., Christie-Blick, N., Kent, D.V., 1999. Paleomagnetic polarity reversals in Marinoan (ca, 600 Ma) glacial deposits of Australia: implications for the duration of low-latitude glaciation in Neoproterozoic time. Geol. Soc. Am. Bull. 111, 1120–1139. 205 Sumner, D.Y., Grotzinger, J.P., 2000. Late Archean aragonite precipitation: petrography, facies associations, and environmental significance. In: Grotzinger, J.P., James, N.P. (Eds.), Carbonate Sedimentation and Diagenesis in the Evolving Precambrian World, 67. SEPM Spec. Publ., pp. 123–144. Swart, Q.D., 1999. Carbonate rocks of the Paleoproterozoic Pretoria and Postmasburg Groups, Transvaal Supergroup. M.S. thesis, Rand Africaans University, Johannesburg, South Africa, 126 pp. Sylvester, G.H., 1973. Depositional environments of diamictites in the Headquarters Formation, Medicine Bow Mountains, southeastern Wyoming. M.S. thesis, University of Wyoming, Laramie, 84 pp. Tsikos, H., Moore, J.M., 1998. The Kalahari manganese field: an enigmatic association of iron and manganese. S. Afr. J. Geol. 101, 287–290. Valley, J.W., 1986. Stable isotope geochemistry of metamorphic rocks. Rev. Mineral 16, 445–489. Veizer, J., Clayton, R.N., Hinton, R.W., 1992. Geochemistry of Precambrian Carbonates: IV. Early Paleoproterozoic (2.25 ± 0.25 Ga) seawater. Geochim. Cosmochim. Acta 56, 875–885. Veizer, J., 1983. Chemical diagenesis of carbonates: Theory and application of trace element technique, Stable isotopes in Sedimentary Geology. SEPM Short Course 10, Dallas, 3-1-3-100. Williams, G.E., Schmidt, P.W., 1997. Paleomagnetism of the Paleoproterozoic Gowganda and Lorrain formations, Ontario; low paleolatitude for Huronian glaciation. Earth Planetary Sci. Lett. 153, 157–169. Williams, H., Hoffman, P.F., Lewry, J.F., Monger, J.W.H., Rivers, T., 1991. Anatomy of North America: thematic geologic portrayals of the continent. Tectonophysics 187, 117–134. Wing, B.A., Bekker, A., Brabson, L., Farquhar, J., Kaufman, A.J., Rumble, D., 2004. Atmospheric chemistry of SO2 during the earliest Huronian glaciation. GSA Abstracts with Programs 36(6), A-339. Wing, B.A., Brabson, E., Farquhar, J., Kaufman, A.J., Rumble III, D., Bekker, A., 2002. 33 S, ␦34 S and ␦13 C constraints on the Paleoproterozoic atmosphere during the earliest Huronian glaciation. Geochim. Cosmochim. Acta 66 (A840) (Abstr.). Wood, J., 1979. No. 18 Regional geology of the Cobalt Embayment, District of Sudbury, Nipissing, and Timiskaming: Summary of fieldwork. Ontario Geological Survey, Misc. paper 90, pp. 9–81. Xiao, S., Bao, H., Wang, H., Kaufman, A.J., Zhou, C., Li, G., Yuan, X., Ling, H., 2004. The Neoproterozoic Quruqtagh Group in eastern Chinese Tianshan: evidence for a post-Marinoan glaciation. Precambrian Res. 130, 1–26. Yang, W., Holland, H.D., 2003. The Hekpoort paleosol profile in Strata 1 at Gaborone, Botswana: soil formation during the Great Oxidation Event. Am. J. Sci. 303, 187–220. Young, G.M., 2002. Stratigraphic and tectonic settings of Proterozoic glaciogenic rocks and banded iron-formations: relevance to the snowball Earth debate. J. Afr. Earth Sci. 35, 451–466. Young, G.M., 1991. Stratigraphy, sedimentology and tectonic setting of the Huronian Supergroup, Geological Association of Canada, Mineralogical Association of Canada, Society of Economic Geologists, Joint Annual Meeting, Toronto ’91, Field trip B5: Guidebook, 34 pp. 206 A. Bekker et al. / Precambrian Research 137 (2005) 167–206 Young, G.M., 1988. Proterozoic plate tectonics, glaciation and ironformations. Sediment. Geol. 58, 127–144. Young, G.M., 1983. Tectono-sedimentary history of Early Proterozoic rocks of the northern Great Lakes region. In: Medaris Jr., L.G. (Ed.), Early Proterozoic Geology of the Great Lake Region. Geol. Soc. Am. Mem. vol. 160, pp. 15–32. Young, G.M., 1973. Tillites and aluminous quartzites as possible time markers for middle Precambrian (Aphebian) rocks of North America. In: Young, G.M. (Ed.), Huronian Stratigraphy and Sedimentation. Geol. Assn. Can. Spec. Pap. vol. 12, pp. 97– 127. Young, G.M., 1973a. Origin of carbonate-rich Early Proterozoic Espanola Formation, Ontario, Canada. Geol. Soc. Am. Bull. 84, 135–160. Young, G.M., 1972. Downward intrusive breccias in the Huronian Espanola Formation, Ontario, Canada. Can. J. Earth Sci. 9, 756–762. Young, G.M., Church, W.R., 1966. The Huronian System in the Sudbury District and adjoining areas of Ontario A review. Geol. Ass. Canada, Proc. Geol. Ass. Canada 17, 65–82. Young, G.M., Long, D.G.F., Fedo, C.M., Nesbitt, H.W., 2001. Paleoproterozoic Huronian basin: product of a Wilson cycle punctuated by glaciations and a meteorite impact. Sed. Geol. 141–142, 233–254. Zolnai, A.I., Price, R.A., Helmstaedt, H., 1984. Regional cross section of the Southern Province adjacent to Lake Huron, Ontario: implications for the tectonic significance of the Murray Fault Zone. Can. J. Earth Sci. 21, 447–456.
© Copyright 2026 Paperzz