Understanding Ferruginous Versus Euxinic Conditions by Simulating Microbial Conditions in Meromictic Lakes.

Ferruginous (iron-rich) conditions have been prominent in oceans throughout the Earth's geologic history but now are reliably found only in a handful of permanently stratified lakes. Microbially mediated iron reduction in such anoxic environments competes with sulfate reduction, which promotes euxinic (sulfide-rich) conditions. Besides the shared demand for organic compounds, the competition is fostered by the produced hydrogen sulfide, which may reduce iron oxides abiotically or co-precipitate with dissolved iron as iron sulfides. Understanding why some environments develop ferruginous rather than euxinic conditions (or vice versa), as well as the attendant effects on methanogenic fermentation, is key to understanding both modern and ancient anoxic ecosystems. Here, we reproduce biogeochemical distributions in multiple anoxic, low-sulfate, meromictic lakes around the world using a biomass-explicit reaction-transport model with a fixed set of metabolism-specific microbial parameters. The results suggest that sulfate reduction and methanogenesis are ubiquitous even in iron-rich systems, and are reflected in microbial surveys. Ferruginous conditions typically develop for surface sulfate concentrations below ≃100 μM. Interestingly, there seems to be a dearth of stably stratified water bodies where sulfate concentrations can persist in the medium-sulfate range of several hundred μM. Rather, when sulfur burial into the sediments becomes iron limited, sulfate tends to accumulate in the water column to much higher (mM) concentrations. A similar mechanism could be suggested to have operated in the variably sulfidic and ferruginous water columns of early oceans. Model simulations also reveal the previously underappreciated role of physical transport in shaping biogeochemical distributions, as minor variations in mixing rates can lead to large variations in microbial abundances. Model applicability across multiple lakes points to an encouraging possibility that geochemical patterns in complex biogeochemical systems may be described from a small number of thermodynamic and kinetic principles using a minimum of fitting parameters.