CO2 Lag as Glacial Cooling Initiates

Introduction

Anyone who examines the Antarctic ice core data for temperature and CO₂ carefully will be forced to noctice that during the Emian (Marine Isotope Stage 5e) temperature peaked and CO₂ reached its maximum for that period and then for 13,000 years CO₂ held that maximum almost constant while temperature fell 7 degrees C - see Figure 1.

This is probably the most dramatic indication in the ice core record showing temperature falling while CO₂ holds high and constant. Can similar behavior be detected in association with other interglacials?

Figure 1: Here CO₂ is seen to hold Constant for 13,000 Years while Temperature Falls 7 Degrees C. Figure taken from [1].

Methodology and Harmonic Fits

The analysis employs harmonic fits to Antarctic ice core proxy temperature and CO₂ data using a greedy algorithm that successively adds periodic components to minimize residual variance [2]. Temperature data are from the EPICA Dome C deuterium-derived temperature proxy on the EDC3 chronology [3,4]. CO₂ data are from the Antarctic composite record, converted from AICC2012 to EDC3 chronology [5]. This approach achieves R² > 0.95 for both temperature and CO₂ over the 350,000-year record. Importantly, the algorithm independently recovers the canonical Milankovitch orbital periods (approximately 100,000, 41,000, and 23,000 years) without prior specification, validating both the methodology and the orbital pacing of glacial cycles [6]. Figure 2 displays the temperature record with harmonic fit and Marine Isotope Stage (MIS) codes. The major interglacial peaks are labeled: MIS 9 (~330,000 years ago), MIS 7e (~240,000 years ago), MIS 5e (the Eemian, ~125,000 years ago), and MIS 1 (the current Holocene). The MIS 7 complex includes several substages representing oscillations within that interglacial period. Figure 3 shows the corresponding CO₂ record with harmonic fit. The CO₂ peaks correlate with temperature peaks but, as we will demonstrate, with systematic temporal offsets that reveal the asymmetric nature of the carbon cycle response.

Figure 2: Temperature Record from date -350,000 to 1911 with Negative dates Indicating BCE. Figure taken from [7]. Marine Isotope Stage (MIS) Codes added

Figure 3: Atmospheric CO₂ Record from date -350,000 to 1911 with Negative dates Indicating BCE. Figure taken from [7].

CO2Versus Temperature Evidence for Asymmetric Response

Rather than examining temperature and CO₂ as separate time series, we can plot CO₂ directly against temperature to reveal the relationship between these variables. Figure 4 presents this composite plot using the evaluation of harmonic fits from date −350,000 to 1850. The plot exhibits pronounced curvature and hysteresis, indicating that CO₂ does not respond symmetrically to warming and cooling.

If CO₂ and temperature were linearly related with a simple phase lag, the parametric plot would trace an ellipse. Instead, the observed pattern shows extended horizontal excursions during cooling phases—periods where temperature decreases substantially while CO₂ remains elevated. This is the graphical signature of the asymmetric absorption mechanism: CO₂ release during warming is relatively rapid, while CO₂ uptake during cooling is rate-limited.

Figure 4: Plot of CO₂ Versus Temperature from Higginbotham 2025b. This Plot is a Composite of the Evaluation of the Fits to CO₂ and Temperature from date -350,000 to 1850. This Plot shows a Curvature in the Composite Plot

Analysis by Interglacial Period

To examine whether the Eemian CO₂ plateau is unique or representative of a general pattern, Figure 5 separates the composite CO₂-temperature plot into four panels corresponding to different interglacial periods. Each panel distinguishes the warming branch (labeled 'a') from the subsequent cooling branch (labeled 'b').

MIS 9 (−350k to −253k)

The upper left panel shows MIS 9, spanning approximately −350,000 to −253,000 years. The warming phase (a) from −350k to −332k shows CO₂ rising from approximately 185 ppm to nearly 300 ppm as temperature increases from −10°C to +4°C relative to present. The cooling phase (b) from −332k to −253k exhibits the characteristic horizontal excursion: temperature begins falling while CO₂ initially holds near its peak before beginning a slower decline. This panel provides clear evidence that the delayed CO₂response during cooling is not unique to the Eemian.

MIS 7e Through MIS 7d (−253k to −221k)

