Nature
Snowball Earth termination by destabilization of equatorial permafrost methane clathrate
DOI: 10.1038/nature06961
ISSN: 0028-0836
Accession: 00006056-200805290-00053
Author(s):

Kennedy, Martin1; Mrofka, David1; Borch, Chris von der2

Issue:
Volume 453(7195), 29 May 2008, pp 642-645
Publication Type:
[Letters]
Publisher:
© 2008 Nature Publishing Group
Institution(s):
1Department of Earth Science, University of California, Riverside, California 92521, USA. 2School of Chemistry, Physics and Earth Sciences, Flinders University, GPO Box 2100 Adelaide, South Australia, 5001 Australia.
doi:10.1038/nature06961
Received 24 September 2007; accepted 18 March 2008.
Supplementary Information is linked to the online version of the paper at www.nature.com/nature .
Author Information Reprints and permissions information is available at www.nature.com/reprints . Correspondence and requests for materials should be addressed to M.K. (martink@mail.ucr.edu).
Abstract

The start of the Ediacaran period is defined by one of the most severe climate change events recorded in Earth history-the recovery from the Marinoan 'snowball' ice age,~ 635 Myr ago (ref. 1). Marinoan glacial-marine deposits occur at equatorial palaeolatitudes2, and are sharply overlain by a thin interval of carbonate that preserves marine carbon and sulphur isotopic excursions of about -5 and +15 parts per thousand, respectively3-5; these deposits are thought to record widespread oceanic carbonate precipitation during postglacial sea level rise1,3,4. This abrupt transition records a climate system in profound disequilibrium3,6 and contrasts sharply with the cyclical stratigraphic signal imparted by the balanced feedbacks modulating Phanerozoic deglaciation. Hypotheses accounting for the abruptness of deglaciation include ice albedo feedback3, deep-ocean out-gassing during post-glacial oceanic overturn7 or methane hydrate destabilization8-10. Here we report the broadest range of oxygen isotope values yet measured in marine sediments (-25% to +12%) in methane seeps in Marinoan deglacial sediments underlying the cap carbonate. This range of values is likely to be the result of mixing between ice-sheet-derived meteoric waters and clathrate-derived fluids during the flushing and destabilization of a clathrate field by glacial meltwater. The equatorial palaeolatitude implies a highly volatile shelf permafrost pool that is an order of magnitude larger than that of the present day. A pool of this size could have provided a massive biogeochemical feedback capable of triggering deglaciation and accounting for the global postglacial marine carbon and sulphur isotopic excursions, abrupt unidirectional warming, cap carbonate deposition, and a marine oxygen crisis. Our findings suggest that methane released from low-latitude permafrost clathrates therefore acted as a trigger and/or strong positive feedback for deglaciation and warming. Methane hydrate destabilization is increasingly suspected as an important positive feedback to climate change11-13 that coincides with critical boundaries in the geological record14,15 and may represent one particularly important mechanism active during conditions of strong climate forcing.



Authigenic carbonates within the Reynella Member of the glacial Elatina Formation occur in coastal cliff exposures of the Marinoan type section in South Australia (Fig. 1 a-d) and provide constraints on the magnitude and timing of methane influence during deglaciation. The Reynella Member is part of a transgressive littoral succession that deepens into the shelfal deposits of the Seacliff sandstone and interbedded cap carbonate of the Nuccaleena Formation 16 (Fig. 1f). The Nuccaleena Formation is the Marinoan cap carbonate that houses the global stratotype section and point for the Ediacaran period, and is widely believed to provide a global tie point shared by the integrated sea level rise driven by global deglaciation 4. On the basis of the low palaeolatitude 2, minimal glacial evidence and trans-gressive sediments, we interpret deposition of the Reynella siltstone to coincide with local ice withdrawal preceding global eustatic rise recorded by the Nuccaleena Formation (Fig. 1f).



