Passive seismic methods using earthquakes can be applied for extracting body waves and obtaining information of subsurface structure. In this study, we retrieve direct and reflected plane waves by applying seismic interferometry to the recorded ground motion from a cluster of regional earthquakes. We apply upgoing/downgoing <it>P/S</it> wavefield decomposition, time windowing, and multidimensional deconvolution to improve the quality of the extraction of reflected waves with seismic interferometry. The wavefield separation and seismic interferometry based on multidimensional deconvolution allow us to reconstruct <it>PP, PS, SP</it> and <it>SS</it> reflected waves without unwanted crosstalk between <it>P</it> and <it>S</it> waves. From earthquake data, we obtain <it>PP, PS</it> and <it>SS</it> reflected plane waves that reflect off the same reflector, and estimate <it>P</it>- and <it>S</it>-wave velocities.