Abstract—We investigate the contribution of gamma radiation of natural radionuclides constituting the Earth crust, radioactive emanations, and their decay product in the ground to the intensity of production of ion pairs in the atmosphere against the background of ionization of the atmosphere by radioactive gases flowing to the atmosphere from the ground and propagating together with their short-lived daughter products. The density of the radon flux to the atmosphere is estimated by three methods: the reservoir method, the integration of altitude profiles of volume activity of radon, based on gamma spectroscopic observation and diffusion model. The distribution of the gamma dose rate dose from the earth radionuclides in the soil and in the atmosphere is calculated using Gleant4 software. The propagation of the radon isotopes and their decay products in the atmosphere is calculated using large eddy simulations supplemented with kinematic simulations of the subgrid transport of a passive scalar. It is shown that depending on the specific activity of nuclides in the ground and the turbulent regime of the atmosphere, the total contribution of the γ-radiation to the ion pair production rate in the atmospheric boundary layer is approximately from 1% to 20% and increases upon a decrease in the penetrability of the upper ground layer for radioactive emanations.