Temperature and mineral dust records serve as valuable palaeoclimatic indicators for studying atmospheric variability across different temporal scales. In this study, we employed Haar fluctuations to analyse global spatiotemporal atmospheric variability over the Last Glacial Cycle, capturing both high- and low-frequency information within the records, regardless of uniform or non-uniform sampling. Furthermore, we utilised Haar fluctuations to compute fluctuation correlations, thereby enhancing our understanding of palaeoclimate dynamics.
Our findings reveal a latitudinal dependency in the transition from macroweather to climate regimes (τc), with polar regions experiencing shorter transitions compared to the tropics and mid-latitudes. These transitions occur at approximately 1/100th of glacial cycle length scales, suggesting a dominant forcing mechanism beyond Milankovitch cycles. Additionally, our analysis shows that polar regions have larger fluctuation amplitudes than lower latitudes as a consequence of the polar amplification effect. Furthermore, fluctuation correlations demonstrate faster synchronisation between the poles themselves compared to lower-latitude sites, achieving high correlation values within 10 kyr.
Therefore, our findings suggest a consistent climate signal propagating from the poles to the Equator, representing the first empirical evidence supporting the hypothesis that the poles play a pivotal role as climate change drivers, influencing the variability in climatic transitions worldwide.
Fondo Nacional de Desarrollo Científico y Tecnológico12316821221220Introduction
The study of palaeoclimate, which involves reconstructing past climate variations using various palaeoclimatic archives, provides crucial insights into understanding Earth's climatic history . Palaeoclimatic archives, such as ice cores and sedimentary records, offer unique opportunities to investigate climate dynamics on various timescales, ranging from annual to millennial and beyond . These archives preserve information about past changes in relevant climate variables such as air and sea temperatures, atmospheric aerosol content, and precipitation regimes. This enables us to discern the patterns and drivers of climate variability over extended periods .
Understanding the temporal complexity of palaeoclimatic records has been a long-standing challenge in climate science. The inherent non-linear and chaotic nature of Earth's climate system often leads to multifaceted interactions and feedback mechanisms acting over a wide range of scales. Conventional statistical methods, although essential, may not fully capture the intricate dynamics within the palaeoclimatic time series, leaving significant patterns undiscovered. Furthermore, because of the many different physical processes interacting within the climate system, its variability spectrum is far from being composed of pure spikes belonging to well-defined forcing processes, such as daily variations or Milankovitch cycles. Instead, the bulk of the spectrum is continuous and scales over a wide range. Indeed, from fractions of a second to hundreds of millions of years, almost all of the variability resides in wide scaling regimes that for a long time had been reduced to an uninteresting “background” . This extends by many orders of magnitude the range of scales considered by , who showed that, over the range of millennia to Milankovitch scales, quasi-periodic signals represent only a small fraction of total variability.
The quantitative importance of the continuum part of the variability is slowly being recognised. A recent ad hoc “solution” proposed by is to redefine the background as a “1/f” noise instead of white noise. This redefinition allowed them to produce a “conceptual landscape” displaying interesting narrow-scale range processes above the new background level. Because the range of scales shown in the conceptual landscape was 1013, the redefinition effectively boosts the background of the high frequencies by the same factor. Although this (partially) offsets the factor of 1015 error in the “educated guess” of , as shown in , the resulting boosted background is not yet representative of the empirical data. Indeed, the ad hoc assumption of a 1/f background constitutes an extreme scaling hypothesis that only approximates empirical measurements. As argued in , one way to improve our analysis of the background spectrum is to investigate real-space fluctuations (instead of Fourier-space spectra), such as the Haar fluctuation.
Classically, atmospheric dynamics is divided into two basic regimes: weather and climate. Generally, the definition of the climate regime is the 30-year average proclaimed by the International Meteorological Organization in 1934, despite abundant data indicating long-term trends, particularly in the context of modern climate change. In a non-steady state, dynamical regimes should be objectively defined, and scaling is one way to achieve this. If we denote by ΔT(Δt) a fluctuation in the temperature over a time interval Δt, then in a scaling regime ΔT(Δt)≈Δtζ. In a general (multifractal) case, ζ is a random exponent. To simplify the discussion, we average this relationship to obtain the average fluctuation 〈ΔT(Δt)〉≈ΔtH, where H is the (non-random) “fluctuation exponent” (“〈〉” denotes averages). If such a relation holds, then, to within the constant factor λH, fluctuations over a long time interval λΔt, (λ>1) are the same as those over a shorter timescale Δt: 〈ΔT(λΔt)〉≈λH〈ΔT(Δt)〉. Through systematic scaling analyses of atmospheric fluctuations over the range of milliseconds to millions of years, and suggested that atmospheric dynamics has five distinct regimes rather than two. They are listed below in order of increasing timescale.
The weather regime: this extends from turbulent dissipation scales (below 1 s) to the lifetime of planetary structures (typically about 10 d in the atmosphere but about 6 months in the ocean). In this regime (after removing the diurnal cycle), average fluctuations 〈ΔT(Δt)〉 have H>0 and hence increase in amplitude with timescale, revealing “unstable” temperature time series that tend to “wander” up and down.
The macroweather regime: beyond the weather regime lies the macroweather regime, in which fluctuations tend to decrease with scale (H<0). In this regime, successive fluctuations tend to cancel, and the temperature anomalies tend to converge to well-defined values. This is the “long-term weather” regime that resembles the classical 30-year average “climate”: by increasing the timescale, the average converges more and more.
The climate regime: the macroweather convergence of anomalies does not extend beyond a critical timescale τc, beyond which the apparent point of convergence starts to wander (again H>0). In the Anthropocene, the situation is fairly straightforward: τc is the scale at which the temperature response to anthropogenic forcing exceeds the internal variability; τc is currently about 15 years. The variability in temperature fluctuations in the Anthropocene clearly shows the importance of τc, a transition between macroweather convergence to a well-defined “climate state” and a new lower-frequency regime of “climate change”. The example of the Anthropocene also provides insight into the causes of the transition: the emergence of new sources of variability, in this case primarily greenhouse gases. While general circulation models (GCMs) can reasonably reproduce the current τc when using anthropogenic forcing, in control runs (fixed forcing), they show macroweather convergence to a well-defined climate state out to their longest simulated timescales . In the pre-industrial epoch, neither the transition time τc nor its origin is known. For example, high-resolution dust and temperature data from Antarctic ice cores suggest that it is of the order of several centuries to several millennia, varying from one interglacial to another . However, these timescales are too short to be consequences of Milankovitch forcing. Is there a new source of external forcing, or is there a new non-linear internal forcing mechanism? There is also reason to expect geographical variation in τc, which is the object of this study.
The macroclimate regime: the wandering (H>0) climate regime extends to Milankovitch scales (≈100 kyr) beyond which, again, H<0 and fluctuations decrease with scale and anomalies converge. Little is known about this regime, although the negative H regime is thought to be a feature of palaeoclimatic data between timescales of about 100 kyr and 1 Myr throughout the Phanerozoic.
The megaclimate regime: this H>0 regime extends from ≈1 Myr to at least several hundred million years and seems to be fundamentally bio-geological. Wandering palaeotemperatures with H>0 apparently coexist with converging (H<0) extinction and origination rates over the entire range. Interestingly, fundamentally important biodiversity seems to be dominated up to about 40 Myr by the destabilising climate but beyond that by stability macrobiology .
Over the past decade, numerous studies have focused on scaling palaeoclimatic data . Some studies have evaluated how well climate models reproduce the observed long-term climate variability e.g.. While climate models and reconstructions show reasonable consistency on seasonal to inter-annual timescales, significant discrepancies emerge on longer timescales and at regional scales, with models underestimating the variability . The best-known variability at millennial timescales is related to Dansgaard–Oeschger (DO) events, which are abrupt temperature changes originating in the North Atlantic region. They influence the climate over the entire Northern Hemisphere and are phase-lagged compared to their Southern Hemisphere counterparts, the Antarctic Isotopic Maximum (AIM) events . They are mechanistically linked to variations in the Atlantic Meridional Ocean Circulation . At multi-millennial timescales, the best-known variability is caused by external orbital forcing, which redistributes incoming solar radiation around the globe. The main periodicities are linked to the precession (∼19–21 kyr), obliquity (∼41 kyr), and eccentricity (∼100 and 400 kyr) of the Earth's orbit, also called the Milankovitch cycles. The weakness of models in capturing low-frequency variability means that, although weather (and to a certain degree macroweather) variability can be analysed in simulation data, the longer climatic timescales and the transition between macroweather and climate can only be reliably investigated in empirical records. However, because of the limited availability of palaeoclimatic records with high resolution and long duration, such studies typically focus on timescales shorter than glacial cycles , and longer spatiotemporal patterns are not yet well understood. As a result, little is known about the geographic differences in τc (the transition timescale between macroweather and climate).
In this study, we compiled a list of palaeoclimatic temperature and dust records that were measured at high resolution and span most of the last glacial–interglacial cycle. Temperature data are typically based on isotope measurements , whereas dust data are typically based on direct particle measurements or elemental data (e.g. calcium; ). Dust particles are usually entrained in continental arid regions and then transported by prevailing winds towards remote regions of the globe, where they are deposited by gravitational settling or precipitation over the surfaces of oceans, lakes, ice sheets, etc. . Unlike isotopes, dust particles measured in palaeoclimatic archives are therefore not representative of a single atmospheric variable but indicate a more holistic state of the atmosphere (hydrological cycle, atmospheric dynamics). Palaeoclimatic proxy records usually have non-equidistant sampling rates and sometimes lose resolution as we move to older time sections, limiting high-frequency analyses. Preprocessing techniques such as linear interpolation can convert non-uniform data resolution to a uniform time series, but this always implies some high-frequency information loss, which is essential for finding scaling transitions. In contrast, Haar fluctuations are convenient for quantifying variability as a function of timescale. This enables us to accurately estimate fluctuations and other scaling exponents , even for data with highly non-uniform chronologies . In this context, Haar fluctuations have been shown to be more robust and less biased than conventional spectral analysis .
Here, we analyse the global spatiotemporal atmospheric variability over the Last Glacial Cycle, using a variety of palaeoclimatic temperature and dust records. We calculated Haar scaling for each site and estimated the macroweather-to-climate transition scale τc for various latitudinal ranges. Because climate involves many intimately coupled processes, understanding variability at any timescale also requires spatiotemporal analysis . Hence, we also used our non-uniform Haar fluctuation to quantify the fluctuation correlations between different datasets as a function of the timescale.
Dataset description
For our analysis, we used palaeoclimatic records of both temperature and dust. Records obtained from marine sediment cores, ice cores, loess deposits, and lake sediments were selected. Most sites have records of either temperature or dust, with only a few sites having both, which may have been measured at different resolutions. Because we are specifically interested in covering the macroweather and climate regimes, our analysis mostly includes palaeoclimatic archives that cover at least the past 50 000 years at centennial resolution or better. Figure and Table show the locations and characteristics of each record included in this study.
The global distribution of the datasets included in this study. The temperature and dust flux records come from ice cores, loess, lake sediments, and marine sediments.
Most sites are located in polar regions, such as Greenland and Antarctica, where measuring proxies at a high temporal resolution is easier to achieve than in the tropics and mid-latitudes. Unfortunately, no datasets that satisfied our selection criteria were available for regions such as Central and South America, central Asia, Australia, the South Atlantic Ocean, and the southern Indian Ocean. Some datasets include multiple glacial cycles, such as the EPICA Dome C records or the Xifeng loess archive, whereas others cover only a single glacial–interglacial cycle or less. To avoid mixing information from multiple glacial cycles, we restricted the analysis to the most recent 130 kyr, roughly the last glacial–interglacial cycle.
Palaeoclimate archives rarely have uniform resolution because each data point corresponds to a thin layer of core material, and each layer typically relates to a different time span depending on depth. Figure shows the cumulative age covered by the sampling for each dataset. Note that a uniform dataset has a linear relationship, but most curves are above the diagonal, indicating a higher resolution in the early samples and a lower resolution in the older samples. We observe that the sampling rates depend on a specific dataset, and most are significantly non-uniform at all scales; that is, the non-uniformity itself is likely to be scaling.
Cumulative age versus cumulative sampling, depicting the proportion of age covered by the proportion of sampling in the dataset, corresponding to the most recent 130 kyr of the datasets. The continuous lines depict temperature datasets, while the dotted lines depict dust datasets.
The datasets contain records associated with various dating scales during the Last Glacial Period. We aimed to use the AICC2012 or GICC05 dating scale on as many records as possible. Table lists the specific dating scales used in this study. In Fig. , we show all the time series over the most recent 130 kyr.
Temperature (a) and dust (b) time series comparison for the last 130 kyr. See Table for the IDs of the records.
Methodology
This study aimed to analyse climate scaling regimes from palaeoclimate data, which tend to have non-equidistant samples and many spikes. Non-uniform spectral estimators, such as Lomb–Scargle, have serious problems with uncontrolled spectral leakage, leading to low accuracy when applied to non-uniformly sampled scaling data. We applied a non-uniform FFT algorithm to the data, but the scaling regimes were unclear from the analysis. Haar fluctuations are more robust on spiky data and better suited for identifying different scaling regimes. Furthermore, no preprocessing of the data, such as normalisation and outlier filtering, is necessary for Haar fluctuation analysis .
Haar fluctuations
Haar fluctuation analysis is a powerful technique for characterising the timescale-dependent variability in a time series F(t). The Haar fluctuation has a simple definition, and the results have a clear interpretation. Furthermore, it can be extended to irregular sampling rates, thus avoiding data preprocessing such as interpolation between time samples. Hence, the algorithm uses all the information available from the data, including shorter timescales represented by high-resolution sections.
The Haar wavelet is a piecewise constant function that is used as a shape function in fluctuation analysis. The Haar fluctuation ΔF over a uniform time interval Δt is defined as the absolute difference in the mean over the first and second halves of an interval:
ΔF(Δt)=1Δt2∫t-Δtt-Δt/2F(t′)dt′-∫t-Δt/2tF(t′)dt′.
When the Haar fluctuation over a range of timescales Δt obeys a power-law relation in the form of ΔF(Δt)∝ΔtH, the mean absolute fluctuation also varies as
〈|ΔF(Δt)|〉∝ΔtH,
where the brackets “〈〉” indicate the sample mean over all available disjointed intervals of size Δt (fluctuation period) in the time series. Here, the exponent H characterises the scaling regime, with H>0 indicating that the average fluctuations increase with scale, while H<0 implies the opposite. This fluctuation exponent estimation is appropriate within the range H∈[-1,1], which is valid for almost all geophysical data analysed to date . A more general qth moment of the fluctuations can be estimated by defining the generalised structure functions Sq(Δt) as
Sq(Δt)=〈|ΔF(Δt)|q〉∝Δtξ(q),
where the exponent
ξ(q)=qH-K(q)
includes a linear term and a non-linear function K(q) that characterises the multifractal scaling properties of the series F(t). This includes both the intermittency and extremes.
The definition of Haar fluctuation in Eq. () assumes a uniform time step. The extension of Haar analysis to non-uniform time series is detailed in Appendix , following a similar approach as demonstrated in Appendix B.
Fluctuation correlation
Thus far, climate variability has been described by scaling relationships in the form of ΔF(Δt)∝ΔtH and calculated using Haar fluctuations for a single time series only. If the governing processes have consistent scaling patterns, we should expect the interrelations between the processes to lead to correlations. The correlations between two time series at different scaling regimes can be calculated using fluctuation correlation analysis. This technique estimates the normalised correlation coefficient (ρ(Δt)) between individual fluctuations at fixed timescales; that is, there is one correlation value at each Δt.
Let us consider two time series denoted as A(t) and B(t). The Haar fluctuations have the property that they are independent of any additive constant; therefore, 〈ΔA〉=〈ΔB〉=0. Furthermore, ΔA(Δt) and ΔB(Δt) denote the Haar fluctuations in A and B, respectively. Then, at a time lag Δt, we can compute the fluctuation correlation as follows:
ρΔA,ΔB(Δt)=〈ΔAΔB〉ΔA21/2ΔB21/2,
with the correlation coefficient ρΔA,ΔB(Δt)∈[-1,1].
Because the time series likely have different sampling rates and durations, their fluctuations will cover timescales that may not match each other, leading to an inaccurate estimation of ρ. To avoid this problem, we use linear interpolation to shift the sampling of the higher-resolution series, say A(t), and fit it with the sampling of the lower-resolution series, say B(t), before computing its fluctuations. Shortening the time series duration to the maximum common Δt is also required. Finally, we fitted a cubic spline to simplify the visualisation because we ended up with a cloud of fluctuations.
Results
This section presents the results of the Haar fluctuation analysis conducted on the selected datasets. Section illustrates the scaling regimes observed at individual sites and provides a comparative analysis of the Haar fluctuation amplitudes. Fluctuation correlations are presented in Sect. .
Haar fluctuation analysis
The Haar analysis yields a fluctuation root mean square (RMS) at each time lag so that the scaling exponent H can be fitted for different regimes, as given by Eq. (). To analyse the spatial differences in the scaling regimes, we clustered the measurement sites into five different groups according to latitude:
Arctic (from 90 to 60°)
northern mid-latitudes (from 60 to 20°)
the tropics (from 20 to -20°)
southern mid-latitudes (from -20 to -60°)
Antarctica (from -60 to -90°).
As explained in Sect. , when the average fluctuations increase with scale (H>0), they represent unstable behaviour, whereas, when they decrease with scale (H<0), they converge towards a mean state. To identify the various scaling regimes, we visually detected different sections with scaling behaviour and fitted a linear segment to estimate the slope (H). These exponents are shown as the slopes of the RMS values from the Haar fluctuation in a log–log plot. Figures to present the Haar fluctuation RMS values for the aforementioned five latitude groups.
Haar fluctuation root-mean-square (RMS) values derived from the Arctic records at various scales, presented in log–log format. Bold letters along each curve denote the corresponding site names, identified by the data IDs in Table . Units for each case vary and can be found in Table . Temperature proxies are displayed as continuous lines, and dust proxies are displayed as dotted lines. Thick dashed lines in black (dust) and red (temperature) are added as slope references, together with the estimation of the H exponent. An arrow labelled with τc points to the estimated climate regime shift (fluctuation minimum). Vertical dashed grey lines indicate the multi-millennial scale range where fluctuations peak.
Haar RMS values derived from the northern mid-latitude records. Temperature proxies are displayed as continuous lines, while dust proxies are displayed as dotted lines. Refer to Fig. for detailed annotations and references.
Haar RMS values derived from the tropical records. Temperature proxies are displayed as continuous lines, while dust proxies are displayed as dotted lines. Refer to Fig. for detailed annotations and references.
Haar RMS values derived from the southern mid-latitude records. Temperature proxies are displayed as continuous lines, while dust proxies are displayed as dotted lines. Refer to Fig. for detailed annotations and references.
Haar RMS values derived from the Antarctic records. Temperature proxies are displayed as continuous lines, while dust proxies are displayed as dotted lines. Refer to Fig. for detailed annotations and references.
The Arctic Haar RMS fluctuations are shown in Fig. . The finer resolution of the dust and temperature NEEM datasets allows us to identify the presence of an H<0 regime (macroweather) and an H>0 regime (climate), with a fluctuation minimum suggesting a climate regime change (τc) between Δt=100 and 300 years. Although the dust GRIP and NGRIP results also show a minimum at higher frequencies, they have a slope greater than 1, falling beyond the accurate range of the Haar analysis. Nevertheless, given the positive value of H at multi-millennial scales and acknowledging that H<0 at sufficiently high frequencies (the macroweather regime), we can deduce that τc is simply below the resolution covered. A similar scenario holds for the RECAP, GRIP, and NGRIP temperature records. In addition, all sites exhibit higher fluctuations towards larger timescales (i.e. H>0), peaking around 17 and 55 kyr. As shown in the following figures, low-frequency peaks can be observed across all latitudes.
The northern mid-latitude RMS fluctuations are shown in Fig. . The terrestrial record from Xifeng (dotted blue) features a fluctuation minimum, suggesting a climate regime transition around 2 kyr. A low-frequency peak appears at timescales ranging between 12 and 30 kyr. In the marine temperature records, a wide millennial-scale peak can also be identified at the Portugal site (Δt≈3 kyr), and a low-frequency peak can be identified in the North Atlantic Ocean (Δt≈30 kyr). In contrast, multi-millennial-scale fluctuation amplitudes appear to plateau at the central California site. None of the marine temperature sites exhibit a clear minimum in their fluctuations. Estimating H for these records is not feasible because of the presence of oscillations at centennial to multi-millennial scales.
Figure shows the RMS fluctuations from the tropics. Lake Tanganyika (LT) exhibits a fluctuation minimum (τc) at timescales of approximately 600–1000 years and a low-frequency peak between 15 and 40 kyr. Conversely, the Indian Ocean features a τc close to 10 kyr, but a clear low-frequency peak cannot be identified. In contrast, marine dust records do not indicate an identifiable climate regime transition. Only 31BB (dotted orange line) shows a symmetry break at approximately Δt=500 years; nonetheless, as its slope on short timescales is less than the minimum detectable by Haar (H<-1), drawing any meaningful conclusions is not possible. However, because H>0 at multi-millennial scales, we may suppose that τc is below the resolution of the data.
The PS75 marine record is the sole dataset located in the southern mid-latitudes. Despite oscillations at higher frequencies, Fig. shows a minimum (τc) in its Haar fluctuation RMS within the timescales of 1.5 to 3 kyr. Beyond this range, the fluctuations increase with scale. A low-frequency peak is not detectable, which is consistent with the results from marine records in the northern mid-latitudes and the tropics.
Absolute amplitude differences in the fluctuations for calibrated temperature datasets in degrees Celsius. Shown are (a) polar regions, (b) northern mid-latitudes, and the (c) tropics. The numbers accompanying the vertical lines estimate the fluctuation amplitude, equivalent to the difference between the peak found at multi-millennial scales minus the minimum found at multi-centennial or millennial scales. Some τc values are shown as reference according to the findings from Figs. –.
Climate regime transition summary. “High-frequency slope” denotes the H exponent found at decadal to centennial timescales (e.g. H14 or H16 in Fig. ), “low-frequency slope” denotes the H exponent found at centennial to millennial timescales (e.g. H15 or H17 in Fig. ), and τc corresponds to the transition timescale from high-frequency slope to low-frequency slope. The LF peak column shows the estimated range in kyr where a fluctuation maximum is attained, whereas the Amplitude column refers to the difference between maximum and minimum fluctuations.
Finally, Fig. shows the Haar fluctuation RMS from the Antarctic records. We observe consistent results with similar scaling patterns across different sites and between dust and temperature proxies. The slopes of the reference lines indicate a minimum fluctuation between timescales of 200 and 750 years in all cases, showing that the macroweather-to-climate transition (our estimate of τc) occurs at multi-centennial scales in the Antarctic. Furthermore, all datasets exhibit a low-frequency peak within the 15 to 45 kyr range.
Figure illustrates the comparison of absolute amplitude differences among the available temperature data, measured in degrees Celsius. In this amplitude comparison, we exclusively utilise temperature data measured in degrees Celsius, thereby avoiding the need for calibration conversions and mitigating potential biases in the analysis. Note that the amplitude calculated here is approximately the total range of fluctuations. When comparing panels a–c, it is evident that the amplitudes gradually decrease from the poles towards lower latitudes, going from 4–6 °C at the poles to 1–2 °C in the tropics, which is in line with the concept of polar amplification . However, it is worth noting that there are no data from the southern mid-latitudes in this analysis. Table summarises the Haar statistics for the different latitudes.
Fluctuation correlation of temperature archives. The light-blue, red, and yellow areas highlight the pole–pole, pole–terrestrial, and pole–marine record correlations, respectively. Both bands cover 1 standard deviation from the mean, in both directions. Panel (a) displays the correlations between series up to Δt=100 kyr, while panel (b) zooms in on only the first 20 kyr of Δt.
Fluctuation correlation analysis
This section aims to quantify the interrelation between the poles and lower latitudes. We can estimate the pairwise time-dependent correlation coefficients using Eq. (). We perform this analysis separately on the temperature and dust records.
Figure shows the fluctuation correlations of the temperature records. Here, we estimate 15 pole–pole correlations (i.e. the fluctuation correlations between Antarctic and Arctic records) and 40 pole–lower-latitude correlations. Lake Tanganyika, the only terrestrial temperature record, is depicted independently of the other marine records. Thus, the light blue, red, and yellow areas highlight the pole–pole, pole–terrestrial, and pole–marine record correlations, respectively. Shaded areas cover 1 standard deviation from the mean, truncated by the upper limit of ρ=1.
Figure shows that all records exhibit lower fluctuation correlation values at high frequencies, which increase with scale. In particular, the poles attain ρ=0.75 at Δt=10 kyr (Fig. b) and remain highly correlated across all scales. On the other hand, the correlation between the poles and lower-latitude sites generally shows lower correlation values with higher standard deviation (wider yellow and red areas), indicating less synchronicity between themselves. It was unexpected that Lake Tanganyika (LT), despite being located in the tropics, correlates as rapidly as the poles, reaching ρ=0.75 at similar scales. However, it shows a sharp drop towards Δt=50 kyr, i.e. at scales of the length of the LT time series.
Similarly, Fig. shows the fluctuation correlations of the dust records. Here, we estimate 8 pole–pole correlations and 30 pole–lower-latitude correlations. Xifeng, the only terrestrial dust record, is depicted independently of the other marine records. Thus, the light-blue, red, and yellow areas highlight the pole–pole, pole–terrestrial, and pole–marine record correlations, respectively.
Fluctuation correlation of dust archives. The light-blue, red, and yellow areas highlight the pole–pole, pole–terrestrial, and pole–marine record correlations, respectively. Both bands cover 1 standard deviation from the mean, in both directions. Panel (a) displays the correlations between series up to Δt=100 kyr, while panel (b) zooms in on only the first 20 kyr of Δt.
In Fig. , we observe a consistent pattern with our previous findings: all records exhibit lower fluctuation correlations at high frequencies, gradually increasing to a nearly perfect positive correlation (ρ=1) around scales of 100 kyr. Moreover, the correlation of polar records rises quickly with scale, reaching ρ=0.75 at Δt≈10 kyr, and they remain strongly correlated with ρ>0.75 over longer timescales. On the other hand, the pole–lower-latitude relations also display stronger fluctuation correlations as Δt increases but at a slower pace (see Fig. b), and they never surpass the pole–pole correlation values.
Marine and terrestrial records seem to correlate similarly with the poles, both reaching a maximum of ρ≈0.8 at approximately Δt=35 kyr then dropping rapidly from Δt=40 to 75 kyr. This drop is visible for both the pole–pole and pole–lower-latitude correlations but is much more pronounced in the latter case, even reaching negative correlations for certain pairs of sites.
Discussion
In this study, we analysed global spatiotemporal atmospheric variability across annual to multi-millennial timescales using Haar fluctuations. Our results suggest that the transition timescale (τc) from macroweather to climate is not consistent globally but varies systematically according to latitude. While τc occurs at centennial to multi-centennial timescales in polar areas (Δt=100–300 years in the Arctic; Δt=200–750 years in Antarctica), it occurs at longer multi-centennial to millennial timescales at lower latitudes (Δt=600–10 000 years in the tropics; Δt=1500–3000 years in the mid-latitudes). Furthermore, almost all our τc estimations are around 1/100th of glacial cycle length scales, which is too short to be attributed to Milankovitch forcing, raising the question of what processes could be behind this transition.
Antarctica was the only group in which consistent results across all sites and variables were observed (Fig. ), probably because of the significant influence of similar atmospheric processes on Antarctic temperature and dust . analysed high-resolution ice core data from the EDC dust series, estimating τc for eight glacial–interglacial phases, averaging approximately 500 years, which is consistent with our Antarctic τc estimation of 200–750 years. For the other latitude groups (Figs. –), the estimation of τc relied on data from only one or two sites, which introduced some uncertainty into the results. However, the overall consistency observed across different latitudes helps mitigate this uncertainty.
Fluctuations from terrestrial records (Figs. and ) display a distinct low-frequency peak that appears to be centred around the precession and obliquity cycles. This peak is most probably related to variability in the East Asian and African monsoons but is absent in marine data, thus highlighting the impact of proxy-specific climate recording processes and underscoring why single-proxy analysis over mixed-proxy approaches should be preferred . The diverse shapes of these peaks across records may be attributed to the varied influence of precession, obliquity, and eccentricity at different locations and to the limitations of the Haar fluctuation analysis, which tends to produce broad peaks rather than clearly identifying specific periodicities.
Another distinct characteristic of terrestrial records is that they do not match the τc calculated in marine records for their respective latitudinal zones. For example, the dust record from Xifeng in China (Fig. ) indicates a transition-scale τc around 2 ka. Such a transition is not corroborated by the marine temperature records of the northern mid-latitudes that suggest a transition scale somewhere below the highest resolution of 0.5 ka. In the tropics, the terrestrial temperature record of Lake Tanganyika in eastern Africa suggests a transition-scale τc around 600 years. In contrast, a marine temperature record from the tropical Indian Ocean suggests a transition around 10 ka, and no conclusive transition could be found in marine dust records from the tropical Pacific within their resolution of about 0.5 ka. These discrepancies suggest that, in addition to the large-scale latitudinal dependency of τc discussed below, some regional differences exist between marine and terrestrial environments.
By analysing the absolute amplitude differences in the fluctuations (Fig. ), we observed that regions with shorter τc values displayed larger fluctuation amplitude differences. Higher fluctuation amplitudes in polar areas and lower values in the mid-latitudes and tropics are a consequence of the polar amplification effect in both temperature and dust . In Fig. we hypothesise that this polar amplification effect may be accompanied by a polar acceleration effect, whereby the transition from macroweather to climate occurs at shorter timescales at higher latitudes. One notable exception is the tropical records from Lake Tanganyika which display a very small τc of about 0.6 ka, comparable with polar values. With the currently available records, this hypothesis is weak because (1) it is partly contradicted by the available data, (2) there are only dust records available for the mid-latitude bands, and (3) there are only temperature records available for the tropical band. Still, the apparent coincidence between the two effects raises the interesting question of potential causality (in either direction) between the two. Additional reconstructions of both temperature (most urgently needed in mid-latitudes) and dust (most urgently needed in the tropics) are now needed to test the hypothesis.
Schematic of the latitude–τc–amplitude relation found considering both dust and temperature archives. Grey markers indicate the estimated range of τc values according to latitude, while the background colour map, ranging from light blue to light red, represents the fluctuation amplitudes across latitudinal bands.
Fluctuation correlation analysis revealed that local variability dominates at high frequencies, whereas global synchronisation appears at longer timescales (Figs. and ). In particular, both the temperature and the dust results show that the poles synchronise within 10 kyr and remain highly correlated at longer timescales. The North Atlantic (where Greenland is located and where all our Arctic records are from) and the South Atlantic Ocean (which strongly influences Antarctic temperatures) are mechanistically linked through the bipolar seesaw . In addition to orbital forcing, these two regions are therefore also linked by this mechanism that produces synchronised multi-millennial variability described as Dansgaard–Oeschger events in the north and Antarctic Isotopic Maxima in the south . In contrast, fluctuation correlations between polar regions and lower latitudes show a slower increase towards multi-millennial scales, indicating no significant synchronising mechanism between the poles and the lower latitudes beyond the orbital forcing.
In temperature correlations, the terrestrial record of Lake Tanganyika (LT) unexpectedly exhibits a high and rapid correlation with the poles, despite its equatorial location. We note that this record also exhibits a very short τc similar in value to the polar records. In contrast, the terrestrial dust record from Xifeng, China, synchronises with the poles at a pace similar to that of marine dust records at lower latitudes. It is therefore not clear whether the continental nature of the LT site plays a role in its synchronisation with the poles. We also note that LT's record is the shortest of all those selected in this study (approximately 60 kyr with 209 samples), reducing the confidence in the low-frequency range of this site's analysis but not in the τc calculation.
Therefore, our results suggest that global synchronisation is primarily driven by the poles, as they exhibit shorter τc values and faster correlation among themselves compared to lower latitudes. This is consistent with the current paradigm that ice sheets at high latitudes were instrumental in driving global climate change . Conversely, as we approach the Equator, the fluctuation amplitude diminishes and correlations decrease, resulting in longer transition timescales for fluctuations to participate in the glacial cycle.
Conclusions
In this study, we utilised Haar fluctuation analysis to explore the scaling properties of various palaeotemperature and palaeodust records. Our work underscores the importance of generating palaeoclimatic datasets with both extensive coverage and high temporal resolution. Through the analysis of uniform and non-uniform time series, we confirmed the efficacy of Haar structure functions in assessing the scaling properties of palaeoclimatic data, offering reliable scaling exponent estimates across a broad range of scales. Moreover, the straightforward interpretation, tied to the behaviour of fluctuations with scale changes, facilitates the characterisation of climate transitions.
Our investigation of the τc parameter reveals latitudinal dependencies in the transition from macroweather to climate regimes. We observed centennial to multi-centennial transitions in polar regions versus longer multi-centennial to millennial transitions at lower latitudes, highlighting the intricate interplay between regional processes and global climatic influences. Furthermore, the τc parameter suggests the existence of a dominant forcing mechanism not aligned with Milankovitch's cycles (1/100th of glacial cycle length scales), prompting inquiries into its underlying causes.
Additionally, we found that the amplitude of fluctuations is greater at higher latitudes and decreases towards the tropics. This observation elucidates how τc shifts from shorter to longer scales and may explain why lower latitudes require more time to synchronise with the poles, as the modulated signal propagating from the poles to the Equator weakens. Moreover, this finding may help reconcile some Holocene results, suggesting that the macroweather climate transition occurs at millennial rather than centennial scales, as Holocene data are often collected at mid-latitudes or at the tropics, where longer transition scales are expected e.g..
Fluctuation correlation analysis revealed that poles synchronise faster than lower latitudes, achieving full correlation before half a cycle (50 kyr), indicating synchronous behaviour during glacial periods. This underscores the significant impact of Dansgaard–Oeschger and/or AIM events and orbital forcing in polar regions. In contrast, the co-variability between polar regions and lower latitudes is less influenced by millennial-scale events and is more concentrated around orbital periodicities, with a distinct drop in dust correlation in the band between obliquity and eccentricity.
In summary, the latitudinal dependency in τc and the observed signal propagation from the poles to the Equator underscore the significance of polar regions in influencing global climatic variability. Thus, our findings represent the first empirical evidence supporting the hypothesis of the pivotal role of the poles as climate change drivers.
Future studies can enhance fluctuation estimation robustness by comparing or averaging results across different glacial cycles. Additionally, assessing lagged cross-correlations between poles, mid-latitudes, and tropics may clarify directional causality, albeit without conclusive evidence. Furthermore, addressing the current scarcity of high-resolution time series from mid-latitudes and the tropics is crucial for improving the characterisation of geographic variability in glacial–interglacial cycles at these latitudes.
Haar fluctuations in non-equidistant data
The definition of Haar fluctuation in Eq. () assumes a uniform time step. Here, we follow an approach similar to Appendix B and extend the Haar analysis to non-uniform time series, as necessary, to accurately analyse the datasets in our study. Firstly, let us assume that there are N measurements of process T(ti), where (ti)i∈[1,N] are irregular time samplings. Then, we define the cumulative sum Si as
Si=∑j≤iTtj.
By taking an index j and an even number k, we can define [tj,tj+k] as the interval where the fluctuation ΔTj,k is estimated. To do so, we determine the sum of T(ti) over the first and second halves of the interval by
ΔS(1)=Sj+k/2-Sj,ΔS(2)=Sj+k-Sj+k/2
and compute the average of the first half minus the average of the second half of the [tj,tj+k] interval as
ΔTj,k=2ΔS(1)tj+k/2-tj-ΔS(2)tj+k-tj+k/2.
Here, the canonical factor 2 multiplied by the difference is a calibration parameter necessary to expand the working range from the anomaly (H<0) and differences (H>0) to the Haar range -1≤H≤1; see .
This definition works for any sampling rate, but its accuracy deteriorates when the first and second halves of the interval cover disproportionate time durations. Hence, we included quality control in the uniformity of the intervals. For this purpose, we define the non-uniformity of the interval Δt as the ratio between the time samples:
ϵ=tj+k/2-tjtj+k-tj.
Then, we define a threshold ϵmin and include only Haar fluctuations on intervals that satisfy ϵmin≤ϵ≤(1-ϵmin). Decreasing ϵmin implies a precision loss and hence smoothens the resulting curve. Differently, taking ϵmin too close to 1/2 (uniform sampling case) can result in the exclusion of too many fluctuations, with the consequence that the statistics will be poor. For this study, we fixed ϵmin at 1/4.
The first-order Haar structural function in Eq. () can be readily applied to the definition of non-uniform intervals. Furthermore, the second-order structure function can be used to estimate the Haar fluctuation root-mean-square (RMS) value as
|ΔF(Δt)|21/2∝Δtξ(2)/2=ΔtH-K(2)/2,
as per . This statistic fully characterises the scaling of Gaussian processes and provides important insights for non-Gaussian datasets. Thus, we used the RMS fluctuation as an approximation to the absolute mean (that is, H≈ξ(2)/2) in this study.
Code and data availability
The data are available as a supplement to the references listed in Table . A repository containing the Python function to estimate Haar fluctuations in non-equidistant datasets is available at https://github.com/NicoAcunaR/Non-equidistant_Haar_fluctuations.
Author contributions
All authors participated in the conceptualisation of the research and the methodology. NA developed the software and visualisation and conducted the formal analysis and investigation. EV, FL, and SL provided supervision. NA prepared the original draft, and all authors contributed to reviewing and editing the final paper.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
Acknowledgements
This study was funded by ANID through FONDECYT grant nos. 1231682 and 1221220.
Financial support
This research has been supported by grants ANID/FONDECYT 1231682 and 1221220.
Review statement
This paper was edited by Laurie Menviel and reviewed by Anders Svensson and two anonymous referees.
ReferencesAbe-Ouchi, A., Saito, F., Kawamura, K., Raymo, M. E., Okuno, J., Takahashi, K., and Blatter, H.: Insolation-driven 100,000-year glacial cycles and hysteresis of ice-sheet volume, Nature, 500, 190–193, 10.1038/nature12374, 2013.Acuña Reyes, N.: Non-equidistant_Haar_fluctuations, GitHub [data set], https://github.com/NicoAcunaR/Non-equidistant_Haar_fluctuations (last access: 22 July 2024), 2024.Buizert, C., Adrian, B., Ahn, J., Albert, M., Alley, R. B., Baggenstos, Danieland Bauska, T. K., Bay, R. C., Bencivengo, B. B., Bentley, C. R., Brook, E. J., Chellman, N. J., Clow, G. D., Cole-Dai, J., Conway, H., Cravens, E., Cuffey, K. M., Dunbar, N. W., Edwards, J. S., Fegyveresi, J. M., Ferris, D. G., Fitzpatrick, J. J., Fudge, T. J., Gibson, C. J., Gkinis, V., Goetz, J. J., Gregory, S., Hargreaves, G. M., Iverson, N., Johnson, J. A., Jones, T. R., Kalk, M. L., Kippenhan, M. J., Koffman, B. G., Kreutz, K., Kuhl, T. W., Lebar, D. A., Lee, J. E., Marcott, S. A., Markle, B. R., Maselli, O. J., McConnell, J. R., McGwire, K. C., Mitchell, L. E., Mortensen, N. B., Neff, P. D., Nishiizumi, K., Nunn, R. M., Orsi, A. J., Pasteris, D. R., Pedro, J. B., Pettit, E. C., Buford Price, P., Priscu, J. C., Rhodes, R. H., Rosen, J. L., Schauer, A. J., Schoenemann, S. W., Sendelbach, P. J., Severinghaus, J. P., Shturmakov, A. J., Sigl, M., Slawny, K. R., Souney, J. M., Sowers, T. A., Spencer, M. K., Steig, E. J., Taylor, K. C., Twickler, M. S., Vaughn, B. H., Voigt, D. E., Waddington, E. D., Welten, K. C., Wendricks, A. W., White, J. W. C., Winstrup, M., Wong, G. J., and Woodruff, T. E.: Precise interpolar phasing of abrupt climate change during the last ice age, Nature, 520, 661–665, 10.1038/nature14401, 2015.Cheng, H., Edwards, R. L., Sinha, A., Spötl, C., Yi, L., Chen, S., Kelly, M., Kathayat, G., Wang, X., Li, X., Kong, X., Wang, Y., Ning, Y., and Zhang, H.: The Asian monsoon over the past 640,000 years and ice age terminations, Nature, 534, 640–646, 10.1038/nature18591, 2016.Davtian, N. and Bard, E.: Biomarker indices and concentrations and biomarker-based temperature estimates from the Iberian Margin core MD95-2042, composite Greenland atmospheric temperature record, and Antarctic stacks, PANGAEA [data set], 10.1594/PANGAEA.957377, 2023.Delmonte, B., Petit, J. R., and Maggi, V.: Glacial to Holocene implications of the new 27 000-year dust record from the EPICA Dome C (East Antarctica) ice core, Clim. Dynam., 18, 647–660, 10.1007/s00382-001-0193-9, 2002.EPICA Community Members, Barbante, C., Barnola, J.-M., Becagli, S., Beer, J., Bigler, M., Boutron, C., Blunier, T., Castellano, E., Cattani, O., Chappellaz, J., Dahl-Jensen, D., Debret, M., Delmonte, B., Dick, D., Falourd, S., Faria, S., Federer, U., Fischer, H., Freitag, J., Frenzel, A., Fritzsche, D., Fundel, F., Gabrielli, P., Gaspari, V., Gersonde, R., Graf, W., Grigoriev, D., Hamann, I., Hansson, M., Hoffmann, G., Hutterli, M. A., Huybrechts, P., Isaksson, E., Johnsen, S., Jouzel, J., Kaczmarska, M., Karlin, T., Kaufmann, P., Kipfstuhl, S., Kohno, M., Lambert, F., Lambrecht, A., Lambrecht, A., Landais, A., Lawer, G., Leuenberger, M., Littot, G., Loulergue, L., Lüthi, D., Maggi, V., Marino, F., Masson-Delmotte, V., Meyer, H., Miller, H., Mulvaney, R., Narcisi, B., Oerlemans, J., Oerter, H., Parrenin, F., Petit, J.-R., Raisbeck, G., Raynaud, D., Röthlisberger, R., Ruth, U., Rybak, O., Severi, M., Schmitt, J., Schwander, J., Siegenthaler, U., Siggaard-Andersen, M.-L., Spahni, R., Steffensen, J. P., Stenni, B., Stocker, T. F., Tison, J.-L., Traversi, R., Udisti, R., Valero-Delgado, F., van den Broeke, M. R., van de Wal, R. S. W., Wagenbach, D., Wegner, A., Weiler, K., Wilhelms, F., Winther, J.-G., and Wolff, E.: One-to-one coupling of glacial climate variability in Greenland and Antarctica, Nature, 444, 195–198, 10.1038/nature05301, 2006.Erhardt, T., Bigler, M., Federer, U., Gfeller, G., Leuenberger, D., Stowasser, O., Röthlisberger, R., Schüpbach, S., Ruth, U., Twarloh, B., Wegner, A., Goto-Azuma, K., Kuramoto, T., Kjær, H. A., Vallelonga, P. T., Siggaard-Andersen, M.-L., Hansson, M. E., Benton, A. K., Fleet, L. G., Mulvaney, R., Thomas, E. R., Abram, N., Stocker, T. F., and Fischer, H.: High-resolution aerosol concentration data from the Greenland NorthGRIP and NEEM deep ice cores, Earth Syst. Sci. Data, 14, 1215–1231, 10.5194/essd-14-1215-2022, 2022.Franzke, C. L. E., Barbosa, S., Blender, R., Fredriksen, H.-B., Laepple, T., Lambert, F., Nilsen, T., Rypdal, K., Rypdal, M., Scotto, M. G., Vannitsem, S., Watkins, N. W., Yang, L., and Yuan, N.: The Structure of Climate Variability Across Scales, Rev. Geophys., 58, e2019RG000657, 10.1029/2019RG000657, 2020.Fredriksen, H.-B. and Rypdal, K.: Spectral Characteristics of Instrumental and Climate Model Surface Temperatures, J. Climate, 29, 1253–1268, 10.1175/JCLI-D-15-0457.1, 2016.Fuhrer, K., Neftel, A., Anklin, M., and Maggi, V.: Continuous measurements of hydrogen peroxide, formaldehyde, calcium and ammonium concentrations along the new grip ice core from summit, Central Greenland, Atmos. Environ. A, 27, 1873–1880, 10.1016/0960-1686(93)90292-7, 1993.Gkinis, V., Vinther, B. M., Popp, T. J., Faber, A.-K., Holme, C. T., Jensen, C. M., Lanzky, M., Lütt, A. M., Mandrakis, V., Ørum, N. O., Pedersen, A.-S., Vaxevani, N., Weng, Y., Capron, E., Dahl-Jensen, D., Hörhold, M., Jones, T. R., Jouzel, J., Landais, A., Masson-Delmotte, V., Oerter, H., Rasmussen, S. O., Steen-Larsen, H. C., Steffensen, J. P., Sveinbjörnsdottir, E. A., Vaughn, B. H., and White, J.: A 120,000-year long climate record from a NW-Greenland deep ice core at ultra-high resolution, Sci. Data, 8, 141, 10.1038/s41597-021-00916-9, 2021.Guo, Z. T., Berger, A., Yin, Q. Z., and Qin, L.: Strong asymmetry of hemispheric climates during MIS-13 inferred from correlating China loess and Antarctica ice records, Clim. Past, 5, 21–31, 10.5194/cp-5-21-2009, 2009.Hébert, R., Rehfeld, K., and Laepple, T.: Comparing estimation techniques for temporal scaling in palaeoclimate time series, Nonlin. Processes Geophys., 28, 311–328, 10.5194/npg-28-311-2021, 2021.Hébert, R., Herzschuh, U., and Laepple, T.: Millennial-scale climate variability over land overprinted by ocean temperature fluctuations, Nat. Geosci., 15, 899–905, 10.1038/s41561-022-01056-4, 2022.Huybers, P. and Curry, W.: Links between annual, Milankovitch and continuum temperature variability, Nature, 441, 329–332, 10.1038/nature04745, 2006.Jacobel, A., McManus, J., Anderson, R., and Winckler, G.: Large deglacial shifts of the Pacific Intertropical Convergence Zone, Nat. Commun., 7, 10449, 10.1038/ncomms10449, 2016.Jouzel, J., Vimeux, F., Caillon, N., Delaygue, G., Hoffmann, G., Masson-Delmotte, V., and Parrenin, F.: Magnitude of isotope/temperature scaling for interpretation of central Antarctic ice cores, J. Geophys. Res.-Atmos., 108, 4361, 10.1029/2002JD002677, 2003.Jouzel, J., Masson-Delmotte, V., Cattani, O., Dreyfus, G., Falourd, S., Hoffmann, G., Minster, B., Nouet, J., Barnola, J. M., Chappellaz, J., Fischer, H., Gallet, J. C., Johnsen, S., Leuenberger, M., Loulergue, L., Luethi, D., Oerter, H., Parrenin, F., Raisbeck, G., Raynaud, D., Schilt, A., Schwander, J., Selmo, E., Souchez, R., Spahni, R., Stauffer, B., Steffensen, J. P., Stenni, B., Stocker, T. F., Tison, J. L., Werner, M., and Wolff, E. W.: Orbital and Millennial Antarctic Climate Variability over the Past 800,000 Years, Science, 317, 793–796, 10.1126/science.1141038, 2007.Kindler, P., Guillevic, M., Baumgartner, M., Schwander, J., Landais, A., and Leuenberger, M.: Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core, Clim. Past, 10, 887–902, 10.5194/cp-10-887-2014, 2014.Laepple, T. and Huybers, P.: Reconciling discrepancies between Uk37 and Mg/Ca reconstructions of Holocene marine temperature variability, Earth Planet. Sc. Lett., 375, 418–429, 10.1016/j.epsl.2013.06.006, 2013.Laepple, T. and Huybers, P.: Ocean surface temperature variability: Large model–data differences at decadal and longer periods, P. Natl. Acad. Sci. USA, 111, 16682–16687, 10.1073/pnas.1412077111, 2014.Lambert, F., Bigler, M., Steffensen, J. P., Hutterli, M., and Fischer, H.: Centennial mineral dust variability in high-resolution ice core data from Dome C, Antarctica, Clim. Past, 8, 609–623, 10.5194/cp-8-609-2012, 2012.Lambert, F., Kug, J.-S., Park, R. J., Mahowald, N., Winckler, G., Abe-Ouchi, A., O'ishi, R., Takemura, T., and Lee, J.-H.: The role of mineral-dust aerosols in polar temperature amplification, Nat. Clim. Change, 3, 487–491, 10.1038/nclimate1785, 2013.Lamy, F., Gersonde, R., Winckler, G., Esper, O., Jaeschke, A., Kuhn, G., Ullermann, J., Martinez-Garcia, A., Lambert, F., and Kilian, R.: Increased Dust Deposition in the Pacific Southern Ocean During Glacial Periods, Science, 343, 403–407, 10.1126/science.1245424, 2014.Lisiecki, L. E. and Raymo, M. E.: A Pliocene-Pleistocene stack of 57 globally distributed benthic δ18O records, Paleoceanography, 20, PA1003, 10.1029/2004PA001071, 2005.Lovejoy, S.: A voyage through scales, a missing quadrillion and why the climate is not what you expect, Clim. Dynam., 44, 3187–3210, 10.1007/s00382-014-2324-0, 2015.Lovejoy, S.: Review article: Scaling, dynamical regimes, and stratification. How long does weather last? How big is a cloud?, Nonlin. Processes Geophys., 30, 311–374, 10.5194/npg-30-311-2023, 2023.Lovejoy, S. and Lambert, F.: Spiky fluctuations and scaling in high-resolution EPICA ice core dust fluxes, Clim. Past, 15, 1999–2017, 10.5194/cp-15-1999-2019, 2019.Lovejoy, S. and Schertzer, D.: Haar wavelets, fluctuations and structure functions: convenient choices for geophysics, Nonlin. Processes Geophys., 19, 513–527, 10.5194/npg-19-513-2012, 2012.Lovejoy, S. and Schertzer, D. The Weather and Climate: Emergent Laws and Multifractal Cascades, Cambridge University Press, ISBN 9781139093811, 10.1017/CBO9781139093811, 2013.Lovejoy, S. and Varotsos, C.: Scaling regimes and linear/nonlinear responses of last millennium climate to volcanic and solar forcings, Earth Syst. Dynam., 7, 133–150, 10.5194/esd-7-133-2016, 2016.Lyle, M., Heusser, L., Ravelo, C., Andreasen, D., Olivarez Lyle, A., and Diffenbaugh, N.: Pleistocene water cycle and eastern boundary current processes along the California continental margin, Paleoceanography, 25, PA4211, 10.1029/2009PA001836, 2010.Maher, B., Prospero, J., Mackie, D., Gaiero, D., Hesse, P., and Balkanski, Y.: Global connections between aeolian dust, climate and ocean biogeochemistry at the present day and at the last glacial maximum, Earth-Sci. Rev., 99, 61–97, 10.1016/j.earscirev.2009.12.001, 2010.Markle, B., Steig, E., Roe, G., Winckler, G., and McConnell, J.: Concomitant variability in high-latitude aerosols, water isotopes and the hydrologic cycle, Na. Geosci., 11, 853–859, 10.1038/s41561-018-0210-9, 2018.Masson-Delmotte, V., Kageyama, M., Braconnot, P., Charbit, S., Krinner, G., Ritz, C., Guilyardi, E., Jouzel, J., Abe-Ouchi, A., Crucifix, M., Gladstone, R. M., Hewitt, C. D., Kitoh, A., LeGrande, A. N., Marti, O., Merkel, U., Motoi, T., Ohgaito, R., Otto-Bliesner, B., Peltier, W. R., Ross, I., Valdes, P. J., Vettoretti, G., Weber, S. L., Wolk, F., and Yu, Y.: Past and future polar amplification of climate change: climate model intercomparisons and ice-core constraints, Clim. Dynam., 26, 513–529, 10.1007/s00382-005-0081-9, 2006.Mitchell, J.: An overview of climatic variability and its causal mechanisms, Quatern. Res., 6, 481–493, 10.1016/0033-5894(76)90021-1, 1976.Mohtadi, M., Lückge, A., Steinke, S., Groeneveld, J., Hebbeln, D., and Westphal, N.: Late Pleistocene surface and thermocline conditions of the eastern tropical Indian Ocean, Quaternary Sci. Rev., 29, 887–896, 10.1016/j.quascirev.2009.12.006, 2010.Naafs, B. D. A., Hefter, J., Grützner, J., and Stein, R.: Warming of surface waters in the mid-latitude North Atlantic during Heinrich events, Paleoceanography, 28, 153–163, 10.1029/2012PA002354, 2013.Nilsen, T., Rypdal, K., and Fredriksen, H.-B.: Are there multiple scaling regimes in Holocene temperature records?, Earth Syst. Dynam., 7, 419–439, 10.5194/esd-7-419-2016, 2016.North Greenland Ice Core Project members, Andersen, K. K., Azuma, N., Barnola, J.-M., Bigler, M., Biscaye, P., Caillon, N., Chappellaz, J., Clausen, H. B., Dahl-Jensen, D., Fischer, H., Flückiger, J., Fritzsche, D., Fujii, Y., Goto-Azuma, K., Grønvold, K., Gundestrup, N. S., Hansson, M., Huber, C., Hvidberg, C. S., Johnsen, S. J., Jonsell, U., Jouzel, J., Kipfstuhl, S., Landais, A., Leuenberger, M., Lorrain, R., Masson-Delmotte, V., Miller, H., Motoyama, H., Narita, H., Popp, T., Rasmussen, S. O., Raynaud, D., Rothlisberger, R., Ruth, U., Samyn, D., Schwander, J., Shoji, H., Siggard-Andersen, M.-L., Steffensen, J. P., Stocker, T., Sveinbjörnsdóttir, A. E., Svensson, A., Takata, M., Tison, J.-L., Thorsteinsson, T., Watanabe, O., Wilhelms, F., and White, J. W. C.: High-resolution record of Northern Hemisphere climate extending into the last interglacial period, Nature, 431, 147–151, 10.1038/nature02805, 2004.Petit, J. R., Jouzel, J., Raynaud, D., Barkov, N. I., Barnola, J.-M., Basile, I., Bender, M., Chappellaz, J., Davis, M., Delaygue, G., Delmotte, M., Kotlyakov, V. M., Legrand, M., Lipenkov, V. Y., Lorius, C., Pépin, L., Ritz, C., Saltzman, E., and Stievenard, M.: Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica, Nature, 399, 429–436, 10.1038/20859, 1999.Rasmussen, S. O., Bigler, M., Blockley, S. P., Blunier, T., Buchardt, S. L., Clausen, H. B., Cvijanovic, I., Dahl-Jensen, D., Johnsen, S. J., Fischer, H., Gkinis, V., Guillevic, M., Hoek, W. Z., Lowe, J. J., Pedro, J. B., Popp, T., Seierstad, I. K., Steffensen, J. P., Svensson, A. M., Vallelonga, P., Vinther, B. M., Walker, M. J., Wheatley, J. J., and Winstrup, M.: A stratigraphic framework for abrupt climatic changes during the Last Glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy, Quaternary Sci. Rev., 106, 14–28, 10.1016/j.quascirev.2014.09.007, 2014.Rehfeld, K., Münch, T., Ho, S. L., and Laepple, T.: Global patterns of declining temperature variability from the Last Glacial Maximum to the Holocene, Nature, 554, 356–359, 10.1038/nature25454, 2018.Reschke, M., Kröner, I., and Laepple, T.: Testing the consistency of Holocene and Last Glacial Maximum spatial correlations in temperature proxy records, J. Quaternary Sci., 36, 20–28, 10.1002/jqs.3245, 2021. Ruddiman, W.: Earth's Climate: Past and Future, W. H. Freeman, ISBN 13:978-0-7167-8490-6, ISBN 10:0-7167-8490-4, 2008.Ruth, U., Barnola, J.-M., Beer, J., Bigler, M., Blunier, T., Castellano, E., Fischer, H., Fundel, F., Huybrechts, P., Kaufmann, P., Kipfstuhl, S., Lambrecht, A., Morganti, A., Oerter, H., Parrenin, F., Rybak, O., Severi, M., Udisti, R., Wilhelms, F., and Wolff, E.: “EDML1”: a chronology for the EPICA deep ice core from Dronning Maud Land, Antarctica, over the last 150 000 years, Clim. Past, 3, 475–484, 10.5194/cp-3-475-2007, 2007.Sato, H., Suzuki, T., Hirabayashi, M., Iizuka, Y., Motoyama, H., and Fujii, Y.: Mineral and Sea-Salt Aerosol Fluxes over the Last 340 kyr Reconstructed from the Total Concentration of Al and Na in the Dome Fuji Ice Core, Atmos. Clim. Sci., 3, 186–192, 10.4236/acs.2013.32020, 2013.Schüpbach, S., Federer, U., Kaufmann, P. R., Albani, S., Barbante, C., Stocker, T. F., and Fischer, H.: High-resolution mineral dust and sea ice proxy records from the Talos Dome ice core, Clim. Past, 9, 2789–2807, 10.5194/cp-9-2789-2013, 2013. Simonsen, M., Baccolo, G., Blunier, T., Borunda, A., Delmonte, B., Frei, R., Goldstein, S., Grinsted, A., Kjær, H., Sowers, T., Svensson, A., Vinther, B., Vladimirova, D., Winckler, G., Winstrup, M., and Vallelonga, P.: East Greenland ice core dust record reveals timing of Greenland ice sheet advance and retreat, Nat. Commun., 10, 4494, 10.1038/s41467-019-12546-2, 2019.Spiridonov, A. and Lovejoy, S.: Life rather than climate influences diversity at scales greater than 40 million years, Nature, 607, 307–312, 10.1038/s41586-022-04867-y, 2022.Stocker, T. F. and Johnsen, S. J.: A minimum thermodynamic model for the bipolar seesaw, Paleoceanography, 18, 1087, 10.1029/2003PA000920, 2003.Svensson, A., Andersen, K. K., Bigler, M., Clausen, H. B., Dahl-Jensen, D., Davies, S. M., Johnsen, S. J., Muscheler, R., Parrenin, F., Rasmussen, S. O., Röthlisberger, R., Seierstad, I., Steffensen, J. P., and Vinther, B. M.: A 60 000 year Greenland stratigraphic ice core chronology, Clim. Past, 4, 47–57, 10.5194/cp-4-47-2008, 2008.Tierney, J. E., Russell, J. M., Huang, Y., Damsté, J. S. S., Hopmans, E. C., and Cohen, A. S.: Northern Hemisphere Controls on Tropical Southeast African Climate During the Past 60,000 Years, Science, 322, 252–255, 10.1126/science.1160485, 2008.Tzedakis, P. C., Margari, V., Skinner, L. C., Menviel, L., Capron, E., Rhodes, R. H., Mleneck-Vautravers, M. J., Ezat, M. M., Martrat, B., Grimalt, J. O., and Hodell, D. A.: Palaeoceanographic, sediment and pollen data from sediment core MD01-2444 on the Portuguese Margin, Marine Isotope Stage 3, PANGAEA [data set], 10.1594/PANGAEA.919846, 2020.Veres, D., Bazin, L., Landais, A., Kele, H. T. M., Lemieux-Dudon, B., Parrenin, F., Martinerie, P., Blayo, E., Blunier, T., Capron, E., Chappellaz, J., Rasmussen, S. O., Severi, M., Svensson, A., Vinther, B., and Wolff, E. W.: The Antarctic ice core chronology (AICC2012): An optimized multi-parameter and multi-site dating approach for the last 120 thousand years, Clim. Past, 9, 1733–1748, 10.5194/cp-9-1733-2013, 2013.von der Heydt, A. S., Ashwin, P., Camp, C. D., Crucifix, M., Dijkstra, H. A., Ditlevsen, P., and Lenton, T. M.: Quantification and interpretation of the climate variability record, Global Planet. Change, 197, 103399, 10.1016/j.gloplacha.2020.103399, 2021.Wunsch, C.: The spectral description of climate change including the 100 ky energy, Clim. Dynam., 20, 353–363, 10.1007/s00382-002-0279-z, 2003.Zhu, F., Emile-Geay, J., McKay, N. P., Hakim, G. J., Khider, D., Ault, T. R., Steig, E. J., Dee, S., and Kirchner, J. W.: Climate models can correctly simulate the continuum of global-average temperature variability, P. Natl. Acad. Sci. USA, 116, 8728–8733, 10.1073/pnas.1809959116, 2019.