Multizone　soundfield　reproduction　in　the　reverberant　room　is　a　complex　problem　in　acoustic　signal　processing.　In　this　paper,　we　introduce　a　method　of　recreating　a　3D　multizone　soundfield　for　reverberant　environments.　Based　on　spherical　addition　theorem　for　Bessel　function,　we　derive　the　desired　global　soundfield　by　translating　local　desired　soundfields　to　a　single　global　coordinate　system.　We　cast　the　global　soundfield　reproduction　into　a　least　squares　framework　and　introduce　the　modal　space　approach　to　derive　the　loudspeaker　weights.　A　practical　method　is　used　to　determining　the　Acoustic　Transfer　Function　(ATF)　between　each　loudspeaker　and　every　point　of　reproduction　region.　The　simulation　results　indicate　that　the　proposed　method　can　reproduce　the　target　3D　multizone　soundfield　accurately.