In this paper, the spatial order of accuracy of multiresolution time-domain methods using basis functions from the Daubechies family are studied. It is observed that MRTD methods using scaling functions from the Daubechies N-vanishing-moment orthonormal family have a spatial order of accuracy equal to 2N. For MRTD methods using scaling functions from the CDF (2, N 2) biorthogonal family, the spatial order of accuracy is 2 + N 2. We also find that the order of accuracy is the same for both staggered and nonstaggered spatial grids. The simulation results are presented which confirm analytical predictions. © 2005 Wiley Periodicals, Inc.