Vertical Distribution of Optical Turbulence at the Peak Terskol Observatory and Mount Kurapdag

,


Introduction
Adaptive optics systems correcting optical distortions in a wide field of view require information about vertical profiles of optical turbulence and wind speed.The key characteristics of the optical turbulence are the structural constant of the air refractive index turbulent fluctuations C 2 n and the atmospheric coherence length, often associated with the Fried parameter r 0 and seeing (Coulman, 1985;Tokovinin, 2023).These quantities determine the possibilities of reaching a high angular resolution in astronomical observations.
There are a number of techniques for measurements of optical turbulence vertical profiles.In astronomical observations, the vertical profiles of optical turbulence above telescope can be determined with a multidirectional Shack-Hartmann wavefront sensor (Ran, 2024).In particular, the vertical profiles of C 2 n are calculated by applying the Solar Differential Image Motion Monitor (S-DIMM+) or Slope Detection and Ranging (SLODAR) techniques (Hickson , 2019(Hickson , , 2020;;Wang Z., 2018;Subramanian, 2023;Shikhovtsev, 2022;Joo, 2024;Bennoui, 2023;Deng, 2023;Wilson R., 2002).The most complete review of methods for measuring optical turbulence vertical profiles is given in a recent paper (Griffiths, 2024).Griffiths R. et al. present the results of a intercomparison of turbulence profiling instruments including Stereo Scintillation Detection and Ranging (S-SCIDAR), DIMM, 24-hour Shack-Hartmann Image Motion Monitor (24hSHIMM), Ring-Image Next Generation Scintillation Sensor (RINGSS), Multi-Aperture Scintillation Sensor-Differential Image Motion Monitor (MASS-DIMM).It is important to emphasize that, the comparison of these turbulence profiling instruments showed that the largest deviations correspond to estimates of the isoplanatic angle θ 0 , the values of which are highly sensitive to turbulence in the upper atmospheric layers.These deviations depicted in the form of scattered points are due to a limited accuracy in determination of the optical turbulence strength at different atmospheric levels and, apparently, differences in radiation propagation paths in measurements.
Measurement of optical turbulence vertical profiles is a hard procedure (Zhang, 2022).Most often, the profile measurements are carried out only for a separate set of nights or days.In order to estimate the optical turbulence characteristics over long time intervals or, in other words, for a statistical ensemble, special atmospheric models are widely used (Bolbasova, 2019;Qing, 2020;Bi , 2023;Yang, 1 2021;Quatresooz, 2023;Wu, 2024).The atmospheric models are applied to estimate the characteristics of mesoscale, micrometeorological and optical turbulence, on spatial scales comparable to the aperture diameters of ground-based telescopes.Below, a few recent papers demonstrating interesting results in this direction are discussed.A very interesting approach to predict optical turbulence was demonstrated (Cuevas, 2024).Cuevas et al applied a combination of models to improve the forecasts of C 2 n profiles and seeing at Paranal using Weather Research and Forecasting model.Also, Macatangay et al. (Macatangay , 2024) utilized advanced numerical simulations with the Weather Research and Forecasting model to predict astronomical seeing.Meteorological model used by the authors reproduced variations very well, especially in comparison with smoothed DIMM data.Using the 21-yr European Centre for Medium-Range Weather Forecasts' fifth set of reanalysis (Era-5) data, the optical turbulence and the wind speed characteristics were simulated at the Lenghu site (Bi , 2024).Analysis of the optical turbulence statistics shows that calculated median value of seeing is slightly lower than measured seeing.In particular, Era-5 derived value is 0.72 ′′ and measured seeing is 0.75 ′′ .Also, in our opinion, significant fundamental results in the optical turbulence diagnostics and forecasting with atmospheric models were obtained by English scientific teams: new methods and tools were created for studying optical turbulence (Masciadri , 2001;Turchi , 2017).
For estimating optical turbulence characteristics with atmospheric models, a key step is the selection of a small-scale turbulence parameterization scheme (Shikhovtsev , 2024).Correct selection and adjustment of the parameterization scheme to the atmospheric conditions of a given site are required to obtain representative estimates of the key characteristics of optical turbulence.Within mountainous regions, including Tibet and the Caucasus, turbulence behavior is complex.The choice of parameterization scheme for mountain regions should be based on analysis of preliminary measurement data.We should emphasize also that modern science still needs to study dependencies in turbulence evolution and determine to what extent these dependencies corresponding to a given mountainous region can be applied for another region.
The present study focuses on the simulation of optical turbulence above the Peak Terskol Observatory (43.27472 °N 42.50083 °E, 3127 m a.s.l.) and Mt.Kurapdag (41.79609 °N,47.37428 °E,3553 m a.s.l.).These sites were chosen as reference locations, for which data of special atmospheric measurements are available (Bolbasova, 2023).The sites of the Terskol Peak Observatory with high astroclimatic characteristics, Mount Kurapdag with low water vapor content (Khaikin, 2022) as well as the Big Telescope Alt-Azimuthal Telescope (BTA), the largest optical telescope with a non-segmented mirror in Eurasia are shown in figure 1. (Shikhovtsev, 2024).Recent studies show that new suitable locations for ground based telescopes include mountain regions of Dagestan.In this regard, the figure shows the location of the mountain village of Gunib as a place with potentially high characteristics.
In March 2023, the special experiments at the Peak Terskol Observatory have been carried to characterize the optical turbulence.Together with the measurements of optical turbulence profiles and meteorological characteristics in the surface layer, variations of C 2 n (z) have been simulated using a gradient method.Below, we discuss the results of estimating optical turbulence profiles using vertical shears of wind speed, mast meteorological measurements and measured height distributions of optical turbulence.

Data
This study aims to estimate the image quality in nighttime at the Terskol Peak Observatory, where optical turbulence monitoring is carried out, and at a new location, on the slope of Mt.Kurapdag, considered as one of the reference sites for methodological tests.In this study, we solved a number of tasks.Firstly, using remote telescopic measurements we determined vertical distributions of the optical turbulence intensity over the Peak Terskol observatory.Vertical distributions of turbulence intensity were obtained under conditions of good and mediocre image quality.Secondly, in order to estimate the vertical profiles of the structural constant of air refractive index turbulent fluctuations C 2 n , we used a gradient method, the foundations of which have been created by Tatarsky.Using the gradient method, the vertical profiles of optical turbulence over the Terskol Peak Observatory and Mt.Kurapdag were determined.
For simulation of C n (z) vertical profiles, we used the European Centre for Medium-Range Weather Forecasts' fifth set of reanalysis (Era-5) data (Hersbach, 2020).The Era-5 reanalysis contains observational data from around the world from 1940 to the present.Compared to previous databases, Era-5  (Huang, 2021;Rao , 2024).The re-analysis contains assimilated meteorological characteristics including air temperature and wind speed components at various pressure levels.The meteorological characteristics are available with a flat horizontal resolution of about 30 km and a continuous 1-hour temporal resolution.In particular, for the grid node closest to the observatory, hourly values of air temperature and horizontal components of wind speed at various pressure surfaces were used.

Optical turbulence vertical distributions from Shack-Hartmann sensor measurements at the Peak Terskol Observatory
In solar astronomy, Shack-Hartmann wavefront sensors with different fields of views are typically employed for estimated wavefront distortions (Townson , 2015).The parameters of standard sensor depend on the optical turbulence characteristics.In the paper, we present the results of an observing campaign of the Shack-Hartmann sensor at the Peak Terskol Observatory in 2023 March.We used the Shack-Hartmann sensor placed in the focal plane of MEADE LX200 Telescope installed inside the Schmidt-Cassegrain Telescope (figure 2).The sensor is equipped with a high frame-rate FLIR Grasshopper GS3-U3-28S5M-C camera.Subaperture size of the Shack-Hartmann sensor is 6x6 mm.Exposure time is 5 ms.The sampling is 65 frame/sec.In experiments, we analyzed 20 full time series and 10000 hartmannograms were obtained during the night of 19 March, 2023.Analysis of these measurement data made it possible to estimate and analyze the character of optical turbulence above the telescope.Obtained data sets correspond to average seeing conditions suitable for astronomical observations and thus were selected for processing.Assuming that the turbulence structure obeys the Kolmogorov model, we can estimate the height distribution of optical turbulence strength.For these calculations, we used an expansion of the spatial spectrum of amplitude perturbations into modes.Model spectrum of the scintillation is related to the height distribution of optical turbulence strength: where variance of the differential displacements of the solar subimages σ 2 α is related to the Fried parameter r 0 by the following equation: where K i is the numerical constant, λ is the light wavelength and D is the telescope diameter.
Analyzing these figures, we can see that pronounced turbulent layers are formed at heights of 7 -8 km and 2 -2.5 km above the ground.The turbulent layer at these heights leads to a significant decrease in image quality.For example, the parameter β determined by the integral C 2 n over height increases from 0.63 to 1.68 ′′ .These changes are due to, in mainly, generation of intensive turbulent fluctuations of air refractive index within these layers.
It is necessary to emphasize that, conditionally, we consider three different cases, including good seeing (β ranges from 0.5 to 1.0 ′′ ), average seeing (β ∼ 1.2 ′′ ) and bad seeing (β ranges from 1.5 to 2.0 ′′ ).Good seeing is associated with a thin lower turbulent layer (∼ 50 -100 m) and the absence of a layer at heights of 2 -2.5 km.With a thickness of the night boundary layer of 150-250 m, the seeing increases to ∼ 1.2 ′′ .Values higher than 1.5 ′′ are often associated with the development of optical turbulence within the surface layer and an additional turbulent layer at heights of 2 -2.5 km.

Analysis of sonic measurements within the atmospheric surface layer.
Estimation of optical turbulence strength C 2 n The lower part of the atmospheric boundary layer is known makes a major contribution to the formation of integral optical turbulence (along line of sight) and, as a consequence, to the resolving power of ground-based telescope (Kornilov , 2014;Panchuk , 2011).Within this layer, variations in the intensity of optical turbulence are most poorly described, both in time and in height.At least, the accuracy of optical turbulence estimating and forecasting based on standard atmospheric models or a small amount of optical remote measurement data is not sufficient for planning an observing time.In order to obtain correct shape of optical turbulence profiles within the lowest atmospheric layers under the different atmospheric conditions we analyzed measurement data carried out with a sonic anemometer.The sonic anemometer "Meteo-2" has been placed on a 7-m meteorological mast at the Peak Terskol Observatory site (figure 5).The observatory is located in the vicinity of the Elbrus, on a flat platform at a height of 3100 m above sea level.
Below, for calculations of atmospheric turbulence characteristics, we used long-period quasi-continuous measurement data of the sonic anemometer (from January 2023 to December 2023).Thanks to the high frequency measurements (10 Hz) of the wind field and the sonic temperature, we calculated a number of the turbulence characteristics including surface values of C 2 n .Based on the 1-year measurements with sonic anemometer, the 3 minute average values of C 2 n were obtained (Odintsov, 2019).In figure 6, we present the number of nights sonic measurements for each month.The fewest number of nights analyzed occurs in February-April.The largest number corresponds to January as well as warm period of observations.
Distributions of measured values of C 2 n within the surface layer of the atmosphere for different seasons at the Peak Terskol Observatory site are depicted in figures 7 and 8.
Table 1 contains the calculated percentiles of C 2 n in the surface layer of the atmosphere.Calculated values of C 2 n obtained from the results of processing sonic measurement data are higher in comparison with estimations derived from optical measurement data under clear sky conditions.This fact is confirmed by previous studies (Lukin, 2015).Higher amplitudes of C 2 n are largely determined by the influence of dynamic processes on the formation of small-scale turbulence in the surface layer of the atmosphere.Due to the fact that astronomical observations are carried out under good astro-optical conditions, when the intensity of optical turbulence along the line of sight is minimal, calculated values in the upper part of table1 differ from statistics corresponding only to clear sky conditions.In this regard, we evaluated the statistics of surface turbulence for medium atmospheric conditions and low intensity of optical turbulence.
Analysis of estimations shows that the intensity of turbulent fluctuations of the air refractive index has a pronounced seasonal dependence.The highest intensity of optical turbulence in the surface layer of the atmosphere is observed in winter.During this period, the median values of C 2 n , estimated for night time (from 0:00 to 6:00 local time), are equal to 8.6 • 10 −15 m −2/3 and 2.6 • 10 −15 m −2/3 for average atmospheric condition and under clear sky, respectively.In spring, we can an increase in the non-dimensional number of atmospheric situations F c with low intensity optical turbulence.At this time, during periods of high image quality, the median of C 2 n is about 8.0 • 10 −16 m −2/3 .In comparison with the winter and spring, in summer and autumn, turbulence in the surface layer is significantly suppressed: the number of atmospheric situations with strong turbulence C 2 n ≥ 10 −14 m −2/3 does not exceed 7%.We should emphasize that the largest number of cases with the best atmospheric conditions (C 2 n ≥ 10 −16 m −2/3 ), occurs in summer and autumn.In this period, about 7 -10 % of the observing time corresponds to weak optical turbulence.

Vertical distribution of optical turbulence strength estimated from Era-5 reanalysis
For simulation of the vertical profiles of optical turbulence, the gradient method was used.Correction of the vertical profiles was performed using sonic measurements and optical observations of image scintillation.The calculations were based on the following equation: where α is the numerical constant which equal to 2.8.The parameter M (z) is the vertical gradient of air refractive index, L 0 is the outer scale of turbulence.The calculation of the air refractive index gradients M (z) at different heights in the atmosphere is performed using the vertical gradients of the natural logarithm of the air potential temperature ∂lnθ ∂z : where P is the atmospheric pressure, T is the air temperature, θ is the air potential air temperature.The potential air temperature is related to the air temperature by the following formula: The key parameter is the outer scale of atmospheric turbulence L 0 Rao (2023).In general, vertical changes of the outer scale L 0 can be determined through vertical gradients of the horizontal component of the wind velocity S: where u and v are the horizontal components of wind velocity.Under conditions of weak turbulence the coefficient a=1.64, b=42, c=0.506 and d=50.Use of standard coefficients a=1.64, b=42, c=0.506 and d=50 gives mediocre results.In particular, we calculated a vertical profile C 2 n (z) for 18 -19 March, 2023 (figure 9 a)).By integrating this profile using formula 9, we obtained the average estimation of β parameter: where H is the height of the optically active atmosphere, α s is the zenith angle, λ is the wavelength of light.In average, the model value of β m is higher than the measured β and and the estimated C 2 n in the upper atmospheric layers are higher than expected values.In order to correct the vertical profile of C 2 n (z) we calculated new parameterization coefficients a new and b new by finding the best fitting model linking measured and model near-surface values of C 2 n .Specifically, new coefficients were calculated by minimizing the standard deviation between the model values C 2 n and the measured, sonic, values of C 2 n (m): where One of the most important characteristics of optical turbulence is the height of atmospheric boundary layer (ABL).For determination of the ABL height we have used a well-known method based on the profile of the bulk Richardson number.The ABL height is defined as the height where the bulk Richardson number exceeds a critical threshold.From figure 9, we can see that the new coefficients used to parameterize the night variations of C 2 n made it possible to reduce the height of ABL.Obtained night-time heights of the atmospheric boundary layer ranges from 100 to 300 m.As a result, the estimated strength of optical turbulent fluctuations decreases more rapidly with height at night above Terskol Peak.We associate this decrease of turbulent fluctuations strength within the lower atmospheric layer above the Terskol Peak with the low heights of the atmospheric boundary layer at night.Considering a large region, the minimum height of the atmospheric boundary layer is observed above Terskol peak not only on individual nights, but also on average, over long time period (figure 10).The high quality of astronomical images at the Terskol Peak is related to the low height of the atmospheric boundary layer.On some nights, its height is limited by 50 -100 m.Mt Kurapdag and Big Telescope Alt-Azimuthal are located in worse atmospheric conditions, the mean night-time heights of the atmospheric boundary layer are ∼ 250 -300 m and 250 -350 m respectively (figure 10 a)).We believe that observed character of vertical changes is in good agreement with theoretical and experimental data on the evolution of the atmospheric boundary layer.
The data of optical and meteorological measurements as well as the fields of wind speed and air temperature extracted from the re-analysis Era-5 , served as the basis for calculation of the average profiles of optical turbulence above the observatory and Mt.Kurapdag.In figure 11 a) we demonstrate calculated vertical profiles of optical turbulence above the Peak Terskol Observatory averaged over January, 2023 and June, 2023.
Analyzing these profiles, it can be noted that the structure of turbulence in the upper atmospheric layers changes significantly in summer.This fact is indicated by deformations of vertical profiles of optical turbulence.On the one hand, these deformations manifest in a greater vertical uniformity of optical turbulence and suppression of turbulent fluctuations in the surface layer of the atmosphere.On the other hand, the turbulence strength in the upper atmospheric layers decreases in summer.
Also, using Era-5 reanalysis data and identified patterns in optical turbulence behaviour at the Peak Terskol Observatory site we calculated averaged vertical profiles of C 2 n (z) above the Mt.Kurapdag.For comparison, the figure 11 b) shows vertical profiles above Terskol Peak and Kurapdag Mountain over 2013 -2023.Analysis of these vertical profiles indicates the presence of an intense lower turbulence layer.It should be noted that within the surface layer, the intensity of turbulence, in terms of C 2 n , at the Mt.Kurapdag site exceeds by more than an order of magnitude the values of C 2 n observed at  4 Estimation of the Fried parameter at the sites with good astroclimatic conditions Figure 12 displays the monthly averaged Fried parameter estimated using Era-5 reanalysis data at the ten astronomical sites.European Southern Observatory and Tibetan Plateau and the Ethiopian Plateau are some of the most good sites.In contrast to the Fried parameter of the European Southern Observatory and Tibetan Plateau (TP) sites, the Ethiopian Plateau has a different pattern, possibly owing to different climate patterns.During the year, the Ethiopian Plateau atmospheric coherence length was high in winter and low in summer.In winter, atmospheric coherence length at the Ethiopian Plateau is better than that at the ESO and TP sites, indicating that the Ethiopian Plateau has high potential conditions that are comparable to some of the best astronomical observatories in the world.
Estimated Fried parameter at the Ali site located in western China and the Tibetan plateau, one of the best astronomical sites, is about 20 cm in the best months.These sites are presently under studies for astronomy Hickson (2020).The Era-5 derived values of the Fried parameter corresponding to spatial coherence length are consistent with measured values at the Maidanak Observatory (Tillayev , 2023).
For comparison, we also compare the variation of the Fried parameter at the site of Peak Terskol. Figure 13 shows variations of the Fried parameter r 0 at the Terskol Peak Observatory.As you can see from the figure, the site of Terskol Peak has a high image quality close to the best places in the world.In particular, in summer β is about 1.0 -1.2 arc sec, in terms of r 0 , these values correspond to 8.3 -10.0 cm.In addition to the average monthly values of the Fried parameter obtained using the gradient method, the figure shows the results of applying the Hufnagel-Vally (HV) climate model to describe vertical changes in C 2 n (Mahmood , 2022).The HV model is based on the processing of experimental data and, on average, describes the structure of optical turbulence well.Nevertheless, detailed information about the vertical distribution of optical turbulence cannot be reflected using this model.Analysis of the figure shows that the gradient model and the model give similar results in terms of r 0 .

