The stable carbon isotopic composition (delta C-13) is an important variable to study the ocean carbon cycle across different timescales. We include a new representation of the stable carbon isotope C-13 into the HAMburg Ocean Carbon Cycle model (HAMOCC), the ocean biogeochemical component of the Max Planck Institute Earth System Model (MPI-ESM). C-13 is explicitly resolved for all oceanic carbon pools considered. We account for fractionation during air-sea gas exchange and for biological fractionation epsilon(p) associated with photosynthetic carbon fixation during phytoplankton growth. We examine two epsilon(p) parameterisations of different complexity: epsilon(Popp)(p) varies with surface dissolved CO2 concentration (Popp et al., 1989), while epsilon(Laws)(p) additionally depends on local phytoplankton growth rates (Laws et al., 1995). When compared to observations of delta C-13 of dissolved inorganic carbon (DIC), both parameterisations yield similar performance. However, with regard to delta C-13 in particulate organic carbon (POC) epsilon(Popp)(p) shows a considerably improved performance compared to epsilon(Laws)(p). This is because epsilon(Laws)(p) produces too strong a preference for C-12, resulting in delta C-13(POC) that is too low in our model. The model also well reproduces the global oceanic anthropogenic CO2 sink and the oceanic C-13 Suess effect, i.e. the intrusion and distribution of the isotopically light anthropogenic CO2 in the ocean. The satisfactory model performance of the present-day oceanic delta C-13 distribution using epsilon(Popp)(p) and of the anthropogenic CO2 uptake allows us to further investigate the potential sources of uncertainty of the Eide et al. (2017a) approach for estimating the oceanic C-13 Suess effect. Eide et al. (2017a) derived the first global oceanic C-13 Suess effect estimate based on observations. They have noted a potential underestimation, but their approach does not provide any insight about the cause. By applying the Eide et al. (2017a) approach to the model data we are able to investigate in detail potential sources of underestimation of the C-13 Suess effect. Based on our model we find underestimations of the C-13 Suess effect at 200 m by 0.24 parts per thousand in the Indian Ocean, 0.21 parts per thousand in the North Pacific, 0.26 parts per thousand in the South Pacific, 0.1 parts per thousand in the North Atlantic and 0.14 parts per thousand in the South Atlantic. We attribute the major sources of underestimation to two assumptions in the Eide et al. (2017a) approach: the spatially uniform preformed component of delta C-13(DIC) in year 1940 and the neglect of processes that are not directly linked to the oceanic uptake and transport of chlorofluorocarbon-12 (CFC-12) such as the decrease in delta C-13(POC) over the industrial period. The new C-13 module in the ocean biogeochemical component of MPI-ESM shows satisfying performance. It is a useful tool to study the ocean carbon sink under the anthropogenic influences, and it will be applied to investigating variations of ocean carbon cycle in the past.