Atmospheric Conditions within Big Telescope Alt-Azimuthal Region and Possibilities of Astronomical Observations

: The paper presents the results of analysis of astroclimatic conditions in the Big Telescope Alt-azimuthal (BTA) region (40 ◦ N–50 ◦ N; 35 ◦ E–55 ◦ E). Using data from the European center for medium-range weather forecast ReAnalysis (ERA-5), we estimated the averaged spatial distributions in total cloud cover, vertical integral of mean kinetic energy, vertical component of wind speed, and wind speed shears, as well as inverse values of Richardson number 1/ Ri . An extensive region with the development of atmospheric ﬂows is formed south and southeast of BTA in winter. High inverse values of the Richardson number, spatial heterogeneities in vertical wind speed, and signiﬁcant wind speed shears in the lower atmosphere are observed in this region. In terms of turbulence development over BTA, the best time for astronomical observations falls in summer, when vertical shears of wind speed are weakened in the lower atmospheric layers. The situation is opposite in the upper troposphere. In winter, BTA is in the region of moderate vertical wind shears. In summer, a region with increased vertical wind speed shears is formed. Taking into account that the intensity of optical turbulence decreases rapidly with height, better image quality can be expected in summer. Such structure of the atmosphere does not allow one to directly apply atmospheric models in order to describe turbulence based on the turbulence strength as function of its ground values, or to use the classical model describing the turbulence velocity as function of air ﬂow velocity at the height corresponding to the 200 hPa level.


