Next Article in Journal
Impact of Aerosol-Cloud Cycling on Aqueous Secondary Organic Aerosol Formation
Next Article in Special Issue
Statistical Analysis of the Spatiotemporal Distribution of Ozone Induced by Cut-Off Lows in the Upper Troposphere and Lower Stratosphere over Northeast Asia
Previous Article in Journal
Activation of sub 2 nm Water Soluble and Insoluble Standard Ions with Saturated Vapors of Butanol in a Boosted TSI ultrafine CPC
Previous Article in Special Issue
Atmospheric Characterization Based on Relative Humidity Control at Optical Turbulence Generator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistics of the Optical Turbulence from the Micrometeorological Measurements at the Baykal Astrophysical Observatory Site

1
Institute of Solar-Terrestrial Physics SB RAS, Irkutsk 664033, Russia
2
V.E. Zuev Institute of Atmospheric Optics SB RAS, Tomsk 634021, Russia
3
Institute on Laser and Information Technologies, Shatura 140700, Russia
4
Federal State Budgetary Institution of Higher Vocational Education Irkutsk State University, Irkutsk 664003, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Atmosphere 2019, 10(11), 661; https://doi.org/10.3390/atmos10110661
Submission received: 26 September 2019 / Revised: 21 October 2019 / Accepted: 27 October 2019 / Published: 30 October 2019
(This article belongs to the Special Issue Atmospheric Turbulence Measurements and Calibration)

Abstract

:
The paper focuses on the investigations of the optical turbulence structure. The distributions of the repeatability of the structure characteristic of the air refractive index are obtained. A scenario of decreasing the intensity of the optical turbulence in the summer is discussed. Numerical estimates of this scenario are given. Using classical methods of wavefront analysis, the results of the first studies of the vertical structure of atmospheric turbulence at the Baykal astrophysical observatory site are presented.

1. Introduction

The degradation of image quality due to the atmospheric turbulence has been a long-known problem. Turbulent fluctuations of both air temperature and wind speed induce variations in the air refractive index. The study of the atmospheric turbulence and its optical properties is important for astroclimatic research, planning time of astronomical observations as well as searching for the physical mechanisms of its formation. In particular, the study of the surface layer over rough terrain is relevant. In rough terrain, the structure of the turbulence is broken. With wave and vortex disturbances arising, mesoscale airflows may induce the intense turbulence and variations in the exchange of heat, momentum, and mass. The work is aimed at studying the structure of optical turbulence at the Baykal astrophysical observatory site. Statistics of the optical turbulence includes a lot of parameters, some of which are analyzed in [1,2]. We focus on the study of the the main characteristic of the optical turbulence that affects the image quality is the structure characteristic of the air refractive index fluctuations C n 2 . The quantity C n 2 is the parameter included in Kolmogorov–Obukhov “ 2 / 3 ”-power law:
D n = C n 2 r 2 / 3 ,
where n is the air refraction index, D n is the structure function, and r is the distance between two given points at which the air refraction index is determined.
As the air refraction index fluctuations are related mainly to the air temperature fluctuations, C n 2 is proportional to the structure characteristic of the air temperature fluctuations C T 2 :
C n 2 = A P < T > 2 2 C T 2 ,
where A is constant for the light wavelength λ ( A = 79.3 × 10 6 at λ = 0.5 μ m), P—the atmospheric pressure, and < T > —averaged air temperature for the given time period. The quantity C n 2 defines the well known Fried radius r 0 , which described the quality of seeing [3]:
r 0 = 0.423 k 2 s e c α 0 H C n 2 ( z ) d z 3 / 5 ,
where α —zenith angle, k = 2 π / λ , λ —the wavelength of the light, H 20 km—the height of the “optically active” atmosphere.
In addition, the isoplanatic angle is determined by the vertical profile C n 2 ( z ) [3]:
θ 0 = 2.91 2 π λ 2 sec 8 / 3 α 0 H C n 2 ( z ) z 5 / 3 d z 3 / 5 .
In order to improve our knowledge about the statistical properties of turbulent mixing in the atmospheric boundary layer, the exchange of turbulent kinetic and potential energy in the atmosphere as well as optical turbulence, long-period measurements of both the air temperature, and wind speed have been performed.

