Geology and geochemistry of paleosols developed on the Hekpoort Basalt, Pretoria Group, South Africa
The Hekpoort paleosols comprise a regional paleoweathering horizon developed on 2.224 +/- 0.021 Ga basaltic andesite lavas at the top of the Hekpoort Formation of the Pretoria Group, Transvaal Supergroup, South Africa. In five separate profiles, from outcrops along road cuts near Waterval Onder and the Daspoort Tunnel and in three drill cores from the Bank Break Area (BB3, BB8, and BB14), the top of the paleosol is a sericite-rich zone. The sericite zone grades downward into a chlorite-rich zone. In core BB8 and in the road cut at the Daspoort Tunnel, we sampled the underlying or parent basaltic andesite into which the chlorite zone grades. We did not obtain samples of the parent material at Waterval Onder and in cores BB3 and BB14, but chemical analyses indicate that the chlorite and sericite zones in these profiles derive from underlying lavas similar to the ones we sampled in core BB8 and at the Daspoort Tunnel. The presence of apparent rip-up clasts of the paleosol in the overlying ironstones of the Strubenkop Formation in the cores from Bank Break makes it very unlikely that most of the alteration was a result of interactions with hydrothermal fluids. Desiccation cracks at the top of the paleosol that were filled with sand during the deposition of the overlying sediments at Waterval Onder point to a subaerial weathering origin. Very little, if any, Al, Ti, Zr, V, or Cr moved a discernible distance during weathering of any of the five profiles. The vertical distribution of Fe, Mg, Mn, Ni, and Co indicates that these elements were largely removed from the top of the soil during weathering. The overall abundance of these elements in each of the profiles indicates that a significant fraction of the complement lost from the top subsequently reprecipitated in the lower portion of the soil as constituents of an Fe2(+) -rich smectite. The loss of Fe from the top of the soil during weathering of the Hekpoort paleosols indicates that atmospheric PO2 was less than 8 x 10(-4) atm about 2.22 Ga. Fe2(+) -rich smectite should only precipitate during soil formation if atmospheric PCO2 is less than or equal to 2 x 10(-2) atm (Rye, Kuo, and Holland, 1995). Ca and Na were largely lost during weathering. Some Na was apparently added to the sericite zone in cores BB3, BB8, and BB14 after weathering. All five profiles are enriched in K and Rb, and most are enriched in Ba. The distribution of these elements indicates that they all were added during post-weathering hydrothermal metasomatism. Rb-Sr analysis of the paleosol at the Daspoort Tunnel indicates that metasomatism last affected that profile 1.925 +/- 0.032 Ga (Macfarlane and Holland, 1991).
Paleosols and the evolution of atmospheric oxygen: a critical review
A number of investigators have used chemical profiles of paleosols to reconstruct the evolution of atmospheric oxygen levels during the course of Earth history (Holland, 1984, 1994; Kirkham and Roscoe, 1993; Ohmoto, 1996). Over the past decade Holland and his co-workers have examined reported paleosols from six localities that formed between 2.75 and 0.45 Ga. They have found that the chemical profiles of these paleosols are consistent with a dramatic change in atmospheric PO2 between 2.2 and 2.0 Ga from < or = 0.002 to > or = 0.03 atm (Holland, 1994). Ohmoto (1996) examined chemical data from twelve reported paleosols ranging in age from 2.9 to 1.8 Ga. He concluded that these chemical profiles indicate that atmospheric PO2 has not changed significantly during the past 3.0 Ga. We seek to resolve the conflict between these reconstructions through a broader examination of the paleosol literature, both to determine which reported paleosols can be definitively identified as such and to determine what these definite paleosols tell us about atmospheric evolution. We here review reports describing over 50 proposed paleosols, all but two are older than 1.7 Ga. Our review indicates that 15 of these reported paleosols can be definitively identified as ancient soils. The behavior of iron uring the formation of these 15 paleosols provides both qualitative and semiquantitative information about the evolution of the redox state of the atmosphere. Every definitely identified pre-2.44 Ga paleosol suffered significant Fe loss during weathering. This loss indicates that atmospheric PO2 was always less than about 5 x l0(-4) atm prior to 2.44 Ga. Analysis of the Hokkalampi paleosol (2.44-2.2 Ga) (Marmo, 1992) and the Ville Marie paleosol (2.38-2.215 Ga) (Rainbird, Nesbitt, and Donaldson, 1990) yield ambiguous results regarding atmospheric PO2. Loss of Fe during the weathering of the 2.245 to 2.203 Ga Hekpoort paleosol (Button, 1979) indicates that atmospheric PO2 was less than 8 x 10(-4) atm shortly before 2.2 Ga. The presence of red beds immediately overlying the Hokkalampi, Ville Marie, and Hekpoort paleosols suggests that by about 2.2 Ga there was an unquantified but substantial amount of oxygen in the atmosphere. Iron loss was negligible during formation of the 2.2 to 2.0 Ga Wolhaarkop (Holland and Beukes, 1990) and Drakenstein (Wiggering and Beukes, 1990) paleosols and during formation of all the later paleosols we previewed. Thus, atmospheric PO2 probably has been > or = 0.03 atm since sometime between 2.2 and 2.0 Ga.
The carbon-isotopic composition of Proterozoic carbonates: Riphean successions from northwestern Siberia (Anabar Massif, Turukhansk Uplift)
Thick carbonate-dominated successions in northwestern Siberia document secular variations in the C-isotopic composition of seawater through Mesoproterozoic and early Neoproterozoic (Early to early Late Riphean) time. Mesoproterozoic dolomites of the Billyakh Group, Anabar Massif, have delta 13C values that fall between 0 and -1.9 permil versus PDB, with values in the upper part of the succession (Yusmastakh Formation) consistently higher than those of the lower (Ust'-Il'ya and Kotuikan formations). Consistent with available biostratigraphic and radiometric data, delta 13C values for Billyakh carbonates compare closely with those characterizing early Mesoproterozoic carbonates (about 1600-1200 Ma) worldwide. In contrast, late Mesoproterozoic to early Neoproterozoic limestones and dolomites in the Turukhansk Uplift exhibit moderate levels of secular variation. Only the lowermost carbonates in the Turukhansk succession (Linok Formation) have delta 13C values that approximate Billyakh values. Higher in the Turukhansk succession, delta 13C values vary from -2.7 to +4.6 permil (with outliers as low as -5.0 permil interpreted as diagentically altered). Again, consistent with paleontological and radiometric data, these values compare well with isotopic values from 1200 to 850 Ma successions elsewhere. Five sections measured in different parts of the Turukhansk basin show nearly identical patterns of variation, confirming that carbonate delta 13C correlates primarily with time and not facies. The Siberian sections illustrate the potential of integrated biostratigraphic and chemostratigraphic data in the intra- and interbasinal correlation of Mesoproterozoic and early Neoproterozoic rocks.
Organic acids in hydrothermal solutions: standard molal thermodynamic properties of carboxylic acids and estimates of dissociation constants at high temperatures and pressures
Experimental standard partial molal volumes, heat capacities, and entropies as well as apparent standard partial molal enthalpies and Gibbs free energies of mono- and dicarboxylic acids and their anions at low temperatures and pressures are used to generate correlations for predicting the same properties at high temperatures and pressures for 59 carboxylic and 18 hydroxyacid species with the revised Helgeson-Kirkham-Flowers (HKF) equation of state. Predicted equilibrium dissociation constants are compared with experimental values from the literature and tabulated as functions of pressure and temperature for 25 carboxylic acids and nine hydroxyacids. Close agreement between independent predictions and experimental data supports the generality of the computational techniques and the accuracy of predicted data. These results allow incorporation of a wide variety of organic acids into quantitative interpretations of geochemical processes.
Modelling the Phanerozoic carbon cycle and climate: constraints from the 87Sr/86Sr isotopic ratio of seawater
A numerical model describing the coupled evolution of the biogeochemical cycles of carbon, sulfur, calcium, magnesium, phosphorus, and strontium has been developed to describe the long-term changes of atmospheric carbon dioxide and climate during the Phanerozoic. The emphasis is on the effects of coupling the cycles of carbon and strontium. Various interpretations of the observed Phanerozoic history of the seawater 87Sr/86Sr ratio are investigated with the model. More specifically, the abilities of continental weathering, volcanism, and surface lithology in generating that signal are tested and compared. It is suggested that the observed fluctuations are mostly due to a changing weatherability over time. It is shown that such a conclusion is very important for the modelling of the carbon cycle. Indeed, it implies that the conventional belief that the evolution of atmospheric carbon dioxide and climate on a long time scale is governed by the balance between the volcanic input of CO2 and the rate of silicate weathering is not true. Rather carbon exchanges between the mantle and the exogenic system are likely to have played a key role too. Further, the increase of the global weathering rates with increasing surface temperature and/or atmospheric CO2 pressure usually postulated in long-term carbon cycle and climate modelling is also inconsistent with the new model. Other factors appear to have modulated the weatherability of the continents through time, such as mountain building and the existence of glaciers and ice sheets. Based on these observations, a history of atmospheric carbon dioxide and climate during Phanerozoic time, consistent with the strontium isotopic data, is reconstructed with the model and is shown to be compatible with paleoclimatic indicators, such as the timing of glaciation and the estimates of Cretaceous paleotemperatures.
Sulfate reduction in deep-sea sediments
Sulfate reduction rates calculated from about 200 DSDP pore water sulfate profiles have been contoured and plotted on a map covering most areas of the world ocean. Rates show a remarkable spatial consistency, with high rates observed near the continental margins, becoming progressively lower toward the central ocean basins. Relatively elevated rates are also found in the eastern equatorial Pacific, a site of upwelling and correspondingly high rates of primary organic production. Overall, the distribution of sulfate reduction in pelagic sediments looks very similar to the distribution of primary organic carbon production. When rates are directly compared, however, the correlation between sulfate reduction and primary production is only moderately strong. Perhaps the most important influence on sulfate reduction is sediment deposition rate and the control this has over the fraction of the sedimentary organic carbon flux that becomes available for sulfate reduction. The slower the rate of sediment deposition the more time for oxic respiration and the less organic carbon that escapes to the zone of sulfate reduction. To predict most accurately sulfate reduction rates, however, the variables of primary production, water depth, and sediment deposition rate must all be integrated.
The Flin Flon paleosol and the composition of the atmosphere 1.8 BYBP
Within the 1800 to 1900 my old Flin Flon-Snow Lake greenstone belt, Amisk Group volcanics are overlain by Missi Group fluvial sediments. Several localities along the Missi-Amisk contact, the volcanics show evidence of subaerial weathering. Field relationships, mineralogical evidence, and chemical analyses confirm that this alteration zone is a paleosol. Pedogenic fabrics and mineralogy were somewhat obscured by greenschist-grade metamorphism associated with the Hudsonian orogeny (1750 my). This is especially true in the upper meter of the paleosol, where metamorphic paragonite and sericitic micas developed in a crenulated fabric. This metamorphism did not, however, obliterate the imprint of weathering on the Amisk volcanics. Features characteristic of well-drained modern soils are evident in the paleosol. Corestones of spheroidally weathered pillow lavas occur at depth within the paleosol (Cr horizon). The corestones decrease in size upward and eventually disappear into a hematite-rich horizon at the top of the paleosol. These macroscopic changes are accompanied by a decrease in CaO and MgO and by an increase in Al2O3, TiO2, and total iron toward the paleosol-Missi contact. Ferrous iron decreases upward toward the contact; FeO was apparently oxidized to ferric iron and retained within the paleosol during weathering. The oxidation and retention of iron within the Flin Flon paleosol indicates that PO2 was probably > or = 10(-2) P.A.L. at the time of weathering. The behavior of iron in the Flin Flon paleosol contrasts sharply with its behavior in the 2200 my Hekpoort paleosol, which is strongly depleted in iron. This difference suggests that a significant increase in the ratio of PO2/PCO2 in the atmosphere took place between 2200 and 1800 mybp.
A new model for atmospheric oxygen over Phanerozoic time
A mathematical model has been constructed that enables calculation of the level of atmospheric O2 over the past 570 my from rates of burial and weathering of organic carbon (C) and pyrite sulfur (S). Burial rates as a function of time are calculated from an assumed constant worldwide clastic sedimentation rate and the relative abundance, and C and S contents, of the three rock types: marine sandstones and shales, coal basin sediments, and other non-marine clastics (red beds, arkoses). By our model, values of O2 versus time, using a constant total sedimentation rate, agree with those for variable sedimentation derived from present-day rock abundances and estimates of erosional losses since deposition. This agreement is the result of our reliance on the idea that any increase in total worldwide sediment burial, with consequently faster burial of C and S and greater O2 production, must be accompanied by a corresponding increase in erosion and increased exposure of C and S on the continents to O2 consumption via weathering. It is the redistribution of sediment between the three different rock types, and not total sedimentation rate, that is important in O2 control. To add stability to the system, negative feedback against excessive O2 fluctuation was provided in the modeling by the geologically reasonable assignment of higher weathering rates to younger rocks, resulting in rapid recycling of C and S. We did not use direct O2 negative feedback on either weathering of C and S or burial of C because weathering rates are assumed to be limited by uplift and erosion, and the burial rate of C limited by the rate of sediment deposition. The latter assumption is the result of modern sediment studies which show that marine organic matter burial occurs mainly in oxygenated shallow water and is limited by the rate of supply of nutrients to the oceans by rivers. Results of the modeling indicate that atmospheric O2 probably has varied appreciably over Phanerozoic time. During the Late Carboniferous and Permian periods O2 was higher than previously because of the rise of vascular land plants and the widespread burial of organic matter in vast coal swamps. A large decrease in O2 during the Late Permian was due probably to the drying-up of the coal swamps and deposition of a large proportion of total sediment in C and S-free continental red beds. Sensitivity study shows that major parameters affecting results are the mean C concentration in coal basins and the relative sizes of the reservoirs of young (rapidly recycled) versus old rocks. Less sensitivity was found for changes over time in total land area undergoing weathering and the use of direct O2 negative feedback on marine carbon burial. Good agreement for rates of C burial calculated via our model and via independent models, which are based on the use of stable carbon isotopes, indicates that the dominant factor that has brought about changes in atmospheric O2 level (and the isotopic composition of dissolved inorganic carbon in seawater) over Phanerozoic time is sedimentation and not weathering or higher temperature phenomena such as basalt-seawater reaction.
The post-Paleozoic chronology and mechanism of 13C depletion in primary marine organic matter
Carbon-isotopic compositions of geoporphyrins have been measured from marine sediments of Mesozoic and Cenozoic age in order to elucidate the timing and extent of depletion of 13C in marine primary producers. These results indicate that the difference in isotopic composition of coeval marine carbonates and marine primary photosynthate was approximately 5 to 7 permil greater during the Mesozoic and early Cenozoic than at present. In contrast to the isotopic record of marine primary producers, isotopic compositions of terrestrial organic materials have remained approximately constant for this same interval of time. This difference in the isotopic records of marine and terrestrial organic matter is considered in terms of the mechanisms controlling the isotopic fractionation associated with photosynthetic fixation of carbon. We show that the decreased isotopic fractionation between marine carbonates and organic matter from the Early to mid-Cenozoic may record variations in the abundance of atmospheric CO2.
A hybrid model of the CO2 geochemical cycle and its application to large impact events
A hybrid model of the carbonate-silicate geochemical cycle is presented which is capable of calculating the partitioning of carbon dioxide between the atmosphere, ocean, and sedimentary rocks. The ocean is subdivided into a shallow, mixed layer, which remains in equilibrium with the atmosphere, and a massive, deep layer which does not. Gradients in dissolved carbon content are established between the mixed layer and the deep ocean as a consequence of downward fluxes of fecal matter and of dead planktonic organisms. The dissolved carbon content and alkalinity of the ocean as a whole are controlled by weathering and metamorphism of sedimentary rocks. Equilibrium solutions are derived for the preindustrial atmosphere/ocean system and for a system that may be similar to that existing during the Late Cretaceous Period. The model is then used to determine how the modern and ancient marine biospheres might be affected by an oceanic impact of a large asteroid or comet. Such an event could perturb the carbon cycle in several different ways. Global darkening caused by stratospheric dust veil could destroy most of the existing phytoplankton in a period of several weeks to several months. At the same time, dissolution of atmospheric NOx compounds synthesized during the impact would lower the pH of ocean surface waters and release CO2 into the atmosphere. Both effects might be enhanced by an influx of CO2 released from upwelling of deep ocean water near the hot impact site, from oxidation of dead organic matter, and from the comet itself. The net result could be to raise surface temperatures by several degrees and to make the surface ocean uninhabitable by calcareous organisms for as much as 20 yrs (the time scale for mixing with deep ocean). It appears unlikely, however, that an impact could create a "Strangelove ocean," in which surface waters remained corrosive to calcium carbonate for thousands or tens of thousands of years. Thus, disruption of the carbon cycle by an impact event cannot by itself explain the scarcity of calcium carbonate in sediments found within the first few centimeters above the K/T boundary.
Comments on the BLAG model: the carbonate-silicate geochemical cycle and its effect on atmospheric carbon dioxide over the past 100 million years
A self-consistent method of determining initial conditions for the model presented by Berner, Lasaga, and Garrels (1983) (henceforth, the BLAG model) is derived, based on the assumption that the CO2 geochemical cycle was in steady state at t = -100 my (million years). This initialization procedure leads to a dissolved magnesium concentration higher than that calculated by Berner, Lasaga, and Garrels and to a low ratio of dissolved calcium to bicarbonate prior to 60 my ago. The latter prediction conflicts with the geologic record of evaporite deposits, which requires that this ratio remain greater than 0.5. The contradiction is probably caused by oversimplifications in the BLAG model, such as the neglect of the cycles of organic carbon and sulfur.
The radium content of varved clay and a possible age of the Hartford, Connecticut, deposits
Evidence for a steepening of geothermal gradients in some deep mines and drill holes
Shell structure and distribution of Cloudina, a potential index fossil for the terminal Proterozoic
Cloudina-bearing biosparites and biomicrites in the lower part of the Nama Group, Namibia, contain a wide morphological diversity of shell fragments that can all be attributed to the two named species C. hartmannae and C. riemkeae. The curved to sinuous tubular shells of Cloudina were multi-layered. Each shell layer was 8 to 50 micrometers thick and in the form of a slightly flaring tube with one end open and the other closed. Growth appears to have been periodic with successive shell layers forming within older layers. Each added layer was slightly elevated from the previous layer at the proximal end and was asymmetrically placed within the older layer so that only a portion of the new shell layer was fused to the previous layer. This type of growth left a relatively large unminerialized area between the shell layers which was often partially or fully occluded by early marine cements. The thin shell layers exhibit both plastic and brittle deformation and were likely formed of a rigid CaCO3-impregnated organic-rich material. Often the shell layers are preferentially dolomitized suggesting an original mineralogy of high-magnesian calcite. Both species in the Nama Group formed thickets, or perhaps bioherms, and this sedentary and gregarious habit suggests that Cloudina was probably a filter-feeding metazoan of at least a cnidarian grade of organization. The unusual shell structure of Cloudina gives rise to a characteristic suite of taphonomic and diagenetic features that can be used to identify Cloudina-bearing deposits within the Nama Group and in other terminal Proterozoic deposits around the world. Species of Cloudina occur in limestones from Brazil, Spain, China, and Oman in sequences consistent with a latest Proterozoic age assignment. In addition, supposed lower Cambrian, pre-trilobitic, shelly fossils from northwest Mexico and the White-Inyo Mountains in California and Nevada, including Sinotubulites, Nevadatubulus, and Wyattia, are all either closely related to or con-generic with Cloudina. Hence, it is probable that these outcrops are latest Proterozoic in age, and that Cloudina or Cloudina-like organisms were widely distributed at that time. It is possible, moreover, to suggest that metazoan biomineralization occurred on a global scale by the latest Proterozoic, at the same time that evidence for complex multicellularity and locomotion in animals appears in siliciclastic "Ediacaran" rocks in the form of body and trace fossils.
A paleoweathering profile from Griqualand West, South Africa: evidence for a dramatic rise in atmospheric oxygen between 2.2 and 1.9 bybp
A core drilled near Wolhaarkop in Griqualand West, South Africa, intersected highly oxidized Kuruman Iron Formation below red beds of the Gamagara Formation. The lateral equivalents of the Kuruman Iron Formation in this drill hole consist largely of siderite, ankerite, magnetite, greenalite, and quartz. The oxidation of the Kuruman Iron Formation in WOL 2 occurred almost certainly during weathering prior to the deposition of the Gamagara Formation. The date of this weathering episode is bracketed between about 2.2 and 1.9 bybp by the age of the Ongeluk lavas in the Transvaal sequence below the unconformity and by the age of the Hartley lavas in the Olifantshoek Group above the unconformity. The ratio of iron to SiO2 in the several facies of the weathered Kuruman Iron Formation in WOL 2 is nearly the same as that in their unweathered equivalents. Since SiO2 loss during weathering was almost certainly minor, the similarity of the Fe/SiO2 ratio in the weathered and unweathered BIF indicates that nearly all the "FeO" in the Kuruman Iron Formation was oxidized and retained as FeO3 during weathering. Such a high degree of iron retention is best explained by an O2 content of the atmosphere > or = 0.03 atm at the time of weathering. Such an O2 pressure is very much greater than that suggested by the composition of paleosols developed on basalt > or = 2.2 bybp but is consistent with the highly oxidized nature of the 1.85 by Flin Flon paleosol. The new data suggest that PO2 rose dramatically from about 1 percent PAL (present atmospheric level) to > or = 15 percent PAL between 2.2 and 1.9 bybp.
Carbonate deposition during the late Proterozoic Era: an example from Spitsbergen
Carbonate sediments reflect the physico-chemical and biological circumstances of their formation; thus, features of limestones and dolomites may provide insights into both environmental and evolutionary change through geological time. The Upper Proterozoic (approx 800-700 Ma) Akademikerbreen Group, Spitsbergen, comprises 2000 m of carbonates, with only minor intercalations of quartz arenite and shale. Although Proterozoic carbonates are often seen as predominantly dolomitic, the Akademikerbreen Group is about 45 percent limestone. Stromatolites are conspicuous in outcrop but constitute only 25 percent of the total section. Micrites and coarser intraclastic carbonates derived mainly from micritric precursors comprise 60 percent of the group, while oolites make up the remaining 15 percent. Distinctive sedimentary features of the group include giant (up to 16 mm) ooids, very early diagenetic calcite nodules and cements, micrites containing subaqueous shrinkage cracks filled with equant microspar cement, and strong 13C enrichment in both carbonates and co-occurring organic matter. The principal features of Akademikerbreen carbonates are widely distributed in coeval successions. However, these rocks appear to differ from older limestones and dolomites in their relative abundance of grainstones and, perhaps, micrites, as well as their paucity of tufa-like laminates and columnar or coniform stromatolites that preserve petrographic evidence of in situ precipitation as a dominant means of carbonate accretion. Upper Proterozoic carbonates also differ from Paleozoic accumulations, but the transition is not abrupt. Most changes accompanying the Proterozoic/Phanerozoic transition can be interpreted in terms of the consequences rather than the causes of metazoan and metaphyte evolution, including the evolution of biomineralization. Carbonate sedimentology reinforces data from other sources which indicate the last 200 to 300 Ma of the Proterozoic Eon was a distinctive interval of Earth history.