The upper right panel captures MIS 7e and the subsequent MIS 7d substage. Unlike MIS 9 and MIS 5e, this panel does not show a pronounced CO₂ plateau during cooling. Instead, CO₂ appears to track temperature downward more closely as the system transitions from MIS 7e into MIS 7d. This anomalous behavior may have a geophysical explanation worth considering. Research on speleothem records from the Austrian Alps indicates that MIS 7d coincided with the lowest Northern Hemisphere summer insolation of the past 800,000 years—387 W/m² at 65°N, caused by a rare combination of minimum obliquity and maximum precession [8]. This extreme orbital configuration drove rapid, hemisphere-wide cooling: atmospheric CO₂ fell below 203 ppm, North Atlantic sea surface temperatures dropped to near-glacial levels, and European pollen records indicate glacial-like conditions.

MIS 7c through MIS 6 (−221k to −137k)

The lower left panel spans the later substages of MIS 7 and the transition into MIS 6 (the Saalian glaciation). This period follows the anomalous MIS 7d interval and shows more complex behavior. Interestingly, the cooling branch (b) from −209k to −137k does appear to show horizontal excursions where CO₂ holds relatively constant while temperature declines—but at a cooler baseline than observed in MIS 9 or MIS 5e. The entire pattern is shifted to temperatures around −1°C to −7°C and CO2 levels in the 230–260 ppm range.

One possible interpretation is that after the extreme MIS 7d cooling event ended, the carbon cycle resumed its normal asymmetric behavior, but operating within the cooler, lower-CO₂ regime characteristic of MIS 7. The "failed interglacial" nature of MIS 7—which never achieved the warmth of MIS 5e or MIS 9 — may have compressed the entire CO₂-temperature relationship into a narrower range while preserving the underlying asymmetric dynamics.

This interpretation remains conjectural. The complexity of the MIS 7 substage oscillations makes it difficult to isolate clean warming and cooling branches comparable to those in MIS 9 and MIS 5e.

MIS 5e through MIS 1 (−137k to +1850)

The lower right panel presents the most recent and best-resolved portion of the record, including the Eemian (MIS 5e) and the current Holocene (MIS 1). The warming phase (a) from −137k to −127k shows the sharp rise into the Eemian peak, with CO₂ and temperature rising together along a steep trajectory.

The subsequent cooling branch (b) from −127k to −110k displays the most dramatic example of the phenomenon under investigation. Temperature falls approximately 7°C (from +4°C to −3°C) while CO₂ barely changes from its peak near 280 ppm. This is the 13,000- year Eemian plateau—the most striking demonstration in the ice core record that CO₂ does not immediately follow temperature during glacial cooling.

The clarity of this signal likely reflects both the magnitude of the Eemian warmth (creating a large CO₂"excess" requiring reabsorption) and the superior temporal resolution of the more recent ice core data. The Eemian thus serves as the type example of asymmetric carbon cycle response, with MIS 9 providing independent confirmation that this behavior is a recurring feature of glacial cycles.

Figure 5: Here the Plot from the Previous Figure is Split into Four Separate Plots. Upper Left Shows the MIS 9 Rise in Temperature (a) and Subsequent fall in Temperature (b). Upper Right Shows the MIS 7e Rise in Temperature (a) followed by fall in Temperature through MIS 7d (b). Lower Left Shows MIS 7c through MIS 6 with Rising Temperature (a) and Falling and Rising Temperature (b). Lower Right Shows MIS 5e through MIS 1 with Rising Temperature of MIS 5e (a) and Falling and Rising Temperature following (b)

Analysis by Interglacial Period

To examine whether the Eemian CO2 plateau is unique or representative of a general pattern, Figure 5 separates the composite CO2-temperature plot into four panels corresponding to different interglacial periods. Each panel distinguishes the warming branch (labeled 'a') from the subsequent cooling branch (labeled 'b').

MIS 9 (−350k to −253k)

The upper left panel shows MIS 9, spanning approximately −350,000 to −253,000 years. The warming phase (a) from −350k to −332k shows CO2 rising from approximately 185 ppm to nearly 300 ppm as temperature increases from −10°C to +4°C relative to present. The cooling phase (b) from −332k to −253k exhibits the characteristic horizontal excursion: temperature begins falling while CO2 initially holds near its peak before beginning a slower decline. This panel provides clear evidence that the delayed CO2 response during cooling is not unique to the Eemian.

MIS 7e Through MIS 7d (−253k to −221k)

