S U M M A R YDeviation of seismic surface waves from the great-circle between source and receiver is illustrated by the anomalies in the arrival angle, that is the difference between the observed backazimuth of the incident waves and the great-circle. Such arrival angle anomalies have been known for decades, but observations remain scattered. We present a systematic study of arrival angle anomalies of fundamental mode Rayleigh waves (20-100 s period interval) from 289 earthquakes and recorded by a broadband network LAPNET, located in northern Finland. These observations are compared with those of full waveform synthetic seismograms for the same events, calculated in a 3-D Earth and also compared with those of seismograms obtained by ambient noise correlation. The arrival angle anomalies for individual events are complex, and have significant variations with period. On average, the mean absolute deviation decreases from ∼9 • at 20 s period to ∼3 • at 100 s period. The synthetic seismograms show the same evolution, albeit with somewhat smaller deviations. While the arrival angle anomalies are fairly well simulated at long periods, the deviations at short periods are very poorly modelled, demonstrating the importance of the continuous improvement of global crustal models. At 20-30 s period, both event data and numerical simulations have strong multipathing, and relative amplitude changes between different waves will induced differences in deviations between very closely located events. The source mechanism has only limited influence on the deviations, demonstrating that they are directly linked to propagation effects, including near-field effects in the source area. This observation is confirmed by the comparison with seismic noise correlation records, that is where the surface waves correspond to those emitted by a point source at the surface, as the two types of observations are remarkably similar in the cases where earthquakes are located close to seismic stations. This agreement additionally confirms that the noise correlations capture the complex surface wave propagation.