Acoustic waves from volcanic eruptions are often observed as infrasound in near fields. Part of them propagate upward and disturb the ionosphere, and can be observed with Total Electron Content (TEC) data from Global Navigation Satellite System (GNSS) receivers. Here we report TEC variations after the 13 February 2014 Plinian eruption of the Kelud volcano, East Java, Indonesia, observed with regional GNSS networks. Significant disturbances in TEC were detected with six GNSS satellites, and wavelet analysis showed that harmonic oscillations started at ~16:25 UT and continued for ~2.5 h. The amplitude spectrum of the TEC time series showed peaks at 3.7 mHz, 4.8 mHz and 6.8 mHz. Long-wavelength standing waves with a wide range of wavelength trapped in the lower atmosphere are excited by the Plinian eruption. Amplitude spectra of the ground motion recorded by seismometers, however, had frequency components at discrete wave-periods. The condition for the resonant oscillations between the atmosphere and the solid Earth is satisfied only at these discrete wave-period and horizontal wavelength pairs, therefore efficient energy transfer from the atmospheric standing waves to the solid Earth Rayleigh waves occurred at discrete periods and resulted in the harmonic ground motion.