The Noah-MP land surface model with modified thermal roughness length parameterization scheme was used for the simulation of land surface process on the Tibetan Plateau (TP) during 1981–2018, using the high-resolution China Meteorological Forcing Data. Characteristics of distributions and changes of surface sensible and latent heat (SH and LH) fluxes on the TP were analyzed using the simulation results after verification against the in situ observations. It is found that the Noah-MP model can well simulate the SH and LH fluxes. The inter-annual variations of SH fluxes on the central and western TP were larger than those on the eastern TP. All the regions showed a decreasing trend in SH flux from mid-1980s, but an increasing trend since about 2002. The LH fluxes on the eastern TP were larger than those on the central and western TP. But the inter-annual variations were smaller than those of central and western TP. The western and eastern TP got an increasing trend in LH flux during the 38 yrs. On the central TP, LH fluxes exhibited an increasing trend from mid-1980s to the beginning of the 21st century, and then decreased. For the whole TP, LH fluxes increased from mid-1980s to the beginning of the 21st century, and showed a decreasing trend from about 2003. The asynchronous changes of wind speed, Ts-Ta (land-air temperature difference) and water vapor condensation may increase the uncertainties of the trends in SH and LH fluxes.