Characteristics and sources of gravity waves observed in NLC over Norway

Introduction Conclusions References


Introduction
Gravity waves contribute strongly to the global dynamics, circulation, structure, variability and thermal balance of the atmosphere.They can be generated in the lower atmosphere through topography, convection and wind shear, and propagate upward into the middle and upper atmosphere (Fritts and Alexander, 2003, and reference therein).Their amplitude grows with altitude due to the exponential decrease of the atmospheric density with height.In the mesosphere the waves can become unstable and break, dissipating their energy and depositing their momentum locally (Lindzen, 1981;Holton, 1983;Garcia and Solomon, 1985).spheric background wind modulates the gravity wave flux into the mesosphere and hence the momentum deposition.In winter time eastward propagating gravity waves are mostly blocked by the strong westerly winds in the stratosphere, and hence a net westward momentum will be deposited in the mesosphere that creates a poleward convergent and downward flow.At the summer pole, the easterly stratospheric winds lead to eastward momentum deposition and hence a divergent and upward flow.This establishes a pole-to-pole circulation in the mesosphere (Lindzen, 1981;Holton, 1982Holton, , 1983;;Garcia and Solomom, 1985) that drives the polar mesospheric temperatures adiabatically away from the radiative equilibrium and creates a cold summer and warm winter mesopause at high latitudes (Haurwitz, 1961;Garcia and Solomon, 1985;Frits and Luo, 1995;Luo et al., 1995;Espy and Stegman, 2002;Fritts and Alexander, 2003).The cold polar summer mesosphere reaches temperatures below 150 K, which is below the frost point for the small amount of water vapor that is present there.The resulting ice clouds that are formed are called Noctilucent Clouds (NLC) when viewed from the ground, or Polar Mesospheric Clouds (PMC) when viewed from space.
Since the thermal balance and the global circulation of the atmosphere are strongly determined by the gravity wave flux into the mesosphere, a correct representation of their effect is critical for global circulation models (GCMs).Longer scale gravity waves can be resolved explicitly in current GCMs.However, shorter-scale waves, which transport the major part of the momentum and energy to the mesopause (Vincent, 1984), have to be parameterized (Hamilton, 1996).For an adequate parameterization sufficient knowledge of seasonal and latitudinal behavior of gravity waves, as well as their source regions are essential (e.g.Rind et al., 1988;Hamilton, 1995), particularly in Polar Regions (Garcia and Boville, 1994;Hamilton, 1996).
Although there are a number of gravity-wave studies at high latitudes, most of these rely upon analysing the wave structures induced by the waves in the night airglow (e.g., Taylor and Henriksen, 1989;Pautet et al., 2005;Nielsen et al., 2006Nielsen et al., , 2009;;Espy et al., 2004Espy et al., , 2006)).However, these observations are not possible during the high latitude polar summer when the mesosphere remains sunlit.Recently, Pautet et al. (2011) rec-Introduction

Conclusions References
Tables Figures

Back Close
Full ognized that even though it is not possible to use the airglow during this period, the structures present in the NLC, known to be caused by gravity waves and their instabilities (Hines, 1960;Thomas, 1991;Fritts et al., 1993;Chandran et al., 2009Chandran et al., , 2010)), could themselves be used to infer the gravity waves present.Analysing NLC images from Stockholm, Sweden (59.5 • N, 18.2 • E), they were able to provide the first climatology of gravity-wave wavelengths, phase speeds and propagation directions from 60 to 64 • N during the polar summer.Similarly, using observations from the Aeronomy of Ice in the Mesosphere (AIM) satellite, Taylor et al. (2011, and references therein) have compiled an impressive climatology of gravity wave wavelengths and wave-front orientation from 75 to 85 • N and S for waves with wavelengths > 20 km.
We present here an extension of the Pautet et al. (2011) climatology to higher latitudes (64 to 74 • N) using NLC images taken from Trondheim, Norway (63.4 • N, 10.5 • E).In doing so, this study also fills the gap between the ground-based and satellite studies.In addition, we seek to expand those studies by characterizing the relative amplitude of the waves in the NLC and utilizing a ray-tracing model to identify the source regions of high latitude gravity waves present during the polar summer.Together, these observations provide a critical extension of the summertime gravity-wave climatology and information as to the efficiency of different sources required for parameterization of gravity-wave effects.

