In this paper we consider the evolution equation ∂tu = ∆µu + f and the corresponding Cauchy problem, where ∆µ represents the Bessel operator ∂ 2x + ( 1 4 − µ 2 )x −2 , for every µ > −1. We establish weighted and mixed weighted Sobolev type inequalities for solutions of Bessel parabolic equations. We use singular integrals techniques in a parabolic setting. Recently, Jones' results have been extended to weighted and mixed weighted L p spaces by Ping, Stinga and Torrea ([42, Theorem 2.4]). In [42], they also considered the equation (1.1) on the whole space R n+1 . In [42, Theorem 2.3], it was proved that if f is a C 2 -function defined on R n+1 with compact support, when n ≥ 2, the function u given byis a solution of (1.1) on R n+1 , where ∂ 2 xi,xj u(t, x) = lim ǫ→0 + Ωǫ