Which global reanalysis dataset represents better in snow cover on the Tibetan Plateau?

The extensive snow cover across the Tibetan Plateau (TP) regions has a major influence on the climate and water supply for over one billion downstream inhabitants. However, an adequate evaluation of snow cover fraction (SCF) variability over the TP simulated by global multiple reanalysis datasets has yet to be undertaken. In this study, we examined eight global reanalysis SCF datasets using the Snow Property Inversion from Remote Sensing (SPIReS) product spanning the period 2001–2020. The results reveal that the HMASR generated the best SCF simulations because of its outstanding spatial and temporal accuracy. The GLDAS and CFSR demonstrated acceptable SCF accuracy with respect to spatial variability, but struggled to reproduce the annual trend. Pronounced SCF overestimations were found when using the ERA5, ERA5L, and JRA55, but SCF was underestimated by MERRA2, and CRAL generated poor spatial pattern. Overall biases were related to the combined effect of precipitation forcing, temperature forcing, snow data assimilation, and SCF parameterization methods, with the dominant factor changing across datasets. In ERA5 and ERA5L, temperature and snowfall bias exhibited significant correlations with SCF bias over most TP areas, therefore having a greater impact on the accuracy of SCF in terms of spatial variability and temporal evolution. On the other hand, the impact of snow assimilation was possibly more pronounced in MERRA2 and CRAL. Although parameterization methods can improve SCF simulation accuracy, their influence was weaker than those of other factors, except for JRA55. To further improve the accuracy of SCF simulation, an ensemble average method was developed. The ensemble based on HMASR and GLDAS generated the most accurate SCF spatial distribution, whereas the ensemble containing ERA5L, CFSR, CRAL, GLDAS, ERA5, and MERRA2 proved optimal for capturing the annual trend. [ABSTRACT FROM AUTHOR]