For the large scale reinforced concrete cooling tower shells, the shape imperfection can be introduced due not only to mistakes in the process of construction but also to the long term behavior of concrete. The shape imperfection evokes the additional responses such as displacements and stresses in addition to the design values. In this study, the statistical behavior of the RC cooling tower shell due to the shape imperfection is investigated using the Monte Carlo simulation. The radius of cooling tower and the shell thickness are adopted as the parameters which cause the shape imperfection. The shape imperfection is modeled as a stochastic field rather than the local one of axisymmetric or bulge type of imperfection. The randomness in the radius is shown to be more affecting the structural responses than the randomness in the shell thickness. In addition to the geometrical randomness, the effect of randomness in the modulus of elasticity on the structural response is also investigated and compared with that of the geometrical ones.