We study the 3-D P wave velocity (Vp) structure of the crust and upper mantle beneath Greenland and surrounding regions using the latest P wave arrival time data. The Greenland Ice Sheet Monitoring Network (GLISN), initiated in 2009, is an international project for seismic observation in these regions using 34 stations. We use a regional seismic tomography method to simultaneously invert P wave arrival times of local earthquakes and P wave relative traveltime residuals of teleseismic events. These data are extracted from the ISC-EHB catalog; however, for the teleseismic data, we picked new arrival times from seismograms using a cross-correlation method. Our tomographic results clearly reveal the Iceland plume, the Jan Mayen plume, and a newly discovered "Svalbard plume," which merge together in the mantle transition zone. A high-Vp body exists beneath the Greenland Sea, which might act as an obstacle against the rising Svalbard plume. Furthermore, our results reveal a remarkable low-Vp anomaly elongated in the NW-SE direction at depths ≤250 km beneath central Greenland, which is connected with the Iceland and Jan Mayen hotspots. Although previous studies have suggested a similar feature, our result is the first to show the low-Vp zone existing at all depths in the Greenland lithosphere, and its spatial distribution coincides with a high heat flux region. These characteristics indicate that the low-Vp zone reflects residual heat from the Iceland plume when the Greenland lithosphere passed over this plume at~80-20 Ma. Our results also indicate the possible existence of residual heat from the Jan Mayen plume. Plain Language Summary Greenland is a stable land mass that has preserved~4 Gyr of Earth's history. In its vicinity, there are the Mid-Atlantic Ridge, the Iceland and Jan Mayen hotspots, and a geothermal area in Svalbard, which indicate the complexity of this region. We apply seismic tomography to analyze the latest data recorded by a new seismograph network, and we obtain detailed 3-D images of the crust and upper mantle. Our results clearly image the conduit-like Iceland plume, the Jan Mayen plume, and a newly discovered "Svalbard plume," which merge together in the mantle transition zone. These new findings will improve our understanding of the tectonic and thermal evolution of this region. We also find a low seismic velocity zone running in the NW-SE direction beneath central Greenland, which is connected with the Iceland and Jan Mayen hotspots. This zone is located within the Greenlandic plate, and its spatial distribution agrees well with an area of high geothermal heat flux. These results suggest that the low-velocity zone reflects residual heat from the Iceland plume when the Greenlandic plate passed over the plume at~80-20 Ma. Our model further suggests a possible heat track left by the Jan Mayen plume.