Observation and processing techniques of the NLC images
Four years of NLC pictures, since the summer of 2007, have been collected using an automatic digital camera (Canon EOS 450D) from Trondheim, Norway by the Leibniz Institute of Atmospheric Physics of Kühlungsborn, Germany (Baumgarten et al., 2009a).The camera takes pictures throughout the twilight period with exposure times of 0.5 s and 0.6 s with a cadence ranging between two to three frames per minute depending on the solar elevation angle.The camera views towards the north with a field of view (FOV) of 54.5 • × 41.5 • .Introduction

Conclusions References
Tables Figures

Back Close
Full  ).At this time, the Sun is in the northnorthwest, towards the bottom and left of the image.The NLC are clearly visible in the image when the solar Mie scattering from the clouds is brighter than the columnintegrated Rayleigh scattered sunlight from the atmosphere.However, in the lower part of the image, towards the Sun, the clouds are not visible due to the loss of contrast against the atmospheric Rayleigh scatter.Similarly, they are not visible in the upper part of the image as the Sun no longer illuminates the clouds.In those areas where the clouds are visible, different types of structures with various forms and sizes are clearly seen in the sunlit NLC.This is the result of short-period gravity waves modulating the ice-particle density, and the resulting intensity fluctuations trace out the waves in the centre and centre right of the image.In addition, very small scale structures are visible that may be the result of instabilities that develop as the gravity waves begin to break (Fritts et al., 1993).
As a first step in the image analysis, the NLC images are mapped to a horizontal plane through the NLC layer "the satellite view" (Witt, 1962;Pautet and Moreels, 2002).
In order to accurately project the NLC images, the viewing geometry of the camera, including the FOV, azimuth and elevation angles as well as rotation of the optical axis, were first calibrated using the visible stars that are manually identified in the images (Garcia et al., 1997).After the images were projected on a linear-scale grid using an optical ray tracer where the refractive index was calculated from air densities taken from MSISE-00.Each pixel of the camera was traced from the camera to the NLC altitude (Baumgarten et al., 2009b).Figure 2 shows image of the NLC picture shown in Fig. 1 projected onto a latitude and longitude grid assuming an NLC average height of 83 km (e.g.Witt, 1962;Fiedler et al., 2009).These images were used to identify gravity-wave events and measure the wavelength, horizontal phase velocity, direction of propagation and relative intensity fluctuation of each event.Introduction

Conclusions References
Tables Figures

Back Close
Full 3 The estimation of wave parameters

Estimation of wavelength and phase
The next step was the wave analysis where individual wave events in the NLC images were visually identified.Within a sequence of images, organized structures propagating in a uniform direction for at least 20 min were identified as gravity-wave events.Due to the short observational periods and challenging scattering conditions, a total of 34 such events were identified in this manner from the four summer observing seasons.
In order to infer the wave parameters, we have developed a robust method using a combination of FFT and a least-squares fitting of a sinusoidal function of the intensity data in order to quantify the horizontal wavelength, phase speed and direction of propagation of the waves in the NLC images.This method is different from the 2-D FFT method used by Pautet et al. (2011), originally developed by Garcia et al. (1997), to analyse the gravity waves from the Stockholm NLC images.Although this fitting method is more labour intensive and more subjective than the Pautet et al. (2011) method, given the low contrast of these high latitude observations it was necessary to manually identify the wave parameters.
Since the images are mapped to an equally spaced latitude and longitude grid, the raw pixel values in the projected image are first converted into kilometre space.Once a wave structure has been visually identified in the NLC image, a line was drawn normal to the wave front.To ensure that the chosen line of sight was truly normal to the wave front, a series of profiles were extracted by slightly changing the angle of the line.The frequency (wavelength) of the wave pattern along each of these lines was determined using a Fourier transform technique.The line with the smallest wavelength is taken to be normal to the wave front (shown in Fig. 2 as an arrow) and is along the propagation direction of the wave.This line was then used for the analysis of the wave structure.In cases where there was more than one wavelength present, the largest wavelength was selected for analysis.In some cases, where the images are blurred or contaminated, Introduction

