In this study, the mechanical properties of the composite plate were considered Gaussian random fields and their effects on the buckling load and corresponding mode shapes were studied by developing a semi-analytical non-intrusive approach. The random fields were decomposed by the Karhunen−Loève method. The strains were defined based on the assumptions of the first-order and higher-order shear-deformation theories. Stochastic equations of motion were extracted using Euler–Lagrange equations. The probabilistic response space was obtained by employing the non-intrusive polynomial chaos method. Finally, the effect of spatially varying stochastic properties on the critical load of the plate and the irregularity of buckling mode shapes and their sequences were studied for the first time. Our findings showed that different shear deformation plate theories could significantly influence the reliability of thicker plates under compressive loading. It is suggested that a linear relationship exists between the mechanical properties’ variation coefficient and critical loads’ variation coefficient. Also, in modeling the plate properties as random fields, a significant stochastic irregularity is obtained in buckling mode shapes, which is crucial in practical applications.