The Tibetan Plateau (TP) region experiences strong land‐atmosphere interactions, and as an elevated heating source, significantly influences the formation of the Asian monsoon. Those interactions are not well represented in current land‐surface models (LSMs), partly due to difficulties in representing heterogeneities in soil structures in LSM. Simulations using the Noah with multiparameterization options (Noah‐MP) LSM are employed to assess the relative importance of parameterizing vertical soil heterogeneity, organic matter, and soil rhizosphere and their impacts on seasonal evolution of soil temperature, soil moisture, and surface energy and water budgets at the sparsely vegetated Amdo site located in central TP. The LSM spin‐up time at the central TP depends on the complexity of the model physics, ranging from 4 years with simplest soil physics to 30 years with the addition of organic matter and spare to dense rhizosphere parameterization in Noah‐MP. Representing layered soil texture and organic matter does not improve low biases in topsoil moisture. Reducing the saturated conductivity from the mucilage in the rhizosphere produces better results. Surface sensible and latent heat fluxes are better simulated in the monsoon season as well. Adding layered soil texture and organic matter in Noah‐MP retard the thawing in deep soil layers, and the rhizosphere effect delays thawing even more in the transient season. Uncertainties in soil initialization significantly affect deep‐soil temperature and moisture, but uncertainties in atmospheric forcing conditions mainly affect topsoil variables and consequently the surface energy fluxes. Differing land‐surface physics cause 36% uncertainty in the accumulated evapotranspiration and subsurface runoff.

Abstract Precipitation is a critical input to land surface and hydrology modeling and prediction. Dynamical downscale modeling has added value to representing precipitation, when compared with the performance of coarse-resolution reanalysis and global climate models, over the Tibetan Plateau (TP). Convection-permitting modeling (CPM) may even outperform dynamical downscale models (DDMs). In this study, 4-km CPM results were compared to 28-km DDM results for a snow season (1 October 2013–31 May 2014) over the TP. The CPM- and DDM-simulated precipitation, as well as three merged gridded precipitation datasets, were evaluated against in situ observations below 4800 m. The five precipitation datasets (CPM, DDM, CMFD, COPRPH, and TRMM) showed large differences over the TP with underestimation of TRMM and overestimation of CPM and DDM compared to observations. The most significant difference occurred in the Brahmaputra Grand Canyon. Given the substantial uncertainty in observed precipitation at high mountains, snow cover simulated by a high-resolution land data assimilation system was used to indirectly evaluate the above precipitation data using MODIS observations. Simulated snow-cover fraction was greatly underestimated using all the merged precipitation datasets. However, simulations using the DDM- and CPM-generated precipitation as input outperformed those using other gridded precipitation data, showing lower biases, higher pattern correlations, and closer probability distribution functions than runs driven by the merged precipitation. The findings of this study generally support the assumption that high-resolution CPM-produced precipitation has an added value for use in land surface and hydrology simulations in high-mountain regions without reliable in situ precipitation observations.

Soil hydrothermal regime of the active layer in the permafrost regions of the Qinghai‐Tibet Plateau (QTP) is important to the underlying permafrost and the climate change dynamics in Asia. However, a large bias still exists in current land surface models in the representation of soil temperature and moisture. This study assessed and augmented the Noah land surface model with multiparameterization options (Noah‐MP) for simulating soil hydrothermal dynamics at the Tanggula (alpine meadow) and Beiluhe (alpine swamp) stations located in the permafrost regions of the QTP. The results showed that the default Noah‐MP tended to underestimate soil temperature and moisture. Specifically, the default model overestimated the snow depth and duration due to the low snow sublimation rate. This resulted in a cold deviation in the soil temperature at two stations. Such underestimation was reduced by introducing a scheme that considered the sublimation loss from wind. Moreover, the remaining cold bias in the soil profiles of two stations was greatly resolved by a combined scheme of roughness length for heat (Z0h) and undercanopy aerodynamic resistance (ra,g). A soil thermal conductivity scheme, which can produce more realistic soil thermal conductivity in frozen soil, further improved the deep soil temperature simulation. The consideration of soil organic matter could mitigate the underestimation of the shallow soil moisture to some extent, but this improvement was more obvious at the Tanggula station, which had coarser mineral soil than the Beiluhe station.