A mature mathematical technique called copula joint function is introduced in this paper, which is commonly used in the financial risk analysis to estimate uncertainty. The joint function is generalized to the n-dimensional Frank’s copula. In addition, we adopt two attenuation models proposed by YU and Boore et al, respectively, and construct a two-dimensional copula joint probabilistic function as an example to illustrate the uncertainty treatment at low probability. The results show that copula joint function gives us a better prediction of peak ground motion than that resultant from the simple linear weight technique which is commonly used in the traditional logic-tree treatment of model uncertainties. In light of widespread application in the risk analysis from financial investment to insurance assessment, we believe that the copula-based technique will have a potential application in the seismic hazard analysis. A mature mathematical technique is introduced in this paper, which is commonly used in the financial risk analysis to estimate uncertainty. The joint function is generalized to the n-dimensional Frank’s copula. YU and Boore et al, respectively, and construct a two-dimensional copula joint probabilistic function as an example to illustrate the uncertainty treatment at low probability. The results show that copula joint function gives us a better prediction of peak ground motion than that resultant from the simple linear weight technique which is commonly used in the traditional logic-tree treatment of model uncertainties. In light of widespread application in the risk analysis from financial investment to insurance assessment, we believe that the copula-based technique will have a potential application in the seismic hazard analysis.