A method for solving elastodynamic problems in radially inhomogeneous elastic materials with spherical anisotropy is presented, i.e. materials having c ijkl = c ijkl (r) in a spherical coordinate system {r, q, f}. The time-harmonic displacement field u(r, q, f) is expanded in a separation of variables form with dependence on q, f described by vector spherical harmonics with r-dependent amplitudes. It is proved that such separation of variables solution is generally possible only if the spherical anisotropy is restricted to transverse isotropy (TI) with the principal axis in the radial direction, in which case the amplitudes are determined by a first-order ordinary differential system. Restricted forms of the displacement field, such as u(r, q), admit this type of separation of variables solution for certain lower material symmetries. These results extend the Stroh formalism of elastodynamics in rectangular and cylindrical systems to spherical coordinates.