In this paper, we first propose a new design strategy of 2D $Z$-complementary array quads (2D-ZCAQs) with feasible array sizes. A 2D-ZCAQ consists of four distinct unimodular arrays satisfying zero 2D auto-correlation sums for non-trivial 2D time-shifts within certain zone. Then, we obtain the upper bounds on the column sequence peak-to-mean envelope power ratio (PMEPR) of the constructed 2D-ZCAQs by using specific auto-correlation properties of some seed sequences. The constructed 2D-ZCAQs with bounded column sequence PMEPR can be used as a potential alternative to 2D Golay complementary array sets for practical applications