This paper develops a two-stage stochastic program and solution procedure to optimize the selection of seismic retrofit strategies to increase the resilience of electric power systems to earthquake hazards. The model explicitly considers the range of earthquake events that are possible and, for each, an approximation to the distribution of damage that is experienced. This is important because electric power systems are spatially distributed; hence their performance is driven by the distribution of damage of the components. We test this solution procedure against the nonlinear integer solver in LINGO 13 and apply the formulation and solution strategy to the Eastern Interconnection where the seismic hazard stems from the New Madrid Seismic Zone.