Abstract:The reaction-transport model is widely used to identify and quantify dissolved chemical species in sediment porewaters. In this paper, a modularized code framework of diagenetic model was proposed as a diagnostic tool to fit the porewater profiles in marine sediments. Based on the conservation principle of the finite volume method, we combined the discretized diagenetic equations with various geochemical reactions, including but not limited to methanogenesis, sulfate reduction, etc.. The code was organized in object-oriented FORTRAN and verified with literature parameters, which proved its robustness and effectiveness. Finally, three different sites（IODP Expedition 311 Site U1327, UBGH2-1 1 , ODP Leg204 Site 1245） are exemplified as case studies.

An Object-Oriented Diagnostic Model for the Quantification of Porewater Geochemistry in Marine Sediments