Iron and sulfur cycling in the cGENIE.muffin Earth system model (v0.9.21)

The coupled biogeochemical cycles of iron and sulfur are central to the long-term biogeochemical evolution of Earth's oceans. For instance, before the development of a persistently oxygenated deep ocean, the ocean interior likely alternated between states buffered by reduced sulfur ("euxinic") and buffered by reduced iron ("ferruginous"), with important implications for the cycles and hence bioavailability of dissolved iron (and phosphate). Even after atmospheric oxygen concentrations rose to modern-like values, the ocean episodically continued to develop regions of euxinic or ferruginous conditions, such as those associated with past key intervals of organic carbon deposition (e.g. during the Cretaceous) and extinction events (e.g. at the Permian-Triassic boundary). A better understanding of the cycling of iron and sulfur in an anoxic ocean, how geochemical patterns in the ocean relate to the available spatially heterogeneous geological observations, and quantification of the feedback strengths between nutrient cycling, biological productivity, and ocean redox requires a spatially resolved representation of ocean circulation together with an extended set of (bio)geochemical reactions.