A zonally averaged, circulation-biogeochemical ocean model is used to explore how the distribution of PO4 and δ13C in the major basins and the atmospheric pCO2 respond to rapid changes in the thermohaline circulation (THC). Different evolutions of the Atlantic THC are simulated by applying surface freshwater pulses typical, for example, of Heinrich events and the last deglaciation. In the model, when the THC completely collapses, PO 4 increased (>0.5 mmol m-3) and δ13C decreased (<0.5‰) in Atlantic bottom waters because of a drop in ventilation by North Atlantic Deep Water (NADW). Although consistent with the traditional interpretation of sedimentary records of benthic foraminiferal Cd/Ca and δ13C rises (>1‰) because of a longer residence time of waters at the surface. These surface anomalies are propogated to the bottom North Atlantic when the THC resumes. When the thermohaline overturning is only partly reduced and at shallower depths, changes in bottom waters are accompanied by a PO4 decrease and δ13C increase at intermediate levels in the mid-latitude Atlantic. This results in enhanced vertical gradients of these properties consistent with chemical and isotopic reconstructions for the last glacial maximum. Finally, the atmospheric pCO2 increases during the cold period in the North Atlatic when the THC is reduced with an amplitude (7-30 μatm) and times scale (~102 to 1-2 × 103 yr) depending on the intensity of the THC change. This is qualitatively consistent with recent data from an Antarctic ice core documenting a pCO2 increase during the Yougner Dryas and after Heinrich events four and five.