ABSTRACT

Iron isotope compositions of sedimentary pyrite are widely used to trace redox and biogeochemical processes, but the isotopic fractionation during FeS-to-pyrite conversion in early diagenesis processes remains poorly constrained under natural sedimentary environments. Here we employ a one-dimensional reaction–transport model that explicitly resolves Fe isotope dynamics during early diagenesis in sediments of the anoxic Santa Barbara Basin, a well-studied modern site of active iron–sulfur cycling. The model reproduces the observed pore-water and solid-phase chemical and isotopic depth profiles and constrains the net isotope fractionation associated with FeS-to-pyrite conversion to approximately −0.5‰, irrespective of the Fe isotope fractionation factors of FeS precipitation/dissolution and dissimilatory iron reduction. This small net fractionation may reflect slow pyrite precipitation in the sediment, which may allow positive equilibrium isotope effects to offset the negative kinetic isotope effect associated with pyrite formation. These results would provide a framework for the broader interpretation of iron isotope records in anoxic sediments.