We investigated the potential capability of the random forest (RF) machine learning (ML) model to estimate snow depth in this work. Four combinations composed of critical predictor variables were used to train the RF model. Then, we utilized three validation datasets from out-of-bag (OOB) samples, a temporal subset, and a spatiotemporal subset to verify the fitted RF algorithms. The results indicated the following: (1) the accuracy of the RF model is greatly influenced by geographic location, elevation, and land cover fractions; (2) however, the redundant predictor variables (if highly correlated) slightly affect the RF model; and (3) the fitted RF algorithms perform better on temporal than spatial scales, with unbiased root-mean-square errors (RMSEs) of ∼4.4 and ∼7.3 cm, respectively. Finally, we used the fitted RF2 algorithm to retrieve a consistent 32-year daily snow depth dataset from 1987 to 2018. This product was evaluated against the independent station observations during the period 1987–2018. The mean unbiased RMSE and bias were 7.1 and −0.05 cm, respectively, indicating better performance than that of the former snow depth dataset (8.4 and −1.20 cm) from the Environmental and Ecological Science Data Center for West China (WESTDC). Although the RF product was superior to the WESTDC dataset, it still underestimated deep snow cover (>20 cm), with biases of −10.4, −8.9, and −34.1 cm for northeast China (NEC), northern Xinjiang (XJ), and the Qinghai–Tibetan Plateau (QTP), respectively. Additionally, the long-term snow depth datasets (station observations, RF estimates, and WESTDC product) were analyzed in terms of temporal and spatial variations over China. On a temporal scale, the ground truth snow depth presented a significant increasing trend from 1987 to 2018, especially in NEC. However, the RF and WESTDC products displayed no significant changing trends except on the QTP. The WESTDC product presented a significant decreasing trend on the QTP, with a correlation coefficient of −0.55, whereas there were no significant trends for ground truth observations and the RF product. For the spatial characteristics, similar trend patterns were observed for RF and WESTDC products over China. These characteristics presented significant decreasing trends in most areas and a significant increasing trend in central NEC.

Accurate long-term daily cloud-gap-filled fractional snow cover products are essential for climate change and snow hydrological studies in the Asian Water Tower (AWT) region, but existing Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover products are not sufficient. In this study, the multiple-endmember spectral mixture analysis algorithm based on automatic endmember extraction (MESMA-AGE) and the multistep spatiotemporal interpolation algorithm (MSTI) are used to produce the MODIS daily cloud-gap-filled fractional snow cover product over the AWT region (AWT MODIS FSC). The AWT MODIS FSC products have a spatial resolution of 0.005° and span from 2000 to 2022. The 2745 scenes of Landsat-8 images are used for the areal-scale accuracy assessment. The fractional snow cover accuracy metrics, including the coefficient of determination (R²), root mean squared error (RMSE) and mean absolute error (MAE), are 0.80, 0.16 and 0.10, respectively. The binarized identification accuracy metrics, including overall accuracy (OA), producer's accuracy (PA) and user's accuracy (UA), are 95.17 %, 97.34 % and 97.59 %, respectively. Snow depth data observed at 175 meteorological stations are used to evaluate accuracy at the point scale, yielding the following accuracy metrics: an OA of 93.26 %, a PA of 84.41 %, a UA of 82.14 % and a Cohen kappa (CK) value of 0.79. Snow depth observations from meteorological stations are also used to assess the fractional snow cover resulting from different weather conditions, with an OA of 95.36 % (88.96 %), a PA of 87.75 % (82.26 %), a UA of 86.86 % (78.86 %) and a CK of 0.84 (0.72) under the MODIS clear-sky observations (spatiotemporal reconstruction based on the MSTI algorithm). The AWT MODIS FSC product can provide quantitative spatial distribution information on snowpacks for mountain hydrological models, land surface models and numerical weather prediction in the Asian Water Tower region. This dataset is freely available from the National Tibetan Plateau Data Center at https://doi.org/10.11888/Cryos.tpdc.272503 (Jiang et al., 2022) or from the Zenodo platform at https://doi.org/10.5281/zenodo.10005826 (Jiang et al., 2023a).

