Frequency-dependent soil-foundation stiffness values are identified with a new method. The method requires only a few (spatially sparse) response signals. The unknown input signal is non-stationary (an earthquake). Timoshenko beam represents the superstructure for computational efficiency.