In this study, we established a wavelet method, based on Haar wavelets and finite difference scheme for two-dimensional time fractional reaction-subdiffusion equation. First by a finite difference approach, time fractional derivative which is defined in Riemann-Liouville sense is discretized. After time discretization, spatial variables are expanded to truncated Haar wavelet series, by doing so a fully discrete scheme obtained whose solution gives wavelet coefficients in wavelet series. Using these wavelet coefficients approximate solution constructed consecutively. Feasibility and accuracy of the proposed method is shown on three test problems by measuring error in norm. Further performance of the method is compared with other methods available in literature such as meshless-based methods and compact alternating direction implicit methods.