The eastern Alpine crust has been shaped by the continental collision of the European and Adriatic plates beginning at 35 Ma and was affected by a major reorganization after 20 Ma. To better understand how the eastern Alpine surface structures link with deep seated processes, we analyze the depth-dependent seismic anisotropy based on Rayleigh wave propagation. Ambient noise recordings are evaluated to extract Rayleigh wave phase dispersion measurements. These are inverted in a two step approach for the azimuthally anisotropic shear velocity structure. Both steps are performed with a reversible jump Markov chain Monte Carlo (rj-McMC) approach that estimates data errors and propagates the modeled uncertainties from the phase velocity maps into the depth inversion. A two layer structure of azimuthal anisotropy is imaged in the Alpine crust, with an orogen-parallel upper crust and approximately orogen-perpendicular layer in the lower crust and the uppermost mantle. In the upper layer, the anisotropy tends to follow major fault lines and may thus be an apparent, structurally driven anisotropy. The main foliation and fold axis orientations might contribute to the anisotropy. In the lower crust, the N-S orientation of the fast axis is mostly confined to regions north of the Periadriatic Fault and may be related to European subduction. Outside the orogen, no clearly layered structure is identified. The anisotropy pattern in the northern Alpine foreland is found to be similar compared to SKS studies which is an indication of very homogeneous fast axis directions throughout the crust and the upper mantle. The formation of the European Alps is due to the continental collision of the Adriatic and the European plates which started approximately 35 million years ago and is still ongoing. To better understand how the plates behaved during this collision process, how the crust was deformed and how distinct features like the Tauern Window were formed, we image the upper 70 km under the Eastern Alps using seismic background noise. We focus on the anisotropic shear-velocity structure of the underground, because it provides valuable additional information that can be related to the fault geometry or the main movement direction of the material. We present a probabilistic approach that attributes errors to the measured data and propagates this error through the two steps of the 3D imaging procedure, to be able to assess the uncertainty of the final model. The results show that there are two anisotropic layers in the crust, an upper crustal layer that is dominated by faster wave propagation parallel to the orientation of the mountain chain, and a lower crustal/uppermost mantle layer that has a dominantly perpendicular fast propagation direction. A new, probabilistic method is presented to obtain an azimuthally anisotropic shear-velocity model from Rayleigh waves The model shows a two-layer anisotropic structure in the eastern Alpine crust