You can also use a cross product to determine the ang rotation for the sphere based on the linear vel at the radius. The problem you will come across is gimble lock which will cause axis to flip when angles exceed 90 degrees. Here you can either do the math fully with quaternions or take the matrix and added euler, conver to quaternions, multiply, convert back to a matrix then extract the up an forward vecs to get the correct rotation.