Time evolution of Marangoni flow instability of molten silicon has been investigated at high Marangoni number region by measuring temperature oscillation using thermocouples and surface oscillation using phase-shift interferometry with a combination of a wavelet transformation. The relationship between azimuthal wave number m and aspect ratio F was experimentally clarified. Through surface oscillation observation, instability with various frequency was found to coexist.