The simulations were recalculated based on the author's corrections Wu et al. (2025). Koehler_NH_iceretreat_over_shales_20ka.dat: Ice retreat over shales in North America is calculated in the area 180°W-75°W, 45°N-90°N based on maps available in 1°x1° spatial resolution of: 1) land ice: ICE6G after Peltier, et al. (2015) lithology from Amiotte Suchet et al. (2003). Koehler_cfluxes_20ka.dat: C fluxes (oxidation of pedrogenic organic carbon) from ice retreat over shales in North America. Koehler_atmD14C_20ka_correction.dat, new simulations: Simulated changes in atmospheric Delta14C (permil) using the BICYCLE-SE mode. (process A): caused by release oxidation of pedrogenic organic carbon, (process B): caused by Arctic permafrost mobilization (Winterfeld et al., 2018). Koehler_atmCO2_20ka_corrected.dat: new simulations: Simulated changes in atmospheric CO2 (ppm) using the BICYCLE-SE model. (process A): caused by release oxidation of pedrogenic organic carbon, (process B): caused by Arctic permafrost mobilization (Winterfeld et al., 2018).