Biogeochemical modelling of the rise in atmospheric oxygen

Understanding the evolution of atmospheric molecular oxygen levels is a fundamental unsolved problem in Earth's history. We develop a quantitative biogeochemical model that simulates the Palaeoproterozoic transition of the Earth's atmosphere from a weakly reducing state to an O-2-rich state. The purpose is to gain an insight into factors that plausibly control the timing and rapidity of the oxic transition. The model uses a simplified atmospheric chemistry (parameterized from complex photochemical models) and evolving redox fluxes in the Earth system. We consider time-dependent fluxes that include organic carbon burial and associated oxygen production, reducing gases from metamorphic and volcanic sources, oxidative weathering, and the escape of hydrogen to space. We find that the oxic transition occurs in a geologically short time when the O-2-consuming flux of reducing gases falls below the flux of organic carbon burial that produces O-2. A short timescale for the oxic transition is enhanced by a positive feedback due to decreasing destruction of O-2 as stratospheric ozone forms, which is captured in our atmospheric chemistry parameterization. We show that one numerically self-consistent solution for the rise of O-2 involves a decline in flux of reducing gases driven by irreversible secular oxidation of the crust caused by time-integrated hydrogen escape to space in the preoxic atmosphere, and that this is compatible with constraints from the geological record. In this model, the timing of the oxic transition is strongly affected by buffers of reduced materials, particularly iron, in the continental crust. An alternative version of the model, where greater fluxes of reduced hydrothermal cations from the Archean seafloor consume O-2, produces a similar history of O-2 and CH4. When climate and biosphere feedbacks are included in our model of the oxic transition, we find that multiple 'Snowball Earth' events are simulated under certain circumstances, as methane collapses and rises repeatedly before reaching a new steady-state.