Each seismic body wave, including quasi compressional, shear, and converted wave modes, carries useful subsurface information. For processing, imaging, amplitude analysis, and forward modeling of each wave mode, we need approximate equations of traveltime, slope (ray-parameter), and curvature as a function of offset. Considering the large offset coverage of modern seismic acquisitions, we propose new approximations designed to be accurate at zero and infinitely large offsets over layered transversely isotropic media with vertical symmetry axis (VTI). The proposed approximation for traveltime is a modified version of the extended generalized moveout approximation that comprises six parameters. The proposed direct approximations for ray-parameter and curvature use new, algebraically simple, equations with three parameters. We define these parameters for each wave mode without ray tracing so that we have similar approximate equations for all wave modes that only change based on the parameter definitions. However, our approximations are unable to reproduce S-wave triplications that may occur in some strongly anisotropic models. Using our direct approximation of traveltime derivatives, we also obtain a new expression for the relative geometrical spreading. We demonstrate the high accuracy of our approximations using numerical tests on a set of randomly generated multilayer models. Using synthetic data, we present simple applications of our approximations for normal moveout correction and relative geometrical spreading compensation of different wave modes.