Complex computer models of real-world processes (or simulators) are an essential ingredient to carry out uncertainty quantification in virtually every field of science and engineer. In coastal emergency risks assessment, storm surge is one of the most severe natural disasters, and it can lead to significant flooding in coastal areas and severe damages to the life and property from a hurricane. Quantifying the risk to storm surge hazard requires large-scale numerical simulations of hurricanes from storm surge modeling systems. The Advanced Circulation model is one of the primary storm surge models used to predict storm surge and control the impact of storm damage in the United States. A crucial need is the development of an emulator to facilitate risk quantification. However, the model outputs are high-dimensional zero-inflated spatial data, and Gaussian processes are unable to handle such data directly. To tackle both modeling and computational challenges, we take a hierarchical modeling approach to develop a zero-inflated emulator. This methodology is applied to storm surge over coastal areas of the United States and coastal emergency risks assessment.