Graphic
[Help with image viewing]
[Email Jumpstart To Image]
Figure 1 Photographic examples and stratigraphic context for Reynella seep carbonates. Distinctive sedimentary facies include vertical cemented chimneys rising through coarse sandstone with sparry cavity cements (a); shallow channels (backpack for scale) with basal lags of authigenic carbonate clasts eroded from exhumed seeps (b); and breccia beds comprising cemented sandstone (5-cm lens cap for scale) resulting from fluidization (c, d). A thin section of tidal rhythmites (e) shows well-developed neap-spring bundles of couplets from samples collected between authigenic carbonate hard grounds. Schematic section (f) shows the association of tidal, seep and stratigraphic units. Section and photos located between Hallet Cove and Marino Rocks around 36° 03' 45" S, 138° 30' 05" E.

These carbonates form an irregular plumbing system containing sub-vertical chimneys 1-2 m in diameter and branching pipes (Fig. 1a) that coalesce into hard grounds. Microcrystalline dolomite and calcite cement fills pores in sandstone and forms up to 60% of the sediment, indicating precipitation before compaction. In plan view, spar-filled tubes are concentrated within the cores of 1-m-diameter chimneys ringed by cemented sandstone. Morphologies include microcrystalline pore-filling calcite and Fe-rich dolomite with pyrite crystals >10 µm indicative of sulphate reduction, acicular cements (formerly aragonite) fringing synsedimentary breccia clasts of carbonate-cemented sandstone, spar-lined vugs and tubes with Mn-cement coatings, sparry calcite-filled veins often forming a honeycomb crust morphology, and isolated centimetre-scale bladed crystal fans of calcite pseudomorphs after barite. Calcite spar is present in early veins, thin vertical tubes and vugs, and within hard grounds. Both dolomite and calcite cements and calcite spar comprise clasts eroded and reworked by currents within channelized, cross-bedded sand and breccia deposits (Fig. 1b). These indicate a dynamic system, including erosion with re-cementation in shallow seafloor sediments rather than a more typical history of cementation during burial. Breccias of similar composition occur as metre-scale, massive, vertically oriented cross-cutting fluidized beds (Fig. 1c, d) that erupted on the sea floor and were reworked by seafloor currents into cross-bedded channel-fill. Sedimentary features indicative of vertical fluid or gas migration, overpressure, seafloor reworking and complex seafloor paragenesis are diagnostic of a methane cold-seep origin 17,18. This complex morphology contrasts with the regionally extensive laminated dolomicrite of the overlying Nuccaleena Formation and the sparry calcite that fills cross-cutting fractures aligned with regional folding.

Authigenic cements within these morphologies show the broadest range of [delta]18O values ever reported from marine sediments 19, ranging from +12‰ to -25‰ (PDB) at centimetre-scale resolution (Fig. 2; Supplementary Table 1). Potential controls on oxygen isotope values include precipitating fluid temperature, fractionation associated with crystal formation, and the composition of the source pore fluid. The heaviest [delta]18O values (Fig. 2; +12‰) occur in pore-filling dolomite cement and are too enriched in 18O to be attributed to conventional mechanisms such as direct dolomite precipitation from sea water at 0° C (ref. 17). Values over 10‰ in marine sediments are known only from methane cold-seep environments where 18O is enriched in pore fluids during methane clathrate formation and dolomite precipitation results from sulphate oxidation of methane 17,20. Values of [delta]13C that range from -10‰ to +10‰; PDB (Fig. 2; Supplementary Table 1) within Reynella cements are also common within cold-seep-associated cements and attributed to concurrent methanotrophy and mixing between isotopically heavy methanogenic CO2 and marine bicarbonate 17,18,20.



