In a volcanic region, unstable conditions of the surface geology result in natural hazards such as landslides. The present study aims to examine the properties of the formation near the surface in the Tatun Volcano Group (TVG), Northern Taiwan. Dense Very Low Frequency magnetotelluric (VLF-MT) surveys were conducted to obtain die near-surface resistivity. A geostatistical approach and stochastic simulation were used to examine the continuity of the data, and to render the spatial distribution of the resistivity. The results show that the estimated range is consistent with the length scale of the geological features such as the hydrothermal phenomena and volcanic edifices. The estimated distribution of the resistivity successfully identified a fresh lava (>= 1000 Omega m), a lineament-bearing, weathered lava (<= a few hundred Omega m), and hydrothermal areas (<= 30 Omega m). On the basis of the present study, we expect further detailed geophysical investigations will be performed from the viewpoint of hazard mitigation in the TVG.