2. The Repeatability of the Optical Turbulence at the Baykal Astrophysical Observatory Site

To study the atmospheric turbulence, large scale air motions, temporal variations of the mean meteorological characteristics, the data of ultrasonic anemometers are widely used [4,5,6,7,8].The results of a field intercomparison experiment of a number of sonic anemometer models point to biases in the measured characteristics. However, biases are generally very small for these sensors [9]. Comparative statistical analysis of the mean meteorological characteristics measured by the ultrasonic anemometers with data of the lidars and mechanical anemometers shows a very good agreement [6,7,10,11].
To estimate the strength of the optical turbulence in the atmospheric boundary layer at the Baykal astrophysical observatory site, we performed the measurements of the air temperature and wind speed using meteorological weather station “AMK-03” (sonic anemometer) mounted on a mast (Figure 1) [8]. The height of the meteorological weather station above the underlying surface was about 30 m. Using the meteorological weather station “AMK-03” mounted on the mast, we carried out the measurements every day during the time period from 1 January to 31 December 2018. The estimations of the air temperature and components of the wind speed are based on the dependence of the group speed of sound in the atmosphere on the air temperature. The data allow us to obtain long period time series of the air temperature and wind speed measured with high sampling rate (10 Hz). Initial averaging of the air temperature, the structure characteristic of the air temperature fluctuations and the structure characteristic of the air refractive index fluctuations was performed every 3 min. The time interval between two measurements was 10 min. The total amount of accumulated data are about 17,000 values of averaged structure characteristic of the air refractive index fluctuations during the year. Using time series recorded, we assess mean air temperature < T > and C n 2 using a standard approach [7].
The distribution of the repeatability of the structure characteristic of the air refractive index fluctuations for the period from from 1 January to 31 December 2018 is shown in Figure 2.
It is important that there is a significant asymmetry in the repeatability of C n 2 . The atmospheric boundary layer at the Baykal astrophysical observatory site is characterized by a moderate turbulence. The values of C n 2 often exceed 10 14 m 2 / 3 .
To analyze the temporal variability of the optical turbulence, we estimated the distributions of the repeatability of the structure characteristic of the air refractive index fluctuations in different seasons. The repeatabilities of the optical turbulence at the Baykal astrophysical observatory site in summer, autumn, winter, and spring are shown in Figure 3, Figure 4, Figure 5 and Figure 6, respectively.
Unfortunately, the measurement data include not only calm weather, clear sky, but also strong winds, cases when the wind speed exceeded 15 m/s. Quantitative estimates of the repeatability of the optical turbulence are shown in Table 1. For convenience, repeatability is given in fractions R:
R = n i N ,
where n i is the number of cases in fixed range of C n 2 changes, and N is the total number of cases in season. The analysis of obtained distributions shows that the structure of the optical turbulence changes significantly during the year. The strongest optical turbulence is observed in the winter, the mean value of C n 2 = 3.3 × 10 13 m 2 / 3 . The most favorable conditions are noted in summer, the mean value of C n 2 = 3.1 × 10 15 m 2 / 3 . It is significant that cases with weak turbulence (∼ 10 17 m 2 / 3 ) are observed throughout the year except for winter. In winter, the range of variations of C n 2 is significantly narrowed. This may be due to the fact that the measurements covered the cold period with open water (the lake has not frozen yet). In addition, compared with spring and summer, the number of cases with moderate turbulence increases in the autumn.
Investigating the C n 2 distributions discussed above, the surface layer turbulence is found to be weaker in summer than in winter. First of all, this is due to the suppression of the surface layer turbulence under the influence of the cold air over the lake in the summer. The high stability of the surface layer in the summer causes a decreasing of the turbulent fluctuations of both wind speed and air temperature. Moreover, significantly increasing the intensity of turbulence and a narrow range of C n 2 changes in winter can be associated with the action of open (unfrozen) water. A warm water surface provides high vertical air temperature gradients inducing the development of convective instability.

