Java, situated in the Pacific Ring of Fire, is one of the most seismically active regions in the world, with frequent earthquakes posing significant risks to its dense population and critical infrastructure. This study aimed to analyze the spatial distribution and intensity patterns of earthquakes in Java from 2022 to 2024 using data from the Meteorology, Climatology, and Geophysics Agency (Badan Meteorologi, Klimatologi, dan Geofisika, BMKG). Spatial interpolation techniques—inverse distance weighted (IDW), nearest neighbor, and Thiessen polygon—were applied to evaluate their effectiveness in mapping earthquake intensity patterns. The dataset included the earthquake magnitude, location, and occurrence time, with performance evaluated using mean absolute percentage error (MAPE) and mean absolute error (MAE). Results showed that the nearest neighbor method achieved the highest accuracy (MAPE of 12.27%, MAE of 0.37), followed by IDW, while the Thiessen polygon method demonstrated limited suitability for continuous seismic phenomena. These findings underscore the importance of selecting appropriate interpolation methods for seismic risk mapping, providing actionable insights for disaster preparedness and urban planning in Java.