The oxygen geochemical cycle: dynamics and stability

The first and possibly only major rise of atmospheric oxygen, from p(O2) no more than 0.1 percent the present atmospheric level (PAL) to p(O2) at least 10 percent PAL, appears to have occurred sometime before 2 Ga ago, although the exact time of and the cause(s) for the rise have been hotly debated. Equally important questions on the atmospheric oxygen concern its stability, especially the mechanisms regulating the atmospheric p(O2) level and the causes and magnitude of p(O2) variations since the first major rise of atmospheric oxygen. Previous efforts to model the p(O2) variation during the Phanerozoic time have typically relied on secondary information, such as the carbon and sulfur isotopic records of sedimentary rocks, and on simple dynamics of the geochemical cycles of O, C, S, and P based on box-type models. As a result, many kinetic questions about the variation and stability of atmospheric oxygen could not have been answered. Here, we quantitatively evaluate the dynamics and stability of atmospheric O2 and CO2, using recent experimental data, field observations, and a new model for the C-O coupled geochemical cycles. We examine the change with time in the fluxes of various compounds (O2, CO2, phosphate, organic C, carbonate C, C-bearing reduced volcanic gases, and C-free reduced volcanic gases) among the various reservoirs (atmosphere, soil, surface ocean, deep ocean, the lower crust and mantle, and upper crust) under a variety of scenarios. Our model does not assume steady-state fluxes for any of the reservoirs. Rather, the model incorporates the kinetic experimental data on oxidation of coal, a proxy for kerogen, the dynamics of soil formation and erosion, the kinetics of the decomposition of organic matter in the Oceans by aerobic and anaerobic bacteria, the equilibrium ocean-atmosphere carbonate model, the observed relationships among the organic burial flux, dissolved O2 content of deep ocean, and sedimentation rates, and the three-box model ocean. The important parameters that strongly influence the dynamics of atmospheric O2 are found to be (1) the total area of soil formation on Earth, (2) the average soil depth, (3) the average rate of physical erosion of soils, which is linked to the average rate of accumulation of clastic sediments in the oceans, (4) the composition and flux of volcanic gas, and (5) the level of atmospheric CO2. We develop kinetic equations linking these parameters to the production and consumption fluxes of atmospheric