The Wasatch fault zone defines the eastern boundary of the actively extending Basin and Range Province (Utah, western United States) and poses a significant seismic hazard to the metropolitan areas along the Wasatch Range. A wealth of paleoseismological data documents ∼24 surface-rupturing Mw ≥ 7 earthquakes along the Wasatch fault during the past 6400 yr. Here, we simulated the Holocene earthquake sequence on the Wasatch, Oquirrh–Great Salt Lake, and West Valley faults using three-dimensional finite-element forward modeling with the goal to calculate coseismic and postseismic Coulomb stress changes and to evaluate the slip and magnitude of hypothetical present-day and future earthquakes. Our results show that a good fit between modeled and observed paleoevents and time-integrated slip rates can be achieved within the uncertainties of the paleoseismological record and model parameters like the fault geometry. The Coulomb stress change analysis for selected paleoearthquakes showed that maximum positive stress changes are induced on faults located along strike of the source fault, while faults parallel to the source fault are generally located in stress shadow zones. Postseismic viscoelastic relaxation considerably modifies the coseismic stress changes; the resulting transient stress changes are recognizable for more than 100 yr after an earthquake. The modeled present-day state of Coulomb stress changes shows that the Brigham City, Salt Lake City, and Provo segments of the Wasatch fault are prone to failure in a Mw ≥ 6.8 earthquake. Our study shows that simulation of an entire earthquake sequence based on a paleoseismological record is feasible and facilitates identification of possible gaps and inconsistencies in the paleoseismological record. Therefore, forward modeling of earthquake sequences may ultimately contribute to improved seismic hazard estimates.