Dynamic properties of the soils at the Lotung test site, Lotung, Taiwan, are estimated from seismic vertical array measurements (input-output data sets) using both time-invariant and time-variant parametric modeling methods (system identification). Soil properties are directly mapped from model parameters to an equivalent lumped mass model of the soil interval. Shear stiffness and damping ratios were calculated for 8 events with M L ranging from 4.5 to 7.0. Shear stiffness ranged between 0.5 and 6MN/m, inversely proportional to PGA. The equivalent viscous damping ratio varied from 2 to 30% of critical damping, proportional to PGA. Degradation of soil behavior, while less pronounced with increasing depth, consistently occurs above a peak input acceleration of 0.07g. Although ''non-linear'' behavior is evident above 0.17g, Event 7 (0.21g) is accurately predicted using a linear constant parameter model estimated from the smaller Event 8 aftershock ground motions.