Speleothem δ18O records provide valuable information about past continental environmental and climatic conditions, although their interpretation is often not straightforward. Here we evaluate a compilation of late Holocene speleothem δ18O records using a Monte Carlo based Principal Component Analysis (MC-PCA) method that accounts for uncertainties in individual speleothem age models and for the variable temporal resolution of each δ18O record. The MC-PCA approach permits not only the identification of temporally coherent changes in speleothem δ18O; it also facilitates their graphical depiction and evaluation of their spatial coherency. The MC-PCA method was applied to 11 Holocene speleothem δ18O records that span most of the European continent (apart from the circum-Mediterranean region). We observe a common (shared) mode of speleothem δ18O variability that suggests millennial-scale coherency and cyclicity during the last 4.5 ka. These changes are likely caused by variability in atmospheric circulation akin to that associated with the North Atlantic Oscillation, reflecting meridionally shifted westerlies. We argue that these common large-scale variations in European speleothem δ18O records are in phase with changes in the North Atlantic Ocean circulation indicated by the vigour of the Iceland Scotland Overflow Water (ISOW), the strength of the subpolar gyre (SPG) and an ocean stacked North Atlantic ice rafted debris (IRD) index. Based on a recent modelling study, we conclude that these changes in the North Atlantic circulation history may be caused by wind stress on the ocean surface driven by shifted westerlies. However, the mechanisms that ultimately force the westerlies remain unclear.