A one-dimensional Particle-In-Cell Monte Carlo Collisions (PIC-MCC) method has been used to model the development and propagation of ionization waves in neon and argon positive columns. Low current conditions are considered, i.e. conditions where stepwise ionization or Coulomb collisions are negligible (linear ionization rate). This self-consistent model describes the development of self-excited moving striations, reproduces many of the well-known experimental characteristics (wavelength, spatial resonances, potential drop over one striation, electron "bunching" effect) of the ionization waves called p, r and s waves in the literature and sheds light on their physical properties and on the mechanisms responsible for their existence. These are the first fully kinetic self-consistent simulations over a large range of conditions reproducing the development of p, r and s ionization waves. Although the spatial resonances and the detailed properties of the striations in the non-linear regime are of kinetic nature, the conditions of existence of the instability can be obtained and understood from a linear stability analysis of a three-moment set of quasi-neutral fluid equations where the electron transport coefficients are expressed as a function of electron temperature and are obtained from solutions of a 0D Boltzman equation. An essential aspect of the instability leading to the development of these striations is the non-Mawellian nature of the electron energy distribution function in the uniform electric field prior to the instability onset, resulting in an electron diffusion coefficient in space much larger than the energy diffusion coefficient.