Sulawesi Island is the active tectonic region, where the tectonic architecture and potential earthquake sources until now remain largely unknown. The worst earthquake, an Mw 7.5 on September 28, 2018, in Palu, Indonesia, was caused catastrophic damage to life and property. The earthquake has highlighted the urgent need to raise knowledge of the cause of possible large future earthquakes and vulnerability. The main objective for this project is to create a thorough earthquake probabilistic hazard analysis map of the region, which is presently unavailable to better prepare for future earthquakes. The neotectonic and structural map was created using was supplemented with the 30-m resolution Shuttle Radar Topography Mission, Centroid Moment Tensor (CMT) solution, and seismologic data. The results demonstrate that faulting controls the geometry and the majority of these faults are active and capable of causing medium to large magnitude earthquakes with moment magnitudes ranging from 6.2 to 7.5 from 44 seismic sources. Our results show Sulawesi's northern deformation regimes have high seismicity risk and vulnerability. This study contributes a realistic seismic source for the Sulawesi neotectonic area particularly at the northwest, north, and east deformation regime, to understand the key large future earthquakes.