It is shown that Lie group analysis of differential equations provides the exact solutions of two-dimensional stratified rotating Boussinesq equations which are a basic model in geophysical fluid dynamics. The exact solutions are obtained as group invariant solutions corresponding to the translation and dilation generators of the group of transformations admitted by the equations. The comparison with the previous analytic studies and experimental observations confirms that the anisotropic nature of the wave motion allows to associate these invariant solutions with uni-directional internal wave beams propagating through the medium. It is also shown that the direction of internal wave beam propagation is in the transverse direction to one of the invariants which corresponds to a linear combination of the translation symmetries. Furthermore, the amplitudes of a linear superposition of wave-like invariant solutions forming the internal gravity wave beams are arbitrary functions of that invariant. Analytic examples of the latitude-dependent invariant solutions associated with internal gravity wave beams that have different general profiles along the obtained invariant and propagating in the transverse direction are considered. The behavior of the invariant solutions near the critical latitude is illustrated. © 2009 Elsevier B.V. All rights reserved.