Thunderstorm is one of the most severe atmospheric disturbances in the world and also in Iran, which is characterized by rapid upward movements, abundant moisture, and climatic instability. Since this phenomenon is usually accompanied with hail, lightning, heavy rain, flood and severe winds, it can cause irreparable damage to the environment. Investigation of spring thunderstorms has a great significance regarding the irreparable damages can cause by them and also because of the higher frequency of this phenomenon in the spring and the necessity for preparedness and disaster mitigation actions. To identify the locations of the major thunderstorm risk areas, the entire country with an area of 1648195 square kilometers, which is located between the 25°-40° north latitude and 44°-63° east longitude is considered. Spatial distribution of the occurrence of hazardous spring thunderstorms was analyzed using a series of monthly thunderstorm frequency data obtained from 25 synoptic stations over a 51-year-long period (1960-2010). Ward's hierarchical clustering and Kriging methods were used for statistical analysis. Initially, total number of thunderstorms in April, May and June were considered as the frequency of occurrence of thunderstorm in different stations in the spring. Measure of central tendency and dispersion which consists of the sum, minimum, maximum, range and coefficient of variation, standard deviation, and skewness were used to clarify the changes of thunderstorms and to determine the spatial and temporal climatic distribution of spring thunderstorms. An appropriate probability distribution function was chosen to determine the distributions of the data. Due to the large volume of data and the uneven distribution of stations, cluster analysis and kriging methods were used to classify different regions into homogeneous groups for zoning and spatial analysis of spring thunderstorms, respectively. The statistical characteristics of spring thunderstorms were reviewed and fitted with a 3-parameter Weibull distribution. Regions considered for this study were classified in four separate clusters according to the simultaneity of thunderstorms in the spring. After zoning, it was found that the highest rates of thunderstorm took place in the northwest and west of country. The northeast of Iran has the second highest number of thunderstorm occurrence. The least number of thunderstorm event had happened in the central and southern half of the country. According to the descriptive statistics parameters, maximum number of thunderstorms occurred in May.. Based on the results of the cluster analysis, there is a similar trend in the central and eastern regions, the rest of the country was clustered into five distinct homogeneous regions, including the northwestern, western, southern, northern, central northern and northeastern regions. Zoning results indicate that the highest number of the occurrence of this phenomenon in the country is concentrated in the northwestern and western regions. Higher frequency of occurrence of thunderstorms in the northwestern and western regions may be attributed to local topographic conditions like high mountains, orientation of the terrain, solar radiation on slopes and existence instability conditions, hillside convection, the presence of water resources and specific climatic conditions in these areas. In addition, as a result of a continuous surface obtained by the method of interpolation with the least amount of systematic error and also the use of correlation functions for recognizing the spatial structure of the data and estimating the model error when using the Kriging method, the weights are chosen in order to have a more optimized interpolation function. Also the cluster analysis may significantly reduce the volume of operation without affecting the results and will help in finding a real band due to more appropriate classification of different geographic areas with greater spatial homogeneity and minimal variance within the group. Based on the results of the spatial analysis, it is clear that Kriging and Ward cluster analysis methods are appropriate for thunderstorm zoning and classification of different regions according to occurrence of thunderstorm, respectively.