Description: Characterizing the variability across timescales is important for understanding the underlying dynamics of the Earth system. It remains challenging to do so from palaeoclimate archives since they are more often than not irregular, and traditional methods for producing timescale-dependent estimates of variability, such as the classical periodogram and the multitaper spectrum, generally require regular time sampling. We have compared those traditional methods using interpolation with interpolation-free methods, namely the Lomb–Scargle periodogram and the first-order Haar structure function. The ability of those methods to produce timescale-dependent estimates of variability when applied to irregular data was evaluated in a comparative framework, using surrogate palaeo-proxy data generated with realistic sampling. The metric we chose to compare them is the scaling exponent, i.e. the linear slope in log-transformed coordinates, since it summarizes the behaviour of the variability across timescales. We found that, for scaling estimates in irregular time series, the interpolation-free methods are to be preferred over the methods requiring interpolation as they allow for the utilization of the information from shorter timescales which are particularly affected by the irregularity. In addition, our results suggest that the Haar structure function is the safer choice of interpolation-free method since the Lomb–Scargle periodogram is unreliable when the underlying process generating the time series is not stationary. Given that we cannot know a priori what kind of scaling behaviour is contained in a palaeoclimate time series, and that it is also possible that this changes as a function of timescale, it is a desirable characteristic for the method to handle both stationary and non-stationary cases alike.

Description: Changes in heat transport associated with fluctuations in the strength of the Atlantic meridional overturning circulation (AMOC) are widely considered to affect the position of the Intertropical Convergence Zone (ITCZ), but the temporal immediacy of this teleconnection has to date not been resolved. Based on a high‐resolution marine sediment sequence over the last deglaciation, we provide evidence for a synchronous and near‐linear link between changes in the Atlantic interhemispheric sea surface temperature difference and continental precipitation over northeast Brazil. The tight coupling between AMOC strength, sea surface temperature difference, and precipitation changes over northeast Brazil unambiguously points to a rapid and proportional adjustment of the ITCZ location to past changes in the Atlantic meridional heat transport.

Description: Changes in climate variability are as important for society to address as are changes in mean climate1. Contrasting temperature variability during the Last Glacial Maximum and the Holocene can provide insights into the relationship between the mean state of the climate and its variability2,3. However, although glacial–interglacial changes in variability have been quantified for Greenland2, a global view remains elusive. Here we use a network of marine and terrestrial temperature proxies to show that temperature variability decreased globally by a factor of four as the climate warmed by 3–8 degrees Celsius from the Last Glacial Maximum (around 21,000 years ago) to the Holocene epoch (the past 11,500 years). This decrease had a clear zonal pattern, with little change in the tropics (by a factor of only 1.6–2.8) and greater change in the mid-latitudes of both hemispheres (by a factor of 3.3–14). By contrast, Greenland ice-core records show a reduction in temperature variability by a factor of 73, suggesting influences beyond local temperature or a decoupling of atmospheric and global surface temperature variability for Greenland. The overall pattern of reduced variability can be explained by changes in the meridional temperature gradient, a mechanism that points to further decreases in temperature variability in a warmer future.

Description: It is virtually certain that the mean surface temperature of the Earth will continue to increase under realistic emission scenarios, yet comparatively little is known about future changes in climate variability. This study explores changes in climate variability over the large range of climates simulated by the Coupled Model Intercomparison Project Phase 5 and 6 (CMIP5/6) and the Paleoclimate Modeling Intercomparison Project Phase 3 (PMIP3), including time slices of the Last Glacial Maximum, the mid-Holocene, and idealized experiments (1 % CO2 and abrupt4×CO2). These states encompass climates within a range of 12 ∘C in global mean temperature change. We examine climate variability from the perspectives of local interannual change, coherent climate modes, and through compositing extremes. The change in the interannual variability of precipitation is strongly dependent upon the local change in the total amount of precipitation. At the global scale, temperature variability is inversely related to mean temperature change on intra-seasonal to multidecadal timescales. This decrease is stronger over the oceans, while there is increased temperature variability over subtropical land areas (40∘ S–40∘ N) in warmer simulations. We systematically investigate changes in the standard deviation of modes of climate variability, including the North Atlantic Oscillation, the El Niño–Southern Oscillation, and the Southern Annular Mode, with global mean temperature change. While several climate modes do show consistent relationships (most notably the Atlantic Zonal Mode), no generalizable pattern emerges. By compositing extreme precipitation years across the ensemble, we demonstrate that the same large-scale modes influencing rainfall variability in Mediterranean climates persist throughout paleoclimate and future simulations. The robust nature of the response of climate variability, between cold and warm climates as well as across multiple timescales, suggests that observations and proxy reconstructions could provide a meaningful constraint on climate variability in future projections.

Description: A comprehensive database of paleoclimate records is needed to place recent warming into the longer-term context of natural climate variability. We present a global compilation of quality-controlled, published, temperature-sensitive proxy records extending back 12,000 years through the Holocene. Data were compiled from 679 sites where time series cover at least 4000 years, are resolved at sub-millennial scale (median spacing of 400 years or finer) and have at least one age control point every 3000 years, with cut-off values slackened in data-sparse regions. The data derive from lake sediment (51%), marine sediment (31%), peat (11%), glacier ice (3%), and other natural archives. The database contains 1319 records, including 157 from the Southern Hemisphere. The multi-proxy database comprises paleotemperature time series based on ecological assemblages, as well as biophysical and geochemical indicators that reflect mean annual or seasonal temperatures, as encoded in the database. This database can be used to reconstruct the spatiotemporal evolution of Holocene temperature at global to regional scales, and is publicly available in Linked Paleo Data (LiPD) format.

