Global tree-ring response and inferred climate variation following the mid-thirteenth century Samalas eruption

The largest explosive volcanic eruption of the Common Era in terms of estimated sulphur yield to the stratosphere was identified in glaciochemical records 40 years ago, and dates to the mid-thirteenth century. Despite eventual attribution to the Samalas (Rinjani) volcano in Indonesia, the eruption date remains uncertain, and the climate response only partially understood. Seeking a more global perspective on summer surface temperature and hydroclimate change following the eruption, we present an analysis of 249 tree-ring chronologies spanning the thirteenth century and representing all continents except Antarctica. Of the 170 predominantly temperature sensitive high-frequency chronologies, the earliest hints of boreal summer cooling are the growth depressions found at sites in the western US and Canada in 1257 CE. If this response is a result of Samalas, it would be consistent with an eruption window of circa May–July 1257 CE. More widespread summer cooling across the mid-latitudes of North America and Eurasia is pronounced in 1258, while records from Scandinavia and Siberia reveal peak cooling in 1259. In contrast to the marked post-Samalas temperature response at high-elevation sites in the Northern Hemisphere, no strong hydroclimatic anomalies emerge from the 79 precipitation-sensitive chronologies. Although our findings remain spatially biased towards the western US and central Europe, and growth-climate response patterns are not always dominated by a single meteorological factor, this study offers a global proxy framework for the evaluation of paleoclimate model simulations.

Seasonality in the relationship between equatorial-mean heat content and interannual eastern equatorial Atlantic sea surface temperature variability

Interannual sea surface temperature (SST) variations in the tropical Atlantic Ocean lead to anomalous atmospheric circulation and precipitation patterns with important ecological and socioeconomic consequences for the semiarid regions of sub-Saharan Africa and northeast Brazil. This interannual SST variability is characterized by three modes: an Atlantic meridional mode featuring an anomalous cross-equatorial SST gradient that peaks in boreal spring; an Atlantic zonal mode (Atlantic Niño mode) with SST anomalies in the eastern equatorial Atlantic cold tongue region that peaks in boreal summer; and a second zonal mode of variability with eastern equatorial SST anomalies peaking in boreal winter. Here we investigate the extent to which there is any seasonality in the relationship between equatorial warm water recharge and the development of eastern equatorial Atlantic SST anomalies. Seasonally stratified cross-correlation analysis between eastern equatorial Atlantic SST anomalies and equatorial heat content anomalies (evaluated using warm water volume and sea surface height) indicate that while equatorial heat content changes do occasionally play a role in the development of boreal summer Atlantic zonal mode events, they contribute more consistently to Atlantic Niño II, boreal winter events. Event and composite analysis of ocean adjustment with a shallow water model suggest that the warm water volume anomalies originate mainly from the off-equatorial northwestern Atlantic, in agreement with previous studies linking them to anomalous wind stress curl associated with the Atlantic meridional mode.

Pacific multidecadal (50–70 year) variability instigated by volcanic forcing during the Little Ice Age (1250–1850)

The Pacific decadal oscillation (PDO) is the leading mode of decadal climate variability over the North Pacific. However, it remains unknown to what extent external forcings can influence the PDO’s periodicity and magnitude over the past 2000 years. We show that the paleo-assimilation products (LMR) and proxy data suggest a 20–40 year PDO occurred during both the Mediaeval Climate Anomaly (MCA, ~ 750–1150) and Little Ice Age (LIA, ~ 1250–1850) while a salient 50–70 year variance peak emerged during the LIA. These results are reproduced well by the CESM simulations in the all-forcing (AF) and single volcanic forcing (Vol) experiments. We show that the 20–40 year PDO is an intrinsic mode caused by internal variability but the 50–70 year PDO during the LIA is a forced mode primarily shaped by volcanic forcing. The intrinsic mode develops in tandem with tropical ENSO-like anomalies, while the forced mode develops from the western Pacific and unrelated to tropical sea surface temperature anomalies. The volcanism-induced land–sea thermal contrast may trigger anomalous northerlies over the western North Pacific (WNP), leading to reduced northward heat transport and the cooling in the Kuroshio–Oyashio Extension (KOE), generating the forced mode. A 50–70 year Atlantic multidecadal oscillation founded during the LIA under volcanic forcing may also contribute to the forced mode. These findings shed light on the interplay between the internal variability and external forcing and the present and future changes of the PDO.