Kobayashi and colleagues present a novel offline simulation conducted with a carbon isotope-enabled biogeochemical model, forced with climate data from a transient MIROC4m simulation of the last deglaciation. Their investigation focuses on the relationship between changes in atmospheric CO₂, stable and radiocarbon, and the varying AMOC during that period. The simulation reveals relatively minor changes in atmospheric CO₂ concentration compared to ice core reconstructions, while demonstrating that changes in water mass ventilation and sourcing align with proxy reconstructions. The authors further analyze the simulated changes in pCO₂, attributing them to physical (temperature and salinity) and biogeochemical (dissolved inorganic carbon and alkalinity) drivers, revealing complex interactions and compensating effects.

The manuscript is well-written and well-illustrated. While similar studies have been previously conducted with intermediate complexity models, this study represents a significant step towards comprehensive transient simulations with an AOGCM, despite not fully achieving this here. The authors transparently acknowledge certain critical factors during the last deglaciation, such as sea-level rise, ice sheets, and Southern Ocean sea surface temperature biases, which were not accounted for in this study. Considering all these processes for a deglacial simulation is very challenging in such a complex model, making it understandable that they are not fully considered here. Therefore, I recommend the manuscript for publication in Climate of the Past after minor revisions. Detailed comments outlining specific areas for improvement are provided below.

L4: Better introduce the abbreviation AMOC already here than in L7.

L5: Here and in the following, I would replace “atmospheric partial pressure of carbon dioxide” just with “atmospheric concentration of CO₂”.

L30: deglacial period not deglaciation period.

L43: investigate not inter.

L80: What is meant by vertical one-dimensional distribution in this context? Previous studies have evaluated the 3D distribution of data.

L87-92: The setup of the model is not fully clear to me. In line 87 it is mentioned that a BGC model was coupled to an ocean model, but in line 90 it is stated that the BGC model was forced with MIROC output. Does this mean that the coupled BGC-ocean model was forced with atmospheric boundary conditions of MIROC? Later on it reads to me as if the BGC model is run entirely in offline mode. Can the authors clarify this in section 2.1?

L91: Can you briefly mention how the MIROC transient simulation was forced, e.g., freshwater fluxes GHG concentrations etc.?

L110: Why were pre-industrial and not LGM values used to initialize these atmospheric values?

L138: Can you give numbers for the AMOC strengths during the LGM and the Holocene?

L146: According to Fig. 2j the Pacific appears to me simulated in rather good agreement with the data. Certainly, much better than the Southern Ocean.

L151: Maybe it is worth mentioning, that indeed the agreement of Kobayashi2021 is much better in the deep Southern Ocean, but quite a bit worse for the mid-depths. It therefore appears to me that it is not as simple as including these processes.

L160: can you give a reason why the SST difference between the LGM and Holocene is so small compared to observations?

L178: There appears to be very little change during the transition from the LGM to HS1, which makes sense, since the AMOC also remains virtually constant. Maybe this needs to be hence slightly rephrased.

L203: As mentioned before, the difference between the LGM and HS1 seems very small in the simulation, because the AMOC also changes very little. In contrast, a larger decrease is observed in the data.

L209: Mention that this is most pronounced in the Atlantic.

L210: It is important to note, that two different things are compared here. The model output has an annual resolution and therefore shows the “true” perturbation magnitude. On the other hand, the reconstructions from marine sediments are smoothed out by processes such as bioturbation, coring artifacts, etc. It is therefore expected that the signal amplitude is bigger in the model than the data for such short perturbations like the YD.

L211: Maybe explicitly mention that this is for the reconstructions.

L224: This is rather surprising to me, as the AMOC strength actually doesn’t change much, but the change in carbon export is rather large in the South Pacific. I’m therefore wondering whether this can really be attributed to an AMOC weakening, or whether other processes dominate this effect? From the pattern and the fact that this negative anomaly persists throughout the deglaciation independent of the AMOC strength, suggest to me that this is primarily a signal of increased iron limitation, which is mentioned in the text.

L227: By which mechanism propagate these changes to the North Pacific?

L304-310: It could also be that the AMOC weakening is too strong in the model.

L334-344: The terrestrial biosphere plays an important role for atmospheric d13C, which I think should be mentioned here as well (see e.g., Jeltsch-Thoemmes et al., 2019, doi:10.5194/cp-15-849-2019).

Fig. 1: I understand the intention to make the model and data timeseries overlap for better comparability. However, I find this in panel c somewhat misleading, as there is a factor of more than two between both y-axes. I would like to see the same increment for both axes as it is done in panel b. Further, can a panel showing global mean surface temperature be added, for instance compared to the data assimilation by Osman et al., 2021 (doi:10.1038/s41586-021-03984-4)?

Fig. 3 and 4: Can you add a similar figure like these two, but for DIC (anomaly?), to better illustrate where the carbon was stored in the course of the deglaciation?

Figs. S4 and S5: The panels are very small and hence hard to read. Can they be made larger, splitting the panels into different rows or columns?

Fig. S7 to S9. I find both these figures very instructive and if possible would like to see them in the main text.

See attached file