3. The First Study of the Vertical Structure of the Surface Layer at the Baykal Astrophysical Observatory

In the previous chapter, we discussed the distributions of the repeatabilities of the structure characteristic of the air refractive index fluctuations that characterize the strength of the optical turbulence at height of 30 m. Although the optical turbulence is the most intense in the lower atmospheric layers, the upper layers can make a significant contribution to variations of Fried radius and the isoplanatic angle. There are a number of techniques for the measurements of the vertical profiles of atmospheric turbulence [12,13,14]. To recover the vertical profiles of the atmospheric turbulence, one may use methods based on the slope detection and the ranging (Slodar) technique. It should be noted that the vertical profiles of the optical turbulence may be restored with a Slodar technique with very high resolution (10 m or less) only in the lower layers of the atmosphere.
To determine the features of the vertical structure of the surface layer at the Baykal astrophysical observatory, we performed measurements of the wavefront local slopes using a Shack–Hartmann sensor of the large solar vacuum telescope for the time period from 26 June 2018 to 28 June 2018. Main parameters of the large solar vacuum telescope and the Shack–Hartmann wavefront sensor used are given in Table 2.
Motion of images in each subaperture is estimated using extended light sources—images of solar edge. According to formula [14], we estimated cross-correlation functions:
C s ( i , j ) = δ x 1 · δ x 2 ,
where
δ x 1 ( i , 0 ) = i = m t s m ( i ) s m ( 0 ) ,
δ x 2 ( i , β ) = i = 1 m s m ( i + β z m ) s m ( β z m ) .
Indices 1 and 2 correspond to light sources (the individual spaced regions of the solar edge), z m is the height of the turbulent layer m, s m ( 0 ) is wavefront slope induced by the turbulent layer m on reference subaperture, i , j —indices of the subapertures, and β —angle between "axial" and "non axial" objects within the field of view of the given subaperture.
In calculations, there is a task to extract the local slopes of the wavefront measured by the Shack–Hartmann sensor using the solar edge images. Multiple centroiding techniques are applied to extract the local slopes including analysis of gravity centers, weighted gravity centers, filtering, and correlation based centroiding techniques. To extract the local slopes of the wavefront, we used a weighted gravity centers technique [15].
Applying the standard technique [13], we estimated the functions C s which characterize the vertical profiles of the optical turbulence in the lower layers. The technique is based on a relation between the wavefronts’ distortions spaced at the telescope pupil and altitude of the turbulent layer through which the wavefronts of a few light sources pass. Positive correlation coefficients of wavefront distortions indicate the existence of a turbulent layer. The examples of calculated functions C s to determine the vertical structure of the atmospheric turbulence are shown in Figure 7.
Analysis of the calculated functions depicted in Figure 7 shows that there are a few peaks. These peaks are induced by the individual turbulent layers. The height of the turbulent layer is determined by angle between the spaced fragments of the solar edge according to formula [14]:
z m = D n θ c o s α .
α is the zenith angle, D / n is the distance between centers of analyzed subapertures, and θ is the angle between individual spaced regions of the solar edge. The peaks correspond the angles of 5–10 and 27–35 arcsec. By triangulating the wavefront distortions and using formula (9), we estimate that the heights corresponding to θ are 1500–3000 m and 400–600 m, respectively.

4. Conclusions