The upper right panel captures MIS 7e and the subsequent MIS 7d substage. Unlike MIS 9 and MIS 5e, this panel does not show a pronounced CO2 plateau during cooling. Instead, CO2 appears to track temperature downward more closely as the system transitions from MIS 7e into MIS 7d. This anomalous behavior may have a geophysical explanation worth considering. Research on speleothem records from the Austrian Alps indicates that MIS 7d coincided with the lowest Northern Hemisphere summer insolation of the past 800,000 years—387 W/m² at 65°N, caused by a rare combination of minimum obliquity and maximum precession [8]. This extreme orbital configuration drove rapid, hemisphere-wide cooling: atmospheric CO2 fell below 203 ppm, North Atlantic sea surface temperatures dropped to near-glacial levels, and European pollen records indicate glacial-like conditions.

MIS 7c through MIS 6 (−221k to −137k)

The lower left panel spans the later substages of MIS 7 and the transition into MIS 6 (the Saalian glaciation). This period follows the anomalous MIS 7d interval and shows more complex behavior. Interestingly, the cooling branch (b) from −209k to −137k does appear to show horizontal excursions where CO2 holds relatively constant while temperature declines—but at a cooler baseline than observed in MIS 9 or MIS 5e. The entire pattern is shifted to temperatures around −1°C to −7°C and CO2 levels in the 230–260 ppm range.

One possible interpretation is that after the extreme MIS 7d cooling event ended, the carbon cycle resumed its normal asymmetric behavior, but operating within the cooler, lower-CO2 regime characteristic of MIS 7. The "failed interglacial" nature of MIS 7—which never achieved the warmth of MIS 5e or MIS 9 — may have compressed the entire CO2-temperature relationship into a narrower range while preserving the underlying asymmetric dynamics.

This interpretation remains conjectural. The complexity of the MIS 7 substage oscillations makes it difficult to isolate clean warming and cooling branches comparable to those in MIS 9 and MIS 5e.

MIS 5e through MIS 1 (−137k to +1850)

The lower right panel presents the most recent and best-resolved portion of the record, including the Eemian (MIS 5e) and the current Holocene (MIS 1). The warming phase (a) from −137k to −127k shows the sharp rise into the Eemian peak, with COâ?? and temperature rising together along a steep trajectory.

The subsequent cooling branch (b) from −127k to −110k displays the most dramatic example of the phenomenon under investigation. Temperature falls approximately 7°C (from +4°C to −3°C) while CO2 barely changes from its peak near 280 ppm. This is the 13,000- year Eemian plateau—the most striking demonstration in the ice core record that CO2 does not immediately follow temperature during glacial cooling.

The clarity of this signal likely reflects both the magnitude of the Eemian warmth (creating a large CO2 "excess" requiring reabsorption) and the superior temporal resolution of the more recent ice core data. The Eemian thus serves as the type example of asymmetric carbon cycle response, with MIS 9 providing independent confirmation that this behavior is a recurring feature of glacial cycles.

Figure 5: Here the Plot from the Previous Figure is Split into Four Separate Plots. Upper Left Shows the MIS 9 Rise in Temperature (a) and Subsequent fall in Temperature (b). Upper Right Shows the MIS 7e Rise in Temperature (a) followed by fall in Temperature through MIS 7d (b). Lower Left Shows MIS 7c through MIS 6 with Rising Temperature (a) and Falling and Rising Temperature (b). Lower Right Shows MIS 5e through MIS 1 with Rising Temperature of MIS 5e (a) and Falling and Rising Temperature following (b)

Discussion

The analysis reveals that the Eemian CO2 plateau, while the most dramatic example, represents an extreme case of a general pattern observable across multiple glacial cycles. In each transition from interglacial warmth to glacial cold, CO2 shows a delayed response, remaining elevated while temperature has already begun its decline. The exception during MIS 7d, if correctly interpreted as a consequence of extreme orbital forcing, actually reinforces the underlying mechanism by demonstrating that sufficiently rapid cooling can overcome the rate limitation.

Mechanisms for Asymmetric Response

Several mechanisms may contribute to the observed asymmetry between CO2 release and absorption:

• Atmospheric transport limitation: When the ocean warms, dissolved CO2 is immediately available at the air-sea interface for release. When the ocean cools, CO2that was released during warming is dispersed throughout the atmosphere and must find its way back to the ocean surface through atmospheric mixing—a rate-limited process.

• Sea ice expansion: Growing sea ice during cooling physically caps the ocean surface, reducing the area available for gas exchange for the coldest ocean. The Southern Ocean and North Atlantic, major regions of deep water formation and CO2 exchange, become progressively ice-covered during glacial descent.

• Plant stress feedback: As CO2 begins to decline during cooling, vegetation—particularly C3 plants that comprise most forests—experiences carbon limitation. This stress reduces photosynthetic uptake and may cause vegetation die- off, releasing stored carbon and creating a negative feedback that resists CO2decline.

