Snowfall and the subsequent evolution of the snowpack play important roles in the cryospheric and hydrospheric processes which occur on the Tibetan Plateau (TP). Current literature provides scarce evidence covering the sensitivity of solid precipitation to land surface physics schemes and initial and boundary conditions on the TP. Six numerical experiments using the Weather Research and Forecasting (WRF) model were conducted to simulate a snow event over the TP in March, 2017. Different land surface physics schemes, i.e. CLM, Noah and Noah‐MP, and initial and boundary conditions provided by atmospheric reanalysis datasets, i.e. the NCEP‐FNL and ERA‐Interim datasets, were applied in sensitivity analyses. The observed near‐surface air temperature, snow depth (SD) and snow water equivalent (SWE) values were used to evaluate each model's performance. The results demonstrate that: (1) the sensitivity of the near‐surface air temperature to land surface physics schemes is greater than it is to both the initial and boundary conditions; (2) the best performance is achieved when applying WRF + CLM with a root mean square error of 8.4°C, a mean absolute deviation of 7.3°C, a correlation coefficient of 0.75 and a spatial correlation coefficient of ~0.5 to air temperature estimates. A potentially important factor appears to be the advanced parametrization of albedo in the CLM scheme; (3) the advanced land surface schemes in the WRF model describes the physics of cryospheric and hydrospheric processes in detail, and the land surface response is determined by multiple variables and parameters in such schemes. The spatial patterns in such variables and parameters determined the detailed spatial variabilities observed in snow cover and amount and its temporal evolution. The WRF model overestimates, however, the intensity and extent of SD and SWE; (4) simulations of solid precipitation are more accurate when applying CLM or Noah‐MP + ERA‐Interim in WRF; and (5) WRF performance with regard to SWE estimates clearly depends upon the discrimination of lighter from heavier snowfall.

Snow albedo is important to the land surface energy balance and to the water cycle. During snowfall and subsequent snowmelt, snow albedo is usually parameterized as functions of snow-related variables in land surface models. However, the default snow albedo scheme in the widely used Noah land surface model shows evident shortcomings in land–atmosphere interaction estimates during snow events on the Tibetan Plateau. Here, we demonstrate that our improved snow albedo scheme performs well after including snow depth as an additional factor. By coupling the Weather Research and Forecasting (WRF) and Noah models, this study comprehensively evaluates the performance of the improved snow albedo scheme in simulating eight snow events on the Tibetan Plateau. The modeling results are compared with WRF run with the default Noah scheme and in situ observations. The improved snow albedo scheme significantly outperforms the default Noah scheme in relation to air temperature, albedo and sensible heat flux estimates by alleviating cold bias estimates, albedo overestimates and sensible heat flux underestimates, respectively. This in turn contributes to more accurate reproductions of snow event evolution. The averaged root mean square error (RMSE) relative reductions (and relative increase in correlation coefficients) for air temperature, albedo, sensible heat flux and snow depth reach 27 % (5 %), 32 % (69 %), 13 % (17 %) and 21 % (108 %), respectively. These results demonstrate the strong potential of our improved snow albedo parameterization scheme for snow event simulations on the Tibetan Plateau. Our study provides a theoretical reference for researchers committed to further improving the snow albedo parameterization scheme.

Meteorological variables (e.g., air temperature (T2), radiation flux, and precipitation) determine the evolution of glacier mass and characteristics. Observations of these variables are not available with adequate spatial coverage and spatiotemporal resolution over the Tibetan Plateau. Albedo is the key factor of net radiation and is determined by the land cover and snow-related variables. This study focuses on evaluating the performance of the albedo parameterization scheme in WRF coupled with Noah-MP in terms of glacio-meteorological variables, by conducting experiments applying the standard surface albedo scheme with the default vegetation and corrected to ice cover and the modified glacial albedo scheme to the Parlung No. 4 Glacier in the 2016 ablation season. In situ glacio-meteorological element observations and MODIS-retrieved albedo are selected to assess the performance of the model. The key results are as follows. First, compared to the air temperature bias of 1.56 °C in WRF applying the standard surface albedo scheme and the default vegetation cover, realistic land-use categories considerably reduce the model warm bias on the glacier. The model using realistic land-use categories yields similar T2 diurnal patterns to the observations, with a mean bias of only −0.5 °C, no matter which glacial albedo scheme is implemented. Second, the default glacial albedo scheme gives a rather high albedo value of 0.68, causing an apparent underestimation of the net shortwave radiation and net radiation; the modified glacial albedo scheme gives a mean albedo value of 0.35, close to the in situ observations, helping to relieve underestimations of net shortwave radiation and net radiation. Compared with the MODIS albedo of the glacier, WRF applying the default glacial albedo scheme apparently overestimates the albedo with a mean error of 0.18, while WRF applying the modified glacial albedo scheme slightly underestimates the albedo with a mean error of only −0.08. Third, the mean net radiation flux (142 W m⁻²) and high ground heat flux (182 W m⁻²) values that were estimated by WRF applying the corrected land cover and the modified glacial albedo scheme result in the heating of the glacier surface and subsurface, causing ice melt and the liquid water content to increase more quickly and preferentially, equating to an estimated ice thickness decrease of 1 m by mid-June in the ablation region. Our study confirms the ability of the WRF model to reproduce glacio-meteorological variables as long as a reasonable glacial albedo scheme and the corrected land cover is applied and provides a theoretical reference for researchers that are committed to further improvement of the glacial albedo scheme.

