Proton emission from a laser ion source Rev. Sci. Instrum. 83, 02B310 (2012) Optical diagnosis and theoretical simulation of laser induced lead plasma spectrum Phys. Plasmas 19, 013302 (2012) Saturation effects in femtosecond laser ablation of silicon-on-insulator Appl. Phys. Lett. 99, 231108 (2011) Effect of air breakdown with a focusing lens on ultrashort laser ablation Appl. Phys. Lett. 99, 234104 (2011) Additional information on J. Appl. Phys. An analytic kinetic model is elaborated to determine the spatial evolution of the concentration of ground and excited Li atoms ͑up to the 3 2 D level͒ present in the plasma generated by pulsed laser ablation of a LiNbO 3 crystal in gas atmospheres ͑at 1 Torr͒ of different nature, both reactive (O 2 ) and inert ͑Ar͒. Whereas the Li atom population densities in the plasma produced in gas are, in general, higher than those found in vacuum, the concentration of ground Li atoms close to the substrate is similar in vacuum and Ar, though lower than in O 2 . It is suggested that the reabsorption of radiation by the plasma is more important within O 2 than in Ar since the absorber ͑ground Li͒ concentration is higher in O 2 environments. A very reasonable qualitative agreement is found when comparing the Li I-670.8 nm spectral line emission intensity predicted by the model with that obtained experimentally in Ar and O 2 atmospheres.