This paper presents a statistical analysis of strong motion data recorded on existing embankment dams. The study aims to enhance our understanding of the seismic response of these structures, providing valuable insights for the design of new dams and the seismic evaluation of existing ones. Utilizing the JCOLD 2014 strong motion database, we examined 190 shaking events from 54 earthquakes recorded on 60 dams. The database, spanning from 1978 to 2012, includes tri-directional strong motion records at both dam foundations and crests. Multivariate analysis and linear/nonlinear fit modeling were employed to investigate seismic response characteristics and their pairwise correlations. The focus was on seismic amplification and fundamental periods, and their relationships with shaking intensity, dam geometry, and input motion frequency content. Nonlinear fit prediction equations and curves were developed, incorporating uncertainty (RMSE) and R-square values in three directions. Despite significant data scatter, these equations offer practical utility in estimating seismic response during preliminary design or seismic assessment phases.