Even though models of water quality have become increasingly detailed over time, their applicability to analyze long-term effect on sedimentary bed formation is yet to be clarified. Hence, an integrated, layer-resolved, process-based, sediment-water coupled, long-term robust, three-dimensional (3D) ecosystem model, including realistic sedimentary and pelagic processes, was developed. The constituents of the integrated model included a multi layered ecosystem model, a quasi-three dimensional hydrodynamic model, a wave hindcasting model, and a sediment deposition and erosion model. Because numerical modeling difficulties arise in setting the initial conditions, especially for a sedimentary bed, this novel modeling approach suggests to initiate the model with the simplest initial conditions of no particulate organic carbon content (POCC) and uniform inorganic sediment distribution. The model was applied to Tokyo Bay and quasi-steady state POCC in the sediment was obtained through long-term computation with realistic sedimentary and pelagic processes. Wave induced bed shear stress (WBSS) is critical for sediment erosion and spatial movement of sediment pollution, while a stable morphology is determined by the balance between the current induced bed shear stress (CBSS) and WBSS in Tokyo Bay. This novel modeling approach, with the simplest sedimentary initial conditions and realistic sedimentary and pelagic processes, provides a great tool for long-term ecosystem modeling in future studies.