Discussion
In this paper, we present the results of studies of the optical turbulence in the atmospheric boundary layer and free atmosphere, up to a height of 20 km.In order to determine the structural constant of turbulent fluctuations of the air refractive index, we refined the parameterization scheme for calculation of this quantity.In the region of interest, this turbulent parameter was calculated based on a comparison of reanalysis data, mast meteorological measurements and observations of image scintillation.In particular, the processing of mast measurement data made it possible not only to estimate C 2 n statistics in the atmospheric surface layer, but also to select optimal parameterization coefficients that determine the contribution of wind speed shear instability to the formation of turbulence.Using these coefficients, we estimated C 2 n values at different heights in the atmospheric boundary layer more accurately.This made it possible to refine the shapes of the vertical profiles of optical turbulence.The correctness of our approach based on the using constant parametrization coefficients (within only atmospheric boundary layer) is supported by the fact that the structure of turbulence in the atmospheric boundary layer should be determined by processes in the surface layer of the atmosphere, at least for long time periods.
Correctness of the calculations is also confirmed by the results of estimating the height of the atmospheric boundary layer.Figures 10 show distributions of the atmospheric boundary layer height (in meters) within the Peak Terskol Observatory region.Analysis of these distributions allows us to notice that a fairly thin boundary layer is formed above the Observatory.This suggests that turbulence develops predominantly in the layer near the ground.During the measurements, the thickness of the atmospheric boundary layer was minimal (compared to average atmospheric conditions (10 a)).
It can be noted that the gradient method outputs are in good agreement with the data of optical measurements of turbulence at different height levels in the atmosphere, at least for the Terskol Peak Observatory.Based on the selected profile calculation scheme for the Terskol Peak Observatory, we estimated the optical turbulence characteristics for other sites as well.Analysis of spatial distribution of the atmospheric boundary layer height, as well as calculated profiles of optical turbulence, can serve as the basis for choosing optimal site for a ground-based astronomical telescope.Taking into account the recent data on the spatial distribution of astroclimatic characteristics within this macroregion, the following areas suitable for optical telescopes are located: in mountain Dagestan, including Gunib, as well as Shalbuzdag and Peak Terskol Observatory.The location of Mt.Kurapdag is not good due to strong surface winds and strong surface optical turbulence, which is the main reason of high values of β (2-3 ′′ ).

