Wetlands can remove organic contaminants, metals and radionuclides from wastewater through various biogeochemical mechanisms. In this study, a mathematical model was developed for simulating the fate and transport of chemical species in marsh wetland sediments. The proposed model is a one-dimensional vertical saturated model which is incorporated advection, hydrodynamic dispersion, biodegradation, oxidative/reductive chemical reactions and the effects from external environments such as the growth of plants and the fluctuation of water level due to periodic tides. The tidal effects causes periodic changes of porewater flow in the sediments and the evapotranspiration and oxygen supply by plant roots affect the porewater flow and redox condition on in the rhizosphere along with seasonal variation. A series of numerical experiments under hypothetical conditions were performed for simulating the temporal and spatial distribution of chemical species of interests using the proposed model. The fate and transport of a trace metal pollutant, chromium, in marsh sediments were also simulated. Results of numerical simulations show that plant roots and tides significantly affect the chemical profiles of different electron acceptors, their reduced species and trace metals in marsh sediments.