Abstract. The modeling of very low-frequency (VLF) electromagnetic (EM) beam propagation in the Earth–ionosphere waveguide (WGEI) is considered. A new tensor impedance method for modeling the propagation of electromagnetic beams in a multi-layered and inhomogeneous waveguide is presented. The waveguide is assumed to possess the gyrotropy and inhomogeneity with a thick cover layer placed above the waveguide. The influence of geomagnetic field inclination and carrier beam frequency on the characteristics of the polarization transformation in the Earth–ionosphere waveguide is determined. The new method for modeling the propagation of electromagnetic beams allows us to study the (i) propagation of the very low-frequency modes in the Earth–ionosphere waveguide and, in perspective, their excitation by the typical Earth–ionosphere waveguide sources, such as radio wave transmitters and lightning discharges, and (ii) leakage of Earth–ionosphere waveguide waves into the upper ionosphere and magnetosphere. The proposed approach can be applied to the variety of problems related to the analysis of the propagation of electromagnetic waves in layered gyrotropic and anisotropic active media in a wide frequency range, e.g., from the Earth–ionosphere waveguide to the optical waveband, for artificial signal propagation such as metamaterial microwave or optical waveguides.