An inverse source problem (ISP) of recovering space‐dependent source term along with diffusion concentration for a two‐dimensional diffusion equation involving integral convolution of arbitrary memory kernel in time is considered. The unique existence of the solution is proved using a bi‐orthogonal system of functions obtained from the associated non‐self‐adjoint spectral problem and its adjoint problem which form Riesz basis in
L2false[false(0,1false)×false(0,1false)false]$$ {L}^2\left[\left(0,1\right)\times \left(0,1\right)\right] $$. In addition, some particular cases of ISP are described as special cases of our results.