In this study, we used MITgcm-REcoM2 to simulate the stepwise glacial-interglacial atmospheric pCO2 change. A general description of the model can be found in the supporting material for Hauck et al. 2013. There are seven simulations included in this dataset: two control runs for interglacial and glacial conditions (IG_ctl, G_ctl); three region-specific sensitivity runs(IG_Gso, IG_Gna, IG_Gns); and two simulations regarding the glacial iron fertilization in the Southern Ocean. Dissolved Inorganic Carbon (DIC) and Export Production in all runs are included in this study, and the potential temperature (THETA), salinity (SALT), meridional velocity (VVEL), sea ice concentration(SIarea), air-sea surface pCO2 (dpCO2surface) difference are available in IG_ctl, IG_Gso, IG_Gna, IG_Gns, and G_ctl.