For sediments encompassed by the rhizosphere, sulfide, sulfate, methane, DOC, and hydrogen ion concentration profiles were determined iteratively prior to model implementation using a central difference scheme: