The subduction interface in northern Chile hosts a large tsunamigenic potential that may trigger devastating tsunamis in the near future. To quantify tsunami hazard, mitigate economic losses, and improve planning and safety protocols for government entities, it is key to understand the behavior of the seismic and tsunami source. In this work, we used a stochastic approach to generate (slip) scenarios in a nonregular domain. All of them satisfy physical laws and show a moment magnitude defined by historical seismicity and recent Global Positioning System records based on seismic coupling. Each scenario sets an initial condition to be propagated, then flow depth and runup are computed in a high‐resolution grid. The junction of these data allows the creation of probability maps for several parameters, as runup, flow depth, and flux moments. This methodology is also valid to obtain the maximum inundation zone, which is compared to the official map of the Chilean government. The results suggest that the studied region has a notable tsunami hazard, which is reflected in the obtained inundation zone.