Abstract. We have developed and employed an Earth system model to explore the forcings of atmospheric pCO 2 change and the chemical and isotopic evolution of seawater over the last glacial cycle. Concentrations of dissolved phosphorus (DP), reactive nitrogen, molecular oxygen, dissolved inorganic carbon (DIC), total alkalinity (TA), 13 C-DIC, and 14 C-DIC were calculated for 24 ocean boxes. The bi-directional water fluxes between these model boxes were derived from a 3-D circulation field of the modern ocean (Opa 8.2, NEMO) and tuned such that tracer distributions calculated by the box model were consistent with observational data from the modern ocean. To model the last 130 kyr, we employed records of past changes in sea-level, ocean circulation, and dust deposition. According to the model, about half of the glacial pCO 2 drawdown may be attributed to marine regressions.