Gravity wave instability structures and turbulence from more than one and a half years of OH* airglow imager observations in Slovenia

We analysed 286 nights of data from the OH* airglow imager FAIM 3 (Fast Airglow IMager) acquired at Otlica Observatory (45.93 °N, 13.91 °E), Slovenia between 26 October 2017 and 6 June 2019. Measurements have been performed 10 with a spatial resolution of 24 m pixel and a temporal resolution of 2.8 s. A two-dimensional Fast Fourier transform is applied to the image data to derive horizontal wavelengths between 48 m and 4.5 km in the upper mesosphere / lower thermosphere (UMLT) region. In contrast to the statistics of larger scale gravity waves (horizontal wavelength up to ca. 50 km, e.g., Hannawald et al., 2019) we find a more isotropic distribution of directions of propagation, pointing to the presence of wave structures created above the stratospheric wind fields. A weak seasonal tendency 15 of a majority of waves propagating eastward during winter may be due instability features from breaking secondary gravity waves that were created in the stratosphere. We also observe an increased southward propagation during summer, which we interpret as an enhanced contribution of secondary gravity waves created as a consequence of primary wave filtering by the meridional mesospheric circulation. We present multiple observations of turbulence episodes captured by our high-resolution airglow imager and estimated the 20 energy dissipation rate in the UMLT from image sequences in 25 cases. Values range around 0.08 and 9.03 W kg and are on average higher than those in recent literature. The values found here would lead to an approximated localized maximum heating of 0.03 3.02 K per turbulence event. These are in the same range as the daily chemical heating rates for the entire atmosphere reported by Marsh (2011), which apparently stresses the importance of dynamical energy conversion in the UMLT.