The research results show that atmospheric turbulence in the surface layer at the Baykal astrophysical observatory site is most developed in the winter. At this time, the number of cases with weak turbulence is minimal. The best conditions for astronomical experiments are observed in the summer. Atmospheric turbulence is suppressed in summer. It was found that in the autumn period the distribution shape is violated: there is mainly moderate turbulence. Increased repeatability of moderate and even weak turbulence is noted in the spring. The nature of the changes of the optical turbulence is consistent with astronomical observations qualitatively. Calculations have enabled quantitative estimates of changes of the optical turbulence. Furthermore, we plan to use an improved algorithm [4] to calculate the statistics of turbulent characteristics. Increased accuracy is ensured by taking into account the shape of the turbulence spectrum at low frequencies.This will allow a detailed study of the structure of atmospheric flows and energy exchange between large and small scales in rough terrain.
In addition, we performed the measurements of the local slopes of the wavefront. By triangulation of local slopes, we reveal the heights of the turbulent layers in the atmospheric boundary layer at the Baykal astrophysical observatory site. The obtained heights (1500–3000 m and 400–600 m) are consistent with the data of lidar observations [16].
The obtained estimations of the structure characteristic of the air refractive index fluctuations as well as revealed features in the vertical distribution of the turbulence will be used for development of the spectral model of turbulent atmosphere [1,17].

Author Contributions

A.S. and M.S. performed the calculations of the distributions of the optical turbulence and developed the method to analyze wavefront distortions. V.L., P.K., and V.N. formulated solutions to the task, D.K. and A.K. provided software processing, and F.A., A.S. and E.K. performed measurements.

Funding

This research of C n 2 distributions was funded by RSF grant “19-79-00061”, the development of the method to analyze wavefront distortions was funded by President grant “075-15-2019-628”. Measurements arepartially supported the base financing of FR program II.16 using the Unique Research Facility Large Solar Vacuum Telescope http://ckp-rf.ru/usu/200615/, project “Fundamentals of breakthrough technologies in the interests of national security.”

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kovadlo, P.G.; Shikhovtsev, A.Y.; Lukin, V.P. Development of the Model of Turbulent Atmosphere at the Large Solar Vacuum Telescope Site as Applied to Image Adaptation. Atmos. Ocean. Opt. 2019, 32, 202–206. [Google Scholar] [CrossRef]
  2. Shikhovtsev, A.Y.; Kovadlo, P.G.; Lukin, V.P.; Kiselev, A.V. Dynamics of the refractive properties of the atmosphere at the Big solar vacuum telescope site: New calculating method. SPIE 2017, 10466, 104660T. [Google Scholar]
  3. Roddier, F. Atmospheric limitations to adaptive image compensation. ASP Conf. Ser. 2002, 266, 546–561. [Google Scholar]
  4. Nosov, V.; Lukin, V.; Nosov, E.; Torgaev, A.; Bogushevich, A. Measurement of atmospheric turbulence characteristics by the ultrasonic anemometers and the calibration processes. Atmosphere 2019, 10, 460. [Google Scholar] [CrossRef]
  5. Bogushevich, A.Y. Ultrasonic methods for estimation of atmospheric meteorological and turbulence parameters. Atmos. Ocean. Opt. 1999, 2, 164–169. [Google Scholar]
  6. Canut, G.; Couvreux, F.; Lothon, M.; Legain, D.; Piguet, B.; Lampert, A.; Maurel, W.; Moulin, E. Turbulence fluxes and variances measured with a sonic anemometer mounted on a tethered balloon. Atmos. Meas. Tech. 2016, 9, 4375–4386. [Google Scholar] [CrossRef] [Green Version]
  7. Gladkikh, V.A.; Nevzorova, I.V.; Oditsov, S.L.; Fedorov, V.A. Structural functions of air temperature over inhomogeneous surface. Part II. Statistics of parameters of structural functions. Atmos. Ocean. Opt. 2013, 26, 955–963. [Google Scholar]
  8. Gladkikh, V.A.; Makienko, A.E. Digital ultrasonic weather station. Pribory 2009, 7, 21–25. [Google Scholar]
  9. Mauder, M.; Zeeman, M.J. Field intercomparision of prevailing sonic anemometers. Atmos. Meas. Tech. 2018, 11, 249–263. [Google Scholar] [CrossRef]
  10. Mauder, M.; Eggert, M.; Gutsmuths, C.; Oertal, S.; Wilhelm, P.; Voelksch, I.; Wanner, L.; Tambke, J.; Bogoev, I. Comparison of turbulence measurements by a CSAT3B sonic anemometer and high-resolution bistatic Doppler lidar. Atmos. Meas. Tech. 2019. in review. [Google Scholar] [CrossRef]
  11. Laso, A.; Manana, M.; Arroyo, A.; Gonzalez, A.; Lecuna, R. A comparision of mechanical and ultrasonic anemometers for ampacity thermal rating in overhead lines. Int. Conf. Renew. Energies Power Qual. 2016, 1, 770–773. [Google Scholar] [CrossRef]
  12. Wilson, R. SLODAR: Measuring optical turbulence altitude with a Shack–Hartmann wavefront sensor. Mon. Not. R. Astron. Soc. 2002, 337, 103–108. [Google Scholar] [CrossRef]
  13. Butterley, T.; Wilson, R.; Sarazin, M. Determination of the profile of atmospheric optical turbulence strength from SLODAR data. Mon. Not. R. Astron. Soc. 2006, 369, 835–845. [Google Scholar] [CrossRef] [Green Version]
  14. Wang, Z.; Zhang, L.; Kong, L.; Bao, H.; Guo, Y.; Rao, X.; Zhong, L.; Zhu, L.; Rao, C. A modified S-DIMM+: Applying additional height grids for characterizing daytime seeing profiles. Mon. Not. R. Astron. Soc. 2018, 478, 1459–1467. [Google Scholar] [CrossRef]
  15. Lavrinov, V.V.; Lavrinova, L.N. Reconstruction of the wavefront distorted by atmospheric turbulence using a Shack–Hartmann sensor. Comput. Opt. 2019, 43, 586–595. [Google Scholar] [CrossRef]
  16. Banakh, V.A.; Smalikho, I.N. Lidar observations of atmospheric internal waves in the boundary layer of the atmosphere on the coast of Lake Baikal. Atmos. Meas. Tech. 2016, 9, 5239–5248. [Google Scholar] [CrossRef] [Green Version]
  17. Shikhovtsev, A.Y.; Kovadlo, P.G. Estimation of mean energy characteristics of atmospheric turbulence at various heights from reanalysis data. IOP Conf. Ser. Earth Environ. Sci. 2018, 211, 012023. [Google Scholar] [CrossRef]
