Description: Proxy reconstructions of tropical Atlantic sea surface temperature (SST) that extend beyond the period of instrumental observations have primarily focused on centennial to millennial variability rather than on seasonal to multidecadal variability. Here we present monthly-resolved records of Sr/Ca (a proxy of SST) from fossil annually-banded Diploria strigosa corals from Bonaire (southern Caribbean Sea). The individual corals provide time-windows of up to 68 years length, and the total number of 295 years of record allows for assessing the natural range of seasonal to multidecadal SST variability in the western tropical Atlantic during snapshots of the mid- to late Holocene. Comparable to modern climate, the coral Sr/Ca records reveal that mid- to late Holocene SST was characterised by clear seasonal cycles, persistent quasi-biennial and prominent interannual as well as inter- to multidecadal-scale variability. However, the magnitude of SST variations on these timescales has varied over the last 6.2 ka. The coral records show increased seasonality during the mid-Holocene consistent with climate model simulations indicating that southern Caribbean SST seasonality is induced by insolation changes on orbital timescales, whereas internal dynamics of the climate system play an important role on shorter timescales. Interannual SST variability is linked to ocean–atmosphere interactions of Atlantic and Pacific origin. Pronounced interannual variability in the western tropical Atlantic is indicated by a 2.35 ka coral, possibly related to a strengthening of the variability of the El Niño/Southern Oscillation throughout the Holocene. Prominent inter- to multidecadal SST variability is evident in the coral records and slightly more pronounced in the mid-Holocene. We finally argue that our coral data provide a target for studying Holocene climate variability on seasonal and interannual to multidecadal timescales, when using further numerical models and high-resolution proxy data.

Description: Previous studies on surface temperature reconstructions for the last 2000 years (2 k) revealed a long-term cooling trend for the last millennium in comparison to the previous millennium. However, knowledge on the decadal- to centennial-scale variability in sea surface temperature and the underlying governing mechanisms throughout the period is limited. We reconstructed high-resolution continuous sea surface temperature changes over the last 2 k in the northwest Pacific margin based on the alkenone unsaturation index. Our alkenone temperature record revealed enhanced and more rapidly changing climate variability during the last millennium (approximately 1200–1850 Common Era) than during the previous millennium. Cold and hot extremes also occurred more frequently during the last millennium. The enhanced and rapidly changing climate variability appears to be associated with frequent volcanic eruptions and grand solar minima. The reconstructed surface temperature variability tends to be associated with variations in the East Asia summer monsoon and the Pacific Decadal Oscillation, implying that these variations are also enhanced in the last millennium than in the previous millennium.

Description: Highlights: • Climate model sensitivity experiments are performed using state-of-the-art ice sheet and freshwater reconstructions • Declining Northern Hemisphere ice sheets increase the sensitivity of the AMOC to North Atlantic meltwater discharge • Deglacial rise in atmospheric CO2 concentration decreases the sensitivity of the AMOC to North Atlantic meltwater discharge • Both effects provide a complementary perspective to existing explanations for abrupt AMOC transitions Abstract: The last deglaciation was characterized by a sequence of abrupt climate events thought to be linked to rapid changes in Atlantic meridional overturning circulation (AMOC). The sequence includes a weakening of the AMOC after the Last Glacial Maximum (LGM) during Heinrich Stadial 1 (HS1), which ends with an abrupt AMOC amplification at the transition to the Bølling/Allerød (B/A). This transition occurs despite persistent deglacial meltwater fluxes that counteract vigorous North Atlantic deep-water formation. Using the Earth system model COSMOS with a range of deglacial boundary conditions and reconstructed deglacial meltwater fluxes, we show that deglacial CO2 rise and ice sheet decline modulate the sensitivity of the AMOC to these fluxes. While declining ice sheets increase the sensitivity, increasing atmospheric CO2 levels tend to counteract this effect. Therefore, the occurrence of a weaker HS1 AMOC and an abrupt AMOC increase in the presence of meltwater, might be explained by these effects, as an alternative to or in combination with changes in the magnitude or routing of meltwater discharge.

