Origin and significance of ultra-slow calcite dissolution rates in deep sea sediments

Abstract

Prediction of mineral-fluid reaction rates in geologic materials is subject to large uncertainties even though the rates are critical for many applications, from studies of past climate change to the engineering of carbon removal and storage. Deep sea carbonate sediments, which have been extensively cored and characterized at dozens of sites by ocean drilling programs, provide unique opportunities to examine how minerals and fluids react over multi-million-year timescales in natural conditions at low temperature. It has already been established that the calcite-fluid reaction in these systems is extremely slow, but the explanation for the ultra-slow rates, and the fact that non-zero rates persist for tens of millions of years, remains elusive. We extend the analysis of pore fluid Sr and Ca concentration data in the literature from ocean drilling archives to 21 drill sites in sediment sections where the fraction of carbonate is larger than 80% to systematically estimate the rates at which calcite dissolves and precipitates as a function of depth and sediment age. Pore fluid Sr/Ca used in this study, and Sr and Ca isotopes in other published analyses, provide estimates of the gross calcite dissolution rates, which are mostly balanced by secondary calcite precipitation. Pore fluid Ca concentration data provide key, previously underappreciated, information on the net calcite dissolution rates. Using these data we evaluate simple models for calcite reaction kinetics and pore fluid sediment evolution and conclude that (1) reaction rates are extremely slow because the reactive site density of the bulk sediment decreases continuously and systematically approximately as age−1, and (2) reaction persists to hundreds of meters burial depth and tens of millions of years because pore fluids never become closed due to diffusive communication with the oceans, and hence remain slightly undersaturated with respect to calcite. The systematic decrease in bulk reactive site density begins immediately upon deposition, greatly slowing the rate of dissolution and hence affecting models of calcite dissolution on the seafloor. The failure to reach equilibrium after millions of years and hundreds of meters of burial, is a consequence of extremely slow net reaction rates that cannot drive the pore fluids to equilibrium against the diffusive loss of Ca and DIC to the ocean. The existence of two or more calcite phases of different solubility is a likely source of a long-lived, but small, driving force for reaction, but its effect decreases with continued recrystallization and pore fluid evolution.