Conclusions References
Tables Figures

Back Close
Full the FFT of the difference image, i.e. the image formed by subtracting two adjacent images, provided the wavelength.
To determine the phase velocity of the wave, several consecutive images, i.e. time steps, of the same wave structure were analysed.First a bias, linear slope and a sinusoid at the fixed, minimum wavelength determined above were fit to the intensity variations along the propagation direction.The bias and linear slope were then removed from the wave structure, and an FFT was applied to the de-trended intensity variations to determine the wavelength and phase of the wave structure along the propagation direction.The same procedure was carried out for all consecutive images to observe the change of the phase of the largest wavelength.It should be noted that the wavelength was determined only from the first image and held constant in all consecutive pictures.This is because other waves can propagate into the analysis regions and become the largest component of the wave structure.However, for the correct calculation of the phase values for each time step and hence the phase speeds of the wave it is important that the same wave is observed in all consecutive pictures.Figures 3 and 4 show an example of the Hovmöller diagrams of the raw data and the extracted wave structure along the propagation direction for consecutive time steps.

Extracting the amplitudes
In order to compare with gravity wave ray-tracing results, it is necessary to have an estimate of the relative amplitude of the waves present in the NLC.Particles can take hours to grow to visible sizes, but evaporation times approach seconds when temperatures are raised (Gadsden and Schröder, 1989).Thus, while the long term effect of short-period waves on NLC will be to reduce their visibility (Rapp et al., 2002), the bright and dark periodic structures in the NLC, with time scales of minutes, are primarily due to the density perturbations caused by the gravity wave rather than changes in the particle radius due to temperature effects.Thus, the relative amplitude of the waves was determined by the amplitude of the NLC intensity variations, determined above, divided by the average Mie scatter of the structured cloud.Introduction

Conclusions References
Tables Figures

Back Close
Full Here the background that was fitted to the raw image above has two components, the Rayleigh scatter from the background atmosphere plus the Mie scatter from the clouds.Due to the limited spatial scales of the NLC, the images were filtered to remove high spatial frequencies in order to determine the Rayleigh scatter component.Tests were done on nights with no apparent NLC to ensure that the pure Rayleigh scatter in these images was not affected by the filtering process.Nights with NLC present were filtered, and the Rayleigh-scattered component was compared to background nights at similar solar angles to ensure that the Mie scattering from the clouds was removed.Dividing the raw NLC image by the filtered "Rayleigh" image resulted in a flat background modulated only by the NLC wave perturbations.Along the direction of propagation, the background fitted to this ratio image represents the ratio of the average Mie plus Rayleigh scatter to Rayleigh scatter.This may be combined with the background from the raw image to obtain the average Mie scatter in the display being analysed.In this way, the relative amplitude of the NLC Mie scatter could be obtained.Since the wave structures are the result of density enhancements associated with the gravity waves, the relative amplitude of the Mie scatter was taken to be directly proportional to the relative density fluctuation caused by the wave.

Ray-tracing analysis
In order to examine the possible source regions for the gravity waves observed in our field of view, we used the Gravity Wave Regional or Global Ray Tracer (GROGRAT) model developed by Eckermann and Marks (1997).GROGRAT is a four-dimensional ray tracing algorithm that traces the evolution of a gravity wave's propagation and amplitude under the effect of background temperature and wind variations in the lower and middle atmosphere (Eckermann and Marks, 1995Marks, , 1997)).Here, we launch a spectrum of gravity waves from an array of latitude and longitudes points, and utilize GROGRAT to determine those waves that enter the FOV of the camera at Trondheim.Introduction

