On 14 December 2023, a heavy snowstorm occurred in the areas including southern Shanxi, southern Hebei, northern Henan, and western Shandong, accompanied by frequent lightning flashes. The ECMWF (European Centre for Medium-Range Weather Forecasts) Reanalysis v5 (ERA5) and the data of lightning, atmospheric electric field, dual-polarization radars are used to analyze the thundersnow process in detail, and the following results are obtained. (1) The thundersnow event is an obvious snowstorm system with elevated convection. In the snow area, there has been a temperature inversion layer at the lower troposphere, and significant warm and wet advection exists between 850 and 600 hPa. The northeasterly below 925 hPa forms a “cold pad” wedging into the conditional symmetric instable layer, where the warm and wet air is lifted to trigger convection and release instable energy, consequently producing relatively strong updraft. (2) The lightning activities during the snowfall are relatively active with the peak flash frequency of 99 times·(10 min)^-1. The ratio of intracloud flashes to cloud-to-ground lightning activities is 1.7∶1 and the positive cloud-to-ground lightning activities account for 19.4%. The spatial distribution and trend of total lightning flashes in the thundersnow event show a good consistency with the snowfall area, and the area with dense lightning activities corresponds exactly to the heavy snowfall area. Lightning flashes are reported within 30 km around 74% of the weather stations with heavy snow (snowfall amount p≥10 mm), while lightning flashes occur within 30 km around 100% of the stations with heavy blizzard (p≥20 mm). (3) The positive K_DP area aloft is closely related to the increase of surface precipitation rate. The spatial consistency between total lightning activities and radar reflectivity at high altitudes such as 6 km is very well because the total lightning activities are basically distributed in the cloud area with reflectivity larger than 20 dBZ, and the echo top height of the lightning falling area is almost above 5 km. It is found that the electrification conditions are good in the cloud where the presence of obviously strong updraft and supercooled liquid water is revealed by positive Z_DR and K_DP in the aloft layer of 5-8 km. Combined with the observations of surface electric field and lightning information, it is speculated that the charge structure of the cloud system during the snowfall may be a normal dipole pattern with positive above negative charge region.