The main result in this paper is an error estimate for interpolation biharmonic polysplines in an annulus A (r1, rN ), with respect to a partition by concentric annular domains A (r1, r2) , ...., A (rN−1, rN ) , for radii 0 < r1 < .... < rN . The biharmonic polysplines interpolate a smooth function on the spheres |x| = rj for j = 1, ..., N and satisfy natural boundary conditions for |x| = r1 and |x| = rN . By analogy with a technique in onedimensional spline theory established by C. de Boor, we base our proof on error estimates for harmonic interpolation splines with respect to the partition by the annuli A (rj−1, rj). For these estimates it is important to determine the smallest constant c d (Ω) , where Ω = A (rj−1, rj) , among all constants c satisfying supvanishing on the boundary of the bounded domain Ω . In this paper we describe c d (Ω) for an annulus Ω = A (r, R) and we will give the estimate min{ 1 2d , 1 8 } (R − r) 2 ≤ c d (A (r, R)) ≤ max{ 1 2d , 1 8 } (R − r) 2where d is the dimension of the underlying space.