Conclusions
In this work we demonstrated that, using the gradient method and additional measurement data, we can estimate the vertical profiles of optical turbulence more accurately.The increase of accuracy of optical turbulence profiles within the lower atmosphere is due to using both measured values of C 2 n and applying new parametrization coefficients between vertical shears of wind speed and model values of C 2 n .Moreover, the gradient method takes into account variations in the height of the atmospheric boundary layer estimated from re-analysis data.Estimations of the atmospheric boundary layer height and the vertical profiles of wind shears are important for obtaining smooth vertical profiles of C 2 n and L 0 .Above the atmospheric boundary layer, vertical profiles of optical turbulence are reconstructed only on the basis of vertical profiles of wind speed shears with parameterization coefficients a=1.64 and b=42, that describe C 2 n variations well for weak turbulence.Thus, we can conclude that the hypotheses of optical turbulence embedded in the gradient method are valid as, in average, used atmospheric models converge to the measurements: the measured and model values of β are practically identical.
For the first time, by linking simulation outputs to measurement data, we have obtained representative optical turbulence vertical profiles for the sites of Peak Terskol Observatory and the Mt.Kurapdag.We have shown that the Kurapdag mountain should not be considered as a potential site for a new optical telescope.

Figure 1 :
Figure 1: Sites of the Terskol Peak Observatory, Big Telescope Alt-Azimuthal Telescope (BTA), Gunib with high astroclimatic indicators and Mount Kurapdag with low water vapor content.Red rectangle shows the BTA region

