Seismic interferometry is a process of generating new seismic data from existing wavefields. This enables us to expand the degree of freedom of source-receiver configuration. Seismic interferometry by multidimensional deconvolution (MDD) is proposed as an alternative to the conventional crosscorrelation method (Wapenaar et al., 2008a; Wapenaar et al., 2008b). We apply MDD to crosswell geometry in order to retrieve crosswell impulse responses from surface sources using numerical modeling and field data. We adopted singular-value decomposition (SVD) for obtaining the pseudoinverse solution to achieve MDD. Since the SVD pseudoinverse is highly dependent on the rank of the MDD matrix, Akaike's information criterion (AIC) is adopted in order to determine the rank of the MDD matrix. We see that the MDD produces higherresolution data compared with the conventional crosscorrelation method. Furthermore, amplitudes of downgoing reflection events are improved in MDD while downgoing reflection events are not recognizable in the crosscorrelation method.