Unavoidable geometric imperfections complicate the load–displacement response of thin-walled cold-formed steel (CFS) structural members. While the inclusion of geometric imperfections in shell finite element (FE) collapse simulations has long been the subject of research, no consensus exits on modeling their full three-dimensional spatial pattern. This paper is the first part in a two-part series and provides a review of the basics and a summary of some efforts geared towards modeling geometric imperfections. We discuss the necessary ingredients for systematic characterization of elastic buckling mode shapes, due to their frequent use in approximating geometric imperfections, and make a conscious effort to express some of the existing geometric imperfection models through three components: imperfection shape, magnitude, and combination coefficient. Different choices and selections of these components, pertaining to either the inherent buckling behavior of the members and/or physical measurements are critically examined. The details required to develop high fidelity nonlinear FE models of CFS members seeded with geometric imperfections are presented. We then perform a comparative study of collapse analysis results obtained from adopting different imperfection modeling strategies. The noticeable variability in the load carrying capacity and the stability behavior highlights an urgent need for reliable probabilistic frameworks for this class of problems. A stochastic framework to generate samples of imperfect shell FE models from limited available data and validate the prescribed or practice-oriented geometric imperfection models is presented in a companion paper (see Part II Farzanian et al., (2023)).