We combine consistently dated benthic carbon isotopic records distributed over the entire Atlantic Ocean with numerical simulations performed by a glacial configuration of the Norwegian Earth System Model with active ocean biogeochemistry, in order to interpret the observed Cibicides δ13C changes at the stadial-interstadial transition corresponding to the end of Heinrich Stadial 4 (HS4) in terms of ocean circulation and remineralization changes. We show that the marked increase in Cibicides δ13C observed at the end of HS4 between ~2000 and 4200 m in the Atlantic can be explained by changes in nutrient concentrations as simulated by the model in response to the halting of freshwater input in the high latitude glacial North Atlantic. Our model results show that this Cibicides δ13C signal is associated with changes in the ratio of southern-sourced (SSW) versus northern-sourced (NSW) water masses at the core sites, whereby SSW is replaced by NSW as a consequence of the resumption of deep water formation in the northern North Atlantic and Nordic Seas after the freshwater input is halted. Our results further suggest that the contribution of ocean circulation changes to this signal increases from ~40 % at 2000 m to ~80 % at 4000 m. Below ~4200 m, the model shows little ocean circulation change but an increase in remineralization across the transition marking the end of HS4. The simulated lower remineralization during stadials than interstadials is particularly pronounced in deep subantarctic sites, in agreement with the decrease in the export production of carbon to the deep Southern Ocean during stadials found in previous studies.