Graphic
[Help with image viewing]
[Email Jumpstart To Image]
Figure 2 Values of [delta]13C and [delta]18O for Nuccaleena Formation and Reynella Member cements and veins. Open triangles, calcite-cemented siltstone; filled triangles, dolomite-cemented siltstones; open squares, spar-filled cavities; filled squares, early calcite spar; asterisks, early spar-filled tubes; dashes, late cross-cutting veins; crosses, Nuccaleena cap carbonate. Upper trend comprises filled symbols and asterisks, lower trendline includes open symbols. Later phases (crosses and dashes) not included in trend lines. Horizontal bars represent ranges of [delta]18O values from (A) isotopically light meltwater-related carbonates 23,24; (B) clathrate enrichment of 18O in seep carbonates 18,19,21; (C) predicted [delta]18O values from carbonate precipitated in 0 °C meteoric water from >70° latitude (-26‰ PDB calcite from -30‰ SMOW water) to low latitude (+14‰ PDB from -2‰ SMOW water)22; and (D) reported range of [delta]18O values for all marine carbonate rocks 30.

Values of [delta]18O for calcite cements within the Reynella Member define a mixing line with a highly depleted meteoric endmember (-25‰ PDB; R2 = 0.84). Although low [delta]18O values can occur during metamorphism by equilibration at >180 °C, a metamorphic origin is incompatible with the presence of temperature-sensitive detrital, dioctahedral-smectite in the laminated mud fraction of the host rock which would have converted to illite above 100 °C (Supplementary Fig. 1; Supplementary Table 2). Thermal equilibration is inconsistent with (1) the microcrystalline morphology of the cement, (2) the isotopic heterogeneity on millimetre scales, and (3) the relatively heavy [delta]18O values of the Nuccaleena cap carbonate and late-stage (tectonic) cross-cutting spar fill fractures (Fig. 2). The isotopic signature of the most depleted cements (-25‰) requires porewater values as low as about -28 to -30‰ (see Supplementary Information), which, in the modern hydrological cycle, is unique to melt-water released from glacial ice 21. In the geological record, these low values are found only in thermally altered carbonate spar 19 or formed from glacial melt water in ancient carbonates 22,23..

The mixing relation connecting enriched and depleted [delta]18O end-members (Fig. 2), sub-centimetre spatial association and synsedimentary timing relations between stable isotopic endmembers indicates formation of seeps in a shallow environment in which meteoric fluid derived from an ice sheet mixed with 18O-enriched clathrate-derived fluids. Such an environment occurs in shelf settings where the hydraulic drive of a melting ice sheet charges the water table, flushing and destabilizing clathrate fields.

Sedimentary evidence indicates that the Reynella seeps formed syndepositionally in high-energy, intertidal, sandstones. The authigenic carbonate occurs within a succession containing superbly preserved tidal bundles (Fig. 1e) including well-developed spring and neap cycles 24, broad (2 m), shallow (1 m) channels (Fig. 1b), reactivation surfaces, bimodal palaeocurrents within cross-bedded sands, and herringbone and flaser lamination, all indicative of a high-energy, macrotidal environment 16,24 unique to coastal or estuarine settings. Reworking of authigenic carbonate in channels followed by subsequent cementation (Fig. 1b) indicates that seep hard grounds formed throughout the tidally influenced deposition of sand.

Methane clathrates form when methane produced by thermogenic or biogenic breakdown of organic matter at depth migrates upwards to stabilize more than 200 m below the sea floor or land surface in a narrow range of pressure and temperature 12,25. To be stable, marine clathrates require both sea water at 0 °C and the additional pressure derived from more than 250 m of water 12. Because the Reynella Member consists of intertidal deposits within a deepening succession midway between the preceding glacial sealevel lowstand and the shelfal deposits of the Seacliff sandstone/Nuccaleena dolomite, the water could not have been deep enough to stabilize clathrates. The stabilization thus must have resulted either from the pressure of a shelf-borne ice sheet or from the permafrost conditions (<0 °C) indicated by local evidence of frost wedging 26. Permafrost clathrates that accumulated beneath the broad exposed shelves of the Arctic during the Pleistocene are now destabilizing as a result of post-glacial flooding and global warming 11,12,25. Similarly, retreat of the Marinoan ice sheet would have destabilized methane clathrates by (1) reducing hydrostatic pressure, (2) warming permafrost shelves with flooding sea water, and (3) flushing clathrate fields with meteoric water derived from meltwater discharge 9,22,23.

