We describe the first steps in the extension of the Symanzik O(a) improveme
nt program for Wilson-type quark actions to anisotropic lattices, with a te
mporal lattice spacing smaller than the spatial one. This provides a fully
relativistic and computationally efficient framework for the study of heavy
quarks. We illustrate our method with accurate results for the quenched ch
armonium spectrum.