Abstract A two-dimensional, anelastic, cloud-resolving numerical model was used to simulate squall systems. Large domain and fine grid resolutions were utilized so that both the convective and mesoscale components of squall lines could be handled adequately. Detailed cloud microphysics including the ice phase and the Coriolis force have been included in the basic model. Both the life cycle and storm structure of observed squall systems have been simulated successfully. Some details in the observed precipitation and kinematic characteristics of squall lines, such as the locations of front-to-rear jet core, the base of the stratiform cloud, the formation of a transition zone, and the organized mesoscale updraft, have been simulated by the model. The storm-generated meso-γ-scale low pressure center located behind the convective updraft has been shown to be instrumental in the initiation and maintenance of the mesoscale circulation and the associated trailing stratiform region. Diagnostically, the horizontal ...