A three-dimensional geodesic finite difference time domain (G-FDTD) algorithm is employed to model the diurnal/seasonal variations of the first five Schumann resonance parameters (SRs), including resonance frequencies and quality factors, in the Earth-ionosphere (E-I) system. In order to make the results more accurate, Prony's method is introduces to obtain these parameters. Through several simulations the SRs shift with time/locations are found and compared with reference results. Besides, different conductivity profiles of the ionosphere including the two-exponential profile and the “knee” profile are also considered in the simulations to assure the results.