The objective of this paper is to investigate the effect of copulas for constructing the bivariate distribution of shear strength parameters on system reliability of geotechnical structures. First, the bivariate distribution of shear strength parameters is constructed using copulas. Second, the implementation procedure of system reliability analysis using direct Monte Carlo simulation (MCS) is developed. Finally, the system reliability of a retaining wall and a rock wedge slope is presented to explore the effect of copula selection on geotechnical system reliability. The results indicate that the system reliability of geotechnical structures under incomplete probability information could not be determined uniquely because the bivariate distribution of cohesion and friction angle with given marginal distributions and correlation coefficient could not be determined uniquely. The copulas for modeling dependence structure between cohesion and friction angle have a significant influence on the system reliability of geotechnical structures. Such an influence includes two phases separately. The first phase is that the dependence structure between shear strength parameters characterized by copulas affects the reliability of single failure mode, depending on the marginal distributions, dependence structure between shear strength parameters, and reliability level of each failure mode. The second phase is that the reliability of each failure mode influences on system reliability, only depending on reliability level of each failure mode and correlations among various failure modes. It is important to distinguish between the effect of copula selection on reliability of each failure mode and that on geotechnical system reliability.