The extension of a volatile permafrost clathrate pool to equatorial palaeolatitudes has important implications for the role of methane in the climate system. Initial warming or, more simply, inherent ice-sheet instability 27, might have been enough to trigger local methane destabilization through meteoric flushing of the permafrost clathrate field on the Marinoan shelf, particularly in this low-latitude setting 2 where it was most sensitive to subtle external forcing. Once initiated, clathrate destabilization provides a positive feedback through greenhouse warming caused by methane and CO2, and is amplified by further destabilization in progressively higher latitudes 9,11 of permafrost and eventually marine clathrates. Extension of the Holocene permafrost reservoir, occurring over the Arctic shelves and estimated 25 at 400 gigatonnes, to total present-day shelf area (including low latitudes) could increase the potential size to ~3,000 Gt (see Supplementary Information). This estimate is conservative because it does not include the extensive intracratonic basinal areas characteristic of the Neoproterozoic (for example 9), nor does it consider destabilization of the more conventional and tenfold larger marine clathrate reservoir 25. It is, however, comparable to past gas hydrate events with estimated releases of ~5,000 Gt in the Toarcian 14 or ~1,000-2,100 Gt in the Palaeocene 15 estimated from marine excursions of ~7‰ and -3‰ [delta]13C, respectively. Permafrost clathrate destabilization could readily provide a trigger for feedbacks to deglaciation such as albedo, marine clathrate destabilization or carbonate reallocation 3,28.

Alkalinity increase capable of accounting for cap carbonate deposition could arise on shorter timescales (<104 yr) from sulphate oxidation of methane following exhaustion of dissolved oxygen in pore and sea water 10 and on longer timescales (>105 yr) from oxidation of methane to CO2 in the atmosphere driving increased continental chemical weathering. The duration of methane release would most likely be controlled by pulsed marine clathrate destabilization similar to that in the early Toarcian 14, which occurred over ~200 kyr. Like other later methane destabilization events 14, release of 3,000 Gt of methane (~2 × 1017 moles of carbon, excluding marine clath rates) into sea water over about 103-104 yr could substantially influence the redox state of parts of the ocean, significantly depleting the dissolved oxygen pool and consuming marine sulphate which maintained a probable stock of less than 4.5 × 1016 mole O2 (ref. 29) and 5 × 1017 mole SO4 (refs 2-5) (see Supplementary Information). Significant sulphate oxidation of methane within pore and sea water predicts the observed 15‰ positive shift in marine sulphate [delta]34S values 5. The release of 3,000 Gt of -60‰ [delta]13C methane from the permafrost reservoir could alone account for a shift of more than -3‰ (see Supplementary Information) of the enigmatic negative [delta]13C excursion in marine carbonates that defines the base of the Ediacaran period 4.

The distinctive features of Marinoan deglaciation that define the base of the Ediacaran period 4 can be attributed to the effects of permafrost methane clathrate destabilization. In contrast to the balanced feedbacks and progressive glacial-interglacial cycles of Cenozoic deglaciation, the violent opening of the highly volatile shelf-permafrost methane clathrate pool could act as a trigger to catastrophic climate and biogeochemical reorganization of the Earth system, abruptly bringing the long-lived and icy Cryogenian period to a close and setting the stage for the appearance of metazoans and dominance of a new Earth system. This event both identifies the range of function of the climate system, and demonstrates a mechanism activated by strong climate forcing not unlike projected future effects of atmospheric CO2.

METHODS SUMMARY

