Models for estimating earthquake ground motions are a key component in seismic hazard analysis. In data-rich regions, these models are mostly empirical, relying on the everincreasing ground-motion databases. However, in areas in which strong-motion data are scarce, other approaches for ground-motion estimates are sought, including, but not limited to, the use of simulations to replace empirical data. In Israel, despite a clear seismic hazard posed by the active plate boundary on its eastern border, the instrumental record is sparse and poor, leading to the use of global models for hazard estimation in the building code and all other engineering applications. In this study, we develop a suite of alternative ground-motion models for Israel, based on an empirical database from Israel as well as on four data-calibrated synthetic databases. Two host models are used to constrain model behavior, such that the epistemic uncertainty is captured and characterized. Despite the lack of empirical data at large magnitudes and short distances, constraints based on the host models or on the physical grounds provided by simulations ensure these models are appropriate for engineering applications. The models presented herein are cast in terms of the Fourier amplitude spectra, which is a linear, physical representation of ground motions. The models are suitable for shallow crustal earthquakes; they include an estimate of the median and the aleatory variability, and are applicable in the magnitude range of 3– 8 and distance range of 1–300 km.