We study properties of spherical Bernstein-Bezier splines. Algorithms for practical implementation of the global splines are presented for a homogeneous case as well as a non-homogeneous. Error bounds are derived for the global splines in terms of Sobolev type spherical semi-norms. Multiple star technique is studied for the minimal energy interpolation problem. Numerical summary supporting theoretical considerations is provided.