Objective: To examine ten years of earthquake data recorded across Indonesia drawing on 5,364 events with magnitudes above M 5.0 between 2016 and 2025. Method: DBSCAN algorithm was run after the optimal neighborhood radius was determined objectively from a k-distance plot. An elbow at about 65 km was identified and the value yielded 16 spatially distinct clusters alongside 460 noise events. K-means algorithm identified four seismic regimes. Results: Of the four regimes, one cluster (Cluster 1) concentrated every major earthquake in the catalog (64 events with M >= 7.0), even though it accounted for fewer than one event in ten. The three remaining clusters captured background seismicity at near-identical mean magnitudes of approximately from 5.33 to 5.35. At the conventional zonal level, Maluku-Sulawesi generated the most events about 40.8% from total events, while Sumatra registered the highest seismic energy output. A Gutenberg-Richter b-value of 0.98 was estimated for the full catalog. Novelty: Introducing earthquake zonation methods based on machine learning for earthquake catalog of Indonesia. These findings support multiple Sustainable Development Goals including the identification of underestimated high-energy rupture corridors informs evidence-based urban risk reduction (SDG 11), strengthens the scientific foundation for earthquake disaster preparedness (SDG 13), introduces an innovative and reproducible machine learning methodology applicable to infrastructure (SDG 9), and contributes a freely transferable workflow that adopt data-driven zonation methods (SDG 17)