Conclusions References
Tables Figures

Back Close
Full In this study, the spectrum of waves suggested by Eckermann (1992), which was obtained from their experimental studies, is used.This spectrum consists of waves with uniform amplitudes of 0.3 m s −1 and discrete wavelengths of 12.5, 25, 50, 100, 200 and 400 km.Each wavelength is launched with phase speeds of 10, 20, 30, and 40 m s −1 at each 45 • of azimuth starting from the north.These waves are launched every 5 • of latitude between 60 • and 80 • N, and every 10 • of longitude between 0 • and 50 • E. The specific parameters entered into GROGRAT for each ray include an initial latitude, longitude, altitude, horizontal wave numbers and initial frequency.The background temperature and pressure fields used in GROGRAT were generated at 2 km altitude intervals from 0 to 80 km using the Mass Spectrometer-Incoherent Scatter Extended empirical model (MSISE-00) (Hedin, 1991) for the time of the display.Similarly, the wind fields were generated using the empirical Horizontal Wind Model HWM-93 (Hedin et al., 1996), for the time of the display.As the waves propagate upward to 80 km, the trajectories and spatial distributions of waves as a function of time is derived due to the interaction with these fields.
Source regions for gravity waves may lie near the tropopause (Fritts and Alexander, 2003) as in the case of convective activity or orographic lifting.However, they may be the result of secondary generation at higher altitudes by either wave breaking (Holton and Alexander, 1999;Vadas and Fritts, 2003) or as a result of the strong zonal wind shears in the mesosphere (Fritts, 2003).To test this, waves were launched from both tropopause heights, 5 km, as well as from 60 km at the same geographic grid points, and the relative amplitudes at NLC heights could be compared with those derived from the NLC image analysis to determine the most probable source altitude.

Results
The characteristic horizontal wave parameters such as wavelength, phase velocity, period, propagation direction and amplitude of the gravity waves were inferred from the NLC images, and a histogram of these wave parameters and the direction of propaga-Introduction

Conclusions References
Tables Figures

Back Close
Full The histogram of wave structures derived from the NLC images (Fig. 5a) shows that there is a strong preference for structures with scales less than 35 km, with an average wavelength of 24 km.However, the characteristic horizontal wavelengths for the ray-traced waves launched from both the 5 km and 60 km altitudes show that while these short wavelength waves do propagate to NLC heights, both distributions are similar and increase with wavelength.Thus the ray tracing shows that most of the waves reaching the mesosphere from either tropospheric or stratospheric sources have wavelengths greater than 40 km.This discrepancy is likely due to the spatial constraints of the camera.Due to the size of the field of view (53 • × 41 • ) of the NLC camera, the maximum observable horizontal wavelength is less than 100 km (Pautet et al., 2011).Thus, very long waves will not be observable as a modulation of the intensity within a single image and will not be identified as waves.Similarly, the high levels of Rayleigh scattering present in our high-latitude images further restrict our usable field of view and will bias us towards shorter period waves.
The observed phase speeds of the waves observed in the NLC (Fig. 5b) and those resulting from the wave tracing (Figs.6b and 7b) all average around 31 m s −1 .However, whereas the ray-traced waves from 5 km and 60 km both show a peak in the distribution near the average value, the waves observed in the NLC have a broader range of phase speeds, with most waves having a phase speed of between 10 and 20 m s −1 .Conversely, most of the waves detected in the NLC and those traced upward from both 5 km and 60 km (Figs.5c, 6c and 7c) have observed periods below 20 min.However, the NLC waves have nearly and equal percentage of waves with observed periods near 30-40 min, resulting in an average observed period of 24 min.In contrast, the waves traced upward from 5 km and 60 km predominantly occur with observed periods less than 5 min, with a steadily decreasing number of waves at longer observed periods, and few waves occurring with periods longer than 20 min.Introduction