Figure 1. Meteorological sensor mounted on the mast, large solar vacuum telescope.
Figure 1. Meteorological sensor mounted on the mast, large solar vacuum telescope.
Atmosphere 10 00661 g001
Figure 2. The repeatability of the optical turbulence at the Baykal astrophysical observatory site for the period from 1 January to 31 December 2018.
Figure 2. The repeatability of the optical turbulence at the Baykal astrophysical observatory site for the period from 1 January to 31 December 2018.
Atmosphere 10 00661 g002
Figure 3. The repeatability of the optical turbulence at the Baykal astrophysical observatory site in summer.
Figure 3. The repeatability of the optical turbulence at the Baykal astrophysical observatory site in summer.
Atmosphere 10 00661 g003
Figure 4. The repeatability of the optical turbulence at the Baykal astrophysical observatory site in autumn.
Figure 4. The repeatability of the optical turbulence at the Baykal astrophysical observatory site in autumn.
Atmosphere 10 00661 g004
Figure 5. The repeatability of the optical turbulence at the Baykal astrophysical observatory site in winter.
Figure 5. The repeatability of the optical turbulence at the Baykal astrophysical observatory site in winter.
Atmosphere 10 00661 g005
Figure 6. The repeatability of the optical turbulence at the Baykal astrophysical observatory site in spring.
Figure 6. The repeatability of the optical turbulence at the Baykal astrophysical observatory site in spring.
Atmosphere 10 00661 g006
Figure 7. The examples of the calculated functions to determine the vertical structure of the atmospheric turbulence. Subaperture diameter is equal to 0.075 m.
Figure 7. The examples of the calculated functions to determine the vertical structure of the atmospheric turbulence. Subaperture diameter is equal to 0.075 m.
Atmosphere 10 00661 g007
Table 1. The repeatability R of the optical turbulence at the Baykal astrophysical observatory site in different seasons.
Table 1. The repeatability R of the optical turbulence at the Baykal astrophysical observatory site in different seasons.
C n 2 R C n 2 R
SummerSpring
10 17 m 2 / 3 0.03 10 17 m 2 / 3 0.06
10 17 10 16 m 2 / 3 0.23 10 17 10 16 m 2 / 3 0.31
10 16 10 15 m 2 / 3 0.42 10 16 10 15 m 2 / 3 0.39
10 15 10 14 m 2 / 3 0.28 10 15 10 14 m 2 / 3 0.21
10 14 10 13 m 2 / 3 0.04 10 14 10 13 m 2 / 3 0.03
10 13 10 12 m 2 / 3 0 10 13 10 12 m 2 / 3 0
AutumnWinter
10 17 m 2 / 3 0.04 10 14 m 2 / 3 0.02
10 17 10 16 m 2 / 3 0.17 10 14 5 × 10 14 m 2 / 3 0.1
10 16 10 15 m 2 / 3 0.37 5 × 10 14 10 13 m 2 / 3 0.33
10 15 10 14 m 2 / 3 0.36 10 13 3 × 10 13 m 2 / 3 0.4
10 14 10 13 m 2 / 3 0.05 3 × 10 13 m 2 / 3 5 × 10 13 m 2 / 3 0.13
10 13 10 12 m 2 / 3 0.01 5 × 10 13 10 12 m 2 / 3 0.02
Table 2. Main parameters of the large solar vacuum telescope and Shack–Hartmann wavefront sensor.
Table 2. Main parameters of the large solar vacuum telescope and Shack–Hartmann wavefront sensor.
ParameterValue
Aperture diameter600 mm
Focal length40 m
Radiation wavelength0.535 μ m
Camera resolution640 × 480 p i x
Pixel size7.4 μ m
Number of subapertures8 × 8
Equivalent size of subaperture7.5 cm
Angular pixel size0.83
Frame frequency70 Hz
Exposure30 ms

