The present paper is concerned with non-isothermal spherical Couette flow of Oldroyd-B fluid in the annular region between two concentric spheres. The inner sphere rotates with a uniform angular velocity while the outer sphere is kept at rest. Moreover, the two spherical boundaries are maintained at fixed temperature values. Hence, the fluid is effect by two heat sources; namely, the viscous heating and the temperature gradient between the two spheres. The viscoelasticity of the fluid is assumed to dominate the inertia such that the latter can be neglected. An approximate analytical solution of the energy and momentum equations isobtained through the expansion of the dynamical fields in power series of Nahme number. The analysis show that, the temperature variation due to the external source appears in the zero order solution and its effect extends to the fluid velocity distribution up to present second order. Viscous heating contributesin the first and second order solutions. In contrast to isothermal case, a first order axial velocity and a second order stream function fields has been appeared. Moreover, at higher orders the temperature distribution depends on the gap width between the two spheres. Finally, there exist a thermal distribution of positive and negative values depend on their positions in the domain region between the two spheres.