Detailed methods are given in the Supplementary Information. We characterized the petrography and mineralogy by X-ray diffraction, petrographic microscopy and back scattered electron microscopy. Isotopic results were obtained on a VG Isogas Prism II, with reproducibility of [delta]18O-enriched samples confirmed on a Thermo Delta-V Plus and at the University of California, Davis, and agree within 0.1‰. [delta]18O = [(18O/16O)sample/(18O/16O)PDB]-1; [delta]13C = [(13C/12C)sample/(13C/12C)PDB]-1.

Acknoledgements

We thank the NSF and NASA Exobiology for funding, A. Derkowski and T. Bristow for assistance with clay mineralogy, D. Winter for replicate isotopic analyses and I. Fairchild for a review.

Author Contributions

M.J.K. collected samples and wrote the manuscript; D.D.M. collected samples, performed geochemical analyses and wrote the supplemental section; C.v.d.B. collected samples. All authors discussed the results and commented on the manuscript.

1. Fairchild, I. J. & Kennedy, M. J. Neoproterozoic glaciation in the Earth system. J. Geol. Soc. Lond. 164, 895-921 (2007). [Context Link]

2. Sohl, L. E., Christie-Blick, N. & Kent, D. V. 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 (1999). Internet Resources Document Delivery Library Holdings [Context Link]

3. Hoffman, P. F., Kaufman, A. J., Halverson, G. P. & Schrag, D. P. A Neoproterozoic snowball Earth. Science 281, 1342-1346 (1998). [Context Link]

4. Knoll, A. H., Walter, M. R., Narbonne, G. M. & Christie-Blick, N. The Ediacaran Period: a new addition to the geologic time scale. Lethaia 39, 13-30 (2006). [Context Link]

5. Halverson, G. P. & Hurtgen, M. T. Ediacaran growth of the marine sulfate reservoir. Earth Planet. Sci. Lett. 263, 32-44 (2007). [Context Link]

6. Allen, P. A. & Hoffman, P. F. Extreme winds and waves in the aftermath of a Neoproterozoic glaciation. Nature 433, 123-127 (2005). [Context Link]

7. Grotzinger, J. P. & Knoll, A. H. Anomalous carbonate precipitates; is the Precambrian the key to the Permian? Palaios 10, 578-596 (1995). Internet Resources Bibliographic Links Document Delivery Library Holdings [Context Link]

8. Jiang, G., Kennedy, M. J. & Christie-Blick, N. Stable isotopic evidence for methane seeps in Neoproterozoic postglacial cap carbonates. Nature 426, 822-826 (2003). [Context Link]

9. Shields, G. A., Deynoux, M., Strauss, H., Paquet, H. & Nahon, D. Barite-bearing cap dolostones of the Taoudeni Basin, northwest Africa; sedimentary and isotopic evidence for methane seepage after a Neoproterozoic glaciation. Precambr. Res. 153, 209-235 (2007). [Context Link]

10. Kennedy, M. J., Christie-Blick, N. & Sohl, L. E. Are Proterozoic cap carbonates and isotopic excursions a record of gas hydrate destabilization following Earth's coldest intervals? Geology 29, 443-446 (2001). Internet Resources Document Delivery Library Holdings [Context Link]

11. Nisbet, E. G. The end of the ice age. Can. J. Earth Sci. 27, 148-157 (1990). [Context Link]

12. MacDonald, G. J. Role of methane clathrates in past and future climates. Clim. Change 16, 247-281 (1990). Internet Resources Bibliographic Links Document Delivery Library Holdings [Context Link]

13. Archer, D. Methane hydrate stability and anthropogenic climate change. Biogeosciences 4, 521-544 (2007). [Context Link]

14. Kemp, D. B., Coe, A. L., Cohen, A. S. & Schwark, L. Astronomical pacing of methane release in the Early Jurassic period. Nature 437, 396-399 (2005). Ovid Full Text Full Text Internet Resources Bibliographic Links Document Delivery Library Holdings [Context Link]

