The satellite and near-ground study demonstrated IR temperature anomalies associated with fault system of the crust in the seismic region of the Earth. The amplitude of anomalies is about ΔT = 3 K. Such change of temperature needs about 100 % increase of CO₂ concentration. Internal gravity waves (IGW) excited by greenhouse effect gases could cause some ionospheric disturbances including density irregularities. The accurate numerical model of 2D «lithospheric greenhouse effect gas antenna* with heat, mass and concentration sources of IGW in the atmosphere is built in the present work. The system of hydrodynamics equations for IGW excitation in the atmosphere by the source in the form of near-ground layer of lithospheric greenhouse effect gases is reduced to the system of two equations for pressure and vertical velocity component. Corresponding effective boundary conditions are obtained by means of limit pass to the case of very thin layer and absolutely rigid lithosphere. Periodical boundary conditions in the horizontal directions are used to avoid computations with continuous spectrum and numerical convergence of the model is checked carefully. It is shown that at the altitude Z = 200 km the value of vertical velocity of IGW with period 1 hour could reach the value of order 4 m/s what is enough for plasma bubbles formation in accordance with previously published data. In distinction to the qualitative model of Gohberg et al., reactive modes in the IGW spectrum are taken into account. It is shown that, although these modes are non-propagating, they influence significantly the IGW excitation and, as a result, the characteristics of propagating modes. In particular, non-taking into account reactive modes could cause under definite conditions overestimating of energy flow by at least two orders of value. This result is analogous to the well known in microwave technique effect of the influence of reactive modes on scattering of electromagnetic waves on inhomogeneities in microwave waveguides.