The mathematics of early diagenesis: from worms to waves

Abstract. The changes that sediments undergo after deposition are collectively known as diagenesis. Diagen-esis is not widely recognized as a source for mathemat-ical ideas; however, the myriad processes responsible for these changes lead to a wide variety of mathematical models. In fact, most of the classical models and meth-ods of applied mathematics emerge naturally from quan-tification of diagenesis. For example, small-scale sedi-ment mixing by bottom-dwelling animals can be described by the diffusion equation; the dissolution of biogenic opal in sediments leads to sets of coupled, nonlinear, ordinary differential equations; and modeling organisms that eat at depth in the sediment and defecate at the surface suggests the one-dimensional wave equa-tion, while the effect of waves on pore waters is governed by the two- or three-dimensional wave equation. Diage-netic modeling, however, is not restricted to classical methods. Diagenetic problems of concern to modern mathematics exist in abundance; these include free-boundary problems that predict the depth of biological mixing or the penetration of O2 into sediments, algebra-ic-differential equations that result from the fast-revers-ible reactions that regulate pH in pore waters, inverse calculations of input functions (histories), and the deter-mination of the optimum choice in a hierarchy of pos-sible diagenetic models. This review highlights and ex-plores these topics with the hope of encouraging further modeling and analysis of diagenetic phenomena. 1