Most climate models show systematic cold biases during snow‐covered period over the Tibetan Plateau (TP), which is associated with snow and surface albedo overestimations. In this work, a snow cover fraction (SCF) scheme and a recently developed albedo scheme for shallow snow are implemented in the Noah‐MP land surface model coupled with the Weather Research and Forecasting (WRF) model. The SCF scheme introduces subgrid orographic variability to reduce the SCF, and the shallow‐snow albedo scheme parameterizes the fresh‐snow albedo as a function of the snow depth (SD). Evaluations by remote sensing data show that both schemes can effectively alleviate the overestimation of the simulated surface albedo, SCF, snow water equivalent, and SD over the TP. The reductions in the modeled SCF and snow albedo directly lead to lower surface albedo values and thus more surface solar radiation absorption, which accelerates snow melting and causes surface warming effects. Further comparisons with Moderate Resolution Imaging Spectroradiometer data and station observations show that both schemes can significantly reduce the cold biases in the surface skin temperature (from −4.39°C to 0.19°C for the TP mean) and 2‐m air temperature (from −4.48°C to −1.05°C for the station mean) during the cold season (October to May of next year) in the study region. This work provides guidance for advancing the snow‐related physics in climate models and the improved WRF model could facilitate weather forecasting and climate prediction for the plateau region.

Microwave emission models are a critical component of snow water equivalent retrieval algorithms applied to passive microwave measurements. Several such emission models exist, but their differences need to be systematically compared. This paper compares the basic theories of two models: the multiple-layer Helsinki University of Technology (HUT) model and the microwave emission model of layered snowpacks (MEMLS). By comparing the mathematical formulation side by side, three major differences were identified: 1) by assuming that the scattered intensity is mostly (96%) in the forward direction, the HUT model simplifies the radiative transfer equation in 4π space into two one-flux equations, whereas MEMLS uses a two-flux theory; 2) the HUT scattering coefficient is much larger than the one of MEMLS; and 3) MEMLS considers the trapped radiation inside snow due to internal reflection by a six-flux model, which is not included in HUT. Simulation experiments indicate that the large scattering coefficient of the HUT model compensates for its large forward scattering ratio to some extent, but the effects of one-flux simplification and the trapped radiation still result in different T _B simulations between the HUT model and MEMLS. The models were compared with observations of natural snow cover at Sodankylä, Finland; Churchill, Canada; and Colorado, USA. No optimization of the snow grain size was performed. It shows that the HUT model tends to underestimate T _B for deep snow. MEMLS with the physically based improved Born approximation performed best among the models, with a bias of -1.4 K and a root-mean-square error of 11.0 K.

Soil freeze-thaw cycles have a profound impact on heat and water fluxes at the land-atmosphere interface and transport in soils. Microwave remote sensing is a widely used technique to detect near-surface soil freeze/thaw states due to significant changes in dielectric properties associated with water phase transitions in soils, where uncertainty remains. This letter proposes a new parameterization scheme for the estimation of unfrozen water content to improve the modeling of soil dielectric properties during freeze-thaw transitions. Predictions from the new model referred to as Zhang-Zhao’s model were compared with dielectric measurements during thawing processes of soil samples collected from Baoding (silty clay soil), Zhangjiakou (loamy sandy soil), and Zhengzhou (clay loam soil) in China. The mean biases of the predictions were 3.25 (4.44 and 2.07 for the thawed value and frozen value, respectively) and 1.54 (2.22 and 0.88 for the thawed value and frozen value, respectively) for the real part and imaginary part, respectively. The model-predicted soil complex relative permittivity (CRP) was highly correlated with measurements, with correlation coefficients ranging from 0.7944 to 0.9865. The normalized root mean square errors of the predictions were 13.72% (real part) and 25.41% (imaginary part).