摘要
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.
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.
基金
supported by the National 973 Project (No. 2013CB835000)
the National Natural Science Foundation of China (Nos. 41230102 and 40903002)
National Key Scientific Project organized by the China Geological Survey (No. GZH201200305-06-04)
