Numerical models for diagenesis and the Neogene Sr isotopic evolution of seawater from DSDP Site 590B

A numerical model for the diagenetic exchange of Sr between carbonates and their pore fluids during sedimentation and compaction has been developed. The model has been applied to data from DSDP Site 590B in order to assess the accuracy with which the Sr isotope record in the carbonate sediment reflects that of seawater. The most important process affecting the Sr in the solid carbonate is exchange with the pore fluid due to solution-reprecipitation, but the concentration or isotopic composition of Sr in the solid itself gives little or no information as to the magnitude of this exchange. The key to determining the rate of exchange is the pore fluid, where the variations of Sr²⁺ and⁸⁷Sr/⁸⁶Sr with depth are very sensitive indicators. The logical structure of applying the model to data from DSDP 590B is to find by successive iteration an ocean history (i.e., the initial⁸⁷Sr/⁸⁶Sr and Sr concentration of each increment of carbonate deposited) and a rate of Sr exchange between pore water and solid carbonate such that the model matches the present Sr concentration and⁸⁷Sr/⁸⁶Sr of both pore water and solid carbonate.Once all the data are matched, the model provides an estimate of the rate of Sr exchange due to solution-reprecipitation and the evolution of⁸⁷Sr/⁸⁶Sr in seawater over the past 20 million years. For DSDP 590B we find that solution-reprecipitation decreases rapidly with depth, from a near surface value of about 10% per million years to about 1% per million years below 200 m. This rate of exchange of Sr results in the carbonates of DSDP 590B preserving an accurate record of the Sr isotopic evolution of the ocean over the past 5 million years, but for ages greater than 5 million years the⁸⁷Sr/⁸⁶Sr ratio of the carbonate is systematically displaced from that of the seawater in which it was deposited. The maximum difference is of order 5 × 10⁻⁵.