With the proliferation of nonlinear loads, high neutral harmonic currents in three-phase four-wire distribution system have been observed. It has been known that the ground impedance has an effect on the neutral currents of a system which operates with harmonics present. On-site measurements of harmonic currents and voltages according to the fall-of-potential method under case study system were made and the ground impedance modeling using the pattern search method for the harmonic analysis was developed The ground impedance model obtained by the proposed method was compared with the frequency characteristics by field tests and has shown appropriate results, and would be applicable to evaluate the harmonic and transient response characteristics of the ground system