This paper introduces system-parameter modelling, commonly known as system identi®cation, to the estimation of dynamic soil properties. These methods are optimized for short and nonstationary signals. The Wildlife Site in California, subject to two large earthquakes (Elmore Ranch and Superstition Hills) on 24 November 1987, was used in an example analysis as the associated data are the only publicly available record of buried and surface motions. During the second, larger, earthquake the site soils lique®ed. For the ®rst time, both stationary and recursive system-parameter modelling methods were applied to the unique input±output data from the two earthquakes. The results show weakening of the soils system as pore pressure increases. Estimates of natural frequency, damping ratio and participation factor are given.