Global Navigation Satellite System Interferometric Reflectometry (GNSS-IR) technology provides a new means of snow depth detection. Multi-satellite and multi-Signal-to-Noise Ratio (SNR) provide more data for daily high-precision snow depth retrieval, but also face the problem of data fusion and effective utilization. Therefore, this study proposes a robust estimation algorithm based on multi-satellite and multi-SNR fusion applied to the observations of a GNSS station in Alaska. This study uses four solutions (Savg, Smed, SRE_avg and SRE_med) to carry out multi-system fusion snow depth inversion and precision comparison research. The Savg has more obvious disadvantages, which is not suitable for snow depth assessment. The SRE_avg and SRE_med have better snow depth retrieval effects in the snowy periods. The correlation coefficient (R), root mean square error (RMSE) and mean error (ME) of the calculated snow depth using the robust estimation algorithm with respect to the nearest in-situ measurements reached 0.759, 3.7 cm and −1.4 cm, respectively. Compared with the Smed, the R is increased by 2.0 %, the RMSE corresponds to an improvement of 2.6 %. Moreover, the ME of the snow depth retrievals, as an indicator of the measurement bias, has significantly decreased by 6.7 %. The result also shows that the snow depth inversion by the robust estimation algorithm is more consistent with the in-situ measurements, further extending and advancing the optimal algorithm for snow depth retrieval.