Introduction
Scheduling of a ground-based astronomical telescope depends strongly on atmospheric conditions. One of the main atmospheric characteristics is the amount of clear sky, which determines the scientific output capacity of the telescope. The amount of clear sky at the telescope site can be estimated either through the sunshine duration or total cloud cover (TCC). The analysis of clear sky and TCC is of special importance not only when choosing sites for new telescopes, but also at the existing astronomical sites [1][2][3].
As a result, on the telescope focal plane, one can observe image motion, scintillation, and spread of fine structure in the images of solar objects or stars. Image motion and spread of fine structure are determined by the sum of intensities of turbulent fluctuations along the telescope line of sight, regardless of the height of turbulent layers. Scintillation and the size of the field of view with the given image quality criteria depend on the height of turbulent layers. As the intensity of turbulence in the upper atmosphere grows, scintillation and so-called isoplanatic angle of the atmosphere increase significantly [4,5]. The isoplanatic angle determines the size of the field of view, within which optical distortions are minimal.
In the article, we consider astroclimatic characteristics within the Big Telescope Alt-Azimuthal region.
The Big Telescope Alt-Azimuthal is one of the key astronomical telescopes in Russia. The Russian 6 m telescope has been successfully operating for almost 48 years [6]. BTA remains among the top twenty largest astronomical instruments. The Special Astrophysical Observatory provides continual operation of the telescope upon requests of both Russian and foreign researchers. Scientific equipment required for astrophysical observations is made and designed by staff members of the Observatory. The telescope is located near Mt. Pastukhova at the altitude of 2070 m above sea level. Scientific equipment, observational methods, and the status of BTA are described in [6,7].
The region of BTA's location is promising in terms of placing new ground-based telescopes operating in the optical range of the electromagnetic spectrum and in the millimeter/submillimeter range. The modern concept of development of ground-based infrastructure in the millimeter/submillimeter range involves construction of three new high-performance millimeter/submillimeter telescopes [8]. One of those is supposed to be built in Russia. This study continues the search for sites suitable for placing a millimeter/submillimeter telescope in Russia and estimation of astroclimatic characteristics of possible sites [9,10].
For the BTA site, we compared variations in observation time with low-frequency variations in TCC. In terms of selecting sites suitable for astronomical telescopes in the BTA region, we analyzed spatial distributions of TCC. To identify the turbulence features within the BTA region, we used the gradient Richardson number Ri associated with the development or suppression of the atmospheric turbulence. We analyzed 1/Ri vertical profiles in the lower turbulent layers of the atmosphere and at the heights corresponding to the 200 hPa level. In a number of studies, the wind speed V 200 on the 200 hPa isobaric surface is considered as a turbulence parameter [11][12][13][14][15][16].
Comparison of ERA-5 data with data measured at different sites shows their good agreement, at least, for the free atmosphere. For example, Virman et al. [18] compared the air temperature vertical profiles derived from radiosounding data with those derived from ERA-5. The authors showed that the largest deviations are observed in the lower layers of the atmosphere for all analyzed stations. The reanalysis data are successfully used in the studies of the structure of atmospheric flows and in order to estimate parameters of large-scale turbulence [11,[19][20][21][22][23][24].
To determine profiles of the 1/Ri parameter and identify atmospheric layers, we also used radiosounding data obtained at the nearest station (Mineralnye Vody (http://weather. uwyo.edu/upperair/sounding.html, accessed on 10 January 2022)). Radiosounding data include atmospheric pressure, height, air temperature, dew point temperature, relative air humidity, wind speed and direction, and potential air temperature. These data are available at different heights in the atmosphere at 00 and 12 h. Figure 1 shows the height changes in vertical displacements of radiosondes.
Vertical displacements were estimated as the difference between the heights of the radiosonde δz = z i − z j (z i > z j ) at different points in time. In winter, the most rapid increase in displacements is observed up to the 1060 m height. In summer, this height reaches 1560 m. Above, vertical displacements continue to increase at a slower rate. Starting from the height of 6000 m, we can assume that the radiosonde moves with approximately constant vertical speed. Variations in the radiosonde speed indicate instabilities in different layers of the atmosphere. It can be noted that high probability of repeatability of minor radiosonde displacements is observed in the lower atmospheric layers, as well as in the layer at 11,200-13,500 m. On average, vertical resolution of the radiosonde data used is 425 m. In the lower atmospheric layer (up to 3000 m), vertical resolution is 235 m.

Total Cloud Cover Variations over BTA
In order to estimate the total cloud cover within the BTA region we used ERA-5 monthly averaged data on single levels from 1979 to 2021. Based on ERA-5 data, we estimated TCC at the BTA site. Figure 2 shows changes in the monthly averaged TCC at the BTA site. From the figure, one can see that TCC decreases from 1979 to 2021. Table 1 presents statistical characteristics of cloudiness for different time intervals. Comparing two intervals (1979-2002; 2012-2021), we can see that at the BTA site, the mean and median values decreased by 0.03 and 0.04, respectively.   Table 2 presents the monthly averaged TCC. Analysis of Table 2 shows that at the BTA site, the monthly averaged TCC decrease, except the period from January to March. The maximum decrease in TCC occurs in November and equals 0.07. At the site of the reference station (Teberda), TCC decreases for all months, except for March and September. At Teberda, the maximum decrease in TCC also falls at November and amounts 0.14.
Comparison of the monthly averaged TCC derived from ERA-5 with the weather station data for available time intervals suggests that ERA-5 overestimates TCC by at least 0.05.  In addition, using ERA-5 data, we compared the observation time at BTA with variations in TCC. Figure 4 shows the amount of observation time at the telescope according to data of the operation service and astronomical observations, as well as low-frequency changes in TCC at BTA. The amount of observation time was calculated from the operation service data, estimated by Panchuk [25]. The observation time begins when the telescope shelter opens and ends when the shelter is closed. The amount of observation time from astronomical observations is the estimate of the clear sky duration from below. It includes only dark-sky optical observation time, excluding telescope downtime and faults.
We assumed that the changes in the amount of observation time indicate only longperiod variations in clear sky. Due to the fact that TCC changes significantly within hours, days and months, we divided the initial time series into low-frequency and high-frequency components. To subtract the high frequency component in the TCC changes, we used the Butterworth bandpass filter with monotonic frequency response. The Butterworth filter frequency response H( f ) is determined by the function: where f is the frequency, and f c is the cutoff frequency. The order of the Butterworth window is given by n r , and H 0 is the amplification coefficient. In the analysis of TCC time series, we used n r =6 and cutoff frequency f c = 2π/T c = 100 months. Figure 5 shows the high-frequency component of TCC changes, which we did not take into account when calculating low-frequency changes in TCC at the BTA site. Figure 6 shows the energy spectrum of the full TCC time series (orange line) and the spectrum of the high-frequency component of TCC changes. Analysis of Figure 6 shows that we did not take into account the total cloudiness fluctuations with periods of less than 24 months for variations in Figure 4.

Spatial Distributions in TCC within the BTA Region
We estimated TCC spatial distributions to find new sites suitable for large groundbased astronomical telescopes. Figure 7 shows TCC spatial distributions averaged in winter and summer for the 2012-2021 period within the BTA region. These distributions are obtained for the region within 40 • N-50 • N and 35 • E-55 • E. TCC values are normalized to the maximum TCC for the region. The analysis of TCC spatial distributions shows that cloudiness is unevenly distributed. There are local regions with increased and decreased TCC.
In winter, BTA is in the transition region, between the local TCC maximum centered southwest of the telescope and the region of low TCC extended in the southeast direction (this region is shown by blue color). The pronounced TCC maximum (corresponds to yellow color) is located north-northeast of BTA. In summer, TCC decreases significantly for this region. However, BTA belongs to the region with high TCC. This pattern of TCC spatial distribution explains low amplitude of seasonal changes in cloudiness at the BTA site [10]. Thus, we can conclude that for 1993-2015, observation time increased by 76 h. The averaged TCC decreased by 0.03.

Structure Constant of Turbulent Fluctuations in the Air Refractive Index at BTA
We used data from a sonic anemometer on a 20 m meteorological mast to estimate the structure constant of turbulent fluctuations in the air refractive index C 2 n . The mast was installed about 24 m southwest of BTA. The sonic anemometer comprises four pairs of ultrasonic transducers capable of measuring air temperature and wind speed in three orthogonal axes. C 2 n is the parameter included in the Kolmogorov-Obukhov "2/3"-power law: where n is the air refraction index, D n is the structure function, and r is the distance between two given points. As fluctuations of the air refraction index are related to fluctuations of air temperature, C 2 n is proportional to the structure constant of air temperature fluctuations C 2 T : where A is the constant, P is the atmospheric pressure, T is the air temperature, <T> is the mean air temperature, V is the mean wind speed, and δt is the time shift. Figure 8 shows The mean value of C 2 n is 1.1 × 10 −15 cm −2/3 for the period from 25 October 2012 to 1 November 2012. In summer, the mean C 2 n decreases and equals 3.2 × 10 −16 cm −2/3 . In addition, we should note that in summer, C 2 n changes within narrower limits compared to autumn. Further, the estimates of C 2 n may be used in the model of turbulent atmosphere for BTA proposed in [26].

Era-5 Data and Richardson Number
As already mentioned, turbulence is the key atmospheric phenomenon limiting the duration of astronomical observations with given image quality criteria. To calculate statistical characteristics of image quality, including the ground-based telescope resolving power in the atmosphere, and predict the contrast and sharpness of images, it is necessary to measure turbulence characteristics at different heights in the atmosphere. These measurements are often carried out within individual programs and for different sites. Due to the lack of data, to describe the structure of atmospheric turbulence within the BTA region, we used the classical approach based on calculation of the Richardson number: where g is the acceleration of gravity, z is the height, and V is the wind speed horizontal component, which is calculated by the formula: where u and v are the wind speed zonal and meridional components, and θ is the potential air temperature, which is reduced to the standard pressure according to the dry adiabatic law: where P 0 is the standard atmospheric pressure, and P is the atmospheric pressure at a certain height above the Earth's surface. Negative values of Ri are observed during unstable thermal stratification, when the air temperature vertical gradient is higher than the dry adiabatic temperature gradient. If the air temperature vertical gradient is lower than the dry adiabatic temperature gradient, the stratification of the atmosphere is stable. In this case, turbulence is suppressed. When Ri = 0, neutral stratification is observed. In neutral stratification, the air density changes with height according to the adiabatic density gradient. The Richardson number does not directly determine the intensity of atmospheric turbulence. However, the analysis of this number indicates the possibility of turbulence formation. The Richardson parameter critical value serves as a criterion for turbulence development. The essential condition for turbulence is Ri < 0.25 [27].
To determine conditions for turbulence formation in the vicinity of BTA, we used ERA-5 data for a long period from 2012 to 2021. We estimated spatial distributions of the 1/Ri i parameter. Due to the fact that turbulent layers formed at different heights in the atmosphere have different effects on image quality, we calculated 1/Ri spatial distributions for the following layers: We estimated Ri using the following formula: where k is the number of atmospheric layers taken into account in calculations of Ri.  In winter, the BTA site is between the region with negative 1/Ri (southeast) associated with the turbulence development and 1/Ri local maximum west of BTA. The critical Richardson number is directly proportional to the degree of spatial inhomogeneity of the turbulence field. In other words, as the spatial inhomogeneity of the turbulence field increases, the critical Richardson number increases. As a result, turbulence can also persist under greater thermal stability of the atmosphere.
In winter, conditions with the transformation of kinetic energy and instability energy into the kinetic energy of turbulent motions are observed in the region south and southeast of BTA. In summer, the 1/Ri field changes significantly. High 1/Ri values are observed in a smaller region south of BTA. In summer, the maximum 1/Ri in this region increases 2.2 times. On average, a slightly stable stratification of the atmosphere is typical for the BTA site. Figure 10 shows  In winter, calculations based on ERA-5 data in the lower atmosphere give high values of wind speed shears in the vicinity of BTA (Figure 11a). In summer, the wind speed shears over BTA decrease significantly (Figure 11b). In the upper troposphere, the situation is opposite. In winter, BTA belongs to the region with moderate wind speed shears (Figure 12a). In summer, a region with increased wind speed shears is formed over BTA (Figure 12b).
The obtained 1/Ri estimates agree with data from other astronomical observatories. For example, to determine 1/Ri profiles, Y. Hach et al. [19] used data from the National Centers for Environmental Prediction (NCEP)/National Centers for Atmospheric Research (NCAR) Reanalysis. The obtained 1/Ri profiles for Oukaimeden and La Palma revealed 1/Ri ranging from ∼0.01 to 0.16 in the lower 20 km atmospheric layer and close to zero values in higher layers.  Mineralniye Vody is the only station located withinthe distance of ∼100 km from BTA. It can be considered as a reference point, at least to describe processes in the middle and upper atmosphere. We believe that the measurements obtained at the Mineralnye Vody station reflect the average atmospheric conditions at the BTA site (at least above 2000 m).
Using the radiosounding data for January 2021 and July 2021, we determined 1/Ri vertical profiles at Mineralnye Vody, shown in Figures 13 and 14. The analysis of vertical profiles at Mineralnye Vody makes it possible to identify atmospheric layers with high 1/Ri. In winter, the atmosphere is predominantly characterized by positive Ri. On average, atmospheric layers, for which the probability of turbulence generation is high, are observed at the heights of ∼250 m and 400 m. For 00 h, the layer is located at 100-280 m (with 1/Ri = 3.3). At 12 h, the layer with 1/Ri = 3.2 corresponds to the heights of ∼500-640 m. In the daytime, negative 1/Ri values are fixed at the height of ∼640 m. Thus, high amplitudes of 1/Ri ∼ 4 indicate the possibility of turbulence generation within the lower atmospheric layer. In January, thickness of this layer varies from 930 m at night to 1170 m in the daytime. In the middle troposphere, the 1/Ri maximum is revealed at 4835 m at night; in the daytime, the height of its formation shifts by 600 m and becomes 5435 m. In addition, small peaks are observed at 10,300-11,400 m below the large-scale jet stream. At night, an additional layer with increased Ri is observed at 8500 m. In the daytime, a turbulent layer is formed at ∼3300 m.
In July, 1/Ri vertical profiles are significantly deformed (Figure 14). In contrast to winter, when negative 1/Ri is predominantly observed in certain time intervals, 1/Ri values are significantly below zero in the lower atmospheric layers in July. In July, the thickness of the most active layer varies from 1360 m at night to 1560 m in the daytime. The maximum 1/Ri = 2.9 in the daytime corresponds to the height of 4860 m; at night the maximum 1/Ri = 1.4 falls at the height of 1360 m. At the BTA site, 1/Ri in summer varies from 0.1 to 0.53. The Richardson parameter Ri varies in the range from 1.9 to 10. Slightly increased 1/Ri amplitudes and the nature of changes over time indicate that in this area, turbulence is formed irregularly, and increased intensity of fluctuations can be observed only in certain time intervals. Such changes allow us to speak about high repeatability of the main parameter of optical turbulence, namely, high values of the Fried parameter r 0 estimated for higher atmospheric layers. Recall that the Fried parameter is determined by the optical turbulence profile: where K = 2π λ , λ is the wavelength of light, α is the zenith angle, H 0 is the height of the lower layer boundary, H A is the height of the upper layer boundary, and C 2 n (z) is the profile of the structure constant of turbulent fluctuations in the air refractive index, which characterizes the intensity of optical turbulence.
Knowledge of C 2 n vertical profile variations is the basis for describing changes in the statistical characteristics of image quality in turbulent atmosphere. For example, the atmospheric resolving power of telescope β is determined by the formula: At low pressure levels, superadiabatic air temperature gradients can be observed only in very thin layers (thickness of several tens or hundreds of meters). Due to low vertical resolution, ERA-5 reanalysis is not sensitive to such air temperature gradients. Low vertical gradients of air temperature lead to a narrow range of 1/Ri changes for BTA.
In the 825-650 hPa layer, the integral 1/Ri estimated from radiosounding data is 0.64 (in January). In July, 1/Ri decreases due to higher repeatability of negative Ri and amounts to 0.25. In the upper layer of the atmosphere (225-150 hPa), 1/Ri varies slightly during the year, from 0.24 in summertime to 0.27 in winter.
It is generally accepted that turbulence is formed when the Richardson number is less than 0.25. This threshold is the result of the linear instability theory. However, turbulence is also observed at Ri > 0. 25. In comparison with the lower atmospheric layers, at heights above 200 hPa, the repeatability of atmospheric situations with Ri < 0.25 is negligible. Using individual 1/Ri(z) profiles, we estimated the repeatability of atmospheric situations with Ri < 0.25 for heights above 200 hPa. The repeatability does not exceed 10-20%. Figure A1 shows that negative Ri values are rarely observed in the upper atmospheric layers.
Taking the account of large vertical wind speed shears leading to the increased turbulent mixing, we can expect that turbulence can be generated and evolve even with the observed range of 1/Ri changes. In this case, typical dynamic Richardson number for these heights should be less than 1 [28]. The estimated 1/Ri changes with height are consistent with data of studies [29], which demonstrate high Ri variability in the lower layers of the atmosphere.

Richardson Number and Vertical Integral of Kinetic Energy
Regions with low and high 1/Ri are associated with the processes of transformation of kinetic energy of the mean motion and instability energy into kinetic energy of turbulent motions. We have obtained spatial distributions of vertical integral of kinetic energy for 2012-2021 within the BTA region. These distributions in winter and summer are shown in Figure 15.
Analysis of Figure 15 shows that in winter, a region with low vertical integral of kinetic energy is formed east of BTA. The minimum vertical integral of kinetic energy is 162.9 × 10 4 J/m 2 . The position of the region with low vertical integral of kinetic energy corresponds well to the atmospheric region with high spatial inhomogeneity in the 1/Ri parameter estimated for the lower layers of the atmosphere. Thus, generation of turbulence at BTA in winter is determined mainly by the lower layers of the atmosphere. During summer, BTA is in a "hollow", determined by low values of the vertical integral of kinetic energy. Compared to winter values, vertical integral of kinetic energy in summer at BTA decreases 1.5-1.7 times.

Vertical Component of Wind Speed
It is known that the wind speed vertical component is a parameter of stability of air masses. Development of vertical movements is the indicator of atmospheric turbulence. For example, improvements in the quality of images near downward air flows are confirmed by acoustic and optical measurements of turbulent fluctuations performed by Nosov V. et al. [9]. In addition, analysis of spatial variability of the wind speed vertical component enables us to judge about the effect of mesoscale disturbances on the structure of air flows on smaller scales. In a stably stratified atmospheric boundary layer, the role of large-scale components [30], mesoscale disturbances, and wave-turbulence interactions in turbulence evolution increases. Large-scale coherent turbulence components can cause mixing across the entire flow thickness [31] and lead to changes in the energy structure of small-scale turbulent motions. The potential significance of modulation effects of mesoscale fluctuations in turbulence was also indicated in [32]. To analyze the structure of flows within the BTA region, we obtained distributions of the wind speed vertical component in winter and summer. The distributions are shown in Figure 16. The analysis of distribution of the wind speed vertical component (Figure 16a) shows that in winter, BTA belongs to the region of alternating downward and upward movements. Such pattern of spatial distribution of the wind speed vertical component suggests that atmospheric flows are deformed significantly above the complex orography. Conditions of the atmosphere stable stratification (positive Ri), a certain air flow orientation, and mountain characteristics contribute to the occurrence and development of mountain waves. Atmospheric turbulence is known to be generated not only by convective instability and Kelvin-Helmholtz instability, but also by mountain-wave breaking [33]. In [34], the authors investigated the relationship between the inertia-gravity wave energy and turbulence. The authors showed that the wave instability parameter correlates with the gravity wave energy and turbulence frequency with correlation coefficients ranging from 0.52 to 0.7. We can expect that in winter, mountain waves are formed south and southeast of BTA. Mountain waves can be an additional source of energy for small-scale turbulence in the spectrum of atmospheric flows within the BTA region. In summer, the spatial heterogeneity in the wind speed vertical component decreases. Compared to winter, the number of regions with ascending currents (with the velocity greater than 0.4 Pa/s) reduces significantly. Apparently, contribution of mountain waves to the formation of turbulence over BTA in summer is insignificant.

Results
The study was carried out to describe astroclimatic conditions in the BTA region (from 35 • E to 55 • E, from 40 • N to 50 • N). Using ERA-5 data and network meteorological data, the following results were obtained: (i) Monthly averaged TCC values were estimated for the BTA site for 1979-2021. Using the Butterworth bandpass filter, changes in the TCC low-frequency component were shown to be in good agreement with the amount of observation time at the telescope. The best time for astronomical observations (in terms of TCC) falls at summer. In particular, for 2002-2021, TCC reaches its minimum in August. The mean TCC is 0.46. (ii) The analysis of TCC long-term changes showed that the mean and median values decreased by 0.03 and 0.04, respectively. The duration of astronomical observations has increased by 76 h. The root mean square deviation of TCC remained unchanged. Thus, we believe that the decrease in TCC over a long time interval has led to increased duration of astronomical observations. (iii) The analysis of TCC spatial distributions showed that BTA was located between the local TCC maximum with the center to the southwest and the low TCC region extending southeastward in winter. In summer, TCC decreases significantly for the BTA region. However, BTA is located in the region with high TCC (compared to background TCC). (iv) The distribution of the 1/ parameter and vertical integral of kinetic energy within the BTA region indicates that in summer, the quality of astronomical images associated with the lower layers of the atmosphere is higher compared to winter. The higher image quality is due to the greater stability of the free atmosphere (low 1/Ri), and to BTA location at the periphery of the region with low vertical integral of kinetic energy. The analysis of winter spatial distributions shows that the vast region of low vertical integral of kinetic energy can be associated with the intense process of transformation of the average kinetic energy into the kinetic energy of turbulence. In general, winter is characterized by the increased spatial non-uniformity of vertical movements, presence of pronounced layers in 1/Ri vertical profiles, and significant wind speed shears in the lower atmosphere. (v) In winter, ERA-5-derived 1/Ri in the lower layers of the atmosphere is close to 0.015. This value indicates that turbulence is formed under weak static stability conditions at the BTA site. Stability of the atmosphere increases with height. (vi) Using radiosounding data, we estimated the 1/Ri integral in the layer from 825 to 650 hPa. The integral is 0.64 in January. The 1/Ri decreases in July due to higher repeatability of negative Ri values and amounts to 0.25. In the upper layer of the atmosphere (225-150 hPa), the average 1/Ri varies slightly during the year and amounts to 0.25. Overestimated values of 1/Ri from radiosounding data in comparison with the reanalysis data are related to insufficient vertical resolution in ERA-5 reanalysis. (vii) One of the main conclusions is that potential sites with high astroclimatic parameters are located east and southeast of BTA. Analyzing the spatial distributions of TCC, 1/Ri, vertical integral of kinetic energy, and the wind speed vertical component, we can conclude that sites suitable for new large astronomical telescopes can be found on the periphery of the low TCC region, in its southeastern part (40.5 • N-42.0 • N; 46.2 • E-48.7 • E). (viii) At the BTA site, the most pronounced atmospheric layer with the wind speed exceeding 15 m/s is observed at the height of ∼12 km [10]. Taking into account the wind speed profiles V(z) [10] and 1/Ri(z), we can conclude that the conditions for turbulence generation are observed mainly in winter under the jet stream (at 10,300-11,400 m). Under the jet stream, high 1/Ri amplitudes and the pattern of changes indicate that turbulence is formed irregularly in this region. Increased strength of fluctuations can be observed only in certain time intervals. These changes allow us to speak about high repeatability of the main parameter of optical turbulence, namely, high values of the Fried parameter r 0 , estimated for higher atmospheric layers. Compared to the lower atmospheric layers, the repeatability of atmospheric situations with Ri < 0 at the heights above 200 hPa is insignificant. Using individual 1/Ri(z) profiles, we estimate the repeatability of atmospheric situations with Ri < 0 for heights above 200 hPa as no more than 10-20% . (ix) In winter, high values of wind speed shears are observed in the lower layers of the atmosphere at BTA. In summer, the wind speed vertical shears over BTA are significantly reduced. In the upper troposphere, the situation is the opposite. In winter, BTA is in the region with moderate wind speed shears. In summer, a region with high wind speed shears is formed. Taking into account that C 2 n decreases rapidly with height, better image quality (associated with higher values of the Fried parameter) can be expected in summer. Such energy structure of the atmosphere does not allow applying atmospheric models to describe turbulence based on C n (z) as a function of surface values C n (0), or to use a model based on turbulence velocity at the height of 200 hPa [35].
The results obtained can be used in the search for new sites suitable for placing new ground-based medium-and large-aperture optical and millimeter telescopes.