Variance-driven decomposition based on the singular spectrum analysis of the European Project for Ice Coring in Antarctica (EPICA) δD, CO2, and CH4 records allowed a novel quantitative structural interpretation of all glacial/interglacial cycles and glacial terminations of the last 800 kyr. This bottom-up approach used the response components of EPICA stacked records to reconstruct the envelope of the thermal response through a physical interference model. The aim was to improve understanding of the intensity, amplitude, and asymmetry features of 73 marine isotope stages/substages (MISs) and seven glacial terminations. The Antarctic stack record can be described by a variance-weighted superposition of ten thermal waves of different origins (mid-term oscillation, orbitals, and suborbitals) that stochastically interfere at a given time according to their relative differences in frequency, amplitude, and polarity. Interglacial/glacial stages resulted from constructive interference and bipolar amplification of warming/cooling responses, respectively. The low-intensity MISs (including 90% of substages) and the unbiased-dated terminations fell in the low-interference regions, where dominant destructive patterns minimize the thermal envelope. The positive skewness of the EPICA stack resulted from constructive interference with a strong bias in the warming direction, especially after the Mid-Brunhes Event. Duration analysis of short eccentricity hemicycles exhibited an intrinsic unexpectedly prolonged mean cooling in the nominal solution (5.8 kyr) and its EPICA response as well (8.6 kyr), along with an interference-induced asymmetry (21.1 kyr). The overall effect has led to the saw-tooth shape of glacial cycles, which was strongly induced by interference.
P. Viaggi (2026) studied this question.