Description: Proxy records from climate archives provide evidence about past climate changes, but the recorded signal is affected by non-climate-related effects as well as time uncertainty. As proxy-based climate reconstructions are frequently used to test climate models and to quantitatively infer past climate, we need to improve our understanding of the proxy record signal content as well as the uncertainties involved. In this study, we empirically estimate signal-to-noise ratios (SNRs) of temperature proxy records used in global compilations of the middle to late Holocene (last 6000 years). This is achieved through a comparison of the correlation of proxy time series from nearby sites of three compilations and model time series extracted at the proxy sites from two transient climate model simulations: a Holocene simulation of the ECHAM5/MPI-OM model and the Holocene part of the TraCE-21ka simulation. In all comparisons, we found the mean correlations of the proxy time series on centennial to millennial timescales to be low (R〈0.2), even for nearby sites, which resulted in low SNR estimates. The estimated SNRs depend on the assumed time uncertainty of the proxy records, the timescale analysed, and the model simulation used. Using the spatial correlation structure of the ECHAM5/MPI-OM simulation, the estimated SNRs on centennial timescales ranged from 0.05 – assuming no time uncertainty – to 0.5 for a time uncertainty of 400 years. On millennial timescales, the estimated SNRs were generally higher. Use of the TraCE-21ka correlation structure generally resulted in lower SNR estimates than for ECHAM5/MPI-OM. As the number of available high-resolution proxy records continues to grow, a more detailed analysis of the signal content of specific proxy types should become feasible in the near future. The estimated low signal content of Holocene temperature compilations should caution against over-interpretation of these multi-proxy and multisite syntheses until further studies are able to facilitate a better characterisation of the signal content in paleoclimate records.

Description: Earth's climate can be understood as a dynamical system that changes due to external forcing and internal couplings. Essential climate variables, such as surface air temperature, describe this dynamics. Our current interglacial, the Holocene (11 700 yr ago to today), has been characterized by small variations in global mean temperature prior to anthropogenic warming. However, the mechanisms and spatiotemporal patterns of fluctuations around this mean, called temperature variability, are poorly understood despite their socioeconomic relevance for climate change mitigation and adaptation. Here we examine discrepancies between temperature variability from model simulations and paleoclimate reconstructions by categorizing the scaling behavior of local and global surface air temperature on the timescale of years to centuries. To this end, we contrast power spectral densities (PSD) and their power-law scaling using simulated and observation-based temperature series of the last 6000 yr. We further introduce the spectral gain to disentangle the externally forced and internally generated variability as a function of timescale. It is based on our estimate of the joint PSD of radiative forcing, which exhibits a scale break around the period of 7 yr. We find that local temperature series from paleoclimate reconstructions show a different scaling behavior than simulated ones, with a tendency towards stronger persistence (i.e., correlation between successive values within a time series) on periods of 10 to 200 yr. Conversely, the PSD and spectral gain of global mean temperature are consistent across data sets. Our results point to the limitation of climate models to fully represent local temperature statistics over decades to centuries. By highlighting the key characteristics of temperature variability, we pave a way to better constrain possible changes in temperature variability with global warming and assess future climate risks.

Description: Perturbations in stratospheric aerosol due to explosive volcanic eruptions are a primary contributor to natural climate variability. Observations of stratospheric aerosol are available for the past decades, and information from ice cores has been used to derive estimates of stratospheric sulfur injections and aerosol optical depth over the Holocene (approximately 10,000 BP to present) and into the last glacial period, extending back to 60,000 BP. Tephra records of past volcanism, compared to ice cores, are less complete, but extend much further into the past. To support model studies of the potential impacts of explosive volcanism on climate variability over across timescales, we present here an ensemble reconstruction of volcanic stratospheric sulfur injection (VSSI) over the last 130,000 years that is based primarily on terrestrial and marine tephra records. VSSI values are computed as a simple function of eruption magnitude, based on VSSI estimates from ice cores and satellite observations for identified eruptions. To correct for the incompleteness of the tephra record we include stochastically generated synthetic eruptions, assuming a constant background eruption frequency from the ice core Holocene record. While the reconstruction often differs from ice core estimates for specific eruptions due to uncertainties in the data used and reconstruction method, it shows good agreement with an ice core based VSSI reconstruction in terms of millennial-scale cumulative VSSI variations over the Holocene. The PalVol reconstruction provides a new basis to test the contributions of forced vs. unforced natural variability to the spectrum of climate, and the mechanisms leading to abrupt transitions in the palaeoclimate record with low-to-high complexity climate models. The PalVol volcanic forcing reconstruction is available at https://doi.org/10.26050/WDCC/PalVolv1 (Toohey, Schindlbeck-Belo, 2023).

Description: The geological record shows that abrupt changes in the Earth system can occur on timescales short enough to challenge the capacity of human societies to adapt to environmental pressures. In many cases, abrupt changes arise from slow changes in one component of the Earth system that eventually pass a critical threshold, or tipping point, after which impacts cascade through coupled climate–ecological–social systems. The chance of detecting abrupt changes and tipping points increases with the length of observations. The geological record provides the only long-term information we have on the conditions and processes that can drive physical, ecological and social systems into new states or organizational structures that may be irreversible within human time frames. Here, we use well-documented abrupt changes of the past 30 kyr to illustrate how their impacts cascade through the Earth system. We review useful indicators of upcoming abrupt changes, or early warning signals, and provide a perspective on the contributions of palaeoclimate science to the understanding of abrupt changes in the Earth system.