Multizone soundfield reproduction with various potential applications has recently been drawn more and more attention in acoustic signal processing. However, the literatures related to multizone soundfield reproduction are mainly limited to 2D multizone scenarios. In this paper, based on spherical harmonic analysis, we seek to recreate a 3D multizone soundfield using a single spherical loudspeaker array. The desired local soundfield coefficients are mapped into the equivalent global soundfield coefficients using spherical addition theorem for Bessel function. The mode matching approach is used to derive loudspeaker weights in the reproduction step. The experiments show that the proposed method can reproduce the target 3D spatial multizone soundfield accurately.Index Terms-Acoustic signal processing, 3D spatial multizone soundfield, spherical addition theorem, spherical loudspeaker array