Description: Previous modelling efforts have investigated climate responses to different Milankovitch forcing during Marine Isotope Stage (MIS) 13. During this time the climate has been highly variable at atmospheric CO2 concentrations of ~240 ppm. As yet, ice sheet-climate feedbacks were missing in previous studies. Therefore we use the state-of-the-art coupled climate-ice sheet model, AWI-ESM-1.2, to investigate the MIS-13 climate and corresponding Northern Hemisphere ice sheet (NHIS) evolution by performing simulations under three different astronomical configurations representing 495, 506 and 517 kyr BP. The simulated excess ice compared to present-day is mainly over the Cordillera, Arctic islands and Tibet. The global mean surface air temperature for the MIS-13 experiments have the same magnitude. At 506 kyr BP with boreal summer at perihelion, the Northern Hemisphere continents are warmer during summer than the other experiments, which could potentially inhibit the development of the ice sheets. The Cordilleran Ice Sheet is found to be especially sensitive to orbital (precession) forcing, at an intermediate CO2 level. This is probably due to its high elevation where the freezing point could be easily maintained. The other ice sheets over northeast America and Eurasia, however, are absent in our simulations. We propose that the alpine-based Cordilleran Ice Sheet is more sensitive and easier to build up than other NHISs in response to the astronomical controlled summer insolation. Dynamic surges are simulated for the Cordilleran Ice Sheet under fixed low orbital forcing. These surges due to internal ice sheet-climate feedbacks could potentially be the mechanism for the millennial scale H-like events.

Description: Highlights: • Multi-centennial oscillation with 100–200 years periods is evident in proxy data and model simulations during the Holocene. • Multi-centennial oscillation is a global signal and is more significant in the Northern Hemisphere high latitudes. • None of the external forcings is found to be the sole driver of the multi-centennial variability. • It indicates the multi-centennial oscillation may be due to potential internal drivers and essential feedbacks. Abstract: Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique opportunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120–130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120–150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemisphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combinations of forcing agents. Significant multi-centennial variability with oscillation periods between 100 and 200 years is found in all forcing scenarios, including those with only orbital forcing. The different forcings induce some variability in the system. Yet, none appear to be the predominant driver based on the spectral analysis. Solar irradiance has long been hypothesized to be a primary driver of multi-centennial variability. However, all the simulations without this forcing have shown significant multi-centennial variability. The results then indicate that internal mechanisms operate on multi-centennial timescales, and the North Atlantic-Arctic is a region of interest for this aspect.

Description: Highlights: • Temperature anomalies for the Mid-Holocene compared to preindustrial are significantly different in the low- and high-resolution versions of the atmospheric model ECHAM5 • For summer, shortwave cloud radiative forcing emerges as an important factor. • For boreal winter, differences are mainly related to circulation changes. • Anomaly differences are regionally as large as the mid-Holocene minus preindustrial temperature signals. Abstract: This study evaluates the dependence of simulated surface air temperatures on model resolution and orography for the mid-Holocene. Sensitivity experiments with the atmospheric general circulation model ECHAM5 are performed with low (∼3.75°, 19 vertical levels) and high (∼1.1°, 31 vertical levels) resolution. Results are compared to the respective preindustrial runs. It is found that the large-scale temperature anomalies for the mid-Holocene (compared to preindustrial) are significantly different in the low- and high-resolution versions. For boreal winter, differences are mainly related to circulation changes caused by the response to thermal forcing in conjunction with orographic resolution. For summer, shortwave cloud radiative forcing emerges as an important factor. The anomaly differences (low minus high resolution version) in the Northern Hemisphere are regionally as large as the anomalous mid-Holocene temperature signals. Furthermore, they depend on the applied surface boundary conditions. We conclude that the resolution matters for the Northern Hemisphere response in mid-Holocene simulations, which should be taken into account in model-model and data-model comparisons.