The snow albedo is a vital component of land–atmosphere coupling models. It plays a critical role in regulating land surface energy exchange by controlling incoming solar radiation absorbed by the land surface and influencing the timing and rate of snowmelt. Accurate snow albedo simulation is essential to obtain surface energy balance and snow-cover estimates. Here, the simulation of albedo and snow cover using the Weather Research and Forecasting Model and an improved snow albedo scheme is verified against satellite-retrieved products during and immediately following eight snowfall events over the Tibetan Plateau. The improved model successfully characterizes the spatial pattern and inverted U-shaped temporal pattern of albedo over the entire Tibetan Plateau. This is attributed to the local optimization of snow-age parameters and explicit consideration of snow depth in the improved scheme. Compared with the previous model, the model proposed herein greatly decreases the overestimated albedo (by 0.13–0.27), yielding a bias range of ±0.08, mean relative bias decrease of 70%, and significant increase in the spatial correlation coefficient of 0.03–0.39 (mean: 0.13). The significant improvements of albedo estimates appear in deep snow-covered regions, largely attributed to parameter optimization related to snow albedo decay, while less improvements appear over the shallow snow-covered regions. Accurate reproduction of the spatiotemporal variation in albedo alleviated snow-cover overestimation by small amounts. For snow-cover estimates, the improved model consistently decreases the false-alarm rate by 0.03, and increases the overall accuracy and equitable threat score by 0.04 and 0.03, respectively. Moreover, the improved scheme shows an equivalent improvement of albedo estimates at both 1- and 5-km grid spacing over the eastern Tibetan Plateau; this is also true for snow-cover estimates. Significance Statement: Snow albedo schemes in widely used numerical weather prediction models show notable shortcomings in complex mountainous regions, hindering accurate surface energy balance and snow-cover prediction. The purpose of this study is to better understand the role of snow albedo on snow-cover estimates and reveal the application potential of an improved snow albedo scheme across the Tibetan Plateau. This is important because snow albedo influences the timing and rate of snowmelt, and in turn snow-cover estimates, through regulating the surface energy budget. Our results highlight the strong application potential of our improved scheme in reducing snow simulation errors, confirm the importance of snow depth on snow albedo, and provide a new perspective for improving the accuracy of snow forecast over the topographically high Tibetan Plateau.

Meteorological variables (e.g. air temperature (T2), radiation flux and precipitation) determine the evolution of glacier mass and characteristics. Observations of these variables are not available at adequate spatial coverage and temporal / spatial resolution over the Tibetan Plateau. This study focuses on evaluating the performance of albedo parameterization scheme in WRF coupled with Noah-MP in terms of glacio-meteorological variables, by conducting sensitivity experiments on Parlung No. 4 Glacier in ablation season in 2016. The control experiment (CTL) uses the model’s default albedo scheme and unrealistic surface type, while sensitivity experiments apply the default albedo scheme (Sens1) and the modified scheme (Sens2) in ice surface. The key results are as follows. First, all experiments yield similar T2 diurnal patterns to the observations and realistic land-use type considerably reduces model warm bias on the glacier. The RMSE of T2 decreases by 1 °C with an improvement of 37 % by sensitivity experiment estimates. Second, Sens1 keeps rather high albedo value of 0.68 causing net shortwave radiation and net radiation apparent underestimation, while Sens2 holds mean albedo value of 0.35 close to observations contributing to relieve net shortwave radiation and net radiation underestimation. Thirdly, compared with Sens1 estimates, more energy from Sens2 is used to heat ice surface resulting in high ground heat flux (182 W m-2), ice melt and liquid water content increase more quickly and preferentially. Our study confirms the ability of WRF model in reproducing glacio-meteorological variables as long as a reasonable albedo parameterization scheme is applied, and provides a theoretical reference for researchers committed to further improving the surface albedo scheme.