Soil temperature changes in cold regions can have great impacts on the land surface energy and water balance, and hence changes in weather and climate, surface and subsurface hydrology and ecosystem. We investigate the spatiotemporal variations of annual soil temperature at depths of 0, 5, 10, 15, 20, and 40 cm during 1983-2013 using observations at 85 stations on the Tibetan Plateau (TP). Our results show that the climatological soil temperatures exhibit a similar spatial pattern among different depths and they are generally higher than surface air temperature at the individual stations. Spatially averaged soil temperature show that the TP has experienced significant warming trends at all six depths during 1983-2013, and the soil at 0-cm depth has the fastest warming rate among all the six layers and the surface air temperature. The first leading mode of joint empirical orthogonal function (EOF) analysis exhibits a spatially prevailing warming pattern across the six depths. This plateau-wide soil warming correlates very well with surface air temperature and sea surface temperature in response to increasing radiative forcing caused by greenhouse gases. The joint EOF2 displays a southeastern-northwestern dipole pattern on the TP in the interannual-decadal variability of soil temperature at all layers, which appears to be related to the warm season precipitation and anomalous atmospheric circulations. The spatial difference of soil warming rates across stations on the TP is associated primarily with the spatial distribution of precipitation (mainly rainfall), with vegetation, snowfall and elevation playing a rather limited role.