AbstractThe beginning and end of the Proterozoic Eon are marked by extreme variations in carbonate carbon isotope values that have been interpreted to record massive perturbations to the global carbon cycle. The lower Proterozoic contains an extended interval of strata characterized by positive carbonate δ13C values. Conversely, uppermost Proterozoic carbonate strata contain thick intervals with extremely negative δ13C values and multiple large swings in carbonate δ13C. Previous attempts to model these pronounced carbon isotope excursions as shifts in the global marine dissolved inorganic carbon (DIC) reservoir have proved to be problematic, as the direction and magnitude of these positive and negative carbon isotope excursions require unrealistic amounts of either organic carbon burial or organic carbon oxidation, respectively. Here we present a modified global carbon cycle model—coupled with oxygen and sulfur cycle mass balances—that includes a parameterization of the recycling of sedimentary isotope anomalies and allows the extent of organic carbon oxidation to vary as a function of atmospheric oxygen levels. Our model is designed to match carbon isotope records while maintaining redox and mass balance with a given set of initial conditions and carbon cycle parameterizations. Using this approach, we demonstrate that there is a range of plausible biogeochemical perturbations that could induce substantial δ13C excursions in the global marine DIC reservoir. However, we also find that there are multiple, nonunique Earth system states for any observed marine δ13C value.