Implications for Causality

The consistent observation of temperature declining while CO2 remains constant poses a challenge for hypotheses in which CO2 is the primary driver of temperature change. If CO2were driving temperature, a 7°C temperature decline while CO2remained constant would require invoking other forcings of comparable magnitude operating in the opposite direction to CO2—forcings that would need to precisely cancel the expected warming effect of elevated CO2 for 13,000 years.

The more parsimonious interpretation is that temperature, driven primarily by orbital (Milankovitch) forcing, leads the system, with CO2 responding as a feedback through ocean degassing, circulation changes, and biosphere responses [9,10]. The asymmetric response times—fast release, slow uptake—create the observed plateaus during the cooling phases of glacial cycles. Climate models continue to have difficulty reproducing glacial inception, suggesting that current understanding of the feedback mechanisms governing interglacial-to-glacial transitions remains incomplete [11,12].

Conclusions

• The Eemian CO₂ plateau—13,000 years of constant CO₂ while temperature fell 7°C—is the most dramatic but not unique example of delayed CO₂ response during glacial cooling.

• Phase plots of CO₂ versus temperature reveal systematic asymmetry: rapid CO₂ rise during warming, delayed CO₂ decline during cooling.

• MIS 9 shows comparable behavior to MIS 5e, confirming this is a recurring feature of glacial cycles spanning at least 350,000 years.

• The MIS 7d interval may represent an instructive exception: extreme orbital forcing (the lowest NH summer insolation in 800,000 years) potentially overwhelmed the rate-limiting mechanisms, allowing CO₂ to track temperature downward.

• The pattern is consistent with rate-limited processes governing CO₂ absorption: atmospheric transport, sea ice formation, and plant stress feedbacks.

• These observations support the interpretation that orbital forcing drives temperature, with CO₂ acting as a delayed feedback rather than the primary driver of glacial-interglacial temperature variations.

Acknowledgement

I would like to thank Renee Hannon for directing my interest toward this subject.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflicts of interest.

References

1. Higginbotham J. Planetary Orbits and Sea-Level. J BiomedRes Environ Sci. 2025 Apr 18; 6(4): 328-339.
2. Higginbotham, J. (2025). Antarctic Ice Core HarmonicAnalysis.
3. Jouzel, J., Masson-Delmotte, V., Cattani, O., Dreyfus, G., Falourd, S., Hoffmann, G., ... & Wolff, E. W. (2007). Orbital and millennial Antarctic climate variability over the past 800,000 years. science, 317(5839), 793-796.
4. Parrenin, F., Barnola, J. M., Beer, J., Blunier, T., Castellano, E., Chappellaz, J., ... & Wolff, E. (2007). The EDC3 chronology for the EPICA Dome C ice core. Climate of the Past, 3(3), 485-497.
5. Bereiter, B., Eggleston, S., Schmitt, J., Nehrbass�Ahles, C.,Stocker, T. F., Fischer, H., ... & Chappellaz, J. (2015). Revision of the EPICA Dome C CO2 record from 800 to 600 kyr before present. Geophysical Research Letters, 42(2), 542-549.
6. Milankovitch, M. (1941). Kanon der Erdbestrahlungen und seine Anwendung auf das Eiszeitenproblem. Royal Serbian Academy. Journal of Natural Sciences and Mathematics, 33.
7. Higginbotham, J. (2026). Climate, Water Vapor, and Volcanic Eruptions. Env Sci Climate Res, 4(1), 01-09.
8. Moseley, G. E., et al. (2021). Precise timing of MIS 7 substages from the Austrian Alps. Climate of the Past, 17, 1443–1454.
9. Humlum, O., Stordahl, K., & Solheim, J. E. (2013). The phase relation between atmospheric carbon dioxide and global temperature. Global and Planetary Change, 100, 51-69.
10. Grabyan, R. (2025). Global Atmospheric COâ?? Lags Temperature: Implications for Short-and Long-Term Climate Predictions. Science of Climate Change, 5(3), 302-326.
11. Vavrus, S. J., He, F., Kutzbach, J. E., Ruddiman, W. F., & Tzedakis, P. C. (2018). Glacial inception in Marine Isotope Stage 19: An orbital analog for a natural Holocene climate. Scientific reports, 8(1), 10213.
12. Willeit, M., Calov, R., Talento, S., Greve, R., Bernales, J., Klemann, V., ... & Ganopolski, A. (2024). Glacial inception through rapid ice area increase driven by albedo and vegetation feedbacks. Climate of the Past, 20(3), 597-623.