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, PO4 increases (>0.5 mmol m-3) and Δ13C decreases (1?) because of a longer residence time of waters at the surface. These surface anomalies are propagated 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 Atlantic when the THC is reduced with an amplitude (7--30 μatm) and timescale (~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 Younger Dryas and after Heinrich events 4 and 5. ¿ 1998 American Geophysical Union