The canyon surrounding a dam can be assumed as an unbounded domain, and the geometry and flexibility of a canyon are parameters that greatly affect the values of natural periods in earth dams. In this thesis, in order to take into account these two effects, canyons are modeled by SBFEM, and earth dams, which have limited geometries, are modeled by FEM. The hybrid FEM-SBFEM technique used for the dynamic three-dimensional analysis of soil-earth dam interactions is validated with results available in the literature. Because the dynamicstiffness matrix of the unbounded domain is complex and frequency-dependent, the classical modesuperposition method is not straightforward for a soil-structure interaction system, and thus, to obtain their fundamental natural frequencies, the modeled dams were excited in the upstream-downstream direction. The natural periods of earth dams in canyons with different geometries shapes and impedance ratios are obtained, and are found to have significant effects on the dams’ natural periods. The results are compared with actual recorded data, and it is found that the graphs put forward in this study may be used by practical engineers for the estimation of natural periods of earth dams in canyons with different shapes and material properties. Several amplification functions corresponding to different canyon conditions are obtained by applying a uniform displacement at the canyons’ boundaries. A comprehensive study is performed to examine the effects of canyon geometry and flexibility on the steady-state responses of the dams, and it is found that these two effects significantly influence the amplification functions. While the flexibility of the canyon does affect the maximum amplification function value, this value does not change for earth dams in canyons that have different shapes but the same length. In addition, the lateral responses of earth dams in the time domain are computed in order to analyze the aforementioned effects under an actual earthquake. The proposed amplification functions are used to compare the recorded response spectra of the El Infiernillo dam under the two 1966 earthquakes with the calculated amplification function, and a reasonable agreement is observed between them. The equivalent linear method (EQL) is implemented into the FEM, and the FEM-SBFEM technique is extended in order to take into consideration the effect of earth dams’ nonlinear behavior. It is observed that such nonlinear behavior greatly affects the natural frequency, the amplification function, and peak crest acceleration of earth dams located in canyons. The effects of canyon geometry and flexibility on the nonlinear behavior are examined, and it is found that by increasing canyon flexibility, the effect of nonlinearity is decreased. The El Infiernillo dam is modeled by the 3D nonlinear FEM-SBFEM, and comparison of the crest amplification function obtained by the proposed method with the recorded data shows the accuracy of the nonlinear FEM-SBFEM