The procedure for the simulation of a time dependent sediment stratigraphy for 4 different sediment species with a global biogeochemical ocean general circulation model is presented. In this feasibility study, a sediment column for the last 131 kyr was accumulated through use of a variable ocean velocity field (mixed from preindustrial and glacial ocean time slices in a kinematic approach). The biogenic components of sediment are entirely a model product. Artificial sediment cores were recovered from the simulated stratigraphy. The approach is the foundation for a systematic comparison of paleoclimate simulations with the sediment record. Through use of adjoint model codes the direct 4-dimensional data assimilation of paleoclimate observations will become possible in the long range. ¿ 2001 American Geophysical Union