Description: Glacial climate is marked by abrupt, millennial-scale climate changes known as Dansgaard–Oeschger cycles. The most pronounced stadial coolings, Heinrich events, are associated with massive iceberg discharges to the North Atlantic. These events have been linked to variations in the strength of the Atlantic meridional overturning circulation. However, the factors that lead to abrupt transitions between strong and weak circulation regimes remain unclear. Here we show that, in a fully coupled atmosphere–ocean model, gradual changes in atmospheric CO2 concentrations can trigger abrupt climate changes, associated with a regime of bi-stability of the Atlantic meridional overturning circulation under intermediate glacial conditions. We find that changes in atmospheric CO2 concentrations alter the transport of atmospheric moisture across Central America, which modulates the freshwater budget of the North Atlantic and hence deep-water formation. In our simulations, a change in atmospheric CO2 levels of about 15 ppmv—comparable to variations during Dansgaard–Oeschger cycles containing Heinrich events—is sufficient to cause transitions between a weak stadial and a strong interstadial circulation mode. Because changes in the Atlantic meridional overturning circulation are thought to alter atmospheric CO2 levels, we infer that atmospheric CO2 may serve as a negative feedback to transitions between strong and weak circulation modes.

Description: Thermodynamic arguments imply that global mean rainfall increases in a warmer atmosphere; however, dynamical effects may result in more significant diversity of regional precipitation change. Here we investigate rainfall changes in the mid-Pliocene Warm Period (~ 3 Ma), a time when temperatures were 2–3ºC warmer than the pre-industrial era, using output from the Pliocene Model Intercomparison Projects phases 1 and 2 and sensitivity climate model experiments. In the Mid-Pliocene simulations, the higher rates of warming in the northern hemisphere create an interhemispheric temperature gradient that enhances the southward cross-equatorial energy flux by up to 48%. This intensified energy flux reorganizes the atmospheric circulation leading to a northward shift of the Inter-Tropical Convergence Zone and a weakened and poleward displaced Southern Hemisphere Subtropical Convergences Zones. These changes result in drier-than-normal Southern Hemisphere tropics and subtropics. The evaluation of the mid-Pliocene adds a constraint to possible future warmer scenarios associated with differing rates of warming between hemispheres.

Description: Abundant hydroclimatic evidence from western Amazonia and the adjacent Andes documents wet conditions during Heinrich Stadial 1 (HS1, 18–15 ka), a cold period in the high latitudes of the North Atlantic. This precipitation anomaly was attributed to a strengthening of the South American summer monsoon due to a change in the Atlantic interhemispheric sea surface temperature (SST) gradient. However, the physical viability of this mechanism has never been rigorously tested. We address this issue by combining a thorough compilation of tropical South American paleorecords and a set of atmosphere model sensitivity experiments. Our results show that the Atlantic SST variations alone, although leading to dry conditions in northern South America and wet conditions in northeastern Brazil, cannot produce increased precipitation over western Amazonia and the adjacent Andes during HS1. Instead, an eastern equatorial Pacific SST increase (i.e., 0.5–1.5 °C), in response to the slowdown of the Atlantic Meridional Overturning Circulation during HS1, is crucial to generate the wet conditions in these regions. The mechanism works via anomalous low sea level pressure over the eastern equatorial Pacific, which promotes a regional easterly low-level wind anomaly and moisture recycling from central Amazonia towards the Andes.

Description: The North Atlantic Oscillation is the dominant atmospheric pressure mode in the North Atlantic region and affects winter temperature and precipitation in the Mediterranean, northwest Europe, Greenland, and Asia1. The index1 that describes the sea-level pressure difference between Iceland and the Azores is correlated with a dipole precipitation pattern over northwest Europe and northwest Africa. How the North Atlantic Oscillation will develop as the Greenland ice sheet melts is unclear2. A potential past analogue is the early Holocene, during which melting ice sheets around the North Atlantic3, 4 freshened surface waters, affecting the strength of the meridional overturning circulation5. Here we present a Holocene rainfall record from northwest Africa based on speleothem δ18O and compare it against a speleothem-based rainfall record from Europe6. The two records are positively correlated during the early Holocene, followed by a shift to an anti-correlation, similar to the modern record, during the mid-Holocene. On the basis of our simulations with an Earth system model, we suggest the shift to the anti-correlation reflects a large-scale atmospheric and oceanic reorganization in response to the demise of the Laurentide ice sheet and a strong reduction of meltwater flux to the North Atlantic, pointing to a potential sensitivity of the North Atlantic Oscillation to the melting of ice sheets.