Abstract Most earthquakes occur at plate boundaries, but some also strike within stable continental interiors. Although dominant causes of such intraplate earthquakes remain elusive, a prevailing hypothesis attributes intraplate stress and seismicity to variations in lithospheric thickness. Here, we test this hypothesis using the Korean Peninsula as a natural laboratory by constructing three‐dimensional numerical models that extend from the surface to a depth of 650 km, incorporating realistic plate boundary configurations and deep thermal and compositional heterogeneities. Our results show that lithospheric thickness variations alone cannot account for the observed seismicity distribution and stress orientations. Instead, models that include weak subduction interfaces, slabs extending into the transition zone, and mantle buoyancy—while excluding shallow lithospheric density contrasts—explain the observed seismicity and stress distribution. In particular, the subducted Pacific slab in the mantle transition zone acts as a gravitational sinker, enhancing crustal compression along the eastern margin of the peninsula.
Lee et al. (Thu,) studied this question.