Earthquakes pose significant natural hazards and impact populations worldwide. Iran is among the most susceptible countries to seismic activity, making comprehensive earthquake risk assessment crucial. This study employs geospatial methods, including integrating satellite, ground-based, and auxiliary data to model earthquake risk across this country. A Fuzzy Inference System (FIS) is used to generate earthquake hazard probability and vulnerability layers, considering factors such as slope, elevation, fault density, building density, proximity to main roads, proximity to buildings, population density, and earthquake epicenter, magnitude, proximity to the epicenter, depth density, peak ground acceleration (PGA). The results highlight high-risk areas in the Alborz and Zagros Mountain ranges and coastal regions. Moreover, the findings indicate that 39.7% (approximately 31.7 million people) of Iran’s population resides in high-risk zones, with 9.6% (approximately 7.7 million) located in coastal areas vulnerable to earthquakes. These findings offer valuable insights for crisis management and urban planning initiatives.