A thermodynamic model is developed to explain low temperature metalorganic chemical vapor deposition of GaN. The two stage epitaxial processes taking into account the physisorbed states on (0001) A and B planes are considered. The physisorbed Ga coverage of both GaN planes is presented as the function of the pressures, temperature, and efficiency of ammonia decomposition. The thermodynamic affinities of each stage of the growth runs are described. The formation conditions of the GaN wurtzite and zinc-blende modifications are discussed. The model developed was applied to estimate the thermodynamic characteristics of low temperature metalorganic chemical vapor deposition of GaN. It is quantitatively shown that the considered growth conditions allows formation of the stable and metastable GaN modifications in the same epitaxial run.