Abstract. The planetary boundary layer (BL) height and stratification are key parameters in determining the exchange of heat, momentum, moisture and trace gases between the surface and the free troposphere. Numerous different methods have been used to quantify the BL height and these methods have been applied to a wide variety of observational data sets obtained from different instruments and to numerical model output. We investigate the BL height at the Hyytiälä SMEAR II station in southern Finland diagnosed from radiosonde observations, a microwave radiometer (MWR) and ERA5 reanalysis. Four different algorithms are used to diagnose the BL height from the radiosondes. The diagnosed BL height is sensitive to the method used and the level of agreement, and the sign of systematic bias, between the 4 different methods depends on the surface-layer stability. For example, for very unstable situations, the median BL height diagnosed from the radiosondes varies from 600 m to 1500 m depending on which method is applied. Good agreement between the BL height in ERA5 and diagnosed from the radiosondes using Richardson number-based methods is found for almost all stability classes, suggesting that ERA5 has adequate vertical resolution near the surface to resolve the BL structure. However, ERA5 overestimates the BL height in very stable conditions highlighting the on-going challenge for numerical models to correctly resolve the stable BL. Furthermore, ERA5 BL height differs most from the radiosondes at 18 UTC suggesting ERA5 does not resolve the evening transition correctly. This study has also shown that BL height estimates from the MWR are reliable in unstable situations but often are inaccurate under stable conditions when, in comparison to ERA5 BL heights, they are much deeper. The errors in the MWR BL height estimates originate from the limitations of the manufacturers algorithm for stable conditions and also the mis-identification of the type of BL. A climatology of the annual and diurnal cycle of BL height and observed surface layer stability was created. The shallowest (353 m) monthly median BL height occurs in February and the deepest (576 m) in June. In winter there is no diurnal cycle in BL height, unstable BLs are rare yet so are very stable BLs. The shallowest BLs occur at night in spring and summer and very stable conditions are most common at night in the warm season. Finally, using ERA5 gridded data we determined that the BL height observed at Hyytiälä is representative of most land areas in southern and central Finland. However, the spatial variability of the BL height is largest during daytime in summer reducing the area over which BL height observations from Hyytiälä would be representative of.