Atmospheric physics is a very complicated natural phenomenon and needs to simplify its basic models for the sea-air oscillator. And it is solved by using the approximate method. The variational iteration method is a simple and valid method. In this paper the coupled system for a sea-air oscillator model of interdecadal climate fluctuations is considered. Firstly, through introducing a set of functions, and computing the variations, the Lagrange multipliers are obtained. And then, the generalized expressions of variational iteration are constructed. Finally, through selecting appropriate initial iteration from the iteration expressions, the approximations of solution for the sea-air oscillator model are solved successively.