<![CDATA[Aceh Province, located in the Sumatra megathrust zone of Indonesia, is one of the most seismically active regions in Southeast Asia. Understanding the spatial distribution of earthquake magnitudes is essential for disaster mitigation and risk management. This study compares two spatial interpolation methods Inverse Distance Weighted (IDW) and Kriging to determine the most accurate approach for mapping earthquake intensity in Aceh Province. A total of 2,255 earthquake events with magnitudes of 2.5 M and above, recorded between 1990 and 2024 by the United States Geological Survey (USGS), were analyzed. IDW was tested using five power parameters (p = 1–5), while Kriging applied three semivariogram models (spherical, exponential, and Gaussian). The interpolation accuracy was assessed through Root Mean Square Error (RMSE), Mean Square Error (MSE), and Mean Absolute Percentage Error (MAPE). Results indicated that Kriging with the exponential semivariogram achieved the highest accuracy, with RMSE = 0.0848, MSE = 0.0072, and MAPE = 1.14%, outperforming IDW (RMSE = 0.2288, MSE = 0.0523, MAPE = 1.24%). The Kriging model effectively represented the gradual spatial decay of seismic energy, identifying Aceh Singkil and northern Simeulue as the most earthquake-prone zones, consistent with regional tectonic patterns. These findings confirm that incorporating spatial autocorrelation enhances interpolation accuracy and geophysical interpretation. The study establishes Kriging as a reliable tool for seismic hazard mapping and provides valuable insights for disaster preparedness, infrastructure planning, and future geostatistical applications in earthquake risk assessment.]]>
