Snow disaster is one of the meteorological disasters with a wide range of impact in winter. The technical parameters for engineering snow disaster prevention include the calculation of snow density, snow pressure, and other snow accumulation parameters. Due to the scarcity of snow observation stations in southern provinces and the lack of data, the calculation of important parameters for snow disaster prevention often uses data from other meteorological stations with snow accumulation observations as a substitute. How to choose representative stations with scientific significance is an urgent problem to be solved. The snowfall climate zoning based on similar climate backgrounds can provide a scientific basis for the selection of representative stations for snow cover parameters in areas without data. In the study of snowfall climate zoning in Hubei Province, due to the significant interannual fluctuations of climate indicators such as the first and last days of snowfall, snowfall amount, we use 12 climate indicators such as the dates of the first and last days of snowfall, the number of snow days, the number of snow cover days, the amount of snowfall and the maximum snow depth, to reflect the fluctuation characteristics of climate indicators in the snowfall area of Hubei Province. At the same time, we draw inspiration from common research methods on snowfall climate in the northern region and adopt the clustering analysis method based on Hilbert Schmidt Independence Criterion (HSIC) kernel function to calculate the overall similarity of multidimensional time series indicators to carry out the classification and zoning of snowfall climate in Hubei Province. The results show that the snowfall climate in Hubei Province can be divided into four climatic zones: southeast, central, northwest, and southwest. The zonal distribution characteristics of the zones are consistent with the climatic background of heavy snowfall in Hubei Province caused by the cold air from the north. The first snow day is delayed from the northwest to the central, southwest, and southeast, and the last snow day is just the opposite. The number of snow days and the number of days with snow cover in the northwest are the greatest. The southeast representative station is Huangshi station; the central representative station is Macheng, Wuhan, and Zhongxiang; the southwest representative station is Xianfeng and Badong; and the northwest representative station is Yunxi and Laohekou. The HSIC kernel function can handle the similarity between sets of indicator sequences with significant interannual fluctuations well, and its clustering method is more reasonable for the climate zoning of Hubei snowfall. The zoning results provide a technical basis for the refined snow disaster prevention in Hubei Province.