The estimation of snow hazard and load faces the small sample size effect because of the short snow depth record at a station. To reduce such an effect, we propose to estimate the return period value of the annual maximum ground snow depth S, s_T, for Canada sites by applying the regional frequency analysis (RFA) and the region of influence approach (ROIA). The use of RFA and ROIA to map Canadian snow hazard is new. The comparison of their performance for snow hazard mapping has not been explored in the literature. We also consider the at-site analysis approach (ASA) for estimating s_T by using three often used probability distributions for S. A comparison of the estimated s_T by using the three approaches (ASA, RFA, ROIA) indicates that there is considerable scatter between the estimated s_T value although the identified overall spatial trends of s_T are similar. It is shown that the two-parameter lognormal distribution for S at most Canadian sites, based on the at-site analysis, is preferred; this differs from the Gumbel distribution used to develop the design snow load in Canadian structural design code. The new findings indicate that it is valuable to consider the lognormal distribution for developing design snow load for Canadian sites.