Lattice Boltzmann models for diffusion equation are generally in Cartesian coordinate system. Very few researchers have attempted to solve diffusion equation in spherical coordinate system. In the lattice Boltzmann based diffusion model in spherical coordinate system extra term, which is due to variation of surface area along radial direction, is modeled as source term. In this study diffusion equation in spherical coordinate system is first converted to diffusion equation which is similar to that in Cartesian coordinate system by using proper variable. The diffusion equation is then solved using standard lattice Boltzmann method. The results obtained for the new variable are again converted to the actual variable. The numerical scheme is verified by comparing the results of the simulation study with analytical solution. A good agreement between the two results is established.