Conclusions References
Tables Figures

Back Close
Full Perhaps the most prominent difference between the waves propagating up from the tropopause and those originating in the stratosphere are their relative amplitudes at NLC altitudes, shown in Figs.6d and 7d, respectively.We see in Fig. 6d that the waves originating near the tropopause attain relative amplitudes ranging between ∼1 and 8 %, whereas those initiated near 60 km, shown in Fig. 7d, only grow to between ∼0.1 and 0.3 %.The relative wave amplitudes observed in the NLC, shown in Fig. 5d, range from ∼1 to 10 %.Thus, the amplitude of the majority of waves observed in the NLC Mie scatter would be more consistent with sources near the tropopause.However, this is not to say that waves generated in the stratosphere do not perturb NLC.Rather, there is likely an observational bias to the large amplitude waves from lower altitudes that present the largest contrast and are easiest to detect.
The propagation directions of the ray-traced waves, shown in Fig. 8b and c, indicate that the waves originating from near the tropopause or from 60 km both have similar eastward directions when they reach NLC altitudes.While the waves observed in the NLC also show a tendency for eastward propagation, Fig. 8a shows that they tend to be directed more towards the north.Additionally a large percentage of the NLC waves are observed moving in a WSW direction (240-270 • azimuth) and towards the ENE (60-90 • azimuth) that the ray tracing has indicated are blocked.A null hypothesis t test was applied on the characteristics of these waves propagating in in these two "blocked" directions to determine whether their wave characteristics differ from the rest of the population.The results, shown in Table 1, indicate that the ENE waves are not significantly different from the general population, but that the WSW waves represent a unique population of waves with shorter wavelengths (λ avg = 14.4 km compared to 24 km) and higher phase speeds (c avg = 42.4 m s −1 compared to 31 m s −1 ).Fritts In order to examine possible source regions for the gravity waves observed in the NLC, the grid points from which the waves were launched at 5 km and 60 km are shown on a map in Fig. 9a and b, respectively.The size of the symbol indicates the percentage of waves originating from a particular grid point that enter the field of view of the NLC imager at Trondheim.The geographical distribution of the wave sources at the two altitudes depict that the coastal regions of Norway at the latitude of 70 • N are the primary source regions of gravity waves propagating upward from 5 km.The waves propagating upward from 60 km represent a broader source region that extends eastward and southward of the camera field of view, but, as indicated above, these would cause smaller perturbations in the NLC.

Discussion
We can compare these climatological results over Norway with the previous work of Pautet et al. (2011), which used imager data from Sweden from 2004 to 2007 that extended to the southern limit of our observations.The comparison of the average horizontal wavelength of the waves shows reasonable agreement, with a mean value of 24 km and 25 km over ∼69 • N and ∼62 • N, respectively.Thus, the wave-like perturbations derived from the NLC images in the summer mesosphere are similar in scale from 60 • to 74 • N. In addition, they are similar to mesospheric short-period GWs observed in the air-glow emissions at middle (Taylor et al., 1998), polar (Nielsen et al., 2006) and equatorial (Nakamura et al., 2003)  studies of gravity waves observed at slightly higher altitude in the airglow emission.For example, Taylor et al. (1997) and Nakamura et al. (2003) found an average value of approximately 48 m s −1 at low latitudes (2 • S), and Hecht et al. (2001) found a mean value of 50 m s −1 at middle latitudes (40 • N).Similarly, the high latitude measurements in the Southern Hemisphere at Halley (76 • S) found a mean value of 48 m s −1 (Nielsen et al., 2009) during winter.Correspondingly, the observed periods found in these studies are typically on the order of 8-10 min.However, Ejiri et al. (2003) found an average value of 35 m s −1 at mid latitude ( 40• N), similar to this work.This wide range of values could be due to different GW sources, for example orographic vs. frontal generation, where the observed phase speeds would be Doppler shifted due to speed of the source relative to the observer.The directions of propagation of the waves observed in the NLC show that the majority of the waves in Norway are propagating towards the north and northeast at azimuths between 0 • and 30 • , but there is a separate population of small-scale fast waves propagating towards the southwest between 240 • and 270 • of azimuth.Similarly, Pautet et al. (2011) found that 60 % of the waves in Sweden propagate towards the north or the northeast directions, but did not observe the waves moving towards the southwest.This matches the propagation directions found for the waves launched from 5 km and 60 km, where the SW direction is blocked by the wind field.Given this blocking, and the fact that the waves we observe moving towards the SW have significantly shorter wavelengths and higher phase speeds than those propagating northward or north-eastward, indicate that these are formed as a result of local, transverse instabilities accompanying wave breaking (Fritts et al., 2002).The shorter observing times available each night at Trondheim as compared to Stockholm, as well as the more restricted field of view resulting from the higher levels of Rayleigh scatter means that these fast moving, short scale instabilities were more likely to have been picked for analysis.Introduction

