The modeling of climate over the Tibetan Plateau (TP) started with the introduction of Global Climate Models (GCMs) in the 1950s. Since then, GCMs have been developed to simulate atmospheric dynamics and eventually the climate system. As the highest and widest international plateau, the strong orographic forcing caused by the TP and its impact on general circulation rather than regional climate was initially the focus. Later, with growing awareness of the incapability of GCMs to depict regional or local-scale atmospheric processes over the heterogeneous ground, coupled with the importance of this information for local decision-making, regional climate models (RCMs) were established in the 1970s. Dynamic and thermodynamic influences of the TP on the East and South Asia summer monsoon have since been widely investigated by model. Besides the heterogeneity in topography, impacts of land cover heterogeneity and change on regional climate were widely modeled through sensitivity experiments. In recent decades, the TP has experienced a greater warming than the global average and those for similar latitudes. GCMs project a global pattern where the wet gets wetter and the dry gets drier. The climate regime over the TP covers the extreme arid regions from the northwest to the semi-humid region in the southeast. The increased warming over the TP compared to the global average raises a number of questions. What are the regional dryness/wetness changes over the TP? What is the mechanism of the responses of regional changes to global warming? To answer these questions, several dynamical downscaling models (DDMs) using RCMs focusing on the TP have recently been conducted and high-resolution data sets generated. All DDM studies demonstrated that this process-based approach, despite its limitations, can improve understandings of the processes that lead to precipitation on the TP. Observation and global land data assimilation systems both present more wetting in the northwestern arid/semi-arid regions than the southeastern humid/semi-humid regions. The DDM was found to better capture the observed elevation dependent warming over the TP. In addition, the long-term high-resolution climate simulation was found to better capture the spatial pattern of precipitation and P-E (precipitation minus evapotranspiration) changes than the best available global reanalysis. This facilitates new and substantial findings regarding the role of dynamical, thermodynamics, and transient eddies in P-E changes reflected in observed changes in major river basins fed by runoff from the TP. The DDM was found to add value regarding snowfall retrieval, precipitation frequency, and orographic precipitation. Although these advantages in the DDM over the TP are evidenced, there are unavoidable facts to be aware of. Firstly, there are still many discrepancies that exist in the up-to-date models. Any uncertainty in the model’s physics or in the land information from remote sensing and the forcing could result in uncertainties in simulation results. Secondly, the question remains of what is the appropriate resolution for resolving the TP’s heterogeneity. Thirdly, it is a challenge to include human activities in the climate models, although this is deemed necessary for future earth science. All-embracing further efforts are expected to improve regional climate models over the TP.