In this article, the notion of orthogonal nonseparable four-dimensional wavelet packs, which is the generalization of orthogonal univariate wavelet packs, is introduced. A new approach for constructing them is presented by iteration method and wavelets as well wavelet frames. The biorthogonality properties of four-dimensi- -onal wavelet packets are discussed. Three biorthogonality formulas concerning these wavelet packs are estabished. A necessary and sufficient condition for the existence of the pyramid decomposition scheme of space is presented.