!!!!!!! DRAFT !!!!!!!
Hydroxyl rotational temperature record from the Auroral Station in Adventdalen, Svalbard (78°N, 15°E): Preliminary results.
Sigernes 1, Kristian Nielsen 2, Charles Deehr 3, Trond Svenøe 1,
Nikita Shumilov 4 and Ove Havnes 4
The University Courses on Svalbard, N-9171 Longyearbyen, Norway.
2 The University Courses on Svalbard and the University of Bergen, Norway.
3 Geophysical Institute, University of Alaska Fairbanks, USA.
4 The Auroral Observatory, University of Tromsø, Norway.
This paper describes the compiled mesospheric winter temperature series derived through over 20 years of ground - based spectral measurements of the hydroxyl airglow layer from the Auroral Station in Adventdalen near Longyearbyen, Svalbard (78°N, 15°E). Diurnal and semi diurnal tides are present. Amplitudes of 20 °K with extremes up to 70 °K within a few hours are usual. The average daily winter temperature for the period was found to be 206 °K with a standard deviation of 16 °K. In general the polar mesosphere has a winter warming, but no clear shapes or symmetry around solstice are observed from year to year. The temperature series is interspersed with periods of 2 - 5 days when the temperatures may drops as much as down to the summer minimum. These minima are probably associated with stratospheric warmings that occur with varying strength all winter. Some years are more quite than others with respect to variation in temperature. The annual mean winter temperatures show nearly zero temperature trend, indicating no observational change in the winter mesospheric temperatures over Svalbard during the last two decades.
The natural thermal structure is believed to be primarily controlled by heating from absorption of short wave solar radiation and as a response balanced by cooling related to emissions in the infrared part of the spectrum. The absorption of UV by ozone (O3) constitutes the principal radiative source of heat in the stratosphere and mesosphere. In addition, molecular oxygen (O2) plays an important role. The absorption of UV by molecular oxygen contribute, especially in the upper mesosphere and lower thermosphere. The absorption processes occur as dissociation, forming atomic oxygen (O). At high altitudes (>80km), the lifetime of atomic oxygen exceeds a day, and the energy may be stored as chemical energy. This energy is released as thermal energy when the atom recombines. Due to both horizontal and vertical transport, much of the stored chemical energy is released in the high latitude winter. This process, along with adiabatic heating caused by vertical velocities, is believed to be the main explanation the warm mesopause temperatures observed in winter [cf. Brasseur and Solomon, 1986]. Radiative cooling of the stratosphere and mesosphere is mainly due to vibrational relaxation in the infrared 15 mm band of carbon dioxide (CO2). In addition, the 9.6 mm band of O3 contribute to the cooling, especially near the stratospause. Water vapor (H2O) also contribute, but to a minor extent compared to CO2 and O3 [cf. Kuhn and London , 1969; London, 1980; Fels and Schwarzkopf, 1981].
The radiative equilibrium described above is not enough in order to describe the complete thermal budget of the mesosphere. Measurements of atmospheric atomic oxygen number densities, crucial to the understanding of the thermal budget, have to be intensified in order to get the complete picture [Rees, 1989]. Other natural minor atmospheric constituents may contribute as well to the coupling between chemistry and radiative transfer. In addition, the general circulation with oscillations due to tidal and planetary waves have to be included. Temporal variations due gravity waves propagating through the mesosphere is also of importance [Viereck and Deehr, 1989; Viereck, 1991; Hamilton, 1996]. External perturbations like variation in the solar flux, energetic particle precipitation, volcanic and anthropogenic emissions are additional sources of concern to the thermal structure of the mesosphere. The increase in greenhouse gases like carbon dioxide and methane in the lower atmosphere are expected to increase the water content and decrease the temperature in the upper mesosphere [Thomas, 1996]. Recent model calculations predict that a doubling of the CO2 content is indeed necessary in order to explain the long term negative temperature trend of the mesosphere [Akmaev and Fomichev, 1998].
The first indirect indications on a negative temperature trend in the mesosphere was reported by Gadsden . From the increasing number of observed noctilucent clouds (NLC), the average temperature of the summer mesosphere was estimated to decrease at a rate of 2.5 °K per decade. Since ice crystals are believed to be the major part of the NLC particles, a decrease in temperature may be the reason why they are more frequently formed and sighted at lower latitudes than earlier. The fact that the F2 peak in electron density in the ionosphere has decreased in altitude, as shown by 39 years of ionosonde observations [Ulich and Turunen, 1997], is consistent with a cooling of the lower parts of the atmosphere. Most of the middle atmosphere is evidently cooling, with values which varies from a few °K at 30-40 km, 10 °K at 50 km and 20 °K at 60-70 km since the mid-60s [Golitsyn et al., 1996]. These results where based on direct in situ rocket measurements. Golitsyn et al.  also reported a record high mesospheric cooling of 30 °K since 1957 using temperature estimates from measurements of the hydroxyl (OH) airglow layer (centered at about 87 km) at Zvenigorod (55.7°N) and Abastumani (41.8°N). These last measurements have not been confirmed by results from other stations. Direct observations by the falling sphere technique [Schmidlin, 1991; Lübken, 1999], which give temperatures up to approximately 95 km, does not give any indications of such a large decrease of the mesopause temperature [Lübken, 2000]. Clearly, there is need for more long term temperature measurements from different locations and by different methods.
In addition to the possibility that the above trends signal climatic changes, measurements of parameters such as the temperature and its long and short term variations are important tools in our effort to understand the processes in the mesosphere and thermosphere A change in the temperature will affect the propagation and dissipation of waves, which have a profound influence on these parts of the atmosphere [Bittner et al., 2000]. Long term observation series such as those reported by Golitsyn et al. , Bittner et al.  and Duck et al.  are therefore important. As a contribution to our understanding on the temperature in the mesosphere and its variations, we will in this paper report mesospheric winter temperatures derived through spectral measurements of the OH airglow since 1980 from the Auroral Station in Adventdalen (78°N, 15°E).
The traditional way and most reliable form of monitoring long term temperatures of the mesosphere is through ground based spectroscopic measurements of the OH airglow layer [Sivjee, 1992; Greet et al., 1998]. The molecular bands of OH was first identified by Meinel  and dominate the airglow spectrum from approximately 5200 Å to 4 mm. The OH - airglow layer is, according to in-situ measurements made with rockets, situated around 87km with a thickness of 8 km regardless of season or latitude [Baker and Stair, 1988]. The principal source of the vibration excited OH radicals, the ozone-hydrogen mechanism, was suggested independently by Bates and Nicolet  and Herzberg 
The exothermicity is sufficient to excite the OH molecules up to a vibrational level of n' = 9. The evidence in favor of the ozone-hydrogen mechanism was obtained in the laboratory by McKinley et al. . Note that both H and O3 are minor constituents in the mesosphere and lower thermosphere. Their production and loss needs to be understood in order to appreciate the importance of the OH emissions [cf. Rees, 1989]. The rotational temperature is inferred from the relative photon emission rates of various lines in a band by constructing synthetic spectra that best match the observations. The obtained rotational temperature is interpreted as kinetic under the assumption that the excited OH molecules are in thermal equilibrium with the atmosphere.
The instrument used in this study, the Ebert-Fastie spectrometer, was constructed by W. G. Fastie at John Hopkins University, Maryland at the beginning of the 70's. Fastie [1952a,b] improved the original design of the monochromator made by Herman Ebert in 1889. In all five of these instruments have been produced. In 1978 two were transferred together with a 1/2 m Ebert-Fastie spectrometer to the Auroral Station in Adventdalen, Svalbard, from the Geophysical Institute, University of Alaska. Since then the photon counting- and computer electronics have been continuously upgraded in order to enhance both the control and performance of the instruments [Sigernes, 1996]. Data from these instruments are widely published and recognized. Throughout the 80's, a series of articles have been published on OH temperature measurements from the station [Myrabø at al., 1983; Myrabø, 1984, 1986; Myrabø and Harang, 1988; Sivjee and Hamwey, 1987]. The work was mainly on tidal variations, the effect of stratospheric warmings and the seasonal variation. In the late 80's, focus turned onto the effect of gravity waves interacting with the OH layer [Myrabø et al., 1987; Viereck, 1988; Viereck and Deehr, 1989].
Fig. 1 shows the experimental setup of the 1/2 m focal length Ebert-Fastie spectrometer used for OH measurement in the early 80's. The instrument was later replaced permanently in 1983 by one of the 1 m focal length Ebert-Fastie's in order to gain higher throughput and spectral resolution. The principal components of the instruments are one large 1m focal length spherical mirror, one plane reflective diffraction grating and a pair of curved slits. The recorded radiance from the sky, is limited by the etendue; the product of the area of the entrance slit and the field of view. Because of the low intensity of the source, the etendue is made as large as possible. The image of the entrance slit is reflected by one part of the spherical mirror onto the grating. The second part of the mirror focuses the diffracted light from the grating onto the exit slit. When the grating turns, the image of the entrance slit is scanned across the exit slit. A collector lens transfers the output of the exit slit to the front of a photo multiplier tube. The tube is mounted in a Peltier cooled housing. Signals from the tube are amplified and discriminated before sent to the computers counting card. The instrument is setup to scan through the first order wavelength range 8312 - 8745Å, which corresponds to the OH(6-2) band of the airglow. The band pass is 4.8Å. A red cut-off filter is used in front of the entrance slit to prevent overlapping spectral orders. The spectrometer is fixed in geographical zenith with a field of view of approximately 5°.
In order to derive the rotational temperatures of the measured OH(6-2) emissions, it is necessary to produce the corresponding synthetic spectra. The basis for calculating these spectra is given by Herzberg  with term values from Krassovsky  and Einstein coefficients from Turnbull and Lowe . The temperatures can then be derived since the upper populated energy states of the OH(6-2) band are distributed according to the Stefan-Boltzmann distribution. Recent work by Conner  explains in detail how the parameters necessary for generating the synthetic spectra as a function of the rotational temperature can be calculated. The temperature is found by choosing the optimal fit between the measured and the synthetic spectrum through iteration, to minimize the least square error. A computer program [Sigernes and Nielsen, 2000] has been created to implement the above procedure. The program is capable of analyzing large sets of data and is together with samples of raw data freely available to the scientific community.
Fig. 2 shows a typical result of a synthetic fit to a measured spectrum from 17. November 1999. In this case, a temperature of 204°K with a fit value of 0.95 was found. The fit value is defined as
where l is the least mean square difference between the measured and the synthetic spectrum. This is the same fit function that was used by Myrabø et al. , Viereck  and Viereck and Deehr . A Fit=0.8 represents a relative uncertainty of approximately ± 2°K. In addition, the covariance between synthetic and measured spectrum has been introduced to select time periods where the measurements are not affected by moon light, clouds or auroral activity. The main problem is scattered moonlight and clouds. That the moon is above the horizon is not a direct problem for measuring the OH(6-2) emissions, as these are measured in the zenith and the moon never rises very far above the horizon at Svalbard. However when the sky is cloudy the moonlight is scattered into the spectrometer making the background intensity rise in both strength and variability, which makes the fit of the synthetic spectra poor. A covariance level of 0.9 turned out to work successfully.
Results and discussion
The OH daily mean rotational winter temperature from 1980 to 2001 is plotted in figure 3. The data set consists of 624 days and it is not uniformly sampled in time. Each daily mean is calculated as the average temperature of hourly deduced temperatures. A minimum of 3 hrs are used to obtained a reliable daily mean. The selected number of hourly obtained temperatures per day depends mainly on sky conditions as mentioned above. Furthermore, it also depends on whether the instrument was operative or used in another configuration / setup. These are the main factors for the large data caps seen in the fig. 3. The instrument was not in operation in 81/82 and during the winter season 98/99 the instrument was setup to study proton aurora instead of airglow. In addition, data from 91/92, 92/93 and 96/97 have not yet been retrieved due to problems with storage media and outsourced tape stations. These temperatures will in the near future be restored by using data from other instruments at the station. During the seasons of 80/81 to 82/83 the temperatures were obtained using the OH(8-3) band. Thus, on average, the 80/81 and 82/83 temperatures have been lowered by 5-6 °K to make them compatible with the OH(6-2) data from the 83/84 to the 89/90 seasons. Furthermore, the temperatures up to the 89/90 season are calculated using the transition probabilities from Mies . These temperatures are approximately 5 °K to low compared to using the transition probabilities of Turnbull and Lowe . Both these two effects have been accounted for to make the whole series compatible.
The average daily winter temperature for 1980 to 2001 is 208 °K. The maximum temperature is 257 °K in January 14 1983. Minimum temperatures down to 150 °K are seen in January of 1983 and 1991. The standard deviation is close to 16°K , indicating that the mesospheric temperature variations over Svalbard in winter are extremely high. The day to day temperatures show large peak like variations lasting for a few days to a week or longer. Note that large short time-scale variations in the temperature may occur within each daily average. Myrabø at al.  detected large amplitude variations of both the OH intensity and the temperature, and were the first to suggest that these variations are actual variations rather than noise. Amplitudes of 20 °K in the internal gravity wave period range (minutes to hours) with extremes up to 70 °K within a few hours are usual.
On a larger time scale, the daily mean temperatures show a wave-like pattern with deviations from the mean, each lasting for about 3-10 days. Myrabø et al.  analyzed the seasonal variation from 80/81, 82/83, 83/84 and 84/85 and reported a semidiurnal tide with an amplitude of 3 °K in the December 1982 data. By early January 1983 the tide increased in amplitude to 10 °K. The increase was correlated with a stratospheric warming [Labitzke, 1981; Labitzke et al., 1985]. A cooling of the mesosphere around New Years with a total amplitude of 20 - 30 °K tends to occur each year. In addition, Walterscheid et al.  reported a semidiurnal tide of 13 °K on day 359 in 1984. The tide was called pseudo - tide, believed to be generated by gravity-wave-mean-flow interactions modulated by the solar-driven tide. Nielsen et al. [2001a] reported a cooling of 30 °K in December 1987, which again was associated with a stratospheric warming. The strength of the diurnal and semi diurnal tides were 30 and 15 °K, respectively.
In general the mesosphere has a winter warming, peaking near solstice. It is interspersed with a period of 2 - 5 days when the temperature drops as much as down to the summer minimum. These minimums are probably associated with stratospheric warmings that occur with varying strength all winter. The 80/81 to 84/85 seasons have a relative low temperatures in late December followed by a warm mesospause in early January. There are exceptions, the 85/86, 89/90 and 99/00 seasons show the opposite tendency. There are also more quite years with smoother seasonal variation. A clear example is the 99/00 season. There is in other words no distinct identifiable shape of the seasonal temperature variation from year to year. Comparison with other OH-measurements to similar measurements performed at mid and low latitude stations show that both short period and seasonal temperature variations are consistently larger in the OH airglow over Svalbard. Similar measurements from Antarctica show smoother temperature variations compared to Svalbard. Both the daily and weekly variations of the temperatures exhibit variations on the scale of 10 °K. Nielsen [2001b] showed that these variations have a statistically significant connection to the variation of the average strength of the stratospheric polar vortex at 10 hPa height, which corresponds well with the general circulation theory and the theory of gravity wave propagation. Some of these variations could also be caused by mesospheric planetary waves as discussed by Nielsen et al. [2001a]. No auroral effects have so far been reported / observed to influence the OH temperature series.
Another thing to consider is the seasonal variation with the direct impact of the solar UV flux on mesosphere. Sahai et al.  reported a correlation of 0.73 between the OH temperature over Calgary (51°N, 114°W) and the 10.7 cm solar flux using 87/88 and 1990 to represent seasons with low and high solar activity, respectively. The 1990 winter temperatures were approximately 30 °K warmer than the 87/88 temperatures, which is in the order of 10 °K higher than model calculations by Huang and Brasseur . At the moment we do not see any obvious connection in our data. Further studies are indeed needed to investigate the correlation between the direct influence of the solar flux and the mesospheric temperatures.
In order to determine the temperature trend with time, we have to correct for variation with solar cycle. As seen above this variation affects the annual variation and is mainly implied by the variation in strength of the stratospheric warmings with solar cycle. We know that the stratospheric warming is associated with mesospheric cooling. A set of criteria is necessary in order to separate the stratospheric - cooled mesosphere from the normal yearly variation:
A first cut would be to average the data for each year independent on whether stratospheric warmings occur or not. Figure 4 show both the monthly and daily temperatures. The overall seasonal or yearly mean winter temperature is calculated as an average between the monthly mean temperatures of December and January. A linear trend is calculated using the yearly averages, expect for the 80/81 season. This season was excluded since a criteria of minimum 14 days per season is applied. The trend is calculated using a polynomial fit of first degree with weights defined as number of days per season divided by the standard deviation for each the yearly average temperatures. This trend is compatible with zero: +0.03 °K/y ± 0.002 °K/y. The uncertainty of the trend estimate is calculated by the bootstrap method [Efron and Tibshirani, 1993]. If we do not apply weights, the trend becomes -0.4 °K/y ± 0.09 °K/y. Note that the errors in each point of the polynomial fits with and without weights are in the order of ±1°K and ±3°K, respectively. Our trend estimate is so far definitely not comparable to the large values up to -1 °K/y at mid latitudes reported by Golitsyn et al . The nearly zero temperature trend in the polar summer mesosphere reported by Lübken  seems to be more compatible with our results. We preliminarily conclude that there is no obvious temperature trend in the polar winter mesosphere over Svalbard.
The Auroral station in Adventdalen is owned by the University of Tromsø. The University Courses on Svalbard (UNIS) operates the station. We deeply appreciate the support of the Optical group at the Geophysical Institute, University of Alaska, who owns most of the instruments. The American program at the observatory is supported by the Atmospheric Science Division of the National Science Foundation.
Akmaev, R. A. and V. I., Fomichev, Cooling of the mesosphere and lower thermosphere due to doubling of CO2, Ann. Geophysicae, Vol. 16, 1501-1512, 1998.
Baker, D. J., and Stairs, Jr., A. T., Rocket measurements of the altitude distribution of the hydroxyl airglow, Physica Scripta, Vol. 37, 611 - 622, 1988.
Bates, D. R, and M., Nicolet, Atmospheric hydrogen, Publ. Ast. Soc. Pacif., Vol. 62, 106-110, 1950.
Bittner, M., Offermann, D., Graef, H. H., Mesopause temperature variability above a midlatitude station in Europe, J. Geophys. Res., Vol. 105, No. D2, 2045 - 2058, 2000.
Brasseur, G. and S., Solomon, Aeronomy of the middle atmosphere, D. Reidel Publishing Company, Dordrecht, Holland, 1986.
Conner, J., An Investigation of the Dynamics of the Mesopause. Fabry-Perot Observations of Winds and Temperatures from Nightglow Emissions, Ph. d. thesis, University of Alaska, Fairbanks, 1993
Duck, T. J., Whiteway, J. A., and Carswell, A. I., A detailed record of High Arctic middle atmospheric temperatures, Geophys. Res., Vol. 105, No. D18, 22.909 - 22.918, 2000.
Efron, B. and Tibshirani, R. T., An Introduction to the Bootstrap, Chapman & Hall, New York, 1993.
Fastie, W. G., A small plane grating monochromator, J. Opt. Soc. Am., Vol. 42, 641, 1952a.
Fastie, W. G., Image forming properties of the Ebert monochromator, J. Opt. Soc. Am., Vol. 42, 647, 1952b.
Fels, S. B., and M. D., Schwarzkopf, An efficient, accurate algorithm for calculation CO2 15 mm - band cooling rates, J. Geophys. Res., Vol. 86, 1205-1232, 1981.
Gadsden, M., A secular change in noctilucent cloud occurrence, J. Atmos. Terr. Phys., Vol. 52, 247-251, 1990.
Greet, P.A., W. J. R., French, G. B., Burns, P. F. B., Williams, R. P., Lowe, and K., Finlayson, OH(6-2) spectra and rotational temperature measurements at Davis, Antarctica, Ann. Geophysicae, Vol. 16, 77-89, 1998.
Golitsyn, G. S., A. L., Semenov, N. N., Shefov, L. M., Fishkova, E. V., Lysenko and S.P. Perov, Long-term temperature trends in the middle and upper atmosphere, Geophys. Res. Lett., Vol. 23, 1741-1744, 1996.
Hamilton, K., Comprehensive Meteorological Modeling of the Middle Atmosphere: a Tutorial Review, J. Atm. Terr. Phys., Vol. 58, No. 14, 1591-1627, 1996.
Herzberg, G., Spectra of Diatomic Molecules, D. Van Nostrad co., inc. New Y.,1950.
Herzberg, G., The atmospheres of the planets, J. Roy. Ast. Canada, Vol 45, 100-123, 1951.
Huang, T. Y. W. and Brasseur, G. P., Effect of long-term solar variability in a two-dimensional interactive model on the middle atmosphere, J. Geophys. Res., Vol. 98, No. 20, 413-420, 1993.
Krassovsky, V. I., N. N., Shefov and V. I., Yarin, Atlas of the Airglow Spectrum 3000-12,400 Å, Planet. Space Sci., Vol. 9, 883-915, 1962.
Kuhn, W. R., and J. London, Infrared radiative cooling in the middle atmosphere (30 - 110 km), J. Atmos. Sci., Vol. 26, 189-201, 1969.
London, J., In Proceedings of the Nato Advanced Institute on Atmospheric Ozone (Portugal), U.S. Dept. of Transportation, FAA - Washington, D.C., USA No. FAA-EE-80-20, 1980.
Labitzke, K., Naujokat, B., Lenschow, R., Petzoldt, K., Rajewski, B. and Wohlfart, R. C., The third winter of MAP*-dynamics, 1984/85; A winter with an extremely intense and early major warming, Beil. Berliner WettKarte, Vol. 67, 1985.
Labitzke, K., Stratospheric-mesospheric midwinter disturbances, A summary of characteristics, J. Geophys. Res., Vol. 86, 9665-9678, 1981.
Lübken, F.-J., K.-H., Fricke and M., Langer, Noctilucent clouds and the thermal structure near the Artic mesopause in summer, J. Geophys. Res., Vol. 101, 9489-9508, 1996.
Lübken, F.-J., Thermal structure of the Arctic summer mesosphere, J. Geophys. Res., 104, 9135 - 9149, 1999.
Lübken, F.-J., Nearly zero temperature trend in the polar summer mesosphere, Geophys. Res. Lett., 27, 3603 - 3606, 2000.
McKinley, J. D., D., Garvin, and M. J., Boudart, Production of excited hydroxyl radicals in the hydrogen atom-ozone reaction, J. Chem. Phys., Vol 23, 784-786, 1955.
Meinel, A. B., Hydride emission band in the spectrum of the night sky, Astrophys. J., Vol. 111, 207, 1950.
Mies, F. H., Calculated vibrational transition probabilities of OH(X 2 P), J. molec. Spectrosc., Vol. 118, 507, 1974.
Myrabø, H. K., C. S., Deehr and G. G., Sivjee, Large-amplitude nightglow OH(6-3) band intensity and rotational temperature variations during a 24-hour period at 78°N, J. Geophys. Res., Vol. 88, 9255, 1983.
Myrabø, H. K., Temperature variation at mesopause levels during winter solstice at 78°N, Planet Space Sci., Vol. 32, 249-255, 1984.
Myrabø, H. K., Winter-season mesopause and lower thermosphere temperatures in the northern polar region, Planet Space Sci., Vol. 34, 1023-1029, 1986.
Myrabø, H. K. and O. E., Harang, Temperatures and tides in the high-latitude mesopause region as observed in the OH night airglow emissions, J. Atmos. Terr. Phys., Vol. 50, 739-748, 1988.
Myrabø, H. K, C. S., Deehr, R., Viereck and K. Henriksen, Polar mesopause gravity wave activity in the sodium and hydroxyl night airglow, J. Geophys. Res., Vol. 92, 2527-2534, 1987.
Nielsen, K. P., Deehr, C. S., Raustein, E., Gjessing, Y. and Sigernes, F., Polar OH-airglow temperature variations in the 87/88 winter, Accepted by Physics and Chemistry of the Earth, 2001a.
Nielsen, K. P., Measurements and evaluation of airglow spectra over Svalbard, Master's thesis in meteorology, University of Bergen, Norway,
Rees, M. H., Physics and chemistry of the upper atmosphere, Cambridge University Press, Cambridge, Great Britain, 1989.
Sahai, Y., Giers, D. H., Cogger, L. L., Fagundes, P. R. and Garbe, G. P., Solar flux and seasonal variations of the mesospause temperatures at 51°N, J. Atm. Terr. Phys., Vol. 58, No. 16, 1927-1934, 1996.
Schmidlin, F. J., The inflatable sphere: A technique for the accurate measurements of the middle atmosphere temperatures, J. Geophys. Res., 96, 22.673 - 22.682, 1991.
Sigernes, F., Optical Studies of Proton Aurora, Dr. thesis, University of Tromsø, 1996.
Sigernes, F. and K. P., Nielsen, Computer Program: A synthetic hydroxyl spectrum generator for Windows, Yet another Arctic Research Consultant Ltd. [YaARCo], 2000. Download: http://fred.unis.no/YaARCo/Software/SyntheticOHZip.exe
Sivjee, G. G., Airglow hydroxyl emissions, Planet. Space Sci., Vol 40, 235-242, 1992.
Sivjee, G. G. and R. M., Hamwey, Temperature and chemistry of the polar mesopause OH, J. Geophys. Res., Vol. 92, 4663-4672, 1987.
Thomas, G. E., Global change in the mesosphere-lower thermosphere region: Has it already arrived?, J. Atmos. Terr. Phys., Vol. 58, 1629-1656, 1996.
Turnbull, D. N. and R. P., Lowe, New Hydroxyl Transition Probabilities and their Importance in Airglow Studies, Planet. Space Sci., Vol. 37, No. 6, 723-738, 1989.
Ulich, T., and Turunen, E., Evidence for log-term cooling of the upper atmosphere in ionosonde data, Geophys. Res. Lett., Vol. 24, No. 9, 1103 - 1106, 1997.
Viereck, R. A., Effects of Gravity Waves on the Polar O2 & OH Airglow, Thesis, University of Alaska Fairbanks, 1988.
Viereck, R. A. and C. S., Deehr, On the Interaction between gravity waves and the OH Meinel (6-2) and O2 atmospheric (0-1) bands in the polar night airglow, J. Geophys. Res., 94, 5397, 1989.
Viereck, R. A., A Review of Mesospheric Dynamics and Chemistry, U. S. National Report to International Union of Geodesy and Geophysics, Reviews of Geophysics 1987-1990, Supplement, 1131-1142, 1991.
Walterscheid, R. L., Sivjee, G. G., Schubert, G. and Hamwey, R. M., Large-amplitude semidiurnal temperature variations in the polar mesopause: evidence of a pseudotide, Nature, Vol. 324, 347 -349, 1986.