This paper presents an application of the bivariate Gaussian copulas in describing—by means of a statistical distribution function— station impinging lightning overvoltages, which could be subsequently employed in the statistical and semi-statistical methods of station insulation coordination. The proposed method further employs a state-of-the-art transmission line model for lightning-surge transient analysis, constructed within the EMTP software. It also makes use of the electrogeometric model of lightning attachment to transmission lines in order to estimate the shielding failure and backflashover rates. The proposed method efficiently generates random lightning currents from these non-normal and correlated bivariate statistical probability distributions and employs them in the Monte– Carlo simulation for the purpose of station equipment insulation coordination.