Conclusions References
Tables Figures

Back Close
Full

Conclusions
The climatology of summertime gravity-waves detected in NLC between 64 and 74 • N is similar to that observed between 60 and 64 • N by Pautet et al. (2011).While the average horizontal wavelengths observed in these two studies agree with previous results from airglow imaging, the observed phase speeds are much slower, and the observed periods correspondingly longer, than those observed in the airglow at lower latitudes and during high-latitude winter.This could be due to differences in the gravitywave sources, or an observational bias brought on by the more restrictive fields of view available to the NLC cameras when compared to the nearly all-sky airglow observations.The gravity-waves observed in the NLC north of 64 • N continue to propagate towards the north or northeast as they have been observed to do south of 64 • N.However, the unique population of fast, short wavelength waves propagating towards the SW is consistent with transverse instabilities generated in-situ by breaking gravity waves.Finally, comparison of the relative amplitude of the waves observed in the NLC Mie-scatter with ray-tracing results show that waves propagating from the tropopause, rather than those resulting from secondary generation in the stratosphere or mesosphere, are more likely to be the sources of the prominent wave structures observed in the NLC.The geographic distribution of the waves generated near the tropopause that can be ray traced into the observational field of view identify the coastal region of Norway along the latitude of 70 • N as the primary source region.As this is the region where many of the polar low pressure systems make landfall, interaction of the frontal systems with the coastal mountainous terrain is a likely generation mechanism for the prominent wave structures observed in the NLC near the mesopause.Introduction

Conclusions References
Tables Figures

Back Close
Full Filtering of wave energy by the strato-Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Figure1is an example of the NLC image taken on the night of16-17 July 2008, at  22:30 UT, which is 23:12 Local Solar Time (LST).At this time, the Sun is in the northnorthwest, towards the bottom and left of the image.The NLC are clearly visible in the image when the solar Mie scattering from the clouds is brighter than the columnintegrated Rayleigh scattered sunlight from the atmosphere.However, in the lower part Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | tion of the waves are shown in Figs. 5 and 8, respectively.For comparison, histograms of the wave parameters from the ray-traced waves that reached the NLC altitude from the source region of 5 km and 60 km are shown in Figs. 6 and 7, respectively.
Discussion Paper | Discussion Paper | Discussion Paper | et al. (2002) have shown that these fast, short period waves are associated with local instabilities generated during wave breaking.These instabilities tend to move along the wave front and transverse to the general wave propagation direction, similar to what is observed here.Discussion Paper | Discussion Paper | Discussion Paper | latitudes, where peak horizontal wavelengths around 25-30 km are found.The average observed phase speed for the Trondheim data was found to be 31 m s −1 , which also compares favourably with the average value of 27 m s −1 found by Pautet et al. (2011) for the summertime data taken from Stockholm.Similarly, there is no significant difference in the average observed periods of the two data sets: 24 min at ∼69 • N and 22 min at ∼62 • N.However, these observed phase speeds measured using NLC data during the high latitude summer are smaller than other previous long-term Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |