Two airglow models, OHCD-90 and MACD-90, are used to simulate the induced variations of O(1S) greenline, O2(0,1) atmospheric band, and OH(8,3) airglow by the CO2 increase, the F10.7 solar cycle variation and geomagnetic activity for a time period from 1960 to 2015. Our results demonstrate that airglow intensity and the peak Volume Emission Rate (VER) variations of the three airglow emissions are strongly correlated, and in phase, with the F10.7 solar cycle variation. In addition, there is a linear trend, be it increasing or decreasing, existing in the airglow intensities and VERs due to the CO2 increase. On other hand, airglow VER peak heights are strongly correlated, and out of phase, with the Ap index variation of geomagnetic activity. The CO2 increase acts to lower the VER peak heights of OH(8,3) airglow and O(1S) greenline by 0.2 km in 55 years and it has no effect on the VER peak height of O2(0,1) atmospheric band.