Introduction
Fully understanding the contribution of gravity waves to atmospheric dynamics is still a major issue when establishing climate models. Due to the various sources and mechanisms of interactions the effects of gravity waves have to be represented in these models using advanced parameterizations (Lindzen, 1981;Holton, 1983;de la Cámara et al., 2016) to cover as many aspects as it is possible given the restricted model resolution. Gravity waves exist on a large span of time scales ranging from several 30 hours down to the Brunt-Väisälä (BV) period, which corresponds to ca. 4 -5 min in the upper mesosphere / lower thermosphere (UMLT) region (Wüst et al., 2017b) and represents the smallest possible period of gravity waves. They show diverse behaviour depending strongly on wave properties like their periodicity (Fritts & Alexander, 2003;Beldon & Mitchell, 2009;Hoffmann et al., 2010;Wüst et al., 2016;Sedlak et al., 2020), which makes it even harder to fully account for them by means of parameterization. Furthermore, gravity wave generation is not restricted to the troposphere but can also take place at higher 35 altitudes, such as secondary wave excitation due to breaking gravity waves (see, e.g., Holton & Alexander, 1999;Satomura & Sato, 1999;Vadas & Fritts, 2001;. As Fritts & Alexander (2003) state, it is necessary to metrologically capture all parts of the gravity wave spectrum. This includes especially dynamics on short scales where gravity wave breaking is induced by the development of instabilities. One of the most prominent features in this context is the formation of Kelvin-Helmholtz instability (KHI), which occurs as a 40 consequence of a dynamically instable atmosphere due to wind shear (Browning, 1971). Gravity wave instability can also be of convective nature when growing wave amplitudes lead to a superadiabatic lapse rate (Fritts & Alexander, 2003). In general, atmospheric instabilities like KHIs often manifest as so-called ripplesperiodic structures with small spatial dimensions and short lifetimes (Peterson, 1979;Adams et al., 1988;Taylor & Hapgood, 1990;Li et al., 2017).
Gravity wave breaking and the conversion of the transported energy into heat takes place in the course of turbulence. Once a 45 wave breaks and motion shifts from laminar to turbulent flow energy is cascaded to smaller and smaller structures until viscosity becomes dominant over inertia and energy is dissipated into the atmosphere by viscous damping (see, e.g., Lübken et al., 1987).
The process of turbulence manifests as formation of vortices, so-called eddies. They cause turbulent mixing of the medium, resulting in the dissipation of turbulent energy at an energy dissipation rate . According to the theory of stratified turbulence, 50 depends on the characteristic length scale and velocity scale of the turbulent features. The energy dissipation rate is then given by = 3 (1) (see, e.g., Chau et al, 2020; who apply this equation to radar observations of KHIs). is a constant which is found to be equal and there are only few values of available at UMLT heights. Lübken et al. (1997) use rocket measurements to retrieve in the height range 65 -120 km. Baumgarten & Fritts (2014) use imaging techniques of mesospheric noctilucent clouds to 60 investigate the formation of KHIs and the onset of turbulence.
At the same height, remote sensing measurements of the OH* airglow are an established access to UMLT dynamics. The OH* airglow is a layer at an average altitude of ca. 86-87 km with a full width at half maximum (FWHM) of ca. 8 km (Baker & Stair, 1988;Liu & Shepherd, 2006;Wüst et al., 2017b). Remote sensing techniques include spectroscopic measurements of strong emission lines and the analysis of temperature time series derived from these (Hines & Tarasick, 1987;Mulligan et al., 65 1995;Bittner et al., 2000;Reisin & Scheer, 2001;Espy & Stegman, 2002;Espy et al., 2003;French & Burns, 2004;Schmidt et al., 2013, Wachter et al., 2015Silber et al., 2016, 2017a, but also twodimensional imaging in the short-wave infrared (SWIR) range (see, e.g., Peterson & Kieffaber, 1973;Hecht et al., 1997;Taylor, 1997;Moreels et al., 2008;Li et al., 2011;Pautet et al., 2014;Hannawald et al., 2016Hannawald et al., , 2019Sedlak et al., 2016;Wüst et al., 2019 and many more). 70 The technology of OH* imaging has undergone a rapid technical progress over the last few decades. Improvements in sensor technology and optics have provided the possibility to observe the signatures of gravity waves that manifest as periodic brightness variations in infrared images of the OH* airglow layer. The observations range from all-sky imaging of large-scale gravity waves (e.g., Taylor et al., 1997;Smith et al., 2009) to high resolution images of smaller gravity waves (Nakamura et al., 1999) and their breaking processes (Hecht et al., 2014;Hannawald et al., 2016). Hannawald et al. (2016) use an airglow 75 imager called FAIM (Fast Airglow IMager) that is well-suited for the observation of small-scale gravity waves with a high temporal resolution of 0.5 s. Based on three years of continuous night-time observations at two different Alpine locations Hannawald et al. (2019) show statistics of gravity wave propagation for waves with horizontal wavelengths smaller than 50 km based on data of the same kind of instrument.
In 2016 we put into operation another FAIM instrument (FAIM 3) which still has a high temporal resolution of 2.8 s, but also 80 a high spatial resolution of up to 17 m pixel -1 (measurements in zenith direction utilizing a 100 mm SWIR objective lens). We were not only able to observe wave patterns on extraordinary small scales (smallest horizontal wavelength 550 m) but also the formation of a vortex which we interprete as the turbulent breakdown of a wave front (Sedlak et al., 2016).
From October 2017 to June 2019 the instrument observed the area around the Gulf of Trieste from Otlica Observatory, Slovenia (45.93 N, 13.91 °E), which is a partner observatory within the context of the Virtual Alpine Observatory (VAO; 85 https://www.vao.bayern.de). This larger data basis includes further observations of small-scale wave features and turbulence which are investigated here.
The focus of this paper is on analysing small-scale dynamics in the UMLT region in FAIM 3 images with regard to two aspects: 4 2. The estimation of dissipated energy by analyzing multiple episodes of turbulence (such as the one exemplarily presented in Sedlak et al., 2016).

Instrumentation 95
FAIM 3 is an OH* airglow imager that has been put into operation in February 2016 at the German Aerospace Center (DLR) in Oberpfaffenhofen, Germany. It consists of the SWIR camera CHEETAH CL manufactured by Xenics nv, which has a thermodynamically cooled 640 x 512 pixels InGaAs sensor array (pixel size 20 µm x 20 µm, operating temperature 233 K). The camera is sensitive to electro-magnetic radiation in the wavelength range from 0.9 to 1.7 µm (for further technical details see Sedlak et al., 2016). equipped with a 100 mm SWIR lens by Edmund Optics ® with aperture angles of 7.3 ° and 5.9 ° in horizontal and vertical direction. Neglecting the curvature of the Earth, this configuration leads to a trapezium-shaped field of view (FOV) with a size of ca. 182 km² (13.1 -14.1 km x 13.4 km) at the mean peak emission height of the OH* layer at ca. 87 km. The mean spatial resolution is therefore 24 m pixel -1 . Due to the abovementioned measurement geometry the FOV is located above the Gulf of Trieste. The integration time of FAIM 3 is 2.8 s, which leads, depending on the season, to the acquisition of ca. 10,000 to 110 18,000 images per night.

Data Basis
All in all, in 477 nights image data were acquired by FAIM 3 at OTL. Since OH* airglow observations are only possible under clear sky conditions, cloudy episodes are sorted out by analyzing keograms. This yields 410 clear sky episodes (durations between 20 min and 13 h) that are distributed over 286 measurement nights. Thus, ca. 60 % of the acquired nights at OTL 115 include suitable OH* observations. Before being analysed, the images undergo the same preprocessing steps as in Hannawald et al. (2016Hannawald et al. ( , 2019 and Sedlak et al. (2016): a flat-field correction is performed and the images are transformed to an equidistant grid, which corresponds to a trapezium-shaped FOV due to the inclination from zenith. For each episode the average image is subtracted to ensure that all remnants of fixed patterns are removed (e.g., reflections of the objective lens in the laboratory window during bright nights). 120 Due to the small FOV of FAIM 3 we renounce the application of a star removal algorithm to avoid an interpolation of too many pixels. In order to extract periodic signatures a two-dimensional Fast Fourier Transform (2d-FFT) is applied to squared cut-outs of each image, so that neither dimension is favored by the analysis. These cut-outs were chosen to have a side length of 406 pixels (equals ca. 9.7 km) as this is the largest possible square fitting into the transformed images. The 2d-FFT is performed on the squared image cut-out as described by Hannawald et al. (2019). A fitted linear intensity gradient is subtracted 125 from the input images and a Hann window is applied during the 2d-FFT to reduce leakage effects. A local maximum filter is applied to automatically find peaks in the spectra and thus plane wave structures, which allows identifying and analyzing single wave events. Zero-padding on the images (to a size of 2160 x 2160 pixels) is used to improve this identification of peaks in the spectra. Hannawald et al. (2019) present a statistical analysis of gravity waves with horizontal wavelengths between 2 and 62 km (with focus on waves with horizontal wavelengths larger than 15 km). With FAIM 3 having a smaller FOV and a higher 130 spatial resolution than the FAIM instrument used therein, we are now able to present statistics of gravity wave parameters that tie in almost seamlessly with the statistics of longer-scale waves of Hannawald et al. (2019): due to the spatial resolution and the FOV size we cover the horizontal wavelength range from 48 m to 4.5 km. Wave structures with horizontal wavelengths of half the FOV size still showed a strong bias toward phases 0 or . Extensive testing showed that this effect disappeared when lowering the upper wavelength limit to 4.5 km. 135 Observed wave structures have to meet several quality criteria in order to be considered a wave event. A wave structure has to be present for at least 20 s and has to be found in at least eight images. This is in contrast to Hannawald et al. (2019) who demand wave signatures to be present for at least 120 s and to appear in at least 100 images within this episode, stating that these restrictions specifically filter out many transient and small-scale wave features as they want to focus on larger persistent waves. 140 Furthermore, FAIM 3 wave events are considered if they have an amplitude of at least 25 % of the maximum observed wave amplitude. Wave structures with this amplitude can just be recognized in the image by the eye. Demanding all the quality criteria mentioned above a total number of 5697 wave events remains. Further restricting these criteria has not significantly altered the distributions of the wave parameters that are presented in the following. An exemplary event and the respective 2dimensional spectrum are shown in Figure 1. 145 We often observe episodes of turbulence in our image series that exhibit the typical dynamics of vortex formation and quasichaotic behavior. While the identification of wave structures is done automatically by the 2d-FFT, finding turbulent vortices is done by hand. Turbulent eddy formation can be well recognized by eye when viewing the episodes in the dynamical course of a video sequence. However, the combined effect of these vortices having a certain variety of shapes and sizes, being almost invisible in single images without comparison to preceding or successive images, and causing (compared to other features such 150 as wave fronts) rather small brightness fluctuations in the images hampers strongly the application of image recognition algorithms. For the given data basis 25 episodes of turbulence with sufficient quality to derive turbulence parameters are found.
The dates along with the respective turbulence parameters are summarized in Table 1.

Statistics of Wave Parameters 155
The wave statistics are presented in Figure 2 and 3. Please note that we a using the word 'wave' for all wave-like structures we find in the images. The question whether these are actual gravity waves is discussed in section 5.
Wave periods range from 21 s to 1498 s (25 min). The median wave period is 359 s (6 min Events with meridional phase speed zero have not been considered for the mean values.

Wave Dissipation
To give an impression of the turbulent dynamics we observe, we present four of our turbulence episodes as video supplement.
On 16 November 2017, 02:16 UTC, the turbulent breakdown of parts of an extended wave field can be observed (video 1). On 170 6 December 2017, 00:26 UTC, several fronts seem to be building up and form rotating vortices (video 2). This can be observed even clearer on 14 October 2018, 17:08 UTC, where the residual movement of turbulent features can be well recognized above the general background movement (video 3). On 4 November 2018, 19:18 UTC, breaking wave fronts seem to form rotating structures of nearly cylindrical shape, while these are accompanied by other turbulently moving eddies (video 4).
We estimate the turbulent energy dissipation rate using equation (1). However, in contrast to Chau et al. (2020) who used 175 radar measurements, we only have horizontal information from our airglow imager. Hecht et al. (2021) demonstrate an approach how to apply equation (1) to purely horizontal airglow imager data, which we adapt to our observations in the following. The characteristic length scale can be read from the images by measuring the size of the turbulent features. The velocity scale is given by the residual velocity of these features. In our observations, they are part of larger instability features, which we assume to be advected by the background wind. We determine by reading the actual velocity of the 180 turbulent features and subtracting the background movement in the resulting direction. This is exemplarily shown in Figure   4. The two patches highlighted therein are both moving to the upper right direction but are approaching each other. This helps distinguishing background and residual movement.
As stated in section 3 we found 25 episodes of turbulence that allowed the derivation of and . Using equation (1), the energy dissipation rate is then calculated by = 3 . The resulting values are shown in Figure . We assume a general read-out 185 error of ±3 pixels, which corresponds to a distance of ±72 m. Velocities are determined by reading the distance a feature covers within an episode of at least ten images, which corresponds to a time span of 28 s. Thus, velocities are estimated with an error of ±2.6 m s -1 . The arising uncertainties of are calculated following the rules of error propagation.
The values of range from 0.08 to 9.03 W kg -1 . The median value is 1.45 W kg -1 .
Assuming the duration of dissipation being equal to the lifetime of the vortex the energy dissipation rate can be converted into 190 the amount of dissipated energy per mass. This is only done for those vortices that both form and decay within the FOV. The time intervals of dissipation are between 241 s and 922 s (4.0 -15.4 min) and can also be found in Table 1. Events are labelled as 'out of FOV' or as 'clouds' if either the formation or the decay of the vortex cannot be observed. No further analysis is performed for these events.
Multiplying energy dissipation rate and duration of dissipation equals the energy per mass that is released in the turbulent 195 process. We retrieve values between 30 and 3015 J kg -1 . Given that the released energy is entirely converted into heat we can make a rough estimate of the resulting temperature change by assuming isobaric conditions (may be approximately fulfilled due to the stable stratification of the atmosphere and small vertical dimension of eddies) and dividing energy per mass by the specific heat capacity of dry air (10 3 J K -1 kg -1 ). The resulting temperature changes in this work are in the range 0.03 -3.02 K (see Figure 6). All values of dissipated energy per mass and maximum temperature change can be found in Table 1. 200 Since we now have a time series of energy dissipation rate we can compare them to gravity wave activity in the UMLT above OTL. Parallel to FAIM 3, SWIR spectrometers called GRIPS (GRound-based Infrared P-branch Spectrometer) instruments deliver time series of OH* rotational temperatures derived from the OH(3-1) P-branch (1.5 µm-1.6 µm) at an initial temporal resolution of 15 s. Unlike the general instrument details discussed by Schmidt et al. (2013), the GRIPS 9 at OTL has a reduced aperture angle of 6.2 ° FWHM increasing its responsivity to smaller structures. As described in Sedlak et al. (2020), gravity 205 wave activitythe so-called significant wavelet intensity (SWI)for the periods 6 -480 min (period resolution 1 min) can be calculated by applying a wavelet analysis to these temperature time series. The FOV of GRIPS 9 is also located above the Gulf of Trieste and in ca. 30 km distance to the FAIM 3 FOV and has a size of approximately 13 km x 19 km. Since the spectroscopic observations are averaged over the entire FOV, GRIPS is most sensitive for gravity waves with horizontal wavelengths of several hundreds of kilometres . The time series of nocturnal SWI is restricted to those nights that exhibited 210 at least one of the turbulence episodes presented above. For each gravity wave period between 6 and 480 min (1 min steps), the correlation between the SWI at the respective period and the energy dissipation rate has been calculated. If there are observations of more than one vortex during one night, the respective energy dissipation rates are averaged to their mean value.
The Pearson correlation coefficient and the P value (significance test) are presented in Figure 7. We find almost no significant correlation for any gravity wave period. Long-periodic SWI (periods > 400 min) shows a slight positive correlation with the As can be seen in Figure 3, the wave structures we observed exhibit multiple directions. The strong tendency to north-eastern direction in summer and to the (south-)west in winter as observed by Hannawald et al. (2019) for medium-scale gravity waves cannot be confirmed for the waves observed here. However, slight tendencies are apparent in Figure 3. The north-western 220 component these authors observed during winter at Mt. Sonnblick in Austria with the FOV being positioned north of the Alps also appears in our data during autumn, winter and spring. During summer we find a conspicuous majority of waves propagating into southern direction.
The number of waves propagating eastward and westward is almost equal for the entire data set. However, as stated in section 4.1, more waves are oriented in eastward direction during winter whereas zonal directions are quite balanced during summer. 225 Although the eastward tendency during winter is quite weak, it contradicts the distribution that is expected for gravity waves being created in the lower atmosphere and propagating upward, being subdued to tropospheric and stratospheric wind filtering.
The eastward oriented mean wind profile during winter would lead to mainly westward propagating gravity waves reaching the UMLT without encountering critical levels. During summer the stratospheric winds reverse to westward direction, so that eastward oriented gravity waves are filtered in the tropopause and westward oriented gravity waves are filtered in the 230 stratosphere (see, e.g., Hoffmann et al., 2010;Hannawald et al., 2019).
As we have no accompanying wind measurements in the height of our observations it is difficult to decide by means of the period whether the wave structures presented in section 4.1 are small-scale gravity waves or instability features. Ca. 63 % of the wave events have an observed period above the BV period (here we used the climatology presented by Wüst et al., 2020), however these could at least in parts also be Doppler-shifted instability features instead of gravity waves. While the distinction 235 between largely extended wave-fields (bands) and small localized wave structures that are related to instability (ripples) is often made at a horizontal wavelength of 10 -20 km (Taylor et al., 1997;Nakamura et al., 1999), Li et al. (2017) remark that even structures with horizontal wavelengths of 5 -10 km may sometimes be gravity waves rather than instability features.
If this would be true for our small-scale wave structures, they might rather be secondary gravity waves (see, e.g., , being generated at greater heights by breaking gravity waves. Secondary gravity waves can either have larger 240 wavelengths and phase speeds than the primary wave if they are created by localized momentum deposition  or smaller wavelengths and phase speeds if they are induced by the nonlinear flow (wave-mean flow and wave-wave interactions; see, e.g., Bacmeister & Schoeberl, 1989;Franke & Robinson, 1999;Bossert et al., 2017). The former type of secondary gravity waves exhibits a rather broad spectrum of wave parameters with horizontal wavelengths longer than 500 km and horizontal phase speeds between 50 and 250 m s -1 , resulting in periods longer than ca. 30 min. The 245 wave structures found in this work have smaller horizontal wavelengths, phase speeds and periods and could therefore be more likely related to the latter type of secondary waves created by nonlinearities. However, these small-scale secondary waves are unlikely to propagate large vertical distances due to their small horizontal phase speeds . They have to be generated at even higher altitudes, i.e. close to the mesopause, to be observable with OH* airglow imagers. Hannawald in the mesosphere might play a vital role in filtering gravity waves. The meridional mesospheric circulation is oriented southward during summer and northward during winter, being much stronger during summer with ca. 10 -14 m / s (Yuan et al., 2008). Simulations by  show that advection by the background wind determines the direction of a newly created secondary wave. Based on these aspects, the accumulation of southward oriented waves we observe during summer could be a hint for gravity waves being filtered by the mesospheric circulation and generating subsequent secondary 255 waves with shorter wavelengths and periods, that are provided with a southward phase speed due to advection. This theory is also in good agreement with our observed meridional phase speeds: in the abovementioned velocity range of the summerly meridional mesospheric circulation (10 -14 m / s) meridional phase speeds are southward in 71 % of cases.
However, regarding the small horizontal wavelengths below 4.5 km, it is more likely that the major part of the observations presented in section 4.1 are related to instability features. The quite slow phase speeds (mean value 13.3 m /s ) are one hint for 260 this as typical gravity wave phase speeds accumulate around 40 m / s (see, e.g., Wachter et al., 2015 andWüst et al., 2018). If Figure 2b was the phase speed distribution of gravity waves, it is likely that a majority of them would encounter critical levels somewhere and would not be observable in the OH* layer. The small spatial scales of the wave structures we observe are typical for ripple structures as they were already observed with FAIM 3 (Sedlak et al., 2016). Their short life spans are not excluded by our quality criteria. Tuan et al. (1979) state that oscillations of this type are usually excited at periods of 4 -10 min, 265 which would explain the large number of wave events we observe in this period range. Observing ripple structures, it would not be surprising to obtain a certain diversity of directions of propagation. In principle, ripples originating from convective instabilities tend to be aligned perpendicular to the wave fronts of the initial wave, whereas ripples arising from dynamic instabilities form parallel to the initial wave fronts (Andreassen et al., 1994;Fritts et al., 1997;Hecht et al., 2000). However, it has been reported that ripples can be rotated by the background wind and that ripples may even be created by a combination 270 of both dynamical and convective instability (Fritts et al., 1996;Hecht, 2004). Considering the fact that the directional peculiarities of our observed wave events fit well with the expected behavior of secondary gravity waves, as discussed above, support the scenario of the wave structures being ripples from dynamic instabilities of secondary gravity waves, that originate from the stratospheric and mesospheric jet. Capturing structures related to instability is not unlikely, considering the numerous observations of turbulent vortices with the FAIM 3 setup. 275 Nevertheless, height-resolved measurements of the horizontal wind would be needed to determine the local wind shear and make a profound statement about atmospheric instability.
It has to be kept in mind that a 2d-FFT was used. Thus, periodic structures are assumed to be stationary, i.e., they extend over the entire image. Faint structures that appear only in small parts of the image (as does for example the 550m wave packet in Sedlak et al., 2016;Fig. 2) would be underrepresented by this analysis. 280 Measuring the energy dissipation rate in the UMLT is still challenging and there are only few studies yet. Rocket measurements of Lübken (1997) deliver energy dissipation rates between ca. 0.01 and 0.1W kg -1 between 85 and 90 km height at high latitudes. Chau et al. (2020) find an energy dissipation rate of 1.125 W kg -1 for their KHI event observed in the summer mesopause and state that this a rather high value compared to the findings of Lübken et al. (2002). Hocking (1999) provides a rescaled overview of earlier values of the energy dissipation rate and these have a maximum magnitude of 0.1 W kg -1 . Hecht et al. (2021) derive 285 a value of 0.97 W kg -1 from airglow images of a KHI event. Ranging from 0.08 up to 9.03 W kg -1 the values of energy dissipation rate derived here are higher than reported by other studies. However, the median value of 1.45 W kg -1 is not too far away from the values of Chau et al. (2020) and Hecht et al. (2021). The vortices we observe do not necessarily mark the smallscale end of the energy cascade. It could be possible that the energy is cascaded further to a larger number of smaller eddies that are no longer visible to our instrument. Parallel in-situ measurements (e.g., lidar, rockets) could be used to estimate the 290 significance of this effect. Additionally, it has to be kept in mind thatexcept for the studies of Hecht et al. (2021), whose value is quite similar to the median value of our datathe values compared here arise from different measurement techniques with different horizontal, vertical and temporal resolutions, so that the accessible scales are not necessarily identical due to the observational filter effect.
The derivation of the turbulence parameters performed here is challenging due to the blurred shape of dynamic signatures in 295 the OH* layer. The length scale and velocity scale of turbulent features have been extracted manually by measuring distances in the images and calculating distances from pixel values. We tried to quantify the read-out error by providing a measurement uncertainty and minimize it by repeating the analysis workflow on the same data multiple times. However, using equation (1) velocity dominates the length scale due to its power of 3, so that strongly depends on a parameter, which is quite difficult to extract from the images. All in all, it seems possible to derive turbulence parameters like the energy dissipation rate from high 300 resolution imager data.
The values of energy dissipation rate derived here show no significant correlation with gravity wave activity in the period range 6 -480 min. Turbulence thus can hardly be related to distinct periods of the gravity wave spectrum with the here-presented data. The slight positive correlation with gravity wave activity at periods larger than 400min may point to a special contribution of long-period gravity waves to the turbulence events we observe. However, this remains speculative at the current stage of 305 research, since this correlation is beyond the level of significance. A larger data basis of turbulence parameters and especially observations of period-resolved gravity wave activity at altitudes below will be needed to answer the question if all parts of the gravity wave spectrum drive turbulence generation in the UMLT equally.
Assuming that the turbulently dissipated energy is entirely converted into heat we find temperature changes of 0.03 -3.02 K that occur within time spans of 4.0 -15.4 min. Marsh (2011) report chemical heating rates in the atmosphere to be around 310 We present an analysis of small-scale dynamics of instability features and turbulence from OH* imager data acquired between 315 26 October 2017 and 6 June 2019 at Otlica Observatory, Slovenia. Measurements have been performed with the imager FAIM 3, which has a spatial resolution of ca. 24 m pixel -1 and a temporal resolution of 2.8 s.
Wave-like structures in the images are systematically identified by applying a 2d-FFT to nocturnal image sequences during clear sky episodes. All events meeting our persistency criteria were used to derive a statistical analysis of wave-like structures with horizontal wavelengths between 48 m and 4.5 km. The small horizontal scales are a strong hint that these are likely 320 instability features of breaking gravity waves like ripples. We generally find variable directions of propagation, which indicates that these wave-like structures may be mostly created above the stratospheric wind fields. However, a weak seasonal dependency is found: zonal directions of propagation are slightly more eastward during winter and westward during summer.
We speculate these to be instability features generated by breaking secondary gravity waves, receiving their zonal direction through advection by the background wind. We find a stronger tendency of southward propagation during summer, which may 325 point to a vital role of gravity wave filtering and excitation of secondary waves and their subsequent instability features by the meridional mesospheric circulation.
Furthermore, we observed and presented OH* imager observations of turbulence with high spatio-temporal resolution. We estimated turbulence parameters from 25 episodes of eddy observations. Following the approach of Hecht et al. (2021) we derived the energy dissipation rates for our observed events by reading the turbulent length and velocity scale from the image 330 series. Our values range between 0.08 and 9.03 W kg −1 and are higher than earlier rocket measurements. The values presented here would cause localized heatings of 0.03 -3.02 K per turbulence event. The largest of these reach the same order of magnitude as the daily chemical heating rates as reported by Marsh (2011). Given that the observed events are representative of typical processes of gravity wave dissipation, this emphasizes the importance of carefully integrating gravity wave turbulence into climate simulations. 335 Being able to derive reasonable values of UMLT turbulence parameters from imager data represents an important progress for measurement techniques of atmospheric dynamics. Airglow imagers are much cheaper and more flexible than rockets or lidars.
Considering the huge amount of data, artificial intelligence could be used in the future to identify and analyse turbulent episodes.

Data availability 340
The data are archived at the WDC-RSAT (World Data Center for Remote Sensing of the Atmosphere). The FAIM and GRIPS instruments are part of the Network for the Detection of Mesospheric Change, NDMC (https://ndmc.dlr.de).     Table 1). b) Histogram of .   Tables   Table 1. Episodes of turbulence observed at OTL and derived parameters from the image sequences. The duration of the turbulence events could not be determined if the vortex was not visible during its entire life span due to being partly outside the FOV ('out of FOV') of FAIM 3 or covered by clouds ('clouds'). In these cases, we noted the dissipated energy per mass and the maximum 565 temperature change as 'not available' (NA