Share and Cite

MDPI and ACS Style

Shikhovtsev, A.; Kovadlo, P.; Lukin, V.; Nosov, V.; Kiselev, A.; Kolobov, D.; Kopylov, E.; Shikhovtsev, M.; Avdeev, F. Statistics of the Optical Turbulence from the Micrometeorological Measurements at the Baykal Astrophysical Observatory Site. Atmosphere 2019, 10, 661. https://doi.org/10.3390/atmos10110661

AMA Style

Shikhovtsev A, Kovadlo P, Lukin V, Nosov V, Kiselev A, Kolobov D, Kopylov E, Shikhovtsev M, Avdeev F. Statistics of the Optical Turbulence from the Micrometeorological Measurements at the Baykal Astrophysical Observatory Site. Atmosphere. 2019; 10(11):661. https://doi.org/10.3390/atmos10110661

Chicago/Turabian Style

Shikhovtsev, Artem, Pavel Kovadlo, Vladimir Lukin, Victor Nosov, Alexander Kiselev, Dmitry Kolobov, Evgeny Kopylov, Maxim Shikhovtsev, and Fedor Avdeev. 2019. "Statistics of the Optical Turbulence from the Micrometeorological Measurements at the Baykal Astrophysical Observatory Site" Atmosphere 10, no. 11: 661. https://doi.org/10.3390/atmos10110661

APA Style

Shikhovtsev, A., Kovadlo, P., Lukin, V., Nosov, V., Kiselev, A., Kolobov, D., Kopylov, E., Shikhovtsev, M., & Avdeev, F. (2019). Statistics of the Optical Turbulence from the Micrometeorological Measurements at the Baykal Astrophysical Observatory Site. Atmosphere, 10(11), 661. https://doi.org/10.3390/atmos10110661

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop