Canonical representations of Sp(1,n) associated with finite dimensional irreducible representations of Sp(1) are defined using vector-valued Berezin kernels. Their decomposition into irreducible representations is determined by decomposing the corresponding reproducing distributions in terms of positive definite trace spherical functions on Sp(1,n). The canonical representatons are also identified with the restriction to Sp(1,n) of certain maximal degenerate representations of SL(n+1,H).