Figure 3 :
Figure 3: Night-time vertical distributions of optical turbulence strength C 2 n ∆z at the Peak Terskol Observatory for the average image quality, 19 March, 2023.In the figure, we used a local time

Figure 4 :
Figure 4: Night-time vertical distributions of optical turbulence strength C 2 n ∆z at the Peak Terskol Observatory for (a)) the best (β = 0.63 ′′ ) and (b)) the worst atmospheric conditions (β = 1.68 ′′ ), 19 March, 2023.In the figure, we used a local time

Figure 5
Figure 5: a) The dome of Schmidt-Cassegrain Telescope and b) Meteo-2 sonic anemometer

Figure 7 :Figure 8 :
Figure 7: Distributions of measured values of C 2 n within the surface layer of the atmosphere at the Peak Terskol Observatory site under average atmospheric conditions, 2023

Figure 10 :
Figure 9: Calculated vertical profile of optical turbulence at the Peak Terskol Observatory, 18 -19 March 2023.Shading corresponds to the interval between the first and the third quartiles estimated for all nights in March 2023

Figure 11 :
Figure 11: Night-time vertical profiles of optical turbulence above the Peak Terskol Observatory and Mt.Kurapdag 15

Table 1 :
Percentiles of C 2 n in the surface layer of the atmosphere at the Peak Terskol Observatory, 2023 is the measured value of C 2 n in the surface layer, C 2 ni is the model (calculated) value of C 2 n (surf ).Summation is performed over the entire time implementation, the length of which is T L .Using this minimization condition, the weight coefficients a new and b new are estimated.New parametrization coefficients at the Terskol Peak Observatory differ significantly from standard values: a new =0.95,b new =10.Figure 9 b) shows corrected vertical profile of C 2 n with the new parametrization coefficients.This profile corresponds to 19 March, 2023.