A simple transformation, for a point on the surface of a spherical model and for the corresponding local velocity, has been used to study the propagation of a disturbance along the surface. Using a Gaussian beam computation scheme, rays are traced and seismograms obtained for a few spherical models with lateral heterogeneities. In the results, deviations from the first order perturbation theory for normal modes are displayed.