The origin of the Steens-Columbia River (SCR) flood basalts, which is presumed to be the onset of Yellowstone volcanism, has remained controversial, with the proposed conceptual models involving either a mantle plume(1-5) or back-arc processes(6-8). Recent tomographic inversions based on the US Array data reveal unprecedented detail of upper-mantle structures of the western USA(9) and tightly constrain geodynamic models simulating Farallon subduction, which has been proposed to influence the Yellowstone volcanism(5,6). Here we show that the best-fitting geodynamic model(10) depicts an episode of slab tearing about 17 million years ago under eastern Oregon, where an associated sub-slab asthenospheric upwelling thermally erodes the Farallon slab, leading to formation of a slab gap at shallow depth. Driven by a gradient of dynamic pressure, the tear ruptured quickly north and south and within about two million years covering a distance of around 900 kilometres along all of eastern Oregon and northern Nevada. This tear would be consistent with the occurrence of major volcanic dikes during the SCR-Northern Nevada Rift flood basalt event both in space and time. The model predicts a petrogenetic sequence for the flood basalt with sources of melt starting from the base of the slab, at first remelting oceanic lithosphere and then evolving upwards, ending with remelting of oceanic crust. Such a progression helps to reconcile the existing controversies on the interpretation of SCR geochemistry and the involvement of the putative Yellowstone plume. Our study suggests a new mechanism for the formation of large igneous provinces.