15. Dickens, G. R., O'Neil, J. R., Rea, D. K. & Owen, R. M. Dissociation of oceanic methane hydrate as a cause of the carbon isotope excursion at the end of the Paleocene. Paleoceanography 10, 965-971 (1995). Internet Resources Document Delivery Library Holdings [Context Link]

16. Preiss, W. V. The Adelaide Geosyncline, late Proterozoic stratigraphy, sedimentation, palaeontology, and tectonics. Bull. Geol. Surv. S. Aust. 53, 1-438 (1987). [Context Link]

17. Greinert, J., Bohrmann, G. & Suess, E. Gas hydrate-associated carbonates and methane-venting at Hydrate Ridge; classification, distribution, and origin of authigenic lithologies. Geophys. Monogr. 124, 99-113 (2000). [Context Link]

18. Aiello, I. W., Garrison, R. E., Moore, J. C., Kastner, M. & Stakes, D. S. Anatomy and origin of carbonate structures in a Miocene cold-seep field. Geology 29, 1111-1114 (2001). Internet Resources Document Delivery Library Holdings [Context Link]

19. Shields, G. & Veizer, J. Precambrian marine carbonate isotope database versions 1.1. Geochem. Geophys. Geosyst. 3, 1-12 (2002). [Context Link]

20. Ritger, S., Carson, B. & Suess, E. Methane-derived authigenic carbonates formed by subduction-induced pore-water expulsion along the Oregon/Washington margin; with Suppl. Data 87-02. Geol. Soc. Am. Bull. 98, 147-156 (1987). [Context Link]

21. Bowen, G. J. & Wilkinson, B. Spatial distribution of [delta]18O in meteoric precipitation. Geology 30, 315-318 (2002). Internet Resources Document Delivery Library Holdings [Context Link]

22. Herbert, C. T. & Compton, J. S. Depositional environments of the Lower Permian Dwyka Diamictite and Prince Albert Shale inferred from the geochemistry of early diagenetic concretions, southwest Karoo Basin, South Africa. Sedim. Geol. 194, 263-277 (2007). [Context Link]

23. Himmler, T., Freiwald, A., Stollhofen, H. & Peckmann, J. Late Carboniferous hydrocarbon-seep carbonates from the glaciomarine Dwyka Group, southern Namibia. Palaeogeogr. Palaeoclimatol. Palaeoecol. 257, 185-197 (2008). [Context Link]

24. Williams, G. E. in Clastic Tidal Sedimentology (eds Smith, D. G., Reinson, G. E., Zaitlin, B. A. & Rahmani, R. A.) 161-178 (Canadian Society of Petroleum Geologists, Calgary, 1991). [Context Link]

25. Kvenvolden, K. A. Methane hydrate; a major reservoir of carbon in the shallow geosphere? Chem. Geol. 71, 41-51 (1988). [Context Link]

26. Williams, G. E. Precambrian permafrost horizons as indicators of palaeoclimate. Precambr. Res. 32, 233-242 (1986). [Context Link]

27. MacAyeal, D. R. Binge/Purge oscillations of the Laurentide ice sheet as a cause of the North Atlantic Heinrich events. Paleoceanography 8, 775-784 (1993). Internet Resources Document Delivery Library Holdings [Context Link]

28. Ridgwell, A. J., Kennedy, M. J. & Caldeira, K. Carbonate deposition, climate stability, and Neoproterozoic ice ages. Science 302, 859-862 (2003). [Context Link]

53. Canfield, D. E., Poulton, S. W. & Narbonne, G. M. Late-Neoproterozoic deep-ocean oxygenation and the rise of animal life. Science 315, 92-95 (2007). [Context Link]

30. Hoefs, J. Stable Isotope Geochemistry 55 (Springer, Berlin, 2004). [Context Link]