In this paper, seismic analyses of single-layer dome structures are carried out, with a particular emphasis on the incorporation of randomness in both geometrical imperfections and ground motions. First, the paper introduces models for random geometrical imperfections and stochastic ground motions, employing the Stochastic Imperfection Mode Superposition Method (SIMSM) for imperfections and utilizing the Spectral Representation Method (SRM) for ground motions. These models are designed to comprehensively account for the inherent uncertainties in real-world structural behavior. Subsequently, the Probability Density Evolution Method (PDEM) is applied to investigate the probabilistic response and seismic reliability of dome structures, which allows for a detailed examination of key performance indicators, such as displacement, plastic ratio, and strength damage, under the influence of random factors. Both deterministic and stochastic analyses are conducted to gain insights into the impact of randomness in geometrical imperfections and ground motions. The computational results unveil a crucial finding: an unfavorable distribution of imperfections can significantly compromise the dynamic performance and safety of single-layer dome structures. Furthermore, the paper delves into parametric analyses, with a specific focus on imperfection amplitude, the rise-span ratio and roof load as pivotal parameters. These analyses provide valuable insights into the sensitivity of the structural response to variations in these critical design factors.