The O(1S) (green line) night airglow emission in response to atmospheric gravity wave (AGW) perturbations was simulated with a linear, one-dimensional model. The results were combined with previously modeled O2(b, 0–1) atmospheric band and OH Meinel band emission response (Liu and Swenson, 2003) to derive amplitude and phase relations among multiple airglow layers in response to gravity waves with various intrinsic parameters and damping rates (β). The simulations show that the vertical profile of the standard deviation of the perturbed green line volume emission rate (VER) has a centroid altitude that is 3 km lower and a full-width-half-maximum 2.1 km smaller than the unperturbed VER profile, similar to findings for the OH and O2(b) band layers. Relative phase differences and amplitudes of vertically propagating waves can be deduced from zenith observations of the layers. Airglow weighted responses to waves are related through a cancellation factor (CF) for both layer intensity and temperature. The vertical wavelength can be deduced from relative phase information of three airglow layers separated in altitude. The vertical flux of horizontal momentum associated with gravity waves is deduced from intrinsic wave parameters. Wave damping versus altitude is used to deduce the flux divergence and local accelerations resulting from dissipative waves. The simulations are useful in calculating wave information and wave effects on the atmosphere from multiwavelength, zenith airglow observations.