Abstract: Current models of the spectrum of X-rays reflected from accretion discs around black holes and other compact objects are commonly calculated assuming that the density of the disc atmosphere is constant within several Thomson depths from the irradiated surface. An important simplifying assumption of these models is that the ionization structure of the gas is completely specified by a single, fixed value of the ionization parameter (xi), which is the ratio of the incident flux to the gas densit…